GPM 半小时降水数据 绘制

1233211234567

第一次写经验帖，请大家多多指教~

       气象数据网上的CMORPH和自动站的融合降水突然下架了，前阵子等来了GPM半小时的降水数据的更新，分辨率0.1°*0.1°，因为是HDF5数据，以前没有使用过，小白摸索了一番，应该是绘制成功了，现在总结一下要注意的地方：

1、数据下载，网上已有很多教程，分享一个:

https://blog.csdn.net/m0_46671865/article/details/121221603?utm_medium=distribute.pc_aggpage_search_result.none-task-blog-2~aggregatepage~first_rank_ecpm_v1~rank_aggregation-18-121221603.pc_agg_rank_aggregation&utm_term=GPM+%E9%99%8D%E6%B0%B4%E6%95%B0%E6%8D%AE%E4%BB%8B%E7%BB%8D&spm=1000.2123.3001.4430
（可能是网速原因，我用python批量下载有点慢，还是一个一个手点下载的快点TAT。)

2、文件名的说明：

比如我下的一个文件命名为：3B-HHR.MS.MRG.3IMERG.20210720-S083000-E085959.0510.V06B

其中20210720开始至V06B前是每个文件不同的信息，20210720则代表年月日，后面的S083000即为从8点30分（UTC）开始，E085959则表示在8点59分59秒结束，即此文件存储的是8点30至8点59分59秒的数据，0510表示从7月20日0点开始的分钟数，510min即8个半小时嘛。

3、绘制变量选择：

在使用NCL绘制脚本前，查阅与数据一起下下来的PDF，点开IMERGE_ATBD_V06.pdf可以发现precipitationCal这一变量是推荐使用的，因此本文使用这个降雨量进行绘图。通过ncl_filedump查看文件描述：

ncl_filedump 3B-HHR.MS.MRG.3IMERG.20210720-S083000-E085959.0510.V06B

4、文件读取：

在README_GPM.pdf中有读取GPM数据的各种示例代码，原文如下：

One excellent source for code snippets in Python, NCL, MATLAB, and IDL can be found at:

https://hdfeos.org/zoo/index_openGESDISC_Examples.php

选择Grid栏中的Python、NCL或MATLAB代码试读取。

其中要注意的是，HDF5文件是0.1°*0.1°的，要将代码中的lon、lat改成如下：

  lon=ispan(0,3600-1,1)*0.1 - 180 + 0.05
  lat=ispan(0,1800-1,1)*0.1 - 90 + 0.05
且原数据是time-lon-lat的摆放顺序，需要通过转置transpose将其转为time-lat-lon才会绘制出正确的区域降水：

data = transpose(data_raw)

5、降水绘制：

下面是轻松愉快的绘图环节，画的是2021年720河南暴雨那一个小时的降水（半小时降水数据相加成一小时）：
因字节限制，图见“一小时降水图”附件
另附两张河南省降水图，感觉比实测要大很多，所以GPM中precipitationCal降水的误差还是值得商讨的：
因字节限制，图见“河南省24h降水图”附件

yq奇

您好 请问下载的逐日的gpm是nc文件应该怎么处理呢？

1233211234567

yq奇 发表于 2021-12-14 10:05
您好 请问下载的逐日的gpm是nc文件应该怎么处理呢？

nc文件除了读取方式不一样，就和hdf5文件一样处理就好了，可以先用ncl_filedump了解一下文件内容

wangmeiyue

{:lxm_24:}{:lxm_24:}

fudong

请问作者下载的是grid的数据吗

fudong

yq奇 发表于 2021-12-14 10:05
您好 请问下载的逐日的gpm是nc文件应该怎么处理呢？

你好 请问你的nc文件读取是怎样的呢

小成都小李亦

yq奇 发表于 2021-12-14 10:05
您好 请问下载的逐日的gpm是nc文件应该怎么处理呢？

你好，想请问一下你的nc文件是怎么下载的？

chinagod

感谢楼主
用ncl绘图时候画不了图，出现了这个问题，看了楼主的贴之后问题解决了：
warning:ContourPlotDraw: out of range coordinates encountered; standard RasterFill rendering method may be unreliable;
看了楼主的贴之后发现，原来原始数据需要转置才能画图。
将数据transpose之后就可以绘图了

其中要注意的是，...原数据是time-lon-lat的摆放顺序，需要通过转置transpose将其转为time-lat-lon才会绘制出正确的区域降水

青青草心

感谢楼主的分享~~请问PrecipitationCal是降水率mm/hr，计算1小时累计降水不应该是前后两个“half-hour"的的降水率求平均然后乘以1hr吗？这样计算的结果肯定是比直接相加小了的

阿天maggie

计算1小时累计降水是前后两个“half-hour"的的降水率求平均然后乘以1hr？