2011 Winter Project Week:Intra-ProceduralProstateMotion - Revision history

Lasso: /* Next steps */

2011-01-25T00:15:33Z

Next steps

Lasso: /* Key Investigators */

2011-01-25T00:12:53Z

Key Investigators

Lasso: /* Next steps */

2011-01-25T00:12:16Z

Next steps

Lasso: /* Bugs */

2011-01-14T17:17:35Z

Bugs

Fedorov: /* Results */

2011-01-14T16:44:41Z

Results

Lasso: /* Queen's registration algo steps: */

2011-01-14T15:13:28Z

Queen's registration algo steps:

Fedorov at 05:54, 14 January 2011

2011-01-14T05:54:05Z

Fedorov at 05:53, 14 January 2011

2011-01-14T05:53:17Z

Fedorov at 05:51, 14 January 2011

2011-01-14T05:51:15Z

Fedorov at 05:50, 14 January 2011

2011-01-14T05:50:43Z

@@ Line 177: / Line 177: @@
 *Slice-to-volume registration compensates large part of the patient motion (by about 50-80%)
 *The registration doesn't compensate prostate motion. To compensate that the registration a second registration step should be performed, that uses a smaller ROI (containing only the prostate). A third step could be to perform non-rigid registration of the prostate.
 ==References==

@@ Line 50: / Line 50: @@
 <h3>Dissemination</h3>
-http://viewvc.slicer.org/viewcvs.cgi/trunk/SLC2011-ProstateRegistration/?root=NAMICSandBox
+http://viewvc.slicer.org/viewcvs.cgi/trunk/IGTLoadableModules/ProstateRegistration/?root=NAMICSandBox
 ==Data==

@@ Line 181: / Line 181: @@
 *Add CTest test cases
 *Add s3ext file to Slicer trunk?
 *Add further registration steps (rigid registration with prostate ROI, deformable registration with prostate ROI)

@@ Line 70: / Line 70: @@
 *TP11: Case011
-==Bugs==
+==Found bugs==
 *itkOptImageToImageMetric gives different result if executed on single vs. multiple threads.
-*ImageSpatialObject< TDimension,  PixelType >::ValueAt evaluates the value at a different position than ImageMaskSpatialObject< TDimension >::IsInside (which determines if the point is inside or not)
+*ImageSpatialObject< TDimension,  PixelType >::ValueAt evaluates the value at a different position than ImageMaskSpatialObject< TDimension >::IsInside (which determines if the point is inside or not) => [http://public.kitware.com/Bug/view.php?id=11699 0011699: ImageMaskSpatialObject IsInside and ValueAt methods give inconsistent results]
 ==BWH slice-volume rigid registration==

@@ Line 102: / Line 102: @@
 |-
 | align="left" |
-| align="left" | '''Tfm magn'''
+| align="left" | '''Mis-registration magnitude recovered by volume-to-volume registration'''
 | align="center | '''Axial slice'''
 | align="center" | '''Axial+Cor slice'''
 | align="center" | '''Axial+Cor+Sag slice'''
-| align="left" | '''Tfm magn'''
+| align="left" | '''Mis-registration magnitude recovered by volume-to-volume registration'''
 | align="center" | '''Axial slice'''
 | align="center" | '''Axial+Cor slice'''

@@ Line 78: / Line 78: @@
 ==Queen's registration algo steps:==
-*rescale the fixed image to have the same min/max intensity as the moving image
+*Rescale the fixed image to have the same min/max intensity as the moving image
-*reads in the GT deformation field
+*Set up fixed image regions: VOI is only the prostate
-*set up fixed image regions based on VOI params
+*Set up all supported metrics: MeanSquaresSparseVolumeToVolumeMetric + NC2SparseVolumeToVolumeMetric + ViolaWellsSparseVolumeToVolumeMetric
-*set up metrics: MeanSquaresSparseVolumeToVolumeMetric + NC2SparseVolumeToVolumeMetric + ViolaWellsSparseVolumeToVolumeMetric
+*Set multiple fixed image regions (one for each slice) for all the metrics
-*set multiple fixed image regions for all the metrics
+*Register: CMAEvolutionStrategyOptimizer, rigid translation, mutual information metric; CenteredTransformInitializer; both volumes smoothed
-*register: CMAEvolutionStrategyOptimizer, rigid translation, mutual information metric; CenteredTransformInitializer; both volumes smoothed
+*Apply determined transform
-*apply determined transform
+*Register: RegularStepGradientDescentOptimizer (constrain translation), Euler 3D transform, mutual information metric; CenteredTransformInitializer; both volumes smoothed
-*register: RegularStepGradientDescentOptimizer (constrain translation), Euler 3D transform, mutual information metric; *CenteredTransformInitializer; both volumes smoothed
+*Apply determined transform
-*apply determined transform
+*Register: MeanSquaresSparseVolumeToVolumeMetric/NC2SparseVolumeToVolumeMetric/ViolaWellsSparseVolumeToVolumeMetric, RegularStepGradientDescentOptimizer, affine transform,
-*register: MeanSquaresSparseVolumeToVolumeMetric/NC2SparseVolumeToVolumeMetric/ViolaWellsSparseVolumeToVolumeMetric, RegularStepGradientDescentOptimizer, affine transform,
+*Apply determined transform
-*apply determined transform
+*Register: bspline transform 5 nodes, same metric, GD optimizer
-*register: bspline transform 5 nodes, same metric, GD optimizer
+*Register: bspline transform 15 nodes, same metric, GD optimizer
 ==Results==

@@ Line 6: / Line 6: @@
 Image:Pelvis_and_prostate_motion_example.gif|Motion of the pelvis and prostate gland over the first 30 minutes of MR-guided transperineal biopsy procedure.
 Image:TransrectalBiopsySagittal.png|Sagittal view of a transrectal MR-guided transrectal biopsy procedure.
-Image:prostate_unreg.gif|[[Unregistered]]
+Image:prostate_unreg.gif|Unregistered
-Image:prostate_rega.gif|[[Registered w/ axial slice]]
+Image:prostate_rega.gif|Registered w/ axial slice
-Image:prostate_regac.gif|[[Registered w/ axial+cor slices]]
+Image:prostate_regac.gif|Registered w/ axial+cor slices
-Image:prostate_regacs.gif|[[Registered w/ axial+cor+sag slices]]
+Image:prostate_regacs.gif|Registered w/ axial+cor+sag slices
-Image:prostate_regvol.gif|[[Registered w/ full volume]]
+Image:prostate_regvol.gif|Registered w/ full volume
 </gallery>