专利名称：神经元信息单位放电分析方法及分析系统的制作方法
技术领域：
本发明涉及生物信息处理技术领域，具体涉及一种神经元单位放电分析方法。
背景技术：
研究神经元的编码与整合，及最终转换成有意义知觉的神经机理，是系统神经科 学的重要内容之一。这一层次的研究目前主要是用微电极核团或神经通路上单个神经元或 者神经元群体的胞外电活动，通过分析特定刺激条件下神经元的放电行为，来解析神经系 统如何进行“神经元”的编码和感知。当前关于神经元信息编码特点解析多采用频率分析，即单位时间内的放电个数， 以及相应分析系统配带的刺激事件前后分析。如文献电刺激初级听皮层对水杨酸耳鸣模 型大鼠下丘外侧核神经元放电的影响(生理学报Acta Physiologica Sinica，April 25， 2009,61(2) 121-126)；以及前扣带皮层对丘脑腹侧基底核神经元自发放电的抑制效应 (生理学报 Acta Physiologica Sinica, June 25，2009，61 (3) :279_284)等，其采用电生理 数据即动作电位分析，主要是用微电极核团或神经通路上单个神经元或者神经元群体的胞 外电活动，通过分析特定刺激条件下神经元的放电行为，来解析神经系统如何进行听觉信 息的编码和感知，按照其文献描述，实现的数据分析指标较少，并且需采用专用数据分析系 统，其通用性和可扩展性不强。
发明内容
本发明针对现有技术中的上述缺陷，利用多导生理仪数字信号powerlab采集控 件，结合Matlab的数据处理功能和图形用户界面⑶I，设计一个用于神经元电生理研究的 实时分析系统并提出一种对神经元单位放电处理的分析方法。本发明解决上述技术问题的技术方案是，提出一种神经元信息单位放电分析方 法，基于神经元动作电位分析、频率分析以及非线性、模式识别技术，通过对采集的原始数 据进行转化，提取相邻动作电位峰电位的时间间隔，进行数据转化，根据不同要求分别得出 其统计值频率、复杂度、优势片段、并绘制出神经节律特点多分析图(放电序列、ISI分布 图、ISI非序列直方图、ISI回归映射图)。具体包括如下步骤生物电信号采集系统将采集到的神经元放电信号送入数据转换单元，将原始采集 的神经元放电信号(原始数据)转化为Matlab数据格式，获得原始间隔序列，峰值选取单 元从神经元优势片段中选取神经元峰电位，确定神经元峰峰间隔；计算单元根据神经元峰 峰间隔计算神经元平均间隔以及平均频率；比较单元将原始间隔序列与平均间隔比较确定 峰电位间隔序列，可将该间隔序列表示为为0、1序列；首先比较单元分别比较上述峰电位 间隔序列的重复性，得到被分段的数目复杂度c (η)，以此直接表达数据串的复杂程度并将 之量化；其次图形单元根据神经元放电的节律性，绘制峰电位间隔ISI时间间隔图；并且根 据神经元的种类设定带宽值bin，根据带宽值bin绘制ISI非序列间隔直方图；按照原始数据的时间顺序相邻两个放电时间差为放电间隔，以相邻两个放电间隔为横坐标，以下一相 邻两个放电间隔为纵坐标，绘制ISI回归映射图，判定神经元的节律特点(可以按照神经元 节律特性首先判断是否具有周期性，如果具有周期性再判断其节律是周期1，还是周期2， 还是周期3等等)。优势片段的获取可采用如下方法，多导生理记录仪中搜索单元可用滑移法在相邻 间隔的范围内寻找各种神经元放电间隔排序模式，根据神经元的间隔模式确定间隔范围 数，在间隔范围数个相邻间隔内寻找神经元放电间隔排序模式，找出重复频率大于重复次 数(可为3次)的间隔组合作为优势片段的原始序列集合；通过平均ISI比对产生替代数 据序列(可计算峰间隔值，如实际ISI >平均ISI则由1替代，小于平均ISI则由0替代，最 后原始数据转化为以0，1为数据的替代数据)替代原始数据序列；搜索单元寻找替代数据 序列中的优势片段作为随机序列集合；比较单元比较原始序列集合与随机序列集合，如果 某一片段在随机序列集合中的出现率大于在原始序列集合中的出现率，则在随机序列集合 中剔除该片段，留下原始序列集合中的片段，如果不是相邻的η片断的重复即为优势片段。本发明还提出一种神经元单位放电分析系统，包括数据采集单元获取优势片段， 峰值选取单元从中选择神经元峰电位，确定神经元峰-峰间隔；根据神经元峰_峰间隔计算 神经元平均间隔，比较单元将神经元原始间隔序列与平均间隔序列比较转化为峰电位间隔 序列；比较单元分别比较上述峰电位间隔序列的重复性，得到被分段的数目复杂度c(n)； 图形用户界面工具以放电时间为横坐标，即时对应的相邻两个放电间隔为纵坐标绘制ISI 时间间隔图；根据神经元的种类设定带宽值bin，根据带宽值bin绘制ISI非序列间隔直方 图；以相邻两个放电间隔为横坐标，以下一相邻两个放电间隔为纵坐标，绘制ISI回归映射 图，判定神经元的节律特点。本发明基于神经元放电信息多样化的特点将非线性、模式分析方法相结合，将神 经元的筛选数据进行充分的利用。通过本发明将数据进行替代处理，提取神经元放电间隔， 应用非线性理论公式计算出放电序列复杂度、频率、序列分布图、序列回归映射图、序列直 方图等能充分反应放电的周期性、节律性、动态性变化等特征。该系统可以实时或后期读取 powerlab记录到的神经元放电信号和刺激参数信息，并给出基本的分析结果，为神经电生 理实验提供了一种方便和定量的工程评测手段，可以帮助实验者更有效率、实时处理并应 用多种线性、非线性方法寻找有特异神经元或者特异性的放电特征。
图1为本发明所述单位放电分析方法流程示意图
具体实施例方式本系统基于多导生理记录仪Powerlab作为生物电信号采集系统，该系统通过 DSP硬件模块实现了数字信号的在线处理，包括放大、滤波，以及神经元的筛选(spike sorting)。生物电信号采集系统采集到的神经元放电信号经前置放大器放大，并进行模数 转换后，经过滤波后的原始波形和硬件分类后的动作电位的时间点(time stamps)通过高 速数据接口输入单位放电分析系统，贮存在专门数据存储区，采图单元Chart图形采集模 块记录神经元，为用户提供一个把专用系统(Tum-Key)和定制系统(Customer)结合起来的平台。应用Matlab编程语言实现了大量数据分析。Matlab单元直接从chart记录数据库 中读取数据，包括原始波形、动作电位时间点、刺激时间点等，并在Matlab单元中给出神经 元初步分析结果。采用Matlab的图形用户界面(⑶I)工具交互显示数据处理结果。用户 只需输入分析的对象，包括本次实验数据的文件名、特定数据段的编号，以及所存储数据的 标注名称，选择需要的分析手段，提供详细的参数信息，如时间窗的大小(放电时间)、计算 发放率(频率)的时间范围等，就直接可以得到相应的分析结果。分析结果直观且便于进 一步分析。针对实验中刺激参数变化的序列，利用Matlab实现了对刺激参数影响神经元电 活动的实时和定量的评价。本发明通过对采集到的神经元放电信号(原始数据)进行转化，提取相邻动作电 位峰电位的时间间隔，进行数据转化，分别得出其统计值放电频率、放电序列复杂度、优势 片段、以及绘制出神经节律特点多分析图(放电序列、ISKinterspike interverval, ISI 峰电位)间隔分布图、ISI非序列直方图、ISI回归映射图)。多导生理记录仪Powerlab作 为生物电信号采集系统，采集到的神经元放电信号送入数据转换单元，将原始数据转换为 Matlab数据；从神经元优势片段中采集数据，峰值选取单元选取神经元峰电位，确定神经 元峰峰间隔；计算单元根据神经元峰电位计算神经元放电个数、放电时间、频率等；根据峰 峰间隔计算平均频率、平均间隔，并将原始间隔序列与平均间隔比较转化为0、1序列；比较 单元分别比较上述峰电位间隔序列的重复性，得到被分段的序列复杂度c(n)。当放电序列中的峰峰间隔片断模式在某一放电序列中出现的次数明显多于在随 机序列中可能出现的次数，这些片段为优势片段(模式)。确定优势片段的具体步骤如下， 根据神经元的间隔模式确定间隔范围数(可选择3 6个相邻间隔)，搜索单元用滑移法在 相邻间隔内寻找各种神经元放电间隔排序模式，找出重复频率大于重复次数(如3次)的 间隔组合作为优势片段的原始集合，产生替代数据序列，应用滑移寻找替代数据序列中的 优势片段作为随机序列集合；比较单元对原始序列和随机序列进行比较，如果原始序列中 的某一随机组合片段在随机序列中的出现率大于在原始序列中的出现率，则在随机序列中 剔除该随机片段；否则保留该原始序列在随机序列中的片段，检验并排除相邻片段的重复 后确定为优势片段。图形用户界面工具计算绘制非序列间隔直方图，找出间隔分布的带宽值(bin)，将 原始序列标准化为相对ISI值；绘制ISI分布图、ISI非序列直方图、ISI回归映射图。具 体为多导生理记录仪Powerlab作为生物电信号采集系统，数据单元采集单元从神经 元放电序列优势片段中采集数据，峰值选取单元从采集的数据中选择有效波幅阈值区分单 一神经元，确定每个神经元的峰电位，并计算每个神经元峰值间隔。为确保间隔有效性并且 避免神经元单一峰的左右峰间隔计算，通过自行输入点个数进行神经元峰点的选择与自我 检查。比如可以自行输入点个数，且如果间隔小于2ms，则提示重新进行点位输入，重新计算 间隔以确保是放电的最高点(峰点)之间的时间距离。计算出平均ISI值，并将原始ISI 值与上述平均ISI值比对，按照整数倍转化为相对ISI值。绘制间隔直方图用于判定神经 元放电序列的节律特点。绘制ISI分布图按照峰电位间隔时间分布图，以放电时间为横坐标，即时对应的 相邻两个放电间隔为纵坐标画散点图；ISI非序列直方图绘制需根据神经元的种类设定带宽值bin，并以带宽值bin宽为横坐标，相应带宽值bin宽的个数为纵坐标，绘制ISI非序 列间隔直方图。根据神经元的种类设定带宽值bin，将原始序列标准化为相对ISI值(如 峰峰间隔在0 1 · bin ms之间为1，在1 · bin ms 2 · bin ms之间为2...依此类推)。 ISI回归映射图相邻两个放电间隔为横坐标，以下一相邻两个放电间隔为纵坐标，绘制ISI 回归映射图。如图1所示为本发明所述单位放电分析流程示意图，具体步骤如下步骤1.将采集到的神经元放电信号转化为Matlab数据，选取神经元峰电位，转换 为峰峰间隔。放电序列中的神经元峰峰间隔片断模式在某一放电序列中出现的次数明显多 于在随机序列中可能出现的次数，这些片段为优势片段(模式)。可采用以下方法获取优势片段。根据神经元的间隔模式确定间隔范围数(可为 3 6个)，可用滑移法在相邻间隔的范围内寻找神经元放电间隔排序模式，找出重复频率 大于重复次数(如3次)的间隔组合作为优势片段的原始序列集合；产生替代数据序列替 代原始数据序列。将matlab数据进行平均ISI计算，再用平均ISI计算获得的数据同原始 数列比对，将之转化为平均ISI的整数倍数字组合。寻找替代数据序列中的出现频率大于3 次的优势片段作为随机序列集合；比较单元将原始序列集合与随机序列集合比较，如果某 一片段在随机序列集合中的出现率大于在原始序列集合中的出现率，则在随机序列集合中 剔除该片段，在原始序列集合中留下该片段，检验确定不是相邻的片断的重复即为优势片 段。还可采用所有的峰电位间隔与第一个间隔比较，获得优势片段。同时按照片段出 现次数的高低进行排序、计算在原始序列中该片段组合出现的总次数，三联片段(例如111 或346)计算出现率最多的三位(如只选择出现次数最高的前三个片段组合)，四联、五联、 六联全部保留。步骤2.根据神经元峰值间隔计算神经元平均频率、平均间隔。峰值提取单元从原 始放电(神经元放电信号)序列中提取峰电位信号，进行相邻间隔峰值时间计算，计算单元 根据神经元峰点的个数，以及放电记录时间计算平均频率，各个间隔与相应神经元峰点的 个数之积总和除以神经元峰点电位个数为平均间隔；计算神经元放电个数、放电时间、频率 等；根据峰峰间隔计算平均频率、平均间隔，并将原始间隔序列与平均间隔比较转化为峰电 位间隔序列(即0、1序列)；通过原始间隔与平均间隔比较，如果即时间隔大于平均间隔， 将序列中的数据置为1，如果小于平均间隔则将序列中的数据置为0，顺次转化从而将原始 数据转化为0，1数据串。步骤3.比较单元分别比较上述峰电位间隔序列的重复性，得到被分段的数目复 杂度c (η)。将上述替代数据(即0，1数字串)输入一个不相关数字集合后进行重复比对， 检测与替代数据的相似性，在初始化过程中自行通过细胞筛检计算峰电位个数，获得放电 脉冲数η (即数据串的数据个数，可趋于无穷，一般η > 1000)，根据长序列的复杂度公式 Limc (Ii)n^00= b (n) = n/l0g2n确定数目复杂度c (η)，以表明数据的复杂程度。该方法确定 复杂度基于电生理数据的直接应用，不需要做其他数据挑选和时间窗的截取，一次得出数 据。还可应用近似熵测度复杂度等方法，反映时间序列的不规则程度，所需数据点数。步骤4.图形用户界面(GUI)工具绘制ISI时间间隔图。按照峰电位间隔时间分 布图，以放电时间为横坐标，即时对应的相邻两个放电间隔为纵坐标画散点图，从中找出神经元放电是否具有周期性样的节律性，根据系统某一运动参量的当前值与前些时刻状态值 之间如果有确定性的函数关系，这种回归可在ISI回归映射图上具有明确的几何结构，从 而可以确定放电串中有无周期性或者其它特异性，由此绘制ISI时间间隔图。步骤5.根据神经元的种类设定带宽值bin，并以带宽值bin宽为横坐标，相应带宽 值bin宽的个数为纵坐标，绘制ISI非序列间隔直方图。根据神经元的种类设定带宽值bin， 将原始序列标准化为相对ISI值(峰峰间隔在0 1 · bin ms之间为1，在1 · bin ms 2-bin ms之间为2...依此类推)。计算相应带宽值bin下不同的峰电位间隔的个数。再将神经元进行分拣，确保神经元为单一神经元后，将所计算的神经元峰电位的 间隔按照即时时间予以排列，并根据优势片段的筛选要求进行重新归类。优势片段的要求 为根据神经元的种类先设定带宽值bin，比方说海马神经元的bin宽设定为35ms，接下来
峰峰间隔(bin宽的个数)设置遵循以下规则0 35ms为l，35ms 70ms为2，......依
此类推，将ISI数值再次转化为以bin宽为电位的相对ISI值。以带宽值bin宽为横坐标， 相应带宽值bin宽的个数为纵坐标，画直方图，连接神经元放电的稀疏以及是否有不同放 电的特征，用于判定神经元的内在放电特性，比如发电特性是burst还是single spike神 经元。步骤6.绘制ISI回归映射图。以相邻两个放电间隔(η和η+1神经元峰点间隔) 为横坐标，以下一相邻两个放电间隔(η+1和η+2神经元峰点间隔)为纵坐标，绘制ISI回 归映射图，用于发现神经元间隔之间的相关性，并可判定神经元的节律特点。该系统基于Matlab环境开发的，与专用数据分析系统相比，具有更强的通用性和 可扩展性。该系统工作稳定，操作简单，可以实时给出神经元基本生理特性的定量分析结 果，为实验数据的进一步分析和实验方案的改进提供了便利的参考信息。本发明能够在实 验过程及离线后提供定量的神经元放电特征信息，提高神经电生理实验的效率，开发了基 于主成分分析方法的实时神经元分类工具箱，为后期的深入数据分析提供参考。
权利要求
神经元单位放电分析方法，其特征在于，包括下列步骤数据采集单元采集神经元放电信号，并转化为Matlab数据，获取优势片段；峰值选取单元从优势片段中选择神经元峰电位，确定神经元峰 峰间隔；根据神经元峰 峰间隔计算神经元平均间隔；比较单元将神经元原始间隔序列与平均间隔序列比较转化为峰电位间隔序列，比较上述峰电位间隔序列的重复性，得到被分段的数目复杂度c(n)；图形用户界面工具绘制ISI时间间隔图，以放电时间为横坐标，对应的相邻两个放电间隔为纵坐标绘制ISI时间间隔图；根据神经元的种类设定带宽值bin，根据带宽值bin绘制ISI非序列间隔直方图；以相邻两个放电间隔为横坐标，以下一相邻两个放电间隔为纵坐标，绘制ISI回归映射图，由此判定神经元的节律特点。
2.根据权利要求1所述的神经元单位放电分析方法，其特征在于，获取优势片段的方 法为，根据神经元的间隔模式确定间隔范围数，在间隔范围数个相邻间隔内寻找神经元放 电间隔排序模式，找出重复频率大于重复次数的间隔组合作为优势片段的原始序列集合; 产生替代数据序列替代原始数据序列；搜索单元寻找替代数据序列中的优势片段作为随机 序列集合；比较单元比较原始序列集合与随机序列集合，如果某一片段在随机序列集合中 的出现率大于在原始序列集合中的出现率，则在随机序列集合中剔除该片段，留下原始序 列集合中的片段，如果不是相邻的η片断的重复即为优势片段。
3.根据权利要求1所述的神经元单位放电分析方法，其特征在于，通过细胞筛检计算 神经元峰电位个数，获得放电脉冲数η，调用公式Limc (η)…=b (n) = n/log2n确定被分段的数目复杂度c (η)。
4.根据权利要求1所述的神经元单位放电分析方法，其特征在于，所述绘制ISI非序列 间隔直方图的方法为，根据神经元的种类设定带宽值bin，以带宽值bin为横坐标，相应带 宽值bin的个数为纵坐标绘制ISI非序列间隔直方图。
5.一种神经元单位放电分析系统，其特征在于，数据采集单元获取优势片段，峰值选取 单元从中选择神经元峰电位，确定神经元峰_峰间隔，根据神经元峰_峰间隔计算神经元 平均间隔；比较单元将神经元原始间隔序列与平均间隔序列比较转化为峰电位间隔序列， 比较上述峰电位间隔序列的重复性，得到被分段的数目复杂度c(n)；图形用户界面工具绘 制ISI时间间隔图，以放电时间为横坐标，即时对应的相邻两个放电间隔为纵坐标绘制ISI 时间间隔图；根据神经元的种类设定带宽值bin，根据带宽值bin绘制ISI非序列间隔直方 图；以相邻两个放电间隔为横坐标，以下一相邻两个放电间隔为纵坐标，绘制ISI回归映射 图，由此判定神经元的节律特点。
6.根据权利要求5所述神经元单位放电分析系统，其特征在于，数据采集单元获取优 势片段具体为，根据神经元的间隔模式确定间隔范围数，在间隔范围数个相邻间隔内寻找 神经元放电间隔排序模式，找出重复频率大于重复次数的间隔组合作为优势片段的原始序 列集合，产生替代数据序列替代原始数据序列；搜索单元寻找替代数据序列中的优势片段 作为随机序列集合；如果某一片段在随机序列集合中的出现率大于在原始序列集合中的出 现率，则在随机序列集合中剔除该片段，留下原始序列集合中的片段，如果不是相邻的η片 断的重复即为优势片段。
7.根据权利要求5所述神经元单位放电分析系统，其特征在于，图形用户界面工具根 据神经元的种类设定带宽值bin，以带宽值bin为横坐标，相应带宽值bin的个数为纵坐标绘制ISI非序列间隔直方图。
全文摘要
本发明请求保护一种单位放电分析系统和分析方法，涉及信息处理技术领域，峰值选取单元从优势片段中选择神经元峰电位，确定神经元峰-峰间隔；通过原始数据的转化，提取相邻动作电位峰电位的时间间隔，进行数据转化，根据不同要求分别得出其统计值频率、复杂度、优势片段、以及绘制出神经节律特点多分析图，判定神经元的节律特点。
文档编号G06F19/00GK101908099SQ20101023304
公开日2010年12月8日 申请日期2010年7月22日 优先权日2010年7月22日
发明者伍亚民, 刘良明, 朱志茹, 隋建峰, 高洁 申请人:中国人民解放军第三军医大学第三附属医院