使用Matplotlib Quiver从NetCDF文件绘制风向异常求助

你遇到的quiver显示网格而非箭头的问题，大概率是数据点过于密集或参数设置不匹配导致的，下面是具体的解决方法：

常见原因及对应解决方案 #

1. 数据点太密集，箭头重叠成网格 #

气象网格化数据的分辨率通常很高（比如0.25°、1°间隔），直接绘制所有点会让箭头挤在一起，看起来像网格。解决方法是抽稀数据，每隔N个点取一个：

# 设置抽稀步长，可根据数据分辨率调整
skip = 5

# 对经纬度和U/V分量进行抽稀
lons_sub = lons[::skip, ::skip]
lats_sub = lats[::skip, ::skip]
ucomp_sub = ucomp[::skip, ::skip]
vcomp_sub = vcomp[::skip, ::skip]

2. 经纬度数组维度不匹配 #

如果你的lons和lats是1D数组（比如从NetCDF中直接读取的一维坐标），需要先转换成2D网格，否则quiver无法正确映射位置：

import numpy as np

# 检查并转换维度
if lons.ndim == 1 and lats.ndim == 1:
    lons_grid, lats_grid = np.meshgrid(lons, lats)
else:
    lons_grid, lats_grid = lons, lats

# 结合抽稀步骤处理网格数据
lons_sub = lons_grid[::skip, ::skip]
lats_sub = lats_grid[::skip, ::skip]

3. 调整quiver显示参数，让箭头更清晰 #

默认的箭头大小、密度不适合高密度数据，手动设置scale（控制箭头长度）、width（箭头宽度）、headwidth（箭头头部大小）参数：

import matplotlib.pyplot as plt

# 绘制抽稀后的风向箭头
plt.quiver(lons_sub, lats_sub, 
           ucomp_sub, vcomp_sub,
           scale=100,  # 值越大箭头越短，根据实际风速调整
           width=0.002,  # 箭头宽度，避免过粗
           headwidth=3)  # 箭头头部宽度，让箭头形状更明显

# 添加坐标轴标签和标题
plt.xlabel('经度')
plt.ylabel('纬度')
plt.title('10米风向分布')
plt.show()

进阶优化：结合风速上色 #

可以把风速作为颜色映射，让可视化更直观：

# 计算抽稀后的风速
net_speed_sub = np.sqrt(ucomp_sub**2 + vcomp_sub**2)

# 绘制带颜色映射的风向箭头
q = plt.quiver(lons_sub, lats_sub, 
               ucomp_sub, vcomp_sub, net_speed_sub,
               scale=100, width=0.002, headwidth=3,
               cmap='viridis')

# 添加颜色条
plt.colorbar(q, label='风速 (m/s)')
plt.xlabel('经度')
plt.ylabel('纬度')
plt.title('10米风向与风速分布')
plt.show()

内容的提问来源于stack exchange，提问作者Samman Amgain