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

From NAMIC Wiki
Jump to: navigation, search
Line 45: Line 45:
 
:#Resample DTI_mask with above nonrigid transform and T1Gd as reference -> yields mask for T1Gd
 
:#Resample DTI_mask with above nonrigid transform and T1Gd as reference -> yields mask for T1Gd
 
:#Register FLAIR and T1Gd to T1_sub (affine only) '''with''' masking
 
:#Register FLAIR and T1Gd to T1_sub (affine only) '''with''' masking
 +
 +
----
 +
 +
[[Projects:RegistrationLibrary:RegLib_C33_v36|For an alternative approach using Slicer 3.6 see here]]
  
 
=== Procedures ===
 
=== Procedures ===

Revision as of 18:58, 16 May 2012

Home < Projects:RegistrationLibrary:RegLib C33

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 Exampe #33: Diffusion Weighted Image Volume: align with structural reference MRI

Input

this is the fixed T1 reference image. All images are aligned into this space lleft this is the FLAIR scan, to be registered with the T1 this is the T1Gd image, to be registered with the T1 lleft this is the DTI tensor image, to be registered with the T1
fixed image/target
T1
moving image
FLAIR
moving image 2a
T1Gd
moving image 2b
DTI tensor

Modules

Objective / Background

This is a typical example of DTI processing. Goal is to align the DTI image with a structural scan that provides accuracte anatomical reference. The DTI contains acquisition-related distortion and insufficient contrast to discern anatomical detail. For treatment planning and evaluation, location of functionally critical fiber tracts relative to the pathology is sought.

Keywords

MRI, brain, head, intra-subject, DTI, DWI

Input Data

  • reference/fixed : T1 axial, 0.5 x 0.5 x 1 , 512 x 512 x 176
  • moving : FLAIR axial
  • moving : T1Gd axial
  • moving: Tensor data of DWI volume, 32 directions, 2.6 mm slice thickness

Overall Strategy

  1. DWI -> DTI, obtain mask +baseline
  2. Subsample T1 to isotropic 1x1x1 mm -> T1_sub
  3. Register DTI_baseline to T1_sub (affine+nonrigid) w/o masking
  4. Resample DTI_mask with above transform to obtain mask for T1_sub
  5. Re-register DTI_baseline to T1_sub (affine+nonrigid) with masking
  6. Resample DTI with result nonrigid transform
  7. Resample DTI_mask with above nonrigid transform and FLAIR as reference -> yields mask for FLAIR
  8. Resample DTI_mask with above nonrigid transform and T1Gd as reference -> yields mask for T1Gd
  9. Register FLAIR and T1Gd to T1_sub (affine only) with masking

For an alternative approach using Slicer 3.6 see here

Procedures

  • Phase I: Load & Convert DICOM
  1. run Dicom2NRRD converter for the DWI Dicom set. Output file: DWI.nrrd
  2. load DWI.nrrd
  3. load T1, FLAIR and T1Gd images via AddVolume...
  4. Subsample T1: The T1 is 512 x 512 x 176 (0.5 x 0.5 x 1) resolution, which is too large as a reference to sample the DTI into. So we generate a smaller, isotropic 1x1x1 version of the T1 that can serve as reference:
    1. Open Filtering/Resample Scalar Volume module
    2. Input: T1, Output: new volume, rename to T1_sub
    3. interpolation: Hamming
    4. Spacing: 1,1,1
  • Phase II: Check initial Alignment and anisotropy
  1. large initial misalignment and high voxel anisotropy can cause problems when registering and resampling the DTI. In such cases it is advisable to resample the DWI first into a more isotropic space and, if necessary, initially realign to a similar orientation as the T1 reference, 'before converting to DTI. In this example DWI is sufficiently aligned and has voxel dimensions of 1 x 1 x 2.6, which is borderline. Anisotropy ratios greater than 1:3 tend to benefit from initial resampling. In this case we proceed without.
  • 'Phase III: convert DWI -> DTI
  1. go to module: Diffusion / Utilities / Diffusion Tensor Estimation
  2. input: DWI, output : DTI, DTI_base and DTI_mask
  3. Input DWI Volume: DWI, Output DTI Volume: create new, rename to "DTI"
  4. Output Baseline: create new, rename to "DTI_base"
  5. Otsu Threshold Mask: create new, rename to "DTI_iso_mask"
  6. check boxes for Remove Islands and Apply Mask
    1. leave defaults
  • Phase IVa: Register DTI-T1 (unmasked)
  1. open Registration : BrainsFit module (presets: Xf1_DTI-T1_unmasked
    1. Registration Phases:
    2. set T1_sub as fixed and DTI_base as moving image
      1. select/check Include Affine registration phase
      2. select/check Include BSpline registration phase
    3. Output Settings:
      1. select a new transform "Slicer BSpline Transform", rename to "Xf1_DTI-T1_unmasked"
      2. select a new volume "Output Image Volume, rename to "DTI_base_Xf1"
    4. Registration Parameters: increase Number Of Samples to 200,000
    5. Registration Parameters: set Number Of Grid Subdivisions to 5,5,3
    6. Leave all other settings at default
    7. click: Apply; (runtime < 1 min 10 sec. on MacPro Quad Core 2.4Ghz)
    8. result: the absence of the skull in the DTI causes significant and incorrect stretch of image in the I-S direction.
    9. we repeat with mask, using the above estimate to produce a mask for the T1
  • Phase IVb: Build mask for T1: resample DTI_mask
  1. go to Filtering/Resample Scalar Vector DWI Volume
    1. Use provided presets or set the following:
    2. input: DTI_mask
    3. reference: T1_sub
    4. output: create new, rename to T1_sub_mask
    5. Transform node: Xf1_DTI-T1_unmasked
    6. Transforms order: check output-to-input
    7. interpolation type: check "nn"
    8. click: Apply
  2. Go to Volumes module
  3. select new "T1_sub_mask: volume, go to Info tab
  4. check the labelmap checkbox
  • Phase IVc: Register DTI-T1 (masked)
  1. open Registration : BrainsFit module (presets: Xf1_DTI-T1_masked)
    1. settings as in Phase IVa above, except:
      1. select a new transform "Slicer BSpline Transform", rename to "Xf2_DTI-T1_masked"
      2. select a new volume "Output Image Volume, rename to "DTI_base_Xf2"
    2. Mask Processing Mode: check "ROI"
    3. Input Fixed Mask: T1_mask
    4. Input Moving Mask: DTI_mask
    5. Apply. Runtime 45 sec.
    6. go to Volumes module to adjust window&level for new result "DTI_base_Xf2"
    7. alignment is now much improved
  • Phase V: 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_sub
    4. Transform Parameters: select transform "Xf2_DTI-T1_masked
    5. check box: output-to-input
    6. Leave all other settings at defaults
    7. Click Apply; runtime ~ 3 min.
  • Phase VIa: Register FLAIR & T1 Gd: obtain masks
  1. we seek affine alignment only for FLAIR and T1Gd. As above, better results are achieved when masking non-brain parenchyma. To obtain masks for the FLAIR and T1Gd we proceed as above by resampling the DT_mask with the Xf2_DTI-T1 transform and the respective image as reference.
  2. go to Filtering/Resample Scalar Vector DWI Volume
    1. input: DTI_mask
    2. reference: FLAIR
    3. output: create new, rename to FLAIR_mask
    4. Transform node: Xf2_DTI-T1_masked
    5. Transforms order: check output-to-input
    6. interpolation type: check "nn"
    7. click: Apply
  3. Go to Volumes module
  4. select new "FLAIR_mask: volume, go to Info tab
  5. check the labelmap checkbox
  • Phase VIb: Mask for T1Gd
  1. ditto: we seek affine alignment only for T1Gd. Identical procedure as for FLAIR above:
  2. go to Filtering/Resample Scalar Vector DWI Volume
    1. input: DTI_mask
    2. reference: T1Gd
    3. output: create new, rename to T1Gd_mask
    4. Transform node: Xf2_DTI-T1_masked
    5. Transforms order: check output-to-input
    6. interpolation type: check "nn"
    7. click: Apply
  3. Go to Volumes module
  4. select new "T1Gd_mask: volume, go to Info tab
  5. check the labelmap checkbox
  • Phase VIc: Register FLAIR & T1Gd
  1. open Registration : BrainsFit module (presets: Xf4_... and Xf5_...)
    1. Registration Phase: check box for "Affine" only
      1. select a new linear transform "Slicer Linear Transform", rename to "Xf4_FLAIR-T1_masked"
    2. Mask Processing Mode: check "ROI"
    3. Input Fixed Mask: T1_mask
    4. Input Moving Mask: FLAIR_mask
    5. Apply. Runtime ~45 sec.
  2. ditto for T1Gd

Registration Results

DTI aligned with T1

Download

Discussion: Registration Challenges

  • The DTI contains acquisition-related distortions (commonly EPI acquisitions) that can make automated registration difficult.

Discussion: Key Strategies

Acknowledgments

This case is the same data as case 40 from the fMRI neurosurgery data set on central.xnat.org