1.

Introduction

Target detection based on synthetic aperture radar (SAR) images has been a hot topic because of its important application in maritime environment monitoring. Since the scattering fields from man-made targets, such as ships and oil platforms, are generally stronger than those from surrounding sea surface, the intensity-based constant false alarm rate (CFAR) algorithms have been widely established to distinguish targets from sea surface.¹2.3.4.5.6.7.8.9.10.11.^–12 However, the detection algorithms based on single-channel SAR data are sometimes ineffective due to various factors, such as speckle noise, the ambiguities of targets, and other strong noise induced by some uncertain reasons. Recently, with the availability of polarimetric SAR sensors, the polarimetric scattering behavior of targets can be described by the scattering matrix, through which features relating to both the intensity and the relative phase can be extracted. Recently, many polarimetric detection algorithms using polarimetric information provided by fully polarimetric SAR are successfully exploited for detecting targets, especially marine targets.⁷8.9.10.11.^–12 In Ref. 12, Novak proposed the polarimetric whitening filter method for target detection. On the other hand, the coherent target decomposition (CTD) theory can also be used to observe ships because the scattering mechanisms of targets are generally different from that of surrounding sea surface. The results in Ref. 13 showed that even in low resolution, SAR imageries marine targets can also be distinguished from sea clutter by CTD theory. Generally, the backscattered field from sea surface is dominated by Bragg scattering with low entropy value. However, owing to the complex metallic structures of man-made targets, the scattering field from a man-made target always contains complex scattering mechanisms with higher entropy value. In Refs. 14 and 15, based on CTD theory, polarimetric SAR images were used to detect ships exploring entropy values. Note that the effectiveness of the entropy-based detector decreases with increasing wind speed because the scattering field from breaking waves also contains complex scattering mechanisms with higher entropy value. Although the methods mentioned above have been widely used to detect targets in SAR images, the strong noises induced by ambiguities and so on, which are often mistaken as real targets and cause false alarms, cannot be effectively removed.

In Ref. 16, the eigenvalues of the coherency scattering matrix are used to distinguish ambiguities from real targets. However, the method in Ref. 16 will no longer perform well in the case of strong depolarized scattering. In Ref. 17, Liu and Gierull found that in $\mathrm{HV}+\mathrm{VH}$ image, the scattering fields of azimuth ambiguities can be suppressed effectively because the HV channel data of the first-order azimuth ambiguities (FOAAs) are approximately equal to VH channel in magnitude but shifted about $\pm \pi$ in phase. In Ref. 18, based on the intensity of $\mathrm{HV}+\mathrm{VH}$ image (it is called as ambiguity free image ${\mathrm{HV}}_{\text{free}}$ by Velotto), a new parameter $\gamma$, which is a simple combination of slope and departure from the uniform distribution of the phase parameters in GK distribution, is proposed for target detection (refer to Ref. 18 for detail). Velotto et al. pointed out that $\gamma$ is always smaller than 10 for real targets but larger than 10 for sea surface and ambiguities. The methods proposed in Refs. 17 and 18 are applicable to detect targets effectively in high-quality SAR images, by this I mean the SAR imageries without exotic noises. Nevertheless, besides ambiguities, many SAR images have also been corrupted by other strong noise induced by uncertain reasons. For the case of the poor-quality SAR images, methods of Liu and Velotto will be no longer in validity.

This paper mainly focuses on detection of azimuth ambiguities in noise severely corrupted fully polarimetric SAR images. We define reciprocal of the mean square value of the relative phase (RMSRP) as a new indicator of marine targets considering the different behaviors of phase difference over real targets and FOAAs. Based on the central limit theorem, we find that the mean value of RMSRP estimated through a $11\times 11$ spatially moving window satisfies the Gaussian distribution well. The theoretical distribution of RMSRP is derived and in terms of this, a RMSRP-based CFAR algorithm is established for marine target detection. The experiment results using Radarsat-2 imageries acquired in fine quad-polarization mode are promising, demonstrating that our detection algorithm is not only robust for removal of azimuth ambiguities in SAR images but also is able to suppress noise induced by unknown reasons. This paper is organized as follows. In Sec. 2, the statistic properties of the relative phase and RMSRP are discussed, and the probability density function (PDF) of RMSRP for surrounding sea surface is derived. Section 3 addresses the RMSRP-based CFAR algorithm. Experimental results are presented in Sec. 4, whereas conclusions are given in Sec. 5.

2.

Definition of the New Parameter RMSRP

In this work, the relative phase is estimated by

For convenience, the new parameter RMSRP is represented by symbol $\mathrm{\Theta }$, and it is defined as

2.2.

For Sea Surface and the Strong Noise

In SAR images, for sea surface and the strong noise, the relative phase between two cross-polarized channels has been widely spread in angle region $[-\pi ,\pi ]$ due to decorrelation effect. The histograms of $\phi$ evaluated by SLC quad-pol RS-2 SAR data are illustrated in Fig. 1. It is clearly shown that the relative phase is equal to 0 rad for real target, and the PDF is very narrow. Just as discussed in Sec. 2.1, the PDF for FOAAs is distributed around $\pm \pi$, but it has been spread because the field intensities of ambiguities are generally much lower than real target and can be easily affected by sea clutter. Figure 1 also shows that the relative phases for sea surface and strong noise are widely spread in angle region $[-\pi ,\pi ]$ due to decorrelation effect.

Fig. 1

Phase differences between cross-polarized channels.

In order to detect target by RMSRP-based CFAR algorithm, the “prior probability” of $\mathrm{\Theta }$ should be obtained first. In Refs. 20 and 21, the PDF of the relative phase $\phi$ has been derived by Lee as

(8)

$$P(\phi )=\frac{\mathrm{\Gamma }(L+1/2){(1-{|}^2)}^L\beta }{2\sqrt{\pi }\mathrm{\Gamma }(L)(1-{\beta }^2)}+\frac{{(1-{|\rho |}^2)}^L}{2\pi }F(L,1;1/2;{\beta }^2),$$

where $F(L,1;1/2;{\beta }^2)$ is a Gauss hypergeometric function, $\beta =|\rho |\cos (\phi -\langle \phi \rangle )$, the correlation coefficient $|\rho |$ and the mean phase $\langle \phi \rangle$ are related to the measured data by

(9)

$$|\rho |\exp (j\langle \phi \rangle )=\frac{\langle S_{\mathrm{HV}}{S}_{\mathrm{VH}}^{*}\rangle }{\sqrt{\langle {|S_{\mathrm{HV}}|}^2\rangle \langle {|S_{\mathrm{VH}}|}^2\rangle }}.$$

Using Eq. (8), the PDF of ${\phi }^2$ is derived as

(10)

$$P({\phi }^2)=\frac{\mathrm{\Gamma }(L+1/2)(1-{|\rho |}^2){}^L\beta }{4\phi \sqrt{\pi }\mathrm{\Gamma }(L)(1-{\beta }^2)}+\frac{{(1-{|\rho |}^2)}^L}{2\pi }F(L,1;1/2;{\beta }^2).$$

Based on the central limit theorem, if the dimension number of a moving window is large enough, the PDF of $\psi$ (i.e., the mean value of ${\phi }^2$ within the moving window) can be described approximately by a Gaussian distribution

(11)

$$P(\psi )=\frac{1}{\sqrt{2\pi {\sigma }_{\psi }^2}}\exp [-\frac{{(\psi -{\mu }_{\psi })}^2}{2{\sigma }_{\psi }^2}],$$

where ${\mu }_{\psi }$ and ${\sigma }_{\psi }^2$ denote the mean value and the mean square deviation of $\psi$, respectively. Then, using the identity equation $P(\psi )|\mathrm{d}\psi |=P(\mathrm{\Theta })|\mathrm{d}\mathrm{\Theta }|$, the PDF of $\mathrm{\Theta }$ is obtained

(12)

$$P(\mathrm{\Theta })=\frac{1}{{\mathrm{\Theta }}^2\sqrt{2\pi {\sigma }_{\psi }^2}}\exp [-\frac{{(1/\mathrm{\Theta }-{\mu }_{\psi })}^2}{2{\sigma }_{\psi }^2}].$$

Using Eqs. (10) and (12), the mean values of $\mathrm{\Theta }$ for different look numbers are evaluated, and the curves are plot in Fig. 2. Just as shown in this figure, the value of $\mathrm{\Theta }$ increases with the increase of look number and correlation coefficient. And the effect of look number is more remarkable when the correlation coefficient is larger. This property implies that multilook processing would improve the ability of $\mathrm{\Theta }$ for target detection because the correlation coefficients for real targets are always larger than those for surrounding sea surface. However, it should be pointed out that the multilook processing would reduce the spatial resolution. On the other hand, to ensure the applicability of the central limit theorem, the dimension of the moving window, in which the mean value of ${\phi }^2$ is evaluated, is set to be $11\times 11\text{pixels}$. In this work, a data set that contains six Radarsat-2 SAR images acquired in fine quad-polarization mode are processed to show the effectiveness of the theoretical PDF. The detailed descriptions of the images are given in Table 1. The comparisons in Fig. 3 illustrate that the theoretical PDFs match the histograms well.

Fig. 2

Plots of $\mathrm{\Theta }$ versus the correlation coefficient for different look numbers.

Table 1

The data set used in this work.

Image’s name	Acquired time (UTC)	Center locations	Inc. angle Nr-Fr (deg)
Image RS_I	July 30, 2008 09:17:22	129°07′16′′E 34°50′51′′N	28.04 to 29.91
Image RS_II	September 18, 2009 10:49:30	109°24′03′′E 18°06′16′′N	32.35 to 34.01
Image RS_III	May 15, 2010 11:56:37	88°18′34′′W 28°33′06′′N	29.17 to 30.92
Image RS_IV	August 19, 2011 10:52:36	120°09′25′′E 38°20′09′′N	36.47 to 38.00
Image RS_V	February 15, 2013 10:13:39	118°14′41′′E 39°01′04′′N	46.02 to 47.25
Image RS_VI	July 13, 2010 09:39:37	122°34′54′′E 31°11′43′′N	48.31 to 49.47
Image RS_VII	August 18, 2012 10:12:01	116°43′27′′E 20°42′55′′N	29.19 to 30.94

Fig. 3

Comparisons of the theoretical probability density functions (PDFs) with the histograms of $\mathrm{\Theta }$.

3.

$\mathrm{\Theta }$-based CFAR Target Detector

Here, a $\mathrm{\Theta }$-based CFAR algorithm is established for marine target detection. The detection threshold $\xi$ is determined by the following integration:

Based on the analysis above, the following rule is used to detect targets, i.e.,

4.

Experimental Results

In this section, three meaningful Radarsat-2 SAR imageries (namely Image RS_V, Image RS_VI and Image RS_VII in Table 1), in which targets, azimuth ambiguities, and strong noise are simultaneously present, have been selected to demonstrate the effectiveness of the $\mathrm{\Theta }$-based target detector. In the following discussions, the false alarm rate $\eta$ is set to be ${10}^{-5}$.

First, the $\mathrm{\Theta }$-based target detector is tested to remove azimuth ambiguities. Figure 4 shows the SAR image of the Dongsha Island acquired at 10:12 (UTC) on August 8, 2012 (namley Image RS_VII in Table 1). From Fig. 4, we can find that the first-order azimuth ambiguity of the Dongsha Island can be seen in $|\mathrm{HH}|$ and $|\mathrm{HV}|$ images. In $|\mathrm{HV}+\mathrm{VH}|$ image (i.e., the so called ${\mathrm{HV}}_{\text{free}}$ image in Ref. 18), the intensity of the azimuth ambiguity has been suppressed due to its special phase difference. From Fig. 4(d), we can find that the value of $\mathrm{\Theta }$ (i.e., RMSRP) for target/azimuth ambiguity is higher/lower than the surrounding sea surface. Thus, the parameter $\mathrm{\Theta }$ can also be used to remove the azimuth ambiguities. Figure 4(e) shows the detection results by the $\mathrm{\Theta }$-based target detector by using threshold 0.77, and it is obvious that this proposed algorithm is good to detect real targets.

Fig. 4

An example of SAR azimuth ambiguity removal using RMSRP. The bright target is Dongsha Island in East China Sea. (a) $|\mathrm{HH}|$, (b) $|\mathrm{HV}|$, (c) $|\mathrm{HV}+\mathrm{VH}|$, (d) RMSRP, (e) Detection result.

By comparison in Fig. 4, it is evident that the azimuth ambiguities can be suppressed effectively within ${\mathrm{HV}}_{\text{free}}$ or RMSRP image. However, in the case of strong noise corrupted SAR images, the ${\mathrm{HV}}_{\text{free}}$-based algorithm will no longer perform well. This behavior is shown in Fig. 5(a), and the SAR image was acquired at 10:13 (UTC) on February 15, 2013 (named Image RS_ V in Table 1) in the Bohai sea. The center longitude and latitude coordinates are about 118°14′41′′E and 39°01′04′′N. The yellow rectangle in Fig. 5(a) indicates the blown-up area, as shown in Figs. 5(b)–5(f). The red circles in Fig. 5(c) highlight the areas with strong noise induced by unknown reasons. Unlike azimuth ambiguities, the strong noise cannot be suppressed effectively in ${\mathrm{HV}}_{\text{free}}$ image [see Fig. 5(d)], while reduced over RMSRP image as given in Fig. 5(e). The reason is attributed to the fact that PDF of the relative phase in the noise-polluted region is always broader than that for real targets, making RMSRP over noise-region much smaller than that of real targets. Figure 5(f) shows the detection results by $\mathrm{\Theta }$-based target detector with threshold 1.49, as we expected that this algorithm is able to suppress both azimuth ambiguities and strong noise. To validate the detection results, the detailed information of targets achieved by Electronic Navigation Chart (ENC) and Automatic Identification System (AIS) are given in Table 2. This ground truth is consistent with the detection results, further supporting our proposed detector. The largest ship in Fig. 5 is T18, with 294 m in length and 32 m in width. There do exist two odd targets circled by red rectangle in Fig. 5(f), and no AIS information indicates these targets.

Fig. 5

An example of target detection in noise-corrupted SAR image (RS_ V in Table 1). (a) $|\mathrm{HH}|$, (b) $|\mathrm{HH}|$, (c) $|\mathrm{HV}|$, (d) $|\mathrm{HV}+\mathrm{VH}|$, (e) RMSRP, (f) Detection result.

Table 2

The information of the targets in Fig. 5.

Note: The information for T27 and T28 are achieved from ENC (electronic navigation chart) C1311700. No AIS information are available for T10 and T26. The information for the other targets is all obtained by AIS (automatic identification system).

Another typical example of target detection in a noise-corrupted SAR imagery is presented in Fig. 6. The SAR image in Fig. 6 was acquired at 09:39 (UTC) on July 13, 2010, in the East China Sea. And the center longitude and latitude coordinates are 122°34′54′′E and 31°11′43′′N, respectively. Besides real targets, some filaments could be seen in Figs. 6(a)–6(c), and we have no idea about the origin of these odd noise. Figure 6(f) demonstrates promising detection result with threshold of 0.67. It is exciting to find that the strong noise has been suppressed in the RMSRP image, highlighting the visually inspected ships. This imagery provides further evidence to the applicability of our proposed algorithm.

Fig. 6

Another example of target detection in noise-corrupted SAR image (named RS_ VI in Table 1). (a) $|\mathrm{HH}|$, (b) $|\mathrm{HH}|$, (c) $|\mathrm{HV}|$, (d) $|\mathrm{HV}+\mathrm{VH}|$, (e) RMSRP, (f) Detection result.

To present a quantitative description about the reduction degree of noise, we introduce the standard deviation of $|\mathrm{HH}|$, $|\mathrm{HV}|$, $|\mathrm{HV}+\mathrm{VH}|$, and RMSRP over the noise-corrupted areas to indicate this improvement. We select two areas in Fig. 5(d) marked as AR1, AR2, and one area AR3 in Fig. 6(b) to demonstrate the improvement of RMSRP in suppressing noise. The detailed estimates for different parameters are given in Table 3, and it is evident that the noise is primarily reduced by RMSRP while others do not exhibit such capability. It also should be noted that the standard deviation of $|\mathrm{HV}+\mathrm{VH}|$ is always the largest value, by this I mean that noise presents the highest variability in $|\mathrm{HV}+\mathrm{VH}|$ images.

Table 3

Standard deviation over noise-corrupted areas.

5.

Conclusions

In this paper, the RMSRP between cross-polarized channels is proposed as a new promising indicator for SAR target detection at sea. Meanwhile, an analytical form of the PDF for RMSRP was derived, and the comparisons between the histograms of the sea clutter and theoretical PDF show the correctness of the statistic model. Based on the statistic model, a RMSRP-based CFAR algorithm is developed to detect targets in noise corrupted SAR images. From the tests over quad-pol RS-2 SAR images, we find that not only the azimuth ambiguities but also strong noise can be reduced successfully by the RMSRP-based algorithm. Standard deviation is introduced to indicate the quantitative improvement in suppressing noise. The experiment results certify that RMSRP is a promising indicator for Pol-SAR target detection purpose, particularly when pol-SAR images have been corrupted by strong noise and azimuth ambiguities. This proposed algorithm would be helpful in detecting targets in SAR imageries and further applying to coastal management.