install.packages("rgbif")

library(dplyr)

library(purrr)

library(readr)  

library(magrittr) 

library(rgbif) 

library(taxize)

#####读取物种信息

gbif_taxon_keys <- 

readr::read_csv("D:/research data/species_list.csv") %>% 

pull("Taxon name") %>%

taxize::get_gbifid_(method="backbone") %>% 

imap(~ .x %>% mutate(original_sciname = .y)) %>%

bind_rows() %T>%

filter(matchtype == "EXACT" & status == "ACCEPTED") %>%

filter(kingdom == "Plantae") %>%

pull(usagekey)

#####筛选条件

occ_download(

pred_in("taxonKey", gbif_taxon_keys),

pred_in("basisOfRecord", c('PRESERVED_SPECIMEN','HUMAN_OBSERVATION','OBSERVATION','MACHINE_OBSERVATION')),

pred("country", "CN"),

pred("hasCoordinate", TRUE),

pred("hasGeospatialIssue", FALSE),

pred_gte("year", 1990),

format = "SIMPLE_CSV",

user="GBIF用户名", pwd="GBIF用户密码", email="GBIF中填写的邮箱地址"

)

######物种信息格式如下