address #375 (#376)

mapme-initiative · Oct 1, 2024 · 9ba6084 · 9ba6084
1 parent 840cbac
commit 9ba6084
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 1 deletion.
diff --git a/NEWS.md b/NEWS.md
@@ -4,6 +4,8 @@
 
 - fixes `portfolio_wide()` throwing an error when single assets with `NULL` 
   values are present 
+- `calc_mangroves_area()` returned NULL if invalid geometries were encountered
+  Now it tries to repair geometries and return the area of valid geometries (#375)
 
 # mapme.biodiversity 0.9.1
 

diff --git a/R/calc_mangroves_area.R b/R/calc_mangroves_area.R
@@ -53,7 +53,10 @@ calc_mangroves_area <- function() {
     }
 
     results <- purrr::map(1:length(gmw), function(j) {
-      intersected <- suppressWarnings(st_intersection(gmw[[j]], x))
+      mangroves <- st_make_valid(gmw[[j]])
+      intersected <- suppressWarnings(st_intersection(mangroves, x))
+      intersected <- st_make_valid(intersected)
+      intersected <- intersected[st_is_valid(intersected), ]
       area <- sum(as.numeric(st_area(intersected)), na.rm = TRUE) / 10000
       year <- strsplit(names(gmw[j]), "_|.gpkg")[[1]][3]