Complex wavelets applied to diffuse optical spectroscopy for brain activity detection

Evaluation of fNIRS signal components elicited by cognitive and hypercapnic stimuli

Functional near infrared spectroscopy (fNIRS) measurements are confounded by signal components originating from multiple physiological causes, whose activities may vary temporally and spatially (across tissue layers, and regions of the cortex). Furthermore, the stimuli can induce evoked effects, which may lead to over or underestimation of the actual effect of interest. Here, we conducted a temporal, spectral, and spatial analysis of fNIRS signals collected during cognitive and hypercapnic stimuli to characterize effects of functional versus systemic responses. We utilized wavelet analysis to discriminate physiological causes and employed long and short source-detector separation (SDS) channels to differentiate tissue layers. Multi-channel measures were analyzed further to distinguish hemispheric differences. The results highlight cardiac, respiratory, myogenic, and very low frequency (VLF) activities within fNIRS signals. Regardless of stimuli, activity within the VLF band had the largest contribution to the overall signal. The systemic activities dominated the measurements from the short SDS channels during cognitive stimulus, but not hypercapnic stimulus. Importantly, results indicate that characteristics of fNIRS signals vary with type of the stimuli administered as cognitive stimulus elicited variable responses between hemispheres in VLF band and task-evoked temporal effect in VLF, myogenic and respiratory bands, while hypercapnic stimulus induced a global response across both hemispheres.

Short-channel regression in functional near-infrared spectroscopy is more effective when considering heterogeneous scalp hemodynamics

Significance: The reliability of functional near-infrared spectroscopy (fNIRS) measurements is reduced by systemic physiology. Short-channel regression algorithms aim at removing systemic "noise" by subtracting the signal measured at a short source-detector separation (mainly scalp hemodynamics) from the one of a long separation (brain and scalp hemodynamics). In literature, incongruent approaches on the selection of the optimal regressor signal are reported based on different assumptions on scalp hemodynamics properties. Aim: We investigated the spatial and temporal distribution of scalp hemodynamics over the sensorimotor cortex and evaluated its influence on the effectiveness of short-channel regressions. Approach: We performed hand-grasping and resting-state experiments with five subjects, measuring with 16 optodes over sensorimotor areas, including eight 8-mm channels. We performed detailed correlation analyses of scalp hemodynamics and evaluated 180 hand-grasping and 270 simulated (overlaid on resting-state measurements) trials. Five short-channel regressor combinations were implemented with general linear models. Three were chosen according to literature, and two were proposed based on additional physiological assumptions [considering multiple short channels and their Mayer wave (MW) oscillations]. Results: We found heterogeneous hemodynamics in the scalp, coming on top of a global close-to-homogeneous behavior (correlation 0.69 to 0.92). The results further demonstrate that short-channel regression always improves brain activity estimates but that better results are obtained when heterogeneity is assumed. In particular, we highlight that short-channel regression is more effective when combining multiple scalp regressors and when MWs are additionally included. Conclusion: We shed light on the selection of optimal regressor signals for improving the removal of systemic physiological artifacts in fNIRS. We conclude that short-channel regression is most effective when assuming heterogeneous hemodynamics, in particular when combining spatial- and frequency-specific information. A better understanding of scalp hemodynamics and more effective short-channel regression will promote more accurate assessments of functional brain activity in clinical and research settings.

Effects of Processing Methods on fNIRS Signals Assessed During Active Walking Tasks in Older Adults

Functional near infrared spectroscopy (fNIRS) is a noninvasive optics-based neuroimaging modality successfully applied to real-life settings. The technology uses light in the near infrared range (650-950nm) to track changes in oxygenated (HbO2) and deoxygenated hemoglobin (Hb) obtained from measured light intensity using light-tissue interaction principles. fNIRS data processing involves artifact removal and hemodynamic signal conversion using modified Beer-Lambert law (MBLL) to obtain Hb and HbO2, reliably. fNIRS signals can get contaminated by various noise sources of physiological and non-physiological origins. Various algorithms have been proposed for the elimination of artifacts from frequency selective filters to blind source separation methods. Hemodynamic signal extraction using raw intensity measurements at multiple wavelengths based on MBLL usually involves apriori knowledge of certain conversion parameters such as molar extinction coefficients ( ε ) and differential path length factor (DPF). Different sets of conversion parameters dependent upon wavelength, chromophores, and age have been reported. Variation in processing algorithms and parameters can cause differences in Hb and HbO2 extraction which can in turn change study outcomes. Using fNIRS, we have previously shown significant increases in oxygenation in the prefrontal cortex from Single-Task-Walking (STW) to Dual-task-Walking (DTW) conditions in older adults due to greater cognitive demands inherent in the latter condition. In the current study, we re-analyzed our data and determined that although using different conversion parameters i.e. ε and age dependent DPF and filter cut-off frequencies at 0.14 and 0.08Hz including those designed to remove confounding effects of Mayer waves had caused some linear increases or decreases on the extracted Hb and HbO2 signals, those effects were minimal in task related comparisons and hence, the overall study outcomes.

Shining a Light on Awareness: A Review of Functional Near-Infrared Spectroscopy for Prolonged Disorders of Consciousness

Qualitative clinical assessments of the recovery of awareness after severe brain injury require an assessor to differentiate purposeful behavior from spontaneous behavior. As many such behaviors are minimal and inconsistent, behavioral assessments are susceptible to diagnostic errors. Advanced neuroimaging tools can bypass behavioral responsiveness and reveal evidence of covert awareness and cognition within the brains of some patients, thus providing a means for more accurate diagnoses, more accurate prognoses, and, in some instances, facilitated communication. The majority of reports to date have employed the neuroimaging methods of functional magnetic resonance imaging, positron emission tomography, and electroencephalography (EEG). However, each neuroimaging method has its own advantages and disadvantages (e.g., signal resolution, accessibility, etc.). Here, we describe a burgeoning technique of non-invasive optical neuroimaging—functional near-infrared spectroscopy (fNIRS)—and review its potential to address the clinical challenges of prolonged disorders of consciousness. We also outline the potential for simultaneous EEG to complement the fNIRS signal and suggest the future directions of research that are required in order to realize its clinical potential.

Cerebral oxygen saturation and peripheral perfusion in the extremely premature infant with intraventricular and/or pulmonary haemorrhage early in life

Extremely preterm infants are at higher risk of pulmonary (PH) and intraventricular (IVH) haemorrhage during the transitioning physiology due to immature cardiovascular system. Monitoring of haemodynamics can detect early abnormal circulation that may lead to these complications. We described time-frequency relationships between near infrared spectroscopy (NIRS) cerebral regional haemoglobin oxygen saturation (CrSO₂) and preductal peripheral perfusion index (PI), capillary oxygen saturation (SpO₂) and heart rate (HR) in extremely preterm infants in the first 72 h of life. Patients were sub-grouped in infants with PH and/or IVH (N _H  = 8) and healthy controls (N _C  = 11). Data were decomposed in wavelets allowing the analysis of localized variations of power. This approach allowed to quantify the percentage of time of significant cross-correlation, semblance, gain (transfer function) and coherence between signals. Ultra-low frequencies (<0.28 mHz) were analyzed as slow and prolonged periods of impaired circulation are considered more detrimental than transient fluctuations. Cross-correlation between CrSO₂ and oximetry (PI, SpO₂ and HR) as well as in-phase semblance and gain between CrSO₂ and HR were significantly lower while anti-phase semblance between CrSO₂ and HR was significantly higher in PH-IVH infants compared to controls. These differences may reflect haemodynamic instability associated with cerebrovascular autoregulation and hemorrhagic complications observed during the transitioning physiology.

Mayer waves reduce the accuracy of estimated hemodynamic response functions in functional near-infrared spectroscopy

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.

Short separation regression improves statistical significance and better localizes the hemodynamic response obtained by near-infrared spectroscopy for tasks with differing autonomic responses

Autonomic nervous system response is known to be highly task-dependent. The sensitivity of near-infrared spectroscopy (NIRS) measurements to superficial layers, particularly to the scalp, makes it highly susceptible to systemic physiological changes. Thus, one critical step in NIRS data processing is to remove the contribution of superficial layers to the NIRS signal and to obtain the actual brain response. This can be achieved using short separation channels that are sensitive only to the hemodynamics in the scalp. We investigated the contribution of hemodynamic fluctuations due to autonomous nervous system activation during various tasks. Our results provide clear demonstrations of the critical role of using short separation channels in NIRS measurements to disentangle differing autonomic responses from the brain activation signal of interest.

Evoked hemodynamic response estimation using ensemble empirical mode decomposition based adaptive algorithm applied to dual channel functional near infrared spectroscopy (fNIRS)

BACKGROUND

The quality of the functional near infrared spectroscopy (fNIRS) recordings is highly degraded by the presence of physiological interferences. It is crucial to efficiently separate the evoked hemodynamic responses (EHRs) from other background hemodynamic activities prior to any further processing.

NEW METHOD

This paper presents a novel algorithm for physiological interferences reduction from the dual channel fNIRS measurements using ensemble empirical mode decomposition (EEMD) technique. The proposed algorithm is comprised of two main steps: (1) decomposing reference signal into its constituents called intrinsic mode functions (IMFs) and (2) adaptively defining appropriate weights of the corresponding IMFs to estimate the proportion of physiological interference in standard channel measurement.

RESULTS

Performance of the proposed algorithm was evaluated using both synthetic and semi-real brain hemodynamic data based on four parameters of relative mean squared error (rMSE), Pearson's correlation coefficient (R(2)), percentage estimation error of peak amplitude (EPA) and peak latency (EL).

COMPARISON WITH EXISTING METHODS

Results obtained from synthetic data revealed that both the EEMD based normalized least mean squares (EEMD-NLMS) and EEMD based recursive least squares (EEMD-RLS) methods could reduce the average rMSE by at least 34% and 49%, respectively, when compared with widely used methods: block averaging, band-pass filtering and principal and/or independent component analysis. Furthermore, the two proposed methods outperform the regression method in reducing rMSE by at least 21% and 35% respectively when applied to semi-real data.

CONCLUSIONS

An effective algorithm for estimating the EHRs from raw fNIRS data was proposed in which no assumption about the amplitude, shape and duration of the responses is considered.

Short separation channel location impacts the performance of short channel regression in NIRS

Near-Infrared Spectroscopy (NIRS) allows the recovery of cortical oxy- and deoxyhemoglobin changes associated with evoked brain activity. NIRS is a back-reflection measurement making it very sensitive to the superficial layers of the head, i.e. the skin and the skull, where systemic interference occurs. As a result, the NIRS signal is strongly contaminated with systemic interference of superficial origin. A recent approach to overcome this problem has been the use of additional short source-detector separation optodes as regressors. Since these additional measurements are mainly sensitive to superficial layers in adult humans, they can be used to remove the systemic interference present in longer separation measurements, improving the recovery of the cortical hemodynamic response function (HRF). One question that remains to answer is whether or not a short separation measurement is required in close proximity to each long separation NIRS channel. Here, we show that the systemic interference occurring in the superficial layers of the human head is inhomogeneous across the surface of the scalp. As a result, the improvement obtained by using a short separation optode decreases as the relative distance between the short and the long measurement is increased. NIRS data was acquired on 6 human subjects both at rest and during a motor task consisting of finger tapping. The effect of distance between the short and the long channel was first quantified by recovering a synthetic hemodynamic response added over the resting-state data. The effect was also observed in the functional data collected during the finger tapping task. Together, these results suggest that the short separation measurement must be located as close as 1.5 cm from the standard NIRS channel in order to provide an improvement which is of practical use. In this case, the improvement in Contrast-to-Noise Ratio (CNR) compared to a standard General Linear Model (GLM) procedure without using any small separation optode reached 50% for HbO and 100% for HbR. Using small separations located farther than 2 cm away resulted in mild or negligible improvements only.

The utility of near-infrared spectroscopy in the regression of low-frequency physiological noise from functional magnetic resonance imaging data

Near-infrared spectroscopy (NIRS) signals have been shown to correlate with resting-state BOLD-fMRI data across the whole brain volume, particularly at frequencies below 0.1Hz. While the physiological origins of this correlation remain unclear, its existence may have a practical application in minimizing the background physiological noise present in BOLD-fMRI recordings. We performed simultaneous, resting-state fMRI and 28-channel NIRS in seven adult subjects in order to assess the utility of NIRS signals in the regression of physiological noise from fMRI data. We calculated the variance of the residual error in a general linear model of the baseline fMRI signal, and the reduction of this variance achieved by including NIRS signals in the model. In addition, we introduced a sequence of simulated hemodynamic response functions (HRFs) into the resting-state fMRI data of each subject in order to quantify the effectiveness of NIRS signals in optimizing the recovery of that HRF. For comparison, these calculations were also performed using a pulse and respiration RETROICOR model. Our results show that the use of 10 or more NIRS channels can reduce variance in the residual error by as much as 36% on average across the whole cortex. However the same number of low-pass filtered white noise regressors is shown to produce a reduction of 19%. The RETROICOR model obtained a variance reduction of 6.4%. Our HRF simulation showed that the mean-squared error (MSE) between the recovered and true HRFs is reduced by 21% on average when 10 NIRS channels are applied and by introducing an optimized time lag between the NIRS and fMRI time series, a single NIRS channel can provide an average MSE reduction of 14%. The RETROICOR model did not provide a significant change in MSE. By each of the metrics calculated, NIRS recording is shown to be of significant benefit to the regression of low-frequency physiological noise from fMRI data.

Improved recovery of the hemodynamic response in diffuse optical imaging using short optode separations and state-space modeling

Diffuse optical imaging (DOI) allows the recovery of the hemodynamic response associated with evoked brain activity. The signal is contaminated with systemic physiological interference which occurs in the superficial layers of the head as well as in the brain tissue. The back-reflection geometry of the measurement makes the DOI signal strongly contaminated by systemic interference occurring in the superficial layers. A recent development has been the use of signals from small source-detector separation (1cm) optodes as regressors. Since those additional measurements are mainly sensitive to superficial layers in adult humans, they help in removing the systemic interference present in longer separation measurements (3 cm). Encouraged by those findings, we developed a dynamic estimation procedure to remove global interference using small optode separations and to estimate simultaneously the hemodynamic response. The algorithm was tested by recovering a simulated synthetic hemodynamic response added over baseline DOI data acquired from 6 human subjects at rest. The performance of the algorithm was quantified by the Pearson R(2) coefficient and the mean square error (MSE) between the recovered and the simulated hemodynamic responses. Our dynamic estimator was also compared with a static estimator and the traditional adaptive filtering method. We observed a significant improvement (two-tailed paired t-test, p<0.05) in both HbO and HbR recovery using our Kalman filter dynamic estimator compared to the traditional adaptive filter, the static estimator and the standard GLM technique.

Quantitative investigation of the effect of the extra-cerebral vasculature in diffuse optical imaging: a simulation study

Diffuse optical imaging (DOI) is a non invasive technique allowing the recovery of hemodynamic changes in the brain. Due to the diffusive nature of photon propagation in turbid media and the fact that cerebral tissues are located around 1.5 cm under the adult human scalp, DOI measurements are subject to partial volume errors. DOI measurements are also sensitive to large pial vessels because oxygenated and deoxygenated hemoglobin are the dominant chromophores in the near infrared window. In this study, the effect of the extra-cerebral vasculature in proximity of the sagittal sinus was investigated for its impact on DOI measurements simulated over the human adult visual cortex. Numerical Monte Carlo simulations were performed on two specific models of the human head derived from magnetic resonance imaging (MRI) scans. The first model included the extra-cerebral vasculature in which constant hemoglobin concentrations were assumed while the second did not. The screening effect of the vasculature was quantified by comparing recovered hemoglobin changes from each model for different optical arrays and regions of activation. A correction factor accounting for the difference between the recovered and the simulated hemoglobin changes was computed in each case. The results show that changes in hemoglobin concentration are better estimated when the extra-cerebral vasculature is modeled and the correction factors obtained in this case were at least 1.4-fold lower. The effect of the vasculature was also examined in a high-density diffuse optical tomography configuration. In this case, the difference between changes in hemoglobin concentration recovered with each model was reduced down to 10%.

1/f noise in diffuse optical imaging and wavelet-based response estimation

In diffuse optical imaging (DOI) data analysis, the functional response is contaminated with physiological noise as in functional magnetic resonance imaging (fMRI). In this work we extend a previously proposed method for fMRI to estimate the parameters of a linear model of DOI time series. The regression is performed in the wavelet domain to infer drift coefficients at different scales and to estimate the strength of the hemodynamic response function (HRF). This multiresolution approach benefits from the whitening property of the discrete wavelet transform (DWT), which approximately decorrelates long-memory noise processes. We also show that a more accurate estimation is obtained by removing some regressors correlating with the protocol. Moreover, we observe that this improvement is related to a quantitative measure of 1/f noise. The performances of the method are first evaluated against a standard spline-cosine drift approach with simulated HRF and real background physiology. Lastly, the technique is applied to experimental event-related data acquired by near-infrared spectroscopy (NIRS).