# Replicating longitudinal fixel-based analysis approach

**URL:** <https://community.mrtrix.org/t/replicating-longitudinal-fixel-based-analysis-approach/2071>\
**Category:** Uncategorized\
**Created:** [November 24, 2018, 4:46pm UTC](https://community.mrtrix.org/t/replicating-longitudinal-fixel-based-analysis-approach/2071 "2018-11-24T16:46:59Z")\
**Posts on this page:** 1\
**Showing post:** 7

<div class="post-metadata">

**Author:** ![maxpietsch](https://community.mrtrix.org/user_avatar/community.mrtrix.org/maxpietsch/32/25_2.png) [@maxpietsch](https://community.mrtrix.org/u/maxpietsch)\
**Post date:** [February 15, 2019, 5:23pm UTC](https://community.mrtrix.org/t/replicating-longitudinal-fixel-based-analysis-approach/2071/7 "2019-02-15T17:23:46Z")

</div>

> [@uclpz](#):
>
> population\_template has a -rigid option for the purpose of creating a longitudinal template but not sure how to implement that to make the intra-subject template.

This gives you the average image (intensity or FODs) of the rigidly aligned images. Registration is performed iteratively in a pairwise manner between each subject image and the current best estimate of the template. After each stage, the template is updated with the average of the rigidly transformed images (see diagram [here](http://community.mrtrix.org/t/fa-template-resolution/1856/4)). This is what happens behind the scenes in `population_template`, nothing you need to worry about.

> [@uclpz](#):
>
> I’m told this needs to be registered to the midpoint of the 3 images but not sure how you work up the midpoint and how that works in practice.

The template is the unweighted average intensity of all images (time points) mapped to the common space via subject to template transformations. In your case, the template works as mid-point so you’d transform each image to the subject-specific template space. You can use population\_tempalte with rigid, affine or nonlinear registration to create this central reference space.

Note that there are other registration approaches to align multiple images to a common central point (based on (pairwise) distances in deformation space for instance) but not within MRtrix. For more information on the topic see for instance [The optimal template effect in hippocampus studies of diseased populations, Avants et al.](https://www.sciencedirect.com/science/article/pii/S1053811909010611). If you wanted to adjust the “centre” point of the subject-specific template, you’d need to get your hands dirty as the MRtrix tools were not designed for this. I guess the easiest would be to adjust population\_template to perform a weighted intensity average based on temporal distance from the target time point. This would bias the template’s representation towards the target time point’s representation and therefore nudge registration towards that time point’s shape as well. It is probably not the best™ approach from a theoretical point of view but I’d guess it would work reasonably well in the absence of large deformations between time-steps.

---

_[View the full topic](https://community.mrtrix.org/t/replicating-longitudinal-fixel-based-analysis-approach/2071)._