|
|
@@ -25,21 +25,25 @@ names(racine) <- racines$codlib |
|
|
|
racine <- racine[racine %in% GHM$GHM] |
|
|
|
rm(racines) |
|
|
|
|
|
|
|
# Lecture codes postaux <-> codes insee et garder ceux existants ---- |
|
|
|
read_csv2("insee.csv", col_types = cols(INSEE = col_character())) %>% |
|
|
|
select(Codepos, insee = INSEE) %>% |
|
|
|
mutate(insee = insee %>% str_pad(5, "left", "0")) %>% |
|
|
|
filter(Codepos %in% GHM$Codepos) -> insee |
|
|
|
# # Lecture codes postaux <-> codes insee et sauvegarde carte pour simplification avec mapshaper |
|
|
|
# read_csv2("insee.csv", col_types = cols(INSEE = col_character())) %>% |
|
|
|
# select(Codepos, insee = INSEE) %>% |
|
|
|
# mutate(insee = insee %>% str_pad(5, "left", "0")) -> insee |
|
|
|
# # Lecture carte et join avec codes postaux ---- |
|
|
|
# readOGR("communes", "communes-20150101-100m") -> communes |
|
|
|
# communes@data %>% |
|
|
|
# left_join(insee) -> communes@data |
|
|
|
# rm(insee) |
|
|
|
# communes <- communes[!is.na(communes$Codepos),] |
|
|
|
# writeOGR(communes, "communes_byCP", "communes_byCP", driver = "ESRI Shapefile") |
|
|
|
|
|
|
|
# Lecture carte ---- |
|
|
|
readOGR("communes_byCP", "communes_byCP") -> communes |
|
|
|
|
|
|
|
# Filtre Codes existants dans la base |
|
|
|
communes <- communes[communes$Codepos %in% GHM$Codepos,] |
|
|
|
communes@data$Codepos <- communes@data$Codepos %>% as.character |
|
|
|
|
|
|
|
# Lecture carte et join avec codes postaux ---- |
|
|
|
readOGR("communes", "communes-20150101-simpl") -> communes |
|
|
|
|
|
|
|
communes@data %>% |
|
|
|
left_join(insee) -> communes@data |
|
|
|
|
|
|
|
communes <- communes[!is.na(communes$Codepos),] |
|
|
|
|
|
|
|
rm(insee) |
|
|
|
|
|
|
|
# long/lat de chaque commune ---- |
|
|
|
|
|
|
|