Research Article  Open Access
Modified Hybrid Freeman/Eigenvalue Decomposition for Polarimetric SAR Data
Abstract
Because of the rapid advancement of the airborne sensors and spaceborne sensors, large volumes of fully polarimetric synthetic aperture radar (PolSAR) data are available, but they are too complex to interpret difficultly. In this paper, a modified hybrid Freeman/eigenvalue decomposition method for the coherency matrix derived from the fully PolSAR sensors is proposed. The proposed modified hybrid Freeman/eigenvalue decomposition uses a real unitary transformation on the coherency matrix to release correlations between the copolarized term and cross polarized term, and the scattering models are derived from eigenvectors of the coherency matrix with reflection symmetry condition. The anisotropy and entropy are used to determine whether the volume scattering component is derived from the manmade structures or not. Moreover, the scattering powers from the proposed hybrid Freeman/eigenvalue decomposition are all nonnegative values. Fully PolSAR data on San Francisco acquired by AIRSAR sensor are used in the experiments to prove the efficacy of the proposed decomposition.
1. Introduction
Since 1985, the first fully polarimetric AIRSAR at Lband was launched by the Jet Propulsion Laboratory (JPL) [1] which started a rapid advancement stage of PolSAR sensors, such as the wellknown airborne sensors: EMISAR [2] by the Technical University of Denmark and ESAR [3] by DLR of Germany, and the famous spaceborne sensors: TerraSARX [4] by DLR of Germany and Radarsat2 [5] by CSA of Canada and so forth. Large volumes of the fully PolSAR data need to be interpreted. Polarimetric target decomposition is an important and useful tool for understanding the PolSAR data [6] by separating received measurements into basic scattering mechanisms. The scattering mechanisms of the PolSAR media are analyzed for the purpose of parameter inversion, terrain classification, and so forth. Currently, eigenvectorbased decompositions and modelbased decomposition methods are commonly used on the secondorder statistics matrix of the PolSAR data. Cloude and Pottier developed the most notable eigenvectorbased decomposition, that is, Entropy/Alpha method [7]. The classical modelbased decomposition was FreemanDurden decomposition (FDD) developed by Freeman and Durden [8], which decomposed the coherency matrix of PolSAR data into three components: surface scattering, doublebounce scattering, and volume scattering with the reflection symmetry condition. The reflection symmetry condition implies that the correlation between copolarized term and cross polarized term is zero. FreemanDurden decomposition is used in various applications since it is easy to understand and accomplish.
When applying FreemanDurden decomposition on the real PolSAR data, some scattering powers are negative those most frequently occur in the doublebounce scattering powers. To solve this deficiency, various modified methods have been developed [9–15]. Cloude improved FreemanDurden decomposition via setting surface scattering model and doublebounce scattering model to be orthogonal [14], that is, wellknown hybrid Freeman/eigenvalue decomposition. After the rotation of the coherency matrix, the scattering powers derived from hybrid Freeman/eigenvalue decomposition are effective at avoiding negative values. Singh improved the hybrid Freeman/eigenvalue decomposition by using different volume scattering models [15] for vegetation areas and oriented structures.
In this paper, a modified version of hybrid Freeman/eigenvalue decomposition for PolSAR data is proposed by using different volume scattering models for the vegetation areas and manmade structures. Moreover, the surface scattering model and doublebounce scattering model are defined as the eigenvectors of the coherency matrix of PolSAR data. The eigenvector with scattering angle greater than denotes the doublebounce scattering model, while the other scattering angle less than represents the surface scattering model. The volume scattering model for the manmade structures is also derived from the eigenvectors of the coherency matrix. We will show how the eigenspace of the coherency matrix enables the proposed Freeman/eigenvalue decomposition and solve the scattering powers as the linear combinations of eigenvalues. In addition, the scattering powers are all nonnegative values.
The rest of this paper is organized as follows. The hybrid Freeman/eigenvalue decomposition is summarized in Section 2. The proposed decomposition is presented in Section 3. Results and discussion of experiments performed on the real PolSAR data are provided in Section 4. Section 5 presents our conclusions.
2. Freeman/Eigenvalue Decomposition
In this section, the Freeman/eigenvalue decomposition is simply introduced for the integrity of this paper. The detail content has been shown in [14].
For monostatic PolSAR sensor with basis, a Pauli vector is used to represent the single look PolSAR data as (1) and the coherency matrix is used to represent the multilook PolSAR data as (2) as follows:where denotes a transposition operator, implies a complex conjugation processing, and denotes the multilook processing.
In the modelbased decomposition, such as FreemanDurden decomposition [8] and hybrid Freeman/eigenvalue decomposition [14], the measured coherency matrix of the PolSAR data is expanded into three components, surface scattering, doublebounce scattering, and volume scattering as follows:where , , and are surface scattering model, doublebounce scattering model, and volume scattering model, respectively.
In hybrid Freeman/eigenvalue decomposition, those three scattering models are defined in (4), (5), and (6). , , and are the corresponding scattering powers in the proper order: and imply the type of the scattering model; namely, denotes the surface scattering, while represents the doublebounce scattering. In [14], the condition is set to reduce the number of unknowns. is the volume parameter. If , hybrid Freeman/eigenvalue decomposition becomes equivalent to FreemanDurden decomposition [8]. The cross polarized term only exists in the volume scattering model, so the volume scattering power is solved as follows:
Since the rank of either surface scattering model or doublebounce scattering model is equal to 1, the corresponding scattering powers and are the eigenvalues of as (10), and and can be solved as (11).
and are obtained as follows: Whichever of and is smaller is always set to zero, and then volume parameter can be estimated as
It enables solving for the scattering powers and by substituting (9) into (11):
3. Proposed Decomposition
3.1. Rotation of the Coherency Matrix
To reduce the power of cross polarized term , a real unitary transformation has been proposed to be implemented on the coherency matrix [12, 16] before decomposition as (12).
is the angle of rotation about the radar line of sight. To minimize , two angles are solved as
After the rotation of the coherency matrix, the imagery element of is equal to zero.
3.2. Eigenvalue Decomposition of the Rotated Coherency Matrix
We assume that the reflection symmetry condition (i.e., ) holds; then each pixel for multilook PolSAR data is represented as
The average coherency matrix after rotation is decomposed [7] aswhere , , and are the eigenvalues of the rotated coherency matrix , and . Because of the reflection symmetry condition and the rotation of the coherency matrix, consists of the corresponding eigenvectors , , and aswhere eigenvectors , , and are orthogonal unit vectors, .
3.3. Scattering Models
Because in hybrid Freeman/eigenvalue decomposition, the condition is set, and, in eigenvectors of the coherency matrix, it can be seen that , we can draw the conclusions: the eigenvectors and contain a surface scattering component and a doublebounce scattering component, and the scattering angles and are equivalent to the angles and of the hybrid Freeman/eigenvalue decomposition. The eigenvector or in with is shown in (17) which denotes the surface scattering targets, and the other eigenvector with is a doublebounce scattering target (18) under the reflection symmetry condition:
If , we obtain , the corresponding eigenvector as (17), and surface scattering model whose form is the same as the one in (4); and doublebounce scattering model is identical to (5). In other cases, if , , surface scattering model as (4), , and doublebounce scattering as (5).
3.3.1. Volume Scattering Model from the ManMade Structures
The entropy and anisotropy are defined in the eigenspace of the coherency matrix to measure the randomness of the targets [7, 17]:
For the manmade structures, usually the value of entropy is large (usually ), and anisotropy is often larger than 0.5 [7, 17]. Based on these, we defined the volume scattering model as
The rank of the volume scattering model is equal to two. The volume scattering model is designed as a distributed target that is diffused from the surface scatterers whose scattering angle and orientation angle and the oriented objects whose scattering angle and orientation angle . For the real PolSAR data, the oriented objects can be modeled as angle and orientation angle because of the reflection symmetry condition. But the surface scatterers usually are not the same as “ and ,” so the surface scattering case is relaxed as (17), and the volume scattering model for the manmade structures is defined as
3.3.2. Volume Scattering Model from the Vegetation Areas
The volume scattering model is used as a unit diagonal matrix whose rank is equal to 3 which has been proposed by An et al. [12]. Because the unit diagonal matrix has the largest entropy () and least anisotropy (), the unit diagonal matrix is also used for the vegetation areas in this paper:
3.4. Freeman/Eigenvalue Decomposition of the Rotated Coherency Matrix
For vegetation areas, the cross polarized term only presents in the volume scattering model, and volume scattering power can be solved as . In the eigenspace of the rotated coherency matrix under the reflection symmetry condition, the cross polarized term is one of the eigenvalues. After the real unitary transformation, reaches the minimum value, usually, . So the volume scattering power is solved as
The surface scattering power and doublebounce scattering power are the eigenvalues of . According to the eigenvalue decomposition of the rotated coherency matrix with reflection symmetry condition, can be expanded asFrom (23), , the surface scattering power , and doublebounce scattering power are solved as (25) while , , and are solved as (26):
For manmade structures, that is, and [7, 17], similar to the case of vegetation areas, the volume scattering power is solved asExpand the coherency matrix after the real unitary transformation into eigenspace:If , it can be seen that, , , then the surface scattering model , and the volume scattering model . In addition, , , and doublebounce scattering model . From the fifth line of (28), we can solve surface scattering powers and doublebounce powers asIf , then , , , and , . The volume scattering model is . From the last line of (28), we can solve surface scattering powers and doublebounce powers as
Because of and because the scattering powers are solved as (23), (25), and (26) or (27), (29), and (30), the scattering powers are all nonnegative values.
4. Experimental Study
To prove the efficacy of proposed Freeman/eigenvalue decomposition, the experiments are conducted on the Lband fully PolSAR data of San Francisco which were acquired by NASA/JPL ARISAR. The spatial resolution and range resolution are both about 10 m and the radar incidence angle is from 5° to 60°. The used PolSAR data is open access from the internet [18]. The original image is shown in Figure 1, with the selected regions, which are used in the later tests. The size of the used PolSAR image in these experiments is 900 × 1024. Before expanding the coherency matrix, the Sigma filter is used to process the speckle of the PolSAR data.
In order to demonstrate the efficacy of the proposed hybrid Freeman/eigenvalue decomposition, the compared methods used FreemanDurden decomposition [8] (FDD 1), FreemanDurden decomposition with the rotation of the coherency matrix [12] (FDD 2), hybrid Freeman/eigenvalue decomposition [14] (HFED 1), and hybrid Freeman/eigenvalue decomposition with extended volume scattering model [15] (HFED 2). The PolSAR data are decomposed into three components: surface scattering power (blue), doublebounce scattering power (red), and volume scattering power (green) in Figure 2.
(a)
(b)
(c)
(d)
(e)
From Figures 2(a)–2(e), three main terrain types in the used data, that is, ocean areas, city blocks, and forests are reconstructed well by the scattering powers derived from these five decompositions. For further analysis, we compare the scattering powers in three selected regions in Figure 1. Those regions are Region 1, Region 2, and Region 3, respectively. The sizes of these regions are 60 × 100, 70 × 70, and 60 × 100 in proper order. The types of ground truth are ocean areas, city blocks, and forests. The mean values of surface scattering power mean_, doublebounce scattering power mean_, and volume scattering power mean_ in the three regions are listed in Tables 1, 2, and 3. The mean scattering powers are all normalized by the total scattering powers (i.e., ). In Region 1, it can be seen that mean_ given by the proposed decomposition is 0.9634, which is about 4.8%, 1.7%, 3.1%, and 3.1% higher than the other four decompositions, respectively. The efficacy of all the five decompositions on Region 1 (sea) provides the excellent performance. In Region 2, the average doublebounce scattering power mean_ of the proposed decomposition also outperforms the other methods. It is about 3.9%, 2.8%, 3.0%, and 0.6% larger than these four decompositions in proper sequence. But, in Region 3, mean_ given by HFED 1 is the best.



We classified the selected regions into three classes. The rule is simple; that is, the max scattering power determines the label of the pixel. If the surface scattering power is the biggest in the three scattering powers, the pixel is labeled “surface scattering class.” It can be seen that the true labels in Region 1, Region 2, and Region 3 are surface scattering class, doublebounce scattering class, and volume scattering class, respectively. The accuracy of the selected zones is listed in Table 4. It can be seen, in Region 1, that surface scattering powers of these decompositions are all the biggest ones, and, in Region 3, these classification accuracies are all higher than 95%. In contrast to the other decompositions, the average scattering powers of the proposed method in Region 3 are not very good, but the classification accuracy is higher than 98%. Moreover, in Region 2, the proposed method obtains the best results, as large as 95.92%, which are about 59.43%, 11.72%, 11.10%, and 9.29% higher than FDD 1, FDD 2, HEFD 1, and HEFD 2, respectively.

5. Conclusions
In this paper, a novel version of hybrid Freeman/eigenvalue decomposition for polarimetric SAR data is proposed. Three conclusions can by drawn about the proposed method. Firstly, the eigenvectors of the rotated coherency matrix with reflection symmetry condition are used as surface scattering or doublebounce scattering models. Secondly, in contrast to conventional FreemanDurden decompositions, the derived scattering powers are all nonnegative values. Thirdly, the volume scattering model is determined by the entropy and anisotropy of the coherency matrix. Experimental results have proved the efficacy of the proposed decomposition.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported in part by the National High Technology Research and Development Program (863 Program) of China (no. 2013CB329402), the National Natural Science Foundation of China (nos. 61173092, 61377011, and 61272282), and the Fund for Foreign Scholars in University Research and Teaching Programs (no. B07048).
References
 J. J. van Zyl, R. Carande, Y. Lou, T. Miller, and K. Wheeler, “The NASA/JPL threefrequency polarimetric AIRSAR system,” in Proceedings of the International Geoscience and Remote Sensing Symposium (IGARSS '92), vol. 1, pp. 649–651, Houston, Tex, USA, May 1992. View at: Publisher Site  Google Scholar
 E. L. Christensen, N. Skou, J. Dall et al., “EMISAR: an absolutely calibrated polarimetric L and Cband SAR,” IEEE Transactions on Geoscience and Remote Sensing, vol. 36, no. 6, pp. 1852–1865, 1998. View at: Publisher Site  Google Scholar
 R. Horn, “The DLR airborne SAR project ESAR,” in Proceedings of the International Geoscience and Remote Sensing Symposium “Remote Sensing for a Sustainable Future” (IGARSS '96), vol. 3, pp. 1624–1628, Lincoln, Neb, USA, May 1996. View at: Publisher Site  Google Scholar
 S. Buckreuss, R. Werninghaus, and W. Pitz, “The German satellite mission TerraSARX,” in Proceedings of the IEEE Radar Conference (RADAR '08), pp. 1–5, Rome, Italy, May 2008. View at: Publisher Site  Google Scholar
 http://www.asccsa.gc.ca/eng/satellites/radarsat2.
 S. R. Cloude and E. Pettier, “A review of target decomposition theorems in radar polarimetry,” IEEE Transactions on Geoscience and Remote Sensing, vol. 34, no. 2, pp. 498–518, 1996. View at: Publisher Site  Google Scholar
 S. R. Cloude and E. Pottier, “An entropy based classification scheme for land applications of polarimetric SAR,” IEEE Transactions on Geoscience and Remote Sensing, vol. 35, no. 1, pp. 68–78, 1997. View at: Publisher Site  Google Scholar
 A. Freeman and S. L. Durden, “A threecomponent scattering model for polarimetric SAR data,” IEEE Transactions on Geoscience and Remote Sensing, vol. 36, no. 3, pp. 963–973, 1998. View at: Publisher Site  Google Scholar
 Y. Yamaguchi, T. Moriyama, M. Ishido, and H. Yamada, “Fourcomponent scattering model for polarimetric SAR image decomposition,” IEEE Transactions on Geoscience and Remote Sensing, vol. 43, no. 8, pp. 1699–1706, 2005. View at: Publisher Site  Google Scholar
 Y. Yamaguchi, A. Sato, W.M. Boerner, R. Sato, and H. Yamada, “Fourcomponent scattering power decomposition with rotation of coherency matrix,” IEEE Transactions on Geoscience and Remote Sensing, vol. 49, no. 6, pp. 2251–2258, 2011. View at: Publisher Site  Google Scholar
 O. Antropov, Y. Rauste, and T. Hame, “Volume scattering modeling in PolSAR decompositions: study of ALOS PALSAR data over boreal forest,” IEEE Transactions on Geoscience and Remote Sensing, vol. 49, no. 10, pp. 3838–3848, 2011. View at: Publisher Site  Google Scholar
 W. An, Y. Cui, and J. Yang, “Threecomponent modelbased decomposition for polarimetric sar data,” IEEE Transactions on Geoscience and Remote Sensing, vol. 48, no. 6, pp. 2732–2739, 2010. View at: Publisher Site  Google Scholar
 G. Singh, Y. Yamaguchi, and S.E. Park, “General fourcomponent scattering power decomposition with unitary transformation of coherency matrix,” IEEE Transactions on Geoscience and Remote Sensing, vol. 51, no. 5, pp. 3014–3022, 2013. View at: Publisher Site  Google Scholar
 S. R. Cloude, Polarisation: Applications in Remote Sensing, Oxford University Press, London, UK, 2009.
 G. Singh, Y. Yamaguchi, S.E. Park, Y. Cui, and H. Kobayashi, “Hybrid freeman/eigenvalue decomposition method with extended volume scattering model,” IEEE Geoscience and Remote Sensing Letters, vol. 10, no. 1, pp. 81–85, 2013. View at: Publisher Site  Google Scholar
 F. Xu and Y.Q. Jin, “Deorientation theory of polarimetric scattering targets and application to terrain surface classification,” IEEE Transactions on Geoscience and Remote Sensing, vol. 43, no. 10, pp. 2351–2364, 2005. View at: Publisher Site  Google Scholar
 J. S. Lee and E. Pottier, Polarimetric Radar Imaging: From Basics to Applications, CRC Press, Boca Raton, Fla, USA, 2009.
 https://earth.esa.int/web/polsarpro/datasets.html.
Copyright
Copyright © 2015 Shuang Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.