Difference between revisions of "DTI Pipeline"

Latest revision as of 16:43, 29 May 2012

Diffusion tensor imaging (DTI) analyzes the tissues that have an internal fibrous structure, which is analogous to the anisotropy of some crystals (e.g.neural axons of white matter or muscle fibers in the heart). In those tissues, the diffusion of water displays anisotropy in certain directions. That means that molecules of water diffuse more rapidly in the direction aligned with the fibrous structure, and more slowly in the direction perpendicular to it. Thus, the measured rate of diffusion will differ depending on the fiber direction and an observer's point of view.

In diffusion-weighted imaging (DWI), at least 3 gradients are applied in 3 different directions, which is sufficient to estimate the trace of the diffusion tensor. From the diffusion tensor, diffusion anisotropy measures such as the fractional anisotropy (FA) eigenvalues and eigenvectors can be computed and displayed. The principal direction of the diffusion tensor can be also used to visualize the white matter connectivity of the brain. DTI has been proven to be very useful to diagnose vascular strokes in the brain and in other clinical applications.

Note: this diffusion model is a rather simplified model of the diffusion process. It assumes that diffusion within each image voxel is linear and homogeneous.

Introduction

The DTI pipeline of MIPAV prepares diffusion weighted images (DWIs) and computes voxel-wise diffusion tensors (DT) for the further analysis of diffusion tensor imaging (DTI) data, see MIPAV DTI Color Display. The pipeline computes maps of diffusion eigenvalues and eigenvectors. It also determines an anatomical correspondence between DTI and structural MRI images (T2) of the same sample.

MIPAV DWI pipeline overview

Error creating thumbnail: Invalid thumbnail parameters

MIPAV DTI Pipeline schematic

Brain MRI introduction

Brain MRI is usually performed with a 1.5- or 3- T MRI machine, with a gradient strength in the range of 20-60 mT/m and a slew rate of 120 T/m/s. Parameters for a single-shot spin-echo echo-planar imaging (EPI) sequence include: a repetition time (TR) of 6000 ms, an echo time (TE) of 100 ms, a field of view (FOV) of 24 cm. These parameters are typically used to obtain 3- to 5-mm axial or coronal sections with a 5-mm intersection gap. The acquisition matrix is usually 96 × 96 with a reconstruction matrix set to 128 × 128. The DWIs are obtained by using 4 linearly increasing b values in 6-7 non-collinear directions (bmax ~ 703-1000 s/mm2). In addition, a T2-weighted (T2W) reference image is obtained without diffusion weighting. Read more: [1], [2].

Determining an anatomical correspondence

between DTI and structural MRI images of the same sample. In MIPAV we use an image-based registration scheme mainly because i) it doesn't require a field map, which is usually not available for DWIs, and ii) it allows one to correct for artifacts produced by a patient motion. For more information, refer to [3], [4].

At least one of the following images is required as a reference in DTI pipeline:

T2 image - MIPAV uses the T2 as a reference image because it is usually less distorted and has a higher signal-to-noise ratio (SNR) than DWIs.

B0 - which is a DWI volume with no gradient applied. It can be used instead of T2, if T2 is not available.

MIPAV aligns all other images to the reference image (either T2 or B0 volume) by optimizing the cost function, which represents the measure of how well the images are spatially aligned.

Note: In MIPAV, we use the term "cost function" to refer to the negative cost function.

MIPAV DTI pipeline outline

The user uploads a DWI image and T2 image to the pipeline using the Import Data panel. A DWI image can be acquired from many different MRI scanners (including Philips, Siemens, GE, etc.) and in various formats. MIPAV reads the gradient information from the image header, or from the B-matrix file uploaded by the user. The gradient (or B-matrix) information is then displayed in the Gradient table. For the list of image types and how MIPAV reads the header information, refer to Image Types section.
In the Pre-processing step, the B0 volume in DWI is detected (or entered by user), and then rigidly aligned to the T2 image. The DW image is then registered to rigidly aligned B0 using the Optimized Automatic Registration 3.5 D algorithm. These steps are necessary to perform a motion correction and eddy current distortion correction.
In the EPI Distortion Correction step, MIPAV calculates deformation vector fields for rigidly aligned B0 and T2, which came from the Pre-processing step. MIPAV then uses both: the transformation matrices obtained in the Pre-processing step, and deformation vector field values to create a corrected DWI image.
MIPAV then calculates a diffusion tensor using pre-processing DWI and the gradient/B-value information and a mask image uploaded by the user.
MIPAV uses the tensor information to create a tensor statistics, including ADC, color map, eigenvalue, eigenvector, FA, RA, and volume ratio images.
MIPAV creates a 3D visualization of fiber bundle tracts in the brain's white matter using the information from the previous step. The user can save fiber tracts information as .vtk and .dat files. See also: Image formats descriptions.

Image types

Image file type	Auto Population of Bvals	Auto Population of Gradients	Auto Population of Bmatrix	Philips Gradient Creator Utility
Philips PAR/REC V3 &V4	Yes	Yes	No	Yes
Philips PAR/REC 4.1 &4.2	Yes	Yes	No	Yes
Philips DCM V3 &V4	Yes	Yes	No	Yes
Philips DCM 4.1 &4.2	Yes	Yes	No	Yes
Siemens Mosaic DCM, see also output files	No	No	Yes	No
Nifti w/ Philips Par File	Yes	Yes	No	Yes
GE DCM	Yes	Yes	No	No
Text File Type	Auto Population of Bvals	Auto Population of Gradients	Auto Population of Bmatrix	Philips Gradient Creator Utility
fslBvalGrad.txt	Yes	Yes	No	No
dtiStudioBvalGrad.txt	Yes	Yes	No	No
mipavStandardBvalGrad.txt	Yes	Yes	No	No
dcm2nii.bvec	Yes	Yes	No	No
fslBmatrix.txt	No	No	Yes	No
mipavStandardBmatrix.txt	No	No	Yes	No

MIPAV DTI pipeline interface

Main article: DTI pipeline interface.

The MIPAV DTI Pipeline interface contains 6 tabs:

In the Import Data the user uploads DW and T2 image (optional), and the gradient information (optional). MIPAV could also read the gradient information from the image header. The Phillips Gradient Creator will then calculate the gradient table based on the information provided.
In the Pre-processing tab, the user is asked to specify the parameters required for aligning B0 to T2, and then DWI to the rigidly aligned B0, and also parameters for motion correction and eddy current distortion correction.
In the EPI Distortion Correction, the user enters parameters needed for calculation of deformation vector fields for rigidly aligned B0 and T2 from the Pre-processing step.
In the Tensor Estimation tab, the user uploads the mask image, selects the tensor estimation algorithm and specifies the output options.
In the Tensor Statistics the user can upload the tensor image and specify which kind of images (e.g. ADC, color map, Eigen value, Eigen vector, FA, RA, and Volume Ratio) he/she would like to have as an output.
In the Visualization, the user can upload images needed for a 3D visualization of fiber bundle tracts in the brain's white matter using the provided information. The user can then save fiber tracts information as .vtk and .dat files.

DTI Pipeline Tutorials

Refer to DTI Pipeline Tutorials.

TBD.

Difference between revisions of "DTI Pipeline"

Latest revision as of 16:43, 29 May 2012

Contents

Introduction

MIPAV DWI pipeline overview

Brain MRI introduction

Determining an anatomical correspondence

MIPAV DTI pipeline outline

Image types

MIPAV DTI pipeline interface

DTI Pipeline Tutorials

See also:

Navigation menu

Views

Personal tools

Navigation

Search

Tools

@@ Line 1: / Line 1: @@
-This article is  '''a stub'''. It needs improvement.
+Diffusion tensor imaging (DTI) analyzes the tissues that have an internal fibrous structure, which is analogous to the anisotropy of some crystals (e.g.neural axons of white matter or muscle fibers in the heart). In those tissues, the diffusion of water displays anisotropy in certain directions. That means that molecules of water diffuse more rapidly in the direction aligned with the fibrous structure, and more slowly in the direction perpendicular to it. Thus, the measured rate of diffusion will differ depending on the fiber direction and an observer's point of view.
+In diffusion-weighted imaging (DWI), at least 3 gradients are applied in 3 different directions, which is sufficient to estimate the trace of the diffusion tensor. From the diffusion tensor, diffusion anisotropy measures such as the fractional anisotropy (FA) eigenvalues and eigenvectors can be computed and displayed. The principal direction of the diffusion tensor can be also used [[DTI pipeline interface#VisualizationTab |to visualize the white matter connectivity of the brain]]. DTI has been proven to be very useful to diagnose vascular strokes in the brain and in other clinical applications.
+'''Note:''' this diffusion model is a rather simplified model of the diffusion process. It assumes that diffusion within each image voxel is linear and homogeneous.
 == Introduction ==
-The DTI pipeline of MIPAV prepares [http://www.mr-tip.com/serv1.php?type=db1&dbs=Diffusion%20Tensor%20Imaging diffusion weighted images (DWIs)] and computes voxel-wise diffusion tensors (DT) for the further analysis of diffusion tensor imaging (DTI) data, see [[DTI Color Display| MIPAV DTI Color Display]]. The pipeline computes maps of diffusion eigenvalues and eigenvectors. It also determines an anatomical correspondence between DTI and structural MRI images of the same sample.
+The DTI pipeline of MIPAV prepares [http://www.mr-tip.com/serv1.php?type=db1&dbs=Diffusion%20Tensor%20Imaging diffusion weighted images (DWIs)] and computes voxel-wise diffusion tensors ([http://en.wikipedia.org/wiki/Diffusion_MRI#Diffusion_tensor_imaging DT]) for the further analysis of diffusion tensor imaging (DTI) data, see [[DTI Color Display| MIPAV DTI Color Display]]. The pipeline computes maps of diffusion [http://en.wikipedia.org/wiki/Diffusion_MRI#Measures_of_anisotropy_and_diffusivity eigenvalues] and [http://en.wikipedia.org/wiki/Diffusion_MRI#Measures_of_anisotropy_and_diffusivity eigenvectors]. It also determines an anatomical correspondence between DTI and structural MRI images (T2) of the same sample.
 == MIPAV DWI pipeline overview ==
-This section is correct.
 [[File:DtiPipelineWikiSmall.jpg‎ |400px|thumb|right|MIPAV DTI Pipeline schematic]]
@@ Line 10: / Line 15: @@
 Brain MRI is usually performed with a 1.5- or 3- T MRI machine, with a gradient strength in the range of 20-60 mT/m and a slew rate of 120 T/m/s. Parameters for a single-shot spin-echo echo-planar imaging (EPI) sequence include: a repetition time (TR) of 6000 ms, an echo time (TE) of 100 ms, a field of view (FOV) of 24 cm. These parameters are typically used to obtain 3- to 5-mm axial or coronal sections with a 5-mm intersection gap. The acquisition matrix is usually 96 × 96 with a reconstruction matrix set to 128 × 128. The DWIs are obtained by using 4 linearly increasing b values in 6-7 non-collinear directions (bmax ~ 703-1000 s/mm2). In addition, a T2-weighted (T2W) reference image is obtained without diffusion weighting. Read more: [http://www.ajnr.org/content/26/6/1455.full], [http://radiology.rsna.org/content/217/2/331.full].
-=== Determining an anatomical correspondence between DTI and structural MRI images of the same sample ===
+=== Determining an anatomical correspondence ===
+between DTI and structural MRI images of the same sample. In MIPAV we use [http://ieeexplore.ieee.org/xpls/abs_all.jsp?arnumber=941749&tag=1 an image-based registration scheme] mainly because i) it doesn't require [http://www.mr-tip.com/serv1.php?type=db1&dbs=High%20Field%20MRI a field map], which is usually not available for DWIs,  and ii) it allows one to correct for artifacts produced by a patient motion. For more information, refer to [http://ee.sharif.edu/~miap/Files/Medical%20Image%20Registration.pdf], [http://www.sciencedirect.com/science/article/pii/S0165027007002270].
-[http://ieeexplore.ieee.org/xpls/abs_all.jsp?arnumber=941749&tag=1 MIPAV uses  an image-based registration scheme] mainly because i) it doesn't require [http://www.mr-tip.com/serv1.php?type=db1&dbs=High%20Field%20MRI a field map], which is usually not available for DWIs,  and ii) it allows one to correct for mis-registration produced by a patient motion. For more information, refer to [http://ee.sharif.edu/~miap/Files/Medical%20Image%20Registration.pdf], [http://www.sciencedirect.com/science/article/pii/S0165027007002270].
-At least one of the following images is required as a reference for image processing:
+At least one of the following images is required as a reference in DTI pipeline:
 *'''[http://www.mr-tip.com/serv1.php?type=db1&dbs=T2%20Weighted%20Image T2 image]''' - MIPAV uses the T2 as a reference image because it is usually less distorted and has a higher signal-to-noise ratio (SNR) than DWIs.
-*'''B0''' - is a DWI volume in which no gradient has been applied.
+*'''B0''' - which is a DWI volume with no gradient applied. It can be used instead of T2, if T2 is not available.
-MIPAV aligns all other images to the target image (either T2  or B0 volume) by optimizing [[Cost functions used in MIPAV algorithms | the cost function]], which represents the measure of  how well the images are spatially aligned.
+MIPAV aligns all other images to the reference image (either T2  or B0 volume) by optimizing [[Cost functions used in MIPAV algorithms | the cost function]], which represents the measure of  how well the images are spatially aligned.
-'''Note:''' In MIPAV, we use the term cost function to refer to the negative cost function.
+'''Note:''' In MIPAV, we use the term "cost function" to refer to [[Cost functions used in MIPAV algorithms |the negative cost function]].
 === MIPAV DTI pipeline outline ===
 <ol>
-<li>A user uploads a DWI image and [http://en.wikipedia.org/wiki/Spin-spin_relaxation_time T2 image] to the pipeline using the Import Data panel. A DWI image can be acquired from many different MRI scanners (including [http://www.healthcare.philips.com/us_en/products/mri/ Philips], [http://www.medical.siemens.com/webapp/wcs/stores/servlet/CategoryDisplay~q_catalogId~e_-1~a_categoryId~e_12754~a_catTree~e_100010,1007660,12754~a_langId~e_-1~a_storeId~e_10001.htm Siemens], [http://www.gehealthcare.com/euen/mri/index.html GE], etc.) and in [[Supported Formats | various formats]]. MIPAV reads the gradient information from the image header, or from the B-matrix file uploaded by the user. The gradient (or B-matrix) information is then displayed in the Gradient table. For the list of image types and how MIPAV reads the header information, refer to [[#ImageTypes | Image Types]] section.</li>
+<li>The user [[DTI pipeline interface#ImportData | uploads a DWI image]] and [http://en.wikipedia.org/wiki/Spin-spin_relaxation_time T2 image] to the pipeline using the Import Data panel. A DWI image can be acquired from many different MRI scanners (including [http://www.healthcare.philips.com/us_en/products/mri/ Philips], [http://www.medical.siemens.com/webapp/wcs/stores/servlet/CategoryDisplay~q_catalogId~e_-1~a_categoryId~e_12754~a_catTree~e_100010,1007660,12754~a_langId~e_-1~a_storeId~e_10001.htm Siemens], [http://www.gehealthcare.com/euen/mri/index.html GE], etc.) and in [[Supported Formats | various formats]]. MIPAV reads the gradient information from the image header, or from the B-matrix file uploaded by the user. The gradient (or B-matrix) information is then displayed in the Gradient table. For the list of image types and how MIPAV reads the header information, refer to [[#ImageTypes | Image Types]] section.</li>
-<li>In the Pre-processing step, the B0 slice in DWI is detected, and then rigidly aligned to the T2 image. The DW image is then registered to rigidly aligned B0 using the[[Optimized automatic registration 3D | Optimized Automatic Registration 3.5 D]] algorithm. These steps are necessary to perform a motion correction and eddy current distortion correction.</li>
+<li>In the Pre-processing step, the B0 volume in DWI is detected (or entered by user), and then rigidly aligned to the T2 image. The DW image is then registered to rigidly aligned B0 using the[[Optimized automatic registration 3D | Optimized Automatic Registration 3.5 D]] algorithm. These steps are necessary to perform a motion correction and eddy current distortion correction.</li>
 <li>In the EPI Distortion Correction step, MIPAV calculates deformation vector fields for rigidly aligned B0 and T2, which came from the Pre-processing step. MIPAV then uses both: the transformation matrices obtained in the Pre-processing step, and deformation vector field values to create a corrected DWI image.</li>
-<li>MIPAV then creates a tensor using pre-processing DWI and the gradient/B-value information  and a mask image uploaded by the user.</li>
-<li>MIPAV uses the tensor information to create a whole bunch of images, including ADC, color map, Eigen value, Eigen vector, FA, RA, and Volume Ratio.</li>
+<li>MIPAV then calculates a diffusion tensor using pre-processing DWI and the gradient/B-value information  and a mask image uploaded by the user.</li>
+<li>MIPAV uses the tensor information to create a tensor statistics, including ADC, color map, eigenvalue, eigenvector, FA, RA, and volume ratio images.</li>
 <li>MIPAV creates a 3D visualization of fiber bundle tracts in the brain's white matter using the information from the previous step. The user can save fiber tracts information as <span style="font-family:courier">[[Image formats descriptions#VtkXml |.vtk]]</span> and <span style="font-family:courier">[http://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/Diffusion .dat]</span> files. See also: [[Image formats descriptions]].
 </ol>
@@ Line 82: / Line 91: @@
 |-
 |GE DCM
-|In progress
+|Yes
-|In progress
+|Yes
-|In progress
+|No
-|In progress
+|No
 |-
 !align="left"|Text File Type
@@ Line 132: / Line 141: @@
 == MIPAV DTI pipeline interface ==
+''Main article: [[DTI pipeline interface]].''
 The MIPAV DTI Pipeline interface contains 6 tabs:
-*In the Import Data tab a user uploads DW and T2 image (optional), and the gradient information (optional). Or, MIPAV reads the gradient information from the image header. The Gradient Creator will then calculate the gradient table based on provided information.
+*In the [[DTI pipeline interface#ImportData |Import Data]] the user uploads DW and T2 image (optional), and the gradient information (optional). MIPAV could also read the gradient information from the image header. The Phillips Gradient Creator will then calculate the gradient table based on the information provided.
-*In the Pre-processing tab, the user can specify the parameters required for aligning B0 to T2, and then DWI to the rigidly aligned B0, and also parameters for motion correction and eddy current distortion correction.
+*In the [[DTI pipeline interface#PreProcessing | Pre-processing tab]], the user is asked to specify the parameters required for aligning B0 to T2, and then DWI to the rigidly aligned B0, and also parameters for motion correction and eddy current distortion correction.
-* In the EPI Distortion Correction tab, the user enters parameters need to calculate deformation vector fields for rigidly aligned B0 and T2 from the Pre-processing step.
+* In the [[DTI pipeline interface #EPIDistortionCorrection | EPI Distortion Correction]], the user enters parameters needed for calculation of deformation vector fields for rigidly aligned B0 and T2 from the Pre-processing step.
-* In the Tensor Estimation tab, the user uploads the mask image, selects the tensor estimation algorithm and specifies the output options.
+* In the [[DTI pipeline interface #TensorEstimation |Tensor Estimation tab]], the user uploads the mask image, selects the tensor estimation algorithm and specifies the output options.
-* In the Tensor Statistics tab the user can upload the tensor image and specifies which kind of images (e.g. ADC, color map, Eigen value, Eigen vector, FA, RA, and Volume Ratio.) he/she would like to have as an output.
+* In the [[DTI pipeline interface#TensorStatistics |Tensor Statistics]] the user can upload the tensor image and specify which kind of images (e.g. ADC, color map, Eigen value, Eigen vector, FA, RA, and Volume Ratio) he/she would like to have as an output.
-* In the Visualization tab, the user can  upload images needed for a 3D visualization of fiber bundle tracts in the brain's white matter using the information. The user can then save fiber tracts information as [[Image formats descriptions#VtkXml |.vtk and .dat]] files.
+* In the [[DTI pipeline interface#VisualizationTab |Visualization]], the user can  upload images needed for a 3D visualization of fiber bundle tracts in the brain's white matter using the provided information. The user can then save fiber tracts information as [[Image formats descriptions#VtkXml |.vtk and .dat]] files.
-=== Import Data tab ===
-This tab is for uploading DW and T2 images, calculating B-values and creating B-matrix (Gradient) table. DTI Pipeline reads raw data in all [[Supported Formats | MIPAV supported formats]] and [http://medical.nema.org/ DICOM] files.
-==== Upload DWI Image box ====
-'''DWI Image Browse''' – check this option if you would like to upload your image of interest from your computer.
-'''Use Active DWI Image''' – check this option to use an active image, which is already opened in MIPAV.
-'''Note:''' in MIPAV, [[Opening and loading image files#ActiveImage | an active image]] is the one that has a red frame.
-<div id="UploadBValue"></div>
-==== Uploading B-Value/Gradient File or B-Matrix file ====
-This options allows a user to manually upload B-Value/Gradient File or B-Matrix file. See also: [[#TableOptions | Table Options]].
-'''Gradient vector file:''' is a 3D set of coordinates representing the direction (on the form of a vector) along which a magnetic gradient has been applied. Each DTI pulse sequence has a set of gradient vectors associated with it. This set is required for subsequent analysis.
-'''B-value:''' is a scaling factor used to describe the attenuating effect of the MR signal on the different gradients. A diffusion gradient can be represented as a 3D vector, where the direction of vector is in the direction of diffusion and its length is proportional to the gradient strength. The gradient strength, or more often used ''diffusion weighting parameter'' term, is expressed in terms of the B-value parameter, which is proportional to the product of the square of the gradient strength (q) and the diffusion time interval (b ~ q2 • Δt). Read more: [http://www.mr-tip.com/serv1.php?type=db1&dbs=B-Value MRI TIP database].</blockquote>
-<div id="GradientCreator"></div>
-=== Philips Gradient Creator ===
-[[File:GradientTable.jpg|200px|thumb|right|The Gradient or B-matrix table]]
-The Gradient Creator computes the correct gradient table which is necessary to calculate diffusion tensors. The table lists gradient vectors that describe the diffusion weighting directions, which are later used for analysis and computing diffusion tensor. The Gradient Creator works with data acquired by Philips MRI scanners.
-It uses the source code from [http://godzilla.kennedykrieger.org/~jfarrell/software_web.htm DTI_gradient_table_creator] created by [http://godzilla.kennedykrieger.org/~jfarrell/index.htm Jonathan Farrell, Ph.D.]
-The code that computes the gradient table based on minimal user input. The list of input parameters is show below, see [[#GradientCreatorInput | Gradient Creator input parameters]]. You can also refer to [http://godzilla.kennedykrieger.org/~jfarrell/software_web.htm Jonathan Farrell's web site for additional information].
-The specifications of the gradient table created for the image  taken on Philips MRI scanner(s) are determined by:
-*The 3 coordinate frames defined on the MR scanner: world, anatomical, and image coordinate system. For more information, refer to [http://www.slicer.org/slicerWiki/index.php/Coordinate_systems Slicer WIKI].
-*The rules on how imaging options (including slice orientation, [[#SliceAngulation | slice angulation]], phase encoding direction, etc.) impact the gradient directions.
-*The scanner software release.
-*Any motion correction, if performed.
-<div id="SliceAngulation"></div>
-==== Slice angulation ====
-[[File:SliceAngulationProblem.jpg|200px|thumb|right|Slice angulation: (a) - the image was acquired orthogonal to scanner bore; (b) -  the image was acquired with slice orientation rotated relative to the scanner bore (a roll has been applied in this example). Image courtesy: Chris Rorden, Michael Harms, Jolinda Smith, Nytavia Wallace.]]
-In diffusion tensor imaging (see [http://unfweb.unf-montreal.ca/jdoyon/cours_6032/Journal%20of%20Mag%20Res%20Imaging%202001.pdf]), tensors are constructed by collecting a series of direction-sensitive diffusion images. MRI scanners save these directions along with the images and later they are used to reconstruct the diffusion properties of the images.
-<div id="ScannerImageOrder"></div>
-Depending to the scanner, the vectors are recorded with reference either to the scanner bore, or to the imaging grid. This should not a problem if the images are acquired precisely orthogonal to the scanner bore, because in that case the image and the scanner share the same frame of reference.
-Problems can arise in oblique acquisitions when the image plane is not aligned with the scanner bore. In this situation, it is important that the gradient vectors used in the imaging software are in the same frame of reference as the image. This requires conversion or re-ordering, otherwise we will get angulation errors.
-'''Note:''' these angulation errors have little influence on the DTI parameters, which are invariant to tensor rotation - ADC, MD and FA. However, they do affect calculation of  the eigenvectors of the tensor.
-==== The Gradient or B-Matrix table ====
-In the Gradient or B-Matrix table the B column contains b values, while X Gradient, Y Gradient and Z gradient columns contain the diffusion gradients applied along the x, y, and z axis  - Gx, Gy, and Gz correspondingly.
-<div id="GradientCreatorInput"></div>
-==== Gradient creator input parameters ====
-{|style="color:black; background-color:#F8F8F8;" border="1" cellpadding="10" class="wikitable"
-! align="left"| Gradient Creator Parameters
-! align="left"| Options
-! align="left"| Philips PAR/REC Version
-|-
-|'''Fatshift'''
-|'''R:''' Right<br> '''L:'''  Left<br> '''A:''' Anterior<br> '''P:'''  Posterior<br> '''H:''' Head<br> '''F:'''  Feet
-|V3 & V4 <br> V4.1 & V4.2
-|-
-|'''Jones30 or Kirby'''
-|Specify which MRI scanner was used to acquire DW images. This option works for [http://mri.kennedykrieger.org/usersupport.html KKI scanners] only.
-|V3 & V4 <br> V4.1 & V4.2
-|-
-|<div id="GradientResolution"></div>'''Gradient Resolution'''
-|'''Low:''' used for 8 DWI volumes <br> '''Medium:''' used for 17 DWI volumes<br> '''High:''' used for 34 DWI volumes
-|V3 & V4 <br> V4.1 & V4.2
-|-
-|'''Gradient Overplus'''
-|Yes <br> No
-|V3 & V4 <br> V4.1 & V4.2
-|-
-|'''Philips Release'''
-|Rel_1.5 <br> Rel_1.7 <br> Rel_2.0 <br> Rel_2.1 <br> Rel_2.5 <br> Rel_11.x <br>
-|V3 & V4 <br> V4.1 & V4.2
-|-
-|'''Patient Position'''
-|Head First <br> Feet First
-|V3 & V4
-|-
-|'''Patient Orientation'''
-|'''SP:''' [http://en.wikipedia.org/wiki/Supine_position Supine] <br> '''PR:''' [http://en.wikipedia.org/wiki/Prone_position Prone]<br> '''RD:''' [http://en.wikipedia.org/wiki/Decubitus Right Decubitus]<br> '''LD:''' [http://en.wikipedia.org/wiki/Decubitus Left Decubitus]<br>
-|V3 & V4
-|-
-|<div id="FoldOver"></div>'''Fold Over'''
-|'''AP:''' Anterior-Posterior <br> '''RL:''' Right-Left <br> '''FH:''' Head-Feet/ Superior-Inferior
-|V3 & V4
-|-
-|'''OS (Operating System)'''
-|Windows <br> VMS
-|V3 & V4 <br> V4.1 & V4.2 (for KKI scanners only)
-|-
-|'''Inverted'''
-|No<br> Yes
-|V3 & V4 <br> V4.1 & V4.2 (for KKI scanners only)
-|}
-'''Compute Gradient Table''' - use this button to compute the gradient table based on data you entered into Gradient Input Parameters.
-==== Error messages ====
-"''Gradient Table Creator Image dimensions 17 are not consistent with gradient table choice - exacted 8 dimensions''"  of a similar message means that [[#GradientResolution | Gradient Resolution]] (e.g. Low, Medium, High) did not correspond to the volume number of the image you uploaded to the pipeline. In this example, the “Gradient Resolution” was set to “Low” which corresponds to 8 volumes/dimensions, while the image had 17 volumes. In order to calculate the gradient table, we needed to change the “Gradient Resolution” parameter to “Medium".
-"''Gradient Table Creator R: is not consistent with Anterior-Posterior fold over''" of a similar message means that [[#FoldOver | the Fold Over option]] selected does not match the image specification. In order to solve that problem, select the correct fold over option.
-<div id="TableOptions"></div>
-==== Table options ====
-Table Options is available when a user uploaded B-value or B-matrix file using [[#UploadBValue | the Upload B-Value Gradient File or B-Matrix File option]]. It allows the user to directly edit the values in the the Gradient table and to save the updated table. To edit the table, press the Edit button. To save the table, press Save Table As.
-==== Applying  the gradient table ====
-When a user entered all parameters needed for the Gradient Table Creator, the Apply Table button becomes available.
-==== Specials notes ====
-For more information regarding the importing DW images, refer to Jonathan Farrell's web site - [http://godzilla.kennedykrieger.org/~jfarrell/software_web.htm http://godzilla.kennedykrieger.org/~jfarrell/software_web.htm].
-*If your DTI data was acquired on a Philips MR Scanner with new software (release version 1.5, 1.7, 2.0 or 2.1), the PAR file version may be V4.1.  You can verify this by looking at the top right hand corner of the file for the following tag:  <span style="font-family:courier">Research image export tool V4.1</span>.  The other known versions are V3 or V4.
-*The V4.1 PAR files list the diffusion weighting directions, however, these directions are provided not in the image space (see [[#SliceAngulation | Slice angulation]]).  Specifically, in the case [[#GradientCreatorInput | the Gradient Overplus option]] is set to YES, the directions need to be rotated and corrected for [[#SliceAngulation |slice angulation]].  This is automatically done by the Gradient Table creator.
-*The Gradient Table Creator does not read gradient directions from the PAR file because this information is not available for all PAR files (V3 or V4).  It uses the directions documented in the Philips source code instead.
-<div id="ImportDataOutput"></div>
-==== Import Data tab output files ====
-The gradient and B-matrix table files are the output files from the Import Data tab. They could be saved in a user specified folder.
-'''Note:''' Siemens  Mosaic 3D DICOM files if used as an input, would be saved as 4D volumes based on DWI volumes.
-=== Pre-processing tab ===
-In this tab, a user can specify parameters for running [[Optimized automatic registration 3D]] (or OAR 3.5D) algorithm that will perform the rigid registration  of B0 to T2 and then DWI to rigidly aligned B0. OAR 3.5D is used to correct the spatial mis-registration of DWI volumes originating from both - subject motion and eddy current-induced distortions. The user has a choice to run OAR 3D or OAR 3.5D. For more information about both algorithms, refer to the [[Optimized automatic registration 3D |OAR documentation]].
-'''Note:''' if a used doesn't have a T2 image, the registration could still be performed based on the B0 volume, but the user choice would be limited by OAR 3.5 D algorithm only.
-==== OAR 3.5D input options ====
-*'''Skip Pre-Processing''' - check this box if you would like to skip the entire step.
-*'''Reference DWI volume number''' - usually MIPAV can figure out the reference volume number (this is B0 volume), but in case MIPAV missed it, a user can enter the reference volume number in the space provided.
-*'''Degrees of freedom''' - in MIPAV, this parameter provides [[Degrees of freedom | the number of independent pieces of information]] that are used to calculate the transformation. The user has an option to select either Motion Correction, or Motion Correction + Eddy current. See also [[Degrees of freedom#AffineTransformations |Affine Transformations]].
-*'''Interpolation''' - in this box, a user can select [[Interpolation methods used in MIPAV | the interpolation method]] to use. The list includes following interpolations: Trilinear, B-spline 3rd order, B-spline 4th order, Cubic Lagrangian, Quintic Lagrangian, Heptic Lagrangian, and Windowed sinc.
-*'''Cost Function''' - in this box, a user can specify [[Cost functions used in MIPAV algorithms | a cost function]] that would be used in the registration. The list of available cost functions is as follows:  [[Cost functions used in MIPAV algorithms#CorrelationRatio | Correlation ratio]], [[Cost functions used in MIPAV algorithms#LeastSquares | Least squares]], [[Cost functions used in MIPAV algorithms#NormalizedCrossCorrelation| Normalized cross-correlation]], [[Cost functions used in MIPAV algorithms#NormalizedMutualInformation |Normalized mutual information]].
-==== OAR 3.5D output options ====
-*'''Display Transformed DWI dataset image''' - check this box if you would like to view the transformed DWI dataset.
-*'''Save trans matrices from transformed DWI dataset to directory''' - check this box, if you would like to save transformation matrices to an chosen Output directory.
-*'''Correct gradients after transformation''' - check this box, if you would like MIPAV to correct the initial gradient information. See [[#GradientCreator | Gradient Creator]].
-==== B0 to structural image OAR 3D output options ====
-*'''Display registered B0 to structural image:''' - if this box is checked, the result image would be displayed in a separate window.
-*'''Save registered B0 to structural image trans matrix to directory:''' - if this box is checked, the transformation matrix would be saved in the user specified directory.
-*'''RUN OAR 3.5D''' - press this button to run the algorithm.
-*'''Output dir:''' - here a user can specify the directory that would be used to save output and intermediate files, such as transformation matrices and registered B0 volumes.
-==== OAR output files ====
-The output files are as follows:
-*the B0 volume rigidly registered to the structural (T2) image
-*the transformation matrix from  B0 rigidly registered to the structural image
-*the re-sampled structural (T2) image from the OAR algorithm
-=== EPI Distortion Correction tab ===
-Clinical diffusion MR studies are mostly based on single-shot echo-planar imaging (EPI) acquisitions. This method is very sensitive to static magnetic field inhomogeneities and artifacts, which appear due to imperfection of the gradient waveforms, and eddy currents during the long readout time. These issues are primarily responsible for creating nonlinear geometric distortion along the phase-encoding direction. Artifacts often appear at air-tissue border and also in the images of the ventral portions of the frontal and temporal lobes. They become even more severe with increasing magnetic field. The mis-registration among a set of DWI volumes due to geometric distortions leads in turn to spatial inaccuracies in the derivation of the diffusion tensor, ADC, and FA since they are typically computed on a pixel-by-pixel basis combining all diffusion-weighting directions.
-To correct geometric distortions due to B0 inhomogeneities, MIPAV uses the special image registration technique, where the distorted EPI image is registered to a corresponding anatomically correct MR image with an intensity based least-squares similarity metric and the VABRA modeled deformation field to improve the sensitivity in areas of low EPI signal.
-===== References =====
-*Chen B., Guo H., Song A.W., Correction for Direction-Dependent Distortions in Diffusion Tensor Imaging Using Matched Magnetic Field Maps, Brain Imaging and Analysis Center, Duke University, [[Media: DTIdistortionCorrection.pdf]]
-*Rohde G.K., Student Member, IEEE, Aldroubi A., and Dawant B.M., Senior Member, IEEE, The Adaptive Bases Algorithm for Intensity-Based Nonrigid Image Registration, IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 22, NO. 11, NOVEMBER 2003, [[Media: VabraRegistration2003.pdf‎]]
-*Wu M., Chang L.-C, Walker L., Lemaitre H., Barnett A.S, Marenco S., and Pierpaoli C., Comparison of EPI Distortion Correction Methods in Diffusion Tensor MRI Using a Novel Framework, 2008, [[Media: EpiDistortionCorrection.pdf]]
-==== Computing the deformation field ====
-All fields except "Compute deformation field" are populated automatically with data from the Pre-processing tab. However, a user has an option to manually upload previously saved files.
-*'''Reference DWI volume number''' usually MIPAV can figure out the reference volume number (this is B0 volume). In case MIPAV missed it, the user can enter the reference volume number in the space provided.
-*'''Resampled T2 image''' - this field refers to the resampled T2 image from the Pre-processing tab. In case the field did not get populated automatically, the user has an option to upload the image using the Browse button.
-*'''Registered DWI image''' - this field refers to the registered DWI volume from the Pre-processing tab. In case the field did not get populated automatically, the user has an option to upload the image using the Browse button.
-*'''Output dir:''' - by default, this field refers to the same Output directory which was indicated in the Pre-processing tab. The user has an option to indicate/change the path to the Output directory using the Browse button.
-*'''Display deformation field''' - if this box is checked, MIPAV will display deformation field after calculation.
-*'''Display B0 registered to T2''' - if this box is checked, MIPAV will display B0 registered to T2 calculation.
-*'''Compute deformation field''' - press this button to start computation. MIPAV uses VABRA registration method to compute deformation field. VABRA  stays for Vectorized Adaptive Bases Registration Algorithm - a deformable registration algorithm, which is an intensity-based nonrigid registration algorithm. As an output, it returns a registered image and the applied deformation field. For more information, refer to [http://jist.projects.nitrc.org/docs/IACL/Registration/Volume/MedicAlgorithmVABRA.html VABRA web site].
-'''MIPAV API note:''' for API information regarding MIPAV VABRA algorithm, refer to [http://mipav.cit.nih.gov/documentation/api/gov/nih/mipav/model/algorithms/registration/vabra/VabraAlgorithm.html VABRA API].
-==== Computing EPI distortion correction ====
-After computing the deformation field, MIPAV creates the following files:
-*Deformation field file (in XML format),
-*B0 to T2 transformation matrix (.mtx), and
-*4D transformation matrix,
-They make the input for the EPI distortion correction algorithm.
-*'''Deformation field''' - this field is populated automatically as soon as MIPAV finishes computing deformation field. However, the user could also manually upload the file using the Browse button.
-*'''B0 to T2 transformation matrix''' -  this field is also populated automatically as soon as MIPAV finishes computing deformation field. However, the user could also manually upload the file using the Browse button.
-*'''4D transformation matrix''' - this field is also populated automatically as soon as MIPAV finishes computing deformation field. However, the user could also manually upload the file using the Browse button.
-*'''Display EPI-distortion corrected image''' - by default, this box is checked.
-*'''Compute EPI-distortion correction''' - press this button to start computation. Note that depending your computer capabilities, it might take a considerable time.
-=== Tensor Estimation tab ===
-This module estimates the diffusion tensor model from a given DWI Volume. The result is a tensor volume that can be used to compute different anisotropy measurements, for example Fractional Anisotropy, and also perform fiber tractography.
-* Mask Image - is an approximated mask of the white matter that can be used to filter out the background. The mask image should be uploaded by the user.
-==== Tensor estimation algorithms ====
-The following algorithms (available via the DTI algorithm drop-down menu) are used in the pipeline:
-'''Weighed, noise-reduction''' - this is the default method used in MIPAV DTI pipeline. It is based on the book by S. Mori,  "Introduction to Diffusion Tensor Imaging".
-'''LLMSE''' - is based on [http://iacl.ece.jhu.edu/~bennett/catnap/faq.shtml CATNAP] pipeline.
-The following 3 methods have been adopted from Camino. Camino is a fully-featured toolkit for Diffusion MR processing and reconstruction. For more information refer to: Cook et al (2006).
-'''CAMINO: Linear''' - this algorithm uses unweighted linear least-squares to fit the diffusion tensor to the log measurements. The algorithm is based on the method explained in the  paper by Basser P.J., Mattielo J., and Lebihan D. (1994).
-'''CAMINO: Non-linear''' - this method provides more accurate noise modelling than Linear fitting. The algorithm is based on the method explained in papers by Jones and Basser (2004) and Alexander and Barker (2005).
-'''CAMINO: Restore''' - this algorithm fits the diffusion tensor robustly in the presence of outliers, for more information refer to the  [http://cmic.cs.ucl.ac.uk/camino/index.php?n=Tutorials.DTI#restore CAMINO Tutorial]. The algorithm is based on the method explained in the paper by Chang, Jones and Pierpaoli (2005).
-'''CAMINO: Weighed linear''' - this algorithm uses a weighted linear fit, as originally proposed by [http://onlinelibrary.wiley.com/doi/10.1002/mrm.20283/pdf Jones and Basser (2004)].
-<div id="ReferencesA"></div>
-===== References =====
-*Alexander D.C. and Barker G.J., "Optimal imaging parameters for fibre-orientation estimation in diffusion MRI", NeuroImage, 27, 357-367, 2005,[http://discovery.ucl.ac.uk/182080/ Abstract].
-*Basser P.J., Mattielo J., and Lebihan D., "Estimation of the effective self-diffusion tensor from the NMR spin echo, Journal of Magnetic Resonance", 103, 247-54, 1994, [[Media:Basser1994.pdf], [http://www.ncbi.nlm.nih.gov/pmc/articles/PMC1275686/ Abstract]
-*Chang L-C, Jones D.K. and Pierpaoli C., "RESTORE: Robust estimation of tensors by outlier rejection, Magnetic Resonance in Medicine", 53(5), 1088-1095, 2005, [http://onlinelibrary.wiley.com/doi/10.1002/mrm.20426/pdf PDF]
-*Cook P. A., Bai Y.,  Nedjati-Gilani S.,  Seunarine K. K., Hall M. G., Parker G. J., Alexander D. C., "Camino: Open-Source Diffusion-MRI Reconstruction and Processing", 14th Scientific Meeting of the International Society for Magnetic Resonance in Medicine, Seattle, WA, USA, p. 2759, May 2006, [http://www.cs.ucl.ac.uk/research/medic/camino/files/camino_2006_abstract.pdf Abstract]
-*Jones D.K. and Basser P.J., "Squashing peanuts and smashing pumpkins: How noise distorts diffusion-weighted MR data", Magnetic Resonance in Medicine, 52(5), 979-993, 2004, [http://onlinelibrary.wiley.com/doi/10.1002/mrm.20283/pdf PDF]
-*S. Mori,  "Introduction to Diffusion Tensor Imaging", ISBN-10: 0444528288 | ISBN-13: 978-0444528285, [http://www.sciencedirect.com/science/book/9780444528285 ScienceDirect], [[Media:SMori2008.pdf]]
-==== Output Options ====
-The following output options available for all tensor estimation algorithms:
-*'''Save Tensor image''' - when this box checked, MIPAV saves the tensor image in the Output directory.
-*'''Display Tensor image''' -  when this box checked, MIPAV displays the tensor image.
-The following output options available only for tensor estimation algorithms based on CAMINO code:
-*'''Save Exit Code image''' -  when this box checked, MIPAV saves the tensor image in the Output directory.
-*'''Display Exit Code image''' - when this box checked, MIPAV displays the exit code image. The exit code is set 0 in brain voxels and -1 in background voxels.
-*'''Save Intensity image''' -  when this box checked, MIPAV saves the tensor image in the Output directory.
-*'''Display Intensity image''' - when this box checked, MIPAV displays the intensity image.
-*'''Output directory''' - the directory that is used to save output and intermediate files.
-=== Tensor Statistics tab ===
-* Tensor image - MIPAV populated this field automatically with a path to the tensor image created in the Tensor Estimation tab. However, a user has an option to upload the file using the Browse button.
-==== Output options ====
-The user could use the following check boxes to specify what kind of images they want to create/display.
-{|style="color:black; background-color:#F8F8F8;" border="1" cellpadding="10" class="wikitable"
-!colspan="6"|Output options
-|-
-|Create ADC image
-|Display ADC image
-|-
-|Create Color image
-|Display Color image
-|-
-|Create Eigen Value image
-|Display  Eigen Value image
-|-
-|Create Eigen Vector image
-|Display Eigen Vector image
-|-
-|Create RA image
-|Display RA image
-|-
-|Create FA image
-|Display FA image
-|-
-|Create Trace image
-|Display Trace image
-|-
-|Create VR image
-|Display VR image
-|-
-|}
-==== Output directory ====
-*'''Output dir:''' - in this field, the user could specify the directory that would be used for output files.
-=== Visualization tab ===
-Tensor Image
-Color Image
-EVector Image
-FA Image
-EValue Image
-T2 Image
-Load
+== DTI Pipeline Tutorials ==
+Refer to [[DTI Pipeline Tutorials]].
+TBD.
-'''See also:'''
+== See also: ==
+*[[Understanding MIPAV capabilities]]
+*[[Developing new tools using the API]]
+*[[Using MIPAV Algorithms | Overview of MIPAV Algorithms]]
+*[[Cost functions used in MIPAV algorithms]]
+*[[Interpolation methods used in MIPAV]]
+*[[DTI Color Display]]
+*[[DTI Estimate tensor]]
 *[[Optimized automatic registration 3D]]
 [[Category:Help:Stub]]
 [[Category:Help:Algorithms]]