tests/testthat/test-RastStackData.R

In ICvectorfields: Vector Fields from Spatial Time Series of Population Abundance

test_that("RastStackData produces a raster stack with appropriate dimensions", {
  xyzdf <- expand.grid(x = c(1:3), y = c(1:3))
  xyzdf$z1 <- runif(9)
  xyzdf$z2 <- runif(9)
  xyzdf$z3 <- runif(9)

  zstack <- RastStackData(xyzdf)
  zstack <- terra::rast(zstack)
  expect_equal(dim(zstack), c(3, 3, 3))
})