Difference between revisions of "DTI"

From Pheno Wiki
Jump to: navigation, search
(Scripts)
(Quality Control)
 
(25 intermediate revisions by 2 users not shown)
Line 1: Line 1:
 
__TOC__
 
__TOC__
=Documentation=
 
All CNP DTI processing and quality assurance is documented in the [https://spreadsheets.google.com/ccc?key=0AhLKRRgAIOCVdG9rY0szNFppcGZ0QmF3ejBYLUY5S0E&hl=en#gid=0 CNP DTI QA Google Document]. To view and edit the document, request access from Katie Karlsgodt or Anna Xu.
 
  
 +
<b> Back to [[HTAC]]</b><br/>
 +
<b> Back to [[LA5C]]</b>
  
=Scripts=
+
=Scan Parameters=
<b>dti_proc.sh</b> is a script written by Russ Poldrack to preprocess the raw CNP DTI data. The path for the script is /space/raid2/data/poldrack/CNP/scripts/dti_proc.sh. The usage is dti_proc <input directory path>.
+
The DTI was acquired using a 64 direction sequence. Parameters were: 2mm slices, TR/TE=9000/93, 1 average, 96x96 matrix, 90 degree flip angle, axial slices, b=1000.
 
+
An alternate version exists called dti_proc_regressor. The difference between the two versions is that the regressor script contains a nuisance regressor based on the Galachan, 2010 paper, that can to some extent correct for vibration artifacts. This artifact exists on BMC data before Fall 2010. After Fall 2010 the table was bolted down, correcting the artifact. The CCN table was bolted at installation, thus avoiding the artifact entirely. Unfortunately it does not entirely correct the issue, so the shareable data uses the original version of the script with people who show the artifact marked for elimination.
+
  
 
=Processing=
 
=Processing=
The dti_proc_regressor script processes the data using FSL Diffusion Toolbox (FDT).  
+
==dti_proc script==
Typically, this is followed by TBSS analyses.
+
dti_proc.sh is a script written by Russ Poldrack to preprocess the raw CNP DTI data. The path for the script is /space/raid2/data/poldrack/CNP/scripts/dti_proc.sh.  
  
 +
To run the script on an entire subject group (CONTROLS, SCHZ, BIPOLAR, ADHD) use the wrapper script run_dti_proc.sh
  
=Quality Control=
+
NOTE: An alternate version exists called dti_proc_regressor. The difference between the two versions is that the regressor script contains a nuisance regressor based on the Galachan, 2010 paper, that can to some extent correct for vibration artifacts. This artifact exists on BMC data before Fall 2010. After Fall 2010 the table was bolted down, correcting the artifact. The CCN table was bolted at installation, thus avoiding the artifact entirely. Unfortunately it does not entirely correct the issue, so the shareable data uses the original version of the script with people who show the artifact marked for elimination.
The checking procedures for CNP follow the Cannon Lab DTI QA Protocol, which was adapted from (((Insert links to documents and papers))). This can be found at: [[Cannon_Lab_Processing_and_QA]]
+
  
The script is located at /space/raid2/data/poldrack/CNP/scripts/DTI_QA.sh
+
===Usage===
 +
For example, to run the script on one CNP subject,<br/>
 +
1. Log on to func and go to the directory /space/raid2/data/poldrack/CNP/scripts:
 +
:> ssh $username@funcserv1
 +
:> cd /space/raid2/data/poldrack/CNP/scripts
  
 +
2. Find the path of the raw DTI directory for one subject, and run the script, e.g.:
 +
:> dti_proc /space/raid2/data/poldrack/CNP/SCHZ/CNP_50006B/raw/DTI_64DIR_3 &
  
 +
===Script Actions===
 +
The dti_proc script processes the data using FSL Diffusion Toolbox (FDT).
  
Return to [[DTI_Resources]]
+
Steps:<br/>
 +
1. B0 image is skull stripped, creating B0_image_brain<br/>
 +
2. the raw data are registered to the first (B0) image using mcflirt- this both corrects for eddy currents and helps motion, the resulting file is dti_mcf<br/>
 +
3. dtifit is run, creating FA, L1, L2, MD and other processed images in subject space<br/>
 +
4. the FA and color maps are registered to MNI space (FA2std and V12std).<br/>
 +
5. Generic ROIs from the JHU atlas are applied to the FA image in standard space. <br/>
 +
:WARNING- the resulting values are NOT to be used as data, rather just as first pass markers for scan integrity. They are not well registered enough to be used for anything, as evidenced by the low FA values.<br/>
 +
6. Motion is calculated<br/>
 +
7. bvecs and bvals for each subject are compared to the CNP standards. <br/>
 +
:NOTE: there are different bvecs and bvals for subjects on the CCN and BMC scanners, as well as for a subset of subjects scanned during the transition and intial set up of the CCN<br/>
 +
7. a pdf is generated (dti_diag.pdf)<br/>
  
Return to [[CNP]]
+
===Output===
 +
For the subject indicated, the script will output the following files in the raw DTI directory:
 +
:B0_image_brain_mask.nii.gz
 +
:B0_image_brain.nii.gz
 +
:B0_image.nii.gz
 +
:dti_diag.pdf
 +
:dtifit_FA2std.nii.gz
 +
:dtifit_FA.nii.gz
 +
:dtifit_L1.nii.gz
 +
:dtifit_L2.nii.gz
 +
:dtifit_L3.nii.gz
 +
:dtifit.log
 +
:dtifit_MD.nii.gz
 +
:dtifit_MO.nii.gz
 +
:dtifit_SO.nii.gz
 +
:dtifit_V12std.nii.gz
 +
:dtifit_V1.nii.gz
 +
:dtifit_V2.nii.gz
 +
:dtifit_V3.nii.gz
 +
:dti_mcf.par
 +
:dti_proc.log
 +
:FA2std.mat
 +
 
 +
=Quality Control=
 +
A complete description of the QA can be found at [[DTI_Quality_Control]]

Latest revision as of 11:41, 14 June 2012

Back to HTAC
Back to LA5C

Scan Parameters

The DTI was acquired using a 64 direction sequence. Parameters were: 2mm slices, TR/TE=9000/93, 1 average, 96x96 matrix, 90 degree flip angle, axial slices, b=1000.

Processing

dti_proc script

dti_proc.sh is a script written by Russ Poldrack to preprocess the raw CNP DTI data. The path for the script is /space/raid2/data/poldrack/CNP/scripts/dti_proc.sh.

To run the script on an entire subject group (CONTROLS, SCHZ, BIPOLAR, ADHD) use the wrapper script run_dti_proc.sh

NOTE: An alternate version exists called dti_proc_regressor. The difference between the two versions is that the regressor script contains a nuisance regressor based on the Galachan, 2010 paper, that can to some extent correct for vibration artifacts. This artifact exists on BMC data before Fall 2010. After Fall 2010 the table was bolted down, correcting the artifact. The CCN table was bolted at installation, thus avoiding the artifact entirely. Unfortunately it does not entirely correct the issue, so the shareable data uses the original version of the script with people who show the artifact marked for elimination.

Usage

For example, to run the script on one CNP subject,
1. Log on to func and go to the directory /space/raid2/data/poldrack/CNP/scripts:

> ssh $username@funcserv1
> cd /space/raid2/data/poldrack/CNP/scripts

2. Find the path of the raw DTI directory for one subject, and run the script, e.g.:

> dti_proc /space/raid2/data/poldrack/CNP/SCHZ/CNP_50006B/raw/DTI_64DIR_3 &

Script Actions

The dti_proc script processes the data using FSL Diffusion Toolbox (FDT).

Steps:
1. B0 image is skull stripped, creating B0_image_brain
2. the raw data are registered to the first (B0) image using mcflirt- this both corrects for eddy currents and helps motion, the resulting file is dti_mcf
3. dtifit is run, creating FA, L1, L2, MD and other processed images in subject space
4. the FA and color maps are registered to MNI space (FA2std and V12std).
5. Generic ROIs from the JHU atlas are applied to the FA image in standard space.

WARNING- the resulting values are NOT to be used as data, rather just as first pass markers for scan integrity. They are not well registered enough to be used for anything, as evidenced by the low FA values.

6. Motion is calculated
7. bvecs and bvals for each subject are compared to the CNP standards.

NOTE: there are different bvecs and bvals for subjects on the CCN and BMC scanners, as well as for a subset of subjects scanned during the transition and intial set up of the CCN

7. a pdf is generated (dti_diag.pdf)

Output

For the subject indicated, the script will output the following files in the raw DTI directory:

B0_image_brain_mask.nii.gz
B0_image_brain.nii.gz
B0_image.nii.gz
dti_diag.pdf
dtifit_FA2std.nii.gz
dtifit_FA.nii.gz
dtifit_L1.nii.gz
dtifit_L2.nii.gz
dtifit_L3.nii.gz
dtifit.log
dtifit_MD.nii.gz
dtifit_MO.nii.gz
dtifit_SO.nii.gz
dtifit_V12std.nii.gz
dtifit_V1.nii.gz
dtifit_V2.nii.gz
dtifit_V3.nii.gz
dti_mcf.par
dti_proc.log
FA2std.mat

Quality Control

A complete description of the QA can be found at DTI_Quality_Control