多模态脑部神经影像特征的处理方法技术

本发明专利技术多模态脑部神经影像特征的处理方法，涉及用于识别图形的图像特征或特性的抽取的图像预处理，先采用样本权重和低秩约束的多模态神经影像特征选择方法对多模态的数据进行特征选择，得到低维特征矩阵，计算每个模态的核矩阵，得到低维特征矩阵，计算每个模态的核矩阵，然后将不同模态的核矩阵融合成一个核矩阵，由此选择出更具有判别性的生物标志物特征，并使用多核支持向量机对阿尔茨海默病新样本病例进行预测分类，克服了现有技术中，所存在的利用的生物标志物特征会对受试者造成伤害，只利用一种脑部影像特征数据或利用不充分的多模态脑部神经影像特征数据不能找出患者的患病脑区，并且所利用的脑部影像中的特征没有医学上的解释性的缺陷。

Processing Method of Multimodal Brain Neuroimaging Features

The processing method of multi-modal brain neuroimaging features of the present invention involves image preprocessing for extracting image features or features for pattern recognition. Firstly, multi-modal neuroimaging feature selection method with sample weights and low rank constraints is used to select features of multi-modal data to obtain low-dimensional feature matrices, calculate the core matrices of each mode, and obtain low-dimensional feature matrices. The core matrix of each mode is calculated, and then the different modes of the core matrix are fused into a core matrix. The more discriminant biomarker features are selected, and the new samples of Alzheimer's disease are predicted and classified by using multi-core support vector machine, which overcomes the existing technology that the existing biomarker features can cause harm to the subjects. Only one biomarker feature is used to predict and classify the new samples of Alzheimer's disease. Brain imaging feature data or inadequate multi-modal neuroimaging feature data can not identify the diseased brain areas of patients, and the features used in brain imaging have no medical explanatory defects.

【技术实现步骤摘要】
多模态脑部神经影像特征的处理方法
本专利技术的技术方案涉及用于识别图形的图像特征或特性的抽取的图像预处理，具体地说是多模态脑部神经影像特征的处理方法。
技术介绍
阿尔茨海默病是一种很难治愈且不可逆转的脑疾病，现阶段治疗阿尔茨海默病的药物效果有限，所以在病变的早期及时对疾病进行干预非常重要。随着神经影像技术发展，现在已经通过不同技术手段，例如磁共振图像或正电子发射计算机断层扫描，来得到多种反映人类脑部状况的影像。医生借助学习得到的知识和实践经验，解读病人脑影像所反映的解剖结构和病理生理信息，但是这种人工解读病人脑影像的方式，是依赖医生的个人经验、知识和情绪，诊断的正确率和效率均较低。现今，利用计算机技术对脑影像进行分析和处理来辅助医生对脑部发生病变的区域进行定量分析，从而提高了诊断的正确率和效率。机器学习方法作为基于数据驱动的预测与分析工具，能够充分利用生物标志数据内在的结构信息构建模型来分析大脑结构或者功能的特性，从而更好地揭示脑认知行为和相关疾病的产生机制。“模态”是指从某一个角度对事物进行描述，那么多模态数据是从不同角度对同一事物进行描述所得到的数据，如果将一种脑部神经影像当作一个模态，那么从不同角度对脑部影像进行描述所得到的就是多模态脑部神经影像。多模态学习技术就是通过探索不同模态之间的关系，来提高多模态的选择和分类的性能。因此，可以采用多模态脑部影像来判断检查者是否患有阿尔茨海默病，或者处于阿尔茨海默病的早期轻度认知障碍，或者处于阿尔茨海默病的晚期轻度认知障碍。脑部影像的维度较高且包含的信息量较大，并不是脑部影像中所有特征都对检测分析阿尔茨海默病有帮助。因此从脑部影像所提供的大量特征中选择出与分类预测任务相关的特征，去掉冗余的或者与分类相关性低的特征是非常重要的研究工作。现有的从脑部影像所提供的大量特征中选择出与分类预测任务相关的特征的方法是基于最小绝对收缩和选择算子方法，该方法大致分为两类：一类是应用单模态数据或者简单级联多模态数据，此类方法不能充分利用多模态数据来挖掘不同模态之间的互补信息。另一类是通过脑脊液中的β-淀粉体(1-42)和t-tau蛋白生物标志物来检测阿尔茨海默病，此类方法需要侵入人体提取生物标志物，对受试者产生较大伤害。CN108198576A公开了基于语音特征非负矩阵分解的阿尔茨海默症初筛方法，该方法利用语音特征对受试者进行是否患病进行判断，其提取出的声学特征没有在临床上证明阿尔茨海默症和健康人存在差异，其提取的特征在医学上没有解释性。CN106062563A公开了用于阿尔茨海默病早期诊断的生物标志物及方法，该方法从脑脊液和血液中提取生物标志物，在对受试者进行这些生物标志物获取时会对受试者造成身体伤害，有侵入性，且特征提取和检验的成本较高。CN108154924A公开了基于支持向量机的阿尔茨海默症特征分类方法及系统，该方法只利用了一种数据，只能从结构上找出脑部变化，获取的特征少，不能充分说明受试者是否患病。CN108682009A公开了一种阿尔茨海默症预测方法、装置、设备及介质，该方法利用磁共振成像数据且使用3D卷积神经网络模型进行运算，该方法只能判断受试者是否患病，不能找出患者的患病脑区，其黑盒特性不能给出判断受试者是否患病的依据。CN106202916A公开了建立一种阿尔茨海默病的分层多流形分析模型，该方法只能对受试者进行分类，并不能找出受试者是否患病的病变脑区特征。CN108682009A公开了一种基于人工神经网络和多模态MRI的阿尔茨海默病智能判别分析方法，该方法只能判断受试者是否患病，不能找出患者的患病脑区，其黑盒特性不能给出判断受试者是否患病的依据。总之，在现有的检测分析阿尔茨海默病技术中，存在利用的生物标志物特征时会对受试者造成伤害，只利用一种脑部影像特征数据或利用不充分的多模态脑部神经影像特征数据不能找出患者的患病脑区，并且所利用的脑部影像中的特征没有医学上的解释性的缺陷。
技术实现思路
本专利技术所要解决的技术问题是：提供多模态脑部神经影像特征的处理方法，是基于样本权重和低秩约束的多模态神经影像特征选择和多核支持向量机分类的方法，即先采用样本权重和低秩约束的多模态神经影像特征选择方法对多模态的数据进行特征选择，得到低维特征矩阵，计算每个模态的核矩阵，得到低维特征矩阵，计算每个模态的核矩阵，然后将不同模态的核矩阵融合成一个核矩阵，由此选择出更具有判别性的生物标志物特征，并使用多核支持向量机对阿尔茨海默病新样本病例进行预测分类。本专利技术方法克服了在现有的检测分析阿尔茨海默病技术中，所存在的利用的生物标志物特征会对受试者造成伤害，只利用一种脑部影像特征数据或利用不充分的多模态脑部神经影像特征数据不能找出患者的患病脑区，并且所利用的脑部影像中的特征没有医学上的解释性的缺陷。上文中，“基于样本权重和低秩约束”的英文为“SampleWeightandLowRankConstraint”，以下简称“SWLRC”，“支持向量机”的英文为“SupportVectorMachine”，以下简称“SVM”，“阿尔茨海默病”的英文为“Alzheimer'sDisease”，以下简称“AD”，“基于最小绝对收缩和选择算子”的英文为“LeastAbsoluteShrinkageandSelectionOperator”，以下简称“LASSO”。本专利技术解决该技术问题所采用的技术方案是：多模态脑部神经影像特征的处理方法，是使用SWLRC的多模态神经影像特征选择方法进行生物标志物的挖掘，再使用多核SVM方法分类，具体步骤如下：第一步，多模态神经影像输入：对基于体素的形态测量法处理的磁共振脑影像和氟脱氧葡萄糖-正电子发射断层扫描脑影像这两种模态脑部影像数据提取感兴趣区模板特征，步骤是，首先使用统计参数映射软件包将基于体素的形态测量法处理的磁共振脑影像和氟脱氧葡萄糖-正电子发射断层扫描脑影像配准到标准空间中，然后采用116个感兴区的SPM软件的工具箱MarsBaR自动解剖标记模板分别将基于体素的形态测量法处理的磁共振脑影像的灰度密度数值和氟脱氧葡萄糖-正电子发射断层扫描脑影像的葡萄糖代谢水平的数值作为特征并进行提取，最后确定去除小脑后的90个感兴趣区的上述两种模态脑影像数据作为使用的特征，进一步对这些使用的特征进行标准化，由此完成多模态神经影像输入；第二步，使用SWLRC的多模态特征选择方法进行特征分析：将上述第一步中得到的每个模态的特征数据用表示，其中上标v表示第v个模态，下标i表示第i个神经影像，将每个神经影像所属的类别用yi表示，此处下标i同样表示第i个神经影像，然后将每个神经影像的数值和所属的类别输入到以下公式(1)所示的目标函数中，公式(1)中，m为神经影像个数，s为模态的个数，βv为第v个模态的特征选择向量，为第i个神经影像的第v个模态所占的权重为多少，且W＝[w1；w2；...；ws]∈Rm×s为训练样本的所有模态的权重矩阵，其中wv表示第v个模态的权重，表示第i个样本每个模态的权重，λS为约束特征稀疏的正则化参数，λR为约束样本多模态关联的正则化参数，R(W)表示W的秩，T表示矩阵和向量的转置，s.t.表示约束条件，由此完成使用SWLRC的多模态特征选择方法进行特征分析；第三本文档来自技高网...

【技术保护点】
1.多模态脑部神经影像特征的处理方法，其特征在于：是使用SWLRC的多模态神经影像特征选择方法进行生物标志物的挖掘，再使用多核SVM方法分类，具体步骤如下：第一步，多模态神经影像输入：对基于体素的形态测量法处理的磁共振脑影像和氟脱氧葡萄糖‑正电子发射断层扫描脑影像这两种模态脑部影像数据提取感兴趣区模板特征，步骤是，首先使用统计参数映射软件包将基于体素的形态测量法处理的磁共振脑影像和氟脱氧葡萄糖‑正电子发射断层扫描脑影像配准到标准空间中，然后采用116个感兴区的SPM软件的工具箱MarsBaR自动解剖标记模板分别将基于体素的形态测量法处理的磁共振脑影像的灰度密度数值和氟脱氧葡萄糖‑正电子发射断层扫描脑影像的葡萄糖代谢水平的数值作为特征并进行提取，最后确定去除小脑后的90个感兴趣区的上述两种模态脑影像数据作为使用的特征，进一步对这些使用的特征进行标准化，由此完成多模态神经影像输入；第二步，使用SWLRC的多模态特征选择方法进行特征分析：将上述第一步中得到的每个模态的特征数据用

【技术特征摘要】
1.多模态脑部神经影像特征的处理方法，其特征在于：是使用SWLRC的多模态神经影像特征选择方法进行生物标志物的挖掘，再使用多核SVM方法分类，具体步骤如下：第一步，多模态神经影像输入：对基于体素的形态测量法处理的磁共振脑影像和氟脱氧葡萄糖-正电子发射断层扫描脑影像这两种模态脑部影像数据提取感兴趣区模板特征，步骤是，首先使用统计参数映射软件包将基于体素的形态测量法处理的磁共振脑影像和氟脱氧葡萄糖-正电子发射断层扫描脑影像配准到标准空间中，然后采用116个感兴区的SPM软件的工具箱MarsBaR自动解剖标记模板分别将基于体素的形态测量法处理的磁共振脑影像的灰度密度数值和氟脱氧葡萄糖-正电子发射断层扫描脑影像的葡萄糖代谢水平的数值作为特征并进行提取，最后确定去除小脑后的90个感兴趣区的上述两种模态脑影像数据作为使用的特征，进一步对这些使用的特征进行标准化，由此完成多模态神经影像输入；第二步，使用SWLRC的多模态特征选择方法进行特征分析：将上述第一步中得到的每个模态的特征数据用表示，其中上标v表示第v个模态，下标i表示第i个神经影像，将每个神经影像所属的类别用yi表示，此处下标i同样表示第i个神经影像，然后将每个神经影像的数值和所属的类别输入到以下公式(1)所示的目标函数中，公式(1)中，m为神经影像个数，s为模态的个数，βv为第v个模态的特征选择向量，为第i个神经影像的第v个模态所占的权重为多少，且W＝[w1；w2；…；ws]∈Rm×s为训练样本的所有模态的权重矩阵，其中wv表示第v个模态的权重，表示第i个样本每个模态的权重，λS为约束特征稀疏的正则化参数，λR为约束样本多模态关联的正则化参数，R(W)表示W的秩，T表示矩阵和向量的转置，s.t.表示约束条件，由此完成使用SWLRC的多模态特征选择方法进行特征分析；第三步，优化目标函数：优化上述第二步中的公式(1)中的目标函数的方法是采用快速收敛的方案，即重复多次搜索公式(1)中的目标函数的局部最优解来近似最优解，并采用交替方向乘子法对上述公式(1)中的W和βv进行优化，目标函数被改写为如下公式(2)所示，s.t.-W+Λ＝0,J＝W,公式(2)中，为和βv的乘积，Λ∈Rm×s为非负松弛变量，J∈Rm×s等于W，‖·‖*表示星范数，其目的是让目标函数符合交替方向乘子法形式并进行求解，在线性规划中，将不等式约束条件变为等式约束条件，由此完成优化目标函数；第四步，目标函数的求解：为了表示简单，将上述第三步中的公式(2)的三项内容分别简记为f1(U)、f2(βv)、f3(J)，由此将公式(2)的增广拉格朗日形式定义为如下公式(3)，公式(3)中，Pv∈Rm×90，Q∈Rm×s，E∈Rm×s均为拉格朗日乘子，μ、ρ、ξ均为增广拉格朗日项的三个参数，‖·‖F表示F范数，Uv为由变量从1到m按行排列而构成的对应于不同模态的矩阵，wv为由构成的列向量，所有模态的矩阵Uv按列排列拼接而构成的跨模态矩阵U，<X1,X2>为矩阵X1TX2的迹，X1和X2表示为两个行数和列数相等的矩阵，因为Λ的解是解析解，即其中Qiv为第i个样本的第v个模态的拉格朗日乘子，将其带入公式(3)，消去Λ，则得到如下公式(4)，并且，将通用的交替方向乘子法框架写成如下公式(5)的形式：且在交替方向乘子法框架中，每一步都是凸优化问题，通过计算得到每一步的最...

【专利技术属性】
技术研发人员：郭迎春，包永进，郝小可，刘依，于洋，朱叶，师硕，阎刚，
申请(专利权)人：河北工业大学，
类型：发明
国别省市：天津,12