updated write_geoms method to check if region is available (#389)

added part to check if region is available
Deltares · Oct 4, 2024 · eb074e5 · eb074e5
1 parent f6932b3
commit eb074e5
Showing 1 changed file with 8 additions and 7 deletions.
diff --git a/hydromt_fiat/fiat.py b/hydromt_fiat/fiat.py
@@ -1328,13 +1328,14 @@ def write_geoms(self):
                     fn.format(name=name),
                 )
                 # At the very last clip based on the region
-                idx = gpd.sjoin(
-                    self.region.to_crs(geom.crs),
-                    geom,
-                    predicate="contains_properly",
-                    how="inner",
-                ).index_right
-                geom = geom.loc[idx]
+                if not self.region.empty:
+                    idx = gpd.sjoin(
+                        self.region.to_crs(geom.crs),
+                        geom,
+                        predicate="contains_properly",
+                        how="inner",
+                    ).index_right
+                    geom = geom.loc[idx]
                 geom.to_file(_fn)
 
         if self.geoms: