MeteoInfoLab脚本示例:AIRS Grid HDF数据

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

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