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.

A Reconstruction Algorithm for Photoacoustic Imaging based on the Nonuniform FFT

Single-stage reconstruction algorithm for quantitative photoacoustic tomography

Multigrid-based Reconstruction Algorithm for Quantitative Photoacoustic Tomography.

A high contrast and resolution reconstruction algorithm in quantitative photoacoustic tomography

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

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

A simple Fourier transform-based reconstruction formula for photoacoustic computed tomography with a circular or spherical measurement geometry

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

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

Effect of Iterative Reconstruction Integrating Sart and Fbp on Photoacoustic Imaging

Fast and Robust Deconvolution-Based Image Reconstruction for Photoacoustic Tomography in Circular Geometry: Experimental Validation

Efficient Block-Sparse Model-Based Algorithm for Photoacoustic Image Reconstruction

Photoacoustic image reconstruction: material detection and acoustical heterogeneities

Photoacoustic Tomography from Weak and Noisy Signals by Using a Pulse Decomposition Algorithm in the Time-Domain.

Model-based reconstructions for quantitative imaging in photoacoustic tomography

Learned backprojection for sparse and limited view photoacoustic tomography

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

An adaptive filtered back-projection for photoacoustic image reconstruction

Compressed sensing and sparsity in photoacoustic tomography

Perspectives In Fourier-Based Image Reconstruction In Photoacoustic Tomography

Research on ADMM Reconstruction Algorithm of Photoacoustic Tomography With Limited Sampling Data