Bug fix to deal with only NA values along one river

Author: ma-z-am

R/extend_river_table.R

@@ -47,7 +47,7 @@ extend_riverTable <- function(
   data_table <- aggregated_data[aggregated_data$qsimVis_ID == river_id &
                                   !is.na(aggregated_data$qsimVis_ID),]
 
-  if(nrow(data_table) > 0L){
+  if(nrow(data_table) > 0L & any(!is.na(data_table[[varName]]))){
     # apply results to closest verknet node, if not already defined
     km_verknet <- river_table$km
     for(i in seq_len(nrow(data_table))){

R/value_to_classes.R

@@ -12,17 +12,22 @@
 #'
 #' @export
 #'
-value_to_classes <- function(river_list, classBreaks, colorVector = NULL){
-
+value_to_classes <- function(
+    river_list, classBreaks, colorVector = NULL
+){
   output_list <- lapply(names(river_list), function(N){
     r <- river_list[[N]]
     x <- r$data
-    x$value_class <- cut(
-      x = x$value,
-      breaks = classBreaks,
-      include.lowest = TRUE,
-      ordered_result = TRUE
-    )
+    x$value_class <- if(any(!is.na(x$value))){
+      cut(
+        x = x$value,
+        breaks = classBreaks,
+        include.lowest = TRUE,
+        ordered_result = TRUE
+      )
+    } else {
+      NA
+    }
     list("data" = x,
          "pp" = r$pp
     )