PostGIS  3.3.9dev-r@@SVN_REVISION@@

◆ hausdorffdistancedensify()

Datum hausdorffdistancedensify ( PG_FUNCTION_ARGS  )

Definition at line 221 of file postgis/lwgeom_geos.c.

222 {
223  GSERIALIZED *geom1;
224  GSERIALIZED *geom2;
225  GEOSGeometry *g1;
226  GEOSGeometry *g2;
227  double densifyFrac;
228  double result;
229  int retcode;
230 
231  geom1 = PG_GETARG_GSERIALIZED_P(0);
232  geom2 = PG_GETARG_GSERIALIZED_P(1);
233  densifyFrac = PG_GETARG_FLOAT8(2);
234 
235  if ( gserialized_is_empty(geom1) || gserialized_is_empty(geom2) )
236  PG_RETURN_NULL();
237 
238  initGEOS(lwpgnotice, lwgeom_geos_error);
239 
240  g1 = POSTGIS2GEOS(geom1);
241  if (!g1)
242  HANDLE_GEOS_ERROR("First argument geometry could not be converted to GEOS");
243 
244  g2 = POSTGIS2GEOS(geom2);
245  if (!g2)
246  {
247  GEOSGeom_destroy(g1);
248  HANDLE_GEOS_ERROR("Second argument geometry could not be converted to GEOS");
249  }
250 
251  retcode = GEOSHausdorffDistanceDensify(g1, g2, densifyFrac, &result);
252  GEOSGeom_destroy(g1);
253  GEOSGeom_destroy(g2);
254 
255  if (retcode == 0) HANDLE_GEOS_ERROR("GEOSHausdorffDistanceDensify");
256 
257  PG_FREE_IF_COPY(geom1, 0);
258  PG_FREE_IF_COPY(geom2, 1);
259 
260  PG_RETURN_FLOAT8(result);
261 }
char result[OUT_DOUBLE_BUFFER_SIZE]
Definition: cu_print.c:267
int gserialized_is_empty(const GSERIALIZED *g)
Check if a GSERIALIZED is empty without deserializing first.
Definition: gserialized.c:152
void lwgeom_geos_error(const char *fmt,...)
#define HANDLE_GEOS_ERROR(label)
GEOSGeometry * POSTGIS2GEOS(const GSERIALIZED *pglwgeom)

References gserialized_is_empty(), HANDLE_GEOS_ERROR, lwgeom_geos_error(), POSTGIS2GEOS(), and result.

Here is the call graph for this function: