一种基于多变量的DTI图像分析方法技术

本发明专利技术公开了一种基于多变量的DTI图像分析方法，用于疾病大脑机制研究中对脑部白质病变区域确认和提取，从而为临床治疗提供影像学证据，具体步骤为：对数据进行预处理，包括去涡流、头动矫正、去脑壳、弥散拟合、构建白质骨架；对预处理后的数据进行特征提取,通过置换检验并以年龄为协变量得到病人组和正常人组的显著差异区域，并分别计算显著差异区域中特定变量的平均值得到特征值；用留一法进行交叉验证，判断是否满足停止条件，若不满足则计算每个特征的平均权值，并去除平均权值最小的特征，直到满足停止条件为止；最后得到的脑部区就是脑部病变区域。本发明专利技术利用DTI成像模式，通过多变量的研究方法，为临床治疗找到了病变区域提供了影像学证据。

【技术实现步骤摘要】
一种基于多变量的DTI图像分析方法
本专利技术涉及一种图像分析方法，具体涉及一种基于多变量的DTI图像分析方法。
技术介绍
扩散张量成像(DTI)是一种非侵入性的能提供活体内水分子扩散运动的成像技术，它能检测传统MRI所不能观察到的组织的微观变化，是MR成像技术的重大突破。基于脑影像信息的模式分类是目前脑影像研究中的热点课题。利用图像分类方法，计算DTI图像具有某种属性的可能性大小，或者自动的判别图像的类别属性，是计算机辅助分析的一个重要应用。文献“AlexanderAL,LeeJE,etal.DiffusiontensorimagingofthecorpuscallosuminAutism.Neuroimage.2007；34(1):61–73.”使用基于感兴趣区域的方法研究胼胝体在孤独症中的作用，但是这种方法需要有关某种病理或者病变区域的先验知识，所以没有很好地推广性。文献“RidgwayGR,HenleySM,etal.Tensimplerulesforreportingvoxel-basedmorphometrystudies.Neuroimage.2008；40(4):1429–1435.”使用基于体素的双样本t检验统计方法研究病人和正常人之间的组间差异。这种方法是在假设变量符合正态分布的情况下仅仅考虑单个变量的影响没有同时考虑多个变量的影响，而且无法判断单个人是否是病人或者是正常人。文献“MadhuraIngalhalikar,etal.DiffusionbasedAbnormalityMarkersofPathology:TowardsLearnedDiagnosticPredictionofASD.Neuroimage.2011；57(3):918–927”将地图集中的各向异性(FA)和平均扩散率(MD)的值作为特征，然后将所得特征加入支持向量机(SVM)中，通过留一法不断地选择可以使支持向量机得到最好的正确率和推广性的特征。然而，这种基于地图集的方法无法提取地图集下子区域的相关变量作为特征，这样就无法找到子区域中病变的区域。同时，这种方法没有考虑年龄因素对白质的影响。
技术实现思路
针对现有技术的不足，本专利技术旨在提供一种推广性更强的基于多变量的DTI图像分析法，本专利技术首先经过预处理和弥散张量计算得到各向异性值(FA)、平均扩散率值(MD)、径向弥散率值(RD)，针对计算得到的三个量进行置换检验后得到病人与正常人间的组间差异，然后将差异区域的各向异性值、平均扩散率值、径向弥散率值分别平均后加入线性支持向量机中用留一法训练支持向量机，直至找到使分类器准确率高、推广性好的特征量，最终得到的这些特征量所在的区域就很有可能是病变区域。为了实现上述目的，本专利技术采用如下技术方案：一种基于多变量的DTI图像分析方法，包括图像数据的采集，其中，所述图像数据分为正常人组和病人组两组，通过3T的磁共振设备采集，再采用弥散加权序列得到弥散图像数据，所述方法还包括如下步骤：步骤1，将所述弥散图像数据进行预处理，具体如下：步骤1.1，用基于Linux的FSL软件对所述弥散图像数据进行去涡流处理和头动矫正；步骤1.2，对经所述步骤1.1处理得出的图像进行去脑壳处理；步骤1.3，对经所述步骤1.2处理得出的图像进行弥散张量拟合得出各向异性图像、平均扩散率图像、径向弥散率图像；步骤1.4，将所述各向异性图像通过非线性配准方法配准到标准空间；步骤1.5，将所有配准到标准空间的各向异性图像进行平均得到平均各向异性图像；步骤1.6，将所述平均各向异性图像进行骨架化，得到白质骨架；步骤1.7，将步骤1.3得出的每个个体的各向异性图像、平均扩散率图像、径向弥散率图像分别投射到所得所述白质骨架上，得到每个个体的各向异性骨架图像、平均扩散率骨架图像、径向弥散率骨架图像；需要说明的是，所述病人组中的图像数据采集自神经退化类的病人。需要进一步说明的是，所述正常人组和病人组的量的最优比例为1:1。步骤2，提取步骤1得出的所述各个个体的各向异性骨架图像、平均扩散率骨架图像、径向弥散率骨架图像的特征：步骤2.1，将所述各个个体的各向异性骨架图像、平均扩散率骨架图像、径向弥散率骨架图像分别进行逐体素置换检验并以年龄作为协变量得到检验后的图像；步骤2.2，将所述检验后的图像设置阈值，得到病人组和正常人组在白质骨架上各向异性值、平均扩散率值、径向弥散率值存在显著差异的区域；步骤3，对步骤2所得的所述存在显著差异的区域进行多变量分析：步骤3.1，分别对所述存在显著差异的区域中的各向异性值、平均扩散率值、径向弥散率值进行平均，得到所述存在显著差异的区域中的平均各向异性值、平均平均扩散率值、平均径向弥散率值；步骤3.2，基于MATLAB软件将步骤3.1中所述平均的各向异性值、平均扩散率值、径向弥散率值作为特征输入到线性支持向量机中，通过留一法对线性支持向量机进行训练，最终得到特征所在区域，从而得到和病变有关的区域。需要说明的是，所述步骤3.2中，通过留一法对线性支持向量机进行训练的具体实施步骤如下：步骤3.2.1，用n表示数据样本中的个体总数，而各个体都有m个特征量，而各个体的类属性都是已知的，即病人或正常人；将得到的数据样本分为两组，一组是测试集，包含一个个体，一组是训练集，包括除测试集中所含个体外的所有人，共n-1个个体；步骤3.2.2，用所述训练集训练所述线性支持向量机，得出所述训练后支持向量机：根据以下公式计算得到权值向量w，w为一个m维列向量，其中的每个元素对应一个特征量；yi(wTxi+b)-1+ξi≥0s.t.ξi≥0；其中，γ是惩罚参数，用来实现算法复杂度和错分样本数的折中；ξi测量错分程度；yi每个人的类属性；xi为每个个体的特征向量；b为常数；需要说明的是，b是通过拉格朗日乘子法将有约束的最小化问题转化为无约束的最小化化问题，求解该无约束优化得到，而ξi大于1时，则第i个人被错分，通过拉格朗日乘子法将有约束的最小化问题转化为无约束的最小化化问题，求解该无约束优化得到。步骤3.2.3，用已知类属性的所述测试集来评估所述训练后的支持向量机的性能：用所述训练后支持向量机来判断所述测试集的类属性，所述训练后支持向量机会给出属性标签1或-1，其中1为病人，-1为正常人，通过所述训练后支持向量机得出的判断结果与所述测试集的实际类属性相比较，两者若一致，则所述训练后支持向量机分类正确，否则，则分类错误；步骤3.2.4，重新将n个个体分为测试集和训练集，所述测试集包含一个个体，且该个体与上一轮的测试集中的个体不相同，余下的所有个体作为训练集，然后按照步骤3.2.2的方法训练所述线性支持向量机，得出训练后支持向量机，然后再按照所述步骤3.2.3的方法评估得出的训练后支持向量机的性能；重复步骤3.2.4n-1次后停止；步骤3.2.5，将每个特征的n次权值求平均权值，并依据平均权值将特征由大到小进行排序，去除排序最低的特征量；步骤3.2.6，重复步骤3.2.1至步骤3.2.4，接着执行步骤3.2.7；步骤3.2.7，依据步骤3.2.6中重复步骤3.2.1至步骤3.2.4后所得出的该轮n次测试中的分类正确率和上一轮n次测试的分类正确率比较的结果判断是本文档来自技高网...

【技术保护点】
一种基于多变量的DTI图像分析方法，包括图像数据的采集，其中，所述图像数据分为正常人组和病人组两组，通过3T的磁共振设备采集，再采用弥散加权序列得到弥散图像数据，其特征在于，所述方法还包括如下步骤：步骤1，将所述弥散图像数据进行预处理：步骤1.1，用基于Linux的FSL软件对所述弥散图像数据进行去涡流处理和头动矫正；步骤1.2，对经所述步骤1.1处理得出的图像进行去脑壳处理；步骤1.3，对经所述步骤1.2处理得出的图像进行弥散张量拟合，得出各向异性图像、平均扩散率图像、径向弥散率图像；步骤1.4，将所述各向异性图像通过非线性配准方法配准到标准空间；步骤1.5，将所有配准到标准空间的所述各向异性图像进行平均得到平均各向异性图像；步骤1.6，将所述平均各向异性图像进行骨架化，得到白质骨架；步骤1.7，将步骤1.3中得出的所述各向异性图像、平均扩散率图像、图像按个体分别投射到所述白质骨架上，得到每个个体的各向异性骨架图像、平均扩散率骨架图像和径向弥散率骨架图像；步骤2，提取步骤1得出的每个个体的所述各向异性骨架图像、平均扩散率骨架图像和径向弥散率骨架图像的特征：步骤2.1，将各个体的各向异性骨架图像、平均扩散率骨架图像、径向弥散率骨架图像分别进行逐体素置换检验并以年龄因素作为协变量得到检验后的图像；步骤2.2，将所述检验后的图像设置阈值，得到所述病人组和所述正常人组在白质骨架上各向异性值、平均扩散率值、径向弥散率值存在显著差异的区域；步骤3，对步骤2所得的存在显著差异的区域进行多变量分析：步骤3.1，分别对所述存在显著差异的区域中的各向异性值、平均扩散率值、径向弥散率值进行平均，得到所述存在显著差异的区域中的平均各向异性值、平均平均扩散率值、平均径向弥散率值；步骤3.2，基于MATLAB软件将步骤3.1中所述的平均各向异性值、平均平均扩散率值、平均径向弥散率值作为特征输入到线性支持向量机中，通过留一法对线性支持向量机进行训练，最终得出特征所在的区域，从而得到与病变有关的区域。...

【技术特征摘要】
1.一种基于多变量的DTI图像分析方法，包括图像数据的采集，其中，所述图像数据分为正常人组和病人组两组，通过3T的磁共振设备采集，再采用弥散加权序列得到弥散图像数据，其特征在于，所述方法还包括如下步骤：步骤1，将所述弥散图像数据进行预处理：步骤1.1，用基于Linux的FSL软件对所述弥散图像数据进行去涡流处理和头动矫正；步骤1.2，对经所述步骤1.1处理得出的图像进行去脑壳处理；步骤1.3，对经所述步骤1.2处理得出的图像进行弥散张量拟合，得出各向异性图像、平均扩散率图像、径向弥散率图像；步骤1.4，将所述各向异性图像通过非线性配准方法配准到标准空间；步骤1.5，将所有配准到标准空间的所述各向异性图像进行平均得到平均各向异性图像；步骤1.6，将所述平均各向异性图像进行骨架化，得到白质骨架；步骤1.7，将步骤1.3中得出的所述各向异性图像、平均扩散率图像、图像按个体分别投射到所述白质骨架上，得到每个个体的各向异性骨架图像、平均扩散率骨架图像和径向弥散率骨架图像；步骤2，提取步骤1得出的每个个体的所述各向异性骨架图像、平均扩散率骨架图像和径向弥散率骨架图像的特征：步骤2.1，将各个体的各向异性骨架图像、平均扩散率骨架图像、径向弥散率骨架图像分别进行逐体素置换检验并以年龄因素作为协变量得到检验后的图像；步骤2.2，将所述检验后的图像设置阈值，得到所述病人组和所述正常人组在白质骨架上各向异性值、平均扩散率值、径向弥散率值存在显著差异的区域；步骤3，对步骤2所得的存在显著差异的区域进行多变量分析：步骤3.1，分别对所述存在显著差异的区域中的各向异性值、平均扩散率值、径向弥散率值进行平均，得到所述存在显著差异的区域中的平均各向异性值、平均平均扩散率值、平均径向弥散率值；步骤3.2，基于MATLAB软件将步骤3.1中所述的平均各向异性值、平均平均扩散率值、平均径向弥散率值作为特征输入到线性支持向量机中，通过留一法对线性支持向量机进行训练，最终得出特征所在的区域，从而得到与病变有关的区域；通过留一法对线性支持向量机进行训练的具体实施步骤如下：步骤3.2.1，用n表示数据样本中的个体总数，而各个体都有m个特征量，而各个体的类属性都...

【专利技术属性】
技术研发人员：刘鹏，杨帆，王强，刘晓明，李军，刘嫣菲，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61