Three-dimensional optoacoustic reconstruction using fast sparse representation

Optoacoustic tomography based on insufficient spatial sampling of ultrasound waves leads to loss of contrast and artifacts on the reconstructed images. Compared to reconstructions based on L2-norm regularization, sparsity-based reconstructions may improve contrast and reduce image artifacts but at a high computational cost, which has so far limited their use to 2D optoacoustic tomography. Here we propose a fast, sparsity-based reconstruction algorithm for 3D optoacoustic tomography, based on gradient descent with Barzilai-Borwein line search (L1-GDBB). Using simulations and experiments, we show that the L1-GDBB offers fourfold faster reconstruction than the previously reported L1-norm regularized reconstruction based on gradient descent with backtracking line search. Moreover, the new algorithm provides higher-quality images with fewer artifacts than the L2-norm regularized reconstruction and the back-projection reconstruction.