geography_segmentize()

Datum geography_segmentize

(

PG_FUNCTION_ARGS

)

Definition at line 1118 of file geography_measurement.c.

{
        LWGEOM *lwgeom1 = NULL;
        LWGEOM *lwgeom2 = NULL;
        GSERIALIZED *g1 = NULL;
        GSERIALIZED *g2 = NULL;
        double max_seg_length;
        uint32_t type1;
 
        /* Get our geometry object loaded into memory. */
        g1 = PG_GETARG_GSERIALIZED_P(0);
        type1 = gserialized_get_type(g1);
 
        /* Convert max_seg_length from metric units to radians */
        max_seg_length = PG_GETARG_FLOAT8(1) / WGS84_RADIUS;
 
        /* We can't densify points or points, reflect them back */
        if ( type1 == POINTTYPE || type1 == MULTIPOINTTYPE || gserialized_is_empty(g1) )
                PG_RETURN_POINTER(g1);
 
        /* Deserialize */
        lwgeom1 = lwgeom_from_gserialized(g1);
 
        /* Calculate the densified geometry */
        lwgeom2 = lwgeom_segmentize_sphere(lwgeom1, max_seg_length);
 
        /*
        ** Set the geodetic flag so subsequent
        ** functions do the right thing.
        */
        lwgeom_set_geodetic(lwgeom2, true);
 
        /* Recalculate the boxes after re-setting the geodetic bit */
        lwgeom_drop_bbox(lwgeom2);
 
        /* We are trusting geography_serialize will add a box if needed */
        g2 = geography_serialize(lwgeom2);
 
        /* Clean up */
        lwgeom_free(lwgeom1);
        lwgeom_free(lwgeom2);
        PG_FREE_IF_COPY(g1, 0);
 
        PG_RETURN_POINTER(g2);
}

References gserialized_get_type(), gserialized_is_empty(), lwgeom_drop_bbox(), lwgeom_free(), lwgeom_from_gserialized(), lwgeom_segmentize_sphere(), lwgeom_set_geodetic(), MULTIPOINTTYPE, POINTTYPE, and WGS84_RADIUS.

Here is the call graph for this function: