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

## ----setup--------------------------------------------------------------------
library(terrainr)

## -----------------------------------------------------------------------------
mt_elbert_points <- data.frame(
  lat = runif(100, min = 39.11144, max = 39.12416),
  lng = runif(100, min = -106.4534, max = -106.437)
)

## -----------------------------------------------------------------------------
mt_elbert_points <- sf::st_as_sf(mt_elbert_points,
                                 coords = c("lng", "lat"))
mt_elbert_points <- sf::st_set_crs(mt_elbert_points, 4326)

## ----eval = FALSE-------------------------------------------------------------
#  library(progressr)
#  handlers("progress")
#  with_progress(
#    output_files <- get_tiles(mt_elbert_points,
#                              output_prefix = tempfile(),
#                              services = c("elevation", "ortho"))
#    )

## ----eval = FALSE-------------------------------------------------------------
#  output_files

## ----echo = FALSE-------------------------------------------------------------
output_files <- list(
  elevation = "/tmp/RtmphTFQvZ/file65e5d859e628_3DEPElevation_1_1.tif",
  ortho = "/tmp/RtmphTFQvZ/file65e5d859e628_USGSNAIPPlus_1_1.tif"
)
output_files

## ----eval=FALSE---------------------------------------------------------------
#  raster::plot(raster::raster(output_files[[1]]))

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

## ----eval=FALSE---------------------------------------------------------------
#  raster::plotRGB(raster::brick(output_files[[2]]), scale = 1)

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

## ----eval=FALSE---------------------------------------------------------------
#  library(ggplot2)
#  
#  elevation_raster <- raster::raster(output_files[[1]])
#  elevation_df <- as.data.frame(elevation_raster, xy = TRUE)
#  elevation_df <- setNames(elevation_df, c("x", "y", "elevation"))
#  
#  ggplot() +
#    geom_raster(data = elevation_df, aes(x = x, y = y, fill = elevation)) +
#    scale_fill_distiller(palette = "BrBG") +
#    coord_sf(crs = 4326)

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

## ----eval = FALSE-------------------------------------------------------------
#  ortho_raster <- raster::stack(output_files[[2]])
#  ortho_df <- as.data.frame(ortho_raster, xy = TRUE)
#  ortho_df <- setNames(ortho_df, c("x", "y", "red", "green", "blue", "alpha"))
#  
#  ggplot() +
#    geom_spatial_rgb(data = ortho_df,
#                     # Required aesthetics r/g/b specify color bands:
#                     aes(x = x, y = y, r = red, g = green, b = blue)) +
#    coord_sf(crs = 4326)

## -----------------------------------------------------------------------------
knitr::include_graphics("ortho_ggplot.jpg")

## ----eval = FALSE-------------------------------------------------------------
#  ggplot() +
#    geom_spatial_rgb(data = ortho_raster,
#                     aes(x = x, y = y, r = red, g = green, b = blue)) +
#    coord_sf(crs = 4326)

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

## ----eval = FALSE-------------------------------------------------------------
#  ggplot() +
#    geom_spatial_rgb(data = output_files[[2]],
#                     aes(x = x, y = y, r = red, g = green, b = blue)) +
#    coord_sf(crs = 4326)

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

## ----eval = FALSE-------------------------------------------------------------
#  ggplot() +
#    geom_spatial_rgb(data = output_files[[2]],
#                     aes(x = x, y = y, r = red, g = green, b = blue)) +
#    geom_sf(data = mt_elbert_points)

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

## ---- eval = FALSE------------------------------------------------------------
#  mt_elbert_overlay <- vector_to_overlay(mt_elbert_points,
#                                         output_files[[2]],
#                                         size = 15,
#                                         color = "red")
#  knitr::include_graphics(mt_elbert_overlay)

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

## ---- eval=FALSE--------------------------------------------------------------
#  ortho_with_points <- combine_overlays(
#    # Overlays are stacked in order, with the first file specified on the bottom
#    output_files[[2]],
#    mt_elbert_overlay,
#    output_file = tempfile(fileext = ".png")
#    )
#  knitr::include_graphics(ortho_with_points)

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

## ----eval = FALSE-------------------------------------------------------------
#  georef_overlay <- georeference_overlay(
#    ortho_with_points,
#    output_files[[2]]
#  )

## ---- eval = FALSE------------------------------------------------------------
#  tile_overlays <- lapply(output_files[[2]],
#                          function(x) vector_to_overlay(mt_elbert_points,
#                                                        x,
#                                                        size = 15,
#                                                        color = "red",
#                                                        na.rm = TRUE))
#  
#  combined_tiles <- mapply(function(x, y) {
#    combine_overlays(x, y, output_file = tempfile(fileext = ".png"))
#    },
#    output_files[[2]],
#    tile_overlays)
#  
#  georef_tiles <- mapply(georeference_overlay, combined_tiles, output_files[[2]])
#  
#  merged_tiles <- merge_rasters(georef_tiles)

## ----eval = FALSE-------------------------------------------------------------
#  elevation_tile <- output_files[[1]]
#  make_manifest(elevation_tile, georef_tiles)

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