pip install numpy matplotlib xarray

• 1.

import xarray as xr

# 假设数据文件名为wind_data.nc
ds = xr.open_dataset('wind_data.nc')

• 1.
• 2.
• 3.
• 4.

import numpy as np

# 假设风速数据存储在ds['wind_speed']，风向数据存储在ds['wind_direction']
wind_speed = ds['wind_speed']
wind_direction = ds['wind_direction']

# 将风向转换为弧度
wind_direction_rad = np.deg2rad(wind_direction)

# 计算U分量
U = -wind_speed * np.sin(wind_direction_rad)

# 计算V分量
V = wind_speed * np.cos(wind_direction_rad)

# 将U、V分量添加到数据集中
ds['U'] = U
ds['V'] = V

• 1.
• 2.
• 3.
• 4.
• 5.
• 6.
• 7.
• 8.
• 9.
• 10.
• 11.
• 12.
• 13.
• 14.
• 15.
• 16.
• 17.
• 18.

import matplotlib.pyplot as plt

# 绘制U分量
plt.figure(figsize=(8, 6))
plt.contourf(ds['longitude'], ds['latitude'], ds['U'], cmap='coolwarm')
plt.colorbar(label='U Component (m/s)')
plt.title('U Component of Wind Field')
plt.xlabel('Longitude')
plt.ylabel('Latitude')
plt.show()

# 绘制V分量
plt.figure(figsize=(8, 6))
plt.contourf(ds['longitude'], ds['latitude'], ds['V'], cmap='coolwarm')
plt.colorbar(label='V Component (m/s)')
plt.title('V Component of Wind Field')
plt.xlabel('Longitude')
plt.ylabel('Latitude')
plt.show()

• 1.
• 2.
• 3.
• 4.
• 5.
• 6.
• 7.
• 8.
• 9.
• 10.
• 11.
• 12.
• 13.
• 14.
• 15.
• 16.
• 17.
• 18.
• 19.

# 导出数据集到NetCDF文件
ds.to_netcdf('uv_wind_field.nc')

• 1.
• 2.