MeteoInfoLab脚本示例：TRMM 2A12 HDF数据

TRMM 2A12 HDF数据是卫星观测的SWATH数据（轨道数据），比格点数据处理起来要麻烦一些。数据的经纬度保存在geolocation变量中，需要先将经纬度数据读出来（均为2维数组），然后读取云水含量数据（cldWater）。虽然都是2维数组，但并不是格点数据，相当于2维的散点数据，点数相当多（3019*208）如果用散点图来绘制的话会非常的慢，需要将其插值为格点数据。插值方法建议选择nearest，该方法速度最快。插值后的格点数据用imshowm函数显示为图像就很快了。

脚本程序：

#Add data file
folder = 'D:/Temp/hdf/'
fns = '2A12_090711_66394_6.HDF'
fn = folder + fns
f = addfile(fn)
#Get data variable
geo = f['geolocation']
lat = geo[:,:,0]
lon = geo[:,:,1]
water = f['cldWater'][:,:,9]
water = water / 1000
#Interpolate to grid data
gwater, x_g, y_g = griddata((lon, lat), water, method='neareast')
#Plot
axesm()
mlayer = shaperead('D:/Temp/map/country1.shp')
geoshow(mlayer)
levs = arange(0, 10, 0.5)
layer = imshowm(x_g, y_g, gwater, levs)
colorbar(layer, orientation='horizontal')
title(fns + ' (cldWater level=9)')
axism([-180, 180, -90, 90])

转载于:https://www.cnblogs.com/yaqiang/p/4622831.html