tests/testthat/test-estimateHaze.R

In RStoolbox: Remote Sensing Data Analysis

context("estimateHaze")

mtlFile  <- system.file("external/landsat/LT52240631988227CUB02_MTL.txt", package="RStoolbox")

test_that("all hazeBand specifications work", {
  hb <- list(single = 2, contiguous = c(1:3), noncontiguous = c(1,3), noncontiguous2 = c(2,4))
  for(i in seq_along(hb)){
      expect_is(hdn <- estimateHaze(lsat, hazeBands = hb[[i]], plot = FALSE), "numeric")
      expect_identical(length(hdn), length(hb[[i]]))
      expect_is(hdn2 <- estimateHaze(lsat, hazeBands = names(lsat)[hb[[i]]], plot = FALSE), "numeric")
      expect_identical(hdn, hdn2)
      expect_equal(names(hdn), names(lsat)[hb[[i]]])
      expect_is(hdn <- estimateHaze(lsat, hazeBands = hb[[i]], maxSlope = FALSE, plot = FALSE), "numeric")
  }
  expect_error(estimateHaze(lsat), "specify the band")
})


vals <- unlist(Map(rep, 1:20, c(1:20)^2))
tera <- rast(vals = vals, ncol = 1, nrow = length(vals))
vals[1]<-NA

test_that("correct haze values are found (and deals with NA)", {
    expect_is(hdn <- estimateHaze(tera, hazeBands = 1, darkProp = .02, maxSlope = FALSE, plot = FALSE),  "numeric")
    expect_equal(hdn, c(lyr.1 = 5))
    expect_is(hdn <- estimateHaze(tera, hazeBands = 1, darkProp = .02, maxSlope = TRUE, plot = FALSE),  "numeric")
    expect_equal(hdn, c(lyr.1 = 4))
})