GEOS2LWGEOM()

LWGEOM * GEOS2LWGEOM	(	const GEOSGeometry *	geom,
		uint8_t	want3d
	)

Definition at line 193 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)
        {
                const GEOSCoordSequence* cs;
                POINTARRAY *pa, **ppaa;
                const GEOSGeometry* g;
                LWGEOM** geoms;
                uint32_t i, ngeoms;
 
        case GEOS_POINT:
                LWDEBUG(4, "lwgeom_from_geometry: it's a Point");
                cs = GEOSGeom_getCoordSeq(geom);
                if (GEOSisEmpty(geom)) return (LWGEOM*)lwpoint_construct_empty(SRID, want3d, 0);
                pa = ptarray_from_GEOSCoordSeq(cs, want3d);
                return (LWGEOM*)lwpoint_construct(SRID, NULL, pa);
 
        case GEOS_LINESTRING:
        case GEOS_LINEARRING:
                LWDEBUG(4, "lwgeom_from_geometry: it's a LineString or LinearRing");
                if (GEOSisEmpty(geom)) return (LWGEOM*)lwline_construct_empty(SRID, want3d, 0);
 
                cs = GEOSGeom_getCoordSeq(geom);
                pa = ptarray_from_GEOSCoordSeq(cs, want3d);
                return (LWGEOM*)lwline_construct(SRID, NULL, pa);
 
        case GEOS_POLYGON:
                LWDEBUG(4, "lwgeom_from_geometry: it's a Polygon");
                if (GEOSisEmpty(geom)) return (LWGEOM*)lwpoly_construct_empty(SRID, want3d, 0);
                ngeoms = GEOSGetNumInteriorRings(geom);
                ppaa = lwalloc(sizeof(POINTARRAY*) * (ngeoms + 1));
                g = GEOSGetExteriorRing(geom);
                cs = GEOSGeom_getCoordSeq(g);
                ppaa[0] = ptarray_from_GEOSCoordSeq(cs, want3d);
                for (i = 0; i < ngeoms; i++)
                {
                        g = GEOSGetInteriorRingN(geom, i);
                        cs = GEOSGeom_getCoordSeq(g);
                        ppaa[i + 1] = ptarray_from_GEOSCoordSeq(cs, want3d);
                }
                return (LWGEOM*)lwpoly_construct(SRID, NULL, ngeoms + 1, ppaa);
 
        case GEOS_MULTIPOINT:
        case GEOS_MULTILINESTRING:
        case GEOS_MULTIPOLYGON:
        case GEOS_GEOMETRYCOLLECTION:
                LWDEBUG(4, "lwgeom_from_geometry: it's a Collection or Multi");
 
                ngeoms = GEOSGetNumGeometries(geom);
                geoms = NULL;
                if (ngeoms)
                {
                        geoms = lwalloc(sizeof(LWGEOM*) * ngeoms);
                        for (i = 0; i < ngeoms; i++)
                        {
                                g = GEOSGetGeometryN(geom, i);
                                geoms[i] = GEOS2LWGEOM(g, want3d);
                        }
                }
                return (LWGEOM*)lwcollection_construct(type, SRID, NULL, ngeoms, geoms);
 
        default:
                lwerror("GEOS2LWGEOM: unknown geometry type: %d", type);
                return NULL;
        }
}

References lwalloc(), lwcollection_construct(), LWDEBUG, lwerror(), lwline_construct(), lwline_construct_empty(), lwpoint_construct(), lwpoint_construct_empty(), lwpoly_construct(), lwpoly_construct_empty(), ptarray_from_GEOSCoordSeq(), and SRID_UNKNOWN.

Here is the call graph for this function: