CN103136726A - 用于恢复图像的深度信息的方法和设备 - Google Patents

用于恢复图像的深度信息的方法和设备 Download PDF

Info

Publication number
CN103136726A
CN103136726A CN201210464952XA CN201210464952A CN103136726A CN 103136726 A CN103136726 A CN 103136726A CN 201210464952X A CN201210464952X A CN 201210464952XA CN 201210464952 A CN201210464952 A CN 201210464952A CN 103136726 A CN103136726 A CN 103136726A
Authority
CN
China
Prior art keywords
view
binocular parallax
pixel
medial view
medial
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
CN201210464952XA
Other languages
English (en)
Other versions
CN103136726B (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.)
Samsung Electronics Co Ltd
Original Assignee
Samsung Electronics Co Ltd
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 Samsung Electronics Co Ltd filed Critical Samsung Electronics Co Ltd
Publication of CN103136726A publication Critical patent/CN103136726A/zh
Application granted granted Critical
Publication of CN103136726B publication Critical patent/CN103136726B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • G06T7/55Depth or shape recovery from multiple images
    • G06T7/593Depth or shape recovery from multiple images from stereo images
    • G06T7/596Depth or shape recovery from multiple images from stereo images from three or more stereo images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/194Segmentation; Edge detection involving foreground-background segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/20Analysis of motion
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/20Analysis of motion
    • G06T7/215Motion-based segmentation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N13/00Stereoscopic video systems; Multi-view video systems; Details thereof
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10016Video; Image sequence
    • G06T2207/10021Stereoscopic video; Stereoscopic image sequence

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Computer Graphics (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

图像处理设备和方法可通过使用多个视图的图像来估计所述多个视图中的中间视图的双目视差图。所述图像处理设备可基于帧的双目视差图来从中间视图中检测运动对象。通过检测运动对象,中间视图中的像素可被分为动态像素和静态像素。所述图像处理设备可将束优化和基于局部三维(3D)线模型的时间优化应用到中间视图以增强静态像素和动态像素的双目视差值。

Description

用于恢复图像的深度信息的方法和设备
本申请要求于2011年11月30日在韩国知识产权局提交的第10-2011-0127153号韩国专利申请的优先权益,该申请的公开通过引用合并于此。
技术领域
这里公开的示例实施例涉及一种用于恢复图像的深度信息的方法和设备,更具体地说,涉及一种用于基于由多个相机捕捉的视频序列,恢复可在时空上一致的一组深度图的方法和设备。
背景技术
在计算机视觉领域中,基于多个图像或视频来恢复深度图像是基本问题。恢复深度图像可被用于三维(3D)重建、3D电影制作、自主机器人导航、视频分割和编辑等。
通过立体匹配,可基于相机参数已知的至少两个的多图像来恢复深度信息。
最近用于恢复深度信息的方法可被划分为两类。第一,若干种方法可使用局部的基于窗口的并当估计视差时可采用局部的赢家通吃所有(Winner-Takes-All)的策略的方案。相反,其它方法可将深度的估计表示为能量最小化问题,并可使用图像切割和置信传播算法以解决能量最小化的问题。
恢复的高质量深度信息可被用于3D视频、立体视频分析、视频编辑和增强现实等。
发明内容
通过提供一种处理图像的方法来实现上述和/或其它方面,所述方法包括:针对一个或多个帧中的每一帧,使用包括在所述帧中的每一帧中的左图像、中间图像和右图像来估计中间视图的双目视差图,所述中间视图包括中间图像和估计的双目视差图;通过使用所述帧的中间视图的深度图来估计指示中间视图的背景的背景视图;基于中间视图和背景视图之差,检测中间视图中的运动对象。
左图像和右图像中的至少一个可相应于多个图像。
估计双目视差图的步骤可包括:当左图像、中间图像和右图像被提供时,通过使用反向传播(BP)算法重复使能量函数的值最小化来计算中间视图的双目视差图,基于与中间视图有关的遮蔽信息和中间视图的双目视差图来确定能量函数的值;通过融合将不同的平面拟合应用到中间视图的双目视差图而获得的结果来增强中间视图的双目视差图。
计算的步骤可包括:基于中间视图的双目视差图来估计遮蔽信息,基于遮蔽信息来估计中间视图的双目视差图。
所述方法还可包括:使用与所述帧中的光流有关的信息来增强中间视图的双目视差图。
估计背景视图的步骤可包括:对所述帧中的目标帧的邻近帧的中间视图进行弯曲以对目标帧的中间视图进行匹配,以产生包括在目标帧的中间视图中的像素的候选双目视差值和候选颜色值;将中值滤波器应用到像素的候选双目视差值和候选颜色值以产生目标帧的背景视图的图像和双目视差图。
检测运动对象的步骤可包括:针对每个像素比较所述帧中的每一帧的中间视图和背景视图,以测量包括在中间视图和背景视图中的每个像素之间的颜色、深度和运动深度的差;基于包括在中间视图中的每个像素的颜色、深度和运动深度的差中的至少一个,从中间视图分割出动态区域。
分割出的动态区域可相应于运动对象。
检测的步骤还可包括:通过使用与所述一个或多个帧中的光流有关的信息来增强分割出的动态区域中的时空一致性。
所述方法还可包括:将束优化和基于局部三维(3D)线模型的时间优化应用到所述帧的中间视图,以增强包括在中间视图中的静态像素和动态像素的双目视差值。
动态像素可相应于与运动对象对应的像素,静态像素可相应于与运动对象不对应的像素。
所述方法还可包括:通过对所述帧的中间视图进行弯曲来产生所述帧的左视图的双目视差图和所述帧的右视图的双目视差图;基于分割,填充包括在左视图的双目视差图和右视图的双目视差图中的孔;将时空优化应用到左视图的双目视差图和右视图的双目视差图。
所述孔可由中间视图中的遮蔽导致的丢失像素而产生。
通过提供一种用于处理图像的设备来实现上述和/或其它方面,所述设备包括:双目视差图估计器,使用包括在一个或多个帧中的每一帧中的左图像、中间图像和右图像来估计所述帧中的每一帧的中间视图的双目视差图,所述中间视图包括中间图像和估计的双目视差图;背景视图估计器,通过使用所述帧的中间视图的深度图来估计指示中间视图的背景的背景视图;运动对象检测器,基于中间视图和背景视图之差,检测中间视图中的运动对象。
实施例的另外方面将在下面的描述中部分地阐明,并且从描述中部分是清楚的,或者可通过本公开的实践可以被理解。
附图说明
专利或申请文件包含至少一个彩色的附图。在请求时以及支付必要的费用时将由我方提供带有彩色附图的本专利或专利申请公开的副本。
从以下结合附图进行的对实施例的描述中,这些和/或其它方面将变得明显并且更易于理解,其中:
图1示出根据示例实施例的图像处理设备;
图2示出根据示例实施例的图像处理方法;
图3示出根据示例实施例的估计双目视差图的方法;
图4示出根据示例实施例的估计背景视图的方法;
图5示出根据示例实施例的检测运动对象的方法;
图6示出根据示例实施例的时空立体匹配的方法;
图7示出根据示例实施例的产生左右双目视差图的方法;
图8示出根据示例实施例的时间融合的结果;
图9示出根据示例实施例的三种差异图和深度误差区域;
图10示出根据示例实施例的从视图分割的运动对象;
图11示出根据示例实施例的执行了时空优化之后的增强的双目视差图。
具体实施方式
现在将详细参考实施例,其示例在附图中示出,其中,相同的标号始终表示相同的元件。以下通过参照附图来描述实施例以解释本公开。
图1示出根据示例实施例的图像处理设备100的框图。
图像处理设备100可包括双目视差图估计器110、背景视图估计器120、运动对象检测器130、时空立体匹配器140和左右双目视差图产生器150。
双目视差图估计器100可包括初始双目视差图计算器111、平面拟合应用器112和光流应用器113。
背景视图估计器120可包括邻近帧弯曲单元121和中值滤波单元122。
运动对象检测器130可包括差测量单元131、动态区域分割单元132和光流应用器133。
时空立体单元140可包括静态像素优化器141、轨迹追踪器142和动态像素优化器143。
左右双目视差图产生器150可包括左右双目视差图初始化单元151、基于分割的孔填充单元152和左右双目视差图时空优化器153。
以下,将描述所述元件的操作和功能。
图2示出根据示例实施例的图像处理方法。
在操作210,双目视差图估计器110可通过使用包括在一个或多个帧中的每一帧中的左图像、中间图像和右图像来估计所述一个或多个帧中的每一帧的中间视图的双目视差图。
中间视图可包括中间图像和估计的双目视差图。
双目视差图估计器100可通过使用每个帧的三目视图来独立地估计每个帧的深度图。这里,每个帧的深度图可相应用于与三目视图中的中间视图的图像相应的深度图。
在操作220,背景视图估计器120可通过使用所述帧的中间视图的深度图来估计指示中间视图的背景的背景视图。
在操作230,运动对象检测器130可基于中间视图和背景视图之间的差来从所述帧的中间视图中检测运动对象。
运动对象检测器130可通过检测运动对象来将中间视图的序列中的静态像素与动态像素分离。动态像素可以是与运动对象相应的像素,静态像素可以是不与运动对象相应的像素。
为了处理无纹理的区域,运动对象检测器130可使用基于分割的算法。可基于具有相似颜色的邻近像素通常具有相似深度值的一般观念来操作基于分割的算法,并可由3D平面来表示通过分割而分离的分段。
在操作240,时空立体匹配器140可将束优化(bundle optimization)和基于局部3D线模型的时空优化应用到所述帧的中间视图,从而增强中间视图中的静态像素和动态像素的双目视差值。
在中间视图的双目视差图被恢复之后,可产生左视图的双目视差图和右视图的双目视差图。
在操作250,左右双目视差图产生器150可基于产生的中间视图的双目视差图来产生左视图的双目视差图和右视图的双目视差图。
图3示出根据示例实施例的估计双目视差图的方法。
在与双目视差图的估计相应的操作210中,可通过例如三目立体匹配来初始化双目视差图。操作210可包括操作310、320和330。
以下,IL、IM和IR分别表示由三目相机捕捉的左图像、中间图像和右图像。这里,由三目相机产生或捕捉的三个图像IL、IM和IR可仅为示例。本实施例可应用于包括图像的帧。例如,IM表示多个图像中的中间图像。IL表示在IM的左侧的多个图像。IR表示在IM的右侧的多个图像。
可由O={OL,OR}来表示中间视图的遮蔽信息(Occlusion information)。OL(x)=1可指示包括在中间图像中的像素x从IL不可见。OR(x)=1可指示像素x在IR中不可见。
Ii(x)表示图像Ii中的像素x的颜色。i可以是L、M和R之一。
D表示中间视图的双目视差图。双目视差可被定义为反深度。当像素x的深度为z时,像素x的双目视差值(D(x))可以是1/z。可由[dmin,dmax]来表示双目视差的范围。
双目视差可被量化为m+1个级。可由等式1来定义量化的双目视差的第k级(dk)。
[等式1]
dk=(m-k)/m·dmin+k/m·dmax
在等式1中,dmin表示量化的双目视差的最小级,dmax表示量化的双目视差的最大级。通常,m可被设置为在50到80之间的范围内。
因此,立体匹配可被表示为标签问题(labeling problem)。
在操作310,初始双目视差图计算器111可通过使用反向传播(BP)算法或融合运动算法来重复最小化能量函数的值,从而计算中间视图的双目视差图。也就是说,初始双目视差图计算器111可最小化图2的能量函数E(D,O;I),从而得到中间视图的双目视差图D。在操作310首先产生的中间视图的双目视差图可被称为中间视图的初始双目视差图。
[等式2]
E(D,O;I)=Ed(D,O;I)+Es(D,O;I)+Ev(D,O;I)
在等式2中,I表示图像,即,IL、IM和IR。E表示能量函数,Ed表示数据项,Es表示平滑项,Ev表示可见度项。Ed可针对给定的I来测量双目视差D的合适性。为了解多标签等式,初始双目视差图计算器111可使用置信传播(belief propagation)或图形切割。
如等式2所示,可基于与中间视图相关的遮蔽信息O以及当IL、IM和IR被提供时中间视图的双目视差图D,确定能量函数E(D,O;I)。
可由等式3来定义Ed
[等式3]
E d ( D , O ; I ) = Σ x ( 1 - O L ( x ) ) ( 1 + O R ( x ) ) ρ d ( F ( x , d x ; I L , I M ) ) +
Σ x ( 1 - O R ( x ) ) ( 1 + O L ( x ) ) ρ d ( F ( x , d x ; I R , I M ) ) +
Σ x O L ( x ) O R ( x ) η
在等式3中,F表示价值函数。η表示惩罚参数。
ρd可以是抗噪声的截断LI正态函数。可由等式4来定义ρd
[等式4]
ρd=-ln((1-ed)exp(-|x|/δ)+ed)
在等式4中,ed和δ可以是用于控制ρd的形状的参数。
如等式5所示,F,即价值函数,可被定义为色差价值的方差和(SSD)。
[等式5]
F(x,d,I)=||I1(x)-I2(I1->2(x,d))||2
I1表示投影之前的图像,即,IM。I2表示投影之后的图像,即,IL或IR。I1(x)表示图像Ii的像素x的颜色值。I1->2表示中间图像中的像素x以双目视差d投影在左图像或右图像的像素。
由等式6来定义Es
E s ( D ) = Σ s , t ∈ C λ ( s , t ) ρ s ( d s , d t )
在等式6中,s,t∈C指示像素s和像素t邻近。也就是说,(s,t)表示中间帧中的一对邻近像素。C表示集。
可如等式7所示来选择关于ds和dt的价值(ρd)。
[等式7]
ρs(ds,dt)=min{|ds-dt|,T}
在等式7中,T表示用于控制价值的限度的参数。
可由等式8来定义关于像素x和像素t的自适应平滑权重(λ(s,t))。
[等式8]
λ ( s , t ) = β u ( s ) | | I ( s ) - I ( t ) + ϵ | |
在等式8中,β表示用于控制平滑项的权重的参数。ε表示用于控制对于色差的敏感度的参数。u(s)表示归一化因子。可由等式9来定义u(s)。
[等式9]
u ( s ) = | N ( s ) | / Σ t ∈ N ( s ) 1 | | I ( s ) - I ( t ) | | + ϵ
在等式9,N(s)表示与像素s邻近的一组像素Ev
为了执行可见度一致性约束,可由等式10来定义Ev
[等式10]
E v ( D , O ; I ) = Σ s β ω | O ( s ) - W ( s ; D ) | + Σ s , t ∈ N β o | O ( s ) - O ( t ) |
在等式10,βω、w(x)和βo可以是用于计算Ev的参数。βω可以是指示针对虽然像素不可见但该像素被设置为无遮蔽的情况的惩罚的参数。βo可以是指示针对在邻近像素中一些像素被遮蔽而其余像素不被遮蔽的情况的惩罚的参数。
这里,将被计算的双目视差图可仅为中间视图的双目视差图。将左视图和中间视图视为示例,可由等式11来定义可见度项。
[等式11]
E v ( D , O L ; I L , I M ) = Σ s β ω | O L ( s ) - W L ( s ; D ) | + Σ s , t ∈ N β o | O L ( s ) - O L ( t ) |
在等式11中,OL(x)表示与左视图有关的遮蔽信息。OL(x)=1指示在左视图中像素x不可见。OL(x)=0指示在左视图中像素x可见。WL(s;D)表示初始遮蔽信息。可基于等式12来获得WL(s;D)。
[等式12]
Figure BDA00002415534600082
等式12被用于基于分割的方法以处理遮蔽。等式12显示像素s被遮蔽。当图像I中的像素t满足以下条件时像素s可被像素t遮蔽:像素t和像素s存在于不同的平面上,在不同的视图(例如,当I=IM时的IL和IR时)中与像素s相应的点和与像素t相应的点相同,并且ds<dt
为了获得更好的结果,初始双目视差图计算器111可在处理遮蔽时使用重复计算可见度图和双目视差图的值的算法。该算法可以是精确地估计双目视差信息以及减轻遮蔽的效果的有效方法。
如上所述,初始双目视差图计算器111可通过迭代优化处理来获得能量函数E(D,O;I)的值。与迭代优化处理相应的操作310可包括以下操作:当双目视差图被提供时,初始双目视差图计算器111基于双目视差图估计遮蔽信息;当遮蔽信息被提供时,初始双目视差图计算器111基于遮蔽信息估计双目视差图。
首先,将描述当双目视差图被提供时估计遮蔽的操作。
当当前估计的双目视差图D被提供时,可由等式13来定义能量函数EO。EO可以是等式2的E(D,O;I)。
[等式13]
E O = Σ x ( 1 - O L ( x ) ) ( 1 + O R ( x ) ) ρ d ( F ( x , d x ; I L , I M ) ) +
Σ x ( 1 - O R ( x ) ) ( 1 + O L ( x ) ) ρ d ( F ( x , d x ; I R , I M ) ) + Σ x O L ( x ) O R ( x ) η +
Σ x β ω ( | O L ( x ) - W L ( x ; D ) | + | O R ( x ) - W R ( x ; D ) | ) +
Σ s , t β o ( | O L ( s ) - O L ( t ) | + | O R ( s ) - O R ( t ) | )
等式13可包括两个变量(即,OL和OR)以及一个二次项OL(x)OR(x)。l因此,EO的值可能不容易被获得。
可分别由等式14和等式15来定义O和W。
[等式14]
O ( p ) = 0 , O L ( p ) = 0 , O R ( p ) = 0 1 , O L ( p ) = 0 , O R ( p ) = 1 2 , O L ( p ) = 1 , O R ( p ) = 0 3 , O L ( p ) = 1 , O R ( p ) = 1
[等式15]
W ( p ) = 0 , W L ( p ) = 0 , W R ( p ) = 0 1 , W L ( p ) = 0 , W R ( p ) = 1 2 , W L ( p ) = 1 , W R ( p ) = 0 3 , W L ( p ) = 1 , W R ( p ) = 1
可由等式16来定义EO
[等式16]
E o = Σ x ( T ( O ( x ) = 0 ) E 0 ( x ; D ) + T ( O ( x ) = 1 ) E 1 ( x ; D ) ) +
Σ x ( T ( O ( x ) = 2 ) E 2 ( x ; D ) + T ( O ( x ) = 3 ) E 3 ( x ; D ) ) +
Σ x β ω ( ψ ( O ( x ) , W ( x ) ) ) + Σ s , t β o ( ψ ( O ( s ) , O ( t ) ) )
在等式16,E0、E1、E2和E3可以是分别针对四种情况的能量函数。T可以是当测试条件为真时返回值1的测试函数。可由等式17来定义ψ(a,b)。
[等式17]
ψ(a,b)=Ceil(|a-b|/2)
初始双目视差图计算器111可通过使用BP算法或融合运动算法获得EO的值。
接下来,将描述当遮蔽信息被提供时估计双目视差的操作。
当当前估计的遮蔽O被提供时,可由等式18来定义能量函数ED。ED可相应于图2的E(D,O;I)。
[等式18]
E D = E d ( D ; O , I ) +
Σ x β ω ( O L ′ ( P M → L ( x , D ( x ) ) ) + O R ′ ( P M → R ( x , D ( x ) ) ) ) +
Σ s , t λ ( s , t ) ρ s ( D ( s ) , D ( t ) )
在等式18中,第一项中的Ed(D;O,I)指示数据价值。
第二项指示可见度项,PM→L(x,D(x))指示投影函数。PM→L(x,D(x))可基于像素x的位置和双目视差值来将中间视图中的像素x投影到左视图上。OL′表示左视图的遮蔽信息。初始双目视差图计算器111可通过基于先前计算的双目视差D弯曲中间视图中的所有像素来获得OL′。然而,OL′是不平的,因此,初始双目视差图计算器111可通过应用中值滤波器来使OL′平滑。
第三项表示上面描述的平滑项。
初始双目视差图计算器111可基于BP算法来计算双目视差图D。
在操作320,平面拟合应用器112可将平面拟合应用到计算的双目视差图D以增强双目视差图。
当在平面拟合中分割无法被完全执行时,计算的双目视差图D可包括一些误差。因此,平面拟合应用器112可使用不同的分割结果来获得平面拟合的各种结果。此外,平面拟合应用器112可将平面拟合的结果融合,从而可获得更高质量的双目视差图。
平面拟合应用器112可融合通过将不同的平面拟合应用到在操作310计算的中间视图的双目视差图D而产生的结果,从而增强中间视图的双目视差图。平面拟合应用器112可融合具有不同的平均位移参数的增强后的双目视差图,从而获得更精确的双目视差图。通过融合,平面拟合应用器112可有效地增强在先前阶段计算的初始双目视差图,特别是针对无纹理区域进行增强。
当k个不同估计的双目视差图{D1′、D2′、...、和Dk′}存在,并且Di=round(Di′)时,k个双目视差建议值{D1、D2、...、和Dk}可存在。因此,每个像素的双目视差值(例如,像素x的双目视差值)可具有k个候选值{D1(x)、D2(x)、...、和Dk(x)}。平面拟合应用器112可将Dk+1(x)=∑iDi(x)/k添加为双目视差值的候选值。在添加之后,k+1个建议值{D1、D2、...、和Dk+1}可存在。
可由等式19来定义能量函数E(D,O;I)。
[等式19]
E ( D , O ; I ) = Σ x ρ ( x , d x ; O , I ) + Σ s , t ∈ N λ ( s , t ) ρ s ( d s , d 1 )
在等式19中,s,t∈N指示像素s和像素t彼此邻近。λ(s,t)表示关于像素s和像素t的自适应平滑权重。
每个像素(例如,像素x)可具有k+1个侯选双目视差值。因此,由平面拟合应用器112估计每个像素的双目视差图可相应于包括k+1个标签的标签问题。平面拟合应用器112可使用平均位移算法来产生不同的双目视差建议值。平面拟合应用器112可使用将不同的参数应用到双目视差建议值的平均位移算法。所述参数可被用于图像分割。
在操作330,光流应用器113可基于一个或多个帧中与光流有关的信息来增强中间视图的双目视差图D。
在一个或多个场景中可存在动态对象。因此,两个不同的帧中的两个点之间的关系可由这两个帧的光流来测量。
为了计算从第k帧到第k+n帧的光流,光流应用器113可使用如下所述的方法。
从第k+i帧到第k+i+1帧的光流相应于Uk+i=(uk+i,vk+i),从第k帧到第k+n帧的光流可以是所描述的方法可被用于计算从第k+n帧到第k帧的反向光流。
通常,尤其是对于被遮蔽的区域或无遮蔽的区域来说,光流是不可靠的。因此,为了减少由于光流误差而包括的现象,光流应用器113可测量光流的可靠性。当x相应于第k帧的像素时,光流应用器113可通过使用光流Uk→k+n来识别第k+n帧中与x相应的点y。此外,光流应用器113可使用光流Uk+n→k,从而将点y映射到第k帧。与点y相应的点可被称为x′。当光流精确时,x和x′可以是相同的点。当光流不精确时,x和x′可相距一定距离。因此,光流应用器113可基于等式20来计算光流的可靠性(Popt)。
[等式20]
P opt = exp ( - | | x - x ′ | | 2 σ )
在等式20中,σ表示用于控制Popt的形状的参数。
光流应用器113可使用光流将第k-1帧的双目视差图和第k+1帧的双目视差图投影在第k帧上。因此,第k帧可具有三个双目视差图
Figure BDA00002415534600123
光流应用器113可基于等式21来针对像素x(即,第k帧的像素)来计算像素x的双目视差值为d的概率P(d|x)。
[等式21]
P ( d | x ) = Σ i = 1 3 Σ y ∈ W ( x ) exp ( - | | d - D k i ( y ) | | 2 2 / T 1 ) · P opt i ( y )
exp ( 1 | | y - x | | 2 2 / T 2 ) .
T 3 T 3 + | I ( x ) - I ( y ) | .
T 4 T 4 + E d ( x , d | I , O )
在等式21中,W(x)可以是中心为x的小窗口。
光流应用器113可基于等式22计算x的双目视差值。
[等式22]
d=arg max P(d|x)
可在图8中描述时间融合的结果。在时间融合之后,一些具有错误深度的区域可被校正。
图4示出根据示例实施例的估计背景视图的方法。
在与背景视图的估计相应的操作220,可使用在操作210被初始化的中间视图的双目视差图或深度图来估计与中间视图相应的背景视图。与预定中间视图相应的背景视图可指示预定中间视图。操作220可包括操作410和420。
背景视图估计器120可推断与中间视图的背景有关的颜色信息和深度信息。推断出的颜色信息和推断出的深度信息可被用于获得高质量的分割结果。这里,分割结果可指示在中间视图中被检测为运动对象的区域。
背景视图估计器120可使用预定帧的邻近帧的中间图像来产生该预定帧的中间视图的背景视图。这里,第k帧的邻近帧可指示第k-k1帧到第k+k2帧。这里,k1和k2可为自然数。也就是说,预定帧的邻近帧可指示包括了该预定帧的时间连续的帧。
可针对每一帧产生背景视图。以下,在一个或多个帧中与由背景视图估计器120产生的背景视图相应的帧可被称为目标帧。
邻近帧弯曲单元121可使用目标帧的邻近帧的中间视图来产生目标帧的背景视图。因此,邻近帧的中间视图的视点可互不相同。因此,邻近帧的中间视图可能需要基于位置、视点等来对目标帧的中间视图进行匹配。
在操作410,邻近帧弯曲单元121可使目标帧的邻近帧的中间视图弯曲以对目标帧的中间视图进行匹配,从而可产生目标帧的中间视图中的像素的候选颜色值和候选双目视差值。邻近帧弯曲单元121可弯曲邻近帧的中间视图以基于邻近帧的中间视图的双目视差图来对目标帧的中间视图进行匹配。
这里,预定像素的侯选颜色值可以是在邻近帧的被弯曲的中间视图中与该预定像素相应的像素的颜色值。侯选双目视差值可以是在邻近帧的被弯曲的中间视图中与预定像素相应的像素的双目视差值。
在操作420,中值滤波单元122可将产生的侯选颜色值的中值推断为侯选颜色值中最可靠的背景值。中值滤波单元122可将中值滤波器施加到目标帧的中间视图中的像素的侯选颜色值和侯选双目视差值,从而产生背景视图的图像和双目视差图。例如,预定像素的侯选颜色值的中值可以是该像素的颜色值。预定像素的侯选双目视差值的中值可以是该像素的双目视差值。
图5示出根据示例实施例的检测运动对象的方法。
在与运动对象的检测相应的操作230,可针对每一像素来对帧的中间视图和背景视图进行比较以从该帧的中间视图检测运动对象。操作230可包括操作510、520和530。
在操作510,差测量单元131可针对每个像素对一个或多个帧中的每一帧的中间视图和背景视图进行比较,从而测量每个像素的色差Dc(x)、深度差Dd(x)和运动深度差Do(x)。
以下将定义将被用于测量像素被包括在动态区域中的概率的三种差。
第一,差测量单元131可基于等式23来计算帧图像I和估计的背景图像IB之间的色差Dc′(x)。
[等式23]
D c ′ ( x ) = | | I ( x ) - I B ( x ) | | 2 | | I ( x ) - I B ( x ) | | 2 + T
在等式23中,T表示用于控制Dc′(x)的形状的参数。差测量单元131可将Dc′(x)归一化至[0,1]。当Dmax表示Dc′(x)的最大值时,差测量单元131可基于等式24来计算I和IB之间的归一化的色差Dc(x)。
[等式24]
Dc(x)=Dc′(x)/Dmax
第二,差测量单元131可基于等式25来计算初始深度和估计的背景深度之间的深度差Dd(x)。
[等式25]
Dd(x)=min{|D(x)-DB(x)|,M}
在等式25中,M表示Dd(x)的最小值的参数。
Dd(x)可被修改为等式26。
[等式26]
Dd(x)=Dd(x)/arg maxdDd(x)
第三,差测量单元131可基于等式27来计算测量初始深度和光流信息的一致性的运动深度差Do′(x)。
[等式27]
D o ′ ( x ) = 1 W Σ i = - n n ( P opt k → k + i ( x ) · | | x o k → k + i - x d k → k + i | | 2 2 )
在等式27中,差测量单元131可基于等式28来计算W。
[等式28]
W = Σ i = - n n P opt k → k + i ( x )
在等式28中,
Figure BDA00002415534600161
表示基于光流信息的在第k+i帧中与x相应的像素。
Figure BDA00002415534600162
表示基于深度信息的在第k+i帧中与x相应的像素。
差测量单元131可基于等式29来计算运动深度差Do(x)。
[等式29]
D o ( x ) = D o ′ ( x ) D o ′ ( x ) + Q
假设与一些区域有关的具有高深度差(Dd)、低色差(Dc)和低运动深度差(Do)的初始双目视差信息存在误差。在该假设下,差测量单元131可通过使能量函数最小化来检测双目视差图中的误差。所述误差的检测可被表示为标签问题。可由等式30来定义标签问题L。
[等式30]
Figure BDA00002415534600164
可由等式31来定义能量函数E(L)。
[等式31]
E(L)=Ed(L)+βEs(L)
可由等式32来定义数据项Ed(L)。
[等式32]
E d ( L ) = Σ x ( T c ( 1 - L ( x ) ) + D d ( x ) max ( D c ( x ) , λ D o ( x ) ) + T d · L ( x ) )
在等式32中,Tc表示门限值,Td表示用于控制L(x)=1的概率的参数。λ表示用于控制运动深度差的权重的参数。
可由等式33定义平滑项。
[等式33]
E s ( L ) = Σ s , t λ ( s , t ) ρ s ( L ( s ) , L ( t ) )
差测量单元131可使用BP算法以使能量函数E(L)最小化。
稍后将参照图9描述三种差异图和深度误差区域。
差测量单元131可在深度误差区域的分割之后计算色差、深度差和运动深度差中的每个的权重。
在操作520,动态区域分割单元132可基于中间视图中的每个像素的色差Dc(x)、深度差Dd(x)和运动深度差Do(x)中的至少一个,通过分割从中间视图获得动态区域。通过分割而获得的动态区域可相应于运动对象。
动态区域分割单元132可基于等式34定义用于动态区域分割的函数C(x)或数据项。
[等式34]
C ( x ) = Σ x D c ( x ) ( 1 - L ( x ) ) + max ( D c ( x ) , D d ( x ) ) 1 + exp ( - ( d ( x ) - M ) / V ) · L ( x )
在等式34中,V表示用于控制C(x)的形状的参数。
动态区域分割单元132可使用基于中间视图中的每个像素的色差Dc(x)、深度差Dd(x)确定值的函数C(x),并可基于从所述函数输出的值来确定中间视图的像素中被包括在动态区域中的像素。当函数C(x)在指定范围内时,动态区域分割单元132可将像素x确定为包括在动态区域中的像素。
在操作530,光流应用器133可通过使用帧中的光流来增强动态区域分割的结果。光流应用器133可通过使用与帧中的光流有关的信息来在分割的动态区域中增强时空一致性。
光流应用器133可将从第k帧到第k+n帧的分割的结果投影到第k帧,从而可产生关于第k帧的2n+1个结果
Figure BDA00002415534600172
可由等式35和等式36来定义关于像素x的分割的结果为l的概率P(l|x)。
[等式35]
P ′ ( d | x ) = Σ i = 1 2 n + 1 Σ y ∈ W ( x ) exp ( - | | d - S k i ( y ) | | 2 2 / T 1 ) · P opt i ( y )
exp ( - | | y - x | | 2 2 / T 2 ) ·
T 3 T 3 + | I ( x ) - I ( y ) |
[等式36]
P(l|x)=P′(l|x)/(P′(l=0|x)+P′(l=1|x))
P(l|x)可用作先验分布,并可被插入到如等式37所示的数据项Ed′(S)。
[等式37]
E d ′ ( S ) = Σ x C ( x ) ( 1 - S ( x ) ) + ( 1 - C ( x ) ) S ( x ) - α log ( P ( S ( x ) | x ) )
在等式37中,α表示用于控制Ed′(S)的形状的参数。
可由如上所述的等式6来定义平滑项Es(L)。
光流应用器133可通过使用BP算法使能量函数最小化。在能量函数被最小化之后,光流应用器133可基于空间信息和时间信息来使用局部分割方法,从而增强动态区域分割的结果。
光流应用器133可在边界上产生重叠的窗口。这里,所述边界是指被分割的动态区域的边界。光流应用器133可在每个窗口内计算前景高斯混合模型(GMM)和背景GMM。光流应用器133可基于前景GMM和背景GMM来增强动态区域分割的结果。
当像素x和边界之间的距离大于门限Q时,像素x包括在前景中的概率
Figure BDA00002415534600185
可以是S(x)。此外,像素x包括在背景中的概率
Figure BDA00002415534600186
可以是
Figure BDA00002415534600187
当像素x与边界之间的距离小于或等于门限Q时,像素x可被视为包括在k个局部窗口内。因此,光流应用器133可基于等式38计算像素x包括在前景中的概率
Figure BDA00002415534600188
[等式38]
P F c ′ ( x ) = 1 k Σ i = 1 k P i ( I ( x ) | s = 1 ) P i ( I ( x ) | s = 0 ) + P i ( I ( x ) | s = 1 )
在等式38中,Pi(I(x)|s=1)可以是像素x包括在前景中的概率,通过使用第i局部窗口的前景GMM来计算所述概率。类似地,光流应用器133可基于等式39来计算像素x包括在背景中的概率
Figure BDA00002415534600192
[等式39]
P B c ′ ( x ) = 1 k Σ i = 1 k P i ( I ( x ) | s = 1 ) P i ( I ( x ) | s = 0 ) + P i ( I ( x ) | s = 1 )
当分别由等式40和等式41定义
Figure BDA00002415534600194
Figure BDA00002415534600195
时,
Figure BDA00002415534600196
Figure BDA00002415534600197
可被视为先验分布。
[等式40]
P F c ( x ) = P F c ′ ( x ) / ( P F c ′ ( x ) + P B c ′ ( x ) )
[等式41]
P B c ( x ) = 1 - P F c ( x )
可从动态区域分割的初始结果获得的形状信息对于减小分割误差会是有用的。中心为像素x的小窗口中的前景像素的数量由N表示,与包括在前景中的像素x相关的形状信息
Figure BDA000024155346001910
可由等式42定义。
[等式42]
P F s ( x ) = N w 2
在等式42,w表示窗口的宽度。此外,与包括在背景中的像素x相关的形状信息
Figure BDA000024155346001912
可由等式43定义。
[等式43]
P B s ( x ) = 1 - P F s ( x )
此外,光流应用器133可使用空间信息和时间信息,从而增强动态区域分割的结果。光流应用器133可基于等式36来增强动态区域分割的结果的时空一致性。
与包括在前景中的像素x有关的时空一致性信息
Figure BDA00002415534600201
可由等式44定义。
[等式44]
P F t ( x ) = P ( l = 1 | x )
与包括在背景中的像素x有关的时空一致性信息
Figure BDA00002415534600203
可由等式45定义。
[等式45]
P B t ( x ) = P ( l = 0 | x )
通过组合三种概率获得的PF(x)和P′B(x)可由等式46和等式47定义。
[等式46]
P F ′ ( x ) = P F c ( x ) · P F s ( x ) · P F t ( x )
[等式47]
P B ′ ( x ) = P B c ( x ) · P B s ( x ) · P B t ( x )
因此,光流应用器133可基于等式48来计算像素x包括在前景中的概率PF(x)。
[等式48]
PF(x)=P′F(x)/(P′F(x)+P′B(x))
此外,光流应用器133可基于等式49来计算像素x包括在背景中的概率PB(x)。
[等式49]
PB(x)=1-PF(x)
如等式50所示,PF(x)和PB(x)可被加到数据项。
[等式50]
E d ( S ) = Σ x ( - log ( P F ( x ) ) ( 1 - S ( x ) ) - log ( P B ( x ) ) S ( x ) )
平滑项Es(D)可由上述的等式6定义。光流应用器133可基于BP算法来获得动态区域分割的结果。
将参照图10描述被分割的运动对象。
图6示出根据示例实施例的时空立体匹配的方法。
在与时空立体匹配相应的操作240中,束优化和时间优化可被应用到中间视图中的静态像素和动态像素。操作240可包括操作610、620和630。
在操作610,静态像素优化器141可将束优化模型应用到中间视图中的静态像素,从而增强静态像素的双目视差值。
假设在短时间内的动态像素的轨迹可被近似为3D线。
在操作620,在增强动态像素的双目视差值之前,轨迹追踪器142可在这样的假设下计算动态像素的3D轨迹。
可由等式51来表示在3D空间中的线(以下被称为线l)。
[等式51]
l:x=x0+tu
在等式51,x0表示线l上的点,u表示方向矢量。当x被假设是第k帧中的二维(2D)点时,与邻近帧中的x相应的位置可被视为x的轨迹。{x1′,x2′,...,和xn′}表示不同帧中的x的位置。轨迹追踪器142可基于相机参数和估计的深度图来将位置映射到3D点{x1,x2,...,和xn}。随后,轨迹追踪器142可基于主成分分析(PCA)来计算动态像素的线性轨迹。
在操作630,动态像素优化器143可通过将使用动态像素的轨迹的时间优化和束优化应用到动态像素,来增强动态像素的双目视差值。
动态像素优化器143可基于等式52来计算通过映射产生的n个点{x1,x2,...,和xn}的均值u。
[等式52]
u = Σ k = 1 n x k / n
动态像素优化器143可基于等式53来计算协方差矩阵S。
[等式53]
S = Σ k = 1 n ( x k - n ) ( x k - u ) T / n
S表示对称矩阵。因此,可如等式54来分解S。
[等式54]
S=Q∑QT
在等式54,Q表示正交矩阵,QT表示对角矩阵。
当λ1表示S的第一最大特征值,λ2表示S的第二最大特征值时,由等式55定义的参数α可指示是否由直线来表示点的轨迹。
[等式55]
α=λ12
当α大于与特征值有关的特征矢量C时,由直线来表示点的轨迹的假设可以是合理的。当α小于特征矢量C时,该假设可以是无效的。
在第一个例子中,在x的双目视差值为d的假设下,当轨迹线与x在d下映射到的3D点之间的距离长时,x的双目视差为d的概率可以为低。当轨迹线与x在d下映射到的3D点之间的距离短时,x的双目视差为d的概率可以为高。因此,为了执行对x的双目视差为d的概率P(d|x)建模,可使用3D空间中的点和线之间的距离。
动态像素优化器143可基于等式56来计算P(d|x)。
[等式56]
P ( d | x ) = exp ( - Q ( d | x ) 2 v ( x ) 2 · σ p )
在等式56,Q(d|x)表示点和线之间的距离。可由等式57来定义v(x)。σp可以是用于控制P(d|x)的形状的参数。
[等式57]
v(x)=mindQ(d|x)
在第二个例子中,动态像素优化器143可使用多个帧中的光流来将第k-l帧和第k+1帧的双目视差图映射到第k帧。通过映射,第k帧可具有三个双目视差图
Figure BDA00002415534600231
针对第k帧中的像素x,动态像素优化器143可基于等式58来计算像素x的双目视差值为d的概率P(d|x)。
[等式58]
P ( d | x ) = Σ i = 1 3 Σ y ∈ W ( x ) exp ( - | | d - D k i ( y ) | | 2 2 / T 1 ) · P opt i ( y )
exp ( - | | y - x | | 2 2 / T 2 ) ·
T 3 T 3 + | I ( x ) - I ( y ) |
在等式58中,W(x)表示中心为x的小窗口。在等式58中,第一项指示d与初始双目视差值Dk(y)之间的相似度。第二项可指示第Ki帧的光流的可靠性。第三项可测量与像素y有关的信息的权重。第四项可类似于等式6的平滑项ES(D)。
动态像素优化器143可基于等式59来计算双目视差似然性P(d|C),从而组合束优化模型和局部3D线模型。
[等式59]
P ( d | C )
= Σ z = 0,1 P ( d , z | C )
= Σ z = 0,1 P ( d | z , C ) P ( z | C )
= Σ Z = 0,1 P ( d | z , C ) P ( C | z ) P ( z ) Σ z = 0,1 P ( C | z ) P ( z )
在等式59中,C表示像素x的颜色值。为了清楚,在等式59中可忽略像素x的位置。z的值为0可指示像素x包括在静态背景中。z的值为1可指示像素x包括在动态背景中。
当像素x与分割边界之间的距离大于门限Q时,可由等式60定义像素x包括在前景中的概率P(C|z=1)。
[等式60]
P(C|z=1)=S(x)
当像素x与分割边界之间的距离大于门限Q时,可由等式61来表示像素x包括在背景中的概率P(C|z=0)。
[等式61]
P(C|z=0)=1-P(C|z=1)
当像素x与分割边界之间的距离小于或等于Q时,动态像素优化器143可基于在局部窗口中计算的前景GMM和背景GMM来获得P(C|z)。
动态像素优化器143可通过以下三个步骤来计算P(z)。
在第一步骤中,动态像素优化器143可通过使用像素x与分割边界之间的距离d,基于等式62和等式63来计算Pd(z=1)和Pd(z=0)。
[等式62]
Pd(z=1)=exp(-d2d)
[等式63]
Pd(z=0)=1-Pd(z=1)
在第二步骤中,形状信息PS(z)可被定义为与等式42的
Figure BDA00002415534600251
相同。
在第三步骤中,时空一致性Pt(z)可被定义为与等式36的P(l|x)相同。
通过以下描述的步骤,动态像素优化器143可基于等式64来计算P(z)。
[等式64]
P(z)=Pd(z)·Ps(z)·Pt(z)
P(d|z=1,C)可被定义为与等式56或等式58的P(d|x)相同。通过使用束优化算法中的几何一致性项Pv,双目视差空间中的几何一致性P(d|z=0,C)可被定义为与P(d|z=1,C)类似。不同帧中的相应像素具有相同的3D位置的现象可以是直观的。微小的差异在于原始束优化模型测量图像空间中的几何一致性。
动态像素优化器143可基于等式65来测量双目视差空间中的几何一致性P(d|z=0,C)。
[等式65]
P ( d | z = 0 , C ) = exp ( - | | h ( d ) - d x ′ | | 2 2 σ v 2 )
在等式65中,x′表示匹配帧i′中的投影的像素。dx表示估计的x′的双目视差值。h(d)表示基于相机参数转换的双目视差值。[R|t]表示从帧i到帧i′的投影矩阵。Ki表示帧i的内在矩阵。x的2D位置可以是(u,v)。
动态像素优化器143可基于等式66计算h(d)。
[等式66]
h ( d ) = 1 ( RK i - 1 ( u , v , 1 / d ) T + t ) [ 3 ]
在等式63中,[3]表示矢量的三个元素。
动态像素优化器143可通过将P(d|C)与由等式3定义的数据项Ed相加,基于等式67来计算动态区域的数据项
[等式67]
E d 1 ( D , O ; I ) = E d ( D , O ; I ) + λ Σ x ( - log P ( D ( x ) | I ( x ) ) )
等式67可应用到动态区域。静态像素优化器141可使用关于静态背景的束优化数据项
Figure BDA00002415534600263
可如等式68所示来组合
[等式68]
E d c ( D , O ; I ) = Σ x ( E d 1 ( D ( x ) , O ; I ) · P ( z = 1 | x , I ( x ) ) +
E d 2 ( D ( x ) , O ; I ) · P ( z = 0 | x , I ( x ) ) )
平滑项
Figure BDA00002415534600268
可被定义为与等式6的ES(D)相同。
为了获得最终结果,时空立体匹配器140可基于BP算法来解能量函数。
将参照图7来描述在时空优化之后增强的双目视差图。
图7示出根据示例实施例的产生左右双目视差图的方法。
除了恢复中间视图的双目视差图之外,左视图和右视图的双目视差图可能需要被估计。
操作250可包括操作710、720和730。
在操作710,左右双目视差图初始化单元151可对一个或多个帧的中间视图进行弯曲,从而产生所述帧的左视图的初始双目视差图和右视图的初始双目视差图。
左右双目视差图初始单元151可将中间视图中的每个帧的估计的运动对象的遮蔽和双目视差图弯曲到左视图和右视图。经过弯曲的双目视差图可具有由于中间视图中的遮蔽导致丢失像素而产生的孔。
在操作720,基于分割的孔填充单元152可按照基于分割的方法来推断丢失像素的深度,从而填充左视图的双目视差图和右视图的双目视差图中的孔。
基于分割的孔填充单元152可组合通过分割运动对象而产生的分割信息,从而推断丢失像素的双目视差值。通常,大多数丢失像素存在于静态区域中。
基于分割的孔填充单元152可使用均值偏移算法来分割左图像和右图像。当分段包括静态像素和动态像素时,基于分割的孔填充单元152可执行进一步的分割。
在每个分段的3D表面可被近似为3D平面的假设下,基于分割的孔填充单元152可针对具有丢失像素的每个分段执行平面拟合。基于分割的孔填充单元152可使用针对每个分段的平面拟合,从而推断关于丢失像素的双目视差值。
由于具有精确的深度的像素可能不足,因此遮蔽区域中的一些分段可能被错误地拟合。如上所述的时空优化方法可被引入,从而增强左双目视差图和右双目视差图。
在操作730,左右双目视差图时空优化单元153可将时空优化应用到左视图的双目视差图和右视图的双目视差图。
图8示出根据示例实施例的时间融合的结果。
第一图像810可相应于在时间融合之前估计的双目视差图。第二图像820可相应于在时间融合之后估计的双目视差图。
图9示出根据示例实施例的三种差异图和深度误差区域。
第一图像910可相应于色差图。第二图像920可相应于深度差图。第三图像930可相应于运动深度差图。第四图像940可相应于基于色差图、深度差图和运动深度差图产生的指示深度误差区域的图像。
图10示出根据示例实施例的从视图分割的运动对象。
第一图像1010可相应于视图的图像。第二图像1020可包括背景和通过分割从视图检测的运动对象。在第二图像1020中用黑色表示的部分可相应于背景。在第二图像1020中用白色表示的部分可相应于通过分割检测的动态对象。
图11示出根据示例实施例的在执行时空优化之后的增强的双目视差图。
第一图像1110可相应于视图的图像。第二图像1120可以是通过时空优化而增强的双目视差图。
根据上述示例实施例的用于恢复图像的深度信息的设备和方法可使用一个或多个处理器,所述一个或多个处理器可包括微处理器、中央处理单元(CPU)、数字信号处理器(DSP)或专用集成电路(ASIC),以及这些处理装置和其它处理装置的部分或组合。
这里使用的术语“模块”和“单元”可以表示但不限于执行特定任务的软件或硬件组件或者装置,诸如现场可编程门阵列(FPGA)或专用集成电路(ASIC)。模块或单元可被配置为位于可寻址的存储介质上并可被配置为在一个或多个处理器上执行。因此,模块或单元可包括(以示例的方式)组件(诸如软件组件、面向对象的软件组件、类组件和任务组件)、处理、功能、属性、程序、子程序、程序代码的片段、驱动器、固件、微码、电路、数据、数据库、数据结构、表、阵列和变量。在组件和模块/单元中提供的功能可被组合为更少的组件和模块/单元,或者可被进一步分离为另外的组件和模块。
根据上述实施例的方法可被记录在非暂时性计算机可读介质中,所述非暂时性计算机可读介质包括用于执行由计算机实施的各种操作的程序指令。所述介质还可包括单独的程序指令、数据文件、数据结构等或者可包括程序指令、数据文件、数据结构等的组合。非暂时性计算机可读介质的示例包括磁介质(诸如硬盘、软盘和磁带);光学介质(诸如CD ROM盘和DVD);磁光介质(诸如光盘)和专门配置以存储和执行程序指令的硬件装置(诸如只读存储器(ROM)、随机存取存储器(RAM)、闪存等)。程序指令的示例包括机器代码(诸如由编译器所产生的)和包含可由计算机使用翻译器执行的高级代码的文件。所述硬件装置可被配置以用作一个或多个软件模块以执行上述实施例的操作,反之亦然。
流程图的每个块可代表包括用于执行指定逻辑功能的一个或多个可执行指令的代码的一个单元、模块、片段或部分。应注意,在一个可选的实施方式中,块中记录的功能可不按顺序发生。例如,根据有关的功能,连续示出的两个块实际上可完全同时被执行,或者所述块有时可以以相反的顺序被执行。
虽然已示出和描述了示例实施例,但是本领域技术人员将理解,在不脱离本公开的原理和精神的情况下,可对这些实施例进行改变,本公开的范围由权利要求及其等同物限定。

Claims (22)

1.一种处理图像的方法,所述方法包括:
使用包括在一个或多个帧中的每一帧中的左图像、中间图像和右图像来估计所述一个或多个帧中的每一帧的多个视图中的中间视图的双目视差图,中间视图包括中间图像和估计的双目视差图;
使用所述一个或多个帧的中间视图的深度图来估计指示中间视图的背景的背景视图;
基于中间视图和背景视图之差,检测中间视图中的运动对象。
2.如权利要求1所述的方法,其中,左图像和右图像中的至少一个相应于多个图像。
3.如权利要求1所述的方法,其中,估计双目视差图的步骤包括:
当左图像、中间图像和右图像被提供时,通过使用反向传播(BP)算法重复使能量函数的值最小化来计算中间视图的双目视差图,基于与中间视图有关的遮蔽信息和中间视图的双目视差图来确定能量函数的值;
通过融合将不同的平面拟合应用到中间视图的双目视差图而获得的结果来增强中间视图的双目视差图。
4.如权利要求3所述的方法,其中,计算的步骤包括:
基于中间视图的双目视差图来估计遮蔽信息;
基于遮蔽信息来估计中间视图的双目视差图。
5.如权利要求3所述的方法,还包括:
使用与所述一个或多个帧中的光流有关的信息来增强中间视图的双目视差图。
6.如权利要求1所述的方法,其中,估计背景视图的步骤包括:
对所述一个或多个帧中的目标帧的邻近帧的中间视图进行弯曲以对目标帧的中间视图进行匹配,从而产生包括在目标帧的中间视图中的像素的候选双目视差值和候选颜色值;
将中值滤波器应用到像素的候选双目视差值和候选颜色值,以产生目标帧的背景视图的图像和双目视差图。
7.如权利要求1所述的方法,其中,检测运动对象的步骤包括:
针对每个像素比较所述一个或多个帧中的每一帧的中间视图和背景视图,以测量包括在中间视图和背景视图中的每个像素之间的颜色、深度和运动深度的差;
基于包括在中间视图中的每个像素的颜色、深度和运动深度的差中的至少一个,从中间视图分割出动态区域,
其中,分割出的动态区域相应于运动对象。
8.如权利要求7所述的方法,其中,检测的步骤还包括:
使用与所述一个或多个帧中的光流有关的信息来增强分割出的动态区域中的时空一致性。
9.如权利要求1所述的方法,还包括:
将束优化和基于局部三维(3D)线模型的时间优化应用到所述一个或多个帧的中间视图,以增强包括在中间视图中的静态像素和动态像素的双目视差值,
其中,动态像素相应于与运动对象对应的像素,静态像素相应于与运动对象不对应的像素。
10.如权利要求1所述的方法,还包括:
通过对所述一个或多个帧的中间视图进行弯曲来产生所述一个或多个帧的左视图的双目视差图和所述一个或多个帧的右视图的双目视差图;
基于分割,填充包括在左视图的双目视差图和右视图的双目视差图中的孔;
将时空优化应用到左视图的双目视差图和右视图的双目视差图,
其中,所述孔相应于由中间视图中的遮蔽而产生的丢失像素。
11.一种用于处理图像的设备,所述设备包括:
双目视差图估计器,使用包括在一个或多个帧中的每一帧中的左图像、中间图像和右图像来估计所述一个或多个帧中的每一帧的多个视图中的中间视图的双目视差图,所述中间视图包括中间图像和估计的双目视差图;
背景视图估计器,使用所述一个或多个帧的中间视图的深度图来估计指示中间视图的背景的背景视图;
运动对象检测器,基于中间视图和背景视图之差,检测中间视图中的运动对象。
12.如权利要求11的设备,其中,双目视差图估计器包括:
初始双目视差图计算器,当左图像、中间图像和右图像被提供时,通过使用反向传播(BP)算法重复使能量函数的值最小化来计算中间视图的双目视差图,基于中间视图的双目视差图和与中间视图有关的遮蔽信息来确定能量函数的值;
平面拟合应用器,通过融合将不同的平面拟合应用到中间视图的双目视差图而获得的结果来增强中间视图的双目视差图。
13.如权利要求12的设备,其中,初始双目视差图计算器基于中间视图的双目视差图来估计遮蔽信息,并基于遮蔽信息来估计中间视图的双目视差图。
14.如权利要求12的设备,还包括:
光流应用器,使用与所述一个或多个帧中的光流有关的信息来增强中间视图的双目视差图。
15.如权利要求11的设备,其中,背景视图估计器包括:
邻近帧弯曲单元,对所述一个或多个帧中的目标帧的邻近帧的中间视图进行弯曲以对目标帧的中间视图进行匹配,从而产生包括在目标帧的中间视图中的像素的候选双目视差值和候选颜色值;
中值滤波单元,将中值滤波器应用到像素的候选双目视差值和候选颜色值以产生目标帧的背景视图的图像和双目视差图。
16.如权利要求11的设备,其中,运动对象检测器包括:
差测量单元,针对每个像素比较所述一个或多个帧中的每一帧的中间视图和背景视图,以测量包括在中间视图和背景视图中的每个像素之间的颜色、深度和运动深度的差;
动态区域分割单元,基于包括在中间视图中的每个像素的颜色、深度和运动深度的差中的至少一个,从中间视图分割出动态区域,
其中,分割出的动态区域相应于运动对象。
17.如权利要求16的设备,其中,运动对象检测器还包括:
光流应用器,使用与所述一个或多个帧中的光流有关的信息来增强分割出的动态区域中的时空一致性。
18.如权利要求11的设备,还包括:
时空立体匹配器,将束优化和基于局部三维(3D)线模型的时间优化应用到所述一个或多个帧的中间视图,以增强包括在中间视图中的静态像素和动态像素的双目视差值,
其中,动态像素相应于与运动对象对应的像素,静态像素相应于与运动对象不对应的像素。
19.如权利要求11的设备,还包括:
左右双目视差图初始化单元,通过对所述一个或多个帧的中间视图进行弯曲来产生所述一个或多个帧的左视图的双目视差图和所述一个或多个帧的右视图的双目视差图;
基于分割的孔填充单元,基于分割填充包括在左视图的双目视差图和右视图的双目视差图中的孔;
左右双目视差图时空优化器,将时空优化应用到左视图的双目视差图和右视图的双目视差图,
其中,所述孔相应于由中间视图中的遮蔽而产生的丢失像素。
20.一种处理图像的方法,所述方法包括:
使用包括在一个或多个帧中的每一帧中的左图像、中间图像和右图像来估计所述一个或多个帧中的每一帧的中间视图的双目视差图;
使用所述一个或多个帧的中间视图的深度图来估计中间视图的背景视图;
基于中间视图与背景视图之差来检测中间视图中的运动对象;
使用重复优化处理来计算中间视图的双目视差图,其中,所述重复优化处理使用反向传播算法来使能量函数的值最小化,
其中,基于与中间视图有关的遮蔽信息确定能量函数的值。
21.如权利要求20的方法,其中,检测的步骤还包括:
针对每个像素比较所述一个或多个帧中的每一帧的中间视图和背景视图,测量包括在中间视图和背景视图中的每个像素之间的颜色、深度和运动深度的差;
基于测量的颜色、深度和运动深度的差中的至少一个,对中间视图进行分割以获得一个或多个分割出的动态区域。
22.如权利要求21所述的方法,其中,检测的步骤还包括:
针对分割出的动态区域的一个或多个边界,计算前景高斯混合模型和背景高斯混合模型。
CN201210464952.XA 2011-11-30 2012-11-16 用于恢复图像的深度信息的方法和设备 Active CN103136726B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
KR1020110127153A KR101918057B1 (ko) 2011-11-30 2011-11-30 이미지의 깊이 정보를 복원하는 방법 및 장치
KR10-2011-0127153 2011-11-30

Publications (2)

Publication Number Publication Date
CN103136726A true CN103136726A (zh) 2013-06-05
CN103136726B CN103136726B (zh) 2017-03-01

Family

ID=48466897

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210464952.XA Active CN103136726B (zh) 2011-11-30 2012-11-16 用于恢复图像的深度信息的方法和设备

Country Status (3)

Country Link
US (1) US8855408B2 (zh)
KR (1) KR101918057B1 (zh)
CN (1) CN103136726B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107209565A (zh) * 2015-01-20 2017-09-26 微软技术许可有限责任公司 固定大小的增强现实对象
CN107610070A (zh) * 2017-09-29 2018-01-19 深圳市佳创视讯技术股份有限公司 基于三摄像采集的自由立体匹配方法
CN110798676A (zh) * 2019-11-29 2020-02-14 苏州新光维医疗科技有限公司 一种利用内镜镜头动态图像形成3d视觉的方法及装置
CN113808185A (zh) * 2021-11-19 2021-12-17 北京的卢深视科技有限公司 图像深度恢复方法、电子设备及存储介质

Families Citing this family (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9098908B2 (en) 2011-10-21 2015-08-04 Microsoft Technology Licensing, Llc Generating a depth map
US9165401B1 (en) * 2011-10-24 2015-10-20 Disney Enterprises, Inc. Multi-perspective stereoscopy from light fields
US9113043B1 (en) * 2011-10-24 2015-08-18 Disney Enterprises, Inc. Multi-perspective stereoscopy from light fields
US9659372B2 (en) * 2012-05-17 2017-05-23 The Regents Of The University Of California Video disparity estimate space-time refinement method and codec
US9235879B2 (en) * 2012-06-29 2016-01-12 Hong Kong Applied Science And Technology Research Institute Co., Ltd. Apparatus, system, and method for temporal domain hole filling based on background modeling for view synthesis
US20140241612A1 (en) * 2013-02-23 2014-08-28 Microsoft Corporation Real time stereo matching
CN104427324A (zh) * 2013-09-02 2015-03-18 联咏科技股份有限公司 视差计算方法及其立体匹配装置
CN104735360B (zh) * 2013-12-18 2017-12-22 华为技术有限公司 光场图像处理方法和装置
CN104616286B (zh) * 2014-12-17 2017-10-31 浙江大学 快速的半自动多视图深度修复方法
WO2017014691A1 (en) 2015-07-17 2017-01-26 Heptagon Micro Optics Pte. Ltd. Generating a distance map based on captured images of a scene
US10699476B2 (en) 2015-08-06 2020-06-30 Ams Sensors Singapore Pte. Ltd. Generating a merged, fused three-dimensional point cloud based on captured images of a scene
WO2017030507A1 (en) * 2015-08-19 2017-02-23 Heptagon Micro Optics Pte. Ltd. Generating a disparity map having reduced over-smoothing
CN105741269A (zh) * 2016-01-25 2016-07-06 中国科学院深圳先进技术研究院 视频切割的方法及装置
CN106123865B (zh) * 2016-06-22 2019-01-29 电子科技大学 面向虚拟图像的机器人导航方法
CN106056622B (zh) * 2016-08-17 2018-11-06 大连理工大学 一种基于Kinect相机的多视点深度视频复原方法
CN106530409B (zh) * 2016-11-03 2019-08-27 浙江大学 立体匹配中局部区域一致性对应方法
CN108961383B (zh) * 2017-05-19 2021-12-14 杭州海康威视数字技术股份有限公司 三维重建方法及装置
CN107220994A (zh) * 2017-06-01 2017-09-29 成都通甲优博科技有限责任公司 一种立体匹配的方法及系统
CN110390645B (zh) 2018-04-23 2024-04-16 康耐视公司 用于立体瞬时图像序列的改善3d数据重构的系统和方法
WO2020006519A1 (en) 2018-06-29 2020-01-02 Eloupes, Inc. Synthesizing an image from a virtual perspective using pixels from a physical imager array
US10803606B2 (en) * 2018-07-19 2020-10-13 National Taiwan University Temporally consistent belief propagation system and method
CN109949397A (zh) * 2019-03-29 2019-06-28 哈尔滨理工大学 一种结合激光点和均值漂移的深度图重建方法
CN110176032B (zh) * 2019-04-28 2021-02-26 暗物智能科技(广州)有限公司 一种三维重建方法及装置
CN110470216B (zh) * 2019-07-10 2022-01-28 湖南交工智能技术有限公司 一种三镜头高精度视觉测量方法及装置
CN115104137A (zh) 2020-02-15 2022-09-23 利蒂夫株式会社 提供基于体育视频的平台服务的服务器的操作方法
CN111369452A (zh) * 2020-02-26 2020-07-03 青海民族大学 大区域图像局部破损点优化提取方法
US11741675B2 (en) * 2020-03-10 2023-08-29 Niantic, Inc. Determining traversable space from single images
US20230245322A1 (en) * 2022-01-28 2023-08-03 Samsung Electronics Co., Ltd. Reconstructing A Three-Dimensional Scene
TWI787141B (zh) * 2022-06-21 2022-12-11 鴻海精密工業股份有限公司 深度估計模型訓練方法、深度估計方法及電子設備

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070031037A1 (en) * 2005-08-02 2007-02-08 Microsoft Corporation Stereo image segmentation
US20090141938A1 (en) * 2007-11-08 2009-06-04 Elctronics And Telecommunications Research Institute Robot vision system and detection method
CN102117486A (zh) * 2009-12-30 2011-07-06 三星电子株式会社 生成三维图像数据的方法和装置
CN101496413B (zh) * 2006-08-01 2011-09-07 高通股份有限公司 用单视场低功率移动装置实时捕获及产生立体图像及视频

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100603603B1 (ko) 2004-12-07 2006-07-24 한국전자통신연구원 변위 후보 및 이중 경로 동적 프로그래밍을 이용한 스테레오 변위 결정 장치 및 그 방법
EP1806697B1 (en) 2006-01-10 2016-08-10 Microsoft Technology Licensing, LLC Segmenting image elements
KR20100008649A (ko) 2008-07-16 2010-01-26 광주과학기술원 참조영상을 이용한 깊이영상 생성방법 및 그 장치, 생성된깊이영상을 부호화/복호화하는 방법 및 이를 위한인코더/디코더, 그리고 상기 방법에 따라 생성되는 영상을기록하는 기록매체
KR101099648B1 (ko) 2009-06-22 2011-12-29 한국과학기술원 입체 영상에서의 폐색을 처리하는 입체 정합 시스템 및 방법, 폐색 처리 알고리즘을 기록한 컴퓨터에 의해 판독가능한 매체
KR101192313B1 (ko) * 2010-04-14 2012-10-17 전자부품연구원 깊이 영상의 시간적 상관도 향상 방법

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070031037A1 (en) * 2005-08-02 2007-02-08 Microsoft Corporation Stereo image segmentation
CN101496413B (zh) * 2006-08-01 2011-09-07 高通股份有限公司 用单视场低功率移动装置实时捕获及产生立体图像及视频
US20090141938A1 (en) * 2007-11-08 2009-06-04 Elctronics And Telecommunications Research Institute Robot vision system and detection method
CN102117486A (zh) * 2009-12-30 2011-07-06 三星电子株式会社 生成三维图像数据的方法和装置

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107209565A (zh) * 2015-01-20 2017-09-26 微软技术许可有限责任公司 固定大小的增强现实对象
CN107209565B (zh) * 2015-01-20 2020-05-05 微软技术许可有限责任公司 用于显示固定大小的增强现实对象的方法和系统
CN107610070A (zh) * 2017-09-29 2018-01-19 深圳市佳创视讯技术股份有限公司 基于三摄像采集的自由立体匹配方法
CN107610070B (zh) * 2017-09-29 2020-12-01 深圳市佳创视讯技术股份有限公司 基于三摄像采集的自由立体匹配方法
CN110798676A (zh) * 2019-11-29 2020-02-14 苏州新光维医疗科技有限公司 一种利用内镜镜头动态图像形成3d视觉的方法及装置
CN113808185A (zh) * 2021-11-19 2021-12-17 北京的卢深视科技有限公司 图像深度恢复方法、电子设备及存储介质
CN113808185B (zh) * 2021-11-19 2022-03-25 北京的卢深视科技有限公司 图像深度恢复方法、电子设备及存储介质

Also Published As

Publication number Publication date
CN103136726B (zh) 2017-03-01
KR101918057B1 (ko) 2019-01-30
US20130136299A1 (en) 2013-05-30
US8855408B2 (en) 2014-10-07
KR20130060865A (ko) 2013-06-10

Similar Documents

Publication Publication Date Title
CN103136726A (zh) 用于恢复图像的深度信息的方法和设备
US11100401B2 (en) Predicting depth from image data using a statistical model
Behl et al. Bounding boxes, segmentations and object coordinates: How important is recognition for 3d scene flow estimation in autonomous driving scenarios?
US8249333B2 (en) Segmenting image elements
Wedel et al. Efficient dense scene flow from sparse or dense stereo data
CN102750711B (zh) 一种基于图像分割和运动估计的双目视频深度图求取方法
Zitnick et al. Stereo for image-based rendering using image over-segmentation
US7755619B2 (en) Automatic 3D face-modeling from video
Yang et al. Near Real-time Stereo for Weakly-Textured Scenes.
Sun et al. Local layering for joint motion estimation and occlusion detection
Yang Fast stereo matching using reliability-based dynamic programming and consistency constraints
US8126273B2 (en) Method for reconstructing three-dimensional images from two-dimensional image data
Lee et al. Silhouette segmentation in multiple views
Lim et al. Generalized background subtraction using superpixels with label integrated motion estimation
US20050286758A1 (en) Color segmentation-based stereo 3D reconstruction system and process employing overlapping images of a scene captured from viewpoints forming either a line or a grid
Huq et al. Occlusion filling in stereo: Theory and experiments
Papadakis et al. Multi-label depth estimation for graph cuts stereo problems
Zhang et al. Recovering consistent video depth maps via bundle optimization
US20220270354A1 (en) Monocular image-based model training method and apparatus, and data processing device
Zhang et al. As-rigid-as-possible stereo under second order smoothness priors
Bulatov et al. Multi-view dense matching supported by triangular meshes
Zhang et al. Simultaneous multi-body stereo and segmentation
Yang et al. Consistent depth maps recovery from a trinocular video sequence
Unger et al. Probabilistic disparity fusion for real-time motion-stereo
Birchfield et al. Correspondence as energy-based segmentation

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant