一种基于二维特征和三维特征融合的阿尔茨海默病及轻度认知功能障碍识别方法

一种基于二维特征和三维特征融合的阿尔茨海默病及轻度认知功能障碍识别方法
【技术领域】
[0001] 本发明属于计算机辅助诊断领域，更为具体地讲，设及一种基于脑部医学图像的 阿尔茨海默病及轻度认知功能障碍计算机辅助识别的方法。
【背景技术】
[0002] 阿尔茨海默病（Alzheimer'SDisease,AD)是一种致死性神经退行性疾病，严重 损害了老年人的身体健康，我国八十岁W上老人发病率在30%，轻度认知功能障碍（Mild CognitiveImpairment,MCI)则被认为是介于正常老化与AD之间的一种临床状态。而随着 影像医学和计算机领域相关技术的发展，如今诊断疾病很大程度上依赖于医学图像的识别 来辅助医生进行诊断，所W通过实现基于医学图像的计算机辅助识别分类来辅助医生诊断 AD和MCI具有重大的作用和意义。其中，图像识别的前期工作和关键环节就是图像特征提 取与选择。
[0003] 传统的特征提取则主要分为底层特征提取和语义特征提取。底层特征主要有灰 度，纹理和形状=大类特征。语义特征则属于高层特征，该类特征需要研究者基于对图像内 容的先验知识W及实际应用需求从语义的角度对其信息进行描述刻画。而与普通图像相 比，医学图像本身具有高精度，信息量巨大，且与临床密切相关等专业性极强的特点。所W 近年来，国内外学者对医学图像的特征提取展开了大量研究。
[0004] 经过对现有文献分析可知，纹理特征非常适用于病变的纵向检测，目前 被广泛应用于医学图像领域，G.Castellano等人在文献（CastellanoG,Boni化a L,LiLM,etal.Textureanalysisofmedicalimages[J].Clinicalradiolo gy，2004, 59 (12) : 1061-1069.)中阐明了医学图像中常用的纹理分析技术，可分为四种形 式；结构型，基于模型，统计分析和基于变换，该四种形式各有优缺点。但是，随着技术的发 展采用单一特征已经难W准确地表达医学图像的特征，从而不利于后期的识别分类。
[0005] 本发明方案基于上述不足，提出了一种基于二维特征和=维特征融合的阿尔茨海 默病及轻度认知功能障碍计算机辅助识别方法，不但保留了纹理特征用于疾病识别的优 势，而且利用S维特征弥补了二维特征的不足，能有效提高疾病的识别率。

【发明内容】

[0006] 本发明的目的在于克服现有辅助识别阿尔茨海默病和轻度认知功能障碍方法的 不足，提出了一种基于二维纹理特征和=维形态学特征相结合的识别方法，突破了传统方 法中单个维度或者单一特征类型在反映医学图像特征上的局限性，有助于全面、准确地表 达医学图像特性，有效地提高疾病识别的准确性，从而更有效的辅助医生的临床诊断。
[0007] 为了实现上述目的，本发明提出了一种基于二维纹理特征和=维形态学特征相结 合的识别方法，包括：
[0008] 对多个样本个体对应的结构磁共振图像进行预处理，采用基于体素的脑形 态学方法（Voxel-basedmo;rphometry，VBM)和通过取幕李代数的微分同胚解剖配准(diffeomorphicanatomicalregistrationthroughexponentiatedliealgebra， DARTED方法，包括将原始图像配准到标准空间，去除非脑组织和预分割为脑白质、脑白质 W及脑脊液=部分等过程；
[0009] 选取每个样本的二维图像横断面序列图像，针对每幅图像，构建不同方向和尺度 的多个Gabor滤波器，然后提取系数幅度序列的均值及其标准方差作为特征值，W及构建 多个灰度共生矩阵，提取能量，对比度，相关度，W及滴共4维特征量；
[0010] 利用二维图像序列及其二维序列的位置信息，按其空间位置关系进行曲面重建， 使用双样本T检验对其进行统计分析，比较样本之间的脑区体积和密度的显著差异区；
[0011] 将显著差异区作为感兴趣区域（regionofinterest,R0I)，提取每个样本的所有 感兴趣区域的体积参数作为=维形态学特征；
[0012] W上提取的特征信息均采用维特征向量表示，然后利用主成分分析法对上述提取 的所有特征向量进行特征级融合，作为每个样本的多维特征向量；
[0013] 将样本分为训练集和测试集，选择径向基核函数的支持向量机（SuppcxrtVector Machine,SVM)，用训练集的特征向量对SVM分类器进行训练，训练完成后，将测试集的特征 向量输入SVM分类器，最后评价测试结果。
[0014] 本发明的特点在于基于=维形态学特征和二维纹理特征融合的阿尔茨海默病及 轻度认知功能障碍识别方法，一方面，灰度共生矩阵思想发展成熟且易于实现，在图像的局 部模式上有较好的鉴别能力，而G油or小波能较好地描述人类的视觉神经元的感受野，具 有时域和频域的联合最佳分辨率，另一方面，提取R0I的形态学特征有助于表现医学图像 特有的生理结构信息，将上述特征相融合可W得到一种具有更强鉴别力的特征向量，能够 有效地提高阿尔茨海默病及轻度认知功能障碍的识别率。
【附图说明】
[0015] 图1是本发明提出的基于=维形态学特征和二维纹理特征的阿尔茨海默病及轻 度认知功能障碍识别方法的总体流程图；
[001引图2是基于VBM-DARTCL的预处理流程图；
[0017] 图3是SVM的示意图。
【具体实施方式】
[0018] 下面结合附图对本发明的【具体实施方式】进行描述，W便本领域的技术人员更好地 理解本发明。需要注意的是，W下描述仅仅用于解释本发明，并不用于限定本发明。
[0019] 本发明提出的基于=维形态学特征和二维纹理特征的阿尔茨海默病及轻度认知 功能障碍识别方法的总体流程图如图1所示，具体如下：
[0020] 步骤（1)是对结构磁共振图像进行预处理，包括患者图像和正常人图像。采用了 基于体素的脑形态学方法和通过取幕李代数的微分同胚解剖配准方法。VBM方法是一种W 体素为单位的形态测量学方法，通过计算单位体积内包含某种体素的密度变化来显示脑组 织的形态学改变，DARTCL算法则是一种VBM分析的算法，通过对流场取幕的方法获得相应 的变形场，从而保证映射的结果是微分同胚的。
[0021] 基于VBM-DARTCL的预处理流程如图2所示；
[0022] 步骤（1. 1)首先是在原始空间上根据大脑的先验概率分布模板对图像进行预分 害d，包括去除非脑组织即顾骨剥离，W及分割出灰质，白质和脑脊液，该一过程保证了后期 统计分析的数据不会受到白质和灰质的相互影响。
[002引步骤（1. 2)利用DARTCL方法将正常人的灰质图像迭代配准到对应的平均图像上， 从而得到多个模板图像，选取最为精确的模板图像作为标准模板，再通过DARTCL配准方法 把所有灰质图像配准到该标准模板上，得到不同图像对应的变形场。
[0024] 步骤（1.扣基于步骤（1.。得到的变形场将所有灰质图像配准到丽I空间进行标 准化，然后对图像进行调制从而保证体素的总数量不变，该过程中包括仿射变换和平滑处 理，平滑处理一般采用高斯平滑核。
[0025] 步骤（2)是对图像进行特征提取，首先选取每个样本的二维图像横断面序列图 像，然后分别对每幅图像提取二维纹理特征和=维形态学特征，具体如下：
[0026] 步骤化1)对图像进行二维纹理特征提取，具体如下：
[0027] 步骤化1. 1)采用灰度共生矩阵统计分析图像中灰度的空间相关特性，提取其二 次统计量作为纹理特征。灰度共生矩阵的定义为从灰度级i的点离开某个固定位置关系d =值X，Dy)达到灰度为j的概率，用Pd(i, _]'）（i,j= 0, 1，2,…，kl)表示，其中L是图像的 灰度级，i和j表示像素的灰度，d表示两个像素间的空间位置关系，包括了距离和方向，一 般方向取0°，45°，90° ,135°四个方向。通过图像构建灰度共生矩阵之后将各个元素 Pd(i，j)除W各元素之和S得到归一化值氏(i，j)，由此得到归一化共生矩阵。考虑到便于计 算W及能够提供较高的分类精度，一般提取如下四个特征值作为纹理特征：
[002引 1)能量；
[0029]
[0030] 能量反映了图像灰度分布均匀程度和纹理粗细程度。
[003。。对比度；
[0032]
[0033] 对比度反映了图像的清晰度和纹理沟纹深浅的程度。
[0034] 如相关；

[0041] 相关用于衡量灰度共生矩阵的元素在行或列方向上的相似程度。
[0042] 4)滴；
[0043]
[0044] 滴反映图像中纹理的复杂程度或非均匀度。
[0045] 步骤化1.2)是基于G油or多通道滤波技术，提取多尺度多方向的纹理特征。 G油or函数是目前唯一能够达到时频测不准关系下届的函数，能最好地兼顾信号在时域和 频域中的分辨能力。二维G油or函数及其傅里叶变换如公式所示：
[0046]
[0047] 但是 由于G油or函数本身不具有小波函数的正交特性，所W需要将g(x，y)作为母 小波对其进行伸缩和变换处理，得到Gabor小波如公式所示：
[0048] g血(X，y) =a_mg(X，y)
[0049] 其中
[0050] X=a_m(xcos日+ysin白）
[0051] y=a-m(ycos日一xsin白）
[005引a^m即是尺度因子，m和n分别是指定小波的尺度和方向，即是说，"1 = 0，1，2-" ，S-1，n= 0, 1，2…，t-1，S和t分别是总的尺度数和总的方向数。即是说，针对WXH的医 学图像H(x，y)的离散G油or小波变换是：
[0053]
[0054] 其中g*_"(i，如是gmn(i.如的共辆复数，i和j即是窗口大小（即滤波器模板大 小）的变量，此时设计的滤波器需要确保Gabor小波滤波器族的响应在频谱半高宽能相互 接触但不重叠。然后由于系数幅度值能够表征不同方向和尺度上的图像能量，所W，计算出 不同方向和尺度上变换后所得的系数幅度序列为：
[00巧]
[0056] 因此，纹理特征分别是系数幅度序列的均值y皿及其标准方差0皿：
[0059] 步骤化2)对图像进行S维形态学特征提取，具体如下：
[0060] 步骤化2. 1)是选取感兴趣区域。利用二维图像序列及其二维序列的位置信息， 按其空间位置关系进行曲面重建，然后将不同组图像中的每个体素作为参数，进行双样本T检验，找出具有显著差异的区域，则是具有统计学意义的簇，将该些簇作为感兴趣区域。其 中，为了使簇具有显著差异，在统计过程中可设置P《0. 05,并且进行抑R校正。
[0061] 步骤化2. 2)是提取每个样本的感兴趣区域的体积特征。将在步骤（2. 2. 1)中得 到的感兴趣区域制作成多个mask,然后用每个mask点乘每个样本的灰质图像，从而计算出 每个样本的感兴趣区域的体积。所W，每个样本的感兴趣区域的体积就是一组特征向量。
[0062] 步骤（3)是对步骤（2)提取所得的特征进行选择融合，W便消除特征间的冗余信 息，且在保留了所需识别信息的基础上降低特征维数，减少计算代价。主成分分析法是一种 把多个实测变量转换为少数几个不相关的综合指标（即主成分）的多元统计分析方法，其 基本思想是基于实测变量之间必定存在一定相关关系，因此，可W用综合各变量信息的互 不重叠的综合指标来代表该些实测变量。
[0063] 每幅图像通过步骤（2)提取了n维特征，一共有k幅图像，那么该个原始特征集就 是：
[0064]
[0065] 对所有特征值采用标准差变换和极差转换进行标准化处理之后，使用XI，X2,…Xk 线性表示综合变量Yi，y2，一yk;
[0066]
[0067]其中，需满足如下条件；（l)?fi+z4H+i4 = 1,/ = ;似y;和yj互不相关； (3)Yi是Xi，X2，"'Xk所有线性组合中的方差最大值，7 2贝是与前者不相关的Xi，X2，…Xk所有线性组合中的方差最大值，W此类推，yk是与之前所有综合变量不相关的X1，X,，…Xk 所有线性组合中的最大方差值。
[0068] 综上所述，得到k个新的综合变量（即主成分），然后根据矩阵代数，把特性向量 Ui，u2，所对应的非零特征根设为…^k，则可W定义前P个主成分yi，y2，… yk的累计贡献率为：
[0069]
[0070] 当P个主成分的累计贡献率超过85%时，即选取该P个主成分作为特征融合的结 果。
[0071] 步骤（4)基于SVM方法对样本图像构建分类器，实现对阿尔茨海默病，轻度认知 功能障碍和正常人之间的自动识别。SVM的基本思想是将输入数据通过内积函数费线性映 射到一个高维空间，然后在该个空间建立一个最大间隔超平面同时求解出最优的线性分类 面，其示意图如图3所示，每个中间节点对应一个支持向量，输出是中间节点的线性组合。
[0072] 构建SVM分类器时，将步骤（3)融合所得的特征向量作为输入，采用径向基核函数 W及leave-one-out交叉验证法设置参数。由于是=分类问题，所W采用一对一的投票策 略，即总共构建3个二分类器，训练集分别对应于（AD，正常人），（MCI，正常人），（AD，MCI)， 测试时将测试样本依次送入该=个分类器，最后得到一组结果，取分类结果投票数最多的 类别作为最终结果，若出现同样多的投票，一般取第一个最大值对应的类别。最后W平均预 测正确率作为评价结果标准。
[0073] W上对本发明所提出的计算机辅助识别阿尔茨海默病W及轻度认知功能障碍方 法进行了详细的介绍，但是【具体实施方式】的描述仅用于解释本发明的方法及其核屯、思想， W便于本技术领的技术人员理解本发明，但应该清楚，本发明不限于【具体实施方式】的范围， 对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明 的精神和范围内，该些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
【主权项】
1. 一种基于二维纹理特征和三维形态学特征结合的阿尔茨海默病和轻度认知功能障 碍的识别方法，包括以下步骤： 步骤一：对多个样本个体对应的结构磁共振图像进行预处理，采用基于体素的脑形态 学方法和通过取幂李代数的微分同胚解剖配准方法，包括配准到标准空间，去除非脑组织 和预分割为脑白质、脑白质以及脑脊液三部分等过程； 步骤二：选取每个样本的二维图像横断面序列图像，针对每幅图像，构建不同方向和尺 度的多个Gabor滤波器，然后提取系数幅度序列的均值及其标准方差作为特征值，以及构 建多个灰度共生矩阵，提取能量，对比度，相关度，以及熵作为特征值； 利用二维图像序列及其二维序列的位置信息，按其空间位置关系进行曲面重建，使用 双样本T检验对其进行统计分析，将样本之间的脑区体积和密度的显著差异区作为感兴趣 区域，提取每个样本的所有感兴趣区域的体积参数作为三维形态学特征； 步骤三：以上提取的特征信息均采用特征向量表示，然后利用主成分分析法对上述提 取的所有特征向量进行特征级融合，作为每个样本的多维特征向量； 步骤四：将样本分为训练集和测试集，选择径向基核函数的支持向量机（SVM)，用训练 集的特征向量对SVM分类器进行训练，训练完成后，将测试集的特征向量输入SVM分类器， 最后评价测试结果。2. 根据权利要求1所述的基于二维纹理特征和三维形态学特征结合的阿尔茨海默病 和轻度认知功能障碍识别方法，其特征在于，所述步骤二和三中将二维纹理特征和三维形 态学特征相结合，具体包括如下步骤： 步骤（1. 1. 1)选取每个样本的二维图像横断面序列图像，为每幅图像构建多个灰度共 生矩阵，提取如下四个特征值作为二维纹理特征： 能量：步骤（1. 1.2)为每幅图像构建不同方向和尺度的多个Gabor滤波器，计算不同尺度和 方向上的系数幅度序列：然后计算如下两个特征值作为二维纹理特征： 系数幅度序列的均值： mmn= E (m, n) /WH 系数幅度序列的标准差：步骤（1.2. 1)利用二维图像序列及其二维序列的位置信息，按其空间位置关系进行曲 面重建，然后将不同组的图像中的每个体素作为参数，进行双样本T检验，找出具有显著差 异的区域作为感兴趣区域，计算出每个样本的感兴趣区域的体积作为三维形态学特征； 步骤（2)对上述得到的特征进行标准化处理，然后采用主成分分析法对其进行融合， 选取其中m个最具判别性的综合指标，即主成分作为每个样本的多维特征向量。3.如权利要求1所述的一种基于二维纹理特征和三维形态学特征结合的阿尔茨海默 病和轻度认知功能障碍识别方法，其特征在于，所述步骤四中的基于SVM的阿尔茨海默病 及轻度认知功能障碍识别的具体步骤如下： 构建三个支持向量机分类器，采用径向基核函数，将经过特征融合得到的特征向量作 为输入，采用一对一投票策略，对阿尔茨海默病，轻度认知功能障碍和正常人进行分类识 另IJ，最后以平均预测正确率作为评价结果标准。
【专利摘要】本发明提出了一种基于三维形态学特征和二维纹理特征的阿尔茨海默病及轻度认知功能障碍计算机辅助识别方法，具体包括：1)对医学图像进行预处理，包括预分割和配准等过程；2)对医学图像进行二维纹理特征提取，包括灰度共生矩阵的二次统计量和Gabor小波变换的多尺度多方向特征值；3)对医学图像进行三维形态学特征提取，即提取感兴趣区域的体积特征；4)对三维形态学特征和二维纹理特征进行特征融合；5)构建支持向量机分类器实现对阿尔茨海默病和轻度认知功能障碍的识别。本发明将三维形态学特征和二维纹理特征结合，能够更加全面准确地表达医学图像的内容，所以该方法能够提高对阿尔茨海默病和轻度认知功能障碍的识别，从而提供更有效的临床辅助诊断。
【IPC分类】G06K9/62, G06K9/46
【公开号】CN104881680
【申请号】CN201510268415
【发明人】蓝天, 王伟, 丁熠, 秦臻, 张聪, 黄若菡, 陈浩, 肖哲, 徐路路, 陈圆
【申请人】电子科技大学
【公开日】2015年9月2日
【申请日】2015年5月25日