optimistic_overlap()

Datum optimistic_overlap

(

PG_FUNCTION_ARGS

)

Definition at line 2712 of file lwgeom_functions_basic.c.

{
        GSERIALIZED *pg_geom1 = PG_GETARG_GSERIALIZED_P(0);
        GSERIALIZED *pg_geom2 = PG_GETARG_GSERIALIZED_P(1);
        double dist = PG_GETARG_FLOAT8(2);
        GBOX g1_bvol;
        double calc_dist;
        LWGEOM *geom1 = lwgeom_from_gserialized(pg_geom1);
        LWGEOM *geom2 = lwgeom_from_gserialized(pg_geom2);
        gserialized_error_if_srid_mismatch(pg_geom1, pg_geom2, __func__);
 
        if (geom1->type != POLYGONTYPE)
        {
                elog(ERROR, "optimistic_overlap: first arg isn't a polygon\n");
                PG_RETURN_NULL();
        }
 
        if (geom2->type != POLYGONTYPE && geom2->type != MULTIPOLYGONTYPE)
        {
                elog(ERROR, "optimistic_overlap: 2nd arg isn't a [multi-]polygon\n");
                PG_RETURN_NULL();
        }
 
        /*bbox check */
        gserialized_get_gbox_p(pg_geom1, &g1_bvol);
 
        g1_bvol.xmin = g1_bvol.xmin - dist;
        g1_bvol.ymin = g1_bvol.ymin - dist;
        g1_bvol.xmax = g1_bvol.xmax + dist;
        g1_bvol.ymax = g1_bvol.ymax + dist;
 
        if ((g1_bvol.xmin > geom2->bbox->xmax) || (g1_bvol.xmax < geom2->bbox->xmin) ||
            (g1_bvol.ymin > geom2->bbox->ymax) || (g1_bvol.ymax < geom2->bbox->ymin))
        {
                PG_RETURN_BOOL(false); /*bbox not overlap */
        }
 
        /*
         * compute distances
         * should be a fast calc if they actually do intersect
         */
        calc_dist =
            DatumGetFloat8(DirectFunctionCall2(ST_Distance, PointerGetDatum(pg_geom1), PointerGetDatum(pg_geom2)));
 
        PG_RETURN_BOOL(calc_dist < dist);
}

References LWGEOM::bbox, gserialized_error_if_srid_mismatch(), gserialized_get_gbox_p(), lwgeom_from_gserialized(), MULTIPOLYGONTYPE, POLYGONTYPE, ST_Distance(), LWGEOM::type, GBOX::xmax, GBOX::xmin, GBOX::ymax, and GBOX::ymin.

Here is the call graph for this function: