tests/testthat/test_to_nii.R

In LabNeuroCogDevel/LNCDR: Laboratory of Neuro-Cognitive Development R Functions

library(LNCDR)
library(oro.nifti)
context("to_nii")
# test with
# cd ../; Rscript testthat.R

test_that("read write save", {
   d<-data.frame(i=20:30,j=20:30,k=20:30,v=100)
   to_nii('../HarOx-sub-2mm.striatumplusthalamus.nii.gz',d,'test')
   expect_true( file.exists('test.nii.gz') )

   nii <- readNIfTI('test.nii.gz')
   nd <- nii@.Data[nii@.Data>0]

   expect_length(nd, 11 )
   expect_true( all(nd==100) )

   file.remove('test.nii.gz')
  
})

test_that("read write save with charater indexies", {
   d<-data.frame(i=20:30,j=20:30,k=20:30,v=100)
   d$i <- as.character(d$i); d$k <- as.character(d$k); d$j <- as.character(d$j) 
   to_nii('../HarOx-sub-2mm.striatumplusthalamus.nii.gz',d,'test')
   expect_true( file.exists('test.nii.gz') )

   nii <- readNIfTI('test.nii.gz')
   nd <- nii@.Data[nii@.Data>0]

   expect_length(nd, 11 )
   expect_true( all(nd==100) )

   file.remove('test.nii.gz')
  
})

test_that("reset afni min/max", {
   skip_if(system("which 3dBrickStat")!=0)
   d <- data.frame(i=20:30, j=20:30, k=20:30, v=100)
   to_nii("../HarOx-sub-2mm.striatumplusthalamus.nii.gz", d, "test")

   nii <- readNIfTI('test.nii.gz')
   niimax <- max(nii@.Data)
   expect_equal(as.numeric(system("3dBrickStat test.nii.gz", intern=T)), niimax)
   file.remove("test.nii.gz")
})