Research Article  Open Access
Distributed Approach for Solving TimeDependent Problems in Multimodal Transport Networks
Abstract
This paper presents an alternative approach for timedependent multimodal transport problem. We describe a new graph structure to abstract multimodal networks, called transfer graph, which adapts to the distributed nature of real information sources of transportation networks. A decomposition of the Shortest Path Problem in transfer graph is proposed to optimize the computation time. This approach was computationally tested in several experimental multimodal networks having different size and complexity. The approach was integrated in the multimodal transport service of the European Carlink platform, where it has been validated in real scenarios. Comparision with other related works is provided.
1. Introduction
Multimodal transport is the combination of two or more means of transport to move passengers or goods from one source to a destination [1]. The traveller can use either public (e.g., bus, taxi, metro, railway, and ship) or private vehicles (e.g., car, motorbike, bicycle, and walking). Multimodal journeys are becoming a real necessity in our society in order to guarantee a high level of mobility both inside cities and at the regional level [2].
From the last decades, many research projects were devoted to develop multimodal transport systems that recommend travellers a combination of transport means for doortodoor journeys [1, 3–5]. These systems must cope with different transport information sources, which are normally distributed in real scenarios, that is, transport data are usually provided by different public and private companies and can be stored in different locations. Another challenge for these systems is to integrate realtime traffic information (e.g., traffic jams, delays, roadworks) in order to adapt proposed routes to real conditions.
Different abstractions have been proposed to represent the timedependent multimodal transport problem. Three main categories of approaches can be identified: the multigraphbased approach [2, 6, 7], the constraint satisfaction problem approach [8], and the gridbased approach [9]. In order to adapt multimodal transport abstraction to the distributed nature of realworld information, the following constrains should be fulfilled: (i) the underlying multimodal network is assumed to be flat and to contain no hierarchical structure, like regional hierarchy; (ii) the involved unimodal networks may be distributed, that is, unimodal networks are stored and accessed separately; (iii) if there are multiple network information sources within a single mode, they may be distributed. The abovementioned approaches relax any of these constraints when modeling and solving the problem.
An alternative representation of the timedependent multimodal transport problem was partially described in [10]. Here we present a formal description of this abstraction as well as the transfer graph approach, which calculates the best paths in this structure. This approach was tested and also validated in the real context of Carlink Platform (Carlink, Wireless Platform for Linking Cars, is a project of the Celtic Cluster Programme Call, which aims to develop an intelligent wireless traffic service platform between cars) [11]. Transfer graph approach provides a solution for the timedependent multimodal transport problem which fulfills the three defined constraints and adapts to the distributed nature of real transport information providers.
The outline of this article is the following. Section 2 describes the transfer graph approach for the timedependent multimodal transport problem. In Section 3 this approach is computationally tested. Section 4 describes Carlink, the real scenario where the approach is validated, mainly the service platform as well as the obtained results. In Section 5 the results are compared with other related works. Finally, this paper ends with some conclusions in Section 6.
2. Transfer Graph Approach for Multimodal Transport Problems
In real scenarios, transport information sources are normally distributed, that is, transport schedules are usually provided by different public and private companies. Besides, this information is frequently modified due to daily unexpected perturbations such as traffic jams, delays, and roadworks. These restrictions complicate the management of transport information and the updates of transport networks.
This section presents the transfer graph approach for solving the Shortest Path Problem (SPP) in discrete timedependent multimodal network. It is able to compute multimodal routes even if all involved networks are kept structurally separated and distributed. Instead of building a solution limited to the unified view of multimodal transport system, our strategy consists in finding a fundamental network representation which covers multimodal networks and then solves the SPP on this abstraction.
The transfer graph approach adapts to the distributed nature of real world transport information sources and allows updates on each transport network independently without requiring any further modification. Our approach is mainly divided into two steps. In the first step, relevant paths are computed within each unimodal network. In the second step, the graph is reduced to a simple timedependent unimodal network, where the classical SPP can be solved.
The next subsections first introduce the basic concepts of graphs and a definition of SPP in timedependent multimodal networks. Then, the transfer graph abstraction is formally described, followed by the presentation of the relevant graph, a simplified structure calculated from transfer graph. Finally, the transfer graph approach to calculate SPP in transfer graph is discussed.
2.1. SPP in Discrete TimeDependent Multimodal Networks
Let denote a directed graph or network, in which is a set of nodes, is a set of transport modes (e.g., train, bus, and car), and is a set of arcs. An arc can be identified by or , where . The means the possibility of going from node to by using transport mode . A value is associated to each arc , indicating the cost of including the arc in the solution (e.g., distance, time).
Definition 2.1. A graph is said to be multimodal if there is at least two transport modes where , , and . If there is only one transport mode in the graph, the graph is said to be unimodal.
In a graph a path or route is a sequence of arcs between a pair of nodes and , , for all , and , for all .
Definition 2.2. Given a graph , a path is said to be multimodal if , , and . If there is only one mode involved in the path, then the path is said to be unimodal.
Let be the set of paths in a graph and a function which represents the cost of a path. An SPP is considered to be multiobjective if maps to a vector of dimension , . Each component , for all , follows the definition of single objective path cost function.
Let denote time in a discrete set . Unlike nontimedependent networks where arc costs are fixed, in timedependent networks the cost of an arc varies according to cost function for all .
Definition 2.3. Given a graph and a discrete time set , a travel for the arc is defined as tuple , , where is the departure time from node , and is the arrival time at node . Tr denotes the set of travels in .
Definition 2.4. Given a graph and a discrete time set , is said to be a timedependent network if for all , where and for all , and .
We conclude this subsection by formulating the SPP in timedependent multimodal graphs. The reader can find a wider description of the SPP in [12].
Definition 2.5. Given a timedependent multimodal graph , two nodes , and time , the SPP is defined as follows: calculate a path from node to that departures from at time and is minimal.
2.2. Transfer Graph
In this subsection, we present the transfer graph, a graphical structure that abstracts the timedependent multimodal transport network.
A transfer graph is described by a set of unimodal networks and a set of arcs connecting them, and it is defined by . Let be an origindestination pair and a set of unimodal networks called components, where each , and . itself is such that and , for all , . Unlike partitioned graph, given two distinct components and , having is not mandatory; however must always hold. Again, let us define two distinct components ; a node is said to be a transfer point when , that is, any node which belongs to more than one component. Finally, a transfer is any arc which connects two components of the transfer graph: , .
Figure 1 illustrates an example of this abstraction, where , , and are components connected by three transfers, and nodes , , and are transfer points. One can observe that source node and target node may belong to more than one component; this is because they can be transfer points. Actually, timedependence and multiobjectivity are not considered in this example, although the approach is able to work with these restrictions.
Transfer graph abstraction separates all transport modes in unimodal networks. So, when a perturbation (e.g., traffic jam, roadworks) is detected for a transport mode, it does not affect globally to the representation, but locally. In those cases, it is only necessary to update those components related to the incident. According to this characteristic we can state Lemma 2.6.
Lemma 2.6. Transfer graph abstraction is more robust against modification of transport network data than classical multimodal network representation.
2.3. Relevant Graph
In the previous subsection we have formally defined the transfer graph. We present now a simplified structure derived from this abstraction, called relevant graph, from which the Shortest Path Problem can be easier computed.
At a high level, paths in a transfer graph can be divided into two groups: intracomponents paths and intercomponent paths. An intercomponent path within a transfer graph is considered as any sequence of arcs which connects two nodes where at least there are two arcs belonging to two distinct components. On the other hand, an intracomponent path within is any sequence of arcs which connects two nodes whose arcs belong only to one component. Intra component paths can be divided into the following categories:
(i)full paths (): the smallest route which starts at source node and ends at target node within the component ;(ii)relevant head paths (): the smallest routes which start at source and end at a transfer point within the component ;(iii)relevant intermediate paths (): the smallest routes which start and end at any transfer point within the component ;(iv)relevant tail paths (): the smallest routes which start at any transfer point and end at target within the component .Given a transfer graph and an origindestination pair , assume that for all components we have computed the best relevant path sets: , , , and (see Table 1). Having these relevant path sets in hand, it is possible to derive a reduced graph from which all possible best inter component paths from to can be computed. We call this graph the relevant graph and use to denote it. The nodes of are , and a subset of all transfer points in the transfer graph. The arc set of is the best relevant intra component paths viewed as edges, that is, paths that are represented only by the initial and final nodes. In fact, is a multigraph, but in general, its scale should be much smaller than the equivalent multigraph of the underlying transfer graph. This is because any node appearing in is a transfer point, except the origin and destination nodes. So, the number of nodes in depends directly on the number of transfer points in the whole transfer graph. Figure 2 shows the relevant graph computed out of Figure 1.

In order to formally describe , let denote the set of relevant nodes for component , and let be the set of relevant paths computed for . From all and we can, respectively, build for all the set of relevant nodes for all component, and for all the set of relevant paths for all component . The relevant graph is formally described by the pair . If we consider timedependent networks with different discrete time , then in order to build the relevant graph, it is necessary to calculate the best relevant paths at any departure time , that is, , , , , for all and for all .
When a relevant graph is built from a transfer graph for nodes and and time , we have a timedependent unimodal multigraph with a reduced number of nodes where all the shortest paths between transfer points, and , are calculated. The computation of the shortest path from to in is basically the classical origindestination SPP on a reduced timedependent unimodal network. This allows us to establish the following lemma.
Lemma 2.7. Let be a transfer graph and its associated relevant graph, then is a less complex structure than in terms of the shortest path computational time.
2.4. Transfer Graph Approach
Once the transfer graph and the relevant graph are formally defined, we present the decomposition of the SPP in these abstractions as well as the transfer graph approach.
Every time a request is fulfilled for obtaining the shortest path from a source and a destination at departure time , we need to calculate the best relevant paths for , , and for all , for all . However, for all is calculated once, since it does not depend on and . Taking advantage of this feature we can decompose the SPP in transfer graphs as follows.
(1)The best paths for , for all , for all are calculated and stored. This is called precalculations.(2)When a request is fulfilled for two specific nodes and , and time , then , , and for all can be calculated.(3)Finally, the relevant graph is constructed and the shortest path is computed.This decomposition reduces time complexity of the problem considerably because precalculations compute the most difficult set of paths in terms of computation time. Precalculations provide an important advantage if multiple and frequent requests are fulfilled (e.g., in a real web service). On the contrary, precalculation must be partly recomputed when a perturbation is detected. Table 2 shows an example of recalculating the relevant graph for Figure 1 if the cost of the arc between nodes and for component increases from two to four. As it can be seen, only few paths are necessary to be recalculated, instead of considering all transport network. This is specially efficient for transfer graphs with multiple components.

Finally, we present the transfer graph approach in Algorithm 1. The procedure Compute calculates the shortest path for a specific component and departure time. This step can be calculated independently for each component, due to the characteristics of transfer graph and the decomposition of the SPP in this abstraction. The procedure StoreDataBase saves in the data base the previous calculated set of the shortest paths. The procedure Compute is called three times to calculate all the specific sets of shortest paths when the user request is submitted. RelevantGraph is the procedure reducing the initial graph to a more compact structure corresponding to a given user request. Finally, the procedure ShortestPath gives the smallest route requested by the user. Obviously all these procedures can be implemented in different ways and giving rise to different complexities. Here, we show the feasibility of our approach by using a variant of Dijsktra's algorithm. Therefore, the transfer graph approach can be applied to scenarios with distributed transport information sources.

3. Simulation Results
So far, we have introduced the transfer graph abstraction and the transfer graph approach. In this section we present a computational test of this approach with several instances of the SPP in timedependent multimodal networks.
We have selected eighteen instances of the SPP with three, four, and five transport modes and increasing number of nodes from 50 to 2000, named . A transfer graph is built for each instance with variable number of arcs and transfers, but keeping high number of transfers and 10 possible random travel times per each arc, in order to simulate real scenarios. Table 3 presents for each problem instance a description of the transfer graph, the characteristic of the relevant graph built from the transfer graph, and the computation time. Instances are grouped by number of nodes.

A variant of Dijkstra's algorithm [13] is used to implement the function of the transfer graph approach (Algorithm 1), which calculates the shortest path in terms of travel time. In these simulation tests we consider that for all , for all are available, that is, the precalculations are already computed. The computation time of precalculations is presented in Table 3.
The transfer graph approach was tested in a computer equipped with an Intel Core2 Duo T7300, 2 GHz, and 2 GB RAM. 100 tests have been done for each instance with random nodes and departure time, calculating the average and maximum CPU time for them all (see Table 3).
When building the relevant graph from the transfer graph, the number of nodes in all instances is reduced around 70%, although the number of arcs increases considerably (see Table 3). Since the complexity of Dijkstra's algorithm depends on the number of nodes , this is an appropriate algorithm to calculate the shortest path in a relevant graph.
Transfer graph approach has been tested in eighteen instances of the problem, providing a solution in a second approximately in average, even for instances with 2000 nodes and 60.000 travels. Precalculation time increases according to the network size and number of travels. Nevertheless, this precalculations must be done once, and considerably reduced the computation time for requests.
Finally, we present the computation cost of precalculations on mediumsized networks with different levels of complexity, according to number of travels. Five transfer graphs are built for the instance with variable number of travels. Figure 3 shows that the computation cost of precalculation is not considerably affected by the number of travels.
4. Experimental Validation
In previous sections the transfer graph approach has been described. In this section, we first present the Carlink Platform (http://carlink.lcc.uma.es/) as well as its services and constraints. Then, it is described how the transfer graph was validated in Carlink Platform for real scenarios.
The Carlink Wireless Traffic Service Platform is designed to provide a basis for a wide range of commercial and governmental traffic, safety, and transport services (see Figure 4). The platform is a wireless ad hoc type communication entity with connectivity to the backbone network via base stations. The communication platform itself is the key element of Carlink, but the various services created to the platform have an important role. The three most important services are the traffic local road weather service (RWS), the traffic management service (TMS), and multimodal transport service (MTS).
The local RWS is derived from the FMI (Finnish Meteorological Institute) road weather model [14] and provides automatic traffic safety enhancement operations by gathering safety related information from vehicles and delivering consolidated data back, automatically without requiring user initiation. The TMS provides the traffic status by using static traffic information offered by the authorities and the dynamic traffic information given by the Carlink floating car application. This application is composed of several cars collecting information about the vehicle density in different locations of the road network and providing information of the traffic status. Finally, the MTS calculates multimodal itineraries according to a travel requested and a set of constraints specified by the user. This service is based in the transfer graph approach which has been extensively described in Section 2.
The platform exploits information of the RWS and TMS, as well as of other services, in order to generate accident or incident warning from traffic and weather conditions (see Figure 4). The MTS uses this information to calculate best routes avoiding bad condition roads and accidental areas to reduce the traveling time, or even recommending to change the transport mode if a blocked road is detected (e.g., to park the car near a train station and taking the train instead).
4.1. Experimental Environment
The Wireless Traffic Service Platform is divided in three different parts: Traffic Service Central Unit (TSCU), Traffic Service Base Station (TSBS), and Mobile End User (MEU) with ad hoc connectivity and (noncontinuous) backbone network connectivity. The whole platform can be seen as data exchange architecture between TSCU and MEU, while TSBS acts as bidirectional data transceivers. MEUs communicate between them in an ad hoc manner when encountering and can also send critical data directly to the TSCU by using a lowercapacity communication method (GRPS). The service cores are external systems connected to TSCU.
TSCU achieves and processes uptodate traffic platform information from the service cores. The TSCU updates regularly TSBS with the more recent traffic platform information related to their area, through the fixed network. Each TSBS disseminates the traffic platform information to every MEU that passes close to it by using a wireless ad hoc connection. Finally, MEUs disseminate this information when encountering (see Figure 5).
MEU and TSBS are equipped with acquisition system that collect their own information (speed, GPSdata, airbag blast, road temperature, etc.) and transmit regularly to TSCU. This data is processed by the TSCU, and if a critical event is detected in an area (e.g., accident, bad road condition), a warning message is delivered in realtime directly to MEUs known to be in the warning area by using GPRS or similar communication method.
The advantage of this architecture is that both users (MEU) and disseminating systems (TSBS) send regularly information of the road condition to the TSCU. This is specially important when an accident occurs, because the TSCU is informed in realtime and can warn drivers in the surrounding of the accident area. Besides, since cars communicate each other when encountering, the number of TSBS necessary in roads can be considerably reduced, and also the system cost.
4.2. Experimental Results in Real Scenario
The validation of the developed system was performed in a real scenario (see Figure 6). The route origin is the Belgian city of Arlon, and the destination is Luxembourg city (Luxembourg). Several thousands of commuters do daily this trip by car, train, bus, or by a combination of these transport modes.
A transfer graph abstracting this scenario was deployed in the Carlink Platform and made available for the MTS. The transfer graph has 73 nodes, 49 transfer points, and 226 arcs and contains three components, associated to train, car, and bus. We consider that changing between two transport modes is done by walking, and this cost is represented in transfers.
The transfer graph approach is implemented in the MTS of the Carlink Platform. The user requests are sent to the MTS through a client application for mobile devices. In this application, the user can select the travel preferences (e.g., transport modes, departure time, origin, and destination), send the request to the MTS, and get the optimal path from the Carlink Platform (Figure 7).
In our scenario, the selected preferences are optimizing travel time to go from Arlon to Luxembourg by using any kind of transport mode and departing at 5:50. The first route recommended by the MTS consists in going directly by car using the highway (Table 4). While we are traveling by car, the TMS is reported at 6:00 with a perturbation in the current itinerary. Since this perturbation modifies the transport data for road transports, precalculations are recomputed for components bus and car, but not for rail. This perturbation increases the travel time for route 1 up to 80 minutes. Automatically, the MTS recalculates a second route taking into account this perturbation and the current GPS position of the traveller. This information is reflected to the user by the client application. The second route is expected to arrive 25 minutes before the first route and its itinerary steps are the following: park near the next train station in Kleinbettingen, take the next train to Luxembourg, and finally the bus to the final destination, as indicated in Table 4.

5. Related Work
The multimodal transport problem under the timedependence constraint is not often studied in literature. As far as we know there are two works related to timedependent multimodal transport problem which present the performance of their approach on different instances of the problem [2, 5].
Authors in [5] abstract the problem with nonhierarchical multimodal graph and propose several data structures and tables, previously precalculated, in order to optimize the computation. The problem is simplified as well as the space search by considering that the optimum paths between two nodes at any departure time are almost always topologically the same. Nevertheless, this constraint cannot be considered in realistic multimodal networks, because perturbations modify the conditions of the networks at different times of the day. Bielli et al. [2] consider the problem of multimodal transport between towns. The problem is geographically divided and abstracted with hierarchical graph: towns are represented by subgraphs, which are interconnected by direct links. The shortest paths inside subgraphs are previously precomputed, reducing the search space. The main goal of this work is the proposed abstraction of the problem, which can be partially parallelized. The proposed approach is tested in networks composed on many smallsized subgraphs.
Table 5 presents a comparison between these works and transfer graph approach on the same instances of the problem. Since different machines are used in the experiments (Pentium II, UltraSPARC III, and T7300), computation time has been normalized by two benchmarks [15, 16]. According to these benchmarks, T7300 and UltraSPARC III are 13.3 and 1.32 times faster than Pentium II, respectively. These ratios are applied to computation time and included in Table 5. After these corrections, our results are still faster and advantageous, even if they were performed in an old Pentium II. The overall results show that transfer graph approach performs better with respect to the normalized computation time. The possible reason of this performance is due to precalculations, and the decomposition of the SPP in transfer graph simplified the problem.

6. Conclusion and Future Work
The main contribution of this research is an approach for the multimodal transport problem which deals with distributed nature of real transport information providers.
The transfer graph can abstract scenarios with distributed transport information sources, since it separates and keeps all transport modes in different unimodal networks. Besides, each unimodal network can be easily updated when unexpected perturbations are detected. This characteristic is very important because the RWS and TMS provide in realtime information about weather and traffic condition, and if an incident or delay is detected, only those unimodal networks related to the warning area must be recalculated, instead of the total network.
The decomposition of the SPP in transfer graphs provides an adequate solution for systems which receive multiple and frequent requests. Nevertheless, precalculations limited the scalability of the transfer graph abstraction in terms of space memory.
The transfer graph approach has been tested with several instances of the SPP, calculating the solution in less than one second in average. This approach provides better performance in terms of computation time than other related works.
The transfer graph approach has been developed as a multimodal transport service for the Carlink Platform and validated in real conditions using several transport means. This service will be deployed and made active and will provide an important solution for thousand of daily travellers in Luxembourg.
Our future work is to take advantage of the inherent parallel and distributed structure exhibited by transfer graph. For this purpose, relevant paths will be computed by using multiobjective Evolutionary Algorithm within each component separately and concurrently and will be also experimented and tested in large real multimodal network spanning many countries and involving several transportation service providers.
Acknowledgments
This work is partially supported by several institutions: the European Union Eureka cluster program Celtic Label , the Luxembourg Government, and the Spanish Regional Government of Andalusia under contract P07TIC03044 (DIRICOM). The authors wish to thank all partners in the Carlink project.
References
 K. G. Zografos and K. N. Androutsopoulos, “Algorithms for itinerary planning in multimodal transportation networks,” IEEE Transactions on Intelligent Transportation System, vol. 9, pp. 175–184, 2008. View at: Google Scholar
 M. Bielli, A. Boulmakoul, and H. Mouncif, “Object modeling and path computation for multimodal travel systems,” European Journal of Operational Research, vol. 175, pp. 1705–1730, 2006. View at: Google Scholar
 T.S. Chang, “Best routes selection in international intermodal networks,” Computers and Operations Research Archive, vol. 35, pp. 2877–2891, 2008. View at: Google Scholar
 Q. Li and C. E. Kurt, “GISbased itinerary planning system for multimodal and fixedroute transit network,” in Proceedings of the MIDContinent Transportation Symposium, pp. 47–50, Iowa State University, Ames, Iowa, USA, May 2000. View at: Google Scholar
 A. Ziliaskopoulos and W. Wardell, “An intermodal optimum path algorithm for multimodal networks with dynamic arc travel times and switching delays,” European Journal of Operational Research, vol. 125, no. 3, pp. 486–502, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. M. Foo, H. W. Leong, Y. Lao, and H. C. Lau, “A multicriteria, multimodal passenger route advisory system,” in Proceedings of the IESCTR, Singapore, July 1999. View at: Google Scholar
 A. Lozano and G. Storchi, “Shortest viable hyperpath in multimodal networks,” Transportation Research Part B, vol. 36, no. 10, pp. 853–874, 2002. View at: Publisher Site  Google Scholar
 D. Chiu, O. Lee, H.F. Leung, E. Au, and M. Wong, “A multimodal agent based mobile route advisory system for public transport network,” in Proceedings of the Annual Hawaii International Conference on System Sciences (HICSS '05), p. 92, Big Island, Hawaii, USA, January 2005. View at: Google Scholar
 M. Fragouli and A. Delis, “Easytransport: an effective navigation and transportation guide for wide geographic areas,” in Proceedings of the 14th IEEE International Conference on Tools with Artificial Intelligence (ICTAI '02), pp. 107–113, Washington, DC, USA, November 2002. View at: Google Scholar
 H. Ayed, D. Khadraoui, Z. Habbas, P. Bouvry, and J. F. Merche, “Transfer graph approach for multimodal transport problems,” in Proceedings of the 2nd International Conference on Modelling, Computation and Optimization in Information Systems and Management Sciences (MCO '08), pp. 538–547, Springer, September 2008. View at: Google Scholar
 “CARLINK: Wireless Platform for Linking Cars,” 2006–2008, http://carlink.lcc.uma.es/. View at: Google Scholar
 S. Pallottino and M. Scutella, “Shortest path algorithms in transportation models: classical and innovative aspects,” in Proceedings of the Equilibrium and Advanced Transportation Modelling Colloquium, Kluwer Academic Publishers, 1997. View at: Google Scholar
 E. W. Dijkstra, “A note on two problems in connexion with graphs,” Numerische Mathematik, vol. 1, pp. 269–271, 1959. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Kangas, M. Hippi, J. Ruotsalainen et al., “The FMI road weather model,” HIRLAM Newsletter, vol. 51, pp. 117–123, October 2006. View at: Google Scholar
 “Standard performance evaluation corporation (SPEC),” http://www.spec.org/. View at: Google Scholar
 “Geekbench benchmark,” http://www.primatelabs.ca/geekbench/. View at: Google Scholar
Copyright
Copyright © 2009 Carlos GalvezFernandez et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.