Nothing
context("wq_calc")
test_that("wq_calc works correctly", {
s2 = terra::rast(system.file("raster/S2_Harsha.tif", package = "waterquality"))
# test number of output layers
expect_equal(terra::nlyr(wq_calc(s2, alg = c("Am092Bsub"), sat = "sentinel2")), 1)
expect_equal(terra::nlyr(wq_calc(s2, alg = c("TurbChip09NIROverGreen", "Am092Bsub"), sat = "sentinel2")), 2)
# warnings
expect_warning(wq_calc(s2, alg = c("Al10SABI", "Am092Bsub", "my_alg"), sat = "sentinel2"))
expect_warning(wq_calc(s2, alg = c("my_alg1", "Al10SABI", "Am092Bsub", "my_alg"), sat = "sentinel2"))
# errors
expect_error(wq_calc(s2[[1]], alg = "all", sat = "sentinel2"))
expect_error(wq_calc(s2, alg = "all", sat = "my_satellite"))
expect_error(wq_calc(s2[[1:4]], alg = "all", sat = "sentinel2"))
expect_error(wq_calc(s2, alg = "many", sat = "sentinel2"))
# wrong alg error
l8 = terra::rast(system.file("raster/L8_Taylorsville.tif", package = "waterquality"))
expect_error(wq_calc(raster_stack = l8, alg = "MM12NDCI", sat = "landsat8"))
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.