Deconvolution algorithms for photoacoustic tomography to reduce blurring caused by finite sized detectors

Most reconstruction algorithms for photoacoustic tomography, like back-projection or time-reversal, work ideally for point-like detectors. For real detectors, which integrate the pressure over their finite size, it was shown that images reconstructed by back-projection or time-reversal show some blurring. Iterative reconstruction algorithms using an imaging matrix can take the finite size of real detectors directly into account, but the numerical effort is significantly higher compared to the use of direct algorithms. For spherical or cylindrical detection surfaces the blurring caused by a finite detector size is proportional to the distance from the rotation center (“spin blur”) and is equal to the detector size at the detection surface. In this work we use deconvolution algorithms to reduce this type of blurring on simulated and on experimental data. Experimental data were obtained on a plastisol cylinder with 6 thin holes filled with an absorbing liquid (OrangeG). The holes were located on a spiral emanating from the center of the cylinder. Data acquisition was done by utilization of a piezoelectric detector which was rotated around the plastisol cylinder.

[1]  Re Gonzalez,et al.  R.C. Eddins, Digital image processing using MATLAB, vol. Gatesmark Publishing Knoxville , 2009 .

[2]  Gene H. Golub,et al.  Generalized cross-validation as a method for choosing a good ridge parameter , 1979, Milestones in Matrix Computation.

[3]  S. Ribaric,et al.  Restoration of images blurred by circular motion , 2000, IWISPA 2000. Proceedings of the First International Workshop on Image and Signal Processing and Analysis. in conjunction with 22nd International Conference on Information Technology Interfaces. (IEEE.

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

[5]  Lihong V. Wang,et al.  Time reversal and its application to tomography with diffracting sources. , 2004, Physical review letters.

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

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

[8]  Leonid Kunyansky,et al.  Fast reconstruction algorithms for the thermoacoustic tomography in certain domains with cylindrical or spherical symmetries , 2011, 1102.1413.

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

[10]  A. R. Mitchell,et al.  The Finite Difference Method in Partial Differential Equations , 1980 .

[11]  G. M.,et al.  Partial Differential Equations I , 2023, Applied Mathematical Sciences.

[12]  Otmar Scherzer,et al.  THERMOACOUSTIC TOMOGRAPHY AND THE CIRCULAR RADON TRANSFORM: EXACT INVERSION FORMULA , 2007 .

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

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

[15]  Minghua Xu,et al.  Analytic explanation of spatial resolution related to bandwidth and detector aperture size in thermoacoustic or photoacoustic reconstruction. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[16]  Markus Haltmeier,et al.  Spatial resolution in photoacoustic tomography: effects of detector size and detector bandwidth , 2010 .

[17]  Hao F. Zhang Functional Photoacoustic Microscopy , 2009 .

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

[19]  L. Trefethen Spectral Methods in MATLAB , 2000 .

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