I Introduction
Speckle reduction of polarimetric synthetic aperture radar (PolSAR) data is often necessary to improve the outcome of applications such as physical parameter estimation, terrain classification and object detection, among others. Polarimetric SAR systems provide information related to the backscattering mechanisms resulting of the interaction between a polarized electromagnetic wave and physical media. However, in the case of many natural media, due to the random distribution of the elementary scatterers present inside a resolution cell, phase and intensity cannot be predicted in a deterministic way
[1]. The effect of the complex interactions between wave and scatterers is a pixeltopixel variability called speckle, that compromises the interpretation of the data. Though resulting of a physical mechanism, speckle has to be described as a signal dependent noise and to be handled in a statistical framework.PolSAR systems provide images of complex vectors called target vectors. Due to the random properties of distributed targets, the coherent information carried by these vectors cannot be used directly and the data has to be processed in an incoherent way. Consequently, the covariance matrix of target vectors has to be estimated by computing an average over a set of independent samples. It is then possible to perform incoherent polarimetric decomposition methods
[2] that allow a physical interpretation in terms of scattering mechanisms. Increasing the number of samples involved in the average allows to reduce speckle provided the statistical properties of the data are spatially stationary. Otherwise, pixels with different underlying scattering mechanisms are mixed, resulting in a loss of spatial resolution. Moreover, image statistics can be described by the fullydeveloped speckle model [1] only for distributed targets. If a strong return from a single scatterer dominates the overall response, the signal is better approximated by a deterministic target. In this case, averaging may be avoided and coherent decompositions can be directly applied (see for example [3], p. 214) to retrieve the polarimetric properties.The challenge faced by speckle reduction methods is then to perform an adaptation to the local structure of the image so as to reduce speckle variability without altering the polarimetric and spatial properties of the image. Thus, homogeneous areas have to be smoothed while boundaries between such regions have to be preserved. As for deterministic targets, they have to be left unfiltered.
Several adaptive methods have already been proposed for speckle reduction of PolSAR data. The popular refined Lee filter [4] is based on a locally linear minimum mean square error assumption and performs matrix averaging so as to preserve the polarimetric properties of the signal. For a better preservation of spatial edges, a set of oriented windows is used. The refined Lee filter is computationally inexpensive and leads to good results as long as the spatial structures in the image are locally aligned with one of the oriented windows.
A region growing based method called intensitydriven adaptiveneighborhood (IDAN) allows an adaptation to the local structure of the image by extending Lee’s minimum mean square approach to spatially adaptive neighborhoods [5]. However, only intensity information is considered in the selection of pixels in the neighborhood.
A model based filter named MBPolSAR [6] allows the independent filtering of the matrix elements by considering an additivemultiplicative multidimensional speckle model. Although taking the full polarimetric information into account, this last method does not take advantage of the spatial information present in edges.
A variational method based on anisotropic diffusion was also proposed [7]
, achieving a high equivalent number of looks while preserving spatial information. Due to the use of partial differential equations, such a method has a high computational cost and is prone to numerical instabilities.
A recent work exploits a multiresolution representation of the image called binary partition tree (BPT) to group pixels with similar polarimetric properties [8]. Though achieving a high number of looks and ensuring the preservation of strong edges, this method only allows the estimation of the mean covariance on a region basis, making it similar to a segmentation procedure.
Other approaches extended the nonlocal means filter to PolSAR data [9, 10]. Nonlocal means computes a weighted average of pixels according to the similarities between local neighborhoods called patches. Although the patches may in theory be located in any part of the image, search is generally restricted to a window to reduce the computational cost. One drawback of this method is that other similar neighborhoods have to be present in the search window.
In this work, we propose to take advantage of the full polarimetric information contained in covariance matrices while preserving the spatial resolution by adapting the bilateral filter [11] to the processing of polarimetric covariance matrices. This filter has been widely used in image processing due to its simplicity and its remarkable performance in terms of edge preservation. The idea of the filter is to compute a weighted mean of pixels contained in a local window according to their spatial and radiometric proximity. The weights are determined by two Gaussian kernels operating in radiometric and spatial domains. Furthermore, it has been shown that the iterative version of the bilateral filter is closely related to anisotropic diffusion [12]. To deal with PolSAR data the bilateral weights have to be adapted to the Hermitian positive definite nature of covariance matrices. One obvious choice would be to replace the traditional Euclidean distance between vectors by the Frobenius norm of matrix differences. However, this would lead to poor results due to the presence of speckle. It is then more suitable to consider distances that are adapted to the statistical nature of speckle. Statistical distances such as the KullbackLeibler distance [13] may then be considered. Moreover, recent studies [14, 15] have shown that a Riemannian affineinvariant metric was a natural choice due to the manifold structure of such matrices. Therefore, we propose to investigate the use of such distances to compute the radiometric proximity between two covariance matrices. We also consider a simpler distance called logEuclidean [16]
that has been successfully applied in the field of diffusion tensor imaging (DTI).
Section II introduces the basic concepts leading to the statistical speckle model for PolSAR data. Section III describes the bilateral filter for graylevel and color images. In section IV we introduce the bilateral filtering of PolSAR data based on distances between covariance matrices. The expression of a matrix based filter is given and several choices of distances are discussed. We also introduce an iterative scheme to progressively refine the estimate. Section V deals with the validation of the filter on synthetic and experimental data. A quantitative performance evaluation allows to understand the effects of parameters in terms of edge preservation, quality of polarimetric reconstruction and amount of smoothing. The new method is then compared to other existing methods using the criteria mentioned above. The method is then validated over two experimental datasets. Finally, the paper is concluded in section VI.
Ii Polarimetric SAR data
Polarimetric SAR systems measure the relation between the transmitted and received electromagnetic wave in two orthogonal polarizations in the form of a scattering matrix [3]
(1) 
where and denote horizontal and vertical polarization, respectively. The reciprocity assumption in the monostatic case leads to . For further analysis, it is convenient to represent the scattering information in the form of a target vector containing the elements of
(2) 
(3) 
where and denote the lexicographic and Pauli bases.
On distributed targets, due to the complex interactions between the wave and multiple scatterers inside the resolution cell, this information is generally considered in a statistical framework. Under the fully developed speckle hypothesis, the target vectors follow a complex circular
variate Normal distribution
[17](4) 
where is the determinant of a matrix and represents the conjugate transpose of a complex vector. This distribution is fully described by its covariance matrix that contains information about power and relative phase between polarimetric channels. If then is called polarimetric covariance, and if then is called polarimetric coherency.
Unfortunately, has to be estimated from the data over several independent samples, resulting in a loss of spatial resolution. The sample covariance is expressed as
(5) 
and is obtained by an operation called multilooking that consists in taking the average of contiguous pixels.
Under this Gaussian hypothesis, the multilook covariance matrix follows a complex Wishart density [17]. The estimation accuracy can be improved by increasing the number of samples. Unfortunately, increasing this number is only possible on homogeneous areas where pixels arise from identical scattering mechanisms.
Iii The bilateral filter for graylevel and color images
Let us consider a 2D gray level image where a pixel is referenced by its 2D location . The original expression of the bilateral filter at is [11]
(6) 
with
(7) 
where the operator is the norm. Weighting functions and have to be chosen such as they attain a maximum in zero and tend to zero if their argument goes to infinity. In this work we adopt the functions defined in the original version of the filter [11] where and with
(8) 
The parameter controls the spatial extent of the filter and is comparable to the window size used in other techniques like the Lee filter [4]. The parameter controls the amount of filtering according to the radiometric proximity between two pixels. The filter thus performs a local average of intensities, assigning weights to pixels depending on their spatial and radiometric similarity to the pixel located at .
The filter is easily generalized to vectorial (e.g. color) images by replacing the absolute value by the norm in equation (6) and (7). In practice, the filter is applied to discrete images and approximated by a weighted sum of the pixel values in a finite local window centered around
(9) 
where the weights are defined by
(10) 
Fig. 1 shows a visual interpretation of the filter behaviour for two examples of spatial structures in grey level images. As it may be observed, the strength of the bilateral filter is that it can adapt not only to straight edges but to any shape of discontinuity provided the edge is characterized by an abrupt intensity transition. For PolSAR data, this should constitute an advantage over other traditional filtering methods that use a discrete set of oriented windows. A drawback of this method is that it relies on the value of the central pixel that is uncertain because of noise. However, this limitation may be overcome by considering an iterative version of the filter as we show in the following. For an extensive survey on applications and variants of the bilateral filter, that is beyond the scope of this paper, the reader may refer to [18].
Iv Extension of the filter to PolSAR data
Iva Bilateral filtering of covariance matrices
To be able to use the bilateral filter on PolSAR data, it has to handle polarimetric covariance matrices instead of vectors. Equation (9) thus becomes a weighted sum of covariances
(11) 
Then the weights given by equation (10) have to be expressed for matrices. Instead of norms of vector differences, matrix distances have to be used so that the weights act like similarities that attain a maximum for identical pixels and tend to zero if the distance goes to infinity. The spatial term of the weights remains unchanged.
The new expression for the weights by replacing the Euclidean norm by a matrix distance is
(12) 
In the following we introduce several distances that are suitable for such a filter.
IvB The symmetrized KullbackLeibler divergence
The KullbackLeibler divergence [13]
has been defined to characterize the discrepancy between two probability distributions
and(13) 
Since this measure is not symmetric, a symmetrized KullbackLeibler divergence may be defined as
(14) 
and verifies . In the case of variate complex circular Normal distributions and , it can be shown that this pseudodistance depends only on the covariance matrices and is expressed as
(15) 
where is the dimension of the multivariate density. Then, provided estimates of the covariances are available, this pseudodistance can be used to compute the weights in equation (12).
IvC Distances on the cone of Hermitian positive definite matrices
Since the covariance matrices lie on the symmetric cone of Hermitian positive definite (HPD) matrices, which is not a vector space but a Riemannian manifold, recent studies showed it was natural to use appropriate distances [14, 15]. An affine invariant metric has been proposed as a replacement for the Euclidean metric to deal with HPD matrices. The corresponding distance between two matrices and is [15]
(16) 
where is the matrix logarithm and the Frobenius norm. This distance corresponds to the length of the geodesic (i.e. the shortest path) between two points on the manifold and is invariant by affine transformation.
IvD Iterative scheme and implementation details
Iterative filtering
The distances introduced in the previous sections assume that the covariance matrices are known for each pixel of the PolSAR image. In practice, this is not the case since the original data is available in the form of the scattering matrix described in equation (1). Moreover, those distances need to be computed on a fullrank matrix, that requires multilooking (defined in equation (5)) of at least independent pixels. Therefore, this multilook image, though leading to a noisy covariance estimate, will be used as an initialization for filtering the image in an iterative way, so as to progressively refine the estimate [12, 19]
(18) 
where the weights (12) are computed on the covariances from iteration .
Weight of the reference pixel
An issue of this type of filter when applied to noisy data is the very high importance given to the central pixel relative to its neighbours, leading to many unfiltered pixels. In fact, since this pixel is taken as a reference, it is always assigned a weight of one (corresponding to a distance of zero). However, its value is uncertain because it is also affected by noise. Consequently, there is no reason to give this pixel such a high weight in the average. To avoid this excessive weighting, authors of [20] propose to set the weight of the central pixel to the maximum weight among its neighbors. We use this recommendation in our experiments.
Rank 1 targets
The use of the statistical and geodesic distances involves matrix inversion and logarithm. In order to avoid numerical problems, the involved matrices have to be fullrank. In theory, in the case of Gaussian distributed targets, a minimal multilooking ensures that those matrices have this property. However, many pixels in PolSAR images cannot be described by distributed targets. If the resolution cell contains a dominant scatterer (due to, for instance, a corner reflector), the target behaviour is deterministic with one dominant scattering mechanism. The ideal model for this type of pixel is a rank 1 matrix, that cannot in theory be handled by the BLF. In practice, due to multilooking and the scattering complexity of natural and made man targets, purely deterministic targets are unlikely to be present in the data. However, it is not excluded that some matrices are close to pure targets and produce numerical instability in the filtering process. For this reason, we propose a modification of the filter that allows to avoid the filtering of illconditioned matrices. For each pixel, before applying the filter, we compute the condition number of the matrix:
(19) 
where and
are the maximal and minimal eigenvalues of the matrix. If the matrix is illconditioned, the condition number goes to infinity and the inverse of this value is close to zero
[21]. Then we compare to a threshold (set to in our experiments) to detect rank deficient matrices. If the pixel is the central pixel in the window, it is left unfiltered. If the pixel is not the central one, its weight is simply set to zero so that it is excluded from the local average.Box  Ref. Lee  IDAN  BLF  BLF  BLF  
6.83  3.43  2.63  1.50  1.15  1.14  
54.5  17.3  6.90  1.71  1.35  1.37  
206  95.3  86.3  492  683  696 




V Experiments
To validate the method, we have divided our experiments into several main steps: we first evaluate the influence of the filter parameters on simulated data generated from sample covariances measured on an experimental dataset. This evaluation helps us to tune the filter parameters. Then, we compare our approach to other filters from the literature. Finally, we apply the method to experimentally acquired data sets.
Va Validation on simulated data
Let us first introduce the strategy employed to generate our synthetic data. To obtain realistic simulations, we have selected four homogeneous areas from the fully polarimetric Lband Oberpfaffenhofen image acquired by the ESAR sensor from DLR in 1998. The mean matrices over these areas have then been computed and used to simulate the 4 region dataset shown in Fig. 2 according to the method described in [3] (pp. 114115) that we briefly recall here:

Compute the square root of the covariance to simulate , where .

Simulate complex random vectors with zero mean and identity covariance matrix.

Compute singlelook vectors .

Compute the multilook covariance by combining the independent samples according to equation (5).
It may be noted that those simulations only allow to generate homogeneous areas according to the fullydeveloped speckle model. The case of nonGaussian and spatially correlated data are beyond the scope of this paper and left for future work.
Additionally, to simulate deterministic targets, we have measured the coherency matrix over a building layover in the ESAR image and added it in the form of parallel lines (close to the top left corner of our synthetic image) without speckle. Assuming that the strong double bounce scattering from the building dominates the overall response, the effects of speckle can be neglected in this case. Note that this is only an approximate way to simulate this kind of target. Nevertheless, we consider it is sufficient to study the spatial properties of our filter in this simulation. In fact, purely deterministic targets are represented by rank 1 matrices (only one scattering mechanism is present). This ideal model is usually not found in experimental data due to the presence of speckle and the effect of multilooking. However, it is important to evaluate the filter in such an extreme case since some targets may be close to rank 1. Therefore, we study this case in a separate simulation, that is described in section VB.
This simulated dataset allows a quantitative evaluation of the bilateral filter (that we refer to as BLF in the following) for the three distances , and described in sections IVB and IVC. The BLF has 3 parameters: , and the number of iterations . Since has the role of a spatial window, which has already been extensively studied for speckle filtering [22], we fix it to and use an window in all experiments. In the following, we focus our study on and .
We have chosen quantitative performance measures according to the following criteria:

closeness to the original signal in the mean square sense,

amount of smoothing in homogeneous areas and

preservation of edges.
To evaluate the performance according to criterion 1, we propose the following quantity:
(20) 
where is the number of pixels in the image, is the dimension of the matrices, is the true covariance at pixel and its estimate by the BLF. This error is the per element root mean square reconstruction error, measuring the average deviation between the true covariances and the ones estimated by the filter.
To measure the amount of smoothing over homogeneous areas (criterion 2), we use the estimated equivalent number of looks (ENL) [3]
(21) 
where and
are the estimated mean and variance over a considered channel. Since our filter performs averages of covariance matrices and does not process the elements independently, all the channels are processed identically. Consequently, the ENL is only measured on the
channel of the filtered image and over a manually selected homogeneous area, assuming the other channels have been applied the same smoothing.Finally, an important feature of speckle filters that is sometimes overlooked is the preservation of edges (criterion 3). Here, an edge is considered to be present if two adjacent pixels belong to a different class in the groundtruth used for simulations. It is difficult to evaluate edge preservation on PolSAR images due to the complex nature of the data. We have found experimentally that the measure as defined previously was not a relevant measure of edge reconstruction. This is due to the fact that pixels from homogeneous areas largely outnumber the ones located around edges and tend to dominate this quantity. Some authors [23, 24] propose to use gradients to evaluate edge preservation. However, gradients can only be computed on intensity and are highly sensitive to noise. To measure the polarimetric reconstruction of edges, we propose to use the same measure as in equation (20) restricted to the pixels that are located on both sides of discontinuities
(22) 
where is the binary mask of such pixels indexed by their location , is defined as in eq. (20) and is equal to if its argument is true and otherwise. The mask is set equal to one if the pixel at location has at least one neighbour belonging to a different class in the ground truth. Thus, we eliminate the influence of pixels located inside homogeneous areas and focus on edge reconstruction.
Fig. 3 represents the curves corresponding to the variation of those values for different values of BLF parameters. We chose the range of values experimentally considering the fact that the curves for have a local minimum for certain combinations of and . We limited the number of displayed curves for visual purpose. However, the evolution of the curves allow the reader to understand the effect of those parameters on the different measures. On the one hand, for low values of the local minimum of is attained for a high number of iterations. On the other hand, higher values of require less iterations to attain the minimum but lead to an optimum with a higher error. This choice of values also allows to show that the global error decreases with the number of iterations and that all the curves converge to the same value. Regarding the , increasing both and leads to higher values. It may be observed that both Riemannian distances and behave in a very similar way. To obtain comparable curves for all the distances, we show the results of the BLF with with a different set of parameters. This difference can be explained in an intuitive way by considering the simpler case of two real positive scalars: varies as the sum of the ratio and its inverse whereas and behave as the logarithm of the ratio. Therefore, the two kinds of distances have different dynamics. The curves representing the ENL show that the BLF with Riemannian distances achieve a higher amount of smoothing in homogeneous areas.
According to this experiment, we select a set of parameters that correspond to a high ENL and good edge preservation while allowing a reasonable computational cost (computation time is mostly linked to the number of iterations.) For and we choose and instead of and because the gain in and is small and running the filter with 9 iterations would dramatically increase the compuational cost. Nevertheless, it is shown in the following that those slightly suboptimal values are sufficient so that the BLF outperforms the other compared methods. Therefore, the selected parameter set for and is and while for we choose and .
Fig. 4 allows a visual comparison of the filtered image for the selected parameters and other parameter values. These images show that for the selected set of parameters, all edges are wellpreserved in the image with a high level of smoothing on homogeneous areas. On the contrary, low values of and result in a poor speckle reduction while too high values lead to blurry edges in low contrast areas. This visually confirms the measurements given in Fig. 3
VB Validation on rank1 targets.
It is important to verify that the method allows to avoid the filtering of purely deterministic targets. As explained in the previous section, ideal deterministic targets are unlikely to be present in experimental data. However, the presence of matrices that are close to rank 1 may result in numerical instability due to the use of matrix inversion and logarithm in the computation of the distances. Therefore, we have simulated a small image containing two ideal deterministic targets that correspond to rank one covariance matrices. A dot representing a purely trihedral target (the ideal representation of a corner reflector) is simulated with a diagonal matrix and a line modelling an ideal dihedral scatterer with matrix . The background is a homogeneous area with speckle simulated by the previously described method (see section VA). We filtered this image with the BLF using the parameters selected in the previous section. Fig. 5 shows the results obtained for the BLF with that involves both matrix inversion and logarithm. Similar results were obtained with the other versions of the filter. It can be noted that, thanks to the modification described in section IVD to handle deterministic targets, the background is well smoothed and both pure targets are preserved. This results in correct parameters for both background and rank 1 targets.
VC Comparison with other filters
We first have visually compared the three versions of the BLF for the previously selected parameters with the boxcar filter with a window, the refined Lee filter with pixel wide oriented windows and the IDAN filter for a neighborhood of pixels. We used implementations available in the RAT [25] free software for the last two filters.
Fig. 6 shows a visual comparison of the results obtained with the different methods on the synthetic image. For a better visualization areas of interest have been extracted from the original image. While the boxcar filter blurs edges, degrading polarimetric information around boundaries, the BLF leads to sharp edges and smoothed estimates in homogeneous areas. Again, we notice that the Riemannian distance based filters perform better than the KullbackLeibler divergence based filter, that produces noisier estimates. This could be explained by the use of matrix logarithms in the Riemannian distances that make these quantities less sensitive to speckle. Although the refined Lee and IDAN filters allow a good preservation of edges, they also result in a noisier image. In the case of the Lee filter this is due to the use of a finite number of oriented windows. In the case of IDAN, this is due to the hard decision of including or not a pixel in the average. This effect is not present in the BLF filtered images thanks to the use of continuously varying weights.
Table I displays the values of the previously defined quality measures for all the filtering methods. All versions of the BLF lead to lower error than the other filters in terms of global reconstruction and edge preservation. Moreover, the BLF outperforms the other methods in terms of ENL. Furthermore, it may be observed that the BLF with Riemannian distances leads to the best estimates.
To evaluate the effects of the BLF on polarimetric information, it is useful to apply a polarimetric decomposition to the filtered data and to compare it to the ideal values. In this work, we choose to consider the entropyalpha () decomposition [2]. Here, denotes the mean alpha angle. For PolSAR data, it is also important that the filter does not lead to a systematic bias in the values of the covariance elements. On table II we display the mean values computed over 4 manually selected homogeneous areas (shown in Fig. 2). We also show the mean of angle and entropy . The estimated values are close to the real ones for all coherency elements, including offdiagonal, as well as for the mean alpha angle and polarimetric entropy. It may be observed that the various versions of the BLF compare well with the other techniques. This study also shows that both Riemannian versions of the BLF outperform the KullbackLeibler based version.
VD Experimental data
Fig. 7 shows the result of the BLF with , applied to the experimental Oberpfaffenhofen dataset introduced at the beginning of the section. For a better visualization, only a cropped area of the original image is shown. This experiment was also done with that showed similar performance. As a preprocessing, a presumming of 2 pixels in both azimuth and range has been applied. The filtering parameters are the ones retained in the previous sections. It can be noted that the filter is able to preserve complex spatial structures in urban and forested areas and strongly reduces speckle on homogeneous areas. Fig. 8 shows a zoom on a subimage with several homogeneous areas and compares the BLF with the boxcar, refined Lee and IDAN filters. The BLF results in sharp edges and very smooth homogeneous areas. As observed in simulations, the refined Lee and IDAN filters also preserve edges but they also lead to a noisier image that can be a problem for applications such as segmentation or edge detection. The drawback of the BLF is a slight oversmoothing of very low contrast edges. A representation of the estimated alpha entropy parameters is shown in the bottom row, using a Hue/Saturation color map. If the parameters look roughly the same for all the filters, the BLF filtered image has a smoother aspect than the other ones. Fig. 9 shows a zoom on a corner reflector. It may be observed that the BLF allows a good reconstruction of the reflector, and performs a better background smoothing than the other methods.
Fig. 10 shows representations of the point density in the plane for boxcar, BLF, Lee and IDAN. Both BLF and boxcar lead to more compact scatter plots than Lee and IDAN filter, suggesting a higher amount of averaging. However, the BLF point cloud contains more very low entropy points than the boxcar one, suggesting its ability to preserve those points whereas boxcar tends to mix them with other higher entropy points.
Finally, we show results obtained with the same set of parameters on the more recent FSAR SBand demo set provided by DLR (Fig. 11). The BLF achieves a highquality reconstruction of the different types of structures in the image. In particular, it may be observed that buildings and roads are visually well preserved, as well as spatial texture in forested areas.
Vi Conclusion
In this paper, we have developed a new speckle filter for PolSAR data based on an iterative version of the bilateral filter. To take advantage of the radiometric similarities between covariance matrices in a local neighbourhood, we have used weights based on suitable distances.
The KullbackLeibler statistical distance has been compared to two Riemannian distances, that are defined on the manifold of Hermitian positive definite matrices.
We have applied this filter to simulated and experimental data. An indepth evaluation has been performed over synthetic data. A study of the results obtained by all versions of the filter with different sets of parameters has been achieved. This study allowed a selection of suitable parameters for the comparison with other filters. The BLF showed good performance for spatial preservation of edges and deterministic targets. Moreover, thanks to an iterative scheme, the BLF outperforms other methods in terms of ENL. The preservation of the polarimetric parameters is also achieved. The best performance has been obtained with the Riemannian distance based filters. Finally, our filter has been applied to experimental data and achieved a highquality restoration of complex spatial structures.
In this work we have considered only fullydeveloped speckle that corresponds to a Gaussian distribution of target vectors. To apply such an approach to very high resolution data, it should also be studied in the case of nonGaussian data. Effects of the filter on spatially correlated texture should also be considered in future work. Since, the use of Riemannian based similarities showed the best performance, it would also be interesting to investigate the relation between Riemannian means and the parameters of probability distributions such as Wishart or Kdistribution. In order to avoid a manual tuning of the parameters, automatic selection should be addressed in future research. Finally, the performance of the filter in terms of sharp edge reconstruction and smoothing of homogeneous areas suggest to investigate its use in applications such as segmentation, edge detection and machine learning.
Acknowledgement
The authors would like to thank the anonymous reviewers for their constructive comments that greatly helped improving the quality of the paper.
The authors also would like to thank DLR for providing the FSAR and ESAR datasets. This work was supported by the German Science Foundation DFG, under project No. HE 2459/151
References
 [1] J. W. Goodman, “Some fundamental properties of speckle,” J. Opt. Soc. Am., vol. 66, no. 11, pp. 1145–1150, Nov 1976.
 [2] S. Cloude and E. Pottier, “A review of target decomposition theorems in radar polarimetry,” IEEE Trans. Geosci. Remote Sens., vol. 34, no. 2, pp. 498–518, 1996.
 [3] J. Lee and E. Pottier, Polarimetric Radar Imaging: From Basics to Applications, ser. Optical Science and Engineering. Taylor & Francis, 2009.
 [4] J. Lee, M. Grunes, and G. de Grandi, “Polarimetric SAR speckle filtering and its implication for classification,” IEEE Trans. Geosci. Remote Sens, vol. 37, no. 5, pp. 2363 –2373, 1999.
 [5] G. Vasile, E. Trouvé, J. Lee, and V. Buzuloiu, “Intensitydriven adaptiveneighborhood technique for polarimetric and interferometric SAR parameters estimation,” IEEE Trans. Geosci. Remote Sens., vol. 44, no. 6, pp. 1609–1621, 2006.
 [6] C. LópezMartínez and X. Fàbregas, “Modelbased polarimetric SAR speckle filter,” IEEE Trans. Geosci. Remote Sens, vol. 46, no. 11, pp. 3894–3907, 2008.
 [7] S. Foucher, G. Farage, and G. Benie, “Speckle filtering of PolSAR and PolInSAR images using tracebased partial differential equations,” in IGARSS, 2006, pp. 2545 –2548.
 [8] A. AlonsoGonzález, C. LópezMartínez, and P. Salembier, “Filtering and segmentation of polarimetric SAR data based on binary partition trees,” Geoscience and Remote Sensing, IEEE Transactions on, vol. PP, no. 99, pp. 597 –609, 2012.
 [9] C. Deledalle, F. Tupin, and L. Denis, “Polarimetric SAR estimation based on nonlocal means,” in IGARSS, 2010, pp. 2515–2518.
 [10] J. Chen, Y. Chen, W. An, Y. Cui, and J. Yang, “Nonlocal filtering for polarimetric SAR data: A pretest approach,” IEEE Trans. Geosci. Remote Sens., vol. 49, no. 5, pp. 1744 –1754, may 2011.
 [11] C. Tomasi and R. Manduchi, “Bilateral filtering for gray and color images,” in ICCV, 1998, p. 839.
 [12] F. Durand and J. Dorsey, “Fast bilateral filtering for the display of highdynamicrange images,” in SIGGRAPH, 2002, pp. 257–266.
 [13] S. Kullback and R. A. Leibler, “On information and sufficiency,” Annals of Mathematical Statistics, vol. 22, pp. 49–86, 1951.
 [14] X. Pennec, P. Fillard, and N. Ayache, “A Riemannian framework for tensor computing,” Int. J. of Computer Vision, vol. 66, no. 1, pp. 41–66, 2006.
 [15] F. Barbaresco, “Interactions between symmetric cone and information geometries: BruhatTits and Siegel spaces models for high resolution autoregressive doppler imagery,” in Emerging Trends in Visual Computing, 2009, pp. 124–163.
 [16] V. Arsigny, P. Fillard, X. Pennec, and N. Ayache, “Logeuclidean metrics for fast and simple calculus on diffusion tensors,” Magnetic Resonance in Medicine, vol. 56, no. 2, pp. 411–421, 2006.
 [17] N. R. Goodman, “Statistical analysis based on a certain multivariate complex gaussian distribution (an introduction),” The Annals of Mathematical Statistics, vol. 34, pp. 152–177, 1963.
 [18] S. Paris, P. Kornprobst, and J. Tumblin, Bilateral Filtering. Hanover, MA, USA: Now Publishers Inc., 2009.
 [19] D. Barash, “A fundamental relationship between bilateral filtering, adaptive smoothing and the nonlinear diffusion equation,” IEEE Trans. Pattern Anal. Mach. Intell, vol. 24, no. 6, pp. 844–847, 2002.
 [20] A. Buades, B. Coll, and J. Morel, “Nonlocal means denoising,” Image Processing On Line, 2011.
 [21] W. Press, S. Teukolsky, W. Vetterling, and B. Flannery, Numerical Recipes in C. Cambridge: Cambridge University Press, 1992.
 [22] R. Touzi, “A review of speckle filtering in the context of estimation theory,” IEEE Trans. Geosci. Remote Sens., vol. 40, no. 11, pp. 2392 – 2404, nov 2002.
 [23] Y. Yu and S. Acton, “Speckle reducing anisotropic diffusion,” IEEE Trans. Image Process., vol. 11, no. 11, pp. 1260 – 1270, nov 2002.
 [24] S. Foucher and C. LopezMartinez, “An evaluation of PolSAR speckle filters,” in IGARSS, vol. 4, july 2009, pp. IV–845 –IV–848.
 [25] A. Reigber, M. Neumann, and T. Dar, “Radar Tools (RAT),” http://radartools.berlios.de/.
Comments
There are no comments yet.