Projects:RegistrationLibrary:RegLib C39

From NAMIC Wiki
Jump to: navigation, search
Home < Projects:RegistrationLibrary:RegLib C39

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

v3.6.3 Slicer3-6Announcement-v1.png Slicer Registration Library Case 39: In Vitro cerebral body CT (thoracic spine)


this is the main fixed reference image. All images are ev. aligned into this space lleft this is the moving image
baseline image follow-up


Objective / Background

This dataset contains in vitro CT of two human vertibral body segments (T10) before and after regional damage. Goal is to align the two images to study in detail the local changes that occurred.


CT, spine cancer, vertebrae

Issues Challenges

  • two images have very different spatial origins (as defined in original DICOM header): recenter first
  • the baseline container contains the vertebra upside down (as apparent from orientation of the pedicle/roots of the vertibral arches coming out of the body, visible in sagittal or coronal view). Hence this image needs to have its z-axis (IS) reversed before registration
  • the containers walls and lid contain dominant image content that would prevent a successful automated registration; cropping is required

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

Procedure / Pipeline

  1. Open case scene file or import image data: RegLib_C28_SlicerScene1.mrml
  2. Place T10intact and T10defect in back- and foreground respectively.
  3. Scroll through sagittal view and note the inverted axis orientation
  4. Flip axis
    1. Go to the Data module
      1. right-click on the Scene node and select "Insert Transform Node" from the menu
      2. select the newly created transform node, then in the MRML Node Inspector below, rename it to "Xf0_AxialFlip"
      3. select the "T10intact" node and drag it inside/on top of the "Xf0_AxialFlip" node
      4. ensure you have T10intact selected for background and the fade slider set on background (B) so you can see changes happening.
    2. Go to the Transforms module
      1. Select "Xf0_AxialFlip" from the Transform Node menu
      2. In the 4x4 matrix displayed, click in the field containing the "1.0000" in row 3 and column 3.
      3. hit Return key to activate editing
      4. change the number to -1. Hit Enter.
      5. you should see the image flip in the sagittal and coronal slice views
  5. Crop Volumes
    1. Go to the Crop Volume module
    2. select "T10intact" as Input Volume; select "Create new MRML ROI Node" in the ROI menu, select "create new volume" in the Output Volume" menu
    3. you should see a small blue rectangle in all 3 views. Click in the views to expand the region to contain only the L4 vertebra (body and processes)
    4. when finished, click on Do ROI resample
    5. go to back to the Data module and change the name of the new "SubVolume" to "T10intact_crop"
    6. delete the ROI nodes (right click: delete node)
    7. repeat the same for the "T10defect" image
  6. Registration
    1. Go to the BRAINSfit module
      1. 'Presets: select the "BRAINSFit_Affine" from the parameter set menu. See here for a guide on how to Load/Save Registration Parameter Presets]]
      2. otherwise select the following parameters:
      3. Fixed Image Volume: :T10intact_crop ; Moving: T10defect_crop
      4. check boxes for "Moments align", Include Rigd registr. phase , Include ScaleVersor3D, include Affine
      5. Slicer Linear Transform: select "create new transform", rename to "Xf1_Affine" or similar
      6. leave rest at defaults. Click Apply
      7. registration should take ~ 10 secs.
      8. use fade slider to verify alignment; compare with result snapshots shown below.

Registration Results

original unregistered unregistered
registered (after flipping axis and cropping volume) registered