HPO-RRT*: a sampling-based algorithm for UAV real-time path planning in a dynamic environment

The real-time path planning of unmanned aerial vehicles (UAVs) in dynamic environments with moving threats is a difficult problem. To solve this problem, this paper proposes a time-based rapidly exploring random tree (time-based RRT*) algorithm, called the hierarchical rapidly exploring random tree algorithm based on potential function lazy planning and low-cost optimization (HPO-RRT*). The HPO-RRT* algorithm can guarantee path homotopy optimality and high planning efficiency. This algorithm uses a hierarchical architecture comprising a UAV perception system, path planner, and path optimizer. After the UAV perception system predicts moving threats and updates world information, the path planner obtains the heuristic path. First, the path planner uses the bias sampling method based on the artificial potential field function proposed in this paper to guide sampling to improve the efficiency and quality of sampling. Then, the tree is efficiently extended by the improved time-based lazy collision checking RRT* algorithm to obtain the heuristic path. Finally, a low-cost path optimizer quickly optimizes the heuristic path directly to optimize the path while avoiding additional calculations. Simulation results show that the proposed algorithm outperforms the three existing advanced algorithms in terms of addressing the real-time path-planning problem of UAVs in a dynamic environment.


Introduction
Recently, unmanned aerial vehicles (UAVs) have gradually played an increasingly critical role in civil fields (such as earthquake relief and vegetation protection) and military fields (such as urban penetration and attacking enemy targets) [1].Path planning, as one of the most fundamental and important problems in autonomous UAV flight, has been widely studied [2].In recent decades, many path-planning algorithms have been proposed, which can be roughly divided into nature-inspired methods, grid-based methods, and samplingbased methods.Nature-inspired algorithms include artificial potential field (APF) methods [3], the interfered fluid dynamical system (IFDS) method [4] and intelligent algorithms.Among them, the APF and IFDS have the advantages of simple calculation and strong real-time performance, but the planned path can easily fall into a local minimum.Intelligent algorithms such as a genetic algorithm [5] can obtain the global optimal solution by iteration, but the calculation is more complicated in high-dimensional spaces.Heuristic search methods include the A* algorithm [6], anytime repairing A* (ARA) algorithm [7], sparse A* algorithm [8] and D* algorithm [9].Although such algorithms can ensure that the optimal solution can be found if it exists, their search time increases sharply with increasing problem scale and spatial dimension.Sampling-based methods such as rapidly exploring random trees (RRT) [10] and probabilistic roadmaps (PRM) [11] usually sacrifice the optimal solution at grid search resolution for the ability to quickly find satisfactory solutions in high-dimensional complex state space and largescale problems.
This paper focuses on solving the real-time path-planning and optimization problem of UAVs in complex dynamic environments with moving threats.Because of the influence of moving threats, planning needs to be constantly revised following the update of environmental information [12].Therefore, the path-planning algorithm must have high planning and optimization efficiency and low storage requirements.Efficient planning and optimization calculations ensure that a UAV can generate the optimal path suitable for flight in real time, and low storage can effectively improve the scale of the UAV planning space.Because RRT and its derived algorithms [13,14] build a tree structure through random sampling to quickly expand the search space, these algorithms still have very low computational complexity and high planning efficiency in high-dimensional environments.Therefore, these algorithms have good prospects in addressing UAV path planning in complex environments with moving threats.Note that most of the early works on RRT algorithms mainly focus on finding the feasible solutions, leaving a great possibility for improving the execution time and path optimization [15].
To deal with the path-planning problem in dynamic environments with moving threats, many improved RRT algorithms have been proposed, which are based on two different ideas.One idea is that algorithms will update the environment information and replan the path immediately after detecting changes in the environment, which requires the algorithms to have faster response speed, higher planning and optimization efficiency.Another idea is that algorithms actively avoid threats by predicting the trajectory of dynamic threats and planning collision-free paths.The algorithms with this idea can reduce the update frequency so that more time can be used to optimize the planning path to deal with the dynamic environment.Recently, many scholars have studied these two categories of RRT algorithms for dynamic environment path planning.The DRRT [16] algorithm uses tree pruning to update the plan.It takes the target position as the root of the tree to extend the tree, which simplifies the extension process.Bryant proposed an online RRT* algorithm [17] that ensures that new targets are always added to the spanning tree through online replanning.RRT X [18] uses a fast rewiring operation to repair damaged branches to rebuild the tree structure after detecting changes in the environment.These three algorithms mainly consider the information of moving threats from the current time to the next step of planning, meaning that path replanning is needed whenever the environment changes.Although the above algorithms optimize the replanning process to improve computing efficiency and ensure real-time performance, they still need a high replanning update frequency to address the movement of threats.In contrast to the above three replanning algorithms, the time constraint is considered by an algorithm based on the partial path-planning method [19].This kind of algorithm adds the nodes without collision with the known threat trajectory to the tree and parallels the path planning by executing the UAV.Influenced by the partial path-planning method, a time-based RRT algorithm named Risk-RRT is proposed in [20].Risk-RRT predicts the movement of a threat through a Gaussian process.Then, RRT is used to plan part of the path and control the robot to advance a certain distance along the planned path.Note that because Risk-RRT uses the traditional RRT as the underlying planner, the algorithm can quickly obtain feasible paths even though it cannot obtain optimized paths.Therefore, to improve the quality of the planned path, Zhang proposed the Risk-RRT* algorithm [21].Risk-RRT* improves the underlying planner-based RRT to RRT*.The rewiring process of RRT* can optimize the tree structure to improve the planning path optimality.However, the rewiring process of RRT* will change the tree structure.Since Risk-RRT* is a time-based RRT* algorithm, the change in tree structure means that each affected node in the tree needs to update the information, including the parent and child node, timestamp and depth of the node.This process often brings many additional calculations, which reduces the percentage of time devoted to the path optimization process within the limited planning time and then affects the quality of the planned path.
This analysis shows that the above algorithms still have space and a possibility to be improved.
(1) The flight environment of UAVs is usually a large threedimensional (3D) space with many moving threats, and the moving frequency of threats is also high.Adopting replanning algorithms requires maintaining a high planning update frequency, which is difficult to achieve.(2) In a large 3D space, if the random sampling method is used to search the whole planning space, many useless sampling nodes may be generated.This result will make the planning inefficient and increase the planning time.
(3) The longer a UAV performs a mission in a battlefield environment, the greater is the probability of being detected and attacked by the enemy.This circumstance necessitates improving the efficiency and quality of the planned path so that the UAV can quickly complete its task on the battlefield and reduce this probability.Therefore, combined with the above motivations, this paper will focus on further improving the partial planning algorithm based on time RRT.The algorithm can improve the efficiency of feasible path planning and quickly optimize the path based on avoiding additional computational consumption to ensure rapid planning efficiency and high path quality in the online dynamic path planning of UAVs with moving threats.
Therefore, in this paper, we propose the HPO-RRT* algorithm to solve the real-time path-planning problem of UAVs in a 3D dynamic environment.Figure 1 shows the basic framework of the algorithm, which adopts a hierarchical architecture.The HPO-RRT* algorithm can quickly plan a path suitable for UAV flight and ensure the optimization of the planned path.First, the HPO-RRT* algorithm uses the UAV perception system to update the world information and predict the moving threats before each plan.Then, an improved time-based RRT* algorithm is used in the path planner to obtain the heuristic path: the sampling bias method based on the APF function is used to guide the sampling towards the target position and avoid threats, the process of selecting the optimal parent node for the new node in RRT* is introduced, and lazy collision checking is adopted to reduce the running time of the algorithm while satisfying the avoidance of dynamic threats.Thus, a partial heuristic path can be quickly obtained.Finally, a low-cost path optimizer is used to quickly optimize the path generated by the path planner to ensure the homotopy optimality of the final planned path.
Our work mainly includes the following innovations: (1) A hierarchical planning framework suitable for realtime fast path planning in a UAV dynamic environment is designed (Fig. 1), which comprises a UAV perception system, a path planner based on the improved timebased RRT* and a low-cost path optimizer, to improve the planning efficiency and ensure path optimization.(2) In the bottom planner, a bias sampling method based on the APF function is used to guide the generation of sampling nodes, which can improve the sampling efficiency and reduce the planning time.(3) In the path planner, the process of neighbourhood optimal parent node selection and lazy collision checking is applied to the tree expansion, which improves the planning efficiency and obtains a path optimized as much as possible.(4) A low-cost path optimizer is designed to quickly optimize the path generated by the path planner without changing the tree structure.(5) The probability completeness and homotopy optimality of the HPO-RRT* algorithm are discussed and proven.
This paper is organized as follows.The following section describes the related work."Problem formulation" describes the problems related to path planning in a dynamic environment."HPO-RRT* algorithm" introduces the HPO-RRT* algorithm in detail in three aspects: the overall framework of HPO-RRT* ("Overall framework"), tree structure optimization and heuristic path planning ("Path planner"), and the low-cost heuristic path optimization method ("Low-cost path optimizer").The analysis and proof of the probability completeness and homotopy optimality of HPO-RRT* are introduced in "Analysis".In "The results of simulation experiments", the simulation experiments are reported.In "Conclusions and future work", we draw conclusions and discuss future work.

Related works
In this section, we first introduce the sampling strategy and related works to improve the RRT convergence speed.Then, Risk-RRT is taken as a typical example to describe the relevant basic knowledge of time-based RRT.

Sampling strategy
As a single-query path-planning method based on sampling, the RRT algorithm has been widely studied in the past three decades because of its high computational efficiency and strong scalability.However, RRT has the problem of slow convergence.Therefore, how to quickly converge has become a research focus.
The results show that changing the sampling strategy can effectively increase the convergence speed of the algorithm.Most RRT-like algorithms use an unbiased and uniform sampling strategy to explore the configuration space, which is pointless and greatly increases computational consumption, resulting in slower convergence speed.Therefore, an importance sampling strategy is proposed, which can extract samples from the preset distribution or configuration space.This method effectively increases the convergence speed [22].Using prior information to guide local sampling is a strategy based on importance sampling.This strategy uses the planned path as a priori information to guide sampling around it, such as informed RRT* [23] and its batch sample iteration improved algorithm BIT* [24].In addition, the ellipsoid heuristic sampling strategy is also adopted in [25,26] to increase the convergence speed.Compared with the prior information-guided sampling strategy, direct bias sampling avoids the uniform sampling process when exploring feasible paths and usually brings a higher convergence speed.This strategy biases random samples to favorable areas through specific methods, such as Theta*-RRT*, to explore the configuration space.Using the APF method as a bias sampling strategy is a good choice.The attractive potential field of the APF is used to guide random sampling towards the target position bias in [27,28].Pharpatara [29] introduced a rotating potential field to the traditional APF to guide offset sampling with direction information.This paper also adopts the offset sampling method based on the APF function to 123 Fig. 1 Basic framework of HPO-RRT* increase the convergence speed.In contrast to [28], we combine attractive and repulsive potential fields to guide random sampling to further improve the sampling efficiency.

Time-based RRT
The difference between time-based RRT and classical RRT mainly lies in the representation of nodes in the tree.The tree node of classical RRT only records the location information of the node, while node x of time-based RRT can be defined as a structure containing six kinds of information, i.e. x : x coor , x parent , x sons , n, t, P collisioncheck (t) .
(1) x coor represents the position information of x and is represented in O − x yz.(2) x parent represents the parent node of x.Note that x parent is also represented by a structure like x. (3)x sons represents any child node of x, similar to x parent .(4) n indicates the depth of x, that is, it needs to pass through n nodes from the root node to the x.(5) t represents the timestamp of x, which can be calculated by depth, i.e. t = t 0 + n t, where t 0 is the timestamp of the root node and t is the time increment required for the UAV to pass through the adjacent nodes.(6) P collisioncheck (t) is the collision risk probability between the node and the obstacles in the environment at time t.In the Risk-RRT algorithm, P collisioncheck (t) is defined as [20]: where P static represents the collision risk probability of the static threat, P moving (t) represents the collision risk probability obtained for the moving threat at time t, and P moving ( p k (t)) represents the collision risk probability between the k-th moving obstacle and the node at time t.Time-based RRT is also an incremental algorithm similar to RRT, i.e. the tree structure continues to expand with increasing sampling nodes.In contrast to RRT, in the tree structure of time-based RRT, tree nodes with high collision risk or timestamps smaller than the current root node will be deleted.

Problem formulation
In this section, we define the path-planning problems in a dynamic environment with moving threats that will be solved in this paper and the notation used to describe them.
The planning environment is abstracted as a configuration space X ⊆ R d .Let X obs ⊆ X be the threat space composed of static threat X sta and moving threat X mov , i.e.X sta , X mov ∈ X obs .The obstacle-free space is defined as X f ree = X \X obs .In a dynamic environment, the configuration space is time varying, so the above two spaces are represented as X obs (t) and X f ree (t).x init , X goal are the initial position and the target area, respectively.In this paper, X goal = x ∈ X f ree |ρ x, x goal ≤ r goal is defined in the dynamic environment to detect whether the goal is reached, where ρ(x 1 , x 2 ) represents the Euclidean distance between x 1 and x 2 .A ball with centre x ∈ X and radius r ∈ R >0 is denoted as B x, r .The APF function is U : R d → R, where the attractive and repulsive potential fields are U att and U rep , respectively.The potential field forces generated by these fields are F att and F rep .
The following paragraphs describe some problems in path planning in a dynamic environment.
In contrast to the static environment, obstacles in a dynamic environment change their positions over time.Therefore, for the dynamic path-planning algorithm, time measurement must be introduced into the nodes so that the collision detector can use the time component to check the feasibility of the planned path.The dynamic environment proposed in this paper is described in Problem 1.
Problem 1 (Dynamic environment).If threats change location over time, the environment is called dynamic, i.e. ( The feasibility of the path is described in Problem 2. Problem 2 (Feasible path planning).In the configuration space X containing X obs (t) and X f ree (t), given x init and X goal x goal , t , a collision-free path σ : [0, T ] → X f ree can be planned, where σ (0) = x init , σ (T ) ∈ X goal x goal , T .
Let be the set of all feasible paths in obstacle-free space.The cost function c(σ ) ∈ R >0 represents the generation value of each feasible path σ .Then, the optimal path-planning problem is described by Problem 3.

Problem 3 (Optimal path planning).
Assuming that there is a solution set for Problem 2, the path in which the cost function c(σ ) is the smallest is the solution σ * of the optimal path planning, i.e. ( In UAV path planning, constraints are placed on the maximum steering angle φ max and the maximum climbing/diving angle γ max [30].Assume that the coordinates of any path node are (x i , y i , z i ).The constraints on φ max and γ max are calculated as follows: where the vector a i is Therefore, in this paper, the cost function between two points x 1 and x 2 in the configuration space is defined as follows: where φ and γ are the steering angle and climbing/diving angle of the UAV flying along the planned path, respectively, which can be calculated through the geometric relationship between the two adjacent path segments.If these angles exceed their limits, it is considered that the planned path cannot meet the needs of the UAV, so the cost function is set to infinity.Thus, the cost function of the planned path can be expressed as:

HPO-RRT* algorithm
In this section, we introduce the HPO-RRT* algorithm proposed in this paper in detail.

Overall framework
To address the path-planning problem of UAVs in dynamic environments with moving threats, the algorithm must be able to solve dynamic problems quickly.However, only obtaining the feasible path is not an optimal choice: the algorithm should also be able to quickly obtain the optimal solution.Note that many time-based RRT* algorithms need to prune the tree structure when optimizing the path, which inevitably results in more wasted branches and additional calculations of time-based nodes in the correction tree.Therefore, the HPO-RRT* algorithm is proposed in this paper to satisfy the efficiency of real-time path planning and the superiority of paths in the dynamic environment of UAVs.The algorithm adopts a hierarchical framework for path planning and optimization.First, the UAV perception system is used for threat prediction and world information updates.Then, an improved time-based RRT* algorithm is used as the path planner.Finally, a low-cost path optimizer is used to optimize the path obtained by the planner.The pseudocode of the algorithm is shown in Algorithm 1, including one main cycle and five main processes.
Similarly, the root node of tree T is set to x init , and the current positions x cur of the UAV and time are initialized.
(2) Termination condition (line 4).After detecting that x cur has reached the target area X goal , the program is immediately terminated.Otherwise, continue running the program until x cur arrives at X goal .(3) UAV perception system (lines 8-15).During any execution, the UAV will fly a certain distance according to the planned path, and then the planner will update the final state x cur of the drone in this execution to the root of the new tree T * new .Next, nodes with timestamps less than x cur in the tree and their descendants need to be deleted.Then, taking zero as the initial depth of x cur , the trajectory of the threat in the future time N l t is predicted through the navigation system, where N l is the maximum depth.Because of the movement of threats and UAVs, all information of maps, trees, UAVs and moving threats must be updated .This process is basi-cally the same as that of Risk-RRT [20] except that the perception system is expanded from 2 to 3D.This paper mainly focuses on improving path planning and optimization, so the description of the perception process refers to [20].(4) Path planner (lines [16][17][18][19][20].The existing time-based RRT algorithm usually uses RRT as the bottom planner [15,31].Obviously, the planned path is not optimal, and a feasible path may not be obtainable within the time interval.In this paper, we use an improved timebased RRT* as the path planner.First, the bias sampling method based on the APF function is used to guide sampling node x ap f rand towards the target region and away from the obstacle.Then, the time-based RRT* process of selecting the optimal parent node in the neighbourhood of the new node is introduced to optimize the tree structure, and the rewiring process in time-based RRT* is removed to avoid damage to the tree structure.Meanwhile, the tree structure is extended in combination with the lazy collision checking process.Finally, the heuristic path is obtained through the tree structure.(5) Path optimizer (line 21).Because all the nodes in the tree have time parameters and the time is irreversible, using the rewiring process in the tree structure to optimize the path will lead to the disorder of node timestamps in the tree, and repairing the tree structure will consume considerable extra time.Therefore, this paper designs a low-cost path optimizer that directly optimizes the heuristic path without changing the tree structure to obtain the optimal path σ opt .

Bias sampling based on the APF function
In the traditional RRT algorithm, the sampling nodes are randomly selected in the collision-free configuration space, and then the extended nodes generate a search tree to the target area.Although this method can ensure the completeness of probability, in the three-dimensional space of UAV flight, the search efficiency is low because of the large random sampling space.In addition, HPO-RRT* is a planning algorithm based on time constraints.The path-planning time of the UAV will be limited to a fixed time interval, so the time of path planning is relatively short.Random sampling is inefficient and may be unable to plan a feasible path.Therefore, to plan the path efficiently, this paper proposes a bias sampling method based on APF functions, which makes the sampling nodes expand towards the target region and away from the threat to plan the path efficiently in a limited time.
The APF constructs the attractive potential field of the target position and the repulsive potential field of the threat.The resultant force F ap f generated by these two potential fields guides the UAV to move towards the target point.In this paper, the following functions are used to define the attractive U att and the repulsive U rep potential fields, as well as their corresponding force vectors F att and F rep , respectively.U att and F att are: U rep and F rep are where ρ o is the maximum impact distance of the threat, and x obs is the location of the threat.According to Eqs. ( 7) and ( 9), the calculation of random sampling through APF bias is described as follows: Then, the bias sampling node can be described as: where k bias ∈ R >0 is the bias step.Note that path planning using the APF function has disadvantages, such as inaccessible targets and easily falling into local minima.However, in this paper, the APF function is only used to bias the random sampling, not to find the path.At the same time, each sampling is an independent random event.Therefore, the shortcomings of the APF function can be effectively avoided.
In this paper, we further improve the approach described in [28].Qureshi and Ayaz [28] calculate the influence of the attractive potential field on x rand and determine whether to reserve bias sampling nodes by determining whether the distance between x rand and x obs is within the range of ρ o .Through the strategy of [28], the random sampling nodes outside the influence range of the repulsive potential field can be biased towards the goal, thereby improving the sampling efficiency.However, if the nodes are in the influence range of the repulsive potential field, they will be discarded.At this time, new nodes will be randomly sampled again to bias.Note that if the random sampling node within the influence range of the repulsive potential field can also be biased, repeated random sampling and collision check of this node will be avoided, thus further improving the sampling efficiency.Therefore, this paper introduces biased sampling of the repulsive potential field to retain all random sampling and expand the expansion space of the tree as much as possible.The pseudocode of the APF bias sampling is shown in Algorithm 2.
First, the sample_random function is used for random sampling in collision-free space X f ree .Second, the attractive force F att of x goal to x rand and the repulsive force F rep of x obs to x rand are calculated.Third, the resultant force F ap f is calculated and then used to bias x rand to x ap f rand .The final bias sample node x ap f rand is obtained.

The expansion of the time-based tree
Similar to RRT*, time-based Risk-RRT* can also optimize the tree structure through the rewiring process to obtain asymptotically optimal partial planned paths when addressing path-planning problems with moving threats.However, because of the irreversibility of time, the rewiring process of Risk-RRT* will destroy the time-based tree structure.The time spent repairing the tree structure increases exponentially with the size of the tree, reducing the planning efficiency.In addition, in the sampling-based path-planning algorithm, collision checking is very time-consuming.Because the planning time using HPO-RRT* is a fixed short time interval, the planning efficiency is particularly important.Therefore, in this paper, the neighbouring parent selection process of Risk-RRT* is retained in the process of expanding the tree to optimize the tree structure, while the rewiring process is eliminated to ensure the efficiency of the path planner.Although this approach will sacrifice the optimality of the path, the path can be quickly optimized by a low-cost path optimizer that is more efficient than the rewiring process.At the same time, to expedite collision checking, a method of lazy collision checking is proposed in this paper, which improves the computational efficiency by ensuring that the planned path is collision free.
In the HPO-RRT* algorithm, the expansion of nodes is closely related to time.Each forward expansion of a node in the tree indicates that the execution time of the UAV increases by the time step t.Therefore, any node x i ∈ T .V in the tree needs to record its timestamp, i.e. x i : x * i , t + N i t , where x * i is the position coordinate of x i , t is the start time of this extension process, and N i is the number of time steps required to extend from the root node x root of the current spanning tree to x i .
Note that the Extend process is improved by using timebased RRT* as the basic planner.It only includes the process of selecting the optimal parent node in the neighbourhood of the new node x new and omits the rewiring process.This approach is used because only after selecting the optimal parent node of x new , x new will be added to the time-based tree structure to obtain the timestamp.However, the time-based tree needs to update the timestamp of the node and all its descendants when rewiring the node and conducting collision checking again, which greatly increases the calculation time.A simple but typical example is shown in Fig. 2. In Fig. 2a, after adding x new to the time-based tree, the rewiring process will be executed in the neighbourhood of x new .Obviously, the cost from x 0 through x new to x 8 is less than the cost from x 0 to x 8 in the original tree.If rewiring is performed, i.e. the new parent node of x 8 is x new (called x 2 in (b)), the time-based 123 Fig. 2 Rewiring process of the time-based tree tree will be updated to the structure shown in Fig. 2b.At this time, the timestamps of x 8 , x 9 , x 10 and x 11 are affected by rewiring and need to be updated to the timestamps of x 3 , x 4 , x 5 and x 6 , respectively.In addition, because of the change in the timestamp, the collision checking of these update nodes must be reperformed to ensure that the nodes will not collide with the moving threats under the new timestamp.If collision checking cannot be passed, all branches formed by this node and all its child nodes need to be deleted.This process brings many additional computing costs.Therefore, in the HPO-RRT* algorithm, we delete the rewiring process of the time-based RRT* to reduce the computational consumption and meet the requirements of efficient and fast tree structure expansion within a fixed time interval.In addition, we keep the process of x new reselecting the new parent node.This process is retained because it is a forwards expansion process, which will only affect the relevant parameters of the new nodex new .Therefore, it will not bring additional computing consumption.Moreover, it can optimize the tree structure to a certain extent.
The planning space in this paper contains static and moving threats.For static threats, the position remains fixed at each time step.However, the coordinates of moving obstacles change with time.During the planning of the current time interval, the UAV can only refer to the threat prediction trajectory before the current planning to expand the tree by N l time steps t (Algorithm 1, line 12).However, as the predicted threat trajectory may not be accurate enough, the path obtained in this partial planning may become infeasible because of the new prediction of the threat trajectory in the next planning.An intuitive idea is to set N l very small to improve the accuracy of obstacle prediction and reduce the number of potentially invalid and repeated collision checking calls to improve the calculation efficiency.However, this approach makes the expansion range of the tree in each cycle very small, which is not suitable for the Fig. 3 Illustration of lazy collision checking large-scale 3D planning space of UAVs.Therefore, to ensure a certain extension range of the forwards tree and reduce the computational loss, this paper proposes a method of lazy collision checking, which delays the collision checking process until it must be called.Figure 3 illustrates lazy collision checking.This method defines a new extended time depth N c < N l .In N c t, all extended branches perform collision checking.The step length of each UAV execution is also set in N c t.However, collision checking will not be performed in the time depth from N c t to N l t.Because of the parallel planning and execution, the planning of the next time interval will start with a new root node, and the time depth for collision checking will continue to expand backwards.Therefore, all nodes on the final expansion tree will complete collision detection.Note that N c greatly influences the performance of HPO-RRT*, so the calculation complexity and planning depth should be fully weighed.
The pseudocode of the Extend process is shown in Algorithm 3.After the SampleAPF process returns the biased sample x ap f rand , the forward expansion process begins.Similar to the time-based RRT*, the process of expanding the new node x * new is realized by the Nearest, Steer and Near functions.In the expansion process of x new , only its position is considered temporarily (i.e.let x new be x * new ).First, the Nearest function returns the closest node x * nearest between x ap f rand and the tree node set T .V , and then Steer generates a new node on x * new .Next, the algorithm searches the neighbourhood formed by the sphere B x * new , r whose centre is x * new and radius is r (where r = γ log n n 1/d [32]) to obtain the node set X near whose distance from x * new in B x * new , r is less than r .After connecting x * new to X near , the timestamp t new of x new is arranged from small to large to form a potential parent node set X parent .First, the algorithm judges whether t new is within the time depth N c t.If this criterion is met, the algorithm performs collision detection to select the optimal parent node of x new ; if t new ∈ (N c t, N l t], it directly selects without collision detection.If t new > N l t, x new exceeds the maximum time depth limit of this planning.Finally, x new and its edge formed with the optimal parent node x parent are added to the tree T .

Planning path finding
Because the growth of the tree can only occur in the fixed planning time of a single cycle, and the maximum time depth of the node is limited to N l t, for the 3D large-scale space X goal in which the UAV flies, it is almost impossible to generate an expansion tree in which the node falls in the target area SD after only one or several cycles.Therefore, this paper introduces a cost function composed of the true distance cost and the estimated distance cost to find the heuristic partial path in the single planning spanning tree, namely, the HeuristicPathFinding process.
In a planning interval, tree expansion is limited by the maximum expansion depth N l .That is, in the nodes of tree T new , the maximum extension time depth is N l t.Based on this limitation, the node set of the true path cost can be defined as X l , including n l nodes x l with a time depth of N l t, i.e.X l := (x * l, 1 , N l t), ...(x * l, i , N l t), ..., x * l, n l , N l t .The set of true costs of all paths that can reach x l, i in the current tree is defined as g(X l ).The set of estimated costs from any node in X l to the target region X goal is defined as h(X l ).Then, the cost f (X l ) of the path is: We aim to find the node x l, i that can minimize f (X l ), i.e.
x l, min = arg min For the time-based tree, the exact value of g(X l ) can be obtained by calculating the path length, while h(X l ) can only be estimated.Therefore, an acceptable estimation cost calculation method should be selected so that the real cost of reaching x goal is not overestimated.The following formula is used for calculating the estimated cost h x l, i of any node x l, i : where k l is the heuristic parameter.
On the basis of the above calculation, we ensure that in each partial planning, we can find a heuristic partial path towards the target from the current time-based tree.Algorithm 4 gives the pseudocode of the HeuristicPathFinding process.First, the algorithm finds all nodes with time depth t l in tree T new and forms a set X l .Then, the node x i, min that minimizes the estimation cost is found by f (X l ).Finally, the heuristic path σ heu generated by the planner is found by x i, min .

Low-cost path optimizer
The path planner using the HPO-RRT* algorithm efficiently generates heuristic paths without guaranteeing the optimality of the paths.Although the time-based Risk-RRT* algorithm can obtain the optimal path, as an online path-planning algorithm used to address moving threats, it takes a long time in the optimization process and sometimes cannot optimize the path well in a short planning time.Qi et al. [33] uses RRT* to find the initial path and then uses the ant colony (ACO) algorithm to reoptimize the path.This hierarchical optimization is interesting, but the ACO algorithm brings high computational complexity.The Bellman-Ford (BF) algorithm [34], as an algorithm for solving the shortest path, provides a method for searching the shortest path in a weighted graph.Therefore, inspired by [33,34], combined with the requirements of online dynamic path planning, this paper proposes a low-cost fast path optimization method.This method directly optimizes the heuristic path σ heu generated by the path planner without changing the tree structure.Therefore, additional complex calculations are not needed to optimize the path more quickly and effectively.The BF algorithm estimates the distance from any node to a specified node by relaxing the trigonometric inequality constraint of the connection between nodes in the weighted graph.In the path optimization process of HPO-RRT*, we use this algorithm to offset the path points to optimize the heuristic path.As shown in Fig. 4, the blue path is the heuristic path σ heu that has not been optimized, and the green translucent path is the path σ * that has been optimized once.In the cycle, bias optimization is performed on the nodes in σ heu (excluding the initial point and the target point) in turn.After all the path nodes are optimized once, path B is obtained.Obviously, cost(σ * ) < cost(σ heu ).
To clearly illustrate the optimization process, the offset to a is taken as an example.For the three nodes x 1 , x 2 and x 3 in the tree, the path cost from x 1 directly to x 3 is minimal.However, if x 1 and x 3 are directly connected (that is, connected by the orange dotted line), the timestamps of x 3 and all its descendants need to be updated, resulting in many extra computations.In addition, if X obs exists nearby, such as the handling of node x 4 , the directly connected edge (x 3 , x 5 ) will collide with the threat, causing the optimization path to fail.Therefore, we try to optimize the path by biasing the position of node x 2 towards edge (x 1 , x 3 ).The specific method is to find the midnode x mid 2 of edge (x 1 , x 3 ).Then, move x 2 along the vector − −−− → x 2 x mid 2 by a step λ opt to obtain the optimized node.For any node x i in the path, the optimization step λ opt is calculated as follows: where k opt ∈ R >0 a is the bias factor.Then, the optimized node x i of any node x new in the path can be expressed as: Algorithm 5 gives the of the Low-costOptimizer process.Define a function to determine whether the optimization process is completed and assign it as false at the beginning of the algorithm.The function CollisionCheck(x i−1 , x new , x i+1 ) means collision checking for the path segment from x i−1 through x new to x i+1 .If the result is true, then (x i−1 , x new , x i+1 ) is a collision-free path segment.The function MaxAngleCheck(x i−1 , x new , x i+1 ) indicates the UAV maximum angle constraint detection for the path segment (x i−1 , x new , x i+1 ).Note that line 18 compares the cost between the path before the current cycle and the optimized path to complete the current cycle.If the cost difference between the two paths before and after optimization is less than the threshold ε, where ε ∈ R, then the path is hardly changed in the optimization process.When this result is obtained, the improvement of path quality in the optimization process is no longer obvious, i.e. the optimal path is essentially determined.Therefore, assign FinishOptimization to true.
In any optimization cycle, first, the initial position is added to the optimization path σ opt .Next, the nodes in heuristic path σ heu are biased optimized.If the optimized path meets the constraints of collision detection and UAV maximum angle, then the optimized path of this cycle planning is feasible.Finally, if the generated path can improve the path quality, start the next cycle with this path as the path to be optimized.The termination conditions of the optimization process shown in Algorithm 5 include the following criteria: (1) the process is terminated when collision detection or the maximum angle constraint of the UAV cannot be met; (2) in the same cycle, if the costs of the optimized path σ opt and the nonoptimized path σ * differ by less than ε, the process terminates.
Obviously, throughout the optimization process, the original tree structure is unchanged, and the timestamps of all nodes are unmodified.We can obtain a better path σ opt than heuristic path σ heu by performing bias optimization on the nodes of the heuristic path.Despite additional computations, the computational cost of our method is negligible relative to methods that optimize tree structures to optimize paths.

Analysis
In this section, we analyse the probability completeness and homotopy optimality of HPO-RRT*.

Probabilistic completeness
Probabilistic completeness is mainly used to analyse the ability of algorithms to find feasible solutions.The HPO-RRT* algorithm finds the path in a single plan through the path planner and then optimizes the path directly generated by the path planner.Therefore, the probabilistic completeness analysis of HPO-RRT* only needs to consider some paths obtained through the path planner in a single plan.
Let V ALG n represent the set of tree nodes generated by algorithm ALG after iteration n.Definition 1 gives the definition of probabilistic completeness.Definition 1 (Probabilistic completeness).Given the initial node x init and the goal region X goal , if algorithm ALG can find the feasible path from x init to x goal ∈ X goal for any path-planning problem with a feasible solution, i.e.
then the algorithm is considered probability complete.
The RRT algorithm has been proven to be probability complete.HPO-RRT* covers all the key processes of RRT in the bottom planner.Similar to RRT, HPO-RRT* has probabilistic completeness, which is described in Theorem 1.

Theorem 1 (Probabilistic completeness of HPO-RRT*).
Given a path-planning problem, if there is a feasible solution, then, Proof of Theorem 1.The proof of Theorem 1 is based on the following three arguments.
(1) As shown in Algorithm 1, the vertex set V PT L−R RT * of the random tree T generated by HPO-RRT* includes node x init , and V PT L−R RT * 0 = x init , which is the same as RRT.
(2) Similar to RRT, the tree generated by HPO-RRT* is connected.In other words, any random sampling node can be connected to the tree.(3) HPO-RRT* only plans the path within a fixed time interval of each cycle, and the time depth limits its expansion range.However, for any cycle i, HPO-RRT* is the same as the RRT algorithm and has probability completeness.That is, in cycle i, when sample n i tends to infinity, the probability of HPO-RRT* finding a path from the current planning starting point (also the subtarget point of the last cycle planning, i.e. x subgoal, i−1 ) to the subgoal node x subgoal, i is one, i.e.
When the number of cycles is infinite, the subgoal node x subgoal, i can surely fall into the goal region X goal , and all partial paths are connected.Therefore, the probability that HPO-RRT* finds a path from x init to X goal can be expressed as: Thus, Theorem 1 is proved.

Homotopy optimality
The asymptotic optimality of RRT* and many of its derivative algorithms has been proven.This section proves that HPO-RRT* has similar properties to these algorithms, which is called homotopy optimality.
Note that HPO-RRT* performs partial path planning within a fixed time interval of each planning.Therefore, the path that ensures the homotopy optimality of HPO-RRT* is obtained by optimizing the heuristic path in a fixed time interval planning, rather than the entire path from the initial position to the goal region.The relevant definitions of homotopy optimality and the proof process are as follows.
Definition 2 (Homotopy class of feasible paths [35]).For two arbitrary paths σ 1 and σ 2 with a fixed initial node and goal node, if one path can be continuously deformed into the other without intersecting any threat, then σ 1 and σ 2 are considered to belong to the same homotopy class [σ ].
When multiple paths belong to the same homotopy class [σ ], there must be a homotopy optimal path σ * opt , which is the least costly path in [σ ].Thus, Lemma 1 is given as follows.

Lemma 1 (Homotopy optimality of σ opt ). Given that the homotopy class containing heuristic path σ heu is
[σ ] PTL−RRT * σ heu ∈ [σ ] PTL−RRT * , the optimization path σ opt returned by the path optimizer of HPO-RRT* is the homotopy optimal path of the homotopy class, which is the path with the lowest cost.
Proof of Lemma 1.In process re-optimization, we use a low-cost path optimization algorithm to quickly optimize heuristic path σ heu .When the path meets the termination conditions of the process, an optimized path σ opt is generated.Therefore, the proof of Lemma 1 can be divided into three arguments according to the conditions of path termination optimization.
(1) After the optimization process, the path cost is almost unchanged, i.e. cost < ε.As we analysed in "Lowcost path optimizer", if cost < ε, then the path has reached optimality, so the generated path σ opt is the homotopy optimal path in the homotopy class.(2) When the return value of CollisionCheck is false, path σ opt is obtained.First, we assume that the termination condition is triggered because edge cannot pass CollisionCheck.As shown in Fig. 5a, there are two kinds of optimized path (i.e.σ opt is the blue path and σ opt is the green path).Suppose the σ opt whose cost is cost σ opt < cost σ opt and cost σ opt − cost σ opt > ε, that is, σ opt must be closer to edge (x 1 , x 4 ) than σ opt .However, σ opt cannot pass the collision checking, and σ opt will collide with obstacles.This causes the CollisionCheck procedure returns a false value.Then, process re-optimization stops and returns the optimized path σ opt before the start of the current cycle as the optimal path.Therefore, no path costs less than σ opt .Then, σ opt can be claimed to be homotopy optimal in the homotopy class.(3) When the return value of Max AngleCheck is false, path σ opt is obtained.First, we assume that the termination condition is triggered because edge cannot pass Max AngleCheck.Similar to Proof 2), Fig. 5b shows two kinds of optimized path (i.e.σ opt is the blue path and σ opt is the green path).The edge (x 0 , x 1 ) is the last segment of the optimized path obtained from the previous cycle planning It is assumed that there is a path σ opt , and cost σ opt < cost σ opt .Then, σ opt must be closer to edge (x 1 , x 4 ) than σ opt .However, the steering angle of σ opt is φ , which is larger than the steering angle φ max (i.e. the maximum steering of UAV) of σ opt .If the angle φ does not meet the UAV angle constraint, Max AngleCheck returns false.Obviously, φ does not meet the constraint.Then,σ opt is infeasible, and re-optimization returns to path σ opt .Therefore, no path costs less than σ opt .It can be proven that σ opt is homotopy optimal.Based on the above three arguments, σ opt is proven to be the homotopy optimal path in the homotopy class containing σ heu .

The results of simulation experiments
In this section, we give the simulation experiments of the proposed algorithm.The main purpose is to verify the feasibility and advantages of the algorithm in a UAV flight environment.Therefore, "Simulation setup" shows various threats in the considered UAV flight dynamic environment and the settings of some basic parameters to make the results as close to the real battlefield environment as possible to increase the reliability.

Dynamic environment design
Constructing a battlefield environment is a prerequisite for path planning.The closer the construction of the battlefield environment is to reality, the better the verified path-planning algorithm will be implemented on UAVs.Therefore, to verify the feasibility and advantages of the HPO-RRT* algorithm in a battlefield environment, this paper simulates the real dynamic environment, including real terrain, radar, air defence missiles, anti-aircraft guns and a tower, as the configuration space for UAV path planning.We have already introduced the modelling and calculation of these threats in our previous work [36].Since this paper mainly focuses on the planning of the HPO-RRT* algorithm, the modelling of threats will not be further described.
Based on the five threats in the above-mentioned real battlefield, this paper designs two 3D dynamic environments, including a threat with known moving trajectories and a random moving threat, to test the performance of the proposed algorithm and compare HPO-RRT* with related algorithms.The terrain threat in the two scenarios uses a unified elevation map.The two scenarios contain static or moving radar, missiles and anti-aircraft artillery threats.
Experiment Scenario 1 Dynamic flight scenario with known moving trajectory threats.
The environment size is 400*400*60 m 3 .Radar 1, Missile 1 and Anti-aircraft gun 1 are moving threats, which start to move from their initial positions.The speed and moving time period are shown in Table 1.The remaining threats are static.When the moving threats touch the boundary of the configuration environment, they will rebound and then move in the opposite direction.There is no interaction between different moving threats.

Experiment Scenario 2 Dynamic flight scenario containing multiple random moving threats
The size of the environment is 400*400*60 m 3 .The initial position of the terrain threat and the other three threats is the same as that of experimental Scenario 1.In Scenario 2, moving threats are randomly selected from all radars, air defence missiles and anti-aircraft guns in the environment.{1, 3, 5, 6} threats can be selected to move at the same time.The movement speed is randomly selected from {1, 3, 5, 7} m/s 2 .The start time of the movement is any time during the planning period, and the time period of each movement is 3 s.All threats except those moving remain stationary.When the moving threats touch the boundary of the configuration environment, they rebound in the opposite direction and then continue to move.There is no interaction between different moving threats.

Basic parameter design
During the preparation of simulation experiments, the performance of the UAV, the location of threats in the environment, and the relevant basic parameters of the HPO-RRT* algorithm must be set.Table 2 shows the specific parameter names and values.

Simulation experiments
In this section, the HPO-RRT* algorithm is applied to the two dynamic scenarios with moving threats designed in "Simulation setup" to verify its feasibility, optimization, efficiency and success rate of the HPO-RRT* algorithm in the dynamic environment of UAV flight.Additionally, to further verify the performance of the HPO-RRT* algorithm, we compared it with three dynamic planning or replanning algorithms derived from RRT, including RRT X , Risk-RRT and Risk-RRT*.In the two sets of scenarios, these algorithms are run many times to reduce the randomness of sampling-based algorithms.We evaluate the performance of each algorithm by taking the navigation time, planning path length and planning success rate as metrics.
Experiment Scenario 1 Figure 6 shows the UAV paths planned by the four algorithms in Scenario 1, and Table 3 shows some performance data in the experiment in Fig. 6.In general, compared with the other three algorithms, the path planned by HPO-RRT* is shorter and smoother.Moreover, in each partial planning, HPO-RRT* can quickly plan a longer path, so the planning efficiency and expansion range of this algorithm are also the best.In addition, HPO-RRT* has the   RRT X is a replanning algorithm.Compared with similar algorithms, RRT X has faster information transmission speed and response speed and has better performance in the problem of path replanning in addressing moving threats.Although the tree structure needs to be pruned in RRT X path planning, which greatly reduces the corresponding speed and planning efficiency of the algorithm, Fig. 6 shows that RRT X can still maintain a good response speed and planning efficiency.In addition, when the tree structure is small at the initial stage of planning, the planned path quality is basically the same as that of the HPO-RRT* algorithm.However, when the sampling increases and the tree structure becomes larger, RRT X takes a long time to prune the tree structure, thus reducing the time to optimize the path.Thus, the quality of the paths planned by RRT X begins to lag behind that of HPO-RRT*.The path planned by RRT X becomes less smooth, and the angle loss increases when the UAV flies along this path.Risk-RRT, Risk-RRT* and HPO-RRT* are time-based RRT algorithms.In other words, little time needs to be spent in tree structure pruning during planning, so they have a fast response speed.However, Risk-RRT is an algorithm whose path planner is time-based RRT, thus using Risk-RRT for planning will generally only obtain one feasible path.Furthermore, the optimality and planning efficiency of the planned path cannot be guaranteed.Therefore, we can see that the paths and navigation times obtained by Risk-RRT and HPO-RRT* are considerably different.Compared with Risk-RRT, the quality of the planned path of Risk-RRT* has been improved to some extent.Risk-RRT* is an algorithm that needs to rewire the tree structure to optimize the path, but Fig. 6 shows that the improvement of the path quality is not large.This result is obtained because part of the planning time is short, and the time allocated to the reconnection process is limited.However, Risk-RRT* requires a large number of rewiring calculations to optimize the path, so it cannot optimize the path better in a limited time.Similar to RRT X , Risk-RRT* has good algorithm performance in the early stage due to its small tree structure.However, with increasing sampling, its performance gradually lags behind that of HPO-RRT*.It can be seen that HPO-RRT* always performs well in planning.HPO-RRT* uses APF bias sampling in the path planner to improve the quality and efficiency of sampling.At the same time, in the expansion of the tree, the lazy collision checking method is used to reduce the computational cost and expand the expansion range of the tree.After obtaining the heuristic path, HPO-RRT* uses the lowcost path optimizer to directly optimize the heuristic path to quickly obtain the optimized path.This optimization method does not change the structure of the time tree; that is, no additional calculation loss occurs in the optimization process.Therefore, compared with the other three algorithms, HPO-RRT* has the smallest navigation time, the shortest planned path and the highest smoothness.
Furthermore, in Scenario 1, each algorithm is repeated 50 times.We use the planned path length, navigation time and success rate to evaluate the performance of each algorithm.Figures 7, 8 and 9 compare the results.We find that the HPO-RRT* algorithm has good stability and algorithm performance.In the length comparison of 50 experiments, HPO-RRT* can use the low-cost path optimizer to fully and quickly optimize the path during planning so that the length of the planned path is always the shortest and the length variance is very small.Additionally, because HPO-RRT* maintains high planning efficiency, the navigation time can be quickly stabilized.Moreover, for the above reasons, the success rate of the algorithm is high.In other words, HPO-RRT* has good performance and stability when addressing UAV path planning in a dynamic environment with known moving threat trajectories.
Experiment Scenario 2 Experiment Scenario 2 focuses on the success rate and optimization of path-planning algorithms in a more complex dynamic environment.The algorithm  performance is tested by setting the number and speed of moving threats in the UAV flight environment to a random generation state.Figure 10 and Table 4 show the performance of the four algorithms in one case with threats of random movement.Furthermore, for each combination of random threat number and speed in Scenario 2, each algorithm is tested in 6 trials.Therefore, for each algorithm, we conducted a total of 96 experiments on 16 combinations of random threats.We use the planned path length, navigation time and success rate to evaluate the performance of each algorithm.The performance of the algorithms is shown in Figs.11 and 12.In general, compared with Scenario 1, the performance of the four algorithms declines in Scenario 2. This decline occurs because the four algorithms need to address more complex moving threats.However, compared with the other three algorithms, HPO-RRT* has the shortest and smoothest planned path length.The navigation time is also the shortest, and the success rate is also guaranteed.Therefore, the efficiency, effectiveness and stability of HPO-RRT* are demonstrated.Figure 10 and Table 4 clearly show that the planning path quality of the HPO-RRT* algorithm remains superlative among the four algorithms.Because Risk-RRT has relatively low planning efficiency and no path optimization process, its planning path and navigation time are relatively poor.In more complex dynamic environments, Risk-RRT* does not have sufficient reconnection time, resulting in less obvious improvement of path quality.Under its ability of rapid replanning and optimization, RRT X can obtain a faster navigation time and higher-quality paths than Risk-RRT and Risk-RRT*.With its efficient planner and optimizer, HPO-RRT* also maintains a high-quality planning path and has the shortest response time in a dynamic environment with random moving threats.
Furthermore, to verify the performance quality of HPO-RRT* more comprehensively, 96 simulations were conducted on 16 combined random moving threats.Figures 11 and 12 show the results.In the combination of a few moving threats and slow moving speed, the four algorithms easily plan the path, and all find the path with basically the same success rate, with the path length of HPO-RRT* being the shortest.However, with an increase in the number and speed of threats, the success rate of the algorithms declines.In addition, it can be seen that HPO-RRT* guarantees almost the same success rate as RRT X , while the success rates of Risk-RRT and Risk-RRT * decrease significantly.Although HPO-RRT* as a time-based RRT algorithm has been proven to be inferior to the replanning algorithm (RRT X ) in addressing path-planning problems in highly dynamic environments (i.e.environments with fast-moving threats and more moving threats), HPO-RRT* can maintain good adaptability to dynamic environments with its high-quality sampling and efficient expansion, and the performance loss of this algorithm is low.In terms of path quality, HPO-RRT* maintains certain advantages in various combinations of random moving threats, which depend on the low-cost path optimizer of HPO-RRT* to quickly and fully complete path optimization at any time.Meanwhile, HPO-RRT* has the shortest and most stable navigation time with the good performance of the efficient path planner and the low-cost path optimizer in its hierarchical planning architecture.In other words, HPO-RRT* can adapt to more complex dynamic environments and has good performance and stability.

Conclusions and future work
In this paper, a time-based RRT named HPO-RRT* is proposed to solve the real-time path-planning problem of UAVs   in a dynamic environment.This algorithm adopts a hierarchical framework including the UAV perception system, path planner and path optimizer, which can optimize the path while improving efficiency.In the path planner, an improved time-based RRT* algorithm is used to obtain heuristic paths.The algorithm of the path planner improves the sampling process, tree expansion process and collision checking process to avoid additional calculation loss and improve planning efficiency.In the low-cost path optimizer, to avoid numerous calculations caused by the optimization tree structure, the heuristic path is directly optimized to quickly and efficiently obtain a homotopy optimal path.Simulations and comparisons show that the proposed algorithm has better performance.
In future work, we will conduct outdoor experimental verification on the independently developed quad-rotor platform to confirm the application value of the algorithm.In addition, we expect to extend HPO-RRT* to the formation and cooperation of UAVs in future work.An intuitive prospect is to use a formation framework to complete the cooperative path 123 planning of UAVs.Of course, a better prospect is to bring HPO-RRT* into the full autonomous cooperative architecture.

Fig. 6
Fig. 6 Experiment Scenario 1. Planned path of different algorithms at different times, in which the blue line represents RRT X ; the green line represents Risk-RRT; the yellow line represents Risk-RRT *; and HPO-RRT* is indicated by a red line.The grey translucent threat

Fig. 7
Fig. 7 Planned path length of 50 experiments of four algorithms

Fig. 10
Fig. 10 Experiment Scenario 2. path diagram of different algorithms at different times, in which the blue line represents RRT X ; the green line represents Risk-RRT; the yellow line represents Risk-RRT*; and HPO-RRT* is indicated by a red line.The grey translucent threat

Fig. 11
Fig. 11 Relationship between the success rate, average path length and obstacle speed and number

Fig. 12
Fig. 12 Comparison of the mean and variance of the navigation time

Table 1
Parameters

Table
Length and navigation time of planned paths with the algorithms in Fig.6477.9represents the shortest path length.21.2 indicates the shortest running time shortest planning and execution time, stable performance, and a high success rate.

Table 4
Length and navigation time of planned paths with the algorithms in Fig.9