Commit b8fadc2b authored by Dorchies David's avatar Dorchies David
Browse files

feat(delta): add maps

Refs #12
parent 2a9a711b
Pipeline #32041 failed with stages
in 36 minutes and 44 seconds
...@@ -244,11 +244,29 @@ Exemple sur tous les indicateurs hydrologiques à Paris: ...@@ -244,11 +244,29 @@ Exemple sur tous les indicateurs hydrologiques à Paris:
deltaIndicators <- calcDelta(Q_indicators, rcp = "rcp4.5", period = "end") deltaIndicators <- calcDelta(Q_indicators, rcp = "rcp4.5", period = "end")
knitr::kable(tableDeltaStation("H8012010", deltaIndicators), digits = 1) knitr::kable(tableDeltaStation("H8012010", deltaIndicators), digits = 1)
``` ```
## Cartes d'évolutions
Exemple sur le QJXA10: ```{r}
plot_map_delta <- function(r, title) {
# Define scale centered around 0
decVal <- r[r < 0]
incVal <- r[!(r < 0)]
decQuant <- quantile(decVal, probs = c(0, 0.1, 0.4, 0.6, 0.9))
incQuant <- quantile(incVal, probs = c(0.1, 0.4, 0.6, 0.9, 1))
breaks <- c(decQuant, incQuant)
breaks[is.na(breaks)] <- 0
plot_seine_map(r, breaks, title)
}
```
```{r}
plot_map_delta(deltaIndicators$med[, "QA"], "QA - évolution RCP4.5 scénario médian (%)")
```
```{r} ```{r}
deltaQJXA10 <- calcDelta(QsimDrias, calcQJXAn, period = "fin", return_period = 2) plot_map_delta(deltaIndicators$min[, "QMNA5"], "QMNA5 - évolution RCP4.5 scénario minimum (%)")
lapply(deltaQJXA10, function(x) x[1:2])
``` ```
```{r}
plot_map_delta(deltaIndicators$max[, "QJXA10"], "QJXA10 - évolution RCP4.5 scénario maximum (%)")
```
Markdown is supported
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