Objectif du tuto

Il s’agit de montrer comment accéder simplement à une masse de “lots de données de référence” disponibles via le flux WFS du Sandre avec le logiciel R et de la fonction wfs_sandre() du package tod.

Activation des `packages` :

library(tod)
library(tidyverse)
library(mapview)

Les données Sandre

Pour s’y retrouver parmi les données diffusées sur le portail du Sandre, il est utile de se reporter à cette page. Ces référentiels géographiques concernent les thèmes suivants :

Url	Thème
https://services.sandre.eaufrance.fr/geo/eth	Hydrographique bd_carthage
https://services.sandre.eaufrance.fr/geo/mdo	Masses d'eau hors meso/mesu
https://services.sandre.eaufrance.fr/geo/MasseDEau_VRAP2010	meso/mesu wise
https://services.sandre.eaufrance.fr/geo/MasseDEau_VEDL2013	meso/mesu sandre
https://services.sandre.eaufrance.fr/geo/MasseDEau_VRAP2016	meso/mesu wise
https://services.sandre.eaufrance.fr/geo/obs	Obstacles à l'ecoulement
https://services.sandre.eaufrance.fr/geo/odp	Ouvrages de dépollution
https://services.sandre.eaufrance.fr/geo/saq	Hydrogéologie
https://services.sandre.eaufrance.fr/geo/sel	Lieux de surveillance
https://services.sandre.eaufrance.fr/geo/stq	Stations de mesure de la qualité des eaux de surface
https://services.sandre.eaufrance.fr/geo/zagri	Zone agricole
https://services.sandre.eaufrance.fr/geo/zppn	Zones de Protection du Patrimoine Naturel
https://services.sandre.eaufrance.fr/geo/zrpe	Zones de Régulations pour l'Environnement
https://services.sandre.eaufrance.fr/geo/zgr	Zones de Gestion de la Ressource
https://services.sandre.eaufrance.fr/geo/hyd	Hydrométrie
https://services.sandre.eaufrance.fr/geo/com	Administratif
https://services.sandre.eaufrance.fr/geo/vic	Vigilance crue

Les possibilités de sélection géographique sont limitées. Elles s’opèrent en rajoutant un suffixe au nom de la couche. Par exemple, pour les obstacles à l’écoulement, qu’on peut trouver par le moteur de recherche du site :

https://www.sandre.eaufrance.fr/Rechercher-un-jeu-de-donnees?keyword=%22Obstacles%22

On peut lire sur la fiche de métadonnées que le nom de la couche est ObstEcoul, ce qui renvoie les données métropole et DROM.

Si l’on ajoute un suffixe on peut ne charger que :

la métropole (ObstEcoul_FXX)
la Guadeloupe (ObstEcoul_GLP)
etc.

Configuration machine

Les étapes de configuration de la machine et d’installation des packages sont considérées comme déjà réalisées. Si ce n’est pas le cas se reporter à ces explications.

Banque Nationale des Prélèvements d’eau (BNPE)

On peut accéder aux couches relatives aux ouvrages et aux points de prélèvement (NB les fichiers sont lourds donc le chargement prend un certain temps).

bnpe_ouvrages <- wfs_sandre(url_wfs = "https://services.sandre.eaufrance.fr/geo/prl",
                            couche = "OuvragePrel")

bnpe_points_prel <- wfs_sandre(url_wfs = "https://services.sandre.eaufrance.fr/geo/prl",
                               couche = "PointPrelevement")

Pour ne charger que les données de la Guadeloupe :

bnpe_ouvrages_glp <- wfs_sandre(url_wfs = "https://services.sandre.eaufrance.fr/geo/prl",
                            couche = "OuvragePrel_GLP")

## Reading layer `OuvragePrel_GLP' from data source 
##   `https://services.sandre.eaufrance.fr/geo/prl?service=wfs&version=2.0.0&request=GetCapabilities' 
##   using driver `WFS'
## Simple feature collection with 93 features and 27 fields
## Geometry type: POINT
## Dimension:     XY
## Bounding box:  xmin: -61.7796 ymin: 15.90608 xmax: -61.22396 ymax: 16.38958
## Geodetic CRS:  WGS 84

bnpe_points_prel_glp <- wfs_sandre(url_wfs = "https://services.sandre.eaufrance.fr/geo/prl",
                               couche = "PointPrelevement_GLP")

## Reading layer `PointPrelevement_GLP' from data source 
##   `https://services.sandre.eaufrance.fr/geo/prl?service=wfs&version=2.0.0&request=GetCapabilities' 
##   using driver `WFS'
## Simple feature collection with 117 features and 26 fields
## Geometry type: POINT
## Dimension:     XY
## Bounding box:  xmin: -61.7796 ymin: 15.90608 xmax: -61.22396 ymax: 16.38958
## Geodetic CRS:  WGS 84

Référentiel des obstacles à l’écoulement (ROE)

obstacles_ecoulement <- wfs_sandre(url_wfs = "https://services.sandre.eaufrance.fr/geo/obs",
                                   couche = "ObstEcoul")

## Reading layer `ObstEcoul' from data source 
##   `https://services.sandre.eaufrance.fr/geo/obs?service=wfs&version=2.0.0&request=GetCapabilities' 
##   using driver `WFS'
## Simple feature collection with 112520 features and 47 fields
## Geometry type: POINT
## Dimension:     XY
## Bounding box:  xmin: -61.79655 ymin: -21.38572 xmax: 55.76965 ymax: 51.05526
## Geodetic CRS:  WGS 84

Hydroécorégions

Chargement :

her2 <-
  wfs_sandre(url_wfs = "https://services.sandre.eaufrance.fr/geo/mdo?",
             couche = "Hydroecoregion2")

## Reading layer `Hydroecoregion2' from data source 
##   `https://services.sandre.eaufrance.fr/geo/mdo?service=wfs&version=2.0.0&request=GetCapabilities' 
##   using driver `WFS'
## Simple feature collection with 114 features and 5 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: -4.79495 ymin: 41.36882 xmax: 9.560037 ymax: 51.08965
## Geodetic CRS:  WGS 84

Affichage :

ggplot(her2) +
  geom_sf()

Les couches de la base TOPAGE

Les bassins hydrographiques :

bassins_hydro <-
  wfs_sandre(url_wfs = "https://services.sandre.eaufrance.fr/geo/topage?",
             couche = "BassinHydrographique_FXX")

## Reading layer `BassinHydrographique_FXX' from data source 
##   `https://services.sandre.eaufrance.fr/geo/topage?service=wfs&version=2.0.0&request=GetCapabilities' 
##   using driver `WFS'
## Simple feature collection with 7 features and 6 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: -5.141592 ymin: 41.36462 xmax: 9.559987 ymax: 51.08945
## Geodetic CRS:  WGS 84

ggplot(bassins_hydro) +
  geom_sf()

D’autres couches (certaines, dépassant 100Mo, sont longues à télécharger si votre connexion est faiblarde).

bv_topo <-
  wfs_sandre(url_wfs = "https://services.sandre.eaufrance.fr/geo/topage?",
             couche = "BassinVersantTopographique_FXX")

## Reading layer `BassinVersantTopographique_FXX' from data source 
##   `https://services.sandre.eaufrance.fr/geo/topage?service=wfs&version=2.0.0&request=GetCapabilities' 
##   using driver `WFS'
## Simple feature collection with 6190 features and 17 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: -5.141592 ymin: 41.36462 xmax: 9.559987 ymax: 51.08945
## Geodetic CRS:  WGS 84

cours_deau <-
  wfs_sandre(url_wfs = "https://services.sandre.eaufrance.fr/geo/topage?",
             couche = "CoursEau")

## Reading layer `CoursEau' from data source 
##   `https://services.sandre.eaufrance.fr/geo/topage?service=wfs&version=2.0.0&request=GetCapabilities' 
##   using driver `WFS'
## Simple feature collection with 146233 features and 10 fields
## Geometry type: MULTILINESTRING
## Dimension:     XY
## Bounding box:  xmin: -61.80598 ymin: -12.99956 xmax: 45.23655 ymax: 51.07512
## Geodetic CRS:  WGS 84

plans_deau <-
  wfs_sandre(url_wfs = "https://services.sandre.eaufrance.fr/geo/topage?",
             couche = "PlanEau")

## Reading layer `PlanEau' from data source 
##   `https://services.sandre.eaufrance.fr/geo/topage?service=wfs&version=2.0.0&request=GetCapabilities' 
##   using driver `WFS'
## Simple feature collection with 34455 features and 13 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: -61.76912 ymin: -12.79715 xmax: 45.29148 ymax: 51.06197
## Geodetic CRS:  WGS 84

BD CarTHAgE

La BD CarTHAgE (Base de Données sur la CARtographie THématique des AGences de l’eau et du ministère chargé de l’environnement) est antérieure à Topage mais contient certaines informations complémentaires :

Par exemple pour le découpage hydrographique en “secteurs hydrographiques”, plus fin que les “bassins hydrographiques” de Topage ci-dessus :

bv <-
  tod::wfs_sandre(url_wfs = "https://services.sandre.eaufrance.fr/geo/sandre",
                  couche = "SecteurHydro_FXX")

## Reading layer `SecteurHydro_FXX' from data source 
##   `https://services.sandre.eaufrance.fr/geo/sandre?service=wfs&version=2.0.0&request=GetCapabilities' 
##   using driver `WFS'
## Simple feature collection with 187 features and 6 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: -5.141592 ymin: 41.36462 xmax: 9.559987 ymax: 51.08945
## Geodetic CRS:  WGS 84

bv %>%
  ggplot() +
    geom_sf()

Les stations et sites hydrométriques

Ce sont les stations où sont effectuées les mesures de débit consultables via la banque HYDRO :

stations_hydro <-
  wfs_sandre(url_wfs = "https://services.sandre.eaufrance.fr/geo/hyd?",
             couche = "StationHydro")

## Reading layer `StationHydro' from data source 
##   `https://services.sandre.eaufrance.fr/geo/hyd?service=wfs&version=2.0.0&request=GetCapabilities' 
##   using driver `WFS'
## Simple feature collection with 5440 features and 15 fields
## Geometry type: POINT
## Dimension:     XY
## Bounding box:  xmin: -61.79058 ymin: -21.38101 xmax: 55.72198 ymax: 50.95467
## Geodetic CRS:  WGS 84

sites_hydro <-
  wfs_sandre(url_wfs = "https://services.sandre.eaufrance.fr/geo/hyd?",
             couche = "SiteHydro")

## Reading layer `SiteHydro' from data source 
##   `https://services.sandre.eaufrance.fr/geo/hyd?service=wfs&version=2.0.0&request=GetCapabilities' 
##   using driver `WFS'
## Simple feature collection with 3396 features and 11 fields
## Geometry type: POINT
## Dimension:     XY
## Bounding box:  xmin: -4.70483 ymin: 41.55289 xmax: 9.517331 ymax: 50.90007
## Geodetic CRS:  WGS 84

Entités hydrogéologiques

Exemple de Mayotte :

entites_hydrogeol_myt <-
  wfs_sandre(url_wfs = "http://services.sandre.eaufrance.fr/geo/saq?",
             couche = "EntiteHydroGeol_MYT")

## Reading layer `EntiteHydroGeol_MYT' from data source 
##   `http://services.sandre.eaufrance.fr/geo/saq?service=wfs&version=2.0.0&request=GetCapabilities' 
##   using driver `WFS'
## Simple feature collection with 19 features and 12 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: 45.01818 ymin: -13.00521 xmax: 45.30006 ymax: -12.63664
## Geodetic CRS:  WGS 84

mapview::mapview(entites_hydrogeol_myt,
                 alpha.regions = 0.01)

Contrats de milieu (contrats de rivière)

cdm <-
  wfs_sandre(url_wfs = "https://services.sandre.eaufrance.fr/geo/zpl?",
             couche = "ContratMilieu_FXX")

## Reading layer `ContratMilieu_FXX' from data source 
##   `https://services.sandre.eaufrance.fr/geo/zpl?service=wfs&version=2.0.0&request=GetCapabilities' 
##   using driver `WFS'
## Simple feature collection with 293 features and 10 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: -61.7655 ymin: 14.36859 xmax: 9.534851 ymax: 50.96395
## Geodetic CRS:  WGS 84

Package {tod} - Flux WFS SANDRE

Pascal Irz

2024-08-31

Objectif du tuto

Activation des `packages` :

Les données Sandre

Configuration machine

Banque Nationale des Prélèvements d’eau (BNPE)

Référentiel des obstacles à l’écoulement (ROE)

Hydroécorégions

Les couches de la base TOPAGE

BD CarTHAgE

Les stations et sites hydrométriques

Entités hydrogéologiques

Contrats de milieu (contrats de rivière)

Package {tod} - Flux WFS SANDRE

Pascal Irz

2024-08-31

Objectif du tuto

Activation des packages :

Les données Sandre

Configuration machine

Banque Nationale des Prélèvements d’eau (BNPE)

Référentiel des obstacles à l’écoulement (ROE)

Hydroécorégions

Les couches de la base TOPAGE

BD CarTHAgE

Les stations et sites hydrométriques

Entités hydrogéologiques

Contrats de milieu (contrats de rivière)

Activation des `packages` :