From f8b8d31c28675af36614f085b5e27cad3ecca625 Mon Sep 17 00:00:00 2001
From: Maxime <maxime.modjeska@inrae.fr>
Date: Mon, 24 Jan 2022 12:09:00 +0100
Subject: [PATCH] cartes pour pauline
---
map/map_entretien.rmd | 798 +++++++++++++++++++++++++++++++++++++++++-
1 file changed, 793 insertions(+), 5 deletions(-)
diff --git a/map/map_entretien.rmd b/map/map_entretien.rmd
index 80a2682..e6755df 100644
--- a/map/map_entretien.rmd
+++ b/map/map_entretien.rmd
@@ -5,13 +5,14 @@ library(sf)
# declaration de parametres
output_path = "data-common/figure/so-ii/map/"
+crs = "WGS84"
# Récupération des données (à faire une fois)
ent_prelim = rio::import("data-common/table/so-ii/membre_roi.ods", which = 1)
```
-```{r entretien_prelim}
+```{r ensemble_personnes_rencontrees}
# traitement des données
## sous-ensemble de données
@@ -19,7 +20,7 @@ dataset = ent_prelim[,c("statut", "latitude", "longitude")]
## géolocalisation des données
dataset = dataset[!is.na(dataset[["latitude"]]), ]
-dataset = sf::st_as_sf(dataset, coords = c("longitude", "latitude"), crs = "WGS84")
+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
@@ -28,8 +29,12 @@ pch_hab = 21
pch = rep(NA, nrow(dataset))
pch[dataset[["statut"]] == "agriculteur"] = pch_agri
pch[dataset[["statut"]] == "habitant"] = pch_hab
-bg = rep("black", nrow(dataset))
-col = rep("white", nrow(dataset))
+bg = rep(NA, nrow(dataset))
+bg_agri = "green"
+bg_hab = "orange"
+bg[dataset[["statut"]] == "agriculteur"] = bg_agri
+bg[dataset[["statut"]] == "habitant"] = bg_hab
+col = rep("black", nrow(dataset))
### synthèse des données
total_agri = length(dataset[["statut"]][dataset[["statut"]] == "agriculteur"])
@@ -43,7 +48,7 @@ dataset_legend = list(
sprintf("Habitant (%s)", total_hab)
),
pch = c(pch_agri,pch_hab),
- pt.bg = "black",
+ pt.bg = c("green","orange"),
pt.cex = 1.4
)
@@ -57,3 +62,786 @@ map_so_ii(dataset,
path = file.path(output_path, "ensemble_personnes_rencontrees.pdf"))
```
+```{r entretien_preliminaire}
+# traitement des données
+
+## sous-ensemble de données
+
+dataset = ent_prelim[ent_prelim[["entretien_preliminaire"]]=="oui",c("statut", "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_agri = "green"
+bg_hab = "orange"
+bg[dataset[["statut"]] == "agriculteur"] = bg_agri
+bg[dataset[["statut"]] == "habitant"] = bg_hab
+col = rep("black", nrow(dataset))
+
+### synthèse des données
+total_agri = length(dataset[["statut"]][dataset[["statut"]] == "agriculteur"])
+total_hab = length(dataset[["statut"]][dataset[["statut"]] == "habitant"])
+
+## Légende Carte en cours
+dataset_legend = list(
+ title = "Entretiens préliminaires",
+ legend = c(
+ sprintf("Agriculteur (%s)", total_agri),
+ sprintf("Habitant (%s)", total_hab)
+ ),
+ pch = c(pch_agri,pch_hab),
+ pt.bg = c("green","orange"),
+ pt.cex = 1.4
+)
+
+# Carte
+map_so_ii(dataset,
+ dataset_legend,
+ theme = "clc",
+ pch = pch,
+ bg = bg,
+ col = col,
+ path = file.path(output_path, "entretiens_preliminaires.pdf"))
+```
+
+```{r }
+
+# traitement des données
+
+## sous-ensemble de données
+
+dataset = ent_prelim[,c(
+ "latitude",
+ "longitude",
+ "statut",
+ "ruissellement",
+ "debordement_cours_eau",
+ "debordement_etang",
+ "remontee_nappe",
+ "rupture_digue",
+ "submersion_marine")]
+
+dataset[["multiple"]] = apply(dataset[,c("ruissellement","debordement_cours_eau","debordement_etang","remontee_nappe","rupture_digue","submersion_marine")],
+ 1,
+ function(x)length(grep("oui",x)))
+
+## 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_ruissellement = "red"
+bg_debordement_cours_eau = "blue"
+bg_debordement_etang = "orange2"
+bg_remontee_nappe = "green"
+bg_rupture_digue = "grey"
+bg_submersion_marine = "darkorange3"
+bg_multiple = "purple"
+bg[dataset[["ruissellement"]] == "oui"] = bg_ruissellement
+bg[dataset[["debordement_cours_eau"]] == "oui"] = bg_debordement_cours_eau
+bg[dataset[["debordement_etang"]] == "oui"] = bg_debordement_etang
+bg[dataset[["remontee_nappe"]] == "oui"] = bg_remontee_nappe
+bg[dataset[["rupture_digue"]] == "oui"] = bg_rupture_digue
+bg[dataset[["submersion_marine"]] == "oui"] = bg_submersion_marine
+bg[dataset[["multiple"]] > 1] = bg_multiple
+col = rep("black", nrow(dataset))
+
+### synthèse des données
+
+total_ruissellement = length(dataset[["ruissellement"]][dataset[["ruissellement"]] == "oui"])
+total_debordement_cours_eau = length(dataset[["debordement_cours_eau"]][dataset[["debordement_cours_eau"]] == "oui"])
+total_debordement_etang = length(dataset[["debordement_etang"]][dataset[["debordement_etang"]] == "oui"])
+total_remontee_nappe = length(dataset[["remontee_nappe"]][dataset[["remontee_nappe"]] == "oui"])
+total_rupture_digue = length(dataset[["rupture_digue"]][dataset[["rupture_digue"]] == "oui"])
+total_submersion_marine = length(dataset[["submersion_marine"]][dataset[["submersion_marine"]] == "oui"])
+total_multiple = length(dataset[["multiple"]][dataset[["multiple"]] > 1])
+
+## Légende Carte en cours
+dataset_legend = list(
+ title = "Types d'aléas",
+ legend = c(
+ sprintf("Ruissellement (%s)",total_ruissellement),
+ sprintf("Débordement cours d'eau (%s)",total_debordement_cours_eau),
+ sprintf("Débordement d'étang (%s)",total_debordement_etang),
+ sprintf("Remontée de nappe (%s)",total_remontee_nappe),
+ sprintf("Rupture de digue (%s)",total_rupture_digue),
+ sprintf("Submersion marine (%s)",total_submersion_marine),
+ sprintf("Multiple (%s)",total_multiple)
+ ),
+ pch = c(pch_agri,pch_hab),
+ pt.bg = c(
+ "red",
+ "blue",
+ "orange2",
+ "green",
+ "grey",
+ "darkorange3",
+ "purple"
+ ),
+ pt.cex = 1.4
+)
+
+# Carte
+map_so_ii(dataset,
+ dataset_legend,
+ theme = "hydro",
+ detail = 3,
+ pch = pch,
+ bg = bg,
+ col = col,
+ path = file.path(output_path, "roi_aleas.pdf"))
+```
+
+```{r inondation_1994}
+
+# traitement des données
+
+## sous-ensemble de données
+
+dataset = ent_prelim[ent_prelim[["inondation_1994"]]=="oui",c(
+ "latitude",
+ "longitude",
+ "statut",
+ "inondation_1994")]
+
+## 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_agri = "green"
+bg_hab = "orange"
+bg[dataset[["statut"]] == "agriculteur"] = bg_agri
+bg[dataset[["statut"]] == "habitant"] = bg_hab
+col = rep("black", nrow(dataset))
+
+### synthèse des données
+total_agri = length(dataset[["statut"]][dataset[["statut"]] == "agriculteur" & dataset[["inondation_1994"]] == "oui"])
+total_hab = length(dataset[["statut"]][dataset[["statut"]] == "habitant" & dataset[["inondation_1994"]] == "oui"])
+
+## Légende Carte en cours
+dataset_legend = list(
+ title = "Inondation 1994",
+ legend = c(
+ sprintf("Agriculteur (%s)", total_agri),
+ sprintf("Habitant (%s)", total_hab)
+ ),
+ pch = c(pch_agri,pch_hab),
+ pt.bg = c("green","orange"),
+ pt.cex = 1.4
+)
+
+# Carte
+map_so_ii(dataset,
+ dataset_legend,
+ legend = TRUE,
+ theme = "catnat",
+ year = "1994",
+ pch = pch,
+ bg = bg,
+ col = col,
+ path = file.path(output_path, "inondation_1994.pdf"))
+```
+
+```{r inondation_1996}
+
+# traitement des données
+
+## sous-ensemble de données
+
+dataset = ent_prelim[ent_prelim[["inondation_1996"]]=="oui",c(
+ "latitude",
+ "longitude",
+ "statut",
+ "inondation_1996")]
+
+## 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_agri = "green"
+bg_hab = "orange"
+bg[dataset[["statut"]] == "agriculteur"] = bg_agri
+bg[dataset[["statut"]] == "habitant"] = bg_hab
+col = rep("black", nrow(dataset))
+
+### synthèse des données
+total_agri = length(dataset[["statut"]][dataset[["statut"]] == "agriculteur" & dataset[["inondation_1996"]] == "oui"])
+total_hab = length(dataset[["statut"]][dataset[["statut"]] == "habitant" & dataset[["inondation_1996"]] == "oui"])
+
+## Légende Carte en cours
+dataset_legend = list(
+ title = "Inondation 1996",
+ legend = c(
+ sprintf("Agriculteur (%s)", total_agri),
+ sprintf("Habitant (%s)", total_hab)
+ ),
+ pch = c(pch_agri,pch_hab),
+ pt.bg = c("green","orange"),
+ pt.cex = 1.4
+)
+
+# Carte
+map_so_ii(dataset,
+ dataset_legend,
+ legend = TRUE,
+ theme = "catnat",
+ year = "1996",
+ pch = pch,
+ bg = bg,
+ col = col,
+ path = file.path(output_path, "inondation_1996.pdf"))
+```
+
+```{r inondation_1999}
+
+# traitement des données
+
+## sous-ensemble de données
+
+dataset = ent_prelim[ent_prelim[["inondation_1999"]]=="oui",c(
+ "latitude",
+ "longitude",
+ "statut",
+ "inondation_1999")]
+
+## 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_agri = "green"
+bg_hab = "orange"
+bg[dataset[["statut"]] == "agriculteur"] = bg_agri
+bg[dataset[["statut"]] == "habitant"] = bg_hab
+col = rep("black", nrow(dataset))
+
+### synthèse des données
+total_agri = length(dataset[["statut"]][dataset[["statut"]] == "agriculteur" & dataset[["inondation_1999"]] == "oui"])
+total_hab = length(dataset[["statut"]][dataset[["statut"]] == "habitant" & dataset[["inondation_1999"]] == "oui"])
+
+## Légende Carte en cours
+dataset_legend = list(
+ title = "Inondation 1999",
+ legend = c(
+ sprintf("Agriculteur (%s)", total_agri),
+ sprintf("Habitant (%s)", total_hab)
+ ),
+ pch = c(pch_agri,pch_hab),
+ pt.bg = c("green","orange"),
+ pt.cex = 1.4
+)
+
+# Carte
+map_so_ii(dataset,
+ dataset_legend,
+ legend = TRUE,
+ theme = "catnat",
+ year = "1999",
+ pch = pch,
+ bg = bg,
+ col = col,
+ path = file.path(output_path, "inondation_1999.pdf"))
+```
+
+```{r inondation_2002}
+
+# traitement des données
+
+## sous-ensemble de données
+
+dataset = ent_prelim[ent_prelim[["inondation_2002"]]=="oui",c(
+ "latitude",
+ "longitude",
+ "statut",
+ "inondation_2002")]
+
+## 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_agri = "green"
+bg_hab = "orange"
+bg[dataset[["statut"]] == "agriculteur"] = bg_agri
+bg[dataset[["statut"]] == "habitant"] = bg_hab
+col = rep("black", nrow(dataset))
+
+### synthèse des données
+total_agri = length(dataset[["statut"]][dataset[["statut"]] == "agriculteur" & dataset[["inondation_2002"]] == "oui"])
+total_hab = length(dataset[["statut"]][dataset[["statut"]] == "habitant" & dataset[["inondation_2002"]] == "oui"])
+
+## Légende Carte en cours
+dataset_legend = list(
+ title = "Inondation 2002",
+ legend = c(
+ sprintf("Agriculteur (%s)", total_agri),
+ sprintf("Habitant (%s)", total_hab)
+ ),
+ pch = c(pch_agri,pch_hab),
+ pt.bg = c("green","orange"),
+ pt.cex = 1.4
+)
+
+# Carte
+map_so_ii(dataset,
+ dataset_legend,
+ legend = TRUE,
+ theme = "catnat",
+ year = "2002",
+ pch = pch,
+ bg = bg,
+ col = col,
+ path = file.path(output_path, "inondation_2002.pdf"))
+```
+
+```{r inondation_2003}
+
+# traitement des données
+
+## sous-ensemble de données
+
+dataset = ent_prelim[ent_prelim[["inondation_2003"]]=="oui",c(
+ "latitude",
+ "longitude",
+ "statut",
+ "inondation_2003")]
+
+## 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_agri = "green"
+bg_hab = "orange"
+bg[dataset[["statut"]] == "agriculteur"] = bg_agri
+bg[dataset[["statut"]] == "habitant"] = bg_hab
+col = rep("black", nrow(dataset))
+
+### synthèse des données
+total_agri = length(dataset[["statut"]][dataset[["statut"]] == "agriculteur" & dataset[["inondation_2003"]] == "oui"])
+total_hab = length(dataset[["statut"]][dataset[["statut"]] == "habitant" & dataset[["inondation_2003"]] == "oui"])
+
+## Légende Carte en cours
+dataset_legend = list(
+ title = "Inondation 2003",
+ legend = c(
+ sprintf("Agriculteur (%s)", total_agri),
+ sprintf("Habitant (%s)", total_hab)
+ ),
+ pch = c(pch_agri,pch_hab),
+ pt.bg = c("green","orange"),
+ pt.cex = 1.4
+)
+
+# Carte
+map_so_ii(dataset,
+ dataset_legend,
+ legend = TRUE,
+ theme = "catnat",
+ year = "2003",
+ pch = pch,
+ bg = bg,
+ col = col,
+ path = file.path(output_path, "inondation_2003.pdf"))
+```
+
+```{r inondation_2004}
+
+# traitement des données
+
+## sous-ensemble de données
+
+dataset = ent_prelim[ent_prelim[["inondation_2004"]]=="oui",c(
+ "latitude",
+ "longitude",
+ "statut",
+ "inondation_2004")]
+
+## 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_agri = "green"
+bg_hab = "orange"
+bg[dataset[["statut"]] == "agriculteur"] = bg_agri
+bg[dataset[["statut"]] == "habitant"] = bg_hab
+col = rep("black", nrow(dataset))
+
+### synthèse des données
+total_agri = length(dataset[["statut"]][dataset[["statut"]] == "agriculteur" & dataset[["inondation_2004"]] == "oui"])
+total_hab = length(dataset[["statut"]][dataset[["statut"]] == "habitant" & dataset[["inondation_2004"]] == "oui"])
+
+## Légende Carte en cours
+dataset_legend = list(
+ title = "Inondation 2004",
+ legend = c(
+ sprintf("Agriculteur (%s)", total_agri),
+ sprintf("Habitant (%s)", total_hab)
+ ),
+ pch = c(pch_agri,pch_hab),
+ pt.bg = c("green","orange"),
+ pt.cex = 1.4
+)
+
+# Carte
+map_so_ii(dataset,
+ dataset_legend,
+ legend = TRUE,
+ theme = "catnat",
+ year = "2004",
+ pch = pch,
+ bg = bg,
+ col = col,
+ path = file.path(output_path, "inondation_2004.pdf"))
+```
+
+```{r inondation_2005}
+
+# traitement des données
+
+## sous-ensemble de données
+
+dataset = ent_prelim[ent_prelim[["inondation_2005"]]=="oui",c(
+ "latitude",
+ "longitude",
+ "statut",
+ "inondation_2005")]
+
+## 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_agri = "green"
+bg_hab = "orange"
+bg[dataset[["statut"]] == "agriculteur"] = bg_agri
+bg[dataset[["statut"]] == "habitant"] = bg_hab
+col = rep("black", nrow(dataset))
+
+### synthèse des données
+total_agri = length(dataset[["statut"]][dataset[["statut"]] == "agriculteur" & dataset[["inondation_2005"]] == "oui"])
+total_hab = length(dataset[["statut"]][dataset[["statut"]] == "habitant" & dataset[["inondation_2005"]] == "oui"])
+
+## Légende Carte en cours
+dataset_legend = list(
+ title = "Inondation 2005",
+ legend = c(
+ sprintf("Agriculteur (%s)", total_agri),
+ sprintf("Habitant (%s)", total_hab)
+ ),
+ pch = c(pch_agri,pch_hab),
+ pt.bg = c("green","orange"),
+ pt.cex = 1.4
+)
+
+# Carte
+map_so_ii(dataset,
+ dataset_legend,
+ legend = TRUE,
+ theme = "catnat",
+ year = "2005",
+ pch = pch,
+ bg = bg,
+ col = col,
+ path = file.path(output_path, "inondation_2005.pdf"))
+```
+
+```{r inondation_2007}
+
+# traitement des données
+
+## sous-ensemble de données
+
+dataset = ent_prelim[ent_prelim[["inondation_2007"]]=="oui",c(
+ "latitude",
+ "longitude",
+ "statut",
+ "inondation_2007")]
+
+## 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_agri = "green"
+bg_hab = "orange"
+bg[dataset[["statut"]] == "agriculteur"] = bg_agri
+bg[dataset[["statut"]] == "habitant"] = bg_hab
+col = rep("black", nrow(dataset))
+
+### synthèse des données
+total_agri = length(dataset[["statut"]][dataset[["statut"]] == "agriculteur" & dataset[["inondation_2007"]] == "oui"])
+total_hab = length(dataset[["statut"]][dataset[["statut"]] == "habitant" & dataset[["inondation_2007"]] == "oui"])
+
+## Légende Carte en cours
+dataset_legend = list(
+ title = "Inondation 2007",
+ legend = c(
+ sprintf("Agriculteur (%s)", total_agri),
+ sprintf("Habitant (%s)", total_hab)
+ ),
+ pch = c(pch_agri,pch_hab),
+ pt.bg = c("green","orange"),
+ pt.cex = 1.4
+)
+
+# Carte
+map_so_ii(dataset,
+ dataset_legend,
+ legend = TRUE,
+ theme = "catnat",
+ year = "2007",
+ pch = pch,
+ bg = bg,
+ col = col,
+ path = file.path(output_path, "inondation_2007.pdf"))
+```
+
+```{r inondation_2014}
+
+# traitement des données
+
+## sous-ensemble de données
+
+dataset = ent_prelim[ent_prelim[["inondation_2014"]]=="oui",c(
+ "latitude",
+ "longitude",
+ "statut",
+ "inondation_2014")]
+
+## 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_agri = "green"
+bg_hab = "orange"
+bg[dataset[["statut"]] == "agriculteur"] = bg_agri
+bg[dataset[["statut"]] == "habitant"] = bg_hab
+col = rep("black", nrow(dataset))
+
+### synthèse des données
+total_agri = length(dataset[["statut"]][dataset[["statut"]] == "agriculteur" & dataset[["inondation_2014"]] == "oui"])
+total_hab = length(dataset[["statut"]][dataset[["statut"]] == "habitant" & dataset[["inondation_2014"]] == "oui"])
+
+## Légende Carte en cours
+dataset_legend = list(
+ title = "Inondation 2014",
+ legend = c(
+ sprintf("Agriculteur (%s)", total_agri),
+ sprintf("Habitant (%s)", total_hab)
+ ),
+ pch = c(pch_agri,pch_hab),
+ pt.bg = c("green","orange"),
+ pt.cex = 1.4
+)
+
+# Carte
+map_so_ii(dataset,
+ dataset_legend,
+ legend = TRUE,
+ theme = "catnat",
+ year = "2014",
+ pch = pch,
+ bg = bg,
+ col = col,
+ path = file.path(output_path, "inondation_2014.pdf"))
+```
+
+```{r inondation_2015}
+
+# traitement des données
+
+## sous-ensemble de données
+
+dataset = ent_prelim[ent_prelim[["inondation_2015"]]=="oui",c(
+ "latitude",
+ "longitude",
+ "statut",
+ "inondation_2015")]
+
+## 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_agri = "green"
+bg_hab = "orange"
+bg[dataset[["statut"]] == "agriculteur"] = bg_agri
+bg[dataset[["statut"]] == "habitant"] = bg_hab
+col = rep("black", nrow(dataset))
+
+### synthèse des données
+total_agri = length(dataset[["statut"]][dataset[["statut"]] == "agriculteur" & dataset[["inondation_2015"]] == "oui"])
+total_hab = length(dataset[["statut"]][dataset[["statut"]] == "habitant" & dataset[["inondation_2015"]] == "oui"])
+
+## Légende Carte en cours
+dataset_legend = list(
+ title = "Inondation 2015",
+ legend = c(
+ sprintf("Agriculteur (%s)", total_agri),
+ sprintf("Habitant (%s)", total_hab)
+ ),
+ pch = c(pch_agri,pch_hab),
+ pt.bg = c("green","orange"),
+ pt.cex = 1.4
+)
+
+# Carte
+map_so_ii(dataset,
+ dataset_legend,
+ legend = TRUE,
+ theme = "catnat",
+ year = "2015",
+ pch = pch,
+ bg = bg,
+ col = col,
+ path = file.path(output_path, "inondation_2015.pdf"))
+```
+
+```{r inondation_2020}
+
+# traitement des données
+
+## sous-ensemble de données
+
+dataset = ent_prelim[ent_prelim[["inondation_2020"]]=="oui",c(
+ "latitude",
+ "longitude",
+ "statut",
+ "inondation_2020")]
+
+## 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_agri = "green"
+bg_hab = "orange"
+bg[dataset[["statut"]] == "agriculteur"] = bg_agri
+bg[dataset[["statut"]] == "habitant"] = bg_hab
+col = rep("black", nrow(dataset))
+
+### synthèse des données
+total_agri = length(dataset[["statut"]][dataset[["statut"]] == "agriculteur" & dataset[["inondation_2020"]] == "oui"])
+total_hab = length(dataset[["statut"]][dataset[["statut"]] == "habitant" & dataset[["inondation_2020"]] == "oui"])
+
+## Légende Carte en cours
+dataset_legend = list(
+ title = "Inondation 2020",
+ legend = c(
+ sprintf("Agriculteur (%s)", total_agri),
+ sprintf("Habitant (%s)", total_hab)
+ ),
+ pch = c(pch_agri,pch_hab),
+ pt.bg = c("green","orange"),
+ pt.cex = 1.4
+)
+
+# Carte
+map_so_ii(dataset,
+ dataset_legend,
+ legend = TRUE,
+ theme = "catnat",
+ year = "2020",
+ pch = pch,
+ bg = bg,
+ col = col,
+ path = file.path(output_path, "inondation_2020.pdf"))
+```
--
GitLab