Next Article in Journal
Assessment of Rangeland Degradation in New Mexico Using Time Series Segmentation and Residual Trend Analysis (TSS-RESTREND)
Next Article in Special Issue
Ambient Seismic Noise and Microseismicity Monitoring of a Prone-To-Fall Quartzite Tower (Ormea, NW Italy)
Previous Article in Journal
Potential of Sentinel-2 Satellite Images for Monitoring Green Waste Compost and Manure Amendments in Temperate Cropland
Previous Article in Special Issue
Long-Term and Emergency Monitoring of Zhongbao Landslide Using Space-Borne and Ground-Based InSAR
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Quality Pixel Selection Applied for Natural Scenes in GB-SAR Interferometry

1
School of Information and Electronics, Beijing Institute of Technology, Beijing 100081, China
2
Chongqing Key Laboratory of Geological Environment Monitoring and Disaster Early-Warning in Three Gorges Reservoir Area, Chongqing Three Gorges University, Chongqing 404000, China
3
Radar Research Lab, School of Information and Electronics, Beijing Institute of Technology, Beijing 100081, China
4
Beijing Institute of Technology Chongqing Innovation Center, Chongqing 401120, China
5
Key Laboratory of Electronic and Information Technology in Satellite Navigation (Beijing Institute of Technology), Ministry of Education, Beijing 100081, China
6
Advanced Technology Research Institute, Beijing Institute of Technology, Jinan 250300, China
*
Author to whom correspondence should be addressed.
Submission received: 23 March 2021 / Revised: 15 April 2021 / Accepted: 19 April 2021 / Published: 21 April 2021
(This article belongs to the Special Issue Remote Sensing for Landslide Monitoring, Mapping and Modeling)

Abstract

:
Phase analysis based on high-quality pixel (HQP) is crucial to ensure the measurement accuracy of ground-based SAR (GB-SAR). The amplitude dispersion (ADI) criterion has been widely applied to identify pixels with high amplitude stability, i.e., permanent scatterers (PSs), which typically are point-wise scatterers such as stones or man-made structures. However, the PS number in natural scenes is few and limits the GB-SAR applications. This paper proposes an improved method to take HQP selection applied for natural scenes in GB-SAR interferometry. In order to increase the spatial density of HQP for phase measurement, three types of HQPs including PS, quasi-permanent scatter (QPS), and distributed scatter (DS), are selected with different criteria. The ADI method is firstly utilized to take PS selection. To select those pixels with high phase stability but moderate amplitude stability, the temporal phase coherence (TPC) is defined. Those pixels with moderate ADI values and high TPC are selected as QPSs. Then the feasibility of the DS technique is explored. To validate the feasibility of the proposed method, 2370 GB-SAR images of a natural slope are processed. Experimental results prove that the HQP number could be significantly increased while slightly sacrificing phase quality.

Graphical Abstract

1. Introduction

Landslide is one of the most frequent and harmful geological hazards all over the world. Its monitoring and early warning plays an important part in disaster prevention and reduction. Surface deformations normally occur before the macro failure of natural and engineered slopes. Deformation measurement is of great significance to forecast and early warn landslide hazards [1,2]. Spaceborne synthetic aperture radar (SAR) is a powerful and effective remote-sensing technique for surveying large areas (km2). Main limitations are due to its long repetition period in the order of several days or weeks [3,4]. Ground-based SAR (GB-SAR) is a good complementary tool of spaceborne SAR suitable for monitoring deformations of small areas (from a single building to a large slope), with the advantages of all-day, all-weather, high-accuracy, and fast image acquisition [5,6].
With two radar images acquired at different times and from the same position, their interferometric phases could be explored to measure deformations based on the differential interferometry technique [7]. Interferometric phase is the most crucial element that determines the accuracy of deformation measurement. However, due to the influence of several decorrelation sources, such as temporal decorrelation, geometric decorrelation, volume decorrelation, and thermal noise, some pixels with low phase quality in an interferogram are not suitable for phase analysis [8]. An adequate pixel selection method to identify high-quality pixels (HQPs) is essential and important for radar interferometry measurement.
A typical selection criterion is the coherence, which is calculated by applying a boxcar window to spatially average two maps belonging to an interferogram. The coherence criterion is suitable for the analysis of natural scenes with distributed scatterers (DSs), but its averaging operation would easily cause low-quality pixels to be selected, and rarely applied in GB-SAR interferometry [9]. The amplitude dispersion index (ADI) criterion is the most classical one to identify pixels with high amplitude stability, i.e., permanent scatterers (PSs) which could keep highly coherent during a long period. PSs typically are point-wise scatterers which are characterized by high reflectivity, such as stones, man-made structures, and outcrops [10]. This method can be well utilized for urban and rocky scenes, if sufficient images (usually no less than 20) are available.
The absence of strong point-wise PSs is the main obstacle for the PS technique to be applied in natural or vegetational scenes. At the present stage, researches on HQP selection are more focused on the spaceborne SAR rather than GB-SAR applications. The DS technique was proposed by Ferretti et al., where DSs usually correspond to debris areas, non-cultivated land with short vegetation, or desert areas [11]. This method needs to identify statistically homogeneous pixels based on the amplitude similarity. Its feasibility should be explored in GB-SAR interferometry. Hooper et al. proposed to use spatial correlation of interferometric phase to find pixels with low phase variance [12,13]. This method demands an iterative spatial filtering to estimate the noise phase of each pixel, and then adopts the temporal coherence to identify those pixels with high phase stability.
This paper focuses on the HQP selection applied for natural scenes in GB-SAR interferometry, and proposes an improved method. Three types of HQPs including PS, QPS (quasi-permanent scatter), and DS, are selected with different criteria. The conventional ADI method which evaluates a pixel’s amplitude stability is firstly utilized to take PS selection. To select those pixels with high phase stability but moderate amplitude stability, the temporal phase coherence (TPC) is defined. Those pixels with moderate ADI values and high TPC are selected as QPSs. Then the feasibility of the DS technique is explored. With the improved method, 2370 GB-SAR images of a natural slope are processed. Experimental results prove that the HQP number could be significantly increased while slightly sacrificing phase quality.

2. High-Quality Pixel

HQPs for natural scenes refer to those scatterers with high phase stability. Based on the PS concept, the idea is to firstly identify those single-point scatterers with moderate amplitude stability and high phase stability, i.e., QPSs. Then the feasibility of the DS technique is explored in GB-SAR interferometry. Although the temporal coherence of those natural scatterers is typically low, the number of pixels sharing the same statistical behavior can be large enough to make it possible for a few of them to exceed the coherence threshold. Together with the PS, QPS, and DS, the spatial density of the HQP utilized for phase analysis could be significantly increased.

2.1. Permanent Scatterer

PSs refer to those point-wise scatterers which are almost not affected by the spatial and temporal decorrelation over a long time period. The ADI criterion is commonly used by analyzing the amplitude series of a pixel and thresholding on the ADI. The ADI D A of a pixel can be expressed as
D A = σ A / m A ,
where σ A and m A are the standard deviation and mean value of a pixel’s amplitude series, respectively.
A same numerical simulation with the published research is repeated to show the ADI criterion [10]. For a deterministic point-wise scatterer, its signal model is built as z k = g + n k ( k = 1 , , K ) , where g is a pixel’s reflectivity. To better analyze the relationship between the dispersion index D A and the phase standard deviation σ V , g is fixed to 1. n is the noise term, and tends to be a complex circular Gaussian distribution with power σ n 2 for both real and imaginary components. The deviation σ n gradually increases from 0.05 to 0.8. For each value of σ n , 5000 estimates of the ADI values of the amplitude series A = | z | are completed. Thirty ( K = 30 ) data are supposed to be available.
The mean values (solid line) and the dispersion (error bars) of the estimated D ^ A are shown in Figure 1a, together with the values of the phase standard deviation (dotted line) σ V . It can be noted that the ADI D A is a good measure of the phase stability for pixels with high SNR (signal-to-noise ratio). Those pixels exhibiting D A value under a given threshold D T (typically lower than 0.25) can be selected as PSs.
Figure 1b shows the scatter map of the dispersion index D A with the phase std σ V . For points located under the red line of D A = 0.25 , their phase std σ V varies from 0.05 rad to 0.33 rad (95% confidence interval). Since the ADI is a criterion in a statistical sense, some PSs with relatively large phase std σ V > 0.25 are inevitably selected, as shown with the cyan points in the set S 2 . Those brown points in the set S 1 with D A 0.25 and σ V 0.25 are ‘true’ PSs which are stable in both amplitude and phase.
For rocky scenes with lots of deterministic point-wise scatterers, the ‘true’ PS number in the Set S 1 is enough for phase analysis. Therefore, based on the ADI criterion, another phase criterion could be further utilized together to reduce and refine the set S 1 + S 2 to S 1 . However, for natural scenes, even if only the ADI criterion is utilized, the PS number in S 1 + S 2 is still few. To increase the point number suitable for differential interferometry, those points with moderate amplitude stability and high phase stability, as shown with the yellow points in the set S 3 , could be utilized. This type of point is defined as QPS.

2.2. Quasi-Permanent Scatterer

QPS refers to a pixel with moderate amplitude stability and high phase stability, i.e., a moderate ADI D A and a small phase std σ V . The interferometric phase φ Int of a pixel in GB-SAR interferometry could be modeled as:
φ Int =   W [ φ defo + φ atm + φ geom + φ noi ] ,
where φ defo is the deformation phase component, φ atm is the phase component due to the atmospheric effects during image acquisitions, φ geom is a geometrical phase term due to the aperture variation or repositioning error. φ noi is the noise term. W [   ] denotes the phase wrapped function.
In order to evaluate the phase stability of a pixel, the TPC γ TPC is defined [14]. N IN denotes the interferogram number, or the length of a pixel’s phase sequence. The value range of γ TPC is from 0 to 1, and a larger value means a better stability of φ noi . A key problem is how to accurately estimate the noise term from the interferometric phase. The basic solution is to design a filtering framework based on the different temporal and spatial properties of those four phase terms in (2).
γ TPC = | n = 1 N IN exp ( j φ noi ( n ) ) | / N IN ,
A further analysis of the numerical simulation in the Section 2.1 is made. Take interferometry process for each data sequence, and the interferometric phase φ Int of any two neighboring data can be regarded as the noise phase φ noi . Based on (3), the TPC γ TPC of each point is calculated. The mean values (solid line) and the dispersion (error bars) of γ TPC , together with the values of the phase std σ V (dotted line), are shown in Figure 2a. Contrast with Figure 1a, γ TPC decreases with the increasement of the noise std σ n . Figure 2b shows the scatter map of the TPC γ TPC with the phase std σ V . For those pixels with phase std σ V lower than 0.25 rad, i.e., left part of the green dotted line, 99.99% of them are with the TPC γ TPC larger than 0.91. To some extent, the TPC γ TPC is an equivalent criterion with the ADI D A .
To better illuminate the relationship between the TPC γ TPC and the ADI D A , Figure 3a shows the scatter map of γ TPC with D A . By setting two proper thresholds of γ T and D T , those pixels with γ TPC γ T and D A > D T in the set S 3 could be selected as QPSs. For those pixels with the ADI D A lower than 0.25, and belonging to S1 + S2, their TPC values are between 0.91 and 0.99 (90% confidence interval). The mean values (solid line) and the dispersion (error bars) of γ TPC with D A , are shown in Figure 3b. When D A equals to 0.25, the variation range of γ TPC is about from 0.91 to 0.96. Moreover, for those pixels with the TPC γ TPC larger than 0.91, and belonging to the set S1 + S3, 99.99% of them are with the ADI D A lower than 0.45. Therefore, those pixels with the TPC γ TPC larger than 0.91, but with the ADI D A between 0.25 and 0.45, could be selected the QPSs.

2.3. Distributed Scatterer

DSs refer to those scatterers which are affected by temporal, geometrical, and volumetric decorrelation, but their phase quality could be improved with their neighboring homogeneous pixels [15].
Ferretti et al. proposed an algorithm named by DespecKS to identify a pair of statistically homogeneous pixels [11]. To evaluate whether two pixels P 1 and P 2 are homogeneous, the two-sample KS test is utilized. By measuring the maximum difference D K between the cumulative probability distribution functions S K P 1 ( x ) and S K P 2 ( x ) of their amplitude sequences, two pixels are homogeneous if D K is lower than a proper threshold.
D K = K / 2 sup | S K P 1 ( x ) S K P 2 ( x ) | ,   S K ( x ) = { 0     i f   x < x i i / K     i f   x i x < x i + 1 1     i f   x x K
where K is the sequence length, and x 1 , x 2 , , x K is the sorted amplitude sequence of a pixel. Figure 4a shows two curves whose amplitudes randomly vary and cannot determine their similarity. Figure 4b shows their sorted amplitude curves. With KS test, both curves satisfy the same normal distribution, and they can be regarded as a pair of homogeneous pixels.
Take an image-pixel P as an example, and then define an estimation window centered on P. Apply two-sample KS test to count the number of its homogeneous pixels. Setting a proper number threshold, P can be regarded as a DS if its homogeneous pixels are enough. For a DS, its phase quality could be enhanced with its surrounding homogeneous pixels. It should be noted that only those DSs whose phase qualities are good enough could be utilized for phase analysis, and a further optimization is essential.

3. Pixel Selection

To select enough HQP for natural scenes, the key problem is how to accurately select QPS, and explore the feasibility of DS selection in GB-SAR interferometry.

3.1. PS Selection

Given N (no less than 20) GB-SAR images, two different ADI thresholds D T 1 (0.1 to 0.25) and D T 2 (0.25 to 0.45) are firstly utilized to take initial selection. With a small ADI threshold to select PS, the possibility of those PSs with large phase errors could be reduced effectively. Pixels with ADI value D A D T 1 are set as the reference PSs. Pixels with ADI value D T 1 < D A D T 2 are set as the candidate set of the QPS ( CS QPS ). Based on those reference PSs, the candidate set should be refined.

3.2. QPS Selection

With N GB-SAR images, every two consecutive images could be explored to construct an interferogram. Then the length of a pixel’s interferometric phase sequence, or the interferogram number is N IN = N 1 . To calculate a pixel’s TPC, the noise phase φ noi in (2) should be estimated. The first three terms φ defo , φ atm , and φ geom in (2) are highly spatially correlated, while φ noi is spatially uncorrelated. Therefore, a spatial filter could be utilized to filter out the spatial phase components and estimate φ noi .
The reference PSs distribute non-uniformly on the interferogram. They can be divided into a certain number of clusters according to their spatial locations with a K-means clustering algorithm, which is a method commonly utilized to automatically partition a point set into K clusters [16]. As shown in Figure 5, those scatter points are partitioned into 5 clusters which are marked with different colors. The sign ‘x’ denotes each cluster center.
For those PSs belonging to a same cluster, calculate their mean location r ¯ and mean phase φ ¯ Clu . Since the interferometric phase φ Int of each pixel is wrapped, the circular period mean filter is utilized, as shown below.
φ ¯ Clu = i = 1 N Clu arg [ exp ( j φ Int Clu ( i ) ) / d Clu ] / N Clu + arg ( d Clu ) d Clu = i = 1 N Clu exp ( j φ Int Clu ( i ) ) / N Clu
For a cluster, N Clu is the total pixel number belonging to it, and φ Int Clu ( i ) denotes the interferometric phase of its ith pixel. d Clu is the mean value of the complex phase of all the pixels belonging to this cluster. arg [    ] is the function to get the phase angle of a complex value.
For every pixel belonging to CS QPS , its spatial phase component φ Spa should be estimated. A 2-D spatial interpolation could be utilized, such as linear or Kriging interpolation. It should be noted that the complex phase exp ( j φ ¯ Clu ) rather than the wrapped phase φ ¯ Clu should be used for interpolation. Subtracting φ Spa from φ Int , the residual phase is φ noi [17].
This method is commonly utilized to compensate the space-variant phase in GB-SAR interferometry, if only those montionless PSs with no deformations are used for K-means cluster. Research has proved that the geometrical term φ geom caused by repeat-track error in GB-SAR interferometry could be similarly modeled with the atmospheric phase. Undoubtedly, φ geom can be well compensated with the interpolation method. Although the deformation term φ defo cannot be compensated thoroughly, its residual component during a short period for a slow-moving scene could be neglected.
Based on (3), the TPC of a pixel could be estimated with (6). φ Int ( n ) is the nth value of a pixel’s interferometric phase sequence. φ Spa ( n ) is its estimated spatial phase component.
γ TPC = | n = 1 N IN exp ( j φ Int ( n ) j φ Spa ( n ) ) | / N IN ,
Setting a proper TPC threshold γ T , those pixels with γ TPC no smaller than γ T are selected as QPSs. Other pixels with γ TPC smaller than γ T are then selected to be the candidate set of the DS ( CS DS ).

3.3. DS Selection

For a pixel belonging to CS DS , it is used as a window center, and then determined its homogeneous pixel according to whether their amplitude curves satisfy a same distribution. If the number of its homogeneous pixels is large enough, it can be selected as a DS. The key problem is how to improve the phase quality of these DSs, since their coherence is low and cannot be directly utilized for phase analysis.
The coherence matrix Γ ( P ) of a DS P could be estimated as,
Γ ( P ) = E [ d ( P Ω ) d H ( P Ω ) ] 1 N HP P Ω d ( P Ω ) d H ( P Ω ) ,
where Ω is the homogeneous pixel set of a DS P , and N HP denotes the number of its homogeneous pixels. d = exp ( j φ ) is the complex vector of a pixel belonging to Ω , and φ = [ φ 1 ,   φ 2 ,   ,   φ N ] H is its phase vector. H denotes the complex conjugate.
Based on the coherence matrix Γ ( P ) , a method of maximum likelihood estimation could be utilized to extract the optimal phase λ = [ ϑ 1 ,   ϑ 2 , , ϑ N ] T of the DS P , as shown in (7). Λ is a column vector and equals to exp ( j λ ) . · denotes the Hadamard product. A newton iteration method can be used to solve this optimization problem.
λ ^ = argmax { Λ H ( | Γ | 1 · Γ ) Λ } ,
Once the optimal solution λ has been obtained, the quality of the estimated phase values ϑ i should be assessed. Another type of coherence γ DS could be utilized to be a ‘goodness of fit’ measure. Select every DS with γ DS higher than the TPC threshold γ T , and then substitute the original phase values with its optimized values.
γ DS = 1 N 2 N n = 1 N k n N exp [ j ( ϑ n ϑ k ) ] exp [ j ( φ n φ k ) ] ,

3.4. HQP Selection

The processing flow developed for the HQP selection is shown in Figure 6, and it can be presented schematically as follows:
(1)
Two different ADI thresholds D T 1 and D T 2 are utilized for initial selection. Pixels with ADI value D A D T 1 are set as the reference PSs. Pixels with D T 1 < D A D T 2 are set as CS QPS .
(2)
Based on the K-means clustering and spatial filtering, γ TPC of every pixel belonging to CS QPS is calculated. Setting a proper γ T , those pixels with γ TPC no smaller than γ T are selected as QPSs, and other pixels belong to CS DS .
(3)
For every pixel belonging to CS DS , its homogeneous pixels are determined. Define those pixels whose homogeneous pixels are enough as DSs. For each DS, phase optimization with its homogeneous pixels is made. Those DSs with coherence γ DS large than γ T are reserved.
(4)
Together with the PS, QPS, and DS as the HQP, the traditional PSInSAR algorithm could be utilized for differential interferometry.

4. Experimental Information

From 10:39 on 20 April 2019 to 8:08 on 30 April 2019 (China Standard Time), a GB-SAR system was utilized to monitor a vegetational scene (N30°50′39.83″, E108°25′3.58″) located in Wanzhou district, Chongqing City, China. The plane form of this slope is a chair-like shape, as shown in Figure 7. The slope is about 80 m to 108 m long, 120 m to 250 m wide, with an average thickness of 13.33 m, and a volume of about 298,000 m3. It is a typical medium-scale middle-layer soil landslide. A village is located at the slope foot and on the slope. Those residential buildings are strong scatterers, i.e., PSs in the GB-SAR images.
Figure 8 shows the photo of the GB-SAR. The instrument was developed by Beijing Institute of Technology, and provided by Suzhou Institute of Technology Reco Sensing Technology Co., Ltd. As a developed type of high-accuracy deformation measurement system, GB-SAR has been widely utilized to monitor surface deformations of natural or geological slopes.
Table 1 show the system parameters. The system works at Ku band with a wavelength of 1.85 cm. The transmitted signal is a frequency-modulated continuous wave. With the transmitting and receiving antennas moving along a high-precision mechanical track, a large aperture of 1.8 m could be synthesized. For a GB-SAR system with a fixed-length track, the Ku-band wave is conducive to a significant improvement in its azimuth resolution. In addition, it benefits increasing the measurement accuracy under the condition of a same level of noise phase, contrast with a lower-band wave. Although the Ku-band electromagnetic wave cannot penetrate vegetations, those bare soils or ground surfaces make it feasible to select QPS and DS [8].
The experiment lasted for about 10 days, a total of 2370 images were acquired during the monitoring period. The average image acquisition time was about 6 min. Figure 9a shows a power image of the monitoring area at the polar coordinate, which was dB processed with respect to the average noise level, and could be regarded as a signal-to-noise ratio (SNR) map. Take the first N = 30 images as an example. Based on (1), the ADI value of every pixel could be calculated with these N amplitude images. Figure 9b shows the ADI map. Pixels with ADI values lower than 0.25 only occupy about 2.3%.
With the 1st image as the master image, and the 2nd and 30th images as the slave images, two interferograms IM 1 2 and IM 1 30 could be acquired. Figure 10a shows the phase interferogram IM 1 2 , whose temporal baseline is just several minutes. As a result of the atmospheric disturbance, for pixels located on the natural slope, their phases vary smoothly. The atmospheric refractivity affects the propagation speed of the radar wave through the troposphere. Different transmission delays could be inevitably caused due to the variant atmospheric condition, and it could even reach centimeter level at one kilometer. Atmospheric phase is one major error source in GB-SAR measurement, and must be well compensated [18]. Figure 10b shows the interferogram IM 1 30 , whose temporal baseline reaches about three hours. Its phase variation is much more noisy in contrast with that of IM 1 2 , since this vegetational scene is affected by the temporal decorrelation over a time period of several hours.

5. Experimental Results

For the first 30 images acquired from 10:40 to 13:34 on 20 April 2019 (China Standard Time), by setting an ADI threshold of D T 1 = 0.25 , 25,575 reference PSs could be selected. Figure 11a shows the phase map of the reference PSs based on the interferogram IM 1 2 , which sparsely distribute on the whole map. By setting an ADI threshold of D T 2 = 0.45 , 159,753 pixels with ADI values D T 1 < D A D T 2 are selected, and they belong to the set CS QPS . Figure 11b shows the phase map of these pixels, which distribute on most parts of the whole map. A regular phase variation caused by the atmospheric disturbance can be observed. In addition, lots of noisy pixels exist, since the threshold D T 2 is larger than the typical threshold D T 1 utilized for PS selection. This phenomenon illuminates the feasibility of further selecting pixels with high phase quality.
These reference PSs are divided into 70 clusters, as shown in Figure 12a. Those black squares denote the cluster centers, and every polygon marks the boundary of every cluster. With an inverse distance weighting interpolation, the spatial phase components of the pixels belonging to CS QPS are estimated, as shown in Figure 12b. It is worth noting that the estimated phases vary regularly both along the range and azimuth direction. For the phase components along the range direction, they are caused by the atmospheric disturbance at different times. While for the components along the azimuth direction, published researches prove that they are caused by the repeat-track error when radar moves along the mechanical track [19].
After compensating the spatial phase components, the TPC of those pixels belonging to CS QPS is calculated based on (6), as shown in Figure 13a. According to the analysis conclusion in Section 2.2, a TPC threshold of γ T = 0.91 equivalent to the ADI threshold of D T = 0.25 is set. A total number of 20,364 QPSs can be selected, as shown in Figure 13b. Contrast with the PS number, the QPS number occupies about 79.6%. Pixels with the TPC γ TPC lower than the threshold γ T = 0.91 belong to the set CS DS . Figure 13c shows the probability density curves of γ TPC of the reference PSs and CS QPS . The shapes of both curves are much different. This illuminates that with a lower ADI threshold, pixels with high amplitude stability are more possible to be also with high phase stability. Figure 13d shows the scatter map of γ TPC with D A for the reference PSs and CS QPS . Its distribution trend is much similar with the simulation results shown in Figure 3a. The PS number belonging to the set S 2 , i.e., γ TPC < 0.91 and D A 0.25 , is 2251 and makes up 8.8%. Although the phase quality of this type of PSs is doubtful, there is low necessity to further refine them for the natural scenes. The QPS belonging to the set S 3 makes up about 12.75% of the set CS QPS . Therefore, for pixels with high ADI values, it is rather hard to squeeze out enough pixels with high phase stability.
After the QPS selection, the remaining pixels are belonging to the set CS DS . For every pixel, it is used as the center of a 5 × 7 rectangular window, where 5 is the pixel number along the azimuth direction, and 7 is along the range number. The statistically homogeneous pixel of every pixel is determined based on the two-sample KS test. Figure 14a shows the homogeneous pixel map of CS DS . By setting a number threshold of N T = 10 , 111,406 DSs could be selected. Based on (7)–(9), γ DS of these DSs are calculated and shown in Figure 14b. With the coherence threshold γ T = 0.91 , 4260 DSs are reserved. Although the DS number is only 16.7% of the PS number, it still benefits increasing the pixel number utilized for differential interferometry. Figure 14c shows the phase map of the selected DSs, which sparsely distribute on the map. Above all, a total of 50,199 HQPs are selected, including 25,575 PSs, 20,364 QPSs, and 4260 DSs. Figure 14d shows the phase map of the HQPs. Contrast with the PSs selected only with the ADI criterion, the HQPs increase by 96.3%. The spatial density of the pixels for deformation analysis has been significantly increased. Other than the increase in the HQP number, analysis on the phase quality should be further made.
For the first 30 images, 29 interferograms IM 1 2 , IM 2 3 , …, and IM 29 30 could be formed by every two consecutive images. Based on the PS-InSAR algorithm, these interferograms are processed, including phase unwrapping, atmospheric compensation, and deformation measurement. Figure 15a shows the cumulative deformation map. Since these images were acquired during three hours, it was impossible for this slow-moving landslide to displace for several millimeters. Then the standard deviation of every HQP’s deformation is calculated, and utilized to quantitively evaluate the phase quality of the PSs, QPSs, and DSs. Figure 15b shows the probability density curves of these three types of pixels. Most of the PSs are with deformation deviations around 0.3 mm, while for most of the QPSs and DSs, their deformation deviations are both around 0.42 mm. It is possibly caused by that a relatively low coherence threshold of γ T = 0.91 is utilized, compared with the ADI threshold of D T 1 = 0.25 .
Every N = 30 images of all the 2370 images are divided into an image group. For each group, the HQP selection is taken with the same processing for the 1st group. Figure 16a shows the number curves of different types of HQPs for all the 79 groups. It should be noted that the number fluctuations of PS, QPS, and DS for time-series image groups are normal phenomena in GB-SAR applications, which is mainly caused by daily temperature variation, and could be solved by an adaptive selection threshold [7]. The PS numbers of different image groups change between 6355 and 46,920, and its average number is 17,702. In addition, the average numbers of QPS and DS are 10,052 and 3389, and their increase ratios contrast with the average PS number are 56.8% and 19.1%, respectively. Therefore, with the improved method to select the QPS and DS, the HQP number could be significantly increased. For the 12th image group, its HQP number is the most, and reaches 78,935, as shown in Figure 16b. The HQP density is much larger than that of the 1st group.
It is worth noting that the DS numbers in these groups are no more than 6958, and the average DS number is only 19.1% of PS. However, the procedures to take DS selection include homogeneous pixel calculation, phase optimization, and coherence thresholding, which need much computation. Although the DS technique in GB-SAR interferometry is feasible, the DS density is rather small. A lower coherence threshold could be set to select more DSs, but larger deformation deviations are unavoidable. Therefore, the necessity of DS selection in GB-SAR applications is low, unless when the PS number is rather few for a natural slope.
Figure 17 shows the cumulative deformation map during the monitoring period of 10 days. Since the monitored slope is a creep landslide, no large deformation area is found during a relatively short period, but some pixels with measurement errors of several mm are unavoidable. No other instruments with deformation measurement ability were available for the creep slope during the monitoring period. Only qualitative explanations can be made to illustrate the effectiveness of the improved method.

6. Discussion

Pixels with high phase equality are essential and important for radar interferometry measurement. A common method is based on the ADI criterion to select PSs, as shown in Figure 1, which typically refer to deterministic point-wise scatterers, such as rocks. However, for natural scenes, the absence of point-wise PSs could severely affect interferometry measurement. Aiming at increasing the spatial density of HQP applied for natural scenes, an improved method to jointly select PS, QPS, and DS has been proposed.
With 30 GB-SAR images of a natural slope, 25,575 pixels with the ADI values D A lower than 0.25 are firstly selected as the PSs, and those pixels with 0.25 < D A 0.45 are selected as CS QPS , as shown in Figure 11. Then based on these PSs, the K-means clustering and spatial interpolation are utilized to calculate the TPC γ TPC of every pixel belonging to CS QPS . A total of 20,364 pixels with γ TPC 0.91 are selected as QPSs, and those pixels with γ TPC < 0.91 are taken as CS DS , as shown in Figure 13. Lastly, the two-sample KS test is utilized to calculate the homogeneous pixels of every pixel belonging to CS DS , and a number threshold is set to select DS. Based on (7) and (8), γ DS of every DS is estimated. With a coherence threshold of γ T = 0.91 , 4260 DSs are reserved.
To validate the effectiveness of the proposed method, analyses are made from two aspects. As for the 30 images acquired during 3 h, the standard deviation of every HQP’s deformation is calculated. As shown in Figure 15b, most PSs are with deformation deviations around 0.3 mm, while for QPSs and DSs, their deformation deviations are both around 0.42 mm. While slightly sacrificing phase quality, the HQP number has been significantly increased. As for the 2370 images acquired during 10 days, 79 image groups are taken HQP selection. Statistical results in Figure 16a show that the average numbers of PSs, QPSs, and DSs are 17,702, 10,052, and 3389, respectively.

7. Conclusions

This paper proposes an improved method to take HQP selection applied for natural scenes in GB-SAR interferometry. In order to increase the spatial density of HQP for phase measurement, three types of HQPs including PS, QPS, and DS, are selected with different criteria. The conventional ADI method which evaluates a pixel’s amplitude stability is firstly utilized to take PS selection. In order to select those pixels with high phase stability but moderate amplitude stability, a TPC criterion is adopted to take QPS selection. In addition, the feasibility of the DS technique in GB-SAR interferometry is explored. With a GB-SAR system, 2370 images of a natural slope were acquired during an experimental period of 10 days. As for the 30 images acquired during 3 h, statistical results of the deformation deviations prove that the HQP number can be significantly increased, while slightly sacrificing phase quality. As for the 2370 images acquired during 10 days, statistical results of the HQP numbers of 79 groups prove that, contrast with the average PS number, an average increase of 56.8% and 19.1% can be got with the QPS and DS selection.
Although experimental datasets validate the effectiveness of the proposed method, this paper needs more experimental verification. More natural scenes, especially those without residential buildings or bare rocks, should be measured to validate the feasibility to select the HQP. In addition, other measurement techniques, such as GPS, spaceborne SAR, and laser scanner, should be also utilized to measure the same scene, and provide comparative verification with the GB-SAR.

Author Contributions

Conceptualization, Y.D.; methodology, Y.D. and W.T.; validation, T.X. and H.Y.; formal analysis, Y.D. and C.H.; investigation, W.T. and T.X.; resources, H.Y.; writing—original draft preparation, Y.D.; writing—review and editing, T.X.; supervision, C.H.; project administration, W.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported in part by Natural Science Foundation of Chongqing, China, under Grants cstc2020jcyj-jqX0008, in part by National Natural Science Foundation of China under Grants 61960206009, 61971037, 31727901, and in part by Chongqing Key Laboratory of Geological Environment Monitoring and Disaster Early-warning in Three Gorges Reservoir Area under Grant ZD2020A0101.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Casagli, N.; Catani, F.; Ventisette, C.D.; Luzi, G. Monitoring, prediction, and early warning using ground-based radar inter-ferometry. Landslides 2010, 7, 291–301. [Google Scholar] [CrossRef]
  2. Kavoura, K.; Konstantopoulou, M.; Depountis, N.; Sabatakakis, N. Slow-moving landslides: Kinematic analysis and move-ment evolution modeling. Environ. Earth Sci. 2020, 79, 1–11. [Google Scholar] [CrossRef]
  3. Hu, C.; Li, Y.; Dong, X.; Wang, R.; Cui, C. Optimal 3D deformation measuring in inclined geosynchronous orbit SAR differential inter-ferometry. Sci. China Inform. Sci. 2017, 60, 060303. [Google Scholar] [CrossRef]
  4. Zheng, W.; Hu, J.; Zhang, W.; Yang, C.; Li, Z.; Zhu, J. Potential of geosynchronous SAR interferometric measurements in estimating three-dimensional surface displacements. Sci. China Inf. Sci. 2017, 60, 060304. [Google Scholar] [CrossRef]
  5. Caduff, R.; Schlunegger, F.; Kos, A.; Wiesmann, A. A review of terrestrial radar interferometry for measuring surface change in the geosciences. Earth Surf. Process. Landforms 2015, 40, 208–228. [Google Scholar] [CrossRef]
  6. Pieraccini, M.; Miccinesi, L. Ground-based radar interferometry: A bibliographic review. Remote Sens. 2019, 11, 1029. [Google Scholar] [CrossRef] [Green Version]
  7. Hu, C.; Deng, Y.; Tian, W.; Wang, J. A PS processing framework for long-term and real-time GB-SAR monitoring. Int. J. Remote Sens. 2019, 40, 6298–6314. [Google Scholar] [CrossRef]
  8. Tofani, V.; Raspini, F.; Catani, F.; Casagli, N. Persistent Scatterer Interferometry (PSI) Technique for Landslide Characterization and Monitoring. Remote Sens. 2013, 5, 1045–1065. [Google Scholar] [CrossRef] [Green Version]
  9. Iglesias, R.; Mallorqui, J.J.; Lopez-Dekker, P. DInSAR Pixel Selection Based on Sublook Spectral Correlation along Time. IEEE Trans. Geosci. Remote Sens. 2013, 52, 3788–3799. [Google Scholar] [CrossRef]
  10. Ferretti, A.; Claudio, P.; Fabio, R. Permanent scatterers in SAR interferometry. IEEE Trans. Geosci. Remote Sens. 2001, 39, 8–20. [Google Scholar] [CrossRef]
  11. Ferretti, A.; Fumagalli, A.; Novali, F.; Prati, C.; Rocca, F.; Rucci, A. A New Algorithm for Processing Interferometric Data-Stacks: SqueeSAR. IEEE Trans. Geosci. Remote Sens. 2011, 49, 3460–3470. [Google Scholar] [CrossRef]
  12. Hooper, A.; Segall, P.; Zebker, H. Persistent scatterer interferometric synthetic aperture radar for crustal deformation analysis, with application to Volcán Alcedo, Galápagos. J. Geophys. Res. Space Phys. 2007, 112, 1–20. [Google Scholar] [CrossRef] [Green Version]
  13. Hooper, A.; Zebker, H.; Segall, P.; Kampes, B. A new method for measuring deformation on volcanoes and other natural terrains using InSAR persistent scatterers. Geophys. Res. Lett. 2004, 31. [Google Scholar] [CrossRef]
  14. Zhao, F.; Mallorqui, J.J. A Temporal Phase Coherence Estimation Algorithm and Its Application on DInSAR Pixel Selection. IEEE Trans. Geosci. Remote Sens. 2019, 57, 8350–8361. [Google Scholar] [CrossRef] [Green Version]
  15. Jiang, M.; Guarnieri, A.M. Distributed scatterer interferometry with the refinement of spatiotemporal coherence. IEEE Trans. Geosci. Remote Sens. 2020, 58, 3977–3987. [Google Scholar] [CrossRef]
  16. Kanungo, T.; Mount, D.; Netanyahu, N.; Piatko, C.; Silverman, R.; Wu, A. An efficient k-means clustering algorithm: Analysis and implementation. IEEE Trans. Pattern Anal. Mach. Intell. 2002, 24, 881–892. [Google Scholar] [CrossRef]
  17. Hu, C.; Deng, Y.; Tian, W.; Zhao, Z. A compensation method for a time–Space variant atmospheric phase applied to time-series GB-SAR images. Remote Sens. 2019, 11, 2350. [Google Scholar] [CrossRef] [Green Version]
  18. Karunathilake, A.; Sato, M. Atmospheric phase compensation in extreme weather conditions for ground-based SAR. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2020, 13, 3806–3815. [Google Scholar] [CrossRef]
  19. Hu, C.; Zhu, M.; Zeng, T.; Tian, W.; Mao, C. High-Precision deformation monitoring algorithm for GBSAR system: Rail de-termination phase error compensation. Sci. China Inf. Sci. 2015, 58, 082307. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Numerical simulation results. (a) Variation curves of the estimated ADI (amplitude dispersion index) D A and phase std σ V with the noise deviation σ n ; (b) scatter map of D A with σ V .
Figure 1. Numerical simulation results. (a) Variation curves of the estimated ADI (amplitude dispersion index) D A and phase std σ V with the noise deviation σ n ; (b) scatter map of D A with σ V .
Remotesensing 13 01617 g001
Figure 2. Numerical simulation results. (a) Variation curves of the temporal phase coherence (TPC) γ TPC and phase std σ V with the noise deviation σ n ; (b) scatter map of γ TPC with σ V .
Figure 2. Numerical simulation results. (a) Variation curves of the temporal phase coherence (TPC) γ TPC and phase std σ V with the noise deviation σ n ; (b) scatter map of γ TPC with σ V .
Remotesensing 13 01617 g002
Figure 3. Numerical simulation results. (a) Scatter map of the TPC γ TPC with the ADI D A ; (b) variation curves of γ TPC with D A .
Figure 3. Numerical simulation results. (a) Scatter map of the TPC γ TPC with the ADI D A ; (b) variation curves of γ TPC with D A .
Remotesensing 13 01617 g003
Figure 4. Numerical simulation. (a) Random amplitude curves; (b) sorted amplitude curves.
Figure 4. Numerical simulation. (a) Random amplitude curves; (b) sorted amplitude curves.
Remotesensing 13 01617 g004
Figure 5. Simulation diagram of cluster partition.
Figure 5. Simulation diagram of cluster partition.
Remotesensing 13 01617 g005
Figure 6. Processing flow for the high-quality pixel (HQP) selection.
Figure 6. Processing flow for the high-quality pixel (HQP) selection.
Remotesensing 13 01617 g006
Figure 7. Scene photo.
Figure 7. Scene photo.
Remotesensing 13 01617 g007
Figure 8. Ground-based synthetic aperture radar (GB-SAR) system.
Figure 8. Ground-based synthetic aperture radar (GB-SAR) system.
Remotesensing 13 01617 g008
Figure 9. (a) Signal-to-noise ratio (SNR) map of the monitoring area; (b) ADI map calculated with the first 30 images.
Figure 9. (a) Signal-to-noise ratio (SNR) map of the monitoring area; (b) ADI map calculated with the first 30 images.
Remotesensing 13 01617 g009
Figure 10. (a) Phase interferogram IM 1 2 ; (b) phase interferogram IM 1 30 .
Figure 10. (a) Phase interferogram IM 1 2 ; (b) phase interferogram IM 1 30 .
Remotesensing 13 01617 g010
Figure 11. (a) Phase map of the reference permanent scatterers (PSs) based on IM 1 2 ; (b) phase map of the pixels belonging to CS QPS (candidate set of the quasi-permanent scatterers).
Figure 11. (a) Phase map of the reference permanent scatterers (PSs) based on IM 1 2 ; (b) phase map of the pixels belonging to CS QPS (candidate set of the quasi-permanent scatterers).
Remotesensing 13 01617 g011
Figure 12. (a) Cluster division map; (b) estimated phase map of the pixels belonging to CS QPS .
Figure 12. (a) Cluster division map; (b) estimated phase map of the pixels belonging to CS QPS .
Remotesensing 13 01617 g012
Figure 13. (a) TPC map of the CS QPS ; (b) phase map of the QPSs (quasi-permanent scatterers) selected with the threshold of γ T = 0.91 . (c) Probability density curves of γ TPC of the reference PSs and CS QPS ; (d) scatter map of γ TPC with D A for the PSs and CS QPS .
Figure 13. (a) TPC map of the CS QPS ; (b) phase map of the QPSs (quasi-permanent scatterers) selected with the threshold of γ T = 0.91 . (c) Probability density curves of γ TPC of the reference PSs and CS QPS ; (d) scatter map of γ TPC with D A for the PSs and CS QPS .
Remotesensing 13 01617 g013
Figure 14. (a) Statistically homogeneous pixel map of CS DS (candidate set of the distributed scatterers); (b) coherence γ DS map of the distributed scatterers (DSs); (c) phase map of the selected DSs after coherence selection; (d) wrapped phase map of the HQPs based on IM 1 2 .
Figure 14. (a) Statistically homogeneous pixel map of CS DS (candidate set of the distributed scatterers); (b) coherence γ DS map of the distributed scatterers (DSs); (c) phase map of the selected DSs after coherence selection; (d) wrapped phase map of the HQPs based on IM 1 2 .
Remotesensing 13 01617 g014
Figure 15. (a) Cumulative deformation map of the first 30 images; (b) probability density curves of the deformation std of three types of pixels: PS, QPS, and DS.
Figure 15. (a) Cumulative deformation map of the first 30 images; (b) probability density curves of the deformation std of three types of pixels: PS, QPS, and DS.
Remotesensing 13 01617 g015
Figure 16. (a) Number curves of different types of HQPs for all the image groups. (b) Interferometric phase map of the 12th group with the most HQPs.
Figure 16. (a) Number curves of different types of HQPs for all the image groups. (b) Interferometric phase map of the 12th group with the most HQPs.
Remotesensing 13 01617 g016
Figure 17. Cumulative deformation map during the monitoring period of 10 days.
Figure 17. Cumulative deformation map during the monitoring period of 10 days.
Remotesensing 13 01617 g017
Table 1. Parameters of the GB-SAR system.
Table 1. Parameters of the GB-SAR system.
ParameterValue
central frequency16.2 GHz
wavelength1.85 cm
synthetic aperture1.8 m
bandwidth600 MHz
resolution0.25 m × 5 mrad @1 km
measuring range50 m to 5 km
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Deng, Y.; Tian, W.; Xiao, T.; Hu, C.; Yang, H. High-Quality Pixel Selection Applied for Natural Scenes in GB-SAR Interferometry. Remote Sens. 2021, 13, 1617. https://doi.org/10.3390/rs13091617

AMA Style

Deng Y, Tian W, Xiao T, Hu C, Yang H. High-Quality Pixel Selection Applied for Natural Scenes in GB-SAR Interferometry. Remote Sensing. 2021; 13(9):1617. https://doi.org/10.3390/rs13091617

Chicago/Turabian Style

Deng, Yunkai, Weiming Tian, Ting Xiao, Cheng Hu, and Hong Yang. 2021. "High-Quality Pixel Selection Applied for Natural Scenes in GB-SAR Interferometry" Remote Sensing 13, no. 9: 1617. https://doi.org/10.3390/rs13091617

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