Non Stationary Texture Analysis from Polarimetric SAR Data D’Hondt, Olivier , Ferro-Famil, Laurent , Pottier, Eric IETR Laboratory, The University of Rennes 1, CNRS UMR 6164, S.A.P.H.I.R Team, Image and Remote Sensing Group Campus de Beaulieu - Bat 11D, 263 Avenue Général Leclerc, CS 74205, 35042 Rennes Cedex, France Tel/Fax : (+33) (0)22323 5716 E-mail : [email protected]

Abstract In this paper, a new model for SAR intensity texture is presented, based on the local spatial structure of the observed image. Texture is represented by a correlated non stationary random process resulting from the local convolution of Gaussian white noise with 2D Gaussian kernels. An estimation method for the texture parameters is then presented, taking into account the presence of speckle and the stochastic multiplicative nature of the intensity texture itself. Next the method is validated by means of simulated speckled intensity texture data. Finally, the texture parameters are estimated over real polarimetric SAR data.

1

Introduction

Data acquired by polarimetric SAR are directly related to physical properties of natural media. A more complete analysis of such data can be obtained by studying the spatial structure variation of the radar scene, caused by the fluctuations of reflectivity between neighbouring resolution cells. Due to the presence of speckle, these variations are usually described by their first and second order statistics. First order statistics were extensively studied with the multiplicative model [1], leading to the well-known K distribution [2]. This distribution was firstly derived for conventional SAR and extended to the multidimensional case [3] in order to deal with polarimetric datasets. Autocorrelation function (ACF) is the common descriptor for the second order statistics of intensity [4]. Generally, ACF is assumed stationary and isotropic. In this paper, a model that takes into account spatial fluctuations of the ACF and local texture orientation is presented and validated on simulated and real polarimetric SAR data.

2

Spatial Fluctuation Model

The spatial information about a stationary random field X (s ) on a discrete grid s = [ x , y ]T is given by its 2D autocorrelation for the lag t : RXX (t ) = E [ X (s )X (s + t ) ] (1) where s and s + t are the two pixel locations considered, and E () denotes the expectation. In general

it is more convenient to deal with the autocovariance, i.e. the autocorrelation of the centred variable: C XX (t ) = E [ ( X (t ) − µX )( X (t + s ) − µX ) ] (2) and the corresponding normalised ACF : C (t ) ρXX (t ) = XX2 (3) σX where µX , σX2 are respectively the mean value and the variance of X .

2.1

Anisotropic Gaussian Process

Gaussian oriented texture can be obtained by convolving a zero-mean Gaussian white noise of unit variance with a bidimensional kernel K (s ) , s = [ x , y ]T . Then the stationary autocorrelation of the resulting field is given by: RXX (t ) = ∑ s ∈Z2 K (s )K (s − t ) (4) In this paper, we assume that SAR data texture correlation is a 2D anisotropic Gaussian function. Relation (4) indicates that this process can be generated by a filter with the following Gaussian impulse response:  1  u(x , y )2 v(x , y )2   A K (x , y ) = exp  −  +   (5) 2  2  σu 2πσu σv σv2   where A is a normalisation constant ensuring that the output texture has unit variance, σ u and σ v are the standard deviations of the kernel along directions u and v defined as:  u   cos θ sin θ   x    =   . (6)  v   − sin θ cos θ   y  Such an ACF can be represented by one of its constant value elliptical contour. By analogy with the one

dimensional case the ellipse resulting from the intersection of the normalised ACF with the plane of constant altitude exp(−1) is chosen. The ellipse axes may then be associated to the principal texture correlation lengths τu = 2σu and τv = 2σv . Texture orientation is defined by θ , the angle formed by the principal axis of the ellipse and a horizontal reference.

2.2

Gaussian Process with NonStationary Covariance

In the present study, non stationary process are represented by a field of kernels whose paralemers vary as a function of the spatial location s . The autocorrelation expression is then more general than the for the stationary case [5]: RXX (s, s ′) = ∑ u ∈Z2 Ks (u )Ks ′ (u ) (7) To simulate such data, we generate a family of Gaussian kernels {Ks (.)} with location dependent parameters. Each pixel of our non-stationary process is the local convolution of a zero mean white noise X with the corresponding kernel Ks : Z (s ) = ∑ u ∈ Z2 Ks (u )X (u ) (8) This method permits to generate a wide class of space varying textures, characterized by three functions of the spatial location τu , τv and θ defined in the above section.

3

Texture Model for SAR Data

In this paper we consider the statistics of single-look intensity data. To capture the intrinsic spatial variability of SAR data we adopt, for a given textured zone of intensity channel, the following model [6]: I (x , y ) = µIT (x , y )F (x , y ) (9) T (x , y ) is the correlated texture field containing the spatial variability of interest. This variable is normalised and µT = E [T ] = 1 . F (x , y ) is the pixel-to-pixel fluctuation due to speckle noise and µI the mean intensity of the area. For model simplicity, speckle contribution is assumed to be spatially uncorrelated, of unit mean µF = E [ F ] = 1 and for single-look data, of unit variance. As explained in the following, it is preferred to work with the centred normalised variable I − µI IN = (10) µI whose statistics are directly related to texture since its standard deviation i.e. its autocovariance at lag zero is equal to the intensity coefficient of variation. Furthermore, second-order statistics of I N are linked to those of T by a simple relation since its covariance can be written as: 1 C NN (t ) = 2 ( E [ I (s )I (s + t ) ] − µI2 ) µI (11) = RFF (t )RTT (t ) − 1

Finally, recalling that speckle is assumed to be uncorrelated of unit mean and variance, texture statistics can be estimated directly from data autocovariance: C (t ) + 1 (12) CTT (t ) = NN −1 1 + δ0,t where δ0,t is the Kronecker delta function. Relation (12) implies that for all non-zero lags CTT is equal to C NN , and for zero lag: G G C NN (0) − 1 2 (13) CTT (0) = σT = . 2 To obtain the estimate of texture covariance, it is only needed to correct the zero lag value of C NN according to (13).

4

Analysis of Simulated and Real Data

4.1

Covariance Estimation method

From the above sections, an estimation procedure for SAR intensity spatial texture covariance is derived. In a first stage the normalised intensity defined in the above section is computed on the region of interest with (10) where the mean intensity value is obtained on the basis of a square sliding window of dimension NS . Then, for each pixel, we compute the estimate of normalised intensity correlation with a spatial lag ranging from 0 to Lmax .The assumption is made that on a small window of size N S this correlation does not vary significantly and can be approximated by the stationary sample covariance estimate: N N 1 S S CˆNN (i, j ) = 2 ∑ ∑ I N (p, q )I N (p + i, q + j ) N S p =1 q = 1 (14) where the indices i and j vary between −Lmax and Lmax . The zero lag value of this estimate is linearly related to the product of texture variance and speckle contribution. Relation (13) is used to correct this value and keep only the covariance corresponding to texture contribution CˆTT (i, j ) . Assuming that this correlation is a 2D Gaussian function, the following stage is the estimation of the local correlation lengths τu,v and local orientation θ . This is done by calculating the two-by-two geometrical covariance matrix of CˆTT (i, j ) in the directions i and j. The eigenvalues of this matrix indicate the lengths of the ellipse axes, and the eigenvector argument corresponding to the higher eigenvalue is the orientation angle. To keep only relevant information, the covariance function is applied a threshold at exp(−1) so as to consider its central part only.

4.2

Estimation on simulated data

To evaluate the method, speckled non-stationary covariance intensity is simulated. Spatial information is modelled by a 128x128 pixels random Gaussian process with non stationary 2D Gaussian covariance generated as described in section 2.2. In order to focus on local orientation estimation, the generated field has fixed correlation lengths τu = 5; τv = 1.5 and θ

process is multiplied by an exponential unit mean noise in order to simulate speckled single-look SAR intensity data. The proposed Gaussian kernel estimation procedure is then applied to a 60x60 zone using a sliding window of size 31, and the covariance is computed for a maximum lag of 7 pixels in both x and y directions. Results are shown in figure 2. The algorithm tends to underestimate the correlation lengths, however the mean value for the estimated axis lengths ratio is close to the real one. Moreover, the trend of the estimated angle along a line is close to the simulated slope. 70

angle (degrees)

60

50

40

30 true value estimate linear regression on estimate

20 Fig1: Simulated non stationary Gaussian texture

10

0

10

20

30 X

40

50

60

Fig 3: Variation of the estimated angle along a horizontal line

4.3

(a)

Estimation on real polarimetric intensity channels

The method is applied to L-band polarimetric SIR-C SAR data acquired over the French Alps during April 1994. As the developed model deals with intensity, the processed data are the three polarimetric intensity channels S HH 2 , SVV 2 and S HV 2 . The texture variance [6] is first computed on each channel in order to

(b) Fig 2: Texture analysis. (a) simulated Gaussian kernels field (b) estimated Gaussian kernels from speckled simulation

angle varies linearly from 0 to 90 degrees. To obtain intensity texture data, the generated random field is squared and normalised so that the resulting process is unit mean and variance. Then this correlated texture Fig 4: Span image of the selected area on a SIR-C image. The arrow shows the 10x10 analysed zone.

find an analysis area with significant spatial variations. Then a 10 by 10 zone is selected for analysis, and the estimation method is applied with the same parameters than in the above section. The results of kernel estimation on the three channels are shown on figure 5. From the analysis, it is obvious that polarimetric intensity channels have different spatial structure. Kernels non stationarity is more pronounced on the VV channel, while HH global orientation is very different from the other channels. Moreover the HV channel shows a most important directionality.

5

Conclusion

A new model is introduced, that extracts the spatial variability information from polarimetric SAR data texture. It has been shown that 2D Gaussian kernels are significant descriptors for local texture orientation and directivity. Subsequently, an algorithm to estimate texture parameters of SAR intensity data has been presented. This process takes into account the presence of speckle and assumes a multiplicative nature for texture variations. The method has been tested and validated over simulated speckled intensity data, and has been employed to retrieve texture parameters of real polarimetric SAR data. This analysis has revealed that the spatial structure of texture was dependent on the polarimetric channel.

6

(a)

References (b)

[1] Lombardo, P. , Oliver, C. J.: Estimation of Texture Parameters in K-distributed Clutter, IEE Proc.-Radar, Sonar Navig., Vol 141, No 4, August 1994 [2] Jakeman E., On the Statistics of K-distributed Noise, Journal of Physics A, Vol 13, pp. 31-48 [3] Yueh S.R., Kong J.A., Jao J.K., Shin R.T., Zebker H.A.,Le Toan T.and Ottl H., KDistribution and polarimetric terrain clutter. Polarimetric remote sensing, edited by J.A. Kong, Amesterdam: Elsevier, 1990 [4] Lombardo, P. , Oliver, C. J.: Estimating the Correlation Properties of K-distributed SAR Clutter, IEE Proc.-Radar, Sonar Navig., Vol 142, No 4, August 1995 [5] Higdon, D. , Swall, J. , Kern, J. ,: NonStationary Spatial Modelling, Bayesian Statistics 6, Oxford University Press, 1998 [6] Ulaby, F. T. , Kouyate, F. , Brisco, B. , Lee Williams, T. H.: Textural Information in SAR Images, IEEE Trans. Geos. Rem. Sensing, Vol GE24,No. 2, March 1986.

(c) Fig 5: Estimated kernels for polarimetric data.(a) HH channel, (b) HV channel, (c) VV channel.

Sample Page

Stationary Covariance. In the present study, non stationary process are repre- sented by a field of kernels whose paralemers vary as a function of the spatial ...

623KB Sizes 1 Downloads 277 Views

Recommend Documents

Sample Title Page
This manuscript has been read and accepted for the. Graduate Faculty in Art ...... Manuals by Albrecht Dürer: The Small Passion and the Engraved Passion: Iconography ...... Mixtec leaders ostentatiously wear Spanish clothes and subserviently pay tri

Sample Title Page
Mixtec Devotional and Doctrinal Books: Format, Content, and Illustrations .... 7. Codex Yanhuitlan, pl. 17. Source: María Teresa Sepúlveda y Herrera, Códice de Yanhuitlán: Estudio preliminar (Mexico City: INAH; Benemérita. Universidad ... Franci

Sample Report Cover Page-1
imply its endorsement, recommendation, or favoring by the LASC, it's Board .... Adopt a policy of annual or bi-annual bicycle and pedestrian counts to collect ...

Sample Report Cover Page-1
The survey was primarily conducted online with additional survey responses gathered at ... the existing network that connect cyclists to public transit, employment centers, ... o Bicycle counts also will provide data to measure the effectiveness of .

Five Sample Armies page 1.pdf
Page 1 of 1. Page 1 of 1. Five Sample Armies page 1.pdf. Five Sample Armies page 1.pdf. Open. Extract. Open with. Sign In. Main menu. Displaying Five Sample Armies page 1.pdf. Page 1 of 1.

Page 1 AIG Sample Criterio Indicators Grades 3-12 STANDARD ...
students whose rotings indicote d need for extensive differentiction. Students must meet four Out of five Criterio listed below. Observable Student Behaviors.

Page 1 Page 2 Page 3 Page 4 Page 5 Page 6 Page 7 Page 8 Page 9 ...
Professional Growth (20%). 1. Education (10%). > Doctoral Degree – 10 pts. > Complete Academic requirement for Doctoral Degree – 7 pts. > Master's Degree ...

Page 1 Page 2 Page 3 Page 4 Page 5 Page 6 Page 7 Page 8 Page 9 ...
content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print,.

Page 1 Page 2 Page 3 Page 4 Page 5 Page 6 Page 7 Page 8 Page 9 ...
8.33 - Alykhiyil risk - (3E)-n)3 systiss's co-sés iss-syss was live 3: 56 lu) is pi 19. (Jessio 1932) r x, Psy-sys - Sssy). -\ss-ny's to oreyss-sys) 3) Sys) -st; it S-P is 3) syst a yook fog \ is 33 sk y) { res's p R&g, ſtyles is is toys - systs S.

Page 1 Page 2 Page 3 Page 4 Page 5 Page 6 Page 7 Page 8 Page 9 ...
Page 1. Page 2. Page 3. Page 4. Page 5. Page 6. Page 7. Page 8. Page 9. Page 10. Page 11. Page 12. Page 13. Page 14. Page 15. Page 16. Page 17. Page 18 ...

sample paper Aptoinn nata sample model question paper - 1.pdf ...
Page 1 of 1. SAMPLE SET - 1. Note:These questions are the collections of student's contributions from different forums & websites. www.aptoinn.in. 917630 1689 / 98847 22837 / 91765 62187. NATA & J.E.E - B.Arch portions covered. We provide the student

Page 1 Page 2 Page 3
Page 1. Page 2. Page 3.

sample paper Aptoinn nata sample model question paper - 10.pdf ...
PART – B MARKS : 100 TIME: 2 Hours. 1 You are visiting a mall where there are food stalls. you are enjoying. your cold drink in the ground floor of the 5 storey atrium space. seeing the people moving on the escalators , shopping , performing. cultu

Page 1 Page 2
Page 1. Page 2.

sample paper Aptoinn nata sample model question paper - 6.pdf ...
Our Branches: Chennai (Tambaram/Saidapet/Anna nagar/Thiruvanmiyur) / Coimbatore / Salem. NATIONAL APTITUDE TEST IN ARCHITECTURE. PART – B ...

SAMPLE - WT.pdf
Explain!the!working!process!of! DNS server. !!!!!!!!!! [1+1+!4]. 8. Suppose!that!you!are!a Web!Programmer!in!one!of!the!leading Web!Development!Company!in!

Sample PDF
Page 2 of 2. Publish to the web. Publishing your docs online is easy – just select Share > Publish as. web page (in documents) or the Publish tab (in ...

Sample Protagonist Group - WordPress.com
Sample Protagonist Group. The radio towers ... A voice crackled over the airwaves. “This… is… ... call him a shrewd businessman, and a few claim he's a mutant with a set of ... Dystopians daily offer their services to Mr. Fandible in exchange .

Trainee Sample
Using the curriculum vitae for policy research: an evaluation of National Institutes of Health center and training support on career trajectories. Monica Gaughan.

Client: Date: Sample: Notes about sample: ALL ... -
Manually select mult factor from table based on average diameter--->. Color. Bc. Bf. Bf. Bf ... Must balance with fungi for the plant desired. Less than 300 ug, will ...

Aakash Institute Sample Paper
(1) Tidal energy. (2) Temperature difference ..... 86. Find the odd one out. (1) Mango. (2) Banana. (3) Apple. (4) Carrot. 87. Complete the series. ? (1). (2). (3). (4) ...

Determining Sample Size
Nov 6, 1992 - This document is Fact Sheet PEOD-6, a series of the Program Evaluation and Organizational Development, Florida Cooperative Extension.

SAMPLE QUOTE.pdf
6:30 - Brew Bus Picks up at French Quarter hotel. DESCRIPTION. *tax & gratuity included in price. *Brew Bus comfortably holds up to a maximum 15 guests.