MeteoInfoLab脚本示例：AIRS Swath HDF数据

MeteoInfo

例子中的AIRS Swath HDF数据在Polar Stereographic（南极）投影中接近矩形，需要先从数据中读出经纬度及相关数据数组，利用surfacem函数绘制Swath数据（散点），在surfacem函数中经纬度数据会被投影到目的投影坐标系（在axesm函数中定义），并从2维散点数据生成surface数据（也就是格点数据）用所谓的surface插值方法。2维散点数据可以组成格网（非矩形格网），遍历目标格点数据的每个格点，该格点落在格网中哪个网格中，就将此网格的散点值赋给该格点。有些拗口，其实算法不是很复杂，和MatLab中的surfacem函数的功能类似。

脚本程序：

1. #Add data file
   
2. folder = 'D:/Temp/hdf/'
   
3. fns = 'AIRS.2002.12.31.001.L2.CC_H.v4.0.21.0.G06100185050.hdf'
   
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 = 'radiances'
    
11. rad_v = f[vname]
    
12. rad = rad_v[:,:,567]
    
13. #Plot
    
14. plot,proj = axesm(proj='stere', lat_0=-90, gridline=True, griddx=30, griddy=30)
    
15. mlayer = shaperead('D:/Temp/map/country1.shp')
    
16. geoshow(mlayer)
    
17. levs = arange(40, 90, 1)
    
18. #slayer = scatterm(lon, lat, rad, levs, edge=True)
    
19. layer = surfacem(lon, lat, rad, levs)
    
20. colorbar(layer, orientation='horizontal')
    
21. title('{0}\n {1}'.format(fns, vname))
    
22. axism()

#Add data file<br /> folder = 'D:/Temp/hdf/'<br /> fns = 'AIRS.2002.12.31.001.L2.CC_H.v4.0.21.0.G06100185050.hdf'<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 = 'radiances'<br /> rad_v = f[vname]<br /> rad = rad_v[:,:,567]<br /> #Plot<br /> plot,proj = axesm(proj='stere', lat_0=-90, gridline=True, griddx=30, griddy=30)<br /> mlayer = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(mlayer)<br /> levs = arange(40, 90, 1)<br /> #slayer = scatterm(lon, lat, rad, levs, edge=True)<br /> layer = surfacem(lon, lat, rad, levs)<br /> colorbar(layer, orientation='horizontal')<br /> title('{0}\n {1}'.format(fns, vname))<br /> axism()