MeteoInfoLab脚本示例：OMI Swath HDF数据

这个例子读取OMI卫星Swath数据中的CloudFaction变量并绘图。

脚本程序：

#Add data file
folder = 'D:/Temp/hdf/'
fns = 'OMI-Aura_L2-OMNO2_2008m0720t2016-o21357_v003-2008m0721t101450.he5'
fn = folder + fns
f = addfile(fn)
lon_v = f['Longitude']
lat_v = f['Latitude']
lon = lon_v[:,:]
lat = lat_v[:,:]
vname = 'CloudFraction'
v = f[vname]
data = v[:,:]*0.001
#Plot
axesm()
mlayer = shaperead('D:/Temp/map/country1.shp')
geoshow(mlayer)
#slayer = scatterm(lon, lat, data, edge=False)
layer = surfacem(lon, lat, data)
colorbar(layer)
title('OMI - ' + vname)
axism()

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