Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Bundled-Optode Method in Functional Near-Infrared Spectroscopy

  • Hoang-Dung Nguyen,

    Affiliation Department of Cogno-Mechatronics Engineering, Pusan National University, 2 Busandaehak-ro, Geumjeong-gu, Busan, 46241, Republic of Korea

  • Keum-Shik Hong ,

    kshong@pusan.ac.kr

    Affiliations Department of Cogno-Mechatronics Engineering, Pusan National University, 2 Busandaehak-ro, Geumjeong-gu, Busan, 46241, Republic of Korea, School of Mechanical Engineering, Pusan National University, 2 Busandaehak-ro, Geumjeong-gu, Busan, 46241, Republic of Korea

  • Yong-Il Shin

    Affiliation Department of Rehabilitation Medicine, School of Medicine, Pusan National University & Research Institute for Convergence of Biomedical Science and Technology, Pusan National University Yangsan Hospital, 20, Geumo-ro, Mulgeum-eup, Yangsan-si, Gyeongsangnam-do, 50612, Republic of Korea

Abstract

In this paper, a theory for detection of the absolute concentrations of oxy-hemoglobin (HbO) and deoxy-hemoglobin (HbR) from hemodynamic responses using a bundled-optode configuration in functional near-infrared spectroscopy (fNIRS) is proposed. The proposed method is then applied to the identification of two fingers (i.e., little and thumb) during their flexion and extension. This experiment involves a continuous-wave-type dual-wavelength (760 and 830 nm) fNIRS and five healthy male subjects. The active brain locations of two finger movements are identified based on the analysis of the t- and p-values of the averaged HbOs, which are quite distinctive. Our experimental results, furthermore, revealed that the hemodynamic responses of two-finger movements are different: The mean, peak, and time-to-peak of little finger movements are higher than those of thumb movements. It is noteworthy that the developed method can be extended to 3-dimensional fNIRS imaging.

Introduction

The objective of this paper is to characterize specific brain regions in the left motor cortex when finger-movements (i.e., flexion and extension) of the right little and thumb fingers are made. For this, the hemodynamic responses (HRs) upon finger movements are measured using functional near-infrared spectroscopy (fNIRS) with a bundled-optode configuration. The activated brain regions (upon 48 channel data) associated with individual finger movements are identified based on t- and p-value analyses of the averaged oxy-hemoglobin (HbO) concentration over five trials.

In the field of fNIRS-based brain-imaging, three techniques are used: continuous-wave (CW), frequency-domain (FD), and time-domain (TD). First, CW-fNIRS, which utilizes constant-intensity light sources of 650~1,000 nm wavelength, is the most popular technique. It measures the scattered light through the brain using light-detectors (optodes) positioned 2~4 cm apart from the light-source on the scalp. Typically, the changes of the hemoglobin concentration can be computed, based on the modified Beer-Lambert law (MBLL), according to the known extinction coefficients and differential path-length factor [114]. Recently, this technique has been integrated into portable equipment using light emitting diodes, electronic circuits, and wireless transceivers [8,15]. Second, FD-fNIRS is a complex technique: It proceeds by modulation of the amplitude of a light source at a very high frequency (i.e., hundreds of megahertz) [1619]. It then measures the absolute values of the HbO and deoxy-hemoglobin (HbR) concentrations based on the changes in the amplitudes and phases of the detected light intensities. The main advantage of the FD approach is the measurement of the absolute values of hemoglobin concentrations. Finally, the TD spectroscopy is an expensive and intricate technique [20,21] that utilizes a very short (on the order of picoseconds) pulsed light to illuminate a specific brain region. This technique allows the measurement of the mean optical path length, absorption coefficient, and scattering coefficient based on the detected back-reflected light. By application of the measured quantitative coefficients, the absolute values of hemoglobin concentrations can then be computed.

The spatially resolved spectroscopy (SRS) approach [6,2225] was used to measure optical properties of tissues. In this approach, one emitter and at least two arranged detectors are used, which means that at least two channels are utilized. The absorption coefficient can then be computed based on the absorbance gradient with respect to the emitter-detector distance. One of the early studies employing this method is Farrell et al. [26]. They measured the diffuse reflectance of light intensities by the distance between an emitter-detector pair (range: 1.0~10.0 mm). Then, from the phantom and muscle tissue data, they estimated the absorption and scattering coefficients using the nonlinear least square fitting approach. Based on the results, the absolute concentration values can be determined by way of the MBLL. Several subsequent studies also focused on the estimation of such coefficients from the measured data of both brain and muscle tissues [20,22,24,27]. In the present work, the SRS method is used with a bundled-optode arrangement (i.e., multiple emitters and multiple detectors arranged side by side). It is noteworthy that short-separation detectors (i.e., 0.5 cm in this paper) as well as long-separation detectors (i.e., 1~4 cm) are used at the same time in our method [28,29]. The superficial noise (i.e., the trend in the scalp) [2830] is reduced by subtraction of the acquired signals between the short-separation detectors from those between the long-separation detectors.

Recently, the HRs in various brain areas (namely, the motor, visual, and somatosensory cortices) have been investigated [3133]. However, finding a specific brain region in a given cortex (e.g., a local region in the motor cortex for a specific finger movement) is difficult. Some previous works have revealed overlapping of the activated sites upon individual finger movements [3438]. One of the early studies on single-finger movement is Schieber and Hibbard [38]. In their work, 10 subjects were asked to extend their right index and middle fingers simultaneously and separately. They found that a separate movement of the fingers was significantly more active than two-finger coincident movement. Moreover, the functional activity increased as the finger-movement rates grew faster [39]. Also, the brain region corresponding to finger extension was more active than that to flexion [40]. Penfield and Rasmussen [41] illustrated the corresponding positions of brain functions in the cerebral cortex (both in the somatosensory and motor areas) by means of an electrical stimulation method (i.e., electrodes placed directly on the patient’s cortices). Recently, the active regions of finger movements in the motor cortex have been located using the functional magnetic resonance imaging technique [36,37,42,43]. Some demonstrated thumb movement to be located in the lateral site and little finger movement in the medial site [35,42,44]. Lately, the fNIRS technique has been utilized to investigate the active locations of the little and thumb fingers in association with the primary somatosensory cortex [45]: In their work, vibro-tactile stimulation was applied to the little and thumb fingers of eight experimental participants. They found that, in five out of eight participants, the active locations of two fingers could be distinguished. Most notably, Yamada et al. [46] established a method for acquisition, by means of a multi-distance probe (one emitter and two detectors), of the brain-functional signals corresponding to left- and right-thumb/index finger tapping. Recently, the active locations of finger movements have been identified by using an invasive technique, that is, an electrocorticography array [47].

In this paper, we present a method for detection of the (absolute) concentrations of HbO and HbR using CW-fNIRS with a bundled-optode configuration (i.e., non-invasive optical methodology). Using the developed method, the HRs upon individual (right) finger movements will be characterized first and the activated brain regions in the (left) motor cortex will be identified. Finally, and perhaps significantly, it is shown that the active locations of the two finger movements are quite distinctive.

Theory: Bundled-Optode

According to the CW-fNIRS technique, the incident light from an emitter will travel into brain tissue and reflect back to the scalp of the subject. Therefore, a portion of the incident light can be measured by a detector placed a few centimeters (2~4 cm) apart from the emitter. This phenomenon is described by the Beer-Lambert equation as follows. (1) where Iin(t, λ) is the intensity of the incident light of wavelength λ at time t, Iout(t, λ) is the intensity of the detected light at time t (of the same wavelength λ), l indicates the distance between an emitter and a detector, and μ(t, λ) expresses the linear attenuation coefficient. Note that μ(t, λ) can be split into two parts [48]: (2) where μa(t, λ) and μs(λ), respectively, stand for the absorption and scattering coefficients. In the present work, the scattering coefficient μs(λ) is assumed constant throughout the brain region, whose value in this paper is used as follows, see [20]. (3) where γ1 = 1.45 mm-1 and γ2 = 4.5×10−4 nm-1. Using (1), the absorbance of the light at the i-th channel, Ai(t, λ), is computed as follows. (4) where the superscript i denotes the channel number, and are the intensities of the incident and detected lights at the i-th channel, respectively, and li indicates the distance between the emitter and the detector at the i-th channel.

Fig 1 illustrates the concept of the bundled-optode method for measurement of absolute HbO/HbR concentrations. Let the italic integer numbers (i.e., 1, 2, 3, 4, 5, 6, 7 and 8) denote the locations of emitters (square) and detectors (circles), and let the plain integer numbers represent the channels formed by the emitter/detector pairs. A short-separation channel is formed from an emitter and a shortly separated detector (i.e., 0.5 cm) whereas long-separation channels are established from an emitter and detectors separated apart (i.e., 1~4 cm). For instance, consider the configuration of one emitter and seven detectors shown in Fig 1(B). It is noted that, from optode 1, optode 2 is separated shortly and optodes 3, 4, 5, 6, 7, and 8 are separated apart. Hence, one short-separation channel (i.e., emitter 1 and detector 2 or, reversely, emitter 8 and detector 7) and six long-separation channels are made (i.e., emitter 1 and six detectors 3~8). In the present work, a 32-optode-configuration operating at a sampling rate of 1.81 Hz was used (see Fig 2). In this case, the time difference between two neighboring channels is 17 msec (i.e., 1/(1.81 × 32)). Therefore, it can be assumed that the linear attenuation coefficient in (4), μ(t, λ), is almost the same in the neighboring channels (for instance, Chs. 1~7).

thumbnail
Fig 1. Bundled-optode arrangement.

(a) Concept of bundled-optode arrangement; (b) 1-emitter and 7-detector configuration for HbO/HbR measurement (see Fig 2 for their assignments): Optode 1 (square) is emitter and optode 2 and optodes 3~8 are short- and long-separation detectors (circles), respectively.

https://doi.org/10.1371/journal.pone.0165146.g001

thumbnail
Fig 2. Optode arrangement on the left motor cortex for detection of right-finger movements.

Squares and circles are emitters and detectors, respectively, and C3 denotes the reference point from the International 10–20 System.

https://doi.org/10.1371/journal.pone.0165146.g002

Then, for two neighboring channels in a local region, the absorbance difference can be defined as (5) where i = 1, j = 2~7. Now, let (6) be the spatial distance between two neighboring channels i and j. Then, from (2) and (5)-(6), the linear attenuation coefficient can be approximated by the spatial absorbance gradient between two neighboring channels (i.e., δAij(t,λ)/δlij) as follows. (7) where aHbO(λ) and aHbR(λ) are the extinction coefficients of HbO and HbR at wavelength λ, respectively, and CHbO(t) and CHbR(t) represent the concentrations of HbO and HbR, respectively. The configuration of one emitter and seven detectors in Fig 1(B), therefore, will make six absorbance gradients.

Finally, for two wavelengths, the following two absorbance spatial gradient equations are obtained.

(8)(9)

Solving for CHbO(t) and CHbR(t), the (absolute) concentration values of HbO and HbR from a pair of i and j channels (i.e., one short-separation channel and one long-separation channel) in a local region can be obtained as follows: (10)

Method

Ethics Statement

This work was approved by the Institutional Review Board of Pusan National University. Written consents were obtained from all subjects prior to the start. All experiments were conducted in accordance with the ethical standards encoded in the latest Declaration of Helsinki.

Equipment

A dual-wavelength (760 and 830 nm) CW-fNIRS instrument (DYNOT; NIRx Medical Technologies, USA) was utilized to acquire the brain signals generated by the movements of two fingers (little, thumb) of the right hand. Forty-eight (48) channels in the left motor cortex at a sampling rate of 1.81 Hz were configured (see Fig 2).

Subjects and Experimental Design

Five healthy male subjects (mean age 36.2; range: 33~37 years) with shaved heads were invited to perform experiments (little and thumb finger movements) involving the left motor cortex. All of the subjects are normal. Four of them are right handed. None of them had any neurological impairment or mental disorder. Prior to the experiments, the nature of the experiment was clearly explained to the subjects. To eliminate any interference from noise, the experiments were conducted in a dark and quiet room. The subjects were asked to sit freely on a comfortable chair and not to move their body during the experiment. Before starting, the subjects were carefully trained how to move their fingers.

Fig 3 describes the experimental paradigm: One experiment comprises ten trials of little and thumb finger movements. One trial consists of a 10 sec task and a 20 sec rest. After a 20 sec initial rest before the first trial, each subject was asked to perform finger movements (i.e., flexion and extension) for 10 sec, which was randomly repeated five times for each finger. Therefore, one experiment took 320 sec. To increase the brain activity during the 10 sec task period, the subjects were asked to move their (right) little/thumb finger as fast as they could, without paying attention to the number of flexions/extensions, and to relax during the 20 sec rest period. A laptop computer with a fifteen-inch screen was utilized to display pictures indicating each finger. The distance between the subject’s eyes and the laptop screen was adjusted to approximately 60 cm, so that the subjects could clearly see the indicated finger. The subjects were also instructed to keep their eyes open during the experiments. During the rest period, a black screen was shown to relax the subjects’ eyes.

thumbnail
Fig 3. The experimental paradigm for ten trials of right finger movements (i.e., one trial is 10 sec flexion/extension of a finger followed by 20 sec rest).

The five trials of little/thumb finger movements are randomly performed after the initial 20 sec rest.

https://doi.org/10.1371/journal.pone.0165146.g003

Data Processing

The measured fNIRS data are the intensities of light, which are to be converted to the absolute concentrations of HbX (i.e., HbO and HbR) using (10) and our own Matlab (the Mathworks, Inc., USA) codes. First, a baseline-correction method is utilized to eliminate the drifting phenomena in the acquired fNIRS data [5,49]. In our work, a fourth order polynomial has been fit to the recorded HbX data of each channel, and the fitted curve has been subtracted from the measured HbX signal [50]. Second, the baseline-corrected HbX signals are filtered by a low-pass filter with a cut-off frequency of 0.15 Hz to remove any respiration and pulse-related physiological noises.

Bundled-Optode Configuration

In the present work, a bundled-optode configuration is used to improve the spatial resolution, which can detect the brain activities in local brain regions. The configuration in Fig 2 consists of 8 emitters and 24 detectors resulting in 48 channels. As in Fig 1(B), a local brain region has been specified by six channels. In each local region, the six absolute HbX concentrations are measured using the developed method in the “Theory” section. The proposed method uses pairs of short- and long-separation channels to compute the absolute values of hemoglobin concentrations (see (10)), while the conventional MBLL computes the variations of HbX concentrations of a single channel.

Contrast-to-Noise Ratio

Contrast-to-noise ratio (CNR) computes the difference of the means of the signals during the task and the rest [5153]. To validate our proposed method, the CNR is used to investigate the signal to noise ratio from the sense that a high CNR value indicates a high ratio of the signal upon the task against that to noise. The CNR is computed as follows. (11) where “task” indicates the task period and “rest” denotes the rest period before starting to move fingers. In the present work, the duration of 6~16 sec was assigned to “task” and that of -6~0 sec from the onset time was set to “rest”.

Regions of Interest

The regions of interest (ROIs) are the brain regions consisting of the channels in which the t-values are higher than the critical t-value (tcrt) [54]. In this work, the tcrt is set to 1.674 according to the degree of freedom (M– 1 = 53, M is number of data samples of a trial) and the significance level (i.e., α = 0.05 for one-tailed test). The t-value of the averaged HbO of each channel is computed using the “robustfit.m” function available in Matlab. With this function, the averaged HbO in each channel is compared to the expected HR model (or the designed HR) by the least squares method, in which the designed HR is produced by convolving a trial period (i.e., 10 sec task and 20 sec rest) with the canonical HR function (i.e., the difference of two gamma functions) [5,10,31]. Then, the t- and p-values are used to determine the active channels: The criteria for “active” are t-value > 0 and p-value < 0.05. If these criteria are not met, the channel is considered inactive.

Local Active Brain Regions

Our objective is to evaluate the brain activities in local brain regions. Fig 4 presents the brain regions related to fingers (blue dashed rectangle) established by electrocorticogram [41]. In our experiment, the movements of little and thumb fingers are randomly performed to enhance the activity strength in the corresponding brain regions [37]. For each experiment, 48 channels with 10 trials per channel (i.e., 48 × 10 = 480 total trials) are considered. The coordinate of each channel is calculated from the known probe geometry in Fig 2. To present the active locations in the three-dimensional (3D) space, the following steps are used: Step i) the HbOs in each channel are averaged across five trials of each finger; Step ii) if the t-value of the averaged HbO is greater than tcrt = 1.674 and the p-value is less than 0.05, the channel is identified as a channel in the ROI. Next, the mean value of the 4~14 sec window of the identified channel is computed. On the other hand, if the t-value of the channel is less than tcrt (i.e., the channel does not belong to the ROI), whose mean value is set to zero to make a clear 3D image, see Fig 5; Step iii) the mean values of 48 channels are normalized within the 0~1 range; Step iv) the normalized mean-value of the averaged HbO of each channel in the ROI is coded in color and plotted at its 3D coordinate using the “fscatter3.m” function available at https://www.mathworks.com/matlabcentral/fileexchange/2993-fscatter3-m. The 3D coordinates assignment procedure of the channels is described below. As seen in Fig 5, the highly active locations are identified with red colors (i.e., large mean values).

thumbnail
Fig 4. Brain regions of fingers directly detected by electrocorticogram in the motor cortex (blue dashed rectangle), see Fig 9 in [41].

https://doi.org/10.1371/journal.pone.0165146.g004

thumbnail
Fig 5. Local activation maps of little (left) and thumb (right) finger movements.

The small circles indicate the optodes’ positions and the color bar at the bottom denotes the activation strength according to the normalized mean value between 0~1 range.

https://doi.org/10.1371/journal.pone.0165146.g005

Three-Dimensional Brain Map

In this work, the probe geometry is known (see Fig 2). Optode 1 is set to the reference coordinate. The coordinates of other optodes are computed based on optode 1. The position of each channel in the 3D space (x, y, z) is calculated from the optodes’ positions. Let the positions of optode 1 (emitter) and optode 17 (detector) be (x1, y1, z1) and (x2, y2, z2), respectively. Then, the coordinate in the 3D space of the target point of this emitter-detector pair is computed as (x, y, z) = ((x1 + x2)/2, (y1 + y2)/2, sqrt((x1x2)2 + (y1y2)2)/2) [55,56]. In this case, the z coordinate indicates the depth in which the brain activity has occurred. According to the optode arrangement in Fig 2, the positions of 48 channels in the 3D space have been computed. The mean HbO value during the 4~14 sec window of each channel is plotted at its computed position with a coded color. Since the depths of channels are different, the active locations at different tissue layers can be detected [57]. Also, our proposed method can effectively remove the superficial noise from the scalp’s layer using (5). Since the banana shaped paths are closely located, it can give a local brain map with much improved spatial resolution [6,58].

Results

We developed a bundled-optode method to measure the concentrations of the HbX of little and thumb finger movements. In most cases, HbOs increased and HbRs decreased during the activation tasks. However, in some channels, HbRs also increased during the tasks. Since some early works [31,59,60] revealed that HbO is more sensitive and more reliable than HbR, HbO data have been mostly analyzed. In addition, some channels showed large baselines drifts [5,49]. In the present paper, a fourth order polynomial (i.e., f(x) = a1x4 + a2x3+ a3x2 + a4x+ a5) was fit to the recorded HbX in each channel [50]. To correct the drifting, the fitted curve has been subtracted from the measured HbX of the corresponding channel. Fig 6(A) shows the HbO (blue thin curve) of a representative channel (Ch. 14, Subject 5) before using the baseline-correction method and the fitted signal (red thick curve), f(x). Fig 6(B) depicts the corrected HbO signal (blue thick curve). As can be seen, the proposed baseline-correction method can effectively remove the baseline drifting in the fNIRS data.

thumbnail
Fig 6. Illustration of the baseline-correction method (Ch. 14, Subject 5).

(a) An HbO signal before applying the baseline-correction method (blue thin curve) and the fitted signal f(x) (red thick curve); (b) the corrected HbO signal (blue thick curve).

https://doi.org/10.1371/journal.pone.0165146.g006

Fig 7 compares the HbOs computed by i) the conventional MBLL method (see Fig 7(A), 7(C) and 7(E)) and ii) the bundled-optode approach (see Fig 7(B), 7(D) and 7(F)) of an active channel (Ch. 20, Subject 5). The blue curves in Fig 7(A) and 7(B), respectively, are the HbO raw data obtained by the MBLL method and by Eq (10) (after the baseline correction was performed). The red thick curves in Fig 7(A) and 7(B), respectively, are the low-pass filtered HbOs from the blue curves. Gray bars in Fig 7(A) and 7(B) illustrate ten stimuli/trials of finger movements. Five trials (i.e., 1, 4, 6, 7, and 9) were assigned to little finger movements and other five trials (i.e., 2, 3, 5, 8, and 10) to thumb finger movements. Fig 7(C) and 7(E) plot the averaged HbOs of little and thumb finger movements, respectively, from Fig 7(A). Fig 7(D) and 7(F) show the averaged HbOs of little and thumb finger movements, respectively, from Fig 7(B). As can be seen, most HbOs increase during the 10 sec task period, and return to the baseline after the stimulation period is over. However, in a few trials (i.e., trials 1 and 7 in Fig 7(A)), it is seen that the HbO can begin to decrease prior to the end of the task: This can happen when the subject loses focus on the experiment.

thumbnail
Fig 7. Comparison of the HbOs computed by the conventional MBLL method and the bundled-optode approach (Ch. 20, Subject 5, blue thin curves: raw data, red thick curves: low-passed filtered).

(a) HbOs computed by the MBLL approach; (b) HbOs computed by the proposed method; (c)-(d) the averaged HbOs of the little finger movements (averaged over trials 1, 4, 6, 7, and 9) in (a)-(b), respectively; (e)-(f) the averaged HbOs of the thumb finger movements (averaged over trials 2, 3, 5, 8, and 10) in (a)-(b), respectively; the boundaries in (c)-(f) denote one standard deviation.

https://doi.org/10.1371/journal.pone.0165146.g007

In our proposed method, a short-separation channel as well as long-separation channels is utilized to compute the absolute concentrations of HbOs. To evaluate the proposed method, the CNRs of HbOs were computed. Fig 8 plots the CNRs of HbOs of 48 channels of a representative subject (Subject 5) for little (Fig 8(A)) and thumb (Fig 8(B)) movements. Red bars denote the CNRs of the HbOs computed by the conventional MBLL method, whereas blue bars indicate the CNRs of those calculated by the proposed approach. As can be seen, the proposed method gives much higher CNRs than those using the conventional MBLL, which demonstrates that the proposed method can effectively remove noises from the scalp. It is interesting to observe that the CNRs of the proposed method in some channels are opposite to those obtained by the MBLL approach (see channels 1, 23, 30, and 47 in both Fig 8(A) and 8(B)), which typically happen in inactive channels. This can be explained by the fact that the noise, when using the proposed method, has been reduced (in both active and inactive cases).

thumbnail
Fig 8. Comparison of the CNRs of HbOs in 48 channels (Subject 5).

(a) The CNR of little movements; (b) the CNR of thumb movements: Red bars are the CNRs of HbOs computed by the conventional MBLL method whereas blue bars indicate the CNRs of those calculated by the proposed approach.

https://doi.org/10.1371/journal.pone.0165146.g008

To illustrate one advantage (noise reduction) of the proposed method, an inactive channel (Ch. 30) is considered. Fig 9(A) and 9(B) present the HbOs computed by the MBLL method and the proposed method, respectively. Like Fig 7(C) and 7(E), Fig 9(C) and 9(E) illustrate the averaged HbOs over 5 trials of little finger movements (i.e., 1, 4, 6, 7, and 9) and other 5 trials of thumb movements (i.e., 2, 3, 5, 8, and 10), respectively, from Fig 9(A). Fig 9(D) and 9(F) depict the averaged HbOs obtained by the proposed method. As seen in Fig 9(C) and 9(E), the HbOs obtained by the MBLL method are more fluctuating than those obtained by the proposed method (particularly in the resting state), which are due to noises. Recall that the mean HbO during the task is calculated between 6~16 sec and that during the rest is computed between -6 to the onset time. Typically, the mean during the task in Fig 9(C) is 2.20 × 10−4 μM and that during the rest is 4.71 × 10−4 μM (see the blue curve). Therefore, the CNR of the little finger movement task in Fig 9(C) is -0.41. Similarly, the mean during the task in Fig 9(E) is -4.47 × 10−5 μM while that during the rest is -1.19 × 10−5 μM (see the blue curve), therefore the CNR of the thumb movement task in Fig 9(E) becomes -0.05. However, if using the proposed method, the mean value during the task in Fig 9(D) is 3.63 μM and that during the rest is 3.59 μM. Therefore, the CNR becomes 1.32 for the little finger movement task in Fig 9(D). Similarly, the mean during the task in Fig 9(F) is 3.64 μM while that during the rest is 3.58 μM. Therefore, the CNR becomes 2.33. This explains why some CNR bar charts in Fig 8 are in opposite directions, and furthermore illustrates the noise removal capability of the proposed method.

thumbnail
Fig 9. Comparison of two approaches for an inactive channel (Ch. 30, Subject 5).

(a) HbOs computed by the MBLL approach (thin blue: raw, thick red: low-pass filtered); (b) HbOs by the proposed method; (c)-(d) the averaged HbO across 5 trials in (a) and (b), respectively, for little finger; (e)-(f) the averaged HbO for thumb finger.

https://doi.org/10.1371/journal.pone.0165146.g009

Additionally, to illustrate the capability of the proposed method, Fig 10 is provided (Ch. 33, Subject 5): Fig 10(A) depicts the HbOs produced by the MBLL method, whereas Fig 10(B) by the proposed method. As seen in Fig 10(A), the fluctuation of the HbO is not consistent with the given stimuli, whereas such fluctuation is very consistent in Fig 10(B). In this case, Ch. 33 would be judged inactive if the conventional MBLL method were used, whereas it would be judged active if the proposed method were used.

thumbnail
Fig 10. An example showing that the MBLL method cannot tell but the proposed method can tell (Ch. 33, Subject 5).

(a) HbOs computed by the conventional MBLL method (blue thin curve: raw data, red thick curve: low-passed filtered); (b) HbOs obtained by the proposed method: In this case, the channel will be judged inactive by the MBLL method, while it will be judged active by the proposed method.

https://doi.org/10.1371/journal.pone.0165146.g010

To locate the active territories upon right finger movements, the ROIs were first identified. That is, if the t-value of a channel was greater than tcrt, the mean value of the averaged HbO of that channel was calculated. In contrast, if the t-value was not greater than tcrt, its mean value was set to zero. Then, colors for the mean values were coded to plot a 3D image. The steps for locating the active brain regions were described in detail in the “Local active brain regions” section. Fig 5 presents the active locations of little and thumb finger movements of five subjects. Maps in the left side represent little finger movements, whereas those in the right side represent thumb movements. The small circles denote the optodes’ positions. The color bar at the bottom indicates the activation strength of the normalized mean value between 0~1 range. As can be seen, the active territories upon little and thumb finger movements are clearly distinguishable for Subjects 1, 4, and 5. However, for Subjects 2 and 3, the active locations of two finger movements are almost indistinguishable.

To assess the representative HR for each finger movement, the HbXs were averaged across all active channels within the ROIs (i.e., t-value > tcrt). Fig 11 describes the averaged HbOs and HbRs of little and thumb finger movements (Subject 4). Fig 11(A) and 11(B) plot the averaged HbO/HbR of little finger movements, respectively. Fig 11(C) and 11(D) illustrate the averaged HbO/HbR of thumb movements, respectively. Fig 12 presents the averaged HRs of little and thumb finger movements averaged across five subjects. Fig 12(A) and 12(B) show the averaged HbO/HbR of little finger movements, respectively. Fig 12(C) and 12(D) picture the averaged HbO/HbR of thumb movements, respectively. As can be seen, at the onset time, the starting HbO values of little and thumb finger movements (see Fig 12(A)–12(C)) are 5.27 μM and 4.97 μM, respectively whereas those for HbRs are 3.44 μM and 3.35 μM (see Fig 12(B)–12(D)). It is noteworthy that the mean, peak, time-to-peak values of little finger movements are higher than those of thumb movements (see Table 1). However, the skewness of little finger movements is smaller than that of thumb movements. Assuming that the HRs of little and thumb finger movements are the same. The two-sample t-test (i.e., the “ttest2.m” function is available in Matlab) was used to test our hypothesis. The experimental results demonstrated that the HR of little and thumb finger movements are different (p-value < 0.05).

thumbnail
Fig 11. Averaged HbOs and HbRs of little and thumb finger movements (Subject 4).

(a)/(b) present the averaged HbO/HbR of little finger movements, (c)/(d) depict the averaged HbO/HbR of thumb movements, (e)/(f) plot the HbO and HbR of little/thumb finger movements (standardized in the 0~1 range) from (a)-(b)/(c)-(d), and the boundaries denote one standard deviation.

https://doi.org/10.1371/journal.pone.0165146.g011

thumbnail
Fig 12. HbOs and HbRs of little and thumb finger movements averaged across five subjects (see S1 Dataset for more details).

(a)/(b) present the averaged HbO/HbR of little finger movements, (c)/(d) plot the averaged HbO/HbR of thumb movements, (e)/(f) plot the HbO and HbR of little/thumb finger movements (standardized in the 0~1 range) from (a)-(b)/(c)-(d), and the boundaries denote one standard deviation.

https://doi.org/10.1371/journal.pone.0165146.g012

thumbnail
Table 1. Comparison of parameters of HRs of little and thumb finger movements.

https://doi.org/10.1371/journal.pone.0165146.t001

Discussion

In this work, a bundled-optode method based on the SRS principle was developed to measure the absolute concentrations of the HbXs upon two finger movements. It is noteworthy that multiple pairs of short- and long-separation channels were utilized in our proposed method for superficial noise removal [28,29]. In the fNIRS area, there is no gold standard to define a short- and/or long-separation detectors yet. However, several researchers demonstrated that the short-separation detector with distance 0.5~1.0 cm could measure the systemic hemodynamic fluctuations from the superficial layers [6,28,29]. In our current work, the short-separation detectors were placed at 0.5 cm to measure the superficial noise. The relative distance for long-separation detectors is normally arranged ~ 3 cm [29]. For our bundled-optode configuration, the long-separation detectors arranged with different distances (1.5, 2.5, 3.0, 3.5, and 4.0 cm) could detect both superficial noises and the brain activities in the different depths. By using (5) in our proposed method, the superficial noises were reduced. Additionally, the patterns in local active brain regions could be precisely detected. Our experimental results demonstrated that the proposed approach can effectively reduce the superficial noise (see Figs 810). The CNRs of HbOs computed by using our proposed method were greater than the CNRs of those using the conventional MBLL approach. Furthermore, our proposed method was applied to detect the HRs of little and thumb finger movements and identified their highly active territories as well. In some recent works, the active locations of individual finger movements were determined by using the center of mass, according to which voxel clusters greater than 4 mm2 were extracted for each finger [37,42]. In our work, the ROIs were identified and plotted in the 3D space in order to present the highly active territories upon individual finger movements. According to our experimental results, the highly active locations upon the movements of the little and thumb fingers were distinguishable. The active locations upon the thumb movements are in the lateral aspect, while those upon the little finger movements are in the medial aspect (see Fig 5, Subjects 1, 4, and 5). Additionally, the highly active locations of little finger movements are deeper than those of thumb movements (see Fig 5, Subjects 1 and 4). Our obtained results are quite consistent with the relevant previous studies [35,37,42,44]. Even though the obtained results are preliminary, they will certainly contribute to the development of future brain-computer interface applications using non-invasive technique [6169]. For example, the movements of the different fingers could be used to formulate commands for control of external devices such as robot arms or wheelchairs [70,71]. Moreover, the method could be extended to 3D fNIRS imaging based on each local brain region.

Reiterating the experimental results, the HRs of little and thumb finger movements are different. Table 1 depicts the parameters of the HRs of two fingers. The mean, peak, and time-to-peak values of the HRs of little finger movements are 5.34 μM, 5.35 μM, and 9.94 sec, respectively, which are quite higher than those of the HRs of thumb movements (i.e., 5.04 μM, 5.06 μM, and 9.39 sec, respectively). However, their skewness values are no significant difference (0.53 ± 0.06).

For Subjects 2 and 3, the active territories of the individual finger movements were indistinguishable (see Fig 5). This phenomenon, which had emerged in some of the relevant previous work [3640], can be explained in two ways: i) the muscle motions of each finger are affected by the others in the same hand; ii) the spatial resolution of our configuration is not yet sufficiently perfect to detect these overlapping regions.

Finally, the current work has certain limitations. Our bundled-optode arrangement comprises 32 optodes that can only record the brain activity in local brain regions with area 1.5 cm × 4 cm. Moreover, we only investigated the active territories of little and thumb finger movements. The active locations of the index, middle and ring finger movements in the same hand were not checked. In future work, we will design a bundled-optode configuration with a larger area and conduct experiments with five-finger movements for the same hand.

Conclusions

In this paper, the bundled-optode method for measuring the absolute concentrations of HRs and identifying the active territories upon two finger (i.e., little and thumb) movements was presented. By identifying the ROIs and plotting them in the 3D space, the active territories upon the two finger movements were distinguishable. Our results showed that the highly active locations of the little and thumb finger movements are different. Additionally, the mean, peak and time-to-peak values of the HRs in active brain regions of little finger movements were higher than those of the HRs of thumb movements, whereas their skewness values were not significantly different. For our future work, further validation is required for all fingers on the same hand and a large subject group.

Supporting Information

S1 Dataset. The data set of little and thumb finger movements (averaged over 5 subjects).

For further understanding our study, recorded fNIRS data were provided as supporting information.

https://doi.org/10.1371/journal.pone.0165146.s001

(XLSX)

Acknowledgments

This research was supported by the Mid-Career Researcher Program (grant no. NRF-2014R1A2A1A10049727) and the Convergence Technology Development Program for Bionic Arm (grant no. 2016M3C1B2912986) through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning, Republic of Korea.

Author Contributions

  1. Conceptualization: KSH.
  2. Data curation: HDN.
  3. Formal analysis: KSH HDN YIS.
  4. Funding acquisition: KSH.
  5. Investigation: KSH HDN.
  6. Methodology: KSH.
  7. Project administration: KSH.
  8. Resources: KSH.
  9. Software: HDN.
  10. Supervision: KSH.
  11. Validation: KSH HDN YIS.
  12. Visualization: HDN.
  13. Writing – original draft: HDN.
  14. Writing – review & editing: KSH.

References

  1. 1. Al-Rawi PG, Smielewski P, Kirkpatrick PJ. Evaluation of a near-infrared spectrometer (NIRO 300) for the detection of intracranial oxygenation changes in the adult head. Stroke. 2001; 32: 2492–2500. pmid:11692006
  2. 2. Delpy DT, Cope M, van der Zee P, Arridge S, Wray S, Wyatt JS. Estimation of optical pathlength through tissue from direct time of flight measurement. Phys Med Biol. 1988; 33: 1433. pmid:3237772
  3. 3. Khan MJ, Hong K-S. Passive BCI based on drowsiness detection: an fNIRS study. Biomed Opt Express. 2015; 6: 4063–4078. pmid:26504654
  4. 4. Naseer N, Hong K-S. Classification of functional near-infrared spectroscopy signals corresponding to the right-and left-wrist motor imagery for development of a brain–computer interface. Neurosci Lett. 2013; 553: 84–89. pmid:23973334
  5. 5. Santosa H, Hong MJ, Hong K-S. Lateralization of music processing with noises in the auditory cortex: an fNIRS study. Front Behav Neurosci. 2014; 8: 418. pmid:25538583
  6. 6. Scholkmann F, Metz AJ, Wolf M. Measuring tissue hemodynamics and oxygenation by continuous-wave functional near-infrared spectroscopy—how robust are the different calculation methods against movement artifacts? Physiol Meas. 2014; 35: 717. pmid:24622337
  7. 7. Teng Y, Ding H, Huang L, Li Y, Shan Q, Ye D, et al. Non-invasive measurement and validation of tissue oxygen saturation covered with overlying tissues. Prog Nat Sci. 2008; 18: 1083–1088.
  8. 8. Bhutta MR, Hong K-S, Kim B-M, Hong MJ, Kim Y-H, Lee S-H. Note: three wavelengths near-infrared spectroscopy system for compensating the light absorbance by water. Rev Sci Instrum. 2014; 85: 026111. pmid:24593411
  9. 9. Hu X-S, Hong K-S, Ge SS. fNIRS-based online deception decoding. J Neural Eng. 2012; 9: 026012. pmid:22337819
  10. 10. Santosa H, Hong MJ, Kim S-P, Hong K-S. Noise reduction in functional near-infrared spectroscopy signals by independent component analysis. Rev Sci Instrum. 2013; 84: 073106. pmid:23902043
  11. 11. Kamran MA, Hong K-S. Linear parameter-varying model and adaptive filtering technique for detecting neuronal activities: an fNIRS study. J Neural Eng. 2013; 10: 056002. pmid:23893789
  12. 12. Kamran MA, Hong K-S. Reduction of physiological effects in fNIRS waveforms for efficient brain-state decoding. Neurosci Lett. 2014; 580: 130–136. pmid:25111978
  13. 13. Haeussinger FB, Heinzel S, Hahn T, Schecklmann M, Ehlis AC, Fallgatter AJ. Simulation of near-infrared light absorption considering individual head and prefrontal cortex anatomy: Implications for optical neuroimaging. Plos One. 2011; 6: e26377. pmid:22039475
  14. 14. Tam ND, Zouridakis G. Temporal decoupling of oxy-and deoxy-hemoglobin hemodynamic responses detected by functional near-infrared spectroscopy (fNIRS). J Biomed Eng Med Imaging. 2014; 1: 18–28.
  15. 15. Everdell N, Airantzis D, Kolvya C, Suzuki T, Elwell C. A portable wireless near-infrared spatially resolved spectroscopy system for use on brain and muscle. Med Eng Phys. 2013; 35: 1692–1697. pmid:23706504
  16. 16. Barker JW, Panigrahy A, Huppert TJ. Accuracy of oxygen saturation and total hemoglobin estimates in the neonatal brain using the semi-infinite slab model for FD-NIRS data analysis. Biomed Opt Express. 2014; 5: 4300–4312. pmid:25574439
  17. 17. Dehaes M, Grant PE, Sliva DD, Roche-Labarbe N, Pienaar R, Boas DA, et al. Assessment of the frequency-domain multi-distance method to evaluate the brain optical properties: Monte Carlo simulations from neonate to adult. Biomed Opt Express. 2011; 2: 552–567. pmid:21412461
  18. 18. Hallacoglu B, Sassaroli A, Wysocki M, Guerrero-Berroa E, Beeri MS, Haroutunian V, et al. Absolute measurement of cerebral optical coefficients, hemoglobin concentration and oxygen saturation in old and young adults with near-infrared spectroscopy. J Biomed Opt. 2012; 17: 0814061–0814068.
  19. 19. Ward LM, Aitchison RT, Tawse M, Simmers AJ, Shahani U. Reduced haemodynamic response in the ageing visual cortex measured by absolute fNIRS. Plos One. 2015; 10: e0125012. pmid:25909849
  20. 20. Matcher S, Cope M, Delpy D. In vivo measurements of the wavelength dependence of tissue-scattering coefficients between 760 and 900 nm measured with time-resolved spectroscopy. Appl Optics. 1997; 36: 386–396.
  21. 21. Ntziachristos V, Ma X, Yodh A, Chance B. Multichannel photon counting instrument for spatially resolved near infrared spectroscopy. Rev Sci Instrum. 1999; 70: 193–201.
  22. 22. Kek KJ, Kibe R, Niwayama M, Kudo N, Yamamoto K. Optical imaging instrument for muscle oxygenation based on spatially resolved spectroscopy. Opt Express. 2008; 16: 18173–18187. pmid:18958095
  23. 23. Suzuki S, Takasaki S, Ozaki T, Kobayashi Y. Tissue oxygenation monitor using NIR spatially resolved spectroscopy. Proc. SPIE. 1999; 3597: 582–592.
  24. 24. Suzuki H, Niwayama M, Yamakawa T, Ohkubo M, Kime R, Katsumura T. Simultaneous determination of absorption coefficients for skin and muscle tissues using spatially resolved measurements. Trans Tech Publ. 2011; 309–312.
  25. 25. Valipour A, McGown A, Makker H, O'Sullivan C, Spiro S. Some factors affecting cerebral tissue saturation during obstructive sleep apnoea. Eur Resp J. 2002; 20: 444–450.
  26. 26. Farrell TJ, Patterson MS, Wilson B. A diffusion theory model of spatially resolved, steady-state diffuse reflectance for the noninvasive determination of tissue optical properties invivo. Med Phys. 1992; 19: 879–888. pmid:1518476
  27. 27. Leung TS, Tachtsidis I, Smith M, Delpy DT, Elwell CE. Measurement of the absolute optical properties and cerebral blood volume of the adult human head with hybrid differential and spatially resolved spectroscopy. Phys Med Biol. 2006; 51: 703. pmid:16424590
  28. 28. Saager RB, Telleri NL, Berger AJ. Two-detector Corrected Near Infrared Spectroscopy (C-NIRS) detects hemodynamic activation responses more robustly than single-detector NIRS. Neuroimage. 2011; 55: 1679–1685. pmid:21256223
  29. 29. Gagnon L, Cooper RJ, Yücel MA, Perdue KL, Greve DN, Boas DA. Short separation channel location impacts the performance of short channel regression in NIRS. Neuroimage. 2012; 59: 2518–2528. pmid:21945793
  30. 30. Hu X-S, Hong K-S, Ge SS. Recognition of stimulus-evoked neuronal optical response by identifying chaos levels of near-infrared spectroscopy time series. Neurosci Lett. 2011; 504: 115–120. pmid:21945547
  31. 31. Hong K-S, Nguyen H-D. State-space models of impulse hemodynamic responses over motor, somatosensory, and visual cortices. Biomed Opt Express. 2014; 5: 1778–1798. pmid:24940540
  32. 32. Jasdzewski G, Strangman G, Wagner J, Kwong K, Poldrack R, Boas DA. Differences in the hemodynamic response to event-related motor and visual paradigms as measured by near-infrared spectroscopy. Neuroimage. 2003; 20: 479–488. pmid:14527608
  33. 33. Aqil M, Hong K-S, Jeong M-Y, Ge SS. Detection of event-related hemodynamic response to neuroactivation by dynamic modeling of brain activity. Neuroimage. 2012; 63: 553–568. pmid:22796989
  34. 34. Dechent P, Frahm J. Functional somatotopy of finger representations in human primary motor cortex. Hum Brain Mapp. 2003; 18: 272–283. pmid:12632465
  35. 35. Kleinschmidt A, Nitschke MF, Frahm J. Somatotopy in the human motor cortex hand area. A high-resolution functional MRI study. Eur J Neurosci. 1997; 9: 2178–2186. pmid:9421177
  36. 36. Meier JD, Aflalo TN, Kastner S, Graziano MS. Complex organization of human primary motor cortex: a high-resolution fMRI study. J Neurophysiol. 2008; 100: 1800–1812. pmid:18684903
  37. 37. Olman CA, Pickett KA, Schallmo MP, Kimberley TJ. Selective BOLD responses to individual finger movement measured with fMRI at 3T. Hum Brain Mapp. 2012; 33: 1594–1606. pmid:21674691
  38. 38. Schieber MH, Hibbard LS. How somatotopic is the motor cortex hand area? Science. 1993; 261: 489–492. pmid:8332915
  39. 39. Rao S, Bandettini P, Binder J, Bobholz J, Hammeke T, Stein EA, et al. Relationship between finger movement rate and functional magnetic resonance signal change in human primary motor cortex. J Cereb Blood Flow Metab. 1996; 16: 1250–1254. pmid:8898698
  40. 40. Yue GH, Liu JZ, Siemionow V, Ranganathan VK, Ng TC, Sahgal V. Brain activation during human finger extension and flexion movements. Brain Res. 2000; 856: 291–300. pmid:10677638
  41. 41. Penfield W, Rasmussen T. The cerebral cortex of man; a clinical study of localization of function. 1st ed. NewYork: The Macmillan company; 1950.
  42. 42. Beisteiner R, Windischberger C, Lanzenberger R, Edward V, Cunnington R, Erdler M, et al. Finger somatotopy in human motor cortex. Neuroimage. 2001; 13: 1016–1026. pmid:11352607
  43. 43. Diedrichsen J, Wiestler T, Krakauer JW. Two distinct ipsilateral cortical representations for individuated finger movements. Cerebral Cortex. 2013; 23: 1362–1377. pmid:22610393
  44. 44. Schieber MH. Constraints on somatotopic organization in the primary motor cortex. J Neurophysiol. 2001; 86: 2125–2143. pmid:11698506
  45. 45. Habermehl C, Holtze S, Steinbrink J, Koch SP, Obrig H, Mehnert J, et al. Somatosensory activation of two fingers can be discriminated with ultrahigh-density diffuse optical tomography. Neuroimage. 2012; 59: 3201–3211. pmid:22155031
  46. 46. Yamada T, Umeyama S, Matsuda K. Multidistance probe arrangement to eliminate artifacts in functional near-infrared spectroscopy. J Biomed Opt. 2009; 14: 064034. pmid:20059272
  47. 47. Hotson G, McMullen DP, Fifer MS, Johannes MS, Katyal KD, Para MP, et al. Individual finger control of a modular prosthetic limb using high-density electrocorticography in a human subject. J Neural Eng. 2016; 13: 026017. pmid:26863276
  48. 48. Tuchin VV. Light scattering study of tissues. Physics-Uspekhi. 1997; 40: 495–515.
  49. 49. Shah A, Seghouane A-K. An Integrated Framework for Joint HRF and Drift Estimation and HbO/HbR Signal Improvement in fNIRS Data. IEEE Trans Med Imaging. 2014; 33: 2086–2097. pmid:24956281
  50. 50. Laughner JI, Ng FS, Sulkin MS, Arthur RM, Efimov IR. Processing and analysis of cardiac optical mapping data obtained with potentiometric dyes. Am J Physiol-Heart Circul Physiol. 2012; 303: H753–H765.
  51. 51. Cui X, Bray S, Reiss AL. Functional near infrared spectroscopy (NIRS) signal improvement based on negative correlation between oxygenated and deoxygenated hemoglobin dynamics. Neuroimage. 2010; 49: 3039–3046. pmid:19945536
  52. 52. Cui X, Bray S, Reiss AL. Speeded near infrared spectroscopy (NIRS) response detection. PloS One. 2010; 5: e15474. pmid:21085607
  53. 53. Welvaert M, Rosseel Y. On the definition of signal-to-noise ratio and contrast-to-noise ratio for fMRI data. PloS One. 2013; 8: e77089. pmid:24223118
  54. 54. Hong KS, Santosa H. Decoding four different sound-categories in the auditory cortex using functional near-infrared spectroscopy. Hear Res. 2016; 333: 157–166. pmid:26828741
  55. 55. Taga G, Homae F, Watanabe H. Effects of source-detector distance of near infrared spectroscopy on the measurement of the cortical hemodynamic response in infants. Neuroimage. 2007; 38: 452–460. pmid:17884584
  56. 56. Nguyen HD, Hong KS. Bundled-optode implementation for 3D imaging in functional near-infrared spectroscopy. Biomed Opt Express. 2016; 7: 3491–3507. pmid:27699115
  57. 57. Strangman GE, Li Z, Zhang Q. Depth sensitivity and source-detector separations for near infrared spectroscopy based on the colin27 brain template. Plos One. 2013; 8: e66319. pmid:23936292
  58. 58. Boas D, Chen K, Grebert D, Franceschini M. Improving the diffuse optical imaging spatial resolution of the cerebral hemodynamic response to brain activation in humans. Opt Lett. 2004; 29: 1506–1508. pmid:15259728
  59. 59. Hu X-S, Hong K-S, Ge SS, Jeong M-Y. Kalman estimator-and general linear model-based on-line brain activation mapping by near-infrared spectroscopy. Biomed Eng Online. 2010; 9: 82. pmid:21138595
  60. 60. Hoshi Y. Functional near-infrared spectroscopy: current status and future prospects. J Biomed Opt. 2007; 12: 062106–062109. pmid:18163809
  61. 61. Naseer N, Hong MJ, Hong K-S. Online binary decision decoding using functional near-infrared spectroscopy for the development of brain–computer interface. Exp Brain Res. 2014; 232: 555–564. pmid:24258529
  62. 62. Khan MJ, Hong MJ, Hong K-S. Decoding of four movement directions using hybrid NIRS-EEG brain-computer interface. Front Hum Neurosci. 2014; 8: 244. pmid:24808844
  63. 63. Naseer N, Hong K-S. fNIRS-based brain-computer interfaces: a review. Front Hum Neurosci. 2015; 9: 3. pmid:25674060
  64. 64. Hong K-S, Naseer N, Kim Y-H. Classification of prefrontal and motor cortex signals for three-class fNIRS–BCI. Neurosci Lett. 2015; 587: 87–92. pmid:25529197
  65. 65. Naseer N, Hong K-S. Decoding answers to four-choice questions using functional near infrared spectroscopy. J Near Infrared Spectrosc. 2015; 23: 23–31.
  66. 66. Bhutta MR, Hong MJ, Kim Y-H, Hong K-S. Single-trial lie detection using a combined fNIRS-polygraph system. Front Psychol. 2015; 6: 709. pmid:26082733
  67. 67. Hong K-S, Naseer N. Reduction of delay in detecting initial dips from functional near-infrared spectroscopy signals using vector-based phase analysis. Int J Neural Syst. 2016; 26: 1650012. pmid:26971785
  68. 68. Tidoni E, Gergondet P, Kheddar A, Aglioti SM. Audio-visual feedback improves the BCI performance in the navigational control of a humanoid robot. Front Neurorobotics. 2014; 8: 20.
  69. 69. Urgen BA, Plank M, Ishiguro H, Poizner H, Saygin AP. EEG theta and Mu oscillations during perception of human and robot actions. Front Neurorobotics. 2013; 7: UNSP 19.
  70. 70. Dura-Bernal S, Zhou XL, Neymotin SA, Przekwas A, Francis JT, Lytton WW. Cortical spiking network interfaced with virtual musculoskeletal arm and robotic arm. Front Neurorobotics. 2015; 9: 13.
  71. 71. Kocaturk M, Gulcur HO, Canbeyli R. Toward building hybrid biological/in silico neural networks for motor neuroprosthetic control. Front Neurorobotics. 2015; 9: 8.