[发明专利]基于脑电棘波频次的预警方法和装置

权利要求书

1.一种基于脑电棘波频次的预警方法，其特征在于，步骤如下：

1)利用脑电放大器和数据采集板采集脑电信号，将采集到的脑电信号输入计算机；

2)计算机对脑电信号滤波，去除伪迹和噪声；

3)计算机对脑电棘波检测：对脑电信号分段，检测每段脑电信号中的棘波；

4)计算机自动记录一段时间内的累计棘波，然后计算棘波频次SR；

5)获得平滑后的平滑棘波频次SR_m：使用平滑滤波器平滑棘波频次SR，得到SR_m；

6)报警：将平滑棘波频次SR_m与相应阈值进行比较，获得预警结果。

2.如权利要求1所述的基于脑电棘波频次的预警方法，其特征在于，所述步骤1)中所述的脑电放大器采用Neurofile NT脑电放大器，所述的数据采集板采用16位A/D转换数据采集板，采样频率为256Hz。

3.如权利要求1所述的基于脑电棘波频次的预警方法，其特征在于，所述步骤2)所述的对脑电信号滤波使用巴特沃斯带通滤波器，巴特沃斯带通滤波器的频率范围为0.5-30Hz。

4.如权利要求1所述的基于脑电棘波频次的预警方法，其特征在于，步骤3)中所述的棘波检测的步骤为：

a.对脑电信号分段，每段脑电信号长度相等，其长度为N点；

b.利用形态学滤波器检测一段脑电信号f(n)(n＝0，1，…，N-1)的棘波，具体步骤为：

①先对该段脑电信号f(n)关于结构元素g(m)＝300(1-|m-L|/L)(其中：m＝0，1，...，M-1；L＝(M-1)/2；M为结构元素的长度)进行形态开-闭(OC)和闭-开(CO)滤波，公式分别为：

oc(f(n))＝fog·g    和

co(f(n))＝f·gog

其中开运算ο和闭运算·公式分别为：

和

(f·g)=(f⊕g)Θg]]>

其中(fΘg)为腐蚀运算，为膨胀运算，公式分别为：

(fΘg)=minm=0,1,...M-1{f(n+m)-g(m)}(n=0,1,...,N-M)]]>和

(f⊕g)=maxm=0,1...M-1{f(n-m)+g(m)}(n=M-1,M,...,N-1)]]>

②利用以下公式得到背景脑电信号y(n)：

y(n)＝[oc(f(n))+co(f(n))]/2

③得到检出的棘波信号为：z(n)＝f(n)-y(n)。