FY2E HDF格式数据处理绘图

圆盘标称投影数据时静止气象卫星常见的数据产品，比如FY2E静止气象卫星就有很多这样的产品（可以从国家卫星气象中心网站上下载）。所谓的圆盘标称投影就是Geostationary投影，主要的投影参数有中央经度和卫星高度。这里给出一个示例脚本程序（MeteoInfoLab环境）可以读取FY2E HDF格式数据、将数据投影为等经纬度并生成数据图层。需要最新的MeteoInfo Java版本（MeteoInfo_Java_1.2.2R1.zip）：http://yun.baidu.com/share/link?shareid=669776748&uk=51062435

脚本程序：

#Add data file
fn = 'D:/Temp/Hdf/FY2E_TBB_IR1_NOM_20150126_1700.hdf'
f = addfile(fn)
#Get data variable
v = f['FY2E_TBB_Hourly_Product']
#Set x/y
x = linspace(-5750000.0,5750000.0, 2288)
y = linspace(-5750000.0,5750000.0, 2288)
#Set x/y dimensions
v.setdim('Y', y)
v.setdim('X', x)
#Get data array
tbb = v[:,:]
#Plot
ax, proj = axesm(proj='geos', origin=(0,104.5,0), h=35785864, gridlabel=False)
mlayer = shaperead('D:/Temp/map/country1.shp')
geoshow(ax, mlayer)
layer = imshowm(ax, tbb, 20, proj=proj)
colorbar(layer)
show()

　　

运行结果：