I am new to Mrtrix3. I thoroughly studied the BATMAN Tutorial and questions asked by other mrtrix users however I am still a bit confused - maybe you can help …
For my project I need to get the streamline count between a seed and a target region (which are both labels from the DKT atlas of Freesurfer). The streamlines do not have to end in the target region, it is enough if they traverse it. The direction, whether the streamlines go from seed to target or from target to seed is also irrelevant.
After reading through all sorts of posts on this topic, I feel like there are several possibilities for this and am unsure which one is the “correct” one (or if there are several correct ones).
One possibility seems to be to make a whole-brain tractography using tckgen, then tcksift, and then construct a connectome using tck2connectome. Using connectome2tck afterwards, I should be able to get a trackfile between the two nodes of the DKT parcellation atlas. Using tckinfo I should get the streamline count, right?
However as I understand, only the streamlines that end in the ROIs are counted here and not the ones that pass through it. Is this correct? Would I get around this by using the -assignment_all_voxels option?
Another option seems to be to get the streamlines between two regions using tckgen by specifying a seed image and inclusion mask, where you don’t just get the streamlines that end in the inclusion mask, but also the ones that go through them.
2a) I read, that you can apply SIFT only to whole-brain tractography, which means it wouldn´t be applicable in this case?
2b) In the tck2connectome command you can specify the scale_invnodevol option. Is there a way I can consider the size of my ROIs when I do not use tck2connectome that has that option??
Then there seems to be a possibility to use tckedit to edit trackfiles already generated with tckgen - is it possible to extract streamlines between two ROIs from a whole-brain tractography using tckedit and use the -include option twice (so -include ROI1 -include ROI2)?
One last related question: doing functional connectivity analysis we mostly use dilated atlas-labels. Would you also recommend to dilated ROIs in structural connectivity analysis? If yes: I guess you cannot dilate a whole parcellation-image, therefore getting the streamline count between two dilated ROIs by using connectome2tck isn´t possible in that case, right?
I would really appreciate some guidance on how you would approach this analysis.
Thanks in advance,