Nothing Special   »   [go: up one dir, main page]

Next Article in Journal
Multisensor Data Fusion by Means of Voxelization: Application to a Construction Element of Historic Heritage
Next Article in Special Issue
Built-Up Area Extraction from GF-3 SAR Data Based on a Dual-Attention Transformer Model
Previous Article in Journal
Tracking Individual Scots Pine (Pinus sylvestris L.) Height Growth Using Multi-Temporal ALS Data from North-Eastern Poland
Previous Article in Special Issue
Object-Oriented Unsupervised Classification of PolSAR Images Based on Image Block
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

Removing InSAR Topography-Dependent Atmospheric Effect Based on Deep Learning

1
College of Earth Science, Chengdu University of Technology, Chengdu 610059, China
2
State Key Laboratory of Geological Disaster Prevention and Geological Environmental Protection, Chengdu University of Technology, Chengdu 610059, China
3
College of Computer Science and Cyber Security (Oxford Brookes College), Chengdu University of Technology, Chengdu 610059, China
4
GeoAI, Smarter Map and LiDAR Lab, Faculty of Geosciences and Environmental Engineering, Southwest Jiaotong University, Chengdu 610097, China
5
Department of Geotechnics and Geomatics, Saveetha School of Engineering, Saveetha Institute of Medical and Technical Sciences, Chennai 602105, India
6
Huadong Engineering Corporation Limited, Hangzhou 311122, China
7
College of Geological Engineering and Geomatics, Chang’an University, Xi’an 710064, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(17), 4171; https://doi.org/10.3390/rs14174171
Submission received: 25 July 2022 / Revised: 17 August 2022 / Accepted: 22 August 2022 / Published: 25 August 2022
(This article belongs to the Special Issue SAR in Big Data Era II)
Graphical abstract
">
Figure 1
<p>(<b>a</b>) Topographic map of the study area; (<b>b</b>) location of the study area; (<b>c</b>,<b>d</b>) images of the study area from Google Earth.</p> ">
Figure 2
<p>Flowchart of the proposed method.</p> ">
Figure 3
<p>(<b>a</b>–<b>f</b>) Original interferometric pair; (<b>g</b>–<b>l</b>) simulated phase interferometric pair by MLP neural network model; (<b>m</b>–<b>r</b>) interferometric pairs after atmospheric correction by MLP neural network model.</p> ">
Figure 4
<p>The correction effects of three topography-dependent atmospheric delay correction methods on InSAR interferograms.</p> ">
Figure 5
<p>Changes in phase standard deviation before and after atmospheric correction by three methods.</p> ">
Figure 6
<p>Interferometric map 20210604–20210610 phase and elevation linear fit analysis:(<b>a</b>) the original phase with elevation; (<b>b</b>) phase corrected by linear model with elevation; (<b>c</b>) phase corrected by GACOS with elevation; (<b>d</b>) phase corrected by MLP neural network model with elevation.</p> ">
Figure 7
<p>(<b>a</b>) Original displacement results; (<b>b</b>) displacement results after linear model correction; (<b>c</b>) displacement results after GACOS correction; (<b>d</b>) displacement results after correction of MLP neural network model.</p> ">
Figure 8
<p>(<b>a</b>) The MLP neural network atmosphere corrects the displacement results; (<b>b</b>) remote interpretation of displacement.</p> ">
Figure 9
<p>(<b>a</b>) Original displacement results; (<b>b</b>) displacement results after linear model correction; (<b>c</b>) displacement results after GACOS correction; (<b>d</b>) displacement results after correction of MLP neural network model.</p> ">
Figure 10
<p>(<b>a</b>) The MLP neural network atmosphere corrects the displacement results; (<b>b</b>) remote interpretation of displacement.</p> ">
Versions Notes

Abstract

:
Atmospheric effects are among the primary error sources affecting the accuracy of interferometric synthetic aperture radar (InSAR). The topography-dependent atmospheric effect is particularly noteworthy in reservoir areas for landslide monitoring utilizing InSAR, which must be effectively corrected to complete the InSAR high-accuracy measurement. This paper proposed a topography-dependent atmospheric correction method based on the Multi-Layer Perceptron (MLP) neural network model combined with topography and spatial data information. We used this proposed approach for the atmospheric correction of the interferometric pairs of Sentinel-1 images in the Baihetan dam. We contrasted the outcomes with those obtained using the generic atmospheric correction online service for InSAR (GACOS) correction and the traditional linear model correction. The results indicated that the MLP neural network model correction reduced the phase standard deviation of the Sentinel-1 interferogram by an average of 64% and nearly eliminated the phase-elevation correlation. Both comparisons outperformed the GACOS correction and the linear model correction. Through two real-world examples, we demonstrated how slopes with displacements, which were previously obscured by a significant topography-dependent atmospheric delay, could be successfully and clearly identified in the interferograms following the correction by the MLP neural network. The topography-dependent atmosphere can be better corrected using the MLP neural network model suggested in this paper. Unlike the previous model, this proposed approach could be adjusted to fit each interferogram, regardless of how much of the topography-dependent atmosphere was present. In order to improve the effectiveness of DInSAR and time-series InSAR solutions, it can be applied immediately to the interferogram to retrieve the effective displacement information that cannot be identified before the correction.

Graphical Abstract">

Graphical Abstract

1. Introduction

InSAR technology is a geodetic measurement technology [1] based on spaceborne Earth observation, which has been developed rapidly in the past three decades. Satellite radar remote sensing has prominent strengths, such as wide-range coverage, high accuracy, and all-day and all-weather applicability [2]. At present, although InSAR technology has been widely used in landslide identification and monitoring [3,4,5,6,7,8], seismic research [9,10,11,12], ground subsidence monitoring [13,14,15], volcano monitoring [16,17,18], glacier movement [19,20] and polar research [21], atmospheric delay error is still one of the major errors limiting the high-accuracy application of InSAR technology [22,23]. When the atmospheric relative humidity changes by 20% in time and space, the atmospheric effect can contribute to a deformation error of 10–14 cm, which is not allowed for the high-accuracy (millimetre level) deformation detection of InSAR [24].
Many scholars have proposed various methods to reduce the influences and errors caused by atmospheric effects. The correction methods are divided into two categories: atmospheric correction with external data or by the model without external data. Atmospheric corrections without external data include atmospheric delay field models [25,26], linear superposition methods [27], topography-dependent turbulence models (referred to as GTTM) [22], and mesoscale meteorological models [28] and linear models [29,30]. The second is the correction method based on external data, which mainly includes Global Navigation Satellite System (GNSS) water vapour data and ground-based meteorological observations correction method [25,31,32], Moderate-resolution Imaging Spectroradiometer (MODIS) and Medium Resolution Imaging Spectrometer (MERIS) water vapour product data correction method [33,34]. Recently, the global medium-resolution numerical atmospheric reanalysis products from the European Centre for Medium-Range Weather Forecasts (ECMWF) have been widely used for the tropospheric delay estimation and correction of individual interferograms and time-series InSAR [35,36,37]. Based on an iterative tropospheric decomposition model that combines continuous GNSS tropospheric estimates with high-resolution ECMWF data, Yu et al. [38] created a generic atmospheric correction online service (GACOS) product with global coverage, all-day, and all-weather data.
The above-generalized method is more appropriate for a broad area, while the study area presents a topography-dependent atmospheric effects with vertical distribution. Among the topography-dependent atmospheric effect correction research, Beauducel et al. (2000) [39] obtained better correction results for topography-dependent atmospheric effects using atmospheric vertical stratification coefficients. Li et al. (2004) [40] proposed a model that integrated continuous GPS (CGPS) and ground meteorological observation data to simulate and correct topography-dependent atmospheric effects in InSAR measurements. Li et al. (2006) [22] established an atmospheric correction model (abbreviated as GTTM) using GPS data, considering the vertical distribution of tropospheric water vapour in terms of topography, which reduces water vapor’s influence on interferograms. External data are often restricted by spatial resolution. Thus, the correction of the topography-dependent atmosphere at a small scale without the help of external data still warrants further study.
This paper proposes a topography-dependent atmospheric correction method based on the MLP neural network model integrated with topographic and spatial data (described in Section 3). The method is applied to the Baihetan dam, comparing the correction results with the traditional linear model and GACOS correction model to evaluate the correction effects of the three models (Section 4). Meanwhile, the method was practically applied to two pairs of interferograms for atmospheric correction in the Baihetan reservoir area, where the topography-dependent atmospheric effects are severe. Potential landslides that could not be detected before the atmospheric correction was extracted using the method (Section 5), which is of practical significance to improve the accuracy of monitoring disaster sites through InSAR technology.

2. Study Area and Data Source

2.1. Study Area

The study area in this research was Baihetan dam, which is located in Ningnan County, Sichuan Province (27.04°N, 102.91°E) and Qiaojia County, Yunnan Province (26.78°N, 103.12°E). It is located at the southeast edge of the Qinghai-Tibet Plateau, with elevations between 1000 and 3000 m, and belongs to the high mountain and plateau geomorphological unit and the Transverse Mountain System in southwest Sichuan and northeast Yunnan. Due to the significant elevation difference, the study area has distinct vertical stereoclimatic features covering all climate types from southern subtropical to cold temperate zones. The area is mainly composed of river erosion landforms, tectonic landforms and ice erosion landforms, with deep river valleys and severe weathering and denudation [41]. Figure 1b shows that the reservoir area is located in the lower reaches of the Jinsha River, and the terrain is high in the east and low in the west. Moreover, the slopes along both sides of the river are steep, the mountains are high, the valleys are deep, and the river valley landforms mostly comprise “V”-shaped development, which is a typical high mountain deep canyon area. Due to impoundment, the water elevation rises dramatically, leading to landslides, mudflows, collapses and other geological hazards in the study area, so a more accurate displacement monitoring of the study area is especially crucial [42]. The real topography of the study area is shown in Figure 1c,d.

2.2. Data Source

A total of 17 Sentinel-1 images, including 10 images of the descending orbit and 7 images of the ascending orbit, were used in this study. The Sentinel-1 satellite incorporates four modes of operation: stripmap mode (SM), interferometric wide swath mode (IW), extra-wide swath mode (EW), and wave mode (WV). In this paper, we selected single-look complex (SLC) images of VV polarization in IW mode, spanning from 21 April 2021 to 4 July 2021, all after the impoundment of the Baihetan Reservoir. The Sentinel-1 single satellite uses topographic observation in progressive scans (TOPS) mode to acquire images on a global scale every 12 days [43]; that is, the satellite has a revisit period of 12 days. Other parameters, such as the incidence angle and resolution of the Sentinel-1 satellite, are shown in Table 1.

3. Methodology

3.1. DInSAR Processing

Differential Interferometry Synthetic Aperture Radar (DInSAR) monitor surface displacement was performed using the phase information of 2 SAR images acquired from the SAR sensor at different times in the same area [44]. As shown in Figure 2, SAR images covering the same area were acquired and co-registered with the precise orbit data to form interferograms. Shuttle radar topography mission digital elevation model (SRTM DEM) data, with a 30 m resolution, were used to remove the topographic and flat effect from the interferogram. The Goldstein adaptive filtering method [45] was used for filtering to remove the noise. Finally, a displacement map was acquired after phase unwrapping using the minimum cost flow (MCF) method [46,47].

3.2. MLP-Based Topography-Dependent Atmospheric Correction Method

In the deep learning community, we have witnessed the rapid development of neural networks. The architecture of neural networks developed from the MLP (Multi-Layer Perceptron), to the CNN (Convolutional Neural Network), to the Transformer, and back to the MLP [48]. In this study, we proposed an MLP with ten layers, referred to as FatMLP, for removing the atmospheric phase in InSAR. FatMLP was designed to learn high-order non-linear interactions of the digital elevation model (DEM), longitude and latitude.
The architecture of FatMLP is shown in Figure 2, which consists of ten layers. In FatMLP, an MLP layer was used to learn the linear correlations of the input variables. We used a wider MLP to capture more linear correlations. We designed deeper layers for higher-order nonlinear interactions of the input variables. The input layer included a fully connected layer with 3 neurons, which represent the DEM, longitude and latitude, respectively. Each pixel of the DEM was treated as a sample. First, the input data were lifted into a high-dimensional space, in which the representation of the original data was limited [49]. To this end, the first two hidden layers were fully connected layers with 4096 neurons. Such a high-dimensional space is useful for learning complex interactions of the input variables. Second, the intermediate sparse data were gradually compressed. The third and fourth hidden layers were two fully connected layers with 2048 neurons. The fifth and sixth layers were two fully connected layers with 1024 neurons. The seventh and eighth layers were two fully connected layers with 512 neurons. Finally, the output layer was fully connected with 1 neuron, representing the atmospheric phase. For each hidden layer, the activation function was the Rectified Linear Unit (ReLU), which was used to learn non-linear interactions of the input variables. The depth of FatMLP was greater than conventional MLPs. Deeper neural networks are useful for learning high-order interactions of the input variables. Equations (1)–(5) show the implementation details of the proposed neural network, where FC represents a fully connected layer. The following hyperparameters and configurations are displayed: 9192 is the batch size; the Adaptive Moment Estimation (Adam) is the optimizer, which is one of the most widely used optimization algorithms in deep learning, with the advantages of fewer hyperparameters and good performance; the rate of learning is 0.001; and the proposed neural network was trained for 50 epochs.
X 4096 = R e l u F C R e l u F C x d e m , x l , x l
X 2048 = R e l u F C R e l u F C X 4096
X 1024 = R e l u F C R e l u F C X 2048
X 512 = R e l u F C R e l u F C X 1024
            Y p = F C R e l u F C X 512

3.3. Linear Model Correction

The distribution of atmospheric water vapour is the primary cause of tropospheric delays, with impacts of atmospheric vertical stratification introducing topography-dependent delayed phases and turbulence introducing local delayed phases. In contrast to delays caused by turbulent mixing, vertical stratification effects are regarded as static stratification in a particular location and time scale. In the case of changes in surface elevation, the vertical stratification tropospheric delay increases as the water vapour concentration and pressure to temperature ratio rise [50]. The concentration of water vapour declines exponentially with altitude, and the theoretical delay curve is an exponential function of altitude [29]. Suppose the interferogram is majorly affected by the static tropospheric delay. In this case, the delay signal is the difference between the corresponding delay curves of the primary and secondary SAR images, and the resulting exponential function is a Taylor-expanded and second-order function. Higher-order terms are neglected, and a linear relationship can be obtained as follows [51]:
Δ Ø = b + K h
where b is a deviation term, h is the topography, Δ Ø is the phase, and K is the transfer function between the topography and the phase. The correction of tropospheric delay using a linear model is applicable to most weather conditions, except for special cases, such as reverse or non-monotonic troposphere. Since the elevation–phase relationship may have a local trade-off with the elevation displacement relationship, the transfer function K is best determined through global linear regression [52].

3.4. GACOS Correction Method

GACOS, an InSAR atmospheric correction model developed by Zhenhong Li’ (2006) [22] at Newcastle University, United Kingdom, is a numerical atmospheric product based on a high-resolution European Centre for Medium-Range Weather Forecasts grid of up to 0.1° × 0.1°. The model utilises an iterative tropospheric decomposition (ITD) model to separate the topography-dependent vertically stratified delay component and the turbulent mixing delay component from the total tropospheric delay in order to generate a high spatially resolved (default 90 m) zenith total delay map for InSAR atmospheric delay correction [53]. InSAR users can obtain ZTD products for the corresponding time period by submitting the location of the study area and the time of SAR image acquisition through (http://www.gacos.net) (accessed on 5 January 2022).

4. Results and Discussion

4.1. Analysis of Correction Effect

As shown by the analysis in Figure 3, the interferograms a–f suffered from more serious atmospheric delay effects associated with elevation, and the atmospheric delay phase was related to the elevation change. The atmospheric simulated interferogram (i.e., MLP-modelled atmospheric correction that needs to be removed from the original interferograms) is shown in Figure 3g,h. The topography-dependent atmospheric delay effects in Figure 3m–r were reduced to different degrees after correction by the MLP neural network method, and the phases of the interferograms were smoother. Overall, the effect of atmospheric delay phase removal was more obvious after correction using the MLP neural network method. Meanwhile, the method is not required to build models, does not introduce any external data, has certain applicability for topography-dependent atmospheric corrections, and the atmospheric correction effect is relatively significant.
Table 2 shows the statistics of the change in the phase standard deviation of the original deconvolution of each interferometric pair. The analysis revealed that the phase standard deviation of all interferometric pairs after correction by the MLP neural network model was significantly reduced compared to the phase standard deviation before correction, with the 20210604–20210610 interferometric pair having the best correction effect and the highest change rate of up to 84%.

4.2. Linear Model and GACOS Correction Method

Figure 4 shows the corrective effect of the Sentinel-1 (21 April 2021 to 4 July 2021) interferogram with the three correction methods. As seen in Figure 4, there was a strong association between the interferogram phase spatial distribution and the topographic relief of the study area, indicating that the delayed vertical stratification of the atmosphere predominates in terms of the overall trend. After correction using the conventional linear model, GACOS, and MLP neural network model, the combination of the second and third rows of Figure 4 demonstrated that the topography-dependent atmospheric delay in the linear model and GACOS correction findings were attenuated to varying degrees. The analysis of the fourth row of Figure 4 revealed, however, that after the correction of the MLP neural network model, the simulated phase of the atmosphere was visibly lowered in comparison to the phase value before correction, and the phase value in the majority of places was close to 0. Thus, the MLP neural network correction had a greater impact than the linear model and GACOS.
To optimize the use of different interferometric SAR atmospheric correction methods and to avoid potential uncertainty, it is necessary to establish appropriate statistical metrics and evaluate the correction performance [54]. The phase standard deviation was applied to measure the effect of atmospheric delay estimation and correction. In particular, the reduction in the interferometric phase was used to demonstrate the reduction in the corrected atmospheric delay phase [36,37]. Using the three aforementioned techniques, atmospheric effects were removed from 27 interferograms, and the standard deviation of the interferometric phase before and after correction was calculated to evaluate the overall performance of the correction results (Figure 5).
Figure 5 shows the magnitude of the reduction in phase standard deviation of the 27 interferometric pairs in the Baihetan reservoir region after atmospheric correction using linear model correction, external data GACOS and MLP neural network, respectively. The analysis revealed that after linear correction and GACOS correction, the phase standard deviation of some interference pairs decreased, which was not significant, while the phase standard deviation of some interference pairs increased instead. Nevertheless, with MLP neural network correction, the phase standard deviation of all 27 interference pairs was reduced, and the reduction was more significant.
Table 3 shows that after correction using the linear model, the interferometric pairs with reduced phase standard deviation accounted for 59% of the 27 interferometric pairs, which is the lowest among the three methods. The GACOS atmospheric correction had 88% of the interferometric pairs with reduced phase standard deviation of the unwinding. The percentage of interferograms with reduced standard deviation was larger than that of the linear model. However, the percentage of interferometric pairs with reduced phase standard deviation after correction using the MLP neural network was 100%, which was significantly higher than that of the linear model and GACOS. Hence, for the topography-dependent atmospheric delay in the Baihetan reservoir area, the atmospheric correction results of the MLP neural network were both superior to those of the linear model and the external data GACOS, and the GACOS was superior to the linear model. Compared with the atmospheric delay correction accuracy, the phase standard deviation of all interferometric pairs was reduced by 11.1% on average after using the linear model for atmospheric delay correction, which was the lowest among the three methods. However, after atmospheric delay correction using GACOS, the average reduction in the standard deviation of the unwrapping phase was 17.4%, which was higher than the linear model and located in the second rank. Finally, after the correction of atmospheric delay using the MLP neural network, the average reduction in phase standard deviation was approximately 64%, which was much higher than the first two correction methods, and the correction effect was more obvious for the topography-dependent atmospheric delay in the Baihetan reservoir area.

4.3. Comparative Analysis of Topographic-Phase Correlations

Jolivet et al. (2014) [36] estimated the correlation coefficient between the interferometric phase and elevation as an alternative to the ratio between turbulent and vertical stratified delays. Here, we used the correlation coefficient to describe the elevation dependence in the interferogram. The smaller difference from zero indicates the lower proportion of signals related to the terrain in the interferometric phases. If the absolute value of the correlation decreases after atmospheric correction, it represents the mitigation of the vertical stratified components [36].
As demonstrated by the example of the analysis of the correlation atmospheric delay correction for the elevation of the Baihetan reservoir area in Figure 6, the correlation between phase and elevation for the original interferograms 20210604–20210610 (Figure 6a) was −0.092. After the linear model and GACOS correction, the correlation between phase and elevation was reduced to −0.082 and −0.002, respectively, as shown in Figure 6b,c. However, there was still some correlation between phase and elevation. Figure 6d illustrates that after rectification using the MLP neural network model, the correlation was reduced to −0.000036, which is quite close to 0. The correlation between phase and elevation was rather weak at this stage.

5. Displacement Identification after MLP Correction

When InSAR technology is applied to early displacement monitoring, some minor displacement information is prone to being obscured due to the atmospheric delay effect in the study area, which poses an invisible threat to human life and property. Therefore, it is essential to effectively identify the deformers in some areas through atmospheric correction. In this paper, the interferometric pair 20210517–20210604 with severe atmospheric delay in the Baihetan reservoir area was atmospherically corrected using the above three methods. The results of the atmospheric correction were processed for phase turn displacement, and it was found that there were signs of displacement in the local area. The specific displacement results are shown in Figure 7. The corrected displacement results of the original displacement results, linear model, GACOS, and MLP neural network model are represented in Figure 7a–d, respectively.
As seen in Figure 7, the original displacement results map can be seen in the slope area without signs of displacement, surrounded by atmospheric phase error. It was seriously impacted, and was unable to determine the exact displacement location and displacement magnitude accurately. Combined with Figure 7b,c, the linear model and GACOS correction were ineffective, and the displacement region was not easily identified. However, from Figure 7d, it is clear that the MLP neural network model had a remarkably effective correction effect. In addition to the obvious displacement in the slope ranges, the displacement volume in the non-displacement area was approximately 0mm, and it was able to analyse the displacement area more accurately.
The displacement area, located near Sanjia Village, Yogu Country, Huidong County, Liangshan Yi Autonomous Prefecture, Sichuan Province, with a central longitude of 103°1′19″ and a central latitude of 26°36′59″, was further analysed. This slope elevation was approximately 772–1028 m, with a differential height of 256 m and a relatively steep slope. Combined with Figure 8a,b, the strong displacement area of this slope was more obvious (as shown by the blue dotted line), which was located in the median–lower part of the slope, with a maximum line-of-sight displacement of up to approximately 28 mm.
Meanwhile, the interferometric pairs 20210604–20210610 with severe atmospheric delay in the Baihetan reservoir area were corrected for atmospheric correction using the three methods mentioned above, and the results of the atmospheric correction were processed for phase turn displacement. From the specific displacement results of Figure 9, we found that the local area showed signs of displacement. The original displacement results, linear model, GACOS, and MLP neural network model corrected displacement results are represented in Figure 9a–d, respectively.
Figure 9a showed that the atmospheric delay effect seriously influenced the original displacement results, and the displacement location and displacement magnitude could not be accurately identified. The linear model correction was not effective, and the atmospheric phase remained serious after the correction (Figure 9b). From Figure 9c, it can be seen that the GACOS correction had a certain effect but failed to accurately identify the location of the landslide. Figure 9d reveals that the MLP neural network model had the optimal correction effect. Moreover, with the exception of the obvious displacement in the slope ranges, the displacement volume in the non-displacement area was approximately 0 mm, allowing for a more accurate analysis of the displacement area.
The further analysis of the displacement area showed that the area is located near Dapingzi village, Nagu Town, Huize County, Qujing City, Yunnan Province, with a central longitude and latitude of 103°4′24″ and 26°31′3″, respectively. The elevation of the slope is approximately 1160–1492 m, with a differential height of 332 m, and the slope is relatively steep. Based on the displacement monitoring results shown in Figure 10a, it is evident that the area of strong displacement in this region was more apparent (as indicated by the blue dotted line), with the maximum line-of-sight displacement reaching approximately 40 mm/year. Figure 10b shows the remote sensing interpretation of the displacement area with clear displacement signs. This active slope with displacement was also validated by the SBAS-InSAR result and our fieldwork.
Among the advantages of GACOS are its versatility and global availability. GACOS has advantages in dealing with long-wave atmospheric delay errors and may have some limitations in local applications with large undulating mountains. Therefore, atmospheric errors can be further and more effectively removed in local applications with the help of more refined models and methods (such as the deep learning methods proposed in this paper). It should be mentioned that in order to maintain a high coherence of interferometry, the temporal baseline range was within one month in this study. The effect of the temporal atmospheric water content variations from wet/dry seasons should be evaluated if interferograms with long temporal baselines were used.

6. Conclusions

This study proposed a topography-dependent atmospheric correction method based on the MLP neural network model integrated with topographic and spatial data. The method was used to correct the atmosphere of an interferometric pair composed of Sentinel-1 ascending and descending orbit SAR in 17 views of the Baihetan dam. Moreover, the correction results were compared with the traditional linear model and GACOS to evaluate the effect of correction and the rate of change of the phase standard deviation of the three models.
The results demonstrated that the MLP neural network method could effectively estimate and correct the topography-dependent atmospheric delay, and the corrected Sentinel-1 wide-field InSAR interferogram phase standard deviation mean reduction showed a 64% improvement compared to that of the traditional linear model and GACOS. Furthermore, through fitting the phase-elevation linear relationship, it was clear that the linear coefficient of the phase and elevation correlation after the MLP neural network model correction was almost 0, which means that the phase and elevation are relatively uncorrelated. Eventually, two practical applications confirmed that the MLP neural network model could successfully identify two landslides in Sanjia and Dapingzi villages in the Baihetan reservoir area. The MLP neural network model proposed in this paper is more credible for topography-dependent atmospheric correction results, which can achieve the high-accuracy application of InSAR technology in reservoir area displacement monitoring.

Author Contributions

Conceptualization, K.D. and X.T.; methodology, C.C. and J.C.; investigation, J.C.; resources, K.D.; software, X.T.; writing—original draft preparation, C.C.; writing—review and editing, K.D., S.P., M.W., X.S. and H.Z.; funding acquisition, Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (grant number 2021YFB3901403), the National Natural Science Foundation of China (grant numbers 41941019, 42072306, and 41801391), the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection Independent Research Project (SKLGP2020Z012) and the fellowship of China Postdoctoral Science Foundation (2020M673322; 2021M690024).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Poreh, D.; Pirasteh, S. InSAR observations and analysis of the Medicina Geodetic Observatory and CosmoSkyMed images. Nat. Hazards 2020, 103, 3145–3161. [Google Scholar] [CrossRef]
  2. Tofani, V.; Segoni, S.; Agostini, A.; Catani, F.; Casagli, N. Use of remote sensing for landslide studies in Europe. Nat. Hazards Earth Syst. Sci. 2013, 13, 299–309. [Google Scholar] [CrossRef]
  3. Dai, K.; Tie, Y.; Xu, Q.; Feng, Y.; Zhuo, G.; Shi, X. Early Identification of Potential Landslide Geohazards inAlpine-canyon Terrain Based on SAR Interferometry—A Case Study of the Middle Section of Yalong River. J. Radars 2020, 9, 554–568. [Google Scholar]
  4. Dai, K.; Li, Z.; Xu, Q.; Bürgmann, R.; Milledge, D.G.; Tomas, R.; Fan, X.; Zhao, C.; Liu, X.; Peng, J.; et al. Entering the era of earth observation-based landslide warning systems: A novel and exciting framework. IEEE Geosci. Remote Sens. Mag. 2020, 8, 136–153. [Google Scholar] [CrossRef]
  5. Fan, X.; Xu, Q.; Scaringi, G.; Dai, L.; Li, W.; Dong, X.; Zhu, X.; Pei, X.; Dai, K.; Havenith, H.B. Failure mechanism and kinematics of the deadly June 24th 2017 Xinmo landslide, Maoxian, Sichuan, China. Landslides 2017, 14, 2129–2146. [Google Scholar] [CrossRef]
  6. Dai, K.; Li, Z.; Tomás, R.; Liu, G.; Yu, B.; Wang, X.; Cheng, H.; Chen, J.; Stockamp, J. Monitoring activity at the Daguangbao mega-landslide (China) using Sentinel-1 TOPS time series interferometry. Remote Sens. Environ. 2016, 186, 501–513. [Google Scholar] [CrossRef]
  7. Dai, K.; Deng, J.; Xu, Q.; Li, Z.; Shi, X.; Wen, N.; Zhang, L.; Zhuo, G. Interpretation and sensitivity analysis of the LOS displacements from InSAR in landslide measurement. GIScience Remote Sens. 2022, 59, 1226–1242. [Google Scholar] [CrossRef]
  8. Pourkhosravani, M.; Mehrabi, A.; Pirasteh, S.; Derakhshani, R. Monitoring of Maskun landslide and determining its quantitative relationship to different climatic conditions using D-InSAR and PSI techniques. Geomat. Nat. Hazards Risk 2022, 13, 1134–1153. [Google Scholar] [CrossRef]
  9. Ding, K.; He, P.; Wen, Y.; Chen, Y.; Wang, D.; Li, S.; Wang, Q. The 2017 M w 7.3 Ezgeleh, Iran earthquake determined from InSAR measurements and teleseismic waveforms. Geophys. J. Int. 2018, 215, 1728–1738. [Google Scholar] [CrossRef]
  10. Zhao, L.; Liang, R.; Shi, X.; Dai, K.; Cheng, J.; Cao, J. Detecting and Analyzing the Displacement of a Small-Magnitude Earthquake Cluster in Rong County, China by the GACOS Based InSAR Technology. Remote Sens. 2021, 13, 4137. [Google Scholar] [CrossRef]
  11. Béjar-Pizarro, M.; Álvarez Gómez, J.A.; Staller, A.; Luna, M.P.; Pérez-López, R.; Monserrat, O.; Chunga, K.; Lima, A.; Galve, J.P.; Martínez Díaz, J.J.; et al. InSAR-based mapping to support decision-making after an earthquake. Remote Sens. 2018, 10, 899. [Google Scholar] [CrossRef] [Green Version]
  12. Mehrabi, A.; Pirasteh, S.; Rashidi, A.; Pourkhosravani, M.; Derakhshani, R.; Liu, G.; Mao, W.; Xiang, W. Observing new evidence of intensified activity of the Anar fault based on the Persistent Scatterer Interferometry and radon anomaly. Remote Sens. 2021, 13, 2072. [Google Scholar] [CrossRef]
  13. Poreh, D.; Pirasteh, S. InSAR and Landsat ETM+ incorporating with CGPS and SVM to determine subsidence rates and effects on Mexico City. Geoenviron. Disasters 2021, 8, 2–19. [Google Scholar] [CrossRef]
  14. Shi, X.; Chen, C.; Dai, K.; Deng, J.; Wen, N.; Yin, Y.; Dong, X. Monitoring and Predicting the Subsidence of Dalian Jinzhou Bay International Airport, China by Integrating InSAR Observation and Terzaghi Consolidation Theory. Remote Sens. 2022, 14, 2332. [Google Scholar] [CrossRef]
  15. Chen, K.; Liu, G.; Xiang, W.; Chen, R.; Sun, T.; Qian, K.; Pirasteh, S.; Chen, X. Assimilation of SBAS-InSAR Vertical Displacement into Land Surface Model to Improve the Estimation of Terrestrial Water Storage. J. Sel. Top. Appl. Earth Obs. Remote Sens. 2022, 15, 2826–2835. [Google Scholar] [CrossRef]
  16. Di Traglia, F.; De Luca, C.; Manzo, M.; Nolesini, T.; Casagli, N.; Lanari, R.; Casu, F. Joint exploitation of spaceborne and ground-based multitemporal InSAR measurements for volcano monitoring: The Stromboli volcano case study. Remote Sens. Environ. 2021, 260, 112441. [Google Scholar] [CrossRef]
  17. Lazecký, M.; Spaans, K.; González, P.J.; Maghsoudi, Y.; Morishita, Y.; Albino, F.; Elliott, J.; Greenall, N.; Hatton, E.; Hooper, A.; et al. LiCSAR: An automatic InSAR tool for measuring and monitoring tectonic and volcanic activity. Remote Sens. 2020, 12, 2430. [Google Scholar] [CrossRef]
  18. Doke, R.; Kikugawa, G.; Itadera, K. Very local subsidence near the hot spring region in Hakone Volcano, Japan, inferred from InSAR time series analysis of ALOS/PALSAR data. Remote Sens. 2020, 12, 2842. [Google Scholar] [CrossRef]
  19. Yan, S.; Ruan, Z.; Liu, G.; Deng, K.; Lv, M.; Perski, Z. Deriving ice motion patterns in mountainous regions by integrating the intensity-based pixel-tracking and phase-based D-InSAR and MAI approaches: A case study of the Chongce glacier. Remote Sens. 2016, 8, 611. [Google Scholar] [CrossRef]
  20. Zhang, X.; Feng, M.; Zhang, H.; Wang, C.; Tang, Y.; Xu, J.; Yan, D.; Wang, C. Detecting Rock Glacier Displacement in the Central Himalayas Using Multi-Temporal InSAR. Remote Sens. 2021, 13, 4738. [Google Scholar] [CrossRef]
  21. Cigna, F.; Bianchini, S.; Casagli, N. How to assess landslide activity and intensity with Persistent Scatterer Interferometry (PSI): The PSI-based matrix approach. Landslides 2013, 10, 267–283. [Google Scholar] [CrossRef] [Green Version]
  22. Li, Z.; Fielding, E.J.; Cross, P.; Muller, J.P. Interferometric synthetic aperture radar atmospheric correction: GPS topography-dependent turbulence model. J. Geophys. Res. Solid Earth 2006, 111, B02404. [Google Scholar] [CrossRef]
  23. Gong, W.; Meyer, F.J.; Liu, S.; Hanssen, R.F. Temporal filtering of InSAR data using statistical parameters from NWP models. IEEE Trans. Geosci. Remote Sens. 2015, 53, 4033–4044. [Google Scholar] [CrossRef]
  24. Zebker, H.A.; Rosen, P.A.; Hensley, S. Atmospheric effects in interferometric synthetic aperture radar surface deformation and topographic maps. J. Geophys. Res. Solid Earth 1977, 102, 7547–7563. [Google Scholar] [CrossRef]
  25. Williams, S.; Bock, Y.; Fang, P. Integrated satellite interferometry: Tropospheric noise, GPS estimates and implications for interferometric synthetic aperture radar products. J. Geophys. Res. Solid Earth 1998, 103, 27051–27067. [Google Scholar] [CrossRef]
  26. Li, Z.; Ding, X.; Huang, C.; Wadge, G.; Zheng, D.; Physics, S.-T. Modeling of atmospheric effects on InSAR measurements by incorporating terrain elevation information. J. Atmos. Sol.-Terr. Phys. 2006, 68, 1189–1194. [Google Scholar] [CrossRef]
  27. Hanssen, R.F. Radar Interferometry: Data Interpretation and Error Analysis; Science Press: Beijing, China; Springer Science & Business Media: Berlin, Germany, 2001; Volume 2. [Google Scholar]
  28. Neu, D.; Petit, D.; Puyssegur, B.; Michel, M.; Grosprêtre, F.; Blum, A. Nouveautes dans le bilan des tumeurs des parties molles en scanner et IRM. J. Radiol. 2007, 88, 1375. [Google Scholar] [CrossRef]
  29. Delacourt, C.; Briole, P.; Achache, J.J.G.R.L. Tropospheric corrections of SAR interferograms with strong topography. Application to Etna. Geophys. Res. Lett. 1998, 25, 2849–2852. [Google Scholar] [CrossRef]
  30. Wicks, C.W., Jr.; Dzurisin, D.; Ingebritsen, S.; Thatcher, W.; Lu, Z.; Iverson, J. Magmatic activity beneath the quiescent Three Sisters volcanic center, central Oregon Cascade Range, USA. Geophys. Res. Lett. 2002, 29, 26-1–26-4. [Google Scholar] [CrossRef]
  31. Bonforte, A.; Ferretti, A.; Prati, C.; Puglisi, G.; Rocca, F.; Physics, S.-T. Calibration of atmospheric effects on SAR interferograms by GPS and local atmosphere models: First results. J. Atmos. Sol.-Terr. Phys. 2021, 63, 1343–1357. [Google Scholar] [CrossRef]
  32. Bock, Y.; Williams, S. Transactions American Geophysical Union. Integrated satellite interferometry in southern California. Eos Trans. Am. Geophys. Union 1997, 78, 293–300. [Google Scholar] [CrossRef]
  33. Li, Z.; Muller, J.P.; Cross, P. Comparison of precipitable water vapor derived from radiosonde, GPS, and Moderate-Resolution Imaging Spectroradiometer measurements. J. Geophys. Res. Atmos. 2003, 108. [Google Scholar] [CrossRef]
  34. Li, Z.; Muller, J.P.; Cross, P.; Albert, P.; Fischer, J.; Bennartz, R. Assessment of the potential of MERIS near-infrared water vapour products to correct ASAR interferometric measurements. Int. J. Remote Sens. 2007, 27, 349–365. [Google Scholar] [CrossRef]
  35. Shi, X.; Yang, C.; Zhang, L.; Jiang, H.; Liao, M.; Zhang, L.; Liu, X. Mapping and characterizing displacements of active loess slopes along the upstream Yellow River with multi-temporal InSAR datasets. Sci. Total Environ. 2019, 674, 200–210. [Google Scholar] [CrossRef] [PubMed]
  36. Jolivet, R.; Agram, P.S.; Lin, N.Y.; Simons, M.; Doin, M.P.; Peltzer, G.; Li, Z. Improving InSAR geodesy using global atmospheric models. J. Geophys. Res. Solid Earth 2014, 119, 2324–2341. [Google Scholar] [CrossRef]
  37. Parker, A.L.; Biggs, J.; Walters, R.J.; Ebmeier, S.K.; Wright, T.J.; Teanby, N.A.; Lu, Z. Systematic assessment of atmospheric uncertainties for InSAR data at volcanic arcs using large-scale atmospheric models: Application to the Cascade volcanoes, United States. Remote Sens. Environ. 2015, 170, 102–114. [Google Scholar] [CrossRef]
  38. Yu, C.; Li, Z.; Penna, N.T. Interferometric synthetic aperture radar atmospheric correction using a GPS-based iterative tropospheric decomposition model. Remote Sens. Environ. 2018, 204, 109–121. [Google Scholar] [CrossRef]
  39. Beauducel, A.; Debener, S.; Brocke, B.; Kayser, J. On the reliability of augmenting/reducing: Peak amplitudes and principal component analysis of auditory evoked potentials. J. Psychophysiol. 2000, 14, 226. [Google Scholar] [CrossRef]
  40. Li, Z.; Ding, X.; Liu, G. Modeling atmospheric effects on InSAR with meteorological and continuous GPS observations: Algorithms and some test results. J. Atmos. Sol.-Terr. Phys. 2004, 66, 907–917. [Google Scholar] [CrossRef]
  41. Zhu, S.; Yin, P.; Wang, M.; Wang, C.; Wang, W.; Li, J.; Zhao, H. Instability mechanism and disaster mitigation measures of long-distancelandslide at high location in Jinsha River junction zone: Case study of Sela landslide in Jinsha River, Tibet. Chin. J. Geotech. Eng. 2021, 43, 688–697. [Google Scholar]
  42. Xiao, R.; Jiang, M.; Li, Z.; He, X. New insights into the 2020 Sardoba dam failure in Uzbekistan from Earth observation. Int. J. Appl. Earth Obs. Geoinf. 2022, 107, 102705. [Google Scholar] [CrossRef]
  43. Torres, R.; Snoeij, P.; Geudtner, D.; Bibby, D.; Davidson, M.; Attema, E.; Potin, P.; Rommen, B.; Floury, N.; Brown, M. GMES Sentinel-1 mission. Remote Sens. Environ. 2012, 120, 9–24. [Google Scholar] [CrossRef]
  44. Bouchra, A.; Mustapha, H.; Mohammed, R. D-InSAR analysis of Sentinel-1 data for landslide detection in northern Morocco, case study: Chefchaouen. J. Geosci. Environ. Protect. 2020, 8, 84–103. [Google Scholar] [CrossRef]
  45. Goldstein, R.M.; Werner, C. Radar interferogram filtering for geophysical applications. Geophys. Res. Lett. 1998, 25, 4035–4038. [Google Scholar] [CrossRef]
  46. Costantini, M. A novel phase unwrapping method based on network programming. IEEE Trans. Geosci. Remote Sens. 1998, 36, 813–821. [Google Scholar] [CrossRef]
  47. Hooper, A.; Zebker, H.A. Phase unwrapping in three dimensions with application to InSAR time series. JOSA A 2007, 24, 2737–2747. [Google Scholar] [CrossRef] [PubMed]
  48. Tolstikhin, I.O.; Houlsby, N.; Kolesnikov, A.; Beyer, L.; Zhai, X.; Unterthiner, T.; Yung, J.; Steiner, A.; Keysers, D.; Uszkoreit, J.; et al. MLP-Mixer: An all-MLP architecture for vision. Adv. Neural Inf. Process. Syst. 2021, 34, 24261–24272. [Google Scholar]
  49. Chan, K.H.R.; Yu, Y.; You, C.; Qi, H.; Wright, J.; Ma, Y.R. A white-box deep network from the principle of maximizing rate reduction. arXiv 2021, arXiv:2105.10446. [Google Scholar]
  50. Doin, M.P.; Lasserre, C.; Peltzer, G.; Cavalié, O.; Doubre, C. Corrections of stratified tropospheric delays in SAR interferometry: Validation with global atmospheric models. J. Appl. Geophys. 2009, 69, 35–50. [Google Scholar] [CrossRef]
  51. Lin, Y.N.N.; Simons, M.; Hetland, E.A.; Muse, P.; DiCaprio, C. A multiscale approach to estimating topographically correlated propagation delays in radar interferograms. Geochem. Geophys. Geosystems 2010, 11. [Google Scholar] [CrossRef]
  52. Cavalié, O.; Doin, M.P.; Lasserre, C.; Briole, P. Ground motion measurement in the Lake Mead area, Nevada, by differential synthetic aperture radar interferometry time series analysis: Probing the lithosphere rheological structure. J. Geophys. Res. Solid Earth 2007, 112. [Google Scholar] [CrossRef]
  53. Yu, C.; Penna, N.T.; Li, Z. Generation of real-time mode high-resolution water vapor fields from GPS observations. J. Geophys. Res. Atmos. 2017, 122, 2008–2025. [Google Scholar] [CrossRef]
  54. Xiao, R.; Yu, C.; Li, Z.; He, X. Statistical assessment metrics for InSAR atmospheric correction: Applications to generic atmospheric correction online service for InSAR (GACOS) in Eastern China. Int. J. Appl. Earth Obs. Geoinf. 2021, 96, 102289. [Google Scholar] [CrossRef]
Figure 1. (a) Topographic map of the study area; (b) location of the study area; (c,d) images of the study area from Google Earth.
Figure 1. (a) Topographic map of the study area; (b) location of the study area; (c,d) images of the study area from Google Earth.
Remotesensing 14 04171 g001
Figure 2. Flowchart of the proposed method.
Figure 2. Flowchart of the proposed method.
Remotesensing 14 04171 g002
Figure 3. (af) Original interferometric pair; (gl) simulated phase interferometric pair by MLP neural network model; (mr) interferometric pairs after atmospheric correction by MLP neural network model.
Figure 3. (af) Original interferometric pair; (gl) simulated phase interferometric pair by MLP neural network model; (mr) interferometric pairs after atmospheric correction by MLP neural network model.
Remotesensing 14 04171 g003
Figure 4. The correction effects of three topography-dependent atmospheric delay correction methods on InSAR interferograms.
Figure 4. The correction effects of three topography-dependent atmospheric delay correction methods on InSAR interferograms.
Remotesensing 14 04171 g004
Figure 5. Changes in phase standard deviation before and after atmospheric correction by three methods.
Figure 5. Changes in phase standard deviation before and after atmospheric correction by three methods.
Remotesensing 14 04171 g005
Figure 6. Interferometric map 20210604–20210610 phase and elevation linear fit analysis:(a) the original phase with elevation; (b) phase corrected by linear model with elevation; (c) phase corrected by GACOS with elevation; (d) phase corrected by MLP neural network model with elevation.
Figure 6. Interferometric map 20210604–20210610 phase and elevation linear fit analysis:(a) the original phase with elevation; (b) phase corrected by linear model with elevation; (c) phase corrected by GACOS with elevation; (d) phase corrected by MLP neural network model with elevation.
Remotesensing 14 04171 g006
Figure 7. (a) Original displacement results; (b) displacement results after linear model correction; (c) displacement results after GACOS correction; (d) displacement results after correction of MLP neural network model.
Figure 7. (a) Original displacement results; (b) displacement results after linear model correction; (c) displacement results after GACOS correction; (d) displacement results after correction of MLP neural network model.
Remotesensing 14 04171 g007
Figure 8. (a) The MLP neural network atmosphere corrects the displacement results; (b) remote interpretation of displacement.
Figure 8. (a) The MLP neural network atmosphere corrects the displacement results; (b) remote interpretation of displacement.
Remotesensing 14 04171 g008
Figure 9. (a) Original displacement results; (b) displacement results after linear model correction; (c) displacement results after GACOS correction; (d) displacement results after correction of MLP neural network model.
Figure 9. (a) Original displacement results; (b) displacement results after linear model correction; (c) displacement results after GACOS correction; (d) displacement results after correction of MLP neural network model.
Remotesensing 14 04171 g009
Figure 10. (a) The MLP neural network atmosphere corrects the displacement results; (b) remote interpretation of displacement.
Figure 10. (a) The MLP neural network atmosphere corrects the displacement results; (b) remote interpretation of displacement.
Remotesensing 14 04171 g010
Table 1. Sentinel-1 satellite SAR sensor major parameters.
Table 1. Sentinel-1 satellite SAR sensor major parameters.
SAR System
Parameters
ValuesSAR System
Parameters
Values
Date of launchApril 2014The angle of incidence29.1–46.0°
Operation bandCResolution5 m × 20 m
Revisit period12dWidth250 km
Proposed shooting modeIWPolarization modeHH + HV/VV + VH/HH/VV
Table 2. The compared standard deviation of the unwrapping phase before and after correction of MLP neural network.
Table 2. The compared standard deviation of the unwrapping phase before and after correction of MLP neural network.
Interferometric PairPhase Standard Deviation (stdDev)
Original PhaseCorrected PhaseRate of Change
20210421–202105032.58471.080758%
20210421–202105152.57630.689673%
20210503–202105152.48710.637974%
20210517–202106042.50690.785669%
20210523–202106043.16980.619280%
20210604–202106104.68200.750884%
Table 3. The number of, percentage of, and average reduction in interference pairs with reduced phase standard deviation and average increment of interference pairs with increased phase standard deviation of the interferogram after the correction using the three correction methods.
Table 3. The number of, percentage of, and average reduction in interference pairs with reduced phase standard deviation and average increment of interference pairs with increased phase standard deviation of the interferogram after the correction using the three correction methods.
The Interferometric Pairs with a Reduced Standard DeviationLinear ModelGACOSMLP Neural Network
Number 162027
Percentage59%88%100%
Average reduction11.1%17.4%64%
The average increment of interferometric pairs with increasing standard deviation.22.1%10.5%0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, C.; Dai, K.; Tang, X.; Cheng, J.; Pirasteh, S.; Wu, M.; Shi, X.; Zhou, H.; Li, Z. Removing InSAR Topography-Dependent Atmospheric Effect Based on Deep Learning. Remote Sens. 2022, 14, 4171. https://doi.org/10.3390/rs14174171

AMA Style

Chen C, Dai K, Tang X, Cheng J, Pirasteh S, Wu M, Shi X, Zhou H, Li Z. Removing InSAR Topography-Dependent Atmospheric Effect Based on Deep Learning. Remote Sensing. 2022; 14(17):4171. https://doi.org/10.3390/rs14174171

Chicago/Turabian Style

Chen, Chen, Keren Dai, Xiaochuan Tang, Jianhua Cheng, Saied Pirasteh, Mingtang Wu, Xianlin Shi, Hao Zhou, and Zhenhong Li. 2022. "Removing InSAR Topography-Dependent Atmospheric Effect Based on Deep Learning" Remote Sensing 14, no. 17: 4171. https://doi.org/10.3390/rs14174171

APA Style

Chen, C., Dai, K., Tang, X., Cheng, J., Pirasteh, S., Wu, M., Shi, X., Zhou, H., & Li, Z. (2022). Removing InSAR Topography-Dependent Atmospheric Effect Based on Deep Learning. Remote Sensing, 14(17), 4171. https://doi.org/10.3390/rs14174171

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

Article Metrics

Back to TopTop