tests/testthat/test_intensityAnalysis.R

In OpenLand: Quantitative Analysis and Visualization of LUCC

context("intensityAnalysis")

set.seed(159)
demo_raster <- .demo_landscape(year =  2000:2005,
                               res = 1,
                               crs = "+proj=utm +zone=21 +south +ellps=GRS80 +units=m +no_defs")
demo_cont <- contingencyTable(demo_raster, pixelresolution = 1)



test_that("Behavior of the intensityAnalysis", {


  demo_int_km2 <-
    intensityAnalysis(
      dataset = demo_cont,
      category_n = "GUP",
      category_m = "SXQ",
      area_km2 = TRUE
    )

  demo_int_pixel <-
    intensityAnalysis(
      dataset = demo_cont,
      category_n = "GUP",
      category_m = "SXQ",
      area_km2 = FALSE
    )


  expect_visible(intensityAnalysis(
    dataset = demo_cont,
    category_n = "GUP",
    category_m = "SXQ",
    area_km2 = FALSE
  ))

  expect_visible(intensityAnalysis(
    dataset = demo_cont,
    category_n = "GUP",
    category_m = "SXQ",
    area_km2 = TRUE
  ))


  expect_equal(demo_int_km2$interval_lvl$intervalData$U,
               demo_int_pixel$interval_lvl$intervalData$U)

  expect_equal(demo_int_km2$interval_lvl$intervalData$PercentChange[3], 80.26)


  expect_equal(demo_int_km2$category_lvlGain$categoryData$GG_km2,
               demo_int_pixel$category_lvlGain$categoryData$GG_pixel/1000000)
  expect_equal(demo_int_km2$category_lvlLoss$categoryData$GL_km2,
               demo_int_pixel$category_lvlLoss$categoryData$GL_pixel/1000000)

  expect_equal(demo_int_pixel$category_lvlGain$categoryData$GG_pixel[8], 1587)
  expect_equal(demo_int_pixel$category_lvlLoss$categoryData$GL_pixel[10], 1684)

  expect_equal(demo_int_km2$transition_lvlGain_n$transitionData$T_i2n_km2,
               demo_int_pixel$transition_lvlGain_n$transitionData$T_i2n_pixel/1000000)
  expect_equal(demo_int_km2$transition_lvlLoss_m$transitionData$T_m2j_km2,
               demo_int_pixel$transition_lvlLoss_m$transitionData$T_m2j_pixel/1000000)

  expect_equal(demo_int_pixel$transition_lvlGain_n$transitionData$T_i2n_pixel[18], 380)
  expect_equal(demo_int_pixel$transition_lvlLoss_m$transitionData$T_m2j_pixel[15], 382)

  #stationarity

  expect_equal(demo_int_km2$category_lvlGain$categoryStationarity$Stationarity,
               demo_int_pixel$category_lvlGain$categoryStationarity$Stationarity)
  expect_equal(demo_int_km2$category_lvlLoss$categoryStationarity$Stationarity,
               demo_int_pixel$category_lvlLoss$categoryStationarity$Stationarity)

  expect_equal(demo_int_km2$transition_lvlGain_n$transitionStationarity$Stationarity,
               demo_int_pixel$transition_lvlGain_n$transitionStationarity$Stationarity)
  expect_equal(demo_int_km2$transition_lvlLoss_m$transitionStationarity$Stationarity,
               demo_int_pixel$transition_lvlLoss_m$transitionStationarity$Stationarity)

  #class

  expect_s4_class(demo_int_km2$interval_lvl, "Interval")
  expect_s4_class(demo_int_pixel$interval_lvl, "Interval")

  expect_s4_class(demo_int_km2$category_lvlGain, "Category")
  expect_s4_class(demo_int_pixel$category_lvlGain, "Category")
  expect_s4_class(demo_int_km2$category_lvlLoss, "Category")
  expect_s4_class(demo_int_pixel$category_lvlLoss, "Category")

  expect_s4_class(demo_int_km2$transition_lvlGain_n, "Transition")
  expect_s4_class(demo_int_pixel$transition_lvlGain_n, "Transition")
  expect_s4_class(demo_int_km2$transition_lvlLoss_m, "Transition")
  expect_s4_class(demo_int_pixel$transition_lvlLoss_m, "Transition")


})