-
Regression Model for Speckled Data with Extremely Variability
Authors:
A. D. C. Nascimento,
J. M. Vasconcelos,
R. J. Cintra,
A. C. Frery
Abstract:
Synthetic aperture radar (SAR) is an efficient and widely used remote sensing tool. However, data extracted from SAR images are contaminated with speckle, which precludes the application of techniques based on the assumption of additive and normally distributed noise. One of the most successful approaches to describing such data is the multiplicative model, where intensities can follow a variety o…
▽ More
Synthetic aperture radar (SAR) is an efficient and widely used remote sensing tool. However, data extracted from SAR images are contaminated with speckle, which precludes the application of techniques based on the assumption of additive and normally distributed noise. One of the most successful approaches to describing such data is the multiplicative model, where intensities can follow a variety of distributions with positive support. The $\mathcal{G}^0_I$ model is among the most successful ones. Although several estimation methods for the $\mathcal{G}^0_I$ parameters have been proposed, there is no work exploring a regression structure for this model. Such a structure could allow us to infer unobserved values from available ones. In this work, we propose a $\mathcal{G}^0_I$ regression model and use it to describe the influence of intensities from other polarimetric channels. We derive some theoretical properties for the new model: Fisher information matrix, residual measures, and influential tools. Maximum likelihood point and interval estimation methods are proposed and evaluated by Monte Carlo experiments. Results from simulated and actual data show that the new model can be helpful for SAR image analysis.
△ Less
Submitted 14 October, 2024;
originally announced October 2024.
-
Statistical Properties of the Entropy from Ordinal Patterns
Authors:
Eduarda T. C. Chagas,
Alejandro. C. Frery,
Juliana Gambini,
Magdalena M. Lucini,
Heitor S. Ramos,
Andrea A. Rey
Abstract:
The ultimate purpose of the statistical analysis of ordinal patterns is to characterize the distribution of the features they induce. In particular, knowing the joint distribution of the pair Entropy-Statistical Complexity for a large class of time series models would allow statistical tests that are unavailable to date. Working in this direction, we characterize the asymptotic distribution of the…
▽ More
The ultimate purpose of the statistical analysis of ordinal patterns is to characterize the distribution of the features they induce. In particular, knowing the joint distribution of the pair Entropy-Statistical Complexity for a large class of time series models would allow statistical tests that are unavailable to date. Working in this direction, we characterize the asymptotic distribution of the empirical Shannon's Entropy for any model under which the true normalized Entropy is neither zero nor one. We obtain the asymptotic distribution from the Central Limit Theorem (assuming large time series), the Multivariate Delta Method, and a third-order correction of its mean value. We discuss the applicability of other results (exact, first-, and second-order corrections) regarding their accuracy and numerical stability. Within a general framework for building test statistics about Shannon's Entropy, we present a bilateral test that verifies if there is enough evidence to reject the hypothesis that two signals produce ordinal patterns with the same Shannon's Entropy. We applied this bilateral test to the daily maximum temperature time series from three cities (Dublin, Edinburgh, and Miami) and obtained sensible results.
△ Less
Submitted 15 September, 2022;
originally announced September 2022.
-
Dual Approaches to Express the Generalized Degree of Polarimetric Purity
Authors:
Avik Bhattacharya,
Subhadip Dey,
Alejandro C. Frery
Abstract:
The degree of polarimetric purity is an invariant dimensionless quantity that characterizes the closeness of a polarization state of a wave to a pure state and is related to the Von Neumann entropy. The polarimetric purity of a plane wave characterized by the second-order statistics (i.e., the covariance matrix) is uniquely described by the degree of polarization. However, the 2D formalism is only…
▽ More
The degree of polarimetric purity is an invariant dimensionless quantity that characterizes the closeness of a polarization state of a wave to a pure state and is related to the Von Neumann entropy. The polarimetric purity of a plane wave characterized by the second-order statistics (i.e., the covariance matrix) is uniquely described by the degree of polarization. However, the 2D formalism is only applicable when the wave propagation direction is fixed. This assumption is typical in optical and radar polarimetric measurements. Therefore, one must consider all the components to describe the general state of wave polarization. Starting from Samson and Barakat, several different concepts have been proposed in the literature to describe the 3D degree of polarization. We discuss two new ways of achieving such description: by the Coefficient of Variation and by a Direct Sum Decomposition.
△ Less
Submitted 12 February, 2022;
originally announced February 2022.
-
A PolSAR Scattering Power Factorization Framework and Novel Roll-Invariant Parameters Based Unsupervised Classification Scheme Using a Geodesic Distance
Authors:
Debanshu Ratha,
Eric Pottier,
Avik Bhattacharya,
Alejandro C. Frery
Abstract:
We propose a generic Scattering Power Factorization Framework (SPFF) for Polarimetric Synthetic Aperture Radar (PolSAR) data to directly obtain $N$ scattering power components along with a residue power component for each pixel. Each scattering power component is factorized into similarity (or dissimilarity) using elementary targets and a generalized random volume model. The similarity measure is…
▽ More
We propose a generic Scattering Power Factorization Framework (SPFF) for Polarimetric Synthetic Aperture Radar (PolSAR) data to directly obtain $N$ scattering power components along with a residue power component for each pixel. Each scattering power component is factorized into similarity (or dissimilarity) using elementary targets and a generalized random volume model. The similarity measure is derived using a geodesic distance between pairs of $4\times4$ real Kennaugh matrices. In standard model-based decomposition schemes, the $3\times3$ Hermitian positive semi-definite covariance (or coherency) matrix is expressed as a weighted linear combination of scattering targets following a fixed hierarchical process. In contrast, under the proposed framework, a convex splitting of unity is performed to obtain the weights while preserving the dominance of the scattering components. The product of the total power (Span) with these weights provides the non-negative scattering power components. Furthermore, the framework along the geodesic distance is effectively used to obtain specific roll-invariant parameters which are then utilized to design an unsupervised classification scheme. The SPFF, the roll invariant parameters, and the classification results are assessed using C-band RADARSAT-2 and L-band ALOS-2 images of San Francisco.
△ Less
Submitted 27 June, 2019;
originally announced June 2019.
-
Comparing Samples from the $\mathcal{G}^0$ Distribution using a Geodesic Distance
Authors:
Alejandro C. Frery,
Juliana Gambini
Abstract:
The $\mathcal{G}^0$ distribution is widely used for monopolarized SAR image modeling because it can characterize regions with different degree of texture accurately. It is indexed by three parameters: the number of looks (which can be estimated for the whole image), a scale parameter and a texture parameter. This paper presents a new proposal for comparing samples from the $\mathcal{G}^0$ distribu…
▽ More
The $\mathcal{G}^0$ distribution is widely used for monopolarized SAR image modeling because it can characterize regions with different degree of texture accurately. It is indexed by three parameters: the number of looks (which can be estimated for the whole image), a scale parameter and a texture parameter. This paper presents a new proposal for comparing samples from the $\mathcal{G}^0$ distribution using a Geodesic Distance (GD) as a measure of dissimilarity between models. The objective is quantifying the difference between pairs of samples from SAR data using both local parameters (scale and texture) of the $\mathcal{G}^0$ distribution. We propose three tests based on the GD which combine the tests presented in~\cite{GeodesicDistanceGI0JSTARS}, and we estimate their probability distributions using permutation methods.
△ Less
Submitted 23 April, 2019;
originally announced April 2019.
-
Parameter Estimation for the Single-Look $\mathcal{G}^0$ Distribution
Authors:
Débora Chan,
Andrea Rey,
Juliana Gambini,
Alejandro C. Frery
Abstract:
The statistical properties of Synthetic Aperture Radar (SAR) image texture reveals useful target characteristics. It is well-known that these images are affected by speckle, and prone to contamination as double bounce and corner reflectors. The $\mathcal{G}^0$ distribution is flexible enough to model different degrees of texture in speckled data. It is indexed by three parameters: $α$, related to…
▽ More
The statistical properties of Synthetic Aperture Radar (SAR) image texture reveals useful target characteristics. It is well-known that these images are affected by speckle, and prone to contamination as double bounce and corner reflectors. The $\mathcal{G}^0$ distribution is flexible enough to model different degrees of texture in speckled data. It is indexed by three parameters: $α$, related to the texture, $γ$, a scale parameter, and $L$, the number of looks which is related to the signal-to-noise ratio. Quality estimation of $α$ is essential due to its immediate interpretability. In this article, we compare the behavior of a number of parameter estimation techniques in the noisiest case, namely single look data. We evaluate them using Monte Carlo methods for non-contaminated and contaminated data, considering convergence rate, bias, mean squared error (MSE) and computational cost. The results are verified with simulated and actual SAR images.
△ Less
Submitted 29 September, 2018;
originally announced October 2018.
-
Revisiting the effect of spatial resolution on information content based on classification results
Authors:
M. G. Palacio,
S. B. Ferrero,
A. C. Frery
Abstract:
Polarimetric Synthetic Aperture Radar (PolSAR) images are an important source of information. Speckle noise gives SAR images a granular appearance that makes interpretation and analysis hard tasks. A major issue is the assessment of information content in these kind of images, and how it is affected by usual processing techniques. Previous works have resulted in various approaches for quantifying…
▽ More
Polarimetric Synthetic Aperture Radar (PolSAR) images are an important source of information. Speckle noise gives SAR images a granular appearance that makes interpretation and analysis hard tasks. A major issue is the assessment of information content in these kind of images, and how it is affected by usual processing techniques. Previous works have resulted in various approaches for quantifying image information content. As Narayanan, Desetty, and Reichenbach(2002) we study this problem from the classification accuracy viewpoint, focusing in the filtering and the classification stages. Thus, through classified images we verify how changing properties of the input data affects their quality. Our input is an actual PolSAR image, the control parameters are the filter (Local Mean or Model Based PolSAR, MBPolSAR), the size of them and the classification method (Maximum Likelihood, ML, or Support Vector Machine, SVM), and the output are the classification precision obtained applying the classification algorithm to the filtered data. To expand the conclusions, this study deals not only with Classification Accuracy, but also with Kappa and Overall Accuracy as measures of map precision. Experiments were conducted on two airborne PolSAR images. Unless Narayanan, Desetty, and Reichenbach(2002) almost all measure values are good and increase with degradation, i.e. the filtering algorithm that we used always improves the classification results at least up to 7x7.
△ Less
Submitted 13 August, 2018; v1 submitted 26 July, 2018;
originally announced July 2018.
-
Fast Matrix Inversion and Determinant Computation for Polarimetric Synthetic Aperture Radar
Authors:
D. F. G. Coelho,
R. J. Cintra,
A. C. Frery,
V. S. Dimitrov
Abstract:
This paper introduces a fast algorithm for simultaneous inversion and determinant computation of small sized matrices in the context of fully Polarimetric Synthetic Aperture Radar (PolSAR) image processing and analysis. The proposed fast algorithm is based on the computation of the adjoint matrix and the symmetry of the input matrix. The algorithm is implemented in a general purpose graphical proc…
▽ More
This paper introduces a fast algorithm for simultaneous inversion and determinant computation of small sized matrices in the context of fully Polarimetric Synthetic Aperture Radar (PolSAR) image processing and analysis. The proposed fast algorithm is based on the computation of the adjoint matrix and the symmetry of the input matrix. The algorithm is implemented in a general purpose graphical processing unit (GPGPU) and compared to the usual approach based on Cholesky factorization. The assessment with simulated observations and data from an actual PolSAR sensor show a speedup factor of about two when compared to the usual Cholesky factorization. Moreover, the expressions provided here can be implemented in any platform.
△ Less
Submitted 21 July, 2018;
originally announced July 2018.
-
Region-Based Classification of PolSAR Data Using Radial Basis Kernel Functions With Stochastic Distances
Authors:
R. G. Negri,
A. C. Frery,
W. B. Silva,
T. S. G. Mendes,
L. V. Dutra
Abstract:
Region-based classification of PolSAR data can be effectively performed by seeking for the assignment that minimizes a distance between prototypes and segments. Silva et al (2013) used stochastic distances between complex multivariate Wishart models which, differently from other measures, are computationally tractable. In this work we assess the robustness of such approach with respect to errors i…
▽ More
Region-based classification of PolSAR data can be effectively performed by seeking for the assignment that minimizes a distance between prototypes and segments. Silva et al (2013) used stochastic distances between complex multivariate Wishart models which, differently from other measures, are computationally tractable. In this work we assess the robustness of such approach with respect to errors in the training stage, and propose an extension that alleviates such problems. We introduce robustness in the process by incorporating a combination of radial basis kernel functions and stochastic distances with Support Vector Machines (SVM). We consider several stochastic distances between Wishart: Bhatacharyya, Kullback-Leibler, Chi-Square, Rényi, and Hellinger. We perform two case studies with PolSAR images, both simulated and from actual sensors, and different classification scenarios to compare the performance of Minimum Distance and SVM classification frameworks. With this, we model the situation of imperfect training samples. We show that SVM with the proposed kernel functions achieves better performance with respect to Minimum Distance, at the expense of more computational resources and the need of parameter tuning. Code and data are provided for reproducibility.
△ Less
Submitted 7 May, 2018;
originally announced May 2018.
-
Modeling Dengue Vector Population Using Remotely Sensed Data and Machine Learning
Authors:
J. M. Scavuzzo,
F. Trucco,
M. Espinosa,
C. B. Tauro,
M. Abril,
C. M. Scavuzzo,
A. C. Frery
Abstract:
Mosquitoes are vectors of many human diseases. In particular, Aedes ægypti (Linnaeus) is the main vector for Chikungunya, Dengue, and Zika viruses in Latin America and it represents a global threat. Public health policies that aim at combating this vector require dependable and timely information, which is usually expensive to obtain with field campaigns. For this reason, several efforts have been…
▽ More
Mosquitoes are vectors of many human diseases. In particular, Aedes ægypti (Linnaeus) is the main vector for Chikungunya, Dengue, and Zika viruses in Latin America and it represents a global threat. Public health policies that aim at combating this vector require dependable and timely information, which is usually expensive to obtain with field campaigns. For this reason, several efforts have been done to use remote sensing due to its reduced cost. The present work includes the temporal modeling of the oviposition activity (measured weekly on 50 ovitraps in a north Argentinean city) of Aedes ægypti (Linnaeus), based on time series of data extracted from operational earth observation satellite images. We use are NDVI, NDWI, LST night, LST day and TRMM-GPM rain from 2012 to 2016 as predictive variables. In contrast to previous works which use linear models, we employ Machine Learning techniques using completely accessible open source toolkits. These models have the advantages of being non-parametric and capable of describing nonlinear relationships between variables. Specifically, in addition to two linear approaches, we assess a Support Vector Machine, an Artificial Neural Networks, a K-nearest neighbors and a Decision Tree Regressor. Considerations are made on parameter tuning and the validation and training approach. The results are compared to linear models used in previous works with similar data sets for generating temporal predictive models. These new tools perform better than linear approaches, in particular Nearest Neighbor Regression (KNNR) performs the best. These results provide better alternatives to be implemented operatively on the Argentine geospatial Risk system that is running since 2012.
△ Less
Submitted 4 May, 2018;
originally announced May 2018.
-
Detecting Changes in Fully Polarimetric SAR Imagery with Statistical Information Theory
Authors:
Abraão D. C. Nascimento,
Alejandro C. Frery,
Renato J. Cintra
Abstract:
Images obtained from coherent illumination processes are contaminated with speckle. A prominent example of such imagery systems is the polarimetric synthetic aperture radar (PolSAR). For such remote sensing tool the speckle interference pattern appears in the form of a positive definite Hermitian matrix, which requires specialized models and makes change detection a hard task. The scaled complex W…
▽ More
Images obtained from coherent illumination processes are contaminated with speckle. A prominent example of such imagery systems is the polarimetric synthetic aperture radar (PolSAR). For such remote sensing tool the speckle interference pattern appears in the form of a positive definite Hermitian matrix, which requires specialized models and makes change detection a hard task. The scaled complex Wishart distribution is a widely used model for PolSAR images. Such distribution is defined by two parameters: the number of looks and the complex covariance matrix. The last parameter contains all the necessary information to characterize the backscattered data and, thus, identifying changes in a sequence of images can be formulated as a problem of verifying whether the complex covariance matrices differ at two or more takes. This paper proposes a comparison between a classical change detection method based on the likelihood ratio and three statistical methods that depend on information-theoretic measures: the Kullback-Leibler distance and two entropies. The performance of these four tests was quantified in terms of their sample test powers and sizes using simulated data. The tests are then applied to actual PolSAR data. The results provide evidence that tests based on entropies may outperform those based on the Kullback-Leibler distance and likelihood ratio statistics.
△ Less
Submitted 16 September, 2018; v1 submitted 26 January, 2018;
originally announced January 2018.
-
Unsupervised Classification of PolSAR Data Using a Scattering Similarity Measure Derived from a Geodesic Distance
Authors:
Debanshu Ratha,
Avik Bhattacharya,
Alejandro C. Frery
Abstract:
In this letter, we propose a novel technique for obtaining scattering components from Polarimetric Synthetic Aperture Radar (PolSAR) data using the geodesic distance on the unit sphere. This geodesic distance is obtained between an elementary target and the observed Kennaugh matrix, and it is further utilized to compute a similarity measure between scattering mechanisms. The normalized similarity…
▽ More
In this letter, we propose a novel technique for obtaining scattering components from Polarimetric Synthetic Aperture Radar (PolSAR) data using the geodesic distance on the unit sphere. This geodesic distance is obtained between an elementary target and the observed Kennaugh matrix, and it is further utilized to compute a similarity measure between scattering mechanisms. The normalized similarity measure for each elementary target is then modulated with the total scattering power (Span). This measure is used to categorize pixels into three categories i.e. odd-bounce, double-bounce and volume, depending on which of the above scattering mechanisms dominate. Then the maximum likelihood classifier of [J.-S. Lee, M. R. Grunes, E. Pottier, and L. Ferro-Famil, Unsupervised terrain classification preserving polarimetric scattering characteristics, IEEE Trans. Geos. Rem. Sens., vol. 42, no. 4, pp. 722731, April 2004.] based on the complex Wishart distribution is iteratively used for each category. Dominant scattering mechanisms are thus preserved in this classification scheme. We show results for L-band AIRSAR and ALOS-2 datasets acquired over San Francisco and Mumbai, respectively. The scattering mechanisms are better preserved using the proposed methodology than the unsupervised classification results using the Freeman-Durden scattering powers on an orientation angle (OA) corrected PolSAR image. Furthermore, (1) the scattering similarity is a completely non-negative quantity unlike the negative powers that might occur in double- bounce and odd-bounce scattering component under Freeman Durden decomposition (FDD), and (2) the methodology can be extended to more canonical targets as well as for bistatic scattering.
△ Less
Submitted 1 December, 2017;
originally announced December 2017.
-
Unassisted Quantitative Evaluation Of Despeckling Filters
Authors:
Luis Gomez,
Raydonal Ospina,
Alejandro C. Frery
Abstract:
SAR (Synthetic Aperture Radar) imaging plays a central role in Remote Sensing due to, among other important features, its ability to provide high-resolution, day-and-night and almost weather-independent images. SAR images are affected from a granular contamination, speckle, that can be described by a multiplicative model. Many despeckling techniques have been proposed in the literature, as well as…
▽ More
SAR (Synthetic Aperture Radar) imaging plays a central role in Remote Sensing due to, among other important features, its ability to provide high-resolution, day-and-night and almost weather-independent images. SAR images are affected from a granular contamination, speckle, that can be described by a multiplicative model. Many despeckling techniques have been proposed in the literature, as well as measures of the quality of the results they provide. Assuming the multiplicative model, the observed image $Z$ is the product of two independent fields: the backscatter $X$ and the speckle $Y$. The result of any speckle filter is $\widehat X$, an estimator of the backscatter $X$, based solely on the observed data $Z$. An ideal estimator would be the one for which the ratio of the observed image to the filtered one $I=Z/\widehat X$ is only speckle: a collection of independent identically distributed samples from Gamma variates. We, then, assess the quality of a filter by the closeness of $I$ to the hypothesis that it is adherent to the statistical properties of pure speckle. We analyze filters through the ratio image they produce with regards to first- and second-order statistics: the former check marginal properties, while the latter verifies lack of structure. A new quantitative image-quality index is then defined, and applied to state-of-the-art despeckling filters. This new measure provides consistent results with commonly used quality measures (equivalent number of looks, PSNR, MSSIM, $β$ edge correlation, and preservation of the mean), and ranks the filters results also in agreement with their visual analysis. We conclude our study showing that the proposed measure can be successfully used to optimize the (often many) parameters that define a speckle filter.
△ Less
Submitted 19 April, 2017;
originally announced April 2017.
-
The Geodesic Distance between $\mathcal{G}_I^0$ Models and its Application to Region Discrimination
Authors:
José Naranjo-Torres,
Juliana Gambini,
Alejandro C. Frery
Abstract:
The $\mathcal{G}_I^0$ distribution is able to characterize different regions in monopolarized SAR imagery. It is indexed by three parameters: the number of looks (which can be estimated in the whole image), a scale parameter and a texture parameter. This paper presents a new proposal for feature extraction and region discrimination in SAR imagery, using the geodesic distance as a measure of dissim…
▽ More
The $\mathcal{G}_I^0$ distribution is able to characterize different regions in monopolarized SAR imagery. It is indexed by three parameters: the number of looks (which can be estimated in the whole image), a scale parameter and a texture parameter. This paper presents a new proposal for feature extraction and region discrimination in SAR imagery, using the geodesic distance as a measure of dissimilarity between $\mathcal{G}_I^0$ models. We derive geodesic distances between models that describe several practical situations, assuming the number of looks known, for same and different texture and for same and different scale. We then apply this new tool to the problems of (i)~identifying edges between regions with different texture, and (ii)~quantify the dissimilarity between pairs of samples in actual SAR data. We analyze the advantages of using the geodesic distance when compared to stochastic distances.
△ Less
Submitted 1 January, 2017;
originally announced January 2017.
-
Classification and Verification of Online Handwritten Signatures with Time Causal Information Theory Quantifiers
Authors:
Osvaldo A. Rosso,
Raydonal Ospina,
Alejandro C. Frery
Abstract:
We present a new approach for online handwritten signature classification and verification based on descriptors stemming from Information Theory. The proposal uses the Shannon Entropy, the Statistical Complexity, and the Fisher Information evaluated over the Bandt and Pompe symbolization of the horizontal and vertical coordinates of signatures. These six features are easy and fast to compute, and…
▽ More
We present a new approach for online handwritten signature classification and verification based on descriptors stemming from Information Theory. The proposal uses the Shannon Entropy, the Statistical Complexity, and the Fisher Information evaluated over the Bandt and Pompe symbolization of the horizontal and vertical coordinates of signatures. These six features are easy and fast to compute, and they are the input to an One-Class Support Vector Machine classifier. The results produced surpass state-of-the-art techniques that employ higher-dimensional feature spaces which often require specialized software and hardware. We assess the consistency of our proposal with respect to the size of the training sample, and we also use it to classify the signatures into meaningful groups.
△ Less
Submitted 26 January, 2016;
originally announced January 2016.
-
Characterization of Vehicle Behavior with Information Theory
Authors:
Andre L. L. Aquino,
Tamer S. G. Cavalcante,
Eliana S. Almeida,
Alejandro C. Frery,
Osvaldo A. Rosso
Abstract:
This work proposes the use of Information Theory for the characterization of vehicles behavior through their velocities. Three public data sets were used: i.Mobile Century data set collected on Highway I-880, near Union City, California; ii.Borlänge GPS data set collected in the Swedish city of Borlänge; and iii.Beijing taxicabs data set collected in Beijing, China, where each vehicle speed is sto…
▽ More
This work proposes the use of Information Theory for the characterization of vehicles behavior through their velocities. Three public data sets were used: i.Mobile Century data set collected on Highway I-880, near Union City, California; ii.Borlänge GPS data set collected in the Swedish city of Borlänge; and iii.Beijing taxicabs data set collected in Beijing, China, where each vehicle speed is stored as a time series. The Bandt-Pompe methodology combined with the Complexity-Entropy plane were used to identify different regimes and behaviors. The global velocity is compatible with a correlated noise with f^{-k} Power Spectrum with k >= 0. With this we identify traffic behaviors as, for instance, random velocities (k aprox. 0) when there is congestion, and more correlated velocities (k aprox. 3) in the presence of free traffic flow.
△ Less
Submitted 14 October, 2015;
originally announced October 2015.
-
Comments on "Detecting Outliers in Gamma Distribution" by M. Jabbari Nooghabi et al. (2010)
Authors:
M. Magdalena Lucini,
Alejandro C. Frery
Abstract:
This note shows that the results presented by Jabbari Nooghabi et al. (2010) do not hold in all expected cases. With this, the technique proposed by Kumar and Lalhita (2012) for detecting upper outliers in Gamma samples is also not valid. Specifically, this note shows that the probability density functions (pdf) under the null hypothesis of the test statistics therein proposed are not always valid…
▽ More
This note shows that the results presented by Jabbari Nooghabi et al. (2010) do not hold in all expected cases. With this, the technique proposed by Kumar and Lalhita (2012) for detecting upper outliers in Gamma samples is also not valid. Specifically, this note shows that the probability density functions (pdf) under the null hypothesis of the test statistics therein proposed are not always valid.
△ Less
Submitted 16 September, 2015;
originally announced September 2015.
-
Optical images-based edge detection in Synthetic Aperture Radar images
Authors:
Gilberto P. Silva Junior,
Alejandro C. Frery,
Sandra Sandri,
Humberto Bustince,
Edurne Barrenechea,
Cédric Marco-Detchart
Abstract:
We address the issue of adapting optical images-based edge detection techniques for use in Polarimetric Synthetic Aperture Radar (PolSAR) imagery. We modify the gravitational edge detection technique (inspired by the Law of Universal Gravity) proposed by Lopez-Molina et al, using the non-standard neighbourhood configuration proposed by Fu et al, to reduce the speckle noise in polarimetric SAR imag…
▽ More
We address the issue of adapting optical images-based edge detection techniques for use in Polarimetric Synthetic Aperture Radar (PolSAR) imagery. We modify the gravitational edge detection technique (inspired by the Law of Universal Gravity) proposed by Lopez-Molina et al, using the non-standard neighbourhood configuration proposed by Fu et al, to reduce the speckle noise in polarimetric SAR imagery. We compare the modified and unmodified versions of the gravitational edge detection technique with the well-established one proposed by Canny, as well as with a recent multiscale fuzzy-based technique proposed by Lopez-Molina et Alejandro We also address the issues of aggregation of gray level images before and after edge detection and of filtering. All techniques addressed here are applied to a mosaic built using class distributions obtained from a real scene, as well as to the true PolSAR image; the mosaic results are assessed using Baddeley's Delta Metric. Our experiments show that modifying the gravitational edge detection technique with a non-standard neighbourhood configuration produces better results than the original technique, as well as the other techniques used for comparison. The experiments show that adapting edge detection methods from Computational Intelligence for use in PolSAR imagery is a new field worthy of exploration.
△ Less
Submitted 24 August, 2015;
originally announced August 2015.
-
Classification of Complex Wishart Matrices with a Diffusion-Reaction System guided by Stochastic Distances
Authors:
Luis Gomez,
Luis Alvarez,
Luis Mazorra,
Alejandro C. Frery
Abstract:
We propose a new method for PolSAR (Polarimetric Synthetic Aperture Radar) imagery classification based on stochastic distances in the space of random matrices obeying complex Wishart distributions. Given a collection of prototypes $\{Z_m\}_{m=1}^M$ and a stochastic distance $d(.,.)$, we classify any random matrix $X$ using two criteria in an iterative setup. Firstly, we associate $X$ to the class…
▽ More
We propose a new method for PolSAR (Polarimetric Synthetic Aperture Radar) imagery classification based on stochastic distances in the space of random matrices obeying complex Wishart distributions. Given a collection of prototypes $\{Z_m\}_{m=1}^M$ and a stochastic distance $d(.,.)$, we classify any random matrix $X$ using two criteria in an iterative setup. Firstly, we associate $X$ to the class which minimizes the weighted stochastic distance $w_md(X,Z_m)$, where the positive weights $w_m$ are computed to maximize the class discrimination power. Secondly, we improve the result by embedding the classification problem into a diffusion-reaction partial differential system where the diffusion term smooths the patches within the image, and the reaction term tends to move the pixel values towards the closest class prototype. In particular, the method inherits the benefits of speckle reduction by diffusion-like methods. Results on synthetic and real PolSAR data show the performance of the method.
△ Less
Submitted 17 July, 2015;
originally announced July 2015.
-
Modifying the Yamaguchi Four-Component Decomposition Scattering Powers Using a Stochastic Distance
Authors:
Avik Bhattacharya,
Arnab Muhuri,
Shaunak De,
Surendar Manickam,
Alejandro C. Frery
Abstract:
Model-based decompositions have gained considerable attention after the initial work of Freeman and Durden. This decomposition which assumes the target to be reflection symmetric was later relaxed in the Yamaguchi et al. decomposition with the addition of the helix parameter. Since then many decomposition have been proposed where either the scattering model was modified to fit the data or the cohe…
▽ More
Model-based decompositions have gained considerable attention after the initial work of Freeman and Durden. This decomposition which assumes the target to be reflection symmetric was later relaxed in the Yamaguchi et al. decomposition with the addition of the helix parameter. Since then many decomposition have been proposed where either the scattering model was modified to fit the data or the coherency matrix representing the second order statistics of the full polarimetric data is rotated to fit the scattering model. In this paper we propose to modify the Yamaguchi four-component decomposition (Y4O) scattering powers using the concept of statistical information theory for matrices. In order to achieve this modification we propose a method to estimate the polarization orientation angle (OA) from full-polarimetric SAR images using the Hellinger distance. In this method, the OA is estimated by maximizing the Hellinger distance between the un-rotated and the rotated $T_{33}$ and the $T_{22}$ components of the coherency matrix $\mathbf{[T]}$. Then, the powers of the Yamaguchi four-component model-based decomposition (Y4O) are modified using the maximum relative stochastic distance between the $T_{33}$ and the $T_{22}$ components of the coherency matrix at the estimated OA. The results show that the overall double-bounce powers over rotated urban areas have significantly improved with the reduction of volume powers. The percentage of pixels with negative powers have also decreased from the Y4O decomposition. The proposed method is both qualitatively and quantitatively compared with the results obtained from the Y4O and the Y4R decompositions for a Radarsat-2 C-band San-Francisco dataset and an UAVSAR L-band Hayward dataset.
△ Less
Submitted 4 April, 2015;
originally announced April 2015.
-
Information Theory Perspective on Network Robustness
Authors:
Tiago A. Schieber,
Laura Carpi,
Alejandro C. Frery,
Osvaldo A Rosso,
Panos M. Pardalos,
Martin G. Ravetti
Abstract:
A crucial challenge in network theory is the study of the robustness of a network after facing a sequence of failures. In this work, we propose a dynamical definition of network's robustness based on Information Theory, that considers measurements of the structural changes caused by failures of the network's components. Failures are defined here, as a temporal process defined in a sequence. The ro…
▽ More
A crucial challenge in network theory is the study of the robustness of a network after facing a sequence of failures. In this work, we propose a dynamical definition of network's robustness based on Information Theory, that considers measurements of the structural changes caused by failures of the network's components. Failures are defined here, as a temporal process defined in a sequence. The robustness of the network is then evaluated by measuring dissimilarities between topologies after each time step of the sequence, providing a dynamical information about the topological damage. We thoroughly analyze the efficiency of the method in capturing small perturbations by considering both, the degree and distance distributions. We found the network's distance distribution more consistent in capturing network structural deviations, as better reflects the consequences of the failures. Theoretical examples and real networks are used to study the performance of this methodology.
△ Less
Submitted 30 October, 2015; v1 submitted 30 October, 2014;
originally announced October 2014.
-
MuSA: Multivariate Sampling Algorithm for Wireless Sensor Networks
Authors:
André L. L. Aquino,
Orlando S. Junior,
Alejandro C. Frery,
Édler Lins de Albuquerque,
Raquel A. F. Mini
Abstract:
A wireless sensor network can be used to collect and process environmental data, which is often of multivariate nature. This work proposes a multivariate sampling algorithm based on component analysis techniques in wireless sensor networks. To improve the sampling, the algorithm uses component analysis techniques to rank the data. Once ranked, the most representative data is retained. Simulation r…
▽ More
A wireless sensor network can be used to collect and process environmental data, which is often of multivariate nature. This work proposes a multivariate sampling algorithm based on component analysis techniques in wireless sensor networks. To improve the sampling, the algorithm uses component analysis techniques to rank the data. Once ranked, the most representative data is retained. Simulation results show that our technique reduces the data keeping its representativeness. In addition, the energy consumption and delay to deliver the data on the network are reduced.
△ Less
Submitted 3 August, 2014;
originally announced August 2014.
-
Parameter Estimation in SAR Imagery using Stochastic Distances and Asymmetric Kernels
Authors:
Juliana Gambini,
Julia Cassetti,
María Magdalena Lucini,
Alejandro C. Frery
Abstract:
In this paper we analyze several strategies for the estimation of the roughness parameter of the $\mathcal G_I^0$ distribution. It has been shown that this distribution is able to characterize a large number of targets in monopolarized SAR imagery, deserving the denomination of "Universal Model" It is indexed by three parameters: the number of looks (which can be estimated in the whole image), a s…
▽ More
In this paper we analyze several strategies for the estimation of the roughness parameter of the $\mathcal G_I^0$ distribution. It has been shown that this distribution is able to characterize a large number of targets in monopolarized SAR imagery, deserving the denomination of "Universal Model" It is indexed by three parameters: the number of looks (which can be estimated in the whole image), a scale parameter, and the roughness or texture parameter. The latter is closely related to the number of elementary backscatters in each pixel, one of the reasons for receiving attention in the literature. Although there are efforts in providing improved and robust estimates for such quantity, its dependable estimation still poses numerical problems in practice. We discuss estimators based on the minimization of stochastic distances between empirical and theoretical densities, and argue in favor of using an estimator based on the Triangular distance and asymmetric kernels built with Inverse Gaussian densities. We also provide new results regarding the heavytailedness of this distribution.
△ Less
Submitted 4 August, 2014; v1 submitted 1 August, 2014;
originally announced August 2014.
-
Bias Correction and Modified Profile Likelihood under the Wishart Complex Distribution
Authors:
Abraão D. C. Nascimento,
Alejandro C. Frery,
Renato J. Cintra
Abstract:
This paper proposes improved methods for the maximum likelihood (ML) estimation of the equivalent number of looks $L$. This parameter has a meaningful interpretation in the context of polarimetric synthetic aperture radar (PolSAR) images. Due to the presence of coherent illumination in their processing, PolSAR systems generate images which present a granular noise called speckle. As a potential so…
▽ More
This paper proposes improved methods for the maximum likelihood (ML) estimation of the equivalent number of looks $L$. This parameter has a meaningful interpretation in the context of polarimetric synthetic aperture radar (PolSAR) images. Due to the presence of coherent illumination in their processing, PolSAR systems generate images which present a granular noise called speckle. As a potential solution for reducing such interference, the parameter $L$ controls the signal-noise ratio. Thus, the proposal of efficient estimation methodologies for $L$ has been sought. To that end, we consider firstly that a PolSAR image is well described by the scaled complex Wishart distribution. In recent years, Anfinsen et al. derived and analyzed estimation methods based on the ML and on trace statistical moments for obtaining the parameter $L$ of the unscaled version of such probability law. This paper generalizes that approach. We present the second-order bias expression proposed by Cox and Snell for the ML estimator of this parameter. Moreover, the formula of the profile likelihood modified by Barndorff-Nielsen in terms of $L$ is discussed. Such derivations yield two new ML estimators for the parameter $L$, which are compared to the estimators proposed by Anfinsen et al. The performance of these estimators is assessed by means of Monte Carlo experiments, adopting three statistical measures as comparison criterion: the mean square error, the bias, and the coefficient of variation. Equivalently to the simulation study, an application to actual PolSAR data concludes that the proposed estimators outperform all the others in homogeneous scenarios.
△ Less
Submitted 18 April, 2014;
originally announced April 2014.
-
Visualization of Skewed Data: A Tool in R
Authors:
R. Ospina,
A. M. Larangeiras,
A. C. Frery
Abstract:
In this work we present a visualization tool specifically tailored to deal with skewed data. The technique is based upon the use of two types of notched boxplots (the usual one, and one which is tuned for the skewness of the data), the violin plot, the histogram and a nonparametric estimate of the density. The data is assumed to lie on the same line, so the plots are compatible. We show that a goo…
▽ More
In this work we present a visualization tool specifically tailored to deal with skewed data. The technique is based upon the use of two types of notched boxplots (the usual one, and one which is tuned for the skewness of the data), the violin plot, the histogram and a nonparametric estimate of the density. The data is assumed to lie on the same line, so the plots are compatible. We show that a good deal of information can be extracted from the inspection of this tool; in particular, we apply the technique to analyze data from synthetic aperture radar images. We provide the implementation in R.
△ Less
Submitted 3 March, 2014;
originally announced March 2014.
-
Information Theory and Image Understanding: An Application to Polarimetric SAR Imagery
Authors:
A. C. Frery,
A. D. C. Nascimento,
R. J. Cintra
Abstract:
This work presents a comprehensive examination of the use of information theory for understanding Polarimetric Synthetic Aperture Radar (PolSAR) images by means of contrast measures that can be used as test statistics. Due to the phenomenon called `speckle', common to all images obtained with coherent illumination such as PolSAR imagery, accurate modelling is required in their processing and analy…
▽ More
This work presents a comprehensive examination of the use of information theory for understanding Polarimetric Synthetic Aperture Radar (PolSAR) images by means of contrast measures that can be used as test statistics. Due to the phenomenon called `speckle', common to all images obtained with coherent illumination such as PolSAR imagery, accurate modelling is required in their processing and analysis. The scaled multilook complex Wishart distribution has proven to be a successful approach for modelling radar backscatter from forest and pasture areas. Classification, segmentation, and image analysis techniques which depend on this model have been devised, and many of them employ some kind of dissimilarity measure. Specifically, we introduce statistical tests for analyzing contrast in such images. These tests are based on the chi-square, Kullback-Leibler, Rényi, Bhattacharyya, and Hellinger distances. Results obtained by Monte Carlo experiments reveal the Kullback-Leibler distance as the best one with respect to the empirical test sizes under several situations which include pure and contaminated data. The proposed methodology was applied to actual data, obtained by an E-SAR sensor over surroundings of We$β$ssling, Bavaria, Germany.
△ Less
Submitted 8 February, 2014;
originally announced February 2014.
-
Contrast Measures based on the Complex Correlation Coefficient for PolSAR Imagery
Authors:
Alejandro C. Frery,
Renato J. Cintra,
Abraão D. C. Nascimento
Abstract:
We derive contrast measures which involve the number of looks and the complex correlation coefficient between polarization channels in PolSAR imagery. Using asymptotic results which characterize the behavior of these measures, we derive statistical regions of confidence which lead to test of hypothesis. An application to real data is performed, confirming the importance of the proposals.
We derive contrast measures which involve the number of looks and the complex correlation coefficient between polarization channels in PolSAR imagery. Using asymptotic results which characterize the behavior of these measures, we derive statistical regions of confidence which lead to test of hypothesis. An application to real data is performed, confirming the importance of the proposals.
△ Less
Submitted 8 February, 2014;
originally announced February 2014.
-
The Generalized Statistical Complexity of PolSAR Data
Authors:
Alejandro C. Frery,
Eliana S. de Almeida,
Osvaldo A. Rosso
Abstract:
This paper presents and discusses the use of a new feature for PolSAR imagery: the Generalized Statistical Complexity. This measure is able to capture the disorder of the data by means of the entropy, as well as its departure from a reference distribution. The latter component is obtained by measuring a stochastic distance between two models: the $\mathcal G^0$ and the Gamma laws. Preliminary resu…
▽ More
This paper presents and discusses the use of a new feature for PolSAR imagery: the Generalized Statistical Complexity. This measure is able to capture the disorder of the data by means of the entropy, as well as its departure from a reference distribution. The latter component is obtained by measuring a stochastic distance between two models: the $\mathcal G^0$ and the Gamma laws. Preliminary results on the intensity components of AIRSAR image of San Francisco are encouraging.
△ Less
Submitted 8 February, 2014;
originally announced February 2014.
-
When Data do not Bring Information: A Case Study in Markov Random Fields Estimation
Authors:
J. Gimenez,
A. C. Frery,
Ana Georgina Flesia
Abstract:
The Potts model is frequently used to describe the behavior of image classes, since it allows to incorporate contextual information linking neighboring pixels in a simple way. Its isotropic version has only one real parameter beta, known as smoothness parameter or inverse temperature, which regulates the classes map homogeneity. The classes are unavailable, and estimating them is central in import…
▽ More
The Potts model is frequently used to describe the behavior of image classes, since it allows to incorporate contextual information linking neighboring pixels in a simple way. Its isotropic version has only one real parameter beta, known as smoothness parameter or inverse temperature, which regulates the classes map homogeneity. The classes are unavailable, and estimating them is central in important image processing procedures as, for instance, image classification. Methods for estimating the classes which stem from a Bayesian approach under the Potts model require to adequately specify a value for beta. The estimation of such parameter can be efficiently made solving the Pseudo Maximum likelihood (PML) equations in two different schemes, using the prior or the posterior model. Having only radiometric data available, the first scheme needs the computation of an initial segmentation, while the second uses both the segmentation and the radiometric data to make the estimation. In this paper, we compare these two PML estimators by computing the mean square error (MSE), bias, and sensitivity to deviations from the hypothesis of the model. We conclude that the use of extra data does not improve the accuracy of the PML, moreover, under gross deviations from the model, this extra information introduces unpredictable distortions and bias.
△ Less
Submitted 7 May, 2014; v1 submitted 7 February, 2014;
originally announced February 2014.
-
Distinguishing noise from chaos: objective versus subjective criteria using Horizontal Visibility Graph
Authors:
Martín Gómez Ravetti,
Laura C. Carpi,
Bruna Amin Gonçalves,
Alejandro C. Frery,
Osvaldo A. Rosso
Abstract:
A recently proposed methodology called the Horizontal Visibility Graph (HVG) [Luque {\it et al.}, Phys. Rev. E., 80, 046103 (2009)] that constitutes a geometrical simplification of the well known Visibility Graph algorithm [Lacasa {\it et al.\/}, Proc. Natl. Sci. U.S.A. 105, 4972 (2008)], has been used to study the distinction between deterministic and stochastic components in time series [L. Laca…
▽ More
A recently proposed methodology called the Horizontal Visibility Graph (HVG) [Luque {\it et al.}, Phys. Rev. E., 80, 046103 (2009)] that constitutes a geometrical simplification of the well known Visibility Graph algorithm [Lacasa {\it et al.\/}, Proc. Natl. Sci. U.S.A. 105, 4972 (2008)], has been used to study the distinction between deterministic and stochastic components in time series [L. Lacasa and R. Toral, Phys. Rev. E., 82, 036120 (2010)]. Specifically, the authors propose that the node degree distribution of these processes follows an exponential functional of the form $P(κ)\sim \exp(-λ~κ)$, in which $κ$ is the node degree and $λ$ is a positive parameter able to distinguish between deterministic (chaotic) and stochastic (uncorrelated and correlated) dynamics. In this work, we investigate the characteristics of the node degree distributions constructed by using HVG, for time series corresponding to $28$ chaotic maps and $3$ different stochastic processes. We thoroughly study the methodology proposed by Lacasa and Toral finding several cases for which their hypothesis is not valid. We propose a methodology that uses the HVG together with Information Theory quantifiers. An extensive and careful analysis of the node degree distributions obtained by applying HVG allow us to conclude that the Fisher-Shannon information plane is a remarkable tool able to graphically represent the different nature, deterministic or stochastic, of the systems under study.
△ Less
Submitted 9 January, 2014;
originally announced January 2014.
-
A New Algorithm of Speckle Filtering using Stochastic Distances
Authors:
Leonardo Torres,
Tamer Cavalcante,
Alejandro C. Frery
Abstract:
This paper presents a new approach for filter design based on stochastic distances and tests between distributions. A window is defined around each pixel, overlapping samples are compared and only those which pass a goodness-of-fit test are used to compute the filtered value. The technique is applied to intensity SAR data with homogeneous regions using the Gamma model. The proposal is compared wit…
▽ More
This paper presents a new approach for filter design based on stochastic distances and tests between distributions. A window is defined around each pixel, overlapping samples are compared and only those which pass a goodness-of-fit test are used to compute the filtered value. The technique is applied to intensity SAR data with homogeneous regions using the Gamma model. The proposal is compared with the Lee's filter using a protocol based on Monte Carlo. Among the criteria used to quantify the quality of filters, we employ the equivalent number of looks, line and edge preservation. Moreover, we also assessed the filters by the Universal Image Quality Index and the Pearson's correlation on edges regions.
△ Less
Submitted 29 August, 2013;
originally announced August 2013.
-
SAR Image Despeckling Algorithms using Stochastic Distances and Nonlocal Means
Authors:
Leonardo Torres,
Alejandro C. Frery
Abstract:
This paper presents two approaches for filter design based on stochastic distances for intensity speckle reduction. A window is defined around each pixel, overlapping samples are compared and only those which pass a goodness-of-fit test are used to compute the filtered value. The tests stem from stochastic divergences within the Information Theory framework. The technique is applied to intensity S…
▽ More
This paper presents two approaches for filter design based on stochastic distances for intensity speckle reduction. A window is defined around each pixel, overlapping samples are compared and only those which pass a goodness-of-fit test are used to compute the filtered value. The tests stem from stochastic divergences within the Information Theory framework. The technique is applied to intensity Synthetic Aperture Radar (SAR) data with homogeneous regions using the Gamma model. The first approach uses a Nagao-Matsuyama-type procedure for setting the overlapping samples, and the second uses the nonlocal method. The proposals are compared with the Improved Sigma filter and with anisotropic diffusion for speckled data (SRAD) using a protocol based on Monte Carlo simulation. Among the criteria used to quantify the quality of filters, we employ the equivalent number of looks, and line and edge preservation. Moreover, we also assessed the filters by the Universal Image Quality Index and by the Pearson correlation between edges. Applications to real images are also discussed. The proposed methods show good results.
△ Less
Submitted 20 August, 2013;
originally announced August 2013.
-
Structural Changes in Data Communication in Wireless Sensor Networks
Authors:
Raquel S. Cabral,
Andre L. L. Aquino,
Alejandro C. Frery,
Osvaldo A. Rosso,
Jaime A. Ramírez
Abstract:
Wireless sensor networks are an important technology for making distributed autonomous measures in hostile or inaccessible environments. Among the challenges they pose, the way data travel among them is a relevant issue since their structure is quite dynamic. The operational topology of such devices can often be described by complex networks. In this work, we assess the variation of measures commo…
▽ More
Wireless sensor networks are an important technology for making distributed autonomous measures in hostile or inaccessible environments. Among the challenges they pose, the way data travel among them is a relevant issue since their structure is quite dynamic. The operational topology of such devices can often be described by complex networks. In this work, we assess the variation of measures commonly employed in the complex networks literature applied to wireless sensor networks. Four data communication strategies were considered: geometric, random, small-world, and scale-free models, along with the shortest path length measure. The sensitivity of this measure was analyzed with respect to the following perturbations: insertion and removal of nodes in the geometric strategy; and insertion, removal and rewiring of links in the other models. The assessment was performed using the normalized Kullback-Leibler divergence and Hellinger distance quantifiers, both deriving from the Information Theory framework. The results reveal that the shortest path length is sensitive to perturbations.
△ Less
Submitted 13 August, 2013;
originally announced August 2013.
-
Comparing Edge Detection Methods based on Stochastic Entropies and Distances for PolSAR Imagery
Authors:
Abraão D. C. Nascimento,
Michelle M. Horta,
Alejandro C. Frery,
Renato J. Cintra
Abstract:
Polarimetric synthetic aperture radar (PolSAR) has achieved a prominent position as a remote imaging method. However, PolSAR images are contaminated by speckle noise due to the coherent illumination employed during the data acquisition. This noise provides a granular aspect to the image, making its processing and analysis (such as in edge detection) hard tasks. This paper discusses seven methods f…
▽ More
Polarimetric synthetic aperture radar (PolSAR) has achieved a prominent position as a remote imaging method. However, PolSAR images are contaminated by speckle noise due to the coherent illumination employed during the data acquisition. This noise provides a granular aspect to the image, making its processing and analysis (such as in edge detection) hard tasks. This paper discusses seven methods for edge detection in multilook PolSAR images. In all methods, the basic idea consists in detecting transition points in the finest possible strip of data which spans two regions. The edge is contoured using the transitions points and a B-spline curve. Four stochastic distances, two differences of entropies, and the maximum likelihood criterion were used under the scaled complex Wishart distribution; the first six stem from the h-phi class of measures. The performance of the discussed detection methods was quantified and analyzed by the computational time and probability of correct edge detection, with respect to the number of looks, the backscatter matrix as a whole, the SPAN, the covariance an the spatial resolution. The detection procedures were applied to three real PolSAR images. Results provide evidence that the methods based on the Bhattacharyya distance and the difference of Shannon entropies outperform the other techniques.
△ Less
Submitted 9 June, 2013;
originally announced June 2013.
-
Speckle Reduction with Adaptive Stack Filters
Authors:
María Elena Buemi,
Alejandro C. Frery,
Heitor S. Ramos
Abstract:
Stack filters are a special case of non-linear filters. They have a good performance for filtering images with different types of noise while preserving edges and details. A stack filter decomposes an input image into stacks of binary images according to a set of thresholds. Each binary image is then filtered by a Boolean function, which characterizes the filter. Adaptive stack filters can be comp…
▽ More
Stack filters are a special case of non-linear filters. They have a good performance for filtering images with different types of noise while preserving edges and details. A stack filter decomposes an input image into stacks of binary images according to a set of thresholds. Each binary image is then filtered by a Boolean function, which characterizes the filter. Adaptive stack filters can be computed by training using a prototype (ideal) image and its corrupted version, leading to optimized filters with respect to a loss function. In this work we propose the use of training with selected samples for the estimation of the optimal Boolean function. We study the performance of adaptive stack filters when they are applied to speckled imagery, in particular to Synthetic Aperture Radar (SAR) images. This is done by evaluating the quality of the filtered images through the use of suitable image quality indexes and by measuring the classification accuracy of the resulting images. We used SAR images as input, since they are affected by speckle noise that makes classification a difficult task.
△ Less
Submitted 8 June, 2013;
originally announced June 2013.
-
An adaptive time integration strategy based on displacement history curvature
Authors:
E. N. Lages,
E. S. S. Silveira,
D. T. Cintra,
A. C. Frery
Abstract:
This work introduces a time-adaptive strategy that uses a refinement estimator based on the first Frenet curvature. In dynamics, a time-adaptive strategy is a mechanism that interactively proposes changes to the time step used in iterative methods of solution. These changes aim to improve the relation between quality of response and computational cost. The method here proposed is suitable for a va…
▽ More
This work introduces a time-adaptive strategy that uses a refinement estimator based on the first Frenet curvature. In dynamics, a time-adaptive strategy is a mechanism that interactively proposes changes to the time step used in iterative methods of solution. These changes aim to improve the relation between quality of response and computational cost. The method here proposed is suitable for a variety of numerical time integration problems, e.g., in the study of bodies subjected to dynamical loads. The motion equation in its space-discrete form is used as reference to derive the formulation presented in this paper. Our method is contrasted with other ones based on local error estimator and apparent frequencies. We check the performance of our proposal when employed with the central difference, the explicit generalized-alpha and the Chung-Lee integration methods. The proposed refinement estimator demands low computational resources, being easily applied to several direct integration methods.
△ Less
Submitted 29 May, 2013;
originally announced May 2013.
-
Analytic Expressions for Stochastic Distances Between Relaxed Complex Wishart Distributions
Authors:
Alejandro C. Frery,
Abraão D. C. Nascimento,
Renato J. Cintra
Abstract:
The scaled complex Wishart distribution is a widely used model for multilook full polarimetric SAR data whose adequacy has been attested in the literature. Classification, segmentation, and image analysis techniques which depend on this model have been devised, and many of them employ some type of dissimilarity measure. In this paper we derive analytic expressions for four stochastic distances bet…
▽ More
The scaled complex Wishart distribution is a widely used model for multilook full polarimetric SAR data whose adequacy has been attested in the literature. Classification, segmentation, and image analysis techniques which depend on this model have been devised, and many of them employ some type of dissimilarity measure. In this paper we derive analytic expressions for four stochastic distances between relaxed scaled complex Wishart distributions in their most general form and in important particular cases. Using these distances, inequalities are obtained which lead to new ways of deriving the Bartlett and revised Wishart distances. The expressiveness of the four analytic distances is assessed with respect to the variation of parameters. Such distances are then used for deriving new tests statistics, which are proved to have asymptotic chi-square distribution. Adopting the test size as a comparison criterion, a sensitivity study is performed by means of Monte Carlo experiments suggesting that the Bhattacharyya statistic outperforms all the others. The power of the tests is also assessed. Applications to actual data illustrate the discrimination and homogeneity identification capabilities of these distances.
△ Less
Submitted 19 April, 2013;
originally announced April 2013.
-
Speckle Reduction in Polarimetric SAR Imagery with Stochastic Distances and Nonlocal Means
Authors:
Leonardo Torres,
Sidnei J. S. Sant'Anna,
Corina C. Freitas,
Alejandro C. Frery
Abstract:
This paper presents a technique for reducing speckle in Polarimetric Synthetic Aperture Radar (PolSAR) imagery using Nonlocal Means and a statistical test based on stochastic divergences. The main objective is to select homogeneous pixels in the filtering area through statistical tests between distributions. This proposal uses the complex Wishart model to describe PolSAR data, but the technique ca…
▽ More
This paper presents a technique for reducing speckle in Polarimetric Synthetic Aperture Radar (PolSAR) imagery using Nonlocal Means and a statistical test based on stochastic divergences. The main objective is to select homogeneous pixels in the filtering area through statistical tests between distributions. This proposal uses the complex Wishart model to describe PolSAR data, but the technique can be extended to other models. The weights of the location-variant linear filter are function of the p-values of tests which verify the hypothesis that two samples come from the same distribution and, therefore, can be used to compute a local mean. The test stems from the family of (h-phi) divergences which originated in Information Theory. This novel technique was compared with the Boxcar, Refined Lee and IDAN filters. Image quality assessment methods on simulated and real data are employed to validate the performance of this approach. We show that the proposed filter also enhances the polarimetric entropy and preserves the scattering information of the targets.
△ Less
Submitted 16 April, 2013;
originally announced April 2013.
-
Classification of Segments in PolSAR Imagery by Minimum Stochastic Distances Between Wishart Distributions
Authors:
Wagner Barreto da Silva,
Corina da Costa Freitas,
Sidnei João Siqueira Sant'Anna,
Alejandro C. Frery
Abstract:
A new classifier for Polarimetric SAR (PolSAR) images is proposed and assessed in this paper. Its input consists of segments, and each one is assigned the class which minimizes a stochastic distance. Assuming the complex Wishart model, several stochastic distances are obtained from the h-phi family of divergences, and they are employed to derive hypothesis test statistics that are also used in the…
▽ More
A new classifier for Polarimetric SAR (PolSAR) images is proposed and assessed in this paper. Its input consists of segments, and each one is assigned the class which minimizes a stochastic distance. Assuming the complex Wishart model, several stochastic distances are obtained from the h-phi family of divergences, and they are employed to derive hypothesis test statistics that are also used in the classification process. This article also presents, as a novelty, analytic expressions for the test statistics based on the following stochastic distances between complex Wishart models: Kullback-Leibler, Bhattacharyya, Hellinger, Rényi, and Chi-Square; also, the test statistic based on the Bhattacharyya distance between multivariate Gaussian distributions is presented. The classifier performance is evaluated using simulated and real PolSAR data. The simulated data are based on the complex Wishart model, aiming at the analysis of the proposal well controlled data. The real data refer to the complex L-band image, acquired during the 1994 SIR-C mission. The results of the proposed classifier are compared with those obtained by a Wishart per-pixel/contextual classifier, and we show the better performance of the region-based classification. The influence of the statistical modeling is assessed by comparing the results using the Bhattacharyya distance between multivariate Gaussian distributions for amplitude data. The results with simulated data indicate that the proposed classification method has a very good performance when the data follow the Wishart model. The proposed classifier also performs better than the per-pixel/contextual classifier and the Bhattacharyya Gaussian distance using SIR-C PolSAR data.
△ Less
Submitted 11 March, 2013;
originally announced March 2013.
-
Entropy-based Statistical Analysis of PolSAR Data
Authors:
Alejandro C. Frery,
Renato J. Cintra,
Abraão D. C. Nascimento
Abstract:
Images obtained from coherent illumination processes are contaminated with speckle noise, with polarimetric synthetic aperture radar (PolSAR) imagery as a prominent example. With an adequacy widely attested in the literature, the scaled complex Wishart distribution is an acceptable model for PolSAR data. In this perspective, we derive analytic expressions for the Shannon, Rényi, and restricted Tsa…
▽ More
Images obtained from coherent illumination processes are contaminated with speckle noise, with polarimetric synthetic aperture radar (PolSAR) imagery as a prominent example. With an adequacy widely attested in the literature, the scaled complex Wishart distribution is an acceptable model for PolSAR data. In this perspective, we derive analytic expressions for the Shannon, Rényi, and restricted Tsallis entropies under this model. Relationships between the derived measures and the parameters of the scaled Wishart law (i.e., the equivalent number of looks and the covariance matrix) are discussed. In addition, we obtain the asymptotic variances of the Shannon and Rényi entropies when replacing distribution parameters by maximum likelihood estimators. As a consequence, confidence intervals based on these two entropies are also derived and proposed as new ways of capturing contrast. New hypothesis tests are additionally proposed using these results, and their performance is assessed using simulated and real data. In general terms, the test based on the Shannon entropy outperforms those based on Rényi's.
△ Less
Submitted 15 October, 2012;
originally announced October 2012.
-
On the Use of Lee's Protocol for Speckle-Reducing Techniques
Authors:
Elsa E. Moschetti,
M. Gabriela Palacio,
Mery Picco,
Oscar H. Bustos,
Alejandro C. Frery
Abstract:
This paper presents two new MAP (Maximum a Posteriori) filters for speckle noise reduction and a Monte Carlo procedure for the assessment of their performance. In order to quantitatively evaluate the results obtained using these new filters, with respect to classical ones, a Monte Carlo extension of Lee's protocol is proposed. This extension of the protocol shows that its original version leads to…
▽ More
This paper presents two new MAP (Maximum a Posteriori) filters for speckle noise reduction and a Monte Carlo procedure for the assessment of their performance. In order to quantitatively evaluate the results obtained using these new filters, with respect to classical ones, a Monte Carlo extension of Lee's protocol is proposed. This extension of the protocol shows that its original version leads to inconsistencies that hamper its use as a general procedure for filter assessment. Some solutions for these inconsistencies are proposed, and a consistent comparison of speckle-reducing filters is provided.
△ Less
Submitted 9 September, 2012;
originally announced September 2012.
-
Assessment of SAR Image Filtering using Adaptive Stack Filters
Authors:
Maria E. Buemi,
Marta Mejail,
Julio Jacobo,
Alejandro C. Frery,
Heitor S. Ramos
Abstract:
Stack filters are a special case of non-linear filters. They have a good performance for filtering images with different types of noise while preserving edges and details. A stack filter decomposes an input image into several binary images according to a set of thresholds. Each binary image is then filtered by a Boolean function, which characterizes the filter. Adaptive stack filters can be design…
▽ More
Stack filters are a special case of non-linear filters. They have a good performance for filtering images with different types of noise while preserving edges and details. A stack filter decomposes an input image into several binary images according to a set of thresholds. Each binary image is then filtered by a Boolean function, which characterizes the filter. Adaptive stack filters can be designed to be optimal; they are computed from a pair of images consisting of an ideal noiseless image and its noisy version. In this work we study the performance of adaptive stack filters when they are applied to Synthetic Aperture Radar (SAR) images. This is done by evaluating the quality of the filtered images through the use of suitable image quality indexes and by measuring the classification accuracy of the resulting images.
△ Less
Submitted 18 July, 2012;
originally announced July 2012.
-
Polarimetric SAR Image Segmentation with B-Splines and a New Statistical Model
Authors:
Alejandro C. Frery,
Julio Jacobo-Berlles,
Juliana Gambini,
Marta Mejail
Abstract:
We present an approach for polarimetric Synthetic Aperture Radar (SAR) image region boundary detection based on the use of B-Spline active contours and a new model for polarimetric SAR data: the GHP distribution. In order to detect the boundary of a region, initial B-Spline curves are specified, either automatically or manually, and the proposed algorithm uses a deformable contours technique to fi…
▽ More
We present an approach for polarimetric Synthetic Aperture Radar (SAR) image region boundary detection based on the use of B-Spline active contours and a new model for polarimetric SAR data: the GHP distribution. In order to detect the boundary of a region, initial B-Spline curves are specified, either automatically or manually, and the proposed algorithm uses a deformable contours technique to find the boundary. In doing this, the parameters of the polarimetric GHP model for the data are estimated, in order to find the transition points between the region being segmented and the surrounding area. This is a local algorithm since it works only on the region to be segmented. Results of its performance are presented.
△ Less
Submitted 17 July, 2012;
originally announced July 2012.
-
Deconvolution of vibroacoustic images using a simulation model based on a three dimensional point spread function
Authors:
Talita Perciano,
Matthew Urban,
Nelson D. A. Mascarenhas,
Mostafa Fatemi,
Alejandro C. Frery,
Glauber T. Silva
Abstract:
Vibro-acoustography (VA) is a medical imaging method based on the difference-frequency generation produced by the mixture of two focused ultrasound beams. VA has been applied to different problems in medical imaging such as imaging bones, microcalcifications in the breast, mass lesions, and calcified arteries. The obtained images may have a resolution of 0.7--0.8 mm. Current VA systems based on co…
▽ More
Vibro-acoustography (VA) is a medical imaging method based on the difference-frequency generation produced by the mixture of two focused ultrasound beams. VA has been applied to different problems in medical imaging such as imaging bones, microcalcifications in the breast, mass lesions, and calcified arteries. The obtained images may have a resolution of 0.7--0.8 mm. Current VA systems based on confocal or linear array transducers generate C-scan images at the beam focal plane. Images on the axial plane are also possible, however the system resolution along depth worsens when compared to the lateral one. Typical axial resolution is about 1.0 cm. Furthermore, the elevation resolution of linear array systems is larger than that in lateral direction. This asymmetry degrades C-scan images obtained using linear arrays. The purpose of this article is to study VA image restoration based on a 3D point spread function (PSF) using classical deconvolution algorithms: Wiener, constrained least-squares (CLSs), and geometric mean filters. To assess the filters' performance, we use an image quality index that accounts for correlation loss, luminance and contrast distortion. Results for simulated VA images show that the quality index achieved with the Wiener filter is 0.9 (1 indicates perfect restoration). This filter yielded the best result in comparison with the other ones. Moreover, the deconvolution algorithms were applied to an experimental VA image of a phantom composed of three stretched 0.5 mm wires. Experiments were performed using transducer driven at two frequencies, 3075 kHz and 3125 kHz, which resulted in the difference-frequency of 50 kHz. Restorations with the theoretical line spread function (LSF) did not recover sufficient information to identify the wires in the images. However, using an estimated LSF the obtained results displayed enough information to spot the wires in the images.
△ Less
Submitted 13 July, 2012;
originally announced July 2012.
-
Stochastic particle packing with specified granulometry and porosity
Authors:
A. C. Frery,
L. Rivarola-Duarte,
V. Carrilho Leão Ramos,
A. Soares Ramos Junior,
W. W. Matos Lira
Abstract:
This work presents a technique for particle size generation and placement in arbitrary closed domains. Its main application is the simulation of granular media described by disks. Particle size generation is based on the statistical analysis of granulometric curves which are used as empirical cumulative distribution functions to sample from mixtures of uniform distributions. The desired porosity i…
▽ More
This work presents a technique for particle size generation and placement in arbitrary closed domains. Its main application is the simulation of granular media described by disks. Particle size generation is based on the statistical analysis of granulometric curves which are used as empirical cumulative distribution functions to sample from mixtures of uniform distributions. The desired porosity is attained by selecting a certain number of particles, and their placement is performed by a stochastic point process. We present an application analyzing different types of sand and clay, where we model the grain size with the gamma, lognormal, Weibull and hyperbolic distributions. The parameters from the resulting best fit are used to generate samples from the theoretical distribution, which are used for filling a finite-size area with non-overlapping disks deployed by a Simple Sequential Inhibition stochastic point process. Such filled areas are relevant as plausible inputs for assessing Discrete Element Method and similar techniques.
△ Less
Submitted 12 July, 2012;
originally announced July 2012.
-
Hypothesis Testing in Speckled Data with Stochastic Distances
Authors:
Abraão D. C. Nascimento,
Renato J. Cintra,
Alejandro C. Frery
Abstract:
Images obtained with coherent illumination, as is the case of sonar, ultrasound-B, laser and Synthetic Aperture Radar -- SAR, are affected by speckle noise which reduces the ability to extract information from the data. Specialized techniques are required to deal with such imagery, which has been modeled by the G0 distribution and under which regions with different degrees of roughness and mean br…
▽ More
Images obtained with coherent illumination, as is the case of sonar, ultrasound-B, laser and Synthetic Aperture Radar -- SAR, are affected by speckle noise which reduces the ability to extract information from the data. Specialized techniques are required to deal with such imagery, which has been modeled by the G0 distribution and under which regions with different degrees of roughness and mean brightness can be characterized by two parameters; a third parameter, the number of looks, is related to the overall signal-to-noise ratio. Assessing distances between samples is an important step in image analysis; they provide grounds of the separability and, therefore, of the performance of classification procedures. This work derives and compares eight stochastic distances and assesses the performance of hypothesis tests that employ them and maximum likelihood estimation. We conclude that tests based on the triangular distance have the closest empirical size to the theoretical one, while those based on the arithmetic-geometric distances have the best power. Since the power of tests based on the triangular distance is close to optimum, we conclude that the safest choice is using this distance for hypothesis testing, even when compared with classical distances as Kullback-Leibler and Bhattacharyya.
△ Less
Submitted 12 July, 2012;
originally announced July 2012.
-
Parametric and Nonparametric Tests for Speckled Imagery
Authors:
Renato J. Cintra,
Abraão D. C. Nascimento,
Alejandro C. Frery
Abstract:
Synthetic aperture radar (SAR) has a pivotal role as a remote imaging method. Obtained by means of coherent illumination, SAR images are contaminated with speckle noise. The statistical modeling of such contamination is well described according with the multiplicative model and its implied G0 distribution. The understanding of SAR imagery and scene element identification is an important objective…
▽ More
Synthetic aperture radar (SAR) has a pivotal role as a remote imaging method. Obtained by means of coherent illumination, SAR images are contaminated with speckle noise. The statistical modeling of such contamination is well described according with the multiplicative model and its implied G0 distribution. The understanding of SAR imagery and scene element identification is an important objective in the field. In particular, reliable image contrast tools are sought. Aiming the proposition of new tools for evaluating SAR image contrast, we investigated new methods based on stochastic divergence. We propose several divergence measures specifically tailored for G0 distributed data. We also introduce a nonparametric approach based on the Kolmogorov-Smirnov distance for G0 data. We devised and assessed tests based on such measures, and their performances were quantified according to their test sizes and powers. Using Monte Carlo simulation, we present a robustness analysis of test statistics and of maximum likelihood estimators for several degrees of innovative contamination. It was identified that the proposed tests based on triangular and arithmetic-geometric measures outperformed the Kolmogorov-Smirnov methodology.
△ Less
Submitted 10 July, 2012;
originally announced July 2012.
-
How good are MatLab, Octave and Scilab for Computational Modelling?
Authors:
Eliana S. de Almeida,
Antonio C. Medeiros,
Alejandro C. Frery
Abstract:
In this article we test the accuracy of three platforms used in computational modelling: MatLab, Octave and Scilab, running on i386 architecture and three operating systems (Windows, Ubuntu and Mac OS). We submitted them to numerical tests using standard data sets and using the functions provided by each platform. A Monte Carlo study was conducted in some of the datasets in order to verify the sta…
▽ More
In this article we test the accuracy of three platforms used in computational modelling: MatLab, Octave and Scilab, running on i386 architecture and three operating systems (Windows, Ubuntu and Mac OS). We submitted them to numerical tests using standard data sets and using the functions provided by each platform. A Monte Carlo study was conducted in some of the datasets in order to verify the stability of the results with respect to small departures from the original input. We propose a set of operations which include the computation of matrix determinants and eigenvalues, whose results are known. We also used data provided by NIST (National Institute of Standards and Technology), a protocol which includes the computation of basic univariate statistics (mean, standard deviation and first-lag correlation), linear regression and extremes of probability distributions. The assessment was made comparing the results computed by the platforms with certified values, that is, known results, computing the number of correct significant digits.
△ Less
Submitted 8 July, 2012;
originally announced July 2012.
-
Nonparametric Edge Detection in Speckled Imagery
Authors:
Edwin Girón,
Alejandro C. Frery,
Francisco Cribari-Neto
Abstract:
We address the issue of edge detection in Synthetic Aperture Radar imagery. In particular, we propose nonparametric methods for edge detection, and numerically compare them to an alternative method that has been recently proposed in the literature. Our results show that some of the proposed methods display superior results and are computationally simpler than the existing method. An application to…
▽ More
We address the issue of edge detection in Synthetic Aperture Radar imagery. In particular, we propose nonparametric methods for edge detection, and numerically compare them to an alternative method that has been recently proposed in the literature. Our results show that some of the proposed methods display superior results and are computationally simpler than the existing method. An application to real (not simulated) data is presented and discussed.
△ Less
Submitted 8 July, 2012;
originally announced July 2012.
-
Polarimetric SAR Image Smoothing with Stochastic Distances
Authors:
Leonardo Torres,
Antonio C. Medeiros,
Alejandro C. Frery
Abstract:
Polarimetric Synthetic Aperture Radar (PolSAR) images are establishing as an important source of information in remote sensing applications. The most complete format this type of imaging produces consists of complex-valued Hermitian matrices in every image coordinate and, as such, their visualization is challenging. They also suffer from speckle noise which reduces the signal-to-noise ratio. Smoot…
▽ More
Polarimetric Synthetic Aperture Radar (PolSAR) images are establishing as an important source of information in remote sensing applications. The most complete format this type of imaging produces consists of complex-valued Hermitian matrices in every image coordinate and, as such, their visualization is challenging. They also suffer from speckle noise which reduces the signal-to-noise ratio. Smoothing techniques have been proposed in the literature aiming at preserving different features and, analogously, projections from the cone of Hermitian positive matrices to different color representation spaces are used for enhancing certain characteristics. In this work we propose the use of stochastic distances between models that describe this type of data in a Nagao-Matsuyama-type of smoothing technique. The resulting images are shown to present good visualization properties (noise reduction with preservation of fine details) in all the considered visualization spaces.
△ Less
Submitted 3 July, 2012;
originally announced July 2012.