读取和绘制海冰运行矢量

MeteoInfo

NSIDC的海冰运动矢量数据是netCDF格式，由于是极地数据，地图投影为lamber azimuthal equal area。MeteoInfoLab的addfile函数打开数据文件后会自动读取数据的投影信息，以及在此投影下数据的x/y坐标值。读取u, v变量的值，用magnitude函数计算风速，用quiver函数绘制风场矢量，绘图时对数据进行了稀疏化处理。

1. fn = 'D:/Temp/nc/icemotion_daily_nh_25km_20220101_20221231_v4.1.nc'
   
2. f = addfile(fn)
   
3. u = f['u'][0]
   
4. v = f['v'][0]
   
5. speed = np.magnitude(u, v)
   
6. 
   
7. axesm(projection=f.proj)
   
8. geoshow('country')
   
9. s = 4
   
10. qv = quiver(u[::s,::s], v[::s,::s], speed[::s,::s], 20, size=15,
    
11.     proj=f.proj)
    
12. quiverkey(qv, 0.93, 0.05, 10, bbox=dict(edge=True, fill=True))
    
13. colorbar()

fn = 'D:/Temp/nc/icemotion_daily_nh_25km_20220101_20221231_v4.1.nc'<br /> f = addfile(fn)<br /> u = f['u'][0]<br /> v = f['v'][0]<br /> speed = np.magnitude(u, v)<br /> <br /> axesm(projection=f.proj)<br /> geoshow('country')<br /> s = 4<br /> qv = quiver(u[::s,::s], v[::s,::s], speed[::s,::s], 20, size=15, <br />     proj=f.proj)<br /> quiverkey(qv, 0.93, 0.05, 10, bbox=dict(edge=True, fill=True))<br /> colorbar()