cinrad 读取雷达 bz2 文件数据

1. 安装
   

   pip install cinrad -i https:mirrors.aliyun.com/pypi/simple

2. 读取文件
   

   import cinrad
   
   filename = r"Z_RADA_C_BABJ_20220720060001_P_DOR_Z9518_BASE_20220821_090303.bin.bz2"
   f = cinrad.io.StandardData(filename)
   print(f)
   print(f.name, "雷达角度：", f.el)
   
   # 查看包含的产品类型
   print(f.available_product(0)) 

   # 获取雷达扫描仰角
   print(f.get_elevation_angles())
   
   # 层数list
   print(f.angleindex_r)

3. 查看数据类型
   

   ele = 0 #选择第1个仰角
   radius = 400 #绘制图像的范围大小，单位km
   r = f.get_data(ele,radius,"REF") #选择反射率数据
   print(r)

4. 雷达回波等
   

   组合反射率  rl = list(f.iter_tilt(230, 'REF'))
   cinrad.easycalc.quick_cr(rl,resolution=[230, 366])
   cr = cinrad.calc.quick_cr([f.get_data(i, 400, 'REF') for i in f.angleindex_r])
   print(cr)
   ppi出图
   fig = cinrad.visualize.PPI(r, dpi=999, add_city_names=True)
   print(fig)
   回波顶高
   et = cinrad.calc.quick_et([f.get_data(i, 400, 'REF') for i in f.angleindex_r])
   print(et)
   
   cr.CR.plot()
   plt.show()
   
   fig = plt.figure()
   ax = Axes3D(fig)  # 创建三维绘图空间
   X = r.longitude.values.flatten()  # 读取ppi中经度纬度高度和反射率数值，并转化成一维
   Y = r.latitude.values.flatten()
   Z = r.height.values.flatten()
   value = r.REF.values.flatten()
   ax.scatter(X, Y, Z, c=value)  # 绘制散点图
   plt.show()

　　　参考： https://cloud.tencent.com/developer/article/1736532