A Reconstruction Algorithm for Photoacoustic Imaging based on the Nonuniform FFT

Abstract

Fourier reconstruction algorithms significantly outperform conventional back-projection algorithms in terms of computation time. In photoacoustic imaging, these methods require interpolation in the Fourier space domain, which creates artifacts in reconstructed images. We propose a novel reconstruction algorithm that applies the one-dimensional nonuniform fast Fourier transform to photoacoustic imaging. It is shown theoretically and numerically that our algorithm avoids artifacts while preserving the computational effectiveness of Fourier reconstruction.

22 pages, 8 figures

References in corpus (1)

• A series solution and a fast algorithm for the inversion of the spherical mean Radon transform

Cited by in corpus (11)

• Fully Dense UNet for 2D Sparse Photoacoustic Tomography Artifact Removal
• Inversion of circular means and the wave equation on convex planar domains
• On Artifacts in Limited Data Spherical Radon Transform: Curved Observation Surface
• Exact Reconstruction Formula for the Spherical Mean Radon Transform on Ellipsoids
• Non-Equispaced Grid Sampling in Photoacoustics with a Non-Uniform FFT
• Artifacts in incomplete data tomography - with applications to photoacoustic tomography and sonar
• Application of a nudging technique to thermoacoustic tomography
• Efficient regularization with wavelet sparsity constraints in PAT
• The spherical mean Radon transform with centers on cylindrical surfaces
• Fast reconstruction algorithms for the thermoacoustic tomography in certain domains with cylindrical or spherical symmetries
• Computational Realization of a Non-Equidistant Grid Sampling in Photoacoustics with a Non-Uniform FFT