从GBIF中获取研究区内物种经纬度信息

#参考网站：https://blog.sciencenet.cn/blog-3432920-1238640.html
install.packages(c("dplyr","purrr","readr","magrittr","taxize"))
library(dplyr)
library(purrr)
library(readr)  
library(magrittr) 
library(rgbif) 
library(taxize)

#####读取物种信息
gbif_taxon_keys <- 
  readr::read_csv("D:/R/plant.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") %>%   #kingdom填物种类型,动物为Animalia，植物为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),
  format = "SIMPLE_CSV",
  user="guohaotian", pwd="GUO1155667", email="18339128227@163.com"
)
###数据下载
##筛选条件运行后，最下面会出现一个网址：https://doi.org/10.15468/dl.zdqzft
  #在网址内下载物种数据