A Morphology Filter-Assisted Extreme-Point Symmetric Mode Decomposition (MF-ESMD) Denoising Method for Bridge Dynamic Deflection Based on Ground-Based Microwave Interferometry

Bridge dynamic deflection is an important indicator of structure safety detection. Ground-based microwave interferometry is widely used in bridge dynamic deflection monitoring because it has the advantages of noncontact measurement and high precision. However, due to the influences of various factors, there are many noises in the obtained dynamic deflection of bridges obtained by ground-based microwave interferometry. To reduce the impacts of noise for bridge dynamic deflection obtained with ground-based microwave interferometry, this paper proposes a morphology filter-assisted extreme-point symmetric mode decomposition (MF-ESMD) for the signal denoising of bridge dynamic deflection obtained by ground-based microwave interferometry. First, the original bridge dynamic deflection obtained with ground-based microwave interferometry was decomposed to obtain a series of intrinsic mode functions (IMFs) with the ESMDmethod. Second, the noise-dominant IMFs were removed according to Spearman’s rho algorithm, and the other decomposed IMFs were reconstructed as a new signal. Finally, the residual noises in the reconstructed signal were further eliminated using the morphological filter method. -e results of both the simulated and on-site experiments showed that the proposed MF-ESMD method had a powerful signal denoising ability.


Introduction
Bridges have become one of the important elements of modern urban transportation. However, due to a combination of various factors such as age, environment, human behavior, overload, and geological activities, corrosion and cracks may occur on the surfaces of bridges and increase the deterioration of the bridges, which may ultimately lead to reducing the performance the bridges or even sudden collapses [1]. Serious damage to the bridges may result in serious property damage and casualties. erefore, it is an essential and extremely important task to carry out the structural health monitoring of the bridges. With the advantage of reflecting the overall stiffness of a bridge structure [2], bridge dynamic deflection has become an important indicator for evaluating the stability of a monitored bridge, which is closely related to its bearing capacity and ability to resist dynamic loads such as strong winds and earthquakes [2][3][4][5]. By analyzing the dynamic deflection of monitored bridges, the internal force distributions can be obtained, and the weak position of the bridge structures can be further determined [6]. erefore, it is necessary to ascertain highprecision bridge dynamic deflections of the monitored bridges.
Compared with the traditional techniques of bridge dynamic deflection measurement such as level gauges and displacement meters, the new technology of ground-based microwave interferometry has the great advantages of noncontact displacement measurement, a wide frequency range of responses, submillimetric displacement sensitivity, and a quick setup [7][8][9][10]. It has been widely used to obtain the dynamic deflection measurements of bridges in recent years [11,12]. Stabile et al. [8] applied the technique to estimate the fundamental dynamic parameters of the reinforced concrete Musmeci Bridge in Basilicata Region (southern Italy), and the results have been validated by the comparison with the ones obtained applying consolidated techniques using data from accelerometers and tromometers. Zhou et al. [12] obtained the subsidence time series of the East Lake High-tech Bridge through groundbased interferometric radar, and a case study was conducted on the subsidence of the East Lake High-tech Bridge during subway shield tunnel crossing underneath this bridge. However, during the data acquisition of bridge dynamic deflection using ground-based microwave interferometry, various factors such as the surrounding environment, thrust of wind, ground-motion, complex traffic activities, artificial operation, and the equipment itself will inevitably generate noise in the collected bridge dynamic deflection. ese factors can cause the loss of precision of the obtained bridge dynamic deflection and further lead to errors in data analysis when evaluating the stability of the monitored bridge [5,13]. erefore, it is of great importance to reduce the influence of noise in the bridge dynamic deflection obtained by groundbased microwave interferometry.
At present, wavelet transform is one of the widely used methods for signal denoising with a good time-frequency localization characteristic [14]. is method can be used to obtain the wavelet coefficients of signals and noises with different properties at each scale [15]. Moreover, with different wavelet threshold criteria, wavelet transform can be used to distinguish different signals with high frequencies and low frequencies. Wavelet transform has proven to be a useful tool for nonstationary signal analysis [16][17][18]. However, the results of wavelet transform can be affected by the selection of the wavelet basis function and wavelet threshold criteria, which may cause the distortion of a signal in the time domain [19,20]. Moreover, wavelet transform is not suitable for nonlinear signals [20], which means it is not suitable for analyzing the dynamic deflection of a monitored bridge because if the monitored bridge were damaged, the obtained time-series displacement would be a nonlinear signal.
Recently, the empirical mode decomposition (EMD) method has been widely applied to analyze nonlinear and nonstationary vibration signals.
is method works to process a signal by performing a timed adaptive decomposition operation that can adaptively decompose a nonstationary signal into a series of intrinsic mode functions (IMFs) [20]. e decomposed components contain the local characteristic signals for different time scales of the source signal [21]. e influence of noise in the obtained vibration signals can be reduced by eliminating a certain number of low-order IMFs with relatively high frequencies [22]. However, the mode mixing effect and the rough trend function can bring out noise with different scales in the decomposed IMFs, which may cause a loss of precision for the EMD denoised method [20,23]. To solve this problem, the ensemble empirical mode decomposition (EEMD) method has been proposed based on the research of the statistical properties of a Gaussian white noise signal, which is an improvement based on the EMD method [21,[24][25][26][27][28].
is method changes the extremum distribution of the original signal by adding white Gaussian noise and further decomposes the original signal to alleviate the mode-mixing effect using the EMD method [21,[24][25][26][27][28][29][30]. However, the EEMD method is easily affected by the amplitude of the added noise and the number of ensemble trials, which can bring out some residual noise in the decomposed IMFs. Moreover, although the EMD and EEMD methods can effectively reduce the influence of periodic noise, it is difficult to reduce the influence of instantaneous noise, such as the instantaneous vibrations of the equipment itself caused by passing vehicles [31]. Extreme-point symmetric mode decomposition (ESMD) is a new alternative time-frequency analysis method for the EMD and EEMD methods [32]. e main advantage of this method is the ability to adaptively determine the optimal global average curve. Any complex raw signal can be decomposed into a series of different physical meaning IMFs, which means that the ESMD method is superior to the EMD and the EEMD methods [7,20,33]. It has been applied in various fields, such as information science, marine and atmospheric science, and seismology [31][32][33][34][35]. Although the ESMD method can effectively reduce the influence of the mode-mixing problems, it is difficult to reduce the influence of noise with different frequencies and scales in decomposed IMFs using the ESMD method.
A mathematical morphology filter is a very typical nonlinear signal processing and analysis technology that has been rapidly developed in recent years [36][37][38][39]. is technology has been widely applied in the analysis of signals such as vibration signals and pulse wave signals [40][41][42]. Morphological filtering based on the theory of mathematical morphology is a new nonlinear filtering method for eliminating noise while protecting signal details. e theory of mathematical morphology is an integrated theory of information theory, statistical estimation theory, adaptive theory, random noise theory, and other theories [36][37][38][39]. Morphological filtering is a noise removal algorithm that is suitable for all types of noise and has a better elimination effect [41]. Compared with traditional linear filtering, morphological filtering is more efficient to a certain extent, and it overcomes the shortcomings of linear filtering in which some primary features of the original signal are harmed or blurred while smoothing the noise. us, the main morphological features of each signal remain more intact [29,36]. Moreover, a morphological filter is more effective in reducing impulse noise and white noise, and it will not cause an abrupt change in the frequency domain or a phase delay in the time domain [29,43]. To address modemixing problems, Zhou et al. [29] proposed a morphological filter-assisted ensemble empirical mode decomposition (MF-EEMD) for signal denoising that could significantly mitigate the mode-mixing problems and achieve a higher decomposition efficiency compared to the EEMD method. erefore, in this study, to effectively eliminate the influences of noise in the dynamic deflection of the monitored bridges obtained by ground-based microwave interferometry, a morphology filter-assisted extreme-point symmetric mode decomposition (MF-ESMD) denoising method for bridge dynamic deflection obtained by ground-based microwave interferometry was proposed. e ESMD method was used to decompose the original obtained bridge dynamic deflection into a series of IMFs, and the noisedominant IMFs were further removed according to Spearman's rho algorithm. A morphological filter was used to further eliminate the residual noises in the reconstructed signal. e MF-ESMD denoising method is introduced in Section 2. Section 3 introduces the results and analysis of the simulation experiments, and Section 4 introduces the results and analysis of the on-site experiment, followed by a conclusion in Section 5. Figure 1 shows the entire workflow of the proposed MF-ESMD denoising method for the bridge dynamic deflection obtained by ground-based microwave interferometry. It consisted of the following two technologies: (1) the original dynamic deflection signal obtained by ground-based microwave interferometry was decomposed into a series of IMFs using the ESMD method. Spearman's rho of each decomposed IMF was calculated to remove the noisedominant IMFs from the decomposed IMFs. (2) A new signal was reconstructed by the remaining IMFs, and the morphological filter method was used to remove the residual noise from the reconstructed signal.

Reduction of Noise-Dominant IMFs.
e ESMD method has a better performance in decomposing a signal into a series of physically meaningful representations compared with the EMD and EEMD methods, which can be more effective to separate noise and useful signals [20]. Moreover, as a nonparametric coefficient correlation analysis method, during the process of the calculation of data correlations, Spearman's rho only depends on ranks, which can guarantee enough robustness of signal noise processing [30,44]. erefore, in this study, the ESMD method and Spearman's rho were integrated in order to remove the noise-dominant IMFs from the decomposed IMFs for the bridge dynamic deflection obtained by ground-based microwave interferometry. For the original dynamic deflection signal y(t), the detailed process of the reduction of noise-dominant IMFs is shown as follows: Step 1: All the maximum and minimum extreme points of the original dynamic deflection signal y(t) were found, and the midpoints of the two adjacent extreme points were defined as F i (i � 1, 2, . . . , n − 1). e midpoints of the left and right borders were set to F 0 and F n using a linear interpolation method. Interpolation curves (named L 1 , L 2 , . . . , L m (m ≥ 1)) were curved using cubic spline interpolation based on the N + 1 points in F i (i � 0, 1, 2, . . . , n), and their average curve L � (L 1 + L 2 + · · · + L m )/m was further calculated.
Step 2: y(t) − L was calculated by repeating Step 1 to obtain IMF 1 until |L| ≤ ε was satisfied or the sifting times reached the preset maximum value K p of 30. In this study, the permitted error ε was set to 0.001 σ 0 , where σ 0 was the standard deviation of the original data y(t) relative to the total mean of y(t).
Step 3: Steps 1 and 2 were repeated for y(t) − IMF i to obtain IMF 2 , . . ., IMF n and a residual R, which is an optimal adaptive global mean (AGM) curve with no more than a certain number of extreme points. en, the variance σ of y(t) − R was calculated.
Step 4: e sifting times K was changed on a finite integer interval [K min , K max ].
en, Steps 1-3 were repeated to obtain the optimal sifting time K 0 , according to the minimum variance σ min on [K min , K max ].
Step 5: e optimal sifting times value K 0 and the permitted error ε were used to perform the optimal signal decomposition for the dynamic deflection signal y(t) to obtain a series of IMFs with an optimal AGM curve.
Step 6: For each decomposed IMF, the noise-dominant IMFs had features that were less similar to the features of the original dynamic deflection signal. Spearman's rho between each decomposed IMF and the original dynamic deflection signal is between [−1, +1]. e bigger the coefficient value, the more similar it was to the original signal [30,44].
rough the coefficient value, noise-dominant IMFs could be removed. Spearman's rho s(i) between the decomposed IMF i and the original signal X(t) could be calculated using equation (1): where IMF i is the i th IMF (i � 1, 2, . . . , n) and N is the number of elements of the original dynamic deflection signal.

Reduction of the Residual Noise.
After the reduction of the noise-dominant IMFs, the other decomposed IMFs were reconstructed as a new signal. ere was only a small amount of residual noise in the reconstructed signal, which may be caused by the surrounding environment, artificial operation, and the equipment itself. A morphology filter is used to perform the reduction of the residual noise in the reconstructed signal. Generally, dilation and erosion are two most basic operations in the algorithm of the mathematical morphology method [14][15][16][17][18]. A was defined as a signal, and B was defined as a structuring element. e dilation operation A ⊕ B was defined as shown in equation (2), and the erosion operation A Θ B was defined as shown in equation (3). Based on the dilation and erosion operations, the opening operation A ∘ B was defined as shown in equation Shock and Vibration (4), and the closing operation A · B was defined as shown in equation (5): However, the opening operation could only change the maxima and maintain the shape of the minima for the reconstructed signal. e closing operation could only change the minima and maintain the shape of the maxima for the reconstructed signal. erefore, a good filtering result could not be obtained if any operation was used on its own. However, the average value of the two operations could be very close to the original signal [36]. To address both the minima and maxima of the signal, a morphological filter was proposed in order to reduce the residual noise in the reconstructed signal, as shown in equation (6), which integrated the opening and closing operations: where s(t) is the final denoised signal, y(t) * is the reconstructed signal, and n is the length of signal y(t) * . g 1 and g 2 are the two structuring elements with different sizes, of which the lengths were less than the length of the reconstructed signal y(t) * . In addition, the structuring element was an important factor affecting the results of the signal denoising. e most frequently used structuring elements include linear structuring elements, triangular structuring elements, and structuring elements [36][37][38][39]. e filtering results were affected by the length and shape of the structuring elements. erefore, in the actual processing of the signal, the structuring elements were designed according to the characteristics of the analysis signal and the computational complexity [36][37][38][39]. In this study, the linear structuring elements were selected due to the simplicity in calculation, and the shape had little effect on the analysis [27].

Simulated Experiments and Analysis
e dynamic deflection signal of bridges obtained by ground-based microwave interferometry is a complicated The original dynamic deflection signal y(t) obtained by ground-based microwave interferometry The remaining IMF components are regarded as useful subsignals to be reconstructed as a new signal y(t) * together with the residual R The final denoised signal s(t) is output The original dynamic deflection signal y(t) is decomposed into a series of IMFs (IMF 1 , IMF 2 , ..., IMF n ) and a residual R using the ESMD method The Spearman coefficientis calculated between each decomposed IMF and the original dynamic deflection signal y(t) Spearman's rho of IMF i (i = 1, 2, ..., i < n) is the minimum value for the first time, and then IMF 1 to IMF i are regarded as the noise dominant IMFs

Reduction of noise dominant IMFs
Reduction of the residual noise The reconstructed signal y(t) * is input Morphology filter is used to perform the reduction of the residual noise in the reconstructed signal y(t) * nonstationary signal mixed with some noise. erefore, in this study, to validate the feasibility and accuracy of the proposed MF-ESMD denoising method, a nonlinear and nonstationary simulated signal Y(t) was generated that consisted of three basic functions with the different frequencies X(t) � X 1 + X 2 + X 3 and a white Gaussian noise signal with an SNR (signal-to-noise ratio) of 20 dB to enhance noise information in the simulated signal Y(t), which is shown as follows: where X 1 � cos(20πt + 2 sin(3πt)), X 2 � sin(6πt), X 3 � cos(10πt), and W(t) is a white Gaussian noise signal. e waveforms of the simulation signal Y(t) and its component signals X 1 , X 2 , and X 3 are shown in Figure 2. As shown in Figure 3, the simulated signal Y(t) was decomposed into eight IMFs with the corresponding meaningful inherent natural frequency from high to low frequency and a residual R using the ESMD method. en Spearman's rho of each decomposed IMFs was calculated, as shown in Table 1. e value of Spearman's rho from the first decomposed IMF 1 was gradually reduced. Spearman's rho of IMF 3 was 0.0291, which was the minimum value for the first time. is indicated that the decomposed IMF 1 to IMF 3 had less similar feathers than the original signal Y(t). erefore, the first three IMF components, IMF 1 to IMF 3 , were regarded as the noise-dominant IMFs, which may have been caused by the surrounding environment, thrust of wind, ground-motion, complex traffic activities, artificial operation, and the equipment itself. e remaining IMF components were regarded as useful subsignals to be reconstructed as a new signal Y(t) * along with the residual R.
For the reconstructed signal Y(t) * , most of the obvious noises were reduced. However, there was a small amount of residual noise in the reconstructed signal Y(t) * . erefore, a morphology filter was used to perform the reduction of the residual noise in the reconstructed signal Y(t) * . In the morphology filter, g 1 � 1 1 1 and g 2 � 1 1 1 1 1 were denoted as shown in equation (6). en, the final denoised signal S(t) could be obtained using the morphology filter. Figure 4 shows the compared results of the curves between the original signal X(t) and the denoised signal S(t). Moreover, to validate the accuracy and robustness of the proposed MF-ESMD method for signal denoising, morphological filter and MF-EEMD methods [27,29] were selected for comparison with the proposed MF-ESMD method. e comparison results of the curves between the original signal X(t) and the denoised signal using the morphological filter and the MF-EEMD methods are shown in Figures 5(a) and 5(d). e inspection of the curves from Figures 4 and 5 highlighted the following aspects: (1) for the middle parts of the curve, as shown in Figures 4(b), 5(b) and 5(f ), the middle parts of the denoised signal were generally consistent with the original signal using the above three methods. e results indicated that the noise could be eliminated effectively using the above three methods for the middle parts of the curve of the simulated signal. (2) For the peak parts of the curve, although there was a maximum deviation of 0.006 mm, the curve of the denoised signal was generally consistent with the original signal using the proposed MF-ESMD method, as shown in Figures 4(c) and 4(d). However, for the denoised signal using the morphological filter and MF-EEMD methods, as shown in Figures 5(c), 5(d), 5(g), and 5(h), there was obvious noise in the peak parts of the curve. e results indicated that the proposed MF-ESMD method had a better performance for peak signal denoising than the morphological filter and the MF-EEMD method.
In order to further evaluate the quality of signal denoising using the proposed MF-ESMD method, two objective evaluation indices, the SNR (signal-to-noise ratio) and RMSE (root-mean-square error), were used in this study. e SNR was a ratio of the useful signal power to the noise power, as shown in equation (8). Generally, the greater the SNR, the better the denoised effect. e SNR could reflect the proportional relationship between the effective part of the original signal and the noise, and it could reflect the denoising ability of the algorithm [7,45].
e RMSE was a deviation value between the observed value and the true value, as shown in equation (9). Generally, the smaller the RMSE value, the better the denoised effect. e RMSE could indicate the difference between the original signal and the denoised signal [7,20]. e bigger the SNR value, the smaller the RMSE value and the better the effect of noise reduction: where y(n) is the n th element of the original signal, s(n) is n th the element of the denoised signal, and N is the number of elements of the signal. e comparison results are shown in Table 2. e inspection from this table highlighted the following aspects: (1) the RMSE was 0.0361 mm for the proposed MF-ESMD denoising method, which was smaller and closer to zero than the RMSE values of the morphological filter method and the MF-EEMD method. (2) e SNR of the simulation signal Y(t) was 20 dB. e SNR was 26.3066 dB for the denoised signal using the morphological filter method, which increased by 31.53% compared to the simulation signal Y(t).
e SNR was 26.4298 dB for the denoised signal using the MF-EEMD method, which increased by 32.15% compared to the simulation signal Y(t). e SNR was 30.6619 dB for the denoised signal using the proposed MF-ESMD denoising method, which increased by 53.3% compared to the simulation signal Y(t). e results showed that the proposed MF-ESMD denoising method had a more powerful denoising ability for a nonlinear and nonstationary signal than the morphological filter method and the MF-EEMD method.

Site Description and Data Acquisition.
In this study, to further validate the accuracy of the proposed MF-ESMD denoising method for the obtained dynamic deflection by ground-based microwave interferometry, the Fengbei Bridge was selected as an experimental bridge. e Fengbei Bridge is one of the most important transportation hubs in Beijing, and it is located in the West Fourth Ring of Beijing, China, as shown in Figure 6(a). Moreover, there is a road crossing under the bridge, as shown in Figure 6(a). An Imaging by Interferometric Survey (IBIS-S) instrument, a typical system based on microwave interferometry, was located on one side of the bridge without passive radar reflectors attached to the bridge, as shown in Figure 6(b). erefore, the instrument inevitably generated noise in the collected bridge dynamic deflection due to the surrounding environment, ground-motion, complex traffic activities, artificial operation, and the equipment itself. e IBIS-S instrument consisted of a radar unit, a control personal computer, a power supply unit, and a tripod. In typical measurement conditions, the sampling rate is up to 200 Hz, the maximal detection distance is up to 1 km, the range resolution is up to 0.50 m, and the displacement measurement accuracy is up to 0.01 mm. In this study, the angle of the altitude of the radar unit was set to 30 so that the two antennas on the radar unit could be aligned to the midspan point of the central span of the bridge. e range resolution was 0.5 m, the sampling frequency was 200 Hz, and the duration of data acquisition was 60 s. e dynamic deflection obtained for the Fengbei Bridge with IBIS-S is shown in Figure 7.  Shock and Vibration

Result Analysis and Discussion.
e original signal, the dynamic deflection of the Fengbei Bridge with IBIS-S, was first decomposed into 9 IMFs and a residual R using the ESMD method, as shown in Figure 8. As described in Section 3, the simulated signal was decomposed into 8 IMFs, which was less than the number of the decomposed IMFs of the dynamic deflection of the Fengbei Bridge with IBIS-S. e reason is that the dynamic deflection signal is more complicated than the simulated signal, due to the influences of surrounding environment, wind thrust, ground motion, complex traffic activities, human operations, and the equipment itself. en, Spearman's rho of each decomposed IMF was calculated, as shown in Table 3. According to Spearman's rho value, Spearman's rho of IMF 3 was 0.0380, which was the minimum value for the first time.
erefore, the first three IMF components, IMF 1 to IMF 3 , were regarded as the noisedominant IMFs, which may have been caused by the surrounding environment, complex traffic activities, artificial operation, and the equipment itself. en, the remaining IMF components were regarded as useful subsignals to be reconstructed as a new signal along with the residual R.
For the reconstructed signal, a mathematical morphology filter was used to perform the reduction of the detailed noise with the parameters g 1 � 1 1 1 and g 2 � 1 1 1 1 1 in equation (6). e resulting denoised signals using the MF-ESMD method and the original signal are shown in Figure 9(a). Figures 9(b) and 9(c) show the comparison results between the original signal and the denoised signal using the morphological filter and the MF-EEMD methods [27,29]. e inspection of the curves from this figure highlighted the following aspects: (1) compared with the original signal (the blue curve in Figure 9(a)), the denoised signal (the red curve in Figure 9(a)) using the MF-ESMD method was smoother and more stable, which indicated that the obvious noise information of the original signal was eliminated very well. (2) e denoised signal obtained by the MF-ESMD method was  smoother at the peaks and troughs of the curve than the denoised signal obtained by the morphological filter method and the MF-EEMD method. is indicated that the proposed MF-ESMD method had a better ability for the noise reduction in the peaks and troughs of the signal than the morphological filter method and the MF-EEMD method.
In this study, due to the complicated affecting factors, there were many kinds of noise in the obtained time series  displacement. It was difficult to evaluate the quality of signal denoising using the indexes of SNR and RMSE. erefore, two new evaluation indexes, RVR and NRR, were used to evaluate the quality of signal denoising. e RVR was the ratio of the variance root of the difference number of the signal after denoising to the variance root of the difference of the original signal, which could reflect the smoothness of the noise signal [7]. e smaller the RVR, the better the denoised signal.
e NRR could reflect the ability to suppress interference and improve the SNR [20]. In contrast, the bigger the NRR, the better the denoised signal: where s(n) is the n th element of the original signal, s(n) is the n th element of the denoised signal, N is the number of elements of the signal, σ 2 1 represents the variance of the measured original signal, and σ 2 2 denotes the variance of the denoised signals. Table 4 shows the values of RVR and NRR for the denoised signals obtained by the morphological filter, the MF-EEMD method, and the proposed MF-ESMD method. e inspection from this table highlighted the following aspects: (1) the RVR was 0.0524 for the proposed MF-ESMD denoising method, which was smaller and closer to zero than the values for the morphological filter method and the MF-EEMD method. (2) e NRR was 3.522 for the proposed MF-ESMD denoising method, which was better than the values of the morphological filter method and the MF-EEMD method. However, the NRR was −4.9465 for the morphological filter method, which was a negative number. As shown in Figure 9(b), when the sampling time was approximately 40 s, there was strong noise in the original signal, which blurred the characteristics of the useful signal using the mathematical morphology. erefore, the morphological filter method could not distinguish the useful signal and noise, and it might treat the useful subsignals as the noise for the instant stronger noise. For the MF-EEMD method and the proposed MF-ESMD method, most of the obvious noises were reduced when the useful subsignal was reconstructed. ere was only a small amount of detailed noise in the reconstructed signal, and the mathematical morphology filter was used to perform a further reduction of the detailed noise in the reconstructed signal. e results showed that the proposed MF-ESMD method had a good denoising ability, which not only efficiently reduced the      effects of noise for dynamic deflection signals but could also retain useful information.

Conclusions
In recent years, ground-based microwave interferometry has been widely used to obtain the dynamic deflection measurements of bridges. Due to the influences of the surrounding environment and the equipment itself, there has been an inevitable increase in the noise information of the obtained dynamic deflection displacement. In this study, in order to improve the accuracy of the bridge dynamic deflection signals obtained by ground-based microwave interferometry, the MF-ESMD denoising method was proposed in order to remove the noise in the original signal.
rough simulated experiments and real data experiments, it was found that the proposed method had the following advantages.
(1) e simulated experiment was performed in order to validate the feasibility and accuracy of the proposed MF-ESMD denoising method for a nonlinear and nonstationary simulated signal. e RMSE was 0.0361 mm and the SNR was 30.6619 dB for the proposed MF-ESMD denoising method. ese values were much better than those of the morphological filter method and the MF-EEMD method. Moreover, the SNR of the proposed MF-ESMD denoising method increased by 53.3% compared with the SNR of the simulation signal. e results showed that the proposed MF-ESMD denoising method was useful for the signal denoising of the nonlinear and nonstationary signals.
(2) Compared with the morphological filter method and the MF-EEMD method, the MF-ESMD denoising method had better values for the indexes of RVR and NRR for the dynamic deflection displacement obtained using ground-based microwave interferometry. e results verified the feasibility and accuracy of the proposed MF-ESMD method for dynamic deflection displacement obtained using groundbased microwave interferometry. (3) is MF-ESMD method integrated the ESMD method and Spearman's rho to remove the noisedominant IMFs from the decomposed IMFs for the bridge. A morphology filter was used to perform the reduction of the detailed noise in the reconstructed signal. Compared with the morphological filter method and the MF-EEMD method, the proposed MF-ESMD method had a good denoising ability, which not only efficiently reduced the effects of noise for the dynamic deflection signals but could also retain the useful information.

Data Availability
All the data included in this study are available upon request by contact with the corresponding author.

Conflicts of Interest
e authors declare no conflicts of interest.