TGV-regularized inversion of the Radon transform for photoacoustic tomography

We propose and study a reconstruction method for photoacoustic tomography (PAT) based on total generalized variation (TGV) regularization for the inversion of the slice-wise 2D-Radon transform in 3D. The latter problem occurs for recently-developed PAT imaging techniques with parallelized integrating ultrasound detection where projection data from various directions is sequentially acquired. As the imaging speed is presently limited to 20 seconds per 3D image, the reconstruction of temporally-resolved 3D sequences of, e.g., one heartbeat or breathing cycle, is very challenging and currently, the presence of motion artifacts in the reconstructions obstructs the applicability for biomedical research. In order to push these techniques forward towards real time, it thus becomes necessary to reconstruct from less measured data such as few-projection data and consequently, to employ sophisticated reconstruction methods in order to avoid typical artifacts. The proposed TGV-regularized Radon inversion is a variational method that is shown to be capable of such artifact-free inversion. It is validated by numerical simulations, compared to filtered back projection (FBP), and performance-tested on real data from phantom as well as in-vivo mouse experiments. The results indicate that a speed-up factor of four is possible without compromising reconstruction quality.

[1]  Kristian Bredies,et al.  Fast quantitative susceptibility mapping using 3D EPI and total generalized variation , 2015, NeuroImage.

[2]  Robert Nuster,et al.  Piezoelectric line detector array for photoacoustic tomography , 2017, Photoacoustics.

[3]  M. Haltmeier,et al.  Temporal back-projection algorithms for photoacoustic tomography with integrating line detectors , 2007 .

[4]  H. Goluba,et al.  Eigenvalue computation in the 20 th century Gene , 2000 .

[5]  Kristian Bredies,et al.  Recovering Piecewise Smooth Multichannel Images by Minimization of Convex Functionals with Total Generalized Variation Penalty , 2011, Efficient Algorithms for Global Optimization Methods in Computer Vision.

[6]  T. Pock,et al.  Second order total generalized variation (TGV) for MRI , 2011, Magnetic resonance in medicine.

[7]  Lihong V. Wang,et al.  A practical guide to photoacoustic tomography in the life sciences , 2016, Nature Methods.

[8]  Kristian Bredies,et al.  Total Generalized Variation in Diffusion Tensor Imaging , 2013, SIAM J. Imaging Sci..

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

[10]  G. Golub,et al.  Eigenvalue computation in the 20th century , 2000 .

[11]  ON THE INVERSE RADON TRANSFORM , 2007 .

[12]  Philip Kollmannsberger,et al.  Architecture of the osteocyte network correlates with bone material quality , 2013, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[13]  Jonas Adler,et al.  Solving ill-posed inverse problems using iterative deep neural networks , 2017, ArXiv.

[14]  Otmar Scherzer,et al.  A Reconstruction Algorithm for Photoacoustic Imaging Based on the Nonuniform FFT , 2009, IEEE Transactions on Medical Imaging.

[15]  Thomas Berer,et al.  All-optical photoacoustic projection imaging. , 2017, Biomedical optics express.

[16]  K. Bredies,et al.  Infimal convolution of total generalized variation functionals for dynamic MRI , 2017, Magnetic resonance in medicine.

[17]  Vasilis Ntziachristos,et al.  Looking at sound: optoacoustics with all-optical ultrasound detection , 2018, Light: Science & Applications.

[18]  Alston S. Householder,et al.  The Theory of Matrices in Numerical Analysis , 1964 .

[19]  THE THEORY OF MATRICES IN NUMERICAL ANALYSIS , 1965 .

[20]  Tony F. Chan,et al.  High-Order Total Variation-Based Image Restoration , 2000, SIAM J. Sci. Comput..

[21]  Marta Betcke,et al.  Accelerated high-resolution photoacoustic tomography via compressed sensing , 2016, Physics in medicine and biology.

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

[23]  Graptor , 2020, Proceedings of the 34th ACM International Conference on Supercomputing.

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

[25]  Robert Nuster,et al.  High resolution three-dimensional photoacoutic tomography with CCD-camera based ultrasound detection , 2014, Biomedical optics express.

[26]  Sanjiv S. Gambhir,et al.  Photoacoustic clinical imaging , 2019, Photoacoustics.

[27]  Lihong V. Wang,et al.  Photoacoustic tomography: principles and advances. , 2016, Electromagnetic waves.

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

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

[30]  Feng Gao,et al.  Enhancing sparse-view photoacoustic tomography with combined virtually parallel projecting and spatially adaptive filtering , 2018, Biomedical optics express.

[31]  Kristian Bredies,et al.  Spatio-temporal TGV denoising for ASL perfusion imaging , 2017, NeuroImage.

[32]  B T Cox,et al.  Fast calculation of pulsed photoacoustic fields in fluids using k-space methods. , 2005, The Journal of the Acoustical Society of America.

[33]  Xiaochuan Pan,et al.  A constrained, total-variation minimization algorithm for low-intensity x-ray CT. , 2010, Medical physics.

[34]  H. Weber,et al.  Temporal backward projection of optoacoustic pressure transients using fourier transform methods. , 2001, Physics in medicine and biology.

[35]  Kristian Bredies,et al.  Total generalized variation regularization for multi-modal electron tomography. , 2019, Nanoscale.

[36]  M. Haltmeier,et al.  Model‐based time reversal method for photoacoustic imaging of heterogeneous media , 2008 .

[37]  K. Bredies,et al.  Joint MR-PET Reconstruction Using a Multi-Channel Image Regularizer. , 2017, IEEE transactions on medical imaging.

[38]  Gaohang Yu,et al.  Sparse-view x-ray CT reconstruction via total generalized variation regularization , 2014, Physics in medicine and biology.

[39]  Yuanyuan Wang,et al.  A photoacoustic image reconstruction method using total variation and nonconvex optimization , 2014, Biomedical engineering online.

[40]  Heather K. Hunt,et al.  Hand-held optoacoustic imaging: A review , 2018, Photoacoustics.

[41]  Stephan Antholzer,et al.  Deep null space learning for inverse problems: convergence analysis and rates , 2018, Inverse Problems.

[42]  P. Lions,et al.  Image recovery via total variation minimization and related problems , 1997 .

[43]  C. S. Sastry,et al.  Reconstruction of sparse-view tomography via preconditioned Radon sensing matrix , 2019 .

[44]  Markus Haltmeier,et al.  Full field detection in photoacoustic tomography. , 2010, Optics express.

[45]  Karl Kunisch,et al.  Total Generalized Variation , 2010, SIAM J. Imaging Sci..

[46]  Huabei Jiang,et al.  Photoacoustic image reconstruction from few-detector and limited-angle data , 2011, Biomedical optics express.

[47]  Lihong V. Wang,et al.  Photoacoustic Tomography: In Vivo Imaging from Organelles to Organs , 2012, Science.

[48]  Günther Paltauf,et al.  Comparison of Piezoelectric and Optical Projection Imaging for Three-Dimensional In Vivo Photoacoustic Tomography , 2019, J. Imaging.

[49]  KEIJO HÄMÄLÄINEN,et al.  Sparse Tomography , 2013, SIAM J. Sci. Comput..

[50]  Kristian Bredies,et al.  Joint MR-PET Reconstruction Using a Multi-Channel Image Regularizer , 2017, IEEE Transactions on Medical Imaging.

[51]  R. Huber.  Variational Regularization for Systems of Inverse Problems , 2019, BestMasters.

[52]  P. Beard.  Biomedical photoacoustic imaging , 2011, Interface Focus.

[53]  P. Burgholzer,et al.  Thermoacoustic tomography with integrating area and line detectors , 2005, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[54]  M. Haltmeier,et al.  Exact and approximative imaging methods for photoacoustic tomography using an arbitrary detection surface. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[55]  G Paltauf,et al.  Weight factors for limited angle photoacoustic tomography , 2009, Physics in medicine and biology.

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

[57]  K. Bredies,et al.  Regularization of linear inverse problems with total generalized variation , 2014 .

[58]  Xiaochuan Pan,et al.  A constrained, total-variation minimization algorithm for low-intensity x-ray CT. , 2011, Medical physics.

[59]  Yiqiu Dong,et al.  An algorithm for total variation regularized photoacoustic imaging , 2015, Adv. Comput. Math..

[60]  Wiendelt Steenbergen,et al.  A framework for directional and higher-order reconstruction in photoacoustic tomography , 2017, Physics in medicine and biology.

[61]  Li Bai,et al.  Denoising optical coherence tomography using second order total generalized variation decomposition , 2016, Biomed. Signal Process. Control..

[62]  Markus Haltmeier,et al.  Experimental evaluation of reconstruction algorithms for limited view photoacoustic tomography with line detectors , 2007 .

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

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