专利摘要

专利摘要

本发明公开了一种基于脑磁共振影像的Aβ蛋白沉积信息检测系统，设置有MRI图像采集装置、图像预处理装置、特征提取装置以及参数反演装置；所述特征提取装置用于提取最优特征子集；所述参数反演装置用于将所述最优特征子集反演出Aβ蛋白沉积含量；该参数反演装置中固化有训练好的SVM模型，所述SVM模型是通过选择多个训练样本图像的多个特征参数，并以每个样本图像对应PET图像的Aβ蛋白沉积含量为评价标准，最终训练出的最优特征子集与Aβ蛋白沉积含量的映射关系。其显著效果是：不仅具有无创、无辐射、安全、自动化程度高等优点，而且可应用于阿尔茨海默病诊断的临床应用。

权利要求

1.一种基于脑磁共振影像的Aβ蛋白沉积信息检测系统，其特征在于：设置有MRI图像采集装置、图像预处理装置、特征提取装置以及参数反演装置；

所述图像预处理装置中依次设置有滤波去噪装置、图像配准装置、头骨剥离装置和组织分割装置；

所述滤波去噪装置用于实现MRI图像的去噪处理；

所述图像配准装置用于实现待测图像头部位置的调整；

所述头骨剥离装置用于去除图像中的脑外组织；

所述组织分割装置用于将图像中的脑组织图像分割为脑新皮层、白质、灰质、丘脑和海马体五个解剖结构；

所述特征提取装置用于提取最优特征子集；

所述参数反演装置用于将所述最优特征子集反演出Aβ蛋白沉积含量；该参数反演装置中固化有训练好的SVM模型，所述SVM模型是通过选择多个训练样本图像的多个特征参数，并以每个样本图像对应PET图像的Aβ蛋白沉积含量为评价标准，最终训练出的最优特征子集与Aβ蛋白沉积含量的映射关系。

2.根据权利要求1所述的基于脑磁共振影像的Aβ蛋白沉积信息检测系统，其特征在于：所述参数反演装置中设置有训练样本图像获取装置、训练样本图像预处理装置、训练样本图像特征提取装置、训练样本特征选择装置、测试样本图像获取装置、测试样本图像预处理装置、测试样本图像特征提取装置以及分类器构建装置；

所述训练样本图像获取装置用于获取训练样本MRI图像和该MRI图像对应的PET图像，并通过PET图像计算Aβ蛋白沉积含量作为分类器的评价标准；

所述训练样本图像预处理装置与所述图像预处理装置的结构和功能相同；

所述训练样本图像特征提取装置用于提取SVM模型所需的多个特征参数；

所述训练样本特征选择装置用于实现多个特征参数的选择；

所述测试样本图像获取装置用于获取测试样本MRI图像和该MRI图像对应的PET图像，并通过PET图像计算Aβ蛋白沉积含量作为分类器的评价标准；

所述测试样本图像预处理装置与所述图像预处理装置的结构和功能相同；

所述测试样本图像特征提取装置用于提取训练样本特征选择装置所选择出的特征参数；

所述分类器构建装置用于构建特征参数与Aβ蛋白沉积含量的反演关系，该分类器的核函数采用径向基函数，当所述测试样本的识别率达到预设目标时，最终训练出的最优特征子集与Aβ蛋白沉积含量的映射关系。

3.根据权利要求2所述的基于脑磁共振影像的Aβ蛋白沉积信息检测系统，其特征在于：所述训练样本图像特征提取装置所提取的多个特征参数包括每个解剖结构的体积、灰度分布的不均匀性、灰度平均、灰度均方差、灰度熵、信息熵、几何矩、对比度、差分矩、自相关、梯度分布的不均匀性、梯度、梯度平均、梯度均方差、梯度熵、逆差矩、平均弥散率和各向异性分数值。

4.根据权利要求2所述的基于脑磁共振影像的Aβ蛋白沉积信息检测系统，其特征在于：所述训练样本图像特征提取装置采用链式智能体遗传算法进行特征选择。

说明书

技术领域

本发明涉及到一种疾病辅助诊断系统，具体地说，是一种基于脑磁共振影像的Aβ蛋白沉积信息检测系统。

背景技术

阿尔茨海默病的影像学和图像特征提取研究表明，提取图像特征来表达阿尔茨海默病的影像学特征是可行和有效的。Aβ蛋白沉积含量是诊断阿尔茨海默病的一个重要指标。脑磁共振图像中灰度特征可以表达明暗等影像学特征，虽然不同序列图像灰度变化较大，但就同一个序列产生的图像来说，灰度对比能较好的分辨Aβ蛋白沉积区域；形状特征可以表达体积、厚度、形态等影像学特征；纹理特征可以表达方向、对比度、粗细、杂乱、密集等影像学特征，图像特征可以有效表征不同成分或解剖结构的影像学特征。因此，提取脑新皮层、白质、灰质、丘脑、海马体等相关解剖结构及纤维束的灰度特征、形状特征、纹理特征可以将与Aβ蛋白沉积相关的信息充分表达出来，即有效显示Aβ蛋白沉积含量信息。

然而，由于图像特征数量通常较多，会增加后续模式分类或参数反演的时间和空间复杂度，也不利于临床应用。此外，特征间具有一定的排斥性和相关性，会影响后续处理的精度。另外，由于待选特征较多，特征与Aβ蛋白沉积含量的关系未明且复杂，因此，Aβ蛋白沉积分布及含量的识别分类精度不高。

发明内容

针对现有技术的不足，本发明的目的是提供一种基于脑磁共振影像的Aβ蛋白沉积信息检测系统，通过对脑磁共振图像进行处理，提取显示相关影像学特征，并通过特征选择和机器学习，利用脑磁共振影像中的图像特征信息来实现定量检测Aβ蛋白沉积分布及含量，识别分类精度高。

为达到上述目的，本发明的技术方案如下：

一种基于脑磁共振影像的Aβ蛋白沉积信息检测系统，设置有MRI图像采集装置、图像预处理装置、特征提取装置以及参数反演装置；

所述图像预处理装置中依次设置有滤波去噪装置、图像配准装置、头骨剥离装置和组织分割装置；

所述滤波去噪装置用于实现MRI图像的去噪处理；

所述图像配准装置用于实现待测图像头部位置的调整；

所述头骨剥离装置用于去除图像中的脑外组织；

所述组织分割装置用于将图像中的脑组织图像分割为脑新皮层、白质、灰质、丘脑和海马体五个解剖结构；

所述特征提取装置用于提取最优特征子集；

所述参数反演装置用于将所述最优特征子集反演出Aβ蛋白沉积含量；该参数反演装置中固化有训练好的SVM模型，所述SVM模型是通过选择多个训练样本图像的多个特征参数，并以每个样本图像对应PET图像的Aβ蛋白沉积含量为评价标准，最终训练出的最优特征子集与Aβ蛋白沉积含量的映射关系。

进一步的技术方案是，所述参数反演装置中设置有训练样本图像获取装置、训练样本图像预处理装置、训练样本图像特征提取装置、训练样本特征选择装置、测试样本图像获取装置、测试样本图像预处理装置、测试样本图像特征提取装置以及分类器构建装置；

所述训练样本图像获取装置用于获取训练样本MRI图像和该MRI图像对应的PET图像，并通过PET图像计算Aβ蛋白沉积含量作为分类器的评价标准；

所述训练样本图像预处理装置与所述图像预处理装置的结构和功能相同；

所述训练样本图像特征提取装置用于提取SVM模型所需的多个特征参数；

所述训练样本特征选择装置用于实现多个特征参数的选择；

所述测试样本图像获取装置用于获取测试样本MRI图像和该MRI图像对应的PET图像，并通过PET图像计算Aβ蛋白沉积含量作为分类器的评价标准；

所述测试样本图像预处理装置与所述图像预处理装置的结构和功能相同；

所述测试样本图像特征提取装置用于提取训练样本特征选择装置所选择出的特征参数；

所述分类器构建装置用于构建特征参数与Aβ蛋白沉积含量的反演关系，该分类器的核函数采用径向基函数，当所述测试样本的识别率达到预设目标时，最终训练出的最优特征子集与Aβ蛋白沉积含量的映射关系。

进一步的技术方案是，所述训练样本图像特征提取装置所提取的多个特征参数包括每个解剖结构的体积、灰度分布的不均匀性、灰度平均、灰度均方差、灰度熵、信息熵、几何矩、对比度、差分矩、自相关、梯度分布的不均匀性、梯度、梯度平均、梯度均方差、梯度熵、逆差矩、平均弥散率和各向异性分数值。

作为优选，所述训练样本图像特征提取装置采用链式智能体遗传算法进行特征选择。

本发明的显著效果是：通过智能遗传算法对SVM模型进行训练，从而得到最优特征子集与Aβ蛋白沉积分布及含量的映射关系，将训练好的SVM模型固化到参数反演装置中，从而构建出本发明所提出的基于脑磁共振影像中磷酸化tau蛋白含量信息的检测系统，不仅具有无创、无辐射、安全、自动化程度高等优点，而且还可应用于阿尔茨海默病诊断的临床应用，可有助于提高诊断的准确性和实时性。

附图说明

图1是本发明的原理结构框图。

具体实施方式

下面结合附图对本发明的具体实施方式以及工作原理作进一步详细说明。

如图1所示，一种基于脑磁共振影像的Aβ蛋白沉积信息检测系统，设置有MRI图像采集装置、图像预处理装置、特征提取装置以及参数反演装置；

所述图像预处理装置中依次设置有滤波去噪装置、图像配准装置、头骨剥离装置和组织分割装置；

所述滤波去噪装置用于实现MRI图像的去噪处理；

所述图像配准装置用于实现待测图像头部位置的调整；

由于每张图像中包括若干张人脑图像，为了使所有的图像的头部在图像中位置和方向上达到一致，需要对图像进行配准。本发明采用基于属性向量的弹性配准(Hierarchical Attribute Matching Mechanism for Elastic Registration，HAMMER)算法对图像样本进行弹性配准，采用互信息值和相关系数作为弹性配准满意的评价准则。

所述头骨剥离装置用于去除图像中的脑外组织；

由于脑MRI图像中脑新皮层、白质、灰质、丘脑、海马体中明显存在磷酸化tau蛋白，亮度较暗，且体积和纹理有变化，因此组织分割装置用于将图像中的脑组织图像分割为脑新皮层、白质、灰质、丘脑和海马体五个解剖结构，具体实施时可以采用ITK程序包(Insight Segmentation and Registration Toolkit)对以上解剖结构进行分割。

所述特征提取装置用于提取最优特征子集；

所述参数反演装置用于将所述最优特征子集反演出Aβ蛋白沉积含量；该参数反演装置中固化有训练好的SVM模型，所述SVM模型是通过选择多个训练样本图像的多个特征参数，并以每个样本图像对应PET图像的Aβ蛋白沉积含量为评价标准，最终训练出的最优特征子集与Aβ蛋白沉积含量的映射关系。

在具体实施过程中，所述参数反演装置中设置有训练样本图像获取装置、训练样本图像预处理装置、训练样本图像特征提取装置、训练样本特征选择装置、测试样本图像获取装置、测试样本图像预处理装置、测试样本图像特征提取装置以及分类器构建装置；

所述训练样本图像获取装置用于获取训练样本MRI图像和该MRI图像对应的PET图像，并通过PET图像计算Aβ蛋白沉积含量作为分类器的评价标准；

所述训练样本图像预处理装置与所述图像预处理装置的结构和功能相同；

所述训练样本图像特征提取装置用于提取SVM模型所需的多个特征参数；

所述训练样本特征选择装置用于实现多个特征参数的选择；

所述测试样本图像获取装置用于获取测试样本MRI图像和该MRI图像对应的PET图像，并通过PET图像计算Aβ蛋白沉积含量作为分类器的评价标准；

所述测试样本图像预处理装置与所述图像预处理装置的结构和功能相同；

所述测试样本图像特征提取装置用于提取训练样本特征选择装置所选择出的特征参数；

所述分类器构建装置用于构建特征参数与Aβ蛋白沉积含量的反演关系，该分类器的核函数采用径向基函数，当所述测试样本的识别率达到预设目标时，最终训练出的最优特征子集与Aβ蛋白沉积含量的映射关系。

所述训练样本图像特征提取装置所提取的多个特征参数包括每个解剖结构的体积、灰度分布的不均匀性、灰度平均、灰度均方差、灰度熵、信息熵、几何矩、对比度、差分矩、自相关、梯度分布的不均匀性、梯度、梯度平均、梯度均方差、梯度熵、逆差矩、平均弥散率和各向异性分数值。其中，第1个特征表明Aβ蛋白沉积对体积的影响；第2-4个特征表明Aβ蛋白沉积对亮度的影响；第6-16个特征表明Aβ蛋白沉积对纹理的影响，第17-18个特征为纤维束特性。因此，同样一幅MR脑组织图像中的特征就有5×18＝90种。将MRI图像纹理特征记为f_t(j＝1,2,p)，形状特征记为f_s(i＝1,2,m)。

在SVM模型训练过程中，可以通过合并纹理特征f_t(j＝1,2,p)和形状特征f_s(i＝1,2,m)，并基于合并后的特征和SVM回归算法，建立反演模型，反演过程中将Aβ蛋白沉积含量记为Count_Aβ，并以此为评价标准，通过反演准确性来进行特征选择，最终获得最优特征子集，记为f_j(j＝1,2,p)，反演得到最接近Count_Aβ的Aβ蛋白沉积含量Count_Aβ’以及回归模型SVM_final，优选的f_j(j＝1,2,p)以及SVM_final，最后以SVM_final作为训练好的SVM模型，通过提取其它MRI图像中的最优特征子集即可反演出该MRI图像对应的Aβ蛋白沉积含量。

在实施过程中，采集40例MR图像样本和对应的PET图像样本(正常人脑的图像和确诊老年痴呆人脑的图像各一半)，选择一幅图像作为参考图像，将其余图像向其配准；然后通过头骨剥离，去除图像中的非脑组织；接着进行组织分割，从脑组织中分离出脑新皮层、白质、灰质、丘脑、海马体等几部分；最后从每个组织中提取出体积、灰度分布的不均匀性、灰度平均、灰度均方差、灰度熵、信息熵、几何矩、对比度、差分矩、自相关、梯度分布的不均匀性、梯度、梯度平均、梯度均方差、梯度熵、逆差矩、平均弥散率和各向异性分数值等18个特征，即每幅图包含5*18＝90个特征。

将MRI图像样本和对应的PET图像样本随机分为A，B两组，每组包含20例MR图像和对应的20例PET图像(正常人脑的图像和确诊老年痴呆人脑的图像各一半)，A组用于参数反演模型训练，B组用于对模型进行性能测试。

基于SVM的参数反演模型，采用A部分作为训练样本，B部分作为测试样本。通过反演准确性来进行特征筛选，最终获得最优特征子集f_j(j＝1,2,p)，以及反演得到最接近Count_Aβ的Aβ蛋白沉积含量Count_Aβ’。特征选择采用链式智能体遗传算法，用二进制编码表示待选特征。由于待选特征为90个，则设置遗传基因长度为90的2进制字符串，数字为‘1’表示选中该特征，‘0’表示未选中，选出的特征作为分类器的输入。种群数量选择150，迭代次数为100次，自适应交叉概率采用如下公式为：

pc=(pc1-pc2)(f′-favg)fmax-favg,f′≥favgpc1,f′<favg]]>

式中，f_avg表示每代种群的平均适应度，f_max表示每代种群的最大适应度，f'表示待交叉的两个个体中较大的适应度值；初始化p_c1＝0.9，p_c2＝0.6,交叉操作采用自适应交叉概率的单点交叉法。

同理，基因变异同样采用自适应的变异概率，其计算表达式如下：

pm=pm1-(pm1-pm2)(fmax-f)fmax-favg,f≥favgpm1,f<favg]]>

式中，f_avg表示每代种群的平均适应度，f_max表示表示每代种群的最大适应度，f表示待变异个体的适应度值；初始化p_m1＝0.1，p_m2＝0.006，变异操作采用自适应变异概率的二进制变异法。

SVM核函数为径向基函数，采用5阶校验法，训练收敛准则为均方误差(Mean squared normalized error，MSE)。利用A部分样本数据进行测试，对SVM进行参数回归，输入向量为图像特征值，输出为Aβ蛋白沉积含量标准值，均方误差满足要求停止，获取参数矩阵，最优特征子集与Aβ蛋白沉积含量的映射关系。

在训练过程中，由于待选特征较多，特征与Aβ蛋白沉积含量的关系未明且复杂，因此，通过同时优化待选图像特征和SVM回归模型参数，来提高图像特征选择和获取与Aβ蛋白沉积含量关系式的精度，结合链式智能体遗传算法对训练样本进行特征提取，以及采用二进制编码与十进制编码混合编码制来表示SVM回归模型参数，保证了后续处理的精度。

基于脑磁共振影像的Aβ蛋白沉积信息检测系统专利购买费用说明

Q：办理专利转让的流程及所需资料

A：专利权人变更需要办理著录项目变更手续，有代理机构的，变更手续应当由代理机构办理。

1：专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。

2：按规定缴纳著录项目变更手续费。

3：同时提交相关证明文件原件。

4：专利权转移的，变更后的专利权人委托新专利代理机构的，应当提交变更后的全体专利申请人签字或者盖章的委托书。

Q:专利著录项目变更费用如何缴交

A：（1）直接到国家知识产权局受理大厅收费窗口缴纳,（2）通过代办处缴纳，（3）通过邮局或者银行汇款,更多缴纳方式

Q：专利转让变更，多久能出结果

A：著录项目变更请求书递交后，一般1-2个月左右就会收到通知，国家知识产权局会下达《转让手续合格通知书》。