Projects:RegistrationLibrary:RegLib C44

Back to ARRA main page
Back to Registration main page
Back to Registration Use-case Inventory

v3.6.3 Slicer Registration Library Case 44: Visible Human Pelvis CT

Input

baseline image		follow-up

Modules

• Slicer 3.6.3 recommended modules: BRAINSFit

Objective / Background

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.

Download

• Data
  • RegLib_C44_Data: Example Dataset (CT images, Presets, Solution, zip file 171 MB)
• Presets
  • RegLib_C44_Presets.mrml: registration parameter presets (.mrml text file 16 kB)

Keywords

CT, pelvis, visible human, inter-subject

Input Data

• fixed: CT , 0.6 x 0.6 x 0.40 mm voxel size, axial; 292 x 292 x 91 unsigned short image
• moving: CT , 0.6 x 0.6 x 0.40 mm voxel size, axial; 512 x 512 x 107 unsigned short image

Overall Strategy

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

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. 'Presets: for either of the registrations below, you can select the "BRAINSFit_...." presets from the parameter set menu instead of setting the parameters manually.
3. Affine Registration
   1. Go to the BRAINSfit module
   2. select the following parameters:
      1. Fixed Image Volume: vhm Moving: vhf
      2. check boxes for Include Rigd registr. phase , Include ScaleVersor3D, include Affine
      3. Slicer Linear Transform: select "create new transform", rename to "Xf1_Affine" or similar
      4. leave rest at defaults. Click Apply
      5. registration should take ~ 10 secs.
      6. use fade slider to verify alignment; compare with result snapshots shown below. Alignment will not be perfect but should be better than before.
4. BSpline Registration
   1. Go to the BRAINSfit module
   2. select the following parameters:
      1. Fixed Image Volume: vhm Moving: vhf
   3. Registration phases: from Initialize with previously generated transform', select "Xf2_..." node created before.
   4. Registration phases: check boxes for Include BSpline registration phase
   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.