AIRS Grid HDF数据是HDF4 EOS格式,数据地理坐标信息可以被MeteoInfo自动识别,脚本程序更为简单。需要注意的是读取数据时Y轴是反向的(卫星数据通常如此)。
脚本程序:
#Add data file folder = 'D:/Temp/hdf/' fns = 'AIRS.2002.08.01.L3.RetStd_H031.v4.0.21.0.G06104133732.hdf' fn = folder + fns f = addfile(fn) vname = 'Temperature_MW_A' t_v = f[vname] t = t_v[11,::-1,:] #Plot axesm() mlayer = shaperead('D:/Temp/map/country1.shp') geoshow(mlayer) layer = imshowm(t, 20) colorbar(layer, orientation='horizontal') title(vname + 'at TempPrsLvls=11') axism()