CN107146229B - 基于元胞自动机模型的结肠息肉图像分割方法 - Google Patents
基于元胞自动机模型的结肠息肉图像分割方法 Download PDFInfo
- Publication number
- CN107146229B CN107146229B CN201710217480.0A CN201710217480A CN107146229B CN 107146229 B CN107146229 B CN 107146229B CN 201710217480 A CN201710217480 A CN 201710217480A CN 107146229 B CN107146229 B CN 107146229B
- Authority
- CN
- China
- Prior art keywords
- cellular
- indicate
- sigmoidoscope
- automation model
- polypiferous
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10068—Endoscopic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30028—Colon; Small intestine
- G06T2207/30032—Colon polyp
Abstract
一种基于元胞自动机模型的结肠息肉图像分割方法,主要解决现有结肠息肉图像分割技术中分割效率低、可重复性差、分割精度低的问题。本发明的技术方案为:(1)读入一幅含有息肉的结肠镜彩色图像;(2)修复图像中的高光区域;(3)初步检测图像中的息肉区域;(4)标记种子像素;(5)构建元胞自动机模型;(6)初始化元胞自动机模型;(7)进行图像分割;(8)输出分割后的图像。本发明利用息肉形状近似椭圆的先验知识,自动标记种子像素,构建元胞自动机模型,通过制定的局部转化规则进行图像分割,充分利用了图像的局部信息,具有分割效率、准确率高的优点,可用于结肠息肉图像的自动分割。
Description
技术领域
本发明属于图像处理技术领域,更进一步涉及医学图像分割技术领域的一种基于元胞自动机CA(Cellular Automata)的结肠息肉图像分割方法。本发明可用于对结肠镜图像中的息肉区域进行分割,以便实现对息肉的分类,判断息肉是否需要切除,进而辅助医生进行结肠镜检查和诊断。
背景技术
结直肠癌是全球第三大最常见的癌症类型。预防结肠癌的关键是检测结肠内潜在的腺瘤性息肉。因此,在结肠镜检查中需要设法只去除结肠上的腺瘤性息肉而保留良性息肉。结肠镜图像中息肉区域的分割是判断息肉类型的基本前提。分割的性能越好,分类的准确性越高,更有利于辅助医生进行结肠镜检查和诊断。
目前,常用的结肠息肉图像分割方法包括人工分割、半自动分割与全自动分割三种方式。人工分割虽然精确度高,但其不仅费时费力,而且分割结果很大程度上取决于医生的经验,一致性和可重复性差。半自动分割方式允许人工交互式操作提供一些有用的信息,然后由计算机进行分割处理。半自动分割通常比较准确,但仍在一定程度上受限于医生的水平。全自动分割完全利用计算机对图像进行分割,不仅能够提高分割的效率,而且分割的准确性更高。因此实现结肠息肉图像分割的主要手段就是设计计算机可执行的算法,进行全自动的分割。在结肠息肉图像的全自动分割算法中,基于知识的分割算法被广泛接受。这类算法在分割时利用了对目标的先验知识,包括专家定义的规则、理想目标的形状、外观模型等高层次知识,适合于对具有特定属性的具体病灶区域的分割。
重庆邮电大学在其申请的专利文献“一种基于云模型和图割相结合的医学图像分割方法”(申请号201510808682.3,公开号CN105389821A)中提出了一种基于云模型和图割相结合的医学图像分割方法。该方法的主要实现步骤是:首先手动标注前景区域;然后对图像进行去噪处理;接着利用逆向云变换分别计算出图像前景和背景的云特征值,利用X条件云发生器计算出每个像素对于前景和背景的隶属度,并计算出数据项和光滑项;最后建立能量函数构造相应网络图,并用最大流/最小割算法实现医学图像分割。该方法存在的不足之处是,需要手动标注前景区域,分割操作在一定程度上依赖于操作者的经验,可重复性差,执行效率低。
Seung-Hwan Bae等人在其发表的论文“Polyp Detection via ImbalancedLearning and Discriminative Feature Learning”(IEEE TRANSACTIONS ON MEDICALIMAGING 34(2015)2379-2393)中提出了一种自动检测内窥镜图像中息肉区域的方法。该方法的主要思想是对大量的训练数据进行特征提取和特征学习,训练出一个检测器,最后利用训练出的检测器来检测内窥镜图像中的息肉区域。该方法存在的不足之处是,检测器对内窥镜图像中息肉区域与非息肉区域的颜色和纹理这些传统特征比较敏感,而且检测结果没有准确地描绘出内窥镜图像中息肉区域的边缘。
发明内容
本发明的目的在于克服上述已有技术的不足,提出了一种基于元胞自动机模型的结肠息肉图像分割方法。
实现本发明目的的思路是,根据结肠息肉图像中拟分割息肉区域的形状近似椭圆的先验知识,在结肠息肉图像的多层级区域树中,选取与椭圆形状最近似的区域,引导种子点的自动标记,使分割操作成为一个全自动的过程,不再需要人工交互,构建的元胞自动机模型充分利用了结肠息肉图像的局部信息,通过元胞之间的局部竞争,实现对结肠息肉图像的分割。
本发明实现的具体步骤包括如下:
(1)读入一幅待处理的含有息肉的结肠镜彩色图像;
(2)修复含有息肉的结肠镜彩色图像中的高光区域:
(2a)利用双阈值法,从含有息肉的结肠镜彩色图像中检测高光区域;
(2b)利用基于样本的图像修复方法,对含有息肉的结肠镜彩色图像中的高光区域,进行修复;
(3)初步检测含有息肉的结肠镜彩色图像中的息肉区域:
(3a)利用gPb-owt-ucm算法,对修复后的含有息肉的结肠镜彩色图像进行层级分割,得到多层级区域树;
(3b)对多层级区域树中的每一个区域进行椭圆拟合操作,得到拟合椭圆;
(3c)计算多层级区域树中每一个区域与其拟合椭圆的重叠度;
(3d)从多层级区域树中的所有区域中选取重叠度最大的区域,作为初级息肉区域;
(4)在初级息肉区域上随机选取5%的像素作为种子像素;
(5)构建元胞自动机模型:
(5a)按照按照每行从左到右、依次逐行将含有息肉的结肠镜彩色图像中的每个像素映射为元胞自动机模型中的每个元胞;
(5b)将摩尔邻域作为元胞自动机模型的邻域;
(5c)将由标号、力量、强度向量三个状态分量组成的状态,作为元胞自动机模型中元胞的状态;
(6)初始化元胞自动机模型:
(6a)对元胞自动机模型中每个种子元胞的状态进行初始化;
(6b)对元胞自动机模型中每个非种子元胞的状态进行初始化:
(7)对修复后的含有息肉的结肠镜彩色图像进行分割:
(7a)按照下式,计算元胞自动机模型中每个元胞属于拟分割息肉区域的概率:
其中,μ(p)表示元胞自动机模型中第p个元胞属于拟分割息肉区域的概率,H表示元胞自动机模型中第p个元胞在多层级区域树中所属的区域,K表示初级息肉区域,D表示H与K取并操作后的区域的拟合椭圆;
(7b)按照下式,计算元胞自动机模型中每个元胞属于非拟分割息肉区域的概率:
δ(p)=1-μ(p)
其中,δ(p)表示元胞自动机模型中第p个元胞属于非拟分割息肉区域的概率,μ(p)表示元胞自动机模型中第p个元胞属于拟分割息肉区域的概率;
(7c)按照下式,计算元胞自动机模型中每个元胞的模糊力量:
其中,表示元胞自动机模型中第p个元胞在第t时刻的模糊力量,δ'(p)表示元胞自动机模型中第p个元胞属于非感兴趣区域的概率,μ'(p)表示元胞自动机模型中第p个元胞属于感兴趣区域的概率,表示元胞自动机模型中第p个元胞在第t时刻的力量;
(7d)按照下式,计算元胞自动机模型中每个元胞的每个邻居元胞的模糊标号:
其中,表示元胞自动机模型中第p个元胞的第q个邻居元胞在第t时刻的模糊标号,表示元胞自动机模型中第p个元胞在第t时刻的标号,δ(p,q)表元胞自动机模型中第p个元胞的第q个邻居元胞属于非感兴趣区域的概率,μ(p,q)表示元胞自动机模型中第p个元胞的第q个邻居元胞属于感兴趣区域的概率,表示元胞自动机模型中第p个元胞的第q个邻居元胞在第t时刻的标号;
(7e)按照下式,计算元胞自动机模型中每个元胞的每个邻居元胞的攻击力量:
其中,表示元胞自动机模型中第p个元胞的第q个邻居元胞在第t时刻的攻击力量,Cp表示元胞自动机模型中第p个元胞的强度向量,Cp,q表示元胞自动机模型中第p个元胞的第q个邻居元胞的强度向量,||·||2表示求欧拉距离操作,表示元胞自动机模型中第p个元胞的第q个邻居元胞在第t时刻的模糊力量;
(7f)按照下式,计算元胞自动机模型中每个元胞被其每个邻居元胞攻击的结果:
其中,r表示元胞自动机模型中第p个元胞被其第q个邻居元胞攻击的结果,表示元胞自动机模型中第p个元胞的第q个邻居元胞在第t时刻的攻击力量,表示元胞自动机模型中第p个元胞在第t时刻的模糊力量;
(7g)将元胞自动机模型中满足r>0的每个元胞中,第t时刻的模糊标号赋值为第t+1时刻的标号,将第t时刻的攻击力量赋值为其第t+1时刻的力量;
(7h)将元胞自动机模型中不满足r>0的每个元胞中,第t时刻的标号赋值为第t+1时刻的标号,将第t时刻的力量赋值为第t+1时刻的力量;
(7i)判断元胞自动机模型中是否存在r>0的元胞,若是,则将时刻t加1后执行步骤(7c),否则,执行步骤(8);
(8)输出分割后的图像。
本发明与现有技术相比具有如下优点:
第一,本发明利用结肠息肉图像中拟分割息肉区域的形状近似椭圆的先验知识,结合结肠息肉图像的多层级区域树,获得初级息肉区域,在初级息肉区域上随机选取5%的种子像素,自动标注前景区域,克服了现有技术需要手动标注前景区域,分割操作在一定程度上依赖于操作者的经验,可重复性差,执行效率低的缺点,使得本发明提高了分割效率和可重复性。
第二,本发明利用元胞自动机对结肠息肉图像分割任务进行建模,将摩尔邻域作为元胞自动机模型的邻域,充分利用结肠息肉图像的局部信息,克服了现有技术没有准确地描绘出结肠息肉图像中息肉区域的边缘的缺点,使得本发明能够得到更加准确的分割结果。
第三,本发明利用元胞自动机对结肠息肉图像分割任务进行建模,元胞自动机模型中的分割规则不仅利用了结肠息肉图像的颜色特征,而且还利用了结肠息肉图像的形状特征,克服了现有技术对颜色传统特征比较敏感的缺点,使得本发明具有更高的鲁棒性。
附图说明
图1是本发明的流程图。
具体实施方式
下面结合附图对本发明作进一步的描述。
参照图1,本发明的具体实现步骤如下。
步骤1,读入一幅待处理的含有息肉的结肠镜彩色图像。
步骤2,修复含有息肉的结肠镜彩色图像中的高光区域。
利用双阈值法,从含有息肉的结肠镜彩色图像中检测高光区域。
所述双阈值法的具体步骤如下:
第1步,按照下式,将输入的含有息肉的结肠镜彩色图像,从红绿蓝RGB颜色空间转换到明度饱和度色调HSV颜色空间:
V←max(R,G,B)
其中,V表示输入的含有息肉的结肠镜彩色图像在HSV颜色空间上的明度分量,←表示赋值操作,max表示求最大值操作,R表示输入的含有息肉的结肠镜彩色图像在RGB颜色空间上的红色分量,G表示输入的含有息肉的结肠镜彩色图像在RGB颜色空间上的绿色分量,B表示输入的含有息肉的结肠镜彩色图像在RGB颜色空间上的蓝色分量,S表示输入的含有息肉的结肠镜彩色图像在HSV颜色空间上的饱和度分量,min表示求最小值操作,otherwise表示条件的补集,H表示输入的含有息肉的结肠镜彩色图像在HSV颜色空间上的色调分量,and表示对条件取交操作。
第2步,按照下式,对输入的含有息肉的结肠镜彩色图像中的每个像素进行高阈值处理,得到高阈值处理后的二值图像:
其中,f(x)表示对输入的含有息肉的结肠镜彩色图像中第x个像素进行高阈值处理后的二值图像,S表示输入的含有息肉的结肠镜彩色图像在HSV颜色空间上的饱和度分量,and表示对条件取交操作,V表示输入的含有息肉的结肠镜彩色图像在HSV颜色空间上的明度分量,otherwise表示条件的补集。
第3步,用半径为5个像素的圆盘形结构元素,对高阈值处理后的二值图像进行膨胀操作,得到膨胀后的二值图像。
第4步,按照下式,对输入的含有息肉的结肠镜彩色图像中的每个像素进行低阈值处理,得到低阈值处理后的二值图像:
其中,f′(y)表示对输入的含有息肉的结肠镜彩色图像中第y个像素进行低阈值处理后的二值图像,S表示输入的含有息肉的结肠镜彩色图像在HSV颜色空间上的饱和度分量,and表示对条件取交操作,V表示输入的含有息肉的结肠镜彩色图像在HSV颜色空间上的明度分量,otherwise表示条件的补集。
第5步,将能使膨胀后的二值图像中的两个不连通的区域变为连通,且值为1的像素称为高像素,将膨胀后的二值图像中与低阈值处理后的二值图像中的高像素位置相同的像素的值置为1,得到双阈值处理后的二值图像,将双阈值处理后的二值图像中值为1的像素组成高光区域。
利用基于样本的图像修复方法,对含有息肉的结肠镜彩色图像中的高光区域,进行修复。
所述基于样本的图像修复方法的具体步骤如下:
第1步,将输入的含有息肉的结肠镜彩色图像中的高光区域作为目标区域,将位于目标区域边缘上所有的像素点组成边缘像素点集合,从输入的含有息肉的结肠镜彩色图像中去除目标区域,得到源区域。
第2步,依次以边缘像素点集合中的每个像素点为中心,在输入的含有息肉的结肠镜彩色图像上选取边长为9个像素的正方形结构,得到待填充块集合。
第3步,按照下式,计算待填充块集合中每个待填充块的信任值:
其中,B(Ψ)表示待填充块集合中第Ψ个待填充块的信任值,∑表示求和操作,q表示叠加块V上的第q个像素,V表示待填充块集合中第Ψ个待填充块与源区域取交操作后的叠加块,∩表示交操作,∈表示属于符号,Cq表示叠加区域上第q个像素的置信度,|·|表示取绝对值操作。
第4步,按照下式,计算待填充块集合中每个待填充块的数据值:
其中,D(Ψ)表示待填充块集合中第Ψ个待填充块的数据值,IΨ表示待填充块集合中第Ψ个待填充块的中心像素的等照度线,·表示点乘操作,nΨ表示待填充块集合中第Ψ个待填充块的中心像素在目标区域边缘上的单位法向量,α表示图像标准化因子,α取值为255。
第5步,按照下式,计算待填充块集合中每个待填充块的优先级:
P(Ψ)=B(Ψ)D(Ψ)
其中,P(Ψ)表示待填充块集合中第Ψ个待填充块的优先级,B(Ψ)表示待填充块集合中第Ψ个待填充块的信任值,D(Ψ)表示待填充块集合中第Ψ个待填充块的数据值。
第6步,从待填充块集合中选出优先级最大的待填充块,得到优先待填充块;
第7步,将与优先待填充块距离最小的块,作为优先待填充块的最佳匹配块,在源区域上查找优先待填充块的最佳匹配块;
第8步,将优先待填充块的最佳匹配块复制到优先待填充块;
第9步,从目标区域中去掉优先待填充块;
第10步,判断目标区域是否为空,若是,则执行本步骤的第1步,否则,得到高光区域修复后的图像。
步骤3,初步检测含有息肉的结肠镜彩色图像中的息肉区域。
利用gPb-owt-ucm算法,对修复后的含有息肉的结肠镜彩色图像进行层级分割,得到多层级区域树。
所述gPb-owt-ucm算法的具体步骤如下:
第1步,对含有息肉的结肠镜彩色图像进行特征分解操作,得到明度通道、颜色A通道、颜色B通道和纹理通道。
第2步,按照下式,计算含有息肉的结肠镜彩色图像中每个像素的边界局部信号,得到边界局部信号矩阵:
其中,mPb(z)表示对含有息肉的结肠镜彩色图像中第z个像素计算边界局部信号后,由所有像素的边界局部信号组成的边界局部信号矩阵,max{·}表示取最大值操作,∑表示求和操作,s表示尺度值,i表示特征通道,αi,s表示直方图权重,表示以含有息肉的结肠镜彩色图像中第z个像素为中心以θ为倾斜角以σi,s为直径的两个半圆的直方图。
第3步,按照下式,计算每一个元素的稀疏值,由所有稀疏值组成稀疏对称相似矩阵:
其中,Wi,j表示稀疏对称相似矩阵中第i行第j列的元素的稀疏值,0≤i<m,0≤j<n,m表示含有息肉的结肠镜彩色图像的长度,n表示含有息肉的结肠镜彩色图像的宽度,max{·}表示取最大值操作,exp表示以自然常数e为底的指数操作,a表示含有息肉的结肠镜彩色图像中第m个像素与第n个像素相连的直线l(u,v)上的任意一点,∈表示属于符号,u的取值与i相等,v的取值与j相等,mPb(z)表示对含有息肉的结肠镜彩色图像中第z个像素计算边界局部信号后,由所有像素的边界局部信号组成的边界局部信号矩阵,ρ为0.1的常数。
第4步,按照下式,计算每个元素的边界全局信号,将所有元素的边界全局信号组成边界全局信号矩阵:
其中,sPb(z)表示对含有息肉的结肠镜彩色图像中第z个像素计算边界全局信号后,由所有像素的边界全局信号组成的边界全局信号矩阵,n取值为16,表示开根号操作,λk表示稀疏对称相似矩阵的第k个特征值,▽表示求梯度操作,V(k)表示稀疏对称相似矩阵的第k个特征向量。
第5步,按照下式,计算含有息肉的结肠镜彩色图像中每个像素的边界概率,得到边界概率矩阵:
其中,gPb(z)表示对含有息肉的结肠镜彩色图像中第z个像素计算边界概率后,由所有像素的边界概率组成的边界概率矩阵,β表示边界局部信号权重,mPb(z)表示边界局部信号矩阵,γ表示边界全局信号权重,sPb(z)表示边界全局信号矩阵;
第6步,对边界概率矩阵做方向分水岭变换,得到最大细节图像。
第7步,按照下式,将最大细节图像映射到图数据结构:
G=(R,E,W)
其中,G表示由最大细节图像映射得到的图数据结构,R表示最大细节图像中所有区域的集合,E表示最大细节图像中所有边缘的集合,W表示边缘度。
第8步,按照边缘度从小到大的顺序,依次将图数据结构中的每个边缘去除,每次去除操作产生一个层级区域,当图数据结构中的所有边缘都被去除时,得到多层级区域树。
对多层级区域树中的每一个区域进行椭圆拟合操作,得到拟合椭圆。
计算多层级区域树中每一个区域与其拟合椭圆的重叠度。
所述的多层级区域树中每一个区域与其拟合椭圆的重叠度是按照下式计算的:
其中,f(A,B)表示多层次区域树中的区域A与拟合椭圆B的重叠度,|·|表示取绝对值操作,∩表示交操作,∪表示并操作。
从多层级区域树中的所有区域中选取重叠度最大的区域,作为初级息肉区域。
步骤4,在初级息肉区域上随机选取5%的像素作为种子像素。
步骤5,构建元胞自动机模型:
按照每行从左到右、依次逐行将含有息肉的结肠镜彩色图像中的每个像素映射为元胞自动机模型中的每个元胞。
将摩尔邻域作为元胞自动机模型的邻域。
将由标号、力量、强度向量三个状态分量组成的状态,作为元胞自动机模型中元胞的状态。
步骤6,初始化元胞自动机模型。
对元胞自动机模型中每个种子元胞的状态进行初始化。
所述的对元胞自动机模型中每个种子元胞状态的初始化是按照下式操作的:
lo=1,θo=1,Co=RGBo
其中,lo表示元胞自动机模型中第o个元胞的标号;θo表示元胞自动机模型中第o个元胞的力量,θo的取值范围是[0,1]上的所有实数,Co表示元胞自动机模型中第o个元胞的强度向量,RGBo表示由含有息肉的结肠镜彩色图像中第o个像素的红色分量R、绿色分量G、蓝色分量B组成的向量。
对元胞自动机模型中每个非种子元胞的状态进行初始化。
所述的对元胞自动机模型中每个非种子元胞状态的初始化是按照下式操作的:
lb=0,θb=0,Cb=RGBb
其中,lb表示元胞自动机模型中第b个元胞的标号;θb表示元胞自动机模型中第b个元胞的力量,θb的取值范围是[0,1]上的所有实数,Cb表示元胞自动机模型中第b个元胞的强度向量,RGBb表示由含有息肉的结肠镜彩色图像中第b个像素的红色分量R、绿色分量G、蓝色分量B组成的向量。
步骤7,对修复后的含有息肉的结肠镜彩色图像进行分割。
第1步,按照下式,计算元胞自动机模型中每个元胞属于拟分割息肉区域的概率:
其中,μ(p)表示元胞自动机模型中第p个元胞属于拟分割息肉区域的概率,H表示元胞自动机模型中第p个元胞在多层级区域树中所属的区域,K表示初级息肉区域,D表示H与K取并操作后的区域的拟合椭圆。
第2步,按照下式,计算元胞自动机模型中每个元胞属于非拟分割息肉区域的概率:
δ(p)=1-μ(p)
其中,δ(p)表示元胞自动机模型中第p个元胞属于非拟分割息肉区域的概率,μ(p)表示元胞自动机模型中第p个元胞属于拟分割息肉区域的概率。
第3步,按照下式,计算元胞自动机模型中每个元胞的模糊力量:
其中,表示元胞自动机模型中第p个元胞在第t时刻的模糊力量,δ'(p)表示元胞自动机模型中第p个元胞属于非感兴趣区域的概率,μ'(p)表示元胞自动机模型中第p个元胞属于感兴趣区域的概率,表示元胞自动机模型中第p个元胞在第t时刻的力量。
第4步,按照下式,计算元胞自动机模型中每个元胞的每个邻居元胞的模糊标号:
其中,表示元胞自动机模型中第p个元胞的第q个邻居元胞在第t时刻的模糊标号,表示元胞自动机模型中第p个元胞在第t时刻的标号,δ(p,q)表元胞自动机模型中第p个元胞的第q个邻居元胞属于非感兴趣区域的概率,μ(p,q)表示元胞自动机模型中第p个元胞的第q个邻居元胞属于感兴趣区域的概率,表示元胞自动机模型中第p个元胞的第q个邻居元胞在第t时刻的标号。
第5步,按照下式,计算元胞自动机模型中每个元胞的每个邻居元胞的攻击力量:
其中,表示元胞自动机模型中第p个元胞的第q个邻居元胞在第t时刻的攻击力量,Cp表示元胞自动机模型中第p个元胞的强度向量,Cp,q表示元胞自动机模型中第p个元胞的第q个邻居元胞的强度向量,||·||2表示求欧拉距离操作,表示元胞自动机模型中第p个元胞的第q个邻居元胞在第t时刻的模糊力量。
第6步,按照下式,计算元胞自动机模型中每个元胞被其每个邻居元胞攻击的结果:
其中,r表示元胞自动机模型中第p个元胞被其第q个邻居元胞攻击的结果,表示元胞自动机模型中第p个元胞的第q个邻居元胞在第t时刻的攻击力量,表示元胞自动机模型中第p个元胞在第t时刻的模糊力量。
第7步,将元胞自动机模型中满足r>0的每个元胞中,第t时刻的模糊标号赋值为第t+1时刻的标号,将第t时刻的攻击力量赋值为其第t+1时刻的力量。
第8步,将元胞自动机模型中不满足r>0的每个元胞中,第t时刻的标号赋值为第t+1时刻的标号,将第t时刻的力量赋值为第t+1时刻的力量。
第9步,判断元胞自动机模型中是否存在r>0的元胞,若是,则将时刻t加1后执行本步骤的第3步,否则,执行步骤8。
步骤8,输出分割后的图像。
Claims (6)
1.一种基于元胞自动机模型的结肠息肉图像分割方法,包括如下步骤:
(1)读入一幅待处理的含有息肉的结肠镜彩色图像;
(2)修复含有息肉的结肠镜彩色图像中的高光区域:
(2a)利用双阈值法,从含有息肉的结肠镜彩色图像中检测高光区域;
(2b)利用基于样本的图像修复方法,对含有息肉的结肠镜彩色图像中的高光区域,进行修复;
(3)初步检测含有息肉的结肠镜彩色图像中的息肉区域:
(3a)利用gPb-owt-ucm算法,对修复后的含有息肉的结肠镜彩色图像进行层级分割,得到多层级区域树;
所述gPb-owt-ucm算法的具体步骤如下:
第一步,对含有息肉的结肠镜彩色图像进行特征分解操作,得到明度通道、颜色A通道、颜色B通道和纹理通道;
第二步,按照下式,计算含有息肉的结肠镜彩色图像中每个像素的边界局部信号,得到边界局部信号矩阵:
其中,mPb(z)表示对含有息肉的结肠镜彩色图像中第z个像素计算边界局部信号后,由所有像素的边界局部信号组成的边界局部信号矩阵,max{·}表示取最大值操作,∑表示求和操作,s表示尺度值,i表示特征通道,αi,s表示直方图权重,表示以含有息肉的结肠镜彩色图像中第z个像素为中心以θ为倾斜角以σi,s为直径的两个半圆的直方图;
第三步,按照下式,计算每一个元素的稀疏值,由所有稀疏值组成稀疏对称相似矩阵:
其中,Wi,j表示稀疏对称相似矩阵中第i行第j列的元素的稀疏值,0≤i<m,0≤j<n,m表示含有息肉的结肠镜彩色图像的长度,n表示含有息肉的结肠镜彩色图像的宽度,max{·}表示取最大值操作,exp表示以自然常数e为底的指数操作,a表示含有息肉的结肠镜彩色图像中第m个像素与第n个像素相连的直线l(u,v)上的任意一点,∈表示属于符号,u的取值与i相等,v的取值与j相等,mPb(z)表示对含有息肉的结肠镜彩色图像中第z个像素计算边界局部信号后,由所有像素的边界局部信号组成的边界局部信号矩阵,ρ为0.1的常数;
第四步,按照下式,计算每个元素的边界全局信号,将所有元素的边界全局信号组成边界全局信号矩阵:
其中,sPb(z)表示对含有息肉的结肠镜彩色图像中第z个像素计算边界全局信号后,由所有像素的边界全局信号组成的边界全局信号矩阵,n取值为16,表示开根号操作,λk表示稀疏对称相似矩阵的第k个特征值,▽表示求梯度操作,V(k)表示稀疏对称相似矩阵的第k个特征向量;
第五步,按照下式,计算含有息肉的结肠镜彩色图像中每个像素的边界概率,得到边界概率矩阵:
其中,gPb(z)表示对含有息肉的结肠镜彩色图像中第z个像素计算边界概率后,由所有像素的边界概率组成的边界概率矩阵,β表示边界局部信号权重,mPb(z)表示边界局部信号矩阵,γ表示边界全局信号权重,sPb(z)表示边界全局信号矩阵;
第六步,对边界概率矩阵做方向分水岭变换,得到最大细节图像;
第七步,按照下式,将最大细节图像映射到图数据结构:
G=(R,E,W)
其中,G表示由最大细节图像映射得到的图数据结构,R表示最大细节图像中所有区域的集合,E表示最大细节图像中所有边缘的集合,W表示边缘度;
第八步,按照边缘度从小到大的顺序,依次将图数据结构中的每个边缘去除,每次去除操作产生一个层级区域,当图数据结构中的所有边缘都被去除时,得到多层级区域树;
(3b)对多层级区域树中的每一个区域进行椭圆拟合操作,得到拟合椭圆;
(3c)计算多层级区域树中每一个区域与其拟合椭圆的重叠度;
(3d)从多层级区域树中的所有区域中选取重叠度最大的区域,作为初级息肉区域;
(4)在初级息肉区域上随机选取5%的像素作为种子像素;
(5)构建元胞自动机模型:
(5a)按照每行从左到右,依次逐行将含有息肉的结肠镜彩色图像中的每个像素映射为元胞自动机模型中的每个元胞;
(5b)将摩尔邻域作为元胞自动机模型的邻域;
(5c)将由标号、力量、强度向量三个状态分量组成的状态,作为元胞自动机模型中元胞的状态;
(6)初始化元胞自动机模型:
(6a)对元胞自动机模型中每个种子元胞的状态进行初始化;
(6b)对元胞自动机模型中每个非种子元胞的状态进行初始化:
(7)对修复后的含有息肉的结肠镜彩色图像进行分割:
(7a)按照下式,计算元胞自动机模型中每个元胞属于拟分割息肉区域的概率:
其中,μ(p)表示元胞自动机模型中第p个元胞属于拟分割息肉区域的概率,H表示元胞自动机模型中第p个元胞在多层级区域树中所属的区域,K表示初级息肉区域,D表示H与K取并操作后的区域的拟合椭圆;
(7b)按照下式,计算元胞自动机模型中每个元胞属于非拟分割息肉区域的概率:
δ(p)=1-μ(p)
其中,δ(p)表示元胞自动机模型中第p个元胞属于非拟分割息肉区域的概率,μ(p)表示元胞自动机模型中第p个元胞属于拟分割息肉区域的概率;
(7c)按照下式,计算元胞自动机模型中每个元胞的模糊力量:
其中,表示元胞自动机模型中第p个元胞在第t时刻的模糊力量,δ'(p)表示元胞自动机模型中第p个元胞属于非感兴趣区域的概率,μ'(p)表示元胞自动机模型中第p个元胞属于感兴趣区域的概率,表示元胞自动机模型中第p个元胞在第t时刻的力量;
(7d)按照下式,计算元胞自动机模型中每个元胞的每个邻居元胞的模糊标号:
其中,表示元胞自动机模型中第p个元胞的第q个邻居元胞在第t时刻的模糊标号,表示元胞自动机模型中第p个元胞在第t时刻的标号,δ(p,q)表元胞自动机模型中第p个元胞的第q个邻居元胞属于非感兴趣区域的概率,μ(p,q)表示元胞自动机模型中第p个元胞的第q个邻居元胞属于感兴趣区域的概率,表示元胞自动机模型中第p个元胞的第q个邻居元胞在第t时刻的标号;
(7e)按照下式,计算元胞自动机模型中每个元胞的每个邻居元胞的攻击力量:
其中,表示元胞自动机模型中第p个元胞的第q个邻居元胞在第t时刻的攻击力量,Cp表示元胞自动机模型中第p个元胞的强度向量,Cp,q表示元胞自动机模型中第p个元胞的第q个邻居元胞的强度向量,||·||2表示求欧拉距离操作,表示元胞自动机模型中第p个元胞的第q个邻居元胞在第t时刻的模糊力量;
(7f)按照下式,计算元胞自动机模型中每个元胞被其每个邻居元胞攻击的结果:
其中,r表示元胞自动机模型中第p个元胞被其第q个邻居元胞攻击的结果,表示元胞自动机模型中第p个元胞的第q个邻居元胞在第t时刻的攻击力量,表示元胞自动机模型中第p个元胞在第t时刻的模糊力量;
(7g)将元胞自动机模型中满足r>0的每个元胞中,第t时刻的模糊标号赋值为第t+1时刻的标号,将第t时刻的攻击力量赋值为其第t+1时刻的力量;
(7h)将元胞自动机模型中不满足r>0的每个元胞中,第t时刻的标号赋值为第t+1时刻的标号,将第t时刻的力量赋值为第t+1时刻的力量;
(7i)判断元胞自动机模型中是否存在r>0的元胞,若是,则将时刻t加1后执行步骤(7c),否则,执行步骤(8);
(8)输出分割后的图像。
2.根据权利要求1所述的基于元胞自动机模型的结肠息肉图像分割方法,其特征在于,步骤(2a)中所述的双阈值法的具体步骤如下:
第一步,按照下式,将输入的含有息肉的结肠镜彩色图像,从红绿蓝RGB颜色空间转换到明度饱和度色调HSV颜色空间:
V←max(R,G,B)
其中,V表示输入的含有息肉的结肠镜彩色图像在HSV颜色空间上的明度分量,←表示赋值操作,max表示求最大值操作,R表示输入的含有息肉的结肠镜彩色图像在RGB颜色空间上的红色分量,G表示输入的含有息肉的结肠镜彩色图像在RGB颜色空间上的绿色分量,B表示输入的含有息肉的结肠镜彩色图像在RGB颜色空间上的蓝色分量,S表示输入的含有息肉的结肠镜彩色图像在HSV颜色空间上的饱和度分量,min表示求最小值操作,otherwise表示条件的补集,H表示输入的含有息肉的结肠镜彩色图像在HSV颜色空间上的色调分量,and表示对条件取交操作;
第二步,按照下式,对输入的含有息肉的结肠镜彩色图像中的每个像素进行高阈值处理,得到高阈值处理后的二值图像:
其中,f(x)表示对输入的含有息肉的结肠镜彩色图像中第x个像素进行高阈值处理后的二值图像,S表示输入的含有息肉的结肠镜彩色图像在HSV颜色空间上的饱和度分量,and表示对条件取交操作,V表示输入的含有息肉的结肠镜彩色图像在HSV颜色空间上的明度分量,otherwise表示条件的补集;
第三步,用半径为5个像素的圆盘形结构元素,对高阈值处理后的二值图像进行膨胀操作,得到膨胀后的二值图像;
第四步,按照下式,对输入的含有息肉的结肠镜彩色图像中的每个像素进行低阈值处理,得到低阈值处理后的二值图像:
其中,f′(y)表示对输入的含有息肉的结肠镜彩色图像中第y个像素进行低阈值处理后的二值图像,S表示输入的含有息肉的结肠镜彩色图像在HSV颜色空间上的饱和度分量,and表示对条件取交操作,V表示输入的含有息肉的结肠镜彩色图像在HSV颜色空间上的明度分量,otherwise表示条件的补集;
第五步,将能使膨胀后的二值图像中的两个不连通的区域变为连通,且值为1的像素称为高像素,将膨胀后的二值图像中与低阈值处理后的二值图像中的高像素位置相同的像素的值置为1,得到双阈值处理后的二值图像,将双阈值处理后的二值图像中值为1的像素组成高光区域。
3.根据权利要求1所述的基于元胞自动机模型的结肠息肉图像分割方法,其特征在于,步骤(2b)中所述基于样本的图像修复方法的具体步骤如下:
第一步,将输入的含有息肉的结肠镜彩色图像中的高光区域作为目标区域,将位于目标区域边缘上所有的像素点组成边缘像素点集合,从输入的含有息肉的结肠镜彩色图像中去除目标区域,得到源区域;
第二步,依次以边缘像素点集合中的每个像素点为中心,在输入的含有息肉的结肠镜彩色图像上选取边长为9个像素的正方形结构,得到待填充块集合;
第三步,按照下式,计算待填充块集合中每个待填充块的信任值:
其中,B(Ψ)表示待填充块集合中第Ψ个待填充块的信任值,∑表示求和操作,q表示叠加块V上的第q个像素,V表示待填充块集合中第Ψ个待填充块与源区域取交操作后的叠加块,∩表示交操作,∈表示属于符号,Cq表示叠加区域上第q个像素的置信度,|·|表示取绝对值操作;
第四步,按照下式,计算待填充块集合中每个待填充块的数据值:
其中,D(Ψ)表示待填充块集合中第Ψ个待填充块的数据值,IΨ表示待填充块集合中第Ψ个待填充块的中心像素的等照度线,·表示点乘操作,nΨ表示待填充块集合中第Ψ个待填充块的中心像素在目标区域边缘上的单位法向量,α表示图像标准化因子,α取值为255;
第五步,按照下式,计算待填充块集合中每个待填充块的优先级:
P(Ψ)=B(Ψ)D(Ψ)
其中,P(Ψ)表示待填充块集合中第Ψ个待填充块的优先级,B(Ψ)表示待填充块集合中第Ψ个待填充块的信任值,D(Ψ)表示待填充块集合中第Ψ个待填充块的数据值;
第六步,从待填充块集合中选出优先级最大的待填充块,得到优先待填充块;
第七步,将与优先待填充块距离最小的块,作为优先待填充块的最佳匹配块,在源区域上查找优先待填充块的最佳匹配块;
第八步,将优先待填充块的最佳匹配块复制到优先待填充块;
第九步,从目标区域中去掉优先待填充块;
第十步,判断目标区域是否为空,若是,则执行第一步,否则,得到高光区域修复后的图像。
4.根据权利要求1所述的基于元胞自动机模型的结肠息肉图像分割方法,其特征在于,步骤(3c)中所述的多层级区域树中每一个区域与其拟合椭圆的重叠度是按照下式计算的:
其中,f(A,B)表示多层次区域树中的区域A与拟合椭圆B的重叠度,|·|表示取绝对值操作,∩表示交操作,∪表示并操作。
5.根据权利要求1所述的基于元胞自动机模型的结肠息肉图像分割方法,其特征在于,步骤(6a)中所述的对元胞自动机模型中每个种子元胞状态的初始化是按照下式操作的:
lo=1,θo=1,Co=RGBo
其中,lo表示元胞自动机模型中第o个元胞的标号;θo表示元胞自动机模型中第o个元胞的力量,θo的取值范围是[0,1]上的所有实数,Co表示元胞自动机模型中第o个元胞的强度向量,RGBo表示由含有息肉的结肠镜彩色图像中第o个像素的红色分量R、绿色分量G、蓝色分量B组成的向量。
6.根据权利要求1所述的基于元胞自动机模型的结肠息肉图像分割方法,其特征在于,步骤(6b)中所述的对元胞自动机模型中每个非种子元胞状态的初始化是按照下式操作的:
lb=0,θb=0,Cb=RGBb
其中,lb表示元胞自动机模型中第b个元胞的标号;θb表示元胞自动机模型中第b个元胞的力量,θb的取值范围是[0,1]上的所有实数,Cb表示元胞自动机模型中第b个元胞的强度向量,RGBb表示由含有息肉的结肠镜彩色图像中第b个像素的红色分量R、绿色分量G、蓝色分量B组成的向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710217480.0A CN107146229B (zh) | 2017-04-05 | 2017-04-05 | 基于元胞自动机模型的结肠息肉图像分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710217480.0A CN107146229B (zh) | 2017-04-05 | 2017-04-05 | 基于元胞自动机模型的结肠息肉图像分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107146229A CN107146229A (zh) | 2017-09-08 |
CN107146229B true CN107146229B (zh) | 2019-08-13 |
Family
ID=59773687
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710217480.0A Active CN107146229B (zh) | 2017-04-05 | 2017-04-05 | 基于元胞自动机模型的结肠息肉图像分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107146229B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107679551B (zh) * | 2017-09-11 | 2020-06-16 | 电子科技大学 | 一种基于分形的涌现现象的识别方法 |
CN108491773B (zh) * | 2018-03-12 | 2022-11-08 | 中国工商银行股份有限公司 | 一种识别方法及系统 |
CN110741410A (zh) * | 2018-04-25 | 2020-01-31 | 北京嘀嘀无限科技发展有限公司 | 用于确定图像中的反光区的系统和方法 |
CN110751605B (zh) * | 2019-10-16 | 2022-12-23 | 深圳开立生物医疗科技股份有限公司 | 一种图像处理方法、装置和电子设备及可读存储介质 |
CN112465766A (zh) * | 2020-11-25 | 2021-03-09 | 武汉楚精灵医疗科技有限公司 | 扁平、微小息肉图像识别方法 |
CN112489061B (zh) * | 2020-12-09 | 2024-04-16 | 浙江工业大学 | 一种基于多尺度信息与并行注意力机制的深度学习肠道息肉分割方法 |
CN112884702B (zh) * | 2020-12-29 | 2023-07-28 | 香港中文大学深圳研究院 | 一种基于内窥镜图像的息肉识别系统和方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102156979A (zh) * | 2010-12-31 | 2011-08-17 | 上海电机学院 | 基于GrowCut的车道线快速检测方法及系统 |
CN102254187A (zh) * | 2011-08-10 | 2011-11-23 | 武汉大学 | 一种基于元胞自动机的遥感影像非监督分类方法 |
CN102789631A (zh) * | 2011-05-17 | 2012-11-21 | 张学 | 基于图像分割和引力模型的遥感图像分辨率提高处理方法 |
CN103793894A (zh) * | 2013-12-04 | 2014-05-14 | 国家电网公司 | 基于云模型元胞自动机角点检测的变电站遥视图像拼接方法 |
CN105894502A (zh) * | 2016-03-30 | 2016-08-24 | 浙江大学 | 基于超图模型的rgbd图像显著性检测方法 |
-
2017
- 2017-04-05 CN CN201710217480.0A patent/CN107146229B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102156979A (zh) * | 2010-12-31 | 2011-08-17 | 上海电机学院 | 基于GrowCut的车道线快速检测方法及系统 |
CN102789631A (zh) * | 2011-05-17 | 2012-11-21 | 张学 | 基于图像分割和引力模型的遥感图像分辨率提高处理方法 |
CN102254187A (zh) * | 2011-08-10 | 2011-11-23 | 武汉大学 | 一种基于元胞自动机的遥感影像非监督分类方法 |
CN103793894A (zh) * | 2013-12-04 | 2014-05-14 | 国家电网公司 | 基于云模型元胞自动机角点检测的变电站遥视图像拼接方法 |
CN105894502A (zh) * | 2016-03-30 | 2016-08-24 | 浙江大学 | 基于超图模型的rgbd图像显著性检测方法 |
Non-Patent Citations (1)
Title |
---|
基于标记提取的GrowCut自动分割算法;杨耀 等;《信息技术》;20150525(第5期);第76-80页 |
Also Published As
Publication number | Publication date |
---|---|
CN107146229A (zh) | 2017-09-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107146229B (zh) | 基于元胞自动机模型的结肠息肉图像分割方法 | |
CN107133616B (zh) | 一种基于深度学习的无分割字符定位与识别方法 | |
CN104392228B (zh) | 基于条件随机场模型的无人机图像目标类检测方法 | |
CN106709517B (zh) | 红树林识别方法及系统 | |
CN108549891B (zh) | 基于背景与目标先验的多尺度扩散显著目标检测方法 | |
KR102102161B1 (ko) | 이미지 내 객체의 대표 특성을 추출하는 방법, 장치 및 컴퓨터 프로그램 | |
CN107274419A (zh) | 一种基于全局先验和局部上下文的深度学习显著性检测方法 | |
CN105893925A (zh) | 基于肤色的人手检测方法及装置 | |
CN109584248A (zh) | 基于特征融合和稠密连接网络的红外面目标实例分割方法 | |
CN108681692A (zh) | 一种基于深度学习的遥感图像内新增建筑识别方法 | |
CN108229288B (zh) | 神经网络训练及衣服颜色检测方法、装置、存储介质、电子设备 | |
CN105279519B (zh) | 基于协同训练半监督学习的遥感影像水体提取方法及系统 | |
CN106611423B (zh) | 基于脊波滤波器和反卷积结构模型的sar图像分割方法 | |
CN106611420B (zh) | 基于反卷积网络和素描图方向约束的sar图像分割方法 | |
CN103914699A (zh) | 一种基于色彩空间的自动唇彩的图像增强的方法 | |
CN108074243A (zh) | 一种细胞定位方法以及细胞分割方法 | |
CN106296695A (zh) | 基于显著性的自适应阈值自然目标图像分割抽取算法 | |
CN106611422B (zh) | 基于素描结构的随机梯度贝叶斯sar图像分割方法 | |
CN102024156B (zh) | 彩色人脸图像中的唇部区域定位方法 | |
CN106683102B (zh) | 基于脊波滤波器和卷积结构学习模型的sar图像分割方法 | |
CN106952271A (zh) | 一种基于超像素分割和em/mpm处理的图像分割方法 | |
CN103177446A (zh) | 基于邻域和非邻域平滑先验的图像前景精准提取方法 | |
CN106846322B (zh) | 基于曲线波滤波器和卷积结构学习的sar图像分割方法 | |
CN109543632A (zh) | 一种基于浅层特征融合引导的深层网络行人检测方法 | |
CN109166102A (zh) | 一种基于判别区域候选对抗网络的图像转图像翻译方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |