read osm data to images and graphs based on admin boundaries

Author: ishaan95

.gitignore

@@ -1,2 +1,2 @@
-data/output/*
+data/*
 __pycache__/*

osm2graphdata.py

@@ -0,0 +1,14 @@
+from pyrosm import OSM, get_data
+import networkx as nx
+
+# Initialize reader
+dir_name = "data/osm/"
+map_name = "SantaBarbara"
+osm = OSM(dir_name+map_name+".osm.pbf")
+for i in range(len(osm.get_boundaries()['osm_type'])):
+    if(osm.get_boundaries()['osm_type'][i]=='relation'): 
+        osm_data = OSM(dir_name+map_name+".osm.pbf", osm.get_boundaries()['geometry'][i])
+        nodes, edges = osm_data.get_network(nodes=True, network_type="driving")
+        nx_graph = osm.to_graph(nodes, edges, graph_type="networkx")
+        nx.write_gpickle(nx_graph, 'data/output/'+map_name+str(i)+'.gpickle')
+        

osm2imgdata.py

@@ -0,0 +1,20 @@
+# Program to return images and graphs given an OSM PBF file.
+from pyrosm import OSM, get_data
+import matplotlib.pyplot as plt
+import osmnx as ox
+
+# Initialize reader
+dir_name = "data/osm/"
+map_name = "SanJose"
+osm = OSM(dir_name+map_name+".osm.pbf")
+for i in range(len(osm.get_boundaries()['osm_type'])):
+    if(osm.get_boundaries()['osm_type'][i]=='relation'): 
+        osm_data = OSM(dir_name+map_name+".osm.pbf", osm.get_boundaries()['geometry'][i])
+        net = osm_data.get_network(network_type="driving")
+        if net is not None:
+            net.plot(linewidth=0.5)
+            plt.axis('off')
+            plt.savefig('data/output/'+map_name+str(i)+'.png', bbox_inches='tight', pad_inches=0, dpi=400)
+
+
+

readmbtiles.py

@@ -0,0 +1,3 @@
+import mbutil
+
+mbutil.mbtiles_to_disk('data/input/mbtiles/PaloAlto-mbtiles-openmaptiles/PaloAlto.mbtiles', 'data/output/tiles')

sliceimages.py

@@ -0,0 +1,15 @@
+#Script to slice up image files as needed.
+
+from PIL import Image
+import numpy as np
+from matplotlib import image as mpimg
+
+np_image = mpimg.imread('data/output/SanJose3.png')
+i=0
+
+
+    
+    
+
+
+