SFCGAL2LWGEOM()

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

Definition at line 278 of file liblwgeom/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(), and SFCGALGeometry2POSTGIS().

Here is the call graph for this function:

Here is the caller graph for this function: