CN106062669A - 用于脑机接口的系统 - Google Patents
用于脑机接口的系统 Download PDFInfo
- Publication number
- CN106062669A CN106062669A CN201580006550.7A CN201580006550A CN106062669A CN 106062669 A CN106062669 A CN 106062669A CN 201580006550 A CN201580006550 A CN 201580006550A CN 106062669 A CN106062669 A CN 106062669A
- Authority
- CN
- China
- Prior art keywords
- spike
- unit
- template
- module
- template matching
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000036982 action potential Effects 0.000 claims abstract description 78
- 210000002569 neuron Anatomy 0.000 claims abstract description 30
- 238000001914 filtration Methods 0.000 claims abstract description 27
- 238000000605 extraction Methods 0.000 claims abstract description 10
- 238000006243 chemical reaction Methods 0.000 claims abstract description 9
- 238000001514 detection method Methods 0.000 claims description 68
- 239000007943 implant Substances 0.000 claims description 13
- 238000010276 construction Methods 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 11
- 238000004422 calculation algorithm Methods 0.000 claims description 9
- 230000000694 effects Effects 0.000 abstract description 29
- 238000012545 processing Methods 0.000 abstract description 14
- 238000000034 method Methods 0.000 abstract description 11
- 238000013459 approach Methods 0.000 abstract description 8
- 230000008901 benefit Effects 0.000 abstract description 4
- 238000012544 monitoring process Methods 0.000 abstract description 3
- 238000011160 research Methods 0.000 abstract description 3
- 238000004891 communication Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 9
- 230000001537 neural effect Effects 0.000 description 9
- 210000004556 brain Anatomy 0.000 description 8
- 230000009467 reduction Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 241001465754 Metazoa Species 0.000 description 4
- 238000012546 transfer Methods 0.000 description 4
- 206010010904 Convulsion Diseases 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 230000007704 transition Effects 0.000 description 3
- 230000001364 causal effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 230000000638 stimulation Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000001225 therapeutic effect Effects 0.000 description 2
- 206010033799 Paralysis Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000003710 cerebral cortex Anatomy 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 206010015037 epilepsy Diseases 0.000 description 1
- 230000001037 epileptic effect Effects 0.000 description 1
- 208000028329 epileptic seizure Diseases 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 238000013021 overheating Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 210000000578 peripheral nerve Anatomy 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 210000000278 spinal cord Anatomy 0.000 description 1
- 238000009966 trimming Methods 0.000 description 1
- 230000002747 voluntary effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/40—Detecting, measuring or recording for evaluating the nervous system
- A61B5/4058—Detecting, measuring or recording for evaluating the nervous system for evaluating the central nervous system
- A61B5/4064—Evaluating the brain
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F3/00—Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
- G06F3/01—Input arrangements or combined input and output arrangements for interaction between user and computer
- G06F3/011—Arrangements for interaction with the human body, e.g. for user immersion in virtual reality
- G06F3/015—Input arrangements based on nervous system activity detection, e.g. brain waves [EEG] detection, electromyograms [EMG] detection, electrodermal response detection
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Neurology (AREA)
- Biomedical Technology (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Neurosurgery (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Dermatology (AREA)
- General Physics & Mathematics (AREA)
- Human Computer Interaction (AREA)
- Physiology (AREA)
- Psychology (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- User Interface Of Digital Computer (AREA)
Abstract
本发明提供一种两步法方案来提供BCI系统。在第一步骤中,本发明提供低功率植入平台以用于放大并过滤胞外记录、执行模数转换(ADC)并实时检测动作电位,该植入平台连接至远程装置,该远程装置能够执行特征提取和锋电位分类的处理器密集型任务,由此为每个神经元生成多个预定模版以用于第二处理步骤。在第二步骤中,该低功率植入平台放大并过滤胞外记录、执行ADC并且检测动作电位,动作电位可在片上被匹配至由外部接收器在第一步骤中生成的预定模版。此两步法方案利用离线处理和在线处理二者的优势,提供了执行单个神经元活动的多个记录的有效并且安全的方法,以用于研究或监测应用或者用于远程装置的控制。
Description
背景技术
神经系统学可以记录由神经元生成的胞外动作电位或锋电位,从而理解信息是如何通过神经元系统呈现和传输的。
直到最近,这些实验通常涉及在数小时的短周期内采样少量神经元(例如,少于15),但随着延久电极阵列的改进,记录数月内来自大量神经元(例如,多于15)的锋电位已成为可能。这些技术激发了脑机接口(BCI)的发展,其直接与神经元系统通信以用于治疗效果或者大脑运行的研究学习。BCI可实现动物或病患大脑与任意外部计算机、数据处理系统或相连接的任意电子装置之间的通信。其中一个示例是使用来自瘫痪病人大脑皮层的神经信号来操作比如机器人假肢的辅助装置的脑机界面(BMI)。
为了获得精确的结果,脑机接口或者用于记录单个神经元活动的任意其他系统应借助于胞外记录来检测动作电位,并且识别哪些动作电位对应于哪些神经元。基于每个神经元生成具有模式化(stereotyped)波形的动作电位这一原理,根据动作电位的形状,对动作电位进行分类。数据被转发给外部计算机,在外部计算机中基于若干独立神经元的动作处理信息或命令。
已提出了具有植入大脑的锋电位分类模块的BCI,锋电位分类模块用于根据模式化的波形的形状对动作电位分类,但这涉及复杂的难以在片上实施的运算,因此当处理来自许多记录电极的数据时,其会消耗过多的功率,这可能导致过热并且还可能会限制植入装置的小型化程度。因此,这对于接受神经元监测的对象的健康而言是异常危险的。为了解决这一问题,设计者们已尝试尽可能地从植入模块中移除处理、将植入模块限于感测技术并且在外部装置中执行所有的信号处理。
虽然这简化了植入装置,但给植入装置和外部处理电路之间的通信开销增加了很大负担。所监测的信号必须全部被传输给外部处理电路,这消耗大量带宽,其中很多带宽是通过传输比如背景噪声和不包含动作电位的信号周期而浪费的。这又延迟了对于治疗应用而言可能很关键的锋电位检测和分类进程。
发明内容
本发明的实施例采取混合方案,其中一些处理功能在外部装置中执行,但一些是在植入装置中以本地方式执行。这减少了通信开销,因为锋电位检测进程和模版匹配进程是在植入装置中以本地方式执行的,并且只需要将比较的结果或者先前未检测到的动作电位传输给外部装置。因此,相较于传统系统,本发明提高了速度而并没有使植入装置复杂化。
根据本发明的一方面,提供了根据权利要求1的用于脑机接口的系统。
根据本发明的一方面,提供了根据权利要求28的植入装置。
根据本发明的一方面,提供了根据权利要求29的远程装置。
可选的特征在从属权利要求中列出。
本发明的实施例提供了两步法方案来提供无线BCI系统。在第一步骤中,低功率植入平台放大并过滤胞外记录,执行模数转换(ADC)并实时检测动作电位。处理器以无线方式连接至能够执行特征提取和锋电位分类的处理器密集型任务的外部接收器,由此为每个神经元生成多个预定模版以用于第二处理步骤。在第二步骤中,低功率植入平台放大并过滤胞外记录,执行ADC并且检测动作电位,动作电位可在片上被匹配至由外部接收器在第一步骤中生成的预定模版。此两步法方案利用了离线处理和在线处理二者的优势,提供了有效并且安全的、用于执行单个神经元活动的多个无线记录的方法。各实施例可用于执行多个单胞记录(single-cell recording),用于研究或监测应用,或者用于控制任意外部电子装置。
附图说明
将参照以下附图以示例方式描述本发明。
图1是本发明的实施例的无线脑机接口的示意图;
图2是本发明的实施例的无线脑机接口的放大示意图;
图3是图1中锋电位检测模块的示意图;
图4是示出了滤波器阶数对锋电位检测性能的影响的图表;
图5是示出了低通截止频率对锋电位检测性能的影响的图表;
图6是示出了高通截止频率对锋电位检测性能的影响的图表;
图7a是示出了ADC采样率对锋电位检测性能的影响的图表;
图7b是示出了ADC采样率对模版构建性能的影响的图表;
图8a是示出了ADC分辨率对锋电位检测性能的影响的图表;
图8b是示出了ADC分辨率对模版构建性能的影响的图表;
图9是示出了根据三种小波系数的动作电位的分布的图表;
图10a是示出了变温相似变量(temperature-analogous variable)对动作电位的群集大小的影响的图表;
图10b是图9a所示动作电位的群集的波迹;
图10c是由图9a所示动作电位的群集生成的多个锋电位模版的波迹;
图11是示出了各种各样的距离度量(distance metric)所需的性能和相对处理器时间的图表;
图12是示出了波峰对准(peak-alignment)误差对模版匹配性能的影响的图表;
图13是示出了窗口尺寸对模版匹配性能的影响的图表。
具体实施方式
图1示出了根据本发明的实施例的一种用于无线脑机接口(BCI)的系统1。BCI包括锋电位检测模块10、模版构建模块20和模版匹配模块30。模版构建模块20以无线方式与锋电位检测模块10和模版匹配模块30通信。锋电位检测模块10可与模版匹配模块30通信。锋电位检测模块10和模版匹配模块30通过有线连接直接通信,并且可在单片上实施。
锋电位检测模块10和模版匹配模块30设置在装置100中,装置100植入大脑中或大脑附近并且联接至用于记录大脑中至少一个神经元处或附近的信号的一个或多个电极。模版构建模块20设置在远程装置200中,远程装置200与植入装置100无线连接。远程装置200可相对于患者在外部使用,或者可植入在与大脑分离的位置处。模版构建模块20通过无线连接以无线方式与锋电位检测模块10和模版构建模块30通信。模版匹配模块30可以以无线方式传输模版匹配进程的结果。
可替换地,系统可通过模版构建模块20和锋电位检测10以及模版匹配模块30之间的电缆通信实现。此实施例不提供无线通信的优势,但仍提供用于在线片上锋电位检测和分类的有效方法。仍可通过减少通过电缆连接的数据传输率来降低植入装置100的功耗。在有线的实施例中,远程装置200可植入身体的其它部分。远程装置200可在不需要时断电,从而进一步降低功耗。
因此,本实施例提供了一种双部件系统,其包括用于检测锋电位并将其与模版相匹配的低功率植入平台100,低功率植入平台100以无线方式与外部接收器200连接,外部接收器200可执行处理器密集型模版构建。
锋电位检测模块10检测由大脑中神经元生成的动作电位,动作电位是通过使用本领域已知的一个或多个电极或者微触压感技术记录的。一个或多个电极中的每个电极记录单个位置,或者可记录多个位置。对于每个记录位置,锋电位检测模块10可检测由靠近该记录位置的多个神经元生成的多个动作电位。
锋电位的形状是由检测到的动作电位的波迹形成的形状。模版构建模块20生成锋电位模版,锋电位模版对应于由具体神经元生成的动作电位的模式化形状。模版构建模块20基于锋电位检测模块10检测到的一个或多个动作电位、为每个模式化锋电位形状生成锋电位模版,以使得可基于将检测到的动作电位与生成的锋电位模版进行匹配而识别未来具有相同或类似形状的动作电位。模版构建模块20可基于检测到的由单个神经元生成的多个动作电位的平均值生成锋电位膜板。模版构建模块20可基于由多个相异的神经元生成的多个相异的动作电位生成多个相异的锋电位模版。模版构建模块20存储所生成的锋电位模版,并且在需要时将模版提供给模版匹配模块30。
模版匹配模块30被设置为根据动作电位的形状群集由不同的神经元生成的动作电位。模版匹配模块30将锋电位检测模块10检测到的动作电位与模版构建模块20生成的一个或多个锋电位模版进行比较。模版匹配模块30根据比较的结果输出信号。输出信号可能表示动作电位与其中一个锋电位模版之间的正匹配。模版匹配模块30可将检测到的一个或多个动作电位与模版构建模块20生成的多个预定锋电位模版进行比较。输出信号可能表示一个或多个动作电位与多个预定锋电位模版之间的一个或多个正匹配。模版匹配模块30可以以无线方式将输出信号传输给外部装置。
在某些情况下,如果出现否定的比较结果,则模版匹配模块30可将检测到的动作电位传输给模版构建模块20,以生成未知的锋电位形状的新模板,其可能代表先前未被设备检测到的明显的神经元活动。然而,否定的比较结果还可能是因为检测到的信号根本不是动作电位,在这种情况下,不期望激活模版构建模块20,因为在这些情况下生成的任何模版都是不需要的。模版构建模块20可基于多个未知的锋电位形状之间的相似性确定是否生成新模版。模版构建模块20能够区分未知的动作电位和不属于动作电位的噪声信号,因为多个噪声信号通常不具有彼此类似的形状。
模版匹配模块30可具有用于缓存从模版构建模块20接收到的、生成的模版的本地存储。可根据用户需要来选择存储的大小。如果期望的神经元活动代表相对少量的相异动作电位类型,则模版匹配模块30处的模版的本地存储是改善本发明的效率的有效途径,因为植入平台可以无监督地有效操作。在更复杂的情况下,模版匹配模块的本地存储能力可能不足以存储所有需要的模版,在这种情况下,与模版构建模块20的定期通信更有效。
在两种情况下可能都需要周期性更新模版以补偿记录条件的变化,这些变化可能对动作电位存在整体影响,比如背景噪声的变化,或者动作电位形状的转变。模版构建模块20可通过取对应于模版的、最新检测到的动作电位的形状的平均值来执行模版的周期性更新。
允许植入装置100无监督地操作的能力部分地取决于锋电位分类和模版匹配进程的优化性质,将在下文中对其进行详细解释。这又使得本发明能够将植入装置100和外部模版构建模块20的功能分开。此外,通过避免大功率和/或大带宽需求,提供了在可记录的频道数量方面的可扩展性。
图2示出了锋电位检测模块10,其包括信号滤波单元11、模数转换器(ADC)15和锋电位检测单元16。信号滤波单元11被配置为使用带通滤波器13过滤大脑内至少一个神经元处或附近记录的信号。胞外记录可包括由靠近的多个神经元生成的一个或多个动作电位。信号滤波单元11应用滤波器13,滤波器13被配置为增加胞外记录的信噪比(SNR)。
在一实施例中,信号滤波单元11应用椭圆滤波器、巴特沃斯滤波器和贝塞尔滤波器中的一种滤波器。可替换地,信号滤波单元11可应用任意其他滤波实施方式。
信号滤波单元11被配置为实时地过滤胞外记录。由信号滤波单元11应用的滤波器13可为因果滤波器,其基于当前讨论的数据点之前的数据点、而不是该数据点之后的那些数据点、计算每个数据点。
如图2所示,ADC 15可位于锋电位检测单元16之前,以使得在数字域中执行锋电位检测。然而,在另一实施例中,ADC 15还可位于锋电位检测单元16之后,以使得在模拟域中执行锋电位检测。
图3是锋电位检测模块10的示意图,其示出信号滤波单元11包括前置放大器12、滤波器13以及可选的增益放大器14。前置放大器12被配置为在应用过滤之前、最大化胞外记录的SNR。在一实施例中,前置放大器12是低噪声低功率前置放大器。增益放大器14可为可编程增益放大器,配置为最大化对ADC 15的输入动态范围。在通过直接对过滤后的数据进行模版匹配而完成检测的情况下,锋电位检测单元16可以是可选的。
在一实施例中,锋电位检测模块10可实施为集成电路。集成电路包括片上的信号滤波单元11、ADC 15和锋电位检测单元16。可减少从原始数据至锋电位包的片外通信,从而显著降低片外通信的数据传输率。片外数据通信的功耗是母线电容、电源电压和数据传输率的函数。通过传输锋电位包,显著地降低了片外数据传输率,减少了系统功耗。
图4示出了在使用椭圆滤波器时,对于具有高、中和低SNR的信号,滤波器13的阶数对锋电位检测单元16性能的影响。信号滤波单元11应用的滤波器13可为1阶、2阶、3阶和4阶滤波器中的一种。较低阶的滤波器引入低水平的相位失真。
减少通带带宽可增加胞外记录的SNR。过窄的通带可降低锋电位检测单元16的性能。在一实施例中,基于合成试验信号的处理、使用本领域已知技术测得的锋电位检测单元16的性能相对于精确检测到的神经元活动的百分比应维持在近似90%或更高。
图5示出了对于具有高、中和低SNR的信号,低通截止频率对锋电位检测单元16的性能的影响。选择滤波器13的低通截止频率以使得锋电位检测单元16的性能最大化。在一实施例中,低通截止频率是3kHz。
图6示出了对于具有高、中和低SNR的信号,高通截止频率对锋电位检测单元16的性能的影响。选择滤波器13的高通截止频率以使得锋电位检测单元16的性能最大化。在本发明一实施例中,高通截止频率是100kHz、200kHz和300kHz中的一个。
ADC 15被配置为将至少一个神经元处或附近记录的信号从模拟信号转换为数字信号。在一实施例中,ADC 15转换由信号滤波单元11输出的过滤后的信号。可替换地,ADC15将胞外记录转换为数字信号,然后通过信号滤波单元11过滤该数字信号。信号滤波单元在数字域中应用的滤波器13可为非因果滤波器,其基于当前讨论的数据点之前的数据点以及当前数据点之后的那些数据点计算该数据点。
可通过调整ADC 15的采样率和分辨率改变胞外记录的数据传输率。减小胞外记录的数据传输率可降低锋电位检测模块10中的功耗水平。过低的采样率或分辨率可降低锋电位检测单元16的性能和模版构建模块20的性能。锋电位检测单元16的性能应维持在80%或更高。如果采样率和ADC分辨率的变化(相较于滤波的变化)是优化锋电位检测的主要参数,则功率代价使得可获得的性能可能小于将滤波作为主要变化参数的情况。然而,通过控制滤波、采样率和ADC分辨率中的每一个,锋电位检测单元16的性能可为90%或更高。相对于由人工合成信号生成模版,信号模版构建模块20的性能应维持在50%或更高。
图7a和图7b示出了对于具有高、中和低SNR的信号,ADC采样率分别对锋电位检测单元16和模版构建模块20的性能的影响。在本发明的一实施例中,ADC 15的采样率大于7kHz。
图8a和图8b示出了对于具有高、中和低SNR的信号,ADC 15的分辨率分别对锋电位检测单元16和模版构建模块20的性能的影响。ADC 15的分辨率在6-12比特的范围内。在一实施例中,ADC 15的分辨率是10比特。
锋电位检测单元16被配置为检测包含在已由信号滤波单元11过滤并且已由ADC15转换过的胞外记录内的一个或多个动作电位。锋电位检测单元16通过将胞外记录的幅度与幅度阈值进行比较而检测该一个或多个动作电位。通过锋电位检测单元16自动选择幅度阈值,但其也可被手动设定,因为锋电位检测单元16的无线配置是作为基于试验或误差的校正步骤以用于实现锋电位检测的预期速率。可替换地,原始数据的一小部分可被传输给远程装置200以用于计算幅度阈值。
在一实施例中,锋电位检测单元16预估胞外记录中背景噪声的标准偏差。锋电位检测单元16可基于预估的标准偏差选择幅度阈值。幅度阈值可为预估的标准偏差的3倍至5倍。
锋电位检测单元16可基于胞外记录的绝对值的中值预估背景噪声的标准偏差。
在一实施例中,幅度阈值Thr被自动设为:
Thr=4σn,
其中,而x是信号。
锋电位检测模块10可包括复用器17,其将来自锋电位检测单元16的多个胞外记录的输出进行组合。来自复用器17的组合输出被发送给模版构建模块20或模版匹配模块30。来自锋电位检测单元16的输出由多个锋电位包构成,每个锋电位包包括检测到的动作电位。可替换地,多个锋电位包可通过串行总线而不是该复用器17输出。锋电位包可按顺序被传输,并且每个锋电位包具有通道ID和时间戳。
模版构建模块20包括特征提取单元21和锋电位分类单元22。特征提取单元21被配置为应用小波变换,以计算检测到的动作电位的多个小波系数。小波变换是信号的时频表示,由信号和小波集合之间的卷积定义。小波是唯一小波函数的扩张(或收敛)型和扰动型。小波系数是通过计算动作电位与小波函数在不同时间移位和扩张(或收敛)比例处的卷积获得的。因此,小波变换仅将为时间函数的信号映射到独立的扩张和移位变量的函数上。小波函数的收敛型映射信号的高频分量,而扩张型映射低频分量。
可替换地,特征提取单元21可被配置为实施主成分分析,以计算检测到的动作电位的多个主成分,或者可使用比如波峰幅度、宽度或峰谷能量之类的锋电位特征。
特征提取单元21从计算出的多个小波系数中选择小波系数优选集。小波系数优选集包括一个或多个系数,其表示多个动作电位之间的区别。特征提取单元21基于计算出的正态偏差(deviation from normality)选择小波系数优选集。特征提取单元21可对计算出的小波系数应用柯尔莫诺夫-斯米尔诺夫正态检验。小波系数优选集可包括具有最低正态分布的一个或多个系数。在一实施例中,小波系数优选集包括选定的10个小波系数。
图9示出了由3个优选的小波系数限定的三维空间中多个动作电位的分布。
锋电位分类单元21从群集中多个动作电位所共享的一个或多个主特征生成锋电位模版,动作电位的群集具有优选的具有相似值的小波系数。动作电位的群集可能位于由小波系数优选集限定的三维空间内的极近处。锋电位分类单元22被配置为通过应用超顺磁聚类(SPC)算法检测动作电位的群集。锋电位分类单元22基于小波系数优选集应用超顺磁聚类算法。
可替换地,锋电位分类单元22可利用K-均值分类器、期望最大化算法、群集的手动界定或其他聚类算法来实施。
SPC算法模拟每个数据点与其最近邻点(nearest neighbour)之间的交互,其中,模拟的交互强度是以增加的间距呈指数跌落的。在使用波特模型的物理模拟中,每个数据点被分配一个随机的初始态,并且极近(即,对应于群集)的数据点将一起改变状态。最近邻点一起改变状态的概率取决于波特模型中的变温相似变量。
高温对应于近邻点一起改变状态的低概率。低温对应于近邻点一起改变状态的高概率,即使是对于较弱的交互强度也是如此。在某些中等的温度范围内存在“超顺磁阶段”,其中,仅那些极近的点将一起改变其状态。
可基于群集大小标准(cluster size criterion)检测动作电位的群集。SPC算法计算:随着变温相似变量的变化、多个群集在尺寸上的变化。在类同于超顺磁阶段的范围内,变温相似变量的增加产生了具有相对大量成员的群集。如果满足群集大小标准,则变温相似变量可为固定的,并且可从检测到的一个或多个动作电位的群集中选择一个或多个主特征。可替换地,可从不同温度中选择具有相对大量成员的群集。
图10a示出了应用于图9中多个动作电位的SPC算法的结果。群集大小由变温相似变量的函数表示,垂直虚线表示进入超顺磁阶段的过渡点。用于生成锋电位模版的该一个或多个主特征是从满足群集大小标准的动作电位的群集中选择的。图10b示出了叠加的多个动作电位。重群集区域可与未群集的外围动作电位一起被识别。图10c示出了从图10c的动作电位的群集的主特征生成的多个锋电位模版。
模版匹配模块30被设置为将检测到的动作电位与锋电位模版进行比较。模版匹配模块30可接收来自锋电位检测模块10的复用器17或串行总线的输出,输出包括多个锋电位包,每个锋电位包包括检测到的动作电位。在一些实施例中,来自锋电位检测模块10的包输出可允许模版匹配模块30无需输入解复用器或输入缓冲器(配置为存储每个通道的检测到的动作电位)即可操作,因为使用可由通道ID和时间戳构成,从而使锋电位的模版匹配可在具体通道中执行。
模版匹配模块30可输出信号,以基于检测动作电位与锋电位模版之间的距离指示正匹配。模版匹配模块30使用距离度量将检测到的动作电位与锋电位模版进行比较,距离度量比如为平方欧几里得距离度量、1-范数距离度量、无穷范数距离度量、马氏距离度量、最近邻距离度量或其他任意适当的距离度量。
图11示出了使用各种各样适当的距离度量的模版匹配模块30的性能和所需的相对处理器时间。在一实施例中,距离度量可为欧几里得距离度量d。欧几里得距离度量是由以下公式计算出的,其中yi是检测到的动作电位的成分,Ti是锋电位模版的成分:
可替换地,模版匹配模块30可使用1-范数距离度量,其计算为:
可选择模版匹配模块30所使用的距离度量,从而最小化模版匹配模块30中的功耗水平或者处理器时间使用,并且同时维持最低性能水平。
模版匹配模块30包括对准单元31,其被配置为将检测到的动作电位与锋电位模版基本对准。对准单元31可将动作电位与锋电位模版对准以使其落入预定波峰对准误差内,但也可基于信号跃迁执行对准,比如上升沿和幅度阈值。
图12示出了对于各种各样适当的距离度量,波峰对准误差对模版匹配模块30性能的影响。不准确的波峰对准可能引起模版匹配模块30性能的降低,而波峰对准误差的减小可引起模版匹配模块30中功耗水平的增加。模版匹配模块30的性能应维持在80%或更高。对准单元31被配置为具有少于约0.125ms的波峰对准误差。
模版匹配模块30被配置为通过修整动作电位而选择窗口尺寸。窗口尺寸的减小会导致模版匹配模块30中存储的数据量的减少以及计算适当距离度量的复杂度的降低。
图13示出了对于各种各样适当的距离度量,窗口尺寸对模版匹配模块30性能的影响。减小的窗口尺寸会导致模版匹配模块30性能的降低。模版匹配模块30的性能应维持在70%或更高。模版匹配模块30被配置为使用如约0.5ms这么小的窗口尺寸。
可使用处理单元实施模版匹配模块30。与使用比如现场可编程门阵列(FPGA)的可配置逻辑装置相比,模版匹配模块30的基于处理器的实施方式可提供更有效的操作。
本发明提供两步法方案来提供无线BCI系统。两步法方案利用离线处理和在线处理二者的优势,提供有效并且安全的神经元研究方法。无线BCI系统允许根据神经元的波形识别具体神经元冲动(neurons firing),并且以经过分类的、便于使用的形式传输神经元活动。
本发明的实施例可用于控制假肢器官装置。例如,由本发明传输的信号可控制一个或多个义肢。可根据检测到的、用于负责失去的四肢的日常活动的、一个或多个神经元的活动来控制义肢。可根据通常用于负责控制身体相关部分的神经元的活动来控制其他假肢器官。
可替换地,本发明的实施例可用于控制传送给电机结构的电刺激,从而在病患的随意控制下产生四肢运动,电机结构包括脊髓、末梢神经或肌肉。
可替换地,本发明的实施例可用于监测一定时间段内一个或多个神经元的活动。例如,本发明可用于记录动物的神经元活动。本发明的低功率无线BCI允许动物在记录过程中自由运动,从而既提供了更多种类的可执行实验、又提供了动物福利方面的改进。
本发明的实施例还可用于为其他临床应用监测神经元活动。对于癫痫病人,本发明可检测表示癫痫发作的一个或多个神经元的活动。医疗装置可被设置为通过使用本发明的无线BCI监测一个或多个神经元的活动来预测可能的癫痫发作,并且可被配置为传送用于防止预测的癫痫发作的电刺激。
本领域技术人员应意识到,可对上述实施例做出多种修改,这些修改落入随附权利要求书限定的本发明的范围内。可将各种互相兼容的实施例的特征进行组合。例如,所描述的任一具体滤波器形状可与上面列出的任一截止频率相兼容,而滤波器的组合又与所描述的任一聚类算法相兼容。根据上述教导,本领域技术人员将能意识到在实现本发明的效果的同时、这些特征可以互相交换的范围。
Claims (29)
1.一种用于脑机接口的系统,包括:
锋电位检测模块,被设置为检测由一个或多个神经元生成的至少一个动作电位;
模版匹配模块,被设置为将检测到的所述至少一个动作电位与一个或多个预定锋电位模版进行比较;以及
模版构建模块,被配置为生成所述一个或多个预定锋电位模版,并且被设置为远离所述锋电位检测模块和所述模版匹配模块。
2.根据权利要求1所述的系统,其中所述锋电位检测模块和所述模版匹配模块被设置在植入装置中,所述植入装置被配置为传输检测到的所述至少一个动作电位并且接收所述一个或多个预定锋电位模版;并且
所述模版构建模块被设置在远程装置中,所述远程装置被配置为接收被传输的所述至少一个动作电位并且传输所述一个或多个预定锋电位模版。
3.根据权利要求1或权利要求2所述的系统,其中所述模版构建模块被设置为与所述锋电位检测模块和所述模版匹配模块无线通信。
4.根据权利要求1-3中任一项所述的系统,其中所述锋电位检测模块包括:
信号滤波单元,被配置为过滤一个或多个胞外记录;
模数转换单元,被配置为转换所述一个或多个胞外记录;以及
锋电位检测单元,被配置为检测所述一个或多个胞外记录内的所述至少一个动作电位。
5.根据权利要求4所述的系统,其中所述信号滤波单元被配置为使用椭圆滤波器、巴特沃斯滤波器和贝塞尔滤波器中的一种滤波器过滤所述一个或多个胞外记录。
6.根据权利要求4或权利要求5所述的系统,其中所述信号滤波单元被配置为使用2阶滤波器过滤所述一个或多个胞外记录。
7.根据权利要求4-6中任一项所述的系统,其中所述信号滤波单元的高通滤波器截止频率被选为使得所述锋电位检测单元的性能值最大化。
8.根据权利要求7所述的系统,其中所述信号滤波单元的高通滤波器截止频率为约300Hz。
9.根据权利要求4-8中任一项所述的系统,其中所述信号滤波单元的低通滤波器截止频率被选为使得所述锋电位检测单元的性能值最大化。
10.根据权利要求9所述的系统,其中所述信号滤波单元的低通滤波器截止频率为约3kHz。
11.根据权利要求4-10中任一项所述的系统,其中所述模数转换单元的采样率被最小化,使得所述锋电位检测单元的性能保持在90%以上。
12.根据权利要求11所述的系统,其中所述模数转换单元的采样率被减小至约7kHz。
13.根据权利要求4-12中任一项所述的系统,其中所述模数转换单元的分辨率被最小化,使得所述锋电位检测单元的性能保持在90%以上。
14.根据权利要求13所述的系统,其中所述模数转换单元的分辨率为6比特。
15.根据权利要求13所述的系统,其中所述模数转换单元的分辨率被最小化,使得所述模版构建单元的性能保持在50%以上。
16.根据权利要求15所述的系统,其中所述模数转换单元的分辨率为10比特。
17.根据权利要求4-16中任一项所述的系统,其中所述锋电位检测单元被配置为基于所述胞外记录中背景噪声的标准偏差的预估,自动选择幅度阈值。
18.根据权利要求4-16中任一项所述的系统,包括用于手动选择幅度阈值的装置。
19.根据权利要求18所述的系统,其中所述背景噪声的标准偏差是基于过滤后的胞外记录的绝对值的中值而预估的。
20.根据前述权利要求中任一项所述的系统,其中所述模版构建模块包括:
特征提取单元,被配置为计算检测到的所述至少一个动作电位中每个动作电位的特征集;以及
锋电位分类单元,被配置为从具有相似特征的动作电位的群集所共享的一个或多个主特征生成一个或多个预定模版。
21.根据权利要求20所述的系统,其中所述特征提取单元被配置为计算包括多个小波系数的小波变换,并且基于计算出的正态偏差,从所述多个小波系数中选择优选集;并且
所述锋电位分类单元被配置为基于小波系数的所述优选集,通过应用超顺磁聚类算法来检测所述动作电位的群集。
22.根据权利要求20或权利要求21所述的系统,其中所述锋电位分类单元被配置为基于群集大小标准检测用于生成所述一个或多个预定模版的所述动作电位的群集。
23.根据前述权利要求中任一项所述的系统,其中所述模版匹配模块被设置为使用平方欧几里得距离度量、1-范数距离度量、无穷范数距离度量、马氏距离度量和最近邻距离度量中的一种距离度量,将检测到的所述至少一个动作电位与所述一个或多个预定锋电位模版进行比较。
24.根据权利要求23所述的系统,其中所述模版匹配模块被配置为使用所述平方欧几里得距离度量或所述1-范数距离度量。
25.根据前述权利要求中任一项所述的系统,其中所述模版匹配模块包括波峰对准单元,所述波峰对准单元被配置为将检测到的所述至少一个动作电位与所述一个或多个预定锋电位模版基本对准,以落入预定波峰对准误差内。
26.根据权利要求25所述的系统,其中所述波峰对准单元被配置为具有小于约0.125ms的波峰对准误差。
27.根据前述权利要求中任一项所述的系统,其中所述模版匹配模块被配置为使用约0.5ms的窗口尺寸。
28.一种用于脑机接口的植入装置,被配置为根据权利要求2所述的系统的一部分。
29.一种用于脑机接口的远程装置,被配置为根据权利要求2所述的系统的一部分。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GBGB1401613.3A GB201401613D0 (en) | 2014-01-30 | 2014-01-30 | System for a brain-computer interface |
GB1401613.3 | 2014-01-30 | ||
PCT/GB2015/050217 WO2015114347A1 (en) | 2014-01-30 | 2015-01-29 | System for a brain-computer interface |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106062669A true CN106062669A (zh) | 2016-10-26 |
CN106062669B CN106062669B (zh) | 2020-03-20 |
Family
ID=50344124
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201580006550.7A Expired - Fee Related CN106062669B (zh) | 2014-01-30 | 2015-01-29 | 用于脑机接口的系统 |
Country Status (6)
Country | Link |
---|---|
US (1) | US10820816B2 (zh) |
EP (1) | EP3100138B1 (zh) |
JP (1) | JP6617108B2 (zh) |
CN (1) | CN106062669B (zh) |
GB (1) | GB201401613D0 (zh) |
WO (1) | WO2015114347A1 (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109376651A (zh) * | 2018-10-24 | 2019-02-22 | 广东工业大学 | 一种基于cuda架构的gpu加速锋电位分类的系统 |
CN109460785A (zh) * | 2018-10-24 | 2019-03-12 | 广东工业大学 | 一种基于gpu加速的锋电位实时分类的方法 |
CN114925734A (zh) * | 2022-07-20 | 2022-08-19 | 浙江大学 | 一种基于神经拟态计算的在线神经元分类方法 |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11504038B2 (en) | 2016-02-12 | 2022-11-22 | Newton Howard | Early detection of neurodegenerative disease |
US11957897B2 (en) | 2016-04-22 | 2024-04-16 | Newton Howard | Biological co-processor (BCP) |
JP7186947B2 (ja) * | 2017-02-10 | 2022-12-12 | 日本光電工業株式会社 | 体内装置からの通信データ量を変更可能なブレインマシンインターフェースシステム、およびその制御方法 |
US20210263589A1 (en) * | 2017-05-26 | 2021-08-26 | Newton Howard | Kinetic intelligent wireless implant/neurons on augmented human |
US10694964B2 (en) | 2017-09-05 | 2020-06-30 | International Business Machines Corporation | Neural spike scanning for high-density implantable neural recording systems |
CN112741598A (zh) * | 2019-10-30 | 2021-05-04 | 株式会社理光 | 波形生成识别方法和计算机可读介质 |
CN113935380B (zh) * | 2021-10-22 | 2022-05-31 | 北京理工大学 | 一种基于模板匹配的自适应运动想象脑机接口方法与系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102096468A (zh) * | 2011-01-20 | 2011-06-15 | 中山大学 | 一种基于脑-机接口的家电遥控装置及方法 |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7136695B2 (en) * | 2001-10-12 | 2006-11-14 | Pless Benjamin D | Patient-specific template development for neurological event detection |
US7187968B2 (en) * | 2003-10-23 | 2007-03-06 | Duke University | Apparatus for acquiring and transmitting neural signals and related methods |
WO2006003662A2 (en) * | 2004-07-06 | 2006-01-12 | Technion Research & Development Foundation Ltd. | Low power hardware algorithms and architectures for spike sorting and detection |
US7957793B2 (en) * | 2004-12-22 | 2011-06-07 | Wisconsin Alumni Research Foundation | Methods for identifying neuronal spikes |
US20060293578A1 (en) * | 2005-02-03 | 2006-12-28 | Rennaker Robert L Ii | Brian machine interface device |
JP4887482B2 (ja) * | 2005-02-28 | 2012-02-29 | 国立大学法人東北大学 | 神経インプラント装置 |
US20090220425A1 (en) * | 2006-01-19 | 2009-09-03 | Drexel University | Method to Quantitatively Measure Effect of Psychotropic Drugs on Sensory Discrimination |
WO2008042900A2 (en) * | 2006-10-02 | 2008-04-10 | University Of Florida Research Foundation, Inc. | Pulse-based feature extraction for neural recordings |
ES2550960T3 (es) * | 2007-02-01 | 2015-11-13 | Boston Scientific Neuromodulation Corporation | Sistema de neuroestimulación para medir la actividad de un paciente |
JP5041584B2 (ja) * | 2007-02-01 | 2012-10-03 | 独立行政法人国立循環器病研究センター | 神経信号用プローバ、神経信号出力装置、神経信号記録装置、神経刺激装置及び神経信号入出力装置 |
US20090082637A1 (en) * | 2007-09-21 | 2009-03-26 | Michael Galperin | Multi-modality fusion classifier with integrated non-imaging factors |
US20090082692A1 (en) * | 2007-09-25 | 2009-03-26 | Hale Kelly S | System And Method For The Real-Time Evaluation Of Time-Locked Physiological Measures |
GB0820629D0 (en) * | 2008-11-11 | 2008-12-17 | Univ Bath | Biocompatible electrode |
US10413203B2 (en) * | 2012-03-27 | 2019-09-17 | Cardiac Pacemakers, Inc. | Baseline determination for phrenic nerve stimulation detection |
-
2014
- 2014-01-30 GB GBGB1401613.3A patent/GB201401613D0/en not_active Ceased
-
2015
- 2015-01-29 EP EP15702568.5A patent/EP3100138B1/en active Active
- 2015-01-29 WO PCT/GB2015/050217 patent/WO2015114347A1/en active Application Filing
- 2015-01-29 JP JP2016549484A patent/JP6617108B2/ja active Active
- 2015-01-29 CN CN201580006550.7A patent/CN106062669B/zh not_active Expired - Fee Related
- 2015-01-29 US US15/115,616 patent/US10820816B2/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102096468A (zh) * | 2011-01-20 | 2011-06-15 | 中山大学 | 一种基于脑-机接口的家电遥控装置及方法 |
Non-Patent Citations (1)
Title |
---|
M. RIZK等: "a fully implantable 96-channel neural data acquisition system", 《JOURNAL OF NEURAL ENGINEERING》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109376651A (zh) * | 2018-10-24 | 2019-02-22 | 广东工业大学 | 一种基于cuda架构的gpu加速锋电位分类的系统 |
CN109460785A (zh) * | 2018-10-24 | 2019-03-12 | 广东工业大学 | 一种基于gpu加速的锋电位实时分类的方法 |
CN109376651B (zh) * | 2018-10-24 | 2021-09-03 | 广东工业大学 | 一种基于cuda架构的gpu加速锋电位分类系统 |
CN109460785B (zh) * | 2018-10-24 | 2021-09-24 | 广东工业大学 | 一种基于gpu加速的锋电位实时分类的方法 |
CN114925734A (zh) * | 2022-07-20 | 2022-08-19 | 浙江大学 | 一种基于神经拟态计算的在线神经元分类方法 |
CN114925734B (zh) * | 2022-07-20 | 2022-11-25 | 浙江大学 | 一种基于神经拟态计算的在线神经元分类方法 |
Also Published As
Publication number | Publication date |
---|---|
EP3100138A1 (en) | 2016-12-07 |
JP6617108B2 (ja) | 2019-12-04 |
EP3100138B1 (en) | 2019-09-11 |
CN106062669B (zh) | 2020-03-20 |
GB201401613D0 (en) | 2014-03-19 |
US10820816B2 (en) | 2020-11-03 |
JP2017509965A (ja) | 2017-04-06 |
US20170164852A1 (en) | 2017-06-15 |
WO2015114347A1 (en) | 2015-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106062669B (zh) | 用于脑机接口的系统 | |
Naderi et al. | Analysis and classification of EEG signals using spectral analysis and recurrent neural networks | |
Yang et al. | A hardware-efficient scalable spike sorting neural signal processor module for implantable high-channel-count brain machine interfaces | |
Mullen et al. | Real-time neuroimaging and cognitive monitoring using wearable dry EEG | |
Bertrand | Distributed signal processing for wireless EEG sensor networks | |
US11992685B2 (en) | Neurostimulator and method for delivering a stimulation in response to a predicted or detected neurophysiological condition | |
US9854985B2 (en) | Brain signal telemetry and seizure prediction | |
US20100081958A1 (en) | Pulse-based feature extraction for neural recordings | |
WO2020187109A1 (zh) | 一种用户睡眠检测方法及系统 | |
Barsakcioglu et al. | An analogue front-end model for developing neural spike sorting systems | |
US20080243022A1 (en) | Seizure prediction using brain signal telemetry | |
Bhattacharya et al. | Time domain multi-feature extraction and classification of human hand movements using surface EMG | |
US9477640B2 (en) | Neural signal processing and/or interface methods, architectures, apparatuses, and devices | |
Delgado-Restituto et al. | System-level design of a 64-channel low power neural spike recording sensor | |
Güngör et al. | A 1.2 nW analog electrocardiogram processor achieving a 99.63% QRS complex detection sensitivity | |
Iniewski | CMOS Biomicrosystems: Where Electronics Meet Biology | |
Uvaydov et al. | Aieeg: Personalized seizure prediction through partially-reconfigurable deep neural networks | |
TWI610267B (zh) | 基於個人化基底的壓縮感知系統及其方法 | |
Ossaba et al. | Open source multichannel EMG armband design | |
Chen et al. | Wearable algorithms: An overview of a truly multi-disciplinary problem | |
Heelan et al. | A mobile embedded platform for high performance neural signal computation and communication | |
Zjajo | Power-efficiency of signal processing circuits in implantable multichannel brain-machine interface | |
Victoria et al. | Integrated circuits and devices application in wearable healthcare systems | |
Wang | Advancing Myoelectric Prothesis Control via EMG-Force with Transfer Learning and via Wireless Sensor Time Synchronization | |
Zhemelev | Electromyography Signals Processing for Gait Phase Recognition |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200320 Termination date: 20220129 |
|
CF01 | Termination of patent right due to non-payment of annual fee |