Simultaneous PET Image Reconstruction and Feature Extraction Method using Non-negative, Smooth, and Sparse Matrix Factorization

Positron emission tomography (PET) is an important imaging technique to visualize a number of functions in the brain or human body. For reconstructing PET images from the sinogram data, an inverse problem has to be solved using numerical optimizations such as expectation-maximization (EM)-based methods. However, the standard EM method suffers from measurement noise added in the sinogram data. In this paper, we propose a new simultaneous PET image reconstruction and parts extraction method using constrained non-negative matrix factorization. In contrast that the many existing methods reconstruct a single PET image independently, we reconstruct the time-series of PET images simultaneously from the time-series of sinograms using non-negative matrix factorization. Furthermore, we impose the smoothness constraint for the temporal feature, and the exclusive LASSO-based sparseness constraint for the spatial feature for robust image reconstruction and physically meaningful feature extraction.

[1]  L. Shepp,et al.  Maximum Likelihood Reconstruction for Emission Tomography , 1983, IEEE Transactions on Medical Imaging.

[2]  M. Raichle Behind the scenes of functional brain imaging: a historical and physiological perspective. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[3]  D. Bailey,et al.  The direct calculation of parametric images from dynamic PET data using maximum-likelihood iterative reconstruction. , 1997 .

[4]  W. Klunk,et al.  Imaging brain amyloid in Alzheimer's disease with Pittsburgh Compound‐B , 2004, Annals of neurology.

[5]  A. Lammertsma,et al.  Simplified Reference Tissue Model for PET Receptor Studies , 1996, NeuroImage.

[6]  A. Ramm,et al.  The RADON TRANSFORM and LOCAL TOMOGRAPHY , 1996 .

[7]  Yukihiko Yamashita,et al.  Smooth nonnegative matrix and tensor factorizations for robust multi-way data analysis , 2015, Signal Process..

[8]  Hidekata Hontani,et al.  A robust PET image reconstruction using constrained non-negative matrix factorization , 2017, 2017 Asia-Pacific Signal and Information Processing Association Annual Summit and Conference (APSIPA ASC).

[9]  H. Sebastian Seung,et al.  Learning the parts of objects by non-negative matrix factorization , 1999, Nature.

[10]  Feiping Nie,et al.  Exclusive Feature Learning on Arbitrary Structures via \ell_{1, 2}-norm , 2014, NIPS.

[11]  H. Sebastian Seung,et al.  Algorithms for Non-negative Matrix Factorization , 2000, NIPS.

[12]  D J Brooks,et al.  Comparison of Methods for Analysis of Clinical [11C]Raclopride Studies , 1996, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[13]  F. Natterer The Mathematics of Computerized Tomography , 1986 .

[14]  K. Lange,et al.  EM reconstruction algorithms for emission and transmission tomography. , 1984, Journal of computer assisted tomography.

[15]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[16]  Robert Y. L. Chu,et al.  Radiological Imaging: The Theory of Image Formation, Detection, and Processing. Volume 2 by H. H. Barrett and W. Swindell , 1983 .

[17]  S. Gambhir Molecular imaging of cancer with positron emission tomography , 2002, Nature Reviews Cancer.

[18]  M. Raichle A Brief History of Human Functional Brain Mapping , 2000 .