Le géocodage dans R avec Google Maps

Avez-vous pensé à utiliser le json appel à la place? En regardant votre code, vous pouvez obtenir la même chose de faire cela (vous aurez besoin d'installer les paquets RCurl et RJSONIO de omegahat.com).

De copier et de les coller dans R:

library(RCurl)
library(RJSONIO)

construct.geocode.url <- function(address, return.call = "json", sensor = "false") {
  root <- "http://maps.google.com/maps/api/geocode/"
  u <- paste(root, return.call, "?address=", address, "&sensor=", sensor, sep = "")
  return(URLencode(u))
}

gGeoCode <- function(address,verbose=FALSE) {
  if(verbose) cat(address,"\n")
  u <- construct.geocode.url(address)
  doc <- getURL(u)
  x <- fromJSON(doc,simplify = FALSE)
  if(x$status=="OK") {
    lat <- x$results[[1]]$geometry$location$lat
    lng <- x$results[[1]]$geometry$location$lng
    return(c(lat, lng))
  } else {
    return(c(NA,NA))
  }
}

Voici comment utiliser les fonctions ci-dessus:

x <- gGeoCode("Philadelphia, PA")

C'est le résultat que vous obtenez. Je pense que dans le code d'origine, lat et de gnl sont mélangés? Mais j'espère que c'est ce que vous voulez:

> x
[1]  39.95233 -75.16379

Espère que ça aide un peu mate,

Tony Breyal

Ceci est une autre option pour le géocodage - il peut être plus facile à analyser:

https://webgis.usc.edu/Services/Geocode/Default.aspx

J'avais besoin pour obtenir tous les anciens d'adresses à partir de geocode pas seulement le premier, j'ai donc écrit une petite fonction pour le faire. Il peut être utilisé pour de géocodage et à de géocodage inversé

geocode <- function(address,reverse=FALSE)  {
  require("RJSONIO")
  baseURL <- "http://maps.google.com/maps/api/geocode/json?sensor=false&"

  # This is not necessary, 
  # because the parameter "address" accepts both formatted address and latlng

  conURL <- ifelse(reverse,paste0(baseURL,'latlng=',URLencode(address)),
                                  paste0(baseURL,'address=',URLencode(address)))  
  con <- url(conURL)  
  data.json <- fromJSON(paste(readLines(con), collapse=""))
  close(con) 
  status <- data.json["status"]
 if(toupper(status) == "OK"){
  t(sapply(data.json$results,function(x) {
      list(address=x$formatted_address,lat=x$geometry$location[1],
                                                 lng=x$geometry$location[2])}))
 } else { 
   warning(status)
   NULL 
 }
}

De géocodage exemple:

de géocodage("Dupont Cir NW, Washington, DC 20036, états-unis")

     address                                                               lat      lng      
[1,] "Dupont Circle Northwest, Washington, DC 20036, USA"                  38.90914 -77.04366
[2,] "Dupont Circle, 1 Dupont Circle Northwest, Washington, DC 20036, USA" 38.90921 -77.04438
[3,] "Washington, DC 20036, USA"                                           38.90808 -77.04061
[4,] "Dupont Circle, Washington, DC 20036, USA"                            38.90958 -77.04344

De Géocodage inversé exemple:

noter que l'adresse peut être mis en forme d'adresse ou de latlng, le reverse paramètre n'est pas utilisé, mais il est, y compris pour une utilisation future avec d'autres services de géocodage

de géocodage("38.910262, -77.043565")

     address                                                    lat      lng      
[1,] "40-58 Dupont Circle Northwest, Washington, DC 20036, USA" 38.91027 -77.04357
[2,] "Washington, DC 20036, USA"                                38.90808 -77.04061
[3,] "Dupont Circle, Washington, DC, USA"                       38.90969 -77.04334
[4,] "Northwest Washington, Washington, DC, USA"                38.94068 -77.06796
[5,] "District of Columbia, USA"                                38.90598 -77.03342
[6,] "Washington, DC, USA"                                      38.90723 -77.03646
[7,] "United States"                                            37.09024 -95.71289