Documentation:Nightly:Registration:RegistrationLibrary:RegLib C09

Back to Registration Library

Slicer Registration Library Case #9: fMRI to structural MRI

Description

This shows registration of an fMRI baseline scan to a structural reference T1. The fMRI suffers from low SNR and bad distortions we seek to correct. Approach: we first compute a rigid alignment of the fMRI to the T1. We then expand this to a nonrigid transform, which we constrain either by explicit limits on the absolute deformation amount or via a mask. Constraints are necessary because the skull present in the T1 would otherwise negatively affect the registration result.

Modules used

• General Registration (BRAINS)
• General Registration (BRAINS)

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.

• RegLib_C09.mrb: input data only, use this to run the tutorial from the start (Slicer mrb file. 40 MB).
• RegLib_C09_full.mrb: includes raw data + all solutions and intermediate files, use to browse/verify (Slicer mrb file. 182 MB).

Keywords

MRI, fMRI, brain, head, intra-subject

Video Screencasts

1. Movie/screencast showing rigid + BSpline registration
2. Movie/screencast showing installation and use of the "SwissSkullStripper" for generating a T1 mask
3. Movie/screencast showing an alternative BSpline registration using masks instead of an explicit deformation limit

Procedure

1. Rigid Registration: register fMRI to T1
   1. open the General Registration (BRAINS) module
   2. Fixed Image Volume: T1 (SPGR)
   3. Moving Image Volume: fMRI
   4. Output Settings:
      1. Slicer BSpline Transform": none
      2. Slicer Linear Transform: create & rename new transform: Xf1_Rigid.tfm
      3. Output Image Volume: none
   5. Registration Phases: check box for Rigid
   6. Main Parameters
      1. Number of Samples: 200,000
   7. click: Apply; runtime < 10 sec (MacPro QuadCore 2.4GHz)
2. Nonrigid BSpline Registration
   1. Fixed Image Volume: T1 SPGR
   2. Moving Image Volume: fMRI
      1. Slicer BSpline Transform": create & rename new transform: Xf2_BSpline.tfm
   3. Slicer Linear Transform: none
      1. Output Image Volume: create & rename new: fMRI_Xf2
   4. Initialization Transform:select "Xf1_Rigid" created above
   5. Initialization Transform Mode: off
   6. Registration Phases: uncheck boxe for Rigid and instead check BSpline
   7. Main Parameters: Number of Samples: 200,000
   8. Advanced Output Settings: Maximum B-Spline Displacement: 3
   9. click: Apply

Alternate Option: Masking

1. if not yet present, install the "SwissSkullStripper" module:

Registration Results

unregistered (click to enlarge)

after rigid registration (click to enlarge)

after nonrigid BSpline registration (click to enlarge) } Acknowledgments