From 2cae20fb0439f908daf516109c4744bc05ea0339 Mon Sep 17 00:00:00 2001
From: Maxime <maxime.modjeska@inrae.fr>
Date: Mon, 24 Jan 2022 15:46:38 +0100
Subject: [PATCH] MAJ cartes Pauline et Fred

---
 map/map_entretien.rmd |  2 +-
 map/map_rex.rmd       | 80 +++++++++++++++++++++++++++++++++++++++++--
 2 files changed, 78 insertions(+), 4 deletions(-)

diff --git a/map/map_entretien.rmd b/map/map_entretien.rmd
index e6755df..8c60ace 100644
--- a/map/map_entretien.rmd
+++ b/map/map_entretien.rmd
@@ -113,7 +113,7 @@ map_so_ii(dataset,
           path = file.path(output_path, "entretiens_preliminaires.pdf"))
 ```
 
-```{r }
+```{r aleas_roi}
 
 # traitement des données
 
diff --git a/map/map_rex.rmd b/map/map_rex.rmd
index c3b8b95..8b6948b 100644
--- a/map/map_rex.rmd
+++ b/map/map_rex.rmd
@@ -8,7 +8,7 @@ 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)
+rex2020 = rio::import("data-common/table/so-ii/rex-2020-09-19-enquete.ods", which = 1)
 
 ```
 
@@ -64,7 +64,7 @@ dataset_legend = list(
 map_so_ii(dataset, 
           dataset_legend, 
           theme = "catnat",
-          theme_legend = TRUE,
+          #theme_legend = TRUE,
           year = "2020", 
           pch = pch, 
           bg = bg, 
@@ -124,10 +124,84 @@ dataset_legend = list(
 map_so_ii(dataset, 
           dataset_legend, 
           theme = "catnat",
-          theme_legend = TRUE,
+          #theme_legend = TRUE,
           year = "2020", 
           pch = pch, 
           bg = bg, 
           col = col, 
           path = file.path(output_path, "rex-2020-09-19-entretien.pdf"))
+```
+
+```{r alea}
+
+# traitement des données
+
+## sous-ensemble de données
+
+dataset = rex2020[rex2020[["statut_rencontre"]] == "realise" & !is.na(rex2020[["debordement"]]),c(
+    "latitude",
+    "longitude",
+    "statut",
+    "activite",
+    "debordement",
+    "ruissellement"
+    )]
+
+dataset[["multiple"]] = apply(dataset[,c("debordement","ruissellement")],
+      1, 
+      function(x)length(grep("TRUE",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 = "blue"
+bg_multiple = "purple"
+bg[dataset[["ruissellement"]] == "TRUE"] = bg_ruissellement
+bg[dataset[["debordement"]] == "TRUE"] = bg_debordement
+bg[dataset[["multiple"]] > 1] = bg_multiple
+col = rep("black", nrow(dataset))
+
+### synthèse des données
+
+total_ruissellement = length(dataset[["ruissellement"]][dataset[["ruissellement"]] == "TRUE"])
+total_debordement = length(dataset[["debordement"]][dataset[["debordement"]] == "TRUE"])
+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 (%s)",total_debordement), 
+        sprintf("Multiple (%s)",total_multiple)
+        ),
+    pch = 22,
+    pt.bg = c(
+        "red",
+        "blue",
+        "purple"
+    ),
+    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-alea.pdf"))
 ```
\ No newline at end of file
-- 
GitLab