lwgeom_make_valid_params()

LWGEOM* lwgeom_make_valid_params	(	LWGEOM *	geom,
		char *	make_valid_params
	)

Definition at line 906 of file liblwgeom/lwgeom_geos_clean.c.

{
        int is3d;
        GEOSGeom geosgeom;
        GEOSGeometry* geosout;
        LWGEOM* lwgeom_out;
 
        LWDEBUG(1, "lwgeom_make_valid enter");
 
        is3d = FLAGS_GET_Z(lwgeom_in->flags);
 
        /*
         * Step 1 : try to convert to GEOS, if impossible, clean that up first
         *          otherwise (adding only duplicates of existing points)
         */
 
        initGEOS(lwgeom_geos_error, lwgeom_geos_error);
 
        lwgeom_out = lwgeom_make_geos_friendly(lwgeom_in);
        if (!lwgeom_out) lwerror("Could not make a geos friendly geometry out of input");
 
        LWDEBUGF(4, "Input geom %p made GEOS-valid as %p", lwgeom_in, lwgeom_out);
 
        geosgeom = LWGEOM2GEOS(lwgeom_out, 1);
        if ( lwgeom_in != lwgeom_out ) {
                lwgeom_free(lwgeom_out);
        }
        if (!geosgeom)
        {
                lwerror("Couldn't convert POSTGIS geom to GEOS: %s", lwgeom_geos_errmsg);
                return NULL;
        }
        else
        {
                LWDEBUG(4, "geom converted to GEOS");
        }
 
#if POSTGIS_GEOS_VERSION < 30800
        geosout = LWGEOM_GEOS_makeValid(geosgeom);
#elif POSTGIS_GEOS_VERSION < 31000
        geosout = GEOSMakeValid(geosgeom);
#else
        if (!make_valid_params) {
                geosout = GEOSMakeValid(geosgeom);
        }
        else {
                /*
                * Set up a parameters object for this
                * make valid operation before calling
                * it
                */
                const char *value;
                char *param_list[OPTION_LIST_SIZE];
                char param_list_text[OPTION_LIST_SIZE];
                strncpy(param_list_text, make_valid_params, OPTION_LIST_SIZE-1);
                param_list_text[OPTION_LIST_SIZE-1] = '\0'; /* ensure null-termination */
                memset(param_list, 0, sizeof(param_list));
                option_list_parse(param_list_text, param_list);
                GEOSMakeValidParams *params = GEOSMakeValidParams_create();
                value = option_list_search(param_list, "method");
                if (value) {
                        if (strcasecmp(value, "linework") == 0) {
                                GEOSMakeValidParams_setMethod(params, GEOS_MAKE_VALID_LINEWORK);
                        }
                        else if (strcasecmp(value, "structure") == 0) {
                                GEOSMakeValidParams_setMethod(params, GEOS_MAKE_VALID_STRUCTURE);
                        }
                        else
                        {
                                GEOSMakeValidParams_destroy(params);
                                lwerror("Unsupported value for 'method', '%s'. Use 'linework' or 'structure'.", value);
                        }
                }
                value = option_list_search(param_list, "keepcollapsed");
                if (value) {
                        if (strcasecmp(value, "true") == 0) {
                                GEOSMakeValidParams_setKeepCollapsed(params, 1);
                        }
                        else if (strcasecmp(value, "false") == 0) {
                                GEOSMakeValidParams_setKeepCollapsed(params, 0);
                        }
                        else
                        {
                                GEOSMakeValidParams_destroy(params);
                                lwerror("Unsupported value for 'keepcollapsed', '%s'. Use 'true' or 'false'", value);
                        }
                }
                geosout = GEOSMakeValidWithParams(geosgeom, params);
                GEOSMakeValidParams_destroy(params);
        }
#endif
        GEOSGeom_destroy(geosgeom);
        if (!geosout) return NULL;
 
        lwgeom_out = GEOS2LWGEOM(geosout, is3d);
        GEOSGeom_destroy(geosout);
 
        if (lwgeom_is_collection(lwgeom_in) && !lwgeom_is_collection(lwgeom_out))
        {
                LWGEOM** ogeoms = lwalloc(sizeof(LWGEOM*));
                LWGEOM* ogeom;
                LWDEBUG(3, "lwgeom_make_valid: forcing multi");
                /* NOTE: this is safe because lwgeom_out is surely not lwgeom_in or
                 * otherwise we couldn't have a collection and a non-collection */
                assert(lwgeom_in != lwgeom_out);
                ogeoms[0] = lwgeom_out;
                ogeom = (LWGEOM*)lwcollection_construct(
                    MULTITYPE[lwgeom_out->type], lwgeom_out->srid, lwgeom_out->bbox, 1, ogeoms);
                lwgeom_out->bbox = NULL;
                lwgeom_out = ogeom;
        }
 
        lwgeom_out->srid = lwgeom_in->srid;
        return lwgeom_out;
}

References LWGEOM::bbox, LWGEOM::flags, FLAGS_GET_Z, GEOS2LWGEOM(), lwalloc(), lwcollection_construct(), LWDEBUG, LWDEBUGF, lwerror(), LWGEOM2GEOS(), lwgeom_free(), lwgeom_geos_errmsg, lwgeom_geos_error(), lwgeom_is_collection(), lwgeom_make_geos_friendly(), MULTITYPE, option_list_parse(), option_list_search(), OPTION_LIST_SIZE, LWGEOM::srid, LWGEOM::type, and genraster::value.

Referenced by lwgeom_make_valid(), and ST_MakeValid().

Here is the call graph for this function:

Here is the caller graph for this function: