1. Introduction
A market of small unmanned aerial vehicles (UAVs)—commonly called drones—is developing very fast [
1] and is stimulated by a very good quality to price ratio—for a relatively small amount a drone can be bought that has excellent performance regarding the possibility of flight and shooting a high-quality video from the air. Unfortunately, the drones could be used not only for commercial, research, surveillance, fun, or any other legal activities but also as a tool of crime. Sometimes intentionally, like in the case of drug smuggling across prison walls or dropping explosive materials [
2]. Another time inappropriate and dangerous use of drones is caused by curiosity or lack of imagination. Camera-equipped UAVs can invade people’s privacy or cause a threat to their life. For example, the report by Bard College’s Center for the Study of the Drone [
3] found that 327 incidents between December 2013 and September 2015 (USA only) posed a "proximity danger" where an unmanned aircraft got within 500 feet (152.4 m) of a plane, helicopter, or other manned aircraft (involved multiengine jet aircraft). The above mentioned, as well as similar accidents, have caused many institutions to start thinking about anti-drone solutions. In the case of objects permanently threatened by drones, for example, airports or buildings crucial for national security, stationary anti-drone systems can be justified. However, there are many cases when using a permanent solution seems to be aimless due to the high cost of the system not corresponding with threat frequency, for example, mass events or meetings with very important persons (VIP). For protecting people and objects in such cases, mobile systems (on vehicle or man-carried) perfectly meet the demands. As an example of such solutions, anti-drone systems Hawk by Hertz Systems Ltd (Poland) [
4] and Radar Backpack Kit by SpotterRF (USA) [
5], are presented in
Figure 1.
Many scientists and engineers involved in the development of anti-drone systems agree that a multi-sensor system could provide a higher probability of drone detection than a single sensor [
6,
7,
8,
9,
10,
11,
12]. Experiments also suggest that one of the most efficient detectors of UAV is a frequency modulated continuous wave (FMCW) radar [
6,
7,
8,
9,
10,
11,
12]. It has a big detection range, transmits low power [
13], and is independent of electromagnetic signals and sounds emitted by drones (UAV operating in a fully autonomous mode and gliding can be detected, too).
Active radars do measurements in a local polar system, where distance and azimuth of the detected object are determined. Three-dimensional radars also provide information about the object elevation. It means that the accuracy of determining the drone localization in a global coordinate system depends not only on the radar performance but also on the accuracy of the radar coordinates and spatial orientation relative to the local horizon and the north. Admittedly, using precise Global Navigation Satellite System (GNSS) positioning, for example, Real Time Kinematic (RTK), can solve the problem of accurate localization. Orientation angles can also be relatively easily determined in the case of a stationary system. However, there would be some difficulties in mobile solutions, especially if the system changes location during the mission. Obviously, some sophisticated techniques, like Moving Base RTK, can provide precise orientation angles; however, using them is not always possible. They also do not ensure expected accuracy in every condition and require additional devices and often specialized knowledge. Although a magnetic compass could be a low-cost and convenient solution, it is not accurate enough if we consider that accuracy better than is required. Due to magnetic variation, magnetic deviation, errors connected with radar and magnetic sensor installation on a vehicle, and method of measurement, accuracy better than is difficult to achieve in a real application.
Long-range FMCW radars used in non-military anti-drone systems have detection ranges of micro-drones (DJI Phantom sizes) not exceeding 3 km, due to the small radar cross section (RCS) of such objects [
14,
15,
16,
17]. However, the classification range has a practical meaning, as anti-drone radars also detect birds and other small objects. Determining if the object is a drone allows relevant actions to be taken. The most common method of classification is the micro-Doppler signature (MDS) analysis. It allows investigating micromotions of the object detected by the radar. The contributions of the rotating turbine or propeller blades in radar backscattering in the form of micro-Doppler contents contain information that there are rotating parts on the object. On this basis, the method distinguishes a drone from a bird [
18,
19,
20,
21,
22,
23]. If the analysis of MDS is a classifier, the classification range is about 50% of the detection range, due to the RCS of the rotating turbine or propeller blades being smaller than the RCS of the drone body. The classification range is the range in which the classifier works accurately. Examples of detection and classification characteristics of the DJI Phantom 4 drone for one of the anti-drone radars are presented in
Figure 2. It should be noted that it is currently one of the radars with the longest range of the micro-drone detection and classification available on the commercial market.
Considering the real range of non-military FMCW radars and the fact that effective protection of the area or object requires detection and classification of the drone several hundred meters before the forbidden zone, it is often necessary to use several mobile systems operating in the network. One of the key conditions for their proper cooperation is the measurement of azimuths relative to the same direction (preferably geographical north, also called the true north, or just, the north). Otherwise, errors will occur in determining the actual location of the drone. In a single-radar system, a large error in determining the coordinates of the detected object may result in incorrect classification as an object “violating” or “not violating” the forbidden area (see
Figure 3) and, consequently, an inadequate response. In network systems, where several radars operate in the same area and their coverages partly overlap, there will be additional uncertainty concerning a location of the detected drone—each of the radars will indicate different coordinates of the same object (see
Figure 4). Moreover, one will not know how many objects are really in the area. A single drone can be classified as a swarm or vice versa. All described cases are very undesirable, and their elimination is necessary to ensure the proper functioning of the anti-drone system.
Due to the fact that non-military anti-drone systems are still developing and many of them are on a low Technology Readiness Level (TRL) - similarly as some of the radar solutions in this area [
3,
24], the market usually offers single radars or radar systems integrated on one platform or assembly set. It ensures that the measurements are in the same local coordinate system (not necessarily oriented to the north). Then, the problem of azimuths misalignment does not exist due to all the devices making the same constant error in azimuths measurements. However, if single radars are integrated into a bigger system, where sensors are distributed over a wide area, bringing all the measurements to a common coordinates system is required. It means that a misalignment reduction of azimuths is necessary. In other cases, after the transformation of local coordinates of the detected object to the global one, different results will be obtained, depending on the azimuth error of the particular radar (see
Figure 4).
The co-author of the paper has participated in the research project granted by the Polish National Center for Research and Development since 2017. One of its aims is to develop an anti-drone system consisting of many mobile sensors (on the vehicles) spread over a wide area. Experiments already performed in a real environment show the problem of azimuth misalignment of particular sensors. It causes improper operation of the whole system and repeats every time the vehicles change location. As an example, the PrintScreen from the developing system console is shown in
Figure 4.
In
Figure 4, it can be seen that one object is identified as three by the system due to azimuth misalignment. The radars one, two, and three made constant azimuths errors of
,
,
, respectively, relative to the north. The errors were determined using GNSS RTK on a drone and the radars.
The azimuth misalignment of
produces the drone coordinate error of about 87.5 m at a distance of 1 km (see
Figure 5). It means that two radars will determine the coordinates of the same object differing by 87.5 m. As already mentioned, in mobile systems (especially on vehicles), it is difficult to obtain the radar azimuths alignment better than
. Then, the coordinates error will be 364 m (see
Figure 5) at a distance of 1 km. Doubtlessly, this will disturb the proper operation of the system. Moreover, variable errors of azimuth measurement of both radars will cause a further increase in the drone position uncertainty. This is a quite serious issue and must be taken into account.
As was mentioned earlier, regarding a permanent system, the azimuths misalignment can be reduced using, for example, geodetic methods. However, radars in a mobile system can change location during the mission. Therefore, in new places, they must be aligned to the rest of the system as quickly as possible. To the best knowledge of the authors, this issue has not been described in publications yet (probably due to the commercial multi-radars mobile systems are still on low TRL). This statement is based on an unsuccessful search for a suitable method description. ‘Suitable’ means giving the desired accuracy at an acceptable time. However, the authors are conscious that this statement may be a kind of abuse. Research done by the authors shows that, at the moment, in practical applications, each of the radars is only roughly set separately using magnetic sensors and more advanced devices, such as gyrocompasses, GNSS RTK, or satellite compasses, are not applied in non-military mobile systems. Such a method does not ensure the desired accuracy.
This work presents a proposal to solve the problem of azimuths misalignment, consisting of the initial/coarse orientation of the radar system to the north using a magnetic compass, and then carrying out a calibration procedure consisting of automatically calculating corrections to azimuths measured by individual radars. The calibration process uses the "friend" drone flying in autonomous mode along a fixed route. The position of the drone does not have to be known, and the drone does not need to communicate with the system. The proposed method is simple, automatic, and requires no specialized knowledge. It gives very good results and can be successfully applied in real applications.
2. Proposed Method of Reducing the Azimuths Misalignment
In an ideal situation (no measurement errors), the coordinates of the drone detected by the radar are expressed by the following system of equations:
where:
| - true coordinates of the drone in the ENU system (East North Up), |
| - true coordinates of the radar in the ENU system, |
| - horizontal distance (2D) to the drone, |
| - azimuth of the drone from the radar position. |
However, in real conditions, all quantities occurring on the right side of equations in the system (1) are affected by errors. Therefore, (1) will take the form
where:
| - radar coordinates determination error, |
| - horizontal distance measurement error, |
| - azimuth measurement error in the ENU system. |
Taking the practical experience in using FMCW radars for drone detection into consideration, it should be noted:
Note 1.: In practice, the radar position is fixed by GNSS before the operation and taken as a constant until the radar is moved to a new location. In such a case, is a constant error. It can be minimized by using Differential GNSS or GNSS RTK measurement techniques and by mounting the satellite system antenna centrally above the center of the radar measuring system. As already mentioned, the GNSS RTK technique is not used in commercial systems. A convenient solution is a Satellite Based Augmentation System (SBAS)—GNSS augmented by the differential system based on satellites. This reduces the radar position error to a sub-meter level. It allows omitting it in further considerations due to a small impact on the final performance of the proposed method.
Note 2.: is a random variable with some small constant component resulting from the leveling of the mobile platform or mounting system of the radar. The constant component can be reduced almost to zero by using a leveling system and then can be modeled as an independent random variable with a normal standard distribution. Then, the highest probability density is concentrated around zero. Therefore, is also omitted in the proposed method.
Note 3.: is a random variable with a small variable component and a very large constant component, resulting from magnetic declination and magnetic deviation, radar and compass assembly on the platform, and the physical possibilities of setting the mobile platform or mounting system of the radar in the direction indicated by the magnetic compass. In practical applications, the variable component varies within
(depending on the class of anti-drone radar), while the constant component can reach
. Such a large error in measuring the azimuth of the detected object will result in an error in determining its location (coordinates). According to Equation (2), it will be larger the greater the distance between the radar and the detected object is (see
Figure 5). This will entail some complications, both in single- and multi-radar systems, which have already been mentioned in the Introduction.
Bearing in mind Notes 1 to 3, and that
the system of Equation (1) can be written as
This system has 3 unknowns: The estimated drone position in the ENU system
and the azimuth measurement error in the ENU system
. The radar measurement results are the horizontal distance
and the azimuth
, where
and:
| - the spatial distance (3D) to the detected drone measured by the radar, |
| - the drone elevation measured by the radar or a priori assumed. |
Estimated radar coordinates (fixed using DGNSS/SBAS measurement) in the ENU system are also known. On the left side of the equations, instead of true coordinates, there are their estimators, because both the position of the radar and the measurements are burdened with the errors. The errors were omitted in the system (4); however, they have an impact and cause some error in determining the coordinates of the detected drone.
From (4), it follows that unambiguous determination of unknowns based on measurements from a single radar is not possible. The next radar measurement will introduce the next 2 unknowns (new 2D coordinates of the drone), and 2 independent equations. Therefore, regardless of the number of measurements done by a single radar, the number of independent equations will always be 1 less than the number of unknowns.
Let us now consider a situation where a few radars in the network are tracking the same object. Then, the following system of independent equations can be formed
where:
| - the estimators of the drone true coordinates calculated on the base of measurements done by i-th radar, |
| - the estimators of the i-th radar true coordinates, |
| - measurements results of horizontal distance and azimuth done by i-th radar, |
| - number of radars in the network. |
Let us mark the position of the detected drone in the ENU system as
. This is the position estimated on the base of measurements done by
i-th radar. Since both estimators of the radar coordinates and radar measurements are affected by errors, then
In this case, system (6) also has no unambiguous solution, because the number of unknowns will be
, and the number of independent equations will be
. However, due to all radars are observing the same object, we can assume a simplification:
and then, system (6) will take the following form:
In this situation, the number of unknowns will be
, and the number of independent equations
, and already for
, all 3 unknowns can be fixed. However, for
, redundant equations will appear. In such a case, we have
numbers of possible equations combinations, which give a solution in the form of
. Unfortunately, due to errors both in the estimated coordinates of the radar position and in the radar measurements, each of these combinations will give a slightly different solution. These differences will depend not only on the radar coordinates errors and the azimuth, elevation, and distance measurement errors but also on the relative location of the radars and the detected object (the system geometry described by the dilution of precision coefficient). It should be noted that the geometry of the system will change during the drone’s flight. Therefore
calculated for the moment
, based on observations from
radars, can be expressed by the following equation:
where:
| - constant component of , |
| - variable in time component of , |
| - variable in time component of resulting from the changing geometry of the radars-detected object arrangement due to drone flight. |
is the sought difference between the radar north and the true north:
Due to
, the determination of
would significantly improve both the accuracy of determining the coordinates of the detected object and minimize the likelihood of classifying a single drone as a swarm and vice versa. Both factors are crucial for the anti-drone system. The calculated
value with the opposite sign would be a correction to the measured azimuth by the radar and could be automatically taken into account by the master system responsible for the multi-sensors data fusion. In the proposed method, determining the difference
, and hence the correction for measuring azimuths consists of eliminating the
component from Equation (10). According to the proposed idea, this is done by multi-radars observations of a flying drone with unknown coordinates. Azimuth measurement error
is estimated in every measurement epoch
, and then the constant component of the azimuth measurement error is calculated as:
where:
| - number of measurement epochs common to all radars, |
| - azimuth measurement error at the moment . |
Therefore, the key issue is to estimate the azimuth measurement error
in every measurement epoch. In the proposed method, instead of the classical form of system of equations describing the coordinates of the detected object (9), the following form is used:
Let us assume the approximate drone position and mark it as
. Let it be a position calculated on the basis of observations from any one radar from the network:
Now let us expand the right sides of the equations of the system (13) in the Taylor series relative to
, limiting to the second order components:
wherein:
Let’s mark the individual derivatives as follows:
Then, the system of Equation (15) will take the form
The system of Equation (18) can be written in matrix form as
is a vector of unknowns. Estimated value of
using least square estimator is:
The solution of (20) is:
Using
, the estimated position of the drone can be calculated as:
If
or
, the next iteration is necessary, taking the estimated drone position calculated from (22) as the new approximate position. The final value of the azimuth measurement error for each of the radars is the value of
from the last iteration. Then, the
value is calculated based on (12). According to the presented idea the correction to measuring azimuths (the difference between the true north and the radar north) for a given radar is:
Then, the azimuth estimator of the detected object for a given radar is:
The practical implementation of the proposed method can be presented in the form of the following algorithm:
BEGIN //drone started a calibration flight and all radars detected a drone |
j := 0; //first measurements epoch |
n := number of radars; |
REPEAT |
FOR i := 1 TO n DO |
Get the radars measurements ; |
Bring to a common moment ; //explanation below the algorithm |
Compute for all radars; //using Equation (20) |
j := j + 1; //next measurement epoch |
UNTIL drone ended a calibration flight |
m := j; //number of measurements epoch |
Computation of ; //using Equation (12) |
END. |
Since observations of the drone by each of the radars are generally not done at the same time, but in some time window (approximately equal to the period of scanning of a particular radar), the values of measured azimuths and distances must be brought to a common moment
. The proposed method uses linear interpolation between the measurements immediately preceding and following
, according to the formulas:
where:
| - the time immediately preceding , when the measurement was done, |
| - the time immediately following , when the measurement was done, |
| - estimated values of azimuth and distance at moment , |
, | - measured azimuths at moments i , |
, | - measured distances at moments i . |
4. Results
According to the assumptions, 95,000 simulations were performed for each of the two variants of the radar location (the 3- and 4-radars variants).
Table 1 and
Table 2 show the calibration accuracy obtained due to the proposed method. The measure used is the difference between the radar north and the true north after calibration.
As it results from the data presented in
Table 1 and
Table 2, the average value of the constant component of the azimuth error after calibration is close to 0 for each of the radars in the network, both in the case of the system consisting of 3 and 4 radars. In the best cases, it has been eliminated. However, in the worst cases, after calibration, it was still around
. Therefore, to better analyze the effectiveness of the proposed method,
Figure 8 presents histograms of the azimuth error constant component for radars R1, R2, and R3 in the 3-radars network (see
Figure 7 (left)) before and
Figure 9 after calibration. On the horizontal axis of the charts, there are ranges of the azimuth error constant component in degrees, and the numbers in the blue rectangles are the numbers of individual ranges.
Figure 8 and
Figure 9 show that for each of the radars, there was a clear improvement—a decrease in the value of the constant component of the azimuth error. Moreover, before calibration, the errors had a uniform distribution (according to the simulation assumptions), while after calibration, it was similar to the Gaussian distribution, with the highest probability density centered around 0. It is worth noting that only about 5% percent of errors after calibration exceeded the value of
, which should be considered a very good result, given that before calibration, such errors were about 65%. It is also apparent from
Figure 8 and
Figure 9 that the method impacts on a constant azimuth measurement error for each of the radars in a similar way. The same is also true in the 4-radars network (shown in
Figure 7 (right)). Histograms of the azimuth error constant component for radars R1, R2, R3, and R4 before calibration are shown in
Figure 10 and after calibration in
Figure 11.
It can be seen in
Figure 8,
Figure 9,
Figure 10 and
Figure 11 that the azimuth error constant component before and after calibration for both 3- and 4-radars networks have similar distributions. Therefore, the results of further analyses are presented based on one selected radar from both tested cases. The solution was applied to limit the volume of the paper. To maintain proper reliability, the results of the analyses are presented for the worst case (least improvement and greatest deterioration after calibration). This was the case of R1 radar from the 4-radars network (see
Figure 7 (right)).
The graphs in
Figure 8,
Figure 9,
Figure 10 and
Figure 11 show that the proposed method significantly improves the accuracy of the azimuth measurement by reducing azimuths misalignment (decreasing the constant error of the azimuths measurements). However, the distribution of the azimuth error constant component after calibration does not yet give a complete picture of the method’s effectiveness. There are no answers to the following important questions:
- a)
Is there radar orientation improvement in each of the cases, or are there cases where the value of the azimuth error constant component increases after calibration?
- b)
If there are cases of radar orientation deterioration relative to the north after calibration, how large are they, and do they preclude the practical application of the proposed method?
- c)
What is the percentage of the radar orientation relative to the north improvement, in comparison to the initial value of the azimuth error constant component?
To answer the above questions, an analysis of the radar north improvement was performed. The measure used was the value by which the difference between the radar north and the true north changed due to calibration. To better show the effect of the proposed method, the radar north improvement was expressed as a percentage of the constant azimuth error before calibration. The results of the analysis for the selected radar (the worst case) are presented in
Figure 12,
Figure 13 and
Figure 14.
Figure 12 shows the percentage histogram of improvement and deterioration of the selected radar orientation relative to the north. The horizontal axis presents the percentages, and the numbers in the blue rectangles are the number of individual ranges. A negative percentage value (
) means that after calibration, the radar orientation has deteriorated.
Figure 12 shows that up to 11.6% of cases, after calibration by the proposed method, the radar orientation deteriorated relative to the true north. Most values of degradation did not exceed 100% of the pre-calibration value, but there were also cases where they reached 200%. At first glance, this seemed very worrying, because even a 20% deterioration may be unacceptable if it concerns large initial values (before calibration). Therefore, it seems to be important to answer the following questions:
- a)
How big is the deterioration in numerical values?
- b)
Does it concern small or large values of initial orientation errors?
- c)
How this might affect the practical application of the proposed method?
To answer them, deterioration cases were thoroughly analyzed.
Figure 13 shows the histogram of the initial (before calibration) constant azimuth errors for the R1 radar in the 4-radars system, where the orientation was deteriorated after the calibration.
Figure 13 shows that the large percentage deterioration of radar orientation relative to the north was related to small initial values (before calibration). The deterioration by more than 50% basically concerned angles not greater than
, and 100 not greater than
. This means that even in the worst case of deterioration, the azimuth misalignment after calibration did not exceed
, which is completely acceptable from a practical point of view. Moreover, it should be emphasized that in all analyzed cases, the deterioration of radar orientation relative to true north after calibration, concerned only one radar in the network, while significantly improving the orientation of the other radars. The deterioration effect occurred when one radar had a small initial error and the rest had large ones.
To show the advantages of the proposed method for the same radar, the results of an analogous analysis regarding the improvement of orientation relative to true north are presented below.
Figure 14 shows the distribution of the initial (before calibration) azimuth constant errors for which the R1 radar orientation improved.
The graphs shown in
Figure 14 show the advantages of the proposed calibration method. It can be seen that the method improves the initial orientation of the radar relative to the north over the entire range of initial errors. It is the most effective for errors above
.