Accelerating image reconstruction in three-dimensional optoacoustic tomography on graphics processing units.

PURPOSE Optoacoustic tomography (OAT) is inherently a three-dimensional (3D) inverse problem. However, most studies of OAT image reconstruction still employ two-dimensional imaging models. One important reason is because 3D image reconstruction is computationally burdensome. The aim of this work is to accelerate existing image reconstruction algorithms for 3D OAT by use of parallel programming techniques. METHODS Parallelization strategies are proposed to accelerate a filtered backprojection (FBP) algorithm and two different pairs of projection/backprojection operations that correspond to two different numerical imaging models. The algorithms are designed to fully exploit the parallel computing power of graphics processing units (GPUs). In order to evaluate the parallelization strategies for the projection/backprojection pairs, an iterative image reconstruction algorithm is implemented. Computer simulation and experimental studies are conducted to investigate the computational efficiency and numerical accuracy of the developed algorithms. RESULTS The GPU implementations improve the computational efficiency by factors of 1000, 125, and 250 for the FBP algorithm and the two pairs of projection/backprojection operators, respectively. Accurate images are reconstructed by use of the FBP and iterative image reconstruction algorithms from both computer-simulated and experimental data. CONCLUSIONS Parallelization strategies for 3D OAT image reconstruction are proposed for the first time. These GPU-based implementations significantly reduce the computational time for 3D image reconstruction, complementing our earlier work on 3D OAT iterative image reconstruction.

[1]  Stephen P. Boyd,et al.  Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers , 2011, Found. Trends Mach. Learn..

[2]  K. Morton,et al.  Numerical Solution of Partial Differential Equations: Introduction , 2005 .

[3]  Vasilis Ntziachristos,et al.  Fast Semi-Analytical Model-Based Acoustic Inversion for Quantitative Optoacoustic Tomography , 2010, IEEE Transactions on Medical Imaging.

[4]  Lihong V. Wang,et al.  Biomedical Optics: Principles and Imaging , 2007 .

[5]  Jeffrey A. Fessler.  Penalized weighted least-squares image reconstruction for positron emission tomography , 1994, IEEE Trans. Medical Imaging.

[6]  Otmar Scherzer,et al.  Reconstruction formulas for photoacoustic sectional imaging , 2011, 1109.0841.

[7]  Avinash C. Kak,et al.  Principles of computerized tomographic imaging , 2001, Classics in applied mathematics.

[8]  J. Shewchuk.  An Introduction to the Conjugate Gradient Method Without the Agonizing Pain , 1994 .

[9]  Lihong V. Wang,et al.  Photoacoustic imaging in biomedicine , 2006 .

[10]  Mark A. Anastasio,et al.  Photoacoustic and Thermoacoustic Tomography: Image Formation Principles , 2015, Handbook of Mathematical Methods in Imaging.

[11]  Lihong V Wang,et al.  Compressed sensing in photoacoustic tomography in vivo. , 2010, Journal of biomedical optics.

[12]  Mark A. Anastasio,et al.  Investigation of iterative image reconstruction in optoacoustic tomography , 2012, Photonics West - Biomedical Optics.

[13]  Fumihiko Ino,et al.  High-performance cone beam reconstruction using CUDA compatible GPUs , 2010, Parallel Comput..

[14]  B T Cox,et al.  k-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields. , 2010, Journal of biomedical optics.

[15]  Alexander A Karabutov,et al.  Optoacoustic imaging of absorbing objects in a turbid medium: ultimate sensitivity and application to breast cancer diagnostics. , 2007, Applied optics.

[16]  Richard Su,et al.  Whole-body three-dimensional optoacoustic tomography system for small animals. , 2009, Journal of biomedical optics.

[17]  Huabei Jiang,et al.  Three-dimensional finite-element-based photoacoustic tomography: reconstruction algorithm and simulations. , 2007, Medical physics.

[18]  Harrison H. Barrett,et al.  Foundations of Image Science , 2003, J. Electronic Imaging.

[19]  Simon R Arridge,et al.  Two-dimensional quantitative photoacoustic image reconstruction of absorption distributions in scattering media by use of a simple iterative method. , 2006, Applied optics.

[20]  Jin Zhang,et al.  Effects of Different Imaging Models on Least-Squares Image Reconstruction Accuracy in Photoacoustic Tomography , 2009, IEEE Transactions on Medical Imaging.

[21]  V. Ntziachristos,et al.  Model-based optoacoustic inversions with incomplete projection data. , 2011, Medical physics.

[22]  Vasilis Ntziachristos,et al.  Accurate Model-Based Reconstruction Algorithm for Three-Dimensional Optoacoustic Tomography , 2012, IEEE Transactions on Medical Imaging.

[23]  Chao Huang,et al.  Aberration correction for transcranial photoacoustic tomography of primates employing adjunct image data. , 2012, Journal of biomedical optics.

[24]  Xing Zhao,et al.  GPU-Based 3D Cone-Beam CT Image Reconstruction for Large Data Volume , 2009, Int. J. Biomed. Imaging.

[25]  Mark A. Anastasio,et al.  Investigation of limited-view image reconstruction in optoacoustic tomography employing a priori structural information , 2010, Optical Engineering + Applications.

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

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

[28]  Huabei Jiang,et al.  Photoacoustic tomography system for noninvasive real-time three-dimensional imaging of epilepsy , 2012, Biomedical optics express.

[29]  Gengsheng Lawrence Zeng,et al.  Unmatched projector/backprojector pairs in an iterative reconstruction algorithm , 2000, IEEE Transactions on Medical Imaging.

[30]  Mark A. Anastasio,et al.  Application of inverse source concepts to photoacoustic tomography , 2007 .

[31]  Weichung Wang,et al.  A fast forward projection using multithreads for multirays on GPUs in medical image reconstruction. , 2011, Medical physics.

[32]  Minghua Xu,et al.  Time-domain reconstruction for thermoacoustic tomography in a spherical geometry , 2002, IEEE Transactions on Medical Imaging.

[33]  M. Anastasio,et al.  Investigation of iterative image reconstruction in three-dimensional optoacoustic tomography , 2012, Physics in medicine and biology.

[34]  L. Kunyansky,et al.  Explicit inversion formulae for the spherical mean Radon transform , 2006, math/0609341.

[35]  Erik Lindholm,et al.  NVIDIA Tesla: A Unified Graphics and Computing Architecture , 2008, IEEE Micro.

[36]  Liang Song,et al.  Compressed-sensing photoacoustic computed tomography in vivo with partially known support , 2012 .

[37]  Quing Zhu,et al.  In vivo photoacoustic tomography of mouse cerebral edema induced by cold injury. , 2011, Journal of biomedical optics.

[38]  Pinhas Ephrat,et al.  Three-dimensional photoacoustic imaging by sparse-array detection and iterative image reconstruction. , 2008, Journal of biomedical optics.

[39]  V Ntziachristos,et al.  Optoacoustic methods for frequency calibration of ultrasonic sensors , 2011, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[40]  Jin Zhang,et al.  Feasibility of half-data image reconstruction in 3-D reflectivity tomography with a spherical aperture , 2005, IEEE Transactions on Medical Imaging.

[41]  V Ntziachristos,et al.  Three-dimensional optoacoustic tomography at video rate. , 2012, Optics express.

[42]  Markus Haltmeier,et al.  Inversion of Spherical Means and the Wave Equation in Even Dimensions , 2007, SIAM J. Appl. Math..

[43]  Lihong V. Wang,et al.  Universal back-projection algorithm for photoacoustic computed tomography. , 2005 .

[44]  Lihong V. Wang,et al.  Tutorial on Photoacoustic Microscopy and Computed Tomography , 2008, IEEE Journal of Selected Topics in Quantum Electronics.

[45]  Jeffrey A. Fessler,et al.  Conjugate-gradient preconditioning methods for shift-variant PET image reconstruction , 1999, IEEE Trans. Image Process..

[46]  R A Kruger,et al.  Thermoacoustic computed tomography--technical considerations. , 1999, Medical physics.

[47]  Makoto Yamakawa,et al.  Model-Based Reconstruction Integrated With Fluence Compensation for Photoacoustic Tomography , 2012, IEEE Transactions on Biomedical Engineering.

[48]  M. Wernick,et al.  Emission Tomography: The Fundamentals of PET and SPECT , 2004 .

[49]  Rakesh,et al.  Determining a Function from Its Mean Values Over a Family of Spheres , 2004, SIAM J. Math. Anal..

[50]  Mark A. Anastasio,et al.  An Imaging Model Incorporating Ultrasonic Transducer Properties for Three-Dimensional Optoacoustic Tomography , 2011, IEEE Transactions on Medical Imaging.

[51]  Lihong V. Wang,et al.  Photoacoustic tomography of water in phantoms and tissue. , 2010, Journal of biomedical optics.

[52]  Jin Zhang,et al.  Half-time image reconstruction in thermoacoustic tomography , 2005, IEEE Transactions on Medical Imaging.

[53]  Lihong V. Wang,et al.  Photoacoustic computed tomography correcting for heterogeneity and attenuation. , 2012, Journal of biomedical optics.

[54]  Yuan Xu,et al.  Exact frequency-domain reconstruction for thermoacoustic tomography. I. Planar geometry , 2002, IEEE Transactions on Medical Imaging.

[55]  Justin P. Haldar,et al.  Accelerating advanced MRI reconstructions on GPUs , 2008, J. Parallel Distributed Comput..

[56]  S. Jacques,et al.  Iterative reconstruction algorithm for optoacoustic imaging. , 2002, The Journal of the Acoustical Society of America.