LWGEOM_GEOS_makeValid()

static GEOSGeometry * LWGEOM_GEOS_makeValid ( const GEOSGeometry *  gin )

static

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

References GEOS2LWGEOM(), LWDEBUGF, lwerror(), lwgeom_geos_errmsg, LWGEOM_GEOS_makeValidCollection(), LWGEOM_GEOS_makeValidLine(), LWGEOM_GEOS_makeValidMultiLine(), LWGEOM_GEOS_makeValidPolygon(), lwgeom_to_ewkt(), and lwnotice().

Referenced by LWGEOM_GEOS_makeValidCollection(), LWGEOM_GEOS_makeValidMultiLine(), and lwgeom_make_valid().

{
        GEOSGeometry* gout;
        char ret_char;

        /*
         * Step 2: return what we got so far if already valid
         */

        ret_char = GEOSisValid(gin);
        if ( ret_char == 2 )
        {
                /* I don't think should ever happen */
                lwerror("GEOSisValid(): %s", lwgeom_geos_errmsg);
                return NULL;
        }
        else if ( ret_char )
        {
                LWDEBUGF(3,
                               "Geometry [%s] is valid. ",
                               lwgeom_to_ewkt(GEOS2LWGEOM(gin, 0)));

                /* It's valid at this step, return what we have */
                return GEOSGeom_clone(gin);
        }

        LWDEBUGF(3,
                       "Geometry [%s] is still not valid: %s. "
                       "Will try to clean up further.",
                       lwgeom_to_ewkt(GEOS2LWGEOM(gin, 0)), lwgeom_geos_errmsg);



        /*
         * Step 3 : make what we got valid
         */

        switch (GEOSGeomTypeId(gin))
        {
        case GEOS_MULTIPOINT:
        case GEOS_POINT:
                /* points are always valid, but we might have invalid ordinate values */
                lwnotice("PUNTUAL geometry resulted invalid to GEOS -- dunno how to clean that up");
                return NULL;
                break;

        case GEOS_LINESTRING:
                gout = LWGEOM_GEOS_makeValidLine(gin);
                if ( ! gout )  /* an exception or something */
                {
                        /* cleanup and throw */
                        lwerror("%s", lwgeom_geos_errmsg);
                        return NULL;
                }
                break; /* we've done */

        case GEOS_MULTILINESTRING:
                gout = LWGEOM_GEOS_makeValidMultiLine(gin);
                if ( ! gout )  /* an exception or something */
                {
                        /* cleanup and throw */
                        lwerror("%s", lwgeom_geos_errmsg);
                        return NULL;
                }
                break; /* we've done */

        case GEOS_POLYGON:
        case GEOS_MULTIPOLYGON:
        {
                gout = LWGEOM_GEOS_makeValidPolygon(gin);
                if ( ! gout )  /* an exception or something */
                {
                        /* cleanup and throw */
                        lwerror("%s", lwgeom_geos_errmsg);
                        return NULL;
                }
                break; /* we've done */
        }

        case GEOS_GEOMETRYCOLLECTION:
        {
                gout = LWGEOM_GEOS_makeValidCollection(gin);
                if ( ! gout )  /* an exception or something */
                {
                        /* cleanup and throw */
                        lwerror("%s", lwgeom_geos_errmsg);
                        return NULL;
                }
                break; /* we've done */
        }

        default:
        {
                char* typname = GEOSGeomType(gin);
                lwnotice("ST_MakeValid: doesn't support geometry type: %s",
                         typname);
                GEOSFree(typname);
                return NULL;
                break;
        }
        }

#if PARANOIA_LEVEL > 1
        /*
         * Now check if every point of input is also found
         * in output, or abort by returning NULL
         *
         * Input geometry was lwgeom_in
         */
        {
                        int loss;
                        GEOSGeometry *pi, *po, *pd;

                        /* TODO: handle some errors here...
                         * Lack of exceptions is annoying indeed,
                         * I'm getting old --strk;
                         */
                        pi = GEOSGeom_extractUniquePoints(gin);
                        po = GEOSGeom_extractUniquePoints(gout);
                        pd = GEOSDifference(pi, po); /* input points - output points */
                        GEOSGeom_destroy(pi);
                        GEOSGeom_destroy(po);
                        loss = pd && !GEOSisEmpty(pd);
                        GEOSGeom_destroy(pd);
                        if ( loss )
                        {
                                lwnotice("%s [%d] Vertices lost in LWGEOM_GEOS_makeValid", __FILE__, __LINE__);
                                /* return NULL */
                        }
        }
#endif /* PARANOIA_LEVEL > 1 */


        return gout;
}

Here is the call graph for this function:

Here is the caller graph for this function: