# Modified Hybrid Freeman/Eigenvalue Decomposition for Polarimetric SAR Data

Authors: Shuang Zhang; Shuang Wang; Bo Chen
Journal: Journal of Sensors
 (2015)

---

## 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 man-made 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.


---

## Body

## 1. Introduction

Since 1985, the first fully polarimetric AIRSAR at L-band was launched by the Jet Propulsion Laboratory (JPL) [1] which started a rapid advancement stage of PolSAR sensors, such as the well-known airborne sensors: EMISAR [2] by the Technical University of Denmark and E-SAR [3] by DLR of Germany, and the famous spaceborne sensors: TerraSAR-X [4] by DLR of Germany and Radarsat-2 [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, eigenvector-based decompositions and model-based decomposition methods are commonly used on the second-order statistics matrix of the PolSAR data. Cloude and Pottier developed the most notable eigenvector-based decomposition, that is, Entropy/Alpha method [7]. The classical model-based decomposition was Freeman-Durden decomposition (FDD) developed by Freeman and Durden [8], which decomposed the coherency matrix of PolSAR data into three components: surface scattering, double-bounce 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. Freeman-Durden decomposition is used in various applications since it is easy to understand and accomplish.When applying Freeman-Durden decomposition on the real PolSAR data, some scattering powers are negative those most frequently occur in the double-bounce scattering powers. To solve this deficiency, various modified methods have been developed [9–15]. Cloude improved Freeman-Durden decomposition via setting surface scattering model and double-bounce scattering model to be orthogonal [14], that is, well-known 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 man-made structures. Moreover, the surface scattering model and double-bounce scattering model are defined as the eigenvectors of the coherency matrix of PolSAR data. The eigenvector with scattering angleα greater than π/4 denotes the double-bounce scattering model, while the other scattering angle α less than π/4 represents the surface scattering model. The volume scattering model for the man-made 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 Section2. 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{H,V} 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:(1)k→p=12SHH+SVVSHH-SVV2SHVt(2)T=k→p·k→p∗t=T11T12T13T12∗T22T23T13∗T23∗T33,where t denotes a transposition operator, ∗ implies a complex conjugation processing, and  denotes the multilook processing.In the model-based decomposition, such as Freeman-Durden decomposition [8] and hybrid Freeman/eigenvalue decomposition [14], the measured coherency matrix of the PolSAR data is expanded into three components, surface scattering, double-bounce scattering, and volume scattering as follows:(3)T=msTs+mdTd+mvTv,where Ts, Td, and Tv are surface scattering model, double-bounce scattering model, and volume scattering model, respectively.In hybrid Freeman/eigenvalue decomposition, those three scattering models are defined in (4), (5), and (6). ms, md, and mv are the corresponding scattering powers in the proper order:(4)Ts=cos2αscosαssinαse-jϕs0cosαssinαsejϕssin2αs0000,αs≤π4,(5)Td=cos2αdcosαdsinαde-jϕd0cosαdsinαdejϕdsin2αd0000,αd>π4,(6)Tv=Fs00010001αs and αd imply the type of the scattering model; namely, αs≤π/4 denotes the surface scattering, while αd>π/4 represents the double-bounce scattering. In [14], the condition αs+αd=π/2 is set to reduce the number of unknowns. Fs is the volume parameter. If Fs=2, hybrid Freeman/eigenvalue decomposition becomes equivalent to Freeman-Durden decomposition [8]. The cross polarized term T33 only exists in the volume scattering model, so the volume scattering power md is solved as follows:(7)mv=T33.Since the rank of either surface scattering modelTs or double-bounce scattering model Td is equal to 1, the corresponding scattering powers ms and md are the eigenvalues of TSD as (10), and ms and md can be solved as (11).α
s and αd are obtained as follows: (8)αd,s=cos-11+T12T22-T33-md,s2-1/2.Whichever of ms and md is smaller is always set to zero, and then volume parameter Fs can be estimated as(9)Fs=T11·T22-T33-T122T33·T22-T33.It enables solving for the scattering powersms and md by substituting (9) into (11):(10)TSD=T-mvTv=cosαsinα0-sinαcosα0001·ms000md0000·cosα-sinα0sinαcosα0001(11)md,s=T11+T22-Fs+1T33T11-T22-Fs-1T332+4T122hhhh±T11-T22-Fs-1T332+4T122·2-1(12)Tθ=1000cos2θsin2θ0-sin2θcos2θT1000cos2θ-sin2θ0sin2θcos2θ.


## 3. Proposed Decomposition

### 3.1. Rotation of the Coherency Matrix

To reduce the power of cross polarized termT33, 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 T33, two angles are solved as(13)θ=14tan-12ReT23T22-T33.After the rotation of the coherency matrix, the imagery element ofT23 is equal to zero.


### 3.2. Eigenvalue Decomposition of the Rotated Coherency Matrix

We assume that the reflection symmetry condition (i.e.,SHHSHV∗≈SVVSHV∗≈0) holds; then each pixel for multilook PolSAR data is represented as(14)Tθ=T11θT12θ0T12∗θT22θ000T33θ.The average coherency matrix after rotation is decomposed [7] as(15)Tθ=U·λ1000λ2000λ3·U∗t,where λ1, λ2, and λ3 are the eigenvalues of the rotated coherency matrix T(θ), and λ1≥λ2≥λ3≥0. Because of the reflection symmetry condition and the rotation of the coherency matrix, U consists of the corresponding eigenvectors k1, k2, and k3 as(16)U=k1k2k3=cosα1cosα20sinα1ejγ1sinα2ejγ20001,where eigenvectors k1, k2, and k3 are orthogonal unit vectors, α1+α2=π/2.


### 3.3. Scattering Models

Because in hybrid Freeman/eigenvalue decomposition, the conditionαs+αd=π/2 is set, and, in eigenvectors of the coherency matrix, it can be seen that α1+α2=π/2, we can draw the conclusions: the eigenvectors k1 and k2 contain a surface scattering component and a double-bounce scattering component, and the scattering angles α1 and α2 are equivalent to the angles αs and αd of the hybrid Freeman/eigenvalue decomposition. The eigenvector k1 or k2 in U with α≤π/4 is shown in (17) which denotes the surface scattering targets, and the other eigenvector with π/4≤α≤π/2 is a double-bounce scattering target (18) under the reflection symmetry condition:(17)ks=cosαssinαsejγs0,αs≤π4,(18)kd=cosαdsinαdejγd0,π4≤αd≤π2.Ifα1≤π/4, we obtain αs=α1, the corresponding eigenvector ks=k1 as (17), and surface scattering model Ts=k1k1∗t whose form is the same as the one in (4); αd=π/2-α1=α2 and double-bounce scattering model Td=k2k2∗t is identical to (5). In other cases, if α2≤π/4, αs=α2, surface scattering model Ts=k2k2∗t as (4), αd=π/2-α2=α1, and double-bounce scattering Td=k1k1∗t as (5).


#### 3.3.1. Volume Scattering Model from the Man-Made Structures

The entropy and anisotropy are defined in the eigenspace of the coherency matrix to measure the randomness of the targets [7, 17]:(19)H=-∑I=1,2,3pilog3pi,withpi=λiλ1+λ2+λ3,A=λ2-λ3λ2+λ3.For the man-made structures, usually the value of entropy(H) is large (usually H>0.7), and anisotropy is often larger than 0.5 [7, 17]. Based on these, we defined the volume scattering model as(20)T=12100000001.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α=0 and orientation angle β=0 and the oriented objects whose scattering angle α=90 and orientation angle β=90. For the real PolSAR data, the oriented objects can be modeled as angle α=90 and orientation angle β=90 because of the reflection symmetry condition. But the surface scatterers usually are not the same as “α=0 and β=0,” so the surface scattering case is relaxed as (17), and the volume scattering model for the man-made structures is defined as(21)Tv_2=12ksks∗t+k3k3∗t.


#### 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 (H=1) and least anisotropy (A=0), the unit diagonal matrix is also used for the vegetation areas in this paper:(22)Tv_1=13100010001.


### 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 asmv=3T33θ. In the eigenspace of the rotated coherency matrix under the reflection symmetry condition, the cross polarized term T33θ is one of the eigenvalues. After the real unitary transformation, T33θ reaches the minimum value, usually, T33θ=λ3. So the volume scattering power is solved as(23)mv=3λ3.The surface scattering powerms and double-bounce scattering power md are the eigenvalues of TSD. According to the eigenvalue decomposition of the rotated coherency matrix with reflection symmetry condition, TSD can be expanded as(24)TθSD=msTs+mdTd=Tθ-mvTv_1=Tθ-3λ3Tv_1=λ1-λ3k1k1∗t+λ2-λ3k2k2∗t.From (23), α1≤π/4, the surface scattering power ms, and double-bounce scattering power md are solved as (25) while α2≤π/4, ms, and md are solved as (26):(25)ms=λ1-λ3,md=λ2-λ3,(26)ms=λ2-λ3,md=λ1-λ3.For man-made structures, that is,H<0.7 and A>0.5 [7, 17], similar to the case of vegetation areas, the volume scattering power is solved as(27)mv=2λ3.Expand the coherency matrix after the real unitary transformation into eigenspace:(28)Tθ=msTs+mdTd+mvTv=TθSD+mvTv_1=TθSD+2·λ3·12ksks∗t+k3k3∗t=λ1k1k1∗t+λ2k2k2∗t+λ3k3k3∗t=λ1-λ2k1k1∗t+λ2k2k2∗t+2·λ3·12k1k1∗t+k3k3∗t=λ1k1k1∗t+λ2-λ3k2k2∗t+2·λ3·12k2k2∗t+k3k3∗t.If α1≤π/4, it can be seen that, αs=α1, ks=k1, then the surface scattering model Ts=k1k1∗t, and the volume scattering model Tv_2=(1/2)k1k1∗t+k3k3∗t. In addition, αd=α2, kd=k2, and double-bounce scattering model Td=k2k2∗t. From the fifth line of (28), we can solve surface scattering powers ms and double-bounce powers md as(29)ms=λ1-λ3,md=λ2.If α2≤π/4, then αs=α2, ks=k2, Ts=k2k2∗t, and αd=α1, Td=k1k1∗t. The volume scattering model is Tv_2=(1/2)k2k2∗t+k3k3∗t. From the last line of (28), we can solve surface scattering powers ms and double-bounce powers md as(30)ms=λ2-λ3,md=λ1.Because ofλ1⩾λ2⩾λ3⩾0 and because the scattering powers are solved as (23), (25), and (26) or (27), (29), and (30), the scattering powers are all nonnegative values.


## 3.1. Rotation of the Coherency Matrix

To reduce the power of cross polarized termT33, 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 T33, two angles are solved as(13)θ=14tan-12ReT23T22-T33.After the rotation of the coherency matrix, the imagery element ofT23 is equal to zero.


## 3.2. Eigenvalue Decomposition of the Rotated Coherency Matrix

We assume that the reflection symmetry condition (i.e.,SHHSHV∗≈SVVSHV∗≈0) holds; then each pixel for multilook PolSAR data is represented as(14)Tθ=T11θT12θ0T12∗θT22θ000T33θ.The average coherency matrix after rotation is decomposed [7] as(15)Tθ=U·λ1000λ2000λ3·U∗t,where λ1, λ2, and λ3 are the eigenvalues of the rotated coherency matrix T(θ), and λ1≥λ2≥λ3≥0. Because of the reflection symmetry condition and the rotation of the coherency matrix, U consists of the corresponding eigenvectors k1, k2, and k3 as(16)U=k1k2k3=cosα1cosα20sinα1ejγ1sinα2ejγ20001,where eigenvectors k1, k2, and k3 are orthogonal unit vectors, α1+α2=π/2.


## 3.3. Scattering Models

Because in hybrid Freeman/eigenvalue decomposition, the conditionαs+αd=π/2 is set, and, in eigenvectors of the coherency matrix, it can be seen that α1+α2=π/2, we can draw the conclusions: the eigenvectors k1 and k2 contain a surface scattering component and a double-bounce scattering component, and the scattering angles α1 and α2 are equivalent to the angles αs and αd of the hybrid Freeman/eigenvalue decomposition. The eigenvector k1 or k2 in U with α≤π/4 is shown in (17) which denotes the surface scattering targets, and the other eigenvector with π/4≤α≤π/2 is a double-bounce scattering target (18) under the reflection symmetry condition:(17)ks=cosαssinαsejγs0,αs≤π4,(18)kd=cosαdsinαdejγd0,π4≤αd≤π2.Ifα1≤π/4, we obtain αs=α1, the corresponding eigenvector ks=k1 as (17), and surface scattering model Ts=k1k1∗t whose form is the same as the one in (4); αd=π/2-α1=α2 and double-bounce scattering model Td=k2k2∗t is identical to (5). In other cases, if α2≤π/4, αs=α2, surface scattering model Ts=k2k2∗t as (4), αd=π/2-α2=α1, and double-bounce scattering Td=k1k1∗t as (5).


### 3.3.1. Volume Scattering Model from the Man-Made Structures

The entropy and anisotropy are defined in the eigenspace of the coherency matrix to measure the randomness of the targets [7, 17]:(19)H=-∑I=1,2,3pilog3pi,withpi=λiλ1+λ2+λ3,A=λ2-λ3λ2+λ3.For the man-made structures, usually the value of entropy(H) is large (usually H>0.7), and anisotropy is often larger than 0.5 [7, 17]. Based on these, we defined the volume scattering model as(20)T=12100000001.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α=0 and orientation angle β=0 and the oriented objects whose scattering angle α=90 and orientation angle β=90. For the real PolSAR data, the oriented objects can be modeled as angle α=90 and orientation angle β=90 because of the reflection symmetry condition. But the surface scatterers usually are not the same as “α=0 and β=0,” so the surface scattering case is relaxed as (17), and the volume scattering model for the man-made structures is defined as(21)Tv_2=12ksks∗t+k3k3∗t.


### 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 (H=1) and least anisotropy (A=0), the unit diagonal matrix is also used for the vegetation areas in this paper:(22)Tv_1=13100010001.


## 3.3.1. Volume Scattering Model from the Man-Made Structures

The entropy and anisotropy are defined in the eigenspace of the coherency matrix to measure the randomness of the targets [7, 17]:(19)H=-∑I=1,2,3pilog3pi,withpi=λiλ1+λ2+λ3,A=λ2-λ3λ2+λ3.For the man-made structures, usually the value of entropy(H) is large (usually H>0.7), and anisotropy is often larger than 0.5 [7, 17]. Based on these, we defined the volume scattering model as(20)T=12100000001.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α=0 and orientation angle β=0 and the oriented objects whose scattering angle α=90 and orientation angle β=90. For the real PolSAR data, the oriented objects can be modeled as angle α=90 and orientation angle β=90 because of the reflection symmetry condition. But the surface scatterers usually are not the same as “α=0 and β=0,” so the surface scattering case is relaxed as (17), and the volume scattering model for the man-made structures is defined as(21)Tv_2=12ksks∗t+k3k3∗t.


## 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 (H=1) and least anisotropy (A=0), the unit diagonal matrix is also used for the vegetation areas in this paper:(22)Tv_1=13100010001.


## 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 asmv=3T33θ. In the eigenspace of the rotated coherency matrix under the reflection symmetry condition, the cross polarized term T33θ is one of the eigenvalues. After the real unitary transformation, T33θ reaches the minimum value, usually, T33θ=λ3. So the volume scattering power is solved as(23)mv=3λ3.The surface scattering powerms and double-bounce scattering power md are the eigenvalues of TSD. According to the eigenvalue decomposition of the rotated coherency matrix with reflection symmetry condition, TSD can be expanded as(24)TθSD=msTs+mdTd=Tθ-mvTv_1=Tθ-3λ3Tv_1=λ1-λ3k1k1∗t+λ2-λ3k2k2∗t.From (23), α1≤π/4, the surface scattering power ms, and double-bounce scattering power md are solved as (25) while α2≤π/4, ms, and md are solved as (26):(25)ms=λ1-λ3,md=λ2-λ3,(26)ms=λ2-λ3,md=λ1-λ3.For man-made structures, that is,H<0.7 and A>0.5 [7, 17], similar to the case of vegetation areas, the volume scattering power is solved as(27)mv=2λ3.Expand the coherency matrix after the real unitary transformation into eigenspace:(28)Tθ=msTs+mdTd+mvTv=TθSD+mvTv_1=TθSD+2·λ3·12ksks∗t+k3k3∗t=λ1k1k1∗t+λ2k2k2∗t+λ3k3k3∗t=λ1-λ2k1k1∗t+λ2k2k2∗t+2·λ3·12k1k1∗t+k3k3∗t=λ1k1k1∗t+λ2-λ3k2k2∗t+2·λ3·12k2k2∗t+k3k3∗t.If α1≤π/4, it can be seen that, αs=α1, ks=k1, then the surface scattering model Ts=k1k1∗t, and the volume scattering model Tv_2=(1/2)k1k1∗t+k3k3∗t. In addition, αd=α2, kd=k2, and double-bounce scattering model Td=k2k2∗t. From the fifth line of (28), we can solve surface scattering powers ms and double-bounce powers md as(29)ms=λ1-λ3,md=λ2.If α2≤π/4, then αs=α2, ks=k2, Ts=k2k2∗t, and αd=α1, Td=k1k1∗t. The volume scattering model is Tv_2=(1/2)k2k2∗t+k3k3∗t. From the last line of (28), we can solve surface scattering powers ms and double-bounce powers md as(30)ms=λ2-λ3,md=λ1.Because ofλ1⩾λ2⩾λ3⩾0 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 L-band 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.Figure 1
Original image of San Francisco.In order to demonstrate the efficacy of the proposed hybrid Freeman/eigenvalue decomposition, the compared methods used Freeman-Durden decomposition [8] (FDD 1), Freeman-Durden 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 ms (blue), double-bounce scattering power md (red), and volume scattering power mv (green) in Figure 2.Figure 2
Decompositions of AIRSAR data on San Francisco withms for blue, md for red, mv for green. (a) FDD 1, (b) FDD 2, (c) HFED 1, (d) HFED 2, and (e) proposed decomposition.
(a)
(b)
(c)
(d)
(e)From Figures2(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 powermean_ms, double-bounce scattering powermean_md, and volume scattering powermean_mv 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., ms+md+mv). In Region 1, it can be seen thatmean_ms 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 double-bounce scattering powermean_md 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_mv given by HFED 1 is the best.Table 1
Means of dominated scattering powers in Region 1.
mean_
m
s
mean_
m
d
mean_
m
v
FDD 1
0.9192
0.0013
0.0794
FDD 2
0.9475
0.0033
0.0492
HFED 1
0.9348
0.0164
0.0488
HFED 2
0.9348
0.0164
0.0488
Proposed method
0.9634
0.0053
0.0313Table 2
Means of dominated scattering powers in Region 2.
mean_
m
s
mean_
m
d
mean_
m
v
FDD 1
0.1082
0.4307
0.4611
FDD 2
0.1676
0.5017
0.3307
HFED 1
0.1686
0.5008
0.3306
HFED 2
0.3245
0.5137
0.1618
Proposed method
0.2404
0.5158
0.2438Table 3
Means of dominated scattering powers in Region 3.
mean_
m
s
mean_
m
d
mean_
m
v
FDD 1
0.0314
0.1535
0.8151
FDD 2
0.0247
0.1555
0.8198
HFED 1
0.0673
0.1098
0.9228
HFED 2
0.0996
0.1097
0.7904
Proposed method
0.2380
0.0754
0.6865We 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, double-bounce scattering class, and volume scattering class, respectively. The accuracy of the selected zones is listed in Table4. 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.Table 4
Classification results (%).
Region 1
Region 2
Region 3
FDD 1
100
36.49
95.97
FDD 2
100
84.65
99.32
HFED 1
100
84.82
99.47
HFED 2
100
86.63
95.50
Proposed method
100
95.92
98.05


## 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 double-bounce scattering models. Secondly, in contrast to conventional Freeman-Durden 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.



---

*Source: 168720-2015-10-20.xml*