lwtin_from_geos()

LWTIN * lwtin_from_geos	(	const GEOSGeometry *	geom,
		uint8_t	want3d
	)

Definition at line 1908 of file liblwgeom/lwgeom_geos.c.

{
        int type = GEOSGeomTypeId(geom);
        int SRID = GEOSGetSRID(geom);
 
        /* GEOS's 0 is equivalent to our unknown as for SRID values */
        if (SRID == 0) SRID = SRID_UNKNOWN;
 
        if (want3d && !GEOSHasZ(geom))
        {
                LWDEBUG(3, "Geometry has no Z, won't provide one");
                want3d = 0;
        }
 
        switch (type)
        {
                LWTRIANGLE** geoms;
                uint32_t i, ngeoms;
        case GEOS_GEOMETRYCOLLECTION:
                LWDEBUG(4, "lwgeom_from_geometry: it's a Collection or Multi");
 
                ngeoms = GEOSGetNumGeometries(geom);
                geoms = NULL;
                if (ngeoms)
                {
                        geoms = lwalloc(ngeoms * sizeof *geoms);
                        if (!geoms)
                        {
                                lwerror("lwtin_from_geos: can't allocate geoms");
                                return NULL;
                        }
                        for (i = 0; i < ngeoms; i++)
                        {
                                const GEOSGeometry *poly, *ring;
                                const GEOSCoordSequence* cs;
                                POINTARRAY* pa;
 
                                poly = GEOSGetGeometryN(geom, i);
                                ring = GEOSGetExteriorRing(poly);
                                cs = GEOSGeom_getCoordSeq(ring);
                                pa = ptarray_from_GEOSCoordSeq(cs, want3d);
 
                                geoms[i] = lwtriangle_construct(SRID, NULL, pa);
                        }
                }
                return (LWTIN*)lwcollection_construct(TINTYPE, SRID, NULL, ngeoms, (LWGEOM**)geoms);
        case GEOS_POLYGON:
        case GEOS_MULTIPOINT:
        case GEOS_MULTILINESTRING:
        case GEOS_MULTIPOLYGON:
        case GEOS_LINESTRING:
        case GEOS_LINEARRING:
        case GEOS_POINT:
                lwerror("lwtin_from_geos: invalid geometry type for tin: %d", type);
                break;
 
        default:
                lwerror("GEOS2LWGEOM: unknown geometry type: %d", type);
                return NULL;
        }
 
        /* shouldn't get here */
        return NULL;
}

References lwalloc(), lwcollection_construct(), LWDEBUG, lwerror(), lwtriangle_construct(), ptarray_from_GEOSCoordSeq(), SRID_UNKNOWN, TINTYPE, and ovdump::type.

Referenced by lwgeom_delaunay_triangulation().

Here is the call graph for this function:

Here is the caller graph for this function: