|||
#R语言根据经纬度获得海拔数据
#编者信息 熊荣川 六盘水师范学院生物信息学实验室 | ||||
Google earth上可以很容易的获得经纬度信息,但是海拔高度却不能直接提供,但据说其数据库中是有相应的海拔高度的数据的,下面是R语言爱好者使用R语言根据经纬度查询google数据库从而获取经纬度的代码: (原文地址 http://stackoverflow.com/questions/21593868/extracting-elevation-from-website-for-lat-lon-points-in-australia-using-r) #预装函数 googEl <- function(locs) { require(RJSONIO) locstring <- paste(do.call(paste, list(locs[, 2], locs[, 1], sep=',')), collapse='|') u <- sprintf('http://maps.googleapis.com/maps/api/elevation/json?locations=%s&sensor=false', locstring) res <- fromJSON(u) out <- t(sapply(res[[1]], function(x) { c(x[['location']]['lat'], x[['location']]['lng'], x['elevation'], x['resolution']) })) rownames(out) <- rownames(locs) return(out) }
| ||||
#使用方法 数据的输入格式如下图,输入为表格(或矩阵)
|
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-20 01:41
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社