PostGIS  3.7.0dev-r@@SVN_REVISION@@

◆ LWGEOM_dwithin()

Datum LWGEOM_dwithin ( PG_FUNCTION_ARGS  )

Definition at line 731 of file lwgeom_functions_basic.c.

732 {
733  double mindist;
734  GSERIALIZED *geom1 = PG_GETARG_GSERIALIZED_P(0);
735  GSERIALIZED *geom2 = PG_GETARG_GSERIALIZED_P(1);
736  double tolerance = PG_GETARG_FLOAT8(2);
737  LWGEOM *lwgeom1 = lwgeom_from_gserialized(geom1);
738  LWGEOM *lwgeom2 = lwgeom_from_gserialized(geom2);
739 
740  if (tolerance < 0)
741  {
742  elog(ERROR, "Tolerance cannot be less than zero\n");
743  PG_RETURN_NULL();
744  }
745 
746  gserialized_error_if_srid_mismatch(geom1, geom2, __func__);
747 
748  if (lwgeom_is_empty(lwgeom1) || lwgeom_is_empty(lwgeom2))
749  {
750  PG_RETURN_BOOL(false);
751  }
752 
753  mindist = lwgeom_mindistance2d_tolerance(lwgeom1, lwgeom2, tolerance);
754 
755  PG_FREE_IF_COPY(geom1, 0);
756  PG_FREE_IF_COPY(geom2, 1);
757  /*empty geometries cases should be right handled since return from underlying
758  functions should be FLT_MAX which causes false as answer*/
759  PG_RETURN_BOOL(tolerance >= mindist);
760 }
void gserialized_error_if_srid_mismatch(const GSERIALIZED *g1, const GSERIALIZED *g2, const char *funcname)
Definition: gserialized.c:432
LWGEOM * lwgeom_from_gserialized(const GSERIALIZED *g)
Allocate a new LWGEOM from a GSERIALIZED.
Definition: gserialized.c:268
double lwgeom_mindistance2d_tolerance(const LWGEOM *lw1, const LWGEOM *lw2, double tolerance)
Function handling min distance calculations and dwithin calculations.
Definition: measures.c:210
static int lwgeom_is_empty(const LWGEOM *geom)
Return true or false depending on whether a geometry is an "empty" geometry (no vertices members)
Definition: lwinline.h:199

References gserialized_error_if_srid_mismatch(), lwgeom_from_gserialized(), lwgeom_is_empty(), and lwgeom_mindistance2d_tolerance().

Here is the call graph for this function: