STRIPSyield: STRIPS yield data

# This is the only script that you should run manully to build the data objects
# to be distributed with the STRIPSYield package. It relies on the shapefiles as
# well as the many functions laid out other .R files stored in the data-raw
# folder. Note that the whole data-raw folder is not part of the user-facing
# part of the STRIPYield package: the raw shapefiles and these scripts are only
# meant to be used locally to build the data objects before the package is
# published. For further questions, see ?whatever.

# 0. Settings -------------------------------------------------------------
yieldPathIn     <- "data-raw/yield_original"    # Folder where the original shapefiles are stored
yieldPathOut    <- "data-raw/yield_curated"     # Folder where the curated shapefiles will be stored
boundaryPathIn  <- "data-raw/boundary_original" # Folder where the original shapefiles are stored
boundaryPathOut <- "data-raw/boundary_curated"  # Folder where the curated shapefiles will be stored
nCores          <- 19                           # Number of cores for point coordinates classification

# 1. Preamble -------------------------------------------------------------
suppressPackageStartupMessages(library(shapefiles))
source('data-raw/R/boundary.R')
source('data-raw/R/build.R')
source('data-raw/R/curate.R')
source('data-raw/R/gis.R')
source('data-raw/R/measurements.R')
source('data-raw/R/shapefiles.R')
source('data-raw/R/yield.R')

# 2. Curate the shapefiles. -----------------------------------------------
unlink(yieldPathOut, recursive = TRUE, force = TRUE)

if (!dir.exists(yieldPathOut))
  dir.create(yieldPathOut)

curate_all_yield_shapefiles(yieldPathIn, yieldPathOut)

# 3. Curate the boundary shapefiles ---------------------------------------
unlink(boundaryPathOut, recursive = TRUE, force = TRUE)

if (!dir.exists(boundaryPathOut))
  dir.create(boundaryPathOut)

curate_all_boundary_shapefiles(boundaryPathIn, boundaryPathOut)

# 3. Build master data.frame ----------------------------------------------
yieldShapes    <- read_all_shapefiles(yieldPathOut)
boundaryShapes <- read_all_shapefiles(boundaryPathOut)
yieldDF        <- build_yield(yieldShapes)
boundaryDF     <- build_boundaries(boundaryShapes)
yieldExtra     <- build_extra(yieldDF, boundaryDF)

warning(
  unique(yieldExtra[yieldExtra$watershed == "OffBounds", ]$year)
)

# 4. Create summary -------------------------------------------------------
yieldSummary   <- "TO DO"

# 5. Rename, subset and export --------------------------------------------
boundaries     <- boundaryDF
yield          <- yieldExtra[
  ,
  c(
    "site", "watershed", "block", "treatment", "year", "crop", "x", "y",
    "moisture", "yieldDryBuAc"
  )
]

usethis::use_data(boundaries, overwrite = TRUE)
usethis::use_data(yield,      overwrite = TRUE)
usethis::use_data(yieldExtra, overwrite = TRUE)

# 6. Build legacy datasets ------------------------------------------------
source('data-raw/R/legacy.R')