Difference between revisions of "Projects:RegistrationLibrary:RegLib C28"

From NAMIC Wiki
Jump to: navigation, search
(Created page with 'Back to ARRA main page <br> Back to Registration main page <br> [[Projects:RegistrationDocumentation:UseCaseInv…')
 
m (Text replacement - "http://www.slicer.org/slicerWiki/index.php/" to "https://www.slicer.org/wiki/")
 
(8 intermediate revisions by one other user not shown)
Line 4: Line 4:
  
  
= <small>v3.6.1</small> [[Image:Slicer3-6Announcement-v1.png‎|150px]]  Slicer Registration Library Case 28: <br> In Vitro cerebral body CT (intra subject)  =
+
= <small>v3.6.1</small> [[Image:Slicer3-6Announcement-v1.png‎|150px]]  Slicer Registration Library Case 28:  In Vitro cerebral body CT (intra subject)  =
 
=== Input ===
 
=== Input ===
 
{| style="color:#bbbbbb; " cellpadding="10" cellspacing="0" border="0"
 
{| style="color:#bbbbbb; " cellpadding="10" cellspacing="0" border="0"
Line 17: Line 17:
  
 
=== Modules ===
 
=== Modules ===
*'''Slicer 3.6.1 recommended modules:  [http://www.slicer.org/slicerWiki/index.php/Modules:BRAINSFit BRAINSFit]
+
*'''Slicer 3.6.1 recommended modules:  [https://www.slicer.org/wiki/Modules:BRAINSFit BRAINSFit]
  
 
===Objective / Background ===
 
===Objective / Background ===
Line 36: Line 36:
 
#Scroll through sagittal view and note the inverted axis orientation
 
#Scroll through sagittal view and note the inverted axis orientation
 
#Flip axis
 
#Flip axis
##Go to the [http://www.slicer.org/slicerWiki/index.php/Modules:Data-Documentation-3.6 Data module]
+
##Go to the [https://www.slicer.org/wiki/Modules:Data-Documentation-3.6 Data module]
 
###right-click on the ''Scene'' node and select "Insert Transform Node" from the menu
 
###right-click on the ''Scene'' node and select "Insert Transform Node" from the menu
 
###select the newly created transform node, then in the ''MRML Node Inspector'' below, rename it to "Xf0_AxialFlip"
 
###select the newly created transform node, then in the ''MRML Node Inspector'' below, rename it to "Xf0_AxialFlip"
###select the "L4_intact" node and drag it inside/on top of the Xf0_AxialFlip" node
+
###select the "L4_intact" node and drag it inside/on top of the "Xf0_AxialFlip" node
##Go to the [http://www.slicer.org/slicerWiki/index.php/Modules:Transforms-Documentation-3.6 Transforms module]
+
###ensure you have L4_intact selected for background and the fade slider set on background (B) so you can see changes happening.
 +
##Go to the [https://www.slicer.org/wiki/Modules:Transforms-Documentation-3.6 Transforms module]
 
###Select "Xf0_AxialFlip" from the ''Transform Node'' menu
 
###Select "Xf0_AxialFlip" from the ''Transform Node'' menu
 
###In the 4x4 matrix displayed, click in the field containing the "1.0000" in row 3 and column 3.  
 
###In the 4x4 matrix displayed, click in the field containing the "1.0000" in row 3 and column 3.  
Line 46: Line 47:
 
###change the number to -1. Hit Enter.
 
###change the number to -1. Hit Enter.
 
###you should see the image flip in the sagittal and coronal slice views
 
###you should see the image flip in the sagittal and coronal slice views
 +
#Crop Volumes
 +
##Go to the [https://www.slicer.org/wiki/Modules:CropVolume-Documentation-3.6 Crop Volume module]
 +
##select "L4_intact" as ''Input Volume''; select "Create new MRML ROI Node" in the ''ROI'' menu, select "create new volume" in the ''Output Volume" menu
 +
##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)
 +
##when finished, click on ''Do ROI resample''
 +
##go to back to the [https://www.slicer.org/wiki/Modules:Data-Documentation-3.6 Data module] and change the name of the new "SubVolume" to "L4_intact_crop"
 +
##delete the ROI nodes (right click: delete node)
 +
##repeat the same for the "L4_defect" image
 +
#Manual Initial Alignment
 +
##in the [https://www.slicer.org/wiki/Modules:Data-Documentation-3.6 Data module] again select "Insert new Transform node" as done above
 +
##select the newly created transform node, then in the ''MRML Node Inspector'' below, rename it to "Xf1_ManualInit" or similar
 +
##select the "L4_defect_crop" node and drag it inside/on top of the "Xf1_ManualInit" node
 +
##ensure you have L4_intact_crop selected for background and  L4_defect_crop  for background
 +
##set the fade slider to halfway between background (B)  and foreground (F) so you can see changes happening.
 +
##Go to the [https://www.slicer.org/wiki/Modules:Transforms-Documentation-3.6 Transforms module]
 +
##Select "Xf1_ManualInit"  from the ''Transform Node'' menu
 +
##adjust the sliders to roughly align the two images
 +
##for more details on manual registration, see the "Manual Registration" tutorial here: [https://www.slicer.org/wiki/Slicer_3.6:Training]
 +
#Run final registration
 +
##Go to the [https://www.slicer.org/wiki/Modules:BRAINSFit BRAINSfit module]
 +
###''''Presets''': select the "BRAINSFit_Affine" from the parameter set menu. See [[Projects:RegistrationDocumentation:ParameterPresetsTutorial|here for a guide]] on how to Load/Save Registration Parameter Presets]]
 +
###otherwise select the following parameters:
 +
###''Fixed Image Volume'': :L4_intact_crop ; ''Moving'': L4_defect_crop
 +
###''Initialize with previously generated transform'':  select "Xf1_ManualInit"  generated above
 +
###check boxes for ''Include Rigd registr. phase'' , ''Include ScaleVersor3D'', ''include Affine''
 +
###''Slicer Linear Transform'': select "create new transform", rename to "Xf1_Affine" or similar
 +
###leave rest at defaults. Click ''Apply''
 +
###registration should take ~ 10 secs.
 +
###use fade slider to verify alignment; compare with result snapshots shown below.
  
##Input Surface: S02_femur
 
##Target Surface: S01_femur
 
##Output Transform: create new transform, rename to "Xf1_surf_S2-S1"
 
##Select/check ''Mean distance mode: RMS''
 
##Leave defaults: iterations: 50, landmarks: 200
 
##Select/check ''Start by matching centroids''
 
##Select/check ''Check mean distance''
 
##Click: ''Apply''; process duration: ~ seconds
 
#Go to the ''Data'' module
 
##find the newly created transform node "Xf1_surf_S2-S1" and drag the "S02_femur" model node (and all other S02 entities) inside the transform.
 
##for the 3D viewer turn on the visibility of the surface models (should be on already)
 
##the surfaces for the two femur models should now align in basic position and orientation
 
#for additional details see tutorial downloads below
 
 
 
=== Registration Results===
 
=== Registration Results===
[[Image:RegLib C05 AGif reg.gif|200px|femur surface models before/after registration]]
+
[[Image:RegLib_C28_unreg.gif|400px|original unregistered]] unregistered <br>
 +
[[Image:RegLib_C28_reg.gif|400px|registered (after flipping axis and cropping volume)]] registered <br>
  
 
===Download ===
 
===Download ===
 
*Data
 
*Data
**[[Media:RegLib_C05_Data.zip‎|'''RegLib_C05_Data''': Example Dataset <small> (MRI images, surface models,Presets, Solution, zip file 9.5 MB) </small>]]
+
**[[Media:RegLib_C28_Data.zip‎|'''RegLib_C28_Data''': Example Dataset <small> (CT images, Presets, Solution, zip file 194 MB) </small>]]
 
*Presets
 
*Presets
**[[Media:RegLib_C05_Presets.mrml‎|'''RegLib_C05_Presets.mrml''': registration parameter presets<small> (.mrml text file 12 kB) </small>]]
+
**[[Media:RegLib_C28_Presets.mrml‎|'''RegLib_C28_Presets.mrml''': registration parameter presets<small> (.mrml text file 12 kB) </small>]]
*[[Projects:RegistrationDocumentation:ParameterPresetsTutorial|Link to User Guide: How to Load/Save Registration Parameter Presets]]
 
*Documentation
 
**[[Media:RegLib_C05_Tutorial.ppt|'''RegLib_C05_Tutorial''': Powerpoint file <small> (.ppt file 1 MB) </small>]]
 
**[[Media:RegLib_C05_Tutorial.ppt|'''RegLib_C05_Tutorial''': PDF file <small> (.pdf file 0.5 MB) </small>]]
 
  
 
=== Acknowledgments ===
 
=== Acknowledgments ===
 +
Thanks to Dr.Ron Alkalay and the Center for Advanced Orthopaedic Studies at Beth Israel Deaconess Medical Center (BIDMC) for providing the example data. Details about their research [http://obl.bidmc.harvard.edu/research.html can be found here].

Latest revision as of 18:07, 10 July 2017

Home < Projects:RegistrationLibrary:RegLib C28

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


v3.6.1 Slicer3-6Announcement-v1.png Slicer Registration Library Case 28: In Vitro cerebral body CT (intra subject)

Input

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

Modules

Objective / Background

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

Keywords

CT, spine cancer, vertebrae

Issues Challenges

  • 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
  • the two vertebrae have strong rotational misalignment; manual initial alignment is recommended

Input Data

  • fixed: CT , 0.342 x 0.342 x 0.20 mm voxel size, axial; 512 x 512 x 193 unsigned short image
  • moving: CT , 0.342 x 0.342 x 0.20 mm voxel size, axial; 512 x 512 x 250 unsigned short image

Procedure / Pipeline

  1. Open case scene file or import image data: RegLib_C28_SlicerScene1.mrml
  2. Place L4_intact and L4_defect 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 "L4_intact" node and drag it inside/on top of the "Xf0_AxialFlip" node
      4. ensure you have L4_intact 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 "L4_intact" 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 "L4_intact_crop"
    6. delete the ROI nodes (right click: delete node)
    7. repeat the same for the "L4_defect" image
  6. Manual Initial Alignment
    1. in the Data module again select "Insert new Transform node" as done above
    2. select the newly created transform node, then in the MRML Node Inspector below, rename it to "Xf1_ManualInit" or similar
    3. select the "L4_defect_crop" node and drag it inside/on top of the "Xf1_ManualInit" node
    4. ensure you have L4_intact_crop selected for background and L4_defect_crop for background
    5. set the fade slider to halfway between background (B) and foreground (F) so you can see changes happening.
    6. Go to the Transforms module
    7. Select "Xf1_ManualInit" from the Transform Node menu
    8. adjust the sliders to roughly align the two images
    9. for more details on manual registration, see the "Manual Registration" tutorial here: [1]
  7. Run final 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: :L4_intact_crop ; Moving: L4_defect_crop
      4. Initialize with previously generated transform: select "Xf1_ManualInit" generated above
      5. check boxes for Include Rigd registr. phase , Include ScaleVersor3D, include Affine
      6. Slicer Linear Transform: select "create new transform", rename to "Xf1_Affine" or similar
      7. leave rest at defaults. Click Apply
      8. registration should take ~ 10 secs.
      9. 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

Download

Acknowledgments

Thanks to Dr.Ron Alkalay and the Center for Advanced Orthopaedic Studies at Beth Israel Deaconess Medical Center (BIDMC) for providing the example data. Details about their research can be found here.