这里演示读取和绘制AVHRR hdf格式数据,以sst(海表面温度)为例。
脚本程序:
#Add data file f = addfile('D:/Temp/hdf/2006001-2006005.s0454pfrt-bsst.hdf') vname = 'bsst' v = f[vname] data = v[::-1,:] data.fill_value = 0.0 #Plot axesm() mlayer = shaperead('D:/Temp/map/country1.shp') geoshow(mlayer) layer = imshowm(data, 20) colorbar(layer) title('AVHRR - ' + vname) axism()