In AnthonyDEBUR/tools4DCE: Tools for rivers water quality monitoring

library(tools4DCE)

Ce document explique comment télécharger des indices hydrobiologiques depuis Hubeau, convertir IBD et IBMR en EQR, affecter une classe de qualité aux résultats selon la typologie des stations de mesures.

Chargement de donnees hydrobiologiques depuis Hubeau

A partir de la plateforme Hubeau, on télécharge les données hydrobiologiques disponibles sur la station 04207400.

donnees<-import_hubeau_indices_hbio(liste_stations = c("04207400"))
str(donnees)

Conversion des notes IBD et IBMR en EQR

Le package tools4DCE fournit une table avec les paramètres nécessaires au calcul des EQR selon la typologie de cours d'eau, conformément à l'arrêté ministériel du 25 janvier 2010 modifié.

str(base_ref_eqr)

Calcul des EQR

Le package contient une fonction pour calculer les EQR.

donnees_IBD<-donnees%>%subset(CdParametre=="5856" & !is.na(RsAna))
donnees_IBD$EQR<-sapply(donnees_IBD$RsAna, function(x){calcule_EQR_hbio(x, CdParametre='5856', typologie="P12-A")})

Graphique des résultats

Le package permet de mettre facilement les résultats sous forme graphique

Création d'un seuil qualité :

seuil_IBD<-makeSeuils(CdParametre = "5856", type_seuil = "DCE", specificites = "<10000 km2")
graphDCE_points(donnees_IBD, col_valeurs = "EQR", seuils=seuil_IBD, ymini=0, ymaxi=1) + ylab("EQR")

Calcul de la classe de qualité correspondant à chaque résultat

Il permet également d'affecter une classe de qualité à chaque résultat.

Création d'un seuil qualité :

donnees_IBD$CLASSE<-sapply(donnees_IBD$EQR, function(x){affecte_une_classe(x, seuil_IBD)})

Chargement des listes floristiques depus Hubeau et recalcule de l'IBD via le SEEE

Import des listes floristiques

liste_flo<-import_hubeau_liste_hbio(liste_stations = "04195400", indice="dia")
str(head(liste_flo, 3))

Mise en forme des données

liste_flo$CODE_OPERATION <-paste0(liste_flo$code_station_hydrobio, "*", liste_flo$date_prelevement)
liste_flo$CODE_STATION <- liste_flo$code_station_hydrobio
liste_flo$DATE <- format(liste_flo$date_prelevement, "%d/%m/%Y")
liste_flo$CODE_TAXON <- liste_flo$code_appel_taxon
liste_flo$RESULTAT <- liste_flo$resultat_taxon
liste_flo <- liste_flo %>% select(CODE_OPERATION, CODE_STATION, DATE, CODE_TAXON, RESULTAT)
liste_flo<-liste_flo%>%subset(!is.na(CODE_TAXON))