Print

Print


Hi FSL Experts,

I used "dcm2nii" to convert a set of PAR/REC (v3) DTI data acquired from 3T Philips 
Gyroscan at the F.M. Kirby Research Center to NIFTI format.  dcm2nii did not 
automatically output the bvecs and bvals data, so I use a free matlab script 
(dti_gradient_table_creator.m) to create the gradient direction data.

Suggested by the previous threads, I also eliminate the last entry of the gradient table 
[100,100,100], which left me with 33 entry (b0+32 directions).  I also used fslroi to 
eliminate the last volume which said to be a diffusion traced image.

However, when I run dtifit, I got the following error message: 
ysheu$ dtifit -k data.nii.gz -o dti -m nodif_brain_mask.nii.gz -r bvecs -b bvals

"terminate called after throwing an instance of 'NEWMAT::IndexException'
Abort trap"

I really don't know what went wrong.  The number of volumes matched to the bvecs and 
bvals, so it should not be the problem.  Is there something wrong with my DTI data?

Any input is appreciated.

Thanks.

Yi-Shin

PS.  I attached my bvecs and bvals and below is the info related to my data.nii.gz 4D DTI 
image:
ysheu$ fslinfo data.nii.gz 
data_type      FLOAT32
dim1           256
dim2           256
dim3           60
dim4           33
datatype       16
pixdim1        0.8281250000
pixdim2        0.8281250000
pixdim3        2.2000000477
pixdim4        6.0850000381
cal_max        0.0000
cal_min        0.0000
file_type      NIFTI-1+