Efficient regularization with wavelet sparsity constraints in PAT

Abstract

In this paper we consider the reconstruction problem of photoacoustic tomography (PAT) with a flat observation surface. We develop a direct reconstruction method that employs regularization with wavelet sparsity constraints. To that end, we derive a wavelet-vaguelette decomposition (WVD) for the PAT forward operator and a corresponding explicit reconstruction formula in the case of exact data. In the case of noisy data, we combine the WVD reconstruction formula with soft-thresholding which yields a spatially adaptive estimation method. We demonstrate that our method is statistically optimal for white random noise if the unknown function is assumed to lie in any Besov-ball. We present generalizations of this approach and, in particular, we discuss the combination of vaguelette soft-thresholding with a TV prior. We also provide an efficient implementation of the vaguelette transform that leads to fast image reconstruction algorithms supported by numerical results.

arXiv comment25 pages, 6 figures

References in corpus (11)

• Simultaneous analysis of Lasso and Dantzig selector
• Investigation of iterative image reconstruction in three-dimensional optoacoustic tomography
• Sparse Regularization with $l^q$ Penalty Term
• Convergence rates and source conditions for Tikhonov regularization with sparsity constraints
• Compressed sensing and sparsity in photoacoustic tomography
• Inversion of circular means and the wave equation on convex planar domains
• A Reconstruction Algorithm for Photoacoustic Imaging based on the Nonuniform FFT
• Discrete Imaging Models for Three-Dimensional Optoacoustic Tomography using Radially Symmetric Expansion Functions
• Mitigation of artifacts due to isolated acoustic heterogeneities in photoacoustic computed tomography using a variable data truncation-based reconstruction method
• On Artifacts in Limited Data Spherical Radon Transform: Curved Observation Surface
• Inversion of the Spherical Mean Transform with Sources on a Hyperplane

Cited by in corpus (1)

• Fully Dense UNet for 2D Sparse Photoacoustic Tomography Artifact Removal