Node coordinates in MNI space?


#1

Hi mrtrix experts!

To go further with connectome analysis and visualization, most of the programs require node coordinates.

What is the easiest way to know them for 84 nodes used in standard “fs_default.txt” I use in mrtrix.

Thank you in advance!


Node Coordinates in MNI space
#2

Hi David,

Definition of the location of the nodes within that parcellation is slightly ambiguous, since the cortical parcels are defined on FreeSurfer’s spherical template rather than within a volumetric template like MNI.

In mrview's connectome tool, this is handled by calculating from the parcellation image the centre of mass of each parcel; the positions are therefore subject-specific.

If it’s only for visualisation and not specifically anatomical reporting, I suppose what you could do is either:

  • Calculate the centre of mass of each parcel in subject space, register your subject’s T1 image to MNI space, then transform each centre of mass according to the estimated transformation;
  • Register your subject’s T1 image to MNI space, transform the subject’s parcellation image to MNI space, then calculate the centre of mass of each parcel.

The calculation of centre-of-mass can be done using this command. To run it for each parcel, you could do something like:

for i in {1..84}; do mrcalc parcellation.mif $i -eq - -quiet | mrcentre - scanner; done

That should give you a 84x3 table at the command-line, containing the XYZ scanner-space position of each parcel.

Rob


#3

Hi Rob,

I too have the same question as David and I’m grateful for the steps you have provided. However I have run into a few issues:
When I run the last step suggested to get MNI coordinates for the parcellated brain I receive this error:
mrcalc: [SYSTEM FATAL CODE: SIGPIPE (13)] Nothing on receiving end of pipe
-bash: mrcentre: command not found

This is after incoporating the ‘mrcentre.cpp’ file in the mrtrix3 script bin. How can I incoporate the mrcentre command into the code for it to work?
Thanks,
Abdul


#4

Hey Abdul,

I ran into the same issue when attempting to run the command specified by Rob. What I did to solve it was move the mrcentre.cpp file into the mrtrix3/cmd directory. Then I ran ./build in the command line in the mrtrix3 directory again to compile the new cpp file. Now you should be able to run Rob’s command from anywhere.

I added an echo -e "\n"; in between the last semicolon and the done so that the output would be easier to copy into a text file and save. Hope this helps!

Rafi


#5
mrcalc: [SYSTEM FATAL CODE: SIGPIPE (13)] Nothing on receiving end of pipe
-bash: mrcentre: command not found

Yes, just placing the file mrcentre.cpp into the cmd/ directory is not sufficient to run the command; it needs to be compiled, just as with all other MRtrix3 commands.

The alternative way to handle cases like this (which can in fact be a lot cleaner) is to set up a separate MRtrix3 “module” / “project”, where additional commands designed to operate against the MRtrix3 libraries can be compiled without interfering with the primary MRtrix3 installation in any way. However I can’t find anywhere in the online documentation where this capability is described (apart from the caveat that is relevant only for Windows). @jdtournier Maybe we need to add this to the “Tips and tricks” section?


#6

Yes, that’s only described in the developer docs, it’s not (yet) mentioned in the user docs.

I’m not sure, it’s quite an unusual thing to need to do – at least at the moment. I guess there’s no harm in adding it in if you have the time?

However, one idea that came up during the ISMRM was that we could add something equivalent to Python’s pip install: just a simple script to download/clone a git repo, set it up as a module, and build its code. That would allow other users to easily contribute their own commands. What do you think?


#7

Hi Rafi,

Very very helpful, thanks!

Abdul