Simultaneous reconstruction and segmentation for dynamic SPECT imaging
- 5 August 2016
- journal article
- research article
- Published by IOP Publishing in Inverse Problems
- Vol. 32 (10) , 104002
Abstract
This work deals with the reconstruction of dynamic images that incorporate characteristic dynamics in certain subregions, as arising for the kinetics of many tracers in emission tomography (SPECT, PET). We make use of a basis function approach for the unknown tracer concentration by assuming that the region of interest can be divided into subregions with spatially constant concentration curves. Applying a regularised variational framework reminiscent of the Chan-Vese model for image segmentation we simultaneously reconstruct both the labelling functions of the subregions as well as the subconcentrations within each region. Our particular focus is on applications in SPECT with the Poisson noise model, resulting in a Kullback–Leibler data fidelity in the variational approach. We present a detailed analysis of the proposed variational model and prove existence of minimisers as well as error estimates. The latter apply to a more general class of problems and generalise existing results in literature since we deal with a nonlinear forward operator and a nonquadratic data fidelity. A computational algorithm based on alternating minimisation and splitting techniques is developed for the solution of the problem and tested on appropriately designed synthetic data sets. For those we compare the results to those of standard EM reconstructions and investigate the effects of Poisson noise in the data.Keywords
All Related Versions
This publication has 18 references indexed in Scilit:
- Simultaneous Reconstruction and Segmentation of Dynamic PET via Low-Rank and Sparse Matrix DecompositionIEEE Transactions on Biomedical Engineering, 2015
- On convergence rates for iteratively regularized Newton-type methods under a Lipschitz-type nonlinearity conditionJournal of Inverse and Ill-posed Problems, 2014
- Convergence rates in expectation for Tikhonov-type regularization of inverse problems with Poisson dataInverse Problems, 2012
- A First-Order Primal-Dual Algorithm for Convex Problems with Applications to ImagingJournal of Mathematical Imaging and Vision, 2010
- Dynamic single photon emission computed tomography—basic principles and cardiac applicationsPhysics in Medicine & Biology, 2010
- A Solver for Dynamic PET Reconstructions based on Forward-Backward-SplittingAIP Conference Proceedings, 2010
- A General Framework for a Class of First Order Primal-Dual Algorithms for Convex Optimization in Imaging ScienceSIAM Journal on Imaging Sciences, 2010
- Error estimates for non-quadratic regularization and the relation to enhancementInverse Problems, 2006
- Mean and covariance properties of dynamic PET reconstructions from list-mode dataIEEE Transactions on Medical Imaging, 2005
- Kinetic parameter estimation from SPECT cone-beam projection measurementsPhysics in Medicine & Biology, 1998