Meteoinfolab Script Example: MODIS sinusoidal projection HDF Data

Many of the terrestrial data of the MODIS satellite are sinusoidal projections, and the data is divided into 10*10 (equatorial area) tiles (http://modis-land.gsfc.nasa.gov/MODLAND_ grid.html), the data format is usually the Hdf-eos format, and as the EOS format Meteoinfo can automatically read the data's coordinate information (including projection, etc.).

Script Program:

f = addfile ('d:/temp/hdf/mod10a1. A2010365.H27V05.005.2011002103013.HDF'snow_albedo_daily_tile' = F[vname][:,:]axesm (projinfo=f.proj, griddx=5, griddy=5, gridline== Shaperead (  ' d:/temp/map/country1.shp '  = Arange (0, 5= IMSHOWM (data, proj=f.proj) title (vname) Colorbar (layer) axism () 

Zoom to data range:

When creating the Axesm coordinate system, the projection parameter is not given, and the data graph is projected automatically by the latitude projection.

Axesm (griddx=30, griddy=30, Gridline=true)

Zoom to data range:

Meteoinfolab Script Example: MODIS sinusoidal projection HDF Data