一种锋电位分离方法技术

本发明专利技术揭示了一种锋电位分离方法，包括如下步骤：首先，通过微电极细胞外记录方式采集神经元电信号；其次，将采集到的信号进行滤波处理；之后，采用阈值法对滤波后的信号进行锋电位检测，并得到锋电位信号；然后，采用离散小波变换对锋电位信号进行四级分解，并得到由小波系数组成的五个分量，且这五个分量组成一个小波系数矩阵；接着，对小波系数矩阵进行进一步特征提取得到特征空间；最后，采用局部加权投票法对特征空间进行聚类，实现锋电位信号分离。本发明专利技术在加权投票的基础上，引入香农熵的概念来评价每个类簇的可靠性，并基于其可靠性给出相应的权值，提出了局部加权投票法，并将其运用于锋电位信号分离，该分离方法实现简单，执行高效。

【技术实现步骤摘要】
一种锋电位分离方法
本专利技术涉及一种分离方法，尤其涉及一种锋电位分离方法，属于锋电位分离领域。
技术介绍
人体共有八大系统，分别是运动系统、神经系统、内分泌系统、循环系统、呼吸系统、消化系统、泌尿系统和生殖系统，上述各系统协调配合，使人体内各种复杂的生命活动能够正常进行。其中，神经系统起到主导作用，它不断接收、分析并处理各个感受器传入的信息。处理后的信息输出到人体各系统器官以调节它们的活动。神经元作为神经系统最基本的结构和功能单位，可以通过突触传递和非突触传递方式进行信息的传递，其中前者为主要传递方式。突触传递又分为化学突触传递和电突触传递。医学上常采集神经元电信号进行医学研究，其中动作电位是重点研究对象。由于动作电位主要由锋电位(spike)组成，因此通常所说的动作电位主要指锋电位。神经元产生的电信号可以通过微电极细胞外记录方式得到，这种方式采集到的信号往往是同一脑区内多个神经元电信号的叠加。神经元电信号可以分别响应不同的生理活动，因此知道哪种电信号对应哪种生理活动至关重要。每种神经元电信号都有自己独特的锋电位信号，因此锋电位信号作为神经元电信号的一部分，可以根据其特征对神经元电信号进行区分。如果能对锋电位信号进行分离，锋电位信号所对应的神经元电信号也随之分离。锋电位分离方法(spikesorting)用以实现从采集到的神经元电信号中检测出锋电位信号，并把锋电位信号划分成不同的类簇，使得每个类簇对应单个的神经元。在锋电位分离方法中，检测锋电位的方式是阈值法，虽然阈值法实现简单，但是阈值的设定是一个问题，阈值过大容易造成漏检，阈值过小容易造成误检。后来人们对阈值法进行了改进，从而产生了双向阈值法、峰值检测法等，其中后者很好地解决了电压漂移对于阈值检测的干扰。虽然双向阈值法和峰值检测法在一定程度上对阈值法进行了改进，但是它们依然没有完美地解决阈值或窗口高度的设定问题。综上所述，如何提供一种锋电位分离方法，实现简单、分离质量高且执行高效，就成为本领域技术人员亟待解决的问题。
技术实现思路
本专利技术结合聚类集成和局部加权的思想，提出了一种基于局部加权投票的锋电位分离方法，该方法采用阈值法检测锋电位，使用四级离散小波分解对锋电位信号进行特征提取，最后使用局部加权投票法对得到的特征空间进行聚类，实现锋电位信号分离。本专利技术技术解决方案是：一种锋电位分离方法，包括如下步骤：S1：通过微电极细胞外记录方式采集神经元电信号；S2：对步骤S1中采集到的信号进行滤波处理；S3：采用阈值法对步骤S2中滤波后的信号进行锋电位检测，并得到锋电位信号；S4：采用离散小波变换对步骤S3中的锋电位信号进行四级分解，并得到由小波系数组成的五个分量，完成初步特征提取，且这五个分量组成一个小波系数矩阵；S5：对步骤S4中的小波系数矩阵进行进一步特征提取得到特征空间，提取的特征包括小波系数的均值、小波系数的标准差、小波系数的斜度、小波系数的能量以及小波系数的熵；S6：采用局部加权投票法对所述特征空间进行聚类，实现锋电位信号分离，所述局部加权投票法包括如下步骤：S61：对所述特征空间进行B次随机采样，并得到B个子特征空间；S62：对B个子特征空间分别使用k均值聚类算法进行聚类，得到B个子特征空间聚类结果；S63：根据k均值聚类算法计算出的聚类中心，对所述特征空间进行再分配，得到B个基于所述特征空间的聚类结果；S64：步骤S63中每个聚类结果均包括多个类簇，计算步骤S63中每个聚类结果中每个类簇的局部权重；S65：计算样本被划分为某个类簇的次数与该类簇的局部权重的乘积，选出乘积最大的类簇作为该样本最终的聚类结果。优选地，所述步骤S1中微电极细胞外记录方式采集神经元电信号是将微电极设置在神经组织的神经元之间，通过微电极尖端检测到的电位变化反映出神经元外介质中的电流。优选地，所述微电极为玻璃微电极或金属微电极。优选地，所述步骤S2中滤波处理的方式为高通滤波，所述高通滤波的截止频率为400Hz。优选地，所述步骤S3中所述阈值法中的阈值为信号幅度最大值的二分之一。优选地，所述步骤S4中的四级分解包括第一级离散小波分解、第二级离散小波分解、第三级离散小波分解及第四级离散小波分解；所述第一级离散小波分解包括如下步骤：首先对锋电位信号分别进行低通滤波与高通滤波，然后对滤波后的信号分别进行二分之一降采样得到第一级近似分量A1与第一级细节分量D1；所述第二级离散小波分解包括如下步骤：首先对第一级近似分量A1分别进行低通滤波与高通滤波，然后对滤波后的信号分别进行二分之一降采样得到第二级近似分量A2与第二级细节分量D2；所述第三级离散小波分解包括如下步骤：首先对第二级近似分量A2分别进行低通滤波与高通滤波，然后对滤波后的信号分别进行二分之一降采样得到第三级近似分量A3与第三级细节分量D3；所述第四级离散小波分解包括如下步骤：首先对第三级近似分量A3分别进行低通滤波与高通滤波，然后对滤波后的信号分别进行二分之一降采样得到第四级近似分量A4与第四级细节分量D4；所述由小波系数组成的五个分量分别为第四近似分量A4、第四细节分量D4、第三细节分量D3、第二细节分量D2以及第一细节分量D1，且五者均是由小波系数Cλ组成的行向量，且1≤λ≤N，其中N为小波系数的个数。优选地，所述步骤S5中小波系数的均值的公式、小波系数的标准差的公式、小波系数的斜度的公式、小波系数的能量的公式以及小波系数的熵的公式分别为：小波系数的均值为：小波系数的标准差为：小波系数的斜度为：小波系数的能量为：小波系数的熵为：所述特征空间的维数为每个分量提取的特征个数与分量的个数的乘积。优选地，所述步骤S63中B个基于所述特征空间的聚类结果为Π＝{π1,π2,...,πi,...,πB}，其中，πi为第i次聚类得到的聚类结果；所述步骤S64中所述类簇的聚类表示定义为：其中，CR(πi)是一个只由0和1组成的矩阵，将B个基于所述特征空间的聚类结果的CR合成一个矩阵T，T＝[CR(π1),CR(π2),...,CR(πB)]，且矩阵T是一个n×b的矩阵，其中n代表样本量的大小，pi表示聚类结果πi中类簇的个数，b表示所有聚类结果的类簇个数的总和；优选地，矩阵T的第f列记作类簇Tf，类簇Tf相对于第i个聚类结果πi的不可靠性定义为：其中，为Tf与的内积，表示B个基于所述特征空间的聚类结果中所有类簇的第f个类簇与第i个聚类结果里的第j个类簇拥有相同样本的个数，[Tf,Tf]表示B个基于所述特征空间的聚类结果中所有类簇中的第f个类簇中样本的个数，故如果所有聚类结果相对独立，则Tf相对于所有聚类结果Π的不可靠性为：其中，B表示所有聚类结果的个数，u(Tf,πi)∈[0,log2pi]，则归一化不可靠性nu(Tf,πi)表达式为：其中，NU(Tf,Π)∈[0,1]，则归一化可靠性ND(Tf,Π)为：ND(Tf,Π)＝1-NU(Tf,Π)。优选地，每个所述类簇的局部权重为：与现有技术相比，本专利技术具有如下优点：1.本专利技术采用了聚类集成，聚类集成比单个聚类算法具有更高的鲁棒性，并且能够提高聚类划分的质量。不同的聚类结果从不同方面反映了数据集的结构，聚类集成通过对多个聚类结果进行合并，可以得到一个能较好地反映数据集内在结构的数据本文档来自技高网...

【技术保护点】
1.一种锋电位分离方法，其特征在于，包括如下步骤：S1：通过微电极细胞外记录方式采集神经元电信号；S2：对步骤S1中采集到的信号进行滤波处理；S3：采用阈值法对步骤S2中滤波后的信号进行锋电位检测，并得到锋电位信号；S4：采用离散小波变换对步骤S3中的锋电位信号进行四级分解，并得到由小波系数组成的五个分量，完成初步特征提取，且这五个分量组成一个小波系数矩阵；S5：对步骤S4中的小波系数矩阵进行进一步特征提取得到特征空间，提取的特征包括小波系数的均值、小波系数的标准差、小波系数的斜度、小波系数的能量以及小波系数的熵；S6：采用局部加权投票法对所述特征空间进行聚类，实现锋电位信号分离，所述局部加权投票法包括如下步骤：S61：对所述特征空间进行B次随机采样，并得到B个子特征空间；S62：对B个子特征空间分别使用k均值聚类算法进行聚类，得到B个子特征空间聚类结果；S63：根据k均值聚类算法计算出的聚类中心，对所述特征空间进行再分配，得到B个基于所述特征空间的聚类结果；S64：步骤S63中每个聚类结果均包括多个类簇，计算步骤S63中每个聚类结果中每个类簇的局部权重；S65：计算样本被划分为某个类簇的次数与该类簇的局部权重的乘积，选出乘积最大的类簇作为该样本最终的聚类结果。...

【技术特征摘要】
1.一种锋电位分离方法，其特征在于，包括如下步骤：S1：通过微电极细胞外记录方式采集神经元电信号；S2：对步骤S1中采集到的信号进行滤波处理；S3：采用阈值法对步骤S2中滤波后的信号进行锋电位检测，并得到锋电位信号；S4：采用离散小波变换对步骤S3中的锋电位信号进行四级分解，并得到由小波系数组成的五个分量，完成初步特征提取，且这五个分量组成一个小波系数矩阵；S5：对步骤S4中的小波系数矩阵进行进一步特征提取得到特征空间，提取的特征包括小波系数的均值、小波系数的标准差、小波系数的斜度、小波系数的能量以及小波系数的熵；S6：采用局部加权投票法对所述特征空间进行聚类，实现锋电位信号分离，所述局部加权投票法包括如下步骤：S61：对所述特征空间进行B次随机采样，并得到B个子特征空间；S62：对B个子特征空间分别使用k均值聚类算法进行聚类，得到B个子特征空间聚类结果；S63：根据k均值聚类算法计算出的聚类中心，对所述特征空间进行再分配，得到B个基于所述特征空间的聚类结果；S64：步骤S63中每个聚类结果均包括多个类簇，计算步骤S63中每个聚类结果中每个类簇的局部权重；S65：计算样本被划分为某个类簇的次数与该类簇的局部权重的乘积，选出乘积最大的类簇作为该样本最终的聚类结果。2.根据权利要求1所述的一种锋电位分离方法，其特征在于：所述步骤S1中微电极细胞外记录方式采集神经元电信号是将微电极设置在神经组织的神经元之间，通过微电极尖端检测到的电位变化反映出神经元外介质中的电流。3.根据权利要求2所述的一种锋电位分离方法，其特征在于：所述微电极为玻璃微电极或金属微电极。4.根据权利要求1所述的一种锋电位分离方法，其特征在于：所述步骤S2中滤波处理的方式为高通滤波，所述高通滤波的截止频率为400Hz。5.根据权利要求1所述的一种锋电位分离方法，其特征在于：所述步骤S3中所述阈值法中的阈值为信号幅度最大值的二分之一。6.根据权利要求1所述的一种锋电位分离方法，其特征在于：所述步骤S4中的四级分解包括第一级离散小波分解、第二级离散小波分解、第三级离散小波分解及第四级离散小波分解；所述第一级离散小波分解包括如下步骤：首先对锋电位信号分别进行低通滤波与高通滤波，然后对滤波后的信号分别进行二分之一降采样得到第一级近似分量A1与第一级细节分量D1；所述第二级离散小波分解包括如下步骤：首先对第一级近似分量A1分别进行低通滤波与高通滤波，然后对滤波后的信号分别进行二分之一降采样得到第二级近似分量A2与第二级细节分量D2；所...

【专利技术属性】
技术研发人员：季薇，李振斌，李云，
申请(专利权)人：南京邮电大学，
类型：发明
国别省市：江苏,32