Add support for parking time function based on the destination zone_data

Scripts/parameters/assignment.py

@@ -2,6 +2,11 @@
 
 from collections import namedtuple
 from typing import Dict, List, Union
+
+import numpy as np
+import pandas as pd
+
+from datahandling.zonedata import ZoneData
 RoadClass = namedtuple(
     "RoadClass",
     (
@@ -171,6 +176,11 @@
 }
 # Default distance unit cost [eur/km]
 dist_unit_cost = 0.12
+# Parking time function
+def parking_time(zone_data: ZoneData) -> pd.Series:
+    density: pd.Series = (zone_data['population'] + zone_data['workplaces']) / zone_data['zone_area']
+    return 0.05993817*np.sqrt(density) + 5.24176150
+
 # Default distance unit time for trucks and trailer trucks [min/km]
 freight_dist_unit_time = 0.2
 # Boarding penalties for different transit modes

Scripts/transform/impedance_transformer.py

@@ -3,6 +3,10 @@
 
 import parameters.impedance_transformation as param
 from parameters.assignment import assignment_classes
+try:
+    from parameters.assignment import parking_time
+except ImportError:
+    parking_time = None
 
 
 class ImpedanceTransformer:
@@ -55,4 +59,11 @@ def transform(self, purpose, impedance):
         for i in range(1, len(purpose.sub_bounds)):
             trips_per_month[purpose.sub_bounds[i], :] = trips_month[i]
         day_imp["transit"]["cost"] /= trips_per_month
+        # Add parking time to car matrices
+        if parking_time is not None:
+            ptime = parking_time(purpose.zone_data).to_numpy()[cols]
+            ptime = numpy.clip(ptime, 0, 30)
+            car_modes = [k for k in day_imp.keys() if k in ('car', 'car_work', 'car_transit')]
+            for k in car_modes:
+                day_imp[k]['time'] += ptime[None,:]
         return day_imp