Next Article in Journal
Stratiform and Convective Rain Classification Using Machine Learning Models and Micro Rain Radar
Next Article in Special Issue
Activity and Kinematics of Two Adjacent Freeze–Thaw-Related Landslides Revealed by Multisource Remote Sensing of Qilian Mountain
Previous Article in Journal
Intercomparison of Planetary Boundary Layer Heights Using Remote Sensing Retrievals and ERA5 Reanalysis over Central Amazonia
Previous Article in Special Issue
Characterizing the Distribution Pattern and a Physically Based Susceptibility Assessment of Shallow Landslides Triggered by the 2019 Heavy Rainfall Event in Longchuan County, Guangdong Province, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An MT-InSAR Data Partition Strategy for Sentinel-1A/B TOPS Data

1
School of Geosciences and Info-Physics, Central South University, Changsha 410083, China
2
Guangzhou Urban Planning and Design Survey Research Institute, Guangzhou 510060, China
3
Guangdong Research Institute of Water Resources and Hydropower, Guangzhou 510635, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(18), 4562; https://doi.org/10.3390/rs14184562
Submission received: 26 July 2022 / Revised: 25 August 2022 / Accepted: 2 September 2022 / Published: 13 September 2022
(This article belongs to the Special Issue Remote Sensing for Natural Hazards Assessment and Control)

Abstract

:
The Sentinel-1A/B satellite launched by European Space Agency (ESA) in 2014 provides a huge amount of free Terrain Observation by Progressive Scans (TOPS) data with global coverage to the public. The TOPS data have a frame width of 250 km and have been widely used in surface deformation monitoring. However, traditional Multi-Temporal Interferometric Synthetic Aperture Radar (MT-InSAR) methods require large computer memory and time when processing full resolution data with large width and long strips. In addition, they hardly correct atmospheric delays and orbital errors accurately over a large area. In order to solve these problems, this study proposes a data partition strategy based on MT-InSAR methods. We first process the partitioned images over a large area by traditional MT-InSAR method, then stitch the deformation results into a complete deformation result by correcting the offsets of adjacent partitioned images. This strategy is validated in a flat urban area (Changzhou City in Jiangsu province, China), and a mountainous region (Qijiang in Chongqing City, China). Compared with traditional MT-InSAR methods, the precision of the results obtained by the new strategy is improved by about 5% for Changzhou city and about 15% for Qijiang because of its advantage in atmospheric delay correction. Furthermore, the proposed strategy needs much less memory and time than traditional methods. The total time needed by the traditional method is about 20 h, and by the proposed method, is about 8.7 h, when the number of parallel processing is 5 in the Changzhou city case. The time will be further reduced when the number of parallel processes increases.

Graphical Abstract

1. Introduction

Due to large coverage and high-precision, Interferometric Synthetic Aperture Radar (InSAR) has been widely used for mapping surface deformation, such as urban surface deformation [1,2], seismic deformation [3,4], landslide displacement [5,6,7,8,9], and mining subsidence [10]. With the fast development of SAR satellite technology [11], the observation range and frequency are both improved [12,13,14], providing cycle monitoring for a large-scale or national wide area. However, the traditional processing strategies for Multi-Temporal Interferometric Synthetic Aperture Radar (MT-InSAR) cannot efficiently process the huge number of images with large spatial and temporal coverage. Furthermore, the possible atmospheric phase screen and orbital errors exist in the SAR images with wide spatial coverage are difficult to be corrected. Therefore, optimizing the InSAR processing strategy and parameters is crucial for the application of wide-area InSAR data.
Using supercomputers or distributed computing systems, such as CASearth Cloud Infrastructure Platform [15], and ESA’s G-POD environment [16], is a way to improve the data processing efficiency, but it is too expensive to be popularized. Another way is to segment a large image into small blocks, which can significantly reduce the computation burden and complexity in one block and improve the efficiency of data processing. Currently, data partition strategies are applied in some steps of SAR data processing, such as phase unwrapping [17,18,19,20], orbital error correction [21], atmospheric correction [22,23], and PS point decomposition [24,25], but not the whole data process. GAMMA software provides a well-known patch-based point target analysis method, Interferometric Point Target Analysis (IPTA) [26]. However, the method of using local reference points between neighboring patches to extend the results to adjacent blocks is highly affected by unstable connections, resulting in errors propagating in the result easily [27]. StaMPS method and software [28] also provide a block strategy to select permanent scatterers, but it is time-consuming for large-scale areas [29]. Data block processing often introduces systematic errors, such as reference basis errors [29,30,31]. To remove these errors, external data (GNSS and leveling) and modeling [32,33] are needed. Furthermore, the correction efficacy and precision strongly depend on the precision and spatial distribution of external data.
To address the above problems, we propose a strategy to divide the original data into small blocks by GAMMA software and process these blocks independently by the traditional MT-InSAR method. Then, we use the least square method to estimate the basis between each block and mosaic the corrected block results to obtain the overall results. To validate our strategy, we selected the Sentinel-1 Terrain Observation by Progressive Scans (TOPS) data of a city in the plain area (Changzhou City, Jiangsu Province) and a city in the mountainous area (Qijiang, Chongqing City) in China for the experiment. The results obtained by the traditional and the proposed methods are compared in terms of precision, memory consumption, and time consumption. We also discuss the optimal overlap ratio of blocks and the application of the proposed method.
This study is organized as follows: Section 2 describes the proposed method in detail. The study area and the datasets are introduced in Section 3. In Section 4, we compare the precision and time consumption of the proposed method and traditional method. The block approach and the applicability scenarios of our method is also discussed in Section 5. Finally, some conclusions are drawn in Section 6.

2. The Block MT-InSAR Data Processing Strategy

In order to solve the great calculation burden caused by Sentinel-1A/B TOPS data of large spatial and temporal coverage, this paper proposes a data partition strategy based on the MT-InSAR data method, referred to as the block MT-InSAR algorithm. First, the TOPS data are co-registered in the study area to obtain registered single look complex (RSLC), and then the RSLC data are partitioned and processed separately by the traditional MT-InSAR algorithm. Then, the results are corrected by the adjustment model based on the spatial consistency of homonymy points (the same ground deformation points located in different blocks within the overlap areas.). Finally, the results are spliced to obtain the continuous overall deformation results. The general flow of the method is shown in Figure 1 and Figure A1.

2.1. Data Partition and Block Processing

Even partition [29], quadtree partition [22], and clustering algorithm partition [19,34] have been used in some parts of data processing, such as atmospheric delay removal, PS network construction, and phase unwrapping. In order to facilitate the splicing process, this paper uses even partition to divide the original data into small blocks, in which the block size and the overlap ratio should be considered.
Block size affects the precision and the processing efficiency of the phase unwrapping, atmospheric delay, and orbit error partition. If the block interferograms are highly coherent and easy to unwrap, the block size has little effect on phase unwrapping precision, but a small block size would lead to high unwrapping efficiency [17]. Additionally, the atmospheric delay in the MT-InSAR consists of a short-scale (few kilometers) and a long-scale (tens of kilometers) component [35], so the block size smaller than these scales is conducive to removing atmospheric delay. However, the too-small block size may remove the long-wavelength deformation signal. The block width and height should be larger than 1/3 ALOS-2 data in range and azimuth for ALOS-2 (70 km) datasets [21]. Therefore, we set the initial block size as ~30 km in length and width for S1A/B TOPS data.
The overlap ratio between blocks also affects the reliability of results and data processing efficiency. The larger the overlap ratio the higher the reliability. For example, if a block area is overlapped with the surrounding blocks in four directions by 10%, 36% of the small block is overlapped with the surrounding blocks, and the overlap area will become 96% when the overlap ratio is 40% in width and height directions. (Figure 2a) However, increasing the overlap ratio will lead to a lot of repeated calculations and reduce the data processing efficiency. In order to improve the result reliability (>50% overlap area), we take off the balance between the reliability of results and choose the overlap ratio of 20~40% for further experiments.

2.2. Results Correction Based on Least Square Estimation

After data partition, we process each block of data to obtain the deformation results using the improved IPTA-InSAR method [36,37]. The obtained deformation results of all the blocks are preprocessed through three steps. (1) Co-registration. Due to the location errors caused by orbital errors and low resolution of DEM; the location of deformation results may have a systematic deviation of about 1–2 pixels after geocoding. Such deviation can be solved by an overall offset correction using some feature points on the ground. (2) Automatic extraction of overlapping regions. The deformation rates of the homonymy points in the overlapping region determine the correction precision, so it is necessary to identify the overlapping regions between the deformation results first. We use the topological relationship between image overlays (quadrangles) to find the coordinates of corner points in the overlap region (Figure 2b). (3) High-quality homonymy point selection. We select homonymy points with high coherence and small uncertainty. After these operations, the deformation results can be corrected by adjustment.
Errors can be removed during data processing. However, data partition makes each block have a local reference point, and the benchmarks of these reference points may be different, resulting in discrepancies between the results of adjacent blocks and affecting the precision of the overall results. The difference in the deformation rates of the homonymy points in the overlapping area is described as
v e l i , k v e l j , k = δ i , k δ j , k
where v e l i , k and v e l j , k denote the deformation rate at deformation point k in image i and j, respectively, and δ i , k , and δ j , k denote the error of the corresponding points. Since the error contains mainly the difference in benchmarks, this value can be assumed as a constant.
The matrix form of Equation (1) is:
V = B X ^ L
where V = [ δ v 1 δ v 2   δ v i   δ v M ] T is the residual of the calculated values and the observations, X ^ = [ x ^ 1 x ^ 2 x ^ M ] T is the difference between the reference points of adjacent images estimated by the least square method. B is the coefficient matrix. L = [ v e l 1 , k v e l 2 , k v e l i , k v e l j , k ] T .
To solve Equation (2), we have to determine the weights of the blocks according to the quality of the data involved in the adjustment.
D ( L ) = σ 0 2 P 1
σ 0 2 denotes the variance of unit weight and P is the weight matrix. Assume that the uncertainty of point i is given by δ. Then, the weight of the point is
p i = C δ i
In this study, the data are partitioned into small blocks, which are processed independently. The Helmert variance component estimation for multiple data classes is applied to optimize the solution weights of each data set.
S θ ^ = W θ
where θ ^ = [ σ ^ 01 2 σ ^ 02 2 σ ^ 0 M 2 ] T is the estimated variance of unit weight. W θ is the square sum of the corrected values, W θ = [ V 1 T P V 1 V 2 T P V 2 V M T P V M ] T . S is the coefficient matrix. After obtaining θ ^ , X ^ is solved using the least square method. Repeat the above process until θ ^ satisfies the given threshold T = 3 δ 0 , and the corresponding solution is the optimal X ^ for each SAR image block. The corrected deformation rate is obtained by Equation (6).
v e l ^ i , k = v e l i , k x ^ i
The posteriori variance of unit weight and the covariance array are used to evaluate the adjustment observation. They can be obtained by
σ ^ 0 2 = V T P V r
Q x ^ i = B i T P B i
In Equation (7), r is the number of redundant observations, and it can be referred to as the number of degrees of freedom, r = N M , with N denoting the number of observations, and M denoting the row number of X ^ . According to the error propagation, the covariance of the estimate of the homonymy points in the overlap region can be obtained by Equation (9).
Q L ^ L ^ = ( B Q x ^ 1 B T P ) Q ( B Q x ^ 1 B T P ) T = B Q x ^ 1 B T
To verify the precision of the adjustment results, the deformation difference of the homonymy points before and after correction are compared. The block processing results are verified by comparing with that of the traditional processing method (the result without partitioning processing).

2.3. Result Mosaicking

The final step is to mosaic the corrected deformation results of all blocks. After geocoding, the block results are horizontally mosaicked. After correction, the deformation of the homonymy points in the overlapping area may still have differences, due to the different errors distribution. We adopt the weighted average method to merge the deformation of the homonymy points.
After correcting the deformation rate, we correct the deformation sequence. Assuming that the deformation is linear, T ( T 1 , T 2 T a T t ) is a deformation sequence, and the corrected deformation sequence at time T a can be obtained by
S ^ a = T 1 T a ( v a + x ^ a ) d T = S a + x ^ a ( T a T 1 )
where S ^ a is the accumulated deformation after correction, v a is the deformation rate. x ^ a is the correction of deformation rate, but cannot be calculated, S a is the accumulated deformation before correction. If the deformation is linear, x ^ a is equal to the correction of average deformation rate in the deformation sequences T . Variant x ^ can be calculated by Equation (2), so the equation can be instead of Equation (11). Figure 3 is the diagram of deformation sequence correction.
S ^ a = S a + x ^ ( T a T 1 )

3. Experiment and Data Processing

3.1. Study Area and Datasets

Two study areas are selected to validate the proposed method. One is in Changzhou City (31°09′–32°04′N, 119°08′–120°12′E), a coastal city in eastern China. This area is a plain with an elevation of about 10 m [38]. It has a highly developed economy and urban industry. The continuous expansion of urban and engineering projects has changed the geological environment and led to frequent geological hazards. So, surface subsidence monitoring in the area is necessary.
The other study area is Qijiang (28°27′–29°11′N, 106°23′–107°03′E), in western China. It is in the transition zone from the southeastern edge of the Sichuan basin to the Yunnan-Guizhou plateau. The topography is undulating. The mountainous area accounts for 67.6% of the total area and the hills account for 32.4%. The average elevation of this area is 254.8 m [39].
These two regions are used to test the applicability of the proposed method under different error conditions.
We collected 110 Sentinel-1A/B TOPS images covering Changzhou City from path T69 and frame 99, between 5 January 2018 and 31 December 2020, and acquired 115 Sentinel-1 images over Qijiang from path T55 and frame 92, between 9 January 2018 and 31 December 2021. Specific image parameters and image acquisition time are shown in Table 1 and Table A1 in Appendix A.

3.2. Data Processing

Using the method described in Section 2, we partitioned the acquired single look complex (SLC) images after co-registration and obtained 30 small blocks with overlapping regions (Figure 4). The block size in Changzhou City is about 7000 × 1400 (pixels), and the overlap rate is about 30%; the block size in Qijiang is about 6400 × 1600 (pixels), and the overlap rate is about 25%.
The spatial baselines of Sentinel-1 images are short, so we connected each image with two (temporally) adjacent images to form a network, only considering the temporal baselines. A multi-look operation (range: azimuth = 5:1) was applied to reduce the noise. After the multi-look operation, the image size of Changzhou city was reduced to 1600 × 1400 (pixels) and that of Qijiang was reduced to 1500 × 1600 (pixels). The data were processed by minimum cost flow (MCF) for phase unwrapping, and Goldstein filtering for noise mitigation. The PS points were selected considering the phase coherence threshold and the amplitude dispersion threshold of the amplitude map. Orbital error phases were removed by polynomial fitting. Most atmospheric phases were removed by differencing between neighboring PS points, and the remaining was removed by spatial-temporal filtering. The topographic residual phases were then removed using linear regression. Finally, the deformation sequence was solved from the remaining phases using Singular Value Decomposition (SVD). The obtained time series of deformation was corrected by the method introduced in Section 2.2. The average deformation of the high-quality homonymy points in the overlapping areas was used for correction. Finally, the result of traditional processing and block processing were obtained. Figure 5 shows the deformation results of Changzhou and Qijiang.

4. Result Analysis

When dealing with the deformation time series of a large area, most conventional algorithms use one reference point for phase unwrapping and solve for PS point deformation rates. If the distance between the PS point and the reference point is large, the precision of the results is low. Reducing the size of image coverage by partition can improve the precision of PS points. However, partition leads to different reference points for different blocks, so the deformation results should be corrected to follow one benchmark. In this section, the partition and the traditional methods are compared in terms of precision and time consumption.

4.1. Precision of the Deformation Rate

We compare the deformation results in Changzhou found by the partition method and the traditional method in Figure 5. The two results show a similar distribution of deformation, but a slight difference in details. We manually selected three regions A, B, and C for analysis. The size of these regions is 1000 × 1000 (pixels) These three regions are stable and outside the deformation region, so we assumed the deformation as 0. The statistical analysis shows that the standard deviation (STD) of the partition results and the traditional results in region A is 3.4 mm/yr and 3.7 mm/yr, respectively, in region B is 4.1 mm/yr and 4.2 mm/yr, respectively, and in region C is 4.6 mm/yr and 4.6 mm/yr, respectively. On the whole, the precision of the results obtained by the partition strategy is slightly higher than that of the traditional processing results.
Figure 5c,d show the deformation results of Qijiang found by the two methods, which generally agreed with each other but some local areas have some slight differences, especially in the circled areas D, E, F (the selection criteria is the same as A, B, C). The traditional results contain a large number of uplift signals, which are not deformation signals but residual errors. These errors are significantly less in the partition results. The STD of the partition results in regions D, E, and F are 3.2 mm/yr, 3.1 mm/yr, and 3.6 mm/yr, respectively, and the correspondence of the traditional results are 3.8 mm/yr, 3.9 mm/yr, and 4.0 mm/yr, respectively. Therefore, the partition method outperforms the traditional method in error removal.
The comparison results in Table 2 show that the partition method has higher precision than the traditional method. In the Changzhou experiment, the former obtained a precision of about 5% higher than the traditional method, and in the Qijiang experiment, the precision improvement is about 15%.

4.2. Precision of the Deformation Sequence

After correcting the deformation rates, we corrected the corresponding deformation sequences. We selected the deformation time series of 3 points in each of Changzhou city (P1, P2, P3) and Qijiang (P4, P5, P6) to test the result precision. These points are in the overlapping regions, and they have different deformation magnitudes. P1 and P4 have large deformation rate, P2 and P5 have medium deformation rate, and P3 and P6 have small deformation rate. The results are shown in Figure 5. In Changzhou City, the difference between the deformation sequences obtained by our method and that obtained by the traditional method is not significant. At P2, the deformation sequences obtained by the two methods almost coincide (Figure 5f), and the deformation rate difference at the three selected points is less than 1 mm/yr. The RMSE between the two results is 5.3 mm for P1, 3.2 mm for P2, and 4.0 mm for P3. We also selected the time series deformation of the three points in Qijiang. The overall deformation trend and deformation magnitude obtained by the two methods are basically the same. The RMSEs between the two results at P4, P5, and P6 are 6.2 mm, 5.0 mm, and 3.0 mm, respectively.
In Figure 5e,h, the annual average deformation rates of these two points are more than 30 mm/yr, the difference between the deformation rates of these two points was about 1.5 mm/yr, and the deformation monitoring precision of InSAR was also basically in this range. A simple proportional function model is used for correcting the time series. If the deformation is nonlinear, the correction of this model might not be appropriate.

4.3. Time and Memory Consumption

The partition and parallel processing strategy can reduce the memory consumption of every single process. Additionally, increasing the number of parallel processing can reduce the time consumed by the whole data processing. The partition strategy can be roughly divided into three stages: image partition, MT-InSAR processing, and correction. The total time consumed by partition depends mainly on the processing times of image blocks, that is, the total number of blocks divided by the number of blocks processed in a single parallel session. The traditional processing spends all its time on MT-InSAR processing.
The program running time and memory consumption of the two strategies are listed in Table 3. The traditional method costs about 20 h. The total time required for block processing is 46.6 h when the number of parallel processing is 1. However, when the number of parallel processing is greater than 2, the block processing needs less time than the traditional processing. Additionally, it only needs 8.7 h when the number of parallel processing is 5. In addition, the computer memory occupied by block processing is much lower than that of traditional processing. In this experiment, the single memory occupied by partitioning is only 1/20 of that occupied by the traditional processing. When the image coverage is large or the computer memory is small, the traditional processing may cause memory overflows and the data cannot be processed successfully, but this problem will not happen to our data partition processing.

5. Discussion

5.1. Space Consistency Correction

The precision of the deformation in each block can be calculated by Formula (9). In the results of Changzhou, the unit weight mean error after adjustment is 0.13 mm/yr, the precision of the block with the highest adjustment precision is 0.27 mm/yr, and the precision of the block with the lowest adjustment precision is 0.59 mm/yr. In the results of Qijiang, the unit weight mean error after adjustment is 0.28 mm/yr, the precision of the block with the highest adjustment precision is 0.40 mm/yr, and the precision of the block with the lowest adjustment precision is 0.81 mm/yr. The adjustment precision is plotted in Figure 6 and Figure 7, which show that the precision of the center blocks is higher than that of the edge blocks.
In theory, after adjustment, the deformation rates of the homonymy points in the overlapping region should be the same. As Figure 6c shows, before the adjustment, the difference between homonymy points in Changzhou is more than 5 mm/yr, with an STD of 5.1 mm/yr, and after the adjustment, the difference almost converges to 0, with an STD of 0.5 mm/yr. In Qijiang, the difference between homonymy points is more than 7 mm before adjustment, with an STD of 3.9 mm/yr, and it is reduced to 0.6 mm/yr after the adjustment. The precision of the block processing results in the two study areas was greatly improved by adjustment, indicating that adjustment can improve the consistency of the block deformation results.
We selected 4 deformation areas (A, B, C, D) in the overlapping area in Changzhou (Figure A2), and compared their deformation results before and after adjustment in Figure 8. The results of the two image blocks in region A have little difference, so the improvement of the result is not significant after the adjustment (Figure 8b,f). However, the results of regions B, C and D are improved obviously after adjustment. The mean values of the differences of deformation in these three regions change from 4.1 mm/yr, −5.6 mm/yr, and 5.6 mm/yr to 0.3 mm/yr, 0.1 mm/yr, and 0.2 mm/yr after correction, and the spatial consistency of the results improves more significantly.

5.2. Effects of Overlap Rate on Result Precision and Time Consumption

In the proposed partition strategy, partitioning the image is the first and most critical step. Different partition strategies provide different precision results, and different data processing efficiency. Because we adopt the even partition strategy, the size and overlap ratio of the blocks have a great impact on the results. To analyze the effects of the overlap ratio on the partition results, we set one block size and obtained the temporal deformation results of Changzhou using the overlap ratios of 10%, 20%, 30%, and 40%, separately. The results are shown in Figure A2. We evaluated the result precision (Table 4).
The precision of the deformation rates obtained by different overlap ratios is similar because the mean values of the high-quality homonymy points in the overlap region are used in adjustment, which is slightly influenced by the overlap region. As long as the block size is the same, the precision of the results obtained by different overlap ratios are similar. However, the time consumed by different overlap ratios is different. The larger the overlap ratio, the larger the number of blocks and the longer the processing. When the number of parallel processing is 5, the total time consumed by the four overlap ratios are 6.6 h, 7.6 h, 9.5 h, and 10.7 h, indicating that the consuming time increases with the increase of overlap ratio. The increase in the overlapping area brings larger double-counted areas and raises the reliability of the results. Considering the precision, reliability, and time consumption, we choose 20% as the best overlap ratio. The experiments show that overlap ratio 20% has similar result precision and time consumption with that of overlap ratio 10%, but it leads to more than 64% overlap area, which contributes to a significantly higher reliability.

5.3. Implications of Data Partition Strategy for MT-InSAR

The administrative boundary of a city is usually an irregular polygon, but the image coverage is a regular quadrangle. Thus, the image coverage contains many data unrelated to the study area. The traditional method will also process these data. If such data accounts for a large proportion of the image, the data processing will waste a lot of time. The proposed method only processes the block data inside the study area, which can improve the data processing efficiency. In Figure 4, the gray blocks do not need processing. Furthermore, the proposed method can refine the data processing for only the blocks with deformation, which further improves the efficiency of data processing.
The difference in the image coverage will definitely lead to the difference in the results. In this paper, we divide data into blocks, and process, correct, and splice the results of all blocks. The atmospheric delay in small range data is easier to remove than that in the data with large range. Studies have shown that the atmospheric phase in InSAR data measurements has a close correlation with spatial scale [35]. The atmospheric phase difference between two PS points with a distance less than 1 km is less than 0.1 rad2 [40], so the smaller the area, the better the atmospheric error removal according to the error propagation law. However, for large deformation areas, long-wavelength deformation may be removed as orbital errors, due to the polynomial fitting [21]. Thus, the proposed method is not fully applicable to the study area with long-wavelength deformation, such as interseismic deformation.
The correction of the partition results is based on the assumption that the deformation rates of the homonymy points are the same. However, the deformation acquired by InSAR is the line-of-sight (LOS) deformation, and the deformation direction at each point is related to the incidence angle. When the deformation of the homonymy point is obtained from the same orbit and has the same incidence angles, the deformation rates should be the same. If the partitioned data are acquired under different imaging geometries, there will be inconsistency in the incidence angles, resulting in different LOS deformation. Therefore, the incidence angle variation of the results should be considered when the partition data are acquired from different orbits.
Finally, the method does not use control points for the adjustment. Although the benchmarks between image blocks are unified, there may be a deviation between the unified benchmark and the real deformation result datum. We only make a simple correction to the result but using external data as control points may improve the correction.

6. Conclusions

In this paper, we propose to partition the data into blocks before obtaining the deformation, to save memory and time for large-scale data processing. To validate this method, we used the Sentinel-1 TOPS data covering Changzhou, a plain area, and Qijiang, a mountainous area in China. The time series deformation results were obtained in these two regions using the traditional processing method (the improved IPTA) and the partition processing method. The latter outperforms the former in precision, time consumption, and memory occupation. Taking Changzhou City as an example, the memory occupation of the traditional processing method is about 27.2 G, and the total time consumed by processing is about 20 h. During partition processing, the memory occupation of each block is only 1.3 G, and the consumed time is 8.7 h when the parallel number is 5. We also compared the precision of the results obtained by the two methods. The results obtained by the partition processing in Changzhou is as about 4.0 mm/yr, while the precision of the traditional processing is about 4.2 mm/yr. The correspondence in Qijiang is about 3.3 mm/yr and 3.9 mm/yr, respectively. The precision of the results obtained by the proposed method is higher than that obtained by traditional processing.
In general, the proposed method can significantly reduce the memory occupation and time consumption of data processing under the condition of sufficient parallelism, and the precision of the results is higher than that obtained by traditional processing. This method is suitable for monitoring the short-wavelength deformation in a large area, such as large-scale urban deformation monitoring and large-scale landslide deformation detection. However, further research is needed for the result splicing and its application to long-wavelength deformation.

Author Contributions

Conceptualization, Y.W. (Yuexin Wang) and G.F.; data curation, S.L.; formal analysis, Z.F. and Y.W. (Yuedong Wang); funding acquisition, Y.W. (Yuexin Wang) and G.F.; methodology, Y.W. (Yuexin Wang) and G.F.; resources, G.F.; Software, Y.W. (Yuexin Wang); supervision, S.L., Y.Z. and H.L.; validation, X.W.; writing—original draft, Y.W. (Yuexin Wang); writing—review and editing, Y.W. (Yuexin Wang), G.F., Z.F., Y.W. (Yuedong Wang), X.W., Y.Z. and H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 42174039), and the Fundamental Research Funds for the Central Universities of Central South University (No. 506021741).

Acknowledgments

The authors would like to thank the European Space Agency (ESA) for providing free Sentinel-1 data. The data were additionally retrieved from the Alaska Satellite Facility Distributed Active Archive Center. We also thank the contributors for the Generic Mapping Tools (Wessel et al., 2013) open-source software.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Parameters of the images used in this study.
Table A1. Parameters of the images used in this study.
Study AreaParametersAcquisition Date (YYYY/MM/DD)
ChangzhouDirectionAscending2018/01/102018/01/222018/02/032018/02/15
PathT692018/02/272018/03/112018/03/232018/04/04
Heading−12.79°2018/04/162018/04/282018/05/102018/05/22
Incidence36.65°2018/06/032018/06/152018/06/272018/07/09
Pixel Spacing
(Rg × Az)
2.33 × 13.982018/07/212018/08/022018/08/142018/09/07
2018/09/192018/10/012018/10/132018/10/25
Number of
images
1102018/11/062018/11/182018/12/122018/12/24
2019/01/052019/01/172019/02/102019/02/16
2019/02/222019/03/062019/03/182019/03/30
2019/04/052019/04/112019/04/232019/04/29
2019/05/052019/05/112019/05/172019/05/23
2019/05/292019/06/042019/06/102019/06/16
2019/06/222019/06/282019/07/042019/07/10
2019/07/162019/07/222019/07/282019/08/03
2019/08/092019/08/152019/08/212019/08/27
2019/09/022019/09/082019/09/202019/09/26
2019/10/022019/10/082019/10/142019/10/20
2019/10/262019/11/012019/11/072019/11/19
2019/11/252019/12/012019/12/072019/12/13
2019/12/192019/12/252019/12/312020/01/12
2020/01/242020/02/052020/02/172020/02/29
2020/03/122020/03/242020/04/052020/04/17
2020/04/292020/05/112020/05/232020/06/04
2020/06/162020/06/282020/07/102020/07/22
2020/07/282020/08/032020/08/152020/08/27
2020/09/082020/09/202020/10/022020/10/14
2020/10/262020/11/072020/11/192020/12/01
2020/12/132020/12/25
QijiangDirectionAscending2018/01/092018/01/212018/02/022018/02/14
PathT552018/02/262018/03/102018/03/222018/04/03
Heading−12.65°2018/04/152018/04/272018/05/092018/05/21
Incidence43.64°2018/06/022018/06/142018/06/262018/07/08
Pixel Spacing
(Rg × Az)
2.33 × 13.962018/07/202018/08/012018/08/252018/09/06
2018/09/182018/09/302018/10/122018/10/24
Number of
images
1142018/11/052018/11/292018/12/112018/12/23
2019/01/042019/01/162019/01/282019/02/09
2019/02/212019/03/052019/03/172019/03/29
2019/04/102019/04/222019/05/042019/05/16
2019/05/282019/06/092019/07/032019/07/15
2019/07/272019/08/082019/08/202019/09/01
2019/09/132019/09/252019/10/072019/10/19
2019/10/312019/11/122019/11/242019/12/06
2019/12/182019/12/302020/01/112020/01/23
2020/02/042020/02/162020/02/282020/03/11
2020/03/232020/04/042020/04/162020/04/28
2020/05/222020/06/032020/06/152020/06/27
2020/07/092020/07/212020/08/022020/08/14
2020/09/072020/09/192020/10/012020/10/13
2020/10/252020/11/062020/11/182020/11/30
2020/12/122020/12/242021/01/052021/01/17
2021/01/292021/02/102021/02/222021/03/06
2021/03/182021/03/302021/04/112021/04/23
2021/05/292021/06/102021/06/222021/07/16
2021/07/282021/08/092021/08/212021/09/02
2021/09/142021/09/262021/10/082021/10/20
2021/11/012021/11/132021/11/252021/12/07
2021/12/192021/12/31
Figure A1. The detailed flowchart of the proposed method.
Figure A1. The detailed flowchart of the proposed method.
Remotesensing 14 04562 g0a1
Figure A2. The deformation results in Changzhou obtained using the overlap ratio of (a) 10%, (b) 20%, (c) 30%, and (d) 40%. A1-C4 are the same areas as A–C described in Section 4.1.
Figure A2. The deformation results in Changzhou obtained using the overlap ratio of (a) 10%, (b) 20%, (c) 30%, and (d) 40%. A1-C4 are the same areas as A–C described in Section 4.1.
Remotesensing 14 04562 g0a2

References

  1. Ng, A.H.-M.; Ge, L.; Li, X.; Zhang, K. Monitoring ground deformation in Beijing, China with persistent scatterer SAR interferometry. J. Geod. 2011, 86, 375–392. [Google Scholar] [CrossRef]
  2. Wang, H.; Wright, T.J.; Yu, Y.; Lin, H.; Jiang, L.; Li, C.; Qiu, G. InSAR reveals coastal subsidence in the Pearl River Delta, China. Geophys. J. Int. 2012, 191, 1119–1128. [Google Scholar] [CrossRef]
  3. Feng, G.; Li, Z.; Shan, X.; Xu, B.; Du, Y. Source parameters of the 2014 Mw 6.1 South Napa earthquake estimated from the Sentinel 1A, COSMO-SkyMed and GPS data. Tectonophysics 2015, 655, 139–146. [Google Scholar] [CrossRef]
  4. Chaussard, E.; Johnson, C.W.; Fattahi, H.; Bürgmann, R. Potential and limits of InSAR to characterize interseismic deformation independently of GPS data: Application to the southern San Andreas Fault system. Geochem. Geophys. Geosystems 2016, 17, 1214–1229. [Google Scholar] [CrossRef]
  5. Dong, J.; Zhang, L.; Tang, M.; Liao, M.; Xu, Q.; Gong, J.; Ao, M. Mapping landslide surface displacements with time series SAR interferometry by combining persistent and distributed scatterers: A case study of Jiaju landslide in Danba, China. Remote Sens. Environ. 2018, 205, 180–198. [Google Scholar] [CrossRef]
  6. Xiong, Z.; Feng, G.; Feng, Z.; Miao, L.; Wang, Y.; Yang, D.; Luo, S. Pre- and post-failure spatial-temporal deformation pattern of the Baige landslide retrieved from multiple radar and optical satellite images. Eng. Geol. 2020, 279, 105580. [Google Scholar] [CrossRef]
  7. Novellino, A.; Cesarano, M.; Cappelletti, P.; di Martire, D.; di Napoli, M.; Ramondini, M.; Sowter, A.; Calcaterra, D. Slow-moving landslide risk assessment combining Machine Learning and InSAR techniques. CATENA 2021, 203, 105317. [Google Scholar] [CrossRef]
  8. Meng, Q.; Confuorto, P.; Peng, Y.; Raspini, F.; Bianchini, S.; Han, S.; Liu, H.; Casagli, N. Regional Recognition and Classification of Active Loess Landslides Using Two-Dimensional Deformation Derived from Sentinel-1 Interferometric Radar Data. Remote Sens. 2020, 12, 1541. [Google Scholar] [CrossRef]
  9. Miele, P.; Di Napoli, M.; Novellino, A.; Calcaterra, D.; Mallorqui, J.J.; Di Martire, D. SAR data and field surveys combination to update rainfall-induced shallow landslide inventory. Remote Sens. Appl. Soc. Environ. 2022, 26, 100755. [Google Scholar] [CrossRef]
  10. Yang, Z.; Li, Z.; Zhu, J.; Preusse, A.; Hu, J.; Feng, G.; Papst, M. High-Resolution Three-Dimensional Displacement Retrieval of Mining Areas From a Single SAR Amplitude Pair Using the SPIKE Algorithm. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 3782–3793. [Google Scholar] [CrossRef]
  11. Ghasemloo, N.; Matkan, A.A.; Alimohammadi, A.; Aghighi, H.; Mirbagheri, B. Estimating the Agricultural Farm Soil Moisture Using Spectral Indices of Landsat 8, and Sentinel-1, and Artificial Neural Networks. J. Geov. Spat. Anal. 2022, 6, 19. [Google Scholar] [CrossRef]
  12. Kellogg, K.; Hoffman, P.; Standley, S.; Shaffer, S.; Rosen, P.; Edelstein, W.; Dunn, C.; Baker, C.; Barela, P.; Shen, Y.; et al. NASA-ISRO Synthetic Aperture Radar (NISAR) Mission. In Proceedings of the 2020 IEEE Aerospace Conference, Big Sky, MT, USA, 7–14 March 2020; pp. 1–21. [Google Scholar] [CrossRef]
  13. Fan, W.; Pan, G.; Wang, L. Development and Application of a Networked Automatic Deformation Monitoring System. J. Geovisualization Spat. Anal. 2020, 4, 11. [Google Scholar] [CrossRef]
  14. Torres, R.; Snoeij, P.; Geudtner, D.; Bibby, D.; Davidson, M.; Attema, E.; Potin, P.; Rommen, B.; Floury, N.; Brown, M.; et al. GMES Sentinel-1 mission. Remote Sens. Environ. 2012, 120, 9–24. [Google Scholar] [CrossRef]
  15. Wang, C.; Tang, Y.; Zhang, H.; You, H.; Zhang, W.; Duan, W.; Wang, J.; Dong, L.; Zhang, B. First mapping of China surface movement using supercomputing interferometric SAR technique. Sci. Bull. 2021, 66, 1608–1610. [Google Scholar] [CrossRef]
  16. De Luca, C.; Cuccu, R.; Elefante, S.; Zinno, I.; Manunta, M.; Casola, V.; Rivolta, G.; Lanari, R.; Casu, F. An On-Demand Web Tool for the Unsupervised Retrieval of Earth’s Surface Deformation from SAR Data: The P-SBAS Service within the ESA G-POD Environment. Remote Sens. 2015, 7, 15630–15650. [Google Scholar] [CrossRef]
  17. Chen, C.W.; Zebker, H.A. Phase unwrapping for large SAR interferograms: Statistical segmentation and generalized network models. IEEE Trans. Geosci. Remote Sens. 2002, 40, 11. [Google Scholar] [CrossRef]
  18. Zhang, K.; Ge, L.; Hu, Z.; Ng, A.H.-M.; Li, X.; Rizos, C. Phase Unwrapping for Very Large Interferometric Data Sets. IEEE Trans. Geosci. Remote Sens. 2011, 49, 4048–4061. [Google Scholar] [CrossRef]
  19. Yu, H.; Xing, M.; Bao, Z. A Fast Phase Unwrapping Method for Large-Scale Interferograms. IEEE Trans. Geosci. Remote Sens. 2013, 51, 4240–4248. [Google Scholar] [CrossRef]
  20. Yuan, Z.; Chen, T.; Xing, X.; Peng, W.; Chen, L. BM3D Denoising for a Cluster-Analysis-Based Multibaseline InSAR Phase-Unwrapping Method. Remote Sens. 2022, 14, 1836. [Google Scholar] [CrossRef]
  21. Du, Y.; Fu, H.; Liu, L.; Feng, G.; Peng, X.; Wen, D. Orbit error removal in InSAR/MTInSAR with a patch-based polynomial model. Int. J. Appl. Earth Obs. Geoinf. 2021, 102, 102438. [Google Scholar] [CrossRef]
  22. Liang, H.; Zhang, L.; Ding, X.; Lu, Z.; Li, X. Toward Mitigating Stratified Tropospheric Delays in Multitemporal InSAR: A Quadtree Aided Joint Model. IEEE Trans. Geosci. Remote Sens. 2019, 57, 291–303. [Google Scholar] [CrossRef]
  23. Shi, M.; Peng, J.; Chen, X.; Zheng, Y.; Yang, H.; Su, Y.; Wang, G.; Wang, W. An Improved Method for InSAR Atmospheric Phase Correction in Mountainous Areas. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2021, 14, 10509–10519. [Google Scholar] [CrossRef]
  24. Goel, K.; Adam, N.; Shau, R.; Rodriguez-Gonzalez, F. Improving the reference network in wide-area Persistent Scatterer Interferometry for non-urban areas. In Proceedings of the 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Beijing, China, 10–15 July 2016; pp. 1448–1451. [Google Scholar] [CrossRef]
  25. Xue, F.; Lv, X.; Dou, F.; Yun, Y. A Review of Time-Series Interferometric SAR Techniques: A Tutorial for Surface Deformation Analysis. IEEE Geosci. Remote Sens. Mag. 2020, 8, 22–42. [Google Scholar] [CrossRef]
  26. Werner, C.; Wegmuller, U.; Strozzi, T.; Wiesmann, A. Interferometric point target analysis for deformation mapping. In Proceedings of the IGARSS 2003 IEEE International Geoscience and Remote Sensing Symposium, Toulouse, France, 21–25 July 2003; Volume 4367, pp. 4362–4364. [Google Scholar] [CrossRef]
  27. Li, T.; Liu, G.; Lin, H.; Jia, H.; Zhang, R.; Yu, B.; Luo, Q. A Hierarchical Multi-Temporal InSAR Method for Increasing the Spatial Density of Deformation Measurements. Remote Sens. 2014, 6, 3349–3368. [Google Scholar] [CrossRef]
  28. Hooper, A. A multi-temporal InSAR method incorporating both persistent scatterer and small baseline approaches. Geophys. Res. Lett. 2008, 35, 16. [Google Scholar] [CrossRef]
  29. Hou, J.; Xu, B.; Li, Z.; Zhu, Y.; Feng, G. Block PS-InSAR ground deformation estimation for large-scale areas based on network adjustment. J. Geod. 2021, 95, 111. [Google Scholar] [CrossRef]
  30. Ge, D. Research on Key Technologies for Regional Ground Subsidence InSAR Monitoring. Ph.D. Thesis, China University of Geosciences, Beijing, China, 2013. [Google Scholar]
  31. Liu, H.; Yu, J.; Zhao, Z.; Jezek, K.C. Calibrating and mosaicking surface velocity measurements from interferometric SAR data with a simultaneous least-squares adjustment approach. Int. J. Remote Sens. 2007, 28, 1217–1230. [Google Scholar] [CrossRef]
  32. Hu, J. Theory and Methodology of InSAR Three-Dimensional Deformation Estimation Based on Modern Measurement Leveling. Ph.D. Thesis, Central South University, Changsha, China, 2013. [Google Scholar]
  33. Kalia, A.C.; Frei, M.; Lege, T. A Copernicus downstream-service for the nationwide monitoring of surface displacements in Germany. Remote Sens. Environ. 2017, 202, 234–249. [Google Scholar] [CrossRef]
  34. Murray, K.D.; Lohman, R.B.; Bekaert, D.P.S. Cluster-Based Empirical Tropospheric Corrections Applied to InSAR Time Series Analysis. IEEE Trans. Geosci. Remote Sens. 2021, 59, 2204–2212. [Google Scholar] [CrossRef]
  35. Li, Z.; Cao, Y.; Wei, J.; Duan, M.; Wu, L.; Hou, J.; Zhu, J. Time-series InSAR ground deformation monitoring: Atmospheric delay modeling and estimating. Earth-Sci. Rev. 2019, 192, 258–284. [Google Scholar] [CrossRef]
  36. Wang, Y.; Feng, G.; Li, Z.; Luo, S.; Wang, H.; Xiong, Z.; Zhu, J.; Hu, J. A Strategy for Variable-Scale InSAR Deformation Monitoring in a Wide Area: A Case Study in the Turpan–Hami Basin, China. Remote Sens. 2022, 14, 3832. [Google Scholar] [CrossRef]
  37. Luo, S.; Feng, G.; Xiong, Z.; Wang, H.; Zhao, Y.; Li, K.; Deng, K.; Wang, Y. An Improved Method for Automatic Identification and Assessment of Potential Geohazards Based on MT-InSAR Measurements. Remote Sens. 2021, 13, 3490. [Google Scholar] [CrossRef]
  38. Zhong, Y.; Dong, S.; Hu, H. Monitoring and Analysis of Ground Settlement in Changzhou City Based on Time-Series InSAR Technology. Geol. J. China Univ. 2019, 25, 131–143. [Google Scholar]
  39. Wang, J.; Wei, Z.; Liu, W.; Zhong, S.; Zhang, G. Functional zoning of land consolidation in mountainous and hilly areas based on “Production-ecological” perspective: A case study of Qijiang District, Chongqing. Areal Res. Dev. 2018, 37, 155–159. [Google Scholar]
  40. 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]
Figure 1. Flowchart of the block MT-InSAR algorithm.
Figure 1. Flowchart of the block MT-InSAR algorithm.
Remotesensing 14 04562 g001
Figure 2. (a) Diagram of overlap ratio and overlap area. w and h are the width and height of the overlap region, respectively, and W and H are the width and height of the image, respectively. So, the overlap ratio is w/W or h/H. The shadow area is the overlap area in the block (red line) with the surrounding blocks, and the overlap area is 1 − (1 − 2 × w/W) × (1 − 2 × h/H). (b) Diagram of the coordinates acquirement of corner points in the overlap region. (c) Diagram of the adjustment model. The circles represent the deformation points. The thicker the circles, the more times the regions are overlapped. v represents the deformation rate of each point.
Figure 2. (a) Diagram of overlap ratio and overlap area. w and h are the width and height of the overlap region, respectively, and W and H are the width and height of the image, respectively. So, the overlap ratio is w/W or h/H. The shadow area is the overlap area in the block (red line) with the surrounding blocks, and the overlap area is 1 − (1 − 2 × w/W) × (1 − 2 × h/H). (b) Diagram of the coordinates acquirement of corner points in the overlap region. (c) Diagram of the adjustment model. The circles represent the deformation points. The thicker the circles, the more times the regions are overlapped. v represents the deformation rate of each point.
Remotesensing 14 04562 g002
Figure 3. Deformation sequence correction diagram.
Figure 3. Deformation sequence correction diagram.
Remotesensing 14 04562 g003
Figure 4. Data coverage for (a) Changzhou City and (b) Qijiang City. The red line is the administrative division boundary. The blue frame is the image coverage after partition and the shaded part is the overlapping area. The gray blocks are not in the study area.
Figure 4. Data coverage for (a) Changzhou City and (b) Qijiang City. The red line is the administrative division boundary. The blue frame is the image coverage after partition and the shaded part is the overlapping area. The gray blocks are not in the study area.
Remotesensing 14 04562 g004
Figure 5. Deformation results of the study areas. The results of Changzhou found by (a) partition method and (b) traditional method. The results of Qijiang found by (c) partition method and (d) traditional method. (ej) Are time series results of the selected points. The red “+” is the reference point. The reference points in (a,c) are virtual reference points after free net adjustment because of there are reference points in each small block before adjustment, and they are the center of gravity of the image coverage. The reference points of (b,d) are the real reference points in data processing.
Figure 5. Deformation results of the study areas. The results of Changzhou found by (a) partition method and (b) traditional method. The results of Qijiang found by (c) partition method and (d) traditional method. (ej) Are time series results of the selected points. The red “+” is the reference point. The reference points in (a,c) are virtual reference points after free net adjustment because of there are reference points in each small block before adjustment, and they are the center of gravity of the image coverage. The reference points of (b,d) are the real reference points in data processing.
Remotesensing 14 04562 g005
Figure 6. (a) Deformation precision of each block after adjustment in Changzhou. (b) The histogram of the adjustment precision of each block after adjustment, (c) the histogram of the difference of the mean value of the overlapping area before and after adjustment of each block.
Figure 6. (a) Deformation precision of each block after adjustment in Changzhou. (b) The histogram of the adjustment precision of each block after adjustment, (c) the histogram of the difference of the mean value of the overlapping area before and after adjustment of each block.
Remotesensing 14 04562 g006
Figure 7. The same as Figure 6, but for Qijiang.
Figure 7. The same as Figure 6, but for Qijiang.
Remotesensing 14 04562 g007
Figure 8. The difference in the homonymy points before and after adjustment in regions A, B, C, and D of Changzhou. (a) shows the location of the four regions. (be) Are the difference of the overlapping areas before adjustment in regions A, B, C, and D, and (fi) are the difference after adjustment. (jm) Shows the statistical histograms of the four regions.
Figure 8. The difference in the homonymy points before and after adjustment in regions A, B, C, and D of Changzhou. (a) shows the location of the four regions. (be) Are the difference of the overlapping areas before adjustment in regions A, B, C, and D, and (fi) are the difference after adjustment. (jm) Shows the statistical histograms of the four regions.
Remotesensing 14 04562 g008
Table 1. The image parameters of the study areas.
Table 1. The image parameters of the study areas.
Study AreaDirectionPathHeadingIncidencePixel Spacing
(Rg × Az)
Num of Images
ChangzhouAscendingT69−12.79°36.65°2.33 × 13.98 m110
QijiangT55−12.65°43.64°2.33 × 13.96 m115
Table 2. Precision of the deformation velocity in Changzhou and Qijiang found by the two methods.
Table 2. Precision of the deformation velocity in Changzhou and Qijiang found by the two methods.
Study AreaStrategyAreaNumber of PointsStd
/(mm/yr)
Mean
/(mm/yr)
Difference
/(mm/yr)
Precision
Improvement
ChangzhouPartitionA564,8673.44.00.25%
B589,4304.1
C593,4794.6
TraditionalA561,0713.74.2
B582,6484.2
C591,1994.6
QijiangPartitionD992,6793.23.30.615%
E916,0983.1
F922,8113.6
TraditionalD986,1633.83.9
E897,8713.9
F923,0174.0
Table 3. Time and memory consumption of traditional processing and block processing with the number of parallel processing 5.
Table 3. Time and memory consumption of traditional processing and block processing with the number of parallel processing 5.
ChangzhouQijiang
TraditionalPartitionTraditionalPartition
Original size (pixels)29,739 × 689229,739 × 689228,104 × 764828,104 × 7648
Partition strategy\6 × 5
~30% overlap
\6 × 5
~25% overlap
Size of block (pixels)\7147 × 1373\6374 × 1574
PlatformCPU: AMD Ryzen 9 5900X 12-Core/RAM:64 G
Multi-look5:1
Average number of points in a block\2,941,200\2,493,200
Total number of points50,941,51251,178,81451,620,12051,883,633
Memory Usage27.2 G1.3 G27.9 G1.3 G
Time of partition\~1 h\~1 h
Time of InSAR processing~20 h~1.2 h~20 h~1.2 h
Time of correction\~0.1 h\~0.1 h
Total time~20 h~8.7 h~20 h~8.7 h
Table 4. Results of using different block overlap ratios.
Table 4. Results of using different block overlap ratios.
Overlap RatioTotal OverlapBlock Size (Amount)Precision
/(mm/yr)
Time of Each BlockTotal Time
Traditional0%29,739 × 6892 (1)A3.720 h20 h
B4.2
C4.6
10%36%7000 × 1500 (20)
4539 × 1500 (5)
A13.51.2 h
0.8 h
6.7 h
B14.2
C14.6
20%64%7000 × 1500 (25)
7000 × 892 (5)
A23.51.2 h
0.7 h
7.8 h
B24.1
C24.6
30%84%7000 × 1500 (30)
5239 × 1500 (6)
A33.41.2 h
0.9 h
10.1 h
B34.1
C34.6
40%96%7000 × 1500 (42)
4539 × 1500 (7)
A43.31.2 h
0.8 h
12.7 h
B44.0
C44.5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Y.; Feng, G.; Feng, Z.; Wang, Y.; Wang, X.; Luo, S.; Zhao, Y.; Lu, H. An MT-InSAR Data Partition Strategy for Sentinel-1A/B TOPS Data. Remote Sens. 2022, 14, 4562. https://doi.org/10.3390/rs14184562

AMA Style

Wang Y, Feng G, Feng Z, Wang Y, Wang X, Luo S, Zhao Y, Lu H. An MT-InSAR Data Partition Strategy for Sentinel-1A/B TOPS Data. Remote Sensing. 2022; 14(18):4562. https://doi.org/10.3390/rs14184562

Chicago/Turabian Style

Wang, Yuexin, Guangcai Feng, Zhixiong Feng, Yuedong Wang, Xiuhua Wang, Shuran Luo, Yinggang Zhao, and Hao Lu. 2022. "An MT-InSAR Data Partition Strategy for Sentinel-1A/B TOPS Data" Remote Sensing 14, no. 18: 4562. https://doi.org/10.3390/rs14184562

APA Style

Wang, Y., Feng, G., Feng, Z., Wang, Y., Wang, X., Luo, S., Zhao, Y., & Lu, H. (2022). An MT-InSAR Data Partition Strategy for Sentinel-1A/B TOPS Data. Remote Sensing, 14(18), 4562. https://doi.org/10.3390/rs14184562

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