LWGEOM_dfullywithin()

Datum LWGEOM_dfullywithin

(

PG_FUNCTION_ARGS

)

Returns boolean true if the second argument is fully contained in a buffer of the first argument.

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

{
        GSERIALIZED *geom1 = PG_GETARG_GSERIALIZED_P(0);
        GSERIALIZED *geom2 = PG_GETARG_GSERIALIZED_P(1);
        LWGEOM *lwg1 = lwgeom_from_gserialized(geom1);
        LWGEOM *lwg2 = lwgeom_from_gserialized(geom2);
        double radius = PG_GETARG_FLOAT8(2);
        GEOSGeometry *buffer1 = NULL;
        GEOSGeometry *geos1 = NULL, *geos2 = NULL;
        char contained;
 
        if (radius < 0.0)
        {
                elog(ERROR, "Tolerance cannot be less than zero\n");
                PG_RETURN_NULL();
        }
 
        gserialized_error_if_srid_mismatch(geom1, geom2, __func__);
 
        if (lwgeom_is_empty(lwg1) || lwgeom_is_empty(lwg2))
                PG_RETURN_BOOL(false);
 
        if (!(lwgeom_isfinite(lwg1) && lwgeom_isfinite(lwg2)))
                PG_RETURN_BOOL(false);
 
        initGEOS(lwpgnotice, lwgeom_geos_error);
 
        geos1 = LWGEOM2GEOS(lwg1, true);
        geos2 = LWGEOM2GEOS(lwg2, true);
        lwgeom_free(lwg1);
        lwgeom_free(lwg2);
        if (!(geos1 && geos2))
                HANDLE_GEOS_ERROR("Geometry could not be converted to GEOS");
 
        buffer1 = GEOSBuffer(geos1, radius, 16);
        GEOSGeom_destroy(geos1);
        if (!(buffer1))
                HANDLE_GEOS_ERROR("Buffer operation failed");
 
        contained = GEOSCovers(buffer1, geos2);
        GEOSGeom_destroy(buffer1);
        GEOSGeom_destroy(geos2);
 
        if (contained == 2) HANDLE_GEOS_ERROR("GEOSContains");
 
        PG_FREE_IF_COPY(geom1, 0);
        PG_FREE_IF_COPY(geom2, 1);
 
        PG_RETURN_BOOL(contained == 1);
}

References gserialized_error_if_srid_mismatch(), HANDLE_GEOS_ERROR, LWGEOM2GEOS(), lwgeom_free(), lwgeom_from_gserialized(), lwgeom_geos_error(), lwgeom_is_empty(), and lwgeom_isfinite().

Here is the call graph for this function: