fix: cell formula again, now final

server/map_handler/add_new_lake.py

@@ -33,17 +33,17 @@ def cut_map(self, cursor, lake_name: str, cell_size_in_km: float = 0.5):
         # List to store all the map tiles
         divided_map = []
 
-        # Convert the cell size to degrees
-        cell_size = cell_size_in_km * 2
-        cell_width = cell_size / 111.32
-        # A slightly more complicated formula is required to calculate the height. This ensures that
-        # the height in km is equal to the width in km regardless of the latitude.
-        cell_height = (cell_size * 360) / (40075 * cos(cell_width))
+        # Select an arbitrary x and y value from within the polygon
+        bounds = polygons[0].bounds
+        start_x, start_y, _, _ = bounds
+
+        cell_width = cell_size_in_km / 111.3200
+        cell_height = cell_width / cos(start_x * 0.01745)
 
         # Process all polygons
         for polygon in polygons:
             # Generate a grid based on the calculated cell size
-            lines = create_grid(polygon, cell_width, cell_height)
+            lines = create_grid(polygon, cell_width*2, cell_height)
             lines.append(polygon.boundary)
             # Merge the grid lines into a single grid object
             lines = unary_union(lines)
@@ -233,7 +233,6 @@ def plot_map(divided_map):
 
     # Configure plot settings
     fig, ax = plt.subplots()
-    ax.set_aspect(1.5)
 
     # Plot each tile
     for tile in tiles:
@@ -241,5 +240,7 @@ def plot_map(divided_map):
         random_color = "#{:06x}".format(random.randint(0, 0xFFFFFF))
         gpd.GeoSeries(tile.geometry).plot(ax=ax, facecolor=random_color, edgecolor='none')
 
+    ax.set_aspect(1.6)
+
     # Display plot
     plt.show()