主权项 |
一种基于局部高斯和玻尔兹曼的水平集图像分割方法,包括如下步骤:(1)输入待分割图像I到局部高斯模型中,得到关于图像I的能量函数E(I);(2)通过Euler‑Lagrange方程得到能量函数E(I)对应的演化方程Y(φ),其中φ为演化曲线;(3)将演化曲线初始化为初始演化曲线φ(X,n),其中X为图像上的点,状态n=0;(4)求解演化曲线为初始演化曲线φ(X,n)的演化方程Y(φ),得到新的演化曲线φ(X,n+1);(5)计算停止函数S,设置阈值T,比较阈值T和停止函数S的大小:(5.1)根据新的演化曲线φ(X,n+1)和初始演化曲线φ(X,n)计算停止函数S:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>S</mi><mo>=</mo><mfrac><mrow><msub><mi>m</mi><mn>1</mn></msub><mo>-</mo><msub><mi>m</mi><mn>3</mn></msub></mrow><msqrt><mrow><msub><mi>σ</mi><mn>1</mn></msub><mo>-</mo><msub><mi>σ</mi><mn>3</mn></msub></mrow></msqrt></mfrac><mo>+</mo><mfrac><mrow><msub><mi>m</mi><mn>2</mn></msub><mo>-</mo><msub><mi>m</mi><mn>4</mn></msub></mrow><msqrt><mrow><msub><mi>σ</mi><mn>2</mn></msub><mo>-</mo><msub><mi>σ</mi><mn>4</mn></msub></mrow></msqrt></mfrac></mrow>]]></math><img file="FDA0000784820740000011.GIF" wi="533" he="137" /></maths>其中,m<sub>1</sub>,m<sub>2</sub>,σ<sub>1</sub>,σ<sub>2</sub>分别为新的演化曲线φ(X,n+1)内外均值和内外方差,m<sub>3</sub>,m<sub>4</sub>,σ<sub>3</sub>,σ<sub>4</sub>分别为初始演化曲线φ(X,n)的内外均值和内外方差;(5.2)设置阈值T,比较阈值T和停止函数S的大小:当T<S时,n=n+1,继续迭代,返回步骤(4);当T≥S时,迭代终止,得到演化曲线φ<sub>r</sub>=φ(X,n+1);(6)对演化曲线φ<sub>r</sub>进行形态学开操作,得到最终分割结果φ<sub>z</sub>。 |