John Snow's cholera clusters, visualized in 3D with rayshader and ggplot2

cholera_snow.R

#theme and ggplot derived from David Kretch, his code: https://github.com/davidkretch/london_cholera_map/blob/master/london_cholera_map.R

library(HistData)
library(ggplot2)
library(ggpointdensity)
library(rayshader)

deaths = Snow.deaths
streets = Snow.streets

themeval = theme(panel.border = element_blank(), 
                 panel.grid.major = element_blank(), 
                 panel.grid.minor = element_blank(), 
                 axis.line = element_blank(), 
                 axis.ticks = element_blank(),
                 axis.text.x = element_blank(), 
                 axis.text.y = element_blank(), 
                 axis.title.x = element_blank(), 
                 axis.title.y = element_blank(),
                 legend.key = element_blank(),
                 plot.margin = unit(c(0.5, 0, 0, 0), "cm"))
js = ggplot() + 
  geom_path(data=streets,aes(x=x, y=y, group = street), color="grey50") +
  geom_pointdensity(data = deaths,aes(x=x, y=y), size=1,adjust=0.1) +
  coord_fixed() + 
  scale_color_viridis_c() +
  theme_bw() +
  themeval 

js2 = ggplot() + 
  geom_path(data = streets, aes(x=x, y=y, group = street), color="white") +
  geom_pointdensity(data = deaths, aes(x=x, y=y), size=1, adjust=0.1) +
  coord_fixed() +
  scale_color_viridis_c() +
  theme_bw() +
  themeval 

ggheight = plot_gg(list(js,js2), multicore = TRUE, raytrace=TRUE,
                   height_aes = "color", shadow_intensity = 0.3,
                   width=8,height=7, soliddepth = -100, save_height_matrix = TRUE,
                   background = "#f5e9dc", shadowcolor= "#4f463c",windowsize=c(1000,1000))

#Add label
render_label(ggheight, "Water Pump", x=1110,y=1020,z=390, textsize = 2)

#Generate animation
for(i in 1:360) {
  render_camera(phi=20,theta=45+i,fov=70,zoom=0.25)
  render_depth(focus=0.77, 
               title_text = "John Snow's London cholera map, # of nearest neighbors", 
               title_size = 35,
               filename = glue::glue("snow{i}"))
}

Getting errors when running script:

Error in if (whichtype %in% c("text", "line")) { :
argument is of length zero
(from ggheight line, followed by)
Error in render_label(ggheight, "Water Pump", x = 1110, y = 1020, z = 390, :
No rgl window currently open.
Error in render_camera(phi = 20, theta = 45 + i, fov = 70, zoom = 0.25) :
No rgl window currently open.

session info:

R version 4.0.0 (2020-04-24)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 10 x64 (build 18362)

Matrix products: default

locale:
[1] LC_COLLATE=English_United States.1252
[2] LC_CTYPE=English_United States.1252
[3] LC_MONETARY=English_United States.1252
[4] LC_NUMERIC=C
[5] LC_TIME=English_United States.1252

attached base packages:
[1] stats graphics grDevices utils datasets
[6] methods base

other attached packages:
[1] rayshader_0.15.1 ggpointdensity_0.1.0
[3] ggplot2_3.3.1 HistData_0.8-6

loaded via a namespace (and not attached):
[1] Rcpp_1.0.4.6 later_1.0.0
[3] pillar_1.4.3 compiler_4.0.0
[5] prettyunits_1.1.1 iterators_1.0.12
[7] tools_4.0.0 progress_1.2.2
[9] digest_0.6.25 rpart_4.1-15
[11] jsonlite_1.6.1 lifecycle_0.2.0
[13] tibble_3.0.1 gtable_0.3.0
[15] viridisLite_0.3.0 pkgconfig_2.0.3
[17] rlang_0.4.5 foreach_1.5.0
[19] shiny_1.4.0.2 rstudioapi_0.11
[21] crosstalk_1.1.0.1 parallel_4.0.0
[23] xfun_0.13 fastmap_1.0.1
[25] knitr_1.28 withr_2.2.0
[27] dplyr_0.8.5 htmlwidgets_1.5.1
[29] vctrs_0.2.4 hms_0.5.3
[31] webshot_0.5.2 manipulateWidget_0.10.1
[33] grid_4.0.0 tidyselect_1.0.0
[35] glue_1.4.0 R6_2.4.1
[37] rgl_0.100.54 purrr_0.3.4
[39] farver_2.0.3 magrittr_1.5
[41] promises_1.1.0 htmltools_0.4.0
[43] scales_1.1.0 codetools_0.2-16
[45] ellipsis_0.3.0 assertthat_0.2.1
[47] xtable_1.8-4 mime_0.9
[49] colorspace_1.4-1 httpuv_1.5.3.1
[51] labeling_0.3 miniUI_0.1.1.1
[53] munsell_0.5.0 doParallel_1.0.15
[55] crayon_1.3.4

For anyone else chasing down a way to reproduce this, a few tips:

1. Overall, easier to do on mac
2. Start with most recent development version (on github) of rayshader
3. lots of dependencies make Windows more challenging, like freetext, then reinstalling rgl
4. Easier to make a mp4 of this with `render_movie("john_snow.mp4", phi=20, fov=70, zoom=0.30, title_text = "John Snow's London cholera map, # of nearest neighbors", title_size = 35). This will not include the bokeh focus effect, however.
5. To make a gif, follow the directions above, which will produce 360 png files, from snow1.png to snow360.png. Unfortunately, these will be ordered in alphanumeric order, so that snow1 is followed by snow10, then snow100, snow101, etc. Not the order that you want. You can put these files into a folder, and extract the vector of files with (fs package)
   files <- fs::dir_ls("folder")
6. Then extract the numeric portion of each file name with (stringr or tidyverse package)
   num <- as.numeric(str_extract(files, "\\d+")
7. then combine these two vectors into a tibble with (tibble package)
   frame <- tibble(files, num) %>% arrange(num)
8. Then turn this into a gif with the gifski package -
   gifski::gifski(frame$files, gif_file = "snow.gif", delay = 0.04, loop = TRUE, progress = TRUE)