safely handle all-NA results, skipping RAT and notify user

ncss-tech · May 2, 2024 · 39d3ab1 · 39d3ab1
1 parent 00598ac
commit 39d3ab1
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 2 deletions.
diff --git a/R/ISSR800.R b/R/ISSR800.R
@@ -49,7 +49,7 @@ ISSR800.wcs <- function(aoi, var, res = 800, quiet = FALSE) {
     stop("package 'terra' is required", call. = FALSE)
   }
 
-  # sanity check: aoi specification
+  # sanity check: AOI specification
   if (!inherits(aoi, c('list', 'Spatial', 'sf', 'sfc', 'bbox', 'RasterLayer', 'SpatRaster', 'SpatVector'))) { 
     stop('invalid `aoi` specification', call. = FALSE)
   }
@@ -176,8 +176,17 @@ ISSR800.wcs <- function(aoi, var, res = 800, quiet = FALSE) {
   # set layer name in object
   names(r) <- var.spec$desc
 
+  # test for all NA
+  # if TRUE, we cannot process RATs
+  .allNA <- terra::global(r, fun = "isNA")$isNA == terra::ncell(r)
+
+  # message when all cells are NA
+  if(.allNA) {
+    message('all cells are NA')
+  }
+
   # optional processing of RAT
-  if (!is.null(var.spec$rat)) {
+  if (!is.null(var.spec$rat) && !.allNA) {
 
     # get rat
     rat <- try(suppressWarnings(read.csv(var.spec$rat, stringsAsFactors = FALSE)), silent = TRUE)

diff --git a/R/WCS-utils.R b/R/WCS-utils.R
@@ -659,6 +659,7 @@ WCS_details <- function(wcs = c('mukey', 'ISSR800', 'soilColor')) {
 )
 
 
+## TODO: dynamically get vintage from ... where?
 .mukey.spec <- list(
 
   'gnatsgo' = list(
-Original file line number
+Diff line change
@@ Expand Up @@
     )
+    ## TODO: dynamically get vintage from ... where?
     .mukey.spec <- list(
       'gnatsgo' = list(
@@ Expand Down @@