An algorithm for total variation regularization in high-dimensional linear problems

This paper describes an iterative algorithm for high-dimensional linear inverse problems, which is regularized by a differentiable discrete approximation of the total variation (TV) penalty. The algorithm is an interlaced iterative method based on optimization transfer with a separable quadratic surrogate for the TV penalty. The surrogate cost function is optimized using the block iterative regularized algebraic reconstruction technique (RSART). A proof of convergence is given and convergence is illustrated by numerical experiments with simulated parallel-beam computerized tomography (CT) data. The proposed method provides a block-iterative and convergent, hence efficient and reliable, algorithm to investigate the effects of TV regularization in applications such as CT.

[1]  Gabor T. Herman,et al.  A Storage-Efficient Algorithm for Finding the Regularized Solution of a Large, Inconsistent System of Equations , 1980 .

[2]  A. Lent,et al.  Iterative algorithms for large partitioned linear systems, with applications to image reconstruction , 1981 .

[3]  P. Joseph.  An Improved Algorithm for Reprojecting Rays through Pixel Images , 1982 .

[4]  A. Kak,et al.  Simultaneous Algebraic Reconstruction Technique (SART): A Superior Implementation of the Art Algorithm , 1984, Ultrasonic imaging.

[5]  M. Daube-Witherspoon,et al.  An Iterative Image Space Reconstruction Algorthm Suitable for Volume ECT , 1986, IEEE Transactions on Medical Imaging.

[6]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[7]  A. R. De Pierro,et al.  On the relation between the ISRA and the EM algorithm for positron emission tomography , 1993, IEEE Trans. Medical Imaging.

[8]  H. Malcolm Hudson,et al.  Accelerated image reconstruction using ordered subsets of projection data , 1994, IEEE Trans. Medical Imaging.

[9]  Klaus Mueller,et al.  The weighted-distance scheme: a globally optimizing projection ordering method for ART , 1997, IEEE Transactions on Medical Imaging.

[10]  G. Kontaxakis,et al.  Optimized image reconstruction for emission tomography using ordered subsets, median root prior and a Web-based interface , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

[11]  Chung-Ming Chen,et al.  Cross-reference weighted least square estimates for positron emission tomography , 1998, IEEE Transactions on Medical Imaging.

[12]  Gengsheng L. Zeng,et al.  Total variation regulated EM algorithm [SPECT reconstruction] , 1999 .

[13]  M. Persson,et al.  Total variation norm for three-dimensional iterative reconstruction in limited view angle tomography , 2001, Physics in medicine and biology.

[14]  I. Daubechies,et al.  An iterative thresholding algorithm for linear inverse problems with a sparsity constraint , 2003, math/0307152.

[15]  Wojciech Zbijewski,et al.  Suppression of intensity transition artifacts in statistical x-ray computer tomography reconstruction through radon inversion initialization. , 2003, Medical physics.

[16]  B. De Man,et al.  Distance-driven projection and backprojection in three dimensions. , 2004, Physics in medicine and biology.

[17]  E. Sidky,et al.  Accurate image reconstruction from few-views and limited-angle data in divergent-beam CT , 2009, 0904.4495.

[18]  Jean-Baptiste Thibault,et al.  A three-dimensional statistical approach to improved image quality for multislice helical CT. , 2007, Medical physics.

[19]  Olaf Dössel,et al.  An Optimal Radial Profile Order Based on the Golden Ratio for Time-Resolved MRI , 2007, IEEE Transactions on Medical Imaging.

[20]  Guang-Hong Chen,et al.  Limited view angle tomographic image reconstruction via total variation minimization , 2007, SPIE Medical Imaging.

[21]  Yair Censor,et al.  On Diagonally Relaxed Orthogonal Projection Methods , 2007, SIAM J. Sci. Comput..

[22]  G T Herman,et al.  Image reconstruction from a small number of projections , 2008, Inverse problems.

[23]  E. Sidky,et al.  Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization , 2008, Physics in medicine and biology.

[24]  José M. Bioucas-Dias,et al.  Adaptive total variation image deblurring: A majorization-minimization approach , 2009, Signal Process..

[25]  Hengyong Yu,et al.  Compressed sensing based interior tomography , 2009, Physics in medicine and biology.

[26]  M. Vannier,et al.  Why do commercial CT scanners still employ traditional, filtered back-projection for image reconstruction? , 2009, Inverse problems.

[27]  Jie Tang,et al.  Performance comparison between total variation (TV)-based compressed sensing and statistical iterative reconstruction algorithms , 2009, Physics in medicine and biology.

[28]  Marc Teboulle,et al.  Fast Gradient-Based Algorithms for Constrained Total Variation Image Denoising and Deblurring Problems , 2009, IEEE Transactions on Image Processing.

[29]  Ruijiang Li,et al.  GPU-based Cone Beam CT Reconstruction via Total Variation Regularization , 2010 .

[30]  Steve B. Jiang,et al.  GPU-based fast cone beam CT reconstruction from undersampled and noisy projection data via total variation , 2010 .

[31]  A B Rosenfeld,et al.  Total variation superiorization schemes in proton computed tomography image reconstruction. , 2010, Medical physics.

[32]  M. Kachelriess,et al.  Improved total variation-based CT image reconstruction applied to clinical data , 2011, Physics in medicine and biology.