Documentation:Nightly:Registration:RegistrationLibrary:RegLib C06

From Slicer Wiki
Revision as of 20:33, 2 September 2013 by Meier (talk | contribs) (→‎Input)
Jump to: navigation, search
Home < Documentation:Nightly:Registration:RegistrationLibrary:RegLib C06

Back to Registration Library

Slicer Registration Library Case #6: Breast MRI Cancer Treatment


this is the main fixed reference image. All images are ev. aligned into this space
RegArrow NonRigid.png
this is the moving image.
exam 1: baseline Pre Rx exam 1: follow-up Post Rx


This is an example of serial/logitudinal imaging for change assessment. We have a pre- and post-treatment MRI. To assess the treatment effect we seek to register the postRx scan to the baseline. Because of the strong changes in image content and contrast between the two exams, some additional work is necessary to obtain a good registration: 1) we will register only the left breast containing the tumor, i.e. we will crop the image to a much smaller FOV, while also obtaining a more isotropic resolution; 2) we have strong intensity bias in the images (originating from MR receiver coil sensitivity) that needs to be corrected first, otherwise it will negatively bias the registration; and 3) we will generate masks that include the breast but not the tumor, to ensure the nonrigid portion of the registration will not try to compensate for the tumor change.

Modules used

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.


MRI, breast, intra-subject, breast cancer, treatment effect

Video Screencasts

  1. Movie/screencast showing importing and cropping the image volume
  2. Movie/screencast showing the bias (intensity inhomogeneity) correction step
  3. Movie/screencast showing building the masks for registration
  4. Movie/screencast showing affine + BSpline registration


  1. Crop:Extract left breast image of PreRx and PostRx scan
    1. Open the Crop Volume module (under Converters menu)
    2. Input volume: "PreRx"
    3. Input ROI: "Create & rename new annotation ROI", rename to "PreRx_CropROI" or similar
    4. Isotropic output voxel: yes (checkbox)
    5. Input spacing scaling constant:increase to 3 (the original image has a 0.4 x 0.4 x 5 mm voxel size. Isotropic resolution is highly desirable for registration, but the full resolution would yield a 784^3 image volume that would use too much memory for purposes of this demonstration).
    6. Interpolator: "Linear" (radio button)
    7. you should see a set of colored dots/beads in at least one of the slice views. If not scroll a bit to locate them. The beads represent sides and corners of an ROI box. Drag the corners of the box to completely enclose the left breast (right side of image).
    8. click on the Crop! button
    9. go to the Data module
      1. several new nodes were created: look for the "PreRx_subvolume-scale_3" entry, double click and rename to "PreRx_Left" or "PreRx_crop" or similar.
      2. Turn off the visibility of the ROI box just created: click on the eye icon next to the "PreRx_CropROI" node. (we do this to more easily distinguish the next ROI box created for the second volume below)
    10. repeat the same for the "PostRx" image. Since this image is smaller, choose 2 for the Input spacing scaling constant (vs. 3 above).
    11. Save intermediate results. You should have now 2 new volumes at ~ 1.3 mm isotropic resolution.
  2. Bias Correction: Correct for local intensity inhomogeneities
    1. Open the N4ITK MRI Bias Correction module (under Filters menu)
      1. Input Image: "PreRx_left"
      2. Mask Image: none
      3. Output Volume: create & rename new: "PreRx_left_n4"
      4. Number of iterations: reduce to 100,80,60
      5. Shrink factor: reduce to 3
      6. leave all parameters at defaults
      7. Apply.
    2. repeat for the Post_Rx image. Same settings
    3. save intermediate results
  3. Build Masks:
    1. open the Foregroud masking (BRAINS) module (under Segmentation:Specialized)
      1. Input Image Volume: PreRx_left_n4
      2. Output Mask: create & rename new: "PreRx_mask"
      3. Apply
    2. (alternatively you can also use the threshold tool in the Editor to obtain a quick mask)
    3. we need to exclude the tumor from the mask, to prevent the registration from trying to match that region. Much should already be excluded from the above segmentation, but we also need to dilate the mask a bit to extend beyon the edge and include some of the background, which will fill in the tumor space again. Hence we make some manual edits to delete the tumor area from the mask:
    4. Go to the Editor module
      1. Click Apply on the popup do confirm colormap selection (choice does not matter)
      2. Master Volume: PreRx_left
      3. Merge Volume: PreRx_mask (generated above)
    5. click on the Brush Icon ("PaintEffect")
    6. Label: 0 (black)
    7. increase radius to ~ 10-15 mm.
    8. navigate to an axial slice showing the tumor
    9. click & drag left mouse to apply, then use arrow keys to move to the next slice.
    10. repeat for the Post_Rx_mask
    11. save intermediate results
  4. Affine pre-Registration:
    1. open the * General Registration (BRAINS) module
      1. Fixed Image Volume: PreRx_left_n4
      2. Moving Image Volume: PostRx_left_n4
      3. Output Settings:
      4. Slicer BSpline Transform": none
      5. Slicer Linear Transform: create & rename new transform, rename to "Xf1_Affine"
      6. Output Image Volume: none
      7. Initialize Transform Mode: select Use Moments Align.
      8. Registration Phases: check boxes for Rigid , Rigid+Scale and Affine
      9. Main Parameters : Number Of Samples: 200,000
      10. Leave all other settings at default
      11. Apply
      12. this should generate a first alignment. note that we need not generate an output volume, because this is an intermediate result, and the BRAINS registration tool will automatically place the moving volume inside the result transform.
      13. Place "PreRx_left_n4" in the background and "PostRx_left_n4" in the foreground. Use fade slider to evaluate the result.
  5. Nonrigid final Registration:
    1. open the General Registration (BRAINS) module
      1. Fixed Image Volume: PreRx_left_n4
      2. Moving Image Volume: PostRx_left_n4
      3. Initialization transform: "Xf1_Affine
      4. Initialize Transform Mode: de-select Use Moments Align, select None
      5. Output Settings:
      6. Slicer BSpline Transform": create & rename new transform, rename to "Xf2_BSpline"
      7. Slicer Linear Transform: none
      8. Output Image Volume: create & rename new transform, rename to "PostRx_left_Xf2" (this time we must request an output volume, since the nonrigid deformations cannot be visualized on the "fly.
      9. Registration Phases: uncheck rigid & affine boxes, check box for BSpline only
      10. Main Parameters:
      11. Number Of Samples: 300,000
      12. B-Spline Grid Size: 7,7,5
      13. Mask Option: select ROI button
        1. ROI Masking input fixed: select "PreRx_mask" generated in Phase 3 above
      14. ROI Masking input moving: select "PostRx_mask" generated in Phase 3 above
      15. Leave all other settings at default
      16. click: Apply

Registration Results

RegLib C04 unregistered1.gif baseline PD & contrast T1Gd before registration (click to enlarge)
RegLib C04 unregistered2.gif baseline and follow-up PD before registration (click to enlarge)
RegLib C04 registered1.gif baseline PD & contrast T1Gd after registration (click to enlarge)
RegLib C04 registered2.gif baseline and follow-up PD after registration (click to enlarge)
RegLib C04 registered3.gif baseline and follow-up T1Gd contrast after registration(click to enlarge)