MeteoInfoLab脚本示例：MODIS Sinusoidal投影HDF数据

MODIS卫星很多陆面数据都是Sinusoidal投影，数据被分为一个个10*10度（赤道地区）的瓦片（http://modis-land.gsfc.nasa.gov/MODLAND_grid.html），数据格式通常是HDF-EOS格式，作为EOS格式MeteoInfo可以自动读出数据的坐标信息（包括投影等）。

脚本程序：

f = addfile('D:/Temp/Hdf/MOD10A1.A2010365.h27v05.005.2011002103013.hdf')
vname = 'snow_Albedo_Daily_Tile'
data = f[vname][:,:]
axesm(projinfo=f.proj, griddx=5, griddy=5, gridline=True)
mlayer = shaperead('D:/Temp/map/country1.shp')
geoshow(mlayer)
levs = arange(0, 100, 5)
layer = imshowm(data, 20, proj=f.proj)
title(vname)
colorbar(layer)
axism()

 

放大到数据范围：

 

在创建axesm坐标系是不给出投影参数，缺省是经纬度投影，数据图形会自动投影。

axesm(griddx=30, griddy=30, gridline=True)

 

放大到数据范围：