Detection Theory for Radar

From Optimal Detection to Adaptive Thresholding

Having formed the range-Doppler map, we must decide which cells contain targets and which are noise. The Neyman-Pearson framework (FSI Ch 01) gives the optimal detector when the noise statistics are known. In practice, the noise level varies with range, clutter, and interference — so the threshold must adapt. This is the role of CFAR (constant false alarm rate) detectors.

For imaging, detection applies per-pixel: each voxel in the reconstructed image is tested against a local noise estimate.

Definition:

Neyman-Pearson Detection in Radar

The binary hypothesis test for a single range-Doppler cell:

H0:r=w(noise only),H_0: r = w \qquad (\text{noise only}), H1:r=αs+w(target present),H_1: r = \alpha s + w \qquad (\text{target present}),

where rr is the matched filter output, α\alpha is the unknown target amplitude, and wCN(0,σ2)w \sim \mathcal{CN}(0, \sigma^2).

The Neyman-Pearson detector compares the test statistic r2|r|^2 to a threshold γ\gamma:

r2H1H0γ,|r|^2 \underset{H_0}{\overset{H_1}{\gtrless}} \gamma,

where γ\gamma is set to achieve a desired false alarm probability: Pfa=eγ/σ2P_{\text{fa}} = e^{-\gamma/\sigma^2} (for complex Gaussian noise with known variance σ2\sigma^2).

Definition:

Cell-Averaging CFAR (CA-CFAR)

When the noise variance σ2\sigma^2 is unknown and varies with range, the CA-CFAR detector estimates it from neighboring cells.

For a cell under test (CUT) at range bin mm, the noise estimate uses 2Nr2N_r reference cells (excluding NgN_g guard cells on each side):

σ2^=12Nrirefri2.\hat{\sigma^2} = \frac{1}{2N_r}\sum_{i \in \text{ref}} |r_i|^2.

The detection rule is:

rm2H1H0αCFARσ2^,|r_m|^2 \underset{H_0}{\overset{H_1}{\gtrless}} \alpha_{\text{CFAR}} \cdot \hat{\sigma^2},

where the threshold multiplier for a desired PfaP_{\text{fa}} is:

αCFAR=2Nr(Pfa1/(2Nr)1).\alpha_{\text{CFAR}} = 2N_r \left(P_{\text{fa}}^{-1/(2N_r)} - 1\right).

The guard cells prevent target energy from leaking into the noise estimate. The reference cells should contain only noise — if a strong target is in the reference window, the noise estimate is inflated and nearby weak targets are masked.

Definition:

Order-Statistic CFAR (OS-CFAR)

The OS-CFAR detector replaces the mean with an order statistic to be robust against interfering targets in the reference window.

Sort the 2Nr2N_r reference cell powers: r(1)2r(2)2r(2Nr)2|r_{(1)}|^2 \leq |r_{(2)}|^2 \leq \cdots \leq |r_{(2N_r)}|^2.

Choose the kk-th order statistic as the noise estimate:

σ2^=r(k)2,\hat{\sigma^2} = |r_{(k)}|^2,

typically with k3Nr/4k \approx 3N_r/4.

The detection rule is rm2>αOSr(k)2|r_m|^2 > \alpha_{\text{OS}} \cdot |r_{(k)}|^2.

OS-CFAR is more robust than CA-CFAR at clutter edges and in multi-target environments, at the cost of slightly higher CFAR loss (reduced PdP_d for the same PfaP_{\text{fa}}).

Theorem: Detection Probability for Swerling I Targets

For a Swerling I target (exponentially distributed RCS, constant across pulses within a CPI) with average SNR\text{SNR} per pulse γˉ\bar{\gamma}, the single-pulse detection probability is:

Pd=Pfa1/(1+γˉ).P_d = P_{\text{fa}}^{1/(1 + \bar{\gamma})}.

For NpN_p pulses with non-coherent integration (square-law combining), the detection probability increases but lacks a closed-form expression — it is computed via:

Pd=1Fχ2Np2(γ1+γˉ),P_d = 1 - F_{\chi^2_{2N_p}}\left(\frac{\gamma}{1 + \bar{\gamma}}\right),

where Fχ2Np2F_{\chi^2_{2N_p}} is the CDF of a chi-squared distribution with 2Np2N_p degrees of freedom and γ=σ2lnPfa\gamma = -\sigma^2 \ln P_{\text{fa}}.

Example: CA-CFAR Threshold Computation

A radar uses CA-CFAR with Nr=16N_r = 16 reference cells on each side and Ng=3N_g = 3 guard cells. Compute the threshold multiplier αCFAR\alpha_{\text{CFAR}} for Pfa=106P_{\text{fa}} = 10^{-6}. What is the CFAR loss compared to the ideal fixed-threshold detector?

CFAR Detection Visualization

Apply CA-CFAR and OS-CFAR to a range profile with targets and clutter. Observe how the adaptive threshold tracks the noise floor and how targets near clutter edges are affected.

Parameters
-4
32

Definition:

M-out-of-N Binary Integration

Binary integration (also called post-detection integration) is a non-coherent technique: apply CFAR detection to each pulse independently, then declare a target present if at least MM out of NN pulses produce a detection.

For independent pulses with single-pulse detection probability Pd,1P_{d,1} and false alarm probability Pf,1P_{f,1}:

Pd=k=MN(Nk)Pd,1k(1Pd,1)Nk,P_d = \sum_{k=M}^{N} \binom{N}{k} P_{d,1}^k (1-P_{d,1})^{N-k}, Pfa=k=MN(Nk)Pf,1k(1Pf,1)Nk.P_{\text{fa}} = \sum_{k=M}^{N} \binom{N}{k} P_{f,1}^k (1-P_{f,1})^{N-k}.

Binary integration is suboptimal compared to coherent integration but does not require pulse-to-pulse phase coherence.

Common Mistake: Target Masking in CA-CFAR

Mistake:

Using CA-CFAR when multiple targets may exist in the reference window, assuming that the noise estimate is unbiased.

Correction:

A strong target in the reference window inflates σ2^\hat{\sigma^2}, raising the threshold and masking nearby weaker targets. This is the target masking problem.

Mitigations: (1) Use OS-CFAR, which selects an order statistic and is robust to a few outliers. (2) Use censored CA-CFAR, which removes outliers before averaging. (3) Use greatest-of (GO-CFAR) which takes the maximum of left and right reference windows, better handling clutter edges.

Historical Note: The Development of CFAR

1960s-1980s

The CFAR concept was developed in the 1960s as radar systems moved from operator-interpreted displays (where the human eye naturally adapted to varying noise levels) to automatic detection. Finn and Johnson (1968) published the first analysis of cell-averaging CFAR for exponential noise. Rohling (1983) introduced OS-CFAR to handle multi-target environments. The progression mirrors a broader theme in radar: replacing human judgment with principled statistical methods, exactly the path from visual interpretation to algorithmic reconstruction in RF imaging.

CFAR (Constant False Alarm Rate)

A detection scheme that adaptively estimates the noise power from reference cells surrounding the cell under test, maintaining a constant false alarm probability despite varying noise levels.

ROC (Receiver Operating Characteristic)

A curve plotting detection probability PdP_d against false alarm probability PfaP_{\text{fa}} as the decision threshold varies. Parameterized by SNR\text{SNR} and target model (Swerling type).

Quick Check

For a Swerling I target with Pfa=106P_{\text{fa}} = 10^{-6} and average SNR=15\text{SNR} = 15 dB (γˉ=31.6\bar{\gamma} = 31.6), what is the single-pulse detection probability?

Pd0.65P_d \approx 0.65

Pd0.90P_d \approx 0.90

Pd0.99P_d \approx 0.99

Pd0.30P_d \approx 0.30

Key Takeaway

CFAR detectors maintain constant PfaP_{\text{fa}} by adaptively estimating the noise power from reference cells. CA-CFAR is optimal for homogeneous noise but suffers from target masking; OS-CFAR and GO-CFAR are more robust. For Swerling I targets, Pd=Pfa1/(1+γˉ)P_d = P_{\text{fa}}^{1/(1+\bar{\gamma})} — fluctuating targets require significantly more SNR\text{SNR} than non-fluctuating ones.

CFAR Detector: Adaptive Threshold in Action

A cell-averaging CFAR detector slides along a range profile containing targets embedded in spatially varying clutter. The guard cells, reference cells, and the adaptive threshold are shown evolving in real time. Detections are flagged where the cell under test exceeds the local threshold, illustrating how CFAR maintains constant false-alarm rate despite non-stationary background power.