ST_TileEnvelope()

Datum ST_TileEnvelope

(

PG_FUNCTION_ARGS

)

Definition at line 2039 of file lwgeom_functions_basic.c.

{
        GSERIALIZED *bounds;
        uint32_t zoomu;
        int32_t x, y, zoom;
        uint32_t worldTileSize;
        double tileGeoSizeX, tileGeoSizeY;
        double boundsWidth, boundsHeight;
        double x1, y1, x2, y2;
        /* This is broken, since 3857 doesn't mean "web mercator", it means
           the contents of the row in spatial_ref_sys with srid = 3857.
           For practical purposes this will work, but in good implementation
           we should de-reference in spatial ref sys to confirm that the
           srid of the object is EPSG:3857. */
        int32_t srid;
        GBOX bbox;
        LWGEOM *g = NULL;
 
        POSTGIS_DEBUG(2, "ST_TileEnvelope called");
 
        zoom = PG_GETARG_INT32(0);
        x = PG_GETARG_INT32(1);
        y = PG_GETARG_INT32(2);
 
        bounds = PG_GETARG_GSERIALIZED_P(3);
        /*
         * We deserialize the geometry and recalculate the bounding box here to get
         * 64b floating point precision. The serialized bbox has 32b float is not
         * precise enough with big numbers such as the ones used in the default
         * parameters, e.g: -20037508.3427892 is transformed into -20037510
         */
        g = lwgeom_from_gserialized(bounds);
        if (lwgeom_calculate_gbox(g, &bbox) != LW_SUCCESS)
                elog(ERROR, "%s: Unable to compute bbox", __func__);
        srid = g->srid;
        lwgeom_free(g);
 
        boundsWidth  = bbox.xmax - bbox.xmin;
        boundsHeight = bbox.ymax - bbox.ymin;
        if (boundsWidth <= 0 || boundsHeight <= 0)
                elog(ERROR, "%s: Geometric bounds are too small", __func__);
 
        if (zoom < 0 || zoom >= 32)
                elog(ERROR, "%s: Invalid tile zoom value, %d", __func__, zoom);
 
        zoomu = (uint32_t)zoom;
        worldTileSize = 0x01u << (zoomu > 31 ? 31 : zoomu);
 
        if (x < 0 || (uint32_t)x >= worldTileSize)
                elog(ERROR, "%s: Invalid tile x value, %d", __func__, x);
        if (y < 0 || (uint32_t)y >= worldTileSize)
                elog(ERROR, "%s: Invalid tile y value, %d", __func__, y);
 
        tileGeoSizeX = boundsWidth / worldTileSize;
        tileGeoSizeY = boundsHeight / worldTileSize;
        x1 = bbox.xmin + tileGeoSizeX * (x);
        x2 = bbox.xmin + tileGeoSizeX * (x+1);
        y1 = bbox.ymax - tileGeoSizeY * (y+1);
        y2 = bbox.ymax - tileGeoSizeY * (y);
 
        PG_RETURN_POINTER(
                geometry_serialize(
                lwpoly_as_lwgeom(
                lwpoly_construct_envelope(
                        srid, x1, y1, x2, y2))));
}

References geometry_serialize(), LW_SUCCESS, lwgeom_calculate_gbox(), lwgeom_free(), lwgeom_from_gserialized(), lwpoly_as_lwgeom(), lwpoly_construct_envelope(), LWGEOM::srid, pixval::x, GBOX::xmax, GBOX::xmin, pixval::y, GBOX::ymax, and GBOX::ymin.

Here is the call graph for this function: