class: center, middle, inverse, title-slide # terrainr: ## Spatial Data Access and Visualization in R ### Mike Mahoney ### 2021-09-07 --- background-image: url("ca.jpg") background-size: contain --- # About Me .pull-left[ - Mike Mahoney - Ph.D. Student at SUNY-ESF - Focusing on environmental visualization as a way to think about large-scale systems - Background in machine learning, remote sensing, and forest ecology - https://mm218.dev ] .pull-right[ <img src="https://github.com/mikemahoney218.png" width="100%" /> ] --- # Outline ### 1. What's terrainr? ### 2. Why visualize in game engines? ### 3. What's next? <br /> <br /> Slides available at https://mm218.dev/esip2021 --- class: inverse center middle # What's terrainr? ??? So, first things first, what is terrainr? --- # What's terrainr? .pull-left[ - New R package with two focuses: 1. Data access and retrieval 2. Spatial data visualization <br /> - https://github.com/ropensci/terrainr ] .pull-right[ <img src="https://github.com/ropensci/terrainr/raw/main/man/figures/logo.png" width="100%" /> ] --- class: title background-image: url("clark_USGSNAIPPlus_1_1.png") # Data access and retrieval --- ```r library(sf) campsites <- "https://opendata.arcgis.com/datasets/06993a2f20bf42d382c8ce6bd54027c4_0.geojson" campsites <- read_sf(campsites) ``` -- ```r library(ggplot2) ggplot(campsites) + geom_sf() ``` <img src="esip_files/figure-html/unnamed-chunk-4-1.png" width="100%" /> --- ```r library(terrainr) downloaded_data <- get_tiles(campsites, output_prefix = "bryce_canyon", services = c("elevation", "ortho"), resolution = 30) ``` -- ```r downloaded_data ``` ``` ## $elevation ## [1] "bryce_canyon_3DEPElevation_1_1.tif" ## ## $ortho ## [1] "bryce_canyon_USGSNAIPPlus_1_1.tif" ``` --- class: inverse title background-image: url("contour.png") # Spatial data visualization --- .pull-left[ ```r library(raster) downloaded_data$elevation |> raster() |> plot() ``` <img src="esip_files/figure-html/unnamed-chunk-7-1.png" width="100%" /> ] .pull-right[ ```r library(raster) downloaded_data$ortho |> stack() |> plotRGB(scale = 1) ``` <img src="esip_files/figure-html/unnamed-chunk-8-1.png" width="100%" /> ] --- ```r library(ggplot2) downloaded_data$elevation |> raster() |> as.data.frame(xy = TRUE) |> setNames(c("x", "y", "z")) |> ggplot(aes(x, y, fill = z)) + geom_raster() + scale_fill_gradientn( colors = (rev(terrain.colors(255))) ) + coord_sf(crs = 4326) ``` <img src="esip_files/figure-html/unnamed-chunk-9-1.png" width="100%" /> --- ```r ggplot() + geom_spatial_rgb(data = downloaded_data$ortho, aes(x, y, r = red, b = blue, g = green)) + coord_sf(crs = 4326) ``` <img src="esip_files/figure-html/unnamed-chunk-10-1.png" width="100%" /> --- ```r ggplot() + geom_spatial_rgb(data = downloaded_data$ortho, aes(x, y, r = red, b = blue, g = green)) + coord_sf(crs = 4326) + geom_sf(data = campsites, color = "red") ``` <img src="esip_files/figure-html/unnamed-chunk-11-1.png" width="100%" /> --- ```r ggplot() + geom_spatial_rgb(data = "img.png", aes(x, y, r = red, b = blue, g = green)) + coord_fixed() + labs(caption = "Artwork by @allison_horst") ``` <img src="esip_files/figure-html/unnamed-chunk-12-1.png" width="100%" /> --- class: center middle ![](usgs.jpeg) --- .pull-left[ ```r library(rayshader) library(tiff) library(raster) downloaded_data$ortho <- readTIFF(downloaded_data$ortho) downloaded_data$elevation <- raster_to_matrix( raster(downloaded_data$elevation) ) rayshader::plot_3d( downloaded_data$ortho, downloaded_data$elevation, zscale = 10, windowsize = 1200, zoom = 0.75, theta = 25, phi = 35, solid = FALSE ) ``` ] .pull-right[ <img src="BryceCanyon.png" width="100%" /> ] --- ```r downloaded_data <- get_tiles(campsites, services = c("elevation", "ortho")) ``` --- ```r downloaded_data <- get_tiles(campsites, services = c("elevation", "ortho")) mapply( merge_rasters, downloaded_data, c("elevation.tif", "orthoimagery.tif") ) ``` --- ```r downloaded_data <- get_tiles(campsites, services = c("elevation", "ortho")) mapply( merge_rasters, downloaded_data, c("elevation.tif", "orthoimagery.tif") ) make_manifest("elevation.tif", "orthoimagery.tif") ``` --- <img src="tiles.png" width="100%" /> --- <img src="vignette.png" width="100%" /> https://docs.ropensci.org/terrainr/articles/unity_instructions.html --- <img src="Capture.JPG" width="100%" /> --- # Thank you! This work was financially supported by the State University of New York via the ESF Pathways to Net Zero Carbon initiative. #### More terrainr: <svg viewBox="0 0 496 512" style="height:1em;position:relative;display:inline-block;top:.1em;" xmlns="http://www.w3.org/2000/svg"> <path d="M165.9 397.4c0 2-2.3 3.6-5.2 3.6-3.3.3-5.6-1.3-5.6-3.6 0-2 2.3-3.6 5.2-3.6 3-.3 5.6 1.3 5.6 3.6zm-31.1-4.5c-.7 2 1.3 4.3 4.3 4.9 2.6 1 5.6 0 6.2-2s-1.3-4.3-4.3-5.2c-2.6-.7-5.5.3-6.2 2.3zm44.2-1.7c-2.9.7-4.9 2.6-4.6 4.9.3 2 2.9 3.3 5.9 2.6 2.9-.7 4.9-2.6 4.6-4.6-.3-1.9-3-3.2-5.9-2.9zM244.8 8C106.1 8 0 113.3 0 252c0 110.9 69.8 205.8 169.5 239.2 12.8 2.3 17.3-5.6 17.3-12.1 0-6.2-.3-40.4-.3-61.4 0 0-70 15-84.7-29.8 0 0-11.4-29.1-27.8-36.6 0 0-22.9-15.7 1.6-15.4 0 0 24.9 2 38.6 25.8 21.9 38.6 58.6 27.5 72.9 20.9 2.3-16 8.8-27.1 16-33.7-55.9-6.2-112.3-14.3-112.3-110.5 0-27.5 7.6-41.3 23.6-58.9-2.6-6.5-11.1-33.3 2.6-67.9 20.9-6.5 69 27 69 27 20-5.6 41.5-8.5 62.8-8.5s42.8 2.9 62.8 8.5c0 0 48.1-33.6 69-27 13.7 34.7 5.2 61.4 2.6 67.9 16 17.7 25.8 31.5 25.8 58.9 0 96.5-58.9 104.2-114.8 110.5 9.2 7.9 17 22.9 17 46.4 0 33.7-.3 75.4-.3 83.6 0 6.5 4.6 14.4 17.3 12.1C428.2 457.8 496 362.9 496 252 496 113.3 383.5 8 244.8 8zM97.2 352.9c-1.3 1-1 3.3.7 5.2 1.6 1.6 3.9 2.3 5.2 1 1.3-1 1-3.3-.7-5.2-1.6-1.6-3.9-2.3-5.2-1zm-10.8-8.1c-.7 1.3.3 2.9 2.3 3.9 1.6 1 3.6.7 4.3-.7.7-1.3-.3-2.9-2.3-3.9-2-.6-3.6-.3-4.3.7zm32.4 35.6c-1.6 1.3-1 4.3 1.3 6.2 2.3 2.3 5.2 2.6 6.5 1 1.3-1.3.7-4.3-1.3-6.2-2.2-2.3-5.2-2.6-6.5-1zm-11.4-14.7c-1.6 1-1.6 3.6 0 5.9 1.6 2.3 4.3 3.3 5.6 2.3 1.6-1.3 1.6-3.9 0-6.2-1.4-2.3-4-3.3-5.6-2z"></path></svg> [ropensci/terrainr](https://github.com/ropensci/terrainr) <svg viewBox="0 0 448 512" style="height:1em;position:relative;display:inline-block;top:.1em;" xmlns="http://www.w3.org/2000/svg"> <path d="M448 360V24c0-13.3-10.7-24-24-24H96C43 0 0 43 0 96v320c0 53 43 96 96 96h328c13.3 0 24-10.7 24-24v-16c0-7.5-3.5-14.3-8.9-18.7-4.2-15.4-4.2-59.3 0-74.7 5.4-4.3 8.9-11.1 8.9-18.6zM128 134c0-3.3 2.7-6 6-6h212c3.3 0 6 2.7 6 6v20c0 3.3-2.7 6-6 6H134c-3.3 0-6-2.7-6-6v-20zm0 64c0-3.3 2.7-6 6-6h212c3.3 0 6 2.7 6 6v20c0 3.3-2.7 6-6 6H134c-3.3 0-6-2.7-6-6v-20zm253.4 250H96c-17.7 0-32-14.3-32-32 0-17.6 14.4-32 32-32h285.4c-1.9 17.1-1.9 46.9 0 64z"></path></svg> [https://docs.ropensci.org/terrainr/](https://docs.ropensci.org/terrainr/) #### Find me online: <svg viewBox="0 0 496 512" style="height:1em;position:relative;display:inline-block;top:.1em;" xmlns="http://www.w3.org/2000/svg"> <path d="M165.9 397.4c0 2-2.3 3.6-5.2 3.6-3.3.3-5.6-1.3-5.6-3.6 0-2 2.3-3.6 5.2-3.6 3-.3 5.6 1.3 5.6 3.6zm-31.1-4.5c-.7 2 1.3 4.3 4.3 4.9 2.6 1 5.6 0 6.2-2s-1.3-4.3-4.3-5.2c-2.6-.7-5.5.3-6.2 2.3zm44.2-1.7c-2.9.7-4.9 2.6-4.6 4.9.3 2 2.9 3.3 5.9 2.6 2.9-.7 4.9-2.6 4.6-4.6-.3-1.9-3-3.2-5.9-2.9zM244.8 8C106.1 8 0 113.3 0 252c0 110.9 69.8 205.8 169.5 239.2 12.8 2.3 17.3-5.6 17.3-12.1 0-6.2-.3-40.4-.3-61.4 0 0-70 15-84.7-29.8 0 0-11.4-29.1-27.8-36.6 0 0-22.9-15.7 1.6-15.4 0 0 24.9 2 38.6 25.8 21.9 38.6 58.6 27.5 72.9 20.9 2.3-16 8.8-27.1 16-33.7-55.9-6.2-112.3-14.3-112.3-110.5 0-27.5 7.6-41.3 23.6-58.9-2.6-6.5-11.1-33.3 2.6-67.9 20.9-6.5 69 27 69 27 20-5.6 41.5-8.5 62.8-8.5s42.8 2.9 62.8 8.5c0 0 48.1-33.6 69-27 13.7 34.7 5.2 61.4 2.6 67.9 16 17.7 25.8 31.5 25.8 58.9 0 96.5-58.9 104.2-114.8 110.5 9.2 7.9 17 22.9 17 46.4 0 33.7-.3 75.4-.3 83.6 0 6.5 4.6 14.4 17.3 12.1C428.2 457.8 496 362.9 496 252 496 113.3 383.5 8 244.8 8zM97.2 352.9c-1.3 1-1 3.3.7 5.2 1.6 1.6 3.9 2.3 5.2 1 1.3-1 1-3.3-.7-5.2-1.6-1.6-3.9-2.3-5.2-1zm-10.8-8.1c-.7 1.3.3 2.9 2.3 3.9 1.6 1 3.6.7 4.3-.7.7-1.3-.3-2.9-2.3-3.9-2-.6-3.6-.3-4.3.7zm32.4 35.6c-1.6 1.3-1 4.3 1.3 6.2 2.3 2.3 5.2 2.6 6.5 1 1.3-1.3.7-4.3-1.3-6.2-2.2-2.3-5.2-2.6-6.5-1zm-11.4-14.7c-1.6 1-1.6 3.6 0 5.9 1.6 2.3 4.3 3.3 5.6 2.3 1.6-1.3 1.6-3.9 0-6.2-1.4-2.3-4-3.3-5.6-2z"></path></svg> [@mikemahoney218](https://github.com/mikemahoney218/) <svg viewBox="0 0 512 512" style="height:1em;position:relative;display:inline-block;top:.1em;" xmlns="http://www.w3.org/2000/svg"> <path d="M459.37 151.716c.325 4.548.325 9.097.325 13.645 0 138.72-105.583 298.558-298.558 298.558-59.452 0-114.68-17.219-161.137-47.106 8.447.974 16.568 1.299 25.34 1.299 49.055 0 94.213-16.568 130.274-44.832-46.132-.975-84.792-31.188-98.112-72.772 6.498.974 12.995 1.624 19.818 1.624 9.421 0 18.843-1.3 27.614-3.573-48.081-9.747-84.143-51.98-84.143-102.985v-1.299c13.969 7.797 30.214 12.67 47.431 13.319-28.264-18.843-46.781-51.005-46.781-87.391 0-19.492 5.197-37.36 14.294-52.954 51.655 63.675 129.3 105.258 216.365 109.807-1.624-7.797-2.599-15.918-2.599-24.04 0-57.828 46.782-104.934 104.934-104.934 30.213 0 57.502 12.67 76.67 33.137 23.715-4.548 46.456-13.32 66.599-25.34-7.798 24.366-24.366 44.833-46.132 57.827 21.117-2.273 41.584-8.122 60.426-16.243-14.292 20.791-32.161 39.308-52.628 54.253z"></path></svg> [@mikemahoney218](https://twitter.com/mikemahoney218/) <svg viewBox="0 0 496 512" style="height:1em;position:relative;display:inline-block;top:.1em;" xmlns="http://www.w3.org/2000/svg"> <path d="M336.5 160C322 70.7 287.8 8 248 8s-74 62.7-88.5 152h177zM152 256c0 22.2 1.2 43.5 3.3 64h185.3c2.1-20.5 3.3-41.8 3.3-64s-1.2-43.5-3.3-64H155.3c-2.1 20.5-3.3 41.8-3.3 64zm324.7-96c-28.6-67.9-86.5-120.4-158-141.6 24.4 33.8 41.2 84.7 50 141.6h108zM177.2 18.4C105.8 39.6 47.8 92.1 19.3 160h108c8.7-56.9 25.5-107.8 49.9-141.6zM487.4 192H372.7c2.1 21 3.3 42.5 3.3 64s-1.2 43-3.3 64h114.6c5.5-20.5 8.6-41.8 8.6-64s-3.1-43.5-8.5-64zM120 256c0-21.5 1.2-43 3.3-64H8.6C3.2 212.5 0 233.8 0 256s3.2 43.5 8.6 64h114.6c-2-21-3.2-42.5-3.2-64zm39.5 96c14.5 89.3 48.7 152 88.5 152s74-62.7 88.5-152h-177zm159.3 141.6c71.4-21.2 129.4-73.7 158-141.6h-108c-8.8 56.9-25.6 107.8-50 141.6zM19.3 352c28.6 67.9 86.5 120.4 158 141.6-24.4-33.8-41.2-84.7-50-141.6h-108z"></path></svg> [mm218.dev](https://mm218.dev) <br /> Slides available at https://mm218.dev/esip2021