postgis/lwgeom_sfcgal.c

Go to the documentation of this file.

/**********************************************************************
 *
 * PostGIS - Spatial Types for PostgreSQL
 * http://postgis.net
 *
 * PostGIS is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 2 of the License, or
 * (at your option) any later version.
 *
 * PostGIS is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with PostGIS.  If not, see <http://www.gnu.org/licenses/>.
 *
 **********************************************************************
 *
 * Copyright 2012-2013 Oslandia <infos@oslandia.com>
 *
 **********************************************************************/


#include "postgres.h"
#include "fmgr.h"
#include "../liblwgeom/liblwgeom.h"

#include "lwgeom_pg.h"
#include "lwgeom_sfcgal.h"
#include "../postgis_config.h"


Datum postgis_sfcgal_version(PG_FUNCTION_ARGS);

Datum sfcgal_from_ewkt(PG_FUNCTION_ARGS);
Datum sfcgal_distance(PG_FUNCTION_ARGS);
Datum sfcgal_distance3D(PG_FUNCTION_ARGS);
Datum sfcgal_area(PG_FUNCTION_ARGS);
Datum sfcgal_area3D(PG_FUNCTION_ARGS);
Datum sfcgal_intersects(PG_FUNCTION_ARGS);
Datum sfcgal_intersects3D(PG_FUNCTION_ARGS);
Datum sfcgal_intersection(PG_FUNCTION_ARGS);
Datum sfcgal_intersection3D(PG_FUNCTION_ARGS);
Datum sfcgal_difference(PG_FUNCTION_ARGS);
Datum sfcgal_difference3D(PG_FUNCTION_ARGS);
Datum sfcgal_union(PG_FUNCTION_ARGS);
Datum sfcgal_union3D(PG_FUNCTION_ARGS);
Datum sfcgal_volume(PG_FUNCTION_ARGS);
Datum sfcgal_extrude(PG_FUNCTION_ARGS);
Datum sfcgal_straight_skeleton(PG_FUNCTION_ARGS);
Datum sfcgal_approximate_medial_axis(PG_FUNCTION_ARGS);
Datum sfcgal_is_planar(PG_FUNCTION_ARGS);
Datum sfcgal_orientation(PG_FUNCTION_ARGS);
Datum sfcgal_force_lhr(PG_FUNCTION_ARGS);
Datum sfcgal_triangulate(PG_FUNCTION_ARGS);
Datum sfcgal_tesselate(PG_FUNCTION_ARGS);
Datum sfcgal_minkowski_sum(PG_FUNCTION_ARGS);
Datum sfcgal_make_solid(PG_FUNCTION_ARGS);
Datum sfcgal_is_solid(PG_FUNCTION_ARGS);


GSERIALIZED *geometry_serialize(LWGEOM *lwgeom);
char* text2cstring(const text *textptr);

static int __sfcgal_init = 0;

void sfcgal_postgis_init(void)
{
    if ( ! __sfcgal_init ) {
        sfcgal_init();
        sfcgal_set_error_handlers((sfcgal_error_handler_t) lwpgnotice, (sfcgal_error_handler_t) lwpgerror);
        sfcgal_set_alloc_handlers(lwalloc, lwfree);
        __sfcgal_init = 1;
    }
}


/* Conversion from GSERIALIZED* to SFCGAL::Geometry */
sfcgal_geometry_t* POSTGIS2SFCGALGeometry(GSERIALIZED *pglwgeom)
{
        sfcgal_geometry_t* g;
        LWGEOM *lwgeom = lwgeom_from_gserialized(pglwgeom);

        if (! lwgeom)
        {
                lwpgerror("POSTGIS2SFCGALGeometry: Unable to deserialize input");
        }
        g = LWGEOM2SFCGAL(lwgeom);
        lwgeom_free(lwgeom);

        return g;
}


/* Conversion from GSERIALIZED* to SFCGAL::PreparedGeometry */
sfcgal_prepared_geometry_t* POSTGIS2SFCGALPreparedGeometry(GSERIALIZED *pglwgeom)
{
        sfcgal_geometry_t* g;
        LWGEOM *lwgeom = lwgeom_from_gserialized(pglwgeom);

        if (!lwgeom)
        {
                lwpgerror("POSTGIS2SFCGALPreparedGeometry: Unable to deserialize input");
        }
        g = LWGEOM2SFCGAL(lwgeom);

        lwgeom_free(lwgeom);

        return sfcgal_prepared_geometry_create_from_geometry(g, gserialized_get_srid(pglwgeom));
}


/* Conversion from SFCGAL::Geometry to GSERIALIZED */
GSERIALIZED* SFCGALGeometry2POSTGIS(const sfcgal_geometry_t* geom, int force3D, int SRID)
{
        GSERIALIZED *result;
        LWGEOM* lwgeom = SFCGAL2LWGEOM(geom, force3D, SRID);

        if (lwgeom_needs_bbox(lwgeom) == LW_TRUE)
                lwgeom_add_bbox(lwgeom);

        result = geometry_serialize(lwgeom);
        lwgeom_free(lwgeom);

        return result;
}


/* Conversion from SFCGAL::PreparedGeometry to GSERIALIZED */
GSERIALIZED* SFCGALPreparedGeometry2POSTGIS(const sfcgal_prepared_geometry_t* geom, int force3D)
{
    return SFCGALGeometry2POSTGIS(sfcgal_prepared_geometry_geometry(geom),
                force3D, sfcgal_prepared_geometry_srid(geom));
}


/* Conversion from EWKT to GSERIALIZED */
PG_FUNCTION_INFO_V1(sfcgal_from_ewkt);
Datum sfcgal_from_ewkt(PG_FUNCTION_ARGS)
{
        GSERIALIZED* result;
        sfcgal_prepared_geometry_t* g;
        text *wkttext = PG_GETARG_TEXT_P(0);
        char *cstring = text2cstring(wkttext);

        sfcgal_postgis_init();

        g = sfcgal_io_read_ewkt( cstring, strlen(cstring) );

        result = SFCGALPreparedGeometry2POSTGIS( g, 0 );
        sfcgal_prepared_geometry_delete( g );
        PG_RETURN_POINTER(result);
}


PG_FUNCTION_INFO_V1(sfcgal_area);
Datum sfcgal_area(PG_FUNCTION_ARGS)
    {
        GSERIALIZED *input;
        sfcgal_geometry_t *geom;
        double result;

        sfcgal_postgis_init();

        input = PG_GETARG_GSERIALIZED_P(0);
        geom = POSTGIS2SFCGALGeometry(input);

        result = sfcgal_geometry_area(geom);
        sfcgal_geometry_delete(geom);

        PG_FREE_IF_COPY(input, 0);

        PG_RETURN_FLOAT8(result);
}


PG_FUNCTION_INFO_V1(sfcgal_area3D);
Datum sfcgal_area3D(PG_FUNCTION_ARGS)
    {
        GSERIALIZED *input;
        sfcgal_geometry_t *geom;
        double result;

        sfcgal_postgis_init();

        input = PG_GETARG_GSERIALIZED_P(0);
        geom = POSTGIS2SFCGALGeometry(input);

        result = sfcgal_geometry_area_3d(geom);
        sfcgal_geometry_delete(geom);

        PG_FREE_IF_COPY(input, 0);

        PG_RETURN_FLOAT8(result);
}


PG_FUNCTION_INFO_V1(sfcgal_is_planar);
Datum sfcgal_is_planar(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input;
        sfcgal_geometry_t *geom;
        int result;

        sfcgal_postgis_init();

        input = PG_GETARG_GSERIALIZED_P(0);
        geom = POSTGIS2SFCGALGeometry(input);

        result = sfcgal_geometry_is_planar(geom);
        sfcgal_geometry_delete(geom);

        PG_FREE_IF_COPY(input, 0);

        PG_RETURN_BOOL(result);
}


PG_FUNCTION_INFO_V1(sfcgal_orientation);
Datum sfcgal_orientation(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input;
        sfcgal_geometry_t *geom;
        int result;

        sfcgal_postgis_init();

        input = PG_GETARG_GSERIALIZED_P(0);
        geom = POSTGIS2SFCGALGeometry(input);

        result = sfcgal_geometry_orientation(geom);
        sfcgal_geometry_delete(geom);

        PG_FREE_IF_COPY(input, 0);

        PG_RETURN_INT32(result);
}


PG_FUNCTION_INFO_V1(sfcgal_intersects);
Datum sfcgal_intersects(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input0, *input1;
        sfcgal_geometry_t *geom0, *geom1;
        int result;

        sfcgal_postgis_init();

        input0 = PG_GETARG_GSERIALIZED_P(0);
        input1 = PG_GETARG_GSERIALIZED_P(1);
        geom0 = POSTGIS2SFCGALGeometry(input0);
        PG_FREE_IF_COPY(input0, 0);
        geom1 = POSTGIS2SFCGALGeometry(input1);
        PG_FREE_IF_COPY(input1, 1);

        result = sfcgal_geometry_intersects(geom0, geom1);
        sfcgal_geometry_delete(geom0);
        sfcgal_geometry_delete(geom1);

        PG_RETURN_BOOL(result);
}


PG_FUNCTION_INFO_V1(sfcgal_intersects3D);
Datum sfcgal_intersects3D(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input0, *input1;
        sfcgal_geometry_t *geom0, *geom1;
        int result;

        sfcgal_postgis_init();

        input0 = PG_GETARG_GSERIALIZED_P(0);
        input1 = PG_GETARG_GSERIALIZED_P(1);
        geom0 = POSTGIS2SFCGALGeometry(input0);
        PG_FREE_IF_COPY(input0, 0);
        geom1 = POSTGIS2SFCGALGeometry(input1);
        PG_FREE_IF_COPY(input1, 1);

        result = sfcgal_geometry_intersects_3d(geom0, geom1);
        sfcgal_geometry_delete(geom0);
        sfcgal_geometry_delete(geom1);

        PG_RETURN_BOOL(result);
}


PG_FUNCTION_INFO_V1(sfcgal_distance);
Datum sfcgal_distance(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input0, *input1;
        sfcgal_geometry_t *geom0, *geom1;
        double result;

        sfcgal_postgis_init();

        input0 = PG_GETARG_GSERIALIZED_P(0);
        input1 = PG_GETARG_GSERIALIZED_P(1);
        geom0 = POSTGIS2SFCGALGeometry(input0);
        PG_FREE_IF_COPY(input0, 0);
        geom1 = POSTGIS2SFCGALGeometry(input1);
        PG_FREE_IF_COPY(input1, 1);

        result = sfcgal_geometry_distance(geom0, geom1);
        sfcgal_geometry_delete(geom0);
        sfcgal_geometry_delete(geom1);

        PG_RETURN_FLOAT8(result);
}


PG_FUNCTION_INFO_V1(sfcgal_distance3D);
Datum sfcgal_distance3D(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input0, *input1;
        sfcgal_geometry_t *geom0, *geom1;
        double result;

        sfcgal_postgis_init();

        input0 = PG_GETARG_GSERIALIZED_P(0);
        input1 = PG_GETARG_GSERIALIZED_P(1);
        geom0 = POSTGIS2SFCGALGeometry(input0);
        PG_FREE_IF_COPY(input0, 0);
        geom1 = POSTGIS2SFCGALGeometry(input1);
        PG_FREE_IF_COPY(input1, 1);

        result = sfcgal_geometry_distance_3d(geom0, geom1);
        sfcgal_geometry_delete(geom0);
        sfcgal_geometry_delete(geom1);

        PG_RETURN_FLOAT8(result);
}


PG_FUNCTION_INFO_V1(sfcgal_tesselate);
Datum sfcgal_tesselate(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input, *output;
        sfcgal_geometry_t *geom;
        sfcgal_geometry_t *result;
        srid_t srid;

        sfcgal_postgis_init();

        input = PG_GETARG_GSERIALIZED_P(0);
        srid = gserialized_get_srid(input);
        geom = POSTGIS2SFCGALGeometry(input);
        PG_FREE_IF_COPY(input, 0);

        result = sfcgal_geometry_tesselate(geom);
        sfcgal_geometry_delete(geom);

        output = SFCGALGeometry2POSTGIS(result, 0, srid);
        sfcgal_geometry_delete(result);

        PG_RETURN_POINTER(output);
}


PG_FUNCTION_INFO_V1(sfcgal_triangulate);
Datum sfcgal_triangulate(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input, *output;
        sfcgal_geometry_t *geom;
        sfcgal_geometry_t *result;
        srid_t srid;

        sfcgal_postgis_init();

        input = PG_GETARG_GSERIALIZED_P(0);
        srid = gserialized_get_srid(input);
        geom = POSTGIS2SFCGALGeometry(input);
        PG_FREE_IF_COPY(input, 0);

        result = sfcgal_geometry_triangulate_2dz(geom);
        sfcgal_geometry_delete(geom);

        output = SFCGALGeometry2POSTGIS(result, 0, srid);
        sfcgal_geometry_delete(result);

        PG_RETURN_POINTER(output);
}


PG_FUNCTION_INFO_V1(sfcgal_force_lhr);
Datum sfcgal_force_lhr(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input, *output;
        sfcgal_geometry_t *geom;
        sfcgal_geometry_t *result;
        srid_t srid;

        sfcgal_postgis_init();

        input = PG_GETARG_GSERIALIZED_P(0);
        srid = gserialized_get_srid(input);
        geom = POSTGIS2SFCGALGeometry(input);
        PG_FREE_IF_COPY(input, 0);

        result = sfcgal_geometry_force_lhr(geom);
        sfcgal_geometry_delete(geom);

        output = SFCGALGeometry2POSTGIS(result, 0, srid);
        sfcgal_geometry_delete(result);

        PG_RETURN_POINTER(output);
}


PG_FUNCTION_INFO_V1(sfcgal_straight_skeleton);
Datum sfcgal_straight_skeleton(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input, *output;
        sfcgal_geometry_t *geom;
        sfcgal_geometry_t *result;
        srid_t srid;

        sfcgal_postgis_init();

        input = PG_GETARG_GSERIALIZED_P(0);
        srid = gserialized_get_srid(input);
        geom = POSTGIS2SFCGALGeometry(input);
        PG_FREE_IF_COPY(input, 0);

        result = sfcgal_geometry_straight_skeleton(geom);
        sfcgal_geometry_delete(geom);

        output = SFCGALGeometry2POSTGIS(result, 0, srid);
        sfcgal_geometry_delete(result);

        PG_RETURN_POINTER(output);
}

PG_FUNCTION_INFO_V1(sfcgal_approximate_medial_axis);
Datum sfcgal_approximate_medial_axis(PG_FUNCTION_ARGS)
{
#if POSTGIS_SFCGAL_VERSION < 12
        lwpgerror("The SFCGAL version this PostGIS binary "
                "was compiled against (%d) doesn't support "
                "'sfcgal_geometry_approximate_medial_axis' function (1.2.0+ required)",
                POSTGIS_SFCGAL_VERSION);
        PG_RETURN_NULL();
#else /* POSTGIS_SFCGAL_VERSION >= 12 */
        GSERIALIZED *input, *output;
        sfcgal_geometry_t *geom;
        sfcgal_geometry_t *result;
        srid_t srid;

        sfcgal_postgis_init();

        input = PG_GETARG_GSERIALIZED_P(0);
        srid = gserialized_get_srid(input);
        geom = POSTGIS2SFCGALGeometry(input);
        PG_FREE_IF_COPY(input, 0);

        result = sfcgal_geometry_approximate_medial_axis(geom);
        sfcgal_geometry_delete(geom);

        output = SFCGALGeometry2POSTGIS(result, 0, srid);
        sfcgal_geometry_delete(result);

        PG_RETURN_POINTER(output);
#endif /* POSTGIS_SFCGAL_VERSION >= 12 */
}


PG_FUNCTION_INFO_V1(sfcgal_intersection);
Datum sfcgal_intersection(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input0, *input1, *output;
        sfcgal_geometry_t *geom0, *geom1;
        sfcgal_geometry_t *result;
        srid_t srid;

        sfcgal_postgis_init();

        input0 = PG_GETARG_GSERIALIZED_P(0);
        srid = gserialized_get_srid(input0);
        input1 = PG_GETARG_GSERIALIZED_P(1);
        geom0 = POSTGIS2SFCGALGeometry(input0);
        PG_FREE_IF_COPY(input0, 0);
        geom1 = POSTGIS2SFCGALGeometry(input1);
        PG_FREE_IF_COPY(input1, 1);

        result = sfcgal_geometry_intersection(geom0, geom1);
        sfcgal_geometry_delete(geom0);
        sfcgal_geometry_delete(geom1);

        output = SFCGALGeometry2POSTGIS(result, 0, srid);
        sfcgal_geometry_delete(result);

        PG_RETURN_POINTER(output);
}


PG_FUNCTION_INFO_V1(sfcgal_intersection3D);
Datum sfcgal_intersection3D(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input0, *input1, *output;
        sfcgal_geometry_t *geom0, *geom1;
        sfcgal_geometry_t *result;
        srid_t srid;

        sfcgal_postgis_init();

        input0 = PG_GETARG_GSERIALIZED_P(0);
        srid = gserialized_get_srid(input0);
        input1 = PG_GETARG_GSERIALIZED_P(1);
        geom0 = POSTGIS2SFCGALGeometry(input0);
        PG_FREE_IF_COPY(input0, 0);
        geom1 = POSTGIS2SFCGALGeometry(input1);
        PG_FREE_IF_COPY(input1, 1);

        result = sfcgal_geometry_intersection_3d(geom0, geom1);
        sfcgal_geometry_delete(geom0);
        sfcgal_geometry_delete(geom1);

        output = SFCGALGeometry2POSTGIS(result, 0, srid);
        sfcgal_geometry_delete(result);

        PG_RETURN_POINTER(output);
}

PG_FUNCTION_INFO_V1(sfcgal_difference);
Datum sfcgal_difference(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input0, *input1, *output;
        sfcgal_geometry_t *geom0, *geom1;
        sfcgal_geometry_t *result;
        srid_t srid;

        sfcgal_postgis_init();

        input0 = (GSERIALIZED*) PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
        srid = gserialized_get_srid(input0);
        input1 = (GSERIALIZED*) PG_DETOAST_DATUM(PG_GETARG_DATUM(1));
        geom0 = POSTGIS2SFCGALGeometry(input0);
        PG_FREE_IF_COPY(input0, 0);
        geom1 = POSTGIS2SFCGALGeometry(input1);
        PG_FREE_IF_COPY(input1, 1);

        result = sfcgal_geometry_difference(geom0, geom1);
        sfcgal_geometry_delete(geom0);
        sfcgal_geometry_delete(geom1);

        output = SFCGALGeometry2POSTGIS(result, 0, srid);
        sfcgal_geometry_delete(result);

        PG_RETURN_POINTER(output);
}


PG_FUNCTION_INFO_V1(sfcgal_difference3D);
Datum sfcgal_difference3D(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input0, *input1, *output;
        sfcgal_geometry_t *geom0, *geom1;
        sfcgal_geometry_t *result;
        srid_t srid;

        sfcgal_postgis_init();

        input0 = (GSERIALIZED*) PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
        srid = gserialized_get_srid(input0);
        input1 = (GSERIALIZED*) PG_DETOAST_DATUM(PG_GETARG_DATUM(1));
        geom0 = POSTGIS2SFCGALGeometry(input0);
        PG_FREE_IF_COPY(input0, 0);
        geom1 = POSTGIS2SFCGALGeometry(input1);
        PG_FREE_IF_COPY(input1, 1);

        result = sfcgal_geometry_difference_3d(geom0, geom1);
        sfcgal_geometry_delete(geom0);
        sfcgal_geometry_delete(geom1);

        output = SFCGALGeometry2POSTGIS(result, 0, srid);
        sfcgal_geometry_delete(result);

        PG_RETURN_POINTER(output);
}

PG_FUNCTION_INFO_V1(sfcgal_union);
Datum sfcgal_union(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input0, *input1, *output;
        sfcgal_geometry_t *geom0, *geom1;
        sfcgal_geometry_t *result;
        srid_t srid;

        sfcgal_postgis_init();

        input0 = (GSERIALIZED*) PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
        srid = gserialized_get_srid(input0);
        input1 = (GSERIALIZED*) PG_DETOAST_DATUM(PG_GETARG_DATUM(1));
        geom0 = POSTGIS2SFCGALGeometry(input0);
        PG_FREE_IF_COPY(input0, 0);
        geom1 = POSTGIS2SFCGALGeometry(input1);
        PG_FREE_IF_COPY(input1, 1);

        result = sfcgal_geometry_union(geom0, geom1);
        sfcgal_geometry_delete(geom0);
        sfcgal_geometry_delete(geom1);

        output = SFCGALGeometry2POSTGIS(result, 0, srid);
        sfcgal_geometry_delete(result);

        PG_RETURN_POINTER(output);
}


PG_FUNCTION_INFO_V1(sfcgal_union3D);
Datum sfcgal_union3D(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input0, *input1, *output;
        sfcgal_geometry_t *geom0, *geom1;
        sfcgal_geometry_t *result;
        srid_t srid;

        sfcgal_postgis_init();

        input0 = (GSERIALIZED*) PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
        srid = gserialized_get_srid(input0);
        input1 = (GSERIALIZED*) PG_DETOAST_DATUM(PG_GETARG_DATUM(1));
        geom0 = POSTGIS2SFCGALGeometry(input0);
        PG_FREE_IF_COPY(input0, 0);
        geom1 = POSTGIS2SFCGALGeometry(input1);
        PG_FREE_IF_COPY(input1, 1);

        result = sfcgal_geometry_union_3d(geom0, geom1);
        sfcgal_geometry_delete(geom0);
        sfcgal_geometry_delete(geom1);

        output = SFCGALGeometry2POSTGIS(result, 0, srid);
        sfcgal_geometry_delete(result);

        PG_RETURN_POINTER(output);
}

PG_FUNCTION_INFO_V1(sfcgal_volume);
Datum sfcgal_volume(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input;
        sfcgal_geometry_t *geom;
        double result;

        sfcgal_postgis_init();

        input = (GSERIALIZED*) PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
        geom = POSTGIS2SFCGALGeometry(input);

        result = sfcgal_geometry_volume(geom);
        sfcgal_geometry_delete(geom);

        PG_FREE_IF_COPY(input, 0);

        PG_RETURN_FLOAT8(result);
}

PG_FUNCTION_INFO_V1(sfcgal_minkowski_sum);
Datum sfcgal_minkowski_sum(PG_FUNCTION_ARGS)
{
        GSERIALIZED *input0, *input1, *output;
        sfcgal_geometry_t *geom0, *geom1;
        sfcgal_geometry_t *result;
        srid_t srid;

        sfcgal_postgis_init();

        input0 = PG_GETARG_GSERIALIZED_P(0);
        srid = gserialized_get_srid(input0);
        input1 = PG_GETARG_GSERIALIZED_P(1);
        geom0 = POSTGIS2SFCGALGeometry(input0);
        PG_FREE_IF_COPY(input0, 0);
        geom1 = POSTGIS2SFCGALGeometry(input1);
        PG_FREE_IF_COPY(input1, 1);

        result = sfcgal_geometry_minkowski_sum(geom0, geom1);
        sfcgal_geometry_delete(geom0);
        sfcgal_geometry_delete(geom1);

        output = SFCGALGeometry2POSTGIS(result, 0, srid);
        sfcgal_geometry_delete(result);

        PG_RETURN_POINTER(output);
}


PG_FUNCTION_INFO_V1(sfcgal_extrude);
Datum sfcgal_extrude(PG_FUNCTION_ARGS)
{
    GSERIALIZED *input, *output;
    sfcgal_geometry_t *geom;
    sfcgal_geometry_t *result;
    double dx, dy, dz;
    srid_t srid;

    sfcgal_postgis_init();

    input = PG_GETARG_GSERIALIZED_P(0);
    srid = gserialized_get_srid(input);

    geom = POSTGIS2SFCGALGeometry(input);
    PG_FREE_IF_COPY(input, 0);

    dx = PG_GETARG_FLOAT8(1);
    dy = PG_GETARG_FLOAT8(2);
    dz = PG_GETARG_FLOAT8(3);

    result = sfcgal_geometry_extrude(geom, dx, dy, dz);
    sfcgal_geometry_delete(geom);

    output = SFCGALGeometry2POSTGIS(result, 0, srid);
    sfcgal_geometry_delete(result);

    PG_RETURN_POINTER(output);
}

PG_FUNCTION_INFO_V1(postgis_sfcgal_version);
Datum postgis_sfcgal_version(PG_FUNCTION_ARGS)
{
        const char *ver = lwgeom_sfcgal_version();
        text *result = cstring2text(ver);
        PG_RETURN_POINTER(result);
}

PG_FUNCTION_INFO_V1(sfcgal_is_solid);
Datum sfcgal_is_solid(PG_FUNCTION_ARGS)
{
        int result;
        GSERIALIZED *input = PG_GETARG_GSERIALIZED_P(0);
        LWGEOM *lwgeom = lwgeom_from_gserialized(input);
        PG_FREE_IF_COPY(input, 0);
        if (! lwgeom)
        {
                elog(ERROR, "sfcgal_is_solid: Unable to deserialize input");
        }
        result = FLAGS_GET_SOLID( lwgeom->flags );

        lwgeom_free(lwgeom);

        PG_RETURN_BOOL(result);
}

PG_FUNCTION_INFO_V1(sfcgal_make_solid);
Datum sfcgal_make_solid(PG_FUNCTION_ARGS)
{
        GSERIALIZED *output;
        GSERIALIZED *input = PG_GETARG_GSERIALIZED_P(0);
        LWGEOM *lwgeom = lwgeom_from_gserialized(input);
        if (! lwgeom)
        {
                elog(ERROR, "sfcgal_make_solid: Unable to deserialize input");
        }

        FLAGS_SET_SOLID( lwgeom->flags, 1);

        output = geometry_serialize( lwgeom );
        lwgeom_free(lwgeom);
        PG_FREE_IF_COPY(input, 0);
        PG_RETURN_POINTER(output);
}