2041 uint32_t worldTileSize;
2042 double tileGeoSizeX, tileGeoSizeY;
2043 double boundsWidth, boundsHeight;
2044 double x1, y1, x2, y2;
2055 POSTGIS_DEBUG(2,
"ST_TileEnvelope called");
2057 zoom = PG_GETARG_INT32(0);
2058 x = PG_GETARG_INT32(1);
2059 y = PG_GETARG_INT32(2);
2061 bounds = PG_GETARG_GSERIALIZED_P(3);
2070 elog(ERROR,
"%s: Unable to compute bbox", __func__);
2075 margin = PG_NARGS() < 4 ? 0 : PG_GETARG_FLOAT8(4);
2078 elog(ERROR,
"%s: Margin must not be less than -50%%, margin=%f", __func__, margin);
2080 boundsWidth = bbox.
xmax - bbox.
xmin;
2081 boundsHeight = bbox.
ymax - bbox.
ymin;
2082 if (boundsWidth <= 0 || boundsHeight <= 0)
2083 elog(ERROR,
"%s: Geometric bounds are too small", __func__);
2085 if (zoom < 0 || zoom >= 32)
2086 elog(ERROR,
"%s: Invalid tile zoom value, %d", __func__, zoom);
2088 zoomu = (uint32_t)zoom;
2089 worldTileSize = 0x01u << (zoomu > 31 ? 31 : zoomu);
2091 if (
x < 0 || (uint32_t)
x >= worldTileSize)
2092 elog(ERROR,
"%s: Invalid tile x value, %d", __func__,
x);
2093 if (
y < 0 || (uint32_t)
y >= worldTileSize)
2094 elog(ERROR,
"%s: Invalid tile y value, %d", __func__,
y);
2096 tileGeoSizeX = boundsWidth / worldTileSize;
2097 tileGeoSizeY = boundsHeight / worldTileSize;
2104 if ((1 + margin * 2) > worldTileSize)
2111 x1 = bbox.
xmin + tileGeoSizeX * (
x - margin);
2112 x2 = bbox.
xmin + tileGeoSizeX * (
x + 1 + margin);
2115 y1 = bbox.
ymax - tileGeoSizeY * (
y + 1 + margin);
2116 y2 = bbox.
ymax - tileGeoSizeY * (
y - margin);
2119 if (y1 < bbox.
ymin) y1 = bbox.
ymin;
2120 if (y2 > bbox.
ymax) y2 = bbox.
ymax;
2126 srid, x1, y1, x2, y2))));
LWGEOM * lwgeom_from_gserialized(const GSERIALIZED *g)
Allocate a new LWGEOM from a GSERIALIZED.
void lwgeom_free(LWGEOM *geom)
LWGEOM * lwpoly_as_lwgeom(const LWPOLY *obj)
LWPOLY * lwpoly_construct_envelope(int32_t srid, double x1, double y1, double x2, double y2)
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...