Skip to content
GitLab
Commit f8b8d31c authored by Maxime Modjeska's avatar Maxime Modjeska
Browse files

cartes pour pauline

parent 69449077
master map_production v1.0.29.1 v1.0.29.0 v1.0.28.0 v1.0.27.0 v1.0.26.0 v1.0.25.0 v1.0.24.0 v1.0.23.0 v1.0.22.0 v1.0.21.0 v1.0.20.0
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 793 additions and 5 deletions
+793 -5
map/map_entretien.rmd
+ 793
− 5
View file @ f8b8d31c
...@@ -5,13 +5,14 @@ library(sf) ...@@ -5,13 +5,14 @@ library(sf)
# declaration de parametres # declaration de parametres
output_path = "data-common/figure/so-ii/map/" output_path = "data-common/figure/so-ii/map/"
crs = "WGS84"
# Récupération des données (à faire une fois) # Récupération des données (à faire une fois)
ent_prelim = rio::import("data-common/table/so-ii/membre_roi.ods", which = 1) 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 # traitement des données
## sous-ensemble de données ## sous-ensemble de données
...@@ -19,7 +20,7 @@ dataset = ent_prelim[,c("statut", "latitude", "longitude")] ...@@ -19,7 +20,7 @@ dataset = ent_prelim[,c("statut", "latitude", "longitude")]
## géolocalisation des données ## géolocalisation des données
dataset = dataset[!is.na(dataset[["latitude"]]), ] 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 ## mise en forme de la visualisation des données
### symboles et couleur des symboles ### symboles et couleur des symboles
...@@ -28,8 +29,12 @@ pch_hab = 21 ...@@ -28,8 +29,12 @@ pch_hab = 21
pch = rep(NA, nrow(dataset)) pch = rep(NA, nrow(dataset))
pch[dataset[["statut"]] == "agriculteur"] = pch_agri pch[dataset[["statut"]] == "agriculteur"] = pch_agri
pch[dataset[["statut"]] == "habitant"] = pch_hab pch[dataset[["statut"]] == "habitant"] = pch_hab
bg = rep("black", nrow(dataset)) bg = rep(NA, nrow(dataset))
col = rep("white", 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 ### synthèse des données
total_agri = length(dataset[["statut"]][dataset[["statut"]] == "agriculteur"]) total_agri = length(dataset[["statut"]][dataset[["statut"]] == "agriculteur"])
...@@ -43,7 +48,7 @@ dataset_legend = list( ...@@ -43,7 +48,7 @@ dataset_legend = list(
sprintf("Habitant (%s)", total_hab) sprintf("Habitant (%s)", total_hab)
), ),
pch = c(pch_agri,pch_hab), pch = c(pch_agri,pch_hab),
pt.bg = "black", pt.bg = c("green","orange"),
pt.cex = 1.4 pt.cex = 1.4
) )
...@@ -57,3 +62,786 @@ map_so_ii(dataset, ...@@ -57,3 +62,786 @@ map_so_ii(dataset,
path = file.path(output_path, "ensemble_personnes_rencontrees.pdf")) 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"))
```
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment