PostGIS  3.4.0dev-r@@SVN_REVISION@@

◆ LWGEOM_FilterByM()

Datum LWGEOM_FilterByM ( PG_FUNCTION_ARGS  )

Definition at line 3358 of file lwgeom_functions_basic.c.

3359 {
3360  GSERIALIZED *geom_in;
3361  GSERIALIZED *geom_out;
3362  LWGEOM *lwgeom_in;
3363  LWGEOM *lwgeom_out;
3364  double min, max;
3365  int returnm;
3366  int hasm;
3367 
3368  if (PG_NARGS() > 0 && !PG_ARGISNULL(0))
3369  {
3370  geom_in = PG_GETARG_GSERIALIZED_P(0);
3371  }
3372  else
3373  {
3374  PG_RETURN_NULL();
3375  }
3376 
3377  if (PG_NARGS() > 1 && !PG_ARGISNULL(1))
3378  min = PG_GETARG_FLOAT8(1);
3379  else
3380  {
3381  min = DBL_MIN;
3382  }
3383  if (PG_NARGS() > 2 && !PG_ARGISNULL(2))
3384  max = PG_GETARG_FLOAT8(2);
3385  else
3386  {
3387  max = DBL_MAX;
3388  }
3389  if (PG_NARGS() > 3 && !PG_ARGISNULL(3) && PG_GETARG_BOOL(3))
3390  returnm = 1;
3391  else
3392  {
3393  returnm = 0;
3394  }
3395 
3396  if (min > max)
3397  {
3398  elog(ERROR, "Min-value cannot be larger than Max value\n");
3399  PG_RETURN_NULL();
3400  }
3401 
3402  lwgeom_in = lwgeom_from_gserialized(geom_in);
3403 
3404  hasm = lwgeom_has_m(lwgeom_in);
3405 
3406  if (!hasm)
3407  {
3408  elog(NOTICE, "No M-value, No vertex removed\n");
3409  PG_RETURN_POINTER(geom_in);
3410  }
3411 
3412  lwgeom_out = lwgeom_filter_m(lwgeom_in, min, max, returnm);
3413 
3414  geom_out = geometry_serialize(lwgeom_out);
3415  lwgeom_free(lwgeom_out);
3416  PG_RETURN_POINTER(geom_out);
3417 }
LWGEOM * lwgeom_from_gserialized(const GSERIALIZED *g)
Allocate a new LWGEOM from a GSERIALIZED.
Definition: gserialized.c:239
void lwgeom_free(LWGEOM *geom)
Definition: lwgeom.c:1155
LWGEOM * lwgeom_filter_m(LWGEOM *geom, double min, double max, int returnm)
Definition: lwmval.c:221
int lwgeom_has_m(const LWGEOM *geom)
Return LW_TRUE if geometry has M ordinates.
Definition: lwgeom.c:941

References lwgeom_filter_m(), lwgeom_free(), lwgeom_from_gserialized(), and lwgeom_has_m().

Here is the call graph for this function: