Bird habitat use analysis

library(tidyverse)
library(mapview)
library(sf)
library(tmap)
library(RColorBrewer)
load("Exploratory_maps.rda")
mapviewOptions(basemaps =c("Esri.WorldImagery",
                             "OpenStreetMap",
                             "Esri.WorldStreetMap"), homebutton = FALSE)
pal<- brewer.pal(12,"Paired")

DEM

library(elevatr)
library(raster)
pal1<-terrain.colors(100)
pal2<-gray(0:100 / 100)
dem<-elevatr::get_elev_raster(birds, z=12)
dslope<-terrain(dem, opt='slope', unit='radians')
daspect<-terrain(dem, opt='aspect',unit='radians')
dhill<-hillShade(dslope,daspect, angle=10)

birds$elevation<-raster::extract(dem,birds)

library(leaflet.extras)
map<-mapview(dem,col.regions=pal1, legend=FALSE, na.color="transparent") + 
  mapview(dhill,col.regions=pal2,na.color="transparent", legend=FALSE) 
map <-map +mapview(birds)
map@map <- map@map  %>%
    addFullscreenControl() 
map

save(list=ls(),file="DEM.rda")