Merge branch 'Axel' of https://github.com/Ahp06/SUMO_Emissions

sumo_project/actions.py

@@ -13,14 +13,19 @@ from model import Area, Vehicle
 
 
 def compute_edge_weight(edge_id):
-    return (traci.edge.getCOEmission(edge_id)
-            + traci.edge.getNOxEmission(edge_id)
-            + traci.edge.getHCEmission(edge_id)
-            + traci.edge.getPMxEmission(edge_id)
-            + traci.edge.getCO2Emission(edge_id))/(traci.edge.getLaneNumber(edge_id))
     
-def adjust_edges_weights():     
-    for edge_id in traci.edge.getIDList():
+    co2 = traci.edge.getCO2Emission(edge_id)
+    co = traci.edge.getCOEmission(edge_id)
+    nox = traci.edge.getNOxEmission(edge_id)
+    hc = traci.edge.getHCEmission(edge_id)
+    pmx = traci.edge.getPMxEmission(edge_id)
+    
+    return (co2 + co + nox + hc + pmx)
+            
+def adjust_edges_weights(area): 
+    area.weight_adjusted = True    
+    for lane in area._lanes:
+        edge_id = traci.lane.getEdgeID(lane.lane_id)
         weight = compute_edge_weight(edge_id)  # by default edges weight = length/mean speed
         traci.edge.setEffort(edge_id, weight)
         

sumo_project/config.py

@@ -8,17 +8,17 @@ import logging
 import os
 import sys
 
+from model import Emission
+
 
 class Config:
-    
     # Total of emissions of all pollutants in mg for n steps of simulation without acting on areas
     # These constants are simulation dependant, you must change them according to your simulation 
-    total_emissions100 = 13615949.148296086
-    total_emissions200 = 43970763.15084738
-    total_emissions300 = 87382632.0821697
+    ref200 = Emission(co2=42816869.05436445, co=1128465.0343051048, nox=18389.648337283958, hc=6154.330914019103,
+                      pmx=885.0829265236318)
 
     def __init__(self):
-        '''Default constructor'''
+        """Default constructor"""
 
     def import_config_file(self, config_file):
         with open(config_file, 'r') as f:
@@ -56,16 +56,17 @@ class Config:
             self.weight_routing_mode = False
             self.lock_area_mode = False
 
-                
     def __repr__(self) -> str:
-        return (str(f'grid : {self.areas_number}x{self.areas_number}\n')
-        + str(f'step number = {self.n_steps}\n')
-        + str(f'window size = {self.window_size}\n')
-        + str(f'weight routing mode = {self.weight_routing_mode}\n')
-        + str(f'lock area mode = {self.lock_area_mode}\n')
-        + str(f'limit speed mode = {self.limit_speed_mode}, RF = {self.speed_rf*100}%\n')
-        + str(f'adjust traffic light mode = {self.adjust_traffic_light_mode} , RF = {self.trafficLights_duration_rf*100}%\n'))
-                
+        return (
+            f'grid : {self.areas_number}x{self.areas_number}\n'
+            f'step number = {self.n_steps}\n'
+            f'window size = {self.window_size}\n'
+            f'weight routing mode = {self.weight_routing_mode}\n'
+            f'lock area mode = {self.lock_area_mode}\n'
+            f'limit speed mode = {self.limit_speed_mode}, RF = {self.speed_rf * 100}%\n'
+            f'adjust traffic light mode = {self.adjust_traffic_light_mode},'
+            f'RF = {self.trafficLights_duration_rf * 100}%\n'
+        )
 
     def init_traci(self):
         if 'SUMO_HOME' in os.environ:
@@ -101,12 +102,6 @@ class Config:
 
         return logger
 
-    def get_basics_emissions(self):
-        if self.n_steps == 100:
-            return self.total_emissions100
+    def get_ref_emissions(self):
         if self.n_steps == 200:
-            return self.total_emissions200
-        if self.n_steps == 300:
-            return self.total_emissions300
-    
-    
+            return self.ref200

sumo_project/emissions.py

@@ -1,19 +1,21 @@
 import argparse
+import csv
+import datetime
+import itertools
+import os
 import sys
 import time
-from traci import trafficlight
-import traci
 from typing import List
 
+import traci
 from parse import search
 from shapely.geometry import LineString
 
 import actions
 from config import Config
-from model import Area, Vehicle, Lane , TrafficLight , Phase , Logic
+from model import Area, Vehicle, Lane, TrafficLight, Phase, Logic, Emission
 
-
-def init_grid(simulation_bounds, areas_number):
+def init_grid(simulation_bounds, areas_number, window_size):
     grid = list()
     width = simulation_bounds[1][0] / areas_number
     height = simulation_bounds[1][1] / areas_number
@@ -22,10 +24,10 @@ def init_grid(simulation_bounds, areas_number):
             # bounds coordinates for the area : (xmin, ymin, xmax, ymax)
             ar_bounds = ((i * width, j * height), (i * width, (j + 1) * height),
                          ((i + 1) * width, (j + 1) * height), ((i + 1) * width, j * height))
-            area = Area(ar_bounds)
-            area.name = 'Area ({},{})'.format(i, j)
+            name = 'Area ({},{})'.format(i, j)
+            area = Area(ar_bounds, name, window_size)
             grid.append(area)
-            traci.polygon.add(area.name, ar_bounds, (0, 255, 0))
+            traci.polygon.add(area.name, ar_bounds, (255, 0, 0))
     return grid
 
 
@@ -40,14 +42,16 @@ def get_all_lanes() -> List[Lane]:
 
 def parse_phase(phase_repr):
     duration = search('duration: {:f}', phase_repr)
-    minDuration = search('minDuration: {:f}', phase_repr)
-    maxDuration = search('maxDuration: {:f}', phase_repr)
-    phaseDef = search('phaseDef: {}\n', phase_repr)
+    min_duration = search('minDuration: {:f}', phase_repr)
+    max_duration = search('maxDuration: {:f}', phase_repr)
+    phase_def = search('phaseDef: {}\n', phase_repr)
 
-    if phaseDef is None: phaseDef = ''
-    else : phaseDef = phaseDef[0]
+    if phase_def is None:
+        phase_def = ''
+    else:
+        phase_def = phase_def[0]
 
-    return Phase(duration[0], minDuration[0], maxDuration[0], phaseDef)
+    return Phase(duration[0], min_duration[0], max_duration[0], phase_def)
 
 
 def add_data_to_areas(areas: List[Area]):
@@ -68,11 +72,13 @@ def add_data_to_areas(areas: List[Area]):
 
 
 def compute_vehicle_emissions(veh_id):
-    return (traci.vehicle.getCOEmission(veh_id)
-            +traci.vehicle.getNOxEmission(veh_id)
-            +traci.vehicle.getHCEmission(veh_id)
-            +traci.vehicle.getPMxEmission(veh_id)
-            +traci.vehicle.getCO2Emission(veh_id))
+    co2 = traci.vehicle.getCO2Emission(veh_id)
+    co = traci.vehicle.getCOEmission(veh_id)
+    nox = traci.vehicle.getNOxEmission(veh_id)
+    hc = traci.vehicle.getHCEmission(veh_id)
+    pmx = traci.vehicle.getPMxEmission(veh_id)
+
+    return Emission(co2, co, nox, hc, pmx)
 
 
 def get_all_vehicles() -> List[Vehicle]:
@@ -87,12 +93,12 @@ def get_all_vehicles() -> List[Vehicle]:
 
 def get_emissions(grid: List[Area], vehicles: List[Vehicle], current_step, config, logger):
     for area in grid:
-        vehicle_emissions = 0
+        total_emissions = Emission()
         for vehicle in vehicles:
             if vehicle.pos in area:
-                vehicle_emissions += vehicle.emissions
+                total_emissions += vehicle.emissions
 
-        area.emissions_by_step.append(vehicle_emissions)
+        area.emissions_by_step.append(total_emissions)
 
         if area.sum_emissions_into_window(current_step, config.window_size) >= config.emissions_threshold:
 
@@ -100,7 +106,8 @@ def get_emissions(grid: List[Area], vehicles: List[Vehicle], current_step, confi
                 logger.info(f'Action - Decreased max speed into {area.name} by {config.speed_rf * 100}%')
                 actions.limit_speed_into_area(area, vehicles, config.speed_rf)
                 if config.adjust_traffic_light_mode and not area.tls_adjusted:
-                    logger.info(f'Action - Decreased traffic lights duration by {config.trafficLights_duration_rf*100}%')
+                    logger.info(
+                        f'Action - Decreased traffic lights duration by {config.trafficLights_duration_rf * 100}%')
                     actions.adjust_traffic_light_phase_duration(area, config.trafficLights_duration_rf)
 
             if config.lock_area_mode and not area.locked:
@@ -108,11 +115,35 @@ def get_emissions(grid: List[Area], vehicles: List[Vehicle], current_step, confi
                     logger.info(f'Action - {area.name} blocked')
                     actions.lock_area(area)
 
-            traci.polygon.setColor(area.name, (255, 0, 0))
+            if config.weight_routing_mode and not area.weight_adjusted:
+                actions.adjust_edges_weights(area)
+
             traci.polygon.setFilled(area.name, True)
 
         else:
             actions.reverse_actions(area)
+            traci.polygon.setFilled(area.name, False)
+
+
+def get_reduction_percentage(ref, total):
+    return (ref - total) / ref * 100
+
+
+def export_data_to_csv(config, grid):
+    csv_dir = os.path.join(SCRIPTDIR, 'csv')
+    if not os.path.exists(csv_dir):
+        os.mkdir(csv_dir)
+        
+    now = datetime.datetime.now().strftime("%Y_%m_%d_%H_%M_%S")
+        
+    with open(f'csv/{now}.csv', 'w') as f:
+        writer = csv.writer(f)
+        # Write CSV headers
+        writer.writerow(itertools.chain(('Step',), (a.name for a in grid)))
+        # Write all areas emission value for each step
+        for step in range(config.n_steps):
+            em_for_step = (f'{a.emissions_by_step[step].value():.3f}' for a in grid)
+            writer.writerow(itertools.chain((step,), em_for_step))
 
 
 def run(config, logger):
@@ -123,66 +154,87 @@ def run(config, logger):
         logger.info('Loading data for the simulation')
         start = time.perf_counter()
 
-        grid = init_grid(traci.simulation.getNetBoundary(), config.areas_number)
+        grid = init_grid(traci.simulation.getNetBoundary(), config.areas_number, config.window_size)
         add_data_to_areas(grid)
 
         loading_time = round(time.perf_counter() - start, 2)
         logger.info(f'Data loaded ({loading_time}s)')
 
-        logger.info('Start of the simulation')
+        logger.info('Simulation started...')
         step = 0
         while step < config.n_steps:  # traci.simulation.getMinExpectedNumber() > 0:
             traci.simulationStep()
 
             vehicles = get_all_vehicles()
             get_emissions(grid, vehicles, step, config, logger)
-
-            if config.weight_routing_mode:
-                actions.adjust_edges_weights()
-
             step += 1
 
+            print(f'step = {step}/{config.n_steps}', end='\r')
+
     finally:
         traci.close(False)
+        export_data_to_csv(config, grid)
+
         simulation_time = round(time.perf_counter() - start, 2)
         logger.info(f'End of the simulation ({simulation_time}s)')
+        logger.info(f'Real-time factor : {config.n_steps / simulation_time}')
 
-        total_emissions = 0
+        total_emissions = Emission()
         for area in grid:
             total_emissions += area.sum_all_emissions()
 
-        logger.info(f'Total emissions = {total_emissions} mg')
+        logger.info(f'Total emissions = {total_emissions.value()} mg')
 
         if not config.without_actions_mode:
-            ref = config.get_basics_emissions()
+            ref = config.get_ref_emissions()
             if not (ref is None):
-                diff_with_actions = (ref - total_emissions) / ref    
-                logger.info(f'Reduction percentage of emissions = {diff_with_actions*100} %')
+                global_diff = (ref.value() - total_emissions.value()) / ref.value()
+
+                logger.info(f'Global reduction percentage of emissions = {global_diff * 100} %')
+                logger.info(f'-> CO2 emissions = {get_reduction_percentage(ref.co2, total_emissions.co2)} %')
+                logger.info(f'-> CO emissions = {get_reduction_percentage(ref.co, total_emissions.co)} %')
+                logger.info(f'-> Nox emissions = {get_reduction_percentage(ref.nox, total_emissions.nox)} %')
+                logger.info(f'-> HC emissions = {get_reduction_percentage(ref.hc, total_emissions.hc)} %')
+                logger.info(f'-> PMx emissions = {get_reduction_percentage(ref.pmx, total_emissions.pmx)} %')
+
+
+def add_options(parser):
+    parser.add_argument("-f", "--configfile", type=str, default='configs/default_config.json', required=False,
+                        help='Choose your configuration file from your working directory')
+    parser.add_argument("-save", "--save", action="store_true",
+                        help='Save the logs into the logs folder')
+    parser.add_argument("-steps", "--steps", type=int, default=200, required=False,
+                        help='Choose the simulated time (in seconds)')
+    parser.add_argument("-ref", "--ref", action="store_true",
+                        help='Launch a reference simulation (without acting on areas)')
+    parser.add_argument("-gui", "--gui", action="store_true",
+                        help="Set GUI mode")
 
 
 def main(args):
     parser = argparse.ArgumentParser(description="")
-    parser.add_argument("-f", "--configfile", type=str, default='configs/default_config.json', required=False)
-    parser.add_argument("-save", "--save", action="store_true")
-    parser.add_argument("-ref", "--ref", action="store_true")
+    add_options(parser)
     args = parser.parse_args(args)
 
-    # > py ./emissions.py -f configs/config1.json -save
-    # will load the configuration file "config1.json" and save logs into the logs directory 
-    
-    # > py ./emissions.py -f configs/config1.json -save -ref & py ./emissions.py -f configs/config1.json -save
-    # same as above but also launches a reference simulation by using -ref option 
-    
     config = Config()
     config.import_config_file(args.configfile)
     config.init_traci()
     logger = config.init_logger(save_logs=args.save)
+
     if args.ref:
         config.without_actions_mode = True
-        config.check_config()
         logger.info(f'Reference simulation')
-    logger.info(f'Loaded configuration file : {args.configfile}')
 
+    if args.steps:
+        config.n_steps = args.steps
+
+    if args.gui:
+        config._SUMOCMD = "sumo-gui"
+
+    config.check_config()
+
+    logger.info(f'Loaded configuration file : {args.configfile}')
+    logger.info(f'Simulated time : {args.steps}s')
     run(config, logger)
 
 

sumo_project/model.py

@@ -1,9 +1,10 @@
-from traci._trafficlight import Logic as SUMO_Logic
+import collections
 from typing import Tuple, Set
 
 from shapely.geometry import Point, LineString
 from shapely.geometry import Polygon
 from shapely.geometry.base import BaseGeometry
+from traci._trafficlight import Logic as SUMO_Logic
 
 
 class Lane:
@@ -17,6 +18,7 @@ class Lane:
         """Overrides the default implementation"""
         return hash(self.lane_id)
 
+
 class Phase:
     def __init__(self, duration: float, minDuration: float, maxDuration: float, phaseDef: str):
         self.duration = duration
@@ -28,11 +30,13 @@ class Phase:
         repr = f'Phase(duration:{self.duration},minDuration:{self.minDuration},maxDuration:{self.maxDuration},phaseDef:{self.phaseDef})'
         return str(repr)
 
+
 class Logic:
     def __init__(self, logic: SUMO_Logic, phases: Set[Phase]):
         self._logic = logic
         self._phases: Set[Phase] = phases
 
+
 class TrafficLight:
 
     def __init__(self, tl_id: str, logics: Set[Logic]):
@@ -43,15 +47,38 @@ class TrafficLight:
         """Overrides the default implementation"""
         return hash(self.tl_id)
 
+
+class Emission:
+    def __init__(self, co2=0, co=0, nox=0, hc=0, pmx=0):
+        self.co2 = co2
+        self.co = co
+        self.nox = nox
+        self.hc = hc
+        self.pmx = pmx
+
+    def __add__(self, other):
+        return Emission(self.co2 + other.co2, self.co + other.co, self.nox + other.nox, self.hc + other.hc,
+                        self.pmx + other.pmx)
+
+    def value(self):
+        return self.co2 + self.co + self.nox + self.hc + self.pmx
+
+    def __repr__(self) -> str:
+        repr = f'Emission(co2={self.co2},co={self.co},nox={self.nox},hc={self.hc},pmx={self.pmx})'
+        return str(repr)
+
+
 class Area:
 
-    def __init__(self, coords, name=''):
+    def __init__(self, coords, name, window_size):
         self.limited_speed = False
         self.locked = False
         self.tls_adjusted = False
+        self.weight_adjusted = False
         self.rectangle = Polygon(coords)
         self.name = name
         self.emissions_by_step = []
+        self.window = collections.deque(maxlen=window_size)
         self._lanes: Set[Lane] = set()
         self._tls: Set[TrafficLight] = set()
 
@@ -78,16 +105,18 @@ class Area:
         self._lanes.remove(lane)
 
     def sum_all_emissions(self):
-        sum = 0 
+        sum = Emission()
         for emission in self.emissions_by_step:
             sum += emission
         return sum
 
     def sum_emissions_into_window(self, current_step, window_size):
+
+        self.window.appendleft(self.emissions_by_step[current_step].value())
+
         sum = 0
-        q = current_step // window_size #Returns the integral part of the quotient
-        for i in range(q*window_size, current_step):
-            sum += self.emissions_by_step[i]
+        for i in range(self.window.__len__()):
+            sum += self.window[i]
         return sum
 
     @classmethod
@@ -98,10 +127,11 @@ class Area:
             (xmax, ymax),
             (xmax, ymin)))
 
+
 class Vehicle:
 
     def __init__(self, veh_id: int, pos: Tuple[float, float]):
-        self.emissions: float = 0.0
+        self.emissions: Emission = Emission()
         self.veh_id = veh_id
         self.pos = Point(pos)