In order to understand the processes of land surface-atmosphere interaction over de-sertification area, it is indispensable to utilize of satellite remote sensing. Two scenes of LandsatTM were used to produce a set of maps of surface reflectance, MSAVI, vegetation coverage, sur-face temperature, net radiation, soil heat flux, sensible heat flux and latent heat flux. Statisticalanalysis based on these maps revealed some quantitative significant land surface characteristics.Future developments of the method are also discussed.