remove openstreetmaps dependency

hub/helpers/geometry_helper.py

@@ -6,9 +6,10 @@ Project Coder Guille Gutierrez guillermo.gutierrezmorote@concordia.ca
 Code contributors: Pilar Monsalvete Alvarez de Uribarri pilar.monsalvete@concordia.ca
 """
 import math
+from pathlib import Path
 
 import numpy as np
-import requests
+from PIL import Image
 from trimesh import Trimesh
 from trimesh import intersections
 
@@ -79,14 +80,15 @@ class GeometryHelper:
     y = int((city.upper_corner[1] - city.lower_corner[1]) * 0.5) + 1
     city_map = [['' for _ in range(y + 1)] for _ in range(x + 1)]
     map_info = [[{} for _ in range(y + 1)] for _ in range(x + 1)]
-    # img = Image.new('RGB', (x + 1, y + 1), "black")  # create a new black image
-    # city_image = img.load()  # create the pixel map
+    img = Image.new('RGB', (x + 1, y + 1), "black")  # create a new black image
+    city_image = img.load()  # create the pixel map
     for building_name in building_names:
       building = city.city_object(building_name)
       line = 0
       for ground in building.grounds:
         length = len(ground.perimeter_polygon.coordinates) - 1
         for i, coordinate in enumerate(ground.perimeter_polygon.coordinates):
+
           j = i + 1
           if i == length:
             j = 0
@@ -106,7 +108,7 @@ class GeometryHelper:
                 'line_start': (coordinate[0], coordinate[1]),
                 'line_end': (next_coordinate[0], next_coordinate[1]),
               }
-              # city_image[x, y] = (100, 0, 0)
+              city_image[x, y] = (100, 0, 0)
             elif city_map[x][y] != building.name:
               neighbour = city.city_object(city_map[x][y])
               neighbour_info = map_info[x][y]
@@ -166,8 +168,51 @@ class GeometryHelper:
               elif building not in neighbour.neighbours:
                 neighbour.neighbours.append(building)
           line += 1
-    # if plot:
-      # img.show()
+
+    if plot:
+      img.show()
+    return lines_information
+
+  @staticmethod
+  def fast_city_mapping(city, building_names=None):
+    lines_information = {}
+    if building_names is None:
+      building_names = [b.name for b in city.buildings]
+    x = int((city.upper_corner[0] - city.lower_corner[0]) * 0.5) + 1
+    y = int((city.upper_corner[1] - city.lower_corner[1]) * 0.5) + 1
+    city_map = [['' for _ in range(y + 1)] for _ in range(x + 1)]
+    for building_name in building_names:
+      building = city.city_object(building_name)
+      line = 0
+      for ground in building.grounds:
+        length = len(ground.perimeter_polygon.coordinates) - 1
+        for i, coordinate in enumerate(ground.perimeter_polygon.coordinates):
+          j = i + 1
+          if i == length:
+            j = 0
+          next_coordinate = ground.perimeter_polygon.coordinates[j]
+          point = GeometryHelper.coordinate_to_map_point(coordinate, city)
+          distance = int(GeometryHelper.distance_between_points(coordinate, next_coordinate))
+          if distance == 0:
+            continue
+          delta_x = (coordinate[0] - next_coordinate[0]) / (distance / 0.5)
+          delta_y = (coordinate[1] - next_coordinate[1]) / (distance / 0.5)
+          for k in range(0, distance):
+            x = MapPoint(point.x + (delta_x * k), point.y + (delta_y * k)).x
+            y = MapPoint(point.x + (delta_x * k), point.y + (delta_y * k)).y
+            if city_map[x][y] == '':
+              city_map[x][y] = building.name
+            elif city_map[x][y] != building.name:
+              neighbour = city.city_object(city_map[x][y])
+              if building.neighbours is None:
+                building.neighbours = [neighbour]
+              elif neighbour not in building.neighbours:
+                building.neighbours.append(neighbour)
+              if neighbour.neighbours is None:
+                neighbour.neighbours = [building]
+              elif building not in neighbour.neighbours:
+                neighbour.neighbours.append(building)
+          line += 1
     return lines_information
 
   @staticmethod
@@ -254,19 +299,24 @@ class GeometryHelper:
     """
     Get Location from latitude and longitude
     """
-    url = 'https://nominatim.openstreetmap.org/reverse?lat={latitude}&lon={longitude}&format=json'
-    response = requests.get(url.format(latitude=latitude, longitude=longitude))
-    if response.status_code != 200:
-      # This means something went wrong.
-      raise Exception('GET /tasks/ {}'.format(response.status_code))
-
-    response = response.json()
+    _data_path = Path(Path(__file__).parent.parent / 'data/geolocation/cities15000.txt').resolve()
+    latitude = float(latitude)
+    longitude = float(longitude)
+    distance = math.inf
+    country = 'Unknown'
     city = 'Unknown'
-    country = 'ca'
-    if 'city' in response['address']:
-      city = response['address']['city']
-    if 'country_code' in response['address']:
-      country = response['address']['country_code']
+    with open(_data_path, 'r') as f:
+      for line_number, line in enumerate(f):
+        fields = line.split('\t')
+        file_city_name = fields[2]
+        file_latitude = float(fields[4])
+        file_longitude = float(fields[5])
+        file_country_code = fields[8]
+        new_distance = math.sqrt(pow((latitude - file_latitude), 2) + pow((longitude - file_longitude), 2))
+        if distance > new_distance:
+          distance = new_distance
+          country = file_country_code
+          city = file_city_name
     return Location(country, city)
 
   @staticmethod