Feature-Based Nonlocal Polarimetric SAR Filtering

Polarimetric synthetic aperture radar (PolSAR) images are inherently contaminated by multiplicative speckle noise, which complicates the image interpretation and image analyses. To reduce the speckle effect, several adaptive speckle filters have been developed based on the weighted average of the similarity measures commonly depending on the model or probability distribution, which are often affected by the distribution parameters and modeling texture components. In this paper, a novel filtering method introduces the coefficient of variance (CV) and Pauli basis (PB) to measure the similarity, and the two features are combined with the framework of the nonlocal mean filtering. The CV is used to describe the complexity of various scenes and distinguish the scene heterogeneity; moreover, the Pauli basis is able to express the polarimetric information in PolSAR image processing. This proposed filtering combines the CV and Pauli basis to improve the estimation accuracy of the similarity weights. Then, the similarity of the features is deduced according to the test statistic. Subsequently, the filtering is proceeded by using the nonlocal weighted estimation. The performance of the proposed filter is tested with the simulated images and real PolSAR images, which are acquired by AIRSAR system and ESAR system. The qualitative and quantitative experiments indicate the validity of the proposed method by comparing with the widely-used despeckling methods.


Introduction
Synthetic aperture radar (SAR) explores the characterization of the Earth's surface in the microwave range, which makes it possible to observe and retrieve data under all weather conditions and all times.Polarimetric SAR (PolSAR) systems further enrich the remote sensing information source, as the appearance of polarization diversity makes the remote sensing system more sensitive to the scattering properties of the imaged scene.However, the coherent imaging system is inherently contaminated by multiplicative speckle noise because of the coherent nature of the scattering.Therefore, the speckle filtering is always a fundamental step for PolSAR image processing.
So far, numerous kinds of speckle filter methods have been developed.The goal of speckle filtering is to suppress speckle while preserve spatial resolution and polarimetric information.Thus, the adaptive filtering methods are developing rapidly by using weighted averages of homogeneous pixels with similar scattering property.These filter algorithms have been widely applied to reduce speckle by using a weighted combination of the homogeneous sample and its mean.The refined Lee filter selected the similar pixels to reinforce homogeneity according to eight non-square windows as the templates [1,2].Similarly, in the improved sigma filter, the two-sigma probability ranges were utilizes to select homogeneous pixels [3][4][5], and the intensity-driven adaptive-neighborhood (IDAN) was proposed by grouping pixels with similar statistical properties [6].Furthermore, the bilateral filter calculates the weighted average based on the similarity between pixels in the local windows of the spatial and radiometric domains [7][8][9].Generally, these filters are performed by the pixels close to the target, which implies these pixels that are similar in the local window.In reality, the notable differences between pixels exist in the heterogeneous regions, and the preservation of the image details is easily affected.
Recently, the nonlocal mean (NLM) filtering method applied to the PolSAR image computes the weights based on the patch similarity instead of the pixels [10].Due to the structure information contained in patches, the patch similarity has been adapted to most state-of-the-art denoising techniques based on the test statistic [11], Bayesian framework [12], etc.Moreover, the nonlocal Lee (NL-Lee) filter combines the Lee filter with the framework of the NLM filter, and the homogeneity similarity is introduced to enhance the patch regularity assumption [13].In [14], the NLM filter is involving in iterative steps and implemented in the form of maximum a posteriori (MAP).These nonlocal algorithms are mostly performed by applying the similarity measure of the matrix statistical estimation, and all the elements of the coherency matrix or covariance matrix are implemented in an identical manner.Therefore, the NLM filters depend on the model or probability distribution to measure the similarity, which are often affected by the distribution parameters and modeling texture components.Thus, few repeated structures may lead to poor results in the case of heterogeneous scenes with lots of details [15].
In the segmentation and classification methods, the extracted features can effectively measure the similarity between the pixels and the class center, and the different polarimetric signatures or features can distinguish the types of terrains from the different perspective [16][17][18][19].Inspired by this idea, the similarity between the central pixel and each neighborhood pixel could be evaluated by the extracted feature.
Given the importance of the structural features for the speckle filters, the coefficient of variance (CV) can reflect the complexity of different scenes, and it has been widely used in SAR and PolSAR image processing [20][21][22][23][24].Moreover, this difference of the textural features can be quantitatively described by CV, which could be used as the feature to measure the similarity.Electromagnetic scattering phenomena play a key role in the SAR image process.Moreover, the analysis of scattering mechanism for PolSAR image is the emphasis for information extraction and interpretation, etc. [25].Scattering-model-based speckle filters have been proposed by combining the classification [26,27].It is well known that the ground targets are much easier to be interpreted by the Pauli images.The Pauli basis (PB) corresponds to the polarization channel and can express most of the polarimetric information, and is usually used as the feature for polarimetric segmentation and classification [28,29].Thus, the Pauli basis is chosen to better describe the scattering mechanisms.
The nonlocal mean filtering method has its advantages, and the pixels could be weighted according to the patch similarity.Therefore, the proposed filter introduces the features to the framework of the NLM filter to improve the estimation accuracy of the similarity weights.In this paper, the selected features are the CV and the Pauli basis (PB) based on the above analysis.Additionally, the similarity weights of the features are deduced according to the test statistic, and the filter is implemented by using weighted average of similar pixels.The combination of the extracted features can effectively enhance the discrimination of the target scattering mechanisms, and the test of the pixel similarity is more robust.The filtering performance is exhibited by using both simulated and real PolSAR data.
The remainder of this paper is organized as follows.The proposed filtering method is explained in detail in Section 2. The experimental results and different filtering methods for comparison analysis are listed in Section 3.The discussion of the results is presented in Section 4. The conclusions are given in Section 5.

Methodology
The proposed filtering method consists of three main parts: (1)

Characteristics of PolSAR Data
The polarimetric characteristics of the ground objects are closely related to their geometric and dielectric properties.The characteristics can reflect the target surface roughness, symmetry, and other features information.PolSAR data offers efficient and reliable information required to extract biophysical and geophysical parameters about the Earth's surface.The polarimetric information of the target can be expressed by the Sinclair scattering matrix [28]: where the elements and produce the power return in the copolarized channels, and the elements and produce the power return in the cross-polarized channels.For a reciprocal target matrix, in the monostatic backscattering case, the Pauli basis vector for PolSAR can be expressed as: where the superscript is the transpose operator.The -look coherence matrix is given as T = 〈 * 〉 = ∑ ( ) ( ) * , where 〈•〉 is the ensemble average, * is the complex conjugate transpose, and is the number of looks.
In Pauli decomposition, the Sinclair scattering matrix is the complex sum of the Pauli matrices, with: where , , , and are all complex and are given by: In the monostatic backscattering case, where = 0 .It is noted that the Pauli basis is an incomplete representation of the polarimetric information of the scene due to only taking intensity into account.However, the expressions indicate that the Pauli spin matrix basis correspond to the matrix of the Pauli decomposition, and the polarimetric observations can be better applied.Thus, the Pauli basis is used to analyze the scattering mechanisms.

Scene Heterogeneity
The CV has been widely used to detect the scene heterogeneity, which was proposed to enhance the SAR filtering adaptively by Lopes et al. [20].It is defined as the ratio between the standard deviation and the mean in a local region as follows:

Characteristics of PolSAR Data
The polarimetric characteristics of the ground objects are closely related to their geometric and dielectric properties.The characteristics can reflect the target surface roughness, symmetry, and other features information.PolSAR data offers efficient and reliable information required to extract biophysical and geophysical parameters about the Earth's surface.The polarimetric information of the target can be expressed by the Sinclair scattering matrix [28]: where the elements S HH and S VV produce the power return in the copolarized channels, and the elements S HV and S V H produce the power return in the cross-polarized channels.
For a reciprocal target matrix, in the monostatic backscattering case, the Pauli basis vector k for PolSAR can be expressed as: where the superscript t is the transpose operator.The L-look coherence matrix is given as , where • is the ensemble average, * t is the complex conjugate transpose, and L is the number of looks.
In Pauli decomposition, the Sinclair scattering matrix is the complex sum of the Pauli matrices, with: where a, b, c, and d are all complex and are given by: In the monostatic backscattering case, where d = 0.It is noted that the Pauli basis is an incomplete representation of the polarimetric information of the scene due to only taking intensity into account.However, the expressions indicate that the Pauli spin matrix basis correspond to the matrix of the Pauli decomposition, and the polarimetric observations can be better applied.Thus, the Pauli basis is used to analyze the scattering mechanisms.

Scene Heterogeneity
The CV has been widely used to detect the scene heterogeneity, which was proposed to enhance the SAR filtering adaptively by Lopes et al. [20].It is defined as the ratio between the standard deviation σ and the mean µ in a local region as follows: For polarimetric SAR data, the scattering vector k usually follows a multivariate complex Gaussian distribution [28].Then, the coherency matrix T can be modeled as having the Wishart distribution: where , p is the dimension of complex scattering vector, L is the number of looks, Σ = E kk * t is Mathematical Expectation of matrix T, tr(•) is the trace of a matrix.However, the CV is based on the SAR image and deal with the intensity or amplitude image.The polarimetric whitening filtering (PWF) synthesizes a filtered real image by combining all elements of the scattering matrix [30,31].Thus, in this paper, the coherency matrix T can be mapped according to the PWF, given by: F has the statistical properties of the original data and is scale invariant, because of the normalized Σ in the local window.Therefore, the CV based on the PWF is used to measure the scene heterogeneity.The mean of F is scale invariant and equal to the dimension of the matrix, , where σ F denotes the standard deviation of F. Since p is a constant, the heterogeneity is determined by σ F .Since this transformation is performed, the CV of F could be described by the standard deviation σ F , i.e., 2.3.The NLM Filtering of CV and PB

Test for Similarity of the Features
For PolSAR image, according to the complex Wishart distribution, Conradsen et al. proposed the likelihood-ratio test for the equality of two complex Wishart matrices [32]: where L 1 , L 2 are the number of looks, p is the dimension of vector k.In the case of monostatic polarimetric SAR on a reciprocal medium, p = 3.For asingle-channel PolSAR image, p = 1, Q can be simplified as [32]: where X and Y are the scalar values for the single-channel HH, HV, or VV data.
Based on the test statistic, for PolSAR speckle filtering, the similarity of the coherence matrix between two patches X and Y is derived [11]: where M is the size of patch size.Thus, the similarity of a single-channel PolSAR image is: We observe that the first term and the equivalent number of looks L areconstantwhen the size of the smoothing window is fixed for a PolSAR image, and do not contribute to the test.Thus, D can be approximately computed as:

Procedure of Weighted Average
In 2005, Buades et al. [10] proposed a nonlocal mean filter algorithm.The idea of this algorithm is to replace the similarity of individual pixels with the structural similarity in the image, and to estimate the pixel value by using the redundant information.Thus, a large number of nonlocal filtering algorithms have been proposed and achieved good performance of speckle suppression, and the form is expressed as [10]: where û(x) and u(y) denote the pixel values in noise free and noisy image, respectively.Ω is the search window.The weight w(x, y) is computed using the noisy patches, which can be defined as: w(x, y) = exp(−D(x, y)/h); h is the smoothing parameter.Z(x, y) is the normalization constant, Z(x, y) = ∑ y∈Ω w(x, y), D(x, y) is the similarity measure between pixel x and y.
The speckle-reduced pixel of the PWF is obtained by the linear combination of intensities of |HH| 2 , |HV| 2 , |VV| 2 and a correlation term.Therefore, the CV based on the PWF is consistent with the statistical properties of a single-channel PolSAR image.Based on the above derivation, the similarity of two patches based on the CV F and the PB can be given as follows: Correspondingly, according to the NLM filter, the weights for the CV F and the PB are defined , the parameter h affects the weight decay.
Homogeneous samples should be close in both the CV F and the PB domains.Accordingly, the filter weight is the product of two feature weights and expressed as: It should be noted that the three values would be calculated according to the three components of the Pauli basis, and the test of the polarization feature is obtained by summing the three values.Moreover, the polarimetric covariance or coherence matrices are Hermitian positive definite matrices based on the test statistic in [32].Hence, when the rank of the matrices is less than three, the undefined channel values are discarded and the similarity measure is determined by other channel values.

Determination of Smoothing Parameters
The parameter h controls the degreeof smoothness in the proposed filter.In addition, the speckle reduction is much better with the parameters increasing.In the two similarity measures, we need to set the smoothing parameters of heterogeneity h CV F and polarimetric scattering h PB .Ideally, the texture information is desirable to be preserved in the higher heterogeneity region, and smaller smoothing parameters are required.In addition, the larger smoothing parameters are capable of suppressing the speckle in the homogeneous region.
In contrast, heterogeneity describes the complexity of the image.Therefore, the smoothing parameter of heterogeneity h CV F can be determined adaptively according to the reciprocal of the heterogeneity: h CV F = 1/CV F .To improve the parameter robustness, the mean value of the center patch CV F is used to determine the heterogeneity parameter.Moreover, the parameters relate to the patch size M and the number of looks L. The differences between patches increase when the patch size grows, but the effect of speckle is reduced with increasing the number of looks.Thus, the heterogeneity parameter can be defined as: The Pauli parameter is tuned experimentally and given as follows: The basic procedure of the proposed method is presented in Algorithm 1 as follows.Output: the filtered image.

Description of the Experimental Datasets
The simulated and real PolSAR datasets are utilized in the experiments to evaluate the performance of the proposed filter, including one-look and multilook datasets, respectively.The simulated PolSAR datasets are generated by using Monte Carlo simulation according to the procedure provided by Lee et al. [33].The given coherent matrices are obtained from the samples of the real PolSAR data, and the regions of different heterogeneity are generated by varying the shape factor in K-distribution [34].
The one-look simulated data with a size of 300 × 300 pixels is mainly composed of field and forest, which is shown in Figure 2a.Additionally, as shown in the red boxes, three isolated points were included to evaluate the preservation of strong targets of the filtered image.The ground truth image is generated by following the method in [35] with the given the coherence matrices and displayed in Figure 2b.The heterogeneity is measured by using the CV F .Given the effect of the patch size, the size is set to 3 × 3. Figure 2c shows the corresponding heterogeneity map.
The one-look simulated data with a size of 300×300 pixels is mainly composed of field and forest, which is shown in Figure 2a.Additionally, as shown in the red boxes, three isolated points were included to evaluate the preservation of strong targets of the filtered image.The ground truth image is generated by following the method in [35] with the given the coherence matrices and displayed in Figure 2b.The heterogeneity is measured by using the .Given the effect of the patch size, the size is set to 3×3.Figure 2c shows the corresponding heterogeneity map.The four-look data includes water body, road, field, forest, and building areas with the dimensions of 300×300 pixels, and the Pauli decomposition image is presented in Figure 3a.The The four-look data includes water body, road, field, forest, and building areas with the dimensions of 300 × 300 pixels, and the Pauli decomposition image is presented in Figure 3a.The corresponding ground truth image and heterogeneity map are shown in Figure 3b,c.As can be seen from the heterogeneity maps in Figures 2c and 3c, the features of the ground targets can be clearly detected, and the boundaries between classes do not spread.The CV F can effectively distinguish between homogeneous and heterogeneous areas and be used to measure the similarity.It should be noted that, except the building area, each ground target is relatively homogeneous with the lower heterogeneity, and the difference between the classes is clearer.
Remote Sens. 2017, 9, x FOR PEER REVIEW 7 of 18 corresponding ground truth image and heterogeneity map are shown in Figure 3b,c.As can be seen from the heterogeneity maps in Figures 2c and 3c, the features of the ground targets can be clearly detected, and the boundaries between classes do not spread.The can effectively distinguish between homogeneous and heterogeneous areas and be used to measure the similarity.It should be noted that, except the building area, each ground target is relatively homogeneous with the lower heterogeneity, and the difference between the classes is clearer.For illustration, two real PolSAR images are employed to validate the performance of the proposed approach: the first real PolSAR dataset is acquired by the AIRSAR system over San Francisco.The region of 400 × 400 pixels is selected from the four-look AIRSAR image, which is mainly composed of water, forest, point targets and urban areas.The Pauli decomposition image of the AIRSAR data is shown in Figure 4a.Correspondingly, the heterogeneity map calculated by is shown in Figure 4b.The second real PolSAR dataset uses a section of the one-look ESAR image over the Oberpfaffenhofen area, where both weak and strong point objects, line, and field areas are included.The Pauli decomposition image and the corresponding heterogeneity map are presented in Figure 5a,b, respectively.It can be seen that the object features such as the building area and the block boundary can be distinguished well.
Central Latitude: 37°46′46″N Central Longitude: 122°28′12″W For illustration, two real PolSAR images are employed to validate the performance of the proposed approach: the first real PolSAR dataset is acquired by the AIRSAR system over San Francisco.The region of 400 × 400 pixels is selected from the four-look AIRSAR image, which is mainly composed of water, forest, point targets and urban areas.The Pauli decomposition image of the AIRSAR data is shown in Figure 4a.Correspondingly, the heterogeneity map calculated by CV F is shown in Figure 4b.The second real PolSAR dataset uses a section of the one-look ESAR image over the Oberpfaffenhofen area, where both weak and strong point objects, line, and field areas are included.The Pauli decomposition image and the corresponding heterogeneity map are presented in Figure 5a,b, respectively.It can be seen that the object features such as the building area and the block boundary can be distinguished well.
the AIRSAR data is shown in Figure 4a.Correspondingly, the heterogeneity map calculated by is shown in Figure 4b.The second real PolSAR dataset uses a section of the one-look ESAR image over the Oberpfaffenhofen area, where both weak and strong point objects, line, and field areas are included.The Pauli decomposition image and the corresponding heterogeneity map are presented in Figure 5a,b, respectively.It can be seen that the object features such as the building area and the block boundary can be distinguished well.

Evaluation and Comparision
The filtering experiments are performed and compared by the local and nonlocal filtering methods, respectively.In this paper, the chosen filtering methods include (a) the refined Lee filter;(b) the iterative weight bilateral (ite-bilateral) filter according to the iterative filter weight of the spatial and polarimetric distances;(c) the nonlocal Lee (NL-Lee) filter based on the combination of the test statistics and the preselection;(d) the speckle reducing anisotropic diffusion (SRAD) filter by utilizing the partial differential equations [35][36][37][38];(e) the MuLog-TV filter of the latest algorithm applying Gaussian denoisers [39];and(f) the proposed CVPB-NLM filter method.
The size of the sliding window for the refined Lee filter is set to be 7 × 7. The ite-bilateral filter with fiveiterations uses an 11× 11 window, and the parameters of the weights sensitivity areℎ = 3, ℎ = 0.6.The nonlocal Lee (NL-Lee) filter uses a 15×15 search window and a 7×7 patch with a threshold coefficientof 20.In addition, a 3×3 neighborhood patch is set for the sigma filter to preselect more homogeneous region.In addition, the NL-Lee filter requires that the coherency or covariance matrix should be not singular, because the NLM weight employed in the NL-Lee filter is based on Hermitian positive definite matrices [13].Thus, in this paper, the NL-Lee filter is only applied to multilook PolSAR images.The SRAD filter is implemented with a maximum of 100 iterations and a time step ∆t = 0.1, including an 11 × 11 window for estimating the noise level.For the proposed CVPB-NLM filter method, the same search patch is set to be 15×15, with the 3×3 heterogeneity and 3×3 Pauli patches size, and the parameters are given as Equations ( 16) and ( 17).

Evaluation and Comparision
The filtering experiments are performed and compared by the local and nonlocal filtering methods, respectively.In this paper, the chosen filtering methods include (a) the refined Lee filter; (b) the iterative weight bilateral (ite-bilateral) filter according to the iterative filter weight of the spatial and polarimetric distances; (c) the nonlocal Lee (NL-Lee) filter based on the combination of the test statistics and the preselection; (d) the speckle reducing anisotropic diffusion (SRAD) filter by utilizing the partial differential equations [35][36][37][38]; (e) the MuLog-TV filter of the latest algorithm applying Gaussian denoisers [39]; and (f) the proposed CVPB-NLM filter method.
The size of the sliding window for the refined Lee filter is set to be 7 × 7. The ite-bilateral filter with fiveiterations uses an 11 × 11 window, and the parameters of the weights sensitivity are h s = 3, h p = 0.6.The nonlocal Lee (NL-Lee) filter uses a 15 × 15 search window and a 7 × 7 patch with a threshold coefficient of 20.In addition, a 3 × 3 neighborhood patch is set for the sigma filter to preselect more homogeneous region.In addition, the NL-Lee filter requires that the coherency or covariance matrix should be not singular, because the NLM weight employed in the NL-Lee filter is based on Hermitian positive definite matrices [13].Thus, in this paper, the NL-Lee filter is only applied to multilook PolSAR images.The SRAD filter is implemented with a maximum of 100 iterations and a time step ∆t = 0.1, including an 11 × 11 window for estimating the noise level.For the proposed CVPB-NLM filter method, the same search patch is set to be 15 × 15, with the 3 × 3 heterogeneity and 3 × 3 Pauli patches size, and the parameters are given as Equations ( 16) and (17).

Experiments on Simulated PolSAR Data
In this section, the simulated PolSAR data was employed to quantitatively analyze and evaluate the performance of the proposed approach.As shown in Figures 6 and 7, the different filtering methods can effectively reduce the speckle.Nevertheless, the edges and the texture details were largely smeared by the refined Lee filtering.Meanwhile, a noticeable pixelation effect was introduced.The ite-bilateral filter made a good performance on the speckle reduction, but some discrete points were introduced near the edges.Similarly, the SRAD tended to oversmooth the point targets, which was conducted according to the iterative nature and diffusion tensor.The edges of the proposed filtering and the MuLoG-TV filter seem more continuous with high contrast while speckle was effectively suppressed.As can be seen, the result of the NL-Lee filter was insufficiently smoothed.Additionally, the performance of the filter results can be quantitatively evaluated according to the reference ground truth image.As mentioned in [35,40], the preservation of the polarimetric and spatial information was analyzed, and the evaluation of the suppression was measured.Specifically, the amplitude and phase of the pixel were denoted by the complex correlation parameters (Ca, Cp).The radiometric parameter was corresponding to the power information (P).The decomposition parameters H/ᾱ/A, including the entropy (H), the anisotropy (A), and the mean alpha angle (ᾱ) were also estimated.These parameters were considered to evaluate the polarimetric information.Then, the absolute relative bias per scattering class was calculated for the filtered images and the ground truth image, and the median of different biases was used to evaluate the preservation of image information.Thus, the smaller value denoted the better ability to retain the information.
The edge preservation degree denoting the spatial information was evaluated according to the ratio of average (EPD-ROA) in the horizontal and vertical directions, respectively [42].The value close to one indicated good edge preservation.The equivalent number of looks (ENL) was applied to measure the degree of suppression, which was defined as the square ratio of the mean to the standard deviation values in a homogeneous region [40][41][42].The larger the ENL was, the better the quality of the speckle reduction was.In this experiment, the homogeneous areas of the simulated data were selected from the field and road.
The quantitative evaluation of the filter results was shown in Table 1, and the best values are marked in bold.For the one-look simulated data, it can be seen that the proposed CVPB-NLM Additionally, the performance of the filter results can be quantitatively evaluated according to the reference ground truth image.As mentioned in [35,40], the preservation of the polarimetric and spatial information was analyzed, and the evaluation of the suppression was measured.Specifically, the amplitude and phase of the pixel were denoted by the complex correlation parameters (Ca, Cp).The radiometric parameter was corresponding to the power information (P).The decomposition parameters H/ᾱ/A, including the entropy (H), the anisotropy (A), and the mean alpha angle (ᾱ) were also estimated.These parameters were considered to evaluate the polarimetric information.Then, the absolute relative bias per scattering class was calculated for the filtered images and the ground truth image, and the median of different biases was used to evaluate the preservation of image information.Thus, the smaller value denoted the better ability to retain the information.
The edge preservation degree denoting the spatial information was evaluated according to the ratio of average (EPD-ROA) in the horizontal and vertical directions, respectively [42].The value close to one indicated good edge preservation.The equivalent number of looks (ENL) was applied to measure the degree of suppression, which was defined as the square ratio of the mean to the standard deviation values in a homogeneous region [40][41][42].The larger the ENL was, the better the quality of the speckle reduction was.In this experiment, the homogeneous areas of the simulated data were selected from the field and road.
The quantitative evaluation of the filter results was shown in Table 1, and the best values are marked in bold.For the one-look simulated data, it can be seen that the proposed CVPB-NLM Additionally, the performance of the filter results can be quantitatively evaluated according to the reference ground truth image.As mentioned in [35,40], the preservation of the polarimetric and spatial information was analyzed, and the evaluation of the suppression was measured.Specifically, the amplitude and phase of the pixel were denoted by the complex correlation parameters (Ca, Cp).The radiometric parameter was corresponding to the power information (P).The decomposition parameters H/ e quantitatively evaluated according to the preservation of the polarimetric and suppression was measured.Specifically, e complex correlation parameters (Ca, wer information (P).The decomposition (A), and the mean alpha angle (ᾱ) were ate the polarimetric information.Then, for the filtered images and the ground to evaluate the preservation of image ty to retain the information.rmation was evaluated according to the directions, respectively [41].The value t number of looks (ENL) was applied to square ratio of the mean to the standard er the ENL was, the better the quality of neous areas of the simulated data were own in Table 1, and the best values are Additionally, the performance of the filter results can be quantitatively evaluated according to the reference ground truth image.As mentioned in [35,40], the preservation of the polarimetric and spatial information was analyzed, and the evaluation of the suppression was measured.Specifically, the amplitude and phase of the pixel were denoted by the complex correlation parameters (Ca, Cp).The radiometric parameter was corresponding to the power information (P).The decomposition parameters H/ᾱ/A, including the entropy (H), the anisotropy (A), and the mean alpha angle (ᾱ) were also estimated.These parameters were considered to evaluate the polarimetric information.Then, the absolute relative bias per scattering class was calculated for the filtered images and the ground truth image, and the median of different biases was used to evaluate the preservation of image information.Thus, the smaller value denoted the better ability to retain the information.
The edge preservation degree denoting the spatial information was evaluated according to the ratio of average (EPD-ROA) in the horizontal and vertical directions, respectively [41].The value close to one indicated good edge preservation.The equivalent number of looks (ENL) was applied to measure the degree of suppression, which was defined as the square ratio of the mean to the standard deviation values in a homogeneous region [40][41][42].The larger the ENL was, the better the quality of the speckle reduction was.In this experiment, the homogeneous areas of the simulated data were selected from the field and road.
The quantitative evaluation of the filter results was shown in Table 1, and the best values are ) were also estimated.These parameters were considered to evaluate the polarimetric information.Then, the absolute relative bias per scattering class was calculated for the filtered images and the ground truth image, and the median of different biases was used to evaluate the preservation of image information.Thus, the smaller value denoted the better ability to retain the information.
The edge preservation degree denoting the spatial information was evaluated according to the ratio of average (EPD-ROA) in the horizontal and vertical directions, respectively [41].The value close to one indicated good edge preservation.The equivalent number of looks (ENL) was applied to measure the degree of suppression, which was defined as the square ratio of the mean to the standard deviation values in a homogeneous region [40][41][42].The larger the ENL was, the better the quality of the speckle reduction was.In this experiment, the homogeneous areas of the simulated data were selected from the field and road.
The quantitative evaluation of the filter results was shown in Table 1, and the best values are marked in bold.For the one-look simulated data, it can be seen that the proposed CVPB-NLM obtained the best filtering result.The SRAD filter can better preserve the edges and the polarimetric information, while the speckle noise was insufficiently smoothed for the one-look simulated image.Meanwhile, the larger ENL and EPD-ROA indicated that the MuLoG-TV filter and the ite-bilateral filter were able to better suppress the speckle and preserve the edges for the one-look simulated image.For the four-look simulated data, because of the effect of discrete points near the edge, the ite-bilateral filter had the smallest ENL.The comprehensive evaluation of the one-look and four-look simulated data certifies that the proposed CVPB-NLM filter made a better tradeoff between speckle removal and information preservation.the edges and some tiny details were smoothed.The refined Lee filter preserved the strong points well; however, the weak ones were blurred.In comparison, the NL-Lee filtering obtained slightly improved visual quality and the edges were preserved better, but the speckle noise in the water was significantly higher than the ite-bilateral filtering and the proposed CVPB-NLM filtering.Visually, the performance of the MuLoG-TV filter was better than the NL-Lee filter, and the SRAD filter obtained a good preservation of spatial details.As can be seen, the CVPB-NLM filtering showed the better preservation of the edges and texture information in the forest and building areas, and the speckle noise was effectively suppressed.For the ESAR image, the corresponding filtered images were shown in Figure 9.We can see that the strong isolated points could be preserved well.Meanwhile, the refined Lee filter blurred the textures and introduced a pixelation effect in the homogeneous regions.The SRAD filtering gave a better detail preservation and the clearer edges.While some details of the ite-bilateral filtered image were removed.As shown, the proposed filter and the MuLoG-TV filter showed better performance of reducing speckle noise while preserving edges.
In the different filtering images, the uniform homogeneous and heterogeneous regions were chosen for computing the ENL and EPD-ROA, which were marked by red rectangles in Figures 4a  and 5a.Specifically, the marked Z1 region was located within the homogeneous area and applied to calculate the ENL, the other heterogeneous areas Z2 and Z3 were calculated for the EPD-ROA.In order to analyze and compare the effects on the polarimetric information of different filtering methods, three 40 × 40 homogeneous areas were selected according to [35].Then, the entire mean of the 40 × 40-pixel areas was assumed to be the reference data.The quantitative evaluation of different filter methods was conducted, which was similar to the evaluation of the simulated data, while the spatial correlation of the real datasets was much higher than the simulated image.The effect of the speckle on the original image was larger, and the entropy (H) of the decomposition parameters could be compared.the edges and some tiny details were smoothed.The refined Lee filter preserved the strong points well; however, the weak ones were blurred.In comparison, the NL-Lee filtering obtained slightly improved visual quality and the edges were preserved better, but the speckle noise in the water was significantly higher than the ite-bilateral filtering and the proposed CVPB-NLM filtering.Visually, the performance of the MuLoG-TV filter was better than the NL-Lee filter, and the SRAD filter obtained a good preservation of spatial details.As can be seen, the CVPB-NLM filtering showed the better preservation of the edges and texture information in the forest and building areas, and the speckle noise was effectively suppressed.For the ESAR image, the corresponding filtered images were shown in Figure 9.We can see that the strong isolated points could be preserved well.Meanwhile, the refined Lee filter blurred the textures and introduced a pixelation effect in the homogeneous regions.The SRAD filtering gave a better detail preservation and the clearer edges.While some details of the ite-bilateral filtered image were removed.As shown, the proposed filter and the MuLoG-TV filter showed better performance of reducing speckle noise while preserving edges.
In the different filtering images, the uniform homogeneous and heterogeneous regions were chosen for computing the ENL and EPD-ROA, which were marked by red rectangles in Figures 4a  and 5a.Specifically, the marked Z1 region was located within the homogeneous area and applied to calculate the ENL, the other heterogeneous areas Z2 and Z3 were calculated for the EPD-ROA.In order to analyze and compare the effects on the polarimetric information of different filtering methods, three 40 × 40 homogeneous areas were selected according to [35].Then, the entire mean of the 40 × 40-pixel areas was assumed to be the reference data.The quantitative evaluation of different filter methods was conducted, which was similar to the evaluation of the simulated data, while the spatial correlation of the real datasets was much higher than the simulated image.The effect of the speckle on the original image was larger, and the entropy (H) of the decomposition parameters could be compared.

Experiments on Real PolSAR Data
The different filtering results of the AIRSAR image were presented in Figure 8, respectively.As shown in Figure 8b, the ite-bilateral filtering smeared the point targets (as shown the red boxes), and the edges and some tiny details were smoothed.The refined Lee filter preserved the strong points well; however, the weak ones were blurred.In comparison, the NL-Lee filtering obtained slightly improved visual quality and the edges were preserved better, but the speckle noise in the water was significantly higher than the ite-bilateral filtering and the proposed CVPB-NLM filtering.Visually, the performance of the MuLoG-TV filter was better than the NL-Lee filter, and the SRAD filter obtained a good preservation of spatial details.As can be seen, the CVPB-NLM filtering showed the better preservation of the edges and texture information in the forest and building areas, and the speckle noise was effectively suppressed.
improved visual quality and the edges were preserved better, but the speckle noise in the water was significantly higher than the ite-bilateral filtering and the proposed CVPB-NLM filtering.Visually,the performance of the MuLoG-TV filter was better than the NL-Lee filter, and the SRAD filter obtained a good preservation of spatial details.As can be seen, the CVPB-NLM filtering showed the better preservation of the edges and texture information in the forest and building areas, and the speckle noise was effectively suppressed.For the ESAR image, the corresponding filtered images were shown in Figure 9.We can see that the strong isolated points could be preserved well.Meanwhile, the refined Lee filter blurred the textures and introduced a pixelation effect in the homogeneous regions.The SRAD filtering gave a better detail preservation and the clearer edges.While some details of the ite-bilateral filtered image were removed.As shown, the proposed filter and the MuLoG-TV filter showed better performance of reducing speckle noise while preserving edges.
In the different filtering images, the uniform homogeneous and heterogeneous regions were chosen for computing the ENL and EPD-ROA, which were marked by red rectangles in Figures 4a  and 5a.Specifically, the marked Z1 region was located within the homogeneous area and applied to calculate the ENL, the other heterogeneous areas Z2 and Z3 were calculated for the EPD-ROA.In order to analyze and compare the effects on the polarimetric information of different filtering methods, three 40 × 40 homogeneous areas were selected according to [35].Then, the entire mean of the 40 × 40-pixel areas was assumed to be the reference data.The quantitative evaluation of different filter methods was conducted, which was similar to the evaluation of the simulated data, while the spatial correlation of the real datasets was much higher than the simulated image.The effect of the speckle on the original image was larger, and the entropy (H) of the decomposition parameters could be compared.For the ESAR image, the corresponding filtered images were shown in Figure 9.We can see that the strong isolated points could be preserved well.Meanwhile, the refined Lee filter blurred the textures and introduced a pixelation effect in the homogeneous regions.The SRAD filtering gave a better detail preservation and the clearer edges.While some details of the ite-bilateral filtered image were removed.As shown, the proposed filter and the MuLoG-TV filter showed better performance of reducing speckle noise while preserving edges.
In the different filtering images, the uniform homogeneous and heterogeneous regions were chosen for computing the ENL and EPD-ROA, which were marked by red rectangles in Figures 4a  and 5a.Specifically, the marked Z1 region was located within the homogeneous area and applied to calculate the ENL, the other heterogeneous areas Z2 and Z3 were calculated for the EPD-ROA.In order to analyze and compare the effects on the polarimetric information of different filtering methods, three 40 × 40 homogeneous areas were selected according to [35].Then, the entire mean of the 40 × 40-pixel areas was assumed to be the reference data.The quantitative evaluation of different filter methods was conducted, which was similar to the evaluation of the simulated data, while the spatial correlation of the real datasets was much higher than the simulated image.The effect of the speckle on the original image was larger, and the entropy (H) of the decomposition parameters could be compared.Table 2 presents the quantitative evaluation with the quantitative indicators of the AIRSAR and ESAR images, moreover, the best results are marked in bold.As shown, ENL of different filters was lower for the AIRSAR image, because the selected relatively homogeneous region in water body was nonstationary in the experimental area.For the AIRSAR image, the proposed CVPB-NLM filter obtained the highest EPD-ROA results.In addition, it can be observed that the ite-bilateral filter and the proposed filter had a much larger ENL and indicated strong speckle reduction ability; however, the EPD-ROA of the ite-bilateral filter was significantly lower than the proposed filter.The ENL of the NL-Lee filter showed that it was not smooth enough.Meanwhile, the MuLoG-TV filter showed the best preservation of polarimetric information and a better speckle reduction ability.The quantitative indicators of the SRAD filter showed a good balance between noise removal and information preservation.Table 2 presents the quantitative evaluation with the quantitative indicators of the AIRSAR and ESAR images, moreover, the best results are marked in bold.As shown, ENL of different filters was lower for the AIRSAR image, because the selected relatively homogeneous region in water body was nonstationary in the experimental area.For the AIRSAR image, the proposed CVPB-NLM filter obtained the highest EPD-ROA results.In addition, it can be observed that the ite-bilateral filter and the proposed filter had a much larger ENL and indicated strong speckle reduction ability; however, the EPD-ROA of the ite-bilateral filter was significantly lower than the proposed filter.The ENL of the NL-Lee filter showed that it was not smooth enough.Meanwhile, the MuLoG-TV filter showed the best preservation of polarimetric information and a better speckle reduction ability.
The quantitative indicators of the SRAD filter showed a good balance between noise removal and information preservation.

Data
Method/Indicators the performance of the MuLoG-TV filter was better than the NL-Lee filter, and the SRAD filter obtained a good preservation of spatial details.As can be seen, the CVPB-NLM filtering showed the better preservation of the edges and texture information in the forest and building areas, and the speckle noise was effectively suppressed.For the ESAR image, the corresponding filtered images were shown in Figure 9.We can see that the strong isolated points could be preserved well.Meanwhile, the refined Lee filter blurred the textures and introduced a pixelation effect in the homogeneous regions.The SRAD filtering gave a better detail preservation and the clearer edges.While some details of the ite-bilateral filtered image were removed.As shown, the proposed filter and the MuLoG-TV filter showed better performance of reducing speckle noise while preserving edges.
In the different filtering images, the uniform homogeneous and heterogeneous regions were chosen for computing the ENL and EPD-ROA, which were marked by red rectangles in Figures 4a  and 5a.Specifically, the marked Z1 region was located within the homogeneous area and applied to calculate the ENL, the other heterogeneous areas Z2 and Z3 were calculated for the EPD-ROA.In order to analyze and compare the effects on the polarimetric information of different filtering methods, three 40 × 40 homogeneous areas were selected according to [35].Then, the entire mean of the 40 × 40-pixel areas was assumed to be the reference data.The quantitative evaluation of different filter methods was conducted, which was similar to the evaluation of the simulated data, while the spatial correlation of the real datasets was much higher than the simulated image.The effect of the speckle on the original image was larger, and the entropy (H) of the decomposition parameters could be compared.the edges and some tiny details were smoothed.The refined Lee filter preserved the strong points well; however, the weak ones were blurred.In comparison, the NL-Lee filtering obtained slightly improved visual quality and the edges were preserved better, but the speckle noise in the water was significantly higher than the ite-bilateral filtering and the proposed CVPB-NLM filtering.Visually, the performance of the MuLoG-TV filter was better than the NL-Lee filter, and the SRAD filter obtained a good preservation of spatial details.As can be seen, the CVPB-NLM filtering showed the better preservation of the edges and texture information in the forest and building areas, and the speckle noise was effectively suppressed.For the ESAR image, the corresponding filtered images were shown in Figure 9.We can see that the strong isolated points could be preserved well.Meanwhile, the refined Lee filter blurred the textures and introduced a pixelation effect in the homogeneous regions.The SRAD filtering gave a better detail preservation and the clearer edges.While some details of the ite-bilateral filtered image were removed.As shown, the proposed filter and the MuLoG-TV filter showed better performance of reducing speckle noise while preserving edges.
In the different filtering images, the uniform homogeneous and heterogeneous regions were chosen for computing the ENL and EPD-ROA, which were marked by red rectangles in Figures 4a  and 5a.Specifically, the marked Z1 region was located within the homogeneous area and applied to calculate the ENL, the other heterogeneous areas Z2 and Z3 were calculated for the EPD-ROA.In order to analyze and compare the effects on the polarimetric information of different filtering methods, three 40 × 40 homogeneous areas were selected according to [35].Then, the entire mean of the 40 × 40-pixel areas was assumed to be the reference data.The quantitative evaluation of different filter methods was conducted, which was similar to the evaluation of the simulated data, while the spatial correlation of the real datasets was much higher than the simulated image.The effect of the speckle on the original image was larger, and the entropy (H) of the decomposition parameters could be compared.The results of the ESAR image exhibited that the CVPB-NLM filter and the MuLoG-TV filter had the better speckle reduction ability.The SRAD filter showed the best preservation of polarimetric information and details information, while the ENL was the lowest.The MuLoG-TV filter showed the best preservation of polarimetric information and speckle reduction ability.Meanwhile, most of the quantitative indicators of the proposed method were better than the other filters.For the real PolSAR data, this proposed filter method and the MuLoG-TV filter achieved a better balance between the speckle suppression and the preservation of the polarimetric information.ormance of the filter results can be quantitatively evaluated according to image.As mentioned in [35,40], the preservation of the polarimetric and alyzed, and the evaluation of the suppression was measured.Specifically, of the pixel were denoted by the complex correlation parameters (Ca, eter was corresponding to the power information (P).The decomposition ng the entropy (H), the anisotropy (A), and the mean alpha angle (ᾱ) were meters were considered to evaluate the polarimetric information.Then, er scattering class was calculated for the filtered images and the ground ian of different biases was used to evaluate the preservation of image ller value denoted the better ability to retain the information.degree denoting the spatial information was evaluated according to the A) in the horizontal and vertical directions, respectively [41].The value edge preservation.The equivalent number of looks (ENL) was applied to ression, which was defined as the square ratio of the mean to the standard geneous region [40][41][42].The larger the ENL was, the better the quality of .In this experiment, the homogeneous areas of the simulated data were road.uation of the filter results was shown in Table 1, and the best values are one-look simulated data, it can be seen that the proposed CVPB-NLM result.The SRAD filter can better preserve the edges and the polarimetric ckle noise was insufficiently smoothed for the one-look simulated image.L and EPD-ROA indicated that the MuLoG-TV filter and the ite-bilateral uppress the speckle and preserve the edges for the one-look simulated mulated data, because of the effect of discrete points near the edge, the itellest ENL.The comprehensive evaluation of the one-look and four-look at the proposed CVPB-NLM filter made a better tradeoff between speckle reservation.
PolSAR Data results of the AIRSAR image were presented in Figure8, respectively.As -bilateral filtering smeared the point targets (as shown the red boxes), and polarimetric decomposition [28].The entropy can describe the polarimetric scattering degree of randomness.As it is shown, the entropy values of the forest areas were high, because of the random scattering process.The homogeneous regions had low degree of scattering randomness.In the building areas, different scattering mechanisms were mixed with the higher entropy values.The refined Lee filter and the SRAD filter preserved the textures in the forest zones, but the built-up regions of the refined Lee filter were blurred for the AIRSAR image.Meanwhile, the NL-Lee filter, the ite-bilateral filter and the MuLoG-TV filter blurred the textures in the forest zones for the AIRSAR image.The proposed method was the one that better preserved the low entropy small spots and lines within both the forest and urban area.The proposed method and the MuLoG-TV filter can distinguish the edges and line targets for the ESAR image, and the homogeneous and heterogeneous regions were clearly observed.The Canny detector [41] was applied to the HH band of the AIRSAR data, and the results were presented in Figure12.As shown, the Gaussian filter was used with the standard deviation of 0.5 in the Canny detector for each filtered image.It can be seen that a few edges and details were lost, even the Canny detector was applied to test the multi-look images.Specifically, the proposed CVPB-NLMfilter and the MuLoG-TV filter preserved the coastline better (red box 1).The refined Lee filter and the NL-Lee filter lose the significant information of forest area, while the detected edges of the SRAD filter were clearer (red box 3).The ite-bilateral filter failed to detect the fine details in the left part (red box 2), although a few linear features were visible in the middle of forest areas.In the building region, the detected edges of the NL-Lee filter and the proposed filter were closer, while the small details were preserved better by the proposed filter.Relatively, the MuLoG-TV filter detected the most details, but there are more false edges.The Canny detector [42] was applied to the HH band of the AIRSAR data, and the results were presented in Figure 12.As shown, the Gaussian filter was used with the standard deviation of 0.5 in the Canny detector for each filtered image.It can be seen that a few edges and details were lost, even the Canny detector was applied to test the multi-look images.Specifically, the proposed CVPB-NLMfilter and the MuLoG-TV filter preserved the coastline better (red box 1).The refined Lee filter and the NL-Lee filter lose the significant information of forest area, while the detected edges of the SRAD filter were clearer (red box 3).The ite-bilateral filter failed to detect the fine details in the left part (red box 2), although a few linear features were visible in the middle of forest areas.In the building region, the detected edges of the NL-Lee filter and the proposed filter were closer, while the small details were preserved better by the proposed filter.Relatively, the MuLoG-TV filter detected the most details, but there are more false edges.The Canny detector [41] was applied to the HH band of the AIRSAR data, and the results were presented in Figure12.As shown, the Gaussian filter was used with the standard deviation of 0.5 in the Canny detector for each filtered image.It can be seen that a few edges and details were lost, even the Canny detector was applied to test the multi-look images.Specifically, the proposed CVPB-NLMfilter and the MuLoG-TV filter preserved the coastline better (red box 1).The refined Lee filter and the NL-Lee filter lose the significant information of forest area, while the detected edges of the SRAD filter were clearer (red box 3).The ite-bilateral filter failed to detect the fine details in the left part (red box 2), although a few linear features were visible in the middle of forest areas.In the building region, the detected edges of the NL-Lee filter and the proposed filter were closer, while the small details were preserved better by the proposed filter.Relatively, the MuLoG-TV filter detected the most details, but there are more false edges.

Main Features of the Proposed Method
In this paper, the features of the heterogeneity and polarimetric scattering are introduced to the proposed filtering method for PolSAR images.The combination of the extracted features effectively improves the estimation accuracy of the similarity weights.In addition, the influence of the distribution model and the probability density parameter are reduced.The filter makes a good balance between speckle removal and detail preservation.
The first feature of this paper is that the pixel similarity is measured based on scene heterogeneity.The heterogeneity is able to quantitatively and objectively measure the complexity of scenes.Moreover, the heterogeneity is an important parameter reflecting image texture characteristics, which has been widely used in SAR(PolSAR) image segmentation, classification, extraction, and filtering processing [29,43,44].The CV is the most commonly used heterogeneity test method.In this paper, the CV is calculated according to the PWF.The accuracy of the heterogeneity can be effectively improved by the CV based on the PWF (CV F ). Thus, the CV F is used as the feature to measure the heterogeneity similarity.The Pauli basis is used as the second feature to test the scattering mechanism.In PolSAR images, the components of the Pauli basis can represent the signatures of the objects and have physical meanings to interpret the scattering mechanisms.The Pauli RGB image is commonly applied to display the PolSAR images.Based on the above analysis, the application of the features in the novel filter reduces the dependence on the distribution model and the probability density.
The relevant filtering methods [7][8][9][11][12][13][14] perform the weighted average by measuring the similarity of the statistical target information or spatial distance.The spatial distance is unsuitable for multiplicative speckle noise of the PolSAR image.In addition, the filter results may be affected with a less repeated structure.In this paper, the similarity of the features is deduced according to the test statistic.The proposed filter is performed by using the weighted average of similar pixels.Compared with other filtering methods, this approach integrates the heterogeneity and polarimetric properties to evaluate the pixels similarity and enhance speckle reduction capabilities.Due to the use of target features, the stability of the proposed algorithm is improved, and it is fit for single-look and multi-look PolSAR images.

Sensitivity Analysis of the Smoothing Parameters
In the proposed method, the two smoothing parameters are defined, including heterogeneity h CV F and the Pauli parameter h PB .According to the reciprocal of the heterogeneity, the parameter h CV F can be determined adaptively; meanwhile, the patch size M and the number of looks L are included.Thus, for a PolSAR image, the map of the parameter h CV F is determined when the patch size is fixed.Similarly, the Pauli parameter h PB is a function of the patch size M and the number of looks L. The Pauli smoothing parameter can be experimentally tuned based on the heterogeneity smoothing parameter.In order to compare with the ground truth, the experimental results based on the simulated PolSAR images are analyzed according to the different Pauli parameter.
Figure 13a presents the maps of the parameter h CV F , and the experimental results of the different Pauli smoothing parameters are shown in Figure 13b-d.As can be seen, the values of the edges, point targets and the building areas are significantly smaller than the homogeneous regions, and the parameter h CV F is able to adaptively control the degree of smoothness according to the map.The filtered results based on the Pauli parameters are affected by applying the different Pauli parameters, and in this paper the determined parameter could achieve better performance.As can be seen, the result shows reasonably good characteristics of preserving edges and scattering properties, andsmall features are much clearer with a larger patch.While, when the smoothing parameters are determined, speckle filtering power is reduced with the increase of the patch.Since, in a larger patch, the difference is greater than in a smaller patch, the pixel targets are regarded as the heterogeneity

Size of the Patch
The influence of the patch sizes is analyzed for the final results, and the experiments are conducted based on the different patch sizes ranging from 3 × 3 to 11 × 11.For comparison, the experimental results of different patch sizes are shown in Figures 14-16.As can be seen, the result shows reasonably good characteristics of preserving edges and scattering properties, andsmall features are much clearer with a larger patch.While, when the smoothing parameters are determined, speckle filtering power is reduced with the increase of the patch.Since, in a larger patch, the difference is greater than in a smaller patch, the pixel targets are regarded as the heterogeneity patches.

Size of the Patch
The influence of the patch sizes is analyzed for the final results, and the experiments are conducted based on the different patch sizes ranging from 3×3 to 11×11.For comparison, the experimental results of different patch sizes are shown in Figures 14,15 and 16.As can be seen, the result shows reasonably good characteristics of preserving edges and scattering properties, andsmall features are much clearer with a larger patch.While, when the smoothing parameters are determined, speckle filtering power is reduced with the increase of the patch.Since, in a larger patch, the difference is greater than in a smaller patch, the pixel targets are regarded as the heterogeneity patches.The quantitative results of different patch sizes are shown in Tables 3 and 4, and the best results are marked in bold.As shown in Table 3, the complex correlation parameters of the simulated data are greatly influenced with the different patch sizes.In addition, there is little difference of the polarimetric information and edge preservation for the simulated image.According to Tables 3 and  4, the ENL of the filtering method decreases by enlarging the size of the patch.The speckle suppression ability is weakened, while the EPD is slightly greater.In conclusion, the choice of the different patch sizes can affect the speckle reduction of the proposed filter, whereas the preservation The quantitative results of different patch sizes are shown in Tables 3 and 4, and the best results are marked in bold.As shown in Table 3, the complex correlation parameters of the simulated data are greatly influenced with the different patch sizes.In addition, there is little difference of the polarimetric information and edge preservation for the simulated image.According to Tables 3 and  4, the ENL of the filtering method decreases by enlarging the size of the patch.The speckle suppression ability is weakened, while the EPD is slightly greater.In conclusion, the choice of the The quantitative results of different patch sizes are shown in Tables 3 and 4, and the best results are marked in bold.As shown in Table 3, the complex correlation parameters of the simulated data are greatly influenced with the different patch sizes.In addition, there is little difference of the polarimetric information and edge preservation for the simulated image.According to Tables 3 and 4, the ENL of the filtering method decreases by enlarging the size of the patch.The speckle suppression ability is weakened, while the EPD is slightly greater.In conclusion, the choice of the different patch sizes can affect the speckle reduction of the proposed filter, whereas the preservation of spatial resolution remains relatively stable.Smaller patches could achieve better balance between noise removal and detail preservation.the edges and some tiny details were smoothed.The refined Lee filter preserved the strong points well; however, the weak ones were blurred.In comparison, the NL-Lee filtering obtained slightly improved visual quality and the edges were preserved better, but the speckle noise in the water was significantly higher than the ite-bilateral filtering and the proposed CVPB-NLM filtering.Visually, the performance of the MuLoG-TV filter was better than the NL-Lee filter, and the SRAD filter obtained a good preservation of spatial details.As can be seen, the CVPB-NLM filtering showed the better preservation of the edges and texture information in the forest and building areas, and the speckle noise was effectively suppressed.
Table 1.Quantitative evaluation of the simulated PolSAR data.
the obtained entropy images, and the Canny detector results were able to more effectively distinguish the various objects than other methods.Future work will integrate more features to enhance the estimation accuracy of the similarity measure.In addition, the metric of the distance between the pixels is important for the proposed filtering.Further research is needed to improve the detection of the weak objects and measure the distance more effectively.
the representation of the PolSAR data characteristics and Pauli decomposition; (2) the heterogeneity CV based on the polarimetric whitening filtering (PWF); (3) the nonlocal mean (NLM) filtering based on the CV and the Pauli basis.

Figure 1 18 Figure 1 .
Figure 1 gives the flowchart of the proposed filtering algorithm.The details of each process are further described below.Remote Sens. 2017, 9, x FOR PEER REVIEW 3 of 18

Figure 1 .
Figure 1.Overview of the proposed algorithm for PolSAR data.

Algorithm 1 :
Nonlocal bilateral filtering based on heterogeneity and scattering.INPUT: PolSAR image coherence matrix (or covariance matrix).(1) For each pixel, the weight is calculated according to the test statistic between two patches.(a) The CV F is calculated from the PWF filtered image, and a heterogeneity map is generated.(b) Calculate the parameter h CV F and h PB according to the CV F , the number of looks L and the patch size M. (c) The test statistic is measured between two patches of the scalar values including the CV F and the Pauli basis of the original image, respectively.(d) The weights are performed based on the test statistic.(2) The filtered coherence matrix (or covariance matrix) is performed by the weighted average based on the combination of the weights for the CV F and the Pauli basis.

Figure 2 .
Figure 2.One-look simulated image used in the experiment: (a) one-look simulated image; (b) onelook ground truth;and (c)one-look heterogeneity map.

Figure 2 .
Figure 2. One-look simulated image used in the experiment: (a) one-look simulated image; (b) one-look ground truth; and (c) one-look heterogeneity map.

Figure 3 .
Figure 3.Four-look simulated image used in the experiment: (a)four-look simulated image; (b) fourlook ground truth; and(c)four-look heterogeneity map.

Figure 3 .
Figure 3. Four-look simulated image used in the experiment: (a) four-look simulated image; (b) four-look ground truth; and(c) four-look heterogeneity map.

Figure4.
Figure4.L-band,multi-look PolSAR image of AIRSAR as the first real dataset used in the experiment: (a)originaldata of AIRSAR; and (b) the heterogeneity image of (a).

Figure 4 . 18 CentralFigure5.
Figure 4. L-band, multi-look PolSAR image of AIRSAR as the first real dataset used in the experiment: (a) originaldata of AIRSAR; and (b) the heterogeneity image of (a).Remote Sens. 2017, 9, x FOR PEER REVIEW 8 of 18

Figure 5 .
Figure 5. Single-look PolSAR image of ESAR as the second real dataset used in the experiment: (a) original data of ESAR; and (b) the heterogeneity image of (a).

Figures 10 and 11
presented the entropy images of the different filtering methods on AIRSAR and ESAR data, which were obtained from H/A/ the one-look simulated data using different filtering algorithms.ite-bilateral filter; (c) SRAD filter; (d) MuLoG-TV filter; and (e) CVPB-NLM (c) (d) (e) (f) ction images of the four-look simulated data using different filtering ee filter; (b) ite-bilateral filter; (c) NL-Lee filter; (d) SRAD filter; (e) MuLoG-NLM filter.

Figure13.
Figure13.The map of the parameter ℎ and the filtered results of simulated PolSAR images for different Pauli parameter based on a 3×3 patch size.(a)The map of the parameter ℎ ; (b) Pauli parameter 0.5 ℎ ;(c) Pauli parameter ℎ size; and (d) Pauli parameter 2 ℎ .4.2.2.Size of the Patch The influence of the patch sizes is analyzed for the final results, and the experiments are conducted based on the different patch sizes ranging from 3×3 to 11×11.For comparison, the experimental results of different patch sizes are shown in Figures 14, 15 and 16.As can be seen, the result shows reasonably good characteristics of preserving edges and scattering properties, andsmall features are much clearer with a larger patch.While, when the smoothing parameters are determined, speckle filtering power is reduced with the increase of the patch.Since, in a larger patch, the difference is greater than in a smaller patch, the pixel targets are regarded as the heterogeneity

Figure 13 .
Figure 13.The map of the parameter h CV F and the filtered results of simulated PolSAR images for different Pauli parameter based on a 3 × 3 patch size.(a) The map of the parameter h CV F ; (b) Pauli parameter 0.5 h PB ; (c) Pauli parameter h PB size; and (d) Pauli parameter 2 h PB .

Figure13.
Figure13.The map of the parameter ℎ and the filtered results of simulated PolSAR images for different Pauli parameter based on a 3×3 patch size.(a)The map of the parameter ℎ ; (b) Pauli parameter 0.5 ℎ ;(c) Pauli parameter ℎ size; and (d) Pauli parameter 2 ℎ .

Figure14.
Figure14.The filtered results of simulated PolSAR images for different patch sizes.(a)3×3 of the patch size; (b) 5×5 of the patch size;(c) 7×7 of the patch size; (d) 9×9 of the patch size; and (e) 11×11of the patch size.

Figure 16 .
Figure 16.The filtered results of ESAR image for different patch sizes.(a) 3 × 3 of the patch size; (b) 5 × 5 of the patch size; (c) 7 × 7 of the patch size;(d) 9 × 9 of the patch size; and (e) 11 × 11 of the patch size.

Table 1 .
Quantitative evaluation of the simulated PolSAR data.

Table 1 .
Quantitative evaluation of the simulated PolSAR data.

Table 1 .
Quantitative evaluation of the simulated PolSAR data.

Table 2 .
Quantitative Evaluation of the real PolSAR data.

Table 1 .
Quantitative evaluation of the simulated PolSAR data.

Table 1 .
Quantitative evaluation of the simulated PolSAR data.

Table 3 .
Filter results of different patch sizes for simulated PolSAR images.