lwtin_from_geos()

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

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

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

Referenced by lwgeom_delaunay_triangulation().

                                                             {
        int type = GEOSGeomTypeId(geom);
        int hasZ;
        int SRID = GEOSGetSRID(geom);

        /* GEOS's 0 is equivalent to our unknown as for SRID values */
        if ( SRID == 0 ) SRID = SRID_UNKNOWN;

        if ( want3d ) {
                hasZ = GEOSHasZ(geom);
                if ( ! hasZ ) {
                        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;
}

Here is the call graph for this function:

Here is the caller graph for this function: