本发明涉及辅助睡眠技术领域,特别是涉及一种睡眠状态分析中去除眼电伪迹的设备。
背景技术:
在睡眠中,人体进行了自我放松及恢复的过程,因此良好的睡眠是保持身体健康的一项基本条件;但是由于工作压力大、生活作息不规律等原因,导致了部分人群的睡眠质量欠佳,表现为失眠、半夜惊醒等。
目前市面上已经有一些设备来帮助人们入睡,提高睡眠质量。例如在某一特定睡眠状态下通过声音、光信号等人工干预,避免在熟睡状态下叫醒用户等。对于辅助睡眠的设备而言,为了真正达到提高用户睡眠质量的目的,正确的识别用户的睡眠状态是非常重要的。
多导睡眠图(Polysomnography,PSG),又称睡眠脑电图,是目前临床上用于睡眠诊断和分析的“金标准”。多导睡眠图利用多种生命体征例如脑电、肌电(颌下)、眼电、呼吸、血氧等对睡眠进行分析。在这些体征信号中,脑电图(electroencephalogram,EOG)处于核心地位。脑电图是利用精密的电子仪器,在头皮上将来自大脑皮层产生的电活动加以记录并放大的波形信号。由于脑电图的信号非常微弱(微伏级),容易被来自其他部位的生物电信号干扰。当眼电信号幅值较低时(即没有较强烈眼球/眼睑活动如眨眼等),眼电信号对脑电信号的干扰比较微弱。而眼电信号幅值较高时,由于眼电信号的频率比正常脑电信号低,高幅值的眼电信号叠加在脑电信号上就形成了一个类似于基线漂移的现象。
为了降低眼电信号所带来的影响,目前有很多去除眼电伪迹的方法。独立成分分析(Indepdent component analysis,ICA)是一种常用的方法。它首先假设输入信号都是统计独立的非高斯的信号的线性组合,然后利用线性变换将来自于信号分离。它的缺点是(1)输入信号的假设条件在实际使用中并不能完全满足;(2)对于分离后的多个信号,还需要进一步判断哪些信号是“纯净的”脑电信号,哪些信号是被分离出的眼电信号。此外,还有方法假设了一个眼电信号对脑电信号的影响因子(如0.2),然后利用脑电信号减去乘以影响因子的眼电信号的方法去除眼电伪迹,如公式:EEGpure=EEGoriginal-0.2*EOG,由于存在个体差异及眼电电极的位置的不同,一个固定的影响因子并不能很好的适应不同的个体。
此外,由于在睡眠状态分析中,脑电信号的波形是一个很重要的睡眠状态指标。例如纺锤波和K复合波的出现,表示进入了非眼快动睡眠的S2期。经过传统方法处理后的脑电信号的波形往往会发生变化,影响了后续对脑电信号的分析效果。
技术实现要素:
基于此,有必要针对上述问题,提供一种睡眠状态分析中去除眼电伪迹的设备,减少去除眼电伪迹过程对脑电信号的波形的影响,确保后续对脑电信号的分析效果。
一种睡眠状态分析中去除眼电伪迹的设备,包括:脑电电极、眼电电极、参考电极、模数转换器、滤波电路以及处理器;
所述脑电电极、眼电电极、参考电极分别连接模数转换器,并依次通过所述模数转换器和滤波电路连接至处理器;
所述脑电电极用于检测用户在睡眠中的原始脑电信号;所述眼电电极用于采集用户在睡眠中的眼电信号;
所述模数转换器将眼电信号和脑电信号转换为数字信号,所述滤波电路对眼电信号和脑电信号进行低频滤波后输入至处理器;
所述处理器,用于对每帧原始脑电信号进行经验模态分解,将其分解成若干个本征模函数,计算各个本征模函数与同一时刻的眼电信号之间的相关系数;查找并删除相关系数大于预设阈值的本征模函数和相关系数最大的本征模函数,利用剩余的本征模函数重建每帧脑电信号。
上述睡眠状态分析中去除眼电伪迹的设备,利用眼电电极采集的用户的眼电信号和脑电电极采集的原始脑电信号,经过模数转换和滤波处理后,由处理器对每帧原始脑电信号进行经验模态分解,将其分解成若干个本征模函数,计算各个本征模函数与同一时刻的眼电信号之间的相关系数;查找并删除相关系数大于预设阈值的本征模函数和相关系数最大的本征模函数,利用剩余的本征模函数重建每帧脑电信号。该设备可以减少去除眼电伪迹过程对脑电信号的波形的影响,保留了原始信号的大部分细节信息,确保后续对脑电信号的分析效果。
附图说明
图1为一个实施例的睡眠状态分析中去除眼电伪迹的设备的结构示意图图;
图2是处理器去除眼电伪迹的算法流程图;
图3是去除眼电伪迹的实验数据结果示意图。
具体实施方式
下面结合附图阐述本发明的睡眠状态分析中去除眼电伪迹的设备的实施例。
参考图1所示,图1为本发明的睡眠状态分析中去除眼电伪迹的设备的结构示意图,包括:脑电电极、眼电电极、参考电极、模数转换器、滤波电路以及处理器;
所述脑电电极、眼电电极、参考电极分别连接模数转换器,并依次通过所述模数转换器和滤波电路连接至处理器;
所述脑电电极用于检测用户在睡眠中的原始脑电信号;所述眼电电极用于采集用户在睡眠中的眼电信号;
所述模数转换器将眼电信号和脑电信号转换为数字信号,所述滤波电路对眼电信号和脑电信号进行低频滤波后输入至处理器;
所述处理器,用于对每帧原始脑电信号进行经验模态分解,将其分解成若干个本征模函数,计算各个本征模函数与同一时刻的眼电信号之间的相关系数;查找并删除相关系数大于预设阈值的本征模函数和相关系数最大的本征模函数,利用剩余的本征模函数重建每帧脑电信号。
上述实施例的睡眠状态分析中去除眼电伪迹的设备,利用眼电电极采集的用户的眼电信号和脑电电极采集的原始脑电信号,经过模数转换和滤波处理后,由处理器对每帧原始脑电信号进行经验模态分解,将其分解成若干个本征模函数,计算各个本征模函数与同一时刻的眼电信号之间的相关系数;查找并删除相关系数大于预设阈值的本征模函数和相关系数最大的本征模函数,利用剩余的本征模函数重建每帧脑电信号。该设备可以减少去除眼电伪迹过程对脑电信号的波形的影响,保留了原始信号的大部分细节信息,确保后续对脑电信号的分析效果。
后续可以利用该设备输出的脑电信号进行睡眠状态监测和分析等,当然,该后续的处理也可以在所述处理器上去实现。
在一个实施例中,所述脑电电极设置在用户的额头位置;所述参考电极设置在用户的耳垂;所述眼电电极设置在眼角位置;如图1所示,图中,脑电电极即图中的“M”,眼电电极包括左右两个电极,即图中的“ROC”和“LOC”,参考电极设置在用户的耳垂,即图中“R”和“L”,加速度传感器即图中“AT”。滤波电路主要是进行低通滤波和滤除工频干扰,为了适应于脑电信号和眼电信号的处理,滤波电路滤波后,输出0-256Hz频段的信号至处理器。
对于去除眼电伪迹功能,主要通过处理器来进行,基于处理器实现的功能,可以在处理器中配置相应的算法模块。
处理器去除眼电伪迹的算法流程包括(1)~(5),具体如下:
(1)处理器控制眼电电极和脑电电极根据设定帧长度采集用户的眼电信号和原始脑电信号;
如在对用户进行辅助睡眠等睡眠状态分析中,处理器可以以设定帧长度,通过用户佩戴的眼电电极和脑电电极,采集用户在睡眠过程中产生的眼电信号和脑电信号。在采集信号时,可以以30s为一帧进行采集,后续对每帧眼电信号和脑电信号进行分析处理。
(2)对该帧原始脑电信号进行经验模态分解,将其分解成若干个本征模函数,得到本征模函数集合;
在此,处理器对脑电信号进行经验模态分解,将其分解成若干个本征模函数(Intrinsic Mode Function,IMF)和残差函数(Redisual,Re)之和的形式。
本征模函数集合包括如下公式:
式中,EEGoriginal表示原始脑电信号,imfi表示第i个本征模函数,Re表示残差函数。
(3)分别计算所述本征模函数集合的各个本征模函数与同一时刻的眼电信号之间的相关系数;
参考图2,图2是处理器去除眼电伪迹的算法流程图,原始脑电信号进行经验模态分解后,得到本征模函数集合,分别计算本征模函数1-n(imf1~imfn)与眼电信号EOG的相关系数1-n(corrcoef1~corrcoefn)。
(4)查找出相关系数大于预设阈值的本征模函数和相关系数最大的本征模函数,并将其从本征模函数集合中删除;
如图2所示,通过设定阈值,在计算完相关系数后,将相关系数大于阈值的本征模函数和相关系数最大的本征模函数删除,剩下的m个本征模函数。
作为一个实施例,在计算完相关系数后,处理器还可以用于在相关系数大于预设阈值的本征模函数中;计算本征模函数与眼电信号的欧氏距离;从欧氏距离最小的本征模函数从本征模函数集合中剔除。
类似的,在计算完相关系数后,处理器还可以用于在相关系数大于预设阈值的本征模函数中;计算本征模函数与眼电信号的余弦距离;从余弦距离最小的本征模函数从本征模函数集合中剔除。
通过上述实施例,在相关系数判断基础上结合了欧氏距离或余弦距离判断,可以将与相关系数判断中无法去除的更多遗留的眼电伪迹去除。
(5)利用本征模函数集合中剩余的本征模函数重建该帧脑电信号;
处理器在去除了眼电伪迹后,利用剩下的m个本征模函数重建去除了眼电伪迹后的脑电信号。作为一个实施例,处理器在重建脑电信号时,根据本征模函数的排列顺序,选择本征模函数集合中位置靠前的若干个本征模函数进行重建脑电信号。
该实施例中,由于本征模函数的排列顺序是按频率由大到小,并且与眼电信号相似度最高的本征模函数一般排列在中间位置,因此在重建脑电信号时,可以仅利用前若干个本征模函数,删除包括相似度最高的本征模函数在内的频率较低的本征模函数后,再重建脑电信号。
其中,处理器可以采用如下公式重建脑电信号:
式中,EEGpure表示重建的脑电信号,corrcoef表示相关系数,imf表示第i个本征模函数,EOG表示眼电信号,corrcoefmax表示最大的相关系数,thre表示预设的相关系数阈值。
在一个实施例中,处理器可以在对原始脑电信号进行经验模态分解前,先将原始脑电信号帧划分为N个时间窗口,并对每个时间窗口的脑电信号进行本征模函数分解;以及在重建脑电信号后,将各个时间窗口重建的脑电信号进行合并,得到脑电信号帧。
上述实施例,通过将采集的原始脑电信号帧划分多个时间窗口并行处理,能够加快信号处理速度,提高睡眠状态分析的效率。
例如,以30s一帧为例,可以以5s或10s为一个时间窗口长度。
本发明的睡眠状态分析中去除眼电伪迹的设备,只去除高幅度眼电所造成的类似于基线漂移的伪迹,并且保留了原始信号的大部分细节信息;后续使用该脑电信号去进行睡眠状态分析时,可以得到更好的效果。
参考图3所示,图3是去除眼电伪迹的实验数据结果示意图。图3(a)为采集的原始脑电信号,图3(b)为采集的眼电信号,图3(c)比较了去除眼电伪迹前后的脑电信号(图中,①为原始脑电信号;②为去除眼电伪迹后的脑电信号),下图是上图截取部分放大图,可以发现,上述区间内的数据点之间,由眼电带来的幅度较大的深V形波动被本方案给消除的同时,并且保留了较多的原始信息。
相对于传统方法(如ICA等)将输入的多路信号视为经过线性组合后的多路源信号,并试图将这些信号彼此分离,能在周期信号上传统方法能获得比较好的效果。而且对于脑电信号而言,由于脑电信号和眼电信号都可以视为随机信号,并且脑电信号容易受到外部干扰,很难将脑电信号和眼电信号彻底分离开,此时脑电信号就会混入额外的噪声信号,加大了后续信号处理分析的难度。而本发明的技术方案,只去除高幅度眼电带来的类似于基线漂移的伪迹,保留了原始信号的大部分细节信息。因此有利于后续的基于时域的脑电信号分析方法的处理。
以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。