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

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

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

x

本帖最后由 高空急流 于 2018-9-21 09:41 编辑

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

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

python setup.py install复制代码设置图片保存路径

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

import cinrad

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雷达数据。例子：

from cinrad.io import CinradReader

f = CinradReader(your_radar_file)

f.reflectivity(elevation_angle_level, data_range) #获取反射率数据(为cinrad.datastruct.R类型)

f.velocity(elevation_angle_level, data_range) #获取速度数据(为cinrad.datastruct.V类型)

f.rhi(azimuth, drange) #获取RHI数据(为cinrad.datastruct.Section类型)复制代码cinrad.utils

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

函数名： composite_reflectivity, echo_tops, vert_integrated_liquidcinrad.easycalc

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

函数名： quick_cr, quick_et, quick_vil

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

列表生成示例：

r_list = [f.reflectivity(i, drange) for i in f.angleindex_r] #SA/SB/CA/CB雷达

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

例子：

from cinrad.visualize.ppi import base_reflectivity

base_reflectivity(R) #绘制基本反射率图片

from cinrad.visualize.rhi import rhi

rhi(Section) #绘制RHI复制代码最新的模块下载：

PyCINRAD-master.zip

(9.36 MB, 下载次数: 658)

2018-7-9 20:41 上传

点击文件名下载附件

PyCINRAD脚本

以下是效果图

Z9270_20180708075600_0.6_230_R.png (867.53 KB, 下载次数: 6)

2018-7-9 20:33 上传

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

Z_RADR_I_Z9551_20170714075600_O_DOR_SA_CAP_0.5_230_V.png (911.41 KB, 下载次数: 1)

2018-7-9 20:33 上传

图2  合肥站雷达基速度

Z_RADR_I_Z9002_20170705075800_O_DOR_SA_CAP_RHI_120_85.png (182.48 KB, 下载次数: 2)

2018-7-9 20:33 上传

图3  青浦站雷达距离高度显示

声明：本程序由HCl开发，最终解释权归HCl所有，切勿用于商业用途。(本人仅参与一小部分工作，经HCl同意在此发布)

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

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