• Blog
  • About
  • Publications
  • Videos
  • Workshops
  • CV
  • Contact Me
  • Hall of Fame
Menu

Andy's Brain Blog

Street Address
City, State, Zip
Phone Number
fMRI, Neuroimaging, and More

fmri, neuroimaging, and more

Andy's Brain Blog

  • Blog
  • About
  • Publications
  • Videos
  • Workshops
  • CV
  • Contact Me
  • Hall of Fame

Creating Masks from the Juelich Atlas

June 2, 2017 Andrew Jahn

The atlases in FSL can be divided into two classes: probabilistic and non-probabilistic. In a probabilistic atlas, the number in each voxel represents the probability that the voxel belongs to a structure. Take the Harvard-Oxford atlas, for example. A voxel at the MNI coordinates 18, -102, 6 has a value of 76. This number means the odds are 76% that the voxel coordinates will be part of the occipital cortex, for a subject taken at random.

In contrast, take a look at a voxel in the frontal areas, which show more variability from subject to subject. At the voxel coordinates 2, 46, 6, the atlas reports that the odds for a random subject are 44% cingulate gyrus and 44% paracingulate gyrus - equal odds. This also shows that structures in a probabilistic atlas can overlap.

A non-probabilistic atlas, on the other hand, does not have any overlap between any of the structures; the voxel values in these atlases represents what structure that voxel belongs to. Take the Juelich Histological DTI-81 atlas for example. Clicking on the MNI coordinates 2, 30, 6 returns a value of 3, which, according to the atlas, means that this voxel belongs to the genu of the corpus callosum. Any voxel with the value 3, as defined by this atlas, belongs to the genu of the corpus callosum; it does not belong to any other structure.

To create masks from a non-probabilistic atlas, add all of the structures from that atlas onto your overlay and note which number corresponds to the structure you want. Assume that we wish to create a mask of the entire corpus callosum from this atlas. We note that the numbers 3, 4, and 5 represent the genu, body, and splenium, respectively. We then highlight the atlas structure in the overlay window and save it as JuelichAtlas.nii.gz.

We then extract only the numbers 3, 4, and 5 by using the -thr and -uthr options of fslmaths:

fslmaths JuelichAtlas.nii.gz -thr 3 -uthr 5 CC.nii.gz

This zeroes out any number below 3 or above 5.

Lastly, we need to binarize the image, again using fslmaths:

fslmaths CC.nii.gz -bin CC_bin.nii.gz

When you overlay this on a template in fslview, you will see that the structures of the corpus callosum all have values of 1. You can then use this mask for an ROI analysis.

← E-Prime Tutorial #3: Duration, Termination, and Pre-ReleaseE-Prime Tutorial #2: Objects and Properties →
Archive
  • September 2024
  • June 2024
  • March 2024
  • January 2024
  • June 2023
  • May 2023
  • April 2023
  • February 2022
  • January 2022
  • January 2021
  • December 2020
  • April 2020
  • December 2018
  • April 2018
  • January 2018
  • December 2017
  • August 2017
  • July 2017
  • June 2017
  • May 2017
  • April 2017
  • March 2017
  • January 2017
  • November 2016
  • August 2016
  • July 2016
  • February 2016
  • January 2016
  • September 2015
  • August 2015
  • July 2015
  • May 2015
  • April 2015
  • December 2014
  • October 2014
  • September 2014
  • August 2014
  • July 2014
  • June 2014
  • May 2014
  • April 2014
  • March 2014
  • February 2014
  • January 2014
  • December 2013
  • November 2013
  • October 2013
  • August 2013
  • July 2013
  • June 2013
  • May 2013
  • April 2013
  • March 2013
  • February 2013
  • January 2013
  • December 2012
  • November 2012
  • October 2012
  • September 2012
  • August 2012
  • July 2012
  • June 2012
  • April 2012
  • March 2012

What's on Andy's Brain this month?

Connect with Andy!

Powered by Squarespace