Theoretical investigation of photon partial pathlengths in multilayered turbid media.

Functional near infrared spectroscopy (fNIRS) is a valuable tool for assessing oxy- and deoxyhemoglobin concentration changes (Δ[HbO] and Δ[HbR], respectively) in the human brain. To this end, photon pathlengths in tissue are needed to convert from light attenuation to Δ[HbO] and Δ[HbR]. Current techniques describe the human head as a homogeneous medium, in which case these pathlengths are easily computed. However, the head is more appropriately described as a layered medium; hence, the partial pathlengths in each layer are required. The current way to do this is by means of Monte Carlo (MC) simulations, which are time-consuming and computationally expensive. In this work, we introduce an approach to theoretically calculate these partial pathlengths, which are computed several times faster than MC simulations. Comparison of our approach with MC simulations show very good agreement. Results also suggest that these analytical expressions give much more specific information about light absorption in each layer than in the homogeneous case.

A Monte Carlo study of near infrared light propagation in the human head with lesions-a time-resolved approach.

Several clinical conditions leading to traumatic brain injury can cause hematomas or edemas inside the cerebral tissue. If these are not properly treated in time, they are prone to produce long-term neurological disabilities, or even death. Low-cost, portable and easy-to-handle devices are desired for continuous monitoring of these conditions and Near Infrared Spectroscopy (NIRS) techniques represent an appropriate choice. In this work, we use Time-Resolved (TR) Monte Carlo simulations to present a study of NIR light propagation over a digital MRI phantom. Healthy and injured (hematoma/edema) situations are considered. TR Diffuse Reflectance simulations for different lesion volumes and interoptode distances are performed in the frontal area and the left parietal area. Results show that mean partial pathlengths, photon measurement density functions and time dependent contrasts are sensitive to the presence of lesions, allowing their detection mainly for intermediate optodes separations, which proves that these metrics represent robust means of diagnose and monitoring. Conventional Continuous Wave (CW) contrasts are also presented as a particular case of the time dependent ones, but they result less sensitive to the lesions, and have higher associated uncertainties.

Implementation of the extended Kalman filter for determining the optical and geometrical properties of turbid layered media by time-resolved single distance measurements.

In this article we propose an implementation of the extended Kalman filter (EKF) for the retrieval of optical and geometrical properties in two-layered turbid media assuming a dynamic setting, where absorption of each layer was changed in different steps. Prior works implemented the EKF in frequency-domain with several pairs of light sources and detectors and for static parameters estimation problems. Here we explore the use of the EKF in single distance, time-domain measurements, together with a corresponding forward model. Results show good agreement between retrieved and nominal values, with rather narrow analytical credibility intervals, indicating that the recovery process has low uncertainty, especially for the absorption coefficients.

A comparison between plausible models in layered turbid media with geometrical variations applying a Bayesian selection criterion.

One possible application of Near Infrared techniques is to analyze human brain metabolic activity. Currently used models take into account the layered structure of the human head but, usually, they do not consider the non-planar surface of some of the boundaries, i.e. gray matter, which results in a much more complex structure, thus leading to more sophisticated models and longer calculation times. The main objective of this work is to determine if it is worth to replace a planar layered structure by a non-planar one. To this end we implement a Bayesian-based quantitative methodology for choosing between two competitive models describing light propagation in layered turbid media. Experiments of time-resolved diffuse reflectance measurements are performed in layered phantoms and complemented with numerical calculations. The resulting Distributions of Time of Flight of both models are compared using Bayesian model selection analysis. The non-planar interface was introduced in the simulations by a simple surface parametrization. Results suggest that, under certain conditions, a multilayer model with planar boundaries is good enough.