我有ncdf格式的时间序列Global MODIS L3SMI数据集, 如何基于经度/纬度范围在特定位置(即印度尼西亚区域)上将此数据集子集化。然后我想将此数据集写入与原始数据集相同的ncdf文件中,而无需更改ncdf结构(即尺寸和变量名称)
可以从此链接https://oceandata.sci.gsfc.nasa.gov/cgi/getfile/A2018064.L3m_DAY_SST_sst_4km.nc下载样本数据集
我尝试使用transNcdfSubset
包中的ncdf.tools
,但给我错误消息:
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'))