圆盘标称投影数据时静止气象卫星常见的数据产品,比如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()
运行结果: