tests/testthat/test-sam.R

In RStoolbox: Tools for Remote Sensing Data Analysis

context("spectral angle mapper")

data(lsat)
lsat[1] <- NA
lsat[2,1] <- NA
pts <- data.frame(x = c(624720, 627480), y = c(-414690, -411090))
endmembers <- extract(lsat, pts)
rownames(endmembers) <- c("water", "vegetation")

test_that("returns correct classes and deals with NA",{
            expect_is(ls <- sam(lsat, endmembers, angles = TRUE), "RasterBrick")
            expect_equal(names(ls), c("water_sa", "vegetation_sa"))
            expect_true(all(is.na(ls[1])))
            expect_true(all(is.na(ls[1,1])))
            expect_is(ls <- sam(lsat, endmembers, angles = FALSE), "RasterLayer")
            expect_equal(names(ls), c("class"))
            expect_true(is.na(ls[1]))
            expect_true(is.na(ls[2,1]))
            
        }
)

sem_mat <- endmembers[1,,drop=FALSE]
sem_vec <- endmembers[1,]
sem_df  <- as.data.frame(endmembers)

test_that("endmember class", {
            expect_is(ls <- sam(lsat, sem_mat, angles = TRUE), "RasterLayer")
            expect_is(ls <- sam(lsat, sem_vec, angles = TRUE), "RasterLayer")
            expect_is(ls <- sam(lsat, sem_df, angles = TRUE), "RasterBrick")
            expect_is(ls <- sam(lsat, sem_df, angles = FALSE), "RasterLayer")
            expect_error(ls <- sam(lsat, sem_vec, angles = FALSE), "only one class")
        })