fix_invalid_geometries.R

fix_invalid_geometries <- function(data) {

  is_valid <- st_is_valid(data)

  invalid_geom <- data[!is_valid, ]

  geom_valid <- data[is_valid, ]

  sf_use_s2(TRUE)

  geom_invalid_fixed <-
    invalid_geom %>%
    s2::as_s2_geography(check = F) %>%
    s2::s2_union() %>%
    s2::s2_rebuild(s2::s2_options(split_crossing_edges = TRUE)) %>%
    st_as_sf() %>%
    bind_cols(invalid_geom %>% st_drop_geometry())

  geom_fixed <-
    bind_rows(geom_valid, geom_invalid_fixed)

  return(geom_fixed)

}