主权项 |
一种基于多特征的灰度不均图像快速分割方法,其特征在于,包括下述步骤:步骤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>φ</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>∈</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>π</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>ϵ</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)接着进行曲线演化;若满足条件,则算法终止,提取零水平集轮廓得到分割结果。 |