A constant false alarm rate (CFAR) detecting method for ships in high-resolution dual-polarization synthetic aperture radar (SAR) amplitude images has been proposed in this paper. First, by the production of amplitude images from two polarimetric channels, a novel detector simply called the PMA detector has been constructed. We testified that the PMA detector could improve the signal-to-clutter ratio (SCR) and make the discrimination of a ship from clutter more easily. Second, the PMA detector’s statistical model has been described by the well-known 𝒢0 distribution when facing complex sea background. The experiments performed on measured dual-polarization TerraSAR-X images demonstrate the good performance of the proposed CFAR detecting method.
1. Introduction
Ship detection in high-resolution synthetic aperture radar (SAR) images has become an increasing interest during the past several years [1]. It plays an important role for various potential applications like marine surveillance, vessel traffic control, military service, and so forth [1–3]. Additionally, with the development of sensor techniques, the advanced polarimetric SAR systems have also been validated that more completed description of target-scattering behavior could be provided than traditional single-channel SAR systems, such as HH, HV, and VV SAR, which is not sufficient for ship detection [2, 4]. Since better discriminating target signal from the surrounding clutter can be obtained by polarimetric SAR systems, ship detection in multipolarimetric or dual-polarimetric SAR images is receiving intensive attentions in present.
Nowadays, some famous contemporary satellite SAR systems, for example, TerraSAR-X, COSMO-SkyMed, and RadarSat-2, support dual-polarization acquisition mode [5]. However, the available data is, mostly, the amplitude data and not the complex-valued data [6]. It is an important data modality because of the image products provided by many satellite SAR systems, for instance, ERS. Thus, we mainly focus on ship detection using only dual-pol SAR amplitude images in this paper.
A wide variety of methods have been proposed for detecting ships in polarimetric SAR images. The basic idea is to reduce three channels of polarimetric SAR data to a single decision criterion [7]. Some popular detectors [7], including the optimal polarimetric detector (OPD), the polarimetric whitening filter (PWF), the span detector, the power maximization synthesis (PMS) detector, and, more recently, the identity likelihood ratio test (ILRT), have been developed and shown that they can perform in a way that targets are more easily discriminated from clutter. Unfortunately, other detectors mentioned earlier except for the span detector are very difficult to implement due to the missing of phase information when the only amplitude is considered [7].
The span detector is a widely used processor which is a noncoherent sum of all polarimetric channels and only makes use of image intensities or amplitudes. It has also been proven that this detector can acquire a better detection performance than HH, HV, or VV individually. Nevertheless, the expected adaptive detection threshold is impossible when facing complex clutter background by this detector, which limits the applications of this detector.
Our goal is to present a powerful detector for ship detection in high-resolution dual-polarization SAR amplitude images. On one hand, this detector can improve the signal-to-clutter ratio (SCR) to enhance the moving targets or restrain the clutter. Meanwhile, we wish that a flexible and adaptive constant false alarm rate (CFAR) threshold could be derived from this detector. Under this consideration, this paper proposed a novel detector similarly with the span detector, simply called the product of multilook amplitudes (PMAs) detector.
2. Dual-Polarization SAR Data Description
The fundamental quantities measured by a polarimetric SAR are the complex scattering matrix with four elements, which, in complete form, can be expressed by [8]
(1)S=[ShhShvSvhSvv],
where Spq is the scattering element with p transmit and q receive polarizations; p and q denote either horizontal (h) or vertical (v) polarization. In a reciprocal medium, the cross-polar elements of the scattering are identical; that is, Shv=Svh. In this matter, the scattering matrix shown in (1) can be organized into a column vector
(2)X=[ShhShvSvv]T,
where the superscript “T” represents transpose. The equation (2) is called single-look complex scattering vector. Based on the coherent nature of SAR, X follows a zero mean multivariate complex Gaussian distribution. The detailed discussion about this distribution can be found in [8] and [10].
2.1. Polarimetric Covariance Matrix for Dual-Polarization Case
For dual-polarization case, the single look scattering vector shown in (2) can be simplified to
(3)u=[S1S2]T.
Herein, for convenience, we use S1 or S2 to indicate one of the scattering elements Shh, Shv, and Svv in any order, as well as S1≠S2. Additionally, to reduce the influence of speckle, SAR data are often multi-look averaged. As the polarimetric information can also be represented by a covariance matrix, the n-look sample covariance matrix is defined as [8, 9]
(4)R=1n∑k=1nu(k)u(k)H=1n∑k=1n[|S1(k)|2S1(k)S2(k)*S1(k)*S2(k)|S2(k)|2],
where the superscript * means complex conjugate and H represents conjugate complex transpose, n is the number of looks, and u(k)=[S1(k),S2(k)]T is the kth single-look image. Assuming statistical ergodicity and constant RCS background, the random matrix R, known as the complex Wishart distribution [10], is with probability density
(5)pR(R)=n2ndet(R)n-2exp[-nTr(C-1R)]πΓ(n)Γ(n-1)det(C)n,
where Γ(·) is the gamma function, and Tr(·) indicates the matrix trace. The symbol det(·) denotes the determinant operator, and the covariance matrix is 2×2 complex, Hermitian, written as
(6)C=E[uuH]=[E(|S1|2)E(|S1|2)E(|S2|2)ρejθE(|S1|2)E(|S2|2)ρe-jθE(|S2|2)],
where ρejθ is the complex correlation coefficient of two components in (3).
2.2. The Joint Distribution of Two Multilook Intensities from Different Polarimetric Channels
When only incomplete polarimetric data are available, for instance, the amplitude or intensity of copolarized components (i.e., HH and VV), the joint distribution of intensity, or amplitude from two correlated polarimetric channels is of importance for constructing signal processing algorithms such as detection and classification in this case. Lee et al. [8] have derived the PDF of joint normalized multilook intensities by integrating (5) with respect to the two off-diagonal elements, which is modeled as
(7)pR1,R2(R1,R2)=nn+1(R1R2)(n-1)/2exp(-n(R1+R2)/(1-ρ2))Γ(n)(1-ρ2)ρn-1×In-1(2nR1R2ρ1-ρ2),R1,R2,n>0,0<ρ<1,
where In-1(·) is the first type modified Bessel function of order n-1, Ri,i=1,2, represents normalized multi-look intensity of ith polarimetric channel with the expression Ri=(1/n)∑k=1n(|Si(k)|2/E(|Si|2)).
Furthermore, to facilitate the posterior derivation, by the transforms with the following forms:
(8)B1=nR1,B2=nR2,
one can obtain the joint density of B1 and B2 as (see [8] for details)
(9)pB1,B2(B1,B2)=(B1B2)(n-1)/2exp(-(B1+B2)/(1-ρ2))Γ(n)(1-ρ2)ρn-1×In-1(2B1B2ρ1-ρ2),B1,B2,n>0,0<ρ<1.
3. The PMA Detector3.1. PMA Detector
In a single-channel SAR image, it is usually assumed that strong backscattering comes from targets. A target point will be lost when the backscattering amplitude or intensity from the target is not large enough compared with the clutter background, generally, due to low signal-to-clutter ratio (SCR). In other words, SCR is an essential factor influencing the detection performance when only amplitude or intensity data are available. Thus, the principle designing a good detector should naturally enable SCR to be improved, that is, enhancing target and restraining clutter.
As we know, the span detector is a widely used processor which is a noncoherent sum of all polarimetric channels and only makes use of image intensities. For dual-polarization multilook case, the span is given by [7]
(10)span=1n∑k=1n|S1(k)|2+1n∑k=1n|S2(k)|2.
This detector can be regarded as the synthetic power of all channels. Consequently, some investigations have shown that a lower noise level and a higher SCR can be obtained by this detector than HH, HV, or VV individually. This conclusion implies that the way of synthetic power can arrive at the purpose of improving SCR, so that the targets can be more easily discriminated from the clutter, compared with that only arbitrary single-channel information is used. Meanwhile, it is also very hard to adaptively give a proper detection threshold by the span detector due to the unknown knowledge of the corresponding statistics. Motivated by these considerations, we construct a detector in this paper by means of another synthetic power, that is, the product of multilook amplitudes from two polarimetric channels, for convenience, called PMA detector which can be defined as
(11)ξ=σ1σ2R1R2=σ1σ2B1B2n,
where σi=E(|Si|2). From an intuitive understanding, in the variable ξ domain, for the targets, their amplitudes of two polarimetric channels are both larger than the surrounding clutter, which results in a much faster cumulative speed of power for targets verse clutter by multiplying the first image amplitude by the second image one, and, hence, targets’ values of ξ are much more prominent and target signal can be enhanced.
3.2. The CFAR Algorithm of PMA Detector
For dual-pol SAR amplitude data, based on the multiplicative model and the assumption that the energy of two channels is balanced, the product of multilook amplitudes from two polarimetric channels can be denoted as
(12)ζ=1n∑k=1n|A1X1(k)|21n∑k=1n|A2X2(k)|2=A12A221n∑k=1n|X1(k)|21n∑k=1n|X2(k)|2=A1A2ξ,
where Ai represents the backscattering RCS amplitude component of ith receiving polarimetric channel, and Xi(k)=|Si(k)|2/E(|Si|2). Recently, Frery et al. [11] have proposed a well-known 𝒢0 distribution to model the clutter regions in SAR images, and the product of multilook amplitudes from two polarimetric channels, that is, the proposed detector in this paper, employs an intensity expression of this distribution with three distinct parameters σ, n, and α, when applying to the modeling of sea background. That is,
(13)pζ(ζ)=σnB(n,-α)ζn-1(1+σζ)n-α,σ,n,-α,ζ>0,
where B(·,·) is the beta function. The estimates α^, σ^, and n^ corresponded, respectively, to the parameters α, σ, and n can be easily obtained with the help of numerical calculation based on the method of log-cumulants (MoLC) [12] as
(14)Ψ(n^)-Ψ(-α^)-ln(σ^)=1N∑i=1N[ln(xi)]=c~^1,Ψ(1,n^)+Ψ(1,-α^)=1N∑i=1N[(ln(xi)-c~^1)2],Ψ(2,n^)-Ψ(2,-α^)=1N∑i=1N[(ln(xi)-)3],
where Ψ(·) represents the digamma function (i.e., the logarithmic derivative of the gamma function), Ψ(r,·) is the rth order polygamma function (i.e., the kth order derivative of the digamma function), and {xi},i∈[1,N], is a given sample set.
Given the density shown in (13), its cumulative distribution function (CDF) is written as [11]
(15)Fζ(x)=σnxnnB(n,-α)F12(n-α,n;n+1;-σx),
where F12(·,·;·;·) is the Gauss hypergeometric function. For a given value of the false alarm probability, denoted by Pfa, the corresponding CFAR threshold T for the distribution shown in (13) can be obtained from
(16)1-Pfa=Fζ(T)=σnTnnB(n,-α)F12(n-α,n;n+1;-σT).
Considering that Fζ(T) is strictly monotonously increasing, the threshold T can be accurately calculated via the numerical solution or a simple bisection method.
4. Experimental Results and Analysis
The test dual-pol SAR amplitude data used in this study are a large TerraSAR-X StripMap mode geocoded scene over Nanjing, China, acquired with high-resolution 6 m × 6 m (azimuth × range) and HH-polarization and VV-polarization. Figure 1(a) provides a fake-color image of this scene. The horizontal and vertical axes are the directions of azimuth and range, respectively. Meanwhile, in order to make a visible comparison, Figure 1(b) gives an optical (ground truth) remote sensing photograph of the test site from SPOT5 satellite.
The image of Nanjing: (a) TerraSAR-X image; (b) SPOT5 optical image.
The red rectangle box region shown in Figure 1(a) consisting of several ships and sea clutter is our investigating area. The product image of multilook amplitudes from two polarimetric channels is shown in Figure 2(a). Figure 2(b) shows the fitting results of the distribution in (13) for the product image of the area indicated in the rectangle box in Figure 2(a). The parameters n, α, and σ are estimated to n^=100.5744, α^=-1.7851, and σ^=120.7748, respectively. As shown in Figure 2(b), it is clear that the fitting result agrees well with the theoretical distribution.
The ship chip: (a) the product image; (b) the fitting result.
Furthermore, giving the theoretical false alarm probability Pfa=10-8, the detection results are shown in Figure 3. It is easy to observe that all ships are detected, whilst a false alarm occurs, which proves the effectiveness of CFAR detection method based on PMA detector for ships.
The ship detecting result.
5. Conclusion
Aiming at the adaptive detection of a ship when only high-resolution dual-polarization SAR amplitude data are available, a CFAR detecting method has been proposed in this paper. We first design a novel PMA detector, which can improve the signal-to-clutter ratio (SCR) and make the discrimination of a ship from clutter more easily. Meanwhile, the PMA detector’s statistical model has been described by the well-known 𝒢0 distribution when facing complex sea background. The experiments performed on measured dual-polarization TerraSAR-X images demonstrate the good performance of the proposed CFAR detecting method.
Acknowledgment
The author would like to appreciate the National Natural Science Foundation of China for the support under Grant no. 41171316.
CrispD. J.2004DSTO, Departement Defence, Australian Government, Canberra, Australia, Public Release Document DSTO-RR-0272OliverC. J.QueganS.1998Norwood, Mass, USAArtech HouseSikanetaI. C.2002University of OttawaUlabyF. T.ElachiC.1990Boston, Mass, USAArtech HouseLivingstoneC. E.ThompsonA. A.The moving object detection experiment on RADARSAT-220043033553682-s2.0-4944225493KrylovV. A.MoserG.SerpicoS. B.ZerubiaJ.Supervised high-resolution dual-polarization SAR image classification by finite mixtures and copulas2011535545662-s2.0-7995744031910.1109/JSTSP.2010.2103925ChaneyR. D.BudM. C.NovakL. M.On the performance of polarimetric target detection algorithms199051110152-s2.0-002551690810.1109/62.63157LeeJ. S.HoppelK. W.MangoS. A.MillerA. R.Intensity and phase statistics of multilook polarimetric and interferometric SAR imagery1994325101710282-s2.0-002849669310.1109/36.312890GierullC. H.Statistics of SAR interferograms with application to moving target detection20012001-045Ottawa, CanadaDefense Research Establishment Ottawa, Department of National DefenseGoodmanN. R.Statistical analysis based on a certain multivariate complex gaussian distribution (an introduction)196334152152180FreryA. C.Jacobo-BerllesJ.GambiniJ.MejailM. E.Polarimetric SAR image segmentation with B-splines and a new statistical model20102143193422-s2.0-7795661151510.1007/s11045-010-0113-4NicolasJ. M.Introduction to second kind statistic: application of log-moments and log-cumulants to SAR image law analysis2002193139167