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

From NAMIC Wiki
Jump to: navigation, search
Line 38: Line 38:
 
##manually edit brain mask with '''Editor''' module. required manual fix at frontal and occipital lobe
 
##manually edit brain mask with '''Editor''' module. required manual fix at frontal and occipital lobe
 
#open '''Expert Automated Registration'' module
 
#open '''Expert Automated Registration'' module
##Settings:  
+
##Settings: Fixed Image: fixed, moving image: moving
::*Fixed Image:
+
##Save Transform: ''Xf1_Affine_wmsk''
::*Moving Image:
+
##Initialization: Centers of Mass, Registration: PipelineAffine
::*Resample Image:
+
##Expected offset: 10 Expected Rotation: 0.2
::*Load Transform:
+
##Expected Scale: 0.1 Expected Skew: 0.05
::*Save Transform: ''RegLib_C10_MRI_AtlasSegmentation_Xform_Affine_wmsk''
+
##Fixed Image Mask: "fixed_mask"
::*Initialization: Centers of Mass,  
+
##Affine Max Iteration: 80 Affine Sampling Ratio: 0.05
::*Registration: PipelineAffine
+
#(alternatively automated Affine Registration:  '''Register Images Multires''' also produces good results
::*Expected offset: 10
+
# run '''Fast Deformable B-spline Registration''' module. Settings:  
::*Expected Rotation: 0.2
+
##Iterations: 20 , Grid Size: 9  
::*Expected Scale: 0.1
+
##fixed image: fixed;  moving image: moving
::*Expected Skew: 0.05
+
##Histogram Bins: 50, Spatial Samples: 50000
::*Fixed Image Mask: ''RegLib_C10_MRI_AtlasSegmentation_fixed_mask''
+
## initial transform:  "Xf1_Affine_wmsk"
::*Affine Max Iteration: 80
+
##Output Transform: Xf2_BSpline1
::*Affine Sampling Ratio: 0.05
 
#(alternatively automated Affine Registration:  '''Register Images Multires''' (Slicer 3.5) also produces good results
 
# run '''Deformable B-spline Registration''' module. Settings:  
 
::*Iterations: 20
 
::*Grid Size: 9  
 
::*Histogram Bins: 50,  
 
::*Spatial Samples: 50000,
 
::*initial transform:  ''RegLib_C10_MRI_AtlasSegmentation_Xform_Affine_wmsk''
 
::*Fixed Image: RegLib_C10_MRI_AtlasSegmentation_fixed
 
::*Fixed Image: RegLib_C10_MRI_AtlasSegmentation_moving
 
::*Output Transform: XForm_BSpline1  -> save an output transform to then apply to other atlas data to be brought into alignment.
 
  
 
=== Registration Results===
 
=== Registration Results===

Revision as of 14:36, 12 October 2010

Home < Projects:RegistrationLibrary:RegLib C10

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 #10: Co-registration of probabilistic tissue atlas for subsequent EM segmentation

Input

this is the fixed T1 reference image. lleft
this is the probabilistic tissue atlas, to be registered to the T1
Target Brain Probabilistic Atlas

Modules

Slicer Registration Library Exampe #10: Co-registration of probabilistic tissue atlas for subsequent EM segmentation

Objective / Background

This is an example of sparse atlas co-registration. Not all atlases have an associated reference image that can be used for registration. Because the atlas represents a map of a particular tissue class probability, its contrast differs significantly from the target image.

Keywords

MRI, brain, head, inter-subject, probabilistic atlas, atlas-based segmentation

Input Data

  • fixed : T1w axial, 0.9375 x 0.9375 x 1.5 mm axial, 256 x 256 x 124
  • moving: Probabilistic Tissue atlas, 0.9375 x 0.9375 x 1.5 mm axial, 256 x 256 x 124

Methods

  • Version 1: Expert Automated + Fast BSpline Modules incl. masking
  1. build brain mask for the fixed image only:
    1. open Skull Stripping module.
    2. Settings: 100 iterations, 20 subdivisions.
    3. Ouput: create new volum, rename to "fixed_mask"
    4. Click: Apply
    5. manually edit brain mask with Editor module. required manual fix at frontal and occipital lobe
  2. open 'Expert Automated Registration module
    1. Settings: Fixed Image: fixed, moving image: moving
    2. Save Transform: Xf1_Affine_wmsk
    3. Initialization: Centers of Mass, Registration: PipelineAffine
    4. Expected offset: 10 , Expected Rotation: 0.2
    5. Expected Scale: 0.1 , Expected Skew: 0.05
    6. Fixed Image Mask: "fixed_mask"
    7. Affine Max Iteration: 80 , Affine Sampling Ratio: 0.05
  3. (alternatively automated Affine Registration: Register Images Multires also produces good results
  4. run Fast Deformable B-spline Registration module. Settings:
    1. Iterations: 20 , Grid Size: 9
    2. fixed image: fixed; moving image: moving
    3. Histogram Bins: 50, Spatial Samples: 50000
    4. initial transform: "Xf1_Affine_wmsk"
    5. Output Transform: Xf2_BSpline1

Registration Results

unregistered
after BSpline non-rigid registration

Download

Link to User Guide: How to Load/Save Registration Parameter Presets

Discussion: Registration Challenges

  • Because the atlas represents a probabilistic image (i.e. contains blurring from combining multiple subjects), its contrast differs significantly from the target image.
  • The atlas has strong rotational misalignment that can cause difficulty for automated affine registration.
  • The two images represent different anatomies, a non-rigid registration is required

Discussion: Key Strategies

  • Because of the strong differences in image contrast, Mutual Information is recommended as the most robust metric.
  • masking (skull stripping) is highly recommended to obtain good results for the initial affine alignment. For the 2nd stage BSpline use the full image (i.e. do NOT use the masked version) unless high-quality masks are available for both fixed & moving image. Using the crude mask created for the initial alignment for the BSpline will likely destabilize.
  • because speed is not that critical, we increase the sampling rate for both affine and BSpline registration
  • we also expect larger differences in scale & distortion than with regular structural scans: so we significantly (2x-3x) increase the expected values for scale and skew from the defaults.
  • a good affine alignment is important before proceeding to non-rigid alignment to further correct for distortions.

Acknowledgments

  • dataset provided by Killian Pohl, Ph.D.