关于Python用cinrad画SA雷达的时间高度剖面图

爱地理不爱物理 · 发表于 2021-12-3 15:19:28

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

您需要登录才可以下载或查看，没有帐号？立即注册

x

想利用python的cinrad库画sa波段雷达反射率的时间高度剖面图，请问有人画过吗，我目前只会画高度-经（纬）度剖面图，网上找了很多资料都没找到相关教程，求大佬们指教！不吝分享

爱地理不爱物理 · 发表于 2021-12-3 15:22:34

目前有一种想法，是通过python将每个时次的bin文件中反射率数据读取出来统一放到一个nc文件中，再利用ncl画图，不知道是否可行

付亚男 · 发表于 2021-12-3 15:25:29

反射率因子这个变量是为了分析风暴的结构而存在的，这有点类似于高度场，单纯考察一个定点的垂直方向上反射率因子的时间变化，是没有意义的，所以没有人画这种图，同样也极少需要分析某个格点上高度场的时间-高度剖面。
我觉得更可取的方式，应该是采用拉格朗日观点，追踪风暴，考察风暴的垂直结构随时间的演变，业务产品中有这类产品。

edwardli · 发表于 2021-12-3 15:54:26

本帖最后由 edwardli 于 2021-12-3 16:02 编辑

其实不难。
本质上雷达体扫到了Python-xarray眼中，就是个三维数组（dataset）。
其次，不论哪种剖面，就是张二维（Z随着x和y的变化）图形。只要“选”到最后x和y维度，都能画出来“剖面图”，x是经度纬度y是高度就是沿着空间的垂直剖面；x是方位角距离y是高度就是雷达的空间剖面（RHI或者VCS）；x是时间y是高度就是时间垂直剖面。。。
对于雷达反射率因子而言，时间垂直剖面可以判断降水相态演变(0度层高度)、辅助识别列车效应等。

好，说回技术层面，读入各个时次的基数据为dataset结构，增加时间维度以合并（merge）所有时次为一个dataset（类似于最古老的NCEP再分析，一个nc文件一整年好多层好多时次）。

绘制的时候，把方位角和距离固定到目标点，剩下的就是contour/contourf/pcolormesh（t，zlev，dBZ）了。当然中间插值可以很多方式，参考pycinrad或者使用其他科学分析包。

爱地理不爱物理 · 发表于 2021-12-6 09:34:41

付亚男发表于 2021-12-3 15:25
反射率因子这个变量是为了分析风暴的结构而存在的，这有点类似于高度场，单纯考察一个定点的垂直方向上反射 ...

分析某一站点强降水出现的原因也没有意义吗

Pickle. · 发表于 2021-12-10 11:21:27

请问cinrad怎么设成白底的

爱地理不爱物理 · 发表于 2021-12-24 08:36:24

Pickle. 发表于 2021-12-10 11:21
请问cinrad怎么设成白底的

这个程序运行会用到utils.py和PPI.py，你把它们找到然后在里面修改设置

		自动登录	找回密码
密码			立即注册