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)
放大到数据范围：
MeteoInfoLab脚本示例：MODIS Sinusoidal投影HDF数据
原文：http://www.cnblogs.com/yaqiang/p/4652389.html