重磅更新！读取CINRAD雷达基数据的Python模块

高空急流

重磅更新！经过不断的开发和完善，此项目的脚本已经变成模块。衍生产品的计算函数可以支持自己构建的numpy.ndarray，而不必须是程序读取出来的CINRAD雷达数据。该模块支持SA/SB/CB三种雷达数据格式，具体安装及函数用法如下。

安装及初始化设置说明：该脚本目前仅能在Python3环境下运行

1. python setup.py install

复制代码

设置图片保存路径

程序默认将图片保存在D:\目录下，如要设置到其他路径，请使用cinrad.set_savepath函数。例子：

1. import cinrad
   
2. cinrad.set_savepath('D:\\1\\')

复制代码

模块介绍
cinrad.datastruct

构建本模块所使用的数据类型

基本类型: cinrad.datastruct.Raw

反射率数据类型: cinrad.datastruct.R (base: cinrad.datastruct.Raw)

速度数据类型: cinrad.datastruct.V (base: cinrad.datastruct.Raw)

剖面数据类型: cinrad.datastruct.Section

该基本类型包含该要素数据，经纬度数据和雷达其他信息（雷达站名，扫描时间等）

cinrad.io

读取CINRAD雷达数据。例子：

1. from cinrad.io import CinradReader
   
2. f = CinradReader(your_radar_file)
   
3. f.reflectivity(elevation_angle_level, data_range) #获取反射率数据（为cinrad.datastruct.R类型）
   
4. f.velocity(elevation_angle_level, data_range) #获取速度数据（为cinrad.datastruct.V类型）
   
5. f.rhi(azimuth, drange) #获取RHI数据（为cinrad.datastruct.Section类型）

复制代码

cinrad.utils

提供雷达衍生产品的计算（接受numpy.ndarray）。将这些功能独立出来的目的是使得计算程序更加通用， 而不仅仅是能计算此程序读取出来的CINRAD雷达数据。

函数名： composite_reflectivity, echo_tops, vert_integrated_liquid

cinrad.easycalc

提供雷达衍生产品的计算（接受list(cinrad.datastruct.Raw)） 使用cinrad.io读取的数据可直接带入该模块下的函数来计算。

函数名： quick_cr, quick_et, quick_vil

传入一个包含每个仰角数据的list即可计算。

列表生成示例：

1. r_list = [f.reflectivity(i, drange) for i in f.angleindex_r] #SA/SB/CA/CB雷达
   
2. r_list = [f.reflectivity(i, drange) for i in range(len(f.elevdeg))] #SC/CC雷达

复制代码

cinrad.visualize

雷达数据可视化，包括ppi和rhi，仅接受cinrad.datastruct.Raw类型。

在ppi下的函数：base_reflectivity, base_velocity, echo_tops, vert_integrated_liquid, composite_reflectivity

在rhi下的函数：rhi

例子：

1. from cinrad.visualize.ppi import base_reflectivity
   
2. base_reflectivity(R) #绘制基本反射率图片
   
3. from cinrad.visualize.rhi import rhi
   
4. rhi(Section) #绘制RHI

复制代码

最新的模块下载： PyCINRAD-master.zip (9.36 MB, 下载次数: 931)

2018-7-9 20:41 上传

点击文件名下载附件
PyCINRAD脚本
以下是效果图

图1  武汉站雷达基本反射率

图2  合肥站雷达基速度

图3  青浦站雷达距离高度显示
声明：本程序由HCl开发，最终解释权归HCl所有，切勿用于商业用途。（本人仅参与一小部分工作，经HCl同意在此发布）

附上github地址，希望大家多多支持哦：https://github.com/CyanideCN/PyCINRAD

有任何问题和想法也欢迎提出，讨论地址：https://github.com/CyanideCN/PyCINRAD/issues

eeeee

更新地区高亮功能
在ppi里面的所有函数里加入了highlight参数。
例子:cinrad.visualize.ppi.base_reflectivity(R, highlight=‘成都市’)
效果图

eeeee

更新一下程序。
新增功能：
1.CC1.0雷达的多普勒速度和基本反射率绘制
2.CA雷达多普勒速度和基本反射率绘制
3.双偏振雷达的基本反射率，差分反射率，协相关系数，差分传播相移率的绘制
4.优化雷达类型的识别
5.一些其他Bug的修复
6.组合反射率绘制


PyCINRAD.zip (9.44 MB, 下载次数: 460)

2018-9-2 19:56 上传

点击文件名下载附件

示例图片

kmdqlb

@高空急流  老师，请问运行出现下列错误，请问是怎么回事？

Traceback (most recent call last):
  File "C:\Users\Administrator\Desktop\PyCINRAD-master\PyCINRAD-master\CINRAD_radar.py", line 26, in <module>
    radarinfo = np.load('RadarStation.npy')
  File "C:\Python27\lib\site-packages\numpy\lib\npyio.py", line 394, in load
    return format.read_array(fid)
  File "C:\Python27\lib\site-packages\numpy\lib\format.py", line 446, in read_array
    array = pickle.load(fp)
TypeError: must be char, not unicode

非对称

佩服楼主的开源精神

kmdqlb

谢谢楼主分享，

kmdqlb

@非对称 老师，上次看了你们江西的GIS系统，好久没见老师发言了

非对称

kmdqlb 发表于 2018-7-10 09:32
@非对称 老师，上次看了你们江西的GIS系统，好久没见老师发言了

我不做预报员好多年了，哈哈

kmdqlb

非对称 发表于 2018-7-10 10:09
我不做预报员好多年了，哈哈

根据你的指导，我结合wcoutour应用在silverlight上生成了一个WebGIS的系统，现已应用于业务中

kmdqlb

非对称 发表于 2018-7-10 10:09
我不做预报员好多年了，哈哈

根据你的指导，我结合wcoutour应用在silverlight上生成了一个WebGIS的系统，现已应用于业务中

kmdqlb

非对称 发表于 2018-7-10 10:09
我不做预报员好多年了，哈哈

根据你的指导，我结合wcoutour应用在silverlight上生成了一个WebGIS的系统，现已应用于业务中

kmdqlb

非对称 发表于 2018-7-10 10:09
我不做预报员好多年了，哈哈

根据你的指导，我结合wcoutour应用在silverlight上生成了一个WebGIS的系统，现已应用于业务中