qualite_esu <- get_qualite_eau(code="04197700", 
                               milieu_preleve="ESU", 
                               date_debut_maj = "2023-01-01", 
                               date_fin_maj = "2024-01-01",
                               date_debut_prelevement="2003-01-08", 
                               date_fin_prelevement="2021-01-25")
#> Page 1 récupérée : 5000 lignes
#> Page 2 récupérée : 5000 lignes
#> Page 3 récupérée : 5000 lignes
#> Page 4 récupérée : 2190 lignes
#> Redémarrage avec date = 2021-01-25
#> Page 1 récupérée : 573 lignes


qualite_eso <- get_qualite_eau(code="BSS001EUDY", 
                               milieu_preleve="ESO", 
                               date_debut_maj = "2023-01-01", 
                               date_fin_maj = "2024-01-01",
                               date_debut_prelevement="2007-06-04", 
                               date_fin_prelevement="2007-10-04")
#> Page 1 récupérée : 770 lignes

library(dplyr)
#> 
#> Attachement du package : 'dplyr'
#> Les objets suivants sont masqués depuis 'package:stats':
#> 
#>     filter, lag
#> Les objets suivants sont masqués depuis 'package:base':
#> 
#>     intersect, setdiff, setequal, union
# Coordonnées approximatives des villes (longitude, latitude en WGS84)
coords_wgs84 <- data.frame(
  ville = c("Rennes", "Redon", "Ploermel"),
  lon = c(-1.6794, -2.0842, -2.3975),
  lat = c(48.1147, 47.6536, 47.9333)
)

# Conversion en objet sf (points)
points_sf <- sf::st_as_sf(coords_wgs84, coords = c("lon", "lat"), crs = 4326)

# Création du polygone (triangle)
triangle <- sf::st_union(points_sf) %>% 
  sf::st_combine() %>% 
  sf::st_cast("POLYGON")

# Création de l'objet sf final
triangle_sf <- sf::st_sf(geometry = sf::st_sfc(triangle))

stations_esu<-get_stations_hubeau(sf_object=triangle_sf,
                    milieu_preleve='ESU')
#> No encoding supplied: defaulting to UTF-8.

stations_eso<-get_stations_hubeau(sf_object=triangle_sf,
                    milieu_preleve='ESO')

library(mapview)
mapview(triangle_sf) + mapview(stations_esu) + mapview(stations_eso, col.regions="red")