1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41
| #!/usr/bin/python3
import matplotlib.pyplot as plt import pyart as pa from sys import exit
if __name__=='__main__': radar=pa.io.read_cfradial('plot_PPI.nc')
#检查扫描方式 if radar.scan_type!='ppi': exit('Error: 请使用扫描方式为PPI的数据文件') display=pa.graph.RadarDisplay(radar) fig=plt.figure(figsize=(7,7))
ax=fig.add_subplot(221) display.plot('reflectivity',0,ax=ax,title='reflectivity',colorbar_label='', vmin=-10,vmax=60, axislabels=('','North South distance from radar (km)')) display.set_limits((-250,250),(-250,250),ax=ax)
ax=fig.add_subplot(222) display.plot('differential_reflectivity',0,ax=ax, title='Differential Reflectivity',vmin=-2,vmax=6,colorbar_label='', axislabels=('','')) display.set_limits((-250,250),(-250,250),ax=ax)
ax=fig.add_subplot(223) display.plot('differential_phase',0,ax=ax, title='Differential Phase',colorbar_label='') display.set_limits((-250,250),(-250,250),ax=ax)
ax=fig.add_subplot(224) display.plot('cross_correlation_ratio',0,ax=ax,colorbar_label='', title='Correlation Coefficient', axislabels=('East West distance from radar (km)','')) display.set_limits((-250,250),(-250,250),ax=ax)
plt.show()
|