Nothing
# devtools::test(".", "calcTOAIrrad")
context("Solar irradiation (ESun) using LUTs")
#-------------------------------------------------------------------------------
test_that("calcTOAIrradModel for 'data.frame' objects works as expected", {
calcTOAIrradModel(lut$L8_RSR, model = "MNewKur")
calcTOAIrradModel(lut$L8_RSR, model = "MNewKur", normalize = FALSE,
esd = calcEarthSunDist("2015-01-01"))
calcTOAIrradModel(lut$L8_RSR, model = "MNewKur", normalize = FALSE,
esd = calcEarthSunDist("2015-07-07"))
calcTOAIrradModel(lut$L7_RSR, model = "MNewKur")
calcTOAIrradModel(lut$L7_RSR, model = "MNewKur", normalize = FALSE,
esd = calcEarthSunDist("2015-01-01"))
calcTOAIrradModel(lut$L7_RSR, model = "MNewKur", normalize = FALSE,
esd = calcEarthSunDist("2015-07-07"))
})
#-------------------------------------------------------------------------------
test_that("calcTOAIrradModel for 'Satellite' objects works as expected", {
## loop over different collections (ie P1L1, C1L1)
path <- c(system.file("extdata", package = "satellite"),
system.file("testdata/LC8", package = "satellite"))
lapply(1:length(path), function(i) {
sid <- ifelse(i == 1, "LC08", "LC8")
files <- list.files(path[i], pattern = glob2rx(paste0(sid, "*.TIF")), full.names = TRUE)
sat <- satellite(files)
test <- calcTOAIrradModel(sat)
expect_equal(as.character(getSatBID(test)[1]), "1")
expect_equal(round(as.numeric(getSatESUN(test)[1]),4), round(1888.4115033, 4))
expect_equal(as.character(getSatBID(test)[2]), "2")
expect_equal(round(as.numeric(getSatESUN(test)[2]),4), round(1974.8429354, 4))
expect_equal(as.character(getSatBID(test)[3]), "3")
expect_equal(round(as.numeric(getSatESUN(test)[3]),4), round(1851.7520559, 4))
expect_equal(as.character(getSatBID(test)[11]), "11")
expect_equal(round(as.numeric(getSatESUN(test)[11]),4), round(0.1068904, 4))
})
})
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.