Pull elevation data with the elevatr package and rayshade with rayshader package

elevatr_rayshade.R

library(raster)
library(elevatr)
library(rayshader)
library(reshape2)
library(ggplot2)

data(lake)
x <- get_elev_raster(lake, z = 12)
x_mat <- as.matrix(x)
x_rsh <- rayshader::rayshade(x_mat, 
   anglebreaks = seq(55,65,10), 
   sunangle = 315, 
   maxsearch = 100)

ggplot(data = reshape2::melt(x_rsh)) + 
  geom_tile(aes(x = Var1, y = Var2, fill = value)) +
  scale_fill_viridis_c() + 
  theme_void() + 
  theme(legend.position = "none")