1#!/bin/sh 2 3if [ $# -lt 2 ] 4then 5echo Missing nifti tool and Binary directory name 6exit 1 7fi 8 9NT=$1 10DATA=$2 11OUT_DATA=$(dirname ${DATA}) #Need to write to separate directory 12cd ${OUT_DATA} 13 14${NT} -quiet -disp_ts 19 36 11 -infiles ${DATA}/stat0.nii | tee o.08.ts.19.36.11 15 16if [ $? -ne 0 ] ; then 17echo disp_ts failed 18exit 1 19fi 20rm -f o.08.ts.19.36.11 21 22${NT} -quiet -disp_ts 19 36 11 -infiles ${DATA}/stat0.nii \ 23 | awk '{print $(NF-1), $NF, $1, $2}' | tee o.09.ts4.1.awk 24 25if [ $? -ne 0 ] ; then 26echo disp_ts failed 27exit 1 28fi 29 30${NT} -quiet -disp_ts 19 36 11 -infiles ${OUT_DATA}/f4.nii | tee o.09.ts4.2.awk 31 32if diff o.09.ts4.1.awk o.09.ts4.2.awk 33then 34echo "" 35else 36echo '** failure, ts4 files differ' 37exit 1 38fi 39rm -f 0.09.ts4.1.awk o.09.ts4.2.awk ${OUT_DATA}/f4.nii 40exit 0 41 42# get the time series for a slice, get the same sub-bricks, 43# compare nifti_images, and display the time series again 44# 45# this should match the previous 46 47if ${NT} -keep_hist -cci 19 36 -1 -1 0 0 0 -prefix ${OUT_DATA}f.19.36 \ 48-infiles ${DATA}/stat0.nii 49then 50echo "" 51else 52echo cci failed 53exit 1 54fi 55 56if ${NT} -keep_hist -cbl -prefix ${OUT_DATA}/f.19.36.t4.nii \ 57 -infiles ${OUT_DATA}/f.19.36.nii'[178..$,0,1]' 58then 59echo "" 60else 61echo cbi failed 62exit 0 63fi 64 65${NT} -diff_nim -infiles ${OUT_DATA}/f.19.36.nii ${OUT_DATA}/f.19.36.t4.nii | tee o.10.diff_nim 66if [ $? -ne 0 ] ; then 67echo f.19.36.nii and f.19.36.t4.nii differ 68else 69echo f.19.36.nii and f.19.36.t4.nii do not differ 70fi 71 72${NT} -quiet -disp_ci 0 0 11 -1 0 0 0 -infiles ${OUT_DATA}/f.19.36.t4.nii \ 73 | tee o.10.dci.4 74#diff o.09.ts4.1.awk o.10.dci.4 75#if( $status ) echo '** failure, o.09 and o.10 ts files differ' 76 77exit 0 78