单通道脑电信号中肌电伪迹的自动去除方法技术

本发明专利技术公开了一种用于单通道脑电信号中肌电伪迹的去除方法，将脑电信号通过SSA分解得到P个信号分量；将P个信号分量按行拼接成一个P维数据矩阵；将P维数据矩阵进行时间延迟处理得到若干个数据矩阵；利用MCCA对若干个数据矩阵进行盲源分离，得到源估计矩阵S和混合矩阵A；识别源估计矩阵中与肌电伪迹相关的源；去除源估计矩阵中的肌电伪迹，将识别为肌电伪迹的源置零，得到消除肌电伪迹后的源估计矩阵S’，并通过重构得到去除肌电伪迹后的多通道脑电信号X′＝A*S′；将多通道脑电信号X’的各行求和，即可最终得到去除了肌电伪迹后的单通道脑电信号x’。本发明专利技术去除肌电伪迹的同时尽可能保留了脑电信息，提高了脑电信号分析的准确性。

Automatic removal of EMG artifact in single channel EEG signals

The invention discloses a method for removing the EMG artifacts in single channel EEG signals, which decomposes the EEG signals through SSA to get P signal components, and splice the P signal components into a P dimensional data matrix according to the line. A number of data matrices are obtained by the time delay processing of the P dimension data matrix, and a number of data are obtained by MCCA. The matrix is separated by blind source, the source estimation matrix S and the mixed matrix A are obtained; the source of the EMG artifacts in the source estimation matrix is identified; the EMG artifacts in the source estimation matrix are removed, the source zeros of the EMG artifacts are identified, the source estimation matrix S after the EMG artifact is obtained, and the EMG artifacts are removed after the reconstruction to remove the EMG artifacts. The multi channel EEG X '= A*S', and the multi channel EEG X' can eventually get the single channel electroencephalogram (x \) after the EMG artifacts are removed. The invention removes EMG artifact while retaining EEG information as much as possible, and improves the accuracy of EEG signal analysis.

【技术实现步骤摘要】
单通道脑电信号中肌电伪迹的自动去除方法
本专利技术涉及脑电信号处理
，特别涉及一种单通道脑电信号中肌电伪迹的自动去除方法。
技术介绍
脑电图(Electroencephalograph,EEG)是大脑神经细胞的电生理活动在头皮上的反映。由于具有高时间分辨率、非侵入式、低成本和适合长期监测等特点，脑电图已被广泛用于研究脑功能和病理性脑机制。但是EEG在记录大脑活动的同时，也记录了大脑以外的其他活动产生的电信号，而这些被记录的不是由大脑活动产生的电信号称为伪迹。常见的伪迹包括眼电(Electrooculogram，EOG)、肌电(Electromyography，EMG)、心电(Electrocardiography，ECG)以及由电路或外部设备带来的噪声等。相比于其他类型的伪迹信号，比如心电(ECG)和眼电(EOG)，头部肌肉收缩(即肌电(EMG)信号)对EEG产生的干扰更难以消除，其主要原因可能在于肌电伪迹具有高振幅、宽频谱以及生理分布更广等特点，导致一些简单的伪迹去除方法失效。伪迹的存在对后续脑电信号分析造成了极大的影响，可能会导致重要信息的丢失。因此，伪迹去除是神经信息处理应用中最重要的预处理步骤之一，对神经科学研究以及实际临床诊断都具有重要意义。目前的肌电伪迹去除方法主要分为两类：伪迹拒绝和伪迹校正。伪迹拒绝是将包含伪迹的脑电时段整个去除，这种方法的主要缺点是它也会消除重要的脑电信息，导致数据丢失；而伪迹校正是指采用各种方法去除伪迹成分的同时尽可能的保留脑电信息。伪迹校正技术主要包括：1)滤波，通过使用特定频段的滤波器来去除脑电信号中的肌电伪迹；然而由于肌电信号频谱分布范围广，和脑电信号有较大的重合，所以在去除肌电伪迹的同时也会损失一部分脑电信号。2)小波变换(wavelettransform，WT)，通过小波变换将观测信号分解为一系列基函数，然后对系数进行阈值处理以去除肌电伪迹，最后通过小波重构得到纯净的EEG信号；但是小波变换方法需要大量的实验来选择合适的小波基函数和分解层次。3)盲源分离(BlindSourceSeparation，BSS)，盲源分离类算法用于从多个观测到的混合信号中分析出未知的信号源以及混合矩阵，然后将伪迹相关的信号源去除，最后通过混合矩阵就可以得到纯净信号；常见的BBS算法有：主成分分析(Principalcomponentanalysis，PCA)、独立成分分析(Independentcomponentanalysis，ICA)和典型相关分析(canonicalcorrelationanalysis，CCA)。4)经验模态分解(empiricalmodedecomposition，EMD)，通过EMD将观测到的EEG信号分解得到若干本征模式分量(IMF)，然后将与肌电伪迹相关的IMF去除再重建信号即可得到纯净EEG。EMD是一种数据驱动的方法，用于对非平稳、非线性、随机过程进行处理，因此非常适用于脑电信号分析和处理，但是EMD计算复杂度高，一般需要耗时数分钟，不太适合实时处理。随着近年来可穿戴设备的快速发展以及未来个人医疗的需求，如何有效去除单通道脑电数据中的伪迹信号变得尤为重要。
技术实现思路
本专利技术是为了克服现有技术的不足，提供一种单通道脑电信号中肌电伪迹的自动去除方法，从观测脑电信号中去除肌电伪迹的同时尽可能保留脑电信息，避免脑电信息的丢失，从而提高脑电信号分析的准确性。本专利技术采用如下技术方案：单通道脑电信号中肌电伪迹的自动去除方法，包括以下步骤：步骤一：将通过单通道脑电电极传感器采集得到的脑电信号通过奇异谱分析算法SSA分解得到P个信号分量；步骤二：将步骤一中得到的P个信号分量按行拼接成一个P维数据矩阵；步骤三：将拼接成的P维数据矩阵进行时间延迟处理，得到若干个数据矩阵；步骤四：利用多重集典型相关分析MCCA对步骤三中得到的若干个数据矩阵进行盲源分离，得到源估计矩阵S和混合矩阵A；步骤五：识别源估计矩阵中与肌电伪迹相关的源；步骤六：去除源估计矩阵中的肌电伪迹，将识别为肌电伪迹的源置零，得到消除肌电伪迹后的源估计矩阵S’，并通过重构得到去除肌电伪迹后的多通道脑电信号X′＝A*S′；步骤七：将多通道脑电信号X’的各行求和，即可最终得到去除了肌电伪迹后的单通道脑电信号x’。从以上技术方案可知，本专利技术首先将混合了肌电伪迹的单通道脑电信号经过奇异谱分析(singularspectrumanalysis，SSA)算法分解得到若干信号分量；接着将得到的所有信号分量拼接成一个数据矩阵并进行延迟处理得到若干个数据矩阵；然后利用多重集典型相关分析(MCCA)进行盲源分离，得到若干典型变量；最后，将自相关系数小于选定阈值的分量当成是肌电伪迹，并重构得到去除肌电伪迹后的干净脑电信号。与已有技术相比，本专利技术所取得的有益效果包括：1、使用基于奇异谱分析(SSA)算法和多重集典型相关分析(MCCA)算法的肌电伪迹去除方法，无需参考信号，只需单通道脑电信号即可有效去除肌电伪迹，同时保留大量有用的脑电信号，适用于近年来快速发展的穿戴式单通道脑电设备，并且去除过程中无需人工参与即可自动去除肌电伪迹。2、相比基于EMD的方法，其运行时间更短，具有较佳的实时性，对日常个人脑电分析技术具有重要意义。附图说明图1为本专利技术方法的流程图；图2为半模拟脑电信号图，其中(a)为真实的干净脑电信号，(b)为模拟的肌电伪迹信号，(c)为两者的合成信号；图3为利用MCCA进行分解后得到的估计源信号；图4为原真实的干净脑电信号与采用本专利技术方法去除肌电伪迹后脑电信号的对比图。具体实施方式下面结合附图对本专利技术的具体实施方式进行描述，以便本领域的技术人员更好地理解本专利技术，但本专利技术的实施方式不限于此。需要特别提醒注意的是，在以下的描述中，当已知功能和设计的详细描述也许会淡化本专利技术的主要内容时，这些详细描述在这里将被忽略。本实施例以半模拟脑电信号为例，其流程图参见图1，具体包括以下步骤：步骤一：将通过单通道脑电电极传感器采集得到的脑电信号通过奇异谱分析(SSA)算法分解得到P个信号分量；具体为：设所采集到的单通道含肌电伪迹的脑电信号为x＝[x(1),x(2),...,x(T)]，其中序列x(t)(t＝1,2,…,T)表示脑电信号x在t时刻的值，脑电信号x经SSA分解后产生P个信号分量，即其中每个分量xi的长度也都为T。本步骤一中利用SSA算法将单通道含肌电伪迹的脑电信号分解为P个信号分量的过程为：(1)计算轨迹矩阵(trajectorymatrix)X对于给定的一维时间序列x＝(x1,x2,…,xN)，根据窗口长度L计算其轨迹矩阵X，轨迹矩阵X为Lx(N‐L+1)阶：(2)奇异值分解(SVD)定义矩阵S＝XXT，计算矩阵S的特征值和特征变量，并将特征值按降序排列：λ1,λ2,…,λN(λ1≥λ2≥,…,≥λN≥0)，则其对应的特征向量为u1,u2,…,uL。我们因此可以将轨迹矩阵X表示为d个初等矩阵的和：X＝Ｘ1+X2+…+Xd，其中d＝argmaxi{λi＞0}，(3)分组将初等矩阵Xi的下标{1,2,…,d}分成p个不相交的子集I1,I2,…,Ip，设I＝{i1,i2,…,im}，那么合成矩阵计算集合I1,I2,本文档来自技高网...

【技术保护点】
1.单通道脑电信号中肌电伪迹的自动去除方法，其特征在于，包括以下步骤：步骤一：将通过单通道脑电电极传感器采集得到的脑电信号通过奇异谱分析算法SSA分解得到P个信号分量；步骤二：将步骤一中得到的P个信号分量按行拼接成一个P维数据矩阵；步骤三：将拼接成的P维数据矩阵进行时间延迟处理，得到若干个数据矩阵；步骤四：利用多重集典型相关分析MCCA对步骤三中得到的若干个数据矩阵进行盲源分离，得到源估计矩阵S和混合矩阵A；步骤五：识别源估计矩阵中与肌电伪迹相关的源；步骤六：去除源估计矩阵中的肌电伪迹，将识别为肌电伪迹的源置零，得到消除肌电伪迹后的源估计矩阵S’，并通过重构得到去除肌电伪迹后的多通道脑电信号X′＝A*S′；步骤七：将多通道脑电信号X’的各行求和，即可最终得到去除了肌电伪迹后的单通道脑电信号x’。

【技术特征摘要】
1.单通道脑电信号中肌电伪迹的自动去除方法，其特征在于，包括以下步骤：步骤一：将通过单通道脑电电极传感器采集得到的脑电信号通过奇异谱分析算法SSA分解得到P个信号分量；步骤二：将步骤一中得到的P个信号分量按行拼接成一个P维数据矩阵；步骤三：将拼接成的P维数据矩阵进行时间延迟处理，得到若干个数据矩阵；步骤四：利用多重集典型相关分析MCCA对步骤三中得到的若干个数据矩阵进行盲源分离，得到源估计矩阵S和混合矩阵A；步骤五：识别源估计矩阵中与肌电伪迹相关的源；步骤六：去除源估计矩阵中的肌电伪迹，将识别为肌电伪迹的源置零，得到消除肌电伪迹后的源估计矩阵S’，并通过重构得到去除肌电伪迹后的多通道脑电信号X′＝A*S′；步骤七：将多通道脑电信号X’的各行求和，即可最终得到去除了肌电伪迹后的单通道脑电信号x’。2.根据权利要求1所述的单通道脑电信号中肌电伪迹的自动去除方法，其特征在于，步骤一中利用SSA将脑电信号分解为P个信号分量的过程为：(1)计算轨迹矩阵X对于给定的一维时间序列x＝(x1,x2,…,xN)，根据窗口长度L计算其轨迹矩阵X，轨迹矩阵X为Lx(N‐L+1)阶：(2)奇异值分解定义矩阵S＝XXT，计算矩阵S的特征值和特征变量，并将特征值按降序排列：λ1,λ2,…,λN，其中λ1≥λ2≥,…,≥λN≥0，则其对应的特征向量为u1,u2,…,uL；将轨迹矩阵X表示为d个初等矩阵的和：X＝X1+X2+…+Xd，其中d＝argmaxi{λi＞0}，(3)分组将初等矩阵Xi的下标{1,2,…,d}分成p个不相交的子...

【专利技术属性】
技术研发人员：陈灿，邢晓芬，徐向民，舒琳，
申请(专利权)人：华南理工大学，
类型：发明
国别省市：广东,44