# Nonlinear registration diffusion to structural

**URL:** <https://community.mrtrix.org/t/nonlinear-registration-diffusion-to-structural/1136>\
**Category:** Uncategorized\
**Created:** [August 29, 2017, 3:06pm UTC](https://community.mrtrix.org/t/nonlinear-registration-diffusion-to-structural/1136 "2017-08-29T15:06:52Z")\
**Posts on this page:** 4\
**Page:** 1

<div class="post-metadata">

**Author:** ![olie](https://community.mrtrix.org/letter_avatar_proxy/v4/letter/o/c4cdca/32.png) [@olie](https://community.mrtrix.org/u/olie)\
**Post date:** [August 29, 2017, 3:06pm UTC](https://community.mrtrix.org/t/nonlinear-registration-diffusion-to-structural/1136/1 "2017-08-29T15:06:53Z")

</div>

Dear experts,

I wanted to know if any of the mrtrix people routinely (confidently) use in their pipeline any nonlinear registration methods between diffusion and structural data, and if so how the particular registration method is integrated with the mrtrix3 code. I am interested in intrasubject registration following denoising/ eddy (topup) corrections/recombinations/ masking/ skull stripping/ bias field correction. For example, DTI-TK has been suggested/published to perform well, but often the preprocessing steps (for tensor calculation in this example) make excursions to camino and other software packages, making things a bit complicated.  
Thank you,

Octavian

---

<div class="post-metadata">

**Author:** ![tbilliet](https://community.mrtrix.org/user_avatar/community.mrtrix.org/tbilliet/32/3859_2.png) [@tbilliet](https://community.mrtrix.org/u/tbilliet)\
**Post date:** [August 30, 2017, 9:00am UTC](https://community.mrtrix.org/t/nonlinear-registration-diffusion-to-structural/1136/2 "2017-08-30T09:00:07Z")

</div>

Hi Octavian,

Interesting topic, especially in the case where no reversed phase encoded data is available (else, an affine registration should be sufficient).  
Maybe this thread can already give you some information of what you could try: [Distortion correction using T1](https://community.mrtrix.org/t/distortion-correction-using-t1/975)

In the past I have experimented by using ANTS with registering the FA image to a WM mask, combined with the registration between the average DWI and the T1 brain image (see Fig. 1 in this paper: [https://www.ncbi.nlm.nih.gov/pubmed/25840837](https://www.ncbi.nlm.nih.gov/pubmed/25840837)). Although it worked well and results seemed very nice, I was told that using the WM mask (which has very strong intensity gradients) can potentially introduce strong local deformations, so perhaps there exist better ways (see the thread I linked above). After that I haven’t experimented with it again.

N.B. I wasn’t yet using mrtrix back then, but ANTS registrations are easily importable into mrtrix (see the `transformconvert` command). So any other intra-subject registration approach based on ANTS can be integrated in your mrtrix pipeline.

Kind regards,

Thibo

---

<div class="post-metadata">

**Author:** ![tbilliet](https://community.mrtrix.org/user_avatar/community.mrtrix.org/tbilliet/32/3859_2.png) [@tbilliet](https://community.mrtrix.org/u/tbilliet)\
**Post date:** [August 30, 2017, 9:03am UTC](https://community.mrtrix.org/t/nonlinear-registration-diffusion-to-structural/1136/3 "2017-08-30T09:03:26Z")

</div>

Oh by the way, I replied although I don’t consider myself one of the experts you refer to 😛

---

<div class="post-metadata">

**Author:** ![rsmith](https://community.mrtrix.org/user_avatar/community.mrtrix.org/rsmith/32/2672_2.png) [@rsmith](https://community.mrtrix.org/u/rsmith)\
**Post date:** [September 5, 2017, 7:45am UTC](https://community.mrtrix.org/t/nonlinear-registration-diffusion-to-structural/1136/4 "2017-09-05T07:45:02Z")

</div>

> I am interested in intrasubject registration following … eddy (topup) corrections …

If you’ve performed image pre-processing including `topup`, then there should be _no need_ to perform _non-linear_ intra-subject registration: the dominant source of non-linear geometric differences between the two modalities is precisely what `eddy` & `topup` correct for. Indeed even affine registration has too many degrees of freedom; it should require a rigid-body registration only.

Trying to perform such correction in the absence of reversed phase-encode data is another topic entirely. I’ve heard many stories of people doing such things, but mostly in ways of which I would be highly sceptical. While it’s not implemented in _MRtrix3_, [this paper](https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4461504/) would be worth reading for the sake of fully understanding the difficulty of the problem.

Cheers  
Rob