used unary_union and convex_hull from shapely to process multipolygon…

… geojson

src/backend/app/projects/project_routes.py

@@ -39,6 +39,8 @@
 from osm_fieldwork.make_data_extract import getChoices
 from osm_fieldwork.xlsforms import xlsforms_path
 from osm_rawdata.postgres import PostgresClient
+from shapely.geometry import mapping, shape
+from shapely.ops import unary_union
 from sqlalchemy.orm import Session
 
 from ..central import central_crud
@@ -665,12 +667,28 @@ async def get_data_extracts(
     # Validatiing Coordinate Reference System
     check_crs(boundary)
     xlsform = f"{xlsforms_path}/{category}.xls"
-
     config_path = f"{data_models_path}/{category}.yaml"
 
+    # Convert each feature into a Shapely geometry
+    geometries = [shape(feature["geometry"]) for feature in boundary["features"]]
+
+    # Merge the geometries into a single geometry (as a MultiPolygon)
+    merged_geometry = unary_union(geometries)
+
+    # Convert the merged MultiPolygon to a single Polygon using convex hull
+    merged_polygon = merged_geometry.convex_hull
+
+    # Convert the merged polygon back to a GeoJSON-like dictionary
+    boundary = {
+        "type": "Feature",
+        "geometry": mapping(merged_polygon),
+        "properties": {},
+    }
+
     # # OSM Extracts using raw data api
     pg = PostgresClient("underpass", config_path)
     data_extract = pg.execQuery(boundary)
+    log.info("Data extracts process completed")
     filter = FilterData(xlsform)
 
     updated_data_extract = {"type": "FeatureCollection", "features": []}