In USDAForestService/FIESTA: Forest Inventory Estimation and Analysis

knitr::opts_chunk$set(echo = TRUE)

library(FIESTA)

WYbhdistfn <- system.file("extdata",
                          "sp_data/WYbighorn_districtbnd.shp",
                          package = "FIESTA")

rastcont <- c("conus_dem_90m.tif",
              "usfs_2016_CONUS_tcc_analytical_90m.tif",
              "conus_tpi_90m.tif")

rastcont_name <- c("dem", "tcc", "tpi")

rastcat <- "evt14_tnt_90m.tif"

rastcat_name <- "tnt"

library(tidyverse)
stunitco_WY <- FIESTA::stunitco %>%
  filter(STATENM == "Wyoming")

get_aux <- spGetAuxiliary(xyplt = FIESTA::WYpltassgn,
               uniqueid = "CN",
               unit_layer = stunitco_WY,
               unitvar = "COUNTYCD",
               rastlst.cont = rastcont,
               rastlst.cont.name = rastcont_name,
               rastlst.cat = rastcat,
               rastlst.cat.name = rastcat_name,
               spMakeSpatial_opts = spMakeSpatial_options(
                 xvar = "LON_PUBLIC",
                 yvar = "LAT_PUBLIC"
               ))

Run this for creating data object Rda's:

get_aux <- readRDS("get_aux.rds")

WYpltassgn <- get_aux$pltassgn
WYunitzonal <- get_aux$unitzonal %>%
  rename(ESTN_UNIT = COUNTYCD) # so that join will be easier and cleaner for examples

# save as RDA to data folder
save(WYunitzonal, file = "../data/WYunitzonal.rda")
save(WYpltassgn, file = "../data/WYpltassgn.rda")