Analysis of Space-Based Observed Infrared Characteristics of Aircraft in the Air

: The space-based infrared observatory of aircraft in the air has the advantages of wide-area, full-time, and passive detection. The optical design parameters for space-based infrared sensors strongly rely on target observed radiation, but there is still a lack of insight into the causes of aircraft observation properties and the impact of instrument performance. A simulation model of space-based observed aircraft infrared characteristics was constructed for this provision, coupling the aircraft radiance with background radiance and instrument performance effects. It was validated by comparing the model predictions to data from both space-based and ground-based measurements. The validation results reveal the alignment between measurements and model predictions and the dependence of overall model accuracy on the background. Based on simulations, the radiance contributions of aircraft and background are quantitatively evaluated, and the detection spectral window for ﬂying aircraft and its causes are discussed in association with instrumental performance effects. The analysis results indicate that the target-background (T-B) contrast is higher in the spectral ranges where aircraft radiation makes an important contribution. The background radiance plays a signiﬁcant role overall, while the observed radiance at 2.5–3 µ m is mainly from skin reﬂection and plume radiance. The skin-reﬂected radiation absence affects the model reliability, and its reduction at nighttime reduces the T-B contrast. The difference in T-B self-radiation and the stronger atmospheric attenuation for background contribute to the higher contrast at 2.7 µ m compared to the other spectral bands.


Introduction
The airplane's invention revolutionized how humans travel, connecting faraway places worldwide.The state of the aircraft during navigation has received extensive attention [1], which is the need for economic development and the focus of national defense construction.In today's complex aviation environment, no single technology can yet track all aircraft types in terms of global coverage [2].Space-based infrared imaging enables the acquisition of aircraft location information on a global scale without time constraints.There is currently no on-orbit infrared instrument designed for flying aircraft detection.The optical design parameters for space-based infrared sensors and the development of detection algorithms strongly rely on target observed radiation.However, there is still a lack of insight into the causes of aircraft observation properties and the effects of instrument performance.
The measurements can provide real infrared data of the aircraft under certain conditions, but the cost is very high.Accordingly, infrared target simulation modeling has been the focus of research in the military field for a long time [3][4][5][6][7].Battlefield requirements change with the development of technology.Mahulikar et al. [8][9][10] discussed the relationship between the infrared radiation level and the locking distance for the infrared modeling of military aircraft and proposed the concept of the infrared cross-section.Coiro [11][12][13] conducted infrared simulation modeling and sensitivity analysis of civil airplanes and unmanned combat air vehicles.The above methods are mainly used for air-to-air or ground-to-air detection scenarios and need to be refined for use in aircraft detectability assessments under space-based infrared observations.Yuan [14] proposed a multispectral integrated model for sea/cloud background radiation characteristics and analyzed the detection performance for the aircraft plume.Zhu [15] established an all-attitude motion characterization and parameter analysis system for aerial targets.These simulation models do not comprehensively account for the aircraft's observed radiation and lack preliminary validation work under the space-based infrared observation.Parts of the radiation were ignored, such as skin or background radiation.Accordingly, it is essential to establish and initially validate a space-based observed aircraft infrared characteristic model, coupling target radiation with background radiation and instrument performance effects.
In the application analysis of the simulation model, Mahulikar et al. [16][17][18][19][20] analyzed the influence of the internal/external radiation sources, the nozzle area, and the skin emissivity on the detection distance and carried out the optimization of the skin emissivity.The sources of infrared characteristic radiation from aircraft in the wide band (3-5 µm and 8-14 µm) are studied [21].The concept of modulating the infrared characteristics of an aircraft based on the skin anisotropic emissive characteristics was proposed and verified by simulation [22].Yuan [23] used the signal-to-noise ratio (SNR) and detection distance to analyze the detection capability of the geostationary infrared imaging system for aircraft plumes.Zhu and Yu et al. [24,25] used the comprehensive signal-to-noise ratio (CSNR) to evaluate the detection performance of infrared systems for the optimization of key parameters.In this literature, metrics such as SNR and CSNR were used to evaluate the detectability of infrared detection systems.Still, they cannot directly tell the influence of infrared radiation characteristics of the target, the background, and the system performance [26].However, there is still a lack of contribution evaluation of aircraft skin emission radiation, reflected radiation, plume radiation, and background radiation to space-based observation radiation.This is not conducive for an insight into the space-based infrared observational properties of the aircraft and the reasons for their formation.
The ground sampling distance (GSD), the modulation transfer function (MTF), the spectral response function (SRF), and the noise equivalent temperature difference (NE∆T) are important performance parameters of the infrared system and also widely known parameters in data applications.In the prior literature [23,24,27], the GSD, detection pixel size, and spectral band of the space-based infrared detection system were studied and analyzed.It was agreed that narrow bands outperformed wide bands [14], with spectral bands such as 2.65-2.9µm and 4.25-4.5 µm considered spectral detection windows for aerial targets.However, these studies focused more on the results or phenomena of optimization than the causes of the characteristic spectral bands and the effects of instrument performance on the relative difference between the target and background.
In view of the above, this paper builds a simulation model of space-based observed aircraft infrared characteristics and uses B7-12 of Gaofen-5 (GF-5) visual and infrared multispectral sensor (VIMS) [28] data and ground-measured plume data [6] for preliminary validation.A radiative contribution evaluation was carried out to discuss the contributions of background radiation, skin emission radiation, skin reflection radiation, and plume radiation at body-leaving and at-sensor radiance, as well as to analyze the effects of diurnal variation and spatial resolution on the radiative contribution.Then, the effect of instrument performance on the T-B contrast was discussed, and the formation of characteristic spectral bands was analyzed in conjunction with the radiance contribution evaluation and atmospheric attenuation effects.In Section 2, we described the simulation methodology, the flow, and the case of validation and analysis.In Section 3, simulated results were compared with space-based and ground-based data for validation.The contribution of aircraft and background radiation to the body-leaving and at-sensor radiance, and the effect of instrument performance on the T-B contrast, were then assessed.Lastly, the discussions and conclusions are provided in Sections 4 and 5, respectively.

Observed Radiance of Aircraft by Space-based Sensor
Within the instrument's linear response range, the measured radiance (also known as restored at-sensor radiance or restored onboard radiance for satellite remote sensing) can be assumed to be the result of a convolution process using an abstract concept of a sensor's equivalent response function and a random noise process, as shown in Formula (1) [29].
where L res_TOA and L TOA are the restored at-sensor radiance and the true top of atmosphere (TOA) radiance, respectively; R sensor is the sensor equivalent response, including the spatial imaging degradation and the spectral response, L noise is the effective instrument noise radiance.Due to the long observation distance, the solid angle of the aircraft is usually smaller than the instantaneous field of view, so the aircraft is a sub-pixel target in the space-based imaging system [14,27].Then, the aircraft radiation signal observed by the space-based sensor includes aircraft and background radiation.The TOA radiance, including the target and background signals, can be expressed in Formula (2).It is generally believed that the infrared radiation of aircraft mainly comes from skin and plume emission radiation [27].Besides, the aircraft skin reflected radiance is also taken into account in this paper.The plume is a block of high-temperature gas mass, and the true observed radiation is the coupling of the gas emission radiance and the background or nozzle radiance, as shown in Formula (3).
where S Skin , S Plume and S NP are the projected areas of the skin, the plume and the nozzle in the observation direction respectively, assuming that S Tar = S Skin + S Plume ; L E_skin and L R_skin are the emitted radiance and reflected radiance of the aircraft skin, respectively; L Plume , L E_plume , L Nozzle and L H Bkg are plume equivalent radiance, plume gas emission radiance, nozzle emission radiance, and background radiance at the flight altitude, respectively; τ Plume is the abstract transmittance of the plume gas.
The restored onboard radiance of the aircraft under the space-based infrared observation can be obtained by substituting Formulas (2) and (3) into Formula (1), as shown in Formula (4).The TOA radiance signal can be simplified into five parts: skin emission radiation, skin reflection radiation, plume radiation, background radiation, and atmospheric path radiation, as shown in Figure 1. is the abstract transmittance of the plume gas.The restored onboard radiance of the aircraft under the space-based infrared observation can be obtained by substituting Formulas (2) and (3) into Formula (1), as shown in Formula (4).The TOA radiance signal can be simplified into five parts: skin emission radiation, skin reflection radiation, plume radiation, background radiation, and atmospheric path radiation, as shown in Figure 1.

The Flow and Metrics of Analysis
Space-based observed aircraft infrared radiation is affected by aircraft radiation, background radiation, atmospheric effects, and instrument performance.To further investigate the causes of space-based observed aircraft infrared radiation, a quantitative assessment of the various impact factors should be carried out.As shown in Figure 2, the research framework in this paper is as follows: Firstly, a model of aircraft space-based infrared observation radiation was developed, coupling aircraft radiation, background radiation, atmospheric effects, and instrument performance characteristics.Preliminary validations were carried out using space-based data and plume static measurement data.Then, based on the simulation model, the contributions of the skin emitted/reflected radiation, plume radiation, background radiation, and path radiation were calculated and evaluated at the aircraft body-leaving radiance and at-sensor radiance.Finally, the effects

The Flow and Metrics of Analysis
Space-based observed aircraft infrared radiation is affected by aircraft radiation, background radiation, atmospheric effects, and instrument performance.To further investigate the causes of space-based observed aircraft infrared radiation, a quantitative assessment of the various impact factors should be carried out.As shown in Figure 2, the research framework in this paper is as follows: Firstly, a model of aircraft space-based infrared observation radiation was developed, coupling aircraft radiation, background radiation, atmospheric effects, and instrument performance characteristics.Preliminary validations were carried out using space-based data and plume static measurement data.Then, based on the simulation model, the contributions of the skin emitted/reflected radiation, plume radiation, background radiation, and path radiation were calculated and evaluated at the aircraft body-leaving radiance and at-sensor radiance.Finally, the effects of GSD, MTF, SRF, and noise on T-B contrast were analyzed regarding the current level of instrument performance.
of GSD, MTF, SRF, and noise on T-B contrast were analyzed regarding the current level of instrument performance.The simulation results were compared with the restored onboard radiance [29].Then the absolute error (AE) and relative error (RE) were adopted to evaluate the simulation accuracy, as shown in Formulas ( 5) and ( 6) [29].T-B contrast was adopted to assess the relative difference and relationship between the aircraft observed radiance and the background radiance.The contrast can be calculated by Formula (7) [30], and a positive number means that the aircraft is brighter than the background.The light-dark relationship in the simulation data was also used as a model reliability assessment metric.Another metric to assess the model's reliability is the consistency of the T-B light-dark relationship between measurements and predictions.The absolute mean values of these metrics were likewise calculated, as shown in Formula (8).The simulation results were compared with the restored onboard radiance [29].Then the absolute error (AE) and relative error (RE) were adopted to evaluate the simulation accuracy, as shown in Formulas ( 5) and ( 6) [29].T-B contrast was adopted to assess the relative difference and relationship between the aircraft observed radiance and the background radiance.The contrast can be calculated by Formula (7) [30], and a positive number means that the aircraft is brighter than the background.The light-dark relationship in the simulation data was also used as a model reliability assessment metric.Another metric to assess the model's reliability is the consistency of the T-B light-dark relationship between measurements and predictions.The absolute mean values of these metrics were likewise calculated, as shown in Formula (8).
where AE i , RE i , CR i , L SI MU ,i and L ONBOARD,i are the absolute error, relative error, T-B contrast, simulated radiance and restored onboard radiance of the i-band, respectively; L Bkg,i is the background radiance; L TB,i is the radiance of the pixel containing the target and background; MRE, MAE and MCR are the mean relative error, the mean absolute error and the mean contrast ratio, respectively; n is the number of bands.
To validate the accuracy of pure aircraft (without background) simulation and its influence on the reliability of the overall simulation model, an attempt was made to separate the influence of background radiance.Assuming that the measured signal of the aircraft is linearly mixed from the pure target signal and the background signal at the pupil, it can be expressed as Formula (9).
where ξ TB is the aircraft observed signal, including both aircraft and background signals; ξ Toa Bkg and ξ Toa Tar are the pure background and aircraft signals at the TOA, respectively; F is the aircraft signal factor, describing the contribution ratio of the target signal to the observed signal.The measured background radiance and the aircraft projected area ratio S Tar /d 2 can be used to estimate ξ Toa Bkg and F for the calculation of ξ Toa Tar , which can be compared with the simulation results of the aircraft without the background.
The spectral relative contribution (RC) of the aircraft skin emission radiance, reflected radiance, plume radiance, background radiance, and path radiance to the total body-leaving radiance and at-sensor radiance for the aircraft are calculated separately by Formula (10).In particular, the path radiance accounts for the flight altitude to the top of the atmosphere, which mainly affects the at-sensor radiance.The analysis of each component's relative contribution provides further insight into the role of each radiation source in the different spectral bands and the effect of atmospheric attenuation.
where L component presents radiance components, including the skin emission/reflection radiance, plume radiance, background radiance and path radiance (only for at-sensor radiance) observed at the body-leaving radiance and at-sensor radiance; L Tot is the total radiance, including the total body-leaving radiance and at-sensor radiance of aircraft.These coefficients can be derived by Formulas (3) and ( 4).The spectral T-B contrast under different instrument performances was calculated to evaluate each performance parameter's effect.Spectral bands with significant and prominent contrast in the target background were selected for further analysis.Finally, in conjunction with the analysis of instrument performance and radiance contributions, the causes of contrast in the characteristic spectral bands were discussed, and the advantages and disadvantages of each of these bands were compared.

Simulation Modeling 2.3.1. Skin Radiance
Aircraft skin radiance included skin emission and reflected radiance.The airframe is usually made of metal and coated, of which emitted radiance can be calculated using Planck's formula, as shown in Formula (11).
where ε is the emissivity of the skin; M BB is the blackbody irradiance; T is the temperature; h is the Planck constant; c is the speed of light; λ is the wavelength, and k is the Boltzmann constant.
During navigation, the aircraft skin temperature is mainly influenced by atmospheric aerodynamic heating, and the heating effect of solar radiation is smaller and can be neglected [17].The stagnation temperature is used to estimate the skin temperature, and the calculation formula can be expressed as in Formula (12) [30].
where T s is the skin temperature; T 0 is the temperature of the atmosphere around the aircraft; r is the recovery coefficient; γ is the specific heat ratio; Ma is the Mach number of flight.
In space-based observation, the skin-reflected radiance mainly considers the direct solar radiance and its scattered radiance, the cloud radiance, and the atmospheric thermal radiance.Defining the latter three as sky radiance, the skin body-leaving radiance can be expressed as: where ρ is the skin reflectivity, and ρ + ε = 1; L sd and L ↓ sky are direct solar radiance and sky downward radiance on the aircraft skin, respectively.Direct solar radiance and sky radiance can be estimated using direct solar irradiance E sd and atmosphere downward diffused irradiance E ↓ sky at a specified horizontal height, which can be derived from "flx" files of MODTRAN [31].Due to the long distance of space-based observation, the aircraft shape can be simplified to calculate the projected area of the aircraft skin [15].

Plume Radiance
Aircraft nozzle radiance was considered grey body radiation with an emissivity of about 0.9 [10], which can be calculated using Planck's formula.The plume, a non-uniformly distributed high-temperature gas, differs from the gray-body radiation characteristics of the surface.Its emission and absorption effects should be considered, and the radiative transfer equation can be expressed as Formula ( 14) [32].
where κ a is the absorption coefficient; L is the local radiance, L BB is the blackbody radiance, s and → s denote the position and optical path vector, respectively.The LOS method [33] was applied in this study to solve the radiative transfer equation.The non-uniform gas in a line of sight (LOS) is uniformly divided into multiple layers, as shown in Figure 3. Plume gas radiation of a LOS can be expressed as Formula (15).The nozzle or background radiance coupled with the plume gas can be calculated as Formulas ( 16) and (17), respectively.Thereby, the plume radiation intensity in Equation ( 3) can be shown as Formula (18).
where L i BB is the blackbody radiance of the ith slab; τ i denotes the transmissivity of the ith slab; (1 − τ i ) is the emissivity of the ith slab, and n represents the number of stratified layers; N is the total number of LOS intersecting the plume; M is the number of LOS intersecting both the nozzle and plume; ∆d is the spatial sampling interval of LOS, and the projected area of the plume in the observation direction can be expressed as where BB L is the blackbody radiance of the ith slab; i τ denotes the transmissivity of the ith slab; (1 ) is the emissivity of the ith slab, and n represents the number of stratified layers; N is the total number of LOS intersecting the plume; M is the number of LOS intersecting both the nozzle and plume; d Δ is the spatial sampling interval of LOS, and the projected area of the plume in the observation direction can be expressed as  The plume fluid field calculation aims to obtain the gas temperature, pressure distribution, and species content.The computational methods have been divided into two categories.One is the simplified model that uses empirical or semi-empirical formulations to obtain the plume fluid field [34,35].The other is the computational fluid dynamics that solves the Navier-Stokes equation to derive the plume fluid field [6,32,36].In contrast to the simplified model, the latter can obtain a fine plume fluid field; however, it requires tedious geometric model construction, meshing, and other manual involvement processes, which costs a large number of computational resources and time.Hence, a simplified model [35] is adopted to complete the calculation of the fluid field distribution.
The absorption coefficients of each species within the specified wavenumber η and temperature intervals can be calculated using the line-by-line method [37] with the aid of the high-temperature database HITEMP [38].

( ) ( , ) ( )
where ( , )  The plume fluid field calculation aims to obtain the gas temperature, pressure distribution, and species content.The computational methods have been divided into two categories.One is the simplified model that uses empirical or semi-empirical formulations to obtain the plume fluid field [34,35].The other is the computational fluid dynamics that solves the Navier-Stokes equation to derive the plume fluid field [6,32,36].In contrast to the simplified model, the latter can obtain a fine plume fluid field; however, it requires tedious geometric model construction, meshing, and other manual involvement processes, which costs a large number of computational resources and time.Hence, a simplified model [35] is adopted to complete the calculation of the fluid field distribution.
The absorption coefficients of each species within the specified wavenumber η and temperature intervals can be calculated using the line-by-line method [37] with the aid of the high-temperature database HITEMP [38].
where S i (η, T) is the line intensity of the ith spectral line at a given wavenumber when the temperature is T; F(η − η 0i ) is the line shape function of the ith spectral line, usually using the Voigt line function, and η 0i is the central wave number of the ith spectral line.The spectral line intensity S i (η, T) can be derived by extrapolation from the reference state line intensity S i (η, T re f ), as shown in Formula (20).
where T re f is the reference temperature; Q represents the partition function; E is the lower-state energy of the transition.

Background Radiation Calculation and Instrument Performance Simulation
The background at-sensor radiance can be expressed as Formula (21) [39], considering the atmospheric adjacency effect with the assumption of a flat subsurface.In the infrared band (>2.5 µm), the atmospheric adjacency effect contributes a relatively small amount of radiation to the at-sensor radiance [29].Thus, the scattered or reflected radiance from ground thermal radiation influenced by adjacency effects is neglected, as shown in Formula (22).The unknown quantities can be calculated by calling MODTRAN several times [40].
where A and A are coefficients describing solar radiance and atmospheric thermal radiance entering the sensor after reflection from the image-pixel surface, respectively; B and B are coefficients describing solar radiation and atmospheric thermal radiation reflected by the area-averaged ground surface into the sensor, respectively; C and C are coefficients describing solar radiation and atmospheric thermal radiation reaching the sensor after scattering in the atmosphere alone, respectively; S represents the atmospheric spherical albedo; ε t and ε b denotes the emissivity of the image-pixel surface and the area-averaged ground surface, respectively; τ d and τ s are the direct transmission along the line of sight and the effective transmission along the scattering path, respectively; T t and T b represent the image-pixel surface and area-averaged ground surface temperatures, respectively; In the case of cloud scenes, the physical parameters of the clouds are considered to be uniformly distributed horizontally within a single pixel, and the cloud radiance is acquired by setting the cloud parameters by "ICLD", "CTHIK", "CALT" and "CEXT" in the MODTRAN.Instrument performance simulations are achieved through spatial imaging degradation, spectral response and noise superimposition.Spatial degradation includes sampling interval degradation and imaging blur, τ h→toa atm L Tar S Tar + L h Bkg S Bkg /d 2 in (2) describes the process of spatial resampling, which can be done by setting the GSD to simulate instrument sampling interval degradation.Imaging blur is caused by factors such as the external imaging environment and the imaging capability of the instrument.The sub-pixel target energy appears distributed into several surrounding pixels [14,27].During instrumental imaging, the imaging blur can be seen as low-pass filtering of the full imaging chain on the ground scene, evaluated by the MTF or the PSF.Accordingly, the spatial domain imaging blur can be expressed as (23).A Gaussian function (24) is usually employed to fit the point spread function, and the MTF is the modulo of the PSF after the Fourier transform.PSF can be calculated by Formula (25).
The spectral response function is commonly applied to describe an instrument's spectral response characteristics [41].The effective spectral radiance obtained by the sensor is considered a weighted average of the continuous radiance spectrum and the spectral response function, as shown in Formula (26).Gaussian functions are often used to fit spectral response functions, where a spectral response curve can be solved for a given central wavelength and full width at half maximum (FWHM) by Formulas ( 27) and (28).
where L SRF_i , SRF i , CW L i , FW HM i are the spectral radiance, spectral response function and its central wavelength and full width at half maximum of the i-band, respectively.Noise superposition is achieved based on noise equivalent radiance [42], which is modeled by adding a normally distributed random number to each band radiance, as shown in Formula ( 29).The NE∆T is applied to gauge the noise level of the instrument.Each band's noise equivalent radiance (NER) can be calculated from the NE∆T given at the instrument's design, as shown in Formula (30).
where L n (i) denotes the i-band noise radiance; T B is the blackbody temperature that NE∆T is defined, and ∆T is the Planck function linear temperature difference.

Materials for Simulation Case Study
The space-based observed aircraft infrared characteristics model obtains the simulated aircraft observation radiance by inputting information such as ground reflectivity, temperature, and aircraft parameters.Two simulation experiments were designed to validate the model's fidelity; one using satellite-based data to validate the space-based observation simulation model for aircraft; and the other using ground-based measurements of the plume to validate the accuracy of the plume modeling, complementing the former.The simulation time and aircraft parameters from the validation of the space-based observation model were also used in the evaluation of the aircraft observed radiance contribution and the impact of instrument performance.
The onboard data gathered by the VIMS of GF-5 satellite were selected as a reference for the simulation validation, mainly to verify the simulation accuracy in the infrared segment.The VIMS provides images in 12 bands from the visible to thermal infrared, six of which were used for validation in this study.The central latitude and longitude of the chosen data are 32.212895• N and 126.392002 • E, located in the East China Sea and imaged on 25 June 2019, with the specific parameters shown in Table 1.
The aircraft's position in the data is shown in Figure 4a, and the aircraft is located above the cirrus, followed by the contrail.The Flightradar24 historical data query shows that the aircraft is a Boeing 777-246, flying from Tokyo, Japan, to Shanghai, China, with an altitude of 12,192 m and a flight speed of 218.64 m/s.The specific parameters are shown in Table 2.The aircraft's position in the data is shown in Figure 4a, and the aircraft is located above the cirrus, followed by the contrail.The Flightradar24 historical data query shows that the aircraft is a Boeing 777-246, flying from Tokyo, Japan, to Shanghai, China, with an altitude of 12,192 m and a flight speed of 218.64 m/s.The specific parameters are shown in Table 2.As shown in Figure 4c, it is not guaranteed that the aircraft is in a given pixel, and its position is inconsistent between different bands due to inter-band offsets.Accordingly, the target signal needs to be extracted band by band.To ensure that the aircraft signals are within a single pixel for simulation validation, the 3 × 3 and 4 × 4 pixel sizes were merged (shown in the black box in Figure 4c) to obtain the mixed signals at 120 m and 160 m spatial resolution.
The sea surface temperature was obtained from the SST CCI data of ECMWF [43] as 295 K.The sea surface reflectance was used from the ECOSTRESS spectral library [44], and the skin emissivity was set to 0.6 [45].The aircraft skin is considered a Lambertian body with emissivity and reflectivity summed into one.The cloud thickness at the aircraft location was inverted using MODTRAN and the measured cloud data (average of the green areas in Figure 4a).Under the cirrus assumption, the cloud base and top altitudes are 8.1 and 9.2 km, at which point the simulation results are closest to the measured results.
The main purpose of the simulation validation experiment was to examine the ability to simulate the space-based observational characteristics of the pixel containing the aircraft, and the spectral response was considered to describe the VIMS instrument performance.The spectral response functions of the B7-B8 were generated based on the spectral ranges given in Table 1 with Gaussian functions, and the B9-B12 spectral response functions were provided by the Numerical Weather Prediction Satellite Application Facility [46], as shown in Figure 4b.
As the VIMS cannot capture the infrared spectral characteristics of the aircraft plume, such as 4.2 µm, the Swedish Defense Research Agency's engine plume measurements [6] were applied to validate the aircraft plume model.The plume simulation considered the effect of CO 2 and H 2 O, with the gas velocity of Mach 0.6, the ambient atmospheric temperature of 290 K, the atmospheric pressure of 101 kPa, air humidity of 35% and detection distance of 20 m perpendicular to the plume.Horizontal path atmospheric attenuation and path radiation were also considered.

Space-Based Simulation Validation
The simulation results of the aircraft-observed radiation for space-based infrared observations were validated by the B7-12 data from VIMS, as shown in Figure 5.The background (green box in Figure 4a) spectral mean and its distribution were calculated and compared with the simulation results, as shown in Figure 5a.The 3 × 3 and 4 × 4 pixel size resampling were selected to extract the aircraft observed radiance spectra (containing both aircraft and background), which were compared with the simulation results as in Figure 5b.The comparison results of aircraft and background spectra curves are given in Figure 5c, d.The results show that the measured radiance and T-B relationship agree with the established model.
In order to quantify the errors between simulations and measurements, the relative and absolute errors at different spatial scales and the T-B contrast were calculated, respectively, as shown in Table 3.The MRE of the aircraft observation characteristics simulations for 3 × 3 and 4 × 4 pixel sizes are 8.32% and 6.42%.The maximum contribution of RE is the B7 band (−28.46%,−20.40%), which has low radiance with sensitivity to errors, as corroborated by the AE.Compared with the simulation accuracy of aircraft observation characteristics, the RE of pure aircraft simulation is larger.The MREs of pure target simulation are 71.22% and 56.71%, and the maximum contribution of RE is B7 (161.95%,123.92%).Besides, the simulated and measured cloud backgrounds were also compared as shown in Figure 5a, with an MRE of 4.52%.

Space-Based Simulation Validation
The simulation results of the aircraft-observed radiation for space-based infrared observations were validated by the B7-12 data from VIMS, as shown in Figure 5.The background (green box in Figure 4a) spectral mean and its distribution were calculated and compared with the simulation results, as shown in Figure 5a.The 3 × 3 and 4 × 4 pixel size resampling were selected to extract the aircraft observed radiance spectra (containing both aircraft and background), which were compared with the simulation results as in Figure 5b.The comparison results of aircraft and background spectra curves are given in Figure 5c, d In order to quantify the errors between simulations and measurements, the relative and absolute errors at different spatial scales and the T-B contrast were calculated, respectively, as shown in Table 3.The MRE of the aircraft observation characteristics simulations for 3 × 3 and 4 × 4 pixel sizes are 8.32% and 6.42%.The maximum contribution of RE is the B7 band (−28.46%,−20.40%), which has low radiance with sensitivity to errors, as corroborated by the AE.Compared with the simulation accuracy of aircraft observation characteristics, the RE of pure aircraft simulation is larger.The MREs of pure target simulation are 71.22% and 56.71%, and the maximum contribution of RE is B7 (161.95%,123.92%).This phenomenon indicates that the accuracy of the overall model is more dependent on the background simulation because the aircraft contribution is smaller.Therefore, the MRE decreases and converges to the background simulation accuracy (MRE of 4.52%), along with the spatial resolution decreases.Meanwhile, it can be found that the simulation accuracy of pure aircraft also changes with the scale (it should not change theoretically), which shows a deviation in the estimation of the aircraft signal factor.There may be two reasons for this deviation.One is that there are unknowns or deviations in aircraft parameters, such as the observation angle of the aircraft, the actual size, etc., resulting in the inability to estimate the aircraft signal factor effectively; the other is that the aircraft projected area ratio may introduce errors.
Objectively, most of the research on infrared signature analysis is controlled by military research institutions with limited details in the open literature [47].The lack of target-related data, especially measured space-based infrared data for the aircraft, is an important factor restricting the validation and improvement of space-based infrared imaging simulation models.Accordingly, the consistency of the T-B relationship between measurements and predictions can also be used to evaluate the reliability of the simulation model.The simulation and the measured T-B contrast results show that the aircraft observed radiance in the B7 band is higher than the background (brighter than the background), while the opposite is true in other bands.This phenomenon indicates alignment between the simulated and measured T-B relative relationships.

Plume Simulation Validation
Plume infrared characteristic measurement experiments were carried out on the engine test stand [6].The results were used to validate the plume simulation accuracy and to supplement the space-based observation simulation validation.The comparison between the measured and simulated results is shown in Figure 6.The MRE in the 4.1-5 µm was calculated to be approximately 61.64% (excluding the position of strong atmospheric absorption).The results show a good agreement between simulation and measurement, with the same spectral characteristics.The partially unknown parameters of the experimental environment have caused errors between the simulation and the experiment.In the reference [4], the spectral radiation intensity is accurate up to 50% after considering all uncertainties associated with the input.By contrast, the simulation accuracy in the paper has been able to achieve a relatively good result with the condition of unknown input parameter uncertainty.

Evaluation of the Aircraft and Background Contribution to the Observed Radiance
It is widely recognized that background radiation affects the characteristics of aircraft space-based observations, especially at low spatial resolutions where the target features are not obvious mixed with the background radiation.However, it still lacks a quantitative analysis of background and aircraft contributions to observed radiation gathered by space-based infrared sensors.The contributions of each radiance component at 2.5-13 μm were calculated for the body-leaving radiance and TOA radiance at different spatial reso- The results show a good agreement between simulation and measurement, with the same spectral characteristics.The partially unknown parameters of the experimental environment have caused errors between the simulation and the experiment.In the reference [4], the spectral radiation intensity is accurate up to 50% after considering all uncertainties associated with the input.By contrast, the simulation accuracy in the paper has been able to achieve a relatively good result with the condition of unknown input parameter uncertainty.

Evaluation of the Aircraft and Background Contribution to the Observed Radiance
It is widely recognized that background radiation affects the characteristics of aircraft space-based observations, especially at low spatial resolutions where the target features are not obvious mixed with the background radiation.However, it still lacks a quantitative analysis of background and aircraft contributions to observed radiation gathered by spacebased infrared sensors.The contributions of each radiance component at 2.5-13 µm were calculated for the body-leaving radiance and TOA radiance at different spatial resolutions and day/night conditions, as shown in Figures 7 and 8.The night scenes mainly considered the absence of solar radiation.
The results show a good agreement between simulation and measurement, with the same spectral characteristics.The partially unknown parameters of the experimental environment have caused errors between the simulation and the experiment.In the reference [4], the spectral radiation intensity is accurate up to 50% after considering all uncertainties associated with the input.By contrast, the simulation accuracy in the paper has been able to achieve a relatively good result with the condition of unknown input parameter uncertainty.

Evaluation of the Aircraft and Background Contribution to the Observed Radiance
It is widely recognized that background radiation affects the characteristics of aircraft space-based observations, especially at low spatial resolutions where the target features are not obvious mixed with the background radiation.However, it still lacks a quantitative analysis of background and aircraft contributions to observed radiation gathered by space-based infrared sensors.The contributions of each radiance component at 2.5-13 μm were calculated for the body-leaving radiance and TOA radiance at different spatial resolutions and day/night conditions, as shown in Figures 7 and 8  As shown in Figure 7a,c, the contribution of skin-reflected radiance to the body-leaving radiance at 2.5-3 µm is up to 98%, while the plume radiance occupies a smaller proportion but is still higher than the background radiance.Moreover, the plume and skin-reflected radiance also dominate near 4.3 µm, as the atmosphere has a strong absorption effect at these two spectral ranges, resulting in the lower background radiation energy.At all other spectrum bands, the background radiation makes up a large radiative contribution, especially in the long-wave infrared band, where it accounts for more than 90% and up to 99%.As shown in Figure 7d-f, the contribution of background radiation shows the same trend between daytime and nighttime.Still, the contribution of skin-reflected radiance at nighttime could be negligible due to the absence of solar radiation.
As seen in Figure 8, the TOA radiance increases the atmosphere path radiance from the aircraft altitude to the sensor.Atmosphere path radiation is an important component around 4.3, 6, and 9.5 µm, especially at 4.3 µm, where it accounts for up to 100%.Figures 7a-c and 8a-c show that the contribution of plume radiation is reduced compared to that of body-leaving radiance.H 2 O and CO 2 are both the main sources of thermal radiation in the plume and the main species of atmospheric attenuation.Therefore, the plume radiation energy is easily attenuated by the atmosphere, which makes it hard to gather signals of the plume from space-based sensors.The skin-reflected radiation is the largest variable in the TOA radiance between daytime and nighttime, which is reduced by the absence of solar radiation, similar to the diurnal variation in the body-leaving radiance.As shown in Figure 7a,c, the contribution of skin-reflected radiance to the body-leaving radiance at 2.5-3 μm is up to 98%, while the plume radiance occupies a smaller proportion but is still higher than the background radiance.Moreover, the plume and skinreflected radiance also dominate near 4.3 μm, as the atmosphere has a strong absorption effect at these two spectral ranges, resulting in the lower background radiation energy.At all other spectrum bands, the background radiation makes up a large radiative contribution, especially in the long-wave infrared band, where it accounts for more than 90% and up to 99%.As shown in Figure 7d-f, the contribution of background radiation shows the same trend between daytime and nighttime.Still, the contribution of skin-reflected radiance at nighttime could be negligible due to the absence of solar radiation.
As seen in Figure 8, the TOA radiance increases the atmosphere path radiance from the aircraft altitude to the sensor.Atmosphere path radiation is an important component around 4.3, 6, and 9.5 μm, especially at 4.3 μm, where it accounts for up to 100%.Figure 7a-c and Figure 8a-c show that the contribution of plume radiation is reduced compared to that of body-leaving radiance.H2O and CO2 are both the main sources of thermal radiation in the plume and the main species of atmospheric attenuation.Therefore, the plume radiation energy is easily attenuated by the atmosphere, which makes it hard to gather signals of the plume from space-based sensors.The skin-reflected radiation is the largest variable in the TOA radiance between daytime and nighttime, which is reduced by the absence of solar radiation, similar to the diurnal variation in the body-leaving radiance.
Jointly, the blue areas of Figures 7 and 8 illustrate that the background radiation has an extremely high radiative contribution overall and increases with GSD.The skin emission radiance is mainly concentrated at 5-13 μm, and its relative contribution decreases  (d-f) are the contributions plots in the nighttime; the blue, red, yellow, purple, and green areas represent the relative contribution of the background radiance, plume radiance, skin reflected radiance, skin emission radiance, and atmosphere path radiance respectively.
Jointly, the blue areas of Figures 7 and 8 illustrate that the background radiation has an extremely high radiative contribution overall and increases with GSD.The skin emission radiance is mainly concentrated at 5-13 µm, and its relative contribution decreases with decreasing spatial resolution, from 5% to 1%.The yellow areas in Figures 7 and 8 together illustrate the importance of skin-reflected radiance, which accounts for a high proportion of both the body-leaving and TOA radiance in the daytime.As shown in Figure 9, it is inconsistent with the real relative relationship (Table 3) that the aircraft observed radiance without the skin reflected radiation in the B7 band is lower than the background radiance.This phenomenon highlights the importance and necessity of considering skin-reflected radiance in the simulation model.The diurnal variation of the skin-reflected radiance is also noteworthy. Figure 10 indicates a significant contrast difference between daytime and nighttime, particularly around 2.5-4.15µm, where the contrast is most significantly reduced and even negative at nighttime.Meanwhile, there is a consistently high level of contrast at 2.7 µm compared with other spectral bands, despite the reduced contrast in the nighttime.
Remote Sens. 2023, 15, x FOR PEER REVIEW 17 of 25 with decreasing spatial resolution, from 5% to 1%.The yellow areas in Figures 7 and 8 together illustrate the importance of skin-reflected radiance, which accounts for a high proportion of both the body-leaving and TOA radiance in the daytime.As shown in Figure 9, it is inconsistent with the real relative relationship (Table 3) that the aircraft observed radiance without the skin reflected radiation in the B7 band is lower than the background radiance.This phenomenon highlights the importance and necessity of considering skin-reflected radiance in the simulation model.The diurnal variation of the skin-reflected radiance is also noteworthy. Figure 10 indicates a significant contrast difference between daytime and nighttime, particularly around 2.5-4.15μm, where the contrast is most significantly reduced and even negative at nighttime.Meanwhile, there is a consistently high level of contrast at 2.7 μm compared with other spectral bands, despite the reduced contrast in the nighttime.

Analysis of the Effect of Instrument Performance on Target-Background Contrast
The instrument performance is an important factor affecting the aircraft observ characteristics.For space-based observations with long observation distances, unid fied objects, and unpredictable instrument performance in orbit, it is required to exa the effect of instrument performance parameters on the relative differences betwee target and background.The impacts of instrument performance parameters such as MTF, SRF, and NEΔT on the T-B contrast were analyzed according to imaging tim mosphere, and aircraft parameters in the space-based simulation validation session.

Analysis of the Effect of Instrument Performance on Target-Background Contrast
The instrument performance is an important factor affecting the aircraft observation characteristics.For space-based observations with long observation distances, unidentified objects, and unpredictable instrument performance in orbit, it is required to examine the effect of instrument performance parameters on the relative differences between the target and background.The impacts of instrument performance parameters such as GSD, MTF, SRF, and NE∆T on the T-B contrast were analyzed according to imaging time, atmosphere, and aircraft parameters in the space-based simulation validation session.
With regard to aircraft space-based observation characteristics, spatial degradation is the most intuitive impact factor.Figure 11 shows the T-B contrast spectra at different GSDs of 70-400 m.The spectral contrasts were shown in three segments due to the large difference between the different spectral ranges.The contrast becomes smaller in absolute terms as GSD increases, with a 97% reduction from 70 to 400 m.Because the lower the spatial resolution is, the higher the background radiance contribution is in the aircraft pixel, which makes the aircraft observed radiance closer to the background radiance.The contrast curves for different MTFs at the GSD of 120 m were presented in Figure 12, demonstrating that T-B relative differences increase with increasing MTF.The results in Figures 11 and 12 show consistently high contrast around 2.7 μm and generally higher contrast at 2.5-3.5 μm than that at 3.5-13 μm.Therefore, 2.7 μm was used in this paper as a general reference for 2.5-3.5 μm.Within 3.5-13 μm, contrast peaks occur around 4.2 μm, 4.4 μm and 5.7 μm.According to Figure 8a-c, it can be found that the skinreflected radiation and the plume radiation at 2.7, 4.2, and 4.4 μm play a dominant role in making the contrast higher.In comparison, the contrast of 5.7 μm is smaller due to the small contribution of aircraft radiation.Besides, it should be noted that 72% of contrasts are negative at 2.5-13 μm, indicating that the aircraft pixel is darker than the pure background.
Apart from spatial degradation, the instrument's spectral response is a non-negligible influence.The T-B contrast and the SRF for VIMS B7-12 at the GSD of 120 m are given in Figure 13.The B7 of VIMS has a central wavelength of 3.68 μm and an FWHM of 0.35 μm.The SRF of B7 does not cover the high contrast region of 2.7 μm, so the T-B contrast The results in Figures 11 and 12 show consistently high contrast around 2.7 µm and generally higher contrast at 2.5-3.5 µm than that at 3.5-13 µm.Therefore, 2.7 µm was used in this paper as a general reference for 2.5-3.5 µm.Within 3.5-13 µm, contrast peaks occur around 4.2 µm, 4.4 µm and 5.7 µm.According to Figure 8a-c, it can be found that the skin-reflected radiation and the plume radiation at 2.7, 4.2, and 4.4 µm play a dominant role in making the contrast higher.In comparison, the contrast of 5.7 µm is smaller due to the small contribution of aircraft radiation.Besides, it should be noted that 72% of contrasts are negative at 2.5-13 µm, indicating that the aircraft pixel is darker than the pure background.
Apart from spatial degradation, the instrument's spectral response is a non-negligible influence.The T-B contrast and the SRF for VIMS B7-12 at the GSD of 120 m are given in Figure 13.The B7 of VIMS has a central wavelength of 3.68 µm and an FWHM of 0.35 µm.The SRF of B7 does not cover the high contrast region of 2.7 µm, so the T-B contrast observed is not high.Meanwhile, the contrast in the spectral response region of the B8-12 is negative, so the aircraft pixel in these bands is darker than the background.observed is not high.Meanwhile, the contrast in the spectral response region of the B8-12 is negative, so the aircraft pixel in these bands is darker than the background.The above discussion illustrates the importance of designing/selecting the appropriate spectral response region for satellite instruments aiming at aerial target detection.The effect of SRF was further analyzed by calculating the T-B contrast for each spectral band at central wavelengths of 2.7, 4.2, 4.4, and 5.7 μm and FWHMs of 0.02, 0.04, 0.08, and 0.16 μm.To account for the impact of spectral calibration error, a random variable with a Gaussian distribution and a standard deviation (STD) of 1 nm was added to the center wavelength.T-B contrast at different central wavelengths and its variation ranges were calculated by several iterations, as shown in Figure 14.The above discussion illustrates the importance of designing/selecting the appropriate spectral response region for satellite instruments aiming at aerial target detection.The effect of SRF was further analyzed by calculating the T-B contrast for each spectral band at central wavelengths of 2.7, 4.2, 4.4, and 5.7 µm and FWHMs of 0.02, 0.04, 0.08, and 0.16 µm.
To account for the impact of spectral calibration error, a random variable with a Gaussian distribution and a standard deviation (STD) of 1 nm was added to the center wavelength.T-B contrast at different central wavelengths and its variation ranges were calculated by several iterations, as shown in Figure 14.
effect of SRF was further analyzed by calculating the T-B contrast for each spectral band at central wavelengths of 2.7, 4.2, 4.4, and 5.7 μm and FWHMs of 0.02, 0.04, 0.08, and 0.16 μm.To account for the impact of spectral calibration error, a random variable with a Gaussian distribution and a standard deviation (STD) of 1 nm was added to the center wavelength.T-B contrast at different central wavelengths and its variation ranges were calculated by several iterations, as shown in Figure 14.The contrast is constantly the highest at the central wavelength of 2.7 μm, as shown in Figure 14.The contrast varies inconsistently with the FWHM, indicating that the optimal FWHM for each spectral band is relatively independent.Furthermore, contrast variations caused by central wavelength shifts should be of attention.According to the error bar length, the 2.7 μm contrast variation range is greater and more sensitive to spectral calibration errors.The error bars of the 2.7 μm are equidistant in length between the positive and negative axes, suggesting that the central wavelength could be further optimized to improve the contrast, with the 5.7 μm position showing a similar phenomenon.It is also found that the length of the error bars decreases with the increase of FWHM, indicating that the wider the spectral response range is, the less the contrast is affected by the central wavelength shift.However, an excessively wide spectral response range also The contrast is constantly the highest at the central wavelength of 2.7 µm, as shown in Figure 14.The contrast varies inconsistently with the FWHM, indicating that the optimal FWHM for each spectral band is relatively independent.Furthermore, contrast variations caused by central wavelength shifts should be of attention.According to the error bar length, the 2.7 µm contrast variation range is greater and more sensitive to spectral calibration errors.The error bars of the 2.7 µm are equidistant in length between the positive and negative axes, suggesting that the central wavelength could be further optimized to improve the contrast, with the 5.7 µm position showing a similar phenomenon.It is also found that the length of the error bars decreases with the increase of FWHM, indicating that the wider the spectral response range is, the less the contrast is affected by the central wavelength shift.However, an excessively wide spectral response range also reduces the contrast.It is thus essential that the appropriate spectral response is selected while properly controlling the spectral calibration error.
Infrared radiation is less energetic, and remote sensing data is more susceptible to instrument noise compared with the visible.It is often expected that the radiance difference between the target and background exceeds the NER, ensuring that the noise does not obscure the target signal.Therefore, this paper compared the T-B radiance difference with the NER of VIMS and ASTER instruments [29,48] and analyzed the effect of noise on contrast.Figure 15 shows the T-B radiance difference at the GSD of 120 m and the corresponding NER at the NE∆T of 0.15K@300K and 0.3K@300K.The results show that the T-B radiance difference at 4.25 µm, 4.57 µm, and 5-7 µm is less than NER with more susceptibility to noise.reduces the contrast.It is thus essential that the appropriate spectral response is selected while properly controlling the spectral calibration error.Infrared radiation is less energetic, and remote sensing data is more susceptible to instrument noise compared with the visible.It is often expected that the radiance difference between the target and background exceeds the NER, ensuring that the noise does not obscure the target signal.Therefore, this paper compared the T-B radiance difference with the NER of VIMS and ASTER instruments [29,48] and analyzed the effect of noise on contrast.Figure 15 shows the T-B radiance difference at the GSD of 120 m and the corresponding NER at the NEΔT of 0.15K@300K and 0.3K@300K.The results show that the T-B radiance difference at 4.25 μm, 4.57 μm, and 5-7 μm is less than NER with more susceptibility to noise.The noise was added in the aircraft observed radiation, assuming that the background radiance is the result of averaging over a uniform scene and is not affected by noise.Figure 16 shows the range of variation in contrast and its standard deviation (STD).The results indicate that the standard deviation is greatest, around 2.7 μm, but its effect is The noise was added in the aircraft observed radiation, assuming that the background radiance is the result of averaging over a uniform scene and is not affected by noise.Figure 16 shows the range of variation in contrast and its standard deviation (STD).The results indicate that the standard deviation is greatest, around 2.7 µm, but its effect is almost negligible due to the large contrast.The range of contrast variation around 4.  The noise was added in the aircraft observed radiation, assuming that the b ground radiance is the result of averaging over a uniform scene and is not affecte noise.Figure 16 shows the range of variation in contrast and its standard deviation (S The results indicate that the standard deviation is greatest, around 2.7 μm, but its eff almost negligible due to the large contrast.The range of contrast variation around μm and 6μm covers the zero axis, indicating a change in the T-B relative relation (light-dark relationships).It is obvious that a change in the T-B light-dark relationsh not expected, which seriously affects the distinction between the target and backgrou

Discussions
The validation and evaluation results illustrate that the proposed model can generate accurate simulation data consistent with the measured data.However, the spectral window and challenges for aircraft detection are still open for further discussion.

Space-Based Infrared Detection Spectral Window for Aircraft
In this paper, aircraft and background radiation contributions are analyzed quantitatively, with a focus on the 2.7 µm, 4.2 µm, 4.4 µm, and 5.7 µm bands as influenced by instrument performance.As seen from Section 3.2, aircraft skin reflections and plume radiance play an important radiance contribution in these spectral bands.This phenomenon indicates that the target-background contrast is higher in the spectral ranges where aircraft radiation makes an important contribution.Because of this, the contrast is consistently higher at 2.5-3.5 µm, where aircraft radiation can contribute up to 98%.
The stronger atmospheric attenuation for the background, in addition to the difference in radiation between the aircraft and the background itself, is an important cause of the higher contrast at 2.7 µm.The atmospheric transmittance from a 0 to 12 km altitude to the TOA and atmospheric profile of CO 2 and H 2 O were given in Figure 17.The whole atmosphere transmittance is almost zero around 2.7, 4.3, and 5-8 µm, which limits the Earth's surface radiation (only path radiance remained) observed by satellite remote sensing.It is the lower atmospheric transmittance and the path radiation at 2.7 µm that results in lower background radiation, meaning that only very small aircraft radiation is required to create high contrast.As seen from Figure 17b, the main gas molecules of atmospheric attenuation are distributed at a 0-10 km altitude, which objectively enables some nonatmospheric windows to be detection spectral windows for aerial targets in the air.
atmosphere transmittance is almost zero around 2.7, 4.3, and 5-8 μm, which limits the Earth's surface radiation (only path radiance remained) observed by satellite remote sensing.It is the lower atmospheric transmittance and the path radiation at 2.7 μm that results in lower background radiation, meaning that only very small aircraft radiation is required to create high contrast.As seen from Figure 17b, the main gas molecules of atmospheric attenuation are distributed at a 0-10 km altitude, which objectively enables some nonatmospheric windows to be detection spectral windows for aerial targets in the air.The degree of T-B contrast affected by instrument performance varies across the characteristic spectral bands.Undeniably, the T-B contrast around 2.7 μm remains consistently high compared to other spectral bands.From Figures 14 and 16, the contrast at 4.2 μm is more sensitive to FWHM and spectral calibration accuracy, while the contrast at 4.4 μm is more susceptible to instrument noise.Moreover, 5.7 μm does not have a significant advantage compared to other spectral bands.

The Challenge of Space-Based Infrared Detection of Aircraft
In previous research, aircraft detection could make use of spatial [49], spectral [50], and motion information [2,51] from the aircraft.However, the low energy and long range of infrared remote sensing limit spatial and spectral information applications in detection.Some researchers [23,24,27] used SNR and CSNR to select the feature spectrum bands, which have the potential to detect and track aircraft based on motion characteristics [52,53].This paper uses radiance contribution analysis to explain why these spectral bands were selected.It then remains doubtful whether aircraft identification can be achieved The degree of T-B contrast affected by instrument performance varies across the characteristic spectral bands.Undeniably, the T-B contrast around 2.7 µm remains consistently high compared to other spectral bands.From Figures 14 and 16, the contrast at 4.2 µm is more sensitive to FWHM and spectral calibration accuracy, while the contrast at 4.4 µm is more susceptible to instrument noise.Moreover, 5.7 µm does not have a significant advantage compared to other spectral bands.

The Challenge of Space-Based Infrared Detection of Aircraft
In previous research, aircraft detection could make use of spatial [49], spectral [50], and motion information [2,51] from the aircraft.However, the low energy and long range of infrared remote sensing limit spatial and spectral information applications in detection.Some researchers [23,24,27] used SNR and CSNR to select the feature spectrum bands, which have the potential to detect and track aircraft based on motion characteristics [52,53].This paper uses radiance contribution analysis to explain why these spectral bands were selected.It then remains doubtful whether aircraft identification can be achieved using a single band.Therefore, the method of aircraft identification for space-based infrared observation is still a challenge.The primary question is what information about the aircraft is used to achieve detection or identification, which affects instrument design and algorithm development.
Both traditional and artificial intelligence algorithms require a certain amount of measured data to achieve feature analysis and training.However, the lack of space-based infrared measurements has limited the study of aircraft detection algorithms.The publicly available infrared datasets are mostly derived from ground-based measurements, and the target size in the images is not representative of space-based observations.It is costly to collect a large-scale dataset with accurate pixel-level annotations [54].The simulation model presented in this paper has the potential to provide space-based infrared data containing aircraft for feature analysis and network training.

Conclusions
In this paper, a simulation model of space-based observed aircraft infrared characteristics is established, coupling target, background radiation, and instrument performance effects.The accuracy of the simulation model was validated by comparing the model predictions with data from space-based and ground-based measurements.Validation results reveal that the measured radiance and T-B relationship agree with the established model.It is also found that the model of space-based observed aircraft infrared characteristics is more dependent on background simulation accuracy.To further understand the causes of the aircraft observed characteristics, the contributions of background radiation, skin reflected/emission radiation, atmosphere path, and plume radiation were evaluated.The evaluation of radiance components indicates that background radiance plays a major role overall, while the observed radiance at 2.5-3 µm is mainly from skin reflection and plume radiance.The lack of skin-reflected radiance decreases the model reliability, and its reduc-tion at nighttime reduces the T-B contrast.The effect of instrument performance parameters on space-based infrared detection was analyzed, and the different changes in contrast on detection windows of 2.7, 4.2, 4.4, and 5.7 µm were highlighted.The results show a consistently high level of contrast at 2.7 µm compared with other spectral bands, although it is susceptible to diurnal variations.The discussions denote that the target-background contrast is higher in the spectral ranges where aircraft radiation makes an important contribution.The difference in T-B self-radiation and the stronger atmospheric attenuation for background contribute to the higher contrast at 2.7 µm.
The model proposed in this paper can provide data for space-based infrared detection algorithm development and onboard instrument performance evaluation.The analysis and discussion based on this model provide insight into the causes of target observation characteristics and the effect of instrument performance on the T-B relative difference.
and path radiance from the aircraft altitude to the TOA, respectively; L Tar and L H Bkg are aircraft and background equivalent radiance at flight altitude H, respectively; S Tar and S Bkg are the projected area of the aircraft and background in the viewing direction, respectively; d is the ground sampling distance, assuming that d 2 = S Tar + S Bkg .
of the skin emission, reflected radiance, plume radiance and background radiance reaching the sensor, respectively.

Figure 1 .
Figure 1.Diagrammatic sketch of aircraft infrared observation using space-based sensors.Schematic CO2 and H2O column density diagrams with altitude are in the upper left corner.

Figure 1 .
Figure 1.Diagrammatic sketch of aircraft infrared observation using space-based sensors.Schematic CO 2 and H 2 O column density diagrams with altitude are in the upper left corner.

Figure 2 .
Figure 2. The flow chart of evaluation and analysis in this study.

Figure 2 .
Figure 2. The flow chart of evaluation and analysis in this study.

Figure 3 .
Figure 3. Observation line of sight uniform division schematic.lΔ is the length of a slab; N is the total number of sight lines intersecting with the plume; μ and ν are the elevation and azimuth angles relative to the aircraft; P, T and X represent the gas pressure, temperature and species content respectively; n represents the number of stratified layers; d Δ is the spatial sampling interval of LOS; n BB L is the blackbody radiance of the n-th slab; n τ denotes the transmissivity of the n-th slab.

−
intensity of the ith spectral line at a given wavenumber when the temperature is T ; is the line shape function of the ith spectral line, usually using the Voigt line function, and 0i η is the central wave number of the ith spectral line.

Figure 3 .
Figure 3. Observation line of sight uniform division schematic.∆l is the length of a slab; N is the total number of sight lines intersecting with the plume; µ and ν are the elevation and azimuth angles relative to the aircraft; P, T and X represent the gas pressure, temperature and species content respectively; n represents the number of stratified layers; ∆d is the spatial sampling interval of LOS; L n BB is the blackbody radiance of the n-th slab; τ n denotes the transmissivity of the n-th slab.

Figure 4 .
Figure 4. Auxiliary data for simulation, (a) aircraft position (red circle), the green box is selected cloud background area; (b) spectral response functions of B7-B12 and sea surface reflectance; (c) aircraft position and pixel aggregation information for B7-12 images.

Figure 4 .
Figure 4. Auxiliary data for simulation, (a) aircraft position (red circle), the green box is selected cloud background area; (b) spectral response functions of B7-B12 and sea surface reflectance; (c) aircraft position and pixel aggregation information for B7-12 images.

Figure 5 .
Figure 5.The comparison of the background and aircraft radiance spectra, (a,b) is the comparison of simulated and measured restored radiance, and (c,d) is the comparison of target and background radiance curves.

Figure 5 .
Figure 5.The comparison of the background and aircraft radiance spectra, (a,b) is the comparison of simulated and measured restored radiance, and (c,d) is the comparison of target and background radiance curves.

Figure 6 .
Figure 6.Comparison of plume measurement and simulation and atmospheric transmittance at 20 m horizontal path.

Figure 6 .
Figure 6.Comparison of plume measurement and simulation and atmospheric transmittance at 20 m horizontal path.

Figure 7 .
Figure 7. Relative contributions for the body-leaving radiance at 2.5-13 µm.(a-c) are the contribution plots of three spatial resolutions in the daytime; (d-f) are the contribution plots in the nighttime; the blue, red, yellow, and purple areas represent the relative contribution of the background radiance, plume radiance, skin reflected radiance and skin emission radiance, respectively.

25 Figure 7 .
Figure 7. Relative contributions for the body-leaving radiance at 2.5-13 μm.(a-c) are the contribution plots of three spatial resolutions in the daytime; (d-f) are the contribution plots in the nighttime; the blue, red, yellow, and purple areas represent the relative contribution of the background radiance, plume radiance, skin reflected radiance and skin emission radiance, respectively.

Figure 8 .
Figure 8. Relative contributions for the TOA radiance at 2.5-13 μm.(a-c) are the contribution plots of three spatial resolutions in the daytime;(d-f) are the contributions plots in the nighttime; the blue, red, yellow, purple, and green areas represent the relative contribution of the background radiance, plume radiance, skin reflected radiance, skin emission radiance, and atmosphere path radiance respectively.

Figure 8 .
Figure 8. Relative contributions for the TOA radiance at 2.5-13 µm.(a-c) are the contribution plots of three spatial resolutions in the daytime;(d-f) are the contributions plots in the nighttime; the blue, red, yellow, purple, and green areas represent the relative contribution of the background radiance, plume radiance, skin reflected radiance, skin emission radiance, and atmosphere path radiance respectively.

Figure 9 .
Figure 9. B7-12 radiance comparison of aircraft and background simulated without the skin reflected radiance.The left axis is the radiance, and the right axis is the contrast, where the negative number means that the target radiance is lower than the background radiance.

Figure 9 .
Figure 9. B7-12 radiance comparison of aircraft and background simulated without the skin reflected radiance.The left axis is the radiance, and the right axis is the contrast, where the negative number means that the target radiance is lower than the background radiance.

Figure 9 .
Figure 9. B7-12 radiance comparison of aircraft and background simulated without the sk flected radiance.The left axis is the radiance, and the right axis is the contrast, where the neg number means that the target radiance is lower than the background radiance.

Figure 10 .
Figure 10.Comparison of target-background contrast ratio in daytime and nighttime.

Figure 10 .
Figure 10.Comparison of target-background contrast ratio in daytime and nighttime.
ote Sens. 2023, 15, x FOR PEER REVIEW 18 ofWith regard to aircraft space-based observation characteristics, spatial degradati is the most intuitive impact factor.Figure11shows the T-B contrast spectra at differ GSDs of 70-400 m.The spectral contrasts were shown in three segments due to the lar difference between the different spectral ranges.The contrast becomes smaller in absol terms as GSD increases, with a 97% reduction from 70 to 400 m.Because the lower spatial resolution is, the higher the background radiance contribution is in the aircr pixel, which makes the aircraft observed radiance closer to the background radiance.T contrast curves for different MTFs at the GSD of 120 m were presented in Figuredemonstratingthat T-B relative differences increase with increasing MTF.

Figure 11 .
Figure 11.Spectral contrast of target and background at different GSDs of 70-400 m.

Figure 11 .
Figure 11.Spectral contrast of target and background at different GSDs of 70-400 m.

Figure 11 .
Figure 11.Spectral contrast of target and background at different GSDs of 70-400 m.

Figure 12 .
Figure 12.Spectral contrast of target and background at different MTFs of 0.1-0.3.

Figure 12 .
Figure 12.Spectral contrast of target and background at different MTFs of 0.1-0.3.

Figure 13 .
Figure 13.Target-background contrast and the spectral response functions of VIMS B7-12 at the GSD of 120 m.

Figure 13 .
Figure 13.Target-background contrast and the spectral response functions of VIMS B7-12 at the GSD of 120 m.

Figure 14 .
Figure 14.Contrast and its variation range with different center wavelengths and FWHMs.

Figure 14 .
Figure 14.Contrast and its variation range with different center wavelengths and FWHMs.

Figure 15 .
Figure 15.Comparison of noise equivalent radiance and radiance difference between target and background.

Figure 15 .
Figure 15.Comparison of noise equivalent radiance and radiance difference between target and background.
25 µm and 6µm covers the zero axis, indicating a change in the T-B relative relationship (light-dark relationships).It is obvious that a change in the T-B light-dark relationship is not expected, which seriously affects the distinction between the target and background.

Figure 15 .
Figure 15.Comparison of noise equivalent radiance and radiance difference between targe background.

Figure 16 .
Figure 16.The variation range of target background contrast and its standard deviation curve sidering the influence of noise.

Figure 16 .
Figure 16.The variation range of target background contrast and its standard deviation curve, considering the influence of noise.

Figure 17 .
Figure 17.Atmospheric transmittance from different altitudes to the top of the atmosphere and atmospheric profile of CO2, H2O, derived from output files of the MODTRAN.

Figure 17 .
Figure 17.Atmospheric transmittance from different altitudes to the top of the atmosphere and atmospheric profile of CO 2 , H 2 O, derived from output files of the MODTRAN.

(
S Skin (L E_skin + L R_skin )ds + S NP (L E_plume + τ Plume L Nozzle )ds + S Plume −S NP (L E_plume + τ Plume L H Bkg )ds + S Bkg L H Bkg ds) + L H→TOA TOA Bkg are the radiance of the skin emission, reflected radiance, plume radiance and background radiance reaching the sensor, respectively.
pth   ⊗ R sensor + L noise TOA Plume and L

Table 1 .
Imaging information and instrument performance parameters of GF-5 VIMS..The VIMS provides images in 12 bands from the visible to thermal infrared, six of which were used for validation in this study.The central latitude and longitude of the chosen data are 32.212895°N and 126.392002°E, located in the East China Sea and imaged on 25 June 2019, with the specific parameters shown in Table1. ment

Table 1 .
Imaging information and instrument performance parameters of GF-5 VIMS.

Table 3 .
Simulation accuracy, target-background contrast calculation results.