R/plotETPMean.R

In paul-carteron/cadastreAnalysis: Analysis of cadastral parcels in France

#' plotETPMean
#'
#' @param ETPMeanRaster Raster des ETP moyennes obtenues avec la fonction importETPMean
#' @param zoneEtude Objet sf contenant le shape de la zoneEtude
#' @param mapBackground Choix du fond de carte : "OpenStreetMap" , "Scan25" , "Ortho"

#' @return Renvoi une carte interactive de la temperature moyennes autour de notre zoneEtude.
#' @export
#'
#' @import leaflet dplyr
#' @importFrom raster values
#'
plotETPMean <- function(ETPMeanRaster, zoneEtude, mapBackground = "OpenStreetMap"){

   if (mapBackground == "Ortho"){
      background = "GeoportailFrance.orthos"
   }else if (mapBackground == "Scan25"){
      background = "GeoportailFrance.ignMaps"
   }else if (mapBackground == "OpenStreetMap"){
      background = "OpenStreetMap.France"
   }else{
      stop("L'argument mapBackground n'est pas rempli. Choisir : \"OpenStreetMap\" , \"Scan25\" ou \"Ortho")
   }

   palette = terrain.colors(length(unique(values(ETPMeanRaster))))
   pal = colorNumeric(palette, domain = values(ETPMeanRaster), na.color = "black")

   res = leaflet() %>%
      addProviderTiles(background) %>%
      addRasterImage(ETPMeanRaster,colors = pal, opacity = 0.6) %>%
      addLegend(
         pal = pal,
         values = values(ETPMeanRaster),
         title = "ETP [mm]")%>%
      addPolylines(data = zoneEtude,
                   opacity = 1,
                   stroke = TRUE,
                   weight = 1,
                   fill = FALSE,
                   color = "black")
   print(res)
   return(res)
}