MeteoInfoLab脚本示例:AVHRR HDF数据

这里演示读取和绘制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()

 

posted on 2015-09-28 16:01  yaqiang  阅读(323)  评论(0)    收藏  举报

导航