Estimation of sparse null space functions for compressed sensing in SPECT

Compressed sensing (CS) [1] is a novel sensing (acquisition) paradigm that applies to discrete-to-discrete system models and asserts exact recovery of a sparse signal from far fewer measurements than the number of unknowns [1- 2]. Successful applications of CS may be found in MRI [3, 4] and optical imaging [5]. Sparse reconstruction methods exploiting CS principles have been investigated for CT [6-8] to reduce radiation dose, and to gain imaging speed and image quality in optical imaging [9]. In this work the objective is to investigate the applicability of compressed sensing principles for a faster brain imaging protocol on a hybrid collimator SPECT system. As a proofof- principle we study the null space of the fan-beam collimator component of our system with regards to a particular imaging object. We illustrate the impact of object sparsity on the null space using pixel and Haar wavelet basis functions to represent a piecewise smooth phantom chosen as our object of interest.

[1]  E.J. Candes,et al.  An Introduction To Compressive Sampling , 2008, IEEE Signal Processing Magazine.

[2]  T. Chan,et al.  Source reconstruction for spectrally-resolved bioluminescence tomography with sparse a priori information. , 2009, Optics express.

[3]  Frédéric Lesage,et al.  The Application of Compressed Sensing for Photo-Acoustic Tomography , 2009, IEEE Transactions on Medical Imaging.

[4]  Grant T Gullberg,et al.  Null-space function estimation for the interior problem , 2012, Physics in medicine and biology.

[5]  Xiaochuan Pan,et al.  A constrained, total-variation minimization algorithm for low-intensity x-ray CT. , 2011, Medical physics.

[6]  D. Donoho,et al.  Sparse MRI: The application of compressed sensing for rapid MR imaging , 2007, Magnetic resonance in medicine.

[7]  Charles K. Chui,et al.  An Introduction to Wavelets , 1992 .

[8]  Xiaochuan Pan,et al.  Image reconstruction exploiting object sparsity in boundary-enhanced X-ray phase-contrast tomography. , 2010, Optics express.

[9]  H Barrett,et al.  Decomposition of images and objects into measurement and null components. , 1998, Optics express.

[10]  Emmanuel J. Candès,et al.  Near-Optimal Signal Recovery From Random Projections: Universal Encoding Strategies? , 2004, IEEE Transactions on Information Theory.

[11]  Jie Tang,et al.  Prior image constrained compressed sensing (PICCS): a method to accurately reconstruct dynamic CT images from highly undersampled projection data sets. , 2008, Medical physics.

[12]  Armando Manduca,et al.  Highly Undersampled Magnetic Resonance Image Reconstruction via Homotopic $\ell_{0}$ -Minimization , 2009, IEEE Transactions on Medical Imaging.

[13]  Emil Y. Sidky,et al.  Algorithm-Enabled Low-Dose Micro-CT Imaging , 2011, IEEE Transactions on Medical Imaging.

[14]  Antonin Chambolle,et al.  A First-Order Primal-Dual Algorithm for Convex Problems with Applications to Imaging , 2011, Journal of Mathematical Imaging and Vision.

[15]  Xiaochuan Pan,et al.  First-order convex feasibility algorithms for x-ray CT. , 2012, Medical physics.