topologypreservesimplify()

Datum topologypreservesimplify

(

PG_FUNCTION_ARGS

)

Definition at line 832 of file postgis/lwgeom_geos.c.

{
        GSERIALIZED     *gs1;
        LWGEOM *lwg1;
        double  tolerance;
        GEOSGeometry *g1, *g3;
        GSERIALIZED *result;
        uint32_t type;
 
        gs1 = PG_GETARG_GSERIALIZED_P(0);
        tolerance = PG_GETARG_FLOAT8(1);
        lwg1 = lwgeom_from_gserialized(gs1);
 
        /* Empty.Simplify() == Empty */
        type = lwgeom_get_type(lwg1);
        if (lwgeom_is_empty(lwg1) || type == TINTYPE || type == TRIANGLETYPE)
                PG_RETURN_POINTER(gs1);
 
        if (!lwgeom_isfinite(lwg1))
        {
                lwpgerror("Geometry contains invalid coordinates");
                PG_RETURN_NULL();
        }
 
        initGEOS(lwpgnotice, lwgeom_geos_error);
 
        g1 = LWGEOM2GEOS(lwg1, LW_TRUE);
        lwgeom_free(lwg1);
        if (!g1)
                HANDLE_GEOS_ERROR("First argument geometry could not be converted to GEOS");
 
        g3 = GEOSTopologyPreserveSimplify(g1,tolerance);
        GEOSGeom_destroy(g1);
 
        if (!g3) HANDLE_GEOS_ERROR("GEOSTopologyPreserveSimplify");
 
        POSTGIS_DEBUGF(3, "result: %s", GEOSGeomToWKT(g3));
 
        GEOSSetSRID(g3, gserialized_get_srid(gs1));
 
        result = GEOS2POSTGIS(g3, gserialized_has_z(gs1));
        GEOSGeom_destroy(g3);
 
        if (!result)
        {
                elog(ERROR,"GEOS topologypreservesimplify() threw an error (result postgis geometry formation)!");
                PG_RETURN_NULL(); /* never get here */
        }
 
        PG_FREE_IF_COPY(gs1, 0);
        PG_RETURN_POINTER(result);
}

References GEOS2POSTGIS(), gserialized_get_srid(), gserialized_has_z(), HANDLE_GEOS_ERROR, LW_TRUE, LWGEOM2GEOS(), lwgeom_free(), lwgeom_from_gserialized(), lwgeom_geos_error(), lwgeom_get_type(), lwgeom_is_empty(), lwgeom_isfinite(), TINTYPE, and TRIANGLETYPE.

Here is the call graph for this function: