Abstract

Two-dimensional pulse-to-pulse canceller (TDPC) of ground clutter can effectively suppress the clutter along the clutter trace, and therefore the moving target detectability of the following space-time adaptive processing (STAP) algorithm can be improved after TDPC as the clutter prefilter. However, TDPC may greatly impair the energy of moving target when inaccurate knowledge is exploited, which is detrimental to target detection. Aiming at this problem, a robust two-dimensional pulse-to-pulse canceller (RTDPC) of ground clutter is proposed. In order to enhance the TDPC’s robustness with inaccurate radar system parameters, which are mainly the platform velocity and crab angle, the errors of estimated platform velocity and crab angle are taken as the prior knowledge and added into the design of the clutter filter coefficient matrix. By exploiting RTDPC as the clutter prefilter, the moving target detectability of the following nonadaptive detection algorithm or STAP algorithm can also be enhanced. The simulated and MCARM data are utilized to verify the clutter suppression performance of RTDPC with inaccurate platform velocity and crab angle.

1. Introduction

For the fast moving platform, the clutter received by the airborne radar needs to be suppressed from the joint spatial-temporal domain, in which the clutter and moving target are separated. Space-time adaptive processing (STAP) has long been viewed as an effective tool for clutter suppression and target detection in airborne radar since it was first proposed by Brennan and Reed [1]. Two main limitations including computational cost and training sample demanding prevent it from being extensively applied and motivate the development of suboptimal dimension/rank-reduced STAP algorithms. By far, many literatures have studied various STAP methods to reduce the computational cost and training sample demanding in adaptive processing and enhance its practicability [24].

The dimension-reduced STAP methods [512], such as auxiliary channel processing (ACP) [5], joint domain localized (JDL) [6], and post-Doppler processing algorithms [7, 8], can reduce the adaptive dimension and hence reduce the computational cost and the number of required training samples. On the other hand, the rank-reduced STAP methods [1320], such as principle components (PC) [13], cross spectrum method (CSM) [16], and diagonal loading (DL) [17], make the eigendecomposition to the clutter plus noise covariance matrix (CPNCM) and try to estimate the clutter subspace. Although the rank-reduced STAP methods exploit the low rank property of clutter and reduce the number of training samples to twice of the clutter rank, they are computational intensive due to the CPNCM eigendecomposition. Besides, the eigenvectors of CPNCM cannot be accurately estimated when the number of training samples is small. Some target detection algorithms in airborne radar based on the reduced-rank STAP are soon proposed [21, 22]. Nowadays, the knowledge-aided (KA) radar, which can use the prior knowledge such as the digital elevation map (DEM) and geospatial databases to enhance STAP performance with small training sample support, has been attracting increasing attentions of researchers and practitioners [2133]. In KA radar, various kinds of information are fused for predicting the clutter distribution. For example, the distribution of the ground clutter in the azimuth-Doppler domain can be determined in advance when the platform velocity and the crab angle, which can be obtained by the inertial navigation system (INS) and global positioning system (GPS), combined with radar operating parameters are known. References [34, 35] brought forward the method of two-dimension pulse-pulse canceller (TDPC) of ground clutter based on this kind of prior knowledge, which is verified to be an effective clutter suppression method for airborne radar.

TDPC is in fact a nonadaptive clutter filter that takes the knowledge of radar operating parameters and platform velocity. However, due to the complex flight environment, the platform velocity and crab angle cannot be always correctly estimated. The inaccurate airborne platform velocity and crab angle used in the design of TDPC lead to the impairment to moving target and degradation of clutter suppression ability. Aiming at this problem, a robust two-dimensional pulse-to-pulse canceller (RTDPC) of ground clutter is proposed. In order to enhance the TDPC’s robustness with inaccurate prior knowledge that are mainly the platform velocity and crab angle, the errors of estimated airborne platform velocity and crab angle are added into the design of RTDPC. The moving target detectability of the following nonadaptive detection algorithm or STAP algorithm can also be enhanced. The simulated and MCARM data are utilized to verify the clutter suppression performance of RTDPC with inaccurate platform velocity and crab angle.

This paper is organized as follows. The ground clutter model and the original TDPC are reviewed in Section 2. Section 3 brings forward the principle of RTDPC. Section 4 illustrates the RTDPC filter cascaded with spatial-temporal matching (STM) or dimension-reduced STAP method. In Section 5, some experiments including simulated data and real data show the performance of RTDPC. Finally, we make a conclusion in Section 6.

2. Ground Clutter Model and the Original TDPC

The configuration of the airborne radar system to be studied is shown in Figure 1, which demonstrates the antenna array moving at a constant speed in the positive x-direction. The antenna array consists of array elements with interelement space . The range cell under test (CUT) is uniformly divided into independent clutter patches without consideration of the earth curvature. The location of the clutter patch is described by its azimuth angle and elevation angle . Assuming that the airborne radar system transmits pulses on the repetition of in one coherent processing interval (CPI), the received clutter vector on the CUT at the pulse is given by [2, 3] the following:where is the normalized Doppler frequency with being the crab angle, is the echo amplitude, and is the spatial steering-vector and can be expressed asin which is the normalized spatial frequency, is the operating wavelength, and is the transpose of a vector. The spatial-temporal snapshot of the clutter can be expressed aswhere is the temporal steering-vector and is the white noise vector. One re-expresses (1) as the matrix-vector form aswhere

In order to design TDPC, a coefficient matrix is established to minimize the residual clutter energy between adjacent pulses, that is, as follows [35]:where is the F-norm. Since the most of the energy contained in each comes from the clutter echo, the clutter can hence be suppressed by minimizing the difference in (6). Let , we have thatwhich can be further described by the Cauthy–Shwarz in equation aswhere is a constant. Now, our goal is to minimize . The following relation can be obtained after expanding :where is the trace of a matrix and is the conjugate transpose of a matrix or a vector. Let the gradient of with respect to be zero. Then, the solution to (9) is

The TDPC coefficient matrix can be obtained if the matrices and are known, where and are only determined by the parameters defined in (5). Thus, if the parameters, such as platform velocity, crab angle, operating wavelength, and so on, are accurately known in advance, the matrices and can be utilized as the prior information for computing . In other words, the TDPC coefficient matrix can be precalculated offline.

A prefiltering matrix of TDPC is devised aswhere is a identity matrix. Then, TDPC can suppress the clutter by the following equation:

3. Design of RTDPC

From the aforementioned designing procedure, the effectiveness of TDPC depends heavily on the accuracy of prior knowledge. In fact, though some parameters, such as operating wavelength and elevation angle of CUT, can be accurately known, some other parameters, such as platform velocity and crab angle, can only be obtained with errors due to the inconstant aircraft control. As a consequence, inaccurate parameters exploited in the design will certainly influence the clutter filtering result of TDPC, and hence in this section, a robust TDPC based on the steering-vector selection will be designed to improve the effectiveness of the original TDPC with inaccurate parameters (mainly platform velocity and crab angle). To start with, the estimated platform velocity error and estimated crab angle error are assumed to be and . Accordingly, the estimated normalized Doppler frequency can be expressed aswhere represents the real angle of a clutter patch with being the azimuth angle of a clutter patch and representing the crab angle. represents the estimated crab angle error, namely, the error generated the estimation. For the first two items and last two items in (13), they can be further simplified as

For a typical airborne radar system, the accuracies for INU and GPS measuring velocity and crab angle are about and , respectively. Generally, the velocity of airborne radar is over hundreds of meters per second. Therefore, reasonably speaking, should be a very small number and can be neglected in the Doppler calculation equation. Meanwhile, by using the facts that and when is very small, equation (14) can be further simplified aswhich becomes a equation related only with estimated crab angle, and the error between the true Doppler frequency and the estimated normalized Doppler frequency is . For verifying the feasibility of equation (15), the following results are demonstrated in Figure 2, where the average estimated Doppler frequency error (defined as ) curves in a range cell for , , and by using equations (13) and (15) are separately plotted. Evidently, as the aforementioned analysis, the variation of with is negligible and the crab angle plays the dominated role in . Consequently, when is determined, the upper bound of can be determined as for , where represents the maximum normalized Doppler frequency.

Let , then the estimated normalized Doppler frequency is located in the range . Then, the normalized Doppler frequency span can be uniformly divided into samples, which are denoted as . Accordingly, a spatial-temporal steering-vector dictionary is defined aswhere is a full-column rank matrix, is the assumed number of clutter patches evenly distributed in the whole azimuth angle, and . After obtaining and the snapshot in the CUT , a search method, which is similar to OMP [36], will be exploited to iteratively find the most relevant steering-vector to from , and this subsection outlines the main steps.

Let , in the first iteration, the most relevant steering-vector is selected as . Then, the corresponding eigenvalue and normalized eigenvector can be expressed as

Then, the residual vector between and the first normalized eigenvector is . For consistently improving the CPNCM estimation accuracy, the iterative procedure will not be stopped until the termination condition or the maximum number of iterations is satisfied at the -th iteration, where is a constant that represents the threshold and the symbol represents the infinity norm. At the -th iteration, the most relevant steering-vector is selected as . As a consequence, the eigenvalue and the normalized eigenvector at the -th iteration are expressed as

It is worth noting that and are mutually orthogonal due to the Schmidt orthogonalization processing in the iteration. For better understanding, the aforementioned algorithm is summarized in Algorithm 1. There exists the contradiction between the performance and computational cost. Intuitively, should be large for maintaining the satisfied performance. However, large will lead to large computational cost. Fortunately, the rule for choosing has been discussed in the previous literatures [4, 27], and will assure the performance.

Input: , , threshold , and maximum iteration number
   For q = 1, 2, 3, ..., do
    Update and
      ,
      and .
    Stop if or is satisfied
Output: , ,..., and , ,..., 90

After obtaining the selected steering-vector from , the robust TDPC can be designed in the similar way as the original TDPC, namely,

Accordingly, in RTDPC, and

4. Cascaded Structure

The proposed RTDPC can be exploited as a prefilter before the STM or the classical STAP method to form the two-stage prefilter and detector. In the first stage, RTDPC as the prefilter can effectively filter out the clutter. The second stage, which is STM or STAP, can further suppress the residual clutter and detect the moving target. Here, we will briefly illustrate the two-stage prefilter and detector.

Assume that the spatial-temporal steering-vector of the target is expressed as , where and , respectively, are the temporal and spatial target steering-vectors. After the clutter prefilter by , the following moving target detection methods as the second stage are illustrated separately:(1)  STM. The RTDPC prefiltering operation on the target vector isCorrespondingly, the weight vector of STM iswhere is the Hadamard product implying the element-by-element operation of two vectors, and and are the temporal and spatial static window vectors adopted for attenuating the sidelobes, respectively.(2)Dimension-Reduced STAP. Suppose that the corresponding dimension-reduced matrix for STAP is , then a dimension-reduced operation is made on the prefiltered data and target steering-vector as follows:As a result, the dimension-reduced space-time filter is as follows:where is a constant and is the correlation matrix of the dimension-reduced data . It is often taken place by the maximum likelihood estimation of itself.

5. Experimental Results

To start with, the simulation experiments are conducted to verify the performance of RTDPC.

In the simulations, the uniform linear array is antennas with interelement . In a CPI, pulses are transmitting at a pulse repetition frequency . The airborne platform is high at a speed of . The CNR (clutter to noise ratio) is . An (signal to noise ratio) moving target is injected into the underdetection range cell with Doppler frequency . The errors for the estimated velocity and crab angle are set to be and , respectively. The MVDR (minimum variance distortionless response) spectral is exploited to testify the ability of RTDPC. It is defined as [3], where means the inverse of a matrix, is the searching steering-vector, and is the covariance of clutter plus noise and target.

In Figure 3, a target with normalized Doppler frequency is injected near the main clutter region. Figure 3(a) shows the original clutter spectral. The clutter ridge is distributed along the diagonal in the angle-Doppler plane. Compared with the strong clutter, the target which is located on the right side of the clutter ridge is much weaker. Figures 3(b) and 3(c) show, with accurate prior knowledge, the clutter spectral after TDPC and RTDPC as clutter prefilters. Obviously, the energy of the moving target is considerably stronger than the main clutter. This result indicates that both TDPC and RTDPC can effectively suppress the clutter by forming the nulling along the clutter ridge and maintain the target constant when the prior knowledge is accurate.

Figures 46 show the clutter filtering performance of TDPC and RTDPC when the prior knowledge is not accurate. It is clearly demonstrated that when crab angle is inaccurate, the SNR will be considerably degraded by TDPC. However, by RTDPC, the SNR almost remains the same as in Figure 3(c) and the clutter is also effectively filtered out. Moreover, as the analysis in Section 3, inaccurate crab angle influences the performances of TDPC more than inaccurate velocity.

In Figures 7 and 8, the robust performance of RTDPC is verified under nonsidelooking scenario. For the variation of crab angle, the clutter ridge is distributed as a half-ellipse in the angle-Doppler plane, which implies the nonlinear relationship between the Doppler frequency and the spatial frequency. Under this condition, RTDPC can still effectively filter the clutter and maintain the SNR constant. However, as in the sidelooking scenario, TDPC considerably degrades the SNR of target when the prior knowledge is inaccurate. From the simulation results, we can conclude that although TDPC can effectively suppress the clutter when the prior knowledge is not accurate, the target is also impaired, which is detrimental to the following target detection.

Real Data. In order to verify a good many achievements in the theory of STAP, MCARM program was established by Rome Laboratory and Northrop Grumman Corporation [37]. The MCARM data used in the experiment come from acquisition 575 on flight 5 (rl050575). The operating parameters used to obtain the data are as follows: the platform velocity , the crab angle , the array element spacing , the radar operating wavelength , and the pulse repetition frequency . As in the simulated experiments, the estimated errors in TDPC and RTDPC for velocity and crab angle are also set to be and . 201–500 range bins are selected to verify the proposed algorithm. The first 11 azimuth channels and the first 32 pulses are exploited here. Two (signal to clutter ratio) moving targets with azimuth angles and Doppler frequency are separately injected near the clutter main beam at the range bins 350 and 400.

Figure 9 shows that the normalized output power (NOP) of each range bin by making use of TDPC + STM and RTDPC + STM. The average NOPs below the target of them are about 19.29 dB, 25.87 dB, and 34.93 dB, respectively. It is clearly demonstrated that after RTDPC filtering, STM can identify the two weak targets from the residual clutter. However, due to the inaccurate prior knowledge, after TDPC filtering, the clutter energy is still strong. Similarly, as shown in Figures 10 and 11, the average NOPs below the target for FA, TDPC + FA, and RTDPC + FA are about 33.48 dB, 41.43 dB, and 43.38 dB and those for EFA, TDPC + EFA, and RTDPC + EFA are about 45.33 dB, 45.86 dB, and 46.71 dB. The results infer that the STAP method followed RTDPC can achieve relatively good performance at the improvement of output SCNR and the corresponding detection performance.

6. Conclusions

The effectiveness of TDPC depends largely on the accuracy of estimated radar parameters. The estimation errors of platform velocity and crab angle in practice caused by the complicated airborne flight environment influence the clutter suppression ability of TDPC. RTDPC that adds the errors of platform velocity and crab angle into the design of filtering coefficient matrix for increasing the TDPC’’s applicability is proposed in this paper. Both the simulated and the MCARM real data are utilized to compare the clutter suppression performance of RTDPC with that of TDPC under the condition of inaccurate platform velocity and crab angle. The moving target detectability of following STAP algorithms is also enhanced after RTDPC as the filter.

Data Availability

All the data used in this work are based on simulations in MATLAB and are available for any research work.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

This work was sponsored by National Natural Science Foundation of China under Grants 61901371, 61901372, and 51974250, in part by the Natural Science Basic Research Program of Shaanxi under Grants 2020JQ-600 and 2020JQ-599, and in part by the Youth Science and Technology Nova Project in Shaanxi Province under Grant 2020KJXX-018.