Deterministic Aided STAP for Target Detection in Heterogeneous Situations

Classical space-time adaptive processing (STAP) detectors are strongly limited when facing highly heterogeneous environments. Indeed, in this case, representative target free data are no longer available. Single dataset algorithms, such as the MLED algorithm, have proved their efficiency in overcoming this problem by only working on primary data. These methods are based on the APES algorithm which removes the useful signal from the covariance matrix. However, a small part of the clutter signal is also removed from the covariancematrix in this operation. Consequently, a degradation of clutter rejection performance is observed.We propose two algorithms that use deterministic aided STAP to overcome this issue of the single dataset APESmethod.The results on realistic simulated data and real data show that these methods outperform traditional single dataset methods in detection and in clutter rejection.


Introduction
In the context of radar signal processing, the purpose of space-time adaptive processing (STAP) is to remove ground clutter returns, in order to enhance slow moving target detection.STAP performs two-dimensional space and time adaptive filtering where different space channels are combined at different times [1].Filter's weights are adaptively computed from training data in the neighborhood of the range cell of interest, called cell under test (CUT).The estimation of these weights is always deducted, more or less directly, from an estimation of the covariance matrix of the received signal, which is the key quantity in the process of adaptation [2].Any implementation of STAP processing must remain absolutely consistent with the strategy of radar processing whose purpose is to obtain a high probability of detection while keeping a very low probability of false alarm.
Classical space-time adaptive processing (STAP) detectors are strongly limited when facing a severe nonstationary environment such as heterogeneous clutter.Indeed, in these cases, representative training data are no longer available.The Maximum Likelihood Estimation Detector (MLED) [3] is a single dataset detector among others [4].It only operates with the data from the cell under test, hence its performance is not impacted by nonstationarity.Of course, no environment is purely heterogeneous or homogeneous and the problem can be addressed by combining primary and secondary data [5].We will here consider the environment to be heterogeneous enough to only use primary data.To make the primary data target-free, the MLED detector removes a thin part of the signal of the Doppler cell under test from the covariance matrix.A slight part of the clutter is removed along the target signal which implies a degradation of clutter rejection, especially if the number of Doppler cells is low.The less Doppler cells, the more the clutter removed from the covariance matrix and the worse the estimation of the covariance matrix.The bad estimation of the matrix can be addressed by using subspace methods [6] but the removal of some clutter is inherent to the APES method.
In this paper, we will show how we can overcome this problem by the use of deterministic aided STAP.Moreover, we will extend this method to the Stop-Band APES which greatly reduces the computational workload of the MLED detector.
Section 2 is devoted to the data model, and Section 3 summarizes the principle of the MLED APES-based detector and the Stop-Band APES algorithm.A deterministic based nonadaptive approach of space-time processing is presented in Section 4. In Section 5, we describe two different approaches for deterministic aided STAP and finally, in Section 6, simulations are given to show that the proposed methods outperform the MLED and Stop-Band algorithms.

Data Model
Consider a radar antenna made of  sensors that acquires   pulse snapshots for each range gate .We will only use the primary data so we will forget the range gate dimension, also called fast-time dimension.Then, the processing algorithm works independently in each range cell.We adopt the following two hypothesis models where  0 and  1 mean that no target or a target is present, respectively: where the received data have been arranged into an  ×   matrix X with   being the number of training data pulse snapshots,  being the number of pulses of the spatio-temporal vector, and  being the complex amplitude of the target.s  is the spatiotemporal steering vector (length ), s  is the temporal steering vector (length   =   −  + 1), and N is the interference (clutter plus noise) matrix.
We make use of a temporal sliding window to work on the temporal dimension; consequently, the estimated covariance matrix R is obtained from X as follows: One classical STAP detector taken as reference uses the Adaptive Matched Filter (AMF) [1,2].The filter w is Detection is achieved by comparing the output SNIR power of the matched filter to a threshold as follows: In case where a strong target is present at this range gate, R contains the target covariance matrix.Consequently, the target is removed with the clutter and it can no longer be detected by ( 4).This happens when many targets are moving at the same speed but are at different distances (roads, highways, convoys, etc.).Another problem with this detector is that the ground clutter has to be homogeneous on the range domain.Otherwise, the clutter used to estimate the covariance matrix will not be representative of the clutter that has to be canceled, leading to a bad clutter rejection.

APES-Based STAP Detectors
3.1.The Maximum Likelihood Estimation Detector.To overcome the previous issues of signal suppression or the none representativeness of secondary data, the MLED detector [7] based on the APES [8] algorithm removes the signal of interest from the covariance matrix.The problem is stated as follows: The obtained solution is where Detection is achieved using the output power normalized by the Adaptive Power Residue (APR = w   = s   Q −1 s  ) as follows: To avoid strong signal loss due to covariance matrix estimation errors [9], one may use in addition diagonal loading [10] or subspace methods [11].

Extension to
where P // is the projector into the target signal subspace: The problem (5) can then be recognized as a minimization of the interference-plus-noise energy outside the subspace spanned by the target as follows: The solution is still w = (Q (−1) s  )/(s   Q (−1) s  ) but with the more general form for Q: This latest formulation not only shows the hyperresolution property along the frequency domain but also allows overcoming one major drawback of the MLED method for our application.The MLED has indeed a high-frequency resolution due to the sharpness of the projection I − P // with P // = (s * t s T t )/(s   s *  ) (dash curve, Figure 1).This is a problem because it requires a strong oversampling to be sure to remove the signal of interest from the covariance matrix and so it leads to an important increase of the computing load.In order to avoid this problem, we propose a new detector called Stop-Band APES.The minimization is using a projector P // on an extended subspace around the Doppler frequency  0 under test.For instance, two adjacent half-cells can be added into the space spanned by P // = S *  (S   S *  ) −1 S   with The sharpness and effectiveness of the cancellation around the target signal are characterized by the frequency response of the projector, which is, for a signal X at frequency  (X = s T t ()), as follows: Figure 1 shows that building a projector with two adjacent half-cells is enough to correctly remove the signal in the cell under test.Nevertheless, compared to the MLED, the Stop-Band APES does not require oversampling of the Doppler resolution for the calculation and the application of the STAP filter.A zero-padding by a factor of 2 will still be required to access the signal that has to be evaluated every half-resolution cells for the creation of the projector [12,13].

Limitations of the MLED and Stop-Band APES.
In order to explore the use of subspace-based methods, we have to go deeper in the formulation of the MLED detector.Indeed, these methods will only work if the clutter subspace of the covariance matrix R remains very close to the clutter subspace of the target-free covariance matrix Q.For a given distance cell, if there is no target at this range, the covariance matrix R only contains interference, that is, clutter and possibly jamming signal and noise, according to (2) as follows: We can demonstrate [11] that the matrix Q is, without approximation: The matrix (NN H )/  is the interference-plus-noise estimated covariance matrix, whereas (Ns *  s   N  )/( 2  ) is the scalar product of interference-plus-noise vectors with their projection on s *  .It follows from ( 16) that the modified covariance matrix Q used for MLED in (8) does no longer contain the target contribution and that the target will not be removed contrarily to the clutter by the MLED STAP filter (6).
The residual clutter-plus-noise covariance matrix is slightly different from the actual covariance matrix (N  )/   (Figure 2).The term (Ns *  s T t N H )/( 2  ) represents the part of the clutter that is removed from the covariance matrix.The number of Doppler cells being usually high, the projector is consequently very sharp; that, the term (Ns *  s T t N H )/( 2  ) is small and both MLED and Stop-Band APES, which removes a wider part of the clutter from the covariance matrix (Ns *  s T t N H )/( 2  ), are all working.This effect can be seen on Figure 2 in a situation with and without target in the Doppler cell tested.
However, in a situation where the number of Doppler cells is low, we will observe a degradation of the clutter rejection performance of the MLED detector, and this degradation will International Journal of Antennas and Propagation be even worse for the Stop-Band APES algorithm.This effect is due to the partitioning which is done only in time domain.If spatio-temporal partitioning is employed, only a single bin of the angle-Doppler plane is removed but the computational cost would hugely increase because of the angle-Doppler scanning.We will present in the next section a deterministic processing and, in Section 5, a new method that makes use of deterministic processing to solve this problem.

Deterministic Space-Time Processing
We will here briefly describe a nonadaptive space-time processing which is the basis of the deterministic aided STAP processing we will introduce in the following section.For a side-mounted antenna, the clutter occupies a one-dimension position in the two-dimensional Doppler-angle domain.The clutter Doppler frequency is a function of the receiving angle as follows: with  being the Doppler frequency of the clutter,  being the receiving angle,  being the platform speed, and  being the wavelength of the radar frequency.Knowing this relation, we can build a filter that will remove all the signal that is in the 1D-domain driven by (17).The general form of the filter, which will be referred in the following to non-adaptive or deterministic processing, has the same form as AMF in (3) as follows: but with where Γ N is the true noise covariance matrix (identity matrix in our case),  is the number of mainlobe clutter patches, and s  (  (  ),   ) is the space-time steering vector of angle   and frequency   obtained with (17).In the same formulation of the filter as MLED and Stop-Band APES in (6), the matrix for each Doppler cell can be written as follows: where the vector s  (  (  ),   ) is the predicted steering vector of the clutter.In this case, to process one Doppler cell, the steering vector s c ( 0 ( 0 ),  0 ) of the Doppler cell under test and the two steering vectors s c ( ±1 ( ±1 ),  ±1 ) from the adjacent Doppler cells are sufficient to correctly remove the clutter.However, the performance of this nonadaptive approach is very limited in practical situations because of the heterogeneity of the clutter (e.g urban or mountainous areas) and because of antenna/receivers calibration errors which make the real steering vector of the antenna slightly different from the actual steering vector used to build the covariance matrix.To illustrate this effect, we compare the non-adaptive processing (18) to the classic adaptive processing on two sets of data.The first data are the simulated data that we build using the true spatial steering vector and the second set of data is made of realistic data simulated by a STAP simulator that emulates phases errors on the receiving channels and randomly adds impulsive echoes in the clutter.In both cases, clutter is Gaussian, homogeneous, and set to 40 dB.No target is present in these data.A side-looking antenna with four uniformly spaced subarrays is used.Aircraft speed is set to 100 ms −1 , radar frequency is 10 GHz, and the pulse frequency (PRF) is 2 kHz.The non-adaptive processing is only applied in the positive speed domain, that is to say that the negative speeds show the sum channel.The adaptive processing is applied on all the Doppler (speed) domains.
As we can see from Figures 3 and 4, non-adaptive processing works well on the error-free simulated data.The clutterto-noise ratio (CNR) is close to 0 dB, like in the adaptive processing as shown on Figure 5, which implies an attenuation of 40 dB.On the realistic simulated data, the non-adaptive processing is not performing well, as it fails to suppress the clutter.Indeed, as we can see on Figure 6, the residual CNR is near 15 dB in the main lobe; the clutter attenuation is limited to 25 dB, implying many false alarms.The full range-Doppler maps also point out this effect in Figures 7 and 8. From these results, we deduce that we cannot use a non-adaptive space-time processing in real situations but we may use the deterministic of the clutter Doppler-angle relation together with adaptive processing to achieve better performance.

Deterministic Aided GMTI STAP.
In GMTI operation, there are two main concerns about heterogeneous environments: clutter heterogeneity (land relief, urban environments) and high-density target area (roads, highways. ..).In many cases, few training data are available and the use of single data set methods is a very helpful alternative (see Section 3.1 and Section 3.2).To overcome the problem of these methods pointed out in Section 3.3, we propose a new method that includes some aspects of the non-adaptive processing.We saw in (16) that the term gg H in (7) removes the interest signal (if any) and also a small part of the clutter.
The idea here is to try to readd this clutter into the covariance matrix.The covariance matrix is then as follows: where g c is the projection of g on the clutter steering vector s c ((), ) as follows: with We can demonstrate that the covariance matrix Q of ( 16) can now be written as follows: If the clutter follows the theoretical Doppler-angle relation of ( 17), then the projection of the signal on the angle-Doppler steering vector will be close to the clutter signal that has been removed from the matrix (P c X ≈ ), and the covariance matrix T will be close to Note that we do not need to set an arbitrary clutter power value because the energy of the clutter is included in g c (cf. ( 22)).
In the case of Stop-Band APES, where the signal notch is wider, we use an extended projector P c as follows: with

Deterministic Aided STAP Processing for Air-to-Air Mode.
In air-to-air situations, the problem is different.The spectral occupation of the mainlobe clutter is much smaller than that of GMTI, whereas clutter sidelobes are much more powerful and have to be cancelled.Moreover, target density is very low, compared to GMTI.As we do not have access to a Dopplerangle relation of the mainlobe clutter, we propose another approach to readd this clutter which is partially removed in the APES-based methods.In air-to-air mode, the mainlobe clutter is pretty homogeneous in the range domain.
We will exploit this property to estimate the matrix g c g H c on the range gates domain.For each Doppler cell, the covariance matrix T is defined by However, the matrix C which was equal to g c g H c in ( 21) is now estimated as follows for each Doppler cells: with g  being the vector g = Xs *  /  of the range cell  and  being the total number of range cells.If the clutter is homogeneous, then we can make the following approximation: This assumption implies that only the homogeneous component of the clutter will be readded in the covariance matrix.
The density of the target has to remain low, otherwise useful signal will be nonnegligible in the matrix and SNR of targets will be attenuated.In the case of Stop-Band APES, ( 12) becomes with X the data of the range cell under test, X  the data of the range cell , and  the total number of range cells.

Results
6.1.GMTI Simulations.We test the GMTI deterministic aided STAP described in Section 5.1 on real airborne data.These data were obtained using the ONERA RAMSES radar system [14], which is a 4-channel ULA antenna.The aircraft speed is   = 85 m ⋅ s −1 , pulse repetition frequency is PRF = 1.5625 kHz, the number of range gates is 300, the number of time taps used to form the space-time data is  = 6, and the total number of time snapshots (radar pulses) is 64.Three targets are present in the scene (see Table 1).The Doppler-range of the sum channel (Figure 9) clearly emphasizes the heterogeneous clutter.The next figures present the results for the classical STAP (estimation on 10 range gates with 2 guard cells), MLED STAP, Stop-Band STAP, and deterministic aided Stop-Band STAP (estimation on 3 range gates with no guard cells for all processing).No oversampling is used for the STAP processors (although a 2x zero-padding is needed to access the data of the half-resolution Doppler cells in the case of Stop-Band) except for the MLED detector, which uses a 4x-oversampling.
The classical STAP processing fails to correctly remove the heterogeneous clutter (Figure 10).The MLED STAP whose signal notch is very sharp also fails to completely remove the clutter.Due to its property of high resolution, the target Doppler extent is very thin and it is difficult to    distinguish the targets on the range-Doppler map (Figure 11).As predicted, the Stop-Bland STAP processing allows even more clutter to be present as shown in Figure 12, whereas the deterministic aided Stop-Band (DA-Stop-Band) effectively cancels the clutter (Figure 13).This is done without any attenuation on target 1 which lies in the clutter.Figure 14 points out the increased clutter attenuation of DA-Stop-Band over classical Stop-Band for range gate number 149, where target 3 is present.Figure 15 shows the superiority of DA-Stop-Band over classical STAP in clutter rejection for range gate 279, an area where the clutter is particularly powerful.

Air-to-Air Simulations.
The air-to-air deterministic aided STAP (see Section 5.2) is tested on realistic synthetic data simulating an air-to-air MTI scenario.A front-looking AMSARlike antenna [15] is used for the simulations.The aircraft speed is   = 300 m ⋅ s −1 , pulse repetition frequency is PRF = 20 kHz, and the number of range gates is 100, corresponding to a physical range of 52.5 km to 59.5 km.The number of time  taps used to form of space-time data is  taps = 8 and the total number of time snapshot (radar pulses) is 128.Five targets are present in the scene (see Table 2).The sum channel (Figure 16) clearly shows that the mainlobe clutter (speeds from 230 to 280 m/s) and the sidelobes clutter occupy a wide part of the range-Doppler map.Only two targets on the left-upper part of the map are detectable without STAP processing.On Figure 17, we can see that the classical STAP processing successfully removes the homogeneous main lobe clutter and does not removes the heterogeneous sidelobes clutter.Classical Stop-Band processing cancels almost all the sidelobes clutter but does not suppress   all the mainlobe clutter (see Figure 18), whereas DA-Stop-Band (Figure 19) totally removes it.
On Figure 20, the effect on clutter attenuation of the DA-Stop-Band is visible through a comparison with classical Stop-Band.We can also observe that both types of Stop-Band processings do not completely remove the sidelobes clutter; this issue can be overcome by using subspace-based algorithms instead of matrix inversion [11].

Conclusion
In this paper, we propose two deterministic aided algorithms both based on the APES method.The first algorithm which relies on the deterministic Doppler-angle relation of the clutter is particularly adapted for GMTI detectors.The results on real data show that it outperforms both classical STAP and APES-based algorithms.The second algorithm, which aims to remove the continuous component of the interference, is on the other hand well adapted to air-to-air modes.In this case, the continuous interference is the main lobe clutter.On realistic simulated data, it totally cancels the mainlobe clutter, whereas classical STAP and traditional APES-based algorithms fail, causing many false alarms.

Figure 1 :
Figure 1: Spectral response of regular MLED S = s  (dash curve) and Stop-Band APES with two adjacent half-cells (solid curve).

Figure 2 :
Figure 2: Angle-Doppler map showing the effect of MLED projector for two different Doppler bins.

Figure 3 :
Figure 3: Comparison between sum channel (bold curve, negative speeds), deterministic (bold curve, positive speeds), and adaptive space-time processing (dash curve) on error-free simulated data.

Figure 4 :Figure 5 :
Figure 4: Range-Doppler map of the nonadaptive processing (positive speeds) and the sum channel (negative speeds) on error-free data simulated data.

InternationalFigure 8 :
Figure 8: Range-Doppler map of the adaptive processing on realistic data.

Figure 9 :Figure 10 :
Figure 9: Range-Doppler map showing the sum channel of the RAMSES data.

Figure 11 :
Figure 11: Range-Doppler map on RAMSES data showing the performance of the MLED detector.

Figure 12 :
Figure 12: Range-Doppler map on RAMSES data showing the performance of the Stop-Band detector.
Figure 13: Range-Doppler map on RAMSES data showing the performance of the deterministic aidedc Stop Band detector.

Figure 14 :
Figure 14: Comparison of classical STAP (dash curve) Stop-Band STAP (dot curve) and deterministic aided Stop-Band (solid curve) on RAMSES data for range gate number 149.

Figure 15 :
Figure 15: Comparison of classical STAP (dash curve) Stop-Band STAP (dot curve) and deterministic aided Stop-Band (solid curve) on RAMSES data for range gate number 279.

Figure 16 :
Figure 16: Range-Doppler map of the air-air realistic data showing the sum channel.

Figure 17 :
Figure 17: Range-Doppler map of the air-air realistic data showing the performance of MLED processing.

Figure 18 :
Figure 18: Range-Doppler map of the air-air simulated data showing the performance of Stop-Band STAP processing.
Figure 19: Range-Doppler map of the air-air simulated data showing the performance of deterministic aided Stop-Band STAP processing.

Figure 20 :
Figure 20: Comparison of Stop-Band STAP (dot curve) and Deterministic-Aided Stop-Band (solid curve) on air-air data at a distance of 58.5 km.