一种心电伪迹在线去除算法制造技术

本发明专利技术涉及核磁环境下脑电信号的心电伪迹在线去除算法。包括：一、基于可变精度Teager能量算子的R峰实时检测算法，二、基于R峰检测的实时平均伪迹去除算法，三、多通道参考的实时自适应噪声消除算法。该方法主要用于在线去除功能磁共振与脑电信号同步采集时所带来的心电伪迹，为功能磁共振与脑电信号联合分析以及实时神经反馈提供了高质量的信号。该算法可以对心电伪迹进行较好的去除，运算速度快。本发明专利技术可以广泛应用在实时功能磁共振与脑电信号的临床与基础研究中，并可产生巨大的社会效益和经济效益。

【技术实现步骤摘要】

本专利技术提供，属于神经信息处理
，具体是指 一种对核磁环境下同步采集的脑电信号中引入的心电伪迹进行在线去除的算法。
技术介绍
实时功能磁共振与脑电联合采集（rtfMRI/EEG)是近些年兴起的一种多模态神经 成像技术，其中脑电信号（EEG)直接反映神经元的电活动，其时间分辨率可以达到毫秒级 另IJ，但空间分辨率较低，而功能磁共振成像（fMRI)的空间分辨率较高，可以达到毫米量级， 但时间分辨率较低。EEG与fMRI在时空分辨率上的互补性，使得rtfMRI/EEG能够为脑功能 研宄提供更为深入和全面的信息。然而，目前传统的EEG/fMRI对采集到的数据还仅限于对 离线数据进行分析，而rtfMRI/EEG则可以对采集到的数据进行在线处理，因此可以在线观 测大脑活动状态或者进行实时神经反馈，以达到改善大脑特定认知功能的目的。rtfMRI/EEG首先面临的是伪迹去除问题，其中EEG采集系统向fMRI信号引入的 伪迹较小，不会显著影响fMRI成像质量；而核磁扫描系统则会向EEG信号中引入严重的伪 迹，其中最为显著的是梯度伪迹和心电伪迹。梯度伪迹由联合采集时磁共振扫描仪不断切 换的梯度磁场产生，具有比较规律的周期重复模式，因而比较容易去除。心电伪迹则是由心 跳引起的微小头动、头部血管随心跳而产生的扩张和收缩等多种因素导致，它与心跳节律 有着明显的相关性，无论功能像的扫描是否进行，心电伪迹都会存在。此外，心电伪迹的形 态因被试、脑电通道和磁共振扫描仪的不同而变化，甚至在同一试次同一通道内，随着时间 的行进，心电伪迹的形态也会发生变化。因此，这些因素的存在使得心电伪迹成为联合采集 中最难以去除的伪迹。目前，心电伪迹去除算法主要分成两类：波形去除法和空间模式去除法。前者按照 一定的方法构造伪迹模板，然后将模板从EEG信号中减除；后者则一般利用独立成分分析 或者主成分分析等方法将EEG信号分解，然后剔除其中的伪迹成分，最后再将信号重构还 原。然而，现有的心电伪迹去除算法方法仍以离线方法为主，无法满足rtfMRI/EEG中对心 电伪迹进行实时去除的要求。
技术实现思路
为此，本专利技术提出了，它可以对核磁环境下实时获取 的EEG信号进行心电伪迹在线去除，在保证伪迹去除效果的前提下，满足rtfMRI/EEG系统 的实时要求。 ，包括： 基于可变精度Teager能量算子的R峰实时检测算法，对心电通道中的R峰位置进 行实时检测； 基于R峰检测的实时平均伪迹去除算法，对脑电通道的心电伪迹进行初步的降噪 处理； 多通道参考的实时自适应噪声消除算法，利用多通道组合参考通道进一步消除实 时平均伪迹去除算法在EEG信号中残余的心电伪迹。 所述的基于可变精度Teager能量算子的R峰实时检测算法，用于实时检测心电 (ECG)信号中R峰的位置。实时传输的ECG信号通过滑动窗限定分析的数据长度（一般为 IOs左右，最优长度可以通过实验来确定），新数据会不断地替换掉窗内已分析的数据。首 先将已去除梯度伪迹的ECG信号送入如下的可变精度Teager能量算子中： ^rk{x(n)}=x(n)2+x(n-k)x(n+k) (I) 通过选择合适的k值，可变精度Teager能量算子可以对ECG信号中的R峰进行增 强，同时抑制其他部分。然后，以当前滑动窗长内最大值的2/3为阈值，在可变精度Teager 能量算子的输出波形中得到R峰的疑似位置。最后通过与R峰模板的相关性校验来确定R 峰的最终位置。 所述的基于R峰检测的实时平均伪迹去除算法，用于初步消除EEG通道中的心电 伪迹。按照前面R峰实时检测算法的结果，将当前滑动窗内的EEG信号进行分段，每一个分 段的中心在R峰之后的210ms，其长度为中线前后各0. 3s，分段的长度由下式计算得出： 其中，L为模板长度，以点数计，n为参加计算的R峰的个数，Ri为第i个R峰间距； 除最新的一个分段外，将其余所有的分段进行平均叠加，计算得到一个心电伪迹模板；从最 新的分段内减去上一步得到的心电伪迹模板。采用实时平均伪迹去除算法后EEG信号中仍 然会存在残余的心电伪迹， 所述的多通道参考的实时自适应噪声消除算法，用于消除EEG通道中的残余心电 伪迹。多通道参考的自适应噪声消除算法（multi-channelreferencedAdaptiveNoise Cancellation，mANC)采用多个心电相关的脑电通道作为参考通道，并组合为复合参考通 道，以便更加准确地捕捉伪迹的动态变化，从而提高心电伪迹的去除效果。参考通道选取的 基本原则是：1)远离脑电的参考电极，2)左右半球对称分布。以四个参考通道为例，组合算 法如下： 其中是用于四个通道进行组合的系数向量，每一个系数对应于 一个参考通道，其具体值由最小二乘算法实时计算得到，以使得P与当前正在处理的EEG通 道之间的差异最小。Ii,j代表第i个参考通道在滑动窗内的第j个点。式（3)中每一行有 k个值，代表滑动窗长内的k个点。p为四个参考通道组合在一起之后的复合参考通道，同 时作为自适应滤波器的参考信号。自适应滤波器中EEG原信号的估计式如下： Si=Xi-WrP^w(0) =0⑷ 其中Xi是正在处理的EEG数据点，w是自适应滤波器的系数矩阵，p是上一步中计 算得出的复合参考通道，元为去除残余心电伪迹之后的EEG数据点。W(O)为w的初始值，在 运算过程中w等自适应滤波器参数都是实时更新的。 本专利技术的优势在于提供了一种快速有效去除心电伪迹的算法，通过可变精度 Teager能量算子和多通道参考的实时自适应噪声消除算法，提高了R峰的检测效果和心电 伪迹的去除效果。相较于传统的心电伪迹去除算法，本专利技术能够更好地估测出EEG通道中 残余的心电伪迹的形状，从而可以达到比传统方法更好的去噪效果。【附图说明】 图1心电伪迹在线去除算法流程图 图2不同k值下可变精度Teager能量算子的输出结果示意图 图3心电伪迹在线去除算法效果图【具体实施方式】 按照图1所示的算法流程图分步骤并结合实验数据来介绍本方法的实施过程。 步骤一、ECG信号中R峰的实时检测： EEG信号由电极帽上的专用电极采集而得，在对ECG信号进行R峰实时检测之前， 与EEG信号处理，首先对ECG信号进行梯度伪迹去除；然后使用可变精度Teager能量算子 处理ECG信号，由于k值决定了算法对于待检信号频率的敏感度，且该值与待检信号的频率 和实验的硬件环境都有关系，需要通过预实验来确定，图2是一次预实验中随着k的变化， 可变精度Teager能量算子输出的变化，随着k的增大，被抑制的其他部分也逐渐增强；之后 利用阈值得到R峰疑似位置，通过与R峰模板的相关性校验确定R峰的最终位置；R峰的初 始模板通过如下方式获得：在启动核磁设备之前先启动EEG采集系统，EEG设备可以记录到 一段不含梯度伪迹的ECG信号；以这一段ECG信号为基础，R峰实时检测算法可以通过平均 叠加的方式计算出一个0. 5s长度的R峰模板；随着滑动窗的移动，R峰模板不断更新。 步骤二、实时平均伪迹去除算法对EEG通道进行初步处理： 根据步骤一检测的R峰最终位置，对滑动窗内的EEG信号进行分段处理，将这些分 段按照R峰的位置对齐，然后叠加平本文档来自技高网...

【技术保护点】
一种心电伪迹在线去除算法，其特征在于，利用可变精度Teager能量算子方法，对心电通道中的R峰位置进行实时检测，并根据R峰检测结果，采用实时平均伪迹去除算法和多通道参考的实时自适应噪声消除算法，对脑电通道中的心电伪迹进行降噪处理，从而完成心电伪迹的去除。

【技术特征摘要】

【专利技术属性】
技术研发人员：赵小杰，姚力，康明轩，
申请(专利权)人：北京师范大学，
类型：发明
国别省市：北京;11