最近30从www.hoelymoley.com 2023 - 07 - 09 - t21:22:12z //www.hoelymoley.com/feeds/question/22612 https://creativecommons.org/licenses/by-sa/4.0/rdf //www.hoelymoley.com/q/22612 3 mpletch1 //www.hoelymoley.com/users/20707 2021 - 07 - 26 - t18:44:24z 2021 - 07 - 28 - t19:02:11z

< p >我有两个数据集,每个都有两个不同的决议。第一个数据集是一个粗分辨率数据集,0.5 x 0.625度分辨率和包含MERRA-2边界层厚度值。第二,细分辨率数据集有0.1 x 0.1度分辨率和包含IMERG降雨率的值。我想找的边界层厚度在每一点上降雨率超过10毫米/小时。自从MERRA-2边界层厚度数据集是粗的,我想,我可以创建一个网格的纬度和经度值MERRA-2数据集和检查如果IMERG数据点位于任何MERRA-2网格框。我很难弄清楚从哪里开始,和任何帮助我将非常感谢。< / p > < p >这是经度和纬度的格式为每个数据集:< / p > < pre > <代码> 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)< /代码> < / pre > < p >编辑:我已经创建了一个从Nemesi regridded数据集的第一个方法(cdo remapcon),下面可以找到他们。上图是regridded数据,底部是原始数据。虽然在新regridded数据集的点的数量从59040增加到5475600,形象似乎比原来的低分辨率图像。 The IMERG data also seems to have higher resolution even though both datasets have the same amount of points. Is this due to an error with the data formatting?

//www.hoelymoley.com/questions/22612/-/22618 # 22618 5 Nemesi //www.hoelymoley.com/users/8087 2021 - 07 - 27 - t08:35:44z 2021 - 07 - 27 - t08:35:44z

< p >注意:我假设NetCDF提供了原始数据格式。< / p > < p >有两个(好,有更多的,但我将重点放在这两个)相对快速处理这个问题的方法。哪一个更喜欢在很大程度上取决于你的偏好和能力使用不同的工具。< / p > < p >我想的最快的方法是使用气候数据操作符(CDO):你可以找到一个详细文档< a href = " https://code.mpimet.mpg.de/projects/cdo/embedded/cdo.pdf " rel = " noreferrer " > < / >和概述< a href = " https://code.mpimet.mpg.de/projects/cdo/wiki/Tutorial " rel = " noreferrer " > < / >。< / p > < p >的方法应该是获得最高分辨率的网格描述文件:< / p > < pre > <代码> cdo -griddes imerg。数控代码< / > < / pre > < p >然后,网格的描述,把它放到<代码>。txt > < /代码文件和重新映射粗分辨率文件来匹配这个网格。在这种情况下,您可能想使用一个保守的重新映射方法来维持原来的值没有改变。< / p > < pre > <代码> cdo remapcon, imerggrid。txt merra_coarse。数控merra_imerggrid。数控代码< / > < / pre > < p >这样你会MERRA-2数据分辨率匹配你IMERG网格。< / p > < p >第二种方法使用<代码> R > < /代码。步骤(很快):< / p > < ol > <李>加载两个NetCDF文件作为光栅;< /李> <李>使光栅点转换(或者只是提取在一个点的坐标和变换空间对象)的imerg光栅(这样你就会获得imerg细胞)的重心,李< / > <李>提取数据从一对象覆盖imerg点。所有这些文章都是记录在< a href = " https://stackoverflow.com/ " > stackoverflow < / >或< a href = " https://gis.stackexchange.com/ " > gis.stackexchange < / >。< /李> < / ol > < p >这最后一种方法,然而,可能会导致一个非常缓慢的过程在大型数据集的情况下。< / p >

map