Next Article in Journal
Lidar-Derived Rockfall Inventory—An Analysis of the Geomorphic Evolution of Rock Slopes and Modifying the Rockfall Activity Index (RAI)
Next Article in Special Issue
High Precision Motion Compensation THz-ISAR Imaging Algorithm Based on KT and ME-MN
Previous Article in Journal
Combining Chlorophyll Fluorescence and Vegetation Reflectance Indices to Estimate Non-Photochemical Quenching (NPQ) of Rice at the Leaf Scale
Previous Article in Special Issue
Joint Direction of Arrival-Polarization Parameter Tracking Algorithm Based on Multi-Target Multi-Bernoulli Filter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

A Coherent Integration and Parameter Estimation Method for Constant Radial Acceleration Weak Target via SOKT-IAR-LVD

School of Electronic and Optical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(17), 4227; https://doi.org/10.3390/rs15174227
Submission received: 20 July 2023 / Revised: 24 August 2023 / Accepted: 27 August 2023 / Published: 28 August 2023
(This article belongs to the Special Issue Advances in Radar Systems for Target Detection and Tracking)

Abstract

:
In order to enhance the detection and parameter estimation capacity to the maneuvering target with complex motions, a low complexity coherent integration and parameter estimation method named SOKT-IAR-LVD is proposed in this paper. In SOKT-IAR-LVD, first, the second-order keystone transform (SOKT) is utilized to eliminate the range curvature induced by target acceleration. Second, improved axis rotation (IAR) is applied to regulate the linear range migration by rotating the fast time axis and the target envelope is aligned along the slow time axis with a quadratic phase characteristic. At last, the target signal is coherently integrated via the Lv’s Distribution (LVD) transform. The target motion parameters, including range, velocity, and acceleration, are estimated by the IAR and LVD results. The integration gain and computational load of SOKT-IAR-LVD are analyzed. Without needing to estimate the Doppler ambiguity number and target acceleration, the computational burden of SOKT-IAR-LVD is three orders of magnitude lower than that of the Radon-Lv’s Distribution (RLVD) method. Simulation results demonstrate that the detection performance of SOKT-IAR-LVD is almost the same as that of RLVD and that the required input SNR of SOKT-IAR-LVD is 17.4 dB lower than that of SOKT–Radon Fourier transform (SOKT-RFT) when the detection threshold is set to 12 dB.

1. Introduction

With the development of stealth technology, modern ground and aerial targets are able to present the characteristics of high speed, strong maneuverability, long range, and low radar cross-section (RCS) [1]. The traditional moving target detection (MTD) algorithms do not have enough capacity to detect these high-speed maneuvering weak targets. The long-time integration technique is an effective way to increase the signal-to-noise ratio (SNR) and improve radar detection performance [2]. However, when the integration time increases, the high speed and acceleration of the target cause range migration (RM) and Doppler frequency migration (DFM) [3], which limit the performance of classical integration algorithms [4]. Therefore, new approaches to eliminate RM and DFM have been investigated [5].
The typical long-time integration techniques are mainly divided into two categories: incoherent integration and coherent integration. The incoherent integration methods only use the amplitude of echoes to accumulate target signal, which causes poor detection performance in low SNR scenarios [6]. Classical incoherent integration methods include the Hough transform, Radon transform, dynamic programming, and particle filtering methods [7].
Coherent integration performs better than incoherent integration by compensating the phase fluctuation among different sampling pulses. The keystone transform (KT) method corrects the range walk by rescaling the slow time for each range frequency [8]. KT and several improved versions of KT have been widely used, as they can correct the RM effectively without any prior knowledge about the target motion parameter. The Radon Fourier transform (RFT) method eliminates the RM via joint searching along the range and velocity directions of the moving target and integrates coherently via Doppler filtering [9]. The axis rotation MTD (AR-MTD) method eliminates the linear RM by rotating the two-dimensional echoes data plane and realizes coherent integration via the classical MTD algorithm [10]. A fast coherent integration method based on sequence reversing transform was presented in [11], providing a good balance between computational cost and detection ability. By employing the symmetric autocorrelation function and the scaled inverse Fourier transform (SCIFT), a coherent detection algorithm was introduced in [12]; this algorithm can detect high-speed targets without brute-force searching of unknown motion parameters. However, the above methods only consider RM correction, and suffer performance degradation when DFM appears.
KT and second-order KT (SOKT) [13]-based methods have been proposed to regulate both RM and DFM. The KT-matched filtering processing (KT-MFP) method corrects the linear RM via KT and then jointly searches in the fold factor and acceleration domain to remove the residual RM and compensate the DFM. KT-MFP achieves coherent integration through the slow-time Fourier transform (FT) [14]. The SKT-DLVD method uses the segmented keystone transform (SKT) to correct the range walks of targets and then applies the Doppler Lv’s transform (DLVD) to estimate the velocities of targets [15]. The KT-LCT method employs KT to eliminate RM. After this, the linear canonical transform (LCT) is applied to compensate DFM and realize coherent integration [16]. In the low-frequency ultra-wideband synthetic aperture radar (SAR), [17] utilizes the first-order KT to correct the range walk and then uses the SOKT to compensate for the range curvature. However, when Doppler ambiguity occurs due to the high speed of the target or limited pulse repetition frequency (PRF), the Doppler ambiguity number has to be estimated before KT and SOKT processing, which increases the computational burden.
In order to integrate coherently under the condition of Doppler ambiguity, the SOKT-RFT method utilizes SOKT for range curvature correction and the improved de-chirping method for DFM compensation. Then, RFT is applied to correct the linear RM. Because SOKT-RFT eliminates RM and DFM in steps, the integration performance can deteriorate due to compensation errors in the previous steps [18]. The SOKT-MFrRT method uses the SOKT to eliminate quadratic range cell migration and the modified fractional Radon transform (MFrRT) to estimate the ambiguity number of Doppler frequency [19].
To eliminate the RM and DFM effects simultaneously, the Radon-fractional Fourier transform (RFRFT) removes the RM effect via three-dimensional searching within the parameter space and realizes integration using the fractional Fourier transform (FRFT) [20]. Inspired by this, the Radon-Lv’s distribution (RLVD) method eliminates the RM via jointly searching in the target motion parameter space and achieves coherent integration via Lv’s distribution. The RLVD obtains better integration and detection performance than RFRFT [21]. The computational complexity of RLVD is quite large due to the need for multi-dimensional joint searching [22]. In this regard, IAR-FRFT eliminates the linear RM via the improved axis rotation (IAR) transform and realizes coherent integration by FRFT [23]. An approach combining the modified axis rotation transform (MART) and Lv’s transform (LVT) is presented in [24]. Compared with RFRFT and RLVD, the computational cost of IAR-FRFT and MART-LVT is decreased; nevertheless, these two methods suffer integration loss because the range curvature induced by the target acceleration is ignored.
For cases of increasing target maneuverability and long observation time, the short-time generalized radon Fourier transform (STGRFT) method was presented in [25] to detect a maneuvering weak target with multiple motion models. The STGRFT was able to eliminate RM and DFM as well as to estimate the model changing-point time and accumulate the target energy distributed in different motion stages. In the wideband radar scenario, a coherent integration algorithm based on the sub-band keystone transform and extended Lv’s distribution (ELVD) was proposed to estimate the motion parameters and reconstruct the high-resolution range profile (HRRP) of a maneuvering weak target [26]. In [27], high-order motion parameter estimation was modeled as an under-estimated linear regression and the complex-field Bayesian compression sensing (BCS) algorithm was designed to resolve the sparse recovery.
In this paper, a low computational complexity coherent integration algorithm named SOKT-IAR-LVD is proposed to eliminate the RM and DFM caused by the constant radial acceleration of the target. First, the SOKT-IAR-LVD employs SOKT to eliminate the range curvature caused by target acceleration and alleviate linear range migration. Second, the IAR is applied to regulate linear range migration by rotating the fast time axis and the target envelope is aligned along the slow time axis with a quadratic phase characteristic. Third, based on the quadratic phase characteristic of the target signal, the LVD is adopted to accumulate the target signal into a well-focused peak in the centroid frequency–chirp rate (CFCR) domain. The motion parameters of the target are estimated by the rotation angle of IAR and the peak position of LVD outputs. The coherent integration gain and computational complexity of SOKT-IAR-LVD are analyzed.
Without needing to estimate the Doppler ambiguity number and target acceleration, SOKT-IAR-LVD possesses a much lower computational cost than RLVD and achieves better integration performance than the IAR-FRFT, SOKT-RFT, and AR-MTD methods.
The rest of this paper is organized as follows. In Section 2, the echo of a constant radial acceleration target is modeled and the phase terms which cause the RM and DFM are analyzed. In Section 3, the procedure of the proposed SOKT-IAR-LVD method is detailed, and the integration gain and computational complexity are analyzed. Section 4 presents simulation results that demonstrate the efficiency of the SOKT-IAR-LVD method.

2. Signal Model

Assume that the radar transmitted waveform is a linear frequency modulated (LFM) signal, as follows:
s T τ , m = rect τ T P exp j π K τ 2 exp j 2 π f c τ + m T r , τ T P / 2 , m = 0 , 1 , , M 1
where f c is the carrier frequency, T r is the pulse repetition interval (PRI), T P is the pulse width, B is the bandwidth of the LFM waveform, K = B / T P is the frequency modulation rate, τ and m T r denote the fast time and slow time, respectively, M is the number of transmitted LFM pulses in a coherent processing interval (CPI), and
rect u = 1 , u 1 / 2 0 , u > 1 / 2 .
For a target moving towards the radar, the instantaneous range between radar and target is
R τ , m = R 0 v m T r 1 2 a m T r 2 v τ
where R 0 is the initial range between the radar and target and v and a represent the radial velocity and acceleration of the target, respectively. For simplicity, we define
R m = R 0 v m T r 1 2 a m T r 2 .
Because v / c 1 , the baseband target echo can be written as
s R τ , m = rect τ 2 R m / c T P exp j π K τ 2 R m c 2 exp j 4 π f c R m c exp j 2 π f d τ
where c is the speed of light, f d = 2 v / λ is the Doppler frequency, and λ = c / f c is the carrier wavelength. The matched filtered output of the received target echo is
s MF τ , m = sinc B f d τ 2 R m c + f d K exp j π f d τ 2 R m c + f d K × exp j 4 π f c f d R m c exp j π f d 2 K
where sinc x denotes the sinc function. In (5), the peak position of the target envelope in the mth PRI is
τ m = 2 R m c f d K
which varies with the increase of slow time. Denote R OFFSET as the range offset during a CPI, which has
R OFFSET = v M 1 T r + 1 2 a M 1 T r 2 .
When R OFFSET is less than half of the range resolution, that is,
R OFFSET < c 4 B ,
the peak position of the target envelope can be approximately regarded as in the same range cell during the integration time, meaning that the traditional MTD method can be employed for coherent integration. However, the condition of (8) is frequently not met for high-speed maneuvering targets, and the integration performance deteriorates when using the MTD method.
The two range migration terms at the right side of (7) should be eliminated before coherent integration. Applying FT on (5) to the fast time τ , the matched filtered output in the range frequency domain is
S MF f , m = rect f f d / 2 B f d exp j 4 π f + f c f d R 0 c exp j 4 π f + f c f d v c m T r × exp j 2 π f + f c f d a c m T r 2 exp j 2 π f f d K exp j π f d 2 K .
Because of the high target speed or low radar PRF, under-sampling induces ambiguity of the measured target velocity. The velocity of the target is written as
v = N a m b v a m b + v 0
where N a m b is the Doppler ambiguity number, v a m b = λ f P R F / 2 is the blind velocity, f P R F = 1 / T r is the PRF, and v 0 = mod v , v a m b is the measured velocity and satisfies v 0 < v a m b / 2 . Because f d / f c 1 , we can ignore the amount of f d in f + f c f d and write (9) as
S MF f , m = rect f f d / 2 B f d exp j 4 π f + f c R 0 c exp j 4 π f + f c v 0 c m T r × exp j 4 π f N a m b v a m b c m T r exp j 2 π N a m b P R F m T r exp j 2 π f + f c a c m T r 2 × exp j 2 π f f d K exp j π f d 2 K
where exp j 2 π N a m b P R F m T r = 1 . Thus, it has
S MF f , m = rect f f d / 2 B f d exp j 4 π f + f c R 0 c exp j 4 π f + f c v 0 c m T r × exp j 4 π f N a m b v a m b c m T r exp j 2 π f + f c a c m T r 2 exp j 2 π f f d K × exp j π f d 2 K .
There are six exponential terms in (12): φ 1 = exp j 4 π f + f c R 0 / c is the initial range term; φ 2 = exp j 4 π f + f c v 0 m T r / c is the Doppler term determined by the measured velocity v 0 , and results in a linear range migration; φ 3 = exp j 4 π f N a m b v a m b m T r / c is the phase term induced by the blind velocity, and also causes a linear range migration; φ 4 = exp j 2 π f + f c a m T r 2 / c is the frequency modulation term induced by the target acceleration, and leads to the range curvature and DFM; φ 5 = exp j 2 π f d f / K indicates the characteristic of the high-speed motion and causes a fixed target envelope offset from its nominal position; and φ 6 = exp j π f d 2 / K is a constant term.
Here, φ 2 and φ 3 suffer from the first-order coupling between f and m, φ 4 suffers from the second-order coupling between f and m, which causes both range curvature and DFM, and φ 2 , φ 3 , and φ 4 should be decoupled prior to coherent integration.

3. SOKT-IAR-LVD Method

In this section, the SOKT-IAR-LVD method to improve coherent integration performance is detailed. After matched filtering in the range frequency domain, SOKT is employed to eliminate the second-order coupling in φ 4 and correct the range curvature caused by the target acceleration. Then, the IAR is used to eliminate the coupling between range frequency and slow time in the terms φ 2 and φ 3 . At last, the LVD is applied to coherently integrate the target echo. The target motion parameters are estimated by the IAR and LVD results. The details of the SOKT-IAR-LVD method are as follows.

3.1. Range Curvature Correction via SOKT

SOKT is utilized to correct the range curvature caused by the target acceleration. SOKT is a process of rescaling the slow time axis for each range frequency. The scaling formula of SOKT is defined as
m T r = f c f + f c m T r
where m T r denotes the new slow-time variable.
Substituting (13) into (12), the SOKT output in the fast time frequency domain is
S SOKT f , m = S MF f , f c f c + f m = rect f f d / 2 B f d exp j 4 π f + f c R 0 c exp j 4 π f c f c + f v 0 c m T r × exp j 4 π f f c f c + f N a m b v a m b c m T r exp j 2 π f c a c m T r 2 × exp j 2 π f f d K exp j π f d 2 K .
As shown in (14), SOKT removes the second-order coupling in φ 4 for range curvature elimination after matched filtering. As f f c , the following approximations hold:
f c f c + f f c + f / 2 f f c / f c + f f .
Therefore, (14) is rewritten as
S SOKT f , m rect f f d / 2 B f d exp j 4 π f + f c R 0 c exp j 4 π f c v 0 2 + N a m b v a m b m T r × exp j 4 π f c v 0 c m T r exp j 2 π f c a c m T r 2 exp j 2 π f f d K exp j π f d 2 K .
As shown in (16), the second-order coupling between f and m in (12) has been eliminated. After performing the inverse Fourier transform (IFT) on (16) into the τ m domain, we have
s SOKT τ , m = sinc B f d τ 2 c R 0 c f d 2 K v 0 2 + N a m b v a m b m T r × exp j π f d τ 2 c R 0 c f d 2 K v 0 2 + N a m b v a m b m T r exp j 4 π f c R 0 c × exp j 4 π f c v 0 c m T r exp j 2 π f c a c m T r 2 exp j π f d 2 K .
The set τ = 2 r / c (17) can now be rewritten as
s SOKT r , m = sinc 2 B f d c r R 0 c f d 2 K v 0 2 + N a m b v a m b m T r × exp 2 j π f d c r R 0 c f d 2 K v 0 2 + N a m b v a m b m T r exp j 4 π f c R 0 c × exp j 4 π f c v 0 c m T r exp j 2 π f c a c m T r 2 exp j π f d 2 K
where r represents the range corresponding to the fast time τ . In (18), the range offset of the target envelope varies linearly with the slow time m T r . In addition, the linear RM caused by the measured velocity v 0 is reduced to half of its value through SOKT processing.
With the fast time sampling frequency f s , the discrete form of (18) is
s SOKT n r , m = sinc B f d f s n r n 0 + f d f s K + 2 f s c v 0 2 + N a m b v a m b m T r × exp j π f d f s n r n 0 + f d f s K + 2 f s c v 0 2 + N a m b v a m b m T r exp j 4 π f c R 0 c × exp j 4 π f c v 0 c m T r exp j 2 π f c a c m T r 2 exp j π f d 2 K
where n r = round 2 r f s / c , n 0 = round ( 2 R 0 f s / c ) , and n r and n 0 represent the range cell numbers of r and R 0 , respectively.

3.2. Range Migration Correction via IAR

The AR-MTD method concentrates the target echoes in a range cell via the axis rotation transform. However, the Doppler resolution may vary with the axis rotation angle. In the SOKT-IAR-LVD method, the IAR regulates the linear range migration in (19) by rotating the fast time axis. The slow time axis remains unchanged in order to maintain a constant Doppler resolution. Figure 1 shows the diagram of the IAR transform. In (19), the target signal envelope after SOKT is distributed along a straight line with a slope of 2 f s T r v 0 / 2 + N a m b v a m b / c in the coordinate system n r m . In Figure 1, the angle between the target signal envelope and slow time axis m is defined as γ , and has
tan γ = 2 f s T r c v 0 2 + N a m b v a m b .
The IAR processing is
n ˜ = round sin β i m + cos β i n r m ˜ = m
where n ˜ m ˜ represents a new coordinate system after axis rotation, β i = i Δ β π / 2 , π / 2 denotes the axis rotation angle, and Δ β is the angle rotation step. We denote N β as the number of searching angles. Because the target velocity is frequently limited to a certain range, the angle rotation searching area β min , β max can be predetermined in order to reduce the computational burden when the target’s velocity is limited in the region v min , v max with
v min = N a m b , min v a m b + v 0 , min
v max = N a m b , max v a m b + v 0 , max .
According to (20), the angle rotation searching area β min , β max is computed as follows:
β min = atan 2 f s T r v 0 , min / 2 + N a m b , min v a m b / c
β max = atan 2 f s T r v 0 , max / 2 + N a m b , max v a m b / c .
Substituting (20) and (21) into (19) gives
s IAR n ˜ , m ˜ = sinc B f d f s cos β i n ˜ n 0 f d f s K cos β i + B f d f s tan β i tan γ m ˜ × exp j π f d f s cos β i n ˜ n 0 f d f s K cos β i + j π f d f s tan β i tan γ m ˜ × exp j 4 π f c R 0 c exp j 4 π f c v 0 c m ˜ T r exp j 2 π f c a c m ˜ T r 2 exp j π f d 2 K .
When β i = γ , (26) is written as
s IAR n ˜ , m ˜ β i = γ = sinc B f d f s cos γ n ˜ n 0 f d f s K cos γ × exp j π f d f s cos γ n ˜ n 0 f d f s K cos γ exp j 4 π f c R 0 c × exp j 4 π f c v 0 c m ˜ T r exp j 2 π f c a c m ˜ T r 2 exp j π f d 2 K .
As shown in (27), the peak position of the target envelope is concentrated at the range cell n ˜ 0 = round n 0 f d f s / K cos γ . The set n ˜ = n ˜ 0 (27) is simplified as follows:
s IAR m ˜ β i = γ , n ˜ = n ˜ 0 = A 1 exp j 2 π 2 f c v 0 c m ˜ T r exp j π 2 f c a c m ˜ T r 2
with
A 1 = exp j 4 π f c R 0 c exp j π f d 2 K .
In (28), the linear range migration has been corrected by IAR and the output of IAR is a chirp signal with a chirp rate of 2 f c a / c . Inspired by the excellent performance of LVD in extracting the parameters of chirp signals [28], the LVD method is adopted to integrate the target signal in (28) and then estimate the target motion parameters.

3.3. Coherent Integration and Parameter Estimation with LVD

In LVD processing, the symmetric instantaneous autocorrelation function (SIAF) of (28) is
R C m ˜ , l = s IAR m ˜ + l s IAR m ˜ l = A 1 2 exp j 4 π f c v 0 c 2 T r l + j 4 π f c a c 2 T r l m ˜ T r
where ∗ denotes the complex conjugation. Because s IAR m ˜ = 0 when m ˜ < 0 or m ˜ M , the number of valid elements in the SIAF matrix is M 2 / 4 .
The variables m ˜ and l in (30) are coupled with each other in the exponential phase term. We set
m ˜ = m b 2 h T r l
where m b is the LVD slow time variable and h is a scaling factor that determines the chirp rate estimation range of s IAR m ˜ . Substituting (31) into (30), we have
R C m b , l = A 1 2 exp j 4 π f c v 0 c 2 T r l + j 4 π f c a h c T r m b
where the coupling between m ˜ and l has been eliminated. Performing two-dimensional (2D) FT on (32), the output of LVD is
L ( p , q ) = A 2 sin c M T r p 2 M T r 2 f c v 0 c sin c M T r q M T r 2 f c a h c
with
A 2 = M T r 2 exp j 4 π f c R 0 c 2 exp j π f d 2 K 2
where p q represents the centroid frequency–chirp rate (CFCR) domain with M / 2 p , q M / 2 1 . The target echo is coherently accumulated at p = 4 M T r v 0 / λ q = 2 M T r a / λ h . Therefore, the velocity and acceleration of target are estimated by the peak position of the LVD results:
p ^ , q ^ = arg max p , q L ( p , q )
v ^ 0 = λ p ^ 4 M T r
a ^ = λ h q ^ 2 M T r .
According to (37), the scaling factor h should be no less than 4 T r a m a x / λ , where a m a x denotes the maximal target radial acceleration.
Finally, based on the IAR and LVD results, the Doppler ambiguity number N a m b and target velocity v are estimated as follows:
N ^ a m b = round c tan γ + v ^ 0 T r f s 2 T r f s v a m b
v ^ = N ^ a m b v a m b + v ^ 0 .

3.4. Procedure of the SOKT-IAR-LVD Method

Based on the above analysis, the flow chart of the SOKT-IAR-LVD method is shown in Figure 2. The procedure of SOKT-IAR-LVD is divided into the following steps.
Step 1. Initialize the parameters of the SOKT-IAR-LVD method, including the rotation angle searching area β min , β max , angle search step Δ β , and scaling factor h.
Step 2. Perform matched filtering on the radar echoes in the frequency domain.
Step 3. Apply SOKT to the matched filtered outputs to eliminate the range curvature, then perform range IFT on the SOKT outputs.
Step 4. For each β i β min , β max , apply IAR to the SOKT outputs to remove the residual linear RM, then perform LVD transform on the IAR results along the slow time.
Step 5. After computing the LVD results for all rotation angles in β min , β max , search for the peak of the LVD results at each range cell and the corresponding values of β i , p ^ , and q ^ .
Step 6. Apply constant false alarm rate (CFAR) detection to the peak of the LVD results. If a target is detected, the range, velocity, and acceleration of the target are estimated using the values of β i , p ^ , and q ^ using (35)–(39), respectively.

3.5. Integration Gain Analysis

The coherent integration gain of SOKT-IAR-LVD mainly depends on the LVD transform. We denote the output SNR of matched filtering as SNR MF , and the output SNR of the IAR transform is the same as SNR MF . Adding the complex Gaussian noise to the IAR output in (28) with SNR MF = A s 2 A s 2 ζ 2 ζ 2 , we have
y m ˜ = s IAR m ˜ + n m ˜
where n m ˜ represents the noise with power ζ 2 and the amplitude s IAR m ˜ is A s . The SIAF of y m ˜ is
R y C m ˜ , l = y m ˜ + l y m ˜ l = s IAR m ˜ + l s IAR m ˜ l + s IAR m ˜ + l n m ˜ l + n m ˜ + l s IAR m ˜ l + n m ˜ + l n m ˜ l .
In (41), the amplitude of the target signal changes to A s 2 , while the noise power is 2 A s 2 ζ 2 + ζ 4 . Because the number of valid elements in the SIAF matrix is M 2 / 4 , the output SNR of LVD is
SNR LVD = M 2 A s 4 4 ( 2 A s 2 ζ 2 + ζ 4 ) .
Therefore, the integration gain of LVD is
G LVD = SNR LVD SNR MF = 10 lo g 10 M 2 / M 2 8 + 4 SNR MF 8 + 4 SNR MF dB .
The integration gain of the RFT- and FRFT-based methods only depends on the pulse number M, which is
G RFT _ FRFT = 10 lo g 10 M dB .
when the output SNR of matched filtering satisfies
SNR MF > 4 4 M 8 M 8 .
The SOKT-IAR-LVD method achieves a higher integration gain than the RFT and FRFT based methods. Condition (45) frequently holds for long coherent integration scenarios.

3.6. Computation Complexity Analysis

Next, the computational load of SOKT-IAR-LVD is analyzed. Because SOKT and LVD can be realized by Chirp-z transform–inverse fast FT (CZT-IFFT) [29] and scaled FT–IFFT (SFT-IFFT) [30], respectively, the computational complexities of SOKT and IAR are O M r M log 2 M and O M r N β , where M r = f s T r is the number of range cells in a PRI. The computation complexity of LVD at a rotation angle β i is O M 2 log 2 M . Therefore, the overall computational cost of SOKT-IAR-LVD is approximately O N β M r M 2 log 2 M .

4. Simulation Results

In this section, the performance of SOKT-IAR-LVT is evaluated and compared with that of the RLVD, SOKT-RFT, IAR-FRFT, and AR-MTD methods. Table 1 shows the parameters of the transmitted radar waveform. Table 2 presents the parameters of the SOKT-IAR-LVD method. As shown in Table 2, the searching rotation angles in SOKT-IAR-LVD is N β = 6573 and the searching number of acceleration is 512. Accordingly, the searching numbers of velocity and acceleration in RLVD, SOKT-RFT, IAR-FRFT, and AR-MTD are set as 6573 and 512, respectively.

4.1. Integration Performance with Different Input SNRs

First, the performance of the SOKT-IAR-LVD method is evaluated in a high SNR condition. Consider a constant radial acceleration target moving towards the radar with initial distance R 0 = 10 km , radial velocity v = 1000 m / s , and acceleration a = 100 m / s 2 . The theoretical numbers of cross-range cells caused by the target velocity and acceleration are 102 and 1, respectively. The theoretical slope of the target envelope after SOKT is −0.19, corresponding to a rotation angle of γ = 10 . 758 .
When the SNR of the received target echo in the time domain is −10 dB, Figure 3a shows that the matched filtered output results in a CPI. The range cell migration in Figure 3a is 103, which is consistent with the theoretical value. Figure 3b shows the SOKT outputs for range curvature correction. The theoretical number of cross-range cells removed by SOKT is f s M 1 T r v 0 + a M 1 T r / c 6 . As shown in Figure 3b, the range cell migration reduces to 97.
When β i = 1157 Δ β = 10.759 , the results of applying the IAR transform to the SOKT results in Figure 3b are shown in Figure 3c. In Figure 3c, the linear range migration has been corrected. The target envelope in every PRI is aligned in the 3930th range cell. The LVD output after applying the LVD transform to the target range cell in Figure 3c is presented in Figure 3d. In Figure 3d, the target signal is coherently integrated at p ^ = 127 , q ^ = 57 . The LVD output SNR is 56.28 dB. Therefore, the overall integration gain, including matched filtering and LVD, is 66.28 dB. The theoretical processing gains in matched filtering and LVD are 21.76 dB and 45.01 dB, respectively. The simulated integration gain accurately coincides with the theoretical value.
The peak value of LVD output at each β i when the rotation angle β i varies from 46 to 0 is shown in Figure 3e. In Figure 3e, the LVD output achieves a maximum at β i = 10 . 759 . Using (35)–(39), the estimated range, velocity, and acceleration of the target are 10.001 km , 999.99 m / s , and 100.19 m / s 2 , respectively. SOKT-IAR-LVD achieves high target motion parameter estimation accuracy.
Next, the performance of SOKT-IAR-LVD in the low SNR condition is evaluated. For a constant radial acceleration weak target moving towards the radar with R 0 = 15 km , v = 4500 m / s and a = 350 m / s 2 , the matched filtered outputs are shown in Figure 4a when the input SNR = 35 dB. In Figure 4a, the target signal is submerged by noise and the matched filtered output SNR is very low. After SOKT and IAR transform at β i = 41 . 986 , the LVD results are shown in Figure 4b. The weak target signal is integrated with the output SNR = 21.10 dB, which is conducive to target detection. The total integration gain is 56.10 dB, which is consistent with the theoretical value of 56.29 dB. The peak position of LVD outputs is at p ^ = 153 , q ^ = 199 . The estimated range, velocity, and acceleration of the target are 15.001 km , 4500.02 m / s , and 349.81 m / s 2 , respectively.
Next, the coherent integration gain of SOKT-IAR-LVD at different input SNRs is evaluated. The target motion parameter is the same as in Figure 3. Figure 5 presents the output SNR curves of five methods when the input SNR increases from −45 dB to −20 dB. In Figure 5, when the output SNR after integration is set to 12 dB the required input SNRs are −39.8 dB, −39.6 dB, −36.3 dB, −22.4 dB, and −20.0 dB for the SOKT-IAR-LVD, RLVD, IAR-FRFT, SOKT-RFT, and AR-MTD methods, respectively. The corresponding input SNR of SOKT-IAR-LVD is 2.9 dB, 17.4 dB, and 19.8 dB lower than those of IAR-FRFT, SOKT-RFT, and AR-MTD, respectively. SOKT-IAR-LVD and RLVD achieve a similar processing gain and possess superior integration ability to SOKT-RFT, IAR-FRFT, and AR-MTD. Because SOKT-RFT estimates the target acceleration by de-chirping and averaging, it inevitably suffers from estimation errors. Therefore, the integration gain of SOKT-RFT is lower than SOKT-IAR-LVD and RLVD. As the target acceleration is not compensated in IAR-FRFT and AR-MTD, the integration performance of these methods is inferior to that of SOKT-IAR-LVD.

4.2. Integration Performance with Different Motion Parameters

In this subsection, the coherent integration performance for different target motion parameters is evaluated. The input SNR is set as −30 dB. With a constant radial acceleration of a = 100 m / s 2 , Figure 6a shows the output SNR when the target velocity increases from 500 m / s to 5000 m / s . Figure 6b presents the output SNR when the target acceleration varies from 50 m / s 2 to 450 m / s 2 and v = 1000 m / s . In Figure 6a,b, SOKT-IAR-LVD and RLVD achieve the highest integration performance. The integration gain of SOKT-IAR-LVD fluctuates slightly around 30.43 dB during the variation of target velocity and acceleration, and is consistent with the theoretical value of 30.55 dB. In Figure 6b, the performance of IAR-FRFT and AR-MTD degrades when the target acceleration increases, which is because the acceleration is not considered in these methods. Suffering from the acceleration estimation error, the performance of SOKT-RFT is stable but inferior to that of SOKT-IAR-LVD.

4.3. Detection Performance

The detection performance of the five methods is evaluated using the cell-averaging CFAR (CA-CFAR) algorithm. The target motion parameter is the same as in Figure 3. The nominal false alarm rate is set to P f a = 10 6 . The detection probability curves when the input SNR varies from −45 dB to −20 dB are plotted in Figure 7 via 1000 Monte Carlo trials for each SNR. As shown in Figure 7, SOKT-IAR-LVD and RLVD achieve almost the same detection probability and are better than IAR-FRFT, SOKT-RFT, and AR-MTD. At P d = 80 % , the input SNR of SOKT-IAR-LVD is about 6.5 dB lower than that of IAR-FRFT and at least 16.5 dB lower than that of SOKT-RFT. The detection ability of SOKT-IAR-LVD for the weak target is superior to IAR-FRFT, SOKT-RFT, and AR-MTD.

4.4. Computational Complexity Comparison

In this subsection, the computational complexity of the five methods is evaluated. The parameters of radar transmit waveform and SOKT-IAR-LVD are the same as in Table 1 and Table 2. As shown in Table 2, the search numbers of rotation angle and acceleration in SOKT-IAR-LVD are selected as N β = 6573 and M = 512 , respectively. Accordingly, the search numbers of velocity and acceleration in RLVD, IAR-FRFT, SOKT-RFT, and AR-MTD are set as N v = 6573 and N a = 512 , respectively. Using the derivation results from Section 3.6, the computation cost of SOKT-IAR-LVD for one range cell is O N β M 2 log 2 M , while that of RLVD is O N v N a M 2 log 2 M [24]. Because N β = N v = 6573 and M = N a = 512 , the computational cost of SOKT-IAR-LVD is nearly three orders of magnitude lower than that of RLVD in the theoretical analysis.
Table 3 shows the number of complex multiplications and additions during one range cell coherent integration. The runtime of one hundred range cells integration on a PC equipped with an Intel Core i5-8250U (1.6 GHz) and 8 GB RAM is shown in Table 4.
From Table 3 and Table 4, it can be seen the computational complexity of SOKT-IAR-LVD is similar to that of IAR-FRFT and nearly three orders of magnitude lower than that of RLVD. Although the computational burdens of SOKT-RFT and AR-MTD are much lower than that of SOKT-IAR-LVD, the detection performance of these methods for the constantly accelerating targets is poor. SOKT-IAR-LVD achieves the best coherent integration performance with a moderate computation cost while estimating the target motion parameter with high accuracy.

5. Conclusions

In order to improve the detection performance for the maneuvering targets through long-time integration, the SOKT-IAR-LVD method utilizes the SOKT to eliminate the range curvature induced by target acceleration. Then, it adopts the IAR to remove the linear range migration. At last, the target signal is coherently integrated via the LVD transform. Because estimation of the Doppler ambiguity number and target acceleration is not required, the computational burden of SOKT-IAR-LVD is three orders of magnitude lower than that of RLVD. The detection performance of SOKT-IAR-LVD is almost the same as that of RLVD, and much better than the IAR-FRFT, SOKT-RFT, and AR-MTD methods in the low SNR scenarios. When the detection threshold is set to 12 dB, the required input SNR of SOKT-IAR-LVD is 2.9 dB, 17.4 dB, and 19.8 dB lower than those of IAR-FRFT, SOKT-RFT, and AR-MTD, respectively. Meanwhile, SOKT-IAR-LVD can estimate the target motion parameter accurately.

Author Contributions

Writing—original draft preparation, R.Z.; investigation, N.X.; writing—review and editing, N.X. and R.Z.; project administration, R.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China under Grants 61971224, 62001227, and 62001232.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank the reviewers for their great help on the article during its review progress.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Huang, X.; Tang, S.; Zhang, L.; Li, S. Ground-Based Radar Detection for High-Speed Maneuvering Target via Fast Discrete Chirp-Fourier Transform. IEEE Access 2019, 7, 12097–12113. [Google Scholar] [CrossRef]
  2. Li, X.; Kong, L.; Cui, G.; Yi, W. Detection and RM correction approach for manoeuvring target with complex motions. IET Signal Process. 2017, 11, 378–386. [Google Scholar] [CrossRef]
  3. Fang, X.; Min, R.; Cao, Z.; Pi, Y. High-order RM and DFM correction method for long-time coherent integration of highly maneuvering target. Signal Process. 2019, 162, 221–233. [Google Scholar] [CrossRef]
  4. Hou, L.; Song, H.; Zheng, M.; Wang, W. Fast moving target imaging and motion parameters estimation based on Radon transform and bi-directional approach. IET Radar Sonar Navig. 2016, 10, 1013–1023. [Google Scholar] [CrossRef]
  5. He, X.; Liao, G.; Zhu, S.; Xu, J.; Guo, Y.; Wei, J. Fast non-searching method for ground moving target refocusing and motion parameters estimation. Digit. Signal Process. 2018, 79, 152–163. [Google Scholar] [CrossRef]
  6. Huang, P.; Xia, X.G.; Liao, G.; Yang, Z.; Zhang, Y. Long-Time Coherent Integration Algorithm for Radar Maneuvering Weak Target with Acceleration Rate. IEEE Trans. Geosci. Remote Sens. 2019, 57, 3528–3542. [Google Scholar] [CrossRef]
  7. Yang, J.; Liu, C.; Wang, Y. Detection and Imaging of Ground Moving Targets with Real SAR Data. IEEE Trans. Geosci. Remote Sens. 2015, 53, 920–932. [Google Scholar] [CrossRef]
  8. Perry, R.; DiPietro, R.; Fante, R. SAR imaging of moving targets. IEEE Trans. Aerosp. Electron. Syst. 1999, 35, 188–200. [Google Scholar] [CrossRef]
  9. Xu, J.; Yu, J.; Peng, Y.N.; Xia, X.G. Radon-Fourier Transform for Radar Target Detection, I: Generalized Doppler Filter Bank. IEEE Trans. Aerosp. Electron. Syst. 2011, 47, 1186–1202. [Google Scholar] [CrossRef]
  10. Rao, X.; Tao, H.; Su, J.; Guo, X.; Zhang, J. Axis rotation MTD algorithm for weak target detection. Digit. Signal Process. 2014, 26, 81–86. [Google Scholar] [CrossRef]
  11. Li, X.; Cui, G.; Yi, W.; Kong, L. Sequence-Reversing Transform-Based Coherent Integration for High-Speed Target Detection. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 1573–1580. [Google Scholar] [CrossRef]
  12. Zheng, J.; Su, T.; Zhu, W.; He, X.; Liu, Q.H. Radar High-Speed Target Detection Based on the Scaled Inverse Fourier Transform. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 1108–1119. [Google Scholar] [CrossRef]
  13. Yang, J.; Zhang, Y. An Airborne SAR Moving Target Imaging and Motion Parameters Estimation Algorithm with Azimuth-Dechirping and the Second-Order Keystone Transform Applied. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 3967–3976. [Google Scholar] [CrossRef]
  14. Sun, Z.; Li, X.; Yi, W.; Cui, G.; Kong, L. Detection of weak maneuvering target based on keystone transform and matched filtering process. Signal Process. 2017, 140, 127–138. [Google Scholar] [CrossRef]
  15. Tian, J.; Cui, W.; Xia, X.G.; Wu, S.L. Parameter Estimation of Ground Moving Targets Based on SKT-DLVT Processing. IEEE Trans. Comput. Imaging 2016, 2, 13–26. [Google Scholar] [CrossRef]
  16. Huang, X.; Zhang, L.; Li, S.; Zhao, Y. Radar high speed small target detection based on keystone transform and linear canonical transform. Digit. Signal Process. 2018, 82, 203–215. [Google Scholar] [CrossRef]
  17. Jungang, Y.; Xiaotao, H.; Thompson, J.; Jin, T.; Zhimin, Z. Low-frequency ultra-wideband synthetic aperture radar ground moving target imaging. IET Radar Sonar Navig. 2011, 5, 994–1001. [Google Scholar] [CrossRef]
  18. Tian, J.; Cui, W.; Shen, Q.; Wei, Z.X.; Wu, S.L. High-speed maneuvering target detection approach based on joint RFT and keystone transform. Sci. China Inf. Sci. 2013, 56, 1–13. [Google Scholar] [CrossRef]
  19. Wang, Y.; Zhou, F.; Wu, Y.; Bao, Z. Improved ambiguity estimation using a modified fractional Radon transform. IET Radar Sonar Navig. 2011, 5, 489–495. [Google Scholar]
  20. 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]
  21. 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]
  22. Wang, C.; Jiu, B.; Liu, H. Maneuvering target detection in random pulse repetition interval radar via resampling-keystone transform. Signal Process. 2021, 181, 107899. [Google Scholar] [CrossRef]
  23. 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]
  24. Li, X.; Sun, Z.; Yi, W.; Cui, G.; Kong, L. Radar Detection and Parameter Estimation of High-Speed Target Based on MART-LVT. IEEE Sensors J. 2019, 19, 1478–1486. [Google Scholar] [CrossRef]
  25. Li, X.; Sun, Z.; Yeo, T.S.; Zhang, T.; Yi, W.; Cui, G.; Kong, L. STGRFT for Detection of Maneuvering Weak Target with Multiple Motion Models. IEEE Trans. Signal Process. 2019, 67, 1902–1917. [Google Scholar] [CrossRef]
  26. Zhang, Q.; Zhang, Y.; Sun, B.; Dong, Z.; Yu, L. Motion Parameters Estimation and HRRP Reconstruction of Maneuvering Weak Target for Wideband Radar Based on SKT-ELVD. IEEE Trans. Aerosp. Electron. Syst. 2023, 59, 2752–2763. [Google Scholar] [CrossRef]
  27. Yu, L.; He, F.; Zhang, Q.; Su, Y.; Zhao, Y. High Maneuvering Target Long-time Coherent Integration and Motion Parameters Estimation Based on Bayesian Compressive Sensing. IEEE Trans. Aerosp. Electron. Syst. 2023, 1–16. [Google Scholar] [CrossRef]
  28. 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]
  29. 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]
  30. Luo, S.; Bi, G.; Lv, X.; Hu, F. Performance analysis on Lv distribution and its applications. Digit. Signal Process. 2013, 23, 797–807. [Google Scholar] [CrossRef]
Figure 1. IAR transform for linear range migration elimination.
Figure 1. IAR transform for linear range migration elimination.
Remotesensing 15 04227 g001
Figure 2. Flowchart of the SOKT-IAR-LVD method.
Figure 2. Flowchart of the SOKT-IAR-LVD method.
Remotesensing 15 04227 g002
Figure 3. Target integration via SOKT-IAR-LVD method: (a) matched filtered outputs, (b) results of SOKT, (c) results of IAR transform, (d) LVD outputs in the target range cell, (e) outputs of LVD peaks at different rotation angles.
Figure 3. Target integration via SOKT-IAR-LVD method: (a) matched filtered outputs, (b) results of SOKT, (c) results of IAR transform, (d) LVD outputs in the target range cell, (e) outputs of LVD peaks at different rotation angles.
Remotesensing 15 04227 g003
Figure 4. Coherent integration of SOKT-IAR-LVD for a weak target: (a) matched filtered output and (b) LVD output in the target range cell.
Figure 4. Coherent integration of SOKT-IAR-LVD for a weak target: (a) matched filtered output and (b) LVD output in the target range cell.
Remotesensing 15 04227 g004
Figure 5. Coherent integration performance of various method when input SNR increases from −45 dB to −20 dB.
Figure 5. Coherent integration performance of various method when input SNR increases from −45 dB to −20 dB.
Remotesensing 15 04227 g005
Figure 6. Output SNR results versus target’s motion parameters: (a) output SNR versus the target’s velocity ( a = 100 m / s 2 ) and (b) output SNR versus the target’s acceleration ( v = 1000 m / s ).
Figure 6. Output SNR results versus target’s motion parameters: (a) output SNR versus the target’s velocity ( a = 100 m / s 2 ) and (b) output SNR versus the target’s acceleration ( v = 1000 m / s ).
Remotesensing 15 04227 g006
Figure 7. Detection probability of the five methods when input SNR increases from −45 dB to −20 dB.
Figure 7. Detection probability of the five methods when input SNR increases from −45 dB to −20 dB.
Remotesensing 15 04227 g007
Table 1. Parameters of transmitted radar waveform.
Table 1. Parameters of transmitted radar waveform.
ParameterValue
Carrier frequency, f c 1 GHz
Bandwidth, B15 MHz
Pulse width, T P 10 μ s
Pulse repetition frequency, f P R F 2000 Hz
Sampling rate, f s 60 MSPS
Number of integrated pulses, M512
Table 2. Parameters of SOKT-IAR-LVD.
Table 2. Parameters of SOKT-IAR-LVD.
ParameterValue
Velocity measured region 0 m / s , 5000 m / s
Acceleration measured region 0 m / s 2 , 450 m / s 2
Rotation angle searching area β min , β max 46 , 0
Rotation angle search step, Δ β 0 . 007
Scaling factor, h3
Table 3. Computational costs of the five methods.
Table 3. Computational costs of the five methods.
MethodComplex MultiplicationsComplex Additions
SOKT-IAR-LVD 3.7908 × 10 10 6.2031 × 10 10
RLVD 1.9409 × 10 13 3.1760 × 10 13
IAR-FRFT 5.5138 × 10 10 9.3046 × 10 10
SOKT-RFT 6.0027 × 10 6 8.0958 × 10 6
AR-MTD 1.5144 × 10 7 3.0288 × 10 7
Table 4. Runtimes of the five methods.
Table 4. Runtimes of the five methods.
MethodRuntime (s)
SOKT-IAR-LVD 15.5545
RLVD 6821.4771
IAR-FRFT 19.0411
SOKT-RFT 0.0098
AR-MTD 0.0167
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, R.; Xu, N. A Coherent Integration and Parameter Estimation Method for Constant Radial Acceleration Weak Target via SOKT-IAR-LVD. Remote Sens. 2023, 15, 4227. https://doi.org/10.3390/rs15174227

AMA Style

Zhang R, Xu N. A Coherent Integration and Parameter Estimation Method for Constant Radial Acceleration Weak Target via SOKT-IAR-LVD. Remote Sensing. 2023; 15(17):4227. https://doi.org/10.3390/rs15174227

Chicago/Turabian Style

Zhang, Renli, and Nan Xu. 2023. "A Coherent Integration and Parameter Estimation Method for Constant Radial Acceleration Weak Target via SOKT-IAR-LVD" Remote Sensing 15, no. 17: 4227. https://doi.org/10.3390/rs15174227

APA Style

Zhang, R., & Xu, N. (2023). A Coherent Integration and Parameter Estimation Method for Constant Radial Acceleration Weak Target via SOKT-IAR-LVD. Remote Sensing, 15(17), 4227. https://doi.org/10.3390/rs15174227

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