A Constrained Variable Projection Reconstruction Method for Photoacoustic Computed Tomography Without Accurate Knowledge of Transducer Responses

Photoacoustic computed tomography (PACT) is an emerging computed imaging modality that exploits optical contrast and ultrasonic detection principles to form images of the absorbed optical energy density within tissue. When the imaging system employs conventional piezoelectric ultrasonic transducers, the ideal photoacoustic (PA) signals are degraded by the transducers' acousto-electric impulse responses (EIRs) during the measurement process. If unaccounted for, this can degrade the accuracy of the reconstructed image. In principle, the effect of the EIRs on the measured PA signals can be ameliorated via deconvolution; images can be reconstructed subsequently by application of a reconstruction method that assumes an idealized EIR. Alternatively, the effect of the EIR can be incorporated into an imaging model and implicitly compensated for during reconstruction. In either case, the efficacy of the correction can be limited by errors in the assumed EIRs. In this work, a joint optimization approach to PACT image reconstruction is proposed for mitigating errors in reconstructed images that are caused by use of an inaccurate EIR. The method exploits the bi-linear nature of the imaging model and seeks to refine the measured EIR during the process of reconstructing the sought-after absorbed optical energy density. Computer-simulation and experimental studies are conducted to investigate the numerical properties of the method and demonstrate its value for mitigating image distortions and enhancing the visibility of fine structures.

[1]  Xu Xiao Photoacoustic imaging in biomedicine , 2008 .

[2]  Gerald R. Harris,et al.  Review of transient field theory for a baffled planar piston , 1981 .

[3]  P. Stepanishen Transient Radiation from Pistons in an Infinite Planar Baffle , 1970 .

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

[5]  James V. Burke,et al.  Algorithmic Differentiation of Implicit Functions and Optimal Values , 2008 .

[6]  Lihong V. Wang,et al.  Enhancement of photoacoustic tomography by ultrasonic computed tomography based on optical excitation of elements of a full-ring transducer array. , 2013, Optics letters.

[7]  Jinyi Qi,et al.  Effect of errors in the system matrix on maximum a posteriori image reconstruction , 2005, Physics in medicine and biology.

[8]  R. Krimholtz,et al.  Equivalent Circuits for Transducers Having Arbitrary Even- Or Odd-Symmetry Piezoelectric Excitation , 1971, IEEE Transactions on Sonics and Ultrasonics.

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

[10]  Ivan Pelivanov,et al.  Focused array transducer for two-dimensional optoacoustic tomography. , 2004, The Journal of the Acoustical Society of America.

[11]  R. Dhillon,et al.  For the safe use of lasers , 1989 .

[12]  Vasilis Ntziachristos,et al.  Model-based optoacoustic inversion with arbitrary-shape detectors. , 2011, Medical physics.

[13]  D. Razansky,et al.  Optoacoustic determination of spatio- temporal responses of ultrasound sensors , 2013, IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control.

[14]  Mark A. Anastasio,et al.  Discrete Imaging Models for Three-Dimensional Optoacoustic Tomography Using Radially Symmetric Expansion Functions , 2013, IEEE Transactions on Medical Imaging.

[15]  V. Ntziachristos,et al.  Acoustic Inversion in Optoacoustic Tomography: A Review , 2013, Current medical imaging reviews.

[16]  A. Gilbert,et al.  A generalization of variable elimination for separable inverse problems beyond least squares , 2013, 1302.0441.

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

[18]  Aleksandr Y. Aravkin,et al.  Estimating nuisance parameters in inverse problems , 2012, 1206.6532.

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

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

[21]  Lihong V. Wang,et al.  Whole-body ring-shaped confocal photoacoustic computed tomography of small animals in vivo. , 2012, Journal of biomedical optics.

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

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

[24]  Lihong V. Wang Photoacoustic imaging and spectroscopy , 2009 .

[25]  Dianne P. O'Leary,et al.  The Use of the L-Curve in the Regularization of Discrete Ill-Posed Problems , 1993, SIAM J. Sci. Comput..

[26]  G. Golub,et al.  Separable nonlinear least squares: the variable projection method and its applications , 2003 .

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

[28]  Richard Su,et al.  Measurement of the spectral directivity of optoacoustic and ultrasonic transducers with a laser ultrasonic source. , 2009, The Review of scientific instruments.

[29]  Quing Zhu,et al.  Curved array photoacoustic tomographic system for small animal imaging. , 2008, Journal of biomedical optics.

[30]  Chao Huang,et al.  Accelerating image reconstruction in three-dimensional optoacoustic tomography on graphics processing units. , 2013, Medical physics.

[31]  Bradley E Treeby,et al.  Automatic sound speed selection in photoacoustic image reconstruction using an autofocus approach. , 2011, Journal of biomedical optics.

[32]  Mark A. Anastasio,et al.  Investigation of the far-field approximation for modeling a transducer's spatial impulse response in photoacoustic computed tomography , 2013, Photoacoustics.

[33]  Dianne P. O'Leary,et al.  Variable projection for nonlinear least squares problems , 2012, Computational Optimization and Applications.

[34]  James G. Nagy,et al.  An Efficient Iterative Approach for Large-Scale Separable Nonlinear Inverse Problems , 2009, SIAM J. Sci. Comput..

[35]  S. Nash,et al.  Linear and Nonlinear Optimization , 2008 .

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

[37]  Gene H. Golub,et al.  The differentiation of pseudo-inverses and non-linear least squares problems whose variables separate , 1972, Milestones in Matrix Computation.

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

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