Evaluation and Bias Removal of Evaluation and Bias Removal of - - PowerPoint PPT Presentation

Evaluation and Bias Removal of Evaluation and Bias Removal of Multi Multi-

• Look Effect on

Look Effect on Entropy/Alpha /Anisotropy Entropy/Alpha /Anisotropy (H/

(H/α α/A /A) )

POLINSAR 2009 WORKSHOP POLINSAR 2009 WORKSHOP

26-29 January 2009 ESA-ESRIN, Frascati (ROME), Italy

Jong Jong-

• Sen Lee*, Thomas Ainsworth

Sen Lee*, Thomas Ainsworth Naval Research Laboratory Naval Research Laboratory Washington DC 20375, USA Washington DC 20375, USA * CSRSR, National Central University, Taiwan * CSRSR, National Central University, Taiwan

J.S. Lee, et al., “Evaluation and bias removal of multi-look effect on Entropy/Alpha/Anisotropy in polarimetric SAR decomposition,” IEEE Transactions on Geoscience and Remote Sensing, October 2008

INTRODUCTION INTRODUCTION

• Entropy/Anisotropy/Alpha (H/A/α): Widely

applied and effective for PolSAR data analysis.

• Geophysical parameter estimation:
• Anisotropy —Surface roughness
• Entropy and Alpha — Soil Moisture
• Entropy — Biomass
• Accurate H/A/α estimation require averaging:
• Underestimate Entropy
• Overestimate Anisotropy
• Alpha?

MOTIVATION MOTIVATION

• Evaluate multi-look asymptotic behavior of

H and A by a simple simulation technique:

• The effect of number of looks on

Averaged α .

• The H/A/α bias problem for L-band and

X-band data

• Devise a bias removal scheme:
• Entropy
• Anisotropy
• Alpha
• C. Lopez-Martinez, E. Pottier and S.R. Cloude, “Statistical assessment of eigenvector based target decomposition

Theorems in radar polarimetry,” IEEE Trans. Geoscuence and Remote Sensing, September 2005.

H / A / H / A / α α DECOMPOSITION DECOMPOSITION

[ ]

U3 = ⎡ ⎣ ⎢ ⎢ ⎢ ⎤ ⎦ ⎥ ⎥ ⎥ cos cos cos sin sin sin sin( sin( sin( ( ) ( ) ( ) ( )cos( )e ( )cos( )e ( )cos( )e )sin ( )e ) sin ( )e )sin( )e

1 2 3 1 1 j 1 2 2 j 2 3 3 j 3 1 1 j 1 2 2 j 2 3 3 j 3

α α α α β α β α β α β α β α β

δ δ δ γ γ γ

T T T

u u u u u u T

* 3 3 3 * 2 2 2 * 1 1 1

] [ λ λ λ + + =

3 SCATTERING PROCESSES

[ ]

k S S S S S

XX YY XX YY XY T

= + − 1 2 2

TARGET VECTOR

[ ] [ ]

T N k k N T

i i T i N i i N

= ⋅ =

= =

∑ ∑

1 1

1 1 *

LOCAL ESTIMATE OF THE COHERENCY MATRIX

S.R. CLOUDE E. POTTIER

3 3 2 2 1 1

P P P α α α α + + =

ENTROPY

) P ( log P H

3 1 i i 3 i

∑

=

− =

3 ROLL INVARIANT PARAMETERS ANISOTROPY

3 2 3 2

A λ λ λ λ + − = α PARAMETER

H / A / H / A / α α DECOMPOSITION DECOMPOSITION

Pi

i k k

=

=

∑

λ λ

1 3

MULTI-LOOK (AVERAGING) EFFECT ON H/A/α:

• UNDERESTIMATE OF H
• OVERESTIMATE OF A
• α DEPENDS ON SCATTERING MECHANISM, BUT

HAS LESS EFFECT.

• C. Lopez-Martinez recommends 9x9 for H and 11x11 for A

H

Original (4 looks) 5x5 9x9

ENTROPY(H) VERSUS MULTI ENTROPY(H) VERSUS MULTI-

• LOOKING

LOOKING

|HH-VV|, |HV|, |HH+VV|

Freeman and Durden Decomposition

Original (4 looks) 5x5 9x9

Anisotropy Anisotropy and and α α VERSUS MULTI VERSUS MULTI-

• LOOKING

LOOKING

A

ALPHA

α

Aniso- tropy

E-SAR L-BAND POLSAR DATA OF OBERPFAFFENHOFEN Freeman/Durden Decomposition

SIMULATION AREA SELECTION SIMULATION AREA SELECTION

Urban (D. B.) Forest (Volume) Grass (Surface)

For a given <T>, simulate single-look complex data:

SIMULATION PROCEDURE SIMULATION PROCEDURE

• 1. Compute
• 3. Form a single-look complex vector

2 / 1

T T T T

T = * 2 / 1 2 / 1

) (

• 2. Simulate a complex random vector, , CN(0,I)

ν

ν

2 / 1

T u =

4. Compute a n look covariance matrix, 5. Compute the n look H/A/α Verification:

T n n

uu n T

* 1

1∑ =

T T vv E T uu E

T T T

= =

* 2 / 1 2 / 1

) ]( [ ] [

EIGENVALUE ESTIMATION EIGENVALUE ESTIMATION

• Mean value of n-look estimation
• Dominant Eigenvalue changes little
• Forest (Volume) eigenvalues within 4 dB

High Entropy, Low Anisotropy

• Urban (D.B.): two dominant scatterings

Medium Entropy, High Anisotropy

• Grass (Surface): One dominant

Low Entropy, Medium Anisotropy

ENTROPY ESTIMATION ENTROPY ESTIMATION

• Entropy is underestimated (1000 samples)
• Rate of increase changes at 5x5 looks.
• 7x7 looks sufficient for Entropy
• 5x5 may severely underestimate entropy.
• Remove bias is possible: 3x3 looks and
• above. 5x5 looks is recommended.
• Boxcar average includes mixed media.

ANISOTROPY ESTIMATION ANISOTROPY ESTIMATION

• Anisotropy is overestimated.
• Rate of increase changes at 5x5 looks.
• For forest (volume) area, impossible to
• btain accurate estimate – very small.
• 9x9 looks sufficient for Anisotropy
• Remove bias is more difficult: 5x5 and
• above. 7x7 is recommended.

AVERAGED ALPHA ESTIMATION AVERAGED ALPHA ESTIMATION

is affected less by multi-looking, except for Surface

• Underestimate or overestimate.
• Peculiar asymptotic behavior for Volume
• Sufficient using 5x5 independent looks
• Bias compensation is required for Surface

α

AVERAGED ALPHA ESTIMATION AVERAGED ALPHA ESTIMATION

1

α

Peculiar asymptotic behavior for Volume

• High entropy – 3 scattering mechanisms
• decreases asymptotically
• and increase asymptotically

3 2 1

λ λ λ ≈ ≈

1

α

2

α

3

α

ENTROPY BIAS REMOVAL ENTROPY BIAS REMOVAL

) ( ) ( ∞ = H n H R

The Ratio Scattering mechanism dependent?

ENTROPY BIAS REMOVAL ENTROPY BIAS REMOVAL

) ( ) ( ∞ = H n H R

) ( ) ( ) ( ˆ n R n H n H

• =

The Ratio Entropy bias removal:

• Ratios are based on Complex

Wishart statistical model

• Linear relation
• Identical linear relation for other

L-Band PolSAR systems

• Identical linear relation for other

frequencies (X-band C-band and P-band).

(Surface) (Urban) (Forest)

ENTROPY BIAS REMOVAL ENTROPY BIAS REMOVAL

AIRSAR L-band

L-Band ALOS/ PALSAR from Tomakomai, Japan

ENTROPY BIAS REMOVAL ENTROPY BIAS REMOVAL

L-band AIRSAR San Francisco

ENTROPY BIAS REMOVAL ENTROPY BIAS REMOVAL

The Ratio

L-band PISAR from Tsukuba, Japan

ENTROPY BIAS REMOVAL ENTROPY BIAS REMOVAL

X-band PISAR from Tsukuba, Japan

ENTROPY BIAS REMOVAL ENTROPY BIAS REMOVAL

AIRSAR L-band

) ( ) ( ∞ = H n H R

ENTROPY BIAS REMOVAL ENTROPY BIAS REMOVAL

) ( ) ( ∞ = H n H R The Ratio

• X-band and L-band have the same ratio

Frequency independent.

• The ratio depends on the number of looks and

) (n H

ENTROPY BIAS REMOVAL ENTROPY BIAS REMOVAL

3x3 Average Bias removed

) 5 . 7 ( H

Entropy bias removal:

) 5 . 7 ( ˆ H

The 3x3 averaged data have ENL=7.5

) ( ) ( ) ( ˆ n R n H n H =

ENTROPY BIAS REMOVAL ENTROPY BIAS REMOVAL

, 5x5 Average Bias removed

) 18 ( ˆ H

) ( ) ( ) ( ˆ n R n H n H =

The 5x5 averaged data have ENL=18

) 18 ( H

ANISOTROPY BIAS REMOVAL ANISOTROPY BIAS REMOVAL

) ( ) ( ∞ = A n A RA

The Ratio Surface and volume scattering requires bias removal. Bias removal requires 49 looks data.

E-SAR L-Band

Anisotropy Anisotropy Bias Bias Removal Removal

Ratio for Volume is very high cause the problem for bias removal. To compensate for bias for Volume class, create a ratio curve for Volume class alone. Do the same for Surface.

ANISOTROPY BIAS REMOVAL ANISOTROPY BIAS REMOVAL

Anisotropy, 7x7 average Bias compensated 13x13 average

Surface and volume scattering requires bias removal.

ALPHA BIAS EVALUATION ALPHA BIAS EVALUATION

) ( ) ( ∞ α α = α n

The Ratio Bias is very small Only the surface category requires bias removal

E-SAR L-Band

ALPHA BIAS EVALUATION ALPHA BIAS EVALUATION

Only the surface scattering category may require bias removal for surface geophysical parameter estimation E-SAR L-Band 5x5 Average 13x13 Average 7x7 Average

SUMMARY SUMMARY

• Evaluated the asymptotic behaviors of H and A as a

function of the number of looks.

• Bias removal:
• Entropy – Robust linear characteristics
• Linear relation is independent radar frequency

and radar systems

• 25 independent looks with bias removal
• 49 without bias removal
• Anisotropy
• Bias removal is required for surface and volume
• 49 independent looks with bias removal
• Alpha
• Bias is small
• Bias removal for entropy and alpha for surface is

required for soil moister estimation

α

PIXEL CORRELATION EFFECT PIXEL CORRELATION EFFECT

Assess the effect of over-sampling at 25%, 50% and 100% in both range and azimuth

• Pixel correlations: 0.234 (25%), 0.415 (50%), 0.636 (100%)
• At 50%, 5x5 looks has the same underestimate in Entropy of 3x5 (0%).
• At 100%, very high number of looks is required to reduce bias.

PIXEL CORRELATION EFFECT PIXEL CORRELATION EFFECT

Alpha Angle of λ1

For Forest (Volume): for forest has the same peculiar effect Over sampling affects

α

1

α

MIXED PIXEL EFFECT MIXED PIXEL EFFECT

Boxcar averaging includes mixed pixels

• Mixed pixels affect surface scattering

pixels:

• Increase Entropy
• Decrease Anisotropy
• Change Averaged Alpha