Abstract
In this paper, we consider 3D Bioluminescence tomography (BLT) source reconstruction from Poisson data in three dimensional space. With a priori information of sources sparsity and MAP estimation of Poisson distribution, we study the minimization of Kullback-Leihbler divergence with ℓ1 and ℓ0 regularization. We show numerically that although several ℓ1 minimization algorithms are efficient for compressive sensing, they fail for BLT reconstruction due to the high coherence of the measurement matrix columns and high nonlinearity of Poisson fitting term. Instead, we propose a novel greedy algorithm for ℓ0 regularization to reconstruct sparse solutions for BLT problem. Numerical experiments on synthetic data obtained by the finite element methods and Monte-Carlo methods show the accuracy and efficiency of the proposed method.
Original language | English (US) |
---|---|
Pages (from-to) | 519-535 |
Number of pages | 17 |
Journal | Journal of Scientific Computing |
Volume | 50 |
Issue number | 3 |
DOIs | |
State | Published - Mar 1 2012 |
Externally published | Yes |
Bibliographical note
Publisher Copyright:© Springer Science+Business Media, LLC 2011.
Keywords
- Bioluminescence tomography
- Orthogonal matching pursuit
- Poisson noise
- Source reconstruction
- ℓ regularization
ASJC Scopus subject areas
- Software
- General Engineering
- Computational Mathematics
- Theoretical Computer Science
- Applied Mathematics
- Numerical Analysis
- Computational Theory and Mathematics