Geocoding example - uses PDOK LocationServer

geocoding_example.R

library(httr)
library(dplyr)
library(readr)
library(purrr)
library(magrittr)

f <- "http://www.landelijkregisterkinderopvang.nl/opendata/export_opendata_lrk.csv"

input <- read_delim(f, delim = ";") %>%
         mutate(adres = paste(opvanglocatie_adres, opvanglocatie_postcode, opvanglocatie_woonplaats)) %>%
         select(lrk_id, adres, bag_id) %>%
         head(100) # Ik pak alleen de eerste 100 rijen in de dataset, anders duurt het zo lang ;)

# Er zit een fout in het inputbestand: voor sommige rijen ontbreken de laatste 1-10 kolommen.

geocode <- function(adres){
  GET(url = "https://geodata.nationaalgeoregister.nl", 
           path = "locatieserver/v3/free", 
           query = list(qf = "exacte_match^1 suggest^0.5 huisnummer^0.5 huisletter^0.5 huisnummertoevoeging^0.5",
                        fq = "type:adres AND bron:BAG",
                        fl = "*",
                        q  = adres,
                        rows = 1)) %>% 
    content() %>% 
    pluck('response', 'docs', 1)
}

output <- map_dfr(input$adres, function(x) {
  attrs <- c("adresseerbaarobject_id", "weergavenaam", "centroide_ll")
  geocode(x) %>% extract(attrs) 
})

matched <- bind_cols(input, output)
print(matched)