TY - JOUR
T1 - Image reconstruction in fluorescence molecular tomography with sparsity-initialized maximum-likelihood expectation maximization
AU - Zhu, Yansong
AU - Jha, Abhinav K.
AU - Wong, Dean F.
AU - Rahmim, Arman
N1 - Funding Information:
NIH BRAIN Initiative Award (R24 MH106083).
Publisher Copyright:
© 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement.
PY - 2018/7/1
Y1 - 2018/7/1
N2 - We present a reconstruction method involving maximum-likelihood expectation maximization (MLEM) to model Poisson noise as applied to fluorescence molecular tomography (FMT). MLEM is initialized with the output from a sparse reconstruction-based approach, which performs truncated singular value decomposition-based preconditioning followed by fast iterative shrinkage-thresholding algorithm (FISTA) to enforce sparsity. The motivation for this approach is that sparsity information could be accounted for within the initialization, while MLEM would accurately model Poisson noise in the FMT system. Simulation experiments show the proposed method significantly improves images qualitatively and quantitatively. The method results in over 20 times faster convergence compared to uniformly initialized MLEM and improves robustness to noise compared to pure sparse reconstruction. We also theoretically justify the ability of the proposed approach to reduce noise in the background region compared to pure sparse reconstruction. Overall, these results provide strong evidence to model Poisson noise in FMT reconstruction and for application of the proposed reconstruction framework to FMT imaging.
AB - We present a reconstruction method involving maximum-likelihood expectation maximization (MLEM) to model Poisson noise as applied to fluorescence molecular tomography (FMT). MLEM is initialized with the output from a sparse reconstruction-based approach, which performs truncated singular value decomposition-based preconditioning followed by fast iterative shrinkage-thresholding algorithm (FISTA) to enforce sparsity. The motivation for this approach is that sparsity information could be accounted for within the initialization, while MLEM would accurately model Poisson noise in the FMT system. Simulation experiments show the proposed method significantly improves images qualitatively and quantitatively. The method results in over 20 times faster convergence compared to uniformly initialized MLEM and improves robustness to noise compared to pure sparse reconstruction. We also theoretically justify the ability of the proposed approach to reduce noise in the background region compared to pure sparse reconstruction. Overall, these results provide strong evidence to model Poisson noise in FMT reconstruction and for application of the proposed reconstruction framework to FMT imaging.
UR - http://www.scopus.com/inward/record.url?scp=85049372976&partnerID=8YFLogxK
U2 - 10.1364/BOE.9.003106
DO - 10.1364/BOE.9.003106
M3 - Article
AN - SCOPUS:85049372976
SN - 2156-7085
VL - 9
SP - 3106
EP - 3121
JO - Biomedical Optics Express
JF - Biomedical Optics Express
IS - 7
M1 - #327334
ER -