Daniel,
I am having some issues with getting this script to work based on header attributes from 3dWarp. I'm using Version ID = AFNI_2008_07_18_1710. Do I need a newer version of afni for 3dWarp to use the same attributes as 3dWarpDrive?
@rotatevectorsallx1 F0001-0708.DTI.1.warp+orig bvecs_afni.1D
if ( F0001-0708.DTI.1.warp+orig == || F0001-0708.DTI.1.warp+orig == -help || F0001-0708.DTI.1.warp+orig == -h ) then
if ( 2 < 2 ) then
set dset = F0001-0708.DTI.1.warp+orig
set tensor1D = bvecs_afni.1D
1dcat bvecs_afni.1D
set ngrads = `cat rotatedgrads.1D|wc -l`
cat rotatedgrads.1D
wc -l
echo Number of gradients is 20
Number of gradients is 20
@ ngrads--
echo Dataset: F0001-0708.DTI.1.warp+orig
Dataset: F0001-0708.DTI.1.warp+orig
echo Gradient vectors: bvecs_afni.1D
Gradient vectors: bvecs_afni.1D
echo Rotating vectors...
Rotating vectors...
rm -f rotatedgrads.1D
set check = `wc tempmatvec.1D`
wc tempmatvec.1D
if ( 0 == 0 ) then
echo Looking for WARPDRIVE attribute instead
Looking for WARPDRIVE attribute instead
set check = `wc tempmatvec0.1D`
wc tempmatvec0.1D
if ( 0 == 0 ) then
echo No rotation matrix in header
No rotation matrix in header
goto END
