Thanks,
I basically have to reclass() all 72 tiles, then aggregate() them
from 30 arc seconds to 5 minutes and eventually put them in a brick
with other rasters
so I can average across the layers.
I thought unionExtent() was going to be the ticket, but that just
effects the extent.
I have 72 tiles in a tif format and each tile is 30X30
tile
class : RasterLayer
dimensions : 3600, 3600, 1296 (nrow, ncol, ncell)
resolution : 0.00833, 0.00833 (x, y)
extent : 60, 90, 60, 90 (xmin, xmax, ymin, ymax)
coord. ref. : NA
values : C:/Users/steve/
unionExtent() looks like the best approach.
On Sun, Oct 16, 2011 at 8:54 PM, steven mosher mosherste...@gmail.com wrote:
I have 72 tiles in a tif format and each tile is 30X30
tile
class : RasterLayer
dimensions : 3600, 3600, 1296 (nrow, ncol, ncell)
resolution :
Steven,
In principle you could give merge all tiles (e.g. make a list of RasterLayer
objects and use do.call ) but merge is quite slow (something that can easily
be be improved; and I expect it will be over the coming months).
What I would do is use run gdalbuildvrt (that comes with the
