专利名称：基于皮质脑电mu节律小波分析的术中运动区功能定位系统的制作方法大脑功能区病变，主要指位于运动、感觉和语言区的肿瘤，血管畸形和癫痫灶，其发病率由世界卫生组织在我国组织的大規模的调查报告仅仅癫痫的患病率就有8%。，我国现有癫痫病人1000多万人，其中药物难治性癫痫占癫痫病人的30%左右，我国目前有300万的难治性癫痫病人需要手术治疗，这还不包括位于功能区的低级别胶质瘤，转移瘤，原发 良性肿瘤，海绵状血管瘤和动静脉畸形等。大脑功能区病变不仅严重威胁人的生命，而且严重影响病人的生存和生活质量，同时造成的个人、社会和经济负担都是长久且巨大的，已成为严重的社会、经济和人文关怀问题。神经外科手术治疗是大脑功能区病变首选治疗方法之一，通过功能区定位确定大脑神经脑功能区边界，帮助医生最大限度地切除病灶而控制肿瘤的生长和复发，同时尽可能地保护病灶周围的正常脑组织，避免神经功能损害，保留正常的神经功能，关系到病人术后的生存质量。如何术中准确实时“脑功能区”定位就是此类手术的关键。目前，神经皮质(运动区)功能定位的方法主要包括显微神经外科技术、神经影像技术、神经电生理技术等方法。经典解剖功能定位对于临床医学具有重要意义，但有一定误差，由于个体差异及肿瘤的占位效应，引起功能区推移和重塑，经典解剖功能定位误差可达20_。依靠影像技术的高分辨率螺旋CT及功能型磁共振(f-MRI)，以及单光子发射计算机断层扫描(SPECT)、正电子发射计算机断层扫描(PET)、脑磁图(MEG)及手术导航系统多可以做到皮质生理解剖定位，但影像学方法存在一定假阳性，尚不能实时监测手术进程以及确定脑功能的状态。功能型磁共振(f-MRI)是依靠脑血流中血氧水平进行功能定位，病变影响脑皮层的血液供应会出现最大可达20mm的误差。正电子发射计算机断层扫描(PET)系统也可以对脑代谢活跃的区域进行定位，但是它与电生理刺激所显示的功能区，仅有65%的符合率。基于电生理技术的术中皮质或皮质下直接电刺激术可实时确定运动、感觉、语言甚至记忆等脑功能的皮质和皮质下功能区定位，是目前最准确、可信的常用脑功能区定位方法，基于电生理技术的术中皮质或皮质下直接电刺激术的精确度可达5 _左右；但是存在电刺激可能损伤大脑皮质、触发癫痫和二次手术等问题，而且操作时间长达0. 5至数小时。上述功能区定位方法的缺陷已表现在神经外科手术治疗实践中，传统手术的功能定位技术不能完全分辨和掌握功能结构与病变的关系，极易在切除病灶时导致大脑功能结构损害，据统计传统手术的永久性神经功能损害并发症为13-27%。另外，由于功能区病变手术容易出现严重并发症，也使得手术医生手术切除不积极，常常进行姑息性切除，如低级别胶质瘤的完全切除和次全切除率仅为43%。这样不仅使病变术后治疗变得困难，而且容易造成疾病的复发或症状难以控制，严重影响治疗预后。由此可见，目前的神经皮质(运动区)功能定位方法在速度、准确和安全性方面不能完全满足脑功能区手术需要。如何能在术中准确、快速、无创，甚至非唤醒状态下定位脑功能区是一直困扰临床和神经医学研究的基础理论问题，亟待解決。目前国内外尚未见有ー种基于皮质脑电mu节律小波分析的术中运动区功能定位系统的报道；同时，国内外也尚无在临床上应用的基于皮质脑电mu节律小波分析的术中运动区功能定位系统。因此研发具有自主知识产权的基于皮质脑电mu节律小波分析的术中运动区功能定位系统，实现准确、快速、无创的脑运动区功能定位，将帮助医生最大限度地切除病灶，同时尽可能地保护正常脑功能，提高患者术后生存质量，对未来大脑外科手术具有重大的应用价值。同吋，为下一歩皮质脑电定位高级认知功能皮质的生物机理科学研究提供新的技术方法手段，对未来大脑高级认知功能科学研究具有重大意义。具有巨大的社会和经济效应前景。发明内容本实用新型的目的在于针对现有技术的缺陷，以运动区特异性脑电mu节律为原理，结合小波变换，公开ー种基于皮质脑电mu节律小波分析的术中运动区功能定位系统。该系统能够准确、快速、无创地检测运动功能区脑电信号输入，并完成脑运动区功能定位图的输出，通过脑运动区皮质脑电信号mu节律的特异性分析，实现人体神经系统大脑皮质运动区功能定位的准确、快速、无创临床应用。mu节律是大脑的感觉运动区皮质的特异性脑电节律,肢体的真实运动或想象运动会在感觉运动皮层区域引起mu和beta节律的事件相关去同步化(ERD)和事件相关同步化(ERS)，而且不同肢体运动的ERD/ERS在皮质运动区上的空间分布也符合躯体特定区域分布的特征。因此，通过检测皮质运动功能区中存在着的特异性脑电mu节律，及其在肢体运动时产生的ERD/ERS在皮质上的空间分布，可以静态和动态检测定位皮质运动功能区的空间分布。小波变换具有多分辨特性，利用Mallat的分解与重构快速算法可以从运动区脑电中提取出mu节律，为运动区特异性脑电信号mu节律检测提供了有力工具。基于上述原理，本实用新型所采用的技术方案如下所述一种基于皮质脑电mu节律小波分析的术中运动区功能定位系统，包括脑电信号采集模块，信号处理模块，功能定位地图输出模块，所述信号处理模块包括脑电信号预处理単元、mu节律特征提取単元和模式分类单元；脑电信号采集模块采集的脑电信号，经由脑电信号预处理単元进行预处理滤波，传送至mu节律特征提取单元提取特异性mu节律特征，再通过模式分类单元进行分类，最后通过功能区定位地图输出模块反馈定位結果。所述脑电信号采集模块包括植入式电极、放大滤波器和A/D转换器，植入式电极采集脑电信号，经由放大滤波器进行放大滤波处理，然后通过A/D转换器将脑电信号转换为数字信号，最后输入到信号处理模块。所 述植入式电极为硬膜钼电极，包括钼6*8或8*8电极阵列，电极直径为4mm，相邻电极间距为10mm。植入式电极安放在人的大脑皮质上。放大滤波器和A/D转换器采用Synamps2放大器，用于电极检测信号的放大和数字化。所述脑电信号预处理単元的预处理滤波包括多尺度分解。所述多尺度分解利用离散db3小波变换进行7层小波分解，采用的小波Mallat算法的分解和重构算法见公式(3)。所述mu节律特征提取单元提取d6单层细节系数，然后进行全点数重构，其重构后的信号Sd6作为mu节律输出。所述模式分类单元以40%为特征阈值对mu节律进行是/否分类，识别特异性电极。 Ximi m=kH /_ M]もkH 丨丨/_JCik)A}啊 il/WI .KJT
(3)；其中，仏G为时域中的小波分解滤波器，A、^ 为时域中的小波重构滤波器；t为离
散时间序列，(=1，2，......，#;プ为分解层数,プ=1,2，......，ゾ，ゾ为分解深度，/(0为原始信
号。も为/⑴在第プ层近似部分的小波系数;沁为/a)在第プ层细节部分的小波系数。重构各单子频带信号时，只提取单层近似或细节系数，其余系数置0，然后对该单层系数进行全点数重构。重构信号特征量(运动事件发生前后ERD时间内能量比ERD)的计算见公式(4)。ERD = (EA-ER)/ER ⑷其中，ER为运动事件前ERD时间窗口内的各子频带重构信号的每个采样点值的平方和，EA为计算运动事件后ERD时间窗口内的各子频带重构信号的每个采样点值的平方和。所述功能区定位地图输出模块输出的运动特异性功能区定位图，是以模式分类单元识别的特异性电极坐标为边界点拟合边界曲线，即运动特异性功能区定位图，闭合曲线围成的区域为运动特异性功能区。本实用新型相对现有技术具有如下优点和效果(I)特异性检出正确率高本实用新型基于运动功能区mu节律的事件ERD特异性，选择了合理的特征频帯、特征值，以及特征提取和分类算法具有可靠的检测原理，从根本上保证了特异性检出正确率，其特异性检出正确率达到78%-100%。(2)电极检测精度高系统采用的电极具有4mm直径和IOmm的植入式电极间距离，具有较高的空间和频率分辨率，能提供电极中心点附近5_半径范围内神经元的电活动信息。特征阈值的实质是在电极中心点附近5mm半径范围内检测到的最小有效特征量。因此计算本系统的定位空间微观精确度可达5 mm。与术中皮质或皮质下直接电刺激术相比，本系统进一步提闻了检测精确度。[0027](3)检测速度快本系统以自发脑电mu节律为运动功能区的特异性脑电为检测对象，自发脑电mu节律的ERD的采样实验周期和时间窗ロ为4秒，因此，理论上采样速度为4秒。考虑可靠性采用10次采样实验结果共同判断的方法，再加上计算机处理的时间，本系统一次功能定位检测时间为60秒。与术中皮质或皮质下直接电刺激术操作时间长达0. 5至数小时相比，本系统极大提高了检测速度，极大減少了医生手术时间和病人的痛苦，节省巨大的人力和物力，具有良好的经济和人文关怀价值。[0028](4)检测无创性本系统提取自发脑电mu节律采用脑电信号的被动检测方式，无主动刺激造成的创伤。避免了术中皮质或皮质下直接电刺激术可能损伤大脑皮质、触发癫痫等问题。极大減少了医生手术时间和病人的痛苦，节省巨大的人力和物力，具有良好的经济和人文关怀价值。

图I为运动区脑功能定位系统结构图。图2为脑电信号采集系统模块结构图。图3为原始信号的小波分解与重构。图4为原始信号的小波变换滤波。图5为运动特异性功能区定位图。

以下结合附图和实施例对本实用新型作进ー步说明，但本实用新型的实施不限于此。一种基于皮质脑电mu节律小波分析的术中运动区功能定位系统，如图I所示，包括脑电信号采集模块，信号处理模块，功能区定位地图输出模块，信号处理模块包括脑电信号预处理単元、mu节律特征提取単元和模式分类单元。脑电信号采集模块的组成如图2所示，包括植入式电极、放大滤波器和A/D转换器。该系统通过植入式电极的硬膜下电极阵列采集皮质脑电信号ECoG，经由放大滤波器进行放大滤波处理，然后通过A/D转换器将脑电信号转换为数字信号，再输入到信号处理模块；信号处理模块的脑电信号预处理単元通过小波分析的分解与重构算法，对各电极采集的数据进行预处理滤波，传送至mu节律特征提取单元提取特异性mu节律特征，再通过模式分类单元进行分类，最后通过功能区定位地图输出模块完成功能区定位图处理输出。在每位病人的大脑皮质上安放硬膜下电极阵列提取ECoG数据，为病人植入硬膜钼 6*8 或 8*8 电极阵列(subdural electrode arrays (Ad-Tech, Racine, WI),姆个电极直径为4mm,相邻电极间的距离为10_ ；Synamps2放大器(Neuroscan, ElPaso, TX)用于电极检测信号的放大和数字化，ECoG数据采样率为1000Hz，通过0. 05-200Hz通带滤波。信号处理模块由综合计算机处理系统提供，向病人提供运动指示并储存指示时间，接收并储存Synamps2放大器的脑电信号数据。毎次实验时，根据计算机显示屏的“运动一休息”指示，病人先运动手指2秒，然后休息2秒；再重复多次上述相同实验。采集位于大脑神经皮质运动功能区的特定皮质区域ECoG数据，共进行10次采集数据，用于计算机系统的处理工作。[0039]根据mu节律和エ频干扰的频带确定分解层次本实用新型所研究的脑电信号包含mu节律(8 - 12 Hz)、一些瞬变信号以及エ频干扰(50 Hz)，经表I频带频率计算公式确定分解层次为7。表I :各单子频带重构信号的ERD指标特征值