geography_segmentize()

Datum geography_segmentize

(

PG_FUNCTION_ARGS

)

Definition at line 1074 of file geography_measurement.c.

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.

Referenced by geography_azimuth().

{
        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);
}

Here is the call graph for this function:

Here is the caller graph for this function: