Basis pursuit deconvolution for improving model-based reconstructed images in photoacoustic tomography.

The model-based image reconstruction approaches in photoacoustic tomography have a distinct advantage compared to traditional analytical methods for cases where limited data is available. These methods typically deploy Tikhonov based regularization scheme to reconstruct the initial pressure from the boundary acoustic data. The model-resolution for these cases represents the blur induced by the regularization scheme. A method that utilizes this blurring model and performs the basis pursuit deconvolution to improve the quantitative accuracy of the reconstructed photoacoustic image is proposed and shown to be superior compared to other traditional methods via three numerical experiments. Moreover, this deconvolution including the building of an approximate blur matrix is achieved via the Lanczos bidagonalization (least-squares QR) making this approach attractive in real-time.

[1]  E.J. Candes,et al.  An Introduction To Compressive Sampling , 2008, IEEE Signal Processing Magazine.

[2]  Zhuang Liu,et al.  Carbon nanotubes as photoacoustic molecular imaging agents in living mice. , 2008, Nature nanotechnology.

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

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

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

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

[7]  Vasilis Ntziachristos,et al.  Acceleration of Optoacoustic Model-Based Reconstruction Using Angular Image Discretization , 2012, IEEE Transactions on Medical Imaging.

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

[9]  Jaya Prakash,et al.  A LSQR-type method provides a computationally efficient automated optimal choice of regularization parameter in diffuse optical tomography. , 2013, Medical physics.

[10]  Manojit Pramanik,et al.  Design and evaluation of a novel breast cancer detection system combining both thermoacoustic (TA) and photoacoustic (PA) tomography. , 2008, Medical physics.

[11]  Phaneendra K. Yalavarthy,et al.  Sparse Recovery Methods Hold Promise for Diffuse Optical Tomographic Image Reconstruction , 2014, IEEE Journal of Selected Topics in Quantum Electronics.

[12]  I. Selesnick Introduction to Sparsity in Signal Processing * , 2012 .

[13]  B. R. Hunt,et al.  The Application of Constrained Least Squares Estimation to Image Restoration by Digital Computer , 1973, IEEE Transactions on Computers.

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

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

[16]  Peter Kuchment,et al.  Mathematics of thermoacoustic and photoacoustic tomography , 2007 .

[17]  Liang Song,et al.  Blind-deconvolution optical-resolution photoacoustic microscopy in vivo. , 2013, Optics express.

[18]  Per Christian Hansen,et al.  Regularization Tools version 4.0 for Matlab 7.3 , 2007, Numerical Algorithms.

[19]  Manojit Pramanik,et al.  Deconvolution-based deblurring of reconstructed images in photoacoustic/thermoacoustic tomography. , 2013, Journal of the Optical Society of America. A, Optics, image science, and vision.

[20]  Manojit Pramanik,et al.  Least squares QR-based decomposition provides an efficient way of computing optimal regularization parameter in photoacoustic tomography , 2013, Journal of biomedical optics.

[21]  Lihong V. Wang,et al.  Noninvasive photoacoustic imaging of the thoracic cavity and the kidney in small and large animals. , 2008, Medical physics.

[22]  Brian W. Pogue,et al.  Model-Resolution-Based Basis Pursuit Deconvolution Improves Diffuse Optical Tomographic Imaging , 2014, IEEE Transactions on Medical Imaging.

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

[24]  K. P. Köstli,et al.  Two-dimensional photoacoustic imaging by use of Fourier-transform image reconstruction and a detector with an anisotropic response. , 2003, Applied optics.

[25]  Mark A. Anastasio,et al.  Full-Wave Iterative Image Reconstruction in Photoacoustic Tomography With Acoustically Inhomogeneous Media , 2013, IEEE Transactions on Medical Imaging.

[26]  Lihong V. Wang,et al.  Photoacoustic tomography of a nanoshell contrast agent in the in vivo rat brain , 2004 .

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

[28]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[29]  Alexander A. Oraevsky,et al.  Optoacoustic tomography: New modality of laser diagnostic systems , 2003 .

[30]  Misha Elena Kilmer,et al.  Choosing Regularization Parameters in Iterative Methods for Ill-Posed Problems , 2000, SIAM J. Matrix Anal. Appl..

[31]  J. Romberg,et al.  Imaging via Compressive Sampling , 2008, IEEE Signal Processing Magazine.

[32]  Jeffrey A. Fessler,et al.  Spatial resolution properties of penalized-likelihood image reconstruction: space-invariant tomographs , 1996, IEEE Trans. Image Process..

[33]  José M. Bioucas-Dias,et al.  Fast Image Recovery Using Variable Splitting and Constrained Optimization , 2009, IEEE Transactions on Image Processing.

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

[35]  Mário A. T. Figueiredo,et al.  Fast frame-based image deconvolution using variable splitting and constrained optimization , 2009, 2009 IEEE/SP 15th Workshop on Statistical Signal Processing.

[36]  Da Xing,et al.  Fast photoacoustic imaging system based on 320-element linear transducer array. , 2004, Physics in medicine and biology.

[37]  Lihong V. Wang,et al.  Functional photoacoustic microscopy for high-resolution and noninvasive in vivo imaging , 2006, Nature Biotechnology.

[38]  Jianhong Shen,et al.  Deblurring images: Matrices, spectra, and filtering , 2007, Math. Comput..

[39]  Linh V. Nguyen,et al.  Reconstruction and time reversal in thermoacoustic tomography in acoustically homogeneous and inhomogeneous media , 2008 .

[40]  Lihong V. Wang,et al.  Thermoacoustic and photoacoustic sensing of temperature. , 2009, Journal of biomedical optics.

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

[42]  L V Wang,et al.  Scanning thermoacoustic tomography in biological tissue. , 2000, Medical physics.