圆盘标称投影数据时静止气象卫星常见的数据产品,比如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()
运行结果:
手机扫一扫
移动阅读更方便
你可能感兴趣的文章