Differences between revisions 33 and 34
Revision 33 as of 2006-05-11 09:32:00
Size: 6956
Editor: JenniPacheco
Comment:
Revision 34 as of 2006-05-11 09:39:30
Size: 7591
Editor: JenniPacheco
Comment:
Deletions are marked like this. Additions are marked like this.
Line 65: Line 65:

== Viewing Surfaces with Tksurfer ==

Now that you've checked out everything in tkmedit you can close it and begin to inspect the surfaces that are output, for this you will use tksurfer. Tksurfer displays one hemisphere at a time, this exercise will go through visualizing things on the left hemisphere only but everything works the same on the right hemisphere (except in the initial command you should specify rh instead of lh if you want to look at the right hemisphere).

To open tksurfer with the left hemisphere inflated surface of your subject use the following command:

{{{
tksurfer good_output lh inflated
}}}
  • To follow this exercise exactly be sure you've downloaded the tutorial data set before you begin. If you choose not to download the data set you can follow these instructions on your own data, but you will have to substitute your own specific paths and subject names.

Inspection of Freesurfer Output

In this exercise you will be able to visualize and inspect a good output so that you can become familiar with what the end product should look like. This will step you through visual inspection of a variety of output, but is not necessarily the recommended procedure to take when trying to verify each subject. Some steps are only necessary to check when there are problems, however it is a good idea for new users to become familiar with what the expected output look like and how to view it.

First you need to set your SUBJECTS_DIR to the appropriate place:

cd $FREESURFER_HOME/subjects/buckner_data/tutorial_subjs
setenv SUBJECTS_DIR ${PWD} 

Viewing Volumes with Tkmedit

The volumes that are output can be loaded into tkmedit, along with surfaces and the subcortical segmentation. With one command line you can load in the brainmask.mgz and wm.mgz volumes, the rh.white and lh.white surfaces, and the subcortical segmentation.

tkmedit good_output brainmask.mgz lh.white -aux wm.mgz -aux-surface rh.white -segmentation aseg.mgz

You should see a tkmedit window open up to this:BR attachment:good_output1.jpg BR

You are currently looking at the brainmask.mgz (loaded as the main volume) with the surfaces displayed and the aseg.mgz (subcortical segmentation) overlayed. The pial (red line), white (yellow line) and orig (green line) surfaces are all shown. You can toggle between the brainmask.mgz (main volume) and the wm.mgz (aux volume) and you can click on and off the aseg.mgz overlay. To become more familar with the buttons in the Tkmedit Toolbox please read the [wiki:FsTutorial/Tools Freesurfer Tools] section of this tutorial.

The things you can look at while this is loaded in tkmedit:

  • talairch transformation
  • intensity normalization
  • skull strip
  • wm.mgz volume
  • the final surfaces
  • subcortical segmentation

For the first run through you might find it easiest to toggle off the aseg with the attachment:icon_seg_volume.gif button.

Talairach transformation

To load in the talairach go to File -> Transforms -> Load Transform for Main Volume and browse to that talairach.xfm file that should be found in the $FREESURFER_HOME/subjects/buckner_data/tutorial_subjs/good_output/mri/transforms directory.

When you do this you should see this:

attachment:talairach.jpg

Scroll through the slices and switch the views to get a good look at the transform. You will see some minor distortions, but you do not want to see any major rotations, translations, twisting or other means to skew the talairach from the orientation that your volume is in. Likewise, when you are looking at a coronal view of your volume you are also looking at a coronal view of the talairach, you do not want to see a different view of the talairach. When you are done looking at the talairach you can turn it off by going to File -> Transforms -> Unload Transform for Main Volume.

Intensity Normalization

Scroll through the brainmask volume and notice that the intensity is all uniform. You should not see any very bright or very dark spots within the cortex. If you click on any voxel that is in the wm you can see that it has been normalized to an intensity of (or very close to) 110. The voxel intensity is shown in the Tkmedit Toolbox. This check and the following, Skull Strip check, can be done simultaneously since they both require you to look at features on the brainmask volume.

Skull Strip

Scroll through the brainmask volume and notice that there is no skull left in your image, notice also that the cerebellum is still included in the brainmask volume. You should not see any large areas of skull left behind, or any areas of cortex or cerebellum removed from this volume. This check and the previous, Intensity Normalization check, can be done simultaneously since they both require you to look at features on the brainmask volume.

White Matter Volume

Switch to view the auxillary volume by using the attachment:icon_aux_volume.gif button, this will show the wm.mgz volume. This volume is comprised of all the voxels that freesurfer is calling white-matter, shown in shades of gray. These are the voxels that were normalized to an intensity of, or very close to, 110 as described above. The bright white voxels are voxels that have been added to the volume during the automatic editting of the wm volume. These edits fill the entire ventricle and basal ganglian defect. You can alternate between the wm.mgz volume ( attachment:icon_aux_volume.gif button) and the brainmask.mgz volume (attachment:icon_main_volume.gif button) to see how well freesurfer has classified the white matter.

Final Surfaces

Switch back to the main volume by using the attachment:icon_main_volume.gif button, this will show the brainmask.mgz volume. The surfaces that are overlaid are the pial surface (red line), white surface (yellow line), and orig surface (green line). You can toggle them on and off with the surface buttons; attachment:icon_pial_surface.gif for the pial surface, attachment:icon_main_surface.gif for the white surface, and attachment:icon_orig_surface.gif for the orig surface. The orig surface is the "first guess" at the boundary between the white matter and gray matter. After topology fixing and some other steps the white surface is generated. The white surface is the best and final estimation at the boundary between the white matter and gray matter. The white surface is the surface used in all calculations of thickness so it is important that this surface follows the boundary of the white matter accurately. The pial surface is should accurately follow the boundary between the gray matter and the CSF. As you scroll through the slices keep in mind that you are looking at a 2-dimensional rendering of a 3-dimensional image, be sure to look at more than just one view too (i.e., sagittal, coronal and horizontal).

Subcortical Segmentation

Toggle on the subcortical segmentation with the attachment:icon_seg_volume.gif button. This will show the complete segmentation of the subcortical structures. Each structure is labeled with a unique color/number distinction. If you click on a voxel the structures name and number label will be shown in the Tkmedit Toolbox. Scrolling through the slices you will be able to see that everything is labeled, and done so accurately. Sometimes it is easier to see the structures and their boundaries looking in either the sagittal or horizontal view, so be sure to check around in all of them.

Viewing Surfaces with Tksurfer

Now that you've checked out everything in tkmedit you can close it and begin to inspect the surfaces that are output, for this you will use tksurfer. Tksurfer displays one hemisphere at a time, this exercise will go through visualizing things on the left hemisphere only but everything works the same on the right hemisphere (except in the initial command you should specify rh instead of lh if you want to look at the right hemisphere).

To open tksurfer with the left hemisphere inflated surface of your subject use the following command:

tksurfer good_output lh inflated

FsTutorial/OutputData_tktools (last edited 2013-11-01 14:23:49 by MaritzaEbling)