买专利,只认龙图腾
首页 专利交易 科技果 科技人才 科技服务 商标交易 会员权益 IP管家助手 需求市场 关于龙图腾
 /  免费注册
到顶部 到底部
清空 搜索

【发明授权】一种面向癫痫患者脑电信号的时频图像分类方法_温州大学_202010419816.3 

申请/专利权人:温州大学

申请日:2020-05-18

公开(公告)日:2021-04-27

公开(公告)号:CN111598003B

主分类号:G06K9/00(20060101)

分类号:G06K9/00(20060101);G06K9/46(20060101);G06T7/45(20170101);G06K9/62(20060101);G06N3/00(20060101)

优先权:

专利状态码:有效-授权

法律状态:2021.04.27#授权;2020.09.22#实质审查的生效;2020.08.28#公开

摘要:本发明公开了一种面向癫痫患者的脑电信号时频图像分类方法,该方法首先对脑电信号进行时频分析得到时频图像用于特征提取,并进行特征选择和分类模型参数优化,从而实现面向癫痫患者的脑电信号时频图像分类。该方法主要包括如下步骤:1时频分析:通过连续小波变换进行频率分解并获得时频图像;2特征提取:结合局部二值模式和灰度共生矩阵描述子提取时频图像的特征。3特征选择与模型参数优化:基于递归消除算法选择有效特征子集以降低计算复杂度,再基于二进制粒子群优化算法同时实现特征子集和参数寻优。本发明方法实现了高精度的癫痫脑电识别,模型泛化性能好。

主权项:1.一种面向癫痫患者的脑电信号时频图像分类方法,其特征在于,包括如下步骤:1时频分析:采用小波分解和连续小波变换的原理将原始脑电信号分解为五个频率子带,并生成整体时频图像和其对应的五个频率子带的时频图;所述步骤1包括如下子步骤:1.1癫痫脑电信号ft的主要频率范围在0~64Hz,采用小波分解的原理将原始脑电信号分解为五个频率子带,分别为δ波、θ波、α波、β波和γ波,其中,δ波的频率范围为0-4Hz,θ波为4-8Hz,α波为8-16Hz,β波为16-32Hz,γ波为32-64Hz;1.2利用步骤1.1得到的五个不同子带的频率范围,选取Morlet作为小波基的连续小波变换方法分别生成对应频率子带的时频图像,再将时频图像按照具有临床意义的脑电信号波段进行有效信息分割获得五个频率子带的时频图;2特征提取:基于步骤1获得的五个频率子带的时频图,分别采用旋转不变的局部二值模式提取局部特征,采用灰度共生矩阵提取全局特征,然后将所述的局部特征和全局特征进行一维向量合并,得到五个频率子带的特征集,进而得到整体时频图像的特征集;所述步骤2包括如下子步骤:2.1对于步骤1.2得到的五个频率子带的时频图,利用局部二值模式进行图像局部纹理特征提取,同时考虑旋转不变性和灰度不变性,根据如下公式得到LBP旋转不变均匀模式和非均匀模式为基于模式得到局部特征向量Li: 其中,r代表邻域半径,P代表在半径r的圆上的邻域点总数,p代表邻域第p个像素点,gi代表图像邻域的灰度值,gc代表图像中心像素的灰度值;ULBPr,p表示以r为半径的圆周上两个邻近值1转换为0或者0转换为1的转换次数;若ULBPr,p≤2,该均匀模式从属于p+1类,反之,全部的非均匀模式从属于1类,整幅图像的纹理直方图的特征向量维数为p+2;2.2对于步骤1.2得到的五个频率子带的时频图,根据灰度共生矩阵通过设置方向θ和距离d两个变量,获得时频图像纹理的数字特征矩阵,并利用两个位置像素之间的联合概率密度Pi,j;d,θ来定义灰度共生矩阵,然后,选择能量、对比度、同质性、相关性四类特征值得到时频图像整体空间的特征表示,进而得到全局特征向量Gi;2.3将步骤2.1的局部特征向量Li和步骤2.2的全局特征向量Gi,进行一维向量合并,得到特征向量Xi,进而得到特征集X=[X1,X2,…,XN],其中,N为脑电信号段总数;3特征选择和模型参数优化:对步骤2获得的整体时频图像特征集选择有效特征子集,并对分类模型进行参数优化,获得最优分类模型;所述步骤3包括如下子步骤:3.1对步骤2.3得到的特征集X,运用十倍交叉验证方法得到十个互不重叠的特征子集,利用递归消除算法选择有效特征集S;3.2对于步骤3.1得到的有效特征集S,基于二进制粒子群优化算法,利用n+2维对整个粒子进行编码得到bool类型的特征向量,改变粒子的位置xij和速度vij寻找全局最优解,根据目标函数计算适度值fit,不断迭代直至达到最大种群数量和最大迭代次数T时,根据内存中的全局最优粒子得到SVM最优参数C,σ和最优特征子集,并通过得到的最优参数C,σ来构建最优分类模型;其中,根据目标函数计算适度值fit的公式如下: 其中,f1表示为训练过程中获得的分类准确率,f2表示为选择的最优特征参数,“1”代表特征参数被选中,“0”表示没有选中;n表示为RFE算法提取有效特征向量的维度;在目标函数f中,α表示为分类准确率的占比,β表示为选择特征子集的权重;4将要识别的脑电信号经过步骤1和2的处理后,输入步骤3得到的分类模型,即可得到准确的分类结果。

全文数据:

权利要求:

百度查询: 温州大学 一种面向癫痫患者脑电信号的时频图像分类方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息,力求客观、公正,但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解,仅供参考使用,不能作为本公司承担任何法律责任的依据或者凭证。