terrain features

Author: mpcsb

.gitignore

@@ -13,12 +13,8 @@
 
 elevation_data/
 
-wfire/pkg/
+pkg/
 
-wfire/pkg/sumdb/
+bin/
 
-wfire/src/bin/
 
-wfire/src/cmd/run.exe
-
-wfire/bin/

src/simulation/terrain/HGT_parser.py

@@ -3,13 +3,13 @@
 import gzip
 import argparse
 import csv 
-# import numpy as np
+
 from tqdm import tqdm
 from gmalthgtparser import HgtParser   
 
 try: 
     parser = argparse.ArgumentParser() 
-    parser.add_argument('-wdir', type=str, default=r'/home/miguel/Documents/projects', help="working dir")
+    parser.add_argument('-wdir', type=str, default=r'/home/miguel/Documents/projects/Wildfire', help="working dir")
     parser.add_argument('lat1', type=float, help="latitude1 for HGT file")
     parser.add_argument('lon1', type=float, help="longitude1 for HGT file")
     parser.add_argument('lat2', type=float, help="latitude2 for HGT file")
@@ -32,6 +32,7 @@
 def linspace(lower, upper, num):  
     return [lower + x*(upper-lower)/num for x in range(num)]
 
+
 def filename_gen(lat, lon):
     ''' receives coordinates used to select hgt file'''
     if lat >= 0: 
@@ -101,6 +102,6 @@ def generate_topology(p1, p2, n_points=100):
 
 terrain = generate_topology(p1, p2, n_points=80)
 
-with open(r'C:\Users\Miguel\Documents\repos\wfire\wfire\src\simulation\terrain\temp\coords.csv', "w", newline="") as f:
+with open(r'/home/miguel/Documents/projects/Wildfire/wfire/src/simulation/terrain/temp/coords.csv', "w", newline="") as f:
     writer = csv.writer(f)
     writer.writerows(terrain)

src/simulation/terrain/example_call.md

@@ -0,0 +1 @@
+/home/miguel/Documents/projects/Wildfire/bin/python3.8 src/simulation/terrain/generate_objects.py 38.798802  -9.332599  38.799060 -9.350210

src/simulation/terrain/generate_objects.py

@@ -0,0 +1,80 @@
+import os
+import argparse
+import csv 
+
+import overpy 
+
+
+try: 
+    parser = argparse.ArgumentParser() 
+    parser.add_argument('-wdir', type=str, default=r'/home/miguel/Documents/projects/Wildfire', help="working dir")
+    parser.add_argument('lat1', type=float, help="latitude1")
+    parser.add_argument('lon1', type=float, help="longitude1")
+    parser.add_argument('lat2', type=float, help="latitude2")
+    parser.add_argument('lon2', type=float, help="longitude2")
+ 
+    args, unknown = parser.parse_known_args()
+
+    os.chdir(args.wdir) 
+except Exception as e:
+    print(e)
+lat1 = args.lat1
+lon1 = args.lon1
+lat2 = args.lat2
+lon2 = args.lon2
+# print(lat1,lat2, lon1, lon2)
+ 
+api = overpy.Overpass()
+
+q = f"""
+    way({lat1}, {lon1}, {lat2}, {lon2}) ["highway"];
+    (._;>;);
+    out body;
+    """ 
+result = api.query(q) 
+
+street_lst = list()
+for way in result.ways: 
+    for node in way.nodes:
+        street_lst.append((float(node.lat), float(node.lon))) 
+  
+with open(r'/home/miguel/Documents/projects/Wildfire/wfire/src/simulation/terrain/temp/street_coordinates.csv', "w", newline="") as f:
+    writer = csv.writer(f)
+    writer.writerows(street_lst)
+
+
+api = overpy.Overpass()
+q = f"""
+    way({lat1}, {lon1}, {lat2}, {lon2}) ["landuse"];
+    (._;>;);
+    out body;
+    """ 
+result = api.query(q) 
+
+natural_lst = list()
+for way in result.ways: 
+    for node in way.nodes:
+        natural_lst.append((float(node.lat), float(node.lon))) 
+  
+with open(r'/home/miguel/Documents/projects/Wildfire/wfire/src/simulation/terrain/temp/landuse_coordinates.csv', "w", newline="") as f:
+    writer = csv.writer(f)
+    writer.writerows(natural_lst)
+
+
+api = overpy.Overpass()
+q = f"""
+    way({lat1}, {lon1}, {lat2}, {lon2}) ["natural"];
+    (._;>;);
+    out body;
+    """ 
+result = api.query(q) 
+
+natural_lst = list()
+for way in result.ways: 
+    for node in way.nodes:
+        natural_lst.append((float(node.lat), float(node.lon))) 
+  
+with open(r'/home/miguel/Documents/projects/Wildfire/wfire/src/simulation/terrain/temp/natural_coordinates.csv', "w", newline="") as f:
+    writer = csv.writer(f)
+    writer.writerows(natural_lst)
+