重磅更新！經(jīng)過不斷的開發(fā)和完善，此項目的腳本已經(jīng)變成模塊。衍生產(chǎn)品的計算函數(shù)可以支持自己構(gòu)建的numpy.ndarray，而不必須是程序讀取出來的CINRAD雷達(dá)數(shù)據(jù)。該模塊支持SA/SB/CB三種雷達(dá)數(shù)據(jù)格式，具體安裝及函數(shù)用法如下。
安裝及初始化設(shè)置****說明：該腳本目前僅能在Python3環(huán)境下運行
- python setup.py install
復(fù)制代碼
設(shè)置圖片保存路徑
程序默認(rèn)將圖片保存在D:\目錄下，如要設(shè)置到其他路徑，請使用cinrad.set_savepath函數(shù)。例子：
import cinrad
cinrad.set_savepath('D:\1\')
復(fù)制代碼
模塊介紹
cinrad.datastruct
構(gòu)建本模塊所使用的數(shù)據(jù)類型
基本類型: cinrad.datastruct.Raw
反射率數(shù)據(jù)類型: cinrad.datastruct.R (base: cinrad.datastruct.Raw)
速度數(shù)據(jù)類型: cinrad.datastruct.V (base: cinrad.datastruct.Raw)
剖面數(shù)據(jù)類型: cinrad.datastruct.Section
該基本類型包含該要素數(shù)據(jù)，經(jīng)緯度數(shù)據(jù)和雷達(dá)其他信息（雷達(dá)站名，掃描時間等）
cinrad.io
讀取CINRAD雷達(dá)數(shù)據(jù)。例子：
from cinrad.io import CinradReader
f = CinradReader(your_radar_file)
f.reflectivity(elevation_angle_level, data_range) #獲取反射率數(shù)據(jù)（為cinrad.datastruct.R類型）
f.velocity(elevation_angle_level, data_range) #獲取速度數(shù)據(jù)（為cinrad.datastruct.V類型）
f.rhi(azimuth, drange) #獲取RHI數(shù)據(jù)（為cinrad.datastruct.Section類型）
復(fù)制代碼
cinrad.utils
提供雷達(dá)衍生產(chǎn)品的計算（接受numpy.ndarray）。將這些功能獨立出來的目的是使得計算程序更加通用， 而不僅僅是能計算此程序讀取出來的CINRAD雷達(dá)數(shù)據(jù)。
函數(shù)名： composite_reflectivity, echo_tops, vert_integrated_liquid
cinrad.easycalc
提供雷達(dá)衍生產(chǎn)品的計算（接受list(cinrad.datastruct.Raw)） 使用cinrad.io讀取的數(shù)據(jù)可直接帶入該模塊下的函數(shù)來計算。
函數(shù)名： quick_cr, quick_et, quick_vil
傳入一個包含每個仰角數(shù)據(jù)的list即可計算。
列表生成示例：
r_list = [f.reflectivity(i, drange) for i in f.angleindex_r] #SA/SB/CA/CB雷達(dá)
r_list = [f.reflectivity(i, drange) for i in range(len(f.elevdeg))] #SC/CC雷達(dá)
復(fù)制代碼
cinrad.visualize
雷達(dá)數(shù)據(jù)可視化，包括ppi和rhi，僅接受cinrad.datastruct.Raw類型。
在ppi下的函數(shù)：base_reflectivity, base_velocity, echo_tops, vert_integrated_liquid, composite_reflectivity
在rhi下的函數(shù)：rhi
例子：
from cinrad.visualize.ppi import base_reflectivity
base_reflectivity(R) #繪制基本反射率圖片
from cinrad.visualize.rhi import rhi
-
rhi(Section) #繪制RHI
image.png
圖1 武漢站雷達(dá)基本反射率
圖2 合肥站雷達(dá)基速度
圖3 青浦站雷達(dá)距離高度顯示
聲明：本程序由HCl開發(fā)，最終解釋權(quán)歸HCl所有，切勿用于商業(yè)用途。（本人僅參與一小部分工作，經(jīng)HCl同意在此發(fā)布）
附上github地址，希望大家多多支持哦：https://github.com/CyanideCN/PyCINRAD
有任何問題和想法也歡迎提出，討論地址：https://github.com/CyanideCN/PyCINRAD/issues
