8 Path_planning_algorithm_under_specific_constraints_in_weighted_directed_graph
8 Path_planning_algorithm_under_specific_constraints_in_weighted_directed_graph
Authorized licensed use limited to: b:on Instituto Politecnico de Coimbra. Downloaded on April 04,2023 at 12:46:56 UTC from IEEE Xplore. Restrictions apply.
In a directed weighted graph, we need to consider the required to pass through the specified intermediate node set
weighted length of the path which starts from a specified {2,3} without repetition, then there are two paths:0→2→3
node, after passing through some specified intermediate →1 and 0→3→2→1,whose LinkID are 1|5|4 and 2|6|3
nodes without repetition (the number of intermediate nodes respectively,with corresponding weighted value 4 and 5.So
is uncertain), and finally ends at the specified termination the shortest path P is obtained as 0→2→3→1, it is an
node.And we aim to find such a feasible path hoping that acyclic directional path.
its weighted length would be as small as possible and
search time would be as short as possible. 2.1. Concept Definition
Most of the current studies of the shortest path
algorithm are just aimed at reaching the end point from the The following are definitions of some related terms:
starting point. Although some researchers considered path LinkID: Index of directed edge;
planning problem with constraints, such as the constraint- Cost: Weight of directed edge;
based multi-robot path planning algorithm given by the SourceID: The index of the starting node of the
M.Ryan[14]; Huang Shuli, Hu Dasha et al.[15] have directed edge;
studied about the shortest path algorithm in the unweighted DestinationID: The index of the terminate node of the
undirected graph with the constraint that the path must directed edgs;
cover a specified intermediate node set. But very few IncludingSet: Node set which must be passed
researchers have considered the complex path planning through, and different node indexes are split with the
problem with various constraints in a weighted directed Separator '|'.
graph as mentioned in the last paragraph. Weight of path: The sum of the weights of all the
For path planning issues raised above, due to the directed edges of the path.
additional constraints, it can not be solved by directly using Acyclic directed path P: A directed path from s
Dijkstra algorithm. Therefore, this paper takes the Dijkstra (starting node) to t (termination node) which is composed
algorithm as the foundation, and draws on the heuristic of a series of directed edges, and is not allowed to
search theory of A* algorithm, and finally gives a feasible repeatedly pass through any of the nodes; The
algorithm consulting the area routing calculation process representation of P is to list the index of the directed edges
of the OSPF routing protocol. In this algorithm, we put constituting the path in turn, and the indexes are split with
forward a new concept of "inspired hop count" ,it is the the Separator '|'.
most important parameter in this algorithm to dynamically
adjust the balance between the weight of path and the time 2.2. Problem Definition
of path searching. The experimental results show that this
algorithm can effectively solve the above path planning Given a weighted directed graph G = (V, E), V for the node
problem.
set, E for the directed edge set, each edge has its own
weight. For a given node s, t, a node set V and its subset ,
2. PROBLEM DESCRIPTION
to find a acyclic directed path P such that the path P can
travel through all nodes in without repitition(no limit to
Here first with a simple example of the problem
node traversal order in ). Note: All weights are integers
description. As shown in Figure 1, assuming the source
within [1,30]; Edgs connecting node A and B are likely to
node number is 0; the intermediate nodes need to pass
be more than one, whose weight may be the same or may
through are 2,3; And the end node number is 1.
be different.
636
Authorized licensed use limited to: b:on Instituto Politecnico de Coimbra. Downloaded on April 04,2023 at 12:46:56 UTC from IEEE Xplore. Restrictions apply.
OSPF is an interior gateway protocols, based on link- nodes there are in the intermediate node set,the more
state algorithms.The process of routing computation by complex the problem will be), and then uses exhaustive
link-state algorithm is shown in Figure 2. RTA, RTB, RTC, search to list all the partial paths and selects one with
RTD are routers in a same area running OSPF, each router minimum path cost, making the chances of getting the
in the area floods Router-LSA (router link-state optimal path as large as possible. The specific algorithm
adervertisement) to form a LSDB(link state database) of the flow chart is shown in Figure 3.
area. Each router uses the SPF algorithm to calculate a
shortest path tree with itself as the root based on the LSDB Start
and this shortest path tree describes the shortest path to all
other routers in the area.The SPF algorithm used here is Load data: s, t, V', V
Dijkstra algorithm. Therefore, learn from the above process
of routing computation, the first step of our algorithm to be
Set inspired hop count: N
given is to calculate the shortest path tree of the nodes in
and the starting node according to the given graph by
using Dijkstra algorithm. To get SPT of nodes in
V' and source node s
LSDB
Is there a node in V' Y
RTA RTB LSA of RTA No solution
unreachable?
RTC LSA of RTB
Flooding
Router-LSA N
LSA of RTC
Calculate the cost of all
RTD LSA of RTD paths with top N hops
N
Computing Path exists?
Routes
IP Routing Table Y
SPT
Figure 2. Routing computation of OSPF Take the first -hop nodes
of the shortest path
Not only the distance to the source node, but also the as the next -hop nodes
distance to the termination node of each node are taken into
consideration in A* algorithm. If g(n) is used to represent N Will reach
the path cost from the starting node to any other node n, the the termination node after N
h(n) represents the estimated path cost(heuristic value) of hops?
the node n to the termination node,define that:
Y
f(n) = g(n) + h(n)
A* algorithm calculates each node’s value of f(n) to find N=1? N=N-1
a minimum one in the path searching before reaching the Y
termination node. So, learn from the above process, our
algorithm to be given will also calculate each node’s value End
of the partial path cost based on the SPT acquired in the Figure 3. Flow chart of the proposed algorithm
first step to find the minimum one as the next node.
Data to be loaded includes a source node s (assume
3.1. Algorithm Description node number is 0), the termination node t (assume node
The algorithm to solve the path planning problem with a number is T-1), an intermediate node set to be passed
specified intermediate node set is mainly divided into two (assume that there are total S nodes whose node number are
steps. The first step is to focus on the source node s and the 1,2, ... S respectively, ie, = {1,2, ... S}), all node set
nodes in the intermediate node set, using the Dijkstra (assume that there are totally T nodes); Assuming that the
algorithm to calculate the shortest path tree of each node inspired hop count N is set to be 3 according to the loaded
above(each node’s SPT describes the minimum path cost data; Firstly, to calculate SPT of the nodes in and the
from itself, the root node, to all other nodes ). The second node s (total +1 nodes)using Dijkstra algorithm, and judge
step is on the basis of the first step, using heuristic search that whether there exist nodes in which are not included
theory, sets the appropriate "inspired hop count" according in the above shortest path tree, if do exist, then the path
to the computational complexity of the problem(the more meeting the requirements can’t be find; If do not exist, then
637
Authorized licensed use limited to: b:on Instituto Politecnico de Coimbra. Downloaded on April 04,2023 at 12:46:56 UTC from IEEE Xplore. Restrictions apply.
start to calculate. Because the inspired hop count N = 3,
that is, the partial optimal path of each calculation covers 3
0 2 3
intermediate nodes in , so the first step of the calculation 0,1
is as follows: 1,2 2,1 3,3 5,1 4,1 6,1
= + +··· + 1
When = 3 :
= + + 2 3 1 3 1 2
Where represents the path weights, 0 、 1 、 2 、 Figure 4. SPT of the node 0,2,3 in the topology graph shown in figure1
638
Authorized licensed use limited to: b:on Instituto Politecnico de Coimbra. Downloaded on April 04,2023 at 12:46:56 UTC from IEEE Xplore. Restrictions apply.
expressed as ( ∗ S!) when N = S, which is easy to see the longest time. However, when there are more nodes in
the time complexity of the algorithm is not sensitive to , the intermediate node set(more than 20), we can no longer
the total number of nodes in V, but very sensitive to S ,the set N = S(S ≥ 20).
number of nodes in intermediate node set V .If we take
N (the insired hop count) into consideration, the time 22
!
complexity should be expressed as ( ∗ ( )!).
11 21
3.4. Program Simulation
639
Authorized licensed use limited to: b:on Instituto Politecnico de Coimbra. Downloaded on April 04,2023 at 12:46:56 UTC from IEEE Xplore. Restrictions apply.
3.4.2 Experiment 2: The effects of parameter on the path network in 3D," 2014 International Conference on Data Mining
weight and search time of the results and Intelligent Computing (ICDMIC), New Delhi, pp. 1-5, 2014.
In experiment 2, we use a more complex topological graph [2]Y. Liu, R. Song and R. Bucknall, "A practical path planning
with 20 intermediate nodes to test our algorithm. The test and navigation algorithm for an unmanned surface vehicle using
results are graphically presented in Figure 7. the fast marching algorithm," OCEANS 2015 - Genova, Genoa,
pp. 1-7, 2015.
[3]Chunbao Wang et al., "Path planning of automated guided
vehicles based on improved A-Star algorithm," 2015 IEEE
International Conference on Information and Automation,
Lijiang, pp. 2071-2076, 2015.
[4]J. Lee, B. y. Kang and D. w. Kim, "Fast genetic algorithm for
robot path planning," in Electronics Letters, vol. 49, no. 23, pp.
1449-1451, Nov. 7 2013.
[5]T.A. Bahtiyar, F. Ardilla, B.S. Marta and D. Pramadihanto,
"Effectiveness of bicycle path planning method and pure pursuit
method on omni-directional mobile robot," 2015 International
Conference on Control, Electronics, Renewable Energy and
Communications (ICCEREC), Bandung, pp. 91-97, 2015.
[6]L. Cheng, C. Liu and B. Yan, "Improved hierarchical A-star
algorithm for optimal parking path planning of the large parking
Figure 7. The effects of parameter N on the results lot," 2014 IEEE International Conference on Information and
Automation (ICIA), Hailar, pp. 695-698, 2014.
From the results of experiment 2, we can see that with the [7]J.R.Jiang, H.W.Huang, J.H.Liao and S.Y.Chen, "Extending
increase of the value of N , the search time increases Dijkstra's shortest path algorithm for software defined
exponentially, while the weight of the path is getting lower networking," 2014 16th Asia-Pacific Network Operations and
and lower. Therefore, the actual value of N depends on Management Symposium (APNOMS), Hsinchu, pp. 1-4, 2014.
your requirements on search time and path weight [8]Mingjun Wei and Yu Meng, "Research on the optimal route
respectively. choice based on improved Dijkstra," 2014 IEEE Workshop on
Advanced Research and Technology in Industry Applications
4. CONCLUSION (WARTIA), Ottawa, ON, pp. 303-306, 2014.
[9]F. H. Tseng, T. T. Liang, C. H. Lee, L. D. Chou and H. C. Chao,
In this paper, we have studied a class of complex path "A Star Search Algorithm for Civil UAV Path Planning with 3G
planning problem: In the weighted directed graph, paths Communication," 2014 Tenth International Conference on
need to be found which start from a specified node, after Intelligent Information Hiding and Multimedia Signal Processing
passing through some specified intermediate nodes without (IIH-MSP), Kitakyushu, pp. 942-945, 2014.
repetition (the number of intermediate nodes is uncertain), [10]Tan Baocheng, Wang Pei, Improvement and Implementation
and finally reach the termination node. We expect to find of A* Path Planning Algorithm[J]. Journal of Xi'an Technological
such a path with a minimum path weight within a certain University, 04:325-329, 2012.
time. To give the solution, this paper took the Dijkstra
[11]J. Yao, C. Lin, X. Xie, A. J. Wang and C. C. Hung, "Path
algorithm as the foundation, and drew on the heuristic
Planning for Virtual Human Motion Using Improved A* Star
search theory of A* algorithm, finally proposed a feasible Algorithm," 2010 Seventh International Conference on
algorithm. In this algorithm, we put forward a new concept Information Technology: New Generations (ITNG), Las Vegas,
of "inspired hop count" which is the most important NV, pp. 1154-1158, 2010.
parameter of the algorithm to dynamically adjust the
[12]Yu Haicong, Lu Feng. A Multi-modal Multi-criteria Route
balance between the path weight and the searching time.
Planning Method Based on Genetic Algorithm[J]. Journal of
Finally, the experimental results were given to show that Surveying and mapping, 01:89-96, 2014.
this algorithm can effectively solve the above path planning
problem and analyse the influence of parameter [13]Shi Enxiu, Chen Minmin, Li Jun, Huang Yumei. Research on
N (‘inspired hop count’) on the performance of the Method of Global Path-planning for Mobile Robot Based on Ant-
colony Algorithm[J]. Transactions of the Chinese Society for
algorithm.
Agricultural Machinery, 06:53-57, 2014.
5. ACKNOWLEDGMENTS [14] M. Ryan, "Constraint-based multi-robot path planning," 2010
IEEE International Conference on Robotics and Automation
This research was partially supported by the National (ICRA), Anchorage, AK, pp. 922-928, 2010.
Nature Science Foundation of China (No.61373084). [15]Huang Shuli, Hu Dasha, Jiang Yuming. Algorithm for finding
shortest path which go through specified intermediate node set[J].
6. REFERENCES Computer Engineering and Applications, 11:41-46, 2015
640
Authorized licensed use limited to: b:on Instituto Politecnico de Coimbra. Downloaded on April 04,2023 at 12:46:56 UTC from IEEE Xplore. Restrictions apply.