确定一个细分辨率的纬度/长度点是否位于一个粗分辨率的纬度/长度网格盒中?

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.3500000e+01 1.3500000e+01 1.3500000e+01 1.3500000e+01] MERRA-2经度:[[-51.875 -51.25 -50.625 -50.]-49.375 -48.75 -48.125 -47.5 -46.875] IMERG纬度:[-0.05 0.05 0.15 0.25 0.25 0.35 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的第一种方法(cdo remapcon)创建了一个受尊重的数据集，它们可以在下面找到。上面的图片是经过处理的数据，下面是原始数据。虽然新数据集中的点数量从59040增加到5475,600，但图像的分辨率似乎比原始图像低。IMERG数据似乎也具有更高的分辨率，即使两个数据集具有相同数量的点。这是因为数据格式错误吗?

Answer 1

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

有两种(好吧，还有更多，但我只关注这两种)相对快速的方法来处理这个问题。选择哪一个在很大程度上取决于您的偏好和使用不同工具的能力。

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

方法应该是从最高分辨率的文件中获取网格描述:

Cdo -griddes imerge .nc

然后，获取网格的描述并将其放入. txt文件并重新映射较粗的分辨率文件以匹配此网格。在这种情况下，您可能希望使用保守的重新映射方法来保持原始值不变。

Cdo remapcon,imerggrid.txt merra_coarse。数控merra_imerggrid.nc

通过这种方式，您将获得与IMERG网格匹配的分辨率的MERRA-2数据。

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

1. 加载两个NetCDF文件作为栅格;
2. 对imerg栅格进行栅格到点的转换(或者只是提取坐标并将它们转换到点空间对象中)(通过这种方式，您将获得imerg单元的质心);
3. 从覆盖imerg点的merra对象中提取数据。所有这些段落都在stackoverflow或gis.stackexchange．

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