使用 shapefile 遮罩栅格，列名错误

Question

在R Raster package 中，当我将netCDF 文件加载到raster 并尝试mask 或crop 带有shapefile 的光栅时，出现错误：

"*Error in `colnames<-`(`*tmp*`, value = names(object)) : 
  attempt to set 'colnames' on an object with less than two dimensions*"

这是示例代码

library(ncdf4)
library(raster)

nc <- brick("in.netcdf.nc", varname = "diff", values=TRUE )
crs(nc) <- "+proj=utm +zone=17 +ellps=GRS80 +datum=NAD83 +units=m +no_defs"
everShp <- shapefile("EVERareas")

nc.sub <- mask(nc, everShp)

绘图显示nc 和everShp 都是二维的并且在空间上重叠。感谢您提供有关此错误的任何建议。

Answer 1

如果我替换

nc <- brick("in.netcdf.nc", varname = "diff", values=TRUE )

与

nc <- raster("in.netcdf.nc", varname = "diff", values=TRUE )

遮罩和裁剪命令按预期工作。