Robust CFAR Detector Based on Truncated Statistics for Polarimetric Synthetic Aperture Radar

Constant false alarm rate (CFAR) algorithms using a local training window are widely used for ship detection with synthetic aperture radar (SAR) imagery. However, when the density of the targets is high, such as in busy shipping lines and crowded harbors, the background statistics may be contaminated by the presence of nearby targets in the training window. Recently, a robust CFAR detector based on truncated statistics (TS) was proposed. However, the truncation of data in the format of polarimetric covariance matrices is much more complicated with respect to the truncation of intensity (single polarization) data. In this article, a polarimetric whitening filter TS CFAR (PWF-TS-CFAR) is proposed to estimate the background parameters accurately in the contaminated sea clutter for PolSAR imagery. The CFAR detector uses a polarimetric whitening filter (PWF) to turn the multidimensional problem to a 1-D case. It uses truncation to exclude possible statistically interfering outliers and uses TS to model the remaining background samples. The algorithm does not require prior knowledge of the interfering targets, and it is performed iteratively and adaptively to derive better estimates of the polarimetric covariance matrix (although this is computationally expensive). The PWF-TS-CFAR detector provides accurate background clutter modeling, a stable false alarm property, and improves the detection performance in high-target-density situations. RadarSat2 data are used to verify our derivations, and the results are in line with the theory

CFAR Ship Detection in Polarimetric Synthetic Aperture Radar Images Based on Whitening Filter

Polarimetric whitening filter (PWF) can be used to filter polarimetric synthetic aperture radar (PolSAR) images to improve the contrast between ships and sea clutter background. For this reason, the output of the filter can be used to detect ships. This paper deals with the setting of the threshold over PolSAR images filtered by the PWF. Two parameter-constant false alarm rate (2P-CFAR) is a common detection method used on whitened polarimetric images. It assumes that the probability density function (PDF) of the filtered image intensity is characterized by a log-normal distribution. However, this assumption does not always hold. In this paper, we propose a systemic analytical framework for CFAR algorithms based on PWF or multi-look PWF (MPWF). The framework covers the entire log-cumulants space in terms of the textural distributions in the product model, including the constant, gamma, inverse gamma, Fisher, beta, inverse beta, and generalized gamma distributions (GΓDs). We derive the analytical forms of the PDF for each of the textural distributions and the probability of false alarm (PFA). Finally, the threshold is derived by fixing the false alarm rate (FAR). Experimental results using both the simulated and real data demonstrate that the derived expressions and CFAR algorithms are valid and robust

Optimal Parameter Estimation in Heterogeneous Clutter for High Resolution Polarimetric SAR Data

International audienceThis letter presents a new estimation scheme for optimally deriving clutter parameters with high-resolution polarimetric synthetic aperture radar (POLSAR) data. The heterogeneous clutter in POLSAR data is described by the spherically invariant random vector model. Three parameters are introduced for the high-resolution POLSAR data clutter: the span, the normalized texture, and the speckle normalized covariance matrix. The asymptotic distribution of the novel span estimator is investigated. A novel heterogeneity test for the POLSAR clutter is also discussed. The proposed method is tested with airborne POLSAR images provided by the Office National d'Études et de Recherches Aerospatiales Radar Aéroporté Multi-spectral d'Etude des Signatures system

Comparison of Target Detectors to Identify Icebergs in Quad-Polarimetric L-Band Synthetic Aperture Radar Data

Icebergs represent hazards to ships and maritime activities and therefore their detection is essential. Synthetic Aperture Radar (SAR) satellites are very useful for this, due to their capability to acquire data under cloud cover and during day and night passes. In this work, we compared six state-of-the-art polarimetric target detectors to test their performance and ability to detect small-sized icebergs 120 m, as they are easier to detect). However, the differences between quad- and dual- or single-polarimetric detectors became much more evident when the PF value was fixed to low detection probabilities 10−6 (i.e., smaller icebergs). In the single-polarimetric mode, the HV channel showed PD values of 0.62 for open ocean and 0.26 for sea ice, compared to values of 0.81 (open ocean) and 0.77 (sea ice) obtained with quad-polarimetric detectors

The Polarimetric Detection Optimization Filter and Its Statistical Test for Ship Detection

Ship detection via synthetic aperture radar (SAR) has been demonstrated to be very useful as polarimetric information helps discriminate between targets and sea clutter. Among the available polarimetric detectors, optimal polarimetric detection (OPD) theoretically provides the best detection performance under the assumption that the fully developed speckle hypothesis stands. This study proposes a polarimetric detection optimization filter (PDOF). The target clutter ratio (TCR) over the speckle variation was maximized using a matrix transform to derive the PDOF. The objective function based on a matrix transform instead of a vector transform is optimized to obtain synthetic effects by combining a polarimetric whitening filter (PWF) and a polarimetric matched filter (PMF). Subspace form of the PDOF (SPDOF) is also proposed, which gives performance comparable to the PDOF. Assuming a Wishart distribution, the exact and approximate expressions of the closed-form probability density function (PDF) of the PDOF are derived. The probability of false alarm (PFA) was derived in a closed-form expression, which allows obtaining the PDOF threshold analytically. Moreover, the gamma model is extended to a generalized gamma distribution (GΓD) to adapt complicated resolutions and sea states. Experiments with simulated and real data validate the correctness and effectiveness of the results. The PDOF detector achieves the best performance in most virtual and real-world environments, especially in cases where the target statistics and clutter are not Wishart-distributed

SCHATTEN MATRIX NORM BASED POLARIMETRIC SAR DATA REGULARIZATION. APPLICATION OVER CHAMONIX MONT-BLANC

International audienceThe paper addresses the filtering of Polarimetry Synthetic Aperture Radar (PolSAR) images. The filtering strategy is based on a regularizing cost function associated with matrix norms called the Schatten p-norms. These norms apply on matrix singular values. The proposed approach is illustrated upon scattering and coherency matrices on RADARSAT-2 PolSAR images over the Chamonix Mont-Blanc site. Several p values of Schatten p-norms are surveyed and their capabilities on filtering PolSAR images is provided in comparison with conventional strategies for filtering PolSAR data

PolSAR Ship Detection Based on Neighborhood Polarimetric Covariance Matrix

The detection of small ships in polarimetric synthetic aperture radar (PolSAR) images is still a topic for further investigation. Recently, patch detection techniques, such as superpixel-level detection, have stimulated wide interest because they can use the information contained in similarities among neighboring pixels. In this article, we propose a novel neighborhood polarimetric covariance matrix (NPCM) to detect the small ships in PolSAR images, leading to a significant improvement in the separability between ship targets and sea clutter. The NPCM utilizes the spatial correlation between neighborhood pixels and maps the representation for a given pixel into a high-dimensional covariance matrix by embedding spatial and polarization information. Using the NPCM formalism, we apply a standard whitening filter, similar to the polarimetric whitening filter (PWF). We show how the inclusion of neighborhood information improves the performance compared with the traditional polarimetric covariance matrix. However, this is at the expense of a higher computation cost. The theory is validated via the simulated and measured data under different sea states and using different radar platforms

Quad polarimetric synthetic aperture radar analysis of icebergs in Greenland and Svalbard

Polarimetric synthetic aperture radar (PolSAR) has been widely used in ocean and cryospheric applications. This is because, PolSAR can be used in all-day operations and in areas of cloud cover, and therefore can provide valuable large-scale monitoring in polar regions, which is very helpful to shipping and offshore maritime operations. In the last decades, attention has turned to the potential of PolSAR to detect icebergs in the Arctic since they are a major hazard to vessels. However, there is a substantial lack of literature exploring the potentialities of PolSAR and the understanding of iceberg scattering mechanisms. Additionally, it is not known if high resolution PolSAR can be used to detect icebergs smaller than 120 metres. This thesis aims to improve the knowledge of the use of PolSAR scattering mechanisms of icebergs, and detection of small icebergs. First, an introduction to PolSAR is outlined in chapter two, and monitoring of icebergs is presented in chapter three. The first data chapter (Chapter 4) is focused on developing a multi-scale analysis of icebergs using parameters from the Cloude-Pottier and the Yamaguchi decompositions, the polarimetric span and the Pauli scattering vector. This method is carried out using ALOS-2 PALSAR quad polarimetric L-band SAR on icebergs in Greenland. This approach outlines the good potential for using PolSAR for future iceberg classification. One of the main important outcomes is that icebergs are composed by a combination of single targets, which therefore may require a more complex way of processing SAR data to properly extract physical information. In chapter five, the problem of detecting icebergs is addressed by introducing six state-of-the-art detectors previously applied to vessel monitoring. These detectors are the Dual Intensity Polarisation Ratio Anomaly Detector (iDPolRAD), Polarimetric Notch Filter (PNF), Polarimetric Matched Filter (PMF), reflection symmetry (sym), Optimal Polarimetric Detector (OPD) and the Polarimetric Whitening Filter (PWF). Cloude-Pottier entropy, and first and third eigenvalues (eig1 and eig3) of the coherency matrix are also utilised as parameters for comparison. This approach uses the same ALOS-2 dataset, but also evaluates detection performance in two scenarios: icebergs in open ocean, and in sea ice. Polarimetric modes (quad-pol, dual-pol, and single intensities) are also considered for comparison. Currently it is very difficult to detect icebergs less than 120 metres in length using this approach, due to the scattering mechanisms of icebergs and sea ice being very similar. However, it was possible to obtain detection performances of the OPD and PWF, which both showed a Probability of Detection (PF) of 0.99 when the Probability of False Alarms (PF) was set to 10-5 in open ocean. Similarly, in dual pol images, the PWF gave the best performance with a PD of 0.90. Results in sea ice found eig3 to be the best detector with a PD of 0.90 while in dual-pol mode, iDPolRAD gave a PD of 0.978. Single intensity detector performance found the HV channel gave the best detection with a PD of 0.99 in open ocean and 0.87 in sea ice. In the previous two approaches, only satellite data is used. However, in chapter six, data from a ground-based Ku-band Gamma Portable Radio Interferometer (GPRI) instrument is introduced, providing images that are synchronised with the satellite acquisitions. In this approach, the same six detectors are applied to three multitemporal RADARSAT-2 quad pol C-band SAR images on icebergs in Kongsfjorden, Svalbard to evaluate the detection performance within a changing fjord environment. As before, we also make use of Cloude-Pottier entropy, eig1 and eig3. Finally, we evaluate the target-to-clutter ratio (TCR) of the icebergs and check for correlation between the backscattering coefficients and the iceberg dimension. The results obtained from this thesis present original additions to the literature that contributes to the understanding of PolSAR in cryospheric applications. Although these methods are applied to PolSAR and ground-based radar on vessels, they have been applied for the first time on icebergs in this thesis. To summarise, the main findings are that icebergs cannot be represented as single or partial targets, but they do exhibit a collection of single targets clustered together. This result leads to the fact that entropy is not sufficient as a parameter to detect icebergs. Detection results show that the OPD and PWF detectors perform best in an open ocean setting and using quad-pol mode. These results are degraded in dual-pol mode, while single intensity detection is best in the HV cross polarisation channel. When these detectors are applied to the RADARSAT-2 in Svalbard, the OPD and PWF detectors also perform best with PD values ranging between 0.5-0.75 for a PF of 0.01-0.05. However, the sea ice present in the fjord degrades performance across all detectors. Correlation plots with iceberg size show that a regression is not straightforward and Computer Vision methodologies may work best for this

The Performance Analysis Based on SAR Sample Covariance Matrix

Multi-channel systems appear in several fields of application in science. In the Synthetic Aperture Radar (SAR) context, multi-channel systems may refer to different domains, as multi-polarization, multi-interferometric or multi-temporal data, or even a combination of them. Due to the inherent speckle phenomenon present in SAR images, the statistical description of the data is almost mandatory for its utilization. The complex images acquired over natural media present in general zero-mean circular Gaussian characteristics. In this case, second order statistics as the multi-channel covariance matrix fully describe the data. For practical situations however, the covariance matrix has to be estimated using a limited number of samples, and this sample covariance matrix follow the complex Wishart distribution. In this context, the eigendecomposition of the multi-channel covariance matrix has been shown in different areas of high relevance regarding the physical properties of the imaged scene. Specifically, the maximum eigenvalue of the covariance matrix has been frequently used in different applications as target or change detection, estimation of the dominant scattering mechanism in polarimetric data, moving target indication, etc. In this paper, the statistical behavior of the maximum eigenvalue derived from the eigendecomposition of the sample multi-channel covariance matrix in terms of multi-channel SAR images is simplified for SAR community. Validation is performed against simulated data and examples of estimation and detection problems using the analytical expressions are as well given