[R-sig-Geo] [raster] conserving dimensions when regridding a 4D Brick or Stack

Tom Roche Tom_Roche at pobox.com
Tue Apr 9 21:27:23 CEST 2013


https://stat.ethz.ch/pipermail/r-sig-geo/2013-April/017947.html
>> Can one

>> + create a 4D RasterBrick or RasterStack
>> + instruct projectRaster

>> so as to get 4D regridded output? Or must one

>> - create 3D RasterLayer's (i.e., one per LAY) from the 4D Brick
>> - regrid each 3D Layer
>> - reassemble the 4D Brick

https://stat.ethz.ch/pipermail/r-sig-geo/2013-April/017973.html
> A RasterBrick is 3D [only].

also RasterStack, no? It Would Be Nice if the docs made this explicit--
or am I missing something?

> all you can do with the raster package is process the data level by
> level

I guessed:

# regridding, including the 4D 3hourlies
https://bitbucket.org/tlroche/gfed-3.1_global_to_aqmeii-na/src/f5e1fd0d607da8bbfe67bb180aaaa21dcf28e1f2/vis_regrid_vis.r?at=master
# reassembling the 4D 3hourlies
https://bitbucket.org/tlroche/gfed-3.1_global_to_aqmeii-na/src/f5e1fd0d607da8bbfe67bb180aaaa21dcf28e1f2/fix_3hourlies.ncl?at=master



More information about the R-sig-Geo mailing list