python ERA5 画水汽通量散度图地图:风速风向矢量图、叠加等高线、色彩分级、添加shp文件、添加位置点及备注




import matplotlib.pyplot as plt
import numpy as np
from metpy.units import units
import xarray as xr
from metpy.calc import dewpoint, theta
from metpy.plots import plot_crosshairs, plot_barb, plot_dove_tail_vorticity, plot_planet
 
# 假设已经有了ERA5的数据,并且存储在了netCDF文件中
era5_file = 'era5_data.nc'
 
# 读取数据
era5 = xr.open_dataset(era5_file)
 
# 选择一个特定的站点
lat = 40.0
lon = -105.0
 
# 选择需要的变量
u = era5['u'][0, :, :].values * units('m/s')  # 东西风
v = era5['v'][0, :, :].values * units('m/s')  # 南北风
t = era5['temperature'][0, :, :].values * units('K')  # 温度
td = era5['dewpoint_temperature'][0, :, :].values * units('K')  # 湿度
hgt = era5['geopotential'][0, :, :].values * units('m^2/s^2')  # 地球气压高度
 
# 计算相对湿度
rh = np.clip(dewpoint(t, td).to('degC'), 0, None)
 
# 计算水汽通量
evaporation = (
    6.112 * np.exp(5419. / (t - rh)) * (hgt / 9.8) ** 2
) * units('mm/day')  # 假设数据单位是mm/day
 
# 创建画图
fig = plt.figure(figsize=(12, 9))
ax = fig.add_subplot(1, 1, 1, projection=era5.metpy.cartopy_crs)
 
# 绘制等高线
ax.contourf(era5['level'].values, era5.metpy.y, era5.metpy.x, hgt.m, 15, cmap='Spectral_r')
 
# 绘制色彩分级的水汽通量
plot_planet(era5.metpy.y, era5.metpy.x, evaporation, cmap='viridis', norm=None, ax=ax)
 
# 绘制风速和风向
plot_barb(era5.metpy.y, era5.metpy.x, u, v, ax=ax)
 
# 绘制笛卡尔流向和笛卡尔强度
plot_dove_tail_vorticity(era5.metpy.y, era5.metpy.x, u, v, ax=ax)
 
# 绘制坐标
plot_crosshairs(lat, lon, ax=ax, linewidth=2)
 
# 设置标题和坐标轴
ax.set_title(f'ERA5 Wind, Vorticity, and Evaporation at {lat:.2f}°N, {lon:.2f}°W')
ax.set_xlabel('Longitude')
ax.set_ylabel('Latitude')
 
# 显示图例
plt.legend()
 
# 保存图片
plt.savefig('era5_map.png', bbox_inches='tight', pad_inches=0.1)
 
# 显示图像
plt.show()

这段代码假设您已经有了ERA5数据,并且存储在了一个名为era5_data.nc的netCDF文件中。代码读取了该文件,并选择了一个特定的站点来绘制风速、风向、笛卡尔流向和笛卡尔强度、水汽通量散度图以及站点坐标。最终保存了一个名为era5_map.png的图像文件。

最后修改于:2024年08月14日 08:47

评论已关闭

推荐阅读

DDPG 模型解析,附Pytorch完整代码
2024年11月24日
DQN 模型解析,附Pytorch完整代码
2024年11月24日
AIGC实战——Transformer模型
2024年12月01日
Socket TCP 和 UDP 编程基础(Python)
2024年11月30日
python , tcp , udp
如何使用 ChatGPT 进行学术润色?你需要这些指令
2024年12月01日
AI
最新 Python 调用 OpenAi 详细教程实现问答、图像合成、图像理解、语音合成、语音识别(详细教程)
2024年11月24日
ChatGPT 和 DALL·E 2 配合生成故事绘本
2024年12月01日
omegaconf,一个超强的 Python 库!
2024年11月24日
【视觉AIGC识别】误差特征、人脸伪造检测、其他类型假图检测
2024年12月01日
[超级详细]如何在深度学习训练模型过程中使用 GPU 加速
2024年11月29日
Python 物理引擎pymunk最完整教程
2024年11月27日
MediaPipe 人体姿态与手指关键点检测教程
2024年11月27日
深入了解 Taipy:Python 打造 Web 应用的全面教程
2024年11月26日
基于Transformer的时间序列预测模型
2024年11月25日
Python在金融大数据分析中的AI应用(股价分析、量化交易)实战
2024年11月25日
AIGC Gradio系列学习教程之Components
2024年12月01日
Python3 `asyncio` — 异步 I/O,事件循环和并发工具
2024年11月30日
llama-factory SFT系列教程:大模型在自定义数据集 LoRA 训练与部署
2024年12月01日
Python 多线程和多进程用法
2024年11月24日
Python socket详解,全网最全教程
2024年11月27日