基于原发性癫痫的静息态功能磁共振的数据处理方法
【技术领域】
[0001]本发明涉及针对疾病检测的数据处理领域，特别涉及一种基于原发性癫痫的静息态功能磁共振的数据处理方法。
【背景技术】
[0002]癫痫(epilepsy)是常见的由于多种病因引起的慢性神经系统疾病，以脑部神经元放电的兴奋性和同步性异常所致过度放电引起的突然、反复和短暂的中枢神经系统功能失常为特征。截止2014年，全球大约有5500万癫痫病人，其发病率高于多发性硬化、帕金森病等常见的神经系统疾病。目前认为，癫痫已经成为神经系统疾病中仅次于脑卒中的第二大常见病，癫痫性猝死也日益引起研究者的关注。在现实生活中，对癫痫患者的歧视和不理解极大地影响着患者及其家人。总之，癫痫是一种致残率高、病程长和以临床反复发作为特点，同时还严重威胁患者身心健康的疾病。
[0003]复杂部分性发作是原发性癫痫(theid1pathic epilepsy) —种常见的类型，指脑部无明显器质性或代谢性异常表现，发病机制尚未完全阐明、可能与遗传因素相关并以不同程度的意识障碍、自主神经症状和精神症状发作为特征的一组癫痫综合征，通常属于药物难治性，约占各类癫痫总数的20 %?30 %。对于CPS患者，很大一部分(约60 %?90 % )可以通过手术切除得以根治，因此术前对致痫灶的准确定位就显得尤为重要。但CPS患者常规的颅脑CT和MRI检查一般没有可见的结构变化。因此探讨致痫灶的无创、精确定位方法一直是癫痫研究的热点。
[0004]头皮脑电图(Scalp Electroencephalogram)检查虽然可以证实癫痫的异常放电，但其空间分辨率较低，与皮层间隔有颅骨、头皮、外界干扰等不利因素的影响，对致痫灶定位的敏感性和特异性低;皮层脑电图或深部脑电图能直接探测脑异常放电区，但系有创性检查且电极所探及的区域相当有限，实际应用受到一定的限制。
[0005]目前，对癫痫的神经影像学研究主要有以下几种:常用结构MRI从形态学方面对海马体积的萎缩程度进行测量及对海马进行波谱分析以确定病变的海马，或者用正电子发射型计算机断层显像和单光子发射计算机断层成像技术从代谢和血流的角度反映癫痫活动相关脑区的功能异常，再就是利用各种刺激诱导局部脑组织血氧水平依赖信号发生变化，间接反映神经元的活动，即“事件/任务相关功能磁共振成像(event/task-relatedfMRI)，，。
[0006]海马体积测量及波谱分析存在主观性;PET和SPECT是侵入性检查，有放射性且空间分辨率较低;事件相关fMRI不但需要设计好合适的任务，而且要求病人高度配合。
[0007]原发性癫痫的静息态功能磁共振:resting-state funct1nal magneticresonance imaging，rfMRI；
[0008]复杂部分性发作:Thecomplex partial seizures，CPS;
[0009]皮层脑电图:Electrocorticogram，ECoG;
[0010]深部脑电图:Deepelectroencephalogram,DEEG；[0011 ]正电子发射型计算机断层显像:Positron Emiss1n Computed Tomography，PET;
[0012]单光子发射计算机断层成像:Singlephoton emiss1n computed tomography,SPECT；
[0013]脑组织血氧水平依赖(bloodoxygenat1n level dependent,BOLD；
[0014]局部一致性:Thereg1nal homogeneity，ReHo;
[0015]低频振幅:Thelow frequency fluctuat1n,ALFF；
[0016]种子相关功能连接:Theseed based funct1nal connectivity ,FC0
【发明内容】
[0017]本发明针对现有技术的不足，提供一种基于原发性癫痫的静息态功能磁共振的数据处理方法，用以解决有创性检查及定位不精确等问题。
[0018]为解决以上问题，本发明采用的技术方案如下:一种基于原发性癫痫的静息态功能磁共振的数据处理方法，包括如下步骤:
[0019]S1.数据采集:常规检查序列包括轴位T2加权成像、Tl加权成像和弥散加权成像，以平行前后联合的水平进行全脑轴位扫描，以排除颅内器质性病变;T1-3D BRAVO序列扫描，作为结构像将数据处理后的功能像叠加其上，用以显示正激活和负激活的脑区;静息态扫描，采用梯度回波单次激发回波平面成像；
[0020]S2.数据预处理:1.对每个被试对象的fMRI原始数据转换为.1mg的格式;2.校正原始图像中层与层获得的时间的不同；3.采用迭代算法对6个参数进行头部刚体模型配准，排除头部沿着X，Y，Z三个方向的平移和围绕这三个方向的旋转超过2mm，或者任意角度旋转大于2°的被试对象;4.将每个被试的图像放入统一的标准空间进行标准化，将其转为朝向和大小都相同的图像;5.空间平滑:使用全宽半高高斯核对fMRI图像进行空间平滑处理;6.对空间平滑后的fMRI图像进行去线性漂移;7.去除高频呼吸心跳及低频噪声；
[0021 ] S3.数据分析:包括ReHo数据分析、ALFF数据分析和FC数据分析，所有的ReHo、ALFF、FC图均单独分析，然后汇总确定有改变的脑区。
[0022]作为优选，ReHo数据分析:计算出脑内每个体素与其周围相邻的3个或3个以上体素在时间序列上的一致性，得出该体素的肯德尔和谐系数，即为该体素的ReHo值，每个体素均有一个ReHo值，所有体素构成了受试者的ReHo脑。
[0023]作为优选，ALFF数据分析:首先保留低频段0.01-0.08Hz内所有频率点上的BOLD信号，将这些数据进行傅立叶变换计算功率谱，功率谱的峰下面积可视为信号的能量，然后对其开方，可得到该信号的振荡幅度，也就是BOLD信号的波动变化强度，进而绘制振幅的统计参数脑图，来描述体素自发活动的强弱。
[0024]本发明的有益效果:本发明的数据处理方法对多灶性癫痫病灶的探索极具价值，避免了复杂的任务刺激导致的脑活动差异，使不同受试者的实验条件基本控制在相同的水平，且易于配合，更容易得到可靠的结论。它从脑功能的角度了解癫痫患者某些临床症状的发病机理，能帮助定位癫痫灶和周围功能区，指导癫痫手术及癫痫灶的切除范围；有助于理解癫痫的病理生理改变和认识脑电图异常波形的意义。随着临床研究的不断扩展rfMRI检查技术、计算方法的改进，可进一步揭示不同脑区更精细的分工及与高级整合之间的相互关系，其对药物难治性癫痫的诊治将具有更大的应用价值。
【附图说明】
[0025]图1为头动校正、空间平滑、低频滤波等预处理后得到的原始图像示意图；
[0026]图2为使用ReHo、ALFF、FC等计算后得到的功能磁共振图像示意图；
[0027]图3为使用Alphasim统计学处理后与正常人脑解剖图融合得到的最终图像示意图。
【具体实施方式】
[0028]为使本发明的目的、技术方案及优点更加清楚明白，以下参照附图并举实施例，对本发明做进一步详细说明。
[0029]实施例:
[0030]1.数据采集工具及对象:
[0031]本实施例使用美国GE公司的3.0T Signa HDx超导型MR扫描仪(GE-Signa HDx，Mi lwaukee，US.)，8NV—Head线圈，以头部线圈作为发射和接受线圈。扫描时受试者清醒、闭目艮、头部固定并最大限度地减少头部及其它部位的主动与被动运动、安静平卧于检查床，均匀呼吸，不执行特定的认知任务，同时要求尽量不行主动思维活动，待受试者熟悉环境后开始扫描。病例组在扫描前进行2-5分钟过度换气实验以增加间期痫样放电的数量;为了保证试验的准确性，正常志愿者也进行2-5分钟的过度换气实验。
[0032]2.数据采集方法:
[0033]相关扫描序列及参数如下:
[0034]2.1常规检查序列包括轴位T2加权成像(T2WI)、T1加权成像(TlWI)和弥散加权成像(DWI)，以平行前后联合的水平进行全脑轴位扫描，以排除颅内器质性病变。
[0035](I) T2WI 扫描参数:TR/TE = 6280/104ms，层厚:5mm，间隔:I.5mm，NEX:1，FOV:24x24cm，Matrix:320x320，扫描层数为19层，扫描时间为82秒。
[0036](2)T1WI扫描参数:TR/TE = 2825/24ms，层厚:5mm，间隔:I.5mm，NEX:1，FOV:24xl8cm，Matrix:320x320，扫描层数为19层，扫描时间为99秒。
[0037](3 )DWI扫描参数:TR/TE = 4450/72.9ms，层厚:5mm，间隔:1.5mm，NEX: 2，FOV:24x24cm,Matrix: 128xl28，b = 1000s/mm，扫描层数为19 层，扫描时间为 40 秒。
[0038]2.2T1-3D BRAVO序列扫描，作为结构像将数据处理后的功能像叠加其上，用以显示正激活和负激活的脑区O参数:TR/TE = 7.8/3.0ms，层厚:Imm，间隔:O，NEX: 0.5，FOV:24x24cm,Matrix:256x256,FA(F1 ip Angle):13°。从右向左采集，扫描范围包括全脑，采集层数为176层，扫描时间为209秒。
[0039]2.3静息态扫描，采用梯度回波(GRE)单次激发回波平面成像(EPI)技术，参数:TR/TE = 2000/30ms，层厚:4mm,间隔:Omm，FOV = 24x24cm，Matrix: 64x64，FA(Fl ip Angle):90°，从下向上采集，扫描范围包括全脑，层数:38，采集256个时间点，扫描时间为