r - 基于经度/纬度的 MODIS L3SMI ncdf 子集

Question

我有 ncdf 格式的时间序列 Global MODIS L3SMI 数据集，如何根据经度/纬度范围在特定位置（即印度尼西亚地区）对该数据集进行子集化。然后我想将此数据集写入与原始数据集相同的 ncdf 文件，而不更改 ncdf 结构上的任何内容（即维度和变量名称）

可以从此链接下载示例数据集https://oceandata.sci.gsfc.nasa.gov/cgi/getfile/A2018064.L3m_DAY_SST_sst_4km.nc

我尝试使用transNcdfSubsetfrom ncdf.toolspackage 但给我错误消息：

transNcdfSubset("A2018064.L3m_DAY_SST_sst_4km.nc", : 文件中的维度名称不是纬度/经度/时间。重命名！

这是我的代码：

transNcdfSubset('A2018064.L3m_DAY_SST_sst_4km.nc', 
            dim.values = list(latitude = c(-25:20), longitude = c(90:145), rgb = c(), eightbitcolor=c()), 
            values.type = c("values", "values", "values", "values")[2], 
            file.output = sub("[.]nc", "_indo.nc", 'A2018064.L3m_DAY_SST_sst_4km.nc'), 
            var.name = readNcdfVarName('A2018064.L3m_DAY_SST_sst_4km.nc'))

标签： rncdf4r-modis