CN110794421A - 伪随机码时延自差分干涉三维成像激光雷达方法及装置 - Google Patents

伪随机码时延自差分干涉三维成像激光雷达方法及装置 Download PDF

Info

Publication number
CN110794421A
CN110794421A CN202010004150.5A CN202010004150A CN110794421A CN 110794421 A CN110794421 A CN 110794421A CN 202010004150 A CN202010004150 A CN 202010004150A CN 110794421 A CN110794421 A CN 110794421A
Authority
CN
China
Prior art keywords
pseudo
random code
differential
phase
random
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.)
Granted
Application number
CN202010004150.5A
Other languages
English (en)
Other versions
CN110794421B (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.)
Sun Jianfeng
Original Assignee
Hangzhou Erae Technology Co 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 Hangzhou Erae Technology Co Ltd filed Critical Hangzhou Erae Technology Co Ltd
Priority to CN202010004150.5A priority Critical patent/CN110794421B/zh
Publication of CN110794421A publication Critical patent/CN110794421A/zh
Application granted granted Critical
Publication of CN110794421B publication Critical patent/CN110794421B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88Lidar systems specially adapted for specific applications
    • G01S17/89Lidar systems specially adapted for specific applications for mapping or imaging

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electromagnetism (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Optical Radar Systems And Details Thereof (AREA)

Abstract

本发明公开了一种伪随机码时延自差分干涉三维成像激光雷达方法,雷达平台中伪随机编码器产生伪随机码经过差分编码器后产生差分伪随机码序列,差分伪随机码序列对输出光束进行相位调制发射至目标,并接收目标的回波信号,回波信号经过非等臂差分干涉仪后进行正交接收,随后经过反正切解相和高通滤波获得位相信息,再进行归一化和二值化处理重建目标时间延迟后的伪随机脉冲序列;将伪随机码与重建得到的目标时间延迟后的伪随机脉冲序列进行互相关处理,获得目标的距离信息,再结合角度信息,经过坐标系变换得到三维位置信息,最终生成三维点云图。本发明可以非常准确获得目标的三维位置信息,具有高抗干扰能力整体系统小型化,运行容易的特点。

Description

伪随机码时延自差分干涉三维成像激光雷达方法及装置
技术领域
本发明涉及激光雷达技术领域,特别涉及一种伪随机码时延自差分干涉三维成像激光雷达方法及装置
背景技术
由于具有高方向性、强抗干扰能力、高测量精度、高时空分辨率等优点,激光雷达技术被广泛应用于航空航天、国民经济和国防等领域,如深空导航、星地同步观测、无人航天器自主交会对接、自主安全着陆、深空探测、车载环境感知和三维地图构建等。传统激光雷达测距方式常采用脉冲激光飞行时间法,即发射较高峰值功率的脉冲,利用激光脉冲往返的时间间隔来获取距离信息,并利用距离的微分来获取速度信息。这种方法结构简单、技术成熟,但在远距离情况下工作时,需要激光器具备极高的峰值功率。为保证激光器的安全运转,出射脉冲的重复频率较低,因而测量的空间分辨率难以提高。
啁啾调频连续波激光雷达对发射激光的频率进行线性调制,回波信号与本振信号进行相干接收,通过获取外差频率实现目标的距离测量,通过波形调制可以实现多普勒速度测量。啁啾调幅连续波激光雷达对发射激光的幅度进行线性调频,将回波强度上调制的延迟啁啾与发射时的初始啁啾混频,得出的差频频率就与回波延迟成正比,通过相干接收可以实现目标的距离和速度的同步探测。这两种技术的优势在于激光器工作在宽脉宽、低峰值功率状态下,可大大提高脉冲重复探测频率,从而提高空间分辨率;另外,采用相干探测可获得更高的灵敏度。但是受激光器硬件条件的限制,成本很高,脉冲重复频率严重受限,啁啾非线性依然是一个未解决的难题。
基于伪随机码和光子计数的测距技术采用直接测量机制,将出射激光按照伪随机码进行强度调制,回波微弱的光子信号序列经过阈值甄别后得到接收的码元序列,该码元序列与原伪随机码序列进行相关运算,相关运算的峰值所对应的距离即为目标距离(在先技术之一,参见 Yufei Zhang, Yan He, Fang Yang, Yuan Luo, and Weibiao Chen,Three-dimensional imaging lidar system based on high speed pseudorandommodulation and photon counting, Chinese Optics Letters, 2016, Vol. 14, No.11, 111101)。该技术的优势在于即使回波信号微弱到光子量级,回波码元出现大量误判,也能够利用回波码元与伪随机码之间的相关性得到正确的目标距离信息,从而大大提高了系统的灵敏度。但是盖革模式的雪崩光电二极管(APD)单光子探测器容易受外界噪声影响,成像速度慢,成本高,限制了这种激光雷达的应用。杨馥等人(在先技术之二,参见邱子胜,杨馥,叶星辰,李姝欣,基于伪随机码相位调制和相干探测的激光测距技术研究,激光与光电子学进展,2018,Vol. 55, No. 5, 052801)提出结合伪随机码相位调制和外差探测技术的方法,能够在宽脉宽情况下获取高的距离分辨率。但是如果激光雷达与测距目标之间存在高速相对运动,会产生较大的多普勒频移,这种多普勒频移会导致相关峰值展宽,造成相关性能的恶化,影响测距精度,所以必须通过一定的方法消除测距过程中目标运动多普勒频移才能保证接收机正常解调。此外,户外大气湍流会对伪随机码相位调制测距产生严重干扰,激光器线宽、激光放大器和位相调制器引入的系统相位噪声都将会严重影响伪随机相位测距的精度。
发明内容
本发明的目的在于,提供一种伪随机码时延自差分干涉三维成像激光雷达方法及装置。本发明可以非常准确获得目标的三维位置信息,具有高抗干扰能力、整体系统小型化,运行容易的特点;此外,本发明不仅可以消除测距过程中目标运动多普勒频移导致的相关峰值展宽,还可以克服大气湍流对伪随机码相位调制测距的干扰,大幅度提高测距精度,同时还降低对激光器和放大器的要求,具有良好的发展前景。
本发明的技术方案:伪随机码时延自差分干涉三维成像激光雷达方法,雷达平台中伪随机编码器产生伪随机码,经过差分编码器后产生差分伪随机码序列,差分伪随机码序列经高速电光相位调制器对激光器的输出光束进行相位调制,得到伪随机码差分相移键控调制激光,将伪随机码差分相移键控调制激光发射至目标并接收目标的回波信号,回波信号经过非等臂差分干涉仪后进行正交接收,得到相干接收信号;将相干接收信号经过反正切解相和高通滤波获得位相信息,再进行归一化和二值化处理重建目标时间延迟后的伪随机脉冲序列;将伪随机码与重建得到的目标时间延迟后的伪随机脉冲序列进行互相关处理,获得远距离目标的高精度距离信息,再结合二维扫描的角度信息,经过坐标系变换得到目标的三维位置信息,最终生成三维点云图。
上述的伪随机码时延自差分干涉三维成像激光雷达方法,所述激光器的输出光束经起偏器起偏,用于提高偏振度,表示为:
Figure DEST_PATH_IMAGE001
其中,为激光器初始相位,
Figure DEST_PATH_IMAGE003
为输出光束主振频率为,
Figure 50590DEST_PATH_IMAGE004
为输出光束线宽引入的 主振频率偏离值,
Figure DEST_PATH_IMAGE005
是时间。
前述的伪随机码时延自差分干涉三维成像激光雷达方法,所述差分伪随机码序列表示为:
Figure 973547DEST_PATH_IMAGE006
其中
Figure DEST_PATH_IMAGE007
是伪随机码,
Figure 952479DEST_PATH_IMAGE008
为伪随机码序列码元宽度,
Figure DEST_PATH_IMAGE009
为码元总数,
Figure 960887DEST_PATH_IMAGE010
是时间延迟为1个比特位的原码,
Figure DEST_PATH_IMAGE011
是差分编码,是 异或运算符。
前述的伪随机码时延自差分干涉三维成像激光雷达方法,所述伪随机码差分相移键控调制激光是通过电光位相调制器把差分编码加载在输出光束的相位上,相邻前后码元的相移为0或π对应要传递的码元0或1。
前述的伪随机码时延自差分干涉三维成像激光雷达方法,所述回波信号经过非等臂差分干涉仪后进行正交接收,得到相干接收信号具体过程是:
所述回波信号表示为:
Figure DEST_PATH_IMAGE013
,
其中,
Figure 526177DEST_PATH_IMAGE014
是振幅,是目标距离,
Figure 944520DEST_PATH_IMAGE016
是激光载频,
Figure DEST_PATH_IMAGE017
为激光器初始相位,
Figure 909065DEST_PATH_IMAGE018
为激光线 宽引入的主振频率偏离值,
Figure DEST_PATH_IMAGE019
是时间,
Figure 500584DEST_PATH_IMAGE020
是目标速度引入的多普勒频移,
Figure DEST_PATH_IMAGE021
是光速,
Figure 132553DEST_PATH_IMAGE022
为大气湍流引入的噪声相位,
Figure DEST_PATH_IMAGE023
为激光放大器引入的噪声相位,
Figure 721798DEST_PATH_IMAGE024
为 位相调制器引入的噪声相位;
回波信号通过非等臂差分干涉仪,首先经过非等臂传输,两个支路的光程延迟对应1比 特位传输数据时间,对应的光程为:
Figure 442148DEST_PATH_IMAGE026
其中
Figure DEST_PATH_IMAGE027
是调制带宽;
两个支路光表示为:
Figure 775040DEST_PATH_IMAGE028
其中
Figure DEST_PATH_IMAGE029
Figure 58254DEST_PATH_IMAGE030
分别是两支路的振幅;
经过非等臂传输后,两个支路光进行干涉,干涉光场表示为:
Figure DEST_PATH_IMAGE031
经过2×4 90º光学桥接器混频后的干涉光场四路输出分别为:
Figure 83979DEST_PATH_IMAGE032
其中是回波信号混频中的噪声位相;
具有正交特性的同相信号和正交信号分别由高速光电平衡探测器进行接收,同相信号和正交信号输出分别为:
其中,
Figure DEST_PATH_IMAGE035
是同相信号的高速光电探测器响应率,
Figure 159699DEST_PATH_IMAGE036
是正交信号的高速光电探测器 响应率。
前述的伪随机码时延自差分干涉三维成像激光雷达方法,将相干接收信号经过反正切解相和高通滤波获得位相信息,再进行归一化和二值化处理重建目标时间延迟后的伪随机脉冲序列的具体过程是:将同相信号和正交信号进行反正切解相重建脉冲位相序列:
Figure DEST_PATH_IMAGE037
其中
Figure 500682DEST_PATH_IMAGE038
为直流量,采用高通滤波器消除;
因此
对重建脉冲位相序列进行二值化处理:
Figure 494045DEST_PATH_IMAGE040
得到重建的目标时间延迟后的伪随机码序列
Figure DEST_PATH_IMAGE041
前述的伪随机码时延自差分干涉三维成像激光雷达方法,所述三维点云图生成的 具体过程是:将原伪随机码序列
Figure 186058DEST_PATH_IMAGE042
经移位延迟
Figure DEST_PATH_IMAGE043
后,与重建的目标时间 延迟后的伪随机码序列进行互相关处理:
Figure 592244DEST_PATH_IMAGE044
当移位延迟
Figure DEST_PATH_IMAGE045
时,相关函数出现峰值,得到远距离目标的视向距离信息
Figure 53312DEST_PATH_IMAGE046
, 距离分辨率为
Figure DEST_PATH_IMAGE047
;在三维空间坐标系
Figure 217577DEST_PATH_IMAGE048
中,,利用二维扫 描得到目标与
Figure 131307DEST_PATH_IMAGE050
轴的夹角
Figure DEST_PATH_IMAGE051
和目标与
Figure 344113DEST_PATH_IMAGE052
轴的夹角
Figure DEST_PATH_IMAGE053
,由此计算出目标点的三维坐标:
最终利用三维坐标得到目标三维点云图。
实现如前述的伪随机码时延自差分干涉三维成像激光雷达方法的装置,包括激光器,激光器的输出端经起偏器连接有高速电光位相调制器,高速电光位相调制器的输出端依次连接有掺饵光纤放大器、光学环形器、二维扫描器和光学望远镜;所述光学环形器的接收端经非等臂差分干涉仪连接有光学桥接器;所述光学桥接器为2×4 90º光学桥接器;所述光学桥接器经高速光电平衡探测器连接有高通滤波器,高通滤波器经高速模数转换器连接有第一高数数据采集器;所述第一高数数据采集器的输出端连接有主控计算机,主控计算机与二维扫描器连接;所述主控计算机还经外触发电路连接有伪随机编码器,伪随机编码器经差分编码器连接有任意波形发生器,任意波形发生器与高速电光位相调制器连接;所述外触发电路还与高速模数转换器连接;所述主控计算机与伪随机编码器之间设有用于采集伪随机码的第二高速数据采集器。
与现有技术相比,本发明通过雷达平台中伪随机编码器产生伪随机码,经过差分编码器后产生差分伪随机码序列,差分伪随机码序列经高速电光相位调制器对激光器的输出光束进行相位调制,得到伪随机码差分相移键控调制激光,将伪随机码差分相移键控调制激光发射至目标并接收目标的回波信号,回波信号经过非等臂差分干涉仪后进行正交接收,得到相干接收信号;将相干接收信号经过反正切解相和高通滤波获得位相信息,再进行归一化和二值化处理重建目标时间延迟后的伪随机脉冲序列;将伪随机码与重建得到的目标时间延迟后的伪随机脉冲序列进行互相关处理,获得远距离目标的高精度距离信息,再结合二维扫描的角度信息,经过坐标系变换得到目标的三维位置信息,最终生成三维点云图;本发明可以非常准确获得目标的三维位置信息,具有高抗干扰能力、整体系统小型化,运行容易的特点;此外,由于输出光束主振频率偏离值和目标速度多普勒频率在1比特位传输数据时间内可以认为是很小的,因此在高速调制的情况下可以消除目标速度多普勒频率和激光器线宽的影响,而且大气湍流引入的噪声相位、激光放大器引入的噪声相位和位相调制器引入的噪声相位都能够完全消除,大幅度提高测距精度,具有非常优秀的测距准确性。经检测本发明调制速率可达到GHz量级,距离分辨率可达到厘米量级,测距精度达到1毫米,系统的重复频率可以达到MHz量级;本发明还降低了对激光器和放大器的要求,有利于减低成本,具有良好的应用前景。
附图说明
图1是本发明的结构原理示意图;
图2是原伪随机码序列、差分伪随机码序列、差分伪随机位相序列、延迟差分伪随机位相序列、自差分干涉位相序列和重建伪随机码序列示意图;
图3是非等臂差分干涉仪示意图;
图4是光学桥接器结构示意图;
图5是空间坐标系中目标二维扫描的角度示意图。
附图中的标记为:1、激光器;2、起偏器;3、高速电光位相调制器;4、掺饵光纤放大器;5、光学环形器;6、二维扫描器;7、光学望远镜;8、非等臂差分干涉仪;9、光学桥接器;10、高速光电平衡探测器;11、高通滤波器;12、高速模数转换器;13、第一高数数据采集器;14、主控计算机;15、外触发电路;16、伪随机编码器;17、差分编码器;18、任意波形发生器;19、第二高速数据采集器。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明,但并不作为对本发明限制的依据。
实施例1:伪随机码时延自差分干涉三维成像激光雷达方法,雷达平台中伪随机编码器产生伪随机码,经过差分编码器后产生差分伪随机码序列,差分伪随机码序列经高速电光相位调制器对激光器的输出光束进行相位调制,得到伪随机码差分相移键控调制激光,将伪随机码差分相移键控调制激光发射至目标并接收目标的回波信号,回波信号经过非等臂差分干涉仪后进行正交接收,得到相干接收信号;将相干接收信号经过反正切解相和高通滤波获得位相信息,再进行归一化和二值化处理重建目标时间延迟后的伪随机脉冲序列;将伪随机码与重建得到的目标时间延迟后的伪随机脉冲序列进行互相关处理,获得远距离目标的高精度距离信息,再结合二维扫描的角度信息,经过坐标系变换得到目标的三维位置信息,最终生成三维点云图。
实施例2:实现伪随机码时延自差分干涉三维成像激光雷达方法的装置,如图1所述,包括激光器1,采用人眼安全的1550nm单频单模连续光纤激光器,激光器线宽200kHz,输出功率为20mW,光纤输出有隔离保护,激光器1的输出端经起偏器2连接有高速电光位相调制器3,高速电光位相调制器3的输出端依次连接有掺饵光纤放大器4、光学环形器5、二维扫描器6和光学望远镜7;所述光学环形器5的接收端经非等臂差分干涉仪8连接有光学桥接器9;所述光学桥接器9为2×4 90º光学桥接器;所述光学桥接器9经高速光电平衡探测器10连接有高通滤波器11,高通滤波器11经高速模数转换器12连接有第一高数数据采集器13;所述第一高数数据采集器13的输出端连接有主控计算机14,主控计算机14与二维扫描器6连接;所述主控计算机14还经外触发电路15连接有伪随机编码器16,伪随机编码器16经差分编码器17连接有任意波形发生器18,任意波形发生器18与高速电光位相调制器3连接;所述外触发电路15还与高速模数转换器12连接;所述主控计算机14与伪随机编码器16之间设有用于采集伪随机码的第二高速数据采集器19。
所述激光器的输出光束经起偏器起偏,用于提高偏振度,确保偏振消光比大于25dB,表示为:
其中,
Figure 723776DEST_PATH_IMAGE002
为激光器初始相位,
Figure 271432DEST_PATH_IMAGE003
为输出光束主振频率为,
Figure 972672DEST_PATH_IMAGE004
为输出光束线宽引入的 主振频率偏离值,
Figure 947581DEST_PATH_IMAGE005
是时间。
伪随机序列既具有确定性的本质,又具有随机信号的某些特征。伪随机相位编码具有调制重频高、保密性强、测距精度与距离无关、峰值旁瓣比高等优点,现在使用最广泛的一种伪随机码序列是最长线性反馈移位寄存器序列,简称m序列。
伪随机码可以记作为
其中,为伪随机码序列码元宽度,
Figure DEST_PATH_IMAGE057
为码元序号,
Figure 733057DEST_PATH_IMAGE058
为码元总数,为伪随机 数,取值为0或1。
伪随机码的主要特性为其自相关函数为
Figure 757644DEST_PATH_IMAGE060
伪随机码的自相关函数是一个尖锐的峰,脉冲宽度越小,波形越尖锐。当
Figure 169034DEST_PATH_IMAGE056
趋近于0 时,上式接近
Figure DEST_PATH_IMAGE061
函数。
伪随机码的互相关特性为
Figure 544652DEST_PATH_IMAGE062
其中,P1和P2为两个独立的伪随机码序列。
因此,利用这个特性则可以提取回波序列位置:用原始伪随机序列与探测器接收到的回波序列进行相关运算,从激光发射到相关峰值对应的时间间隔代表了激光脉冲序列从发出到由目标反射并返回探测器的时间间隔,由此则可给出目标的距离。本实施例中,伪随机码序列由伪随机编码器产生,再经过差分编码器后,产生差分伪随机码序列,通过任意波形发生器驱动电光相位调制器产生伪随机码差分相移键控(DPSK)调制激光。
差分相移键控(Differential Phase Shift Keying, DPSK)调制是一种相干调制方式,通过调制载波的相对相位来加载信息,即把要传递的数字信息加载在载波上相邻码元的相对相位变化的调制,称为相对调相。因此差分编码可以把绝对码转换为相对码(或差分码),编码规则遵循异或原理如图2所示,所述差分伪随机码序列可以表示为:
其中
Figure 699690DEST_PATH_IMAGE064
是伪随机码,
Figure DEST_PATH_IMAGE065
为伪随机码序列码元宽度,
Figure 109942DEST_PATH_IMAGE009
为码元总数,
Figure 426654DEST_PATH_IMAGE010
是时间延迟为1个比特位的原码,
Figure 555147DEST_PATH_IMAGE011
是差分编码,
Figure 248297DEST_PATH_IMAGE012
是 异或运算符。
通过电光位相调制器把差分编码加载在输出光束的相位上,相邻前后码元的相移(即相位改变量)为0或π对应要传递的码元0或1。
所述回波信号表示为:
Figure 44214DEST_PATH_IMAGE013
,
其中,
Figure 797407DEST_PATH_IMAGE014
是振幅,
Figure 475513DEST_PATH_IMAGE066
是目标距离,
Figure 706774DEST_PATH_IMAGE016
是激光载频,为激光器初始相位,
Figure DEST_PATH_IMAGE067
为激光线 宽引入的主振频率偏离值,
Figure 481624DEST_PATH_IMAGE019
是时间,
Figure 850289DEST_PATH_IMAGE020
是目标速度引入的多普勒频移,
Figure 885241DEST_PATH_IMAGE021
是光速,为大气湍流引入的噪声相位,
Figure 281904DEST_PATH_IMAGE023
为激光放大器引入的噪声相位,
Figure 137865DEST_PATH_IMAGE024
为 位相调制器引入的噪声相位;
如图3所示,所述的非等臂差分干涉仪两臂间的光程差为0.3米,对应1个比特位的时间 延迟,回波信号通过非等臂差分干涉仪,首先经过非等臂传输,两个支路的光程延迟对应1 比特位传输数据时间
Figure 710929DEST_PATH_IMAGE025
,对应的光程为:
Figure DEST_PATH_IMAGE069
其中
Figure 539207DEST_PATH_IMAGE070
是调制带宽;
两个支路光表示为:
其中
Figure DEST_PATH_IMAGE071
Figure 148360DEST_PATH_IMAGE030
分别是两支路的振幅;
经过非等臂传输后,两个支路光进行干涉,干涉光场表示为:
Figure 790694DEST_PATH_IMAGE031
再输入至光学桥接器,如图3-4所示,采用自由空间结构2×4 90º光学桥接器,由半波片#1、偏振分束棱镜#1、半波片#2、偏振分束棱镜#2、半波片#5、偏振分束组合棱镜#2、全反射棱镜#1、全反射棱镜#2、半波片#3、四分之一波片、半波片#4和偏振分束组合棱镜#1组成;经过2×4 90º光学桥接器混频后的干涉光场四路输出分别为:
其中
Figure 441435DEST_PATH_IMAGE033
是回波信号混频中的噪声位相;
具有正交特性的同相信号和正交信号分别由高速光电平衡探测器进行接收,高速光电平衡探测器的带宽为1GHz,同相信号和正交信号输出分别为:
Figure DEST_PATH_IMAGE073
其中,
Figure 478180DEST_PATH_IMAGE035
是同相信号的高速光电探测器响应率,
Figure 393047DEST_PATH_IMAGE036
是正交信号的高速光电探测器 响应率。
由高速光电探测器采集同相信号和正交信号,高速光电探测器每个通道采样率5GHz,采集后的同相信号和正交信号经过高通滤波器后由高速模数转换器转换成数字信号,并由第一高速数据采集器采集,与此同时伪随机编码器产生的原伪随机码也有第二高速数据采集器采集,三通道数据输入主控计算机。为保证同步性,采用主控计算机控制的外触发电路同时给高速模数转换器和伪随机编码器提供触发信号,实现伪随机码的同步接收。
在主控计算机中,如果同相信号和正交信号的高速光电探测器响应率保持一致,通过反正切解相重建脉冲位相序列:
Figure 789393DEST_PATH_IMAGE074
其中
Figure DEST_PATH_IMAGE075
为直流量,采用高通滤波器消除;
因此
Figure 865933DEST_PATH_IMAGE076
对重建脉冲位相序列进行二值化处理:
Figure DEST_PATH_IMAGE077
得到重建的目标时间延迟后的伪随机码序列
Figure 121465DEST_PATH_IMAGE078
将伪随机码序列经移位延迟
Figure 105602DEST_PATH_IMAGE043
后,与重建的目标时间延迟后的伪 随机码序列进行互相关处理:
Figure 559717DEST_PATH_IMAGE080
当移位延迟
Figure 338317DEST_PATH_IMAGE045
时,相关函数出现峰值,得到远距离目标的视向距离信息, 距离分辨率为
Figure DEST_PATH_IMAGE081
;在三维空间坐标系
Figure 603393DEST_PATH_IMAGE048
中,
Figure 912015DEST_PATH_IMAGE049
,利用二维 扫描得到目标与
Figure 658254DEST_PATH_IMAGE050
轴的夹角和目标与
Figure 745476DEST_PATH_IMAGE052
轴的夹角
Figure 908604DEST_PATH_IMAGE053
,由此计算出目标点的三维坐标:
Figure 26077DEST_PATH_IMAGE054
最终利用三维坐标得到目标三维点云图。
经测试,本发明的测距分辨率为15cm,测距精度1mm,伪随机码序列的长度1024比特,重复频率为977kHz。
综上所述,本发明可以非常准确获得目标的三维位置信息,具有高抗干扰能力、整体系统小型化,运行容易的特点;此外,本发明不仅可以消除测距过程中目标运动多普勒频移导致的相关峰值展宽,还可以克服大气湍流对伪随机码相位调制测距的干扰,大幅度提高测距精度,同时还降低对激光器和放大器的要求,具有良好的发展前景。

Claims (8)

1.伪随机码时延自差分干涉三维成像激光雷达方法,其特征在于:雷达平台中伪随机编码器产生伪随机码,经过差分编码器后产生差分伪随机码序列,差分伪随机码序列经高速电光相位调制器对激光器的输出光束进行相位调制,得到伪随机码差分相移键控调制激光,将伪随机码差分相移键控调制激光发射至目标并接收目标的回波信号,回波信号经过非等臂差分干涉仪后进行正交接收,得到相干接收信号;将相干接收信号经过反正切解相和高通滤波获得位相信息,再进行归一化和二值化处理重建目标时间延迟后的伪随机脉冲序列;将原伪随机码与重建得到的目标时间延迟后的伪随机脉冲序列进行互相关处理,获得远距离目标的高精度距离信息,再结合二维扫描的角度信息,经过坐标系变换得到目标的三维位置信息,最终生成三维点云图。
2.根据权利要求1所述的伪随机码时延自差分干涉三维成像激光雷达方法,其特征在于:所述激光器的输出光束经起偏器起偏,用于提高偏振度,表示为:
Figure 825695DEST_PATH_IMAGE001
其中,为激光器初始相位,
Figure 460256DEST_PATH_IMAGE003
为输出光束主振频率为,为输出光束线宽引入的 主振频率偏离值,
Figure 404258DEST_PATH_IMAGE005
是时间。
3.根据权利要求1所述的伪随机码时延自差分干涉三维成像激光雷达方法,其特征在于:所述差分伪随机码序列表示为:
其中
Figure 431437DEST_PATH_IMAGE007
是伪随机码,
Figure 102108DEST_PATH_IMAGE008
为伪随机码序列码元宽度,为码元总数,
Figure 436455DEST_PATH_IMAGE010
是时间延迟为1个比特位的原码,是差分编码,是 异或运算符。
4.根据权利要求3所述的伪随机码时延自差分干涉三维成像激光雷达方法,其特征在于:所述伪随机码差分相移键控调制激光是通过电光位相调制器把差分编码加载在输出光束的相位上,相邻前后码元的相移为0或π对应要传递的码元0或1。
5.根据权利要求1所述的伪随机码时延自差分干涉三维成像激光雷达方法,其特征在于:所述回波信号经过非等臂差分干涉仪后进行正交接收,得到相干接收信号具体过程是:
所述回波信号表示为:
Figure 654444DEST_PATH_IMAGE013
,
其中,
Figure 553130DEST_PATH_IMAGE014
是振幅,
Figure 263597DEST_PATH_IMAGE015
是目标距离,是激光载频,
Figure 854295DEST_PATH_IMAGE017
为激光器初始相位,
Figure 127144DEST_PATH_IMAGE018
为激光 线宽引入的主振频率偏离值,
Figure 324908DEST_PATH_IMAGE019
是时间,是目标速度引入的多普勒频移,
Figure 636120DEST_PATH_IMAGE021
是光速,
Figure 876609DEST_PATH_IMAGE022
为大气湍流引入的噪声相位,
Figure 296089DEST_PATH_IMAGE023
为激光放大器引入的噪声相位,为 位相调制器引入的噪声相位;
回波信号通过非等臂差分干涉仪,首先经过非等臂传输,两个支路的光程延迟对应1比 特位传输数据时间
Figure 200252DEST_PATH_IMAGE025
,对应的光程为:
其中
Figure 518418DEST_PATH_IMAGE027
是调制带宽;
两个支路光表示为:
其中
Figure 146026DEST_PATH_IMAGE029
Figure 728317DEST_PATH_IMAGE030
分别是两支路的振幅;
经过非等臂传输后,两个支路光进行干涉,干涉光场表示为:
Figure 122389DEST_PATH_IMAGE031
经过2×4 90º光学桥接器混频后的干涉光场四路输出分别为:
Figure 815538DEST_PATH_IMAGE032
其中
Figure 80298DEST_PATH_IMAGE033
是回波信号混频中的噪声位相;
具有正交特性的同相信号和正交信号分别由高速光电平衡探测器进行接收,同相信号和正交信号输出分别为:
Figure 833490DEST_PATH_IMAGE034
其中,
Figure 714858DEST_PATH_IMAGE035
是同相信号的高速光电探测器响应率,
Figure 211699DEST_PATH_IMAGE036
是正交信号的高速光电探测器 响应率。
6.根据权利要求5所述的伪随机码时延自差分干涉三维成像激光雷达方法,其特征在于:将相干接收信号经过反正切解相和高通滤波获得位相信息,再进行归一化和二值化处理重建目标时间延迟后的伪随机脉冲序列的具体过程是:将同相信号和正交信号进行反正切解相重建脉冲位相序列:
Figure 658861DEST_PATH_IMAGE037
其中为直流量,采用高通滤波器消除;
因此
Figure 686040DEST_PATH_IMAGE039
对重建脉冲位相序列进行二值化处理:
Figure 720992DEST_PATH_IMAGE040
得到重建的目标时间延迟后的伪随机码序列
7.根据权利要求6所述的伪随机码时延自差分干涉三维成像激光雷达方法,其特征在 于:所述三维点云图生成的具体过程是:将原伪随机码序列
Figure 55338DEST_PATH_IMAGE042
经移位延迟
Figure 902510DEST_PATH_IMAGE043
后,与重建的目标时间延迟后的伪随机码序列进行互相关处理:
Figure 741153DEST_PATH_IMAGE044
当移位延迟时,相关函数出现峰值,得到远距离目标的视向距离信息
Figure 163224DEST_PATH_IMAGE046
, 距离分辨率为
Figure 240901DEST_PATH_IMAGE047
;在三维空间坐标系
Figure 883235DEST_PATH_IMAGE048
中,
Figure 97179DEST_PATH_IMAGE049
,利用二维 扫描得到目标与
Figure 533976DEST_PATH_IMAGE050
轴的夹角
Figure 98950DEST_PATH_IMAGE051
和目标与轴的夹角,由此计算出目标点的三维坐标:
Figure 17861DEST_PATH_IMAGE054
最终利用三维坐标得到目标三维点云图。
8.实现如权利要求1-7任一项所述的伪随机码时延自差分干涉三维成像激光雷达方法的装置,其特征在于:包括激光器(1),激光器(1)的输出端经起偏器(2)连接有高速电光位相调制器(3),高速电光位相调制器(3)的输出端依次连接有掺饵光纤放大器(4)、光学环形器(5)、二维扫描器(6)和光学望远镜(7);所述光学环形器(5)的接收端经非等臂差分干涉仪(8)连接有光学桥接器(9);所述光学桥接器(9)为2×4 90º光学桥接器;所述光学桥接器(9)经高速光电平衡探测器(10)连接有高通滤波器(11),高通滤波器(11)经高速模数转换器(12)连接有第一高数数据采集器(13);所述第一高数数据采集器(13)的输出端连接有主控计算机(14),主控计算机(14)与二维扫描器(6)连接;所述主控计算机(14)还经外触发电路(15)连接有伪随机编码器(16),伪随机编码器(16)经差分编码器(17)连接有任意波形发生器(18),任意波形发生器(18)与高速电光位相调制器(3)连接;所述外触发电路(15)还与高速模数转换器(12)连接;所述主控计算机(14)与伪随机编码器(16)之间设有用于采集伪随机码的第二高速数据采集器(19)。
CN202010004150.5A 2020-01-03 2020-01-03 伪随机码时延自差分干涉三维成像激光雷达方法及装置 Active CN110794421B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010004150.5A CN110794421B (zh) 2020-01-03 2020-01-03 伪随机码时延自差分干涉三维成像激光雷达方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010004150.5A CN110794421B (zh) 2020-01-03 2020-01-03 伪随机码时延自差分干涉三维成像激光雷达方法及装置

Publications (2)

Publication Number Publication Date
CN110794421A true CN110794421A (zh) 2020-02-14
CN110794421B CN110794421B (zh) 2020-05-15

Family

ID=69448674

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010004150.5A Active CN110794421B (zh) 2020-01-03 2020-01-03 伪随机码时延自差分干涉三维成像激光雷达方法及装置

Country Status (1)

Country Link
CN (1) CN110794421B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111308489A (zh) * 2020-03-16 2020-06-19 杭州爱莱达科技有限公司 双通道伪随机码相位调制光学域解析聚焦相干激光方法及装置
CN111308488A (zh) * 2020-03-16 2020-06-19 杭州爱莱达科技有限公司 对称伪随机码相位调制全光聚焦相干激光方法及装置
CN111650601A (zh) * 2020-06-01 2020-09-11 杭州爱莱达科技有限公司 车载相干激光雷达高分辨3d成像方法及装置
CN115015886A (zh) * 2022-06-01 2022-09-06 浙江科技学院 双通道伪随机码调相测距测速激光雷达方法及装置
CN115494512B (zh) * 2022-11-15 2023-04-11 中国科学院西安光学精密机械研究所 一种基于伪随机编码的多重频单光子测距方法及系统
CN116338649A (zh) * 2021-12-23 2023-06-27 深圳市速腾聚创科技有限公司 激光雷达抗干扰方法、激光雷达抗干扰装置及电子设备
CN117949927A (zh) * 2024-03-27 2024-04-30 中国科学院西安光学精密机械研究所 基于单光子的天基空间碎片定位方法、系统、介质及设备
CN118011416A (zh) * 2024-04-09 2024-05-10 北京中科国光量子科技有限公司 一种基于稳定相位编解码的激光雷达及相位补偿方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101236253A (zh) * 2008-03-07 2008-08-06 中国科学院上海光学精密机械研究所 高精度测速测距激光雷达系统及测速测距方法
CN104486279A (zh) * 2014-11-27 2015-04-01 西安空间无线电技术研究所 一种uqpsk微波直接调制器调制特性测试方法
CN107991681A (zh) * 2017-11-22 2018-05-04 杭州爱莱达科技有限公司 基于衍射光学的激光雷达及其扫描方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101236253A (zh) * 2008-03-07 2008-08-06 中国科学院上海光学精密机械研究所 高精度测速测距激光雷达系统及测速测距方法
CN104486279A (zh) * 2014-11-27 2015-04-01 西安空间无线电技术研究所 一种uqpsk微波直接调制器调制特性测试方法
CN107991681A (zh) * 2017-11-22 2018-05-04 杭州爱莱达科技有限公司 基于衍射光学的激光雷达及其扫描方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
YANAN ZHI: ""2.5- and 5-Gbps time-delay self-homodyne interference", 《OPTICAL ENGINEERING》 *
YUFEI ZHANG: ""Three-dimension imaging lidar system based on high speed pseudorandom modulation and photon counting"", 《CHINESE OPTICS LETTERS》 *
邱子胜: ""基于伪随机码相位调制和相干探测的激光测距技术研究"", 《激光与光电子学进展》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111308488B (zh) * 2020-03-16 2022-02-11 杭州爱莱达科技有限公司 对称伪随机码相位调制全光聚焦相干激光雷达方法及装置
CN111308488A (zh) * 2020-03-16 2020-06-19 杭州爱莱达科技有限公司 对称伪随机码相位调制全光聚焦相干激光方法及装置
CN111308489A (zh) * 2020-03-16 2020-06-19 杭州爱莱达科技有限公司 双通道伪随机码相位调制光学域解析聚焦相干激光方法及装置
CN111308489B (zh) * 2020-03-16 2022-02-11 杭州爱莱达科技有限公司 伪随机码相位调制激光雷达的光学域解析聚焦方法及装置
CN111650601B (zh) * 2020-06-01 2022-10-25 孙建锋 车载相干激光雷达高分辨3d成像方法及装置
CN111650601A (zh) * 2020-06-01 2020-09-11 杭州爱莱达科技有限公司 车载相干激光雷达高分辨3d成像方法及装置
CN116338649A (zh) * 2021-12-23 2023-06-27 深圳市速腾聚创科技有限公司 激光雷达抗干扰方法、激光雷达抗干扰装置及电子设备
CN115015886A (zh) * 2022-06-01 2022-09-06 浙江科技学院 双通道伪随机码调相测距测速激光雷达方法及装置
CN115015886B (zh) * 2022-06-01 2024-05-24 浙江科技学院 双通道伪随机码调相测距测速激光雷达方法及装置
CN115494512B (zh) * 2022-11-15 2023-04-11 中国科学院西安光学精密机械研究所 一种基于伪随机编码的多重频单光子测距方法及系统
CN117949927A (zh) * 2024-03-27 2024-04-30 中国科学院西安光学精密机械研究所 基于单光子的天基空间碎片定位方法、系统、介质及设备
CN117949927B (zh) * 2024-03-27 2024-08-16 中国科学院西安光学精密机械研究所 基于单光子的天基空间碎片定位方法、系统、介质及设备
CN118011416A (zh) * 2024-04-09 2024-05-10 北京中科国光量子科技有限公司 一种基于稳定相位编解码的激光雷达及相位补偿方法
CN118011416B (zh) * 2024-04-09 2024-06-04 北京中科国光量子科技有限公司 一种基于稳定相位编解码的激光雷达及相位补偿方法

Also Published As

Publication number Publication date
CN110794421B (zh) 2020-05-15

Similar Documents

Publication Publication Date Title
CN110794421B (zh) 伪随机码时延自差分干涉三维成像激光雷达方法及装置
CN110780310B (zh) 偏振分集双通道测速及测距相干激光雷达测量方法及装置
CN110806586B (zh) 无扫描线性调频连续波测速测距激光三维成像方法及装置
CN111337902B (zh) 多通道高重频大动态范围测距测速激光雷达方法及装置
CN115015886B (zh) 双通道伪随机码调相测距测速激光雷达方法及装置
CN106707291B (zh) 一种双频线性调频相干测风激光雷达
CN106932785B (zh) 一种时分复用的偏振相干多普勒测风激光雷达
US5889490A (en) Method and apparatus for improved ranging
CN105629258B (zh) 基于伪随机码相位调制和外差探测的测速测距系统及方法
CN111999739A (zh) 相位调制测距测速相干激光雷达方法及装置
CN111650601B (zh) 车载相干激光雷达高分辨3d成像方法及装置
US20020071109A1 (en) Laser radar system and method
CN111308489B (zh) 伪随机码相位调制激光雷达的光学域解析聚焦方法及装置
CN114035174B (zh) 双通道双啁啾线性调频连续波激光雷达方法及装置
CN105137446A (zh) 基于光学正交解调的相干零差多普勒测速激光雷达系统
CN102928832A (zh) 基于高速伪随机码调制和光子计数的远程激光测距系统
CN103048662B (zh) 一种三波束全光纤相干调频连续波激光雷达
CN105372670A (zh) 基于光学正交解调的相干零差多普勒测速激光雷达系统
CN105204030B (zh) 光学正交解调的相干零差测速激光雷达的数据处理方法
Spollard et al. Mitigation of phase noise and Doppler-induced frequency offsets in coherent random amplitude modulated continuous-wave LiDAR
CN116381644B (zh) 一种偏振不敏感的抗欺骗干扰量子激光雷达
CN116400324B (zh) 一种偏振无关的抗欺骗干扰量子激光雷达
CN115552284A (zh) LiDAR设备和处理
CN110646779A (zh) 相位编码不饱和调制方法、装置、激光雷达测距测速方法及激光雷达系统
Ding et al. Dual-chirp photonics-based radar for distance and velocity measurement based on compressive sensing

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
TR01 Transfer of patent right

Effective date of registration: 20220927

Address after: Room 1803, No. 12, Lane 333, Hongshi Road, Juyuan New District, Jiading District, Shanghai, 201800

Patentee after: Sun Jianfeng

Address before: 310053 room 701, block B, 581 torch Avenue, Puyan street, Binjiang District, Hangzhou City, Zhejiang Province

Patentee before: HANGZHOU ILIDAR TECHNOLOGY CO.,LTD.

TR01 Transfer of patent right