MeteoInfoLab脚本示例：AVHRR HDF数据

MeteoInfo

这里演示读取和绘制AVHRR hdf格式数据，以sst（海表面温度）为例。

脚本程序：

1. f = addfile('D:/Temp/hdf/2006001-2006005.s0454pfrt-bsst.hdf')
   
2. vname = 'bsst'
   
3. v = f[vname]
   
4. data = v[::-1,:]
   
5. data.fill_value = 0.0
   
6. #Plot
   
7. axesm()
   
8. mlayer = shaperead('D:/Temp/map/country1.shp')
   
9. geoshow(mlayer)
   
10. layer = imshowm(data, 20)
    
11. colorbar(layer)
    
12. title('AVHRR - ' + vname)
    
13. axism()

f = addfile('D:/Temp/hdf/2006001-2006005.s0454pfrt-bsst.hdf')<br /> vname = 'bsst'<br /> v = f[vname]<br /> data = v[::-1,:]<br /> data.fill_value = 0.0<br /> #Plot<br /> axesm()<br /> mlayer = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(mlayer)<br /> layer = imshowm(data, 20)<br /> colorbar(layer)<br /> title('AVHRR - ' + vname)<br /> axism()

梁的丰

真好，给楼主点赞...

qdsonglinyy

MARK，绘制近海海域闪电密度分布时应该用得上