In pokyah/defHydWal: Déficit hydrique en Wallonia : cartes et graphiques

library(dplyr)
library(ggsn) #from github !!
devtools::load_all()

## show grouped code output instead of single lines
knitr::opts_chunk$set(collapse = FALSE)
knitr::knit_hooks$set(document = function(x){
  gsub("```\n*```r*\n*", "", x)
})

knitr::opts_chunk$set(warning = FALSE, echo = FALSE, message = FALSE, fig.width = 8, fig.height = 6, results = "asis")

data(plu.obs)
data(cum.norm)
data(st_info)
# checking if the input date goes beyond 30-06-2018 (in case TRUE, needs API call to get the data not available as .Rda object)
if (params$dto > as.Date("2018-06-30")) {
  plu.obs <- dhw_add_days(dto = params$dto,
    data = plu.obs)
} 
# filtering the data to only keep relevant stations and the period of interest
plu.obs = dhw_filter(data = plu.obs,
  dfrom = params$dfrom,
  dto = params$dto)

cum.norm = dhw_filter(data = cum.norm,
  dfrom = params$dfrom,
  dto = params$dto)

cum.obs = dhw_obs.cums(plu.obs)
cum = dhw_join_cums(obs = cum.obs, norm = cum.norm)

Comparaison avec les données historiques : Situation à Ernage

data(ernage.all_y)
all_y.plot = dhw_plot(ernage.all_y, int = params$int)

data(ernage.two_y)
two_y.plot = dhw_plot(ernage.two_y, int = params$int)

all_y.plot
two_y.plot

Cartographie pour l'ensemble de la Région Wallonne

dhw_summary = dhw_summarize(cum, st_info)

data(grid.grid)
dh.spatial = dhw_idw(data = dhw_as.spatial(dhw_summary), param = "defExHyd", grid = grid.grid)
ip.spatial = dhw_idw(data = dhw_as.spatial(dhw_summary), param = "ind_plu", grid = grid.grid)

Indice pluviométrique

data(wallonia.sf)
ip.map = dhw_map(gridded.data.df = dhw_2_ggmap(ip.spatial),
  boundaries.sf = wallonia.sf,
  title.chr = paste0("Carte de l'indice pluviométrique calculé pour la période du ",
   format(as.Date(params$dfrom, "%Y-%m-%d"), "%d %b %Y"), " au ", format(as.Date(params$dto, "%Y-%m-%d"), "%d %b %Y")),
  legend.chr = "",
  target.chr = "var1.pred",
  nb_classes.num = params$ip_cl,
  reverse_pal.bool = FALSE,
  pretty_breaks.bool = TRUE
)
ip.map

Déficit hydrique

dh.map = dhw_map(gridded.data.df = dhw_2_ggmap(dh.spatial),
  boundaries.sf = wallonia.sf,
  title.chr = paste0("Carte du déficit hydrique calculé pour la période du ",
     format(as.Date(params$dfrom, "%Y-%m-%d"), "%d %b %Y"), " au ", format(as.Date(params$dto, "%Y-%m-%d"), "%d %b %Y")),
  legend.chr = "mm",
  target.chr = "var1.pred",
  nb_classes.num = params$dh_cl,
  reverse_pal.bool = FALSE,
  pretty_breaks.bool = TRUE
)
dh.map