PostGIS  3.4.0dev-r@@SVN_REVISION@@

◆ LWGEOM_dwithin()

Datum LWGEOM_dwithin ( PG_FUNCTION_ARGS  )

Definition at line 732 of file lwgeom_functions_basic.c.

733 {
734  double mindist;
735  GSERIALIZED *geom1 = PG_GETARG_GSERIALIZED_P(0);
736  GSERIALIZED *geom2 = PG_GETARG_GSERIALIZED_P(1);
737  double tolerance = PG_GETARG_FLOAT8(2);
738  LWGEOM *lwgeom1 = lwgeom_from_gserialized(geom1);
739  LWGEOM *lwgeom2 = lwgeom_from_gserialized(geom2);
740 
741  if (tolerance < 0)
742  {
743  elog(ERROR, "Tolerance cannot be less than zero\n");
744  PG_RETURN_NULL();
745  }
746 
747  gserialized_error_if_srid_mismatch(geom1, geom2, __func__);
748 
749  if (lwgeom_is_empty(lwgeom1) || lwgeom_is_empty(lwgeom2))
750  {
751  PG_RETURN_BOOL(false);
752  }
753 
754  mindist = lwgeom_mindistance2d_tolerance(lwgeom1, lwgeom2, tolerance);
755 
756  PG_FREE_IF_COPY(geom1, 0);
757  PG_FREE_IF_COPY(geom2, 1);
758  /*empty geometries cases should be right handled since return from underlying
759  functions should be FLT_MAX which causes false as answer*/
760  PG_RETURN_BOOL(tolerance >= mindist);
761 }
void gserialized_error_if_srid_mismatch(const GSERIALIZED *g1, const GSERIALIZED *g2, const char *funcname)
Definition: gserialized.c:403
LWGEOM * lwgeom_from_gserialized(const GSERIALIZED *g)
Allocate a new LWGEOM from a GSERIALIZED.
Definition: gserialized.c:239
double lwgeom_mindistance2d_tolerance(const LWGEOM *lw1, const LWGEOM *lw2, double tolerance)
Function handling min distance calculations and dwithin calculations.
Definition: measures.c:207
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:203

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

Here is the call graph for this function: