1. Introduction
Quite often, ground vehicles operate in GNSS-denied environments. In such cases, methodologies are required to enable the localisation of ground vehicles. This paper presents a method to accurately localise the ground vehicles using UAVs. It is assumed that sufficiently accurate GNSS locations of the UAVs are available, and that the UAVs always operate above the tree canopies; for example, localising a UGV that travels under a forest tree canopy can be localised using some UAVs which are flying above the forest tree canopy. In such situations, UAVs can be deployed to collaboratively estimate the location of the UGVs in real time.
The motivations behind developing the research outputs presented in this paper are (1) to develop a UGV localisation method that uses a minimum number of UAVs for UGV localisation; (2) to avoid unobservabilities, which may arise when UAVs are used for UGV localisation; (3) to use the proposed UAV–UGV collaborative system in adverse environment/field conditions. Such UAV–UGV collaborative systems can be used to localise UGVs in battlefields and disaster zones. For example, a UGV can be sent to a high-priority rescue mission where the GNSS reception is weak, and the environment has thick smoke and flames. Moreover, in a rescue scenario during bushfires, such a system can be used to rescue people who are surrounded by bushfires. In such an unfortunate situation, the firetrucks must be operating autonomously, since the firefighters cannot be sent to rescue and the field conditions are adverse. Similarly, UGVs that are used for farming may need the assistance of UAVs for their localisation when their GNSS reception is poor. For example, if a farm is covered by a tall tree canopy, the farming UGVs operating under such a tall tree canopy will not have sufficient GNSS reception. Thereupon, UAVs with sufficient GNSS reception, which hover/fly over the tree canopy, can be used to collaboratively localise farming UGVs operating under the tall tree canopy. Henceforth, a robust localisation method, which has been presented in this paper, is a necessity to navigate such UGVs.
Since there is a motivation to use the proposed UAV–UGV collaborative system in adverse environment/field conditions, sensing methods must be robust against dust, smoke, darkness, high heat, glare, etc. RADAR, LIDAR, Vision Camera, IR camera, and ultrasonic ranging are prospective remote sensing methods. The vision cameras are unable to perform the localisation properly in smoke and also during the night time. Due to the high heat fluxes present, IR cameras will jam the localisation if the localisation is supposed to be performed on a hostile battlefront with frequent glares. Ultrasonic ranging can have accuracy problems due to the impracticality of distinguishing the UGVs clearly from the other sound disturbances. LIDAR seems to show a promising ranging solution, even in vegetation clutter [
1], but the inability to penetrate through the dense smoke makes it unsuitable for sensing during a situation like a forest fire [
2].
However, RADAR technology shows promising results, as it is not affected by adverse environmental conditions such as bad weather [
3] or smoke [
4]. RADAR image processing is very cumbersome, and requires expert human intervention to interpret the RADAR scan images. Ultra Wide Band (UWB) sensing is also a RADAR range-finding technique. Unlike RADAR, it does not have to choose a location where surroundings induce minimal clutter, since UWB signals can sense through clutter [
5]. In UWB sensing, the large bandwidth enhances reliability as the sensing signal contains different frequencies, which increases the possibility that at least a handful of the emitted signal can go through/around obstacles, and the high bandwidth offers improved ranging accuracy [
6]. Moreover, UWB sensors are capable of delivering range measurements at Non-Line-of-Sight (NLOS) situations (e.g., ranging through a forest canopy) without significant degradation of the range measurements [
5,
7]. Often, UWB sensors can be operated in RADAR mode or range sensing mode. Again, if the UWB sensors are used in the RADAR mode, sensing data processing is very complex, e.g., identifying a UGV travelling under a tree canopy using RADAR mode UWB images, which are acquired from a UAV that is flying above the tree canopy. On the contrary, when UWB sensors are used in the range sensing mode, such data processing complexities do not arise, while the sensing robustness is also safeguarded. Therefore, the proposed UAV–UGV collaborative localisation method is a UWB range-only localisation method.
In UWB ranging, Time of Arrival (TOA) techniques are providing less complex, reliable and cost-effective solutions. There are three commonly used TOA techniques, namely: (i) basic two-way ranging TOA, (ii) synchronous two-way TOA, and (iii) asynchronous two-way TOA. The basic two-way TOA expects ideal instrumentation, which results in low accuracy [
6]. In the synchronous two-way TOA method, the time delay in returning back the response to the initial signal sender has been compensated [
8]. In addition to the advantages of the synchronous two-way TOA method, the asynchronous two-way TOA method has compensated for frequency and/or phase mismatches between the UWB transceivers [
9]. Concerning the aforementioned advantages, UWB range sensing is assumed to be performed by an asynchronous two-way TOA ranging algorithm. Most modern UWB sensors that use the asynchronous two-way ranging method have achieved ranging accuracy up to
cm [
7,
10]. Henceforth, the observer model in the localisation algorithm does not have to account for either the time delay in range sensing or the frequency and/or the phase mismatches between the UWB transceivers [
9].
In a practical application of collaboratively localising a UGV using UAVs, reducing the number of UAVs that have to be utilised for the task is equally important. By reducing the number of UAVs, the capital cost that has to be spent can be reduced. Furthermore, the operational costs can also be reduced since the electricity power cost is low when a lesser number of UAVs are to be airborne. Due to the limited flight time of UAVs, additional UAVs are kept by users to run UAV operations without interruption. The additional UAVs are utilised for the operation while the battery swapping is performed. In that regard, if the number of UAVs required for a UAV operation is minimised, the additional number of UAVs that have to be purchased for the application can also be reduced.
However, as the number of drones is reduced to a minimum of two, processing of the range data for UGV localisation faces substantial challenges. The main problem to be addressed is the ambiguity of the localisation due to the loss of system observability. This paper addresses this problem, and shows the successful localisation of ground vehicles using the proposed method. In the literature, the aforementioned problem is known as the “flip ambiguity phenomenon”, and it has been researched in Wireless Sensor Networks (WSN) and in tracking/localisation. In [
11], flip ambiguity has been overcome by using a high number of location anchors/nodes in the WSN so that the ambiguously localised nodes can be identified, and their localisation is supposed to be refined to avoid the flip ambiguity. Since two UAVs are used as anchors in this research, identification of the ambiguous localisation using many UAVs is not possible. In [
12], flip ambiguity in intra-localisation of UAVs in a UAV swarm has been addressed, along with the measurement errors. The solution is based on geometric constraints in a 2D plane like in a WSN, which are based on the range measurement constraints, communication range constraints and kinematic information constraints. In [
13], an Extended Kalman Filter (EKF) has been designed to localise a GNSS unavailable UAV in a UAV swarm, and the flip ambiguity in localisation has been overcome by estimating the angular velocity of the UAV. However, UAVs have low process noise in their motion. Nevertheless, in a noisy process situation, such as in a UGV motion on a farm/forest ground, angular velocity estimations will have significant deviations, so that granting the angular velocity estimation as crisp information to address the flip ambiguity in localisation will not be a reliable solution to a UGV localisation. Therefore, in this research, a constrained state estimation-based method is developed to address the localisation issues arising from the flip ambiguity.
The system observability was analysed in a deterministic approach to identify the unobservable situations in the proposed localisation method. Based on that, a methodology was developed to successfully avoid the localisation errors caused by the unobservability. In this research, constrained stochastic estimation has been used for localisation. The constraints mitigate the challenges that arise when only two UAVs are utilised for localisation. In order to check the ambiguity aversion performance while using the constrained stochastic location estimation, stochastic observability has been analysed during temporary unobservable scenarios using simulations and experiments.
Due to the strictly/narrowly focused operational scenario considered in this research, the authors did not find comparable past research works/methods that possess similar system implementations. Therefore, the authors believe that the system design of this UAV–UGV collaborative localisation is novel. In the constrained state estimation-based localisation method, all of the kinematics-based constraints are newly formulated for localisation ambiguity aversion. Moreover, when the UGV localisation is performed by UAVs in real time, a method has to be developed to validate the efficacy of the CEKF-based localisation method presented. Thereupon, a novel analytical method is formulated to show the efficacy of the CEKF-based localisation method using the Constrained Posterior Cramér–Rao Bound (CPCRB).
The rest of this paper is organised as follows:
Section 2 describes the preliminaries of the motion model, observation model and the EKF-based localisation.
Section 3 describes the unobservability identification method of the EKF-based localisation.
Section 4 and
Section 5 explain the proposed method of overcoming localisation unobservabilities and how the presented unobservability aversion techniques assure the unobservability aversion, respectively. Finally,
Section 6 and
Section 7 present simulation results and experiment results, respectively.
4. Aversion of Ambiguities
Based on the identified singularities in
Section 3, avoiding ambiguity will be of interest in accurately localising the UGVs in real time. Nevertheless, the motion of the UGV and UAVs is subjective to the mission objectives that cannot be altered due to observability issues. Therefore, developing a method which enables a UGV to successfully pass through the ambiguity region in
Figure 2 is the key to eradicating the ambiguity errors in localisation. Thereupon, CEKF has been utilised by imposing constraints based on heuristics, covariance-based accuracy margins and inter-dependency of state variables.
Even if the system model depicts the real system to a greater extent, KFs with constraints [
21,
22,
23,
24] have been used in estimation problems to improve the estimation accuracy by avoiding the unrealistic state estimates irrespective of the system nonlinearities [
21].
In this CEKF-based UAV–UGV collaborative localisation method, the estimation projection method has been used for the state estimation, where the unconstrained
a posteriori estimate of the EKF
is projected into the constrained space and obtaining the constrained estimate
by:
such that
are the linear equality constraints and linear inequality constraints, respectively.
The nonlinear constraints can also be handled in the same way by linearising the nonlinear constraints using the first-order Taylor expansion of the constraints [
21,
25]. In (
13), the matrix
is the weighting matrix. The value of
is equal to the inverse of the
a posteriori covariance matrix of estimation if the projection is based on the maximum probability approach. If the projection is based on the least-squares approach, the matrix
is equal to the identity matrix. In the proposed localisation method, the maximum probability state projection method is followed (
). Due to the maximum probability constraining approach, and also because the EKF is unbiased, the overall CEKF state estimation process is still a minimum variance estimation approach.
Figure 3 illustrates the difference between the two projection approaches intuitively.
In
Figure 3, the least squares method has been followed to calculate
by projecting the estimate
onto the constraint space. In that method, there is no concern about safeguarding the minimum variance objective of the estimator when imposing constraints.
Imposing constraints to an EKF (i.e., CEKF) is a quadratic optimisation problem as formulated in (
13). Active set-based quadratic programming can be utilised to find the constrained estimate
, by identifying the active constraints in each step of the optimisation. Hildreth’s quadratic programming procedure, which is simple and reliable in real-time implementation can be utilised to solve this quadratic optimisation problem [
26,
27,
28].
In the following three sections, the novel constraints, which are the main contributions of this paper, have been formulated with derivations. The symbol depicts the standard deviation of the a posteriori estimation of a variable .
Note: The constraints of the CEKF, explained in this section, are to be imposed only after the EKF’s covariance matrix of estimation is sufficiently converged. All of the constraints are formulated based on uncertainty (i.e., 99.73% confidence level).
4.1. The Position Constraint
Assuming that the UGV motion can be successfully modelled by the CV motion model, a position constraint can be imposed on the UGV state estimation based on the
uncertainty of the position caused by the uncertainty of the velocity. For instance, if the estimated position along
x-direction
and the estimated velocity along
x-direction
of the UGV’s front UWB tag at the previous time step is known, then the position of that tag at this time step can be kinematically anticipated as
. However,
and
are stochastic variables. Therefore, deterministic anticipation is inappropriate, and have to take
to stochastically anticipate the position of that tag. Hence, a compound inequality can be written for that tag as:
and a kinematic constraint can be formulated at the
estimation step based on this compound inequality. In such a kinematically constrained EKF-based estimator,
and
are constrained estimates, and have to be written as
and
, respectively. Then,
will be the free variable for the optimisation in the CEKF, which was explained in
Section 4. Hence, (
14) can be modified as:
The constraints in (
15) can be extended to all the other position variables (
and
) of the UGV tag positions’ state vector (
) by writing another three compound inequalities analogous to (
15). Altogether, all four compound inequalities should be imposed simultaneously in the CEKF, and also, the compound inequalities have to be re-arranged as inequalities such that the consolidated inequality becomes:
where the
diagonal matrices
and
independently have the matrices
and
as their block diagonal matrices, respectively, where
and
is the element-wise square root of a matrix.
Further derivations have to be performed to (
16) in order to implement its constraints in the form of
. First, (
16), must be re-arranged such that (
16) can be obtained as a nonlinear inequality in the form of
, where the
vector consists of all the variables in (
16), which are independent of
. In order to linearise (
16), a first-order Taylor expansion at the
a priori estimate can be used [
25]. Hence, the Jacobian of (
16) can be written as
Using the techniques presented in [
21,
24,
25], (
16) can be written as a linear inequality as
in the form of
.
4.2. The Heading Constraint
Assuming that the UGV motion can be modelled by the CV motion model, a heading constraint can be imposed based on the 3
uncertainty of the heading of a UGV range sensing tag, based on a UGV tag’s velocity triangle. This is a constraint which is based on the kinematics of the UGV on a planar terrain. For simplicity, first, the subsequent derivation was performed by considering a given UWB tag on the UGV as a particle moving on a Cartesian plane. In
Figure 4, the position plane is the Cartesian coordinate frame, relative to which a UGV tag’s positions are denoted. Moreover, the velocity plane is the Cartesian coordinate frame, relative to which a UGV tag’s velocity vectors are denoted. The velocity plane is an instantaneous coordinate frame, the origin of which is placed on the respective UGV tag’s estimated position
at the previous time step. The velocity of a UGV tag is drawn in a solid black arrow, as shown in
Figure 4 on the velocity plane.
When
and
denote the
x-position of a tag,
y-position of a tag, the velocity of a tag along the
x-axis and the velocity of a tag along the
y-axis, respectively; if an assumption is made (for simplicity) that the direction of a UGV tag’s velocity is in the first quadrant of the velocity plane as shown in
Figure 4, then the compound inequality (
19) can be obtained. In the compound inequality (
19), the expression
denotes the direction angle of the relative position vector of the current tag position estimate with respect to the previous tag position estimate. Nevertheless, this direction angle cannot be a deterministic variable due to the uncertainty of the velocity estimates. Referring to the velocity plane in
Figure 4, the maximum and the minimum angle of the estimated velocity vector at the previous time step can be easily identified based on the velocity estimation uncertainties’
edge limits. Hence, the left-hand side expression and the right-hand side expression of the compound inequality (
19) can be formulated. Furthermore, (
19) can be simplified to (
20), where a compound inequality can be obtained, which can be used to impose as a constraint on the UGV state’s velocity variables in the EKF-based localisation.
If the compound inequality constraint in (
20) is generalised such that the constraint can be imposed while the UGV tag’s velocity is in any quadrant of the velocity plane, then a common pattern can be observed. Hence, a consolidated inequality expression for the front UGV tag can be written as
where
and the elements of
are:
If (
21) is extended to both tags of the UGV, a linear inequality in the form of
can be obtained, which can be written as
where
is the equivalent weight matrix in (
21) written for the rear tag of the UGV.
4.3. Node Separation Constraint
Since the UGV is a rigid body and the localisation sensor nodes are rigidly attached to the UGV in a specific separation distance, a heuristic equality constraint can be imposed based on this aspect. If the localisation tag positions’ coordinates are taken from the UGV state, the equality constraint can be formulated as
in the form of
, where
is the sensor node separation distance. In this research, it is assumed that the vehicle length is equal to the sensor node separation. In order to linearise (
23), first-order Taylor expansion at the
a priori estimate can be used, such that
Hence, the linearised equality constraint can be obtained in the form of
as follows:
In order to successfully localise the UGV without undergoing any localisation singularities/ambiguities, the three constraints which are discussed have to be imposed concurrently in the CEKF-based UAV–UGV collaborative localisation, where the final expression encompassing all the constraints will be:
8. Conclusions and Future Works
In a UAV–UGV collaborative localisation with range-only observations, severe localisation errors can occur due to the unavailability of UAVs to be deployed to make sufficient range observations. It is well known that a unique global localisation to track a vehicle on a 2D terrain is possible when three or more drones are deployed to make range-only observations. However, when only two UAVs are deployed, there exists an unobservable region within which it is difficult to localise the UGVs using range-only observations when EKF is used.
The presented CEKF-based localisation method successfully avoids the ambiguity in localisation due to the unobservability. Successful numerical simulations and field experiments show the efficacy in localisation despite the localisation unobservability/ambiguity arising when only two UAVs are deployed for the collaborative localisation task.
The presented CEKF-based localisation method successfully localise a UGV travelling on a 2D terrain, using only two UAVs if an initial position guess (fairly close enough to the UGV) is provided. Nevertheless, the initial position guess must be in the same side where the actual UGV is at; with respect to the unobservability boundary (in
Figure 2). After this initial position guess is provided, the subsequent location estimations will converge to the UGV location. Localising a UGV travelling on a 2D terrain, using less than three UAVs cannot be done due to the estimation singularities/ambiguities which may occur during localisation. While the location estimation is ongoing, the ambiguity aversion features of the presented localisation algorithm will rectify any localisation errors whenever the UGV encounters a localisation ambiguity. Because of that, localising a UGV travelling on a 2D terrain, using less than three UAVs is possible with the CEKF-based localisation method presented.
The stochastic observability analysis, which has been presented in this paper, gives solid evidence of how the CEKF-based UAV–UGV collaborative localisation method can overcome estimation singularities. CEKF-based localisation increases the stochastic observability of the localisation more than the EKF-based localisation method. Furthermore, the kinematic constraints, which take the state variables’ variances, give an assurance of successful localisation despite the localisation singularity almost all the time.
However, there can be quite a few instances where UAVs must be slightly re-positioned to ensure a successful localisation. The CPCRB-based stochastic observability analysis framework, discussed in
Section 5.3 can be used to identify the instances when the UAVs must be slightly re-positioned.
In conclusion, a UGV which is travelling on a horizontal plane can be localised (and the heading can also be estimated) by only two UAVs using the presented CEKF-based range-only localisation method.
In future works, the proposed system can be further researched to analyse localisation efficacy (1) for different range measurement sensitivities, (2) when occasional range sensing occlusions occur, and (3) in localising UGVs travelling in hilly terrains.