使用shapefile屏蔽栅格，使用colnames中的错误

Question

在R Raster package中，当我将netCDF文件加载到raster并尝试mask或crop带有{{1}的栅格时我收到错误：

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都是2D和空间重叠。感谢您对此错误提出的任何建议。

Answer 1

如果我替换

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

与

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

mask和crop命令按预期工作。