Graph-based sinogram denoising for tomographic reconstructions

Limited data and low-dose constraints are common problems in a variety of tomographic reconstruction paradigms, leading to noisy and incomplete data. Over the past few years, sinogram denoising has become an essential preprocessing step for low-dose Computed Tomographic (CT) reconstructions. We propose a novel sinogram denoising algorithm inspired by signal processing on graphs. Graph-based methods often perform better than standard filtering operations since they can exploit the signal structure. This makes the sinogram an ideal candidate for graph based denoising since it generally has a piecewise smooth structure. We test our method with a variety of phantoms using different reconstruction methods. Our numerical study shows that the proposed algorithm improves the performance of analytical filtered back-projection (FBP) and iterative methods such as ART (Kaczmarz), and SIRT (Cimmino). We observed that graph denoised sinograms always minimize the error measure and improve the accuracy of the solution, compared to regular reconstructions.

[1]  Per Christian Hansen,et al.  AIR Tools - A MATLAB package of algebraic iterative reconstruction methods , 2012, J. Comput. Appl. Math..

[2]  F. Natterer.  The Mathematics of Computerized Tomography , 1986 .

[3]  P. Gilbert.  Iterative methods for the three-dimensional reconstruction of an object from projections. , 1972, Journal of theoretical biology.

[4]  Å. Björck,et al.  Accelerated projection methods for computing pseudoinverse solutions of systems of linear equations , 1979 .

[5]  Patrick L. Combettes,et al.  Proximal Splitting Methods in Signal Processing , 2009, Fixed-Point Algorithms for Inverse Problems in Science and Engineering.

[6]  B. F. Logan,et al.  The Fourier reconstruction of a head section , 1974 .

[7]  Patrick J. La Riviere.  Penalized‐likelihood sinogram smoothing for low‐dose CT , 2005 .

[8]  Daniel Kolditz,et al.  Iterative reconstruction methods in X-ray CT. , 2012, Physica medica : PM : an international journal devoted to the applications of physics to medicine and biology : official journal of the Italian Association of Biomedical Physics.

[9]  A. Tsybakov,et al.  Minimax theory of image reconstruction , 1993 .

[10]  Jose-Jesus Fernandez,et al.  Computational methods for electron tomography. , 2012, Micron.

[11]  Pascal Frossard,et al.  The emerging field of signal processing on graphs: Extending high-dimensional data analysis to networks and other irregular domains , 2012, IEEE Signal Processing Magazine.

[12]  J. Fessler.  Statistical Image Reconstruction Methods for Transmission Tomography , 2000 .

[13]  Patrick J. La Rivière,et al.  Penalized-likelihood sinogram restoration for computed tomography , 2006, IEEE Transactions on Medical Imaging.

[14]  Zhengrong Liang,et al.  An experimental study on the noise properties of x-ray CT sinogram data in Radon space , 2008, Physics in medicine and biology.

[15]  G. Herman,et al.  Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and x-ray photography. , 1970, Journal of theoretical biology.

[16]  R A Brooks,et al.  Theory of image reconstruction in computed tomography. , 1975, Radiology.

[17]  Jiang Hsieh,et al.  Computed Tomography: Principles, Design, Artifacts, and Recent Advances, Fourth Edition , 2022 .

[18]  Zhengrong Liang,et al.  Sinogram noise reduction for low-dose CT by statistics-based nonlinear filters , 2005, SPIE Medical Imaging.

[19]  Rabab Kreidieh Ward,et al.  A sinogram denoising algorithm for low-dose computed tomography , 2016, BMC Medical Imaging.

[20]  Rebecca S Lewis,et al.  Projected cancer risks from computed tomographic scans performed in the United States in 2007. , 2009, Archives of internal medicine.

[21]  P. L. La Riviere.  Penalized-likelihood sinogram smoothing for low-dose CT. , 2005, Medical physics.

[22]  G. Moneta,et al.  Projected Cancer Risks From Computed Tomographic Scans Performed in the United States in 2007 , 2010 .