将R光栅堆栈写入NetCDF

Question

我获得了一个包含1981年每月温度数据的R网格文件，我阅读了该数据，并试图使用以下代码将其写入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网格文件中的数据？尤其是每个月的数据！

编辑：

新的工作守则：

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

主要是将raster替换为brick。另外，这三个点.../是没有意义的。它可能是一两个点(或不必要的)，而package = "raster"参数是没有意义的。