EROS and TR-fNIRS report different kinds of data in regards to brain activity. EROS detects changes in the light scattering properties of neural tissue, thought to be related to individual neuronal activity. The latter uses the same light signal and photon arrival time measurement technology (NIRS) to track oxygenated and deoxygenated hemoglobin changes. This report is a brief, cursory proposal of combining EROS and time-resolved fNIRS into a working hybrid system as a literature review. The table below compares the three relevant neuroimaging techniques on a wide array of dimensions, from spatial and temporal resolution to real-world accessibility that will be discussed in the rest of this report.
Combining event-related optical signal (EROS) and time-resolved functional infrared spectroscopy (TR-fNIRS)
Introduction
Comparison of Techniques
| Techniques | Measures | Temporal resolution | Signal origin | Relatively low-cost | Portability |
|---|---|---|---|---|---|
| EROS | Individual photon arrival times to determine changes in light scattering | 10-100 ms | Neuronal activation | Yes | Yes |
| fNIRS | Changes in oxygenated and deoxygenated hemoglobin concentrations | Scale of seconds | Neurovascular coupling | Yes | Yes |
| Time-resolved techniques | Individual photon arrival times to determine absolute hemoglobin concentrations | 10-100 ms | Neurovascular coupling | No | No |
Discussion
fNIRS has good spatial resolution, although not as powerful as fMRI, but poorer temporal resolution due to the detection of blood-flow responses being slower (on the scale of seconds). With its strengths in both temporal (on the scale of milliseconds) and spatial recording, EROS is prone to signal contamination and produces a lower signal-to-noise ratio due to measurements being limited to the cerebral cortex (Gratton & Fabiana, 2001). fNIRS detection also has limited depth, causing it to suffer from the same problem. Time-resolved NIRS enhances the depth sensitivity of NIRS by recording the time it took for single photons to arrive at the detector (light that travels farther into brain tissue takes longer, or has a longer “time of flight,” than light travelling closer to the surface) (Abdalmalak et al., 2020). This greater depth penetration reduces scalp and superficial tissue contamination to the signal.
Benefits and Challenges
The main benefit of combining these specific techniques is the combination of hemodynamic data with neuronal activation responses. Detection of light scattering in EROS, although limited to only structures up to around 3 cm of the brain, is still better than EEG and EROS also compensates for spatial resolution lacking by ERPs (Tse et al., 2007). Secondly, time-resolved techniques can provide absolute values (in hemoglobin content) as opposed to fNIRS alone. These technologies also make up for the poorer temporal resolution in fMRI.
The devices used in time-resolved techniques are more complex than fNIRS. Furthermore, in the integration of EROS with TR-fNIRS, the main obstacle is that light scattering itself limits depth sensitivity, as diffusion of photons in a medium decreases photons detected in one focused direction (Torricelli et al., 2013). This problem would have to be explored and solved to make this combined system an effective one. EROS and fNIRS can be used to measure responses concurrently, modulating light beams to gain time-resolved responses. However, building in the separate measurement systems for both EROS and TR-fNIRS requires sophistication, for instance in synchronizing temporal measurements. The more millisecond scale of neural activity detected by EROS is to be synchronized with the slower hemodynamic responses detected by (TR-)fNIRS. One solution is using a single clock to time-lock measurements to specific events. Thus, we synchronize data at the occurrence or onset of specific events, such as elicited responses to experimental stimulus.
Technical Considerations
The same optodes can be used to emit and detect light for this multimodal system to be placed on the scalp. Although, again, the measurement systems for detecting EROS and for TR-fNIRS may need more sophistication. To produce anatomical images containing spatial data, we can scan cerebral images using MRI, then, after simultaneously acquiring date of light scattering properties and time-resolved hemodynamic changes from subjects, optode locations can be projected onto the MRI reconstructions (Chen et al., 2017). There are many steps to this procedure that can be time-consuming and computationally costly.
The detection and measurement devices of this multimodal system are ones already employed in their independent uses. To obtain the separate time-resolved data vs. EROS data—oxygenated hemoglobin concentrations and changes in light scattering/absorption properties in neural tissue respectively—two data can be collected: 1) the intensity of the light reaching the detector and 2) the time-of-flight of individual photons reaching the detectors. Mapping the resulting EROS and hemoglobin concentration readings should be displayed in combined activation topographic maps of the target brain areas across multiple snapshots in time for comparison, especially because EROS data is differential. Photon delay or time-of-flight should give spatial information about what area of the brain the signal is coming from. Oxygenated and deoxygenated hemoglobin concentrations at differing layers of the brain provides spatial depth information. Meanwhile, calculations of EROS derived from the light reaching the detector all-in-all is a better tell of temporal information. Activation maps and time course graphs would likely be based on elicited or target events in experimental situations.
Applications and Future Directions
The figures above display example sketches of time course graphs for 1) EROS events based on photon delay and 2) oxygenated blood concentration levels time stamps. We can overlap oxygenated hemoglobin concentration and deoxygenated hemoglobin concentration. The timing is specific to a specific event, for example, a stimulus presentation or a simple motor task. By comparing these graphs, information that can be attained from EROS and TR-fNIRS, we can get temporal and spatial information about changes in activity in neural tissue and blood oxygen. In essence, data from this multimodal data will likely be richer when focused on responses to specific events.
It is especially because these individual systems are still emerging technologies that there are still challenges to be faced. Individually, EROS and fNIRS are relatively low cost and both are attractive options because of their portability. On the other hand, time-resolved NIRS technology is bulky, limiting dynamic range and portability, and is costly (Gratton & Fabiani, 2009). This would be an obstacle in widely introducing this multimodal system into labs and clinical settings. It also puts us far from more widespread, commercial uses. For instance, among detection systems in TR-NIRS such as non-linear optic gating, time-correlated single-photon counting (TCSPC), and the use of time-gated intensified charge coupled device (ICCD), the streak camera is one of significant consideration due to its high temporal resolution (1-10 picoseconds) and usability as a multi-channel detector (Torricelli et al., 2013) (more channels improve spatial resolution). Its greatest drawbacks are its costs and complexity, hindering it from being more widely used in clinical settings. Availability and costs for light sources and data acquisition for different techniques also pose a challenge, especially when the choices of devices to use narrow when we select for faster temporal detection (Contini et al., 2006). In terms of light sources, “there is always a trade-off between output power and pulse duration” (Torricelli et al., 2013). As an example, for pulsed diode lasers, to obtain a less than 100 picosecond pulse duration, an output (capacity) power is preferred at less than 1mW. This choice affects its clinical usability. All of these considerations put together explain the lack of commercial devices for TR-NIRS systems. It would be useful for TR-fNIRS to be further tested in clinical and BCI applications to strengthen our understanding of how it could work in combination with EROS.
Conclusions
These obstacles don’t diminish the usefulness of integrating these systems together. There are a number of existing studies exploring research questions and clinical applications where time-resolved techniques combined with EROS measurements would be useful. One study was mentioned above: whether it would be effective to detect hemodynamic shifts earlier using EEG combined with fNIRS (Jawad Khan et al., 2018). The study’s practical purpose in exploring this was in an attempt to reduce the time it took for mental commands to be detected for a quicker BCI response by having the detection devices predict them ahead of time. The strain that bulky time-resolved NIRS technology would put on clients, however, remains a drawback. Still, the idea of using different types of neuroimaging systems to create a smarter, more user-responsive single multimodal system is compelling in the world of BCIs. Similarly, some recent research (e.x., Abdalmalak et al., 2017; Abdalmalak et al., 2020) studied TR-fNIRS as a base for BCIs, particularly medical BCIs geared to patients with disorders of consciousness (DOC), because of the advantage of its depth sensitivity. Is a multimodal EROS and TR-fNIRS approach toward BCIs for patients with disorders of consciousness (DOC) not only a feasible one, but one that enhances current technologies? Can this technique generate BCI actions faster and more accurately from sharper data from brain regions, produced by statistically comparing neuronal changes and hemodynamic data? Another area of interest is measuring and studying cognitive load. fNIRS has been used in studies of cognitive load. For example, a recent study by Karmakar and colleagues (2023) generated topographical brain images in a fNIRS-based experiment exploring a deep learning model convolutional neural network (CNN) to classify different brain responses into different cognitive states. Stepping back to look at the bigger picture, we can ask how the proposed multimodal system in this report, with its greater depth sensitivity, can provide richer data to help us test and develop deep learning models in clinical and neuroscientific studies? Can this system extract clearer data that can further differentiate mental states—a problem addressed in the study discussed above—to improve classification accuracy? How can this help our understanding of cognitive load? As we move towards enhancing the portability and accessibility of TR-fNIRS and consider it being used in tandem with EROS detection, there is a potential for greater insight into these questions and real-world, clinical matters.
References
- Abdalmalak A., Milej D., Diop M., Shokouhi M., Naci L., Owen A.M., & St. Lawrence K. (2017). Can time-resolved NIRS provide the sensitivity to detect brain activity during motor imagery consistently? Biomedical Optics Express, 8 (4), 2162-2172. https://doi.org/10.1364/BOE.8.002162.
- Abdalmalak A., Milej D., Yip L.C.M., Khan A.R., Diop M., Owen A.M., & St. Lawrence K. (2020). Assessing time-resolved fNIRS for brain-computer interface applications of mental communication. Front. Neurosci., 14 (105). https://doi.org/10.3389/fnins.2020.00105.
- Chen M., Blumen H.M., Izzatoglu M., & Holtzer R. (2017). Spatial co-registration of functional near-infrared spectroscopy to brain MRI. J Neuroimaging, 27 (5), 453-460. https://doi.org/10.1111/jon.12432.
- Contini D., Torricelli A., Pifferi A., Spinelli L., Paglia F., & Cubeddu R. (2006). Multi-channel time-resolved system for functional near infrared spectroscopy. Optics Express, 14 (12), 5418-5432. https://doi.org/10.1364/OE.14.005418.
- Gratton G., & Fabiani M. (2001). Shedding light on brain function: The event-related optical signal. Trends in Cognitive Science, 5 (8), 357-363. https://doi.org/10.1016/S1364-6613(00)01701-0.
- Gratton G., & Fabiana M. (2009). Fast optical signals: Principles, methods, and Experimental Results. In Frostig R.D. (Ed.), In Vivo Optical Imaging of Brain Function. 2nd edition. Boca Raton (FL): CRC Press/Taylor & Francis. https://doi.org/10.1080/00107514.2010.487572.
- Jawad Khan J., Ghafoor U., & Hong K-S. (2018). Early detection of hemodynamic responses using EEG: A hybrid EEG-fNIRS study. Front. Hum. Neurosci., 12 (479). https://doi.org/10.3389/fnhum.2018.00479.
- Karkamar S., Kamilya S., Dey P., Guhathakurta P.K., Dalui M., Bera T.K., Halder S., Koley C., Pal T., & Basu A. (2022). Real time detection of cognitive load using fNIRS: A deep learning approach. Biomedical Signal Processing and Control, 80 (1). https://doi.org/10.1016/j.bspc.2022.104227.
- Torricelli A., Contini D., Pifferi A., Caffini M., Re R., Zucchelli L., & Spinelli L. (2013). Time domain functional NIRS imaging for human brain mapping. NeuroImage, 85 (1), 28-50. https://doi.org/10.1016/j.neuroimage.2013.05.106.
- Tse C-Y., Lee C-L., Sullivan J., Garnsey S.M., Dell G.S., Fabiana M., & Gratton G. (2007). Imaging cortical dynamics of language processing with the event-related optical signal. Proc. Natl. Acad. Sci. U.S.A., 104 (43), 17157-17162. https://doi.org/10.1073/pnas.0707901104.
Other Course Work
📊 Assignment 1: EEG Analysis
📓 assignment1 - Kyle Cordez.ipynb
Jupyter Notebook💡 Opens in Google Colab for interactive execution. Requires Google account.
🎨 Assignment 2: BrainImation
📄 Cordez_assignment 2 write-up.pdf
View Original PDFASSIGNMENT 2 WRITE-UP: “Attention-Controlled Rain”
Kyle Cordez
PSYCH 403
Using the “Attention Meter” as the base, I added rain as a visual that changes in speed as attention
shifts, specifically as the beta waves shifts.
The biggest addition is the Raindrop class that creates a single falling raindrop, with properties horizontal
position (x), vertical position (y), velocity (vel), height (h), and width (w), all of which are randomized
values within set ranges. A fixed number of raindrops are created in setup() and no new raindrops are
generated after. Raindrop has a method show() that draws it. I take advantage of the looping nature of
draw() to repeatedly set each drop to the top of the screen (randomized from 0 to negative height) when
it reaches the bottom of the screen (y=height), handled within the Raindrop method fall() that handles
the movement of each raindrop, by calling another method update(). update() also re-randomizes the
raindrop’s other properties. There is a loop that goes through every raindrop created, running show()
and fall() for each.
The most relevant property to this assignment is the velocity. The velocity is randomized in a range,
defined by global variables low_v and high_v, which have default values assigned in setup(). The function
updateVelocity() checks eegData.beta, then increases or decreases low_v and high_v. When
eegData.beta falls above 0.5, it is added to both low_v and high_v. When it falls below, it is subtracted
from both. At eegData.beta == 0.5, the values are set back to the default. The velocity range caps at 0.5
at the lowest and 15 at most. Since fall() is repeatedly called for each raindrop, the velocity is regularly
updated after the range is changed.
Other changes
To draw more attention to the rain (and simply for a more interesting visual experience), I made changes
to the formatting of the attention bar and the text labels, including:
– moving the labels and the bar away from the centre of the screen
– font sizes
– the text displaying the percentage of attention changes in hue the same as the attention bar
– change properties of the attention bar (size dimensions + no stroke)
🎥 Cordez_assignment 2 attention-controlled rain.mp4
💡 Videos require Google Drive access. Open in new tab if it doesn't load.
🎮 Cordez_assignment 2 attention-example.js
💡 Code is embedded in this portfolio - opens instantly in the live BrainImation editor (no internet required!)
🎯 Midterm Project
📄 midterm writeup - Kyle Cordez.pdf
View Original PDFMIDTERM PART 1: BCI INNOVATION
ADAPTIVE DIFFICULTY SYSTEM
The visuals are a simple mimicry of a vessel flying through space. As the ship appears to soar
through the stars (the movement to the bottom of the screen of the stars, defined by the Star
object, makes it appear so), the speed and the colour of the stars might change. This change is
dictated by the collected alpha and beta waves, which in turn are dictated by level of attention
or stress.
When stress is high (or attention levels are high), beta waves go up. Past a specific threshold,
eegData.beta is subtracted from the speed of the stars. The stars become desaturated, reverting
to all white. When the user is relaxed, alpha waves go up. At a threshold, the stars speed up by
increments of values of eegData.alpha. The stars become a range of more saturated colours.
The visuals of the stars become calm but remain interesting when neither of the above
conditions are met. The speed is somewhere in-between. The stars are a unified colour.
There is a scoring system that runs for a set amount of time; when the countdown hits 0, the
final score is recorded. One point is added to the score for every second the user remains in the
target values.
The visuals, controlled by levels of attention/stress, serve as direct visible feedback to the user
about their current brain state. Its practical use is providing the user information about their
levels of focus, and serving as a learning tool to understand how to adjust their mental state.
The scoring system and timer serves as further motivation. They track the user’s alertness levels
over a fixed period of time, allowing them to see how long they’ve been in a specific range of
brain waves or “on track.”
MIDTERM PART 2: ERP EXPERIMENT
AUDITORY N100
In the experiment, I chose to measure auditory N100. This component peaks after the onset of
an unpredictable auditory stimulus, elicited by the perception of it. Since it is stronger when the
stimulus is unpredictable, it gets weaker when a presented stimulus is repetitive and, clearly,
predicted. For unpredictability, in the design of the stimulus presentation, I have set three
different frequencies for each category manipulated for–high (>=2000Hz) vs. low(<=300)
frequency—to be randomly selected from and there is a 50/50 chance of it actually playing
every passed interval.
The system records only from two channels (TP9 and TP10, two lateral channels near heschel’s
gyrus) and displays them twice separately for each category. While EEG data is continuously
being added for averaging, it also adds specifically after a stimulus (trial) runs to the graphs of
the corresponding frequency category.
If I was testing with real EEG, I would expect to see peaks ~100ms after the stimulus has been
played. I would also expect to see, specifically, that the amplitude of the averaged EEG data
peaks greater when the most recently played stimulus has a frequency different from the
frequency of the stimulus played previously. This peak increases proportionally to the difference
between the frequency of the recent stimulus and the frequency of the last played stimulus.
This solidifies the findings that auditory N100 peaks are greater when the perceived stimulus is
unpredictable. As an example, we expect bigger amplitudes when the frequency of the most
recent stimulus is 2000Hz and the previous frequency was 250Hz, and smaller when the last
two frequencies were 200Hz and 250Hz as the difference is smaller.
🎥 midtermpart2 - Kyle Cordez.mp4
💡 Videos require Google Drive access. Open in new tab if it doesn't load.
🎥 midtermpart1 - Kyle Cordez.mp4
💡 Videos require Google Drive access. Open in new tab if it doesn't load.
📝 cordez_midterm_part2 - Kyle Cordez.txt
💡 Code is embedded in this portfolio - opens instantly in the live BrainImation editor (no internet required!)
📝 cordez_midterm_part1 - Kyle Cordez.txt
💡 Code is embedded in this portfolio - opens instantly in the live BrainImation editor (no internet required!)