A parallel adaptive finite element simplified spherical harmonics approximation solver for frequency domain fluorescence molecular imaging.

Fluorescence molecular imaging/tomography may play an important future role in preclinical research and clinical diagnostics. Time- and frequency-domain fluorescence imaging can acquire more measurement information than the continuous wave (CW) counterpart, improving the image quality of fluorescence molecular tomography. Although diffusion approximation (DA) theory has been extensively applied in optical molecular imaging, high-order photon migration models need to be further investigated to match quantitation provided by nuclear imaging. In this paper, a frequency-domain parallel adaptive finite element solver is developed with simplified spherical harmonics (SP(N)) approximations. To fully evaluate the performance of the SP(N) approximations, a fast time-resolved tetrahedron-based Monte Carlo fluorescence simulator suitable for complex heterogeneous geometries is developed using a convolution strategy to realize the simulation of the fluorescence excitation and emission. The validation results show that high-order SP(N) can effectively correct the modeling errors of the diffusion equation, especially when the tissues have high absorption characteristics or when high modulation frequency measurements are used. Furthermore, the parallel adaptive mesh evolution strategy improves the modeling precision and the simulation speed significantly on a realistic digital mouse phantom. This solver is a promising platform for fluorescence molecular tomography using high-order approximations to the radiative transfer equation.

[1]  Haiou Shen,et al.  A tetrahedron-based inhomogeneous Monte Carlo optical simulator , 2010, Physics in medicine and biology.

[2]  L Wang,et al.  MCML--Monte Carlo modeling of light transport in multi-layered tissues. , 1995, Computer methods and programs in biomedicine.

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

[4]  Benjamin S. Kirk,et al.  Library for Parallel Adaptive Mesh Refinement / Coarsening Simulations , 2006 .

[5]  Ivo Babuška,et al.  A posteriori error analysis and adaptive processes in the finite element method: Part I—error analysis , 1983 .

[6]  C. L. Hutchinson,et al.  Fluorescence and Absorption Contrast Mechanisms for Biomedical Optical Imaging Using Frequency‐Domain Techniques , 1997, Photochemistry and photobiology.

[7]  Eva M. Sevick-Muraca,et al.  Non-contact fluorescence optical tomography with scanning area illumination , 2006, 3rd IEEE International Symposium on Biomedical Imaging: Nano to Macro, 2006..

[8]  Benjamin S. Kirk,et al.  Adaptive finite element simulation of flow and transport applications on parallel computers , 2007 .

[9]  Alan Perkins,et al.  Diagnostic imaging , 2012, The Lancet.

[10]  John McGhee,et al.  Radiative transport in fluorescence-enhanced frequency domain photon migration. , 2006, Medical physics.

[11]  Simon R Cherry,et al.  In vivo molecular and genomic imaging: new challenges for imaging physics. , 2004, Physics in medicine and biology.

[12]  Michael Chu,et al.  Light transport in biological tissue using three-dimensional frequency-domain simplified spherical harmonics equations , 2009, Physics in medicine and biology.

[13]  Harvey R Herschman,et al.  Molecular Imaging: Looking at Problems, Seeing Solutions , 2003, Science.

[14]  Greg Turk,et al.  Generating random points in triangles , 1990 .

[15]  Eric C Frey,et al.  Development of a 4-D digital mouse phantom for molecular imaging research. , 2004, Molecular imaging and biology : MIB : the official publication of the Academy of Molecular Imaging.

[16]  S R Cherry,et al.  Optical imaging of Cerenkov light generation from positron-emitting radiotracers , 2009, Physics in medicine and biology.

[17]  Singiresu S. Rao.  The finite element method in engineering , 1982 .

[18]  Hamid Dehghani,et al.  In vivo bioluminescence tomography with a blocking-off finite-difference SP3 method and MRI/CT coregistration. , 2009, Medical physics.

[19]  Huabei Jiang,et al.  A higher order diffusion model for three-dimensional photon migration and image reconstruction in optical tomography , 2009, Physics in medicine and biology.

[20]  Milton V. Marshall,et al.  Imaging of lymph flow in breast cancer patients after microdose administration of a near-infrared fluorophore: feasibility study. , 2008, Radiology.

[21]  Eva M Sevick-Muraca,et al.  Fluorescence-enhanced, near infrared diagnostic imaging with contrast agents. , 2002, Current opinion in chemical biology.

[22]  J. Culver,et al.  Time-dependent whole-body fluorescence tomography of probe bio-distributions in mice. , 2005, Optics express.

[23]  Y Liu,et al.  Flux vector formulation for photon propagation in the biological tissue. , 2007, Optics letters.

[24]  Graham F. Carey,et al.  Computational grids : generation, adaptation, and solution strategies , 1997 .

[25]  S R Arridge,et al.  Recent advances in diffuse optical imaging , 2005, Physics in medicine and biology.

[26]  Yujie Lu,et al.  Experimental bioluminescence tomography with fully parallel radiative-transfer-based reconstruction framework. , 2009, Optics express.

[27]  Vasilis Ntziachristos,et al.  The inverse source problem based on the radiative transfer equation in optical molecular imaging , 2005 .

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

[29]  U. Whitcher.  Boundary Conditions , 2011 .

[30]  Jie Tian,et al.  A multilevel adaptive finite element algorithm for bioluminescence tomography. , 2006, Optics express.

[31]  R. Weissleder.  Molecular imaging: exploring the next frontier. , 1999, Radiology.

[32]  Eva M Sevick-Muraca,et al.  Near-infrared fluorescence contrast-enhanced imaging with intensified charge-coupled device homodyne detection: measurement precision and accuracy. , 2003, Journal of biomedical optics.

[33]  Xuejun Gu,et al.  A parallel reduced-space sequential-quadratic programming algorithm for frequency-domain small animal optical tomography , 2009, BiOS.

[34]  L. O. Svaasand,et al.  Boundary conditions for the diffusion equation in radiative transfer. , 1994, Journal of the Optical Society of America. A, Optics, image science, and vision.

[35]  Monte Carlo simulation of time-dependent, transport-limited fluorescent boundary measurements in frequency domain. , 2007, Medical physics.

[36]  M. Jiang,et al.  Uniqueness theorems in bioluminescence tomography. , 2004, Medical physics.

[37]  B. Rice,et al.  Three-dimensional reconstruction of in vivo bioluminescent sources based on multispectral imaging. , 2007, Journal of biomedical optics.

[38]  Geoffrey McLennan,et al.  Practical reconstruction method for bioluminescence tomography. , 2005, Optics express.

[39]  John McGhee,et al.  Radiative transport-based frequency-domain fluorescence tomography , 2008, Physics in medicine and biology.

[40]  R. Leahy,et al.  Hyperspectral and multispectral bioluminescence optical tomography for small animal imaging , 2005, Physics in medicine and biology.

[41]  Jie Tian,et al.  Spectrally resolved bioluminescence tomography with the third-order simplified spherical harmonics approximation. , 2009, Physics in medicine and biology.

[42]  Jae Hoon Lee,et al.  Fast intersections on nested tetrahedrons (FINT): An algorithm for adaptive finite element based distributed parameter estimation , 2008, J. Comput. Phys..

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

[44]  M. Eppstein,et al.  Fluorescence-enhanced optical imaging in large tissue volumes using a gain-modulated ICCD camera. , 2003, Physics in medicine and biology.

[45]  I-Chih Tan,et al.  Lymphatic imaging in humans with near-infrared fluorescence. , 2009, Current opinion in biotechnology.

[46]  S. Gambhir,et al.  Molecular imaging in living subjects: seeing fundamental biological processes in a new light. , 2003, Genes & development.

[47]  Akira Ishimaru,et al.  Wave propagation and scattering in random media , 1997 .

[48]  B. Pogue,et al.  Spectrally resolved bioluminescence optical tomography. , 2006, Optics letters.

[49]  Arion F. Chatziioannou,et al.  Čerenkov radiation imaging as a method for quantitative measurements of beta particles in a microfluidic chip , 2008 .

[50]  E. Lewis,et al.  Computational Methods of Neutron Transport , 1993 .

[51]  C. Contag,et al.  Advances in in vivo bioluminescence imaging of gene expression. , 2002, Annual review of biomedical engineering.

[52]  A G Yodh,et al.  Fluorescent heterogeneities in turbid media: limits for detection, characterization, and comparison with absorption. , 1998, Applied optics.

[53]  Antonio Pifferi,et al.  Accelerated Monte Carlo models to simulate fluorescence spectra from layered tissues. , 2003, Journal of the Optical Society of America. A, Optics, image science, and vision.

[54]  Michael Z. Lin,et al.  Mammalian Expression of Infrared Fluorescent Proteins Engineered from a Bacterial Phytochrome , 2009, Science.

[55]  Edward W. Larsen,et al.  Light transport in biological tissue based on the simplified spherical harmonics equations , 2006, J. Comput. Phys..

[56]  Paolo Cignoni,et al.  Generating Random Points in a Tetrahedron , 2000, J. Graphics, GPU, & Game Tools.

[57]  Yujie Lu,et al.  A Parallel Adaptive Finite Element Method for the Simulation of Photon Migration with the Radiative-Transfer-Based Model. , 2009, Communications in numerical methods in engineering.

[58]  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.

[59]  D. Boas,et al.  Three dimensional Monte Carlo code for photon migration through complex heterogeneous media including the adult human head. , 2002, Optics express.

[60]  George Karypis,et al.  Multilevel k-way Partitioning Scheme for Irregular Graphs , 1998, J. Parallel Distributed Comput..