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()

原文地址:https://www.cnblogs.com/yaqiang/p/4844228.html