https://github.com/CyanideCN/PyCINRAD/blob/master/README_zh.md
https://gitee.com/CyanideCN/PyCINRAD
import cinrad
file = 'Z_RADR_I_Z9535_20210508095956_O_DOR_SAD_CAP_FMT.bin.bz2'
f = cinrad.io.StandardData(file)
print(f.name, "雷达角度:", f.el)
print(f.name, f.code, f.scantime, f.stationlat, f.stationlon, f.radarheight)
rl = list(f.iter_tilt(230, 'REF'))
# f.get_data(elevation, range, datatype)方法,datatype是数据类型,包括'REF', 'VEL', 'SW'
# 以及一些双偏振要素,range是数据的半径,单位是公里,elevation是仰角序号,从零开始计数。
# rl = [f.get_data(i, 230, 'REF') for i in f.angleindex_r]
# 组合反射率
cr = cinrad.calc.quick_cr(rl)
# ppi出图
fig = cinrad.visualize.PPI(cr)
fig('test.png')
vil = cinrad.calc.quick_vil(rl) # 计算VIL
et = cinrad.calc.quick_et(rl) # 计算ET