参考 https://mne.tools/stable/auto_tutorials/index.html
伪影检测
伪影:
- 低频漂移
- 电源线噪声
- 心电(ECG)
- 眼电(EOG)
低频漂移
mne示例数据包括ssp投影,所以用del_proj()方法删除,plot()方法参数remove_dc为True是绘制数据时删除直流分量。
import os
import numpy as np
import mne
raw_name = 'H:/PycharmProjects/untitled1/simple/sample_audvis_raw.fif'
raw = mne.io.read_raw_fif(raw_name)
raw.crop(0, 60).load_data()
eeg_ch = mne.pick_types(raw.info, meg=False, eeg=True)
raw.plot(duration=60, order=meg_ch, n_channels=len(eeg_ch))
raw.copy().del_proj().plot(duration=60, order=meg_ch, n_channels=len(eeg_ch),remove_dc=False)
由上图可观察到漂移的半个周期差不多为10s,即频率为1/20Hz,分别用0.1Hz和0.2Hz进行高通滤波;