R/BV.ComByCoordL93.R

In jbfagotfede39/aquatools: Tools for freshwater managers

#' Recherche de caractéristiques administratives d'un point
#'
#' Cette fonction permet de fournir les caractéristiques géographiques d'un point donné (commune, département, région, population, etc.) à partir de la GeoAPI
#' @name BV.ComByCoordL93
#' @param X Coordonnée X à rechercher (Lambert 93)
#' @param Y Coordonnée Y à rechercher (Lambert 93)
#' @import rgeoapi
#' @import sf
#' @import tidyverse
#' @export
#' @examples
#' BV.ComByCoordL93(936722,6603736)
#' data %>% rowwise() %>% mutate(Commune = BV.ComByCoordL93(X,Y) %>% select(name) %>% as.character()) %>% ungroup()
#' data %>% rowwise() %>% mutate(codeInsee = BV.ComByCoordL93(X,Y) %>% select(codeInsee) %>% as.character()) %>% ungroup()
#' data %>% rowwise() %>% mutate(Departement = BV.ComByCoordL93(X,Y) %>% select(codeDepartement) %>% as.character()) %>% ungroup()
#' data %>% rowwise() %>% mutate(Region = BV.ComByCoordL93(X,Y) %>% select(codeRegion) %>% as.character()) %>% ungroup()
#' data %>% rowwise() %>% mutate(Population = BV.ComByCoordL93(X,Y) %>% select(population) %>% as.character()) %>% ungroup()
#' data %>% rowwise() %>% mutate(Surface = BV.ComByCoordL93(X,Y) %>% select(surface) %>% as.character()) %>% ungroup()

BV.ComByCoordL93 <- function(
  X = as.numeric(0),
  Y = as.numeric(0)
    )
{

  #### Conversion des coordonnées en Lambert WGS84 ####
  Point <- sf::st_as_sf(data.frame(X,Y), coords = c("X","Y")) %>% st_set_crs(2154) %>% st_transform(4269) %>% st_coordinates()
  
  #### Recherche en tant que telle ####
data <- rgeoapi::ComByCoord(lat = Point[,2], lon = Point[,1])

  #### Sortie des données ####
return(data)
  
} # Fin de la fonction