Présentation

"Cette note est un essai de la réalisation de la carte thématique avec le package mapsf"
[1] "Cette note est un essai de la réalisation de la carte thématique avec le package mapsf"

Chargement des packages nécessaires

library(sf)
library(s2)
library(mapsf)
library(raster)
library(tmap)
library(tidyverse)

Définition de l’espace de travail

setwd("E:/Logiciel R/CARTOGRAPHIE AVEC R/TP CARTO")
cg <- st_read("Departements_du_CongoV.shp", quiet=TRUE)

Importation des données

library(readxl)
centres <- read_excel("Centres d'état civil VF 29_06_2022.xlsx")
cg1 <- st_read("Departements_du_CongoV.shp", quiet=TRUE) 
d <- select(centres,"Département","Centres_principaux" = "Centres principaux", "Centres_secondaires" = "Centres secondaires")
cg1 <- left_join(cg,d, by = c("Noms_Dép" = "Département"))

Suppression de certaines variables

#Suppresion de certaines variables----
cg1$OBJECTID <- NULL
cg1$Id <- NULL
cg1$Ojectid <- NULL
cg1$Centres_secondaires
 [1] 219 452 260 419 496 517 154 170 123 225 628 113
names(cg1)
sf::sf_use_s2 (FALSE) 
Spherical geometry (s2) switched off

Production de la carte

# Ma carte réussi----

#mf_init(x = d2, theme = "candy")
#expandBB valeurs fractionnaires pour développer le cadre de sélection avec, dans chaque direction (bas, gauche, haut, droite)
mf_init(x = cg1, theme = "candy", expandBB = c(0,0,0,0))
#mf_shadow(d2, add = TRUE)
mf_shadow(cg1, add = TRUE)
mf_map(cg1, add = TRUE)
mf_map(cg1, var = "Centres_secondaires",type = "prop",
      leg_pos = "right",
      add = TRUE)
Warning in st_centroid.sfc(x = st_geometry(x), of_largest_polygon = TRUE) :
  st_centroid does not give correct centroids for longitude/latitude data
mf_label(cg1, var ="Noms_Dép",halo = FALSE, cex = 0.8,
         overlap = FALSE,
         col = "black",
         lines = TRUE)
Warning in st_centroid.sfc(x = sf::st_geometry(x), of_largest_polygon = max(sf::st_is(sf::st_as_sf(x),  :
  st_centroid does not give correct centroids for longitude/latitude data
mf_arrow(pos = "topright")
#mf_legend(type = "prop", pos = "topright")#,val = c(1, 5, 10), inches = .3)
mf_layout(title = "Cartographie des centres d'état civil du Congo",
          arrow = FALSE,
          credits = paste0("Sources: INS, 2022\n", 
                           "mapsf ", 
                           packageVersion("mapsf")))

#mf_legend(type = "prop", pos = "topright", val = c(1, 5, 10), inches = .3)
mf_label(cg1, var ="Noms_Dép",halo = FALSE, cex = 0.8,
         overlap = FALSE,
         col = "black",
         lines = TRUE, add = TRUE)
Warning in st_centroid.sfc(x = sf::st_geometry(x), of_largest_polygon = max(sf::st_is(sf::st_as_sf(x),  :
  st_centroid does not give correct centroids for longitude/latitude data
Error in strwidth(words[i], cex = cex[i], ...) : 
  plot.new has not been called yet
LS0tDQp0aXRsZTogIkNBUlRPR1JBUEhJRSBERVMgQ0VOVFJFUyBEJ0VUQVQgQ0lWSUwgRFUgQ09OR08iDQpvdXRwdXQ6DQogIGh0bWxfbm90ZWJvb2s6IA0KICAgIHRvYzogeWVzDQogICAgZmlnX3dpZHRoOiAyNg0KICAgIGZpZ19oZWlnaHQ6IDE5DQotLS0NCg0KIyBQcsOpc2VudGF0aW9uDQoNCmBgYHtyfQ0KIkNldHRlIG5vdGUgZXN0IHVuIGVzc2FpIGRlIGxhIHLDqWFsaXNhdGlvbiBkZSBsYSBjYXJ0ZSB0aMOpbWF0aXF1ZSBhdmVjIGxlIHBhY2thZ2UgbWFwc2YiDQpgYGANCg0KIyBDaGFyZ2VtZW50IGRlcyBwYWNrYWdlcyBuw6ljZXNzYWlyZXMNCg0KYGBge3J9DQpsaWJyYXJ5KHNmKQ0KbGlicmFyeShzMikNCmxpYnJhcnkobWFwc2YpDQpsaWJyYXJ5KHJhc3RlcikNCmxpYnJhcnkodG1hcCkNCmxpYnJhcnkodGlkeXZlcnNlKQ0KYGBgDQoNCiMgRMOpZmluaXRpb24gZGUgbCdlc3BhY2UgZGUgdHJhdmFpbA0KDQpgYGB7cn0NCnNldHdkKCJFOi9Mb2dpY2llbCBSL0NBUlRPR1JBUEhJRSBBVkVDIFIvVFAgQ0FSVE8iKQ0KY2cgPC0gc3RfcmVhZCgiRGVwYXJ0ZW1lbnRzX2R1X0NvbmdvVi5zaHAiLCBxdWlldD1UUlVFKQ0KYGBgDQoNCiMgSW1wb3J0YXRpb24gZGVzIGRvbm7DqWVzDQoNCmBgYHtyfQ0KbGlicmFyeShyZWFkeGwpDQpjZW50cmVzIDwtIHJlYWRfZXhjZWwoIkNlbnRyZXMgZCfDqXRhdCBjaXZpbCBWRiAyOV8wNl8yMDIyLnhsc3giKQ0KY2cxIDwtIHN0X3JlYWQoIkRlcGFydGVtZW50c19kdV9Db25nb1Yuc2hwIiwgcXVpZXQ9VFJVRSkgDQpkIDwtIHNlbGVjdChjZW50cmVzLCJEw6lwYXJ0ZW1lbnQiLCJDZW50cmVzX3ByaW5jaXBhdXgiID0gIkNlbnRyZXMgcHJpbmNpcGF1eCIsICJDZW50cmVzX3NlY29uZGFpcmVzIiA9ICJDZW50cmVzIHNlY29uZGFpcmVzIikNCmNnMSA8LSBsZWZ0X2pvaW4oY2csZCwgYnkgPSBjKCJOb21zX0TDqXAiID0gIkTDqXBhcnRlbWVudCIpKQ0KYGBgDQoNCiMgU3VwcHJlc3Npb24gZGUgY2VydGFpbmVzIHZhcmlhYmxlcw0KDQpgYGB7cn0NCiNTdXBwcmVzaW9uIGRlIGNlcnRhaW5lcyB2YXJpYWJsZXMtLS0tDQpjZzEkT0JKRUNUSUQgPC0gTlVMTA0KY2cxJElkIDwtIE5VTEwNCmNnMSRPamVjdGlkIDwtIE5VTEwNCmBgYA0KDQpgYGB7cn0NCmNnMSRDZW50cmVzX3NlY29uZGFpcmVzDQpgYGANCmBgYHtyfQ0KbmFtZXMoY2cxKQ0KYGBgDQoNCmBgYHtyfQ0Kc2Y6OnNmX3VzZV9zMiAoRkFMU0UpIA0KYGBgDQoNCg0KIyBQcm9kdWN0aW9uIGRlIGxhIGNhcnRlDQoNCmBgYHtyfQ0KIyBNYSBjYXJ0ZSByw6l1c3NpLS0tLQ0KDQojbWZfaW5pdCh4ID0gZDIsIHRoZW1lID0gImNhbmR5IikNCiNleHBhbmRCQiB2YWxldXJzIGZyYWN0aW9ubmFpcmVzIHBvdXIgZMOpdmVsb3BwZXIgbGUgY2FkcmUgZGUgc8OpbGVjdGlvbiBhdmVjLCBkYW5zIGNoYXF1ZSBkaXJlY3Rpb24gKGJhcywgZ2F1Y2hlLCBoYXV0LCBkcm9pdGUpDQptZl9pbml0KHggPSBjZzEsIHRoZW1lID0gImNhbmR5IiwgZXhwYW5kQkIgPSBjKDAsMCwwLDApKQ0KI21mX3NoYWRvdyhkMiwgYWRkID0gVFJVRSkNCm1mX3NoYWRvdyhjZzEsIGFkZCA9IFRSVUUpDQptZl9tYXAoY2cxLCBhZGQgPSBUUlVFKQ0KbWZfbWFwKGNnMSwgdmFyID0gIkNlbnRyZXNfc2Vjb25kYWlyZXMiLHR5cGUgPSAicHJvcCIsDQogICAgICBsZWdfcG9zID0gInJpZ2h0IiwNCiAgICAgIGFkZCA9IFRSVUUpDQptZl9sYWJlbChjZzEsIHZhciA9Ik5vbXNfRMOpcCIsaGFsbyA9IEZBTFNFLCBjZXggPSAwLjgsDQogICAgICAgICBvdmVybGFwID0gRkFMU0UsDQogICAgICAgICBjb2wgPSAiYmxhY2siLA0KICAgICAgICAgbGluZXMgPSBUUlVFKQ0KbWZfYXJyb3cocG9zID0gInRvcHJpZ2h0IikNCiNtZl9sZWdlbmQodHlwZSA9ICJwcm9wIiwgcG9zID0gInRvcHJpZ2h0IikjLHZhbCA9IGMoMSwgNSwgMTApLCBpbmNoZXMgPSAuMykNCm1mX2xheW91dCh0aXRsZSA9ICJDYXJ0b2dyYXBoaWUgZGVzIGNlbnRyZXMgZCfDqXRhdCBjaXZpbCBkdSBDb25nbyIsDQogICAgICAgICAgYXJyb3cgPSBGQUxTRSwNCiAgICAgICAgICBjcmVkaXRzID0gcGFzdGUwKCJTb3VyY2VzOiBJTlMsIDIwMjJcbiIsIA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgIm1hcHNmICIsIA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgcGFja2FnZVZlcnNpb24oIm1hcHNmIikpKQ0KDQpgYGANCg0KYGBge3J9DQojbWZfbGVnZW5kKHR5cGUgPSAicHJvcCIsIHBvcyA9ICJ0b3ByaWdodCIsIHZhbCA9IGMoMSwgNSwgMTApLCBpbmNoZXMgPSAuMykNCm1mX2xhYmVsKGNnMSwgdmFyID0iTm9tc19Ew6lwIixoYWxvID0gRkFMU0UsIGNleCA9IDAuOCwNCiAgICAgICAgIG92ZXJsYXAgPSBGQUxTRSwNCiAgICAgICAgIGNvbCA9ICJibGFjayIsDQogICAgICAgICBsaW5lcyA9IFRSVUUpDQptZl9hcnJvdyhwb3MgPSAidG9wcmlnaHQiKQ0KI21mX2xlZ2VuZCh0eXBlID0gInByb3AiLCBwb3MgPSAidG9wcmlnaHQiKSMsdmFsID0gYygxLCA1LCAxMCksIGluY2hlcyA9IC4zKQ0KbWZfbGF5b3V0KHRpdGxlID0gIkNhcnRvZ3JhcGhpZSBkZXMgY2VudHJlcyBkJ8OpdGF0IGNpdmlsIGR1IENvbmdvIiwNCiAgICAgICAgICBhcnJvdyA9IEZBTFNFLA0KICAgICAgICAgIGNyZWRpdHMgPSBwYXN0ZTAoIlNvdXJjZXM6IElOUywgMjAyMlxuIiwgDQogICAgICAgICAgICAgICAgICAgICAgICAgICAibWFwc2YgIiwgDQogICAgICAgICAgICAgICAgICAgICAgICAgICBwYWNrYWdlVmVyc2lvbigibWFwc2YiKSkpDQpgYGANCg0K