Projects:RegistrationLibrary:RegLib C34

From NAMIC Wiki
Revision as of 17:32, 1 December 2010 by Meier (talk | contribs) (→‎Modules)
Jump to: navigation, search
Home < Projects:RegistrationLibrary:RegLib C34

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 #34: Mouse Brain MRI

Input

MRI baseline lleft MRI follow-up
fixed image 1/target
MRI mouse 1
moving image
MRI mouse 2

Modules

Objective / Background

Aim is to register the brain of a 5 days old mouse (Post Natal Day 5, PND5) with an atlas.

Download

Keywords

MRI, brain, mouse, masking, non-human, non-rigid,

Input Data

  • reference/fixed: DTI baseline: 2.5 x 2.5 x 2.5 mm, 128 x 112 x 44
  • reference/fixed: 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
  • moving : T1 , 1x1x1.1 mm voxel size, 256 x 256 x 193

Registration Challenges

  • The original images have very unusual dynamic range which can cause problems with modules that expect a minimal dynamic range in absolute terms
  • target image content is of lower intensity than the surrounding tissue, making masking essential

Key Strategies

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. 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 IVa: Manual Initial Alignment
  1. the two datasets in this case have very large initial misalignment and also a clipped FOV, which makes it difficult for an automated initializer to find the proper position. This increases the likelihood for the automated registration to go awry. To avoid this we adjust the main (translational) misalignment manually. You may skip this step if your data has closer initial alignment
  2. Go to the Data module, right-click on the "Scene" node and select "Insert Transform Node"
  3. select the new transform node and then rename it to "Xf0_ManualInit_T1-DTI", in the MRML Node Inspector" (press return after typing)
  4. place the DTI_base image in the background and the T1_bc image in the foreground
  5. move the fade slider halfway to the middle so that you can see both volumes
  6. go to the Transforms module, move the translation sliders to move the T1 image on top of the DTI, or enter -20 for both LR and IS translation
  7. in the BRAINSfit module below, select this initialization as starting transform.
  • 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. enter the following settings:
    2. Fixed Image: DTI_base, Moving Image: T1_bc
    3. Registration phases: from Initialize with previously generated transform', select "Xf0_ManualInit" node created before.
    4. Registration phases: check boxes for Include Rigid registration phase, Include Scale Versor3D registration phase Include Affine registration phase. To combine with Phase V below also check the Include BSpline phase box.
    5. Output: Slicer Linear transform: create new, rename to "Xf1_T1-DTI_masked"
    6. Registration Parameters: increase Number Of Samples to 200,000
    7. Control Of Mask Processing Tab: check ROI box, Input Fixed Mask : DTI_mask; Input Moving Mask: T1_mask
    8. Leave all other settings at default
    9. click apply; runtime ca. 1-2 min.
  3. go to Data module. You should see the "T1_bc " 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. enter the following settings:
    2. Fixed Image: DTI_base, Moving Image: T1_bc
    3. Registration phases: from Initialize with previously generated transform', select "Xf1_..." node created before.
    4. Registration phases: check boxes for Include BSpline registration phase
    5. Output: Slicer BSpline transform: create new, rename to "Xf2_T1-DTI_BSpl"
    6. Output Image Volume: create new, rename to "T1_bc_Xf2"; Pixel Type: "short"
    7. Registration Parameters: increase Number Of Samples to 200,000; Number of Grid Subdivisions: 3,3,3
    8. Control Of Mask Processing Tab: check ROI box, Input Fixed Mask : DTI_mask; Input Moving Mask: T1_mask
    9. Leave all other settings at default
    10. click apply; runtime ca. 1-2 min.
  • The result volume "T1_bc_Xf2" is the one we want, since the DTI remains stationary, there is no additional resampling required. Note that the final transform "Xf2_T1-DTI_BSpl" contains both the affine and nonrigid BSpline portions.
  1. To check alignment, go to the Volumes module, select the newly produced DTI volume
  2. under the Display tab, select Color Orientation from the Scalar Mode menu
  3. set DTI as background and the new "T1_bc_Xf2" volume as foreground
  4. Set fade slider to see DTI overlay onto the SPGR image, you should see alignment similar to the result shown below.

Registration Results

Surface model of original unregistered brains original, not registered
original unregistered brains original, not registered
affine registered brains affine registered
registered brains after 5x5x5 nonrigid BSpline alignment registered brains after 5x5x5 nonrigid BSpline alignment]