Anthony Paul Doulgeris

Estimation of the Equivalent Number of Looks in Polarimetric Synthetic Aperture Radar Imagery

This paper addresses estimation of the equivalent number of looks (ENL), an important parameter in statistical modeling of multilook synthetic aperture radar (SAR) images. Two new ENL estimators are discovered by looking at certain moments of the multilook polarimetric covariance matrix, which is commonly used to represent multilook polarimetric SAR (PolSAR) data, and assuming that the covariance matrix is complex Wishart distributed. First, a second-order trace moment provides a polarimetric extension of the ENL definition and also a matrix-variate version of the conventional ENL estimator. The second estimator is obtained from the log-determinant matrix moment and is also shown to be the maximum likelihood estimator under the Wishart model. It proves to have much lower variance than any other known ENL estimator, whether applied to single-polarization or PolSAR data. Moreover, this estimator is less affected by texture and thus provides more accurate results than other estimators should the assumption of Gaussian statistics for the complex scattering coefficients be violated. These are the first known estimators to use the full covariance matrix as input, rather than individual intensity channels, and therefore to utilize all the statistical information available. We finally demonstrate how an ENL estimate can be computed automatically from the empirical density of small sample estimates calculated over a whole scene. We show that this method is more robust than procedures where the estimate is calculated in a manually selected region of interest.

Classification With a Non-Gaussian Model for PolSAR Data

In this paper, we present a generalized Wishart classifier derived from a non-Gaussian model for polarimetric synthetic aperture radar (PolSAR) data. Our starting point is to demonstrate that the scale mixture of Gaussian (SMoG) distribution model is suitable for modeling PolSAR data. We show that the distribution of the sample covariance matrix for the SMoG model is given as a generalization of the Wishart distribution and present this expression in integral form. We then derive the closed-form solution for one particular SMoG distribution, which is known as the multivariate K-distribution. Based on this new distribution for the sample covariance matrix, termed as the K -Wishart distribution, we propose a Bayesian classification scheme, which can be used in both supervised and unsupervised modes. To demonstrate the effect of including non-Gaussianity, we present a detailed comparison with the standard Wishart classifier using airborne EMISAR data.

Automated Non-Gaussian Clustering of Polarimetric Synthetic Aperture Radar Images

This paper presents an automatic image segmentation method for polarimetric synthetic aperture radar data. It utilizes the full polarimetric information and incorporates texture by modeling with a non-Gaussian distribution for the complex scattering coefficients. The modeling is based upon the well-known product model, with a Gamma-distributed texture parameter leading to the K-Wishart model for the covariance matrix. The automatic clustering is achieved through a finite mixture model estimated with a modified expectation maximization algorithm. We include an additional goodness-of-fit test stage that allows for splitting and merging of clusters. This not only improves the model fit of the clusters, but also dynamically selects the appropriate number of clusters. The resulting image segmentation depicts the statistically significant clusters within the image. A key feature is that the degree of sub-sampling of the input image will affect the detail level of the clustering, revealing only the major classes or a variable level of detail. Real-world examples are shown to demonstrate the technique.

Leads in Arctic pack ice enable early phytoplankton blooms below snow-covered sea ice

The Arctic icescape is rapidly transforming from a thicker multiyear ice cover to a thinner and largely seasonal first-year ice cover with significant consequences for Arctic primary production. One critical challenge is to understand how productivity will change within the next decades. Recent studies have reported extensive phytoplankton blooms beneath ponded sea ice during summer, indicating that satellite-based Arctic annual primary production estimates may be significantly underestimated. Here we present a unique time-series of a phytoplankton spring bloom observed beneath snow-covered Arctic pack ice. The bloom, dominated by the haptophyte algae Phaeocystis pouchetii, caused near depletion of the surface nitrate inventory and a decline in dissolved inorganic carbon by 16 ± 6 g C m−2. Ocean circulation characteristics in the area indicated that the bloom developed in situ despite the snow-covered sea ice. Leads in the dynamic ice cover provided added sunlight necessary to initiate and sustain the bloom. Phytoplankton blooms beneath snow-covered ice might become more common and widespread in the future Arctic Ocean with frequent lead formation due to thinner and more dynamic sea ice despite projected increases in high-Arctic snowfall. This could alter productivity, marine food webs and carbon sequestration in the Arctic Ocean.

Scale mixture of Gaussian modelling of polarimetric SAR data

This paper describes a flexible non-Gaussian statistical method used to model polarimetric synthetic aperture radar (POLSAR) data. We outline the theoretical basis of the well-know product model as described by the class of Scale Mixture models and discuss their appropriateness for modelling radar data. The statistical distributions of several Scale mixture models are then described, including the commonly used Gaussian model, and techniques for model parameter estimation are given. Real data evaluations are made using airborne fully polarimetric SAR studies for several distinct land cover types. Generic scale mixture of Gaussian features is extracted from the model parameters and a simple clustering example presented.

Goodness-of-Fit Tests for Multilook Polarimetric Radar Data Based on the Mellin Transform

The advent of polarimetric synthetic aperture radar has spurred a growing interest in statistical models for complex-valued covariance matrices, which is the common representation of multilook polarimetric radar images. In this paper, we respond to an emergent need by proposing statistical tests for the simple and composite goodness-of-fit (GoF) problem for a class of compound matrix distributions. The tests are based on Mellin-kind matrix cumulants. These are derived from a novel characteristic function for positive definite Hermitian random matrices, defined in terms of a matrix-variate Mellin transform instead of the conventional Fouriér transform, and belong to a new framework for statistical analysis of multilook polarimetric radar data recently introduced by the authors. The cumulant-based tests are easy to compute, and the asymptotic sampling distribution of the test statistic is chi-square distributed in the simple hypothesis case. Under the composite hypothesis, the sampling distribution is obtained by Monte Carlo simulations. We evaluate the power of the proposed GoF tests with simulated data. We also use them to assess the fit of several matrix distributions to real data acquired by Radarsat-2 in fine-quad polarization mode.

Correction and Characterization of Radio Frequency Interference Signatures in L-Band Synthetic Aperture Radar Data

Radio frequency interference (RFI) is a known issue in low-frequency radar remote sensing. In synthetic aperture radar (SAR) image processing, RFI can cause severe degradation of image quality, distortion of polarimetric signatures, and an increase of the SAR phase noise level. To address this issue, a processing system was developed that is capable of reliably detecting, characterizing, and mitigating RFI signatures in SAR observations. In addition to being the basis for image correction, the robust RFI-detection algorithms developed in this paper are used to retrieve a wealth of RFI-related information that allows for mapping, characterizing, and classifying RFI signatures across large spatial scales. The extracted RFI information is expected to be valuable input for SAR-system design, sensor operations, and the development of effective RFI-mitigation strategies. The concepts of RFI detection, analysis, and mapping are outlined. Large-scale RFI mapping results are shown. In case studies, the benefit of detailed RFI information for customized RFI filtering and sensor operations is exemplified.

A Textural–Contextual Model for Unsupervised Segmentation of Multipolarization Synthetic Aperture Radar Images

This paper proposes a novel unsupervised, non-Gaussian, and contextual segmentation method that combines an advanced statistical distribution with spatial contextual information for multilook polarimetric synthetic aperture radar (PolSAR) data. This extends on previous studies that have shown the added value of both non-Gaussian modeling and contextual smoothing individually or for intensity channels only. The method is based on a Markov random field (MRF) model that integrates a K-Wishart distribution for the PolSAR data statistics conditioned to each image cluster and a Potts model for the spatial context. Specifically, the proposed algorithm is constructed based upon the stochastic expectation maximization (SEM) algorithm. A new formulation of SEM is developed to jointly perform clustering of the data and parameter estimation of the K-Wishart distribution and the MRF model. Experiments on simulated and real PolSAR data demonstrate the added value of using an appropriate statistical representation, in combination with contextual smoothing.

An Automatic

We have recently presented a novel unsupervised, non-Gaussian, and contextual clustering algorithm for segmentation of polarimetric synthetic aperture radar (PolSAR) images. This represents one of the most advanced PolSAR unsupervised statistical segmentation algorithms and uses the doubly flexible two-parameter U-distribution model for the PolSAR statistics and includes a Markov random field (MRF) approach for contextual smoothing. A goodness-of-fit testing stage adds a statistically rigorous approach to determine the significant number of classes. The fully automatic algorithm was demonstrated with good results for both simulated and real data sets. This paper discusses a rethinking of the overall strategy and leads to some simplifications. The primary issue was that the MRF optimization depends on the number of classes and did not behave well under the split-and-merge environment. We explain the reasons behind a separation of the cluster evaluation from the contextual smoothing and a modified rationale for the adaptive number of classes. Both aspects have simplified the overall algorithm while maintaining good visual results.

{\cal U}

This paper presents a processing chain for the change detection of Arctic glaciers from multitemporal multipolarization synthetic aperture radar (SAR) images. We produce terrain-corrected multilook complex covariance data by including the effects of topography on both geolocation and SAR radiometry as well as azimuth slope variations on polarization signature. An unsupervised contextual non-Gaussian clustering algorithm is employed for the segmentation of each terrain-corrected polarimetric SAR image and subsequently labeled with the aid of ground-truth data into glacier facies. We demonstrate the consistency of the segmentation algorithm by characterizing the expected random error level for different SAR acquisition conditions. This allows us to determine whether an observed variation is statistically significant and therefore can be used for the postclassification change detection of Arctic glaciers. Subsequently, the average classified images of succeeding years are compared, and changes are identified as the detected differences in the location of boundaries between glacier facies. In the current analysis, a series of dual-polarization C-band ENVISAT ASAR images over the Kongsvegen glacier, Svalbard, is used for demonstration.