如何从R中的城市名称和国家/地区获取经度和纬度坐标？

Question

我有一长串的城市名称和国家/地区,我想在地图上绘制它们.为了做到这一点,我需要每个城市的经度和纬度信息.

我的表被调用test并具有以下结构:

Cityname  CountryCode
New York  US
Hamburg   DE
Amsterdam NL

Answer 1

使用以下代码我已成功解决了问题.

library(RJSONIO)
nrow <- nrow(test)
counter <- 1
test$lon[counter] <- 0
test$lat[counter] <- 0
while (counter <= nrow){
  CityName <- gsub(' ','%20',test$CityLong[counter]) #remove space for URLs
  CountryCode <- test$Country[counter]
  url <- paste(
    "http://nominatim.openstreetmap.org/search?city="
    , CityName
    , "&countrycodes="
    , CountryCode
    , "&limit=9&format=json"
    , sep="")
  x <- fromJSON(url)
  if(is.vector(x)){
    test$lon[counter] <- x[[1]]$lon
    test$lat[counter] <- x[[1]]$lat    
  }
  counter <- counter + 1
}

由于这是调用外部服务(openstreetmaps.org),因此对于较大的数据集可能需要一段时间.但是,您可能只在新城市添加到列表中时偶尔执行此操作.

Answer 2

还有其他一些选择.

ggmaps

ggmaps有一个geocode使用谷歌地图进行地理编码的功能.这限制你每天2,500.

taRifx.geo

taRifx.geo的最新版本具有geocode使用Google或Bing Maps进行地理编码的功能.Bing版本要求您使用(免费)Bing帐户,但作为回报,您可以对更多条目进行地理编码.此版本的功能:

• 服务选择(支持Bing和Google Maps)
• 登录支持(特别是对于Bing,它需要一个帐户密钥,但在交换中允许每天更多的请求)
• 一次对整个data.frame进行地理编码,包括一些节省时间的操作,例如忽略已经进行过地理编码的任何行
• 强大的批量地理编码(这样任何错误都不会导致整个data.frame的地理编码丢失,对于更大的工作)
• 路线寻找(从A点到B点的旅行时间)

Answer 3

试试这个我认为这将是解决这个问题的更好方法

> library(ggmap) 
Loading required package: ggplot2
Google Maps API Terms of Service: http://developers.google.com/maps/terms.
Please cite ggmap if you use it: see citation('ggmap') for details.

#Now you can give city name or country name individually

> geocode("hamburg")
Information from URL : http://maps.googleapis.com/maps/api/geocode/json?address=hamburg&sensor=false
       lon      lat
1 9.993682 53.55108

geocode("amsterdam")
Information from URL : http://maps.googleapis.com/maps/api/geocode/json?address=amsterdam&sensor=false
       lon      lat
1 4.895168 52.37022

> geocode("new york")
Information from URL : http://maps.googleapis.com/maps/api/geocode/json?address=new+york&sensor=false
        lon      lat
1 -74.00594 40.71278