tests/testthat/test-utils.R

In ecmwf/caliver: Calibration and Verification of Gridded Model Outputs

context("utils")

test_that("Testing the file utils.R - .round.choose", {

  # Check whether the results are correct
  expect_equal(.round.choose(41, 5, 1), 45)
  expect_equal(.round.choose(41, 5, 0), 40)

})

test_that("Testing the file utils.R - .background_map_fun", {

  # Check whether the result is correct
  example_box <- as(raster::extent(-180, +180, -90, +90), "SpatialPolygons")
  p <- {
    raster::plot(example_box)
    .background_map_fun()
  }
  expect_equal(length(p$layers), 0)

})

test_that("Testing the file utils.R - .quant_function", {

  # Check whether the result is correct
  x <- structure(c(65.4, 84.4, 78, 89.1, 78.8, 77.9, 52.4, 65.3, 58.3,
                   53, 68.2, 84.9, 82.6, 89.6, 82.4, 84.5, 68.5, 72.8,
                   66.3, 54.3, 69.6, 85.1, 84.9, 89.8, 84.1, 87.8, 76.6,
                   76.6, 70.3, 55, 71, 85.4, 87.2, 90.1, 85.9, 91.1, 84.7,
                   80.3, 74.3, 55.6, 71.8, 85.5, 88.6, 90.2, 87, 93.1, 89.5,
                   82.6, 76.7, 56), .Dim = c(10L, 5L),
                 .Dimnames = list(c("X1993.01.01", "X1993.01.02", "X1993.01.03",
                                    "X1993.01.04", "X1993.01.05", "X1993.01.06",
                                    "X1993.01.07", "X1993.01.08", "X1993.01.09",
                                    "X1993.01.10"),
                                  c("75%", "85%", "90%", "95%", "98%")))
  expect_equal(x, round(.quant_function(b[[1:10]]), 1))

})

test_that("Testing the file utils.R - .convert_long_from_180_to_360", {

  x <- .convert_long_from_180_to_360(-10)
  expect_equal(x, 350)

})

test_that("Testing the file utils.R - .transform_raster", {

  x <- .transform_raster(raster_in = r, variable_name = "FWI")
  expect_equal(raster::cellStats(x, "sum"), raster::cellStats(r, "sum"))

  x <- .transform_raster(raster_in = r, variable_name = "BasisRegions")
  expect_equal(raster::cellStats(x, "sum"), raster::cellStats(r, "sum"))

})

test_that("Testing the file utils.R - .create_rat", {

  x2compare <- structure(list(ID = structure(1:3, .Label = c("1", "2", "3"),
                                             class = "factor"),
                              Class = structure(1:3, .Label = c("A", "B", "C"),
                                                class = "factor")),
                         row.names = c(NA, 3L), class = "data.frame")

  x <- .create_rat(ids = c(1, 2, 3), classes = c("A", "B", "C"))
  expect_equal(x, x2compare)

})

test_that("Testing the file utils.R - .get_layers_for_clima", {

  x <- .get_layers_for_clima(b = b, raster_date = as.Date("2016-01-01"),
                             expand = FALSE)
  expect_equal(nlayers(x), 3)

})

test_that("Testing the file utils.R - .classify_hazard", {


  # Check whether the result is correct
  x <- .classify_hazard(b[[1]][], b[[2]][], b[[3]][],
                        b[[4]][], b[[5]][], b[[6]][])

  # Check the function does not generate a dummy layer
  expect_equal(sum(x), 6)

})

test_that("Testing the file utils.R - .utci_classification", {

  # Check whether the result is correct
  x <- .utci_classification(rtp = r)

  # Check the function does not generate a dummy layer
  expect_equal(raster::cellStats(x, sum)[[1]], 120)

})