library(raster)
library(sf)
library(tidyverse)
source("data-raw/01_borders_download.R")
source("data-raw/02_landsat_download.R")
source("data-raw/03_srtm_download.R")
source("data-raw/04_nlcd_download.R")
## park borders
"inst/vector/nps_boundary.zip" %>%
download_borders(.) %>%
preprocess_borders(.)
## landsat8
filenames <- c("inst/raster/landsat_b2.tif",
"inst/raster/landsat_b3.tif",
"inst/raster/landsat_b4.tif",
"inst/raster/landsat_b5.tif")
bands <- c(2, 3, 4, 5)
map2(filenames, bands, download_landsat8) %>%
map(~preprocess_lansat8(., "inst/vector/zion.gpkg"))
stack(filenames) %>%
writeRaster(., filename = "inst/raster/landsat.tif", overwrite=TRUE, datatype="INT2U", options=c("COMPRESS=DEFLATE"))
## srtm
"inst/raster/srtm_14_05.zip" %>%
download_srtm(.) %>%
preprocess_srtm(., "inst/vector/zion.gpkg")
## nlcd
"inst/raster/nlcd2011" %>%
download_nlcd(.) %>%
preprocess_nlcd(., "inst/vector/zion.gpkg")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.