发明名称 一种基于多特征的灰度不均图像快速分割方法
摘要 本发明公开了一种基于多特征的灰度不均图像快速分割方法,通过引入相似性理论快速估计偏压场信息,从而简化了局部信息模型,不但在运行速度上有了很大的提高,而且降低了对初始化轮廓信息的敏感度;本方法与经典算法相比,通过构造与阶跃函数更近似的Heaviside函数,在分割曲线的光滑度上表现更优;本方法通过引入双重终止条件,根据图像内容的不同自适应的结束曲线演化进程,从而提高了分割算法的速度。
申请公布号 CN105139398A 申请公布日期 2015.12.09
申请号 CN201510531170.7 申请日期 2015.08.26
申请人 武汉大学 发明人 何发智;于海平
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 武汉科皓知识产权代理事务所(特殊普通合伙) 42222 代理人 张火春
主权项 一种基于多特征的灰度不均图像快速分割方法,其特征在于,包括下述步骤:步骤1,输入待分割图像:I0;步骤2,设置初始轮廓C<sub>0</sub>,使用公式(2)初始化水平集函数φ<sub>0</sub>,设置时间步长:Δt=0.1,用来控制曲线光滑度函数Heaviside中的参数设置为:ε=1.5,长度惩罚项参数:μ=λ×255<sup>2</sup>,λ∈(0,1);公式(3)H<sup>new</sup>(x)是曲线光滑度Heaviside函数;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&phi;</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>2</mn></mtd><mtd><mrow><mi>x</mi><mo>&Element;</mo><msub><mi>C</mi><mn>0</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mn>2</mn></mrow></mtd><mtd><mrow><mi>o</mi><mi>t</mi><mi>h</mi><mi>e</mi><mi>r</mi><mi>w</mi><mi>i</mi><mi>s</mi><mi>e</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000789648630000011.GIF" wi="1180" he="178" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mi>H</mi><mrow><mi>n</mi><mi>e</mi><mi>w</mi></mrow></msup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><mfrac><mi>&pi;</mi><mn>2</mn></mfrac><mo>-</mo><mi>a</mi><mi>t</mi><mi>a</mi><mi>n</mi><mo>(</mo><mfrac><mi>x</mi><mi>&epsiv;</mi></mfrac><mo>)</mo></mrow><mo>)</mo><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000789648630000012.GIF" wi="1178" he="128" /></maths>其中,π为圆周率;步骤3,根据步骤2中初始化信息建立局部信息模型、全局信息模型以及正则化能量模型,将三类模型嵌入到水平集框架中,得到总能量信息模型并进行曲线演化;所述的的总能量信息模型表示如下:E(φ)=α·E<sup>L</sup>(φ)+β·E<sup>G</sup>(φ)+R<sub>p</sub>(φ)         (4)其中E(φ)表示总的能量模型;E<sup>L</sup>(φ)表示局部统计信息模型;E<sup>G</sup>(φ)表示全局统计信息模型;R<sub>p</sub>(φ)用来表示距离正则化能量模型;α,β是分别用来控制局部和全局信息模型的非负常数;步骤4,将步骤3中三类信息模型嵌入到水平集框架进行曲线演化,使用双重终止条件进行判断曲线演化是否终止:若条件不满足,则跳转到步骤(3)接着进行曲线演化;若满足条件,则算法终止,提取零水平集轮廓得到分割结果。
地址 430072 湖北省武汉市武昌区珞珈山武汉大学