Generating fibre-to-field angle from DTI

[lukevano]

Dear all,

I think I may have found a simple solution based on the advice given in this thread

Fibre orientation wrt B0

Short answer: yes it’s perfectly possible. If you want to use DTI to do this, something like this should do the trick: dwi2tensor dwi.mif dt.mif tensor2metric dt.mif -modulate none -vec vec.mif mrconvert vec.mif -coord 3 2 z-component.mif mrcalc z-component.mif -acos angle.mif Or as a one-liner using Unix pipes: dwi2tensor dwi.mif - | tensor2metric - -modulate none -vec - | mrconvert - -coord 3 2 - | mrcalc - -acos angle.mif If you’re using fixel data, you’ll need to take a close look throug…

Would the following code work to generate an image where each voxel value is the fibre orientation relative to B0?

mrconvert data.nii.gz dwi.mif -fslgrad bvecs bvals
dwi2response tournier dwi.mif response.txt
dwi2fod csd dwi.mif response.txt fod.mif
sh2peaks -num 1 fod.mif first_peak.mif
mrconvert first_peak.mif -coord 3 2 z-component.mif
mrcalc z-component.mif -acos angle.mif
mrconvert angle.mif angle.nii.gz