2012 Winter Project Week:CTLiverRegistration

From NAMIC Wiki
Jump to: navigation, search
Home < 2012 Winter Project Week:CTLiverRegistration

Key Investigators

  • AZE Research and Development, Karl Diedrich
  • Brigham and Women's Hospital, Nobuhiko Hata
  • Brigham and Women's Hospital, Atsushi Yamada


Register CT images of livers taken at different time points for the purpose of detecting changes in liver tumors.

Approach, Plan

  1. Acquire unregistered 4D liver data sets.
  2. Segment bone using upper and lower Hounsfield unit threshold.
  3. Create a distance map from the threshold segmented bones. Look into using vtkITKDistanceTransform for distance map. The bones are very sparse making poor registration subjects.
  4. Conduct a fixed registration of the fixed and moving bones near the liver using BRAINSfit in slicer.
  5. Apply the resulting transform to the full moving liver data set.


  1. Acquired additional unregistered 4D liver time series data sets for development and testing from http://midas.kitware.com/community/view/47.
    • Rotationally different data. LiverFullA B unreg.png
    • Different point in respiration to test non-rigid registration. Stanford2LiverPatient0 0 90saggital.png
  2. Used Slicer to segment the bones using a lower threshold of 252. Segmented the bones with a label of 255 to allow blurring. See upper left images below.
  3. ITK contains a filter class, itkSignedDanielssonDistanceMapImageFilter, http://www.itk.org/Doxygen/html/classitk_1_1SignedDanielssonDistanceMapImageFilter.html for creating a distance map. Look into using vtkITKDistanceTransform for distance map. Need to develop a Slicer CLI module for generating distance maps of the fixed and moving bone images.
  4. Registered two time points of the same patient as fixed and moving images in BRAINSfit using Initialize Transform Mode: useGeometryAlign Registration Phase: Rigid (6 DOF) AND Rigid+Scale (7DOF) AND Rigid+Scale+Skew (10 DOF). See lower left images below.
  5. Applied the resulting transform to the full moving liver data set in the Slicer Data module Nodes panel by dragging the moving liver data set onto the output transform name. See lower right images below.

LiverSceneViewAxial.png LiverSceneViewSaggital.png LiverFullTransformedFixed.png Liver1 2coronal.png

Delivery mechanism

Develop a Slicer CLI module using the ITK distance map itkSignedDanielssonDistanceMapImageFilter, http://www.itk.org/Doxygen/html/classitk_1_1SignedDanielssonDistanceMapImageFilter.html. Then chain the steps together in Python.


Hans Johnson, Greg Harris, Kent Williams. BRAINSFit: Mutual Information Rigid Registrations of Whole-Brain 3D Images, Using the Insight Toolkit. The Insight Journal. 2009.