NetCDF 到光栅砖“无法为‘ncdf4’找到函数‘砖’的继承方法”

Question

光栅包的问题非常简单，也使用 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，将其推入数组）时，我没有收到相同的错误，但 raster 包说它直接支持 .nc 文件（这么久因为它们与 cf-1 兼容，而 Era-Interim .nc 是）

非常感谢任何帮助，已经尝试过许多 Netcdf 文件（非 Era Interim 也是如此）。

Answer 1

感谢Renaud Lancelot，他给出了清晰的源代码。我已经修改了他的代码以适应你的数据

 # 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)