將柵格堆棧寫入NetCDF

Question

我有一個R網格文件，其中包含1981年的每月溫度數據，我已閱讀並嘗試使用以下代碼將其寫入NetCDF：

library(raster)
library(ncdf4)
library(RNetCDF)

test <- raster('.../TavgM_1981.gri', package = "raster")
rstack = stack(test)

writeRaster(rstack, "rstack.nc", overwrite=TRUE, format="CDF",     varname="Temperature", varunit="degC", 
        longname="Temperature -- raster stack to netCDF", xname="X",   yname="Y",zname="nbands",
        zunit="numeric")

這將寫入NetCDF文件 ，但似乎只有一個月（我不確定是哪一個），而不是一年中使用全景查看時的12個月。

是否可以寫入NetCDF文件並保留R-grid文件中盡可能多的數據？ 特別是每個月的數據！

編輯：

新的工作代碼：

test <- brick('/TavgM_1981.gri')
writeRaster(test, "rstack.nc", overwrite=TRUE, format="CDF",     varname="Temperature", varunit="degC", 
        longname="Temperature -- raster stack to netCDF", xname="Longitude",   yname="Latitude", zname="Time (Month)")

Answer 1

正如dww指出的那樣，要獲取所有圖層，

test <- raster('.../TavgM_1981.gri', package = "raster")

應該

test <- brick('TavgM_1981.grd')

最主要的是用brick代替raster 。 而且，三個點.../毫無意義。 它可以是一個或兩個點（或不必要），並且package = "raster"參數是沒有意義的。