我发现了栅格砖的意外问题。我创建了一个光栅砖并将其写入文件,但是当我再次加载它并尝试应用一个简单的函数(例如:plot或summary)时,我收到以下错误:
Error in .readCellsRaster(x, uniquecells, layers) : replacement has length zero
代码如下:
rb <- brick(list) # create brick from a list of rasters
masklay <- raster("mask_layer.grd")
rbm <- mask(rb, masklay, inverse=FALSE, maskvalue=NA) # apply a mask
writeRaster(rbm,filename="Brick.grd" ,format="raster", overwrite=FALSE,datatype='INT2U')
但是,在R中加载它不起作用:
rbcheck <- brick("Brick.grd")
plot(rbcheck[[1]], ext=e)
Error in .readCellsRaster(x, uniquecells, layers) : replacement has length zero
我不明白这个错误的含义或如何解决它,我读到它是一个错误,但我不知道该怎么做来处理它。如果它是相关的,砖中的大部分单元都是NA,砖本身很大(4 193 KB)
欢迎任何评论或帮助,
提前致谢!