CN101582168A - 基于模糊连接度的抠图样本集构造方法 - Google Patents

基于模糊连接度的抠图样本集构造方法 Download PDF

Info

Publication number
CN101582168A
CN101582168A CNA2009100627318A CN200910062731A CN101582168A CN 101582168 A CN101582168 A CN 101582168A CN A2009100627318 A CNA2009100627318 A CN A2009100627318A CN 200910062731 A CN200910062731 A CN 200910062731A CN 101582168 A CN101582168 A CN 101582168A
Authority
CN
China
Prior art keywords
pixel
prime
background
prospect
border
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.)
Granted
Application number
CNA2009100627318A
Other languages
English (en)
Other versions
CN101582168B (zh
Inventor
何发智
吴玉娥
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changshu Zijin Intellectual Property Service Co ltd
Original Assignee
Wuhan University WHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN2009100627318A priority Critical patent/CN101582168B/zh
Publication of CN101582168A publication Critical patent/CN101582168A/zh
Application granted granted Critical
Publication of CN101582168B publication Critical patent/CN101582168B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)

Abstract

本发明公开了一种基于模糊连接度的抠图样本集构造方法:根据用户提供的区域划分,分别构造像素相对于已知前景、已知背景的相似权重图;求解未知像素与已知前景和背景边界的模糊连接度,并使用树形结构记录最强路径;遍历树形结构,获取每个未知像素通过最强路径关联的前景和背景边界根节点;以根节点为中心搜索邻近已知像素,构造未知像素对应的前景和背景样本集。应用本方法可以有效地处理抠图中未知像素的样本集的构造问题,与已有的抠图采样方法相比,本发明的方法大大降低了对用户提供的区域划分的依赖,样本集构造更加精确、鲁棒,适用于不同的用户交互情形,同时,能够对样本集进行快速有效的更新,具有很好的实用价值。

Description

基于模糊连接度的抠图样本集构造方法
技术领域
本发明属于数字图像处理技术领域,特别是涉及基于模糊连接度的抠图样本集构造方法。
背景技术
所谓抠图,就是通过恢复充分或部分的像素覆盖,将前景对象从背景中精确地分离出来。它是图像、视频编辑和电影制作中的关键技术,其结果可用于与其他背景合成,产生新的影像或特殊效果。图像抠图方法假设一幅输入图像I是前景图像If和背景图像Ib的合成。输入图像I中像素i的观测颜色Ii是其对应的前景颜色Ii f和背景颜色Ii b按照如下合成方程的线性组合
I i = α i I i f + ( 1 - α i ) I i b
αi是i的前景不透明度,αi∈[0,1]。由于仅仅已知图像的观测颜色Ii,Ii f、Ii b和αi均为未知,抠图问题本质上是一个欠约束的问题,因此,几乎所有的抠图方法都要求用户提供一定的约束信息,将图像划分为三个区域:前景区域F,背景区域B和未知区域U。
近几年,国内外研究者从不同的角度提出了大量的解决抠图问题的方法,其中一类重要的方法是通过使用与未知区域U中的未知像素邻近的已知前景与背景区域的统计信息估算未知像素的前景颜色,背景颜色及前景不透明度值,该方法的潜在假设是未知像素的真实前景、背景颜色包含在其对应的前景、背景样本集中。颜色样本被进一步处理,以估算出未知像素的最优前景、背景颜色。一旦前景颜色Ii f,背景颜色Ii b被确定,未知像素的前景不透明度值αi能够很容易地从合成公式中求解出。因此如何准确构造未知像素的样本集成为影响这类方法性能优劣的至关重要的一步。
在抠图问题中,未知像素的最优样本应满足两个重要条件:1.前景样本与背景样本的颜色及未知像素的观测颜色Ii符合(1)中的线性关系;2.前景样本、背景样本的颜色应尽可能接近于未知像素的观测颜色Ii。条件1的合理性显而易见。条件2则是构造样本集时必须满足的要求,由于未知区域中大部分像素为充分的前景像素(即前景不透明度αi为1的像素)和充分的背景像素(即前景不透明度αi为0的像素),少量的混合像素(即前景不透明度αi介于0,1之间的像素)的颜色则由邻近的前景和背景颜色混合而成,因此未知像素的前景样本和背景样本应尽可能选取与其颜色最相似的像素。然而,以前的大部分实现抠图采样的方法其主要思想是邻接最近优先原则,即总是选取在空间位置上离待考察的未知像素近的已知像素作为其样本,其潜在假设是两个像素越邻近则它们所对应的抠图参数越相似,该策略在用户提供的区域划分足够精确的情况下是有效的,然而,一旦用户提供的区域划分较粗糙时,就容易错误地选择到与未知像素相对位置距离近但颜色差异却非常大的样本,这类样本容易导致在后续抠图处理中产生较大的误差。
早期的基于采样的抠图系统,如贝叶斯方法采取构造以未知像素为中心的滑动窗口搜集样本集,将滑动窗口覆盖的已知像素作为对应未知像素的样本,然而这种方法由于受到窗口大小的限制,容易导致低效的采样处理,且可能无法充分地跨越前景和背景的颜色变化。Robust抠图因此提出了沿已知前景和背景边界采样的策略,与局部窗口采样方法相比,该方法能够更好地跨越已知区域的颜色变化。但是这两种方法均不同程度地受到用户提供的区域划分的精确度的约束,当用户提供的区域划分比较粗糙时,对远离前景边界(即前景区域与未知区域之间的分界线)或背景边界(即背景区域与未知区域之间的分界线)的未知像素,两种方法都可能采样到空间位置上距离近的错误的样本。最近,Improving color modeling for alpha matting方法假设前景对象在空间上连通的,通过求解与测地距离相对应的最大权重路径搜索与未知像素最相似的前景样本,该方法在一定程度上降低了空间位置对前景样本集正确构造的限制,但是由于测地距离的计算本身的局限性,该方法容易为包含在具有细长结构的区域中的未知像素构造出错误的样本集。此外,该方法认为背景区域由于前景对象的遮挡通常是不连通的,因此对背景样本集的构造仍然简单地采取了沿背景边界采样的策略,事实上,对背景颜色分布模式复杂的情形,这种处理是不足够的,往往导致采样错误的发生。
发明内容
本发明目的在于提供一种基于模糊连接度的抠图样本集构造方法,以克服在用户提供的区域划分粗糙的情况下容易产生采样错误的问题。
本发明提供的一种基于模糊连接度的抠图样本集构造方法,包含以下步骤:
步骤1,根据用户提供的区域划分,分别构造图像中像素相对于已知前景、已知背景的相似权重图Af、Ab
步骤2,根据相似权重图Af、Ab求解未知像素与已知的前景边界和背景边界的模糊连接度,并使用树形结构记录最强路径;
步骤3,遍历树形结构,获取每个未知像素通过最强路径关联的前景边界像素和背景边界像素,并分别记为前景边界根节点和背景边界根节点;
步骤4,以前景边界根节点为中心搜索邻近已知前景像素,构造每个未知像素对应的前景样本集;以背景边界根节点为中心搜索邻近已知背景像素,构造每个未知像素对应的背景样本集。
而且,像素相对于已知前景、已知背景的相似权重图Af、Ab按照如下步骤构造,
步骤1.1,根据用户通过交互信息提供的区域划分,将图像划分为三个区域:前景区域F,背景区域B和未知区域U,前景区域F即为图像的已知前景,背景区域B即为图像的已知背景,直接与未知区域U中的未知像素相邻接的前景/背景像素构成了前景/背景边界,位于前景/背景边界上的像素即为前景/背景边界像素;
步骤1.2,分别对已知前景的颜色和已知背景的颜色建模,定义图像中所有像素本身及邻接像素间相对于已知前景和背景的颜色相似度的权重函数,构造出相似权重图Af及Ab,定义Ao(p,q),o∈{f,b},表示相似权重图中像素p,q之间的边的权重值。
而且,所述模糊连接度定义如下:
(1)定义两个像素间路径的连接强度,
在图像内任意两个像素s和t之间存在着大量的路径,构成路径集合P,每条路径γ∈P由一系列连续邻接的节点p0,p1,p2,…,pn-1,pn连接而成,表示为γ(s→t)={p0,p1,p2,…,pn-1,pn},其中起始的节点p0为s,结束的节点pn为t,且 A o i = 0 · · · , n - 1 ( p i , p i + 1 ) > 0 , o∈{f,b},
Figure A20091006273100142
是节点pi和pi+1之间的边的相似权重,由相似权重图Af或Ab中取得;
路径γ(s→t)的连接强度记为
str o ( γ ( s → t ) ) min i = 0 · · · , n - 1 A o ( p i , p i + 1 ) , o ∈ { f , b }
(2)像素s和t的模糊连接度定义为两者间所有路径中最大的连接强度,记为
FC o ( s , t ) max γ ∈ P str o ( γ ( s → t ) ) , o ∈ { f , b } ;
(3)图像中像素到前景边界或背景边界的模糊连接度,根据该像素与前景边界像素或背景边界像素的模糊连接度得到。
而且,步骤2具体实施方式如下,
(1)求解未知像素与前景边界的模糊连接度和最强路径,步骤如下,
步骤2.11,设置相关变量并初始化,包括有:
设置像素与前景边界间的模糊连接度FCf,初始化时将已知前景中所有的前景像素与前景边界间的FCf值设置为1,图像中所有非前景像素与前景边界间的FCf值设为0;
设置以模糊连接度FCf值为优先权的优先队列Qf,初始化时将所有前景边界像素加入到优先队列Qf中;
设置标识位Flagf,当Flagf值为“1”时表示对应像素可以被插入优先队列Qf中,为“0”则不能插入优先队列Qf中;初始化时前景边界像素与未知像素的Flagf值设为1,图像中剩余像素的Flagf值设为0;
设置树形结构Tf,用于记录前景边界像素和未知像素的父节点和孩子节点;初始化时把所有前景边界像素加入到树形结构Tf中,且对应的父节点设为-1,孩子节点设为空;并且把所有未知像素加入到树形结构Tf中,且对应的父节点和孩子节点设为空;
步骤2.12,当优先队列Qf不为空时,提取出优先队列Qf中具有最大的FCf值的像素p,继续执行步骤2.13、2.14;当Qf为空时,求解终止;
步骤2.13,对像素p的每个邻接像素q:
1)如果由相似权重图Af中取得的像素p,q之间的边的相似权重Af(p,q)>0,且Flagf(q)等于1,则执行后续操作;
2)比较像素p与前景边界的模糊连接度FCf(p)和Af(p,q)的大小;
3)令 FC min f { FC f ( p ) , A f ( p , q ) } ;
FCf(q)为像素q与前景边界的模糊连接度,如果 FC min f > FC f ( q ) ) , 则执行后续操作;
4)设置 FC f ( q ) = FC min f ;
5)如果q在树形结构Tf中,则
更新q的父节点指针为p;
把q加为p的孩子节点;
如果q不在树形结构Tf中,则
把q插入到树形结构Tf中,q的父节点指针指向p;
把q加为p的孩子节点;
把节点q插入到队列Qf中;
步骤2.14,将p从Qf中删除,返回步骤2.12;
(2)求解未知像素与背景边界的模糊连接度和最强路径,步骤如下,
步骤2.21,设置相关变量并初始化,包括有:
设置像素与背景边界间的模糊连接度FCb,初始化时将已知背景中的所有背景像素与背景边界间的FCb值设置为1,图像中所有非背景像素与背景边界间的FCb值设为0;
设置以模糊连接度FCb值为优先权的优先队列Qb,初始化时将所有背景边界像素加入到优先队列Qb中;
设置标识位Flagb,当Flagb值为“1”时表示对应像素可以被插入优先队列Qb中,为“0”则不能插入优先队列Qb中;初始化时背景边界像素与未知像素的Flagb值设为1,图像中剩余像素的Flagb值设为0;
设置树形结构Tb,用于记录背景边界像素和未知像素的父节点和孩子节点;初始化时把所有背景边界像素加入到树形结构Tb中,且对应的父节点设为-1,孩子节点设为空;并且把所有未知像素加入到树形结构Tb中,且对应的父节点和孩子节点设为空;
步骤2.22,当优先队列Qb不为空时,提取出优先队列Qb中具有最大的FCf值的像素p,继续执行步骤2.23、2.24;当Qb为空时,求解终止;
步骤2.23,对像素p的每个邻接像素q:
1)如果由相似权重图Ab中取得的像素p,q之间的边的相似权重Ab(p,q)>0,且Flagb(q)等于1,则执行后续操作;
2)比较p与背景边界的模糊连接度FCb(p)和Ab(p,q)的大小;
3)令 FC min b = { FC b ( p ) , A b ( p , q ) } ,
FCb(q)为像素q与背景边界的模糊连接度,如果 FC min b > FC b ( q ) ) , 则执行后续操作;
4)设置 FC b ( q ) = FC min b ;
5)如果q在树形结构Tb中,则
更新q的父节点指针为p;
把q加为p的孩子节点;
如果q不在树形结构Tb中,则
把q插入到树形结构Tb中,q的父节点指针指向p;
把q加为p的孩子节点;
把节点q插入到队列Qb中;
步骤2.24,将p从Qb中删除,返回步骤2.22。
而且,用户通过交互提供新的笔划以补充区域划分时,只对与新的笔划所产生边界之间存在连通路径的未知像素重新处理,更新前景和背景样本集;
(1)当新的前景笔划加入后,具体处理方式如下,
1)如果新加入的前景笔划包含了新的前景颜色,则对新的已知前景的颜色重新建模,定义图像中所有像素本身及邻接像素间相对于已知前景的颜色相似度的权重函数,更新相似权重图Af,然后转向步骤2);如果新加入的前景笔划未包含新的前景颜色,直接转向步骤2);
2)计算新加入的前景笔划所产生前景边界F”与原有的前景边界F’间的模糊连接度值FC(F’,F”);
3)根据FC(F’,F”)将新的前景笔划加入后所产生新的未知区域U’划分为三个子集合:
Ω U ′ 1 = { i | i ∈ U ′ , FC f ( F ′ , i ) > FC ( F ′ , F ′ ′ ) }
Ω U ′ 2 = { i | i ∈ U ′ , FC f ( F ′ , i ) = FC ( F ′ , F ′ ′ ) }
&Omega; U &prime; 3 = { i | i &Element; U &prime; , FC f ( F &prime; , i ) < FC ( F &prime; , F &prime; &prime; ) }
其中FCf(F′,i)表示新的未知区域U’中任意像素i与原有的前景边界F’的模糊连接度;
4)包含在集合ΩU′ 1和ΩU′ 3中的未知像素,其对应的FCf及前景样本集仍使用原值;对包含在集合ΩU′ 2中的未知像素则按照步骤5)做进一步处理;
5)如果集合ΩU′ 2不为空,则:
首先,初始化相关变量:设置被新的前景笔划所覆盖的未知像素的FCf为1,ΩU′ 2中未知像素的FCf为0;设置ΩU′ 2中未知像素的Flagf位为1,图像中剩余像素的Flagf位为0;把新的前景边界F”中的像素加入到队列Qf中;
然后,重新计算ΩU′ 2中未知像素与新的前景边界F”间的FCf和最强路径;
最后,根据计算出的结果做如下处理:如果未知像素的FCf值为0,则该未知像素的FCf及前景样本集仍使用原值;否则,该未知像素与F”之间通过区域ΩU′ 2是连通的,其FCf设置为重新计算出的FCf值,并根据新的最强路径更新该像素的前景样本集;
(2)当新的背景笔划加入后,具体处理方式如下,
1)如果新加入的背景笔划包含了新的背景颜色,则对新的已知背景的颜色重新建模,定义图像中所有像素本身及邻接像素间相对于已知背景的颜色相似度的权重函数,更新相似权重图Ab,然后转向步骤2);如果新加入的背景笔划未包含新的背景颜色,直接转向步骤2);
2)计算新加入的背景笔划所产生背景边界B”与原有的背景边界B’间的模糊连接度值FC(B’,B”);
3)根据FC(B’,B”)将新的背景笔划加入后所产生新的未知区域U’划分为三个子集合:
&Omega; U &prime; 1 = { i | i &Element; U &prime; , FC b ( B &prime; , i ) > FC ( B &prime; , B &prime; &prime; ) }
&Omega; U &prime; 2 = { i | i &Element; U &prime; , FC b ( B &prime; , i ) = FC ( B &prime; , B &prime; &prime; ) }
&Omega; U &prime; 3 = { i | i &Element; U &prime; , FC b ( B &prime; , i ) < FC ( B &prime; , B &prime; &prime; ) }
其中FCb(B′,i)表示新的未知区域U’中任意像素i与原有的背景边界B’的模糊连接度;
4)包含在集合ΩU′ 1和ΩU′ 3中的未知像素,其对应的FCb及背景样本集仍使用原值;对包含在集合ΩU′ 2中的未知像素则按照步骤5)做进一步处理;
5)如果集合ΩU′ 2不为空,则:
首先,初始化相关变量:设置被新的背景笔划所覆盖的未知像素的FCb为1,ΩU′ 2中未知像素的FCb为0;设置ΩU′ 2中未知像素的Flagb位为1,图像中剩余像素的Flagb位为0;把新的背景边界B”中的像素加入到队列Qb中;
然后,重新计算ΩU′ 2中未知像素与新的背景边界B”间的FCb和最强路径;
最后,根据计算出的结果做如下处理:如果未知像素的FCb值为0,则该未知像素的FCb及背景样本集仍使用原值;否则,该未知像素与F”之间通过区域ΩU′ 2是连通的,其FCb设置为重新计算出的FCb值,并根据新的最强路径更新该像素的背景样本集。
而且,(1)计算新加入的前景笔划所产生的前景边界F”与原有的前景边界F’间的模糊连接度FC(F’,F”),具体实现方式如下
前景边界F”与原有的前景边界F’间的模糊连接度为两个像素集合间的模糊连接度,定义为:
FC ( F &prime; , F &prime; &prime; ) = max f &prime; &Element; F &prime; max f &prime; &prime; = F &prime; &prime; FC ( f &prime; , f &prime; &prime; )
考察F”中每个像素i与F’的模糊连接度FCf(F’,i),并从中挑选出最大的FCf(F’,i)值作为FC(F’,F”),即
FC ( F &prime; , F &prime; &prime; ) = max i &Element; F &prime; &prime; FC f ( F &prime; , i )
(2)计算新加入的背景笔划所产生的背景边界B”与原有的背景边界B’间的模糊连接度FC(B’,B”),具体实现方式如下
背景边界B”与原有的背景边界B’间的模糊连接度为两个像素集合间的模糊连接度,定义为:
FC ( B &prime; , B &prime; &prime; ) = max b &prime; &Element; B &prime; max b &prime; &prime; &Element; B &prime; &prime; FC ( b &prime; , b &prime; &prime; )
考察B”中每个像素i与B’的模糊连接度FCb(B’,i),并从中挑选出最大的FCb(B’,i)值作为FC(B’,B”),即
FC ( B &prime; , B &prime; &prime; ) = max i &Element; B &prime; &prime; FC b ( B &prime; , i ) .
而且,步骤3具体实施方式为,遍历树形结构,通过孩子节点或父节点指针追溯未知像素到前景或背景边界的最强路径,获取最强路径所关联的前景或背景边界像素。
而且,步骤4具体实施方式为,获取某个未知像素对应的前景边界根节点后,以该根节点为中心沿前景边界搜索m个前景像素,所得的m个前景像素即前景样本构成了前景样本集,m为用户指定的每个未知像素的前景样本的个数;获取某个未知像素对应的背景边界根节点后,以该根节点为中心沿背景边界搜索n个背景像素,所得的n个背景像素即背景样本构成了背景样本集,其中n为用户指定的每个未知像素的背景样本的个数。
本发明还要求保护所述抠图样本集构造方法在图像处理领域或视频处理领域中的应用。
本发明通过计算未知像素到前景和背景边界的最强路径,搜集与未知像素在颜色空间最相似的样本集,大大降低了对用户提供的区域划分的依赖,样本集构造过程更加精确、鲁棒。并且当新的用户笔划加入时,能够对样本集进行快速有效的更新。由于其鲁棒性和灵活性,本发明的方法可以被应用到多种抠图系统中,产生更加先进的抠图方法,也可以在视频技术领域发挥作用。
附图说明
图1是本发明实施例中前景样本集构造流程图。
图2是本发明实施例中前景样本集更新的流程示意图。
图3是本发明实施例中用户提供的较精确的区域划分。
图4是本发明实施例中用户使用少量笔划提供的较粗糙的区域划分。
具体实施方式
本发明提供的基于模糊连接度的抠图样本集构造方法,其本质思想是选取与未知像素在颜色空间最相似的已知像素作为其样本,本方法通过计算未知像素到前景和背景边界的模糊连接度及对应的最强路径,以最强路径所关联的已知边界像素为中心搜集未知像素的样本集。
本发明的基于模糊连接度的抠图样本集构造方法,包括以下步骤:
步骤1,根据用户提供的区域划分,分别构造图像中像素相对于已知前景、已知背景的相似权重图Af、Ab。准确地,是像素相对于已知前景的相似权重图Af,像素相对于已知背景的相似权重图Ab
步骤2,根据相似权重图Af、Ab求解未知像素与前景边界和背景边界的模糊连接度,并使用树形结构记录最强路径。
具体是根据相似权重图Af求解未知像素与前景边界的模糊连接度,根据相似权重图Ab求解未知像素与背景边界的模糊连接度。
步骤3,遍历树形结构,获取每个未知像素通过最强路径关联的前景边界像素和背景边界像素,并分别记为前景边界根节点和背景边界根节点。
步骤4,以前景边界根节点为中心搜索邻近已知前景像素,构造每个未知像素对应的前景样本集;以背景边界根节点为中心搜索邻近已知背景像素,构造每个未知像素对应的背景样本集。
以下结合实施例详细说明各步骤实现方式:
步骤1
用户提供的区域划分,具体实施时可以是由用户提供较多的交互信息指定大部分前景和背景实现区域划分,也可以是由用户提供少量笔划标识部分前景和背景像素进行区域划分。“由用户提供较多的交互信息指定大部分前景和背景”是指用户交互量较多的情况,指定了图像中的大部分前景和背景像素,“仅使用少量笔划标识部分前景和背景像素”对应于粗糙的用户交互,在该情形下用户仅通过简单的几笔指定少量的前景和背景像素。本发明既支持充分交互的情况,也支持交互粗糙的情况。将图像划分为三个区域,即前景区域F,背景区域B和未知区域U之后,前景和背景确定为已知信息,前景边界和背景边界也就确定成为已知信息了。因此,前景区域F即为图像的已知前景,背景区域B即为图像的已知背景,直接与未知区域U中的未知像素相邻接的前景/背景像素构成了前景/背景边界,位于前景/背景边界上的像素即为前景/背景边界像素。
构造相似权重图属于现有技术,本发明提出具体实现方式为:分别对已知前景的颜色和已知背景的颜色建模,定义图像中所有像素本身及邻接像素间相对于已知前景和背景的颜色相似度的权重函数,构造出相似权重图Af及Ab,定义Ao(p,q),o∈{f,b},表示相似权重图中像素p,q之间的边的权重值。具体实施时,也可以采用其它方式得到出相似权重图Af及Ab,例如通过考察像素颜色属性和空间位置属性上的距离得到。
实施例中,将一幅输入图像抽象为4-邻接的图,为前景和背景分别构建颜色分布模型——高斯混合模型GMM,并使用FuzzyMatte方法构造Af、Ab。高斯混合模型GMM和FuzzyMatte方法为现有技术,为便于实施参考,本发明提供实施例相关说明,包括使用FuzzyMatte方法中定义的如下形式的权重函数,构造Af、Ab
Figure A20091006273100241
其中
Figure A20091006273100242
表示两个像素之间的颜色相似度,其定义为:
Figure A20091006273100243
公式中:I(p1)、I(p2)分别表示两个邻接像素p1,p2的颜色向量,∑max o是对GMM中的每个高斯分量平均三个颜色通道的方差值后,具有最大的平均方差的前景或背景高斯分量的协方差矩阵,exp()表示取指数,上标T表示矩阵转置。本发明采用f标示前景,b标示背景。
Aφ o(p1,p2)像素与已知前景或背景颜色间的相似度,以计算像素到已知前景的颜色相似度Aφ f(p1,p2)为例,其步骤为:
首先定义一个像素p与前景颜色间的相似度:
S f ( p ) max i exp ( - 1 2 [ I ( p ) - m i f ] T ( &Sigma; i f ) - 1 [ I ( p ) - m i f ] )
公式中:i表示前景高斯混合模型GMM中的高斯分量的索引,mi f,∑i f分别是第i个前景高斯分量所对应的均值向量和协方差矩阵。
Figure A20091006273100245
符号表示在i取得不同值时,所有对应的exp()的取值中的最大值。则
Figure A20091006273100246
这里,
W min f ( p 1 , p 2 ) = min [ S f ( p 1 ) , S f ( p 2 ) ]
W max f ( p 1 , p 2 ) = max [ S b ( p 1 ) , S b ( p 2 ) ]
像素与已知背景颜色间的相似度Aφ b(p1,p2)的计算与Aφ f(p1,p2)的定义相似。
Ao(p1,p2)的计算公式中的符号λ是权衡两种相似度衡量
Figure A20091006273100251
与Aφ o(p1,p2)之间的相对重要性的参数,λ∈[0,1]。
步骤2
本发明为了求解未知像素与已知前景和背景边界间的模糊连接度,首先引入两个像素间路径强度的定义,然后设定两个像素间的模糊连接度的定义。在一幅图像中像素s和t之间存在着大量的路径P,每条路径γ∈P由一系列连续的邻接节点连接而成,γ(s→t)={p0,p1,p2,…,pn-1,pn},p0=r,pn=t,且 A o i = 0 &CenterDot; &CenterDot; &CenterDot; , n - 1 ( p i , p i + 1 ) > 0 . 路径γ(s→t)的强度被定义为:
str o ( &gamma; ( s &RightArrow; t ) ) min i = 0 &CenterDot; &CenterDot; &CenterDot; , n - 1 A o ( p i , p i + 1 ) , o &Element; { f , b }
则s、t的模糊连接度定义为两者间所有路径中最大的连接强度,即
FC o ( s , t ) max &gamma; &Element; P str o ( &gamma; ( s &RightArrow; t ) ) , o &Element; { f , b } ;
计算两个像素间的模糊连接度本质上是一个单源最短路径问题。
具体求解未知像素到已知前景和背景边界的模糊连接度,并使用树形结构记录最强路径,包括前景和背景两部分,需要分别处理。
实施例中求解未知像素到前景边界的模糊连接度和最强路径,其步骤为:
步骤2.11,设置相关变量并初始化,包括有:
设置像素与前景边界间的模糊连接度FCf,初始化时将已知前景中所有的前景像素与前景边界间的FCf值设置为1,图像中所有非前景像素与前景边界间的FCf值设为0;
设置以模糊连接度FCf值为优先权的优先队列Qf,初始化时将所有前景边界像素加入到优先队列Qf中;
设置标识位Flagf,当Flagf值为“1”时表示对应像素可以被插入优先队列Qf中,为“0”则不能插入优先队列Qf中;初始化时前景边界像素与未知像素的Flagf值设为1,图像中剩余像素的Flagf值设为0;
设置树形结构Tf,用于记录前景边界像素和未知像素的父节点和孩子节点;初始化时把所有前景边界像素加入到树形结构Tf中,且对应的父节点设为-1,孩子节点设为空;并且把所有未知像素加入到树形结构Tf中,且对应的父节点和孩子节点设为空;
步骤2.12,当优先队列Qf不为空时,提取出优先队列Qf中具有最大的FCf值的像素p,继续执行步骤2.13、2.14;当Qf为空时,求解终止;
步骤2.13,对像素p的每个邻接像素q:
1)如果由相似权重图Af中取得的像素p,q之间的边的相似权重Af(p,q)>0,且Flagf(q)等于1,则执行后续操作;
2)比较像素p与前景边界的模糊连接度FCf(p)和Af(p,q)的大小;
由于像素p与前景边界的模糊连接度根据像素与前景边界像素的模糊连接度得出,FCf(p)可定义为 FC f ( p ) = FC f ( F &prime; , p ) = max p &prime; &Element; F &prime; FC f ( p &prime; , p ) , 其中F’标记前景边界,p’标记前景边界F’中的任一个像素。即考察像素p与前景边界F’中所有前景边界像素的模糊连接度,取其中最大者为像素p与前景边界F’的模糊连接度。
3)令 FC min f { FC f ( p ) , A f ( p , q ) } ;
FCf(q)为像素q与前景边界的模糊连接度,如果 FC min f > FC f ( q ) ) , 则执行后续操作;
与FCf(p)获取过程相同,
因此FCf(q)可定义为 FC f ( q ) = FC f ( F &prime; , q ) = max q &prime; &Element; F &prime; FC f ( q &prime; , q ) , 其中F’标记前景边界,q’标记前景边界F′中的任一个像素。
4)设置 FC f ( q ) = FC min f ;
5)如果q在树形结构Tf中,则
更新q的父节点指针为p;
把q加为p的孩子节点;
如果q不在树形结构Tf中,则
把q插入到树形结构Tf中,q的父节点指针指向p;
把q加为p的孩子节点;
把节点q插入到队列Qf中;
步骤2.14,将p从Qf中删除,返回步骤2.12。
实施例中求解未知像素到背景边界的模糊连接度和最强路径,步骤如下:
步骤2.21,设置相关变量并初始化,包括有:
设置像素与背景边界间的模糊连接度FCb,初始化时将已知背景中的所有背景像素与背景边界间的FCb值设置为1,图像中所有非背景像素与背景边界间的FCb值设为0;
设置以模糊连接度FCb值为优先权的优先队列Qb,初始化时将所有背景边界像素加入到优先队列Qb中;
设置标识位Flagb,当Flagb值为“1”时表示对应像素可以被插入优先队列Qb中,为“0”则不能插入优先队列Qb中;初始化时背景边界像素与未知像素的Flagb值设为1,图像中剩余像素的Flagb值设为0;
设置树形结构Tb,用于记录背景边界像素和未知像素的父节点和孩子节点;初始化时把所有背景边界像素加入到树形结构Tb中,且对应的父节点设为-1,孩子节点设为空;并且把所有未知像素加入到树形结构Tb中,且对应的父节点和孩子节点设为空;
步骤2.22,当优先队列Qb不为空时,提取出优先队列Qb中具有最大的FCb值的像素p,继续执行步骤2.23、2.24;当Qb为空时,求解终止;
步骤2.23,对像素p的每个邻接像素q:
1)如果由相似权重图Ab中取得的像素p,q之间的边的相似权重(p,q)>0,且Flagb(q)等于1,则执行后续操作;
2)比较p与背景边界的模糊连接度FCb(p)和Ab(p,q)的大小;
由于像素p与背景边界的模糊连接度根据像素与背景边界像素的模糊连接度得出,FCb(p)可定义为 FC b ( p ) = FC b ( B &prime; , p ) = max p &prime; &Element; B &prime; FC b ( p &prime; , p ) , 其中B’标记背景边界,p’标记背景边界中B’中的任一个像素。即考察像素p与背景边界B’中所有前景边界像素的模糊连接度,取其中最大者为像素p与背景边界B’的模糊连接度。
3)令 FC min b = { FC b ( p ) , A b ( p , q ) } ,
FCb(q)为像素q与背景边界的模糊连接度,如果 FC min b > FC b ( q ) ) , 则执行后续操作;
与FCb(p)获取过程相同,
因此FCb(q)可定义为 FC b ( q ) = FC b ( B &prime; , q ) = max q &prime; &Element; B &prime; FC b ( q &prime; , q ) 其中B’标记背景边界,q’标记背景边界B’中的任一个像素。
4)设置 FC b ( q ) = FC min b ;
5)如果q在树形结构Tb中,则
更新q的父节点指针为p;
把q加为p的孩子节点;
如果q不在树形结构Tb中,则
把q插入到树形结构Tb中,q的父节点指针指向p;
把q加为p的孩子节点;
把节点q插入到队列Qb中;
步骤2.24,将p从Qb中删除,返回步骤2.22。
步骤3
本发明提供进一步技术方案:当所有未知像素与前景边界和背景边界间的模糊连接度FCf,FCb求解出后,遍历对应的树形结构,通过孩子节点或父节点指针追溯每个未知像素到前景或背景边界的最强路径,获取最强路径所关联的前景或背景边界像素,即为最强路径所关联的前景或背景边界根节点。根节点即与对应未知像素在颜色空间最相似的已知像素。
求解模糊连接度的过程实际上也是确定最强路径的过程,因此,为了记录路径本发明给对应的像素分配两个变量:即步骤2中所说的“父节点”“孩子节点”指针变量,从起始节点到目的节点各个邻接节点相互链接构成了一条路径,对一个具体的像素i,其“父节点”是这条路径上i的上一个邻接节点,其“孩子节点”是这条路径上i的下一个邻接节点。由于在求解模糊连接度的过程中,已经记录了路径中每个节点的“父节点”和“孩子节点”信息(起始节点即前景/背景边界像素节点的“父节点”为-1,目的节点即当前被考察的未知像素的“孩子节点”为空)。因此,模糊连接度求解过程结束后,从起始点(“目的节点”)开始,通过“孩子节点”(“父节点“)指针找到下(上)一个节点,直到“孩子节点”(“父节点”)为空(-1),这样就很容易地确定了一条最强路径,并获取最强路径所关联的前景或背景边界根节点,根节点即与对应未知像素在颜色空间最相似的已知像素。
步骤4
本发明提供进一步技术方案:获取某个未知像素对应的前景边界根节点后,以该根节点为中心沿前景边界搜索前景像素,作为该未知像素的前景样本,搜索所得所有前景样本构成前景样本集。对背景样本集的构造过程与之类似,获取某个未知像素对应的背景边界根节点后,以该根节点为中心沿背景边界搜索背景像素,作为该未知像素的背景样本,搜索所得所有背景样本构成背景样本集。具体实施时,也可以根节点为中心在局部邻域内搜索前景像素,局部邻域意指以该根节点为中心所构成的区域,例如矩形或圆形区域。具体实施时,可以事先由用户根据需要指定每个未知像素的前景样本的个数m和背景样本的个数n。一般应用时建议m和n取值范围为3~10。
当用户提供新的交互信息时,可以按照交互信息重新进行区域划分,然后重复步骤1~4。但为了提高效率,本发明还提供了简便的样本集更新方案:用户通过交互信息提供新的笔划以补充区域划分时,只对与新的笔划所产生边界之间存在连通路径的未知像素重新处理,更新前景和背景样本集。用户提供的笔划可能为前景笔划或背景笔划。
实施例中,当新的前景笔划加入后,具体处理方式如下,
1)如果新加入的前景笔划包含了新的前景颜色,则对新的已知前景的颜色重新建模,定义图像中所有像素本身及邻接像素间相对于已知前景的颜色相似度的权重函数,更新相似权重图Af,然后转向步骤2);如果新加入的前景笔划未包含新的前景颜色,直接转向步骤2);
2)计算新加入的前景笔划所产生前景边界F”与原有的前景边界F’间的模糊连接度值FC(F’,F”);
3)根据FC(F’,F”)将新的前景笔划加入后所产生新的未知区域U’划分为三个子集合:
&Omega; U &prime; 1 = { i | i &Element; U &prime; , FC f ( F &prime; , i ) > FC ( F &prime; , F &prime; &prime; ) }
&Omega; U &prime; 2 = { i | i &Element; U &prime; , FC f ( F &prime; , i ) = FC ( F &prime; , F &prime; &prime; ) }
&Omega; U &prime; 3 = { i | i &Element; U &prime; , FC f ( F &prime; , i ) < FC ( F &prime; , F &prime; &prime; ) }
其中FCf(F′,i)表示新的未知区域U’中任意像素i与原有的前景边界F’的模糊连接度,可定义为 FC f ( F &prime; , i ) = max p &prime; &Element; F &prime; FC f ( p &prime; , i ) , 即考察像素i与的前景边界F’中所有前景边界像素的模糊连接度,取其中最大者为像素i与原有的前景边界F’的模糊连接度。
4)包含在集合ΩU′ 1和ΩU′ 3中的未知像素,其对应的FCf及前景样本集仍使用原值;对包含在集合ΩU′ 2中的未知像素则按照步骤5)做进一步处理;
5)如果集合ΩU′ 2不为空,则:
首先,初始化相关变量:设置被新的前景笔划所覆盖的未知像素的FCf为1,ΩU′ 2中未知像素的FCf为0;设置ΩU′ 2中未知像素的Flag位为1,图像中剩余像素的Flagf位为0;把新的前景边界F”中的像素加入到队列Qf中;
然后,重新计算ΩU′ 2中未知像素与新的前景边界F”间的FCf和最强路径;
最后,根据计算出的结果做如下处理:如果未知像素的FCf值为0,则该未知像素的FCf及前景样本集仍使用原值;否则,该未知像素与F”之间通过区域ΩU′ 2是连通的,其FCf设置为重新计算出的FCf值,并根据新的最强路径更新该像素的前景样本集。
实施例中,当新的背景笔划加入后,具体处理方式如下,
1)如果新加入的背景笔划包含了新的背景颜色,则对新的已知背景的颜色重新建模,定义图像中所有像素本身及邻接像素间相对于已知背景的颜色相似度的权重函数,更新相似权重图Ab,然后转向步骤2);如果新加入的背景笔划未包含新的背景颜色,直接转向步骤2);
2)计算新加入的背景笔划所产生背景边界B”与原有的背景边界B’间的模糊连接度值FC(B’,B”);
3)根据FC(B’,B”)将新的背景笔划加入后所产生新的未知区域U’划分为三个子集合:
&Omega; U &prime; 1 = { i | i &Element; U &prime; , FC b ( B &prime; , i ) > FC ( B &prime; , B &prime; &prime; ) }
&Omega; U &prime; 2 = { i | i &Element; U &prime; , FC b ( B &prime; , i ) = FC ( B &prime; , B &prime; &prime; ) }
&Omega; U &prime; 3 = { i | i &Element; U &prime; , FC b ( B &prime; , i ) < FC ( B &prime; , B &prime; &prime; ) }
其中FCb(B′,i)表示新的未知区域U’中任意像素i与原有的背景边界B’的模糊连接度,可定义为 FC b ( B &prime; , i ) = max p &prime; &Element; B &prime; FC b ( p &prime; , i ) 即考察像素i与的背景边界B’中所有前景边界像素的模糊连接度,取其中最大者像素i与原有的背景边界B’的模糊连接度。
4)包含在集合ΩU′ 1和ΩU′ 3中的未知像素,其对应的FCf及背景样本集仍使用原值;对包含在集合ΩU′ 2中的未知像素则按照步骤5)做进一步处理;
5)如果集合ΩU′ 2不为空,则:
首先,初始化相关变量:设置被新的背景笔划所覆盖的未知像素的FCb为1,ΩU′ 2中未知像素的FCb为0;设置ΩU′ 2中未知像素的Flagb位为1,图像中剩余像素的Flagb位为0;把新的背景边界B”中的像素加入到队列Qb中;
然后,重新计算ΩU′ 2中未知像素与新的背景边界B”间的FCb和最强路径;
最后,根据计算出的结果做如下处理:如果未知像素的FCb值为0,则该未知像素的FCb及背景样本集仍使用原值;否则,该未知像素与F”之间通过区域ΩU′ 2是连通的,其FCb设置为重新计算出的FCb值,并根据新的最强路径更新该像素的背景样本集。
为了便于实施,提供实施例中新旧边界间模糊连接度的具体过程如下:
(1)计算新加入的前景笔划所产生前景边界F”与原有的前景边界F’间的模糊连接度FC(F’,F”),具体实现方式如下
前景边界F”与前景边界F’间的模糊连接度为两个像素集合间的模糊连接度,定义为:
FC ( F &prime; , F &prime; &prime; ) = max f &prime; &Element; F &prime; max f &prime; &prime; = F &prime; &prime; FC ( f &prime; , f &prime; &prime; )
考察F”中每个像素i与F’的模糊连接度FCf(F’,i),并从中挑选出最大的FCf(F’,i)值作为FC(F’,F”),即
FC ( F &prime; , F &prime; &prime; ) = max i &Element; F &prime; &prime; FC f ( F &prime; , i )
(2)计算新加入的背景笔划所产生背景边界B”与原有的背景边界B’间的模糊连接度FC(B’,B”),具体实现方式如下
背景边界B”与原有的背景边界B’间的模糊连接度为两个像素集合间的模糊连接度,定义为:
FC ( B &prime; , B &prime; &prime; ) = max b &prime; &Element; B &prime; max b &prime; &prime; &Element; B &prime; &prime; FC ( b &prime; , b &prime; &prime; )
考察B”中每个像素i与B’的模糊连接度FCb(B’,i),并从中挑选出最大的FCb(B’,i)值作为FC(B’,B”),即
FC ( B &prime; , B &prime; &prime; ) = max i &Element; B &prime; &prime; FC b ( B &prime; , i ) .
本发明的提供的基于模糊连接度的抠图样本集构造方法,能够被应用到多种涉及对已知像素采样的抠图系统中,在图像处理领域发挥重要作用。并且,由于视频本质上是一系列时间上连续的图像的组合,抠图是在视频编辑中的常用技术,本发明所提供基于模糊连接度的抠图样本集构造方法可以很自然地推广到视频处理领域。
本领域技术人员可以按照以上技术方案,采用计算机软件技术实现自动化样本集生成。为了便于实施参考,本发明提供了实施例细节说明:
实施例利用VC++8.0软件平台实现。对前景和背景像素均构建包含5个高斯分量的高斯混合模型,参数λ的设置则根据不同的输入图像而有所不同。高斯混合模型中所包含的高斯分量的个数及λ的值可以做适当调整,以期取得更好的效果。对输入图像,由用户对该图像进行类似于如图3所示的较精确的区域划分,或者进行类似于如图4所示用少量笔划提供的较粗糙的区域划分。以前景样本集的构造为例,总结前景样本集生成的基本流程,参见图1:输入图像后,根据交互信息将图像划分为前景、背景、未知区域;构造像素相对于已知前景的相似权重图;求解未知像素与前景边界的模糊连接度和最强路径;遍历树形结构,获取未知像素通过最强路径关联的前景边界根节点;以根节点为中心,搜索邻近已知前景像素,构造前景样本集。以前景样本集的更新为例,总结在用户增加新的前景笔划后对前景样本集进行更新的基本流程,参见图2:输入新的前景笔划后,判断新的前景笔划中是否包含了新的前景颜色,若是则重构像素相对于前景的相似权重图然后计算新的前景边界F”与已有的前景边界F’的模糊连接度FC(F’,F”),否则直接计算FC(F’,F”);然后根据FC(F’,F”)将新的未知像素集合划分为三个子集合,计算子集合2(即集合ΩU′ 2)中的未知像素与F”的模糊连接度和最强路径;判断计算出的模糊连接度值是否为0,是则该像素的模糊连接度值及前景样本集不变,否则该像素的前景模糊连接度及最强路径设置为重新计算出的值,然后根据新的最强路径关联的前景边界根节点更新该未知像素的前景样本集。
为了评估样本集构造结果的优劣,可以采用Robust抠图中基于样本对抠图参数估算的方法求解出输入图像的前景不透明度,该方法的性能表现强依赖于样本集的构造,因此,所提取出的抠图结果能够很好地反映出样本集构造的精确程度。而且,由于沿边界采样的方法在采样效率和精确度上明显优于局部窗口采样方法,可以只将本发明的方法与沿边界采样的方法进行比较。比较结果可以证明,沿边界采样方法容易错误地搜集到空间位置上距离未知像素近的错误样本,基于该方法所得的样本集估算出的前景不透明度中包含了大量显见的错误。相比之下,基于模糊连接度的采样方法沿到前景/背景边界的最强路径,搜索到与该未知像素在颜色空间最相似的前景/背景边界像素,然后以前景/背景边界像素为中心沿边界构造样本集,在此基础上提取出的前景不透明度精确完整。对远离前景/背景边界的未知像素,在给定相同的区域划分的情况下,基于沿边界采样的方法产生的抠图结果中会出现更多的错误,主要集中在背景或前景颜色较复杂,用户笔划没有充分覆盖到邻近已知颜色的未知区域中。而基于模糊连接度的方法则通过搜索未知像素到前景和背景边界的最强路径构造未知像素的样本集,而不是仅仅简单地将空间位置上距离近的已知像素作为样本,能够得到更高质量的抠图结果,在样本集构造的正确性与鲁棒性上,本发明的方法有显著改进。
本发明不限于实施例,还可以有许多变形。本领域的普通技术人员能从本发明公开的内容直接导出或联想到的所有变形,均应认为是本发明的保护范围。

Claims (9)

1.基于模糊连接度的抠图样本集构造方法,其特征在于:包括以下步骤,
步骤1,根据用户提供的区域划分,分别构造图像中像素相对于已知前景、已知背景的相似权重图Af、Ab
步骤2,根据相似权重图Af、Ab求解未知像素与已知的前景边界和背景边界的模糊连接度,并使用树形结构记录最强路径;
步骤3,遍历树形结构,获取每个未知像素通过最强路径关联的前景边界像素和背景边界像素,并分别记为前景边界根节点和背景边界根节点;
步骤4,以前景边界根节点为中心搜索邻近已知前景像素,构造每个未知像素对应的前景样本集;以背景边界根节点为中心搜索邻近已知背景像素,构造每个未知像素对应的背景样本集。
2.根据权利要求1所述的抠图样本集构造方法,其特征在于:像素相对于已知前景、已知背景的相似权重图Af、Ab按照如下步骤构造,
步骤1.1,根据用户通过交互信息提供的区域划分,将图像划分为三个区域:前景区域F,背景区域B和未知区域U,前景区域F即为图像的已知前景,背景区域B即为图像的已知背景,直接与未知区域U中的未知像素相邻接的前景/背景像素构成了前景/背景边界,位于前景/背景边界上的像素即为前景/背景边界像素;
步骤1.2,分别对已知前景的颜色和已知背景的颜色建模,定义图像中所有像素本身及邻接像素间相对于已知前景和背景的颜色相似度的权重函数,构造出相似权重图Af及Ab,定义Ao(p,q),o∈{f,b},表示相似权重图中像素p,q之间的边的权重值。
3.根据权利要求2所述的抠图样本集构造方法,其特征在于:所述模糊连接度定义如下:
(1)定义两个像素间路径的连接强度,
在图像内任意两个像素s和t之间存在着大量的路径,构成路径集合P,每条路径γ∈P由一系列连续邻接的节点p0,p1,p2,…,pn-1,pn连接而成,表示为γ(s→t)={p0,p1,p2,…,pn-1,pn},其中起始的节点p0为s,结束的节点pn为t,且 A o i = 0 &CenterDot; &CenterDot; &CenterDot; , n - 1 ( p i , p i + 1 ) > 0 , o∈{f,b},
Figure A2009100627310003C2
是节点pi和pi+1之间的边的相似权重,由相似权重图Af或Ab中取得;
路径γ(s→t)的连接强度记为
str o ( &gamma; ( s &RightArrow; t ) ) = min t = 0 &CenterDot; &CenterDot; &CenterDot; , n - 1 A o ( p i , p i + 1 ) , o &Element; { f , b }
(2)像素s和t的模糊连接度定义为两者间所有路径中最大的连接强度,记为
FC o ( s , t ) = max &gamma; &Element; P str o ( &gamma; ( s &RightArrow; t ) ) , o &Element; { f , b } ;
(3)图像中像素到前景边界或背景边界的模糊连接度,根据该像素与前景边界像素或背景边界像素的模糊连接度得到。
4.根据权利要求3所述的抠图样本集构造方法,其特征在于:步骤2具体实施方式如下,
(1)求解未知像素与前景边界的模糊连接度和最强路径,步骤如下,
步骤2.11,设置相关变量并初始化,包括有:
设置像素与前景边界间的模糊连接度FCf,初始化时将已知前景中所有的前景像素与前景边界间的FCf值设置为1,图像中所有非前景像素与前景边界间的FCf值设为0;
设置以模糊连接度FCf值为优先权的优先队列Qf,初始化时将所有前景边界像素加入到优先队列Qf中;
设置标识位Flagf,当Flagf值为“1”时表示对应像素可以被插入优先队列Qf中,为“0”则不能插入优先队列Qf中;初始化时前景边界像素与未知像素的Flagf值设为1,图像中剩余像素的Flagf值设为0;
设置树形结构Tf,用于记录前景边界像素和未知像素的父节点和孩子节点;初始化时把所有前景边界像素加入到树形结构Tf中,且对应的父节点设为-1,孩子节点设为空;并且把所有未知像素加入到树形结构Tf中,且对应的父节点和孩子节点设为空;
步骤2.12,当优先队列Qf不为空时,提取出优先队列Qf中具有最大的FCf值的像素p,继续执行步骤2.13、2.14;当Qf为空时,求解终止;
步骤2.13,对像素p的每个邻接像素q:
1)如果由相似权重图Af中取得的像素p,q之间的边的相似权重Af(p,q)>0,且Flagf(q)等于1,则执行后续操作;
2)比较像素p与前景边界的模糊连接度FCf(p)和Af(p,q)的大小;
3)令 FC min f = { FC f ( p ) , A f ( p , q ) } ;
FCf(q)为像素q与前景边界的模糊连接度,如果 FC min f > FC f ( q ) ) , 则执行后续操作;
4)设置 FC f ( q ) = FC min f ;
5)如果q在树形结构Tf中,则
更新q的父节点指针为p;
把q加为p的孩子节点;
如果q不在树形结构Tf中,则
把q插入到树形结构Tf中,q的父节点指针指向p;
把q加为p的孩子节点;
把节点q插入到队列Qf中;
步骤2.14,将p从Qf中删除,返回步骤2.12;
(2)求解未知像素与背景边界的模糊连接度和最强路径,步骤如下,
步骤2.21,设置相关变量并初始化,包括有:
设置像素与背景边界间的模糊连接度FCb,初始化时将已知背景中的所有背景像素与背景边界间的FCb值设置为1,图像中所有非背景像素与背景边界间的FCb值设为0;
设置以模糊连接度FCb值为优先权的优先队列Qb,初始化时将所有背景边界像素加入到优先队列Qb中;
设置标识位Flagb,当Flagb值为“1”时表示对应像素可以被插入优先队列Qb中,为“0”则不能插入优先队列Qb中;初始化时背景边界像素与未知像素的Flagb值设为1,图像中剩余像素的Flagb值设为0;
设置树形结构Tb,用于记录背景边界像素和未知像素的父节点和孩子节点;初始化时把所有背景边界像素加入到树形结构Tb中,且对应的父节点设为-1,孩子节点设为空;并且把所有未知像素加入到树形结构Tb中,且对应的父节点和孩子节点设为空;
步骤2.22,当优先队列Qb不为空时,提取出优先队列Qb中具有最大的FCf值的像素p,继续执行步骤2.23、2.24;当Qb为空时,求解终止;
步骤2.23,对像素p的每个邻接像素q:
1)如果由相似权重图Ab中取得的像素p,q之间的边的相似权重Ab(p,q)>0,且Flagb(q)等于1,则执行后续操作;
2)比较p与背景边界的模糊连接度FCb(p)和Ab(p,q)的大小;
3)令 FC min b = { FC b ( p ) , A b ( p , q ) } ,
FCb(q)为像素q与背景边界的模糊连接度,如果 FC min b > FC b ( q ) ) , 则执行后续操作;
4)设置FCb(q)=FCmin b
5)如果q在树形结构Tb中,则
更新q的父节点指针为p;
把q加为p的孩子节点;
如果q不在树形结构Tb中,则
把q插入到树形结构Tb中,q的父节点指针指向p;
把q加为p的孩子节点;
把节点q插入到队列Qb中;
步骤2.24,将p从Qb中删除,返回步骤2.22。
5.根据权利要求4所述的抠图样本集构造方法,其特征在于:用户通过交互提供新的笔划以补充区域划分时,只对与新的笔划所产生边界之间存在连通路径的未知像素重新处理,更新前景和背景样本集;
(1)当新的前景笔划加入后,具体处理方式如下,
1)如果新加入的前景笔划包含了新的前景颜色,则对新的已知前景的颜色重新建模,定义图像中所有像素本身及邻接像素间相对于已知前景的颜色相似度的权重函数,更新相似权重图Af,然后转向步骤2);如果新加入的前景笔划未包含新的前景颜色,直接转向步骤2);
2)计算新加入的前景笔划所产生前景边界F”与原有的前景边界F’间的模糊连接度值FC(F’,F”);
3)根据FC(F’,F”)将新的前景笔划加入后所产生新的未知区域U’划分为三个子集合:
&Omega; U &prime; 1 = { i | i &Element; U &prime; , FC f ( F &prime; , i ) > FC ( F &prime; , F &prime; &prime; ) }
&Omega; U &prime; 2 = { i | i &Element; U &prime; , FC f ( F &prime; , i ) > FC ( F &prime; , F &prime; &prime; ) }
&Omega; U &prime; 3 = { i | i &Element; U &prime; , FC f ( F &prime; , i ) > FC ( F &prime; , F &prime; &prime; ) }
其中FCf(F′,i)表示新的未知区域U’中任意像素i与原有的前景边界F’的模糊连接度;
4)包含在集合ΩU′ 1和ΩU′ 3中的未知像素,其对应的FCf及前景样本集仍使用原值;对包含在集合ΩU′ 2中的未知像素则按照步骤5)做进一步处理;
5)如果集合ΩU′ 2不为空,则:
首先,初始化相关变量:设置被新的前景笔划所覆盖的未知像素的FCf为1,ΩU′ 2中未知像素的FCf为0;设置ΩU′ 2中未知像素的Flagf位为1,图像中剩余像素的Flagf位为0;把新的前景边界F”中的像素加入到队列Qf中;
然后,重新计算ΩU′ 2中未知像素与新的前景边界F”间的FCf和最强路径;
最后,根据计算出的结果做如下处理:如果未知像素的FCf值为0,则该未知像素的FCf及前景样本集仍使用原值;否则,该未知像素与F”之间通过区域ΩU′ 2是连通的,其FCf设置为重新计算出的FCf值,并根据新的最强路径更新该像素的前景样本集;
(2)当新的背景笔划加入后,具体处理方式如下,
1)如果新加入的背景笔划包含了新的背景颜色,则对新的已知背景的颜色重新建模,定义图像中所有像素本身及邻接像素间相对于已知背景的颜色相似度的权重函数,更新相似权重图Ab,然后转向步骤2);如果新加入的背景笔划未包含新的背景颜色,直接转向步骤2);
2)计算新加入的背景笔划所产生背景边界B”与原有的背景边界B’间的模糊连接度值FC(B’,B”);
3)根据FC(B’,B”)将新的背景笔划加入后所产生新的未知区域U’划分为三个子集合:
&Omega; U &prime; 1 = { i | i &Element; U &prime; , FC b ( B &prime; , i ) > FC ( B &prime; , B &prime; &prime; ) }
&Omega; U &prime; 2 = { i | i &Element; U &prime; , FC b ( B &prime; , i ) > FC ( B &prime; , B &prime; &prime; ) }
&Omega; U &prime; 3 = { i | i &Element; U &prime; , FC b ( B &prime; , i ) > FC ( B &prime; , B &prime; &prime; ) }
其中FCb(B′,i)表示新的未知区域U’中任意像素i与原有的背景边界B’的模糊连接度;
4)包含在集合ΩU′ 1和ΩU′ 3中的未知像素,其对应的FCb及背景样本集仍使用原值;对包含在集合ΩU′ 2中的未知像素则按照步骤5)做进一步处理;
5)如果集合ΩU′ 2不为空,则:
首先,初始化相关变量:设置被新的背景笔划所覆盖的未知像素的FCb为1,ΩU′ 2中未知像素的FCb为0;设置ΩU′ 2中未知像素的Flagb位为1,图像中剩余像素的Flagb位为0;把新的背景边界B”中的像素加入到队列Qb中;
然后,重新计算ΩU′ 2中未知像素与新的背景边界B”间的FCb和最强路径;
最后,根据计算出的结果做如下处理:如果未知像素的FCb值为0,则该未知像素的FCb及背景样本集仍使用原值;否则,该未知像素与F”之间通过区域ΩU′ 2是连通的,其FCb设置为重新计算出的FCb值,并根据新的最强路径更新该像素的背景样本集。
6.根据权利要求5所述的抠图样本集构造方法,其特征在于:
(1)计算新加入的前景笔划所产生的前景边界F”与原有的前景边界F’间的模糊连接度FC(F’,F”),具体实现方式如下
前景边界F”与原有的前景边界F’间的模糊连接度为两个像素集合间的模糊连接度,定义为:
FC ( F &prime; , F &prime; &prime; ) = max f &prime; &Element; F &prime; max f &prime; &prime; &Element; F &prime; &prime; FC ( f &prime; , f &prime; &prime; )
考察F”中每个像素i与F’的模糊连接度FCf(F’,i),并从中挑选出最大的FCf(F’,i)值作为FC(F’,F”),即
FC ( F &prime; , F &prime; &prime; ) = max i &Element; F &prime; &prime; FC f ( F &prime; , i )
(2)计算新加入的背景笔划所产生的背景边界B”与原有的背景边界B’间的模糊连接度FC(B’,B”),具体实现方式如下
背景边界B”与原有的背景边界B’间的模糊连接度为两个像素集合间的模糊连接度,定义为:
FC ( B &prime; , B &prime; &prime; ) = max b &prime; &Element; B &prime; max b &prime; &prime; &Element; B &prime; &prime; FC ( b &prime; , b &prime; &prime; )
考察B”中每个像素i与B’的模糊连接度FCb(B’,i),并从中挑选出最大的FCb(B’,i)值作为FC(B’,B”),即
FC ( B &prime; , B &prime; &prime; ) = max i &Element; B &prime; &prime; FC b ( B &prime; , i ) .
7.根据权利要求4或5或6所述的抠图样本集构造方法,其特征在于:步骤3具体实施方式为,遍历树形结构,通过孩子节点或父节点指针追溯未知像素到前景或背景边界的最强路径,获取最强路径所关联的前景或背景边界像素。
8.根据权利要求1或2或3或4或5或6所述的抠图样本集构造方法,其特征在于:步骤4具体实施方式为,获取某个未知像素对应的前景边界根节点后,以该根节点为中心沿前景边界搜索m个前景像素,所得的m个前景像素即前景样本构成了前景样本集,m为用户指定的每个未知像素的前景样本的个数;获取某个未知像素对应的背景边界根节点后,以该根节点为中心沿背景边界搜索n个背景像素,所得的n个背景像素即背景样本构成了背景样本集,其中n为用户指定的每个未知像素的背景样本的个数。
9.根据权利要求1或2或3或4或5或6所述的抠图样本集构造方法,其特征是所述抠图样本集构造方法在图像处理领域或视频处理领域中的应用。
CN2009100627318A 2009-06-16 2009-06-16 基于模糊连接度的抠图样本集构造方法 Expired - Fee Related CN101582168B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100627318A CN101582168B (zh) 2009-06-16 2009-06-16 基于模糊连接度的抠图样本集构造方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100627318A CN101582168B (zh) 2009-06-16 2009-06-16 基于模糊连接度的抠图样本集构造方法

Publications (2)

Publication Number Publication Date
CN101582168A true CN101582168A (zh) 2009-11-18
CN101582168B CN101582168B (zh) 2011-06-15

Family

ID=41364303

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100627318A Expired - Fee Related CN101582168B (zh) 2009-06-16 2009-06-16 基于模糊连接度的抠图样本集构造方法

Country Status (1)

Country Link
CN (1) CN101582168B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101989353A (zh) * 2010-12-10 2011-03-23 中国科学院深圳先进技术研究院 图像抠图方法
CN103206957A (zh) * 2012-12-26 2013-07-17 辽宁师范大学 车辆自主导航的车道线检测与跟踪方法
CN104468831A (zh) * 2014-12-26 2015-03-25 山东中创软件商用中间件股份有限公司 一种基于网络拓扑结构的消息服务网络集群
WO2018141232A1 (zh) * 2017-02-06 2018-08-09 腾讯科技(深圳)有限公司 图像处理方法、计算机存储介质及计算机设备
CN109325946A (zh) * 2018-09-14 2019-02-12 北京石油化工学院 一种危险化学品堆垛监测方法和系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1224242C (zh) * 2003-04-14 2005-10-19 浙江大学 数字图像处理中以任意图像为背景的快速抠图方法
CN1260682C (zh) * 2003-05-18 2006-06-21 浙江大学 数字图像处理中基于hvs的自然图像抠图方法
CN1564198A (zh) * 2004-04-13 2005-01-12 浙江大学 基于感知颜色空间的自然图像抠图方法
CN100539698C (zh) * 2007-07-10 2009-09-09 浙江大学 一种数字视频处理中交互式时空一致的视频抠图方法
CN100565584C (zh) * 2008-03-25 2009-12-02 武汉大学 一种具有纠偏性的自然图像抠图的全局优化方法

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101989353A (zh) * 2010-12-10 2011-03-23 中国科学院深圳先进技术研究院 图像抠图方法
CN101989353B (zh) * 2010-12-10 2012-10-17 中国科学院深圳先进技术研究院 图像抠图方法
CN103206957A (zh) * 2012-12-26 2013-07-17 辽宁师范大学 车辆自主导航的车道线检测与跟踪方法
CN103206957B (zh) * 2012-12-26 2016-03-02 辽宁师范大学 车辆自主导航的车道线检测与跟踪方法
CN104468831A (zh) * 2014-12-26 2015-03-25 山东中创软件商用中间件股份有限公司 一种基于网络拓扑结构的消息服务网络集群
WO2018141232A1 (zh) * 2017-02-06 2018-08-09 腾讯科技(深圳)有限公司 图像处理方法、计算机存储介质及计算机设备
CN109325946A (zh) * 2018-09-14 2019-02-12 北京石油化工学院 一种危险化学品堆垛监测方法和系统
CN109325946B (zh) * 2018-09-14 2021-08-24 北京石油化工学院 一种危险化学品堆垛监测方法和系统

Also Published As

Publication number Publication date
CN101582168B (zh) 2011-06-15

Similar Documents

Publication Publication Date Title
CN104732506A (zh) 一种基于人脸语义分析的人物照片颜色风格转换方法
CN110335277A (zh) 图像处理方法、装置、计算机可读存储介质和计算机设备
Yang et al. User-friendly interactive image segmentation through unified combinatorial user inputs
CN103942794B (zh) 一种基于置信度的图像协同抠图方法
CN102034247B (zh) 一种基于背景建模对双目视觉图像的运动捕捉方法
CN101582168B (zh) 基于模糊连接度的抠图样本集构造方法
CN107146201A (zh) 一种基于改进图像融合的图像拼接方法
US20070003154A1 (en) Video object cut and paste
CN100448271C (zh) 基于全景图拼接的视频编辑方法
CN107862698A (zh) 基于k均值聚类的光场前景分割方法及装置
CN107657625A (zh) 融合时空多特征表示的无监督视频分割方法
CN103400386A (zh) 一种用于视频中的交互式图像处理方法
CN106056606A (zh) 图像处理方法及装置
CN106952276A (zh) 一种图像抠图方法及装置
CN105957078A (zh) 一种基于图割的多视点视频分割方法
CN108564528A (zh) 一种基于显著性检测的肖像照片自动背景虚化方法
CN103578107B (zh) 一种交互式图像分割方法
CN109410243A (zh) 基于区域交叠与模板匹配的多目标跟踪方法
Zhang et al. GAIN: Gradient augmented inpainting network for irregular holes
CN102254325A (zh) 一种运动模糊场景的分割及前景提取方法和系统
CN107909079A (zh) 一种协同显著性检测方法
JP2015041254A (ja) Trimap生成装置、Trimap生成方法、及びプログラム
CN109887008A (zh) 基于前后向平滑和o(1)复杂度视差立体匹配方法、装置和设备
Kataoka et al. Ten-million-order human database for world-wide fashion culture analysis
CN106447718A (zh) 一种2d转3d深度估计方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
ASS Succession or assignment of patent right

Owner name: CHANGSHU ZIJIN INTELLECTUAL PROPERTY SERVICE CO.,

Free format text: FORMER OWNER: WUHAN UNIVERSITY

Effective date: 20121211

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 430072 WUHAN, HUBEI PROVINCE TO: 215500 SUZHOU, JIANGSU PROVINCE

TR01 Transfer of patent right

Effective date of registration: 20121211

Address after: 215500 Changshou City South East Economic Development Zone, Jiangsu, Jin Road, No. 8

Patentee after: Changshu Zijin Intellectual Property Service Co.,Ltd.

Address before: 430072 Hubei city of Wuhan province Wuchang Luojiashan

Patentee before: Wuhan University

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20160316

Address after: 215500 No. 8, Jin Du Road, Changshou City hi tech Industrial Development Zone, Jiangsu, China

Patentee after: Changshu Nanjing Normal University Development Research Institute Co.,Ltd.

Address before: 215500 Changshou City South East Economic Development Zone, Jiangsu, Jin Road, No. 8

Patentee before: Changshu Zijin Intellectual Property Service Co.,Ltd.

TR01 Transfer of patent right

Effective date of registration: 20201222

Address after: Building 3, No.88, Xianshi Road, high tech Zone, Changshu City, Suzhou City, Jiangsu Province

Patentee after: Changshu Zijin Intellectual Property Service Co.,Ltd.

Address before: 215500 No.8, Jindu Road, Changshu hi tech Industrial Development Zone, Jiangsu Province

Patentee before: CHANGSHU NANJING NORMAL UNIVERSITY DEVELOPMENT RESEARCH INSTITUTE Co.,Ltd.

TR01 Transfer of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110615