update bike speed and network filter for new bike network

input/v6.4/berlin-v6.4-vehicleTypes.xml

@@ -24,9 +24,9 @@
 		<width meter="1.0"/>
 
 		<!-- The speed was set such that the obtained avg. speed per trip matches the one reported in the SrV survey -->
-		<!-- Note that this speed includes all the effects of delays that might occur, such as waiting traffic lights -->
-		<!-- 11.5 kmh -->
-		<maximumVelocity meterPerSecond="3.28"/>
+		<!-- Note that this speed includes all the effects of delays that might occur, such as waiting for traffic lights -->
+		<!-- The reported data is 10.29 km/h, it needs to set higher here (+0.5) because this includes access/egress as well -->
+		<maximumVelocity meterPerSecond="3"/>
 		<passengerCarEquivalents pce="0.2"/>
 		<networkMode networkMode="bike"/>
 	</vehicleType>

input/v6.4/berlin-v6.4.config.xml

@@ -80,11 +80,12 @@
 	-->
 
 	<module name="simwrapper" >
+		<param name="sampleSize" value="0.1" />
+
 		<parameterset type="params" >
 			<param name="context" value="" />
 			<param name="mapCenter" value="13.39,52.51" />
 			<param name="mapZoomLevel" value="9.1" />
-			<param name="sampleSize" value="0.1" />
 			<param name="shp" value="./area/area.shp" />
 		</parameterset>
 	</module>

src/main/java/org/matsim/prepare/counts/CreateCountsFromVMZ.java

@@ -94,6 +94,8 @@ private void matchWithNetwork(Path network) throws TransformException, IOExcepti
 			manager.addLinkFilter(l -> !l.getId().toString().startsWith("pt_"));
 			manager.addLinkFilter(l -> !NetworkUtils.getHighwayType(l).equals("service"));
 			manager.addLinkFilter(l -> !NetworkUtils.getHighwayType(l).equals("living_street"));
+			manager.addLinkFilter(l -> !NetworkUtils.getHighwayType(l).equals("path"));
+			manager.addLinkFilter(l -> !NetworkUtils.getHighwayType(l).equals("cycleway"));
 
 			net = manager.applyFilters();
 		}

src/main/python/extract_bike_speeds.py

@@ -0,0 +1,31 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+
+import numpy as np
+import os
+from matsim.scenariogen.data import EconomicStatus, TripMode, preparation, run_create_ref_data
+
+
+from extract_ref_data import trip_filter
+
+def person_filter(df):
+    df = df[df.reporting_day <= 4]
+    df = df[df.location == "Berlin"]
+
+    return df
+
+if __name__ == "__main__":
+    d = os.path.expanduser("~/Development/matsim-scenarios/shared-svn/projects/matsim-berlin/data/SrV/")
+
+    result = run_create_ref_data.create(
+        d + "Berlin+Umland",
+        person_filter, trip_filter,
+        run_create_ref_data.InvalidHandling.REMOVE_TRIPS,
+    )
+
+    t = result.trips
+
+    t["speed"] = (t.gis_length * 3600) / (t.duration * 60)
+    t = t[t.main_mode == TripMode.BIKE]
+
+    t.to_csv('bike_speeds.csv', columns=["t_weight", "age", "speed"], index=False)