data-raw/00_data_prep.R

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")
Nowosad/spDataLarge documentation built on Oct. 13, 2023, 11:45 p.m.