LWGEOM_GEOS_makeValidMultiLine()

static GEOSGeometry* LWGEOM_GEOS_makeValidMultiLine ( const GEOSGeometry *  gin )

static

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

References lwalloc(), lwerror(), lwfree(), LWGEOM_GEOS_makeValid(), LWGEOM_GEOS_makeValidLine(), and lwrealloc().

Referenced by LWGEOM_GEOS_makeValid().

{
        GEOSGeometry** lines;
        GEOSGeometry** points;
        GEOSGeometry* mline_out=0;
        GEOSGeometry* mpoint_out=0;
        GEOSGeometry* gout=0;
        uint32_t nlines=0, nlines_alloc;
        uint32_t npoints=0;
        uint32_t ngeoms=0, nsubgeoms;
        uint32_t i, j;

        ngeoms = GEOSGetNumGeometries(gin);

        nlines_alloc = ngeoms;
        lines = lwalloc(sizeof(GEOSGeometry*)*nlines_alloc);
        points = lwalloc(sizeof(GEOSGeometry*)*ngeoms);

        for (i=0; i<ngeoms; ++i)
        {
                const GEOSGeometry* g = GEOSGetGeometryN(gin, i);
                GEOSGeometry* vg;
                vg = LWGEOM_GEOS_makeValidLine(g);
                /* Drop any invalid or empty geometry */
                if (!vg)
                        continue;
                if (GEOSisEmpty(vg))
                {
                        GEOSGeom_destroy(vg);
                        continue;
                }
                if ( GEOSGeomTypeId(vg) == GEOS_POINT )
                {
                        points[npoints++] = vg;
                }
                else if ( GEOSGeomTypeId(vg) == GEOS_LINESTRING )
                {
                        lines[nlines++] = vg;
                }
                else if ( GEOSGeomTypeId(vg) == GEOS_MULTILINESTRING )
                {
                        nsubgeoms=GEOSGetNumGeometries(vg);
                        nlines_alloc += nsubgeoms;
                        lines = lwrealloc(lines, sizeof(GEOSGeometry*)*nlines_alloc);
                        for (j=0; j<nsubgeoms; ++j)
                        {
                                const GEOSGeometry* gc = GEOSGetGeometryN(vg, j);
                                /* NOTE: ownership of the cloned geoms will be
                                 *       taken by final collection */
                                lines[nlines++] = GEOSGeom_clone(gc);
                        }
                }
                else
                {
                        /* NOTE: return from GEOSGeomType will leak
                         * but we really don't expect this to happen */
                        lwerror("unexpected geom type returned "
                                "by LWGEOM_GEOS_makeValid: %s",
                                GEOSGeomType(vg));
                }
        }

        if ( npoints )
        {
                if ( npoints > 1 )
                {
                        mpoint_out = GEOSGeom_createCollection(GEOS_MULTIPOINT,
                                                               points, npoints);
                }
                else
                {
                        mpoint_out = points[0];
                }
        }

        if ( nlines )
        {
                if ( nlines > 1 )
                {
                        mline_out = GEOSGeom_createCollection(
                                        GEOS_MULTILINESTRING, lines, nlines);
                }
                else
                {
                        mline_out = lines[0];
                }
        }

        lwfree(lines);

        if ( mline_out && mpoint_out )
        {
                points[0] = mline_out;
                points[1] = mpoint_out;
                gout = GEOSGeom_createCollection(GEOS_GEOMETRYCOLLECTION,
                                                 points, 2);
        }
        else if ( mline_out )
        {
                gout = mline_out;
        }
        else if ( mpoint_out )
        {
                gout = mpoint_out;
        }

        lwfree(points);

        return gout;
}

Here is the call graph for this function:

Here is the caller graph for this function: