3 ways to check NIRS data quality

1 min read

Before performing any data analysis, we should check the data quality first. Below are 3 ways to do so.

1. Visual check of the time series

The best pattern detector is our eyes and brains! In many cases, if we visually see the data, we know what is wrong. You may use the method and program in this post to plot the time courses of all channels (not just one).

Visual data quality check of NIRS time courses
Visual data quality check of NIRS time courses

In the plot above, all 48 channels are plotted together (the y-axis). We can clearly see two types of noise:

  1. The spikes which occur in most channels after time point 7000. These spikes are caused by head motion.
  2. The high noise level in the “red” channels (channels 39, 34, 31). This is more evident if we plot the variance of each channel (figure below). As we can easily see, the variance of channels 39, 34 and 31 is much higher than other channels.

    Variance vs Channel
    Variance vs Channel

2. Existence of the “heart-beat” band

If the NIRS signal was acquired well, then the heart beating signal should be captured, leaving a bright brand in the frequency ~1Hz in the wavelet transform plot, just like the left plot in the figure below (the band close to period 8). If there is no such band, it does not necessarily mean the signal is trash, but you need to be cautious. To use the wavelet transform toolbox, please download here: https://www.alivelearn.net/?p=1561


NIRS wavelet
NIRS wavelet
Example:
figure;wt(hbo(:,1))

3. Correlation between hbo and hbr

The third way is to check the correlation between hbo and hbr. They are supposed to have negative correlation, at least in young healthy subjects. If not, or if they have perfect negative correlation (-1), then they might contain too much noise. We have a separate article on this method. Please check out https://www.alivelearn.net/?p=1767

Correlation between oxy and deoxy-Hb
Correlation between oxy and deoxy-Hb

Below is the scripts used for the 3 methods.

[hbo,hbr,mark]=readHitachData('SA06_MES_Probe1.csv');

figure;plotTraces(hbr,1:52,mark)

figure;wt(hbo(:,1))

for ii=1:52; wt(hbo(:,ii)); pause; end

[badchannels] = checkDataQuality(hbo,hbr);

Do you have other ways to check data quality? Please let me know!



写作助手,把中式英语变成专业英文


Want to receive new post notification? 有新文章通知我

第五十四期fNIRS Journal Club通知2024/07/27, 10am 黄鑫

电子烟作为一种新型的尼古丁输送装置,常以安全健康、有助于戒烟的形象出现在公众视野中,使得大批消费者转而投入到此类产品的怀抱中。然而,相关研究却仍存较大短板,深入了解电子烟对行为与脑功能活动的影响将更好
Wanling Zhu
10 sec read

第五十三期fNIRS Journal Club视频 李洪

Youtube: https://youtu.be/czyJTrXhees优酷:https://v.youku.com/v_show/id_XNjQwNjc5OTk0OA==.html 个体在处理不同
Wanling Zhu
12 sec read

第五十三期fNIRS Journal Club通知2024/06/22, 10am 李洪

个体在处理不同记忆负荷信息时会表现出一定的行为差异。作为一项新兴指标,瞬时脑信号变异性能够揭示个体内部因任务需求不断变化而进行的神经资源分配,从而为了解大脑如何适应和处理不同复杂程度的信息提供了新的见
Wanling Zhu
8 sec read

Leave a Reply

Your email address will not be published. Required fields are marked *