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

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…')
 
Line 4: Line 4:
  
 
==Slicer Registration Use Case Exampe #5: Inter-subject Knee MRI ==
 
==Slicer Registration Use Case Exampe #5: Inter-subject Knee MRI ==
 +
 +
[[Image:RegLib C05 KneeMRI1.png|70px|lleft|RegLib 05: Knee MRI subject 1]] [[Image:RegLib C05 KneeMRI2.png|70px|lleft|RegLib 05: Knee MRI subject 2]]
 +
 +
 
{| style="color:#bbbbbb; background-color:#333333;" cellpadding="10" cellspacing="0" border="0"
 
{| style="color:#bbbbbb; background-color:#333333;" cellpadding="10" cellspacing="0" border="0"
|[[Image:RegUC2_thumb_t0.png|200px|lleft|this is the fixed reference image. All images are aligned into this space]]  
+
|[[Image:RegLib C05 KneeMRI1.png|200px|lleft|this is the fixed reference image. All images are aligned into this space]]  
 
|[[Image:Arrow_left_gray.jpg|100px|lleft]]  
 
|[[Image:Arrow_left_gray.jpg|100px|lleft]]  
|[[Image:RegUC2_thumb_t1.png|200px|lleft|this is the moving image. The transform is calculated by matching this to the reference image]]
+
|[[Image:RegLib C05 KneeMRI2.png|200px|lleft|this is the moving image. The transform is calculated by matching this to the reference image]]
 
|align="left"|LEGEND<br>
 
|align="left"|LEGEND<br>
 
[[Image:Button_red_fixed.jpg|20px|lleft]]  this indicates the reference image that is fixed and does not move. All other images are aligned into this space and resolution<br>
 
[[Image:Button_red_fixed.jpg|20px|lleft]]  this indicates the reference image that is fixed and does not move. All other images are aligned into this space and resolution<br>
Line 39: Line 43:
 
===Download ===
 
===Download ===
 
*'''[[Media:RegLib_Case_01_TumorGrowth.zip‎|download entire package  <small> (Data,Presets,Tutorial, Solution, zip file 33.7 MB) </small>]]'''
 
*'''[[Media:RegLib_Case_01_TumorGrowth.zip‎|download entire package  <small> (Data,Presets,Tutorial, Solution, zip file 33.7 MB) </small>]]'''
**[[Media:RegPreset_RegUC-001.txt|download registration parameter presets file  <small> (MRML file, import as scene) </small>]]
 
**[[Media:RegLib_C01_Data_TumorGrowth.zip|download image dataset only  <small>(NRRD, 10.7 MB, filename: RegLib_C01_Data_TumorGrowth.zip) </small>]]
 
**[[Media:RegLib_Case_01_NRRD_TumorGrowth.zip|download image dataset only  <small>(NRRD, 10.7 MB, filename: RegLib_Case_01_NRRD_TumorGrowth.zip) </small> ]]
 
**[[Media:RegLib_C01_DataNIFTI_TumorGrowth.zip|download image dataset in NIFTI format <small>(NIFTI / nii, 10.7 MB, filename: RegLib_C01_DataNIFTI_TumorGrowth.zip) </small> ]]
 
**[[Media:RegXForm_RegUC-001.tfm.txt|result transform file <small>(ITK .tfm file, load into slicer and apply to the target volume)</small>]]
 
**Tutorials (step-by -step walk through):
 
***[[Media:RegLib_C01_VideoTutorial_TumorGrowth.mov|download/play video tutorial <small>(quicktime, 15.9 MB, filename: RegLib_C01_VideoTutorial_TumorGrowth.mov) </small>]]
 
***[[Media:RegLib_C01_PPTTutorial_TumorGrowth.ppt.zip‎|download power point tutorial <small>(zip file, 2.8 MB, filename: RegLib_C01_PPTTutorial_TumorGrowth.ppt.zip) </small>]]
 
***[[Media:RegInstr_RegUC-001.txt‎|download step-by step text instructions <small>(rtf text file) </small>]]
 
*'''[[Media:RegLib_C01_TumorGrowth_MultiresSolution_Dec09.zip‎|Multiresolution testresult package  <small> (Data,Xform, Solution, zip file 16.5 MB) </small>]]'''
 
*'''[http://www.insight-journal.org/midas/item/bitstream/2332 Download package from MIDAS server<small> (Data,Xform) </small>]'''
 
  
 
<!--
 
<!--

Revision as of 21:55, 3 February 2010

Home < Projects:RegistrationLibrary:RegLib C05

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

Slicer Registration Use Case Exampe #5: Inter-subject Knee MRI

RegLib 05: Knee MRI subject 1 RegLib 05: Knee MRI subject 2


this is the fixed reference image. All images are aligned into this space lleft this is the moving image. The transform is calculated by matching this to the reference image LEGEND

lleft this indicates the reference image that is fixed and does not move. All other images are aligned into this space and resolution
lleft this indicates the moving image that determines the registration transform.

lleft T1 SPGR lleft T1 SPGR
0.9375 x 0.9375 x 1.4 mm
256 x 256 x 112
RAS
0.9375 x 0.9375 x 1.2 mm
256 x 256 x 130
RAS

Objective / Background

The final goal is to align a segmentation prior model to aid in cartilage segmentation.

Keywords

MRI, knee, inter-subject, segmentation

Input Data

  • Button red fixed white.jpgreference/fixed : T1 SPGR , 0.9375 x 0.9375 x 1.4 mm voxel size, axial, RAS orientation.
  • Button green moving white.jpg moving: T1 SPGR , 0.9375 x 0.9375 x 1.2 mm voxel size, sagittal, RAS orientation.
  • Content preview: Have a quick look before downloading: Does your data look like this? Media:RegUC2_lightbox.png
  • download dataset to load into slicer (~17 MB zip archive)

Registration Results

after orientation
before/after registration, indiv. images in red/cyan, match in gray

Download


Discussion: Registration Challenges

  • accuracy is the critical criterion here. We need the registration error (residual misalignment) to be smaller than the change we want to measure/detect. Agreement on what constitutes good alignment can therefore vary greatly.
  • the two images have strong differences in coil inhomogeneity. This affects less the registration quality but hampers evaluation. Most of the difference does not become apparent until after registration in direct juxtaposition. Bias field correction beforehand is recommended.
  • we have slightly different voxel sizes
  • if the pathology change is substantial it might affect the registration.

Discussion: Key Strategies

  • the two images have identical contrast, hence we consider "sharper" cost functions, such as NormCorr or MeanSqrd
  • general practice is to register the follow-up to the baseline. However here the follow-up has slightly higher resolution. From an image quality/data perspective it would be better to use the highest resolution image as your fixed/reference. But here follow the most common convention, i.e. fixed image is the baseline.
  • because we seek to assess/quantify regional size change, we must use a rigid (6DOF) scheme, i.e. we must exclude scaling.
  • if the pathology change is soo large that it might affect the registration, we should mask it out. The simplest way to do this is to build a box ROI from the ROItool and feed it as input to the registration. Remember that masking does not mean that masked areas aren't matched, they just do not contribute to the cost function driving the registration, but move along passively. Next more involved level would be to outline the tumor. If a segmentation is available, you can use that.
  • these two images are not too far apart initially, so we reduce the default of expected translational misalignment
  • because accuracy is more important than speed here, we increase the sampling rate from the default 2% to 15%.
  • we also expect minimal differences in scale & distortion: so we can either set the expected values to 0 or run a rigid registration
  • we test the result in areas with good anatomical detail and contrast, far away from the pathology. With rigid body motion a local measure of registration accuracy is representative and can give us a valid limit of detectable change.