1. Introduction
Simulation is an indispensable process in the calibration and validation of instruments or analysis of algorithms. It is also an important way to analyze the characteristics of a target in specific circumstances. Over the last two decades, digital simulation has become the major approach to simulating hyperspectral imaging along the sun–target–observer image chain [
1,
2,
3,
4]. In the process of digital simulation, solar irradiance, skylight, and reflected background radiance are considered to illuminate the target. Then, the radiance to the sensor is calculated as a combination of directly reflected radiance and upwelled radiance of atmosphere [
5,
6,
7]. Lastly, the sensor characteristics are described using common models. This is a good approach to simulating images in all kinds of radiation and imaging geometric conditions. However, its performance is limited by the knowledge of characteristics of the target, environment, and sensor.
Although the precision reflectance model of buildings and canopies [
8,
9] and the mixing model of minerals (such as the Hapke model [
10,
11]) have been developed, and the precision model of a spectrometer can also be specifically built [
12], it is still important to measure reflected radiance in a hardware-in-the-loop (HWIL) simulation when researching the influence of environment, such as heavy metal toxicity in plants, nonlinear mixing of minerals, or weathered minerals [
13,
14,
15,
16,
17,
18,
19,
20,
21,
22]. The traditional approaches include conducting flight campaigns over experimental sites [
13,
14,
15] and experiments in labs [
16,
17,
18,
19,
20,
21,
22]. The flight campaign approach has several shortcomings such as the great expense, time consumption, the complexity and difficulty of making simultaneous in-situ measurements of atmospheric parameters and ground reflective characteristics. The result of using a flight campaign is also limited by many environmental conditions such as solar zenith, weather, and visibility.
Therefore, most experiments for researching or instrument testing are accomplished in laboratories by measuring the spectral reflectance of the target. The common light source illuminating targets are halogen lamps [
16,
17,
18,
19] or halogen lamps combined with a collimator or integrating sphere [
20,
21,
22]. In this context, the indoor experiment faces some challenges. Firstly, the irradiance illuminated on the target can hardly simulate the geometric characteristics of solar irradiance and skylights. Secondly, the spectral characteristic of halogen lamps is quite different from that of solar irradiance or skylights [
23].
In order to simulate the remote multispectral imaging process indoors, an HWIL simulation facility named Imagery Simulation Facility (ISF) was developed by ITEK Optical System (Lexington, MA, USA) [
24], which includes both solar and skylight simulators. However, the spectral characteristic of the tungsten lamps selected for the skylight simulator was quite different from real skylights. Therefore, they are respectively filtered with several different filters for the spectral range from 400 nm to 2500 nm, or replaced by a different kind of lamp with accurate spectral characteristics from 400 nm to 800 nm for color film imaging simulation [
24]; this results in either of the following two problems. First, the spectral isotropy over the hemisphere is decreased by different filters. Second, the facility can only work in the visible region. Except for the above problems, the change in skylight irradiance with different solar zeniths was not reported. Therefore, the ISF facility cannot be used to test a hyperspectral system or algorithm, especially for studies on mineral identification, which are mostly based on spectral features in the short-wave infrared (SWIR) region.
In this paper, a hemispherical lamp array with 80 identical lamps is designed as a skylight simulator for the indoor simulation of hyperspectral remote sensing. A metal-halide lamp is chosen as the light source in the simulator and designed to point to the center of the array. The spatial distribution of lamps is modeled and tested considering the slightly different performances among individual lamps. A series of experiments are conducted to test the performance of the simulator. As the result, the spatial uniformity of irradiance is greater than 0.91. The spectral mismatch is about 1/243 that of the ISF. The spatial distribution of radiance can be adjusted, and the accuracy of the adjustment is greater than 0.895.
Section 2 analyzes the characteristics of the skylight and draws the requisite performance of the skylight simulator.
Section 3 indicates the development of the skylight simulator. In
Section 4, the accuracy of the simulator is tested. Conclusions are drawn in
Section 5.
2. Theory
The spectral irradiance onto the surface,
, can be divided into three parts including direct solar irradiance
, skylight irradiance
, and background radiance
[
25]:
where
is the fraction of the hemispherical sky that could be seen from the target,
is the zenith of the target, and
is the transmission from solar to the scattering volume. As shown in
Figure 1 [
25], the downwelled irradiance (skylight)
could be modeled as the integration of the directional radiance from the scattering of solar light by a small unit volume over the hemisphere area above the target [
25]:
where
is the angle between the normal to the target and the ray from the volume
,
r is the distance from the target to the volume which will be integrated from the target to the top of atmosphere (TOA), and
and
are the transmissions along the paths
and
. This reveals that
is attenuated through the sun–volume path
, scattered in a deflection of
onto the target, in the fraction of the angular scattering coefficient
, and attenuated again through the volume–target path
.
Accordingly, the skylight generally keeps three characteristics. (a) The target is illuminated by a hemispherical skylight and the spatial radiation distribution is influenced by the solar direction. (b) The total irradiance varies with the atmospheric status. (c) The total spectral irradiance arriving at the ground is spatially uniform.
In order to achieve a group of standard spectra of irradiance, MODTRAN 4.1 (Spectral Sciences Inc., United States Air Force, USA) was used to calculate the radiance from different directions and irradiance onto the ground surface for different dates and times [
26]. Ground locations were set at Hami, Xinjiang, China, which is a typical zone with many porphyry copper–gold mineralization subzones [
27]. The parameters used in MODTRAN calculation are listed in
Table 1.
As shown in
Figure 2, the skylight irradiates the target from every direction and varies with incident angle, holding the same spectral characteristic. The radiance comes to the maximum in the solar direction (zenith 17°, azimuth 180°) and the minimum appears in the opposite direction of solar incidence, as shown in
Figure 2a,b. The spectral isotropy of radiance is tested by calculating the spectral correlation coefficients between azimuthal-average spectral radiance with different zenith angles and between zenithal-average spectral radiance with different azimuth angles, respectively, as shown in
Table 2 and
Table 3. The Pearson correlation coefficient is calculated according to Equation (3), where
and
represent the radiances in band
with different zenith/azimuth angles respectively. As the correlation in
Table 2 and
Table 3 is greater than 0.980, the spectral characteristic of radiance can be summarized as stable with the change of zenith and azimuth. The radiance at 16:00 shows the same characteristic in
Figure 2c,d (the solar zenith is 55° and azimuth is 180°). Caused by the change in solar direction, the spatial distribution of the radiance changed. The correlation of the max radiance, min radiance, and middle one at 12:00 and 16:00 is also calculated, as shown in
Table 4. The spectral characteristic at different times can be defined as stable, because the correlation is greater than 0.957.
The stability of the spectral characteristic ensures the probability of adjusting the irradiance of the skylight simulator. It is obvious that the spatial uniformity of irradiance and spatial distribution of radiance onto the ground are two very important characteristics of the skylight. When adjusting the zenith distribution of the radiance, the difference in radiance between each target could be negligible. However, when adjusting the azimuth distribution of radiance, lamps need to point to different targets, respectively, to maintain the uniformity of irradiance, which will cause great differences in radiance on targets. As a result, the simulation of azimuth distribution of radiance is almost impossible. Therefore, the skylight simulator is supposed to be adjusted to simulate the zenith distribution of radiance, and maintain the spatial uniformity of spectral irradiance at the same time. The irradiance and the zenith distribution of the radiance would be similar to the MODTRAN-simulated result shown in
Figure 2.
5. Conclusions
The skylight simulator is an important facility in the HWIL simulation of remote sensing, especially for hyperspectral remote sensing. The previous facilities, including the ISF developed by ITEK Optical System, are not applicable in a hyperspectral remote sensing experiment. The two main reasons for this are the mismatch of spectral characteristics and lack of ability of adjustment for the simulation of different scenes, because of the limitations of the light source and structure. A new facility was designed and developed with a new wide-beam metal-halide lamp and hemisphere structure to spatially and spectrally simulate skylight illumination.
The performance of the proposed skylight simulator was tested using a spectrometer with different accessories in different situations. The spatial uniformity of spectral irradiance is greater than 0.91 in almost all bands from 350–2200 nm. The spectral match of the simulator with the real skylight improves by about 243 times in the visible and near infrared (VNIR) and SWIR regions compared with a halogen-lamp-based system such as ISF. The spectral isotropy over the hemisphere is greater than 0.9995. The accuracy of adjusting the irradiance of the simulator is greater than 0.895 for three modes. The accuracy of irradiance distribution along different zeniths is greater than 0.741.
To validate the spectral performance of the simulation, the spectra of irradiance, radiance reflected by targets, and reflectance of targets were collected with the skylight simulator and compared with those collected in the field. The relative deviation of irradiance is less than 0.021. The RMSRE of spectral radiance is 0.6–0.7 with the relative deviation of radiance less than 0.184. The RMSRE of spectral reflectance is less than 0.233, and the correlation is greater than 0.966. The spectral isotropy over the hemisphere is greater than 0.997.
All the above spectral accuracies in the 350–2500 nm region are about 243 times those of the latest facility. Simultaneously, the isotropy of the spectral characteristic is ensured, and the adjusting accuracies of irradiance and radiance are greater than 0.895 and 0.741, respectively.
Future work will be focused on developing new light sources working in 350–2500 nm to improve the accuracy of the spectral simulation, and, in particular, to decrease the mismatch of the spectral characteristic in the 350–600 nm region.