add check that polygons dont overlap

hydromt_fiat/workflows/aggregation_areas.py

@@ -1,6 +1,7 @@
 import geopandas as gpd
 from typing import List, Union
 from pathlib import Path
+import sys 
 
 def process_value(value):
     if isinstance(value, list) and len(value) == 1:
@@ -43,6 +44,7 @@ def join_exposure_aggregation_multiple_areas(
 
 
         assert attribute_name in aggregation_gdf.columns, f"Attribute {attribute_name} not found in {file_path}"
+
 
         # If you overwrite the exposure_gdf with the joined data, you can append all 
         # aggregation areas to the same exposure_gdf
@@ -52,19 +54,9 @@ def join_exposure_aggregation_multiple_areas(
             op="intersects",
             how="left",
         )
+
+        assert exposure_gdf["Object ID"].is_unique, "Error! Polygons overlap! Please clean your data from overlaping features."
 
-        # aggregate the data if duplicates exist
-        aggregated = (
-            exposure_gdf.groupby("Object ID")[attribute_name].agg(list).reset_index()
-        )
-        exposure_gdf.drop_duplicates(subset="Object ID", keep="first", inplace=True)
-        exposure_gdf.drop(columns=attribute_name, inplace=True)
-        exposure_gdf = exposure_gdf.merge(aggregated, on="Object ID")
-
-        # Create a string from the list of values in the duplicated aggregation area 
-        # column
-        exposure_gdf[attribute_name] = exposure_gdf[attribute_name].apply(process_value)
-
         # Rename the 'aggregation_attribute' column to 'new_column_name'. Put in 
         # Documentation that the order the user put the label name must be the order of the gdf
         exposure_gdf.rename(columns={attribute_name: f"Aggregation Label: {label_name}"}, inplace=True)