tests/testthat/test-load_and_reshape.R

In frycast/rsar: What the Package Does (One Line, Title Case)

filename <- system.file(
  "extdata", "MG_VH_sub_norm_testclip.tif", package="rsar")

vh_m1 <- load_SAR_matrix( filename )
vh_b1 <- raster::brick( filename )
vh_m2 <- brick_to_matrix( vh_b1 )
vh_b2 <- matrix_to_brick( vh_m1 )
m1 <- SAR_matrix(m = matrix(0, 4, 3),
           brick_nrow = 4, brick_ncol = 3)


test_that('Tif data loads correctly', {

  x <- load_SAR_matrix(filename)
  testthat::expect_equal(class(x), c('SAR_matrix', 'matrix'))
  testthat::expect_equal(typeof(x), 'double')
  testthat::expect_equal(dim(x), c(2150, 30))
  testthat::expect_equal(attributes(x)$brick_dim, c(43, 50, 30))
  testthat::expect_equal(attributes(x)$crs@projargs,
     "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0")
  testthat::expect_equal(round(attributes(x)$extent@ymax, 3), -37.259)
  testthat::expect_equal(round(attributes(x)$extent@xmax, 3), 140.497)
})


test_that("data loads and reshapes", {

  vh_m1 <- load_SAR_matrix(filename)
  vh_b1 <- raster::brick(filename)
  vh_m2 <- brick_to_matrix(vh_b1)
  vh_b2 <- matrix_to_brick(vh_m1)

  testthat::expect_equal(dim(vh_m1), dim(vh_m2))
  testthat::expect_equal(dim(vh_b1), dim(vh_b2))
  testthat::expect_equal(summary(vh_m1)[1], summary(vh_m2)[1])
  testthat::expect_equal(summary(vh_b1)[1], summary(vh_b2)[1])

  testthat::expect_equal(attr(vh_m1,"brick_dim"), c(43, 50, 30))

  testthat::expect_gt(dim(vh_m2)[1], 1)
  testthat::expect_gt(dim(vh_m2)[2], 1)

  testthat::expect_gt(dim(vh_b2)[1], 1)
  testthat::expect_gt(dim(vh_b2)[2], 1)
  testthat::expect_gt(dim(vh_b2)[3], 1)

  testthat::expect_gt(dim(m1)[1], 1)
  testthat::expect_gt(dim(m1)[2], 1)
})


test_that("Time series data is processed correctly", {

  ts1 <-  choose_order(vh_m1)
  ts2 <-  choose_order(vh_m2)

  testthat::expect_true(ts1 >= 1)
  testthat::expect_true(ts2 >= 1)
})


test_that("NA values are tracked in reshaping", {

  b <- raster::brick( filename )
  b[ 2, 2 ] <- NA

  m1 <- brick_to_matrix( b )
  m2 <- load_SAR_matrix( filename, set_na_indices = 2 )
  expect_equal( nrow( m1 ), 2149 )
  expect_equal( attr( m1, "brick_na_indices" ), 52 )
  expect_equal( nrow( m2 ), 2149 )
  expect_equal( attr( m2, "brick_na_indices" ), 2 )

  b <- matrix_to_brick( m1 )
  expect_equal( dim( b ), c( 43, 50, 30 ) )
})

test_that("Fit function is processed correctly", {

  fit1 <-  fit_AR_to_SAR(vh_m1)
  fit2 <-  fit_AR_to_SAR(vh_m2)

  testthat::expect_equal(dim(fit1)[1], dim(vh_m1)[1])
  testthat::expect_equal(dim(fit2)[1], dim(vh_m2)[1])
})