CN113534179B - 一种激光雷达测距方法及装置 - Google Patents

一种激光雷达测距方法及装置 Download PDF

Info

Publication number
CN113534179B
CN113534179B CN202110678058.1A CN202110678058A CN113534179B CN 113534179 B CN113534179 B CN 113534179B CN 202110678058 A CN202110678058 A CN 202110678058A CN 113534179 B CN113534179 B CN 113534179B
Authority
CN
China
Prior art keywords
signal
sequence
echo
points
linear interpolation
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
Application number
CN202110678058.1A
Other languages
English (en)
Other versions
CN113534179A (zh
Inventor
程坤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CETHIK Group Ltd
Original Assignee
CETHIK Group Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CETHIK Group Ltd filed Critical CETHIK Group Ltd
Priority to CN202110678058.1A priority Critical patent/CN113534179B/zh
Publication of CN113534179A publication Critical patent/CN113534179A/zh
Application granted granted Critical
Publication of CN113534179B publication Critical patent/CN113534179B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/89—Lidar systems specially adapted for specific applications for mapping or imaging
    • G01S17/894—3D imaging with simultaneous measurement of time-of-flight at a 2D array of receiver pixels, e.g. time-of-flight cameras or flash lidar
    • 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/483—Details of pulse systems
    • G01S7/486—Receivers
    • G01S7/487—Extracting wanted echo signals, e.g. pulse detection
    • 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/491—Details of non-pulse systems
    • G01S7/493—Extracting wanted echo signals
    • 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
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Abstract

本发明公开了一种激光雷达测距方法及装置，方法包括：获取调制信号以及回波信号，回波信号为待测对象对发射光脉冲进行反射得到，发射光脉冲为根据调制信号调制得到；对调制信号进行线性插值后做傅里叶变换，得到调制信号的频域信号；对回波信号进行线性插值后做傅里叶变换，得到回波信号的频域信号；根据调制信号的频域信号和回波信号的频域信号，基于傅里叶逆变换得到待测对象与激光雷达之间的距离。本发明通过线性插值方法转换到频率域计算，降低了ADC采样频率的要求。

Description

一种激光雷达测距方法及装置
技术领域
本申请属于激光测距技术领域，具体涉及一种激光雷达测距方法及装置。
背景技术
全固态FLASH激光雷达无机械旋转部件，可以进行像素级实时三维成像，为ADAS(汽车高级驾驶辅助系统)提供超高的距离分辨和空间分辨能力，因此FLASH激光雷达越来越受到市场的青睐。其中测距范围、空间分辨率、测距精度是FLASH激光雷达最主要的性能指标。
FLASH激光雷达系统主要包括发射模块、接收模块、控制模块和计算模块。在接收模块端，目标物的反射回波经过光电转换、前置信号放大、滤波处理和模数转换后进入信号处理，解调出目标物的距离、目标物上的反射率、角度等信息。但现有技术中为了得到较高精度的测距值，需要对激光雷达数据进行高频采集，这不仅增加了数据采集硬件的成本，还容易因数据采集丢帧而影响测距精度。
发明内容
本申请的目的在于提供一种激光雷达测距方法及装置，通过线性插值方法转换到频率域计算，降低了ADC采样频率的要求。
为实现上述目的，本申请所采取的技术方案为：
以下还提供了若干可选方式，但并不作为对上述总体方案的额外限定，仅仅是进一步的增补或优选，在没有技术或逻辑矛盾的前提下，各可选方式可单独针对上述总体方案进行组合，还可以是多个可选方式之间进行组合。
一种激光雷达测距方法，所述激光雷达测距方法，包括：
步骤S1、获取调制信号以及回波信号，所述回波信号为待测对象对发射光脉冲进行反射得到，所述发射光脉冲为根据所述调制信号调制得到；
步骤S2、对所述调制信号进行线性插值后做傅里叶变换，得到调制信号的频域信号；
步骤S3、对所述回波信号进行线性插值后做傅里叶变换，得到回波信号的频域信号；
步骤S4、根据调制信号的频域信号和回波信号的频域信号，基于傅里叶逆变换得到待测对象与激光雷达之间的距离，包括：
步骤S41、取调制信号的频域信号对应的序列D(m)和回波信号的频域信号对应的序列X(m)，将序列D(m)和X(m)相乘得到序列Y(m)，m＝1,2,3,…,N+k(N-1)，N为调制信号中信号点的总个数，k为线性插值时在每相邻两个信号点之间插入点的个数；
步骤S42、对序列Y(m)进行傅里叶逆变换得到序列y(m)；
步骤S43、取序列y(m)中实部最大的元素y(i)，1≤i≤N+k(N-1)，则该元素y(i)对应的时间刻度为i/(N+k(N-1))；
步骤S44、取调制信号的周期为T，则回波信号在激光雷达与待测对象之间的往返时间t为t＝T×i/(N+k(N-1))+Δt，Δt为发射光脉冲发射前对应的随机延迟时间；
步骤S45、得到待测对象与激光雷达之间的距离d为式中c为光速。
作为优选，所述调制信号为随机幅度调制信号，所述发射光脉冲为根据所述调制信号调制得到，包括：
生成随机延迟时间和随机幅度调制信号；
根据所述随机延迟时间延迟间隔时间后，基于所述随机幅度调制信号产生对应的发射光脉冲，其中随机幅度调制信号的周期为T，随机延迟时间的间隔时间为0～1/2T。
作为优选，所述对所述调制信号进行线性插值后做傅里叶变换，得到调制信号的频域信号，包括：
令调制信号为b(n)，其中n表示调制信号中信号点的时刻序号，且n＝1,2,3,…,N，N为调制信号中信号点的总个数；
取线性插值的个数k，在调制信号b(n)的每相邻两个信号点之间进行线性插值生成新的序列c(m)，m＝1,2,3,…,N+k(N-1)，即在b(q)和b(q+1)之间均匀插入k个点，得到b(q),b(q,1),b(q,2),…,b(q,k),b(q+1)，q＝1,2,3,…,N-1；
对序列c(m)进行傅里叶变换生成序列C(m)：
取序列C(m)的共轭序列得到序列D(m)作为调制信号的频域信号。
作为优选，所述对所述回波信号进行线性插值后进行频谱变换，得到回波信号的频谱信号，包括：
令回波信号为r(n)，其中n表示回波信号中信号点的时刻序号，即为调制信号中信号点的时刻序号，且n＝1,2,3,…,N，N为回波信号中信号点的总个数，即为调制信号中信号点的总个数；
取线性插值的个数k，在回波信号r(n)的每相邻两个信号点之间进行线性插值生成新的序列x(m)，m＝1,2,3,…,N+k(N-1)，即在r(q)和r(q+1)之间均匀插入k个点，得到r(q),r(q,1),r(q,2),…,r(q,k),r(q+1)，q＝1,2,3,…,N-1；
对序列x(m)进行傅里叶变换得到序列X(m)如下：
将序列X(m)作为回波信号的频谱信号。
一种激光雷达测距装置，所述激光雷达测距装置包括控制模块、发射模块、接收模块和计算模块，其中：
所述控制模块，用于生成调制信号；
所述发射模块，用于根据所述调制信号调制得到发射光脉冲并发出；
所述接收模块，用于接收回波信号，所述回波信号为待测对象对发射光脉冲进行反射得到；
所述计算模块，用于获取调制信号以及回波信号，对所述调制信号进行线性插值后做傅里叶变换，得到调制信号的频域信号，对所述回波信号进行线性插值后做傅里叶变换，得到回波信号的频域信号，根据调制信号的频域信号和回波信号的频域信号，基于傅里叶逆变换得到待测对象与激光雷达之间的距离；
所述计算模块基于傅里叶逆变换得到待测对象与激光雷达之间的距离，执行如下操作：
取调制信号的频域信号对应的序列D(m)和回波信号的频域信号对应的序列X(m)，将序列D(m)和X(m)相乘得到序列Y(m)，m＝1,2,3,…,N+k(N-1)，N为调制信号中信号点的总个数，k为线性插值时在每相邻两个信号点之间插入点的个数；
对序列Y(m)进行傅里叶逆变换得到序列y(m)；
取序列y(m)中实部最大的元素y(i)，1≤i≤N+k(N-1)，则该元素y(i)对应的时间刻度为i/(N+k(N-1))；
取调制信号的周期为T，则回波信号在激光雷达与待测对象之间的往返时间t为t＝T×i/(N+k(N-1))+Δt，Δt为发射光脉冲发射前对应的随机延迟时间；
得到待测对象与激光雷达之间的距离d为式中c为光速。
作为优选，所述调制信号为随机幅度调制信号，所述发射模块根据所述调制信号调制得到发射光脉冲，执行如下操作：
生成随机延迟时间和随机幅度调制信号；
根据所述随机延迟时间延迟间隔时间后，基于所述随机幅度调制信号产生对应的发射光脉冲，其中随机幅度调制信号的周期为T，随机延迟时间的间隔时间为0～1/2T。
作为优选，所述计算模块对所述调制信号进行线性插值后做傅里叶变换，得到调制信号的频域信号，执行如下操作：
令调制信号为b(n)，其中n表示调制信号中信号点的时刻序号，且n＝1,2,3,…,N，N为调制信号中信号点的总个数；
取线性插值的个数k，在调制信号b(n)的每相邻两个信号点之间进行线性插值生成新的序列c(m)，m＝1,2,3,…,N+k(N-1)，即在b(q)和b(q+1)之间均匀插入k个点，得到b(q),b(q,1),b(q,2),…,b(q,k),b(q+1)，q＝1,2,3,…,N-1；
对序列c(m)进行傅里叶变换生成序列C(m)：
取序列C(m)的共轭序列得到序列D(m)作为调制信号的频域信号。
作为优选，所述计算模块对所述回波信号进行线性插值后进行频谱变换，得到回波信号的频谱信号，执行如下操作：
令回波信号为r(n)，其中n表示回波信号中信号点的时刻序号，即为调制信号中信号点的时刻序号，且n＝1,2,3,…,N，N为回波信号中信号点的总个数，即为调制信号中信号点的总个数；
取线性插值的个数k，在回波信号r(n)的每相邻两个信号点之间进行线性插值生成新的序列x(m)，m＝1,2,3,…,N+k(N-1)，即在r(q)和r(q+1)之间均匀插入k个点，得到r(q),r(q,1),r(q,2),…,r(q,k),r(q+1)，q＝1,2,3,…,N-1；
对序列x(m)进行傅里叶变换得到序列X(m)如下：
将序列X(m)作为回波信号的频谱信号。
本申请提供的激光雷达测距方法及装置，对调制信号和回波信号进行线性插值，降低对回波信号采样频率的要求，提高激光雷达测距精度。
附图说明
图1为本申请的激光雷达测距方法的流程图；
图2为本申请的激光雷达测距装置的结构示意图；
图3为本申请发射模块和接收模块的一个实施例结构示意图。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
除非另有定义，本文所使用的所有的技术和科学术语与属于本申请的技术领域的技术人员通常理解的含义相同。本文中在本申请的说明书中所使用的术语只是为了描述具体的实施例的目的，不是在于限制本申请。
其中一个实施例中，提供一种激光雷达测距方法，该方法适用于任意场景下的激光测距，例如自动驾驶系统、测绘系统、机器视觉系统等等，本实施例提供的测距方法可有效降低数据采集频率，因此可有效克服现有技术中测距范围和成像距离过近的问题，实现在远、近距离情况下均能够具备高精度测距的目的。
具体的，如图1所示，本实施例中的激光雷达测距方法，包括：
步骤1、获取调制信号以及回波信号，其中，回波信号为待测对象对发射光脉冲进行反射得到，发射光脉冲为根据所述调制信号调制得到。
对于本申请而言，待测对象可以是移动对象，也可以是固定对象，并且对象可以是人、动物、物体等，本实施例中不进行限制。
其中调制信号可以是固定幅度的周期信号，也可以是采用随机幅度的周期信号。为了提高本申请测距过程中的抗干扰能力，在一个实施例中，调制信号为随机幅度调制信号，同时根据调制信号调制得到发射光脉冲，具体包括以下步骤：
生成随机延迟时间和随机幅度调制信号；根据所述随机延迟时间延迟间隔时间后，基于所述随机幅度调制信号产生对应的发射光脉冲，其中随机幅度调制信号的周期为T，随机延迟时间的间隔时间为0～1/2T。
其中，随机延迟时间和随机幅度调制信号可以由随机数发生器产生，将随机数发生器针对随机幅度调制信号产生的数据经过数字模拟变换器转化为对应的随机幅度调制信号。
步骤2、对调制信号进行线性插值后做傅里叶变换，得到调制信号的频域信号。
具体为，令调制信号为b(n)，其中n表示调制信号中信号点的时刻序号，且n＝1,2,3,…,N，N为调制信号中信号点的总个数。
取线性插值的个数k，在调制信号b(n)的每相邻两个信号点之间进行线性插值生成新的序列c(m)，m＝1,2,3,…,N+k(N-1)，即在b(q)和b(q+1)之间均匀插入k个点，得到b(q),b(q,1),b(q,2),…,b(q,k),b(q+1)，q＝1,2,3,…,N-1。
对序列c(m)进行傅里叶变换生成序列C(m)：
取序列C(m)的共轭序列得到序列D(m)作为调制信号的频域信号。
本实施例对调制信号先进行线性插值以扩充信号点，并在线性插值后进行傅里叶变换，建立离散的信号点之间的关联，使得处理后的调制信号仍具备相应的物理意义。
步骤3、对回波信号进行线性插值后做傅里叶变换，得到回波信号的频域信号。
具体为，令回波信号为r(n)，其中n表示回波信号中信号点的时刻序号，即为调制信号中信号点的时刻序号，且n＝1,2,3,…,N，N为回波信号中信号点的总个数，即为调制信号中信号点的总个数。
取线性插值的个数k，在回波信号r(n)的每相邻两个信号点之间进行线性插值生成新的序列x(m)，m＝1,2,3,…,N+k(N-1)，即在r(q)和r(q+1)之间均匀插入k个点，得到r(q),r(q,1),r(q,2),…,r(q,k),r(q+1)，q＝1,2,3,…,N-1。
对序列x(m)进行傅里叶变换得到序列X(m)如下：
将序列X(m)作为回波信号的频谱信号。
本实施例同时对调制信号和回波信号进行线性插值以及傅里叶变换，得到具有可比性的两组数据，以便于后续进行测距计算，避免仅基于回波信号进行计算，既提高了测距计算精度，又降低了测距计算复杂度。
需要说明的是，本实施例的激光雷达测距方法中对步骤2和步骤3的执行顺序没有严格限制，即可以先执行步骤2再执行步骤3，或先执行步骤3再执行步骤3，或者同步执行步骤2和步骤3。
步骤4、根据调制信号的频域信号和回波信号的频域信号，基于傅里叶逆变换得到待测对象与激光雷达之间的距离。
步骤S41、取调制信号的频域信号对应的序列D(m)和回波信号的频域信号对应的序列X(m)，将序列D(m)和X(m)相乘得到序列Y(m)，m＝1,2,3,…,N+k(N-1)，N为调制信号中信号点的总个数，k为线性插值时在每相邻两个信号点之间插入点的个数。
步骤S42、对序列Y(m)进行傅里叶逆变换得到序列y(m)。
步骤S43、取序列y(m)中实部最大的元素y(i)，1≤i≤N+k(N-1)，则该元素y(i)对应的时间刻度为i/(N+k(N-1))。
步骤S44、取调制信号的周期为T，则回波信号在激光雷达与待测对象之间的往返时间t为t＝T×i/(N+k(N-1))+Δt，Δt为发射光脉冲发射前对应的随机延迟时间。
步骤S45、得到待测对象与激光雷达之间的距离d为式中c为光速。
如图2所示，在另一个实施例中，还提供了一种激光雷达测距装置，激光雷达测距装置包括控制模块、发射模块、接收模块和计算模块，其中：
控制模块，用于生成调制信号；
发射模块，用于根据所述调制信号调制得到发射光脉冲并发出；
接收模块，用于接收回波信号，所述回波信号为待测对象对发射光脉冲进行反射得到；
计算模块，用于获取调制信号以及回波信号，对所述调制信号进行线性插值后做傅里叶变换，得到调制信号的频域信号，对所述回波信号进行线性插值后做傅里叶变换，得到回波信号的频域信号，根据调制信号的频域信号和回波信号的频域信号，基于傅里叶逆变换得到待测对象与激光雷达之间的距离。
其中，计算模块基于傅里叶逆变换得到待测对象与激光雷达之间的距离，执行如下操作：
取调制信号的频域信号对应的序列D(m)和回波信号的频域信号对应的序列X(m)，将序列D(m)和X(m)相乘得到序列Y(m)，m＝1,2,3,…,N+k(N-1)，N为调制信号中信号点的总个数，k为线性插值时在每相邻两个信号点之间插入点的个数；
对序列Y(m)进行傅里叶逆变换得到序列y(m)；
取序列y(m)中实部最大的元素y(i)，1≤i≤N+k(N-1)，则该元素y(i)对应的时间刻度为i/(N+k(N-1))；
取调制信号的周期为T，则回波信号在激光雷达与待测对象之间的往返时间t为t＝T×i/(N+k(N-1))+Δt，Δt为发射光脉冲发射前对应的随机延迟时间；
得到待测对象与激光雷达之间的距离d为式中c为光速。
为了提高本申请数据的有效性，如图3所示，本申请的发射模块包括随机数发生器、数字模拟变换器和激光器，接收模块包括光电传感器、前置放大滤波电路和模拟数据变换器。
在发射模块工作时，由随机数发生器生成针对随机延迟时间和随机幅度调制信号的随机数，并且由数字模拟变换器将针对随机幅度调制信号的随机数进行数模转换得到随机幅度调制信号，最后由激光器接收随机幅度调制信号生成发射光脉冲发出。
在接收模块工作时，由光电传感器接收由被测物体返回的回波信号，并将回波信号经过前置放大滤波电路，对回波信号进行滤波以提高信噪比，后将滤波后的回波信号输入模拟数据变换器，经过模数转换后变为数字信号以便于后续处理。并且为了控制数据的一致性，本实施例设置数字模拟变换器和模拟数字变换器的调制频率一致。
在另一个实施例中，调制信号为随机幅度调制信号，发射模块根据所述调制信号调制得到发射光脉冲，执行如下操作：
生成随机延迟时间和随机幅度调制信号；
根据所述随机延迟时间延迟间隔时间后，基于所述随机幅度调制信号产生对应的发射光脉冲，其中随机幅度调制信号的周期为T，随机延迟时间的间隔时间为0～1/2T。
在另一个实施例中，计算模块对所述调制信号进行线性插值后做傅里叶变换，得到调制信号的频域信号，执行如下操作：
令调制信号为b(n)，其中n表示调制信号中信号点的时刻序号，且n＝1,2,3,…,N，N为调制信号中信号点的总个数；
取线性插值的个数k，在调制信号b(n)的每相邻两个信号点之间进行线性插值生成新的序列c(m)，m＝1,2,3,…,N+k(N-1)，即在b(q)和b(q+1)之间均匀插入k个点，得到b(q),b(q,1),b(q,2),…,b(q,j),b(q+1)，q＝1,2,3,…,N-1；
对序列c(m)进行傅里叶变换生成序列C(m)：
取序列C(m)的共轭序列得到序列D(m)作为调制信号的频域信号。
在另一个实施例中，计算模块对所述回波信号进行线性插值后进行频谱变换，得到回波信号的频谱信号，执行如下操作：
令回波信号为r(n)，其中n表示回波信号中信号点的时刻序号，即为调制信号中信号点的时刻序号，且n＝1,2,3,…,N，N为回波信号中信号点的总个数，即为调制信号中信号点的总个数；
取线性插值的个数k，在回波信号r(n)的每相邻两个信号点之间进行线性插值生成新的序列x(m)，m＝1,2,3,…,N+k(N-1)，即在r(q)和r(q+1)之间均匀插入k个点，得到rr(q),r(q,1),r(q,2),…,r(q,k),r(q+1)，q＝1,2,3,…,N-1；
对序列x(m)进行傅里叶变换得到序列X(m)如下：
将序列X(m)作为回波信号的频谱信号。
关于一种激光雷达测距装置的具体限定，可参见上述对于一种激光雷达测距方法的限定，这里不再进行赘述。
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。

Claims (8)

1.一种激光雷达测距方法，其特征在于，所述激光雷达测距方法，包括：
步骤S1、获取调制信号以及回波信号，所述回波信号为待测对象对发射光脉冲进行反射得到，所述发射光脉冲为根据所述调制信号调制得到；
步骤S2、对所述调制信号进行线性插值后做傅里叶变换，得到调制信号的频域信号；
步骤S3、对所述回波信号进行线性插值后做傅里叶变换，得到回波信号的频域信号；
步骤S4、根据调制信号的频域信号和回波信号的频域信号，基于傅里叶逆变换得到待测对象与激光雷达之间的距离，包括：
步骤S41、取调制信号的频域信号对应的序列D(m)和回波信号的频域信号对应的序列X(m)，将序列D(m)和X(m)相乘得到序列Y(m)，m＝1,2,3,…,N+k(N-1)，N为调制信号中信号点的总个数，k为线性插值时在每相邻两个信号点之间插入点的个数；
步骤S42、对序列Y(m)进行傅里叶逆变换得到序列y(m)；
步骤S43、取序列y(m)中实部最大的元素y(i)，1≤i≤N+k(N-1)，则该元素y(i)对应的时间刻度为i/(N+k(N-1))；
步骤S44、取调制信号的周期为T，则回波信号在激光雷达与待测对象之间的往返时间t为t＝T×i/(N+k(N-1))+Δt，Δt为发射光脉冲发射前对应的随机延迟时间；
步骤S45、得到待测对象与激光雷达之间的距离d为式中c为光速。
2.如权利要求1所述的激光雷达测距方法，其特征在于，所述调制信号为随机幅度调制信号，所述发射光脉冲为根据所述调制信号调制得到，包括：
生成随机延迟时间和随机幅度调制信号；
根据所述随机延迟时间延迟间隔时间后，基于所述随机幅度调制信号产生对应的发射光脉冲，其中随机幅度调制信号的周期为T，随机延迟时间的间隔时间为0～1/2T。
3.如权利要求1所述的激光雷达测距方法，其特征在于，所述对所述调制信号进行线性插值后做傅里叶变换，得到调制信号的频域信号，包括：
令调制信号为b(n)，其中n表示调制信号中信号点的时刻序号，且n＝1,2,3,…,N，N为调制信号中信号点的总个数；
取线性插值的个数k，在调制信号b(n)的每相邻两个信号点之间进行线性插值生成新的序列c(m)，m＝1,2,3,…,N+k(N-1)，即在b(q)和b(q+1)之间均匀插入k个点，得到b(q),b(q,1),b(q,2),…,b(q,k),b(q+1)，q＝1,2,3,…,N-1；
对序列c(m)进行傅里叶变换生成序列C(m)：
取序列C(m)的共轭序列得到序列D(m)作为调制信号的频域信号。
4.如权利要求3所述的激光雷达测距方法，其特征在于，所述对所述回波信号进行线性插值后进行频谱变换，得到回波信号的频谱信号，包括：
令回波信号为r(n)，其中n表示回波信号中信号点的时刻序号，即为调制信号中信号点的时刻序号，且n＝1,2,3,…,N，N为回波信号中信号点的总个数，即为调制信号中信号点的总个数；
取线性插值的个数k，在回波信号r(n)的每相邻两个信号点之间进行线性插值生成新的序列x(m)，m＝1,2,3,…,N+k(N-1)，即在r(q)和r(q+1)之间均匀插入k个点，得到r(q),r(q,1),r(q,2),…,r(q,k),r(q+1)，q＝1,2,3,…,N-1；
对序列x(m)进行傅里叶变换得到序列X(m)如下：
将序列X(m)作为回波信号的频谱信号。
5.一种激光雷达测距装置，其特征在于，所述激光雷达测距装置包括控制模块、发射模块、接收模块和计算模块，其中：
所述控制模块，用于生成调制信号；
所述发射模块，用于根据所述调制信号调制得到发射光脉冲并发出；
所述接收模块，用于接收回波信号，所述回波信号为待测对象对发射光脉冲进行反射得到；
所述计算模块，用于获取调制信号以及回波信号，对所述调制信号进行线性插值后做傅里叶变换，得到调制信号的频域信号，对所述回波信号进行线性插值后做傅里叶变换，得到回波信号的频域信号，根据调制信号的频域信号和回波信号的频域信号，基于傅里叶逆变换得到待测对象与激光雷达之间的距离；
所述计算模块基于傅里叶逆变换得到待测对象与激光雷达之间的距离，执行如下操作：
取调制信号的频域信号对应的序列D(m)和回波信号的频域信号对应的序列X(m)，将序列D(m)和X(m)相乘得到序列Y(m)，m＝1,2,3,…,N+k(N-1)，N为调制信号中信号点的总个数，k为线性插值时在每相邻两个信号点之间插入点的个数；
对序列Y(m)进行傅里叶逆变换得到序列y(m)；
取序列y(m)中实部最大的元素y(i)，1≤i≤N+k(N-1)，则该元素y(i)对应的时间刻度为i/(N+k(N-1))；
取调制信号的周期为T，则回波信号在激光雷达与待测对象之间的往返时间t为t＝T×i/(N+k(N-1))+Δt，Δt为发射光脉冲发射前对应的随机延迟时间；
得到待测对象与激光雷达之间的距离d为式中c为光速。
6.如权利要求5所述的激光雷达测距装置，其特征在于，所述调制信号为随机幅度调制信号，所述发射模块根据所述调制信号调制得到发射光脉冲，执行如下操作：
生成随机延迟时间和随机幅度调制信号；
根据所述随机延迟时间延迟间隔时间后，基于所述随机幅度调制信号产生对应的发射光脉冲，其中随机幅度调制信号的周期为T，随机延迟时间的间隔时间为0～1/2T。
7.如权利要求5所述的激光雷达测距装置，其特征在于，所述计算模块对所述调制信号进行线性插值后做傅里叶变换，得到调制信号的频域信号，执行如下操作：
令调制信号为b(n)，其中n表示调制信号中信号点的时刻序号，且n＝1,2,3,…,N，N为调制信号中信号点的总个数；
取线性插值的个数k，在调制信号b(n)的每相邻两个信号点之间进行线性插值生成新的序列c(m)，m＝1,2,3,…,N+k(N-1)，即在b(q)和b(q+1)之间均匀插入k个点，得到b(q),b(q,1),b(q,2),…,b(q,k),b(q+1)，q＝1,2,3,…,N-1；
对序列c(m)进行傅里叶变换生成序列C(m)：
取序列C(m)的共轭序列得到序列D(m)作为调制信号的频域信号。
8.如权利要求7所述的激光雷达测距装置，其特征在于，所述计算模块对所述回波信号进行线性插值后进行频谱变换，得到回波信号的频谱信号，执行如下操作：
令回波信号为r(n)，其中n表示回波信号中信号点的时刻序号，即为调制信号中信号点的时刻序号，且n＝1,2,3,…,N，N为回波信号中信号点的总个数，即为调制信号中信号点的总个数；
取线性插值的个数k，在回波信号r(n)的每相邻两个信号点之间进行线性插值生成新的序列x(m)，m＝1,2,3,…,N+k(N-1)，即在r(q)和r(q+1)之间均匀插入k个点，得到r(q),r(q,1),r(q,2),…,r(q,k),r(q+1)，q＝1,2,3,…,N-1；
对序列x(m)进行傅里叶变换得到序列X(m)如下：
将序列X(m)作为回波信号的频谱信号。
CN202110678058.1A 2021-06-18 2021-06-18 一种激光雷达测距方法及装置 Active CN113534179B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110678058.1A CN113534179B (zh) 2021-06-18 2021-06-18 一种激光雷达测距方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110678058.1A CN113534179B (zh) 2021-06-18 2021-06-18 一种激光雷达测距方法及装置

Publications (2)

Publication Number Publication Date
CN113534179A CN113534179A (zh) 2021-10-22
CN113534179B true CN113534179B (zh) 2023-08-01

Family

ID=78125192

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110678058.1A Active CN113534179B (zh) 2021-06-18 2021-06-18 一种激光雷达测距方法及装置

Country Status (1)

Country Link
CN (1) CN113534179B (zh)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101582870A (zh) * 2008-05-15 2009-11-18 中兴通讯股份有限公司 同步实现方法和装置
KR20100071360A (ko) * 2008-12-19 2010-06-29 한국전자통신연구원 Ｍｉｍｏ 및 ｏｆｄｍ 전송 방식을 지원하는 단말 복조 시스템 및 그 제어방법
EP2296362A1 (en) * 2009-09-09 2011-03-16 European Central Bank A method for generating a security bi-level image for a banknote
JP2011058976A (ja) * 2009-09-10 2011-03-24 Fuji Heavy Ind Ltd レーダ装置
EP2860549A1 (en) * 2013-10-12 2015-04-15 Institute of Electronics, Chinese Academy of Sciences Method and device for direction of arrival (DOA) detection of SAR signals
CN106405514A (zh) * 2016-08-23 2017-02-15 中国人民解放军国防科学技术大学 一种非直线航迹条件下合成孔径雷达回波信号的仿真方法
CN107942340A (zh) * 2017-09-28 2018-04-20 北京华航无线电测量研究所 一种太赫兹sar振动补偿成像方法
CN109581350A (zh) * 2018-11-22 2019-04-05 西安八阵图电子科技有限公司 基于时频积分插值的雷达测距测速方法及装置
CN109932718A (zh) * 2019-03-11 2019-06-25 南京航空航天大学 多旋翼无人机载的圆周轨迹环视sar全孔径成像方法
CN110764067A (zh) * 2019-10-15 2020-02-07 哈尔滨工程大学 傅里叶变换的lfm-bpsk复合调制雷达信号参数估计方法
CN111239705A (zh) * 2020-02-12 2020-06-05 北京未感科技有限公司 激光雷达的信号处理方法、装置、设备及存储介质
KR20210022401A (ko) * 2019-08-20 2021-03-03 삼성전자주식회사 라이다 장치 및 그 동작 방법

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101582870A (zh) * 2008-05-15 2009-11-18 中兴通讯股份有限公司 同步实现方法和装置
KR20100071360A (ko) * 2008-12-19 2010-06-29 한국전자통신연구원 Ｍｉｍｏ 및 ｏｆｄｍ 전송 방식을 지원하는 단말 복조 시스템 및 그 제어방법
EP2296362A1 (en) * 2009-09-09 2011-03-16 European Central Bank A method for generating a security bi-level image for a banknote
JP2011058976A (ja) * 2009-09-10 2011-03-24 Fuji Heavy Ind Ltd レーダ装置
EP2860549A1 (en) * 2013-10-12 2015-04-15 Institute of Electronics, Chinese Academy of Sciences Method and device for direction of arrival (DOA) detection of SAR signals
CN106405514A (zh) * 2016-08-23 2017-02-15 中国人民解放军国防科学技术大学 一种非直线航迹条件下合成孔径雷达回波信号的仿真方法
CN107942340A (zh) * 2017-09-28 2018-04-20 北京华航无线电测量研究所 一种太赫兹sar振动补偿成像方法
CN109581350A (zh) * 2018-11-22 2019-04-05 西安八阵图电子科技有限公司 基于时频积分插值的雷达测距测速方法及装置
CN109932718A (zh) * 2019-03-11 2019-06-25 南京航空航天大学 多旋翼无人机载的圆周轨迹环视sar全孔径成像方法
KR20210022401A (ko) * 2019-08-20 2021-03-03 삼성전자주식회사 라이다 장치 및 그 동작 방법
CN110764067A (zh) * 2019-10-15 2020-02-07 哈尔滨工程大学 傅里叶变换的lfm-bpsk复合调制雷达信号参数估计方法
CN111239705A (zh) * 2020-02-12 2020-06-05 北京未感科技有限公司 激光雷达的信号处理方法、装置、设备及存储介质

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于互相关时延估计的激光引信远距离定距;徐伟等;探测与控制学报(第03期);全文 *
基于单个宽带回波提取目标距离和速度的方法;刘明敬;电光与控制;第25卷(第12期);全文 *

Also Published As

Publication number Publication date
CN113534179A (zh) 2021-10-22

Similar Documents

Publication Publication Date Title
AU2010340812B2 (en) Pulse radar range profile motion compensation
JP2006510236A (ja) ディジタル濾波方法
CN111239705B (zh) 激光雷达的信号处理方法、装置、设备及存储介质
CN113156396B (zh) 一种优化干扰源对激光雷达影响的方法与装置
JP2005085167A (ja) 相関処理装置、相関処理方法、パルス圧縮処理装置及びパルス圧縮処理方法
CN110596671A (zh) 一种lfmcw测速测距雷达的优化处理方法及系统
CN115508820A (zh) 一种线性调频脉冲雷达的目标探测方法
CN113805166A (zh) 一种雷达物位计的目标跟踪测距方法及系统
CN113534179B (zh) 一种激光雷达测距方法及装置
CN114035163A (zh) 一种基于算法模型处理雷达回波信号的方法
CN113514842A (zh) 一种距离测量方法、系统及装置
CN111198355B (zh) 基于fpga的pcr回波信号处理系统及方法
CN113835077B (zh) 基于变脉冲重复频率的搜索雷达目标检测方法及系统
CN115561728A (zh) 一种fmcw雷达的单帧速度解模糊方法和装置
CN111736146B (zh) 一种基于速度滤波的双基地检测前跟踪方法及装置
CN114924246A (zh) 目标检测方法、装置、毫米波雷达及存储介质
CN1298107A (zh) 通过雷达频率带宽复用提高分辨率的方法及电路
CN117730264A (zh) 激光雷达抗干扰方法、装置、可读存储介质及终端设备
CN114545353A (zh) 一种基于fri采样和中国余数定理的脉冲多普勒信号欠采样与参数估计方法
CN102680949A (zh) 一种正弦波嵌入式混沌雷达探测装置
CN112630739B (zh) 一种梳状谱信号的宽带条件任意速度副本产生方法
CN108132460A (zh) 基于频域通道均衡的脉冲压缩补偿算法
CN116430354B (zh) 一种fmcw激光雷达目标信息解算方法及系统
CN113567983B (zh) 一种利用滑轨的雷达合成孔径稀疏成像方法及系统
CN110415287B (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