HY-2B散射计风场数据读取及绘图

MeteoInfo

HY-2B卫星搭载了一台Ku波段微波散射计（HSCAT-B），主要用于全球海面风场测量。数据是hdf5格式，下面的示例脚本用到了经纬度变量（wvc_lon, wvc_lat）和风速风向变量（wind_speed_selection, wind_dir_selection）。

1. fn = 'D:/Temp/FY/H2B_OPER_SCA_L2B_OR_20210819T225905_20210820T004328_14133_pwp_250_07_owv.h5'
   
2. f = addfile(fn)
   
3. wdir = f['wind_dir_selection'][:]
   
4. wspeed = f['wind_speed_selection'][:]
   
5. wdir[wdir<0] = nan
   
6. wspeed[wspeed<0] = nan
   
7. lon = f['wvc_lon'][:]
   
8. lat = f['wvc_lat'][:]
   
9. lon[lon>1000] = nan
   
10. lat[lat>1000] = nan
    
11. 
    
12. geoshow('country')
    
13. barbs(lon, lat, wdir, wspeed, wspeed, isuv=False)
    
14. colorbar(shrink=0.8)
    
15. xlim(-96, -80)
    
16. ylim(18, 30)

fn = 'D:/Temp/FY/H2B_OPER_SCA_L2B_OR_20210819T225905_20210820T004328_14133_pwp_250_07_owv.h5'<br /> f = addfile(fn)<br /> wdir = f['wind_dir_selection'][:]<br /> wspeed = f['wind_speed_selection'][:]<br /> wdir[wdir<0] = nan<br /> wspeed[wspeed<0] = nan<br /> lon = f['wvc_lon'][:]<br /> lat = f['wvc_lat'][:]<br /> lon[lon>1000] = nan<br /> lat[lat>1000] = nan<br /> <br /> geoshow('country')<br /> barbs(lon, lat, wdir, wspeed, wspeed, isuv=False)<br /> colorbar(shrink=0.8)<br /> xlim(-96, -80)<br /> ylim(18, 30)

贫道敬孔

不明觉厉，给王老师赞一个

莫莫泡沫

感谢，还在学习中