基于特征学习的ct图像肝脏分割方法及系统的制作方法
【专利摘要】本发明公开一种基于特征学习的CT图像肝脏分割方法及系统，能有效提高CT图像中肝脏的分割精度。所述方法包括：S101、读取训练图像集和待分割图像，其中，所述训练图像集中的训练图像和待分割图像为腹部的CT图像；S102、提取所述训练图像和待分割图像的Haar特征、局部二进制模式特征、方向梯度直方图特征和共生矩阵特征；S103、利用主成分分析方法将提取的所有特征进行特征融合，获取更有效的特征；S104、利用分类器对所述待分割图像每个像素的特征进行分类，得到肝脏概率图；S105、结合所述肝脏概率图与待分割图像，修改随机游走分割算法的图模型权值，实现肝脏的分割。
【专利说明】
基于特征学习的CT图像肝脏分割方法及系统
技术领域
[0001] 本发明涉及机器学习技术领域，具体涉及一种基于特征学习的CT图像肝脏分割方 法及系统。
【背景技术】
[0002] 医学图像分割辅助医生识别病人的内部组织器官及病灶区域，在计算机辅助治疗 及手术规划中发挥至关重要的作用。所以，肝脏的自动分割是医生诊治如肝硬化、肝脏肿 瘤、肝移植等肝脏疾病的基础。在腹部CT图像中，肝脏与邻近器官的灰度值差异较小，肝脏 本身灰度不均匀且其形状各异，自动、精确的分割出肝脏难度较大。所以，临床医生迫切需 要一种简单、快速、准确的肝脏分割方法。
[0003] 现有的随机游走分割方法具有快速简单等优点，但它对CT图像中对比度低的区域 分割效果较差，特别是肝脏与大血管、胃等邻近器官的连接处，单纯地依赖灰度值难以有效 地实现肝脏的分割。
【发明内容】
[0004] 有鉴于此，本发明实施例提供一种基于特征学习的CT图像肝脏分割方法及系统， 能够有效提高CT图像中肝脏的分割精度。
[0005] -方面，本发明实施例提出一种基于特征学习的CT图像肝脏分割方法，包括：
[0006] S101、读取训练图像集和待分割图像，其中，所述训练图像集中的训练图像和待分 割图像为腹部的CT图像；
[0007] S102、提取所述训练图像和待分割图像的Haar特征、局部二进制模式特征、方向梯 度直方图特征和共生矩阵特征；
[0008] S103、利用主成分分析方法将提取的所有特征进行特征融合，获取更有效的特征；
[0009] S104、利用分类器对所述待分割图像每个像素的特征进行分类，得到肝脏概率图；
[0010] S105、结合所述肝脏概率图与待分割图像，修改随机游走分割算法的图模型权值， 实现肝脏的分割。
[0011] 另一方面，本发明实施例提出一种基于特征学习的CT图像肝脏分割系统，包括：
[0012] 读取模块，用于读取训练图像集和待分害_像，其中，所述训练图像集中的训练图 像和待分割图像为腹部的CT图像；
[0013] 提取模块，用于提取所述训练图像和待分割图像的Haar特征、局部二进制模式特 征、方向梯度直方图特征和共生矩阵特征；
[0014] 融合模块，用于利用主成分分析方法将提取的所有特征进行特征融合，获取更有 效的特征；
[0015] 分类模块，用于利用分类器对所述待分割图像每个像素的特征进行分类，得到肝 脏概率图；
[0016] 分割模块，用于结合所述肝脏概率图与待分割图像，修改随机游走分割算法的图 模型权值，实现肝脏的分割。
[0017] 本发明实施例提供的基于特征学习的CT图像肝脏分割方法及系统，通过融合图像 的不同特征并分类，实现图像中像素点的学习，获得图像对应的肝脏概率图，在此基础上， 利用随机游走算法自动地分割出肝脏区域。实验结果证明，本发明能够获得较好的分割结 果，在对图像进行特征学习的基础上结合随机游走算法，相比于单独地使用随机游走算法， 更加有效地提高了 CT图像中肝脏的分割精度。
【附图说明】
[0018] 图1为本发明基于特征学习的CT图像肝脏分割方法一实施例的流程示意图；
[0019] 图2为本发明基于特征学习的CT图像肝脏分割系统一实施例的结构示意图。
【具体实施方式】
[0020] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例 中的附图，对本发明实施例中的技术方案进行清楚地描述，显然，所描述的实施例是本发明 一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有 做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0021] 如图1所示，本实施例公开一种基于特征学习的CT图像肝脏分割方法，包括：
[0022] S101、读取训练图像集和待分割图像，其中，所述训练图像集中的训练图像和待分 割图像为腹部的CT图像；
[0023] S102、提取所述训练图像和待分割图像的Haar特征、局部二进制模式特征、方向梯 度直方图特征和共生矩阵特征；
[0024] S103、利用主成分分析方法将提取的所有特征进行特征融合，获取更有效的特征；
[0025] S104、利用分类器对所述待分割图像每个像素的特征进行分类，得到肝脏概率图；
[0026] S105、结合所述肝脏概率图与待分割图像，修改随机游走分割算法的图模型权值， 实现肝脏的分割。
[0027] 本实施例提供的基于特征学习的CT图像肝脏分割方法，通过融合图像的不同特征 并分类，实现图像中像素点的学习，获得图像对应的肝脏概率图，在此基础上，利用随机游 走算法自动地分割出肝脏区域。实验结果证明，本发明能够获得较好的分割结果，在对图像 进行特征学习的基础上结合随机游走算法，相比于单独地使用随机游走算法，更加有效地 提高了 CT图像中肝脏的分割精度。
[0028] 可选地，在本发明基于特征学习的CT图像肝脏分割方法的另一实施例中，所述 S103中具体做法为：
[0029]对图像中的任一样本点，用以上四种特征组合表示其周围邻域的局部纹理信息：
[0030] f =[fHaar；fLBP；fHOG；fGLM],
[0031 ] 其中，fHaar为Haar特征，fLBP为局部二进制模式特征，fHQG为方向梯度直方图特征， fcLM为共生矩阵特征；
[0032] 对组合后的特征进行降维，以去除特征中相互关系密切的部分，获取降维后的像 素点的特征。
[0033] 组合后的特征维数较高，特征间可能存在信息冗余，造成计算内存的浪费。本发明 通过对高维特征的降维，去除特征中相互关系密切的部分，减小噪声干扰，提取出最能反映 其特征的数据。
[0034] 可选地，在本发明基于特征学习的CT图像肝脏分割方法的另一实施例中，所述降 维的方法为主成分分析方法、独立成分分析方法、线性判别分析方法、局部线性嵌入方法、 拉普拉斯特征映射方法或互信息提取方法。
[0035] 可选地，在本发明基于特征学习的CT图像肝脏分割方法的另一实施例中，所述 S104具体为：
[0036]对所提取的图像纹理特征进行分类，由于所提取的腹部CT图像纹理特征具有较强 的非线性特点，且样本量较大，本发明采用AdaBoost算法对提取样本点集的特征信息进行 分类，构造区分肝脏与非肝脏纹理信息的分类器，利用分类器计算测试图像像素点x的肝脏 概率p(x)。
[0037] 在训练过程中，利用分类回归树作为弱分类器，通过训练将T个弱分类器集成为一 个强分类器。N个训练样本的初始样本权值为：
[0038] Di(n) = 1/N,
[0039] 其中，N为样本点数量。当一个弱分类器训练完成后，根据该弱分类器对训练样本 集的分类错误率：
[0040] Y ￡>丨(")[1( n=l
[0041] 将其权重调整为：
[0043]同时更新样本权值：
[0044] D.,,(/;) - Dt(n)cxp(-c/yJi(/,)) / Z ,
[0045] 其中，D(n)为样本权值，弱分类器1!(匕'）￡{-1，1}，1为分类结果属于肝脏，-1为分 类结果不属于肝脏，表示训练图像中像素点本来属于肝脏但分类结果不属于肝脏 时结果等于l，h表示弱分类器，f' n表示第n个样本点的降维后特征向量，h(f'n)为弱分类器 对特征f'n进行分类，D t(n)为第t次迭代更新后的样本权值，.0,为第t次迭代的权重，yn表示 第n个样本点的类别; Z为归一化因子。被错误分类的样本相应的新权值较高，使得下一次训 练的弱分类器更倾向于对这些样本的分类。
[0046] 在测试过程中，对于测试样本X，其特征向量为f'x，集成后的强分类器的判决结果 为所有弱分类器判决结果的加权和：
[0047] 峨卜坌哪'；)，
[0048] 利用此加权和，计算测试像素点x属于肝脏区域的概率估计：
[0050]其中，ht(f'x)为弱分类器对特征f'x进行第t次分类，y=l表示像素点属于肝脏，p (y=l |x)为测试像素点x属于肝脏区域的概率估计，y为属于肝脏区域的指引，H为所有弱分 类器。
[0051]所使用的分类器可以是利用分类回归树作为弱分类器的AdaBoost分类器，以支持 向量机作为弱分类器的AdaBoost分类器，支持向量机分类器，决策树分类器，人工神经网络 分类器，朴素贝叶斯分类器，随机森林分类器等，且不局限于这几种分类器。对于测试图像 上的每个像素点，都可得到相应的肝脏概率估计。对于任意一幅待测试的腹部CT图像，提取 纹理特征并分类，便可得到对应的肝脏概率图，自动实现测试图像像素点的模糊分类。 [0052]可选地，在本发明基于特征学习的CT图像肝脏分割方法的另一实施例中，所述 S105具体为：
[0053] 对于待分割图像Itest，由统计可知，肝脏区域的灰度值对应腹部图像灰度值统计 的某一固定峰值附近。基于此特征，针对肝脏区域所占面积较大的第Z0层切片ItMtUo)，选 取待分割图像灰度统计图中，此峰值两边的波谷 S1、S2作为选取肝脏区域的阈值，阈值化的 二值图像包括包含肝脏区域的二值图像gl和包含背景区域的二值图像g 2:
[0054] gi(k，l，zo) = l(si<I(Itest(k，l，Z0)Xs2)，
[0055] g2(k，l，zo) = l(I(Itest(k，l，zo)Xsi)，
[0056] 其中，k，l代表像素点在切片中的位置，Itest(k，l，ZQ)表示切片zo中k，l对应位置处 的像素点，K*)表示像素点*的灰度值，对图像g 2做腐蚀运算，即获取第Z0层切片非肝脏区域 的确信像素点;对图像gi做腐蚀运算并找到最大连通区域n，再做腐蚀运算后的像素点即为 第zo层切片肝脏区域的确信像素点，
[0057]对于除zo层的其他切片，为防止当肝脏所占面积较小时，无法由最大连通域找到 肝脏区域，本发明利用相邻切片间的连接信息确定肝脏区域，首先，利用相邻切片的最大连 通区域n作为先验信息，找出包含n的最小矩形区域并向外延拓八个像素;然后，在此区域 上重复上述步骤，便可自动选取所有切片中肝脏与非肝脏区域的确信像素点作为随机游走 的种子点 seedin、seedout，
[0058] 原始的随机游走算法只根据种子点提取CT图像的先验信息，并没有有效地利用图 像特征空间的信息，本发明将随机游走推广至特征空间，结合图像空间的先验信息并用于 图像的特征空间中，基于待分割图像对应的肝脏概率图，利用改进的随机游走算法确定肝 脏的边界，考虑到在原始待分割图像中，对于肝脏边缘模糊区域，像素点的对比度较低，难 以识别肝脏区域，本发明通过利用肝脏概率图的信息，首先将其映射为带权的无向图G = (V，E)，结点集V由肝脏概率图P的像素点集构成，边集E由邻接像素关系构成，其权值由像素 点的先验概率决定：
[0059] = QKp(~P(P(i) -P(j)Y) y
[0060] 其中，i，j为待分割图像Ites冲的像素索引，ei伪两像素间的相邻关系，0为调整参 数，p(i)和P(j)分别为索引为i的像素的先验概率和索引为j的像素的先验概率，然后计算 每个未标记像素点到达肝脏区域种子点seecU的概率q:
[0061 ] V =呢min [(叫 | P(7)- P(./) I)- p e^E 1
[0062] s. t.q(seedin) = l，q( seedout )=0，
[0063] 其中，q(seedin)表示seedin的先验概率，q( seed〇ut)表示seed〇ut的先验概率。
[0064] 最后判断此概率，大于或者等于1/2时，像素点的类别label被判为1，即属于肝脏 区域;否则，属于非肝脏，最终实现肝脏的分割，即： fl. (/(/)> 1/2
[0065] label(i)= <
[0, (/(/')< 1/2
[0066] 可选地，在本发明基于特征学习的CT图像肝脏分割方法的另一实施例中，对图像 g2做腐蚀运算所使用的腐蚀参数为6,对图像 gl做的两次腐蚀运算所使用的腐蚀参数为4。 [0067]如图2所示，本实施例公开一种基于特征学习的CT图像肝脏分割系统，包括：
[0068] 读取模块1，用于读取训练图像集和待分割图像，其中，所述训练图像集中的训练 图像和待分割图像为腹部的CT图像；
[0069] 提取模块2,用于提取所述训练图像和待分割图像的Haar特征、局部二进制模式特 征、方向梯度直方图特征和共生矩阵特征；
[0070] 融合模块3,用于利用主成分分析方法将提取的所有特征进行特征融合，获取更有 效的特征；
[0071] 分类模块4,用于利用分类器对所述待分割图像每个像素的特征进行分类，得到肝 脏概率图；
[0072] 分割模块5,用于结合所述肝脏概率图与待分割图像，修改随机游走分割算法的图 模型权值，实现肝脏的分割。
[0073] 本实施例提供的基于特征学习的CT图像肝脏分割系统，通过融合图像的不同特征 并分类，实现图像中像素点的学习，获得图像对应的肝脏概率图，在此基础上，利用随机游 走算法自动地分割出肝脏区域。实验结果证明，本发明能够获得较好的分割结果，在对图像 进行特征学习的基础上结合随机游走算法，相比于单独地使用随机游走算法，更加有效地 提高了 CT图像中肝脏的分割精度。
[0074]虽然结合附图描述了本发明的实施方式，但是本领域技术人员可以在不脱离本发 明的精神和范围的情况下做出各种修改和变型，这样的修改和变型均落入由所附权利要求 所限定的范围之内。
【主权项】
1. 一种基于特征学习的CT图像肝脏分割方法，其特征在于，包括： 5101、 读取训练图像集和待分割图像，其中，所述训练图像集中的训练图像和待分割图 像为腹部的CT图像； 5102、 提取所述训练图像和待分割图像的Haar特征、局部二进制模式特征、方向梯度直 方图特征和共生矩阵特征； 5103、 利用主成分分析方法将提取的所有特征进行特征融合，获取更有效的特征； 5104、 利用分类器对所述待分割图像每个像素的特征进行分类，得到肝脏概率图； 5105、 结合所述肝脏概率图与待分割图像，修改随机游走分割算法的图模型权值，实现 肝脏的分割。2. 根据权利要求1所述的基于特征学习的CT图像肝脏分割方法，其特征在于，所述S103 中具体做法为： 对图像中的任一样本点，用以上四种特征组合表示其周围邻域的局部纹理信息： f = [fHaar；fLBP；fHOG；fGLM], 其中，fHaar为Haar特征，fLBP为局部二进制模式特征，fHQC为方向梯度直方图特征，fcLM为 共生矩阵特征； 对组合后的特征进行降维，以去除特征中相互关系密切的部分，获取降维后的像素点 的特征。3. 根据权利要求2所述的基于特征学习的CT图像肝脏分割方法，其特征在于，所述降维 的方法为主成分分析方法、独立成分分析方法、线性判别分析方法、局部线性嵌入方法、拉 普拉斯特征映射方法或互信息提取方法。4. 根据权利要求1述的基于特征学习的CT图像肝脏分割方法，其特征在于，所述S104中 的分类器为决策树分类器、人工神经网络分类器、朴素贝叶斯分类器、随机森林分类器、支 持向量机分类器、利用分类回归树作为弱分类器的AdaBoost分类器或者支持向量机作为弱 分类器的AdaBoost分类器。5. 根据权利要求1述的基于特征学习的CT图像肝脏分割方法，其特征在于，所述S104具 体为： 对于所述待分割图像上的每个像素点X，利用分类器对该像素点X的特征进行分类，都 可得到相应的肝脏概率估计P(X)，从而得到对应的肝脏概率图P，自动实现所述待分割图像 像素点的模糊分类。6. 根据权利要求1所述的基于特征学习的CT图像肝脏分割方法，其特征在于，所述S105 具体为： 对于待分割图像ltd，针对肝脏区域所占面积较大的第ZO层切片Ite3st(ZO),选取待分割 图像灰度统计图中，此峰值两边的波谷S1、s2作为选取肝脏区域的阈值，阈值化的二值图像 包括包含肝脏区域的二值图像 81和包含背景区域的二值图像g2: gl ( k , I , ZO ) 一 I ( Sl^ I ( Itest ( k , I , ZO ) ) ^ S2 ) 7 g2 ( k , I , ZO ) 一 1(1( I test (k, I , ZO ) ) € SI )， 其中，k，l代表像素点在切片中的位置，Ite3St(k，l，ZQ)表示切片ZQ中k，l对应位置处的像 素点，K*)表示像素点*的灰度值，对图像g2做腐蚀运算，即获取第ZO层切片非肝脏区域的确 信像素点;对图像gl·做腐蚀运算并找到最大连通区域Π ，再做腐蚀运算后的像素点即为第Z0 层切片肝脏区域的确信像素点， 对于除ZO层的其他切片，利用相邻切片间的连接信息确定肝脏区域，首先，利用相邻切 片的最大连通区域rl作为先验信息，找出包含rl的最小矩形区域并向外延拓八个像素;然 后，在此区域上重复上述步骤，便可自动选取所有切片中肝脏与非肝脏区域的确信像素点 作为随机游走的种子点seedin、seed out， 将随机游走推广至特征空间，结合图像空间的先验信息并用于图像的特征空间中，基 于待分割图像对应的肝脏概率图，利用改进的随机游走算法确定肝脏的边界，考虑到在原 始待分割图像中，对于肝脏边缘模糊区域，像素点的对比度较低，难以识别肝脏区域，通过 利用肝脏概率图的信息，首先将其映射为带权的无向图G=(V，E)，结点集V由肝脏概率图P 的像素点集构成，边集E由邻接像素关系构成，其权值由像素点的先验概率决定：其中，i，j为待分割图像ltd中的像素索引，为两像素间的相邻关系，β为调整参数，P (i)和P(j)分别为索引为i的像素的先验概率和索引为j的像素的先验概率，然后计算每个 未标记像素点到达肝脏区域却；占。。。丄-的擗索η . s.?.q(seedin) - I，q(seed〇ut) - 0， 其中，q (seedin)表示seedin的先验概率，q (seedout)表示seedout的先验概率。 最后判断此概率，大于或者等于1/2时，像素点的类别label被判为1，即属于肝脏区域； 否则，属于非肝脏，最终实现肝_&#_ ΗΠ ·7. 根据权利要求6所述的基于特征学习的CT图像肝脏分割方法，其特征在于，对图像沿 做腐蚀运算所使用的腐蚀参数为6,对图像 gl做的两次腐蚀运算所使用的腐蚀参数为4。8. -种基于特征学习的CT图像肝脏分割系统，其特征在于，包括： 读取模块，用于读取训练图像集和待分割图像，其中，所述训练图像集中的训练图像和 待分割图像为腹部的CT图像； 提取模块，用于提取所述训练图像和待分割图像的Haar特征、局部二进制模式特征、方 向梯度直方图特征和共生矩阵特征； 融合模块，用于利用主成分分析方法将提取的所有特征进行特征融合，获取更有效的 特征； 分类模块，用于利用分类器对所述待分割图像每个像素的特征进行分类，得到肝脏概 率图； 分割模块，用于结合所述肝脏概率图与待分割图像，修改随机游走分割算法的图模型 权值，实现肝脏的分割。
【文档编号】G06T7/00GK105894517SQ201610256397
【公开日】2016年8月24日
【申请日】2016年4月22日
【发明人】艾丹妮, 杨健, 王涌天, 丛伟建, 张盼
【申请人】北京理工大学