《睡眠状态分析中去除眼电伪迹的设备.pdf》由会员分享，可在线阅读，更多相关《睡眠状态分析中去除眼电伪迹的设备.pdf（11页珍藏版）》请在专利查询网上搜索。
1、(19)中华人民共和国国家知识产权局 (12)发明专利 (10)授权公告号 (45)授权公告日 (21)申请号 201610840431.8 (22)申请日 2016.09.21 (65)同一申请的已公布的文献号 申请公布号 CN 106236083 A (43)申请公布日 2016.12.21 (73)专利权人 广州视源电子科技股份有限公司 地址 510530 广东省广州市广州黄埔区云 埔四路6号 (72)发明人 赵巍 胡静 韩志 (74)专利代理机构 广州华进联合专利商标代理 有限公司 44224 代理人 潘桂生 (51)Int.Cl. A61B 5/0476(2006.01) A61B 。
2、5/00(2006.01) (56)对比文件 CN 105342605 A,2016.02.24,全文. CN 104688220 A,2015.06.10,全文. CN 101869477 A,2010.10.27,全文. US 2011/0257556 A1,2011.10.20,全文. CN 103720471 A,2014.04.16,全文. 审查员 李慧 (54)发明名称 睡眠状态分析中去除眼电伪迹的设备 (57)摘要 本发明涉及一种睡眠状态分析中去除眼电 伪迹的设备， 包括： 脑电电极、 眼电电极、 参考电 极及其连接的模数转换器， 以及通过模数转换器 和滤波电路连接的处理器； 脑。
3、电电极用于检测原 始脑电信号； 眼电电极用于采集眼电信号； 模数 转换器用于模数转换， 滤波电路用于低频滤波后 输入至处理器； 处理器， 用于对每帧滤波后的脑 电信号进行经验模态分解， 将其分解成若干个本 征模函数， 计算各个本征模函数与同一时刻的眼 电信号之间的相关系数； 查找并删除相关系数大 于预设阈值的本征模函数和相关系数最大的本 征模函数， 利用剩余的本征模函数重建每帧脑电 信号。 本发明可以减少去除眼电伪迹过程对脑电 信号的波形的影响， 保留了原始信号的大部分细 节信息。 权利要求书2页 说明书5页 附图3页 CN 106236083 B 2018.02.16 CN 10623608。
4、3 B 1.一种睡眠状态分析中去除眼电伪迹的设备， 其特征在于， 包括： 脑电电极、 眼电电极、 参考电极、 模数转换器、 滤波电路以及处理器； 所述脑电电极、 眼电电极、 参考电极分别连接模数转换器， 并依次通过所述模数转换器 和滤波电路连接至处理器； 所述脑电电极用于检测用户在睡眠中的脑电信号； 所述眼电电极用于采集用户在睡眠 中的眼电信号； 所述模数转换器将眼电信号和脑电信号转换为数字信号， 所述滤波电路对眼电信号和 脑电信号进行低频滤波后输入至处理器； 所述处理器， 用于对每帧滤波后的脑电信号进行经验模态分解， 将其分解成若干个本 征模函数， 计算各个本征模函数与同一时刻的眼电信号之间。
5、的相关系数； 查找并删除相关 系数大于预设阈值的本征模函数和相关系数最大的本征模函数， 利用剩余的本征模函数重 建每帧脑电信号。 2.根据权利要求1所述的睡眠状态分析中去除眼电伪迹的设备， 其特征在于， 所述脑电 电极设置在用户的额头位置； 所述参考电极设置在用户的耳垂； 所述眼电电极设置在眼角 位置； 所述滤波电路输出0-256Hz频段的信号。 3.根据权利要求1所述的睡眠状态分析中去除眼电伪迹的设备， 其特征在于， 所述处理 器， 用于根据设定帧长度采集用户的眼电信号和脑电信号； 对该帧脑电信号进行经验模态 分解， 将其分解成若干个本征模函数， 得到本征模函数集合； 分别计算所述本征模函数。
6、集合 的各个本征模函数与同一时刻的眼电信号之间的相关系数； 查找出相关系数大于预设阈值 的本征模函数和相关系数最大的本征模函数， 并将其从本征模函数集合中删除； 利用本征 模函数集合中剩余的本征模函数重建该帧脑电信号。 4.根据权利要求3所述的睡眠状态分析中去除眼电伪迹的设备， 其特征在于， 所述本征 模函数集合包括如下公式： 式中， EEGoriginal表示脑电信号， imfi表示第i个本征模函数， Re表示残差函数。 5.根据权利要求1所述的睡眠状态分析中去除眼电伪迹的设备， 其特征在于， 所述处理 器采用如下公式重建脑电信号： 式中， EEGpure表示重建的脑电信号， corrcoe。
7、f表示相关系数， imfi表示第i个本征模函 数， EOG表示眼电信号， corrcoefmax表示最大的相关系数， thre表示预设的相关系数阈值。 6.根据权利要求3所述的睡眠状态分析中去除眼电伪迹的设备， 其特征在于， 所述处理 器， 还用于在相关系数大于预设阈值的本征模函数中， 计算本征模函数与眼电信号的欧氏 距离， 将欧氏距离最小的本征模函数从本征模函数集合中剔除。 7.根据权利要求3所述的睡眠状态分析中去除眼电伪迹的设备， 其特征在于， 所述处理 器， 还用于在相关系数大于预设阈值的本征模函数中， 计算本征模函数与眼电信号的余弦 距离， 将余弦距离最小的本征模函数从本征模函数集合中。
8、剔除。 8.根据权利要求3所述的睡眠状态分析中去除眼电伪迹的设备， 其特征在于， 所述处理 权 利 要 求 书 1/2 页 2 CN 106236083 B 2 器， 还用于在对脑电信号进行经验模态分解前， 将脑电信号帧划分为若干个时间窗口， 并对 每个时间窗口的脑电信号进行本征模函数分解； 以及在重建脑电信号后将各个时间窗口重 建的脑电信号进行合并， 得到脑电信号帧。 9.根据权利要求3所述的睡眠状态分析中去除眼电伪迹的设备， 其特征在 于， 所述处 理器， 还用于在重建脑电信号时， 根据本征模函数的排列顺序， 选择本征模函数集合中位置 靠前的若干个本征模函数进行重建脑电信号。 10.根据权。
9、利要求8所述的睡眠状态分析中去除眼电伪迹的设备， 其特征在于， 所述预 设阈值为0.5， 所述帧长度为30s， 所述时间窗口长度为5s或10s。 权 利 要 求 书 2/2 页 3 CN 106236083 B 3 睡眠状态分析中去除眼电伪迹的设备 技术领域 0001 本发明涉及辅助睡眠技术领域， 特别是涉及一种睡眠状态分析中去除眼电伪迹的 设备。 背景技术 0002 在睡眠中， 人体进行了自我放松及恢复的过程， 因此良好的睡眠是保持身体健康 的一项基本条件； 但是由于工作压力大、 生活作息不规律等原因， 导致了部分人群的睡眠质 量欠佳， 表现为失眠、 半夜惊醒等。 0003 目前市面上已经有。
10、一些设备来帮助人们入睡， 提高睡眠质量。 例如在某一特定睡 眠状态下通过声音、 光信号等人工干预， 避免在熟睡状态下叫醒用户等。 对于辅助睡眠的设 备而言， 为了真正达到提高用户睡眠质量的目的， 正确的识别用户的睡眠状态是非常重要 的。 0004 多导睡眠图(Polysomnography,PSG)， 又称睡眠脑电图， 是目前临床上用于睡眠诊 断和分析的 “金标准” 。 多导睡眠图利用多种生命体征例如脑电、 肌电(颌下)、 眼电、 呼吸、 血 氧等对睡眠进行分析。 在这些体征信号中， 脑电图(electroencephalogram， EOG)处于核心 地位。 脑电图是利用精密的电子仪器， 在。
11、头皮上将来自大脑皮层产生的电活动加以记录并 放大的波形信号。 由于脑电图的信号非常微弱(微伏级)， 容易被来自其他部位的生物电信 号干扰。 当眼电信号幅值较低时(即没有较强烈眼球/眼睑活动如眨眼等)， 眼电信号对脑电 信号的干扰比较微弱。 而眼电信号幅值较高时， 由于眼电信号的频率比正常脑电信号低， 高 幅值的眼电信号叠加在脑电信号上就形成了一个类似于基线漂移的现象。 0005 为了降低眼电信号所带来的影响， 目前有很多去除眼电伪迹的方法。 独立成分分 析(Indepdent component analysis,ICA)是一种常用的方法。 它首先假设输入信号都是统 计独立的非高斯的信号的线性。
12、组合， 然后利用线性变换将来自于信号分离。 它的缺点是(1) 输入信号的假设条件在实际使用中并不能完全满足； (2)对于分离后的多个信号， 还需要进 一步判断哪些信号是 “纯净的” 脑电信号， 哪些信号是被分离出的眼电信号。 此外， 还有方法 假设了一个眼电信号对脑电信号的影响因子(如0.2)， 然后利用脑电信号减去乘以影响因 子的眼电信号的方法去除眼电伪迹， 如公式： EEGpureEEGoriginal-0.2*EOG， 由于存在个体差 异及眼电电极的位置的不同， 一个固定的影响因子并不能很好的适应不同的个体。 0006 此外， 由于在睡眠状态分析中， 脑电信号的波形是一个很重要的睡眠状态。
13、指标。 例 如纺锤波和K复合波的出现， 表示进入了非眼快动睡眠的S2期。 经过传统方法处理后的脑电 信号的波形往往会发生变化， 影响了后续对脑电信号的分析效果。 发明内容 0007 基于此， 有必要针对上述问题， 提供一种睡眠状态分析中去除眼电伪迹的设备， 减 少去除眼电伪迹过程对脑电信号的波形的影响， 确保后续对脑电信号的分析效果。 0008 一种睡眠状态分析中去除眼电伪迹的设备， 包括： 脑电电极、 眼电电极、 参考电极、 说 明 书 1/5 页 4 CN 106236083 B 4 模数转换器、 滤波电路以及处理器； 0009 所述脑电电极、 眼电电极、 参考电极分别连接模数转换器， 并。
14、依次通过所述模数转 换器和滤波电路连接至处理器； 0010 所述脑电电极用于检测用户在睡眠中的脑电信号； 所述眼电电极用于采集用户在 睡眠中的眼电信号； 0011 所述模数转换器将眼电信号和脑电信号转换为数字信号， 所述滤波电路对眼电信 号和脑电信号进行低频滤波后输入至处理器； 0012 所述处理器， 用于对每帧滤波后的脑电信号进行经验模态分解， 将其分解成若干 个本征模函数， 计算各个本征模函数与同一时刻的眼电信号之间的相关系数； 查找并删除 相关系数大于预设阈值的本征模函数和相关系数最大的本征模函数， 利用剩余的本征模函 数重建每帧脑电信号。 0013 上述睡眠状态分析中去除眼电伪迹的设备。
15、， 利用眼电电极采集的用户的眼电信号 和脑电电极采集的脑电信号， 经过模数转换和滤波处理后， 由处理器对每帧滤波后的脑电 信号进行经验模态分解， 将其分解成若干个本征模函数， 计算各个本征模函数与同一时刻 的眼电信号之间的相关系数； 查找并删除相关系数大于预设阈值的本征模函数和相关系数 最大的本征模函数， 利用剩余的本征模函数重建每帧脑电信号。 该设备可以减少去除眼电 伪迹过程对脑电信号的波形的影响， 保留了信号的大部分细节信息， 确保后续对脑电信号 的分析效果。 附图说明 0014 图1为一个实施例的睡眠状态分析中去除眼电伪迹的设备的结构示意图图； 0015 图2是处理器去除眼电伪迹的算法流。
16、程图； 0016 图3是去除眼电伪迹的实验数据结果示意图。 具体实施方式 0017 下面结合附图阐述本发明的睡眠状态分析中去除眼电伪迹的设备的实施例。 0018 参考图1所示， 图1为本发明的睡眠状态分析中去除眼电伪迹的设备的结构示意 图， 包括： 脑电电极、 眼电电极、 参考电极、 模数转换器、 滤波电路以及处理器； 0019 所述脑电电极、 眼电电极、 参考电极分别连接模数转换器， 并依次通过所述模数转 换器和滤波电路连接至处理器； 0020 所述脑电电极用于检测用户在睡眠中的脑电信号； 所述眼电电极用于采集用户在 睡眠中的眼电信号； 0021 所述模数转换器将眼电信号和脑电信号转换为数字。
17、信号， 所述滤波电路对眼电信 号和脑电信号进行低频滤波后输入至处理器； 0022 所述处理器， 用于对每帧滤波后的脑电信号进行经验模态分解， 将其分解成若干 个本征模函数， 计算各个本征模函数与同一时刻的眼电信号之间的相关系数； 查找并删除 相关系数大于预设阈值的本征模函数和相关系数最大的本征模函数， 利用剩余的本征模函 数重建每帧脑电信号。 0023 上述实施例的睡眠状态分析中去除眼电伪迹的设备， 利用眼电电极采集的用户的 说 明 书 2/5 页 5 CN 106236083 B 5 眼电信号和脑电电极采集的脑电信号， 经过模数转换和滤波处理后， 由处理器对每帧滤波 后的脑电信号进行经验模态。
18、分解， 将其分解成若干个本征模函数， 计算各个本征模函数与 同一时刻的眼电信号之间的相关系数； 查找并删除相关系数大于预设阈值的本征模函数和 相关系数最大的本征模函数， 利用剩余的本征模函数重建每帧脑电信号。 该设备可以减少 去除眼电伪迹过程对脑电信号的波形的影响， 保留了信号的大部分细节信息， 确保后续对 脑电信号的分析效果。 0024 后续可以利用该设备输出的脑电信号进行睡眠状态监测和分析等， 当然， 该后续 的处理也可以在所述处理器上去实现。 0025 在一个实施例中， 所述脑电电极设置在用户的额头位置； 所述参考电极设置在用 户的耳垂； 所述眼电电极设置在眼角位置； 如图1所示， 图中。
19、， 脑电电极即图中的 “M” ， 眼电电 极包括左右两个电极， 即图中的 “ROC” 和 “LOC” ， 参考电极设置在用户的耳垂， 即图中 “R” 和 “L” ， 加速度传感器即图中 “AT” 。 滤波电路主要是进行低通滤波和滤除工频干扰， 为了适应 于脑电信号和眼电信号的处理， 滤波电路滤波后， 输出0-256Hz频段的信号至处理器。 0026 对于去除眼电伪迹功能， 主要通过处理器来进行， 基于处理器实现的功能， 可以在 处理器中配置相应的算法模块。 0027 处理器去除眼电伪迹的算法流程包括(1)(5)， 具体如下： 0028 (1)处理器控制眼电电极和脑电电极根据设定帧长度采集用户的。
20、眼电信号和脑电 信号； 0029 如在对用户进行辅助睡眠等睡眠状态分析中， 处理器可以以设定帧长度， 通过用 户佩戴的眼电电极和脑电电极， 采集用户在睡眠过程中产生的眼电信号和脑电信号。 在采 集信号时， 可以以30s为一帧进行采集， 后续对每帧眼电信号和脑电信号进行分析处理。 0030 (2)对该帧脑电信号进行经验模态分解， 将其分解成若干个本征模函数， 得到本征 模函数集合； 0031 在此， 处理器对脑电信号进行经验模态分解， 将其分解成若干个本征模函数 (Intrinsic Mode Function， IMF)和残差函数(Redisual， Re)之和的形式。 0032 本征模函数集。
21、合包括如下公式： 0033 0034 式中， EEGoriginal表示脑电信号， imfi表示第i个本征模函数， Re表示残差函数。 0035 (3)分别计算所述本征模函数集合的各个本征模函数与同一时刻的眼电信号之间 的相关系数； 0036 参考图2， 图2是处理器去除眼电伪迹的算法流程图， 脑电信号进行经验模态分解 后， 得到本征模函数集合， 分别计算本征模函数1-n(imf1imfn)与眼电信号EOG的相关系数 1-n(corrcoef1corrcoefn)。 0037 (4)查找出相关系数大于预设阈值的本征模函数和相关系数最大的本征模函数， 并将其从本征模函数集合中删除； 0038 如。
22、图2所示， 通过设定阈值， 在计算完相关系数后， 将相关系数大于阈值的本征模 函数和相关系数最大的本征模函数删除， 剩下的m个本征模函数。 0039 作为一个实施例， 在计算完相关系数后， 处理器还可以用于在相关系数大于预设 说 明 书 3/5 页 6 CN 106236083 B 6 阈值的本征模函数中； 计算本征模函数与眼电信号的欧氏距离； 从欧氏距离最小的本征模 函数从本征模函数集合中剔除。 0040 类似的， 在计算完相关系数后， 处理器还可以用于在相关系数大于预设阈值的本 征模函数中； 计算本征模函数与眼电信号的余弦距离； 从余弦距离最小的本征模函数从本 征模函数集合中剔除。 004。
23、1 通过上述实施例， 在相关系数判断基础上结合了欧氏距离或余弦距离判断， 可以 将与相关系数判断中无法去除的更多遗留的眼电伪迹去除。 0042 (5)利用本征模函数集合中剩余的本征模函数重建该帧脑电信号； 0043 处理器在去除了眼电伪迹后， 利用剩下的m个本征模函数重建去除了眼电伪迹后 的脑电信号。 作为一个实施例， 处理器在重建脑电信号时， 根据本征模函数的排列顺序， 选 择本征模函数集合中位置靠前的若干个本征模函数进行重建脑电信号。 0044 该实施例中， 由于本征模函数的排列顺序是按频率由大到小， 并且与眼电信号相 似度最高的本征模函数一般排列在中间位置， 因此在重建脑电信号时， 可以。
24、仅利用前若干 个本征模函数， 删除包括相似度最高的本征模函数在内的频率较低的本征模函数后， 再重 建脑电信号。 0045 其中， 处理器可以采用如下公式重建脑电信号： 0046 0047 式中， EEGpure表示重建的脑电信号， corrcoef表示相关系数， imfi表示第i个本征模 函数， EOG表示眼电信号， corrcoefmax表示最大的相关系数， thre表示预设的相关系数阈值。 0048 在一个实施例中， 处理器可以在对脑电信号进行经验模态分解前， 先将脑电信号 帧划分为N个时间窗口， 并对每个时间窗口的脑电信号进行本征模函数分解； 以及在重建脑 电信号后， 将各个时间窗口重建。
25、的脑电信号进行合并， 得到脑电信号帧。 0049 上述实施例， 通过将采集的脑电信号帧划分多个时间窗口并行处理， 能够加快信 号处理速度， 提高睡眠状态分析的效率。 0050 例如， 以30s一帧为例， 可以以5s或10s为一个时间窗口长度。 0051 本发明的睡眠状态分析中去除眼电伪迹的设备， 只去除高幅度眼电所造成的类似 于基线漂移的伪迹， 并且保留了信号的大部分细节信息； 后续使用该脑电信号去进行睡眠 状态分析时， 可以得到更好的效果。 0052 参考图3所示， 图3是去除眼电伪迹的实验数据结果示意图。 图3(a)为采集的脑电 信号， 图3(b)为采集的眼电信号， 图3(c)比较了去除眼。
26、电伪迹前后的脑电信号(图中， 为 脑电信号； 为去除眼电伪迹后的脑电信号)， 下图是上图截取部分放大图， 可以发现， 上述 区间内的数据点之间， 由眼电带来的幅度较大的深V形波动被本方案给消除的同时， 并且保 留了较多的信息。 0053 相对于传统方法(如ICA等)将输入的多路信号视为经过线性组合后的多路源信 号， 并试图将这些信号彼此分离， 能在周期信号上传统方法能获得比较好的效果。 而且对于 脑电信号而言， 由于脑电信号和眼电信号都可以视为随机信号， 并且脑电信号容易受到外 部干扰， 很难将脑电信号和眼电信号彻底分离开， 此时脑电信号就会混入额外的噪声信号， 加大了后续信号处理分析的难度。。
27、 而本发明的技术方案， 只去除高幅度眼电带来的类似于 说 明 书 4/5 页 7 CN 106236083 B 7 基线漂移的伪迹， 保留了信号的大部分细节信息。 因此有利于后续的基于时域的脑电信号 分析方法的处理。 0054 以上所述实施例的各技术特征可以进行任意的组合， 为使描述简洁， 未对上述实 施例中的各个技术特征所有可能的组合都进行描述， 然而， 只要这些技术特征的组合不存 在矛盾， 都应当认为是本说明书记载的范围。 0055 以上所述实施例仅表达了本发明的几种实施方式， 其描述较为具体和详细， 但并 不能因此而理解为对发明专利范围的限制。 应当指出的是， 对于本领域的普通技术人员来 说， 在不脱离本发明构思的前提下， 还可以做出若干变形和改进， 这些都属于本发明的保护 范围。 因此， 本发明专利的保护范围应以所附权利要求为准。 说 明 书 5/5 页 8 CN 106236083 B 8 图1 说 明 书 附 图 1/3 页 9 CN 106236083 B 9 图2 说 明 书 附 图 2/3 页 10 CN 106236083 B 10 图3 说 明 书 附 图 3/3 页 11 CN 106236083 B 11 。