1. Introduction
Unfazed by the inherent constraints of light and barometric altitude sensors, radars employ EM waves as an optimal means to attain precise altitude information during the autonomous landing stage of UASs [
1]. Amongst the radar waveforms capable of measuring altitude, pulse Doppler radars, with their low-duty cycle, are limited by higher peak power demand, raising concerns from a public health standpoint [
2]. Additionally, the inherent rapid switching requirement impedes the measurement of the very low altitudes necessary in the touchdown stage. Conversely, the simultaneous transmission and reception in FMCW radars make them more suitable for a dynamic range of altitude requirements [
3]. Over the past three decades, mmWave automotive radars have advanced significantly, incorporating compact, highly flexible, and cost-effective 45 nm chipsets that integrate control, signal processing, and RFFE. The maturation of these platforms, coupled with researchers’ yearning for exploration, has expanded potential use cases beyond the automotive realm. The current landscape of widespread 5G communication network deployment has led to interference with universally allocated bands for RAs [
4]. Aptly, mmWave bands, particularly 77 GHz, hold considerable potential as a feasible avenue for migrating RAs in commercial aviation and a broad range of UASs.
This article is a continuation of ongoing work with an initial investigation covering theoretical and mathematical discourse concerning the potential use of mmWave automotive radars for UAS altimetry [
5]. In commercial aviation, MOPS serve as a reference for RA development [
6,
7]. The focal point of the preceding article was to align the MOPS for RAs in commercial aviation with the operational requirements of UASs. In the absence of dedicated MOPS for UASs, an adaptation from existing standards for commercial aviation was a reasonable compromise given the largely similar requirements, particularly in the landing stage. This required a systematic approach, starting with the fundamentals of FMCW radars and then delving into the complexities of performance metrics due to their interconnectedness. The rationale for choosing an automotive radar was based on commercial availability and cost-effectiveness in the mmWave band. This selection was further motivated by the promising SWaP characteristics. Moreover, the maturity of said technology was appraised with a state-of-the-art review. A compelling case for the exploration of automotive radars beyond their conventional domain was supported by comparative works. Lastly, detailed explanations for deriving waveform specifications from the operational requirements were furnished using a realistic test case.
The process of altitude measurement in UASs using FMCW radars can be categorized into two main domains. The first involves the generation of the waveform by the radar sensor in alignment with the operational requirements of the UAS. The second concerns the signal processing stage, where the waveform is converted into usable altitude data. Given the scarcity of reference literature and the novel nature of mmWave altimetry, a comprehensive and simplified discussion was essential to establish both feasibility and a foundation for this new direction. With this in mind, a detailed discussion covering state-of-the-art developments, the adaptation of MOPS for UASs, and the theoretical framework for performance metrics was provided in the previous article. Therefore, it was deemed necessary to dedicate this article specifically to the signal processing domain.
This study provides a seamless extension of the preceding work, with the framework for waveform design following accordingly. While the MOPS for commercial aviation are a good reference, there are no restrictions to stretching the performance metrics even further. These improvements are poised to serve a wide range of existing and future requirements of UASs. Furthermore, there is tremendous potential for capability enhancement. Fittingly, this work encompasses three unique stages of UAS flight: cruise, landing approach, and touchdown. The ensuing sections rationalize the unique operational requirements for each stage followed by derivation of resultant waveform specifications from a theoretical and mathematical standpoint. The premise of defining these stages is to exhibit the versatility of an adaptive waveform designed to cater to a dynamic range of operational requirements. The cruise stage is inferred from the previous work, with the addition of signal processing aspects. The term “cruise” refers to the normal flight operation of a drone while it is carrying out its intended task. This stage does not entail landing and the goal is to maintain constant altitude AGL in surveillance applications [
8]. Appropriately, high altitude estimation coupled with the finest possible range resolution is targeted for this stage.
RAs are most commonly utilized in aviation during the landing approach [
9]. In the existing landscape of industrial expansion, there is a growing trend toward the deployment of small-sized UASs in a plethora of applications, including agriculture [
10], law enforcement [
11], and traffic management [
12]. At the same time, the safety of these UASs is equally critical, particularly during the autonomous landing. Conventional RAs having a single Tx and Rx antenna pair with a wide HPBW lack the spatial information necessary for estimating true altitude along the radar boresight directly under the airframe. For instance, a highly reflective off-boresight object may cause an inaccurate altitude measurement [
13]. In light of this limitation, it is proposed that the AoA may be leveraged for the estimation of true altitude. It is envisaged that the addition of this capability will augment situational awareness significantly and improve the overall safety standards for autonomous landing of UASs. To this end, a landing scenario comprising a VTOL drone on a ship in a water body is simulated. The aim is to augment other onboard sensors during the landing stage.
WoW systems in aircraft are used to estimate the instant when its weight starts to rest on its wheels. In the context of fixed-wing commercial UASs, such a system informs the FCC that the aircraft has safely landed, and necessary control processes are engaged accordingly [
14]. By providing mm-accuracy altitude with a high level of confidence, the RA can serve as a reliable backup to the WoW system in UASs.
Primarily, considering the operational advantages of enhanced range accuracy and precise altitude estimation, the primary objective of this research endeavor is to develop an adaptive waveform. Consequently, it is imperative to elucidate in sufficient detail the signal-processing aspects related to the estimation of high altitude, AoA, and high accuracy during the three stages of flight, all within the constraints imposed by the hardware of automotive radars.
The software-defined architecture of automotive radar chosen for this endeavor allows for the realization of the intended objective. However, rationalization of operational requirements is mandatory to avoid setting too ambitious benchmarks. The solution lies in prioritizing the relevant performance metric for the respective flight stage while balancing the other metrics with reasonable trade-offs. For instance, the cruise stage does not necessitate a quick update rate and fine range resolution, but it requires a higher altitude. Similarly, during the landing stage, update rate and range resolution take precedence over high altitude. Lastly, the touchdown stage necessitates the highest possible update rate and range resolution with a very lenient altitude requirement. From an SWaP standpoint, a single radar capable of adaptively switching waveforms aligns well with the objective of improving operational endurance on a wide range of UASs.
To achieve the said objectives, legacy and contemporary signal processing techniques in the literature were reviewed. Each stage of the flight was treated as a radar scenario and simulated using the radar toolbox in MATLAB [
15]. The simulation environment provided plentiful options for generating FMCW waveforms as well as propagation configurations using phased array antennae. Typically, the task of propagating FMCW waveforms in the radar LoS and acquiring signals through ADC is performed in hardware. Subsequently, ADC samples are transferred to simulation software for post-processing. Administratively, this approach is inefficient, as emulating multiple scenarios for each flight stage consumes significant time, energy, and resources. To address this, the study presented in [
16] integrated the entire process into the simulation software, utilizing the hardware specifications of the commercial radar platform, IWR1843, operating within the 77–81 GHz frequency band. This method optimized the waveform design process while significantly conserving resources. Afterward, range profiles were subject to CFAR detection for altitude reporting. However, the touchdown stage, employing the ZFFT algorithm, was implemented directly on the radar due to the relative ease of emulating the target, the runway surface. The radar was mounted on a tripod stand and positioned against a flat concrete wall to simulate the runway surface at close range.
TDM-MIMO radars have been prevalent in automotive applications for a while now [
17]. They utilize the concept of virtual antennas to estimate the AoA of the target in the radar FOV. Conventionally, fast-time and slow-time matrices are used for detecting targets in range and Doppler bins. The use of multiple virtual antennas extends these matrices in the third dimension to form radar cubes [
18]. Signal processing techniques are then applied to these radar cubes to simultaneously estimate the range, velocity, and AoA. Within the scope of the proposed adaptive mmWave FMCW waveform for UAS altimetry, all three aspects are important. Most conventional RAs provide range information only [
19]. Academic discussions on the use of radar cubes in radar applications are largely uncharted in the literature.
The radial velocity resulting from the ROD and variations in terrain profile present a substantial difficulty that must be carefully considered during the waveform design process. Fittingly, explanations in theory and mathematics are provided for characterization and subsequent compensation. Until now, a detailed discussion on this topic has been lacking. The following sections demonstrate that, in the absence of a direct requirement to estimate the ROD, the resultant radial velocity component does not influence altitude measurement despite variations in terrain profile. Moreover, when properly accounted for, it does not impact the AoA, even with the application of TDM-MIMO.
Three widely utilized variants of CFAR employ specific parameters to balance sensitivity and the probability of false alarms [
20]. Hence, it was imperative to discuss the selection of an appropriate CFAR variant tailored to the requirements of each flight stage. This study seeks to address this gap by evaluating popular CFAR variants to optimize detection performance for ground surfaces.
The use of mmWave automotive radars has been reported for mm-level accuracy in liquid-level sensing [
21]. Traditionally, ZFFT has been used for the conservation of resources by targeting a specific portion of interest from the available spectrum. Legacy approaches aim to achieve the same spectral resolution while significantly reducing the order of FFT. Although used in mm-level accurate sensing of liquids, ZFFT for UAS altimetry has not been investigated. Contrary to the conventional approach, the objective is to target a small portion of the spectrum and perform a high-order FFT for increased spectral resolution. Consequently, very high accuracy can be leveraged for the estimation of touchdown instant as a software redundancy in the WoW system of UASs. In this work, the specifics of ZFFT are appraised for potential use in the intended application. A description of the improvised algorithm with a methodology for efficient implementation is furnished, considering the constraints of radar hardware. A coarse altitude estimation method is provided for identifying the ROI. The mathematical expressions and theoretical aspects of cited concepts may seem daunting to readers new to this field or from other backgrounds. To address this concern and in continuation with the tone of the previous article, a tutorial theme is adopted. The entire discussion is broken down into fundamental and relatable concepts.
Thus, the main contributions of this body of work are the following:
Signal processing flow in a radar cube for mmWave altimetry.
Design of an adaptive FMCW waveform for varied operational requirements.
Appraisal and compensation of radial velocity due to UAS ROD and terrain profile.
Use of TDM-MIMO for situational awareness in autonomous landing of UASs.
mmWave altimetry as a software redundancy for WoW systems using ZFFT.
Simulation of radar scenarios using hardware specifications of a commercial radar platform, IWR1843, operating at 77 GHz.
Characterization of CFAR variants for optimal detection performance.
The remainder of the article is structured as follows:
Section 2 provides an overview of the reference literature.
Section 3 summarizes the fundamental concepts and aspects common to the entirety of the article. To seamlessly bridge the discussion, an overview of methodology is provided in
Section 4. The ensuing
Section 5 covers the specifics of the cruise stage.
Section 6 entails AoA estimation using TDM-MIMO for enhanced situational awareness.
Section 7 summarizes the details of ZFFT for use in WoW systems during the touchdown stage. The authors’ discussion, future research direction, and challenges are penned in
Section 8. Finally,
Section 9 concludes the article.
5. Cruise
The cruise stage for an airborne platform is defined as the level flight segment after arrival at the initial cruise altitude until the start of the descent to the destination [
44]. In the context of UASs, this is the stage of flight where a designated task is performed. For applications requiring the UAS to maintain a constant altitude AGL, RAs are mandatory. As already explained in the preceding text, the main purpose of an RA is the real-time estimation of altitude AGL. The authors’ preliminary study orchestrated a realistic test case to maximize the range with the finest possible resolution. It is pertinent to note that although velocity estimation is not an inherent objective, the variations in aircraft altitude and terrain profile introduce a radial velocity component [
45]. It is crucial to characterize and remediate the impact of this component on altitude estimation. To this end, terrain models were adopted from MATLAB [
46].
Figure 5 shows the standard deviation in terrain elevation for various land types.
Due to the unavailability of step size, it was considered appropriate to consider it equal to the lateral motion covered by the drone in one second. Accordingly, for a lateral velocity of 20 m/s, the terrain exhibits an altitude variation of 10 m. Consequently, the radial velocity component equates to 10 m/s. This represents a stringent requirement aimed at catering to a worst-case scenario. The operational requirements for the cruise phase are listed in
Table 2, with
Table 3 containing the resultant waveform specifications.
These limitations laid the foundation for waveform design in the preliminary study as well. It was argued in the proceeding discussion that these limitations affect the performance of RAs from a signal-processing standpoint as well. It was exhibited in previous work that the choice of a longer chirp duration aids in reducing noise BW while simultaneously allowing for an IF filter with relatively narrower BW to cater to a high-altitude requirement. It is made possible by using a gradual chirp slope to achieve the same chirp BW with a longer chirp duration. This allows for a higher value of maximum range without compromising the range resolution. However, this upper bounds the
to a very small value, leading to Doppler folding with target velocity aliasing into neighboring Doppler bins. The scenarios were simulated using the Radar Toolbox in MATLAB (see
Figure 6) and it was observed that for a radial velocity component of 10 m/s, the range profile predominantly appeared in the 0th Doppler bin. This approach of empirical validation aligns well with the scope of the application, since 10 m/s represents a worst-case scenario, since actual radial velocities are expected to be significantly lower. Consequently, the 0th Doppler bin was selected as the range profile for the subsequent stage of CFAR detection, ensuring consistent detection performance.
Nevertheless, a long chirp duration has its own inherent limitations. While the range profile resides in the 0th Doppler bin, the radial velocity component causes the target to be distributed across multiple range bins. This phenomenon is illustrated in
Figure 7 for a single target, where the peak is spread over multiple neighboring bins. For a velocity of 10 m/s, the target moves 1 cm over a chirp duration of 1 ms. Point cloud data are exported from the radar at the end of each FMCW chirp frame, governing the update rate. Consequently, over the span of 16 chirps with a total duration of 16 ms, the cumulative distance covered is 16 cm. Given a range resolution of 10 cm, the target will spend a duration of 10 ms in the initial range bin and 6 ms in the neighboring bin either to the right or left, depending on the direction of relative radial motion. The only drawback of a very long chirp duration in the cruise stage is, therefore, measurement accuracy being compromised by a maximum value of 10 cm. Considering a range accuracy requirement of 45 cm, as specified in
Table 2, this trade-off is acceptable within the scope of the hardware limitations since the range resolution only degrades by a maximum cumulative value of 20 cm.
5.1. CFAR for Detection of Ground Surface
Once the range profile is acquired, the next step is to determine whether the peaks in the FFT range profile correspond to the ground surface. The detection technique employed for this study is the CFAR algorithm. It maintains a constant false alarm rate by dynamically adjusting the detection threshold based on the noise level. This approach helps in the identification of ground surface in the presence of clutter due to specular reflection. This is particularly important for the intended application since the target is the surface of the ground with a varying terrain elevation profile.
There are three major variants of the CFAR algorithm based on how the noise threshold is calculated: CFAR-CA, CFAR-CASO, and CFAR-CAGO. CFAR-CA uses the average noise level from both sides of the target cell to evaluate the detection threshold, while CFAR-CASO and CFAR-CAGO modify this approach by using the smallest and greatest average noise levels from the training cells, respectively, from either side of the target cell. Guard cells, on the other hand, are the ones adjacent to the target cell that help prevent signal leakage from the target cell into the noise estimation process. The guard cells are deliberately excluded from the noise estimation to ensure that the presence of the target signal does not bias the noise estimation, leading to a more accurate approximation. The mathematical expression for the noise level in CFAR-CA is given by
where
is the number of training cells on either side of the
kth target cell,
. While straightforward from an implementation standpoint, it can be less effective in environments with significant variation in terrain profile. Appropriately, CFAR-CASO is better suited for rough terrains. The mathematical expressions for noise level on the left and right side of the target cell in CFAR-CASO are given by [
3]
The smallest of the two averages is eventually chosen to be the noise level as follows:
The CFAR algorithm, by keeping a constant threshold, calculates the detection threshold as a function of the PFA [
47]. The scaling factor for estimating the threshold factor is expressed as
It is evident from (13) that a small value of PFA signifies a higher value of the scaling factor. Subsequently, the threshold for the
kth cell is given by
For all the target cells with signal amplitudes above the threshold, a target matrix containing binary elements is generated as
where
is the
kth target cell of the range profile with a size equal to the order of the range FFT,
. Like any detection method, CFAR-CASO is employed to validate potential target returns and to ensure that statistically significant targets are considered, effectively reducing false alarms.
Figure 8 shows a side-by-side comparison of the CFAR variants under discussion for a ground surface, offering five unique targets indicative of rough terrain. It is apparent that CFAR-CASO has a higher detection threshold but is still able to detect more targets. This indicates a lower PFA while still being able to detect more targets with a better approximation of noise. In conclusion, while CFAR-CA provides a straightforward approach to target detection, the ability of CFAR-CASO to handle clutter more effectively makes it the preferred choice for RAs operating over terrains with varying elevation profiles. Discussion on CFAR-CAGO has been deliberately skimmed due to its limited applicability to the intended use case.
The altimetry application requires that only a single valid target is reported as the altitude of the drone. Fittingly, it is essential to identify the strongest target return from a set of validated targets as the ground surface. Given that the radiation pattern of the IWR1843BOOST features a single grating lobe, the strongest return is expected from the ground surface at the radar boresight [
16]. The peak grouping methodology coupled with the antenna characteristics favors accurate altitude reporting.
5.2. Utility of SIMO for High Altitude
Albeit having high altitude requirements, having multiple receive antennae may be leveraged for AoA estimation. It can be particularly useful in filtering ground returns from off-boresight angles and better approximation of true altitude. Another inherent benefit of multiple Rx antennae is that even if the AoA capability is not utilized, integration of FMCW chirps along the antenna dimension improves SNR. Increased SNR is a desirable scenario for high-altitude estimation.
Since there is considerable signal processing overhead, it is appropriate to appraise the benefit of incorporating SIMO in the cruise stage.
Figure 9 illustrates the radar platform, IWR1843, onboard a drone cruising at an altitude of 180 m. A simplified scenario is elucidated, with a hilltop and flat terrain being point A and B, respectively. Given the minimum altitude requirement of 150 m, an arbitrary value of
is chosen to be 180 m. Appropriately, the Pythagorean theorem can be applied for the evaluation of the base length,
. For both points to be identified as separate targets in the angular domain,
is considered equal to angular resolution,
, 28.64°. Consequently, the requirement of minimum separation,
, between points A and B amounts to approximately 98.3 m. This shows that the 1 × 4 SIMO configuration offers little to no value in terms of angular information. The only benefit of having multiple antennae is, therefore, the SNR improvement through noncoherent integration along the antenna dimension. The 2D range-Doppler map in (4) was provided for a single antenna. The same expression can be extended for multiple antennae. The SNR-enhanced version commonly referred to as the detection matrix is mathematically expressed as
where
corresponds to range-Doppler maps for antenna index
, and
is the absolute value of the complex entries in the respective maps. The detection matrix is better suited for CFAR detection owing to the improved SNR, mathematically expressed as
where
and
signify the signal power before and after integration. As the phase is not summed coherently, the gain is scaled only by a factor of
.
8. Discussion
Having established the foundation for migration to mmWave altimetry for UASs in the preliminary work, this article aims to advance the discussion by providing a signal processing framework. The article maintains a tutorial approach to engage a broad readers, ranging from application engineers to seasoned academicians.
The framework for developing a comprehensive mmWave altimetry solution for UASs using automotive radars is focused on overcoming inherent hardware limitations, including Tx power and IF filter bandwidth. Both the preliminary work and this article focus on optimizing performance metrics for each stage of flight within these hardware constraints. Tailoring these metrics to meet the requirements of an mmWave altimeter with a stage adaptive waveform benefits greatly from the software-defined architecture. The contribution of this body of work is the meticulous documentation of engineering involved in waveform design by prioritizing metrics according to their relevance in the respective stage of flight. Contrary to the concept of legacy RA, all nine performance metrics listed in
Table 1 have been involved in the waveform design.
The signal processing aspects were elucidated, with each stage being simulated as a radar scenario. As highlighted in this study as well as the preliminary work, the intertwined nature of performance metrics requires careful consideration. For instance, the radial velocity component in the cruise stage does not have a significant impact as the target was empirically estimated to be always present in the 0th Doppler bin. However, it was observed that the target migrates to neighboring range bins due to the motion of the platform relative to the ground surface. Despite the radial velocity component arising from variations in the terrain profile and a chirp with a very long duration, the resultant range accuracy remained within the limits specified by operational requirements.
A novel contribution of this work was proposing TDM-MIMO for maximizing angular resolution as long as ROD remains within the limit in the landing approach stage. Owing to the relatively relaxed maximum range requirement, the chirp duration was reduced to cater to the maximum ROD. Fundamental concepts involved in Doppler compensation for a TDM-MIMO radar were presented mathematically for resolving phase ambiguity due to simultaneous motion and off-boresight angular position of the target. Lastly, in the touchdown stage, range accuracy was prioritized above all in the waveform design process since , AoA, and maximum range were largely irrelevant.
For the detection methodology, the suitability of three CFAR variants was appraised, with CFAR-CASO offering a reasonable compromise over others in the cruise and landing approach stage. The touchdown stage does not warrant CFAR detection owing to the straightforward characteristics of the runway. The underpinnings of CFAR, associated signal processing aspects, and the flow of operations were covered in reasonable detail from a mathematical standpoint.
8.1. Challenges
Having covered mmWave altimetry for UASs in great detail with supporting arguments, mathematical illustrations, and insights, it is now imperative to lay the groundwork for the culmination of the work. It is envisioned that the framework for deriving waveform specifications from operational requirements and simulation results must be validated in a real-world scenario. A two-stage approach for the execution of this yearning is to simulate a radar scenario with the target being the ground surface using the DTED of an actual runway. In the existing study, the targets were treated as uniform bodies, and a free space propagation model was employed. However, the true emulation of an actual scenario must consider the effect of specular reflections from the ground terrain. Afterward, it is planned to mount the radar on a small-sized drone and emulate all three stages of flight. Furthermore, a complete simulation of the landing scenario is planned using clutter generated from the actual DTED of runway surroundings. The advantage of employing TDM-MIMO in the determination of actual altitude, as opposed to a traditional RA, shall also be investigated in the upcoming research phase. This direction has groundbreaking potential to improve landing safety and maintain constant altitude in surveillance operations. The promising simulation results and alignment with practical operational requirements demonstrate the viability of mmWave altimetry using automotive radars, achieved through the optimization of performance metrics. Given the growing applications of UASs in modern life, the potential of mmWave altimetry for UASs warrants a dedicated MOPS tailored to specific operational needs.
8.2. Future Work
Despite promising outcomes, there are inherent obstacles to the pursuit of future goals. Firstly, the simulation scenarios covered in this work must be implemented on actual hardware. The zoomed-in range profile using the ZFFT algorithm illustrated in
Section 7 was generated using actual hardware, but the update rate was not part of the equation. During the touchdown estimation stage, the update rate is equally as critical as the range accuracy. The real-time reporting of altitude to FCC is necessary to engage subsequent control actions. Similarly, there is a need to implement requisite waveforms for all three flight stages. It is safe to assume that the cruise and landing approach entail relatively slow update rates and the real challenge resides in the touchdown stage. Accordingly, documenting the implementation details and the characterization of performance profiles comprising range accuracy and processing latency is mandatory. Another daunting challenge is managing the seamless transition between flight stages or, in other words, the engineering design of the hysteresis loop [
62]. The constant switching between waveform configurations at the borderline of altitude limits must be handled gracefully. An important goal from a systems engineering view is to ensure that failures and abnormal limits are well catered for at all times. One particular scenario that requires careful consideration is the drone exceeding the maximum limit of ROD in the landing approach stage. Consequently, Doppler compensation shall not be performed correctly, resulting in erroneous AoA. For a safety-critical application such as an emergency landing, this could lead to a catastrophic outcome if the landing relies solely on a single radar sensor. It is crucial to guarantee accurate reporting or no reporting at all, without any possibility of incorrect readings. One potential way out is to extend the
and resort to methods aimed at resolving the velocity ambiguity, employing methodology cited in a related dissertation [
63]. However, in the existing scope, there is a limit to extending the
, and a scenario of
being exceeded may arise regardless. In this context, a consolidated solution is mandatory that outlines the direction for future research. Lastly, it is envisioned that the summary of cited challenges, proposed solutions, and succinct details entailing this endeavor and preliminary work culminating into clutter simulation and experimental validation shall be documented in a comprehensive letter.