1
|
Lewis AV, Fang Q. Revisiting equivalent optical properties for cerebrospinal fluid to improve diffusion-based modeling accuracy in the brain. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.08.20.608859. [PMID: 39229084 PMCID: PMC11370459 DOI: 10.1101/2024.08.20.608859] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Indexed: 09/05/2024]
Abstract
Significance The diffusion approximation (DA) is used in functional near-infrared spectroscopy (fNIRS) studies despite its known limitations due to the presence of cerebrospinal fluid (CSF). Nearly all of these studies rely on a set of empirical CSF optical properties, recommended by a previous simulation study, that were not selected for the purpose of minimizing DA modeling errors. Aim We aim to directly quantify the accuracy of DA solutions in brain models by comparing those with the goldstandard solutions produced by the mesh-based Monte Carlo (MMC), based on which we derive updated recommendations. Approach For both a 5-layer head and Colin27 atlas models, we obtain DA solutions by independently sweeping the CSF absorption ( μ a ) and reduced scattering coefficients. Using an MMC solution with literature CSF optical properties as reference, we compute the errors for surface fluence, total brain sensitivity and brain energy-deposition, and identify the optimized settings where the such error is minimized. Results Our results suggest that previously recommended CSF properties can cause significant errors (8.7% to 52%) in multiple tested metrics. By simultaneously sweeping μ a and , we can identify infinite numbers of solutions that can exactly match DA with MMC solutions for any single tested metric. Furthermore, it is also possible to simultaneously minimize multiple metrics at multiple source/detector separations, leading to our new recommendation of setting while maintaining physiological μ a for CSF in DA simulations. Conclusion Our new recommendation of CSF equivalent optical properties can greatly reduce the model mismatches between DA and MMC solutions at multiple metrics without sacrificing computational speed. We also show that it is possible to eliminate such a mismatch for a single or a pair of metrics of interest.
Collapse
|
2
|
Srinivasan S, Acharya D, Butters E, Collins-Jones L, Mancini F, Bale G. Subject-specific information enhances spatial accuracy of high-density diffuse optical tomography. FRONTIERS IN NEUROERGONOMICS 2024; 5:1283290. [PMID: 38444841 PMCID: PMC10910052 DOI: 10.3389/fnrgo.2024.1283290] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/25/2023] [Accepted: 02/02/2024] [Indexed: 03/07/2024]
Abstract
Functional near-infrared spectroscopy (fNIRS) is a widely used imaging method for mapping brain activation based on cerebral hemodynamics. The accurate quantification of cortical activation using fNIRS data is highly dependent on the ability to correctly localize the positions of light sources and photodetectors on the scalp surface. Variations in head size and shape across participants greatly impact the precise locations of these optodes and consequently, the regions of the cortical surface being reached. Such variations can therefore influence the conclusions drawn in NIRS studies that attempt to explore specific cortical regions. In order to preserve the spatial identity of each NIRS channel, subject-specific differences in NIRS array registration must be considered. Using high-density diffuse optical tomography (HD-DOT), we have demonstrated the inter-subject variability of the same HD-DOT array applied to ten participants recorded in the resting state. We have also compared three-dimensional image reconstruction results obtained using subject-specific positioning information to those obtained using generic optode locations. To mitigate the error introduced by using generic information for all participants, photogrammetry was used to identify specific optode locations per-participant. The present work demonstrates the large variation between subjects in terms of which cortical parcels are sampled by equivalent channels in the HD-DOT array. In particular, motor cortex recordings suffered from the largest optode localization errors, with a median localization error of 27.4 mm between generic and subject-specific optodes, leading to large differences in parcel sensitivity. These results illustrate the importance of collecting subject-specific optode locations for all wearable NIRS experiments, in order to perform accurate group-level analysis using cortical parcellation.
Collapse
Affiliation(s)
- Sruthi Srinivasan
- Department of Engineering, University of Cambridge, Cambridge, United Kingdom
| | - Deepshikha Acharya
- Department of Engineering, University of Cambridge, Cambridge, United Kingdom
| | - Emilia Butters
- Department of Engineering, University of Cambridge, Cambridge, United Kingdom
- Department of Psychiatry, University of Cambridge, Cambridge, United Kingdom
| | - Liam Collins-Jones
- Department of Medical Physics and Biomedical Engineering, University College London, London, United Kingdom
- Department of Clinical Neurosciences, University of Cambridge, Cambridge, United Kingdom
| | - Flavia Mancini
- Department of Engineering, University of Cambridge, Cambridge, United Kingdom
| | - Gemma Bale
- Department of Engineering, University of Cambridge, Cambridge, United Kingdom
- Department of Physics, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
3
|
Cockx H, Oostenveld R, Tabor M, Savenco E, van Setten A, Cameron I, van Wezel R. fNIRS is sensitive to leg activity in the primary motor cortex after systemic artifact correction. Neuroimage 2023; 269:119880. [PMID: 36693595 DOI: 10.1016/j.neuroimage.2023.119880] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2022] [Revised: 11/17/2022] [Accepted: 01/13/2023] [Indexed: 01/22/2023] Open
Abstract
BACKGROUND functional near-infrared spectroscopy (fNIRS) is an increasingly popular tool to study cortical activity during movement and gait that requires further validation. This study aimed to assess (1) whether fNIRS can detect the difficult-to-measure leg area of the primary motor cortex (M1) and distinguish it from the hand area; and (2) whether fNIRS can differentiate between automatic (i.e., not requiring one's attention) and non-automatic movement processes. Special attention was attributed to systemic artifacts (i.e., changes in blood pressure, heart rate, breathing) which were assessed and corrected by short channels, i.e., fNIRS channels which are mainly sensitive to superficial scalp hemodynamics. METHODS Twenty-three seated, healthy participants tapped four fingers on a keyboard or tapped the right foot on four squares on the floor in a specific order given by a 12-digit sequence (e.g., 434141243212). Two different sequences were executed: a beforehand learned (i.e., automatic) version and a newly learned (i.e., non-automatic) version. A 36-channel fNIRS device including 12 short channels covered multiple motor-related cortical areas including M1. The fNIRS data were analyzed with a general linear model (GLM). Correlation between the expected functional hemodynamic responses (i.e. task regressor) and the short channels (i.e. nuisance regressors), necessitated performing a separate short channel regression instead of integrating them in the GLM. RESULTS Consistent with the M1 somatotopy, we found significant HbO increases of very large effect size in the lateral M1 channels during finger tapping (Cohen's d = 1.35, p<0.001) and significant HbO increases of moderate effect size in the medial M1 channels during foot tapping (Cohen's d = 0.8, p<0.05). The cortical activity differences between automatic and non-automatic tasks were not significantly different. Importantly, leg movements produced large systemic fluctuations, which were adequately removed by the use of all available short channels. DISCUSSION Our results indicate that fNIRS is sensitive to leg activity in M1, though the sensitivity is lower than for finger activity and requires rigorous correction for systemic fluctuations. We furthermore highlight that systemic artifacts may result in an unreliable GLM analysis when short channels show signals that are similar to the expected hemodynamic responses.
Collapse
Affiliation(s)
- Helena Cockx
- Donders Institute for Brain, Cognition and Behaviour, Biophysics Department, Faculty of Science, Radboud University, Heyendaalseweg 135, 6525AJ Nijmegen, the Netherlands.
| | - Robert Oostenveld
- Donders Institute for Brain Cognition and Behaviour, Donders Center for Cognitive Neuroimaging, Radboud University, Kapittelweg 29, 6525EN Nijmegen, the Netherlands; NatMEG, Karolinska Institutet, Nobels Väg 9, D2:D235, 17177 Stockholm, Sweden.
| | - Merel Tabor
- Donders Institute for Brain, Cognition and Behaviour, Biophysics Department, Faculty of Science, Radboud University, Heyendaalseweg 135, 6525AJ Nijmegen, the Netherlands
| | - Ecaterina Savenco
- Donders Institute for Brain, Cognition and Behaviour, Biophysics Department, Faculty of Science, Radboud University, Heyendaalseweg 135, 6525AJ Nijmegen, the Netherlands
| | - Arne van Setten
- Donders Institute for Brain, Cognition and Behaviour, Biophysics Department, Faculty of Science, Radboud University, Heyendaalseweg 135, 6525AJ Nijmegen, the Netherlands
| | - Ian Cameron
- Donders Institute for Brain, Cognition and Behaviour, Biophysics Department, Faculty of Science, Radboud University, Heyendaalseweg 135, 6525AJ Nijmegen, the Netherlands; OnePlanet Research Center, Toernooiveld 300, 6525EC Nijmegen, the Netherlands; Biomedical Signals and Systems Group, Faculty of Electrical Engineering, Mathematics and Computer Science (EEMCS), University of Twente, Drienerlolaan 5, 7522NB Enschede, the Netherlands.
| | - Richard van Wezel
- Donders Institute for Brain, Cognition and Behaviour, Biophysics Department, Faculty of Science, Radboud University, Heyendaalseweg 135, 6525AJ Nijmegen, the Netherlands; Biomedical Signals and Systems Group, Faculty of Electrical Engineering, Mathematics and Computer Science (EEMCS), University of Twente, Drienerlolaan 5, 7522NB Enschede, the Netherlands.
| |
Collapse
|
4
|
Cai Z, Machado A, Chowdhury RA, Spilkin A, Vincent T, Aydin Ü, Pellegrino G, Lina JM, Grova C. Diffuse optical reconstructions of functional near infrared spectroscopy data using maximum entropy on the mean. Sci Rep 2022; 12:2316. [PMID: 35145148 PMCID: PMC8831678 DOI: 10.1038/s41598-022-06082-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2021] [Accepted: 01/24/2022] [Indexed: 02/07/2023] Open
Abstract
Functional near-infrared spectroscopy (fNIRS) measures the hemoglobin concentration changes associated with neuronal activity. Diffuse optical tomography (DOT) consists of reconstructing the optical density changes measured from scalp channels to the oxy-/deoxy-hemoglobin concentration changes within the cortical regions. In the present study, we adapted a nonlinear source localization method developed and validated in the context of Electro- and Magneto-Encephalography (EEG/MEG): the Maximum Entropy on the Mean (MEM), to solve the inverse problem of DOT reconstruction. We first introduced depth weighting strategy within the MEM framework for DOT reconstruction to avoid biasing the reconstruction results of DOT towards superficial regions. We also proposed a new initialization of the MEM model improving the temporal accuracy of the original MEM framework. To evaluate MEM performance and compare with widely used depth weighted Minimum Norm Estimate (MNE) inverse solution, we applied a realistic simulation scheme which contained 4000 simulations generated by 250 different seeds at different locations and 4 spatial extents ranging from 3 to 40[Formula: see text] along the cortical surface. Our results showed that overall MEM provided more accurate DOT reconstructions than MNE. Moreover, we found that MEM was remained particularly robust in low signal-to-noise ratio (SNR) conditions. The proposed method was further illustrated by comparing to functional Magnetic Resonance Imaging (fMRI) activation maps, on real data involving finger tapping tasks with two different montages. The results showed that MEM provided more accurate HbO and HbR reconstructions in spatial agreement with the main fMRI cluster, when compared to MNE.
Collapse
Affiliation(s)
- Zhengchen Cai
- Department of Physics and PERFORM Centre, Concordia University, Montreal, Canada.
| | - Alexis Machado
- Multimodal Functional Imaging Lab, Biomedical Engineering Department, McGill University, Montreal, Canada
| | - Rasheda Arman Chowdhury
- Multimodal Functional Imaging Lab, Biomedical Engineering Department, McGill University, Montreal, Canada
| | - Amanda Spilkin
- Department of Physics and PERFORM Centre, Concordia University, Montreal, Canada
| | - Thomas Vincent
- Department of Physics and PERFORM Centre, Concordia University, Montreal, Canada
- Neurology and Neurosurgery Department, Montreal Neurological Institute, McGill University, Montreal, Canada
- Centre de médecine préventive et d'activité physique, Montréal Heart Institute, Montréal, Canada
| | - Ümit Aydin
- Department of Physics and PERFORM Centre, Concordia University, Montreal, Canada
- MRC Social, Genetic and Developmental Psychiatry Centre, Institute of Psychiatry, Psychology and Neuroscience, King's College London, London, UK
| | - Giovanni Pellegrino
- Neurology and Neurosurgery Department, Montreal Neurological Institute, McGill University, Montreal, Canada
| | - Jean-Marc Lina
- École de technologie supérieure de l'Université du Québec, Montréal, Canada
- Centre de Recherches Mathématiques, Université de Montréal, Montréal, Canada
| | - Christophe Grova
- Department of Physics and PERFORM Centre, Concordia University, Montreal, Canada
- Multimodal Functional Imaging Lab, Biomedical Engineering Department, McGill University, Montreal, Canada
- Neurology and Neurosurgery Department, Montreal Neurological Institute, McGill University, Montreal, Canada
- Centre de Recherches Mathématiques, Université de Montréal, Montréal, Canada
| |
Collapse
|
5
|
Cai Z, Uji M, Aydin Ü, Pellegrino G, Spilkin A, Delaire É, Abdallah C, Lina J, Grova C. Evaluation of a personalized functional near infra-red optical tomography workflow using maximum entropy on the mean. Hum Brain Mapp 2021; 42:4823-4843. [PMID: 34342073 PMCID: PMC8449120 DOI: 10.1002/hbm.25566] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2021] [Revised: 05/25/2021] [Accepted: 06/11/2021] [Indexed: 11/20/2022] Open
Abstract
In the present study, we proposed and evaluated a workflow of personalized near infra-red optical tomography (NIROT) using functional near-infrared spectroscopy (fNIRS) for spatiotemporal imaging of cortical hemodynamic fluctuations. The proposed workflow from fNIRS data acquisition to local 3D reconstruction consists of: (a) the personalized optimal montage maximizing fNIRS channel sensitivity to a predefined targeted brain region; (b) the optimized fNIRS data acquisition involving installation of optodes and digitalization of their positions using a neuronavigation system; and (c) the 3D local reconstruction using maximum entropy on the mean (MEM) to accurately estimate the location and spatial extent of fNIRS hemodynamic fluctuations along the cortical surface. The workflow was evaluated on finger-tapping fNIRS data acquired from 10 healthy subjects for whom we estimated the reconstructed NIROT spatiotemporal images and compared with functional magnetic resonance imaging (fMRI) results from the same individuals. Using the fMRI activation maps as our reference, we quantitatively compared the performance of two NIROT approaches, the MEM framework and the conventional minimum norm estimation (MNE) method. Quantitative comparisons were performed at both single subject and group-level. Overall, our results suggested that MEM provided better spatial accuracy than MNE, while both methods offered similar temporal accuracy when reconstructing oxygenated (HbO) and deoxygenated hemoglobin (HbR) concentration changes evoked by finger-tapping. Our proposed complete workflow was made available in the brainstorm fNIRS processing plugin-NIRSTORM, thus providing the opportunity for other researchers to further apply it to other tasks and on larger populations.
Collapse
Affiliation(s)
- Zhengchen Cai
- Multimodal Functional Imaging Lab, Department of Physics and PERFORM CentreConcordia UniversityMontréalQuébecCanada
| | - Makoto Uji
- Multimodal Functional Imaging Lab, Department of Physics and PERFORM CentreConcordia UniversityMontréalQuébecCanada
| | - Ümit Aydin
- Multimodal Functional Imaging Lab, Department of Physics and PERFORM CentreConcordia UniversityMontréalQuébecCanada
- Social, Genetic and Developmental Psychiatry CentreInstitute of Psychiatry, Psychology and Neuroscience, King's College LondonLondonUK
| | - Giovanni Pellegrino
- Neurology and Neurosurgery Department, Montreal Neurological InstituteMcGill UniversityMontréalQuébecCanada
- Multimodal Functional Imaging Lab, Biomedical Engineering DepartmentMcGill UniversityMontréalQuébecCanada
| | - Amanda Spilkin
- Multimodal Functional Imaging Lab, Department of Physics and PERFORM CentreConcordia UniversityMontréalQuébecCanada
| | - Édouard Delaire
- Multimodal Functional Imaging Lab, Department of Physics and PERFORM CentreConcordia UniversityMontréalQuébecCanada
| | - Chifaou Abdallah
- Neurology and Neurosurgery Department, Montreal Neurological InstituteMcGill UniversityMontréalQuébecCanada
- Multimodal Functional Imaging Lab, Biomedical Engineering DepartmentMcGill UniversityMontréalQuébecCanada
| | - Jean‐Marc Lina
- Département de Génie ElectriqueÉcole de Technologie SupérieureMontréalQuébecCanada
- Centre De Recherches En MathématiquesMontréalQuébecCanada
| | - Christophe Grova
- Multimodal Functional Imaging Lab, Department of Physics and PERFORM CentreConcordia UniversityMontréalQuébecCanada
- Neurology and Neurosurgery Department, Montreal Neurological InstituteMcGill UniversityMontréalQuébecCanada
- Multimodal Functional Imaging Lab, Biomedical Engineering DepartmentMcGill UniversityMontréalQuébecCanada
- Centre De Recherches En MathématiquesMontréalQuébecCanada
| |
Collapse
|
6
|
Forbes SH, Wijeakumar S, Eggebrecht AT, Magnotta VA, Spencer JP. Processing pipeline for image reconstructed fNIRS analysis using both MRI templates and individual anatomy. NEUROPHOTONICS 2021; 8:025010. [PMID: 35106319 PMCID: PMC8786393 DOI: 10.1117/1.nph.8.2.025010] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/14/2021] [Accepted: 05/18/2021] [Indexed: 05/29/2023]
Abstract
Significance: Image reconstruction of fNIRS data is a useful technique for transforming channel-based fNIRS into a volumetric representation and managing spatial variance based on optode location. We present an innovative integrated pipeline for image reconstruction of fNIRS data using either MRI templates or individual anatomy. Aim: We demonstrate a pipeline with accompanying code to allow users to clean and prepare optode location information, prepare and standardize individual anatomical images, create the light model, run the 3D image reconstruction, and analyze data in group space. Approach: We synthesize a combination of new and existing software packages to create a complete pipeline, from raw data to analysis. Results: This pipeline has been tested using both templates and individual anatomy, and on data from different fNIRS data collection systems. We show high temporal correlations between channel-based and image-based fNIRS data. In addition, we demonstrate the reliability of this pipeline with a sample dataset that included 74 children as part of a longitudinal study taking place in Scotland. We demonstrate good correspondence between data in channel space and image reconstructed data. Conclusions: The pipeline presented here makes a unique contribution by integrating multiple tools to assemble a complete pipeline for image reconstruction in fNIRS. We highlight further issues that may be of interest to future software developers in the field.
Collapse
Affiliation(s)
- Samuel H. Forbes
- University of East Anglia, School of Psychology, Lawrence Stenhouse Building, Norwich, United Kingdom
| | | | - Adam T. Eggebrecht
- Washington University, Mallinckrodt Institute of Radiology, St Louis, Missouri, United States
| | | | - John P. Spencer
- University of East Anglia, School of Psychology, Lawrence Stenhouse Building, Norwich, United Kingdom
| |
Collapse
|
7
|
Yücel MA, Lühmann AV, Scholkmann F, Gervain J, Dan I, Ayaz H, Boas D, Cooper RJ, Culver J, Elwell CE, Eggebrecht A, Franceschini MA, Grova C, Homae F, Lesage F, Obrig H, Tachtsidis I, Tak S, Tong Y, Torricelli A, Wabnitz H, Wolf M. Best practices for fNIRS publications. NEUROPHOTONICS 2021; 8:012101. [PMID: 33442557 PMCID: PMC7793571 DOI: 10.1117/1.nph.8.1.012101] [Citation(s) in RCA: 126] [Impact Index Per Article: 42.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Accepted: 12/02/2020] [Indexed: 05/09/2023]
Abstract
The application of functional near-infrared spectroscopy (fNIRS) in the neurosciences has been expanding over the last 40 years. Today, it is addressing a wide range of applications within different populations and utilizes a great variety of experimental paradigms. With the rapid growth and the diversification of research methods, some inconsistencies are appearing in the way in which methods are presented, which can make the interpretation and replication of studies unnecessarily challenging. The Society for Functional Near-Infrared Spectroscopy has thus been motivated to organize a representative (but not exhaustive) group of leaders in the field to build a consensus on the best practices for describing the methods utilized in fNIRS studies. Our paper has been designed to provide guidelines to help enhance the reliability, repeatability, and traceability of reported fNIRS studies and encourage best practices throughout the community. A checklist is provided to guide authors in the preparation of their manuscripts and to assist reviewers when evaluating fNIRS papers.
Collapse
Affiliation(s)
- Meryem A. Yücel
- Boston University, Neurophotonics Center, Biomedical Engineering, Boston, Massachusetts, United States
- Massachusetts General Hospital, Harvard Medical School, MGH/HST Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Charlestown, Massachusetts, United States
| | - Alexander v. Lühmann
- Boston University, Neurophotonics Center, Biomedical Engineering, Boston, Massachusetts, United States
- Massachusetts General Hospital, Harvard Medical School, MGH/HST Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Charlestown, Massachusetts, United States
| | - Felix Scholkmann
- University Hospital Zurich, University of Zurich, Department of Neonatology, Biomedical Optics Research Laboratory, Neonatology Research, Zurich, Switzerland
- University of Bern, Institute for Complementary and Integrative Medicine, Bern, Switzerland
| | - Judit Gervain
- Université de Paris, CNRS, Integrative Neuroscience and Cognition Center, Paris, France
- Università di Padova, Department of Social and Developmental Psychology, Padua, Italy
| | - Ippeita Dan
- Chuo University, Faculty of Science and Engineering, Applied Cognitive Neuroscience Laboratory, Tokyo, Japan
| | - Hasan Ayaz
- Drexel University, School of Biomedical Engineering, Science and Health Systems, Philadelphia, Pennsylvania, United States
- Drexel University, College of Arts and Sciences, Department of Psychology, Philadelphia, Pennsylvania, United States
- Drexel University, Drexel Solutions Institute, Philadelphia, Pennsylvania, United States
- University of Pennsylvania, Department of Family and Community Health, Philadelphia, Pennsylvania, United States
- Children’s Hospital of Philadelphia, Center for Injury Research and Prevention, Philadelphia, Pennsylvania, United States
| | - David Boas
- Boston University, Neurophotonics Center, Biomedical Engineering, Boston, Massachusetts, United States
| | - Robert J. Cooper
- University College London, DOT-HUB, Department of Medical Physics and Biomedical Engineering, Biomedical Optics Research Laboratory, London, United Kingdom
| | - Joseph Culver
- Washington University School of Medicine, Department of Radiology, St. Louis, Missouri, United States
| | - Clare E. Elwell
- University College London, Department of Medical Physics and Biomedical Engineering, London, United Kingdom
| | - Adam Eggebrecht
- Washington University School of Medicine, Mallinckrodt Institute of Radiology, St. Louis, Missouri, United States
| | - Maria A. Franceschini
- Massachusetts General Hospital, Harvard Medical School, MGH/HST Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Charlestown, Massachusetts, United States
| | - Christophe Grova
- Concordia University, Department of Physics and PERFORM Centre, Multimodal Functional Imaging Lab, Montreal, Québec, Canada
- McGill University, Biomedical Engineering Department, Multimodal Functional Imaging Lab, Montreal, Québec, Canada
| | - Fumitaka Homae
- Tokyo Metropolitan University, Department of Language Sciences, Tokyo, Japan
| | - Frédéric Lesage
- Polytechnique Montréal, Department Electrical Engineering, Montreal, Canada
| | - Hellmuth Obrig
- University Hospital Leipzig, Max-Planck-Institute for Human Cognitive and Brain Sciences and Clinic for Cognitive Neurology, Leipzig, Germany
| | - Ilias Tachtsidis
- University College London, Department of Medical Physics and Biomedical Engineering, London, United Kingdom
| | - Sungho Tak
- Korea Basic Science Institute, Research Center for Bioconvergence Analysis, Ochang, Cheongju, Republic of Korea
| | - Yunjie Tong
- Weldon School of Biomedical Engineering Purdue University, West Lafayette, Indiana, United States
| | - Alessandro Torricelli
- Politecnico di Milano, Dipartimento di Fisica, Milan, Italy
- Consiglio Nazionale delle Ricerche, Istituto di Fotonica e Nanotecnologie, Milan, Italy
| | | | - Martin Wolf
- University Hospital Zurich, University of Zurich, Department of Neonatology, Biomedical Optics Research Laboratory, Neonatology Research, Zurich, Switzerland
| |
Collapse
|
8
|
Ebrahimpour A, Zakariaee SS, Hejazi M. Sensitivity Laplacian Ratio-Based Optimization of the Projection Selection for Diffuse Optical Tomography. JOURNAL OF MEDICAL SIGNALS AND SENSORS 2020; 10:119-124. [PMID: 32676448 PMCID: PMC7359954 DOI: 10.4103/jmss.jmss_67_18] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/25/2018] [Revised: 10/28/2019] [Accepted: 01/06/2020] [Indexed: 11/17/2022]
Abstract
BACKGROUND In diffuse optical tomography, determining the optimal angle between the source and detector is an effective method to reduce the number of projections while maintaining the quality of the reconstructed images. In this study, a new parameter is introduced to evaluate the source-detector geometries. METHODS A two-dimensional mesh with the radius of 20 mm and 7987 nodes were built. In each reconstruction, 0.5 mm heterogeneity with the absorption coefficient of 0.06 mm-1 and the dispersion coefficient of 0.6 mm-1 was added in different parts of the sample randomly. The relationship between the mean square error (MSE), sensitivity Laplacian ratio (SLR), and sensitivity standard deviation ratio (SSR) was evaluated based on their correlation coefficients. The quality of the images achieved using the optimized projections were compared with that of the full projections for the same depths. RESULTS MSE decreases by increasing the SLR magnitudes which indicate that the parameter could be used to evaluate the scanning geometries. There was a negative correlation coefficient (R = -0.76) with the inverse relationship between the SLR and MSE indices. SSR does not have a significant relationship with the quality of the reconstructed images. For each scanning depth, the comparison of the images obtained using the full and optimized-selective projections did not show any considerable difference despite the decrease of the projection numbers in scanning geometry with the optimized-selective projections. CONCLUSION The unnecessary projections could be eliminated by placing the detectors at the specific angles, which were determined using the SLR. Thus, a proper compromise between the quality of the reconstructed images and reconstruction time might establish.
Collapse
Affiliation(s)
- Anita Ebrahimpour
- Department of Medical Physics and Biomedical Engineering, Faculty of Medicine, Tehran University of Medical Sciences, Tehran, Iran
| | - Seyed Salman Zakariaee
- Department of Medical Physics, Faculty of Paramedical Sciences, Ilam University of Medical Sciences, Ilam, Iran
| | - Marjaneh Hejazi
- Department of Medical Physics and Biomedical Engineering, Faculty of Medicine, Tehran University of Medical Sciences, Tehran, Iran,Address for correspondence: Dr. Marjaneh Hejazi, Department of Medical Physics and Biomedical Engineering, Faculty of Medicine, Tehran University of Medical Sciences, Tehran, Iran. E-mail:
| |
Collapse
|
9
|
Dragojević T, Vidal Rosas EE, Hollmann JL, Culver JP, Justicia C, Durduran T. High-density speckle contrast optical tomography of cerebral blood flow response to functional stimuli in the rodent brain. NEUROPHOTONICS 2019; 6:045001. [PMID: 31620545 PMCID: PMC6782685 DOI: 10.1117/1.nph.6.4.045001] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/08/2019] [Accepted: 09/12/2019] [Indexed: 05/20/2023]
Abstract
Noninvasive, three-dimensional, and longitudinal imaging of cerebral blood flow (CBF) in small animal models and ultimately in humans has implications for fundamental research and clinical applications. It enables the study of phenomena such as brain development and learning and the effects of pathologies, with a clear vision for translation to humans. Speckle contrast optical tomography (SCOT) is an emerging optical method that aims to achieve this goal by directly measuring three-dimensional blood flow maps in deep tissue with a relatively inexpensive and simple system. High-density SCOT is developed to follow CBF changes in response to somatosensory cortex stimulation. Measurements are carried out through the intact skull on the rat brain. SCOT is able to follow individual trials in each brain hemisphere, where signal averaging resulted in comparable, cortical images to those of functional magnetic resonance images in spatial extent, location, and depth. Sham stimuli are utilized to demonstrate that the observed response is indeed due to local changes in the brain induced by forepaw stimulation. In developing and demonstrating the method, algorithms and analysis methods are developed. The results pave the way for longitudinal, nondestructive imaging in preclinical rodent models that can readily be translated to the human brain.
Collapse
Affiliation(s)
- Tanja Dragojević
- Institut de Ciències Fotòniques, The Barcelona Institute of Science and Technology, Barcelona, Spain
- Address all correspondence to Tanja Dragojević, E-mail:
| | - Ernesto E. Vidal Rosas
- Institut de Ciències Fotòniques, The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Joseph L. Hollmann
- Institut de Ciències Fotòniques, The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Joseph P. Culver
- Washington University, School of Medicine, Department of Radiology, St. Louis, Missouri, United States
- Washington University, Department of Physics, St. Louis, Missouri, United States
| | - Carles Justicia
- Institut d’Investigacions Biomèdiques de Barcelona, Consejo Superior de Investigaciones Científicas, Department of Brain Ischemia and Neurodegeneration, Barcelona, Spain
- Institut d’Investigacions Biomèdiques August Pi i Sunyer, Àrea de Neurociències, Barcelona, Spain
| | - Turgut Durduran
- Institut de Ciències Fotòniques, The Barcelona Institute of Science and Technology, Barcelona, Spain
- Institució Catalana de Recerca i Estudis Avançats, Barcelona, Spain
| |
Collapse
|
10
|
Daly MJ, Wilson BC, Irish JC, Jaffray DA. Navigated non-contact fluorescence tomography. ACTA ACUST UNITED AC 2019; 64:135021. [DOI: 10.1088/1361-6560/ab1f33] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
|
11
|
Doulgerakis M, Eggebrecht AT, Dehghani H. High-density functional diffuse optical tomography based on frequency-domain measurements improves image quality and spatial resolution. NEUROPHOTONICS 2019; 6:035007. [PMID: 31482102 PMCID: PMC6702521 DOI: 10.1117/1.nph.6.3.035007] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/09/2019] [Accepted: 07/30/2019] [Indexed: 05/18/2023]
Abstract
Measurements of dynamic near-infrared (NIR) light attenuation across the human head together with model-based image reconstruction algorithms allow the recovery of three-dimensional spatial brain activation maps. Previous studies using high-density diffuse optical tomography (HD-DOT) systems have reported improved image quality over sparse arrays. These HD-DOT systems incorporated multidistance overlapping continuous wave measurements that only recover differential intensity attenuation. We investigate the potential improvement in reconstructed image quality due to the additional incorporation of phase shift measurements, which reflect the time-of-flight of the measured NIR light, within the tomographic reconstruction from high-density measurements. To evaluate image reconstruction with and without the additional phase information, we simulated point spread functions across a whole-scalp field of view in 24 subject-specific anatomical models using an experimentally derived noise model. The addition of phase information improves the image quality by reducing localization error by up to 59% and effective resolution by up to 21% as compared to using the intensity attenuation measurements alone. Furthermore, we demonstrate that the phase data enable images to be resolved at deeper brain regions where intensity data fail, which is further supported by utilizing experimental data from a single subject measurement during a retinotopic experiment.
Collapse
Affiliation(s)
- Matthaios Doulgerakis
- University of Birmingham, School of Computer Science, Birmingham, England, United Kingdom
| | - Adam T. Eggebrecht
- Washington University School of Medicine, Mallinckrodt Institute of Radiology, St. Louis, Missouri, United States
| | - Hamid Dehghani
- University of Birmingham, School of Computer Science, Birmingham, England, United Kingdom
- Address all correspondence to Hamid Dehghani, E-mail:
| |
Collapse
|
12
|
Wheelock MD, Culver JP, Eggebrecht AT. High-density diffuse optical tomography for imaging human brain function. THE REVIEW OF SCIENTIFIC INSTRUMENTS 2019; 90:051101. [PMID: 31153254 PMCID: PMC6533110 DOI: 10.1063/1.5086809] [Citation(s) in RCA: 64] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/23/2018] [Accepted: 04/14/2019] [Indexed: 05/08/2023]
Abstract
This review describes the unique opportunities and challenges for noninvasive optical mapping of human brain function. Diffuse optical methods offer safe, portable, and radiation free alternatives to traditional technologies like positron emission tomography or functional magnetic resonance imaging (fMRI). Recent developments in high-density diffuse optical tomography (HD-DOT) have demonstrated capabilities for mapping human cortical brain function over an extended field of view with image quality approaching that of fMRI. In this review, we cover fundamental principles of the diffusion of near infrared light in biological tissue. We discuss the challenges involved in the HD-DOT system design and implementation that must be overcome to acquire the signal-to-noise necessary to measure and locate brain function at the depth of the cortex. We discuss strategies for validation of the sensitivity, specificity, and reliability of HD-DOT acquired maps of cortical brain function. We then provide a brief overview of some clinical applications of HD-DOT. Though diffuse optical measurements of neurophysiology have existed for several decades, tremendous opportunity remains to advance optical imaging of brain function to address a crucial niche in basic and clinical neuroscience: that of bedside and minimally constrained high fidelity imaging of brain function.
Collapse
Affiliation(s)
- Muriah D. Wheelock
- Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri 63110, USA
| | | | - Adam T. Eggebrecht
- Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri 63110, USA
| |
Collapse
|
13
|
Wang B, Pan T, Zhang Y, Liu D, Jiang J, Zhao H, Gao F. A Kalman-based tomographic scheme for directly reconstructing activation levels of brain function. OPTICS EXPRESS 2019; 27:3229-3246. [PMID: 30732347 DOI: 10.1364/oe.27.003229] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
In functional near-infrared spectroscopy (fNIRS), the conventional indirect approaches first separately recover the spatial distribution of the changes in the optical properties at every time point, and then extract the activation levels by a time-course analysis process at every site. In the tomographic implementation of fNIRS, i.e., diffuse optical tomography (DOT), these approaches not only suffer from the ill-posedness of the optical inversions and error propagation between the two successive steps, but also fail to achieve satisfactory temporal resolution due to the requirement for a complete data set. To cope with the above adversities of the indirect approaches, we propose herein a direct approach to tomographically reconstructing the activation levels by incorporating a Kalman scheme. Dynamic simulative and phantom experiments were conducted for the performance validation of the proposed approach, demonstrating its potentials to improve the calculated images and to relax the speed limitation of the instruments.
Collapse
|
14
|
Machado A, Cai Z, Pellegrino G, Marcotte O, Vincent T, Lina JM, Kobayashi E, Grova C. Optimal positioning of optodes on the scalp for personalized functional near-infrared spectroscopy investigations. J Neurosci Methods 2018; 309:91-108. [PMID: 30107210 DOI: 10.1016/j.jneumeth.2018.08.006] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2017] [Revised: 08/07/2018] [Accepted: 08/08/2018] [Indexed: 12/21/2022]
Abstract
BACKGROUND Application of functional Near InfraRed Spectroscopy (fNIRS) in neurology is still limited as a good optical coupling and optimized optode coverage of specific brain regions remains challenging, notably for prolonged monitoring. METHODS We propose to evaluate a new procedure allowing accurate investigation of specific brain regions. The procedure consists in: (i) A priori maximization of spatial sensitivity of fNIRS measurements targeting specific brain regions, while reducing the number of applied optodes in order to decrease installation time and improve subject comfort. (ii) Utilization of a 3D neuronavigation device and usage of collodion to glue optodes on the scalp, ensuring good optical contact for prolonged investigations. (iii) Local reconstruction of the hemodynamic activity along the cortical surface using inverse modelling. RESULTS Using realistic simulations, we demonstrated that maps derived from optimal montage acquisitions showed, after reconstruction, spatial resolution only slightly lower to that of ultra high density montages while significantly reducing the number of optodes. The optimal montages provided overall good quantitative accuracy especially at the peak of the spatially reconstructed map. We also evaluated real motor responses in two healthy subjects and obtained reproducible motor responses over different sessions. COMPARISON WITH EXISTING METHODS We are among the first to propose a mathematical optimization strategy, allowing high sensitivity measurements. CONCLUSIONS Our results support that using personalized optimal montages should allow to conduct accurate fNIRS studies in clinical settings and realistic lifestyle conditions.
Collapse
Affiliation(s)
- A Machado
- Multimodal Functional Imaging Laboratory, Biomedical Engineering Department, McGill University, Canada.
| | - Z Cai
- Physics Department and PERFORM center, Concordia University, Montreal, Canada
| | - G Pellegrino
- Multimodal Functional Imaging Laboratory, Biomedical Engineering Department, McGill University, Canada; IRCCS Fondazione Ospedale San Camillo Via Alberoni, Venice, Italy
| | - O Marcotte
- GERAD, École des HEC, Montréal, Canada; Département d'informatique, Université du Québec à Montréal, Canada; Centre de Recherches Mathématiques, Université de Montréal, Québec, Canada
| | - T Vincent
- Physics Department and PERFORM center, Concordia University, Montreal, Canada
| | - J-M Lina
- École de technologie supérieure de l'Université du Québec, Canada; Centre de Recherches Mathématiques, Université de Montréal, Québec, Canada
| | - E Kobayashi
- Montreal Neurological Institute, Department of Neurology and Neurosurgery, McGill University, Canada
| | - C Grova
- Multimodal Functional Imaging Laboratory, Biomedical Engineering Department, McGill University, Canada; Physics Department and PERFORM center, Concordia University, Montreal, Canada; Montreal Neurological Institute, Department of Neurology and Neurosurgery, McGill University, Canada; Centre de Recherches Mathématiques, Université de Montréal, Québec, Canada
| |
Collapse
|
15
|
Doulgerakis M, Eggebrecht AT, Wojtkiewicz S, Culver JP, Dehghani H. Toward real-time diffuse optical tomography: accelerating light propagation modeling employing parallel computing on GPU and CPU. JOURNAL OF BIOMEDICAL OPTICS 2017; 22:1-11. [PMID: 29197176 PMCID: PMC5709934 DOI: 10.1117/1.jbo.22.12.125001] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/28/2017] [Accepted: 11/06/2017] [Indexed: 05/18/2023]
Abstract
Parameter recovery in diffuse optical tomography is a computationally expensive algorithm, especially when used for large and complex volumes, as in the case of human brain functional imaging. The modeling of light propagation, also known as the forward problem, is the computational bottleneck of the recovery algorithm, whereby the lack of a real-time solution is impeding practical and clinical applications. The objective of this work is the acceleration of the forward model, within a diffusion approximation-based finite-element modeling framework, employing parallelization to expedite the calculation of light propagation in realistic adult head models. The proposed methodology is applicable for modeling both continuous wave and frequency-domain systems with the results demonstrating a 10-fold speed increase when GPU architectures are available, while maintaining high accuracy. It is shown that, for a very high-resolution finite-element model of the adult human head with ∼600,000 nodes, consisting of heterogeneous layers, light propagation can be calculated at ∼0.25 s/excitation source.
Collapse
Affiliation(s)
- Matthaios Doulgerakis
- University of Birmingham, School of Computer Science, Birmingham, United Kingdom
- Address all correspondence to: Matthaios Doulgerakis, E-mail:
| | - Adam T. Eggebrecht
- Washington University School of Medicine, Department of Radiology, St. Louis, Missouri, United States
| | | | - Joseph P. Culver
- Washington University School of Medicine, Department of Radiology, St. Louis, Missouri, United States
- Washington University in St. Louis, Department of Biomedical Engineering, St. Louis, Missouri, United States
- Washington University School of Medicine, Division of Biology and Biomedical Sciences, St. Louis, Missouri, United States
| | - Hamid Dehghani
- University of Birmingham, School of Computer Science, Birmingham, United Kingdom
| |
Collapse
|
16
|
Chiarelli AM, Maclin EL, Low KA, Fantini S, Fabiani M, Gratton G. Low-resolution mapping of the effective attenuation coefficient of the human head: a multidistance approach applied to high-density optical recordings. NEUROPHOTONICS 2017; 4:021103. [PMID: 28466026 PMCID: PMC5400126 DOI: 10.1117/1.nph.4.2.021103] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/03/2017] [Accepted: 03/30/2017] [Indexed: 05/29/2023]
Abstract
Near infrared (NIR) light has been widely used for measuring changes in hemoglobin concentration in the human brain (functional NIR spectroscopy, fNIRS). fNIRS is based on the differential measurement and estimation of absorption perturbations, which, in turn, are based on correctly estimating the absolute parameters of light propagation. To do so, it is essential to accurately characterize the baseline optical properties of tissue (absorption and reduced scattering coefficients). However, because of the diffusive properties of the medium, separate determination of absorption and scattering across the head is challenging. The effective attenuation coefficient (EAC), which is proportional to the geometric mean of absorption and reduced scattering coefficients, can be estimated in a simpler fashion by multidistance light decay measurements. EAC mapping could be of interest for the scientific community because of its absolute information content, and because light propagation is governed by the EAC for source-detector distances exceeding 1 cm, which sense depths extending beyond the scalp and skull layers. Here, we report an EAC mapping procedure that can be applied to standard fNIRS recordings, yielding topographic maps with 2- to 3-cm resolution. Application to human data indicates the importance of venous sinuses in determining regional EAC variations, a factor often overlooked.
Collapse
Affiliation(s)
| | - Edward L. Maclin
- University of Illinois, Beckman Institute, Urbana, Illinois, United States
| | - Kathy A. Low
- University of Illinois, Beckman Institute, Urbana, Illinois, United States
| | - Sergio Fantini
- Tufts University, Department of Biomedical Engineering, Medford, Massachusetts, United States
| | - Monica Fabiani
- University of Illinois, Beckman Institute, Urbana, Illinois, United States
- University of Illinois, Department of Psychology, Champaign, Illinois, United States
| | - Gabriele Gratton
- University of Illinois, Beckman Institute, Urbana, Illinois, United States
- University of Illinois, Department of Psychology, Champaign, Illinois, United States
| |
Collapse
|
17
|
Abstract
In this review I introduce the historical context and methods of optical neuroimaging, leading to the modern use of functional near-infrared spectroscopy (fNIRS) and high-density diffuse optical tomography (HD-DOT) to study human brain function. In its most frequent application, optical neuroimaging measures a hemodynamically-mediated signal indirectly related to neural processing, similar to that captured by fMRI. Compared to other approaches to measuring human brain function, optical imaging has many advantages: it is noninvasive, frequently portable, acoustically silent, robust to motion and muscle movement, and appropriate in many situations in which fMRI is not possible (for example, due to implanted medical devices). Challenges include producing a full-brain field of view, homogenous spatial resolution, and accurate source localization. Experimentally, optical neuroimaging has been used to study phoneme, word, and sentence processing in a variety of paradigms. With continuing technical and methodological improvements the future of optical neuroimaging is increasingly bright.
Collapse
Affiliation(s)
- Jonathan E Peelle
- Department of Otolaryngology, Washington University in St. Louis, St. Louis MO USA
| |
Collapse
|
18
|
Li L, Cazzell M, Babawale O, Liu H. Automated voxel classification used with atlas-guided diffuse optical tomography for assessment of functional brain networks in young and older adults. NEUROPHOTONICS 2016; 3:045002. [PMID: 27752518 PMCID: PMC5052324 DOI: 10.1117/1.nph.3.4.045002] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/16/2016] [Accepted: 09/12/2016] [Indexed: 05/12/2023]
Abstract
Atlas-guided diffuse optical tomography (atlas-DOT) is a computational means to image changes in cortical hemodynamic signals during human brain activities. Graph theory analysis (GTA) is a network analysis tool commonly used in functional neuroimaging to study brain networks. Atlas-DOT has not been analyzed with GTA to derive large-scale brain connectivity/networks based on near-infrared spectroscopy (NIRS) measurements. We introduced an automated voxel classification (AVC) method that facilitated the use of GTA with atlas-DOT images by grouping unequal-sized finite element voxels into anatomically meaningful regions of interest within the human brain. The overall approach included volume segmentation, AVC, and cross-correlation. To demonstrate the usefulness of AVC, we applied reproducibility analysis to resting-state functional connectivity measurements conducted from 15 young adults in a two-week period. We also quantified and compared changes in several brain network metrics between young and older adults, which were in agreement with those reported by a previous positron emission tomography study. Overall, this study demonstrated that AVC is a useful means for facilitating integration or combination of atlas-DOT with GTA and thus for quantifying NIRS-based, voxel-wise resting-state functional brain networks.
Collapse
Affiliation(s)
- Lin Li
- University of Texas at Arlington, Department of Bioengineering and Joint Graduate Program Between University of Texas at Arlington and University of Texas Southwestern Medical Center, Arlington, Texas 76019, United States
- University of California, Department of Neurology, David Geffen School of Medicine, Los Angeles, California 90095, United States
| | - Mary Cazzell
- Cook Children’s Medical Center, 801 Seventh Avenue, Fort Worth, Texas 76104, United States
| | - Olajide Babawale
- University of Texas at Arlington, Department of Bioengineering and Joint Graduate Program Between University of Texas at Arlington and University of Texas Southwestern Medical Center, Arlington, Texas 76019, United States
| | - Hanli Liu
- University of Texas at Arlington, Department of Bioengineering and Joint Graduate Program Between University of Texas at Arlington and University of Texas Southwestern Medical Center, Arlington, Texas 76019, United States
- Address all correspondence to: Hanli Liu, E-mail:
| |
Collapse
|
19
|
Yücel MA, Selb J, Aasted CM, Lin PY, Borsook D, Becerra L, Boas DA. Mayer waves reduce the accuracy of estimated hemodynamic response functions in functional near-infrared spectroscopy. BIOMEDICAL OPTICS EXPRESS 2016; 7:3078-88. [PMID: 27570699 PMCID: PMC4986815 DOI: 10.1364/boe.7.003078] [Citation(s) in RCA: 98] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/23/2016] [Revised: 06/30/2016] [Accepted: 07/05/2016] [Indexed: 05/07/2023]
Abstract
Analysis of cerebral hemodynamics reveals a wide spectrum of oscillations ranging from 0.0095 to 2 Hz. While most of these oscillations can be filtered out during analysis of functional near-infrared spectroscopy (fNIRS) signals when estimating stimulus evoked hemodynamic responses, oscillations around 0.1 Hz are an exception. This is due to the fact that they share a common spectral range with typical stimulus evoked hemodynamic responses from the brain. Here we investigate the effect of hemodynamic oscillations around 0.1 Hz on the estimation of hemodynamic response functions from fNIRS data. Our results show that for an expected response of ~1 µM in oxygenated hemoglobin concentration (HbO), Mayer wave oscillations with an amplitude > ~1 µM at 0.1 Hz reduce the accuracy of the estimated response as quantified by a 3 fold increase in the mean squared error and decrease in correlation (R(2) below 0.78) when compared to the true HRF. These results indicate that the amplitude of oscillations at 0.1 Hz can serve as an objective metric of the expected HRF estimation accuracy. In addition, we investigated the effect of short separation regression on the recovered HRF, and found that this improves the recovered HRF when large amplitude 0.1 Hz oscillations are present in fNIRS data. We suspect that the development of other filtering strategies may provide even further improvement.
Collapse
Affiliation(s)
- Meryem A. Yücel
- MGH/HST Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital (MGH), Harvard Medical School, Charlestown, MA, USA
| | - Juliette Selb
- MGH/HST Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital (MGH), Harvard Medical School, Charlestown, MA, USA
| | - Christopher M. Aasted
- Center for Pain and the Brain, Departments of Anaesthesia and Radiology, Boston Children’s Hospital, Boston, MA, USA
| | - Pei-Yi Lin
- MGH/HST Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital (MGH), Harvard Medical School, Charlestown, MA, USA
| | - David Borsook
- Center for Pain and the Brain, Departments of Anaesthesia and Radiology, Boston Children’s Hospital, Boston, MA, USA
| | - Lino Becerra
- Center for Pain and the Brain, Departments of Anaesthesia and Radiology, Boston Children’s Hospital, Boston, MA, USA
| | - David A. Boas
- MGH/HST Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital (MGH), Harvard Medical School, Charlestown, MA, USA
| |
Collapse
|
20
|
Chiarelli AM, Maclin EL, Low KA, Mathewson KE, Fabiani M, Gratton G. Combining energy and Laplacian regularization to accurately retrieve the depth of brain activity of diffuse optical tomographic data. JOURNAL OF BIOMEDICAL OPTICS 2016; 21:36008. [PMID: 26987429 PMCID: PMC4796096 DOI: 10.1117/1.jbo.21.3.036008] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/02/2015] [Accepted: 02/22/2016] [Indexed: 05/31/2023]
Abstract
Diffuse optical tomography (DOT) provides data about brain function using surface recordings. Despite recent advancements, an unbiased method for estimating the depth of absorption changes and for providing an accurate three-dimensional (3-D) reconstruction remains elusive. DOT involves solving an ill-posed inverse problem, requiring additional criteria for finding unique solutions. The most commonly used criterion is energy minimization (energy constraint). However, as measurements are taken from only one side of the medium (the scalp) and sensitivity is greater at shallow depths, the energy constraint leads to solutions that tend to be small and superficial. To correct for this bias, we combine the energy constraint with another criterion, minimization of spatial derivatives (Laplacian constraint, also used in low resolution electromagnetic tomography, LORETA). Used in isolation, the Laplacian constraint leads to solutions that tend to be large and deep. Using simulated, phantom, and actual brain activation data, we show that combining these two criteria results in accurate (error <2 mm) absorption depth estimates, while maintaining a two-point spatial resolution of <24 mm up to a depth of 30 mm. This indicates that accurate 3-D reconstruction of brain activity up to 30 mm from the scalp can be obtained with DOT.
Collapse
Affiliation(s)
- Antonio M. Chiarelli
- University of Illinois, Beckman Institute, 405 North Mathews Avenue, Urbana, Illinois 61801, United States
| | - Edward L. Maclin
- University of Illinois, Beckman Institute, 405 North Mathews Avenue, Urbana, Illinois 61801, United States
| | - Kathy A. Low
- University of Illinois, Beckman Institute, 405 North Mathews Avenue, Urbana, Illinois 61801, United States
| | - Kyle E. Mathewson
- University of Alberta, Department of Psychology, P217 Biological Sciences Building, Edmonton, Alberta T6G 2E9, Canada
| | - Monica Fabiani
- University of Illinois, Beckman Institute, 405 North Mathews Avenue, Urbana, Illinois 61801, United States
| | - Gabriele Gratton
- University of Illinois, Beckman Institute, 405 North Mathews Avenue, Urbana, Illinois 61801, United States
| |
Collapse
|
21
|
Hejazi SM, Sarkar S, Darezereshki Z. Fast multislice fluorescence molecular tomography using sparsity-inducing regularization. JOURNAL OF BIOMEDICAL OPTICS 2016; 21:26012. [PMID: 26927222 DOI: 10.1117/1.jbo.21.2.026012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/30/2015] [Accepted: 02/04/2016] [Indexed: 05/05/2023]
Abstract
Fluorescence molecular tomography (FMT) is a rapidly growing imaging method that facilitates the recovery of small fluorescent targets within biological tissue. The major challenge facing the FMT reconstruction method is the ill-posed nature of the inverse problem. In order to overcome this problem, the acquisition of large FMT datasets and the utilization of a fast FMT reconstruction algorithm with sparsity regularization have been suggested recently. Therefore, the use of a joint L1/total-variation (TV) regularization as a means of solving the ill-posed FMT inverse problem is proposed. A comparative quantified analysis of regularization methods based on L1-norm and TV are performed using simulated datasets, and the results show that the fast composite splitting algorithm regularization method can ensure the accuracy and robustness of the FMT reconstruction. The feasibility of the proposed method is evaluated in an in vivo scenario for the subcutaneous implantation of a fluorescent-dye-filled capillary tube in a mouse, and also using hybrid FMT and x-ray computed tomography data. The results show that the proposed regularization overcomes the difficulties created by the ill-posed inverse problem.
Collapse
Affiliation(s)
- Sedigheh Marjaneh Hejazi
- Tehran University of Medical Sciences, Medical Physics and Biomedical Engineering Department, School of Medicine, Tehran 1417613151, IranbTehran University of Medical Sciences, Research Center for Molecular and Cellular in Imaging, Bio-optical Imaging Gro
| | - Saeed Sarkar
- Tehran University of Medical Sciences, Medical Physics and Biomedical Engineering Department, School of Medicine, Tehran 1417613151, IrancTehran University of Medical Sciences, Research Center for Science and Technology in Medicine, Imam Khomeini Hospital
| | - Ziba Darezereshki
- Tehran University of Medical Sciences, Medical Physics and Biomedical Engineering Department, School of Medicine, Tehran 1417613151, Iran
| |
Collapse
|
22
|
Wu X, Eggebrecht AT, Ferradal SL, Culver JP, Dehghani H. Fast and efficient image reconstruction for high density diffuse optical imaging of the human brain. BIOMEDICAL OPTICS EXPRESS 2015; 6:4567-84. [PMID: 26601019 PMCID: PMC4646563 DOI: 10.1364/boe.6.004567] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/31/2015] [Revised: 10/14/2015] [Accepted: 10/16/2015] [Indexed: 05/18/2023]
Abstract
Real-time imaging of human brain has become an important technique within neuroimaging. In this study, a fast and efficient sensitivity map generation based on Finite Element Models (FEM) is developed which utilises a reduced sensitivitys matrix taking advantage of sparsity and parallelisation processes. Time and memory efficiency of these processes are evaluated and compared with conventional method showing that for a range of mesh densities from 50000 to 320000 nodes, the required memory is reduced over tenfold and computational time fourfold allowing for near real-time image recovery.
Collapse
Affiliation(s)
- Xue Wu
- School of Computer Science, University of Birmingham, Birmingham, B15 2TT, UK
| | - Adam T. Eggebrecht
- Department of Radiology, Washington University School of Medicine, 4525 Scott Avenue, St. Louis, MO, 63110, USA
| | - Silvina L. Ferradal
- Division of Newborn Medicine, Boston Children's Hospital, Harvard Medical School, Boston, MA, 02115, USA
| | - Joseph P. Culver
- Department of Radiology, Washington University School of Medicine, 4525 Scott Avenue, St. Louis, MO, 63110, USA
- Department of Biomedical Engineering, Washington University, One Brookings Drive, St. Louis, MO, 63130, USA
| | - Hamid Dehghani
- School of Computer Science, University of Birmingham, Birmingham, B15 2TT, UK
| |
Collapse
|
23
|
Mozumder M, Tarvainen T, Seppänen A, Nissilä I, Arridge SR, Kolehmainen V. Nonlinear approach to difference imaging in diffuse optical tomography. JOURNAL OF BIOMEDICAL OPTICS 2015; 20:105001. [PMID: 26440615 DOI: 10.1117/1.jbo.20.10.105001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/03/2015] [Accepted: 09/02/2015] [Indexed: 06/05/2023]
Abstract
Difference imaging aims at recovery of the change in the optical properties of a body based on measurements before and after the change. Conventionally, the image reconstruction is based on using difference of the measurements and a linear approximation of the observation model. One of the main benefits of the linearized difference reconstruction is that the approach has a good tolerance to modeling errors, which cancel out partially in the subtraction of the measurements. However, a drawback of the approach is that the difference images are usually only qualitative in nature and their spatial resolution can be weak because they rely on the global linearization of the nonlinear observation model. To overcome the limitations of the linear approach, we investigate a nonlinear approach for difference imaging where the images of the optical parameters before and after the change are reconstructed simultaneously based on the two datasets. We tested the feasibility of the method with simulations and experimental data from a phantom and studied how the approach tolerates modeling errors like domain truncation, optode coupling errors, and domain shape errors.
Collapse
Affiliation(s)
- Meghdoot Mozumder
- University of Eastern Finland, Department of Applied Physics, P.O. Box 1627, Kuopio 70211, Finland
| | - Tanja Tarvainen
- University of Eastern Finland, Department of Applied Physics, P.O. Box 1627, Kuopio 70211, FinlandbUniversity College London, Department of Computer Science, Gower Street, London WC1E 6BT, United Kingdom
| | - Aku Seppänen
- University of Eastern Finland, Department of Applied Physics, P.O. Box 1627, Kuopio 70211, Finland
| | - Ilkka Nissilä
- Aalto University School of Science, Department of Neuroscience and Biomedical Engineering, P.O. Box 12200, Aalto 00076, FinlanddHelsinki University Central Hospital, HUS Medical Imaging Center, BioMag Laboratory, P.O. Box 340, HUS 00029, Finland
| | - Simon R Arridge
- University College London, Department of Computer Science, Gower Street, London WC1E 6BT, United Kingdom
| | - Ville Kolehmainen
- University of Eastern Finland, Department of Applied Physics, P.O. Box 1627, Kuopio 70211, Finland
| |
Collapse
|
24
|
Kurihara K, Kawaguchi H, Obata T, Ito H, Okada E. Magnetic resonance imaging appropriate for construction of subject-specific head models for diffuse optical tomography. BIOMEDICAL OPTICS EXPRESS 2015; 6:3197-3209. [PMID: 26417492 PMCID: PMC4574648 DOI: 10.1364/boe.6.003197] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/17/2015] [Revised: 07/24/2015] [Accepted: 07/24/2015] [Indexed: 06/05/2023]
Abstract
Subject-specific head models of which their geometry is based on structural magnetic resonance images are essential to accurately estimate the spatial sensitivity profiles for image reconstruction in diffuse optical tomography. T1-weighted magnetic resonance images, which are commonly used for structural imaging, are not sufficient for the threshold-based segmentation of the superficial tissues. Two types of pulse sequences, which provide a high contrast among the superficial tissues, are introduced to complement the segmentation to construct the subject-specific head models. The magnetic resonance images acquired by the proposed pulse sequences are robust to the threshold level and adequate for the threshold-based segmentation of the superficial tissues compared to the T1- and T2-weighted images. The total scan time of the proposed pulse sequences is less than one-fourth of that for the T2-weighted pulse sequence.
Collapse
Affiliation(s)
- Kazuki Kurihara
- Department of Electronics and Electrical Engineering, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama, 223-8522, Japan
| | - Hiroshi Kawaguchi
- Human Informatics Research Institute, National Institute of Advanced Industrial Science and Technology, Central 6, 1-1-1 Higashi, Tsukuba, Ibaraki, 305-8566, Japan
- Molecular Imaging Center, National Institute of Radiological Sciences, 4-9-1 Anagawa, Inage-ku, Chiba, 263-8555, Japan
| | - Takayuki Obata
- Molecular Imaging Center, National Institute of Radiological Sciences, 4-9-1 Anagawa, Inage-ku, Chiba, 263-8555, Japan
- Research Center for Charged Particle Therapy, National Institute of Radiological Sciences, 4-9-1 Anagawa, Inage-ku, Chiba, 263-8555, Japan
| | - Hiroshi Ito
- Molecular Imaging Center, National Institute of Radiological Sciences, 4-9-1 Anagawa, Inage-ku, Chiba, 263-8555, Japan
- Advanced Clinical Research Center, Fukushima Medical University, 1 Hikarigaoka, Fukushima, 960-1295, Japan
| | - Eiji Okada
- Department of Electronics and Electrical Engineering, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama, 223-8522, Japan
| |
Collapse
|
25
|
Wu X, Eggebrecht AT, Ferradal SL, Culver JP, Dehghani H. Evaluation of rigid registration methods for whole head imaging in diffuse optical tomography. NEUROPHOTONICS 2015; 2:035002. [PMID: 26217675 PMCID: PMC4509792 DOI: 10.1117/1.nph.2.3.035002] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/20/2015] [Accepted: 06/18/2015] [Indexed: 05/22/2023]
Abstract
Functional brain imaging has become an important neuroimaging technique for the study of brain organization and development. Compared to other imaging techniques, diffuse optical tomography (DOT) is a portable and low-cost technique that can be applied to infants and hospitalized patients using an atlas-based light model. For DOT imaging, the accuracy of the forward model has a direct effect on the resulting recovered brain function within a field of view and so the accuracy of the spatially normalized atlas-based forward models must be evaluated. Herein, the accuracy of atlas-based DOT is evaluated on models that are spatially normalized via a number of different rigid registration methods on 24 subjects. A multileveled approach is developed to evaluate the correlation of the geometrical and sensitivity accuracies across the full field of view as well as within specific functional subregions. Results demonstrate that different registration methods are optimal for recovery of different sets of functional brain regions. However, the "nearest point to point" registration method, based on the EEG 19 landmark system, is shown to be the most appropriate registration method for image quality throughout the field of view of the high-density cap that covers the whole of the optically accessible cortex.
Collapse
Affiliation(s)
- Xue Wu
- University of Birmingham, School of Computer Science, Edgbaston, Birmingham B15 2TT, United Kingdom
| | - Adam T. Eggebrecht
- Washington University School of Medicine, Department of Radiology, 4525 Scott Avenue, St. Louis, Missouri 63110, United States
| | - Silvina L. Ferradal
- Fetal-Neonatal Neuroimaging and Developmental Science Center, Boston Children’s Hospital, 300 Longwood Avenue, Boston, Massachusetts 02115, United States
| | - Joseph P. Culver
- Washington University School of Medicine, Department of Radiology, 4525 Scott Avenue, St. Louis, Missouri 63110, United States
- Washington University, Department of Biomedical Engineering, One Brookings Drive, St. Louis, Missouri 63130, United States
| | - Hamid Dehghani
- University of Birmingham, School of Computer Science, Edgbaston, Birmingham B15 2TT, United Kingdom
- Address all correspondence to: Hamid Dehghani, E-mail:
| |
Collapse
|
26
|
Davies DJ, Su Z, Clancy MT, Lucas SJE, Dehghani H, Logan A, Belli A. Near-Infrared Spectroscopy in the Monitoring of Adult Traumatic Brain Injury: A Review. J Neurotrauma 2015; 32:933-41. [PMID: 25603012 DOI: 10.1089/neu.2014.3748] [Citation(s) in RCA: 94] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Cerebral near-infrared spectroscopy (NIRS) has long represented an exciting prospect for the noninvasive monitoring of cerebral tissue oxygenation and perfusion in the context of traumatic brain injury (TBI), although uncertainty still exists regarding the reliability of this technology specifically within this field. We have undertaken a review of the existing literature relating to the application of NIRS within TBI. We discuss current "state-of-the-art" NIRS monitoring, provide a brief background of the technology, and discuss the evidence regarding the ability of NIRS to substitute for established invasive monitoring in TBI.
Collapse
Affiliation(s)
- David J Davies
- 1 Department of Neurosurgery Clinical Research, Queen Elizabeth Hospital , Edgbaston, Birmingham, United Kingdom
| | - Zhangjie Su
- 1 Department of Neurosurgery Clinical Research, Queen Elizabeth Hospital , Edgbaston, Birmingham, United Kingdom
| | - Michael T Clancy
- 2 School of Computational Science Medical Imaging Group, University of Birmingham , Edgbaston, Birmingham, United Kingdom
| | - Samuel J E Lucas
- 3 Department of Exercise Physiology, School of Sport, Exercise, and Rehabilitation Sciences, University of Birmingham , Edgbaston, Birmingham, United Kingdom
| | - Hamid Dehghani
- 4 Department of Medical Imaging, School of Computer Science, University of Birmingham , Edgbaston, Birmingham, United Kingdom
| | - Ann Logan
- 5 Department of Molecular Neuroscience, School of Clinical and Experimental Medicine, University of Birmingham , Edgbaston, Birmingham, United Kingdom
| | - Antonio Belli
- 6 Department of Surgical Neurology, National Institute for Health Research, Queen Elizabeth Hospital , Edgbaston, Birmingham, United Kingdom
| |
Collapse
|
27
|
Aasted CM, Yücel MA, Cooper RJ, Dubb J, Tsuzuki D, Becerra L, Petkov MP, Borsook D, Dan I, Boas DA. Anatomical guidance for functional near-infrared spectroscopy: AtlasViewer tutorial. NEUROPHOTONICS 2015; 2:020801. [PMID: 26157991 PMCID: PMC4478785 DOI: 10.1117/1.nph.2.2.020801] [Citation(s) in RCA: 235] [Impact Index Per Article: 26.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/20/2014] [Accepted: 04/02/2015] [Indexed: 05/18/2023]
Abstract
Functional near-infrared spectroscopy (fNIRS) is an optical imaging method that is used to noninvasively measure cerebral hemoglobin concentration changes induced by brain activation. Using structural guidance in fNIRS research enhances interpretation of results and facilitates making comparisons between studies. AtlasViewer is an open-source software package we have developed that incorporates multiple spatial registration tools to enable structural guidance in the interpretation of fNIRS studies. We introduce the reader to the layout of the AtlasViewer graphical user interface, the folder structure, and user files required in the creation of fNIRS probes containing sources and detectors registered to desired locations on the head, evaluating probe fabrication error and intersubject probe placement variability, and different procedures for estimating measurement sensitivity to different brain regions as well as image reconstruction performance. Further, we detail how AtlasViewer provides a generic head atlas for guiding interpretation of fNIRS results, but also permits users to provide subject-specific head anatomies to interpret their results. We anticipate that AtlasViewer will be a valuable tool in improving the anatomical interpretation of fNIRS studies.
Collapse
Affiliation(s)
- Christopher M. Aasted
- Center for Pain and the Brain, Harvard Medical School, 1 Autumn Street, Boston, Massachusetts 02215, United States
- Department of Anesthesiology, Perioperative, and Pain Medicine, Boston Children’s Hospital, 300 Longwood Avenue, Boston, Massachusetts 02115, United States
| | - Meryem A. Yücel
- Department of Radiology, Athinoula Martinos Center for Biomedical Imaging, 149 13th Street, Charlestown, Massachusetts 02129, United States
| | - Robert J. Cooper
- University College London, Department of Medical Physics and Biomedical Engineering, Gower Street, London WC1E 6BT, United Kingdom
| | - Jay Dubb
- Department of Radiology, Athinoula Martinos Center for Biomedical Imaging, 149 13th Street, Charlestown, Massachusetts 02129, United States
| | - Daisuke Tsuzuki
- Chuo University, Faculty of Science and Engineering, 1-13-27 Kasuga, Bunkyo, Tokyo 112-8551, Japan
| | - Lino Becerra
- Center for Pain and the Brain, Harvard Medical School, 1 Autumn Street, Boston, Massachusetts 02215, United States
- Department of Anesthesiology, Perioperative, and Pain Medicine, Boston Children’s Hospital, 300 Longwood Avenue, Boston, Massachusetts 02115, United States
- Department of Radiology, Boston Children’s Hospital, 300 Longwood Avenue, Boston, Massachusetts 02115, United States
- Department of Psychiatry, McLean Hospital, 115 Mill Street, Belmont, Massachusetts 02478, United States
| | - Mike P. Petkov
- Department of Anesthesiology, Perioperative, and Pain Medicine, Boston Children’s Hospital, 300 Longwood Avenue, Boston, Massachusetts 02115, United States
| | - David Borsook
- Center for Pain and the Brain, Harvard Medical School, 1 Autumn Street, Boston, Massachusetts 02215, United States
- Department of Anesthesiology, Perioperative, and Pain Medicine, Boston Children’s Hospital, 300 Longwood Avenue, Boston, Massachusetts 02115, United States
- Department of Radiology, Boston Children’s Hospital, 300 Longwood Avenue, Boston, Massachusetts 02115, United States
- Department of Psychiatry, McLean Hospital, 115 Mill Street, Belmont, Massachusetts 02478, United States
| | - Ippeita Dan
- Chuo University, Faculty of Science and Engineering, 1-13-27 Kasuga, Bunkyo, Tokyo 112-8551, Japan
| | - David A. Boas
- Department of Radiology, Athinoula Martinos Center for Biomedical Imaging, 149 13th Street, Charlestown, Massachusetts 02129, United States
| |
Collapse
|
28
|
Wu X, Eggebrecht AT, Ferradal SL, Culver JP, Dehghani H. Quantitative evaluation of atlas-based high-density diffuse optical tomography for imaging of the human visual cortex. BIOMEDICAL OPTICS EXPRESS 2014; 5:3882-900. [PMID: 25426318 PMCID: PMC4242025 DOI: 10.1364/boe.5.003882] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/13/2014] [Revised: 08/21/2014] [Accepted: 09/25/2014] [Indexed: 05/24/2023]
Abstract
Image recovery in diffuse optical tomography (DOT) of the human brain often relies on accurate models of light propagation within the head. In the absence of subject specific models for image reconstruction, the use of atlas based models are showing strong promise. Although there exists some understanding in the use of some limited rigid model registrations in DOT, there has been a lack of a detailed analysis between errors in geometrical accuracy, light propagation in tissue and subsequent errors in dynamic imaging of recovered focal activations in the brain. In this work 11 different rigid registration algorithms, across 24 simulated subjects, are evaluated for DOT studies in the visual cortex. Although there exists a strong correlation (R(2) = 0.97) between geometrical surface error and internal light propagation errors, the overall variation is minimal when analysing recovered focal activations in the visual cortex. While a subject specific mesh gives the best results with a 1.2 mm average location error, no single algorithm provides errors greater than 4.5 mm. This work demonstrates that the use of rigid algorithms for atlas based imaging is a promising route when subject specific models are not available.
Collapse
Affiliation(s)
- Xue Wu
- School of Computer Science, University of Birmingham, Birmingham, B15 2TT,
UK
| | - Adam T. Eggebrecht
- Department of Radiology, Washington University School of Medicine, 4525 Scott Avenue, St Louis, MO, 63110,
USA
| | - Silvina L Ferradal
- Department of Radiology, Washington University School of Medicine, 4525 Scott Avenue, St Louis, MO, 63110,
USA
- Department of Biomedical Engineering, Washington University, One Brookings Drive, St. Louis, MO, 63130,
USA
| | - Joseph P. Culver
- Department of Radiology, Washington University School of Medicine, 4525 Scott Avenue, St Louis, MO, 63110,
USA
- Department of Biomedical Engineering, Washington University, One Brookings Drive, St. Louis, MO, 63130,
USA
| | - Hamid Dehghani
- School of Computer Science, University of Birmingham, Birmingham, B15 2TT,
UK
| |
Collapse
|
29
|
Giacometti P, Diamond SG. Correspondence of electroencephalography and near-infrared spectroscopy sensitivities to the cerebral cortex using a high-density layout. NEUROPHOTONICS 2014; 1:025001. [PMID: 25558462 PMCID: PMC4280681 DOI: 10.1117/1.nph.1.2.025001] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
This study investigates the correspondence of the cortical sensitivity of electroencephalography (EEG) and near-infrared spectroscopy (NIRS). EEG forward model sensitivity to the cerebral cortex was calculated for 329 EEG electrodes following the 10-5 EEG positioning system using a segmented structural magnetic resonance imaging scan of a human subject. NIRS forward model sensitivity was calculated for the same subject using 156 NIRS source-detector pairs selected from 32 source and 32 detector optodes positioned on the scalp using a subset of the 10-5 EEG positioning system. Sensitivity correlations between colocalized NIRS source-detector pair groups and EEG channels yielded R = 0.46 ± 0.08. Groups of NIRS source-detector pairs with maximum correlations to EEG electrode sensitivities are tabulated. The mean correlation between the point spread functions for EEG and NIRS regions of interest (ROI) was R = 0.43 ± 0.07. Spherical ROIs with radii of 26 mm yielded the maximum correlation between EEG and NIRS averaged across all cortical mesh nodes. These sensitivity correlations between EEG and NIRS should be taken into account when designing multimodal studies of neurovascular coupling and when using NIRS as a statistical prior for EEG source localization.
Collapse
Affiliation(s)
- Paolo Giacometti
- Thayer School of Engineering at Dartmouth, 14 Engineering Drive, Hanover, New Hampshire 03755, United States
- Address all correspondence to: Paolo Giacometti, E-mail:
| | - Solomon G. Diamond
- Thayer School of Engineering at Dartmouth, 14 Engineering Drive, Hanover, New Hampshire 03755, United States
| |
Collapse
|
30
|
Varma HM, Valdes CP, Kristoffersen AK, Culver JP, Durduran T. Speckle contrast optical tomography: A new method for deep tissue three-dimensional tomography of blood flow. BIOMEDICAL OPTICS EXPRESS 2014; 5:1275-89. [PMID: 24761306 PMCID: PMC3986001 DOI: 10.1364/boe.5.001275] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/31/2014] [Revised: 03/10/2014] [Accepted: 03/11/2014] [Indexed: 05/18/2023]
Abstract
A novel tomographic method based on the laser speckle contrast, speckle contrast optical tomography (SCOT) is introduced that allows us to reconstruct three dimensional distribution of blood flow in deep tissues. This method is analogous to the diffuse optical tomography (DOT) but for deep tissue blood flow. We develop a reconstruction algorithm based on first Born approximation to generate three dimensional distribution of flow using the experimental data obtained from tissue simulating phantoms.
Collapse
Affiliation(s)
- Hari M. Varma
- ICFO-Institut de Ciències Fotòniques, Av. Carl Friedrich Gauss, 3, 08860, Castelldefels (Barcelona),
Spain
- These authors contributed equally to this work
| | - Claudia P. Valdes
- ICFO-Institut de Ciències Fotòniques, Av. Carl Friedrich Gauss, 3, 08860, Castelldefels (Barcelona),
Spain
- These authors contributed equally to this work
| | - Anna K. Kristoffersen
- ICFO-Institut de Ciències Fotòniques, Av. Carl Friedrich Gauss, 3, 08860, Castelldefels (Barcelona),
Spain
| | - Joseph P. Culver
- Department of Radiology, Washington University School of Medicine, St. Louis, MO 63110,
USA
- Department of Physics, Washington University, St. Louis, MO 63130,
USA
| | - Turgut Durduran
- ICFO-Institut de Ciències Fotòniques, Av. Carl Friedrich Gauss, 3, 08860, Castelldefels (Barcelona),
Spain
| |
Collapse
|
31
|
Lin ZJ, Li L, Cazzell M, Liu H. Atlas-guided volumetric diffuse optical tomography enhanced by generalized linear model analysis to image risk decision-making responses in young adults. Hum Brain Mapp 2014; 35:4249-66. [PMID: 24619964 PMCID: PMC4282392 DOI: 10.1002/hbm.22459] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2013] [Revised: 11/25/2013] [Accepted: 12/18/2013] [Indexed: 11/06/2022] Open
Abstract
Diffuse optical tomography (DOT) is a variant of functional near infrared spectroscopy and has the capability of mapping or reconstructing three dimensional (3D) hemodynamic changes due to brain activity. Common methods used in DOT image analysis to define brain activation have limitations because the selection of activation period is relatively subjective. General linear model (GLM)-based analysis can overcome this limitation. In this study, we combine the atlas-guided 3D DOT image reconstruction with GLM-based analysis (i.e., voxel-wise GLM analysis) to investigate the brain activity that is associated with risk decision-making processes. Risk decision-making is an important cognitive process and thus is an essential topic in the field of neuroscience. The Balloon Analog Risk Task (BART) is a valid experimental model and has been commonly used to assess human risk-taking actions and tendencies while facing risks. We have used the BART paradigm with a blocked design to investigate brain activations in the prefrontal and frontal cortical areas during decision-making from 37 human participants (22 males and 15 females). Voxel-wise GLM analysis was performed after a human brain atlas template and a depth compensation algorithm were combined to form atlas-guided DOT images. In this work, we wish to demonstrate the excellence of using voxel-wise GLM analysis with DOT to image and study cognitive functions in response to risk decision-making. Results have shown significant hemodynamic changes in the dorsal lateral prefrontal cortex (DLPFC) during the active-choice mode and a different activation pattern between genders; these findings correlate well with published literature in functional magnetic resonance imaging (fMRI) and fNIRS studies.
Collapse
Affiliation(s)
- Zi-Jing Lin
- Department of Bioengineering, Joint Program of Biomedical Engineering between University of Texas at Arlington and University of Texas Southwestern Medical Center at Dallas, University of Texas at Arlington, Arlington, Texas; National Synchrotron Radiation Research Center, Hsinchu, Taiwan
| | | | | | | |
Collapse
|
32
|
Lin ZJ, Ren M, Li L, Liu Y, Su J, Yang SH, Liu H. Interleaved imaging of cerebral hemodynamics and blood flow index to monitor ischemic stroke and treatment in rat by volumetric diffuse optical tomography. Neuroimage 2013; 85 Pt 1:566-82. [PMID: 23872158 DOI: 10.1016/j.neuroimage.2013.07.020] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2013] [Revised: 07/05/2013] [Accepted: 07/06/2013] [Indexed: 12/29/2022] Open
Abstract
Diffuse optical tomography (DOT) has been used by several groups to assess cerebral hemodynamics of cerebral ischemia in humans and animals. In this study, we combined DOT with an indocyanine green (ICG)-tracking method to achieve interleaved images of cerebral hemodynamics and blood flow index (BFI) using two middle cerebral artery occlusion (MCAO) rat models. To achieve volumetric images with high-spatial resolution, we first integrated a depth compensation algorithm (DCA) with a volumetric mesh-based rat head model to generate three-dimensional (3D) DOT on a rat brain atlas. Then, the experimental DOT data from two rat models were collected using interleaved strategy for cerebral hemodynamics and BFI during and after ischemic stroke, with and without a thrombolytic therapy for the embolic MCAO model. The acquired animal data were further analyzed using the integrated rat-atlas-guided DOT method to form time-evolving 3D images of both cerebral hemodynamics and BFI. In particular, we were able to show and identify therapeutic outcomes of a thrombolytic treatment applied to the embolism-induced ischemic model. This paper demonstrates that volumetric DOT is capable of providing high-quality, interleaved images of cerebral hemodynamics and blood perfusion in small animals during and after ischemic stroke, with excellent 3D visualization and quantifications.
Collapse
Affiliation(s)
- Zi-Jing Lin
- Department of Bioengineering, Joint Graduate Program between University of Texas at Arlington and University of Texas Southwestern Medical Center, University of Texas at Arlington, TX 76019, USA
| | | | | | | | | | | | | |
Collapse
|
33
|
Tian F, Liu H. Depth-compensated diffuse optical tomography enhanced by general linear model analysis and an anatomical atlas of human head. Neuroimage 2013; 85 Pt 1:166-80. [PMID: 23859922 DOI: 10.1016/j.neuroimage.2013.07.016] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2013] [Revised: 05/24/2013] [Accepted: 07/04/2013] [Indexed: 11/17/2022] Open
Abstract
One of the main challenges in functional diffuse optical tomography (DOT) is to accurately recover the depth of brain activation, which is even more essential when differentiating true brain signals from task-evoked artifacts in the scalp. Recently, we developed a depth-compensated algorithm (DCA) to minimize the depth localization error in DOT. However, the semi-infinite model that was used in DCA deviated significantly from the realistic human head anatomy. In the present work, we incorporated depth-compensated DOT (DC-DOT) with a standard anatomical atlas of human head. Computer simulations and human measurements of sensorimotor activation were conducted to examine and prove the depth specificity and quantification accuracy of brain atlas-based DC-DOT. In addition, node-wise statistical analysis based on the general linear model (GLM) was also implemented and performed in this study, showing the robustness of DC-DOT that can accurately identify brain activation at the correct depth for functional brain imaging, even when co-existing with superficial artifacts.
Collapse
Affiliation(s)
- Fenghua Tian
- Department of Bioengineering, Joint Program in Biomedical Engineering between UT Arlington and UT Southwestern Medical Center at Dallas, University of Texas at Arlington, Arlington, TX, USA
| | | |
Collapse
|
34
|
Ferradal SL, Eggebrecht AT, Hassanpour M, Snyder AZ, Culver JP. Atlas-based head modeling and spatial normalization for high-density diffuse optical tomography: in vivo validation against fMRI. Neuroimage 2013; 85 Pt 1:117-26. [PMID: 23578579 DOI: 10.1016/j.neuroimage.2013.03.069] [Citation(s) in RCA: 65] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2013] [Revised: 03/27/2013] [Accepted: 03/28/2013] [Indexed: 10/27/2022] Open
Abstract
Diffuse optical imaging (DOI) is increasingly becoming a valuable neuroimaging tool when fMRI is precluded. Recent developments in high-density diffuse optical tomography (HD-DOT) overcome previous limitations of sparse DOI systems, providing improved image quality and brain specificity. These improvements in instrumentation prompt the need for advancements in both i) realistic forward light modeling for accurate HD-DOT image reconstruction, and ii) spatial normalization for voxel-wise comparisons across subjects. Individualized forward light models derived from subject-specific anatomical images provide the optimal inverse solutions, but such modeling may not be feasible in all situations. In the absence of subject-specific anatomical images, atlas-based head models registered to the subject's head using cranial fiducials provide an alternative solution. In addition, a standard atlas is attractive because it defines a common coordinate space in which to compare results across subjects. The question therefore arises as to whether atlas-based forward light modeling ensures adequate HD-DOT image quality at the individual and group level. Herein, we demonstrate the feasibility of using atlas-based forward light modeling and spatial normalization methods. Both techniques are validated using subject-matched HD-DOT and fMRI data sets for visual evoked responses measured in five healthy adult subjects. HD-DOT reconstructions obtained with the registered atlas anatomy (i.e. atlas DOT) had an average localization error of 2.7mm relative to reconstructions obtained with the subject-specific anatomical images (i.e. subject-MRI DOT), and 6.6mm relative to fMRI data. At the group level, the localization error of atlas DOT reconstruction was 4.2mm relative to subject-MRI DOT reconstruction, and 6.1mm relative to fMRI. These results show that atlas-based image reconstruction provides a viable approach to individual head modeling for HD-DOT when anatomical imaging is not available.
Collapse
Affiliation(s)
- Silvina L Ferradal
- Department of Biomedical Engineering, Washington University, Whitaker Hall, One Brookings Dr., St. Louis, MO, 63130, USA; Department of Radiology, Washington University School of Medicine, East Bldg., 4525 Scott Ave, St. Louis, MO, 63110, USA
| | | | | | | | | |
Collapse
|
35
|
Jha AK, Clarkson E, Kupinski MA. An ideal-observer framework to investigate signal detectability in diffuse optical imaging. BIOMEDICAL OPTICS EXPRESS 2013; 4:2107-23. [PMID: 24156068 PMCID: PMC3799670 DOI: 10.1364/boe.4.002107] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/12/2013] [Revised: 08/21/2013] [Accepted: 08/26/2013] [Indexed: 05/12/2023]
Abstract
With the emergence of diffuse optical tomography (DOT) as a non-invasive imaging modality, there is a requirement to evaluate the performance of the developed DOT systems on clinically relevant tasks. One such important task is the detection of high-absorption signals in the tissue. To investigate signal detectability in DOT systems for system optimization, an appropriate approach is to use the Bayesian ideal observer, but this observer is computationally very intensive. It has been shown that the Fisher information can be used as a surrogate figure of merit (SFoM) that approximates the ideal observer performance. In this paper, we present a theoretical framework to use the Fisher information for investigating signal detectability in DOT systems. The usage of Fisher information requires evaluating the gradient of the photon distribution function with respect to the absorption coefficients. We derive the expressions to compute the gradient of the photon distribution function with respect to the scattering and absorption coefficients. We find that computing these gradients simply requires executing the radiative transport equation with a different source term. We then demonstrate the application of the SFoM to investigate signal detectability in DOT by performing various simulation studies, which help to validate the proposed framework and also present some insights on signal detectability in DOT.
Collapse
Affiliation(s)
- Abhinav K. Jha
- College of Optical Sciences, University of Arizona, Tucson, AZ,
USA
| | - Eric Clarkson
- College of Optical Sciences, University of Arizona, Tucson, AZ,
USA
- Department of Medical Imaging, University of Arizona, Tucson, AZ,
USA
| | - Matthew A. Kupinski
- College of Optical Sciences, University of Arizona, Tucson, AZ,
USA
- Department of Medical Imaging, University of Arizona, Tucson, AZ,
USA
| |
Collapse
|
36
|
Zhan Y, Eggebrecht AT, Culver JP, Dehghani H. Singular value decomposition based regularization prior to spectral mixing improves crosstalk in dynamic imaging using spectral diffuse optical tomography. BIOMEDICAL OPTICS EXPRESS 2012; 3:2036-49. [PMID: 23024899 PMCID: PMC3447547 DOI: 10.1364/boe.3.002036] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/18/2012] [Revised: 08/02/2012] [Accepted: 08/04/2012] [Indexed: 05/24/2023]
Abstract
The spectrally constrained diffuse optical tomography (DOT) method relies on incorporating spectral prior information directly into the image reconstruction algorithm, thereby correlating the underlying optical properties across multiple wavelengths. Although this method has been shown to provide a solution that is stable, the use of conventional Tikhonov-type regularization techniques can lead to additional crosstalk between parameters, particularly in linear, single-step dynamic imaging applications. This is due mainly to the suboptimal regularization of the spectral Jacobian matrix, which smoothes not only the image-data space, but also the spectral mapping space. In this work a novel regularization technique based on the singular value decomposition (SVD) is presented that preserves the spectral prior information while regularizing the Jacobian matrix, leading to dramatically reduced crosstalk between the recovered parameters. Using simulated data, images of changes in oxygenated and deoxygenated hemoglobin concentrations are reconstructed via the SVD-based approach and compared with images reconstructed by using non-spectral and conventional spectral methods. In a 2D, two wavelength example, it is shown that the proposed approach provides a 98% reduction in crosstalk between recovered parameters as compared with conventional spectral reconstruction algorithms, and 60% as compared with non-spectrally constrained algorithms. Using a subject specific multilayered model of the human head, a noiseless dynamic simulation of cortical activation is performed to further demonstrate such improvement in crosstalk. However, with the addition of realistic noise in the data, both non-spectral and proposed algorithms perform similarly, indicating that the use of spectrally constrained reconstruction algorithms in dynamic DOT may be limited by the contrast of the signal as well as the noise characteristics of the system.
Collapse
Affiliation(s)
- Yuxuan Zhan
- School of Computer Science, University of Birmingham, Birmingham, B15 2TT, UK
| | - Adam T. Eggebrecht
- Department of Radiology, Washington University School of Medicine, 4525 Scott Avenue, St Louis, MO, 63110, USA
| | - Joseph P. Culver
- Department of Radiology, Washington University School of Medicine, 4525 Scott Avenue, St Louis, MO, 63110, USA
| | - Hamid Dehghani
- School of Computer Science, University of Birmingham, Birmingham, B15 2TT, UK
| |
Collapse
|