在 R 中使用 cdfextract 时出错(clim.pact 包)

error using cdfextract in R (clim.pact package)

我对使用 R 还很陌生,目前我正在尝试使用它来提取 netCDF 文件特定区域的数据并将其转换为 csv。最后,我想为我感兴趣的区域中的每个变量创建一个时间序列。

我已经为 R 下载了软件包 ncdfclim.pact,并且由于我的 netCDF 文件非常大(25 年的每日数据)我想我应该能够提取一个使用 cdfextract 获取我感兴趣区域的数据,然后我会将那个较小的文件转换为 csv。但是,当我尝试 运行 时,我总是遇到同样的错误。这是我一直在使用的脚本和出现的错误:

cdfextract("file.nc","tas", x.rng=c(-93, -85), y.rng=c(16, 23),t.rng=c(1, 9125))

[1] "time.origin= 01 Jan 0850" Error in system("rm -f cdfcont.txt", intern = TRUE) : 'rm' not found

我的文件是 file.nc,感兴趣的变量是 tas。我认为这可能与时间格式有关,但当我检查文件的时间单位时,它说它在 "days since 01 Jan 0850" 中,所以我认为单位是正确的。

无论如何,欢迎就此问题提出任何意见!谢谢!

如果您不必使用 R,可以使用 NCO 包中的 nckshttp://nco.sourceforge.net/nco.html#ncks-netCDF-Kitchen-Sink

ncks 是一个很棒的命令行工具,用于从 netCDF 文件中提取数据。通常,您可以在某些 lat/lon 波段上提取变量,然后使用类似以下内容将该数据输出到 .csv 文件中:

ncks -v tas -d lat,16.,23. -d lon,-93.,-85. -s '%f\n' file.nc > output.csv

我无法判断您是否首先需要 tas 的面积平均值,然后是 .csv 文件中该面积平均值的时间序列。如果是这样,您需要先使用 ncwa 进行面积平均:http://nco.sourceforge.net/nco.html#Averaging