CN112946675B - 一种基于时间融合的距离测量方法、系统和设备 - Google Patents
一种基于时间融合的距离测量方法、系统和设备 Download PDFInfo
- Publication number
- CN112946675B CN112946675B CN202110120017.0A CN202110120017A CN112946675B CN 112946675 B CN112946675 B CN 112946675B CN 202110120017 A CN202110120017 A CN 202110120017A CN 112946675 B CN112946675 B CN 112946675B
- Authority
- CN
- China
- Prior art keywords
- distance
- error
- time
- photon
- histogram
- 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.)
- Active
Links
- 230000004927 fusion Effects 0.000 title claims abstract description 28
- 238000000691 measurement method Methods 0.000 title claims abstract description 21
- 238000005259 measurement Methods 0.000 claims abstract description 53
- 238000000034 method Methods 0.000 claims abstract description 18
- 238000012545 processing Methods 0.000 claims description 15
- 238000004590 computer program Methods 0.000 claims description 8
- 238000013139 quantization Methods 0.000 claims description 8
- 230000004044 response Effects 0.000 claims description 5
- 230000003068 static effect Effects 0.000 claims description 2
- 238000001514 detection method Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000002829 reductive effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007500 overflow downdraw method Methods 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/06—Systems determining position data of a target
- G01S17/08—Systems determining position data of a target for measuring distance only
- G01S17/10—Systems determining position data of a target for measuring distance only using transmission of interrupted, pulse-modulated waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/495—Counter-measures or counter-counter-measures using electronic or electro-optical means
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明公开了一种基于时间融合的距离测量方法、系统和设备,方法包括:根据直方图计算光子从发射到被接收的飞行时间,并根据飞行时间计算目标的距离;其中,所述直方图包括有连续的时间间隔,所述时间间隔内包含由发射器发出的脉冲光束被目标反射后、经由采集器采集的脉冲光束中光子的计数值;计算测量目标距离时产生的测量误差,并根据所述测量误差设定误差阈值;计算当前帧的测量距离与历史帧的测量距离的距离差值;对比所述距离差值与误差阈值,并根据对比结果输出最终的目标距离。
Description
技术领域
本发明涉及测距技术领域,具体涉及一种基于时间融合的距离测量方法、系统和设备。
背景技术
利用飞行时间原理(TOF,Time of Flight)可以对目标进行距离测量以获取包含目标的深度值的深度图像,而基于飞行时间原理的距离测量系统已被广泛应用于消费电子、无人架驶、AR/VR等领域。基于飞行时间原理的距离测量系统通常包括发射器和采集器,利用发射器发射脉冲光束照射目标视场并利用采集器采集反射光束中的光子,计算光子由发射到反射接收的飞行时间来计算物体的距离。其中,时间数字转换器(TDC)用于记录光子从发射到被采集的飞行时间并生成光子信号,利用该光子信号寻找直方图电路中相应的时间bin(时间间隔),使得该时间间隔内的光子计数值加1,当进行大量重复脉冲探测后,可以统计得到时间信号对应的光子计数的直方图,根据直方图计算出光子从发射到接收的飞行时间进一步计算出物体的距离。
对于基于光子计数飞行时间的距离测量系统,其优点是具有高灵敏度,即使单光子量级的回波也能响应引起雪崩,但是当回波信号为单光子量级时,由于光子的量子性,回波中的信号光子随机地分布在信号脉冲内,导致单次光子探测的时刻分布具有随机性,因此需要通过光子计数(TCSPC)来减少随机误差(测量误差)。
当回波信号较弱时,实际接收到的光子数有限,信号光子的泊松噪声抖动成为影响测距精度的主要因素;当环境光较强时,环境光造成的大量环境光子的泊松噪声抖动也会影响测距精度。而泊松噪声区别与常见的高斯噪声,常用的滤波算法效果有限,降低泊松噪声抖动的方法包括增加重复探测次数、多次测量距离求均值等,以实现降低测量误差从而提高测距精度,但是目前的这些方法均会对系统的帧频会造成影响。
发明内容
有鉴于此,本发明提出一种基于时间融合的距离测量方法、系统和设备,以解决上述现有技术所存在的至少一种技术问题。
为了解决上述技术问题,本发明的一种技术方案为:
一种基于时间融合的距离测量方法,包括:根据直方图计算光子从发射到被接收的飞行时间,并根据飞行时间计算目标的距离;其中,所述直方图包括有连续的时间间隔,所述时间间隔内包含由发射器发出的脉冲光束被目标反射后、经由采集器采集的脉冲光束中光子的计数值;计算测量所述目标的距离时产生的测量误差,并根据所述测量误差设定误差阈值;计算当前帧的测量距离与历史帧的测量距离的距离差值;对比所述距离差值与误差阈值,并根据对比结果输出最终的目标距离。
在一些实施例中,所述测量误差根据单光子分布的抖动误差和探测到的信号光子数计算得到,具体公式为:
其中,σt为测距标准差,用于表征所述测量误差;σeach为单光子分布的抖动误差;Ns表征探测到的信号光子数,根据直方图中与脉冲光束对应的回波区间内的光子计数总和确定。
在一些实施例中,所述信号光子数等于所述回波区间内的光子计数总和减去所述回波区间内的噪声光子总数;其中,所述回波区间内的噪声光子总数等于直方图中每个时间间隔内包含的噪声光子数乘以所述回波区间内包含的时间间隔数量。
在一些实施例中,所述噪声光子数通过如下方式计算:从直方图中选取远离脉冲峰值位置的一局部区域;将所述局部区域内的光子计数总值按照该局部区域内时间间隔的数量求均值,记为所述噪声光子数;或,所述噪声光子数通过如下方式计算:选取直方图中除脉冲位置处以外的区域,并将该区域内的光子计数总值按照该区域内时间间隔的数量求均值,记为所述噪声光子数。
在一些实施例中,所述单光子分布的抖动误差通过如下方式计算:
其中,Tpulse表示脉冲持续时间,即脉宽;K表示量化参数;
或者,所述单光子分布的抖动误差通过误差叠加的方法计算得到;其中,进行误差叠加的误差项包括信号光子分布的抖动误差、发射脉冲波形的抖动误差、单光子响应抖动误差、噪声光子计数误差、雪崩脉冲波形抖动误差、TDC计时抖动误差或直方图累计的量化误差。
在一些实施例中,根据所述测量误差设定误差阈值Th,令Th=εσt;其中,σt表示测量误差,ε为预先设定的线性控制参数。
在一些实施例中,所述计算当前帧的测量距离与历史帧的测量距离的距离差值,包括:
ΔD=Dcurrent-Dprevious
其中,ΔD表示所述距离差值,Dcurrent为当前帧的测量距离,Dprevious为历史帧的测量距离;优选地,历史帧即为当前帧的前一帧。
在一些实施例中,所述对比所述距离差值与误差阈值,并根据对比结果输出最终的目标距离,包括:当所述距离差值小于或等于所述误差阈值时,对当前帧的测量距离与历史帧的测量距离进行加权融合后输出最终的目标距离;当所述距离差值大于所述误差阈值时,直接输出当前帧的测量距离作为最终的目标距离,并将历史帧的测量距离更新为当前帧的测量距离。
本发明的另一技术方案为:
一种基于时间融合的距离测量系统,包括:发射器,用于朝向目标发射脉冲光束;采集器,用于采集被目标反射回的所述脉冲光束中的光子并形成光子信号;处理电路,与所述发射器以及所述采集器连接,用于对所述光子信号进行处理以形成直方图,并执行前述的基于时间融合的距离测量方法的步骤。
本发明的另一技术方案为:
一种基于时间融合的距离测量设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现前述的基于时间融合的距离测量方法。
本发明技术方案的有益效果是:根据距离差值与误差阈值的对比来灵活输出最终的目标距离,有效的减少了在连续帧测量过程中由光子泊松噪声抖动产生的测量误差,提升测距精度的同时不会对系统的帧频产生影响。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例的距离测量系统的示意图;
图2是本发明实施例的距离测量方法的流程图;
图3是本发明实施例的一种示例性的直方图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或成一体;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
另需要理解的是,术语“长度”、“宽度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”和“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
图1所示为本发明一个实施例基于时间融合的距离测量系统的示意图。距离测量系统10包括发射器11、采集器12以及处理电路13。其中,发射器11包括由一个或多个激光器组成的光源111,用于向目标20发射脉冲光束30,至少部分脉冲光束经过目标20反射形成反射光束40回到采集器12;采集器12包括由多个像素组成的像素阵列121,用于采集反射光束40中的光子并输出光子信号,处理电路13同步发射器11与采集器12的触发信号以计算光束中的光子从发射到接收所需要的飞行时间。
发射器11包括光源111、发射光学元件112以及驱动器113等。在一个实施例中,光源111是在单块半导体基底上生成多个VCSEL光源以形成的VCSEL阵列光源芯片。其中,光源111可以在驱动器113的驱动控制下以一定频率(脉冲周期)向外发射脉冲光束,脉冲光束经过发射光学元件112投射到目标场景上形成照明斑点,其中,发射频率根据测量距离进行设定。
采集器12包括像素阵列121、过滤单元122和接收光学元件123等。接收光学元件123将目标反射的斑点光束成像到像素阵列121上,像素阵列121包括多个采集光子的像素,所述像素可以是APD、SPAD、SiPM等采集光子的单光子器件中的一种,像素阵列121采集到光子的情况被视为光子检测事件发生并输出光子信号。在一个实施例中,像素阵列121包括由多个SPAD组成,SPAD可以对入射的单个光子进行响应并输出指示所接收光子在每个SPAD处相应到达时间的光子信号。一般地,还包括有与像素阵列连接的信号放大器、时数转换器(TDC)、数模转换器(ADC)等器件中的一种或多种组成的读出电路(图中未示出)。这些电路既可以与像素整合在一起,作为采集器的一部分,也可以作为处理电路13的一部分,后面为便于描述,将统一视作处理电路13的一部分。
处理电路13用于接收光子信号并进行处理计算出光子从发射到接收的飞行时间,并进一步计算出目标的距离信息。在一个实施例中,处理电路13包括TDC电路以及直方图存储器,TDC电路接收光子信号用于确定光子从发射到采集的飞行时间,并生成表征飞行时间信息的时间码,利用时间码寻找直方图存储器中的对应位置,并使得直方图存储器的对应位置处存储的数值加“1”,根据直方图存储器的位置作为时间bin(时间间隔)构造直方图。
图2所示是本发明一个实施例一种基于时间融合的距离测量方法的流程图。该距离测量方法包括如下步骤:
S1、根据直方图计算光子从发射到接收的飞行时间,并根据飞行时间计算目标距离;其中,所述直方图包括连续的时间间隔,所述时间间隔内包含由发射器发出的脉冲光束被目标反射后、经由采集器采集的所述脉冲光束中光子的计数值。
结合图1所示,在距离测量系统中,处理电路13控制发射器11朝向目标区域发射脉冲光束,被目标反射的部分脉冲光束入射到采集器12,采集器12采集反射的脉冲光束中的光子并生成包含光子的飞行时间的光子信号;处理电路13接收光子信号并进行处理形成直方图,直方图包括连续的时间间隔,每个时间间隔用于表示在检测时段内采集器采集到光子的计数值。处理电路13对直方图进行处理后,根据直方图计算出光子从发射到被接收的飞行时间,其中采用的计算方法包括寻峰法、质心法、模板匹配法等。在本发明的实施例中,根据直方图计算光子飞行时间的方法不做具体限制。
图3所示是本发明一个实施例的直方图的示意图。直方图可以被称为检测数据,用于表示在检测时段内采集器采集到光子的时间分布。一般地,时间间隔301的大小为十几到几十皮秒,则一个脉冲光束的光子信号在直方图中对应分布在一个回波区间内,该回波区间包括多个时间间隔,时间间隔的数量根据脉冲光束的脉冲持续时间(脉宽)确定,比如一个脉冲光束在直方图中对应的脉冲位置为回波区间302。则可以选取脉冲峰值位置所处时间间隔的时间作为脉冲光束的飞行时间,一般选择时间间隔的中间量为该时间间隔的时间。在本发明实施例中,假设脉冲光束的脉宽为2ns,直方图中一个时间间隔301的大小为100ps,则一个脉冲的光子信号在直方图中对应分布在连续的20个时间间隔内,即在直方图中脉冲位置处的时间间隔的数量为20。
S2、计算测量目标距离时产生的测量误差,并根据所述测量误差设定误差阈值。
在距离测量系统进行目标距离测量时会产生测量误差,测量误差主要受到单光子分布的抖动误差影响。由于距离测量过程中不同光子之间是相互独立,则采用TCSPC(时间相关单光子计数)技术能够使得抖动误差减小,因此在计算测量误差时,可以先计算出单光子分布的抖动误差,再根据单光子分布的抖动误差计算出测量误差。其中,单光子分布的抖动误差受到多种因素的影响,包括系统参数、环境光以及信号光强度等等。
在一个实施例中,单光子分布的抖动误差主要受到信号光子分布的抖动误差影响,信号光子分布的抖动误差根据发射脉冲光束的宽度(简称脉宽,即脉冲持续时间)计算,则将信号光子分布的抖动误差记为单光子分布的抖动误差,具体公式为:
其中,σeach为单光子分布的抖动标准差,用于表征单光子分布的抖动误差;Tpulse表征脉冲持续时间。在一个实施例中,根据脉冲半峰宽确定脉冲持续时间,K表征量化参数,可通过理论分析或者实验测量计算确定。
进一步的,再根据单光子分布的抖动误差计算测量误差,计算公式为:
其中,σt为测距标准差,用于表征测量误差;Ns表征探测到的信号光子数,可以根据直方图中与脉冲光束对应的回波区间内全部时间间隔内光子计数总和确定;优选地,所述信号光子数等于该回波区间内的光子计数总和减去该回波区间内的噪声光子总数。
回波区间内的噪声光子总数根据直方图中每个时间间隔内包含的噪声光子数确定,具体而言,回波区间内的噪声光子总数等于直方图中每个时间间隔内包含的噪声光子数乘以所述回波区间内包含的时间间隔数量。在一个实施例中,从直方图中截取局部区域计算噪声光子数,具体的,根据直方图中脉冲峰值位置选择远离脉冲峰值位置的局部区域用于计算噪声光子数。比如以直方图中间位置的时间间隔为分界线,若脉冲峰值位置处于直方图的后半部分,则在直方图的前半部分选择局部区域计算噪声光子数,即选定局部区域后,将该局部区域内的光子计数总值按照该局部区域内时间间隔的数量求均值,记为所述噪声光子数;同理,若位于前半部分,则从后半部分选择局部区域计算噪声光子数。在一个实施例中,根据初始直方图全部的时间间隔计算噪声光子数,即选取直方图中除脉冲位置处以外的区域,并将该区域内的光子计数总值按照该区域内时间间隔的数量求均值,记为所述噪声光子数。
在一个实施例中,影响单光子的抖动误差的影响因素还包括发射脉冲波形的抖动误差、单光子响应抖动误差、噪声光子计数误差、雪崩脉冲波形抖动误差、TDC计时抖动误差、直方图累计的量化误差等等,因此,采用误差叠加方法计算出单个光子的抖动误差,具体计算公式为:
其中,对应一种存在的误差情况,即利用多个误差的平方和开根号计算出单个光子的抖动误差。在误差叠加公式中,至少包括两种误差,其中一种误差为信号光子分布的抖动误差。在距离测量系统中,上述的多种误差值可以根据实际测量的情况确定或者可以根据系统参数确定。
具体的,发射脉冲波形的抖动误差受到光源及驱动器性能的影响,可根据系统中光源及驱动器的参数直接确定出发射脉冲波形的抖动误差。单光子响应抖动误差受到光子响应上升沿速率以及斯密特触发器阈值的影响,与采集器的设计参数相关。噪声光子计数误差受到环境光噪声强度的影响,根据直方图计算噪声光子数,并根据噪声光子数计算出对应的噪声光子计数误差。TDC计时抖动误差主要受到TDC电路本身的性能参数、TDC的计时精度等特性的影响,可根据TDC电路的具体设计通过理论分析或者实验测试确定TDC计时抖动误差。直方图累计的量化误差与直方图中时间bin的宽度有关,当时间bin的宽度大于TDC的时间分辨率时,在将每个光子的计数结果量化到每个具体的时间bin内时会产生非线性误差,可以根据直方图中时间bin的宽度计算出直方图累计的量化误差。
根据计算出的测量误差σt设定误差阈值Th,令Th=εσt,其中,σt用于表征测量误差;ε为线性控制参数,可以根据不同场景进行设定。
S3、计算当前帧的测量距离与历史帧的测量距离的距离差值;
在距离测量过程中,每测量一帧都可以得到一个距离值,但由于测量误差的存在,则可能导致多帧距离值之间存在测量误差,因此基于测量时间的连续性采用时间融合的方法修正测量误差。
具体的,计算当前帧的测量距离与历史帧的测量距离的距离差值ΔD,计算公式为:
ΔD=Dcurrent-Dprevious
其中,Dcurrent为当前帧的测量距离,Dprevious为历史帧的测量距离,优选地,历史帧即为当前帧的前一帧。
S4:对比所述距离差值与误差阈值,并根据对比结果输出目标距离。
根据步骤S3中计算出来的距离差值ΔD与误差阈值Th进行对比,并根据对比结果输出目标距离。当距离差值ΔD小于或等于误差阈值Th时,则可以认为距离测量系统与目标间相对静止,没有产生明显的空间移动,偏差主要来自于测量误差,则对当前帧测量距离与历史帧测量距离进行加权融合后输出最终的目标距离Doutput:
Doutput=α*Dcurrent+(1-α)*Dprevious
其中,Doutput表征目标距离,α为进行融合的权重。并且将历史帧的测量距离更新为Doutput,用于作为下一帧数据融合的参考。
在实际应用中,也会存在距离差值ΔD大于误差阈值Th的情况,此时则可以认为目标产生了明显的空间移动,包括沿着深度方向或者移出视场范围的方向移动,则此时距离差值不再主要受到测量误差的影响,此时直接输出当前帧的测量距离Dcurrent作为目标距离,并将历史帧的测量距离更新为Dcurrent。
作为本发明另一实施例,还提供一种基于时间融合的距离测量设备,包括:存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序;其中,所述处理器执行所述计算机程序时实现前述实施例所述的距离测量方法。
本发明的实施例可以包括或利用包括计算机硬件的专用或通用计算机,如下面更详细讨论的。在本发明的范围内的实施例还包括用于携带或存储计算机可执行指令和/或数据结构的物理和其他计算机可读介质。这样的计算机可读介质可以是可以被通用或专用计算机系统访问的任何可用介质。存储计算机可执行指令的计算机可读介质是物理存储介质。携带计算机可执行指令的计算机可读介质是传输介质。因此,作为示例而非限制,本发明的实施例可以包括至少两种截然不同的计算机可读介质:物理计算机可读存储介质和传输计算机可读介质。
本申请实施例还提供一种计算机设备,所述计算机设备包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序,其中,所述处理器执行所述计算机程序时至少实现前述实施例方案中所述的距离测量方法。
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述或记载的部分,可以参见其它实施例的相关描述。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于时间融合的距离测量方法,其特征在于,包括:
根据直方图计算光子从发射到被接收的飞行时间,并根据飞行时间计算目标的距离;其中,所述直方图包括有连续的时间间隔,所述时间间隔内包含由发射器发出的脉冲光束被目标反射后、经由采集器采集的脉冲光束中光子的计数值;
基于单光子分布的抖动误差计算测量所述目标的距离时产生的测量误差,并根据所述测量误差设定误差阈值;
计算当前帧的测量距离与历史帧的测量距离的距离差值;
对比所述距离差值与所述误差阈值,并根据对比结果确定与所述目标之间相对静止或产生了明显的空间移动,输出最终的目标距离。
2.如权利要求1所述的基于时间融合的距离测量方法,其特征在于,所述测量误差根据单光子分布的抖动误差和探测到的信号光子数计算得到,具体公式为:
其中,σt为测距标准差,用于表征所述测量误差;σeach为单光子分布的抖动误差;Ns表征探测到的信号光子数,根据直方图中与脉冲光束对应的回波区间内的光子计数总和确定。
3.如权利要求2所述的基于时间融合的距离测量方法,其特征在于,所述信号光子数等于所述回波区间内的光子计数总和减去所述回波区间内的噪声光子总数;其中,所述回波区间内的噪声光子总数等于直方图中每个时间间隔内包含的噪声光子数乘以所述回波区间内包含的时间间隔数量。
4.如权利要求3所述的基于时间融合的距离测量方法,其特征在于,所述噪声光子数通过如下方式计算:
从直方图中选取远离脉冲峰值位置的一局部区域;
将所述局部区域内的光子计数总值按照该局部区域内时间间隔的数量求均值,记为所述噪声光子数;
或,所述噪声光子数通过如下方式计算:
选取直方图中除脉冲位置处以外的区域,并将该区域内的光子计数总值按照该区域内时间间隔的数量求均值,记为所述噪声光子数。
5.如权利要求2所述的基于时间融合的距离测量方法,其特征在于,所述单光子分布的抖动误差通过如下方式计算:
其中,Tpulse表示脉冲持续时间,即脉宽;K表示量化参数;
或者,
所述单光子分布的抖动误差通过误差叠加的方法计算得到;其中,进行误差叠加的误差项包括信号光子分布的抖动误差、发射脉冲波形的抖动误差、单光子响应抖动误差、噪声光子计数误差、雪崩脉冲波形抖动误差、TDC计时抖动误差或直方图累计的量化误差。
6.如权利要求1所述的基于时间融合的距离测量方法,其特征在于,根据所述测量误差设定误差阈值Th,令Th=εσt;其中,σt表示测量误差,ε为预先设定的线性控制参数。
7.如权利要求1所述的基于时间融合的距离测量方法,其特征在于,所述计算当前帧的测量距离与历史帧的测量距离的距离差值,包括:
ΔD=Dcurrent-Dprevious
其中,ΔD表示所述距离差值,Dcurrent为当前帧的测量距离,Dprevious为历史帧的测量距离;优选地,历史帧即为当前帧的前一帧。
8.如权利要求1所述的基于时间融合的距离测量方法,其特征在于,所述对比所述距离差值与误差阈值,并根据对比结果输出最终的目标距离,包括:
当所述距离差值小于或等于所述误差阈值时,对当前帧的测量距离与历史帧的测量距离进行加权融合后输出最终的目标距离;
当所述距离差值大于所述误差阈值时,直接输出当前帧的测量距离作为最终的目标距离,并将历史帧的测量距离更新为当前帧的测量距离。
9.一种基于时间融合的距离测量系统,其特征在于,包括:
发射器,用于朝向目标发射脉冲光束;
采集器,用于采集被目标反射回的所述脉冲光束中的光子并形成光子信号;
处理电路,与所述发射器以及所述采集器连接,用于对所述光子信号进行处理以形成直方图,并执行权利要求1-8任一项所述的基于时间融合的距离测量方法的步骤。
10.一种基于时间融合的距离测量设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于:所述处理器执行所述计算机程序时实现权利要求1-8任一项所述的基于时间融合的距离测量方法。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110120017.0A CN112946675B (zh) | 2021-01-28 | 2021-01-28 | 一种基于时间融合的距离测量方法、系统和设备 |
PCT/CN2021/106027 WO2022160611A1 (zh) | 2021-01-28 | 2021-07-13 | 一种基于时间融合的距离测量方法、系统和设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110120017.0A CN112946675B (zh) | 2021-01-28 | 2021-01-28 | 一种基于时间融合的距离测量方法、系统和设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112946675A CN112946675A (zh) | 2021-06-11 |
CN112946675B true CN112946675B (zh) | 2023-12-01 |
Family
ID=76238837
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110120017.0A Active CN112946675B (zh) | 2021-01-28 | 2021-01-28 | 一种基于时间融合的距离测量方法、系统和设备 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN112946675B (zh) |
WO (1) | WO2022160611A1 (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112946675B (zh) * | 2021-01-28 | 2023-12-01 | 深圳奥锐达科技有限公司 | 一种基于时间融合的距离测量方法、系统和设备 |
CN115657056A (zh) * | 2021-07-07 | 2023-01-31 | 奥比中光科技集团股份有限公司 | 一种距离测量系统及其相对精度的确定方法、装置、设备 |
CN116009010A (zh) * | 2021-10-21 | 2023-04-25 | 华为技术有限公司 | 一种探测方法及装置 |
CN114236504A (zh) * | 2021-12-17 | 2022-03-25 | 奥诚信息科技(上海)有限公司 | 一种基于dToF的探测系统及其光源调整方法 |
CN115128623A (zh) * | 2022-04-13 | 2022-09-30 | 深圳市灵明光子科技有限公司 | 激光测距方法和装置 |
CN114509740B (zh) * | 2022-04-18 | 2022-08-09 | 深圳阜时科技有限公司 | 飞行时间偏移的修正方法、ToF装置、电子设备及存储介质 |
CN114637019B (zh) * | 2022-05-18 | 2022-08-30 | 杭州宇称电子技术有限公司 | 一种基于时间分段自适应计数量化的抗环境光方法 |
CN115546079B (zh) * | 2022-11-25 | 2023-03-10 | 杭州宇称电子技术有限公司 | Tof直方图动态范围拓展方法及其应用 |
CN116626650B (zh) * | 2023-07-26 | 2023-11-10 | 中国人民解放军国防科技大学 | 一种单光子探测的运动目标全波形恢复方法与系统 |
CN118294972A (zh) * | 2024-03-22 | 2024-07-05 | 中国电子科技集团公司第三十八研究所 | 基于光子计数激光雷达的自适应测距方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008120532A1 (ja) * | 2007-03-29 | 2008-10-09 | Brother Kogyo Kabushiki Kaisha | 移動局の静止判定を利用した位置検出システム |
CN109696654A (zh) * | 2017-10-23 | 2019-04-30 | 深圳市优必选科技有限公司 | 定位方法及终端设备 |
WO2019136612A1 (zh) * | 2018-01-09 | 2019-07-18 | 深圳市沃特沃德股份有限公司 | 基于距离度量的视觉词典闭环检测方法与装置 |
CN110082772A (zh) * | 2019-05-05 | 2019-08-02 | 中国科学院国家天文台长春人造卫星观测站 | 一种信号回波率实时可控的卫星激光测距系统、方法及装置 |
WO2021000313A1 (zh) * | 2019-07-04 | 2021-01-07 | 深圳市大疆创新科技有限公司 | 基于侧向毫米波雷达的侧向静止物体和移动速度检测方法 |
CN112255636A (zh) * | 2020-09-04 | 2021-01-22 | 奥诚信息科技(上海)有限公司 | 一种距离测量方法、系统及设备 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DK3206045T3 (da) * | 2016-02-15 | 2021-06-21 | Airborne Hydrography Ab | Enkelt foton lidar scanner |
CN105954733B (zh) * | 2016-06-17 | 2018-11-13 | 南京理工大学 | 基于光子飞行时间相关性的时域滤波方法 |
EP3370080B1 (en) * | 2017-03-01 | 2021-04-28 | STMicroelectronics (Grenoble 2) SAS | Range and parameter extraction using processed histograms generated from a time of flight sensor - parameter extraction |
WO2020150574A1 (en) * | 2019-01-18 | 2020-07-23 | Sense Photonics, Inc. | Digital pixels and operating methods thereof |
US20200341144A1 (en) * | 2019-04-26 | 2020-10-29 | Ouster, Inc. | Independent per-pixel integration registers for lidar measurements |
CN110954919B (zh) * | 2019-12-13 | 2022-03-18 | 华中科技大学 | 一种面阵激光探测器的固定值噪声确定方法及去除方法 |
CN112946675B (zh) * | 2021-01-28 | 2023-12-01 | 深圳奥锐达科技有限公司 | 一种基于时间融合的距离测量方法、系统和设备 |
-
2021
- 2021-01-28 CN CN202110120017.0A patent/CN112946675B/zh active Active
- 2021-07-13 WO PCT/CN2021/106027 patent/WO2022160611A1/zh active Application Filing
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008120532A1 (ja) * | 2007-03-29 | 2008-10-09 | Brother Kogyo Kabushiki Kaisha | 移動局の静止判定を利用した位置検出システム |
CN109696654A (zh) * | 2017-10-23 | 2019-04-30 | 深圳市优必选科技有限公司 | 定位方法及终端设备 |
WO2019136612A1 (zh) * | 2018-01-09 | 2019-07-18 | 深圳市沃特沃德股份有限公司 | 基于距离度量的视觉词典闭环检测方法与装置 |
CN110082772A (zh) * | 2019-05-05 | 2019-08-02 | 中国科学院国家天文台长春人造卫星观测站 | 一种信号回波率实时可控的卫星激光测距系统、方法及装置 |
WO2021000313A1 (zh) * | 2019-07-04 | 2021-01-07 | 深圳市大疆创新科技有限公司 | 基于侧向毫米波雷达的侧向静止物体和移动速度检测方法 |
CN112255636A (zh) * | 2020-09-04 | 2021-01-22 | 奥诚信息科技(上海)有限公司 | 一种距离测量方法、系统及设备 |
Non-Patent Citations (1)
Title |
---|
光子计数激光雷达游走误差与时间抖动的研究;罗佳承;刘博;华康健;;半导体光电(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
WO2022160611A1 (zh) | 2022-08-04 |
CN112946675A (zh) | 2021-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112946675B (zh) | 一种基于时间融合的距离测量方法、系统和设备 | |
CN112924981B (zh) | 一种飞行时间测距方法、系统和设备 | |
CN108802753B (zh) | 用于确定距对象的距离的设备以及相应的方法 | |
CN110073244B (zh) | 用于确定光子的飞行时间的直方图读出方法和电路 | |
CN109100702B (zh) | 用于测量到对象的距离的光电传感器和方法 | |
CN111679290B (zh) | 光子计数校正方法、激光雷达以及计算机可读介质 | |
US20220196812A1 (en) | Time of flight sensor | |
CN112817001B (zh) | 一种飞行时间测距方法、系统和设备 | |
US11754686B2 (en) | Digital pixel | |
US11668826B2 (en) | Photodetector | |
US20220120872A1 (en) | Methods for dynamically adjusting threshold of sipm receiver and laser radar, and laser radar | |
CN112731425B (zh) | 一种处理直方图的方法、距离测量系统及距离测量设备 | |
EP4016124A1 (en) | Time of flight calculation with inter-bin delta estimation | |
EP4421528A1 (en) | Detection method and apparatus, and terminal | |
US20240248181A1 (en) | Methods and devices for peak signal detection | |
US11885878B2 (en) | Method for gating an output of a single photon avalanche diode (SPAD) cluster based on an ambient count of ambient light events generated by the SPAD cluster in a region of interest of a time of flight (TOF) ranging array | |
JP2020134224A (ja) | 光学的測距装置 | |
JP2024526378A (ja) | レーザーレーダーの測距方法、レーザーレーダー及びコンピュータ可読記憶媒体 | |
US20240159903A1 (en) | Data processing method for lidar and lidar | |
CN115616533A (zh) | 激光雷达的探测方法、发射单元以及激光雷达 | |
US20230273304A1 (en) | Efficient Fault Detection For Lidar Sensors | |
US20230052228A1 (en) | Time of flight sensor | |
CN115184946A (zh) | 基于雷达的距离检测方法、装置、雷达及终端 | |
WO2024131914A1 (en) | Method for dynamically adjusting single-photon detector and idar | |
CN113608230A (zh) | 测距方法、装置及设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |