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

From NAMIC Wiki
Jump to: navigation, search
Line 47: Line 47:
  
 
=== Procedures ===
 
=== Procedures ===
'''with  BrainsFit''' (ca. 3 min):
 
 
*'''Phase I: LOAD DATA'''
 
*'''Phase I: LOAD DATA'''
 
#download example dataset
 
#download example dataset

Revision as of 18:57, 30 August 2010

Home < Projects:RegistrationLibrary:RegLib C29

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 #29: Intra-subject Brain DTI

Input

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

Modules

Objective / Background

This is a classic case of a multi-sequence MRI exam we wish to spatially align to the anatomical reference scan (T1-SPGR). The scan of interest is the DTI image to be aligned for surgical planning/reference.

Keywords

MRI, brain, head, intra-subject, DTI, T1, T2, non-rigid, tumor, surgical planning

Input Data

  • reference/fixed : T1 SPGR , 0.5x0.5x1 mm voxel size, 512 x 512 x 176
  • moving 1: T2 0.5x0.5x1.5 mm voxel size, 512 x 512 x 92
  • moving 2a: DTI baseline: 1 x 1 x 3 mm, 256 x 256 x 41
  • moving 2b: 1 x 1 x 3 mm, 256 x 256 x 41 x 9 (tensor), original: DWI 256 x 256 x 41 x 32 directions

Registration Challenges

  • The DTI sequence (EPI) contains string distortions we seek to correct via non-rigid alignment
  • the DTI baseline is similar in contrast to a T2, albeit at much lower resolution
  • we have different amounts of voxel-anisotropy

Key Strategies

  • Slicer 3.6.1 recommended modules: BrainsFit, Robust Multiresolution Affine
  • we use an 2-step approach: 1) we first co-register the T2 with the SPGR T1. 2) we register the DTI baseline to the T2 (the T1 provides insufficient similarity), 3) combine this with the Affine1 transform obtained from registering the T2 to the T1. affine transform with 12 DOF (rather than a rigid one) to address distortion differences between the two protocols; 4) resample the DTI volume with the new transform


Procedures

  • Phase I: LOAD DATA
  1. download example dataset
  2. load into 3DSlicer 3.6.1 (Load Scene)
  3. To convert the DWI into a DTI: use the Converters / DICOM to NRRD Converter module
  • Phase II:REGISTER DTI TO T2
  1. open Registration : BrainsFit module
    1. Registration Phases:
    2. set T2 as fixed and DTI_baseline as moving image
    3. select/check Include Rigid registration phase
    4. select/check Include Affine registration phase
    5. select/check Include BSpline registration phase
    6. select Initialize with CenterofHeadAlign
    7. select Include Rigid registration phase
    8. select "Include Affine registration phase"
  2. Output Settings:
    1. select a new transform "Slicer BSpline Transform"
    2. select a new volume "Output Image Volume
  3. Registration Parameters: increase Number Of Samples to 200,000
  4. Leave all other settings at default
  5. click apply; runtime ca. 3 min.
  • Phase III: Register T2 to SPGR
  1. open Registration : BrainsFit module
    1. Registration Phases:
    2. select/check Include Rigid registration phase
    3. select/check Include Affine registration phase
    4. select a new transform Output Transform
  2. Registration Parameters: increase Number Of Samples to 200,000
  3. Leave all other settings at default
  4. click apply; runtime ca. 1-2 min.
  • Phase IV: Merge transforms
  1. We now have an affine transform to align the T2 with the SPGR and an BSpline+Affine to align the DTIBaseline with the T2
  2. there are 2 alternatives:
    1. Resample the DTI twice consecutively with each transform
    2. Combine the transforms and resample only once. This is the preferred approach since it reduces interpolation blurring. However this requires editing the transform file as described below:
      1. Select Save, select the BSpline and the Affine transforms to be saved
      2. in a Text Editor of your choice (outside Slicer), open the BSpline transform file and delete the BSpline portion Transform 0. Then Save As.. 'under a new name.
      3. Reload the newly edited transform into slicer, named e.g. as DTIBaseline_Affine0 (via Add Data)
      4. Move it inside the other affine transform that aligns the T2 to the SPGR.
      5. Select Harden Transforms'
      6. Save the newly "hardened/combined" transform to a file, e.g. named DTI2T1_AffineOnly.tfm
      7. in a Text Editor of your choice (outside Slicer), open again the BSpline transform file and also the above DTI2T1_AffineOnly.tfm
      8. replace the affine portion in the BSpline file with the affine portion of DTI2T1_AffineOnly.tfm and save under a new name,e.g. DTI2T1_combined.tfm
      9. This is now the new transform that will resample the DTI directly to the SPGR (see phase V below)
  • Phase V: Resample DTI
  1. Load the combined transform (Add Data)
  2. 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_reg2T1
    3. Reference Volume: select SPGR
    4. Transform Parameters: select merged transform DTI2T1_combined.tfm
    5. check box: output-to-input
    6. Leave all other settings at defaults
    7. Click Apply; runtime ca. 3-4 min.
  3. Go to the Volumes module, select the newly produced DTI_reg2T1 volume
  4. under the Display tab, select Color Orientation from the Scalar Mode menu
  5. set SPGR as background and new DTI_reg2T1 volume as foreground
  6. Set fade slider to see DTI overlay onto the SPGR image

for more details see the tutorial(s) under Downloads

Registration Results

Unregistered Data + segmentation labelmap unregistered

Conclusion: BrainsFit provides good registration but requires intermediate step of registering to T2.

Download