pyCINRAD 雷达数据

重磅更新！经过不断的开发和完善，此项目的脚本已经变成模块。衍生产品的计算函数可以支持自己构建的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

   
   
   image.png

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

image.png

图2 合肥站雷达基速度

image.png

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

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

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