C波段双偏振雷达数据处理和可视化

bugsuse

发布于 2020-04-21 16:20:27

2.7K0

发布于 2020-04-21 16:20:27

关于双偏振雷达数据处理和可视化之前在github发过matlab版的程序，以前的推送也专门说过气象数据处理：气象雷达数据II。之所以想要再次更新是因为Python中有了更好的处理雷达数据的库--PyART，相较于之前发布的matlab程序而言，整体的设计都要好太多，所以就有了加入国内雷达数据到此库的想法。

国内S波段雷达数据读取的API已经添加了，而C波段双偏振多普勒雷达数据的读取API一直搁浅，其实整个程序在去年夏天已经完成，但因为存在一些小问题，一放就是差不多一年时间，这两天抽个时间把问题解决了。算是完成了今年的一个小任务了。

先附上链接：https://github.com/bugsuse/pyart

测试代码如下：

import pyart

from pyart.io import c98dfile_archive

import matplotlib.pyplot as plt

radar, nfile = c98dfile_archive('NUIST.20140928.070704.AR2')

display = pyart.graph.RadarDisplay(radar)
fig = plt.figure(figsize=(12, 9))

# plot super resolution reflectivity
ax = fig.add_subplot(111)
display.plot('dBZ', 0, title='WSR-98D Reflectivity',
             vmin=0, vmax=70, colorbar_label='', ax=ax)

所得图形如下：