Python 中 MNE 读取EEG竞赛数据绘图和提取epoch（gdf格式）

上一篇文章中介绍了研究读取gdf格式数据的心路历程，绕了半天竟然直接操作就可以。本来以为会很容易，毕竟数据都已经成功读入了，可是我又错了，弄了好久读取的数据绘图都不对（打印看数据是对的），感觉是mne的一个bug，在Github上提issue可是网络不给力，图片死活传不上，最后给开发团队发了封邮件，也不知道能不能回我。。。。好在经过不懈努力，试了多种方法，总算有一种方法可以成功绘制数据图像了。

问题描述，读取数据绘图，可以明显看出EEG数据绘制的图像有问题。

rawDataGDF = mne.io.read_raw_gdf("B0101T.gdf",preload=True,eog=['EOG:ch01', 'EOG:ch02', 'EOG:ch03'])

rawDataGDF.plot()

打印输出数据查看，并没有发现异常。

又尝试了N多的方法，也没有定位到问题的根本原因，虽然严重怀疑是描述信息（info）有问题，可是比对了好久也没有发现异常。最后没办法了尝试自己构造RawArray数据结构，虽然麻烦一些，也是最后的办法了，好在这个方法可行。

sfreq = 250

data = np.squeeze(np.array([rawDataGDF['EEG:Cz'][0], rawDataGDF['EEG:C3'][0], rawDataGDF['EEG:C4'][0]]))

ch_types = ['eeg', 'eeg', 'eeg']

ch_names = ['EEG:Cz', 'EEG:C3', 'EEG:C4',]

info = mne.create_info(ch_names=ch_names,sfreq=sfreq,ch_types=ch_types)

raw = mne.io.RawArray(data, info)

raw.plot(n_channels=3,show=True)

最后将代码整理了下，用于处理EEG竞赛的数据，嗯嗯，真香~~

import mne
import matplotlib.pyplot as plt
import numpy as np
%matplotlib auto
#文件路径和文件名称
path = "F:\\BaiduNetdiskDownload\\BCICompetition\\BCICIV_2b_gdf\\"
fileName = "B0101T"
#事件的对应关系
eventDescription = {'276': "eyesOpen", '277': "eyesClosed", '768': "startTrail", '769':"cueLeft", '770':"cueRight", '781':"feedback", '783':"cueUnknown",'1023': "rejected", '1077': 'horizonEyeMove', '1078': "verticalEyeMove", '1079':"eyeRotation",'1081':"eyeBlinks",'32766':"startRun"}
# 读取原始数据
rawDataGDF = mne.io.read_raw_gdf(path+fileName+".gdf",preload=True,eog=['EOG:ch01', 'EOG:ch02', 'EOG:ch03'])
# 定义数据通道和通道类型
ch_types = ['eeg', 'eeg', 'eeg','eog','eog','eog']
ch_names = ['EEG:Cz', 'EEG:C3', 'EEG:C4','EOG:ch01', 'EOG:ch02', 'EOG:ch03']
# 创建数据的描述信息
info = mne.create_info(ch_names=ch_names, sfreq=rawDataGDF.info['sfreq'], ch_types=ch_types)
# 创建数据结构体
data = np.squeeze(np.array([rawDataGDF['EEG:Cz'][0], rawDataGDF['EEG:C3'][0], rawDataGDF['EEG:C4'][0], rawDataGDF['EOG:ch01'][0] , rawDataGDF['EOG:ch02'][0] ,rawDataGDF['EOG:ch03'][0]]))
# 创建RawArray类型的数据
rawData = mne.io.RawArray(data, info)
rawData.plot()
# 获取事件
event,_= mne.events_from_annotations(rawDataGDF)
print(event)
event_id = {}
for i in _:#整理event_idevent_id[eventDescription[i]] = _[i]
# 提取epoch
epochs = mne.Epochs(rawData,event,event_id, tmax=1.5 ,event_repeated = 'merge')
epochs.plot()
ev_left = epochs['cueLeft']
ev_left.plot()

原始数据图像

epochs数据图像

左手对应得epoch

Python 中 MNE 读取EEG竞赛数据绘图和提取epoch（gdf格式）相关推荐

Python 中 MNE 读取EEG竞赛数据（gdf格式）
最近在研究运动想象相关的算法,发现网上有一些公开的脑电数据,发现之前EEG还有这方面的竞赛,不知道为啥现在没有了,或许还有只不过我到目前还没找到,哈哈.前半部分都是说如何获取数据,后面才是标题的重点内 ...
Python中ArcPy读取Excel表格数据创建矢量要素图层并生成属性表字段与内容
1 任务需求首先,我们来明确一下本文所需实现的需求. 现有一个记录北京市部分PM2.5浓度监测站点信息的Excel表格数据,格式为.xls:文件内包含站点编号.X与Y坐标.站点名称等四列数据 ...
Python中如何读取xml的数据
<?xml version="1.0" encoding="utf-8" ?> - <catalog><maxid>4< ...
Python 中MNE库去伪迹(ICA)案例的逐句解析
本文是在阅读博主zhoudapeng01的文章Python 中MNE库去伪迹(ICA)_zhoudapeng01的博客-CSDN博客_ica mne后做的一个翻译工作. 主要介绍的还是ICA,需要了解 ...
python应用中调用spark_在python中使用pyspark读写Hive数据操作
1.读Hive表数据 pyspark读取hive数据非常简单,因为它有专门的接口来读取,完全不需要像hbase那样,需要做很多配置,pyspark提供的操作hive的接口,使得程序可以直接使用SQL语 ...
python处理netcdf_在python中高效读取netcdf变量
我需要能够在 python中快速读取大量的netCDF变量(每个文件1个变量).我发现netCDF4库中的数据集功能与其他语言(例如IDL)中的读取实用程序相比相当慢. 我的变量的形状为(2600,5 ...
从Mysql某一表中随机读取n条数据的SQL查询语句
若要在i ≤ R ≤ j 这个范围得到一个随机整数R ,需要用到表达式 FLOOR(i + RAND() * (j – i + 1)).例如, 若要在7 到 12 的范围(包括7和12)内得到一个随机 ...
python读取ansi编码文件,如何在Python中同时读取ANSI和Unicode txt文件？
我是python新手,遇到了一个奇怪的问题: 当一个目录中有50个txt文件时,我希望读取每个.txt文件并将其内容保存在一个唯一的变量中,例如:**file = open(fcf[i], 'r') ...
【机器学习】python使用matplotlib进行二维数据绘图并保存为png图片
端到端机器学习导航: [机器学习]python借助pandas加载并显示csv数据文件,并绘制直方图 [机器学习]python使用matplotlib进行二维数据绘图并保存为png图片 [机器学习]p ...

Python 中 MNE 读取EEG竞赛数据绘图和提取epoch（gdf格式）

Python 中 MNE 读取EEG竞赛数据绘图和提取epoch（gdf格式）相关推荐

最新文章

热门文章