#编者信息

熊荣川

六盘水师范学院生物信息学实验室

xiongrongchuan@126.com

http://blog.sciencenet.cn/u/Bearjazz

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)

}

#使用方法

数据的输入格式如下图，输入为表格（或矩阵）