## ---- include = FALSE---------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ---- eval = FALSE------------------------------------------------------------
#  zion <- tmaptools::geocode_OSM("Zion National Park")$coords

## ---- eval = FALSE------------------------------------------------------------
#  library(terrainr)
#  library(sf)
#  library(magrittr)
#  
#  zion <- data.frame(x = zion[["x"]], y = zion[["y"]]) %>%
#    st_as_sf(coords = c("x", "y"), crs = 4326) %>%
#    set_bbox_side_length(8000)

## ---- eval = FALSE------------------------------------------------------------
#  merged_tiles <- zion %>%
#    get_tiles(services = c("elevation", "ortho")) %>%
#    lapply(merge_rasters)

## ----eval = FALSE-------------------------------------------------------------
#  make_unity(
#    project = "zion",
#    heightmap = merged_tiles$elevation,
#    overlay = merged_tiles$ortho
#  )

## ---- echo = FALSE------------------------------------------------------------
knitr::include_graphics("terrain_surface.jpg")