爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 46952|回复: 138

[经验总结] Python绘制GPM DPR二级降水和潜热

  [复制链接]

新浪微博达人勋

发表于 2022-4-10 12:27:51 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 灭火器 于 2023-9-24 17:16 编辑

更新:最近整理了代码和数据上传到github。
https://github.com/ZhaJiMan/gpm-dpr-visualization
-------------------------------------------------------------------------------

对老帖子python绘制GPM CSH潜热垂直剖面的内容进行了更新,这次完整地画出DPR降水的水平分布和垂直分布,同时比较CSH和SLH潜热的垂直分布。
在选取任意剖面时,由于DPR swath数据是非规则网格,所以我也没想出怎么通过插值得到结果,只好采用取小框框做平均的方法,剖面上每个节点用节点方框内数据点的平均值表示,如下图所示
diagram.png
另一种算法是使用最邻近插值,但我试过之后觉得出图效果没有均值的方法好看。读者可以自己实现看看。
代码使用h5py读取HDF5文件,使用Cartopy包画地图。填色图通过contourf绘制,但也因此无法显示出缺测部分。
游客,如果您要查看本帖隐藏内容请回复

下面以2021年9月29日的台风“蒲公英”为例,具体用到的文件为:
2A.GPM.GMI.GPROF2017v1.20210929-S142530-E155804.043106.V05B.HDF5
2A.GPM.DPR.V8-20180723.20210929-S142530-E155804.043106.V06A.HDF5
2B.GPM.DPRGMI.2HCSHv4-1.20210929-S142530-E155804.043106.V06A.HDF5
2A.GPM.DPR.GPM-SLH.20210929-S142530-E155804.043106.V06B.HDF5
加入GMI是因为DPR的扫描宽度太窄,用更宽的GMI完整表示降水的水平分布。使用的2ADPR产品仍然是V06版,官方供应仅到2021年12月1日,之后的产品均为V07版。因为我对V07的文件结构不熟,并且文件大小暴增,所以这里还是用的老版本。选取了三条剖面,剖面中叠加了零度层高度和地形(可惜洋面看不出地形效果)。
horizontal.png
merged.png





密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2022-4-10 14:34:23 | 显示全部楼层
h5py读取HDF5文件
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-4-10 16:20:10 | 显示全部楼层
以2021年9月29日的台风“蒲公英”为例
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-4-10 17:08:53 | 显示全部楼层
厉害了!!
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2022-4-11 16:11:30 来自手机 | 显示全部楼层
df5文件是从哪里下载的
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2022-4-11 17:15:27 | 显示全部楼层
ZP盛夏光年 发表于 2022-4-11 16:11
df5文件是从哪里下载的

NASA的Earthdata Search网站
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-4-11 17:16:36 | 显示全部楼层
以2021年9月29日的台风“蒲公英”为例
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-4-12 20:59:30 | 显示全部楼层
灭火器 发表于 2022-4-11 17:15
NASA的Earthdata Search网站

好的谢谢楼主,学习了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-4-14 09:18:41 | 显示全部楼层
{:eb348:}{:eb348:}{:eb348:}{:eb348:}
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-4-14 10:05:38 | 显示全部楼层
谢谢楼主 学习了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表