3 August 2016 Radar high-speed maneuvering target detection based on joint second-order keystone transform and modified integrated cubic phase function
Author Affiliations +
J. of Applied Remote Sensing, 10(3), 035009 (2016). doi:10.1117/1.JRS.10.035009
Abstract
Long-time coherent integration is an effective means to improve the radar detection ability of high-speed maneuvering targets with jerk motion. However, the range migration (RM) and Doppler frequency migration (DFM) have a great impact on the integration performance. To overcome these problems, a unique method, called the second-order keystone transform modified integrated cubic phase function (SKT-MICPF), is proposed. In this method, the velocity compensation and SKT are jointly employed to correct the RM. After the RM correction, the azimuth echoes of a range cell where a target is located can be modeled as a cubic phase signal (CPS), whose chirp rate (CR) and quadratic CR are related to the target’s radial acceleration and jerk, respectively. Thereafter, an effective parameters’ estimation algorithm for CPS, called MICPF, is proposed and applied to compensate the DFM. After that, coherent integration and target detection are accomplished via the fast Fourier transform and constant false alarm rate technique, successively. Compared with the improved axis rotation discrete chirp Fourier transform, the SKT-MICPF achieves close detection performance, but greatly reduces the computational complexity. The results of simulation and real radar data demonstrate the validity of the proposed algorithm.
Zhang, Su, Li, and Zheng: Radar high-speed maneuvering target detection based on joint second-order keystone transform and modified integrated cubic phase function

1.

Introduction

Modern national defenses have higher requirements for the detection of high-speed maneuvering targets.1,2 Prolonging the integration time is an effective means to improve the detection performance. However, with the increasing of the integration time, the integration performance will degrade due to the range migration (RM) and Doppler frequency migration (DFM). To address these problems, various successful algorithms12.3.4.5.6.7.8.9.10.11.12.13.14.15 have been developed.

The Radon Fourier transform (RFT),2 the keystone transform (KT),3 the axis rotation moving target detection,4 and the scaled inverse Fourier transform-based method5 can be applied to correct the RM and accomplish coherent integration for a target with a uniform radial velocity. The Radon fractional FT,6 the Radon Lv’s distribution (RLVD),7 the KT Lv’s transform,8 the two-dimensional (2-D) match filter-based method,9 and the KT-RFT10 are the effective methods for a high-speed maneuvering target with a uniform radial acceleration. However, for a highly maneuvering target with a jerky motion, the aforementioned algorithms23.4.5.6.7.8.9.10 will suffer from integration loss due to the RM and DFM induced by the target’s jerky motion.

For a high-speed maneuvering target with a uniform jerk, an algorithm based on the generalized KT and second-order dechirp process is proposed in Ref. 11, which employs the range frequency domain form of the compressed echoes, i.e., Eq. (5), to accomplish the target’s motion parameters estimation. However, because the target’s energy spreads in the range frequency domain, this method suffers from great signal-to-noise ratio (SNR) loss. The cubic phase function (CPF)-based methods12,13 and the adjacent correlation function LVT14 also face similar problems. An improved axis rotation discrete chirp Fourier transform (IAR-DCFT) method is presented in Ref. 15, which can correct the RM by the IAR, and achieve coherent integration via the DCFT. Nevertheless, the implementation of the DCFT is computationally demanding.16,17

After the RM correction, the azimuth echoes of a range cell where a target is located can be considered as a cubic phase signal (CPS), whose chirp rate (CR) and quadratic CR (QCR) are related to the target’s radial acceleration and jerk, respectively. Therefore, the estimation algorithms for the CPS can be applied to estimate the target’s radial acceleration and jerk to further compensate the DFM. The scaled FT (SCFT)-based method,18 the keystone time chirp rate distribution,19 and the modified chirp rate quadratic chirp rate distribution20 are the successful algorithms for CPS. However, due to the fourth-order kernel functions that are applied, the antinoise performances of those methods are still low. The integrated cubic phase function (ICPF)21 which is developed from the CPF22,23 can obtain a higher antinoise performance than the aforementioned approaches1819.20 due to its bilinearity. In Ref. 24, a space switching-based ICPF (SS-ICPF) algorithm is proposed for CPS. Nevertheless, the space switch method which includes the circular shifting operation and map operation would burden the complexity of the estimator and increase the computational cost.

Motivated by the previous work, an algorithm, known as second-order keystone transform modified integrated cubic phase function (SKT-MICPF), is presented for the high-speed maneuvering targets with jerky motions. First, the velocity compensation and SKT are jointly employed to correct the RM. Second, an estimation method for CPS, named MICPF, is proposed and applied to compensate the DFM. Finally, the fast Fourier transform (FFT) and constant false alarm rate (CFAR) technique are employed to accomplish coherent integration and target detection, successively. Comparisons with the IAR-DCFT show that the SKT-MICPF can reduce the computational cost greatly with little detection performance loss.

The rest of this paper is organized as follows. In Sec. 2, the signal model for high-speed maneuvering targets with jerky motion is built. The SKT-based RM correction method is presented in Sec. 3. In Sec. 4, an MICPF method is proposed for CPS and its performances are also analyzed. A detection approach for high-speed maneuvering targets with jerky motion is illustrated in Sec. 5. In Sec. 6, the experiments with synthetic data and real data are performed to demonstrate the effectiveness of the proposed algorithm. Section 7 gives the conclusions.

2.

Signal Model for High-Speed Maneuvering Targets

Suppose that the radar transmits a linear frequency modulation signal, which can be expressed as

(1)

st(t)=rect(t/Tr)exp[jπ(2fct+γt2)],
where rect(tTr)={1,|t|Tr/20,|t|>Tr/2, Tr is the pulse duration, and t, fc, and γ denote the fast time, the carrier frequency, and the modulation rate, respectively.

Assume that there are P high-speed maneuvering targets with a constant radial jerk in the scene, thus the instantaneous slant range between the radar and the p’th target satisfies

(2)

Rp(tm)=rp+vTptm+aptm2/2+kptm3/6,
where tm represents the slow time, and rp,vp,ap, and kp denote the initial slant range, the radial velocity, acceleration, and jerk of the p’th target, respectively.

The received baseband echoes can be expressed as2

(3)

Sr(tm,t)=p=1PAprect(t2Rp(tm)/cTr)exp{jπγ[t2Rp(tm)c]2}×exp[j2πfc2Rp(tm)c]+N(tm,t),tm[T/2,T/2],
where Ap denotes the p’th target reflectivity, and c and T represent the speed of light and the integration time, respectively. N(tm,t) is the complex white Gaussian noise.

After the pulse compression, the echoes in the slow time-fast time domain can be stated as

(4)

Sc(tm,t)=p=1PAcpsinc{B[t2(rp+vTptm+aptm2/2+kptm3/6)c]}×exp[j2π2(rp+vTptm+aptm2/2+kptm3/6)λ]+Nc(tm,t),
where B denotes the bandwidth, Nc(tm,t) is the compressed noise, and λ=c/fc is the wavelength.

3.

Range Migration Correction

Performing the range FFT on Eq. (4), we get the compressed signal in the slow time range frequency domain.

(5)

Sc(tm,fr)=p=1PA1prect(frB)exp[j2π(fc+fr)2(rp+vTptm+aptm2/2+kptm3/6)c],
where fr denotes the range frequency.

In order to correct the RMs induced by the targets’ radial velocities, a compensation function is constructed, which takes the form

(6)

H1(tm,fr,vs)=exp[j2π(fc+fr)2vstmc],
where vs is the searching velocity.

Multiplying Eq. (6) with Eq. (5) yields

(7)

S1(tm,fr)=H1(tm,fr,vs)×Sc(tm,fr)=p=1PA1prect(frB)exp{j2π(fc+fr)2[rp+(vTpvs)tm+aptm2/2+kptm3/6]c}.
To correct the RMs caused by the targets’ radial accelerations, we perform the SKT on Eq. (7). The scaling equation of the SKT takes the form10

(8)

tm=(fcfc+fr)1/2tn.
After the SKT, we get

(9)

S2(tn,fr)=p=1PA1prect(frB)exp[j2π(fc+fr)2rpc]exp[j2πfc(1+frfc)1/22(vTpvs)tnc]×exp(j2πfcaptn2c)exp[j2πfc(1+frfc)1/2kptn3/3c].
The SKT can be implemented via the FFT-based chirp-z transform with the computational cost of O(MNlog2M),3 where M and N denote the number of integration pulses and range cells, respectively. For a narrowband radar, fcfr is true.11 By using the first-order approximations (1+fr/fc)1/21+fr/(2fc) and (1+fr/fc)1/21fr/(2fc), Eq. (9) can be rewritten as

(10)

S2(tn,fr)p=1PA1prect(frB)exp[j2π(fc+fr)2rpc]exp[j2π(fc+12fr)2(vTpvs)tnc]×exp(j2πfcaptn2c)exp[j2π(fc12fr)kptn3/3c].
Applying the inverse FFT (IFFT) on Eq. (10) along the fr-axis, we get

(11)

S3(tn,t)p=1PA2psinc{B[t2[rp+(vTpvs)tn/2kptn3/12]c]}×exp{j2π2[rp+(vTpvs)tn+aptn2/2+kptn3/6]λ}.
From Eq. (11), it can be easily seen that the RMs induced by the targets’ radial accelerations and part of the RMs caused by the targets’ radial jerks have been corrected.

Assume that the radial velocity of the p’th target is matched by the searching velocity, i.e., |vTpvs|Δvs/2, where Δvs denotes the interval of the searching velocity. Therefore, the residual RM of the p’th target satisfies

(12)

(|(vTpvs)2tn||kp|12tn3)|tn=±T/2Δvs4T+|kp|max48T3.
The RMs induced by the targets’ radial jerks are usually small for a narrowband radar.12,15 We limit the residual RMs induced by the targets’ radial jerks to be half a range resolution in this paper. If the interval of searching velocity satisfies

(13)

Δvs2ΔR/T,
where ΔR=c/(2B) is the range resolution, the residual RM of the p’th target will be less than a range resolution, which indicates that the RM effect can be ignored.2 Thus, the interval of searching velocity can be selected in terms of Eq. (13).

With the searching procedure, all the targets can be corrected successively. However, the DFMs induced by the targets’ radial accelerations and jerks still remain, which have a great impact on the integration performance.

4.

Parameters Estimation of CPS Based on Modified Integrated Cubic Phase Function

In order to compensate the DFM, a modified ICPF is proposed for CPS and applied to compensate the DFM.

4.1.

Modified Integrated Cubic Phase Function

Consider a monocomponent CPS expressed as follows:

(14)

s(tm)=ρexp[j2π(a1tm+12a2tm2+16a3tm3)],tm[T/2,T/2],
where ρ, a1, a2, and a3 denote the amplitude, the centroid frequency (CF), the CR, and the QCR, respectively. T represents the integration time.

The instantaneous frequency rate (IFR)22,23 of s(tm) can be stated as

(15)

IFR=d2Φ(tm)dtm2=2π(a2+a3tm),
where Φ(tm) denotes the phase of s(tm).

The CPF22,23 is defined as

(16)

CPF(tm,f)=τm2R(tm,τm)exp(j2πfτm2)dτm2=ρ1exp[j2π(2a1tm+a2tm2+13a3tm3)]δ(fa2a3tm),
where τm is a lag-time variable with respect to tm, f denotes the IFR domain, and δ(·) is the Dirac delta function. R(tm,τm) is the kernel function defined as

(17)

R(tm,τm)=s(tm+τm)s(tmτm)=ρ2exp[j2π(2a1tm+a2tm2+13a3tm3)]exp[j2π(a2+a3tm)τm2].

After the Fourier transform, the signal energy peaks along the inclined line fa2a3tm=0 in the time-IFR domain. The FT performed on Eq. (16) is with respect to τm2, which indicates that the unambiguous range of f is

(18)

|f|fs2/2,
where fs denotes the sampling frequency.

In realistic applications, the IFR of a CPS is usually much smaller than fs2/2.1819.20,2223.24 To obtain a higher estimation accuracy, we limit f to be |f|fs/2, and introduce a zoom factor to balance the resolution and range of the IFR. Thus, a modified CPF is defined as

(19)

MCPF(tm,fτ)=τm2R(tm,τm)exp(j2πfτητm2)dτm2=ρ1exp[j2π(2a1tm+a2tm2+13a3tm3)]δ(fτa2ηa3ηtm),fτ[fs2,fs2],
where η is the introduced zoom factor. The selection criterion will be discussed in Sec. 4.2.

Because the sampling corresponds to τm2 is nonuniform, the Fourier transform performed on Eq. (19) cannot be implemented by the FFT. Fortunately, the nonuniform FFT (NUFFT) can be applied to speed up the FT with the computational cost of O(M2log2M),24,25 where M is the number of integration pulses.

After the NUFFT, the signal energy peaks along the inclined line

(20)

fτa2/ηa3tm/η=0.
Due to the coupling between tm and τm in the second exponential phase term of R(tm,τm), the inclined line is not parallel to the tm-axis. Due to the first exponential phase term of R(tm,τm), the SCFT18,26 cannot remove this coupling. It is observed from Eq. (19) that, after the NUFFT operation, the CR and QCR have been packaged into the inclined line, i.e., Eq. (20).

The Radon transform,27 the Hough transform,28 and the modified Hough transform29 are the effective methods to accomplish energy accumulation along a straight line. Nevertheless, a 2-D searching procedure is inevitable which will burden the system complexity and increase the computational cost.

To realize energy accumulation along the inclined line without brute-force searching, the IFFT is performed on |MCPF(tm,fτ)|2 with respect to fτ

(21)

R1(tm,τm)=IFFTfτ[|MCPF(tm,fτ)|2]=ρ2exp[j2π(a2ητm+a3ητmtm)],
where τm denotes the lag-time variable with respect to fτ. IFFTfτ[·] denotes the IFFT operation along the fτ-axis. The first exponential phase term of R(tm,τm) does not appear in R1(tm,τm), which indicates that the SCFT becomes applicable to remove the coupling between tm and τm. Applying the SCFT on Eq. (21) along the tm-axis, we obtain

(22)

R2(ft,τm)=SCFT(ξτmtm)[R1(tm,τm)]=ρ3exp(j2πa2ητm)δ(fta3ηξ),
where ft is the frequency domain with respect to τmtm, SCFT(ξτmtm)[·] denotes the SCFT operation along the tm-axis, ξ is a zoom factor which is introduced to balance the resolution and range of the QCR, and its selection criterion will be analyzed in Sec. 4.2.

After the SCFT, the signal energy distributes along the beeline fta3/ηξ=0 in the τmft domain. Thereafter, applying the FFT along the τm-axis to achieve further energy accumulation, we obtain

(23)

MICPF(ft,fτ)=FFTτm[R2(ft,τm)]=ρ4δ(fτa2η)δ(fta3ηξ),
where FFTτm[·] denotes the FFT operation along the τm-axis.

A sole peak has been accumulated at the point [a2/η,a3/(ηξ)] in the fτft domain. Thereafter, a peak detection technique1819.20,24 is employed to estimate the CR and QCR of the CPS.

From the above discussion, the abbreviated expression of the proposed MICPF can be expressed as

(24)

MICPF(ft,fτ)=FFTτm{SCFT(ξτmtm){IFFTfτ|NUFFT(ητm2)[KF(s(tm))]|2}},
where KF[·] denotes the kernel function defined in Eq. (17).

4.2.

Selection Criteria of Zoom Factors

Equation (23) indicates that the CR and QCR of the CPS can be estimated as

(25)

a2=ηf^τ,|f^τ|fs/2,

(26)

a3=ηξf^t,|f^t|fs/2,
where (f^t,f^τ) denotes the peak location. Furthermore, the estimated ranges of the CR and QCR are

(27)

|a2|ηfs/2,

(28)

|a3|ηξfs/2.
To avoid spectrum aliasing, it is assumed that |a2|fs/T, |a3|(9fs)/(2T2),22,23 thus we have

(29)

|a2+a3tm|13fs/(4T).
In order to avoid energy loss, Eq. (19) indicates that η should satisfy

(30)

|(a2+a3tm)/η|fs/2.
Therefore, η should be

(31)

η6.5/T.

Equations (27) and (31) indicate that the estimated range of the CR is larger than its unambiguous range, so it is suitable for the CR estimation. Similarly, the estimated range of the QCR should be larger than its unambiguous range, which indicates that ξ should be

(32)

ξ9/(ηT2).

However, a bigger zoom factor will result in resolution degradation. Therefore, we suggest η=6.5/T, ξ=18/(13T) for the MICPF.

In realistic applications, the CR and QCR of a CPS may be limited to a smaller range. Under this situation, the zoom factors can be selected in terms of the real ranges of the parameters.

4.3.

Performance Analyses

The cross-term suppression performance, the computational complexity, and the antinoise performance are analyzed to evaluate the performance of the proposed MICPF. The SCFT-based method,18 the SS-ICPF,24 and the method in Ref. 30 are also presented for comparison. These are chosen due to the following: (1) these three approaches are the typical methods for the parameters estimation of the CPS, (2) no or few redundant data is necessary, and (3) the CR and QCR are estimated simultaneously, which can suppress the error propagation.

4.3.1.

Cross-term suppression performance

The bilinearity of the MICPF indicates that cross-terms will appear for multi-CPS. After the NUFFT, the cross-terms in the time-IFR domain take the cosine oscillation structure.24 Thus, the cross-terms cannot be accumulated as the self-terms, because: (1) MICPF can achieve energy accumulation along the inclined line with where signal energy peaks, and (2) the cosine oscillation structure will spread the cross-terms energy. This conclusion can be proven by the experiment performed in Sec. 6.1.

However, if the amplitudes of different CPSs differ significantly, the weak CPSs may be submerged by the cross-terms generated by the strong CPSs. Under this situation, similar to other estimation methods,1819.20 the CLEAN technique31 can be employed to separate the weak and strong CPSs, then the strong CPSs and weak CPSs can be estimated iteratively.

4.3.2.

Computational complexity

In what follows, the computational complexity is analyzed in terms of the number of operations, i.e., complex multiplications (Mc) and complex additions (Ac). Denote the sampling points of signal by M. The main procedures of the MICPF include: the NUFFT-based FT performed on Eq. (19) [O(M2log2M) Mc and O(M2log2M) Ac], the IFFT applied on Eq. (21) [O(M2log2M) Mc and O(M2log2M) Ac], the chirp-z-based SCFT taken of Eq. (22) [O(M2log2M) Mc and O(M2log2M) Ac], and the FFT performed on Eq. (23) [O(M2log2M) Mc and O(M2log2M) Ac]. On the other hand, O(M2log2M) Mc and O(M3) Ac are needed for the SS-ICPF,24 and O(M2log2M) Mc and O(M2log2M) Ac are needed for the SCFT-based method18 and the method in Ref. 30. The computational costs of the MICPF, the SS-ICPF, the SCFT-based method, and the method in Ref. 30 are listed in Table 1.

Table 1

Computational cost.

MethodsMultiplicationsAdditions
MICPFO(M2 log2 M)O(M2 log2 M)
SS-ICPFO(M2 log2 M)O(M3)
SCFT-basedO(M2 log2 M)O(M2 log2 M)
Ref. O(M2 log2 M)O(M2 log2 M)

It can be seen from Table 1 that the computational complexity of the proposed method is lower than that of the SS-ICPF. Moreover, for the SS-ICPF, the space switch method which includes a circular shifting operation and map operation will burden the complexity of the estimator.24,29 The computational cost of the MICPF, the SCFT-based method, and the method in Ref. 30 are in the same order. However, the antinoise performance of the MICPF is better than that of the SCFT-based method and the method in Ref. 30, which will be analyzed in Sec. 4.3.3.

4.3.3.

Antinoise performance

First-order perturbation analysis as used in Ref. 21 is an effective method to analyze the antinoise performance of estimation algorithms. However, due to the NUFFT operation, the modulus square operation, the IFFT operation, the SCFT operation, and the FFT operation that are performed on R(tm,τm), the first-order perturbation analysis becomes extremely difficult. Similar to Refs. 1819.20, the mean square error (MSE) is adopted to analyze the antinoise performance of MICPF in this section.

Example: Consider a mono-CPS contaminated with a zero-mean complex white Gaussian noise, where the sampling frequency and sampling points of the signal are 128 Hz and 256, respectively. The CF, CR, and QCR of the CPS are set as 40 Hz, 60  Hz/s, and 50  Hz/s2, respectively, and the test input SNRs are [1010] (dB). 100 iterations of Monte-Carlo experiments are performed for each input SNR value. The zoom factors η and ξ are set as 2 and 0.5, respectively. It is worthwhile noting that the zoom factors are set in terms of the parameters’ actual values. However, it is not in conflict with the selection criteria described in Sec. 4.2. Moreover, in Sec. 6, the zoom factors are also set in terms of the real ranges of the targets’ motion parameters.

The simulation results and the corresponding Cramer–Rao bounds (CRBs) are shown in Fig. 1, and the expressions of the CRBs can be found in Ref. 32. In Fig. 1, the threshold SNR of the MICPF is 8  dB, because: (1) the bilinear kernel function is applied in the MICPF and (2) signal energy is well accumulated by the proposed nonsearching method. The SS-ICPF has the same antinoise performance as the MICPF, but with a greater computational cost. However, the threshold SNR of the SCFT-based method is 3  dB due to the fourth-order nonlinear kernel function being applied. Compared with the SCFT-based method, the method in Ref. 30 introduces a constant delay to suppress noise, and achieves a 1 dB antinoise performance gain. Nevertheless, it still has a 4 dB loss compared with the MICPF.

Fig. 1

Parameters estimation performance: (a) MSE of the CR estimation and (b) MSE of the quadratic CR estimation.

JARS_10_3_035009_f001.png

From the above analyses and experiments, we can conclude that the proposed method is more suitable for CPS than the SS-ICPF, the SCFT-based method, and the method in Ref. 30.

5.

Target Detection Based on Second-Order Keystone Transform Modified Integrated Cubic Phase Function

In this section, by employing the SKT-based RM correction method and the MICPF, a coherent detection algorithm, known as SKT-MICPF, is presented to realize the detection of a high-speed maneuvering target with a constant radial jerk.

5.1.

Target Detection Based on Second-Order Keystone Transform Modified Integrated Cubic Phase Function

For simplicity, a single target is considered in the scene. After the RM correction with the proposed method in Sec. 3, the azimuth echoes of a range cell where a target is located can be stated as

(33)

s(tn)=A1exp{j2π2[r+(vTvs)tn+atn2/2+ktn3/6]λ}.

Performing the MICPF on s(tn), we get

(34)

MICPF(ft,fτ)=A2δ(fτ+2aλη)δ(ft+2kληξ).

Target’s radial acceleration and jerk can be estimated via peak detection technique. Thus, a compensation function is constructed to compensate the DFM, which takes the form

(35)

H2(a^,k^,tn)=exp[j2π(a^λtn2+k^3λtn3)],
where a^ and k^ denote the estimated acceleration and jerk, respectively.

After the DFM compensation, the FFT is applied to achieve coherent integration

(36)

s(fd)=FFTtn{s(tn)×H2(a^,k^,tn)}=A3exp(j2π2rλ)δ(fd+2Vλ),
where V=vTvs.

The processing procedures of the proposed method, i.e., velocity compensation, SKT, and DFM compensation, will not change the distributions. Therefore, the CFAR technique is utilized to target detection, and the threshold can be set as33

(37)

T^=Nw(PFA1/Nw1)β^2,
where Nw, PFA, and β^2 denote the length of the reference windows, the probability of false alarm, and the estimated noise power, respectively.

It is worth noting that the strong target may be detected even when the RM has not been completely corrected, i.e., different search velocities yield positive detections. Under this situation, we confirm that the target’s radial velocity has been matched by the searching velocity, if both of the following two criteria are satisfied: (1) the searching velocity satisfies |V^|Δvs/2, where V^ denotes the estimated value of V; and (2) the searching velocity corresponds to the best integration result (the maximum integrated peak).

By the criterion 1, the range of searching velocity has been confined to vsn=vs0+nva(n=0,±1,±2,), where vs0 denotes the searching velocity which matches the target’s radial velocity vT, i.e., |vTvs0|Δvs/2. Thereafter, the criterion 2 is introduced to select vs0 among vsn=vs0+nva(n=0,±1,±2,). On one hand, the RM can be completely corrected when the searching velocity is vs0, then optimal integration can be obtained. On the other hand, due to the long-time integration, vaT>ΔR can be easily satisfied, which indicates that the RM cannot be completely removed when the searching velocities are vsn=vs0+nva(n=±1,±2,), and will result in integration loss. Therefore, criterion 2 implies that the final matched velocity can be selected from the set of vsn.

Furthermore, the target’s radial velocity can be estimated as

(38)

v^T=vs0+V^.

Although a single target is considered in the above analysis, SKT-MICPF has a good detection performance for multitargets since it can distinguish multitargets from the initial range, the radial velocity, acceleration, and jerk.

5.2.

Procedure of Second-Order Keystone Transform Modified Integrated Cubic Phase Function

The flowchart of the SKT-MICPF is shown in Fig. 2, which is mainly composed of seven steps. I and N denote the number of searching velocities and range cells, respectively.

Fig. 2

Flowchart of the SKT-MICPF.

JARS_10_3_035009_f002.png
  • Step 1: Complete the downconversion, pulse compression, and range FFT. Initialize the range of the searching velocity, let i=1.

  • Step 2: Construct the i’th compensation function H1(tm,fr,vsi) as Eq. (6), and multiply H1(tm,fr,vsi) with Sc(tm,fr). Thereafter, the SKT and range IFFT are performed successively. Initialize the range of searching range cell, let n=1.

  • Step 3: Get the data s(tn) of the n’th range cell, and perform the MICPF on s(tn). Carry out the peak detection in the ftfτ domain. If a peak has been detected, estimate the acceleration and jerk, and then go to step 4. Otherwise, go to step 5.

  • Step 4: Construct a compensation function H2(a^,k^,tn) as Eq. (35), and multiply H2(a^,k^,tn) with s(tn). Accomplish coherent integration and target detection via the FFT and CFAR detection, successively. If a target has been detected, estimate V, and note down the coordinate, then go to step 5. Otherwise, go directly to step 5.

  • Step 5: If n<N, set n=n+1, and repeat step 3 and step 4 until n=N.

  • Step 6: If i<I, set i=i+1, and repeat step 2 to step 5 until i=I.

  • Step 7: After all the searching velocities and range cells have been detected, estimate targets’ radial velocities by the two criteria in Sec. 5.1.

5.3.

Analysis of Computational Complexity

Denote the number of searching velocities, integration pulses, and range cells by I, M, and N, respectively. Assume that the number of searching angles of IAR-DCFT is I. The main procedures of the SKT-MICPF include: I times chirp-z-based SKT [O(MNIlog2M)], MI times N points IFFT [O(MNIlog2N)], and NI times MICPF [O(M2NIlog2M)]. Therefore, the computational cost of the SKT-MICPF is O(M2NIlog2M). On the other hand, the computational cost of the IAR-DCFT is O(M3NIlog2M).15

From the above analysis, we can see that the computational complexity of the proposed method is much less than that of the IAR-DCFT. In real applications, the initial slant range and radial velocity of a target are limited to some fixed region via prior knowledge,6 and the computational cost of the proposed method can be further decreased.

Remark: Although the SKT-MICPF is proposed for high-speed maneuvering targets’ detection, it can also be applied to inverse synthetic aperture radar (ISAR) imaging. For ISAR imaging of a target with complex motion, after the range alignment and the phase adjustment, the azimuth echoes of a range cell can be modeled as multi-CPS.1718.19.20,24,30 If the parameters of those CPSs can be estimated, high-resolution imaging can be obtained. Moreover, the performance of parameters’ estimation has a great impact on the imaging quality. The proposed MICPF can achieve high-precision estimation for multi-CPS. Comparisons with the SS-ICPF, the SCFT-based method, and the method in Ref. 30 show that the proposed MICPF outperforms those methods, which indicates that the proposed MICPF can also be applied to ISAR imaging, and may achieve a better imaging quality than those methods.

6.

Verification of the Proposed Algorithm

We are to validate the effectiveness of the proposed detection algorithm with synthetic radar data and real radar data.

6.1.

Verification with Synthetic Radar Data

Two high-speed maneuvering targets with a constant radial jerk are contained in the radar echoes. The radar parameters and targets’ motion parameters are listed in Tables 2 and 3, respectively. The echoes are contaminated with a zero-mean complex white Gaussian noise.

Table 2

Radar parameters of the synthetic data.

Parameter (unit)Values
Carrier frequency (GHz)1
Pulse width (50
Bandwidth (MHz)10
Sample frequency (MHz)10
PRF (Hz)1000
Number of integrated pulses2000

Table 3

Targets’ motion parameters.

TargetRange cell indexRadial velocity (m/s)Radial acceleration (m2/s)Radial jerk (m3/s)
A200209745−36
B2002106−4231.5

Assume the ranges of the targets’ radial acceleration and jerk are limited to [50,50]  m/s2 and [40,40]  m/s3, respectively. The zoom factor of the MICPF can be set as η=1.2, ξ=0.5 in terms of the radar parameters and the ranges of the target’s motion parameters. According to the above parameters, the estimated resolutions of acceleration and jerk are 0.09  m/s2 and 0.045  m/s3, respectively. The SNRs (after the pulse compression) of the two targets are SNRA=3  dB and SNRB=6  dB. The interval of the searching velocity is set as 15  m/s for the RM correction in terms of Eq. (13).

Figure 3 shows the simulation results of the SKT-MICPF. The result after the pulse compression is given in Fig. 3(a). Because of the low SNR, the two targets are buried in the noise. Figure 3(b) shows the result after the pulse compression without noise. Due to the high-speed, serious RMs occur to the two targets. The results after the velocity compensation (2100  m/s) and the SKT are shown in Figs. 3(c) and 3(d), respectively. The RMs of the two targets have been corrected. Although the RMs have been corrected, the DFMs still remain, which make the targets’ energy spread in the Doppler frequency domain as shown in Fig. 3(e). The distribution of the targets’ energy in the slow time-IFR domain is shown in Fig. 3(f). The two targets peak along two inclined lines in the slow time-IFR domain, respectively. Figure 3(g) gives the accumulation result via the proposed method. Two peaks have been accumulated and well distinguished in the acceleration-jerk domain, whereas the cross-term is not accumulated. Thus, the targets’ radial accelerations and jerks can be estimated by the peak detection technique. Figures 3(h) and 3(i) show the coherent integration results after the RM and DFM compensation. Targets are well-concentrated and can be easily detected via the CFAR technique.

Fig. 3

Simulation results of the SKT-MICPF [(b)–(e) without noise]: (a) result after the pulse compression, (b) result after the pulse compression without noise, (c) result after the velocity compensation (the searching velocity is 2100  m/s), (d) result after the SKT, (e) result after the FFT performed on the 200th range cell, (f) the distribution of targets energy in the slow time-IFR domain, (g) result of the MICPF, (h) coherent integration result of target A after the RM and DFM compensation, and (i) coherent integration result of target B after the RM and DFM compensation.

JARS_10_3_035009_f003.png

If there is another target (target C) with a radial velocity, e.g., 2000  m/s, it is also located in the 200th range cell. Assume that the searching velocities are set as [1800152400]  m/s. When the searching velocity is 1995  m/s, the RM of target C can be completely compensated, whereas serious RMs still remain for targets A and B. Therefore, target C can be detected first. When the searching velocity is 2100  m/s, the RMs of targets A and B can be eliminated further to parameters estimation and targets detection.

The integration results of the RFT,2 RLVD,7 IAR-DCFT,15 and the proposed SKT-MICPF are further performed for comparison. Consider target A in the scene, where the SNR after the pulse compression is set as 13  dB.

Figure 4(a) shows the integration result via the RFT. Because of the RM and DFM induced by the target’s radial acceleration and jerk, RFT becomes invalid. The integration result of the RLVD is given in Fig. 4(b). The RM induced by the target’s radial velocity and acceleration and the DFM caused by the target’s radial acceleration can be eliminated by the RLVD. However, the DFM induced by the target’s radial jerk affects the integration performance. The integration results of the IAR-DCFT and SKT-MICPF are shown in Figs. 4(c) and 4(d), respectively. The target is well-focused via the two approaches, which indicates that the IAR-DCFT and SKT-MICPF are suitable for the detection of a high-speed maneuvering target with a jerky motion. However, the SKT-MICPF has a lower computational complexity, which has been analyzed in Sec. 5.3. The time costs of the SKT-MICPF and IAR-DCFT are also given in Table 4.

Fig. 4

Integration results of the RFT, RLVD, IAR-DCFT, and SKT-MICPF: (a) integration result of the RFT, (b) integration result of the RLVD, (c) integration result of the IAR-DCFT, and (d) integration result of the SKT-MICPF.

JARS_10_3_035009_f004.png

Table 4

Time cost of the SKT-MICPF and IAR-DCFT.

AlgorithmsTime cost (s)
SKT-MICPF34.4294
IAR-DCFT1682.7359

Main configuration of the computer: CPU: Intel Core i7-4770 3.4 GHz;

RAM: 8.00 GHz; operating system: Windows 7; software: MATLAB® 2012a.

The detection performances of the RFT, RLVD, IAR-DCFT, and SKT-MICPF are further evaluated by Monte-Carlo experiments, and target A is considered in the scene. The false alarm rate is set to be 106, the test input SNRs are [25120] (dB, after the pulse compression), and 100 Monte-Carlo experiments are performed for each input SNR value. The simulation results are shown in Fig. 5.

Fig. 5

Detection probability curves against input SNR.

JARS_10_3_035009_f005.png

It is observed from Fig. 5 that, due to its ability to deal with the DFM effect induced by the target’s radial jerk, the detection performance of the proposed method is superior to those of the RFT and RLVD. The IAR-DCFT has a better detection performance (2 dB) than the SKT-MICPF. This is because: (1) the DCFT has the optimal antinoise performance16,17 and (2) the MICPF suffers from SNR loss due to its bilinearity. However, the computational cost of the IAR-DCFT is much greater than that of the SKT-MICPF.

From the above experiments, we can see that the SKT-MICPF offers a much better performance than the RFT and RLVD, while offering a much lower computational cost than the higher performing IAR-DCFT. This makes the proposed approach a stronger candidate for practical adoptions.

6.2.

Verification with Real Radar Data

The real radar data were collected in February 2010 by the National University of Defense Technology, Hunan, China. The parameters of the radar system are listed in Table 5. The 2000 integration pulses and 200 range cells are extracted from the collected data.

Table 5

Radar parameters of the real data.

Parameter (unit)Values
Wave bandL
Pulse width (75
Bandwidth (MHz)2
Sample frequency (MHz)3.25
PRT (910

According to Eq. (13), the interval of the searching velocity is set as 30  m/s for the RM correction. The processing results are shown in Fig. 6.

Fig. 6

Processing results of real radar data via the SKT-MICPF: (a) radar echoes, (b) pulse compression result, (c) RM correction result (the searching velocity is 180  m/s), (d) MICPF result, (e) integration result via the MTD (the zoomed-in area: range cell: [117:125], Doppler cell: [877:857]), and (f) coherent integration result via the SKT-MICPF.

JARS_10_3_035009_f006.png

Figures 6(a) and 6(b) show the original radar echoes and the result after the pulse compression, respectively. It is obvious from Fig. 6(b) that the target envelope has been shifted away from its original position due to the RM. Figure 6(c) gives the result of RM correction via the proposed method (the searching velocity is 180  m/s). After the velocity compensation and SKT, the RM has been corrected. The result after the MICPF performed on the 121th range cell is shown in Fig. 6(d). We can see that a peak has been accumulated, and then the target’s radial acceleration and jerk can be estimated by the peak detection technique. Figure 6(e) gives the integration result via the MTD, and the nearby area of the integrated peak is also zoomed-in to evaluate the integration performance. Due to the RM and DFM, target energy distributes in both the range and Doppler directions. The integration result of the proposed method is shown in Fig. 6(f). Because the proposed method can compensate the RM and DFM, the integrated peak of the proposed SKT-MICPF is larger than that of the MTD.

Although the maneuverability and velocity of the passenger plane are low, the RM and DFM affect the integration performance greatly. The advantages of the SKT-MICPF will be more obvious for a target moving at a higher speed and maneuverability.

7.

Conclusions

In this paper, an SKT-based RM correction method is first proposed then an MICPF is presented for the parameters’ estimation of CPS. Thereafter, by using the proposed RM correction method and MICPF, a coherent target detection algorithm, known as SKT-MICPF, is proposed for a high-speed maneuvering target with a constant jerkiness. The advantages of the proposed method include the following: (1) the estimations of a target’s radial acceleration and jerk are not searched which greatly reducing the computational cost; (2) the RM correction and DFM compensation can be accomplished in a low SNR scene, which indicates the proposed method can get a high antinoise performance; and 3) it has a good detection performance for multitargets since it can distinguish multitargets from the initial slant range, radial velocity, radial acceleration, and radial jerk. Comparisons with several methods show that the SKT-MICPF offers a much better performance than the RFT and RLVD, while offering a much lower computational cost than the higher performing IAR-DCFT. This makes the proposed approach a stronger candidate in realistic applications.

Acknowledgments

This work was partially supported by the National Natural Science Foundation of China (61001204) and the Fundamental Research Funds for the Central Universities (JY0000902020). Jiancheng Zhang and Tao Su proposed the target detection algorithms and wrote this paper; Jiancheng Zhang and Yanyan Li performed the experiments; Jibin Zheng analyzed the data and revised the article.

References

1. 

M. Xing et al., “New parameter estimation and detection algorithm for high speed small target,” IEEE Trans. Aerosp. Electron. Syst. 47(1), 214–224 (2011).IEARAX0018-9251http://dx.doi.org/10.1109/TAES.2011.5705671Google Scholar

2. 

J. Xu et al., “Radon-Fourier transform for radar target detection I: generalized Doppler filter bank,” IEEE Trans. Aerosp. Electron. Syst. 47(2), 1186–1202 (2011).IEARAX0018-9251http://dx.doi.org/10.1109/TAES.2011.5751251Google Scholar

3. 

D. Zhu, Y. Li and Z. Zhu, “A keystone transform without interpolation for SAR ground moving-target imaging,” IEEE Geosci. Remote Sens. Lett. 4(1), 18–22 (2007).http://dx.doi.org/10.1109/LGRS.2006.882147Google Scholar

4. 

X. Rao et al., “Axis rotation MTD algorithm for weak target detection,” Digital Signal Process. 26, 81–86 (2014).DSPREJ1051-2004http://dx.doi.org/10.1016/j.dsp.2013.12.003Google Scholar

5. 

J. Zheng et al., “Radar high-speed target detection based on the scaled inverse Fourier transform,” IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 8(3), 1108–1119 (2015).IJSQEN1077-260Xhttp://dx.doi.org/10.1109/JSTARS.2014.2368174Google Scholar

6. 

X. Chen et al., “Maneuvering target detection via Radon-fractional Fourier transform-based long-time coherent integration,” IEEE Trans. Signal Process. 62(4), 939–953 (2014).ITPRED1053-587Xhttp://dx.doi.org/10.1109/TSP.2013.2297682Google Scholar

7. 

X. Li et al., “Coherent integration for maneuvering target detection based on Radon-Lv’s distribution,” IEEE Signal Process. Lett. 22(9), 1467–1471 (2015).IESPEJ1070-9908http://dx.doi.org/10.1109/LSP.2015.2390777Google Scholar

8. 

J. Tian et al., “Joint estimation algorithm for multi-targets’ motion parameters,” IET Radar Sonar Navig. 8(8), 939–945 (2014).http://dx.doi.org/10.1049/iet-rsn.2013.0346Google Scholar

9. 

S. Zhu et al., “A new method for radar high-speed maneuvering weak target detection and imaging,” IEEE Geosci. Remote Sens. Lett. 11(7), 1175–1179 (2014).http://dx.doi.org/10.1109/LGRS.2013.2283887Google Scholar

10. 

J. Tian, W. Cui and S. Wu, “A novel method for parameter estimation of space moving targets,” IEEE Geosci. Remote Sens. Lett. 11(2), 389–393 (2014).http://dx.doi.org/10.1109/LGRS.2013.2263332Google Scholar

11. 

L. Kong et al., “Coherent integration algorithm for a maneuvering target with high-order range migration,” IEEE Trans. Signal Process. 63(17), 4474–4486 (2015).ITPRED1053-587Xhttp://dx.doi.org/10.1109/TSP.2015.2437844Google Scholar

12. 

P. Huang et al., “Approach for space-based radar manoeuvring target detection and high-order motion parameter estimation,” IET Radar Sonar Navig. 9(6), 732–741 (2015).http://dx.doi.org/10.1049/iet-rsn.2014.0192Google Scholar

13. 

X. Li et al., “A low complexity coherent integration method for maneuvering target detection,” Digital Signal Process. 49, 137–147 (2016).DSPREJ1051-2004http://dx.doi.org/10.1016/j.dsp.2015.10.008Google Scholar

14. 

X. Li et al., “A fast detection method for maneuvering target in coherent radar,” IEEE Sensors J. 15(11), 6722–6729 (2015).ISJEAZ1530-437Xhttp://dx.doi.org/10.1109/JSEN.2015.2464788Google Scholar

15. 

X. Rao et al., “Long-time coherent integration detection of weak manoeuvring target via integration algorithm, improved axis rotation discrete chirp-Fourier transform,” IET Radar Sonar Navig. 9(7), 917–926 (2015).http://dx.doi.org/10.1049/iet-rsn.2014.0344Google Scholar

16. 

X. Xia, “Discrete chirp-Fourier transform and its application to chirp rate estimation,” IEEE Trans. Signal Process. 48(11), 3122–3133 (2000).ITPRED1053-587Xhttp://dx.doi.org/10.1109/78.875469Google Scholar

17. 

L. Wu et al., “ISAR imaging of target with complex motion based on discrete chirp Fourier transform for cubic chirp,” IEEE Trans. Geosci. Remote Sens. 50(10), 4201–4212 (2012).IGRSD20196-2892http://dx.doi.org/10.1109/TGRS.2012.2189220Google Scholar

18. 

X. Bai et al., “ISAR imaging of a ship target based on parameter estimation of multi-component quadratic frequency-modulated signals,” IEEE Trans. Geosci. Remote Sens. 52(2), 1418–1429 (2014).IGRSD20196-2892http://dx.doi.org/10.1109/TGRS.2013.2251348Google Scholar

19. 

J. Zheng et al., “ISAR imaging of targets with complex motions based on the keystone time-chirp rate distribution,” IEEE Geosci. Remote Sens. Lett. 11(7), 1275–1279 (2014).http://dx.doi.org/10.1109/LGRS.2013.2291992Google Scholar

20. 

Y. Li, T. Su and J. Zheng, “Inverse synthetic aperture radar imaging of targets with complex motion based on modified chirp rate-quadratic chirp rate for cubic phase signal,” J. Appl. Remote Sens. 9, 095036 (2015).http://dx.doi.org/10.1117/1.JRS.9.095036Google Scholar

21. 

P. Wang et al., “Integrated cubic phase function for Linear FM signal analysis,” IEEE Trans. Aerosp. Electron. Syst. 46(3), 963–977 (2010).IEARAX0018-9251http://dx.doi.org/10.1109/TAES.2010.5545167Google Scholar

22. 

P. O’shea, “A fast algorithm for estimating the parameters of a quadratic FM signal,” IEEE Trans. Signal Process. 52(2), 385–393 (2004).ITPRED1053-587Xhttp://dx.doi.org/10.1109/TSP.2003.821097Google Scholar

23. 

P. O’shea, “A new technique for instantaneous frequency rate estimation,” IEEE Signal Process. Lett. 9(8), 251–252 (2002).IESPEJ1070-9908http://dx.doi.org/10.1109/LSP.2002.803003Google Scholar

24. 

J. Zheng et al., “ISAR imaging for fluctuating ships based on a fast bilinear parameter estimation algorithm,” IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 8(8), 3954–3966 (2015).IJSQEN1077-260Xhttp://dx.doi.org/10.1109/JSTARS.2015.2440911Google Scholar

25. 

Q. Liu and N. Nguyen, “An accurate algorithm for non uniform fast Fourier transforms (NUFFT’s),” IEEE Microwave Guided Wave Lett. 8(1), 18–20 (1998).IMGLE31051-8207http://dx.doi.org/10.1109/75.650975Google Scholar

26. 

Y. Li et al., “Inverse synthetic aperture radar imaging of targets with nonsevere maneuverability based on the centroid frequency chirp rate distribution,” J. Appl. Remote Sens. 8, 095065 (2015).http://dx.doi.org/10.1117/1.JRS.9.095065Google Scholar

27. 

J. Carretero-moya et al., “Application of the Radon transform to detect small-targets in sea clutter,” IET Radar Sonar Navig. 3(2), 155–166 (2009).http://dx.doi.org/10.1049/iet-rsn:20080123Google Scholar

28. 

B. Carlso, E. D. Evans and S. L. Wilson, “Search radar detection and track with the Hough transform,” IEEE Trans. Aerosp. Electron. Syst. 30(1), 102–108 (1992).IEARAX0018-9251http://dx.doi.org/10.1109/7.250410Google Scholar

29. 

J. Zeng et al., “Modified Hough transform for searching radar detection,” IEEE Geosci. Remote Sens. Lett. 5(4), 683–686 (2008).http://dx.doi.org/10.1109/LGRS.2008.2002574Google Scholar

30. 

J. Zheng et al., “ISAR imaging of targets with complex motions based on a noise-resistant parameter estimation algorithm without nonuniform axis,” IEEE Sensors J. 8(16), 2509–2518 (2016).ISJEAZ1530-437Xhttp://dx.doi.org/10.1109/JSEN.2016.2516040Google Scholar

31. 

J. Misiurewicz et al., “Radar detection of helicopters with application of CLEAN method,” IEEE Trans. Aerosp. Electron. Syst. 48(4), 3525–3537 (2012).IEARAX0018-9251http://dx.doi.org/10.1109/TAES.2012.7324734Google Scholar

32. 

B. Ristic and B. Boashash, “Comments on “The Cramer-Rao lower bounds for signals with constant amplitude and polynomial phase”,” IEEE Trans. Signal Process. 46(6), 1708–1709 (1998).ITPRED1053-587Xhttp://dx.doi.org/10.1109/78.678501Google Scholar

33. 

M. Richards, Fundamentals of Radar Signal Processing, pp. 352–354, McGraw-Hill, New York (2005).Google Scholar

Biography

Jiancheng Zhang received his BS degree in measurement and control technology and instrumentation from Xidian University, Shaanxi, China, in 2011. He is currently pursuing his PhD at the National Laboratory of Radar Signal Processing, Xidian University. His research interests include maneuvering target detection, parameter estimation, and time-frequency analysis.

Tao Su received his BS in information theory, MS degree in mobile communication, and PhD in signal and information processing from Xidian University, Xi’an, China, in 1990, 1993, and 1999, respectively. He is currently a professor at the National Laboratory of Radar Signal Processing, Xidian University. His research interests include high-speed real-time signal processing on radar, sonar and telecommunications, digital signal processing, parallel processing system design, and FPGA IP design.

Yanyan Li received her BS and MS degrees in industrial design and signal and information processing from Xidian University, Xi’an, China, in 2008 and 2011, respectively. She is currently pursuing her PhD at the National Laboratory of Radar Signal Processing, Xidian University. Her research interests include ISAR imaging, target detection, and target tracking.

Jibin Zheng received his BS degree in electronic information science and technology from Shandong Normal University, Shandong, China, in 2009, and PhD in signal and information processing from Xidian University in 2015. Since 2015, he has been on the faculty at National Laboratory of Radar Signal Processing at Xidian University. His research interests include target detection, synthetic aperture radar (SAR) and inverse SAR signal processing, cognitive radar, and time-frequency analysis.

Jiancheng Zhang, Tao Su, Yanyan Li, Jibin Zheng, "Radar high-speed maneuvering target detection based on joint second-order keystone transform and modified integrated cubic phase function," Journal of Applied Remote Sensing 10(3), 035009 (3 August 2016). http://dx.doi.org/10.1117/1.JRS.10.035009
JOURNAL ARTICLE
17 PAGES


SHARE
KEYWORDS
Target detection

Radar

Design for manufacturing

Signal to noise ratio

Chromium

Detection and tracking algorithms

Integration

Back to Top