Try running bedpostx_datacheck on your data folder.
Peace,
Matt.
-----Original Message-----
From: FSL - FMRIB's Software Library [mailto:[log in to unmask]] On Behalf
Of Yi-Shin Sheu
Sent: Thursday, July 09, 2009 3:16 PM
To: [log in to unmask]
Subject: [FSL] DTIFIT gives error message (data from Philips 3T scanner)
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+
|