Closed tashrifbillah closed 4 years ago
Issue illustration by @yrathi
Only one word has been changed in the commit:
Please use temporary branch: https://github.com/pnlbwh/multi-shell-dMRIharmonization/tree/fix-warp
Again, it's already set up in our platforms.
It has been observed that lower bvalue-shell diffusion and rish measures are being directly warped to MNI space instead of being warped to template space first and then to MNI space.
Reason:
In the first pass, for bmax-shell, only
prefix*FA*1Warp.nii.gz
transforms are present that properly warps subject space to template space.Later during debugging, template space is being warped to MNI space using
prefix*FA*ToMNI1Warp.nii.gz
However, in the second pass, with pattern
prefix*FA*1Warp.nii.gz
, there are two transforms present:prefix*FA*1Warp.nii.gz
andprefix*FA*ToMNI1Warp.nii.gz
. Thus, for lower bvalue-shells, scale maps are indeed created from MNI space data, not the template space data. This error propagated down the pipeline and resulted in bad result, again for lower bvalue-shells only.Refer to codes: https://github.com/pnlbwh/multi-shell-dMRIharmonization/blob/c7ce6a4412bd7e91332949fcf170788d13796809/lib/buildTemplate.py#L50
and
https://github.com/pnlbwh/multi-shell-dMRIharmonization/blob/c7ce6a4412bd7e91332949fcf170788d13796809/lib/debug_fa.py#L64