st_coordinates()
sf
inla()
library('spocc') df <- occ(query = "Bradypus variegatus", from = "gbif", date = c("2000-01-01", "2019-12-31"), gbifopts = list(country = "CR"), has_coords = TRUE, limit = 1000) d <- occ2df(df)
library(sf) d <- st_as_sf(d, coords = c("longitude", "latitude")) st_crs(d) <- 4326 mapview(d)