1.一种清醒状态下的脑电信号数据类型的标注方法,其特征在于,包括:
在用户开始睡眠过程后,采集用户的实时眼电信号和脑电信号样本;
分别对所述实时眼电信号和脑电信号样本进行小波分解,并根据设定低频段的小波系数进行信号重建得到眼电信号和脑电信号;
根据同一时刻脑电信号和眼电信号的相关性以及眨眼眼电波形的特征,在眼电信号上检测眨眼活动;
当检测到眨眼活动时,计算所述脑电信号的样本熵,将该样本熵与预先计算的样本熵阈值进行比较,若所述样本熵大于所述样本熵阈值,将所述脑电信号样本的信号类型标注为清醒状态;
所述样本熵阈值的计算公式如下:
sampen_vali=sampen(y[p_start:p_end])
p_start=(i-1)*time_length*fs+1
p_end=p_start+time_length*fs-1
p_end<T·fs
式中,其中sampen_thre为样本熵阈值,sampen_vali为样本熵集合中第i个样本的样本熵,sampen为求样本熵的运算,其输入y[p_start:p_end]为脑电信号y在第p_start点开始到第p_end点为止的部分,time_length为计算样本熵的每个样本的时间长度,fs为脑电信号的采样率,T为开始采集脑电信号后的设定时间,v为设定参数。
2.根据权利要求1所述的清醒状态下的脑电信号数据类型的标注方法,其特征在于,所述根据同一时刻脑电信号和眼电信号的相关性以及眨眼眼电波形的特征,在眼电信号上检测眨眼活动的步骤包括:
利用具有设定信号幅度和时间长度的滑动窗口截取低频眼电信号;
分别计算滑动窗口内眼电信号与同一时刻脑电信号的波形的相关系数,滑动窗口内眼电信号波形尖峰的尖锐程度参数以及尖峰的持续时间;
若所述相关系数、尖锐程度参数和持续时间均分别满足预设的相关系数阈值、尖锐程度参数阈值和持续时间阈值,判断该滑动窗口内眼电信号存在眨眼活动。
3.根据权利要求2所述的清醒状态下的脑电信号数据类型的标注方法,其特征在于,所述计算滑动窗口内眼电信号波形尖峰的尖锐程度参数的步骤包括:
分别计算眼电信号波形在滑动窗口内的上部区域面积和下部区域面积,计算公式如下:
式中,pi为滑动窗口内的眼电信号,pmax为滑动窗口内眼电信号的最大值,pmin为滑动窗口内眼电信号的最小值,areaup表示上部区域面积,areadown表示下部区域面积;
根据所述上部区域面积和下部区域面积计算所述眼电信号波形尖峰的面积,计算公式如下:
式中,blinkarea表示尖峰的面积,if表示满足条件;
根据尖峰面积计算尖锐程度参数,计算公式如下:
blinkratio=blinkarea/in-blinkarea
式中,blinkratio表示尖锐程度参数,in-blinkarea表示非尖峰部分的面积。
4.根据权利要求3所述的清醒状态下的脑电信号数据类型的标注方法,其特征在于,所述计算滑动窗口内眼电信号波形尖峰的持续时间的步骤包括:
根据所述上部区域面积和下部区域面积计算所述眼电信号波形尖峰的方向,计算公式如下:
式中,blinkdirection表示尖峰方向,标记为1表示尖峰向下,标记为-1表示尖峰向上;
在尖峰向下时,根据尖峰两侧的局部极大值点计算尖峰的持续时间;或者在尖峰向上时,根据尖峰两侧的局部极小值点计算尖峰的持续时间;计算公式如下:
blinkts=vertexright-vertexleft,
式中,blinkts表示尖峰的持续时间,vertexright表示尖峰的右侧顶点时刻,表示尖峰的左侧顶点时刻,right_min_loca表示尖峰右侧的局部极小值点,right_max_loca表示尖峰右侧的局部极大值点,left_min_loca表示尖峰左侧的局部极小值点,left_max_loca表示尖峰左侧的局部极大值点,s.t.表示约束条件。
5.根据权利要求4所述的清醒状态下的脑电信号数据类型的标注方法,其特征在于,所述滑动窗口的幅度为75微伏至300微伏,所述眨眼时眼电信号与脑电信号之间的相关系数阈值为0.9,所述尖锐程度参数阈值为0.3,所述持续时间阈值为0.3秒。
6.根据权利要求1至5任一项所述的清醒状态下的脑电信号数据类型的标注方法,其特征在于,还包括:
从开始采集脑电信号样本和实时眼电信号后的设定时间T内,将所有脑电信号样本标注为清醒状态。
7.根据权利要求1所述的清醒状态下的脑电信号数据类型的标注方法,其特征在于,所述参数v的取值如下:
其中,
x=v
式中,X表示样本熵集合,Xi表示样本熵集合X中第i个元素,Φ(x)表示标准正态分布函数的积分,P(X≤x)表示样本熵的集合X中的取值小于x的概率。
8.一种清醒状态下的脑电信号数据类型的标注系统,其特征在于,包括:
信号采集模块,用于在用户开始睡眠过程后,采集用户的实时眼电信号和脑电信号样本;
信号重建模块,用于分别对所述实时眼电信号和脑电信号样本进行小波分解,并根据设定低频段的小波系数进行信号重建得到眼电信号和脑电信号;
眨眼检测模块,用于根据同一时刻脑电信号和眼电信号的相关性以及眨眼眼电波形的特征,在眼电信号上检测眨眼活动;
样本标注模块,用于当检测到眨眼活动时,计算所述脑电信号的样本熵,将该样本熵与预先计算的样本熵阈值进行比较,若所述样本熵大于所述样本熵阈值,将所述脑电信号样本的信号类型标注为清醒状态;
所述样本熵阈值的计算公式如下:
sampen_vali=sampen(y[p_start:p_end])
p_start=(i-1)*time_length*fs+1
p_end=p_start+time_length*fs-1
p_end<T·fs
式中,其中sampen_thre为样本熵阈值,sampen_vali为样本熵集合中第i个样本的样本熵,sampen为求样本熵的运算,其输入y[p_start:p_end]为脑电信号y在第p_start点开始到第p_end点为止的部分,time_length为计算样本熵的每个样本的时间长度,fs为脑电信号的采样率,T为开始采集脑电信号后的设定时间,v为设定参数。