After having analyzed the requirement on the aerodynamic earth's surface roughness in two-dimensional distribution in the research field of interaction between land surface and atmosphere, this paper presents a new way to calculate the aerodynamic roughness using the earth's surface geometric roughness retrieved from SAR (Synthetic Aperture Radar) and TM thermal infrared image data. On the one hand, the SPM (Small Perturbation Model) was used as a theoretical SAR backscattering model to describe the relationship between the SAR backscattering coefficient and the earth's surface geometric roughness and its dielectric constant retrieved from the physical model between the soil thermal inertia and the soil surface moisture with the simultaneous TM thermal infrared image data and the ground microclimate data. On the basis of the SAR image matching with the TM image, the non-volume scattering surface geometric information was obtained from the SPM model at the TM image pixel scale, and the ground pixel surface's equivalent geometric roughness-height standard RMS (Root Mean Square) was achieved from the geometric information by the transformation of the typical topographic factors. The vegetation (wheat, tree) height retrieved from spectrum model was also transferred into its equivalent geometric roughness. A completely two-dimensional distribution map of the equivalent geometric roughness over the experimental area was produced by the data mosaic technique. On the other hand, according to the atmospheric eddy currents theory, the aerodynamic surface roughness was iterated out with the atmosphere stability correction method using the wind and the temperature profiles data measured at several typical fields such as bare soil field and vegetation field. After having analyzed the effect of surface equivalent geometric roughness together with dynamic and thermodynamic factors on the aerodynamic surface roughness within the working area, this paper first establishes a scale transformation model to calculate the aerodynamic surface roughness from surface equivalent geometric roughness. The final result retrieved from above series of models was validated using the filed measurements. It is concluded that the proposed approach is operational and feasible to derive the aerodynamic surface roughness at the pixel scale from a combination of the SAR and TM images.