Traffic Assignment Procedure#

Along with a network data model, traffic assignment is the most technically challenging portion to develop in a modeling platform, especially if you want it to be fast. In AequilibraE, we aim to make it as fast as possible, without making it overly complex to use, develop and maintain, although we know that complex is subjective.

Running traffic assignment in AequilibraE consists in creating the traffic classes that are going to be assigned, add them to a traffic assignment object, set the traffic assignment parameters, and run the assignment.

TrafficClass#

The TrafficClass object holds all the information pertaining to a specific traffic class to be assigned. There are three pieces of information that are required in the instantiation of this class:

  • name: name of the class. It has to be unique among all classes used in a multi-class traffic assignment

  • graph: it is the Graph object corresponding to that particular traffic class/mode

  • matrix: it is the AequilibraE matrix with the demand for that traffic class, which can have an arbitrary number of user-classes setup as different layers (cores) of the matrix object.

>>> from aequilibrae.paths import TrafficClass

>>> project = create_example(project_path)
>>> project.network.build_graphs()

# We get the graphs for cars and trucks
>>> graph_car = project.network.graphs['c']
>>> graph_truck = project.network.graphs['T']

# And also get the matrices for cars and trucks
>>> matrix_car = project.matrices.get_matrix("demand_mc")
>>> matrix_car.computational_view("car")

>>> matrix_truck = project.matrices.get_matrix("demand_mc")
>>> matrix_truck.computational_view("trucks")

# We create the Traffic Classes
>>> tc_car = TrafficClass("car", graph_car, matrix_car)
>>> tc_truck = TrafficClass("truck", graph_truck, matrix_truck)

It is also possible to modify the default values for the following parameters of a traffic classe by using a method call:

  • Passenger-car equivalent (PCE) is the standard way of modeling multi-class traffic assignment equilibrium in a consistent manner (see [3] for the technical detail), and its value is set to 1.0 by default.

>>> tc_truck.set_pce(2.5)

  • Fixed costs: in case there are fixed costs associated with the traversal of links in the network, the user can provide the name of the field in the graph that contains that network.

>>> tc_truck.set_fixed_cost("distance")

  • Value-of-Time (VoT) is the mechanism to bring time and monetary costs into a consistent basis within a generalized cost function. In the event that fixed cost is measured in the same unit as free-flow travel time, then vot must be set to 1.0.

>>> tc_truck.set_vot(0.35)

TrafficAssignment#

>>> from aequilibrae.paths import TrafficAssignment

>>> assig = TrafficAssignment()

AequilibraE’s traffic assignment is organized within an object with the same name which contains a series of member variables that should be populated by the user, providing thus a complete specification of the assignment procedure.

  • classes: list of completely specified traffic classes

# You can add one or more traffic classes to the assignment instance
>>> assig.add_class(tc_truck) 

>>> assig.set_classes([tc_car, tc_truck])

  • vdf: the volume-delay function (VDF) to be used, being one of BPR, BPR2, CONICAL, or INRETS

>>> assig.set_vdf('BPR')

  • vdf_parameters: the parameters to be used in the volume-delay function, other than volume, capacity and free-flow time. VDF parameters must be consistent across all graphs.

    Because AequilibraE supports different parameters for each link, its implementation is the most general possible while still preserving the desired properties for multi-class assignment, but the user needs to provide individual values for each link OR a single value for the entire network.

    Setting the VDF parameters should be done AFTER setting the VDF function of choice and adding traffic classes to the assignment, or it will fail.

# The VDF parameters can be either an existing field in the graph, passed as a parameter:
>>> assig.set_vdf_parameters({"alpha": "b", "beta": "power"}) 

# Or as a global value:
>>> assig.set_vdf_parameters({"alpha": 0.15, "beta": 4})

  • time_field: the field of the graph that corresponds to free-flow travel time. The procedure will collect this information from the graph associated with the first traffic class provided, but will check if all graphs have the same information on free-flow travel time

>>> assig.set_time_field("free_flow_time")

  • capacity_field: the field of the graph that corresponds to the link capacity. The procedure will collect this information from the graph associated with the first traffic class provided, but will check if all graphs have the same information on capacity

>>> assig.set_capacity_field("capacity")

  • algorithm: the assignment algorithm to be used, being one of all-or-nothing, bfw, cfw, fw, franke-wolfe, or msa.

>>> assig.set_algorithm("bfw")

Volume-delay function#

For now, the only VDF functions available in AequilibraE are

\[CongestedTime_{i} = FreeFlowTime_{i} * (1 + \alpha * (\frac{Volume_{i}}{Capacity_{i}})^\beta)\]

  • Spiess’ conical [2]

\[CongestedTime_{i} = FreeFlowTime_{i} * (2 + \sqrt[2][\alpha^2*(1- \frac{Volume_{i}}{Capacity_{i}})^2 + \beta^2] - \alpha *(1-\frac{Volume_{i}}{Capacity_{i}})-\beta)\]

  • and French INRETS (alpha < 1)

Before capacity

\[CongestedTime_{i} = FreeFlowTime_{i} * \frac{1.1- (\alpha *\frac{Volume_{i}}{Capacity_{i}})}{1.1-\frac{Volume_{i}}{Capacity_{i}}}\]

and after capacity

\[CongestedTime_{i} = FreeFlowTime_{i} * \frac{1.1- \alpha}{0.1} * (\frac{Volume_{i}}{Capacity_{i}})^2\]

More functions will be added as needed/requested/possible.

Setting Preloads#

We can also optionally include a preload vector for constant flows which are not being otherwise modelled. For example, this can be used to account for scheduled public transport vehicles, adding an equivalent load to each link along the route accordingly. AequilibraE supports various conditions for which PT trips to include in the preload, and allows the user to specify the PCE for each type of vehicle in the public transport network.

To create a preload for public transport vehicles operating between 8 AM to 10 AM, do the following:

>>> from aequilibrae.transit import Transit

# Times are specified in seconds from midnight
>>> transit = Transit(project)
>>> preload = transit.build_pt_preload(start=8*3600, end=10*3600)

# Add the preload to the assignment
>>> assig.add_preload(preload, 'PT_vehicles')

Executing an Assignment#

Finally, run traffic assignment!

>>> assig.execute()

References#