## 1.

## Introduction

Recently, the detection of moving small dim targets under heavy IR background clutter has been an active research area, and many algorithms for this problem have been developed.^{1, 2, 3, 4} The algorithms used in IR searching and tracking systems are adequate for applications with bright targets against background clutter, and make use of only the spatial information of the targets and clutter without considering their temporal behavior. Some new approaches^{2, 3} incorporated temporal and spatial information have good performance for the detection of small dim targets in IR image sequences, but heavy computational complexity is a frustration in practical applications.

Since that the temporal profile of a pixel through which a target passes is different from those through which clutter pass, Tzannes and Brooks^{3} developed the models of temporal behavior of clutter and targets on a single-pixel basis. These are used to develop a temporal likelihood-ratio test for target detection. Being different from the mentioned algorithm based on 1-D signal procession, our algorithm views the fluctuation of the gray-scale value of a pixel, while a target passes, as a weak transient signal. Owing to the good performance for transient signal detection in Gaussian noise, the power-law detector^{5, 6} is employed in our algorithm to detect small dim targets.

The rest of the paper is organized as follows. Section 2 describes the power-law detector and the motivation behind the use of the power-law detector for small dim target detection in an IR image. Section 3 describes our algorithm. Section 4 presents experimental results with corresponding analyses. Finally, conclusions are given in Sec. 5.

## 2.

## Power-Law Detector

The power-law detector, developed by Nuttall,^{5} performs well for transient signal detection without knowledge of the signal, hence its wide applications. In a Gaussian background, Nuttall considered that the transient signal detection can be transformed into the problem of arbitrary
$M$
bits signal detection among the
$N$
discrete Fourier transform (DFT) bins, where
$M$
is the number of signal-present bins. For a 1-D signal
$x\left(t\right)$
including a transient signal, the power-law statistic can be formulated as

## 1

$$\sum _{j=1}^{N}{\mid X\left(j\right)\mid}^{v}\{\begin{array}{ll}\u2a7d\lambda & \text{signal}\phantom{\rule{0.3em}{0ex}}\text{absent}\\ >\lambda & \text{signal}\phantom{\rule{0.3em}{0ex}}\text{present},\end{array}\phantom{\}}$$Given the contiguity tendency of a transient signal in frequency, Willett and Wang^{6} modified
$X\left(j\right)$
in Eq. 1 and obtained new random variables, i.e.,
${W}_{j}={X}_{j-1}+{X}_{j}$
,
$j=1,\dots ,N$
by combining two contiguous frequency bins. A new power-law detector is defined as

## 2

$$T\left(W\right)=\sum _{j=1}^{N}{W}_{j}^{v}=\sum _{j=1}^{N}{({X}_{j-1}+{X}_{j})}^{v}\{\begin{array}{ll}\u2a7d\lambda & \text{signal}\phantom{\rule{0.3em}{0ex}}\text{absent}\\ >\lambda & \text{signal}\phantom{\rule{0.3em}{0ex}}\text{present}.\end{array}\phantom{\}}$$The power-law detector and its improved version are nonparametric detectors, and prior knowledge of the transient signal such as the model or the frequency spectrum is not necessary. According to the assumption that the temporal behavior of the background noise is additive Gaussian noise,^{3} the fluctuation of the gray-scale value of a pixel in temporal profiles caused by a passing target can be viewed as a transient signal, and its detection is transformed into the problem of transient signal detection in Gaussian noise.

## 3.

## Algorithm Implementation

The algorithm consists of two parts, clutter suppression in the spatial domain and power-law detection in the time domain, as shown in Fig. 1.

First, a time window with length of
$L$
frames is introduced. The temporal profile of each pixel can be viewed as a 1-D signal in the time window. To decrease the computational burden, the current image can be processed by a max-median filter^{7} in the spatial domain to suppress the clutter and obtain the candidate target pixels. Then the power-law detector is used to detect real target on the temporal profiles of the candidate target pixels. Nevertheless, to let the “pulse” fully present in the time window, the current image is put at the center of the time window. Then, the length of the window is related to the size and the speed of the transient signal.^{3} Generally,
$L$
is set as 7 to 15.

With the candidate pixels by the max-median filter^{7} just obtained, the small dim targets detection problem can be formulated as a binary hypothesis test:

## 3.

## 4

$$T\left(W\right)=\sum _{j=1}^{N}{W}_{j}^{v}=\sum _{j=1}^{N}{({X}_{j-1}+{X}_{j})}^{v}\{\begin{array}{ll}\u2a7d\lambda & {H}_{0}\\ >\lambda & {H}_{1},\end{array}\phantom{\}}$$Similar to the description in the work of Wang and Willet,^{6}
$T\left(W\right)$
is also rewritten as

## 5

$$T\left(W\right)=\sum _{j=1}^{N}{W}_{j}^{v}=\sum _{j=1}^{N}{({X}_{j-1}+{X}_{j})}^{v}=\sum _{j=1}^{N\u22152}{W}_{2j-1}+\sum _{j=1}^{N\u22152}{W}_{2j}={T}_{o}+{T}_{e}.$$^{6}${X}_{j}$ follows an independent and identically distributed (i.i.d) exponential distribution under hypothesis ${H}_{0}$ , and $F(x,y,n)$ is approximately Gaussian noise. Because that ${W}_{2j-1}$ and ${W}_{2j}$ , $j=1,\dots ,N\u22152$ , are i.i.d exponential variables, ${T}_{o}$ and ${T}_{e}$ converge to the distribution $N(\mu ,{\sigma}^{2})$ via the central limit theorem, with

## 6.

## 4.

## Experiments and Analysis

The signal-to-clutter ratio (SCR) is defined as

where ${a}^{2}$ denotes the target intensity, and ${\sigma}^{2}$ is the variance of background clutter.## 4.1.

### Experiment 1

A target of size
$2\times 2$
is embedded in an IR image sequence of 95 frames, which goes through the cloud along a straight line. The parameters of the algorithm are selected as follows: the false alarm probability
${P}_{\mathrm{FA}}={10}^{-4}$
, thus the resulting threshold
$\lambda =3200$
, the length of sliding window
$L=9$
, and the length of the DFT
$N=128$
. Note that the effective detection begins from the fifth frame. Figure 2 shows the detection results with
$\mathrm{SCR}=1$
for our algorithm. Figure 3 illustrates the detection performances with different SCRs of our algorithm and the third-order-cumulant^{4} based algorithm.

When
$\mathrm{SCR}=1.0$
, only a few of possible target pixels are detected after max-median filtering and self-adapted valve segmentation,^{7} as shown in Fig. 2c. Then, for these pixels, our algorithm is used for further detection and the results are shown in Fig. 2d. Note that our algorithm can detect the small gray change caused by a target passing through the clouds. From Fig. 3, it is clear that our algorithm has a better performance than the third-order-cumulant-based algorithm for
$\mathrm{SCR}\ge 0.7$
. However, due to the limit of clutter suppression, the detection probabilities decline sharply when
$\mathrm{SCR}<0.7$
, while the third-order-cumulant-based algorithm also has good performance.

## 4.2.

### Experiment 2

The second sequence is acquired by a PtSi IR cameras with focal plane arrays^{3} of
$320\times 244\phantom{\rule{0.3em}{0ex}}\text{pixels}$
. The selected sequence of 95 consecutive frames is used for the algorithm evaluation. Figure 4a shows a single image of this sequence, which is a daytime scene that includes two point targets and many drifting and evolving clouds. The left target is detected accurately by our algorithm in 85 frames, while the Tzannes and Brooks algorithm detects the target in just 68 frames, as shown in Figs. 4b and 4c. Due to the higher SCR of the right target, there is little difference of detection performance between the two algorithms.

## 5.

## Conclusions

We presented a novel algorithm based on the improved power-law detector for the detection of moving small dim targets in image sequences with heavy clutter. It was heuristically demonstrated that the proposed algorithm is very successful for the detection of pixel-sized target. The effectiveness of the algorithm was demonstrated with two experiments. The experimental results illustrate that our algorithm can promote the performance of small dim target detection under heavy clutter background.

## Acknowledgments

This work is supported by the National Natural Science Foundation of China (No. 60677040).

## References

*Signal and Data Processing of Small Targets*, Proc. SPIE0277-786X 10.1117/12.364049 3809, 74–83 (1999). Google Scholar