與shape文件，錯誤colnames

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) 

繪圖示出了兩個nceverShp和是2D和在空間上重疊。感謝有關此錯誤的任何建議。

Answer 1

如果我更換

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

與

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

如預期的面具和作物命令的工作。