请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 7416|回复: 7

[求助] python 处理ERA5数据并画风场

[复制链接]

新浪微博达人勋

发表于 2022-7-4 20:51:28 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 menseye 于 2022-7-5 10:44 编辑

我用python处理ERA5数据,想画出风场,该数据是2021年10月30日18时950百帕的u和v数据,我感觉思路没问题,分别提取u,v分量,再利用bars画出风羽图,运行程序后没有报错,但是只出来空白的地图,请教各位大神,应该怎么修改?
十分感激各位大佬相助!!!

以下为程序
import numpy as np
from datetime import datetime
import xarray as xr
from scipy.ndimage import gaussian_filter
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cfeature
import metpy.calc as mpcalc
from metpy.units import units
import  netCDF4 as nc
from cartopy.mpl.gridliner import LATITUDE_FORMATTER, LONGITUDE_FORMATTER


ds=nc.Dataset(r'D:\Python_total\DATA\EC_data\adaptor.mars.internal-1656849067.1887913-23322-3-5e39ab33-abfe-4988-abe4-5b45d8bad5c4.nc')
time=ds.variables['time'][:]

lat=ds.variables['latitude'][:]
lon=ds.variables['latitude'][:]           #经评论区大佬提醒,这里经度提取错误,再次感谢评论区各位大佬批评指正!
lons,lats=np.meshgrid(lon,lat)
u=ds.variables["u"][0][:][:]#仅有一个时次,所以用0代表选用的时次
v=ds.variables["v"][0][:][:]
print(u)

print(ds.variables['time'][:])

uwind=u.data*units('m/s')
vwind=v.data*units('m/s')
sped=mpcalc.wind_speed(uwind,vwind)
#sped0800=sped.sel(time=1067888)

datacrs=ccrs.PlateCarree()
fig=plt.figure(figsize=(9,9))
ax=fig.add_axes([0.05,0.05,0.9,0.9],projection=datacrs)
rejion=[110,130,30,50]
ax.set_extent(rejion,datacrs)


ax.barbs(lons, lats, uwind, vwind,
         pivot='middle',color='black', regrid_shape=20, transform=datacrs)
ax.add_feature(cfeature.COASTLINE.with_scale('50m'))
ax.add_feature(cfeature.OCEAN.with_scale('50m'))
ax.add_feature(cfeature.LAND.with_scale('50m'))

plt.show()

程序出图

程序出图


                               
登录/注册后可看大图


程序出图,无报错

程序出图,无报错
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2022-7-5 15:43:16 | 显示全部楼层
请问楼主下载的是ERA5 hourly data on pressure levels from 1959 to present在950hpa上的资料吗?
密码修改失败请联系微信:mofangbao
回复 支持 0 反对 1

使用道具 举报

新浪微博达人勋

发表于 2022-7-5 08:50:00 | 显示全部楼层
先画海岸线,陆地特征那些,之后再画风羽图,你那个add_feature直接覆盖了风羽
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-7-5 09:12:02 | 显示全部楼层
lat=ds.variables['latitude'][:]
lon=ds.variables['latitude'][:]
lon没有正确赋值
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2022-7-5 10:28:18 | 显示全部楼层
cookie-o-o 发表于 2022-7-5 09:12
lat=ds.variables['latitude'][:]
lon=ds.variables['latitude'][:]
lon没有正确赋值

感谢大佬,修改之后就画出风场了。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2022-7-5 10:31:22 | 显示全部楼层
千与千叶 发表于 2022-7-5 08:50
先画海岸线,陆地特征那些,之后再画风羽图,你那个add_feature直接覆盖了风羽

我把lon的赋值修改了就出来风场了,所以应该不是这个原因
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2022-7-7 12:41:21 | 显示全部楼层
dj11235 发表于 2022-7-5 15:43
请问楼主下载的是ERA5 hourly data on pressure levels from 1959 to present在950hpa上的资料吗?

是的,就是这个资料;
另,鼠标误触,点了反对,抱歉!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-9-1 17:27:05 | 显示全部楼层
menseye 发表于 2022-7-7 12:41
是的,就是这个资料;
另,鼠标误触,点了反对,抱歉!

好的谢谢啦~
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表