MeteoInfoLab脚本示例：AIRS Grid HDF数据

MeteoInfo

AIRS Grid HDF数据是HDF4 EOS格式，数据地理坐标信息可以被MeteoInfo自动识别，脚本程序更为简单。需要注意的是读取数据时Y轴是反向的（卫星数据通常如此）。

脚本程序：

1. #Add data file
   
2. folder = 'D:/Temp/hdf/'
   
3. fns = 'AIRS.2002.08.01.L3.RetStd_H031.v4.0.21.0.G06104133732.hdf'
   
4. fn = folder + fns
   
5. f = addfile(fn)
   
6. vname = 'Temperature_MW_A'
   
7. t_v = f[vname]
   
8. t = t_v[11,::-1,:]
   
9. #Plot
   
10. axesm()
    
11. mlayer = shaperead('D:/Temp/map/country1.shp')
    
12. geoshow(mlayer)
    
13. layer = imshowm(t, 20)
    
14. colorbar(layer, orientation='horizontal')
    
15. title(vname + 'at TempPrsLvls=11')
    
16. axism()

#Add data file <br /> folder = 'D:/Temp/hdf/' <br /> fns = 'AIRS.2002.08.01.L3.RetStd_H031.v4.0.21.0.G06104133732.hdf' <br /> fn = folder + fns <br /> f = addfile(fn) <br /> vname = 'Temperature_MW_A' <br /> t_v = f[vname] <br /> t = t_v[11,::-1,:] <br /> #Plot <br /> axesm() <br /> mlayer = shaperead('D:/Temp/map/country1.shp') <br /> geoshow(mlayer) <br /> layer = imshowm(t, 20) <br /> colorbar(layer, orientation='horizontal') <br /> title(vname + 'at TempPrsLvls=11') <br /> axism()

eeenene

王老师这个怎么提取对应站点数据