Skip to content

Instantly share code, notes, and snippets.

@datagistips
Last active January 25, 2023 16:07
Show Gist options
  • Save datagistips/8f55425396e93af6df5a897df70a92c3 to your computer and use it in GitHub Desktop.
Save datagistips/8f55425396e93af6df5a897df70a92c3 to your computer and use it in GitHub Desktop.
Récupérer les mairies de France avec geo.api.gouv.fr
# Auparavant, on trouvait facilement les emplacements des mairies depuis Le Répertoire Géographique des Communes
# https://geo.data.gouv.fr/fr/datasets/2d67a8d0da4ba99004a88cbe93fc0e0fceb43bb8
# La donnée n'est plus disponible
# Mais l'API geo.api.gouv.fr permet de retrouver l'emplacement des mairies.
# Cependant, il y a des soucis avec la requête https://geo.api.gouv.fr/communes?geometry=mairie&format=geojson
# qui retourne 'Bad request'
# On utilise ici une itération le long des différents départements
# On exporte le résultat en GéoPackage
library(sf)
library(tidyverse)
library(jsonlite)
library(glue)
# Fonction pour récupérer les mairies d'un département
get_mairies_dep <- function(dep) {
# Read API
url <- glue("https://geo.api.gouv.fr/departements/{dep}/communes?geometry=mairie&format=geojson")
res <- read_json(url)
# Get communes
feats <- res$features
out <- vector(mode = "list")
for(i in 1:length(feats)) {
feat <- res$features[[i]]
# Informations attributaires
nom <- feat$properties$nom
code <- feat$properties$code
siren <- feat$properties$siren
codeRegion <- feat$properties$codeRegion
population <- feat$properties$population
# Géométrie
coords <- feat$geometry$coordinates
long <- coords[[1]]
lat <- coords[[2]]
geom <- st_point(c(long, lat)) %>% st_sfc
# Data frame
df <- data.frame(nom, code, siren, codeRegion, population)
st_geometry(df) <- geom
out[[i]] <- df
}
res <- do.call(rbind, out) %>% st_set_crs(4326)
return(res)
}
# Fonction pour récupérer les mairies de plusieurs départements
get_mairies <- function(deps) {
out <- vector(mode = "list")
for(i in 1:length(deps)) {
dep <- deps[i]
print(dep)
out[[i]] <- get_mairies_dep(dep)
}
res <- do.call(rbind, out)
return(res)
}
# Récupérer les mairies des Bouches-du-Rhône
mairies_d13 <- get_mairies_dep("13")
mairies_d13
# Récupérer les codes départements nationaux
url <- "https://geo.api.gouv.fr/departements"
j <- read_json(url)
deps <- sapply(j, function(x) x$code)
# On récupère les mairies sur les départements concernés
mairies_fr <- get_mairies(deps)
# On exporte la donnée en GPKG
st_write(mairies_fr, "mairies_fr.gpkg", delete_dsn = T)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment