本发明属医学图像处理及应用领域，涉及意识障碍者意识恢复预测的功能影像分析方法，具体涉及一种基于意识障碍者功能的磁共振检测的数据分析方法。
背景技术：
现有技术公开了昏迷(coma)是在各种致病因素作用下,脑干网状上行激动系统或大脑皮质内神经元广泛受损引起的严重意识障碍,是多种疾病发展过程中的一个阶段。随着交通事故、心脑血管疾病等意外事件的逐年增多及重症医学的快速发展，存活的严重脑损伤患者不断增多，致使意识障碍患者的数量不断增加。
常见的意识障碍类型主要包括昏迷(coma)，植物状态(vegetativestate,vs)和微意识状态(minimallyconsciousstate,mcs)，临床实践显示，维持此类患者生存的支持治疗及日常护理的费用非常昂贵，给社会造成极大的负担，给患者家庭带来沉重的精神打击和巨大的经济压力，并伴随引发一系列伦理、法律等问题。业内认为，意识障碍患者的预后预测是目前医学领域亟待解决的重大难题，准确可靠的预后判断将为患者家属及临床医师制定合理有效的医疗决策提供有益的信息。
以临床量表为主的行为学判断是目前最为常用的评估方法，wijdicks等进行了一项汇总分析，格拉斯哥昏迷评分(glasgowcomascale,gcs)可很好地预测混合病因昏迷患者的预后，其受试者工作曲线(roc)下面积达0.87；estraneo等发现反应性恢复的患者具有更高的昏迷恢复量表(comarecoveryscale-revised，crs-r)得分和低残障等级量表(disabilityratingscale，drs)得分，正中神经sep存在和crs-r>6分是反应性恢复的有效预测因子，但，该类方法主观性极强，其结果易受评估人经验及接受培训程度的影响。
神经电生理具有极好的时间分辨率，并可连续床旁监测，gütling等发现脑电反应性(eegreactivity，eeg-r)区分预后良好和预后不良的正确率达92％，而中枢传导时间(centralconductiontime，cct)的正确率仅82％，gcs的正确率为72％，但神经电生理其空间分辨率差随着病情进展，该类指标的敏感性逐渐下降。
近年神经影像学在意识障碍脑功能判断中日益发挥重要作用，受到业内极大的关注，tollard等将弥散张量成像(diffusiontensorimaging,dti)和磁共振波谱分析(mrs)结合用于预测脑外伤患者的长期预后，其中观察到将dti和mrs结合后对1年后未恢复预测的敏感性达86％；lancet的一项研究采用pet和fmri，pet对102例患者75例的预测结局正确(74％)，而fmri对65例患者中36例预测正确(56％)。
基于现有技术中对于意识障碍病人中具有被唤醒可能性预测的研究目前仍在摸索阶段的现状，本申请的发明人拟提供一种意识障碍者意识恢复预测的功能影像分析方法，拟通过对患者的神经功能影像数据的分析，获得能用于临床参考的初步预判方法，进一步更好地为临床及社会医疗实践提供有效的指导，有助于更好的配置医疗资源。
与本发明相关的现有技术有：
1.wijdicksef,rabinsteinaa,bamletwr,mandrekarjn:fourscoreandglasgowcomascaleinpredictingoutcomeofcomatosepatients:apooledanalysis.neurology2011,77(1):84-85.
2.estraneoa,morettap,termet,trojanol:predictorsofrecoveryofresponsivenessinprolongedanoxicvegetativestate.authorreply.neurology2013,81(14):1274-1275.
3.gutlinge,gonsera,imhofhg,landist:eegreactivityintheprognosisofsevereheadinjury.neurology1995,45(5):915-918.
4.tollarde,galanaudd,perlbargv,sanchez-penap,lefury,abdennourl,cozzonep,lehericys,chirasj,puybassetl:experienceofdiffusiontensorimagingand1hspectroscopyforoutcomepredictioninseveretraumaticbraininjury:preliminaryresults.criticalcaremedicine2009,37(4):1448-1455.
5.stenderj,gosserieso,brunoma,charland-vervillev,vanhaudenhuysea,demertzia,chatellec,thonnardm,thibauta,heineletal:diagnosticprecisionofpetimagingandfunctionalmriindisordersofconsciousness:aclinicalvalidationstudy.lancet2014,384(9942):514-522.。
技术实现要素：
本发明要解决的技术问题是提供一种意识障碍者意识恢复预测的功能影像分析方法，具体涉及一种基于意识障碍者功能的磁共振检测的数据分析方法。
本发明基于人脑功能网络构建使用的传统方法是皮尔逊相关系数(如图2传统方法1，图4a和图5a所示)，但是该方法没有考虑到多个脑区之间的关系，仅考虑了任意两个脑区之间的功能连接，因此无法揭示复杂的大脑功能连接网络以及利用稀疏表征方法(如图2传统方法2，图4b和图5b所示)则采用机器学习的思想，考虑到多个脑区之间的关系，因此能更好的刻画这种复杂网络的特点，同时，传统稀疏表征方法考虑了功能连接的稀疏性，符合人脑网络经济高效的假设，可是，该方法对所有脑区等同对待，忽视了强连接和弱连接的区别，其优化结果往往错误地压低了强连接，或者错误地提升了弱连接，使得脑功能连接网络矩阵失去了结构性和模块化特征，而这一特征是人脑网络的重要特性等基础现状，
本发明提供了一种准确率高、重复性好的预测脑损伤长期无意识患者能否恢复意识的新的算法-机器学习的数据分析方法，其包括，基于静息态功能磁共振数据(rs-fmri)，采用“加权组稀疏”算法构建人脑功能连接矩阵，从该矩阵中，用稀疏表征的特征筛选方法挑选对分类贡献高的功能连接特征用于自动预测；采用并采用线性支撑向量机构建预测模型，得到最终的对意识恢复与否的预测结果。
更具体的，本发明的一种基于意识障碍者功能的磁共振检测的数据分析方法，其包括步骤：
1)对rs-fmri数据进行预处理，
所述rs-fmri数据包括层获取时间校正，头动校正，对齐到标准空间，空间平滑，时域带通滤波，从数据中去除白质、脑脊液平均信号和头动曲线；
本方法中，采用脑区分割图谱(如图1a所示)提取不同脑区(或称感兴趣区，regionofinterest/roi)的平均时间序列；该时间序列反映区域平均血氧水平依赖信号(如图1b上所示)；
2)对步骤1)获得的数据，采用“加权稀疏表征”算法，优化地计算任意一个脑区和其他脑区信号之间的表征关系，采用脑区信号之间的相关系数加权的l-1范数作为约束项(如图1b下所示)，获得在同时考虑其他脑区的影响的情况下，任意两个脑区之间的关系，即一个方阵，该矩阵代表人脑功能连接网络(如图1c所示)；
所述的步骤2)包括如下分步骤：
2)-1，对脑区i，其时间序列由其他脑区的线性组合表达，此为图1b下的模型的第一项，通过计算其线性组合系数w的l-1范数(图1b下的模型的第二项)进行约束；本发明的实施例中，特别地，在该约束项上进行加权，加权系数通过脑区i和另外任一脑区的皮尔逊相关系数的e负指数变形得到；通过优化该目标函数，即最小化该目标函数取值，得到最佳w，即为脑区i和所有其他脑区之间的功能连接系数；
2)-2对各个脑区分别进行步骤2)-1，获得每个脑区和所有其他脑区之间功能连接系数，由于步骤2.1得到的为一个列向量，所有脑区所对应的列向量叠在一起形成一个方阵；
2)-3对上述方阵进行转置，然后与转置前的方阵相加并除以2，实现方阵的对称化；因为传统意义上的脑功能连接矩阵为对称阵，即脑区i和脑区j的功能连接等于脑区j和脑区i的功能连接；
3)对所有n个被试采用步骤1)-2)计算人脑功能连接网络，
从所有n个意识丧失被试中，选择n-1个被试作为训练数据(如图2所示)，将这些被试的意识恢复结果作为标号，采用稀疏表征的方法对所有被试所有功能连接系数(人脑功能连接网络方阵里的上三角矩阵)进行特征筛选，筛选的原则是选择较少的对分类目标贡献较大的特征，
通过下述方法选择较少的对分类目标贡献较大的特征：
建立一个以特征权重的l-1范数为约束项的线型表征优化模型，其中每个功能连接强度乘以一个特征权重，用所有功能连接强度的加权和来逼近训练数据的标号，l-1范数可以得到稀疏的优化结果，即只有少数特征被选择，该特征权重用于特征筛选，非零权重所对应的特征即被选中；
4)用被选中的特征，在特征所构成的高维空间中构建线型支撑向量机(supportvectormachine/svm)，学习最优分类面，该分类面由所有进入svm中的特征的权重系数决定；
5)将剩下的一个意识丧失被试，作为测试数据，对该被试的人脑功能连接网络矩阵，用步骤3)所得到的有用特征序号，提取相同特征(相同位置的功能连接强度)，将其放入步骤4)中svm所学习到的分类模型中(应用步骤4学习到的权重系数)得到预测结果；图2表示了步骤3-5所示的过程；
6)更换测试被试，将剩下的被试作为训练数据，重复步骤3-5.每次均能得到对测试被试的预测结果，将所有预测结果和该测试被试的标号相比较，得到预测模型的准确率、敏感度、特异度，制作受试者操作曲线(receiveroperativecurve/roc)并得到线下面积；
7)评估模型，经过利用大数据集的实验，以及步骤1-6，构建出一个最优分类模型，实验结果显示，包括了最优的脑区间功能连接特征的选择以及这些特征的加权方式预测苏醒与否的准确率超过89％，敏感度和特异度分别超过了88％和90％；
实验结果证实，本发明的方法的性能明显优于传统方法(如图3所示)；
8)新被试的预测，将所有n个被试作为训练样本，利用步骤3所述的特征选择确定对预测有用的特征，用步骤4重新构建分类模型；当有新被试时，可利用该模型预测新被试的苏醒结果，用于医生参考，具体的预测方法是，对该新被试的rs-fmri数据，采用步骤1所述方法进行预处理，并采用步骤2所述方法进行脑网络构建，最后，采用上述n个被试重新构建的大数据预测模型进行特征挑选和结果预测。
本发明基于现有技术的现状，提供了一种基于皮尔逊相关系数所得的“功能连接先验”对稀疏表征优化函数中稀疏性的约束项进行加权的一种基于意识障碍者功能的磁共振检测的数据分析方法，其中的加权方式使得强连接被尽可能保留，同时保持了稀疏表征的优点。本方法构建的人脑功能连接网络，同时具有模块化和稀疏性两个符合神经科学假设的重要特性，所提出的加权稀疏表征方法能够得到具有模块化和稀疏性特征的脑功能连接网络(如图4c所示)，结果优于传统方法(如图4a和4b所示)。
附图说明
图1、利用加权稀疏表征方法构建人脑功能连接网络。
图2、基于人脑功能连接网络的模式判别流程图。
图3、模型评估结果。
图4显示了实验结果，证明所提出的加权稀疏表征方法，其中(c)本发明方法，(a、b)传统方法。
图5、显示了传统方法(a、b)和本方法(c)在两组患者功能连接网络构建上的差异。
具体实施方式
实施例1意识障碍者的预后预测临床检测试验
1)对rs-fmri数据进行预处理，
所述rs-fmri数据包括层获取时间校正，头动校正，对齐到标准空间，空间平滑，时域带通滤波，从数据中去除白质、脑脊液平均信号和头动曲线；
本方法中，采用脑区分割图谱(如图1a所示)提取不同脑区(或称感兴趣区，regionofinterest/roi)的平均时间序列；该时间序列反映区域平均血氧水平依赖信号(如图1b上所示)；
2)对步骤1)获得的数据，采用“加权稀疏表征”算法，优化地计算任意一个脑区和其他脑区信号之间的表征关系，采用脑区信号之间的相关系数加权的l-1范数作为约束项(如图1b下所示)，获得在同时考虑其他脑区的影响的情况下，任意两个脑区之间的关系，即一个方阵，该矩阵代表人脑功能连接网络(如图1c所示)；
所述的步骤2)包括如下分步骤：
2)-1，对脑区i，其时间序列由其他脑区的线性组合表达，此为图1b下的模型的第一项，通过计算其线性组合系数w的l-1范数(图1b下的模型的第二项)进行约束；本发明的实施例中，特别地，在该约束项上进行加权，加权系数通过脑区i和另外任一脑区的皮尔逊相关系数的e负指数变形得到；通过优化该目标函数，即最小化该目标函数取值，得到最佳w，即为脑区i和所有其他脑区之间的功能连接系数；
2)-2对各个脑区分别进行步骤2)-1，获得每个脑区和所有其他脑区之间功能连接系数，由于步骤2.1得到的为一个列向量，所有脑区所对应的列向量叠在一起形成一个方阵；
2)-3对上述方阵进行转置，然后与转置前的方阵相加并除以2，实现方阵的对称化；因为传统意义上的脑功能连接矩阵为对称阵，即脑区i和脑区j的功能连接等于脑区j和脑区i的功能连接；
3)对所有n个被试采用步骤1)-2)计算人脑功能连接网络，
从所有n个意识丧失被试中，选择n-1个被试作为训练数据(如图2所示)，将这些被试的意识恢复结果作为标号，采用稀疏表征的方法对所有被试所有功能连接系数(人脑功能连接网络方阵里的上三角矩阵)进行特征筛选，筛选的原则是选择较少的对分类目标贡献较大的特征，
通过下述方法选择较少的对分类目标贡献较大的特征：
建立一个以特征权重的l-1范数为约束项的线型表征优化模型，其中每个功能连接强度乘以一个特征权重，用所有功能连接强度的加权和来逼近训练数据的标号，l-1范数可以得到稀疏的优化结果，即只有少数特征被选择，该特征权重用于特征筛选，非零权重所对应的特征即被选中；
4)用被选中的特征，在特征所构成的高维空间中构建线型支撑向量机(supportvectormachine/svm)，学习最优分类面，该分类面由所有进入svm中的特征的权重系数决定；
5)将剩下的一个意识丧失被试，作为测试数据，对该被试的人脑功能连接网络矩阵，用步骤3)所得到的有用特征序号，提取相同特征(相同位置的功能连接强度)，将其放入步骤4)中svm所学习到的分类模型中(应用步骤4学习到的权重系数)得到预测结果；图2表示了步骤3-5所示的过程；
6)更换测试被试，将剩下的被试作为训练数据，重复步骤3-5.每次均能得到对测试被试的预测结果，将所有预测结果和该测试被试的标号相比较，得到预测模型的准确率、敏感度、特异度，制作受试者操作曲线(receiveroperativecurve/roc)并得到线下面积；
7)评估模型，经过利用大数据集的实验，以及步骤1-6，构建一个最优分类模型，实验结果显示，包括了最优的脑区间功能连接特征的选择以及这些特征的加权方式预测苏醒与否的准确率超过89％，敏感度和特异度分别超过了88％和90％；
实验结果证实，本发明的方法的性能明显优于传统方法(如图3所示)；
8)新被试的预测，将所有n个被试作为训练样本，利用步骤3所述的特征选择确定对预测有用的特征，用步骤4重新构建分类模型；当有新被试时，可利用该模型预测新被试的苏醒结果，用于医生参考，具体的预测方法是，对该新被试的rs-fmri数据，采用步骤1所述方法进行预处理，并采用步骤2所述方法进行脑网络构建，最后，采用上述n个被试重新构建的大数据预测模型进行特征挑选和结果预测。