Updated Klimabäume Data for Würzburg

python_scripts/add_trees_to_open_street_map/add_trees.py

@@ -40,7 +40,12 @@ from export_OSM_roads_to_csv import write_roads_to_csv
 # From RegionChooser, or https://wktmap.com/
 # Case-study for Grombühl Paper. (4 Blocks):
 # WKT = "POLYGON((9.942183 49.804010, 9.942301 49.802740, 9.945477 49.802698, 9.945563 49.804014, 9.942167 49.804055, 9.942183 49.804010))"
-WKT = "POLYGON ((9.940352 49.804078, 9.942176 49.804646, 9.949343 49.804508, 9.953356 49.804383, 9.954686 49.804355, 9.955094 49.801932, 9.951661 49.800422, 9.949 49.800796, 9.94591 49.800713, 9.941339 49.801724, 9.94061 49.802527, 9.940352 49.804078))"
+# WKT = "POLYGON ((9.940352 49.804078, 9.942176 49.804646, 9.949343 49.804508, 9.953356 49.804383, 9.954686 49.804355, 9.955094 49.801932, 9.951661 49.800422, 9.949 49.800796, 9.94591 49.800713, 9.941339 49.801724, 9.94061 49.802527, 9.940352 49.804078))"
+# Grombühl Mini:
+# WKT = "POLYGON((9.946389 49.803571, 9.946378 49.803007, 9.948400 49.802990, 9.948433 49.803557, 9.946389 49.803571))"
+# Grafenbühl Mini:
+WKT = "POLYGON((9.151047 48.908034, 9.151348 48.907276, 9.153536 48.907417, 9.153413 48.908263, 9.151047 48.908034))"
+# WKT = "POLYGON((9.146925 48.909689, 9.148727 48.904838, 9.154628 48.904302, 9.153963 48.909844, 9.146925 48.909689))"
 # Grafenbühl:
 # WKT = "POLYGON((9.147551 48.908059, 9.148635 48.907953, 9.149525 48.907819, 9.151177 48.907819, 9.151413 48.907840, 9.153226 48.908087, 9.153387 48.906705, 9.149160 48.906634, 9.148999 48.906620, 9.147551 48.908059))"
 
@@ -62,7 +67,7 @@ EXISTING_TREES = 'existing_trees/baumkataster/Baum.shp'
 EPSG_ID = 25832
 
 # Trees will be planted every TREE_DISTANCE along roads:
-TREE_DISTANCE = 25  # [m]
+TREE_DISTANCE = 15  # [m]
 # Unless there's already another tree closer than MIN_DISTANCE away:
 MIN_DISTANCE = 10  # [m]
 
@@ -78,12 +83,13 @@ TREE_SPECIES = {
     'Acer platanoides': 40
 }
 
-TREE_CROWN_DIAMETER = [4, 6, 8]  # [m]
-TREE_HEIGHT = [8, 10, 12] # [m]
+TREE_CROWN_DIAMETER = [8]  # [m]
+TREE_HEIGHT = [12] # [m]
 TREE_TRUNK_DIAMETER = 90  # [cm]
 
 # Percentage of roads to which trees should be added
-ROAD_PERCENTAGE = 50
+# TODO: Length weighted?
+ROAD_PERCENTAGE = 100
 
 ########################
 #       STREETS        #
@@ -95,6 +101,7 @@ ROAD_PERCENTAGE = 50
 DEFAULT_WIDTHS = {
     'unclassified': 0,
     #   'residential': 8,
+    'footway': -1,
     'motorway': -1,
     'trunk': -1,
 }