CFAR Edge Detector for Polarimetric SAR Images

Abstract—Finding the edges between different regions in an image is one of the fundamental steps of image analysis, and several edge detectors suitable for the special statistics of synthetic aperture radar (SAR) intensity images have previously been developed. In this paper, a new edge detector for polarimetric SAR images is presented using a newly developed test statistic in the complex Wishart distribution to test for equality of covariance matrices. The new edge detector can be applied to a wide range of SAR data from single-channel intensity data to multifrequency and/or multitemporal polarimetric SAR data. By simply changing the parameters characterizing the test statistic according to the applied SAR data, constant false-alarm rate detection is always obtained. An adaptive filtering scheme is presented, and the distributions of the detector are verified using simulated polarimetric SAR images. Using SAR data from the Danish airborne polarimetric SAR, EMISAR, it is demonstrated that superior edge detection results are obtained using polarimetric and/or multifrequency data compared to using only intensity data. Index Terms—Complex Wishart distribution, edge detection, polarimetry, synthetic aperture radar (SAR), Wishart likelihoodrati

Signal subspace change detection in structured covariance matrices

International audienceTesting common properties between covariance matricesis a relevant approach in a plethora of applications. In thispaper, we derive a new statistical test in the context of structuredcovariance matrices. Specifically, we consider low rank signalcomponent plus white Gaussian noise structure. Our aim is totest the equality of the principal subspace, i.e., subspace spannedby the principal eigenvectors of a group of covariance matrices. Adecision statistic is derived using the generalized likelihood ratiotest. As the formulation of the proposed test implies a non-trivialoptimization problem, we derive an appropriate majorizationminimizationalgorithm. Finally, numerical simulations illustratethe properties of the newly proposed detector compared to thestate of the art

A multi-family GLRT-based algorithm for oil spill detection

This paper deals with detection of oil spills from multi-polarization SAR images. The problem is cast in terms of a composite hypothesis test aimed at discriminating between the Polarimetric Covariance Matrix (PCM) equality (absence of oil spills in the tested region) and the situation where the region under test exhibits a PCM with at least an ordered eigenvalue smaller than that of a reference covariance. This last setup reflects the physical condition where the back scattering associated with the oil spills leads to a signal, in some eigen-directions, weaker than the one gathered from a reference area where it is a-priori known the absence of any oil slicks. A Multi-family Generalized Likelihood Ratio Test (MGLRT) approach is pursued to come up with an adaptive detector ensuring the Constant Alarm False Rate (CFAR) property. At the analysis stage, the behavior of the new architecture is investigated in comparison with a benchmark (but non-implementable) structure and some other sub-optimum adaptive detectors available in open literature. The study, conducted in the presence of both simulated and real data, confirms the practical effectiveness of the new approach

Trace Coherence: A New Operator for Polarimetric and Interferometric SAR images

Quadratic forms play an important role in the development of several Polarimetric and Interferometric Synthetic Aperture Radar (Pol-InSAR) methodologies, which are very powerful tools for Earth Observation. This work investigates integrals of Pol-InSAR operators based on quadratic forms, with special interest for the Pol-InSAR coherence. A new operator is introduced, namely Trace Coherence, that provides an approximation for the center of mass of the Coherence Region (CoRe). The latter is the locus of points on the polar plot containing all the possible coherence values. Such center of mass can be calculated as the integral of Pol-InSAR coherences over the scattering mechanisms. The Trace Coherence provides a synthetic information regarding the partial target as one single entity. Therefore, it provides a representation, which is not dependent on the selection of one specific polarization channel. It may find application in change detection (e.g. Coherent Change Detection and differential DEM), classification (e.g. building structure parameters) and modeling (e.g. for the retrieval of forest height). In calculating the integral of the Pol-InSAR coherences, an approximate Trace Coherence expression is derived and shown to improve the calculation speed by several orders of magnitude. The Trace Coherence approximation is investigated using Monte Carlo simulations and validated ESA (DLR) L-band quad-polarimetric data acquired during the AGRISAR 2006 campaign. The result of the analysis using simulated and real data is that the average error in approximating the integral of the Coherence Region is 0.025 in magnitude and 3° in phase (in scenarios with sufficiently high coherence)

Two-Stage Change Detection for Synthetic Aperture Radar

Coherent change detection using paired synthetic aperture radar (SAR) images is often performed using a classical coherence estimator that is invariant to the true variances of the populations underlying each paired sample. While attractive, this estimator is biased and requires a significant number of samples to yield good performance. Increasing sample size often results in decreased image resolution. Thus, we propose the use of Berger's coherence estimate because, with the same number of pixels, the estimator effectively doubles the sample support without sacrificing resolution when the underlying population variances are equal or near equal. A potential drawback of this approach is that it is not invariant since its distribution depends on the pixel pair population variances. While Berger's estimator is inherently sensitive to the inequality of population variances, we propose a method of insulating the detector from this acuity. A two-stage change statistic is introduced to combine a noncoherent intensity change statistic given by the sample variance ratio, followed by the alternative Berger estimator, which assumes equal population variances. The first-stage detector identifies pixel pairs that have nonequal variances as changes caused by the displacement of sizeable object. The pixel pairs that are identified to have equal or near-equal variances in the first stage are used as an input to the second stage. The second-stage test uses the alternative Berger coherence estimator to detect subtle changes such as tire tracks and footprints. We show experimentally that the proposed method yields higher contrast SAR change detection images than the classical coherent change detector (state of the art), the alternative coherent change detector, and the intensity change detector. Experimental results are presented to show the effectiveness and robustness of the proposed algorithm for SAR change detection

Polarization Optimization for the Detection of Multiple Persistent Scatterers Using SAR Tomography

The detection of multiple interfering persistent scatterers (PSs) using Synthetic Aperture Radar (SAR) tomography is an efficient tool for generating point clouds of urban areas. In this context, detection methods based upon the polarization information of SAR data are effective at increasing the number of PSs and producing high-density point clouds. This paper presents a comparative study on the effects of the polarization design of a radar antenna on further improving the probability of detecting persistent scatterers. For this purpose, we introduce an extension of the existing scattering property-based generalized likelihood ratio test (GLRT) with realistic dependence on the transmitted/received polarizations. The test is based upon polarization basis optimization by synthesizing all possible polarimetric responses of a given scatterer from its measurements on a linear orthonormal basis. Experiments on both simulated and real data show, by means of objective metrics (probability of detection, false alarm rate, and signal-to-noise ratio), that polarization waveform optimization can provide a significant performance gain in the detection of multiple scatterers compared to the existing full-polarization-based detection method. In particular, the increased density of detected PSs at the studied test sites demonstrates the main contribution of the proposed method

Forcing scale invariance in multipolarization SAR change detection

This paper considers the problem of coherent (in the sense that both amplitudes and relative phases of the polarimetric returns are used to construct the decision statistic) multi-polarization SAR change detec- tion starting from the availability of image pairs exhibiting possible power mismatches/miscalibrations. The principle of invariance is used to characterize the class of scale-invariant decision rules which are insensitive to power mismatches and ensure the Constant False Alarm Rate (CFAR) property. A maximal invariant statistic is derived together with the induced maximal invariant in the parameter space which significantly compress the data/parameters domain. A Generalized Likelihood Ratio Test (GLRT) is synthesized both for the cases of two- and three-polarimetric channels. Interestingly, for the two-channel case, it is based on the comparison of the condition number of a data-dependent matrix with a suitable threshold. Some additional invariant decision rules are also proposed. The performance of the considered scale-invariant structures is compared to those from two non- invariant counterparts using both simulated and real radar data. The results highlight the robustness of the proposed method and the performance tradeoff involve