在R中进行反向地理编码的循环

Question

在R中进行反向地理编码的循环

11

我正在尝试对一个大型数据集（约10万条）进行反向地理编码。我使用了来自 ggmap 软件包的 revgeocode 函数。我已经获得了1个条目的结果。

48 Grand View Terrace, San Francisco, 
CA 94114, USA            
48 Grand View Terrace Eureka Valley San Francisco        
San Francisco County                  California United States
postal_code postal_code_suffix

但我需要自动化这个过程，并将其用于整个数据集。

我尝试过

r <- lapply(revgeocode(location = (c(z$lon),c(z$lat)),
             output = "more",
            messaging = FALSE, sensor = FALSE, override_limit = FALSE,
            client = "", signature = ""))

并且在每一步中都遇到了意外的','错误。

我也尝试编写了以下循环：

r <- for(i in 1:10){
  revgeocode(location = ("z$lon", "z$lat"),output = "more", messaging =      FALSE, sensor = FALSE, override_limit = FALSE,client = "", signature = "")}

并且收到了类似的错误。

请提供一些材料或有用的链接，帮助我编写反向地理编码循环。如何验证数据的真实性？

- marine8115

1

Google Maps API 每天的查询限制为 2500 次。因此，您可能希望采取一些措施来考虑这一点。 - Rick Arko

@RickArko 是的，我知道这个事实！有没有办法使用这个 ggmap 代码调用API？我想每1000个查询需要支付0.5美元。 - marine8115

2个回答

0

我曾经遇到过一个类似的问题，就是集成API密钥的问题。基本上，这是一个将API密钥集成到R调用的URL中的问题。如果这个不能帮助你，你需要改变核心代码（在Github上查找），以允许调用密钥的参数。

- pdx

getGeoData <- function(latlng, api_key){ geo_data <- getURL(paste("https://maps.googleapis.com/maps/api/geocode/json?","latlng=",latlng,"&key=",sep="")) geo_data <- fromJSON(geo_data) return(geo_data$results[[1]])}根据帖子的建议，我编写了一个反向地理编码函数。我还尝试应用以下循环，以便可以处理多个查询 for (i in 1:10) { geo_data[i] = getGeoData(unique(y1[i,4]))} 。我收到错误消息number of items to replace is not a multiple of replacement length。有什么建议来解决这个问题吗？ - marine8115

我是一名R语言的新手，但应该写成geo_data[[i]]才对吧？我现在工作很忙，但我很快会去看一下。 - pdx

网页内容由stack overflow 提供, 点击上面的

可以查看英文原文，
原文链接

- Steven Beaupré · Accepted Answer

根据这个答案, 你可以在你的数据框中创建新变量。

我们使用 mapply() 处理你的坐标并将结果返回到列表 res 中。

res <- mapply(FUN = function(lon, lat) { 
  revgeocode(c(lon, lat), output = "more") 
  }, 
  df$lon, df$lat
  )

然后，我们使用data.table中的rbindlist()将列表转换为data.frame（使用fill = TRUE，因为res的所有元素长度不相同，即某些结果没有返回street_number和postal_code），并将其与原始数据cbind()。

cbind(df, data.table::rbindlist(res, fill = TRUE))

更新

回复您的评论，如果您想处理超过2500个查询，您可以订阅Google Maps API高级版以解锁更高的配额。然后，您可以使用signature和client参数将您的凭据传递给revgeocode()。

如文档所述：

购买Google Maps API高级版许可证后，您将收到一封来自Google的欢迎电子邮件，其中包含您的客户端ID和私有加密密钥。

您的客户端ID用于访问Google Maps API高级版的特殊功能。所有客户端ID都以gme-前缀开头。将您的客户端ID作为客户端参数的值传递。使用您的私有加密密钥生成唯一的数字签名。将此签名作为签名参数的值传递。

你可以通过检查revgeocode()源代码来了解它是如何在后台运作的，同时可以查看URL的构造方式：

sensor4url <- paste("&sensor=", sensor, sep = "")
client4url <- paste("&client=", client, sep = "")
signature4url <- paste("&signature=", signature, sep = "")
url_string <- paste("http://maps.googleapis.com/maps/api/geocode/json?latlng=", 
        loc4url, sensor4url, sep = "")
    if (userType == "business") {
        url_string <- paste(url_string, client4url, signature4url, 
            sep = "")
    }

数据

df <- structure(list(lat = c(32.31, 32.19, 34.75, 35.09, 35.35, 34.74 ), lon = 
c(119.827, 119.637, 119.381, 119.364, 119.534, 119.421 )), .Names = 
c("lat", "lon"), row.names = c(21L, 32L, 37L, 48L, 50L, 89L), class = "data.frame")