用于利用使用幅度调制的连续光的飞行时间相机进行深度测量的方法

文档序号:74678 发布日期:2021-10-01 浏览:47次 >En<

阅读说明:本技术 用于利用使用幅度调制的连续光的飞行时间相机进行深度测量的方法 (Method for depth measurement with a time-of-flight camera using amplitude modulated continuous light ) 是由 B·米尔巴赫 M·博古斯劳斯基 T·索利尼亚克 于 2020-02-14 设计创作,主要内容包括:本发明涉及一种利用使用幅度调制的连续光的飞行时间相机(1)进行深度测量的方法,所述方法包括以下步骤:针对相机(1)的传感器阵列(2)的多个像素(3)中的每个像素,以高于所述幅度调制的连续光的调制频率的采样频率获取(120)至少一个样本序列,所述样本序列包括至少四个幅度样本(A-0,A-1,A-2,A-3)。为了能够使用飞行时间相机进行准确并且有效的深度测量,本发明规定该方法还包括:—针对每个像素(3)的每个样本序列确定(130)置信度值(C),所述置信度值指示幅度样本(A-0,A-1,A-2,A-3)与幅度的正弦时间演变的对应程度;以及—针对多个分箱区域(4)中的每个分箱区域基于来自所述分箱区域(4)的像素(3)的样本序列的幅度样本(A-0,A-1,A-2,A-3)来确定(260)分箱的深度值(D-b),所述每个分箱区域包括多个像素(3),其中,所述样本序列对所述分箱的深度值(D-b)的贡献取决于其置信度值(C)。(The invention relates to a method for depth measurement with a time-of-flight camera (1) using amplitude modulated continuous light, comprising the following steps: sensor array for a camera (1)(2) At least one sample sequence comprising at least four amplitude samples (a) at a sampling frequency higher than the modulation frequency of the amplitude modulated continuous light is acquired (120) 0 ,A 1 ,A 2 ,A 3 ). In order to enable an accurate and efficient depth measurement using a time-of-flight camera, the invention provides that the method further comprises: -determining (130), for each sequence of samples of each pixel (3), a confidence value (C) indicative of the amplitude samples (a) 0 ,A 1 ,A 2 ,A 3 ) Degree of correspondence with sinusoidal temporal evolution of amplitude; and-for each of a plurality of binning regions (4), based on amplitude samples (a) of a sequence of samples from pixels (3) of the binning region (4) 0 ,A 1 ,A 2 ,A 3 ) To determine (260) a binned depth value (D) b ) Each binning region comprising a plurality of pixels (3), wherein the depth values (D) of the sequence of samples for the binning b ) The contribution of (C) depends on its confidence value (C).)

用于利用使用幅度调制的连续光的飞行时间相机进行深度测 量的方法

技术领域

概括地说,本公开内容涉及一种利用使用幅度调制的连续光的飞行时间相机进行深度测量的方法。

背景技术

飞行时间相机用于提供在三维对象或风景的图像中的像素深度信息。相机包括具有多个像素的(通常为二维的)传感器阵列。每个像素都提供信息,从中可以得出空间中记录点的深度(即距相机的距离)。除了使用光脉冲的TOF相机外,另一种类型的TOF相机使用幅度调制的连续光。换句话说,相机发出连续场的幅度调制的光,该光被从相机视场中的对象反射。反射光被各个像素接收。由于幅度调制,可以从幅度推导出接收的光的相位,并可以通过相对相位差导出飞行时间,并且从而确定到反射对象的距离。根据公知的方法,采用锁定像素,其中每个像素的读出与光的调制频率同步。具体而言,每个像素的读出频率可以是调制频率的4倍。这也被称为4抽头(4-tap)方法。其基于在四个连续时间点接收和评估四个连续幅度样本,每个时间间隔对应于90°相移。每个幅度测量都可以被称为抽头。

4抽头方法的一个潜在问题是对象在测量期间可能移动。因此,针对连续抽头检测到的幅度样本可以对应于不同的实际深度,例如前景中移动对象的深度和背景的深度。例如,如果对象的运动具有垂直于相机光轴的分量,则给定像素可以在一个抽头处对应于对象的一部分,并且在下一个抽头处可能对应于背景的一部分,或反之亦然。图像中靠近对象边缘的像素可能会出现此问题,并且通常会导致错误和不期望的深度。如果对象仅远离或朝向相机移动,由于相机感知的拉伸或收缩,也会发生这种情况,因此边缘附近的像素也可能在对象和背景之间发生变化。这种效果在文献中也被称为“飞行像素”。

EP 2 966 475 A1公开了一种对来自场景的TOF数据进行分箱(binning)的方法,以用于提高TOF测量的准确性并降低其中的噪声,TOF数据包括相位数据和置信度数据。根据该方法,通过用多个调制的信号照射场景来获取多个TOF数据,并且每个调制的信号分别与由相位和置信度数据定义的向量相关联。将多个向量相加以获得分箱的向量,并且对分箱的向量的相位和置信度数据进行处理以获得场景的深度数据。根据描述,“置信度”对应于反射信号的幅度。

DE 10 2015 195 161 A1公开了一种用于检测目标空间中对象的运动的装置,其中该对象位于距图像捕获设备一定距离处,该图像捕获设备被配置为:测量该距离并提供指示该距离的传感器信号,如果对象静止,则传感器信号可以在包括奇次谐波的分解中进行分解。该装置包括:确定电路,其被配置为:接收传感器信号并基于传感器信号的分解的至少一个偶次谐波生成至少一个运动信号;以及检测电路,其被配置为:基于至少一个运动信号来检测对象的运动并提供指示对象的运动的检测信号。图像捕获设备被配置为:捕获包括多个像素的图形,并且其中,确定电路被配置为:在不依赖于多个像素的相邻像素的情况下,接收传感器信号并且生成针对多个像素中的每个像素的运动信号。

发明目的

因此,本发明的目的是利用飞行时间相机来实现准确且有效的深度测量。

该问题通过根据权利要求1所述的方法来解决。

发明内容

本发明提供了一种利用使用幅度调制的连续光的飞行时间相机来进行深度测量的方法。本文中的深度测量当然是指测量距相机的距离,以便得到3D图像。使用幅度调制的连续光的飞行时间(TOF)相机的原理是众所周知的,并且已经在上文进行了解释。虽然术语“光”可以指可见光,但是应当理解也可以使用红外光或紫外光。

在第一步骤中,该方法包括:针对相机的传感器阵列的多个像素中的每个像素,以高于幅度调制的连续光的调制频率的采样频率获取至少一个样本序列,样本序列包括至少四个幅度样本。像素具体而言可以是锁定像素。传感器阵列包括多个(通常在几百到几千之间)像素,通常以二维图案布置,但也可以设想一维布置。样本序列包括至少四个幅度样本,该幅度样本以高于幅度调制的连续光的调制频率的采样频率而被采样。特别地,采样频率可以是调制频率的整数倍。幅度样本通常与发射信号和接收信号之间的相关函数的幅度相对应。对于与静止对象的表面相对应的像素,相关函数应该是正弦的,即它应该对应于具有正常非零相移的正弦函数(或为余弦函数,分别地)。

在该方法的另一个步骤中,针对每个像素的每个样本序列确定置信度值,该置信度值指示幅度样本与幅度的正弦时间演变的对应程度。确定置信度值的方法有多种,下面将讨论其中的一些方法。通常,置信度被如下定义,使得高对应度的置信度高,而低对应度的置信度低。该方法步骤基于以下假设:对于静止对象,幅度样本应对应于正弦函数。根据诸如测量误差等各种因素,幅度样本通常不会与正弦时间演变完全对应,而仅在一定程度上对应。然而,如果对象在运动并且一些幅度样本实际上对应于从对象的表面接收到的信号,而另一些对应于从背景接收到的信号,则幅度样本和任何正弦函数之间通常会有很大的差异。换句话说,幅度样本将根本不会对应于正弦时间演变,这会影响置信度值。需要说明的是,一般至少需要4个幅度样本来确定对应程度。正弦函数可以用4个参数来描述,即幅度、偏移、相位和频率。由于在当前情况下正弦函数的频率是已知的,因此保留3个参数。因此,总是有可能找到针对3个(或更少)幅度样本的“拟合”正弦函数。另一方面,如果有4个或更多幅度样本,则可以确定与正弦时间演变的任何偏差。

在该方法的另一个步骤中,针对多个分箱区域中的每个分箱区域,每个分箱区域包括多个像素,基于来自分箱区域的像素的样本序列的幅度样本来确定分箱的深度值,其中,样本序列对分箱的深度值的贡献取决于其置信度值。该方法步骤可以被描述为“分箱步骤”。这样的分箱在本领域中是已知的,但是本发明的方法应用可以被称为“智能分箱”等的先前未知的变体。定义多个分箱区域,其中每个分箱区域包括多个像素。每个分箱区域可以是矩形的,例如包括m x n像素或n x n像素。也可以将分箱区域称为像素组。通常,每个分箱区域是相干的,即它对应于传感器阵列的相干区域。可以想象,两个分箱区域重叠,使得给定的像素属于多于一个分箱区域。但是,通常情况下,不同的分箱区域是分开的。在这种情况下,分箱区域共同可以被视为低分辨率图像的单位或“像素”,而传感器阵列的像素对应于高分辨率图像。来自分箱区域的像素样本序列的幅度样本对该分箱区域的分箱的深度值有贡献。也可以说对来自分箱区域的像素的样本序列的幅度样本的信息进行组合以获得分箱的深度值。如何计算分箱的深度值有多种可能性,其中一些将在下面进一步讨论。

如果针对单个像素(和单个样本序列)确定深度值,这将通过确定正弦函数的相对相位来完成,该相位将对应于深度值。然而,如上所述,由于正弦函数的相位是根据幅度样本确定的,如果一些幅度样本对应于对象而一些对应于其背景,则相位信息会受到影响。从具有这种幅度样本的样本序列得出的任何深度信息大多是无用的。然而,根据本发明的方法,针对每个像素的每个样本序列确定置信度值,由此可以评估相位信息的可靠性。然后,当使用来自给定像素的幅度样本的信息来确定分箱的深度值时,并非分箱区域中的所有像素的样本序列都被同等对待(如本领域已知的分箱方法),而是样本序列对分箱的深度值的贡献取决于其置信度值。定性地,具有高置信度值的样本序列的贡献通常大于具有低置信度值的样本序列的贡献。后者的贡献甚至可能为零,即这个样本序列可以被完全忽略。

本发明的概念允许增加分箱的深度值的可靠性,因为幅度样本的“被破坏的”序列被较少考虑或甚至根本不考虑。此外,置信度值可以单独基于各个像素的信息来确定,即不需要考虑其他像素。正如下面将变得明显的那样,也可以用最少的处理功率和存储器来计算置信度值。这也意味着即使通过简单、低成本的处理单元也可以实时进行计算。

优选地,该方法包括:以比幅度调制的连续光的调制频率高四倍的采样频率获取四个幅度样本。这对应于公知的4抽头方法。幅度为A且相位为φ的正弦波调制接收信号可以用二维向量表示:

r(A,φ)=(A·cosφ,A·sinφ)=(d13,d02) (公式1)

其中,d02=A0-A2,以下称为第一差,以及d13=A1-A3,以下称为第二差,它们分别是两个幅度样本Ak、k=0…3的两两差值。其中第一幅度样本A0对应于0°的相位角,第二幅度样本A1对应于90°的相位角,第三幅度样本A2对应于180°的相位角,而第四幅度样本A3对应于270°的相位角。因此,接收到的信号的幅度和相位可以计算为

φ=atan2(d02,d13) (公式2)

虽然信号的幅度A与接收到的光子数量成比例,但相位φ与由相对应的像素看到的对象的深度D成比例。

其中D是测量的距离相机的对象深度,c是光速,而fmod是信号的调制频率。

因此,深度可以计算为

关于可以如何定义置信度值有多种选择。置信度值C的一种可能定义如下:

其中幅度A根据公式3计算,但可以近似为

在一个简单的实施例中,针对每个像素仅获取一个样本序列,因此对于置信度值并且对于分箱的深度值,可以仅考虑该样本序列的幅度样本。在这样的实施例中,还可以说针对每个像素确定了(单个)置信度值,因为每个像素对应于单个样本序列。在另一实施例中,该方法包括:针对至少一个像素,获取多个样本序列。这些也可以被称为多次曝光。一方面,必须按顺序执行的多次曝光增加了某些像素受对象的移动影响的可能性。另一方面,考虑多个样本序列可以有助于减少噪声或其他测量误差的影响。具体而言,可以针对不同的样本序列执行抽头测量以用于获取具有不同积分时间的幅度样本,即每个样本序列对应于不同的积分时间。除了更改曝光时间外,还可以更改其他参数,例如调制频率。

即使针对给定的像素获取多个样本序列,每个置信度值也可以通过个体样本序列的幅度样本彼此之间的关系来确定。换言之,在不考虑其他样本序列(如果存在)的幅度样本的情况下考虑样本序列的幅度样本。公式6是这种实施例的示例。在大多数情况下,这种方式是可靠的,因为单个样本序列的幅度样本不受例如不同的曝光时间的影响。但是应该注意,如果已经确定了多个样本序列,则应该针对每个样本序列确定单独的置信度值,例如因为针对给定像素的一个样本序列可能不受对象运动的影响,而另一个样本序列受到影响,因此产生不可靠的数据。

对于不同的样本序列可以如何基于它们各自的置信度值对分箱的深度值做出贡献存在多种可能性。例如,可以存在连续范围的贡献或加权因子。另一种可能性可以被称为“二进制”分类。在这种情况下,该方法包括:如果置信度值满足预定义标准则将每个样本序列分类为有效,否则分类为无效。换句话说,每个样本序列只能为有效或无效。在这种情况下,该方法还包括:仅当像素有效时,才使用样本序列的幅度样本来确定分箱的深度值。换句话说,如果样本序列被认为是无效的,则该样本序列的幅度样本将被完全忽略,并且不影响分箱的深度值。如果针对每个像素只考虑一个样本序列,也可以说每个像素分别被分类为有效或无效。如果仅单独考虑特定曝光的样本序列,情况也是如此。

具体而言,可以基于置信度值与第一阈值的关系来对样本序列进行分类。换言之,确定置信度值是高于还是低于第一阈值,并取决于结果将样本序列分类为有效或无效。通常,对于置信度是程度高的,则置信度值高,并且如果置信度值高于第一阈值,则样本序列被分类为有效。第一阈值(也可以被称为运动参数)通常是预定义的。其可以被估计、计算或可以通过校准来确定,例如在没有移动对象的情况下使用静止的风景。从图形上讲,这种方式可以用置信度掩码来描述,该置信度掩码掩蔽了具有的置信度值低于第一阈值的所有像素。

根据一个实施例,分箱的深度值是基于来自分箱区域的像素的样本序列的幅度样本的线性组合确定的,其中,每个样本序列对线性组合的贡献取决于各个样本序列的置信度值。换句话说,在线性组合中,向每个幅度样本分配系数,该系数取决于相应样本序列的置信度值。一般来说,当然,对于高置信度值,系数会较高,而对于低置信度值,系数会较低。具体而言,如果置信度值高于第一阈值,则系数可以是1,而如果置信度值低于第一阈值,则系数可以为零。在这种情况下,线性组合对应于所有有效样本序列的总和,而忽略所有无效样本序列。例如,可以针对每个有效样本序列定义上述第一和第二差值,并且然后定义与所有有效样本序列的各个差值之和相对应的“分箱的”第一和第二差值。可以理解,这等效于首先取幅度样本的总和来确定“分箱的”幅度样本,而然后计算分箱的第一和第二差值。当已确定分箱的插值时,可根据公式2来计算“分箱的”相位。对差值进行求和也可以看作是向量加法,第一差值和第二差值是向量的分量,这个向量加法得到分箱的向量。因此,这种方式也可以被称为“向量分箱”或“加权的向量分箱”。

根据另一实施例,分箱的深度值Db是通过对来自分箱区域的像素的样本序列的像素深度值进行平均来确定的,其中,每个像素深度值的权重取决于相应像素的相应样本序列的置信度值,并且其中,像素深度值是基于像素的样本序列的幅度样本来确定的。换言之,像素深度值D是针对分箱区域中像素的每个样本序列单独确定的,或者可能仅针对分箱区域中像素的有效样本序列而确定。可以确定这些像素深度值D,例如使用公式2和公式5来确定。然后,对像素深度值进行平均以确定分箱的深度值Db,但以加权的方式,使得权重取决于相应像素的相应样本序列的置信度值。这种方式也可以被称为“深度分箱”或“加权的深度分箱”。同样,对于高置信度值,像素深度值的权重较高,而对于低置信度值,像素深度值的权重较低。具体而言,如果置信度值低于第一阈值,即仅考虑针对有效样本序列的像素深度值,则像素深度值的权重可以为零。如果针对每个像素只存在一个样本序列,计算可以如下进行:

在这样的实施例中,分箱的深度值与所有有效像素(即具有有效样本序列的像素)的像素深度值的(算术)平均值相对应。如果针对每个像素获取多个样本序列,则分别针对每个样本序列确定像素深度值,并且必须修改公式8以对所有像素的所有有效样本序列(或所有积分时间的所有有效像素)求平均。

在一个实施例中,该方法包括:确定像素的样本序列的第一幅度样本和第三幅度样本之间的第一差值,并且向第一组分配具有正第一差值的样本序列,以及向第二组分配具有负第一差值的样本序列。具体而言,这可以仅指有效样本序列,而无效样本序列不包括在上述两组中的任一组中。已经在上文提到了第一差值d02=A0-A2。当将第一差值和第二差值视为向量的分量时,具有正第一差值的任意两个向量的相位差小于180°,而具有负第一差值的任意两个向量的相位差小于180°。当来自第一组(或第二组,分别地)的两个向量相加时,所得向量的相位介于相加向量的相位之间。相应地,根据公式5,与结果向量相对应的深度也在与两个相加向量相对应的深度之间。需要说明的是,一般而言,可能所有样本序列都具有正的第一差值,或者也可以是所有样本序列分别具有负的第一差值。如果是这种情况,当然没有必要将样本序列划分成组。

此外,该方法可以包括:定义向量,其具有第二幅度样本和第四幅度样本之间的第二差值作为第一分量,以及具有第一差值作为第二分量。上文已经解释了第二差值d13=A1-A3。此外,该方法可以包括:定义第一组向量rP=[xP,yP],它是第一组中的向量基于相应样本序列的置信度值的线性组合,以及第二组向量rM=[xM,yM],它是第二组中的向量基于相应样本序列的置信度值的线性组合。在第一组向量中,基于各个样本序列的置信度值对第一组中的向量进行线性组合。换句话说,向量的系数或权重取决于置信度值。具体而言,如果置信度值低于第一阈值,则权重可以为零,并且如果置信度值高于第一阈值,则权重可以为一,在这种情况下,仅添加有效样本序列。这同样适用于第二组中的向量的线性组合。更具体地,每个组向量可以是相应组中所有向量的总和,在这种情况下,第一和第二组向量rP,rM的分量计算如下,其中,多个积分时间的总和是可选的:

在这里,每个样本序列对应单独的积分时间。在公式9a-公式9d中,假设针对每个像素有相同数量的样本序列(或积分时间,分别地)。“有效像素”的总和应理解为具有有效第It(It=0…n)个样本序列的所有像素的总和。替换地,可以对所有像素和相应像素的所有有效样本序列求和。类似于公式2,组向量的相位计算如下:

以及

如上所述,与第一组向量相对应的相位(因此深度)在第一组的各个向量的区间内。并且,与第二组向量相对应的相位在第二组的各个向量的区间内。在另一步骤中,该方法包括:基于第二组向量和第一组向量之间的相位差来确定分箱的深度值。换句话说,确定第二组向量和第一组向量之间的相位差(或取决于相位差的量)以及分箱的深度值的确定(或计算)取决于该相位差。通常,假设两个相位都在0°和360°之间。

根据一个实施例,该方法还包括:如果相位差低于第二阈值,则基于第一组向量和第二组向量二者来确定分箱的深度值;以及如果相位差高于第二阈值,则仅基于第一组向量和第二组向量中的一个来确定分箱的深度值。与第一阈值一样,也可以通过使用静止场景的校准来估计、计算或确定第二阈值。如果相位差低于第二阈值,这通常对应于所有(或大多数)第一差值具有相对较低的绝对值的情况,其中,第一差值中的一些为正,并且一些为负,而第二差值大部分为负。在这种情况下,可以假设第一组向量和第二组向量仅在较小程度上不同,并且可以将第一组向量与第二组向量相加,

rb=rP+rM (公式10)

之后,可以再次使用公式5,基于得到的分箱的向量来确定分箱的深度值。另一方面,如果相位差高于第二阈值,则这表明第一组向量和第二组向量可以对应于不同的深度,例如前景对象和背景的深度。在这种情况下,完全丢弃一个组向量并仅基于另一个向量来确定分箱的深度值更为合适。

根据一个实施例,第二阈值为180°。这可以看作是用于避免在由第一组向量和第二组向量给定的区间之外的不切实际的深度值的最小条件。然而,第二阈值可以更小,例如小于90°或小于60°。如果第二阈值为180°,则不必明确计算相位差。相反,可以示出在以下条件下相位差低于180°:

xMyP<xPyM

测试该条件只需要两次乘法,而因此可以节省处理功率。替换地,也可以基于第一组向量和第二组向量的相位来明确地计算相位差。

如果相位差高于第二阈值,则存在多种可能的标准来确定第一和第二组向量中的哪一个被认为更可靠。一般来说,可以合理假设:如果一个组向量基于更多数量的有效样本序列,则该组向量是更可靠的。因此,如果相位差高于第二阈值,则分箱的深度值可以是基于具有更多有效样本序列的组的组向量来确定的。换句话说,如果第一组具有的有效样本序列比第二组具有的更多,则基于第一组向量来确定分箱的深度值,反之亦然:

其中,NP,NM分别是第一组和第二组中有效样本序列的数量。

根据另一个实施例,该方法包括:如果两个组向量的第一分量xP,xM都为负,则基于第一组向量和第二组向量二者来确定分箱的深度值;以及如果至少一个第一分量为正,则仅基于第一组向量和第二组向量中的一个来确定分箱的深度值。如果第一组向量的第一分量和第二组向量的第一分量都为负,则意味着:组向量位于第二和第三象限中,对应于接近180°的相位和接近模糊深度一半的对象深度。但是,如果至少一个第一分量为正,则向量要么位于相反的象限(分别为第一和第三象限或第二和第四象限)中,要么位于第一和第四象限中,其中,一个向量对应于接近360°的相位或接近模糊深度的深度。

附图说明

本发明的进一步细节和优点将从下文参考附图的非限制性实施例的详细描述中显而易见,其中:

图1是可用于本发明方法的TOF相机和对象的示意图;

图2是示出函数和四个幅度样本的时间演变图;

图3是向量图;

图4是图示传感器阵列上的幅度值的图;

图5是示出函数和多个幅度样本的时间演变的另一个图;

图6是图示根据现有技术的深度测量的结果的第一图;

图7是图示根据现有技术的深度测量的结果的第二图;

图8是图示向量相加的向量图;

图9是示出根据现有技术的深度测量的结果的第三图;

图10是图示本发明方法的第一实施例的流程图。

图11是图示二元置信度掩码的图;

图12是图示二元置信度掩码的构造和该置信度掩码的应用的图;

图13是图示向量相加的另一个向量图;

图14是图示两组向量的位置的向量图;

图15是示出根据本发明的深度测量的结果的第一图;

图16是示出根据本发明的深度测量的结果的第二图;

图17是示出根据本发明的深度测量的结果的第三图;

图18是图示本发明方法的第二实施例的流程图;以及

图19是示出根据本发明的深度测量的结果的第四图。

具体实施方式

图1示意性地示出了适用于使用幅度调制的连续光进行深度测量的TOF相机1。它包括具有多个(例如几千或几万)像素3的矩形传感器阵列2。此外,它可以包括存储器和处理单元,为简单起见未示出。相机1被配置为:使用一个或多个光发射器5发射幅度调制的连续光10。光10被相机1视野中的3D对象20或景物反射,并且反射光11被传感器阵列2的像素3接收。具有相位延迟τ的原始调制函数s(t)与接收函数q(t)进行相关,从而产生相关函数c(τ)。接收函数q(t)的幅度以比光10的调制频率fmod高四倍的频率被采样。换句话说,四个幅度样本A0…3,也被称为抽头,用于取回调制光的相位φ,如图2所示。每四个幅度样本A0…3是针对相应像素3的样本序列的一部分。

具有幅度A和相位φ的正弦信号可以由2D向量表示,该向量可以从抽头测量中确定的4个幅度样本中确定,即

r(A,φ)=(A·cosφ,A·sinφ)=(d13,d02) (公式1)

其中,d02=A0-A2,以下称为第一差,以及d13=A1-A3,以下称为第二差,它们分别是两个幅度样本Ak,k=0…3的两两差值。因此,接收到的信号的幅度和相位可以计算为

φ=atan2(d02,d13) (公式2)

虽然信号的幅度A与接收到的光子数量成比例,但相位φ与由相对应的像素看到的对象的深度D成比例。

其中D是像素深度值,即距相机的像素的距离,c是光速,而fmod是信号的调制频率。因此,深度可以计算为

图3是当调制频率为fmod=20MHz时从深度D=2m的对象接收到的信号的向量表示。

然而,运动伪影可能沿着在相机1的场景中移动的对象20的边缘发生。由于随后执行抽头测量,因此靠近对象边缘的像素3在一个幅度样本的获取期间可能“看到”对象表面,而在随后的获取中它可能看到背景。

图4通过示例示出了运动伪影的发生。在获取每个幅度样本期间,具有‘O’形状的对象20向上移动1个像素且向左移动1个像素。灰度值表示对象存在于相应像素情况下的抽头测量次数。黑色像素代表对象存在情况下的零抽头测量,而白色像素代表对象存在情况下的(四个中的)四个抽头测量。右侧是完整图像,左侧是围绕‘O’左上角外边缘的缩放。

图5图示了由对象20的运动引入的误差。取决于像素3在获取期间看到的内容,测量的幅度样本Ak位于前景或背景信号的正弦互相关曲线上。如果将一个像素3的4个后续幅度样本Ak相加,则可以识别出沿对象20边缘的模糊效果。在4x4像素3的邻域中,可以找到在所有抽头中“看到”前景对象20的像素3以及部分看到前景对象20或背景的像素3。如果随后执行了具有不同积分时间的获取,则对应于额外的样本序列,它们也可能对应于不同的深度。

由于前景和背景对象的深度和反射率不同,幅度样本Ak可能会有很大差异。因此,根据公式2和公式5计算的相位和深度可能是错误的。结果由图6中的图说明,它是在7m背景前2m深度处的“O”形目标的高分辨率深度图像,在水平和垂直两个方向上以每次抽头获取移位1个像素的方式移动。由于运动伪影,沿对象边缘,计算的深度在1.65m和5.41m之间变化。应该注意的是,测量的深度不仅位于前景和背景深度之间,也可能在这个深度范围之外。相对应的像素可以被称为飞行像素。

根据现有技术,有两种主要的方式来缓解该问题,这两种方法都使用分箱方法。多个像素,例如4x4像素被视为分箱区域,针对其确定单个深度值。在第一种方式中,对分箱区域中所有像素的幅度样本Ak进行相加,并使用公式2(用总和而不是个体幅度样本Ak)和公式5来计算单个深度值。这种方法可以被称为“抽头分箱”。其结果如图7所示。人们认识到,在位于0.31m和6.93m之间的范围内测量的深度存在异常值。换句话说,仍然存在对象20和背景的深度范围之外的深度值,并且与图6的高分辨率图像相比,甚至增加了飞行像素效应。这种增加的一个原因可以通过图8的向量图来理解。增加幅度样本Ak对应于如图8针对表示对象20的深度的第一向量和表示背景的第二向量所示的向量相加。由于两个向量的相位相差超过180°,因此得到的向量的相位小于任一向量的相位。这导致深度值超出深度范围。

根据另一种方式,针对分箱区域中的每个单独像素确定像素深度值并且对这些像素深度值求平均以确定针对分箱区域的深度值。这种方式可以被称为“像素分箱”。结果如图9所示。平均导致位于1.81m到6.9m范围内的深度值变得模糊。

通过本发明方法减少或消除了上述问题。图10是图示本发明方法的第一实施例的流程图。

在该方法开始后,在100处选择分箱区域4。这可以是例如包括4x 4像素3的区域(也参见图12)。接下来,在110处选择分箱区域4内的像素3。在120处,针对该像素3的样本序列确定幅度样本Ak。在130处,基于幅度样本Ak来计算置信度值C。针对相应像素3的每个样本序列计算个体置信度值C,即如果只有一个样本序列,则针对每个像素计算一个置信度值C。置信度值C的一种可能定义如下:

其中幅度A根据公式3计算,但可以近似为

根据这个定义,置信度值C总是在0和1之间的范围内,其中,最高可能值1代表完美的正弦函数。在140处,将该置信度值C与第一阈值Cmin进行比较,该第一阈值可以通过使用静止场景的校准来计算、估计或确定。在以下示例中,第一阈值Cmin可以是0.25。第一阈值Cmin也可以被称为“运动参数”,因为它可能适合区分受对象运动影响的样本序列和不受对象运动影响的样本序列。如果置信度值C小于第一阈值Cmin,则相应的样本序列在190处被分类为无效,并且基本不作进一步考虑。另一方面,如果置信度值C大于第一阈值Cmin,则相应的样本序列在150处被分类为有效。幅度值或第一和第二差值d02,d13分别被视为向量的第二和第一分量,被保留以供进一步处理。

这个过程可以被看作是二元置信度掩码的创建,这在图11中以图形方式说明。图11的上半部分与图4相对应,而下半部分示出了相对应的置信度掩码,其中,左半部分是对象20边缘附近部分的放大图。黑色表示被视为无效的像素,而白色表示具有被视为有效的样本序列的像素。人们认识到,抽头模糊的区域被置信度掩码所掩蔽。

图12进一步图示了用于4×4像素3的分箱区域4的置信度掩码的构造和分箱过程,其中为了简单起见,假设了针对每个像素3的单个样本序列。首先,如a)所示,针对每个像素确定单独的幅度样本(其中,不同的阴影分别代表采样编号或时间点)。然后,如b)所示,针对每个像素确定置信度值(其中,暗色调代表高置信度值)。在c)处显示置信度掩码,其中,黑色代表无效像素(或样本序列,分别地),而白色代表有效像素。将此置信度掩码与个体幅度样本一起使用有效地产生针对整个分箱区域4的分箱的幅度样本,如d)所示,其中不同的阴影再次分别表示采样编号或时间点)。如果考虑与若干积分时间相对应的若干本序列,则可以针对每个积分时间构建置信度掩码。

在160处,确定第一差值d02的符号。如果符号为正,则在170处将样本序列及其向量分配到第一组,并且如果符号为负,则在180处将样本序列及其向量分配到第二组。如虚线箭头所示,在该方法的简化版本中也可以跳过步骤160、170和180。

对分箱区域4中的所有像素3以及在适用的情况下对每个像素3的所有样本序列重复目前为止提到的步骤。当在200处确定最后一个像素3已被处理时,该方法在210处继续通过分别将第一组和第二组中的向量相加来计算第一组向量rP=[xP,yP]和第二组向量rM=[xM,yM]。换句话说,对具有正第一差值d02的所有向量进行求和,并且对具有负第一差值d02的所有向量进行求和。因此,第一和第二组向量rP,rM的分量计算如下,其中多个积分时间上的总和是可选的:

需要注意的是,在第一组和第二组中的任一组中,仅对有效样本序列的向量进行求和,而针对分箱过程,忽略无效样本序列。第一组中任意两个向量的相位相差小于180°,因此这些向量的相加不会导致飞行像素。这同样适用于第二组中的向量。求和向量的相位相差小于180°这一事实保证了得到的组向量不受分箱效应的影响。

在220处,计算第二组向量和第一组向量的相位差ΔΦ(假设这两个相位都在0°和360°之间)并将其与第二阈值Φmax进行比较。特别地,第二阈值Φmax可以等于180°。如果相位差ΔΦ较小,如图13的示例所示,则第一和第二组向量rP,rM在230处简单相加以计算分箱的向量rb=[xb,yb],即:

rb=rP+rM (公式10)

如果相位差ΔΦ较大,如图14的示例所示,则这可以表明这些组分别对应于背景和前景对象的像素3。无论哪种方式,将两个组向量rP,rM相加都会导致飞行像素。由于这些原因,在240处选择一个组向量rP,rM作为分箱的向量rb,即较大组(即具有的有效样本序列的数量较多的组)的组向量:

其中,NP,NM分别是第一组和第二组中有效样本序列的数量。最后,使用公式2和公式5基于分箱的向量rb来确定分箱的深度值Db

如果记录了若干具有若干积分时间It=1,2,..n的样本序列,则分箱值,例如分量xb,yb可以标准化为:

其中,NIt是在特定积分时间期间具有有效样本序列的像素3的数量,而TIt是积分时间的长度。这产生不显示人为跳跃的归一化的幅度,因为它独立于分箱中考虑的像素3的数量,并且因此允许应用标准图像处理方法,例如对分箱抽头或幅度进行杂散光补偿。

在虚线所指示的方法的简化版本中,在250处,将有效像素3的所有向量相加以确定分箱的向量rb

之后,基于分箱的向量rb确定分箱的深度值Db。图15示出了该简化版本的结果。与图7相比,可以看到显著的改进,但仍然存在具有的深度超出分箱的像素3的深度范围的异常“飞行像素”。如果将向量分配给第一组和第二组并分别处理,这个问题就会减少,如图16中可见。虽然图16示出了单个积分时间的结果,图17示出了两个积分时间的结果,但第一积分时间比第二积分时间长4倍。在这种情况下,飞行像素的数量减少到几乎可以忽略不计的量。

有两种可能的替代方法来检查220处相位差ΔΦ的不等式。首先,可以检查以下关系:

xMyP<xPyM (公式13)

如果是,则该方法在230处继续,如果不是,则在240处继续。该条件与向量(x,y)的斜率有关,其与y/x成比例。对于区分两个向量之间的夹角是小于还是大于180°的关键情况,其中一个向量必须在象限1中,而另一个在象限中,或者一个向量位于象限2中,而另一个在象限4中。对于任何其他情况,区别是微不足道的。如果一个向量在象限1中,而另一个向量在象限4中,则两个向量之间的夹角明显大于180°。如果一个向量在象限2中,另一个向量在象限3中,则两个向量之间的夹角明显小于180°。

其次,可以判断xP和xM二者是否都是负数,这意味着对象深度接近模糊深度的一半。如果是,则该方法在230处继续,如果不是,则在240处继续。

图18是图示本发明方法的第二实施例的流程图。步骤100、110、120、130和140与第一实施例相同,而为简洁起见不再赘述。如果样本序列在150处被分类为有效,则在155处根据公式2和公式5确定像素深度值D。在处理完分箱区域4中的所有像素3后,通过对像素深度值D求平均值来确定分箱的深度值Db

如果针对每个像素3获取多个样本序列,则分别针对每个样本序列确定像素深度值,并且必须修改公式8以对所有像素3的所有有效样本序列或所有积分时间的所有有效像素3求平均。计算分箱的深度值Db的结果为低分辨率深度图像,其中分箱的深度值Db表示各个分箱区域4中有效像素的算术均值。图19示出了本实施例中计算的深度图像的示例。与示出了没有区分有效像素和无效像素的平均过程的结果的图9相比,减少了飞行像素的影响。然而,应当注意,本发明方法的该第二实施例需要计算全高分辨率图像上每个有效像素(和每个样本序列,在适用的情况下)的像素深度值D,这会导致计算工作量增加和/或存储器要求增加。

附图标记列表

1 TOF相机

2 传感器阵列

3 像素

4 分箱区域

5 存储器

10 光

11 反射光

20 对象

26页详细技术资料下载
上一篇:一种医用注射器针头装配设备
下一篇:对人-产品交互的无线测量

网友询问留言

已有0条留言

还没有人留言评论。精彩留言会获得点赞!

精彩留言,会给你点赞!