PostGIS  3.2.2dev-r@@SVN_REVISION@@

◆ ST_ClipByBox2d()

Datum ST_ClipByBox2d ( PG_FUNCTION_ARGS  )

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

1500 {
1501  static const uint32_t geom_idx = 0;
1502  static const uint32_t box2d_idx = 1;
1504  LWGEOM *lwgeom1, *lwresult ;
1505  GBOX bbox1 = {0};
1506  GBOX *bbox2;
1507  uint8_t type;
1508  int32_t srid;
1509  lwflags_t flags;
1510 
1511  if (!gserialized_datum_get_internals_p(PG_GETARG_DATUM(geom_idx), &bbox1, &flags, &type, &srid))
1512  {
1513  /* empty clips to empty, no matter rect */
1514  PG_RETURN_DATUM(PG_GETARG_DATUM(geom_idx));
1515  }
1516 
1517  /* WARNING: this is really a BOX2DF, use only xmin and ymin fields */
1518  bbox2 = (GBOX *)PG_GETARG_POINTER(box2d_idx);
1519  bbox2->flags = 0;
1520 
1521  /* if bbox1 is covered by bbox2, return lwgeom1 */
1522  if (gbox_contains_2d(bbox2, &bbox1))
1523  {
1524  PG_RETURN_DATUM(PG_GETARG_DATUM(geom_idx));
1525  }
1526 
1527  /* If bbox1 outside of bbox2, return empty */
1528  if (!gbox_overlaps_2d(&bbox1, bbox2))
1529  {
1530  /* Get type and srid from datum */
1531  lwresult = lwgeom_construct_empty(type, srid, 0, 0);
1532  result = geometry_serialize(lwresult) ;
1533  lwgeom_free(lwresult) ;
1534  PG_RETURN_POINTER(result);
1535  }
1536 
1537  lwgeom1 = lwgeom_from_gserialized(PG_GETARG_GSERIALIZED_P(geom_idx));
1538  lwresult = lwgeom_clip_by_rect(lwgeom1, bbox2->xmin, bbox2->ymin,
1539  bbox2->xmax, bbox2->ymax);
1540 
1541  lwgeom_free(lwgeom1);
1542 
1543  if (!lwresult)
1544  PG_RETURN_NULL();
1545 
1546  result = geometry_serialize(lwresult) ;
1547  PG_RETURN_POINTER(result);
1548 }
char result[OUT_DOUBLE_BUFFER_SIZE]
Definition: cu_print.c:267
int gbox_overlaps_2d(const GBOX *g1, const GBOX *g2)
Return LW_TRUE if the GBOX overlaps on the 2d plane, LW_FALSE otherwise.
Definition: gbox.c:323
int gbox_contains_2d(const GBOX *g1, const GBOX *g2)
Return LW_TRUE if the first GBOX contains the second on the 2d plane, LW_FALSE otherwise.
Definition: gbox.c:339
LWGEOM * lwgeom_from_gserialized(const GSERIALIZED *g)
Allocate a new LWGEOM from a GSERIALIZED.
Definition: gserialized.c:239
int gserialized_datum_get_internals_p(Datum gsdatum, GBOX *gbox, lwflags_t *flags, uint8_t *type, int32_t *srid)
Peak into a GSERIALIZED datum to find its bounding box and some other metadata.
void lwgeom_free(LWGEOM *geom)
Definition: lwgeom.c:1138
uint16_t lwflags_t
Definition: liblwgeom.h:313
LWGEOM * lwgeom_clip_by_rect(const LWGEOM *geom1, double x0, double y0, double x1, double y1)
LWGEOM * lwgeom_construct_empty(uint8_t type, int32_t srid, char hasz, char hasm)
Definition: lwgeom.c:2111
type
Definition: ovdump.py:42
double ymax
Definition: liblwgeom.h:371
double xmax
Definition: liblwgeom.h:369
double ymin
Definition: liblwgeom.h:370
double xmin
Definition: liblwgeom.h:368
lwflags_t flags
Definition: liblwgeom.h:367

References GBOX::flags, gbox_contains_2d(), gbox_overlaps_2d(), gserialized_datum_get_internals_p(), lwgeom_clip_by_rect(), lwgeom_construct_empty(), lwgeom_free(), lwgeom_from_gserialized(), result, ovdump::type, GBOX::xmax, GBOX::xmin, GBOX::ymax, and GBOX::ymin.

Here is the call graph for this function: