Documentation:Nightly:Registration:RegistrationLibrary:RegLib C39 - Revision history

Meier: /* Download (from NAMIC MIDAS) */

2013-10-03T15:41:11Z

Download (from NAMIC MIDAS)

Meier: /* Download (from NAMIC MIDAS) */

2013-10-03T15:41:03Z

Download (from NAMIC MIDAS)

Meier at 15:39, 3 October 2013

2013-10-03T15:39:42Z

Meier at 15:37, 3 October 2013

2013-10-03T15:37:34Z

Meier: /* Procedure */

2013-10-03T15:35:27Z

Procedure

Meier: Created page with 'Back to Registration Library = Slicer Registration Library Case #39: In Vitro cerebral body CT (thoracic spine) = ==…'

2013-10-03T14:05:10Z

Created page with 'Back to Registration Library = Slicer Registration Library Case #39: In Vitro cerebral body CT (thoracic spine) = ==…'

New page

[[Documentation/Nightly/Registration/RegistrationLibrary|Back to Registration Library]]
= Slicer Registration Library Case #39: In Vitro cerebral body CT (thoracic spine) =
== Input ==
{| style="color:#bbbbbb; " cellpadding="10" cellspacing="0" border="0"
|[[Image:RegLib_C39_Thumb1.png|150px|lleft|this is the fixed reference image. All images are aligned into this space]]
|[[Image:RegArrow_Rigid.png|100px|lleft]]
|[[Image:RegLib_C39_Thumb2.png|150px|lleft|this is the moving image. The transform is calculated by matching this to the reference image]]
|-
|fixed image/target
|
|moving image
|}

== Modules used ==
*[[Documentation/Nightly/Modules/BRAINSFit| ''General Registration (BRAINS)'']]

== 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. ''
*[http://slicer.kitware.com/midas3/download/?items=xx '''RegLib_C39.mrb'''] (input data only, Slicer mrb file. 40 MB) 
*[http://slicer.kitware.com/midas3/download/?items=xx '''RegLib_C39_full.mrb'''] (input data + results, Slicer mrb file. 100 MB). 

== Description ==
We have 2 in vitro CT scans of a human thoracic spine, before and after application of a controlled structural damage. We apply a rigid registration to bring the two images into alignment. Because the container tube and fixtures add image content that is "distracting" to the registration, we crop both images to contain just the vertebra of interest. Also the spines were not imaged in the same orientation, and as consequence we have an inverted z-axis in one of the images we have to correct first. Finally the two images have very different spatial origins (as defined in original DICOM header), so we need to recenter first before running the registration.

=== Keywords ===
CT, spine cancer, vertebrae

== Video Screencasts ==
#[[Media:RegLib_C39_1_WindowLevelCenter.mov|Movie/screencast showing window & level and re-centering steps]]
#[[Media:RegLib_C39_2_OrientationFix.mov|Movie/screencast showing steps for fixing the orientation issue (inverting z-axis)]]
#[[Media:RegLib_C39_3_CropVolumes.mov|Movie/screencast showing cropping of both volumes to the region of interest]]
#[[Media:RegLib_C39_4_Register.mov|Movie/screencast showing rigid registration and applying to image '''without''' resampling]]

== Procedure ==
#'''Window & Level , Recenter''': open the [[Documentation/Nightly/Modules/Volumes| ''Volumes'']] module.
##in the ''Display'' tab, click on "Auto W/L" to get a preliminary adjustment, then drag left mouse key in the slice views horizontal & vertical to adjust contrast
##once you have acceptable levels, note the window & level entries in the display tab, then change the active volume and apply the same settings to the other image by typing them directly into the window and level fields. When judging alignment with very high precision, it is good practice to ensure window & level settings are identical for both volumes
##for both volumes, click the "Center Volume" button to reset the origin
##place "T10_intact" and "T10_defect" into foreground and background, respectively. Use the fade slider to assess misalignment
#'''Orientation Fix''': open the [[Documentation/Nightly/Modules/Data| ''Data'']] module.
##you will notice that the "T10_intact" spine is upside down (sagittal view) compared to "T10_defect". This is because the tube was inverted during the 2nd acquisition. We will correct this with a transform.
##in the Data module, right click on the "Scene" node and select "Insert Transform" from the pop-up menu
##rename the transform to "Xf0_OrientationFix" or similar
##drag the "T10_intact" volume node into/onto the transform node.
##adjust sagittal view so that "T10_intact" is visible
##switch to the [[Documentation/Nightly/Modules/Transforms| ''Transforms'']] module.
##the newly created "Xf0_OrientationFix" transform should already be selected.
##Note the 3x3 matrix showing mostly zeros and three 1.0 along the diagonal. Double click on the last 1.0 in the bottom right corner. It represents the direction of the z-axis which we want to invert. Change the sign, i.e. replace "1.0" with "-1.0". Then hit ENTER.
##you should see the view update and the "T10_intact" image flip.
##return to the [[Documentation/Nightly/Modules/Data| ''Data'']] module.
##right-click on the "T10_intact" node inside the transform and select "Harden Transform". This will apply the new orientation to the image header. The node will move back out t the main level.
##double click on "T10_intact" and rename it to "T10_intact_flip" or similar. Note that your original image has been modified. To avoid accidentally overwriting the original on disk with the new version it is important to immediately rename the node upon applying the transform.
##save new "T10_intact_flip" image under a new name.
#'''Crop''': open the [[Documentation/Nightly/Modules/Crop| ''Crop'']] module.
#
#'''Rigid Registration''': open the [[Documentation/Nightly/Modules/BRAINSFit| ''General Registration (BRAINS)'']] module.
##''Fixed Image Volume'': "T10_intact"
##''Moving Volume'': "T10_defect"
##Registration phases: check box for ''Rigid" only
##''Slicer Linear Transform'': select "create new transform", rename to "Xf1_Rigid" or similar
##increase sample points to 200,000
##leave rest at defaults. Click ''Apply''
##registration should take ~ 10 secs.
##use fade slider to verify alignment; compare with result snapshots shown below.
##note: you can also change the colormaps for the fixed and moving volumes to better judge the alignment: go to the [[Documentation/Nightly/Modules/Volumes| ''Volumes'']] module and in the ''Display'' tab, select "green" and "magenta" as the respective colormaps for the two volumes.

== Registration Results==
[[Image:RegLib_C39_unregistered.gif|300px]] before registration (click to enlarge) 
[[Image:RegLib_C39_registered.gif|300px]] after rigid alignment. (click to enlarge) 

==Acknowledgments ==
Registration problem and dataset kindly provided by Dr. Ron Alkalay from the Beth Israel Deaconess Medical Center.

@@ Line 18: / Line 18: @@
 <small>''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. ''</small>
 *[http://slicer.kitware.com/midas3/download/?items=119824 '''RegLib_C39.mrb''']  <small>(input data only, Slicer mrb file. 25 MB) </small><br>
-*[http://slicer.kitware.com/midas3/download/?items=119823 '''RegLib_C39_full.mrb''']<br><small>(input data + results, Slicer mrb file. 45 MB). </small>
+*[http://slicer.kitware.com/midas3/download/?items=119823 '''RegLib_C39_full.mrb'''] <small>(input data + results, Slicer mrb file. 45 MB). </small>
 == Description ==

@@ Line 17: / Line 17: @@
 == Download (from NAMIC MIDAS) ==
 <small>''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. ''</small>
-*[http://slicer.kitware.com/midas3/download/?items=xx '''RegLib_C39.mrb'''] <br><small>(input data only, Slicer mrb file. 40 MB) </small><br>
+*[http://slicer.kitware.com/midas3/download/?items=119824 '''RegLib_C39.mrb''']  <small>(input data only, Slicer mrb file. 25 MB) </small><br>
-*[http://slicer.kitware.com/midas3/download/?items=xx '''RegLib_C39_full.mrb''']<br><small>(input data + results, Slicer mrb file. 100 MB). </small>
+*[http://slicer.kitware.com/midas3/download/?items=119823 '''RegLib_C39_full.mrb''']<br><small>(input data + results, Slicer mrb file. 45 MB). </small>
 == Description ==

@@ Line 33: / Line 33: @@
 == Procedure ==
-#'''Window & Level , Recenter''': open the  [[Documentation/Nightly/Modules/Volumes| ''Volumes'']] module. ([[Media:RegLib_C39_1_WindowLevelCenter.mov|Movie/screencast]])
+#'''Window & Level , Recenter''': open the  [[Documentation/Nightly/Modules/Volumes| ''Volumes'']] module. '''( [[Media:RegLib_C39_1_WindowLevelCenter.mov|Movie/screencast]] )'''
 ##in the ''Display'' tab, click on "Auto W/L" to get a preliminary adjustment, then drag left mouse key in the slice views horizontal & vertical to adjust contrast
 ##once you have acceptable levels, note the window & level entries in the display tab, then change the active volume and apply the same settings to the other image by typing them directly into the window and level fields. When judging alignment with very high precision, it is good practice to ensure window & level settings are identical for both volumes
 ##for both volumes, click the "Center Volume" button to reset the origin
 ##place "T10_intact" and "T10_defect" into foreground and background, respectively. Use the fade slider to assess misalignment
-#'''Orientation Fix''': open the  [[Documentation/Nightly/Modules/Data| ''Data'']] module.  ([[Media:RegLib_C39_2_OrientationFix.mov|Movie/screencast])
+#'''Orientation Fix''': open the  [[Documentation/Nightly/Modules/Data| ''Data'']] module.  '''( [[Media:RegLib_C39_2_OrientationFix.mov|Movie/screencast]] )'''
 ##you will notice that the "T10_intact" spine is upside down (sagittal view) compared to "T10_defect". This is because the tube was inverted during the 2nd acquisition. We will correct this with a transform.
 ##in the Data module, right click on the "Scene" node and select "Insert Transform" from the pop-up menu
@@ Line 52: / Line 52: @@
 ##double click on "T10_intact" and rename it to "T10_intact_flip" or similar. Note that your original image has been modified. To avoid accidentally overwriting the original on disk with the new version it is important to immediately rename the node upon applying the transform.
 ##save new "T10_intact_flip" image under a new name.
-#'''Crop''':open the  [[Documentation/Nightly/Modules/Crop_Volume| ''Crop Volume'']]  module (under ''Converters'' menu)([[Media:RegLib_C39_3_CropVolumes.mov|Movie/screencast]])
+#'''Crop''':open the  [[Documentation/Nightly/Modules/Crop_Volume| ''Crop Volume'']]  module (under ''Converters'' menu) . '''([[Media:RegLib_C39_3_CropVolumes.mov|Movie/screencast]] )'''
 ##''Input volume'':  "T10_intact_flip"
 ##''Input ROI'': "Create & rename new annotation ROI", rename to "intact_CropROI" or similar
@@ Line 64: / Line 64: @@
 ##repeat the same for the "T10_defect" image.
 ##Save intermediate results.
-#'''Rigid Registration''': open the  [[Documentation/Nightly/Modules/BRAINSFit| ''General Registration (BRAINS)'']] module.
+#'''Rigid Registration''': open the  [[Documentation/Nightly/Modules/BRAINSFit| ''General Registration (BRAINS)'']] module.   '''([[Media:RegLib_C39_4_Register.mov|Movie/screencast]])'''
 ##''Fixed Image Volume'': "T10_intact"
 ##''Moving Volume'': "T10_defect"