class: center, middle, inverse, title-slide # Données géographiques et cartographie sous R ## Manipulation de données communales en historique et réalisation de cartes ### Kim Antunez ### 06/03/2018 --- background-image: url("img/giphy1.gif") background-size: 100% class: inverse, center, middle # 0. Introduction sur les géographies communales --- # Code Officiel Géographique (COG) Le **Code officiel géographique** (COG) est un document de référence publié par l’**Insee** qui rassemble les **codifications (numérotations et libellés) des communes**, des cantons, des arrondissements, des départements, des régions... Ce code est **révisé annuellement**, en fonction notamment des **fusions et associations** de communes ou de territoires, des **défusions** et des **changements de dénomination**. --- # Types de modifications communales (1/7)
--- # Types de modifications communales (2/7)
--- # Types de modifications communales (3/7) - les populations d'une même commune varient d'une année sur l'autre - la liste des communes change également chaque année <center><img src="img/ex_COG.png" height="320px"></center> Comment fait-on pour intégrer des données de plusieurs millésimes dans un même data.frame ? → **Je dois savoir transformer les pop2011 en COG2012 ou les pop2012 en COG2011** --- # Types de modifications communales (4/7) - Variables numériques et fusions = **FACILE** - *Question : Quelle situation pour la population (2011) en COG 2012 ?* <center><img src="img/varnum_fusion1.png" height="400px"></center> --- # Types de modifications communales (4/7) - Variables numériques et fusions = **FACILE** - *Question : Quelle situation pour la population (2011) en COG 2012 ?* <center><img src="img/varnum_fusion2.png" height="400px"></center> --- # Types de modifications communales (5/7) - Variables numériques et défusions = **DIFFICILE** - *Question : Quelle situation pour la population (2013) en COG 2014 ?* <center><img src="img/varnum_defusion1.png" height="400px"></center> --- # Types de modifications communales (5/7) - Variables numériques et défusions = **DIFFICILE** - *Question : Quelle situation pour la population (2013) en COG 2014 ?* <center><img src="img/varnum_defusion2.png" height="400px"></center> --- # Types de modifications communales (6/7) - Typologie et défusions = **FACILE** - *Question : Quelles typologies pour les communes défusionnées en COG 2014 ?* <center><img src="img/typo_defusion1.png" height="300px"></center> --- # Types de modifications communales (6/7) - Typologie et défusions = **FACILE** - *Question : Quelles typologies pour les communes défusionnées en COG 2014 ?* <center><img src="img/typo_defusion2.png" height="300px"></center> --- # Types de modifications communales (7/7) - Typologie et fusions = **DIFFICILE** - *Question : Quelle typologie pour la commune fusionnée en COG 2012 ?* <center><img src="img/typo_fusion1.png" height="300px"></center> --- # Types de modifications communales (7/7) - Typologie et fusions = **DIFFICILE** - *Question : Quelle typologie pour la commune fusionnée en COG 2012 ?* <center><img src="img/typo_fusion2.png" height="300px"></center> --- # Les différents évènements depuis 1968
--- # Le package COGugaison (1/2) - détecter le millésime du code officiel géographique d'une table de données communales : fonction **`COG_akinator`** et apparier une liste de communes d'une base de données avec celle du COG d'un millésime donné avec `apparier_COG` - visualiser les modifications communales (fusions, défusions, changements de codes ou de noms) qui ont eu lieu entre deux dates : `modifications_communales` - visualiser sur un graphique interactif la trajectoire d'une commune donnée, c'est-à-dire toutes les modifications qui ont eu lieu depuis 1968 : **`trajectoire_commune`** et `trajectoire_commune_shiny` - transformer des tables de données numériques en géographie au premier janvier d'une année souhaitée : **`changement_COG_varNum`** --- # Le package COGugaison (2/2) - transformer des typologies de communes en géographie au premier janvier d'une année souhaitée en ayant le choix entre plusieurs hypothèses de classement en cas de fusion de communes de classes différentes (attribuer une unique classe à toutes les communes fusionnées, attribuer la classe qui contient le plus de population, définir une classe absorbante, une classe absorbée ou une classe spécifique aux regroupements de plusieurs communes de classes différentes) : **`changement_COG_typo`**. Il est également possible d'isoler dans une table les communes fusionnées appartenant à des classes différentes : `changement_COG_typo_details` - permettre d'agréger les tables de données communales à de nombreux échelons supra-communaux administratifs (EPCI, arrondissements, cantons-villes, départements, régions) ou d'étude (bassins de vie, zones d'emploi, unités urbaines, aires urbaines) : **`nivsupra`** - gérer des cas particuliers comme les codes Insee des communes corses (`modification_Corse`) ou des arrondissements municipaux de Paris, Lyon, et Marseille (**`enlever_PLM`**) ou encore l'ancienne commune de l'Oudon (`modification_Oudon`) --- # Les données utiles dans COGugaison (1/4) > - COG ```r head(COG2017) ``` ``` ## CODGEO LIBGEO POP ## 1 01001 L'Abergement-Clémenciat 767 ## 2 01002 L'Abergement-de-Varey 239 ## 3 01004 Ambérieu-en-Bugey 14022 ## 4 01005 Ambérieux-en-Dombes 1627 ## 5 01006 Ambléon 109 ## 6 01007 Ambronay 2570 ``` --- # Les données utiles dans COGugaison (2/4) > - TABLES_PASSAGE ```r head(PASSAGE_2015_2016_insee) ``` ``` ## cod2015 cod2016 annee typemodif ratio ## 1 14697 14472 2016-01-01 c 1 ## 2 01015 01015 2016-01-01 f 1 ## 3 01340 01015 2016-01-01 f 1 ## 4 01080 01080 2016-01-01 f 1 ## 5 01119 01080 2016-01-01 f 1 ## 6 01176 01187 2016-01-01 f 1 ``` --- # Les données utiles dans COGugaison (3/4) ```r head(table_supracom_2017) ``` ``` ## CODGEO LIBGEO DEP REG EPCI NATURE_EPCI ARR CV ## 1 01001 L'Abergement-Clémenciat 01 84 200069193 CC 012 0108 ## 2 01002 L'Abergement-de-Varey 01 84 240100883 CC 011 0101 ## 3 01004 Ambérieu-en-Bugey 01 84 240100883 CC 011 0101 ## 4 01005 Ambérieux-en-Dombes 01 84 200042497 CC 012 0122 ## 5 01006 Ambléon 01 84 200040350 CC 011 0104 ## 6 01007 Ambronay 01 84 240100883 CC 011 0101 ## ZE2010 UU2010 TUU2014 TDUU2014 AU2010 TAU2014 CATAEU2010 BV2012 ## 1 8213 01000 0 05 997 00 120 01093 ## 2 8201 01000 0 04 002 09 112 01004 ## 3 8201 01302 3 32 002 09 112 01004 ## 4 8213 01000 0 06 002 09 112 69123 ## 5 8216 01000 0 03 998 00 300 01034 ## 6 8201 01000 0 07 002 09 112 01004 ``` --- # Les données utiles dans COGugaison (4/4) ```r head(libelles_supracom_2017) ``` ``` ## NIVGEO CODGEO LIBGEO NB_COM ## 1 ARR 011 Belley 115 ## 2 ARR 012 Bourg-en-Bresse 201 ## 3 ARR 013 Gex 27 ## 4 ARR 014 Nantua 65 ## 5 ARR 021 Château-Thierry 108 ## 6 ARR 022 Laon 244 ``` --- # Exemple de fonction COGugaison (1/2) <center><img src="img/typo_fusion2.png" height="120px"></center> **changement_COG_typo** : Transformer des typologies - `methode_classe_fusion` : pour toutes les communes qui ont fusionné entre 2 dates indiquer comme classe la valeur inscrite dans "mot_fusion" y compris pour les fusions de communes de mêmes classes - `methode_difference` : créer une classe spécifique dont le nom est contenu dans mot_difference (par défaut) - `methode_max_pop` : attribuer la classe contenant le maximum de population des communes fusionnées - `methode_classe_absorbante` : attribuer la classe dite absorbante à toute commune fusionnée contenant au moins une ancienne commune appartenant à cette classe absorbante - `methode_classe_absorbee` : ne pas tenir compte de cette classe dite "absorbée" pour toute commune fusionnée contenant au moins une ancienne commune appartenant à cette classe absorbée --- # Exemple de fonction COGugaison (2/2) ```r table_modif_com <- cbind(COG1968, COG1968[, 1])[, c(1,4)] %>% setNames(c("CODGEO", "typo")) table_modif_com[1, ] ``` ``` ## CODGEO typo ## 1 01001 01001 ``` ```r table_modif_com <- changement_COG_typo(table_modif_com, methode_fusion = "methode_classe_fusion", mot_fusion = "communes fusionnées", annees = 1968:2017, codgeo_entree = "CODGEO", typos = "typo", donnees_insee = F) table_modif_com$typo[table_modif_com$typo %in% ( table_modif_com$typo[duplicated(table_modif_com$typo) & table_modif_com$typo != "communes fusionnées"])] <- "communes defusionnées" table_modif_com$typo[which(!table_modif_com$typo %in% c("communes fusionnées", "communes defusionnées"))] <- "communes inchangées" table_modif_com[1, ] ``` ``` ## CODGEO typo ## 1 01001 communes inchangées ``` --- # Modifications communales depuis 1968 <img src="180306_RLadies_COGugaison_carto_files/figure-html/unnamed-chunk-12-1.svg" style="display: block; margin: auto;" /> --- background-image: url("img/fond_carte.png") background-size: 100% # Démarche **Objectif** : Réaliser des cartes du nombre de femmes et nombre de femmes pour 100 hommes de 1968 à aujourd'hui - Récupérer les données sur insee.fr et observer à quoi elles ressemblent - Transformer les données (changements de COG + agrégation à une échelle supra-communale adaptée) - Réaliser les cartes --- background-image: url("img/giphy2.gif") background-size: 100% class: inverse, center, middle # 1. Récupération des données --- # Téléchargeons les données... <center><a href="https://www.insee.fr/fr/statistiques/1893204" target="_self"><img src="img/screenshot_insee.jpg" height="450px"></a></center> --- # Observons à quoi elles ressemblent... <center><img src="img/screenshot_fichier.jpg"></center> --- # Chargement des tables com1968...com2014 ```r library(readxl) recup_bdd <- function(annee) { bdd <- read_excel("data/pop-sexe-age-quinquennal6814.xls", sheet = paste0("COM_", annee), skip = 13, col_names = TRUE, na = "") %>% mutate(codgeo = paste0(substr(DR, 1, 2), CR)) %>% filter(!is.na(get(paste0("ageq_rec01s1rpop", annee)))) %>% mutate(total_s1 = rowSums(select(., contains("s1"))), total_s2 = rowSums(select(., contains("s2")))) %>% select(codgeo, total_s1, total_s2) return(bdd) } bdd_1975 <- recup_bdd(1975) bdd_1999 <- recup_bdd(1999) bdd_2014 <- recup_bdd(2014) c(nrow(bdd_1975), nrow(bdd_1999), nrow(bdd_2014)) ``` ``` ## [1] 36413 36696 35906 ``` --- background-image: url("img/giphy3.gif") background-size: 100% class: inverse, center, middle # 2. Transformation des données --- # En quel COG sont les données ? ```r COG_akinator(bdd_1975$codgeo, donnees_insee = TRUE) ``` ``` # [1] "COG1975" ``` ```r COG_akinator(bdd_1999$codgeo, donnees_insee = TRUE) ``` ``` # [1] "COG1999" ``` ```r COG_akinator(bdd_2014$codgeo, donnees_insee = TRUE) ``` ``` # [1] "COG2016" ``` ```r # REMARQUE :les recensements récents sont en COG n+2 ``` --- # Changements de COG ```r # Transformer les tables en bons COG # idem pour bdd_1999 et bdd_2014 bdd_1975_COG2016 <- bdd_1975 %>% enlever_PLM() %>% changement_COG_varNum(1975:2016, donnees_insee = TRUE) ``` ```r # Fusionner les tables bdd <- plyr::join_all(list(COG2016 %>% rename(codgeo = CODGEO) %>% select(-POP), bdd_1975_COG2016, bdd_1999_COG2016, bdd_2014_COG2016), type = "left", by = "codgeo") %>% setNames(nm = c("CODGEO", "LIBGEO", "total_s1_1975", "total_s2_1975", "total_s1_1999", "total_s2_1999", "total_s1_2014", "total_s2_2014")) # vérifier les totaux. Par exemple... sum(bdd_2014$total_s1, na.rm = TRUE) - sum(bdd$total_s1_2014, na.rm = TRUE) ``` ``` # [1] 0 ``` --- # Echelles géo supra-communales (1/2) <img src="180306_RLadies_COGugaison_carto_files/figure-html/unnamed-chunk-19-1.svg" style="display: block; margin: auto;" /> --- # Echelles géo supra-communales (2/2) <img src="180306_RLadies_COGugaison_carto_files/figure-html/unnamed-chunk-20-1.svg" style="display: block; margin: auto;" /> --- # Agrégation par Zone d'Emploi ```r bdd[1, ] ``` ``` # CODGEO LIBGEO total_s1_1975 total_s2_1975 total_s1_1999 # 1 01001 L'Abergement-Clémenciat 210 240 384 # total_s2_1999 total_s1_2014 total_s2_2014 # 1 332 385 380 ``` ```r bdd_ze <- bdd %>% nivsupra(COG = 2016, nivsupra = "ZE2010") %>% mutate(rapport_FH_1975 = 100 * total_s2_1975/total_s1_1975, rapport_FH_1999 = 100 * total_s2_1999/total_s1_1999, rapport_FH_2014 = 100 * total_s2_2014/total_s1_2014) ``` --- background-image: url("img/giphy4.gif") background-size: 100% class: inverse, center, middle # 3. Cartographie --- # Format des couches cartographiques (1/3) **Avant = spdf (SpatialPolygonsDataFrame)** - {rgdal} : import/export des données et les opérations de projections - {sp} : gestion et l’affichage des données : *coordinates*, *SpatialPolygonsDataFrame* - {rgeos} : géotraitements : *gUnion*, *gSimplify*, *gBuffer* - format : polygons*data.frame imbriqués dans des listes **Maintenant = sf** - un seul package : sf (nov. 2016) - auteur (Edzer Pebesma) = contributeur / maintainer des packages précédents - format : data.frame enrichi d’une « colonne-liste » de géométries - avantage : format considéré plus pratique « tidy » et compatible avec %>% --- # Format des couches cartographiques (2/3) ```r library(rgdal) ``` ```r comm_2016 <- readOGR( dsn = "data/couches_geo", layer = "COM_2016_CARTElette", stringsAsFactors = FALSE, encoding = "UTF-8", verbose = FALSE, ) ``` ```r str(comm_2016) ``` ``` # Formal class 'SpatialPolygonsDataFrame' [package "sp"] with 5 slots # ..@ data :'data.frame': 35885 obs. of 4 variables: # .. ..$ insee : chr [1:35885] "01001" "01002" "01004" "01005" ... # .. ..$ nom : chr [1:35885] "L'Abergement-Clémenciat" "L'Abergement-de-Varey" "Ambérieu-en-Bugey" "Ambérieux-en-Dombes" ... # .. ..$ wikipedia: chr [1:35885] "fr:L'Abergement-Clémenciat" "fr:L'Abergement-de-Varey" "fr:Ambérieu-en-Bugey" "fr:Ambérieux-en-Dombes" ... # .. ..$ surf_ha : num [1:35885] 1567 909 2449 1602 597 ... # ..@ polygons :List of 35885 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 848615 6563353 # .. .. .. .. .. .. ..@ area : num 15703242 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:44, 1:2] 846839 847056 847011 847281 847063 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 848615 6563353 # .. .. .. ..@ ID : chr "0" # .. .. .. ..@ area : num 15703242 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 887861 6548427 # .. .. .. .. .. .. ..@ area : num 9261360 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:19, 1:2] 886003 886699 886996 887375 887769 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 887861 6548427 # .. .. .. ..@ ID : chr "1" # .. .. .. ..@ area : num 9261360 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 883678 6542925 # .. .. .. .. .. .. ..@ area : num 24363731 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:44, 1:2] 880225 880443 880936 881747 882368 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 883678 6542925 # .. .. .. ..@ ID : chr "2" # .. .. .. ..@ area : num 24363731 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 847960 6545832 # .. .. .. .. .. .. ..@ area : num 1.6e+07 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:50, 1:2] 845172 845410 845616 845559 845817 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 847960 6545832 # .. .. .. ..@ ID : chr "3" # .. .. .. ..@ area : num 1.6e+07 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 901675 6520012 # .. .. .. .. .. .. ..@ area : num 6028013 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:21, 1:2] 899808 900067 900467 900567 900777 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 901675 6520012 # .. .. .. ..@ ID : chr "4" # .. .. .. ..@ area : num 6028013 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 882413 6547832 # .. .. .. .. .. .. ..@ area : num 33522314 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:49, 1:2] 877923 878436 878295 878722 878761 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 882413 6547832 # .. .. .. ..@ ID : chr "5" # .. .. .. ..@ area : num 33522314 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 880717 6540140 # .. .. .. .. .. .. ..@ area : num 5253294 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:17, 1:2] 878877 880563 881019 882492 882978 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 880717 6540140 # .. .. .. ..@ ID : chr "6" # .. .. .. ..@ area : num 5253294 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 906468 6524295 # .. .. .. .. .. .. ..@ area : num 7252612 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:23, 1:2] 905162 905472 906102 906569 906684 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 906468 6524295 # .. .. .. ..@ ID : chr "7" # .. .. .. ..@ area : num 7252612 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 916619 6538230 # .. .. .. .. .. .. ..@ area : num 29774898 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:18, 1:2] 913977 914015 914427 914499 914276 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 916619 6538230 # .. .. .. ..@ ID : chr "8" # .. .. .. ..@ area : num 29774898 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 904910 6570773 # .. .. .. .. .. .. ..@ area : num 15157770 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:21, 1:2] 902752 903120 903291 904241 904993 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 904910 6570773 # .. .. .. ..@ ID : chr "9" # .. .. .. ..@ area : num 15157770 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 894320 6547718 # .. .. .. .. .. .. ..@ area : num 21867079 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:27, 1:2] 891411 891720 891924 891878 892088 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 894320 6547718 # .. .. .. ..@ ID : chr "10" # .. .. .. ..@ area : num 21867079 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 893767 6535342 # .. .. .. .. .. .. ..@ area : num 14094114 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:35, 1:2] 890996 891883 891738 892029 893050 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 893767 6535342 # .. .. .. ..@ ID : chr "11" # .. .. .. ..@ area : num 14094114 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 907115 6579494 # .. .. .. .. .. .. ..@ area : num 23524342 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:19, 1:2] 904080 904927 905749 905659 905847 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 907115 6579494 # .. .. .. ..@ ID : chr "12" # .. .. .. ..@ area : num 23524342 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 905481 6516034 # .. .. .. .. .. .. ..@ area : num 22725883 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:48, 1:2] 902607 902720 902846 903903 903748 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 905481 6516034 # .. .. .. ..@ ID : chr "13" # .. .. .. ..@ area : num 22725883 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 850365 6599370 # .. .. .. .. .. .. ..@ area : num 17468645 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:37, 1:2] 846835 847047 847631 847915 848041 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 850365 6599370 # .. .. .. ..@ ID : chr "14" # .. .. .. ..@ area : num 17468645 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 892351 6540112 # .. .. .. .. .. .. ..@ area : num 7797419 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:20, 1:2] 890503 890559 890918 890978 892090 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 892351 6540112 # .. .. .. ..@ ID : chr "15" # .. .. .. ..@ area : num 7797419 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 900385 6531564 # .. .. .. .. .. .. ..@ area : num 6955999 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:19, 1:2] 898151 898590 899444 899615 900246 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 900385 6531564 # .. .. .. ..@ ID : chr "16" # .. .. .. ..@ area : num 6955999 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 841011 6545410 # .. .. .. .. .. .. ..@ area : num 5619680 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:10, 1:2] 839376 839412 839785 841829 842173 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 841011 6545410 # .. .. .. ..@ ID : chr "17" # .. .. .. ..@ area : num 5619680 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 908620 6533558 # .. .. .. .. .. .. ..@ area : num 3271716 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:19, 1:2] 907642 907746 908385 908548 908820 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 908620 6533558 # .. .. .. ..@ ID : chr "18" # .. .. .. ..@ area : num 3271716 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 844739 6588971 # .. .. .. .. .. .. ..@ area : num 4657708 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:16, 1:2] 843409 845031 845657 845838 845560 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 844739 6588971 # .. .. .. ..@ ID : chr "19" # .. .. .. ..@ area : num 4657708 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 867790 6578577 # .. .. .. .. .. .. ..@ area : num 1.9e+07 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:47, 1:2] 865028 865505 865604 866181 866383 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 867790 6578577 # .. .. .. ..@ ID : chr "20" # .. .. .. ..@ area : num 1.9e+07 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 850270 6582255 # .. .. .. .. .. .. ..@ area : num 4e+07 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:79, 1:2] 846604 846839 847141 847724 847790 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 850270 6582255 # .. .. .. ..@ ID : chr "21" # .. .. .. ..@ area : num 4e+07 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 848490 6580543 # .. .. .. .. .. .. ..@ area : num 818057 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:9, 1:2] 847822 848045 848262 849012 849280 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 848490 6580543 # .. .. .. ..@ ID : chr "22" # .. .. .. ..@ area : num 818057 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 863631 6527896 # .. .. .. .. .. .. ..@ area : num 1.8e+07 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:32, 1:2] 860978 861122 861562 861639 862194 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 863631 6527896 # .. .. .. ..@ ID : chr "23" # .. .. .. ..@ area : num 1.8e+07 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 846871 6558909 # .. .. .. .. .. .. ..@ area : num 9038886 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:19, 1:2] 845203 845984 845984 846850 847023 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 846871 6558909 # .. .. .. ..@ ID : chr "24" # .. .. .. ..@ area : num 9038886 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 873546 6594570 # .. .. .. .. .. .. ..@ area : num 14126842 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:57, 1:2] 871390 871693 871652 872007 873024 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 873546 6594570 # .. .. .. ..@ ID : chr "25" # .. .. .. ..@ area : num 14126842 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 835960 6546058 # .. .. .. .. .. .. ..@ area : num 953698 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:8, 1:2] 835316 835342 835766 836344 836572 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 835960 6546058 # .. .. .. ..@ ID : chr "26" # .. .. .. ..@ area : num 953698 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 903132 6574280 # .. .. .. .. .. .. ..@ area : num 7556646 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:14, 1:2] 901326 901770 903029 902980 903279 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 903132 6574280 # .. .. .. ..@ ID : chr "27" # .. .. .. ..@ area : num 7556646 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 866046 6531102 # .. .. .. .. .. .. ..@ area : num 13278058 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:26, 1:2] 863723 863821 864145 863984 864166 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 866046 6531102 # .. .. .. ..@ ID : chr "28" # .. .. .. ..@ area : num 13278058 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 918754 6560254 # .. .. .. .. .. .. ..@ area : num 15371679 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:52, 1:2] 915464 916178 916196 916062 916269 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 918754 6560254 # .. .. .. ..@ ID : chr "29" # .. .. .. ..@ area : num 15371679 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 908609 6520235 # .. .. .. .. .. .. ..@ area : num 22547714 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:41, 1:2] 906380 906571 906447 906638 906534 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 908609 6520235 # .. .. .. ..@ ID : chr "30" # .. .. .. ..@ area : num 22547714 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 914520 6576390 # .. .. .. .. .. .. ..@ area : num 17554346 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:27, 1:2] 910520 911279 911865 913120 913823 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 914520 6576390 # .. .. .. ..@ ID : chr "31" # .. .. .. ..@ area : num 17554346 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 905528 6535180 # .. .. .. .. .. .. ..@ area : num 19822895 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:35, 1:2] 903381 903748 903753 903475 903467 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 905528 6535180 # .. .. .. ..@ ID : chr "32" # .. .. .. ..@ area : num 19822895 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 892682 6529541 # .. .. .. .. .. .. ..@ area : num 15421825 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:32, 1:2] 890142 890972 891329 891628 892199 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 892682 6529541 # .. .. .. ..@ ID : chr "33" # .. .. .. ..@ area : num 15421825 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 875604 6582485 # .. .. .. .. .. .. ..@ area : num 18304400 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:38, 1:2] 872761 872773 872984 872965 873136 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 875604 6582485 # .. .. .. ..@ ID : chr "34" # .. .. .. ..@ area : num 18304400 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 913517 6532262 # .. .. .. .. .. .. ..@ area : num 10438617 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:19, 1:2] 911715 912718 912487 912573 913242 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 913517 6532262 # .. .. .. ..@ ID : chr "35" # .. .. .. ..@ area : num 10438617 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 856802 6587333 # .. .. .. .. .. .. ..@ area : num 11184800 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:29, 1:2] 854813 854886 855199 855385 856140 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 856802 6587333 # .. .. .. ..@ ID : chr "36" # .. .. .. ..@ area : num 11184800 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 883335 6540263 # .. .. .. .. .. .. ..@ area : num 3442777 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:30, 1:2] 881332 881522 881794 881881 882401 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 883335 6540263 # .. .. .. ..@ ID : chr "37" # .. .. .. ..@ area : num 3442777 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 842289 6570521 # .. .. .. .. .. .. ..@ area : num 2734913 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:9, 1:2] 841143 842100 842417 843043 842791 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 842289 6570521 # .. .. .. ..@ ID : chr "38" # .. .. .. ..@ area : num 2734913 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 855030 6528617 # .. .. .. .. .. .. ..@ area : num 10564831 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:23, 1:2] 853057 853701 853772 854080 854669 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 855030 6528617 # .. .. .. ..@ ID : chr "39" # .. .. .. ..@ area : num 10564831 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 913718 6557021 # .. .. .. .. .. .. ..@ area : num 13667258 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:27, 1:2] 909719 910030 910455 911056 911567 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 913718 6557021 # .. .. .. ..@ ID : chr "40" # .. .. .. ..@ area : num 13667258 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 857788 6541298 # .. .. .. .. .. .. ..@ area : num 15856696 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:37, 1:2] 855261 855375 855776 855858 856834 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 857788 6541298 # .. .. .. ..@ ID : chr "41" # .. .. .. ..@ area : num 15856696 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 849591 6571460 # .. .. .. .. .. .. ..@ area : num 11553936 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:27, 1:2] 847187 847400 847672 848025 847942 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 849591 6571460 # .. .. .. ..@ ID : chr "42" # .. .. .. ..@ area : num 11553936 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 874711 6529978 # .. .. .. .. .. .. ..@ area : num 9447399 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:11, 1:2] 872539 875288 876763 875892 875242 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 874711 6529978 # .. .. .. ..@ ID : chr "43" # .. .. .. ..@ area : num 9447399 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 857385 6529257 # .. .. .. .. .. .. ..@ area : num 9274135 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:30, 1:2] 854285 854581 856377 857015 857512 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 857385 6529257 # .. .. .. ..@ ID : chr "44" # .. .. .. ..@ area : num 9274135 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 853627 6587850 # .. .. .. .. .. .. ..@ area : num 9568910 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:21, 1:2] 851742 852642 852719 852566 852783 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 853627 6587850 # .. .. .. ..@ ID : chr "45" # .. .. .. ..@ area : num 9568910 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 890508 6569225 # .. .. .. .. .. .. ..@ area : num 4885442 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:29, 1:2] 888833 888952 889325 889655 889657 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 890508 6569225 # .. .. .. ..@ ID : chr "46" # .. .. .. ..@ area : num 4885442 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 854110 6548980 # .. .. .. .. .. .. ..@ area : num 26171882 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:50, 1:2] 850387 850430 850665 850647 851077 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 854110 6548980 # .. .. .. ..@ ID : chr "47" # .. .. .. ..@ area : num 26171882 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 873126 6569704 # .. .. .. .. .. .. ..@ area : num 24117892 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:42, 1:2] 870128 870302 870138 870252 870536 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 873126 6569704 # .. .. .. ..@ ID : chr "48" # .. .. .. ..@ area : num 24117892 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 866216 6534107 # .. .. .. .. .. .. ..@ area : num 8931113 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:20, 1:2] 863752 864418 864492 864876 865282 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 866216 6534107 # .. .. .. ..@ ID : chr "49" # .. .. .. ..@ area : num 8931113 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 890492 6551443 # .. .. .. .. .. .. ..@ area : num 16899801 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:40, 1:2] 888247 888439 888302 888424 888757 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 890492 6551443 # .. .. .. ..@ ID : chr "50" # .. .. .. ..@ area : num 16899801 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 847006 6591282 # .. .. .. .. .. .. ..@ area : num 7721185 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:27, 1:2] 845004 845121 845325 845640 845614 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 847006 6591282 # .. .. .. ..@ ID : chr "51" # .. .. .. ..@ area : num 7721185 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 903802 6508049 # .. .. .. .. .. .. ..@ area : num 11340432 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:30, 1:2] 901360 902159 902342 902829 903872 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 903802 6508049 # .. .. .. ..@ ID : chr "52" # .. .. .. ..@ area : num 11340432 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 911978 6542601 # .. .. .. .. .. .. ..@ area : num 9690440 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:23, 1:2] 910071 910118 910515 910886 910925 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 911978 6542601 # .. .. .. ..@ ID : chr "53" # .. .. .. ..@ area : num 9690440 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 902018 6556214 # .. .. .. .. .. .. ..@ area : num 24073611 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:44, 1:2] 899548 899823 900627 900782 901187 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 902018 6556214 # .. .. .. ..@ ID : chr "54" # .. .. .. ..@ area : num 24073611 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 909178 6516269 # .. .. .. .. .. .. ..@ area : num 6876333 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:19, 1:2] 908230 908186 908836 908886 909277 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 909178 6516269 # .. .. .. ..@ ID : chr "55" # .. .. .. ..@ area : num 6876333 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 862998 6532148 # .. .. .. .. .. .. ..@ area : num 7525123 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:26, 1:2] 861343 862194 862193 862482 862506 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 862998 6532148 # .. .. .. ..@ ID : chr "56" # .. .. .. ..@ area : num 7525123 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 896452 6566320 # .. .. .. .. .. .. ..@ area : num 4484666 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:15, 1:2] 895254 895419 895986 895872 897836 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 896452 6566320 # .. .. .. ..@ ID : chr "57" # .. .. .. ..@ area : num 4484666 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 893104 6522041 # .. .. .. .. .. .. ..@ area : num 12169331 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:27, 1:2] 890081 890494 890747 890830 891249 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 893104 6522041 # .. .. .. ..@ ID : chr "58" # .. .. .. ..@ area : num 12169331 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 865380 6570105 # .. .. .. .. .. .. ..@ area : num 10230206 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:33, 1:2] 863594 863642 863775 863761 863942 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 865380 6570105 # .. .. .. ..@ ID : chr "59" # .. .. .. ..@ area : num 10230206 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 897504 6532087 # .. .. .. .. .. .. ..@ area : num 10589536 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:34, 1:2] 895578 895700 896386 895834 896547 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 897504 6532087 # .. .. .. ..@ ID : chr "60" # .. .. .. ..@ area : num 10589536 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 892628 6561488 # .. .. .. .. .. .. ..@ area : num 1e+07 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:21, 1:2] 890952 891188 891122 891361 891529 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 892628 6561488 # .. .. .. ..@ ID : chr "61" # .. .. .. ..@ area : num 1e+07 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 891345 6555450 # .. .. .. .. .. .. ..@ area : num 12168498 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:27, 1:2] 888717 889063 889441 890139 890668 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 891345 6555450 # .. .. .. ..@ ID : chr "62" # .. .. .. ..@ area : num 12168498 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 874659 6561803 # .. .. .. .. .. .. ..@ area : num 15793113 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:42, 1:2] 872392 872962 872602 872737 872451 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 874659 6561803 # .. .. .. ..@ ID : chr "63" # .. .. .. ..@ area : num 15793113 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 936783 6584223 # .. .. .. .. .. .. ..@ area : num 6469807 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:27, 1:2] 935011 935628 935288 935600 935782 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 936783 6584223 # .. .. .. ..@ ID : chr "64" # .. .. .. ..@ area : num 6469807 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 878873 6567473 # .. .. .. .. .. .. ..@ area : num 9228563 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:23, 1:2] 876869 876918 877869 880311 880124 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 878873 6567473 # .. .. .. ..@ ID : chr "65" # .. .. .. ..@ area : num 9228563 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 911024 6529628 # .. .. .. .. .. .. ..@ area : num 19738462 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:33, 1:2] 907641 907990 908406 909199 909096 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 911024 6529628 # .. .. .. ..@ ID : chr "66" # .. .. .. ..@ area : num 19738462 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 867860 6547607 # .. .. .. .. .. .. ..@ area : num 3.3e+07 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:56, 1:2] 864240 864496 864382 864587 865137 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 867860 6547607 # .. .. .. ..@ ID : chr "67" # .. .. .. ..@ area : num 3.3e+07 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 839102 6549825 # .. .. .. .. .. .. ..@ area : num 1.7e+07 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:31, 1:2] 836874 836909 837255 837350 837609 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 839102 6549825 # .. .. .. ..@ ID : chr "68" # .. .. .. ..@ area : num 1.7e+07 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 896452 6542556 # .. .. .. .. .. .. ..@ area : num 4749291 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:23, 1:2] 895140 895332 895975 896564 897194 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 896452 6542556 # .. .. .. ..@ ID : chr "69" # .. .. .. ..@ area : num 4749291 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 889740 6561040 # .. .. .. .. .. .. ..@ area : num 7969393 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:17, 1:2] 888054 888497 889010 889141 889637 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 889740 6561040 # .. .. .. ..@ ID : chr "70" # .. .. .. ..@ area : num 7969393 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 929183 6568420 # .. .. .. .. .. .. ..@ area : num 8790767 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:20, 1:2] 927736 928086 928130 927911 927808 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 929183 6568420 # .. .. .. ..@ ID : chr "71" # .. .. .. ..@ area : num 8790767 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 908648 6540466 # .. .. .. .. .. .. ..@ area : num 18211017 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:32, 1:2] 906658 906933 907392 907555 907853 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 908648 6540466 # .. .. .. ..@ ID : chr "72" # .. .. .. ..@ area : num 18211017 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 900142 6550923 # .. .. .. .. .. .. ..@ area : num 31469958 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:26, 1:2] 896889 896893 897182 897683 897986 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 900142 6550923 # .. .. .. ..@ ID : chr "73" # .. .. .. ..@ area : num 31469958 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 917496 6572840 # .. .. .. .. .. .. ..@ area : num 31904867 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:44, 1:2] 913883 914417 914554 914441 915326 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 917496 6572840 # .. .. .. ..@ ID : chr "74" # .. .. .. ..@ area : num 31904867 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 914235 6548862 # .. .. .. .. .. .. ..@ area : num 18216485 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:19, 1:2] 910709 914545 915039 916014 916765 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 914235 6548862 # .. .. .. ..@ ID : chr "75" # .. .. .. ..@ area : num 18216485 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 843273 6556539 # .. .. .. .. .. .. ..@ area : num 12744172 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:25, 1:2] 840484 841058 841426 841428 842112 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 843273 6556539 # .. .. .. ..@ ID : chr "76" # .. .. .. ..@ area : num 12744172 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 856480 6567107 # .. .. .. .. .. .. ..@ area : num 13274413 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:47, 1:2] 854193 854376 854285 854526 854453 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 856480 6567107 # .. .. .. ..@ ID : chr "77" # .. .. .. ..@ area : num 13274413 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 856262 6554614 # .. .. .. .. .. .. ..@ area : num 13632486 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:61, 1:2] 853739 853941 854384 854409 854700 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 856262 6554614 # .. .. .. ..@ ID : chr "78" # .. .. .. ..@ area : num 13632486 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 906727 6568592 # .. .. .. .. .. .. ..@ area : num 18263627 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:28, 1:2] 903748 903901 904779 905127 905845 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 906727 6568592 # .. .. .. ..@ ID : chr "79" # .. .. .. ..@ area : num 18263627 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 872334 6532948 # .. .. .. .. .. .. ..@ area : num 6799025 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:18, 1:2] 870605 871342 871484 871271 871551 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 872334 6532948 # .. .. .. ..@ ID : chr "80" # .. .. .. ..@ area : num 6799025 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 878638 6544077 # .. .. .. .. .. .. ..@ area : num 1.6e+07 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:33, 1:2] 876550 877448 878317 878752 878718 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 878638 6544077 # .. .. .. ..@ ID : chr "81" # .. .. .. ..@ area : num 1.6e+07 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 870514 6549727 # .. .. .. .. .. .. ..@ area : num 15224890 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:35, 1:2] 868537 868733 869026 869156 869294 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 870514 6549727 # .. .. .. ..@ ID : chr "82" # .. .. .. ..@ area : num 15224890 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 913885 6562587 # .. .. .. .. .. .. ..@ area : num 37797358 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:37, 1:2] 910094 910350 910452 910688 911043 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 913885 6562587 # .. .. .. ..@ ID : chr "83" # .. .. .. ..@ area : num 37797358 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 873534 6544153 # .. .. .. .. .. .. ..@ area : num 1.4e+07 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:36, 1:2] 871339 871497 871448 871770 871884 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 873534 6544153 # .. .. .. ..@ ID : chr "84" # .. .. .. ..@ area : num 1.4e+07 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 851186 6559996 # .. .. .. .. .. .. ..@ area : num 1.8e+07 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:44, 1:2] 848717 848902 848743 849116 849128 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 851186 6559996 # .. .. .. ..@ ID : chr "85" # .. .. .. ..@ area : num 1.8e+07 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 854043 6595076 # .. .. .. .. .. .. ..@ area : num 16260714 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:33, 1:2] 851774 851917 851971 852252 852558 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 854043 6595076 # .. .. .. ..@ ID : chr "86" # .. .. .. ..@ area : num 16260714 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 886067 6577794 # .. .. .. .. .. .. ..@ area : num 21699332 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:30, 1:2] 884051 884460 884198 884225 884825 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 886067 6577794 # .. .. .. ..@ ID : chr "87" # .. .. .. ..@ area : num 21699332 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 858535 6568680 # .. .. .. .. .. .. ..@ area : num 16693716 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:44, 1:2] 855288 855556 855558 855803 855812 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 858535 6568680 # .. .. .. ..@ ID : chr "88" # .. .. .. ..@ area : num 16693716 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 911270 6535036 # .. .. .. .. .. .. ..@ area : num 7636816 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:43, 1:2] 909405 909528 909450 909723 910030 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 911270 6535036 # .. .. .. ..@ ID : chr "89" # .. .. .. ..@ area : num 7636816 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 908543 6525001 # .. .. .. .. .. .. ..@ area : num 10500792 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:27, 1:2] 906580 906954 908680 909742 910192 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 908543 6525001 # .. .. .. ..@ ID : chr "90" # .. .. .. ..@ area : num 10500792 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 874857 6535104 # .. .. .. .. .. .. ..@ area : num 2.2e+07 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:33, 1:2] 872432 872651 872985 873954 874572 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 874857 6535104 # .. .. .. ..@ ID : chr "91" # .. .. .. ..@ area : num 2.2e+07 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 902225 6528142 # .. .. .. .. .. .. ..@ area : num 6014209 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:41, 1:2] 899193 899394 899725 899919 900170 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 902225 6528142 # .. .. .. ..@ ID : chr "92" # .. .. .. ..@ area : num 6014209 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 899338 6560022 # .. .. .. .. .. .. ..@ area : num 6656574 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:18, 1:2] 897590 897707 898681 899296 900523 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 899338 6560022 # .. .. .. ..@ ID : chr "93" # .. .. .. ..@ area : num 6656574 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 850833 6588824 # .. .. .. .. .. .. ..@ area : num 17144675 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:35, 1:2] 848373 848568 849400 848820 849030 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 850833 6588824 # .. .. .. ..@ ID : chr "94" # .. .. .. ..@ area : num 17144675 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 934571 6580435 # .. .. .. .. .. .. ..@ area : num 5773268 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:28, 1:2] 933045 933449 933571 933795 933835 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 934571 6580435 # .. .. .. ..@ ID : chr "95" # .. .. .. ..@ area : num 5773268 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 921344 6573226 # .. .. .. .. .. .. ..@ area : num 46311647 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:55, 1:2] 917045 917214 917873 918654 919094 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 921344 6573226 # .. .. .. ..@ ID : chr "96" # .. .. .. ..@ area : num 46311647 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 846189 6537637 # .. .. .. .. .. .. ..@ area : num 20063580 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:30, 1:2] 843257 843566 843466 843664 844305 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 846189 6537637 # .. .. .. ..@ ID : chr "97" # .. .. .. ..@ area : num 20063580 # .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots # .. .. .. ..@ Polygons :List of 1 # .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots # .. .. .. .. .. .. ..@ labpt : num [1:2] 888945 6569690 # .. .. .. .. .. .. ..@ area : num 4605718 # .. .. .. .. .. .. ..@ hole : logi FALSE # .. .. .. .. .. .. ..@ ringDir: int 1 # .. .. .. .. .. .. ..@ coords : num [1:12, 1:2] 887664 888044 888563 888901 888959 ... # .. .. .. ..@ plotOrder: int 1 # .. .. .. ..@ labpt : num [1:2] 888945 6569690 # .. .. .. ..@ ID : chr "98" # .. .. .. ..@ area : num 4605718 # .. .. [list output truncated] # ..@ plotOrder : int [1:35885] 35838 4348 35823 35839 35859 35854 35866 35856 35863 35828 ... # ..@ bbox : num [1:2, 1:2] 98983 5988966 1242433 7110574 # .. ..- attr(*, "dimnames")=List of 2 # .. .. ..$ : chr [1:2] "x" "y" # .. .. ..$ : chr [1:2] "min" "max" # ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slot # .. .. ..@ projargs: chr "+proj=lcc +lat_1=49 +lat_2=44 +lat_0=46.5 +lon_0=3 +x_0=700000 +y_0=6600000 +ellps=GRS80 +units=m +no_defs" ``` --- # Format des couches cartographiques (3/3) ```r library(sf) ``` ```r comm_2016 <- st_read(paste0(getwd(), "/data/couches_geo/COM_2016_CARTElette.shp"), stringsAsFactors = FALSE, quiet = TRUE) ``` ```r str(comm_2016) ``` ``` # Classes 'sf' and 'data.frame': 35885 obs. of 5 variables: # $ insee : chr "01001" "01002" "01004" "01005" ... # $ nom : chr "L'Abergement-Clémenciat" "L'Abergement-de-Varey" "Ambérieu-en-Bugey" "Ambérieux-en-Dombes" ... # $ wikipedia: chr "fr:L'Abergement-Clémenciat" "fr:L'Abergement-de-Varey" "fr:Ambérieu-en-Bugey" "fr:Ambérieux-en-Dombes" ... # $ surf_ha : num 1567 909 2449 1602 597 ... # $ geometry :sfc_MULTIPOLYGON of length 35885; first list element: List of 1 # ..$ :List of 1 # .. ..$ : num [1:44, 1:2] 846839 847056 847011 847281 847063 ... # ..- attr(*, "class")= chr "XY" "MULTIPOLYGON" "sfg" # - attr(*, "sf_column")= chr "geometry" # - attr(*, "agr")= Factor w/ 3 levels "constant","aggregate",..: NA NA NA NA # ..- attr(*, "names")= chr "insee" "nom" "wikipedia" "surf_ha" ``` --- # Créer des couches supra-communales (1/2) <center><img src="img/agreger.png"></center> --- # Créer des couches supra-communales (2/2) ```r # chargement de la couche comm_2016 comm_2016 <- st_read(paste0(getwd(),"/data/couches_geo/COM_2016_CARTElette.shp"), stringsAsFactors = FALSE) %>% merge(table_supracom_2016[, -2], by.x = "insee", by.y = "CODGEO", all.x = T) # fonction pour créer les couches supra-communales creer_couche <- function(maille) { couche <- comm_2016 %>% dplyr::select(maille) %>% group_by(get(maille)) %>% summarize() %>% setNames(nm = c(maille, colnames(.)[-1])) %>% st_simplify(preserveTopology = FALSE, dTolerance = 10) st_write(couche, paste0(getwd(), "/data/couches_geo/", maille, "_2016_CARTElette.shp")) } # appliquer la fonction for (maille in c("DEP", "REG", "CV", "EPCI", "ZE2010", "AU2010", "BV2012", "ARR", "UU2010")) { creer_couche(maille) } ``` --- # Cartes avec cartography (1/4) ```r couche_ze <- st_read(paste0(getwd(), "/data/couches_geo/ZE2010_2016_CARTElette.shp"), stringsAsFactors = FALSE) %>% merge(bdd_ze, by="ZE2010", all.x = TRUE) library(cartography) cartography_rapportHF <- function(annee, breaks, col){ par(mar=c(1,1,1,1)) choroLayer(x = couche_ze, #sf data var = paste0("rapport_FH_", annee), # var in df breaks = breaks, # list of breaks col = col, # colors border = "black", # color of the polygons borders legend.pos = "topleft", # position of the legend legend.title.txt = paste0("Nombre de femmes pour\n100 hommes en ", annee), # title of the legend legend.values.rnd = 0, # number of decimal in the legend values colNA = "#303030", #color of NA values add = FALSE) # add the layer to the current plot ``` --- # Cartes avec cartography (2/4) ```r propSymbolsLayer(x = couche_ze, var = paste0("total_s2_", annee), # population symbols = "circle", # type of symbol border = "grey", # color of the symbols borders col= NA, lwd = 2, # width of the symbols borders legend.pos = "left", legend.title.txt = paste0("Nombre de\nfemmes en ", annee)) layoutLayer(title = paste0("Répartition des femmes en France en ", annee), sources = paste0("Insee, RP",annee,". Fond antuki/CARTElette"), author = "antuki with cartography", scale = 0, #automatic scale bar north = TRUE, frame = TRUE, coltitle = "black", #color of the title col = "#688994" # color of the frame ) } ``` --- # Cartes avec cartography (3/4) ```r cartography_rapportHF(annee = 2014, breaks = quantile(couche_ze$rapport_FH_2014, seq(0, 1, 0.2), na.rm = TRUE) %>% replace(2, 100), col = carto.pal(pal1 = "blue.pal", n1 = 1, pal2 = "red.pal", n2 = 4)) ``` --- # Cartes avec cartography (4/4) <img src="180306_RLadies_COGugaison_carto_files/figure-html/unnamed-chunk-33-1.svg" style="display: block; margin: auto;" /> --- # Cartes avec ggplot2 (1/4) ```r # Créer thème personnalisé theme_map2 <- function() { cowplot::theme_map() %+replace% theme(legend.position = c(0.95, 0.95), legend.justification = c(0, 1), legend.key.width = unit(1, "lines"), legend.text = element_text(size = 10), legend.title = element_text(size = 11), panel.grid.major = element_line(color = "transparent"), plot.title = element_text(size = 12, face = "bold"), plot.caption = element_text(size = 10, colour = "grey20"), plot.margin = unit(c(0, 0, 0, 0), "cm")) } ``` --- # Cartes avec ggplot2 (2/4) ```r ggplot2_rapportHF <- function(annee, breaks, col){ carto <- ggplot( bind_cols(couche_ze %>% mutate(typo = cut(get(paste0("rapport_FH_", annee)), breaks = breaks, dig.lab = 4, include.lowest = TRUE)), couche_ze %>% st_centroid() %>% st_coordinates() %>% as.data.frame() )) + geom_sf(aes(fill = typo), color = "black") + geom_point(aes(x = X, y = Y, size= with(couche_ze, get(paste0("total_s2_", annee)))), shape = 1, fill = NA, color = "grey", stroke = 1.5) + scale_size(name = paste0("Nombre de\nfemmes en ", annee), range = c(0, 18), breaks = c(50000, 100000, 500000),labels = c("50 000", "100 000", "500 000"))+ scale_fill_manual(name = paste0("Nombre de femmes pour\n100 hommes en ", annee), values = col, na.value = "#303030")+ theme_map2() + labs( title = paste0("Répartition des femmes en France en ", annee), caption = paste0("Insee, RP", annee, ". Fond antuki/CARTElette") ) + guides(size = guide_legend(order = 1)) return(carto) } ``` --- # Cartes avec ggplot2 (3/4) ```r ggplot2_rapportHF(annee = 2014, breaks = quantile(couche_ze$rapport_FH_2014, seq(0, 1, 0.2), na.rm = TRUE) %>% replace(2, 100), col = carto.pal(pal1 = "blue.pal", n1 = 1, pal2 = "red.pal", n2 = 4)) ``` --- # Cartes avec ggplot2 (4/4) <img src="180306_RLadies_COGugaison_carto_files/figure-html/unnamed-chunk-37-1.svg" style="display: block; margin: auto;" /> --- # Cartes animées avec magick (1/2) ```r library(magick) img <- image_graph(800, 500, res = 96) out <- lapply(c(1975, 1999, 2014), function(annee) { p <- ggplot2_rapportHF(annee = annee, breaks = quantile(c(couche_ze$rapport_FH_1975, couche_ze$rapport_FH_1999, couche_ze$rapport_FH_2014), seq(0, 1, 0.2), na.rm = TRUE) %>% replace(2, 100), col = carto.pal(pal1 = "blue.pal", n1 = 1, pal2 = "red.pal", n2 = 4)) print(p) }) animation <- image_animate(img, fps = 0.5) print(animation) image_write(animation, "carte_animee.gif") ``` --- # Cartes animées avec magick (2/2) <center><img src="img/carte_animee.gif"></center> --- # Cartes interactives avec ggiraph (1/4) ```r library(rgdal) library(ggiraph) couche_ze_rgdal <- merge(fortify( rgeos::gBuffer(readOGR(dsn = "data/couches_geo", layer = "ZE2010_2016_CARTElette", stringsAsFactors = FALSE), width = 0, byid = TRUE), region = "ZE2010"), bdd_ze, by.x = "id", by.y = "ZE2010", all.x = TRUE) couche_ze_rgdal <- couche_ze_rgdal %>% mutate(typo = cut(get( paste0("rapport_FH_", 2014)), breaks = quantile(couche_ze$rapport_FH_2014, seq(0, 1, 0.2), na.rm = TRUE) %>% replace(2, 100), dig.lab = 4, include.lowest = TRUE)) ``` --- # Cartes interactives avec ggiraph (2/4) ```r ggiraph_rapportHF <- function(annee, breaks, col) { carto <- ggplot(couche_ze_rgdal, aes(x = long, y = lat)) + geom_map_interactive(map = couche_ze_rgdal, mapping = aes(map_id = id, tooltip = paste0("Ratio : ", round(get(paste0("rapport_FH_", annee)), 0)), fill = typo), color = NA) + coord_fixed() + scale_fill_manual(name = paste0("Nombre de femmes pour\n100 hommes en ", annee), values = col, na.value = "#303030") + cowplot::theme_map() + labs(title = paste0("Répartition des femmes en France en ", annee), caption = paste0("Insee, RP", annee, ". Fond antuki/CARTElette")) ggiraph(code = { print(carto) }, tooltip_extra_css = "background-color:black;color:white;font-style: italic;padding:10px;border-radius:10px 20px 10px 20px;", tooltip_opacity = 0.75, width = 0.9, width_svg = 8, height_svg = 6) } ``` --- # Cartes interactives avec ggiraph (3/4) ```r ggiraph_rapportHF(annee = 2014, breaks = quantile(couche_ze$rapport_FH_2014, seq(0, 1, 0.2), na.rm = TRUE) %>% replace(2, 100), col = carto.pal(pal1 = "blue.pal", n1 = 1, pal2 = "red.pal", n2 = 4)) ``` --- # Cartes interactives avec ggiraph (4/4)
--- # Cartes interactives avec leaflet (1/3) ```r library(leaflet) leaflet_rapportHF <- function(){ carto <- leaflet(bind_cols(couche_ze, couche_ze %>% st_centroid() %>% st_transform('+proj=longlat +ellps=GRS80 +no_defs') %>% st_coordinates() %>% as.data.frame()) %>% filter(!substr(ZE2010, 1, 2) %in% c("01", "02", "03", "04"))) %>% addTiles('http://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png') %>% addCircles(lng = ~X, lat = ~Y, weight = 1, radius = ~sqrt(total_s2_2014)*50, popup = ~paste(nom," : ", round(total_s2_2014, 0)), color = ~colorQuantile(carto.pal(pal1 = "blue.pal", n1 = 1, pal2 = "red.pal", n2 = 4), domain=couche_ze$rapport_FH_2014, n = 5) (rapport_FH_2014), fillOpacity = 0.7) ``` --- # Cartes interactives avec leaflet (2/3) ```r carto <- carto %>% addLegend("topright", colors = carto.pal(pal1 = "blue.pal", n1 = 1, pal2 = "red.pal", n2 = 4), labels = paste0("[", round(quantile(couche_ze$rapport_FH_2014, seq(0,1,0.20),na.rm=TRUE), 1), ";", c(round(quantile( couche_ze$rapport_FH_2014,seq(0, 1, 0.20),na.rm=TRUE), 1)[-1], NA),")")[-6], title = "Nombre de<br>femmes pour<br>100 hommes<br>en 2014", opacity = 1) return(carto) } ``` --- # Cartes interactives avec leaflet (3/3)
--- class: center, middle <center><img src="img/avatar.gif" height="300px"></center> <link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/4.7.0/css/font-awesome.min.css"> <a href="http://twitter.com/antuki13" class="social"><i class="fa fa-twitter fa-2x" title="Mon twitter"></i>@antuki13</a> <a href="http://github.com/antuki" class="social"><i class="fa fa-github fa-2x" title="Mon github"></i>antuki</a> <a href="http://antuki.github.io" class="social"><i class="fa fa-bold fa-2x" title="Mon blog"></i>antuki.github.io</a> Diaporama créé avec le package R [**xaringan**](https://github.com/yihui/xaringan) avec le thème [rladies](https://alison.rbind.io/post/r-ladies-slides/) et quelques gifs rigolos de [giphy](https://giphy.com/). Avec également l'utilisation de [remark.js](https://remarkjs.com), [knitr](http://yihui.name/knitr), et [R Markdown](https://rmarkdown.rstudio.com).