Time–Frequency Attribute Analysis of Channel 1 Data of Lunar Penetrating Radar

: The Lunar Penetrating Radar (LPR) carried by the Chang’E-3 (CE-3) and Chang’E-4 (CE-4) mission plays a very important role in lunar exploration. The dual-frequency radar on the rover (DFR) provides a meaningful opportunity to detect the underground structure of the CE-3 landing site. The low-frequency channel (channel 1) maps the underground structure to a depth of several hundred meters, while the high-frequency channel (channel 2) can observe the stratigraphic structure of gravel near the surface. As the low-frequency radar image is troubled by unknown noise, time–frequency analysis of a single trace is applied. Then, a method named complete ensemble empirical mode decomposition (CEEMD) is conducted to decompose the channel 1 data, and the Hilbert transform gives us the chance for further data analysis. Finally, combined with regional geology, previous studies, and channel 2 data, a usability analysis of LPR channel 1 data provides a reference for the availability of the CE-4 LPR data.


Introduction
Exploration of the lunar subsurface structure is a powerful force for the analysis of the stratigraphic thickness, time, and contact relationships, which exert a vital role not only on theorizing the evolution of the subsurface structure, but also on estimating the total resources of the lunar subsurface.
In order to detect the shallow geological structure of the moon, spaceborne radar and synthetic aperture radar (SAR) are mainly used. As early as 1972, Apollo 17 carried out on-board radar detection [1]. It operated at 5, 15, and 150 MHz. In addition, it mainly detected the underground geological structure of the moon. Since 20 November 2007, the Lunar-Radar Sounder (LRS) on Japan's Kaguya probes (SELENE) has been surveying the moon underground. The LRS uses the high-frequency (HF) band (5 MHz) to enable the acquisition of subsurface data at depths of several kilometers [2]. The small SAR was launched on the moon in September 2007 and it carried out the Indian Chandrayaan mission. It mapped the scattering characteristics of the moon's pole and determined the existence and extent of polar ice [3,4]. The Mini-RF instrument on the Lunar Reconnaissance Orbiter (LRO) is a SAR, and it works with bistatic radar from the Arecibo Observatory to conduct data analysis of the region. The results show that there may be near-surface sediments of water ice [4,5].
Unlike the spaceborne radar described above, the lunar radar on Chang'E-3 (CE-3) is a field instrument for underground structure detection [6]. The dual-frequency radar (DFR) offers a great The DFR on the YR offers a meaningful chance to map the underground structure about several hundred meters deep from the CH-1 (60 MHz) and its stratigraphic structure near the surface from the CH-2 (500 MHz). For the underground structure, the CH-1 data are selected. According to the recording parameters, real conditions, and data condition, a set of LPR data pretreatment procedures (Table 1) was projected. The data can be preprocessed to access the CH-1 image ( Figure 2).

Processing
Interpretative Statement i Data reading The data is decomposed into 9 parts. According to their storage format (*.psd, which is the standard storage format for aerospace and aerospace), data and location information will be read one by one.
ii Data registration The data is split into 9 parts and ought to be stitched together. iii

Time lag adjustment
The arrival time of the radar echo is delayed by 69.664 ns, which corresponds to the start time of the recorded data. iv Deep useless data deleting As the SNRs of Part 1 and Part 2 are low and the rest of the data are not collected, the last 10,000 ns of data without the study value needs to be deleted.
v Low SNR traces deleting Some parts with low SNR should be removed. The DFR on the YR offers a meaningful chance to map the underground structure about several hundred meters deep from the CH-1 (60 MHz) and its stratigraphic structure near the surface from the CH-2 (500 MHz). For the underground structure, the CH-1 data are selected. According to the recording parameters, real conditions, and data condition, a set of LPR data pretreatment procedures (Table 1) was projected. The data can be preprocessed to access the CH-1 image ( Figure 2).

Processing
Interpretative Statement i Data reading The data is decomposed into 9 parts. According to their storage format (*.psd, which is the standard storage format for aerospace and aerospace), data and location information will be read one by one.
ii Data registration The data is split into 9 parts and ought to be stitched together. iii

Time lag adjustment
The arrival time of the radar echo is delayed by 69.664 ns, which corresponds to the start time of the recorded data. iv Deep useless data deleting As the SNRs of Part 1 and Part 2 are low and the rest of the data are not collected, the last 10,000 ns of data without the study value needs to be deleted.
v Low SNR traces deleting Some parts with low SNR should be removed.
vi Traces selecting Rover patrols are uneven because the rover may stop at a waypoint to collect other scientific data, for example, APXS data, VNIS data, and so forth. However, LPR never stops collecting data, which results in multiple acquisitions of multiple traces at the same place. The repeating traces are stacked and averaged. Two sets of data are obtained. One is all data without a copy track, and the other is a repeating track that has been stacked and averaged at the waypoint.
vii Location Location information is added to the image.
Appl. Sci. 2020, 10, 535 4 of 17 vi Traces selecting Rover patrols are uneven because the rover may stop at a waypoint to collect other scientific data, for example, APXS data, VNIS data, and so forth. However, LPR never stops collecting data, which results in multiple acquisitions of multiple traces at the same place. The repeating traces are stacked and averaged. Two sets of data are obtained. One is all data without a copy track, and the other is a repeating track that has been stacked and averaged at the waypoint. vii Location Location information is added to the image. There are many kinds of space electromagnetic waves on the moon; the terrain of the moon is very complex under the adverse circumstances, so the signal-to-noise ratio (SNR) of LPR data is low. The repeated gathers of each stop point are superimposed, the purpose of which is to improve the SNR of the single trace at the waypoint [24], and provide more accurate and reasonable onedimensional data for the subsequent analysis of the frequency-domain (FD) and the time-frequency domain (TFD).
In this paper, several analysis and processing methods are utilized, for example, the Fourier transform, short-time Fourier transform, CEEMD, and Hilbert transform. The Fourier transform (FT) and short-time Fourier transform (STFT) are commonly used as data analysis methods, the principle of which will not be repeated here. This section mainly introduces the CEEMD and Hilbert transform.

CEEMD
CEEMD comes from EMD and EEMD. The EMD method resolves a given signal into a series of IMFs. These IMFs are components that satisfy certain conditions. The reconstructed original signal is [22] where i c is the ith IMF and n r is the residue of data x after n IMFs are extracted. The decomposition flowchart of EMD is shown in Figure 3. The termination conditions of the EMD in the figure are: (1) The number of extreme points is the same as or, at most, one of the zero crossings; (2) the average value of the envelope defined by the local extremum at any point is zero. There are many kinds of space electromagnetic waves on the moon; the terrain of the moon is very complex under the adverse circumstances, so the signal-to-noise ratio (SNR) of LPR data is low. The repeated gathers of each stop point are superimposed, the purpose of which is to improve the SNR of the single trace at the waypoint [24], and provide more accurate and reasonable one-dimensional data for the subsequent analysis of the frequency-domain (FD) and the time-frequency domain (TFD).
In this paper, several analysis and processing methods are utilized, for example, the Fourier transform, short-time Fourier transform, CEEMD, and Hilbert transform. The Fourier transform (FT) and short-time Fourier transform (STFT) are commonly used as data analysis methods, the principle of which will not be repeated here. This section mainly introduces the CEEMD and Hilbert transform.

CEEMD
CEEMD comes from EMD and EEMD. The EMD method resolves a given signal into a series of IMFs. These IMFs are components that satisfy certain conditions. The reconstructed original signal is [22] where c i is the ith IMF and r n is the residue of data x after n IMFs are extracted. The decomposition flowchart of EMD is shown in Figure 3. The termination conditions of the EMD in the figure are: (1) The number of extreme points is the same as or, at most, one of the zero crossings; (2) the average value of the envelope defined by the local extremum at any point is zero. As mentioned above, the main disadvantage of EMD is the frequent occurrence of mode overlap, which is defined as a single IMF composed of signals of different scales, or a signal of likeness scale in varying IMF components. In order to overcome the frequent appearance of mode overlap, several new noise-assisted data analysis methods are proposed, such as variational mode decomposition (VMD) [25] and ensemble EMD (EEMD) [22]. The EEMD method can be expressed as follows [23]:

1.
Generate where w i [n] (i = 1, . . . , I) are different implementations of white Gaussian noise, Set up IMF k to be the k-th pattern of x[n], obtained as the corresponding average IMF i k [23]: The resulting decomposition of EEMD solves the problem of mode mixing; however, it introduces new ones. In the technique of CEEMD, the decomposition mode will be recorded as IMF k and be proposed to count a unique first residue as (4) is obtained in the same way as EEMD. Then, calculate the first EMD mode over aggregate r 1 [n] plus various realizations of a given noise, obtaining IMF 2 through averaging. The next residue can be defined as As mentioned above, the main disadvantage of EMD is the frequent occurrence of mode overlap, which is defined as a single IMF composed of signals of different scales, or a signal of likeness scale in varying IMF components. In order to overcome the frequent appearance of mode overlap, several new noise-assisted data analysis methods are proposed, such as variational mode decomposition (VMD) [25] and ensemble EMD (EEMD) [22]. The EEMD method can be expressed as follows [23]: The resulting decomposition of EEMD solves the problem of mode mixing; however, it introduces new ones. In the technique of CEEMD, the decomposition mode will be recorded as  k IMF and be proposed to count a unique first residue as First, we define an operator E j (·), which, given a signal, produces the j-th mode from EMD [23]. w i is the white noise with N(0, 1). If x[n] is the target datastore, the method can be described by the following algorithm [23]:

1.
Disintegrate by EMD I realizations x[n] + ε 0 w i [n] to obtain their first modes and calculate 2.
In the first phase (k = 1), obtain the first residue as in Equation (4).

5.
Decompose realizations r k [n] + ε k E k w i [n] , i = 1, . . . , I, until their first EMD mode and define the (k + 1)-th mode as [23] 6. Go to step 4 to continue to get the next k.
Repeat steps 4-6. The direct residual part cannot be further decomposed, and the final residue satisfies [23] R K is the total number of modes. Therefore, the given signal x[n] can be written as Equation (10) indicates that the data decomposition has been completed and the original data can be reconstructed.

Hilbert Transform
After introducing the algorithm of CEEMD, we describe the principle of the Hilbert transform. For any time series X(t), we can get its Hilbert transform result: The Hilbert inverse transformation is For a real signal X(t), there are many kinds of complex methods, and it is the simplest and most effective to obtain the analytical signal by the Hilbert transform method. Assuming that the Hilbert transform of X(t) is Y(t), the expression of the parsed signal is where a(t) and θ(t) are the instantaneous amplitude and phase of the signal, respectively, and they can be written in the following form: The signal's instantaneous frequency can be obtained by deriving the instantaneous phase: Unlike the Fourier transform, a and ω in the equation are not constant, but a function of time, so the Hilbert spectrum can characterize the temporal variation of a data sequence; the instantaneous frequency is relative to the analytical signal, which is the analytical signal phase.
The frequency obtained by the signal FT reflects the distribution characteristics of the signal with time, and the frequency of the signal with time cannot be obtained. The instantaneous frequency reflects the frequency of the signal at a certain time, which can better describe the characteristics of the signal. This is also the main difference between them. The Hilbert transform can only deal with the stationary signal of a single frequency, which has great limitations for non-stationary signals with multiple frequencies at a certain time. Therefore, the Hilbert transform must meet some constraints in the application to obtain a physical instantaneous frequency.

Analysis in Frequency Domain and the Time-Frequency Domain
The FT and the STFT are performed to obtain the frequency spectrum ( Figure 4) and the time-frequency spectrum ( Figure 5) (et al., N104, N106, N202, and N206, for examples), and, thus, the time-frequency spectra will expose the issue and the kind of processing that should be done.
Unlike the Fourier transform, a and ω in the equation are not constant, but a function of time, so the Hilbert spectrum can characterize the temporal variation of a data sequence; the instantaneous frequency is relative to the analytical signal, which is the analytical signal phase.
The frequency obtained by the signal FT reflects the distribution characteristics of the signal with time, and the frequency of the signal with time cannot be obtained. The instantaneous frequency reflects the frequency of the signal at a certain time, which can better describe the characteristics of the signal. This is also the main difference between them. The Hilbert transform can only deal with the stationary signal of a single frequency, which has great limitations for non-stationary signals with multiple frequencies at a certain time. Therefore, the Hilbert transform must meet some constraints in the application to obtain a physical instantaneous frequency.

Analysis in Frequency Domain and the Time-Frequency Domain
The FT and the STFT are performed to obtain the frequency spectrum ( Figure 4) and the timefrequency spectrum ( Figure 5) (et al., N104, N106, N202, and N206, for examples), and, thus, the timefrequency spectra will expose the issue and the kind of processing that should be done. It can be seen from the frequency spectra of the sample traces ( Figure 4) that in the frequency domain, the maximum energy of all sample gathers is around 12 MHz (as indicated by the arrow), whereas at around 60 MHz (the center frequency of CH-1), the energy is not prominent.  It can be seen from the frequency spectra of the sample traces ( Figure 4) that in the frequency domain, the maximum energy of all sample gathers is around 12 MHz (as indicated by the arrow), whereas at around 60 MHz (the center frequency of CH-1), the energy is not prominent.
Take a look at the time-frequency spectra of the sample traces: (1) As the arrows in the N104 time-frequency spectrum of STFT show, there are 150, 100, and 60 MHz electromagnetic energies in the data throughout the entire time delay. Obviously, the energies at 150 and 100 MHz are not the transmitting signals. The electromagnetic energies at those frequencies may have come from the lunar space or the instrument (needs further research). (2) Compared to the shallow signal, the energy of the deep electromagnetic signal is very weak. Gain should be processed to highlight the deep information ( Figure 6). (3) Both the energies of the high-frequency and low-frequency information in the shallow part are strong. As the time delay increases, high-frequency energy is absorbed, and only some low-frequency information is left. For this kind of nonlinear frequency aliasing phenomenon, a nonlinear decomposition method must be processed. Take a look at the time-frequency spectra of the sample traces: (1) As the arrows in the N104 time-frequency spectrum of STFT show, there are 150, 100, and 60 MHz electromagnetic energies in the data throughout the entire time delay. Obviously, the energies at 150 and 100 MHz are not the transmitting signals. The electromagnetic energies at those frequencies may have come from the lunar space or the instrument (needs further research). (2) Compared to the shallow signal, the energy of the deep electromagnetic signal is very weak. Gain should be processed to highlight the deep information ( Figure 6). (3) Both the energies of the high-frequency and low-frequency information in the shallow part are strong. As the time delay increases, high-frequency energy is absorbed, and only some lowfrequency information is left. For this kind of nonlinear frequency aliasing phenomenon, a nonlinear decomposition method must be processed.

CEEMD of the Trace at N104
In terms of the N104 trace, CEEMD is performed (Figure 7). The noise standard deviation is set to be 0.2. The number of realizations is 200 and the maximum number of sifting iterations allowed is 2000. The N104 trace is divided into seven IMFs and one residual.
After decomposing the trace at N104, in order to observe the time-frequency information of the data in each component, we performed a short-time Fourier transform for each component (Figure 8). Analyze each component one by one: (1) In the STFT spectrum of the IMF1 component, energy is concentrated mainly at 150 MHz, and a small amount of energy is at 100 and 60 MHz; this component can be considered a high-frequency noise component. (2) In the STFT spectrum of the IMF2 and IMF3 components, the main energy is collected near 60 MHz, and some of the energy scattered around it is considered to be the component containing the useful signal.  Figure 7) and its STFT spectrum (Figure 8e), which are around 3500 and 5800 ns, respectively. The reflection of these two positions was used for reflection in the articles by Xiao et al. [6] and Zhang et al. [10], but it can be seen from the STFT spectrum that the frequency is between 10 and 15 MHz, and it is quite different from the center frequency of the transmitter. Li et al. also used the S transform and LPR ground data to verify the trap of these two pitfalls [20]. In summary, the IMF5 component is considered to be the unusable component. (5) The STFT spectra of IMF5, IMF6, and IMF7 and the residual component show that the energy is collected below 10 MHz, which is a low-frequency unwanted component.

CEEMD of the Trace at N104
In terms of the N104 trace, CEEMD is performed (Figure 7). The noise standard deviation is set to be 0.2. The number of realizations is 200 and the maximum number of sifting iterations allowed is 2000. The N104 trace is divided into seven IMFs and one residual.  Figure 7) and its STFT spectrum (Figure 8e), which are around 3500 and 5800 ns, respectively. The reflection of these two positions was used for reflection in the articles by Xiao et al. [6] and Zhang et al. [10], but it can be seen from the STFT spectrum that the frequency is between 10 and 15 MHz, and it is quite different from the center frequency of the transmitter. Li et al. also used the S transform and LPR ground data to verify the trap of these two pitfalls [20]. In summary, the IMF5 component is considered to be the unusable component. (5) The STFT spectra of IMF5, IMF6, and IMF7 and the residual component show that the energy is collected below 10 MHz, which is a low-frequency unwanted component.
On the other hand, from the Hilbert spectrum of IMFs at N104 (Figure 9), it can be seen that the energy in the shallow part is mainly distributed in the low-frequency IMF components below 20 MHz. The deep part mainly includes three parts, a high-frequency IMF component of 150 MHz, an IMF component of 60 MHz, which is consistent with the transmission frequency, and a low-frequency pitfall component of 20 MHz or less. Here again, an aliasing component around 3500 and 5800 ns can be found.  On the other hand, from the Hilbert spectrum of IMFs at N104 (Figure 9), it can be seen that the energy in the shallow part is mainly distributed in the low-frequency IMF components below 20 MHz. The deep part mainly includes three parts, a high-frequency IMF component of 150 MHz, an IMF component of 60 MHz, which is consistent with the transmission frequency, and a low-frequency pitfall component of 20 MHz or less. Here again, an aliasing component around 3500 and 5800 ns can be found.   After decomposing and analyzing the single trace at N104 as an example, we perform CEEMD on all CH-1 data ( Figure 10) and the traces of stagnation points ( Figure 11) in order to obtain a deeper understanding and utilization of the CH-1 data.
After comprehensively observing the CEEMD components of the overall CH-1 data and the CEEMD components of the stagnation points' traces, it can be found that in IMF5 (shown by the arrows in Figures 10 and 11), two pitfalls at 3500 and 5800 ns can be found (IMF5 in the stagnation points' traces is more obvious than IMF5 in the overall CH-1 data).

Data Interpretation
As can be known from the above analysis, only the frequencies of the IMF2 and IMF3 components are in accordance with the radar's transmission frequency and are available. These two components are superimposed for further analytical interpretation. Figure 12a shows the data after IMF3 is superimposed with IMF4, and Figure 12b plots all the traces. From these two figures, especially in Figure 12b, it is found that the SNR of the data below 3000 ns is very low, making it difficult to pick up effective and convincing events. Therefore, we will focus on the shallow part. Figure 13 shows the information within the shallow part of 3000 ns, but the events below 300 ns are not continuous. It is necessary to be very careful to interpret such information. This paper aims to decompose and analyze the CH-1 data, and the deep information no longer needs deeper explanation. Instead, information within the first 200 ns corresponding to the CH-2 data is selected for analysis. According to the hierarchical information of Xiao et al. [6], Zhang et al. [10], Fa et al. [11], and Zhang et al. [14], there are three well-recognized interfaces (Figure 14), namely: 1.
The interface between the ejected materials of the Chang'E-3 crater and the lunar regolith; 2.
The interface between the lunar regolith and paleoregolith; 3.
The bedrock interface. on all CH-1 data ( Figure 10) and the traces of stagnation points ( Figure 11) in order to obtain a deeper understanding and utilization of the CH-1 data. After comprehensively observing the CEEMD components of the overall CH-1 data and the CEEMD components of the stagnation points' traces, it can be found that in IMF5 (shown by the arrows in Figures 10 and 11), two pitfalls at 3500 and 5800 ns can be found (IMF5 in the stagnation points' traces is more obvious than IMF5 in the overall CH-1 data).

Data Interpretation
As can be known from the above analysis, only the frequencies of the IMF2 and IMF3 components are in accordance with the radar's transmission frequency and are available. These two components are superimposed for further analytical interpretation. Figure 12a shows the data after IMF3 is superimposed with IMF4, and Figure 12b plots all the traces. From these two figures, especially in Figure 12b, it is found that the SNR of the data below 3000 ns is very low, making it difficult to pick up effective and convincing events. Therefore, we will focus on the shallow part.  Figure 13 shows the information within the shallow part of 3000 ns, but the events below 300 ns are not continuous. It is necessary to be very careful to interpret such information. This paper aims to decompose and analyze the CH-1 data, and the deep information no longer needs deeper explanation. Instead, information within the first 200 ns corresponding to the CH-2 data is selected for analysis. According to the hierarchical information of Xiao et al. [6], Zhang et al. [10], Fa et al. [11], and Zhang et al. [14], there are three well-recognized interfaces (Figure 14), namely: 1 The interface between the ejected materials of the Chang'E-3 crater and the lunar regolith; 2 The interface between the lunar regolith and paleoregolith; 3 The bedrock interface.    Figure 13 shows the information within the shallow part of 3000 ns, but the events below 300 ns are not continuous. It is necessary to be very careful to interpret such information. This paper aims to decompose and analyze the CH-1 data, and the deep information no longer needs deeper explanation. Instead, information within the first 200 ns corresponding to the CH-2 data is selected for analysis. According to the hierarchical information of Xiao et al. [6], Zhang et al. [10], Fa et al. [11], and Zhang et al. [14], there are three well-recognized interfaces (Figure 14), namely: 1 The interface between the ejected materials of the Chang'E-3 crater and the lunar regolith; 2 The interface between the lunar regolith and paleoregolith; 3 The bedrock interface.

Discussion
In the use of CEEMD to process radar data, Li et al. used the CEEMD decomposition method not only to perform time-frequency analysis but also to perform denoising and target recognition on GPR data [26,27]. They mainly used CEEMD to decompose GPR data, used the Hilbert transform to analyze the data, and then performed denoising accordingly. The method in this paper not only uses CEEMD to decompose the data, but also uses the STFT and Hilbert transform to analyze the single component and the overall data comprehensively, which improves the application of CEEMD in ground-penetrating-radar processing.
On the other hand, regarding the use of CH-1 data from the lunar radar, although Li et al. [20] believed that Xiao et al. [6] and Zhang et al. [10] had objections to the use of CH-1 data, and used the ground experimental data and S transform to support their view, they did not perform in-depth decomposition and analysis of the overall CH-1 data. In this paper, the non-linear decomposition method is used to deeply analyze the false reflectors and available information of CH-1 data.   Figure 15 shows the information within 200 ns of the CH-1 data. There are three distinct events at the same position corresponding to the CH-2 data. It proves the rationality of stratification at the shallow part.

Discussion
In the use of CEEMD to process radar data, Li et al. used the CEEMD decomposition method not only to perform time-frequency analysis but also to perform denoising and target recognition on GPR data [26,27]. They mainly used CEEMD to decompose GPR data, used the Hilbert transform to analyze the data, and then performed denoising accordingly. The method in this paper not only uses CEEMD to decompose the data, but also uses the STFT and Hilbert transform to analyze the single component and the overall data comprehensively, which improves the application of CEEMD in ground-penetrating-radar processing.
On the other hand, regarding the use of CH-1 data from the lunar radar, although Li et al. [20] believed that Xiao et al. [6] and Zhang et al. [10] had objections to the use of CH-1 data, and used the ground experimental data and S transform to support their view, they did not perform in-depth decomposition and analysis of the overall CH-1 data. In this paper, the non-linear decomposition method is used to deeply analyze the false reflectors and available information of CH-1 data.

Discussion
In the use of CEEMD to process radar data, Li et al. used the CEEMD decomposition method not only to perform time-frequency analysis but also to perform denoising and target recognition on GPR data [26,27]. They mainly used CEEMD to decompose GPR data, used the Hilbert transform to analyze the data, and then performed denoising accordingly. The method in this paper not only uses CEEMD to decompose the data, but also uses the STFT and Hilbert transform to analyze the single component and the overall data comprehensively, which improves the application of CEEMD in ground-penetrating-radar processing.
On the other hand, regarding the use of CH-1 data from the lunar radar, although Li et al. [20] believed that Xiao et al. [6] and Zhang et al. [10] had objections to the use of CH-1 data, and used the ground experimental data and S transform to support their view, they did not perform in-depth decomposition and analysis of the overall CH-1 data. In this paper, the non-linear decomposition method is used to deeply analyze the false reflectors and available information of CH-1 data.
The most important significance of this paper is that the LPR data of CE-4 is being collected on the moon, and the CH-1 data collected by CE-4 can be used by the CEEMD combined with the short-time Fourier transform and the Hilbert in the case of an extremely low SNR. The idea of special transformation is analyzed to improve the credibility of the CE-4 data.

Conclusions
The LPR equipped on the YR measured the lunar geological structure in Northern Imbrium. As the data were complicated, a group of data preprocessing was designed. When the data were still distinct after preprocessing, the frequency spectrum and time-frequency spectrum of LPR CH-1 data were carried out for analysis. Subsequently, CEEMD helped to analyze the data, and the IMFs indicated the information distribution. Finally, combining with the previous studies and LPR CH-2 data, we proved the credibility of the CH-1 data at the shallow part.
CEEMD has been put forward as an adaptive time-frequency LPR data analysis method. It has been verified as quite multifunctional in extracting signals. The IMFs of the overall CH-1 data and the stagnation points' traces describe the frequency information distribution. After comprehensively observing the CEEMD components of the overall CH-1 data and the CEEMD components of the stagnation points' traces, it can be found that in IMF5, two pitfalls at 3500 and 5800 ns can be found. Only the frequencies of the IMF2 and IMF3 components are in accordance with the radar's transmission frequency and are available. These two components are superimposed for further analytical interpretation. Superimposed components provide us with the chance to analyze geological information at the shallow part.
This paper aims to decompose and analyze the CH-1 data, and the deep information no longer needs a deeper explanation. Instead, information within the first 200 ns corresponding to the CH-2 data is selected for analysis. According to the hierarchical information of former researchers, there are three well-recognized interfaces. As the CE-4 LPR data are going to be published soon, CEEMD, which can prove the availability of the data, provides solid help and support for the articles published on the CE-4 radar data.
Those results obtained in this paper provide valuable reference for understanding LPR data and valuable information for future lunar exploration missions.

Conflicts of Interest:
The authors declare no conflict of interest.