Next Article in Journal
Generalized Vision-Based Detection, Identification and Pose Estimation of Lamps for BIM Integration
Next Article in Special Issue
Comparison of CO2 Vertical Profiles in the Lower Troposphere between 1.6 µm Differential Absorption Lidar and Aircraft Measurements Over Tsukuba
Previous Article in Journal
Activities of Daily Living Ontology for Ubiquitous Systems: Development and Evaluation
Previous Article in Special Issue
Atmospheric Pollution Monitoring in Urban Area by Employing a 450-nm Lidar System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improvement of CO2-DIAL Signal-to-Noise Ratio Using Lifting Wavelet Transform

1
School of Remote Sensing & Geomatics Engineering, Nanjing University of Information Science and Technology, Nanjing 210044, China
2
School of Remote Sensing and Information Engineering, Wuhan University, Wuhan 430079, China
3
School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China
4
Electronic Information School, Wuhan University, Wuhan 430072, China
5
State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing, Wuhan University, Wuhan 430079, China
*
Author to whom correspondence should be addressed.
Submission received: 25 May 2018 / Revised: 2 July 2018 / Accepted: 18 July 2018 / Published: 20 July 2018

Abstract

:
Atmospheric CO2 plays an important role in controlling climate change and its effect on the carbon cycle. However, detailed information on the dynamics of CO2 vertical mixing remains lacking, which hinders the accurate understanding of certain key features of the carbon cycle. Differential absorption lidar (DIAL) is a promising technology for CO2 detection due to its characteristics of high precision, high time resolution, and high spatial resolution. Ground-based CO2-DIAL can provide the continuous observations of the vertical profile of CO2 concentration, which can be highly significant to gaining deeper insights into the rectification effect of CO2, the ratio of respiration photosynthesis, and the CO2 dome in urban areas. A set of ground-based CO2-DIAL systems were developed by our team and highly accurate long-term laboratory experiments were conducted. Nonetheless, the performance suffered from low signal-to-noise ratio (SNR) in field explorations because of decreasing aerosol concentrations with increasing altitude and surrounding interference according to the results of our experiments in Wuhan and Huainan. The concentration of atmospheric CO2 is derived from the difference of signals between on-line and off-line wavelengths; thus, low SNR will cause the superimposition of the final inversion error. In such a situation, an efficient and accurate denoising algorithm is critical for a ground-based CO2-DIAL system, particularly in field experiments. In this study, a method based on lifting wavelet transform (LWT) for CO2-DIAL signal denoising was proposed. This method, which is an improvement of the traditional wavelet transform, can select different predictive and update functions according to the characteristics of lidar signals, thereby making it suitable for the signal denoising of CO2-DIAL. Experiment analyses were conducted to evaluate the denoising effect of LWT. For comparison, ensemble empirical mode decomposition denoising was also performed on the same lidar signal. In addition, this study calculated the coefficient of variation (CV) at the same altitude among multiple original signals within 10 min and then performed the same calculation on the denoised signal. Finally, high-quality signal of ground-based CO2-DIAL was obtained using the LWT denoising method. The differential absorption optical depths of the denoised signals obtained via LWT were calculated, and the profile distribution information of CO2 concentration was acquired during field detection by using our developed CO2-DIAL systems.

1. Introduction

At present, the occurrence of global warming has been further confirmed by the fifth working report of the International Panel of Climate Change [1]. In response to this global challenge, many countries signed the Paris Agreement at the World Climate Conference on 12 December 2015, thereby setting the long-term goal that global temperature will be controlled to rise within 2 °C compared with the industrialized period and realizing balance between carbon emission and sequestration by 2050. The continued increase of CO2 concentration, as the major positive radiative forcing factor, has exerted a huge impact on global climate. However, the temporal and spatial distributions of global CO2 remain largely uncertain. In particular, the vertical profile information of CO2 concentration is still lacking, which restricts an in-depth understanding of the carbon cycle. The precise monitoring of the vertical distribution of CO2 concentration is currently an important and challenging scientific issue [2,3,4]. As a representative of active detection technology, differential absorption lidar (DIAL) is widely used in probing atmospheric CO2 given that it cannot only measure CO2 column but can also provide vertical profiles of CO2 concentration distribution. To date, multiple sets of CO2-DIAL systems have been established by developed countries, such as the United States, Germany, and Japan [5,6,7,8,9,10].
A set of ground-based CO2-DIAL systems has developed by our team since 2013, and highly precise long-term laboratory experiments have been conducted. Results with a random error of less than 1% were achieved in horizontal detection experiments [11]. However, the performance suffered from low signal-to-noise ratio (SNR) in the vertical detection experiments due to decreasing aerosol concentrations with respect to altitude according to the result of our experiments in Wuhan and Huainan. In 2017, the CO2-DIAL system was transported to Yangbajing, a small town in Tibet with an altitude of 4300 m, and SNR decreased sharply because of the clean air and low atmospheric pressure in the area. Moreover, two beam detection lasers are transmitted for the CO2-DIAL system and the inversion process requires a logarithmic ratio of on-line and off-line wavelength signal intensities. Inversion accuracy depends directly on the denoising quality of on-line and off-line wavelength signals, and low SNR causes the superimposition of the final inversion error. In such cases, an efficient signal denoising algorithm is essential for improving the detection distance and inversion accuracy of the ground-based CO2-DIAL system.
Numerous methods have been improved to reduce the noise level [12,13,14,15,16]. Multiple-pulse accumulation and averaging is one of the most widely used denoising methods. However, the denoising effect of this method is poor and useful information can be easily lost in the signal. Other low-pass filters, including Fourier transform (FT), empirical mode decomposition (EMD), and traditional wavelet transform (TWT), are widely used in signal denoising and have achieved good effects. Nevertheless, these methods focus on the denoising of echo signals without considering the subsequent application of signals. Although the denoising of on-line and off-line wavelengths improves the SNR of lidar signals, CO2-DIAL always loses important information required for differential absorption detection, thereby reducing the accuracy of CO2 concentration inversion. The signals of on-line and off-line wavelengths are received simultaneously by CO2-DIAL; thus, the difference information produced by CO2 absorption between two signals should be retained while improving SNR. Consequently, the aforementioned methods are ineffective for denoising CO2-DIAL signals.
This study introduced a denoising method based on lifting wavelet transform (LWT) for the lidar signals of our ground-based CO2-DIAL system. The proposed denoising method is an improvement of TWT and a spatial filter applied to each signal as a function of range. It selects different predictive and update functions according to signal characteristics. Consequently, the LWT denoising method can preserve the characteristic information of signals while improving SNR, which makes it suitable for the signal denoising of CO2-DIAL. The proposed method plays a significant role in the accurate field detection of CO2.

2. Principle and System

Figure 1 shows the diagram of the ground-based CO2-DIAL [17] system. In the system, two laser pulses with similar wavelengths are transmitted. One laser, called the on-line wavelength (λon), is located at the absorption peak of the detected component to obtain maximum absorption. The other laser, called the offline wavelength (λoff), is located near the valley of absorption, which results in minimum absorption.
The lidar equations of the on-line and off-line wavelengths can be written as
P ( λ , r ) = K P 0 A c τ 2 r 2 β ( λ , r ) exp { 2 0 r [ α 0 ( λ , r ) + N g ( r ) σ g ( λ ) ] d r } ,
where r is the detection range, P(λ,r) is the received power of range r (λ can be on-line and off-line wavelengths), P0 is the laser output power, K is the calibration constant for lidar, A is the light area of the receiving telescope, c is the speed of light, τ is the laser pulse duration, β(λ,r) is the backscatter coefficient of the atmosphere, α0(λ,r) is the extinction coefficient of the atmosphere (excluding the trace gas under study), Ng(r) is the number density of trace gas, and σg(r) is the absorption cross section of trace gas [18,19].
In the CO2-DIAL system, Equation (1) with the on-line wavelength is divided by that with the off-line wavelength. The number density of the range-resolved CO2 can then be derived using the following equation:
N g ( r ) = 1 2 Δ r [ σ g ( λ o n ) σ g ( λ o f f ) ] ln [ P ( λ o f f , r 2 ) P ( λ o n , r 1 ) P ( λ o n , r 2 ) P ( λ o f f , r 1 ) ] ,
where r1 and r2 are the beginning and end of the integration interval, respectively; and Δr = r2r1 is the range resolution.
In addition, the concentration of CO2 (unit: parts per million, ppm) can be obtained by using the following equation with number density:
C = N R T N A P × 10 6 ,
where NA is the Avogadro constant of approximately 6.022 × 10 23 mol 1 ; P and T are the pressure and temperature of gas, respectively; and R is the ideal gas constant of approximately 8.314 J/(K·mol).
Equation (2) can also be expressed as the following equation of optical depth (OD) and differential absorption OD (DAOD):
N g ( r ) = D A O D 2 Δ r Δ σ = O D ( o n ) O D ( o f f ) 2 Δ r [ σ g ( λ o n ) σ g ( λ o f f ) ] .
When the laser wavelength range is small and the time frame is short, several atmospheric parameters change slightly with wavelength, which can be regarded as constant in a DIAL system.
The configuration of our ground-based CO2-DIAL system is shown in Figure 2. It is composed of a laser launch system, an optical reception detection system, and a data acquisition and processing system. A pulsed dye laser is used for wavelength modulation in our CO2-DIAL system. The target lasers with a wavelength of 1572 nm and a laser repetition frequency of 20 Hz are generated via difference frequency mixing between the fundamentals of the Nd:YAG laser (1064 nm) and the dye laser (~634 nm), which is pumped by the second harmonic of the Nd:YAG laser (532 nm). Consequently, the wavelength of the output laser can be changed by tuning the red laser, the wavelength of which is changed using a stepper motor. In addition, 8% of the pulsed laser light transmitted from the laser system is used as the laser source of the wavelength control unit, whereas the rest is vertically launched into the air. The optical reception detection system is composed of a receiving telescope with an aperture of 1000 mm, a filter bank, and near-infrared photomultiplier tube for echo signal reception, background noise removal, and photoelectric conversion. The data acquisition and processing system is based on a Licel transient recorder for realizing data acquisition and processing functions [11,20,21,22]. The resolution and the sampling frequency of the digitizer is 12-Bit and 20 MHz, respectively. Figure 3 presents the physical map of our ground-based CO2-DIAL system, and the exterior view of the Yangbajing site on the Tibetan Plateau is shown in Figure 4.
On the basis of the lidar equations, the quality of a CO2 signal will seriously affect the inversion results, particularly for our ground-based CO2-DIAL system. In a CO2-DIAL system, the concentration of atmospheric CO2 is derived from the difference of signals between on-line and off-line wavelengths, and thus, inversion accuracy depends directly on denoising quality. The high-accuracy detection of CO2 concentration does not only require an echo signal with high SNR but also the maximum retention of the characteristic information of lidar signals. To improve the detection accuracy of CO2-DIAL, a flexible and accurate denoising method based on LWT for the lidar signals of our ground-based CO2-DIAL system was proposed.

3. Lidar Signal Denoising Based on LWT

The concept of wavelet denoising is derived from stretching and translation, which are conducted on a wavelet basis. First, a signal is transformed into a multiscale signal, and then the wavelet coefficients of the effective signal are extracted at each scale as much as possible while removing the wavelet coefficients that belong to the noise signal. Finally, the signal is reconstructed via inverse wavelet transform. Although TWT is widely used in signal denoising, it requires selecting an appropriate wavelet basis to adapt to different signals. Once a wavelet basis is selected, this basis can no longer be changed during the noise reduction process. This nonadaptive selection of a wavelet basis is a difficult point in TWT application [23].
Although the LWT proposed in this study is based on TWT, it differs from TWT. LWT does not rely on FT, and all its operations are conducted in the time domain. Its wavelet basis is no longer fixed because of the adoption of two independent filters (i.e., predictors and updators) in the lifting wavelet scheme. Each layer decomposition can select the appropriate wavelet basis according to signal characteristics. Furthermore, each step is completely reversible in LWT, thereby providing it with strong adaptability. Consequently, the LWT denoising method is suitable for changeable lidar signals and can effectively preserve the characteristic information of signals [24,25,26].
In the current study, LWT is applied to the denoising of lidar signals based on the signal and noise characteristics of a CO2-DIAL system. The basic principle of LWT is to divide a wavelet transform into several steps, with each step conducted according to the characteristics of the signal for an appropriate transformation process. The main steps of lidar signal denoising via LWT are as follows [27,28,29,30,31].

3.1. Lifting Wavelet Decomposition (Split–Predict–Update)

(1) Split
This step splits data into two sequences according to odd and even numbers. The two sequences are not intersecting, i.e.,
S = ( S o ( n ) , S e ( n ) ) , S o ( n ) = x ( 2 k 1 ) , S e ( n ) = x ( 2 k ) , k z ,
where So(n) and Se(n) correspond to the odd and even sequence signals, respectively.
(2) Predict
P(Se(n)) is obtained by selecting an appropriate filter function (marked as P) and applying it to the even signal. So(n) can be replaced by P(Se(n)) because Se(n) and So(n) are parity adjacent, and a change between them should be minimal in case of a smooth signal. The advantage of prediction considerably reduces the number of signals for processing. The difference between So(n) and P(Se(n)) [marked as d(n)] also represents the approximation degree between the parity signals. In this case, d(n) can be regarded as the high-frequency (HF) component of the signal.
d ( n ) = S o ( n ) P ( S e ( n ) )
(3) Update
The main purpose of this step is to enable the predicted Se(n) to maintain the same global features as the original signal and to obtain the low-frequency (LF) components of the original signal. In the update process, a function U is introduced to modify d(n), and then the updated d(n) is used to update Se(n) to obtain a new subset of S(K).
S ( k ) = S e ( n ) + U ( d ( n ) )
A layer decomposition of S(n) is achieved through three steps. Furthermore, the same processes of split, predict, and update are performed on the subset. After j times lifting wavelet decomposition, S(n) can be expressed as [Snj, dj, dj−1, , d1], where Snj represents the LF of the signal, i.e., the useful information, whereas [dj, , d1] represents the HF of the signal. The principle diagram of signal decomposition is shown in Figure 5. With the assumption that P = [p1, p2] and U = [u1, u2, u3, u4], the decomposition process of a lifting wavelet is shown in Figure 6.

3.2. Threshold Processing

After wavelet decomposition, an appropriate threshold processing method is necessary to deal with the wavelet coefficients obtained through decomposition. At present, hard-field and soft-field value processing are the primary threshold processing methods. In threshold processing, the wavelet coefficients are compared with a given threshold. When the absolute value of the wavelet coefficient is less than that of the threshold, the two processing methods make all the wavelet coefficients zero [32,33,34,35]. However, when the absolute value of the wavelet coefficients is higher than the threshold, the hard-field threshold processing method keeps the coefficients unchanged, whereas the soft-field threshold processing method regards the difference between the wavelet coefficients and the threshold as new wavelet coefficients. In this work, the soft-field threshold processing method is adopted for lidar signal denoising due to its better smoothing effect.
d τ = sgn ( d ( n ) ) ( | d ( n ) τ | ) = { 0 , | x | τ d ( n ) τ , | x | > τ ,
τ = σ 2 ln N , σ = 1 0.6745 M e d ( | d | ) ,
where τ is the threshold, σ is the standard deviation estimate of the noise, N is the sampling number of the detail signal, and Med() is the median function.

3.3. Lifting Wavelet Reconstruction (Resume Update–Resume Prediction–Merge)

The last step of denoising is the reconstruction of the wavelet coefficients that were processed by the threshold. In LWT denoising, signal reconstruction is the opposite of signal decomposition. Therefore, a reconstructed signal can be achieved by simply replacing the calculation order and the sign, thereby preventing the failure of signal reconstruction. This feature is another advantage of LWT denoising. The principle diagram of signal reconstruction is shown in Figure 7, and the main processes are described as follows.
(1) Resume update
The even signal Se(n) is resumed from the LF signal S(k) and the HF signal d(n).
S e ( n ) = S ( k ) U ( d ( n ) )
(2) Resume prediction
The odd signal So(n) is resumed by the even signal Se(n) and the HF signal d(n).
S o ( n ) = P ( S e ( n ) ) + d ( n )
(3) Merge
The original signal is resumed by the odd signal So(n) and the even signal Se(n).
S = ( S o ( n ) , S e ( n ) )
The lifting wavelet reconstruction process is illustrated in Figure 8.
Finally, the flowchart of lidar signal denoising via LWT for atmospheric CO2 detection is shown in Figure 9.

4. Experiment Analysis

LWT is applied to the processing of CO2-DIAL signals to evaluate the effect of denoising. The lidar signals, which are collected from 8:00 p.m. to 11:30 p.m. on 30 July 2016 in Huainan, China, are used as the experiment data. The time resolution of our CO2-DIAL system is 30 s, the laser transmission frequency is 20 Hz, the resolution and sampling frequency of Licel is 12-Bit and 20 MHz, respectively. The wavelengths of on-line and off-line lasers in our system are 1572.018 nm and 1572.150 nm, respectively. For simplicity, this paper presents only off-line wavelength signals because the echo signals of on-line and off-line wavelengths are similar.

4.1. Background Noise Deduction

The background noises of lidar signals are mainly composed of two parts: system background noise and sky background light noise. These noises can be measured and removed roughly in a lidar detection system. In our experiments, the process of system background noise removal are as follows: The CO2-DIAL system was first operating normally for 5 min with the telescope covered to ensure that there are no scattering signals received by the infrared detector. Then, lidar signals were collected and used to estimate a system noise for subtraction from the subsequent live signal. As to the sky background noise, the signals received by the telescope at high detection area can be considered as the sky background light noise because the laser detector cannot reach this height. Therefore, in our experiment, we calculated the mean value of signal intensity of about 20–22 km altitude range as sky background noise and subtracted them from lidar signals. Figure 10 illustrates the lidar signal after the subtraction of background noise, whereas Figure 11 shows the signal corrected by distance.
It is undeniable that taking these measures cannot remove all the noise, and the random noise still exists. In order to obtain high accuracy vertical profile information of CO2 concentration, an effective signal denoising method is still needed to improve the quality of the signal in the subsequent data processing.

4.2. Selection of Wavelet Function and Decomposition Scale

Before denoising, the wavelet function and the decomposition scale are first selected. This selection directly determines the quality of denoising. A comparison of various wavelet functions with lidar signals shows that the db5 wavelet exhibits a high degree of similarity with lidar signals. Figure 12 shows the graph of the db5 wavelet scaling function and the wavelet function. Consequently, this study adopts the db5 wavelet lifting algorithm for lidar signal denoising.
A simulation experiment was conducted to select a suitable wavelet decomposition scale for signal denoising via LWT, which can enable the denoised signal obtained through LWT to preserve the characteristic information of signals while improving SNR. In the experiment, the concentration of CO2 was assumed as a fixed value, and then the theoretical lidar signals were obtained by the inversion equation, i.e., Equation (2). The db5 wavelet was used to denoise analog signals at different decomposition levels, and the decomposition scales are one to four layers in turn. After signal denoising is completed, the linear relationship between the original and denoised signals was calculated to evaluate the denoising effect.
The results are provided in Figure 13, where the abscissa represents the original signal and the ordinate represents the denoised signal. A good linear relationship exists between the two signals, thereby indicating that the lifting wavelet algorithm exerts a good denoising effect on lidar signals and can maintain the characteristic of the original signal. In addition, the best denoising effect occurs when the decomposition level has three layers according to the linear fitting parameters after the denoising of different decomposition levels. Such effect can effectively preserve the characteristic information of signals.
Table 1 lists the linear fitting parameters between the original signal and the denoised signal under different wavelet decomposition scales. The table shows that when the decomposition scale has three layers, the degree of linear fitting is the highest, which indicates the best denoising effect. After the decomposition of layer three, the denoising effect deteriorates mainly because the signal itself is smoothed away with considerably more useful information as the decomposition scale increases. Therefore, the db5 wavelet function decomposition scale used in the lifting wavelet denoising process in this work has three levels.

4.3. Lifting Wavelet Signal Denoising

In this study, lidar signals with an altitude range of 1000–3000 m are selected to conduct wavelet denoising. Figure 14 shows the lidar signals in this region, which is above the atmospheric boundary layer and has a considerably weakened signal intensity, thereby causing low SNR. Therefore, lidar signal denoising is important in this region. On the basis of a pre-period analysis and the selection of a lifting wavelet function and a decomposition level, the mother wavelet of the denoising algorithm in this study is a db5 wavelet, and the decomposition scale has three layers. The denoising effects of the online and off-line wavelength signals are shown in Figure 15.
The figure shows that compared with the original signal, the signal denoised using the lifting wavelet algorithm is considerably smoother while still retaining the basic characteristics of the original signal.

4.4. Evaluation of Denoising Effect

For the measured signal, the relative error among echo signals is generally used to evaluate signal quality. The time interval between each signal is extremely short and CO2 concentration does not change considerably within such a short time. Thus, the relative error among the signal intensities of echo signals in the adjacent time slot should be minimal at the same detection altitude. However, the relative error increases due to the influence of noise. Therefore, signal quality can be evaluated by calculating the relative error of multiple echo signals at the same altitude over a certain period, i.e., the coefficient of variation (CV).
CV is calculated as follows:
C V = σ / μ ,
where σ is the standard deviation of the echo signal intensity at the same altitude, and u is the mean value of data.
To evaluate the denoising effect of lifting wavelet, this study calculated the CV at the same altitude among multiple original signals within 10 min and then performed the same calculation on the denoised signal.
The CVs at different altitudes are shown in Figure 16. The CVs of the on-line and off-line wavelengths are considerably reduced after denoising using the lifting wavelet algorithm, which indicates that the algorithm exhibits good capability for lidar signal denoising. The figure also illustrates that SNR decreases while CV increases at a high altitude.
To prove the signal denoising effect of lifting wavelet, the ensemble empirical mode decomposition (EEMD) denoising is conducted on the same lidar signal in this study. EEMD, which is an improvement of EMD, was proposed by Huang et al. to solve the modal aliasing problem of EMD. It decomposes a signal into a series of intrinsic mode functions, which are arranged from high to low frequencies at different time scales. Noise generally exists in high-frequency components; thus, signal noise reduction can be achieved by removing high-frequency components [13,36,37]. EEMD does not require pre-selection of the basis function. This method decomposes based on the signal itself, which makes it adaptable; it also exhibits good denoising effect on lidar signals [38,39]. The denoising effect of the lifting wavelet method proposed in this work can be measured well by comparing with the result of EEMD denoising. The results are presented in Figure 17.
Figure 18 shows the CVs of the intensity of multiple signals at the same height. The average values of CV calculated from the original signal and the denoised signals using the two methods are listed in Table 2. The lifting wavelet algorithm proposed in this study is demonstrated to be capable of achieving better denoising effect than EEMD.
In addition, the DAOD of lidar signal after denoising via LWT is calculated. The results are presented in Figure 19, and the linear fitting parameters are listed in Table 3.
The DAOD singular value of the signal caused by noise is eliminated after denoising via LWT, and the linear fitting effect is also considerably improved. These results show that the LWT algorithm proposed in this study can effectively improve the quality of echo signals. The LWT denoising method cannot only improve SNR but can also preserve the characteristic information of signals, which is favorable for the subsequent detection of CO2 using our CO2-DIAL system.
Finally, long-term field experiments were conducted in Huainan, China, and the profile distribution information of CO2 concentration was obtained using our CO2-DIAL system, as shown in Figure 20.

5. Conclusions

In this study, a denoising method based on LWT was proposed to acquire high-quality echo signals of our ground-based CO2-DIAL system in field exploration. The LWT denoising method can preserve the characteristic information of signals while improving SNR, which is highly significant for the accurate field detection of CO2. Experiment analyses were conducted to evaluate the denoising effect of LWT. For comparison, EEMD denoising was also conducted on the same lidar signal. Moreover, this study calculated the CVs at the same altitude among multiple original signals within 10 min, and then performed the same calculation on the denoised signal. The respective average values of CV calculated from the denoised signals of LWT and EEMD were 0.1637 and 0.1751 for on-line wavelength and 0.1508 and 0.1645 for off-line wavelength, respectively. The results show that the CVs of multiple signals at the same height were reduced after denoising via LWT, and the average CV value calculated from the denoised signals of LWT was smaller than that calculated from the denoised signals of EEMD, thereby indicating that LWT achieved better denoising effect than EEMD. In addition, the DAOD values of the original and denoised signals via LWT were calculated. The DAOD singular value of the signal caused by noise was eliminated after denoising by LWT, and the linear fitting effect was considerably improved, with R2 increased from 0.619 to 0.884, thereby implying that the detection accuracy of CO2 was significantly improved. Finally, long-term field experiments were conducted in Huainan, China, and the profile distribution information of CO2 concentration was obtained using our CO2-DIAL system with the LWT denoising method.

Author Contributions

Conceptualization, C.X. and W.G.; Methodology, C.X.; Software, X.M. and Y.Z; Validation, C.X. and G.H.; Formal Analysis, C.X. and G.H.; Investigation, G.H.; Resources, C.X. and G.H.; Data Curation, C.X.; Writing—Original Draft Preparation, C.X.; Writing—Review & Editing, C.X. and G.H.; Visualization, C.X.; Supervision, G.H.; Project Administration, X.M.; Funding Acquisition, C.X., G.H.,X.M. and W.G.

Funding

This work was supported by the National Key R&D Program of China (Grant No. 2017YFC0212600), the National Natural Science Foundation of China (Grant Nos. 41601351, 41627804 and U1531122), the China Postdoctoral Science Foundation funded project (Grant No. 2017T100580), and the Talent Launch Fund of Nanjing University of Information Science and Technology (Grant No. 2017r066).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Barros, V.; Field, C.; Dokke, D.; Mastrandrea, M.; Mach, K.; Bilir, T.; Chatterjee, M.; Ebi, K.; Estrada, Y.; Genova, R. Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part B: Regional Aspects. In Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; IPCC: Stanford, CA, USA, 2015. [Google Scholar]
  2. Tans, P.P.; Fung, I.Y.; Takahashi, T. Observational contrains on the global atmospheric CO2 budget. Science 1990, 247, 1431–1438. [Google Scholar] [CrossRef] [PubMed]
  3. Stephens, B.B.; Gurney, K.R.; Tans, P.P.; Sweeney, C.; Peters, W.; Bruhwiler, L.; Ciais, P.; Ramonet, M.; Bousquet, P.; Nakazawa, T. Weak northern and strong tropical land carbon uptake from vertical profiles of atmospheric CO2. Science 2007, 316, 1732–1735. [Google Scholar] [CrossRef] [PubMed]
  4. Gatti, L.V.; Miller, J.B.; D’amelio, M.T.S.; Martinewski, A.; Basso, L.S. Vertical profiles of CO2 above eastern Amazonia suggest a net carbon flux to the atmosphere and balanced biosphere between 2000 and 2009. Tellus B 2010, 62, 581–594. [Google Scholar] [CrossRef]
  5. Gibert, F.; Edouart, D.; Cénac, C.; Le Mounier, F.; Dumas, A. 2-μm Ho emitter-based coherent DIAL for CO2 profiling in the atmosphere. Opt. Lett. 2015, 40, 3093–3096. [Google Scholar] [CrossRef] [PubMed]
  6. Ismail, S.; Koch, G.; Abedin, N.; Refaat, T.; Rubio, M.; Davis, K.; Miller, C.; Vay, S.; Singh, U. Development of a ground-based 2-micron differential absorption lidar system to profile tropospheric CO2. In Proceedings of the NASA Earth Science Technology Conference 2006, College Park, MD, USA, 27–29 June 2006. [Google Scholar]
  7. Koch, G.J.; Barnes, B.W.; Petros, M.; Beyon, J.Y.; Amzajerdian, F.; Yu, J.; Davis, R.E.; Ismail, S.; Vay, S.; Kavaya, M.J.; et al. Coherent differential absorption lidar measurements of CO2. Appl. Opt. 2004, 43, 5092–5099. [Google Scholar] [CrossRef] [PubMed]
  8. Fiorani, L.; Santoro, S.; Parracino, S.; Nuvoli, M.; Minopoli, C.; Aiuppa, A. Volcanic CO2 detection with a DFM/OPA-based lidar. Opt. Lett. 2015, 40, 1034–1036. [Google Scholar] [CrossRef] [PubMed]
  9. Queiβer, M.; Granieri, D.; Burton, M. A new frontier in CO2 flux measurements using a highly portable DIAL laser system. Sci. Rep. 2016, 6. [Google Scholar] [CrossRef] [PubMed]
  10. Shibata, Y.; Nagasawa, C.; Abo, M. Development of 1.6 μm DIAL using an OPG/OPA transmitter for measuring atmospheric CO2 concentration profiles. Appl. Opt. 2017, 56, 1194–1201. [Google Scholar] [CrossRef] [PubMed]
  11. Han, G.; Gong, W.; Ma, X.; Xiang, C.Z.; Liang, A.L.; Zheng, Y.X. A ground-based differential absorption lidar for atmospheric vertical CO2 profiling. Acta Phys. Sin. 2015, 64, 0244206. [Google Scholar]
  12. Milton, M.; Woods, P.; Jolliffe, B.; Swann, N.; McIlveen, T. Measurements of toluene and other aromatic hydrocarbons by differential-absorption LIDAR in the near-ultraviolet. Appl. Phys. B 1992, 55, 41–45. [Google Scholar] [CrossRef]
  13. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.; Yen, N.-C.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. Lond. A 1998, 454, 903–995. [Google Scholar] [CrossRef]
  14. Zhou, Z.; Hua, D.; Wang, Y.; Yan, Q.; Li, S.; Li, Y.; Wang, H. Improvement of the signal to noise ratio of Lidar echo signal based on wavelet de-noising technique. Opt. Lasers Eng. 2013, 51, 961–966. [Google Scholar] [CrossRef]
  15. Mao, F.; Gong, W.; Li, C. Anti-noise algorithm of lidar data retrieval by combining the ensemble Kalman filter and the Fernald method. Opt. Express 2013, 21, 8286–8297. [Google Scholar] [CrossRef] [PubMed]
  16. Wu, S.; Liu, Z.; Liu, B. Enhancement of lidar backscatters signal-to-noise ratio using empirical mode decomposition method. Opt. Commun. 2006, 267, 137–144. [Google Scholar] [CrossRef]
  17. Xiang, C.; Ma, X.; Liang, A.; Han, G.; Gong, W.; Yan, F. Feasibility Study of Multi-Wavelength Differential Absorption LIDAR for CO2 Monitoring. Atmosphere 2016, 7, 89. [Google Scholar] [CrossRef]
  18. Gong, W.; Liang, A.; Han, G.; Ma, X.; Xiang, C. Sensitivity of on-line wavelength during retrieval of atmospheric CO2 vertical profile. Photonics Res. 2015, 3, 146–152. [Google Scholar] [CrossRef]
  19. Chantjaroen, C.; Repasky, K.S. Development of a Singly-Resonant Optical Parametric Oscillation for a Laser Transmitter of an Eye-Safe Ground-Based Differential Absorption Lidar (DIAL) for Carbon Dioxide (CO2) and Methane (CH4) Studies. In Proceedings of the AGU Fall Meeting Abstracts, San Francisco, CA, USA, 12–16 December 2016. [Google Scholar]
  20. Qi, Z.; Zhang, T.; Han, G.; Li, D.; Ma, X.; Gong, W. A nonlinear merging method of analog and photon signals for CO2 detection in lower altitudes using differential absorption lidar. Opt. Commun. 2017, 388, 68–76. [Google Scholar] [CrossRef]
  21. Gong, W.; Ma, X.; Han, G.; Xiang, C.; Liang, A.; Fu, W. Method for wavelength stabilization of pulsed difference frequency laser at 1572 nm for CO2 detection lidar. Opt. Express 2015, 23, 6151–6170. [Google Scholar] [CrossRef] [PubMed]
  22. Han, G.; Gong, W.; Lin, H.; Ma, X.; Xiang, C. On-line wavelength calibration of pulsed laser for CO2 DIAL sensing. Appl. Phys. B 2014, 117, 1041–1053. [Google Scholar] [CrossRef]
  23. Fang, H.-T.; Huang, D.-S.; Wu, Y.-H. Antinoise approximation of the lidar signal with wavelet neural networks. Appl. Opt. 2005, 44, 1077–1083. [Google Scholar] [CrossRef] [PubMed]
  24. Ghael, S.; Sayeed, A.M.; Baraniuk, R.G. Improved wavelet denoising via empirical Wiener filtering. Proc. SPIE 1997, 3169, 389–399. [Google Scholar] [Green Version]
  25. Sweldens, W. The lifting scheme: A construction of second generation wavelets. SIAM J. Math. Anal. 1998, 29, 511–546. [Google Scholar] [CrossRef]
  26. Jansen, M.H.; Oonincx, P.J. Second Generation Wavelets and Applications; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2005. [Google Scholar]
  27. Ercelebi, E. Electrocardiogram signals de-noising using lifting-based discrete wavelet transform. Comput. Biol. Med. 2004, 34, 479–493. [Google Scholar] [CrossRef]
  28. Daubechies, I.; Sweldens, W. Factoring wavelet transforms into lifting steps. J. Fourier Anal. Appl. 1998, 4, 247–269. [Google Scholar] [CrossRef]
  29. Sree, P.S.J.; Kumar, P.; Siddavatam, R.; Verma, R. Salt-and-pepper noise removal by adaptive median-based lifting filter using second-generation wavelets. Signal Image Video Proc. 2013, 7, 1–8. [Google Scholar] [CrossRef]
  30. Szczesna, A. Designing lifting scheme for second generation wavelet-based multiresolution processing of irregular surface meshes. Proc. Comput. Graph. Vis. 2008, 103–111. [Google Scholar]
  31. Quellec, G.; Lamard, M.; Cazuguel, G.; Cochener, B.; Roux, C. Adaptive nonseparable wavelet transform via lifting and its application to content-based image retrieval. IEEE Trans. Image Proc. 2010, 19, 25–35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Donoho, D.L.; Johnstone, I.M. Threshold selection for wavelet shrinkage of noisy data. In Proceedings of the 16th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Baltimore, MD, USA, 3–6 November 1994; Volume 1, pp. A24–A25. [Google Scholar]
  33. Kaur, L.; Gupta, S.; Chauhan, R.C. Image Denoising Using Wavelet Thresholding. ICVGIP 2002, 2, 16–18. [Google Scholar]
  34. Han, M.; Liu, Y.; Xi, J.; Guo, W. Noise smoothing for nonlinear time series using wavelet soft threshold. IEEE Signal Process. Lett. 2007, 14, 62–65. [Google Scholar] [CrossRef]
  35. Berkner, K.; Wells, R.O. Smoothness estimates for soft-threshold denoising via translation-invariant wavelet transforms. Appl. Comput. Harmon. Anal. 2002, 12, 1–24. [Google Scholar] [CrossRef]
  36. Barnhart, B.L. The Hilbert-Huang Transform: Theory, Applications, Development; The University of Iowa: Iowa City, IA, USA, 2011. [Google Scholar]
  37. Wu, Z.; Huang, N.E. Ensemble empirical mode decomposition: A noise-assisted data analysis method. Adv. Adapt. Data Anal. 2009, 1, 1–41. [Google Scholar] [CrossRef]
  38. Li, J.; Gong, W.; Ma, Y. Atmospheric lidar noise reduction based on ensemble empirical mode decomposition. In Proceedings of the International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Melbourne, Australia, 5 August–1 September 2012; Volume XXXIX-B8, pp. 127–129. [Google Scholar]
  39. He, J.; Liu, W.; Zhang, Y. New method of lidar ceilometer backscatter signal processing based on Hilbert-Huang transform. Infrared Laser Eng. 2012, 41, 397–403. [Google Scholar]
Figure 1. Diagram of ground-based CO2-DIAL system.
Figure 1. Diagram of ground-based CO2-DIAL system.
Sensors 18 02362 g001
Figure 2. Configuration of our ground-based CO2-DIAL system.
Figure 2. Configuration of our ground-based CO2-DIAL system.
Sensors 18 02362 g002
Figure 3. Physical map of our ground-based CO2-DIAL system.
Figure 3. Physical map of our ground-based CO2-DIAL system.
Sensors 18 02362 g003
Figure 4. Exterior view of the Yangbajing probing site.
Figure 4. Exterior view of the Yangbajing probing site.
Sensors 18 02362 g004
Figure 5. Diagram of lifting wavelet decomposition.
Figure 5. Diagram of lifting wavelet decomposition.
Sensors 18 02362 g005
Figure 6. Rendering of lifting wavelet decomposition.
Figure 6. Rendering of lifting wavelet decomposition.
Sensors 18 02362 g006
Figure 7. Diagram of lifting wavelet reconstruction.
Figure 7. Diagram of lifting wavelet reconstruction.
Sensors 18 02362 g007
Figure 8. Rendering of lifting wavelet reconstruction.
Figure 8. Rendering of lifting wavelet reconstruction.
Sensors 18 02362 g008
Figure 9. Flowchart of lidar signal denoising via LWT.
Figure 9. Flowchart of lidar signal denoising via LWT.
Sensors 18 02362 g009
Figure 10. Lidar signal after the subtraction of background noise (30 July 2016).
Figure 10. Lidar signal after the subtraction of background noise (30 July 2016).
Sensors 18 02362 g010
Figure 11. Lidar signal corrected by distance (30 July 2016).
Figure 11. Lidar signal corrected by distance (30 July 2016).
Sensors 18 02362 g011
Figure 12. Graph of the db5 wavelet scaling function and wavelet function.
Figure 12. Graph of the db5 wavelet scaling function and wavelet function.
Sensors 18 02362 g012
Figure 13. Denoising effect of db5 LWT under different decomposition scales (the decomposition scales are one to four layers in turn).
Figure 13. Denoising effect of db5 LWT under different decomposition scales (the decomposition scales are one to four layers in turn).
Sensors 18 02362 g013aSensors 18 02362 g013b
Figure 14. Range-corrected signal at an altitude range of 1000–3000 m (30 July 2016).
Figure 14. Range-corrected signal at an altitude range of 1000–3000 m (30 July 2016).
Sensors 18 02362 g014
Figure 15. Original and denoised signals via LWT ((left) on-line wavelength; (right) off-line wavelength).
Figure 15. Original and denoised signals via LWT ((left) on-line wavelength; (right) off-line wavelength).
Sensors 18 02362 g015
Figure 16. CVs of original signals and denoised signals via LWT ((left) on-line wavelength; (right) off-line wavelength).
Figure 16. CVs of original signals and denoised signals via LWT ((left) on-line wavelength; (right) off-line wavelength).
Sensors 18 02362 g016
Figure 17. Original and denoised signals using EEMD ((left) online wavelength; (right) off-line wavelength).
Figure 17. Original and denoised signals using EEMD ((left) online wavelength; (right) off-line wavelength).
Sensors 18 02362 g017
Figure 18. CVs of the original signals and denoised signals using EEMD ((left) on-line wavelength; (right) off-line wavelength).
Figure 18. CVs of the original signals and denoised signals using EEMD ((left) on-line wavelength; (right) off-line wavelength).
Sensors 18 02362 g018
Figure 19. Calculated DAOD of the signal after denoising via LWT. The above is the DAOD of original signal, and the below is the DAOD of denoised signal via LWT.
Figure 19. Calculated DAOD of the signal after denoising via LWT. The above is the DAOD of original signal, and the below is the DAOD of denoised signal via LWT.
Sensors 18 02362 g019
Figure 20. Pseudo-color map of CO2 vertical concentration detected using our CO2-DIAL system.
Figure 20. Pseudo-color map of CO2 vertical concentration detected using our CO2-DIAL system.
Sensors 18 02362 g020
Table 1. Denoising evaluation parameters of db5 LWT under different decomposition scales.
Table 1. Denoising evaluation parameters of db5 LWT under different decomposition scales.
n = 1n = 2n = 3n = 4
k0.9750.9870.9910.931
R20.8900.9350.9400.854
Table 2. Average values of CV calculated by the original and denoised signals.
Table 2. Average values of CV calculated by the original and denoised signals.
Original SignalEEMD Denoised SignalLWT Denoised Signal
On-line0.21150.17510.1637
Off-line0.19400.16450.1508
Table 3. Linear fitting parameters of the original and denoised signals.
Table 3. Linear fitting parameters of the original and denoised signals.
Original SignalLWT Denoised Signal
k8.286E-58.529 E-5
R20.6190.884

Share and Cite

MDPI and ACS Style

Xiang, C.; Han, G.; Zheng, Y.; Ma, X.; Gong, W. Improvement of CO2-DIAL Signal-to-Noise Ratio Using Lifting Wavelet Transform. Sensors 2018, 18, 2362. https://0-doi-org.brum.beds.ac.uk/10.3390/s18072362

AMA Style

Xiang C, Han G, Zheng Y, Ma X, Gong W. Improvement of CO2-DIAL Signal-to-Noise Ratio Using Lifting Wavelet Transform. Sensors. 2018; 18(7):2362. https://0-doi-org.brum.beds.ac.uk/10.3390/s18072362

Chicago/Turabian Style

Xiang, Chengzhi, Ge Han, Yuxin Zheng, Xin Ma, and Wei Gong. 2018. "Improvement of CO2-DIAL Signal-to-Noise Ratio Using Lifting Wavelet Transform" Sensors 18, no. 7: 2362. https://0-doi-org.brum.beds.ac.uk/10.3390/s18072362

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