tests/testthat/test-cog.R

In neuroconductor/neurobase: 'Neuroconductor' Base Package with Helper Functions for 'nifti' Objects

test_that("cog", {
  set.seed(5)
  dims = rep(10, 3)
  
  arr = array(rnorm(prod(dims)), dim = dims)
  nim = oro.nifti::nifti(arr)
  mystr = "FLOAT32"
  datatype(nim) <- convert.datatype()[[mystr]]
  bitpix(nim) <- convert.bitpix()[[mystr]]
  nim = cal_img(nim)
  
  ee = structure(c(5.45383104125737, 5.44597249508841, 5.62671905697446
  ), .Names = c("dim1", "dim2", "dim3"))
  expect_equal(cog(nim), ee)
  
  expect_equal(xyz(nim), ceiling(ee))
  
})