On the efficiency of proximal methods in CBCT and PET

Cone Beam Computerized Tomography (CBCT) and Positron Emission Tomography (PET) Scans are medical imaging devices that require solving ill-posed inverse problems. The models considered come directly from the physics of the acquisition devices, and take into account the specificity of the (Poisson) noise. We propose various fast numerical schemes to compute the solution. In particular, we show that a new algorithm recently introduced by A. Chambolle and T. Pock is well suited in the PET case when considering non differentiable regularizations such as total variation or wavelet ℓ1-regularization. Numerical experiments indicate that the proposed algorithms compare favorably with respect to well-established methods in tomography.

[1]  J. Pesquet,et al.  Parallel algorithm and hybrid regularization for dynamic PET reconstruction , 2010, IEEE Nuclear Science Symposuim & Medical Imaging Conference.

[2]  Jeffrey A. Fessler,et al.  Fast Monotonic Algorithms for Transmission Tomography , 1999, IEEE Trans. Medical Imaging.

[3]  Alin Achim,et al.  18th IEEE International Conference on Image Processing, ICIP 2011, Brussels, Belgium, September 11-14, 2011 , 2011, ICIP.

[4]  Patrick L. Combettes,et al.  Signal Recovery by Proximal Forward-Backward Splitting , 2005, Multiscale Model. Simul..

[5]  Mohamed-Jalal Fadili,et al.  A Proximal Iteration for Deconvolving Poisson Noisy Images Using Sparse Representations , 2008, IEEE Transactions on Image Processing.

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

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

[8]  Nelly Pustelnik,et al.  Nested Iterative Algorithms for Convex Constrained Image Recovery Problems , 2008, SIAM J. Imaging Sci..

[9]  Antonin Chambolle,et al.  A First-Order Primal-Dual Algorithm for Convex Problems with Applications to Imaging , 2011, Journal of Mathematical Imaging and Vision.

[10]  丸山 徹 Convex Analysisの二,三の進展について , 1977 .

[11]  Hakan Erdogan,et al.  Monotonic algorithms for transmission tomography , 1999, IEEE Transactions on Medical Imaging.

[12]  Alfred O. Hero,et al.  Ieee Transactions on Image Processing: to Appear Penalized Maximum-likelihood Image Reconstruction Using Space-alternating Generalized Em Algorithms , 2022 .

[13]  Eero P. Simoncelli,et al.  Image quality assessment: from error visibility to structural similarity , 2004, IEEE Transactions on Image Processing.

[14]  L. Feldkamp,et al.  Practical cone-beam algorithm , 1984 .

[15]  Rebecca Willett,et al.  This is SPIRAL-TAP: Sparse Poisson Intensity Reconstruction ALgorithms—Theory and Practice , 2010, IEEE Transactions on Image Processing.