1. Introduction
With the development of aerospace technology, the detection of space targets has attracted much attention in recent years. Small space targets, including missile shrapnel and precision-guided weapon seekers, pose a considerable threat to defense security. If these targets are detected, it will significantly reduce threats to our safety. Thus, detecting small space targets greatly challenges our early warning radar. Space targets are usually characterized by a low reflection cross-section (RCS), which results in a low SNR. Therefore, coherent or incoherent integration in a long CPI is an effective way to obtain a high SNR gain. However, space targets usually have high-speed and maneuverable motion, leading to range and Doppler migration and causing severe performance loss of integration gain [
1,
2,
3]. Estimating motion parameters and compensating for the motion phase is an effective way to obtain fully coherent integration gain and detect space targets successfully.
In recent years, many methods have been proposed to realize coherent or incoherent integration of high-speed and maneuverable targets. Generally, they can be divided into three categories: (1) Incoherent integration method. A Hough transform (HT) [
4,
5,
6] is a typical incoherent integration method, which finds the slot of the envelope exceeding the detection threshold and integrates the signal along tracks to realize long time incoherent integration. However, the non-coherent SNR gain loss leads to performance degradation in very low SNR conditions. (2) Decoupling of range frequency and slow time. After down-conversion and pulse compression, the targets’ motion related to slow time couples with range frequency in the range frequency-slow time domain. As for uniform motions, only the linear coupling phase needs to be considered. The typical keystone transformation method [
7,
8,
9] modifies the slow time to decouple the range frequency and slow time for high-speed targets. After decoupling, the range frequency and slow time are independent, and thus, the signal can be coherently integrated into the fast time and Doppler domains independently. For maneuverable target detection, many modified keystone transformation methods have been developed, such as the high-order keystone transformation [
10,
11,
12] and deramp-keystone transformation [
13,
14], which aim to be a specific motion model to realize decoupling and coherent integration. Another decoupling method is based on Radon transformation [
15,
16,
17], which jointly searches range and velocity parameters to realize decoupling. These methods can estimate high-order motion model parameters such as acceleration and jerk [
18,
19]. (3) Time-frequency analysis method. Unlike methods (1) and (2) that integrate in the fast time and Doppler domain, this method uses specific time-frequency distribution to realize parameter estimation and energy accumulation. Lv’s distribution (LVD) [
20,
21,
22,
23] uses a scaled correlation between the adjacent slow times. This way, the quadratic phase can be converted to a linear phase, and the echo can be integrated into the frequency-chirp rate domain. The method breaks through the tradeoff between resolution and cross-term [
24], and parameter searching is unnecessary. The method has good performance in the low SNR scene. The Wigner–Ville distribution (WVD) [
25,
26,
27] is a typical time-frequency analysis method to detect constant acceleration motion. With the changing of slow time, the acceleration produces a chirp rate modulation, and the velocity has a centroid frequency related to slow time in the echo phase. In this case, the echo can be treated as a linear frequency modulation signal (LFM) in slow time. Therefore, the time-frequency analysis method points to this feature to integrate echo in the time-frequency domain. By using the WVD time-frequency analysis method, the target echo can be well integrated. The modified WVD method has been studied in [
28,
29,
30]. Though the time-frequency analysis method can estimate and accumulate low SNR target signals, some limitations should be noticed. The time-frequency transform is non-orthogonal, and the side lobes cannot be ignored. In conclusion, these traditional methods share problems of model limitations and high computational complexity.
Except for signal domain processing, recently, the tracking-before-detection (TBD) methods [
4,
31,
32,
33,
34,
35,
36] have been studied, which combine target detection with tracking. Unlike traditional tracking-after-detection methods that use target estimation parameters to track, TBD methods use original signal data. Combined with the signal and target motion transition formulation, the energy can be integrated along the target’s motion tracks. After integration, the tracks can be outputted when the integration amplitude exceeds the CFAR threshold. There are four categories of TBD methods: dynamic programming (DP), the recursive Bayesian approach, the finite-set statistics theory (FISST) method, and the histogram-probabilistic multi-hypothesis tracker (H-PMHT). However, the TBD method usually uses incoherent integration along the target tracks, and thus, the technique loses performance when the single pulse is a low SNR condition.
Based on the weak target detection consideration, there are two critical problems: model mismatching and high computational complexity. The problem of model mismatch arises because the target cannot keep the same motion model in a long CPI, and it needs a higher-order model to fit the motion model. Given these problems, we propose a robust and efficient layer integration algorithm to detect weak and maneuvering radar targets with a long CPI for narrowband radar. This paper uses the piecewise way to simplify the complex motion over a long time. First, the echoes over a long time are divided into several echo segments over a short time. The complex motion over a long time can be regarded as combining several different uniform motions in a short time. Then, the coherent integration for every signal segment is carried out, and a low detection threshold is used to obtain the rough detection results in every signal segment. Because of every segment’s short coherent integration time, many false alarms are in the detection results. A target association and signal integration between segments using the H-PMHT are proposed to increase SNR gain and exclude false alarms. The high-order motion parameter is estimated based on the velocity variation between target association results. Because the target can be associated and the target’s signal can be integrated between segments, and as the false alarms are hardly integrated, the false alarms can be excluded. The target can be detected with the segment fusion layer by layer. In addition, to avoid the target glint problem, the target information transmission policy is also designed to prevent the loss of target information. Compared with traditional methods, the proposed method mainly has three advantages:
(1) It uses piecewise integration and TBD between slow-time segments to increase model adaption and decrease computational complexity;
(2) The layer integration detection mode, updated CFAR threshold, and target-alarm transmission association mechanism are designed so that the glint targets can be detected effectively;
(3) It uses coherent integration processing and compensates for the motion parameters which are estimated by associating between segments, helping to detect the low SNR targets effectively.
The rest of paper is organized as follows: In
Section 2, the algorithm framework is built, and in
Section 3, the method is introduced extensively, including the signal model, threshold detection, target association and target-alarm information base design. In
Section 4, the algorithm’s complexity and limitations are discussed. In
Section 5, some simulation experiments are designed. Finally, some conclusions are summarized in
Section 6.
2. Algorithm Framework
Figure 1 shows the flow chart of the method of layered target detection based on information transmission and the updated threshold, and
Figure 2 shows the layer integration mode. The algorithm can be summarized as the segment division that decreases the order of motion model to realize the coherent integration in short periods, which offers enough of an SNR to ensure the correctness of the TBD association. Furthermore, since the low detection threshold results in many false alarms, updating the increased threshold excludes the non-integration alarms. Moreover, the target-alarm transmission machine ensures that the glint target information can be transmitted with each layer increasing. The process is divided into five steps as follows:
(1). Preprocessing: The pulse compression and segment division are based on radar parameters in slow time. To detect all targets, we must select a certain number of pulses to obtain enough of an SNR and ensure that the target cannot have Doppler migration in every segment. The specific method of segment determination is discussed in the following sections.
(2). Noise power estimation and low-threshold detection: Range migration is corrected using the keystone transform in every slow-time segment. The noise power is estimated, and the CFAR threshold is set to detect targets based on noise power. To make sure that all targets can be detected, the threshold should be low, which means that many false alarms may be involved; thus, the increased threshold with layers is necessary to exclude false alarms and preserve the target.
(3). Target association and interaction: Combined with the TBD method, targets can be associated between segments. The target usually has the case of glint RCS, which means that the targets cannot be detected in a specific segment, and the target information cannot be transmitted into the following layer. The problem can be solved by designing a target-alarm information transmission mechanism to preserve targets and alarm information in every layer so that the information on the glint target would not be discarded as the layers increase.
(4). Compensation and integration: Compensate for the range and Doppler phase of the associated targets using the estimated parameters and coherently integrate them. Then, the integrated segments are transmitted to the next layer and processed using the same steps.
(5). CFAR threshold updating: With the increasing layers, the targets are integrated, and the alarms are not. Thus, the CFAR threshold needs to be increased based on the integration gain.
The algorithm was summarized above, and the specific algorithm is described in the next section.
3. Layer Integration and Detection Method with Information Transmission and Threshold Updating Algorithm Description
In
Section 2, we described our overall algorithm framework. In this section, we discuss the algorithm specifically using mathematical derivation from four aspects, which are the layer signal model, integration and detection, tracking and association, and the information transmission mechanism.
3.1. Layer Signal Model
Suppose that radar transmits an LFM signal:
where
is the rectangle window function,
is the pulse width,
is the radar’s carrier frequency, and
is the chirp rate. There are
targets with a high-order motion in the scenario, and the
target motion formulation can be expressed as:
where
is the slow time. The received signal of targets can be written as follows after down-conversion:
The echo in the range frequency-slow time domain
can be written as:
where
is the signal’s bandwidth, and
is the Fourier transform of the LFM signal. By multiplying Equation (4) by
(
denotes conjugate), the signal model is given by [
14]:
From Equation (5), we can see that the phase is nonlinear due to maneuverable motion. We use the piecewise concept to adapt to high-order complex motion models, which use several uniform motions in continuous slow time to approximate high-order motions. This way, we divide the long CPI into several segments, thus the phase can be approximately regarded as linear and integrated using the keystone transform. Moreover, the assumption of the piecewise segmental method is reasonable since narrowband radar usually has a low pulse repetition frequency (PRF) to reduce range ambiguity. Combined with dividing several segments of slow time, every segment’s velocity resolution is low. Thus, it can be regarded as a constant velocity with a low-velocity resolution. Moreover, the acceleration error can be compensated for by using the estimation parameter of the segment association. We use the long CPI to integrate signals to get a high SNR gain. However, it will bring substantial computational costs if we process the signal during the slow time. For example, suppose the keystone transform is used to correct the range migration. In that case, the computational complexity is proportional to the square of the slow-time number using the discrete Fourier transform (DFT) method. Therefore, if we use the long CPI to obtain coherent integration results, the computational burden usually is unacceptable. On the contrary, the segment division decreases the complexity affected by the number of pulses, and parallel computing can be used in every segment to increase calculation efficiency. Then, we derive the signal model based on the segment division.
Suppose that there are
pulses in one CPI, and the signal transmits to the
layer and has
segments. In this case, the number of pulses in the
layer is
, and the
segment’s slow time can be represented as:
where
is the pulse repetition period. Based on the piecewise consideration, the segment division should ensure that the acceleration motion cannot cause the Doppler migration in every segment. In this way, the additional phase caused by acceleration should be less than
. For the practical detection scene, the maximum acceleration has a fixed bound. For example, the acceleration of flight is less than 10
. Therefore, we make the target’s maximum acceleration value
, and the maximum phase caused by the acceleration can be represented as:
where
is the calculator to calculate the angle of the complex value and
is the pulse number in one segment. Thus, it can be equivalent to:
Further, in order to ensure the maximum energy integration, the pulse number of one segment can be determined by:
where
is rounded down to an integer. Substituting Equation (6) into Equation (5), we have:
Suppose that the radar maximum non-ambiguous velocity is
, the number of ambiguity is
, and the velocity without ambiguity is
. Then, the motion equation of each segment can be written as:
Thus, the signal model in the
segment of the
layer can be written as:
In this way, we use the segment division method to realize model matching and reduce the computational burden problem. Thus, the signal can be decoupled and integrated in every segment using the keystone transformation as subsection B.
3.2. Keystone Uncoupling and Updating CFAR Threshold
After segmentation division and piecewise approximation using the method in subsection A, the nonlinear phase can be regarded as a linear phase coupled with the range frequency in a short time period. The keystone transformation modifies the slow time axis to decouple the range frequency and slow time from the coupling phase; this specific method references [
11]. Note that the keystone transformation is performed without the ambiguous velocity; however, the velocity of the target is usually ambiguous in low frequency radar, thus the ambiguous velocity should be compensated in the different ambiguity channels. The number of ambiguity can be determined by the energy of the ambiguous channel output.
The ambiguous channel of the keystone signal model unfolding Equation (10) is:
where
is the velocity without ambiguity. The compensation function of the ambiguous velocity can be expressed as:
The signal after the ambiguity compensation is:
By applying the keystone transformation, and substituting
into Equation (13), the signal can be written as:
The keystone transformation can be processed by the chirp Z-transform-inverse Fourier transform (CZT-IFFT) or scaled Fourier transform-inverse Fourier transform (SFT-IFFT) method to improve operational efficiency. Therefore, the signal can be integrated by moving the target detection (MTD) after the keystone transformation.
After integrating the first layer, the CFAR detection threshold should be determined based on noise power. The high speed and maneuvering target detection are usually applied in detecting space targets such as missiles and aircraft. The space target detection usually has these characteristics: (1) the clutter has little effect on echo; (2) the echo is mainly affected by the thermal noise of the transmitter and receiver; (3) the distribution of noise with different ranges and Doppler cells is usually homogeneous. In this case, the influence of clutter on the detection is small, and only Gaussian white noise is considered. Assuming that all the measurement cells share the same noise distribution, based on the CFAR detection theory, the adaptive threshold corresponding with the constant false alarm rate and statistical noise power should be applied to the detector. As for this condition, the white Gaussian noise obeys
, and the amplitude obeys the Rayleigh distribution after envelope detection. The possibility density function is:
where
is the hypothesis that there are no targets in the detection cell and
is the amplitude of echo. Based on the Neyman–Pearson criterion, the alarm rate can be deduced as [
37]:
where
is the detection threshold. Thus, the detection threshold can be represented as
. As for the noise power estimation, we can choose the range cells in the far distance to estimate it, since the echo power is proportional to
, where
is the range between the cell and radar. The power of the target in the far distance is much lower than the noise and can be ignored. The noise power can be estimated by these cells as:
where
is the number of reference cells needed in order to estimate noise. Thus, the threshold can be determined by noise power estimation. It is noted that a low threshold should be used in the first layer to ensure that all the targets can be detected. However, the low threshold means a high false alarm rate and many false alarms involved. Updating the increased threshold is necessary to exclude false alarms. Assume that the threshold in the
layer is
, and the number of pulses in one segment of this layer is
. The target can obtain a
dB gain of power by coherent integration and the noise obtains a
dB gain of power. Based on the difference in the integration gain, the threshold power can be increased by
dB to distinguish the target and false alarm. Therefore, the threshold in
can be determined by:
In this way, with the layer increasing, the targets can exceed the increased threshold and false alarms cannot since the false alarms have a low probability of being associated with the layer detector passing by; therefore, the false alarms will be excluded and the target will be integrated and detected in the last layer.
3.3. H-PMHT Association and Integration
After segment division, integration, and target detection in every segment, the signals need to be integrated between segments. The targets have a high-speed and maneuverable motion, resulting in the range and velocity change in the adjacent segment. Based on the detection and integration result, the problem becomes associating targets between segments, compensating for the phase caused by the maneuverable motion, and integrating them. Combined with the TBD concept, we use the H-PMHT method to associate targets. The H-PMHT method quantizes the amplitude of measurement cells as a synthetic histogram with several shots. The shot is treated as a measurement, and the sum of shots is treated as the total number of acquiring measurements. The shot distribution of a segment obeys multinomial distribution, where the probability of each cell can be considered as a superposition of targets and noise. Moreover, successive scans are associated with the state transformation function. Thus, the joint probability distribution can be obtained, and parameters can be estimated using the expectation maximization (EM) method. Moreover, [
37,
38] derived that the Kalman filter can be used to estimate parameters, including motion parameters and amplitude, in the subsequent scans, whereas [
39,
40] revealed that the estimated intensity can indicate whether or not the association is successful. The target will be treated as an alarm if the estimated intensity drops below 0 dB.
Suppose that targets have constant acceleration between segments; then, the state transformation function in the
layer can be expressed as:
where the
and
target states in the segment are
, and
, which represent the range, velocity, and acceleration observation value, respectively.
is the state error. The acceleration can be denoted by
. The measurement function can be stated as:
Based on the H-PMHT, the distribution of the mean cell value is given by [
40]:
where
is the state distribution of the
target, and when
, it is noise distribution. We assume that the noise distribution is uniform in the detection cells, thus the pdf of noise is:
where
and
are the number of range cells and number of Doppler cells. As for target distribution, it is reasonable that we assume the range and velocity distribution obeys the normal density function with the variance
and
. Thus, the
target contribution can be expressed as [
39]:
The state and proportion estimations can be updated by using recursive implementation. By combining the intensity and state estimations, we can confirm whether or not the targets are associated, and the specific algorithm is summarized as Algorithm [
40].
Targets have been associated between segments
using the TBD method. The next step is to compensate for the range and velocity migration, and to integrate the signal between segments. Suppose that the target state in two segments can be represented as
and
, the signal in two segments can be derived as follows:
Thus, the compensation function can be written as:
After compensating for the motion phase of targets based on Equation (23), the signal can be coherently integrated.
3.4. Target Transmission Mechanism
The signal can be integrated using the above method in long coherent periods; however, there is an additional question to consider regarding glint target integration and detection. The glint target usually has a low RCS in a specific segment, and the echo cannot be integrated effectively. Based on the layer detection framework, the false alarms which cannot exceed the increased threshold will be discarded and result in missing the detection of the glint target, which is not integrated and treated as a false alarm. Considering this question, if we preserve the false alarm information and design the mechanism of information transmission and interaction, the glint targets could be detected. From the target detection analysis,
Figure 3 shows all the possible situations after association and integration.
According to the above analysis, false alarms are unlikely to be associated between segments. Even though they can be successfully associated with a specific layer, they are not likely to pass by the final layer and exceed the threshold. Therefore, only the cases of the targets’ association and detection should be considered to avoid missing detection. The association and detection of targets can be divided into four issues as follows:
(1) Targets are well associated and exceed the threshold (normal targets). It is an ideal condition, and the target can be detected in the last layer;
(2) Targets are well associated and do not exceed the threshold (weak targets). The target is unstable and the integration gain is less than the increased threshold, thus the target’s information should be preserved and transmitted into the following layers;
(3) Targets are not associated and exceed the threshold (strong glint targets). Although the target is glint, the power is strong enough in a specific segment, and the target also can be detected in the following layers;
(4) Targets are not associated and do not exceed the threshold (glint targets). If the target is glint and weak, the target will not be associated and exceeds the threshold. Thus, it cannot be detected in the following layer, and the information should be preserved.
As for these four types of targets, the target and false alarm pursuit mechanism is designed as shown in
Figure 4. In short, the alarm information will be preserved, and the association steps between targets and false alarms will be added to complete the association of type (2) and (4) targets. In a nutshell, the association method is designed with three stages: (a) associate between detected targets, which ensures that normal targets can be associated; (b) non-associated targets associate with false alarms, which makes sure that weak or glint targets can be integrated; and (c) associate between false alarms, which makes sure that the glint target which is not integrated and detected can be associated. This way, the four kinds of targets are associated and transmitted into the final layer to be detected. The specific processing method refers to
Figure 4.