一种基于复杂网络的脑功能磁共振图像分类方法

文档序号:6439798阅读:393来源:国知局
专利名称:一种基于复杂网络的脑功能磁共振图像分类方法
技术领域
本发明属于图像处理技术领域,具体涉及一种基于复杂网络的脑功能磁共振图像分类方法。
背景技术
功能磁共振成像(functional Magnetic Resonance Imaging, fMRI)以其高时空分辨率,非侵入式等特点在神经疾病诊断治疗方面得到了广泛应用。fMRI —般指基于血氧水平依赖(blood oxygen level-dependent, BOLD)的磁共振成像,它通过测量由神经活动引起的脑血流和脑血氧等成分变化而造成的磁共振信号变化来反应脑活动。脑是一个复杂的系统,在受到刺激条件或经历病变时脑的磁共振图像会发生相应的变化。利用图像分类方法,计算脑功能磁共振图像具有某种属性的可能性大小,或者自动判别图像的类别属性, 是计算机辅助分析的一个重要应用。传统的功能磁共振图像分类方法主要有感兴趣区域(ROI)方式和体素(voxel)方式两种分类方法。感兴趣区域方式的分类方法依据目标结构的先验知识,将样本和目标分割成多个目标区域,并据此对目标进行分类;体素方式的分类方法采用复杂的非线性配准, 以最大限度地实现个体间的精确对应,然后以图像的每一个空间单位(体素)作为分类依据。这两种方法都假设目标与样本的内部组织结构是一一对应的。前者认为先验的图像区域存在于每一个目标图像当中,并且能够准确分割;后者假定非线性配准后的体素是一一对应的。然而,这样的假设在很多情况下并不合理。人在不同状态下的脑功能磁共振图像会受到多方面因素的干扰,传统的分类方法都不是根据脑的固有属性对脑功能磁共振图像进行分类的,因此都会导致分类性能的下降。

发明内容
(一 )要解决的技术问题为了克服已有技术的不足,本发明所要解决的技术问题是设计一种分类准确率高、泛化性能强的脑功能磁共振图像分类方法。( 二 )技术方案为实现上述目的,本发明提出一种基于复杂网络的脑功能磁共振图像分类方法, 包括以下步骤步骤M 对训练样本图像和测试样本图像进行预处理,然后进行脑区分割,并提取各个脑区的平均时间序列;步骤Sb 计算各个平均时间序列之间的偏相关系数,得到偏相关系数矩阵;步骤& 将所述偏相关系数矩阵二值化,得到复杂网络模型;步骤Sd 计算该复杂网络模型的特征路径长度、成本和集群度作为功能磁共振图像的特征;步骤Se 利用训练样本图像的网络参数作为该功能磁共振图像的特征中的训练样本图像的特征,来训练一自适应提高(adaboost)分类器;步骤Sf 利用训练好的该自适应提高(adaboost)分类器对测试样本图像进行分类。(三)有益效果本发明针对脑功能磁共振图像分类问题,通过构建脑网络模型、计算网络特征参数、训练自适应提高(adaboost)分类器等方法有效提高了图像分类的准确性和稳定性。本发明能够利用脑功能磁共振图像中尽可能多的信息,脑网络参数能够从本质上反应脑的活动,弥补了传统分类方法不能体现脑活动固有属性的不足,能够精确的对脑功能磁共振图像进行分类。


图1是本发明提供的基于复杂网络的脑功能磁共振图像分类的方法流程图;图2是依照本发明实施例使用本发明所述分类方法(方法A)对比现有的基于局部特征的分类方法(方法B),分类受试者操作特性(ROC)的对比曲线。
具体实施例方式为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。基于复杂网络的脑功能磁共振图像分类是一种全新的脑功能磁共振图像分类方法。该方法首先建立复杂脑网络模型,计算脑网络的特征路径长度和集群度,用以表征不同的图像模式;然后利用该特征路径长度和集群度来训练一个自适应提高(adaboost)分类器;最后利用训练好的该自适应提高(adaboost)分类器对测试样本图像进行分类。参照图1,根据本发明所述的一种人脑功能磁共振成像图像分类方法,能够依据训练样本图像来确定测试样本图像的类别,具体实施步骤如下步骤Μ,对训练样本图像和测试样本图像进行预处理,然后进行脑区分割,并提取各个脑区的平均时间序列;1.脑功能磁共振图像的预处理由于磁共振扫描过程中各种各样的噪声的影响,被试个体自身存在尺度和位置上的差异,非常有必要在分析数据之前对数据做一定的预处理。在整个的实验的数据获取中, 主要的噪声信息来源有(1)物理头动;(2)图像内层间扫描时间差别;(3)外在磁场的不均勻性等。脑功能磁共振图像预处理的常见步骤有切片扫描时间对齐,图像序列对齐,联合配准,标准化(或称均一化),空间平滑滤波和时间平滑滤波等。2.脑功能磁共振图像的分割采用国际通用的结构标记模板(AAL),将全脑分为90个脑区。结构标记模板是脑功能磁共振图像研究领域使用最为广泛的脑结构模板。3.提取各脑区的平均时间序列依据预处理后的脑功能磁共振图像的数据,提取包含于相应脑区内部的各个体素在不同时间点上激活值的时间序列Y(矩阵维数DXN),其中D为包含于球体内部的体素数目,N为时间点数。所述激活值是指各个体素在不同时间点上的血氧水平依赖(BOLD)强度。
步骤Sb 计算平均各个时间序列之间的偏相关系数。该步骤Sb具体包括如下步骤1.计算平均时间序列间的协方差系数依据步骤提取的各个脑区的时间序列,计算各个平均时间序列之间的协方差矩阵S,S的每个元素Sy为第i个和第j个时间序列之间的协方差系数,
权利要求
1.一种基于复杂网络的脑功能磁共振图像分类方法,其特征在于,包括以下步骤 步骤Μ 对训练样本图像和测试样本图像进行预处理,然后进行脑区分割,并提取各个脑区的平均时间序列;步骤Sb 计算各个平均时间序列之间的偏相关系数,得到偏相关系数矩阵; 步骤& 将所述偏相关系数矩阵二值化,得到复杂网络模型; 步骤Sd 计算该复杂网络模型的特征路径长度、成本和集群度作为功能磁共振图像的特征;步骤利用训练样本图像的网络参数作为该功能磁共振图像的特征中的训练样本图像的特征,来训练一自适应提高(adaboost)分类器;步骤Sf 利用训练好的该自适应提高(adaboost)分类器对测试样本图像进行分类。
2.如权利要求1所述的基于复杂网络的脑功能磁共振图像分类方法,其特征在于,所述对训练样本图像和测试样本图像进行预处理,是在保留脑功能图像细节的同时,使用脑功能图像与标准模板进行仿射配准变换方式的预处理,并提高脑功能图像的信噪比。
3.如权利要求1所述的基于复杂网络的脑功能磁共振图像分类方法,其特征在于,所述提取脑区平均时间序列的步骤为首先按照标准脑结构模板将全脑分为90个脑区,分别提取每个脑区内部各个体素在不同时间点上的激活值,再将各个体素的激活值进行平均,得到脑区平均时间序列。
4.如权利要求1所述的基于复杂网络的脑功能磁共振图像分类方法,其特征在于,计算各个平均时间序列之间的偏相关系数的方法为首先计算各个平均时间序列之间的协方差矩阵S,该协方差矩阵维度为90X90,S的每个元素Sy为第i个和第j个时间序列之间的协方差系数,
5.如权利要求1所述的基于复杂网络的脑功能磁共振图像分类方法,其特征在于,所述将所述偏相关系数矩阵二值化得到复杂网络模型的步骤,包括选取阈值将经过Fisher变换的偏相关系数矩阵F 二值化,该变换后的偏相关系数矩阵维度为90X90,二值化后1表示两个脑区之间有连接,即网络中两个节点之间的边存在,0则表示两个脑区之间没有连接,即网络中的两个节点之间没有边;阈值选取的方法为使选用此阈值进行了二值化后的网络中实际存在的边的数量是网络中可能存在的边的数量的十分之一。
6.如权利要求1所述的基于复杂网络的脑功能磁共振图像分类方法,其中所述计算该复杂网络模型的特征路径长度的步骤为用特征路径长度矩阵描述网络中任意两个节点i,j的特征路径长度Iij,网络平均特征路径长度L描述了网络中任意两个节点的特征路径长度的平均值,即L=~1~ y /,N(N-I)1^mj v其中,N为网络中节点的个数,即分割的脑区数90 ;Iij为节点i,j之间的特征路径长度,V为网络中所有节点的集合。
7.如权利要求1所述的基于复杂网络的脑功能磁共振图像分类方法,其中所述计算该复杂网络模型的成本的步骤为用网络中实际存在的所有边的数量比上网络中最多可能存在的边的数量,即κ _ 2 _ 1 y κ “NH聊-\、乙 2其中,N为网络中节点的个数,Ki为网络中连接到节点i的边的数量。
8.如权利要求1所述的基于复杂网络的脑功能磁共振图像分类方法,其中所述计算该复杂网络模型的集群度的步骤为某一节点i的集群度Ci的值等于它的相邻节点之间存在的边的数目与它们之间所有可能的边数的比值,即C - e' - 2e'! MVi1) ^ikl-I) 2其中,e,表示节点i的邻点之间存在的边数,Iii表示节点i的邻点的数目,^1 表示节点i的邻点之间可能存在的边数。
全文摘要
本发明是一种基于复杂网络的脑功能磁共振图像分类方法,该方法包括以下步骤对训练样本图像和测试样本图像进行预处理,并进行区域分割,提取每个区域的平均时间序列;计算该平均时间序列间的偏相关系数,将所述偏相关系数矩阵二值化,得到复杂网络模型,计算该复杂网络模型的特征路径长度、成本和集群度,分别得到训练样本图像和测试样本图像的网络特征;利用训练得到的该自适应提高(adaboost)分类器;利用训练得到的该自适应提高(adaboost)分类器对测试样本图像进行分类。本发明利用了脑功能磁共振图像中尽可能多的信息,能够精确的对脑功能磁共振图像进行分类。
文档编号G06K9/66GK102509123SQ201110392269
公开日2012年6月20日 申请日期2011年12月1日 优先权日2011年12月1日
发明者刘振宇, 田捷, 白丽君 申请人:中国科学院自动化研究所
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1