Difference between revisions of "Projects:SlicerFAQ:DWI Registration"

From NAMIC Wiki
Jump to: navigation, search
m (Text replacement - "http://www.slicer.org/slicerWiki/index.php/" to "https://www.slicer.org/wiki/")
 
(2 intermediate revisions by one other user not shown)
Line 9: Line 9:
 
=== Overall Strategy ===
 
=== Overall Strategy ===
 
#if voxel size ratios are 1.5 or greater:
 
#if voxel size ratios are 1.5 or greater:
##Resample the T2 reference to isotropic voxel size:[http://www.slicer.org/slicerWiki/index.php/Modules:ResampleVolume-Documentation-3.6 ''Resample Volume'' module] or [http://www.slicer.org/slicerWiki/index.php/Modules:ResampleScalarVectorDWIVolume-Documentation-3.6 ''Resample Scalar/Vector/DWI Volume'' module]
+
##Resample the T2 reference to isotropic voxel size:[https://www.slicer.org/wiki/Modules:ResampleVolume-Documentation-3.6 ''Resample Volume'' module] or [https://www.slicer.org/wiki/Modules:ResampleScalarVectorDWIVolume-Documentation-3.6 ''Resample Scalar/Vector/DWI Volume'' module]
##Resample the DWI (not DTI)  to isotropic voxel size: [http://www.slicer.org/slicerWiki/index.php/Modules:ResampleScalarVectorDWIVolume-Documentation-3.6 ''Resample Scalar/Vector/DWI Volume'' module]. Note that depending on the size of the original DWI, the isotropic volumes can become very large files. Make sure you have sufficient RAM and disk space available.
+
##Resample the DWI (not DTI)  to isotropic voxel size: [https://www.slicer.org/wiki/Modules:ResampleScalarVectorDWIVolume-Documentation-3.6 ''Resample Scalar/Vector/DWI Volume'' module]. Note that depending on the size of the original DWI, the isotropic volumes can become very large files. Make sure you have sufficient RAM and disk space available.
#Extract the baseline and DTI from the isotropic DWI:  [http://www.slicer.org/slicerWiki/index.php/Modules:DiffusionTensorEstimation-Documentation-3.6 ''Diffusion Tensor Estimation'' module]. Also select to generate a mask image.
+
#Extract the baseline and DTI from the isotropic DWI:  [https://www.slicer.org/wiki/Modules:DiffusionTensorEstimation-Documentation-3.6 ''Diffusion Tensor Estimation'' module]. Also select to generate a mask image.
 
#Register the DWI baseline scan to the reference. Preferred structural reference is a T2-weighted scan. If no T2 is available use the next closest contrast.  
 
#Register the DWI baseline scan to the reference. Preferred structural reference is a T2-weighted scan. If no T2 is available use the next closest contrast.  
##If initial misalignment (rotation/translation) is significant (more than ~20% of the FOV), consider a preliminary manual alignment first, using the [http://www.slicer.org/slicerWiki/index.php/Modules:Transforms-Documentation-3.6 ''Transforms'' module]. If the final target is a T1 but a T2 is available, use the T2 as reference for the DWI, but then perform a second T2-T1 or T1-T2 registration.
+
##If initial misalignment (rotation/translation) is significant (more than ~20% of the FOV), consider a preliminary manual alignment first, using the [https://www.slicer.org/wiki/Modules:Transforms-Documentation-3.6 ''Transforms'' module]. If the final target is a T1 but a T2 is available, use the T2 as reference for the DWI, but then perform a second T2-T1 or T1-T2 registration.
##try [http://www.slicer.org/slicerWiki/index.php/Modules:BRAINSFit ''BRAINSfit'' module] as first method of choice
+
##try [https://www.slicer.org/wiki/Modules:BRAINSFit ''BRAINSfit'' module] as first method of choice
 
##masking is recommended for better results. We have a mask for the DWI but not the T2 reference. To obtain a mask for the reference
 
##masking is recommended for better results. We have a mask for the DWI but not the T2 reference. To obtain a mask for the reference
 
##*you can either perform a prelim. registration without masking, then apply this transform to the DWI mask
 
##*you can either perform a prelim. registration without masking, then apply this transform to the DWI mask
 
##*you can obtain a mask from another scan (e.g.T1, via the skull stripping or other segmentation tool) and perform a registration to the reference
 
##*you can obtain a mask from another scan (e.g.T1, via the skull stripping or other segmentation tool) and perform a registration to the reference
 
##*you can try and obtain a mask for the reference directly, via  
 
##*you can try and obtain a mask for the reference directly, via  
#Resample the DTI tensor with above transform, using the dedicated  [http://www.slicer.org/slicerWiki/index.php/Modules:ResampleDTIVolume-Documentation-3.6 ''Resample DTIVolume'' module].
+
#Resample the DTI tensor with above transform, using the dedicated  [https://www.slicer.org/wiki/Modules:ResampleDTIVolume-Documentation-3.6 ''Resample DTIVolume'' module].
  
 
=== Procedures ===
 
=== Procedures ===
Line 110: Line 110:
 
#set ''T1'' as background and new  ''DTI_Xf2'' volume as foreground
 
#set ''T1'' as background and new  ''DTI_Xf2'' volume as foreground
 
#Set fade slider to see DTI overlay onto the T2 image
 
#Set fade slider to see DTI overlay onto the T2 image
 +
=== Caveats ===
 +
*In the BRAINSfit registration step
 +
**if images are not already closely aligned, check boxes to include "Rigid" & "ScaleVersor3D" phases
 +
**always check boxes for "Affine" and "BSpline" phases
 +
**remember to select "Slicer BSpline Transform" and "Output Image Volume" in Output Settings
 +
**reduce default grid size settings to 5,5,5 or similar, depending on level of distortion
 +
**increase sample points to 200,000
 +
**registration will take 2-4 minutes, depending on CPU and grid size setting
 +
**upon completion: drag the result volume back out from the result transform node to avoid confusion: since the transform contains a BSpline, moving a volume inside the node will not show an effect
 +
**when using masks: remember that 2 masks for both fixed and moving are required. The moving mask usually comes from the DTI estimation, the fixed mask can be obtained separately or via a prior (unmasked) registration that then allows to resample the DTI mask and thus produce a fixed mask.
 +
*DTI resampling
 +
**if the DTI has voxel anisotropy of ratio 3 or greater, consider isotropic resampling of DWI first (see *[[Projects:RegistrationLibrary:RegLib_C03| case 3]] for details.
 +
**use the dedicated ''Resample DTI Volume'' node (under ''All Modules'')
 +
**in the 'Resample DTI Volume'', make sure to check the "output-to-input" box in the ''Transform Parameters'' box. This will execute the affine transform first and the BSpline second.
 +
**always set a ''Reference Volume''. This will be the T1 or T2 baseline to which the DTI is registered. Check the resolution of the reference: if the image is very large (e.g. 512 x 512 x 100 etc.) you may get memory problems trying to resample the 4D tensor volume. In that case first subsample the reference scan (e.g. to 1x1x1) and the resample the DTI into this new space
 +
**validate the resulting resampled DTI using ''Color Orientation'' in the Volumes/Display settings. Make sure the FOV is not clipped and there is no "color cast" indicating blurring artifacts from anisotropic DTI voxels.

Latest revision as of 17:44, 10 July 2017

Home < Projects:SlicerFAQ:DWI Registration

How do I register a DWI/DTI image to a reference image ?

  • Problem: The DWI/DTI image is not in the same orientation as the reference image that I would like to use to locate particular anatomy; the DWI image is distorted and does not line up well with the structural images
  • Explanation: DWI images are often acquired as EPI sequences that contain significant distortions, particularly in the frontal areas. Also because the image is acquired before or after the structural scans, the subject may have moved in between and the position is no longer the same.
  • Fix: obtain a baseline image from the DWI sequence, register that with the structural image and then apply the obtained transform to the DTI tensor. The two chief issues with this procedure deal with the difference in image contrast between the DWI and the structural scan, and with the common anisotropy of DWI data.
  • you can find example cases of DWI registration in the Slicer Registration Case Library, which includes example datasets and step-by-step instructions. Find an example closest to your scenario and perform the registration steps recommended there.

Overall Strategy

  1. if voxel size ratios are 1.5 or greater:
    1. Resample the T2 reference to isotropic voxel size:Resample Volume module or Resample Scalar/Vector/DWI Volume module
    2. Resample the DWI (not DTI) to isotropic voxel size: Resample Scalar/Vector/DWI Volume module. Note that depending on the size of the original DWI, the isotropic volumes can become very large files. Make sure you have sufficient RAM and disk space available.
  2. Extract the baseline and DTI from the isotropic DWI: Diffusion Tensor Estimation module. Also select to generate a mask image.
  3. Register the DWI baseline scan to the reference. Preferred structural reference is a T2-weighted scan. If no T2 is available use the next closest contrast.
    1. If initial misalignment (rotation/translation) is significant (more than ~20% of the FOV), consider a preliminary manual alignment first, using the Transforms module. If the final target is a T1 but a T2 is available, use the T2 as reference for the DWI, but then perform a second T2-T1 or T1-T2 registration.
    2. try BRAINSfit module as first method of choice
    3. masking is recommended for better results. We have a mask for the DWI but not the T2 reference. To obtain a mask for the reference
      • you can either perform a prelim. registration without masking, then apply this transform to the DWI mask
      • you can obtain a mask from another scan (e.g.T1, via the skull stripping or other segmentation tool) and perform a registration to the reference
      • you can try and obtain a mask for the reference directly, via
  4. Resample the DTI tensor with above transform, using the dedicated Resample DTIVolume module.

Procedures

  • Note: below the step-by-step cookbook from a typical DWI case in the Slicer Registration Case Library. This assumes a T2 reference. If voxel anisotropy is not a concern, you may skip phases I-III.
  • Phase I: resample T2 to isotropic voxel size
    1. original T2 has 0.5 x 0.5 x 3 mm voxel size, which we resample to 1 x 1 x 1 mm
    2. open Filtering / Resample Scalar Volume module
    3. Spacing: enter 1,1,1
    4. Interpolation: check hamming
    5. Input Volume: select T2raw, Output Volume: create new, rename to "T2"
    6. click: Apply
  • Phase II: Manually align DWI_raw to T2
  1. open Data module
  2. in MRML tree, select the "Scene" node, then right click & select "Insert Transform Node"
  3. click on MRML Node Inspector and rename Xform node to "Xf0_DWI_init"
  4. drag DWI_raw into the "Xf0_DWI_init" node
  5. go to "Transforms" module
  6. select T2 as background and DWI_raw as foreground, set fading slider halfway
  7. manually adjust transform to align the two volumes, particularly LR rotation and IS translation
  8. Save
  9. go to Filtering / Resample Scalar/Vector/DWI module
  10. Input Volume: "DWI_raw", Reference Volume: "T2", Output Volume: "create new Diffusion Weighted Volume", rename to "DWI_Xf0"
  11. Transform Node: select "Xf0_DWI_Init"
  12. Click: Apply
  13. this will correct both the initial misalignment and the anisotropy
  14. save new DWI_Xf0 to disk
  15. depending on RAM of your machine, consider deleting the original DWI_raw node
  • Phase III: DWI -> DTI
  1. open "Diffusion Tensor Estimation" module
  2. Input DWI Volume: DWI_Xf0, Output DTI Volume: create new, rename to "DTI"
  3. Output Baseline: create new, rename to "DTI_base"
  4. Otsu Threshold Mask: create new, rename to "DTI_mask"
  5. Estimation Parameters: reduce Otsu Omega Threshold" to 0.3
  6. check boxes for Remove Islands and Apply Mask
  7. Click: Apply
  • Phase IV: Register DTI (unmasked)
  1. open Registration : BrainsFit module (presets: Xf1_DTI-T2_unmasked)
    1. Registration Phases:
    2. set T2 as fixed and DTI_base as moving image
      1. select/check Initialize GeometryCenter Align
      2. select/check Include Rigid registration phase
      3. select/check Include ScaleVersor3D registration phase
      4. select/check Include Affine registration phase
      5. select/check Include BSpline registration phase
    3. Output Settings:
      1. select a new transform "Slicer BSpline Transform", rename to "Xf1_DTI-T2_unmasked"
      2. select a new volume "Output Image Volume, rename to "DT_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.
  • Phase V: Resample DTI_mask (presets: DTI_mask_Xf1)
    • we use the above Xform to produce a mask for the T2.
  1. Open Resample Scalar/Vector/DWI Volume module
    1. Input Volume: DTI_mask; Output volume: create new volume, rename to "DTI_mask_Xf1"
    2. Transform Node: "Xf1_DTI-T2_unmasked"
    3. select/check: output-to-input
    4. Interpolation Type: select: nn
    5. click: Apply
    6. Go to Volumes module, select the new "DTI_mask_Xf1", in the Info tab, check the Labelmap box
  • Phase VI:REGISTER DTI TO T2 with masking
  1. open Registration : BrainsFit module (presets: Xf2_DTI-T2_masked)
    1. set T2_Xf1 as fixed and DTI_baseline as moving image
      1. select/check Include Rigid registration phase
      2. select/check Include ScaleVersor3D registration phase
      3. select/check Include Affine registration phase
      4. select/check Include BSpline registration phase
    2. Output BSpline Transform: create new , rename to "Xf2_DTI-T2_masked"
    3. Output Volume: create new, rename to "DTI_base_Xf2"
    4. Registration Parameters: increase Number Of Samples to 200,000
    5. Registration Parameters: set Number Of Grid Subdivisions to 7,7,5
    6. Control of Mask Processing
      1. select/check: ROI (rightmost box)
      2. Input Fixed Mask: select "DTI_mask_Xf1"
      3. Input Moving Mask: select "DTI_mask"
    7. Leave all other settings at default
    8. click: Apply; runtime < 1 min.
  • Phase VII: 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_Xf2
    3. Reference Volume: select T2
    4. Transform Parameters: select transform "Xf2_DTI-T2_masked
    5. check box: output-to-input
    6. Leave all other settings at defaults
    7. Click Apply; runtime ~ 2 min.
  3. Go to the Volumes module, select the newly produced DTI_Xf2 volume
  4. under the Display tab, select Color Orientation from the Scalar Mode menu
  5. set T1 as background and new DTI_Xf2 volume as foreground
  6. Set fade slider to see DTI overlay onto the T2 image

Caveats

  • In the BRAINSfit registration step
    • if images are not already closely aligned, check boxes to include "Rigid" & "ScaleVersor3D" phases
    • always check boxes for "Affine" and "BSpline" phases
    • remember to select "Slicer BSpline Transform" and "Output Image Volume" in Output Settings
    • reduce default grid size settings to 5,5,5 or similar, depending on level of distortion
    • increase sample points to 200,000
    • registration will take 2-4 minutes, depending on CPU and grid size setting
    • upon completion: drag the result volume back out from the result transform node to avoid confusion: since the transform contains a BSpline, moving a volume inside the node will not show an effect
    • when using masks: remember that 2 masks for both fixed and moving are required. The moving mask usually comes from the DTI estimation, the fixed mask can be obtained separately or via a prior (unmasked) registration that then allows to resample the DTI mask and thus produce a fixed mask.
  • DTI resampling
    • if the DTI has voxel anisotropy of ratio 3 or greater, consider isotropic resampling of DWI first (see * case 3 for details.
    • use the dedicated Resample DTI Volume node (under All Modules)
    • in the 'Resample DTI Volume, make sure to check the "output-to-input" box in the Transform Parameters box. This will execute the affine transform first and the BSpline second.
    • always set a Reference Volume. This will be the T1 or T2 baseline to which the DTI is registered. Check the resolution of the reference: if the image is very large (e.g. 512 x 512 x 100 etc.) you may get memory problems trying to resample the 4D tensor volume. In that case first subsample the reference scan (e.g. to 1x1x1) and the resample the DTI into this new space
    • validate the resulting resampled DTI using Color Orientation in the Volumes/Display settings. Make sure the FOV is not clipped and there is no "color cast" indicating blurring artifacts from anisotropic DTI voxels.