<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Sig. Proc.</journal-id>
<journal-title>Frontiers in Signal Processing</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Sig. Proc.</abbrev-journal-title>
<issn pub-type="epub">2673-8198</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">781777</article-id>
<article-id pub-id-type="doi">10.3389/frsip.2021.781777</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Signal Processing</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Multi-Frequency Radar Micro-Doppler Based Classification of Micro-Drone Payload Weight</article-title>
<alt-title alt-title-type="left-running-head">Dhulashia et&#x20;al.</alt-title>
<alt-title alt-title-type="right-running-head">Micro-Drone Payload Weight Based Classification</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Dhulashia</surname>
<given-names>Dilan</given-names>
</name>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1546552/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Peters</surname>
<given-names>Nial</given-names>
</name>
<uri xlink:href="https://loop.frontiersin.org/people/1546500/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Horne</surname>
<given-names>Colin</given-names>
</name>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Beasley</surname>
<given-names>Piers</given-names>
</name>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ritchie</surname>
<given-names>Matthew</given-names>
</name>
<uri xlink:href="https://loop.frontiersin.org/people/969906/overview"/>
</contrib>
</contrib-group>
<aff>
<institution>Department of Electronic and Electrical Engineering, University College London</institution>, <addr-line>London</addr-line>, <country>United&#x20;Kingdom</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1227373/overview">Luca Pallotta</ext-link>, Roma Tre University, Italy</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1506808/overview">Tao Shan</ext-link>, Beijing Institute of Technology, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1167638/overview">Christos Ilioudis</ext-link>, University of Strathclyde, United&#x20;Kingdom</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Dilan Dhulashia, <email>dilan.dhulashia.15@ucl.ac.uk</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Radar Signal Processing, a section of the journal Frontiers in Signal Processing</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>22</day>
<month>12</month>
<year>2021</year>
</pub-date>
<pub-date pub-type="collection">
<year>2021</year>
</pub-date>
<volume>1</volume>
<elocation-id>781777</elocation-id>
<history>
<date date-type="received">
<day>23</day>
<month>09</month>
<year>2021</year>
</date>
<date date-type="accepted">
<day>29</day>
<month>11</month>
<year>2021</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2021 Dhulashia, Peters, Horne, Beasley and Ritchie.</copyright-statement>
<copyright-year>2021</copyright-year>
<copyright-holder>Dhulashia, Peters, Horne, Beasley and Ritchie</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these&#x20;terms.</p>
</license>
</permissions>
<abstract>
<p>The use of drones for recreational, commercial and military purposes has seen a rapid increase in recent years. The ability of counter-drone detection systems to sense whether a drone is carrying a payload is of strategic importance as this can help determine the potential threat level posed by a detected drone. This paper presents the use of micro-Doppler signatures collected using radar systems operating at three different frequency bands for the classification of carried payload of two different micro-drones performing two different motions. Use of a KNN classifier with six features extracted from micro-Doppler signatures enabled mean payload classification accuracies of 80.95, 72.50 and 86.05%, for data collected at S-band, C-band and W-band, respectively, when the drone type and motion type are unknown. The impact on classification performance of different amounts of situational information is also evaluated in this&#x20;paper.</p>
</abstract>
<kwd-group>
<kwd>micro-drone</kwd>
<kwd>multi-frequency radar</kwd>
<kwd>micro-Doppler</kwd>
<kwd>payload classification</kwd>
<kwd>counter-UAV</kwd>
<kwd>Doppler signatures</kwd>
<kwd>multi-functional radar</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Over the past decade there has been a substantial increase in the number of micro-drones, i.e.,&#x20;Unmanned Aerial Vehicles (UAVs), available on both the commercial and consumer markets. This increase can be attributed to their continually improving capabilities, ease of use, and low cost, making them attractive to businesses and hobbyists alike. There have been countless positive applications of micro-drone technology, e.g., agricultural and environmental surveying, disaster response and search and rescue. However, micro-drones have become infamous for their misuse, in performing anti-social, and often criminal activities, such as illegal filming of private/sensitive locations or dangerous flying through restricted airspace creating a collision hazard to other aircraft. Of particular interest in this work is the ability for micro-drones to be utilized as vehicles for transporting payloads of dangerous substances e.g. explosives, chemical/biological weapons and illicit drugs. This capability is of concern to security agencies worldwide (<xref ref-type="bibr" rid="B9">Forbes, 2019</xref>), with calls being made to limit commercial drone payload capacities in order to limit the risk posed (<xref ref-type="bibr" rid="B1">BBC, 2016</xref>).</p>
<p>Current conventional radar based air surveillance systems often perform poorly when tasked with detecting and tracking micro-drones, due to their considerably smaller Radar Cross Section (RCS) when compared with standard aircraft, thus establishing the need for a dedicated system with adequate sensitivity for detection of these small RCS targets (<xref ref-type="bibr" rid="B20">Rahman and Robertson, 2018</xref>). Additionally, birds present an issue with drone detection systems due to their comparable RCS, altitude and speed of flight, all of which make discrimination between these two targets challenging using RCS information alone (<xref ref-type="bibr" rid="B20">Rahman and Robertson, 2018</xref>). In order to reduce the false alarm rate and thus mature radar based drone detection technology, there is a widening body of research into methods for discriminating between birds and micro-drones. The micro-Doppler signatures of drones have been identified as one method for discriminating micro-drones from other targets and general clutter (<xref ref-type="bibr" rid="B24">Tahmoush, 2014</xref>; <xref ref-type="bibr" rid="B22">Ritchie et&#x20;al., 2016</xref>; <xref ref-type="bibr" rid="B19">Samiur and Robertson, 2019</xref>; <xref ref-type="bibr" rid="B2">Beasley et&#x20;al., 2020</xref>). Micro-Doppler signatures are produced by periodic micro-motions of a target or object, such as flapping of a bird&#x2019;s wings, the gait of a human walk or rotation of a micro-drone&#x2019;s propellers (<xref ref-type="bibr" rid="B3">Chen, 2019</xref>). Micro-Doppler is the respective Doppler shift induced by these target micro-motions which appear as side bands around the target&#x2019;s bulk Doppler. Further details on the mechanisms that produce micro-Doppler have been extensively documented for a range of radar targets in <xref ref-type="bibr" rid="B25">Miceli (2014)</xref> and <xref ref-type="bibr" rid="B3">Chen (2019)</xref>. There are several publications detailing micro-Doppler signatures of drones, where comparisons are made with bird micro-Doppler e.g., <xref ref-type="bibr" rid="B24">Tahmoush (2014)</xref>, <xref ref-type="bibr" rid="B22">Ritchie et&#x20;al. (2016)</xref>, <xref ref-type="bibr" rid="B20">Rahman and Robertson (2018)</xref>, and <xref ref-type="bibr" rid="B2">Beasley et&#x20;al. (2020)</xref>. This work has been taken further in <xref ref-type="bibr" rid="B13">Molchanov et&#x20;al. (2014)</xref>, <xref ref-type="bibr" rid="B12">Kim et&#x20;al. (2017)</xref>, <xref ref-type="bibr" rid="B19">Samiur and Robertson (2019)</xref>, and <xref ref-type="bibr" rid="B17">Rahman and Robertson (2020)</xref>, where drone micro-Doppler signatures are fed into classifiers in order to discriminate between both different classes of drones and between drones and birds. For a more comprehensive review of the current literature on drone detection and classification using machine learning, please refer to a recent paper published on the topic by <xref ref-type="bibr" rid="B23">Taha and Shoufan (2019)</xref>.</p>
<p>There have been a number of works published on multi-frequency micro-Doppler based classification of targets, with gesture recognition accounting for a large portion. Such works look to compare the suitability of particular frequency bands in the role of micro-Doppler signature generation for classification. A recent paper by <xref ref-type="bibr" rid="B10">Gurbuz et&#x20;al. (2021)</xref>, investigated the use of three frequencies S-bands, namely 77&#xa0;GHz, 24 and 10&#xa0;GHz, for use in sign language recognition. It was reported that, in nearly all cases, using the 77&#xa0;GHz radar data resulted in the highest classification accuracy, where the lowest classification accuracy was reported to result from using the 10&#xa0;GHz radar&#x20;data.</p>
<p>A comparatively more limited body of literature is published on the topic of micro-Doppler based classification of unloaded and loaded drones. A single data-set collected by University College London&#x2019;s (UCL) multistatic S-Band radar system, NETRAD (<xref ref-type="bibr" rid="B5">Doughty, 2008</xref>), has been used in a number of publications (<xref ref-type="bibr" rid="B7">Fioranelli et&#x20;al., 2015</xref>; <xref ref-type="bibr" rid="B21">Ritchie et&#x20;al., 2017</xref>; <xref ref-type="bibr" rid="B15">Patel et&#x20;al., 2019</xref>; <xref ref-type="bibr" rid="B14">Pallotta et&#x20;al., 2020</xref>), to test the accuracy of a variety of payload classification algorithms. In this data-set, a DJI Phantom 2 was loaded with a variety of payloads between 0 and 600&#xa0;g and made to either hover or fly towards the central radar node. <xref ref-type="bibr" rid="B7">Fioranelli et&#x20;al. (2015)</xref> and <xref ref-type="bibr" rid="B21">Ritchie et&#x20;al. (2017)</xref> report high classification accuracies using a feature extraction based approach. In the hovering scenario, accuracies of 96&#x2013;97% were achieved using centroid based micro-Doppler features. In contrast, for payload classification of the drone in flight, a maximum classification accuracy of 95% was achieved when using singular value decomposition (SVD) derived features. A similar feature based approach was used on the same data-set in a more recent publication by <xref ref-type="bibr" rid="B14">Pallotta et&#x20;al. (2020)</xref>, differing mainly through the use of spectral kurtosis (<xref ref-type="bibr" rid="B6">Dwyer, 1983</xref>) to extract features from the spectrograms. The algorithm produced by <xref ref-type="bibr" rid="B14">Pallotta et&#x20;al. (2020)</xref> achieved classification accuracies for the hovering drone of 92.61%, marginally lower than the result reported by <xref ref-type="bibr" rid="B21">Ritchie et&#x20;al. (2017)</xref>. In the case of the moving drone, their classifier achieved an accuracy of 70.83%, considerably lower than the 95% accuracy reported in the work by <xref ref-type="bibr" rid="B21">Ritchie et&#x20;al. (2017)</xref>. The UCL data-set is again used by <xref ref-type="bibr" rid="B15">Patel et&#x20;al. (2019)</xref>; however, a different approach is used to tackle the classification problem. In <xref ref-type="bibr" rid="B7">Fioranelli et&#x20;al. (2015)</xref>, <xref ref-type="bibr" rid="B21">Ritchie et&#x20;al. (2017)</xref>, and <xref ref-type="bibr" rid="B14">Pallotta et&#x20;al. (2020)</xref> work, a single spectrogram time-frequency representation was used, whereas in the work by <xref ref-type="bibr" rid="B15">Patel et&#x20;al. (2019)</xref>, two additional transforms are used, namely the cepstrogram and the cadence velocity diagram (CVD). Instead of selecting individual features to feed a classifier, a Convolutional Neural Network (CNN), AlexNet, is provided with the three time-frequency representation images. In both the hovering and moving drone scenarios, the AlexNet based image classifier had an overall accuracy of around 95&#x2013;96%, comparable to the results reported by <xref ref-type="bibr" rid="B7">Fioranelli et&#x20;al. (2015)</xref> and <xref ref-type="bibr" rid="B21">Ritchie et&#x20;al. (2017)</xref>.</p>
<p>
<xref ref-type="bibr" rid="B18">Rahman et&#x20;al. (2020)</xref> from The University of St Andrews have since released a substantial piece of work in which a 24&#xa0;GHz Frequency Modulated Continuous Wave (FMCW) radar and 96&#xa0;GHz Continuous Wave (CW) radar were simultaneously used to sample the micro-Doppler signatures of two large drones (DJI S900 hexacopter and a Joyance JT5L-404 crop spraying quadcopter) with a variety of payloads up to a maximum of 5&#xa0;kg. <xref ref-type="bibr" rid="B18">Rahman et&#x20;al. (2020)</xref> concluded through visual inspection of the spectrograms derived from the FMCW radar data that there were no consistent unique signatures that could solely be attributed to a drone carrying a payload. However, more success was had in the case of the CW radar data, though a conclusion was still reached that observed features could not be used as the sole basis of a classification algorithm. Nevertheless, an attempt was made to classify between unloaded and loaded drones using two machine learning approaches, namely feature extraction and a CNN. However, both approaches were found to have limited success.</p>
<p>The work presented in this paper looks to expand on previous work by comparing the effectiveness of radar data captured at three different frequency bands to be used for determining the payload weight which different drones are carrying. Radar systems operating at C-band, S-band and W-band (mm-wave) are used to measure two different types of small quadcopter type drones carrying three different payload weights and performing two different movement styles each. A series of investigations with increasing complexity are then carried out using the micro-Doppler signatures obtained from the radar measurements to be able to classify the payload type. The rest of this paper is organized as follows. <xref ref-type="sec" rid="s2">Section 2</xref> provides an overview of the experimental setup and the equipment used, including information about the radar systems and drones. <xref ref-type="sec" rid="s3">Section 3</xref> presents the theory and the methodology of the signal processing and the classification procedures used. <xref ref-type="sec" rid="s4">Section 4</xref> reports the results obtained and provides some analysis from these. Discussion on the observations are then summarized in <xref ref-type="sec" rid="s5">Section 5</xref>, before conclusion are finally presented in <xref ref-type="sec" rid="s6">Section&#x20;6</xref>.</p>
</sec>
<sec id="s2">
<title>2 Experimental Setup</title>
<p>In this section, the two radar platforms used in this research to capture radar data, the C-band and S-band ARESTOR RF sensor system and the W-band Texas Instruments radar system, are described. Details of the drone models flown as targets and the experimental setup employed to gather the radar data are then provided.</p>
<sec id="s2-1">
<title>2.1&#x20;C-Band and S-Band RF Sensor</title>
<p>The UCL developed ARESTOR RF sensor platform is used to collect S-band and C-band radar data from the UAV under test. ARESTOR is a highly flexible custom-built RF sensor system, based on the Xilinx RF System on a Chip (RFSoC) <xref ref-type="bibr" rid="B27">Xilinx (2021b)</xref>, which integrates multiple RF sampling rate analogue to digital converters and digital to analogue converters with a large field-programmable gate array (FPGA) fabric and multiple ARM processing&#x20;cores.</p>
<p>ARESTOR is an eight channel transmit/eight channel receive RF sensor capable of direct sampling of RF signals up to S-band frequencies. The system consists of a Xilinx ZCU111 RFSoC evaluation board <xref ref-type="bibr" rid="B26">Xilinx (2021a)</xref>, shown in <xref ref-type="fig" rid="F1">Figure&#x20;1</xref>, plus FPGA images and software infrastructure which allows the customization of the ZCU111 such that it implements a range of radar and RF sensor components, including FMCW radar channels, passive radar channels and wideband electronic surveillance channels.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Xilinx ZCU111 RFSoC evaluation board with attached RF interface card.</p>
</caption>
<graphic xlink:href="frsip-01-781777-g001.tif"/>
</fig>
<p>Synchronization between multiple evaluation boards allows the construction of systems with multiple Tx/Rx coherent channels beyond the eight provided by a single board, and, with the addition of our in-house GPSDO solution, allows for multiple, spatially separated, ARESTOR subsystems to be synchronized for coherent multistatic operation.</p>
<p>For a more detailed overview of the ARESTOR RF sensor platform, see <xref ref-type="bibr" rid="B16">Peters et&#x20;al. (2021)</xref>.</p>
<p>In this study, the sensor is configured for dual-frequency FMCW radar operation, with one channel operating within the native frequency range of the RFSoC device for the S-band capture, and for C-band by extending the frequency range using a prototype RF front end which extends the upper limits of the frequency coverage to approximately 12&#xa0;GHz.</p>
<p>The S-band channel uses a pair of 18&#xa0;dBi dish antennas, and the C-band channel uses a pair of ultra-wideband tapered slot antennas, each providing 12 dBi of&#x20;gain.</p>
<p>Radar operating parameters employed in the data captures are shown in <xref ref-type="table" rid="T1">Table&#x20;1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Operating parameters of radars and waveforms&#x20;used.</p>
</caption>
<table>
<thead>
<tr>
<td align="left">Parameter</td>
<td align="center">W-Band</td>
<td align="center">C-Band</td>
<td align="center">S-Band</td>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Operating mode</td>
<td align="center">FMCW</td>
<td align="center">FMCW</td>
<td align="center">FMCW</td>
</tr>
<tr>
<td align="left">Waveform</td>
<td align="center">Rising sawtooth</td>
<td align="center">Rising sawtooth</td>
<td align="center">Rising sawtooth</td>
</tr>
<tr>
<td align="left">PRF (Hz)</td>
<td align="center">5,000</td>
<td align="center">5,000</td>
<td align="center">5,000</td>
</tr>
<tr>
<td align="left">Chirp length (us)</td>
<td align="center">160</td>
<td align="center">100</td>
<td align="center">100</td>
</tr>
<tr>
<td align="left">Carrier frequency (GHz)</td>
<td align="center">77.00</td>
<td align="center">5.80</td>
<td align="center">2.44</td>
</tr>
<tr>
<td align="left">RF bandwidth (MHz)</td>
<td align="center">2,560</td>
<td align="center">80</td>
<td align="center">80</td>
</tr>
<tr>
<td align="left">Duty cycle (%)</td>
<td align="center">80</td>
<td align="center">50</td>
<td align="center">50</td>
</tr>
<tr>
<td align="left">Tx power (dBm)</td>
<td align="center">12</td>
<td align="center">20</td>
<td align="center">20</td>
</tr>
<tr>
<td align="left">Tx gain (dBi)</td>
<td align="center">9</td>
<td align="center">12</td>
<td align="center">18</td>
</tr>
<tr>
<td align="left">Rx gain (dBi)</td>
<td align="center">30</td>
<td align="center">12</td>
<td align="center">18</td>
</tr>
<tr>
<td align="left">Sampling rate (ksps)</td>
<td align="center">6,250</td>
<td align="center">1875 (post mixing and decimation)</td>
<td align="center">1875 (post mixing and decimation)</td>
</tr>
<tr>
<td align="left">Polarisation</td>
<td align="center">Horizontal</td>
<td align="center">Horizontal</td>
<td align="center">Horizontal</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s2-2">
<title>2.2&#x20;W-Band Radar</title>
<p>A Texas Instruments W-band radar was used to collect the W-band data during the drone captures. This system is comprised of an AWR1642BOOST radar evaluation board and a Texas Instruments DCA1000EVM capture card together being controlled from a PC using the Texas Instruments mmWave Studio application. The main radar operating parameters used and the parameters used for the chirped waveform design are shown in the corresponding columns within <xref ref-type="table" rid="T1">Table&#x20;1</xref>. A picture of the system, including the radar board and the capture card is shown in <xref ref-type="fig" rid="F2">Figure&#x20;2</xref>.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Texas Instruments W-band radar system including AWR1642BOOST radar evaluation board and DCA1000EVM capture&#x20;card.</p>
</caption>
<graphic xlink:href="frsip-01-781777-g002.tif"/>
</fig>
<p>This radar is relatively simple to use out of the box, requiring no bespoke design from the user end beyond waveform design. While the radar includes two transmit channels and four receive channels, only a single transmit channel was used with the four receive channels in the collection of the data used in this study. The system can be used in either a vertical or a horizontal polarization by physical manipulation of the antenna orientations. The transmit and receive antennas included on the system and used are formed of etched patch designs. In the horizontal polarization used, each receive antenna is comprised of a series of connected patches in a horizontal alignment. The use of four receive antennas which are vertically stacked allows collection of four data streams which can be combined prior to Doppler-time processing in order to achieve an improved SNR with minimal distortion within resultant Doppler signatures. This is due to the relatively small angle differences between each receive antenna and a target. The transmit channel of the radar provides a total gain of 9&#xa0;dBi, while each receive channel provides a gain of 30&#xa0;dBi. A similar version to this system has been proposed and used within <xref ref-type="bibr" rid="B11">Gusland et&#x20;al. (2021)</xref>, where an initiative has been formed in order to collect and share radar data within the radar community.</p>
</sec>
<sec id="s2-3">
<title>2.3 Drone Types</title>
<p>Two different quadcopter style drones were used in the captures performed for this research. These are the DJI Phantom 2 and the DJI Spark. Both drones have a target user base of hobbyists and are two of the most common types of low cost quadcopter style drones which are typically used for recreational purposes.</p>
<p>The DJI Spark is the smaller of the two drones used, with main body dimensions of 143&#x20;&#xd7; 143&#x20;&#xd7; 55&#xa0;mm, a diagonal length (excluding blades) of 175&#xa0;mm, blade length of 110&#xa0;mm and a mass (including blades and battery) of 300&#xa0;g. It has a maximum velocity of 14&#xa0;ms<sup>&#x2212;1</sup> and control communications are at a frequency of 2.400&#x2013;2.483&#xa0;GHz or 5.725&#x2013;5.825&#xa0;GHz. This drone will be referred to as the &#x201c;small drone&#x201d; in the following sections of this paper. An image of the drone can be seen in <xref ref-type="fig" rid="F3">Figure&#x20;3A</xref>. The payload masses used with this drone are: 0, 50 and 100&#xa0;g. As proportions of the mass of the drone, these payloads equate to: 0, 16.67 and 33.33%, respectively.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>The two drones used within tests in this research: <bold>(A)</bold> DJI Spark, smaller drone <bold>(B)</bold> DJI Phantom 2, larger&#x20;drone.</p>
</caption>
<graphic xlink:href="frsip-01-781777-g003.tif"/>
</fig>
<p>The DJI Phantom 2 is the larger of the two drones used, with dimensions of 300&#x20;&#xd7; 300&#x20;&#xd7; 190&#xa0;mm, a diagonal length (excluding blades) of 425&#xa0;mm, blade length of 240&#xa0;mm and a mass (including blades and battery) of 1,250&#xa0;g. It has a maximum velocity of 15&#xa0;ms<sup>&#x2212;1</sup> and control communications are at a frequency of 2.400&#xa0;GHz. Unlike the Spark, the Phantom includes two lower rails beneath the main body of the drone which act as landing gear, providing it with a distinctive difference in shape from the small drone. This drone will be referred to as the &#x201c;big drone&#x201d; in the following sections of this paper. An image of the drone can be seen in <xref ref-type="fig" rid="F3">Figure&#x20;3B</xref>. The payload masses used with this drone are: 0, 250 and 500&#xa0;g. As proportions of the mass of the drone, these payloads equate to: 0, 20 and 40%, respectively.</p>
<p>It should be noted that the communications signal used for controlling the smaller drone is always within the same frequency range as either the C-band or S-band radar systems. Likewise, the communications signal used for controlling the big drone is always within the same frequency range as the S-band radar system. No interference from either of the drone communications signals were observed within the data collected.</p>
</sec>
<sec id="s2-4">
<title>2.4 Data Captures</title>
<p>The drone captures were carried out on the same day in a flat-terrain outdoor setting over short grass with minimal other electromagnetic (EM) or non-drone target interference. A diagram depicting the approximate positions of the drone relative to the antennas of the monostatic radars and the approximate distances over which the drone was allowed to fly can be seen in <xref ref-type="fig" rid="F4">Figure&#x20;4</xref>.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Depiction of layout for experimental setup used and distances over which drone flights were carried&#x20;out.</p>
</caption>
<graphic xlink:href="frsip-01-781777-g004.tif"/>
</fig>
<p>Two captures, each of 1&#xa0;min duration, were carried out for each drone, motion and payload combination, meaning a total of 24 captures per frequency band (two drone types, two motions and three different payloads, each recorded twice). Radar data was captured simultaneously at the three different frequency bands during each capture. A final set of captures were carried out to record a static target with no drone present in the environment at S-band using the same operating parameters as the S-band parameters used for the drone data collection, shown in <xref ref-type="table" rid="T1">Table&#x20;1</xref>.</p>
</sec>
</sec>
<sec id="s3">
<title>3 Signal Processing and Classification</title>
<p>In this section, the signal processing methodology used to facilitate the generation of Doppler-time spectrograms is presented. The extracted features and machine learning model used for classification are described, and details of the classification investigations for which results will be reported are&#x20;given.</p>
<sec id="s3-1">
<title>3.1 Signal Processing</title>
<p>An overview of the signal processing performed on the captured data sets is shown in the block diagram provided in <xref ref-type="fig" rid="F5">Figure&#x20;5</xref>, in which data structures are contained within blue blocks and operational stages are contained within red blocks.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Block diagram showing stages involved in data processing from collected raw data to features used for classification. Data structures are contained within blue blocks and operational stages are contained within red blocks.</p>
</caption>
<graphic xlink:href="frsip-01-781777-g005.tif"/>
</fig>
<p>For each frequency band, the data obtained within a capture is initially organized into samples recorded for each returned chirp, commonly referred to by fast and slow time samples. A range fast Fourier transform (FFT) is performed on the samples within each chirp to obtain the range-time data. A moving target indication (MTI) algorithm is used on this range-time data to remove returns from static targets within the vicinity and to enable a greater concentration onto the response from the drone. This stage is desirable as returns from static targets can produce a prominent component at the 0&#xa0;Hz Doppler frequency within the spectrograms which are used for classification. A strong component at this Doppler frequency has the potential to dominate over intricate micro-Doppler components resulting from reflections from the drone within the spectrogram and would cause the features used for classification to be based upon the return from the stationary objects instead of the drone. A range-time sample synchronization process is then used in order to remove capture initiation offsets between data collected from the three frequency bands to ensure features within the generated spectrograms will arise from the same observed physical phenomena. Following this, the range bins in which the drone target is determined to be located in at any given time are collected and formed into a single dimensional time spanning vector, <italic>R</italic>
<sub>&#x3a3;</sub>. A short time Fourier transform (STFT) can be performed across these bins to generate a Doppler-time spectrogram containing the Doppler signature of the drone. The value of the STFT at frequency bin <italic>k</italic> at time slice <italic>m</italic> is then given by <xref ref-type="disp-formula" rid="e1">Eq. 1</xref>, where <italic>L</italic> is the number of samples in a window length, <italic>x</italic> [<italic>i</italic>] is the <italic>i</italic>th sample within <italic>R</italic>
<sub>&#x3a3;</sub>, and <italic>h</italic>
<sub>
<italic>i</italic>
</sub> is the <italic>i</italic>th sample within the weighting window used for tapering the short time samples.<disp-formula id="e1">
<mml:math id="m1">
<mml:mi>F</mml:mi>
<mml:mfenced open="[" close="]">
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mo>&#x3d;</mml:mo>
<mml:munderover accentunder="false" accent="false">
<mml:mrow>
<mml:mo>&#x2211;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:munderover>
<mml:msub>
<mml:mrow>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>x</mml:mi>
<mml:mfenced open="[" close="]">
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>j</mml:mi>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>L</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
<p>The parameters used to perform the STFT on the range-time data for all the datasets used in this research are shown in <xref ref-type="table" rid="T2">Table&#x20;2</xref>.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Short time Fourier transform parameters used to generate Doppler-time spectrograms.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Parameter</th>
<th align="center">Value</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Window type</td>
<td align="center">Hamming</td>
</tr>
<tr>
<td align="left">Window length (samples)</td>
<td align="center">500</td>
</tr>
<tr>
<td align="left">Window length (% of PRI)</td>
<td align="center">10</td>
</tr>
<tr>
<td align="left">Overlap length (samples)</td>
<td align="center">450</td>
</tr>
<tr>
<td align="left">Overlap factor (% of window length)</td>
<td align="center">90</td>
</tr>
<tr>
<td align="left">FFT Points (samples)</td>
<td align="center">4,000</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The resultant Doppler-time spectrogram for the drone motion throughout the capture is then cut into a series of one second long segments.</p>
</sec>
<sec id="s3-2">
<title>3.2 Feature Extraction</title>
<p>Feature extraction is carried out on each of the Doppler-time segments and the obtained values for each of the features are recorded into feature vectors along with a class label of the weight of the payload. It should be noted here that a Doppler-time segment (and the use of the word &#x201c;segment&#x201d;) refers to a single, one second long spectrogram slice containing values in decibels. The use of one second long spectrogram samples has been shown to be an effective length for drone micro-Doppler classification in <xref ref-type="bibr" rid="B21">Ritchie et&#x20;al. (2017)</xref>. This short sample length ensures that feature values are based on micro-movements intrinsic to the drone instead of events such as: variations in flight path, the amount the drone moves within the 1&#xa0;min capture, or other unrelated phenomena occurring within the scene during the capture period. If a one second long spectrogram slice is represented as matrix <bold>D</bold>
<sub>complex</sub>, then the Doppler-time segment upon which feature extraction is carried out is given by <bold>D</bold>
<sub>dB</sub>:<disp-formula id="e2">
<mml:math id="m2">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="bold">D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mtext>dB</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>20</mml:mn>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>log</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>10</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mfenced open="|" close="|">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="bold">D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mtext>complex</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:math>
<label>(2)</label>
</disp-formula>
</p>
<p>A comprehensive list of the features which are extracted from Doppler-time segments and used for classification are presented in <xref ref-type="table" rid="T3">Table&#x20;3</xref>. The features extracted from the Doppler-time segments are based on statistical values and empirical characteristics obtained from the segments. These provide measures of the power distributions and variations within spectrogram samples. The feature set is comprised of some features which have previously been used in <xref ref-type="bibr" rid="B4">Dhulashia et&#x20;al. (2021)</xref> for the purpose of human activity micro-Doppler signature classification, simple statistical measures of the segment, and an application specific feature well suited to finding the mean bulk Doppler frequency for drone targets.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Descriptions of features extracted from Doppler-time segments which comprise feature vectors used for classification.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Feature</th>
<th align="center">Feature description</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">1</td>
<td align="left">Mean of mean values from each time sample within segment</td>
</tr>
<tr>
<td align="left">2</td>
<td align="left">Variance of all values within segment</td>
</tr>
<tr>
<td align="left">3</td>
<td align="left">Median power of all values within segment</td>
</tr>
<tr>
<td align="left">4</td>
<td align="left">Maximum power of all values within segment</td>
</tr>
<tr>
<td align="left">5</td>
<td align="left">Mean of standard deviation of values in each time sample within segment</td>
</tr>
<tr>
<td align="left">6</td>
<td align="left">Mean bulk Doppler frequency within segment</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The feature vectors are stored within categorized tables based on the drone type and movement type. All of the feature vectors for each of the four combinations (two drone types and two movement types) formed from these categories from all payload captures are stored for each of the frequencies examined in the work. The feature vector collections for each frequency can then be used to train a classifier upon choice of movement type and/or drone&#x20;type.</p>
</sec>
<sec id="s3-3">
<title>3.3&#x20;K-Nearest Neighbors Classifier</title>
<p>The classifier model used in this research is the K-Nearest Neighbor (KNN) algorithm. This model is trained by mapping the training data into a six-dimensional feature space based on the extracted feature values for each data segment. Test samples are similarly mapped into the feature space. The class types of the &#x201c;<italic>k</italic>&#x201d; nearest training data points to the test data point within the feature space are known and a vote is made amongst these to find the majority class type. The test data point is then classified as being a member belonging to this class type. In all classification tests included in this research, the number of neighbors over which voting is carried out (i.e.,&#x20;<italic>k</italic> value) is three. The KNN method for classification is suitable as it is non-parametric and therefore makes no assumptions regarding the probabilistic distribution of the data, as well as being well suited to problems where the data sets consist of smaller amounts of data. The choice to use three nearest neighbors for class label voting is to reduce computation times while minimizing the impact on classification of erroneous or low quality individual training data samples as well as potential overfitting.</p>
</sec>
<sec id="s3-4">
<title>3.4 Classification Methodology</title>
<p>In a fully deployable radar system, for usage outside of the experimental research context, a hierarchical classification procedure would typically be used. The initial stage of such a classification methodology would perform target detection and recognition. Here, a distinction would be made regarding the type of target which is being observed. This would enable the system to determine whether the observed target is a drone or some other similar RCS target of non-interest, for example, a bird. Such target determination methods are presented in <xref ref-type="bibr" rid="B22">Ritchie et&#x20;al. (2016)</xref> and <xref ref-type="bibr" rid="B20">Rahman and Robertson (2018)</xref>. If the system determines that the target is a drone, a second classification stage would then be employed in order to classify the payload status of the drone. In the context of the work presented in this paper, a single classification stage is used. It is assumed that the target is present and that the target is already known to be a&#x20;drone.</p>
<p>Since the classification methodology employed in this work did not use a hierarchical classification scheme, the data samples from the static target captures allows for the inclusion of a class of type, &#x201c;No Drone,&#x201d; within the training set data when constructing classifiers used for the investigations within this research. The purpose of this is to ensure that if a provided test sample is of sufficiently low quality, for example, if the drone was not observed, the classifier will be able to determine this as being the case instead of forcefully attempting to classify the sample into one of the three payload based classes pertaining to the drone being present.</p>
<p>The research presented in the results section of this paper is broken down into four classification investigations based on supervised learning problems of differing complexities. The first classification investigation aims to classify the payload carried by a drone, given that the drone type and drone motion is known. This classification is found using data from each frequency band separately. This is achieved by creating classifier training data sets using one second spectrogram segments obtained from the captures for all three payloads but only for a particular drone, motion, and frequency band combination. 80 samples per payload mass are used to construct these training data sets, as well as an additional 40 samples of type &#x201c;No Drone,&#x201d; which correspond to samples taken when no drone was flying, resulting in a total training set size of 280 samples. Unseen test samples corresponding to the same drone and motion combination but of varying payload are then provided to the classifier for classification. A total of 2,000 test samples are used in each classification test in a Monte-Carlo fashion. Within a single classification test type, the classifier model will be retrained with a new selection of randomly chosen training data several times in order to allow the number of Monte-Carlo tests to be reached. The same format of Monte-Carlo tests along with the same number of repetitions are also used in investigations two, three and four. <xref ref-type="fig" rid="F6">Figure&#x20;6</xref> shows the feature table banks from which data is randomly selected depending on the classification test type and how data is provided to the classifier. It should be noted that the algorithm used for randomized file selection ensures test data remains unseen by the classifier during training and that only samples corresponding to the test type being carried out are&#x20;used.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Diagram depicting stages involved in selection of feature vectors from banks to results storage.</p>
</caption>
<graphic xlink:href="frsip-01-781777-g006.tif"/>
</fig>
<p>The second classification investigation aims to classify the payload being carried by a known drone type but with unknown motion. As such, training and test sets are both created using a mixture of data from both &#x201c;Hovering&#x201d; and &#x201c;Moving&#x201d; motion types for each of the frequency bands being investigated, as well as the training set including &#x201c;No Drone&#x201d; samples. Similarly to the first investigation, 80 samples per motion type per drone type are included in each training set, along with 40 samples of type &#x201c;No Drone,&#x201d; meaning training sets have a total size of 520 samples.</p>
<p>The third classification investigation looks to classify the payload being carried by an unknown drone type but with known motion. This is performed for a given frequency band by creating classifier training sets comprised of samples for all payloads from both the big and small drone but for only a single motion type, as well as 40 type &#x201c;No Drone&#x201d; samples. The total training set size is then 520 samples, and the classification tests are repeated for each of the two motion&#x20;types.</p>
<p>The fourth classification investigation is concerned with the classification of the payload given that both the drone and motion types are unknown. While the numerical value of the payload weights is different for both drones, the use of the payload labels: &#x201c;Empty,&#x201d; &#x201c;Light&#x201d; and &#x201c;Heavy,&#x201d; ensures that the problem remains a four class problem (also including the &#x2018;No Drone class). The training set for this investigation for a given frequency band consists of 80 samples from each of the payload, drone type and motion type combinations, as well as the 40 samples from the &#x201c;No Drone&#x201d; class, making the total training set contain 1,000 samples.</p>
<p>In all four investigations, it is ensured that the test sets are balanced such that an equal number of samples from each of the true payload classes are provided, as well as equal numbers from each motion type and/or drone type, in the cases where either (or both) of these variables are considered unknown.</p>
</sec>
</sec>
<sec id="s4">
<title>4 Results</title>
<p>This section first presents examples of Doppler-time spectrograms obtained from captures using two different payload weights at each of the three frequency bands before presenting the results from the different classification tests in the order introduced in the preceding section.</p>
<p>A series of example Doppler-time spectrograms corresponding to 1&#xa0;min of radar data are shown in <xref ref-type="fig" rid="F7">Figure&#x20;7</xref>. The spectrograms show the differences between data collected at the three different frequency bands being used to generate a Doppler signature corresponding to the same capture. Spectrograms are also shown for two different payloads: No payload and heavy payload (500&#xa0;g). All the example spectrograms are generated from data collected using the big drone and the &#x201c;Moving&#x201d; motion&#x20;type.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Example Doppler signatures from recordings of large drone moving: <bold>(A)</bold> S-band, 0&#xa0;g payload <bold>(B)</bold> S-band, 500&#xa0;g payload <bold>(C)</bold> C-band, 0&#xa0;g payload <bold>(D)</bold> C-band, 500&#xa0;g payload <bold>(E)</bold> W-band, 0&#xa0;g payload <bold>(F)</bold> W-band, 500&#xa0;g payload.</p>
</caption>
<graphic xlink:href="frsip-01-781777-g007.tif"/>
</fig>
<p>The initial classification investigation aimed to classify the weight of the carried payload by the drone given that the training set data was formed only from samples corresponding to the same drone and motion type as test samples being used, along with samples from a static target. <xref ref-type="fig" rid="F8">Figure&#x20;8</xref> shows the results from this classification investigation for both drone and motion types at each of the three frequency&#x20;bands.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Payload classification accuracy from data collected at each frequency band for a known drone and motion&#x20;type.</p>
</caption>
<graphic xlink:href="frsip-01-781777-g008.tif"/>
</fig>
<p>It can be seen from <xref ref-type="fig" rid="F8">Figure&#x20;8</xref> that when the motion type is &#x201c;Hovering,&#x201d; higher payload classification accuracies were achieved using S-band and C-band when the small drone was used, whereas a higher classification accuracy was achieved using W-band when the big drone was used. When the big drone is used, the performance ranking of the three frequency bands, from best to worst is: W-band, C-band, S-band. When the small drone is used, the ranking is reversed: S-band, C-band, W-band. This shows that despite the smaller payload weight to body weight ratio, when the smaller drone is used, the payload classification problem is still easier for the small drone when the motion type is &#x201c;Hovering&#x201d; and either S-band or C-band are being&#x20;used.</p>
<p>When the motion type is &#x201c;Moving,&#x201d; it is seen that the performance ranking of the three frequency bands, from best to worst is: W-band, S-band, C-band. This is the same ranking order for both drone types; however, higher classification accuracies are observed for the big drone compared to small across all frequency bands. The classification accuracies for samples from the big drone were higher than that from samples from the small drone by the following amounts: S-band &#x3d; &#x2b;28.90%, C-band &#x3d; &#x2b;20.60%, W-band &#x3d; &#x2b;6.50%, showing that the smallest performance spread due to different drone types was seen when using W-band. When the motion type is &#x201c;Moving,&#x201d; payload classification accuracy for both drones at all frequencies is lower than for the &#x201c;Hovering&#x201d; motion type, showing that classification of payload weight using signatures from the drone hovering is more accurate than classification using signatures from moving drones.</p>
<p>By taking into account the graphs for both motion types in <xref ref-type="fig" rid="F8">Figure&#x20;8</xref>, it can be seen that W-band data leads to the best classification accuracy in three out of the four cases and S-band data has the best payload classification accuracy for the remaining case. C-band ranks second best when the motion type is &#x201c;Hovering,&#x201d; for both drone types. S-band ranks second best for both drones when the motion type is &#x201c;Moving.&#x201d;</p>
<p>Out of the four cases for which data is shown in <xref ref-type="fig" rid="F8">Figure&#x20;8</xref>, the largest performance spread between any two frequency bands occurs in the case of small drone and a motion type of &#x201c;Moving,&#x201d; between W-band and C-band, where the classification accuracy difference is 47.65%. Across the four cases included in the first investigation, the best overall payload classification accuracy is achieved at W-band when using the big drone in a &#x201c;Hovering&#x201d; motion, while the worst overall payload classification accuracy is achieved at C-band when using the small drone in a &#x201c;Moving&#x201d; motion.</p>
<p>
<xref ref-type="fig" rid="F9">Figure&#x20;9</xref> shows the results from the second classification investigation, in which training sets are constructed using samples from a specific drone type but for both moving and hovering motions, as well as static non-drone target samples.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Payload classification accuracy from data collected at each frequency band for unknown motion but known drone&#x20;type.</p>
</caption>
<graphic xlink:href="frsip-01-781777-g009.tif"/>
</fig>
<p>From <xref ref-type="fig" rid="F9">Figure&#x20;9</xref> it can be seen that the ranking order of the frequency bands from best to worst payload classification accuracy is the same for both drone types: W-band, S-band, C-band. The payload classification accuracy is seen to be higher for the small drone compared to the big drone for S-band and C-band, though this difference is very small (&#x3c;2%). The opposite is seen for results obtained using the W-band data. The classification accuracy spread across the three frequency bands for the big drone is 23.85%, while for the small drone this spread is 19.50%, meaning less variation in performance across frequency bands is seen when using the smaller drone compared to the bigger&#x20;drone.</p>
<p>
<xref ref-type="fig" rid="F10">Figure&#x20;10</xref> shows the results from the third classification investigation where payload classification accuracies are found using mixed drone data but known movements (the left hand set of bars correspond to the hovering movement and the central set of bars correspond to the moving drone movement), as well as for the fourth classification investigation where mixed drone data and mixed movement data is used (shown in the set of bars on the right side of the graph).</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Payload classification accuracy from data collected at each frequency band for unknown drone&#x20;type.</p>
</caption>
<graphic xlink:href="frsip-01-781777-g010.tif"/>
</fig>
<p>From the first two sets of bars (bars for &#x201c;Hovering&#x201d; and &#x201c;Moving&#x201d;) in <xref ref-type="fig" rid="F10">Figure&#x20;10</xref>, it can be determined that when the drone type is unknown, the highest payload classification accuracy is achieved for motion type &#x201c;Hovering&#x201d; for all three frequency bands. The smallest variation in performance across S-bands also occurs for hovering data, where all three frequency bands performed within 4.35% of each other. While the performance differences are small for hovering data, the frequency band ranking from best to worst is seen to be: W-band, C-band, S-band. In the case where the data type is from the drones performing the &#x201c;Moving&#x201d; motion, the frequency band performance ranking is: W-band, S-band, C-band. The performance differences seen between the &#x201c;Hovering&#x201d; and &#x201c;Moving&#x201d; signatures for each band are: S-band &#x3d; &#x2212;17.55%, C-band &#x3d; &#x2212;38.25% and W-band &#x3d; &#x2212;6.90%, showing the smallest variation in performance across different motion types to be achieved by W-band.</p>
<p>When the signatures used in an investigation are of unknown motion and unknown drone (corresponding to the set of bars on the right hand side of <xref ref-type="fig" rid="F10">Figure&#x20;10</xref>), the payload classification accuracy of the frequency bands from best to worst is: W-band, S-band, C-band. It is also seen that S-band and C-band perform better in this investigation than when moving only data is used, but worse than when hovering only data (i.e.,&#x20;the results from investigation three), while W-band shows a minor performance reduction compared to the case where the motion type is known to be &#x201c;Moving&#x201d; (&#x3c;1%).</p>
<p>
<xref ref-type="fig" rid="F11">Figure&#x20;11</xref> includes three confusion matrices pertaining to each of the three frequency bands. Each confusion matrix is constructed by recording the total number of classification attempt outcomes for each investigation carried out within the research presented. The values in the confusion matrices shown the proportional breakdown of outcome class predictions given a particular true payload&#x20;class.</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Confusion matrices from all classification attempts using samples generated from data obtained at each frequency band: <bold>(A)</bold> S-band <bold>(B)</bold> C-band <bold>(C)</bold> W-band.</p>
</caption>
<graphic xlink:href="frsip-01-781777-g011.tif"/>
</fig>
<p>It can be seen from the confusion matrices that the highest proportion of correct classifications for each payload type occurred when using W-band data, followed by S-band data, and the worst performance is seen from C-band data. When using C-band data, the highest misclassification situation occurs when the correct payload is &#x201c;Heavy&#x201d; but the &#x201c;Empty&#x201d; class is predicted.</p>
<p>The second highest misclassification situation when using C-band data occurs jointly when the true payload is &#x201c;Light&#x201d; but the &#x201c;Empty&#x201d; class is predicted, and when the true payload is &#x201c;Empty&#x201d; but the &#x201c;Heavy&#x201d; class is predicted. Three misclassification scenarios occur with the same rate as these two variations when using S-band. For all three frequency bands, no instances were observed where a situation was classed as having no drone present. For all three frequency bands, the most likely true payload to be misclassified is the &#x201c;Heavy&#x201d; payload. Similarly, for all three frequency bands, the most likely true payload to be correctly classified is the &#x2018;Empty&#x2019; payload, though the true class of &#x201c;Light&#x201d; payload has an equal misclassification rate as this in the case of S-band.</p>
<p>All classification tests reported in this section of were repeated using a set of eight conventional SVD features, previously used for micro-Doppler classification in <xref ref-type="bibr" rid="B8">Fioranelli et&#x20;al. (2017)</xref> and <xref ref-type="bibr" rid="B4">Dhulashia et&#x20;al. (2021)</xref>. The four investigations reported include a total of 27 different reported classification accuracy results. A comparison of the classification accuracy achieved in each of the tests using the SVD features with that achieved using the features proposed in this work can be seen in <xref ref-type="fig" rid="F12">Figure&#x20;12</xref>.</p>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>Comparison of classification accuracy achieved in each test using feature set proposed in this paper and using SVD feature&#x20;set.</p>
</caption>
<graphic xlink:href="frsip-01-781777-g012.tif"/>
</fig>
<p>As shown in <xref ref-type="fig" rid="F12">Figure&#x20;12</xref>, the feature set reported in this paper outperformed the SVD feature set in 25 out of the 27 classification tests, with classification accuracies being an average of 26.67% higher when using the feature set proposed in this&#x20;paper.</p>
</sec>
<sec id="s5">
<title>5 Discussion</title>
<p>The results shown within this paper provide evidence that it is possible to classify the payload of different quadcopter micro-drone platforms. The tests included used payloads varying from 0% up to 40% of the drone body weight and achieved a peak classification success of 98.55%. From the results shown in <xref ref-type="fig" rid="F8">Figures 8</xref>, <xref ref-type="fig" rid="F10">10</xref>, it can be seen that higher payload classification accuracies were achieved when using signatures from the &#x201c;Hovering&#x201d; scenarios, compared to those from when the drone is performing the &#x201c;Moving&#x201d; motion. This may be due to a greater number of degrees of freedom in the drone movement for the &#x201c;Moving&#x201d; motion, causing greater variations within the signatures for this motion type. Examples of these include more deviations in: Flight path, drone altitude and drone&#x20;tilt.</p>
<p>It can be determined that the best payload classification accuracy can be achieved using signatures obtained at W-band frequencies. However, the W-band radar system is significantly range limited compared to the C-band and S-band systems, and suffers from rapid SNR reduction as the range to the target increases. For the W-band data, the lack of observable Helicopter Rotation Modulation (HERM) lines within Doppler signatures means that the features extracted were dominated by the bulk Doppler contribution which appears to be sufficiently different in each payload example to allow high classification success rates despite the lack of HERM lines. This may be due to the changed reflectivity of the target from the added metallic weights creating the payload. As such, this conclusion is limited to the short range examples tested in this research and the W-band system is likely to have a significant drop off in classification success rate with increases in range, unless highly directive antennas coupled with significant power amplification is built into the radar system.</p>
<p>In the ideal case, a single system capable of simultaneous measurements of all three frequency bands would be used along with transmission and reception hardware to ensure no difference in gain values. Such a system would allow for a comparison of the classifiability of micro-Doppler signatures for the three frequency bands while ensuring no other operating variables have any bearing on the results. The measured SNR of the S-band and C-band signals collected in this work were seen to be similar, thus legitimizing a comparison of data obtained using different hardware architectures.</p>
<p>In comparing the C-band and S-band results, the S-band was found to generally outperform the C-band measurements, particularly for the &#x201c;Moving&#x201d; and combined movement cases. Therefore, the recommendation based on the findings of this work would be to utilize the S-band frequencies for the task of micro-drone payload classification, due to this frequency band&#x2019;s long range capabilities and high classification success rates. While the superior classification accuracy of the W band is likely due to the greater magnitude of Doppler frequencies spanned by the bulk Doppler signature, the S-band and C-band signature classification results from a greater contribution of the HERM lines observed within the spectrograms. The bulk Doppler contribution will be lower for these bands since this component of the signature spans only a small proportion of the overall Doppler frequency span. The S-band is shown to result in a higher overall classification accuracy than the C-band. Two factors are thought to contribute to this. The first of these results from the differences in the distribution of the HERM lines within the signatures from the two bands. In the S-band spectrograms, it can be observed that the HERM lines are unambiguously captured, while in the C-band spectrograms HERM lines may extend beyond the Doppler frequency range and result in overlap. Varying the payload on the drone results in changes in the rate at which the motor rotates which induces changes in the distribution of the HERM lines. For a higher payload and a higher rate of motor rotation this would appear only as the spreading of HERM line spacing in the S-band spectrograms. Conversely, this change would result in some increase in HERM line spacing and others becoming closer in the C-band spectrograms, due to the overlap wrapping effect. The S-band and C-band performance differences are also partially likely due to the higher SNR of the HERM lines within the S-band spectrograms compared to the C-band spectrograms. Since the features used for classification are based on measuring the distribution of energy within the signature, a more sensitive distinction between signatures of different payload classes is achievable for higher variations between the energy distributions in the spectrograms from different payload classes. The stronger HERM lines observed in the S-band spectrograms would result in higher variations between the energy distributions within segments from different payload classes, thereby resulting in a higher overall classification accuracy compared to the use of C-band spectrograms, which exhibit less energy distribution variations across segments from different payloads due to the lower SNR HERM line contributions. The disadvantage of S-band is that a system operating at this frequency band requires a larger antenna due to the longer wavelengths involved. This may present additional logistical challenges if deployment is required on a lightweight or small platform. Another potential reason for the higher classification success rates achieved when using S-band could be that the features being extracted have a greater propensity towards characteristics more readily observed within the S-band spectrograms compared to the C-band spectrograms.</p>
</sec>
<sec id="s6">
<title>6 Conclusion</title>
<p>This paper has presented results of using micro-Doppler features extracted from radar data obtained at three frequency bands to classify between three different payloads carried by two types of quadcopter micro-drones performing two different motions. It has been shown that the features and classification methodology used in this work can successfully determine the payload being carried by a drone in problems chosen to have varying complexities and situational knowledge. The reported results show that W-band is the preferable choice of radar frequency band to use for micro-drone payload classification. Out of the three frequency bands investigated, S-band is the second most preferable and C-band is the least preferable. It has also been shown that micro-Doppler signatures collected from hovering drones lead to more accurate payload classification than signatures from moving drones.</p>
<p>Further work will aim at collecting additional data in different conditions to validate these preliminary results, including the use of additional models of micro-drones and selecting different locations in which to perform&#x20;tests.</p>
<p>The results shown within this paper are all based on data captured using horizontal polarization (H-pol) geometries. Future work may look to investigate the benefit of a multiple polarized sensing system when looking to classify a drone payload. H-pol was selected to maximize the micro-Doppler signatures but vertical polarization (V-pol) may also give information on flight dynamics, particularly when the drone is moving at higher velocities and travelling at a significantly tilted orientation.</p>
<p>In all four of the classification investigations presented in this work, classification decisions were made solely on the data from a selected frequency band at any given decision instance. If a multi-frequency system were to be used with a voting system to make classification decisions, it is likely that the payload classification accuracy would be greatly increased. This could be achieved by using the data collected from all three frequency bands and applying separate classifiers, before voting on the outcomes to obtain a joint decision. It is expected that this fused decision would provide an improved overall success&#x20;rate.</p>
</sec>
</body>
<back>
<sec id="s7">
<title>Data Availability Statement</title>
<p>The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="s8">
<title>Author Contributions</title>
<p>NP and CH developed and operated the ARESTOR platform in order to capture the S-band and C-band data, wrote the description of the ARESTOR platform and provided inputs on data analysis. DD operated the W-band radar and completed the data analysis and wrote the results and discussion. PB provided the literature review and analysis of the drone data. MR flew the drone, provided inputs on the analysis and wrote the discussion and conclusion.</p>
</sec>
<sec id="s9">
<title>Funding</title>
<p>Funding for this research was provided by: Defence and Security Accelerator (DASA), Defence Science and Technology Laboratory (DSTL), Engineering and Physical Sciences Research Council (EPSRC) grant: EP/R513143/1 and Thales&#x20;UK.</p>
</sec>
<sec sec-type="COI-statement" id="s10">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s11">
<title>Publisher&#x2019;s Note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ack>
<p>The authors would like to thank Sarah and David Tutt for kindly allowing use of their land, on which the experiments performed as part of this work were carried out. The authors would also like to thank Thales UK for their support of this research.</p>
</ack>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="web">
<collab>BBC</collab> (<year>2016</year>). <article-title>Warning over Drones Use by Terrorists</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.bbc.co.uk/news/technology-35280402">https://www.bbc.co.uk/news/technology-35280402</ext-link>
</comment>. [<comment>Accessed January 12, 2019</comment>]. </citation>
</ref>
<ref id="B2">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Beasley</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Ritchie</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Griffiths</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Miceli</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Inggs</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Lewis</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). &#x201c;<article-title>Multistatic Radar Measurements of Uavs at X-Band and L-Band</article-title>,&#x201d; in <conf-name>2020 IEEE Radar Conference (RadarConf20)</conf-name>, <fpage>1</fpage>&#x2013;<lpage>6</lpage>. <pub-id pub-id-type="doi">10.1109/RadarConf2043947.2020.9266444</pub-id> </citation>
</ref>
<ref id="B3">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>V.</given-names>
</name>
</person-group> (<year>2019</year>). <source>The Micro-doppler Effect in Radar</source>. <edition>Second Edition</edition>. <publisher-loc>Norwood, MA</publisher-loc>: <publisher-name>Artech House</publisher-name>. </citation>
</ref>
<ref id="B4">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Dhulashia</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Ritchie</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Vishwakarma</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Chetty</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2021</year>). &#x201c;<article-title>Human Micro-doppler Signature Classification in the Presence of a Selection of Jamming Signals</article-title>,&#x201d; in <conf-name>2021 IEEE Radar Conference (RadarConf21)</conf-name>, <fpage>1</fpage>&#x2013;<lpage>6</lpage>. <pub-id pub-id-type="doi">10.1109/RadarConf2147009.2021.9455278</pub-id> </citation>
</ref>
<ref id="B5">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Doughty</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2008</year>). <source>Development and Performance Evaluation of a Multistatic Radar System</source>. <comment>Ph.D. thesis</comment>. <publisher-loc>London</publisher-loc>: <publisher-name>University College London</publisher-name>. </citation>
</ref>
<ref id="B6">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Dwyer</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>1983</year>). &#x201c;<article-title>Detection of Non-gaussian Signals by Frequency Domain Kurtosis Estimation</article-title>,&#x201d; in <conf-name>ICASSP &#x2019;83. IEEE International Conference on Acoustics, Speech, and Signal Processing</conf-name>, <fpage>607</fpage>&#x2013;<lpage>610</lpage>. <pub-id pub-id-type="doi">10.1109/ICASSP.1983.1172264</pub-id>
<volume>8</volume> </citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fioranelli</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Ritchie</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Griffiths</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Borrion</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Classification of Loaded/unloaded Micro&#x2010;drones Using Multistatic Radar</article-title>. <source>Electron. Lett.</source> <volume>51</volume>, <fpage>1813</fpage>&#x2013;<lpage>1815</lpage>. <pub-id pub-id-type="doi">10.1049/el.2015.3038</pub-id> </citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fioranelli</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Ritchie</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>G&#xfc;rb&#xfc;z</surname>
<given-names>S. Z.</given-names>
</name>
<name>
<surname>Griffiths</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Feature Diversity for Optimized Human Micro-doppler Classification Using Multistatic Radar</article-title>.&#x20;<source>IEEE Trans. Aerosp. Electron. Syst.</source> <volume>53</volume>, <fpage>640</fpage>&#x2013;<lpage>654</lpage>. <pub-id pub-id-type="doi">10.1109/TAES.2017.2651678</pub-id> </citation>
</ref>
<ref id="B9">
<citation citation-type="web">
<collab>Forbes</collab> (<year>2019</year>). <article-title>Warning over Terrorist Attacks Using Drones Given by Eu Security Chief</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.forbes.com/sites/zakdoffman/2019/08/04/europes-security-chief-issues-dire-warning-on-terrorist-threat-from-drones/?sh=10ee551d7ae4">https://www.forbes.com/sites/zakdoffman/2019/08/04/europes-security-chief-issues-dire-warning-on-terrorist-threat-from-drones/?sh&#x3d;10ee551d7ae4</ext-link>
</comment>. [<comment>Accessed January 12, 2019</comment>] <comment>[Dataset].</comment> </citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gurbuz</surname>
<given-names>S. Z.</given-names>
</name>
<name>
<surname>Mahbubur Rahman</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Kurtoglu</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Malaia</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Gurbuz</surname>
<given-names>A. C.</given-names>
</name>
<name>
<surname>Griffin</surname>
<given-names>D. J.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Multi-frequency Rf Sensor Fusion for Word-Level Fluent Asl Recognition</article-title>. <source>IEEE Sensors J.</source> <volume>1</volume>, <fpage>1</fpage>. <pub-id pub-id-type="doi">10.1109/JSEN.2021.3078339</pub-id> </citation>
</ref>
<ref id="B11">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Gusland</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Christiansen</surname>
<given-names>J.&#x20;M.</given-names>
</name>
<name>
<surname>Torvik</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Fioranelli</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Gurbuz</surname>
<given-names>S. Z.</given-names>
</name>
<name>
<surname>Ritchie</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2021</year>). &#x201c;<article-title>Open Radar Initiative: Large Scale Dataset for Benchmarking of Micro-doppler Recognition Algorithms</article-title>,&#x201d; in <conf-name>2021 IEEE Radar Conference (RadarConf21)</conf-name>, <fpage>1</fpage>&#x2013;<lpage>6</lpage>. <pub-id pub-id-type="doi">10.1109/RadarConf2147009.2021.9455239</pub-id> </citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname>
<given-names>B. K.</given-names>
</name>
<name>
<surname>Kang</surname>
<given-names>H.-S.</given-names>
</name>
<name>
<surname>Park</surname>
<given-names>S.-O.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Drone Classification Using Convolutional Neural Networks with Merged Doppler Images</article-title>. <source>IEEE Geosci. Remote Sensing Lett.</source> <volume>14</volume>, <fpage>38</fpage>&#x2013;<lpage>42</lpage>. <pub-id pub-id-type="doi">10.1109/lgrs.2016.2624820</pub-id> </citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Molchanov</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Harmanny</surname>
<given-names>R. I. A.</given-names>
</name>
<name>
<surname>de Wit</surname>
<given-names>J.&#x20;J.&#x20;M.</given-names>
</name>
<name>
<surname>Egiazarian</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Astola</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Classification of Small Uavs and Birds by Micro-doppler Signatures</article-title>. <source>Int. J.&#x20;Microw. Wireless Technol.</source> <volume>6</volume>, <fpage>435</fpage>&#x2013;<lpage>444</lpage>. <pub-id pub-id-type="doi">10.1017/s1759078714000282</pub-id> </citation>
</ref>
<ref id="B14">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Pallotta</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Clemente</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Raddi</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Giunta</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2020</year>). &#x201c;<article-title>A Feature-Based Approach for Loaded/unloaded Drones Classification Exploiting Micro-doppler Signatures</article-title>,&#x201d; in <conf-name>2020 IEEE Radar Conference (RadarConf20)</conf-name>, <fpage>1</fpage>&#x2013;<lpage>6</lpage>. <pub-id pub-id-type="doi">10.1109/radarconf2043947.2020.9266458</pub-id> </citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Patel</surname>
<given-names>J.&#x20;S.</given-names>
</name>
<name>
<surname>Al&#x2010;Ameri</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Fioranelli</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Anderson</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Multi&#x2010;time Frequency Analysis and Classification of a Micro&#x2010;drone Carrying Payloads Using Multistatic Radar</article-title>. <source>J.&#x20;Eng.</source> <volume>2019</volume>, <fpage>7047</fpage>&#x2013;<lpage>7051</lpage>. <pub-id pub-id-type="doi">10.1049/joe.2019.0551</pub-id> </citation>
</ref>
<ref id="B16">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Peters</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Horne</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Ritchie</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2021</year>). &#x201c;<article-title>Arestor: A Multi-Role Rf Sensor Based on the Xilinx Rfsoc</article-title>,&#x201d; in <conf-name>2021 EuRad conference</conf-name>, <conf-loc>London</conf-loc>. </citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rahman</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Robertson</surname>
<given-names>D. A.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Classification of Drones and Birds Using Convolutional Neural Networks Applied to Radar micro&#x2010;Doppler Spectrogram Images</article-title>. <source>IET radar, sonar &#x26; navigation</source> <volume>14</volume>, <fpage>653</fpage>&#x2013;<lpage>661</lpage>. <pub-id pub-id-type="doi">10.1049/iet-rsn.2019.0493</pub-id> </citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rahman</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Robertson</surname>
<given-names>D. A.</given-names>
</name>
<name>
<surname>Govoni</surname>
<given-names>M. A.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Radar Signatures of Drones Equipped with Heavy Payloads and Dynamic Payloads Generating Inertial Forces</article-title>. <source>IEEE Access</source> <volume>8</volume>, <fpage>220542</fpage>&#x2013;<lpage>220556</lpage>. <pub-id pub-id-type="doi">10.1109/ACCESS.2020.3042798</pub-id> </citation>
</ref>
<ref id="B19">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Rahman</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Robertson</surname>
<given-names>D. A.</given-names>
</name>
</person-group> (<year>2019</year>). &#x201c;<article-title>Millimeter-wave Radar Micro-doppler Feature Extraction of Consumer Drones and Birds for Target Discrimination</article-title>,&#x201d; in <source>Radar Sensor Technology XXIII</source>. Editors <person-group person-group-type="editor">
<name>
<surname>Ranney</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Doerry</surname>
<given-names>A.</given-names>
</name>
</person-group>. <publisher-loc>United&#x20;States: SPIE</publisher-loc>. <publisher-name>Proceedings of SPIE</publisher-name>. <pub-id pub-id-type="doi">10.1117/12.2518846</pub-id> </citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rahman</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Robertson</surname>
<given-names>D. A.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Radar Micro-doppler Signatures of Drones and Birds at K-Band and W-Band</article-title>. <source>Sci. Rep.</source> <volume>8</volume>, <fpage>1</fpage>&#x2013;<lpage>11</lpage>. <pub-id pub-id-type="doi">10.1038/s41598-018-35880-9</pub-id> </citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ritchie</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Fioranelli</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Borrion</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Griffiths</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>20172017</year>). <article-title>Multistatic micro&#x2010;Doppler Radar Feature Extraction for Classification of Unloaded/loaded Micro&#x2010;drones</article-title>. <source>IET Radar, Sonar &#x26; Navigation</source> <volume>11</volume> (<issue>1</issue>), <fpage>116</fpage>&#x2013;<lpage>124</lpage>. <pub-id pub-id-type="doi">10.1049/iet-rsn.2016.0063</pub-id> </citation>
</ref>
<ref id="B22">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Ritchie</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Fioranelli</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Griffiths</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Torvik</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2016</year>). &#x201c;<article-title>Monostatic and Bistatic Radar Measurements of Birds and Micro-drone</article-title>,&#x201d; in <conf-name>2016 IEEE Radar Conference (RadarConf)</conf-name>, <fpage>1</fpage>&#x2013;<lpage>5</lpage>. <pub-id pub-id-type="doi">10.1109/RADAR.2016.7485181</pub-id> </citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Taha</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Shoufan</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Machine Learning-Based Drone Detection and Classification: State-Of-The-Art in Research</article-title>. <source>IEEE Access</source> <volume>7</volume>, <fpage>138669</fpage>&#x2013;<lpage>138682</lpage>. <pub-id pub-id-type="doi">10.1109/access.2019.2942944</pub-id> </citation>
</ref>
<ref id="B24">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Tahmoush</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2014</year>). &#x201c;<article-title>Detection of Small Uav Helicopters Using Micro-doppler</article-title>,&#x201d; in <conf-name>Proceedings of SPIE - The International Society for Optical Engineering</conf-name>, <conf-loc>Baltimore, MD</conf-loc> (<publisher-name>SPIE</publisher-name>), <fpage>907717</fpage>&#x2013;<lpage>9077176</lpage>. <pub-id pub-id-type="doi">10.1117/12.2054329</pub-id> </citation>
</ref>
<ref id="B25">
<citation citation-type="book">
<person-group person-group-type="editor">
<name>
<surname>Miceli</surname>
<given-names>W. J.</given-names>
</name>
</person-group> (Editor) (<year>2014</year>). <source>Radar Micro-doppler Signatures: Processing And Applications. Radar, Sonar &#x26; Navigation</source>. <publisher-loc>London</publisher-loc>: <publisher-name>Institution of Engineering and Technology</publisher-name>. </citation>
</ref>
<ref id="B26">
<citation citation-type="web">
<collab>Xilinx</collab> (<year>2021a</year>). <article-title>Zcu111 Evaluation Kit</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.xilinx.com/products/boards-and-kits/zcu111.html">https://www.xilinx.com/products/boards-and-kits/zcu111.html</ext-link>
</comment>. [<comment>Accessed September 1, 2021</comment>]. </citation>
</ref>
<ref id="B27">
<citation citation-type="web">
<collab>Xilinx</collab> (<year>2021b</year>). <article-title>Zynq Ultrascale&#x2b; Rfsoc</article-title>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.xilinx.com/products/silicon-devices/soc/rfsoc.html">https://www.xilinx.com/products/silicon-devices/soc/rfsoc.html</ext-link>
</comment>. [<comment>Accessed September 1, 2021</comment>]. </citation>
</ref>
</ref-list>
</back>
</article>