Accelerating spatially non-uniform update for sparse target recovery in fluorescence molecular tomography by ordered subsets and momentum methods

Fluorescence molecular tomography (FMT) is a significant preclinical imaging modality that has been actively studied in the past two decades. However, it remains a challenging task to obtain fast and accurate reconstruction of fluorescent probe distribution in small animals due to the large computational burden and the ill-posed nature of the inverse problem. We have recently studied a non-uniform multiplicative updating algorithm, and obtained some further speed gain with the ordered subsets (OS) method. However, increasing the number of OS leads to larger approximation errors and the speed gain from larger number of OS is marginal. In this paper, we propose to further enhance the convergence speed by incorporating the first order momentum method that uses previous iterations to achieve a quadratic convergence rate. Using cubic phantom experiment, we have shown that the proposed method indeed leads to a much faster convergence.

[1]  Simon R Cherry,et al.  In vivo molecular and genomic imaging: new challenges for imaging physics. , 2004, Physics in medicine and biology.

[2]  Vasilis Ntziachristos,et al.  Looking and listening to light: the evolution of whole-body photonic imaging , 2005, Nature Biotechnology.

[3]  Jeffrey A. Fessler,et al.  Combining Ordered Subsets and Momentum for Accelerated X-Ray CT Image Reconstruction , 2015, IEEE Transactions on Medical Imaging.

[4]  Changqing Li,et al.  Nonconvex regularizations in fluorescence molecular tomography for sparsity enhancement , 2014, Physics in medicine and biology.

[5]  R. Leahy,et al.  Joint L1 and total variation regularization for fluorescence molecular tomography , 2012, Physics in medicine and biology.

[6]  Emmanuel J. Candès,et al.  NESTA: A Fast and Accurate First-Order Method for Sparse Recovery , 2009, SIAM J. Imaging Sci..

[7]  Hakan Erdogan,et al.  Ordered subsets algorithms for transmission tomography. , 1999, Physics in medicine and biology.

[8]  Hanli Liu,et al.  Optimization of probe geometry for diffuse optical brain imaging based on measurement density and distribution. , 2009, Applied optics.

[9]  Yue Zhao,et al.  Comparison of Regularization Methods in Fluorescence Molecular Tomography , 2014 .

[10]  Yurii Nesterov,et al.  Smooth minimization of non-smooth functions , 2005, Math. Program..

[11]  Changqing Li,et al.  Nonuniform update for sparse target recovery in fluorescence molecular tomography accelerated by ordered subsets. , 2014, Biomedical optics express.

[12]  K. Lange The MM Algorithm , 2013 .

[13]  Changqing Li,et al.  Simultaneous PET and Multispectral 3-Dimensional Fluorescence Optical Tomography Imaging System , 2011, The Journal of Nuclear Medicine.

[14]  Jeffrey A. Fessler,et al.  Accelerating Ordered Subsets Image Reconstruction for X-ray CT Using Spatially Nonuniform Optimization Transfer , 2013, IEEE Transactions on Medical Imaging.

[15]  Jeffrey A. Fessler,et al.  Grouped-coordinate ascent algorithms for penalized-likelihood transmission image reconstruction , 1997, IEEE Transactions on Medical Imaging.

[16]  Zhou Yu,et al.  Fast Model-Based X-Ray CT Reconstruction Using Spatially Nonhomogeneous ICD Optimization , 2011, IEEE Transactions on Image Processing.

[17]  Alvaro R. De Pierro,et al.  A modified expectation maximization algorithm for penalized likelihood estimation in emission tomography , 1995, IEEE Trans. Medical Imaging.

[18]  R. Leahy,et al.  A three-dimensional multispectral fluorescence optical tomography imaging system for small animals based on a conical mirror design. , 2009, Optics express.

[19]  Simon R. Arridge,et al.  Inverse Methods for Optical Tomography , 1993, IPMI.

[20]  J. Laible,et al.  Coupled complex adjoint sensitivities for frequency-domain fluorescence tomography: theory and vectorized implementation , 2003 .

[21]  Michael Unser,et al.  An Efficient Numerical Method for General $L_{p}$ Regularization in Fluorescence Molecular Tomography , 2010, IEEE Transactions on Medical Imaging.

[22]  Guobao Wang,et al.  Three-dimensional fluorescence optical tomography in small-animal imaging using simultaneous positron-emission-tomography priors. , 2009, Optics letters.

[23]  Y. Nesterov A method for solving the convex programming problem with convergence rate O(1/k^2) , 1983 .

[24]  L. R. Dice Measures of the Amount of Ecologic Association Between Species , 1945 .