微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

使用ARM_Pyart访问2级雷达数据中的其他扫描

如何解决使用ARM_Pyart访问2级雷达数据中的其他扫描

全部

我正在尝试使用arm_pyart从2级雷达数据文件访问中翼龙扫描。我的脚本当前仅绘制以文件命名的基本扫描,但是我知道文件包含其他扫描。我已经阅读了所有pyart文档,但找不到任何说明如何访问这些其他扫描的内容。任何在正确方向上的帮助/指向都会有所帮助。下面是我的python代码的示例。

# Import Libraries
import matplotlib.pyplot as plt 
import cartopy.crs as ccrs 
import pyart 
import cartopy 

# Level 2 file name.
file = 'KGRR20160820_173712_V06'

# Takes the date and time from the file name. 
# This is used for saving a figure (Not shown).
name=str(file[4:12])
name1=str(file[13:19])

# Reads the data.
radar=pyart.io.read_nexrad_archive(filename)
display = pyart.graph.RadarMapdisplay(radar)

# Creates map of SW Michigan.
fig=plt.figure(figsize=[16,16]) 
ax = plt.subplot(221,projection=ccrs.PlateCarree())
ax.set_extent([-86.7,-84.8,41.5,43.5])
ax.add_feature(cartopy.feature.STATES,linewidth=0.5,zorder=10)
ax.coastlines(linewidth=0.5,resolution='10m',zorder=0)
ax.add_feature(cartopy.feature.BORDERS,zorder=0)


# Plots the radar reflectivity factor for the 0.5 tilt.   
display.plot_ppi_map(radar.fields['reflectivity']['data'][sweep_1],vmin=10,vmax=70,title='Base \n Reflectivity',colorbar_label='dBZ',cmap=pyart.graph.cm.NWSRef,embelish=False)

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。