LWGEOM2GEOS()

GEOSGeometry* LWGEOM2GEOS	(	const LWGEOM *	lwgeom,
		int	autofix
	)

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

References COLLECTIONTYPE, POINTARRAY::flags, FLAGS_NDIMS, free(), LWCOLLECTION::geoms, getPoint_internal(), LINETYPE, LWDEBUGF, lwerror(), LWGEOM2GEOS(), lwgeom_free(), lwgeom_has_arc(), lwgeom_is_empty(), lwgeom_stroke(), lwtype_name(), malloc(), MULTILINETYPE, MULTIPOINTTYPE, MULTIPOLYGONTYPE, LWCOLLECTION::ngeoms, POINTARRAY::npoints, LWPOLY::nrings, LWPOINT::point, LWLINE::points, POINTTYPE, POLYGONTYPE, ptarray_addPoint(), ptarray_to_GEOSCoordSeq(), ptarray_to_GEOSLinearRing(), LWPOLY::rings, LWGEOM::srid, and LWGEOM::type.

Referenced by _lwt_AddFaceSplit(), _lwt_CheckEdgeCrossing(), _lwt_EdgeMotionArea(), _lwt_GetEqualEdge(), boundary(), isvalid(), LWGEOM2GEOS(), lwgeom_buildarea(), lwgeom_centroid(), lwgeom_clip_by_rect(), lwgeom_delaunay_triangulation(), lwgeom_difference(), lwgeom_extract_unique_endpoints(), lwgeom_geos_noop(), lwgeom_intersection(), lwgeom_is_simple(), lwgeom_linemerge(), lwgeom_make_valid(), lwgeom_node(), lwgeom_normalize(), lwgeom_offsetcurve(), lwgeom_sharedpaths(), lwgeom_snap(), lwgeom_symdifference(), lwgeom_unaryunion(), lwgeom_union(), LWGEOMARRAY2GEOS(), lwline_split_by_line(), lwpoly_split_by_line(), lwpoly_to_points(), lwt_AddPoint(), lwt_AddPolygon(), lwt_ChangeEdgeGeom(), POSTGIS2GEOS(), PrepGeomCacheBuilder(), rt_raster_compute_skewed_raster(), rt_raster_gdal_polygonize(), rt_raster_gdal_rasterize(), rt_raster_geos_spatial_relationship(), rt_raster_intersects(), and rt_raster_surface().

{
        GEOSCoordSeq sq;
        GEOSGeom g, shell;
        GEOSGeom *geoms = NULL;
        /*
        LWGEOM *tmp;
        */
        uint32_t ngeoms, i, j;
        int geostype;
#if LWDEBUG_LEVEL >= 4
        char *wkt;
#endif

        LWDEBUGF(4, "LWGEOM2GEOS got a %s", lwtype_name(lwgeom->type));

        if (lwgeom_has_arc(lwgeom))
        {
                LWGEOM *lwgeom_stroked = lwgeom_stroke(lwgeom, 32);
                GEOSGeometry *g = LWGEOM2GEOS(lwgeom_stroked, autofix);
                lwgeom_free(lwgeom_stroked);
                return g;
        }

        switch (lwgeom->type)
        {
                LWPOINT *lwp = NULL;
                LWPOLY *lwpoly = NULL;
                LWLINE *lwl = NULL;
                LWCOLLECTION *lwc = NULL;

        case POINTTYPE:
                lwp = (LWPOINT *)lwgeom;

                if ( lwgeom_is_empty(lwgeom) )
                {
                        g = GEOSGeom_createEmptyPolygon();
                }
                else
                {
                        sq = ptarray_to_GEOSCoordSeq(lwp->point, 0);
                        g = GEOSGeom_createPoint(sq);
                }
                if ( ! g )
                {
                        /* lwnotice("Exception in LWGEOM2GEOS"); */
                        return NULL;
                }
                break;
        case LINETYPE:
                lwl = (LWLINE *)lwgeom;
                /* TODO: if (autofix) */
                if ( lwl->points->npoints == 1 ) {
                        /* Duplicate point, to make geos-friendly */
                        lwl->points = ptarray_addPoint(lwl->points,
                                               getPoint_internal(lwl->points, 0),
                                               FLAGS_NDIMS(lwl->points->flags),
                                               lwl->points->npoints);
                }
                sq = ptarray_to_GEOSCoordSeq(lwl->points, 0);
                g = GEOSGeom_createLineString(sq);
                if ( ! g )
                {
                        /* lwnotice("Exception in LWGEOM2GEOS"); */
                        return NULL;
                }
                break;

        case POLYGONTYPE:
                lwpoly = (LWPOLY *)lwgeom;
                if ( lwgeom_is_empty(lwgeom) )
                {
                        g = GEOSGeom_createEmptyPolygon();
                }
                else
                {
                        shell = ptarray_to_GEOSLinearRing(lwpoly->rings[0], autofix);
                        if ( ! shell ) return NULL;
                        /*lwerror("LWGEOM2GEOS: exception during polygon shell conversion"); */
                        ngeoms = lwpoly->nrings-1;
                        if ( ngeoms > 0 )
                                geoms = malloc(sizeof(GEOSGeom)*ngeoms);

                        for (i=1; i<lwpoly->nrings; ++i)
                        {
                                geoms[i-1] = ptarray_to_GEOSLinearRing(lwpoly->rings[i], autofix);
                                if ( ! geoms[i-1] )
                                {
                                        --i;
                                        while (i) GEOSGeom_destroy(geoms[--i]);
                                        free(geoms);
                                        GEOSGeom_destroy(shell);
                                        return NULL;
                                }
                                /*lwerror("LWGEOM2GEOS: exception during polygon hole conversion"); */
                        }
                        g = GEOSGeom_createPolygon(shell, geoms, ngeoms);
                        if (geoms) free(geoms);
                }
                if ( ! g ) return NULL;
                break;
        case MULTIPOINTTYPE:
        case MULTILINETYPE:
        case MULTIPOLYGONTYPE:
        case COLLECTIONTYPE:
                if ( lwgeom->type == MULTIPOINTTYPE )
                        geostype = GEOS_MULTIPOINT;
                else if ( lwgeom->type == MULTILINETYPE )
                        geostype = GEOS_MULTILINESTRING;
                else if ( lwgeom->type == MULTIPOLYGONTYPE )
                        geostype = GEOS_MULTIPOLYGON;
                else
                        geostype = GEOS_GEOMETRYCOLLECTION;

                lwc = (LWCOLLECTION *)lwgeom;

                ngeoms = lwc->ngeoms;
                if ( ngeoms > 0 )
                        geoms = malloc(sizeof(GEOSGeom)*ngeoms);

                j = 0;
                for (i=0; i<ngeoms; ++i)
                {
                        GEOSGeometry* g;

                        if( lwgeom_is_empty(lwc->geoms[i]) )
                                continue;

                        g = LWGEOM2GEOS(lwc->geoms[i], 0);
                        if ( ! g )
                        {
                                while (j) GEOSGeom_destroy(geoms[--j]);
                                free(geoms);
                                return NULL;
                        }
                        geoms[j++] = g;
                }
                g = GEOSGeom_createCollection(geostype, geoms, j);
                if ( geoms ) free(geoms);
                if ( ! g ) return NULL;
                break;

        default:
                lwerror("Unknown geometry type: %d - %s", lwgeom->type, lwtype_name(lwgeom->type));
                return NULL;
        }

        GEOSSetSRID(g, lwgeom->srid);

#if LWDEBUG_LEVEL >= 4
        wkt = GEOSGeomToWKT(g);
        LWDEBUGF(4, "LWGEOM2GEOS: GEOSGeom: %s", wkt);
        free(wkt);
#endif

        return g;
}

Here is the call graph for this function:

Here is the caller graph for this function: