一种基于二维特征和三维特征融合的阿尔茨海默病及轻度认知功能障碍识别方法
【技术领域】
[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] 但是