Improving model-based functional near-infrared spectroscopy analysis using mesh-based anatomical and light-transport models

Abstract. Significance: Functional near-infrared spectroscopy (fNIRS) has become an important research tool in studying human brains. Accurate quantification of brain activities via fNIRS relies upon solving computational models that simulate the transport of photons through complex anatomy. Aim: We aim to highlight the importance of accurate anatomical modeling in the context of fNIRS and propose a robust method for creating high-quality brain/full-head tetrahedral mesh models for neuroimaging analysis. Approach: We have developed a surface-based brain meshing pipeline that can produce significantly better brain mesh models, compared to conventional meshing techniques. It can convert segmented volumetric brain scans into multilayered surfaces and tetrahedral mesh models, with typical processing times of only a few minutes and broad utilities, such as in Monte Carlo or finite-element-based photon simulations for fNIRS studies. Results: A variety of high-quality brain mesh models have been successfully generated by processing publicly available brain atlases. In addition, we compare three brain anatomical models—the voxel-based brain segmentation, tetrahedral brain mesh, and layered-slab brain model—and demonstrate noticeable discrepancies in brain partial pathlengths when using approximated brain anatomies, ranging between −1.5  %   to 23% with the voxelated brain and 36% to 166% with the layered-slab brain. Conclusion: The generation and utility of high-quality brain meshes can lead to more accurate brain quantification in fNIRS studies. Our open-source meshing toolboxes “Brain2Mesh” and “Iso2Mesh” are freely available at http://mcx.space/brain2mesh.

[1]  J. Haueisen,et al.  Influence of head models on EEG simulations and inverse source localizations , 2006, Biomedical engineering online.

[2]  Isabelle Bloch,et al.  Realistic models of children heads from 3D-MRI segmentation and tetrahedral mesh construction , 2004, Proceedings. 2nd International Symposium on 3D Data Processing, Visualization and Transmission, 2004. 3DPVT 2004..

[3]  David R. Kaeli,et al.  Scalable and massively parallel Monte Carlo photon transport simulations for heterogeneous computing platforms , 2017, Journal of biomedical optics.

[4]  J. M. Cascón,et al.  An automatic strategy for adaptive tetrahedral mesh generation , 2009 .

[5]  John E. Richards,et al.  Age-specific MRI brain and head templates for healthy adults from 20 through 89 years of age , 2015, Front. Aging Neurosci..

[6]  Myung Yung Jeong,et al.  Cortical Signal Analysis and Advances in Functional Near-Infrared Spectroscopy Signal: A Review , 2016, Front. Hum. Neurosci..

[7]  Marco Attene,et al.  A lightweight approach to repairing digitized polygon meshes , 2010, The Visual Computer.

[8]  Bharat K. Soni,et al.  Handbook of Grid Generation , 1998 .

[9]  Michael Garland,et al.  Surface simplification using quadric error metrics , 1997, SIGGRAPH.

[10]  David A Boas,et al.  Comparison of a layered slab and an atlas head model for Monte Carlo fitting of time-domain near-infrared spectroscopy data of the adult head , 2014, Journal of biomedical optics.

[11]  Eric L. Miller,et al.  A Multi-Modality Image Reconstruction Platform for Diffuse Optical Tomography , 2008 .

[12]  D. Louis Collins,et al.  Twenty New Digital Brain Phantoms for Creation of Validation Image Data Bases , 2006, IEEE Transactions on Medical Imaging.

[13]  Ziji Wu,et al.  Multiple material marching cubes algorithm , 2003 .

[14]  Bernhard Preim,et al.  Comparison of Fundamental Mesh Smoothing Algorithms for Medical Surface Models , 2006, SimVis.

[15]  Valer Jurcak,et al.  10/20, 10/10, and 10/5 systems revisited: Their validity as relative head-surface-based positioning systems , 2007, NeuroImage.

[16]  Qianqian Fang,et al.  Quantitative assessment of diffuse optical tomography sensitivity to the cerebral cortex using a whole-head probe. , 2012, Physics in medicine and biology.

[17]  Bart Vanrumste,et al.  Review on solving the forward problem in EEG source analysis , 2007, Journal of NeuroEngineering and Rehabilitation.

[18]  David A. Boas,et al.  Validating atlas-guided DOT: A comparison of diffuse optical tomography informed by atlas and subject-specific anatomies , 2012, NeuroImage.

[19]  Xavier Tricoche,et al.  Influence of tissue conductivity anisotropy on EEG/MEG field and return current computation in a realistic head model: A simulation and visualization study using high-resolution finite element modeling , 2006, NeuroImage.

[20]  Pieter Verboven,et al.  Modeling the propagation of light in realistic tissue structures with MMC-fpf: a meshed Monte Carlo method with free phase function. , 2013, Optics express.

[21]  Qianqian Fang,et al.  Dual-grid mesh-based Monte Carlo algorithm for efficient photon transport simulations in complex three-dimensional media , 2019, Journal of biomedical optics.

[22]  Peter Lindstrom,et al.  Fast and memory efficient polygonal simplification , 1998, Proceedings Visualization '98 (Cat. No.98CB36276).

[23]  Arthur W. Toga,et al.  The generation of tetrahedral mesh models for neuroanatomical MRI , 2011, NeuroImage.

[24]  M. Fox,et al.  Connectivity Predicts deep brain stimulation outcome in Parkinson disease , 2017, Annals of neurology.

[25]  Ron Kikinis,et al.  Real-time registration of volumetric brain MRI by biomechanical simulation of deformation during image guided neurosurgery , 2002 .

[26]  Laurent Rineau,et al.  High-Quality Consistent Meshing of Multi-label Datasets , 2007, IPMI.

[27]  Yu Huang,et al.  Automated MRI segmentation for individualized modeling of current flow in the human head , 2013, Journal of neural engineering.

[28]  Steven J. Owen,et al.  A Survey of Unstructured Mesh Generation Technology , 1998, IMR.

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

[30]  Hamid Dehghani,et al.  Near infrared optical tomography using NIRFAST: Algorithm for numerical model and image reconstruction. , 2009, Communications in numerical methods in engineering.

[31]  Qianqian Fang,et al.  Selective photobiomodulation for emotion regulation: model-based dosimetry study , 2019, Neurophotonics.

[32]  Vladislav Toronov,et al.  Measurement of the optical properties of a two-layer model of the human head using broadband near-infrared spectroscopy. , 2010, Applied optics.

[33]  Marco Ferrari,et al.  A brief review on the history of human functional near-infrared spectroscopy (fNIRS) development and fields of application , 2012, NeuroImage.

[34]  Clifford C. Chou,et al.  Modeling of the Brain for Injury Prevention , 2011 .

[35]  Carmen E Sanchez,et al.  Age-Specific MRI Templates for Pediatric Neuroimaging , 2012, Developmental neuropsychology.

[36]  William E. Lorensen,et al.  Marching cubes: A high resolution 3D surface construction algorithm , 1987, SIGGRAPH.

[37]  D. Stegeman,et al.  Simulating Transcranial Direct Current Stimulation With a Detailed Anisotropic Human Head Model , 2014, IEEE Transactions on Neural Systems and Rehabilitation Engineering.

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

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

[40]  Solomon G. Diamond,et al.  Algorithm to find high density EEG scalp coordinates and analysis of their correspondence to structural and functional regions of the brain , 2014, Journal of Neuroscience Methods.

[41]  Simon R. Arridge,et al.  A 4D neonatal head model for diffuse optical imaging of pre-term to term infants , 2014, NeuroImage.

[42]  Lothar Lilge,et al.  High-performance, robustly verified Monte Carlo simulation with FullMonte , 2018, Journal of biomedical optics.

[43]  Richard M. Leahy,et al.  Brainstorm: A User-Friendly Application for MEG/EEG Analysis , 2011, Comput. Intell. Neurosci..

[44]  Daniel A. Rüfenacht,et al.  Light transport in tissue by 3D Monte Carlo: Influence of boundary voxelization , 2008, Comput. Methods Programs Biomed..

[45]  David A. Boas,et al.  Tetrahedral mesh generation from volumetric binary and grayscale images , 2009, 2009 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[46]  Paolo Cignoni,et al.  Metro: Measuring Error on Simplified Surfaces , 1998, Comput. Graph. Forum.

[47]  Ross T. Whitaker,et al.  Lattice Cleaving: Conforming Tetrahedral Meshes of Multimaterial Domains with Bounded Quality , 2012, IMR.

[48]  Dan Liu,et al.  Influence of extracerebral layers on estimates of optical properties with continuous wave near infrared spectroscopy: analysis based on multi-layered brain tissue architecture and Monte Carlo simulation , 2019, Computer assisted surgery.

[49]  B. Joe,et al.  Relationship between tetrahedron shape measures , 1994 .

[50]  D. Boas,et al.  Effective scattering coefficient of the cerebral spinal fluid in adult head models for diffuse optical imaging. , 2006, Applied optics.

[51]  Ashraf El-Hamalawi,et al.  Mesh Generation – Application to Finite Elements , 2001 .

[52]  Jason F. Shepherd,et al.  A meshing pipeline for biomedical computing , 2009, Engineering with Computers.

[53]  Hang Si,et al.  TetGen, a Delaunay-Based Quality Tetrahedral Mesh Generator , 2015, ACM Trans. Math. Softw..

[54]  Hamid Dehghani,et al.  Fast segmentation and high-quality three-dimensional volume mesh creation from medical images for diffuse optical tomography , 2013, Journal of biomedical optics.

[55]  Yu Shang,et al.  A Nth-order linear algorithm for extracting diffuse correlation spectroscopy blood flow indices in heterogeneous tissues. , 2014, Applied physics letters.

[56]  Alexander Opitz,et al.  Electric field calculations in brain stimulation based on finite elements: An optimized processing pipeline for the generation and usage of accurate individual head models , 2013, Human brain mapping.

[57]  Richard S. Frackowiak,et al.  New tissue priors for improved automated classification of subcortical brain structures on MRI☆ , 2016, NeuroImage.

[58]  Mariette Yvinec,et al.  Mesh Generation from 3D Multi-material Images , 2009, MICCAI.

[59]  Mahlega S. Hassanpour,et al.  Mapping distributed brain function and networks with diffuse optical tomography , 2014, Nature Photonics.

[60]  Steve Oudot,et al.  Provably good sampling and meshing of surfaces , 2005, Graph. Model..