BRAINSFit to simpleITK
Contents
Goal
We are developing registration module in Slicer version 4 (which is using ITKv4) for deformable registration of prostate MRI.
We want to develop a SimpleITK prostate registration tool in Slicer4/ITKv4. This Slicer4/ITKv4 registration tool should be functional, accurate and fast (i.e., comparable with the the functionality we had in Slicer3/ITKv3, which has been evaluated previously [1]).
Details on the registration approach
Registration is applied to align preprocedural and intraoperational MR T2 image volumes. We are using masks the prostate for both image data sets. Registration is done using MMI metric with rigid, affine and B-spline stages applied in sequence. In Slicer3/BRAINSFit we use gradient descent for rigid/affine, and LBFGS for B-spline.
Parameters we are using to call BRAINSFit in Slicer 3.6 can be found here.
Sample data can be found here.
A (very simplified) graphical overview of BRAINSFit registration can be found here.
Current status
Registration code on github: here
Latest update: April 28, 2015
Things that are working
- exhaustive search based initialization procedure implemented and produces good results (cross-correlation metric)
Things that are implemented, but are not working
- SimpleITK code for rigid registration
Things that are not implemented, but need to be implemented
- SimpleITK code for affine registration
- SimpleITK code for BSpline registration
Unresolved issues out of our direct control
Standing SimpleITK issues
- no API for consistent initialization of the metric seed
Standing ITKv4 issues
- capability to sample voxels only within the mask (see topic under review on Gerrit: http://review.source.kitware.com/#/c/19684/)
Standing Slicer4 issues
References
[1] Fedorov A, Tuncali K, Fennessy FM, Tokuda J, Hata N, et al. (2012) Image registration for targeted MRI-guided transperineal prostate biopsy. J Magn Reson Imaging 36: 987–992. Available: http://dx.doi.org/10.1002/jmri.23688.