Difference between revisions of "Projects:RegistrationLibrary:RegLib C30"

From NAMIC Wiki
Jump to: navigation, search
Line 120: Line 120:
  
 
=== Registration Results===
 
=== Registration Results===
[[Image:RegLib_C29_AGif_full.gif|500px|Registered DTI superimposed on SPGR and T2]] registered (cycles show T1 and T2 and color DTI overlay)<br>
+
[[Image:RegLib_C30_AGif_reg.gif‎ |500px|Registered DTI superimposed on T1]] registered<br>
 
<br>
 
<br>

Revision as of 16:44, 29 September 2010

Home < Projects:RegistrationLibrary:RegLib C30

Back to ARRA main page
Back to Registration main page
Back to Registration Use-case Inventory

v3.6.1 Slicer3-6Announcement-v1.png Slicer Registration Library Case #30: Intra-subject Brain DTI

Input

this is the fixed reference image. All images are aligned into this space lleft this is the DTI Baseline scan, to be registered with the T1 this is the DTI tensor image, in the same orientation as the DTI Baseline
fixed image/target
T1
moving image 2a
DTI baseline
moving image 2b
DTI tensor

Modules

Objective / Background

This is a common case of a DTI exam with no T2 available as structural reference and a T1 that has strong field inhomogeneity. We wish to spatially align the DTI to the anatomical reference scan (T1-SPGR).

Download

Keywords

MRI, brain, head, intra-subject, DTI, T1, non-rigid,

Input Data

  • reference/fixed : T1 , 1x1x1.1 mm voxel size, 256 x 256 x 193
  • moving: DTI baseline: 2.5 x 2.5 x 2.5 mm, 128 x 112 x 44
  • moving DTI tensor: : 2.5 x 2.5 x 2.5 mm, 128 x 112 x 44 x 9 (tensor), original: DWI 256 x 256 x 41 x 36 directions

Registration Challenges

  • The DTI sequence (EPI) has a low resolution, a clipped FOV and distortions we seek to correct via non-rigid alignment
  • the DTI baseline is similar in contrast to a T2, but we have only a T1 as structural reference

Key Strategies

  • Slicer 3.6.1 recommended modules: BrainsFit
  • to align the DTI with the T1 we need 2 preprocessing steps: 1. reduce the bias field inhomogeneity in the reference T1 and 2. obtain a skull-stripping / brain mask for the T1
  • The DWI is already isotropic and hence no resampling is required before obtaining the DTI
  • the DTI-T2 registration includes non-rigid deformation to correct for the strong distortions from the EPI acquisition. Because of the nonrigid component a mask of the brain parenchyma helps greatly in obtaining a meaningful transform.
  • The DTI estimation provides an automated mask for the DTI_baseline scan, but we have no mask for the T1. We can either obtain one through separate segmentation or by sending the DTI_mask through an additional registration step. We use the former here.
  • thus the full pipeline is this:
  1. Bias Field Correction of T1 -> T1_bc
  2. Skull Stripping of T1_bc
  3. DWI -> DTI estimation (incl. DTI_base and DT_mask output)
  4. Affine registration of DTI_baseline to T1_bc, unmasked
  5. non-rigid (BSpline) registration of DTI_baseline to T1_bc, masked, using above affine as starting pose
  6. resample DTI with result Affine+BSpline transform

Procedures

  • Phase I: LOAD DATA
  1. download example dataset
  2. load into 3DSlicer 3.6.1 (Load Scene)
Note that the main data package does not contain the DWI. Please load separately from the second Data download package.
  • Phase II: Bias Correction of T1
  1. open Filtering : N4ITK MRI Bias Field Correction module
    1. Input: T1
    2. Output: create new volume, rename to T1_bc
    3. Use presets or set parameters: BSpline grid: 2,2,2 ; Shrink Factor: 4
    4. Click: Apply
    5. the bias-corrected image will be of type 'float'; to save memory convert to 'uint':
  2. Open Cast Image module
    1. Input: T1_bc, Output: T1_bc
    2. Output Type: unsigned short
  • Phase III: Skull Stripping
  1. Open "Skull Skull Stripper for Structural MR: This is an extension module. If you have never used it before you need to load it first via the extension manager
    1. Click on the gear wheel icon to the right of the Modules menu. Select the Skull Stripper module and click download and install. Save prelim. work and then restart Slicer to load
    2. use presets or choose the following settings: Input Volume: T1_bc ; Brain Mask: create new, rename to "T1_mask"; Subdivisions: 20, Dilation Radius: 1
  • Phase IV: Affine Registration
  1. you may be able to skip this step and perform affine + nonrigid registration in one step, however the stepwise approach gives you more insight and control
  2. Go to the BRAINSFit module
    1. use presets Xf1_DTI-T1_masked_Aff1 or enter the following settings:
    2. Fixed Image: T1_bc, Moving Image: DTI_base
    3. Registration phases: check boxes for Include Rigid registration phase, Include Scale Versor3D registration phase Include Affine registration phase
    4. Output: Slicer Linear transform: create new, rename to "Xf1_DTI-T1_masked"
    5. Registration Parameters: increase Number Of Samples to 200,000
    6. Control Of Mask Processing Tab: check ROI box, Input Fixed Mask : T1_mask; Input Moving Mask: DTI_mask
    7. Leave all other settings at default
    8. click apply; runtime ca. 1-2 min.
  3. go to Data module. You will see the "DTI_base" volume moved inside the "Xf1.." transform node
  4. select T1_bc as background and DTI_base as foreground and check the alignment via the toggle switch
  • Phase V: Nonrigid Registration
  1. Go to the BRAINSFit module
    1. use presets Xf2_DTI-T1_BSpl or enter the following settings:
    2. Fixed Image: T1_bc, Moving Image: DTI_base
    3. Registration phases: check boxes for Include BSpline registration phase
    4. Output: Slicer BSpline transform: create new, rename to "Xf2_DTI-T1_BSpl"
    5. Output Image Volume: create new, rename to "DTI_base_Xf2"; Pixel Type: "short"
    6. Registration Parameters: increase Number Of Samples to 200,000; Number of Grid Subdivisions: 3,3,3
    7. Control Of Mask Processing Tab: check ROI box, Input Fixed Mask : T1_mask; Input Moving Mask: DTI_mask
    8. Leave all other settings at default
    9. click apply; runtime ca. 1-2 min.
  • Phase VI: Resample DTI
  1. Open the Resample DTI Volume module (found under: All Modules)
    1. Input Volume: select DTI
    2. Output Volume: select New DTI Volume, rename to DTI_Xf2
    3. Reference Volume: select T1
    4. Transform Parameters: select transform "Xf2_DTI-T1
    5. check box: output-to-input
    6. Leave all other settings at defaults
    7. Click Apply; runtime 3-4 min.
  2. Go to the Volumes module, select the newly produced DTI_Xf2 volume
  3. under the Display tab, select Color Orientation from the Scalar Mode menu
  4. set T1 as background and new DTI_Xf2 volume as foreground
  5. Set fade slider to see DTI overlay onto the SPGR image

for more details see the tutorial(s) under Downloads

Registration Results

Registered DTI superimposed on T1 registered