Right now I'm not interested in the eddy currents. I have DTI data corrected for eddy currents and head motion, and would like to add simulated motion to test how well different algorithms can correct for the motion.
And will this “added” motion then be considered ground truth? That would of course hinge crucially on the initial correction being absolutely perfect, which makes it all a little circular.
I would recommend using the simulations by Graham, Drobnjak and Zhang as your ground truth instead. It is my understanding that they are planning to release simulated diffusion data along with ground truth. If any of them are reading this they might be
able to confirm?
Compared to fMRI data, I think that motion correction of DTI data is more challenging, since DTI volumes have different intensity and noise levels (depending on measurement direction and b-factor). I would therefore like to only use the motion
correction part of eddy or eddy_correct, but maybe that's impossible.
As long as your initial manipulations hasn’t cause areas of NaN (or zeros) at the edges of the FOV it should be fine with either eddy or eddy_correct.
My fMRI software BROCCOLI uses an image registration algorithm that is not based on image intensity, instead it uses local phase from quadrature filters (and the local phase is invariant to the intensity level). My hypothesis is that this algorithm
may be better for images of different contrast/intensity/noise levels, compared to methods based on intensity.
Good luck Jesper