2.845

2023影响因子

(CJCR)

  • 中文核心
  • EI
  • 中国科技核心
  • Scopus
  • CSCD
  • 英国科学文摘

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

单通道脑电信号眼电伪迹去除算法研究

刘志勇 孙金玮 卜宪庚

尹培丽, 王建华, 陈阳泉, 王飞跃. 平行测量:复杂测量系统的一个新型理论框架及案例研究. 自动化学报, 2018, 44(3): 425-433. doi: 10.16383/j.aas.2018.c170192
引用本文: 刘志勇, 孙金玮, 卜宪庚. 单通道脑电信号眼电伪迹去除算法研究. 自动化学报, 2017, 43(10): 1726-1735. doi: 10.16383/j.aas.2017.c160191
YIN Pei-Li, WANG Jian-Hua, CHEN Yang-Quan, WANG Fei-Yue. Parallel Measurements: A New Theory and Framework for Complex Measurement System and A Case Study. ACTA AUTOMATICA SINICA, 2018, 44(3): 425-433. doi: 10.16383/j.aas.2018.c170192
Citation: LIU Zhi-Yong, SUN Jin-Wei, BU Xian-Geng. EOG Artifact Removing Method for Single-channel EEG Signal. ACTA AUTOMATICA SINICA, 2017, 43(10): 1726-1735. doi: 10.16383/j.aas.2017.c160191

单通道脑电信号眼电伪迹去除算法研究

doi: 10.16383/j.aas.2017.c160191
基金项目: 

哈尔滨市科技创新人才研究专项资金 2015RA XXJ038

中央高校基本科研业务费专项资金 2013004

国家自然科学基金 61301012

中央高校基本科研业务费专项资金 2013005

详细信息
    作者简介:

    刘志勇 哈尔滨工业大学电气工程及自动化学院博士研究生.2008年获得哈尔滨工业大学电气工程及自动化学院学士学位.主要研究方向为生物医学信号处理.E-mail:liuzhiyong563@hit.edu.cn

    卜宪庚 哈尔滨医科大学基础医学院教授.主要研究方向为物联网及其应用.E-mail:bxg@ems.hrbmu.edu.cn

    通讯作者:

    孙金玮 哈尔滨工业大学电气工程及自动化学院教授.主要研究方向为生物医学传感器, 主动噪声控制理论.本文通信作者, E-mail:jwsun@hit.edu.cn

EOG Artifact Removing Method for Single-channel EEG Signal

Funds: 

Sci-tech Innovation Foundation of Harbin 2015RA XXJ038

Fundamental Research Funds for the Central Universities 2013004

National Natural Science Foundation of China 61301012

Fundamental Research Funds for the Central Universities 2013005

More Information
    Author Bio:

    Ph. D. candidate at the School of Electrical Engineering and Automation, Harbin Institute of Technology. He received his bachelor degree from Harbin Institute of Technology in 2008. His main research interest is biomedical signal processing

    Professor at the School of Basic Medical, Harbin Medical University. His research interest covers internet of things and its application

    Corresponding author: SUN Jin-Wei  Professor at the School of Electrical Engineering and Automation, Harbin Institute of Technology. His research interest covers biomedical sensors and active noise control theory. Corresponding author of this paper, E-mail:jwsun@hit.edu.cn
  • 摘要: 由眨眼和眼动产生的眼电伪迹(Electrooculography,EOG)信号是脑电信号(Electroencephalography,EEG)中的主要噪声信号之一.目前,多通道脑电信号中眼电伪迹的去除算法已经较为成熟.而在单通道脑电信号的眼电伪迹去除中,由于采集通道数量较少且缺乏参考眼电信号,目前尚无十分有效的去除方法.本文提出一种基于小波变换(Wavelet transform,WT)、集合经验模态分解(Ensemble empirical mode decomposition,EEMD)和独立成分分析(Independent component analysis,ICA)的WT-EEMD-ICA单通道脑电信号眼电伪迹去除算法.实验表明:WT-EEMD-ICA算法有效地解决了单通道WT-ICA算法中的超完备问题,能够有效去除单通道脑电信号中的眼电伪迹,并且分离出的眼电伪迹成分与参考通道采集的眼电信号相关性较强.
  • 机器是改造世界的工具, 仪器是认识世界的工具, 而认识世界往往是改造世界的先导.仪器是测量技术的载体, 测量技术是机械制造的"眼睛".随着机械制造业智能数字化、交叉综合化、柔性集成化以及精密微型化的发展总趋势, 使得测量技术的范围越来越广, 精度要求越来越高, 面临的挑战越来越大[1-2].测量参数的定义更加复杂, 被测对象的尺寸小到纳米、原子量级, 大到几十米甚至几百米, 形貌也从简单型面扩展到复杂型面.计量模式从非现场离线计量发展到现场在线计量, 从传统的静态计量到现在的静态与动态同步计量.计量形态也从简单信息到多信息融合.传统的计量方式和量值溯源方法受到极大挑战.

    传统的量值溯源方式是利用高精度标准计量器具实现测量系统的量值传递.譬如, 齿轮测量系统通采用渐开线样板、螺旋线样板, 来验证相关测量任务的正确性[3], 利用球面样板[4], 楔形样板[5]等实现齿轮齿距、齿形、齿向的测量溯源.然而对于弧齿锥齿轮、环面蜗杆、齿轮刀具等复杂工件, 其高精度实物样板难以制造, 目前工程上仍然采用上述简单样板作为替代物来校准测量系统.这种替代作法, 难以证明测量系统对复杂工件的测量结果是实际工件误差的真实反映, 不能满足复杂工件的测量溯源要求.由于缺少实物标准器, 复杂测量系统的测量结果溯源问题亟待解决.

    针对复杂系统, 王飞跃提出了平行系统的思想, 通过构造实际系统与人工系统并行互动的平行系统, 借助人工系统使复杂问题简单化, 以此实现复杂系统的控制与管理[6].平行系统主要通过整合人工社会, 计算实验和平行执行等方法, 形成新的计算研究体系[7-11].目前, 平行系统已经在经济、能源、交通、机器学习、视觉计算等领域取得了大量的成果, 并不断丰富和完善起来[12-34].我们尝试将平行系统的思想扩展并引入到测量领域建立复杂测量系统的一种新型理论框架, 构建与物理测量系统平行执行的人工测量系统, 将物理测量系统能够溯源至人工测量系统.用人工测量系统验证物理测量系统, 以复杂验证复杂, 而不是简单替代复杂的方式, 解决传统方法无法解决的复杂测量系统的测量溯源问题.

    量值溯源是通过一条具有规定不确定度的不间断的比较链, 使测量结果或测量标准的值能够与规定的参考标准(通常是国家计量基准或国际计量基准)联系起来的特性.溯源的过程也是测量结果不确定度的评定过程.目前, 复杂测量系统的测量不确定度评定主要采用以下三种方法:校准工件法, 灵敏系数分析法(GUM法)和计算机仿真法[35].

    对于一般形位误差的测量, 通常采用量块、量规、标准球等标准器具实现测量结果的溯源.通过采用与被测工件形状相似、尺寸接近的校准工件, 校准测量系统, 从而确定该测量任务的精度.在同一测量条件下, 利用该测量系统对与校准工件相似的被测工件进行测量.因为校准工件已溯源至国际单位, 所以通过校准工件法可以实现被测工件测量结果的溯源[36].该方法要求被测件与校准工件几乎相同, 校准工件的精度要远远高于被测工件和测量系统自身的精度.

    对于无校准工件的测量任务, 可利用高精度仪器对复杂测量系统的各个误差成分进行检测或经验估计, 建立不确定度传递函数的数学表达式, 利用GUM法合成实现测量结果的不确定度评定[37].该方法适用于所有具有明确解析表达式且一阶可微的测量模型, 并且要求尽可能完全地罗列出所有的不确定度来源[38].此时, 溯源依赖于误差成分的检测方法以及不确定度传递函数的正确性.

    针对无法建立解析表达式的复杂测量模型, 研究学者们用计算机仿真法实现测量结果的不确定度评定.该方法是根据输入量与输出量的关系建立测量过程的数学模型, 比如德国联邦物理技术研究院(Physikalisch-Technische Bundesanstalt, PTB)建立的虚拟坐标测量机(Virtual coordinate measuring machine, VCMM)[39]、日本建立的虚拟齿轮检查仪(Virtual gear checker, VGC)[40]、北京工业大学建立的虚拟坐标测量机[41]等.通过对复杂测量系统的各个误差分量进行检测, 经验估计其误差概率分布, 将实际测量点或工件理论模型、真实测量过程中获取的误差分布作为输入; 采用蒙特卡洛法产生特定误差分布的随机数, 通过测量过程的数学模型进行大量重复仿真, 由仿真结果的统计特性评定测量结果的不确定度.该方法实现的溯源依赖于各个误差分量的检测, 人为概率分布的假设以及测量过程数学模型等的准确性.

    复杂测量系统的溯源研究现状说明, 校准工件法, 灵敏系数分析法无法满足对无解析表达式、无校准工件的复杂工件测量任务的溯源要求, 计算机仿真法相对于其他方法具有一定的优势.然而, 复杂测量系统的测量过程是一个复杂的动态过程, 不确定度来源囊括测量系统硬件、软件、被测工件、安装、环境、人员等方面的误差.这些误差非常复杂且相互耦合, 各项误差的确定、人为概率分布假设以及不确定度传递函数的正确性难以评价, 从而使复杂测量系统无校准工件时的测量溯源的正确性难以保证.

    另外, 计算机仿真法是以离线、静态、辅助的形式应用于实际测量系统的评估.该方法涵盖大量测量误差原始数据的获取及统计概率分布的假定, 并通过虚拟系统产生大量误差数据, 虽然有数据驱动的影子, 但其属于确定论或概率方法、需要事先设置、是"程式化"的设计, 并不能满足目前测量系统的实时动态的运行需求.

    若要实现复杂测量系统的溯源, 必须将测量系统及测量过程都溯源至标准国际单位制SI上.我们尝试借助平行系统的思想, 构建与物理测量系统"等价"的人工测量系统, 人工测量系统与物理测量系统实现测量过程的平行执行, 将物理测量系统及其测量过程溯源至人工测量系统的理论模型上, 用人工测量系统的计算实验和平行执行特性验证物理测量系统, 解决传统方法无法解决的复杂测量系统的测量溯源问题.

    ACP (Artificial societies, computational experiments, and parallel execution)的核心就是把复杂系统中"虚"的部分, 分解成可定量、可计算、可执行的过程.针对复杂的物理测量系统, 构造人工测量系统, 使来自物理、社会及信息社会的知识形式化、计算化、可视化, 以在线嵌入和实时反馈的方式实现描述解析、预测解析和诱导解析的功能.目标就是促使实际测量过程趋向人工测量过程, 即主动控制技术.通过计算和比较发现更优化的运行状态, 引导实际系统逼近人工系统, 从而借助人工系统减少实际坐标测量系统相关目标的不确定性, 化多样为归一, 使复杂变简单, 以此实现复杂测量系统的溯源要求.

    ACP步骤:利用人工测量系统(A)对复杂测量问题建模; 利用计算实验(C)对复杂测量任务进行各种情况下的大量模拟、分析和评估, 优化得到有效的解决方案; 将人工测量系统和实际测量系统并举, 通过一定的方式进行虚实互动, 以平行执行(P), 引导和管理物理测量过程.平行测量系统的理论框架如图 1所示.

    图 1  平行测量系统的理论框架
    Fig. 1  Theoretical framework for parallel measurement systems

    平行测量系统是由物理测量系统与人工测量系统构成的一个虚实互动、平行执行的系统.人工系统能反映并引导物理测量系统优化运行, 同时人工系统基于数据驱动, 引导人工系统运行.人工测量系统(虚拟空间)是实际测量系统(物理空间)的伴生系统.它集成了外部物理世界的大数据资源和虚拟人工世界的内部信息系统的数据和模型, 实时跟踪记录信息, 按需描述系统的状态.

    通过认识实际系统各要素间的联系和演化规律, 抽象出数学模型, 建立与实际系统测量特性"等价"的人工系统, 根据实际测量信息数据驱动人工测量系统进行计算实验, 然后对具体的测量方案进行全面、准确、及时的评估和修正.最后将人工系统与实际系统相结合, 通过平行执行实施优化控制和管理.平行测量系统的基本架构如图 2所示.

    图 2  平行测量系统的基本架构图
    Fig. 2  Basic architecture of parallel measurement systems

    1) 人工测量系统

    人工测量系统是物理测量系统在虚拟空间的映射, 是通过对物理测量系统的探索与研究, 抽象出来的一个或多个模型组成的系统.一个完整的物理测量系统主要包括四大要素:测量设备、测量方法、测量人员和被测对象.对上述要素进行抽象, 采用数据驱动和语义建模, 构建信息和行为之间的反馈.人工测量系统的架构主要由基础架构层, 数据、知识层(A), 计算实验层(C)和平行执行层(P)四层组成.

    建立人工测量系统的前提是对实际测量系统的运行数据进行获取、存储和建模.建模初期, 人工测量系统在时间上落后于实际测量系统的运行, 处于被动地位; 通过不断的数据更新、模型完善, 追赶实际测量系统, 直至与实际测量系统对应起来; 最后处于主导地位, 使实际测量系统逼近人工测量系统, 并以人工测量系统来引导实际测量系统的运行.

    由于复杂测量系统"不可分不可知"[10]的本质, 因此不存在"一劳永逸"的解决方案, 必须要引入一个会学习且可自适应变化的人工测量系统, 不断地学习处在变化和发展之中的物理测量系统.学习人工的经验知识, "感知"环境和测量系统的动态变化, 反映测量系统特性的变化, 存储与分析历史大数据, 形成能够平行地反映物理测量系统状态发展和变化的人工测量系统.

    人工测量系统的构建采用计算机图形学、虚拟现实、机器学习、知识自动化等技术.其测量场景可借鉴"平行视觉"[16]的技术方法, 被测对象的三维重构采用光学图像采集和数据处理[42-43]、高精度重构技术[44]等, 人员的建模可基于代理控制的智能体"平行学习"[13]技术, 误差数据库的建立采用激光干涉仪等静态误差测量仪器或装置[45-46]、动态误差采集装置[47-48]、温度湿度传感器等完成误差数据的采集, 通过蒙特卡洛方法模拟随机误差的分布[49-50].

    2) 计算实验

    评价一个测量系统的特性主要考察其稳定性、偏倚、线性、重复性和再现性, 且评价是在统计条件下进行[51].利用人工测量系统计算实验的可设计性以及可反复进行的特点, 我们可以对实际测量系统进行各种关于性能可靠性、方案鲁棒性、测量结果可靠性及正确性等同一条件下短时间内的大量反复测量实验、"极限"实验等, 而这些"实验"一般无法用于实际测量系统, 特别是正在运行的实际测量系统, 或无条件, 或经济成本巨大.

    利用人工测量系统的计算实验优化测量策略.研究某一测量策略对测量精度的影响时, 完全没有必要建立针对某一实际测量系统实际测量任务的精确测量模型, 因为测量系统模型相当复杂, 理论上一个测量系统的测量结果难以直接推广到另一测量系统, 而且一个合适的测量模型, 尽管并未代表任何一个实际测量系统的测量过程, 但它却能描述某一"可能测量系统"的测量过程.因此, 一项合理有效的测量策略, 必须在一个合理的测量系统模型中产生良好的效果, 否则该项策略就需改进.

    利用人工测量系统实现实时动态的误差评估.人工测量系统中误差数据需要从实际中采集小数据, 然后根据误差的特性, 利用蒙特卡洛法或神经网络法产生误差的大数据, 产生与物理测量系统相同的误差特性, 并不断和物理测量系统的运行状态进行对比, 修正人工测量系统的差异, 从而实现人工测量系统和物理测量系统在误差特性上的表现的一致性.最后, 根据人工测量系统的误差仿真对实际测量系统进行误差分离、误差补偿和测量不确定度评定.这比目前的仿真系统多了修正和反馈环节, 真正实现物理测量系统的实时动态化模拟运行.

    3) 平行执行

    复杂测量系统属于工程复杂性较高的系统, 我们采用平行控制方法, 并以控制为主.

    物理测量系统与人工测量系统平行执行.通过物理测量系统与人工测量系统的相互连接, 借助误差反馈实时修正人工测量系统, 进行离线和在线的调整, 使二者测量特性等价, 达到人工测量系统能够验证物理测量系统的目的.平行执行过程主要有两部分:离线状态更新和在线优化及评估.

    a) 离线状态更新:人工测量系统执行虚拟测量时必须根据实际情况对许多信息进行更新, 其中包括测量系统中测量设备部件的磨损、维修和更换情况, 测量人员知识水平和技巧能力, 被测工件的表面数据信息, 测量方案等.

    b) 在线优化及评估:人工测量系统试图尽可能地模拟物理测量系统, 同一水平的测量人员在相同的环境设置下采用相同的测量策略基于代理控制分别同时执行实际测量和虚拟测量, 对二者的测量过程和测量结果进行对比和分析, 通过观察实际测量系统与人工测量系统评估的状态之间的不同, 产生误差反馈信号, 对人工测量系统的评估方式或参数进行修正, 减少差别, 并开始分析新一轮的优化和评估.优化和评估的极限即可达到二者测量特性的完全一致, 也可根据一致性评价标准界定终止条件.

    总之, 人工测量系统平行执行过程以计算实验为基础, 一方面通过接收物理测量系统感知的环境、测量设备、工件等各类信息, 对人工测量系统进行离线和在线的状态更新, 实现对物理测量系统的在线监控和评价; 另一方面可根据对物理测量系统感知信息的处理结果分析, 实现在线测量策略的制定、选择、实施及效果评估, 并通过误差反馈实现二者测量特性的一致, 最终使得物理测量系统能够溯源至人工测量系统.

    齿轮在位测量系统是典型的复杂测量系统.随着工业技术的发展, 齿轮测量在测量原理方面实现了从"比较测量"到"啮合运动测量"再到"模型化测量"; 在测量技术上, 经历了"以机械为主"到"机电结合"的过程后, 发展为现在的"光-机-电"与"信息技术"的综合集成[52]; 测量功能也远不限于齿轮, 还扩展到齿轮刀具、环面蜗杆、弧齿锥齿轮等复杂工件.被测工件的多样化, 仪器设备的复杂化, 信息的融合化, 测量的动态化等使得齿轮在位系统发展具有显著的复杂系统的特征.因此, 必须借助人工测量系统进行齿轮测量系统的设计、控制、评估和管理.

    人工齿轮测量系统与物理齿轮测量系统是平行系统的两大组成部分.人工齿轮测量系统是实际齿轮测量系统在虚拟空间的映射, 包括齿轮测量设备中的控制系统硬件、软件和机械系统的仿真, 各种误差源的数据, 工件数据等.平行齿轮测量系统如图 3所示.

    图 3  平行齿轮测量系统
    Fig. 3  Parallel gear measurement system

    图 3图 2的实例化说明.人工测量系统测量过程的仿真要求和真实测量过程具有一致的功能和特性.人工测量系统能够仿真真实齿轮测量系统的测量全过程.首先采用计算机图形学、虚拟现实技术、平行视觉技术对物理齿轮测量系统进行功能抽象, 构造人工齿轮测量系统的几何模块、运动模块、误差模块和虚拟工件模块, 并建立各模块的联系[53-57].利用传感器技术、高精度测量技术对物理齿轮测量系统的几何误差、运动误差、测头误差等误差数据进行采集, 存储至云端.通过误差大数据发现误差规律, 再利用蒙特卡洛法产生大量的误差数据, 用于测量不确定度的评定.采用光学图像采集与处理技术对被测工件表面特征数据进行获取, 利用高精度重构技术对工件模型进行重建, 确保被测工件的一致性.根据真实测量过程中的实时误差采集数据对人工测量系统进行修正, 使二者拥有一致的误差特性.

    人工齿轮测量系统与物理齿轮测量系统采用同一测量软件系统, 即实际测量软件既能驱动物理齿轮测量系统, 也能不做修改地驱动人工齿轮测量系统.人工齿轮测量系统在数据驱动与测量软件控制下与物理测量系统平行执行, 根据计算实验产生大量数据, 根据数据挖掘技术, 产生有效决策, 优化物理齿轮测量系统.

    1) 测量方案优化应用

    当有新的测量任务到达时, 人工齿轮测量系统首先模拟多个齿轮测量方案, 通过分析模拟运行结果, 将最优的测量方案发送给物理齿轮测量系统, 物理齿轮测量系统在测量齿轮的过程中, 将运行的状况实时反馈给人工齿轮测量系统, 人工齿轮测量系统以物理齿轮测量系统的实时状态为初始条件, 模拟运行测量方案, 再将模拟运行结果发送给物理齿轮测量系统, 物理齿轮测量系统根据模拟运行结果实时调整测量参数, 实现物理齿轮测量系统的测量方案优化.

    2) 软件误差溯源应用

    实际齿轮测量软件可同时驱动物理测量系统和人工齿轮测量系统, 根据这一特性通过人工齿轮测量系统的仿真运行和计算实验来测试和分析齿轮测量软件的误差.通过构建误差已知的数字化工件[58], 导入至人工齿轮测量系统中进行仿真测量, 根据仿真测量结果与已知的工件误差, 检验齿轮测量软件的误差.将检验结果反馈回物理齿轮测量系统, 指导物理齿轮测量系统测量软件的修正.

    3) 复杂工件的测量溯源应用

    对于无校准工件的复杂工件, 利用人工齿轮测量系统和数字化工件的仿真测量实验指导复杂工件的校准工件设计与研制, 并在物理齿轮测量系统中进行测量, 使物理齿轮测量系统逼近人工齿轮测量系统.譬如, 弧齿锥齿轮、环面蜗杆、齿轮刀具等难以制造高精度校准工件的复杂工件, 寻求利用简单曲面如球面、柱面、锥面和平面等组合的方式代替复杂曲面, 实现复杂曲面的溯源.并建立该替代曲面的数字化模型及测量策略, 在人工齿轮测量系统中进行仿真测量.并根据测量结果优化工件模型和测量策略, 从而指导复杂工件的校准工件设计、研制与测量, 最终实现复杂工件的测量溯源.

    本文分析了制造领域内测量技术随工业技术的发展现状和趋势, 总结了复杂测量系统的溯源研究现状, 指出现有的研究无法满足目前复杂测量系统的动态化测量溯源需求, 并指出了人工测量系统的必要性和传统仿真系统的不同之处.基于ACP方法提出了平行测量系统的理论框架, 研究了平行测量系统的学习培训、实验评估和管理控制的工作过程.最后将平行测量系统的理论应用于齿轮在位测量系统案例, 实现齿轮测量的优化、控制和溯源.

    平行测量系统的提出, 为制造业的智能数字化、交叉综合化、柔性集成化等提供理论技术基础.未来可提出平行制造系统, 并将平行制造系统与平行测量系统进行有效的无缝连接, 实现制造与测量的真正一体化, 让测量更好地为制造服务, 使制造质量实时有标准可依, 使产品决策更加经济高效.

    本文首次提出平行测量系统, 其应用研究仍处于探索阶段, 一些理论和技术仍待进一步的研究和深入.然而, 毫无疑问, 平行测量是迈向智能测量的一条有效途径.


  • 本文责任编委 田捷
  • 图  1  WT-EEMD-ICA算法框图

    Fig.  1  The block of WT-EEMD-ICA algorithm

    图  2  实验设计框图

    Fig.  2  The block of experiment design

    图  3  EEG$_{\rm ac}$和EOG$_{\rm ac}$波形

    Fig.  3  Wave of EEG$_{\rm ac}$和EOG$_{\rm ac}$

    图  4  $P1$段EEG$_{\rm ac}$和EOG$_{\rm ac}$波形

    Fig.  4  Wave of EEG$_{\rm ac}$和EOG$_{\rm ac}$ in $P1$

    图  5  带通滤波前后EEG与EOG信号对比

    Fig.  5  EEG and EOG processed by band pass filter

    图  6  带阻滤波前后EEG与EOG信号对比

    Fig.  6  EEG and EOG processed by band stop filter

    图  7  EEG信号小波分解图

    Fig.  7  Wavelet coefficients of EEG

    图  8  $d4$小波系数EEMD结果图

    Fig.  8  The EEMD results of $d4$

    图  9  $d4$小波系数ICA结果图

    Fig.  9  The ICA results of $d4$

    图  10  $d4$小波系数WT-EEMD-ICA结果

    Fig.  10  The WT-EEMD-ICA results of $d4$

    图  11  $d3$小波系数WT-EEMD-ICA结果

    Fig.  11  The WT-EEMD-ICA results of $d3$

    图  12  $P1$段信号WT-EEMD-ICA结果

    Fig.  12  The WT-EEMD-ICA results of signals in $P1$

    图  13  $P2$段信号WT-EEMD-ICA结果

    Fig.  13  The WT-EEMD-ICA results of signals in $P2$

    图  14  $P1+P2$段信号WT-EEMD-ICA结果

    Fig.  14  The WT-EEMD-ICA results of signals in $P1+P2$

    图  15  WT-EEMD-ICA算法稳定性实验结果

    Fig.  15  Stability results of WT-EEMD-ICA

    图  16  $P1$段信号WT算法结果

    Fig.  16  The WT algorithm results of signals in $P1$

    图  17  $P1$段信号WT-ICA算法结果

    Fig.  17  The WT-ICA algorithm results of signals in $P1$

    图  18  $P1$段信号EEMD-ICA算法结果

    Fig.  18  The EEMD-ICA algorithm results of signals in $P1$

    图  19  $P1$段信号DE-ICA算法结果

    Fig.  19  The DE-ICA algorithm results of signals in $P1$

    图  20  $P1$段信号HHT算法结果

    Fig.  20  The HHT algorithm results of signals in $P1$

    表  1  脑电分量与眼电分量近似熵对比

    Table  1  Comparison of $ApEn$ between EEG and EOG

    $ApEn$脑电分量眼电分量
    均值1.350.55
    标准差0.260.23
    下载: 导出CSV

    表  2  重构信号与原始信号相关系数

    Table  2  The correlation coefficients between reconstructed and original signal

    对比算法 $S1$ $S2$EOG
    WT-EEMDICA0.820.740.54
    WT0.620.670.56
    WT-ICA0.540.590.57
    EEMD-ICA0.700.590.31
    HHT0.640.570.32
    下载: 导出CSV
  • [1] Marker R J, Maluf K S. Effects of electrocardiography contamination and comparison of ECG removal methods on upper trapezius electromyography recordings. Journal of Electromyography and Kinesiology, 2014, 24(6):902-909 doi: 10.1016/j.jelekin.2014.08.003
    [2] 胡进, 侯增广, 陈翼雄, 张峰, 王卫群.下肢康复机器人及其交互控制方法.自动化学报, 2014, 40(11):2377-2390 http://www.aas.net.cn/CN/abstract/abstract18514.shtml

    Hu Jin, Hou Zeng-Guang, Chen Yi-Xiong, Zhang Feng, Wang Wei-Qun. Lower limb rehabilitation robots and interactive control methods. Acta Automatica Sinica, 2014, 40(11):2377-2390 http://www.aas.net.cn/CN/abstract/abstract18514.shtml
    [3] Zeng H, Song A G, Yan R Q, Qin H Y. EOG artifact correction from EEG recording using stationary subspace analysis and empirical mode decomposition. Sensors, 2013, 13(11):14839-14859 doi: 10.3390/s131114839
    [4] 王行愚, 金晶, 张宇, 王蓓.脑控:基于脑-机接口的人机融合控制.自动化学报, 2013, 39(3):208-221 http://www.aas.net.cn/CN/abstract/abstract17800.shtml

    Wang Xing-Yu, Jin Jing, Zhang Yu, Wang Bei. Brain control:human-computer integration control based on braincomputer interface. Acta Automatica Sinica, 2013, 39(3):208-221 http://www.aas.net.cn/CN/abstract/abstract17800.shtml
    [5] Jie D, Tao W, Zhang A T, Dai H Y. The removal of blink and saccade artifact in EEG recordings by independent component analysis. In:Proceedings of the 3rd International Conference on Biomedical Engineering and Informatics. Yantai, China:IEEE, 2010, 3:1071-1075 http://ieeexplore.ieee.org/xpl/tocresult.jsp?sortType%3Dasc_p_Sequence%26filter%3DAND%28p_IS_Number%3A5639235%29&pageNumber=3
    [6] Daly I, Scherer R, Billinger M, Müller-Putz G. FORCe:fully online and automated artifact removal for brain-computer interfacing. IEEE Transactions on Neural Systems and Rehabilitation Engineering, 2015, 23(5):725-736 doi: 10.1109/TNSRE.2014.2346621
    [7] Mannan M M N, Jeong M Y, Kamran M A. Hybrid ICAregression:automatic identification and removal of ocular artifacts from electroencephalographic signals. Frontiers in Human Neuroscience, 2016, 10:193-193 https://core.ac.uk/display/43879890
    [8] Qi W. EOG artifacts removal in EEG measurements for affective interaction with brain computer interface. In:Proceedings of the 8th International Conference on Intelligent Information Hiding and Multimedia Signal Processing. Piraeus-Athens, Greece:IEEE, 2012. 471-475
    [9] Chang W D, Cha H S, Kim K, Im C H. Detection of eye blink artifacts from single prefrontal channel electroencephalogram. Computer Methods and Programs in Biomedicine, 2016, 124:19-30 doi: 10.1016/j.cmpb.2015.10.011
    [10] Vigon L, Saatchi M R, Mayhew J E W, Fernandes R. Quantitative evaluation of techniques for ocular artefact filtering of EEG waveforms. IEE Proceedings-Science Measurement and Technology, 2000, 147(5):219-229 doi: 10.1049/ip-smt:20000475
    [11] Fatourechi M, Bashashatia A, Ward R K, Birch G E. EMG and EOG artifacts in brain computer interface systems:a survey. Clinical Neurophysiology, 2007, 118(3):480-494 doi: 10.1016/j.clinph.2006.10.019
    [12] Gomez-Herrero G, Clercq W D, Anwar H, Kara O, Egiazarian K, van Huffel S, van Paesschen W. Automatic removal of ocular artifacts in the EEG without an EOG reference channel. In:Proceedings of the 7th Nordic Signal Processing Symposium. Reykjavik, Iceland:IEEE, 2006. 130-133
    [13] 陈永强, 王宏霞.基于新型混合模型的欠定盲分离方法.自动化学报, 2014, 40(7):1412-1420 http://www.aas.net.cn/CN/abstract/abstract18412.shtml

    Chen Yong-Qiang, Wang Hong-Xia. A method for underdetermined blind source separation based on new mixture model. Acta Automatica Sinica, 2014, 40(7):1412-1420 http://www.aas.net.cn/CN/abstract/abstract18412.shtml
    [14] Jung T P, Makeig S, Westerfield M, Townsend J, Courchesne E, Sejnowski T J. Removal of eye activity artifacts from visual event-related potentials in normal and clinical subjects. Clinical Neurophysiology, 2000, 111(10):1745-1758 doi: 10.1016/S1388-2457(00)00386-2
    [15] del R Millan J, Mourino J, Franze M, Cincotti F, Varsta M, Heikkonen J, Babiloni F. A local neural classifier for the recognition of EEG patterns associated to mental tasks. IEEE Transactions on Neural Networks, 2002, 13(3):678-686 doi: 10.1109/TNN.2002.1000132
    [16] Croft R J, Barry R J. Removal of ocular artifact from the EEG:a review. Neurophysiologie Clinique, 2000, 30(1):5-19 doi: 10.1016/S0987-7053(00)00055-1
    [17] Zhou W D, Gotman J. Removing eye-movement artifacts from the EEG during the intracarotid amobarbital procedure. Epilepsia, 2005, 46(3):409-414 doi: 10.1111/epi.2005.46.issue-3
    [18] Sarossy M G, Lee M H A, Bach M. A fast automated method for calculating the EOG arden ratio. Documenta Ophthalmologica, 2014, 128(3):169-178 doi: 10.1007/s10633-014-9430-5
    [19] 吴明权, 李海峰, 马琳.单通道脑电信号中眼电干扰的自动分离方法.电子与信息学报, 2015, 37(2):367-372 doi: 10.11999/JEIT140602

    Wu Ming-Quan, Li Hai-Feng, Ma Lin. Automatic electrooculogram separation method for single channel electroencephalogram signals. Journal of Electronics and Information Technology, 2015, 37(2):367-372 doi: 10.11999/JEIT140602
    [20] James C J, Lowe D. Single channel analysis of electromagnetic brain signals through ICA in a dynamical systems framework. In:Proceedings of the 23rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society. Istanbul, Turkey:IEEE, 2001. 1974-1977
    [21] Davies M E, James C J. Source separation using single channel ICA. Signal Processing, 2007, 87(8):1819-1832 doi: 10.1016/j.sigpro.2007.01.011
    [22] Mammone N, La Foresta F, Morabito F C. Automatic artifact rejection from multichannel scalp EEG by wavelet ICA. IEEE Sensors Journal, 2012, 12(3):533-542 doi: 10.1109/JSEN.2011.2115236
    [23] Wu Z H, Huang N E. Ensemble empirical mode decomposition:a noise-assisted data analysis method. Advances in Adaptive Data Analysis, 2009, 1(1):1-41 doi: 10.1142/S1793536909000047
    [24] Tang Y, Li J P, Wu H. A simple and accurate ICA algorithm for separating mixtures of up to four independent components. Acta Automatica Sinica, 2011, 37(7):794-799 doi: 10.1016/S1874-1029(11)60210-3
    [25] Li H Q, Feng X L, Cao L, Li E B, Liang H A, Chen X L. A new ECG signal classification based on WPD and ApEn feature extraction. Circuits Systems and Processing, 2016, 35(1):339-352 doi: 10.1007/s00034-015-0068-7
    [26] Li M A, Yang L B, Yang J F. Separation of EOG artifacts from EEG signals using Hilbert-Huang transform. In:Proceedings of the 2011 International Conference on Electric Information and Control Engineering. Wuhan, China:IEEE, 2011. 4453-4456
  • 期刊类型引用(11)

    1. 王东庆,周建华,伏云发. EWT算法在单通道脑电信号眼电伪迹自动去除中的研究. 电子测量与仪器学报. 2023(02): 121-129 . 百度学术
    2. 李沛洋,高晓辉,朱鹏程,黄伟杰,李存波,司亚静,徐鹏,田银. 一种基于局部密度的自适应眼电伪迹去除方法. 电子与信息学报. 2022(02): 464-476 . 百度学术
    3. 彭双凌,郭泉成,陆超. 基于短时能量谱的眼电伪迹去除算法研究. 科技与创新. 2022(06): 57-60 . 百度学术
    4. 姚健康,熊根良. 小波包分析和FastICA相结合对单通道脑电信号的去噪研究. 现代电子技术. 2021(07): 60-65 . 百度学术
    5. 张丹,张帅,李小俚,康健楠. 基于小波相干性算法的孤独症儿童脑电评估. 自动化学报. 2021(03): 678-684 . 本站查看
    6. 熊馨,杨秋红,周建华,徐保磊,李永程,尹旭贤,伏云发. 脑机融合控制中脑电伪迹处理方法. 昆明理工大学学报(自然科学版). 2021(03): 56-70 . 百度学术
    7. 李昕,安占周,李秋月,蔡二娟,王欣. 基于多重多尺度熵的孤独症静息态脑电信号分析. 自动化学报. 2020(06): 1255-1263 . 本站查看
    8. 吴全玉,张文强,潘玲佼,陶为戈,刘晓杰. 一种结合自适应噪声完备经验模态分解和盲反卷积去除脑电中眼电伪迹的新方法. 数据采集与处理. 2020(04): 720-729 . 百度学术
    9. 孙晓娟,王利. 基于CEEMD小波包算法的降噪方法研究. 计算机与现代化. 2020(09): 73-76+82 . 百度学术
    10. 牛群峰,周季冬,王莉,惠延波. 便携式BCI设备快速自动去眼电伪迹算法的研究. 自动化与仪表. 2019(02): 103-108 . 百度学术
    11. 马煜. 多通道网络传输数据断点区域检测方法研究. 电子设计工程. 2019(13): 79-83 . 百度学术

    其他类型引用(28)

  • 加载中
  • 图(20) / 表(2)
    计量
    • 文章访问数:  2396
    • HTML全文浏览量:  582
    • PDF下载量:  956
    • 被引次数: 39
    出版历程
    • 收稿日期:  2016-03-09
    • 录用日期:  2016-07-18
    • 刊出日期:  2017-10-20

    目录

    /

    返回文章
    返回