Open Access
5 October 2023 scatterBrains: an open database of human head models and companion optode locations for realistic Monte Carlo photon simulations
Melissa M. Wu, Roarke W. Horstmeyer, Stefan A. Carp
Author Affiliations +
Abstract

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.

1.

Introduction

In the past several decades, Monte Carlo (MC) light simulations have emerged as the gold standard technique for simulating photon transport through tissue.1,2 Advances in these algorithms have allowed for modeling increasingly complex transport geometries or tissue types39 while drastically reducing the runtime through parallel or cloud computing architectures.1014 Such developments have spawned new studies that use MC simulations to, among many other efforts, investigate potential tissue areas to image,1518 determine the promise of new imaging modalities,1922 optimize hardware in existing modalities,2328 improve recovery of desired biomarkers,29,30 and further elucidate the physical principles behind light-tissue interaction.3135 This list is far from exhaustive—the number of distinctive studies exploiting MC simulations has truly burgeoned across subfields in the past couple of decades.2 For details on methodology, the interested reader can peruse Zhu and Liu’s review paper36 on MC simulations in biological tissue, although much work has been developed since then as well.

Many works in the near-infrared and diffuse correlation spectroscopy (NIRS/DCS) communities have taken advantage of the Monte Carlo eXtreme (MCX)37 or mesh-based Monte Carlo (MMC)38 software packages written by Qianqian Fang and his team.26,30,32,3941 Both packages have the capability to simulate arbitrary optode locations on any given complex, heterogeneous tissue model. However, many recent works using MCX/MMC have still used atlas-based head models (or even simpler slab or sphere geometries) while restricting the simulations to a fixed probe configuration at a single location.24,30,32,4244 This is driven in part by difficulties in manually placing optodes and potentially limits the works’ investigative scope. Recently, both Nizam et al.45 and Bürmen et al.46 have acknowledged the need in the field for open reference datasets, as well as user-friendly multiple optode placement, for both machine learning training and validation of custom MC software. Both groups have provided datasets based on slab geometries with embedded absorbing targets of various shapes meant for training tomographic reconstruction algorithms. However, these are unrepresentative of realistic head, skull, and brain geometries and would not be as useful for the growing segment of the NIRS/DCS community targeting noninvasive detection of brain hemodynamics. Although a myriad of NIRS toolboxes exist,47 almost all are focused on signal or data analysis as their primary functionality and thus may often provide a single atlas-based head model. For users looking to run large-scale simulations considering many optode locations on diverse, realistic head models, such software may not be suitable.

In this work, we present scatterBrains, an open database of human head models along with companion optode locations of interest, and a toolbox to create specifications for running MC simulations of light transport, including the code to create input files compatible with MMC. The head models were derived from T1-weighted magnetic resonance imaging (MRI) scans acquired from healthy adult volunteers in the course of a previous study30 and subsequently segmented into five tissue layers. Each model was then “polished,” i.e., gaps or unwanted protrusions as a result of imperfect segmentation were manually smoothed over. The models were then converted into tetrahedral meshes using the iso2mesh toolbox.48 For the straightforward selection of a surface location of interest to run MMC simulations, 805 predefined scalp locations per model, comprising full head coverage, are provided for the user.16 Also included is a pipeline to automatically place user-specified optodes at any of these (or other) scalp locations and to write these and other parameters to MMC input files. Example post-processing for DCS is also given. The head models, their companion surface locations, and the processing pipeline comprise the scatterBrains open database.

2.

Database Description and Methods

As mentioned above, the MRI structural scans that formed the basis for the head models in the toolbox were acquired in a previous research study approved by the Mass General Brigham Institutional Review Board (MGB IRB). Further, the MGB IRB approved the de-identification procedures and the sharing of the de-identified models in this open dataset.

The repository consists of 16 separate folders, 1 for each subject, as well as some functions and example scripts. Each subject folder contains files for the head model (both mesh- and volume-based) and the companion surface locations. Dependencies are described in Sec. 3.

2.1.

Subject Head Models

The methodology for head model generation is described in our previous paper.16 Each model is segmented into five tissue layers representing the scalp, skull, cerebrospinal fluid, gray matter, and white matter. The head models were de-identified by removing the face area. The ears were also removed for more accurate optode placement across the head (users can contact the authors for models with the ears if needed). The details of the subject information is given in Table 1. Additional subject information regarding extracerebral thickness and cerebral sensitivity statistics are provided in Supplementary Material.

Table 1

Subject demographics.

Subject IDSexAgeRaceHispanic/Latino
1M24More than oneUnknown/not reported
2M29WhiteNo
3F24WhiteYes
4F24AsianNo
5F25WhiteNo
6F21Black or AfricanNo
7F39WhiteNo
8F32WhiteNo
9F23Black or AfricanNo
10F24WhiteNo
11M30AsianNo
12M48WhiteNo
13F23More than oneYes
14F24WhiteNo
15F38WhiteNo
16F25WhiteNo

2.2.

Companion Surface Locations

The methodology for surface location placement is also described in our previous work cited above. For this dataset, an additional 6 points near each ear were removed, resulting in a total of 805 surface locations for each model. A visualization of a subject head model and its companion surface locations is shown in Fig. 1(a).

Fig. 1

Pipeline demonstrating usage of scatterBrains functionality. (a) An example subject head with surface locations and a cross-section of the internal anatomy is shown. The subject head is lying supine with the de-identified face directed upward. The tan, gray, blue, pink, and red surfaces are the scalp, skull, cerebrospinal fluid, gray matter, and white matter, respectively. The user can select any one of the predefined surface locations or choose their own. (b) The optode placement algorithm will automatically generate source and detector positions at the surface location based on user-inputted distances. (c) A MC light propagation simulation can be run with the MMC software. scatterBrains provides all of the input information necessary for the simulation, and can also generate the input file directly passed to MMC. (d) The user can process the output photon trajectory data from the simulation. Analysis steps for DCS are also provided in the toolbox.

JBO_28_10_100501_f001.png

2.3.

Associated Functionality

Included in the toolbox is a pipeline for automatically placing sources and detectors at a given pre-defined head surface location [Fig. 1(b)]. We briefly describe the methodology here and denote the selected head surface location as the “point-of-interest” (POI). First, the intersection of a plane, defined by the POI and a point near each ear (also provided for the user), and the head mesh is taken. A segment of the intersection perimeter surrounding the POI is extracted, interpolated, and smoothed. The sources and detectors are then defined as points along this segment, per user-specified distances, with the POI centered between the source and furthest detector. Of note, the source–detector distances are calculated along the curved arc of the segment (not by the Euclidean distance)—we believe that this is more physically accurate for optical probes that are designed on a flat strip first and then wrapped around the head. Figure 2 shows the different steps of the optode placement pipeline. Users need only to input their desired POI (which can be visualized and selected from 1 of the provided 805 surface points) and the desired source–detector distances. The placement in the example below takes less than a second to run.

Fig. 2

Optode placement pipeline for an example POI is shown. Here, the subject head is lying supine, with the facial features removed. The chosen POI is located above the subject eyebrow, slightly toward the subject’s right side. (a) A mesh-plane intersection is made with the plane defined by the POI and two points near the ears (not shown, but provided for the user). Its perimeter is extracted and shown. (b) An outer segment of this intersection perimeter surrounding the POI is taken. (c) Optodes are placed along this segment according to the desired distances inputted by the user, with the POI centered between the source and maximum detector distance. The red dot is the source, and the blue dots are detectors at 5 and 50 mm (the detector distance has been exaggerated for improved visualization).

JBO_28_10_100501_f002.png

scatterBrains provides the code to write the optode locations, as well as other simulation parameters (e.g., optical properties for each tissue layer), directly to an input file, which can be passed to MMC along with the provided subject model. After running the simulation [Fig. 1(c)], the user can process the photon trajectory results as they wish [Fig. 1(d)]. The details of these steps and other remaining functionality are straightforward and presented in Sec. 3.

3.

Usage

The repository hosting scatterBrains is publicly available on Github at https://github.com/wumelissa/scatterBrains and can be run out-of-the-box. We have tested the code on the Windows and Linux platforms—the repository and dependencies are cross-platform and should readily run on MacOS as well. The data and functions are written in Matlab and require the iso2mesh toolbox48 to be installed to run. The MMC package38 is needed to actually execute the MC simulation and provide photon history files for post-processing. The algorithm shown in Fig. 2 can be run on any of the surface locations in Fig. 1(a), for any combination of source–detector distances.

An end-to-end example script is provided for the user to demonstrate loading a model, selecting a POI, and placing optodes. The main steps of the complete pipeline demonstrated in the script are as follows:

  • Step 1. Select subject model

  • Step 2. Visualize model, and select surface location [Fig. 1(a)]

  • Step 3. Specify desired detector distances, and place optodes [automated, Fig. 1(b)]

  • Step 4. Specify additional MC parameters

  • Step 5. Specify desired tissue optical properties

  • Step 6. Run MMC [Fig. 1(c)]

  • Step 7. (DCS processing): Generate photon autocorrelations from MMC output [Fig. 1(d)]

  • Step 8. (DCS processing): Fit autocorrelations to semi-infinite model [Fig. 1(d)].

Steps 1–3 are described in Sec. 2. For Step 4, the code to generate the MMC input file (.json format) is provided, i.e., users only need to fill in their desired parameters in the Matlab script. The same applies for Step 5, for which the code is provided to generate an MMC optical properties file (.dat format). During Step 6, users can then run a simulation with their own MMC installation—MMC and iso2mesh are not redistributed in the scatterBrains repository because they are freely available. In Steps 7 and 8, post-processing code for generating diffuse correlation spectroscopy simulated measurements is also provided.

Aside from this end-to-end demo, additional example scripts are included to calculate the distance of the POI to the nearest brain tissue and to generate volume and surface meshes from the original segmented volume file. A Google groups forum is provided for community discussions and author support, although MMC-specific questions will be directed toward the MMC community support forum.

4.

Discussion

In this work, we present an open dataset of realistic head models, companion surface locations, and code for user-friendly input file generation suitable for running MC simulations using the MMC package. As mentioned above, open datasets of slab tissue models with embedded blood vessels have recently been released for ease of large-scale MC simulations.45,46 However, these are primarily targeted toward diffuse optical tomography applications or smaller-scale spectroscopy measurements. For NIRS/DCS brain imaging, in which vessel reconstruction is not the main goal, modeling other tissue features (e.g., extracerebral layers, gyri and sulci of the brain) would be more valuable. Moreover, as measurements with larger source–detector separations (and thus a probed tissue area) are increasingly attained,49 larger and more comprehensive models are needed to capture photon paths.32 As mentioned previously, the NIRS community has many open-source toolboxes,47 primarily for signal and data analysis (as opposed to being database repositories); toolboxes such as AtlasViewer,50 NIRSTORM,51 and NIRS-SPM52 will generally provide a single model atlas, e.g., Colin27.53 If custom probe placement functionality is available in such toolboxes, it is generally restricted to graphical user interface (GUI)-level interaction, which may make preparing and running batch simulations tedious. We anticipate that our open dataset, surface locations, and optode placement functionality can be used in conjunction with these existing toolboxes for analyzing diverse subject models and surface locations. For DCS, toolboxes or datasets such as these have yet to be released—recent works involving MC simulations have often used layered slab or spherical models.26,30,32,44 Thus, for future DCS studies involving MC simulations, we expect that scatterBrains will be particularly useful. For both NIRS and DCS, we expect that our database will also be of use for large-scale machine learning training, in which thousands of optode configurations across the entire head can be automatically generated on a diverse set of subjects. Finally, we hope that scatterBrains can also be of use for cross-validating results from custom MC simulation software.

There are several means by which scatterBrains can be expanded upon in the near future. The current models are derived from predominantly female subjects in their 20s and 30s (driven by the demographics of those who responded to the study advertisement). A more diverse database in age, sex, and race/ethnicity would be beneficial for studying light transport in a wider range of subjects. Structural MRI scans are freely available to the public,54 so we anticipate no barrier in achieving such an improvement. With potential future expansion of the database, we additionally envision compiling “averaged” head models based on sex, age, etc., for use as references. On the pre-processing side, our group plans to eventually include both the image segmentation55 and surface location placement code (Secs. 2.1 and 2.2) available as part of the toolbox, so users can import their own MRI scans and subsequently use them as model inputs to MMC. For optode placement, we also aim to incorporate more complex optode placement schemes in the automated pipeline, as recent work has exploited non-linear setups in source–detector placement.56 With regard to DCS post-processing, we also plan to include MC fitting methods30 for DCS data, as well as other potential fitting schemes that have demonstrated improvement over the semi-infinite diffuse correlation equation. 41,57 Finally, because NIRS/DCS have also been used for blood flow monitoring in other tissue areas,49,58,59 in the long-term, we hope that scatterBrains will be a launching pad for open database generation of various organ models to use for MC simulations.

As an open database, we encourage users to contribute models and/or processing code to scatterBrains. Our goal is for this database to be easily accessible, user-friendly, and representative in its models. Although there is still work to be done in the future, we hope that collaboration with the NIRS/DCS community will help us toward realizing these aims.

5.

Conclusion

We have presented an open database of human head models and companion surface locations derived from 16 adult subjects. For each model, the code for automatic optode placement at any of the 805 pre-defined surface locations can allow for straightforward and fast input file generation to MMC simulation software. We hope to expand this database and associated functionalities in the future.

Disclosures

The authors have no relevant financial interests in the paper at this time and no other potential conflicts of interest.

Code, Data, and Materials Availability

All code and data in this paper are publicly available at https://github.com/wumelissa/scatterBrains

Acknowledgments

This work could not have been done without access to iso2mesh and MMC. We would like to thank Fang and his group for developing and releasing these tools and for his advice on preparing this open source database. This work was supported by the National Institutes of Health (Grant No. R01NS100750) (SAC).

References

1. 

L. Wang, S. L. Jacques and L. Zheng, “MCML—Monte Carlo modeling of light transport in multi-layered tissues,” Comput. Methods Programs Biomed., 47 (2), 131 –146 https://doi.org/10.1016/0169-2607(95)01640-F CMPBEK 0169-2607 (1995). Google Scholar

2. 

Q. Fang, F. Martelli and L. Lilge, “Special section guest editorial: introduction to the special section celebrating 30 years of Open Source Monte Carlo Codes in biomedical optics,” J. Biomed. Opt., 27 (8), 083001 https://doi.org/10.1117/1.JBO.27.8.083001 JBOPFO 1083-3668 (2022). Google Scholar

3. 

Y. Zhang and Q. Fang, “BlenderPhotonics: an integrated open-source software environment for three-dimensional meshing and photon simulations in complex tissues,” J. Biomed. Opt., 27 (8), 083014 https://doi.org/10.1117/1.JBO.27.8.083014 JBOPFO 1083-3668 (2022). Google Scholar

4. 

S. Yan and Q. Fang, “Hybrid mesh and voxel based Monte Carlo algorithm for accurate and efficient photon transport modeling in complex bio-tissues,” Biomed. Opt. Express, 11 (11), 6262 –6270 https://doi.org/10.1364/BOE.409468 BOEICL 2156-7085 (2020). Google Scholar

5. 

S. Yan et al., “Graphics-processing-unit-accelerated Monte Carlo simulation of polarized light in complex three-dimensional media,” J. Biomed. Opt., 27 (8), 083015 https://doi.org/10.1117/1.JBO.27.8.083015 JBOPFO 1083-3668 (2022). Google Scholar

6. 

J. Gröhl et al., “SIMPA: an open-source toolkit for simulation and image processing for photonics and acoustics,” J. Biomed. Opt., 27 (8), 083010 https://doi.org/10.1117/1.JBO.27.8.083010 JBOPFO 1083-3668 (2022). Google Scholar

7. 

C. K. Hayakawa et al., “MCCL: an open-source software application for Monte Carlo simulations of radiative transport,” J. Biomed. Opt., 27 (8), 083005 https://doi.org/10.1117/1.JBO.27.8.083005 JBOPFO 1083-3668 (2022). Google Scholar

8. 

L. McMillan, G. D. Bruce and K. Dholakia, “Meshless Monte Carlo radiation transfer method for curved geometries using signed distance functions,” J. Biomed. Opt., 27 (8), 083003 https://doi.org/10.1117/1.JBO.27.8.083003 JBOPFO 1083-3668 (2022). Google Scholar

9. 

J. Jönsson and E. Berrocal, “Multi-scattering software: part I: online accelerated Monte Carlo simulation of light transport through scattering media,” Opt. Express, 28 (25), 37612 –37638 https://doi.org/10.1364/OE.404005 OPEXFF 1094-4087 (2020). Google Scholar

10. 

L. Yu et al., “Scalable and massively parallel Monte Carlo photon transport simulations for heterogeneous computing platforms,” J. Biomed. Opt., 23 (1), 010504 https://doi.org/10.1117/1.JBO.23.1.010504 JBOPFO 1083-3668 (2018). Google Scholar

11. 

Q. Fang and S. Yan, “MCX Cloud—a modern, scalable, high-performance and in-browser Monte Carlo simulation platform with cloud computing,” J. Biomed. Opt., 27 (8), 083008 https://doi.org/10.1117/1.JBO.27.8.083008 JBOPFO 1083-3668 (2022). Google Scholar

12. 

Q. Fang and S. Yan, “Graphics processing unit-accelerated mesh-based Monte Carlo photon transport simulations,” J. Biomed. Opt., 24 (11), 115002 https://doi.org/10.1117/1.JBO.24.11.115002 JBOPFO 1083-3668 (2019). Google Scholar

13. 

M. Raayai Ardakani et al., “Framework for denoising Monte Carlo photon transport simulations using deep learning,” J. Biomed. Opt., 27 (8), 083019 https://doi.org/10.1117/1.JBO.27.8.083019 JBOPFO 1083-3668 (2022). Google Scholar

14. 

S. Wang et al., “Scalable and accessible personalized photodynamic therapy optimization with FullMonte and PDT-SPACE,” J. Biomed. Opt., 27 (8), 083006 https://doi.org/10.1117/1.JBO.27.8.083006 JBOPFO 1083-3668 (2022). Google Scholar

15. 

J. Guo et al., “Non-invasive optical monitoring of human lungs: Monte Carlo modeling of photon migration in Visible Chinese Human and an experimental test on a human,” Biomed. Opt. Express, 13 (12), 6389 –6403 https://doi.org/10.1364/BOE.472530 BOEICL 2156-7085 (2022). Google Scholar

16. 

M. M. Wu et al., “Complete head cerebral sensitivity mapping for diffuse correlation spectroscopy using subject-specific magnetic resonance imaging models,” Biomed. Opt. Express, 13 (3), 1131 –1151 https://doi.org/10.1364/BOE.449046 BOEICL 2156-7085 (2022). Google Scholar

17. 

M. Jermyn et al., “Intraoperative brain cancer detection with Raman spectroscopy in humans,” Sci. Transl. Med., 7 (274), 274ra19 https://doi.org/10.1126/scitranslmed.aaa2384 STMCBQ 1946-6234 (2015). Google Scholar

18. 

M. Moradi and Y. Chen, “Monte Carlo simulation of diffuse optical spectroscopy for 3D modeling of dental tissues,” Sensors, 23 (11), 5118 https://doi.org/10.3390/s23115118 SNSRES 0746-9462 (2023). Google Scholar

19. 

X. Cheng et al., “Development of a Monte Carlo-wave model to simulate time domain diffuse correlation spectroscopy measurements from first principles,” J. Biomed. Opt., 27 (8), 083009 https://doi.org/10.1117/1.JBO.27.8.083009 JBOPFO 1083-3668 (2022). Google Scholar

20. 

J. Mao et al., “Monte Carlo-based full-wavelength simulator of Fourier-domain optical coherence tomography,” Biomed. Opt. Express, 13 (12), 6317 –6334 https://doi.org/10.1364/BOE.475428 BOEICL 2156-7085 (2022). Google Scholar

21. 

T. Boonya-ananta et al., “Synthetic photoplethysmography (PPG) of the radial artery through parallelized Monte Carlo and its correlation to body mass index (BMI),” Sci. Rep., 11 (1), 2570 https://doi.org/10.1038/s41598-021-82124-4 SRCEC3 2045-2322 (2021). Google Scholar

22. 

J. Fine, M. J. McShane and G. L. Coté, “Monte Carlo method for assessment of a multimodal insertable biosensor,” J. Biomed. Opt., 27 (8), 083017 https://doi.org/10.1117/1.JBO.27.8.083017 JBOPFO 1083-3668 (2022). Google Scholar

23. 

N. Okui and E. Okada, “Wavelength dependence of crosstalk in dual-wavelength measurement of oxy- and deoxy-hemoglobin,” J. Biomed. Opt., 10 (1), 011015 https://doi.org/10.1117/1.1846076 JBOPFO 1083-3668 (2005). Google Scholar

24. 

D. Mazumder et al., “Optimization of time domain diffuse correlation spectroscopy parameters for measuring brain blood flow,” Neurophotonics, 8 (3), 035005 https://doi.org/10.1117/1.NPh.8.3.035005 (2021). Google Scholar

25. 

G. Strangman, M. A. Franceschini and D. A. Boas, “Factors affecting the accuracy of near-infrared spectroscopy concentration calculations for focal changes in oxygenation parameters,” NeuroImage, 18 (4), 865 –879 https://doi.org/10.1016/S1053-8119(03)00021-1 NEIMEF 1053-8119 (2003). Google Scholar

26. 

H. Zhao and E. M. Buckley, “Influence of source–detector separation on diffuse correlation spectroscopy measurements of cerebral blood flow with a multilayered analytical model,” Neurophotonics, 9 (3), 035002 https://doi.org/10.1117/1.NPh.9.3.035002 (2022). Google Scholar

27. 

S. Brigadoi and R. J. Cooper, “How short is short? Optimum source–detector distance for short-separation channels in functional near-infrared spectroscopy,” Neurophotonics, 2 (2), 025005 https://doi.org/10.1117/1.NPh.2.2.025005 (2015). Google Scholar

28. 

Ajmal et al., “Monte Carlo analysis of optical heart rate sensors in commercial wearables: the effect of skin tone and obesity on the photoplethysmography (PPG) signal,” Biomed. Opt. Express, 12 (12), 7445 –7457 https://doi.org/10.1364/BOE.439893 BOEICL 2156-7085 (2021). Google Scholar

29. 

K. B. Sung and H. H. Chen, “Enhancing the sensitivity to scattering coefficient of the epithelium in a two-layered tissue model by oblique optical fibers: Monte Carlo study,” J. Biomed. Opt., 17 (10), 107003 https://doi.org/10.1117/1.JBO.17.10.107003 JBOPFO 1083-3668 (2012). Google Scholar

30. 

M. M. Wu et al., “Improved accuracy of cerebral blood flow quantification in the presence of systemic physiology cross-talk using multi-layer Monte Carlo modeling,” Neurophotonics, 8 (1), 015001 https://doi.org/10.1117/1.NPh.8.1.015001 (2021). Google Scholar

31. 

T. C. Kao and K. B. Sung, “Quantifying tissue optical properties of human heads in vivo using continuous-wave near-infrared spectroscopy and subject-specific three-dimensional Monte Carlo models,” J. Biomed. Opt., 27 (8), 083021 https://doi.org/10.1117/1.JBO.27.8.083021 JBOPFO 1083-3668 (2022). Google Scholar

32. 

H. Zhao and E. M. Buckley, “Influence of oversimplifying the head anatomy on cerebral blood flow measurements with diffuse correlation spectroscopy,” Neurophotonics, 10 (1), 015010 https://doi.org/10.1117/1.NPh.10.1.015010 (2023). Google Scholar

33. 

F. B. Haeussinger et al., “Simulation of near-infrared light absorption considering individual head and prefrontal cortex anatomy: implications for optical neuroimaging,” PLoS One, 6 (10), e26377 https://doi.org/10.1371/journal.pone.0026377 POLNCL 1932-6203 (2011). Google Scholar

34. 

L. Pattelli and G. Mazzamuto, “Experimental imaging and Monte Carlo modeling of ultrafast pulse propagation in thin scattering slabs,” J. Biomed. Opt., 27 (8), 083020 https://doi.org/10.1117/1.JBO.27.8.083020 JBOPFO 1083-3668 (2022). Google Scholar

35. 

S. You et al., “High cell density and high-resolution 3D bioprinting for fabricating vascularized tissues,” Sci. Adv., 9 (8), eade7923 https://doi.org/10.1126/sciadv.ade7923 STAMCV 1468-6996 (2023). Google Scholar

36. 

C. Zhu and Q. Liu, “Review of Monte Carlo modeling of light transport in tissues,” J. Biomed. Opt., 18 (5), 050902 https://doi.org/10.1117/1.JBO.18.5.050902 JBOPFO 1083-3668 (2013). Google Scholar

37. 

Q. Fang and D. A. Boas, “Monte Carlo simulation of photon migration in 3D turbid media accelerated by graphics processing units,” Opt. Express, 17 (22), 20178 https://doi.org/10.1364/OE.17.020178 OPEXFF 1094-4087 (2009). Google Scholar

38. 

Q. Fang, “Mesh-based Monte Carlo method using fast ray-tracing in Plücker coordinates,” Biomed. Opt. Express, 1 (1), 165 –175 https://doi.org/10.1364/BOE.1.000165 BOEICL 2156-7085 (2010). Google Scholar

39. 

F. Scholkmann et al., “A review on continuous wave functional near-infrared spectroscopy and imaging instrumentation and methodology,” NeuroImage, 85 6 –27 https://doi.org/10.1016/j.neuroimage.2013.05.004 NEIMEF 1053-8119 (2014). Google Scholar

40. 

H. Ayaz et al., “Optical imaging and spectroscopy for the study of the human brain: status report,” Neurophotonics, 9 (S2), S24001 https://doi.org/10.1117/1.NPh.9.S2.S24001 (2022). Google Scholar

41. 

H. Zhao, E. Sathialingam and E. M. Buckley, “Accuracy of diffuse correlation spectroscopy measurements of cerebral blood flow when using a three-layer analytical model,” Biomed. Opt. Express, 12 (11), 7149 –7161 https://doi.org/10.1364/BOE.438303 BOEICL 2156-7085 (2021). Google Scholar

42. 

Y. Gao et al., “Short-separation regression incorporated diffuse optical tomography image reconstruction modeling for high-density functional near-infrared spectroscopy,” Neurophotonics, 10 (2), 025007 https://doi.org/10.1117/1.NPh.10.2.025007 (2023). Google Scholar

43. 

X. Zhai et al., “Brain space image reconstruction of functional near-infrared spectroscopy using a Bayesian adaptive fused sparse overlapping group lasso model,” Neurophotonics, 10 (2), 023516 https://doi.org/10.1117/1.NPh.10.2.023516 (2023). Google Scholar

44. 

M. Helton et al., “Numerical approach to quantify depth-dependent blood flow changes in real-time using the diffusion equation with continuous-wave and time-domain diffuse correlation spectroscopy,” Biomed. Opt. Express, 14 (1), 367 –384 https://doi.org/10.1364/BOE.469419 BOEICL 2156-7085 (2023). Google Scholar

45. 

N. I. Nizam et al., “Monte Carlo-based data generation for efficient deep learning reconstruction of macroscopic diffuse optical tomography and topography applications,” J. Biomed. Opt., 27 (8), 083016 https://doi.org/10.1117/1.JBO.27.8.083016 JBOPFO 1083-3668 (2022). Google Scholar

46. 

M. Bürmen, F. Pernuš and P. Naglič, “MCDataset: a public reference dataset of Monte Carlo simulated quantities for multilayered and voxelated tissues computed by massively parallel PyXOpto Python package,” J. Biomed. Opt., 27 (8), 083012 https://doi.org/10.1117/1.JBO.27.8.083012 JBOPFO 1083-3668 (2022). Google Scholar

47. 

“OpenFNIRS – an fNIRS hardware and software ecosystem,” https://openfnirs.org/software/ (). Google Scholar

48. 

Q. Fang and D. A. Boas, “Tetrahedral mesh generation from volumetric binary and grayscale images,” in IEEE Int. Symp. on Biomed. Imaging: From Nano to Macro, 1142 –1145 (2009). https://doi.org/10.1109/ISBI.2009.5193259 Google Scholar

49. 

A. Pifferi et al., “Initial non-invasive in vivo sensing of the lung using time domain diffuse optics,” http://arxiv.org/abs/2205.08211 (2022). Google Scholar

50. 

C. M. Aasted et al., “Anatomical guidance for functional near-infrared spectroscopy: AtlasViewer tutorial,” Neurophotonics, 2 (2), 020801 https://doi.org/10.1117/1.NPh.2.2.020801 (2015). Google Scholar

51. 

, published online March 21, 2023. https://github.com/Nirstorm/nirstorm (). Google Scholar

52. 

J. C. Ye et al., “NIRS-SPM: statistical parametric mapping for near-infrared spectroscopy,” NeuroImage, 44 (2), 428 –447 https://doi.org/10.1016/j.neuroimage.2008.08.036 NEIMEF 1053-8119 (2009). Google Scholar

53. 

C. J. Holmes et al., “Enhancement of MR images using registration for signal averaging,” J. Comput. Assist. Tomogr., 22 (2), 324 https://doi.org/10.1097/00004728-199803000-00032 JCATD5 0363-8715 (1998). Google Scholar

54. 

C. Markiewicz et al., “The OpenNeuro resource for sharing of neuroscience data,” eLife, 10 e71774 https://doi.org/10.7554/eLife.71774 (2021). Google Scholar

55. 

K. L. Perdue and S. G. Diamond, “T1 magnetic resonance imaging head segmentation for diffuse optical tomography and electroencephalography,” J. Biomed. Opt., 19 (2), 026011 https://doi.org/10.1117/1.JBO.19.2.026011 JBOPFO 1083-3668 (2014). Google Scholar

56. 

S. Xu et al., “Imaging dynamics beneath turbid media via parallelized single-photon detection,” Adv. Sci., 9 (24), 2201885 https://doi.org/10.1002/advs.202201885 (2022). Google Scholar

57. 

W. B. Baker et al., “Pressure modulation algorithm to separate cerebral hemodynamic signals from extracerebral artifacts,” Neurophotonics, 2 (3), 035004 https://doi.org/10.1117/1.NPh.2.3.035004 (2015). Google Scholar

58. 

Y. Shang, K. Gurley and G. Yu, “Diffuse correlation spectroscopy (DCS) for assessment of tissue blood flow in skeletal muscle: recent progress,” Anat. Physiol., 3 (2), 128 https://doi.org/10.4172/2161-0940.1000128 (2013). Google Scholar

59. 

M. Ferrari and V. Quaresima, “A brief review on the history of human functional near-infrared spectroscopy (fNIRS) development and fields of application,” NeuroImage, 63 (2), 921 –935 https://doi.org/10.1016/j.neuroimage.2012.03.049 NEIMEF 1053-8119 (2012). Google Scholar
CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 International License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Melissa M. Wu, Roarke W. Horstmeyer, and Stefan A. Carp "scatterBrains: an open database of human head models and companion optode locations for realistic Monte Carlo photon simulations," Journal of Biomedical Optics 28(10), 100501 (5 October 2023). https://doi.org/10.1117/1.JBO.28.10.100501
Received: 24 August 2023; Accepted: 23 September 2023; Published: 5 October 2023
Lens.org Logo
CITATIONS
Cited by 1 scholarly publication.
Advertisement
Advertisement
KEYWORDS
Head

Simulations

Data modeling

Databases

Tissues

Monte Carlo methods

Sensors

Back to Top