Next Article in Journal
A Tunable Hyperspectral Imager for Detection and Quantification of Marine Biofouling on Coated Surfaces
Next Article in Special Issue
Signal Property Information-Based Target Detection with Dual-Output Neural Network in Complex Environments
Previous Article in Journal
A Lightweight YOLOv5-MNE Algorithm for SAR Ship Detection
Previous Article in Special Issue
Radar Target Detection Algorithm Using Convolutional Neural Network to Process Graphically Expressed Range Time Series Signals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Weak and Maneuvering Target Detection with Long Observation Time Based on Segment Fusion for Narrowband Radar

1
College of Oceanography and Space Informatics, China University of Petroleum (East China), Qingdao 257099, China
2
Beijing Institute of Radio Measurement, Beijing 100854, China
3
Beijing Institute of Tracking Telecommunications Technology, Beijing 100094, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(18), 7086; https://doi.org/10.3390/s22187086
Submission received: 8 August 2022 / Revised: 12 September 2022 / Accepted: 14 September 2022 / Published: 19 September 2022
(This article belongs to the Special Issue Airborne Distributed Radar Technology)

Abstract

:
Detecting high-speed and maneuvering targets is challenging in early warning radar applications. Modern early warning radar has many functions such as detection, tracking, imaging, and recognition which need a high signal-to-noise ratio (SNR). Thus, long-time coherent integration is a necessary method to realize high SNR requirements. However, high-speed and maneuverable motion cause range and Doppler migration, which brings about serious coherent integration loss. Traditional integration methods usually have the drawbacks of model mismatching and high computational complexity. This paper establishes a novel long coherent processing interval (CPI) integration algorithm that detects maneuvering and weak targets which have a low reflection cross-section (RCS) and low echo SNR. The range and Doppler migration problems are solved via a layer integration by blending the association in a tracking-before-detection (TBD) technique. Compact SNR gain is achieved with a target information transmission mechanism and an updated constant false alarm ratio (CFAR) threshold. The algorithm is applicable in multiple target scenarios by considering different velocity ambiguities and maneuvers. A simulation and real-measured experiments confirm the effectiveness of the algorithm.

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:
s ( t ) = rect t T p exp j 2 π f c t + π γ t 2
where r e c t [ ] is the rectangle window function, T P is the pulse width, f c is the radar’s carrier frequency, and γ is the chirp rate. There are K targets with a high-order motion in the scenario, and the k - t h target motion formulation can be expressed as:
R k ( τ ) = R 0 k v k τ a k τ 2 h k τ 3 ,
where τ is the slow time. The received signal of targets can be written as follows after down-conversion:
s r ( t , τ ) = k = 1 K σ k rect t R k ( τ ) T p exp j 4 π f c R k ( τ ) c exp j π γ t 2 R k ( τ ) c 2 .
The echo in the range frequency-slow time domain ( f , τ ) can be written as:
s r , c ( f , τ ) = P ( f ) k = 1 K σ k rect f B exp j 4 π ( f + f c ) R k ( τ ) c
where B is the signal’s bandwidth, and P ( f ) is the Fourier transform of the LFM signal. By multiplying Equation (4) by P f * ( * denotes conjugate), the signal model is given by [14]:
s r , c ( f , τ ) = P ( f ) 2 k = 1 K σ k rect f B exp j 4 π ( f + f c ) R k ( τ ) c .
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 N pulses in one CPI, and the signal transmits to the l - th layer and has n l segments. In this case, the number of pulses in the l layer is N l = N n l , and the q - th segment’s slow time can be represented as:
τ ( l , q , p ) = ( q 1 ) N l T r + p T r , p = 1 , 2 , 3 N l ,
where T r 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 π / 4 . For the practical detection scene, the maximum acceleration has a fixed bound. For example, the acceleration of flight is less than 10 m / s 2 . Therefore, we make the target’s maximum acceleration value a max , and the maximum phase caused by the acceleration can be represented as:
a n g l e exp j 4 π c f c 1 2 a max N 2 Δ T 2 < π 4
where a n g l e is the calculator to calculate the angle of the complex value and N is the pulse number in one segment. Thus, it can be equivalent to:
1 2 a max N 2 Δ T 2 < λ 16 N < λ 2 a max Δ T 2
Further, in order to ensure the maximum energy integration, the pulse number of one segment can be determined by:
N = λ 2 a max Δ T 2
where is rounded down to an integer. Substituting Equation (6) into Equation (5), we have:
s r , c f , l , q , p = P ( f ) 2 k = 1 K σ k rect f B exp j 4 π f + f c R k τ l , p , q c .
Suppose that the radar maximum non-ambiguous velocity is v max , the number of ambiguity is n v , and the velocity without ambiguity is v 0 . Then, the motion equation of each segment can be written as:
v k ( l , q ) = n v v max + v 0 k + a k ( q 1 ) N l T r .
R k l , q , p = R 0 k l , q v k l , q p T r .
Thus, the signal model in the q - th segment of the l - th layer can be written as:
s r , c ( f , l , p , q ) = P ( f ) 2 k = 1 K σ k rect f B exp j 4 π ( f + f c ) R k τ ( l , p , q ) c .
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:
s r , c f , l , p , q = P ( f ) 2 k = 1 K σ k rect f B exp j 4 π ( f + f c ) R 0 k , l , q n v v max + v 0 k l , q p T r λ ,
where v 0 k l , q is the velocity without ambiguity. The compensation function of the ambiguous velocity can be expressed as:
H c ( f , n v ) = exp j 4 π λ ( 1 + f f c ) ( n v v max p T r ) .
The signal after the ambiguity compensation is:
s ( 1 ) r , c ( f , l , p , q ) = s r , c ( f , l , p , q ) H c ( f , n v ) = k = 1 K σ k rect f B exp j 4 π ( f + f c ) R 0 k , l , q v 0 k ( l , q ) p T r c .
By applying the keystone transformation, and substituting ( f + f c ) p T r = f c p T r into Equation (13), the signal can be written as:
s ( 1 ) r , c f , l , p , q = P ( f ) 2 k = 1 K σ k rect f B exp j 4 π ( f + f c ) R 0 k , l , q c exp j 4 π v 0 k ( l , q ) p T r λ
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 n ( t ) ~ N ( 0 , σ 2 ) , and the amplitude obeys the Rayleigh distribution after envelope detection. The possibility density function is:
f A | H 0 = A σ 2 exp A 2 2 σ 2 ,
where H 0 is the hypothesis that there are no targets in the detection cell and A is the amplitude of echo. Based on the Neyman–Pearson criterion, the alarm rate can be deduced as [37]:
P f a = V T A σ 2 exp A 2 2 σ 2 d A = exp V T 2 2 σ 2 ,
where V T is the detection threshold. Thus, the detection threshold can be represented as V T = 2 σ 2 ln 1 P f a . 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 1 / r 4 , where r 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:
σ ^ 2 = 1 W w = 1 W s w 2 ,
where W 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 l t h layer is V T l , and the number of pulses in one segment of this layer is N l . The target can obtain a 10 log 10 N l 2 dB gain of power by coherent integration and the noise obtains a 10 log 10 N l dB gain of power. Based on the difference in the integration gain, the threshold power can be increased by 10 log 10 N l 2 dB to distinguish the target and false alarm. Therefore, the threshold in l + 1 t h can be determined by:
V T l + 1 = 10 10 log 10 V T l 2 + 10 log 10 N l 2 .
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 l - th layer can be expressed as:
x p + 1 = 1 N l T r ( N l T r ) 2 2 0 1 N l T r 0 0 1 x p + ε ,
where the p - th and p + 1 - th target states in the segment are x p = r p , v p , a p T , x p + 1 = r p + 1 , v p + 1 , a p + 1 T , and r , v , a , which represent the range, velocity, and acceleration observation value, respectively. ϵ is the state error. The acceleration can be denoted by a p = a p + 1 = v p + 1 v p N l T r . The measurement function can be stated as:
y p + 1 = 1 0 0 0 1 0 x p + 1 .
Based on the H-PMHT, the distribution of the mean cell value is given by [40]:
P ( r , v ) = k = 0 K π k P k ( r , v ) ,
where P k ( r , v ) is the state distribution of the k t h target, and when k = 0 , it is noise distribution. We assume that the noise distribution is uniform in the detection cells, thus the pdf of noise is:
P 0 = 1 N R N V ,
where N R and N V 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 σ r and σ v . Thus, the k - th target contribution can be expressed as [39]:
P k ( r , v ) = N ( R k , σ r ) N ( V k , σ v ) , k = 1 , 2 K .
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 ( q , q + 1 ) 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 ( r 1 , v 1 ) and ( r 2 , v 2 ) , the signal in two segments can be derived as follows:
s r 1 , v 1 = σ k rect f B exp j 4 π f + f c r 1 c exp j 4 π v 1 ( q n l T r + p T r ) λ s r 2 , v 2 = σ k rect f B exp j 4 π ( f + f c ) r 2 c exp j 4 π v 2 ( q + 1 ) n l T r + p T r λ .
Thus, the compensation function can be written as:
H r = exp j 4 π ( f + f c ) ( r 1 r 2 ) c H v = exp j 4 π ( v 1 v 2 ) ( q + 1 ) n l T r + p T r λ H a = exp j 4 π a q + 1 n l T r + p T r 2 λ .
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.

4. Algorithm Complexity and Parameter Determination

4.1. Computational Complexity

As for computational complexity, the LVD method is an effective way to solve the maneuvering motion integration problem. It needs a correlation process and keystone transformation in every range cell. Suppose that N r is the number of range cells, N v is the number of velocity cells, and M is the number of ambiguous channels. The computational complexity is O N r M N v 2 log N v [20]. The proposed method needs to implement one keystone transformation in the first layer. Moreover, it does not add the signal dimension, and the number of targets and alarms determines the complexity of the association. Suppose the number of targets is K , the number of segments is n l , the number of false alarms has an expectation N v N r n l P f a and N a is the number of iterations; thus, the proposed method’s computational complexity is less than:
O M N r N v log N v + K + N r N v n l P f a N a log 2 n l .
Therefore, our method is more efficient compared with the LVD distribution method. Compared with parameter searching methods such as the Radon Fourier transform [16], our approach only needs to search ambiguous parameters. Thus, many repetition complements are omitted to increase calculation efficiency. Based on the complexity analysis, the radar parameters are set as follows: carrier frequency—1 GHz; bandwidth—15 MHz; pulse repetition frequency—2000 kHz; the number of pulses—960; the number of range cells—100; and false alarm rate— 10 7 . In this case, we set a single target scenario and used an Intel® CoreTM i7-7700 CPU and Matlab 2017 to compare the proposed method with the LVD method for calculation efficiency. The LVD method uses 31.779885 s on average to calculate one range cell. However, the total calculation time of the proposed method is 5.344031 s. The results show that the proposed method is more efficient than the correlation method, such as LVD.

4.2. Parameter Determination

A. The threshold of the first layer determination. The first layer threshold should be determined by the lowest signal power in a single pulse. Assume that the lowest SNR of a target after pulse compression to be detected is P 0 min dB, and the number of integration pulses in the first layer is N 1 . In this case, the target’s power after coherent integration of the first layer is P 0 min + 10 log 10 N 1 . If the target with minimum power can be detected in the first layer detector, the threshold V T 1 should meet the requirement as:
V T 1 10 P 0 min + 10 log 10 N 1   .
B. The number of pulses in the first layer segment determination. The number of pulses in the first layer is determined by two aspects: Doppler migration and the SNR requirement in the first layer detector. The Doppler migration limitation makes sure that the target echo satisfies the piecewise assumption. The SNR requirement makes sure that the target which has the lowest power requirement can be detected. Based on these two criterions, the number of pulses can be determined as follows:
(1) Doppler migration limitation. Based on the uniform motion assumption, the target cannot cross the Doppler cell in the segment. The Doppler migration is caused by an equivalent acceleration, and the maximum equivalent acceleration is defined as a max in the realistic application. Assuming that the number of pulses is N 1 in the first layer, the pulse repetition frequency is f r , and the wave length is λ , the target with the maximum acceleration should satisfy the condition:
a max N 1 f r f r λ 2 N 1
Simplify it, and the limitation can be derived as:
N 1 λ f r 2 2 a max .
(2) The SNR requirement in the first layer detection. The SNR can be increased by 10 log 10 N 1 dB after N 1 coherent integration pulses in the first layer. Assuming that the lowest SNR of a target that can be detected in the first layer is Q 1 min dB, then the minimum SNR of the single pulse Q 0 min should meet the requirement:
Q 0 min Q 1 min 10 log 10 N 1 .
Therefore, the N 1 can be determined as:
N 1 10 Q 1 min Q 0 min / 10 .
The number of pulses N 1 in the first layer segment can be determined by the above two limitations. However, there is a situation where the N 1 based on the SNR criterion is bigger than the Doppler criterion. In this case, we should increase the motion order in the segment, such that we assume that the target has constant acceleration in every segment and has constant jerk between segments. In this way, the number of pulses N 1 can be increased to make sure that it satisfies the Doppler and SNR criterions at the same time.

5. Experiment Analysis

The experiment comprises three parts to prove that our method is effective for multiple targets, especially glint target integration and detection. First, a scenario of multiple targets with different range velocity accelerations and ambiguity numbers is set to prove that our method is valid. Second, a glint target scenario is set to verify that our approach has better performance for glint target integration compared with the traditional method. An integration comparison experiment uses our technique and traditional techniques to analyze integration and detection performance. Third, Monte Carlo experiments are designed to quantitatively compare the detection probability and integration SNR.

5.1. Integration SNR Upper Limit

In order to make a reference to integration performance, first we derive the upper limit of the pulse compression and coherent integration signal. Suppose that the length of transmitting sampling is L s , and the RCS of the target is A t ; thus, the signal’s power P s after pulse compression is:
P s = L s A t 2 .
We assume that the noise power spectral density is P n and the SNR after pulse compression can be derived.
S N R P C = P s P n = 10 log 10 L s A t 2 P n .
If we coherently integrate the N pulse compression signal, the ideal integration SNR can be written as:
S N R o u t i d e a l = S N R P C + 10 log 10 ( N )
Whatever integration method we use, the integration SNR output cannot exceed Equation (34); thus, we can use the SNR upper limit to evaluate the integration method performance.

5.2. Case of Multiple Target Integration and Detection

To prove that our algorithm is effective in the aspects of multi-targets, multi-ambiguity channels, glint targets, and low SNR detection, the scenario that includes three targets is set to evaluate the multi-target detection performance. The ranges of target1, target2, and target3 are in the 2000th, 1000th, and 300th cells, and the velocities are 600 m/s, 1600 m/s, and 2300 m/s, respectively. The accelerates of the targets are 60 m/s2, 60 m/s2 and 20 m/s2, respectively. Based on the radar parameters shown in Table 1, target1 is in the zeroth ambiguity channel, and target2 and target3 are in the first ambiguity channel. In this way, we can verify the ability of the multi-ambiguity channel processing. Third, target1 is the glint target, whose RCS is too low to be detected in the segments. As for the target glint, it is assumed that the glint segment obeys binomial distribution. The glint probability is set as 0.25. The glint target integration and detection ability can be tested in this case. Finally, we set the single-pulse SNR to be −33 dB to verify the low SNR detection performance. Using the above target scenario to test the integration and detection performance, this experiment can test the performance of multi-targets, multi-ambiguity channels, glint targets, and low SNR conditions. Figure 5 shows the results of echo pulse compression. Before the coherent integration, the tracks of targets cannot be seen in Figure 5 under the low SNR condition.
Using the proposed method of integration and detection to obtain the visualized integration and detection result in every layer, a three layer detector is set. In other words, the original signal can be divided into four segments. Figure 6, Figure 7 and Figure 8 show the integration and detection results of every layer in the zeroth ambiguity channel, whereas Figure 9, Figure 10 and Figure 11 show the integration and detection results of every layer in the first ambiguity channel. From Figure 6, only target1 energy can be effectively integrated after using an ambiguity match filter. Moreover, due to the target glint, the target1 cannot be combined and detected in the second segment of the first layer, as shown in Figure 6b, but the target1 in the rest of the segments of the first layer can be detected, as shown in Figure 6a,c,d. Therefore, the target cannot be associated with the first and second segments in the first layer. From Figure 7a, as the threshold increases, the target cannot be detected in the first segment of the second layer. Thus, the target is treated as a false alarm. Although target1 can be detected in the second segment of the second layer, the information of target1 in the first segment of the second layer will be lost. However, due to the design of the target and alarm information transmission mechanism, the information of target1 can be transmitted successfully. The target1 in the first segment can be associated with the second layer detector. After parameter estimation and compensation, the target1 is integrated with the last layer and exceeds the previous detection threshold to be detected successfully, as shown in Figure 8. If the transmission mechanism is not designed, the information of target1 will be discarded in the second layer, which results in missing detection. As for the integration and detection result in the first ambiguity channel, from Figure 9, Figure 10 and Figure 11, target2 and target3 without glint are in the same ambiguity channel. After ambiguity matching, they can be well integrated and detected in the first layer, as shown in Figure 9. The motion formulation can be determined to realize the association based on the detection result. The two targets are well associated due to precise parameter estimation and are transmitted to the following layers, as shown in Figure 10. In this case, the targets without glint can be integrated and detected in the last layer from Figure 11. In conclusion, based on the experiment results, the ambiguity channel matching filter solves the problem of high-speed ambiguity. The piecewise segment division realizes integration in a short period, and the association between segments ensures that motion parameters can be determined and the acceleration motion can be compensated for. By transmitting with layers, the targets can be well integrated. Moreover, the information transmission and perseveration mechanism can ensure that the glint target is associated with the segment in all layers. Due to one experiment having a high random error, the quantitative integration and detection performance is not discussed in this subsection. The specific integration performance of glint and non-glint targets will be analyzed in Section 3 and Section 4.

5.3. Case of Glint Target Integration Performance Analysis

The target-alarm information transmission mechanism is designed as above. In this subsection, we design a comparison experiment with the LVD distribution method to prove that the proposed method has good integration performance with the glint target. Suppose that the target has a weak RCS in a specific segment, the glint case obeys binomial distribution with a probability of 0.25, and the radar parameters are set as shown in Table 1. The target is in the 2000th range cell. The velocity is 450 m/s, and the acceleration is 50 m/s2. A random simulation shows the target glint in the 251st–500th slow-time period. The proposed method and LVD distribution method integration results are shown in Figure 12 and Table 2.
From Figure 12, it can be seen that the glint target can be well integrated using the proposed method. The integration case with the LVD method can also realize coherent integration, but the integration SNR is much lower than the proposed method. Moreover, the LVD method has higher side lobes than the proposed method, resulting in the missing detection of a weak target in the side lobe area. As for the quantitative analysis shown in Table 2, the result shows that the integration SNR loss is 0.91 dB using the proposed method compared with LVD, which has a 11.8 dB SNR loss. This experiment proves that the proposed method performs better in glint target integration and detection. The reasons might be as follows: (1) The proposed method uses segment processing, and the signal is divided into several segments to be coherently integrated into every segment. It makes sure that the segment without glint can be sufficiently integrated, but LVD uses a correlation operation and decreases the SNR. (2) The proposed method uses an information transmission mechanism to be adapted to the target glint model, which makes sure that the target in every segment can be transmitted into all layers. (3) As for the LVD method, it uses a WVD and keystone transform to complete the correlation operation. The missing signal causes a significant effect on the WVD operation due to model mismatching. Based on the above reasons, the integration SNR has a considerable loss by using LVD, and the proposed method has a more robust ability for model matching.

5.4. Case of Quantitative Performance

In this subsection, we analyze a quantitative experiment’s integration and detection performance because the proposed method mainly solves the problem of SNR integration loss and the high complexity of the correlation method. In addition, the proposed method uses the TBD concept to realize association and integration. Therefore, it is necessary to compare the proposed method with traditional correlation and TBD detection methods. This subsection uses the traditional correlation method of LVD and the traditional TBD method of DP as comparative methods. The radar parameters are set as shown in Table 1. Suppose that the target velocity is 450 m/s, the acceleration is 50 m/s2, and 100 Monte Carlo trails are implemented. The results of output SNRs with different input SNRs are shown in Figure 13. The integration up bound can be easily determined using Equations (32)–(34), shown as the blue line in the figure. The SNRs of a single pulse before the pulse compression are set from −40 to 0 dB to test the output signal SNR using different methods. From the experiment results, the proposed method has the highest SNR output compared with other methods. The SNR integration loss is less than 1 dB compared with the upper limit. The LVD method has an integration SNR loss of about 3 dB compared with the integration up bound. The DP method has the worst performance because it uses the incoherent integration method. Our method has a better performance compared with other methods due to the reasons as follows: The LVD method uses the correlation between slow times, and the noise correlation causes the SNR loss. Moreover, the cross-terms in the correlation calculation bring about a high side lobe. In this case, the output SNR gain has a performance loss. On the other hand, the LVD method ignores the range migration caused by acceleration. Thus, it will cause model mismatching and result in SNR loss; DP uses the motion formulation to search the estimation area based on transmission states, and it uses incoherent integration to explore and determine the targets’ tracks when the final integration result is compared with the detection threshold. In conclusion, the estimation precision of motion parameters and incoherent integration limits the performance of the DP method. Therefore, the incoherent integration gain is much lower than a coherent integration method such as LVD and the proposed method, as the proposed method uses a piecewise concept to minimize model error. Moreover, the association and integration between segments solve the shortage of SNR loss, which comes from correlation. On the other hand, the proposed method combines TBD tracking with coherent integration and solves the problem of low SNR gain using the traditional DP method.

5.5. Real-Measured Data Analysis

The real-measured data of the space target are used to verify that the proposed method effectively completes signal integration for weak targets. The radar parameters are listed in Table 3.
The recorded echoes are shown in Figure 14. The original echoes have a very low SNR and we almost cannot find the target’s tracks in Figure 14.
For this target, we use the proposed detection method of segment division to realize coherent integration. The detection false alarm rate in the first layer is set as 10 2 and the false alarm rate in the last layer is set as 10 6 , and we output the integration results of every layer, which is shown in Figure 15.
From the results, we can see that the energy of the target is integrated layer by layer, and the false alarm cannot be integrated and excluded. In order to verify that the target is integrated gradually in every layer, the statistical SNR is calculated in Figure 16. From the results, we can see that the proposed method successfully realizes target integration and detection for the space target in real-measured data.
Finally, the estimation results of range, velocity, and acceleration are also given in Table 4. In all, based on the real-measured results, the proposed method can be used in practical applications.

6. Conclusions

Using the piecewise approximation mindset and segment integration method, the proposed method has advantages in increasing SNR gain and computational efficiency for small and weak targets. As seen in the simulated and real-measured experiments, the proposed method has more than 1 dB of SNR integration gain compared with traditional methods such as LVD. In addition, the computational speed increases almost 6 times more than in LVD. However, the proposed method relies on the robust detection results in the first layer, and a more complex motion will lead to the performance degradation of the proposed method. In addition, the target association is not considered good; thus, we will focus on robust initial detection and target association in future works.

Author Contributions

Conceptualization, S.W., Y.D. and Q.Z.; methodology, S.W., Y.D. and Q.Z.; software, S.W., Y.D. and Q.Z.; validation, S.W., Y.D. and Q.Z.; formal analysis, S.W., Y.D. and Q.Z.; writing—original draft preparation, S.W. and Q.Z.; writing—review and editing, S.W. and Q.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rao, X.; Tao, H.; Su, J.; Xie, J.; Zhang, X. Detection of constant radial acceleration weak target via IAR-FRFT. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 3242–3253. [Google Scholar] [CrossRef]
  2. Xu, J.; Xia, X.; Peng, S.; Yu, J.; Peng, Y.; Qian, L. Radar Maneuvering Target Motion Estimation Based on Generalized Radon-Fourier Transform. IEEE Trans. Signal Process. 2012, 60, 6190–6201. [Google Scholar]
  3. Li, Y.; Zhang, X.; Wang, G.; Li, L. A Novel HT-TBD Detection Approach for Near-Space Target. In Proceedings of the 2017 3rd IEEE International Conference on Computer and Communications (ICCC), Chengdu, China, 13–16 December 2017; pp. 1720–1724. [Google Scholar]
  4. Su, J.; Xing, M.; Wang, G.; Bao, Z. High-speed multi-target detection with narrowband radar. IET Radar Sonar Navig. 2010, 4, 595–603. [Google Scholar] [CrossRef]
  5. Carlson, B.D.; Evans, E.D.; Wilson, S.L. Errata: Search radar detection and track with the Hough transform. IEEE Trans. Aerosp. Electron. Syst. 2003, 39, 382–383. [Google Scholar] [CrossRef]
  6. Du, P.; Wang, Y.; Sun, W. Radar Weak Target Detection and Location Based on Hough Transform. In Proceedings of the 6th International Conference on Signal Processing, Beijing, China, 26–30 August 2002; pp. 1480–1483. [Google Scholar]
  7. Perry, R.P.; DiPietro, R.C.; Fante, R.L. SAR imaging of moving targets. IEEE Trans. Aerosp. Electron. Syst. 1999, 35, 188–200. [Google Scholar] [CrossRef]
  8. Zhu, D.; Li, Y.; Zhu, Z. A Keystone Transform Without Interpolation for SAR Ground Moving-Target Imaging. IEEE Geosci. Remote Sens. Lett. 2007, 4, 18–22. [Google Scholar] [CrossRef]
  9. Li, Y.; Zeng, T.; Long, T.; Wang, Z. Range Migration Compensation and Doppler Ambiguity Resolution by Keystone Transform. In Proceedings of the 2006 CIE International Conference on Radar, Shanghai, China, 16–19 October 2006; pp. 1–4. [Google Scholar]
  10. Huang, P.; Liao, G.; Yang, Z.; Xia, X.G.; Ma, J.; Zheng, J. Ground Maneuvering Target Imaging and High-Order Motion Parameter Estimation Based on Second-Order Keystone and Generalized Hough-HAF Transform. IEEE Trans. Geosci. Remote Sens. 2017, 55, 320–335. [Google Scholar] [CrossRef]
  11. Huang, P.; Liao, G.; Yang, Z.; Xia, X.; Ma, J.; Ma, J. Long-Time Coherent Integration for Weak Maneuvering Target Detection and High-Order Motion Parameter Estimation Based on Keystone Transform. IEEE Trans. Signal Process. 2016, 64, 4013–4026. [Google Scholar] [CrossRef]
  12. Kong, L.; Li, X.; Cui, G.; Yi, W.; Yang, Y. Coherent Integration Algorithm for a Maneuvering Target with High-Order Range Migration. IEEE Trans. Signal Process. 2015, 63, 4474–4486. [Google Scholar] [CrossRef]
  13. Zheng, J.; Su, T.; Liu, H.; Liao, G.; Liu, Z.; Liu, Q.H. Radar High-Speed Target Detection Based on the Frequency-Domain Deramp-Keystone Transform. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 285–294. [Google Scholar] [CrossRef]
  14. Sun, G.; Xing, M.; Xia, X.; Wu, Y.; Bao, Z. Robust Ground Moving-Target Imaging Using Deramp Keystone Processing. IEEE Trans. Geosci. Remote Sens. 2013, 51, 966–982. [Google Scholar] [CrossRef]
  15. Carretero-Moya, J.; Gismero-Menoyo, J.; Asensio-Lopez, A.; Blanco-del-Campo, A. A Coherent Radon Transform for Small Target Detection. In Proceedings of the 2009 IEEE Radar Conference, Pasadena, CA, USA, 4–8 May 2009; pp. 1–4. [Google Scholar]
  16. Chen, X.; Guan, J.; Zhao, Z.; Yonghua, X. Radar Maneuvering Target Detection Based on Scaling Processing and Radon-Fourier Transform. In Proceedings of the IET International Radar Conference 2015, Hangzhou, China, 14–16 October 2015; pp. 1–4. [Google Scholar]
  17. Yu, J.; Xu, J.; Peng, Y.; Xia, X. Radon-Fourier Transform for Radar Target Detection (III): Optimality and Fast Implementations. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 991–1004. [Google Scholar] [CrossRef]
  18. Qian, L.; Xu, J.; Xia, X.; Sun, W.; Long, T.; Peng, Y. Fast implementation of generalised Radon-Fourier transform for manoeuvring radar target detection. Electron. Lett. 2012, 48, 1427–1428. [Google Scholar] [CrossRef]
  19. Chen, X.; Guan, J.; Liu, N.; He, Y. Maneuvering Target Detection via Radon-Fractional Fourier Transform-Based Long-Time Coherent Integration. IEEE Trans. Signal Process. 2014, 62, 939–953. [Google Scholar] [CrossRef]
  20. Lv, X.; Bi, G.; Wan, C.; Xing, M. Lv’s Distribution: Principle, Implementation, Properties, and Performance. IEEE Trans. Signal Process. 2011, 59, 3576–3591. [Google Scholar] [CrossRef]
  21. Li, Y.; Su, T.; Zheng, J.; He, X. ISAR Imaging of Targets with Complex Motions Based on Modified Lv’s Distribution for Cubic Phase Signal. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 4775–4784. [Google Scholar] [CrossRef]
  22. Li, X.; Cui, G.; Yi, W.; Kong, L. Coherent Integration for Maneuvering Target Detection Based on Radon-Lv’s Distribution. IEEE Signal Process. Lett. 2015, 22, 1467–1471. [Google Scholar] [CrossRef]
  23. Zhang, Y.; Cheng, Q.; Sun, J.; Zhang, H. Near Space Maneuvering Target Detection Based on Polynomial Fourier Transform and Lv’s Distribution. In Proceedings of the 2017 10th International Congress on Image and Signal Processing, BioMedical Engineering and Informatics (CISP-BMEI), Shanghai, China, 14–16 October 2017; pp. 1–5. [Google Scholar]
  24. Tian, J.; Cui, W.; Lv, X.L.; Wu, S.; Hou, J.G.; Wu, S.L. Joint estimation algorithm for multi-targets’ motion parameters. IET Radar Sonar Navig. 2014, 8, 939–945. [Google Scholar] [CrossRef]
  25. Yinfang, M.; Guoan, C.; Junfeng, W. SAR/ISAR Imaging of Multiple Moving Targets Based on Combination of WVD and HT. In Proceedings of the International Radar Conference, Beijing, China, 8–10 October 1996; pp. 342–345. [Google Scholar]
  26. Wang, H.; Jiang, Y. Real-time parameter estimation for SAR moving target based on WVD slice and FrFT. Electron. Lett. 2018, 54, 47–49. [Google Scholar] [CrossRef]
  27. Huang, P.; Liao, G.; Yang, Z.; Xia, X.; Ma, J. Ground Moving Target Imaging and Motion Parameter Estimation Using Radon-Second-Order WVD Transform. In Proceedings of the 2015 IEEE 5th Asia-Pacific Conference on Synthetic Aperture Radar (APSAR), Singapore, 1–4 September 2015. [Google Scholar]
  28. Xing, M.; Wu, R.; Li, Y.; Bao, Z. New ISAR imaging algorithm based on modified Wigner-Ville distribution. IET Radar Sonar Navig. 2009, 3, 70–80. [Google Scholar] [CrossRef]
  29. Lan, Y.; Li, Z.; Li, W.; Yang, Q.; Wu, J.; Yang, J.; Yang, H.; Yang, X. Maritime Moving Target Detection with Space-based Passive Radar Based on Multi-frame Modified Wigner-Ville Distribution. In Proceedings of the EUSAR 2018, 12th European Conference on Synthetic Aperture Radar, Aachen, Germany, 4–7 June 2018; pp. 1–4. [Google Scholar]
  30. Huang, P.; Liao, G.; Yang, Z.; Xia, X.; Ma, J.; Zhang, X. A Fast SAR Imaging Method for Ground Moving Target Using a Second-Order WVD Transform. IEEE Trans. Geosci. Remote Sens. 2016, 54, 1940–1956. [Google Scholar] [CrossRef]
  31. Hongbo, Y.; Guohong, W.; Qian, C. A Novel Unscented Filter TBD Algorithm FOR Weak Radar Target. In Proceedings of the 2013 IEEE International Conference on Signal Processing, Communication and Computing (ICSPCC 2013), Kunming, China, 5–8 August 2013; pp. 1–5. [Google Scholar]
  32. Li, W.; Yi, W.; Wang, J. A Method for Resolving the Merit Function Expansion of Dynamic Programming TBD. In Proceedings of the 2018 21st International Conference on Information Fusion (FUSION), Cambridge, UK, 10–13 July 2018; pp. 1–7. [Google Scholar]
  33. Yu, H.; Wang, G.; Cao, Q.; Sun, Y. A Fusion Based Particle Filter TBD Algorithm for Dim Targets. Chin. J. Electron. 2015, 24, 590–595. [Google Scholar] [CrossRef]
  34. Chongyang, L.; Jing, J.; Xun, F.; Yangi, W. PF-TBD Algorithm Based on PSO. In Proceedings of the 2016 CIE International Conference on Radar (RADAR), Guangzhou, China, 10–13 October 2016; pp. 1–5. [Google Scholar]
  35. Hongpeng, Z.; Yong, H.; Baichen, J.; Jian, G. A Gaussian Mixture probability Hypothesis Density Smoothing Algorithm for Multi-Target Track-Before-Detect. In Proceedings of the 2016 IEEE International Conference on Signal and Image Processing (ICSIP), Beijing, China, 13–15 August 2016; pp. 667–671. [Google Scholar]
  36. Wang, Z.; Li, M.; Lu, Y.; Chen, H.; Wu, Y. Efficient TR-TBD algorithm for slow-moving weak multi-targets in heavy clutter environment. IET Signal Process. 2017, 11, 422–428. [Google Scholar] [CrossRef]
  37. Skolnik, M. Radar Handbook, 2nd ed.; McGrawHill: New York, NY, USA, 1990. [Google Scholar]
  38. Pakfiliz, A.G.; Efe, M. Multi-Target Tracking in Clutter with Histogram Probabilistic Multi-Hypothesis Tracker. In Proceedings of the 18th International Conference on Systems Engineering (ICSEng’05), Las Vegas, NV, USA, 16–18 August 2005; pp. 137–142. [Google Scholar]
  39. Davey, S.J.; Rutten, M.G.; Cheung, B.; Cheung, B. A Comparison of Detection Performance for Several Track-Before-Detect Algorithms. In Proceedings of the 2008 11th International Conference on Information Fusion, Cologne, Germany, 30 June–3 July 2008; pp. 1–8. [Google Scholar]
  40. Wieneke, M.; Davey, S. Histogram-PMHT for extended targets and target groups in images. IEEE Trans. Aerosp. Electron. Syst. 2014, 50, 2199–2217. [Google Scholar] [CrossRef]
Figure 1. Overall processing flow chart.
Figure 1. Overall processing flow chart.
Sensors 22 07086 g001
Figure 2. Layer integration mode.
Figure 2. Layer integration mode.
Sensors 22 07086 g002
Figure 3. Cases of target association and detection.
Figure 3. Cases of target association and detection.
Sensors 22 07086 g003
Figure 4. The transmission and interaction mechanism based on detection result. (Y is yes, and N is no).
Figure 4. The transmission and interaction mechanism based on detection result. (Y is yes, and N is no).
Sensors 22 07086 g004
Figure 5. The pulse compression result.
Figure 5. The pulse compression result.
Sensors 22 07086 g005
Figure 6. The 0th ambiguity channel and the 1st layer output result. (a) is the result of integration and detection in the 1st segment; (b) is the result of integration and detection in the 2nd segment; (c) is the result of integration and detection in the 3rd segment; (d) is the result of integration and detection in the 4th segment.
Figure 6. The 0th ambiguity channel and the 1st layer output result. (a) is the result of integration and detection in the 1st segment; (b) is the result of integration and detection in the 2nd segment; (c) is the result of integration and detection in the 3rd segment; (d) is the result of integration and detection in the 4th segment.
Sensors 22 07086 g006
Figure 7. The 0th ambiguity channel and the 2nd layer output result. (a) is the result of integration and detection in the 1st segment; (b) is the result of integration and detection in the 2nd segment.
Figure 7. The 0th ambiguity channel and the 2nd layer output result. (a) is the result of integration and detection in the 1st segment; (b) is the result of integration and detection in the 2nd segment.
Sensors 22 07086 g007
Figure 8. The 0th ambiguity channel and the 3rd layer integration and detection result.
Figure 8. The 0th ambiguity channel and the 3rd layer integration and detection result.
Sensors 22 07086 g008
Figure 9. The 1st ambiguity channel and the 1st layer output result. (a) is the result of integration and detection in the 1st segment; (b) is the result of integration and detection in the 2nd segment; (c) is the result of integration and detection in the 3rd segment; (d) is the result of integration and detection in the 4th segment.
Figure 9. The 1st ambiguity channel and the 1st layer output result. (a) is the result of integration and detection in the 1st segment; (b) is the result of integration and detection in the 2nd segment; (c) is the result of integration and detection in the 3rd segment; (d) is the result of integration and detection in the 4th segment.
Sensors 22 07086 g009
Figure 10. The 1st ambiguity channel and the 2nd layer output result. (a) is the result of integration and detection in the 1st segment; (b) is the result of integration and detection in the 2nd segment.
Figure 10. The 1st ambiguity channel and the 2nd layer output result. (a) is the result of integration and detection in the 1st segment; (b) is the result of integration and detection in the 2nd segment.
Sensors 22 07086 g010
Figure 11. The 1st ambiguity channel and the 3rd layer integration and detection result.
Figure 11. The 1st ambiguity channel and the 3rd layer integration and detection result.
Sensors 22 07086 g011
Figure 12. The result of glint target integration; the left figure is proposed method, the right figure is LVD method.
Figure 12. The result of glint target integration; the left figure is proposed method, the right figure is LVD method.
Sensors 22 07086 g012
Figure 13. The output SNRs with different input SNRs.
Figure 13. The output SNRs with different input SNRs.
Sensors 22 07086 g013
Figure 14. Real-measured data for space target.
Figure 14. Real-measured data for space target.
Sensors 22 07086 g014
Figure 15. The integration results of every layer. (af) is the integration result of 1 to 6 layer, respectively.
Figure 15. The integration results of every layer. (af) is the integration result of 1 to 6 layer, respectively.
Sensors 22 07086 g015
Figure 16. Integration SNR of every layer. The red line is the integration SNR of every layer.
Figure 16. Integration SNR of every layer. The red line is the integration SNR of every layer.
Sensors 22 07086 g016
Table 1. Radar parameters in experiments.
Table 1. Radar parameters in experiments.
Radar ParametersValue
Carrier frequency1 GHz
Bandwidth5 MHz
Pulse width 75   μ s
Pulse repetition frequency500 Hz
Table 2. Integration SNR of glint target using LVD and proposed method.
Table 2. Integration SNR of glint target using LVD and proposed method.
MethodSNR InputSNR Output
Ideal−30 dB34.63 dB
Proposed−30 dB33.72 dB
LVD−30 dB22.88 dB
Table 3. The real-measured radar parameters.
Table 3. The real-measured radar parameters.
Radar ParametersValue
Waveband1 GHz
Bandwidth1 MHz
PRF 800 Hz
Integration Time1 s
Table 4. The results of target parameter estimation.
Table 4. The results of target parameter estimation.
RangeVelocityAcceleration
7770 m1249.875 m/s10.5469 m/s2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wei, S.; Dai, Y.; Zhang, Q. Weak and Maneuvering Target Detection with Long Observation Time Based on Segment Fusion for Narrowband Radar. Sensors 2022, 22, 7086. https://doi.org/10.3390/s22187086

AMA Style

Wei S, Dai Y, Zhang Q. Weak and Maneuvering Target Detection with Long Observation Time Based on Segment Fusion for Narrowband Radar. Sensors. 2022; 22(18):7086. https://doi.org/10.3390/s22187086

Chicago/Turabian Style

Wei, Shaopeng, Yan Dai, and Qiang Zhang. 2022. "Weak and Maneuvering Target Detection with Long Observation Time Based on Segment Fusion for Narrowband Radar" Sensors 22, no. 18: 7086. https://doi.org/10.3390/s22187086

APA Style

Wei, S., Dai, Y., & Zhang, Q. (2022). Weak and Maneuvering Target Detection with Long Observation Time Based on Segment Fusion for Narrowband Radar. Sensors, 22(18), 7086. https://doi.org/10.3390/s22187086

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop