We should have a new table with distribution per Municipio.

[dlizcano]

We should have a new table with distribution per Municipio.

library (mammalcol)
library (rgbif)
library (sf)

get map of municipios

munimapcol <- sf::st_as_sf(gadm(country="COL", level=2, path=tempdir()))

Crea lista vacia

gbif_data_sp <- list()

Descargar datos de ocurrencia de GBIF para cada especie; Este proceso puede tomar tiempo si existen muchos puntos!

loop for(i in 1:length(taxon$scientificName)){

gbif_data_sp[[i]] <- occ_data(scientificName = taxon$scientificName[i],#"Mazama sanctaemartae",
hasCoordinate = TRUE,
country ="CO",
#basisOfRecord= c("PRESERVED_SPECIMEN", "MACHINE_OBSERVATION"),
limit = 20000)
}

removes the metadata

gbif_data_sp_list <- list() # empty list
for (i in 1:length(taxon$scientificName)){
gbif_data_sp_list[[i]] <- gbif_data_sp[[i]]$data
}

convert to dataframe

df <- ldply(gbif_data_sp_list, data.frame)