ST_Voronoi()

Datum ST_Voronoi

(

PG_FUNCTION_ARGS

)

Definition at line 3604 of file postgis/lwgeom_geos.c.

References geometry_serialize(), gserialized_get_gbox_p(), lwgeom_free(), lwgeom_from_gserialized(), lwgeom_voronoi_diagram(), PG_FUNCTION_INFO_V1(), POSTGIS_GEOS_VERSION, and ST_MinimumClearance().

Referenced by ST_Node().

{
#if POSTGIS_GEOS_VERSION < 35
        lwpgerror("The GEOS version this PostGIS binary "
                "was compiled against (%d) doesn't support "
                "'ST_Voronoi' function (3.5.0+ required)",
                POSTGIS_GEOS_VERSION);
        PG_RETURN_NULL();
#else /* POSTGIS_GEOS_VERSION >= 35 */
        GSERIALIZED* input;
        GSERIALIZED* clip;
        GSERIALIZED* result;
        LWGEOM* lwgeom_input;
        LWGEOM* lwgeom_result;
        double tolerance;
        GBOX clip_envelope;
        int custom_clip_envelope;
        int return_polygons;

        /* Return NULL on NULL geometry */
        if (PG_ARGISNULL(0))
                PG_RETURN_NULL();

        /* Read our tolerance value */
        if (PG_ARGISNULL(2))
        {
                lwpgerror("Tolerance must be a positive number.");
                PG_RETURN_NULL();
        }

        tolerance = PG_GETARG_FLOAT8(2);

        if (tolerance < 0)
        {
                lwpgerror("Tolerance must be a positive number.");
                PG_RETURN_NULL();
        }

        /* Are we returning lines or polygons? */
        if (PG_ARGISNULL(3))
        {
                lwpgerror("return_polygons must be true or false.");
                PG_RETURN_NULL();
        }
        return_polygons = PG_GETARG_BOOL(3);

        /* Read our clipping envelope, if applicable. */
        custom_clip_envelope = !PG_ARGISNULL(1);
        if (custom_clip_envelope) {
                clip = PG_GETARG_GSERIALIZED_P(1);
                if (!gserialized_get_gbox_p(clip, &clip_envelope))
                {
                        lwpgerror("Could not determine envelope of clipping geometry.");
                        PG_FREE_IF_COPY(clip, 1);
                        PG_RETURN_NULL();
                }
                PG_FREE_IF_COPY(clip, 1);
        }

        /* Read our input geometry */
        input = PG_GETARG_GSERIALIZED_P(0);

        lwgeom_input = lwgeom_from_gserialized(input);

        if(!lwgeom_input)
        {
                lwpgerror("Could not read input geometry.");
                PG_FREE_IF_COPY(input, 0);
                PG_RETURN_NULL();
        }

        lwgeom_result = lwgeom_voronoi_diagram(lwgeom_input, custom_clip_envelope ? &clip_envelope : NULL, tolerance, !return_polygons);
        lwgeom_free(lwgeom_input);

        if (!lwgeom_result)
        {
                lwpgerror("Error computing Voronoi diagram.");
                PG_FREE_IF_COPY(input, 0);
                PG_RETURN_NULL();
        }

        result = geometry_serialize(lwgeom_result);
        lwgeom_free(lwgeom_result);

        PG_FREE_IF_COPY(input, 0);
        PG_RETURN_POINTER(result);

#endif /* POSTGIS_GEOS_VERSION >= 35 */
}

Here is the call graph for this function:

Here is the caller graph for this function: