PostGIS  2.3.8dev-r@@SVN_REVISION@@

◆ optimistic_overlap()

Datum optimistic_overlap ( PG_FUNCTION_ARGS  )

Definition at line 2460 of file lwgeom_functions_basic.c.

References LWGEOM::bbox, error_if_srid_mismatch(), FALSE, gserialized_get_gbox_p(), LWGEOM_affine(), lwgeom_from_gserialized(), LWGEOM_mindistance2d(), MULTIPOLYGONTYPE, PG_FUNCTION_INFO_V1(), POLYGONTYPE, LWGEOM::srid, LWGEOM::type, GBOX::xmax, GBOX::xmin, GBOX::ymax, and GBOX::ymin.

Referenced by LWGEOM_azimuth().

2461 {
2462  GSERIALIZED *pg_geom1 = PG_GETARG_GSERIALIZED_P(0);
2463  GSERIALIZED *pg_geom2 = PG_GETARG_GSERIALIZED_P(1);
2464  double dist = PG_GETARG_FLOAT8(2);
2465  GBOX g1_bvol;
2466  double calc_dist;
2467 
2468  LWGEOM *geom1 = lwgeom_from_gserialized(pg_geom1);
2469  LWGEOM *geom2 = lwgeom_from_gserialized(pg_geom2);
2470 
2471  error_if_srid_mismatch(geom1->srid, geom2->srid);
2472 
2473  if (geom1->type != POLYGONTYPE)
2474  {
2475  elog(ERROR,"optimistic_overlap: first arg isn't a polygon\n");
2476  PG_RETURN_NULL();
2477  }
2478 
2479  if (geom2->type != POLYGONTYPE && geom2->type != MULTIPOLYGONTYPE)
2480  {
2481  elog(ERROR,"optimistic_overlap: 2nd arg isn't a [multi-]polygon\n");
2482  PG_RETURN_NULL();
2483  }
2484 
2485  /*bbox check */
2486  gserialized_get_gbox_p(pg_geom1, &g1_bvol );
2487 
2488  g1_bvol.xmin = g1_bvol.xmin - dist;
2489  g1_bvol.ymin = g1_bvol.ymin - dist;
2490  g1_bvol.xmax = g1_bvol.xmax + dist;
2491  g1_bvol.ymax = g1_bvol.ymax + dist;
2492 
2493  if ( (g1_bvol.xmin > geom2->bbox->xmax) ||
2494  (g1_bvol.xmax < geom2->bbox->xmin) ||
2495  (g1_bvol.ymin > geom2->bbox->ymax) ||
2496  (g1_bvol.ymax < geom2->bbox->ymin) )
2497  {
2498  PG_RETURN_BOOL(FALSE); /*bbox not overlap */
2499  }
2500 
2501  /*
2502  * compute distances
2503  * should be a fast calc if they actually do intersect
2504  */
2505  calc_dist = DatumGetFloat8 ( DirectFunctionCall2(LWGEOM_mindistance2d, PointerGetDatum( pg_geom1 ), PointerGetDatum( pg_geom2 )));
2506 
2507  PG_RETURN_BOOL(calc_dist < dist);
2508 }
int gserialized_get_gbox_p(const GSERIALIZED *g, GBOX *box)
Read the bounding box off a serialization and calculate one if it is not already there.
Definition: g_serialized.c:398
GBOX * bbox
Definition: liblwgeom.h:397
LWGEOM * lwgeom_from_gserialized(const GSERIALIZED *g)
Allocate a new LWGEOM from a GSERIALIZED.
#define POLYGONTYPE
Definition: liblwgeom.h:86
double xmax
Definition: liblwgeom.h:292
void error_if_srid_mismatch(int srid1, int srid2)
Definition: lwutil.c:369
int32_t srid
Definition: liblwgeom.h:398
double ymin
Definition: liblwgeom.h:293
double xmin
Definition: liblwgeom.h:291
double ymax
Definition: liblwgeom.h:294
#define MULTIPOLYGONTYPE
Definition: liblwgeom.h:89
#define FALSE
Definition: dbfopen.c:168
Datum LWGEOM_mindistance2d(PG_FUNCTION_ARGS)
uint8_t type
Definition: liblwgeom.h:395
Here is the call graph for this function:
Here is the caller graph for this function: