Home > brain, nirs > 3 ways to check NIRS data quality

3 ways to check NIRS data quality

August 31st, 2017

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: http://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 http://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!

Author: Xu Cui Categories: brain, nirs Tags:
Try Stork, a research tool we developed

Stork is a publication alert app developed by us at Stanford. As a researcher we often forget to follow up important publications - and it's practically impossible to search many keywords or researchers' names everyday. Stork can help us to search everyday and notifies us when there are new publications/grants. How Stork helped me?

About the author:

Xu Cui is a human brain research scientist in Stanford University. He lives in the Bay Area in the United States. He is also the founder of Stork (smart publication alert app), PaperBox and BizGenius.

 

He was born in He'nan province, China. He received education in Beijing University(BS), University of Tennessee (Knoxville) (MS), Baylor College of Medicine (PhD) and Stanford University (PostDoc). Read more ...
  1. No comments yet.
CAPTCHA Image