Validating atlas-guided DOT: a comparison of diffuse optical tomography informed by atlas and subject-specific anatomies

Evaluation of rigid registration methods for whole head imaging in diffuse optical tomography

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.

Anatomical guidance for functional near-infrared spectroscopy: AtlasViewer tutorial

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.

Dynamic causal modelling for functional near-infrared spectroscopy

Functional near-infrared spectroscopy (fNIRS) is an emerging technique for measuring changes in cerebral hemoglobin concentration via optical absorption changes. Although there is great interest in using fNIRS to study brain connectivity, current methods are unable to infer the directionality of neuronal connections. In this paper, we apply Dynamic Causal Modelling (DCM) to fNIRS data. Specifically, we present a generative model of how observed fNIRS data are caused by interactions among hidden neuronal states. Inversion of this generative model, using an established Bayesian framework (variational Laplace), then enables inference about changes in directed connectivity at the neuronal level. Using experimental data acquired during motor imagery and motor execution tasks, we show that directed (i.e., effective) connectivity from the supplementary motor area to the primary motor cortex is negatively modulated by motor imagery, and this suppressive influence causes reduced activity in the primary motor cortex during motor imagery. These results are consistent with findings of previous functional magnetic resonance imaging (fMRI) studies, suggesting that the proposed method enables one to infer directed interactions in the brain mediated by neuronal dynamics from measurements of optical density changes.

Data-driven approach to optimum wavelength selection for diffuse optical imaging

The production of accurate and independent images of the changes in concentration of oxyhemoglobin and deoxyhemoglobin by diffuse optical imaging is heavily dependent on which wavelengths of near-infrared light are chosen to interrogate the target tissue. Although wavelengths can be selected by theoretical methods, in practice the accuracy of reconstructed images will be affected by wavelength-specific and system-specific factors such as laser source power and detector sensitivity. We describe the application of a data-driven approach to optimum wavelength selection for the second generation of University College London's multichannel, time-domain optical tomography system (MONSTIR II). By performing a functional activation experiment using 12 different wavelengths between 690 and 870 nm, we were able to identify the combinations of 2, 3, and 4 wavelengths which most accurately reproduced the results obtained using all 12 wavelengths via an imaging approach. Our results show that the set of 2, 3, and 4 wavelengths which produce the most accurate images of functional activation are [770, 810], [770, 790, 850], and [730, 770, 810, 850] respectively, but also that the system is relatively robust to wavelength selection within certain limits. Although these results are specific to MONSTIR II, the approach we developed can be applied to other multispectral near-infrared spectroscopy and optical imaging systems.

Quantitative evaluation of atlas-based high-density diffuse optical tomography for imaging of the human visual cortex

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.

Sensitivity of near-infrared spectroscopy and diffuse correlation spectroscopy to brain hemodynamics: simulations and experimental findings during hypercapnia

Near-infrared spectroscopy (NIRS) and diffuse correlation spectroscopy (DCS) are two diffuse optical technologies for brain imaging that are sensitive to changes in hemoglobin concentrations and blood flow, respectively. Measurements for both modalities are acquired on the scalp, and therefore hemodynamic processes in the extracerebral vasculature confound the interpretation of cortical hemodynamic signals. The sensitivity of NIRS to the brain versus the extracerebral tissue and the contrast-to-noise ratio (CNR) of NIRS to cerebral hemodynamic responses have been well characterized, but the same has not been evaluated for DCS. This is important to assess in order to understand their relative capabilities in measuring cerebral physiological changes. We present Monte Carlo simulations on a head model that demonstrate that the relative brain-to-scalp sensitivity is about three times higher for DCS (0.3 at 3 cm) than for NIRS (0.1 at 3 cm). However, because DCS has higher levels of noise due to photon-counting detection, the CNR is similar for both modalities in response to a physiologically realistic simulation of brain activation. Even so, we also observed higher CNR of the hemodynamic response during graded hypercapnia in adult subjects with DCS than with NIRS.

A 4D neonatal head model for diffuse optical imaging of pre-term to term infants

Diffuse optical tomography is most accurate when an individual's MRI data can be used as a spatial prior for image reconstruction and for visualization of the resulting images of changes in oxy- and deoxy-hemoglobin concentration. As this necessitates an MRI scan to be performed for each study, which undermines many of the advantages of diffuse optical methods, the use of registered atlases to model the individual's anatomy is becoming commonplace. Infant studies require carefully age-matched atlases because of the rapid growth and maturation of the infant brain. In this paper, we present a 4D neonatal head model which, for each week from 29 to 44 weeks post-menstrual age, includes: 1) a multi-layered tissue mask which identifies extra-cerebral layers, cerebrospinal fluid, gray matter, white matter, cerebellum and brainstem, 2) a high-density tetrahedral head mesh, 3) surface meshes for the scalp, gray-matter and white matter layers and 4) cranial landmarks and 10-5 locations on the scalp surface. This package, freely available online at www.ucl.ac.uk/medphys/research/4dneonatalmodel can be applied by users of near-infrared spectroscopy and diffuse optical tomography to optimize probe locations, optimize image reconstruction, register data to cortical locations and ultimately improve the accuracy and interpretation of diffuse optical techniques in newborn populations.

Atlas-guided volumetric diffuse optical tomography enhanced by generalized linear model analysis to image risk decision-making responses in young adults

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.

T1 magnetic resonance imaging head segmentation for diffuse optical tomography and electroencephalography

Accurate segmentation of structural magnetic resonance images is critical for creating subject-specific forward models for functional neuroimaging source localization. In this work, we present an innovative segmentation algorithm that generates accurate head tissue layer thicknesses that are needed for diffuse optical tomography (DOT) data analysis. The presented algorithm is compared against other publicly available head segmentation methods. The proposed algorithm has a root mean square scalp thickness error of 1.60 mm, skull thickness error of 1.96 mm, and summed scalp and skull error of 1.49 mm. We also introduce a segmentation evaluation metric that evaluates the accuracy of tissue layer thicknesses in regions of the head where optodes are typically placed. The presented segmentation algorithm and evaluation metric are tools for improving the localization accuracy of neuroimaging with DOT, and also multimodal neuroimaging such as combined electroencephalography and DOT.

Probing the early development of visual working memory capacity with functional near-infrared spectroscopy

Visual working memory (VWM) is a core cognitive system with a highly limited capacity. The present study is the first to examine VWM capacity limits in early development using functional neuroimaging. We recorded optical neuroimaging data while 3- and 4-year-olds completed a change detection task where they detected changes in the shapes of objects after a brief delay. Near-infrared sources and detectors were placed over the following 10-20 positions: F3 and F5 in left frontal cortex, F4 and F6 in right frontal cortex, P3 and P5 in left parietal cortex, and P4 and P6 in right parietal cortex. The first question was whether we would see robust task-specific activation of the frontal-parietal network identified in the adult fMRI literature. This was indeed the case: three left frontal channels and 11 of 12 parietal channels showed a statistically robust difference between the concentration of oxygenated and deoxygenated hemoglobin following the presentation of the sample array. Moreover, four channels in the left hemisphere near P3, P5, and F5 showed a robust increase as the working memory load increased from 1 to 3 items. Notably, the hemodynamic response did not asymptote at 1-2 items as expected from previous fMRI studies with adults. Finally, 4-year-olds showed a more robust parietal response relative to 3-year-olds, and an increasing sensitivity to the memory load manipulation. These results demonstrate that fNIRS is an effective tool to study the neural processes that underlie the early development of VWM capacity.

Comparison of a layered slab and an atlas head model for Monte Carlo fitting of time-domain near-infrared spectroscopy data of the adult head

Near-infrared spectroscopy (NIRS) estimations of the adult brain baseline optical properties based on a homogeneous model of the head are known to introduce significant contamination from extracerebral layers. More complex models have been proposed and occasionally applied to in vivo data, but their performances have never been characterized on realistic head structures. Here we implement a flexible fitting routine of time-domain NIRS data using graphics processing unit based Monte Carlo simulations. We compare the results for two different geometries: a two-layer slab with variable thickness of the first layer and a template atlas head registered to the subject's head surface. We characterize the performance of the Monte Carlo approaches for fitting the optical properties from simulated time-resolved data of the adult head. We show that both geometries provide better results than the commonly used homogeneous model, and we quantify the improvement in terms of accuracy, linearity, and cross-talk from extracerebral layers.

Effects of spatial pattern scale of brain activity on the sensitivity of DOT, fMRI, EEG and MEG

The objective of this work is to quantify how patterns of cortical activity at different spatial scales are measured by noninvasive functional neuroimaging sensors. We simulated cortical activation patterns at nine different spatial scales in a realistic head model and propagated this activity to magnetoencephalography (MEG), electroencephalography (EEG), diffuse optical tomography (DOT), and functional magnetic resonance imaging (fMRI) sensors in arrangements that are typically used in functional neuroimaging studies. We estimated contrast transfer functions (CTF), correlation distances in sensor space, and the minimum resolvable spatial scale of cortical activity for each modality. We found that CTF decreases as the spatial extent of cortical activity decreases, and that correlations between nearby sensors depend on the spatial extent of cortical activity. For cortical activity on the intermediate spatial scale of 6.7 cm(2), the correlation distances (r>0.5) were 1.0 cm for fMRI, 2.0 cm for DOT, 12.8 for EEG, 9.5 cm for MEG magnetometers and 9.7 cm for MEG gradiometers. The resolvable spatial pattern scale was found to be 1.43 cm(2) for MEG magnetometers, 0.88 cm(2) for MEG gradiometers, 376 cm(2) for EEG, 0.75 cm(2) for DOT, and 0.072 cm(2) for fMRI. These findings show that sensitivity to cortical activity varies substantially as a function of spatial scale within and between the different imaging modalities. This information should be taken into account when interpreting neuroimaging data and when choosing the number of nodes for network analyses in sensor space.

Spatial registration for functional near-infrared spectroscopy: from channel position on the scalp to cortical location in individual and group analyses

Functional near-infrared spectroscopy (fNIRS) has now become widely accepted as a common functional imaging modality. In order for fNIRS to achieve genuine neuroimaging citizenship, it would ideally be equipped with functional and structural image analyses. However, fNIRS measures cortical activities from the head surface without anatomical information of the object being measured. In this review article, we will present a methodological overview of spatial registration of fNIRS data to overcome this technical drawback of fNIRS. We first introduce and explore the use of standard stereotaxic space and anatomical labeling. Second, we explain different ways of describing scalp landmarks using 10-20 based systems. Third, we describe the simplest case of fNIRS data co-registration to a subject's own MRI. Fourth, we extend the concept to fNIRS data registration of group data. Fifth, we describe probabilistic registration methods, which use a reference-MRI database instead of a subject's own MRIs, and thus enable MRI-free registration for standalone fNIRS data. Sixth, we further extend the concept of probabilistic registration to three-dimensional image reconstruction in diffuse optical tomography. Seventh, we describe a 3D-digitizer-free method for the virtual registration of fNIRS data. Eighth, we provide practical guidance on how these techniques are implemented in software. Finally, we provide information on current resources and limitations for spatial registration of child and infant data. Through these technical descriptions, we stress the importance of presenting fNIRS data on a common platform to facilitate both intra- and inter-modal data sharing among the neuroimaging community.

Depth-compensated diffuse optical tomography enhanced by general linear model analysis and an anatomical atlas of human head

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.

Reducing motion artifacts for long-term clinical NIRS monitoring using collodion-fixed prism-based optical fibers

As the applications of near-infrared spectroscopy (NIRS) continue to broaden and long-term clinical monitoring becomes more common, minimizing signal artifacts due to patient movement becomes more pressing. This is particularly true in applications where clinically and physiologically interesting events are intrinsically linked to patient movement, as is the case in the study of epileptic seizures. In this study, we apply an approach common in the application of EEG electrodes to the application of specialized NIRS optical fibers. The method provides improved optode-scalp coupling through the use of miniaturized optical fiber tips fixed to the scalp using collodion, a clinical adhesive. We investigate and quantify the performance of this new method in minimizing motion artifacts in healthy subjects, and apply the technique to allow continuous NIRS monitoring throughout epileptic seizures in two epileptic in-patients. Using collodion-fixed fibers reduces the percent signal change of motion artifacts by 90% and increases the SNR by 6 and 3 fold at 690 and 830 nm wavelengths respectively when compared to a standard Velcro-based array of optical fibers. The SNR has also increased by 2 fold during rest conditions without motion with the new probe design because of better light coupling between the fiber and scalp. The change in both HbO and HbR during motion artifacts is found to be statistically lower for the collodion-fixed fiber probe. The collodion-fixed optical fiber approach has also allowed us to obtain good quality NIRS recording of three epileptic seizures in two patients despite excessive motion in each case.

Atlas-based head modeling and spatial normalization for high-density diffuse optical tomography: in vivo validation against fMRI

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.