基于全局凸优化变分模型的快速多相图像分割方法技术

一种基于全局凸优化变分模型的快速多相图像分割方法，主要包括多标记特征函数的定义、非凸能量泛函的构、建能量泛函的转换、全局凸优化的能量泛函构建和基于能量最小化的对偶方法求解过程。在基于区域竞争模型分割模型和多标记特征函数定义的基础上，通过构造非凸的能量泛函和能量泛函的凸化表示，并用对偶方法计算能量泛函的最小化解。本发明专利技术所提出的分割方法既能解决非凸目标泛函的局部极优问题，使得分割结果与初始条件无关，又能极大地提高算法的计算效率。

【技术实现步骤摘要】
基于全局凸优化变分模型的快速多相图像分割方法
本专利技术涉及的是一种图像处理
图像分割的方法，具体是一种基于全局凸优化变分模型(GlobalConvexVariational，GCV)的快速多相图像分割方法。技术背景图像分割是图像特征提取和分类的重要环节，图像分割的目的就是将图像中的灰度同质区域分离开，并通过各个同质区域的边界来表达。基于变分偏微分方程的图像分割方法，由于具有将模型的初始估计和图像数据先验知识统一于特征提取过程中，并利用分割过程中获得的先验知识来指导分割过程等优势，因此成为近年来图像分割的研究热点。传统的两相Chan-Vese分割模型无法分割多区域同质目标，其非凸能量泛函使得目标函数陷入局部最优，导致分割结果依赖于演化曲线的初始化，而且还使得一些快速数值计算方法无法应用。基于变分水平集的多相图像分割方法，采用差分格式求解复杂的目标函数导致计算效率非常低，极大地限制了该算法的实际应用价值。因此，如何建立有效分割模型的能量泛函，以消除非凸能量泛函局部极优问题使得分割结果依赖于初始条件的影响，是多相图像分割领域的研究热点之一。经对现有技术文献的检索发现，基于变分法理论的多相水平集方法最初通过集成多种信息来构建分割模型(ZhaoH.K.“Avariationallevelsetapproachtomultiphasemotion”(多相运动的变分水平集方法)，JournalofComputationalPhysics(1996)7:79-195)。多相分段恒定的变分水平集分割模型(VeseL.,ChanT.“AmultiphaselevelsetframeworkforimagesegmentationusingtheMumfordandShahmodel”(基于Mumford-Shah模型的多相水平集图像分割框架)，InternationalJournalofComputerVision)(2002)50(3):271-293)，进一步简化了Mumford-Shah模型，用N个水平集函数表达2N个相，避免多个水平集函数的重叠。方江雄等提出的多分辨率多区域变分水平集图像分割方法(专利号：CN102044077B)，通过多分辨率技术执行分割图像的演化曲线来解决初始化水平集能量函数陷入局部能量最小值，降低了噪声的干扰、减小了搜索的空间。但是，这些多相图像分割方法所建立的能量泛函均是非凸函数，导致局部极优问题的产生，而且还使得现有的快速算法无法直接应用。
技术实现思路
本专利技术的目的是，通过提出一种基于全局凸优化变分模型的快速多相图像分割方法，来解决非凸目标泛函的局部极优问题，使得分割结果与初始条件无关，并能用快速的数值计算方法求解。本专利技术的技术方案：在基于区域竞争模型分割模型和多标记特征函数定义的基础上，通过构造非凸的能量泛函和能量泛函的凸化表示，并用对偶方法计算能量泛函的最小化问题，从而避免能量泛函的局部极优问题。具体步骤如下：步骤1：多标记特征函数的定义。在不相交的图像子域Ωi中(整个图像域为Ω)的点x，多标记特征函数ui(x)∈[0,1]，其中i＝1,……,N，其表达式如下：多标记约束项满足步骤2：非凸能量泛函的构建。本专利技术采用了FangJiangxiong(方江雄)等在OpticalEngineering上“Statisticalapproachestoautomaticlevelsetimagesegmentationwithmultipleregions”中定义的能量泛函。假设给定的图像子域Ω中点坐标为x，用N-1个水平集函数φi表达N个区域(i＝1,…,N-1)，其能量泛函表达式如下：其中r＝(r1,…,rN)，λi(i＝1,2,…,N)是大于0的正数，第i区域像素均值为ci，I(x)为x坐标点的像素值，表示图像梯度，φ＝(φ1,…,φN-1)，H(x)是Heaviside函数，H(φ0)≡0，δ(x)为平滑函数H(x)的导数，λ1＞0,i＝1,…,N。根据欧拉拉格朗日等式最小化能量泛函，其演化方程如下：其中t时间变量，假设Hi＝H(φi(x))，变量Φi(φi)和像素均值ci表达式如下：为了找到全局最小化分割模型，通过去掉函数δ(x)，其分割模型的能量泛函如下：其中c＝(c1,…,cN)。步骤3：能量泛函的转换。为求解步骤2中的非凸能量泛函，多标记特征函数重定义如下：其中λ0＞0。通过定义变量ui,i＝1,…,N-1，本专利技术采用凸松弛方法来解决非凸问题，假设其凸集η定义如下：用图像变量来替代标记函数，其能量泛函表达式可重写为：其中U＝{u1,…,uN-1}。该能量泛函由数据项和规则项构成，规则项中用全变分范数表示，函数g(x)＝1/(1+x2)，能量泛函表示可改写为：步骤4：全局凸优化的能量泛函构建。针对任意ci∈R+,i＝1,…,N-1，函数g(x)∈[0,1]，步骤3中非凸约束能量泛函转化为凸优化非约束最小化能量泛函如下：其中常量α＞0远大于λ2，惩罚因子步骤5：基于能量最小化的对偶方法求解过程。为最小化能量泛函E(U,c)，本专利技术用对偶方法增加辅助变量来解决最小化问题。针对每相凸问题变分问题，通过用图像变量进行参数交替，其表达式如下：各子区域Ωi像素均值在迭代过程中更新方程式如下：其中i＝1,…,N-1。依据Chan等在“SIAMJournalonAppliedMathematics”上提出的“Algorithmsforfindingglobalminimizersofimagesegmentationanddenoisingmodels”中方法，通过在能量泛函中增加对偶变量(ui,vi)，其能量泛函可改写如下：其中θ是非常小值。因为能量泛函是凸函数，最小化能量泛函就能得到全局最小值。为解决ui和vi的凸优化问题，采用两步迭代算法计算最小化问题。第一步：当vi固定，求解ui，其表达可写成：假定ui(x)＝vi(x)-θidivpi，向量pi满足等式通过固定点方法我们通过条件来终止固定点迭代问题。第二步：求解vi(x)，与vi(x)能量泛函表达式如下：其中vi(x)＝min{max{ui(x)-θpi(x),0},1}。本专利技术通过在基于区域竞争模型分割模型基础上，构造非凸的能量泛函和能量泛函的凸化表示，并用对偶方法计算能量泛函的最小化解。本专利技术所提出的方法既能解决非凸目标泛函的局部极优问题，使得分割结果与初始条件无关，又能极大地提高算法的计算效率。附图说明图1表示本专利技术实施例中基于全局凸优化变分模型的快速多相图像分割方法流程图；图2为在不同的初始位置条件下，基于全局凸优化变分模型的快速方法分割医学图像结果。其中：第一行表示分割图像的初始轮廓图；第二行和第三行分别表示在不同初始条件下演化曲线的最终停止位置和分割结果；图3比较了本专利技术所提出的GCV分割模型和多区域竞争分割(Multi-regionCompetitionSegmentation，MCS)模型分割遥感图像结果；其中：图3(a)和图3(e)分别显示了MCS模型和GCV模型分割遥感图像的初始轮廓；图3(b)和图3(c)分别显示了MCS模型分割结果的两目标区域；图3(f)和图3(g)分别显示了GCV模型分割结果本文档来自技高网...

【技术保护点】
一种基于全局凸优化变分模型的快速多相图像分割方法，其能量泛函构建过程如下：假设给定的图像子域Ω中点坐标为x，用N‑1个水平集函数φi表达N个区域(i＝1,…,N‑1)，其能量泛函表达式如下：min(r,φ)E(r,φ)=min(c,φ)Σi=1N-1∫Ω|▿H(φi(x))|dx+Σi=1N-1λi∫Ω(I(x)-ci)2H(φi(x))Πk=0k=i-1(1-H(φi(x)))dx+λN∫Ω(I(x)-cN)2Πk=1N-1(1-H(φi(x)))dx]]>其中r＝(r1,…,rN),I为图像的灰度均值，表示图像梯度，φ＝(φ1,…,φN‑1)，H(x)是Heaviside函数，H(φ0)≡0，δ(x)为平滑函数H(x)的导数，λ1＞0,i＝1,…,N，其分割方法的具体步骤如下：步骤1：输入分割图像，定义多标记特征函数；在不相交的图像子域Ωi中，给每个区域定义一个多标记特征函数ui(x)∈[0,1](i＝1,……,N‑1)；步骤2：分割模型能量泛函的构建，根据图像概率密度函数分布规律，用N‑1个水平集函数φi表达N个区域(i＝1,…,N‑1)；步骤3：能量泛函的转换，重新定义多标记特征函数，将规则项中用全变分范数TVg(ui)=∫Ωg(|▿I(x)|)|▿ui|dx]]>表示，来简化能量泛函形式；步骤4：全局凸优化能量泛函的构建，将非凸约束能量泛函转化为凸优化非约束最小化能量泛函；步骤5：基于能量最小化的对偶方法求解过程，通过增加辅助变量vi(x)，用对偶方法来解决能量泛函的最小化问题。...

【技术特征摘要】
1.一种基于全局凸优化变分模型的快速多相图像分割方法，其能量泛函构建过程如下：假设给定的图像子域Ω中点坐标为x，用N-1个水平集函数φi表达N个区域，其中i＝1，2，…，N-1，其能量泛函表达式如下：其中r＝(r1,…,rN)，λj是正数，j＝1,2,…,N，ri＝r(ci,Φi,x)＝(λi(I-ci)2-Φi(φi))，第i区域像素均值为ci，I(x)为x坐标点的像素值，表示图像梯度，φ＝(φ1,…,φN-1)，H(x)是Heaviside函数，H(φ0)≡0，δ(x)为平滑函数H(x)的导数，U＝{u1,…,uN-1}，Φi(φi)定义如下：其分割方法的具体步骤如下：步骤1：输入分割图像，定义多标记特征函数；在不相交的图像子域Ωi中，给每个区域定义一个多标记特征函数ui(x)∈[0,1]，其中i＝1,……,N-1；步骤2：分割模型能量泛函的构建；根据图像概率密度函数分布规律，用N-1个水平集函数φi表达N个区域，其中i＝1，…，N-1；步骤3：能量泛函的转换；重新定义多标记特征函数，将规则项中用全变分范数表示，来简化能量泛函形式，其中g(x)＝1/(1+x2)；步骤4：全局凸优化能量泛函的构建，将非凸约束能量泛函转化为凸优化非约束最小化能量泛函；步骤5：基于能量最小化的对偶方法求解过程；通过增加辅助变量vi(x)，用对偶方法来解决能量泛函的最小化问题。2.根据权利要求1所述的基于全局凸优化变分模型的快速多相图像分割方法，其特征在于：能量泛函的转换，多标记特征函数ui(x)重定义如下：其中λ0＞0；通过已定义的变量ui，i＝1，…，N-1，采用凸松弛方法来解决非凸问题，其凸集η定义如下：用图像变量来替代标记函数，其能量泛函表达式可重写为：该能量泛...

【专利技术属性】
技术研发人员：刘花香，方江雄，肖静，刘彬，郭文，
申请(专利权)人：东华理工大学，
类型：发明
国别省市：江西;36