ST_Voronoi()

Datum ST_Voronoi

(

PG_FUNCTION_ARGS

)

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

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

References geometry_serialize(), gserialized_get_gbox_p(), lwgeom_free(), lwgeom_from_gserialized(), and lwgeom_voronoi_diagram().

Here is the call graph for this function: