Problem with adding a raster and a brick

On Thu, 24 May 2018 at 18:41 Mark R Payne <markpayneatwork at gmail.com>
wrote:

Hi,

I seem to be having a problem adding the following two raster objects
together - one is a rasterLayer, the other is a rasterBrick. The extent,
resolution, and origin are the same, so according to my understand it
should work. The objects look like so:

obs.clim

class : RasterLayer
dimensions : 60, 200, 12000 (nrow, ncol, ncell)
resolution : 0.5, 0.5 (x, y)
extent : -70, 30, 50, 80 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
data source :
/home/mpayne/Documents/Predictability_engine/scratch/Bluefin/HadISST/
obs_climatology.nc
names : sst
z-value : 1988-02-15
zvar : sst

mdl.anom

class : RasterBrick
dimensions : 60, 200, 12000, 1 (nrow, ncol, ncell, nlayers)
resolution : 0.5, 0.5 (x, y)
extent : -70, 30, 50, 80 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
data source :

/home/mpayne/Documents/Predictability_engine/scratch/Bluefin/DCPP-hindcasts/IPSL-CM5A-LR/B.realmean/IPSL-CM5A-LR_S19610101_L000000_realmean_
anom.nc

names : X1961.08.16
Date : 1961-08-16
varname : tos

When I try to add them together, I get the follow error:

mdl.val <- obs.clim + mdl.anom

Error in .readRowsNetCDF(x = x, row = row, nrows = nrows, col = col, ncols
= ncols) :
no slot of name "band" for this object of class ".MultipleRasterData"

I can reproduce, with netcf-file-backed bricks. I would try

obs.clim[[1]] + mdl.anom[[1]]

or

readAll(obs.clim) + readAll(mdl.anom)

those differ mainly in what will happen to multi-layer versions of the
files, the first will give just the first layer. The key is breaking the
file-read-delay, both objects have no data in memory but the add operation
invalidates and causes read to occur - but with a bug - (I will check in
dev raster and report if necessary).

FYI, R-Sig-Geo is generally a better forum for spatial stuff: