A Framework for Directional and Higher-Order Reconstruction in Photoacoustic Tomography

Abstract

Photoacoustic tomography is a hybrid imaging technique that combines high optical tissue contrast with high ultrasound resolution. Direct reconstruction methods such as filtered backprojection, time reversal and least squares suffer from curved line artefacts and blurring, especially in case of limited angles or strong noise. In recent years, there has been great interest in regularised iterative methods. These methods employ prior knowledge on the image to provide higher quality reconstructions. However, easy comparisons between regularisers and their properties are limited, since many tomography implementations heavily rely on the specific regulariser chosen. To overcome this bottleneck, we present a modular reconstruction framework for photoacoustic tomography. It enables easy comparisons between regularisers with different properties, e.g. nonlinear, higher-order or directional. We solve the underlying minimisation problem with an efficient first-order primal-dual algorithm. Convergence rates are optimised by choosing an operator dependent preconditioning strategy. Our reconstruction methods are tested on challenging 2D synthetic and experimental data sets. They outperform direct reconstruction approaches for strong noise levels and limited angle measurements, offering immediate benefits in terms of acquisition time and quality. This work provides a basic platform for the investigation of future advanced regularisation methods in photoacoustic tomography.

arXiv commentsubmitted to "Physics in Medicine and Biology". Changes from v1 to v2: regularisation with directional wavelet has been added; new experimental tests have been included

References in corpus (8)

• An inertial forward-backward algorithm for monotone inclusions
• Convex optimization problem prototyping for image reconstruction in computed tomography with the Chambolle-Pock algorithm
• Accelerated High-Resolution Photoacoustic Tomography via Compressed Sensing
• On the Adjoint Operator in Photoacoustic Tomography
• Initial results of finger imaging using Photoacoustic Computed Tomography
• Universal inversion formulas for recovering a function from spherical means
• Compressed sensing and sparsity in photoacoustic tomography
• Inversion of the spherical means transform in corner-like domains by reduction to the classical Radon transform

Cited by in corpus (7)

• Model based learning for accelerated, limited-view 3D photoacoustic tomography
• Deep Learning in Photoacoustic Tomography: Current approaches and future directions
• A Partially Learned Algorithm for Joint Photoacoustic Reconstruction and Segmentation
• Faster PET Reconstruction with Non-Smooth Priors by Randomization and Preconditioning
• Enhancing Compressed Sensing 4D Photoacoustic Tomography by Simultaneous Motion Estimation
• Approximate k-space models and Deep Learning for fast photoacoustic reconstruction
• Photo-acoustic tomographic image reconstruction from reduced data using physically inspired regularization