MeteoInfoLab脚本示例：AMSR-E卫星数据投影

AMSR-E（http://nsidc.org/data/amsre/index.html）数据中的Land3数据是HDF-EOS4格式，投影是Cylindrical_Equal_Area。这里示例读取数据并投影至等经纬度投影。

脚本程序：

#Add data file
fn = '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] = -9999
data.fill_value = -9999
data = data * 0.1
yn = data.dimlen(0)
xn = data.dimlen(1)
#Project data
toproj = projinfo()    #longlat projection
x = linspace(-180, 180, xn)
y = linspace(-90, 90, yn)
lon, lat = meshgrid(x, y)
data = data.project(lon, lat, toproj=toproj)
#Plot
axesm(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])