EEG信号——时域分析、频域分析、时频分析-程序员宅基地

技术标签: 算法  脑电数据处理  生物医学工程  脑机接口  

1.EEG信号的数据本质
1)原始的EEG数据是由很多个样本点数所构成的一个有限的离散的时间序列数据。至于样本点数的多少,则由采样率所决定,比如采样率为1000Hz,那么每秒就有1000个数据样本点。
2)每个样本点数据代表的是脑电波幅的大小,物理学上称为电压值,单位为伏特(V),由于脑电信号通常较弱,所以更常使用的单位为微伏(uV)。
2.主要出现区域
脑电波主要出现在颞叶和枕叶。 熟睡、 婴儿及严重器质性脑患者中都能观察到δ节律。波形为正弦波, 30~50μV。波幅沿枕-顶-额-颞降低; 正常情况下有明显调幅现象。 与视觉活动有关, 在觉醒安静闭目状态时出现的数量最多且振幅也最高。 当进入睡眠时, α节律完全消失。 α节律的频率、 振幅和空域分布等因素是反映大脑的机能状态的重要指标。 随着脑发育成熟或年龄的变化, α节律的数量逐渐增多,频率逐渐提高, 至成年期趋于稳定, 到老年期α节律则逐渐变慢。睡眠状态下, 浅睡眠期的背景节律为5~6Hz( θ波); 深睡眠期的背景节律为2~3Hz( δ波)。

3.脑电的产生机制

头表的EEG电位是由皮层椎体细胞群产生的, 它们排列整齐、 相互平行, 有利于总和形成较强的电场。 椎体细胞的有效方向垂直于皮层表面( Niedermeyer and Da Silva, 2005)
此外由于人体组织的低通滤波作用, 高频动作电位在传至头皮时将会极大衰减。 持续时间较长的突触后电位被认为是EEG信号产生的基础。动物实验研究发现, 局部场电位主要与皮层椎体细胞的突触后电位变化有关。 利用微电极记录, 可以观察到脑波的节律与兴奋性或抑制性突触后电位的变化一周。 与动作电位相比, 突触后电位振幅较低,但出现时间较长, 有利于通过多细胞的总和过程满足头表观测的需要。

3.信号的两个视角:

时域信号:幅值随时间变化的信号。
频域信号:幅值随频率变化的信号。

· 示例:
· 正弦/余弦波的时域表示、频域表示。
在这里插入图片描述

时域vs频域

时域分析: 脑电波幅随时间的变化(常 刺激后 - 刺激前)。
频域分析:分析脑电信号各频段的频谱能量

1)时域分析
因此,原始的EEG信号就是一个时域上的数据。时域分析关注的是EEG波幅随时间进程的变化情况,**事件相关电位(ERP)分析就是最常用的时域分析方法,能够快速得到由某个事件(刺激)所引起的波幅值的变化,如常见的事件相关电位成分P300,MMN、N400等。**时域分析的优点在于其计算简单和快速,而且由于不需要进行滤波处理,相较于频域分析/时频分析具有更高的时间精度和准确性。

2)频域分析

一般信号的时域、频域视角:

时域 -> 频域: 傅里叶级数分解\傅里叶变换**
频域->时域: 傅里叶逆变换

然而,时域分析并不足以全面反映EEG信号中所蕴藏的信息。我们常说的a波、B波、v波等就是根据不同频率的EEG信号来进行划分的,反映的是脑电信号随时间变化的快慢。
在这里插入图片描述

因此,频域分析实现的是EEG信号从时域到频域的转换,频域分析的结果为各个频率上的能量值分布,也就是我们常说的power值,可以简单理解为振幅的平方,单位为uV2.
在这里插入图片描述

频域分析使用的是傅里叶变换。根据傅里叶定理,任何连续测量的时序或信号,都可以表示为不同频率的正弦波信号的无限叠加。脑电信号可以看成就是由不同的正弦信号混壘而成的混合信号,通过傅里叶变换,就能够将这个混合信号重新分解成具有不同频率的正弦波,从而获得频域上的信息。
频域分析不仅可以用于分析任务态的数据,还常用于分析静息态的数据。
3)时频分析
但是
傅里叶变换有一个局限性,就是其只适用于稳态数据,而脑电数据就属于非稳态数据。
频域分析不能反映频率随时间的变化。因此,单一维度的时域分析或频域分析都不能全面反映信号特征,这时候就需要时频分析。

简而言之,脑电信号的频域分析

EEG频域分析方法:对各被试、各通道、各分段的脑电信号做傅里叶变换,得到各个频率点上的能量,再做统计分析。
主要适用场景:静息态数据等。
统计方式:
(1)沿着频率点统计:选择感兴趣的空间区域(通道),沿着频率点进行统计检验。
(2)沿着通道统计:选择感兴趣的频段,沿着通道进行统计检验
(3)提取特定通道、特定频段的能量数值,进行统计检验。

时域分析的缺点:仅能体现幅值随时间的变化。(缺少频段信息)
频域分析的缺点:仅能查看频段能量的分布,并不能体现特定频段能量随时间的变化。(缺少时间信息,仅适用于稳态数据分析)。

顾名思义,时频分析既包含时域,又包含了频域的信息,其方法是通过对脑电数据进行加窗处理,并假设在该时间窗内数据是稳态的,从而进行傅里叶变换,提取该时间窗内的频域信息。将窗口沿着时间轴向前滑动,并对每个时间窗内的数据进行同样的处理,这样就能得到随时间变化的频率的信息,所得到的结果就是时频图。其横轴代表的是时间,纵轴代表的是频率,每个时间频率所对应的点代表的就是power值。

在这里插入图片描述

窗口的大小会影响到时间精度和频率精度。窗口越大,时间精度越低,频率精度越高,适合分析低频慢波,窗口越小,时间精度越高,频率精度越低,适合分析高频快波。当然,具体还是根据研究需要进行设定。时间窗可以是固定大小,也可以具有自适应性。如短时傅里叶变换的时间窗就是大小固定的,而小波变换的时间窗则可以随着频率变化而伸缩,使用更灵活。

简而言之,,脑电信号的时频分析,在进行频谱分析时,保留时间信息。
、、
主要适用于:事件相关实验; 非稳态数据等。
主要实现方式:短时傅里叶变换、小波变换等。
分析方法:对各被试、各通道、各分段的数据做时频变换,得到各时频点的能量,再做统计分析。
统计方式:
(1)沿着频率点统计:选择感兴趣的空间区域(通道),沿着时频点进行统计检验。
(2)沿着通道统计:选择感兴趣的时频区域,沿着通道进行统计检验。
(3)提取特定通道、特定时频区域的能量数值,进行统计检验。

以下是现有文献做出的时频分析图示例:

在这里插入图片描述

在这里插入图片描述

在这里插入图片描述

在这里插入图片描述

总之,不同的分析关注的是脑电信号的不同维度,都有各自的优势与局限性,需要根据自己的研究需求选择合适的方法。

声明:转载自知乎:https://zhuanlan.zhihu.com/p/401681076
仅是自己整理理解,以及大家分享交流之用,无他意。
侵权删。

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/qq_38341276/article/details/124140782

智能推荐

稀疏编码的数学基础与理论分析-程序员宅基地

文章浏览阅读290次,点赞8次,收藏10次。1.背景介绍稀疏编码是一种用于处理稀疏数据的编码技术,其主要应用于信息传输、存储和处理等领域。稀疏数据是指数据中大部分元素为零或近似于零的数据,例如文本、图像、音频、视频等。稀疏编码的核心思想是将稀疏数据表示为非零元素和它们对应的位置信息,从而减少存储空间和计算复杂度。稀疏编码的研究起源于1990年代,随着大数据时代的到来,稀疏编码技术的应用范围和影响力不断扩大。目前,稀疏编码已经成为计算...

EasyGBS国标流媒体服务器GB28181国标方案安装使用文档-程序员宅基地

文章浏览阅读217次。EasyGBS - GB28181 国标方案安装使用文档下载安装包下载,正式使用需商业授权, 功能一致在线演示在线API架构图EasySIPCMSSIP 中心信令服务, 单节点, 自带一个 Redis Server, 随 EasySIPCMS 自启动, 不需要手动运行EasySIPSMSSIP 流媒体服务, 根..._easygbs-windows-2.6.0-23042316使用文档

【Web】记录巅峰极客2023 BabyURL题目复现——Jackson原生链_原生jackson 反序列化链子-程序员宅基地

文章浏览阅读1.2k次,点赞27次,收藏7次。2023巅峰极客 BabyURL之前AliyunCTF Bypassit I这题考查了这样一条链子:其实就是Jackson的原生反序列化利用今天复现的这题也是大同小异,一起来整一下。_原生jackson 反序列化链子

一文搞懂SpringCloud,详解干货,做好笔记_spring cloud-程序员宅基地

文章浏览阅读734次,点赞9次,收藏7次。微服务架构简单的说就是将单体应用进一步拆分,拆分成更小的服务,每个服务都是一个可以独立运行的项目。这么多小服务,如何管理他们?(服务治理 注册中心[服务注册 发现 剔除])这么多小服务,他们之间如何通讯?这么多小服务,客户端怎么访问他们?(网关)这么多小服务,一旦出现问题了,应该如何自处理?(容错)这么多小服务,一旦出现问题了,应该如何排错?(链路追踪)对于上面的问题,是任何一个微服务设计者都不能绕过去的,因此大部分的微服务产品都针对每一个问题提供了相应的组件来解决它们。_spring cloud

Js实现图片点击切换与轮播-程序员宅基地

文章浏览阅读5.9k次,点赞6次,收藏20次。Js实现图片点击切换与轮播图片点击切换<!DOCTYPE html><html> <head> <meta charset="UTF-8"> <title></title> <script type="text/ja..._点击图片进行轮播图切换

tensorflow-gpu版本安装教程(过程详细)_tensorflow gpu版本安装-程序员宅基地

文章浏览阅读10w+次,点赞245次,收藏1.5k次。在开始安装前,如果你的电脑装过tensorflow,请先把他们卸载干净,包括依赖的包(tensorflow-estimator、tensorboard、tensorflow、keras-applications、keras-preprocessing),不然后续安装了tensorflow-gpu可能会出现找不到cuda的问题。cuda、cudnn。..._tensorflow gpu版本安装

随便推点

物联网时代 权限滥用漏洞的攻击及防御-程序员宅基地

文章浏览阅读243次。0x00 简介权限滥用漏洞一般归类于逻辑问题,是指服务端功能开放过多或权限限制不严格,导致攻击者可以通过直接或间接调用的方式达到攻击效果。随着物联网时代的到来,这种漏洞已经屡见不鲜,各种漏洞组合利用也是千奇百怪、五花八门,这里总结漏洞是为了更好地应对和预防,如有不妥之处还请业内人士多多指教。0x01 背景2014年4月,在比特币飞涨的时代某网站曾经..._使用物联网漏洞的使用者

Visual Odometry and Depth Calculation--Epipolar Geometry--Direct Method--PnP_normalized plane coordinates-程序员宅基地

文章浏览阅读786次。A. Epipolar geometry and triangulationThe epipolar geometry mainly adopts the feature point method, such as SIFT, SURF and ORB, etc. to obtain the feature points corresponding to two frames of images. As shown in Figure 1, let the first image be ​ and th_normalized plane coordinates

开放信息抽取(OIE)系统(三)-- 第二代开放信息抽取系统(人工规则, rule-based, 先抽取关系)_语义角色增强的关系抽取-程序员宅基地

文章浏览阅读708次,点赞2次,收藏3次。开放信息抽取(OIE)系统(三)-- 第二代开放信息抽取系统(人工规则, rule-based, 先关系再实体)一.第二代开放信息抽取系统背景​ 第一代开放信息抽取系统(Open Information Extraction, OIE, learning-based, 自学习, 先抽取实体)通常抽取大量冗余信息,为了消除这些冗余信息,诞生了第二代开放信息抽取系统。二.第二代开放信息抽取系统历史第二代开放信息抽取系统着眼于解决第一代系统的三大问题: 大量非信息性提取(即省略关键信息的提取)、_语义角色增强的关系抽取

10个顶尖响应式HTML5网页_html欢迎页面-程序员宅基地

文章浏览阅读1.1w次,点赞6次,收藏51次。快速完成网页设计,10个顶尖响应式HTML5网页模板助你一臂之力为了寻找一个优质的网页模板,网页设计师和开发者往往可能会花上大半天的时间。不过幸运的是,现在的网页设计师和开发人员已经开始共享HTML5,Bootstrap和CSS3中的免费网页模板资源。鉴于网站模板的灵活性和强大的功能,现在广大设计师和开发者对html5网站的实际需求日益增长。为了造福大众,Mockplus的小伙伴整理了2018年最..._html欢迎页面

计算机二级 考试科目,2018全国计算机等级考试调整,一、二级都增加了考试科目...-程序员宅基地

文章浏览阅读282次。原标题:2018全国计算机等级考试调整,一、二级都增加了考试科目全国计算机等级考试将于9月15-17日举行。在备考的最后冲刺阶段,小编为大家整理了今年新公布的全国计算机等级考试调整方案,希望对备考的小伙伴有所帮助,快随小编往下看吧!从2018年3月开始,全国计算机等级考试实施2018版考试大纲,并按新体系开考各个考试级别。具体调整内容如下:一、考试级别及科目1.一级新增“网络安全素质教育”科目(代..._计算机二级增报科目什么意思

conan简单使用_apt install conan-程序员宅基地

文章浏览阅读240次。conan简单使用。_apt install conan