Effect of probe arrangement on reproducibility of images by near-infrared topography evaluated by a virtual head phantom

Diffuse optical reconstructions of functional near infrared spectroscopy data using maximum entropy on the mean

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\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {cm}^2$$\end{document}cm2 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.

Diffuse Optical Tomography Using fNIRS Signals Measured from the Skull Surface of the Macaque Monkey

Diffuse optical tomography (DOT), as a functional near-infrared spectroscopy (fNIRS) technique, can estimate three-dimensional (3D) images of the functional hemodynamic response in brain volume from measured optical signals. In this study, we applied DOT algorithms to the fNIRS data recorded from the surface of macaque monkeys’ skulls when the animals performed food retrieval tasks using either the left- or right-hand under head-free conditions. The hemodynamic response images, reconstructed by DOT with a high sampling rate and fine voxel size, demonstrated significant activations at the upper limb regions of the primary motor area in the central sulcus and premotor, and parietal areas contralateral to the hands used in the tasks. The results were also reliable in terms of consistency across different recording dates. Time-series analyses of each brain area revealed preceding activity of premotor area to primary motor area consistent with previous physiological studies. Therefore, the fNIRS–DOT protocol demonstrated in this study provides reliable 3D functional brain images over a period of days under head-free conditions for region-of-interest–based time-series analysis.

Exclusive detection of cerebral hemodynamics in functional near-infrared spectroscopy by reflectance modulation of the scalp surface

SIGNIFICANCE

Functional near-infrared spectroscopy (fNIRS) is a technique for detecting regional hemodynamic responses associated with neural activation in the cerebral cortex. The absorption changes due to hemodynamic changes in the scalp cause considerable signal contamination in the fNIRS measurement. A method for extracting hemodynamic changes in the cerebral tissue is required for reliable fNIRS measurement.

AIM

To exclusively detect cerebral functional hemodynamic changes, we developed an fNIRS technique using reflectance modulation of the scalp surface.

APPROACH

The theoretical feasibility of the proposed method was proven by a simulation calculation of light propagation. Its practical feasibility was evaluated by a phantom experiment and brain activation simulation mimicking human fNIRS experiments.

RESULTS

The simulation calculation revealed that the partial path length of the scalp was changed by reflectance modulation of the scalp surface. The influence of absorption change in the superficial layer was successfully reduced by the proposed method, using only measurement data, in the phantom experiment. The proposed method was applicable to human experiments of standard designs, achieving statistical significance within an acceptable experimental time-frame.

CONCLUSIONS

Removal of the scalp hemodynamic effect by the proposed technique will increase the quality of fNIRS data, particularly in measurements in neonates and infants that typically would require a dense optode arrangement.

Optimal positioning of optodes on the scalp for personalized functional near-infrared spectroscopy investigations

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.

Functional near-infrared spectroscopy for monitoring macaque cerebral motor activity during voluntary movements without head fixation

We developed an fNIRS system for monitoring macaque cerebral motor activity during voluntary movements without head fixation. fNIRS data at 27 channels in 7.5 mm spatial interval were calibrated by simulating light propagation through the macaque cranial tissues. The subject was instructed to repeatedly (75 times) retrieve a food pellet with alternating left or right hands from a food well for each session. We detected significant increases in oxygenated hemoglobin (Hb) and decrease in deoxygenated Hb in the primary motor area (M1) contralateral to the hand used. In more rostral and ventral regions in both hemispheres, the hemodynamic similarly changed regardless of used hand. Direct feeding to the mouth eliminated activity in the hand M1 whereas that at bilateral ventral regions (mouth M1 area) remained. Statistical analyses for the hemodynamics between left/right-hand use revealed the location of each hand M1 in either hemisphere. In these regions, the maximum amplitude and time of the maximum amplitude in the hemodynamic response evoked by food retrieval were highly correlated with the time associated with food retrieval. We could assign each channel to an appropriate functional motor area, providing proof of principle for future studies involving brain damage models in freely moving macaque monkeys.

Magnetic resonance imaging appropriate for construction of subject-specific head models for diffuse optical tomography

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.

Optimal optode montage on electroencephalography/functional near-infrared spectroscopy caps dedicated to study epileptic discharges

Functional near-infrared spectroscopy (fNIRS), acquired simultaneously with electroencephalography (EEG), allows the investigation of hemodynamic brain responses to epileptic activity. Because the presumed epileptogenic focus is patient-specific, an appropriate source/detector (SD) montage has to be reconfigured for each patient. The combination of EEG and fNIRS, however, entails several constraints on montages, and finding an optimal arrangement of optodes on the cap is an important issue. We present a method for computing an optimal SD montage on an EEG/fNIRS cap that focuses on one or several specific brain regions; the montage maximizes the spatial sensitivity. We formulate this optimization problem as a linear integer programming problem. The method was evaluated on two EEG/fNIRS caps. We simulated absorbers at different locations on a head model and generated realistic optical density maps on the scalp. We found that the maps of optimal SD montages had spatial resolution properties comparable to those of regular SD arrangements for the whole head with significantly fewer sensors than regular SD arrangements. In addition, we observed that optimal montages yielded improved spatial density of fNIRS measurements over the targeted regions together with an increase in signal-to-noise ratio.

Direct cortical hemodynamic mapping of somatotopy of pig nostril sensation by functional near-infrared cortical imaging (fNCI)

Functional near-infrared spectroscopy (fNIRS) is a neuroimaging technique for the noninvasive monitoring of human brain activation states utilizing the coupling between neural activity and regional cerebral hemodynamics. Illuminators and detectors, together constituting optodes, are placed on the scalp, but due to the presence of head tissues, an inter-optode distance of more than 2.5cm is necessary to detect cortical signals. Although direct cortical monitoring with fNIRS has been pursued, a high-resolution visualization of hemodynamic changes associated with sensory, motor and cognitive neural responses directly from the cortical surface has yet to be realized. To acquire robust information on the hemodynamics of the cortex, devoid of signal complications in transcranial measurement, we devised a functional near-infrared cortical imaging (fNCI) technique. Here we demonstrate the first direct functional measurement of temporal and spatial patterns of cortical hemodynamics using the fNCI technique. For fNCI, inter-optode distance was set at 5mm, and light leakage from illuminators was prevented by a special optode holder made of a light-shielding rubber sheet. fNCI successfully detected the somatotopy of pig nostril sensation, as assessed in comparison with concurrent and sequential somatosensory-evoked potential (SEP) measurements on the same stimulation sites. Accordingly, the fNCI system realized a direct cortical hemodynamic measurement with a spatial resolution comparable to that of SEP mapping on the rostral region of the pig brain. This study provides an important initial step toward realizing functional cortical hemodynamic monitoring during neurosurgery of human brains.

Exploration of cerebral activation using hemodynamic modality separation method in high-density multichannel fNIRS

Hemodynamic modality separation (HMS) is a method for separating the functional near infrared spectroscopy (fNIRS) signal into the cerebral functional and systemic physiological components based on their difference in hemodynamic modalities: 1) Changes in oxyhemoglobin and deoxyhemoglobin (ΔHbO and ΔHbR) in the cerebral capillaries during neural activation negatively correlate with each other; 2) Other physiological hemodynamic changes originating from major vessels cause a positive correlation in ΔHbO and ΔHbR. We applied this simple method to a high-density multichannel (HDM) fNIRS measurement. In the case of functional signal detection in the parietal area of human adults during a single-sided finger-tapping task, conventional fNIRS data showed very unclear signal laterality, while the functional components separated by the HMS method highly localized at the contralateral area of the tapping side. Using the HMS method for HDM NIRS, we successfully explored cerebral activation in the parietal area. This is the first report that HMS method was utilized for the exploratory detection of cerebral activity.

Separation of fNIRS signals into functional and systemic components based on differences in hemodynamic modalities

In conventional functional near-infrared spectroscopy (fNIRS), systemic physiological fluctuations evoked by a body's motion and psychophysiological changes often contaminate fNIRS signals. We propose a novel method for separating functional and systemic signals based on their hemodynamic differences. Considering their physiological origins, we assumed a negative and positive linear relationship between oxy- and deoxyhemoglobin changes of functional and systemic signals, respectively. Their coefficients are determined by an empirical procedure. The proposed method was compared to conventional and multi-distance NIRS. The results were as follows: (1) Nonfunctional tasks evoked substantial oxyhemoglobin changes, and comparatively smaller deoxyhemoglobin changes, in the same direction by conventional NIRS. The systemic components estimated by the proposed method were similar to the above finding. The estimated functional components were very small. (2) During finger-tapping tasks, laterality in the functional component was more distinctive using our proposed method than that by conventional fNIRS. The systemic component indicated task-evoked changes, regardless of the finger used to perform the task. (3) For all tasks, the functional components were highly coincident with signals estimated by multi-distance NIRS. These results strongly suggest that the functional component obtained by the proposed method originates in the cerebral cortical layer. We believe that the proposed method could improve the reliability of fNIRS measurements without any modification in commercially available instruments.

The influence of frontal sinus in brain activation measurements by near-infrared spectroscopy analyzed by realistic head models

Adequate modeling of light propagation in the head is important to predict the sensitivity of NIRS signal and the spatial sensitivity profile of source-detector pairs. The 3D realistic head models of which the geometry is based upon the anatomical images acquired by magnetic resonance imaging and x-ray computed tomography are constructed to investigate the influence of the frontal sinus on the NIRS signal and spatial sensitivity. Light propagation in the head is strongly affected by the presence of the frontal sinus. The light tends to propagate around the frontal sinus. The influence of the frontal sinus on the sensitivity of the NIRS signal to the brain activation is not consistent and depends on the depth of the frontal sinus, the optical properties of the superficial tissues and the relative position between the source-detector pair and the frontal sinus. The frontal sinus located in the shallow region of the skull tends to reduce the sensitivity of the NIRS signal while the deep frontal sinus can increase the sensitivity of the NIRS signal.

Development of double density whole brain fNIRS with EEG system for brain machine interface

Brain-machine interfaces (BMI) are expected as new man-machine interfaces. Non-invasive BMI have the potential to improve the quality of life of many disabled individuals with safer operation. The non-invasive BMI using the functional functional near-infrared spectroscopy (fNIRS) with the electroencephalogram (EEG) has potential applicability beyond the restoration of lost movement and rehabilitation in paraplegics and would enable normal individuals to have direct brain control of external devices in their daily lives. To shift stage of the non-invasive BMI from laboratory to clinical, the key factor is to develop high-accuracy signal decoding technology and highly restrictive of the measurement area. In this article, we present the development of a high-accuracy brain activity measurement system by combining fNIRS and EEG. The new fNIRS had high performances with high spatial resolution using double density technique and a large number of measurement channels to cover a whole human brain.

Simulation of near-infrared light absorption considering individual head and prefrontal cortex anatomy: implications for optical neuroimaging

Functional near-infrared spectroscopy (fNIRS) is an established optical neuroimaging method for measuring functional hemodynamic responses to infer neural activation. However, the impact of individual anatomy on the sensitivity of fNIRS measuring hemodynamics within cortical gray matter is still unknown. By means of Monte Carlo simulations and structural MRI of 23 healthy subjects (mean age: 25.0±2.8 years), we characterized the individual distribution of tissue-specific NIR-light absorption underneath 24 prefrontal fNIRS channels. We, thereby, investigated the impact of scalp-cortex distance (SCD), frontal sinus volume as well as sulcal morphology on gray matter volumes (V(gray)) traversed by NIR-light, i.e. anatomy-dependent fNIRS sensitivity. The NIR-light absorption between optodes was distributed describing a rotational ellipsoid with a mean penetration depth of (23.6±0.7) mm considering the deepest 5% of light. Of the detected photon packages scalp and bone absorbed (96.4±9.7)% and V(gray) absorbed (3.1±1.8)% of the energy. The mean V(gray) volume (1.1±0.4) cm3 was negatively correlated (r=-.76) with the SCD and frontal sinus volume (r=-.57) and was reduced by 41.5% in subjects with relatively large compared to small frontal sinus. Head circumference was significantly positively correlated with the mean SCD (r=.46) and the traversed frontal sinus volume (r=.43). Sulcal morphology had no significant impact on V(gray). Our findings suggest to consider individual SCD and frontal sinus volume as anatomical factors impacting fNIRS sensitivity. Head circumference may represent a practical measure to partly control for these sources of error variance.

Structural atlas-based spatial registration for functional near-infrared spectroscopy enabling inter-study data integration

OBJECTIVE

The use of functional near-infrared spectroscopy (fNIRS) is growing, leading to a need for methods to summarise data from multiple studies. However, this is difficult using the current channel-based methods when experiments do not share a common channel (CH) arrangement. Thus, we proposed and implemented a CH-independent analysis method for summarising fNIRS data.

METHODS

We defined sub-regions as spatial bins to organise fNIRS data. Sub-regions were defined on the standard brain surface based on macro- and micro-structural information. After probabilistically estimating CH location in standard stereotaxic brain space, the CH-based data were reorganised into these spatial bins to evaluate sub-region-based activation.

RESULTS

Sub-regions with sizes corresponding to fNIRS spatial resolution were defined. We demonstrated this method by integrating data from two of our fNIRS studies that shared the same region of interest but used different channel arrangements.

CONCLUSIONS

Using this method, data from multiple fNIRS studies with different CH arrangements can be integrated in standard brain space, while keeping in mind the brain structure-function relationship.

SIGNIFICANCE

The current method will facilitate an effective use of accumulating fNIRS data by allowing integration of data from multiple studies.

Significance of background optical properties, time-resolved information and optode arrangement in diffuse optical imaging of term neonates

The significance of accurate knowledge of background optical properties and time-resolved information in reconstructing images of hemodynamic changes in the neonatal brain from diffuse optical imaging data was studied using Monte Carlo (MC) simulation. A segmented anatomical magnetic resonance (MR) image and literature-derived optical properties for each tissue type were used to create a voxel-based anatomical model. Small absorbing perturbations were introduced into the anatomical model to simulate localized hemodynamic responses related to brain activation. Perturbation MC (pMC) was used as the primary method of image reconstruction. For comparison, reconstructions were also performed using the finite element method (FEM) to solve the diffusion approximation (DA) to the radiative transfer equation (RTE). The effect of optode layout was investigated using three different grids. Of the factors studied, the density of the optode grid was found to have the greatest effect on image quality. The use of time-resolved information significantly improved the spatial accuracy with all optode grids. Adequate knowledge and modeling of the optical properties of the background was found to significantly improve the spatial accuracy of the reconstructed images and make the recovery of contrast of absorption changes more consistent over simplified modeling. Localization accuracy of small perturbations was found to be 2-3 mm with accurate a priori knowledge of the background optical properties, when a grid with high optode density (>1 optode cm(-2)) was used.