如何在R中使用双线性插值来降低netCDF的分辨率(Regrid)?

时间:2016-02-10 20:04:02

标签: r geospatial netcdf cdo-climate

我有从here下载的netCDF文件。它们的分辨率为0.5 * 0.5。我想以1 * 1的较粗分辨率重新格式化这些文件。我找到了一些链接。 First link讨论了R中的重新网格,但没有使用双线性插值。 second link处理双线性插值,但使用气候数据运算符(我不太熟悉)。然后我遇到了一个R包HiClimR。在此程序包中,命令coarseR会降低数据的分辨率。我将netCDF文件转换为excel文件并使用coarseR。但在得到结果后,我发现这个命令实际上以某种方式跳过经度纬度,并将分辨率降低到1 * 1。在坚果壳中,我的问题是

(1)使用coarseR来降低分辨率是否正确? (2)如何在R?

中使用双线性变换来解决我的具体问题

非常感谢提前。

1 个答案:

答案 0 :(得分:2)

这个解决方案不在R中,只是指出CDO内置了一些非常好的重新嵌入函数。例如,使用双线性插值重新编程为1x1规则网格,您可以这样做:

cdo remapbil,r720x360 in.nc out.nc

但是,如果要转换为较粗糙的网格,则建议使用保守的重新映射技术,否则在重映射期间可能会遗漏点。这对于诸如降水等高度异质的领域尤为重要。在这种情况下,CDO提供一阶和二阶保守重映射技术。使用一阶技术

cdo remapcon,r720x360 in.nc out.nc

(请注意,偶尔您可能会发现CDO因转换过程中的精度损失而摇摆不定,在这种情况下,它会建议您使用选项" -b32"或" -b64&# 34;。除非被迫增加文件大小,否则你不想这样做。)