Search MarineGazetteer for polygon overlap

gazetteerSearch

library(sf)
library(dplyr)

#get WKT
bb_points <- event %>% 
  st_as_sf(coords = c("decimalLongitude", "decimalLatitude"), 
           crs = "WGS84") %>% 
  st_bbox() %>% 
  st_as_sfc() %>% 
  st_cast(to = "POINT")

lapply(bb_points, mregions2::gaz_search) %>% 
  data.table::rbindlist() %>% 
  distinct() %>% 
  filter(preferredGazetteerNameLang =="English" & status == "standard")

Righteous. Thanks for sharing!