发明名称 一种基于广义数据场和Ncut算法的图像分割方法
摘要 本发明提供一种基于广义数据场和Ncut算法的图像分割方法,首先,将特征空间划分为层次网格,第一层网格和第二层网格,同时形成与之对应的网格特征空间Ω<sub>s</sub>和Ω<sub>b</sub>,第一层每8个相邻小网格形成第二层的一个大网格;然后通过运用GDF算法基于第二层网格计算得到第一层网格的势值分布。基于势值分布,再将第一层网格进行聚类,聚类结果映射到图像上,从而实现对于一幅图像的初始分割操作,将其划分为不同的彼此不相交的区域;最后,基于图像初始分割的结果构建无向加权图后,运用基于区域的Ncut算法来合并性质相同的区域,直至达到最佳图像分割结果。本发明在图像分割上具有快速、简单、准确的优点。
申请公布号 CN102855624B 申请公布日期 2015.09.16
申请号 CN201210265614.3 申请日期 2012.07.23
申请人 武汉大学 发明人 王树良;李英;尹进飞;陈其良;李伟
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 武汉科皓知识产权代理事务所(特殊普通合伙) 42222 代理人 薛玲
主权项 一种基于广义数据场和Ncut算法的图像分割方法,其特征在于:包括以下步骤,步骤1、层次网格的划分以及势值估计,具体包括以下步骤,步骤1.1、将图像的RGB特征空间转换为L*u*v*或L*a*b*特征空间,将L*u*v*或L*a*b*特征空间Ω划分为2N×2N×2N个小网格作为第一层网格,计算每个小网格内数据点的均值,并以此作为该小网格的特征值,形成一个新的特征空间Ω<sub>s</sub>;步骤1.2将八邻域小网格合并成为一个大网格作为第二层网格,形成一个新的特征空间Ω<sub>b</sub>及其相应的网格空间坐标值;步骤1.3根据势值估计公式计算第一层网格内每个小网格的势值<img file="FDA0000746302730000011.GIF" wi="128" he="53" /><img file="FDA0000746302730000012.GIF" wi="895" he="304" />在特征空间Ω<sub>b</sub>中,<img file="FDA0000746302730000013.GIF" wi="677" he="103" />表示坐标为(i<sub>th</sub>,j<sub>th</sub>,k<sub>th</sub>)的网格,<img file="FDA0000746302730000014.GIF" wi="874" he="132" />是(i<sub>th</sub>,j<sub>th</sub>,k<sub>th</sub>)网格的质量,x,y是空间坐标的变量,<img file="FDA0000746302730000015.GIF" wi="105" he="59" />是位于网格(i<sub>th</sub>,j<sub>th</sub>,k<sub>th</sub>)内数据点的数量,<img file="FDA0000746302730000016.GIF" wi="329" he="88" />是对应于网格(i<sub>th</sub>,j<sub>th</sub>,k<sub>th</sub>)的空间坐标值,σ<sub>X</sub>、σ<sub>Y</sub>是空间影响因子,用于控制对象间的相互作用力,σ=ch=c(h<sub>1</sub>,h<sub>2</sub>,h<sub>3</sub>)<sup>T</sup>,σ<sub>t</sub>=ch<sub>t</sub>,t=1、2、3,c是比例常数,h=(h<sub>1</sub>,h<sub>2</sub>,h<sub>3</sub>)<sup>T</sup>是核密度估计的窗宽,K(x)为单位势函数,满足:∫K(x)dx=1,∫xK(x)dx=0,0<R(K)=∫K(x)<sup>2</sup>dx<∞;步骤2、采用下山法对小网格聚类,根据小网格的势值分布,基于分层网格结构,聚类过程是从任意一个顶点出发,按照梯度值增加的方向不断搜索,顶点以及被搜索的小网格归为一个聚类,直到梯度不再增加为止,将沿路搜索到的小网格划分到聚类中,将聚类结果映射到图像,将图像划分为不同的彼此不相交的区域,具体步骤为:步骤2.1.对步骤1.3中的公式求偏导,得到公式:<img file="FDA0000746302730000021.GIF" wi="1712" he="743" />,利用其计算第一层网格中每个小网格的偏导,并以此来确定所有的顶点网格,通过六邻域模式组合顶点网格来描绘聚类{C<sub>k</sub>},k=1,2,...,v,v是图像的顶点,其中C<sub>k</sub>至少包含一个顶点网格;步骤2.2对于每一个k=1,2,...,v,聚类C<sub>k</sub>中的顶点网格作为初始数据点,沿着梯度值上升的方向搜索网格,直到梯度值不再上升即<img file="FDA0000746302730000022.GIF" wi="215" he="53" />为止,将沿路搜索到的小网格划分到聚类C<sub>k</sub>中;步骤2.3搜索完毕后,对于每一个聚类C<sub>k</sub>,k=1,2,...,v,将每一聚类内所有的数据点映射到图像上,并合并在空间上数据点个数少于M个点的图像碎块,20≤M≤100;一幅图像被划分为R块不重合的初始区域Ω<sub>i</sub>,i=1,2,...,R;步骤3、运用基于区域的Ncut算法合并过分分割的区域;其中所用到的计算权重矩阵W的公式为:<img file="FDA0000746302730000023.GIF" wi="1643" he="282" />n<sub>i</sub>、n<sub>j</sub>分别是Ω<sub>i</sub>、Ω<sub>j</sub>包含的特征数据点的数量,i,j=1,2,...,R且i≠j,F<sub>i</sub>、F<sub>j</sub>分别对应于Ω<sub>i</sub>、Ω<sub>j</sub>中数据点的平均值,i,j=1,2,...,R且i≠j,σ<sub>I</sub>是固定的影响因子。
地址 430072 湖北省武汉市武昌区珞珈山武汉大学