Abstract:Purpose: The purpose of this study is to develop an improved filtered-back-projection (FBP) algorithm for photoacoustic tomography (PAT), which allows image reconstruction with higher quality compared to images reconstructed through traditional algorithms. Methods: A rigorous expression of a weighting function has been derived directly from a photoacoustic wave equation and used as a ramp filter in Fourier domain. The authors' new algorithm utilizes this weighting function to precisely calculate each photoacoustic signal's contribution and then reconstructs the image based on the retarded potential generated from the photoacoustic sources. In addition, an adaptive criterion has been derived for selecting the cutoff frequency of a low pass filter. Two computational phantoms were created to test the algorithm. The first phantom contained five spheres with each sphere having different absorbances. The phantom was used to test the capability for correctly representing both the geometry and the relative absorbed energy in a planar measurement system. The authors also used another phantom containing absorbers of different sizes with overlapping geometry to evaluate the performance of the new method for complicated geometry. In addition, random noise background was added to the simulated data, which were obtained by using an arc-shaped array of 50 evenly distributed transducers that spanned 160° over a circle with a radius of 65 mm. A normalized factor between the neighbored transducers was applied for correcting measurement signals in PAT simulations. The authors assumed that the scanned object was mounted on a holder that rotated over the full 360° and the scans were set to a sampling rate of 20.48 MHz. Results: The authors have obtained reconstructed images of the computerized phantoms by utilizing the new FBP algorithm. From the reconstructed image of the first phantom, one can see that this new approach allows not only obtaining a sharp image but also showing the correct signal strength of the absorbers. The reconstructed image of the second phantom further demonstrates the capability to form clear images of the spheres with sharp borders in the overlapping geometry. The smallest sphere is clearly visible and distinguishable, even though it is surrounded by two big spheres. In addition, image reconstructions were conducted with randomized noise added to the observed signals to mimic realistic experimental conditions. Conclusions: The authors have developed a new FBP algorithm that is capable for reconstructing high quality images with correct relative intensities and sharp borders for PAT. The results demonstrate that the weighting function serves as a precise ramp filter for processing the observed signals in the Fourier domain. In addition, this algorithm allows an adaptive determination of the cutoff frequency for the applied low pass filter.

Implementation of Medical Image Reconstruction Algorithm for Photoacoustic Imaging Using k-wave Toolkit

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

S-Wave Accelerates Optimization-based Photoacoustic Image Reconstruction in vivo

Single-stage reconstruction algorithm for quantitative photoacoustic tomography

Real-time reconstruction of the medical complex frequency domain photoacoustic image

Effect of Iterative Reconstruction Integrating Sart and Fbp on Photoacoustic Imaging

A Reconstruction Algorithm for Photoacoustic Imaging based on the Nonuniform FFT

Photoacoustic computed tomography for joint reconstruction of initial pressure and sound speed in vivo using a feature coupling method

Multigrid-based Reconstruction Algorithm for Quantitative Photoacoustic Tomography.

Efficient Discrete Cosine Transform Model–based Algorithm for Photoacoustic Image Reconstruction

Biomedical Photoacoustic Imaging with Unknown Spatially Distributed Ultrasound Sensor Array

Feature Coupling Photoacoustic Computed Tomography for Joint Reconstruction of Initial Pressure and Sound Speed in Vivo

An adaptive filtered back-projection for photoacoustic image reconstruction

Research on the photoacoustic spectrum analysis using k-Wave

Simultaneous Algebraic Reconstruction Technique Based on Total Variation for Photoacoustic Image

A high contrast and resolution reconstruction algorithm in quantitative photoacoustic tomography

Total Variation-Based Limited-View Photoacoustic Reconstruction Method with Segmentation-Based Regularization

Conjugate Gradient Preconditioning Methods with Symmetric Algebraic Reconstruction Technique in Photoacoustic Imaging

Multiple Delay and Sum with Enveloping Beamforming Algorithm for Photoacoustic Imaging

An Iterative Reconstruction Method for Photoacoustic Tomography in Heterogeneous Media Based on Point Source Response

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