如何在python中读取多普勒雷达数据？（代码）反射率因子、垂直累计液态水、回波顶高等

import cinrad

f = cinrad.io.StandardData(
    r"G:\Radar20230601\Z9417Z_RADR_I_Z9417_20230601132355_O_DOR_SAD_CAP_FMT.bin.bz2")
rl = list(f.iter_tilt(230, 'REF'))  # 数据类型（'REF'，即反射率）的数据读取出来。

print(f.available_product(0))
print('*' * 50)

tilt_number = 0
radius = 230
data_dtype = 'REF'
r = f.get_data(tilt_number, radius, data_dtype)  # 获取反射率数据
print(r)
rl = list(f.iter_tilt(radius, 'REF'))
print("-" * 50)

# 绘制某一仰角的反射率图像
# start1 = time.clock()
x = f.angleindex_r  # 获取反射率对应仰角层[ 0  2  4  5  6  7  8  9 10]
r = f.get_data(x[0], 230, 'REF')
fig = cinrad.visualize.PPI(r, dpi=999, add_city_names=True)

# 绘制VIL图像           # 绘制垂直累计液态水
vil = cinrad.calc.quick_vil(rl)
fig = cinrad.visualize.PPI(vil, dpi=999, add_city_names=True)
fig("VIL.png")

# 绘制ET图像             # 绘制回波顶高
# start3 = time.clock()
et = cinrad.calc.quick_et(rl)
fig = cinrad.visualize.PPI(et, dpi=999, add_city_names=True)
fig("ET.png")

# 绘制CR图像             # 绘制组合反射率
cr = cinrad.calc.quick_cr(rl)
fig = cinrad.visualize.PPI(cr, dpi=999, add_city_names=True)
fig("CR.png")

# 绘制某一仰角的