全局和局部信息自适应调整的图像分割活动轮廓方法技术

技术编号:13230438 阅读:45 留言:0更新日期:2016-05-13 12:43
本发明专利技术提供一种基于全局和局部信息自适应调整的图像分割活动轮廓方法,(1)定义了一个新颖的自适应平衡函数,其能够根据图像自身特性自动调整各部分的权重,进而驱动曲线演化;(2)在权重函数中,加入了高斯滤过程去正则化水平集函数,同时增加了一项下降因子,加快了曲线的演化速度;(3)惩罚项的引入保证了模型的精确计算和平稳演化;本发明专利技术无论在分割精度还是处理速度上都达到了较好的分割效果,可实现对对灰度分布不均匀的异质图像的分割。

【技术实现步骤摘要】

本专利技术涉及图像处理领域,尤其是一种计算速度快、对初始轮廓曲线位置及图像 噪声具有鲁棒性且能够分割灰度分布不均匀的异质图像的基于全局和局部信息自适应调 整的图像分割活动轮廓方法。
技术介绍
目前,图像分割已在日常生活、医疗技术、军事等领域都得到了广泛的应用,例如 目标提取、边缘检测、目标跟踪等。基于活动轮廓模型的图像分割是该领域中一次重要的理 论方法革新,其基本思想是利用图像的几何特性建立一个能量泛函,在变分法下求能量函 数极小值,得到相应的Euler-Lagrange方程,然后利用泛函分析和数值分析等领域的相关 知识对模型的合理性进行分析,最终提取出感兴趣的图像区域。 现有的活动轮廓模型进行图像分割的模型主要包括以下几类:基于边缘的分割模 型、基于区域的分割模型、基于边缘和区域信息相结合的分割模型等。其中,测地活动轮廓 模型利用图像自身的特性在黎曼空间中寻找能够刻画其特性的测地线,从而实现对目标物 体的分割。水平集函数的引用使模型能够很好地适应演化曲线的拓扑结构变化,其实质就 是将轮廓曲线演化过程转化为偏微分方程求解问题。与参数活动轮廓相比,该模型的数值 计算稳定性得到了增强,曲线的几何特性更容易被捕捉。同时,该模型中还可以添加某些附 加?目息,使其鲁棒性得到提尚,但边缘泄露问题会时常发生;C-V模型是由Chan等人在 Mumford-Shan能量泛函基础上提出的,它假定待分割图像只包含目标和背景两个区域,且 两区域的平均灰度值差别较大。因此,C-V模型只要根据目标和背景区域的平均灰度值的差 别即可完成分割,而不需要利用图像的梯度信息。C-V模型的出现具有重大意义,不仅解决 了 GAC模型中存在的大量缺陷问题,而且能够较好地处理边缘噪声严重以及内部区域出现 遮挡的图像。但是,C-V模型的适用范围亦被严重限制,而且演化速度慢,对异质图像的处理 不是很理想。为此,研究人员提出了很多改进模型,其中,基于边缘信息和区域信息相融合 的方案成为了主要的研究趋势和研究方向,作为该思想的一个主要代表,SBGFRLS模型将C-V模型的全局信息引入到GAC模型中,不仅利用了图像的边缘信息,也结合了图像的区域信 息,进而解决了 GAC模型容易陷入能量函数极小值的问题。但是,该模型依然没能解决GAC和 C-V模型不能分割灰度分布不均匀的异质图像的问题。
技术实现思路
本专利技术是为了解决现有技术所存在的上述技术问题,提供一种计算速度快、对初 始轮廓曲线位置及图像噪声具有鲁棒性且能够分割灰度分布不均匀的异质图像的基于全 局和局部信息自适应调整的图像分割活动轮廓方法。 本专利技术的技术解决方案是:一种基于全局和局部信息自适应调整的图像分割活动 轮廓方法,其特征在于:定义能量泛函形式如下,其中,离项为全局信息项;·项为局部信息项;·项是惩罚项;^Ιν?,???ε: |〇|是关于 图像:?的自适应平衡函数,|7|是图像梯度模值,衡是常数,用以自适应平衡图像的全局信 息和局部信息;#是正参数,其取值范围是%和€是两个常数,分别近似等于轮廓 内部和轮廓外部的平均灰度值,分别为在龙点附近的逼近曲线内外的图像 局部灰度值的光滑函数;模型中自适应平衡函数的选取如下:其中|穿丨>丨¥爾_:,%表示方差为逐的二维高斯滤波器,用来平滑初始图像; 进一步由Euler-Lagrange方程,可得到模型的水平演化方程: i为用来控制演化曲线的平稳演化 和精确计算的惩罚项;利用平方差公式对模型进行因式分解,根据改进CV模型的思想,对于 模型进行变形,得到新的水平集演化方程:步骤1.初始化水平集函数權:??錄; 步骤2.对初始图像进行高斯滤波,利用公式(3)计算权重函数步骤3 ·利用公式(1)、(2)分别计算%,%,: 4,爲;步骤5.利用有限差分法,根据公式(5)更新水平集函数:其中:r:为迭代步长,令:愁I; 步骤6.使用"停止准则"检查演化曲线是否稳定收敛,若稳定收敛,则停止迭代;否则, 转入步骤2。 与现有技术相比,本专利技术具有以下优点:第一,引入的自适应平衡函数可根据曲线 演化的当前状态自适应调整权重信息控制演化的总体演化趋势。而且,全局信息保证模型 的对于噪声和弱边缘图像的处理效果,细节信息保证了模型对复杂背景区域的目标分割精 度。第二,模型通过引入惩罚项,确保了模型演化的平稳性,提高了演化的速度。第三,模型 在演化过程中通过使用高斯卷积来规范水平集函数,避免了传统活动轮廓模型由于重新初 始化和正则化而引起的计算复杂,通过对多种类型图像分割的仿真实验表明,本专利技术具有 分割精度高、速度快和对初始轮廓曲线位置及图像噪声具有鲁棒性的特点,可实现对灰度 分布不均匀的异质图像的分割。【附图说明】 图1是当爾时,自适应平衡函数斤?的曲线图。 图2是本专利技术方法与C-V方法、LBF方法、自适应拟合方法对Ri ce图像的分割结果对 比。 图3是本专利技术方法与C-V方法、LBF方法、自适应拟合方法对血管图像的分割结果对 比。 图4是本专利技术方法与C-V方法、LBF方法、自适应拟合方法对脑部核磁共振图像的分 割结果对比。 图5是本专利技术方法与C-ν方法、LBF方法、自适应拟合方法对脑部CT图像的分割结果 对比。【具体实施方式】 -种基于,其特征在于:定 义能量泛函·%形式如下,其中,焉项为全局信息项,该项使模型在灰度均匀区域快速演化;笔:项为局部信息项, 该项使模型在具有细节信息的区域进行更加细腻的演化;项是惩罚项,用以纠正水平集 函数的误差;壞1矣,1是关于图像?的自适应平衡函数,Μ是图像梯度模值,是常 数,用以自适应平衡图像的全局信息和局部信息;1是正参数,其取值范围是|5ι|; _和 ?2是两个常数,分别近似等于轮廓内部和轮廓外部的平均灰度值,$1?)分别为 在《点附近的逼近曲线内外的图像局部灰度值的光滑函数;其中|表示方差为逐的二维高斯滤波器,用来平滑初始图像;对名 的选取,如果邊过大,可能会发生边缘泄露现象,使边缘检测不够精确;反之,如果泛过小, 将会使模型对噪声过于敏感。:为是下降因子参数,用以控制函数的下降速度;,用以控制演化曲线的平稳演化和精确计算。 进一步由Euler-Lagrange方程,可得到模型的水平演化方程:其中,厶为Laplacian算子;利用平方差公式对模型进行因式分解,根据改进C-V模型 的思想,对于模型进行变形,得到新的水平集演化方程:步骤1 ·初始化水平集函数; 步骤2.对初始图像进行高斯滤波,利用公式(3)计算权重函数歹|;抱|; 当-时,自适应平衡函数的曲线图如图1所示。 步骤3.利用公式(1)、(2)分别计算太,殺;步骤5.利用有限差分法,根据公式(5)更新水平集函数:其中T为迭代步长,说为常数,为了实现水平集函数的正则化并加速演化过程:令步骤6.使用"停止准则"(即演化曲线趋近目标边界时,水平集函数值逐渐变小,当到达 目标边界时,函数值达到最小值,停止演化)。检查演化曲线是否稳定收敛,若稳定收敛,则 停止迭代;否则,转入步骤2。 将本专利技术实施例与其它方法的分割效果进行对比,结果如图2-5所示,可以看出, 本专利技术方法对初始轮廓曲线位置及图像噪声具有更高的鲁棒性,且分割精度更高、速度更 快,可实现对灰度分布不均本文档来自技高网
...

【技术保护点】
一种基于全局和局部信息自适应调整的图像分割活动轮廓方法,其特征在于:定义能量泛函形式如下,其中,项为全局信息项;项为局部信息项;项是惩罚项;是关于图像的自适应平衡函数,是图像梯度模值,是常数,用以自适应平衡图像的全局信息和局部信息;是正参数,其取值范围是;和是两个常数,分别近似等于轮廓内部和轮廓外部的平均灰度值,分别为在x点附近的逼近曲线内外的图像局部灰度值的光滑函数;、及的表达式分别如下:模型中自适应平衡函数的选取如下:(3)其中,表示方差为的二维高斯滤波器;进一步由Euler‑Lagrange方程,可得到模型的水平演化方程:(4)其中,为Laplacian算子,为用来控制演化曲线的平稳演化和精确计算的惩罚项;利用平方差公式对模型进行因式分解,根据改进C‑V模型的思想,对于模型进行变形,得到新的水平集演化方程:(5)步骤1.初始化水平集函数;步骤2.对初始图像进行高斯滤波,利用公式(3)计算权重函数;步骤3.利用公式(1)、(2)分别计算;步骤4.计算惩罚项的值;步骤5.利用有限差分法,根据公式(5)更新水平集函数:其中为迭代步长,令;步骤6.使用“停止准则”检查演化曲线是否稳定收敛,若稳定收敛,则停止迭代;否则,转入步骤2。...

【技术特征摘要】

【专利技术属性】
技术研发人员:王相海方玲玲张冲宋若曦
申请(专利权)人:辽宁师范大学
类型:发明
国别省市:辽宁;21

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1