Simulation of spatially correlated PolSAR images using inverse transform method

Abstract. This paper proposes an algorithm of simulating spatially correlated polarimetric synthetic aperture radar (PolSAR) images based on the inverse transform method (ITM). Three flexible non-Gaussian models are employed as the underlying distributions of PolSAR images, including the KummerU, W and M models. Additionally, the spatial correlation of the texture component is considered, which is described by a parametric model called the anisotropic Gaussian function. In the algorithm, PolSAR images are simulated by multiplying two independent components, the speckle and texture, that are generated separately. There are two main contributions referring to two important aspects of the ITM. First, the inverse cumulative distribution functions of all the considered texture distributions are mathematically derived, including the Fisher, Beta, and inverse Beta models. Second, considering the high computational complexities the implicitly expressed correlation transfer functions of these texture distributions have, we develop an alternative fast scheme for their computation by using piecewise linear functions. The effectiveness of the proposed simulation algorithm is demonstrated with respect to both the probability density function and spatial correlation.

[1]  S. Quegan,et al.  Understanding Synthetic Aperture Radar Images , 1998 .

[2]  Huanxin Zou,et al.  SAR Image Segmentation via Hierarchical Region Merging and Edge Evolving With Generalized Gamma Distribution , 2014, IEEE Geoscience and Remote Sensing Letters.

[3]  David A. Clausi,et al.  Unsupervised Polarimetric SAR Image Segmentation and Classification Using Region Growing With Edge Penalty , 2012, IEEE Transactions on Geoscience and Remote Sensing.

[4]  Gui Gao,et al.  Statistical Modeling of SAR Images: A Survey , 2010, Sensors.

[5]  Laurent Ferro-Famil,et al.  Non Stationary Texture Analysis from Polarimetric SAR Data , 2004 .

[6]  Jong-Sen Lee,et al.  Polarimetric SAR speckle filtering and its implication for classification , 1999, IEEE Trans. Geosci. Remote. Sens..

[7]  Sidnei J. S. Sant'Anna,et al.  Classification of Segments in PolSAR Imagery by Minimum Stochastic Distances Between Wishart Distributions , 2013, IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing.

[8]  Stian Normann Anfinsen,et al.  Statistical Analysis of Multilook Polarimetric Radar Images with the Mellin Transform , 2010 .

[9]  Lionel Bombrun,et al.  Fisher Distribution for Texture Modeling of Polarimetric SAR Data , 2008, IEEE Geoscience and Remote Sensing Letters.

[10]  Torbjørn Eltoft,et al.  Goodness-of-Fit Tests for Multilook Polarimetric Radar Data Based on the Mellin Transform , 2011, IEEE Transactions on Geoscience and Remote Sensing.

[11]  Alejandro C. Frery,et al.  The polarimetric 𝒢 distribution for SAR data analysis , 2005 .

[12]  Henri Maître,et al.  A new statistical model for Markovian classification of urban areas in high-resolution SAR images , 2004, IEEE Transactions on Geoscience and Remote Sensing.

[13]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[14]  Stian Normann Anfinsen,et al.  A Complete Coverage of Log-Cumulant Space in Terms of Distributions for Polarimetric SAR Data , 2011 .

[15]  Oscar H. Bustos,et al.  Simulation of Spatially Correlated Clutter Fields , 2009, Commun. Stat. Simul. Comput..

[16]  P. Odell,et al.  A Numerical Procedure to Generate a Sample Covariance Matrix , 1966 .

[17]  Armand Lopes,et al.  Optimal speckle reduction for the product model in multilook polarimetric SAR imagery and the Wishart distribution , 1997, IEEE Trans. Geosci. Remote. Sens..

[18]  Joseph A. O'Sullivan,et al.  Quantitative statistical assessment of conditional models for synthetic aperture radar , 2004, IEEE Transactions on Image Processing.

[19]  Gabriel Vasile,et al.  Hierarchical segmentation of Polarimetric SAR images using heterogeneous clutter models , 2009, 2009 IEEE International Geoscience and Remote Sensing Symposium.

[20]  Kevin Barraclough,et al.  I and i , 2001, BMJ : British Medical Journal.

[21]  Ya-Qiu Jin,et al.  Imaging Simulation of Polarimetric SAR for a Comprehensive Terrain Scene Using the Mapping and Projection Algorithm , 2006, IEEE Transactions on Geoscience and Remote Sensing.

[22]  J. Kong,et al.  K-Distribution and Polarimetric Terrain Radar Clutter , 1990, Progress In Electromagnetics Research.

[23]  Torbjørn Eltoft,et al.  Application of the Matrix-Variate Mellin Transform to Analysis of Polarimetric Radar Images , 2011, IEEE Transactions on Geoscience and Remote Sensing.

[24]  Oscar H. Bustos,et al.  Generalized Method for Sampling Spatially Correlated Heterogeneous Speckled Imagery , 2001, EURASIP J. Adv. Signal Process..

[25]  Laurent Ferro-Famil,et al.  Spatially Nonstationary Anisotropic Texture Analysis in SAR Images , 2007, IEEE Transactions on Geoscience and Remote Sensing.

[26]  Jong-Sen Lee,et al.  Intensity and phase statistics of multilook polarimetric and interferometric SAR imagery , 1994, IEEE Trans. Geosci. Remote. Sens..

[27]  Michael J. Collins,et al.  Modeling and Simulation of SAR Image Texture , 2009, IEEE Transactions on Geoscience and Remote Sensing.

[28]  Leslie M. Novak,et al.  Studies of target detection algorithms that use polarimetric radar data , 1988 .

[29]  E. Pottier,et al.  Polarimetric Radar Imaging: From Basics to Applications , 2009 .

[30]  C. D. Vale,et al.  Simulating multivariate nonnormal distributions , 1983 .

[31]  N. R. Goodman Statistical analysis based on a certain multivariate complex Gaussian distribution , 1963 .

[32]  Leslie M. Novak,et al.  Optimal speckle reduction in polarimetric SAR imagery , 1990 .