一种基于相关函数的多通道脑电信号压缩方法与流程

本发明涉及脑电信号处理技术，尤其涉及一种基于相关函数的多通道脑电信号压缩方法。
背景技术：
：随着微电子、生物医学、脑机接口等领域的发展，以无线体域网(wirelessbodyareanetwork，wban)为中心，对目标对象进行远程生理信号监测，包括脑电信号(eeg)、心电信号(ecg)等逐渐成为常态。eeg信号能够反映出大脑的神经电活动，并且不具有伪装性，常常用于癫痫等脑疾病的诊断以及人体的情绪识别和身份识别。在上述领域中，通常需要对eeg信号进行多通道长时间的测量以得到大量的数据，并从中提取相应的有用信息。然而，由于脑电信号的长时间测量和记录会产生大量的数据，导致目前的测量设备需要设计足够容量的存储空间才能对大量数据进行全部存储，而这则会导致随身携带的测量设备的体积较大，降低其携带便利性。因此，需要设计一种对多通道脑电信号进行压缩的方法，实现减少传输的数据量，延长无线设备的使用时间，降低设备的传输带宽。技术实现要素：为了解决上述技术问题，本发明提供一种基于相关函数的多通道脑电信号压缩方法，通过去相关差分运算消除临近通道的相关性，以减少脑电信号数据的存储空间。本发明采用如下技术方案实现：一种基于相关函数的多通道脑电信号压缩方法，包括：计算每个通道的自相关函数和任意两个通道间的互相关函数；利用自相关函数和互相关函数，计算任意两个通道的差分信号和原信号的能量差值；基于能量差值，得到参考/编码通道顺序表；根据参考/编码通道顺序表，为每个通道指定唯一的参考通道，构造差分编码树；基于差分编码树对各通道进行差分运算得到各通道的差分信号；对差分信号中所包含的数据进行霍夫曼编码后得到压缩数据。优选地，通道的自相关函数和两个通道间的互相关函数的计算方法如下：对于通道数为m，每通道样本数为l的多通道脑电信号，先求每个通道中l个样本值的平方和，再对样本数求平均即得到其自相关函数；两个通道对应时刻样本值之积求和，再对样本数求平均得到互相关函数；设脑电信号通道i的l个样本值为xi(1)、xi(2)、...、xi(l)，通道j的l个样本值为xj(1)、xj(2)、...、xj(l)，则通道i和j自相关函数公式为：通道i和j之间的互相关函数公式为：优选地，任意两个通道的差分信号和原信号的能量差值计算步骤包括：设编码通道j的参考通道为p(j)，则编码通道j和参考通道p(j)各自的自相关函数乘以通道信号的样本数求得编码通道j和参考通道p(j)原信号的能量；编码通道j与通道参考p(j)对应时刻样本值之差的平方和，即为差分信号的能量，差分信号能量与参考通道p(j)原信号的能量之和，即为差分编码信号总能量，再减去编码通道j和参考通道p(j)的总能量，得到编码通道j和通道参考p(j)两个通道的差分信号和原信号的能量差值。优选地，参考/编码通道顺序表通过以下步骤得到：把所有能量差值进行升序排序，并对每个能量差值所对应的参考通道和编码通道的顺序也做同样的调整，得到参考/编码通道顺序表。优选地，对能量差值进行升序排序采用的方法为快速排序算法，以第一个能量差值作为枢轴变量。优选地，构造差分编码树的步骤包括：设待编码通道总数为m，首先从参考/编码通道顺序表选出第一对参考通道p(j1)和编码通道j1，若j1尚未指定参考通道，且p(j1)指向j1的连接与已有参考通道-编码通道连接不形成回路，则把p(j1)指定为j1的参考通道，并把待编码通道总数减1；否则丢弃该通道对，顺序选择下一对通道；重复直到待编码通道总数减少为1，表示只剩一个通道还未指定参考通道，把剩下的这一个通道指定为初始参考通道，差分编码树构造完成。优选地，基于差分编码树对各通道进行差分运算得到各通道差分信号的步骤包知：根据差分编码树，除初始参考通道外，把每个通道的信号样本值减去对应参考通道的信号样本值，得到各通道的差分信号。优选地，对差分信号所包含的数据进行霍夫曼编码后得到压缩数据的步骤包括：通过统计排序单元采用边输入边统计排序的方式来对差分信号中的字符出现频率进行统计，并且按照字符频率对字符进行排序；通过编码单元根据差分信号中字符的频率，对排序后的字符进行霍夫曼编码后，通过数据输出单元将编码后得到的数据输出。优选地，数据输出单元先输出霍夫曼编码表，然后再输出霍夫曼编码后得到的数据。优选地，统计排序单元将输入的差分信号按8bits为单位进行划分后进行字符的统计排序。与现有技术相比，本发明至少具有以下有益效果：本发明基于多通道脑电信号临近通道间的相关性，采用差分编码的方式消除通道间的相关性。与原信号相比，得到的差分信号的动态范围和能量均有减小，有效地压缩了信号值的范围。然后再对差分编码后的数据进行霍夫曼编码，减少表示脑电信号所需的数据量，从而降低传输功耗和存储消耗。附图说明图1为本发明一个实施例中利用相关函数对多通道脑电信号进行压缩的流程图；图2为本发明一个实施例中构造差分编码树流程图。具体实施方式下面结合附图和具体实施例对本发明做进一步的详细说明。对于以下实施例中的步骤编号，其仅为了便于阐述说明而设置，对步骤之间的顺序不做任何限定，实施例中的各步骤的执行顺序均可根据本领域技术人员的理解来进行适应性调整。本发明的基本思路包括：对于多通道脑电信号，首先计算各个通道的自相关函数和任意两个通道间的互相关函数，然后根据自相关和互相关函数计算任意两个通道间的差分信号与原信号的能量差值，然后对所有能量差值进行升序排序，并调整对应的参考通道和编码通道顺序，得到参考/编码通道顺序表；根据通道顺序构造差分编码树，并按照差分编码树对各通道进行差分运算得到差分信号；最后对差分信号进行霍夫曼编码，得到压缩数据。本发明以基于国际脑电图学会建议标准的10-10系统和bci2000脑机接口采集的多通道脑电信号作为原始数据库，该脑电信号共64通道，采样频率为160hz，记录时间为2分钟。利用图1所示的压缩方法对其进行数据压缩，包括：s1、计算每个通道的自相关函数和任意两个通道间的互相关函数，对于通道数为m＝64，每通道样本数为l＝160的多通道脑电信号，先求每个通道中l个样本值的平方和，再对样本数求平均即得到其自相关函数；两个通道对应时刻样本值之积求和，再对样本数求平均记得到互相关函数。设脑电信号通道i的l个样本值为xi(1)、xi(2)、...、xi(l)，通道j的l个样本值为xj(1)、xj(2)、...、xj(l)，则通道i和j自相关函数公式为：通道i和j之间的互相关函数公式为：s2、利用步骤s1得到的64个自相关函数和个互相关函数，求出每个通道原信号的能量和任意两个通道间差分信号的能量差值。样本数l＝160乘以通道自相关函数即为该通道的原信号能量，以通道p(j)作为参考通道，通道j作为编码通道(p(j)≠j)，那么编码通道j与参考通道p(j)对应时刻样本值之差的平方和，即为差分信号的能量差分信号的能量与参考通道p(j)能量ep(j)之和，即为差分编码信号总能量再减去编码通道j和参考通道p(j)的总能量即为二者能量差值δej,p(j)。s3、根据步骤s2的所有共4032个能量差值δej,p(j)，求出参考/编码通道顺序表。首先以第一个能量差值作为枢轴变量，对所有δej,p(j)采用快速排序算法进行升序排序，再把每个能量差值δej,p(j)对应的编码通道j和参考通道p(j)按照相同的顺序进行调整，得到参考/编码通道顺序表。s4、根据参考/编码通道顺序表，为每个通道指定唯一的参考通道，构造差分编码树。如图2所示，设待编码通道总数为m＝64，首先选出顺序表中的第一对参考通道p(j1)和编码通道j1，若j1尚未指定参考通道，且p(j1)指向j1的连接与已有参考通道-编码通道连接不形成回路，则把p(j1)指定为j1的参考通道，并把待编码通道总数减1；否则为j1指定参考通道，待编码通道总数保持不变，丢弃能量差值并选择下一对参考/编码通道。s5、重复执行步骤s4，直到待编码通道总数减少为1，表示只剩一个通道还未指定参考通道，把该通道指定为初始参考通道，共有63对参考/编码通道对:j1--p(j1)，j2--p(j2)，……，j63--p(j63)，则差分编码树构造完成。s6、初始参考通道信号保持不变，其它通道的信号样本值xj(n)减去其参考通道的信号样本值xp(j)(n),，得到每个通道的差分信号：s7、对差分信号中所包含的数据进行霍夫曼编码后得到压缩数据。通过统计排序单元采用边输入边统计排序的方式对差分信号中的字符出现频率进行统计，并且按照字符频率从大到小对字符进行排序。具体地，将差分信号数据按顺序输入霍夫曼编码模块，数据输入有两个通道，一方面将输入的数据暂存，以用于后面编码输出，另一方面则输入到统计排序单元进行统计排序。统计排序单元将输入的数据按8bits为单位进行划分后进行字符的统计排序，其中的排序模块会根据输入的数据把对应的频率计数器数值加1，即相当于对输入的相同数据进行个数统计。直至差分信号数据输入完毕，完成字符统计排序。s8、通过编码单元根据差分信号中字符的出现频率，对排序后的字符进行霍夫曼编码后，通过数据输出单元先输出霍夫曼编码表，然后再输出编码，得到压缩数据。为了便于理解，以下将以5个字符为例说明霍夫曼编码的处理过程。假设一个信号包含的所有字符组成的集合x＝{x1,x2,x3,x4,x5}，其概率为p(x1)＝0.1、p(x2)＝0.2、p(x3)＝0.25、p(x4)＝0.4、p(x5)＝0.05。根据霍夫曼编码概率高码字短，概率低码字长的原则，首先按照概率的降序顺序对所有符号进行排列，结果如表1所示，表中w(xi)表示xi的码字。表1符号频率码字x40.4w(x4)x30.25w(x3)x20.2w(x2)x10.1w(x1)x50.05w(x5)每次把频率最低的两个字符合并，将其频率相加，重新为新字符分配码字，对新字符集按频率从大到小排序，直到仅剩一个字符，此时频率为1。最终这一信号的霍夫曼编码由表2给出：表2符号频率码字x40.41x30.2501x20.2000x10.10010x50.050011编码结果表明，出现可能性最大的字符x4分配的码字长度最短，并且所有的码字都符合前缀编码的特点，这就保证了给定任何一个编码序列，都只能够通过唯一的方式解码。由于差分信号的动态范围小于原信号，因此可以用更少的比特数表示每一个信号样本值，实现数据压缩目的。实验采用一组14个完成动作任务的脑电信号作为样本，信号通道数量为64，长度为160，进行对比测试。信号的压缩效果以压缩比进行衡量。实验结果如表3所示。采用未去相关压缩方法的平均压缩率约为2.20，采用本发明去相关压缩方法的平均压缩率约为2.40，平均压缩率提升了约9.12％，证明本发明多通道去相关压缩方法可以进一步提高脑电信号的压缩率。表3任务12345压缩率(未去相关)2.282.152.202.152.08压缩率(去相关)2.412.382.412.402.35任务678910压缩率(未去相关)2.212.202.202.262.22压缩率(去相关)2.402.422.392.402.40任务11121314-压缩率(未去相关)2.232.232.152.22-压缩率(去相关)2.442.402.362.41-上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。当前第1页12