1. Introduction
The land surface temperature (LST) and land surface emissivity (LSE) are key parameters for many areas such as mineral identification [
1,
2], gas plume detection [
3], plant species [
4], soil moisture retrieval [
5]. The emergence of many hyperspectral thermal infrared (HTIR) imaging spectrometers, which have hundreds of contiguous spectral bands in the atmospheric window of 8–14 μm with spectral resolution higher than 100 nm, especially the commercial ones, has helped move thermal infrared hyperspectral remote sensing from concept to reality. We can retrieve simultaneously the LST and LSE of one pixel from the HTIR data with a temperature and emissivity separation (TES) algorithm. With the development of thermal infrared technology, many HTIR imaging spectrometers have appeared, for example, SEBASS [
6], MAKO [
7], Hyper-Cam LW [
8], Sieleters B3 [
9], hyperspectral thermal emission spectrometer (HyTES) [
10], AisaOWL [
11] and the airborne thermal infrared hyperspectral imaging system (ATHIS) [
12].
Different from single-band, dual-band or multispectral thermal infrared remote sensing which focus on the retrieval of LST, HTIR remote sensing is able to retrieve the LST and LSE simultaneously with higher quantitative inversion accuracy, which a way for an important application mode of HTIR data, namely using the retrieval LSE spectrum to identify ground objects. The spectral information in the LSE reflects the pattern characteristics of an object, which can also be understood as the changing characteristics of the emissivity curve. However, most ground objects have high emissivity and their emissivity changes are relatively small (for example, the average emissivity of water is about 0.98, and its fluctuation is only 0.02.) Such small changes are general characteristics of low-reflection ground objects in reflective optical remote sensing. On the one hand, it naturally comes to mind to improve the performance of the HTIR instrument when we want to obtain high-quality emissivity curves from HTIR data to describe the emission features of ground objects, or even to build connection between the emissivity curve of an object and its physical properties. On the other, from the principle of thermal infrared imaging, it is known that the spectral response accuracy and noise level of the instrument directly affect the quality of the acquired data and further affect HTIR inversion results. Therefore, it is necessary to study the quantitative influence of spectral response change (namely center wavelength shift and FWHM change) and instrument noise level on the LST and LSE inversion, which will provide significant support for the instrument design and development and the HTIR applications.
The difficulty of inversing the LST and LSE simultaneously from the hyperspectral or multispectral thermal infrared data observed in the same phase, namely separating the temperature and emissivity, lies in solving
N equations (one sensor output radiance spectrum with
N bands) to get
N + 1 unknowns (
N emissivities and one temperature), which is an underdetermined equation problem. It is necessary to introduce additional conditions in order to solve the problem, usually increasing the number of equations by proposing assumptions about the target emissivity spectrum shape or target temperature. There are several existing TES algorithms provided for multispectral thermal infrared data, such as reference band method (MBR) [
13], normalized emissivity method [
14,
15], alpha residues method [
16], temperature-independent spectral indices (TISI) [
17], spectral ratio method [
18], maximum-minimum difference (MMD) [
19], split-window [
20], or grey body emissivity [
21]. Compared with single-band or multispectral thermal infrared data, the hyperspectral thermal infrared data can provide more hypothetical conditions in line with real physical properties of objects for TES, effectively increasing the stability of the underdetermined equations to be solved. Most of existing TES algorithms for hyperspectral thermal infrared data are based on the well-known observation that the emissivity spectra of most materials are much smoother than the atmospheric transmission features (smoothness assumption). Based on cost function representations, there are spectral smoothing index class algorithms, e.g., the iterative spectrally smooth temperature-emissivity separation (ISSTES) [
22] and its improved version, the automatic retrieval of temperature and emissivity using spectral smoothness (ARTEMISS) [
23]; downwelling radiance residual class algorithms, e.g., the downwelling radiance residual index (DRRI) [
24], the correlation-based temperature and emissivity separation (CBTES) [
25], stepwise refining algorithm of temperature and emissivity separation (SRATES) [
26]; linear constraint class, e.g., linear spectral emissivity constraint (LSEC) [
27], improved LSEC (I-LSEC) [
28], pre-estimate shape LSEC (PES-LSEC) [
29]; and wavelet class algorithms, e.g., wavelet transform method for separating temperature and emissivity (WTTES) [
30], multi-scale wavelet-based temperature and emissivity separation (MSWTES)[
31].
As many TES more algorithms are described, more articles are devoted to the sensitivity analysis of TES algorithms. The impact factors of a TES algorithm performance are: assumptions reasonableness, robustness, atmospheric correction accuracy, spectral calibration accuracy, instrument noise and others [
32,
33,
34,
35,
36]. Among the literature [
33] and [
36] involve the influence analysis of spectral calibration accuracy and instrument noise on the TES algorithms. In [
33], an error analysis for the ARTEMISS algorithm shows that spectral shifts of more than 1/20th of the spectral sampling can produce fitting errors of the order of the sensor noise and errors in the FWHM of up to 10% still yield useful results. In [
36], the researchers studied the influence of the instrument spectral properties (the center wavenumber shifts were 4%, 20%, 40%, 60% and FWHM were widen by 5%, 10%, 20%) and noise (Gaussian random, 0.1 K, 0.2 K) on the retrieval results for five TES algorithms, and the results shows that even if some spectral response changes cause only a slight change in the observed radiance or the bright temperature spectral data they have a great impact on the LST and LSE retrieval. The above studies both have shown that the changes in spectral response characteristics have a great impact on the TES inversion results, however these related studies have paid more attention to exploring the application boundary conditions of the newly proposed TES algorithm or to prove the performance of the new method, moreover the number of samples are small, for example 54 in [
36], and the emissivity values of the samples in their experiments are high. Although TES algorithm improvements are helpful to improve the application of HTIR data under certain conditions, the advancement in the manufacturing level of HTIR instruments and quantification level of HTIR data can better promote the development of HTIR remote sensing from the source. Therefore, from the perspective of instrument manufacturing, we analyzed quantitatively the influence of the spectral parameters change and noise of an instrument on the TES inversion results, in order to provide a reference for HTIR instruments manufacturing and provide guidance on the application scope for a HTIR data with a given spectral parameters.
In terms of the instrument development, this article selected three instruments with currently common spectral resolutions (100 nm, 50 nm, and 35 nm) as representative examples, simulated tens of millions of radiances of 1525 emissivity curves with 17 center wavelength shift ratios, 6 FWHM change ratios and 6 noise equivalent differential temperatures (NEDTs) under 15 atmospheric conditions with 6 object temperatures, inverted them with two TES methods (ISSTES [
22] and ARTEMISS [
23]) respectively, and quantitatively analyzed the relationship between the LST and LSE inversion errors and each of the impact factors. After this Introduction section,
Section 2 will introduce the data simulation process, the TES algorithms and the evaluation metrics for the inversion results, and
Section 3,
Section 4 and
Section 5 will present the results, discussion and conclusions, respectively.
4. Discussion
The research in this paper focused on the influence of the core performance indicators, in terms of spectra and noise, of the common HTIR instruments on the LST and LSE inversion errors. It provided valuable references for the performance improvement and data processing of HTIR instruments. However, further research is needed on the relationship between the inversion errors and the performance of the HTIR instruments with higher spectral resolution and the validation of the modelled effects with real thermal infrared hyperspectral data. Moreover, during the result analysis, we found that ISSTES method has its own limitations. It is intelligible that after all, every TES method has certain assumptions. This also results in some TES methods (for example, ISSTES) being strong sensitive to data distortion while some TES methods (for example, ARTEMISS) being weakly sensitive to data distortion. Therefore, while committed to improving the performance of the instrument, we should also pay attention to how to further improve the anti-data-distortion ability of the TES method, especially noise. Instrument noise cannot be avoided or precisely corrected with the data itself, while center wavelength shift and FWHM change can usually be corrected by some methods to reduce the errors. Besides, atmospheric correction error is also affected by factors such as center wavelength shift, FWHM change and NEDT, which in turn affects the error of the TES algorithms. It needs further research, especially, whether the superposition effect will amplify the influence of instrument spectral parameter changes and noise on the TES algorithms.