MeteoInfoLab脚本示例：OMI Swath HDF数据

MeteoInfo

这个例子读取OMI卫星Swath数据中的CloudFaction变量并绘图。

脚本程序：

1. #Add data file
   
2. folder = 'D:/Temp/hdf/'
   
3. fns = 'OMI-Aura_L2-OMNO2_2008m0720t2016-o21357_v003-2008m0721t101450.he5'
   
4. fn = folder + fns
   
5. f = addfile(fn)
   
6. lon_v = f['Longitude']
   
7. lat_v = f['Latitude']
   
8. lon = lon_v[:,:]
   
9. lat = lat_v[:,:]
   
10. vname = 'CloudFraction'
    
11. v = f[vname]
    
12. data = v[:,:]*0.001
    
13. #Plot
    
14. axesm()
    
15. mlayer = shaperead('D:/Temp/map/country1.shp')
    
16. geoshow(mlayer)
    
17. #layer = scatterm(lon, lat, data, edge=False)
    
18. layer = surfacem(lon, lat, data)
    
19. colorbar(layer)
    
20. title('OMI - ' + vname)
    
21. axism()

#Add data file <br /> folder = 'D:/Temp/hdf/' <br /> fns = 'OMI-Aura_L2-OMNO2_2008m0720t2016-o21357_v003-2008m0721t101450.he5' <br /> fn = folder + fns <br /> f = addfile(fn) <br /> lon_v = f['Longitude'] <br /> lat_v = f['Latitude'] <br /> lon = lon_v[:,:] <br /> lat = lat_v[:,:] <br /> vname = 'CloudFraction' <br /> v = f[vname] <br /> data = v[:,:]*0.001 <br /> #Plot <br /> axesm() <br /> mlayer = shaperead('D:/Temp/map/country1.shp') <br /> geoshow(mlayer) <br /> #layer = scatterm(lon, lat, data, edge=False) <br /> layer = surfacem(lon, lat, data) <br /> colorbar(layer)<br /> title('OMI - ' + vname) <br /> axism()

shan

很好，大牛啊!

nowayback

好漂亮的图！

狼牙李

感谢楼主分享