tests/testthat/test-process_otsu_regenera.R

In OtsuFire: Fire Scars, Severity and Regeneration Mapping Using 'Otsu' Thresholding

test_that("process_otsu_regenera runs on synthetic RBR raster with and without fixed threshold", {
  skip_on_cran()

  # Set paths
  python <- get_env_tool("PYTHON_EXE_PATH")
  gdal_polygonize <- get_env_tool("GDAL_POLYGONIZE")
  skip_if(is.null(python) || !file.exists(python), "Python not available")
  skip_if(is.null(gdal_polygonize) || !file.exists(gdal_polygonize), "gdal_polygonize.py not found")

  tmpdir <- file.path(tempdir(), "regenera_test")
  dir.create(tmpdir, showWarnings = FALSE, recursive = TRUE)

  # Create synthetic raster
  r <- terra::rast(nrows = 10, ncols = 10, xmin = 0, xmax = 1000, ymin = 0, ymax = 1000)
  terra::values(r) <- matrix(sample(-200:200, 100, replace = TRUE), nrow = 10)
  terra::crs(r) <- "EPSG:3035"
  rbr_path <- file.path(tmpdir, "RBR_P2.tif")
  terra::writeRaster(r, rbr_path, overwrite = TRUE)

  # --- Test with fixed threshold ---
  result_fixed <- process_otsu_regenera(
    rbr_post = list(P2 = rbr_path),
    output_dir = tmpdir,
    python_exe = python,
    gdal_polygonize_script = gdal_polygonize,
    fire_year = 2000,
    use_fixed_threshold = TRUE,
    fixed_threshold_value = -100
  )

  expect_true("P2" %in% names(result_fixed))
  p2_fixed <- result_fixed$P2[[1]]
  expect_true(all(c("raster", "shapefile") %in% names(p2_fixed)))
  expect_true(file.exists(p2_fixed$raster))
  expect_true(file.exists(p2_fixed$shapefile))
  shp_fixed <- sf::st_read(p2_fixed$shapefile, quiet = TRUE)
  expect_s3_class(shp_fixed, "sf")
  expect_gt(nrow(shp_fixed), 0)

  # --- Test with default threshold (< 0) ---
  result_auto <- process_otsu_regenera(
    rbr_post = list(P2 = rbr_path),
    output_dir = tmpdir,
    python_exe = python,
    gdal_polygonize_script = gdal_polygonize,
    fire_year = 2000,
    use_fixed_threshold = FALSE,
    bind_all=FALSE
  )

  expect_true("P2" %in% names(result_auto))
  p2_auto <- result_auto$P2[[1]]
  expect_true(all(c("raster", "shapefile") %in% names(p2_auto)))
  expect_true(file.exists(p2_auto$raster))
  expect_true(file.exists(p2_auto$shapefile))
  shp_auto <- sf::st_read(p2_auto$shapefile, quiet = TRUE)
  expect_s3_class(shp_auto, "sf")
  expect_gt(nrow(shp_auto), 0)
})