Adaptive Graph-Based Total Variation for Tomographic Reconstructions

Sparsity exploiting image reconstruction (SER) methods have been extensively used with total variation (TV) regularization for tomographic reconstructions. Local TV methods fail to preserve texture details and often create additional artifacts due to over-smoothing. Nonlocal TV (NLTV) methods have been proposed as a solution to this but they either lack continuous updates due to computational constraints or limit the locality to a small region. In this letter, we propose adaptive graph-based TV. The proposed method goes beyond spatial similarity between different regions of an image being reconstructed by establishing a connection between similar regions in the entire image regardless of spatial distance. As compared to NLTV, the proposed method is computationally efficient and involves updating the graph prior during every iteration making the connection between similar regions stronger. Moreover, it promotes sparsity in the wavelet and graph gradient domains. Since TV is a special case of graph TV, the proposed method can also be seen as a generalization of SER and TV methods.

[1]  D. Brenner,et al.  Computed tomography--an increasing source of radiation exposure. , 2007, The New England journal of medicine.

[2]  Steve B. Jiang,et al.  4D Computed Tomography Reconstruction from Few-Projection Data via Temporal Non-local Regularization , 2010, MICCAI.

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

[4]  R. Vershynin,et al.  A Randomized Kaczmarz Algorithm with Exponential Convergence , 2007, math/0702226.

[5]  I. M. Otivation Playing with Duality: An Overview of Recent Primal-Dual Approaches for Solving Large-Scale Optimization Problems , 2018 .

[6]  Y. Censor Row-Action Methods for Huge and Sparse Systems and Their Applications , 1981 .

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

[8]  Jinyi Qi,et al.  Iterative reconstruction techniques in emission computed tomography , 2006, Physics in medicine and biology.

[9]  K. P. Kim,et al.  Radiation exposure from CT scans in childhood and subsequent risk of leukaemia and brain tumours: a retrospective cohort study , 2012, The Lancet.

[10]  Steve B. Jiang,et al.  Low-dose CT reconstruction via edge-preserving total variation regularization. , 2010, Physics in medicine and biology.

[11]  P. P. Vaidyanathan,et al.  Extending Classical Multirate Signal Processing Theory to Graphs—Part I: Fundamentals , 2017, IEEE Transactions on Signal Processing.

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

[13]  Heinz H. Bauschke,et al.  Convex Analysis and Monotone Operator Theory in Hilbert Spaces , 2011, CMS Books in Mathematics.

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

[15]  Frédéric Lesage,et al.  The Application of Compressed Sensing for Photo-Acoustic Tomography , 2009, IEEE Transactions on Medical Imaging.

[16]  Reiner Hegerl,et al.  Cryo-electron tomography of biological specimens , 2006, IEEE Signal Processing Magazine.

[17]  Fei Yang,et al.  Compressed magnetic resonance imaging based on wavelet sparsity and nonlocal total variation , 2012, 2012 9th IEEE International Symposium on Biomedical Imaging (ISBI).

[18]  Pierre Vandergheynst,et al.  Graph-based sinogram denoising for tomographic reconstructions , 2016, 2016 38th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC).

[19]  G Bricogne,et al.  Maximum-entropy three-dimensional reconstruction with deconvolution of the contrast transfer function: a test application with adenovirus. , 1996, Journal of structural biology.

[20]  Xiaoqun Zhang,et al.  TICMR: Total Image Constrained Material Reconstruction via Nonlocal Total Variation Regularization for Spectral CT , 2016, IEEE Transactions on Medical Imaging.

[21]  L. Landweber An iteration formula for Fredholm integral equations of the first kind , 1951 .

[22]  Jeffrey A. Fessler,et al.  Fast X-Ray CT Image Reconstruction Using a Linearized Augmented Lagrangian Method With Ordered Subsets , 2014, IEEE Transactions on Medical Imaging.

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

[24]  Eric Todd Quinto,et al.  Electron lambda-tomography , 2009, Proceedings of the National Academy of Sciences.

[25]  D. L. Donoho,et al.  Compressed sensing , 2006, IEEE Trans. Inf. Theory.

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

[27]  J. Frank Electron tomography : methods for three-dimensional visualization of structures in the cell , 2005 .

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

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

[30]  Ozan Öktem,et al.  A componentwise iterated relative entropy regularization method with updated prior and regularization parameter , 2007 .

[31]  A. Brandt Algebraic multigrid theory: The symmetric case , 1986 .

[32]  David G. Lowe,et al.  Scalable Nearest Neighbor Algorithms for High Dimensional Data , 2014, IEEE Transactions on Pattern Analysis and Machine Intelligence.

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

[34]  Pierre Vandergheynst,et al.  Stationary Signal Processing on Graphs , 2016, IEEE Transactions on Signal Processing.

[35]  Rowan Leary,et al.  Compressed sensing electron tomography. , 2013, Ultramicroscopy.

[36]  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.

[37]  L. Tanoue Computed Tomography — An Increasing Source of Radiation Exposure , 2009 .

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

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

[40]  José M. F. Moura,et al.  Big Data Analysis with Signal Processing on Graphs: Representation and processing of massive data sets with irregular structure , 2014, IEEE Signal Processing Magazine.

[41]  Laurent D. Cohen,et al.  Non-local Regularization of Inverse Problems , 2008, ECCV.

[42]  Jie Tang,et al.  Prior image constrained compressed sensing (PICCS): a method to accurately reconstruct dynamic CT images from highly undersampled projection data sets. , 2008, Medical physics.

[43]  Guy Gilboa,et al.  Nonlocal Operators with Applications to Image Processing , 2008, Multiscale Model. Simul..

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

[45]  Stanley Osher,et al.  Image Recovery via Nonlocal Operators , 2010, J. Sci. Comput..

[46]  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.

[47]  Pierre Vandergheynst,et al.  UNLocBoX A matlab convex optimization toolbox using proximal splitting methods , 2014, ArXiv.

[48]  Emil Y Sidky,et al.  Compressive sensing in medical imaging. , 2015, Applied optics.

[49]  Pierre Vandergheynst,et al.  GSPBOX: A toolbox for signal processing on graphs , 2014, ArXiv.

[50]  José M. F. Moura,et al.  Discrete Signal Processing on Graphs , 2012, IEEE Transactions on Signal Processing.

[51]  Jiayu Song,et al.  Sparseness prior based iterative image reconstruction for retrospectively gated cardiac micro-CT. , 2007, Medical physics.

[52]  Y. Censor Finite series-expansion reconstruction methods , 1983, Proceedings of the IEEE.