A methodology to improve estimation of stimulus-evoked hemodynamic response from fNIRS measurements.

Functional near-infrared spectroscopy (fNIRS) is a non-invasive optical neuroimaging method used to investigate functional activity of the cerebral cortex evoked by cognitive, visual, auditory and motor tasks, detecting regional changes of oxy- and deoxy-hemoglobin concentration. Accurate estimation of the stimulus-evoked hemodynamic response (HR) from fNIRS signals in order to quantitatively investigate cognitive functions requires to cope with several noise components. Some of them appear as random disturbances (typically tackled through averaging techniques), while others are due to physiological sources, such as heart beat, respiration, vasomotor waves, and are particularly challenging to be dealt with because they lie in the same frequency band of HR. In this work we present a new two-steps methodology for the HR estimation from fNIRS data. The first step is a pre-processing stage where physiological trends in fNIRS data are reduced by exploiting a mathematical model identified from the signal of a reference channel. In the second step, the pre-processed data of the other channels are filtered with a recently presented non-parametric Bayesian approach (Scarpa et al., Optics Express, 2010). The presented method for HR estimation is compared with widely used methods: conventional averaging, band-pass filtering and principal component analysis (PCA). Results on simulated data reveal the ability of the proposed method to improve the accuracy of the estimates of the functional hemodynamic response, as well as the estimate of peak amplitude and latency. Encouraging preliminary results in a representative real data set showing an improvement of contrast to noise ratio are also reported.

Bayesian filtering of human brain hemodynamic activity elicited by visual short-term maintenance recorded through functional near-infrared spectroscopy (fNIRS).

Functional near-infrared spectroscopy (fNIRS) is a neuroimaging technique that measures changes in oxy-hemoglobin (ΔHbO) and deoxy-hemoglobin (ΔHbR) concentration associated with brain activity. The signal acquired with fNIRS is naturally affected by disturbances engendering from ongoing physiological activity (e.g., cardiac, respiratory, Mayer wave) and random measurement noise. Despite its several drawbacks, the so-called conventional averaging (CA) is still widely used to estimate the hemodynamic response function (HRF) from noisy signal. One such drawback is related to the number of trials necessary to derive stable HRF functions adopting the CA approach, which must be substantial (N >> 50). In this work, a pre-processing procedure to remove artifacts followed by the application of a non-parametric Bayesian approach is proposed that capitalizes on a priori available knowledge about HRF and noise. Results with the proposed Bayesian approach were compared with CA and with a straightforward band-pass filtering approach. On simulated data, a five times lower estimation error on HRF was obtained with respect to that obtained by CA, and 2.5 times lower than that obtained by band pass filtering. On real data, the improvement achieved by the present method was attested by an increase in the contrast to noise ratio (CNR) and by a reduced variability in single trial estimation. An application of the present Bayesian approach is illustrated that was optimized to monitor changes in hemodynamic activity reflecting variations in visual short-term memory load in humans, which are notoriously hard to detect using functional magnetic resonance imaging (fMRI). In particular, statistical analyses of HRFs recorded during a memory task established with high reliability the crucial role of the intraparietal sulcus and the intra-occipital sulcus in posterior areas of the human brain in visual short-term memory maintenance.