The NaNs should be masked out at the time rast is defined, and not in each classification round.

R-package/R/inner_procedure.r

@@ -71,6 +71,7 @@ sample_nb <- function(raster,
   } else {
     rast <- raster
   }
+  rast <- raster::mask(rast, raster::calc(rast, fun = sum))
   ###
   l <- 1        ###6. opt=260
   pbtn1 <- matrix(1, nrow = 1, ncol = 1)

R-package/R/model_opt.r

@@ -42,7 +42,6 @@ model_opt_r <- function(k,
   points <- NULL
   models <- NULL
   oobe <- matrix(NA, nrow = n, ncol = 1)
-  rast <- raster::mask(raster, raster::calc(raster, fun = sum))
   for (j in 1:n) {
     ###Vorbereitung Klassifizierung
     if (j == 1) {