首页 / 专利库 / 诊断设备和程序 / 医学影像学 / 分子成像 / 磁共振成像 / 弥散张量成像 / 一种基于弥散张量成像技术的疾病诊断辅助系统及方法

一种基于弥散张量成像技术的疾病诊断辅助系统及方法

阅读:657发布:2020-05-19

专利汇可以提供一种基于弥散张量成像技术的疾病诊断辅助系统及方法专利检索,专利查询,专利分析的服务。并且本 发明 属于医学图像辅助诊断技术领域,尤其涉及一种基于 弥散张量成像 技术的 疾病 诊断辅助系统及方法。所述基于弥散张量成像技术的疾病诊断辅助系统包括图像预处理模 块 、专家 知识库 模块、张量学习模块和诊断结果输出模块;所述图像预处理模块用于弥散张量图像的采集、弥散张量图像重配准以及弥散张量图像特征提取;所述专家知识库模块用于建立大脑和脊髓相关疾病的专家知识库;所述张量学习模块用于利用专家知识库和张量学习 算法 对弥散张量图像进行张量分类模型训练以及张量分类模型优化和测试;所述诊断结果输出模块用于输出诊断结果。本发明充分挖掘图像的原有信息,提高模式分类 精度 ,极大地降低计算量,并有效地保整了疾病诊断的实时性。,下面是一种基于弥散张量成像技术的疾病诊断辅助系统及方法专利的具体信息内容。

1.一种基于弥散张量成像技术的疾病诊断辅助系统,其特征在于:包括图像预处理模、专家知识库模块、张量学习模块和诊断结果输出模块;所述图像预处理模块用于弥散张量图像的采集、弥散张量图像重配准以及弥散张量图像特征提取;所述专家知识库模块用于建立大脑和脊髓相关疾病的专家知识库;所述张量学习模块用于利用专家知识库和张量学习算法对弥散张量图像进行张量分类模型训练以及张量分类模型优化和测试;所述诊断结果输出模块用于输出诊断结果。
2.根据权利要求1所述的基于弥散张量成像技术的疾病诊断辅助系统,其特征在于,所述图像预处理模块包括图像采集单元、图像配准单元和特征提取单元;
所述图像采集单元用于采集弥散张量图像;
所述图像配准单元用于进行基于张量相似度的弥散张量图像配准,再进行基于标量相似度的弥散张量图像配准;
所述特征提取单元用于基于张量代数的多线性核主成分分析方法进行弥散张量图像的特征提取和降维
3.根据权利要求2所述的基于弥散张量成像技术的疾病诊断辅助系统,其特征在于,所述弥散张量图像的特征提取和降维具体为:首先选取适当的核函数把原始数据映射到各自的特征空间;再通过张量积将所有的特征映射到多重线性子空间,使得各子空间能捕捉到大部分正交多维的变异量;按照最小交替平方原理,计算求解新特征。
4.根据权利要求1所述的基于弥散张量成像技术的疾病诊断辅助系统,其特征在于,所述张量学习模块包括张量分类模型训练单元和张量分类模型优化单元;
所述张量分类模型训练单元用于利用专家知识库和最优投影支持张量机对弥散张量图像进行张量分类模型训练;
所述张量分类模型优化单元用于进行张量分类模型优化和测试。
5.根据权利要求4所述的基于弥散张量成像技术的疾病诊断辅助系统,其特征在于,所述张量分类模型训练具体为:基于最优投影支持张量机对弥散张量图像进行模式识别,寻找图像在边界方向上的最优投影,并通过最大化投影后样本的类间矩阵和类内矩阵的散布之比来计算最优的投影矢量,所述计算最优的投影矢量具体包括:计算类间矩阵Sb和类T T
内矩阵Sw散布之比的极大值max[(VSbV)/(VSwV)],求解第一个边界方向上的最优投影矢量v;把训练样本通过v投影,求解二次规划问题,同理求得第二个边界方向上的最优投影T T
矢量;用类间散布矩阵和类内散布矩阵定义投影离散系数R=Tr(Sb×Sb)/Tr(Sw×Sw),R代表训练样本经过最优投影矢量投影后样本间的分离度,R值越大表示训练样本经过投影后样本间距离就越大。
6.根据权利要求5所述的基于弥散张量成像技术的疾病诊断辅助系统,其特征在于,所述张量分类模型训练还包括:基于最优投影支持张量机分析两种投影方案为:每个子分类器的训练样本都向同一个方向投影,或者每个子分类器确定不同的投影矢量,然后训练样本分别向各自的投影矢量方向上投影。
7.一种基于弥散张量成像技术的疾病诊断辅助方法,包括:
步骤a:采集弥散张量图像,并进行弥散张量图像重配准以及弥散张量图像特征提取;
步骤b:建立大脑和脊髓相关疾病的专家知识库;
步骤c:利用专家知识库和张量学习算法对弥散张量图像进行张量分类模型训练以及张量分类模型优化和测试;
步骤d:输出诊断结果。
8.根据权利要求7所述的基于弥散张量成像技术的疾病诊断辅助方法,其特征在于,在所述步骤a中,所述弥散张量图像重配准包括基于张量相似度的弥散张量图像配准及基于标量相似度的弥散张量图像配准。
9.根据权利要求7所述的基于弥散张量成像技术的疾病诊断辅助方法,其特征在于,在所述步骤a中,所述弥散张量图像特征提取为:基于张量代数的多线性核主成分分析方法进行弥散张量图像的特征提取和降维;具体包括:首先选取适当的核函数把原始数据映射到各自的特征空间,再通过张量积将所有的特征映射到多重线性子空间,使得各子空间能捕捉到大部分正交多维的变异量,按照最小交替平方原理,计算求解新特征。
10.根据权利要求7所述的基于弥散张量成像技术的疾病诊断辅助方法,其特征在于,在所述步骤c中,所述张量分类模型训练具体为:基于最优投影支持张量机对弥散张量图像进行模式识别,寻找图像在边界方向上的最优投影,并通过最大化投影后样本的类间矩阵和类内矩阵的散布之比来计算最优的投影矢量,所述计算最优的投影矢量具体包T T
括:计算类间矩阵Sb和类内矩阵Sw散布之比的极大值max[(VSbV)/(VSwV)],求解第一个边界方向上的最优投影矢量v;把训练样本通过v投影,求解二次规划问题,同理求得第二个边界方向上的最优投影矢量;用类间散布矩阵和类内散布矩阵定义投影离散系数R=T T
Tr(Sb×Sb)/Tr(Sw×Sw),R代表训练样本经过最优投影矢量投影后样本间的分离度,R值越大表示训练样本经过投影后样本间距离就越大。

说明书全文

一种基于弥散张量成像技术的疾病诊断辅助系统及方法

技术领域

[0001] 本发明属于医学图像辅助诊断技术领域,尤其涉及一种基于弥散张量成像技术的疾病诊断辅助系统及方法。

背景技术

[0002] 弥散张量成像(Diffusion Tensor Imaging,DTI)技术是当今唯一无创伤的白质神经纤维束活体成像方法,是核磁共振成像(MRI)的特殊形式。不同于核磁共振成像是追踪分子中的氢原子,弥散张量成像是依据水分子移动方向制图。弥散张量成像图可以揭示脑瘤如何影响神经细胞连接,引导医疗人员进行大脑手术。它还可以揭示同中、多发性硬化症、精神分裂症阅读障碍等有关大脑和脊髓的细微反常变化。弥散张量成像数据本质上是二阶张量结构,它的每一个体素包含了水分子在白质神经纤维束内弥散的三维空间信息。
[0003] 通过机器学习方法能从原始弥散张量图像中提取出判断人体疾病的有效信息,从而为预测分析大脑和脊髓相关的疾病提供有的帮助。但是传统的机器学习方法都是基于向量模式的算法,例如支持向量机、线性判别分析以及神经网络等,这些算法要么就只处理弥散张量成像数据的一些标量指标,而不能充分利用弥散张量图像的结构空间信息;要么在分析处理之前,先将张量展开为向量。然而,这种做法会带来以下问题:1、破坏原始数据的结构和张量结构数据的内在相关性。2、大脑弥散张量成像数据在向量化后,其生成的统计参数(典型的为协方差阵)的维数极大,会破坏原始数据的结构,丢失数据的内在相关性,从而导致极高的计算复杂度和存储代价。

发明内容

[0004] 本发明提供了一种基于弥散张量成像技术的疾病诊断辅助系统及方法,旨在解决现有的向量模式学习算法不能充分利用弥散张量图像数据的结构空间信息,且在把张量数据向量化的过程中,会破坏原始数据的结构,丢失数据的内在相关性,增加计算复杂度和存储代价的技术问题。
[0005] 本发明是这样实现的,一种基于弥散张量成像技术的疾病诊断辅助系统,包括图像预处理模、专家知识库模块、张量学习模块和诊断结果输出模块;所述图像预处理模块用于弥散张量图像的采集、弥散张量图像重配准以及弥散张量图像特征提取;所述专家知识库模块用于建立大脑和脊髓相关疾病的专家知识库;所述张量学习模块用于利用专家知识库和张量学习算法对弥散张量图像进行张量分类模型训练以及张量分类模型优化和测试;所述诊断结果输出模块用于输出诊断结果。
[0006] 本发明实施例采取的技术方案还包括:所述图像预处理模块包括图像采集单元、图像配准单元和特征提取单元;
[0007] 所述图像采集单元用于采集弥散张量图像;
[0008] 所述图像配准单元用于进行基于张量相似度的弥散张量图像配准,再进行基于标量相似度的弥散张量图像配准;
[0009] 所述特征提取单元用于基于张量代数的多线性核主成分分析方法进行弥散张量图像的特征提取和降维
[0010] 本发明实施例采取的技术方案还包括:所述弥散张量图像的特征提取和降维具体为:首先选取适当的核函数把原始数据映射到各自的特征空间;再通过张量积将所有的特征映射到多重线性子空间,使得各子空间能捕捉到大部分正交多维的变异量;按照最小交替平方原理,计算求解新特征。
[0011] 本发明实施例采取的技术方案还包括:所述张量学习模块包括张量分类模型训练单元和张量分类模型优化单元;
[0012] 所述张量分类模型训练单元用于利用专家知识库和最优投影支持张量机对弥散张量图像进行张量分类模型训练;
[0013] 所述张量分类模型优化单元用于进行张量分类模型优化和测试。
[0014] 本发明实施例采取的技术方案还包括:所述张量分类模型训练具体为:基于最优投影支持张量机对弥散张量图像进行模式识别,寻找图像在边界方向上的最优投影,并通过最大化投影后样本的类间矩阵和类内矩阵的散布之比来计算最优的投影矢量,所述计算T最优的投影矢量具体包括:计算类间矩阵Sb和类内矩阵Sw散布之比的极大值max[(VSbV)/T
(VSwV)],求解第一个边界方向上的最优投影矢量v;把训练样本通过v投影,求解二次规划问题,同理求得第二个边界方向上的最优投影矢量;用类间散布矩阵和类内散布矩阵定义T T
投影离散系数R=Tr(Sb×Sb)/Tr(Sw×Sw),R代表训练样本经过最优投影矢量投影后样本间的分离度,R值越大表示训练样本经过投影后样本间距离就越大。
[0015] 本发明实施例采取的技术方案还包括:所述张量分类模型训练还包括:基于最优投影支持张量机分析两种投影方案为:每个子分类器的训练样本都向同一个方向投影,或者每个子分类器确定不同的投影矢量,然后训练样本分别向各自的投影矢量方向上投影。
[0016] 本发明实施例采取的另一技术方案为:一种基于弥散张量成像技术的疾病诊断辅助方法,包括:
[0017] 步骤a:采集弥散张量图像,并进行弥散张量图像重配准以及弥散张量图像特征提取;
[0018] 步骤b:建立大脑和脊髓相关疾病的专家知识库;
[0019] 步骤c:利用专家知识库和张量学习算法对弥散张量图像进行张量分类模型训练以及张量分类模型优化和测试;
[0020] 步骤d:输出诊断结果。
[0021] 本发明实施例采取的技术方案还包括:在所述步骤a中,所述弥散张量图像重配准包括基于张量相似度的弥散张量图像配准及基于标量相似度的弥散张量图像配准。
[0022] 本发明实施例采取的技术方案还包括:在所述步骤a中,所述弥散张量图像特征提取为:基于张量代数的多线性核主成分分析方法进行弥散张量图像的特征提取和降维;具体包括:首先选取适当的核函数把原始数据映射到各自的特征空间,再通过张量积将所有的特征映射到多重线性子空间,使得各子空间能捕捉到大部分正交多维的变异量,按照最小交替平方原理,计算求解新特征。
[0023] 本发明实施例采取的技术方案还包括:在所述步骤c中,所述张量分类模型训练具体为:基于最优投影支持张量机对弥散张量图像进行模式识别,寻找图像在边界方向上的最优投影,并通过最大化投影后样本的类间矩阵和类内矩阵的散布之比来计算最优的投影矢量,所述计算最优的投影矢量具体包括:计算类间矩阵Sb和类内矩阵Sw散布之比的极T T大值max[(VSbV)/(VSwV)],求解第一个边界方向上的最优投影矢量v;把训练样本通过v投影,求解二次规划问题,同理求得第二个边界方向上的最优投影矢量;用类间散布矩阵和T T
类内散布矩阵定义投影离散系数R=Tr(Sb×Sb)/Tr(Sw×Sw),R代表训练样本经过最优投影矢量投影后样本间的分离度,R值越大表示训练样本经过投影后样本间距离就越大。
[0024] 本发明实施例的基于弥散张量成像技术的疾病诊断辅助系统及方法针对弥散张量图像中方向信息比灰度信息更敏感的特点,提出二重度量配准方法;其次,提出一种多线性核主成分分析方法实现弥散张量图像的特征提取和降维;最后,针对向量模式学习算法在处理张量结构数据方面的局限性,结合分类器的结构风险最小化原理和现有的张量模式学习算法,并基于最优投影矢量准则,提出基于最优投影支持张量机,实现基于张量模式算法和弥散张量图像的疾病诊断实时识别,可以在不破坏原始数据结构的前提下,充分挖掘图像的原有信息,提高模式分类精度,可以极大地降低计算量和计算成本,并有效地保整了疾病诊断的实时性。本发明可用于在星形细胞瘤间变型、精神脱髓鞘性病变、中枢神经系统发育异常、精神分裂症、抑郁症、脊髓型颈椎病等与脑和脊髓相关联的疾病分析和预测方面,同时本发明也可以用于正常人的大脑和脊髓相关疾病的预防和健康护理方面。附图说明
[0025] 图1是本发明实施例的基于弥散张量成像技术的疾病诊断辅助系统的结构示意图;
[0026] 图2是本发明实施例的基于弥散张量成像技术的疾病诊断辅助方法的流程图

具体实施方式

[0027] 为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
[0028] 请参阅图1,是本发明实施例的基于弥散张量成像技术的疾病诊断辅助系统的结构示意图。本发明实施例的基于弥散张量成像技术的疾病诊断辅助系统包括图像预处理模块、专家知识库模块、张量学习模块和诊断结果输出模块;其中,图像预处理模块用于弥散张量图像的采集、基于张量模式的弥散张量图像重配准以及基于多线性核主成分分析方法的弥散张量图像特征提取;专家知识库模块用于建立大脑和脊髓相关疾病的专家知识库;张量学习模块用于利用专家知识库和张量学习算法对弥散张量图像进行张量分类模型训练及张量分类模型优化和测试;诊断结果输出模块用于输出诊断结果。
[0029] 具体地,图像预处理模块包括图像采集单元、图像配准单元和特征提取单元,[0030] 图像采集单元用于采集弥散张量图像;
[0031] 图像配准单元用于针对弥散张量图像中方向信息比灰度信息更敏感的特点,首先进行基于张量相似度的弥散张量图像配准,再进行基于标量相似度的弥散张量图像配准。
[0032] 特征提取单元用于基于张量代数的多线性核主成分分析方法实现弥散张量图像的特征提取和降维,具体为:首先选取适当的核函数把原始数据映射到各自的特征空间,再通过张量积将所有的特征映射到多重线性子空间,使得各子空间能捕捉到大部分正交多维的变异量,按照最小交替平方原理,计算求解新特征。其中,主成分分析(Principal Component Analysis,PCA)方法是将多个变量通过线性变换以选出较少个数重要变量的一种多元统计分析方法,又称主分量分析。
[0033] 张量学习模块包括张量分类模型训练单元和张量分类模型优化单元,[0034] 张量分类模型训练单元用于利用专家知识库和最优投影支持张量机(OPSTM)对弥散张量图像进行张量分类模型训练;具体为:利用Fisher准则,基于最优投影支持张量机对弥散张量图像进行模式识别,寻找图像在边界方向上的最优投影,并通过最大化投影后样本的类间矩阵和类内矩阵的散布之比来计算最优的投影矢量;算法分为两步来实现:T T
(1)、计算类间矩阵Sb和类内矩阵Sw散布之比的极大值max[(VSbV)/(VSwV)],求解第一个边界方向上的最优投影矢量v。(2)、把训练样本通过v投影,求解二次规划问题,同理可求得第二个边界方向上的最优投影矢量。本发明利用类间散布矩阵和类内散布矩阵定义投影T T
离散系数R=Tr(Sb×Sb)/Tr(Sw×Sw),R代表训练样本经过最优投影矢量投影后样本间的分离度,R值越大表示训练样本经过投影后样本间距离就越大。将基于最优投影支持张量机从二分类推广到多分类中,分析两种投影方案:1、每个子分类器的训练样本都向同一个方向投影。2、每个子分类器确定不同的投影矢量,然后训练样本分别向各自的投影矢量方向上投影。由于把训练样本投影到张量子空间中,不同类的训练样本会相互分离,同类的训练样本会相互靠近。这使得选取同一投影矢量的方案识别率高于选取非同一投影矢量的方案。
[0035] 张量分类模型优化单元用于进行张量分类模型优化和测试;具体为:结合多线性代数运算规则,针对流数据特点,把最优投影支持张量机推广到在线形式,并将优化问题中的不等式约束转化为等式约束,对弥散张量成像数据(二阶张量)来说,仅需要求解两个线性方程组,进行简单的矩阵逆运算,可以大幅降低计算量。在目标函数优化方法上,充分考虑流数据的特点,利用随机梯度下降法完成拉格朗日乘子迭代。为了能保持流数据的原始结构信息并降低计算成本,利用张量的CP分解辅助张量内积运算,这可以在保持原始张量的自然结构信息的同时,减少计算量并节省存储空间。
[0036] 请参阅图2,是本发明实施例的基于弥散张量成像技术的疾病诊断辅助方法的流程图。本发明实施例的基于弥散张量成像技术的疾病诊断辅助方法包括以下步骤:
[0037] 步骤100:采集弥散张量图像;
[0038] 步骤200:对所述采集的弥散张量图像进行配准。具体而言:
[0039] 针对弥散张量图像中方向信息比灰度信息更敏感的特点,首先进行基于张量相似度的弥散张量图像配准,再进行基于标量相似度的弥散张量图像配准
[0040] 步骤300:基于张量代数的多线性核主成分分析方法进行弥散张量图像的特征提取和降维;
[0041] 在步骤300中,弥散张量图像的特征提取和降维具体为:首先选取适当的核函数把原始数据映射到各自的特征空间,再通过张量积将所有的特征映射到多重线性子空间,使得各子空间能捕捉到大部分正交多维的变异量,按照最小交替平方原理,计算求解新特征。
[0042] 步骤400:建立大脑和脊髓相关疾病的专家知识库;
[0043] 步骤500:寻找上述弥散张量图像在边界方向上的最优投影,并通过最大化投影后样本的类间矩阵和类内矩阵的散布之比来计算最优的投影矢量。具体而言:
[0044] 分为两步来实现:(1)、计算类间矩阵Sb和类内矩阵Sw散布之比的极大值T Tmax[(VSbV)/(VSwV)],求解第一个边界方向上的最优投影矢量v。(2)、把训练样本通过v投影,求解二次规划问题,同理可求得第二个边界方向上的最优投影矢量。
[0045] 步骤600:利用专家知识库和最优投影支持张量机进行张量分类模型训练;
[0046] 在步骤600中,张量分类模型训练具体为:利用Fisher准则,基于最优投影支持张量机对弥散张量图像进行模式识别。本发明利用类间散布矩阵和类内散布矩阵定义投影离T T散系数R=Tr(Sb×Sb)/Tr(Sw×Sw),R代表训练样本经过最优投影矢量投影后样本间的分离度,R值越大表示训练样本经过投影后样本间距离就越大。将基于最优投影支持张量机从二分类推广到多分类中,分析两种投影方案:1、每个子分类器的训练样本都向同一个方向投影。2、每个子分类器确定不同的投影矢量,然后训练样本分别向各自的投影矢量方向上投影。由于把训练样本投影到张量子空间中,不同类的训练样本会相互分离,同类的训练样本会相互靠近。这使得选取同一投影矢量的方案识别率高于选取非同一投影矢量的方案。
[0047] 步骤700:进行张量分类模型优化和测试;
[0048] 在步骤700中,张量分类模型优化和测试具体为:结合多线性代数运算规则,针对流数据特点,把最优投影支持张量机推广到在线形式,并将优化问题中的不等式约束转化为等式约束,对弥散张量成像数据(二阶张量)来说,仅需要求解两个线性方程组,进行简单的矩阵逆运算,可以大幅降低计算量。在目标函数优化方法上,充分考虑流数据的特点,利用随机梯度下降法完成拉格朗日乘子迭代。为了能保持流数据的原始结构信息并降低计算成本,利用张量的CP分解辅助张量内积运算,这可以在保持原始张量的自然结构信息的同时,减少计算量并节省存储空间。
[0049] 步骤800:输出诊断结果。
[0050] 本发明实施例的基于弥散张量成像技术的疾病诊断辅助系统及方法针对弥散张量图像中方向信息比灰度信息更敏感的特点,提出二重度量配准方法;其次,提出一种多线性核主成分分析方法实现弥散张量图像的特征提取和降维;最后,针对向量模式学习算法在处理张量结构数据方面的局限性,结合分类器的结构风险最小化原理和现有的张量模式学习算法,并基于最优投影矢量准则,提出基于最优投影支持张量机,实现基于张量模式算法和弥散张量图像的疾病诊断实时识别,可以在不破坏原始数据结构的前提下,充分挖掘图像的原有信息,提高模式分类精度,可以极大地降低计算量和计算成本,并有效地保整了疾病诊断的实时性。本发明可用于在星形细胞瘤间变型、精神脱髓鞘性病变、中枢神经系统发育异常、精神分裂症、抑郁症、脊髓型颈椎病等与脑和脊髓相关联的疾病分析和预测方面,同时本发明也可以用于正常人的大脑和脊髓相关疾病的预防和健康护理方面。
[0051] 以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈