Temperature and emissivity retrieval from remotely sensed images using the ''grey body emissivity'' method

The problem of extracting temperature and emissivity information from thermal infrared multispectral radiance data is of great importance in remote sensing applications. The Earth surface temperature is a relevant parameter for climatology, meteorology, oceanography, fire mapping, wind field determination, and energy balance at the atmosphere-Earth interface. In addition, ground emissivity is very important for lithological mapping, geological inspection, petrography, the dating of lava hows, etc. A new method for emissivity and temperature retrieval has been developed in our laboratory which is based on the assumption that emissivity is a slowly varying function of wavelength. It is shown that this assumption holds true for many kinds of natural targets. The resulting model is called ''grey body emissivity,'' and may be used to retrieve the emissivity of soils and other natural materials. The problem with using an effective wavelength or a bandwidth integral in order to represent the sampled radiance is discussed, thus leading to the implementation of two algorithms for the grey body emissivity. It is demonstrated that one of them, extensively discussed throughout the paper, provides the maximum likelihood values for temperature and emissivity. The other reveals two drawbacks: it produces noisy processed images and works well only for target temperatures greater than 274 K. Comparisons with other models are carried out utilizing remotely sensed images, laboratory measurements, and numerical simulations. The main approximations and the resulting errors, contained in the considered models, are briefly discussed. The influence of the atmospheric effects on utilized remote sensing data is taken into account. Finally, we investigate the emissivity dependence on the viewing angle. Theoretical predictions which rely on Fresnel's and Kirchhoff's laws are compared with the emissivity maps extracted from TIMS images over sea. Starting from the image data alone, the procedure allows us to compute the refraction index and the absorption coefficient of the concerned marine water.