CN111012341B - 基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法 - Google Patents

基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法 Download PDF

Info

Publication number
CN111012341B
CN111012341B CN202010016383.7A CN202010016383A CN111012341B CN 111012341 B CN111012341 B CN 111012341B CN 202010016383 A CN202010016383 A CN 202010016383A CN 111012341 B CN111012341 B CN 111012341B
Authority
CN
China
Prior art keywords
artifact
electroencephalogram
signal
electroencephalogram signal
artifacts
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.)
Active
Application number
CN202010016383.7A
Other languages
English (en)
Other versions
CN111012341A (zh
Inventor
崔兴然
高之琳
顾忠泽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN202010016383.7A priority Critical patent/CN111012341B/zh
Publication of CN111012341A publication Critical patent/CN111012341A/zh
Application granted granted Critical
Publication of CN111012341B publication Critical patent/CN111012341B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/369—Electroencephalography [EEG]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6802—Sensor mounted on worn items
    • A61B5/6803—Head-worn items, e.g. helmets, masks, headphones or goggles
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7253—Details of waveform analysis characterised by using transforms
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Psychology (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明公开了基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法，步骤：通过可穿戴脑电设备采集头皮表层的原始脑电信号；对原始脑电信号中的设备噪声进行滤波；对原始脑电信号中的各种伪迹噪声进行识别，包括运动伪迹、毛刺伪迹、眼动伪迹和肌电伪迹；基于伪迹识别结果，进行自适应伪迹裁剪和信号拼接，得到干净的脑电信号；采用神经网络分类和指标参数对干净脑电信号进行综合质量评估。本发明解决了可穿戴脑电设备由于通道数限制和实时性造成的伪迹去除困难以及脑电信号质量评估方法缺失的问题。

Description

基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法
技术领域
本发明属于脑电信号监测领域，特别涉及了一种脑电信号质量评估方法。
背景技术
脑电信号是由大脑神经元放电产生的，反映大脑活动的微弱生物电信号。由于脑电具有容易获取，非侵入性，以及高时间分辨率的特点，其在科学研究和疾病诊断等方面发挥着越来越重要的作用。随着脑电监测技术的不断发展，市场上出现了一大批可穿戴脑电设备，主要应用于商业，医疗和教育等领域。但是这些可穿戴脑电设备在实际应用过程中面临着很多问题，其中最主要的两大问题是：噪声难以实时去除和缺少信号质量。脑电信号是一种低幅值非稳态信号，在采集过程中受到各种噪声的干扰，比如：低频趋势噪声，高频噪声，工频噪声，运动伪迹，毛刺伪迹，眼动伪迹和肌电伪迹等等。这些干扰的存在严重影响了脑电的基本节律，因此需要提出能够有效识别并去除这些伪迹的算法。由于可穿戴脑电设备在实际应用过程中通常是实时分析的，这就增加了伪迹识别和滤除的难度。同时由于失去了人工筛选和评估的支持，信号质量低的无效脑电信号就无法被提前剔除，这严重影响了脑电信号后续分析的准确性。因此一个实时伪迹自动去除和信号质量评估***可以确保后续分析结果的准确性。
在噪声干扰去除方面，主要目的是在尽可能保留脑电信号的基础上最大程度上的去除干扰。针对低频趋势、工频干扰和高频噪声可以利用巴特沃斯滤波器有效滤除。针对运动，毛刺，眼动和肌电等伪迹，由于其与脑电信号的基本节律存在很多交叠的部分，去伪迹过程中势必会去除部分有用的脑电信号。对于多通道脑电，普遍运用独立成分分析(ICA)盲源分离方法，可以提取出各种伪迹，将伪迹从原始脑电中减去，即可滤除。但是针对目前的便携可穿戴脑电设备，其通道数不足以支持使用ICA方法。目前针对通道数较少的脑电信号的伪迹去除方法主要有快速傅里叶变换法和小波变换法。快速傅里叶变换法和小波变换法由于其本质是基于一个基波分解脑电，这些伪迹的形状和频率各不相同，需要人为的设定基波和滤波阈值，且很难在保留信号的同时有效去除伪迹。这些方法都是不适用于可穿戴脑电采集设备的实时应用的。
通常情况下，经过滤波后的脑电信号即被认为是干净的脑电信号，可用于后续分析，但是实际应用过程中存在一些特殊情况，例如原始信号的信噪比过低，信号幅值低于正常值，脑电信号波段能量占比与正常值不符等等。这些特殊情况的脑电信号质量一般较差，不适合应用于后续分析的，应该被剔除。在科研中，所有脑电数据都会经过研究人员的人为筛选，剔除信号质量差的数据段，但是实际应用过程中无法做到人为筛选，因此迫切需要一个自动评估脑电信号质量的技术。
发明内容
为了解决上述背景技术提到的技术问题，本发明提出了基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法。
为了实现上述技术目的，本发明的技术方案为：
基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法，包括以下步骤：
(1)通过可穿戴脑电设备采集头皮表层的原始脑电信号，原始脑电信号中包含了设备噪声和各种伪迹噪声；
(2)对原始脑电信号中的设备噪声进行滤波；
(3)对原始脑电信号中的各种伪迹噪声进行识别，包括运动伪迹、毛刺伪迹、眼动伪迹和肌电伪迹；
(4)基于步骤(3)的伪迹识别结果，进行自适应伪迹裁剪和信号拼接，得到干净的脑电信号；
(5)采用神经网络分类和指标参数对步骤(4)得到的脑电信号进行综合质量评估。
进一步地，在步骤(2)中，利用巴特沃斯低通滤波器滤出高频噪声；利用巴特沃斯高通滤波器滤出低频噪声；利用巴特沃斯带阻滤波器作为陷波滤出工频干扰。
进一步地，所述巴特沃斯低通滤波器的截止频率小于100Hz，所述巴特沃斯高通滤波器的截止频率小于0.3Hz。
进一步地，在步骤(3)中，所述运动伪迹的识别方法如下：
(1a)利用巴特沃斯带通滤波器提取出0.3-2Hz之间的脑电信号，其中就包含了运动伪迹；
(1b)利用Hilbert方法提取运动伪迹的包络线，包络线能够保留运动伪迹的大致形状并忽略其他信号的影响；
(1c)对上一步获得的包络线做平滑处理，进一步降低其他信号的影响；
(1d)对平滑后的信号做标准分数Z变换，判断超过预设阈值a的数据为运动伪迹，保存其在脑电信号上的索引位置，a为正常干净脑电信号幅值的10倍。
进一步地，在步骤(3)中，所述毛刺伪迹的识别方法如下：
(2a)计算信号的差分和近似导数，识别出有瞬间电位变化的信号；
(2b)采用中值滤波平滑差分和近似导数的绝对值；
(2c)对平滑后的信号做标准分数Z变换，判断超过预设阈值b的数据为毛刺伪迹，保存其在脑电信号上的索引位置。
进一步地，在步骤(3)中，所述眼动伪迹的识别方法如下：
(3a)利用巴特沃斯带通滤波器提取出0.3-10Hz之间的脑电信号，其中就包含了眼动伪迹；
(3b)利用Hilbert方法提取眼动伪迹的包络线，包络线能够保留眼动伪迹的大致形状并忽略其他信号的影响；
(3c)对上一步获得的包络线做平滑处理，进一步降低其他信号的影响；
(3d)对平滑后的信号做标准分数Z变换，判断超过预设阈值c的数据为眼动伪迹，保存其在脑电信号上的索引位置，c为正常干净脑电信号幅值的5倍。
进一步地，在步骤(3)中，所述肌电伪迹的识别方法如下：
(4a)利用巴特沃斯带通滤波器提取出100-120Hz之间的脑电信号，其中就包含了肌电伪迹；
(4b)利用Hilbert方法提取肌电伪迹的包络线，包络线能够保留肌电伪迹的大致形状并忽略其他信号的影响；
(4c)对上一步获得的包络线做平滑处理，进一步降低其他信号的影响；
(4d)对平滑后的信号做标准分数Z变换，判断超过预设阈值d的数据为肌电伪迹，保存其在脑电信号上的索引位置。
进一步地，步骤(4)的具体过程如下：
(401)取可穿戴脑电设备所有通道上识别出的4种伪迹索引位置的并集；
(402)对任意伪迹，寻找该伪迹索引位置前后幅值在[-5,5]μV范围内点的索引位置，并对此取可穿戴脑电设备所有通道上的交集；
(403)在上述交集中找到伪迹前后最近的点，作为伪迹的开始和结束位置，重新定义准确的伪迹位置；
(404)裁剪掉准确的伪迹位置，并将裁剪后剩余的信号拼接在一起，得到干净的脑电信号。
进一步地，在步骤(5)中，计算脑电信号质量评价参数：
(A)信噪比SNRhigh和SNRnotch：
Figure GDA0003547373210000051
Figure GDA0003547373210000052
上式中，S为去除设备噪声的脑电信号，Nhigh为高频噪声，Nnotch为工频干扰；
(B)有效信号率E：
Figure GDA0003547373210000053
上式中，Lleft为伪迹裁剪和信号拼接后脑电信号的数据长度，Loriginal为原始脑电信号的数据长度；
(C)运动伪迹占比率Emove：
Figure GDA0003547373210000054
上式中，Lmove为运动伪迹部分脑电信号的数据长度；
(D)毛刺伪迹占比率Ejump：
Figure GDA0003547373210000055
上式中，Ljump为毛刺伪迹部分脑电信号的数据长度；
(E)眼动伪迹占比率Eeog：
Figure GDA0003547373210000056
上式中，Leog为眼动伪迹部分脑电信号的数据长度；
(F)肌电伪迹占比率Eemg：
Figure GDA0003547373210000057
上式中，Lemg为肌电伪迹部分脑电信号的数据长度；
(G)干净脑电信号的幅值；
(H)Alpha波段功率谱密度占比：利用小波将Alpha波段从干净脑电信号中提取出来，分别计算Alpha波段和干净脑电信号的功率谱密度，最终获得Alpha波段在脑电中的功率谱密度占比；
(I)干净脑电信号的模糊熵。
进一步地，在步骤(5)中，利用小波时频分析方法生成干净脑电信号的时频图，从而将一维信号转换为二维图像，并作为深度神经网络的输入，深度神经网络的输出为输入数据所属于不同类别的概率，类别包括低质量不可用、中等质量可用和高等质量可用。
采用上述技术方案带来的有益效果：
(1)本发明解决了目前可穿戴脑电设备由于通道数限制和实时性造成的伪迹去除困难的问题，实现了伪迹自动识别、裁剪和拼接，该方法在保留脑电信息完整性的同时能够较为彻底的去除设备噪声、运动伪迹、毛刺伪迹、眼动伪迹和肌电伪迹的干扰，相比于其他方法，该方法不受通道数影响并适用于实时***。
(2)本发明解决了目前脑电信号质量评估方法缺失的问题，基于真实脑电信号的特性，总结出信噪比，有效信号率，各伪迹占比率，干净脑电信号的幅值、Alpha波段功率谱密度占比和模糊熵作为信号质量评估的指标参数。同时基于深度神经网络模型，直接判断出该段数据的信号质量等级，决定是否用于后续分析。评估结果可进一步提高后续脑电分析的可靠性、准确性和可解释性。
附图说明
图1是本发明方法的整体流程图；
图2是本发明设备噪声滤波的流程图和过程结果图；
图3是本发明运动伪迹识别的流程图和过程结果图；
图4是本发明毛刺伪迹识别的流程图和两种典型毛刺伪迹的过程结果图；
图5是本发明眼动伪迹识别的流程图和过程结果图；
图6是本发明肌电伪迹识别的流程图和过程结果图；
图7是本发明伪迹裁剪和拼接的示意图；
图8是本发明深度神经网络的结构示意图。
具体实施方式
以下将结合附图，对本发明的技术方案进行详细说明。
如图1所示，本发明主要包括5个步骤：信号采集，设备噪声滤波，四种伪迹(运动伪迹，毛刺伪迹，眼动伪迹，肌电伪迹)识别，伪迹裁剪和信号拼接，脑电信号质量评估。每个步骤的具体操作将在下文进行详细介绍。
步骤1——信号采集：便携可穿戴脑电设备采集头皮表层的原始脑电信号，原始脑电信号中除了脑电信号还包含了设备噪声和各种伪迹噪声。设备噪声是设备本身所造成的噪声干扰，包括：低频噪声，工频干扰和高频噪声。基于这些设备噪声，得到的设备信噪比，可以用于评估该设备本身的信号质量。各种伪迹噪声通常是设备和人体接触不良或者人为造成的一些伪迹干扰，由于这些伪迹的幅值通常远大于脑电信号本身，因此它们的存在会影响分析的准确性，需要进行去除。伪迹噪声包括：被试运动造成的运动伪迹，电极和被试接触不稳定造成的毛刺伪迹，被试眨眼造成的眼动伪迹，肌电伪迹。
步骤2——设备噪声滤波：一般脑电的预处理包括：高通滤波，低通滤波和工频陷波，其目的在于去除设备噪声。脑电信号是头皮层的微弱电信号，幅值通常为50μV，有效频率范围[0.3-100]Hz。如图2中的(a)所示，对原始脑电信号，利用巴特沃斯低通滤波器(截止频率的选取取决于后续分析所关注的脑电频段，但是应当低于头皮脑电的最高有效频率100Hz)滤出的高频噪声；利用巴特沃斯高通滤波器(截止频率的选取一般建议在0.3Hz以下，过高的截止频率将影响脑电信号的完整性)滤出的低频噪声；49-51Hz的巴特沃斯带阻滤波器滤出的50Hz工频干扰(不同地区的工频噪声可能是60Hz)。如图2中的(b)所示得到的滤除设备噪声后的脑电信号，相对原始信号已经没有了设备噪声的干扰，可以看得到脑电本身的信息了。
步骤3——四种伪迹识别：
运动伪迹识别：运动伪迹通常来自于被试的头部运动或者走路等运动，幅值远高于50μV，频率范围通常为[0.3,2]Hz。要消除运动伪迹干扰主要有三种方法：
1、通过高通滤波器滤出低频的运动伪迹的，但是这种方法往往会导致大部分有用的低频脑电信号的丢失；2、根据被试头部的加速度传感器，从脑电信号中分离出脑电中的运动趋势，但是并不是所有脑电设备都同时装有加速度传感器的，且该方法算法尚不成熟，因此该方法的使用受到了很大的限制；3、识别到运动伪迹的地方，直接从脑电中减去该段信号。
本发明能够准确的识别出运动伪迹的所在位置，如图3中的(a)所示具体识别步骤如下：
1.运动伪迹的主要频段一般为0.3-2Hz之间，利用巴特沃斯带通滤波器提取出0.3-2Hz之间的脑电信号，其中就包含了运动伪迹；
2.由于运动伪迹的幅值通常是脑电信号的十几倍甚至几十倍，因此利用Hilbert方法提取其包络线，包络线可以保留运动伪迹的大致形状并忽略其他信号的影响；
3.对上一步获得包络线做平滑处理，进一步降低其他信号的影响；
4.对平滑后的信号做标准分数Z变换，判断超过阈值范围(正常干净脑电信号幅值的十倍)的数据为运动伪迹，保存其在脑电信号上的索引位置；
5.所有通道重复上述步骤。
运动伪迹识别结果如图3中的(b)所示。
毛刺伪迹识别：毛刺伪迹通常是由于电极和皮肤接触不良或者信号传输过程中存在问题，造成电位瞬间变化。为了消除这类伪迹的干扰，首先要识别毛刺伪迹，如图4中的(a)所示采用以下步骤：
1.计算信号的差分和近似导数，可以识别到有瞬间电位变化的信号；
2.中值滤波平滑差分和近似导数的绝对值；
3.对平滑后的信号做标准分数Z变换，判断超过阈值范围的数据为毛刺伪迹，保存其在脑电信号上的索引位置；
4.所有通道重复上述步骤。
毛刺伪迹识别结果如图4中(b)、(c)所示，(b)、(c)分别表示两种典型的毛刺伪迹。
眼动伪迹识别：眼动伪迹是脑电数据中普遍存在的一种伪迹，且出现频繁影响巨大。幅值范围为[50,200]μV，频率范围通常为[0.3,10]Hz。目前常用的眼动伪迹去除方法是盲源分离方法，但是该方法要求足够多的通道数，否则无法有效分离出眼动伪迹信号，因此该方法对于便携脑电设备显然是不适用。而其他类似于基于小波变换或者经验模态分解方法，存在操作复杂和易滤除有用频率段信号的缺点。对于需要应用到实际生活中的便携脑电设备，本发明识别眼动伪迹的方法，如图5中的(a)所示具体步骤如下：
1.利用巴特沃斯带通滤波器提取出0.3-10Hz之间的脑电信号，其中就包含了眼动伪迹；
2.由于眼动伪迹的幅值通常是脑电信号的几倍，因此利用Hilbert方法提取其包络线，包络线可以保留眼动伪迹的大致形状并忽略其他信号的影响；
3.对上一步获得包络线做平滑处理，进一步降低其他信号的影响；
4.对平滑后的信号做标准分数Z变换，判断超过阈值范围(正常干净脑电信号幅值的五倍)的数据为眼动伪迹，保存其在脑电信号上的索引位置；
5.所有通道重复上述步骤。
眼动伪迹识别结果如图5中的(b)所示。
肌电伪迹识别：肌电伪迹普遍存在于脑电信号的采集过程中，其幅值范围和频率范围都非常广泛，很难从脑电中分离和提取出来。但是一般肌电伪迹的幅值要高于50μV且频率较高。虽然上文已经利用低通滤波，滤除了高频噪声，但是由于采用的3阶巴特沃斯滤波器，仍有部分高频信号遗留。基于以上原因，本发明如图6中的(a)所示采用如下步骤识别肌电伪迹：
1.利用巴特沃斯带通滤波器提取出100-120Hz之间的脑电信号，其中就包含了高频肌电伪迹；
2.由于肌电伪迹的幅值通常是脑电信号的几倍，因此利用Hilbert方法提取其包络线，包络线可以强调肌电伪迹的大致位置；
3.对上一步获得包络线做平滑处理，进一步明确肌电伪迹的位置；
4.对平滑后的信号做标准分数Z变换，判断超过阈值范围的数据为肌电伪迹，保存其在脑电信号上的索引位置；
5.所有通道重复上述步骤。
肌电伪迹识别结果如图6中的(b)所示。
步骤4——自适应伪迹裁剪和脑电拼接：通过以上步骤识别出运动、毛刺、眼动和肌电伪迹后，需要将其从脑电信号中裁剪，并将脑电信号重新拼接。未必要裁剪和拼接过程中出现伪迹残留或者拼接毛刺的问题。本发明设计了一种自适应裁剪和拼接方法，有的伪迹在所有通道上都存在，而有的伪迹可能仅在于某一特定通道，但是为确保通道数据长度的一致性，所有通道都要裁剪同样长的伪迹段后进行拼接。因此，如图7所示具体步骤如下：
1.取所有通道上识别到的四种伪迹的索引位置的并集；
2.对任意伪迹，寻找该伪迹索引位置前后幅值在[-5,5]μV范围内点的索引位置，并对此取所有可穿戴脑电设备通道上的交集；；
3.在这一交集中找到伪迹前后最近的点，作为伪迹的开始和结束位置；
4.裁剪掉重新定义的伪迹位置并将裁剪后剩余的信号拼接在一起；
5.对所有伪迹重复步骤2-4，完成伪迹裁剪和脑电拼接。
步骤5——脑电信号质量评估：脑电信号质量的评估主要用于评估该段脑电数据是否可以用于分析，评估方法包括指标参数和神经网络分类。
1.信号质量的指标参数
(A)信噪比SNRhigh和SNRnotch：
Figure GDA0003547373210000111
Figure GDA0003547373210000112
上式中，S为去除设备噪声的脑电信号，Nhigh为高频噪声，Nnotch为工频干扰；
(B)有效信号率E：
Figure GDA0003547373210000113
上式中，Lleft为伪迹裁剪和信号拼接后脑电信号的数据长度，Loriginal为原始脑电信号的数据长度；
(C)运动伪迹占比率Emove：
Figure GDA0003547373210000114
上式中，Lmove为运动伪迹部分脑电信号的数据长度；
(D)毛刺伪迹占比率Ejump：
Figure GDA0003547373210000115
上式中，Ljump为毛刺伪迹部分脑电信号的数据长度；
(E)眼动伪迹占比率Eeog：
Figure GDA0003547373210000116
上式中，Leog为眼动伪迹部分脑电信号的数据长度；
(F)肌电伪迹占比率Eemg：
Figure GDA0003547373210000117
上式中，Lemg为肌电伪迹部分脑电信号的数据长度；
(G)干净脑电信号的幅值：正常的脑电信号其幅值在50μV左右，范围为[20,100]μV。
(H)Alpha波段功率谱密度占比：利用小波将Alpha波段从干净脑电信号中提取出来，分别计算Alpha波段和干净脑电信号的功率谱密度，最终获得Alpha波段在脑电中的功率谱密度占比；
(I)干净脑电信号的模糊熵：利用模糊熵计算经过伪迹去除后的干净脑电信号的熵，如果信号中有噪声残留会导致熵值降低。
以上各指标可以从不同角度详细描述信号的特征。
2.基于深度神经网络的质量分类
利用小波时频分析方法生成伪迹去除后脑电信号的时频图，实现将一维信号转换为二维图像，作为深度神经网络模型的输入。脑电的信号质量被分为三个等级：低质量不可用、中等质量可用和高质量可用，这三个等级作为三个类别用于模型的输出。如图8所示，采用了类似于残差网络体系结构中的连接，通过允许信息在深度神经网络中很好地传播来优化训练。该网络由1个基本块和16个残差块组成。基本块具有两个卷积层(convolutionallayer，简称conv)，残差块具有一个卷积层，并采用最大池化(max pool)。在每个卷积层之后，都添加批量归一化(batch normalization，简称BN)和整流线性单位(rectifiedlinear unit,简称ReLU)，采用预激活块设计。由于采用预激活块结构，网络的第一层和最后一层是特殊情况。我们还在卷积层之间和非线性之后应用Dropout。最终的完全连接层(fully connected layers，简称FC)和归一化(softmax)输出时频图属于各个类别的概率。由于信号质量评估对于数据长度和实验环境等没有特殊要求，因此可以在短时间内收集大量的数据作为训练集，训练集脑电生成的时频图可以用于训练模型。测试集脑电生成的时频图作为模型的输入，得到该段脑电数据属于不同类别的概率。
3.综合评估
上述两个部分分别得到了用于描述原始脑电和去伪迹后脑电的指标参数，以及去伪迹后脑电的信号质量所属类别。信号质量所属类别直接决定该段数据是否用于后续分析，但是深度神经网络的最大缺点在于其不可解释性。当遇到低质量脑电时，无法快速找到造成信号低的原因，而指标参数的辅助则可以很好的弥补这一缺憾。两种方法的结合和互补，在保证其高准确率和实时性的同时能够快速解释影响信号质量高低的原因并及时调整。例如，当检测到连续“低质量不可用”数据段时，查看指标参数，发现运动伪迹成分很大，那么可以提醒用户减少头部运动，提高数据可用性。
实施例仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明保护范围之内。

Claims (9)

1.基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法，其特征在于，包括以下步骤：
(1)通过可穿戴脑电设备采集头皮表层的原始脑电信号，原始脑电信号中包含了设备噪声和各种伪迹噪声；
(2)对原始脑电信号中的设备噪声进行滤波；
(3)对原始脑电信号中的各种伪迹噪声进行识别，包括运动伪迹、毛刺伪迹、眼动伪迹和肌电伪迹；
(4)基于步骤(3)的伪迹识别结果，进行自适应伪迹裁剪和信号拼接，得到干净的脑电信号；具体过程如下：
(401)取可穿戴脑电设备所有通道上识别出的4种伪迹索引位置的并集；
(402)对任意伪迹，寻找该伪迹索引位置前后幅值在[-5,5]μV范围内点的索引位置，并对此取可穿戴脑电设备所有通道上的交集；
(403)在上述交集中找到伪迹前后最近的点，作为伪迹的开始和结束位置，重新定义准确的伪迹位置；
(404)裁剪掉准确的伪迹位置，并将裁剪后剩余的信号拼接在一起，得到干净的脑电信号；
(5)采用神经网络分类和指标参数对步骤(4)得到的脑电信号进行综合质量评估。
2.根据权利要求1所述基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法，其特征在于，在步骤(2)中，利用巴特沃斯低通滤波器滤出高频噪声；利用巴特沃斯高通滤波器滤出低频噪声；利用巴特沃斯带阻滤波器作为陷波滤出工频干扰。
3.根据权利要求2所述基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法，其特征在于，所述巴特沃斯低通滤波器的截止频率小于100Hz，所述巴特沃斯高通滤波器的截止频率小于0.3Hz。
4.根据权利要求1所述基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法，其特征在于，在步骤(3)中，所述运动伪迹的识别方法如下：
(1a)利用巴特沃斯带通滤波器提取出0.3-2Hz之间的脑电信号，其中就包含了运动伪迹；
(1b)利用Hilbert方法提取运动伪迹的包络线，包络线能够保留运动伪迹的大致形状并忽略其他信号的影响；
(1c)对上一步获得的包络线做平滑处理，进一步降低其他信号的影响；
(1d)对平滑后的信号做标准分数Z变换，判断超过预设阈值a的数据为运动伪迹，保存其在脑电信号上的索引位置，a为正常干净脑电信号幅值的10倍。
5.根据权利要求1所述基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法，其特征在于，在步骤(3)中，所述毛刺伪迹的识别方法如下：
(2a)计算信号的差分和近似导数，识别出有瞬间电位变化的信号；
(2b)采用中值滤波平滑差分和近似导数的绝对值；
(2c)对平滑后的信号做标准分数Z变换，判断超过预设阈值b的数据为毛刺伪迹，保存其在脑电信号上的索引位置。
6.根据权利要求1所述基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法，其特征在于，在步骤(3)中，所述眼动伪迹的识别方法如下：
(3a)利用巴特沃斯带通滤波器提取出0.3-10Hz之间的脑电信号，其中就包含了眼动伪迹；
(3b)利用Hilbert方法提取眼动伪迹的包络线，包络线能够保留眼动伪迹的大致形状并忽略其他信号的影响；
(3c)对上一步获得的包络线做平滑处理，进一步降低其他信号的影响；
(3d)对平滑后的信号做标准分数Z变换，判断超过预设阈值c的数据为眼动伪迹，保存其在脑电信号上的索引位置，c为正常干净脑电信号幅值的5倍。
7.根据权利要求1所述基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法，其特征在于，在步骤(3)中，所述肌电伪迹的识别方法如下：
(4a)利用巴特沃斯带通滤波器提取出100-120Hz之间的脑电信号，其中就包含了肌电伪迹；
(4b)利用Hilbert方法提取肌电伪迹的包络线，包络线能够保留肌电伪迹的大致形状并忽略其他信号的影响；
(4c)对上一步获得的包络线做平滑处理，进一步降低其他信号的影响；
(4d)对平滑后的信号做标准分数Z变换，判断超过预设阈值d的数据为肌电伪迹，保存其在脑电信号上的索引位置。
8.根据权利要求1所述基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法，其特征在于，在步骤(5)中，计算脑电信号质量评价参数：
(A)信噪比SNRhigh和SNRnotch：
Figure FDA0003547373200000031
Figure FDA0003547373200000032
上式中，S为去除设备噪声的脑电信号，Nhigh为高频噪声，Nnotch为工频干扰；
(B)有效信号率E：
Figure FDA0003547373200000033
上式中，Lleft为伪迹裁剪和信号拼接后脑电信号的数据长度，Loriginal为原始脑电信号的数据长度；
(C)运动伪迹占比率Emove：
Figure FDA0003547373200000034
上式中，Lmove为运动伪迹部分脑电信号的数据长度；
(D)毛刺伪迹占比率Ejump：
Figure FDA0003547373200000041
上式中，Ljump为毛刺伪迹部分脑电信号的数据长度；
(E)眼动伪迹占比率Eeog：
Figure FDA0003547373200000042
上式中，Leog为眼动伪迹部分脑电信号的数据长度；
(F)肌电伪迹占比率Eemg：
Figure FDA0003547373200000043
上式中，Lemg为肌电伪迹部分脑电信号的数据长度；
(G)干净脑电信号的幅值；
(H)Alpha波段功率谱密度占比：利用小波将Alpha波段从干净脑电信号中提取出来，分别计算Alpha波段和干净脑电信号的功率谱密度，最终获得Alpha波段在脑电中的功率谱密度占比；
(I)干净脑电信号的模糊熵。
9.根据权利要求1所述基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法，其特征在于，在步骤(5)中，利用小波时频分析方法生成干净脑电信号的时频图，从而将一维信号转换为二维图像，并作为深度神经网络的输入，深度神经网络的输出为输入数据所属于不同类别的概率，类别包括低质量不可用、中等质量可用和高等质量可用。
CN202010016383.7A 2020-01-08 2020-01-08 基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法 Active CN111012341B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010016383.7A CN111012341B (zh) 2020-01-08 2020-01-08 基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010016383.7A CN111012341B (zh) 2020-01-08 2020-01-08 基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法

Publications (2)

Publication Number Publication Date
CN111012341A CN111012341A (zh) 2020-04-17
CN111012341B true CN111012341B (zh) 2022-04-22

Family

ID=70198537

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010016383.7A Active CN111012341B (zh) 2020-01-08 2020-01-08 基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法

Country Status (1)

Country Link
CN (1) CN111012341B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112115856B (zh) * 2020-09-17 2024-06-14 浙江迈联医疗科技有限公司 脑电质量评估方法、存储介质及***
CN112043262B (zh) * 2020-09-28 2023-01-31 南京伟思医疗科技股份有限公司 一种新生儿脑电信号伪差的识别***和方法
CN112826510B (zh) * 2021-03-11 2023-02-07 清华大学 电生理信号采集***
CN113040788B (zh) * 2021-06-02 2021-09-21 博睿康科技（常州）股份有限公司 基于频谱分析的脑电信号质量的检测方法
CN113812965B (zh) * 2021-08-19 2024-04-09 杭州回车电子科技有限公司 睡眠状态识别方法、装置、电子装置和存储介质
CN114041803B (zh) * 2021-11-16 2022-10-14 电子科技大学 一种基于区间中值滤波的tms脉冲伪迹的去除装置
CN114189795B (zh) * 2022-02-15 2022-05-17 武汉左点科技有限公司 助听器尖峰噪声自适应调节方法及设备
CN114886388B (zh) * 2022-07-12 2022-11-22 浙江普可医疗科技有限公司 一种麻醉深度监测过程中脑电信号质量的评估方法及装置
CN117064409B (zh) * 2023-10-12 2024-01-30 深圳般意科技有限公司 经颅直流电干预刺激效果实时评估方法、装置及终端
CN117257327A (zh) * 2023-10-24 2023-12-22 沈阳工业大学 一种脑电信号质量评估方法及***
CN117643474B (zh) * 2023-12-24 2024-05-28 沈阳工业大学 基于深度神经网络的脑电信号处理方法及***

Citations (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5603088A (en) * 1995-02-28 1997-02-11 Motorola, Inc. Method and apparatus for determining a quality level of an analog signal in a radio communication system
CA2369163A1 (en) * 2001-01-31 2002-07-31 Markus Naf Automatic colour defect correction
EP1237621A1 (en) * 1999-12-07 2002-09-11 Krasnow Institute Adaptive electric field modulation of neural systems
CN105249961A (zh) * 2015-11-02 2016-01-20 东南大学 基于蓝牙脑电耳机的实时驾驶疲劳检测***及检测方法
CN105361880A (zh) * 2015-11-30 2016-03-02 上海乃欣电子科技有限公司 肌肉运动事件的识别***及其方法
WO2016110804A1 (en) * 2015-01-06 2016-07-14 David Burton Mobile wearable monitoring systems
CN106537901A (zh) * 2014-03-26 2017-03-22 马克·W·帕布利科弗 用于提供定制的娱乐内容的计算机处理方法和***
CN106568589A (zh) * 2016-11-04 2017-04-19 东南大学 一种基于经验小波变换碰摩声发射消噪方法
CN107095672A (zh) * 2017-06-10 2017-08-29 深圳创达云睿智能科技有限公司 一种去除脑电信号中伪迹信号的方法
CN107468243A (zh) * 2017-09-01 2017-12-15 北京大智商医疗器械有限公司 一种评估昏迷程度的方法及装置
EP3323742A1 (en) * 2016-11-22 2018-05-23 Tetra Laval Holdings & Finance S.A. A method of operating an apparatus for applying components to packaging containers
CN108175406A (zh) * 2017-12-14 2018-06-19 中山大学附属第医院 一种治疗中枢性睡眠呼吸暂停症的方法及装置
CN108478218A (zh) * 2018-04-24 2018-09-04 武汉理工大学 基于高性能电路与虚拟仪器的脑电信号采集***
CN108960182A (zh) * 2018-07-19 2018-12-07 大连理工大学 一种基于深度学习的p300事件相关电位分类识别方法
CN109009092A (zh) * 2018-06-15 2018-12-18 东华大学 一种去除脑电信号噪声伪迹的方法
CN109318207A (zh) * 2018-11-07 2019-02-12 西安交通大学 一种利用肌电定时的下肢运动准备电位检测***及方法
CN109497998A (zh) * 2018-12-14 2019-03-22 深圳和而泰数据资源与云技术有限公司 脑电信号质量评估方法、脑电信号采集器和头戴设备
CN110269608A (zh) * 2019-06-10 2019-09-24 浙江普可医疗科技有限公司 一种去除信号干扰的方法、装置以及可读存储介质
CN110353672A (zh) * 2019-07-15 2019-10-22 西安邮电大学 一种脑电信号中眼部伪迹去除***及去除方法
CN110568483A (zh) * 2019-07-22 2019-12-13 中国石油化工股份有限公司 基于卷积神经网络的地震线性噪声压制效果自动评估方法
WO2019239143A1 (en) * 2018-06-15 2019-12-19 CLARK, David Julian An eeg signal processor, and associated methods

Patent Citations (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5603088A (en) * 1995-02-28 1997-02-11 Motorola, Inc. Method and apparatus for determining a quality level of an analog signal in a radio communication system
EP1237621A1 (en) * 1999-12-07 2002-09-11 Krasnow Institute Adaptive electric field modulation of neural systems
CA2369163A1 (en) * 2001-01-31 2002-07-31 Markus Naf Automatic colour defect correction
CN106537901A (zh) * 2014-03-26 2017-03-22 马克·W·帕布利科弗 用于提供定制的娱乐内容的计算机处理方法和***
WO2016110804A1 (en) * 2015-01-06 2016-07-14 David Burton Mobile wearable monitoring systems
CN105249961A (zh) * 2015-11-02 2016-01-20 东南大学 基于蓝牙脑电耳机的实时驾驶疲劳检测***及检测方法
CN105361880A (zh) * 2015-11-30 2016-03-02 上海乃欣电子科技有限公司 肌肉运动事件的识别***及其方法
CN106568589A (zh) * 2016-11-04 2017-04-19 东南大学 一种基于经验小波变换碰摩声发射消噪方法
EP3323742A1 (en) * 2016-11-22 2018-05-23 Tetra Laval Holdings & Finance S.A. A method of operating an apparatus for applying components to packaging containers
CN107095672A (zh) * 2017-06-10 2017-08-29 深圳创达云睿智能科技有限公司 一种去除脑电信号中伪迹信号的方法
CN107468243A (zh) * 2017-09-01 2017-12-15 北京大智商医疗器械有限公司 一种评估昏迷程度的方法及装置
CN108175406A (zh) * 2017-12-14 2018-06-19 中山大学附属第医院 一种治疗中枢性睡眠呼吸暂停症的方法及装置
CN108478218A (zh) * 2018-04-24 2018-09-04 武汉理工大学 基于高性能电路与虚拟仪器的脑电信号采集***
CN109009092A (zh) * 2018-06-15 2018-12-18 东华大学 一种去除脑电信号噪声伪迹的方法
WO2019239143A1 (en) * 2018-06-15 2019-12-19 CLARK, David Julian An eeg signal processor, and associated methods
CN108960182A (zh) * 2018-07-19 2018-12-07 大连理工大学 一种基于深度学习的p300事件相关电位分类识别方法
CN109318207A (zh) * 2018-11-07 2019-02-12 西安交通大学 一种利用肌电定时的下肢运动准备电位检测***及方法
CN109497998A (zh) * 2018-12-14 2019-03-22 深圳和而泰数据资源与云技术有限公司 脑电信号质量评估方法、脑电信号采集器和头戴设备
CN110269608A (zh) * 2019-06-10 2019-09-24 浙江普可医疗科技有限公司 一种去除信号干扰的方法、装置以及可读存储介质
CN110353672A (zh) * 2019-07-15 2019-10-22 西安邮电大学 一种脑电信号中眼部伪迹去除***及去除方法
CN110568483A (zh) * 2019-07-22 2019-12-13 中国石油化工股份有限公司 基于卷积神经网络的地震线性噪声压制效果自动评估方法

Non-Patent Citations (10)

* Cited by examiner, † Cited by third party
Title
Recognition of Emotional States Using Multiscale Information Analysis of High Frequency EEG Oscillations;Gao, Zhilin；Cui, Xingran；Wan, Wang;《ENTROPY》;20190630;第21卷(第6期);全文 *
便携式BCI设备快速自动去眼电伪迹算法的研究;牛群峰,周季冬,王莉,惠延波;《自动化与仪表》;20191231;第34卷(第2期);全文 *
儿童早期脑的发育与医教融合大健康;杨元魁,崔兴然.;《教育家》;20191231(第28期);全文 *
基于伪迹去除的单通道脑机接口技术研究;张璐;《中国优秀硕士学位论文全文数据库 信息科技辑》;20190115;正文第10-11页 *
柔性通信对抗模拟训练***研究与实现;常宏斌.;《万方》;20180830;全文 *
***成瘾戒断群体在MID任务中的奖赏期待与加工机制研究;周昱;《万方》;20191125;全文 *
经颅磁刺激同步干预的头皮脑电信号伪迹离线去除方法综述;银珊,李颖洁;《生物医学工程学杂志》;20191231;第36卷(第1期);全文 *
脑血管病后的康复训练现况和临床应用;崔宝娟;《万方》;20160623;全文 *
针刺对脑缺血后神经元损伤保护作用的研究;许能贵;《万方》;20011231;全文 *
面向转化医学的非线性信号分析及其应用研究;崔兴然.;《万方》;20131129;全文 *

Also Published As

Publication number Publication date
CN111012341A (zh) 2020-04-17

Similar Documents

Publication Publication Date Title
CN111012341B (zh) 基于可穿戴脑电设备的伪迹去除和脑电信号质量评估方法
Hamad et al. Feature extraction of epilepsy EEG using discrete wavelet transform
Islam et al. A wavelet-based artifact reduction from scalp EEG for epileptic seizure detection
Fathima et al. Wavelet based features for epileptic seizure detection
CN105496363A (zh) 基于检测睡眠脑电信号对睡眠阶段进行分类的方法
CN108577834A (zh) 一种用于癫痫间期棘波自动检测的方法
CN114190944B (zh) 基于脑电信号的鲁棒情绪识别方法
Wang et al. Wavelet transform based ECG denoising using adaptive thresholding
CN114757236B (zh) 基于tqwt与svmd的脑电信号去噪优化方法及***
CN113723557A (zh) 一种基于多频带时空卷积网络的抑郁脑电分类***
Shahbakhti et al. SWT-kurtosis based algorithm for elimination of electrical shift and linear trend from EEG signals
Tibdewal et al. Detection and classification of eye blink artifact in electroencephalogram through discrete wavelet transform and neural network
CN114533089A (zh) 一种基于表面肌电信号的下肢动作识别分类方法
CN117918863A (zh) 一种脑电信号实时伪迹处理与特征提取的方法和***
CN116784860B (zh) 基于形态学心拍模板聚类的心电信号特征提取***
Alturki et al. A DWT-band power-SVM based architecture for neurological brain disorders diagnosis using EEG signals
Demirel et al. Single-channel EEG based arousal level estimation using multitaper spectrum estimation at low-power wearable devices
Li et al. Ocular Artifact Removal Algorithm of Single Channel EEG Based on CEEMDAN-ICA-WTD
CN113855047A (zh) 一种基于异常片段识别的fvep信号增强方法
Wei et al. Automatic recognition of chewing noises in epileptic EEG based on period segmentation
CN113208633A (zh) 一种基于eeg脑波的情绪识别方法、***
CN113180663A (zh) 一种基于卷积神经网络的情绪识别方法、***
CN112971786A (zh) 一种基于脑肌电信号小波相干系数的卒中康复评估方法
CN118078305B (zh) 用于可穿戴式血流动力监测仪的动态心电数据处理方法
CN116035594B (zh) 一种基于分割-降噪网络的脑电伪迹去除方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant