Biomedical Photoacoustic Imaging Optimization with Deconvolution and EMD Reconstruction

A photoacoustic (PA) signal of an ideal optical absorbing particle is a single N-shape wave. PA signals are a combination of several individual N-shape waves. However, the N-shape wave basis leads to aliasing between adjacent micro-structures, which deteriorates the quality of final PA images. In this paper, we propose an image optimization method by processing raw PA signals with deconvolution and empirical mode decomposition (EMD). During the deconvolution procedure, the raw PA signals are de-convolved with a system dependent deconvolution kernel, which is measured in advance. EMD is subsequently adopted to further process the PA signals adaptively with two restrictive conditions: positive polarity and spectrum consistency. With this method, signal aliasing is alleviated, and the micro-structures and detail information, previously buried in the reconstructing images, can now be revealed. To validate our proposed method, numerical simulations and phantom studies are implemented, and reconstructed images are used for illustration.

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

[2]  Jie Yuan,et al.  Spread Spectrum Photoacoustic Tomography With Image Optimization , 2017, IEEE Transactions on Biomedical Circuits and Systems.

[3]  Lihong V. Wang,et al.  Noninvasive laser-induced photoacoustic tomography for structural and functional in vivo imaging of the brain , 2003, Nature Biotechnology.

[4]  N. Huang,et al.  The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis , 1998, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

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

[6]  A. Oraevsky,et al.  Laser optoacoustic imaging system for detection of breast cancer. , 2009, Journal of biomedical optics.

[7]  Quing Zhu,et al.  A real-time photoacoustic tomography system for small animals. , 2009, Optics express.

[8]  Joanna Brunker,et al.  Pulsed photoacoustic Doppler flowmetry using time-domain cross-correlation: accuracy, resolution and scalability. , 2012, The Journal of the Acoustical Society of America.

[9]  Minghua Xu,et al.  Thermoacoustic and Photoacoustic Tomography of Thick Biological Tissues toward Breast Imaging , 2005, Technology in cancer research & treatment.

[10]  Robert A Kruger,et al.  Thermoacoustic computed tomography using a conventional linear transducer array. , 2003, Medical physics.

[11]  De Cai,et al.  In vivo deconvolution acoustic-resolution photoacoustic microscopy in three dimensions. , 2016, Biomedical optics express.

[12]  Geng Ku,et al.  Noninvasive imaging of hemoglobin concentration and oxygenation in the rat brain using high-resolution photoacoustic tomography. , 2006, Journal of biomedical optics.

[13]  Lihong V Wang,et al.  Photoacoustic tomography and sensing in biomedicine , 2009, Physics in medicine and biology.

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

[15]  Xiaojun Liu,et al.  Adaptive optimization on ultrasonic transmission tomography-based temperature image for biomedical treatment , 2017 .

[16]  Y. Saijo,et al.  Basic study of improvement of axial resolution and suppression of time side lobe by phase-corrected Wiener filtering in photoacoustic tomography , 2018, Japanese Journal of Applied Physics.