[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