Light Illumination and Detection Patterns for Fluorescence Diffuse Optical Tomography Based on Compressive Sensing

Fluorescence diffuse optical tomography (FDOT) is an emerging molecular imaging modality that uses near infrared light to excite the fluorophore injected into tissue; and to reconstruct the fluorophore concentration from boundary measurements. The FDOT image reconstruction is a highly ill-posed inverse problem due to a large number of unknowns and limited number of measurements. However, the fluorophore distribution is often very sparse in the imaging domain since fluorophores are typically designed to accumulate in relatively small regions. In this paper, we use compressive sensing (CS) framework to design light illumination and detection patterns to improve the reconstruction of sparse fluorophore concentration. Unlike the conventional FDOT imaging where spatially distributed light sources illuminate the imaging domain one at a time and the corresponding boundary measurements are used for image reconstruction, we assume that the light sources illuminate the imaging domain simultaneously several times and the corresponding boundary measurements are linearly filtered prior to image reconstruction. We design a set of optical intensities (illumination patterns) and a linear filter (detection pattern) applied to the boundary measurements to improve the reconstruction of sparse fluorophore concentration maps. We show that the FDOT sensing matrix can be expressed as a columnwise Kronecker product of two matrices determined by the excitation and emission light fields. We derive relationships between the incoherence of the FDOT forward matrix and these two matrices, and use these results to reduce the incoherence of the FDOT forward matrix. We present extensive numerical simulation and the results of a real phantom experiment to demonstrate the improvements in image reconstruction due to the CS-based light illumination and detection patterns in conjunction with relaxation and greedy-type reconstruction algorithms.

[1]  Birsen Yazici,et al.  Discretization Error Analysis and Adaptive Meshing Algorithms for Fluorescence Diffuse Optical Tomography: Part II , 2010, IEEE Transactions on Medical Imaging.

[2]  Birsen Yazici,et al.  Direct Reconstruction of Pharmacokinetic-Rate Images of Optical Fluorophores From NIR Measurements , 2009, IEEE Transactions on Medical Imaging.

[3]  A. Robert Calderbank,et al.  Why Gabor frames? Two fundamental measures of coherence and their role in model selection , 2010, Journal of Communications and Networks.

[4]  Michael Elad,et al.  Sparse and Redundant Representations - From Theory to Applications in Signal and Image Processing , 2010 .

[5]  Wolfgang Bangerth,et al.  Plane-wave fluorescence tomography with adaptive finite elements. , 2006, Optics letters.

[6]  A. Chatziioannou,et al.  Tomographic bioluminescence imaging by use of a combined optical-PET (OPET) system: a computer simulation feasibility study , 2005, Physics in medicine and biology.

[7]  Turgut Durduran,et al.  Compressed sensing in diffuse optical tomography. , 2010, Optics express.

[8]  Jean-Luc Starck,et al.  Sparse Solution of Underdetermined Systems of Linear Equations by Stagewise Orthogonal Matching Pursuit , 2012, IEEE Transactions on Information Theory.

[9]  Vasilis Ntziachristos,et al.  Free-space fluorescence molecular tomography utilizing 360° geometry projections , 2007 .

[10]  Joel A. Tropp,et al.  Greed is good: algorithmic results for sparse approximation , 2004, IEEE Transactions on Information Theory.

[11]  Birsen Yazici,et al.  Discretization Error Analysis and Adaptive Meshing Algorithms for Fluorescence Diffuse Optical Tomography in the Presence of Measurement Noise , 2011, IEEE Transactions on Image Processing.

[12]  Pierre Vandergheynst,et al.  Dictionary Preconditioning for Greedy Algorithms , 2008, IEEE Transactions on Signal Processing.

[13]  Mário A. T. Figueiredo,et al.  Gradient Projection for Sparse Reconstruction: Application to Compressed Sensing and Other Inverse Problems , 2007, IEEE Journal of Selected Topics in Signal Processing.

[14]  Koon-Pong Wong,et al.  An investigation of coded aperture imaging for small animal SPECT , 2000 .

[15]  Deanna Needell,et al.  Uniform Uncertainty Principle and Signal Recovery via Regularized Orthogonal Matching Pursuit , 2007, Found. Comput. Math..

[16]  Michael Elad,et al.  Optimized Projections for Compressed Sensing , 2007, IEEE Transactions on Signal Processing.

[17]  I. Daubechies,et al.  An iterative thresholding algorithm for linear inverse problems with a sparsity constraint , 2003, math/0307152.

[18]  Arye Nehorai,et al.  Image reconstruction for diffuse optical tomography using sparsity regularization and expectation-maximization algorithm. , 2007, Optics express.

[19]  A B Thompson,et al.  Diagnostic imaging of breast cancer using fluorescence-enhanced optical tomography: phantom studies. , 2004, Journal of biomedical optics.

[20]  Xavier Intes,et al.  Time-resolved diffuse optical tomography with patterned-light illumination and detection. , 2010, Optics letters.

[21]  B. Yazıcı,et al.  Two-level domain decomposition methods for diffuse optical tomography , 2006 .

[22]  Turgut Durduran,et al.  Sparse Image Reconstruction in Diffuse Optical Tomography: An Application of Compressed Sensing , 2010 .

[23]  Deanna Needell,et al.  CoSaMP: Iterative signal recovery from incomplete and inaccurate samples , 2008, ArXiv.

[24]  N. Pullman Matrix theory and its applications , 1976 .

[25]  A. Adibi,et al.  Optimal sparse solution for fluorescent diffuse optical tomography: theory and phantom experimental results. , 2007, Applied optics.

[26]  Wolfgang Bangerth,et al.  Non-contact fluorescence optical tomography with scanning patterned illumination. , 2006, Optics express.

[27]  Britton Chance,et al.  An adaptive multigrid algorithm for region of interest diffuse optical tomography , 2003, Proceedings 2003 International Conference on Image Processing (Cat. No.03CH37429).

[28]  David L. Donoho,et al.  Sparse Solution Of Underdetermined Linear Equations By Stagewise Orthogonal Matching Pursuit , 2006 .

[29]  Gongguo Tang,et al.  Verifiable and Computable Performance Analysis of Sparsity Recovery , 2011 .

[30]  Birsen Yazici,et al.  Optimum source design for detection of heterogeneities in diffuse optical imaging , 2006, SPIE BiOS.

[31]  David A Boas,et al.  Assessing the future of diffuse optical imaging technologies for breast cancer management. , 2008, Medical physics.

[32]  Roummel F. Marcia,et al.  Compressive coded apertures for high-resolution imaging , 2010, Photonics Europe.

[33]  Yonina C. Eldar,et al.  Sensing Matrix Optimization for Block-Sparse Decoding , 2010, IEEE Transactions on Signal Processing.

[34]  Jong Chul Ye,et al.  Compressive Diffuse Optical Tomography: Noniterative Exact Reconstruction Using Joint Sparsity , 2011, IEEE Transactions on Medical Imaging.

[35]  Michael Elad,et al.  Stable recovery of sparse overcomplete representations in the presence of noise , 2006, IEEE Transactions on Information Theory.

[36]  Guillaume Bal,et al.  Transport- and diffusion-based optical tomography in small domains: a comparative study. , 2007, Applied optics.

[37]  Vasilis Ntziachristos,et al.  Free-space fluorescence molecular tomography utilizing 360 degrees geometry projections. , 2007, Optics letters.

[38]  Vadim A. Markel,et al.  Optical tomography with structured illumination. , 2009, Optics letters.

[39]  Britton Chance,et al.  Extended Kalman Filtering for the Modeling and Analysis of ICG Pharmacokinetics in Cancerous Tumors Using NIR Optical Methods , 2006, IEEE Transactions on Biomedical Engineering.

[40]  Vasilis Ntziachristos,et al.  Preconditioning of the fluorescence diffuse optical tomography sensing matrix based on compressive sensing. , 2012, Optics letters.

[41]  Henry Arguello,et al.  Code aperture optimization for spectrally agile compressive imaging. , 2011, Journal of the Optical Society of America. A, Optics, image science, and vision.

[42]  Gongguo Tang,et al.  Performance Analysis for Sparse Support Recovery , 2009, IEEE Transactions on Information Theory.

[43]  A. Joshi,et al.  Adaptive finite element based tomography for fluorescence optical imaging in tissue. , 2004, Optics express.

[44]  Jong Chul Ye,et al.  Exact reconstruction formula for diffuse optical tomography using simultaneous sparse representation , 2008, 2008 5th IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[45]  Arkadi Nemirovski,et al.  On verifiable sufficient conditions for sparse signal recovery via ℓ1 minimization , 2008, Math. Program..

[46]  R. Leahy,et al.  Digimouse: a 3D whole body mouse atlas from CT and cryosection data , 2007, Physics in medicine and biology.

[47]  Guillermo Sapiro,et al.  Learning to Sense Sparse Signals: Simultaneous Sensing Matrix and Sparsifying Dictionary Optimization , 2009, IEEE Transactions on Image Processing.

[48]  Birsen Yazici,et al.  Discretization Error Analysis and Adaptive Meshing Algorithms for Fluorescence Diffuse Optical Tomography: Part I , 2010, IEEE Transactions on Medical Imaging.

[49]  Michael A. Saunders,et al.  Atomic Decomposition by Basis Pursuit , 1998, SIAM J. Sci. Comput..

[50]  V. Ntziachristos Fluorescence molecular imaging. , 2006, Annual review of biomedical engineering.