CN104835168B - 基于全局凸优化变分模型的快速多相图像分割方法 - Google Patents
基于全局凸优化变分模型的快速多相图像分割方法 Download PDFInfo
- Publication number
- CN104835168B CN104835168B CN201510250167.8A CN201510250167A CN104835168B CN 104835168 B CN104835168 B CN 104835168B CN 201510250167 A CN201510250167 A CN 201510250167A CN 104835168 B CN104835168 B CN 104835168B
- Authority
- CN
- China
- Prior art keywords
- energy functional
- convex
- energy
- image
- functional
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
Abstract
一种基于全局凸优化变分模型的快速多相图像分割方法,主要包括多标记特征函数的定义、非凸能量泛函的构、建能量泛函的转换、全局凸优化的能量泛函构建和基于能量最小化的对偶方法求解过程。在基于区域竞争模型分割模型和多标记特征函数定义的基础上,通过构造非凸的能量泛函和能量泛函的凸化表示,并用对偶方法计算能量泛函的最小化解。本发明所提出的分割方法既能解决非凸目标泛函的局部极优问题,使得分割结果与初始条件无关,又能极大地提高算法的计算效率。
Description
技术领域
本发明涉及的是一种图像处理技术领域图像分割的方法,具体是一种基于全局凸优化变分模型(Global Convex Variational,GCV)的快速多相图像分割方法。
技术背景
图像分割是图像特征提取和分类的重要环节,图像分割的目的就是将图像中的灰度同质区域分离开,并通过各个同质区域的边界来表达。基于变分偏微分方程的图像分割方法,由于具有将模型的初始估计和图像数据先验知识统一于特征提取过程中,并利用分割过程中获得的先验知识来指导分割过程等优势,因此成为近年来图像分割的研究热点。传统的两相Chan-Vese分割模型无法分割多区域同质目标,其非凸能量泛函使得目标函数陷入局部最优,导致分割结果依赖于演化曲线的初始化,而且还使得一些快速数值计算方法无法应用。基于变分水平集的多相图像分割方法,采用差分格式求解复杂的目标函数导致计算效率非常低,极大地限制了该算法的实际应用价值。因此,如何建立有效分割模型的能量泛函,以消除非凸能量泛函局部极优问题使得分割结果依赖于初始条件的影响,是多相图像分割领域的研究热点之一。
经对现有技术文献的检索发现,基于变分法理论的多相水平集方法最初通过集成多种信息来构建分割模型(Zhao H.K.“A variational level set approach tomultiphase motion”(多相运动的变分水平集方法),Journal of Computational Physics(1996)7:79-195)。多相分段恒定的变分水平集分割模型(Vese L.,Chan T.“A multiphaselevel set framework for image segmentation using the Mumford and Shah model”(基于Mumford-Shah模型的多相水平集图像分割框架),International Journal ofComputer Vision)(2002)50(3):271-293),进一步简化了Mumford-Shah模型,用N个水平集函数表达2N个相,避免多个水平集函数的重叠。方江雄等提出的多分辨率多区域变分水平集图像分割方法(专利号:CN102044077B),通过多分辨率技术执行分割图像的演化曲线来解决初始化水平集能量函数陷入局部能量最小值,降低了噪声的干扰、减小了搜索的空间。但是,这些多相图像分割方法所建立的能量泛函均是非凸函数,导致局部极优问题的产生,而且还使得现有的快速算法无法直接应用。
发明内容
本发明的目的是,通过提出一种基于全局凸优化变分模型的快速多相图像分割方法,来解决非凸目标泛函的局部极优问题,使得分割结果与初始条件无关,并能用快速的数值计算方法求解。
本发明的技术方案:在基于区域竞争模型分割模型和多标记特征函数定义的基础上,通过构造非凸的能量泛函和能量泛函的凸化表示,并用对偶方法计算能量泛函的最小化问题,从而避免能量泛函的局部极优问题。具体步骤如下:
步骤1:多标记特征函数的定义。在不相交的图像子域Ωi中(整个图像域为Ω)的点x,多标记特征函数ui(x)∈[0,1],其中i=1,……,N,其表达式如下:
多标记约束项满足
步骤2:非凸能量泛函的构建。本发明采用了Fang Jiangxiong(方江雄)等在Optical Engineering上“Statistical approaches to automatic level set imagesegmentation with multiple regions”中定义的能量泛函。假设给定的图像子域Ω中点坐标为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等在“SIAM Journal on Applied Mathematics”上提出的“Algorithmsfor finding global minimizers of image segmentation and denoising models”中方法,通过在能量泛函中增加对偶变量(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-regionCompetition Segmentation,MCS)模型分割遥感图像结果;
其中:图3(a)和图3(e)分别显示了MCS模型和GCV模型分割遥感图像的初始轮廓;图3(b)和图3(c)分别显示了MCS模型分割结果的两目标区域;图3(f)和图3(g)分别显示了GCV模型分割结果的两目标区域;图3(d)和图3(h)分别显示了MCS模型和GCV模型分割遥感图像结果。
具体实施方式
在基于区域竞争模型分割模型和多标记特征函数定义的基础上,通过构造非凸的能量泛函和能量泛函的凸化表示,并用对偶方法求解能量泛函的最小解。本发明所提出的分割方法既能解决非凸目标泛函的局部极优问题,使得分割结果与初始条件无关,又能极大地提高分割算法的计算效率。本发明具体实施步骤包括如下:
(1)输入分割图像,设置初始化参数:分割区域数目N,两个权重系数附加θ,最大迭代次数;
(2)各区域参量初始化:计算各区域的像素均值ci,特征函数ui和v(ui),初始化其中i=1,…N-1;
(3)运行如下的递归过程当满足条件或循环次数不大于最大迭代次数,重复如下操作:
①通过固定点方法所求解的公式,计算各向量pi,i=1,…N-1;
②计算各区域像素均值ci,i=1,…N;
③根据公式v(ui)=min{max{ui(x)-θipi(x),0},1},更新参数v(ui),i=1,…N-1。
(4)输出各分割区域的图像和分割结果。
图2显示了在不同的初始位置条件下,基于全局凸优化变分模型的快速算法分割医学图像结果。实验中,分割区域数为3,试验中区域标记参数λi=1,i=1,…,N。图中第一行分别显示了初始化轮廓形状分别为两个圆形、矩形、直线和两个三角形。第二行和第三行分别显示了基于GCV模型的多相分割方法分割同一幅医学图像后,演化曲线最终停止位置和分割结果。
图3比较了本发明所提出的GCV分割模型和多区域竞争分割(MCS)模型分割遥感图像结果。实验中,分割区域数为3,区域标记参数λ0=0.5。在两组图像试验中,图像初始化轮廓位置相同,图(a)和图(e)显示了演化曲线最终停靠位置轮廓。图(b)和图(c)显示MCS模型对应的两目标区域。图(f)和图(g)显示GCV模型对应的两目标区域。图(d)和图(h)显示了MCS模型和GCV模型分割的最终结果。从分割的结果看,基于GCV模型的分割方法比基于MCS模型的分割方法有更好的效果。
Claims (4)
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,采用凸松弛方法来解决非凸问题,其凸集η定义如下:
用图像变量来替代标记函数,其能量泛函表达式可重写为:
该能量泛函由数据项和规则项构成,规则项中用全变分范数表示,函数g(x)=1/(1+x2),能量泛函表示可改写为:
。
3.根据权利要求1所述的基于全局凸优化变分模型的快速多相图像分割方法,其特征在于:全局凸优化的能量泛函构建,对任意ci∈R+,i=1,…,N-1,函数g(x)∈[0,1],步骤3中非凸约束能量泛函转化为凸优化非约束最小化能量泛函如下:
其中κ(ui)=max{min{2ui-1},1},惩罚因子α>1/2||R(x)||L∞。
4.根据权利要求1所述的基于全局凸优化变分模型的快速多相图像分割方法,其特征在于:基于能量最小化的对偶方法求解过程,通过用对偶方法增加辅助变量来解决最小化问题,其表达式如下:
其中,α为大于0的常数,κ(x)=max{min{2|x-1|},1};各子区域Ωi像素均值在迭代过程中更新方程式如下:
依据Chan等在“SIAM Journal on Applied Mathematics”上提出的“Algorithms forfinding global minimizers of image segmentation and denoising models”中方法,通过在能量泛函中增加对偶变量(ui,vi),其能量泛函可改写如下:
其中θ取值为0.01;因为能量泛函是凸函数,最小化能量泛函就能得到全局最小值;为解决ui和vi的凸优化问题,采用两步迭代算法计算最小化问题;
第一步:当vi固定,求解ui,其表达可写成:
假定ui(x)=vi(x)-θidivpi,向量pi满足等式通过固定点方法
其中t为时间变量;我们通过条件来终止固定点迭代问题;
第二步:求解vi(x),与vi(x)能量泛函表达式如下:
其中vi(x)=min{max{ui(x)-θpi(x),0},1}。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510250167.8A CN104835168B (zh) | 2015-05-15 | 2015-05-15 | 基于全局凸优化变分模型的快速多相图像分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510250167.8A CN104835168B (zh) | 2015-05-15 | 2015-05-15 | 基于全局凸优化变分模型的快速多相图像分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104835168A CN104835168A (zh) | 2015-08-12 |
CN104835168B true CN104835168B (zh) | 2017-08-18 |
Family
ID=53813033
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510250167.8A Expired - Fee Related CN104835168B (zh) | 2015-05-15 | 2015-05-15 | 基于全局凸优化变分模型的快速多相图像分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104835168B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106803260B (zh) * | 2016-12-28 | 2019-08-09 | 辽宁师范大学 | 基于局部熵凸优化的红外船舶图像活动轮廓分割方法 |
CN107665494A (zh) * | 2017-10-11 | 2018-02-06 | 青岛大学 | 一种自适应含噪sar图像全变分分割方法 |
CN108898611B (zh) * | 2018-06-19 | 2021-09-24 | 东华理工大学 | 基于显著感知先验的模糊区域活动轮廓分割模型 |
CN110415262B (zh) * | 2019-08-07 | 2022-01-14 | 青岛大学 | 实现三维图像多相分割计算机装置及三维分割方法、设备 |
CN111105430B (zh) * | 2019-11-28 | 2023-08-08 | 青岛大学 | 基于Landmark单纯形约束的变分水平集图像分割方法 |
CN114580212A (zh) * | 2022-04-28 | 2022-06-03 | 支付宝(杭州)信息技术有限公司 | 一种资源的优化配置方法、装置以及设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103544709A (zh) * | 2013-11-07 | 2014-01-29 | 东北大学 | 基于显著活动轮廓模型的硅钢板表面微小缺陷检测方法 |
CN104318569A (zh) * | 2014-10-27 | 2015-01-28 | 北京工业大学 | 基于深度变分模型的空间显著性区域提取方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7706610B2 (en) * | 2005-11-29 | 2010-04-27 | Microsoft Corporation | Segmentation of objects by minimizing global-local variational energy |
-
2015
- 2015-05-15 CN CN201510250167.8A patent/CN104835168B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103544709A (zh) * | 2013-11-07 | 2014-01-29 | 东北大学 | 基于显著活动轮廓模型的硅钢板表面微小缺陷检测方法 |
CN104318569A (zh) * | 2014-10-27 | 2015-01-28 | 北京工业大学 | 基于深度变分模型的空间显著性区域提取方法 |
Non-Patent Citations (3)
Title |
---|
A new continuous max-flow algorithm for multiphase image segmentation using super-level set functions;Jun Liu 等;《Journal of Visual Communication and Image Representation》;20140831;第25卷(第6期);第1472-1488页 * |
Fast and globally convex multiphase active contours for brain MRI segmentation;Juan C.Moreno 等;《Computer Vision & Image Understanding》;20140509;第125卷(第8期);第237-250页 * |
多相图像分割的变分模型及其Split Bbregman迭代算法;王琦;《中国优秀硕士学位论文全文数据库信息科技辑》;20120615;第I138-2070页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104835168A (zh) | 2015-08-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104835168B (zh) | 基于全局凸优化变分模型的快速多相图像分割方法 | |
Xu et al. | Object segmentation using graph cuts based active contours | |
Bae et al. | Global minimization for continuous multiphase partitioning problems using a dual approach | |
Kim et al. | Improved simple linear iterative clustering superpixels | |
Jeon et al. | Unsupervised hierarchical image segmentation with level set and additive operator splitting | |
CN110458172A (zh) | 一种基于区域对比度检测的弱监督图像语义分割方法 | |
CN110189339A (zh) | 深度图辅助的主动轮廓抠图方法及系统 | |
CN101504767B (zh) | 基于水平集接力的图像分割方法 | |
CN102509327B (zh) | 一种图像空洞填补的多尺度全局采样方法 | |
CN108629809B (zh) | 一种精确高效的立体匹配方法 | |
CN104299263A (zh) | 一种基于单幅图像建模云场景的方法 | |
CN104732545A (zh) | 结合稀疏近邻传播和快速谱聚类的纹理图像分割方法 | |
Schuster et al. | FlowFields++: Accurate optical flow correspondences meet robust interpolation | |
Wang et al. | Interactive multilabel image segmentation via robust multilayer graph constraints | |
CN102446357B (zh) | 基于自适应有限元的水平集sar图像分割方法 | |
Yang | Local smoothness enforced cost volume regularization for fast stereo correspondence | |
CN104867143A (zh) | 基于局部引导核拟合能量模型的水平集图像分割方法 | |
CN107909597B (zh) | 一种具有边缘保持的多尺度mrf模型图像分割方法 | |
CN103413332B (zh) | 基于两通道纹理分割主动轮廓模型的图像分割方法 | |
Singh et al. | Automatic brain MRI image segmentation using FCM and LSM | |
Pan et al. | A novel level set approach for image segmentation with landmark constraints | |
CN108898611B (zh) | 基于显著感知先验的模糊区域活动轮廓分割模型 | |
CN104574379B (zh) | 一种基于目标多部件学习的视频分割方法 | |
CN111105430B (zh) | 基于Landmark单纯形约束的变分水平集图像分割方法 | |
CN105374024B (zh) | 高分辨率卫星影像水上桥梁提取的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170818 Termination date: 20180515 |