[发明专利]一种脑电爆发抑制指数计算方法及系统

权利要求书

1.一种脑电爆发抑制指数计算方法，其特征在于，包括：

利用离散小波变换，将预设时间段内的脑电信号分解为多层；

根据脑电信号主要能量集的分布区间，对包含脑电信号的频带小波层数系数进行处理，通过自适应的能量阈值提取出每层的噪声系数和信号系数，并重构出噪声波形序列与信号波形序列，所述自适应的能量阈值通过以下公式计算：

λ_dk＝C_dk*PD(k)

式中，k取2、3、4、5，PD(k)表示第k层小波对应的能量，dk表示第k层小波，λ_dk是第k层小波的自适应能量阈值，C_dk是第k层小波的自适应能量阈值计算比例；

对噪声波形序列和信号波形序列分别进行差分处理，并将差分后的信号波形序列和噪声波形序列进行逐点平方，以及利用预设传递函数进行滑动平均，得到信号水平序列和噪声水平序列，所述预设传递函数为

其中，N＝FS，FS为脑电信号的采样率；

基于信号水平序列和噪声水平序列利用预设爆发抑制判定方式，得到抑制信号点，累加抑制信号点对应的抑制信号大于预设持续时间的抑制段，并计算抑制段占总爆发抑制信号段的比例，得到脑电爆发抑制指数，所述预设爆发抑制判定方式包括：

将计算得到的信号水平和噪声水平逐点做比较，当信号水平序列大于噪声水平序列时，判定此点为爆发信号点，否则，则判定为抑制信号点；或，判断信号水平序列是否大于自适应阈值，如果大于等于自适应阈值，则判定为爆发信号点，否则，则判定为抑制信号点。

2.根据权利要求1所述的脑电爆发抑制指数计算方法，其特征在于，所述利用离散小波变换，将预设时间段内的脑电信号分解为多层的步骤，包括：

采用预设母小波并根据脑电信号频带分类，使用离散小波变换将采样率为FS，时长为T秒的脑电信号分解为N层，得到对应的组小波系数和近似系数。

3.根据权利要求2所述的脑电爆发抑制指数计算方法，其特征在于，所述根据脑电信号主要能量集的分布区间，对包含脑电信号的频带小波层数系数进行处理，通过自适应的能量阈值提取出每层的噪声系数和信号系数，并重构出噪声波形序列与信号波形序列，重构出噪声波形序列与信号波形序列的步骤，包括：使用离散小波变换将采样率为256Hz，时长为60s内的脑电图信号进行6层分解，得到6组小波系数D1-D6和最后一组近似系数A6；

将A6和D1、D2的系数置零去除有低频干扰与高频干扰的数据，计算D2-D5的系数对应的能量，则脑电波的能量为D2-D5层对应的能量之和；

对D2层的系数取绝对值后进行从小到大的排序，得到排序后的系数序列SD2，对排序后的系数SD2从小到大逐个进行累加，当累加到第n个数时的值刚好大于等于自适应的能量阈值，把当前的值则作为D2层的系数阈值，并以同样的方式得到D3、D4、D5层的系数阈值；

将各个层系数的绝对值大于对应阈值的系数作为信号系数，小于对应阈值的系数作为信号系数作为噪声系数；

将信号系数进行重构得到信号波形序列，将噪声系数进行重构得到噪声波形序列。