[wiki:FsTutorial top] | [wiki:FsTutorial/FslFeatFreeSurfer previous]

Overlaying FSL Feat statistical maps

In this exercise, the data set of subject bert is used for demonstration. To begin the exercises, first enter the following directory, and then set the current directory to be the default subjects directory using this command:

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

1.0 Overlaying the statistical map onto the bert's orig volume

Use the following command to display the zmap (zstat1.img) overlaid onto the bert's orig volume. It will display the automatic segmentation, and will also set the threshold at z = 2.3:

 tkmedit bert orig.mgz lh.white -aux brain.mgz \
   -overlay ./fbert.feat/stats/zstat1.nii.gz \
   -overlay-reg ./fbert.feat/reg/freesurfer/anat2exf.register.dat \
   -segmentation ${SUBJECTS_DIR}/bert/mri/aparc+aseg.mgz ${FREESURFER_HOME}/FreeSurferColorLUT.txt \
   -fthresh 2.3 -fmid 3.3 -fslope 1

You should see the image below:BR attachment:tkm-zstat1-cor-128.th23.small.jpg When you click or mouse over a voxel, the cortical or subcortical structure that that voxel belongs to will be displayed in the control panel. You can view any of the volumes in the stats dir in this way as well as the clustered maps in the feat directory.

2.0 View statistical maps on bert's surface

To view any of the statistical maps on bert's surface, close the tkmedit GUI (or open a new terminal window) and run:

tksurfer bert lh inflated \
  -overlay ./fbert.feat/stats/zstat1.nii.gz \
  -overlay-reg ./fbert.feat/reg/freesurfer/anat2exf.register.dat \
  -fthresh 2.3 -fmid 3.3 -fslope 1

You should see the image below:BR attachment:tks-zstat1-rh-lat.th23.small.jpg BR To load the cortical parcellation, goto File->Label->Import Annotation, and enter lh.aparc.annot. The parcellation will be opaque. To show only the ouline (as above), View->LabelStyle->Outline.When you click or mouse over a vertex, the control panel will display the name of the cortical structure. You can view any of the volumes in the stats dir in this way as well as the clustered maps in the feat directory. You can also run the tkmedit and tksurfer commands above in separate shells and use the Save-Point/Goto-Point functionality to navigate through the volume and surface.

3.0 Resampling COPEs to common surface

When performing surface-based group analysis of functional data, the COPEs must be sampled to the common surface space (similar to how you must reslice the COPEs from the native functional space into standard space prior to group analysis). The surface sampling can be done in two ways.

First, you can use feat2surf as in:

feat2surf --feat fbert.feat --cope-only

This will create fbert.feat/reg_surf-lh-fsaverage/stats for the left hemi (and another for the right). There will be cope1.nii.gz. This looks like a volume because it is in nifti format, but it is really a surface stored in a volume format (note it's dimensions are 1974 x 1 x 83 = 163842 = number of vertices in fsaverage's surface). You can then concatenate these files from different subjects to perform group analysis with FreeSurfer's mri_glmfit or FSL's randomise or flame.

Alternatively, when you have multiple runs/subjects, you can use mris_preproc, something like:

mris_preproc --out lh.cope1.nii.gz --target fsaverage --hemi lh \
  --iv fbert.feat/stats/cope1.nii.gz fbert.feat/reg/freesurfer/anat2exf.register.dat

This will create lh.cope1.nii.gz. In this case it will have only one frame, but if you had multiple runs or subjects, then you would add more --iv arguments pointing to the other volumes. This can then be used as input mri_glmfit, randomise, or flame.