发明名称 图像分割方法与系统
摘要 本发明公开了一种图像分割方法与系统,先采用基于局部对比度和颜色一致性的显著性分析方法,得到一个初始的前景和背景,再通过迭代分析,得到更精确的前景和背景。整个过程可以由计算机完成,不需人工介入,具有较高的适用性,所得的初始前景相比人工设定的前景,更为准确,迭代之后得到的最终的前景也具有更高的准确性,同时还减少了迭代次数,从而提高了分割效率。
申请公布号 CN103208123B 申请公布日期 2016.03.02
申请号 CN201310139396.3 申请日期 2013.04.19
申请人 广东图图搜网络科技有限公司 发明人 钟海兰
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 广州华进联合专利商标代理有限公司 44224 代理人 王茹;曾旻辉
主权项 一种图像分割方法,其特征在于,包括如下步骤:抽取图像的局部对比度和颜色一致性特征进行显著性分析,得到图像的初始前景区域和初始背景区域;对所述初始前景区域和所述初始背景区域进行迭代分析,获得精确前景区域,所述局部对比度为图像上预定范围的区域与其周边区域颜色直方图的卡方距离;所述颜色一致性为图像上预定范围的区域内颜色的一致性;所述初始前景区域为局部对比度和颜色一致性综合结果高于第一前景阈值,且颜色一致性高于第二前景阈值的区域;所述初始背景区域为局部对比度和颜色一致性综合结果低于第一背景阈值,或局部对比度低于第二背景阈值的区域;其中,对所述初始前景区域和所述初始背景区域采用颜色聚类和最小割的迭代方法进行迭代分析,所述步骤抽取图像的局部对比度和颜色一致性特征进行显著性分析包括:在图像上均匀采样c<sub>1</sub>×c<sub>2</sub>个像素点,按如下公式计算每个像素点x的局部对比度特征LC,并对LC进行归一化:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>L</mi><mi>C</mi><mo>=</mo><munder><mo>&Sigma;</mo><mrow><mo>{</mo><msup><mi>x</mi><mo>&prime;</mo></msup><mo>|</mo><msup><mi>x</mi><mo>&prime;</mo></msup><mo>&Element;</mo><msup><mi>R</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>}</mo></mrow></munder><msub><mi>w</mi><mrow><msup><mi>x</mi><mo>&prime;</mo></msup><mi>x</mi></mrow></msub><mi>d</mi><mrow><mo>(</mo><msup><mi>R</mi><mo>*</mo></msup><mo>(</mo><mi>x</mi><mo>)</mo><mo>,</mo><msubsup><mi>R</mi><mi>s</mi><mo>*</mo></msubsup><mo>(</mo><mi>x</mi><mo>)</mo><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000788852950000011.GIF" wi="900" he="176" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>w</mi><mrow><msup><mi>x</mi><mo>&prime;</mo></msup><mi>x</mi></mrow></msub><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mn>0.5</mn><msubsup><mi>&sigma;</mi><mi>x</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>|</mo><mo>|</mo><mi>x</mi><mo>-</mo><msup><mi>x</mi><mo>&prime;</mo></msup><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000788852950000012.GIF" wi="780" he="95" /></maths>σ<sub>x</sub>=s<sub>σ</sub>·|R<sup>*</sup>(x)|<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><mi>R</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mi>arg</mi><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><mi>R</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></munder><mi>d</mi><mrow><mo>(</mo><mi>R</mi><mo>(</mo><mi>x</mi><mo>)</mo><mo>,</mo><msub><mi>R</mi><mi>s</mi></msub><mo>(</mo><mi>x</mi><mo>)</mo><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000788852950000013.GIF" wi="782" he="103" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>d</mi><mrow><mo>(</mo><mi>R</mi><mo>(</mo><mi>x</mi><mo>)</mo><mo>,</mo><msub><mi>R</mi><mi>s</mi></msub><mo>(</mo><mi>x</mi><mo>)</mo><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>&Sigma;</mi><mi>i</mi></msub><mfrac><msup><mrow><mo>(</mo><msup><mi>R</mi><mi>i</mi></msup><mo>-</mo><msubsup><mi>R</mi><mi>s</mi><mi>i</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mrow><msup><mi>R</mi><mi>i</mi></msup><mo>+</mo><msubsup><mi>R</mi><mi>s</mi><mi>i</mi></msubsup></mrow></mfrac></mrow>]]></math><img file="FDA0000788852950000014.GIF" wi="851" he="183" /></maths>其中,R(x)表示以像素点x为中心的长边为预定长度的区域,R<sub>s</sub>(x)表示R(x)周围预定范围内的区域,<img file="FDA0000788852950000015.GIF" wi="148" he="79" />表示R<sup>*</sup>(x)周围预定范围内的区域,s<sub>σ</sub>为经验值,R<sup>i</sup>表示R(x)对应的颜色直方图的第i个分量,<img file="FDA0000788852950000016.GIF" wi="78" he="83" />表示R<sub>s</sub>(x)对应的颜色直方图的第i个分量;采用K均值聚类法将图像聚成k类,将聚类结果作为初始解,再采用高斯混合模型将图像聚成k类,得到图像中每个像素点y属于聚类c的概率p(c|y),c=1,2,…,k,按如下公式计算图像中每个像素点y的颜色一致性特征CS,并对其中的V(c)和D(c)进行归一化,得到归一化的CS:CS=∑<sub>c</sub>p(c|y)·(1‑V(c))·(1‑D(c))<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><mi>c</mi><mo>|</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>w</mi><mi>c</mi></msub><mi>N</mi><mrow><mo>(</mo><mi>y</mi><mo>|</mo><msub><mi>&mu;</mi><mi>c</mi></msub><mo>,</mo><msub><mi>&sigma;</mi><mi>c</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>&Sigma;</mi><mi>c</mi></msub><msub><mi>w</mi><mi>c</mi></msub><mi>N</mi><mrow><mo>(</mo><mi>y</mi><mo>|</mo><msub><mi>&mu;</mi><mi>c</mi></msub><mo>,</mo><msub><mi>&sigma;</mi><mi>c</mi></msub><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000788852950000021.GIF" wi="658" he="159" /></maths>V(c)=V<sub>h</sub>(c)+V<sub>V</sub>(c)<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>V</mi><mi>h</mi></msub><mrow><mo>(</mo><mi>c</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>&Sigma;</mi><mi>y</mi></msub><mi>p</mi><mrow><mo>(</mo><mi>c</mi><mo>|</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac><msub><mi>&Sigma;</mi><mi>y</mi></msub><mi>p</mi><mrow><mo>(</mo><mi>c</mi><mo>|</mo><mi>y</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mo>|</mo><msub><mi>y</mi><mi>h</mi></msub><mo>-</mo><msub><mi>M</mi><mi>h</mi></msub><mrow><mo>(</mo><mi>c</mi><mo>)</mo></mrow><msup><mo>|</mo><mn>2</mn></msup></mrow>]]></math><img file="FDA0000788852950000022.GIF" wi="1093" he="163" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>M</mi><mi>h</mi></msub><mrow><mo>(</mo><mi>c</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>&Sigma;</mi><mi>y</mi></msub><mi>p</mi><mrow><mo>(</mo><mi>c</mi><mo>|</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac><msub><mi>&Sigma;</mi><mi>y</mi></msub><mi>p</mi><mrow><mo>(</mo><mi>c</mi><mo>|</mo><mi>y</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>y</mi><mi>h</mi></msub></mrow>]]></math><img file="FDA0000788852950000023.GIF" wi="829" he="162" /></maths><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>V</mi><mi>v</mi></msub><mrow><mo>(</mo><mi>c</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>&Sigma;</mi><mi>y</mi></msub><mi>p</mi><mrow><mo>(</mo><mi>c</mi><mo>|</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac><msub><mi>&Sigma;</mi><mi>y</mi></msub><mi>p</mi><mrow><mo>(</mo><mi>c</mi><mo>|</mo><mi>y</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mo>|</mo><msub><mi>y</mi><mi>v</mi></msub><mo>-</mo><msub><mi>M</mi><mi>v</mi></msub><mrow><mo>(</mo><mi>c</mi><mo>)</mo></mrow><msup><mo>|</mo><mn>2</mn></msup></mrow>]]></math><img file="FDA0000788852950000024.GIF" wi="1092" he="167" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>M</mi><mi>v</mi></msub><mrow><mo>(</mo><mi>c</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>&Sigma;</mi><mi>y</mi></msub><mi>p</mi><mrow><mo>(</mo><mi>c</mi><mo>|</mo><mi>y</mi><mo>)</mo></mrow></mrow></mfrac><msub><mi>&Sigma;</mi><mi>y</mi></msub><mi>p</mi><mrow><mo>(</mo><mi>c</mi><mo>|</mo><mi>y</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>y</mi><mi>v</mi></msub></mrow>]]></math><img file="FDA0000788852950000025.GIF" wi="830" he="167" /></maths>D(c)=∑<sub>y</sub>p(c|y)d<sub>y</sub>其中,V(c)表示聚类c的方差,D(c)表示聚类c的距离权重,w<sub>c</sub>、μ<sub>c</sub>、σ<sub>c</sub>分别表示聚类c的权重、均值和方差,N(y|μ<sub>c</sub>,σ<sub>c</sub>)表示y服从均值为μ<sub>c</sub>、方差为σ<sub>c</sub>的高斯分布,y<sub>h</sub>、y<sub>v</sub>分别表示像素点y的横向和纵向坐标,d<sub>y</sub>表示像素点y距离图像中心点的距离;按如下公式确定所述初始背景区域Background和所述初始前景区域Foreground:Background=(LC+CS&lt;f<sub>1B</sub>)|(LC&lt;f<sub>2B</sub>)Foreground=(LC+CS&gt;f<sub>1F</sub>)&amp;(CS&gt;f<sub>2F</sub>)其中,|表示逻辑“或”,&amp;表示逻辑“和”,f<sub>1B</sub>、f<sub>2B</sub>、f<sub>1F</sub>、f<sub>2F</sub>为经验值,分别表示所述第一背景阈值、第二背景阈值、第一前景阈值和第二前景阈值。
地址 528311 广东省佛山市顺德区北滘镇三洪奇居委会三乐路北1号M栋4楼