日志
python读取雷达基数据代码并绘图
|
import matplotlib.pyplot as plt
import pyart
from glob import glob
import os
filedir = 'C:/radar/eg/'
flist = sorted(glob(os.path.join(filedir,'*.nc')))
for filename in flist:
#create the plot using radardisplay
radar = pyart.io.read(filename)
display = pyart.graph.RadarDisplay(radar)
fig = plt.figure(figsize=[5,5])
ax = fig.add_subplot(111)
display.plot('dBZ',2,vmin=-16,vmax=64.0,title='reflectivity')
ax.set_xlim([-200,200])
ax.set_ylim([-200,200])
plt.savefig(os.path.join(filedir,filename.replace('.nc','.png')),dpi=300)
plt.close()