Freesurfer Subparcellation

From CCN Wiki
Revision as of 14:49, 21 December 2016 by Chris (talk | contribs) (Created page with "Freesurfer annotations and labels describe swaths of cortical regions, defined anatomically during the recon-all process, or functionally or statistically (e.g., significant G...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Freesurfer annotations and labels describe swaths of cortical regions, defined anatomically during the recon-all process, or functionally or statistically (e.g., significant GLM clusters). Tools exist to further subdivide these regions to the desired level of granularity.

Using mris_divide_parcellation

This program divides one or more parcellations into divisions perpendicular to the long axis of the label. The number of divisions can be specified in one of two ways, depending upon the nature of the fourth argument.

First, a splitfile can be specified as the fourth argument. The splitfile is a text file with two columns. The first column is the name of the label in the source annotation, and the second column is the number of units that label should be divided into. The names of the labels depends upon the source parcellation. For aparc.annot and aparc.a2005.annot, the names can be found in $FREESURFER_HOME/FreeSurferColorLUT.txt. For aparc.annot, the labels are between the ranges of 1000-1034. For aparc.a2005s.annot, the labels are between the ranges of 1100-1181. The name for the label is the name of the segmentation without the 'ctx-lh'. Note that the name included in the splitfile does not indicate the hemisphere. For example, 1023 is 'ctx-lh-posteriorcingulate'. You should put 'posteriorcingulate' in the splitfile. Eg, to divide it into three segments, the following line should appear in the splitfile:

posteriorcingulate 3

Only labels that should be split need be specified in the splitfile.

The second method is to specify an area threshold (in mm^2) as the fourth argument, in which case each label is divided until each subdivision is smaller than the area threshold.

The output label name will be the original name with _divN appended, where N is the division number. N will go from 2 to the number of divisions. The first division has the same name as the original label.

Sample Split File

In this example, I have carried out a group analysis and did a cluster-size threshold. This saved a number of clusters in two .annot files, the smallest of which being about 100 mm2