[raster] Implement ST_MakeEmptyCoverage

[dzwarg]

Make a helper method that creates an empty coverage.

An empty coverage is a set of rasters that are regularly tiled that cover the entire pixel area of the coverage, tiling based on the tilewidth/tileheight.

[dzwarg]

Implementation of non-rotating ST_MakeEmptyCoverage()

[dzwarg]

To handle rotation:

Create a single fake raster with as many cells as tiles in the coverage. Iterate over the x,y dimensions of that raster, and use the ULX and ULY of each fake raster pixel to generate a new empty raster with proper ULX and ULY values for the rotated component tiles in the coverage.

[dzwarg]

Implementation of rotatable ST_MakeEmptyCoverage()

[strk]

This is needed for ST_CreateOverviews (#2247)

[strk]

I've taken a look at the function signature in the latest patch and got a bit scared :) I imagine a coverage creation as two steps:

1. Creation of the table that would represent the coverage configuration
2. Population of the coverage

The first step could be similar to AddRasterColumn, that is it only creates a structure, w/out inserting the data.

The way I'd call such a coverage-creation function would be by passing it:

Mandatory:

1. The spatial extent to be covered by the coverage
2. The pixel size in geographical units (the scale)

Optional:

3. SRID. Unconstrained if not given. [ok, maybe this could be mandatory]
4. The max size of coverage tiles in pixels. Unconstrained if not given.

I see this function instead creates a whole set of empty tiles, which is possibly a big waste of time as you probably want to later modify those tiles to have actual data in them...

Also for a big coverage it could be very easy to end up out of memory.

[strk]

Comments about this @dustymugs ?

[strk]

A pointer to the manual page explaining that the whole resultset from a SET-RETURNING plpgsql function is hold in memory (possibly flushed to disk temporarely) until completed: ​http://www.postgresql.org/docs/current/static/plpgsql-control-structures.html#PLPGSQL-STATEMENTS-RETURNING

[pracine]

Replying to strk:

I've taken a look at the function signature in the latest patch and got a bit scared :) I imagine a coverage creation as two steps:

1. Creation of the table that would represent the coverage configuration
2. Population of the coverage

The first step could be similar to AddRasterColumn, that is it only creates a structure, w/out inserting the data.

The way I'd call such a coverage-creation function would be by passing it:

Mandatory:

1. The spatial extent to be covered by the coverage
2. The pixel size in geographical units (the scale)

Optional:

3. SRID. Unconstrained if not given. [ok, maybe this could be mandatory]
4. The max size of coverage tiles in pixels. Unconstrained if not given.

I see this function instead creates a whole set of empty tiles, which is possibly a big waste of time as you probably want to later modify those tiles to have actual data in them...

Also for a big coverage it could be very easy to end up out of memory.

ST_MakeEmptyCoverage() should create EMPTY rasters tiles, thus having no bands and hence almost no data. It is useful for any query that want a set of reference tiles to recreate a raster coverage from an existing one or from a vector coverage (rasterizing). You use, for example, to create on the fly an empty raster table in the FROM geomtable, (SELECT ST_MakeEmptyCoverage() rast) that is used only as reference by function in the SELECT. Many functions have a reference raster as parameter from which they borrow alignment information.

[strk]

Good point, I hadn't considered the rasters were empty ! Then I guess the only thing missing here is documentation in the manual.

[ainomieli]

Here is boilerplate for creating tile set with bounds (:x0, :y0, :x1, :y1) and grid dimensions (:width, :height):

SELECT cell FROM ST_MakeEmptyCoverage(
  1, 1, -- tile width, tile height
  :width, :height, -- grid columns, grid rows
  :x0, :y0, -- upperleftx, upperlefty
  (:x1 - :x0)/(:width)::float, -- scalex
  (:y1 - :y0)/(:height)::float, -- scaley
  0., 0., -- skewx, skewy
  4326 -- srid
  ) cell;

• (:x0, :y0) is the upper left corner
• (:x1, :y1) is the lower right corner
• :width is the number of columns in the grid
• :height is the number of rows in the grid

[ainomieli]

preliminary docs

[robe]

I hate pushing this but I don't have time to test it out and giving the patch is so old, not sure how it plays with our current code.

[robe]

ainomieli,

If you are still around, can you give me your name as you want to appear in credits?

Thanks, Regina

[robe]

In 15562:

ST_MakeEmptyCoverage (David Zwarg, ainomieli)
Closes #2249 for PostGIS 2.4.0