Light propagation through weakly scattering media: a study of Monte Carlo vs. diffusion theory with application to neuroimaging

One of the major challenges within Optical Imaging, photon propagation through clear layers embedded between scattering tissues, can be now efficiently modelled in real-time thanks to the Monte Carlo approach based on GPU. Because of its nature, the photon propagation problem can be very easily parallelized and ran on low cost hardware, avoiding the need for expensive Super Computers. A comparison between Diffusion and MC photon propagation theory is presented in this work with application to neuroimaging, investigating low scattering regions in a mouse-like phantom. Regions such as the Cerebral Spinal Fluid, are currently not taken into account in the classical computational models because of the impossibility to accurately simulate light propagation using fast Diffusive Equation approaches, leading to inaccuracies during the reconstruction process. The goal of the study presented here, is to reduce and further improve the computation accuracy of the reconstructed solution in a highly realistic scenario in the case of neuroimaging in preclinical mouse models.

[1]  Martin Schweiger,et al.  The Toast++ software suite for forward and inverse modeling in optical tomography , 2014, Journal of biomedical optics.

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

[3]  E. Okada,et al.  Hybrid Monte Carlo-diffusion method for light propagation in tissue with a low-scattering region. , 2003, Applied optics.

[4]  Simon R. Arridge,et al.  Development of in-vivo fluorescence imaging with the Matrix-Free method , 2010 .

[5]  S. Arridge Optical tomography in medical imaging , 1999 .

[6]  Dario Fasino,et al.  An inverse Robin problem for Laplace's equation: theoretical results and numerical methods , 1999 .

[7]  M. Schweiger,et al.  Linear and nonlinear reconstruction for optical tomography of phantoms with nonscattering regions. , 2005, Applied optics.

[8]  David A. Boas,et al.  Factors affecting the accuracy of near-infrared spectroscopy concentration calculations for focal changes in oxygenation parameters , 2003, NeuroImage.

[9]  Michael C. Kolios,et al.  Speckle statistics in OCT images: Monte Carlo simulations and experimental studies. , 2014, Optics letters.

[10]  David A Boas,et al.  Monte Carlo simulation of photon migration in 3D turbid media accelerated by graphics processing units. , 2009, Optics express.

[11]  V. Ntziachristos Going deeper than microscopy: the optical imaging frontier in biology , 2010, Nature Methods.

[12]  Stefan Andersson-Engels,et al.  A matrix-free algorithm for multiple wavelength fluorescence tomography. , 2009, Optics express.

[13]  Qianqian Fang,et al.  Mesh-based Monte Carlo method using fast ray-tracing in Plücker coordinates , 2010, Biomedical optics express.