发明名称 基于局部引导核拟合能量模型的水平集图像分割方法
摘要 一种基于局部引导核拟合能量模型的水平集图像分割方法,主要包括水平集函数的定义、分割模型能量泛函的构建、能量泛函的形式简化、水平集函数的演化和水平集函数的光滑处理。本发明通过将基于局部引导核拟合能量泛函提取图像局部信息,并在每次迭代过程中进行高斯滤波,以避免周期性重新初始化水平集函数。本发明所提出的分割方法不仅提高了灰度不均匀场景中弱边界目标的分割精度,而且完全避免了周期性的重新初始化问题,同时还降低了算法的计算复杂度。
申请公布号 CN104867143A 申请公布日期 2015.08.26
申请号 CN201510249115.9 申请日期 2015.05.15
申请人 东华理工大学 发明人 方江雄;刘花香;刘军;曾正军;饶利民
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 江西省专利事务所 36100 代理人 胡里程
主权项 一种基于局部引导核拟合能量模型的水平集图像分割方法,其分割模型的能量泛函表达式如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><mi>E</mi><mi>LKF</mi></msup><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><msub><mo>&Integral;</mo><mi>&Omega;</mi></msub><msup><mrow><mo>|</mo><mo>|</mo><mi>&theta;</mi><mrow><mo>(</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mi>&theta;</mi><mrow><mo>(</mo><msup><mi>I</mi><mi>LKF</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0000717764630000011.GIF" wi="770" he="127" /></maths>θ(·)表示从观察空间到高维特征空间的一个映射,I(x)和I<sup>LKF</sup>(x)分别表示原始图像灰度与局部引导核拟合图像灰度,其中I<sup>LKF</sup>=u<sub>1</sub>H(φ)+u<sub>2</sub>(1‑H(φ)),u<sub>1</sub>和u<sub>2</sub>分别表示曲线C内部R<sub>C(s)</sub>和外部区域<img file="FDA0000717764630000013.GIF" wi="99" he="74" />的拟合函数,u<sub>1</sub>和u<sub>2</sub>其表达式如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>u</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>K</mi><mi>&alpha;</mi></msub><mo>*</mo><mo>[</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mi>H</mi><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>]</mo></mrow><mrow><msub><mi>K</mi><mi>&alpha;</mi></msub><mo>*</mo><msub><mi>H</mi><mi>&epsiv;</mi></msub><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><msub><mi>u</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>K</mi><mi>&alpha;</mi></msub><mo>*</mo><mo>[</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>H</mi><mi>&epsiv;</mi></msub><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>]</mo></mrow><mrow><msub><mi>K</mi><mi>&alpha;</mi></msub><mo>*</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>H</mi><mi>&epsiv;</mi></msub><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000717764630000012.GIF" wi="1038" he="143" /></maths>其中K<sub>α</sub>是的高斯核函数,方差α是尺度参数;其分割方法的具体步骤如下:步骤1:输入分割图像,定义水平集函数:用符号距离函数来表示水平集函数φ(x,t),即表示在t时刻点x到轮廓曲线的最短距离;步骤2:分割模型能量泛函的构建:通过考虑图像的局部拟合信息,来建立图像分割模型的能量泛函;步骤3:能量泛函的形式简化:通过将能量泛函中的数据项核函数用原始数据空间的非欧氏距离向量来替换,使二维的数据空间数据转化为一维的数据空间,以便能量泛函能用梯度下降法进行求解;步骤4:水平集函数的演化:采用梯度下降流方法,求取演化曲线的水平集演化方程式;步骤5:水平集函数的光滑处理:在水平集函数每次迭代后对符号距离函数进行高斯滤波操作,使演化曲线能变光滑。
地址 330013 江西省南昌市经开区广兰大道418号
您可能感兴趣的专利