tests/testthat/test-35-coerce.R

In tractor.base: Read, Manipulate and Visualise Magnetic Resonance Images

context("Compatibility")

reportr::setOutputLevel(Warning)
options(reportrStderrLevel=reportr::OL$Fatal)

test_that("we can interoperate with RNifti", {
    skip_if_not_installed("RNifti")
    
    path <- system.file("extdata", "nifti", "maskedb0.nii.gz", package="tractor.base")
    
    image <- readImageFile(path)
    expect_equal(RNifti::pixdim(image), rep(2.5,3L))
    expect_equal(RNifti::pixunits(image), c(spatial="mm",temporal="s"))
    
    image <- RNifti::retrieveNifti(image)
    expect_equal(RNifti::pixdim(image), rep(2.5,3L))
    expect_equal(RNifti::pixunits(image), c("mm","s"))
    
    image <- as(RNifti::readNifti(path,internal=FALSE), "MriImage")
    expect_equal(image$getVoxelDimensions(), rep(2.5,3L))
    expect_equal(image$getVoxelUnits(), c(spatial="mm",temporal="s"))
    
    image <- as(RNifti::readNifti(path,internal=TRUE), "MriImage")
    expect_equal(image$getVoxelDimensions(), rep(2.5,3L))
    expect_equal(image$getVoxelUnits(), c(spatial="mm",temporal="s"))
})

test_that("we can interoperate with oro.nifti", {
    skip_if_not_installed("oro.nifti")
    
    path <- system.file("extdata", "nifti", "maskedb0.nii.gz", package="tractor.base")
    
    image <- readImageFile(path)
    expect_true(validObject(as(image, "nifti")))
    
    image <- suppressWarnings(as(oro.nifti::readNIfTI(path), "MriImage"))
    expect_equal(image$getVoxelDimensions(), rep(2.5,3L))
    expect_equal(image$getVoxelUnits(), c(spatial="mm",temporal="s"))
})