Difference between revisions of "Documentation:Nightly:Registration:RegistrationLibrary:RegLib C03"

From Slicer Wiki
Jump to: navigation, search
m (Text replacement - "\[http:\/\/www\.slicer\.org\/slicerWiki\/index\.php\/([^ ]+) ([^]]+)]" to "$2")
 
(22 intermediate revisions by 3 users not shown)
Line 3: Line 3:
 
== Input ==
 
== Input ==
 
{| style="color:#bbbbbb; " cellpadding="10" cellspacing="0" border="0"
 
{| style="color:#bbbbbb; " cellpadding="10" cellspacing="0" border="0"
|[[Image:RegLib_C03_thumb1.png|150px|lleft|this is the fixed T2 reference image. All images are aligned into this space]]  
+
|[[Image:RegLib_C03_thumb1.png|150px|left|this is the fixed T2 reference image. All images are aligned into this space]]  
|[[Image:RegArrow_NonRigid.png|100px|lleft]]  
+
|[[Image:RegArrow_NonRigid.png|100px|left]]  
|[[Image:RegLib_C03_thumb2.png|150px|lleft|this is the DTI Baseline scan, to be registered with the T2]]
+
|[[Image:RegLib_C03_thumb2.png|150px|left|this is the DTI Baseline scan, to be registered with the T2]]
|[[Image:RegLib_C03_thumb3.png|150px|lleft|this is the DTI tensor image, in the same orientation as the DTI Baseline]]
+
|[[Image:RegLib_C03_thumb3.png|150px|left|this is the DTI tensor image, in the same orientation as the DTI Baseline]]
 
|-
 
|-
 
|fixed image 1/target<br>T2
 
|fixed image 1/target<br>T2
Line 14: Line 14:
 
|}
 
|}
  
==Objective / Background ==
+
== Description ==
Goal is to align the DTI image with the structural reference T2 scan that provides accuracte anatomical reference.
+
This is a simple case of diffusion MRI. Goal is to align the Diffusion Tensor Image (DTI) with the structural reference T2 scan that provides accuracte anatomical reference.<br>
 +
'''Approach''': we compute the registration transform using a scalar DTI_baseline scan and then apply this transform to the DTI tensor image. Because reorienting tensor data requires a different form of resampling, we must use a dedicated module for this purpose. Because DTI acquisitions tend to contain strong geometric distortions, we apply a nonrigid transform to the DTI, which makes resampling necessary.
  
 
== Modules used ==
 
== Modules used ==
Line 22: Line 23:
  
 
== Alternate Versions ==
 
== Alternate Versions ==
*this example covers the most basic form of directly registering a DTI + baseline to a T2. There is another (more advanced) version that show how to address additional issues of a strong initial rotation and strong voxel-anisotropy for the raw DWI image acquired.  [[Documentation:Nightly:Registration:RegistrationLibrary:RegLib_C03B|You will find the advanced version here]].
+
*this example covers the most basic form of directly registering a DTI + baseline to a T2. There is another (more advanced) version that show how to address additional issues of a strong initial rotation and strong voxel-anisotropy for the raw DWI image acquired:
 
*[http://na-mic.org/Wiki/index.php/Projects:RegistrationLibrary:RegLib_C03 for the Slicer '''4.1''' version of this case see here]
 
*[http://na-mic.org/Wiki/index.php/Projects:RegistrationLibrary:RegLib_C03 for the Slicer '''4.1''' version of this case see here]
 
*[http://na-mic.org/Wiki/index.php/Projects:RegistrationLibrary:RegLib_C03_v3 for the Slicer '''3.6.3''' version of this case see here]
 
*[http://na-mic.org/Wiki/index.php/Projects:RegistrationLibrary:RegLib_C03_v3 for the Slicer '''3.6.3''' version of this case see here]
  
 
== Download (from NAMIC MIDAS) ==
 
== Download (from NAMIC MIDAS) ==
*[http://slicer.kitware.com/midas3/download/?items=95265,1 '''RegLib_C03_raw.mrb''': raw data only, use this to run the tutorial from the start <small>(Slicer mrb file. 75 MB). </small>]
+
<small>''Why 2 sets of files?  The "input data" mrb includes only the unregistered data to try the method yourself from start to finish. The full dataset includes intermediate files and results (transforms, resampled images etc.). If you use the full dataset we recommend to choose different names for the images/results you create yourself to distinguish the old data from the new one you generated yourself. ''</small>
*[http://slicer.kitware.com/midas3/download/?items=95057,1 '''RegLib_C03.mrb''': includes raw data + all solutions and intermediate files, use to browse/verify <small>(Slicer mrb file. 130 MB). </small>]
+
*[http://slicer.kitware.com/midas3/download/?items=95490 '''RegLib_C03.mrb''': input data only, use this to run the tutorial from the start <small>(Slicer mrb file. 50 MB). </small>]
 +
*[http://slicer.kitware.com/midas3/download/?items=95491 '''RegLib_C03_full.mrb''': includes raw data + all solutions and intermediate files, use to browse/verify <small>(Slicer mrb file. 108 MB). </small>]
 +
 
 +
== Keywords ==
 +
MRI, brain, head, intra-subject, DTI, DWI
  
 
== Video Screencasts ==
 
== Video Screencasts ==
Line 34: Line 39:
  
 
== Procedure ==
 
== Procedure ==
This assumes you have the following: 1) a T2 reference image, 2) a DTI baseline image and  3) the DTI volume (both obtained from the  [http://www.slicer.org/slicerWiki/index.php/Documentation/4.1/Modules/DiffusionTensorEstimation Diffusion Tensor Estimation module]). If you do not have a baseline image, generate a scalar Trace image from the DTI, using the ''Diffusion Tensor Scalar Measurements'' module:
+
This assumes you have the following: 1) a T2 reference image, 2) a DTI baseline image and  3) the DTI volume (both obtained from the  [[Documentation/4.1/Modules/DiffusionTensorEstimation|Diffusion Tensor Estimation module]]). If you do not have a baseline image, generate a scalar Trace image from the DTI, using the ''Diffusion Tensor Scalar Measurements'' module:
#open the [[Documentation/Nightly/Modules/BRAINSFit|''General Registration (BRAINS)'' module]]
+
#'''Compute Registration''': open the [[Documentation/Nightly/Modules/BRAINSFit|''General Registration (BRAINS)'' module]]
 
##''Input Images: Fixed Image Volume'': T2
 
##''Input Images: Fixed Image Volume'': T2
##'''Input Images: Moving Image Volume'': DTI_baseline
+
##''Input Images: Moving Image Volume'': DTI_baseline
 
##''Output Settings'':  
 
##''Output Settings'':  
 
###''Slicer BSpline Transform'' (create new transform, rename to: "Xf1_DTbase-T2_BSpline")
 
###''Slicer BSpline Transform'' (create new transform, rename to: "Xf1_DTbase-T2_BSpline")
Line 44: Line 49:
 
##''Registration Phases'':  select/check ''Rigid'' , ''Rigid+Scale'', ''Affine'', ''BSpline''
 
##''Registration Phases'':  select/check ''Rigid'' , ''Rigid+Scale'', ''Affine'', ''BSpline''
 
##''Main Parameters'':  
 
##''Main Parameters'':  
###increase ''Number Of Samples'' to 300,000
+
###increase ''Percentage of Samples'' to 0.006
 
###set  ''B-Spline Grid Size'' to 7,7,5 (we have lower resolution in the IS-direction (z), hence we set a smaller (5) grid size there)
 
###set  ''B-Spline Grid Size'' to 7,7,5 (we have lower resolution in the IS-direction (z), hence we set a smaller (5) grid size there)
 
##Leave all other settings at default
 
##Leave all other settings at default
 
##click: ''Apply''; runtime < 1 min.
 
##click: ''Apply''; runtime < 1 min.
#Resample DTI: We have now computed the registration transform, but the output volume produced above is a registered version of the baseline, which we need for validation only. To get the actual DTI registered we now apply this transform to the tensor image.
+
#'''Resample DTI:''' We have now computed the registration transform, but the output volume produced above is a registered version of the baseline, which we need for validation only. To get the actual DTI registered we now apply this transform to the tensor image.
#Open the [[Documentation/Nightly/Modules/ResampleDTIVolume|''Resample DTI Volume'']] module (found under: All Modules) (Do not resample a tensor image with other modules, this one must be used to correctly transform the tensor data)
+
##Open the [[Documentation/Nightly/Modules/ResampleDTIVolume|''Resample DTI Volume'']] module (found under: All Modules) (Do not resample a tensor image with other modules, this one must be used to correctly transform the tensor data)
##Input Volume: select DTI
+
###''Input Volume'': select DTI
##Output Volume: select ''create new Diffusion Tensor Volume'',and rename it to ''DTI_Xf1''
+
###''Output Volume'': select ''create new Diffusion Tensor Volume'',and rename it to ''DTI_Xf1''
##Reference Volume: select ''T2''
+
###''Reference Volume'': select ''T2''
##Transform Parameters: select transform node "Xf1_DTI-T2_BSpline", for  ''Deformation Field'': none ; '''check the ''displacement'' checkbox'''
+
###''Transform Parameters'': select transform node "Xf1_DTI-T2_BSpline", for  ''Deformation Field'': none ; '''check the ''displacement'' checkbox'''
##Leave all other settings at defaults
+
###Leave all other settings at defaults
##Click Apply; runtime ~ 2 min.
+
###Click Apply; runtime ~ 2 min.
#set ''T2'' as background and new  ''DTI_Xf1'' volume as foreground
+
##set ''T2'' as background and new  ''DTI_Xf1'' volume as foreground
#fade between back- and foreground to see DTI overlay onto the T2 image. Note that you can also fade via holding the OPTION+CMD keys (mac) + dragging left mouse.
+
##fade between back- and foreground to see DTI overlay onto the T2 image. Note that you can also fade via holding the OPTION+CMD keys (mac) + dragging left mouse.
  
== Registration Results (animated gifs, click to enlarge & play) ==
+
==Registration Results==
[[Image:RegLib_C03_baseline_unregistered.gif|400px]] baseline & T2 before registration <br>
 
[[Image:RegLib_C03_baseline_registered.gif|400px]]  baseline to T2 after affine+nonrigid alignment <br>
 
[[Image:RegLib_C03_DTI_registered.gif|400px]] DTI and T2 before & after registration <br>
 
  
== Keywords ==
+
{|cellpadding="10" cellspacing="0" border="0"
MRI, brain, head, intra-subject, DTI, DWI
+
|[[Image:RegLib_C03_baseline_unregistered.gif|300px]] ||baseline & T2 before registration (click to enlarge)
 +
|-
 +
|[[Image:RegLib_C03_baseline_registered.gif|300px]]||baseline to T2 after affine+nonrigid alignment (click to enlarge)
 +
|-  
 +
|[[Image:RegLib_C03_DTI_registered.gif|300px]]||DTI and T2 before & after registration (click to enlarge)
 +
|}
  
 
== Discussion: Key Strategies ==
 
== Discussion: Key Strategies ==
 
*the strong EPI-based distortions of the DTI image make nonrigid registration necessary
 
*the strong EPI-based distortions of the DTI image make nonrigid registration necessary
 
*initial alignment & overlap is sufficient so that no "initialization" methods are necessary and registration can succeed without.
 
*initial alignment & overlap is sufficient so that no "initialization" methods are necessary and registration can succeed without.
*contrast & initial pose are similar enough for registration to succeed without any masking. However the DTI estimation procedure '''does''' provide an optional mask that is usually very helpful in registering cases with more "distracting" image content.   [[Projects:RegistrationLibrary:RegLib_C03_2| For an example see the extended version of this case here.]]
+
*contrast & initial pose are similar enough for registration to succeed without any masking. However the DTI estimation procedure '''does''' provide an optional mask that is usually very helpful in registering cases with more "distracting" image content.
 
*the DTI in this example is isotropic and hence can be resampled directly. If the DTI contains strong anisotropy of ratios 1:3 or greater, reorienting the DTI can lead to strong artifacts (e.g. in axial direction appear as blue cast in the color orientation view). In that case it is necessary to resample the DWI in the original orientation to an isotropic size before reorienting. It may also be advisable to first reorient the DWI and perform the DTI estimation afterwards.
 
*the DTI in this example is isotropic and hence can be resampled directly. If the DTI contains strong anisotropy of ratios 1:3 or greater, reorienting the DTI can lead to strong artifacts (e.g. in axial direction appear as blue cast in the color orientation view). In that case it is necessary to resample the DWI in the original orientation to an isotropic size before reorienting. It may also be advisable to first reorient the DWI and perform the DTI estimation afterwards.
 
 
=== Acknowledgments ===
 
=== Acknowledgments ===

Latest revision as of 02:21, 27 November 2019

Home < Documentation:Nightly:Registration:RegistrationLibrary:RegLib C03

Back to Registration Library

Slicer Registration Library Case #3: Diffusion Weighted Image Volume: align with structural reference MRI

Input

this is the fixed T2 reference image. All images are aligned into this space
RegArrow NonRigid.png
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 1/target
T2
moving image 2a
DTI baseline
moving image 2b
DTI tensor

Description

This is a simple case of diffusion MRI. Goal is to align the Diffusion Tensor Image (DTI) with the structural reference T2 scan that provides accuracte anatomical reference.
Approach: we compute the registration transform using a scalar DTI_baseline scan and then apply this transform to the DTI tensor image. Because reorienting tensor data requires a different form of resampling, we must use a dedicated module for this purpose. Because DTI acquisitions tend to contain strong geometric distortions, we apply a nonrigid transform to the DTI, which makes resampling necessary.

Modules used

Alternate Versions

Download (from NAMIC MIDAS)

Why 2 sets of files? The "input data" mrb includes only the unregistered data to try the method yourself from start to finish. The full dataset includes intermediate files and results (transforms, resampled images etc.). If you use the full dataset we recommend to choose different names for the images/results you create yourself to distinguish the old data from the new one you generated yourself.

Keywords

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

Video Screencasts

  1. Movie/screencast showing registration+resampling for Case #03

Procedure

This assumes you have the following: 1) a T2 reference image, 2) a DTI baseline image and 3) the DTI volume (both obtained from the Diffusion Tensor Estimation module). If you do not have a baseline image, generate a scalar Trace image from the DTI, using the Diffusion Tensor Scalar Measurements module:

  1. Compute Registration: open the General Registration (BRAINS) module
    1. Input Images: Fixed Image Volume: T2
    2. Input Images: Moving Image Volume: DTI_baseline
    3. Output Settings:
      1. Slicer BSpline Transform (create new transform, rename to: "Xf1_DTbase-T2_BSpline")
      2. Slicer Linear Transform none
      3. Output Image Volume (create new volume, rename to: "DTIbaseline_Xf1"
    4. Registration Phases: select/check Rigid , Rigid+Scale, Affine, BSpline
    5. Main Parameters:
      1. increase Percentage of Samples to 0.006
      2. set B-Spline Grid Size to 7,7,5 (we have lower resolution in the IS-direction (z), hence we set a smaller (5) grid size there)
    6. Leave all other settings at default
    7. click: Apply; runtime < 1 min.
  2. Resample DTI: We have now computed the registration transform, but the output volume produced above is a registered version of the baseline, which we need for validation only. To get the actual DTI registered we now apply this transform to the tensor image.
    1. Open the Resample DTI Volume module (found under: All Modules) (Do not resample a tensor image with other modules, this one must be used to correctly transform the tensor data)
      1. Input Volume: select DTI
      2. Output Volume: select create new Diffusion Tensor Volume,and rename it to DTI_Xf1
      3. Reference Volume: select T2
      4. Transform Parameters: select transform node "Xf1_DTI-T2_BSpline", for Deformation Field: none ; check the displacement checkbox
      5. Leave all other settings at defaults
      6. Click Apply; runtime ~ 2 min.
    2. set T2 as background and new DTI_Xf1 volume as foreground
    3. fade between back- and foreground to see DTI overlay onto the T2 image. Note that you can also fade via holding the OPTION+CMD keys (mac) + dragging left mouse.

Registration Results

RegLib C03 baseline unregistered.gif baseline & T2 before registration (click to enlarge)
RegLib C03 baseline registered.gif baseline to T2 after affine+nonrigid alignment (click to enlarge)
RegLib C03 DTI registered.gif DTI and T2 before & after registration (click to enlarge)

Discussion: Key Strategies

  • the strong EPI-based distortions of the DTI image make nonrigid registration necessary
  • initial alignment & overlap is sufficient so that no "initialization" methods are necessary and registration can succeed without.
  • contrast & initial pose are similar enough for registration to succeed without any masking. However the DTI estimation procedure does provide an optional mask that is usually very helpful in registering cases with more "distracting" image content.
  • the DTI in this example is isotropic and hence can be resampled directly. If the DTI contains strong anisotropy of ratios 1:3 or greater, reorienting the DTI can lead to strong artifacts (e.g. in axial direction appear as blue cast in the color orientation view). In that case it is necessary to resample the DWI in the original orientation to an isotropic size before reorienting. It may also be advisable to first reorient the DWI and perform the DTI estimation afterwards.

Acknowledgments