确定如果一个高分辨率纬度/经度点在于一个粗网格分辨率纬度/经度箱子吗?

Question

我有两个数据集,每个都有两个不同的决议。第一个数据集是一个粗分辨率数据集,0.5 x 0.625度分辨率和包含MERRA-2边界层厚度值。第二,细分辨率数据集有0.1 x 0.1度分辨率和包含IMERG降雨率的值。我想找的边界层厚度在每一点上降雨率超过10毫米/小时。自从MERRA-2边界层厚度数据集是粗的,我想,我可以创建一个网格的纬度和经度值MERRA-2数据集和检查如果IMERG数据点位于任何MERRA-2网格框。我很难弄清楚从哪里开始,和任何帮助我将非常感谢。

这是每个数据集的纬度和经度的格式:

MERRA-2纬度:[[1.3500000 e + 01 1.3500000 e + 01 1.3500000 e + 01 1.3500000 e + 01] MERRA-2经度:[[-51.875 -51.25 -50.625 -50。IMERG纬度:-49.375 -48.75 -48.125 -47.5 -46.875 -0.05 0.05 0.15 0.25 0.35 0.45 0.55 0.65 0.75 0.85 0.95 1.05]IMERG经度:[-49.95 -49.85 -49.75 -49.65 -49.55 -49.45 -49.35 - -49.25 -49.15 - -49.05)

编辑:我已经创建了一个从Nemesi regridded数据集的第一个方法(cdo remapcon),下面可以找到他们。上图是regridded数据,底部是原始数据。虽然在新regridded数据集的点的数量从59040增加到5475600,形象似乎比原来的低分辨率图像。IMERG数据也似乎更高分辨率虽然数据集都有相同数量的点。这是由于一个错误的数据格式?

Answer 1

注意:我假设NetCDF提供了原始数据格式。

有两个(好,有更多的,但我将重点放在这两个)相对快速处理这个问题的方法。哪一个更喜欢在很大程度上取决于你的偏好和能力使用不同的工具。

我想的最快的方法是使用气候数据操作符(CDO):你可以找到详细的文档在这里和概述在这里。

这种方法应该是获得最高分辨率的网格描述文件:

cdo -griddes imerg.nc

然后,把网格的描述它. txt文件和重新映射粗分辨率文件来匹配这个网格。在这种情况下,您可能想使用一个保守的重新映射方法来维持原来的值没有改变。

cdo remapcon imerggrid。txt merra_coarse。数控merra_imerggrid.nc

这样你会MERRA-2数据分辨率匹配你IMERG网格。

第二种方法使用R。步骤(很快):

1. 加载两个NetCDF文件光栅;
2. 光栅点转换(或者只是提取在一个点的坐标和变换空间对象)的imerg光栅(这样你会获得的质心imerg细胞);
3. 提取数据从一对象覆盖imerg点。所有这些文章都是有据可查的stackoverflow或gis.stackexchange。

然而,这最后一种方法可能会导致一个非常缓慢的过程在大型数据集的情况下。