如何从每日 netcdf 计算每月数据答案

【问题标题】：How to calculate monthly data from a daily netcdf如何从每日 netcdf 计算每月数据
【发布时间】：2015-04-14 10:35:06
【问题描述】：

我能够使用包“ncdf4”在 R 中打开一个 netcdf 文件并提取降雨量（使用 ncvar_get）。但是，我找不到有关如何从每日值计算每月值的信息。在保持 netcdf 格式时这可能吗，还是应该将其转换为其他格式？ - 但是什么格式可以保留与每个降雨值相关的纬度、经度和时间信息？我对 netcdfs 不熟悉，因此我们将不胜感激。

【问题讨论】：

对于这种数据分析，我强烈建议您使用Climate Data Operators (CDO)。它既快速又简单。如果要执行高级分析，CDO 不是很好，但在这种情况下效果很好。您要查找的运算符是 monmean

标签： r netcdf cdo-climate

【解决方案1】：

cdo 的确切命令是

cdo monmean in.nc out.nc

我在此处的 cdo 中发布了 video tutorial on temporal statistics。

编辑：

在 J W 的评论中，您也可以计算“平均值”而不是“平均值”

cdo monavg in.nc out.nc

两者之间的区别在于它们处理缺失数据的方式。 monmean 忽略缺失值，而 monavg 包括缺失值。引用手册：“1、2、miss 和 3 的平均值是 (1+2+3)/3 = 2，而平均值是 (1+2+miss+3)/4 = miss/4 = miss” - see section 2.8 of the manual 了解更多详情。

我已经发布了video tutorial on the difference between mean and avg here。

最后但同样重要的是，您还可以直接从 R using this package which provides R wrappers for CDO 使用 CDO 命令

【讨论】：

还有cdo monavg in.nc out.nc
是的，很好的补充，谢谢——但读者应该注意，如果数据中有缺失值，这些计算是不一样的——monmean 忽略缺失值，而 monavg 包括缺失值。引用手册：“1、2、miss 和 3 的平均值为 (1+2+3)/3 = 2，而平均值为 (1+2+miss+3)/4 = miss/4 = miss” - 有关详细信息，请参阅第 2.8 节。

【解决方案2】：

我不是 R 程序员，但我可以推荐命令行工具 nco 和/或 cdo 进行此类操作

http://nco.sourceforge.net/nco.html -- 见平均值 nces、ncra 和 ncwa。你可能想要记录平均器ncra

https://code.zmaw.de/projects/cdo -- 查看https://code.zmaw.de/projects/cdo/embedded/1.6.4/cdo_refcard.pdf 获取每月、每日、季节性、统计的命令

【讨论】：