Article citation information:

Ślesicki, B., Ślesicka, A., Kawalec, A. Improve the safety of air transport, especially in militarized terrain, by use of side looking airborne radar and space time adaptive processing. Scientific Journal of Silesian University of Technology. Series Transport. 2024, 123, 335-346. ISSN: 0209-3324. DOI: https://doi.org/10.20858/sjsutst.2024.123.17.

 

 

Błażej ŚLESICKI[1], Anna ŚLESICKA[2], Adam KAWALEC[3]

 

 

 

IMPROVE THE SAFETY OF AIR TRANSPORT, ESPECIALLY IN MILITARIZED TERRAIN, BY USE OF SIDE LOOKING AIRBORNE RADAR AND SPACE TIME ADAPTIVE PROCESSING

 

Summary. The paper explores the potential to enhance aviation safety, particularly in militarized regions, by outfitting aircraft with Side Looking Airborne Radar (SLAR) and employing space-time adaptive processing (STAP) algorithms. The research objective revolves around implementing a model of side-looking airborne radar and the corresponding STAP algorithms. This technology enables the detection of slow-moving targets amidst strong interference, encompassing both passive (clutter) and active (jammer) elements. Slow-moving targets relative to the aircraft's speed include tanks, combat vehicles, command vehicles, artillery, and logistical assets of enemy forces. The theoretical framework of space-time adaptive processing is presented, elucidating the sequential steps of the classical Sample Matrix Inversion Space-Time Adaptive Processing (SMI STAP) algorithm. The paper underscores the significance of characteristic parameters delineating a linear STAP processor. The proposed solution facilitates the detection of enemy combat measures and enhances aviation safety. It outlines a radar model installed beneath the aircraft's fuselage and elucidates algorithms for space-time adaptive processing of radar signals. The simulations conducted within the article were executed using the MATLAB environment. The simulation results indeed suggest that the proposed solution holds promise for deployment in equipping aircraft of one's own military and those engaged in operations within conflict zones. This paper stands as one of the few contributions in the literature addressing the augmentation of aircraft safety through radar and space-time adaptive processing.

Keywords: radar, safety, space-time adaptive processing, signal processing, airborne radar, air transport

 

 

1. INTRODUCTION

 

Today, radars can be categorized based on their function, operational mode, range, probing signal type, or signal processing methodology. Nevertheless, this merely provides an overview of the wide range of radar systems customized for diverse applications. Additionally, radars can be classified based on the platform on which they are deployed.

Deploying radar on an aircraft surface inevitably introduces several technical challenges. Firstly, the high velocity of the aircraft leads to moving received echoes, necessitating specialized algorithms for their mitigation. Secondly, achieving effective target detection amidst significant clutter poses another substantial challenge. Consequently, this study explores the implementation of Side Looking Airborne Radar (SLAR) installed beneath the aircraft fuselage, coupled with Space-Time Adaptive Processing (STAP) algorithms to counteract interfering signals.

STAP represents a contemporary signal processing approach applied in radar systems. This technique is instrumental in detecting ground-based moving targets via radar systems mounted on airborne platforms. In the literature, radar interference is categorized into passive interference (referred to as “clutter”) and active interference (known as “jammer”). Passive interference encompasses echoes from signals reflected off buildings, vegetation, and other terrain obstacles, while active interference includes radio interference transmitters hindering the operation of one's own radar systems [1].

Aircraft equipped with active radar sensors and state-of-the-art signal processing capabilities possess real-time control and surveillance capabilities over their operational airspace. Enhanced operational awareness for pilots mitigates risks and threats posed by enemy combat assets. Therefore, the proposed solution holds potential for deployment in outfitting aircraft for military forces and those engaged in operations within conflict zones.

 

 

2. A LITERATURE REVIEW

 

The history of the beginning of STAP technology dates back to the early 1970s [2].
As a result of the work done by researchers, several monographs were published summarizing the current state of knowledge on STAP [3]. The publications of the aforementioned authors laid the foundation for future STAP research [4]. Until the late 1990s, STAP algorithms were based on statistical methods for estimating the disturbance covariance matrix. A very popular method, which later often served as a reference for new methods, was the Sample Matrix Inversion (SMI) method [3].

The publication [5] marked a significant advancement in scientific inquiry regarding the evolution of methods for estimating interference covariance matrices. It introduced a novel form of non-statistical interference covariance matrix estimation termed the Direct Data Domain method, commonly abbreviated as D3, making its debut appearance in the literature.

Subsequently, research on statistical methods for estimating clutter covariance matrices was discontinued, which seems appropriate, given the drawbacks with which they were characterized. The main ones include the need to access a huge amount of training data contained in distance bins, which was not infrequently difficult to fulfil. In the following years, efforts were made to introduce improvements to D3 methods by eliminating the disadvantages that occurred, which caused phenomena such as difficulties in detecting a target against a background of inhomogeneous interference or confusing target detection [6].

Another approach to estimate the clutter covariance matrix was through the Knowledge-Aided STAP (KA-STAP) method [7]. The concept behind this technique was to derive the clutter covariance matrix using information about the radar-scanned area or the present targets [8]. It is important to highlight that during the early 21st century, one of the primary research objectives of the US Defence Advanced Research Projects Agency (DARPA) was the development of the KA-STAP method [9]. Consequently, it is noteworthy that, based on the existing literature, it can be inferred that assessing the accuracy of prior knowledge about the scanned terrain compared to the actual radar data received. This constituted the primary focus of scientific research on KA-STAP algorithms.

In recent years, a fresh avenue of research has emerged in the advancement of non-statistical techniques for clutter covariance matrix estimation [10]. These approaches rely on sparse recovery algorithms within STAP [11]. Consequently, according to reports from the literature, there has been a surge in research publications over the past few years focusing on the refinement of methods for clutter covariance matrix estimation in STAP, particularly emphasizing the utilization of sparse recovery algorithms [12].

 

 

 

Fig. 2. STAP signal processing [3]


 

In order to realize the essence of the issues discussed in the publication, it seems very helpful to present a diagram of STAP processing, on which the process of estimating the clutter covariance matrix can be clearly located [3]. Hence, the above figure highlights the most important steps in the processing of the radar signal in STAP technology.

 

 

3. RESEARCH PROBLEM

 

The problem of this research is the implementation of a model of side looking airborne radar and the corresponding algorithms for space-time adaptive signal processing.

 

 

 

Fig. 2. Geometry of radar [4]

 

 

The paper adopts a model of the localization of the radar and the target according to the diagram in Figure 2. The letter P was used to denote the target, the distance of the radar to the target was denoted Rs, the elevation angle between the radar and the target was denoted θ, and the azimuth angle between the radar and the target was denoted φ. The flight altitude of the aircraft is denoted H. The aircraft is moving in a straight line at a constant speed Va. Located under the fuselage, the radar has a uniform linear antenna (ULA) radiating the target P at an angle α. The radar emits electromagnetic waves at a wavelength λ, the radar signal is emitted as a sequence of M coherent pulses with a pulse repetition rate fr. The pulses are emitted by N antennas spaced at a fixed distance d from each other. Accordingly, the received echo forms a radar data cube. The radar data cube consists of composite signal samples collected for M pulses by N antenna array elements for distance intervals from 1 to K.

Every STAP algorithm analyses the unprocessed data, utilizing a particular snapshot of the radar cube at a predetermined distance bin k [3]. Following this, assessments regarding the presence or absence of the target within the specified distance bin are commonly conducted.
To accomplish this, a filter is devised, distinguished by its strong amplification of the pertinent signal from the target, alongside substantial suppression of all other signals (including interference from stationary objects and jammers). STAP is engineered to eliminate echoes originating from interference sources while preserving the signal originating from the target of interest [1].


 

4. RESEARCH METHOD

 

The data contained for the k-th snapshot of the raw data cube is a matrix [1]:

 

                                      .                                        (1)

 

For further processing, the matrix should be regrouped into a vector:

 

            .                         (2)

 

 denotes the received echo.   consists of the sum of signals from the object , clutter , jammer  and noise  [1]:

 

                                               .                                                      (3)

 

According to the literature, the first step in STAP processing is to determine the control vector  for each target. This vector is the Kronecker product of the temporal steering vector  and the spatial steering vector .  The steering vector in the time domain is shown as [1]:

 

                                                   (4)

 

where  is the Doppler frequency shift:

 

                                                              (5)

 

where Vr is the radial velocity between the target and the radar. The temporal steering vector is given as [1]:

 

                                               (6)

 

where  is the spatial frequency, which is expressed as:

 

.                                                            (7)

 

Finally, the steering vector takes the form [1]:

                                             (8)

 

where the symbol   denotes Kronecker product.

 

Subsequently, the following stage involves the calculation of the clutter covariance matrix. To achieve this, the distance bin under examination, denoted as k, is partitioned into  bins
of clutter. It is postulated that the clutter component within a specific distance bin results from the amalgamation of signals emanating from each clutter patch. Consequently, the covariance matrix encompassing noise, clutter, and jammer for a particular snapshot of the radar data cube is delineated as [1], respectively [13]:

 

                                                                                                    (9)

 

                                                        (10)

 

                                                                       (11)

 

where  is the power of the clutter source,  and  denote the unit matrix of dimension  and ,  is the steering vector of an individual clutter bin,  is the control vector of a given active interference source, CNR is the clutter to noise ratio measured in decibels and JNR denotes the jammer to noise ratio measured in decibels, respectively. Thus, the covariance matrix representing the combined effects of clutter, jammer, and noise can be expressed by the equation:

 

.                                                (12)

 

By establishing both the steering vector and the covariance matrix encompassing clutter, jammer, and noise, ascertain the weight vector [14]:

 

                                          (13)

 

where  is a scalar,  is the inverse covariance matrix of the sum of clutter, jammer and noise.  is the steering vector for one target to be detected. In practice, both  and   are unknown [15].

 

As previously noted, the STAP processor functions as a linear filter, with its primary objective being the elimination of clutter, jammer, and noise to facilitate target detection. The relationship governing these processes is detailed in [16]:

 

                                                         (14)

 

where  is the resulting scalar. The symbol H denotes the Hermitian transpose of the given matrix.

 

At this point it is worth quoting a very important parameter describing the degree of clutter in the received echo signal. This parameter is the clutter ridge slope β. In general, the higher its value, the more difficult it is to detect the target. Its relationship is given by the formula [17]:

 

.                                                           (15)

 

This equation defines the position of the clutter in the space-time plane. It is worth noting that the distance between antennas d takes a constant value at the stage of radar antenna design. Hence, when operating the radar at a specific pulse repetition frequency , the value of the β parameter can be directly influenced by selecting the value of the flight speed of the platform  [3].

Another crucial measure indicating the clutter level in the received signal is the rank of the clutter covariance matrix. The rank of a matrix signifies the highest count of linearly independent vectors that form its rows (or columns). The rank of a matrix is expounded upon by Brennan's equation in the referenced paper [3]. The authors in question approximated the rank of the clutter covariance matrix as expressed by [18]:

 

           .                                     (16)

 

The parentheses denote rounding to the nearest integer.

 

 

5. RESEARCH RESULTS

 

The simulations conducted in this paper were executed using the MATLAB platform. Initially, simulation analyses were undertaken to assess the scope of application and practical viability of the proposed solution. In this context, the impact of the platform's velocity on radar performance under fixed parameters, as commonly referenced in literature [19], was investigated. Furthermore, the investigations were enhanced by examining the eigendecomposition of the clutter covariance matrix across various radar and platform configurations.

The last section presents the results of a complex simulation of a heterogeneous environment in which a target was randomly placed. The task of the radar and the STAP processing used was to eliminate clutter and correctly detect the target.

The radar and platform models as specified in paragraphs 3 and 4 were adopted for the simulation, with the parameters included in the corresponding tables.


 

Tab. 1

Parameters adopted to simulation.

 

Parameter

Value

Antennas

10

Pulses

10

Radar operating frequency

8 GHz

Wavelength

0.0375 m

Distance between antennas

0.01875 m

Platform flight altitude

2000 m

Clutter ridge β

0.4…2.0

Pulse repetition frequency

12 000 Hz

Clutter to noise ratio

30 dB

Signal-to-noise ratio

10 dB

 

5.1. Scope of applicability and operational use of the developed solution

 

Figure 3 depicts a graph illustrating the clutter's location in the received signal for various velocities of the platform's flight. As previously established, the velocity of the platform's flight directly influences the β parameter. The β parameter values are presented within the range of 0.4 to 2.0. For the parameters outlined in Table 1 and a β parameter value of 1, a velocity of 90 m/s was determined, representing the optimal scenario in terms of target detection precision. Figure 3 illustrates that with an increase in the β parameter, clutter (represented by more and more lines) emerges across a wider range of azimuth angles observed by the radar.

 

5.2. Analysis of the eigendecomposition of the clutter covariance matrix for individual radar and platform parameters

 

Additionally, the computational complexity, which directly impacts the processing time and target detection efficiency, was calculated for the parameters listed in Table 1. The figure below uses vertical lines to represent the rows of the clutter covariance matrix determined directly from Brennan's equation for various values of the parameter β. For the given simulation data, the row of the clutter covariance matrix is , as shown by the red vertical line in Figure 4. It is also evident that an increase in the parameter β leads to a higher order of the clutter covariance matrix, thereby increasing the computational complexity of STAP processing and making target detection more challenging.

 

5.5. Correct target detection

 

To verify the accuracy of the proposed STAP method using the SMI algorithm, an additional simulation was conducted. A radar system consisting of 10 antennas operating at 10 GHz was assumed. The distance between the antennas was set to half the wavelength, which in this case is 0.015 m. The radar was mounted on an airborne platform moving along the axis of the antennas at a constant speed of 225 m/s at an altitude of 2000 m. For the given speed and pulse repetition frequency, the parameter β is equal to 1. The radar-cross section (RCS) of the target was set to 1 .

 

Fig. 3. Graphs of the location of clutter in the received signal for
different values of the platform's flight velocity

 

 

 

Fig. 4. Plot of the eigendecomposition of the clutter covariance matrix against the value of
the parameter
β

 

A commonly utilized model of heterogeneous forest-covered terrain, known in the MATLAB environment as the gamma model, was employed as clutter. Additionally, a jammer was positioned in the field. The parameters used for the simulation are presented in the following tables.

 

                                                                                                                                         Tab. 2

Parameters adopted to simulation

 

Parameter

Value

Velocity of target [ x, y, z]

[40 m/s, 40 m/s, 0 m/s,]

Location of target

x = 1500; y = 1500; z = 0;

Power of jammer

100 W

Location of jammer [x, y, z]

x = 1100; y = 1200; z = 0;

Noise

Gauss noise

 

As a result of the simulation, a string of 10 pulses was transmitted through the antenna array, with a pulse repetition frequency of 30 kHz and a pulse duration of 33 µs. Then the signal reflected from the target but also the clutter from the ground surface, the interference signal from the transmitter and the noise on the receiving side were received.

Figure 5a illustrates the values of signals received by the radar's antenna array as a function of range, following the transmission of the first pulse. At this point, due to the presence of significant clutter, the radar is unable to accurately determine the target's position. It is evident that the radar incorrectly suggests the target is located 1000 meters away.

Figure 5b depicts the values of signals received by the radar's antenna array as a function of range, after the first pulse transmission. This time, however, the raw data has undergone STAP processing using the proposed SMI-STAP algorithm implemented in the MATLAB environment. It is clear that the radar correctly identifies the target as being approximately 2200 meters away in a straight line.

 

Fig. 6. Space time spectrum of interference before (a) and after (b) SMI-STAP processing

6. SUMMARY AND CONCLUSION

 

The article explores the potential for improving the safety of air transport, particularly in militarized areas, by outfitting aircraft with SLAR and the STAP algorithm. It outlines the theoretical foundations of STAP processing and describes the specific steps of the classical SMI STAP processing algorithm. The theoretical analyses, calculations, and simulation results conducted during the study lead to the following conclusions:

-       in STAP processing, several approaches to estimating the clutter covariance matrix are suitable, with a shift away from statistical methods. At present, the most advanced techniques for estimating the clutter covariance matrix are non-statistical methods;

-       the influence of the flight parameters of the flying platform on which the radar is mounted, such as speed and the angle of deviation from the axis of the antennas through the parameter β, which uniquely determines the degree of clutter and the frequency of its occurrence, is emphasized;

-       the paper is one of the few items in the literature on enhancing aircraft safety with radar and STAP processing. It represents an original contribution to the development of knowledge and radar technology, and at the same time, this item, through the publication of some of the results in world journals, represents an important point on the map of the development of the field of STAP;

-       the theoretical analyses presented in the article, the calculations performed, and the simulation results obtained are of great practical importance. The use of the STAP technique in radar systems contributes to increasing their usability and reducing computational complexity requirements.

 

In summary, this paper demonstrates that it is feasible to detect slow-moving objects (such as tanks, artillery, and drones) and enhance pilots' situational awareness during flight missions, especially in air transport operations.

 

 

References

 

1.        Klemm Richard. 1998. Space-time Adaptive Processing: Principles and Applications. The Institution of Electrical Engineers. ISBN: 978-0-8529-6946-5.

2.        Reed Irving, John Mallett, Lawrence Brennan. 1974. Rapid convergence rate in adaptive arrays”. IEEE Trans. Aerosp. Electron. Syst. 10(6): 853-863. DOI: 10.1109/TAES.1974.307893.

3.        Ward James. 1994. Space-Time Adaptive Processing for Airborne Radar. Lincoln Laboratory Technical Report 1015.

4.        Guerci Joseph. 2014. Space-Time Adaptive Processing for Radar. Artech House. ISBN: 978-1-6080-7820-2.

5.        Sarkar Tapan, Srikanth Nagraja, Michael Wicks. 1998. A deterministic direct data domain approach to signal estimation utilizing non uniform and uniform 2D arrays”. Dig. Sig. Proc. 8: 114-125.

6.        Carlo Jeffrey, Tapan Sarkar, Michael Wicks. 2003. A Least Squares Multiple Constraint Direct Data Domain Approach for STAP”. In: 2003 IEEE Radar Conference: 431-438. 5-8 May 2003, Huntsville, AL, United States.

7.        Melvin Wiliam, Gregory Showman. 2006. An approach to knowledge-aided covariance estimation”. IEEE Transactions on Aerospace and Electronic Systems 42(3): 1021-1042. DOI: 10.1109/TAES.2006.248216.

8.        Zhu Xumin, Peter Stoica. 2011. Knowledge-aided space-time adaptive processing”. IEEE Trans. Aerosp. Electron. Syst. 47(2): 1325-1333. DOI: 10.1109/TAES.2011.5751261.

9.        Peng Hao, Yuze Sun, Yang Xiaopeng. 2019. Robust knowledge-aided sparse recovery STAP method for non-homogeneity clutter suppression”. The Journal of Engineering 20: 6373-6376. DOI: 10.1049/joe.2019.0273.

10.    Chen Jie, Xiaoming Huo. 2006. Theoretical results on sparse representations of multiple-measurement vectors”. IEEE Trans. on Signal Processing 54(12): 4634-4643. DOI: 10.1109/TSP.2006.881263.

11.    Duan Keqing, Zetao Wang, Wenchong Xie. 2017. Sparsity-based STAP algorithm with multiple measurement vectors via sparse Bayesian learning strategy for airborne radar”. IET Signal Processing 11(5): 544-553. DOI: https://doi.org/10.1049/iet-spr.2016.0183.

12.    Zang Wei. 2019. Reduced dimension STAP based on sparse recovery in heterogeneous clutter environments”. IEEE Trans. on Aerospace and Electronics Systems 56(1): 785-795. DOI: 10.1109/TAES.2019.2921141.

13.    Cristallini Diego, Wolfram Bürger. 2012. A robust direct data domain approach for STAP”. IEEE Trans. on Sig. Proc. 60(3):1283-1294. DOI: 10.1109/TSP.2011.2176335.

14.    Cristallini Diego. 2012. Exploiting robust direct data domain STAP for GMTI in very high resolution SAR”. In: 2012 IEEE Radar Conference: 0348-0353. 7-11 May 2012, Atlanta, GA, United States.

15.    Guo Yiduo, Guisheng Liao, Weike Feng. 2017. Sparse representation-based algorithm for airborne radar in beam-space post-Doppler reduced-dimension space-time adaptive processing”. IEEE Access 5: 5896-5903. DOI: 10.1109/ACCESS.2017.2689325.

16.    Li Ming, Guohao Sun, Zishu He. 2019. Direct Data Domain STAP Based on Atomic Norm Minimization”. In: 2019 IEEE Radar Conference: 1-6. 22-26 April 2019, Boston, Massachusetts, United States.

17.    Ma Zeqiang, Yimin  Liu, Huadong Meng. 2013. Jointly sparse recovery of multiple snapshots in STAP”. In: 2013 IEEE Radar Conference. 1-4. 29 April - 3 May 2013, Ottawa, Ontario, Canada.

18.    Satyabrata Sen. 2015. Low-rank matrix decomposition and spatio-temporal sparse recovery for STAP radar”. IEEE Journal of Selected Topics in Signal Processing 9(8): 1510-1523. DOI: 10.1109/JSTSP.2015.2464187.

19.    Knee Peter. 2012. Sparse representations for Radar with MATLAB. Examples. Morgan & Claypool. ISBN: 978-1-6270-5034-0.

 

 

Received 05.12.2024; accepted in revised form 10.03.2024

 

 

 

Scientific Journal of Silesian University of Technology. Series Transport is licensed under a Creative Commons Attribution 4.0 International License



[1] Faculty of Aviation Division, Polish Air Force University, Dywizjonu 303 no 35 Street, 08-520 Dęblin, Poland. Email: b.slesicki@law.mil.pl. ORCID: https://orcid.org/ 0000-0002-0857-1081

[2] Institute of Navigation, Polish Air Force University, Dywizjonu 303 no 35 Street, 08-520 Dęblin, Poland. Email: a.slesicka@law.mil.pl. ORCID: https://orcid.org/ 0000-0002-6313-030X

[3] Faculty of Mechatronics, Armament and Aerospace, Military University of Technology, gen. Sylwestra Kaliskiego 2 Street, 00 -908 Warszawa, Poland. Email: adam.kawalec@wat.edu.pl. ORCID: https://orcid.org/0000-0003-3930-7504