R/load_env_data.R

In Doi90/bushfireSOS: Fit Species Distribution Models In Response To Australian Bushfire Crisis

#' Load Environmental Data
#'
#' @param stack_dir Character. The directory to the raster_tile folder
#' @param region
#'
#' @return
#' @export
#'
#' @examples
#'

load_env_data <- function(stack_dir,
                          region){

  ## Add ACT to NSW
  if("NSW" %in% region){
    if(!"ACT" %in% region)
      region <- c(region, "ACT")
  }

  ## make a temp file for the vrt
  tmp <- tempdir(check = FALSE)
  outfile <- file.path(tmp, "bushfire_raster.vrt")

  ## select the regions
  infile <- list.files(stack_dir, pattern = ".tif$")
  inname <- substr(infile, 1, nchar(infile) - 4)
  infile <- infile[which(inname %in% region)]

  ## make the virtual tiles
  gdalUtils::gdalbuildvrt(gdalfile = file.path(stack_dir, infile),
                          output.vrt = outfile,
                          vrtnodata = -9999,
                          overwrite = TRUE)

  ## read the tiles
  tempraster <- raster::brick(outfile)

  return(tempraster)

}