char *msudf_simplifyPreserveTopology(UDF_INIT *initid,UDF_ARGS *args, char *buf,
unsigned long *length, char *is_null, char *error)
{
char *result;
GEOSGeom geom2;
GEOSGeom geom1;
double tolerance;
DEBUG("msudf_simplifyPreserveTopology");
geom1 = msudf_getGeometry((unsigned char *)args->args[0],args->lengths[0]);
if (geom1 == NULL) {
strcpy(error,"Invalid geometry.");
*is_null = 1;
return NULL;
}
tolerance = *(double*) args->args[1] ;
geom2 = GEOSTopologyPreserveSimplify(geom1,tolerance);
if (geom2 != NULL) {
GEOSSetSRID(geom2,GEOSGetSRID(geom1));
result = msudf_returnGeometry(initid,length,geom2);
GEOSGeom_destroy(geom1);
GEOSGeom_destroy(geom2);
return result;
} else {
GEOSGeom_destroy(geom1);
*is_null = 1;
return NULL;
}
}
shapeObj *msGEOSTopologyPreservingSimplify(shapeObj *shape, double tolerance)
{
#ifdef USE_GEOS
GEOSGeom g1, g2;
if(!shape)
return NULL;
if(!shape->geometry) /* if no geometry for the shape then build one */
shape->geometry = (GEOSGeom) msGEOSShape2Geometry(shape);
g1 = (GEOSGeom) shape->geometry;
if(!g1) return NULL;
g2 = GEOSTopologyPreserveSimplify(g1, tolerance);
return msGEOSGeometry2Shape(g2);
#else
msSetError(MS_GEOSERR, "GEOS Simplify support is not available.", "msGEOSTopologyPreservingSimplify()");
return NULL;
#endif
}
static ERL_NIF_TERM
topology_preserve_simplify(ErlNifEnv* env, int argc, const ERL_NIF_TERM argv[])
{
GEOSGeometry **geom;
GEOSGeometry *simpler_geom;
double dbl_tol;
int int_tol;
ERL_NIF_TERM eterm;
if(!enif_get_resource(env, argv[0], GEOSGEOM_RESOURCE, (void**)&geom)) {
return 0;
}
if (enif_get_int(env, argv[1], &int_tol)) {
dbl_tol = int_tol;
}
else if (!enif_get_double(env, argv[1], &dbl_tol)) {
return 0;
}
simpler_geom = GEOSTopologyPreserveSimplify(*geom, dbl_tol);
eterm = geom_to_eterm(env, simpler_geom);
GEOSGeom_destroy(simpler_geom);
return eterm;
}
