diff --git a/map/map_rex.rmd b/map/map_rex.rmd new file mode 100644 index 0000000000000000000000000000000000000000..c3b8b95b624a1bfa03c9acf3bbafa2ff61ae64a3 --- /dev/null +++ b/map/map_rex.rmd @@ -0,0 +1,133 @@ +```{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