编辑-地球科学江南体育网页版堆栈交换江南电子竞技平台

你不是登录。你的编辑将被放置在一个队列,直到同行评议。

我们欢迎编辑,使读者更容易理解和更有价值。因为社区成员审查编辑,请尽量使文章更好比你怎么发现,例如,通过修正语法或添加额外的资源和超链接。

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

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

回答

注意:我假设NetCDF提供了原始数据格式。有两个(好,有更多的,但我将重点放在这两个)相对快速处理这个问题的方法。哪一个更喜欢在很大程度上取决于你的偏好和能力使用不同的工具。我想的最快的方法是使用气候数据操作符(CDO):你可以找到详细的文档(在这里)[1]和[这][2]概述。这种方法应该是获得最高分辨率的网格描述文件:cdo -griddes imerg。数控,把网格的描述”。txt文件和重新映射粗决议文件来匹配这个网格。在这种情况下,您可能想使用一个保守的重新映射方法来维持原来的值没有改变。cdo remapcon imerggrid。txt merra_coarse。数控merra_imerggrid。数控这样你会MERRA-2数据分辨率匹配你IMERG网格。 The second approach uses `R`. Steps (very quickly): 1. load the two NetCDF files as raster; 2. make a raster to point transformation (or just extract the coordinates and transform them in a Point Spatial object) of the imerg raster (in this way you would obtain the centroids of the imerg cells); 3. extract the data from the merra object overlaying the imerg points. All these passages are well documented in [stackoverflow][3] or [gis.stackexchange][4]. This last methodology, however, may result in a very slow process in case of large datasets. [1]: https://code.mpimet.mpg.de/projects/cdo/embedded/cdo.pdf [2]: https://code.mpimet.mpg.de/projects/cdo/wiki/Tutorial [3]: https://stackoverflow.com/ [4]: https://gis.stackexchange.com/

编辑总结

取消

\ 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日,在大礼堂开幕

| 显示3更多的评论

正确的小错误或错误
澄清的意思,而不更改它
添加相关资源或链接
总是尊重作者的意图
不使用编辑与作者回复

创建代码用引号”或式子~篱笆
' ' '
像这样
' ' '
添加语言标识符来突出显示的代码
”“python
def函数(foo):
打印(foo)
' ' '
把段落之间的回报
在结束linebreak添加2个空间
_italic_或* *的* *
报价将>开始的
做链接(尽量使用https)
< https://example.com >
(例子)(https://example.com)
< a href = " https://example.com " > < / >示例
MathJax方程罪\ \ ^ 2θ美元

格式帮助»
回答帮助»

MathJax帮助»