Documentation:Nightly:Registration:RegistrationLibrary:RegLib C44

Back to Registration Library

Slicer Registration Library Case 44: Visible Human Pelvis CT

Input

baseline image		follow-up

Modules used

• General Registration (BRAINS)

Description

This dataset contains CT of the visible human male and female pelvis. This serves as a test example for exploring non-rigid registration for inter-subject comparison from CT.

1. Open case scene file or import image data: RegLib_C44_SlicerScene.mrml
2. Overall strategy will be
   1. use vhm as fixed, vhf as moving volume
   2. because we're interested in the pelvic bone structure differences, we create a mask that will focus the registration on the skeletal region only.
   3. perform affine alignment
   4. using above affine alignment as starting point, perform low-level BSpline alignment
   5. To study the deformation only, extract the BSpline deformation from the transform, i.e. remove the affine portion. This can be done in a text editor, or by first resampling the vhf volume with the Affine before computing the BSpline.
   6. Visualize the deformation by applying the pure BSpline to a auxiliary grid image
3. Change colormap of deformed grid (e.g. hot), window and level to view the gridlines of interest, and overlay on the grayscale images

Download

RegLib_C44.mrb
(input data only, Slicer mrb file. 74 MB)
RegLib_C44_full.mrb
(input data + results, Slicer mrb file. 82 MB).

Video Screencasts

1. Movie/screencast showing generating a registration mask
2. Movie/screencast showing affine and nonrigid BSpline registration
3. Movie/screencast showing visualization of the deformation via the Transform Visualizer module

Keywords

CT, pelvis, visible human, inter-subject

Procedure / Pipeline

1. Mask generation
   1. Go to the Editor module
   2. "Master Volume": select vhm
   3. A new labelmap "vhm-label" will be created
   4. Select "vhm" to be visible in the slice viewer
   5. Select the Threshold tool from the editor toolbar
   6. Adjust the lower threshold (slider bar) until most of the bone is highlighted, e.g. somewhere around an intensity value of 80. Leave the upper threshold unchanged at the max.
   7. Click Apply
   8. Morphologically clean the segmentation:
   9. Run a Median Filter with a 3x3x3 neighborhood to remove speckle noise
   10. Return to the Editor module and use the Change Island tool to remove the segmentation of the arms.
   11. Apply 2-3 rounds of morphological dilation to expand the region to include surrounding area.
   12. Rename the labelmap to "vhm_mask" or similar. You will find a "vhm_mask_dil" in the example dataset for comparison
   13. Repeat the above segmentation procedure for "vhf"
2. Registration:
   1. Go to the General Registration (BRAINS)
      1. Fixed Image Volume: vhm
      2. Moving Volume: vhf
      3. check boxes for Include Rigd registr. phase , Include ScaleVersor3D, include Affine
      4. Slicer Linear Transform: select "create new transform", rename to "Xf1_Affine" or similar
      5. leave rest at defaults. Click Apply
      6. registration should take ~ 10 secs.
      7. use fade slider to verify alignment; compare with result snapshots shown below. Alignment will not be perfect but should be better than before.
      8. note: you can also change the colormaps for the fixed and moving volumes to better judge the alignment: go to the Volumes module and in the Display tab, select "green" and "magenta" as the respective colormaps for the two volumes (vhf, vhm)
3. BSpline Registration
   1. Go to the General Registration (BRAINS)
      1. Fixed Image Volume: vhm Moving: vhf
   2. Registration phases: from Initialize with previously generated transform', select "Xf1_Affine" node created before.
   3. Registration phases: uncheck boxes for rigid, scale and affine and check box for BSpline
   4. Output: Slicer Linear transform: set to None
   5. Output: Slicer BSpline transform: create new, rename to "Xf2_BSpline_msk" or similar
   6. Output Image Volume: create new, rename to "vhf_Xf2"; Pixel Type: "short"
   7. Registration Parameters: increase Number Of Samples to 200,000; Number of Grid Subdivisions: 7,7,7
   8. Control Of Mask Processing Tab: check ROI box, for Input Fixed Mask and Input Moving Mask select the two dilated labelmaps from above
   9. Leave all other settings at default
   10. click apply

Registration Results

unregistered
registered (affine)]
registered (nonrigid w/o masking)
registered (nonrigid+masking)
registered deformation only of vhf
deformation visualized by grid image overlay

Acknowledgments

Original CT from the Visible Human Project shared by the University of Iowa.