PostGIS  3.1.6dev-r@@SVN_REVISION@@

◆ optimistic_overlap()

Datum optimistic_overlap ( PG_FUNCTION_ARGS  )

Definition at line 2600 of file lwgeom_functions_basic.c.

2601 {
2602  GSERIALIZED *pg_geom1 = PG_GETARG_GSERIALIZED_P(0);
2603  GSERIALIZED *pg_geom2 = PG_GETARG_GSERIALIZED_P(1);
2604  double dist = PG_GETARG_FLOAT8(2);
2605  GBOX g1_bvol;
2606  double calc_dist;
2607  LWGEOM *geom1 = lwgeom_from_gserialized(pg_geom1);
2608  LWGEOM *geom2 = lwgeom_from_gserialized(pg_geom2);
2609  gserialized_error_if_srid_mismatch(pg_geom1, pg_geom2, __func__);
2610 
2611  if (geom1->type != POLYGONTYPE)
2612  {
2613  elog(ERROR, "optimistic_overlap: first arg isn't a polygon\n");
2614  PG_RETURN_NULL();
2615  }
2616 
2617  if (geom2->type != POLYGONTYPE && geom2->type != MULTIPOLYGONTYPE)
2618  {
2619  elog(ERROR, "optimistic_overlap: 2nd arg isn't a [multi-]polygon\n");
2620  PG_RETURN_NULL();
2621  }
2622 
2623  /*bbox check */
2624  gserialized_get_gbox_p(pg_geom1, &g1_bvol);
2625 
2626  g1_bvol.xmin = g1_bvol.xmin - dist;
2627  g1_bvol.ymin = g1_bvol.ymin - dist;
2628  g1_bvol.xmax = g1_bvol.xmax + dist;
2629  g1_bvol.ymax = g1_bvol.ymax + dist;
2630 
2631  if ((g1_bvol.xmin > geom2->bbox->xmax) || (g1_bvol.xmax < geom2->bbox->xmin) ||
2632  (g1_bvol.ymin > geom2->bbox->ymax) || (g1_bvol.ymax < geom2->bbox->ymin))
2633  {
2634  PG_RETURN_BOOL(false); /*bbox not overlap */
2635  }
2636 
2637  /*
2638  * compute distances
2639  * should be a fast calc if they actually do intersect
2640  */
2641  calc_dist =
2642  DatumGetFloat8(DirectFunctionCall2(ST_Distance, PointerGetDatum(pg_geom1), PointerGetDatum(pg_geom2)));
2643 
2644  PG_RETURN_BOOL(calc_dist < dist);
2645 }
void gserialized_error_if_srid_mismatch(const GSERIALIZED *g1, const GSERIALIZED *g2, const char *funcname)
Definition: gserialized.c:404
int gserialized_get_gbox_p(const GSERIALIZED *g, GBOX *gbox)
Read the box from the GSERIALIZED or calculate it if necessary.
Definition: gserialized.c:65
LWGEOM * lwgeom_from_gserialized(const GSERIALIZED *g)
Allocate a new LWGEOM from a GSERIALIZED.
Definition: gserialized.c:239
#define MULTIPOLYGONTYPE
Definition: liblwgeom.h:121
#define POLYGONTYPE
Definition: liblwgeom.h:118
Datum ST_Distance(PG_FUNCTION_ARGS)
double ymax
Definition: liblwgeom.h:371
double xmax
Definition: liblwgeom.h:369
double ymin
Definition: liblwgeom.h:370
double xmin
Definition: liblwgeom.h:368
uint8_t type
Definition: liblwgeom.h:476
GBOX * bbox
Definition: liblwgeom.h:472

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: