Quantifying tissue optical properties of human heads in vivo using continuous-wave near-infrared spectroscopy and subject-specific three-dimensional Monte Carlo models

Monte Carlo-Based Optical Simulation of Optical Distribution in Deep Brain Tissues Using Sixteen Optical Sources

Optical-based imaging has improved from early single-location research to further sophisticated imaging in 2D topography and 3D tomography. These techniques have the benefit of high specificity and non-radiative safety for brain detection and therapy. However, their performance is limited by complex tissue structures. To overcome the difficulty in successful brain imaging applications, we conducted a simulation using 16 optical source types within a brain model that is based on the Monte Carlo method. In addition, we propose an evaluation method of the optical propagating depth and resolution, specifically one based on the optical distribution for brain applications. Based on the results, the best optical source types were determined in each layer. The maximum propagating depth and corresponding source were extracted. The optical source propagating field width was acquired in different depths. The maximum and minimum widths, as well as the corresponding source, were determined. This paper provides a reference for evaluating the optical propagating depth and resolution from an optical simulation aspect, and it has the potential to optimize the performance of optical-based techniques.

Retrieval of chromophore concentration changes in a digital human head model using analytical mean partial pathlengths of photons

Significance

Continuous-wave functional near-infrared spectroscopy has proved to be a valuable tool for assessing hemodynamic activity in the human brain in a non-invasively and inexpensive way. However, most of the current processing/analysis methods assume the head is a homogeneous medium, and hence do not appropriately correct for the signal coming from the scalp. This effect can be reduced by considering light propagation in a layered model of the human head, being the Monte Carlo (MC) simulations the gold standard to this end. However, this implies large computation times and demanding hardware capabilities.

Aim

In this work, we study the feasibility of replacing the homogeneous model and the MC simulations by means of analytical multilayered models, combining in this way, the speed and simplicity of implementation of the former with the robustness and accuracy of the latter.

Approach

Oxy- and deoxyhemoglobin (HbO and HbR, respectively) concentration changes were proposed in two different layers of a magnetic resonance imaging (MRI)-based meshed model of the human head, and then these changes were retrieved by means of (i) a typical homogeneous reconstruction and (ii) a theoretical layered reconstruction.

Results

Results suggest that the use of analytical models of light propagation in layered models outperforms the results obtained using traditional homogeneous reconstruction algorithms, providing much more accurate results for both, the extra- and the cerebral tissues. We also compare the analytical layered reconstruction with MC-based reconstructions, achieving similar degrees of accuracy, especially in the gray matter layer, but much faster (between 4 and 5 orders of magnitude).

Conclusions

We have successfully developed, implemented, and validated a method for retrieving chromophore concentration changes in the human brain, combining the simplicity and speed of the traditional homogeneous reconstruction algorithms with robustness and accuracy much more similar to those provided by MC simulations.

scatterBrains: an open database of human head models and companion optode locations for realistic Monte Carlo photon simulations

Significance

Monte Carlo (MC) simulations are currently the gold standard in the near-infrared and diffuse correlation spectroscopy (NIRS/DCS) communities for generating light transport paths through tissue. However, realistic and diverse models that capture complex tissue layers are not easily available to all; moreover, manually placing optodes on such models can be tedious and time consuming. Such limitations may hinder the adoption of representative models for basic simulations and the use of these models for large-scale simulations, e.g., for training machine learning algorithms.

Aim

We aim to provide the NIRS/DCS communities with an open-source, user-friendly database of morphologically and optically realistic head models, as well as a succinct software pipeline to prepare these models for mesh-based Monte Carlo simulations of light transport.

Approach

Sixteen anatomical models were created from segmented T1-weighted magnetic resonance imaging (MRI) head scans and converted to tetrahedral mesh volumes. Approximately 800 companion scalp surface locations were distributed on each model, comprising full head coverage. A pipeline was created to place custom source and optical detectors at each location, and guidance is provided on how to use these parameters to set up MC simulations.

Results

The models, head surface locations, and all associated code are freely available under the scatterBrains project on Github.

Conclusions

The NIRS/DCS community benefits from having shared resources for conducting MC simulations on realistic head geometries. We hope this will make MRI-based head models and virtual optode placement easily accessible to all. Contributions to the database are welcome and encouraged.

Performance estimation of optical skin probe in short wavelength infrared spectroscopy based on Monte-Carlo simulation

Optical throughput and optical path length are key parameters to obtain high signal to noise ratio and sensor sensitivity for the detection of skin tissue components based on short wavelength infrared (SWIR) spectroscopy. These parameters should be taken into account at the stage of optical system design. We aim to develop a method to estimate the optical efficiency and the effective water path length of a newly designed SWIR spectroscopy skin measurement system using Monte-Carlo photon migration simulation. To estimate the optical efficiency and the effective water path length, we investigated the characteristics of Monte-Carlo photon migration simulation utilizing one layered simple skin model. Simulation of photon transport in skin was conducted for transmission, transflection, and reflection optical configurations in both first overtone (1540 ~ 1820 nm) and combination (2040 ~ 2380 nm) wavelength ranges. Experimental measurement of skin spectrum was done using Fourier transform infrared spectroscopy based system to validate the estimation performance. Overall, the simulated results for optical efficiency and effective water path length are in good agreements with the experimental measurements, which shows the suggested method can be used as a means for the performance estimation and the design optimization of various in-vivo SWIR spectroscopic system.

Special Section Guest Editorial: Introduction to the Special Section Celebrating 30 years of Open Source Monte Carlo Codes in Biomedical Optics

The editorial introduces the JBO Special Section Celebrating 30 Years of Open Source Monte Carlo Codes in Biomedical Optics for Volume 27, Issue 8.