Added option to download 10 m DEM in addition to 30 m DEM. Added missing processing of subbasins. Bug fixes.

Author: ryanrbart

DESCRIPTION

@@ -11,7 +11,7 @@ Description: RHESSys Preprocessing is a set of functions that will produce the f
 License: MIT + file LICENSE
 Encoding: UTF-8
 LazyData: true
-RoxygenNote: 7.2.3
+RoxygenNote: 7.3.0
 Imports:
   terra,
   sf,

R/clip_landscape_rasters_by_basin.R

@@ -31,7 +31,7 @@ clip_landscape_rasters_by_basin <- function(work_folder,
 
     raster_to_clip <- terra::rast(file.path(work_folder, rasters_to_clip[x]))
     clipped_raster <- terra::mask(raster_to_clip, mask = basin)
-    terra::writeRaster(clipped_raster, filename = file.path(work_folder, names_output[x]), overwrite=TRUE)
+    terra::writeRaster(clipped_raster, filename = file.path(work_folder, paste0(names_output[x], ".tif")), overwrite=TRUE)
 
     return(clipped_raster)
   })

R/download_terrain_variables.R

@@ -12,12 +12,15 @@
 #' @param variable Name of the desired variable. Current options include "dem",
 #'   and "nlcd".
 #' @param label The name of the study area.
+#' @param res_ned The resolution of the NED product. '1' indicates the 1 arc-second
+#'   (~30 m) NED (the default). A '13' indicates the 1/3 arc-second (~10 m)
+#'   dataset.
 #' @param nlcd_year Year of nlcd data. Acceptable values are 2019 (default),
 #'   2016, 2011, 2008, 2006, 2004, and 2001.
-#' @param res Desired resolution in meters.
-#' @param res_example An spatraster with desired resolution for the downloaded
+#' @param res_final Desired final resolution in meters.
+#' @param res_final_example An spatraster with desired resolution for the downloaded
 #'   data to be resampled to. In many cases, this will be same raster as bbox.
-#'   If both `res` and `res_example` are not null, `res_example` takes
+#'   If both `res_final` and `res_final_example` are not null, `res_final_example` takes
 #'   precedent.
 #' @param output_file Path and file name of raster to be saved.
 #'
@@ -30,9 +33,10 @@ download_terrain_variables <- function(bbox,
                                        data_source,
                                        variable,
                                        label,
+                                       res_ned = "1",
                                        nlcd_year = NULL,
-                                       res = NULL,
-                                       res_example = NULL,
+                                       res_final = NULL,
+                                       res_final_example = NULL,
                                        output_file){
 
   # -----------------
@@ -47,6 +51,7 @@ download_terrain_variables <- function(bbox,
     if (variable == "dem"){
       terrain_raster <- FedData::get_ned(template = bbox,
                                          label = label,
+                                         res = res_ned,
                                          force.redo = TRUE)
       method <- "bilinear"
     }
@@ -62,24 +67,24 @@ download_terrain_variables <- function(bbox,
   # -----------------
   # Reproject
 
-  terrain_raster_reproject <- terra::project(terrain_raster, y = terra::crs(paste0("epsg:", proj_epsg)))
+  terrain_raster_reproject <- terra::project(terrain_raster, y = terra::crs(paste0("epsg:", proj_epsg)), method = method)
 
   # -----------------
   # Change resolution
 
   # Resample function in Terra requires an example of the raster with the
-  # correct new resolution. If 'res_example' is available, we use it. Otherwise,
-  # we create an empty 'res_example' raster with the new resolution. Note that
+  # correct new resolution. If 'res_final_example' is available, we use it. Otherwise,
+  # we create an empty 'res_final_example' raster with the new resolution. Note that
   # the terra::res function keeps the xmin and ymin extents the same, but
   # slightly adjusts the xmax and ymax extents to permit the new resolution. We
-  # then resample from the original raster to the 'res_example' raster.
+  # then resample from the original raster to the 'res_final_example' raster.
 
   # Create raster with new resolution
-  if (!is.null(res_example)){
-    terrain_raster_empty <- res_example
+  if (!is.null(res_final_example)){
+    terrain_raster_empty <- res_final_example
   } else {
     terrain_raster_empty <- terrain_raster_reproject
-    terra::res(terrain_raster_empty) <- res
+    terra::res(terrain_raster_empty) <- res_final
   }
 
   # Resample
@@ -93,6 +98,7 @@ download_terrain_variables <- function(bbox,
   # -----------------
   # Write output
 
+  setGDALconfig("GDAL_PAM_ENABLED", "FALSE")     # Prevent xml files from being generated
   terra::writeRaster(terrain_raster_res, filename = output_file, overwrite=TRUE)
 
   return(terrain_raster_res)

R/make_landscape_rasters.R

@@ -101,6 +101,11 @@ make_landscape_rasters <- function(dem,
                                 output = file.path(work_folder, "streams_landscape.tif"),
                                 threshold = stream_threshold)
 
+  # Subbasins
+  whitebox::wbt_subbasins(d8_pntr = file.path(work_folder, "d8_pointer_landscape.tif"),
+                          streams = file.path(work_folder, "streams_landscape.tif"),
+                          output = file.path(work_folder, "subbasin_landscape.tif"))
+
   # -----------------
   # Create patch maps