Abstract

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 ESMD method. 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. The results of both the simulated and on-site experiments showed that the proposed MF-ESMD method had a powerful signal denoising ability.

1. 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. Therefore, 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 [25]. 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]. Therefore, it is necessary to ascertain high-precision 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 [710]. 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 ground-based 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. These 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]. Therefore, it is of great importance to reduce the influence of noise in the bridge dynamic deflection obtained by ground-based microwave interferometry.

At present, wavelet transform is one of the widely used methods for signal denoising with a good time-frequency localization characteristic [14]. This 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 [1618]. 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. This 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]. The decomposed components contain the local characteristic signals for different time scales of the source signal [21]. The 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, 2428]. This 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, 2430]. 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]. The 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 [3135]. 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 [3639]. This technology has been widely applied in the analysis of signals such as vibration signals and pulse wave signals [4042]. Morphological filtering based on the theory of mathematical morphology is a new nonlinear filtering method for eliminating noise while protecting signal details. The theory of mathematical morphology is an integrated theory of information theory, statistical estimation theory, adaptive theory, random noise theory, and other theories [3639]. 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. Thus, 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 mode-mixing 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.

Therefore, 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. The ESMD method was used to decompose the original obtained bridge dynamic deflection into a series of IMFs, and the noise-dominant 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. The 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.

2. Methods

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 noise-dominant 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.

2.1. Reduction of Noise-Dominant IMFs

The 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]. Therefore, 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 , 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 were found, and the midpoints of the two adjacent extreme points were defined as . The midpoints of the left and right borders were set to and using a linear interpolation method. Interpolation curves (named ) were curved using cubic spline interpolation based on the points in , and their average curve was further calculated.Step 2: was calculated by repeating Step 1 to obtain until was satisfied or the sifting times reached the preset maximum value of 30. In this study, the permitted error was set to 0.001 , where was the standard deviation of the original data relative to the total mean of .Step 3: Steps 1 and 2 were repeated for to obtain , …, and a residual R, which is an optimal adaptive global mean (AGM) curve with no more than a certain number of extreme points. Then, the variance of was calculated.Step 4: The sifting times was changed on a finite integer interval . Then, Steps 1–3 were repeated to obtain the optimal sifting time , according to the minimum variance on .Step 5: The optimal sifting times value and the permitted error were used to perform the optimal signal decomposition for the dynamic deflection signal 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]. The bigger the coefficient value, the more similar it was to the original signal [30, 44]. Through the coefficient value, noise-dominant IMFs could be removed. Spearman’s rho between the decomposed and the original signal could be calculated using equation (1):where is the IMF and is the number of elements of the original dynamic deflection signal.

2.2. Reduction of the Residual Noise

After the reduction of the noise-dominant IMFs, the other decomposed IMFs were reconstructed as a new signal. There 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 [1418]. was defined as a signal, and was defined as a structuring element. The dilation operation was defined as shown in equation (2), and the erosion operation was defined as shown in equation (3). Based on the dilation and erosion operations, the opening operation was defined as shown in equation (4), and the closing operation was defined as shown in equation (5):where , , , is the length of , and is the width of .

However, the opening operation could only change the maxima and maintain the shape of the minima for the reconstructed signal. The closing operation could only change the minima and maintain the shape of the maxima for the reconstructed signal. Therefore, 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 is the final denoised signal, is the reconstructed signal, and is the length of signal . and are the two structuring elements with different sizes, of which the lengths were less than the length of the reconstructed signal . In addition, the structuring element was an important factor affecting the results of the signal denoising. The most frequently used structuring elements include linear structuring elements, triangular structuring elements, and structuring elements [3639]. The filtering results were affected by the length and shape of the structuring elements. Therefore, in the actual processing of the signal, the structuring elements were designed according to the characteristics of the analysis signal and the computational complexity [3639]. 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].

3. Simulated Experiments and Analysis

The dynamic deflection signal of bridges obtained by ground-based microwave interferometry is a complicated nonstationary signal mixed with some noise. Therefore, in this study, to validate the feasibility and accuracy of the proposed MF-ESMD denoising method, a nonlinear and nonstationary simulated signal was generated that consisted of three basic functions with the different frequencies and a white Gaussian noise signal with an SNR (signal-to-noise ratio) of 20 dB to enhance noise information in the simulated signal , which is shown as follows:where , , , and is a white Gaussian noise signal. The waveforms of the simulation signal and its component signals , , and are shown in Figure 2.

As shown in Figure 3, the simulated signal was decomposed into eight IMFs with the corresponding meaningful inherent natural frequency from high to low frequency and a residual using the ESMD method. Then Spearman’s rho of each decomposed IMFs was calculated, as shown in Table 1. The value of Spearman’s rho from the first decomposed IMF1 was gradually reduced. Spearman’s rho of IMF3 was 0.0291, which was the minimum value for the first time. This indicated that the decomposed IMF1 to IMF3 had less similar feathers than the original signal . Therefore, the first three IMF components, IMF1 to IMF3, 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. The remaining IMF components were regarded as useful subsignals to be reconstructed as a new signal along with the residual .

For the reconstructed signal , most of the obvious noises were reduced. However, there was a small amount of residual noise in the reconstructed signal . Therefore, a morphology filter was used to perform the reduction of the residual noise in the reconstructed signal . In the morphology filter, and were denoted as shown in equation (6). Then, the final denoised signal could be obtained using the morphology filter. Figure 4 shows the compared results of the curves between the original signal and the denoised signal . 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. The comparison results of the curves between the original signal and the denoised signal using the morphological filter and the MF-EEMD methods are shown in Figures 5(a) and 5(d). The 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. The 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. The 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. The 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. The 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]. The 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. The RMSE could indicate the difference between the original signal and the denoised signal [7, 20]. The bigger the SNR value, the smaller the RMSE value and the better the effect of noise reduction:where is the element of the original signal, is the element of the denoised signal, and N is the number of elements of the signal.

The comparison results are shown in Table 2. The 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) The SNR of the simulation signal was 20 dB. The SNR was 26.3066 dB for the denoised signal using the morphological filter method, which increased by 31.53% compared to the simulation signal . The SNR was 26.4298 dB for the denoised signal using the MF-EEMD method, which increased by 32.15% compared to the simulation signal . The 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 . The 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.

4. On-Site Experiment and Analysis

4.1. 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. The 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). Therefore, 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. The 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. The range resolution was 0.5 m, the sampling frequency was 200 Hz, and the duration of data acquisition was 60 s. The dynamic deflection obtained for the Fengbei Bridge with IBIS-S is shown in Figure 7.

4.2. Result Analysis and Discussion

The 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. The 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. Then, Spearman’s rho of each decomposed IMF was calculated, as shown in Table 3. According to Spearman’s rho value, Spearman’s rho of IMF3 was 0.0380, which was the minimum value for the first time. Therefore, the first three IMF components, IMF1 to IMF3, were regarded as the noise-dominant IMFs, which may have been caused by the surrounding environment, complex traffic activities, artificial operation, and the equipment itself. Then, 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 and in equation (6). The 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]. The 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) The 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. This 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. Therefore, two new evaluation indexes, RVR and NRR, were used to evaluate the quality of signal denoising. The 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]. The smaller the RVR, the better the denoised signal. The 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 is the element of the original signal, is the element of the denoised signal, N is the number of elements of the signal, represents the variance of the measured original signal, and 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. The 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) The 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. Therefore, 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. There 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. The 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.

5. 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. Through simulated experiments and real data experiments, it was found that the proposed method had the following advantages.(1)The 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. The RMSE was 0.0361 mm and the SNR was 30.6619 dB for the proposed MF-ESMD denoising method. These 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. The 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. The results verified the feasibility and accuracy of the proposed MF-ESMD method for dynamic deflection displacement obtained using ground-based microwave interferometry.(3)This MF-ESMD method integrated the ESMD method and Spearman’s rho to remove the noise-dominant 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

The authors declare no conflicts of interest.

Authors’ Contributions

Xianglei Liu and Mengzhuo Jiang conceptualized the study. Mengzhuo Jiang and Ziqi Liu developed methodology. Mengzhuo Jiang and Hui Wang performed formal analysis. Ziqi Liu and Hui Wang were involved in data curation. Xianglei Liu performed English editing. Xianglei Liu was responsible for financial support.

Acknowledgments

This study was sponsored by the National Natural Science Foundation of China (grant no. 41871367), the Ministry of Science and Technology of the People’s Republic of China (grant no. 2018YFE0206100), the Importation and Development of High-Caliber Talents Project of Beijing Municipal Institutions (grant no. CIT&TCD201704053), the Science and Technology Project of Ministry of Housing and Urban-Rural Development of the People’s Republic of China (grant no. 2017-K4-002), the Scientific Research Project of Beijing Educational Committee (grant no. KM201910016007), the Major Projects of Beijing Advanced Innovation Center for Future Urban Design (grant no. UDC2018031321), and the BUCEA Post Graduate Innovation Project.