Fix missing clip in ElevationPoints function

rosepearson · Jan 26, 2025 · fb06004 · fb06004
1 parent 022122c
commit fb06004
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/src/geofabrics/geometry.py b/src/geofabrics/geometry.py
@@ -644,6 +644,7 @@ def _set_up(
             crs=polygon_list[0].crs,
         )
         polygon = polygon.to_crs(self.catchment_geometry.crs["horizontal"])
+        polygon = polygon.clip(self.catchment_geometry.catchment, keep_geom_type=True)
         points = points.clip(polygon.buffer(0), keep_geom_type=True)
         points = points.clip(self.catchment_geometry.catchment, keep_geom_type=True)
         points = points.reset_index(drop=True)