MNE学习笔记3

ICA修复伪影

导入存在明显ECG、EOG伪影的数据

import os
import mne
from mne.preprocessing import (ICA, create_eog_epochs,create_ecg_epochs, corrmap)
raw_name = 'H:/PycharmProjects/untitled1/simple/sample_audvis_raw.fif'

raw = mne.io.read_raw_fif(raw_name)
raw.crop(0, 60).load_data() 

#可视化存在伪影的数据
regexp = r'(MEG [12][45][123]1|EEG 00.)'
artifact_picks = mne.pick_channels_regexp(raw.ch_names, regexp=regexp)
raw.plot(order=artifact_picks, n_channels=len(artifact_picks))

eog_evoked = create_eog_epochs(raw).average()
eog_evoked.apply_baseline(baseline=(None, -0.2))
eog_evoked.plot_joint()
ecg_evoked = create_ecg_epochs(raw).average()
ecg_evoked