R工具经纬度提取海拔 rgbif

一 前期准备

1 参考网址 https://rdrr.io/cran/rgbif/man/elevation.html

2 Rverson 4.0.2

3 提前注册  http://www.geonames.org/login，得到用户名XXX

二 安装需要的包

install.packages(c("here","sp","rgbif"))

library(rgbif)

library(here)

library(sp)

三 读取经纬度数据

geo=read.csv("Farjon.csv",header=T)

geo=subset(geo,!is.na(Latitude))

head(geo)

#   taxon Latitude Longitude Elevation

#1 A_alba 50.35000 16.883333         0

# 2 A_alba 49.13333 22.150000         0

 #3 A_alba 50.88333 18.433333         0

四 提取海拔

options(geonamesUsername = "XXX")

file.edit(here::here(".Rprofile"))

user <- "XXX" #需注册，见前期准备

lats <- geo$Latitude

lons <- geo$Longitude

abie_ele=elevation(latitude = lats, longitude = lons, elevation_model = "srtm1",username = user)  #srtm1=30m*30m，最高分辨率。