# Target detection in synthetic aperture radar imagery: a state-of-the-art survey

**Khalid El-Darymli**

C-CORE, Captain Robert A. Bartlett Building, Morrissey Road, St. John’s, Newfoundland, A1C 3X5, Canada

Memorial University, St. John’s, Newfoundland, A1B 3X5, Canada

**Peter McGuire**

C-CORE, Captain Robert A. Bartlett Building, Morrissey Road, St. John’s, Newfoundland, A1C 3X5, Canada

**Desmond Power**

C-CORE, Captain Robert A. Bartlett Building, Morrissey Road, St. John’s, Newfoundland, A1C 3X5, Canada

**Cecilia Moloney**

Memorial University, St. John’s, Newfoundland, A1B 3X5, Canada

*J. Appl. Remote Sens*. 7(1), 071598 (Mar 18, 2013). doi:10.1117/1.JRS.7.071598

#### Open Access

**This article was corrected** | View correction

**Abstract.**
Target detection is the front-end stage in any automatic target recognition system for synthetic aperture radar (SAR) imagery (SAR-ATR). The efficacy of the detector directly impacts the succeeding stages in the SAR-ATR processing chain. There are numerous methods reported in the literature for implementing the detector. We offer an umbrella under which the various research activities in the field are broadly probed and taxonomized. First, a taxonomy for the various detection methods is proposed. Second, the underlying assumptions for different implementation strategies are overviewed. Third, a tabular comparison between careful selections of representative examples is introduced. Finally, a novel discussion is presented, wherein the issues covered include suitability of SAR data models, understanding the multiplicative SAR data models, and two unique perspectives on constant false alarm rate (CFAR) detection: signal processing and pattern recognition. From a signal processing perspective, CFAR is shown to be a finite impulse response band-pass filter. From a statistical pattern recognition perspective, CFAR is shown to be a suboptimal one-class classifier: a Euclidean distance classifier and a quadratic discriminant with a missing term for one-parameter and two-parameter CFAR, respectively. We make a contribution toward enabling an objective design and implementation for target detection in SAR imagery.

Synthetic aperture radar (SAR) offers distinctive active remote sensing capabilities for both military and civilian applications. Target, clutter, and noise are three terms of military origins associated with automatic target recognition (ATR), and their definition depends on the application of interest. In the case of SAR imagery, target refers to the object(s) of interest in the imaged scene. Clutter refers to manmade (building, vehicles, etc.) and/or natural objects (trees, topological features, etc.) that tend to dominate the imaged scene. Noise refers to imperfections in the SAR image which are a result of electronic noise in the SAR sensor, as well as computational inaccuracies introduced by the SAR signal processor. The general structure of an end-to-end ATR system for SAR imagery (SAR-ATR), as reported in the literature, is depicted in Fig. 1. To account for the prohibitive amounts of processing pertaining to the input SAR imagery, the strategy is to divide and conquer. Accordingly, the SAR-ATR processing is split into three distinctive stages: detector (also known as prescreener), low-level classifier (LLC, also known as discriminator), and high-level classifier (HLC).^{1}^{–}^{9} The first two stages together are commonly known as the focus-of-attention module. While this is the most common structure reported in the literature, it should be highlighted that (theoretically) there is no restriction on the number of stages.

As depicted in Fig. 1, the input SAR image creates an extremely high computational load due to its high resolution and/or the presence of various clutter types and objects. As the SAR data progresses throughout the SAR-ATR processing chain, its load is reduced. The HLC stage deals with SAR data that has relatively lower computational load. To the contrary, the computational complexity of the SAR-ATR chain increases as the SAR data progresses from the *front-end* stage toward the *back-end* stage.

Detection is the front-end stage in any SAR-ATR processing chain. The detector interfaces with the input SAR image to identify all regions of interest (ROIs), thus ROIs can be passed-in to the LLC stage for further analysis. One may think of the detector as a dimensionality reduction scheme that properly reduces the dimensionality of the SAR data. The detector should be designed to balance the tradeoff between computational complexity, detection efficacy, and outlier rejection. On the one hand, it is required that the detector is relatively computationally simple, thus it can operate in real-time or near-real-time. On the other hand, it is required that the detector enjoys a low probability of false alarm (PFA), and a high probability of detection (PD). Indeed, these often conflicting factors distinguish one detector from another.

There are numerous strategies for implementing the detector. This is evident in the overwhelming number of research articles published on the topic in the open literature. Different researchers tend to approach the topic from various perspectives. This makes it even more challenging and time consuming to relate the various research findings and to grasp the relationship between these various approaches. This shows a dire need for a survey that offers an umbrella under which various research activities can be broadly probed and taxonomized. This is precisely the goal of this paper.

In this paper, we restrict our attention to *single-channel* SAR imagery (i.e., single polarization). This is because the development of our survey is motivated by our endeavor to develop SAR-ATR algorithms for Spotlight mode Radarsat-2 data, which is a single channel.^{10} However, many of the topics described in this survey are either applicable or extendable to multichannel SAR imagery. For readers interested in multichannel SAR image processing, please refer to pertinent references.^{11}^{–}^{17}

The remainder of this paper is organized as follows. In Sec. 2, taxonomy of the detection methods is introduced. Primarily, three taxa are proposed: single-feature-based, multifeature-based, and expert-system-oriented. In Sec. 3, the various methods for implementing the detection module are comprehensively surveyed. First, a classification methodology for the various strategies under each taxon is introduced. Next, a description of the different sub-taxa is elaborated. The description commences with the Bayesian approach being the optimal approach, and it ends with the various detection strategies that fall under the expert-system-oriented taxon. Representative examples pertinent to SAR imagery are carefully chosen, and relevant comments are pinpointed throughout this section. The issues approached include proper choice of the size of the constant false alarm rate (CFAR) sliding window, CFAR as a suboptimal one-class classifier, CFAR loss, and cell averaging CFAR (CA-CFAR) as the baseline algorithm for comparison with other CFAR techniques. In Sec. 4, a compact tabular comparison between the surveyed methods is offered. In Sec. 5, a novel discussion is presented. The issues covered under the discussion include suitability of SAR data models, understanding the multiplicative SAR data models, and two unique perspectives on CFAR detection: a signal processing perspective and a statistical pattern recognition perspective. In Sec. 6, the paper is concluded.

The detection module takes the entire SAR image and identifies the ROIs. Ultimately, the detected regions in the image are passed-in to the next stage in the SAR-ATR chain for further analysis. The goodness of any detection module is typically judged based upon three aspects of significance: computational complexity, PD, and false alarm rate (also known as PFA). The detection module should enjoy a low computational complexity such that it operates in real-time or near-real-time. This is in contrast to the succeeding stages in the SAR-ATR chain, which are relatively more computationally expensive. Further, a good detection module should provide a means to refine detections, reduce clutter false alarms, and pass ROIs; thus the detection method enjoys a reasonable PFA and acceptable PD.

We broadly taxonomize the detection algorithms reported in the open literature into three major taxa: single-feature-based, multifeature-based, and expert-system-oriented. This taxonomy is depicted in Fig. 2.

The single-feature-based taxon bases the detection in the SAR image on a single feature; typically the brightness in the pixel intensity commonly known as the radar cross-section (RCS). Various methods in the literature fall under this taxon. The single-feature-based approach is placed at the base of the pyramid in Fig. 2 because it is the most common and widely used in the literature. Further, the single-feature-based approach is the building block for the other two taxa.

The multifeature-based taxon bases the detection decision on a fusion of two or more features extracted from the input SAR image. Besides RCS, additional features that can be inferred and fused include multiresolution RCS and fractal dimension. Obviously, this taxon builds on the previous taxon and is expected to provide relatively improved detection performance, along with fewer false alarms. Multiple methods in the literature fall under this taxon.

Finally, the expert-system-oriented taxon is the most sophisticated. It extends the two aforementioned taxa and utilizes a multistage (two or more stages) artificial intelligence (AI) approach that bases the detection process in the SAR image on exploitation of prior knowledge about the imaged scene, clutter, and/or target(s). Prior knowledge is exploited through various means such as image segmentation, scene maps, previously gathered data, etc.

As the sophistication of the detection taxon increases, the complexity-performance tradeoff arises. Caution should be taken when opting for a certain approach in order to balance this tradeoff carefully.

Based on the aforementioned taxonomy, we broadly classify the various detection schemes and relevant methods reported in the literature in Fig. 3.

Primarily, under the single-feature-based taxon, the sliding window CFAR (CFAR-based) sub-taxon is the most popular. The various CFAR methods can be viewed through three perspectives. First, based on the specifications of the sliding window, there is fixed-size versus adaptive, as well as rectangle-shaped versus nonrectangle-shaped. Second, based on the strategy used to implement the CFAR technique, there are various strategies, including cell-averaging CFAR (CA-CFAR), smallest of CA-CFAR (SOCA-CFAR), greatest of CA-CFAR (GOCA-CFAR), and order statistics CFAR (OS-CFAR). Third, based on the method used to estimate the threshold (for a desired PFA) in the boundary ring and/or the approach utilized for estimating the target signature (for a desired PD), two subclasses emerge: parametric and nonparametric. Under the parametric subclass, two approaches are recognized: only background modeling and background and target modeling. A choice of the parametric model that best represents the SAR data in use has to be made among various parametric models. Unlike the parametric approach, the nonparametric approach does not assume any form for the background/target model(s). Rather, it directly infers an approximate model from the training data. One such method to perform the model inference is the kernel-density-estimation (KDE) method.

Less popular non-CFAR-based methods, such as those that rely on a coherence image, represent the other sub-taxon of single-feature-based methods. The single-feature-based taxon has the limitation that it bases the detection decision solely on RCS, and thus it can become overwhelmed in regions in the SAR image where there is heterogeneous clutter and/or a high density of targets.

Methods under the multifeature-based taxon try to circumvent this drawback by basing the detection decision on a fusion of two or more features. Obviously, this taxon can utilize a suitable method among those presented under the single-feature-based taxon and incorporate additional features besides RCS, such as multiresolution RCS analysis, fractal dimension, etc. Multiresolution methods can be either space-scale-based or space-frequency-based. Prime examples of methods that utilize space-scale features are those based on the wavelet transform, including the discrete wavelet transform (DWT), and the continuous wavelet transform (CWT). Prime examples of methods that utilize space-frequency features include linear space-frequency methods such as the Gabor transform and the S-transform, along with bilinear (also known as quadratic) space-frequency methods such as Cohen’s class distributions (Wigner distribution, Wigner-Ville distribution, pseudo-Wigner-Ville distribution, etc.).

Finally, a more robust taxon is the expert-system-oriented approach, which incorporates intelligence into the process to guide the decision making. In its simplest form, detection decisions can be guided by a structure map of the imaged scene generated from properly segmenting the SAR image. Further, methods of AI can be appropriately integrated to achieve near-optimal context utilization.

Next we review the various methods introduced above under each taxon. Further, representative examples pertaining to SAR imagery under each method are carefully chosen and presented.

Single-feature-based detection algorithms base their search for target detection in the SAR image on a single feature. CFAR is the most popular single-feature-based detection algorithm. Despite the many variations of CFAR under this category, they are considered single-feature-based because they base the search for ROIs on RCS alone. Indeed, as it is evident from the numerous works published in the literature, CA-CFAR is the baseline approach for target detection in SAR imagery. To understand the limitations of the single-feature-based-CFAR approach, it is important to review its underlying assumptions.

An optimal detector (theoretically speaking) should utilize the Bayesian approach, which for a zero-one cost, reduces to the maximum *a posteriori* (MAP) criterion^{18} as

This is simply a binary (i.e., two-class; dichotomizer) classification problem where $x$ is a feature vector that represents the pixel values, typically obtained from the boundary ring in a sliding window with suitable guard cells centered on the ROI. This window is typically called a CFAR stencil and is depicted in Fig. 4 for a size of $9\xd79\u2009\u2009pixels$. The boundary ring is shown in bright green with suitable pixel labels.

Note that the choice of this stencil size here is for illustration purposes only. Proper choice of the stencil size will depend on the type of SAR image and the target size. More precisely, the size (and shape) of the guard ring will be governed by the geometrical size of the target.^{19} However, it should be noted that choosing the stencil dimension using only the prior knowledge of target size yields a detection loss (i.e., CFAR loss that leads to a suboptimal performance), because the backscatter of the target in SAR imagery is dependent on the operating conditions and is weakly linked with the target’s geometric shape.^{20} Despite these challenges, it is recommended^{21}^{,}^{22} that the target window size (i.e., pixels under test, PUTs) should be about the size of the smallest object that one wishes to detect, the guard ring window size should be about the size of the largest object, and the boundary ring window size should be large enough to estimate the local clutter statistics accurately.

Posterior probabilities can be expressed by the Bayes rule as

Accordingly, the MAP criterion can be expressed as a likelihood ratio test (LRT) as

Equation (4) is the main formula motivating the design for the CFAR algorithm and its variants. Indeed, many popular CFAR algorithms in the literature assume that only the background-class is characterized without characterizing the target-class, and thus they adopt a suboptimal anomaly detection (AD). This converts the optimal binary classification problem into a one-class classification problem^{23} as

^{24}

^{,}

^{25}as shown in Eq. (5).

In the former approach (i.e., the AD approach), the CFAR algorithm assigns PUTs to the background if it finds that the PUTs are consistent with the background distribution; otherwise, the PUTs are labeled as detected. For a desired PFA, the scaling factor $\alpha $ is adaptively estimated throughout the image from the boundary ring in the sliding window focused on the ROI as

Be it a fixed-size or adaptive sliding window, the various CFAR methods can be viewed through two perspectives intermingled. The first is based on the method used to estimate the threshold scaling factor (for a desired PFA) in the boundary ring and/or the approach utilized for estimating the target signature (for a desired PD). There are two strategies: parametric CFAR and nonparametric CFAR. The second is based on the method used to implement the CFAR technique. There are various strategies, including CA-CFAR, SOCA-CFAR, GOCA-CFAR, and OS-CFAR. Thus, any CFAR detector can be viewed as a combination of these two perspectives: one strategy pertinent to the estimation of the threshold scaling factor and one strategy for implementing the CFAR technique. To visualize the interrelation between the various methods and strategies, refer to Fig. 3.

Under this section, we review parametric CFAR and nonparametric CFAR. Under parametric CFAR, we review one-parameter CFAR and two-parameter CFAR. Under one-parameter CFAR, we review various implementation strategies, including CA-CFAR, SOCA-CFAR, GOCA-CFAR, and OS-CFAR. Under two-parameter CFAR, we discuss the most common implementation strategy. We then briefly discuss CFAR loss. Then an interesting remark that addresses an important issue pertaining to CFAR usage is presented. Finally, the topic of nonparametric CFAR is briefly approached.

Parametric CFAR methods can be classified into two classes: methods based only on background modeling (i.e., AD) and methods based on both background and target modeling. We briefly review these two methods here.

All the parametric CFAR algorithms that perform AD (i.e., only models the background clutter) have one thing in common. They all assume that the background clutter can be roughly modeled by a certain probability distribution, i.e., $p(x|\omega B)$. Then, to perform the CFAR detection, they estimate the model distribution parameters from the boundary ring in the CFAR stencil. This PDF model is used to estimate the threshold scaling factor $\alpha $, for a desired CFAR (i.e., PFA), as the focusing window is systematically slid over the SAR image. However, variant classes of CFAR algorithms primarily differ in two aspects.

First, there is the probability distribution chosen, i.e., $p(x|\omega B)$, to model the background clutter. For example, some CFAR algorithms assume a homogeneous clutter and model the background clutter with an exponential distribution (i.e., for SAR image in the power-domain; magnitude-squared), or a Rayleigh distribution (i.e., for SAR image in the magnitude-domain). This class of distribution models is characterized by one-parameter (i.e., the mean) and thus is referred to in the literature as a one-parameter CFAR. Other CFAR algorithms model the background clutter in the SAR image with the more realistic but more complex Weibull distribution,^{26}^{,}^{27}$K$-distribution,^{28} alpha-stable distribution,^{29}^{,}^{30} or beta-prime ($\beta \u2032$) distribution,^{31} among other models. This class of distribution models is characterized by two parameters (mean and variance, scale and shape parameters, etc.), and because of this, the CFAR algorithm is referred to in the literature as a two-parameter CFAR.

Second, there is the method used to estimate the model parameters pertaining to the detection Threshold from the boundary ring. For example, there is CA-CFAR, GOCA-CFAR, SOCA-CFAR, OS-CFAR,^{32}^{–}^{34} and best linear unbiased estimator CFAR (BLUE-CFAR),^{26}^{,}^{27} among others.

The CFAR algorithms that model both the background and the target perform procedures similar to those mentioned above. However, besides estimating the background model parameters in the focusing window, they also estimate the target model parameters. Thus, the detection threshold in the sliding window is determined based on the NP criterion as shown in Eq. (4). Gan and Wang^{35} and Rong-Bing and Jian-Guo^{36} offer examples on this approach.

Under this section, we briefly review various one-parameter CFAR implementation strategies, including CA-CFAR, SOCA-CFAR, GOCA-CFAR, and OS-CFAR. Obviously, these same strategies are also utilized for implementing two-parameter CFAR. Accordingly, a proper understanding of these implementation strategies paves the way for better understanding two-parameter CFAR.

CA-CFAR was the first CFAR test proposed in 1968 by Fin and Johnson.^{37}^{–}^{39} The adaptive Threshold is comprised of two parts. One is estimated from the boundary ring called $Z$, and the other is found from the corresponding PDF distribution for a desired PFA. This explains why $\alpha $ is referred to as the threshold scaling factor. Thus, Threshold is given by

Assuming a power image, and assuming the clutter is independent and identically distributed (iid), and that the PDF of a pixel $xi$ is exponential distributed, whereby both the $I$ and $Q$ channels with power $Z/2$ in each channel (total $Z$), as follows:

^{40}by

^{40}

In reference to Fig. 4, CA-CFAR computes the arithmetic average of the pixels in the boundary ring and then compares it with the PUT. The decision to conclude either that a PUT is detected (i.e., target-class) or not detected (i.e., background-class) is contingent upon the threshold scaling factor $\alpha $, as illustrated below. The SAR image is in the (noncomplex) power domain:

Note that if a log-detector CFAR is used (i.e., the SAR image is assumed to be in the log-domain), one simply takes the logarithm for both sides in Eq. (15). Accordingly, the log-detector one-parameter CFAR is governed by

Additionally, note that it is possible to consider more than one pixel for the PUT in the inequalities above. However, in this case, $XPUT$ should be replaced with the ML estimate of the arithmetic average pertaining to the $M$ PUTs given by

^{40}(also known as CFAR loss).

SOCA-CFAR and GOCA-CFAR are variants of CA-CFAR that split the boundary ring in the sliding stencil into separate leading and lagging windows (indeed, there are four such windows, one on each side of the PUT) that are used to estimate separate statistics.^{41} This is depicted in Fig. 5.

In SOCA-CFAR, the smallest of the four mean estimates is used to perform the test, while in GOCA-CFAR, the greatest of the four mean estimates is used. Assume there are a total of $N$ pixels in each window. Thus, there are four mean estimates:

Accordingly, for SOCA-CFAR, the detection decision for power-domain (i.e., square-law) detection is

^{40}.

SOCA-CFAR is designed to handle strong clutter returns in the boundary ring, but it is susceptible to clutter edges. On the other hand, GOCA-CFAR is capable of performing better than CA-CFAR and SOCA-CFAR at clutter edges, but its performance degrades when strong returns appear in the boundary ring. Further, compared to CA-CFAR, both SOCA and GOCA suffer from additional CFAR loss, due to the consideration of only a part of the boundary ring.

OS-CFAR was first proposed by Rohling^{42} to counter multiple targets in the CFAR stencil. OS-CFAR rank orders the pixels in the boundary ring according to their value. Assuming that there are N pixels in the boundary ring of Fig. 4, OS-CFAR orders them in an ascending order:

^{42}$Q=3/4$ is found to work best. However, in later work,

^{43}

^{,}

^{44}$Q=4/5$ is found to give better results. Obviously, the best value of $Q$ is dependent on the type of the SAR data used. Thus, bearing these values in mind, it is desirable to empirically check the value of $Q$ that best suits the data in use. For derivation of the threshold scaling factor $\alpha $ for a desired PFA under the Gaussian assumption, the reader is referred to page 372 in Richards.

^{40}

In heterogeneous/nonhomogeneous clutter backgrounds, and for contiguous targets, OS-CFAR is known to outperform CA-CFAR.^{38}^{,}^{42}^{,}^{45}^{,}^{46} However, the performance of OS-CFAR degrades during clutter transitions.^{45} This motivated researchers to develop this method further to handle various scenarios. Although only tested on 1-D radar data (i.e., range profiles), a variant called switched order statistics (SWOS) CFAR is designed for detecting targets in nonhomogenous clutter and/or multiple interfering target scenarios.^{47} This algorithm builds on a relevant method known as selection and estimation^{48} (SE) and is able to determine whether the cells in the boundary ring belong to homogeneous or nonhomogeneous clutter and thus adaptively adjust the detection threshold for a desired PFA. SWOS-CFAR is shown to outperform standard OS-CFAR^{49} and switching CFAR (S-CFAR). A generalized version of S-CFAR termed GS-CFAR is proposed, wherein it is shown that GS-CFAR yields some improvement in the detection probability in the case of interfering targets and clutter edges.^{50} Many additional variations of CFAR exist in ^{40}.

Gandhi and Kassam^{51} show that the performance of the CA-CFAR processor approaches that of the NP detector in that the limit of the number of pixels in the boundary ring approaches infinity under the condition that homogeneity is maintained in the boundary ring. Thus, CA-CFAR achieves asymptotically optimal performance under these (theoretical) assumptions. Accordingly, CA-CFAR can be used as the baseline algorithm for comparison with other CFAR techniques.

We summarize the four basic (square-law) CFAR detectors (i.e., CA-CFAR, SOCA-CFAR, GOCA-CFAR, and OS-CFAR) in Table 1.

Method | Formula | Advantages | Disadvantages |
---|---|---|---|

CA-CFAR | $XPUT1N\u2211i=1Nxi\u2276\omega T\omega B\alpha $ | Optimal in homogeneous clutter (the baseline detector) | Susceptible to nonhomogeneous clutter in the boundary ring |

SOCA-CFAR | $XPUTmin{.}\u2276\omega T\omega B\alpha $ | Designed to handle strong clutter returns in boundary ring | Susceptible to clutter edges in the boundary ring |

GOCA-FAR | $XPUTmax{.}\u2276\omega T\omega B\alpha $ | Perform well on clutter edges | Susceptible to strong returns appear in the boundary ring |

OS-CFAR | $XPUTx(Q)\u2276\omega T\omega B\alpha $ | Perform well in heterogeneous/nonhomogeneous clutter backgrounds | Susceptible to clutter transitions |

Unlike one-parameter distribution models discussed earlier, more realistic two-parameter distribution models characterize the clutter in the boundary ring of the CFAR stencil by two parameters (mean and variance, scale and shape, etc.). Examples of two-parameter distribution models include log-normal distribution and Weibull distribution. For high-resolution SAR imagery, compound two-parameter distributions such as $K$-distribution, $Go$-distribution, and $\beta $’-distribution are typically used.

A conventional two-parameter CA-CFAR algorithm based on the log detector has the form^{14}^{,}^{52}^{,}^{53}

Accordingly, rearranging Eq. (23) yields

PUT is assumed to be a single pixel. If more than one pixel is assumed, then PUT is the ML estimate of the arithmetic mean such that

However, note that this will entail replacing^{54}$\sigma ^B\u2009log$ in Eqs. (23) and (24) with $\sigma ^m\u2009log$:

^{19}

^{,}

^{21}

^{,}

^{22}

^{,}

^{54}

If we keep the assumption implicit in Eqs. (23) and (24) (i.e., magnitude image is Rayleigh distributed, and power intensity image is exponential distributed), the two-parameter CFAR applied to the image in the (nonlog) domain is given^{55} by

As it is explained thus far, the CFAR approach aims at maintaining a CFAR by locally adapting the detection threshold to the background clutter in the SAR image. However, a detection loss, commonly referred to as CFAR loss, is the price paid for this threshold adaptation. This is due to the fact that, in real-world applications, the noise level in the boundary ring is not constant and/or the number of reference pixels used in the estimation is not large enough. Further, CFAR loss can be viewed as the required increase in the signal to noise ratio (SNR) in order to maintain the desired PD.^{56} The value of this CFAR loss is dependent upon a number of factors, including CFAR method used (e.g., CA, GOCA, SOCA, etc.), number of pixels in the boundary ring, PFA, target statistics, clutter statistics, and noise. Antipov and Baldwinson^{57} offer a notable work on this topic, though applied to 1-D radar data.

A universal curve for CFAR loss in a single-hit detection for steady or Rayleigh target^{56} is depicted in Fig. 6. The parameter $\chi $ is given by

^{56}).

CFAR Method | SAR image type | $k$ |
---|---|---|

CA-CFAR | Square-law detector | $k=0$ |

Linear envelope detector | $k=0.09$ | |

Log detector | $k=0.65$ | |

GOCA-CFAR | Square-law detector | $k=0.37$ |

Linear envelope detector | $k=0.5$ | |

Log detector | $k=1.26$ |

We now present a final word on the log-detector CFAR. Indeed, caution should be exercised when dealing with log-CFAR. As reported in ^{40} log-CFAR introduces an additional CFAR loss into the process. To circumvent this CFAR loss, the size of the CFAR stencil, more precisely the boundary ring in the stencil, needs to be increased by 65 percent^{58} by following

Indeed, despite this additional hurdle, log-detection CFAR is shown to be preferred over linear threshold CFAR processing for typical nonhomogeneous clutter conditions where background clutter surrounding the target is contaminated by other targets, bright clutter pixels, etc.^{59}^{,}^{60}

This section highlights the fact that the CFAR detectors described earlier that are originally designed under the Gaussian assumption (i.e., the background clutter distribution in the $I$ and $Q$ channels) are typically applied in the literature to nonGaussian clutter. Typically, a suitable clutter model (Weibull distribution, $K$-distribution, $\beta $’-distribution, etc.) is adopted, wherein the distribution model parameters are estimated from the boundary ring in the sliding stencil and are used to estimate the scaling factor ($\alpha $) for the desired PFA.

In ^{19}, the two-parameter CFAR scheme is originally designed under Rayleigh/Gaussian assumptions (similar to the one presented earlier) and applied to nonGaussian data in high-resolution SAR images. In Armstrong and Griffiths,^{61} the performance of two-parameter CA-CFAR, GOCA-CFAR, and OS-CFAR (originally designed under Gaussian assumption) is evaluated under $K$-distributed clutter. In Refs. ^{31} and ^{62} two-parameter CA-CFAR (originally designed based on Gaussian assumption) is applied to $\beta $’-distributed high-resolution SAR data. Additional examples are provided in Refs. ^{53} and ^{54}. Obviously, such CFAR detection schemes are applied to clutter distributions that are different from those on which the CFAR detector was originally designed. Subsequently, the conclusion is reached that CFAR techniques prove to be effective regardless of this fact, which explains their popularity.^{63} However, CFAR loss is the price paid for this process.^{40}

In the parametric CFAR methods explained in the preceding sections, the background clutter and/or the target pixels are assumed to obey a certain distribution, and the pixels in the reference window (i.e., boundary ring) are used to estimate the corresponding model parameters. Nonparametric CFAR algorithms do not assume any prior model for the background or the target. Rather, they use nonparametric methods^{18}^{,}^{64} to directly infer the model from the SAR data. An example on a nonparametric CFAR using kernel density estimation (KDE) for inferring the background and the target models is proposed in Gao.^{65}

While most detection techniques reported in the literature rely on a rectangle-shaped hollow stencil with a suitable size and a guard ring, there are some other methods that replace the rectangle-shaped window with another shape. One such example is provided in ^{66}, where the rectangle-shaped CFAR window is replaced with a 2-D gamma kernel. The method utilizes a CFAR detection strategy and is referred to as gamma-CFAR. Similar to traditional rectangle-shaped-window CFAR, the CFAR threshold in gamma-CFAR is estimated from the radial pixel intensity information around a PUT. Further, gamma-CFAR has a free parameter which can be used to estimate the size of its region of support and its standard deviation required for the CFAR test. Although it is not implemented in ^{66} it is stated that this free parameter can be set adaptively. The gamma stencil is also applied in a non-CFAR fashion based on a quadratic Gamma discriminant (QGD) that has eight free parameters.^{66} Test results show that QGD outperforms gamma-CFAR. Extension of QGD to an artificial neural network (ANN) based on a multilayer perceptron (MLP) is provided in ^{67}.

Besides CFAR-based methods, there are other approaches that do not use CFAR. For example, in Ouchi et al.,^{68} the detection is based on a coherence image produced from the multilook SAR image via cross correlation between two SAR images extracted by moving windows of a small size over the original image. The method can detect objects buried in a speckle. In Howard, Roberts, and Brankin,^{69} the detection is based on genetic programming. Relevant work is provided in Conte, Lops, and Ricci,^{70} although applied to 1-D radar data, where the detection is based on a generalized likelihood ratio test (GLRT) in a $K$-distributed clutter. The null hypothesis ($Ho$) represents the $K$-distributed clutter, and the alternative hypothesis ($H1$) is modeled as being a compound of target signal and a $K$-distributed noise.

All the aforementioned methods for target detection are single-feature-based in that they distinguish the target pixels from the background only on the basis of the level of pixel brightness (i.e., RCS). This poses a limitation on these methods, more significantly, in ROIs with heterogeneous clutter and/or contiguous targets.

Methods under the multifeature-based taxon try to circumvent this drawback by basing the detection decision on a fusion of two or more features. Obviously, this taxon can utilize a suitable method among those presented under the single-feature-based taxon and incorporate additional features besides RCS. Methods that fall under the multifeature-based taxon can be broadly classified into two major sub-taxa: those that utilize arbitrary user-chosen features and those that rely on systematic multiresolution analysis. Examples of arbitrary user-chosen features are provided in Refs. ^{9} and ^{71}, wherein decisions on ROIs in the SAR image are based on a fusion of three multistage features extracted in parallel from the SAR image, namely CFAR features, variance features, and extended fractal (EF) features. Obviously, this approach is not purely CFAR. Another relevant example of the multifeature method is provided in ^{72}. Further, in Subotic et al.,^{73} parametric circular complex Gaussian models are utilized for both target and clutter. However, unlike the traditional CFAR approach, which works on a single resolution, the resolution of the SAR image in this approach is varied to produce multiple features. It is based on the conjecture that target signatures, when viewed as a function of resolution, are sufficiently different from clutter. This allows detection performance gains over single-resolution detection algorithms.

This perspective motivates the more systematic multiresolution analysis methods, which can be broadly classified into space-scale-based and space-frequency-based. Space-scale methods produce space-scale features based on the wavelet transform, including the DWT, and CWT. In ^{74} the detection is based on DWT that yields a spatial correlation of the sub-bands that attenuates the background noise and enhances the structured pattern of the target signature. A description of the more relevant detection strategies based on *CWT* can be found in Antoine et al.^{75}

Prime examples of methods that utilize space-frequency features include linear space-frequency methods, such as the S-transform,^{76} and bilinear (also known as quadratic) space-frequency methods, such as Cohen’s class distributions^{77}^{,}^{78} (Wigner distribution, Wigner-Ville distribution, pseudo-Wigner-Ville distribution, etc.).

Expert-system-oriented target detection is a multistage (two or more stages) AI approach that bases the detection process on exploitation of prior knowledge about the imaged scene and/or target(s). Prior knowledge is captured via context utilization. In the broader sense, context here refers to all the available means that can help to convey information about the background clutter and/or target(s) of interest in the scene. Such means include image segmentation, scene maps, digital elevation model (DEM), previously gathered data, and geographical information system (GIS).

In its simplest form, context utilization in CFAR can be realized via unsupervised/semisupervised SAR image segmentation. Typically, prior to applying CFAR to the SAR image, the SAR image is segmented to extract an image structure map. Then, the conventional CFAR method of choice is aided with this map to enable it to adaptively select the suitable area over which the background statistics can be properly estimated. Further, smaller segments can be labeled as potential targets, while larger ones can be labeled as background. An example for utilizing annealed segmentation with one-parameter CFAR for SAR target detection is provided in McConnell and Oliver.^{79}

A relevant work is termed variability index CFAR (VI-CFAR), although applied on 1-D range data.^{80}^{,}^{81} Based on the estimations of the mean in the boundary ring of the stencil, VI-CFAR switches between one of the following three CFAR methods: CA-CFAR, SOCA-CFAR, and GOCA-CFAR. This approach tries to benefit from the strength of each CFAR method through deployment of the switching scheme.

Another relevant algorithm is reported in Gao et al.^{62} The algorithm builds on the beta-prime CFAR ($\beta $’-CFAR) algorithm^{31} reported earlier in this paper. A binary index map is created based on globally thresholding the input SAR image. The index map is comprised of zeros assigned to pixels in the input SAR image that are found to be less than a predetermined global threshold, and ones assigned to pixels found to be greater than the global threshold. Then, a sliding window stencil is placed over the image (i.e., pixel by pixel) where the parameters of the $Go$-distribution (i.e., $Go$-distribution reduces to $\beta $’-distribution for a single look SAR image; more on this is given in Sec. 5.2) are estimated from the pixel values in the boundary ring correspond to $0$ in the index map. The size of the stencil is chosen based on the prior knowledge of the target size as described in Salazar.^{31} Then, a local threshold in the stencil is determined based on the parameters estimated, and a decision on the PUT is taken. The window is then systematically slid to subsequent pixels in the image until the whole image is scanned. Further refinements on detections are achieved by placing a constraint on the size and allowable distance between detections.

More robust expert-system-oriented approaches utilize a mixture of multiple different CFAR experts (CA-CFAR, OS-CFAR, etc.), each of which is specialized to handle a suitable kind of clutter. The expert-system-oriented system uses available context information, extracted by one or more means (as explained earlier in this section), to assign the CFAR experts to suitable regions in the SAR image. Unfortunately, there is a lack of work published on this approach for SAR imagery. However, there are works published on 1-D radar data. One such example^{82}^{–}^{84} was implemented by the U.S. Air Force Research Laboratory.

Rimbert and Bell^{43}^{,}^{44} present another interesting work. It is motivated by the observation that, in homogeneous regions (i.e., locally stationary clutter), a larger reference window provides a clutter estimate that has a mean with smaller variance than a similar estimate based on a smaller reference window. Conversely, a smaller reference window provides a more reliable mean estimate in nonhomogeneous regions. Accordingly, a target detection scheme that adapts the reference window size and selects one of two CFAR detectors (i.e., CA-CFAR, and OS-CFAR) based on the type of the reference region is proposed. The proposed detection scheme checks the homogeneity of the reference region using a goodness-of-fit statistical test with an assumed parametric model for the clutter. One can think of this process as producing a structure map for the reference window (i.e., not for the whole image, as discussed in this section for the segmentation-based method). A CFAR method termed cell under test inclusive CFAR (CI-CFAR) is also introduced, which operates in a manner similar to OS-CFAR but combines the PUT with the clutter pixels in the reference window. Although the authors reported in one work^{43} that CI-CFAR outperforms OS-CFAR, in another work,^{44} the authors acknowledge that this conclusion is incorrect, due to the use of a simulation with errant detection threshold parameters. It should be noted that the results presented are not based on 2-D SAR data. Further, the analysis was based on the assumption that both the clutter and the target obey the central limit theorem (CLT) (i.e., Gaussian), which is a nonrealistic assumption for high-resolution SAR imagery. Additionally, reference stencils used are not hollow, and they do not consider any guard regions.

Finally, any relevant detection method that utilizes any form of intelligence/inference is of particular interest and to fit under this taxon.

A comparison between selected examples, most of which are from amongst those cited in Sec. 3, on the detection module pertaining to SAR imagery is provided in Table 3. The choice of the examples was carefully made to cover the different methods under each taxon. Comparison aspects include SAR image type, feature(s), clutter/target type, clutter model (if applicable), and target model (if applicable). The comparison does not attempt to assess the algorithm performance based on reported PD and PFA, as this is infeasible, given the variant types of methods and data used, as well as the different sensor characteristics and operating conditions. Rather, the aim of this comparison is to show the major differences between the various methods and depict their applicability to certain scenarios.

$T$ | ST | Refs. | Image type | Feature(s) | Clutter/target | Clutter model | Target model | Comments |
---|---|---|---|---|---|---|---|---|

Single-feature-based | CFAR-Based | Parametric methods: based only on background modeling | ||||||

^{85} | SIR-C/X SAR data | One/RCS | Sea/ship | Joint log-normal | N/A | Uses sliding window CA-CFAR. Tested on isolated targets in a homogeneous clutter. | ||

^{33} | HH/airborne/high resolution | One/RCS | Homogeneous/glinting area targets | Exponential | N/A | Uses sliding window OS-CFAR. KDE is used to estimate PFA and PD. Logarithm of intensity is considered for calculations. | ||

^{26}, ^{27} | Spotlight/HH/airborne/high resolution | One/RCS | Land/extended objects | Weibull | N/A | Uses sliding window CFAR. Homogeneous clutter in local windows is assumed. Location-scale type (BLUE estimate)^{86} is used for parameter estimation. | ||

^{28}, ^{86} | Lincoln Lab/HH and HV | One/RCS | Land/bridge, and power-line tower | K-distribution | N/A | Uses sliding window CFAR. OS-CFAR and CA-CFAR are compared under Weibull and $K$-distribution clutters. | ||

Parametric methods: based on background and target modeling | ||||||||

^{35}, ^{36} | High resolution (no info provided) | One/RCS | Land/vehicles | Log-normal | Log-normal | Uses sliding window-CFAR. (not AD). NP criterion is utilized. | ||

Nonparametric methods: based on kernel density estimation (KDE) | ||||||||

^{65} | Radarsat-1/space-borne/high resolution | One/RCS | Sea/ship (homogeneous clutter, and isolated targets) | Nonparametric. Estimated using KDE | N/A | Uses sliding window CFAR, for AD. Gaussian kernel is used for the KDE of the background model. | ||

Non-CFAR-based | Based on a coherence Image | |||||||

^{68} | Radarsat-1 | One/RCS (correlation-based) | Sea/ship (homogeneous clutter, and isolated targets) | N/A | N/A | A coherence image produced from the multilook image via cross correlation between two images extracted by moving windows of a small size over the original image. Can detect objects buried in speckle noise. | ||

Based on genetic programming | ||||||||

^{69} | Low resolution SAR imagery (ERS data) | One/RCS | Sea/ship (homogeneous clutter, and isolated targets) | N/A | N/A | Two-stage evolution strategy. | ||

Nonrectangle/nonsquare-shaped-stencil methods | ||||||||

CFAR | ^{66} | Fully polarimetric MIT high resolution SAR data. PWF is used to generate a single image. | One/RCS | Natural and cultural clutter (including a parking lot)/vehicles | Gamma kernel with three parameter | N/A | The method is termed gamma-CFAR, and uses two 2-D gamma kernels to form the sliding stencil. PWF stands for polarimetric whitening filtering. | |

Non-CFAR | ^{66} | Fully polarimetric MIT high resolution SAR data. PWF is used to generate a single image. | One/RCS | Natural and cultural clutter/vehicles | QGD with eight parameters | N/A | The method is based on 2-D gamma kernels but replaces CFAR with QGD. QGD is shown to outperform gamma-CFAR. | |

Multifeature-based | CFAR Utilized | Based on multifeature fusion | ||||||

^{9}, ^{71} | MSTAR, and TESAR emulated imagery. | Three features (RCS, variance, & extended fractals) | Heterogeneous clutter/vehicles (isolated targets). | Gaussian | N/A | Features generation process is multistaged (four stages). Among other features, uses sliding window CFAR. Then detection of ROIs is based on fusion of the features. | ||

CFAR-Based | Based on multiresolution | |||||||

^{73} | Synthetic X-band SAR, and real SAR imagery | RCS at multiresolutions | Land/vehicles (quite homogeneous clutter, and isolated target) | Zero-mean circular complex Gaussian | Circular Complex Gaussian | Sliding window-CFAR (Not AD). NP criterion is utilized. It is based on that target signatures, when viewed as a function of resolution, are sufficiently different from clutter. This allows detection performance gains over single-resolution detection algorithms. | ||

Non-CFAR-based | ^{74} | Synthetic (emulated), and Radarsat-1 | One/RCS | Sea/ship (homogeneous clutter, and isolated targets) | N/A | N/A | (Space-scale analysis) Application of DWT that yields a spatial correlation of the sub-bands attenuates background noise and enhances the structured pattern of ship signature. | |

Expert-System-Oriented | CFAR-Based | Based on segmentation | ||||||

^{79} | Simulated targets manually inserted in real and synthetic SAR images | Image structure, and RCS | Non-homogeneous clutter/simulated targets on land inserted in the clutter. | Constant background is assumed. | N/A | Two-stage process including, segmentation and CFAR. CFAR criterion is based on AD. No model is assumed for the target. | ||

Based on index matrix | ||||||||

^{62} | Real SAR imagery, X-band, HH | Global thresholding, RCS, and target size | Nonhomogeneous clutter/vehicle targets | $\beta \u2032$-Distribution | N/A | The CFAR detection process is guided via a binary index matrix. Detections are refined via checking the distance between detection segments. The algorithm design utilizes the beta-prime CFAR algorithm reported in ^{31}. |

Obviously, target detection based on parametric modeling of SAR imagery is the most popular in the literature. This discussion focuses on parametric CFAR methods that utilize stochastic models for modeling the background clutter. Primarily, the issue of the suitability of these models to represent the SAR data is briefly discussed. This is followed by a concise summary of the popular multiplicative (also known as compound) SAR data models, conditions of their applicability, and the interrelation between them. Finally, our focus is shifted to parametric CFAR. We tame CFAR and present some novel discussions from two different perspectives: the signal processing perspective and the statistical pattern recognition perspective.

In Table 3, there are various parametric clutter models used, including log-normal-distribution, $K$-distribution, and exponential distribution. Accordingly, in Table 4, we briefly summarize some of the major statistical distributions, along with the backscatter types in SAR images they are typically used to model.

Distribution | Backscatter type/comments |
---|---|

Normal, and Rayleigh^{87} | Homogenous (i.e., bare ground surfaces, dense forest canopies, snow covered ground). |

Weibull, and log-normal^{87} | Other clutter types such as sea surface. |

Modified beta^{88} | Different ice types. |

$K$-distribution^{89} | Models heterogeneous backgrounds. |

It offers a multiplicative model (compound distribution). | |

$K$-distribution originally proposed in ^{90} for modeling microwave sea echoes. | |

It then became popular for modeling multilook^{91} and polarimetric SAR signature.^{92} | |

It has much poorer performance in extremely heterogonous clutter such as urban clutter.^{93} | |

$Go$-distribution^{93} | Models extremely heterogeneous clutter background. Thus, it has better performance than $K$-distribution.^{93}^{,}^{94} |

It is a compound distribution. | |

$G$-distribution^{93} | It can model extremely heterogeneous clutter such as urban regions that $K$-distribution cannot.^{93}^{,}^{94} |

It is a compound distribution. | |

$K$ and $Go$ distributions are special cases of this class. |

As is evident in Table 3 (and relevant works published in the literature), in many cases, similar distributions are being randomly chosen to model the clutter at various sensor characteristics, such as frequency, polarization, imaging mode (e.g., Spotlight, Stripmap, ScanSAR, etc.), and resolution. Obviously, the choice of a proper model for the clutter backscatter depends, not only on the clutter type, but also on these sensor characteristics. Several works^{89}^{,}^{95}^{,}^{96} have noted that the suitability of a certain distribution to model a certain clutter depends on the data being used, as well as the corresponding sensor characteristics and operating conditions. Surprisingly, this issue is generally overlooked in many works published in the literature. In many such works, the justifications for choosing some model (e.g., $K$-distribution) is merely based on the assumption that the distribution model is found suitable for a certain clutter type (e.g., ocean clutter) in some published work, and thus it can be automatically adopted for a similar clutter. Indeed, such a conclusion can be misleading.

Accordingly, prior to randomly opting for a popular parametric distribution to model the clutter, one should consider validating the applicability of the distribution on the data using a suitable goodness-of-fit technique.^{97}^{,}^{98} One such validation approach for SAR data is based on the Cramer-Von Mises (CVM) distance^{97} as presented in di Bisceglie and Galdi.^{26} The normally used Kolmogorov-Smirnov test is discarded, because the independency assumptions of its usage are violated, due to the dependency of the SAR data being non-Gaussian distributed and generally correlated.^{97} Thus, the CVM method is used to measure the distributional distance between the design cumulative distribution function (CDF) model and the empirical CDF estimated from the available high-resolution SAR image. The design distribution that scores a minimum distance within some threshold is typically chosen. This procedure is demonstrated in di Bisceglie and Galdi^{26} on Spotlight (9.6 GHz, HH polarized, geometric resolution of $0.4\u2009\u2009m\xd70.4\u2009\u2009m$) SAR data (with a specific clutter type) pertaining to the Rayleigh and Weibull distributions. On the SAR image used in di Bisceglie and Galdi,^{26} the Weibull distribution is found to be a more suitable model for the clutter.

One final observation on this matter is that, though the CVM test was used in the work mentioned above, CVM is not the best technique when it comes to characterizing tailed SAR data. Indeed, the $k$-sample Anderson-Darling test^{99} offers a better procedure, as it places more weight on observations in the tail of the CDF distribution. Both the CVM test and Anderson-Darling test belong to the quadratic class of empirical distribution functions (EDF) for statistical tests.^{99} Finally, information-theoretic approaches for characterizing the goodness-of-fit can be found in ^{100}

In any target detection scheme that depends on parametric modeling (e.g., the popular CFAR-based), the selection of an appropriate probability distribution to model the pixels in the SAR image (i.e., radar backscatter) is a must, because the thresholding operation in any such detector is dependent on the clutter distribution. In cases where the random scatterers in a resolution cell in the SAR image have sizes on the order of the wavelength of the radar signal, the total backscatter can be modeled as the sum of isolated returns in the cell.^{101} This invokes the CLT, wherein the $I$ and $Q$ components of the total complex-valued backscatter can naturally be assumed to be normally distributed. This implies that the total backscattered amplitude and phase can be modeled as Rayleigh and uniform distributions, respectively. Thus, the power in each resolution cell is modeled as an exponential distribution. Conversely, in high-resolution SAR, the above mentioned assumptions are violated, because the number of random scatters in a resolution cell is not large, and thus the CLT cannot apply. This renders the clutter non-normally distributed, which motivates the need for a suitable model.

The multiplicative model (also known as the compound model) for SAR image formation has been popularly used in the literature to model the clutter background. The model is based on the hypothesis that the SAR image is formed from the product of a backscatter and speckle random processes as

Typically, for a single look (i.e., $n=1$) SAR image, the complex speckle $Y$ is characterized as bivariate normal density for complex imagery, which reduces to the exponential distribution in the power or intensity domain. Further, for multilook imagery (i.e., $n>1$), the two-parameter gamma distribution $\Gamma (\alpha ,\lambda )$ characterizes the speckle in the power domain, and this reduces to the square root gamma distribution $\Gamma (\alpha ,\lambda )$ in the magnitude domain. Depending on the type of the background clutter (i.e., homogeneous, heterogeneous, or extremely heterogeneous), and the pertinent sensor characteristics (i.e., operating conditions) such as frequency, polarization, and gazing angles, several different distributions are used in the literature to model the backscatter $X$. For each case, the manner in which $X$ (and subsequently $Z$) manifests itself depends^{31}^{,}^{94} on whether the SAR image is single-look (i.e., $n=1$) or multilook (i.e., $n>1$).

First, for homogeneous regions and a single look SAR image, $X$ is typically modeled as a constant that equals the average power in the homogeneous region (i.e., $C=2\sigma G2$). Accordingly, the power-domain SAR image $Z$ is modeled as exponential distributed: