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

如何修复我的风倒钩以便它们正确绘制?

如何解决如何修复我的风倒钩以便它们正确绘制?

我一直在尝试使用 Python 在地图上绘制风倒钩以进行研究。在浏览了一些在线指南后,我得出了这个结果,当风倒钩应该在它们的坐标位置绘制时,它们绘制在地图的中心(参见图片链接)。

Wind barb issues

我正在使用来自 WRF 模型运行的数据。这是我的代码如下:

# Get WRF variables
wspd_10m = getvar(ds,'wspd_wdir10',units='ms-1')[0,:]
u10 = getvar(ds,'U10')
v10 = getvar(ds,'V10')
wrf_lats,wrf_lons = latlon_coords(wspd_10m)
wrf_lons = to_np(wrf_lons)
wrf_lats = to_np(wrf_lats)
barb_lons,barb_lats = latlon_coords(u10)
barb_lons = to_np(barb_lons)
barb_lats = to_np(barb_lats)

# Timestamp
timestamp = to_np(wspd_10m.Time).astype('M8[s]').astype('O').isoformat()
time = wspd_10m.Time
time_str = str(time.values)

cart_proj = get_cartopy(wspd_10m)
fig = plt.figure(figsize=(12,6))
ax = plt.axes(projection=cart_proj)
plt.contour(wrf_lons,wrf_lats,wspd_10m,colors='black',transform=ccrs.PlateCarree())
plt.contourf(wrf_lons,cmap=get_cmap('rainbow'),transform=ccrs.PlateCarree())
plt.barbs(barb_lons[::50,::50],barb_lats[::50,to_np(u10[::50,::50]),to_np(v10[::50,length=6)

plot_background(ax)

plt.colorbar(ax=ax,shrink=0.98)

ax.set_extent([-104.35,-94.45,36.37,44.78])
ax.set_title('10m Wind Speed and Direction (m/s) ' + time_str[:19])

plt.show()

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