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")
# 绘制某一仰角的
如何在python中读取多普勒雷达数据?(代码)反射率因子、垂直累计液态水、回波顶高等
于 2024-02-19 16:33:46 首次发布
本文介绍如何在Python中读取以.bin结尾的多普勒雷达数据,包括读取反射率因子和垂直累计液态水等信息。同时,概述了多普勒天气雷达的工作原理,中国新一代天气雷达网的构成,如雷达数据采集子系统(RDA)、雷达产品生成系统(RPG)和主用户终端子系统(PUP)。
摘要由CSDN通过智能技术生成