博文
获取某一点的海拔高程的R函数
|
刚审稿发现有个可以获取一个北美经纬度点的海拔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
https://blog.sciencenet.cn/blog-267448-1243981.html
上一篇：《数量生态学-R语言的应用》中文版第二版已经出版
下一篇：R4.0的版本下定义泛函数的注意事项
