tests/testthat/test_spct_utils.R

In aphalo/photobiology: Photobiological Calculations

context("thin_wl")

test_that("source_spct", {

  my.spct <- sun.spct

  energy_as_default()
  expect_equal(nrow(thin_wl(my.spct)), 495L)
  expect_known_value(thin_wl(my.spct), "./data/wl-thin-default-value-e")
  expect_known_value(thin_wl(my.spct, max.wl.step = 40), "./data/wl-thin-default-value-step-e")
  expect_known_value(thin_wl(my.spct, max.slope.delta = 0.003), "./data/wl-thin-default-value-slope-e")

  photon_as_default()
  expect_equal(nrow(thin_wl(my.spct)), 496L)
  expect_known_value(thin_wl(my.spct), "./data/wl-thin-default-value-q")
  expect_known_value(thin_wl(my.spct, max.wl.step = 40), "./data/wl-thin-default-value-step-q")
  expect_known_value(thin_wl(my.spct, max.slope.delta = 0.003), "./data/wl-thin-default-value-slope-q")

  unset_radiation_unit_default()
})

test_that("response_spct", {

  my.spct <- ccd.spct

  energy_as_default()
  expect_equal(nrow(thin_wl(my.spct)), 105L)
  expect_known_value(thin_wl(my.spct), "./data/wl-thin-default-value-re")
  expect_known_value(thin_wl(my.spct, max.wl.step = 40), "./data/wl-thin-default-value-step-re")
  expect_known_value(thin_wl(my.spct, max.slope.delta = 0.003), "./data/wl-thin-default-value-slope-re")

  photon_as_default()
  expect_equal(nrow(thin_wl(my.spct)), 119L)
  expect_known_value(thin_wl(my.spct), "./data/wl-thin-default-value-rq")
  expect_known_value(thin_wl(my.spct, max.wl.step = 40), "./data/wl-thin-default-value-step-rq")
  expect_known_value(thin_wl(my.spct, max.slope.delta = 0.003), "./data/wl-thin-default-value-slope-rq")

  unset_radiation_unit_default()
})

test_that("filter_spct", {

  my.spct <- yellow_gel.spct

  Tfr_as_default()
  expect_equal(nrow(thin_wl(my.spct)), 233L)
  expect_known_value(thin_wl(my.spct), "./data/wl-thin-default-value-tfr", update = FALSE)
  expect_known_value(thin_wl(my.spct, max.wl.step = 40), "./data/wl-thin-default-value-step-tfr", update = FALSE)
  expect_known_value(thin_wl(my.spct, max.slope.delta = 0.003), "./data/wl-thin-default-value-slope-tfr", update = FALSE)

  unset_filter_qty_default()
})

test_that("reflector_spct", {

  my.spct <- Ler_leaf_rflt.spct

  expect_equal(nrow(thin_wl(my.spct)), 316L)
  expect_known_value(thin_wl(my.spct), "./data/wl-thin-default-value-rfr")
  expect_known_value(thin_wl(my.spct, max.wl.step = 40), "./data/wl-thin-default-value-step-rfr")
  expect_known_value(thin_wl(my.spct, max.slope.delta = 0.003), "./data/wl-thin-default-value-slope-rfr")
})

test_that("object_spct", {

  my.spct <- Ler_leaf.spct

  expect_error(thin_wl(my.spct))
  expect_equal(nrow(thin_wl(my.spct, col.names = "Tfr")), 353L)
  expect_known_value(thin_wl(my.spct, col.names = "Tfr"), "./data/wl-thin-default-value-otfr")
  expect_known_value(thin_wl(my.spct, col.names = "Rfr"), "./data/wl-thin-default-value-orfr")
})

test_that("chroma_spct", {

  my.spct <- ciexyzCC10.spct

  expect_warning(thin_wl(my.spct))
  expect_equal(suppressWarnings(thin_wl(my.spct)), my.spct)
})

test_that("calibration_spct", {

  my.spct <- calibration_spct(w.length = 400:450, irrad.mult = 1)

  expect_warning(thin_wl(my.spct))
  expect_equal(suppressWarnings(thin_wl(my.spct)), my.spct)
})

test_that("generic_mspct", {

  my.mspct <- source_mspct(list(sun.spct))

  energy_as_default()
  expect_equal(nrow(thin_wl(my.mspct)[[1]]), 495L)
  expect_known_value(thin_wl(my.mspct), "./data/wl-thin-default-value-mspct-e")

  photon_as_default()
  expect_equal(nrow(thin_wl(my.mspct)[[1]]), 496L)
  expect_known_value(thin_wl(my.mspct), "./data/wl-thin-default-value-mspct-q")

  unset_radiation_unit_default()
})