map rex pour Fred

map/map_rex.rmd

+```{r data-recup}
+# chargement des bibliothèques
+library(geau)
+library(sf)
+
+# declaration de parametres
+output_path = "data-common/figure/so-ii/map/"
+crs = "WGS84"
+
+# Récupération des données (à faire une fois)
+rex2020 = rio::import("data-common/table/so-ii/rex-2020-09-19.ods", which = 1)
+
+```
+
+```{r rex_contact_pris}
+# traitement des données
+
+## sous-ensemble de données
+dataset = rex2020[,c("statut","activite", "latitude", "longitude")]
+
+## géolocalisation des données
+dataset = dataset[!is.na(dataset[["latitude"]]), ]
+dataset = sf::st_as_sf(dataset, coords = c("longitude", "latitude"), crs = crs)
+
+## mise en forme de la visualisation des données
+### symboles et couleur des symboles
+pch_agri = 24 
+pch_hab = 21
+pch = rep(NA, nrow(dataset))
+pch[dataset[["statut"]] == "agriculteur"] = pch_agri
+pch[dataset[["statut"]] == "habitant"] = pch_hab
+bg = rep(NA, nrow(dataset))
+bg_viti = "deeppink4"
+bg_marai = "green"
+bg_cereale = "yellow"
+bg_hab = "orange"
+bg[dataset[["activite"]] == "viticulture"] = bg_viti
+bg[dataset[["activite"]] == "maraichage"] = bg_marai
+bg[dataset[["activite"]] == "cereale"] = bg_cereale
+bg[dataset[["statut"]] == "habitant"] = bg_hab
+col = rep("black", nrow(dataset))
+
+### synthèse des données
+total_viti = length(dataset[["activite"]][dataset[["activite"]] == "viticulture"])
+total_marai = length(dataset[["activite"]][dataset[["activite"]] == "maraichage"])
+total_cereale = length(dataset[["activite"]][dataset[["activite"]] == "cereale"])
+total_hab = length(dataset[["statut"]][dataset[["statut"]] == "habitant"])
+
+## Légende Carte en cours
+dataset_legend = list(
+    title = "Contacts pris",
+    legend = c(
+        sprintf("Viticulture (%s)", total_viti), 
+        sprintf("Maraichage (%s)", total_marai),
+        sprintf("Grandes cultures (%s)", total_cereale), 
+        sprintf("Habitant (%s)", total_hab)
+        ),
+    pch = c(pch_agri,pch_agri,pch_agri,pch_hab),
+    pt.bg = c("deeppink4","green","yellow","orange"),
+    pt.cex = 1.4
+)
+
+# Carte
+map_so_ii(dataset, 
+          dataset_legend, 
+          theme = "catnat",
+          theme_legend = TRUE,
+          year = "2020", 
+          pch = pch, 
+          bg = bg, 
+          col = col, 
+          path = file.path(output_path, "rex-2020-09-19-contact.pdf"))
+```
+
+```{r rex_entretiens_realises}
+# traitement des données
+
+## sous-ensemble de données
+dataset = rex2020[rex2020[["statut_rencontre"]] == "realise",c("statut","activite", "latitude", "longitude")]
+
+## géolocalisation des données
+dataset = dataset[!is.na(dataset[["latitude"]]), ]
+dataset = sf::st_as_sf(dataset, coords = c("longitude", "latitude"), crs = crs)
+
+## mise en forme de la visualisation des données
+### symboles et couleur des symboles
+pch_agri = 24 
+pch_hab = 21
+pch = rep(NA, nrow(dataset))
+pch[dataset[["statut"]] == "agriculteur"] = pch_agri
+pch[dataset[["statut"]] == "habitant"] = pch_hab
+bg = rep(NA, nrow(dataset))
+bg_viti = "deeppink4"
+bg_marai = "green"
+bg_cereale = "yellow"
+bg_hab = "orange"
+bg[dataset[["activite"]] == "viticulture"] = bg_viti
+bg[dataset[["activite"]] == "maraichage"] = bg_marai
+bg[dataset[["activite"]] == "cereale"] = bg_cereale
+bg[dataset[["statut"]] == "habitant"] = bg_hab
+col = rep("black", nrow(dataset))
+
+### synthèse des données
+total_viti = length(dataset[["activite"]][dataset[["activite"]] == "viticulture"])
+total_marai = length(dataset[["activite"]][dataset[["activite"]] == "maraichage"])
+total_cereale = length(dataset[["activite"]][dataset[["activite"]] == "cereale"])
+total_hab = length(dataset[["statut"]][dataset[["statut"]] == "habitant"])
+
+## Légende Carte en cours
+dataset_legend = list(
+    title = "Contacts pris",
+    legend = c(
+        sprintf("Viticulture (%s)", total_viti), 
+        sprintf("Maraichage (%s)", total_marai),
+        sprintf("Grandes cultures (%s)", total_cereale), 
+        sprintf("Habitant (%s)", total_hab)
+        ),
+    pch = c(pch_agri,pch_agri,pch_agri,pch_hab),
+    pt.bg = c("deeppink4","green","yellow","orange"),
+    pt.cex = 1.4
+)
+
+# Carte
+map_so_ii(dataset, 
+          dataset_legend, 
+          theme = "catnat",
+          theme_legend = TRUE,
+          year = "2020", 
+          pch = pch, 
+          bg = bg, 
+          col = col, 
+          path = file.path(output_path, "rex-2020-09-19-entretien.pdf"))
+```
\ No newline at end of file