用于确定到场景的距离的方法和设备与流程

文档序号:30711991发布日期:2022-07-10 20:23阅读:128来源:国知局
用于确定到场景的距离的方法和设备与流程
用于确定到场景的距离的方法和设备
1.公开领域
2.本公开涉及一种用于确定到场景的距离的系统。更具体而言,本发明涉及一种使用激光来照亮场景并检测反射激光的基于飞行时间的感测方法。本发明还涉及用于检测到场景中的一个或多个对象的距离的lidar设备。


背景技术:

3.基于主动照明的用于确定场景距离的设备也称为lidar(光检测和测距)设备。lidar设备通过用激光照亮场景并通过在检测器中检测反射激光来测量到场景的距离,检测器通常位于发射激光的激光源附近。因此,激光的发射和对反射激光的检测之间的时间间隔与到反射该激光的场景对象或场景部分的距离的两倍成正比。
4.大多数已知的lidar设备都使用直接tof(dtof)检测方法。这些系统包括在纳秒脉冲范围内操作的强大脉冲激光器、扫描脉冲激光束的机械扫描系统、以及脉冲检测器。这一类型的系统目前可从包括加利福尼亚州摩根山velodyne lidar在内的供应商处获得。作为现有系统的示例,velodyne hdl-64e在机械旋转结构中以每秒5到15转的速度使用64个高功率激光器和64个雪崩二极管检测器。
5.这些dtof系统具有一些缺点。例如,这些系统需要具有无法用当前可用的半导体激光器获得的过高功率电平的激光器,当前可用的半导体激光器的功率电平低几个数量级。另外,用于扫描目的的机械旋转元件的使用进一步限制了此类系统的小型化、可靠性和成本降低的前景。
6.lidar设备的紧凑性是汽车行业中应用的重要因素,其中lidar设备例如将耦合到汽车的前挡风玻璃或耦合到汽车的保险杠上。事实上,lidar设备是开发自动驾驶或驾驶员辅助系统的关键因素。在这一上下文中,lidar设备被用于检测障碍物,诸如车辆环境中的其他车辆或对象。
7.在wo2017/068199中,提出了一种使用半导体激光器和基于cmos的多阱像素检测器的lidar设备。这允许制造紧凑型lidar设备,其中照明系统和检测系统一起放置在紧凑的壳体中。wo2017/068199中公开的设备使用投影仪,该投影仪被配置成用于用由以光斑图案同时发射的多个在空间上分开的脉冲激光形成的照明图案来照亮场景。激光脉冲作为脉冲链发射,即以给定的脉冲频率发射脉冲序列。基于cmos的多阱像素检测器检测代表场景所反射的在空间上分开的脉冲或在空间上分开的脉冲的一部分的反射激光的光斑。以此方式,可以在不使用任何旋转或机械扫描元件的情况下获得场景的深度图。
8.wo2017/068199中使用的照明和检测方法的原理在如图1上示出,其中所发射的激光脉冲和所反射的脉冲分别用附图标记11和12标识。该图1仅是示意图,且与所发射的脉冲11的振幅相比,所反射的脉冲12的振幅通常较小。这一方法也被称为距离选通方法。检测器被配置成用于检测代表在第一时间窗口tw1期间获得的反射激光的第一量的第一电荷量和代表在第二时间窗口tw2期间获得的反射激光的第二量的第二电荷量。第二时间窗口发生在第一时间窗口之后,并且这两个时间窗口的宽度等于激光脉冲的脉宽pw。处理装置允许
基于第一电荷量和第二电荷量来计算到场景的对象的距离。
9.在wo2017/068199中,场景是用脉冲序列来照亮的。以此方式,可以在第一tw1和第二tw2时间窗口中累积足够的电荷,以获得允许计算到场景的对象的距离的信噪比。用脉冲序列照亮场景、累积反射激光、读出电荷以及基于所累积的电荷来计算距离通常称为帧或帧测量。然而,当执行单个帧测量时,获得的距离的精度(也称为时间精度或时间误差)较低。事实上,当生成多个帧时,逐帧的测得距离存在很大差异(spread)。换言之,时间误差(即测量分布的sigma值σ)很大。因此,总是拍摄多个帧,并确定平均对象距离。以此方式,与单个帧测量的误差相比,平均对象距离值的误差用因子来减小,其中nf是帧数。对于经完美校准的lidar设备,从多个帧获得的所计算的平均距离等于标准偏差所确定的置信区间内的实际距离。
10.当与上述dtof系统相比,wo2017/068199中公开的lidar设备具有许多优点。事实上,此类固态lidar设备紧凑、稳健,并且可以以成本高效的方式生产。然而,这些类型的lidar设备也有缺点,即需要拍摄大量帧以获得可接受的空间准确度的事实。这大大缩短了lidar检测器的总体响应时间,并限制了在给定时间段内可以测量的点的数量。
11.因此,需要改进此类紧凑型lidar设备,更确切而言,需要使这些类型的lidar设备可在广泛的距离范围内操作,包括等于或大于100米的距离,并以高精准度进行距离测量。
12.概要
13.本发明的目的是提供一种用于确定到场景中的一个或多个对象的距离的方法,更具体而言,用于以高精准度确定对象距离,以使得即使使用单个帧测量,也可以在有限的帧数下获得可接受的空间准确度。换言之,目的是减小上述时间误差σ。另一目的是提供一种包括检测和处理装置的lidar设备,该检测和处理装置用于以高精准度确定对象距离,使得能够极大地减少要拍摄的帧数或甚至限制到单个帧测量。
14.本发明在所附独立权利要求中定义。从属权利要求定义了有利的实施例。
15.根据本发明的第一方面,一种用于通过使用脉冲激光照亮场景并检测与激光飞行时间相关的反射激光来确定到场景中的一个或多个对象的距离的方法,并且其中场景中的一个或多个对象位于最小距离和最大距离之间的距离范围内,并且其中0≤d
min
≤0.6
×dmax
,优选地0≤d
min
≤0.4
×dmax
,更优选地0≤d
min
≤0.2
×dmax
,其中d
min
和d
max
分别是限定距离范围的最小距离和最大距离。
16.根据本发明的方法包括以下步骤:
17.a)确定脉冲激光的脉宽,该脉宽小于与最大距离相关联的最大飞行时间,以使得
18.pw=(tof
max

t
dl
)/n,其中tof
max
=2
×dmax
/c,
19.其中pw是脉宽,d
max
是最大距离,c是光速,tof
max
是最大飞行时间,t
dl
是预定延迟时间窗口,其中0≤t
dl
≤2
×dmin
/c,且n是预定义脉宽缩减因子,
20.b)确定脉冲频率,使得f
p
≤1/((n+1)
×
pw+t
dl
),其中f
p
是脉冲频率,并且其中f
p
=1/p
p
,其中p
p
是脉冲周期,
21.c)确定用于在脉冲激光的脉冲周期内检测反射激光的检测时间段,并且其中td=(n+1)
×
pw,其中td是检测时间段,
22.以及将检测时间段拆分成数个连贯检测时间窗口tw[i],其中i=1到m,且m是检测
时间窗口的数量,并且其中:
[0023]
且m=α
×
(n+1)
[0024]
其中α是表示每脉宽所选择的检测时间窗口的数量的整数,其中α≥
[0025]
1,优选地α≤10,优选地α≤5,并且其中α和n被选择成使得α
×
n≥
[0026]
2,优选地α
×
n≥3,更优选地α
×
n≥4,
[0027]
d)用由多个空间上分开的脉冲激光束形成的光斑图案来照亮场景,并且其中每一脉冲激光束包括脉冲序列,该脉冲序列具有在步骤a)中确定的脉宽和在步骤b)中确定的脉冲频率,并且其中多个脉冲激光束中的每一者的脉冲是同时发射的,
[0028]
e)在如步骤d)中的所定义地照亮场景时,在每次同时发射脉冲之后,在步骤c)中定义的检测时间段(td)的连贯检测时间窗口tw[i]期间检测反射激光的光斑,并且其中这些光斑表示被场景中的一个或多个对象反射的在空间上分开的脉冲激光束,
[0029]
以及对于诸脉冲的每一同步发射,与脉冲激光束的诸脉冲的发射同步地执行对反射激光的光斑的所述检测:
[0030]
a)如果延迟时间窗口为零,则执行所述同步以使得第一检测时间窗口tw[1]与脉冲的发射时间窗口交叠或至少部分交叠,以及
[0031]
b)如果延迟时间窗口不为零,则执行所述同步以使得第一检测时间窗tw[1]相对于脉冲的发射时间窗口延迟且延迟等于该延迟时间窗口,
[0032]
f)对于每一所检测的反射激光的光斑,通过累积在检测时间窗口tw[i]中检测到的反射激光的量来获得与检测时间窗口tw[i]相关联的曝光值,并且其中对于脉冲序列的所有脉冲执行该累积,
[0033]
g)对于每一所检测的反射激光的光斑:
[0034]
i)标识在连贯检测时间窗口tw[i]中的哪些检测时间窗口中检测到反射激光,
[0035]
ii)基于所标识的检测时间窗口来计算对象距离,并且其中该对象距离对应于检测反射激光的光斑的检测器与导致反射激光被检测到的场景中的对象或该场景中的对象的一部分之间的距离。
[0036]
有利地,通过相对于最大飞行时间tof
max
(即对应于要测量的最大距离的tof
max
)来减小激光脉冲的脉宽,并减小用于检测反射激光的检测时间窗口,距离确定的时间精度(上述定义为测量分布的sigma值σ)得到了改善,即sigma值σ降低。换言之,时间精度得到了改善。这允许减少要拍摄的帧数,或甚至允许执行单个帧测量。如果lidar设备经过完美校准并且如果仅执行单个帧测量,则时间精度确定测得距离值与实际距离值相等的概率。
[0037]
有利地,定义检测时间段和脉冲频率,使得在发射后续激光脉冲之前的在检测时间段内检测来自位于最大距离处的对象的反射激光。这避免了混叠问题。
[0038]
有利地,利用根据本发明的方法,利用具有所述的脉宽的相同激光脉冲序列和所述的检测时间窗口,可以在较宽的距离范围内(即从短距离d
min
到最大距离d
max
)检测对象。
[0039]
在检测时间段期间因变于时间来执行对反射激光的检测具有在最多两个检测时间窗口内检测到反射激光的额外优点,其中检测时间段被划分成m个连贯检测时间窗口,其中每一检测时间窗口的时间宽度基本上等于脉宽,即α=1。因此,因为m≥3,所以存在可被用作背景检测时间窗口的至少一个检测时间窗口,使得执行背景测量不需要额外的时隙。
[0040]
同样,如果α》1,则在最多α+1个检测时间窗口中检测反射激光,并且因此脉宽因子可被选择成使得检测时间窗口可供用于执行背景测量。
[0041]
有利地,由于距离测量的经改进精度,可以减少帧数,并且因此帧数可以降低并且因此改进了lidar设备的整体响应时间。
[0042]
最小距离d
min
的下限“0”必须解释为0米。
[0043]
在各实施例中,最小距离d
min
等于0米,并且因此延迟时间窗口t
dl
被设置成0纳秒。
[0044]
在各实施例中,0.005
×dmax
≤d
min
≤0.6
×dmax
。在其他实施例中,0.01
×dmax
≤d
min
≤0.6
×dmax
。在又一些实施例中,0.01
×dmax
≤d
min
≤0.4
×dmax

[0045]
在各实施例中,n和α的积等于或小于1000,优选地等于或小于500,更优选地等于或小于100。
[0046]
在各实施例中,预定义脉宽缩减因子n等于或大于3,而在其他实施例中,n等于或大于4。
[0047]
在一些实施例中,预定义脉宽缩减因子n包含在以下范围内:2≤n≤20,优选地3≤n≤20。
[0048]
在各实施例中,最大距离d
max
包括在以下范围内:30米≤d
max
≤1000米,优选地在范围50米≤d
max
≤500米内。然而,本发明不限于待确定的特定最大距离。
[0049]
如上所定义并在权利要求书中定义的检测时间段td必须被解释为检测到反射激光的检测时间段,因此不排除在检测时间段td之后的附加检测时间窗口被进一步用于检测例如背景或执行其他测量。
[0050]
在各实施例中,对于i=1到m,连贯检测时间窗口tw[i]的持续时间基本相等。
[0051]
在各实施例中,对于每一连贯检测时间窗口tw[i],其中i=1到m,(|(pw/α)-tw[i]|/(pw/α))≤0.10,优选地(|(pw/α)-tw[i]|/(pw/α))≤0.05,更优选地(|(pw/α)-tw[i]|/(pw/α))≤0.02。
[0052]
在各实施例中,根据本发明的方法是用于使用lidar设备确定到场景中的一个或多个对象的距离的方法。
[0053]
注意,根据本发明的方法不限于所定义的步骤的特定次序,并且用字母a)到g)标识的步骤不应被解释为指示特定的时间次序。例如,步骤a)到c)可以以不同次序执行或并行执行。字母a)至g)仅用于标识由本发明方法所包括的不同步骤,并且当例如在文本中进一步详细讨论这些步骤时,这些字母用作参考。
[0054]
在各实施例中,对反射激光的检测是用激光接收设备来执行的,该激光接收设备包括具有多个像素的多阱像素检测器,并且其中每一像素包括光敏区和多个电荷存储阱w[i],其中i=1到nw,其中nw是电荷存储阱的数量,且nw≤m,并且其中该多个电荷存储阱被配置成用于存储在检测时间段td期间由光敏区检测到的电荷。在优选实施例中,nw《m且m≥4。
[0055]
在根据本发明的实施例中,在步骤c)中,如果整数α被选择成大于一,则α被选择成低于nw,优选地α=round(nw/2),其中round是将比值nw/2四舍五入为最接近整数的函数。如果nw是奇数,则除以2后的四舍五入可以设置成最低或最高整数。
[0056]
在各实施例中,α+1个连贯检测时间窗口中的每一子序列tw[k]到tw[k+α]与多个电荷存储阱w[i]中的α+1个电荷存储阱的唯一性序列相关联,使得在其中检测到反射激光的检测时间窗口是能毫无疑义地标识的,其中k=1到m-α。
[0057]
在其中α=1的各实施例中,每一对连贯检测时间窗口tw[k]到tw[k+1]与多个电荷存储阱w[i]中的唯一性的一对电荷存储阱相关联,使得在其中检测到反射激光的检测时间窗口是能毫无疑义地标识的,其中k=1到m-1。
[0058]
在其中α=1的实施例中,在步骤a)中,脉宽缩减因子n被定义成使得在电荷存储阱的数量是奇数的情况下,n≤nw!/((nw-2)!
×
2),或者在电荷存储阱的数量是偶数的情况下,n≤[nw!/((nw-2)!
×
2)]

[(nw!/((nw-1)!
×
2))-1]。
[0059]
在各实施例中,在步骤g)ii)中,对象距离的计算基于所标识的检测时间窗口和/或基于至少针对其中标识了反射激光的最早和最后一个检测时间窗口获得的曝光值。在一些实施例中,对象距离的计算基于在标识反射激光的最早检测时间窗口之前的检测时间窗口的数量的计数。根据本公开的第二方面,提供了一种用于确定到场景的一个或多个对象的距离的lidar设备。此类lidar设备包括:投影仪,其被配置成用于用由多个在空间上分开的脉冲激光束形成的光斑图案照亮该场景,其中每一脉冲激光束包括具有脉宽的脉冲序列,并且其中该多个脉冲激光束中的每一脉冲是同时发射的;包括基于cmos的多阱像素检测器的光接收设备,该多阱像素检测器被配置成用于检测反射激光的光斑,该反射激光的光斑表示由该场景的一个或多个对象反射的光斑图案;以及控制器,该控制器用于控制投影仪和光接收设备,以便与脉冲激光束的脉冲的同时发射同步地检测和积累反射激光。
[0060]
根据本发明的lidar设备的特征在于,基于cmos的多阱像素检测器被配置并可操作用于在落在脉冲激光束的脉冲周期内的检测时间段期间检测反射激光,并且其中该检测时间段被划分成多个连贯检测时间窗口tw[i],其中i=1到m,m是检测时间窗口的数量,其中m≥3,优选地m≥4,更优选地m≥5,并且其中
[0061]
pw是脉宽,td是检测时间段,且α是整数,其中α≥1,优选地α≤10,更优选地α≤5。
[0062]
对于i=1到m,连贯检测时间窗口tw[i]的持续时间基本相等。在各实施例中,对于每一检测时间窗口tw[i],其中i=1到m,(|pw-tw[i]|/pw)≤0.10,优选地(|pw-tw[i]|/pw)≤0.05,更优选地(|pw-tw[i]|/pw)≤0.02。
[0063]
基于cmos的多阱像素检测区还被配置且可操作用于对于每一所检测的反射激光的光斑,通过累积表示在检测时间窗口tw[i]中检测到的反射激光的电荷量来获得与检测时间窗口tw[i]相关联的曝光值,并且其中对于脉冲序列的所有脉冲执行该累积。
[0064]
lidar的特征还在于,控制器被配置成用于执行同步,以使得i)第一检测时间窗口tw[1]与脉冲的脉冲发射时间窗口交叠或至少部分交叠,或者另选地使得ii)第一检测时间窗口tw[1]相对于脉冲的发射时间窗口被延迟,并且该延迟等于预定义延迟时间窗口。
[0065]
优选地,预定延迟时间因子等于或小于2
×dmin
/c,d
min
是最小距离且c是光速,使得可以检测位于最小距离d
min
之处或附近的对象并确定距离。
[0066]
根据本发明的lidar设备的特征进一步在于,其包括处理装置,该处理装置被配置成用于对于每一检测到的反射激光的光斑,i)标识在连贯检测时间窗口的哪些检测时间窗口中获得代表反射激光的曝光值,ii)基于所标识的检测时间窗口和/或在这些所标识的检测时间窗口中累积的电荷,计算到场景的对象的对象距离。
[0067]
在各实施例中,基于cmos的多阱像素检测器包括多个像素,并且其中每一像素包
括光敏检测区和多个电荷存储阱w[i],其中i=1到nw且nw是电荷存储阱的数量,其中nw≤m,并且其中多个电荷存储阱被配置成用于存储在检测时间段期间由光敏区检测到的电荷。在优选实施例中,其中nw《m且m≥4。
[0068]
在包括此类基于cmos的多阱像素检测区并且其中nw《m且m≥4的各实施例中,α+1个连贯检测时间窗口tw[k]到tw[k+α]中的每一子序列与多个电荷存储阱w[i]中的α+1个电荷存储阱的唯一性序列相关联,使得在其中检测到反射激光的检测时间窗口是能毫无疑义地标识的,其中k=1到m-α。
[0069]
在m≥4且α=1的各实施例中,对于多阱像素检测器nw《m中的每一像素,检测时间窗口中的每一对连贯检测时间窗口tw[k]和tw[k+1]与唯一性的一对电荷存储阱w[m]和w[n]相关联,使得在其中检测到反射激光的检测时间窗口是能毫无疑义地标识的。
[0070]
在其中α=1的实施例中,电荷存储阱的数量nw和检测时间窗口的数量m定义如下:nw=3且m=4,或者nw=4且5≤m≤6,或者nw=5且7≤m≤11。
[0071]
在优选实施例中,lidar设备被配置成使得检测时间窗口的数量m和脉宽pw是相关的,并且该关系被定义为:pw=(tof
max

t
dl
)/((m-α)/α),其中tof
max
=2
×dmax
/c,其中tof
max
是最大飞行时间且d
max
是与最大飞行时间相关联的最大距离,c是光速,且t
dl
是预定义延迟时间窗口,其中0≤t
dl
≤(0.6
×
tof
max
),优选地0≤t
dl
≤(0.4
×
tof
max
),更优选地0≤t
dl
≤(0.2
×
tof
max
)。最大距离是lidar设备的工作范围,一般而言是lidar设备的最大工作范围。
[0072]
在各实施例中,与lidar设备的最大工作范围相对应的最大距离d
max
在以下范围内:30米≤d
max
≤1000米,优选地在范围50米≤d
max
≤500米。
[0073]
在各实施例中,投影仪被配置成用于以脉冲频率f
p
提供脉冲激光束,使得f
p
≤1/((m/α)
×
pw+t
dl
)。
[0074]
在各实施例中,由脉冲激光束产生的激光具有在800nm和1600nm之间的波长。
[0075]
在根据本发明的实施例中,α≥2。
[0076]
本发明还涉及包括根据本发明的lidar设备的交通工具。
[0077]
附图的简短描述
[0078]
将通过示例并参考附图更详细地解释本公开的这些和其他方面,附图中:
[0079]
图1示意性地解说了用于光照和检测的现有技术时序图,
[0080]
图2示意性地解说了形成脉冲链的脉冲序列的示例,
[0081]
图3示意性地解说了多个帧的重复,
[0082]
图4示意性地解说了根据本发明的实施例的脉冲发射和检测的各时序图,
[0083]
图5示意性地解说了根据本发明的实施例的脉冲发射和检测的时序图的又一些示例,
[0084]
图6示意性地解说了脉冲发射和检测的时序图,其中脉冲发射时间窗口和用于检测反射激光的第一检测时间窗口之间存在延迟,
[0085]
图7示意性地解说了与给定数量的检测时间窗口相关联的多个检测器阱的两个示例,
[0086]
图8示意性地解说了脉冲发射和检测的时序图的示例,其中α=2,
[0087]
图9示意性地解说了使用根据本发明的方法的距离测量的时间误差的减小,
[0088]
图10示意性地解说了根据本发明的实施例的lidar设备。
[0089]
这些附图的绘制既不是按比例也是不成比例的。一般而言,在附图中,相同组件由相同附图标记来表示。
[0090]
实施例的详细描述
[0091]
将根据具体实施例来描述本公开,这些实施例是本公开的解说性实施例而不应被解释为是限制性的。本领域技术人员将领会,本公开不受已经具体示出和/或描述的内容的限制,并且可以根据本公开的整体教学开发出替换或修改的实施例。所描述的附图仅是示意性的且是非限制性的。
[0092]
动词“包括”以及相应词形变化的使用并不排除所陈述的元素以外的其他元素的存在。在元素前面使用冠词“一”、“一个”或“该”并不排除存在多个这样的元素。
[0093]
此外,说明书中和权利要求中的术语第一、第二等等用于在类似的要素之间进行区分,并且不一定用于在时间上、空间上、以排名或以任何其他方式描述顺序。应当理解,如此使用的术语在适当的环境下是可互换的,并且本文中所描述的本公开的实施例能够以不同于本文中所描述或所解说的顺序操作。
[0094]
贯穿本说明书,对“一个实施例”或“一实施例”的引用意味着结合该实施例描述的特定特征、结构或特性被包括在本公开的一个或多个实施例中。因此,短语“在一个实施例中”或“在一实施例中”贯穿本说明书在各个地方的出现并不一定全部指代同一实施例,而是可以指代同一实施例。此外,在一个或多个实施例中,如通过本公开将对本领域普通技术人员显而易见的,特定的特征、结构或特性能以任何合适的方式进行组合。
[0095]
当诸如在脉冲激光束的脉宽中使用词语“脉宽”时,它必须被解释为时间脉宽,以时间单位表示,例如以纳秒表示。
[0096]
用于确定到场景的距离的方法,概述
[0097]
根据本公开的第一方面,提供了一种用于确定到场景的一个或多个对象的距离的方法。场景是例如一区域,例如安装在汽车挡风玻璃或保险杠上的lidar设备所观察的区域。取决于lidar设备的视野,场景可以覆盖大面积或较小面积。汽车应用的视野是例如30
°
x10
°
、120
°×
30
°
或任何其他合适的视野。场景可以包括例如处于距lidar设备不同距离处的各种对象,或几个对象或仅一个对象。该方法旨在对场景进行距离绘图,从而标识到诸对象的不同距离或到场景的诸部分的距离。根据本发明的方法不仅限于用于汽车应用的lidar,而且该方法还可应用于其中lidar例如安装在飞机或卫星上的其他领域。
[0098]
该方法包括用由多个在空间上分开的脉冲激光束形成的光斑图案来照亮场景的步骤。每一脉冲激光束包括具有脉宽pw和脉冲频率f
p
的脉冲序列。多个脉冲激光束中的每一者的诸脉冲是同时发射的。这种脉冲序列也被称为脉冲链(train),并且脉冲链50的示例在图2上示意性地示出。在这一解说性示例中,仅示出了5个脉冲,但实际上,脉冲链中的脉冲数通常要大得多。例如,在一些实施例中,脉冲链中的脉冲数在50到500个脉冲之间。脉冲11的脉宽pw和脉冲周期p
p
与脉冲频率f
p
成反比,如图2上所示。
[0099]
在各实施例中,如图2上所示,脉冲是块脉冲,即矩形脉冲。然而,本发明不限于脉冲的特定形状。在其他实施例中,脉冲可以具有另一形状。
[0100]
在各实施例中,由形成离散光斑图案的激光束产生的激光的波长在800nm和1600nm之间。
[0101]
在又一步骤中,检测反射激光的光斑。这些光斑表示由场景中的一个或多个对象
或场景中的一个或多个对象的一部分反射的光斑图案。该检测在脉冲激光束的每一脉冲周期p
p
内并与脉冲激光束的诸脉冲的同时发射同步地执行。
[0102]
反射激光的检测例如用光接收设备执行,光接收设备包括基于cmos的多阱像素检测器和光学元件,诸如一个或多个光学透镜和窄带滤光器。下文将更详细地讨论可被用于应用根据本发明的方法的光接收设备。
[0103]
根据本发明的方法包括确定由用于形成照亮场景的光斑图案的脉冲激光束生成的脉冲激光的最佳脉宽pw的步骤。
[0104]
根据本发明的方法,脉宽pw是相关于最大距离d
max
确定的。实际上,场景中的一个或多个对象位于最小距离d
min
和最大距离d
max
之间的距离范围中的距离ds处。一般而言,0≤d
min
≤0.6
×dmax
,优选地0≤d
min
≤0.4
×dmax
,更优选地0≤d
min
≤0.2
×dmax
。换言之,对象可以位于不同的距离处,并且该方法允许确定到一个或多个对象的距离,无论该一个或多个对象实际上位于在d
min
和d
max
之间的任何距离处。
[0105]
该最大距离可解释为所需的最大工作范围,即最大距离对应于仍应是能毫无疑义地测量的到场景的最大距离。例如,如果将d
max
定义为100米,则意味着根据本发明的方法,场景中的位于100米或以下的对象应该是可检测的,并且距离是可测量的,而到场景中位于最大距离之外的对象的距离不需要能够被确定。对于根据本发明的方法,该最大距离通常等于或大于50米,并且通常等于或小于500米。但取决于应用,待测的最大距离也可更大,例如1000米或甚至更大。根据本发明的方法不受具体最大距离的限制。
[0106]
根据本发明的方法,脉宽pw被设置为小于与最大距离d
max
相关联的最大飞行时间tof
max
。当激光在已经被发射并被对象反射回来后被检测到时,它所经过的距离是到对象的距离的两倍。因此,tof
max
等于(2
×dmax
/c),其中c是光速。脉宽pw被设置成等于[tof
max

t
dl
]/n,且n是预定义脉宽缩减因子,且t
dl
是预定延迟时间窗口。
[0107]
预定延迟时间因子通常等于或小于2
×dmin
/c,使得可以检测位于最小距离d
min
之处或附近的对象并确定距离。
[0108]
在图4中,示出了其中预定延迟时间窗口为零的各实施例。如图4示意性地解说的,在这一情形中下,脉宽缩减因子n对应于最大可测量飞行时间tof
max
与所发射激光脉冲11的脉宽pw之间的比率。
[0109]
根据本发明的方法,对于每一检测到的反射激光的光斑,获得作为时间的函数的曝光值。
[0110]
根据本发明的方法包括确定落在脉冲激光束的脉冲周期p
p
内与用于检测反射激光的检测时间段td的步骤,并且还包括将检测时间段td划分成m个连贯检测时间窗口tw[i],其中i=1到m,m是检测时间窗口的数量,并且其中m=α
×
(n+1)。检测时间段td表示为:
[0111][0112]
其中α是整数,其中α≥1。因子α表示每脉宽的时间窗口的数量。下面将进一步讨论其中α=1和α》1的两个实施例。
[0113]
在其中α=1的各实施例中,预定义脉宽缩减因子n通常等于或大于2,优选地等于或大于3,且更优选地等于或大于4。
[0114]
在各实施例中,α和n被选择成使得α
×
n≥2,优选地α
×
n≥3,更优选地α
×
n≥4。
[0115]
每一检测窗口tw[i]必须被解释为一时间段。在各实施例中,对于i=1到m,连贯检测时间窗口tw[i]的持续时间基本相等。
[0116]
在图4中,示出了对于n从2到5以及α=1的值而言的检测时间窗口的数量。图4上示出的各实施例是其中预定延迟时间窗口为零的实施例。
[0117]
在图6中,示出了实施例的示例,其中预定延迟时间窗口t
dl
不为零。在各实施例中,0≤t
dl
≤0.6
×
tof
max
,优选地0≤t
dl
≤0.4
×
tof
max
,更优选地0≤t
dl
≤0.2
×
tof
max
。本领域技术人员将取决于所使用的检测装备以及取决于需要使用同一装备测量的对象距离的范围来选择应用或不应用延迟时间窗t
dl
。例如,如果目的是测量在最小距离d
min
等于d
max
的10%和最大距离d
max
之间的范围内的对象距离,则可以设置等于tof
max
的10%的时间延迟,以减少用于检测反射激光的检测器的像素阱的数量,这将在下文中进一步讨论。另一方面,如果对象能定位在0米(即d
min
=0)和最大距离d
max
之间的任何距离处并需要在此被检测,则本领域技术人员可以选择不实现延迟时间窗口。
[0118]
如图4、图6和图8上所示,每一检测时间窗口tw[i]几乎等于比率pw/α,其中pw是所发射的激光脉冲11且α是阿尔法因子,对于图4和图5的示例而言α等于1且对于图8的示例而言α等于2。在一些实施例中,每一检测时间窗口tw[i]不一定完全等于比率pw/α,且只要检测时间段td等于m乘以比值pw/α,就会出现一定量的偏差。一般而言,对于每一时间窗口tw[i],(|(pw/α)-tw[i]|/(pw/α))≤0.10,优选地(|(pw/α)-tw[i]|/(pw/α))≤0.05,更优选地(|(pw/α)-tw[i]|/(pw/α))≤0.02。本领域技术人员将取决于所使用的装备和用于确定到对象的距离所需的最终准确度来选择最佳值。在一些实施例中,(|(pw/α-tw[i]|/(pw/α))≤0.01。
[0119]
对于每一检测到的反射激光的光斑,该方法包括通过累积在检测时间窗口tw[i]中检测到的反射激光的量来获得与检测时间窗口tw[i]相关联的曝光值qi的步骤。这一积累是在产生光斑的脉冲激光束的脉冲序列的所有脉冲上进行的。事实上,取决于到场景中的反射激光的对象的距离,会出现不同的飞行时间tof,并且因此取决于要测量的距离,可以在不同时间窗口中检测到激光。
[0120]
如果α=1,并且由于脉宽pw几乎与检测时间窗口的宽度相同,反射脉冲12将在两个连贯检测时间窗口或单个检测时间窗口中被检测到。在图4中,示出了一些任意tof示例,其解说了反射脉冲通常与两个连贯时间窗口交叠。
[0121]
对于某些对象距离,并且如果α=1,则反射激光也可以在单个时间窗口内被检测到。例如,如果第一检测时间窗与脉冲11的发射一致,即如果没有应用延迟时间窗口t
dl
,并且如果到待测量对象的距离对应于等于或多个脉宽pw的飞行时间,则反射激光将在单个检测时间窗口中被检测到。类似地,如果应用了延迟时间窗口t
dl
,并且如果到待测量对象的距离对应于等于延迟时间窗口t
dl
加上多个脉宽pw的飞行时间,则反射激光在单个检测时间窗口中被检测到。
[0122]
更一般而言,对于α等于或大于1的任何值,反射激光在α个或α+1个连贯检测时间窗口中被检测到。
[0123]
根据本发明的方法,确定图2上解说的脉冲频率f
p
,使得f
p
≤1/((n+1)
×
pw+t
dl
)。脉冲频率的上限被定义成使得从位于最大距离d
max
处的对象反射的反射激光在发射后续脉
冲之前被检测到。一般而言,脉冲频率低于上限1/(((n+1)
×
pw)+t
dl
),以将平均激光功率保持在例如出于眼睛安全原因或为了避免混叠问题或为了执行2d测量而定义的限度内。在各实施例中,如图2上所解说的,脉冲链50的脉冲频率f
p
通常在khz范围内,例如在10khz和500khz之间。
[0124]
在一些实施例中,频率f
p
被确定为使得1/((2
×n×
pw)+t
dl
)≤f
p
≤1/(((n+1)
×
pw)+t
dl
)。
[0125]
根据本发明的方法提供了针对每一检测到的反射激光的光斑来确定对象距离的又一步骤。通过首先标识在其中检测到反射激光的相应时间窗口来确定对象距离。在第二步骤中,基于对所标识的相应检测时间窗口的标识和/或在所标识的检测时间窗口中获得的曝光值,来计算到场景的对象的对象距离ds,其中ds≤d
max
,如下文更详细地讨论的。
[0126]
对于现有技术的时序图,例如如图1上所示,仅使用两个检测时间窗口,并且因此可以基于在两个时间窗口中检测到的激光量来确定对象距离。另一方面,根据本发明的时序图,例如如图4上所示,没有提前知悉将在什么检测时间窗口内检测到反射激光。因此,当在图案照明序列期间已经执行了对反射激光的检测之后,根据本发明的方法包括标识在其中已检测到激光的检测时间窗口的步骤。针对检测到的每一反射激光的光斑执行该步骤。如上所述,利用根据本发明的方法,反射激光在α个或者在α+1个连贯检测时间窗口中被检测到。
[0127]
在读出在检测时间窗口中获得的数据后,可以使用各种软件算法来执行对已检测到反射激光的检测时间窗口的标识。事实上,反射激光在α个或者α+1个连贯检测时间窗口被检测到,并且因此通过比较每一检测时间窗口中的计数数量,具有最高计数数量的α或α+1个检测时间窗口可以与具有仅检测到背景计数的检测时间窗口区分开。另外,该算法还可以将α+1个检测时间窗口中的预期强度曲线纳入考虑,其中在检测到反射激光的最早和最后一个检测时间窗口中预期强度较低,而在最早和最后一个检测时间窗口之间预期强度较高且基本相同。对于在α个检测时间窗口中检测到反射激光的情况,在α个连贯检测时间窗口中预期基本上相同的强度。注意,如果例如α=1,则反射激光在最多两个检测时间窗口中被检测到,并且因此始终至少一个检测时间窗口可供用于背景检测。
[0128]
在各实施例中,例如使用块函数的拟合算法可被应用来确定具有最高计数数量的检测时间窗口的位置。
[0129]
在标识出检测到反射激光的检测时间窗口之后,可以基于所标识的检测时间窗口来计算到场景的对象的距离。
[0130]
在各实施例中,如更详细地进一步讨论的,对象距离的计算基于所标识的检测时间窗口和/或基于至少针对其中标识了反射激光的最早和最后一个检测时间窗口获得的曝光值。在其他实施例中,对象距离的计算基于在标识反射激光的最早检测时间窗口之前的检测时间窗口的数量的计数。对于α=1,脉宽缩减因子n的上限主要由可供用于检测反射激光的检测器技术且尤其是检测器响应时间(例如cmos相机的速度)来确定。例如,对于d
max
=100米并且使用等于10的脉宽缩减因子,脉宽pw等于约66纳秒。这样的脉宽仍然与例如当前基于cmos的像素传感器技术兼容。
[0131]
对于α》1,具有缩减因子n的脉冲是相关于缩减因子α来确定的。事实上,如上所讨论的,检测时间窗口取决于n和α。如下文将讨论的,如果使用的检测器技术包括需要打开和
关闭的电荷存储阱,则检测时间窗口应保持在与cmos技术兼容的时间窗口内。在各实施例中,n和α的积等于或小于1000,优选地等于或小于500,更优选地等于或小于100。本领域技术人员将基于可用的检测技术并且还基于最大距离d
max
来选择α和n的值。d
max
越大,积α
×
n要取得越大,以获得具有给定经缩减宽度的检测时间窗口。
[0132]
用于确定到对象的距离的公式
[0133]
由于对反射激光的检测是与脉冲激光束的诸脉冲的同时发射同步地执行的,基于对检测到激光的检测时间窗口的标识,可以确定激光的飞行时间并且因此确定到对象的距离。
[0134]
本领域技术人员将定义用于基于已检测到激光的所标识的检测时间窗口来计算到对象的距离的公式。可以使用各种选项来计算对象距离。在这里下面将概述用于计算距离的公式的示例。如将要讨论的,在一些实施例中,这些公式将一些所标识的检测时间窗口的曝光值纳入考虑,尤其是其中标识了反射激光的最早和最后检测时间窗口的曝光值。在其他实施例中,仅基于对已检测到反射激光的检测时间窗口的标识来确定对象距离。在这些实施例中,为了确定对象距离,对检测到反射激光的最早检测时间窗口之前的检测时间窗口的数量进行计数。
[0135]
首先讨论在α=1(即脉宽等于检测时间段)时用于计算距离ds的公式。如果在两个连贯时间窗口i和i+1中检测到反射激光,则可以使用以下公式计算距离ds:
[0136]ds
=(t
dl
×
c/2)+[(i-1)
×c×
pw+c
×
pw
×qi+1
/(qi+q
i+1
)]/2,
[0137]
其中qi和q
i+1
是分别在时间窗口i和时间窗口i+1中获得的曝光值。通常,曝光值会针对例如由到达检测器的环境光或其他背景噪声源产生的背景事件进行校正。
[0138]
或者,如果在两个连贯时间窗口i和i+1中检测到反射激光,则也可以使用以下替换公式计算距离ds:ds=(t
dl
×
c/2)+[i
×c×
pw-c
×
pw
×
qi/(qi+q
i+1
)]/2。如果时间窗口i中的反射激光的量大于时间窗口i+1中的反射激光的量,则优选地使用该替换公式。
[0139]
另一方面,如果仅在一个时间窗口i中检测到反射激光,其中i≥2,则可以使用以下公式计算到场景的距离:ds=(t
dl
×
c/2)+[(i-1)
×c×
pw]/2。
[0140]
图8上示意性地示出了检测时序图的示例,其中α=2。在该示例中,脉宽缩减因子n等于15,即所发射的脉冲11的脉宽pw等于最大飞行时间除以15。在α=2时,检测时间窗口的总数m等于32=α
×
(n+1)。对于如图8中所示的检测时序图,其中α》1,还可以推导各公式以确定距离ds。
[0141]
如果α》1,则在α+1个检测时间窗口或α个检测时间窗口中检测到反射激光。例如,如果α=2,即tw[i]=pw/2,则在最多达三个连贯时间窗口中检测到反射激光。一般而言,对于α=2,反射脉冲部分地在时间窗口tw[i]中观察到,完全在时间窗口tw[i+1]中检测到,且部分地在时间窗口tw[i+2]中检测到。在某些情况下,取决于对象在场景中的位置,还可以在两个连贯时间窗口tw[i]和tw[i+1]中检测到反射激光。对于图8中所示的示例,在时间窗口14到16中检测到反射激光脉冲12。
[0142]
如果α≥1且如果在α+1个连贯检测时间窗口中检测到反射激光,则根据本发明的方法确定距离的更一般公式可以如下表示为:
[0143]ds
=t
dl
×
c/2+z
×
(tof
max-t
dl
)
×
c/2,
[0144]
其中z=(r+z)/(α
×
n)且z=qb/(qa+qb),
[0145]
在该公式中,ds是对象距离,qa是针对连贯检测时间窗口中的在其中标识反射激光的最早检测时间窗口tw[a]获得的曝光值,qb是针对连贯检测时间窗口中的在其中标识反射激光的最后检测时间窗口tw[b]获得的曝光值,且r是对最早时间窗口tw[a]之前的检测时间窗口的数量进行计数的整数。在图8上所示的示例中,反射激光在时间窗口14到16中被检测,且最早的检测时间窗口tw[a]是其中检测到反射激光的检测时间窗口14,且检测时间窗口16是其中检测到反射激光的最后检测时间窗口tw[b]。本示例中的整数r等于13,对应于最早检测时间窗口14之前的检测时间窗口的数量。这些曝光值qa和qb必须被解释为已针对背景计数(例如由背景光产生的背景计数)进行校正的曝光值。因此,如果在α+1个连贯检测时间窗口检测到反射激光,并且当使用上述公式时,通过将在所标识的检测时间窗口中获得的曝光值纳入考虑来确定距离。
[0146]
如果α≥1且在α个连贯检测时间窗口中而不是在α+1个连贯检测时间窗口检测到反射激光,则可以使用以下表达式计算对象距离ds:
[0147]ds
=t
dl
×
c/2+(r/(α
×
n))
×
(tof
max-t
dl
)
×
c/2,
[0148]
其中r是对在其中标识反射激光的最早检测时间窗口之前的检测时间窗口的数量进行计数的整数。因此,如果在α个连贯检测时间窗口中检测到反射激光,则在不使用在所标识的检测时间窗口中获得的曝光值的情况下确定距离。
[0149]
在一些实施例中,对于在α或者α+1个连贯检测时间窗口中检测到反射激光的两种情形,对象距离ds可在不考虑在所标识的检测时间窗口中获得的曝光值的情况下确定。实际上,如果检测时间窗口很短,例如通过为α和/或脉宽因子n选择很大的值,则在距离计算中可以忽略最早检测时间窗口,其中仅在该检测时间窗口的一部分期间检测到反射激光。在这些实施例中,可以使用与上述相同的表达式来计算对象距离ds,即:ds=t
dl
×
c/2+(r/(α
×
n)x(tof
max-t
dl
)x c/2,其中r是对在其中标识反射激光的最早检测时间窗口之前的检测时间窗口的数量进行计数的整数。例如,这种用于计算对象距离的办法可以应用于以下情况:α和/或n很大或者待确定距离很大,例如百米或更大。
[0150]
激光的脉冲频率
[0151]
根据本发明的方法可达到的最高脉冲频率f
p
等于1/((n+1)
×
pw+t
dl
)。然而,脉冲频率可以更低并且等于例如1/((n+1)
×
pw+t
dl
+(q
×
pw)),其中q≥1且其中q
×
pw形成延迟时隙,该延迟时隙被定义成使得对从场景中的位于大于d
max
的距离处的对象的反射产生的虚假反射的检测被消除以避免混叠。在其他实施例中,还可以出于眼睛安全的原因实现延迟时隙,以便例如降低所发射激光的平均功率。当使用如图1所示的仅具有两个检测时间窗口的现有技术时序图时,即其中n=1的图,通常在第一和第二时间窗口之后提供附加时隙以用于背景测量。如上所讨论的,使用根据本发明的时序图,用于背景测量的附加时隙不是强制性的,因为背景估计总是可以从检测时间窗口tw[i]中的在其中没有检测到反射激光并且只检测到背景的一个检测时间窗口获得。
[0152]
然而,在根据本发明的一些实施例中,上述附加延迟时隙q
×
pw可以被用于或可以部分地被用于执行附加背景测量。因此,检测时间段td必须被解释为检测反射激光的检测时间段,它不排除检测时间段td之后的附加检测时间窗口被用于检测例如背景,或者附加时间窗口被用于任何其他合适目的。
[0153]
在图5中,示出了根据本发明的时间图的示例,其中在检测时段td之后提供了附加
延迟时隙。检测时间段之后的延迟时隙在图5上用附图标记21来标识。在图5上所示的这些示例中,没有实现如上所讨论的预定延迟时间窗口t
dl
,并且如图所解说地执行同步,使得所发射的脉冲11与第一检测时间窗口交叠。在其他实施例中,可以另外地实现预定延迟时间窗t
dl
,以相对于脉冲的发射时间窗口来延迟第一检测时间窗口,如上所述。如上所讨论的,该附加延迟时隙21可被用于多个目的,减少虚假反射、出于眼睛调节目的来降低平均功率、或甚至用于执行附加背景测量。
[0154]
在各实施例中,脉冲频率f
p
等于或小于1/((n+1)
×
pw+t
dl
+im),其中im是提供用于使用环境光或照亮场景的光源对场景进行成像以获得场景的2d摄影图像的时隙,使得用tof方法确定的距离可以与该摄影图像相组合以形成场景的3d图像。如上所述,不强制提供附加时隙来测量用于校正曝光值所需的背景,因为总是存在在其中没有激光被反射的可用时间窗口。然而,在一些实施例中,提供用于测量背景以校正和改善2d图像的附加时隙。
[0155]
脉冲检测和脉冲发射之间的同步如上所讨论的,反射激光的光斑是与脉冲激光束的诸脉冲的同时发射同步地检测的。
[0156]
在其中延迟时间窗t
dl
如上所规定为零的实施例中,执行激光束脉冲的发射与反射激光的检测的同步,使得第一检测时间窗口tw[1]与脉冲的发射时间窗口交叠或至少部分交叠。在图4和图5上所示的其中α=1的实施例中,第一时间窗口tw[1]被解说为与激光脉冲11的发射一致,即脉冲的脉宽与第一时间窗口完全交叠。
[0157]
在其他实施例中,第一检测时间窗口tw[1]与脉冲的发射时间窗口部分地交叠。这可例如在第一检测时间窗口tw[1]相对于脉冲的发射被延迟时发生。如果该延迟小于第一检测时间窗口,且因为在α=1的情况下检测时间窗口具有与脉冲相同或几乎相同的宽度,则在第一时间窗口和脉冲的发射之间始终存在部分交叠。在各实施例中,如果脉宽和检测时间窗口不完全相同,则第一检测时间窗口也可以与脉冲的发射部分交叠,即不完全交叠。更具体而言,如果α》1,则脉宽pw大于检测时间窗口,并且因此在同步时,激光脉冲的发射和第一检测时间窗口之间可发生部分交叠。第一检测时间窗口tw[1]与脉冲的发射时间窗口交叠或至少部分交叠。
[0158]
在各实施例中,通过在第一检测时间窗口tw[1]期间在激光脉冲的发射之前开始检测,第一检测时间窗口与脉冲的发射时间窗口部分交叠,即激光脉冲的发射相对于第一检测时间窗口被延迟。因此,即使t
dl
如上所规定为零,仍有一些余量用于同步激光脉冲的发射和第一检测时间窗口,并且例如通过同步引入正或负延迟,使得第一检测时间窗口和激光脉冲的发射之间仅存在部分交叠。
[0159]
更一般而言,如果预定义延迟时间窗口t
dl
如上所定义地不为零,例如在图6上所示,则执行所述同步以使得第一检测时间窗tw[1]相对于脉冲的发射时间窗口被延迟且延迟等于延迟时间窗口t
dl

[0160]
距离准确度,时间精度
[0161]
如上所述,通过拍摄多个帧并根据在每一单个帧测量期间获得的各对象距离ds来计算平均对象距离,计算得到的平均对象距离值上的误差在与单个帧测量的时间误差相比时以因子减小,其中nf是帧的数量。使用根据本发明的方法,目的是通过单个帧测量获得高精度,以便减少确定距离所需的帧数,且甚至使用单个帧执行距离测量。帧可以以
帧速率ff重复,其通常远低于脉冲频率f
p
。在图3中,示意性地解说了帧60的重复,并示出了帧率ff。在该示例中,示出了三个帧的重复,在实践中,帧重复的数量通常要大得多,因为一般而言,例如对于安装在汽车或无人机上的lidar设备,距离是以连续的方式测量的。如图3示意性地解说的,在每一脉冲链50之后,需要处理时间65来读出曝光值并处理所采集的数据。可达到的帧率ff通常在hz范围内,在各实施例中,帧率例如在5hz和50hz之间。
[0162]
例如在图4、图6和图8上所示的根据本发明的时序图的主要优点是,与图1上所示的现有技术时序图相比,单个帧中的距离测量的时间精度大大提高。这是以下事实的结果:如果检测时间窗口变得更小,则时间精度会提高。当将使用当前方法获得的精度与使用图1上所示的现有技术时序图(其中使用了两个检测窗口)获得的精度进行比较时,获得因子等于的改进。因此,如果α=1,改进因子与分割因子n成正比。
[0163]
利用现有技术时序图,需要拍摄大量帧以获得在可接受值内的计算得到的平均距离上的标准偏差单个帧测量的时间误差在这里表示为标准偏差σ。单个帧测量是基于脉冲序列50的测量,如图2上所示。随着使用根据本发明的方法的单个帧测量的精度的提高,可以减少帧的数量nf,或者甚至可以在单个帧内执行测量。以此方式,可以快得多地获得到场景的对象的距离。下文将进一步讨论精度的提高以及因此帧数的减少。
[0164]
具有位于给定检测时间窗口tw[i]中的相应飞行时间的对象距离的时间误差,用标准偏差σ表示,可以写成如下:
[0165][0166]
其中σ
ζ
是独立于脉宽pw且仅取决于一个帧内的总信号强度s、在一个帧期间的测得背景b、像素噪声σ
px
和等于qb/(qa+qb)的计数比因子ζ的精度因子,其中qa和qb是在其中检测到反射激光的最早和最后检测时间窗口中检测到的曝光值。
[0167]
更准确而言,对于α≥1,上述公式中提到的σ的精度因子σ
ζ
与以下因素成正比:
[0168]
其中z=qb/(qa+qb),其中如上所讨论的,qa和qb分别是在其中检测到反射激光的最早和最后检测时间窗口的曝光值。
[0169]
时间误差σ的示例在图9上针对划分因子n和因子α的不同值示出。图9上的σ值是相对值,以解说对于所示的各种示例而言σ的相对差异。各示例在图9上以字母a到h标识,并且图9上的图例指示所选择的n和α的对应值。sigma的值在此表示为比率d/d
max
=tof/tofmax的函数,从0到1,即覆盖用于确定距离的整个距离范围。在该示例中,延迟时间窗口t
dl
假定为零。现有技术脉冲时序图,其中α=1且n=1,如图1上所示,被用作基准时间误差σ。使用现有技术脉冲时序图,值σ在tof/tof
max
=0.5时具有最大值。如用示例b到h所解说的,其中α
×
n≥2,对于这些示例中的每一示例,作为tof/tof
max
的函数的σ值表示重复的弧形图案。每一弧代表在一检测时间窗口内的时间误差σ,该检测时间窗口与1/(α
×
n)成正比,如上所讨论的。在图9中,对于α=1与n=2、3、5和10相组合,以及对于α=2与n=1.5、n=6和n=15相组合,示出了根据本发明的示例。如图9上所解说的,与现有技术时序图相比,根据本发明的新
时序图,σ的值以因子减少。
[0170]
例如,对于α=1,通过使用脉宽缩减因子n=4,时间精度提高了四倍,并且为获得计算出的平均对象距离的相同标准偏差而拍摄的帧的数量可以减少两倍。然而,这是在脉冲链期间发射的光子总数保持不变并且背景和噪声贡献保持不变的条件下进行的。例如,如果脉冲链中的脉冲数保持不变,则如果脉宽缩减因子是4且如果目的是将帧数减少2倍,则脉冲内发射的光子数需要增加4倍。
[0171]
如果α》1,则在与α=1相比时,时间检测窗口被进一步减少,并且因此单个帧测量的精度被进一步提高。但是,精度改进因子不与α乘以n成线性关系,而是改进因子与成正比,如上所述。这一附加因子得自如下事实:单个检测时间窗口中的最大可检测激光与1/α成正比。例如,如果α=2且n等于15,则与图1上所示的现有技术相比,精度的提高是倍。如果α=1且n等于15,则精度的提高是15倍。
[0172]
注意,图9上所示的曲线是假设不存在背景计数的理论曲线。在实践中,要针对背景计数来校正曝光值,并且因此取决于背景的量,在考虑背景校正后,上述定义的改进因子可能略有不同。
[0173]
多阱像素检测器的使用以及n和α的确定
[0174]
在各实施例中,为了执行根据本发明的方法,反射激光检测是用激光接收设备来执行的,该激光接收设备包括具有多个像素的多阱像素检测器,并且其中每一像素包括光敏区和多个电荷存储阱w[i],其中i=1到nw,其中nw是电荷存储阱的数量,其中nw≤m。 这些多个电荷存储阱被配置成用于存储在检测时间段td期间由光敏区域检测到的电荷。
[0175]
在各实施例中,对于每一像素,一电荷存储阱与每一检测时间窗口tw[i]相关联,并且因此在这些实施例中,电荷存储阱的数量nw等于检测时间窗口的数量m。
[0176]
在优选实施例中,电荷存储阱的数量nw低于检测时间窗口的数量m。因此,在这些实施例中,一个或多个电荷存储阱存储一个以上检测时间窗口的电荷。这允许在不增加电荷存储阱数量的情况下使用更高的划分因子n。确立了一个关系,该关系定义哪个电荷存储阱存储在什么检测时间窗口期间检测到的电荷。为了确立检测时间窗口和电荷存储阱之间的关系,需要定义电荷存储阱的唯一性序列,如下所述。
[0177]
在包括多阱像素检测器的实施例中,α+1个连贯时间窗口tw[k]到tw[k]+α]的每一子序列与α+1个电荷存储阱的唯一性序列w[i]相关联,其中k=1到m-α。因此,定义了α+1个电荷存储阱的m-α个唯一性序列。电荷存储阱的这些唯一性序列被配置成使得当在α或α+1个电荷存储阱中检测到反射激光时,可标识在其中检测到反射激光的相应检测时间窗口,并且对象距离是使用以上讨论的公式确定的。电荷存储阱的唯一性序列以及与检测时间窗口的关联定义了什么电荷存储阱被用于存储在什么检测时间段期间的电荷。
[0178]
α+1个电荷存储阱的唯一性序列必须被解释为电荷存储阱序列,该电荷存储阱序列被配置成允许毫无疑义地标识在什么最早检测时间窗口和什么最后检测时间窗口中检测到反射激光。
[0179]
例如,如果α=1,则电荷存储阱序列α+1是一对,并且如果例如nw等于3,则如果对
1-2被定义为与一对检测时间窗口相关联的电荷存储阱的唯一性对,那么对2-1不是另一唯一性对。事实上,如果将使用对1-2和2-1,则无法确定在其中检测到激光的相应检测时间窗口,因为无法区分开这两对。另一方面,对3-1和/或对3-2可被用作又一唯一性对。
[0180]
类似地,如果α=2,则例如序列1-2-3可被定义为与三个连贯检测时间窗口的序列相关联的电荷存储阱的唯一性序列,其中阱1和阱3可检测例如分别对应于最早和最后检测时间窗口的激光。然而,序列3-2-1不能用作又一唯一性序列,因为如果将在阱1、2和3中检测到反射激光,则不能确定在什么检测时间段内有效地检测到反射激光。例如,可以使用第四存储阱4配置又一唯一性序列,其中组合4-2-3可以形成三个电荷存储阱的又一唯一性序列。
[0181]
如果例如α=4且nw等于5,序列1-2-3-4-5可被选择成与用于检测反射激光的5个连贯检测时间窗口的序列相关联的5个阱的唯一性序列。在该序列中,在阱1和阱5中可以检测到对应于例如最早和最后检测时间段的反射激光。如果这一序列1-2-3-4-5被用作唯一性序列,则序列5-2-3-4-1、序列1-3-4-2-5、序列5-3-4-2-1、序列1-4-2-3-5、序列5-4-2-3-1、序列1-4-3-2-5和序列5-4-3-2-1都可被用作又一唯一性序列。事实上,如果阱序列的任一此类又一组合将与5个连贯检测时间窗口的又一序列相关联,则与序列1-2-3-4-5相对应的检测时间窗口期间检测到的反射激光是无法区分开的。
[0182]
在图8中,示出了一示例,其中n=15且α=2,且因此检测时间窗口的数量m等于32,并且其中电荷存储阱的数量nw等于5,即《m。在该示例中,如图8上所示,电荷存储阱tw[1]例如可用于在检测时间窗口1、4、7、10、13、16和31期间对电荷进行积分。
[0183]
对于其中使用多阱像素检测器的这些实施例,其中每一像素包括给定数量nw个电荷存储阱,划分因子n和定义检测时间窗口的数量m=α
×
(n+1)的阿尔法因子被选择成使得检测时间窗口的数量m不大于对应于可以由给定数量nw个电荷存储阱形成的、α+1个连贯时间窗口的唯一性序列的最大数量的最大值m
max

[0184]
对于图9上所示的示例,图例中示出了所使用的电荷存储阱的数量nw。对于示例b到h,所使用的电荷存储阱的数量在3到5之间变化。
[0185]
各对电荷存储阱,α=1
[0186]
如上所述,因子α可以等于或高于1。首先,将讨论具有α=1的多个实施例。
[0187]
在n≥3且对于多阱像素检测器nw《m中的每一像素的各实施例中,检测时间窗口中的每一对连贯检测时间窗口tw[k]和tw[k+1]与唯一性的一对电荷存储阱w[m]和w[n]相关联。例如,在图7中,示出了n=10的示例,并且因此检测时间窗口的数量等于11。在图7上,示出了两个示例,其中一个nw=11,而另一个nw=5。对于其中nw=5,因此nw<m的实施例,每一对连贯时间窗口与唯一性的一对电荷存储阱相关联。例如,如图7上所解说的,连贯检测时间窗口tw1和tw2例如与唯一性阱对w1和w2相关联,连贯检测时间窗口tw7和tw8例如与唯一性阱对w1和w4相关联,连贯检测时间窗口tw10和tw11例如与唯一性阱对w5和w1相关联。如上所述,对于唯一性的一对电荷存储阱,将理解,例如阱对w1-w2与阱对w2-w1是同一对。
[0188]
在各实施例中,其中n≥3并且其中对于多阱像素检测器nw《m的每一像素,多个电荷存储阱中的至少一个电荷存储阱被配置成用于累积并存储在检测时间窗口的两个不同时间窗口tw[k]和tw[m]期间检测到的电荷,并且其中两个不同时间窗口是两个非连贯时间窗口,使得|m-k|》1。在图7中,对于nw=5的情况,例如电荷存储阱w5正在累积并存储来自检
测时间窗口tw5和tw10这两者的电荷,例如电荷存储阱w3在检测时间窗口tw3和tw6期间累积并存储电荷,例如电荷存储阱w1在检测时间窗口tw1、tw7和tw11期间累积并存储电荷。
[0189]
在根据本发明的方法的其中α=1的实施例中,在应用确定脉宽pw并应用脉宽缩减因子n的步骤中,脉宽缩减因子n被定义成使得n≤n
max
,其中在电荷存储阱的数量是奇数的情况下,n
max
=nw!/((nw-2)!
×
2),或者在电荷存储阱的数量是偶数的情况下,n
max
=[nw!/((nw-2)!
×
2)]

[(nw!/((nw-1)!
×
2))-1]。这一最大数量n
max
定义了给定数量的电荷存储阱可以形成的唯一性对的最大数量。当对n采用如上述公式中定义的这一上限时,对于给定数量的阱,时间检测窗口的最大数量m
max
等于n
max
+1。
[0190]
α+1个电荷存储阱的唯一性序列,α》1
[0191]
将进一步讨论其中α大于1的实施例。使用大于1的因子α的优点是对于给定数量nw的电荷存储阱,可以形成的α+1个电荷存储阱的唯一性序列的数量更高,并且因此对于相同数量的电荷存储阱,可以使用较短的检测时间窗口,从而得到经改进的精度。
[0192]
例如,如果α=2,即tw[i]=pw/2,则在最多达三个连贯时间窗口中检测到反射激光。因此,需要至少三个电荷存储阱来检测反射的激光,并需要第四电荷存储阱用于测量背景。一般而言,对于α=2,反射脉冲部分地在时间窗口tw[i]中观察到,完全在时间窗口tw[i+1]中检测到,且部分地在时间窗口tw[i+2]中检测到。在某些情况下,取决于对象在场景中的位置,还可以在两个连贯时间窗口tw[i]和tw[i+1]中检测到反射激光。
[0193]
对于α=2且对于给定数量nw的电荷存储阱,可以计算出3个电荷存储阱的唯一性序列的最大数量c
max
。对于从4到10的nw,c
max
分别等于12、30、60、105、168、252和360。
[0194]
更一般而言,对于任何数α,可由nw个电荷存储阱形成的、α+1个电荷存储阱的唯一性序列的最大数量c
max
可被计算出。
[0195]
本领域已知用于计算唯一性序列的最大数量c
max
和用于确立电荷存储阱的唯一性序列的数学算法。例如,当使用例如5个电荷存储阱时,为了计算3个唯一性电荷存储阱的序列,如图8上所示,使用迭代程序,将电荷阱序列必须允许在最早和最后检测时间窗口中毫无疑义地检测到的激光的边界条件纳入考虑。在图8上所示的该示例中,其中α=2且nw=5时,可以形成的3个电荷存储阱的唯一性序列的最大数量c
max
等于30。因此,在该示例中,可使用的最大脉冲划分因子n
max
等于15。
[0196]
唯一性序列的最大数量c
max
定义了最大脉宽缩减因子n
max
,其中α
×nmax
=c
max
。因此,当使用每像素具有给定数量nw个电荷存储阱的基于cmos的多阱像素检测器时,脉宽缩减因子n被选择成使得n等于或小于n
max

[0197]
对于选择相对小的α和相对小数量nw的电荷存储阱,唯一性序列的最大数量c
max
可变大。例如,对于α=2且nw=4,c
max
=12,对于α=2且nw=6,c
max
=60。通过取α=3或更大并取nw=7或更大,可以形成数百个唯一性组合,且c
max
变得大于100。如果取α=4或更大并取nw=9或更大,可以形成数千个唯一性组合,且c
max
变得大于1000。
[0198]
用于确定到场景的距离的lidar设备
[0199]
根据本发明的第二方面,提供了一种用于确定到场景的一个或多个对象的距离的lidar设备。图10上示意性地解说了lidar设备1的实施例。
[0200]
此类lidar设备1包括配置成用于用照明图案150来照亮场景99的投影仪100。该图案由多个在空间上分开的脉冲激光束形成,其中每一脉冲激光束包括脉冲序列。多个脉冲
激光束的诸脉冲是同时发射的。如上所讨论的,脉冲序列,也称为脉冲链在图2上示意性地解说。
[0201]
优选地,激光由半导体激光器产生。在各实施例中,激光由包含垂直腔面发射(vcsel)激光器阵列的半导体光源产生,以产生所需图案。为了使该系统即使在长距离和高环境光水平(例如,在日光下)下也能最优地工作,本发明的实施例中使用的vcsel优选地被布置成每单位面积每光斑发射最大光学功率。更优选地,激光器应具有最小的波长差异;用单模激光器可以实现特别低的波长差异。因此,可以以必要的空间和时间准确度来可再现地产生基本相同的脉冲。
[0202]
除了半导体激光器之外,投影仪还包括投影透镜,该投影透镜被配置成在所定义的视野内用光图案照亮场景。
[0203]
在各实施例中,由激光器产生的激光的波长在800nm和1600nm之间。
[0204]
该lidar设备还包括光接收设备300,该光接收设备300包括配置用于检测反射激光的光斑的基于cmos的多阱像素检测器。多阱像素检测器包括多个像素,其中每一像素包括用于检测激光的光敏区。如上所讨论的,这些反射激光的光斑代表由场景中的一个或多个对象反射的在空间上分开的光束。例如,在wo2017068199中公开了基于cmos的多阱像素检测器。反射激光在图10上以附图标记350示意性地示出,并且形成包含反射激光的脉冲的反射图案。反射图案作为多个光斑在基于cmos的多阱像素检测器上被观察到。此类光斑必须被解释为激光脉冲与检测器的检测表面的交点。
[0205]
根据本发明的基于cmos的多阱像素检测器被配置并可操作用于在落在脉冲激光束的脉冲周期p
p
内的检测时间段td期间检测反射激光,并且其中检测时间段被划分成m个连贯检测时间窗口tw[i],使得其中m是检测时间窗口的数量。这一数量m等于α
×
(n+1)且m≥3个,优选地m≥4,更优选地m≥5。在各实施例中,对于i=1到m,连贯检测时间窗口tw[i]的持续时间基本相等。因子α是上文在讨论根据本发明的方法时讨论的整数,其中α≥1,优选地α≤10,更优选地α≤5。
[0206]
在各实施例中,检测时间窗口几乎等于脉宽,即对于每一检测时间窗口tw[i],其中i=1到m:(|pw-tw[i]|/pw)≤0.10,优选地(|pw-tw[i]|/pw)≤0.05,更优选地(|pw-tw[i]|/pw)≤0.02。
[0207]
根据本发明的基于cmos的多阱像素检测器被配置并可操作用于针对每一检测到的反射激光的光斑,获得与检测时间窗口tw[i]相关联的曝光值qi。通过累积代表在检测时间窗口tw[i]中检测到的反射激光的电荷量来获得曝光值,使得在检测到激光与在其中检测到激光的时间窗口之间建立相关性。对脉冲序列中的所有脉冲执行电荷的累积。
[0208]
对反射激光的检测要与激光脉冲序列的发射同步地执行。因此,lidar设备包括控制器200,该控制器用于控制投影仪100和光接收设备300,以便与所发射的脉冲序列同步地检测并累积反射激光。在各实施例中,控制器包括同步装置,其可包括常规时钟电路或振荡器。
[0209]
在各实施例中,控制器200被配置成用于执行同步,使得第一检测时间窗口tw[1]与脉冲的脉冲发射时间窗口交叠或至少部分交叠。脉冲的发射时间窗口必须被解释为在其中发射脉冲的时间窗口,并且因此该发射时间窗口具有脉宽的宽度。上文已经在讨论根据
本发明的方法时讨论了与脉冲的发射时间窗口交叠或部分交叠的第一检测时间窗口的示例。
[0210]
在其他实施例中,控制器200被配置成用于执行脉冲的发射与检测器之间的同步,使得第一检测时间窗口tw[1]相对于脉冲发射时间窗口被延迟,并且该延迟等于预定义延迟时间窗口t
dl

[0211]
优选地,预定延迟时间因子等于或小于2
×dmin
/c,其中d
min
是最小距离,使得可以检测位于最小距离d
min
之处或附近的对象并确定距离。
[0212]
基于cmos的多阱像素检测器必须被解释为像素的2d阵列,每一像素形成用于检测激光的个体检测器。事实上,每一像素包括个体光敏区和至少三个电荷存储阱。电荷存储阱的数量取决于检测时间窗口的数量m。
[0213]
更一般而言,基于cmos的多阱像素检测器的每一像素包括多个电荷存储阱w[i],其中i=1到nw,nw是电荷存储阱的数量,其中nw≤m。
[0214]
电荷存储阱是用于累积得自落在像素的光敏区上的激光的电荷的储存库。每个阱都具有由外部信号控制的单独转移栅,该外部信号与激光脉冲同步并与所指定的n+1个检测时间窗口同步。以此方式,通过控制转移栅,在给定时间窗口tw[i]中由像素的光敏区检测到的反射激光将被存储在相关联的阱w[k]中。
[0215]
光接收设备还包括用于将反射激光投射到多像素检测器上的透镜系统。例如,如果照明图案的所有脉冲被场景的对象反射,则在光接收设备的焦平面中的、多像素检测器的光敏区所在的位置处获得图案的图像。
[0216]
根据本发明的lidar设备包括处理装置400,其被配置成对于每一检测到的反射激光的光斑,标识在连贯检测时间窗口中的哪些检测时间窗口中检测到反射激光,即已获得曝光值。反射激光在α或者α+1个连贯检测时间窗口中被检测到。处理装置还基于对检测时间窗口的标识和/或基于在所标识的检测时间窗口中获得的最早和最后曝光值来计算对象距离ds。上面讨论的用于执行标识的方法通常被实现为用于执行该标识的计算机算法。除该标识算法外,又一计算机算法可被用于如上讨论地使用该公式和用于执行这一计算的各种选项来计算距离ds。在各实施例中,处理装置400通常包括处理器、微处理器或计算机。
[0217]
注意,当多像素检测器检测到从场景反射的激光的光斑时,该反射激光是在多像素检测器中的多个像素中被检测到的,并且为了确定反射激光的量,通常对由在多个像素中检测到的所检测激光进行求和。在各实施例中,当执行求和时,与像素相关联的权重因子被纳入考虑。当使用在空间上分开的光斑的图案执行照明时,在多像素检测器上检测到的各种光斑也在空间上分开,使得在检测到的各种光斑之间没有交叠或极少交叠。
[0218]
当确定反射激光的曝光值时,处理装置被配置成减去背景。该背景通常从其中未检测到反射激光的检测时间窗口获得。
[0219]
在一些实施例中,在读出数据之前,在cmos检测器的水平上执行已检测到激光的时间窗口的确定。在这些实施例中,基于cmos的多阱像素检测器包括例如比较器,该比较器被配置成用于区分背景计数和激光。
[0220]
像素检测器中每一像素的阱w[k]的数量nw不一定等于时间窗口tw[i]的数量m。在优选实施例中,nw《m,并且因此在这些实施例中,由于存储阱少于时间窗口,所以多个电荷存储阱中的至少一个电荷存储阱被配置成用于存储在检测时间窗口tw[i]的两个不同时间
窗口期间检测到的电荷。如上所讨论的,这两个不同时间窗口是两个非连贯时间窗口。
[0221]
在其中m≥4且α=1的实施例中,对于多阱像素检测器nw《m的每一像素,每对两个连贯检测时间窗口与唯一性的一对电荷存储阱相关联。
[0222]
更一般而言,对于α≥1,α+1个连贯检测时间窗口tw[k]到tw[k+α]中的每一子序列与该多个电荷存储阱w[i]中的α+1个电荷存储阱的唯一性序列相关联,使得在其中检测到反射激光的检测时间窗口是能毫无疑义地标识的,其中k=1到m-α。上文已经讨论了电荷储存阱的唯一性序列的定义。
[0223]
如果α=1,则对于给定数量nw的电荷存储阱,可以作出的唯一性的对组合p的数量可以表示为排列p的数量,并且可以如下计算:p=nw!/((n
w-2)!
×
2)。
[0224]
符号“!”是排列符号。例如,对于三个阱,可以形成三个唯一性的对组合,即1-2、1-3和2-3。对于例如4个阱,可以形成6个唯一性组合:即1-2、1-3、1-4、2-3、2-4和3-4。对于五个阱,可以形成10个唯一性的组合,并且因此当使用例如11个检测时间窗口时,只需要5个阱。在图7中,示出了一示例,其中脉宽缩减因子n等于10,并且因此根据本发明提供了11个检测时间窗口。在图4的顶部,示出了一示例,其中每一检测时间窗口具有其适当的阱,并且因此每一像素将需要11个阱。在图7的底部,只有5个阱被用于形成11个检测时间窗口的10个唯一性的对组合。有利地,通过形成唯一性的阱对,所需的阱的数量大大减少,并且基于cmos的多阱像素检测器的复杂度降低。
[0225]
对于α=1且对于给定数量nw的存储阱,还可以计算可被应用于减小脉宽的最大脉宽缩减因子n
max
,如上所讨论的。如果阱的数量是奇数,则可以应用以下公式来计算能与给定数量的阱一起应用的最大缩减因子:n
max
=nw!/((n
w-2)!
×
2)。
[0226]
例如,如果阱数是5,则n
max
=10。这也意味着可被用于给定数量的阱的检测时间窗口的最大数量m
max
可表示为:m
max
=1+nw!/((n
w-2)!
×
2)。
[0227]
如果α=1且阱数nw是偶数,则用于计算最大脉宽缩减因子n
max
的公式是不同的且可如下表示:n
max
=nw!/((n
w-2)!
×
2)

((nw!/((n
w-1)!
×
2))-1)。
[0228]
例如,如果阱数是4,则n
max
=5。类似地,在阱数是偶数的情形中,检测时间窗口的最大数量m
max
可以表示为:m
max
=1+nw!/((n
w-2)!
×
2)

((nw!/((n
w-1)!
×
2))-1)。
[0229]
在各实施例中,5≤m≤6且nw=4。在其他实施例中,7≤m≤11且nw=5。在又一实施例中,m=4且nw=3。
[0230]
更一般而言,如果α》1,已知的数学算法可被用于定义可被用于毫无疑义地能够标识在其中检测到反射激光的检测时间窗口的、电荷存储阱的唯一性序列的最大数量c
max

[0231]
在优选实施例中,根据本发明的lidar设备被配置成使得脉宽pw与检测时间窗口的数量m相关,并且该关系可以表示为:
[0232]
pw=(tof
max

t
dl
)/((m-α)/α),其中tof
max
=2
×dmax
/c,其中tof
max
是最大飞行时间tof
max
且d
max
是与最大飞行时间相关联的最大距离,c是光速,且t
dl
是上述预定义延迟时间窗口,其中0≤t
dl
≤(0.6
×
tof
max
),优选地0≤t
dl
≤(0.4
×
tof
max
),更优选地0≤t
dl
≤(0.2
×
tof
max
)。
[0233]
如上所述,d
max
将被解释为lidar设备的工作范围,通常是最大工作范围。在各实施例中,根据本发明的lidar设备的最大工作范围d
max
例如在以下范围内:30米≤d
max
≤1000米。
[0234]
遵循上述定义的脉宽pw通常预期在纳秒到微秒范围内。假设例如最大范围是100米和300米,则激光来回的相应飞行时间分别是666纳秒和2微秒,并且因此取决于所定义的脉宽缩减因子n,或时间检测窗口的等效数量,并取决于lidar设备的所需最大工作距离d
max
,脉冲长度pw预期在33纳秒到1微秒之间。具有缩减因子的脉冲通常等于或小于20。
[0235]
在各实施例中,lidar设备的投影仪被配置成用于脉冲频率f
p
提供脉冲激光束,使得f
p
≤1/((m/α)
×
pw+t
dl
)。
[0236]
根据本发明的lidar设备适于集成到交通工具中。集成到交通工具中的lidar设备被布置成可操作地覆盖该交通工具周围区域的至少一部分。区域的该至少一部分对应于需要确定距离的场景。所覆盖的区域取决于lidar设备的视野(fov),并且在各实施例中,fov是例如30
°×
10
°
或120
°×
30
°
或任何其他合适视野。
当前第1页1 2 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1