GEOS2LWGEOM()

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

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

References GEOS2LWGEOM(), lwalloc(), lwcollection_construct(), LWDEBUG, lwerror(), lwline_construct(), lwline_construct_empty(), lwpoint_construct(), lwpoint_construct_empty(), lwpoly_construct(), lwpoly_construct_empty(), ptarray_from_GEOSCoordSeq(), ptarray_to_GEOSCoordSeq(), SRID_UNKNOWN, and ovdump::type.

Referenced by convexhull(), GEOS2LWGEOM(), GEOS2POSTGIS(), GEOSARRAY2LWGEOM(), isvaliddetail(), lwgeom_buildarea(), lwgeom_centroid(), lwgeom_clip_by_rect(), lwgeom_delaunay_triangulation(), lwgeom_difference(), lwgeom_extract_unique_endpoints(), LWGEOM_GEOS_buildArea(), LWGEOM_GEOS_makeValid(), LWGEOM_GEOS_makeValidPolygon(), LWGEOM_GEOS_nodeLines(), lwgeom_geos_noop(), lwgeom_intersection(), lwgeom_linemerge(), lwgeom_make_valid(), lwgeom_node(), lwgeom_normalize(), lwgeom_offsetcurve(), lwgeom_sharedpaths(), lwgeom_snap(), lwgeom_symdifference(), lwgeom_unaryunion(), lwgeom_union(), lwgeom_voronoi_diagram(), lwline_split_by_line(), lwpoly_split_by_line(), and rt_raster_surface().

{
        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;
                }
        }

/*
        if ( GEOSisEmpty(geom) )
        {
                return (LWGEOM*)lwcollection_construct_empty(COLLECTIONTYPE, SRID, 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;

        }

}

Here is the call graph for this function:

Here is the caller graph for this function: