PostGIS  3.0.0dev-r@@SVN_REVISION@@

◆ LWGEOM_expand()

Datum LWGEOM_expand ( PG_FUNCTION_ARGS  )

Definition at line 1601 of file lwgeom_functions_basic.c.

References gbox_expand(), gbox_expand_xyzm(), dumpnode::geom, geometry_serialize(), LW_FAILURE, lwgeom_add_bbox(), lwgeom_calculate_gbox(), lwgeom_free(), lwgeom_from_gserialized(), lwgeom_get_srid(), lwgeom_has_m(), lwgeom_has_z(), lwgeom_is_empty(), lwgeom_set_srid(), LWGEOM_to_BOX(), lwpoly_as_lwgeom(), lwpoly_construct_rectangle(), GBOX::mmax, GBOX::mmin, PG_FUNCTION_INFO_V1(), GBOX::xmax, GBOX::xmin, GBOX::ymax, GBOX::ymin, GBOX::zmax, and GBOX::zmin.

Referenced by LWGEOM_makepoly().

1602 {
1603  GSERIALIZED *geom = PG_GETARG_GSERIALIZED_P(0);
1604  LWGEOM *lwgeom = lwgeom_from_gserialized(geom);
1605  int srid = lwgeom_get_srid(lwgeom);
1606  LWPOLY *poly;
1607  GSERIALIZED *result;
1608  GBOX gbox;
1609 
1610  POSTGIS_DEBUG(2, "LWGEOM_expand called.");
1611 
1612  /* Can't expand an empty */
1613  if ( lwgeom_is_empty(lwgeom) )
1614  {
1615  lwgeom_free(lwgeom);
1616  PG_RETURN_POINTER(geom);
1617  }
1618 
1619  /* Can't expand something with no gbox! */
1620  if ( LW_FAILURE == lwgeom_calculate_gbox(lwgeom, &gbox) )
1621  {
1622  lwgeom_free(lwgeom);
1623  PG_RETURN_POINTER(geom);
1624  }
1625 
1626  if (PG_NARGS() == 2)
1627  {
1628  /* Expand the box the same amount in all directions */
1629  double d = PG_GETARG_FLOAT8(1);
1630  gbox_expand(&gbox, d);
1631  }
1632  else
1633  {
1634  double dx = PG_GETARG_FLOAT8(1);
1635  double dy = PG_GETARG_FLOAT8(2);
1636  double dz = PG_GETARG_FLOAT8(3);
1637  double dm = PG_GETARG_FLOAT8(4);
1638 
1639  gbox_expand_xyzm(&gbox, dx, dy, dz, dm);
1640  }
1641 
1642  {
1643  POINT4D p1 = { gbox.xmin, gbox.ymin, gbox.zmin, gbox.mmin };
1644  POINT4D p2 = { gbox.xmin, gbox.ymax, gbox.zmin, gbox.mmin };
1645  POINT4D p3 = { gbox.xmax, gbox.ymax, gbox.zmax, gbox.mmax };
1646  POINT4D p4 = { gbox.xmax, gbox.ymin, gbox.zmax, gbox.mmax };
1647 
1648  poly = lwpoly_construct_rectangle(lwgeom_has_z(lwgeom), lwgeom_has_m(lwgeom), &p1, &p2, &p3, &p4);
1649  }
1650 
1652  lwgeom_set_srid(lwpoly_as_lwgeom(poly), srid);
1653 
1654  /* Construct GSERIALIZED */
1655  result = geometry_serialize(lwpoly_as_lwgeom(poly));
1656 
1658  lwgeom_free(lwgeom);
1659  PG_FREE_IF_COPY(geom, 0);
1660 
1661  PG_RETURN_POINTER(result);
1662 }
void gbox_expand_xyzm(GBOX *g, double dx, double dy, double dz, double dm)
Move the box minimums down and the maximums up by the distances provided.
Definition: g_box.c:122
LWGEOM * lwgeom_from_gserialized(const GSERIALIZED *g)
Allocate a new LWGEOM from a GSERIALIZED.
void gbox_expand(GBOX *g, double d)
Move the box minimums down and the maximums up by the distance provided.
Definition: g_box.c:104
double xmax
Definition: liblwgeom.h:295
void lwgeom_free(LWGEOM *geom)
Definition: lwgeom.c:1128
int32_t lwgeom_get_srid(const LWGEOM *geom)
Return SRID number.
Definition: lwgeom.c:907
int lwgeom_has_z(const LWGEOM *geom)
Return LW_TRUE if geometry has Z ordinates.
Definition: lwgeom.c:914
LWGEOM * lwpoly_as_lwgeom(const LWPOLY *obj)
Definition: lwgeom.c:311
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:172
#define LW_FAILURE
Definition: liblwgeom.h:78
int lwgeom_calculate_gbox(const LWGEOM *lwgeom, GBOX *gbox)
Calculate bounding box of a geometry, automatically taking into account whether it is cartesian or ge...
Definition: lwgeom.c:737
double zmax
Definition: liblwgeom.h:299
double ymin
Definition: liblwgeom.h:296
LWGEOM * geom
double xmin
Definition: liblwgeom.h:294
double ymax
Definition: liblwgeom.h:297
GSERIALIZED * geometry_serialize(LWGEOM *lwgeom)
double mmin
Definition: liblwgeom.h:300
double zmin
Definition: liblwgeom.h:298
void lwgeom_add_bbox(LWGEOM *lwgeom)
Compute a bbox if not already computed.
Definition: lwgeom.c:677
double mmax
Definition: liblwgeom.h:301
void lwgeom_set_srid(LWGEOM *geom, int srid)
Set the SRID on an LWGEOM For collections, only the parent gets an SRID, all the children get SRID_UN...
int lwgeom_has_m(const LWGEOM *geom)
Return LW_TRUE if geometry has M ordinates.
Definition: lwgeom.c:921
LWPOLY * lwpoly_construct_rectangle(char hasz, char hasm, POINT4D *p1, POINT4D *p2, POINT4D *p3, POINT4D *p4)
Definition: lwpoly.c:80
Here is the call graph for this function:
Here is the caller graph for this function: