Leaflet for R

load libraries

library(virtualspecies)
library(ggplot2)
library(plyr)
library(reshape2)
library(maptools)
library(vegan)
library(dplyr)
library(leaflet)
library(RColorBrewer)

load bioclim

worldclim <- getData("worldclim", var = "bio", res = 10)

map a raster in leaflet

r <- worldclim[[2]]

pal <- colorNumeric("RdYlBu", values(r), na.color = "transparent")

leaflet() %>% addTiles() %>%
  addRasterImage(r, colors = pal, opacity = 0.8) %>%
  addLegend(pal = pal, values = values(r), title = "worldclim[[2]]")

## Warning in rgdal::rawTransform(projfrom, projto, nrow(xy), xy[, 1], xy[, :
## 54 projected point(s) not finite

make a map with different baselayers

r <- worldclim[[1]]
r2 <- worldclim[[2]]

pal <- colorNumeric("RdYlBu", values(r), na.color = "transparent")
pal2 <- colorNumeric("Blues", values(r2), na.color = "transparent")

mapbox <- "http://api.tiles.mapbox.com/v4/mapbox.outdoors/{z}/{x}/{y}.png?access_token=pk.eyJ1Ijoiam9zaHBlcHBlciIsImEiOiJuTWdrY2k4In0.HCCXtgU04scrTB_-ON4kjA"

leaflet() %>% addTiles() %>%
  addRasterImage(r, colors = pal, opacity = 0.6, group="one") %>%
  addRasterImage(r2, colors = pal2, opacity = 0.6, group="two") %>%
  addLegend(pal = pal, values = values(r), title = "worldclim[[2]]", position = "bottomleft") %>%
  addProviderTiles("CartoDB.Positron", group="Simple") %>%
  addTiles(urlTemplate = mapbox, group="Outdoors") %>%
  addProviderTiles("Esri.WorldImagery", group="Satellite") %>%
  addLayersControl(
    baseGroups = c("Satellite", "Simple", "Outdoors"),
    overlayGroups = c("one","two")
  ) %>%
  hideGroup("two")

## Warning in rgdal::rawTransform(projfrom, projto, nrow(xy), xy[, 1], xy[, :
## 54 projected point(s) not finite

## Warning in rgdal::rawTransform(projfrom, projto, nrow(xy), xy[, 1], xy[, :
## 54 projected point(s) not finite