GPU-accelerated Double-stage Delay-multiply-and-sum Algorithm for Fast Photoacoustic Tomography Using LED Excitation and Linear Arrays

Double-stage delay-multiply-and-sum (DS-DMAS) is an algorithm proposed for photoacoustic image reconstruction. The DS-DMAS algorithm offers a higher contrast than conventional delay-and-sum and delay-multiply and-sum but at the expense of higher computational complexity. Here, we utilized a compute unified device architecture (CUDA) graphics processing unit (GPU) parallel computation approach to address the high complexity of the DS-DMAS for photoacoustic image reconstruction generated from a commercial light-emitting diode (LED)-based photoacoustic scanner. In comparison with a single-threaded central processing unit (CPU), the GPU approach increased speeds by nearly 140-fold for 1024 × 1024 pixel image; there was no decrease in accuracy. The proposed implementation makes it possible to reconstruct photoacoustic images with frame rates of 250, 125, and 83.3 when the images are 64 × 64, 128 × 128, and 256 × 256, respectively. Thus, DS-DMAS can be efficiently used in clinical devices when coupled with CUDA GPU parallel computation.

Acoustic dipole radiation based electrical impedance contrast imaging approach of magnetoacoustic tomography with magnetic induction

PURPOSE

Different from the theory of acoustic monopole spherical radiation, the acoustic dipole radiation based theory introduces the radiation pattern of Lorentz force induced dipole sources to describe the principle of magnetoacoustic tomography with magnetic induction (MAT-MI). Although two-dimensional (2D) simulations have been studied for cylindrical phantom models, layer effects of the dipole sources within the entire object along the z direction still need to be investigated to evaluate the performance of MAT-MI for different geometric specifications. The purpose of this work is further verifying the validity and generality of acoustic dipole radiation based theory for MAT-MI with two new models in different shapes, dimensions, and conductivities.

METHODS

Based on the theory of acoustic dipole radiation, the principles of MAT-MI were analyzed with derived analytic formulae. 2D and 3D numerical studies for two new models of aluminum foil and cooked egg were conducted to simulate acoustic pressures and corresponding waveforms, and 2D images of the scanned layers were reconstructed with the simplified back projection algorithm for the waveforms collected around the models. The spatial resolution for conductivity boundary differentiation was also analyzed with different foil thickness. For comparison, two experimental measurements were conducted for a cylindrical aluminum foil phantom and a shell-peeled cooked egg. The collected waveforms and the reconstructed images of the scanned layers were achieved to verify the validation of the acoustic dipole radiation based theory for MAT-MI.

RESULTS

Despite the difference between the 2D and 3D simulated pressures, good consistence of the collected waveforms proves that wave clusters are generated by the abrupt pressure changes with bipolar vibration phases, representing the opposite polarities of the conductivity changes along the measurement direction. The configuration of the scanned layer can be reconstructed in terms of shape and size, and the conductivity boundaries are displayed in stripes with different contrast and bipolar intensities. Layer effects are demonstrated to have little influence on the collected waveforms and the reconstructed images of the scanned layers for the two new models. The experimental results have good agreements with numerical simulations, and the reconstructed 2D images provide conductivity configurations in the scanned layers of the aluminum foil and the egg models.

CONCLUSIONS

It can be concluded that the acoustic pressure of MAT-MI is produced by the divergence of the induced Lorentz force, and the collected waveforms comprise wave clusters with bipolar vibration phases and different amplitudes, providing the information of conductivity boundaries in the scanned layer. With the simplified back projection algorithm for diffraction sources, collected waveforms can be used to reconstruct 2D conductivity contrast image and the conductivity configuration in the scanned layer can be obtained in terms of shape and size in stripes with the spatial resolution of the acoustic wavelength. The favorable results further verify the validity and generality of the acoustic dipole radiation based theory and suggest the feasibility of MAT-MI as an effective electrical impedance contrast imaging approach for medical imaging.

Biomedical photoacoustics in China

During the last decade, along with its explosive growth globally, biomedical photoacoustics has become a rapidly growing research field in China as well. In particular, photoacoustic tomography (PAT), capable of imaging intact biological tissue in vivo at great depths, has generated intense interest among Chinese researchers. This review briefly summarizes the current status and recent progress of the research in PAT in China. The focus is on the technology development and biomedical applications of three representative embodiments of PAT: photoacoustic microscopy, photoacoustic computed tomography, and photoacoustic endoscopy. In addition, recent development and studies in other related areas are also reviewed shortly.

4-D photoacoustic tomography

Photoacoustic tomography (PAT) offers three-dimensional (3D) structural and functional imaging of living biological tissue with label-free, optical absorption contrast. These attributes lend PAT imaging to a wide variety of applications in clinical medicine and preclinical research. Despite advances in live animal imaging with PAT, there is still a need for 3D imaging at centimeter depths in real-time. We report the development of four dimensional (4D) PAT, which integrates time resolutions with 3D spatial resolution, obtained using spherical arrays of ultrasonic detectors. The 4D PAT technique generates motion pictures of imaged tissue, enabling real time tracking of dynamic physiological and pathological processes at hundred micrometer-millisecond resolutions. The 4D PAT technique is used here to image needle-based drug delivery and pharmacokinetics. We also use this technique to monitor 1) fast hemodynamic changes during inter-ictal epileptic seizures and 2) temperature variations during tumor thermal therapy.