Metal artifact reduction in CT-based attenuation correction of PET using sobolev sinogram restoration

In patients bearing metallic implants, CT images reconstructed by the filtered back-projection algorithm usually suffer substantially from streaking metal artifacts. The CT-based attenuation correction of PET using such images can lead to pseudo-uptakes and thus equivocal findings. In this paper, we introduce a new metal artifact reduction (MAR) algorithm based on Bayesian iterative restoration techniques applied in the sinogram space. We proposed a Sobolev prior in the maximum a posteriori (MAP) estimation of the projections corrupted by metallic implants. The Sobolev prior was invoked in order to impose the a priori knowledge that a CT sinogram is a smooth dataset in which it is highly probable that neighboring projections have similar photon counts. We compared the proposed prior with a total variation (TV) one which imposes piece-wise smoothness on the sinograms being restored. We also compared it with a smoothed TV (TVs) prior which ranks between the Sobolev and TV ones. We formulated the MAP estimation as a convex constrained optimization problem and solved it for the Sobolev and TVs priors by a projected gradient descent algorithm and for the TV/s priors by a sophisticated primal-dual projected gradient algorithm. The results of artifact simulations in a real CT image showed that the Sobolev-based MAR algorithm outperforms its TV and TVs-based counterparts in terms of convergence rate and is comparable with the TVs in projection recovery. It was demonstrated that the proposed MAR algorithm has high applicability in fast and efficient CT-based attenuation correction of PET data.

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

[2]  Wenyuan Xu,et al.  Behavioral analysis of anisotropic diffusion in image processing , 1996, IEEE Trans. Image Process..

[3]  O. Nelles,et al.  An Introduction to Optimization , 1996, IEEE Antennas and Propagation Magazine.

[4]  Shiying Zhao,et al.  A wavelet method for metal artifact reduction with multiple metallic objects in the field of view , 2002 .

[5]  ANTONIN CHAMBOLLE,et al.  An Algorithm for Total Variation Minimization and Applications , 2004, Journal of Mathematical Imaging and Vision.

[6]  Yurii Nesterov,et al.  Introductory Lectures on Convex Optimization - A Basic Course , 2014, Applied Optimization.

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

[8]  Stephen P. Boyd,et al.  Convex Optimization , 2004, Algorithms and Theory of Computation Handbook.

[9]  Li Zhang,et al.  Metal artifact reduction in CT images by sinogram TV inpainting , 2008 .

[10]  Zhiqiang Chen,et al.  Metal artifact reduction in dual energy CT by sinogram segmentation based on active contour model and TV inpainting , 2009, 2009 IEEE Nuclear Science Symposium Conference Record (NSS/MIC).

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

[12]  Denis Laurendeau,et al.  39 oral: An Algorithm for Efficient Metal Artifact Reductions in Permanent Seed Implants , 2009 .

[13]  H. Zaidi,et al.  Qualitative and Quantitative Assessment of Metal Artifacts Arising from Implantable Cardiac Pacing Devices in Oncological PET/CT Studies: A Phantom Study , 2011, Molecular Imaging and Biology.

[14]  Jacob Geleijns,et al.  Development and validation of segmentation and interpolation techniques in sinograms for metal artifact suppression in CT. , 2010, Medical physics.

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

[16]  Marc Teboulle,et al.  Gradient-based algorithms with applications to signal-recovery problems , 2010, Convex Optimization in Signal Processing and Communications.

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

[18]  Denis Laurendeau,et al.  An algorithm for efficient metal artifact reductions in permanent seed. , 2010, Medical physics.

[19]  L. Xing,et al.  Metal artifact reduction in x-ray computed tomography (CT) by constrained optimization. , 2011, Medical physics.