Differences between revisions 17 and 18
Deletions are marked like this. Additions are marked like this.
Line 1: Line 1:
[[FsTutorial|Back to Top Tutorial Page]]<<BR>>
[[FsTutorial/MultiModal|Back to Multimodal Top]]
Other multimodal tutorials: <<BR>>
[[FsTutorial/MultiModalRegistration|A. Multimodal Registration]],
[[FsTutorial/MultiModalFmriIndividual|B. Individual fMRI Integration]],
[[FsTutorial/MultiModalFmriGroup|C. Surface-based Group fMRI Analysis]]<<BR>>
[[FsTutorial|Back to Top Tutorial Page]]<<BR>> [[FsTutorial/MultiModal|Back to Multimodal Top]] Other multimodal tutorials: <<BR>> [[FsTutorial/MultiModalRegistration|A. Multimodal Registration]], [[FsTutorial/MultiModalFmriIndividual|B. Individual fMRI Integration]], [[FsTutorial/MultiModalFmriGroup|C. Surface-based Group fMRI Analysis]]<<BR>>
Line 8: Line 3:
The purpose of this tutorial is to give you experience with the
integration of Diffusion Tensor Imaging (DTI) with FreeSurfer. The
data were collected at MGH as part of the MIND.
The purpose of this tutorial is to give you experience with the integration of Diffusion Tensor Imaging (DTI) with FreeSurfer. The data were collected at MGH as part of the MIND.
Line 13: Line 6:
It is assumed that you already have knowledge of DTI and its analysis. This paragraph is supplied as a summary for completeness. DTI attempts to measure the diffusion of water in the brain. White matter tracts tend to be like little straws which constrain the direction diffusion. The DTI analysis measures the orientation and the "strength" of the orientation. Eg, CSF will have no strong orientation. This strength is often measured as the fractional anisotropy (FA). The diffusivity is often measured as the Apparent Diffusion Coefficient (ADC).
Line 14: Line 8:
It is assumed that you already have knowledge of DTI and its
analysis. This paragraph is supplied as a summary for completeness.
DTI attempts to measure the diffusion of water in the brain. White
matter tracts tend to be like little straws which constrain the
direction diffusion. The DTI analysis measures the orientation and the
"strength" of the orientation. Eg, CSF will have no strong
orientation. This strength is often measured as the fractional
anisotropy (FA). The diffusivity is often measured as the Apparent
Diffusion Coefficient (ADC).

Prior to DTI analysis, it is customary to motion correct, and
sometimes eddy current correct, the diffusion weighted images. This is
done by selecting one of the images as a template (usually a low-b
volume). This forces all the DTI-related maps to be in-line with the
template. Consequently, this template should be used as the
registration target.
Prior to DTI analysis, it is customary to motion correct, and sometimes eddy current correct, the diffusion weighted images. This is done by selecting one of the images as a template (usually a low-b volume). This forces all the DTI-related maps to be in-line with the template. Consequently, this template should be used as the registration target.
Line 32: Line 11:

This data was acquired at MGH as part of a MIND Consortium study. The
raw data are in dwi.nii.gz. The FreeSurfer anatomical analysis for
this subject is M87102113.v4. Seventy images were acquired (10 low-b
and 60 diffusion weighted). The bvalues and gradient directions are
stored in bvals.dat and bvects.dat, respectively. These data were
analyzed with FreeSurfer's [[dt_recon|dt_recon]] program, which
created the fa.nii, adc.nii, lowb.nii, and register.dat. The data are
located in:
This data was acquired at MGH as part of a MIND Consortium study. The raw data are in dwi.nii.gz. The FreeSurfer anatomical analysis for this subject is M87102113.v4. Seventy images were acquired (10 low-b and 60 diffusion weighted). The bvalues and gradient directions are stored in bvals.dat and bvects.dat, respectively. These data were analyzed with FreeSurfer's [[dt_recon]] program, which created the fa.nii, adc.nii, lowb.nii, and register.dat. The data are located in:
Line 51: Line 22:

As always, check the registration with the template. In this case, the
template is the low-b image:
As always, check the registration with the template. In this case, the template is the low-b image:
Line 60: Line 29:

The registration should be accurate. For more information on
registration, see the [[FsTutorial/MultiModalRegistration|Registration]]
Tutorial. 
The registration should be accurate. For more information on registration, see the [[FsTutorial/MultiModalRegistration|Registration]] Tutorial.
Line 66: Line 32:

View the FA on the subject's anatomical with the white matter
parcellation (wmparc.mgz) as the segmentation:
View the FA on the subject's anatomical with the white matter parcellation (wmparc.mgz) as the segmentation:
Line 72: Line 36:
tkmedit M87102113.v4 orig.mgz -aux brain.mgz \ tkmedit xxxxxxxxxxx.v4 orig.mgz -aux brain.mgz \
Line 75: Line 39:
  -fthresh 0.2 -fmax 1    -fthresh 0.2 -fmax 1
Line 78: Line 42:
Notes:
Line 79: Line 44:
Notes:
 
1. The FA is a values between 0 and 1, so the thresholds are set to 0.2 and 1.
  1. When the window first comes up, there will be a lot of activity
  outside the brain. To remove this, set the 'View->Mask Functional
  Overlay To Aux Volume' flag.
 
1. In the image below, you can see that white matter tends to have a higher FA.
 1. The FA is a values between 0 and 1, so the thresholds are set to 0.2 and 1.
 1. When the window first comes up, there will be a lot of activity
 outside the brain. To remove this, set the 'View->Mask Functional Overlay To Aux Volume' flag.
 1. In the image below, you can see that white matter tends to have a higher FA.
Line 86: Line 49:
{{attachment:fa.cor128.gif}} {{attachment:fa.cor128.jpg|fa.cor128.gif}}
Line 89: Line 52:

As with the [[FsTutorial/MultiModalFmriIndividual|Individual fMRI Analysis]], we
will first resample the FA into the subject's anatomical space:
As with the [[FsTutorial/MultiModalFmriIndividual|Individual fMRI Analysis]], we will first resample the FA into the subject's anatomical space:
Line 101: Line 62:
Notes:
Line 102: Line 64:
Notes:
 
1. Nearest neighbor is used to avoid averaging FA across voxels.
  2. The output fa.anat.mgh will be 256^3, 1mm3.
 1. Nearest neighbor is used to avoid averaging FA across voxels.
 1. The output fa.anat.mgh will be 256^3, 1mm3.
Line 106: Line 67:
Now run mri_segstats getting a report for each of the following ROIs
(index/name pairs are found in the LUT):
Now run mri_segstats getting a report for each of the following ROIs (index/name pairs are found in the LUT):
Line 110: Line 71:
3024 wm-lh-precentral  3024 wm-lh-precentral
Line 112: Line 73:
3021 wm-lh-pericalcarine  3021 wm-lh-pericalcarine
Line 114: Line 75:
   4 Left-Lateral-Ventricle     4 Left-Lateral-Ventricle
Line 116: Line 77:

The wm-lh-precentral, wm-lh-superiortemporal, wm-lh-pericalcarine are
created by the gyral white matter parcellation.
The wm-lh-precentral, wm-lh-superiortemporal, wm-lh-pericalcarine are created by the gyral white matter parcellation.
Line 125: Line 84:
  --seg $SUBJECTS_DIR/M87102113.v4/mri/wmparc.mgz \   --seg $SUBJECTS_DIR/xxxxxxxxx.v4/mri/wmparc.mgz \
Line 131: Line 90:
Click [[FsTutorial/MultiModalFaStats|HERE]] to see the output. The CC has an average FA of about 0.75, gyral parcellations are about 0.4, the left putamen is 0.27, and the ventricle is 0.2. This is as expected because the CC is highly directional with no crossing fibers so we would expect the CC to have the highest FA. The gyral white matter is also directional but has fibers crossing in them, so one expects the FA to be lower than CC. The gray matter (putamen) is still lower. The ventricle has no fibers, so we expect it to have the lowest FA.
Line 132: Line 92:
Click [[FsTutorial/MultiModalFaStats|HERE]] to see the output. The CC
has an average FA of about 0.75, gyral parcellations are about 0.4,
the left putamen is 0.27, and the ventricle is 0.2. This is as
expected because the CC is highly directional with no crossing fibers
so we would expect the CC to have the highest FA. The gyral white
matter is also directional but has fibers crossing in them, so one
expects the FA to be lower than CC. The gray matter (putamen) is still
lower. The ventricle has no fibers, so we expect it to have the lowest
FA.

Other multimodal tutorials:
[[FsTutorial/MultiModalRegistration|A. Multimodal Registration]],
[[FsTutorial/MultiModalFmriIndividual|B. Individual fMRI Integration]],
[[FsTutorial/MultiModalFmriGroup|C. Surface-based Group fMRI Analysis]]<<BR>>
Other multimodal tutorials: [[FsTutorial/MultiModalRegistration|A. Multimodal Registration]], [[FsTutorial/MultiModalFmriIndividual|B. Individual fMRI Integration]], [[FsTutorial/MultiModalFmriGroup|C. Surface-based Group fMRI Analysis]]<<BR>>

Back to Top Tutorial Page
Back to Multimodal Top Other multimodal tutorials:
A. Multimodal Registration, B. Individual fMRI Integration, C. Surface-based Group fMRI Analysis

The purpose of this tutorial is to give you experience with the integration of Diffusion Tensor Imaging (DTI) with FreeSurfer. The data were collected at MGH as part of the MIND.

DTI Basics

It is assumed that you already have knowledge of DTI and its analysis. This paragraph is supplied as a summary for completeness. DTI attempts to measure the diffusion of water in the brain. White matter tracts tend to be like little straws which constrain the direction diffusion. The DTI analysis measures the orientation and the "strength" of the orientation. Eg, CSF will have no strong orientation. This strength is often measured as the fractional anisotropy (FA). The diffusivity is often measured as the Apparent Diffusion Coefficient (ADC).

Prior to DTI analysis, it is customary to motion correct, and sometimes eddy current correct, the diffusion weighted images. This is done by selecting one of the images as a template (usually a low-b volume). This forces all the DTI-related maps to be in-line with the template. Consequently, this template should be used as the registration target.

This Data Set

This data was acquired at MGH as part of a MIND Consortium study. The raw data are in dwi.nii.gz. The FreeSurfer anatomical analysis for this subject is M87102113.v4. Seventy images were acquired (10 low-b and 60 diffusion weighted). The bvalues and gradient directions are stored in bvals.dat and bvects.dat, respectively. These data were analyzed with FreeSurfer's dt_recon program, which created the fa.nii, adc.nii, lowb.nii, and register.dat. The data are located in:


setenv TUTORIAL_DATA $FREESURFER_HOME/subjects/buckner_data/tutorial_subjs
setenv SUBJECTS_DIR $TUTORIAL_DATA
cd $TUTORIAL_DATA/multimodal/dti


Check the Registration

As always, check the registration with the template. In this case, the template is the low-b image:


tkregister2 --mov lowb.nii --reg register.dat --surf


The registration should be accurate. For more information on registration, see the Registration Tutorial.

FA Viewing

View the FA on the subject's anatomical with the white matter parcellation (wmparc.mgz) as the segmentation:


tkmedit xxxxxxxxxxx.v4 orig.mgz -aux brain.mgz \
  -seg wmparc.mgz \
  -reg register.dat -overlay fa.nii \
  -fthresh 0.2 -fmax 1


Notes:

  1. The FA is a values between 0 and 1, so the thresholds are set to 0.2 and 1.
  2. When the window first comes up, there will be a lot of activity

    outside the brain. To remove this, set the 'View->Mask Functional Overlay To Aux Volume' flag.

  3. In the image below, you can see that white matter tends to have a higher FA.

fa.cor128.gif

FA ROI Analysis

As with the Individual fMRI Analysis, we will first resample the FA into the subject's anatomical space:


mri_vol2vol --mov fa.nii \
  --reg register.dat \
  --fstarg --interp nearest \
  --o fa.anat.mgh


Notes:

  1. Nearest neighbor is used to avoid averaging FA across voxels.
  2. The output fa.anat.mgh will be 256^3, 1mm3.

Now run mri_segstats getting a report for each of the following ROIs (index/name pairs are found in the LUT):

 251 CC_Posterior (Posterior Corpus Callosum)
3024 wm-lh-precentral
3030 wm-lh-superiortemporal
3021 wm-lh-pericalcarine
  12 Left-Putamen
   4 Left-Lateral-Ventricle

The wm-lh-precentral, wm-lh-superiortemporal, wm-lh-pericalcarine are created by the gyral white matter parcellation.

Now run mri_segstats:


mri_segstats \
  --seg $SUBJECTS_DIR/xxxxxxxxx.v4/mri/wmparc.mgz \
  --ctab $FREESURFER_HOME/FreeSurferColorLUT.txt \
  --id 251 --id 3021 --id 3024 --id 3030 --id 12 --id 4 \
  --i fa.anat.mgh --sum fa.stats


Click HERE to see the output. The CC has an average FA of about 0.75, gyral parcellations are about 0.4, the left putamen is 0.27, and the ventricle is 0.2. This is as expected because the CC is highly directional with no crossing fibers so we would expect the CC to have the highest FA. The gyral white matter is also directional but has fibers crossing in them, so one expects the FA to be lower than CC. The gray matter (putamen) is still lower. The ventricle has no fibers, so we expect it to have the lowest FA.

Other multimodal tutorials: A. Multimodal Registration, B. Individual fMRI Integration, C. Surface-based Group fMRI Analysis

FsTutorial/MultiModalDtiIndividual_tktools (last edited 2014-01-28 16:27:11 by LouisVinke)