基于nc文件画风场

zhangmeiru

一大碗年糕

今天也碰到这个报错，把xarray的U,V的值要先提取出来才行
lon = mveof['lon']
lat = mveof['lat']
u_mveof = mveof.sel(mode=1)['U'].drop_vars('mode').data
v_mveof = mveof.sel(mode=1)['V'].drop_vars('mode').data
加个.data

咸墨

读取文件创建画布后使用这个函数就可以了
ax.quiver(lon, lat, u, v)

dafenqi

咸墨 发表于 2021-2-6 21:20
读取文件创建画布后使用这个函数就可以了
ax.quiver(lon, lat, u, v)

请问你使用ax.quiver(lon, lat, u, v)这个函数有没有报错：IndexError: 2-dimensional boolean indexing is not supported. ？
我是用的NCEP数据
代码如下：
dsu=xr.open_dataset(uwind_dpath)
dsv=xr.open_dataset(vwind_dpath)

uwnd=dsu['uwnd']
vwnd=dsv['vwnd']


u=uwnd.loc['2019-01-01',850,:,:]
v=vwnd.loc['2019-01-01',850,:,:]

lat=dsu['lat']
lon=dsu['lon']

fig=plt.figure()
ax=fig.add_subplot(111,projection=ccrs.PlateCarree(180))
ax.quiver(lon,lat,u,v,transform=ccrs.PlateCarree())

咸墨

dafenqi 发表于 2021-2-14 19:43
请问你使用ax.quiver(lon, lat, u, v)这个函数有没有报错：IndexError: 2-dimensional boolean indexing  ...

ax.quiver(lon,lat,u,v,transform=ccrs.PlateCarree())
把这一句的transform参数去掉就行，改为
ax.quiver(lon,lat,u,v)

dafenqi

咸墨 发表于 2021-2-15 22:32
ax.quiver(lon,lat,u,v,transform=ccrs.PlateCarree())
把这一句的transform参数去掉就行，改为
ax.qui ...

嗯 谢谢您回复！ 我去试试！！

YanPenG

你可以看下我之前那个 http://bbs.06climate.com/forum.php?mod=viewthread&tid=60499   基于nc数据画的 温度 湿度 风向图

zhangmeiru

一大碗年糕 发表于 2020-12-12 16:38
今天也碰到这个报错，把xarray的U,V的值要先提取出来才行
lon = mveof['lon']
lat = mveof['lat']

对的，谢谢分享

hanwangyu

咸墨 发表于 2021-2-15 22:32
ax.quiver(lon,lat,u,v,transform=ccrs.PlateCarree())
把这一句的transform参数去掉就行，改为
ax.qui ...

非常感谢！正好遇到这个问题