Meteoinfolab Script Example: AMSR-E satellite data projection

The LAND3 data in the AMSR-E (http://nsidc.org/data/amsre/index.html) data is HDF-EOS4 format, and the projection is Cylindrical_equal_area. Here the sample reads the data and projects it to the latitude and longitude projection.

Script Program:

#Add Data Filefn ='AMSR_E_L3_DAILYLAND_V06_20091231.HDF'F= AddFile (Os.path.join ('D:/TEMP/HDF', FN))#vname = ' d_soil_moisture 'VName ='A_tb36.5h_ (res_1)'Data= F[vname][::-1,:]data[data==0] = 9999Data.fill_value=-9999Data= Data * 0.1yn=data.dimlen (0) Xn= Data.dimlen (1)#Project DataToproj = Projinfo ()#Longlat Projectionx = Linspace (-180, 180, xn) y= Linspace (-90, 90, yn) lon, lat=Meshgrid (x, y) data= Data.project (lon, LAT, toproj=toproj)#PlotAXESM (tickfontsize=12) Mlayer= Shaperead ('d:/temp/map/country1.shp') geoshow (Mlayer, Edgecolor='k') Levs= Arange (140, 320, 10) Layer=imshowm (x, y, data, Levs) Colorbar (layer, FontSize=12) Xticks (Arange (-180, 181, 30)) Yticks (Arange (-90, 91, 30) title ([FN, VName])

Meteoinfolab Script Example: AMSR-E satellite data projection