CN117665741A - 激光雷达回波信号处理方法及装置、激光雷达探测系统 - Google Patents
激光雷达回波信号处理方法及装置、激光雷达探测系统 Download PDFInfo
- Publication number
- CN117665741A CN117665741A CN202211035396.4A CN202211035396A CN117665741A CN 117665741 A CN117665741 A CN 117665741A CN 202211035396 A CN202211035396 A CN 202211035396A CN 117665741 A CN117665741 A CN 117665741A
- Authority
- CN
- China
- Prior art keywords
- pulse
- time
- signal
- intensity
- pulses
- 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.)
- Pending
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 84
- 238000003672 processing method Methods 0.000 title abstract description 11
- 238000000034 method Methods 0.000 claims abstract description 79
- 238000012216 screening Methods 0.000 claims abstract description 77
- 238000012545 processing Methods 0.000 claims abstract description 39
- 238000006243 chemical reaction Methods 0.000 claims description 27
- 238000004364 calculation method Methods 0.000 claims description 21
- 238000001914 filtration Methods 0.000 claims description 19
- 230000003287 optical effect Effects 0.000 claims description 14
- 238000004458 analytical method Methods 0.000 claims description 11
- 238000012952 Resampling Methods 0.000 claims description 6
- 239000000306 component Substances 0.000 description 78
- 230000000875 corresponding effect Effects 0.000 description 42
- 238000010586 diagram Methods 0.000 description 17
- 230000001360 synchronised effect Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000004590 computer program Methods 0.000 description 3
- 238000002592 echocardiography Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000009499 grossing Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 2
- 230000003190 augmentative effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000002310 reflectometry Methods 0.000 description 2
- 229910052710 silicon Inorganic materials 0.000 description 2
- 239000010703 silicon Substances 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 102100034112 Alkyldihydroxyacetonephosphate synthase, peroxisomal Human genes 0.000 description 1
- 241000252073 Anguilliformes Species 0.000 description 1
- 101000799143 Homo sapiens Alkyldihydroxyacetonephosphate synthase, peroxisomal Proteins 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000000848 angular dependent Auger electron spectroscopy Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000035622 drinking Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 238000007634 remodeling Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000717 retained effect 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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
- G01S17/93—Lidar systems specially adapted for specific applications for anti-collision purposes
- G01S17/931—Lidar systems specially adapted for specific applications for anti-collision purposes of land vehicles
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C3/00—Measuring distances in line of sight; Optical rangefinders
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
本发明公开了一种激光雷达回波信号处理方法及装置、激光雷达探测系统,该方法包括:获取接收信号,接收信号包括多个脉冲及其对应的时间信息和强度信息;将接收信号在每个统计单元的强度值分解为直流分量和交流分量,统计单元为时刻点或时钟节拍;根据接收信号的脉冲信息筛选多个脉冲得到候选脉冲,脉冲信息包括直流分量和交流分量;从候选脉冲中选择一个脉冲的到达时间作为TOF值。本发明方案可以提高基于dTOF的回波处理能力,进而提升激光雷达系统的测距性能和精度。
Description
技术领域
本发明涉及信号处理技术领域,具体涉及一种激光雷达回波信号处理方法及装置、激光雷达探测系统。
背景技术
随着高级驾驶辅助系统(Advanced Driving Assistance System,ADAS)到自动驾驶(Automated Driving,AD)技术的不断探索和发展,用于道路及障碍物的人工智能(Artificial Intelligence,AI)识别技术不断提高。对道路及障碍物的感知识别技术包括感知障碍物的精确距离、形状、以及道路类型等感知特征。这些特征感知获取的主流方法之一便是基于直接测量飞行时间(direct Time-Of-Flight,dTOF)原理的激光雷达探测,形成由场景深度信息描绘成的物体形状点云图。
dToF是通过发出短脉冲光然后测量发射的光返回所需的时间来检测与物体的距离,适用于工业机器人等需快速进行测距避障检测的应用、增强现实(Augmented Reality,AR)应用等多种领域。但现有的基于dToF原理的激光雷达探测装置还普遍存在测距能力不足、功耗开销较大的问题。
发明内容
本发明提供一种激光雷达回波信号处理方法及装置、激光雷达探测系统,以提高基于dTOF的回波处理能力,进而提升激光雷达系统的测距性能和精度。
为此,本发明实施例提供如下技术方案:
一方面,本发明实施例提供一种激光雷达回波信号处理方法,所述方法包括:
获取接收信号,所述接收信号包括多个脉冲及其对应的时间信息和强度信息;
将所述接收信号在每个统计单元的强度值分解为直流分量和交流分量,所述统计单元为时刻点或时钟节拍;
根据所述接收信号的脉冲信息筛选所述多个脉冲得到候选脉冲,所述脉冲信息包括所述直流分量和所述交流分量;
从所述候选脉冲中选择一个脉冲的到达时间作为TOF值。
可选地,所述每个统计单元的直流分量为所述接收信号的当前统计单元之前的历史统计单元的强度统计值;所述每个统计单元的交流分量为所述当前统计单元的强度值减去所述直流分量。
可选地,所述方法还包括对所述脉冲进行波形重塑,得到重塑后的接收信号,所述波形重塑包括:判断当前统计单元的信号强度值是否处于波谷,且所述当前统计单元的信号强度与前一脉冲的峰值强度之差是否大于所述当前统计单元的交流分量的关联摆幅,所述交流分量的关联摆幅为所述当前统计单元的交流分量与预设系数的乘积;若是,则将当前统计单元的强度值归为零;若否,则保持当前统计单元的强度值不变。
可选地,所述根据所述接收信号的脉冲信息筛选所述多个脉冲得到候选脉冲包括:
从所述多个脉冲中选择峰值强度最高的设定数量的脉冲为第一筛选脉冲;
计算每个所述第一筛选脉冲的脉冲宽度,并根据所述直流分量和所述交流分量计算每个所述第一筛选脉冲的信噪比;
如果所述第一筛选脉冲的脉冲宽度小于第一阈值,且所述第一筛选脉冲的信噪比小于第二阈值,则删除所述第一筛选脉冲,否则保留所述第一筛选脉冲为候选脉冲。
可选地,所述方法还包括:计算每个所述第一筛选脉冲或每个所述候选脉冲的到达时间。
可选地,所述从所述候选脉冲中选择一个脉冲的到达时间作为TOF值包括:
从所述候选脉冲中选择强度值最大的脉冲的到达时间作为TOF值;或者
从所述候选脉冲中选择到达时间最早并且强度值大于设定值的脉冲的到达时间作为TOF值;或者
从所述候选脉冲中选择到达时间最晚并且强度值大于设定值的脉冲的到达时间作为TOF值。
可选地,所述计算每个所述第一筛选脉冲或每个所述候选脉冲的到达时间包括:
根据所述脉冲的峰值时刻确定所述脉冲的边沿时刻;
根据所述边沿时刻至所述峰值时刻的时间段内脉冲的强度计算到达强度;
根据所述到达强度计算所述脉冲的到达时间。
可选地,所述根据所述脉冲的峰值时刻确定所述脉冲的边沿时刻包括:
确定所述脉冲的峰值时刻和所述脉冲的开始时刻;
确定距离所述峰值时刻前预设时间差的时刻为阈值时刻;
判断所述开始时刻是否在所述阈值时刻至所述峰值时刻之间;
若是,则确定所述开始时刻为所述脉冲的边沿时刻;
若否,则确定所述阈值时刻为所述脉冲的边沿时刻。
可选地,所述根据所述边沿时刻至所述峰值时刻的时间段内脉冲的强度计算到达强度包括:计算所述边沿时刻至所述峰值时刻之间的脉冲强度的中位值或均值为到达强度。
可选地,所述根据所述到达强度计算所述脉冲的到达时间包括:
确定与所述到达强度最近的两个强度为第一强度和第二强度,确定所述第一强度对应的时刻为第一时刻,所述第二强度对应的时刻为第二时刻;
根据所述到达强度、所述第一强度和所述第二强度计算得到时间权重;
根据所述第一时刻、所述第二时刻和所述时间权重计算得到所述脉冲的到达时间。
可选地,所述方法还包括:对所述接收信号进行滤波,以平滑所述接收信号。
可选地,所述方法还包括:对所述TOF值进行校准。
可选地,所述方法还包括:获取关联时间及关联空间内探测单元产生的电信号;将所述电信号转换为数字信号,所述数字信号为所述接收信号。
可选地,所述关联空间包括:回波信号在探测单元阵列中对应的理论区域以及相邻区域,所述相邻区域与所述理论区域不重叠或部分重叠;所述关联时间与所述统计单元对应。
可选地,所述将所述电信号转换为数字信号包括:利用时间数字转换方法将所述电信号转换为数字信号;或者利用模数转换方法将所述电信号转换为数字信号。
可选地,所述利用时间数字转换方法将所述电信号转换为数字信号包括:基于关联时间及关联空间内点亮的所述探测单元的数量,生成一组直方图;根据所述直方图得到所述数字信号。
可选地,所述基于关联时间及关联空间内点亮的所述探测单元的数量,生成一组直方图包括:计算每个所述统计单元中关联空间内点亮的所述探测单元总数,根据每个所述统计单元的所述总数,生成一组直方图;或者统计每个所述统计单元中所述理论区域内点亮的探测单元的数量,以及,统计每个所述统计单元中所述相邻区域内点亮的探测单元的数量,根据所述理论区域内点亮的探测单元的数量和所述相邻区域内点亮的探测单元的数量,生成一组直方图。
可选地,所述利用模数转换方法将所述电信号转换为数字信号包括:获取关联时间及关联空间内探测单元对应的电流信号或电压信号,将所述电流信号或电压信号转换为数字信号。
另一方面,本发明实施例还提供一种激光雷达回波信号处理装置,所述装置包括:
信号获取模块,用于获取接收信号,所述接收信号包括多个脉冲及其对应的时间信息和强度信息;
直交分量计算模块,用于将所述接收信号在每个统计单元的强度值分解为直流分量和交流分量,所述统计单元为时刻点或时钟节拍;
脉冲检测模块,用于根据所述接收信号的脉冲信息筛选所述多个脉冲得到候选脉冲,所述脉冲信息包括所述直流分量和所述交流分量;
TOF分析模块,用于计算各脉冲的到达时间;
TOF决策模块,用于从所述候选脉冲中选择一个脉冲的到达时间作为TOF值。
可选地,所述装置还包括:重采样滤波模块,用于对所述接收信号进行滤波,以平滑所述接收信号。
可选地,所述装置还包括:波形重塑模块,用于对所述脉冲进行波形重塑,得到重塑后的接收信号;所述波形重塑包括:判断当前统计单元的信号强度值是否处于波谷,且所述当前统计单元的信号强度与前一脉冲的峰值强度之差是否大于所述当前统计单元的交流分量的关联摆幅,所述交流分量的关联摆幅为所述当前统计单元的交流分量与预设系数的乘积;若是,则将当前统计单元的强度值归为零;若否,则保持当前统计单元的强度值不变;
可选地,所述脉冲检测模块,具体用于从所述多个脉冲中选择峰值强度最高的设定数量的脉冲为第一筛选脉冲;计算每个所述第一筛选脉冲的脉冲宽度,并根据所述直流分量和所述交流分量计算每个所述第一筛选脉冲的信噪比;如果所述第一筛选脉冲的脉冲宽度小于第一阈值,且所述第一筛选脉冲的信噪比小于第二阈值,则删除所述第一筛选脉冲,否则保留所述第一筛选脉冲为候选脉冲。
可选地,所述TOF分析模块,具体用于计算每个所述第一筛选脉冲,或者计算每个所述候选脉冲的到达时间。
可选地,所述装置还包括:TOF校准模块,用于对所述TOF值进行校准。
可选地,所述信号获取模块包括:信号采集模块,用于获取关联时间及关联空间内探测单元产生的电信号,将所述电信号转换为数字信号,所述数字信号为所述接收信号。
可选地,所述信号采集模块包括TDC模块或者ADC模块;
所述TDC模块,用于利用时间数字转换方法将所述电信号转换为数字信号;
所述ADC模块,用于利用模数转换方法将所述电信号转换为数字信号。
另一方面,本发明实施例还提供一种激光雷达探测系统,用于根据TOF时间确定探测物体的距离,所述系统包括:发射装置、接收装置以及前面所述的激光回波信号处理装置;
所述发射装置,用于发射光脉冲信号;
所述接收装置,包括多个探测单元,用于接收光信号生成接收信号;
所述激光回波信号处理装置,与所述接收装置相连,用于根据接收信号确定发射光脉冲信号的TOF时间,根据所述TOF时间确定探测物体的距离。
本发明提供的激光雷达回波信号处理方法及装置、激光雷达探测系统,获取接收信号,将接收信号中每个统计单元的强度值分解为直流分量和交流分量,根据这些信息对接收信号中的多个脉冲进行筛选,筛选出最有可能的候选脉冲,最后从候选脉冲中选择一个脉冲的到达时间作为TOF值。本发明方案在提高了dTOF回波处理能力的同时,降低了硬件实现的功耗及开销,降低了应用的复杂度,可以极大提升产品的工程化与产品化。
进一步地,在确定每个统计单元的直流分量时,只需采集前向历史值作为参考,即只考虑本统计单元和之前统计单元的数据,而非采集所有统计单元的数据后再进行计算,从而可以逐个统计单元计算统计值,不仅可以极大提高信号处理的实时性,而且还可节省额外的存储开销。本发明方案易于芯片化实现或FPGA实现,且硬化实现逻辑简单,不需要大量的片上RAM作缓存。
进一步地,在对多个脉冲进行筛选时,综合考虑了脉冲的峰值强度、脉冲宽度、信噪比,排除了由于干扰而产生的不符合要求的脉冲,使最终的TOF值更准确。
附图说明
图1是现有技术中激光雷达探测装置的工作原理示意图;
图2是本发明实施例激光雷达回波信号处理方法的一种流程图;
图3是本发明实施例中记录的最强脉冲的示意图;
图4是本发明实施例中近距离回波脉冲的示意图。
图5是本发明实施例中计算脉冲的到达时间的流程图;
图6是本发明实施例中利用时间数字转换方法将电信号转换为数字信号的流程图;
图7是本发明实施例中基于单光子计数机理形成的直方图示意图;
图8是本发明实施例中SPAD阵列中与光阑尺寸形状相匹配的宏像素示意图;
图9是本发明实施例中对波形峰值受噪声影响产生的抖动滤波前后对比示意图;
图10是对图7所示直方图滤波后的直方图示意图;
图11是本发明实施例激光雷达回波信号处理装置的一种结构示意图;
图12是本发明实施例激光雷达回波信号处理装置的另一种结构示意图;
图13是本发明实施例激光雷达回波信号处理装置的另一种结构示意图;
图14是本发明实施例激光雷达回波信号处理装置的另一种结构示意图;
图15是本发明实施例激光雷达探测系统的结构示意图。
具体实施方式
为使本发明的上述目的、特征和有益效果能够更为明显易懂,下面结合附图对本发明的具体实施例做详细的说明。
下面首先结合图1对激光雷达探测装置的工作原理进行简要说明。
参照图1,现有的激光雷达测量装置主要由两部分组成,即发射装置TX和接收装置RX。对应视场中某一个角度(极坐标的水平角与垂直角)位置上的物体10,雷达发射一个很窄的红外光脉冲打到物体10表面,通过物体10反射回来的脉冲回波被探测器接收,例如,探测器为单光子雪崩二极管(single photon avalanche diode,SPAD)阵列,反射的脉冲回波轰击SPAD阵列,SPAD阵列中被点亮的点便可参与计数或者采集电流的变化,即反射回波返回的这段关联时间内形成单光子计数或采集电流变化。假设从激光发射时刻点开始计算,捕捉回波的最大时间片段为2048纳秒,那么最远探测距离即为光速乘上1024纳秒(激光往返时间的一半)。当雷达以一定的扫描方式(逐行、逐列、或其他特定的顺序)对视场中不同角度做完整遍历扫描时,便形成一帧完整的具有一定分辨率的场景深度点云图,该点云图正是后端自动驾驶主控计算平台所需的包含场景深度特征感知的数据。
一种dToF的核心组件包括:激光器、单光子雪崩二极管(Single PhotonAvalanche Diode,SPAD)和时间数字转换器(Time to Digital Convert,TDC)。SPAD是一种具有单光子探测能力的光电探测雪崩二极管,只要有微弱的光信号就能产生电流。其中,激光器可以是垂直腔面发射激光器(vertical cavity surface emitting laser,VCSEL)、边发射激光器(edge emitting laser,EEL),或者其他类型的激光器。dToF模组的激光器向场景中发射脉冲波,SPAD接收从目标物体反射回来的脉冲波。TDC能够记录每次接收到的光信号的飞行时间,也就是发射脉冲和接收脉冲之间的时间间隔。dToF会在一次测量时间内(即测量一个像素点对应的时间)发射和接收N次光信号,每次发射和接收光信号对应的时间为一个时间窗口,然后对记录的N次飞行时间做直方图统计,其中出现频率最高的飞行时间t用来计算待测物体的深度。对于dToF而言,由于雪崩过程存在量子噪声和放大器噪声,以及dToF中TDC模块存在的固有噪声,导致目前dToF的实际精度(真实深度值和相机的测量值之间的差)只能达到cm级。
针对上述问题,本发明实施例提供一种激光雷达回波信号处理方法回波信号处理方法及装置、激光雷达探测系统,以提高基于dTOF的回波处理能力,进而提升激光雷达装置的测距性能和精度。
如图2所示,是本发明实施例激光雷达回波信号处理方法的一种流程图,包括以下步骤:
步骤201,获取接收信号,所述接收信号包括多个脉冲及其对应的时间信息和强度信息。
其中,所述时间信息和强度信息包括:脉冲的开始时间、结束时间、峰值对应的时间点和强度值,还包括每个统计单元的开始时间、结束时间、以及强度值。
需要说明的是,所述接收信号是指数字信号,或者对激光雷达探测单元产生的模拟信号经过模数转换后的数字信号。
步骤202,将所述接收信号在每个统计单元的强度值分解为直流分量和交流分量,所述统计单元为时刻点或时钟节拍。
为了提取出有效回波脉冲,需要进行噪声统计。为此,可以将每个统计单元的强度值分解为直流分量和交流分量。
在本发明实施例的接收信号是指,对应激光雷达视场中某一角度的探测器,其中所述某一角度的探测器可以包括至少一个探测单元(即接收信号是激光雷达视场中对应于一个像素点的探测器产生的信号)。所述统计单元可以以时间点为单位,即统计每个时间点的强度信息;或者所述统计单元以时钟节拍为单位,即统计每个时钟节拍内的强度信息,其中时钟节拍是指一个固定的时间间隔。
例如,当探测器选择SPAD阵列时(即对应于一个像素点的探测器包括多个SPAD),接收信号包括一次测量时间内(发射和接收N次光信号)探测器中点亮的SPAD数量,统计单元为时钟节拍,所述一次测量时间包括多个连续的时钟节拍,每个统计单元计算该时间段内探测器中点亮的SPAD数量作为该统计单元下接收到的光信号的强度值。又例如,当探测器选择SiPM时,探测器输出的接收信号为模拟信号,其中模拟信号的大小作为接收到的光信号的强度值(为便于后续计算需要将模拟信号转换为数字信号),统计单元为时间点或者时间节拍;当以时间点为统计单元时,该时间点模拟信号的大小即为该统计单元的强度值;当以时钟节拍为统计单元时,该时钟节拍对应的时间段内模拟信号大小的统计值(如均值或者中位值或者其他统计方法)即为该统计单元的强度值。
与传统的计算均值、方差时的统计方式的区别在于,在本发明实施例中,只需采集前向历史值作为参考,即只考虑本统计单元和之前统计单元的数据,而非采集所有统计单元的数据后再进行计算。
具体地,对于每个统计单元,将所述接收信号的当前统计单元之前的历史统计单元的强度统计值作为该统计单元的直流分量,将当前统计单元的强度值减去所述直流分量作为该统计单元的交流分量。
统计单元的直流分量和交流分量可以表征为:
DC=Funcdi=[1,n](F2(i)) (1)
AC=Funcai=[1,n](DC(i)) (2)
上式(1)、(2)中DC、AC分别表示直流分量、交流分量当前统计单元的统计值,即直流分量和交流分量。其中,i表示统计单元的序号,F2(i)表示第i个统计单元的强度。Funcd()为自定义函数(例如平均、中位值、方差统计等均可),Funca()为该统计单元的强度值减去直流分量。
需要说明的是,对于初始统计单元,可以认为其直流分量DC为0。
这样,可以逐个统计单元计算统计值,也就是说,在每个统计单元所属时间内,即可获得当前统计单元的统计值,这种方式不仅可以极大提高信号处理的实时性,而且还可节省额外的存储开销。
步骤203,根据所述接收信号的脉冲信息筛选所述多个脉冲得到候选脉冲,所述脉冲信息包括所述直流分量和所述交流分量。
在对所述多个脉冲进行筛选时,可以先选出一定数量峰值强度最高的脉冲,即所述第一筛选脉冲,然后再根据这些脉冲的脉冲宽度和信噪比确定候选脉冲。
在实际应用中,可以记录前N个最强脉冲,并记录其时间信息,若新脉冲的强度大于N个记录中强度最小的脉冲,则用新脉冲替换掉原最小强度的脉冲,以此获得前N个最强脉冲。
比如,记录4个最强脉冲,其对应的时间信息和强度信息如下:
PD0:PD0_TOF=[1780(开始时间),1806(峰值时刻),1826(结束时间)]
PD1:PD1_TOF=[1038(开始时间),1049(峰值时刻),1058(结束时间)]
PD2:PD2_TOF=[1990(开始时间),2007(峰值时刻),2017(结束时间)]
PD3:PD3_TOF=[1914(开始时间),1924(峰值时刻),1935(结束时间)]
PD0_peak=[24.826171875(强度值)]
PD1_peak=[23.866015625(强度值)]
PD2_peak=[36.39296875(强度值)]
PD3_peak=[25.03359375(强度值)]
上述4个最强脉冲对应的波形图如图3所示。
然后,再从记录的4个最强脉冲中选出所述候选脉冲。具体地,计算每个所述第一筛选脉冲的脉冲宽度和信噪比,如果第一筛选脉冲的脉冲宽度小于第一阈值,且所述第一筛选脉冲的信噪比小于第二阈值,则删除所述第一筛选脉冲,否则保留所述第一筛选脉冲为候选脉冲。
脉冲宽度可以根据脉冲对应的时间信息即脉冲的开始时间、结束时间计算得到,脉冲的信噪比可以根据脉冲强度值中的直流分量和交流分量计算得到。
脉冲的信噪比和脉冲宽度可表征如下:
SNR(n)=(peak(n)-dc)/ac (3)
PW(n)=tof_end(n)-tof_start(n) (4)
其中,SNR(n)表示第n个脉冲的信噪比,peak(n)表示该脉冲的强度,dc和ac分别表示该脉冲内包括的统计单元的平均直流分量和平均交流分量;
其中,PW(n)表示第n个脉冲的脉冲宽度,tof_end(n)表示该脉冲的结束时间,tof_start(n)表示该脉冲的开始时间。
在具体应用中,对所述第一筛选脉冲可以首先判断其脉冲宽度,若脉冲宽度大于等于第一阈值,则保留该第一筛选脉冲;如果脉冲的宽度小于第一阈值,则再次对其进行信噪比的判断,如果信噪比小于第二阈值,则删除该第一筛选脉冲(判断为噪声或假性回波),否则保留该第一筛选脉冲。
上述对脉冲宽度和信噪比的判断主要是考虑到近距离(如10米以内)回波拖尾较长且接近饱和以及脉冲宽度较大导致的信噪比过低的情况。由于统计单元的直流分量是当前统计单元之前的所有单元的强度的平均值,因此在近距离的场景下,脉冲的强度很高且宽度很宽,计算得到的直流分量较大,交流分量较小,信噪比比较低。如图4所示,由于脉冲的宽度较宽,每个统计单元的强度基本差异很小,所以计算得到的dc约等于该脉冲的强度值,即dc很大,ac会很小,进而计算得到的信噪比很小。如果仅依据信噪比筛选确定所述候选脉冲,则会漏掉这部分脉冲。通过增加对脉冲宽度的判断,则可避免这种情况,将有可能是近距离的回波对应的脉冲保留下来。
步骤204,从所述候选脉冲中选择一个脉冲的到达时间作为TOF值。
从候选脉冲中选择一个脉冲的到达时间作为TOF值时,可以基于不同的决策原则进行选择,比如可以是但不限于以下几种原则:
(1)最强回波原则:从所述候选脉冲中选择强度值最大的脉冲的到达时间作为TOF值;
(2)第一较强回波原则:从所述候选脉冲中选择到达时间最早并且强度值大于设定值的脉冲的到达时间作为TOF值;
(3)最远距离回波原则:从所述候选脉冲中选择到达时间最晚并且强度值大于设定值的脉冲的到达时间作为TOF值。
因为点云图每个视场角度位置上通常只用一个深度或距离值来表征,因此需要挑选出一个最可能的脉冲的到达时间作为测量得到的TOF值。为此,在上述过程中,还需要计算各脉冲的到达时间。需要说明的是,在具体应用中,可以在得到各候选脉冲后,再计算每个候选脉冲的到达时间;也可以先计算各脉冲的到达时间,再进行筛选,比如对所述多个脉冲中一定数量峰值强度最高的脉冲,即上述第一筛选脉冲,计算所述第一筛选脉冲的到达时间,然后再确定所述候选脉冲。
在本发明实施例中,每个脉冲的到达时间T的计算可与信号检测形成前后级流水运算,无需等到所有脉冲检测完成再启动该计算。即当脉冲0的特征值获取完,便可进行脉冲0的到达时间T的估算,直至第N个最强脉冲分析完。
具体地,如图5所示,示出了本发明实施例中计算每个第一筛选脉冲或每个候选脉冲的到达时间的流程图,包括以下步骤:
在步骤501,根据脉冲的峰值时刻确定所述脉冲的边沿时刻。
在本发明实施例中,可以通过选取脉冲峰值时刻点之前有限时间差内的高频分量进行分析,该有限时间差的大小可以根据需要设定,相当于一个TOF方向的摆动阈值,考虑到计算资源的限制以及计算需要一定的时间,该阈值不宜过大或过小,具体可根据计算资源及设备所需的TOF误差抖动要求来设定,比如可以将该阈值设置为12个统计单元。另外,对于不同的TOF距离,可以设置不同的阈值,以细化适应不同远近距离的TOF误差抖动摆幅。
具体地,确定所述脉冲的峰值时刻和所述脉冲的开始时刻,并确定距离所述峰值时刻前预设时间差的时刻为阈值时刻;判断所述开始时刻是否在所述阈值时刻至所述峰值时刻之间;若是,则确定所述开始时刻为所述脉冲的边沿时刻;若否,则确定所述阈值时刻为所述脉冲的边沿时刻。
在步骤502,根据所述边沿时刻至所述峰值时刻的时间段内脉冲的强度计算到达强度。
也就是说,将脉冲的开始时刻至边沿时刻之间的信号看作噪声,不参与后续到达时间T的计算,将脉冲的边沿时刻至峰值时刻之间的信号看作该脉冲,参与后续到达时间T的计算。
具体地,可以计算所述边沿时刻至所述峰值时刻之间的脉冲强度的中位值或均值或峰值的某一百分比(比如40%)作为到达强度Ft。
在步骤503,根据所述到达强度计算所述脉冲的到达时间。
具体地,可以确定与所述到达强度最近的两个强度为第一强度和第二强度,确定所述第一强度对应的时刻为第一时刻,所述第二强度对应的时刻为第二时刻;根据所述到达强度、所述第一强度和所述第二强度计算得到时间权重;根据所述第一时刻、所述第二时刻和所述时间权重计算得到所述脉冲的到达时间。
比如,寻找与到达强度最近的前后时刻两个强度为F3(i-1)、F3(i),通过计算到达强度Ft与F3(i-1)、F3(i)之间的强度差作为距离差获取权重,并插值得到强度Ft的对应对达时间T。
插值算式可以表征如下:
其中,(i-1)和(i)指该点的时间值,绝对值表示的部分为各时间值的权重。
当然,本发明实施例并不局限于上式(5)的TOF估算插值算子,也可是线性插值等。
本发明提供的激光雷达回波信号处理方法,将接收信号中的每个统计单元的强度值分解为直流分量和交流分量,根据这些信息对接收信号中的多个脉冲进行筛选,筛选出最有可能的候选脉冲,最后从候选脉冲中选择一个脉冲的到达时间作为TOF值。本发明方案在噪声计算、以及有效回波探测方面具有很好的性能表现,通过将每个统计单元的强度值分解成直流分量和交流分量,进而根据每个统计单元的直流分量和交流分量进行波形重塑、信噪比计算以及候选脉冲的筛选等步骤,在低信噪比的情况下也可以筛选出有效回波进行到达时间的计算,如远距离(回波信号弱)低信噪比有效回波探测;并且增加了通过脉宽筛选候选脉冲的条件,可以进一步的在低信噪比的情况下也可以筛选出有效回波进行到达时间的计算,如超近距离(回波信号强)宽饱和有效回波探测。相较于传统的通过固定强度筛选有效脉冲非方法,本申请的方案大大提高了dTOF回波处理能力。另外,与现有技术相比,本发明方案的处理流程处理为单纯的流水线处理结构,易于芯片化实现或FPGA(FieldProgrammable Gate Array,可编程阵列逻辑)实现,且硬化实现逻辑简单,不需要过多的存储单元,降低了硬件实现的功耗及开销,降低了应用的复杂度,可以极大提升产品的工程化与产品化。
在另一种非限制性实施例中,本发明激光雷达回波信号处理方法还可进一步对上述步骤204中得到的TOF值进行校准,比如,可以将对应各TOF值的补偿值写入相应的列表中,通过查表方式来得到校准后的TOF值。具体地,预先设置TOF补偿值对应表;根据所述TOF值查找所述TOF补偿值对应表,得到与所述TOF值最接近的两个TOF值对应的补偿值;根据所述最接近的两个TOF值对应的补偿值对所述TOF值进行校准。
前面提到,在上述步骤201中所述的接收信号是指对激光雷达探测单元产生的电信号经过模数转换后得到的数字信号。基于不同的探测器,所述探测单元也会有所不同,比如可以是:SPAD(Single Photon Avalanche Diode,单光子雪崩二极管)、SiPM(Siliconphotomultiplier,硅光电倍增管)、APD(Avalanche Photo Diode,雪崩光电二极管)等,对此本发明实施例不做限定。
在本发明实施例中,所述电信号具体可以是关联时间及关联空间内探测单元产生的电信号,也就是说,将关联时间及关联空间内探测单元产生的电信号转换后得到的数字信号,所述数字信号即为步骤201中所述的接收信号。
其中,所述关联空间包括:回波信号在探测单元阵列中对应的理论区域以及相邻区域,所述相邻区域与所述理论区域不重叠或部分重叠,而且相邻区域相对于理论区域的位置不做限定,比如可以是上、下、左、右、斜角方向等均可。另外,相邻区域和理论区域可以是任意形状的探测器阵列,比如矩形、圆形、多边形等,对此本发明实施例不做限定。需要说明的是,所述回波信号为激光雷达在某一视场角度下对应的回波信号,所述理论区域为光阑未发生偏移时该视场角度下在探测单元阵列中对应的区域。
其中,所述关联时间与所述统计单元相对应。
在具体应用中,将上述电信号转换为数字信号可以采用多种方式来实现,比如:
(1)利用时间数字转换方法将所述电信号转换为数字信号;
(2)利用模数转换方法将所述电信号转换为数字信号。
下面对上述两种方式转换方式分别进行详细说明。
如图6所示,是本发明实施例中利用时间数字转换方法将电信号转换为数字信号的流程图,包括以下步骤:
步骤601,基于关联时间及关联空间内点亮的探测单元数量,生成一组直方图。
比如,在一种非限制性实现方式中,以探测器采用SPAD接收阵列为例,可以基于关联时间及关联空间内单光子计数机理,用统计单元时间段内激活的单光子计数来表征该统计单元接收的信号强度,若干个这样的统计单元接收的信号强度在时间轴上连续关联起来便形成一组直方图。也就是说,计算每个所述统计单元中关联空间内点亮的所述探测单元总数,根据每个所述统计单元的所述总数,生成一组直方图。如图7所示,为2048个纳秒连续关联时间上形成的直方图示意图。其中,关联空间为理想情况下回波光斑覆盖的SPAD接收阵列的理论区域。考虑到由于设备装调工艺偏差,可能会导致通过光阑的光斑偏离SPAD阵列上正确的接收位置,关联空间还包括与该理论区域相邻的相邻区域的SPAD阵列。计算每个所述统计单元中关联空间内点亮的所述探测单元总数,即计算每个统计单元对饮的时间内,包含理论区域和相邻区域内的关联空间中点亮的探测单元总数,并生成直方图。
考虑到前述光斑偏离SPAD阵列上正确的接收位置,在另一种非限制性实现方式中,还可以采用搜索区域匹配累加的光子计数方法生成所述直方图。具体地,统计每个统计单元中理论区域内点亮的探测单元的数量,以及统计每个统计单元中相邻区域内点亮的探测单元的数量,根据所述理论区域内点亮的探测单元的数量和所述相邻区域内点亮的探测单元的数量,生成一组直方图。即分别统计理论区域和相邻区域内点亮的探测单元总数,分别将统计的数据生成在一组直方图中对应的统计单元中。
在另一种非限制性实现方式中,分别统计每个统计单元中理论区域内点亮的探测单元的数量,以及统计每个统计单元中相邻区域内点亮的探测单元的数量,根据所述理论区域内点亮的探测单元的数量形成一张初始直方图,根据每个所述相邻区域内点亮的探测单元的数量各自形成对应的初始直方图,将各个所述初始直方图分配不同的权重,加权形成理论区域对应视场的直方图,其中理论区域对应的初始直方图权重高于相邻区域的初始直方图权重。例如,如图8所示,假设位于SPAD阵列中心区域的Macro-Pixel(图8中九宫格中心,对应光斑C位置的区域)为理论区域,其环绕其周围的8个Macro-Pixel为相邻区域,理论区域和8个相邻区域各自生成对应初始直方图,将理论区域和8个相邻区域对应初始直方图叠加形成理论区域的直方图,并根据该直方图进行TOF计算,其中理论区域的初始直方图权重为50%,每个相邻区域的初始直方图的权重为6.25%。采用上述方法,可以在光阑发生偏移等情况下保证探测器端能准确获取到回波信号,此外由于通常情况下回波信号落在理论区域的概率更大,因此为理论区域的初始直方图分配更高的权重比,此方法还可以抵消因为扩大探测器的区域范围(由原始的理论区域扩大到理论区域加相邻区域)而造成的环境光噪声增大的问题,可以提高回波信号处理的鲁棒性。
比如图8所示,假设一个宏像素(MP,Macro Pixel)由20*12个SPAD组成,与光阑尺寸形状相匹配。理想情况下回波光斑正好覆盖一个MP,即所述理论区域,但由于接收偏差等因素,光斑覆盖范围可能产生偏移。若假设水平方向搜索最大左右偏移3/4宏像素水平尺寸,垂直方向搜索最大上下偏移1/2宏像素垂直尺寸,即所述相邻区域,则光斑搜索范围为50*24个SPAD组成的阵列范围,即所述关联空间。当然本发明实施例不局限于以上搜索尺寸与精度大小。
需要说明的是,本发明实施例并不局限参与搜索的区域块数,可由工程实现中做定义。
上述搜索区域匹配累加的光子计数方法,以某一空间搜索精度考虑多种偏移情况,可以有效实现关联时间内关联空间的单光子计数。步骤602,根据所述直方图得到所述数字信号。
利用上述第(2)种模数转换方法将所述电信号转换为数字信号时,需要获取关联时间及关联空间内探测单元对应的电流信号或电压信号,将所述电流信号或电压信号转换为数字信号。
需要说明的是,为了解决光斑偏移的问题,同样可以采用上述扩大搜索范围的方法,不同的是,采用模数转换方法将所述电信号转换为数字信号时,需要将关联区域内所有探测单元输出的电信号合并采集,再将采集的信号通过模数转换成数字信号。
进一步地,在本发明激光雷达回波信号处理方法另一非限制性实施例中,还可对上述得到的数字信号进行滤波,以平滑信号、消减噪声,同时改善信号波形数据的关联性。此处的滤波主要是滤除波形峰值受噪声影响抖动的情况。如图9a波形所示,图9a中圆圈部分所示为受噪声影响产生的抖动(该脉冲的峰值处被分裂为两个小峰),上述滤波可以滤除该抖动,形成平滑的脉冲信号。
进一步地,还能够对急剧变化的信号脉冲边沿产生一个小高频过冲,从而有助于区分不同反射率的回波表征情况,尤其对于探测距离较近的脉冲回波,因为近距离脉冲回波信号强度接近饱和,脉冲高度几乎都相同,辨别的难度会增加。如图9b波形所示,图9a中圆圈部分所示为平滑后并增加了小高频过冲的脉冲信号。
因此,通过对数字信号进行滤波,可以进一步提高后续TOF分析计算的准确性。
需要说明的是,本发明实施例对上述滤波的具体运算公式不做限定,只要能起到平滑信号作用即可。
比如,在上述滤波的一种非限制性实现中,滤波运算可以表征如下:
/>
其中,F1表示信号I经过FIR(Finite Impulse Response,有限冲激响应)滤波器滤波后的信号波形,该滤波器卷积系数C1皆为正值,起到低通平滑效果。
其中,F2表示信号F1经过同频半相位重采样插值后的信号,例如,采样第n个点,则采样插值的中心点在第n+0.5的位置,若阶数半径s2设为1,那么参与运算的点有F1(n-1)、F1(n)、F1(n+1)、F1(n+2),其中与中心点相邻的奇数位置FIR系数为正,而偶数位置系数为负,如此随阶数多少交替下去。
为了提高滤波器的带宽与输出效果,可以使用多层重采样滤波。另外,对于上述公式(8)、(9)中的滤波器阶数、系数、层数(公式(8)和公式(9)重复的次数)不做限定。
比如,对于图7所示的直方图,滤波后的直方图如图10所示。可以看出,与滤波前的直方图(图7)相比,滤波后的直方图(图10)更加平滑,滤除了因噪声而导致的不规则抖动,便于后续的处理和计算。
本发明实施例提供的激光雷达回波信号处理方法,在噪声、假性回波滤除、以及有效回波探测方面都有强劲的性能表现。考虑到接收信号包括了经物体反射的回波信号和干扰信号,因此,对接收信号进行处理,将接收信号在每个统计单元的强度值分解为直流分量和交流分量,利用这些信息筛选出多个最有可能的回波进行TOF计算,进而可以得到探测距离。本发明方案在信噪比较低的两种情况:远距离回波信号弱、以及超近距离宽饱和有效回波,仍然可以筛选出有效回波进行TOF计算。而且,本发明方案具有较强的鲁棒性与简易实操性,而且实现功耗及开销、复杂度都很低。
相应地,本发明实施例还提供一种激光雷达回波信号处理装置,如图11所示,是该装置的一种结构示意图。
该实施例的激光雷达回波信号处理装置包括以下各模块:
信号获取模块111,用于获取接收信号,所述接收信号包括多个脉冲及其对应的时间信息和强度信息;
直交分量计算模块112,用于将所述多个脉冲在每个统计单元的强度值分解为直流分量和交流分量,所述统计单元为时刻点或时钟节拍;
脉冲检测模块113,用于根据所述接收信号的脉冲信息筛选所述多个脉冲得到候选脉冲,所述脉冲信息包括所述直流分量和所述交流分量;
TOF分析模块114,用于计算各脉冲的到达时间;
TOF决策模块115,用于从所述候选脉冲中选择一个脉冲的到达时间作为TOF值。
需要说明的是,所述接收信号是指对激光雷达探测单元产生的电信号经过模数转换后的数字信号。
其中,所述脉冲检测模块113具体可以从所述多个脉冲中选择峰值强度最高的设定数量的脉冲为第一筛选脉冲;计算每个所述第一筛选脉冲的脉冲宽度,并根据所述直流分量和所述交流分量计算每个所述第一筛选脉冲的信噪比;如果所述第一筛选脉冲的脉冲宽度小于第一阈值,且所述第一筛选脉冲的信噪比小于第二阈值,则删除所述第一筛选脉冲,否则保留所述第一筛选脉冲为候选脉冲。
在具体应用中,所述TOF分析模块114可以计算每个第一筛选脉冲,或者计算每个候选脉冲的到达时间,具体计算方式可参照前面本发明方法实施例中的描述,在此不再赘述。即TOF分析模块可以在筛选出第一筛选脉冲后对每个第一筛选脉冲的到达时间进行计算,或者在经过第二次筛选得到候选脉冲后再对候选脉冲的到达时间进行计算。
上述TOF决策模块115在从所述候选脉冲中选择一个脉冲的到达时间作为TOF值时,可以基于不同的决策原则进行选择,具体可参照前面本发明方法实施例中的描述,在此不再赘述。
如图12所示,是本发明实施例激光雷达回波信号处理装置的另一种结构示意图。
与图11所示实施例的区别在于,在该实施例中,所述激光雷达回波信号处理装置还包括:重采样滤波模块121,用于对所述接收信号进行滤波,以平滑所述接收信号。
需要说明的是,本发明实施例中对所述重采样滤波模块121具体采用何种滤波算法不做限定,只要能起到平滑信号的作用即可。
如图13所示,是本发明实施例激光雷达回波信号处理装置的另一种结构示意图。
与图12所示实施例相比,激光雷达回波信号处理装置还包括:波形重塑模块131,用于对所述脉冲进行波形重塑,得到重塑后的接收信号。
波形重塑的目的是将每个脉冲的起始与结束点(即脉冲强度)归到零点,以便于后续计算。
为了使脉冲的起始与结束点回归到零点(便于后续运算),传统的做法是减去统一固定直流电平(即信号均值),这种方法有一个缺陷是统计的直流电平不一定很理想,尤其相对于每一个当前时刻点,这样就会导致各时刻点减去的直流电平有些是过大、有些会过小,从而会影响TOF插值计算时的信号脉冲强度的准确性。
为此,在本发明实施例中,不需要将每个时刻点的信号强度减去直流电平(或信号均值)来使信号起始点和结束点趋于归零,而是以前面所述的统计单元为单位,保留每个统计单元的直流分量,从而保留了统计单元的原始特性,只需在TOF分析计算时考虑每个统计单元的交流分量或高频分量即可。
本发明实施例中的波形重塑算法主要有两点:一是通过判断当前统计单元的强度值是否处于波谷位置,二是判断当前统计单元的强度值与前一脉冲的峰值强度之差是否高于当前统计单元的交流分量的关联摆幅;若满足这两个条件,则将该统计单元的强度值归为零点。在具体实现时,所述回波重塑包括:判断当前统计单元的信号强度值是否处于波谷,且所述当前统计单元的信号强度与前一脉冲的峰值强度之差是否大于所述当前统计单元的交流分量的关联摆幅,所述交流分量的关联摆幅为所述当前统计单元的交流分量与预设系数的乘积;若是,则将当前统计单元的强度值归为零;若否,则保持当前统计单元的强度值不变。
需要说明的是,前一脉冲的峰值强度可通过前后统计单元作极值递归比较记录下来。统计单元的信号强度值更新可表征为:
其中,alpha为所述预设系数,F3为重塑后的波形。
如图14所示,是本发明实施例激光雷达回波信号处理装置的另一种结构示意图。
与图11所示实施例的区别在于,在该实施例中,所述激光雷达回波信号处理装置还包括:TOF校准模块141,用于对所述TOF值进行校准。
因为通过TOF决策模块,虽然每个视场角位置获得一个TOF精确估算值,但由于设备等固有因素导致该值与理论值之间往往存在一定的误差,这就需要对每台雷达作一次出厂校准,即对不同距离上的TOF估算值作一次补偿,以使其逼近理论值。
对TOF值的补偿算式可表征如下:
TOF_after_calibration(d)=TOF_before_calibration(d)+delta(d) (11)
其中,delta(d)表示在探测距离d上的校准补偿值。
在工程实现中,可将该距离换算成TOF值并量化成一个离散的TOF表,这样对待校准的TOF值,通过查表寻找与其最接近的两个TOF项所对应的校准补偿值delta作插值即可。当然,本发明实施例并不局限于使用何种插值方法,具体可由工程实现时定义。
具体的,实际应用中误差的大小与被测物体的反射率大小有关,当物体反射率不同时,在相同距离下对应的校准补偿值不同,因此可以通过校准前的距离以及接收信号的强度计算物体反射率,并在校准表中确定在校准前的距离时与计算的反射率最接近的两个反射率对应的补偿值,通过插值法计算出校准补偿值。
需要说明的是,上述TOF校准模块141同样可应用于图12和图13所示实施例。
进一步地,在本发明装置另一实施例中,还可包括:信号采集模块(未图示),用于获取关联时间及关联空间内探测单元产生的电信号,将所述电信号转换为数字信号,所述数字信号即为所述信号获取模块获取的接收信号。
在实际应用中,所述信号采集模块可以采用多种方式将所述电信号转换为数字信号,比如,所述信号采集模块可以包括TDC模块或者ADC模块。
其中:
所述TDC模块,用于利用时间数字转换方法将所述电信号转换为数字信号;
所述ADC模块,用于利用模数转换方法将所述电信号转换为数字信号。
上述不同转换方法的具体实现方式在前面本发明方法实施例中有详细说明,在此不再赘述。
本发明提供的激光雷达回波信号处理装置,针对包括经物体反射的回波信号和干扰信号的接收信号中的多个脉冲,将每个统计单元的强度值分解为直流分量和交流分量,根据这些信息对接收信号中的多个脉冲进行筛选,筛选出最有可能的候选脉冲,最后从候选脉冲中选择一个脉冲的到达时间作为TOF值。本发明方案在保障dTOF回波处理能力的同时,降低了硬件实现的功耗及开销,降低了应用的复杂度,可以极大提升产品的工程化与产品化。
相应地,本发明实施例还提供一种激光雷达探测系统,用于根据所述TOF时间确定探测物体的距离。如图15所示,该激光雷达探测系统150包括:发射装置151、接收装置152以及前面各实施例所述的激光回波信号处理装置153。其中:
所述发射装置151,发射光脉冲信号;
所述接收装置152,包括多个探测单元,用于接收光信号生成接收信号;
所述激光回波信号处理装置153,与所述接收装置152相连,用于根据接收信号确定发射光脉冲信号的TOF时间,根据所述TOF时间确定探测物体的距离。
在具体实施中,关于上述实施例中描述的各个装置、产品包含的各个模块/单元,其可以是软件模块/单元,也可以是硬件模块/单元,或者也可以部分是软件模块/单元,部分是硬件模块/单元。
例如,对于应用于或集成于芯片的各个装置、产品,其包含的各个模块/单元可以都采用电路等硬件的方式实现,或者,至少部分模块/单元可以采用软件程序的方式实现,该软件程序运行于芯片内部集成的处理器,剩余的(如果有)部分模块/单元可以采用电路等硬件方式实现;对于应用于或集成于芯片模组的各个装置、产品,其包含的各个模块/单元可以都采用电路等硬件的方式实现,不同的模块/单元可以位于芯片模组的同一组件(例如芯片、电路模块等)或者不同组件中,或者,至少部分模块/单元可以采用软件程序的方式实现,该软件程序运行于芯片模组内部集成的处理器,剩余的(如果有)部分模块/单元可以采用电路等硬件方式实现;对于应用于或集成于终端的各个装置、产品,其包含的各个模块/单元可以都采用电路等硬件的方式实现,不同的模块/单元可以位于终端内同一组件(例如,芯片、电路模块等)或者不同组件中,或者,至少部分模块/单元可以采用软件程序的方式实现,该软件程序运行于终端内部集成的处理器,剩余的(如果有)部分模块/单元可以采用电路等硬件方式实现。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质为非易失性存储介质或非瞬态存储介质,其上存储有计算机程序,所述计算机程序被处理器运行时执行上述图2或图5或图6对应实施例提供的方法的步骤。
本发明实施例还提供了一种电子设备,包括存储器和处理器,所述存储器上存储有可在所述处理器上运行的计算机程序,所述处理器运行所述计算机程序时执行上述图2或图5或图6对应实施例所提供的方法的步骤。
在本发明实施例中,所述处理器可以为中央处理单元(central processingunit,简称CPU),该处理器还可以是其他通用处理器、数字信号处理器(digital signalprocessor,简称DSP)、专用集成电路(application specific integrated circuit,简称ASIC)、现成可编程门阵列(field programmable gate array,简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
还应理解,本发明实施例中的存储器可以是易失性存储器或非易失性存储器,或可包括易失性和非易失性存储器两者。其中,非易失性存储器可以是只读存储器(read-only memory,简称ROM)、可编程只读存储器(programmable ROM,简称PROM)、可擦除可编程只读存储器(erasable PROM,简称EPROM)、电可擦除可编程只读存储器(electricallyEPROM,简称EEPROM)或闪存。易失性存储器可以是随机存取存储器(random accessmemory,简称RAM),其用作外部高速缓存。通过示例性但不是限制性说明,许多形式的随机存取存储器(random access memory,简称RAM)可用,例如静态随机存取存储器(staticRAM,简称SRAM)、动态随机存取存储器(DRAM)、同步动态随机存取存储器(synchronousDRAM,简称SDRAM)、双倍数据速率同步动态随机存取存储器(double data rate SDRAM,简称DDR SDRAM)、增强型同步动态随机存取存储器(enhanced SDRAM,简称ESDRAM)、同步连接动态随机存取存储器(synchlink DRAM,简称SLDRAM)和直接内存总线随机存取存储器(direct rambus RAM,简称DR RAM)。
应理解,本文中术语“和/或”,仅仅是一种描述关联对象的关联关系,表示可以存在三种关系,例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B这三种情况。另外,本文中字符“/”,表示前后关联对象是一种“或”的关系。
本发明实施例中出现的“多个”是指两个或两个以上。
本发明实施例中出现的第一、第二等描述,仅作示意与区分描述对象之用,没有次序之分,也不表示本发明实施例中对设备个数的特别限定,不能构成对本发明实施例的任何限制。
本发明实施例中出现的“连接”是指直接连接或者间接连接等各种连接方式,以实现设备间的通信,本发明实施例对此不做任何限定。
虽然本发明披露如上,但本发明并非限定于此。任何本领域技术人员,在不脱离本发明的精神和范围内,均可作各种更动与修改,因此本发明的保护范围应当以权利要求所限定的范围为准。
Claims (27)
1.一种激光雷达回波信号处理方法,其特征在于,所述方法包括:
获取接收信号,所述接收信号包括多个脉冲及其对应的时间信息和强度信息;
将所述接收信号在每个统计单元的强度值分解为直流分量和交流分量,所述统计单元为时刻点或时钟节拍;
根据所述接收信号的脉冲信息筛选所述多个脉冲得到候选脉冲,所述脉冲信息包括所述直流分量和所述交流分量;
从所述候选脉冲中选择一个脉冲的到达时间作为TOF值。
2.根据权利要求1所述的方法,其特征在于:
所述每个统计单元的直流分量为所述接收信号的当前统计单元之前的历史统计单元的强度统计值;
所述每个统计单元的交流分量为所述当前统计单元的强度值减去所述直流分量。
3.根据权利要求2所述的方法,其特征在于,所述方法还包括对所述脉冲进行波形重塑,得到重塑后的接收信号,所述波形重塑包括:
判断当前统计单元的信号强度值是否处于波谷,且所述当前统计单元的信号强度与前一脉冲的峰值强度之差是否大于所述当前统计单元的交流分量的关联摆幅,所述交流分量的关联摆幅为所述当前统计单元的交流分量与预设系数的乘积;
若是,则将当前统计单元的强度值归为零;
若否,则保持当前统计单元的强度值不变。
4.根据权利要求2所述的方法,其特征在于,所述根据所述接收信号的脉冲信息筛选所述多个脉冲得到候选脉冲包括:
从所述多个脉冲中选择峰值强度最高的设定数量的脉冲为第一筛选脉冲;
计算每个所述第一筛选脉冲的脉冲宽度,并根据所述直流分量和所述交流分量计算每个所述第一筛选脉冲的信噪比;
如果所述第一筛选脉冲的脉冲宽度小于第一阈值,且所述第一筛选脉冲的信噪比小于第二阈值,则删除所述第一筛选脉冲,否则保留所述第一筛选脉冲为候选脉冲。
5.根据权利要求4所述的方法,其特征在于,所述方法还包括:
计算每个所述第一筛选脉冲或每个所述候选脉冲的到达时间。
6.根据权利要求5所述的方法,其特征在于,所述从所述候选脉冲中选择一个脉冲的到达时间作为TOF值包括:
从所述候选脉冲中选择强度值最大的脉冲的到达时间作为TOF值;或者
从所述候选脉冲中选择到达时间最早并且强度值大于设定值的脉冲的到达时间作为TOF值;或者
从所述候选脉冲中选择到达时间最晚并且强度值大于设定值的脉冲的到达时间作为TOF值。
7.根据权利要求5所述的方法,其特征在于,所述计算每个所述第一筛选脉冲或每个所述候选脉冲的到达时间包括:
根据所述脉冲的峰值时刻确定所述脉冲的边沿时刻;
根据所述边沿时刻至所述峰值时刻的时间段内脉冲的强度计算到达强度;
根据所述到达强度计算所述脉冲的到达时间。
8.根据权利要求7所述的方法,其特征在于,所述根据所述脉冲的峰值时刻确定所述脉冲的边沿时刻包括:
确定所述脉冲的峰值时刻和所述脉冲的开始时刻;
确定距离所述峰值时刻前预设时间差的时刻为阈值时刻;
判断所述开始时刻是否在所述阈值时刻至所述峰值时刻之间;
若是,则确定所述开始时刻为所述脉冲的边沿时刻;
若否,则确定所述阈值时刻为所述脉冲的边沿时刻。
9.根据权利要求7所述的方法,其特征在于,所述根据所述边沿时刻至所述峰值时刻的时间段内脉冲的强度计算到达强度包括:
计算所述边沿时刻至所述峰值时刻之间的脉冲强度的中位值或均值为到达强度。
10.根据权利要求7所述的方法,其特征在于,所述根据所述到达强度计算所述脉冲的到达时间包括:
确定与所述到达强度最近的两个强度为第一强度和第二强度,确定所述第一强度对应的时刻为第一时刻,所述第二强度对应的时刻为第二时刻;
根据所述到达强度、所述第一强度和所述第二强度计算得到时间权重;
根据所述第一时刻、所述第二时刻和所述时间权重计算得到所述脉冲的到达时间。
11.根据权利要求1所述的方法,其特征在于,所述方法还包括:
对所述接收信号进行滤波,以平滑所述接收信号。
12.根据权利要求1所述的方法,其特征在于,所述方法还包括:
对所述TOF值进行校准。
13.根据权利要求1至12任一项所述的方法,其特征在于,所述方法还包括:
获取关联时间及关联空间内探测单元产生的电信号;
将所述电信号转换为数字信号,所述数字信号为所述接收信号。
14.根据权利要求13所述的方法,其特征在于,所述关联空间包括:回波信号在探测单元阵列中对应的理论区域以及相邻区域,所述相邻区域与所述理论区域不重叠或部分重叠;
所述关联时间与所述统计单元对应。
15.根据权利要求14所述的方法,其特征在于,所述将所述电信号转换为数字信号包括:
利用时间数字转换方法将所述电信号转换为数字信号;或者
利用模数转换方法将所述电信号转换为数字信号。
16.根据权利要求15所述的方法,其特征在于,所述利用时间数字转换方法将所述电信号转换为数字信号包括:
基于关联时间及关联空间内点亮的所述探测单元的数量,生成一组直方图;
根据所述直方图得到所述数字信号。
17.根据权利要求16所述的方法,其特征在于,所述基于关联时间及关联空间内点亮的所述探测单元的数量,生成一组直方图包括:
计算每个所述统计单元中关联空间内点亮的所述探测单元总数,根据每个所述统计单元的所述总数,生成一组直方图;或者
统计每个所述统计单元中所述理论区域内点亮的探测单元的数量,以及,统计每个所述统计单元中所述相邻区域内点亮的探测单元的数量,根据所述理论区域内点亮的探测单元的数量和所述相邻区域内点亮的探测单元的数量,生成一组直方图。
18.根据权利要求15所述的方法,其特征在于,所述利用模数转换方法将所述电信号转换为数字信号包括:
获取关联时间及关联空间内探测单元对应的电流信号或电压信号,将所述电流信号或电压信号转换为数字信号。
19.一种激光雷达回波信号处理装置,其特征在于,所述装置包括:
信号获取模块,用于获取接收信号,所述接收信号包括多个脉冲及其对应的时间信息和强度信息;
直交分量计算模块,用于将所述接收信号在每个统计单元的强度值分解为直流分量和交流分量,所述统计单元为时刻点或时钟节拍;
脉冲检测模块,用于根据所述接收信号的脉冲信息筛选所述多个脉冲得到候选脉冲,所述脉冲信息包括所述直流分量和所述交流分量;
TOF分析模块,用于计算各脉冲的到达时间;
TOF决策模块,用于从所述候选脉冲中选择一个脉冲的到达时间作为TOF值。
20.根据权利要求19所述的装置,其特征在于,所述装置还包括:
重采样滤波模块,用于对所述接收信号进行滤波,以平滑所述接收信号。
21.根据权利要求19所述的装置,其特征在于,所述装置还包括:
波形重塑模块,用于对所述脉冲进行波形重塑,得到重塑后的接收信号;
所述波形重塑包括:判断当前统计单元的信号强度值是否处于波谷,且所述当前统计单元的信号强度与前一脉冲的峰值强度之差是否大于所述当前统计单元的交流分量的关联摆幅,所述交流分量的关联摆幅为所述当前统计单元的交流分量与预设系数的乘积;若是,则将当前统计单元的强度值归为零;若否,则保持当前统计单元的强度值不变。
22.根据权利要求19所述的装置,其特征在于:
所述脉冲检测模块,具体用于从所述多个脉冲中选择峰值强度最高的设定数量的脉冲为第一筛选脉冲;计算每个所述第一筛选脉冲的脉冲宽度,并根据所述直流分量和所述交流分量计算每个所述第一筛选脉冲的信噪比;如果所述第一筛选脉冲的脉冲宽度是否小于第一阈值,且所述第一筛选脉冲的信噪比小于第二阈值,则删除所述第一筛选脉冲,否则保留所述第一筛选脉冲为候选脉冲。
23.根据权利要求22所述的装置,其特征在于:
所述TOF分析模块,具体用于计算每个所述第一筛选脉冲,或者计算每个所述候选脉冲的到达时间。
24.根据权利要求19所述的装置,其特征在于,所述装置还包括:
TOF校准模块,用于对所述TOF值进行校准。
25.根据权利要求19至24任一项所述的装置,其特征在于,所述信号获取模块包括:
信号采集模块,用于获取关联时间及关联空间内探测单元产生的电信号,将所述电信号转换为数字信号,所述数字信号为所述接收信号。
26.根据权利要求25所述的装置,其特征在于,所述信号采集模块包括TDC模块或者ADC模块;
所述TDC模块,用于利用时间数字转换方法将所述电信号转换为数字信号;
所述ADC模块,用于利用模数转换方法将所述电信号转换为数字信号。
27.一种激光雷达探测系统,用于根据TOF时间确定探测物体的距离,其特征在于,所述系统包括:发射装置、接收装置以及如权利要求19至26任一项所述的激光回波信号处理装置;
所述发射装置,用于发射光脉冲信号;
所述接收装置,包括多个探测单元,用于接收光信号生成接收信号;
所述激光回波信号处理装置,与所述接收装置相连,用于根据接收信号确定发射光脉冲信号的TOF时间,根据所述TOF时间确定探测物体的距离。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211035396.4A CN117665741A (zh) | 2022-08-26 | 2022-08-26 | 激光雷达回波信号处理方法及装置、激光雷达探测系统 |
PCT/CN2023/079470 WO2024040912A1 (zh) | 2022-08-26 | 2023-03-03 | 激光雷达回波信号处理方法及装置、激光雷达探测系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211035396.4A CN117665741A (zh) | 2022-08-26 | 2022-08-26 | 激光雷达回波信号处理方法及装置、激光雷达探测系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117665741A true CN117665741A (zh) | 2024-03-08 |
Family
ID=90012268
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211035396.4A Pending CN117665741A (zh) | 2022-08-26 | 2022-08-26 | 激光雷达回波信号处理方法及装置、激光雷达探测系统 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN117665741A (zh) |
WO (1) | WO2024040912A1 (zh) |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4129547B2 (ja) * | 2004-02-13 | 2008-08-06 | 三菱電機株式会社 | 震動変位計算精度の向上手段 |
JP5554688B2 (ja) * | 2010-11-19 | 2014-07-23 | 株式会社デンソー | レーダ装置 |
JP6497911B2 (ja) * | 2014-11-28 | 2019-04-10 | キヤノン株式会社 | 信号処理方法、信号処理装置、音響波処理方法および音響波処理装置 |
EP3370079B1 (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 - pulse detection |
WO2018181013A1 (ja) * | 2017-03-29 | 2018-10-04 | 株式会社デンソー | 光検出器 |
EP3460509A1 (en) * | 2017-09-22 | 2019-03-27 | ams AG | Method for calibrating a time-of-flight system and time-of-flight system |
US10830880B2 (en) * | 2018-03-20 | 2020-11-10 | Panosense Inc. | Selecting LIDAR pulse detector depending on pulse type |
CN114402225A (zh) * | 2020-06-03 | 2022-04-26 | 深圳市大疆创新科技有限公司 | 测距方法、测距装置和可移动平台 |
CN111679290B (zh) * | 2020-06-04 | 2023-05-05 | 上海禾赛科技有限公司 | 光子计数校正方法、激光雷达以及计算机可读介质 |
CN112924981B (zh) * | 2021-01-28 | 2023-10-31 | 深圳奥锐达科技有限公司 | 一种飞行时间测距方法、系统和设备 |
CN112817001B (zh) * | 2021-01-28 | 2023-12-01 | 深圳奥锐达科技有限公司 | 一种飞行时间测距方法、系统和设备 |
CN114488175A (zh) * | 2022-01-21 | 2022-05-13 | 深圳市灵明光子科技有限公司 | 一种提高信噪比的直方图构造方法及激光测距芯片 |
-
2022
- 2022-08-26 CN CN202211035396.4A patent/CN117665741A/zh active Pending
-
2023
- 2023-03-03 WO PCT/CN2023/079470 patent/WO2024040912A1/zh unknown
Also Published As
Publication number | Publication date |
---|---|
WO2024040912A1 (zh) | 2024-02-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20240028667A1 (en) | Method and apparatus for processing a histogram output from a detector sensor | |
US20210181317A1 (en) | Time-of-flight-based distance measurement system and method | |
US20210116572A1 (en) | Light ranging apparatus | |
US20220196812A1 (en) | Time of flight sensor | |
KR20140057625A (ko) | 이동거리시간차 신호들을 프로세싱하는데 있어서 또는 프로세싱과 관련된 개선들 | |
WO2022160611A1 (zh) | 一种基于时间融合的距离测量方法、系统和设备 | |
CN110361751B (zh) | 时间飞行深度相机及单频调制解调的降低噪声的距离测量方法 | |
CN112817001B (zh) | 一种飞行时间测距方法、系统和设备 | |
JP6261681B2 (ja) | タイムオブフライト信号の処理における又はこれに関する改良 | |
JP2004538491A (ja) | 三次元距離画像を記録するための方法及び装置 | |
CN113466836A (zh) | 一种测距方法、装置及激光雷达 | |
JP2015219120A (ja) | 距離測定装置 | |
WO2020237448A1 (zh) | 回波信号处理方法、装置、系统及存储介质 | |
US20220187430A1 (en) | Time of flight calculation with inter-bin delta estimation | |
CN112470026A (zh) | 激光雷达及其探测方法、存储介质和探测系统 | |
WO2022206031A1 (zh) | 确定噪声水平的方法、激光雷达以及测距方法 | |
CN113424077A (zh) | 光学测距装置 | |
CN115453497A (zh) | 一种低信噪比稀疏雷达信号的滤波方法及系统 | |
EP3721261B1 (en) | Distance time-of-flight modules | |
CN117665741A (zh) | 激光雷达回波信号处理方法及装置、激光雷达探测系统 | |
US20210396876A1 (en) | Optical distance measurement apparatus | |
US20210389462A1 (en) | System and method for histogram binning for depth detection | |
KR20190092957A (ko) | 거리 측정 방법 | |
US20240036208A1 (en) | Time of flight sensing | |
Kufner et al. | A Digital Signal Processing System for Terrestrial Laser Range Finding Systems |
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 |