基于混合区域的活动轮廓模型图像分割方法技术

本发明专利技术公开了一种基于混合区域的活动轮廓模型图像分割算法，该方法包括：通过图像增强算法对原始图像进行增强处理；同时使用原始图像的局部灰度拟合项和增强图像的全局灰度拟合项，构建基于混合区域的活动轮廓模型的能量泛函；通过变分法求解构建能量泛函，使用变分水平集的求解框架进行数学计算，即轮廓曲线以隐式水平集的形式表示，然后通过梯度下降流和欧拉一拉格朗日方程将构建的能量泛函转化为一个偏微分方程，并通过迭代逼近的方式求出所述偏微分方程的最优解，进而获得最终的目标边界轮廓。本发明专利技术对初始轮廓曲线的位置以及图像的灰度不均匀性不敏感，分割效率和精度都比较高，且受到噪声的影响较小，得到的分割结果鲁棒性较强。

【技术实现步骤摘要】
基于混合区域的活动轮廓模型图像分割方法
本专利技术涉及一种图像处理领域的图像分割方法，具体涉及一种基于混合区域的活动轮廓模型图像分割方法。技术背景人类大脑从外界获取的各种信息中，绝大部分信息是通过眼睛获得的，物体或者场景发出的光或者反射的光通过眼睛中的晶状体折射到视网膜上，在视网膜上形成其倒立的图像，再由视觉感知神经传达给大脑，经大脑处理后即可看到该物体或场景的图像，而图像中包含其所表达物体或场景丰富的描述信息，因此图像信息是人类日常生活中获取信息的主要途径。随着科学技术的飞速发展，特别是成像技术、计算机技术和信号处理理论的发展，人们将从各种成像器件中获得的图像信息，转变成数字信号，并在计算机中对其进行各种形式的加工和处理，如降噪、压缩、分割、增强、复原、特征抽取、分析和理解等，这些加工和处理即为人们常说的数字图像处理过程。图像分割是指按照图像中各个区域特有的性质，把图像分成几个区域，然后再把感兴趣的区域从这些区域中提取出来的技术。在数字图像处理的过程中，许多加工和处理都是建立在好的图像分割的基础上，只有先从图像中准确分割出目标物体，才能更好得进行下一步的图像处理及分析的工作。因此，图像分割是数字图像处理中一个非常基本、非常关键的技术。虽然图像分割技术很重要，而且到目前为止相关的科研人员也提出了种类多样的图像分割算法，但是由于成像系统的多样性和复杂性，造成了所成图像的多样性和复杂性，所以即使到目前为止科研人员已经提出了许多优秀的图像分割方法，但是还没有一个通用的图像分割算法，各种分割方法也只能用于针对某一类图像的分割问题，也因为这个原因，使得图像分割技术一直是数字图像处理技术当中研究的重点和难点。传统的图像分割方法包括阈值法、边缘检测法以及基于区域提取的方法。而随着待解决问题复杂程度的提高，需要越来越多的高层知识融入到图像分割当中。此时出现了基于活动轮廓模型的图像分割方法，该方法依据轮廓曲线的表达形式的不同而分为参数活动轮廓模型和几何活动轮廓模型。而几何活动轮廓模型具有能有效处理曲线拓扑结构变化的优点，使得研究人员对其的研究越来越多。几何活动轮廓模型中，利用图像全局灰度信息的模型会对图像灰度不均匀性敏感；利用图像局部灰度信息的模型会对初始轮廓位置敏感；而结合图像全局和局部灰度信息的模型会存在一个全局项权重系数设置的问题。
技术实现思路
本专利技术的目的在于提供一种基于混合区域的活动轮廓模型图像分割方法。实现本专利技术目的的技术方案为：一种基于混合区域的活动轮廓模型图像分割方法，包括以下步骤：第一步，原始图像的增强处理：通过图像增强算法对原始图像进行增强处理；第二步，构建能量泛函：同时使用原始图像的局部灰度拟合项和增强图像的全局灰度拟合项，构建基于混合区域的活动轮廓模型的能量泛函；第三步，通过变分法求解构建能量泛函，使用变分水平集的求解框架进行数学计算，即轮廓曲线以隐式水平集的形式表示，然后通过梯度下降流和欧拉一拉格朗日方程将构建的能量泛函转化为一个偏微分方程，并通过迭代逼近的方式求出所述偏微分方程的最优解，进而获得最终的目标边界轮廓。与现有技术相比，本专利技术的显著优点为：(1)本专利技术对初始轮廓曲线的位置以及图像的灰度不均匀性不敏感，分割效率和精度都比较高，且受到噪声的影响较小，得到的分割结果鲁棒性较强；(2)本专利技术能够克服现有技术的不足，将图像中感兴趣的目标从背景中清除、完整地分割出来，并且降低了算法的复杂程度。附图说明图1是本专利技术基于混合区域的活动轮廓模型的图像分割方法的流程图。图2(a)和图2(b)分别是本专利技术仿真实验中原始图像和增强图像示意图。图3(a)～图3(c)是本专利技术仿真实验中选取不同初始轮廓位置时的分割结果图。图4(a)～图4(c)是本专利技术仿真试验中对待分割图像加入不同噪声后的分割结果图。图5(a)～图5(c)是仿真实验中CV模型的分割结果示意图。图6(a)～图6(c)是仿真实验中LRCV模型的分割结果示意图。图7(a)～图7(c)是本专利技术仿真实验中分割结果示意图。图8(a)～图8(c)是本专利技术仿真实验中与手工分割结果的对比图。具体实施方式结合图1，本专利技术的一种基于混合区域的活动轮廓模型图像分割方法，包括以下步骤：第一步，原始图像的增强处理：通过图像增强算法对原始图像进行增强处理，即使用各项同性非抽取小波变换算法改善原始图像中存在的灰度不均匀现象，同时凸显图像中的目标边界；第二步，构建能量泛函：同时使用原始图像的局部灰度拟合项和增强图像的全局灰度拟合项，构建基于混合区域的活动轮廓模型的能量泛函；第三步：通过变分法求解构建能量泛函，使用变分水平集的求解框架进行数学计算，即轮廓曲线以隐式水平集的形式表示，然后通过梯度下降流和欧拉一拉格朗日方程将构建的能量泛函转化为一个偏微分方程，并通过迭代逼近的方式求出所述偏微分方程的最优解，进而获得最终的目标边界轮廓。进一步的，所述的图像增强算法包括以下步骤：S11，通过低通滤波器处理原始图像I0，得到一个尺度系数aj，其中a0＝I0；S12，将相邻两个尺度系数aj和aj-1相减获得小波系数wj；aj＝hj↑*aj-1，hj↑为低通滤波器；S13，选择小波系数进行相加得到增强图像Ie＝w1+w2+w3。进一步的，第二步中构建的能量泛函为：其中，λ1、λ2、η1和η2皆为非负常数，ci(x)为图像局部灰度拟合项，是与图像大小一致的矢量，用于近似表示轮廓曲线内外局部区域的灰度；ci为图像全局灰度拟合项，为轮廓曲线内外的灰度均值，R为水平集正则项，其定义如下：其中μ和ν分别为非负常数，为梯度运算符；R的前半部分为距离惩罚项，可以保证φ在演化过程中与符号距离函数的一致性；后半部分为长度约束项，可以保证轮廓曲线C的光滑性。由于轮廓曲线难以显示表示，因此在实际运算中常用水平集函数φ表示轮廓曲线，从而得到：其中，M1＝H(φ)和M2＝1-H(φ)分别用于描述曲线C的内部和外部区域，H(φ)为正则化Heaviside函数，ε为大于零的常数。进一步的，根据变分法将能量泛函转化为偏微分方程的形式为：F1＝-λ1(I0(x)-c1(x))2+λ2(I0(x)-c2(x))2F2＝-η1(Ie(x)-c1)2+η2(Ie(x)-c2)2其中，div(·)表示曲率运算符，为Dirac函数；根据变分法，得到常数c1和c2满足如下欧拉-拉格朗日方程：∫(Ie(x)-ci)Midx＝0经求解上式得到：根据LRCV(Localregion-basedChan-Vese)模型中局部灰度拟合项的求解得到：式中，Kσ(x-y)为高斯核函数，x为图像中的某一点，y为x邻域范围内的某一点。下面将结合附图，对本专利技术的实施例进行详细描述。实施例如图1所示，本专利技术的基于混合区域的活动轮廓模型的图像分割方法，通过以下实施步骤实现：步骤1，对输入的待分割图像I0进行增强处理：首先对图像I0进行低通滤波，滤波器为上采样滤波算子hj↑，是通过对h0中相邻系数间插入2j-1个零值获得，而h0＝[1,4,6,4,1]/16则由三次B样条曲线衍生获得。经低通滤波获得各级尺度系数cj，然后将相邻的两个尺度系数cj和cj-1相减可计算获得小波系数wj，再将1-3级小波系数相加得到增强图像，即Ie＝w1+w2+w3。步骤2，本文档来自技高网...

【技术保护点】
一种基于混合区域的活动轮廓模型图像分割方法，其特征在于，包括以下步骤：第一步，原始图像的增强处理：通过图像增强算法对原始图像进行增强处理；第二步，构建能量泛函：同时使用原始图像的局部灰度拟合项和增强图像的全局灰度拟合项，构建基于混合区域的活动轮廓模型的能量泛函；第三步，通过变分法求解构建能量泛函，使用变分水平集的求解框架进行数学计算，即轮廓曲线以隐式水平集的形式表示，然后通过梯度下降流和欧拉一拉格朗日方程将构建的能量泛函转化为一个偏微分方程，并通过迭代逼近的方式求出所述偏微分方程的最优解，进而获得最终的目标边界轮廓。

【技术特征摘要】
1.一种基于混合区域的活动轮廓模型图像分割方法，其特征在于，包括以下步骤：第一步，原始图像的增强处理：通过图像增强算法对原始图像进行增强处理；第二步，构建能量泛函：同时使用原始图像的局部灰度拟合项和增强图像的全局灰度拟合项，构建基于混合区域的活动轮廓模型的能量泛函；第三步，通过变分法求解构建能量泛函，使用变分水平集的求解框架进行数学计算，即轮廓曲线以隐式水平集的形式表示，然后通过梯度下降流和欧拉一拉格朗日方程将构建的能量泛函转化为一个偏微分方程，并通过迭代逼近的方式求出所述偏微分方程的最优解，进而获得最终的目标边界轮廓。2.根据权利要求1所述的基于混合区域的活动轮廓模型图像分割方法，其特征在于，所述的图像增强算法包括以下步骤：S11，通过低通滤波器处理原始图像I0，得到一个尺度系数aj，其中a0＝I0；S12，将相邻两个尺度系数aj和aj-1相减获得小波系数wj，aj＝hj↑*aj-1，hj↑为低通滤波器；S13，选择小波系数进行相加得到增强图像Ie＝w1+w2+w3。3.根据权利要求1所述的基于混合区域的活动轮廓模型图像分割方法，其特征在于，第二步中构建的能量泛函为：其中，λ1、λ2、η1和η2皆为非负常数，ci(x)为图像局部灰度拟合项，是与图像大小一致的矢量，用于近似表示轮廓曲线内外局部区域的灰度；ci为图像全局灰度拟合项，为轮廓曲线内外的灰度均值，R为水平集正则项，其定义如下：

【专利技术属性】
技术研发人员：胡茂海，李诗宇，
申请(专利权)人：南京理工大学，
类型：发明
国别省市：江苏,32