用pyart绘制雷达图像并进行反射率修正

前言

这是一篇关于使用 Python 库 Py-ART 读取雷达数据文件、绘制雷达图并进行反射率修正的教程。

这篇文章是我在2020年8月做大创项目的时候写的，由于中间修改过代码，并不确定代码 100% 正确。

原地址：Py-ART 简易中文教程 随着gitee关闭了pages的功能，现在已经无法访问了。

反射率衰减修正

下面是从nc文件中读取数据，计算修正反射率衰减，绘制出图像的例子。

#!/usr/bin/python3

#引入所依赖的第三方库  

import matplotlib.pyplot as plt
import pyart
from sys import exit

if __name__=='__main__':  
	
    #读取数据文件  
	
    radar=pyart.io.read_cfradial('CfRadial_1.nc')  
    
    #检查扫描方式
    
    if radar.scan_type!='ppi':
    	exit('Error: 请使用扫描方式为PPI的数据文件')
    
    #计算修正，返回计算出的衰减(spec_at)和修正后的反射率(cor_z)  
    #radar 是读取数据文件返回的对象，0 是水平反射率的偏移量