NetCDF 到 Raster Brick“无法找到‘ncdf4’的函数‘brick’的继承方法”

光栅包的问题非常简单，还使用 ​​ncdf4 加载 ECMWF Era-Interim Netcdf 文件。

只需这样做：

a <- nc_open("SSTs.nc")
B <- brick(a, varname="sst")

返回这个：

    Error in (function (classes, fdef, mtable)  : 
  unable to find an inherited method for function ‘brick’ for signature ‘"ncdf4"’

该文件只是全球 1 个月（2016 年 1 月）的 SST 数据。当我将其转换为数组时（即提取维度/变量，并将时间转换为 UTC，将其推入数组中），我没有得到相同的错误，但光栅包表示它直接支持 .nc 文件（所以只要因为它们与 cf-1 兼容，而 Era-Interim .nc 也是如此）

非常感谢任何帮助，已尝试使用许多 Netcdf 文件（也包括非 Era Interim）。

thank 雷诺·兰斯洛特，谁给出明确的源代码。我修改了他的代码以适应您的数据

 # load package
 library(sp)
 library(raster)
 library(ncdf4)

 # read ncdf file
 nc<-nc_open('D:/SSTs.nc')

 # extract variable name, size and dimension
 v <- nc$var[[1]]
 size <- v$varsize
 dims <- v$ndims
 nt <- size[dims]              # length of time dimension
 lat <- nc$dim$latitude$vals   # latitude position
 lon <- nc$dim$longitude$vals  # longitude position

 # read sst variable
 r<-list()
 for (i in 1:nt) {
   start <- rep(1,dims)     # begin with start=(1,1,...,1)
   start[dims] <- i             # change to start=(1,1,...,i) to read    timestep i
   count <- size                # begin with count=(nx,ny,...,nt), reads entire var
   count[dims] <- 1             # change to count=(nx,ny,...,1) to read 1 tstep

   dt<-ncvar_get(nc, varid = 'sst', start = start, count = count)

   # convert to raster
   r[i]<-raster(dt)
 }

 # create layer stack with time dimension
 r<-stack(r)

 # transpose the raster to have correct orientation
 rt<-t(r)
 extent(rt)<-extent(c(range(lon), range(lat)))

 # plot the result
 spplot(rt)