Next Article in Journal
Parameterized Modeling and Calibration for Orbital Error in TanDEM-X Bistatic SAR Interferometry over Complex Terrain Areas
Next Article in Special Issue
Target Classification of Similar Spatial Characteristics in Complex Urban Areas by Using Multispectral LiDAR
Previous Article in Journal
Multi-Source Hydrological Data Products to Monitor High Asian River Basins and Regional Water Security
Previous Article in Special Issue
Land Cover Classification with Multispectral LiDAR Based on Multi-Scale Spatial and Spectral Feature Selection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parameter Simulation and Design of an Airborne Hyperspectral Imaging LiDAR System

1
College of Optical and Electronic Technology, China Jiliang University, Hangzhou 310018, China
2
Anhui Institute of Optics and Fine Mechanics, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei 230031, China
3
Innovation Academy for Institute of Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China
4
State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing, Wuhan University, 129, Luoyu Road, Wuhan 430072, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(24), 5123; https://doi.org/10.3390/rs13245123
Submission received: 8 November 2021 / Revised: 8 December 2021 / Accepted: 10 December 2021 / Published: 17 December 2021
(This article belongs to the Special Issue Land Cover Classification Using Multispectral LiDAR Data)

Abstract

:
With continuous technological development, the future development trend of LiDAR in the field of remote sensing and mapping is to obtain the elevation and spectral information of ground targets simultaneously. Airborne hyperspectral imaging LiDAR inherits the advantages of active and passive remote sensing detection. This paper presents a simulation method to determine the design parameters of an airborne hyperspectral imaging LiDAR system. In accordance with the hyperspectral imaging LiDAR equation and optical design principles, the atmospheric transmission model and the reflectance spectrum of specific ground targets are utilized. The design parameters and laser emission spectrum of the hyperspectral LiDAR system are considered, and the signal-to-noise ratio of the system is obtained through simulation. Without considering the effect of detector gain and electronic amplification on the signal-to-noise ratio, three optical fibers are coupled into a detection channel, and the power spectral density emitted by the supercontinuum laser is simulated by assuming that the signal-to-noise ratio is equal to 1. The power spectral density emitted by the laser must not be less than 15 mW/nm in the shortwave direction. During the simulation process, the design parameters of the hyperspectral LiDAR system are preliminarily demonstrated, and the feasibility of the hyperspectral imaging LiDAR system design is theoretically guaranteed in combination with the design requirements of the supercontinuum laser. The spectral resolution of a single optical fiber of the hyperspectral LiDAR system is set to 2.5 nm. In the actual prototype system, multiple optical fibers can be coupled into a detection channel in accordance with application needs to further improve the signal-to-noise ratio of hyperspectral LiDAR system detection.

1. Introduction

After nearly half a century of development, Earth observation technology has played an increasingly important role in environmental monitoring, resource exploration, urban planning, forestry investigation, and other fields [1,2,3,4,5,6,7]. As the core method of Earth observation, remote sensing has gradually developed into a comprehensive technical system that covers passive multispectral imaging, hyperspectral imaging, active and passive microwave remote sensing, LiDAR remote sensing, and other methods of information acquisition. Remote sensing can be performed for aviation and Earth observation on multiple platforms, such as aerospace and ground [8,9]. Realizing full-time, full-band, high-resolution Earth observation is important in the development of remote sensing technology because it can address many problems in environmental monitoring, urban construction, and military applications. In the process of Earth observation, the simultaneous acquisition and identification of spatial elevation and spectral information is a frontier scientific issue for the development of Earth observation technology in the future [10,11,12].
Among the main remote sensing technology methods available at present, active LiDAR ranging can acquire the elevation information of ground objects, and passive hyperspectral imaging can obtain rich spectral information [13,14]. The two technical methods show unique advantages in the acquisition of ground target information, but neither can achieve the simultaneous acquisition of elevation and spectral information [15,16]. Moreover, the acquisition of the spectral information of ground objects is limited due to the principle of active single-wavelength LiDAR ranging. Passive hyperspectral imaging is greatly affected by light intensity and weather changes during the acquisition of ground object information, and its capability to acquire spatial elevation information is limited. To combine the technical advantages of active and passive detection, domestic and foreign researchers have conducted extensive work and have made many attempts [17,18,19,20], all of which have achieved good results. Relevant research has also fully verified that increasing the detection bands is crucial to improve the capability of LiDAR spectral information acquisition, especially in vegetation remote sensing monitoring, and thus shows great application potential [21,22,23]. However, for the detection requirements of different ground targets and different growth states of the same plant, other characteristic wavelengths that characterize physical properties need to be added to achieve enhanced spectral detection effects. However, when the wavelength that characterizes the properties of the ground object is increased, the corresponding wavelength band of the supercontinuum laser and the detection channel of the hyperspectral LiDAR system also need to be increased, thereby increasing the complexity and cost of the hyperspectral LiDAR system. The next step of related research should be to incorporate the radar detection effect when conducting a comprehensive evaluation to determine when to increase the number of optimal detection bands and to identify the best characteristic wavelength.
Hyperspectral LiDAR technology for Earth observation was introduced to address these issues. Researchers agree that the spectral band of Earth observation LiDAR ranges from a single wavelength to multiple spectra, which is the development direction of hyperspectral imaging [24,25,26,27]. The capability to obtain elevation and spectral information can be improved by increasing the detection spectral band of hyperspectral LiDAR. Compared with the traditional passive hyperspectral imaging technology, the active hyperspectral imaging LiDAR can obtain spectral information and elevation information at the same time with one laser pulse, realizing the perfect matching of information. However, due to the energy limitation of supercontinuum lasers and technical bottlenecks, such as the simultaneous detection and acquisition of multi-channel weak echo signals, the development of hyperspectral imaging LiDAR has been limited to a certain extent [28,29].
This study introduces the structure of an airborne hyperspectral imaging LiDAR system. The spectrum range of the system is 400–900 nm, and it can realize 50+ channels of synchronous detection and acquisition. The hyperspectral LiDAR system adopts the optical fiber array focal plane splitting method to realize the efficient coupling of the echo signals of all channels. At the same time, the weak energy of the supercontinuum laser in the shortwave band leads to a poor signal-to-noise ratio (SNR) of the echo signal. Thus, the system uses single-tube detectors. Compared with traditional line and area array detectors, the advantage of single-tube detectors is that the responsivity and gain of each channel can be adjusted independently. This study simulates the parameters of the system before it is built and examines the effects of the power of the supercontinuum laser and the design parameters of the LiDAR system on the SNR of the echo signal. The main theoretical basis of this research work is carried out around the LiDAR equation. We combined the principles of optical design and finally determined the design parameters of the system components and the transmitting power spectrum of the supercontinuum laser.

2. Description of Airborne Hyperspectral LiDAR

The development of multi-spectral LiDAR systems has verified that the spectral detection capability for ground objects can be enhanced by increasing the number of LiDAR detection bands. Therefore, further development of LiDAR systems with increased numbers of bands can achieve the accurate acquisition of spectral and elevation information of ground objects. Compared with traditional multi-spectral LiDAR, the hyperspectral LiDAR system has stronger spectral detection capabilities and can achieve spectral coverage from visible light to near-infrared by increasing the number of bands. Hyperspectral LiDAR technology based on the basic structure of the supercontinuum laser and grating spectroscopy is expected to be the main development direction in the next stage. The system will gradually develop from a theoretical prototype to an experimental prototype and finally to an engineering prototype. The application platform of the hyperspectral LiDAR system will gradually expand from the ground platform to the vehicle-mounted platform and finally to the airborne multi-platform.
As shown in Figure 1, the hyperspectral imaging LiDAR system mainly includes lasers, scanning mirrors, receiving telescopes, grating spectrometers, detector arrays, an acquisition system, and a storage system. The scanning mirror reflects the laser beam to the ground target for line scanning, and the echo signal is received by the telescope along the coaxial optical path [30]. The off-axis parabolic mirror is used as the optical receiving telescope to receive the full-band laser echo signal scattered by the ground object, and the receiving field-of-view angle of the telescope is limited by the small aperture diaphragm to suppress the background light signal entering the grating spectrometer. The echo signal passing through the aperture diaphragm is split by the grating, and the echo signals of different bands are coupled to the corresponding detection channel. The hyperspectral LiDAR system can couple the echo signals of different bands to the corresponding detectors through the optical fiber in accordance with the detection needs to realize the flexible selection of the center wavelength and bandwidth of different detection channels. The spectrum range of the echo signals is 400–900 nm. The acquisition system realizes rapid parallel acquisition of data from more than 50 detection channels and obtains the spatial elevation and spectral information of the target to be measured.
The basic working principles of hyperspectral imaging LiDAR and multi-spectral imaging LiDAR for Earth observation are similar to those of ordinary LiDAR systems. The essence of LiDAR is the interaction between light and matter. The laser emission system sends out a light signal, the laser pulse interacts with the target, and the echo signal is received and processed by the detection and acquisition system to obtain the required information. The echo signal received by LiDAR may be a reflected signal or an elastic scattering signal received by the system. In hyperspectral and multispectral LiDAR, the useful echo signals obtained are all reflected signals of ground objects.

3. Materials and Methods

The parameter simulation of the hyperspectral imaging LiDAR system is based on the theoretical model established with the traditional LiDAR system and combined with the hyperspectral imaging LiDAR equation. Then, the simulation model of the hyperspectral imaging LiDAR system is constructed. Using the design parameters of the hyperspectral imaging LiDAR system and the performance parameters of the laser and in combination with the characteristic spectrum of the ground object, the detection results of LiDAR are simulated under various conditions. The obtained results can enable LiDAR system workers to select the optimized design approach when examining the system, and they can have a complete understanding of system performance before the system prototype is completed. Therefore, simulation of the LiDAR system is an indispensable part of current LiDAR system design and research. A complete theory is available for the simulation of radar systems [31]. Regarding the simulation of the LiDAR system, current domestic and foreign scholars have also carried out extensive research and established a complete theoretical system [32,33].
Figure 2 shows a parametric simulation model of hyperspectral LiDAR. Before designing the airborne hyperspectral imaging LiDAR, we must preliminarily determine the theoretical parameters of each part of the optical components. On the basis of the traditional LiDAR model and the SNR results, the system design parameters are further determined and optimized, and the performance of the supercontinuum laser is evaluated. In the simulation of the design parameters of the hyperspectral LiDAR system, the ground object is assumed to be soil and red leaves, and the echo signal intensity after the laser pulse signal transmitted by the hyperspectral lidar passes through the ground object is studied. At the same time, the influence of the atmosphere and sky background light on SNR is considered.
The theoretical model of hyperspectral LiDAR is similar to the traditional single-wavelength LiDAR, so the feasibility demonstration of the hyperspectral LiDAR design scheme is similar to that of single-wavelength LiDAR. In the simulation process, we use the SNR of hyperspectral LiDAR as the basis to demonstrate the feasibility of the design parameters. On the basis of the simulation results of SNR, we derive the performance indicators of the laser and further optimize the design parameters of hyperspectral LiDAR.

3.1. Hyperspectral Waveform Modeling

Airborne hyperspectral imaging LiDAR uses a supercontinuum laser as a light source to transmit laser pulse signals in a wide spectrum to the measured target. An optical receiving and detecting system is adopted to receive and detect laser echo signals in a wide spectrum scattered by the target before obtaining the elevation and spectral information of the target. Generally, the intensity of the echo signal received by hyperspectral LiDAR is mainly affected by three aspects. First, the echo signal strength of LiDAR is determined by the interaction between the laser beam and the ground target. For hyperspectral imaging LiDAR, the reflectivity of the ground target affects the strength of the echo signal. The second aspect is the influence of the medium between hyperspectral LiDAR and the ground target on the echo signal, such as the attenuation of the echo signal by the atmosphere. Lastly, the strength of the echo signal of hyperspectral LiDAR is related to the configuration of the LiDAR system and the technical factors of the system, such as the influence of the transmittance of the LiDAR system on the strength of the echo signal. The strength of the echo signal directly determines whether hyperspectral LiDAR can detect ground objects effectively.
As shown in Figure 3, the single-pulse energy emitted by the laser is E T , and the distance between the receiving system and the ground target is z . The laser emits a pulse signal that is incident on the surface of the ground target, and the pulse signal uniformly scatters in the hemispherical space in a Lambertian form. Without considering the laser energy loss, the energy density of the scattered laser pulse is defined as D s .
D s = E T 2 π z 2
The intensity of the echo signal received by the hyperspectral LiDAR system is defined as P R , the reflectivity of the ground target is ρ , and the receiving aperture of the telescope is D R .
P R = E T π 0.5 D R 2 2 π z 2 ρ
During the airborne LiDAR flight test, the optical components of the LiDAR system experience a certain loss of laser pulse energy. The coupling efficiency among various parts of the optical components is also important. In the process of data processing, the attenuation effect of the laser pulse energy by the atmosphere needs to be considered. In the theoretical derivation process of the LiDAR equation, the system attenuation coefficient is recorded as η s y s , and the atmospheric transmission factor is recorded as η a t m . Equation (3) is then obtained.
P R = E T π 0.5 D R 2 2 π z 2 ρ η s y s η a t m
In the design process of airborne multi-channel hyperspectral imaging LiDAR, the working mechanism of single-shot laser pulse emission and multi-channel synchronous acquisition of echo signals is usually adopted. First, after the single pulse emitted by the laser is scattered by the ground object, the LiDAR receiving system splits the wide-spectrum echo signal through the beam splitter element. The beam splitter element is generally a grating or a Fabry–Perot interferometer. Then, the echo signal is coupled to the corresponding detection channel through light transmission. In the hyperspectral LiDAR system discussed in this article, we use optical fiber as a light transmission device to couple the echo signals of different spectral bands to the corresponding single-tube detectors.
In accordance with the actual working mode of the hyperspectral imaging LiDAR system discussed in this article and in combination with Equation (3), we obtain the echo signal intensity of the multi-channel hyperspectral LiDAR equation, which can be further defined by Equation (4) [34].
P R λ , z = ρ 0 η λ Δ λ β G λ D R 2 8 z 2 ε λ T a t m λ , 0 , z 2 ,
where λ is the central wavelength of each channel obtained during spectral calibration (the determination of center wavelength λ is one of the core contents of the calibration of the radar system), P R λ , z is the optical power of the return signal received by the channel with a center wavelength of λ , ρ 0 is the average spectral power density produced by the laser, η λ is the spectral power density distribution function normalized by the average spectral power density of the laser, Δ λ is the spectral bandwidth corresponding to one channel, and β T λ is the reflectivity of the target. Moreover, D R is the effective aperture of the receiving telescope, z is the distance between the LiDAR and the measured object and can be measured in real time through the ranging channel, ε λ is the optical efficiency of the LiDAR system, and T a t m λ , 0 , z is the transmittance of the atmosphere between the LiDAR and the measured surface at wavelength λ .
In the design process of the LiDAR system, SNR is an issue that cannot be ignored. Generally, the process includes sky background light noise, stray light noise of the LiDAR system itself, and dark current noise and thermal noise in the photoelectric detection system. If the noise cannot be effectively suppressed, the SNR of the echo signal will be reduced, which will greatly affect the measurement accuracy of the hyperspectral LiDAR system. We focus on the influence of sky background light noise on SNR.
When the LiDAR system is working, the detection performance at night is stronger than that in the day due to the influence of sky background light. To ensure that the hyperspectral LiDAR system can detect ground objects throughout the day, its design needs to ensure that the radar can work normally during the day. Assuming that sky background light is uniformly scattered in all directions in space, the actual sky background light intensity received by the LiDAR system can be understood as the sum of the sky background light intensity within the range of the telescope projected to the ground.
Therefore, the daytime background light power received by the LiDAR system can be defined by Equation (5).
B R λ , z = L λ β G λ π 2 θ 2 D R 2 16 Δ λ ε λ T λ , 0 , z ,
where B R λ , z is the sky background light power scattered by the ground surface and received by each channel of the hyperspectral LiDAR, L λ is the monochromatic sky background light radiance received by each detection channel, and θ is the field of view of the receiving telescope.
By referring to the traditional LiDAR SNR calculation method [35] combined with Equations (4) and (5), we can use Equation (6) to define the SNR of the echo signal received by the hyperspectral LiDAR from the surface target object scattering.
S N R λ , z = N P R λ , z P R λ , z + B R λ , z
In the actual flight operation of hyperspectral LiDAR, when the flying altitude is Z from the ground target, the SNR of the surface-scattered echo signal received by the channel corresponding to the center wavelength of λ is defined as S N R λ , z . N is the cumulative laser pulse number.

3.2. Parameter Simulation of Hyperspectral LiDAR

In accordance with Equation (6), we can simulate the SNR of the hyperspectral LiDAR system, which can be used as an important theoretical basis for determining the design parameters of the LiDAR system. Table 1 shows the design parameters to be determined for the hyperspectral LiDAR system. The spectral resolution of the hyperspectral LiDAR system is set to 2.5 nm, and the clear aperture of the telescope is set to 0.1 m. The working band of hyperspectral LiDAR is set to 400–900 nm. The optical efficiency of the hyperspectral LiDAR system is also an important design index. We expect the optical efficiency of the hyperspectral LiDAR system to be 0.3. We optimize and confirm the main technical parameters of the LiDAR system through simulation to ensure that the design parameters meet the pre-designed performance requirements of hyperspectral LiDAR.
The spectral power density in the 400–900 nm range is given in Figure 4. The spectral power density distribution function is obtained by normalizing the power spectrum of the laser output. The power density spectrum of the laser is provided by the manufacturer, and the power density spectrum corresponding to different wavelengths in the simulation process is weaker than the power density spectrum of the actual laser.
Sky brightness usually refers to the part of solar radiation energy absorbed and scattered by the atmosphere and the spatial distribution of the radiance of scattered light formed within the atmosphere [36]. The brightness of the sky is generally determined by the position of the sun and the characteristics of the scatterers in the atmosphere, and it usually involves atmospheric molecules and aerosol particles in the atmosphere. Sky brightness is the radiation caused by these factors. Sky brightness is an important direction of research on the detection and recognition of targets in space [37], the tracking and imaging of targets in space [38], and the inversion of the optical properties of aerosols and clouds [39]. Figure 5 shows the spectral curve of sky background light radiance. The sky background light radiance mentioned in this article is based on the American ASTM G173-03 model [40].
In the process of simulating the parameters of the LiDAR system to obtain SNR, vegetation and soil are used as the main ground objects for analysis, and vegetation is mainly represented by red leaves. As shown in Figure 6, the reflectance of the ground targets in the 400–900 nm spectral band is selected. The reflectance spectrum curves of the red leaves and soil in Figure 6 are all measured results provided by the State Key Laboratory of Information Engineering for Surveying, Mapping, and Remote Sensing. In actual flight operations, hyperspectral imaging LiDAR is mainly used for the remote sensing observation of soil and trees.
For airborne hyperspectral LiDAR, the influence of atmospheric transmittance on the measurement results of the LiDAR system is crucial. In the field of Earth observation and remote sensing measurement, the absorption of broad-spectrum laser energy by atmospheric molecules directly affects the final measurement results of the LiDAR system. As shown in Figure 7, the atmospheric transmittance curve corresponding to 400–900 nm is selected during the simulation process. This curve is based on the CART mode provided by the Hefei Institute of Material Science, Chinese Academy of Sciences [41].

4. Results and Discussion

We simulate the design parameters of the system on the basis of the hyperspectral LiDAR equation. For the measurement of the spectral data of each channel, only the total energy value received by each channel needs to be considered in the simulation process. The spectral range of the hyperspectral LiDAR system is 400–900 nm, the detection channel is set to 50+, and the spectral resolution of a single fiber is set to 2.5 nm.
According to Equation (5), without considering the influence of atmospheric transmittance on the measurement results, the brightness spectrum of the sky background light received by the hyperspectral LiDAR system is mainly affected by sky background light radiance and ground target reflectivity. As shown in Figure 8, assuming that the ground objects are red leaves and soil, the background light intensity spectra of the red leaves and soil are obtained with Equation (5). In the shortwave direction, because the reflectivity of red leaves is lower than that of soil, the sky background light intensity obtained by the hyperspectral LiDAR system when the ground target is red leaves is lower than that of soil. In the long-wave direction, due to the sharp increase in the reflectivity of red leaves, the corresponding sky background light intensity is considerably higher than that of soil.
In accordance with Equation (4), we obtain the power spectrum of the echo signal received by the LiDAR system through simulation. As shown in Figure 9, without considering the optical efficiency of the hyperspectral LiDAR system, the intensity of the echo signal is mainly affected by the laser power. Atmospheric transmittance is another influencing factor. Given that atmospheric transmittance has a sharply attenuated peak at 700–800 nm, the corresponding laser echo signal power also shows a tendency to decrease sharply in this spectral band. Comparison of the simulation results in Figure 8 and Figure 9 indicates that in the daytime, the intensity of the echo signal received by the LiDAR system is considerably weaker than the intensity of the sky background light.
Analysis of the simulation results in Figure 8 shows that the theoretical calculation value of sky background light power is in the order of 10−12 (PW). However, in actual operation, the suppression of daytime background light by the hyperspectral LiDAR prototype is generally in the order of 10−9 (nW). It is a difficult problem that needs to be solved in the design process of the hyperspectral LiDAR system, and emphasis should be placed on the interference of background light during the day.
In accordance with Equation (6), the laser power received by the LiDAR system and the intensity of sky background light are considered comprehensively. Figure 10 presents the numerical simulation results of the SNR of each channel when hyperspectral LiDAR works in the daytime. Under the influence of laser emission power, the SNR of the hyperspectral LiDAR system shows an upward trend in the shortwave direction and reaches a peak in the 600–650 nm spectral band. In the long-wave direction, the overall SNR of the hyperspectral LiDAR system shows a downward trend. The simulation results of the hyperspectral LiDAR system confirm that the design parameters are feasible. When the LiDAR system is flying at an airborne altitude of 500 m, the hyperspectral LiDAR system can complete the measurement of the spectral information of ground targets during daytime.
However, in the SNR simulation process of the hyperspectral LiDAR system, the influence of detector gains and electronic amplification on SNR is not considered. In the actual working process of the LiDAR system, the sky background light signal and the scattered echo signal undergo the same photoelectric conversion and signal amplification process. Therefore, the relationship between the root mean square of the scattered echo signal and its standard deviation must be the same as that of the sky background light signal. Moreover, when the LiDAR system is working at daytime, the strong sky background light signal should be one of the main considerations.
As presented in Figure 8, the theoretically calculated value of the background light power received by LiDAR is in the order of 10−12. In actual operation, the suppression of daytime background light is generally in the order of 10−9. Therefore, in the SNR simulation process of the LiDAR system, the SNR must be calculated when the background light power received by the LiDAR is increased by 1000 times.
Figure 11 presents the numerical simulation results of the SNR of each channel of hyperspectral LiDAR after the sky background light signal is amplified by 1000 times. The SNR results of all channels show a significant downward trend. The SNR of several detection channels is lower than 1, and the maximum SNR of detection channels is close to 4. In the data processing of the echo signal, a too low SNR increases the difficulty of data processing. Therefore, improving the SNR of the LiDAR system is one of the problems that all scientific researchers need to solve in the process of designing and processing the LiDAR system.
During the preliminary design of the LiDAR system, the spectral resolution of a single optical fiber is set to 2.5 nm. A higher spectral resolution can enhance the LiDAR system’s capability to distinguish ground objects, but the SNR of each detection channel cannot be guaranteed. Therefore, to improve the SNR of the receiving system as much as possible while satisfying the requirements of the LiDAR system for spectral resolution, multiple optical fibers can be coupled into one detection channel in the production process of the LiDAR prototype system. The SNR of the echo signal is improved in the process of reducing the spectral resolution of the LiDAR system.
In the design process of the hyperspectral LiDAR prototype, three optical fibers are coupled into one detection channel under the precondition that the spectral resolution can be distinguished from the ground objects. Therefore, the simulation process of the SNR of a single detection channel can be regarded as the superposition of the echo signals of three laser pulses. Figure 12 shows the power spectrum density of the laser emission obtained by simulation when the SNR of the echo optical signal of three laser pulses received by a single fiber is equal to 1. Under the premise of ensuring that the SNR of the LiDAR system meets the requirements, we need to increase the laser’s emission power, especially in the shortwave direction, and we need to ensure that the laser’s emission power spectral density is not less than 15 mW/nm.
The emission power of supercontinuum lasers is one of the important factors that determine the development of airborne hyperspectral imaging LiDAR. A laser with high emission power can make the structure of the hyperspectral LiDAR system compact, and a lightweight, compact LiDAR system is crucial for the airborne platform. If the emission power of a supercontinuum laser is not improved, then the receiving aperture of the hyperspectral LiDAR system must be increased. An extremely large receiving aperture increases the volume of the LiDAR system, which is unfavorable for the narrow space of the airborne platform. Moreover, an extremely large receiving aperture is not conducive to the suppression of sky background light by the hyperspectral LiDAR system, which has a negative impact on the SNR of the system. Therefore, when optimizing the optical–mechanical structure design of hyperspectral LiDAR, how to increase the emission power of the supercontinuum laser is another problem that researchers in this field have been trying to solve.
According to the simulation results of the SNR of the hyperspectral LiDAR system, the energy spectral distribution of the supercontinuum laser presents a nonlinear effect, especially in the shortwave direction that requires high laser energy. On the one hand, we need to increase the emission power of the supercontinuum laser. On the other hand, we need to improve the detection capabilities of LiDAR systems. In the design process of the LiDAR prototype, we refer to grating spectrometer technology and the high-efficiency coupling technology of optical fiber arrays. At the same time, high-sensitivity APD and PMT detectors are used to detect echo signals. Compared with traditional line and area array detectors, single-tube detectors present outstanding advantages. Single-tube detectors can independently control the bandwidth, sensitivity, and gain, thereby further improving the detection capability of the hyperspectral LiDAR system for wide-spectrum echo signals. In accordance with the spectral scattering characteristics of the ground target, the center wavelength and bandwidth of different detection channels can be optimized, and multiple optical fibers can be coupled into one detection channel to realize the flexible selection of the center wavelength and bandwidth of different detection channels.
Furthermore, under the premise of meeting the spectral resolution requirements of the hyperspectral LiDAR system, coupling multiple optical fibers into a detection channel can ensure that the hyperspectral LiDAR system can detect ground targets under harsh environmental conditions. In particular, it improves the detection capability of the hyperspectral LiDAR system for weak echo signals during daytime.
Background light suppression technology, including the narrow field of view receiving technology of the LiDAR system and the numerical aperture matching of the entire LiDAR system, must be adopted to further improve the detection capability of the hyperspectral LiDAR system in the daytime. In the LiDAR optical receiving system, all optical devices, such as telescopes, grating spectrometers, and optical fibers, have corresponding numerical apertures. The numerical aperture determines the size of the incident angle that the LiDAR system can receive. An extremely large numerical aperture causes excess background light signals to enter the LiDAR optical receiving system, and an extremely small numerical aperture causes part of the echo signals to not fully enter the LiDAR optical receiving system. Therefore, in the optical design of a hyperspectral LiDAR system, we need to ensure that the numerical apertures of all optical components match. At the same time, under the premise of ensuring the effective reception of laser echo signals, sky background light and stray light signals should be received as little as possible (or not at all) to improve the SNR of the hyperspectral LiDAR system for echo signal detection.

5. Conclusions

This work introduces the basic structure of airborne hyperspectral imaging LiDAR. A core structure based on optical fiber array focal plane spectroscopy is proposed to examine the characteristics of the full waveform sampling of hyperspectral imaging LiDAR. The hyperspectral imaging LiDAR equation is derived, and the design parameters of the LiDAR system are simulated based on the equation. Regardless of the influence of detector gain and electronic amplification on SNR, the simulation process uses a single fiber to receive three laser pulse echo optical signals at the same time as a precondition. When the obtained SNR is equal to 1, the simulation attains the power spectral density of the laser emission. The simulation results show that in the design process of a hyperspectral LiDAR system, we need to increase the emission power of the laser, especially in the shortwave direction, and we need to ensure that the power spectral density emitted by the laser is not less than 15 mW/nm. At the same time, in the design process of a hyperspectral LiDAR system, we should focus on the suppression of stray light, especially the suppression of sky background light during daytime. In consideration of the difficulty of detecting echo signals during daytime in the field of Earth observation, background light combined suppression technology is proposed. The intensity of sky background light can be effectively suppressed on the premise of ensuring high optical efficiency of the LiDAR optical system, thereby improving the SNR of the echo signal of the receiving system.
In this study, a method to demonstrate the design parameters of hyperspectral imaging LiDAR is proposed. The parameters obtained from the simulation are applied to the design of a hyperspectral imaging LiDAR prototype, which is about to enter the flight test stage. For our next step, flight test data will be analyzed to further demonstrate the reliability of the design parameters. We believe that this method can perform a preliminary demonstration of all LiDAR design parameters. The development of this kind of hyperspectral imaging LiDAR is expected to usher in a new era in the field of Earth observation and remote sensing.

Author Contributions

All authors made significant contributions to the work. Conceptualization, L.Q., D.W. and D.L.; formal analysis, L.Q., D.W., D.L. and S.S. (Shalei Song); funding acquisition, D.L. and W.G.; methodology, L.Q., D.W., S.S. (Shalei Song), and S.S. (Shuo Shi); project administration, D.L., D.W., S.S. (Shalei Song), and S.S. (Shuo Shi); software, D.W. and L.Q.; validation, L.Q., D.W., S.S. (Shalei Song), and S.S. (Shuo Shi); writing—original draft, L.Q.; writing—review and editing, L.Q., D.W., D.L. and L.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key R&D Program of China (Grant No. 2018YFB0504503).

Acknowledgments

The authors thank the Lidar Studio of Anhui Institute of Optics and Fine Mechanics and the State Key Laboratory of Surveying, Mapping and Remote Sensing Information Engineering of Wuhan University for their project support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, L.X.; Coops, N.C.; Aven, N.W.; Pang, Y. Mapping urban tree species using integrated airborne hyperspectral and LiDAR remote sensing data. Remote Sens. Environ. 2017, 200, 170–182. [Google Scholar] [CrossRef]
  2. Kuras, A.; Brell, M.; Rizzi, J.; Burud, I. Hyperspectral and Lidar Data Applied to the Urban Land Cover Machine Learning and Neural-Network-Based Classification: A Review. Remote Sens. 2021, 13, 3393. [Google Scholar] [CrossRef]
  3. Li, W.; Niu, Z.; Sun, G.; Gao, S.; Wu, M.Q. Deriving backscatter reflective factors from 32-channel full-waveform LiDAR data for the estimation of leaf biochemical contents. Opt. Express 2016, 24, 4771–4785. [Google Scholar] [CrossRef]
  4. Marrs, J.; Ni-Meister, W. Machine Learning Techniques for Tree Species Classification Using Co-Registered LiDAR and Hyperspectral Data. Remote Sens. 2019, 11, 819. [Google Scholar] [CrossRef] [Green Version]
  5. Lee, Z.P.; Churnside, J.; Mao, Z.H.; Wu, S.H.; Zibordi, G. Active and passive optical remote sensing of the aquatic environment: Introduction to the feature issue. Appl. Opt. 2020, 59, APS1–APS2. [Google Scholar] [CrossRef] [PubMed]
  6. Qian, L.Y. Focal Plane Spectroscopic Technique of Optical Fiber Array for Hyperspectral Imaging Lidar. PhD’s Thesis, University of Science and Technology of China, Hefei, China, 2021.
  7. Sankey, T.; Donager, J.; McVay, J.; Sankey, J.B. UAV lidar and hyperspectral fusion for forest monitoring in the southwestern USA. Remote Sens. Environ. 2017, 195, 30–43. [Google Scholar] [CrossRef]
  8. Kaasalainen, S.; Lindroos, T.; Hyyppa, J. Toward hyperspectral lidar: Measurement of spectral backscatter intensity with a supercontinuum laser source. IEEE Geosci. Remote Sens. Lett. 2007, 4, 211–215. [Google Scholar] [CrossRef]
  9. Manzo, C.; Valentini, E.; Taramelli, A.; Filipponi, F.; Disperati, L. Spectral characterization of coastal sediments using Field Spectral Libraries, Airborne Hyperspectral Images and Topographic LiDAR Data (FHyL). Int. J. Appl. Earth Obs. Geoinf. 2015, 36, 54–68. [Google Scholar] [CrossRef]
  10. Chen, B.W.; Shi, S.; Sun, J.; Gong, W.; Yang, J.; Du, L.; Guo, K.H.; Wang, B.H.; Chen, B.W. Hyperspectral lidar point cloud segmentation based on geometric and spectral information. Opt. Express 2019, 27, 24043–24059. [Google Scholar] [CrossRef]
  11. Gong, W.; Sun, J.; Shi, S.; Yang, J.; Du, L.; Zhu, B.; Song, S.L. Investigating the Potential of Using the Spatial and Spectral Information of Multispectral LiDAR for Object Classification. Sensors 2015, 15, 21989–22002. [Google Scholar] [CrossRef] [Green Version]
  12. Qian, L.Y.; Wu, D.C.; Zhou, X.J.; Zhong, L.J.; Wei, W.; Wang, Y.J.; Shi, S.; Song, S.L.; Gong, W.; Liu, D. Optical system design for a hyperspectral imaging lidar using supercontinuum laser and its preliminary performance. Opt. Express 2021, 29, 17542–17553. [Google Scholar] [CrossRef]
  13. Tyo, J.S.; Goldstein, D.L.; Chenault, D.B.; Shaw, J.A. Review of passive imaging polarimetry for remote sensing applications. Appl. Opt. 2006, 45, 5453–5469. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Hakala, T.; Suomalainen, J.; Kaasalainen, S.; Chen, Y.W. Full waveform hyperspectral LiDAR for terrestrial laser scanning. Opt. Express 2012, 20, 7119–7127. [Google Scholar] [CrossRef]
  15. Chen, B.W.; Shi, S.; Gong, W.; Sun, J.; Chen, B.W.; Du, L.; Yang, J.; Guo, K.H.; Zhao, X.M. True-Color Three-Dimensional Imaging and Target Classification BASED on Hyperspectral LiDAR. Remote Sens. 2019, 11, 1541. [Google Scholar] [CrossRef] [Green Version]
  16. Pan, Z.G.; Glennie, C.; Fernandez-Diaz, J.C.; Starek, M. Comparison of bathymetry and seagrass mapping with hyperspectral imagery and airborne bathymetric lidar in a shallow estuarine environment. Int. J. Remote Sens. 2016, 37, 516–536. [Google Scholar] [CrossRef]
  17. Brell, M.; Segl, K.; Guanter, L.; Bookhagen, B. Hyperspectral and Lidar Intensity Data Fusion: A Framework for the Rigorous Correction of Illumination, Anisotropic Effects, and Cross Calibration. IEEE Trans. Geosci. Remote Sens. 2017, 55, 2799–2810. [Google Scholar] [CrossRef] [Green Version]
  18. Anderson, J.E.; Plourde, L.C.; Martin, M.E.; Braswell, B.H.; Smith, M.L.; Dubayah, R.O.; Hofton, M.A.; Blair, J.B. Integrating waveform lidar with hyperspectral imagery for inventory of a northern temperate forest. Remote Sens. Environ. 2008, 112, 1856–1870. [Google Scholar] [CrossRef]
  19. Heiskanen, J.; Adhikari, H.; Piiroinen, R.; Packalen, P.; Pellikka, P.K.E. Do airborne laser scanning biomass prediction models benefit from Landsat time series, hyperspectral data or forest classification in tropical mosaic landscapes? Int. J. Appl. Earth Obs. Geoinf. 2019, 81, 176–185. [Google Scholar] [CrossRef]
  20. Hui, Z.Y.; Wu, B.P.; Hu, Y.J.; Ziggah, Y.Y. Improved progressive morphological filter for digital terrain model generation from airborne lidar data. Appl. Opt. 2017, 56, 9359–9367. [Google Scholar] [CrossRef]
  21. Shoot, C.; Andersen, H.E.; Moskal, L.M.; Babcock, C.; Morton, D.C.J.R.S. Classifying Forest Type in the National Forest Inventory Context with Airborne Hyperspectral and Lidar Data. Remote Sens. 2021, 13, 1863. [Google Scholar] [CrossRef]
  22. Alonzo, M.; Bookhagen, B.; Roberts, D.A. Urban tree species mapping using hyperspectral and lidar data fusion. Remote Sens. Environ. 2014, 148, 70–83. [Google Scholar] [CrossRef]
  23. Qu, Y.; Shaker, A.; Korhonen, L.; Silva, C.A.; Song, J.J.R.S. Direct Estimation of Forest Leaf Area Index based on Spectrally Corrected Airborne LiDAR Pulse Penetration Ratio. Remote Sens. 2020, 12, 217. [Google Scholar] [CrossRef] [Green Version]
  24. Chen, B.W.; Shi, S.; Gong, W.; Zhang, Q.J.; Yang, J.; Du, L.; Sun, J.; Zhang, Z.B.; Song, S.L. Multispectral LiDAR Point Cloud Classification: A Two-Step Approach. Remote Sens. 2017, 9, 373. [Google Scholar] [CrossRef] [Green Version]
  25. Chen, Y.W.; Raikkonen, E.; Kaasalainen, S.; Suomalainen, J.; Hakala, T.; Hyyppa, J.; Chen, R.Z. Two-channel Hyperspectral LiDAR with a Supercontinuum Laser Source. Sensors 2010, 10, 7057–7066. [Google Scholar] [CrossRef] [Green Version]
  26. Dalponte, M.; Bruzzone, L.; Gianelle, D. Tree species classification in the Southern Alps based on the fusion of very high geometrical resolution multispectral/hyperspectral images and LiDAR data. Remote Sens. Environ. 2012, 123, 258–270. [Google Scholar] [CrossRef]
  27. Gong, W.; Song, S.; Zhu, B.; Shi, S.; Li, F.; Cheng, X. Multi-wavelength canopy LiDAR for remote sensing of vegetation: Design and system performance. J. Photogramm. Remote. Sens. 2012, 69, 1–9. [Google Scholar] [CrossRef]
  28. Zhang, C.L.; Sun, X.J.; Zhang, R.W.; Zhao, S.J.; Lu, W.; Liu, Y.W.; Fan, Z.Q. Impact of solar background radiation on the accuracy of wind observations of spaceborne Doppler wind lidars based on their orbits and optical parameters. Opt. Express 2019, 27, A936–A952. [Google Scholar] [CrossRef] [PubMed]
  29. Ceolato, R.; Riviere, N.; Hespel, L. Reflectances from a supercontinuum laser-based instrument: Hyperspectral, polarimetric and angular measurements. Opt. Express 2012, 20, 29413–29425. [Google Scholar] [CrossRef]
  30. Qian, L.Y.; Wu, D.C.; Liu, D.; Zhou, X.J.; Wei, W.; Zhong, L.J.; Wang, W.J.; Wang, Y.J.; Gong, W. Analysis and Design of Hyperspectral Imaging Lidar Scanning Mirror. Acta Opt. Sin. 2021, 41, 232–238. [Google Scholar]
  31. Deng, B.; Li, X.A.; Wang, H.Q.; Qin, Y.L.; Wang, J.T. Fast Raw-Signal Simulation of Extended Scenes for Missile-Borne SAR With Constant Acceleration. IEEE Geosci. Remote Sens. Lett. 2011, 8, 44–48. [Google Scholar] [CrossRef]
  32. Heaps, W.S.; Burris, J.; French, J.A. Lidar technique for remote measurement of temperature by use of vibrational-rotational Raman spectroscopy. Appl. Opt. 1997, 36, 9402–9405. [Google Scholar] [CrossRef]
  33. Muhadi, N.A.; Abdullah, A.F.; Bejo, S.K.; Mahadi, M.R.; Mijic, A. The Use of LiDAR-Derived DEM in Flood Applications: A Review. Remote Sens. 2020, 12, 2308. [Google Scholar] [CrossRef]
  34. Qian, L.Y.; Wu, D.C.; Zhou, X.J.; Zhong, L.J.; Wei, W.; Wang, W.J.; Wang, Y.J.; Gong, W.; Liu, D. Calibration of a Hyperspectral-Imaging Lidar for earth Surface Identification and Classification. Acta Opt. Sin. 2020, 40, 1601001. [Google Scholar]
  35. Zhang, Z.Y. Research on the control and inversion algorithm of the scanning lidar. Master’s Thesis, University of Science and Technology of China, Hefei, China, 2019. [Google Scholar]
  36. Arshinov, Y.; Bobrovnikov, S.; Serikov, I.; Ansmann, A.; Wandinger, U.; Althausen, D.; Mattis, I.; Optics, D.M.J.A. Daytime operation of a pure rotational Raman lidar by use of a Fabry-Perot interferometer. Appl. Opt. 2005, 44, 3593–3603. [Google Scholar] [CrossRef] [PubMed]
  37. Hassebo, Y.Y.; Gross, B.; Oo, M.; Moshary, F.; Ahmed, S. Polarization-discrimination technique to maximize the lidar signal-to-noise ratio for daylight operations. Appl. Opt. 2006, 45, 5521–5531. [Google Scholar] [CrossRef]
  38. Ma, Y.; Anderson, J.; Crouch, S.; Shan, J. Moving Object Detection and Tracking with Doppler LiDAR. Remote Sens. 2019, 11, 1154. [Google Scholar] [CrossRef] [Green Version]
  39. Chiu, J.C.; Marshak, A.; Wiscombe, W.J.; Valencia, S.C.; Welton, E.J. Cloud optical depth retrievals from solar background “signals” of micropulse lidars. IEEE Geosci. Remote Sens. Lett. 2007, 4, 456–460. [Google Scholar] [CrossRef]
  40. Baguckis, A.; Novickovas, A.; Mekys, A.; Tamosiunas, V. Compact hybrid solar simulator with the spectral match beyond class A. J. Photonics Energy 2016, 6, 10. [Google Scholar] [CrossRef]
  41. Cui, S.C.; Zhen, X.B.; Wang, Z.; Yang, S.Z.; Zhu, W.Y.; Li, X.B.; Huang, H.H.; Wei, H.L. Toward a new radiative-transfer-based model for remote sensing of terrestrial surface albedo. Opt. Lett. 2015, 40, 3842–3845. [Google Scholar] [CrossRef]
Figure 1. Theoretical model of the hyperspectral LiDAR system.
Figure 1. Theoretical model of the hyperspectral LiDAR system.
Remotesensing 13 05123 g001
Figure 2. Simulation model of hyperspectral LiDAR parameters.
Figure 2. Simulation model of hyperspectral LiDAR parameters.
Remotesensing 13 05123 g002
Figure 3. Schematic of laser pulse detection.
Figure 3. Schematic of laser pulse detection.
Remotesensing 13 05123 g003
Figure 4. Normalized power density spectrum distribution function of the laser output.
Figure 4. Normalized power density spectrum distribution function of the laser output.
Remotesensing 13 05123 g004
Figure 5. Spectral curve of the radiance of sky background light.
Figure 5. Spectral curve of the radiance of sky background light.
Remotesensing 13 05123 g005
Figure 6. (a) Reflectance spectrum of red leaves. (b) Reflectance spectrum of soil.
Figure 6. (a) Reflectance spectrum of red leaves. (b) Reflectance spectrum of soil.
Remotesensing 13 05123 g006
Figure 7. Atmospheric transmittance spectrum curve at an altitude of 500 m obtained in CART mode.
Figure 7. Atmospheric transmittance spectrum curve at an altitude of 500 m obtained in CART mode.
Remotesensing 13 05123 g007
Figure 8. Power spectrum of sky background light obtained by simulation.
Figure 8. Power spectrum of sky background light obtained by simulation.
Remotesensing 13 05123 g008
Figure 9. Power spectrum of the echo signal obtained by simulation.
Figure 9. Power spectrum of the echo signal obtained by simulation.
Remotesensing 13 05123 g009
Figure 10. Signal-to-noise ratio of hyperspectral LiDAR.
Figure 10. Signal-to-noise ratio of hyperspectral LiDAR.
Remotesensing 13 05123 g010
Figure 11. Signal-to-noise ratio obtained after amplifying the sky background light.
Figure 11. Signal-to-noise ratio obtained after amplifying the sky background light.
Remotesensing 13 05123 g011
Figure 12. Emission power spectrum of the supercontinuum laser.
Figure 12. Emission power spectrum of the supercontinuum laser.
Remotesensing 13 05123 g012
Table 1. Parameters of SNR simulation calculation.
Table 1. Parameters of SNR simulation calculation.
Model Inputs IncludedIndexSimulation Input Value
λ Width of the spectrum received by the hyperspectral LiDAR400–900 nm
z Flight altitude of the hyperspectral LiDAR500 m
ρ λ Power density spectrum distribution curve of the laser outputSimulation calculation curve
Δ λ Spectral resolution of a single fiber2.5 nm
I B G λ Radiance of monochromatic sky background lightASTM G173-03
θ Full angle of the receiving field1 mrad
β G λ Reflectivity of the ground targetReflectance of soil and red leaves
D R Clear aperture of the receiving telescope0.1 m
ε λ Optical efficiency of the hyperspectral LiDAR system0.3
T λ , 0 , z Atmospheric transmittance between the hyperspectral LiDAR and the ground targetModel of CART radiation transmission
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Qian, L.; Wu, D.; Liu, D.; Song, S.; Shi, S.; Gong, W.; Wang, L. Parameter Simulation and Design of an Airborne Hyperspectral Imaging LiDAR System. Remote Sens. 2021, 13, 5123. https://doi.org/10.3390/rs13245123

AMA Style

Qian L, Wu D, Liu D, Song S, Shi S, Gong W, Wang L. Parameter Simulation and Design of an Airborne Hyperspectral Imaging LiDAR System. Remote Sensing. 2021; 13(24):5123. https://doi.org/10.3390/rs13245123

Chicago/Turabian Style

Qian, Liyong, Decheng Wu, Dong Liu, Shalei Song, Shuo Shi, Wei Gong, and Le Wang. 2021. "Parameter Simulation and Design of an Airborne Hyperspectral Imaging LiDAR System" Remote Sensing 13, no. 24: 5123. https://doi.org/10.3390/rs13245123

APA Style

Qian, L., Wu, D., Liu, D., Song, S., Shi, S., Gong, W., & Wang, L. (2021). Parameter Simulation and Design of an Airborne Hyperspectral Imaging LiDAR System. Remote Sensing, 13(24), 5123. https://doi.org/10.3390/rs13245123

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop