微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

在 R 或 CDO 中,随着时间的推移,多个 netcdf 层之间的每日时间步长插值

如何解决在 R 或 CDO 中,随着时间的推移,多个 netcdf 层之间的每日时间步长插值

我对空间变量(假设 var_A_1985.nc、var_A_2000.nc 和 var_A_2010.nc)随着时间的推移进行了多次观察,只有一个时间层。是否可以通过线性插值这些值来有效地重建从 01/01/1985 到 01/01/2010 的每日值?

我目前探索了两种选择:

  • 气候数据运算符 (CDO):我看到了运算符 inttime 和 intyear,但这些似乎不适用于我的情况。

  • R 作为 RasterBrick 元素操作

类似:

from PyInstaller.utils.hooks import exec_statement

mpl_data_dir = exec_statement(
    "import matplotlib; print(matplotlib._get_data_path())")

datas = [
    (mpl_data_dir,"matplotlib/mpl-data"),]

第二种方法确实有效,但需要很长时间。有没有更有效的方法来做到这一点?

解决方法

您可以使用 CDO 执行此操作,如下所示:

cdo -inttime,1985-01-01,12:00:00,1day -mergetime var*.nc outfile

当然要调整开始时间。

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。