NetCDF to Raster Brick「無法找到'ncdf4'的函數'brick'的繼承方法」

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"’ 

的文件只是在全球各地SST數據，1個月（Jan2016）。當我將它轉換成一個數組時（即提取維數/變量，並將時間轉換爲UTC，將其轉換爲數組），我沒有得到相同的錯誤，但是柵格包說它直接支持.nc文件因爲它們與cf-1兼容，這是Era-Interim .nc的）

任何幫助非常感謝，試圖用很多Netcdf文件（非時代臨時）。

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)