Signed Real-Time Delay Multiply and Sum Beamforming for Multispectral Photoacoustic Imaging

Reconstruction of photoacoustic (PA) images acquired with clinical ultrasound transducers is usually performed using the Delay and Sum (DAS) beamforming algorithm. Recently, a variant of DAS, referred to as Delay Multiply and Sum (DMAS) beamforming has been shown to provide increased contrast, signal-to-noise ratio (SNR) and resolution in PA imaging. The main reasons for the use of DAS beamforming in photoacoustics are its simple implementation, real-time capability, and the linearity of the beamformed image to the PA signal. This is crucial for the identification of different chromophores in multispectral PA applications. In contrast, current DMAS implementations are not responsive to the full spectrum of sound frequencies from a photoacoustic source and have not been shown to provide a reconstruction linear to the PA signal. Furthermore, due to its increased computational complexity, DMAS has not been shown yet to work in real-time. Here, we present an open-source real-time variant of the DMAS algorithm, signed DMAS (sDMAS), that ensures linearity in the original PA signal response while providing the increased image quality of DMAS. We show the applicability of sDMAS for multispectral PA applications, in vitro and in vivo. The sDMAS and reference DAS algorithms were integrated in the open-source Medical Imaging Interaction Toolkit (MITK) and are available as real-time capable implementations.

[1]  Cheri X Deng,et al.  Photoacoustic spectrum analysis for microstructure characterization in biological tissue: A feasibility study. , 2012, Applied physics letters.

[2]  Liang Song,et al.  Delay-multiply-and-sum-based synthetic aperture focusing in photoacoustic microscopy , 2016, Journal of biomedical optics.

[3]  Benjamin A. Rockwell,et al.  A procedure for multiple-pulse maximum permissible exposure determination under the Z136.1-2000 American National Standard for Safe Use of Lasers , 2001 .

[4]  F. Gran,et al.  Broadband minimum variance beamforming for ultrasound imaging , 2009, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[5]  Ali Mahloojifar,et al.  Double-Stage Delay Multiply and Sum Beamforming Algorithm: Application to Linear-Array Photoacoustic Imaging , 2018, IEEE Transactions on Biomedical Engineering.

[6]  S L Jacques,et al.  Optical properties of intralipid: A phantom medium for light propagation studies , 1992, Lasers in surgery and medicine.

[7]  Charles L. Lawson,et al.  Solving least squares problems , 1976, Classics in applied mathematics.

[8]  F. Harris On the use of windows for harmonic analysis with the discrete Fourier transform , 1978, Proceedings of the IEEE.

[9]  Ali Mahloojifar,et al.  Linear-array photoacoustic imaging using minimum variance-based delay multiply and sum adaptive beamforming algorithm , 2018, Journal of biomedical optics.

[10]  S. Jacques Corrigendum: Optical properties of biological tissues: a review , 2013 .

[11]  B. Cox,et al.  Photoacoustic tomography in absorbing acoustic media using time reversal , 2010 .

[12]  Klaus H. Maier-Hein,et al.  The Medical Imaging Interaction Toolkit: challenges and advances , 2013, International Journal of Computer Assisted Radiology and Surgery.

[13]  Piero Tortoli,et al.  High Frame-Rate, High Resolution Ultrasound Imaging With Multi-Line Transmission and Filtered-Delay Multiply And Sum Beamforming , 2017, IEEE Trans. Medical Imaging.

[14]  L. J. Griffiths,et al.  An alternative approach to linearly constrained adaptive beamforming , 1982 .

[15]  S. Jacques Optical properties of biological tissues: a review , 2013, Physics in medicine and biology.

[16]  Chulhong Kim,et al.  Programmable Real-time Clinical Photoacoustic and Ultrasound Imaging System , 2016, Scientific Reports.

[17]  Stanislav Y. Emelianov,et al.  Optical wavelength selection for improved spectroscopic photoacoustic imaging☆ , 2013, Photoacoustics.

[18]  R Zander,et al.  The oxygen status of arterial human blood. , 1990, Scandinavian journal of clinical and laboratory investigation. Supplementum.

[19]  Yves Rosseel,et al.  On the Definition of Signal-To-Noise Ratio and Contrast-To-Noise Ratio for fMRI Data , 2013, PloS one.

[20]  Sevan Harput,et al.  Spatial resolution and contrast enhancement in photoacoustic imaging with filter delay multiply and sum beamforming technique , 2016, 2016 IEEE International Ultrasonics Symposium (IUS).

[21]  Lena Maier-Hein,et al.  Context encoding enables machine learning-based quantitative photoacoustics , 2017, Journal of biomedical optics.

[22]  Giovanni Magenes,et al.  The Delay Multiply and Sum Beamforming Algorithm in Ultrasound B-Mode Medical Imaging , 2015, IEEE Transactions on Medical Imaging.

[23]  Kristofer Schlachter,et al.  An Introduction to the OpenCL Programming Model , 2012 .