tests/testthat/test-nii.R

In muschellij2/NiftiArray: HDF5 Delayed Array for Nifti Objects

testthat::context("Trying to make NiftiArray objects")

nii_fname = system.file("extdata", "example.nii.gz",
                        package = "RNifti")
h5_fname = tempfile(fileext = ".h5")
img = RNifti::readNifti(nii_fname)
img_hdr = nifti_header(img)

check_array = function(x) {
  testthat::expect_is(x, "NiftiArray")
  testthat::expect_is(x, "HDF5Array")
  testthat::expect_is(x, "DelayedArray")
}
testthat::test_that("Writing an Array", {
  res = writeNiftiArray(img)
  check_array(res)
})

testthat::test_that("Writing and Reading an Array", {

  res = writeNiftiArray(img, filepath = h5_fname)
  testthat::expect_true(file.exists(h5_fname))
  check_array(res)

  hdr = nifti_header(res)
  testthat::expect_is(hdr, "niftiHeader")
  rm(res)

  res = NiftiArray(h5_fname)
  check_array(res)
})

testthat::test_that("Checking equivalent headers", {

  res = writeNiftiArray(img)
  hdr = nifti_header(res)
  rm(res)

  testthat::expect_true(all(names(img_hdr) == names(hdr)))
  testthat::expect_true(length(img_hdr) == length(hdr))
  xx = mapply(function(x, y) {
    testthat::expect_equal(x, y)
    NULL
  }, img_hdr, hdr)

  a_img_hdr = attributes(img_hdr)
  a_hdr = attributes(hdr)
  print(a_img_hdr)
  a_img_hdr$strings = NULL
  print(a_hdr)
  xx = mapply(function(x, y) {
    testthat::expect_equal(x, y)
    NULL
  }, a_img_hdr, a_hdr)
  rm(xx)

  pd = RNifti::pixdim(hdr)
  testthat::expect_true(all(pd == 2.5))
  testthat::expect_equal(hdr$qoffset_x, 122.033897399902)

})


testthat::test_that("Writing and Reading just nifti file on disk", {

  res = writeNiftiArray(nii_fname)
  check_array(res)

  hdr = nifti_header(res)
  testthat::expect_is(hdr, "niftiHeader")
  rm(res)

  res = NiftiArray(nii_fname)
  check_array(res)
})