CN105249962B - 头皮脑电信号回顾性癫痫发作点检测方法及系统 - Google Patents
头皮脑电信号回顾性癫痫发作点检测方法及系统 Download PDFInfo
- Publication number
- CN105249962B CN105249962B CN201510736832.4A CN201510736832A CN105249962B CN 105249962 B CN105249962 B CN 105249962B CN 201510736832 A CN201510736832 A CN 201510736832A CN 105249962 B CN105249962 B CN 105249962B
- Authority
- CN
- China
- Prior art keywords
- eeg signals
- epileptic attack
- retrospective
- attack point
- sample entropy
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明属于头皮脑电信号技术领域,提出了一种头皮脑电信号回顾性癫痫发作点检测方法及系统。本发明方法是对去除了各种伪迹脑电信号,通过非线性动力学样本熵阈值检测法,进行回顾性分析确定癫痫发作点的。本发明的头皮脑电信号回顾性癫痫发作点检测系统,包括脑电信号接收模块、癫痫发作点确定模块,信息输出模块。其中,脑电信号接收模块用于接收临床采集到的原始脑电信号。癫痫发作点确定模块用于通过脑电信号接收模块接收的脑电信号分析确定回顾性癫痫发作点。信息输出模块用于将癫痫发作点确定模块确定的回顾性癫痫发作点输出。采用本发明方法或者系统脑电信号数据可在10秒内完成解混,快速确定癫痫发作点,效果显著。
Description
技术领域
本发明涉及头皮脑电信号技术领域,尤其涉及头皮脑电信号回顾性癫痫发作点检测方法及系统。
背景技术
癫痫是一种神经系统紊乱的疾病,据文献报道该病在人群中的发病率约为0.5%~2%,以大脑神经细胞群反复超同步放电引起的发作性、突然性、短暂性脑功能紊乱为特征,为患者的生活带来极大的不便。脑电图(EEG)检查是目前临床上癫痫诊断和病灶定位的常用技术,通过脑电图分析癫痫,容易为医生和患者所接受。脑电图信号的时间分辨率高,能在毫秒级水平上精确地反映大脑的时变特性。而24小时长程脑电信号是癫痫诊断的重要依据,但面对该海量数据,目前依然靠人工目测搜寻发作点进行下一步的分析,工作量巨大且误判率高。
但由于人体是一个复杂的网络系统,脑电图中不可避免的会受到一些无关的电生理信号如眼电、肌电和心电的干扰,因此临床采集到的脑电信号中经常会含有大量的伪迹影响,而且由于大脑神经网络是互联网络,不同通道的脑电信号也会存在相互串扰,使得每一通道的脑电最终都表现为一种混叠信号。串扰现象给脑电的分析带来不便,甚至有可能可出错误的结果。目前已有的一些癫痫定位检测方法会因为信号质量的影响而降低准确性。
对于脑电信号中信息的研究方法有时域频域分析与概率统计分析,但时域的方法如棘波、尖波的检测,概率统计的方法如神经网络、主成分分析,都无法涵盖癫痫脑电的大部分特征。
发明内容
本发明需要解决的问题是:针对上述缺陷,本发明提出了一种头皮脑电信号回顾性癫痫发作点检测方法,包括采集原始脑电信号,获得去除各种伪迹脑电信号;
对去除伪迹的脑电信号通过非线性动力学样本熵阈值检测法,进行回顾性分析确定癫痫发作点。
优选的是,所述获得去除各种伪迹脑电信号方法如下:
(1)对脑电信号进行去均值化处理;
(2)对去均值化处理后的脑电信号进行白化处理;
(3)求解最佳分离矩阵;
(4)将步骤(2)获得的电脑信号与步骤(3)得到的最佳分离矩阵相乘,获得去除各种伪迹脑电信号。
上述任意实施方式中优选的是,所述脑电信号均值化处理如下:
X=(Xj-mj)/Sj
其中,X为脑电信号均值化处理后的脑电信号,Xj为原始脑电信号,mj和Sj分别为每个脑电信号样本的均值和标准差。
上述任意实施方式中优选的是,所述白化处理是对任意多维信号施加一个线性变换,使多维信号变为白色信号的处理过程。
上述任意实施方式中优选的是,所述求解最佳分离矩阵方法如下:建立优化算法与求解判别函数,用判别函数判断优化算法何时收敛,求解出分离矩阵。
上述任意实施方式中优选的是,所述判别函数如下:
其中,其中W是待求的解混矩阵,H(Y)随机变量Y的熵,J(Y)随机变量Y的负熵,Ygauss是n个高斯随机变量组成的随机矢量,与随机变量Y有相同的均值和协方差阵。
上述任意实施方式中优选的是,所述优化算法是混合编码的遗传算法。
上述任意实施方式中优选的是,所述混合编码的遗传算法具体如下:
(1)随机产生矩阵的初始种群,将随机生成的种群按实数编码的方式进行编码,
计算初始种群的适应度函数即负熵最大化的代价函数,
对适应度值排序,进行种群的实数交叉和变异操作,进入迭代运算;
(2)达到指定的遗传代数或者负熵的变化达到阈值,停止遗传操作,获得种群用格雷码;
(3)对得到的种群用格雷码进行重新编码;
(4)对适应度值排序,进行种群的二进制交叉和变异操作,进入迭代运算;
(5)达到指定的遗传代数或者负熵的变化达到阈值,停止遗传操作,获得最佳分离矩阵。
上述任意实施方式中优选的是,所述步骤(2)与步骤(5)中的阈值均是经验值。
上述任意实施方式中优选的是,所述确定癫痫发作点具体如下:
取样本熵SampEn(m,r,N)为-ln(A/B),当本熵达到指定阈值时,确定该点为癫痫发作点;
其中,B={[(N-m-1)(N-m)]/2}Bm(r),A={[(N-m-1)(N-m)]/2}Am(r),N为计算样本熵所取用的窗口长度值,m为组成矢量矩阵时所使用的矩阵维度,r为相似性容限,Bm(r)为某个窗口中符合相似性容限的数据平均值,Bm(r)为某个窗口中符合相似性容限的数据平均值,Am(r)为相邻窗口中符合相似性容限的数据平均值。
上述任意实施方式中优选的是,所述阈值为15。
上述任意实施方式中优选的是,所述相似性容限r取信号0.15倍的标准差。
为解决上述技术问题,本发明又提出了一种头皮脑电信号回顾性癫痫发作点检测系统,包括脑电信号接收模块、癫痫发作点确定模块,信息输出模块,其中,脑电信号接收模块用于接收临床采集到的原始脑电信号,癫痫发作点确定模块用于通过脑电信号接收模块接收的脑电信号分析确定回顾性癫痫发作点,信息输出模块用于将癫痫发作点确定模块确定的回顾性癫痫发作点输出。
优选的是,所述癫痫发作点确定模块包括脑电信号伪迹去除模块与发作点分析确定模块。
上述任意实施方式中优选的是,所述脑电信号伪迹去除模块包括去均值化处理模块、白化处理模块、求解最佳分离矩阵模块、获得无伪迹脑电信号模块。
上述任意实施方式中优选的是,所述信息输出模块包括显示接口输出、音频接口输出、网络接口输出。
本发明的头皮脑电信号回顾性癫痫发作点检测方法及系统,通过脑电信号的负熵及两阶段遗传优化算法快速实现盲源分离,由于脑电信号数据量及其庞大,一般算法均无法做到实时运算,本发明由于使用了遗传算法,优化效果显著,进过测试20秒24导联的脑电数据可在10秒内完成解混,快速确定癫痫发作点,效果显著。
本发明的头皮脑电信号回顾性癫痫发作点检测方法及系统,通过脑电的非线性动力学样本熵阈值检测法,对患者海量脑电信号进行回顾性分析,实现癫痫监测的全程自动化。
本发明的头皮脑电信号回顾性癫痫发作点检测方法及系统中,所使用的样本熵是条件概率的负平均自然对数的精确值,因此它的计算不依赖数据长度,具有良好的广泛适应性。
附图说明
图1为按照本发明的头皮脑电信号回顾性癫痫发作点检测方法在实施例中的流程图。
图2为按照本发明的头皮脑电信号回顾性癫痫发作点检测系统在实施例中的组成示意图。
具体实施方式
下面结合附图对本发明作进一步详细描述,有必要在此指出的是,以下具体实施方式只用于对本进行进一步的说明,不能理解为对本发明保护范围的限制,该领域的技术人员可以根据上述发明内容对本发明作出一些非本质的改进和调整。
如图1所示,本发明的头皮脑电信号回顾性癫痫发作点检测方法,包括如下步骤:
步骤一、获取去除各种伪迹脑电信号
通过脑电信号的负熵及两阶段遗传优化算法实现快速盲源分离,去除脑电信号
各种伪迹,尤其是心电噪声的影响。盲源快速分离算法可大致描述为:
S-->混合矩阵A-->X=AS-->解混矩阵W-->Y=WX
其中,S为理想的独立源信号,X为实际可以采集到的观测信号(混合信号),可以看作是独立源信号和混合矩阵A的乘积,盲源快速分离算法就是要找到一个解混矩阵W,让采集信号通过后得到分离信号Y,以期望Y能尽量接近S。因为盲源快速分离算法无法区分噪声和独立源信号,所以采用该方法对信号进行盲源分离时,需要对观测数据进行去均值化和白化两项预处理,用以简化盲源快速分离算法,而且对比试验中表明,经过上述两项预处理的数据,在优化迭代的时候,能体现出更加优良的收敛性能。
如图1所示,获取去除各种伪迹脑电信号具体过程如下:
1、对脑电信号进行去均值化处理。
对采集到的脑点信号Xj=[xj1(t),xj2(t),…,xjn(t)],进行去均值化处理,即X=(Xj-mj)/Sj,其中,mj和Sj分别为每个脑电信号样本的均值和标准差。去均值化的目的是从脑电信号中减去信号的均值向量,使得脑电信号成为零均值矢量。
2、对去均值化处理后的脑电信号进行白化处理
白化就是对任意多维信号施加一个线性变换,使其变为白色信号的处理过程,白化也称为球化或者归一化解相关,对应的变换矩阵称为白化矩阵,白化后的矩阵各分量互不相关,具有独立性且具有单位方差。
令Q为X的白化矩阵,则有cov(QX)=I,其中,I为单位矩阵。再将X=AS代入上式并令矩阵B=QA,则cov(BS)=I。由于S是白色的随机向量,因而矩阵B一定是正交矩阵。白化就是使原来的混合矩阵A简化为一个正交矩阵B。可见白化可以有效地降低问题的复杂度。
3、求解判别函数
快速盲源分离算法的基本目标就是要找到一个线性变换矩阵W,W也称为解混矩阵或分离矩阵,使变换后的各信号之间尽可能统计独立。可以分解为两个基本问题:建立最优化算法和判别函数,用判别函数判断优化算法何时收敛并求解出分离矩阵W。
由信息论理论可知,在所有等方差的随机变量中,高斯变量的熵最大。因而利用熵来度量非高斯性,为得到一种对高斯分布为零、总是非负的非高斯性的度量,常用熵的修正形式,即负熵。
设随机变量Y的概率密度函数为p(Y),则它的熵为:H(Y)=-∫pY(y)lgPY(y)dy,
信号的负熵定义如下:
J(Y)=H(Ygauss)-H(Y),
其中Ygauss是n个高斯随机变量组成的随机矢量,与Y有相同的均值和协方差阵。可见负熵总是非负。负熵完全可以度量信号之间的独立性。
输出的互信息可以用负熵表示为:最小化互信息就等于最大化负熵。因而基于负熵最大化的代价函数为:
其中,其中W是待求的解混矩阵,H(Y)随机变量Y的熵,J(Y)随机变量Y的负熵,Ygauss是n个高斯随机变量组成的随机矢量,与随机变量Y有相同的均值和协方差阵。
上述的代价函数即确定为判别函数。
4、求解最佳分离矩阵
用基于负熵最大化的代价函数作为遗传算法的适应度函数,可用于判断是否达到收敛条件。
混合编码的遗传算法:是本发明的头皮脑电信号回顾性癫痫发作点检测方法采用的优化算法,具体步骤如下
(1)随机产生矩阵W=[w1,w2,...wn]T的初始种群;
将随机生成的种群按实数编码的方式进行编码。
计算初始种群的适应度函数-负熵最大化的代价函数;
对适应度值排序,进行种群的实数交叉和变异操作,进入迭代运算;
(2)达到指定的遗传代数或者负熵的变化达到阈值,停止遗传操作,获得种群用格雷码。阈值是一个经验值,一般为一个很小的数,代表几乎没有变化。
(3)对得到的种群用格雷码进行重新编码。
(4)对适应度值排序,进行种群的二进制交叉和变异操作,进入迭代运算;
(5)达到指定的遗传代数或者负熵的变化达到阈值,停止遗传操作,获得最佳分离矩阵。
上述步骤中的阈值是一个经验值,一般为一个很小的数,代表几乎没有变化。
5、获得去除各种伪迹脑电信号
将去均值化和白化后的脑电信号X=[x1(t),x2(t),…,xn(t)]和最佳矩阵W=[w1,w2,...wn]T相乘,解混获得无伪迹信号的干扰的脑电信号X'=[x1(t)',x2(t)',…,xn(t)']。
步骤二、对去除伪迹的脑电信号通过非线性动力学样本熵阈值检测法,进行回顾性分析确定癫痫发作点
将去处各种伪迹的脑电信号X'=[x1(t)',x2(t)',…,xn(t)']通过非线性动力学样本熵阈值检测法进行回顾性分析,在长程脑电数据中找到癫痫发作的数据段,提示可能出现癫痫发作,以辅助医生的完成海量数据筛选和分析,实现癫痫监测的全程自动化。具体实施步骤如下:
选取去处各种伪迹的脑电信号计算样本熵,该数据要求包括发作前期和发作期两部分数据。样本熵曲线是由多个样本熵点绘制而成,数据长度D=样本熵个数M*窗口大小N。下面介绍单个样本熵的计算。
设定窗口的大小N,即这该段信号即为[x1(t)',x2(t)',…,xN(t)']x(1),x(2),x(3)……x(N),按顺序组成一组m维矢量,Xm,i'(t)=[x'm(t),x'm+1(t),……,x'm+i-1(t)],1<i<N-m。参数m的设定在本实验中以经验值来确定。
定义矢量X'm,i(t)与X'm,j(t)之间的距离d[X'm,i(t),X'm,j(t)],为两者对应元素中差值最大的一个,即:d[X'm,i(t),X'm,j(t)]=max[|x(i+k)-x(j+k)|],0<k<m-1,i<j,1<i,j<N-m.
定义阈值r为脑电信号序列的相似容限,常用序列的标准差SD来表示,则r=0.15*std(N)。对于给定阈值r,对每一个i值统计d[X'm,i(t),X'm,j(t)]小于r的数目(模板匹配数),并计算该数目与距离总数的比值,用表示,即:
的数目,i,j}。同时用Bm(r)表示其对于所
有的i的平均值,即:
维数加1并不断重复上述步骤,计算得到与Am(r)。这该序列的样本熵为:SampEn(m,r)=limN{-ln[Am(r)/Bm(r)]},当窗口N确定时,即为有限值时,SampEn(m,r,N)=-ln[Am(r)/Bm(r)]。对取定的模板向量参数,当模板长度为m,阈值为r时,设B为此时模板匹配的总数,即:B={[(N-m-1)(N-m)]/2}Bm(r)。当维数加1时,对于此时的阈值r,设A为此时的模板匹配的总数,即:A={[(N-m-1)(N-m)]/2}Am(r)。取样本熵SampEn(m,r,N)为-ln(A/B)。当样本熵达到指定阈值时,记录该点为发作点。具体是:由于本算法为回顾性算法,经多次对比实验,本发明的头皮脑电信号回顾性癫痫发作点检测方法中此处的阈值设置为15时,预测准确率最高时,记录该点为发作点。
如图2所示,本发明的头皮脑电信号回顾性癫痫发作点检测系统,包括脑电信号接收模块、癫痫发作点确定模块,信息输出模块。其中,脑电信号接收模块用于接收临床采集到的原始脑电信号。癫痫发作点确定模块用于通过脑电信号接收模块接收的脑电信号分析确定回顾性癫痫发作点。信息输出模块用于将癫痫发作点确定模块确定的回顾性癫痫发作点输出。
癫痫发作点确定模块包括脑电信号伪迹去除模块、发作点分析确定模块,脑电信号伪迹去除模块包括去均值化处理模块、白化处理模块、求解最佳分离矩阵模块、获得无伪迹脑电信号模块。其中,均值化处理模块是对原始脑电信号进行均值处理。白化处理模块是对均值化处理模块处理后的脑电信号进行白化处理。求解最佳分离矩阵模块用于求解最佳分离矩阵,首先求解代价函数,求解代价函数作为遗传算法的适应度函数,可用于判断是否达到收敛条件,最终求得最佳分离矩阵。获得无伪迹脑电信号模块是将白化处理模块的到的脑电信号与求解最佳分离矩阵模块获得的最佳分离矩阵相乘获得无伪迹脑电信号。发作点分析确定模块是对去除伪迹的脑电信号通过非线性动力学样本熵阈值检测法,进行回顾性分析确定癫痫发作点。
信息输出模块将确定的癫痫发作点可以以不同方式输出。包括显示接口输出、音频接口输出、网络接口输出。当确定为发作点时,在显示上以红色竖线标记方式通过显示接口输出至显示器,同时通过音频接口输出声音报警,还可通过网络接口通过3G网络或wifi无线网发送报警信息。
Claims (16)
1.一种头皮脑电信号回顾性癫痫发作点检测方法,包括采集原始脑电信号,其特征在于:
获得去除各种伪迹脑电信号;
对去除伪迹的脑电信号通过非线性动力学样本熵阈值检测法,进行回顾性分析确定癫痫发作点;
所述非线性动力学样本熵阈值检测法具体包括如下步骤:
(1)、选取去除各种伪迹的脑电信号,所述脑电信号要求包括发作前期和发作期两部分数据;
(2)、根据所述脑电信号计算单个的样本熵;
(3)、多个样本熵点绘制成样本熵曲线,数据长度D=样本熵个数M*窗口大小N;
(4)、定义阈值,当样本熵达到指定阈值时,记录该点为发作点。
2.根据权利要求1所述的一种头皮脑电信号回顾性癫痫发作点检测方法,其特征在于:所述获得去除各种伪迹脑电信号方法如下:
(1)对脑电信号进行去均值化处理;
(2)对去均值化处理后的脑电信号进行白化处理;
(3)求解最佳分离矩阵;
(4)将步骤(2)获得的脑电信号与步骤(3)得到的最佳分离矩阵相乘,获得去除各种伪迹脑电信号。
3.根据权利要求2所述的一种头皮脑电信号回顾性癫痫发作点检测方法,其特征在于:所述脑电信号均值化处理如下:
X=(Xj-mj)/Sj
其中,X为脑电信号均值化处理后的脑电信号,Xj为原始脑电信号,mj和Sj分别为每个脑电信号样本的均值和标准差。
4.根据权利要求2所述的一种头皮脑电信号回顾性癫痫发作点检测方法,其特征在于:所述白化处理是对任意多维信号施加一个线性变换,使多维信号变为白色信号的处理过程。
5.根据权利要求2所述的一种头皮脑电信号回顾性癫痫发作点检测方法,其特征在于:所述求解最佳分离矩阵方法如下:建立优化算法与求解判别函数,用判别函数判断优化算法何时收敛,求解出分离矩阵。
6.根据权利要求5所述的一种头皮脑电信号回顾性癫痫发作点检测方法,其特征在于:所述判别函数如下:
其中,其中W是待求的解混矩阵,H(Y)随机变量Y的熵,J(Y)随机变量Y的负熵,Ygauss是n个高斯随机变量组成的随机矢量,与随机变量Y有相同的均值和协方差阵,N为窗口大小,M为样本熵个数。
7.根据权利要求5所述的一种头皮脑电信号回顾性癫痫发作点检测方法,其特征在于:所述优化算法是混合编码的遗传算法。
8.根据权利要求7所述的一种头皮脑电信号回顾性癫痫发作点检测方法,其特征在于:所述混合编码的遗传算法具体如下:
(1)随机产生矩阵的初始种群,将随机生成的种群按实数编码的方式进行编码,
计算初始种群的适应度函数即负熵最大化的代价函数,
对适应度值排序,进行种群的实数交叉和变异操作,进入迭代运算;
(2)达到指定的遗传代数或者负熵的变化达到阈值,停止遗传操作,获得种群用格雷码;
(3)对得到的种群用格雷码进行重新编码;
(4)对适应度值排序,进行种群的二进制交叉和变异操作,进入迭代运算;
(5)达到指定的遗传代数或者负熵的变化达到阈值,停止遗传操作,获得最佳分离矩阵。
9.根据权利要求8所述的一种头皮脑电信号回顾性癫痫发作点检测方法,其特征在于:
所述混合编码的遗传算法中的步骤(2)与步骤(5)中的阈值均是经验值。
10.根据权利要求1所述的一种头皮脑电信号回顾性癫痫发作点检测方法,其特征在于:单个样本熵的计算包括如下步骤:
(1)设定窗口的大小N,该段信号即为[x1(t)',x2(t)',…,xN(t)']x(1),x(2),x(3)……x(N),按顺序组成一组m维矢量,Xm,i'(t)=[x'm(t),x'm+1(t),......,x'm+i-1(t)],其中:1<i<N-m;
(2)定义矢量X'm,i(t)与X'm,j(t)之间的距离d[X'm,i(t),X'm,j(t)],为两者对应元素中差值最大的一个,即:d[X'm,i(t),X'm,j(t)]=max[|x(i+k)-x(j+k)|],其中,0<k<m-1,i<j,1<i,j<N-m;
(3)定义阈值r为脑电信号序列的相似容限,用序列的标准差SD来表示,则r=0.15*std(N);对于给定阈值r,对每一个i值,统计d[X'm,i(t),X'm,j(t)]小于r的数目,并计算该数目与距离总数的比值,用表示,即:
同时用Bm(r)表示其对于所有的i的平均值,即:
(4)维数加1并不断重复上述步骤,计算得到与Am(r);则该序列的样本熵为:SampEn(m,r)=limN{-ln[Am(r)/Bm(r)]},当窗口N确定时,即为有限值时,SampEn(m,r,N)=-ln[Am(r)/Bm(r)];
其中N为计算样本熵所取用的窗口长度值,m为组成矢量矩阵时所使用的矩阵维度,r为相似性容限,Bm(r)为某个窗口中符合相似性容限的数据平均值,Bm(r)为某个窗口中符合相似性容限的数据平均值,Am(r)为相邻窗口中符合相似性容限的数据平均值。
11.根据权利要求10所述的一种头皮脑电信号回顾性癫痫发作点检测方法,其特征在于:所述确定癫痫发作点具体如下:
取样本熵SampEn(m,r,N)为-ln(A/B),当所述样本熵达到指定阈值时,确定该点为癫痫发作点;
其中,B={[(N-m-1)(N-m)]/2}Bm(r),A={[(N-m-1)(N-m)]/2}Am(r),
其中N为计算样本熵所取用的窗口长度值,m为组成矢量矩阵时所使用的矩阵维度,r为相似性容限,Bm(r)为某个窗口中符合相似性容限的数据平均值,Bm(r)为某个窗口中符合相似性容限的数据平均值,Am(r)为相邻窗口中符合相似性容限的数据平均值。
12.根据权利要求11所述的一种头皮脑电信号回顾性癫痫发作点检测方法,其特征在于:所述指定阈值为15。
13.根据权利要求11所述的一种头皮脑电信号回顾性癫痫发作点检测方法,其特征在于:所述相似性容限r取信号0.15倍的标准差。
14.一种头皮脑电信号回顾性癫痫发作点检测系统,其特征在于:包括脑电信号接收模块、癫痫发作点确定模块,信息输出模块,其中,脑电信号接收模块用于接收临床采集到的原始脑电信号,癫痫发作点确定模块用于通过脑电信号接收模块接收的脑电信号分析确定回顾性癫痫发作点,信息输出模块用于将癫痫发作点确定模块确定的回顾性癫痫发作点输出;
所述癫痫发作点确定模块包括脑电信号伪迹去除模块与发作点分析确定模块;
癫痫发作点确定模块对去除伪迹的脑电信号通过非线性动力学样本熵阈值检测法,进行回顾性分析确定癫痫发作点;
所述非线性动力学样本熵阈值检测法具体包括如下步骤:
(1)、选取去除各种伪迹的脑电信号,所述脑电信号要求包括发作前期和发作期两部分数据;
(2)、根据所述脑电信号计算单个的样本熵;
(3)、多个样本熵点绘制成样本熵曲线,数据长度D=样本熵个数M*窗口大小N;;
(4)、定义阈值,当样本熵达到指定阈值时,记录该点为发作点。
15.根据权利要求14所述的头皮脑电信号回顾性癫痫发作点检测系统,其特征在于:所述脑电信号伪迹去除模块包括去均值化处理模块、白化处理模块、求解最佳分离矩阵模块、获得无伪迹脑电信号模块。
16.根据权利要求14所述的头皮脑电信号回顾性癫痫发作点检测系统,其特征在于:所述信息输出模块包括显示接口输出、音频接口输出、网络接口输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510736832.4A CN105249962B (zh) | 2015-11-03 | 2015-11-03 | 头皮脑电信号回顾性癫痫发作点检测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510736832.4A CN105249962B (zh) | 2015-11-03 | 2015-11-03 | 头皮脑电信号回顾性癫痫发作点检测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105249962A CN105249962A (zh) | 2016-01-20 |
CN105249962B true CN105249962B (zh) | 2019-04-30 |
Family
ID=55090130
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510736832.4A Expired - Fee Related CN105249962B (zh) | 2015-11-03 | 2015-11-03 | 头皮脑电信号回顾性癫痫发作点检测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105249962B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106361324B (zh) * | 2016-10-27 | 2020-07-14 | 河南科技大学 | 基于遗传约束盲源分离的胎儿心电精确提取方法 |
CN108021873B (zh) * | 2017-11-22 | 2022-02-15 | 湖北师范大学 | 一种聚类非对称互信息的脑电信号癫痫分类方法及系统 |
EP3920781A4 (en) * | 2019-02-08 | 2022-10-12 | Nanyang Technological University | METHODS AND SYSTEM FOR SEIZURE DETECTION |
CN111870241B (zh) * | 2020-07-03 | 2022-05-17 | 杭州电子科技大学 | 基于优化的多维样本熵的癫痫发作信号检测方法 |
CN113499083B (zh) * | 2021-05-18 | 2023-10-20 | 生物岛实验室 | 一种头皮脑电高频振荡自动检测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102429657A (zh) * | 2011-09-22 | 2012-05-02 | 上海师范大学 | 一种癫痫脑电信号分类检测装置及方法 |
CN102835955A (zh) * | 2012-09-08 | 2012-12-26 | 北京工业大学 | 一种无需设定阈值的脑电信号中眼电伪迹自动去除的方法 |
CN103690160A (zh) * | 2013-11-18 | 2014-04-02 | 浙江大学 | 一种基于非高斯时序模型的脑电特征提取方法 |
WO2014152565A1 (en) * | 2013-03-15 | 2014-09-25 | Simon Adam J | Wavelet analysis in neuro diagnostics |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6631291B2 (en) * | 2001-05-18 | 2003-10-07 | Instrumentarium Corp. | Closed loop drug administration method and apparatus using EEG complexity for control purposes |
-
2015
- 2015-11-03 CN CN201510736832.4A patent/CN105249962B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102429657A (zh) * | 2011-09-22 | 2012-05-02 | 上海师范大学 | 一种癫痫脑电信号分类检测装置及方法 |
CN102835955A (zh) * | 2012-09-08 | 2012-12-26 | 北京工业大学 | 一种无需设定阈值的脑电信号中眼电伪迹自动去除的方法 |
WO2014152565A1 (en) * | 2013-03-15 | 2014-09-25 | Simon Adam J | Wavelet analysis in neuro diagnostics |
CN103690160A (zh) * | 2013-11-18 | 2014-04-02 | 浙江大学 | 一种基于非高斯时序模型的脑电特征提取方法 |
Non-Patent Citations (1)
Title |
---|
基于遗传算法的多通道癫痫脑电信号盲源分离;沈晋慧 等;《中国医学物理学杂志》;20131130;第30卷(第6期);第4547-4552页 |
Also Published As
Publication number | Publication date |
---|---|
CN105249962A (zh) | 2016-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105249962B (zh) | 头皮脑电信号回顾性癫痫发作点检测方法及系统 | |
CN110811609B (zh) | 基于自适应模板匹配与机器学习算法融合的癫痫棘波智能检测装置 | |
Riaz et al. | EMD-based temporal and spectral features for the classification of EEG signals using supervised learning | |
CN110338786B (zh) | 一种癫痫样放电的识别与分类方法、系统、装置和介质 | |
CN108577834B (zh) | 一种用于癫痫间期棘波自动检测的方法 | |
CN109077715A (zh) | 一种基于单导联的心电信号自动分类方法 | |
Kołodziej et al. | A new method of EEG classification for BCI with feature extraction based on higher order statistics of wavelet components and selection with genetic algorithms | |
CN109497997A (zh) | 基于多数据采集的癫痫发作检测设备及预警系统 | |
US20070225932A1 (en) | Methods, systems and computer program products for extracting paroxysmal events from signal data using multitaper blind signal source separation analysis | |
CN106236080B (zh) | 基于多通道的脑电信号中肌电噪声的消除方法 | |
CN104840186A (zh) | 一种充血性心力衰竭患者自主神经功能的评估方法 | |
CN111000555B (zh) | 一种癫痫脑电信号的训练数据生成方法、自动识别模型建模方法和自动识别方法 | |
CN114145754B (zh) | 基于eeg交叉频率耦合的卒中脑功能评估装置 | |
CN109858537A (zh) | 改进的eemd与csp结合的脑电信号特征提取方法 | |
CN105046273A (zh) | 基于多尺度样本熵的癫痫脑皮层电图信号分类方法 | |
CN108958489A (zh) | 一种基于脑电和眼动仪的图像感兴趣区域快速检出方法 | |
Song et al. | Epileptic seizure detection using brain-rhythmic recurrence biomarkers and onasnet-based transfer learning | |
Botcharova et al. | Power-law distribution of phase-locking intervals does not imply critical interaction | |
CN112336369B (zh) | 一种多通道心音信号的冠心病风险指数评估系统 | |
Prabhakar et al. | EM based non-linear regression and singular value decomposition for epilepsy classification | |
CN113208631A (zh) | 一种基于eeg脑波的眨眼检测方法及系统 | |
CN116491956A (zh) | 一种基于多域特征融合的异常脑电自动检测方法 | |
CN113907770A (zh) | 基于特征融合的棘慢复合波检测识别方法及系统 | |
CN114048784A (zh) | 基于共空间模式模糊小波包的脑电信号疲劳特征提取方法 | |
Hussain et al. | Electroencephalography (EEG) analysis of alcoholic and control subjects using multiscale permutation entropy |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190430 Termination date: 20211103 |