Lausanne Parcellation: Difference between revisions

From CCN Wiki
Jump to navigation Jump to search
(Created page with "The Lausanne Parcellation is a workflow that subdivides the 30+ general brain regions into many smaller partitions of approximately equal size. This gives us a finer-grained p...")
 
No edit summary
Line 2: Line 2:


The Lausanne 2008 parcellation scheme refers to the 2008 PLoS Biology Paper [http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.0060159 Mapping the Structural Core of Human Cerebral Cortex]
The Lausanne 2008 parcellation scheme refers to the 2008 PLoS Biology Paper [http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.0060159 Mapping the Structural Core of Human Cerebral Cortex]
== Lausanne Parcellation Workflow ==
=== Recon-all ===
A precondition for subparcellating the cortex in surface space is that you will first require a parcellated cortical surface mesh on which to operate. You will need to take a T1 (i.e., anatomical; MPRAGE) volume to the end of the recon-all pipeline, by running the autorecon1, autorecon2 and autorecon3 steps, as outlined elsewhere.
=== ConnectomeMapper in Python ===
The Lausanne 2008 subparcellation is most easily carried out using the ConnectomeMapper Python library. This should already be installed on your workstation. I believe these libraries use calls to FreeSurfer to carry out the parcellation, but I am not totally sure that this is the case, nor what actually goes on under the hood. Nonetheless, the connectome mapper toolkit (cmtk) Python interface provides an easy way to get this done using the Python scripting language.
First, in your Linux terminal, navigate to your $SUBJECTS_DIR
cd $SUBJECTS_DIR
Hopefully this is indeed the directory that contains the folder for the subject whose surface mesh you wish to sub-parcellate. If not, you will need to navigate to the correct directory, and reset your SUBJECTS_DIR environment variable as required.
Next, launch ipython in your terminal window:
ipython
You will now find yourself in a python interactive command shell that looks something like this:
Python 2.7.6 (default, Jun 22 2015, 17:58:13)
Type "copyright", "credits" or "license" for more information.
IPython 1.2.1 -- An enhanced Interactive Python.
?        -> Introduction and overview of IPython's features.
%quickref -> Quick reference.
help      -> Python's own help system.
object?  -> Details about 'object', use 'object??' for extra details.
In [1]:
Suppose you wish to sub-parcellate the subject with ID 605_001. You would type the following commands in the interactive Python shell:
import os
import nipype.interfaces.cmtk as cmtk
parcellate = cmtk.Parcellate()
parcellate.inputs.subject_id="605_001"
parcellate.inputs.subjects_dir=os.environ['SUBJECTS_DIR'];
parcellate.run();
That's it. If you had a different participant in mind, you would simply provide a different folder name for the <code>parcellate.inputs.subject_id</code> value.
The parcellation step will run for quite a while.

Revision as of 10:58, 29 April 2016

The Lausanne Parcellation is a workflow that subdivides the 30+ general brain regions into many smaller partitions of approximately equal size. This gives us a finer-grained partitioning of the cortex. Because they are smaller, these partitions are more likely to behave relatively more homogeneously during a particular task. Becaues they are of comparable size, it is more justifiable to weight them equally in your analyses.

The Lausanne 2008 parcellation scheme refers to the 2008 PLoS Biology Paper Mapping the Structural Core of Human Cerebral Cortex

Lausanne Parcellation Workflow

Recon-all

A precondition for subparcellating the cortex in surface space is that you will first require a parcellated cortical surface mesh on which to operate. You will need to take a T1 (i.e., anatomical; MPRAGE) volume to the end of the recon-all pipeline, by running the autorecon1, autorecon2 and autorecon3 steps, as outlined elsewhere.

ConnectomeMapper in Python

The Lausanne 2008 subparcellation is most easily carried out using the ConnectomeMapper Python library. This should already be installed on your workstation. I believe these libraries use calls to FreeSurfer to carry out the parcellation, but I am not totally sure that this is the case, nor what actually goes on under the hood. Nonetheless, the connectome mapper toolkit (cmtk) Python interface provides an easy way to get this done using the Python scripting language.

First, in your Linux terminal, navigate to your $SUBJECTS_DIR

cd $SUBJECTS_DIR

Hopefully this is indeed the directory that contains the folder for the subject whose surface mesh you wish to sub-parcellate. If not, you will need to navigate to the correct directory, and reset your SUBJECTS_DIR environment variable as required.

Next, launch ipython in your terminal window:

ipython

You will now find yourself in a python interactive command shell that looks something like this:

Python 2.7.6 (default, Jun 22 2015, 17:58:13) 
Type "copyright", "credits" or "license" for more information. 

IPython 1.2.1 -- An enhanced Interactive Python.
?         -> Introduction and overview of IPython's features.
%quickref -> Quick reference.
help      -> Python's own help system.
object?   -> Details about 'object', use 'object??' for extra details.

In [1]: 

Suppose you wish to sub-parcellate the subject with ID 605_001. You would type the following commands in the interactive Python shell:

import os
import nipype.interfaces.cmtk as cmtk
parcellate = cmtk.Parcellate()
parcellate.inputs.subject_id="605_001"
parcellate.inputs.subjects_dir=os.environ['SUBJECTS_DIR'];
parcellate.run();

That's it. If you had a different participant in mind, you would simply provide a different folder name for the parcellate.inputs.subject_id value. The parcellation step will run for quite a while.