Abstract:
In providing a method for finding a shortest path from a starting place to a destination place in a traffic network including turn restrictions, U-turns and P-turns using a Dijkstra algorithm and a Floyd-warshall algorithm, the Dijkstra algorithm includes the steps of: assigning a virtual arc connection value from a starting node to a destination node based on traffic information in order to do the turn restriction, the U-turn and the P-turn, wherein the starting node indicates the starting place and the destination node indicates the destination place; selecting a smallest travel time value out of total travel time value for a temporary label node from the starting node to all nodes except for the starting node and assigning the smallest travel time value to a permanent label node; and determining the shortest path by tracing a permanent node stating from the destination node. The Floyd-warshall algorithm includes the steps of: computing a travel time value between two nodes for all pairs of nodes; if a left-turn restriction is included in continuous three nodes, adding a virtual node L′ in the continuous three nodes, wherein let the three nodes be nodes V, L and W; computing matrix of the pairs of nodes based on the virtual node L′; and finding a smallest travel time using the matrix of the pairs of nodes.

Description:
FIELD OF THE INVENTION  
         [0001]    This invention relates to a method for finding a shortest path and a computer-readable record media storing instructions for performing the method; and more particularly, to a method for finding a shortest path in a traffic network using a Dijkstra algorithm and a Floyd-warshall algorithm and a computer-readable record media storing instructions for performing the method.  
         DESCRIPTION OF THE PRIOR ART  
         [0002]    A Dijkstra algorithm and a Floyd-warshall algorithm are ones, which are frequently used for finding a shortest path.  
           [0003]    The conventional Dijkstra algorithm makes it possible to find a shortest travel time value from a starting node to all other nodes in a network with nonnegative arc values. According to this algorithm, all nodes are classified into two groups of node. One is a permanent node that represents a shortest travel time value from a starting node to a relevant node. The other one is a temporary node that represents an upper bound of a real shortest travel time value from the starting node to the relevant node. Also, a label of a node refers to a shortest travel time value from the starting node along a path whose internal nodes are all permanently labeled. Initially, the starting node is labeled with the permanent node whose travel time value is zero.  
           [0004]    According to the Dijkstra algorithm, a node having a minimum temporary label is selected to thereby become the permanent node. The Dijkstra algorithm terminates when all nodes including from the relevant node to neighboring nodes become the permanent node.  
           [0005]    The Floyd-warshall algorithm makes it possible to find a shortest path between all pairs of nodes in a network.  
           [0006]    According to the Floyd-warshall algorithm, a matrix of π ij  is computed. The π ij  denotes a travel time value for each of all pairs of nodes, i and j. A π ij   (m)  ( denotes a matrix of a shortest travel time value from the node i and to the node j using 1, 2, . . . , m−1 as the internal node. A π ij   (k)  denotes a matrix of a travel time value for each of all pairs of nodes i and j within a starting node and an arrival node. A π ij   (k+1)  is computed based on the π ij   (k) . Here, a Triple operation needs to be performed in order to compute the π ij   (k+1) .  
           [0007]    [0007]FIG. 7 shows one embodiment of a Triple operation according to a Floyd-warshall algorithm.  
           [0008]    Here, it is supposed that a π ij   (m)  denotes a matrix for a shortest travel time value from the node i and to the node j using 1, 2, . . . , m−1 as an internal node. A π ij   (k+1)  is generated by calculating a given π ij   (k)  on a following Triple operation.  
           [0009]    That is to say, as for a shortest travel time using nodes 1, 2, . . . , k as an internal node, if π ij   (k) &lt;π ik   (k) +π kj   (k) , π ij   (k) =π lj   (k+1)  without passing a node k , however if π ij   (k) &gt;π ik   (k) +π kj   (k) , π ij   (k) =π ik   (k) +π kj   (k)  with passing the node k.  
           [0010]    When representing the above by a mathematical equation, the equation is as follows:  
           π ij   (k+1) =min {π ij   (k) , π ik   (k) +π kj   (k) } 
           [0011]    In other words, a less one of π ij   (k)  and π lk   (k) +π kj   (k)  is assigned to π ij   (k+1) .  
           [0012]    As described above, the conventional Dijkstra algorithm and Floyd-warshall algorithm are used to find a shortest path, however there is a problem that the same node should not be passed twice. Also, there is another problem that the above two algorithm cannot be applied to a real traffic network, in which a Left-turn restriction, a U-turn and a P-turn are included.  
         SUMMARY OF THE INVENTION  
         [0013]    It is an object of the present invention to provide a method for finding a shortest path from a starting place to a destination place in a traffic network to thereby reach the destination place in shorter time.  
           [0014]    In accordance with an aspect of the present invention, there is provided a method for finding a shortest path from a starting place to a destination place in a traffic network including one or more turn restrictions, one or more U-turns and one or more P-turns using a Dijkstra algorithm, the method including the steps of: assigning a virtual arc connection value from a starting node to a destination node based on traffic information in order to do the turn restriction, the U-turn and the P-turn, wherein the starting node indicates the starting place and the destination node indicates the destination place; selecting a smallest travel time value out of total travel time value for a temporary label node from the starting node to all nodes except for the starting node and assigning the smallest travel time value to a permanent label node; and determining the shortest path by tracing a permanent node stating from the destination node. 
       
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS  
       [0015]    Other objects and aspects of the invention will become apparent from the following description of the embodiments with reference to the accompanying drawings, in which:  
         [0016]    [0016]FIG. 1 shows a configuration of a DRGS (dynamic route guidance system) in accordance with the present invention;  
         [0017]    [0017]FIG. 2 shows a configuration of a traffic network including a turn-restriction;  
         [0018]    [0018]FIG. 3 is a configuration of a traffic network including a turn-restriction and a U-turn;  
         [0019]    [0019]FIG. 4 is a configuration of a traffic network including a turn-restriction, a P-turn and a U-turn to which a Diijkstra algorithm is applied;  
         [0020]    [0020]FIGS. 5A and 5B are flow charts illustrating one embodiment of an improved Dijkstra algorithm in accordance with the present invention;  
         [0021]    [0021]FIG. 6 is a configuration of a traffic network including a turn-restriction, a U-turn and a P-turn to which a Floyd-warshall algorithm is applied;  
         [0022]    [0022]FIG. 7 shows one embodiment of a Triple operation according to a Floyd-warshall algorithm; and  
         [0023]    [0023]FIG. 8 is a flow chart illustrating one embodiment of an improved Floyd-warshall algorithm in accordance with the present invention. 
     
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS  
       [0024]    [0024]FIG. 1 shows a configuration of a DRGS (dynamic route guidance system) in accordance with the present invention.  
         [0025]    Reference numerals  100 ,  110 ,  120 ,  130 ,  140  and  150  denote respectively a real-time movement management system, an integrated logistics information network (ILIN), traffic information system, an operation center, a freight carrier and a traveler information system.  
         [0026]    The traffic information system  120  includes, as a traffic information source, police, a road administrator, road traffic information center or the like. The freight carrier  140  includes a vehicle having an OBU (on board unit). Each vehicle can find a desired location using the equipment for the GPS and an operator of the freight carrier can perform monitoring information about a location of the freight and the vehicle, a travel time, a vehicle status, some freight loaded on a specific vehicle using the same.  
         [0027]    The real-time movement management system  100 , the traffic information system  120 , the operation center  130 , the freight carrier  140  and the traveler information system  150  are all coupled with one another through the integrated logistics information network (ILIN)  110 .  
         [0028]    The DRGS performs guiding a vehicle to an appropriate route, wherein a driver can escape from a traffic congestion to thereby reach his/her destination earlier. Accordingly, a traffic flow can be smooth and vehicle movement efficiency can be increased.  
         [0029]    The DRGS needs to estimate a travel time along each path between crossroads and the estimation is based on data that are measured at each of a variety of time slots in the past time, wherein the data are collected using equipment for a GPS (global positioning system). In other words, the travel time estimation for next minutes is required because the traffic flow varies with short time. Also, when this estimation is compared with a next estimation for a traffic variable, it is possible to perform a proper traffic signal control and find an optimal route according to a status of the traffic flow.  
         [0030]    As shown in FIG. 1, the real-time movement management system (RTMMS)  100  provides the operation center  130  with a vehicle location, a travel time, vehicle information or the like via an uplink and receives traffic information such as a new shortest path corresponding to a current traffic status via a downlink.  
         [0031]    The traffic information system  120  receives information about an overall traffic from the road traffic information center, the police, the road administrator, a weather station or the like and provides the operation center  130  with the information.  
         [0032]    The freight carrier  140  provides the operation center  130  with information about a location of the freight and the vehicle, the travel time, a vehicle status, some freight loaded on a specific vehicle or the like via the integrated logistics information network (ILIN)  110  and receives information analyzed in the operation center  130  via the same.  
         [0033]    The traveler information system  150  provides the operation center  130  with information related to a route guidance service, for example, traffic information, regulations, a weather in the road and a traffic congestion via the integrated logistics information network (ILIN)  110  and receives information analyzed in the operation center  130  via the same.  
         [0034]    [0034]FIG. 2 shows a configuration of a traffic network including a turn-restriction.  
         [0035]    As shown in FIG. 2, a value on a corresponding path across two nodes represents a travel time value from one node to the other node. The Dijkstra algorithm in accordance with the present invention uses continuous three nodes including a previous permanent node, a current permanent node and a temporary node adjacent to the current permanent node in order to exclude the turn-restriction in the course of finding a shortest path. In other words, if a path along some continuous three nodes is identical to the turn-restriction path, the path cannot be the shortest path.  
         [0036]    For example, it is assumed that a vehicle is going to travel from a starting node ( 3 )  200  to a destination node ( 11 )  201 . If there are turn-restriction paths including a path  210  along nodes  6 - 5 - 8  and a path  211  along nodes  3 - 4 - 7 , a shortest path is a path along nodes  3 - 6 - 7 - 10 - 9 - 8 - 11  and an optimal total travel time along the shortest path is 14, wherein 14 equals a sum of travel times on each path, that is, 3+2+4+2+1+2. Consequently, the turn-restriction paths including the path  210  along nodes  6 - 5 - 8  and the path  211  along nodes  3 - 4 - 7  cannot be paths as a comparison object in order to find a shortest path.  
         [0037]    [0037]FIG. 3 is a configuration of a traffic network including a turn-restriction and a U-turn.  
         [0038]    An improved Dijkstra algorithm in accordance with the present invention is applied to the traffic network in FIG. 3.  
         [0039]    A virtual arc connection is used in consideration of the U-turn in applying the improved Dijkstra algorithm to the traffic network.  
         [0040]    As shown in FIG. 3, each value between nodes represents a travel time value from one node to the other node and the traffic network includes two turn-restriction paths including a path  310  along nodes  6 - 5 - 8  and a path  311  along nodes  3 - 4 - 7 . It is possible for a driver to do a U-turn along nodes  6 - 5 - 2 - 5 - 8  in order to travel from a node  6  to a node  8  and a U-turn along nodes  3 - 4 - 1 - 4 - 7  in order to travel from a node  3  to a node  7 .  
         [0041]    For the above two U-turns, the driver should pass the nodes  5  and  4  one more time, which the driver already passed. However, the U-turn cannot be used in the Dijkstra algorithm because the Dijkstra algorithm doesn&#39;t consider the U-turn path containing a node passed twice.  
         [0042]    In accordance with the present invention, virtual arcs including a path  300  from a node  2  to a node  8  and a path  301  from a node  1  to a node  7  and virtual arc values (the virtual arc value=travel time) including a value of 5 for the path  300  and a value of 3 for the path  301  wherein 5=2+3 and 3=2+1, are used to find a shortest path according to the improved Dijkstra algorithm.  
         [0043]    If each of the two virtual arc values 5 or 3 is respectively less than each of existing values from the node  2  to the node  8  or from the node  1  to the node  7 , each virtual arc value is used to find the shortest path. Otherwise, each existing value is used to find the shortest path instead of the virtual arc value.  
         [0044]    [0044]FIG. 4 is a configuration of a traffic network including a turn-restriction, a P-turn and a U-turn to which a Diijkstra algorithm is applied.  
         [0045]    As shown in FIG. 4, a P-turn path  401  is formed along nodes H-E-B-C-F-E-D and a U-turn path  402  is formed along nodes H-E-F-E-D.  
         [0046]    A virtual arc  403  from a node F to a node D is used to find a shortest path, wherein a node E is in the center of the virtual arc  403 .  
         [0047]    If there exists already an arc from the node F to the node D in the traffic network, a following consideration needs to be taken into. If a virtual arc value from the node F to the node D is less than an existing arc value, the virtual arc value is used to find the shortest path. In other words, if the travel time of the path from the node F to the node D via the node E is shorter than that of the other paths from the node F to the node D, the path  403  is selected to find the shortest path.  
         [0048]    [0048]FIGS. 5A and 5B are flow charts illustrating one embodiment of an improved Dijkstra algorithm in accordance with the present invention.  
         [0049]    As shown in FIGS. 5A and 5B, at the step S 501 , a starting node and a destination node are selected in a traffic network.  
         [0050]    At the step S 502 , it is performed to assign a virtual arc value based on traffic information for a U-turn. In this case, if there exists already an arc in the traffic network, followings should be considered. Comparing the virtual arc value with an existing arc value, if the virtual arc value is less than an existing arc value, the virtual arc value is used to find a short path. Otherwise, the existing arc value is used.  
         [0051]    At the step S 503 , it is performed to assign a d lj  to a π j   b . Here, the d lj  represents a travel time value from a starting node, usually a node  1  to a node j. The π j   b  represents a total travel time of a temporary label node from the starting node to the node j.  
         [0052]    At the step S 504 , it is performed to assign the starting node  1  to a b and if there exists no arc between the starting node  1  and the node j, d 1j =∞ (an infinite value). Here, the b represents a previous permanent node for the node j that is a temporary label node.  
         [0053]    At the step S 505 , it is performed to assign 0 to a π 1   * .  
         [0054]    A π j   *  represents an optimal total travel time value from the starting node to the node j. Accordingly, the π 1   *  represents a total travel time value from the starting node to the node  1 .  
         [0055]    At the step S 506 , it is performed to select a shortest travel time for every temporary label π j   b  assign a permanent label to the same. If a permanent label node is an i, a minimum value (min [π j   b ]) out of the total travel time of a temporary label node from the starting node to the node j (π j   b ) is selected to thereby be designated as a π i   b .  
         [0056]    At the step S 507 , it is performed to assign the π l   b  to a π l   * . The π i   *  represents a total travel time value from the starting node to the node i.  
         [0057]    At the step S 508 , it is performed determine if there exists no more temporary label, and if not, the logic flow proceeds to step S 513 , otherwise proceeds to step S 509 .  
         [0058]    At the step S 509 , it is performed to modify all temporary nodes adjacent to the permanent node i.  
         [0059]    At the step S 510 , if a path along nodes b-i-j is not the turn-restriction path, it is performed to compare the total travel time of the temporary label node from the starting node to the node j (π j   b ) with a sum of the total travel time value from the starting node to the permanent node i (π l   * ) and a travel time value from the permanent node i to a temporary node j adjacent to the i (d lj ) and then as a result of the comparison, select less travel time to thereby assign the less travel time (min [π j   b , π i   * +d lj ]) to a new total travel time of the temporary label node from the starting node to the node j (π j   b ). Here, if the path along nodes b-i-j is the turn-restriction path, the b-i-j path is excluded from the shortest path.  
         [0060]    At the step S 511 , it is performed to determine if the new total travel time of the temporary label node from the starting node to the node j (π j   b ) is larger than the sum of the total travel time value from the starting node to the permanent node i (π i   * ) and a travel time value from the permanent node i to a temporary node j adjacent to the i (d ij ), and if not, the logic flow proceeds to step S 513 , otherwise proceeds to step S 512 .  
         [0061]    If π j   b &gt;π i   * +d ij , at the step S 512 , it is performed to assign the previous permanent node i to b for π j   b  and then the logic flow returns to the step S 506 .  
         [0062]    At the step S 513 , it is performed to find a shortest path by tracing b representing a previous permanent node for the current permanent node, starting from the destination node.  
         [0063]    At the step S 514 , it is performed to trace an actual path after eliminating a virtual arc for the U-turn and the P-turn.  
         [0064]    [0064]FIG. 6 is a configuration of a traffic network including a turn-restriction, a U-turn and a P-turn to which a Floyd-warshall algorithm is applied.  
         [0065]    According to an improved Floyd-warshall algorithm, in case a turn restriction arc is included in a traffic network, the turn restriction arc is excluded from objects for a shortest path by assigning a significantly large value of a travel time, for example, an infinite value (∞) to the arc.  
         [0066]    As shown in FIG. 6, the traffic network includes a Left-turn restriction path  601  along nodes V-L-W, two U-turn paths including a path along nodes V-L-a 3 -L-W and a path along nodes V-L-a 5 -L-W and a P-turn path along nodes V-L-a 3 -a 4 -a 5 -L-W. The node L  602  is in the center of the traffic network.  
         [0067]    When assuming that a starting node is a node S  604  and a destination node is a node T  605 , the Left-turn restriction path  601  along the nodes V-L-W needs to be considered.  
         [0068]    To find a shortest path, it is required to obtain a first shortest travel time matrix (π ij   (k) ) by considering two paths including a path along nodes V-a 6 -a 5 -L-W and a path along nodes V-a 6 -a 5 -a 4 -a 3 -L-W as if an arc along nodes V-L doesn&#39;t exist, that is, π VL   (0) =∞ and also to obtain a second shortest travel time matrix (π ij   (k) ) by considering two paths including a path along nodes V-L-a 3 -a 2 -W and a path along nodes V-L-a 5 -a 4 -a 3 -a 2 -W as if an arc along nodes L-W doesn&#39;t exist, that is, π LW   (0) =∞.  
         [0069]    It is required to compare the first shortest travel time matrix (π ij   (k) ) with the second shortest travel time matrix (π ij   (k) ) to find the shortest path in the traffic network including the Left-turn restriction path  601  along the nodes V-L-W.  
         [0070]    A major disadvantage on the above approach is that a computation time and a memory capacity are increased 2 n  times as late and big, wherein the n of the 2 n  denotes a number of the Left-turn restriction. Accordingly, the improved Floyd-warshall algorithm should be applied to the traffic network including a limited number of the Left-turn restriction.  
         [0071]    The U-turn and the P-turn cannot be used to find an optimal path in a conventional shortest path algorithm because the conventional algorithm doesn&#39;t consider the U-turn and the P-turn containing a node passed twice.  
         [0072]    As shown in FIG. 6, it is required to add a virtual node L′  603  in the traffic network. Consequently, it is possible to consider the P-turn along the nodes V-L-a 3 -a 4 -a 5 -L′-W and the two U-turns including the path along the nodes V-L-a 3 -L′-W and the path along the nodes V-L-a 5 -L′-W, each of which don&#39;t pass the same node twice.  
         [0073]    Any virtual node can be added in a traffic network including any P-turns or U-turns, wherein a node needs to be in the center of the P-turns or the U-turns.  
         [0074]    [0074]FIG. 8 is a flow chart illustrating one embodiment of an improved Floyd-warshall algorithm in accordance with the present invention.  
         [0075]    At the step S 800 , a starting node and a destination node are selected in a traffic network.  
         [0076]    At the step S 801 , it is performed to compute a travel time value between two nodes (π lj ) for all pairs of nodes. If i=j, a π ii  equals 0. That is to say, the travel time value between two nodes (π lj ) at the same node is 0. Also, if an arc doesn&#39;t exist between two nodes, the travel time value between two nodes (π ij ) equals an infinite value (∞).  
         [0077]    At the step S 802 , if a path along nodes i-L-j is a Left-turn restriction when i=V and j=W, it is performed to add a virtual node L′ in an initialization matrix π (0) .  
         [0078]    At the step S 803 , it is performed to compute a first matrix π (0) , assuming that π VL   (0) =∞, π VL′   (0) =∞ and compute a second matrix π (0) , assuming that π LW   (0) =∞, π L′W   (0) =∞.  
         [0079]    At the step S 804 , it is performed to indicate a row m and a column m in each matrix π (m−1) .  
         [0080]    At the step S 805 , it is performed to compute each π (m)  for the first matrix and the second matrix using a Triple operation.  
         [0081]    At the step S 806 , if π VW   (m−1) &gt;π Vm   (m−1) +π mW   (m−1) , (V≠m and W≠m), then π VW   (m) =π Vm   (m−1) +π mW   (m−1) .  
         [0082]    At the step S 807 , it is performed to determine if the above m equals N+1. The m represents a sum of the number of existing nodes and the number of virtual nodes, and if not, the logic flow returns to the step S 804 , otherwise proceeds to step S 808 .  
         [0083]    If the above m equals N+1, at the step S 808 , it is performed to compare each π (m)  to find a π ST   (m) . The π ST   (m)  represents a shortest travel time between the starting node S and the destination node T for the first matrix and the second matrix.  
         [0084]    At the step S 809 , it is performed to find the π ST   (m) .  
         [0085]    It is possible to store instructions for performing the present invention in accordance with the present invention in a computer-readable record media.  
         [0086]    There is an advantageous effect that a driver can reach his/her destination in shorter time along an optimal shortest route in a real-time dynamic route guidance system (DRGS) since techniques of a virtual arc connection, a virtual node or the like are applied to a plurality of road situations of a real traffic network including a turn restriction, a U-turn and a P-turn based on improved Dijkstra algorithm and Floyd-warshall algorithm in accordance with the present invention.  
         [0087]    Also, there is another advantageous effect that it is possible to lower a computation time and memory capacity that are required to find a shortest route based on the improved Dijkstra algorithm in accordance with the present invention.  
         [0088]    Although the preferred embodiments of the invention have been disclosed for illustrative purpose, those skilled in the art will appreciate that various modifications, additions, and substitutions are possible, without departing from the scope and sprit of the invention as disclosed in the accompanying claims.