3
\ begingroup美元

我有两个数据集,每个都有两个不同的决议。第一个数据集是一个粗分辨率数据集,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数据也似乎更高分辨率虽然数据集都有相同数量的点。这是由于一个错误的数据格式?

原始数据集新regridded数据集

\ endgroup美元
4
  • 2
    \ begingroup美元 你应该使用GIS方法。你甚至可以使用python(例如geopandas)很快。根据我的经验这是更快的加载地理空间数据(如本机NetCDF ?),因为它是耗时的重新创建网格的纬度/经度文本输出。 \ endgroup美元
    - - - - - -f.thorpe
    2021年7月26日19:50
  • \ begingroup美元 通过GIS方法最近的插值点最近的邻居或你的意思是? \ endgroup美元
    - - - - - -mpletch1
    2021年7月26日20:35
  • \ begingroup美元 如果它只是一个数据集,你可以做什么通常被称为一个“相交”,这是一个类型的空间加入。需要从光栅转换多边形,但这可能太耗费时间如果你处理timeseries数据。如果您使用的是时间序列数据,我建议你使用适当的光栅工具转换从一个网格,让他们直接可比。 \ endgroup美元
    - - - - - -f.thorpe
    2021年7月26日20:47
  • \ begingroup美元 不幸的是,我没有访问ArcMap软件到目前为止,我有困难geopandas Python库安装到我的Python环境。 \ endgroup美元
    - - - - - -mpletch1
    2021年7月29日14:26

1回答1

5
\ begingroup美元

注意:我假设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点。所有这些文章都是有据可查的stackoverflowgis.stackexchange

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

\ endgroup美元
8
  • \ begingroup美元 执行上述方法之后,似乎比以前新regridded图像的分辨率较低。这意味着发生或我某个地方出错了吗? \ endgroup美元
    - - - - - -mpletch1
    2021年7月27日在19:13
  • \ begingroup美元 regridded形象现在应该有相同的分辨率imerg数据 \ endgroup美元
    - - - - - -Nemesi
    2021年7月28日7点15分
  • \ begingroup美元 奇怪是因为点regridded图像的数量从5475600增加到59040,形象似乎不太详细。我要一些图片链接到我最初的问题。 \ endgroup美元
    - - - - - -mpletch1
    2021年7月28日,在福音18:36
  • 1
    \ begingroup美元 图像看起来分散的原因,是由于保守派重新映射。如果你想要一个平滑图像,您应该使用双线性插值重新映射命令(remapbil)。我最初建议remapcon的原因是因为你要求“我想找到边界层厚度每一点降雨率超过10毫米/小时”,我以为你想要原始数据。重新映射与biliniear插值会改变数据平滑。 \ endgroup美元
    - - - - - -Nemesi
    2021年7月29日,8分到达
  • \ begingroup美元 @mpletch1在我看来你也绘制轮廓/插值,这就是为什么粗IMERG数据较低的阴谋出现光滑、更详细。如果你是,那么我建议开关,因为它很难评估regridding。 \ endgroup美元
    - - - - - -Deditos
    2021年7月29日,在大礼堂开幕

你的答案

通过点击“发布你的答案”,你同意我们服务条款并承认您已阅读并理解我们的隐私政策的行为准则

不是你要找的答案?浏览其他问题标记问你自己的问题