Abstract:SignificancePhotoacoustic computed tomography (PACT) is a promising non-invasive imaging technique for both life science and clinical implementations. To achieve fast imaging speed, modern PACT systems have equipped arrays that have hundreds to thousands of ultrasound transducer (UST) elements, and the element number continues to increase. However, large number of UST elements with parallel data acquisition could generate a massive data size, making it very challenging to realize fast image reconstruction. Although several research groups have developed GPU-accelerated method for PACT, there lacks an explicit and feasible step-by-step description of GPU-based algorithms for various hardware platforms.AimIn this study, we propose a comprehensive framework for developing GPU-accelerated PACT image reconstruction (GPU-accelerated photoacoustic computed tomography), to help the research community to grasp this advanced image reconstruction method.ApproachWe leverage widely accessible open-source parallel computing tools, including Python multiprocessing-based parallelism, Taichi Lang for Python, CUDA, and possible other backends. We demonstrate that our framework promotes significant performance of PACT reconstruction, enabling faster analysis and real-time applications. Besides, we also described how to realize parallel computing on various hardware configurations, including multicore CPU, single GPU, and multiple GPUs platform.ResultsNotably, our framework can achieve an effective rate of ∼871 times when reconstructing extremely large-scale three-dimensional PACT images on a dual-GPU platform compared to a 24-core workstation CPU. In this paper, we share example codes via GitHub.ConclusionsOur approach allows for easy adoption and adaptation by the research community, fostering implementations of PACT for both life science and medicine.

Three-dimensional Photoacoustic Tomography Based on Graphics-Processing-unit-accelerated Finite Element Method.

3-D Ultrasonic Image Reconstruction in Frequency Domain Using a Virtual Transducer Model.

The Matrix Transform for Reconstruction on Finite-Element-Based Photoacoustic Tomography

Accelerating Image Reconstruction in Three-Dimensional Optoacoustic Tomography on Graphics Processing Units

Comprehensive framework of GPU-accelerated image reconstruction for photoacoustic computed tomography

GPU-Based Parallel Computing Methods for Medical Photoacoustic Image Reconstruction

A fast marching method based back projection algorithm for photoacoustic tomography in heterogeneous media

Pet Image Reconstruction: GPU-accelerated Particle Filter Framework

Finite-Element-Based Photoacoustic Tomography: Phantom And Chicken Bone Experiments

Fast and High-Resolution Three-Dimensional Hybrid-Domain Photoacoustic Imaging Incorporating Analytical-Focused Transducer Beam Amplitude

FPGA Acceleration of Image Reconstruction for Real-Time Photoacoustic Tomography

Graphics processing unit (GPU)-accelerated particle filter framework for positron emission tomography image reconstruction.

Accelerated 3-D T-FDK Reconstruction Algorithm Using the Commordity Graphics Hardware

Multigrid-based Reconstruction Algorithm for Quantitative Photoacoustic Tomography.

Reconstruction Algorithm in Photoacoustic Tomography

Finite-element-based photoacoustic tomography in time domain

Rapid Three-Dimensional Photoacoustic Imaging Reconstruction for Irregularly Layered Heterogeneous Media.

Parallel Computing For Quantitative Blood Flow Imaging In Photoacoustic Microscopy

Fluorescence Molecular Tomography Using a Two-Step Three-Dimensional Shape-Based Reconstruction with Graphics Processing Unit Acceleration.

The Cublas and Cula Based Gpu Acceleration of Adaptive Finite Element Framework for Bioluminescence Tomography

GPU-accelerated Two Dimensional Synthetic Aperture Focusing for Photoacoustic Microscopy