Graph Cuts算法实现过程中关键问题研究.docVIP

Graph Cuts算法实现过程中关键问题研究 　　摘 要：Graph Cuts（GC）算法是一个基于图论的交互式目标提取算法，该方法将图像建模为一个区域和边缘的约束模型，通过求解该模型的最小割获得一个优化的目标分割边界。在使用程序实现该算法时会遇到种种问题，详细叙述了GC算法的实现过程，为进一步研究GC相关方法奠定基础。 　　关键词：目标提取;Graph Cuts;优化分割 　　DOI：10.11907/rjdk.143584 　　中图分类号：TP303 　　文献标识码：A 文章编号：1672-7800（2014）012-0044-02 　　作者简介：张晓峰（1979-），男，江苏南通人，博士，南通大学计算机科学与技术学院讲师，研究方向为数字图像处理。 　　0 引言 　　 Graph Cuts（GC）lt;supgt;[1]lt;/supgt;是一种广泛应用的的优化算法，在目标提取领域中该算法通过简单地交互标记出部分前景和背景区域，该方法就可以将目标从复杂的背景中分离出来。由于人工交互的存在，该类算法获得了比较可靠的先验知识，随后的自动分割又使得整个过程不需要太多的人工干预就能获得较好的分割结果，因此算法获得了广泛的应用并吸引了很多学者研究并改进[2-4]。 　　 GC算法基于图论，将图像分割问题与图的最小割（min cuts）联系起来，在给定的图像中求解一个最优分割位置。从图像到图论，再到最优分割，各种复杂的概念夹杂在一起，使得算法实现过程中有很多细节问题需要解决。本文详细叙述了GC算法实现细节，为进一步研究GC及其相关方法奠定基础。 　　1 优化模型建立 　　模型建立是GC算法的重要部分，该算法建立了类似MAP-MRF的估计方程，既考虑节点集合P，也包含每个节点邻域中的边缘{p，q}。先定义一个分割A=（Alt;subgt;1lt;/subgt;，…，Ap，…，A|P|），A是一个二值向量，其中每个元素属于前景或背景。然后可以通过A的代价方程E（A）表示区域和边缘的约束特性： 　　E（A）=λ×R（A）+B（A）（1） 　　 其中，R（A）是区域特性，B（A）是边缘特性。λ是一个大于0的值，它调节了R（A）和B（A）在整个表达式中的重要性。 　　R（A）=∑p∈NRp（Ap） 　　（2） 　　 它表示了区域项R（A）是每一个节点属性的总和。一般在目标提取任务中，Rp（?）可以认为是节点p和前景、背景联系的紧密程度，此时Rp（?）包含了两个值。 　　B（A）=∑{p，q}∈NB{p，q}?δ（Ap，aq）（3） 　　 其中，N是邻域系统，如四邻域。 　　δ（Ap，Aq）=1 if A”p≠Aq 　　0 otherwize 　　（4） 　　 即：如果Ap和Aq分类相同则忽略它们之间的边缘，它的直观理解是相同性质的节点应该处于相同区域。边缘项B（A）是所有边性质B{p，q}的总和。B{p，q}是一个大于零的项，它惩罚了节点p和q的不连续性。 　　2 从图像到图的转换 　　 假设待分割的图像为G=（V，E），其中V为顶点集合，一个像素为一个顶点，E为边的集合，相邻两个像素之间有一条边。图1（a）表示的是一个图像，其中椭圆代表顶点，实线代表边。图像中一般只看到像素，边是一个隐含的关系，图1（a）中的边是节点和其四邻域中点之间的关系。 　　图1 图像到图的转换 　　 GC算法将图像转换为图时保留了所有的节点和边，另外还添加了两个节点T和S。T表示背景，而S表示前景。这两个节点和原来所有节点之间都有联系，表示原来的节点和这两节点的相似程度。 　　 因此，将图像转变成图有这样几个部分：①保留原有的像素节点;②所有节点与其相邻节点（如四邻域）之间为边，文献中将这种边称为n-link。值得注意的是，n-link的边是无向的，若转变为有向图则可以添加双向相同的边;③添加背景节点T和前景节点S，并与原有的像素节点建立边，文献中将这种边称为t-link。 　　 图1（b）表示的是由一个3×3的图像转化而成的图。 　　3 优化模型的具体函数实现 　　图建立好后，就可以将图和式（1）建立的优化模型联系起来，最后通过优化算法求解图1（b）中标记为“cut”的“最小割”（min cut），将前景和背景分开。下面是优化模型式（1）中具体使用的函数，这一步能够将理论和实际联系起来。 　　3.1 区域项 　　设“obj”表示前景，“bkg”表示背景，则Rp（?）可以表示为一种对数形式： 　　Rp（”obj”）=-lnPr（Ip|obj）（5） 　　Rp（”bkg”）=-lnPr（Ip|bkj）（6） 　　 Pr（?）表示概率，式（5）和（6）分别表示的是节点属于前景、背景的概率，并