Disjoint sparsity for signal separation and applications to hybrid inverse problems in medical imaging

Abstract

The main focus of this work is the reconstruction of the signals $f$ and $g_i$, $i=1,...,N$, from the knowledge of their sums $h_i=f+g_i$, under the assumption that $f$ and the $g_i$'s can be sparsely represented with respect to two different dictionaries $A_f$ and $A_g$. This generalizes the well-known "morphological component analysis" to a multi-measurement setting. The main result of the paper states that $f$ and the $g_i$'s can be uniquely and stably reconstructed by finding sparse representations of $h_i$ for every $i$ with respect to the concatenated dictionary $[A_f,A_g]$, provided that enough incoherent measurements $g_i$ are available. The incoherence is measured in terms of their mutual disjoint sparsity. This method finds applications in the reconstruction procedures of several hybrid imaging inverse problems, where internal data are measured. These measurements usually consist of the main unknown multiplied by other unknown quantities, and so the disjoint sparsity approach can be directly applied. As an example, we show how to apply the method to the reconstruction in quantitative photoacoustic tomography, also in the case when the Grüneisen parameter, the optical absorption and the diffusion coefficient are all unknown.

33 pages, 6 figures

References in corpus (2)

• Single-stage reconstruction algorithm for quantitative photoacoustic tomography
• A one-step reconstruction algorithm for quantitative photoacoustic imaging

Cited by in corpus (6)

• Stochastic Proximal Gradient Algorithms for Multi-Source Quantitative Photoacoustic Tomography
• Reconstruction of a compactly supported sound profile in the presence of a random background medium
• Linearized Reconstruction for Diffuse Optical Spectroscopic Imaging
• A signal separation technique for sub-cellular imaging using dynamic optical coherence tomography
• Non-zero constraints in elliptic PDE with random boundary values and applications to hybrid inverse problems
• Efficient inclusion of total variation type priors in quantitative photoacoustic tomography