lwgeom_segmentize_sphere()

LWGEOM* lwgeom_segmentize_sphere	(	const LWGEOM *	lwg_in,
		double	max_seg_length
	)

Create a new, densified geometry where no segment is longer than max_seg_length.

Derive a new geometry with vertices added to ensure no vertex is more than max_seg_length (in radians) from any other vertex.

Input geometry is not altered, output geometry must be freed by caller.

Parameters

lwg_in	= input geometry
max_seg_length	= maximum segment length in radians

Definition at line 1749 of file lwgeodetic.c.

{
        POINTARRAY *pa_out;
        LWLINE *lwline;
        LWPOLY *lwpoly_in, *lwpoly_out;
        LWCOLLECTION *lwcol_in, *lwcol_out;
        uint32_t i;
 
        /* Reflect NULL */
        if ( ! lwg_in )
                return NULL;
 
        /* Clone empty */
        if ( lwgeom_is_empty(lwg_in) )
                return lwgeom_clone(lwg_in);
 
        switch (lwg_in->type)
        {
        case MULTIPOINTTYPE:
        case POINTTYPE:
                return lwgeom_clone_deep(lwg_in);
                break;
        case LINETYPE:
                lwline = lwgeom_as_lwline(lwg_in);
                pa_out = ptarray_segmentize_sphere(lwline->points, max_seg_length);
                return lwline_as_lwgeom(lwline_construct(lwg_in->srid, NULL, pa_out));
                break;
        case POLYGONTYPE:
                lwpoly_in = lwgeom_as_lwpoly(lwg_in);
                lwpoly_out = lwpoly_construct_empty(lwg_in->srid, lwgeom_has_z(lwg_in), lwgeom_has_m(lwg_in));
                for ( i = 0; i < lwpoly_in->nrings; i++ )
                {
                        pa_out = ptarray_segmentize_sphere(lwpoly_in->rings[i], max_seg_length);
                        lwpoly_add_ring(lwpoly_out, pa_out);
                }
                return lwpoly_as_lwgeom(lwpoly_out);
                break;
        case MULTILINETYPE:
        case MULTIPOLYGONTYPE:
        case COLLECTIONTYPE:
                lwcol_in = lwgeom_as_lwcollection(lwg_in);
                lwcol_out = lwcollection_construct_empty(lwg_in->type, lwg_in->srid, lwgeom_has_z(lwg_in), lwgeom_has_m(lwg_in));
                for ( i = 0; i < lwcol_in->ngeoms; i++ )
                {
                        lwcollection_add_lwgeom(lwcol_out, lwgeom_segmentize_sphere(lwcol_in->geoms[i], max_seg_length));
                }
                return lwcollection_as_lwgeom(lwcol_out);
                break;
        default:
                lwerror("lwgeom_segmentize_sphere: unsupported input geometry type: %d - %s",
                        lwg_in->type, lwtype_name(lwg_in->type));
                break;
        }
 
        lwerror("lwgeom_segmentize_sphere got to the end of the function, should not happen");
        return NULL;
}

References COLLECTIONTYPE, LWCOLLECTION::geoms, LINETYPE, lwcollection_add_lwgeom(), lwcollection_as_lwgeom(), lwcollection_construct_empty(), lwerror(), lwgeom_as_lwcollection(), lwgeom_as_lwline(), lwgeom_as_lwpoly(), lwgeom_clone(), lwgeom_clone_deep(), lwgeom_has_m(), lwgeom_has_z(), lwgeom_is_empty(), lwgeom_segmentize_sphere(), lwline_as_lwgeom(), lwline_construct(), lwpoly_add_ring(), lwpoly_as_lwgeom(), lwpoly_construct_empty(), lwtype_name(), MULTILINETYPE, MULTIPOINTTYPE, MULTIPOLYGONTYPE, LWCOLLECTION::ngeoms, LWPOLY::nrings, LWLINE::points, POINTTYPE, POLYGONTYPE, ptarray_segmentize_sphere(), LWPOLY::rings, LWGEOM::srid, and LWGEOM::type.

Referenced by geography_segmentize(), lwgeom_segmentize_sphere(), test_geos_subdivide(), and test_lwgeom_segmentize_sphere().

Here is the call graph for this function:

Here is the caller graph for this function: