|
刚审稿发现有个可以获取一个北美经纬度点的海拔R包“elevatr”
install.packages("elevatr")
library(elevatr)
mt_wash <- data.frame(x = -71.3036, y = 44.2700) mt_mans <- data.frame(x = -72.8145, y = 44.5438) mts <- rbind(mt_wash,mt_mans) ll_prj <- "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs" mts_sp <- sp::SpatialPoints(sp::coordinates(mts), proj4string = sp::CRS(ll_prj)) get_elev_point(locations = mt_wash, prj = ll_prj) get_elev_point(locations = mt_wash, units="feet", prj = ll_prj) get_elev_point(locations = mt_wash, units="meters", prj = ll_prj) get_elev_point(locations = mts_sp)
请查看
https://cran.r-project.org/web/packages/elevatr/vignettes/introduction_to_elevatr.html
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-22 04:36
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社