Most of the pixels in thermal infrared remote sensing images are three-dimensional non-isothermal pixel, especially for the pixels with the size of meters, tens of meters or hundreds of meters which have received widespread attention in geoscience and remote sensing. Even though the sizes of some pixels reach centimeters, the three-dimensional non-isothermal phenomenon may still arise. So, it is very important to accurately determine the component temperatures in one pixel for the related researches in geoscience. The remote sensing data used to carry out the related inversion experiments in this paper was the airborne remote sensing data obtained by WSIS (Wide Spectrum Imaging Spectrometer) the imaging wave bands of which include VNIR (visible light and near infrared), SWIR (short wave infrared) and TIR (thermal infrared). Firstly, the components of all the pixels in the image were determined through the VNIR images using linear mixing spectral model. Secondly, the emissivity of each component in every pixel in the image was determined according to a prior knowledge base of emissivity of many surface features. Thirdly, the so called average temperature of every pixel was retrieved using the TES (temperature and emissivity separation) algorithm. The retrieved temperature was regarded as initial value. The multi-band equations were established after the linearization of Planck function, and the component temperatures of every pixel in the image were inversed. The results show that the accuracy of the component temperatures inversion in one pixel can be improved obviously, with the combination of the VNIR, SWIR and TIR images.