如何將 R 中的列表轉換為柵格？

Question

我想將我必須的列表轉換為柵格，以便允許我在特定位置提取數據。

#create example data
tmintest=array(1:100, c(30,189,420))

#create example list
Variable <- list(Varname="Tmin")
Data = tmintest
xyCoords <- list(x = seq(-40.37,64.37,length.out=420), y = seq(25.37,72.37,length.out=189))
Dates <- list(start = seq(as.Date("2036-01-01"), as.Date("2065-01-01"), by="years"), end=seq(as.Date("2036-12-31"), as.Date("2065-12-31"), by="years"))
All <- list(Variable = Variable, Data=Data, xyCoords=xyCoords,Dates=Dates)

我從光柵文檔中看到它不接受列表，所以我嘗試取消列出All

unlistAll<-unlist(All)
raster(unlistAll)

返回：

Error in .local(.Object, ...) 

因此，我將使用上面列表中的信息從頭開始創建它：

myraster<-raster(nrows=108, ncols=21, xmn=-40.375, xmx=64.375,ymn=25.375,ymx=72.375,vals=All$Data)

但我得到了錯誤：

Error in setValues(r, vals) : values must be a vector

我會以錯誤的方式解決這個問題嗎？ 我只想從我的列表或柵格中提取特定 x/y 坐標的數據。

Answer 1

要獲得程度：

xy <- cbind(x=rep(x, length(y)), rep(y, each=length(x)))
r <- rasterFromXYZ(xy)

創建一個 RasterBrick 並設置范圍

b <- brick(tmintest)
extent(b) <- extent(r)

你是如何得到這些值的？ 如果你從 ncdf 文件中讀取它們，你可以只用brick(filename)