From the view of statistical signal processing, the phase denoising is to retrieve the true value from the noisy observations, i.e. interferometric pair of SAR images. The statistical model of interferometric data and noise model of interferometric phase construct the basic foundation of InSAR technology. In order to describe the existing methods, it is necessary to introduce the signal model beforehand.
A. Statistical model
Let an interferometric pair of single look complex (SLC)
SAR images be defined in a vector form as
where 1z and 2z are the obtained SAR images of two
channels, and T denotes the transpose operator. Usually, k
can be assumed to follow a bivariate complex Gaussian distribution with zero means and its probability density
function (PDF) is expressed as [25]
where H is the transpose and conjugate operator, and is
a 2 2
where
reflectivities, and is the correlation coefficient between
and
where is the so-called coherence with a real value, and
is the actual/noise-free interferometric phase. In (4),
multiplication between two SAR images.
From (2) and (3), the InSAR observations, i.e.
under the constraint of and the PDF is expressed
as [27]
where . Then, the phase distribution can be obtained by integrating 1a and 2a out, which is expressed
as [28]
The curve of versus different values of coherence is plotted in Fig.1 (a). As the increase of
, the distribution of
is much more compact with smaller noise turbulence, showing a higher quality of interferometric phase. In fact, all the noise sources aforementioned in the introduction can be treated as the decorrelation issues,
reflecting from the decrease of coherence. As a result, lower coherence means presence of stronger noise while higher coherence indicates smaller noise level. It is well known that the phase noise level of SAR interferogram is spatial-variant, which is highly related to reflectivity and tomography of the observed scenes. Fig.2 (a) and (b) show the interferometric phase and estimated coherence using the simulated interferometric data of Mount Etna, which is used to illustrate this phenomenon.
With the InSAR development, the distribution of interferometric phase has been further completely investigated. Several researchers have independently derived the PDF of the multilook interferometric phase, which is
expressed as [11]
where L is the look number, is the Gamma function, and
the interferometric phase distribution is unimodal, symmetric
and modulus 2 with its peak at the location of actual
interferometric phase
noise between coherence and look number is discovered and the curve of phase standard deviation with respect to coherence for specific look numbers is shown as Fig.1 (b) [11]. Moreover, the Cramer-Rao Bound (CRB) of the standard
deviation of phase is also conducted as [11, 28]
and look number L is helpful to reduce the phase noise. For multilook filter, the increase of look number means increasing the used window size to include more available samples, assuming that all the samples in the window are stationary and homogeneous to follow i.i.d. However, this assumption is not always true for the scenes with high topography and heterogeneity. Therefore, the implementation of phase denoising needs to consider the InSAR phase characteristics.
Fig.1 Statistics of interferometric phase. (a) Probability density function ( ), (b) standard deviation of noisy phase.
Fig.2 Simulated interferometric data of Mount Etna. (a) Noisy interferometric phase, (b) the coherence of (a).
B. Phase noise model
In the real phase domain, the interferometric phase noise
can be modelled as additive noise [13]
where v is the additive noise with zero mean. In (9),
important characteristic for phase denoising. In this case, the noise standard deviation is denoted as and can be estimated from using a lookup table shown as Fig.1 (b). As has also been aforementioned in section II-A, v depends on the factors of coherence
and look number L .
Due to the phase warped property,
phase unwrapping is a necessary step for the noise reduction in the real phase domain. Alternatively, phase denoising using the complex data is more effective and powerful. The noisy phase
in the complex domain is modelled as [29]
where
phase, noise-free phase and phase noise in the complex
domain, respectively. In (10),
determined by coherence , which can also be used to
indicate the phase quality. The real and imaginary parts of
are treated as zero-mean additive noise, which can be modelled as independence from the actual phase . In a similar way,
also depends on the terms of coherence
and look number L .
Last but not the least, phase denoising can also be realized from interferogram estimation/filtering, together with the de-speckling processing. As a foundation of interferogram
denoising, the signal model of interferogram
[21]
where
signal-dependent noise (the dependence arises from the speckle noise model, which is not focused in this paper). Under the constraint of , reduces to be e
while the corresponding likelihood function is denoted as (6).
Based on the aforementioned interferometric data model, a
variety of InSAR phase denoising methods have been proposed in the literature. According to the history of algorithm development and technology feature, the InSAR phase denoising methods can mainly be categorized into four categories: traditional local spatial-domain filters, transformed-domain filters, NL filters and newly advanced methods. The first two categories are the groups of traditional and mature methods in this field. Fig.3 shows the statistics for the journal and conference publications of the two group methods (from web of science). It can be seen from Fig.3 that there is no significant growth in the number of publications since 2002, indicating the relatively mature research status. The NL filter can be treated as a new generation of the phase denoising technology in this field. Fig.4 shows the publication statistics of NL filters that there is significant growth since 2009 (from web of science). Moreover, as one typical advanced technique, the publication statistics of sparse methods are also shown as Fig.5 (from web of science). As a very new technology, the study on sparse technique is very limited, which needs more attention to evaluate the potentials. In the following, a non-exhaustive review is attempted to classify the existing methods, analyzing the strength and limitation of different technologies.
Fig.3 The statistics for total publications of local filters and transformed-domain filters.
Fig.4 The statistics for publications of nonlocal filters.
Fig.5 The statistics for publications of sparse methods.
As a mature research field, the local filter methods have a long history and the multilook filter [11, 30-32] can be treated as one of the very earliest methods. The basic principle of this group methods employs local window with pixels connected to select stationary and homogenous samples for phase estimation. In this section, the introduction of local filters is divided into four groups. For clarity, the processing flowchart is shown as Fig.6.
A. Median/mean filters
This group methods utilize the mathematical morphologies of mean [19, 33-36] and median [16-18, 37-39] functions to estimate the phase, assuming that the local phase in a small window is statistically stationary and homogenous to follow the i.i.d. assumption. In this view, the expectation of pivoting mean/median filter is to obtain optimal filtered results from local statistics. The formulation of pivoting mean filter is
expressed as
where
is an average operator. Meanwhile, the formulation
of median filter is given by
where is a median operator. The preconditions of both the filters highly depends on the stationary and homogenous property of local phase that the topography fluctuation in the local region is slow for the sampling rate and the adjacent pixels have high correlation. When the noise is additive Gaussian noise, the pivoting mean filter is statistically optimal from a view of ML estimation. However, the stationary and homogenous assumption conflicts with the fast variation of topography. Similar to the multilook filter, the mean filter unfortunately has over-smoothing effect and can hardly deal well with the phase details of high topography variation. Compared with the mean filter, the pivoting median form has a better preservation of phase fringe, but the power of noise suppression reduces in some degree. Besides, one major limitation of this group methods is non-adaptive regardless of the local spatial-variation of the selected samples, such as noise level, using an equal average. Several modifications have been designed with improvement by distinguishing the contributions between different samples, such as weighted average [36, 39].
B. Lee filter and its refined forms
The basic principle of this group methods uses locally directional window and applies Lee filter [13] or refined version [14, 40, 41] with adaptive local noise adjustment. The processing of Lee filter can be divided into two steps. The first step is the stationary and homogenous sample selection by evaluating from candidates of different directional windows. To adaptively capture the phase fringe, Lee filter has predesigned different windows of 16 directions and tries to choose the best association one. The next step is statistical parameter estimation and minimum mean square error (MMSE) estimation. The generalized mathematical formulation of Lee
filter is shown as below
where
coefficient of Lee filter and as aforementioned that can be estimated from
using a lookup table. To calculate
and var( )
needs to be unwrapped in the filtering window.
As an alternative way, the complex formulation of Lee filter is
given by
where
value of
coherence, the value of is very close to 1 and Lee filter
area of low coherence, Lee filter reduces to be a mean filter as
depends on the accuracy of directional window selection and parameter estimation. The major drawback of Lee filter is the number limitation of predefined window and sensitivity to strong noise, which has limitation when facing with the phase fringe of complicated structures and textures. In the area of high topography, Lee filter does not work well and tends to introduce some undesired artifacts, such as discontinuity. Several modifications have been designed to Lee filter with improvements, such as local adaptive filter [15, 40] and refined Lee filter [14]. To overcome the number limitation of window, the strategy of adaptively estimating the directional window is employed by retrieving the local frequency to determine the direction of phase fringe [15]. The refiled Lee filter proposed in [14] has done some modifications to improve the robustness of directional window selection and statistical parameter estimation.
C. Local frequency estimators
The local frequency estimators are based on the polynomial phase model, such as linear [42-45] and nonlinear [46-49]. The basic principle of this group methods estimates the local frequency of the interferometric phase to derive the phase fringe. Usually, the local phase in a small window is modelled
as sine wave with one major component, which is shown as
where
at site ( , )k l , ,and ( , )
are the coefficient and two-dimensional (2-D) frequency, respectively. Then, ( , )
can be estimated using the spectrum estimators, such as ML [42], vector covariance matrix (VCM) [42], and multiple-signal classification (MUSIC) [45]. Based on the estimated fringe frequencies, the filtered phase can be
estimated as
where
the linear frequency model is not accurate enough by neglecting other frequency components. To overcome this limitation, several innovative strategies, such as iteration with varying the window size [18] and nonlinear phase modelling [47-49], are applied to preserve the phase details of fringe and texture. The order of model selection is crucial to this group methods. The higher the order choose, the more accurate the phase estimation. However, the cost is the increase of computational load. As a result, it needs to balance between order selection and computational complexity. Besides, the integration methods of phase frequency estimation and unwrapping have also been studied [51-57], which is not introduced here.
D. Hybrid methods
The procedure of this group methods is usually decomposed into two step estimations. The basic idea is that the interferometric phase can be separated into the principal and residual components, such as low-resolution (LR) fringe and high-resolution (HR) pattern [18, 41, 60, 61]. The principal component can usually be extracted using the technique of local frequency estimation and so on. After principal phase compensation, the residual component of interferometric phase can be estimated more precisely. It is obvious that the phase is
estimated in a hybrid way, which is shown as
where
residual phase, respectively. For the local filters, the sample selection strategy, such as the window form, is critical to the phase/frequency estimation. Until now, the developed window can be roughly classified as oriented/directional window [13-15], adaptive size-varying window [18, 43] and region
growing mask [59, 60]. It is highlighted that the region
Fig.6 The processing flowchart of the local filters.
Fig.7 Filtered phase images using local filters. (a) Pivoting mean filter [19] (residue number 2831), (b) pivoting median filter [39] (residue number 2826), (c) Lee filter [13] (residue number 6161), (d) local frequency estimator using ML method [42] (residue number 2409), (e) region growing method [60] (residue number 2189).
Fig.8 The processing flowchart of the transformed-domain filters.
Fig.9 Filtered phase images using the transformed-domain filters. The first and second rows are Goldstein filter [62] and wavelet filter [29], respectively.
Fig.10 The processing flowchart of the nonlocal filters.
growing (named as intensity-driven adaptive neighborhood, IDAN) method [60] employs local window with adaptive shape and size to search for the similar pixels of being stationary and homogenous. Among the traditional local filters, the region growing method can select the most relevant samples with high quality, which can be treated as one best of the established sample selectors.
E. Comparison of local filters
In this subsection, numerical experiments using the aforementioned interferometric data are performed to compare some typical local filters, i.e. pivoting mean filter [19], pivoting median filter [39], Lee filter [13], ML method [42], and region growing method [60]. Fig.7 shows the filtered images that the rows are the filtered phase, phase difference (between filtered and noisy ones), the residue and coherence, respectively. It can be seen that the pivoting mean and median filters tend to introduce phase discontinuity in the
highly-sloped area, producing many residues. The ML method
has a better performance in the high topography regions. However, this method will produce excessive filtering, which will lead to loss of phase details. Besides, the size of the filter window highly affects the ML results. In comparison with the three aforementioned filters, Lee filter has a better performance on noise reduction and detail preservation. However, Lee filter still does not work well by introducing artifacts in the area of very low coherence. Among these filters, the performance of region growing method is nearly the best. However, the improvement of phase estimation in the highly-sloped regions is still very limited. It can be concluded that the traditional local filters are not very effective to deal with the high topography variation, which is an inherent limitation of the local filters.
Instead of filtering in the spatial domain, the InSAR phase
denoising can also be performed in the transformed domain. The transformed-domain filters have some superiorities because the phase signatures can be enhanced to be more separable from the noise in the transformed domain. In this section, we introduce the transformed-domain filters and divide them into two groups. For clarity, the processing flowchart of the transformed-domain filters is shown as Fig. 8.
A. Frequency domain methods
The frequency domain usually refers to the Fourier transform and the discrete cosine transform. The most typical and effective one is the Goldstein filter [62] and its refined versions [63-74]. The basic idea lies into the fact that the interferometric phase is locally stationary and homogeneous and most of the useful components of local phase is usually limited to a narrow band while the noise is distributed in the entire frequency band. The Goldstein filter transfers the interferometric phase image into frequency domain and weights the phase spectrum for noise reduction. The
generalized formulation of Goldstein filter is expressed as
where
phase image in a local window centered at site
represents the smoothing operation, and is the filtered parameter. In fact, proper selection of
is very important for the Goldstein filter. The performance of Goldstein filter is greatly affected by the patch size of local window and filter parameter
. When , it means no filter power without any processing. When 1
, the Goldstein filter has the strongest filtering power. To improve the robustness, Baran [63] has defined the filtering parameter as 1
, meaning that the coherence reflects the noise level to adaptively control the filtering strength. Furthermore, several improvements have been achieved by providing more moderate estimation of the filtering parameter [67], adaptive window size [65], topography compensation [64] and etc [70-74].
B. Wavelet domain filters
The wavelet transform can enhance the signatures of interferometric phase using the multiresolution representation. The multiresolution property has the advantages of reliable
noise separation and perfect nonstationary signal analysis. Lopez and Fabregas propose to apply wavelet transform for phase denoising by modelling the phase noise in the wavelet domain [29], which provides a foundation of this group methods. Based on the assumption of linear wavelet transform, the additive phase noise model in (10) can be easily generated to the wavelet domain. The wavelet transform of (10) can be
expressed as
where and
are the wavelet transforms of
,
and
, respectively. In the i th scale of 2-D wavelet
transform, the detailed form of
where represents the phase information at site
is not change before and after the wavelet transform. When
one deep scale of wavelet transform is applied, the wavelet transform multiplies the phase term by a factor of two without altering the noise level [29]. The behavior of wavelet transform can effectively locate, preserve and amplify the actual phase in the wavelet domain. In this way, the phase and noise can be more easily separated. Then, the phase noise reduction is performed by means of appropriate nonlinear processing, such as wavelet shrinkage or filtering, which can
be expressed in a generalized formulation as
where
the inverse wavelet transform is applied to
image can be obtained with noise reduction. The performance of wavelet filter depends on the scale of wavelet decomposition and the shrinkage/filter function used in the wavelet domain. For better multiresolution capability, the undecimated wavelet transform [76] and wavelet packets [77, 78] are used to improve the phase denoising performance. It is commonly accepted that the wavelet-transformed filters have good preservation of spatial resolution due to the perfect representation of nonstationary signals from multiresolution property.
C. Comparison of transformed-domain filters
In this subsection, the experiments using the same dataset are carried out to show the performances of transformed-domain filers. Fig.9 shows the results using the Goldstein and wavelet filters (residue numbers are 18620 and 25289), respectively. For the Goldstein filter [62], the parameter is set as with a strong filtering effect. As Fig.9 shows that the Goldstein filter has a good phase denoising performance in the high coherence area. However, there are so many residues and large phase differences in the low coherence area, indicating the phase distortion in the high topography areas. For the wavelet filter [29], the dual-tree complex wavelet transform is used for excellent phase preservation. It can be seen from Fig.9 that, in comparison with the Goldstein filter, the wavelet filter has a better performance of phase detail preservation because the wavelet coefficients of noise are still remained with only amplifying the signal components.
As a new generation, the technology of NL filter has emerged high success in InSAR phase denoising, which has been a hot research topic now. The paradigm of NL filter lies into the sample selection strategy of measuring the patch similarity without restricting the search window in a local region, overcoming the limitation of neighbor connection. In the image processing, the concept of NL filter can be tracked to the pioneering work of Lee [79] and the following seminal work of Buades [80, 81] has developed the patch-based filter framework. With the technology development, the NL filter has been attempted into InSAR by establishing new technology, known as NL-InSAR [21], which can be treated as one of the most typical and effective methods. The NL filters employ the patch-wise method to measure the patch similarity in the interferogram domain and weighted average is applied to estimate the filtered interferogram together with phase denoising. In this framework, the most relevant samples can be selected with high quality and phase denoising is achieved by distinguishing the contribution between different samples. It is commonly accepted that the NL filters have the advantage of phase fringe/texture preservation by effectively capturing the phase structure and extracting the phase redundancy in the patch space. In comparison with the traditional local filters, the
NL filters are more power with some unique superiorities, such as better performance in dealing with the high topography and high heterogeneity (such as isolated scenes). For clarity, the processing flowchart of NL filters is shown as Fig. 10.
[80, 81] has established the framework of NL filter by means
of patch-based method and weighted average strategy. To
estimate the noise free interferogram at site
generalized formulation of NL filter is expressed as
where is the weighted coefficient by considering the contribution from the pixel at site t . The key idea of NL filter is that the weight is measured from the similarity
between the patches
(centered at site t ), known as patch-based method. Usually, ,is designed using an exponential kernel of the patch
similarity parameter as
where a large value of means stronger patch similarity to place large weight , and h controls the decay of the exponential function. Essentially, the choice of h is a tradeoff between filtering strength and detail preservation. From (24), a choice of large value h leads to less discriminative weight with stronger filtering strength while a smaller one enhances the patch dissimilarity for better detail preservation.
The design of the weights is an important issue of the NL filter, which is associated to the patch similarity. So it is necessary to properly model the similarity criteria for the NL filters of InSAR, which is crucial to the phase denoising performance. Besides, the transform-domain filter is also helpful for phase denoising. It is very interesting to combine the NL principle in combination with the transform-domain filter, such as wavelet transform. Accordingly, the review of NL filters is divided into two main groups in this section.
A. Probabilistic patch-based filters
The probabilistic patch-based (PPB) filter [21] extends the NL filter to the application of SAR/InSAR using the weighted maximum likelihood estimator (WMLE) [81]. Based on the statistical model of interferometric data in (5), the generalized formulation of PPB filter for interferogram denoising is shown
as [21]
where
pixel at site u and
pixel at site t . The detailed form of InSAR parameters
estimation is given by
For the PPB filter in [21], the patch similarity is determined as: the patch similarity can be treated as measuring how likely
the two patches
Based on the independence assumption of the pixels in each patch, the patch similarity can be computed pixel-wise. The
weight is defined as [21]
where
Accordingly, it also means that the similarity parameter is
computed as [21]
It is also important that the PPB filter employs a strategy of iterative manner. Using the intermediate results, the estimation accuracy of weight computation can be enhanced by using a
posterior probability as
In this way, the weights can be refined at each iteration. The performance of PPB filter relies on the accuracy of parameter estimation, i.e.
topography and scene heterogeneity conflicts with the collection of relevant samples. In this case, the PPB filter may produce over-smoothing in the areas of high topography and high heterogeneity. Several modifications [22, 82-84], such as topography compensation and heterogeneity measurement, have been designed to improve the reliability of the PPB filter in interferogram denoising. Besides, due to the iterative manner of PPB approach, the fast solution of NL filter has also been studied on the GPU platform to consider the practical application [85].
B. Block matching 3-D filters
The InSAR block matching 3-D (BM3D) filter [86] modifies the traditional BM3D filter into the application of InSAR domain. The concept of InSAR BM3D filter can be decomposed into three main steps: 1) use patch-based method to collect the similar patches, and group/stack these similar patches; 2) wavelet transform on the 3-D block data, phase filtering in the wavelet domain and then inverse wavelet transform back to the patch domain; 3) return all the filtered patches to the original locations and aggregate them with overlapped data fusion. The basic principle of noise reduction lies into the fact that the 3-D wavelet transform can enhance the signatures of interferometric phase, providing more reliable noise separation. The formulation of 3-D wavelet
transform on patch block data is expressed as
The phase noise reduction in the wavelet domain is decomposed into two steps, where the first step is the hard
thresholding formulated as
where
threshold, and the estimation of
step using Wiener filtering as
where is the noise variance in the wavelet domain. Therefore, the InSAR BM3D filters have both the superiorities of NL filter and wavelet-domain filter. In addition to the wavelet-domain filter, some other advanced methods have also been proposed, such as pyramidal representation [87], higher order SVD [24] and etc [88-92].
Experimental analysis: As a typical method, the NL-InSAR algorithm proposed in [21] is applied to illustrate the NL filter performance. First, the experiments are performed using the aforementioned data and the filtered images are shown as Fig.11. The search window seize is 15 15, the patch size is
deals well with the smooth area with little disturbance, however, there is artifacts for the filtered phase in the high heterogeneity and high-slopped areas. Then, more experiments are implemented using the RADARSAT-2 repeat-pass interferometric sample data. Due to the resolution limitation,
the search window seize is 11 11, the patch size is 3 3
filter is very powerful of phase and amplitude denoising in the application of building areas.
With the recent development of new concept of signal processing, such as sparse signal processing and machine/deep learning [93, 94], some advanced methods have been proposed to do attempt in InSAR phase denoising. Among these methods, the group of sparse methods [95, 96] has drawn a lot of attentions to show potentials in some degree. In a high probability, this group methods may be treated as good candidate of next technology generation of InSAR phase denoising. In this section, a review of the sparse methods is introduced.
The sparse methods can retrieve the noise-free image by
applying a sparse constraint in a generalized formulation as
where 0 is the 0L -norm of a vector, x / s is the noise-free
image to be recovered, containing the actual phase ,
is
a dictionary or operator to sparsely represent x ,
are the vectors of interferometric pair of SAR images
function by mapping
as a uniform formulation of (9), (10) and (11). Accordingly, the sparse technology can perform noise reduction on real phase, complex phase and interferogram, respectively. In fact, the direction solution of (34) is very difficult because it is a NP-hard problem. Meanwhile, it also depends on the forms of
groups of sparse methods are discussed, considering their
contributions.
A. Sparse regularization
This group of sparse methods employs the Bayesian rules to restore the filtered image through a maximum a posteriori (MAP) estimator [97-100]. Usually, the MAP estimation is transferred to solve an energy minimum problem, which can
be expressed in a generalized formulation as
where is the data fidelity term from the
likelihood function, is the
sparse prior, and is a regularized coefficient to balance the two terms. In (35), can usually have three different formulations according to the noise reduction in different domains. Take the interferogram regularization as an
example,
of (5) as
where x is composed of R and
shown as (4) and (5)). On the other hand, the format of
many different sparse regularization approaches, such as total variation (TV) and wavelet-domain regularizations, which are
listed as follows
(a) (b) (c) (d)
Fig.11 Filtered phase images using the NL-InSAR filter [21] (Mount Etna data). (a) Filtered phase, (b) phase difference, (c) residues, (d) coherence of filtered images.
(a) (b) (c) (d)
Fig.12 Filtered images using the NL-InSAR filter [21] (RADARSAT-2 data). (a) Original phase, (b) filtered phase, (b) phase difference, (d) filtered amplitude.
(a) (b) (c) (d)
Fig.13 Filtered images using the sparse methods (RADARSAT-2 data). TV regularization [23]: (a) filtered phase, (b) filtered amplitude. Wavelet regularization [100]: (c) filtered phase, (d) filtered amplitude.
(a) (b) (c) (d)
Fig.14 Filtered phase images using the sparse wavelet regularization [100].(a) Filtered phase, (b) phase difference, (c) residues (3773), (d) coherence of filtered images.
Fig.15 Filter phase images using methods of different categories (Mount Etna data) [100]. (a) Original, (b) Lee filter [13]. (c) wavelet filter [29], (d) region growing [60], (e) sparse wavelet regularization [100].
1) TV regularization [23, 98]. Usually, the TV regularization applies the gradient operation on the 2-D phase image for
noise reduction, which is expressed as
where
local smooth effect using gradient operation. To improve the robustness of phase reduction, joint TV regularization of phase and amplitude of InSAR data has also been studied, where the
regularization is formulated as
where denotes the pixel of amplitude image R at size u , and
is used to balance the regularization between amplitude and phase. In fact, the sparse TV regularization is not very suited to phase denoising when dealing with the scenes of high topography and high heterogeneity, not satisfying the sparsity assumption on phase gradient.
2) Wavelet-domain regularization [100, 101]. The generalized
formulation of sparse wavelet regularization is given by
where
In this framework, the phase noise reduction is achieved by sparse regularization of phase in the wavelet domain. Besides, the amplitude and phase can be jointly regularized in the wavelet domain to improve the sparsity degree, which is helpful to noise reduction. Accordingly, both the amplitude and phase noise reduction can be achieved using this approach.
B. Sparse coding
Different from sparse regularization using a fixed dictionary, the group methods of sparse coding involve the dictionary learning of to improve the sparse representation of phase
image. The generalized formulation is shown as [102-104]
where is a vector of one patch of InSAR observations (i.e.
of the noise-free image. The processing flowchart of sparse
coding methods can be summarized into three steps: 1) for each pixel, collect the square window patch and arrange it in a vector form; 2) all the patches are used for sparse coding to obtain the filtered results with perfect sparse signal estimation; 3) aggregate the filtered patches to the original locations with data fusion. Finally, the filtered phase image can be obtained
in this manner.
Experimental analysis: Here, two sparse regularization
algorithms are employed to illustrate the performance of
aforementioned RADARSAT-2 data are performed using the
sparse TV [23] and wavelet [100] regularizations. Fig.13
shows the filtered results that both the algorithms can
effectively reduce the noise in the urban areas. Next, the
wavelet regularization is applied on the aforementioned
interferometric data of Mount Etna and the filtered results are
shown as Fig.14.
For clarity, Fig.15 shows more details of the results of the
aforementioned algorithms. As Fig.15 shows, Lee filter
introduces some artifacts of fringes ambiguous and broken in
the areas of high topography and low coherence. In
comparison of wavelet filter, region growing and sparse
interferometric phase noise reduction in smooth areas and
effectively suppress the phase residues of grainy noise. The
sparse wavelet method almost does the best performance in
dealing with blurring and discontinuities of phase fringes.
More details can be referred to [100].
Comments: In addition to the sparse methods, some other advanced technologies have also been tried in InSAR phase denoising, such as random Markov Random Field (MRF) [105-108], tensor decomposition [109], convolutional neural networks (CNN) [110] and optimization integration with phase unwrapping [111-115]. Roughly, most of them can be classified as the machine learning or deep learning methods and the details are not introduced in this paper.
In this paper, we have reviewed the InSAR phase denoising technology. In the past few decades, a large number of methods have been proposed with great achievements, which are classified into four categories in this paper: 1) traditional local filters; 2) transformed-domain filters; 3) NL filters; and (4) newly advanced methods. More importantly, we have introduced the concepts and surveyed the strengths and limitations of these methods, both theoretically and experimentally. At present, the research achievements on InSAR phase denoising look exciting, which have been
successfully applied in the current radar sensors. However, there is still improvement space for the current research in some fields, such as the incoming technology of NL filters and upcoming technology of newly advanced methods. To the best of our knowledge, we try to provide some suggestions on the research tendency, indicating the potential directions in the future.
A. High efficiency and high precision algorithms
The group of NL filters uses the patch-wise based strategy for sample selection and iterative manner to improve the estimation performance. The filter parameters, such as search window size, patch size, similarity measure criteria and iteration number, are crucial to the efficiency and precision of algorithm implementation in practice. In particular, the complicate texture of interferometric images, known as high topography variation and high heterogeneity, is a high challenge in the face of the NL filters, hindering their practical application [116]. Therefore, how to improve the performances of robustness and efficiency is a very important issue for NL filters, even for all the phase denoising methods.
B. Sparse methods
The performance of sparse methods is highly dependent on the model of a sparse prior. The sparse representation of interferometric images is still an open question. It is hard to evaluate the sparse solution of interferometric phase: 1) the sparse results depend on the sparse signal recovery technique, such as sparse regularization and sparse coding; 2) it is difficult to distinguish the interferometric phase from noise in the sparse space because there is no explicit theoretical support. Therefore, how to design a well-established paradigm of sparse phase denoiser is significantly important to promote the practical development of this class technology. Meanwhile, high efficiency is also an important factor.
C. Evaluating phase denoising performance
It is commonly accepted that it is still an open question to objectively evaluate the phase denoising performance. The purpose of phase denoising is to reduce the noise and preserve the phase details as much as possible. In the flat area, noise reduction is the most important issue. In the area of high topography variation, the preservation of phase details without introducing artifacts becomes more important. In addition, the computational complexity is also an important issue. Therefore, it needs further study on how best to comprehensively evaluate the phase denoising performance, which is very helpful to the
architecture development of InSAR signal processing.
Phase denoising technology plays an important role to be a mandatory step in InSAR signal processing. Its ultimate objective is to satisfy the practical demands along with the development of new InSAR systems. The purpose of this paper is to provide necessary guideline to the researchers with some novel inspirations. It is our hope that this paper is helpful to promote the further development of InSAR technology.
We would like to thank the editors and anonymous reviewers for their valuable comments, which have helped to improve the quality of this paper. This work was supported in part by the National Natural Science Foundation of China (NSFC) under Grant 61701106 and Grant 61861136002, in part by the Natural Science Foundation of Jiangsu Province under Grant BK20170698, in part by Shanghai Aerospace Science and Technology Innovation Fund under Grant SAST2018-044. The RADARSAT-2 data were provided by MacDonald, Dettwiler and Associates Ltd online in a public share.
AUTHORS
Gang Xu (M’16) received the B.S. degree and the Ph.D degree in electrical engineering from Xidian University, Xi’an China, in 2009 and 2015, respectively. From 2015 to 2016, he was a full-time Postdoctoral Research Fellow with the School of Electrical and Electronic Engineering (EEE), Nanyang Technological University (NTU), Singapore. He is currently an Associate Professor with the State Key Laboratory of Millimeter Waves, School of Information Science and Engineering, Southeast University, Nanjing, China. His major research interests include high-resolution radar imaging and SAR/ISAR interferometric processing.
Yandong Gao received the B.S. and M.S. degrees in survey and mapping engineering from University of Science and
Technology Liaoning, Anshan, China, in 2013 and 2016, respectively. He received the Ph.D. degree in geodesy and surveying engineering from China University of Mining and Technology, Xuzhou, China, in 2019. He is currently a Postdoctoral Research Fellow with the School of Environment Science and Spatial Informatics, China University of Mining and Technology, Xuzhou, China. His major research interests include the interferometric phase filtering, phase unwrapping, and synthetic aperture radar interferometry signal processing.
Jinwei Li received the M.S. degree and the Ph.D degree in electrical engineering from Xidian University, Xi’an China, in 2009 and 2015, respectively. From 2016, he working as a member at Xi’an institute of space radio technology, Xi’an, China. His major research interests include the fields of interferometric phase filtering, phase unwrapping, and synthetic aperture radar interferometry signal processing.
Mengdao Xing (M’04) was born in Zhejiang, China, in November, 1975. He received the Bachelor degree and the Ph.D degree in electrical engineering from Xidian University, Xi’an, China, in 1997 and 2002, respectively. He is currently a full professor with the National Laboratory of Radar Signal Processing, Xidian University, Xi’an, China. His major research interests include SAR, ISAR and over the horizon radar (OTHR).
[1] L. C. Graham, “Synthetic interferometric radar for topographic mapping,” Proc. IEEE, vol. 62, pp. 763-768, June 1974.
[2] R. Lanari et al., “Generation of digital elevation models by using SIR-C/X-SAR multifrequency two-pass interferometry: The Etna case study,” IEEE Trans. Geosci. Remote Sensing, vol. 34, pp. 1097-1113, Sep. 1996.
[3] J. Moreira et al., “X-SAR interferometry: First results,” IEEE Trans. Geosci. Remote Sensing, vol. 33, pp. 950-956, July 1995.
[4] D. Masonnet, M. Rossi, C. Carmona, F. Adragna, G. Peltzer, K. Feigl, T. Rabaute, “The displacement field of the landers earthquake mapped by radar interferometry,” Nature, vol. 364, no. 6433, pp. 138-142, 1993.
[5] G. Krieger, A. Moreira, H. Fiedler, I. Hajnsek, M. Werner, M. Younis, M. Zink, “TanDEM-X: a satellite formation for high-resolution SAR interferometry,” IEEE Trans. Geosci. Remote Sensing, vol. 45, no. 11, pp. 3317-3341,
2007.
[6] D. Massonnet, P. Briole, A. Arnaud, “Deflation of mount Etna monitored by spaceborne radar interferometry,” Nature, vol. 375, no. 6532, pp. 567-570, 1995.
[7] R. Bamler, P. Hartl. “Synthetic aperture radar interferometry,” Inverse Problems, vol. 14, no.4. pp. R1-R54, 1998.
[8] H. A. Zebker and J. Villasenor, “Decorrelation in interferometric radar echos,” IEEE Trans. Geosci. Remote Sensing, vol. 30, pp. 950-959, July 1992.
[9] R. Goldstein, H. A. Zebker, and C. L. Werner, “Satellite radar interferometry: Two dimensional phase unwrapping,” Radio Sci., vol. 23, pp. 713-720, 1988.
[10] D. C. Ghiglia and M. D. Pritt, Two-dimensional Phase Un-wrapping: Theory, Algorithms, and Software. New York: Wiley, 1998.
[11] J. S. Lee, K. W. Hoppel, S. A. Mango, and A. R. Miller, “Intensity and phase statistics of multilook polarimetric and interferometric SAR imagery,” IEEE Trans. Geosci. Remote Sens., vol. 32, no. 5, pp. 1017-1028, Sep. 1994.
[12] M. S. Seymour and I. G. Cumming, “Maximum likelihood estimation for SAR interferometry,” in Proc. IGARSS, Aug. 8–12, 1994, pp. 2272-2275.
[13] J. S. Lee, K. P. Papathanassiou, T. L. Ainsworth, M. R. Grunes, and A. Reigber, “A new technique for noise filtering of SAR interferometric phase images,” IEEE Trans. Geosci. Remote Sens., vol. 36, no. 5, pp. 1456-1465, Sep.1998.
[14] C. F. Chao, K. S. Chen, and J. S. Lee. “Refined filtering of interferometric phase from InSAR data,” IEEE Trans. Geosci. Remote Sens., vol. 51, no. 12, pp. 5315-5323.
[15] N. Wu, D. Z. Feng, and J. Li, “A locally adaptive filter of interferometric phase images,” IEEE Geosci. Remote Sens. Lett., vol. 3, no. 1, pp. 73-77, Jan. 2006.
[16] A. L. B. Candeias, J. C. Mura, L. V. Dutra, J. R. Moreira, and P. P. Santos, “Interferogram phase noise reduction using morphological and modified median filters,” in Proc. IEEE IGARSS’95, vol. 1, pp. 166-168.
[17] W. B. Abdallah and R. Abdelfattah, “An enhanced weighted median filter for noise reduction in SAR interferograms,” Adv. Concepts Intell. Vis. Syst., vol. 8192, pp. 49-59, 2013.
[18] D. Meng, V. Sethu, E. Ambikairajah, and L. Ge, “A novel technique for noise reduction in InSAR images,” IEEE Geosci. Remote Sens. Lett., vol. 4, no. 2, pp. 226-230, Apr. 2007.
[19] J. Xue, K. Deng, and F. Zhu, “Gauss weighted periodic pivoting filtering algorithm for InSAR interferogram,” Sci. Surv. Mapp., vol. 3, pp. 40-45, 2011.
[20] O. O. Bezvesiniy, I. M. Gorovyi, and Vavriv, D. M, “Effects of local phase errors in multi-look SAR images,” Prog. Electromagn. Res B, vol. 53, no. 53, pp. 1-24, Jan.2013.
[21] C. A. Deledalle, L. Denis, and F. Tupin, “NL-InSAR: Nonlocal interferogram estimation,” IEEE Trans. Geosci. Remote Sens., vol. 49, no. 4, pp. 1441-1452, Apr. 2011.
[22] G. Baier, C. Rossi, M. Lachaise , X. X. Zhu, and R. Bamler, “A nonlocal InSAR filter for high-resolution DEM generation from TanDEM-X interferograms,” IEEE Trans. Geosci. Remote Sens., vol. 56, no. 11, pp. 6469-6483, 2018.
[23] L. Denis, F. Tupin, J. Darbon, and M. Sigelle, “SAR image regularization with fast approximate discrete
minimization,” IEEE Trans. Image Process., vol. 18, no. 7, pp. 1588-1600, Jul. 2009.
[24] X. Lin, F. Li, D. Meng, D. Hu, and C. Ding, “Nonlocal SAR interferometric phase filtering through higher order singular value decomposition,” IEEE Geosci. Remote Sens. Lett., vol. 12, no. 4, pp. 806-810, Apr. 2015
[25] R. Touzi, A. Lopes, J. Bruniquel, and et al, “Coherence estimation for SAR imagery,” IEEE Trans. Geosci. Remote Sens., vol. 37, no. 1, pp. 135-149, 1999.
[26] N. Goodman, “Statistical analysis based on a certain multivariate complex Gaussian distribution (an introduction),” Ann. Math. Stat., vol. 3, no. 1, pp. 152-177, Mar. 1963.
[27] J. Goodman, Speckle Phenomena in Optics: Theory and Applications. Greenwood, CO: Roberts Co. Publishers, 2006.
[28] D. Middleton, Introduction to Statistical Communication Theory. New York: McGraw-Hill, 1960, ch. 9.
[29] C. Lopez-Martinez and X. Fabregas, “Modeling and reduction of SAR interferometric phase noise in the wavelet domain,” IEEE Trans. Geosci. Remote Sens., vol. 40, no. 12, pp. 2553–2566, Dec. 2002.
[30] M. Schmitt, U. Stilla, “Adaptive multilooking of airborne single-pass multi-baseline InSAR stacks,” IEEE Trans Geosci Remote Sens., vol. 52, no. 1, pp. 305–312, Jan.2014.
[31] R. Touzi, A. Lopes, “Statistics of the Stokes parameters and of the complex coherence parameters in one-look and multilook speckle fields,” IEEE Trans. Geosci. Remote Sens., vol. 34, no. 2, pp. 519–531, Apr. 1996.
[32] M. Jiang, X. Ding, R. F. Hanssen, R. Malhotra, and L. Chang, “Fast statistically homogeneous pixel selection for covariance matrix estimation for multitemporal InSAR,” IEEE Trans. Geosci. Remote Sens, vol. 53, no. 3, pp. 1213-1224, 2015.
[33] G. Fornaro and A. Monti Guarnieri, “Minimum mean square error space varying filtering of interferometric SAR data,” IEEE Trans. Geosci. Remote Sens., vol. 40, no. 1, pp. 11–21, Jan. 2002.
[34] T. Zhi, L. Jingwen, and Z. Yinging, “Analysis on noise reduction method for interferometric SAR image,” in Proc. IGARSS, 2004, vol. 6, pp. 4243–4346
[35] X. Ding, Z. Li, D. Zheng, C. Huang, and W. Zou, “SAR interferogram filtering with 2-D Vondrak filter,” in Proc. IGARSS, 2005, vol. 7, pp. 4825–4828.
[36] J. V. Lorenzo-Ginori. K. Plataniotis, A.N. Venetsanopoulos, “Nonlinear filtering for phase image denoising,” IEE Proc-Vision Image and Signal Processing, 2002.
[37] E. Rodriguez, J. M. Martin. Theory and Design of Interferometric Synthetic Aperture Radar [J]. Radar and Signal Processing IEE Proceedings-F, 1992, 139(2): 147-159.
[38] W. B. Abdallah and R. Abdelfattah, “A modification to the ASM filter for improving SAR interferograms,” in Proc. 21st Eur. Signal Process. Conf. EURASIP, 2013, pp. 1–5.
[39] K. Duan, J. Xiang, and F. Wang, “An algorithm of weighted periodic pivoting median filtering for InSAR phase fringe,” J. Air Force Radar Acad., vol. 1, pp. 4–6, 2005.
[40] S. H. Fu, X. J. Long, X. Yang, and Q.F. Yu, “Directionally adaptive filter for synthetic aperture radar interferometric
phase images,” IEEE Trans. Geosci. Remote Sens., vol. 51, no. 1, pp. 552-559, 2013.
[41] Q. S. Wang, H. F. Huang, A. X. Yu, and Z. Dong, “An efficient and adaptive approach for noise filtering of SAR interferometric phase images,” IEEE Geosci. Remote Sens. Lett., vol. 8, no. 6, pp. 1140–1144, 2011.
[42] E. Trouvé, M. Caramma, and H. Maître, “Fringe detection in noisy complex interferograms,” Appl. Opt., vol. 35, no. 20, pp. 3799–3806, Jul. 1996.
[43] Z. Suo, Z. Li, and Z. Bao, “A new strategy to estimate local fringe frequencies for InSAR phase noise reduction,” IEEE Geosci. Remote Sens. Lett., vol. 7, no. 4, pp. 771–775, Oct. 2010.
[44] Y. D. Gao, S. B. Zhang, T. Li, L. Guo, Q. F. Chen, and S. J. Li, “A novel two-step noise reduction approach for interferometric phase images,” Opt. Las. Engine., vol. 121, pp. 1-10, 2019.
[45] Z. F. Li, Z. Bao, H. Li, and G. S. Liao, “Image autocoregistration and InSAR interferogram estimation using joint subspace projection,” IEEE Trans. Geosci. Remote Sens., vol. 44, no. 2, pp. 288–297, Feb. 2006.
[46] Y. D. Gao, S. B. Zhang, K. F. Zhang, and S. J. Li, “Frequency domain filtering SAR interferometric phase noise using the amended matrix pencil model,” Comp Modeling in Eng & Sci.,vol.119, no.2, pp.349-363,2019.
[47] B. Friedlander, J. M. Francos, “An estimation algorithm for 2-D polynomial phase signals,” IEEE Transactions on Image Processing, vol.5, no.6, pp.1084-1087, 1996.
[48] D. B. Ma, M. Liu, Y. Q. Deng, Y. Lin, “A piece-wise polynomial fitting method to filter the Interferogram phase noise,” In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, IGARSS’02, Toronto, Canada, 2002, 3459-3461.
[49] A. Pepe, Y. Yang, M. Manzo, R. Lanari., “Improved EMCF-SBAS processing chain based on advanced technique for the noise-filtering and selection of small baseline multi-look DInSAR interferograms,” IEEE Trans. Geosci. Remote Sensing, 53(8): 4394-4417, 2015.
[50] U. Spagnolini, “2-D phase unwrapping and instantaneous frequency estimation,” IEEE Trans. Geosci. Remote Sensing, vol.33, no.3, 579-589, 1995.
[51] E. Trouve, J. M. Nicolas, and H. Maitre, “Improving phase unwrapping techniques by the use of local frequency estimates,” IEEE Trans. Geosci. Remote Sens., vol. 36, no. 6, pp. 1963–1972, Nov. 1998.
[52] J. Bioucas-Dias, V. Katkovnik, J. Astola, and K. Egiazarian, “Absolute phase estimation: adaptive local denoising and global unwrapping,” Appl. Opt., vol. 47, no. 29, pp. 5358–5369, Oct. 2008.
[53] V. Katkovnik, J. Astola, and K. Egiazarian, “Phase local approximation (PhaseLa) technique for phase unwrap from noisy data,” IEEE Trans. Image Process., vol. 17, no. 6, pp. 833–846, Jun. 2008.
[54] D. Chirico and G. Schirinzi, “Multichannel interferometric SAR phase unwrapping using extended Kalman smoother,” Int. J. Microw. Wireless Technol., vol. 5, pp. 429–436, Jun. 2013.
[55] O. Loffeld, H. Nies, S. Knedlik, and W. Yu, “Phase unwrapping for SAR interferometry-A data fusion approach by Kalman filtering,” IEEE Trans. Geosci. Remote Sens., vol. 46, no. 1, pp. 47–58, Jan. 2008.
[56] X. Xie and Y. Pi, “Phase noise filtering and phase unwrapping method based on unscented Kalman filter,” J.
Syst. Eng. Electron., vol. 22, no. 3, pp. 365–372, 2011.
[57] R. Ambrosino, F. Baselice, G. Ferraioli, and G. Schirinzi, “Extended Kalman filter for multichannel InSAR height reconstruction,” IEEE Trans. Geosci. Remote Sens., vol. 55, no. 10, pp. 5854-5863, 2017.
[58] J. J. Martinez-Espla, T. Martinez-Marin, and J. M. Lopez-Sanchez, “A particle filter approach for InSAR phase filtering and unwrapping,” IEEE Trans. Geosci. Remote Sens., vol. 47, no. 4, pp. 1197–1211, Apr. 2009.
[59] G. Vasile, E. Trouvé, J. S. Lee, and V. Buzuloiu, “Intensity-driven adaptive-neighborhood technique for polarimetric and interferometric SAR parameters estimation,” IEEE Trans. Geosci. Remote Sens., vol. 44, no. 6, pp. 1609–1621, Jun. 2006.
[60] G. Vasile et al., “High-resolution SAR interferometry: Estimation of local frequencies in the context of Alpine glaciers,” IEEE Trans. Geosci. Remote Sens., vol. 46, no. 4, pp. 1079–1090, Apr. 2008.
[61] Q. Yu, X. Yang, S. Fu, X. Liu, and X. Sun, “An adaptive contoured window filter for interferometric synthetic aperture radar,” IEEE Geosci. Remote Sens. Lett., vol. 4, no. 1, pp. 23–26, 2007.
[62] R. M. Goldstein and C. L. Werner, “Radar interferogram filtering for geophysical applications,” Geophys. Res. Lett., vol. 25, no. 21, pp. 4035–4038, Nov. 1998.
[63] I. Baran, M. Stewart, B. Kampes, and P. Lilly, “A modification to the Goldstein radar interferogram filter,” IEEE Trans. Geosci. Remote Sens., vol. 41, no. 9, pp. 2114–2118, Sep. 2003.
[64] Z. G. Ding, Z. Wang, S. Lin, T. D. Liu, Q. Zhang, and T. Long, “Local fringe frequency estimation based on multifrequency InSAR for phase-noise reduction in highly sloped terrain,” IEEE Geosci. Remote Sens. Lett., vol. 14, no. 9, pp. 1527–1531, 2017.
[65] Z. Y. Suo, J. q. Zhang, M. Li, Q. J. Zhang, and C. Fang, “Improved InSAR phase noise filter in frequency domain,” IEEE Trans. Geosci. Remote Sens., vol. 54, no. 2, pp. 1185-1195, 2016.
[66] A. M.-Quereda, J. M. Lopez-Sanchez, J. Selva, and P. J. Gonzalez. “An improved phase filter for differential SAR interferometry based on an iterative method,” IEEE Trans. Geosci. Remote Sens., vol. 56, no. 8, pp. 4477-4491, 2018.
[67] Z. W. Li, X.L. Ding, C. W. M. Huang, Z. J. Jun, Y. L. Chen, “Improved filtering parameter determination for the Goldstein radar interferogram filter,” ISPRS Journal of Photogrammetry and Remote Sensing., vol. 63, no. 6, pp. 621–634, 2008.
[68] C.Y. Zhao, Q. Zhang, X.L. Ding and J. Zhang. “An iterative Goldstein SAR interferogram filter,” International Journal of Remote Sensing., vol.33, no.11, pp.3443-3455, 2012.
[69] M. Jiang, X.L. Ding, X, Tian, R. Malhotra, W. X. Kong. “A hybrid method for optimization of the adaptive Goldstein filter,” ISPRS Journal of Photogrammetry and Remote Sensing., vol. 98, pp. 29-43, 2014.
[70] M. Jiang, X. L. Ding, Z. W. Li, X. Tian, W. Zhu, C. S. Wang, B. Xu, “The improvement for Baran phase filter derived from unbiased InSAR coherence,” IEEE Journal of Selected Topics in Applied Earth Observation and Remote Sensing, vil.7, no.7, pp. 3002-3010, 2014.
[71] R. Song, H. D. Guo, G. Liu, Z. Perski, H. Y. Yue, C. M. Han, J. H. Fan, “Improved Goldstein SAR interferogram
filter based on adaptive-neighborhood technique,” IEEE Geoscience and Remote Sensing Letters, vol.12, no.1, pp.40-144, 2015.
[72] K. M. Qian, “Two-dimensional windowed Fourier transform for fringe pattern analysis: Principles, applications and implementations,” Optics and Lasers Engineering, vol.45, no.2, pp.304-317, 2007.
[73] H. Fattahi, M. J. V. Zoej, M. R. Mobasheri, M. Dehghani, M. R. Sahebi, “Windowed Fourier transform for noise reduction of SAR interferograms,” IEEE Geoscience and Remote Sensing Letters, vol.6, no.3, pp.418-422, 2009.
[74] S. Hensley, "An analytic expression for the phase noise of the goldstein-werner filter," Accepted by IEEE Trans. Geosci. Remote Sensing (in pross.).
[75] W. B. Abdallah, R. Abdelfattah, “Two-dimensional wavelet algorithm for interferometric synthetic aperture radar phase filtering enhancement,” Journal of Applied Remote Sensing, vol.9, no.1, pp. 096061, 2015.
[76] Y. Bian and B. Mercer, “Interferometric SAR phase filtering in the wavelet domain using simultaneous detection and estimation,” IEEE Trans. Geosci. Remote Sens., vol. 49, no. 4, pp. 1396–1416, Apr. 2011.
[77] H. Ru-Yun and W. Yao-Nan, “InSAR interferogram filtering based on wavelet transform,” Acta Geod. Cartographica Sin., vol. 2, pp. 6–16, 2006.
[78] Z. Xianjie, F. Rongshan, D. Zhiyang, and L. Bin, “Noise reduction in interferograms using the wavelet packet transform and Wiener filtering,” IEEE Geosci. Remote Sens. Lett., vol. 5, no. 3, pp. 404–408, Jul. 2008.
[79] J. S. Lee, “Digital image smoothing and the sigma filter,” Comput. Vis. Graph. Image Process., vol. 24, no. 2, pp. 255–269, 1983.
[80] A. Buades, B. Coll, and J. Morel, “A review of image denoising algorithms, with a new one,” Multiscale Model. Simul., vol. 4, no. 2, pp. 490-530, 2005.
[81] A. Buades, B. Coll, and J. Morel, “A non-local algorithm for image denoising,” in Proc. IEEE Comput. Soc. Conf. CVPR, 2005, vol. 2, pp. 60-65.
[82] J. W. Li, Z. F. Li , Z. Bao, et al, “Noise filtering of high-resolution interferograms over vegetation and urban areas with a refined nonlocal filter,” IEEE Geoscience and Remote Sensing Letters, 2015, 12(1):77-81.
[83] F. Sica, , et al. "The offset-compensated nonlocal filtering of Interferometric phase,” Remote Sensing, vol.10., no.9, pp.1359, 2018.
[84] C. A. Deledalle, L. Denis, F. Tupin, A. Reigber, and M. Jäger, “NL-SAR: A unified nonlocal framework for resolution-preserving (Pol)(In)SAR denoising,” IEEE Trans. Geosci. Remote Sens., vol. 53, no. 4, pp. 2021–2038, Apr. 2015.
[85] G. Baier and X. X. Zhu, “GPU-based nonlocal filtering for large scale SAR processing,” in Proc. IEEE Int. Geosci. Remote Sens. Symp. (IGARSS), Jul. 2016, pp. 7608-7611.
[86] F. Sica, D. Cozzolino, X. X. Zhu , L. Verdoliva , and G. Poggi, “InSAR-BM3D: a nonlocal filter for SAR interferometric phase restoration,” IEEE Trans. Geosci. Remote Sens., vol. 56, no. 6, pp. 3456-3467, 2018.
[87] R. Chen, W. Yu, R. Wang, G. Liu, and Y. Shao, “Interferometric phase denoising by pyramid nonlocal means filter,” IEEE Geosci. Remote Sens. Lett., vol. 10, no. 4, pp. 826–830, Jul. 2013.
[88] Y. J. Wang, Y. K. Deng, R. Wang, and J. L. Wang,
“Adaptive multilooking based on complex patch for multitemporal interferometry,” IEEE J. Sel. Topics Appl. Earth Observ. Remote Sens., vol. 11, no. 3, pp. 907-918, 2018.
[89] F. Sica, D. Reale, G. Poggi, L. Verdoliva, and G. Fornaro, “Nonlocal adaptive multilooking in SAR multipass differential interferometry,” IEEE J. Sel. Topics Appl. Earth Observ. Remote Sens., vol. 8, no. 4, pp. 1727-1742, Apr. 2015.
[90] O. D’Hondt, C. López-Martínez, S. Guillaso, and O. Hellwich, “Nonlocal filtering applied to 3-D reconstruction of tomographic SAR data,” IEEE Trans. Geosci. Remote Sens., vol. 56, no. 1, pp. 272–285, Jan, 2018.
[91] X. X. Zhu, R. Bamler, M. Lachaise, F. Adam, Y. Shi, and M. Eineder, “Improving TanDEM-X DEMs by non-local InSAR filtering,” in Proc. EUSAR, Jun. 2014, pp. 1-4.
[92] Y. L. Shi, X. X. Zhu, R. Bamler, “Optimized parallelization of non-local means filter for image noise reduction of InSAR image,” In Proceedings of the 2015 IEEE International Conference on Information and Automation, 2015, Lijiang, China, 1515-1518.
[93] X. X. Zhu, et al, "Deep learning in remote sensing: A comprehensive review and list of resources," IEEE Geoscience and Remote Sensing Magazine, 5.4 (2017): 8-36.
[94] Li. P. Zhang, L. F. Zhang, and B. Du, "Deep learning for remote sensing data: A technical tutorial on the state of the art," IEEE Geoscience and Remote Sensing Magazine, 4.2 (2016): 22-40.
[95] L.C. Potter, E. Ertin, J.T. Parker, M. Çetin, "Sparsity and compressed sensing in radar imaging", Proc. IEEE, vol. 98, no. 6, pp. 1006-1020, 2010.
[96] M. Çetin, I. Stojanovic, N. O. Önhon, K. R. Varshney, S. Samadi, W. C. Karl, A. S. Willsky, "Sparsity-driven synthetic aperture radar imaging", IEEE Signal Process. Mag., vol. 31, no. 4, pp. 27-40, Jul. 2014.
[97] G. Ferraiulo and G. Poggi, “A Bayesian filtering technique for SAR interferometric phase fields,” IEEE Trans. Image Process., vol. 13, no. 10, pp. 1368–1378, Oct. 2004.
[98] L. Denis, F. Tupin, J. Darbon, and M. Sigelle, “Joint regularization of phase and amplitude of InSAR data: Application to 3-D reconstruction,” IEEE Trans. Geosci. Remote Sens., vol. 47, no. 11, pp. 3774-3785, Nov. 2009.
[99] G. Ferraiuolo, V. Pascazio, and G. Schirinzi, “Maximum a posteriori estimation of height profiles in InSAR imaging,” IEEE Geosci. Remote Sens. Lett., vol. 1, no. 2, pp. 66–70, Apr. 2004.
[100] G. Xu, M.-D. Xing, X.-G. Xia, L. Zhang, Y.-Y. Liu, and Z. Bao, “Sparse regularization of interferometric phase and amplitude for InSAR image formation based on Bayesian representation,” IEEE Trans. Geosci. Remote Sens., vol. 53, no. 4, pp. 2123-2136, Apr. 2015.
[101] Fang, Dongsheng, et al. “A sparsity-based InSAR phase denoising algorithm using nonlocal wavelet shrinkage,” Remote Sensing, vol.8, no.10, pp.830, 2016.
[102] H. Hongxing, J. M. Bioucas-Dias, and V. Katkovnik, “Interferometric phase image estimation via sparse coding in the complex domain,” IEEE Trans. Geosci. Remote Sens., vol. 53, no. 5, pp. 2587-2602, May 2015.
[103] X. M. Luo , X. F. Wang , Y. H. Wang, and S. Q. Zhu, “Efficient InSAR phase noise reduction via compressive sensing in the complex domain,” IEEE J. Sel. Topics Appl.
Earth Observ. Remote Sens., vol. 11, no. 5, pp. 1615-1632, 2018.
[104] C. Ojha, A. Fusco, I. M. Pinto, “Interferometric SAR phase denoising using proximity-based K-SVD technique,” Sensors, vol. 19, no.12, pp. 2684-2703, 2019.
[105] R. Chen, W. Yu, R. Wang, G. Liu, and Y. Shao, “Integrated denoising and unwrapping of InSAR phase based on Markov random fields,” IEEE Trans. Geosci. Remote Sens., vol. 51, no. 8, pp. 4473–4485, Aug. 2013.
[106] W. B. Abdallah and R. Abdelfattah, “A Joint Markov Random Field Approach for SAR Interferogram Filtering and Unwrapping,” IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing., vol. 9, no. 7, pp. 1–10, April. 2016.
[107] A. B. Suksmono, A. Hirose, “Adaptive noise reduction of InSAR images based on a complex-valued MRF model and its application to phase unwrapping,” IEEE Trans. Geosci. Remote Sensing, 2002, 40(3): 699-709.
[108] H. Y. Li, H. J. Song, R. Wang, H. Wang, G. Liu, R. P. Chen, X. L. Li, Y. K. Deng, T. Balz, “A modification to the complex-valued MRF modeling filter of interferometric SAR phase,” IEEE Geoscience and Remote Sensing Letters, 2014, 12(3): 681-685.
[109] J. Kang, Y. Y. Wang, M. Körner, and X. X. Zhu, “Robust object-based multipass InSAR deformation reconstruction,” IEEE Trans. Geosci. Remote Sens., vol. 55, no. 8, pp. 4239 - 4251, 2018.
[110] S. Mukherjee, et al, “CNN-based InSAR denoising and coherence metric,” IEEE sensors, 28-31, 2018.
[111] H. Yu, Y. Lan, Z. Yuan, J. Xu, and H. Lee, “Phase Unwrapping in InSAR: A Review,” IEEE Geosci. Remote Sens. Mag., vol. 7, no. 1, pp. 40-58, Mar, 2019.
[112] H. Yu, H. Lee, N. Cao, and Y. Lan, “Optimal baseline design for multi-baseline InSAR phase unwrapping,” IEEE Trans. Geosci. Remote Sens., vol. 57, no. 8, pp.5231–5243, Aug. 2019.
[113] H. Yu, Y. Lan, H. Lee, and N. Cao, “2-D phase unwrapping using minimum infinity-norm,” IEEE Geosci. Remote Sens. Lett., vol.15, no. 12, pp. 1887–1891, Dec. 2018.
[114] H. Yu, H. Lee, T. Yuan, and N. Cao, “A novel method for deformation estimation based on multi-baseline InSAR phase unwrapping,” IEEE Trans. Geosci. Remote Sens., vol. 56, no. 9, pp.5231–5243, Sept. 2018.
[115] M. Jiang, Z. Miao, P. Gamba, and B. Yong, “Application of multitemporal InSAR covariance and information fusion to robust road extraction,” IEEE Trans. Geosci. Remote Sens., vol. 55, no. 6, pp. 3611-3622, 2017.
[116] X. X. Zhu, G. Baier, M. Lachaise, Y. L. Shi, F. Adam, and R. Bamler, “Potential and limits of non-local means InSAR filtering for TanDEM-X high-resolution DEM generation,” Remote Sensing of Environment, vol. 218, pp. 148-161, 2018.