Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

New script dwicat #1626

merged 2 commits into from Jun 18, 2019


None yet
3 participants
Copy link

commented Jun 8, 2019

Hey Robert,
Could you get it approved by some other developer.

Lena & Jakub

New script dwicat
Co-authored-by: Lena Dorfschmidt <>
Co-authored-by: Robert E. Smith <>

@Lestropie Lestropie added the scripts label Jun 8, 2019

@Lestropie Lestropie self-assigned this Jun 8, 2019

@Lestropie Lestropie marked this pull request as ready for review Jun 8, 2019

dwicat: Give 4D b=0 images to mrhistmatch
Enabled by changes in #1628

Provide all b=0 volumes of each input series to mrhistmatch, rather than explicitly calculating the mean b=0 image. This provides a greater amount of data for the least squares determination of scaling factor between series. This was however not possible prior to #1628, as previously mrhistmatch worked explicitly on 3D images only.

@Lestropie Lestropie requested a review from jdtournier Jun 8, 2019


This comment has been minimized.

Copy link

commented Jun 8, 2019

Closes #1428.

Copy link

left a comment

Looks great, thanks!

My one comment (based on a recent discussion where this functionality will be most welcome) would be to add an option to set the average b=0 intensity - useful in cases where the data might come from different scanners where the scaling might differ wildly between subjects. But it's nothing that can't be done in a subsequent step, so don't let that stop this from merging.

@jdtournier jdtournier merged commit 5cb2397 into MRtrix3:dev Jun 18, 2019

2 checks passed

continuous-integration/appveyor/pr AppVeyor build succeeded
continuous-integration/travis-ci/pr The Travis CI build passed

This comment has been minimized.

Copy link

commented Jun 25, 2019

Manually setting average b=0 intensity would be pretty easy. What I actually wanted to do here, but couldn't convince myself that my math was correct, was to intensity-match all protocols to all other protocols, producing a scaling factor matrix, and then scale each protocol individually in such a way as to converge their intensities toward a common target, rather than arbitrarily using the first scan as the "target" intensity.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
You can’t perform that action at this time.