Optimization of signalized network configurations using the Lane-based method

Vehicle movements at signalized intersections should follow the guidance of lane marking arrows. Turns are permitted or banned depending on the existence of lane marking arrows establishing network link connectivity. Lane markings are the interface that joins consecutive upstream and downstream intersections. Traffic flows from origins to destinations across intersections should be governed by lane markings. In this study, conventionally fixed lane markings are relaxed as binary variables to be optimized by maximizing the green bandwidths. The proposed methodology is innovative in that it extends the lane-based design framework to incorporate green band maximization for enhancing traffic signal coordination. Path flows are controlled to satisfy flow conservations and to compile turning flows at intersections. With these turning flows as inputs, lane markings can be optimized together with the coordinated traffic signal settings. Path flows, path travel times, and path choices are evaluated through new linear constraints. For path travel times, cruise times along lanes and total delays at the ends of lanes are evaluated. The non-linear total delays are approximated by the proposed linearized delay function. The model coefficients are calibrated by network data as in a previous study. The problem is formulated as a binary-mixed-integer-linear-program and solved by standard branch-and-bound techniques using a CPLEX solver. To avoid non-linearity in the problem formulation, the bandwidth maximization approach is adopted instead of minimizing the total network delay in the design of the signalized network. A network with four intersections is provided to demonstrate how optimized lane markings can create efficient network link configurations. The numerical results are promising compared with those obtained in previous studies and show that the overall network performance can be improved.


Introduction
In road networks, links with shorter travel times are preferred by users over travel between origins and destinations (OD). Route choices depend on total travel times along links and delays at intersections. In signalized networks, signal settings can affect route choices [1][2][3][4][5]. Examples of routing optimization problems for moving hazardous materials are based on multi-objective functions by probabilistic search methods such as those of neural networks and genetic algorithms [6][7][8][9]. Coordinated signal controls in ring road systems can be improved by proper speed guidance [10]. Network flow patterns depend on signal settings that are updated until consistent patterns are achieved [11][12][13]. A comprehensive review of this approach can be found in the literature [14]. A bi-level programing formulation integrated a lower-level user equilibrium assignment model and an upper-level control optimization model. Signal operations were optimized in the upper level and users' equilibrium was achieved in the lower level [15][16][17][18][19][20].
The coordination of linked traffic signals is important when considering signalized networks. The TRANSYT traffic model simulates network traffic considering platoon dispersions and signal coordination [21]. A path-based assignment routine was developed to enhance the TRANSYT model so that traffic flows could become responsive to signal settings [22]. Groupbased signal optimization was conducted for this enhanced TRANSYT platform [23]. A sensitivity analysis, based on the derivatives of the performance index with respect to signal settings, was conducted to optimize the network signal settings [24]. By displaying different signal timings, shorter or longer paths inside networks are generated as OD route choices for users. Linked connections also determine the availability of those paths. Such network configurations have normally been given as users' fixed inputs. However, lane marking arrows show users which turns are permitted when approaching intersection levels are variables and should be designed according to the demand flow patterns [25][26][27][28]. Practical demonstrations of the lanebased methods for signalized intersections can be found in the literature [29,30]. More recent studies applying the concept of lane-based design for signalized intersections are noted below. The lane assignments for multiple left-turn and right-turn lanes were determined by optimizing a multi-objective function with microscopic traffic simulations. Transportation efficiency, energy consumption, and road safety were remarkably improved [31]. The lane-based concept, varying the approach lanes at signalized intersections, was used in Dalian. The designs of timevarying traffic signals and intersection settings at different times of the day were produced by minimizing the total delay. Four scenarios with field data verifications demonstrated promising results [32]. Lane allocations and signal timings were optimized to design symmetric intersection layouts. A binary-mixed-integer-linear-program was formulated and solved by the standard branch-and-bound method. Instead of the usual four phases of traffic signal settings at intersections, three were found to be sufficient for separating all conflicting movements [33]. A new stage-based modeling approach was developed to prohibit left-turn traffic in urban road networks by minimizing total travel time [34]. The underlying concept is similar to the one proposed in this study. The proposed lane-based optimization framework is more computationally attractive, generating lane markings that configure an urban road network in which the proper locations to permit or prohibit turns may be optimized. A modeling interface was established between the connectivity of network links and lane marking designs at intersections [35]. Maximizing the common flow multiplier to scale the whole OD (input) demand flow matrix, as done in a previous study [36], can maximize the capacities of all individual signalized intersections within a network and can generate lane markings for configuring the network connections [37]. However, equilibrium flow conditions and proper signal coordination were not considered.
The MAXBAND method was developed to maximize the bandwidths across adjacent intersections and ensure smooth arterial progression across intersections [38]. This concept was extended to address multi-arterial systems [39]. Systematic traffic-dependent criteria were developed for optimizing a weighted sum of bandwidths to form a new MULTIBAND model [40,41]. MULTIBAND was applied to design asymmetrical traffic signal coordination for arterials [42] and model arterial networks [43]. In the MULTIBAND model, individual bandwidths are weighted based on ratios of flow volumes to saturation flow rates (i.e. flow factors). In their formulation, flow volumes and saturation flow rates are known as given parameters. In the proposed lane-based optimization framework, both flow factors and bandwidths are variables and their products in MULTIBAND become a nonlinear term. Vehicles from the upstream intersections enter downstream links and take the link travel times to the downstream intersections. Average travel times depend on both the physical separation between upstream and downstream intersections and average vehicle speeds. The slopes of the vehicle trajectories in the time-space diagram represent average vehicle speeds [38]. Ideally, vehicles can move without waiting if proper green signal phases and offsets are displayed. Offsets are the time differences between two green start times at the upstream and downstream intersections. The signalized network is well coordinated if the offsets are properly designed.
To optimize the lane marking patterns at individual intersections, 13 sets of linear constraints are required, taking the turning flow patterns as model inputs for defining a feasible solution region for the traffic signal settings [25,27,37]. The constraints are set to (i) ensure that the number of lane marking arrows for turns from upstream do not exceed their respective numbers of downstream exit lanes; (ii) satisfy the requirements for flow conservation, equalizing the given turning flows and the sum of assigned lane flows; (iii) establish lane markings to provide consistency in the assigned turning flows; (iv) avoid (internal) conflicts of lane markings across adjacent approach lanes from the same approaches; (v) establish single lane markings (for a single turn) or shared lane markings (for multiple turns) on approach lanes, thus ensuring the use of all approach lanes; (vi) restrict the optimized cycle time to be within the users' specified allowable range; (vii) ensure that the optimized green duration times exceed the users' specified minimum duration of green times; (viii) restrain the starts of green times within the traffic signal cycle; (ix) provide minimum clearance times (or intergreen times) to separate the right-of-way of conflicting traffic movements; (x) introduce successor functions to regulate the signal display orders for conflicting pairs of traffic movements; (xi) synchronize the traffic signal settings for turns and lanes according to the optimized lane marking patterns; (xii) equalize the flow factors on adjacent approach lanes with identical lane marking designs on both approach lanes; and (xiii) restrict the maximum degree of saturation to below the users' specified figure for all approach lanes. These 13 sets of governing constraints are standard formulation for the lane-based design for individual signalized intersections and have been applied for reserve capacity maximization [25], total delay minimization [26], and intersection layouts with short lanes [30]. Other applications can also be found for U-turn traffic control [44], designing the layout of a displaced left-turn intersection [45], handling fluctuating traffic demand [46], urban arterials with reversible lane usages [47], for exclusive bus lanes [48], and designing pre-signals with sorting areas before entering the common areas of intersections [49].
As reviewed above, in most research methods relating to signalized network designs the network configurations are fixed users' inputs, including the connections of network links and their link saturation flow rates, indicating that their discharge abilities are excluded from the optimization framework. The given OD flows might be forced to use the available but congested paths provided by inefficient link connections leading to sub-optimal flow distributions. In a recent study [37], lane markings were optimized by maximizing the common flow multiplier to the given OD matrix, thus generating network configurations. However, equilibrium flow conditions and the coordination of traffic signal settings were not modeled, leading to unrealistic results. In this study, users' movement paths are free to be generated depending on the optimization results of lane markings and signal timings that serve all OD flows from the given matrix. This study attempts to relax lane markings as design variables and optimize the network configurations such that the overall link connectivity of a signalized network is enhanced. All network links or lanes can be used effectively to let through all OD demand flows. Offsets and other signal settings are optimized in a unified framework by maximizing the green bandwidths along the selected major routes. Path travel times are evaluated containing lane cruise times and delays at the ends of lanes to present users' equilibrium flow conditions. A network consisting of four signalized intersections is given as a case study to show the effectiveness of the proposed algorithms. Five scientific measures are evaluated for their assessment of the performance of optimized network settings: (I) the lengths of green bands, (II) the numbers and patterns of lane markings to enable feasible paths for connecting OD pairs, (III) the effectiveness of signal timings to coordinate upstream and downstream intersections, (IV) the equilibrium flow conditions for using different paths to connect the OD pairs, and (V) total network travel times X all paths path flow � path travel time The mathematical problem is formulated as a binary-mixed-integer-linear-program (BMILP) and is solved by a standard branch-and-bound method through the CPLEX solver. Network settings optimized by the proposed algorithms are then compared with those obtained in a previous study to show their enhancement.

Linear governing constraint sets
Flow conservation to equalize path flows and given OD flows. Eq (1) is an equality constraint set for equalizing the sum of available path flows and the given OD demand flows, thus satisfying the flow conservation requirements: Turning flows at intersections. OD demand flows are usually given by users as model inputs that are to be assigned onto different paths. Path flows may pass through different intersections in the network. Turning flows from arm i to arm j at intersection n can be deduced from the path flows. Fig 1 shows  Required lane markings for a path to serve OD pairs. In the present formulation, lane markings are optimized to connect network links and ensure that the required paths in Eq (7) serve a network's OD flows. i to arm j at intersection n is required. Eq (9) is required in the design framework for optimizing lane markings.
Classifying the sequences and patterns of green duration times by transforming green signal phases for bandwidth evaluations. The shapes typical of the green bands between two consecutive intersections are illustrated in the time-space diagram given in Fig 2. To evaluate the green bandwidth in the proposed formulation, a linear transformation is performed to shift green phases from the upstream intersections horizontally along the time axis with the duration of average travel times from the upstream to downstream intersections. After the transformation, linear constraint sets can be developed based on simple rules to compare the start of green and end of green times for the upstream and downstream intersections. Auxiliary binary variables are introduced to monitor different outcomes of the comparisons for measuring the actual green bandwidths across pairs of consecutive upstream and downstream intersections.
As Fig 2 shows, there are two green signal phases for upstream intersection n and downstream intersection m, respectively. To evaluate the green bandwidth between the two green signal phases (i.e. thickened solid lines), a linear transformation is applied to shift the original green signal phase at upstream intersection n horizontally along the time axis with the amount of t n;m ¼ S n;m =V n;m , which is the average travel time for vehicles leaving from upstream intersection n to reach downstream intersection m. S n;m is the distance separation between upstream n and downstream m and V n;m is the average speed of vehicles moving from upstream n to downstream m. After performing the transformation, the revised green start time is Y n;i;k þ t n;m . Similarly, the end of green time also moves to Y n;i;k þ F n;i;k þ t n;m . Two vertical dashed lines [(i) and (ii)] may be drawn: (i) from the end of green at upstream intersection n and (ii) from the start of green at downstream intersection m. Furthermore, bandwidth W n;m can be measured, representing the overlapping length of the two green signal phases [i.e. between the two vertical dashed lines (i) and (ii)] after the transformation. The practical transformation process is operated as follows.
The transformation process enables the evaluations of bandwidth found in the proposed programming approach. In the optimization process, traffic signal settings and green bandwidths can be varied and optimized with different and overlapping patterns used to generate the actual bandwidth W n,m . Linear constraint sets are developed to differentiate all of these cases and evaluate the corresponding bandwidth W n,m . In modeling, four conditions (I-IV) in the signal settings are considered. By comparing the transformed start of green time of upstream intersection n, which is Y n;i;k þ t n;m , and the start of green time of downstream intersection m, Y m;i 0 ;k 0 þ x n;m (where x n,m is an integer variable = 0, 1, 2, . . .), we may subtract the two green start times in Eq (10) to establish Condition I. ML 1;n;i;k;m;i 0 ;k 0 � Y m;i 0 ;k 0 þ x n;m À Y n;i;k À t n;m � MðL 1;n;i;k;m;i 0 ;k 0 À 1Þ; 8ðm; i 0 Þ 2 Dðn; i; jÞ; n 2 f1; 2; . . . ; Ng; i 6 ¼ j 2 f1; 2; . . . ; I n g; k 2 f1; 2; . . . ; L n;i g; k 0 2 f1; 2; . . . ; L m;i 0 g ð10Þ where L 1;n;i;k;m;i 0 ;k 0 is a binary variable (= 1 if the transformed start of green time for the vehicle movement of the upstream intersection ðn; i; kÞ is earlier than the start of green time for the vehicle movement of the downstream intersection ðm; i 0 ; k 0 Þ, or = 0 if not). L n;i is the total number of approach lanes on arm i at intersection n. L m;i 0 is the total number of approach lanes on arm i 0 at intersection m. A user's defined mathematical function Dðn; i; jÞ is to identify arm i' at downstream intersection m, (m, i'), taking the traffic from arm i to arm j at upstream intersection n.
To establish Condition II, a similar comparison as in Eq (10) is required for the two ends of the green times of upstream intersection n and downstream intersection m. In Eq (11), L 2;n;i;k;m;i 0 ;k 0 is another binary variable to state the comparison result (= 1 if the transformed end of green time for the vehicle movement of the upstream intersection ðn; i; kÞ is earlier than the end of green time for the vehicle movement of the downstream intersection ðm; i 0 ; k 0 Þ, or = 0 if not). ML 2;n;i;k;m;i 0 ;k 0 � ðY m;i 0 ;k 0 þ F m;i 0 ;k 0 Þ þ x n;m À ðY n;i;k þ F n;i;k Þ À t n;m � MðL 2;n;i;k;m;i 0 ;k 0 À 1Þ; 8ðm; i 0 Þ 2 Dðn; i; jÞ; n 2 f1; 2; . . . ; Ng; i 6 ¼ j 2 f1; 2; . . . ; I n g; k 2 f1; 2; . . . ; L n;i g; k 0 2 f1; 2; . . . ; L m;i 0 g ð11Þ To establish Condition III, the comparison continues for the transformed end of green time for the vehicle movement of the upstream intersection ðn; i; kÞ and the start of green time for the vehicle movement of the downstream intersection ðm; i 0 ; k 0 Þ in Eq (12). Again, binary variable L 3;n;i;k;m;i 0 ;k 0 states whether the comparison result (= 1 if the start of green time for vehicle movement of the downstream intersection ðm; i 0 ; k 0 Þ is earlier than the transformed end of green time for the vehicle movement of the upstream intersection ðn; i; kÞ, or = 0 if not).
Maximum length of green bandwidth. Apart from mathematical constraint sets developed in previous sections, which relate the adjacent upstream and downstream intersections, it must be noted that the actual (green) bandwidth should always be shorter than or equal to the green duration times operating at individual intersections for the traffic moving along the selected major routes and OD pairs. Eqs (20) and (21) are given to restrict the maximum (green) bandwidth in the solution process. 0 � W n;m � F n;i;k ; 8ðm; i 0 Þ 2 Dðn; i; jÞ; n 2 f1; 2; . . . ; Ng; i 6 ¼ j 2 f1; 2; . . . ; I n g; k 2 f1; 2; . . . ; L n;i g ð20Þ 0 � W n;m � F m;i 0 ;k 0 ; 8ðm; i 0 Þ 2 Dðn; i; jÞ; n 2 f1; 2; . . . ; Ng; i 6 ¼ j 2 f1; 2; . . . ; I n g; k 0 2 f1; 2; . . . ; L m;i 0 g ð21Þ Bandwidths for coordinating upstream and downstream intersections. To put forward an effective coordination between upstream and downstream traffic signal settings, maximizing the green bandwidth is set as the objective in the present formulation. Before bandwidth can be optimized in the objective function, we must ensure that the bandwidth exists in which all transformed upstream green signal phases can overlap the downstream green signal phases. There are two possibilities, governed by Eqs (22) and (23).
If the transformed start of green time for the vehicle movement of the upstream intersection ðn; i; kÞ is earlier than the start of green time for the vehicle movement of the downstream intersection ðm; i 0 ; k 0 Þ, then to ensure that the two signal phases can overlap to produce effective green bandwidth, as in Cases 1 and 3, we need to force the transformed end of green time at the upstream intersection ðn; i; kÞ to be later than the start of green time at the downstream intersection ðm; i 0 ; k 0 Þ. The required constraint set is given in Eq (22). 1 À L 1;n;i;k;m;i 0 ;k 0 � L 3;n;i;k;m;i 0 ;k 0 À 1 � L 1;n;i;k;m;i 0 ;k 0 À 1; 8ðm; i 0 Þ 2 Dðn; i; jÞ; n 2 f1; 2; . . . ; Ng; i 6 ¼ j 2 f1; 2; . . . ; I n g; k 2 f1; 2; . . . ; L n;i g; k 0 2 f1; 2; . . . ; L m;i 0 g ð22Þ If the transformed start of green time at the upstream intersection ðn; i; kÞ is later than the start of green time at the downstream intersection ðm; i 0 ; k 0 Þ, as in Cases 2 and 5, then the end of green time at the downstream intersection ðm; i 0 ; k 0 Þ must be later than the transformed start of green at the upstream intersection ðn; i; kÞ. The related constraint set is given in Eq (23). L 1;n;i;k;m;i 0 ;k 0 � L 4;n:i;k;m;i 0 ;k 0 À 1 � À L 1;n;i;k;m;i 0 ;k 0 ; 8ðm; i 0 Þ 2 Dðn; i; jÞ; n 2 f1; 2; . . . ; Ng; i 6 ¼ j 2 f1; 2; . . . ; I n g; k 2 f1; 2; . . . ; L n;i g; k 0 2 f1; 2; . . . ; L m;i 0 g ð23Þ In this study, maximum bandwidths are designed to be the maximum length that can be optimized "in common" from the overlapping bandwidths shown in Total intersection delay by approximated delay function. A 3-term Webster's delay formula was developed (but could be simplified and approximated by using the first two terms of the original formula) for the uniform and random delays at signalized intersections [50]. Webster's delay formula grows infinitely and unrealistically large when the degree of saturation approaches 100%. Then, a proposal emerged for developing a coordinate transformation technique to predict the intersection delays that are subject to time-varying traffic demand and capacity for both under-and over-saturated conditions [51]. According to the TRANSYT traffic model [21,22,35,52], link travel times include two components: link cruise time and end of link delay. The cruise time on a link has been considered to be the free flow travel time along the link. The delay at the end of link consists of uniform delay due to alternations of green and red traffic signals and random and oversaturation delay due to random traffic arrivals and degree of saturation > 100%, respectively [53]. Uniform delay can be calculated directly using the first term of Webster's equation that depends on cycle length, the proportion of a signal cycle, which is effectively green, and the degree of saturation. Random and over-saturation delay can be evaluated based on a sheared delay formula [24], depending on link capacity, initial queue length, degree of saturation, analysis period, and the constant in the Pollaczek-Khintchine (P-K) formula [54]. The P-K constant in the sheared delay formula was calibrated for use in a coordinated network in Hong Kong [55].
Obviously, the shear delay formula calculating delays at the ends of links are non-linear functions. Adding those uniform delay and random and oversaturation delay (formula) directly to the proposed formulation may turn the entire design problem into a binary-mixedinteger-non-linear-program, which is very difficult to solve [26,28]. The linear approximation technique and simple functional form were successfully applied to predict total delays in a signalized system [56]. In this study, a similar concept and linearization technique is applied to evaluate the delay time at the end of a lane that is related to its lane flow (including turning and straight-ahead movements), proportion of a signal cycle which is effectively green, and degree of saturation. For slowly varying non-linear delay function, it is possible to decouple the original product terms into individual linear terms for approximations. In the present formulation, total lane flow, ratio of lane flow to saturation flow, effective green duration time, and cycle length are extracted from the lane-based model to be the linearized function variables for predicting delays at ends of lanes. Eq (25) is the approximated function.
8n 2 f1; 2; . . . ; Ng; i 2 f1; 2; . . . ; I n g; k 2 f1; 2; . . . ; L n;i g ð25Þ Optimization of signalized network configurations using the Lane-based method traveling along the approach and the delay at the downstream intersection. In practice, a simplified Bureau of Public Roads function depending on flow-to-capacity ratio is commonly used. Obviously, such link performance function is non-linear in nature due to a power term applied to the flow-to-capacity ratio resulting with asymptotic link travel time values when link flows approach their capacity values (i.e. degree of saturation is near 100.0%). In the present linear formulation, directly applying the link performance function to calculate path travel times may turn the entire formulation into a non-linear problem that is difficult to solve.
In the proposed formulation, the lane travel time, g n;i;k , in Eq (26) is the sum of (1) vehicle traversing time on the used lane between upstream intersection n and downstream intersection m, t n;m ¼ S n;m =V n;m (i.e. the cruise time) and (2) total delay at the end of the lane, φ n;i;k . g n;i;k ¼ t n;m þ � n;i;k ; 8ðm; i 0 Þ 2 Dðn; i; jÞ; n 2 f1; 2; . . .; Ng; i 6 ¼ j 2 f1; 2; . . .; Ng; k 2 f1; 2; . . .; L n;i g ð26Þ In this section, linear constraint sets are developed to evaluate the travel times along lanes using piecewise linear functions. A two-piece linearized function is proposed as a replacement to maintain the linearity in the formulation, as shown in Fig 6. To match the original shape of the non-linear link performance function in Fig 5, we have to ensure that a sharp increase in the travel time along a lane if the lane flow exceeds the lane capacity given in Fig 6. The threshold position is near the flow-to-capacity ratio = 100% (i.e. demand flow = capacity) as given by the original link performance function. In traffic signal controls, the degree of saturation is known as the flow-to-capacity ratio, which is equal to the lane flow factor multiplying cycle length and dividing by the effective green duration time. The corresponding threshold position appears when the flow factor equals the proportion of a signal cycle that is effectively green. In Fig 6, two linear lines, based on two different increasing slopes, simulate the original curve shape of the link performance function for evaluating lane travel time g n;i;k . In Eq (27), an auxiliary binary variable P n;i;k for lane k from arm i at intersection n is forced to be zero when the flow factor y n,i,k is smaller than or equal to the proportion of a signal cycle which is effectively green, F n;i;k for unsaturation conditions. Or P n;i;k ¼ 1 when y n;i;k > F n;i;k for over-saturation conditions. MP n;i;k � y n;i;k À F n;i;k > MðP n;i;k À 1Þ; 8n 2 f1; 2; . . . ; Ng; i 2 f1; 2; . . . ; I n g; k 2 f1; 2; . . . ; L n;i g ð27Þ Based on the binary results of P n;i;k , Eqs (28) and (29) are developed to control the numerical coefficients used to weight the lane flow factor in the proposed linear delay function in Eq (25) to evaluate the total delays at the ends of lanes. Whenever P n;i;k ¼ 0, Eq (28) forces the numerical coefficient Z y ¼ ε, which is the first segment of the 2-piece linearized function with a gentle increasing slope ε. Or P n;i;j ¼ 1, Eq (29) forces Z y ¼ ε 0 which is the second segment of the 2-piece linearized function with a sharper increasing slope ε 0 . MP n;i;k � Z y À ε � À MP n;i;k ; 8n 2 f1; . . . ; Ng; i 2 f1; . . . ; I n g; k 2 f1; . . . ; L n;i g ð28Þ Mð1 À P n;i;k Þ � Z y À ε 0 � À Mð1 À P n;i;k Þ; 8n 2 f1; 2; . . . ; Ng; i 2 f1; 2; . . . ; I n g; k 2 f1; 2; . . . ; L n;i g ð29Þ where ε < ε 0 numerically. Ensuring identical path travel times for used paths connecting origins and destinations. In a network, users tend to use the shortest path(s) to travel between OD pairs. In the present formulation, paths are composed by linking up series of intersections on different approach lanes based on lane marking patterns. To assign path flows onto different available paths, path travel times, and associated total delays at intersections should be considered. In Eq (30) (2). To ensure that all longer paths are unfavorable in route choices, Eq (31) can be used to guarantee that path travel times are longer when paths are not chosen by users. Path flows along those longer paths must be zero. With Eq (31) Optimization of signalized network configurations using the Lane-based method in the proposed formulation, the resultant traffic flow pattern satisfies the Wardropian userequilibrium principle in which no user in the network can unilaterally change routes to reduce travel times [58].

Objective function for maximizing green bandwidths
Green bandwidth maximization enables the coordination of traffic signal settings in a network system. With such an approach, it is possible to determine the best offsets of green start times for all involved signal phases along a selected major path. In this study, bandwidth across all pairs of consecutive intersections along selected major paths can be evaluated linearly. Other linear constraint sets for configuring network link connections, described above in an effort to outline a solution region for linear optimization, are setting up individual signal-controlled intersections with lane marking designs and safe traffic signal timings [20,22], evaluating total delays at ends of lanes, estimating lane travel times, and ensuring identical travel times for used paths. In the present formulation, lane markings are the key binary type variables to be optimized for the design of network configuration. To take full advantage of the design features of current lane-based methods, we need to maintain a linear platform in the optimization framework. Maximizing the green bandwidth is thus chosen as the objective for optimization in the proposed study. The problem can then be formulated as a BMILP problem encapsulating the lane markings, link connections for defining network configurations, traffic signal settings including offsets for coordinating signal timings across adjacent intersections, and green bandwidths for a signalized network design problem. A new approach is introduced to enable the design of network configurations, relaxing lane markings as design variables to support network link connections and link up adjacent intersections. As for the objective function, the overall common bandwidth is applied for establishing the linear objective function in Eq (32) for optimization. By maximizing the green bandwidths, it is possible to obtain optimized signal timings, lane marking patterns, lane flows, and path flows. To solve the proposed optimization problem, users aiming to receive green bands should provide their origin nodes, destination nodes, and selected major paths with their preferences. The selected major paths become more attractive for users when connecting their OD pairs. In a study network, OD demand flows are given as users' input data. The resultant lane markings on approach lanes at various intersections should be compiled to configure the network's link connections. Assigned lane flows that are compatible with the optimized set of lane markings can be optimized to match their respective path flows. In addition, all path flows can be combined to equalize the given OD demand flow patterns, thus satisfying the requirements for flow conservation. Optimized signal timings are effective and efficient for practical operations. Better linkages between pairs of upstream and downstream intersections can be established by providing proper offsets. Mathematically, the bandwidth maximization problem can be formulated as a BMILP: Optimization of signalized network configurations using the Lane-based method

Optimization results from a case study network with four intersections
In this section, a four-intersection signal-controlled network is modeled with N = 4. Lane markings for all intersections are to be optimized to configure network links and establish link connectivity. Fig 7 presents   Optimization of signalized network configurations using the Lane-based method In Fig 7, for some OD pairs, such as from origin o = 1 to destination d = 2, there are shorter paths that only pass through two intersections to reach their destinations, while the longer paths need to pass through all four intersections in the network. The turning radius is set to be 12.0m for all left-and right-turn movements from both nearside and non-nearside lanes. The arbitrary large positive number M takes on 10,000. Clearance times and minimum green duration times are set to be 6.0 and 5.0 seconds, respectively. Effective green duration times are 1.0 seconds longer than the actual (or displayed) green duration times i.e. e ¼ 1:0 second. To optimize two green bandwidths (for inbound and outbound traffic in opposite directions) in the present numerical example, traffic flow from origin o = 2 to destination d = 4 through intersections n = 2, n = 3, and n = 4 is defined as the first selected major path, and another selected major path is from origin o = 4 to destination d = 2 through the intersections n = 4, n = 3, and n = 2.  (32). The distances of all lanes connecting upstream and downstream intersections are assumed to be the same with t n;m ¼ S n;m =V n;m ¼ 26 seconds.
Based on the sheared delay function, which was calibrated and validated through the NET-SIM model using the Hong Kong coordinated network in Sha Tin and their traffic data as inputs, to predict the random-and-oversaturation delays [48] and the Webster's delay function for calculating the uniform delays (adding them together forming the total delays), the total delays at the ends of links are evaluated using the numerical data results from a signalized network in a previous study [30]. These calculated total delays are set as the target values to be approximated by computing the linearized delay formula in Eq (25). By minimizing the least square error between the two sets of total delays through genetic algorithms, the numerical coefficients of the delay formula in Eq (25) are determined with the following results, � n;i;k ¼ 59:25000 � z À 99:00000 � ðF n;i;k þ ezÞ þ 0:39583 � y n;i;k À 0:00004 � X j q n;i;j;k . The predictor variables in Eq (25) are also extracted numerically from the present lane-based network model for predicting the total delays at the ends of lanes.
Based on all the above model inputs, the proposed lane-based model can be applied to optimize the green bandwidth adopted in objective function Eq (32). Fig 8 shows  the offsets of the green start times are also modeled to control upstream and downstream traffic. Their offsets are also given in Fig 9. Due to the equal importance from the given numerical weights for the green bands with respect to the inbound and outbound traffic directions, another green band along the opposite traffic direction has the same green bandwidth, W o¼4;d¼2;h o;d ¼1 =z = 32.9s as given in Fig 10. To serve the OD flows T o;d in the example network, different paths are generated and used.  Optimization of signalized network configurations using the Lane-based method flows could then be revised, taking into consideration turning radius and turning flow proportions. Then, flow factors can be deduced from dividing the total lane flows by the revised saturation flows. These flow factors can also be used to allocate available green duration times to various traffic lanes, and used as model parameters for computing the total delays at ends of lanes using the proposed approximated linear function. Green start and effective green duration times are given in Tables 2-5. The last column presents the degree of saturation for all approach traffic lanes, all of which are unsaturated. In one of the equality constraint sets in the formulation, total delays at ends of lanes comprising uniform and random-and-oversaturation delays are used to predict the total delays from origins to destinations across different lanes and intersections. The linearized delay formula in Eq (25) has been proposed to approximate the non-linear uniform and sheared delays. Based on the optimization results, the two delays have been computed and are given in Table 6  Optimization of signalized network configurations using the Lane-based method for the purposes of comparison and verification. It is found that the two calculated delays for all of the traffic lanes are consistent, with discrepancies found to be within ±2.5%. These findings, which result from using a linear approximation delay formula to replace the non-linear sheared delay function, thus avoiding non-linearity, seem to indicate a possible alternative approach.

Sensitivity analysis and result comparisons
To illustrate how the proposed algorithms enhance the performance of the optimized network settings, we first optimize the same network by taking the same OD flow matrix and other network settings as inputs, specifically by maximizing the common flow multiplier according to the method developed from a recent study [31] and compare their results. In addition, to conduct a sensitivity analysis, all shared lane markings in the optimized network configuration are removed for evaluation. Five performance indicators are retrieved from the optimization results for comparison: (I) the length of green bands along major paths, (II) the numbers and patterns of lane marking arrows on approach traffic lanes to enable feasible paths serving the given OD flows, (III) the effectiveness of signal timing for coordinating upstream and downstream traffic signals, (IV) existence of equilibrium flow conditions for using different paths to connect the same OD pairs, and (V) total network travel times obtained from the latest method used to establish network link connections by varying and optimizing lane markings at the intersection levels. The method maximizes the common flow multiplier of the given OD flow matrix. The scaled OD flows that enter the study network are subject to a simple requirement for flow conservation based on equalizing the given OD flows and the sum of path flows in the used paths. The method optimizes the traffic signal settings for all individual intersections without modeling traffic signal coordination, which is expected to be refined by the proposed study method. The proposed method enhances the five performance indicators mentioned above, allowing them to obtain more convincing and realistic results.
For performance indicator (I) the widths of green bands (across 3 intersections), the previous method give 15.4 s for connecting origin o = 2 to destination d = 4 through intersection n = 3 and zero in the opposite direction. It also records 2.1s for connecting origin o = 3 to destination d = 1 through intersection n = 2 and zero in its opposite direction, according to the green start and duration times given in Fig 13. The proposed algorithms are applied to two major paths defined by users as inputs connecting origin o = 2 to destination d = 4 through intersection n = 3 (inbound direction) and its opposite direction connecting origin o = 4 to destination d = 2 through intersection n = 3 (outbound direction). The sum of the OD flows for this OD pair is the largest (= 600 + 400 = 1,100 pcu/h). The proposed MAXBAND For performance indicator (III), the proposed algorithm's optimized signal timings are more effective in generating longer green bands along users' specified major paths, as found in performance indicator (I). Both sets of signal timing results, as given in Figs 8 and 13, are able to separate conflicting traffic movements at the intersection level. Hence, the proposed optimization algorithms lead to better coordination of optimized traffic signal settings. As for performance indicator (IV), users from the same OD pair, under equilibrium conditions, may travel using different paths and no user should be able to save travel time by switching to another path. In addition, the resultant path travel times for the same OD pairs should be identical. The proposed algorithms develop the necessary governing constraint sets to control path flows, satisfying this path travel time requirement. From origin o = 1 to destination d = 3 or 4 and o = 2 to d = 4, there are two available paths optimized by the previous method, as given in Fig 14 and Table 7. Their respective path travel times are not identical, violating the equilibrium requirement. This problem is completely resolved by the proposed design algorithms presented in Fig 12 and Table 8. The two path travel times are now equalized: the path from o = 1 to d = 3 is now equalized to 169.5s and the path from o = 2 to d = 4 is now equalized to 162.1s. For path 1 and path 2, path travel times from o = 1 to d = 4 are now 235.6 seconds and 81.7 seconds, respectively. The OD flow (= 400 pcu/h) uses path 2 with a shorter travel time (= 81.7s) and the flow along path 1 is zero due to a longer travel time (= 235.6s). For performance indicator (V), both methods model 15 paths in the network, using them to serve the four OD pairs shown in Figs 11 and 14. By comparing the 15 path travel times optimized by the proposed algorithms, given in Fig 12 and Table 8, with those optimized by the previous method, given in Table 7, the algorithms reduce the total network travel times of all modeled paths (= 156.2 pcu-h/h in Table 8) more than the previous method (= 159.0 pcuh/h in Table 7).
With respect to the optimized network configuration shown in Fig 8, the four shared lane markings demonstrate that the proposed lane-based optimization method is able to use all approach lanes effectively and optimize the lane turning patterns by fixing proper shared lane locations. In practical designs, engineers and network planners may design simple network configurations without involving shared lane markings. A sensitivity analysis compares the five performance indicators with respect to the optimized network configuration in Fig 8 and the simple network design (by removing the four shared lane markings from the optimized network configuration). The summary of those comparisons is given in Table 9. The proposed lane-based optimization method, which optimizes the network configuration with shared lane marking designs, outperforms the simple network design.

Conclusions
This study's main innovations are the extension of the lane-based design concept to the optimization of individual lane usage, which includes lane markings, assigned lane flows, and the Lengths of green bands are longer in the optimized network settings # Four optimized shared lane markings are given in the optimized network settings � Degree of saturation is below 0.9 for all approach traffic lanes in the optimized network settings https://doi.org/10.1371/journal.pone.0216958.t009 Optimization of signalized network configurations using the Lane-based method traffic signal settings used to configure a signalized network. The lane markings at individual intersections are defined as binary type variables and new governing constraint sets are developed them as network links. Available paths can then be established to serve the users' input OD flows. Efficient network link connections can obtained if proper lane markings are optimized. Multiple paths can also be generated to connect the origin and destination pairs through different turns at different intersections, according to the optimized lane marking patterns. Link connectivity and network configuration are enhanced and optimized in a unified framework that combines assigned lane flows and traffic signal settings. The proposed bandwidth optimization method facilitates the coordination of the traffic signal settings that control individual signal-controlled intersections. The offsets in regulating green start times across upstream and downstream intersections are also optimized. All of these are new features for extending the lane-based method to the design of signalized networks. The problem is formulated as a BMILP that can solved by a standard branch-and-bound routine. A numerical example with four intersections in a network with two approach lanes and two exit lanes is given to demonstrate the effectiveness of the proposed algorithms. Five performance indicators have been extracted from the optimization results and show that the proposed algorithms outperform the previous method that was developed to optimize the common flow multiplier, but without considering the bandwidth optimization for traffic signal coordination. To summarize, the optimized bandwidth is found to be 32.9 seconds and the total network travel time is 156.2 pcu-h/h. Future research on network configuration designs should incorporate safety performance to guard against shared lane markings on non-nearside lanes permitting three movement turns that may lead to statistically higher accident rates. In designing large-size networks, such restrictions should be imposed. The potential tradeoffs between engineering performance and safety considerations also merit future research attention.