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])
