SFCGAL2LWGEOM()

LWGEOM* SFCGAL2LWGEOM	(	const sfcgal_geometry_t *	geom,
		int	force3D,
		int32_t	srid
	)

Definition at line 323 of file lwgeom_sfcgal.c.

{
        uint32_t ngeoms, nshells, nsolids;
        uint32_t i, j, k;
        int want3d;
 
        assert(geom);
 
        want3d = force3D || sfcgal_geometry_is_3d(geom);
 
        switch (sfcgal_geometry_type_id(geom))
        {
        case SFCGAL_TYPE_POINT: {
                if (sfcgal_geometry_is_empty(geom))
                        return (LWGEOM *)lwpoint_construct_empty(srid, want3d, 0);
 
                POINTARRAY *pa = ptarray_from_SFCGAL(geom, want3d);
                return (LWGEOM *)lwpoint_construct(srid, NULL, pa);
        }
 
        case SFCGAL_TYPE_LINESTRING: {
                if (sfcgal_geometry_is_empty(geom))
                        return (LWGEOM *)lwline_construct_empty(srid, want3d, 0);
 
                POINTARRAY *pa = ptarray_from_SFCGAL(geom, want3d);
                return (LWGEOM *)lwline_construct(srid, NULL, pa);
        }
 
        case SFCGAL_TYPE_TRIANGLE: {
                if (sfcgal_geometry_is_empty(geom))
                        return (LWGEOM *)lwtriangle_construct_empty(srid, want3d, 0);
 
                POINTARRAY *pa = ptarray_from_SFCGAL(geom, want3d);
                return (LWGEOM *)lwtriangle_construct(srid, NULL, pa);
        }
 
        case SFCGAL_TYPE_POLYGON: {
                if (sfcgal_geometry_is_empty(geom))
                        return (LWGEOM *)lwpoly_construct_empty(srid, want3d, 0);
 
                uint32_t nrings = sfcgal_polygon_num_interior_rings(geom) + 1;
                POINTARRAY **pa = (POINTARRAY **)lwalloc(sizeof(POINTARRAY *) * nrings);
 
                pa[0] = ptarray_from_SFCGAL(sfcgal_polygon_exterior_ring(geom), want3d);
                for (i = 1; i < nrings; i++)
                        pa[i] = ptarray_from_SFCGAL(sfcgal_polygon_interior_ring_n(geom, i - 1), want3d);
 
                return (LWGEOM *)lwpoly_construct(srid, NULL, nrings, pa);
        }
 
        case SFCGAL_TYPE_MULTIPOINT:
        case SFCGAL_TYPE_MULTILINESTRING:
        case SFCGAL_TYPE_MULTIPOLYGON:
        case SFCGAL_TYPE_MULTISOLID:
        case SFCGAL_TYPE_GEOMETRYCOLLECTION: {
                ngeoms = sfcgal_geometry_collection_num_geometries(geom);
                LWGEOM **geoms = NULL;
                if (ngeoms)
                {
                        nsolids = 0;
                        geoms = (LWGEOM **)lwalloc(sizeof(LWGEOM *) * ngeoms);
                        for (i = 0; i < ngeoms; i++)
                        {
                                const sfcgal_geometry_t *g = sfcgal_geometry_collection_geometry_n(geom, i);
                                geoms[i] = SFCGAL2LWGEOM(g, 0, srid);
                                if (FLAGS_GET_SOLID(geoms[i]->flags))
                                        ++nsolids;
                        }
                        geoms = (LWGEOM **)lwrealloc(geoms, sizeof(LWGEOM *) * ngeoms);
                }
                LWGEOM *rgeom = (LWGEOM *)lwcollection_construct(
                    SFCGAL_type_to_lwgeom_type(sfcgal_geometry_type_id(geom)), srid, NULL, ngeoms, geoms);
                if (ngeoms)
                {
                        if (ngeoms == nsolids)
                                FLAGS_SET_SOLID(rgeom->flags, 1);
                        else if (nsolids)
                                lwnotice(
                                    "SFCGAL2LWGEOM: SOLID in heterogeneous collection will be treated as a POLYHEDRALSURFACETYPE");
                }
                return rgeom;
        }
 
#if 0
        case SFCGAL_TYPE_CIRCULARSTRING:
        case SFCGAL_TYPE_COMPOUNDCURVE:
        case SFCGAL_TYPE_CURVEPOLYGON:
        case SFCGAL_TYPE_MULTICURVE:
        case SFCGAL_TYPE_MULTISURFACE:
        case SFCGAL_TYPE_CURVE:
        case SFCGAL_TYPE_SURFACE:
 
        /* TODO curve types handling */
#endif
 
        case SFCGAL_TYPE_POLYHEDRALSURFACE: {
                ngeoms = sfcgal_polyhedral_surface_num_polygons(geom);
 
                LWGEOM **geoms = NULL;
                if (ngeoms)
                {
                        geoms = (LWGEOM **)lwalloc(sizeof(LWGEOM *) * ngeoms);
                        for (i = 0; i < ngeoms; i++)
                        {
                                const sfcgal_geometry_t *g = sfcgal_polyhedral_surface_polygon_n(geom, i);
                                geoms[i] = SFCGAL2LWGEOM(g, 0, srid);
                        }
                }
                return (LWGEOM *)lwcollection_construct(POLYHEDRALSURFACETYPE, srid, NULL, ngeoms, geoms);
        }
 
        /* Solid is map as a closed PolyhedralSurface (for now) */
        case SFCGAL_TYPE_SOLID: {
                nshells = sfcgal_solid_num_shells(geom);
 
                for (ngeoms = 0, i = 0; i < nshells; i++)
                        ngeoms += sfcgal_polyhedral_surface_num_polygons(sfcgal_solid_shell_n(geom, i));
 
                LWGEOM **geoms = 0;
                if (ngeoms)
                {
                        geoms = (LWGEOM **)lwalloc(sizeof(LWGEOM *) * ngeoms);
                        for (i = 0, k = 0; i < nshells; i++)
                        {
                                const sfcgal_geometry_t *shell = sfcgal_solid_shell_n(geom, i);
                                ngeoms = sfcgal_polyhedral_surface_num_polygons(shell);
 
                                for (j = 0; j < ngeoms; j++)
                                {
                                        const sfcgal_geometry_t *g = sfcgal_polyhedral_surface_polygon_n(shell, j);
                                        geoms[k] = SFCGAL2LWGEOM(g, 1, srid);
                                        k++;
                                }
                        }
                }
                LWGEOM *rgeom = (LWGEOM *)lwcollection_construct(POLYHEDRALSURFACETYPE, srid, NULL, ngeoms, geoms);
                if (ngeoms)
                        FLAGS_SET_SOLID(rgeom->flags, 1);
                return rgeom;
        }
 
        case SFCGAL_TYPE_TRIANGULATEDSURFACE: {
                ngeoms = sfcgal_triangulated_surface_num_triangles(geom);
                LWGEOM **geoms = NULL;
                if (ngeoms)
                {
                        geoms = (LWGEOM **)lwalloc(sizeof(LWGEOM *) * ngeoms);
                        for (i = 0; i < ngeoms; i++)
                        {
                                const sfcgal_geometry_t *g = sfcgal_triangulated_surface_triangle_n(geom, i);
                                geoms[i] = SFCGAL2LWGEOM(g, 0, srid);
                        }
                }
                return (LWGEOM *)lwcollection_construct(TINTYPE, srid, NULL, ngeoms, geoms);
        }
 
        default:
                lwerror("SFCGAL2LWGEOM: Unknown Type");
                return NULL;
        }
}

References LWGEOM::flags, FLAGS_GET_SOLID, FLAGS_SET_SOLID, lwalloc(), lwcollection_construct(), lwerror(), lwline_construct(), lwline_construct_empty(), lwnotice(), lwpoint_construct(), lwpoint_construct_empty(), lwpoly_construct(), lwpoly_construct_empty(), lwrealloc(), lwtriangle_construct(), lwtriangle_construct_empty(), POLYHEDRALSURFACETYPE, ptarray_from_SFCGAL(), SFCGAL_type_to_lwgeom_type(), and TINTYPE.

Referenced by lwgeom_sfcgal_noop().

Here is the call graph for this function:

Here is the caller graph for this function: