Обратное геокодирование для извлечения компонентов адреса

Я пытаюсь изменить геокодирование с помощью R. Сначала я использовал ggmap, но не смог заставить его работать с моим ключом API. Теперь я пытаюсь это с GoogleWay.

newframe[,c("Front.lat","Front.long")]

  Front.lat Front.long
1 -37.82681   144.9592
2 -37.82681   145.9592

newframe$address <- apply(newframe, 1, function(x){
  google_reverse_geocode(location = as.numeric(c(x["Front.lat"], 
x["Front.long"])),
                         key = "xxxx")
})

Это извлекает переменные в виде списка, но я не могу понять структуру.

Я изо всех сил пытаюсь выяснить, как извлечь компоненты адреса, перечисленные ниже, как переменные в newframe

postal_code, administrative_area_level_1, administrative_area_level_2, locality, route, street_number

Я бы предпочел каждый компонент адреса как отдельную переменную.

Источник

3 ответа

Решение

После обратного геокодирования в newframe$address компоненты адреса могут быть извлечены следующим образом:

# Make a boolean array of the valid ("OK" status) responses (other statuses may be "NO_RESULTS", "REQUEST_DENIED" etc).
sel <- sapply(c(1: nrow(newframe)), function(x){
  newframe$address[[x]]$status == 'OK'
})

# Get the address_components of the first result (i.e. best match) returned per geocoded coordinate.
address.components <- sapply(c(1: nrow(newframe[sel,])), function(x){
  newframe$address[[x]]$results[1,]$address_components
})

# Get all possible component types.
all.types <- unique(unlist(sapply(c(1: length(address.components)), function(x){
  unlist(lapply(address.components[[x]]$types, function(l) l[[1]]))
})))

# Get "long_name" values of the address_components for each type present (the other option is "short_name").
all.values <- lapply(c(1: length(address.components)), function(x){
  types <- unlist(lapply(address.components[[x]]$types, function(l) l[[1]]))
  matches <- match(all.types, types)
  values <- address.components[[x]]$long_name[matches]
})

# Bind results into a dataframe.
all.values <- do.call("rbind", all.values)
all.values <- as.data.frame(all.values)
names(all.values) <- all.types

# Add columns and update original data frame.
newframe[, all.types] <- NA
newframe[sel,][, all.types] <- all.values

Обратите внимание, что я сохранил только первый тип, заданный для каждого компонента, фактически пропустив "политический" тип, поскольку он появляется в нескольких компонентах и, вероятно, является излишним, например, "administrator_area_level_1, political".

Источник

API Google возвращает ответ в формате JSON. Который при переводе на R естественно формирует вложенные списки. Внутренне в googleway это делается через jsonlite::fromJSON()

В googleway Я дал вам возможность вернуть необработанный JSON или список, используя simplify аргумент.

Я намеренно вернул ВСЕ данные из ответа Google и оставил их на усмотрение пользователя, чтобы извлечь интересующие их элементы с помощью обычных операций по подстановке списков.

Сказав все это, в версии разработки googleway я написал несколько функций, помогающих получить доступ к элементам различных вызовов API. Вот три из них, которые могут быть полезны для вас

## Install the development version
# devtools::install_github("SymbolixAU/googleway")

res <- google_reverse_geocode(
  location = c(df[1, 'Front.lat'], df[1, 'Front.long']), 
  key = apiKey
  )

geocode_address(res)
# [1] "45 Clarke St, Southbank VIC 3006, Australia"                    
# [2] "Bank Apartments, 275-283 City Rd, Southbank VIC 3006, Australia"
# [3] "Southbank VIC 3006, Australia"                                  
# [4] "Melbourne VIC, Australia"                                       
# [5] "South Wharf VIC 3006, Australia"                                
# [6] "Melbourne, VIC, Australia"                                      
# [7] "CBD & South Melbourne, VIC, Australia"                          
# [8] "Melbourne Metropolitan Area, VIC, Australia"                    
# [9] "Victoria, Australia"                                            
# [10] "Australia"

geocode_address_components(res)
#        long_name short_name                                  types
# 1             45         45                          street_number
# 2  Clarke Street  Clarke St                                  route
# 3      Southbank  Southbank                    locality, political
# 4 Melbourne City  Melbourne administrative_area_level_2, political
# 5       Victoria        VIC administrative_area_level_1, political
# 6      Australia         AU                     country, political
# 7           3006       3006                            postal_code

geocode_type(res)
# [[1]]
# [1] "street_address"
# 
# [[2]]
# [1] "establishment"      "general_contractor" "point_of_interest" 
# 
# [[3]]
# [1] "locality"  "political"
# 
# [[4]]
# [1] "colloquial_area" "locality"        "political"
Источник

Ты можешь использовать ggmap:revgeocode без труда; Смотри ниже:

library(ggmap)
df <- cbind(df,do.call(rbind,
        lapply(1:nrow(df),
          function(i) 
            revgeocode(as.numeric(
              df[i,2:1]), output = "more")      
                [c("administrative_area_level_1","locality","postal_code","address")])))

#output:
df
#   Front.lat Front.long administrative_area_level_1  locality
#   1 -37.82681   144.9592                    Victoria Southbank
#   2 -37.82681   145.9592                    Victoria    Noojee
#     postal_code                                     address
#   1        3006 45 Clarke St, Southbank VIC 3006, Australia
#   2        3833 Cec Dunns Track, Noojee VIC 3833, Australia

Можете добавить "route" а также "street_number" переменным, которые вы хотите извлечь, но, как вы можете видеть, второй адрес не имеет номера улицы, и это приведет к ошибке.

Примечание: вы также можете использовать sub и извлекать информацию из адреса.

Данные:

df <- structure(list(Front.lat = c(-37.82681, -37.82681), Front.long = 
      c(144.9592, 145.9592)), .Names = c("Front.lat", "Front.long"), class = "data.frame", 
      row.names = c(NA, -2L))
Источник
Другие вопросы по тегам