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

时间:2012-12-16 20:16:02

标签: r maps openstreetmap rjsonio

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

我的表名为test,具有以下结构:

Cityname  CountryCode
New York  US
Hamburg   DE
Amsterdam NL

4 个答案:

答案 0 :(得分:15)

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

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),因此对于较大的数据集可能需要一段时间。但是,您可能只在新城市添加到列表中时偶尔执行此操作。

答案 1 :(得分:11)

为您提供其他一些选择。

<强> ggmaps

ggmaps有一个功能geocode,它使用Google地图进行地理编码。这限制你每天2,500。

<强> taRifx.geo

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

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

答案 2 :(得分:6)

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

&#13;
&#13;
> 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
&#13;
&#13;
&#13;

答案 3 :(得分:-1)

试试这个...

function geocodeAddress(geocoder, resultsMap) {
    var address = document.getElementById('address').value;
    geocoder.geocode({'address': address}, function(results, status) {
      if (status === 'OK') {
        resultsMap.setCenter(results[0].geometry.location);
        var marker = new google.maps.Marker({
          map: resultsMap,
          position: results[0].geometry.location,
          draggable:true
        });
        var infowindow = new google.maps.InfoWindow({
            content: "Please drag this marker to your position.."
          });
          infowindow.open(resultsMap,marker);
         document.getElementById('lat').value=marker.getPosition().lat();
        document.getElementById('lng').value=marker.getPosition().lng();
         marker.addListener('drag', handleEvent);
marker.addListener('dragend', handleEvent);
      } else {
        alert('Geocode was not successful for the following reason: ' + status);
      }
    });
  }

get full code from here..