CN106066481A - 一种基于dmd的四维激光雷达成像装置及方法 - Google Patents
一种基于dmd的四维激光雷达成像装置及方法 Download PDFInfo
- Publication number
- CN106066481A CN106066481A CN201610348806.9A CN201610348806A CN106066481A CN 106066481 A CN106066481 A CN 106066481A CN 201610348806 A CN201610348806 A CN 201610348806A CN 106066481 A CN106066481 A CN 106066481A
- Authority
- CN
- China
- Prior art keywords
- signal
- dmd
- frequency
- target
- centerdot
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 35
- 238000000034 method Methods 0.000 title claims abstract description 10
- 238000012545 processing Methods 0.000 claims abstract description 36
- 238000005070 sampling Methods 0.000 claims description 34
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 claims description 26
- 229910052782 aluminium Inorganic materials 0.000 claims description 26
- 239000011159 matrix material Substances 0.000 claims description 22
- 230000003287 optical effect Effects 0.000 claims description 19
- 230000001427 coherent effect Effects 0.000 claims description 16
- 230000003321 amplification Effects 0.000 claims description 9
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 9
- 230000006835 compression Effects 0.000 claims description 7
- 238000007906 compression Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 7
- 238000005286 illumination Methods 0.000 claims description 6
- 230000010355 oscillation Effects 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 6
- 239000000284 extract Substances 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 230000008447 perception Effects 0.000 claims description 4
- 230000011514 reflex Effects 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 7
- 238000001514 detection method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000033001 locomotion Effects 0.000 description 5
- 238000013459 approach Methods 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 230000009017 pursuit movement Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 229910000530 Gallium indium arsenide Inorganic materials 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 230000035945 sensitivity 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/88—Lidar systems specially adapted for specific applications
- G01S17/89—Lidar systems specially adapted for specific applications for mapping or imaging
-
- 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/66—Tracking systems using electromagnetic waves other than radio 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/4802—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- 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/4861—Circuits for detection, sampling, integration or read-out
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Optical Radar Systems And Details Thereof (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种基于DMD的四维激光雷达成像装置及方法,该装置包括发射装置,固定在发射装置出光端的目标,固定于目标反射光路上的第一汇聚透镜,固定于第一汇聚透镜汇聚光路上的数字微镜器件DMD,固定于数字微镜器件DMD反射光路上的第二汇聚透镜,固定于第二汇聚透镜汇聚光路上的接收装置;与接收装置连接的混频器,与混频器连接的中频放大器,与中频放大器连接的正交相位检波器,与正交相位检波器连接的信号处理装置;还包括与所述发射装置、混频器和正交相位检波器连接的高稳定信号源,与发射装置、数字微镜器件DMD和信号处理装置连接的同步装置;本发明还公开了基于DMD的四维激光雷达成像方法;本发明对于距离离散分布和连续分布的目标均具有很好的探测效果,成像速度快,图像分辨率高。
Description
技术领域
本发明涉及激光主动成像技术领域,具体涉及一种基于DMD的四维激光雷达成像装置及方法。
背景技术
激光作为远距离探测的常用光源,具有能量高、方向性好、传播距离远、单色性好等优点,用激光作为光源的成像系统具有成像距离远、抗干扰能力强等优点,被广泛用于战场、考古、测量、检测和城市建模等多个领域。
目前,激光主动成像技术主要包括:扫描式激光主动成像和无扫描激光主动成像两种方案,但这两种方案均存在某些不足之处。
扫描式激光主动成像系统主要利用机械扫描,声(电)光扫描,MEMS扫描等扫描方式实现对目标的三维激光成像。该方案具有探测距离远,信号处理电路简单,成本低廉等优势,但是,扫描式激光主动成像系统的成像速度受扫描速度和激光器重复频率的制约,难以获得高速的激光图像。目前,扫描式激光主动成像的图像更新率一般低于10Hz,难以满足精确制导技术的需求。
无扫描激光主动成像探测系统主要采用基于APD阵列探测器的无扫描激光成像技术,(如MIT的Gen-III系统,我国哈工大,航天二院研制的激光主动成像雷达)。该方案具有结构简单,成像速度快的优点,但是,无扫描激光主动成像系统中的核心器件是APD阵列探测器,这种探测器价格非常昂贵,以Si型32×32APD阵列为例,每只大约需要十几万元人民币,若换做灵敏度更高的军用级的InGaAs型APD阵列探测器,单只价格则会接近百万,高昂的价格严重制约了无扫描激光主动成像技术的研究和应用;此外,受制造工艺和成本的制约,APD阵列的分辨率比较低,一般在32×32左右,难以获得清晰的目标图像。
近几年出现的压缩感知理论为解决上述问题提供了新思路。应用压缩感知理论,可以用单元探测器代替传感器阵列,在采集信号的同时进行信号的压缩,虽然采样频率低于奈奎斯特采样频率,仍旧可以以极大的精确性恢复出原始信号。基于压缩感知的激光成像技术,是近几年随着压缩感知理论的发展与应用而出现的一种新型激光主动成像方式,具有成像速度快,图像分辨率高,成本低等优势,恰恰弥补了传统激光主动成像方案的不足,极有可能为激光主动成像制导提供一条新的技术途径。
发明内容
为了解决上述问题,本发明的目的在于提供一种基于DMD的四维激光雷达成像装置及方法,能够实时地获得目标场景中每一个分辨单元的强度信息和距离信息,不同于目前已有的基于压缩感知的激光主动成像装置对于距离连续分布的目标无法获得较高的距离分辨率,本发明对于距离离散分布和连续分布的目标均具有很好的探测效果,成像速度快,图像分辨率高,成本低,体积小,可以大大提高后续识别、跟踪运动目标的实时性和准确性,增强雷达适应复杂战场的能力。
实现本发明目的的技术思路是:高稳定信号源为发射装置、混频器以及正交相位检波器提供载频一定的正弦信号。同步装置控制发射装置、数字微镜器件DMD以及信号处理装置按周期T同步工作。在每一个周期中,发射装置发射脉冲宽度为t、载频为f0的相参脉冲串信号照在目标上,目标返回来的信号经由第一汇聚透镜传给数字微镜器件DMD,数字微镜器件DMD对照在其上的入射光按照压缩感知二值随机采样矩阵Φ的形式进行调制后传给第二汇聚透镜。第二汇聚透镜将数字微镜器件DMD调制过的光信号传给接收装置,接收装置将载频为f0+fdi(i=1,2,3,…m,m为具有不同运动状态的物体的个数)的光信号转换为电信号,并作为主振信号传给混频器,其中fdi为反应第i个目标运动信息的多普勒频率。在混频器中,主振信号与载频为f0-fI的本振信号进行混频,其中fI为中频信号的频率,并通过选频回路选出频率为fdi+fI的信号输出到中频放大器。中频放大器对其进行功率放大后,进入正交相位检波器,将载频为fdi的信号检出来送入信号处理装置。在信号处理装置中对信号做FFT变换,得到信号的频谱图。信号的频谱图为频率fdi(i=1,2,3,…m,m为具有不同运动状态的物体的个数)处的sinc函数。从幅度谱图中分别提取出w=fdi处的幅值Wi,并计算出所有多普勒频率为fdi的回波合成信号的强度信息Ei。从相位谱图w=fdi处的幅值计算出相应的相位Ψi,得到强度距离混合信息的一次压缩测量值经过M个采样周期,可以生成M维强度/距离测量向量Y,对复数信号恢复效果很好的压缩感知重构算法(如SpaRSA,不限于SpaRSA算法)恢复出四维图像中每一个分辨单元的强度距离混合信息,并据此分别解算出相应的强度信息和距离信息。
为达到以上目的,本发明采用如下技术方案:
一种基于DMD的四维激光雷达成像装置,包括发射装置1,固定在发射装置1出光端的目标10,固定于目标10反射光路上的第一汇聚透镜3,固定于第一汇聚透镜3汇聚光路上的数字微镜器件DMD2,固定于数字微镜器件DMD2反射光路上的第二汇聚透镜4,固定于第二汇聚透镜4汇聚光路上的接收装置5;与接收装置5连接的混频器6,与混频器6连接的中频放大器7,与中频放大器7连接的正交相位检波器8,与正交相位检波器8连接的信号处理装置9;还包括与所述发射装置1、混频器6和正交相位检波器8连接的高稳定信号源11,与所述发射装置1、数字微镜器件DMD2和信号处理装置9连接的同步装置12;
所述发射装置1,发射周期为T、脉冲宽度为t、初相为0的相参脉冲串光信号照在目标10上;
所述数字微镜器件DMD2,依照压缩感知二值随机采样矩阵Φ的形式变换其上的铝镜状态,对目标10的反射光进行调制,压缩感知二值随机采样矩阵Φ的元素值即对应数字微镜器件DMD2的铝镜状态,当矩阵元素值为+1时,数字微镜器件DMD2对应位置的铝镜偏转+12°,将光信号反射到第二汇聚透镜4并传给接收装置5;当矩阵元素值为0时,数字微镜器件DMD 2对应位置的铝镜偏转-12°,将其上的入射光反射到吸收平面进行入射光吸收;
所述接收装置5,采用单点探测器,接收第二汇聚透镜4传来的光信号并转换为电信号以便后续处理;
所述混频器6,是一种典型的频谱搬移电路,它将接收装置5传来的信号作为主振信号,混频器6的主振信号的频率为f0+fdi,与高稳定信号源11提供的频率为f0-fI的本振信号进行混频,采用选频回路选出频率为fdi+fI的中频信号;
所述中频放大器7,对混频器6传来的中频信号进行功率放大,以便进行后续处理;
所述正交相位检波器8,从中频放大器7传来的放大的中频信号中提取出多普勒频率fdi,得到的信号是频率为多普勒频率fdi的复信号;
所述信号处理装置9,对正交相位检波器8传来的复信号进行FFT变换,并运用压缩感知得到目标10的强度信息和距离信息,即可到目标10的四维图像;
所述高稳定信号源11,为发射装置1提供载频为f0的信号,为混频器6提供频率为f0-fI的正弦本振信号,为正交相位检波器8提供频率为fI的相参基准信号;
所述同步装置12,保证发射装置1发射光信号、数字微镜器件DMD2变换其上铝镜状态以及信号处理装置9对这一次采样得到的数据进行处理三者之间的同步进行;
所述正交相位检波器8包括第一乘法器30,与第一乘法器30连接的第一低通滤波器21;第二乘法器31,与第二乘法器31连接的第二低通滤波器22;与第一乘法器30和第二乘法器31连接的移相器20。
本发明提出的基于DMD的四维激光雷达成像方法包括如下步骤:
步骤一,利用信号处理装置9生成一个M×N维的压缩感知二值随机采样矩阵Φ和N×N维的正交稀疏表示矩阵Ψ,其中N=p×q,p为四维图像的纵向分辨单元的个数,q为四维图像的横向分辨单元的个数,M为压缩感知的采样次数,也即相参脉冲串信号的周期数;
步骤二,信号处理装置9由压缩感知二值采样矩阵Φ生成M个测量子矩阵,并编号为1,2,3,…M;其中:M为压缩感知的采样次数;
步骤三,根据四维图像的横向分辨单元的个数p和四维图像的纵向分辨单元的个数q对数字微镜器件DMD2上的铝镜进行分块;
步骤四,发射装置1发射载频为f0、周期为T、脉冲宽度为t、初相为0的相参脉冲串光信号照在目标10上,设当前为第i个周期,并令i=1;
步骤五,目标10的反射光信号经第一汇聚透镜3到达数字微镜器件DMD 2上,同步装置12控制数字微镜器件DMD2按照编号为i的测量子矩阵设置其上的铝镜状态,完成对目标10的反射光的随机采样;
步骤六,数字微镜器件DMD2上的反射光经第二汇聚透镜4到达接收装置5,接收装置5将接收到的光信号转化为电信号,并作为主振信号传给混频器6;
步骤七,混频器6将主振信号与高稳定信号源11提供的频率为f0-fI的正弦本振信号进行混频,采用选频回路选出频率为fdi+fI的中频信号,并传给中频放大器7;
步骤八,中频放大器7对传来的含有多普勒频率fdi的中频信号进行功率放大,并传给正交相位检波器8;
步骤九,正交相位检波器8将中频放大器7传来的中频信号分别与两个相互正交的相参基准信号相乘后,通过低通滤波器21滤除高频分量后,得到频率为多普勒频率fdi的复信号,并传给信号处理装置9;
步骤十,信号处理装置9对第i个周期得到的复信号进行FFT变换,得到第i个幅度频谱图和相位谱图;
步骤十一,令周期序号数i=i+1,并判断i是否小于压缩感知的采样次数M,如果i小于等于M,则转至步骤五,如果i大于M,则转至步骤十二;
步骤十二,信号处理装置9根据步骤十一得到的M个幅度频谱图和相位谱图,通过压缩感知计算得出目标10的强度信息和距离信息,即得到目标10的四维图像。
与现有技术相比,本发明具有以下优点:
第一,由于本发明摒弃了传统的扫描式成像,也没有采用价格昂贵且分辨率较低的APD阵列探测器,转而采用分辨率极高且状态转换速率也极高的数字微镜器件DMD,所以成像速度快,图像分辨率高,成本低,恰恰弥补了传统扫描式以及无扫描式激光主动成像方案的不足。
第二,由于本发明采用FFT变换对接收到的回波信号进行分析,能够实时地获得目标场景中每一个分辨单元的强度信息和距离信息,利用多维信息可以提取出抗干扰性更强的特征,从而大大提高导引头识别、跟踪运动目标的实时性和准确性,显著增强导弹应对复杂战场环境的能力。
第三,由于本发明并未像现有的基于压缩感知的激光主动成像技术依据脉冲飞行时间测量法来判断目标的距离,所以连续分布的目标回波相互混叠的现象对本发明的测距效果无影响,能够显著提高激光图像的距离分辨率,对于距离离散分布和连续分布的目标均能获得良好的成像效果。
附图说明
图1为本发明系统的结构示意图。
图2为本发明系统中正交相位检波器的结构示意图。
图3为本发明方法的实现流程图。
具体实施方式
以下结合附图及具体实施例对本发明作进一步的详细描述。
如图1所示,本发明基于DMD的四维激光雷达成像装置,包括发射装置1,固定在发射装置1出光端的目标10,固定于目标10反射光路上的第一汇聚透镜3,固定于第一汇聚透镜3汇聚光路上的数字微镜器件DMD2,固定于数字微镜器件DMD2反射光路上的第二汇聚透镜4,固定于第二汇聚透镜4汇聚光路上的接收装置5;与接收装置5连接的混频器6,与混频器6连接的中频放大器7,与中频放大器7连接的正交相位检波器8,与正交相位检波器8连接的信号处理装置9;还包括与所述发射装置1、混频器6和正交相位检波器8连接的高稳定信号源11,与所述发射装置1、数字微镜器件DMD2和信号处理装置9连接的同步装置12。
如图2所示,正交相位检波器8,包括第一乘法器30,与第一乘法器30连接的第一低通滤波器21;第二乘法器31,与第二乘法器31连接的第二低通滤波器22;与第一乘法器30和第二乘法器31连接的移相器20。
本发明基于DMD的四维激光雷达成像装置的工作原理如下:
高稳定信号源11为发射装置1、混频器6以及正交相位检波器8提供载频一定的正弦信号。同步装置12控制发射装置1、数字微镜器件DMD2以及信号处理装置9按周期T同步工作。在每一个周期中,发射装置发射脉冲宽度为t、载频为f0的相参脉冲串信号照在目标10上,目标10返回来的信号经由第一汇聚透镜3传给数字微镜器件DMD2,数字微镜器件DMD2对照在其上的入射光按照压缩感知二值采样矩阵Φ的形式进行调制后传给第二汇聚透镜4。第二汇聚透镜4将数字微镜器件DMD2调制过的光信号传给接收装置5,接收装置5将载频为f0+fdi(i=1,2,3,…m,m为具有不同运动状态的物体的个数)的光信号转换为电信号,并作为主振信号传给混频器6,其中fdi为反应第i个目标运动信息的多普勒频率。在混频器6中,主振信号与载频为f0-fI的本振信号进行混频,其中fI为中频信号的频率,并通过选频回路选出频率为fdi+fI的信号输出到中频放大器7。中频放大器7对其进行功率放大后,进入正交相位检波器8,将载频为fdi的信号检出来送入信号处理装置9。在信号处理装置9中对信号做FFT变换,得到信号的频谱图。信号的频谱图为频率fdi(i=1,2,3,…m,m为具有不同运动状态的物体的个数)处的sinc函数。从相位谱图中提取各个不同运动目标的强度信息Ei,从相位谱图中提取相应的相位信息Ψi,并计算一次压缩测量过程中的强度距离混合信息经过M个采样周期,可以得到M维强度/距离测量向量Y,再由对复数信号恢复效果很好的压缩感知重构算法(如SpaRSA,不限于SpaRSA算法)恢复出四维图像中每一个分辨单元的强度距离混合信息,并通过相应算法分别解算出每一个分辨单元对应的强度信息和距离信息,即可得到四维图像。
如图3所示,本发明提出的基于DMD的四维激光雷达成像方法,包括如下步骤:
步骤一,利用信号处理装置9生成一个M×N维的压缩感知二值采样矩阵Φ和N×N维的正交稀疏表示矩阵Ψ,其中N=p×q,p为四维图像的纵向分辨单元的个数,q为四维图像的横向分辨单元的个数,M为压缩感知的采样次数,也即相参脉冲串信号的周期数;
其中bist为Φ中第i行,第st列的元素,1≤s≤p,1≤t≤q,bist的值取+1或0;
步骤二,信号处理装置9由压缩感知二值采样矩阵Φ生成M个测量子矩阵Φi,并令编号i为1,2,3,…M,其中:M为压缩感知的采样次数;每个随机测量子矩阵Φi的形式如下所示:
第i个测量子矩阵Φi即对应原压缩感知二值采样矩阵Φ的第i行元素;
步骤三,根据四维图像的横向分辨单元个数p和四维图像的纵向分辨单元个数q对数字微镜器件DMD2上的铝镜进行分块;
设数字微镜器件DMD2上的铝镜个数为m×n,m为数字微镜器件中横向铝镜的个数,n为数字微镜器件中纵向铝镜的个数;将数字微镜器件DMD2上的铝镜分为p×q块,则每块有[min{m,n}/max{p,q}]2个铝镜。以数字微镜器件DMD2上铝镜个数为1024*768和图像分辨率为64*64为例,将数字微镜器件DMD2上的铝镜分为64*64块,每块有144个铝镜,对应于目标场景中的一个分辨单元;
步骤四,发射装置1发射载频为f0、周期为T、脉冲宽度为t、初相为0的相参脉冲串光信号照在目标10上,设当前为第i个周期,并令i=1;
步骤五,目标10的反射光信号经第一汇聚透镜3到达数字微镜器件DMD2上,同步装置12控制数字微镜器件DMD2按照编号为i的测量子矩阵设置其上的铝镜状态,完成对目标10的反射光的随机采样,当测量子矩阵相应元素值为+1时,将数字微镜器件DMD2对应位置的铝镜偏转+12°,将其上的入射光反射给第二汇聚透镜4并传给接收装置5;当随机矩阵相应元素值为0时,将数字微镜器件DMD2对应位置的铝镜偏转-12°,将其上的入射光反射到吸收平面进行入射光吸收;
步骤六,数字微镜器件DMD2上的反射光经第二汇聚透镜4到达接收装置5,接收装置5将接收到的光信号转化为电信号,并作为主振信号传给混频器6;
步骤七,混频器6将主振信号与高稳定信号源11提供的频率为f0-fI的本振信号进行混频,采用选频回路选出频率为fdi+fI的中频信号,并传给中频放大器7,具体如下:
接收装置5传来的主振信号可以表示为其中cij代表第i个运动目标上的第j点反射信号的强度信息,代表第i个运动目标上的第j点反射信号的距离信息。高稳定信号源11提供的频率为f0-fI的本振信号可以表示为cos(2π(f0-fI)t),两者混频后可得到如下信号:
选频回路的低通滤波器滤除掉频率为2f0+fdi-fI的信号后,保留频率为fdi+fI的中频信号,传给中频放大器7;
步骤八,中频放大器7对传来的含有多普勒频率fdi的中频信号进行功率放大,增益为K,并传给正交相位检波器8;
步骤九,正交相位检波器8将中频放大器7传来的中频信号分别与两个相互正交的相参基准信号相乘后,通过第一低通滤波器21和第二低通滤波器22滤除高频分量后,得到频率为多普勒频率fdi的复信号,并传给信号处理装置9,具体如下:
设相参基准信号sc(t)=cos(2πfIt),其经过移相器20,可得到与其正交的另一个相参基准信号sc'(t)=cos(2πfIt+π/2)=-sin(2πfIt)。将中频信号与相参基准信号sc(t)=cos(2πfIt)相乘,可得信号为:
第一低通滤波器21滤除掉其中的高频分量,可得同相通道信号为
将中频信号与相参基准信号sc'(t)=-sin(2πfIt)相乘,可得信号为:
第二低通滤波器22滤除掉其中的高频分量,可得正交通道信号为:
由此可获得正交双通道信号X(t)为:
步骤十,信号处理装置9对第i个周期得到的复信号进行FFT变换,得到第i个幅度频谱图和相位谱图;
步骤十一,令周期序号数i=i+1,并判断i是否小于压缩感知的采样次数M,如果i小于等于M,则转至步骤五,如果i大于M,则转至步骤十二;
步骤十二,信号处理装置9根据步骤十一得到的M个幅度频谱图和相位谱图,通过压缩感知计算得出目标10的四维图像,具体如下:
12.1)选取M个幅度频谱图中频率为fdi处的幅值Wi,并用公式Ei=Wi/τ分别计算出相应的Ei,其中τ已知,为相参脉冲串光信号的脉冲宽度;
12.2)计算一次采样过程中强度距离混合信息的一次压缩测量值同时有Aij反映了第i个运动目标第j点反射信号的强度信息,反映了第i个运动目标第j点反射信号的距离信息。其中aij和bij由一次采样过程中第i个运动目标第j点反射信号的距离信息和强度信息共同决定;
12.3)经过M个采样周期,即可得到M个这样的测量值,并由此生成M维强度/距离测量向量Y。同时计算压缩感知的全息算子A=ΦΨ,其中Φ为M×N维的压缩感知二值采样矩阵,Ψ为N×N维的正交稀疏表示矩阵;
12.4)根据公式Y=Ah,通过对复数信号恢复效果很好的压缩感知重构算法SpaRSA计算出h,进而由公式C=Ψh计算出目标10每一点处的距离信息和强度信息,C为N×1维的列向量,其中N=p×q:
C=[C11 C12...C1q C21 C22...C2q....Cp1..Cpq]T
变换得到:
C中的每一个元素都是由目标10中每一个分辨单元对应的强度信息Dij和距离信息ψij共同决定的,即
12.5)根据如下公式计算出目标10中每一个分辨单元对应的距离信息ψij:
12.6)根据如下公式计算出目标10中每一个分辨单元对应的强度信息Dij:
12.7)得到了目标10中每一个分辨单元对应的距离信息和强度信息,即得到了目标10的四维图像。
Claims (2)
1.一种基于DMD的四维激光雷达成像装置,其特征在于:包括发射装置(1),固定在发射装置(1)出光端的目标(10),固定于目标(10)反射光路上的第一汇聚透镜(3),固定于第一汇聚透镜(3)汇聚光路上的数字微镜器件DMD(2),固定于数字微镜器件DMD(2)反射光路上的第二汇聚透镜(4),固定于第二汇聚透镜(4)汇聚光路上的接收装置(5);与接收装置(5)连接的混频器(6),与混频器(6)连接的中频放大器(7),与中频放大器(7)连接的正交相位检波器(8),与正交相位检波器(8)连接的信号处理装置(9);还包括与所述发射装置(1)、混频器(6)和正交相位检波器(8)连接的高稳定信号源(11),与所述发射装置(1)、数字微镜器件DMD(2)和信号处理装置(9)连接的同步装置(12);
所述发射装置(1),发射周期为T、脉冲宽度为t、初相为0的相参脉冲串光信号照在目标(10)上;
所述数字微镜器件DMD(2),依照压缩感知二值随机采样矩阵Φ的形式变换其上的铝镜状态,对目标(10)的反射光进行调制,压缩感知二值采样矩阵Φ的元素值即对应数字微镜器件DMD(2)的铝镜状态,当矩阵元素值为+1时,数字微镜器件DMD(2)对应位置的铝镜偏转+12°,将光信号反射到第二汇聚透镜4并传给接收装置(5);当矩阵元素值为0时,数字微镜器件DMD(2)对应位置的铝镜偏转-12°,将其上的入射光反射到吸收平面进行入射光吸收;
所述接收装置(5),采用单点探测器,接收第二汇聚透镜(4)传来的光信号并转换为电信号以便后续处理;
所述混频器(6),是一种典型的频谱搬移电路,它将接收装置(5)传来的信号作为主振信号,混频器(6)的主振信号的频率为f0+fdi,与高稳定信号源(11)提供的频率为f0-fI的本振信号进行混频,采用选频回路选出频率为fdi+fI的中频信号;
所述中频放大器(7),对混频器(6)传来的中频信号进行功率放大,以便进行后续处理;
所述正交相位检波器(8),从中频放大器(7)传来的放大的中频信号中提取出多普勒频率fdi,得到的信号是频率为多普勒频率fdi的复信号;
所述信号处理装置(9),对正交相位检波器(8)传来的复信号进行FFT变换,并运用压缩感知得到目标(10)的强度信息和距离信息,即可到目标(10)的四维图像;
所述高稳定信号源(11),为发射装置(1)提供载频为f0的信号,为混频器6提供频率为f0-fI的正弦本振信号,为正交相位检波器(8)提供频率为fI的相参基准信号;
所述同步装置(12),保证发射装置(1)发射光信号、数字微镜器件DMD(2)变换其上铝镜状态以及信号处理装置(9)对这一次采样得到的数据进行处理三者之间的同步进行。
所述正交相位检波器(8)包括第一乘法器(30),与第一乘法器(30)连接的第一低通滤波器(21);第二乘法器(31),与第二乘法器(31)连接的第二低通滤波器(22);与第一乘法器(30)和第二乘法器(31)连接的移相器(20)。
2.权利要求1基于DMD的四维激光雷达成像装置的成像方法,其特征在于:包括如下步骤:
步骤一,利用信号处理装置(9)生成一个M×N维的压缩感知二值采样矩阵Φ和N×N维的正交稀疏表示矩阵Ψ,其中N=p×q,p为四维图像的纵向分辨单元的个数,q为四维图像的横向分辨单元的个数,M为压缩感知的采样次数,也即相参脉冲串信号的周期数;
步骤二,信号处理装置(9)由压缩感知二值采样矩阵Φ生成M个测量子矩阵,并编号为1,2,3,…M;其中:M为压缩感知的采样次数;
步骤三,根据四维图像的横向分辨单元的个数p和四维图像的纵向分辨单元的个数q对数字微镜器件DMD(2)上的铝镜进行分块;
步骤四,发射装置1发射载频为f0、周期为T、脉冲宽度为t、初相为0的相参脉冲串光信号照在目标(10)上,设当前为第i个周期,并令i=1;
步骤五,目标(10)的反射光信号经第一汇聚透镜(3)到达数字微镜器件DMD(2)上,同步装置(12)控制数字微镜器件DMD(2)按照编号为i的测量子矩阵设置其上的铝镜状态,完成对目标(10)的反射光的随机采样;
步骤六,数字微镜器件DMD(2)上的反射光经第二汇聚透镜(4)到达接收装置(5),接收装置(5)将接收到的光信号转化为电信号,并作为主振信号传给混频器6;
步骤七,混频器(6)将主振信号与高稳定信号源(11)提供的频率为f0-fI的正弦本振信号进行混频,采用选频回路选出频率为fdi+fI的中频信号,并传给中频放大器(7);
步骤八,中频放大器(7)对传来的含有多普勒频率fdi的中频信号进行功率放大,并传给正交相位检波器(8);
步骤九,正交相位检波器(8)将中频放大器(7)传来的中频信号分别与两个相互正交的相参基准信号相乘后,通过低通滤波器(21)滤除高频分量后,得到频率为多普勒频率fdi的复信号,并传给信号处理装置(9);
步骤十,信号处理装置(9)对第i个周期得到的复信号进行FFT变换,得到第i个幅度频谱图和相位谱图;
步骤十一,令周期序号数i=i+1,并判断i是否小于压缩感知的采样次数M,如果i小于等于M,则转至步骤五,如果i大于M,则转至步骤十二;
步骤十二,信号处理装置(9)根据步骤十一得到的M个幅度频谱图和相位谱图,通过压缩感知计算得出目标(10)的强度信息和距离信息,即得到目标(10)的四维图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610348806.9A CN106066481B (zh) | 2016-05-24 | 2016-05-24 | 一种基于dmd的四维激光雷达成像装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610348806.9A CN106066481B (zh) | 2016-05-24 | 2016-05-24 | 一种基于dmd的四维激光雷达成像装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106066481A true CN106066481A (zh) | 2016-11-02 |
CN106066481B CN106066481B (zh) | 2019-01-18 |
Family
ID=57421045
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610348806.9A Active CN106066481B (zh) | 2016-05-24 | 2016-05-24 | 一种基于dmd的四维激光雷达成像装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106066481B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107402385A (zh) * | 2017-07-31 | 2017-11-28 | 岭纬公司 | 适用于激光雷达的dlp信号接收系统及其自适应调节方法 |
WO2019010679A1 (zh) * | 2017-07-13 | 2019-01-17 | 中山大学 | 一种基于降维稀疏表示的主动式场景三维信息获取方法 |
CN110471082A (zh) * | 2019-08-13 | 2019-11-19 | 西安电子科技大学 | 单像素激光计算成像装置及方法 |
CN113050118A (zh) * | 2021-05-31 | 2021-06-29 | 之江实验室 | 单臂激光外差关联成像雷达空间后端调制方法及装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5734347A (en) * | 1996-06-10 | 1998-03-31 | Mceligot; E. Lee | Digital holographic radar |
CN202485671U (zh) * | 2010-11-26 | 2012-10-10 | 北京空间机电研究所 | 基于压缩感知理论的航天相机成像装置 |
CN103234479A (zh) * | 2013-04-10 | 2013-08-07 | 浙江大学 | 基于压缩感知脉冲编码调制的三维信息获取装置和方法 |
CN104142506A (zh) * | 2014-08-15 | 2014-11-12 | 中国科学院上海技术物理研究所 | 一种基于压缩感知的激光雷达成像系统 |
CN105223582A (zh) * | 2015-09-01 | 2016-01-06 | 西安交通大学 | 一种基于压缩感知的激光雷达成像装置及成像方法 |
CN105242281A (zh) * | 2015-09-01 | 2016-01-13 | 西安交通大学 | 基于apd阵列的三维激光成像系统及方法 |
-
2016
- 2016-05-24 CN CN201610348806.9A patent/CN106066481B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5734347A (en) * | 1996-06-10 | 1998-03-31 | Mceligot; E. Lee | Digital holographic radar |
CN202485671U (zh) * | 2010-11-26 | 2012-10-10 | 北京空间机电研究所 | 基于压缩感知理论的航天相机成像装置 |
CN103234479A (zh) * | 2013-04-10 | 2013-08-07 | 浙江大学 | 基于压缩感知脉冲编码调制的三维信息获取装置和方法 |
CN104142506A (zh) * | 2014-08-15 | 2014-11-12 | 中国科学院上海技术物理研究所 | 一种基于压缩感知的激光雷达成像系统 |
CN105223582A (zh) * | 2015-09-01 | 2016-01-06 | 西安交通大学 | 一种基于压缩感知的激光雷达成像装置及成像方法 |
CN105242281A (zh) * | 2015-09-01 | 2016-01-13 | 西安交通大学 | 基于apd阵列的三维激光成像系统及方法 |
Non-Patent Citations (4)
Title |
---|
凡守涛: ""单脉冲气象雷达系统分析与设计"", 《中国博士学位论文全文数据库 信息科技辑》 * |
刘丽华等: ""连续波多普勒测速雷达射频前端电路设计与仿真"", 《电子科技》 * |
迪可新等: ""基于压缩感知的空间稀疏目标成像方法研究"", 《微电子学与计算机》 * |
邓正芳等: ""微镜阵列激光成像技术及其性能分析"", 《激光与红外》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019010679A1 (zh) * | 2017-07-13 | 2019-01-17 | 中山大学 | 一种基于降维稀疏表示的主动式场景三维信息获取方法 |
CN107402385A (zh) * | 2017-07-31 | 2017-11-28 | 岭纬公司 | 适用于激光雷达的dlp信号接收系统及其自适应调节方法 |
CN110471082A (zh) * | 2019-08-13 | 2019-11-19 | 西安电子科技大学 | 单像素激光计算成像装置及方法 |
CN113050118A (zh) * | 2021-05-31 | 2021-06-29 | 之江实验室 | 单臂激光外差关联成像雷达空间后端调制方法及装置 |
CN113050118B (zh) * | 2021-05-31 | 2021-08-13 | 之江实验室 | 单臂激光外差关联成像雷达空间后端调制方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106066481B (zh) | 2019-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106066481A (zh) | 一种基于dmd的四维激光雷达成像装置及方法 | |
Kusy et al. | Radio interferometric tracking of mobile wireless nodes | |
CN103592026B (zh) | 基于压缩感知和编码变换的时间飞行成像光谱系统及方法 | |
CN101587500B (zh) | 双站合成孔径雷达海面成像的计算机仿真方法 | |
CN112558066B (zh) | 双极化sar图像系统 | |
CN105242281B (zh) | 基于apd阵列的三维激光成像系统及方法 | |
CN107192355B (zh) | 一种双光梳精密测角方法及测角系统 | |
CN103994765B (zh) | 一种惯性传感器的定位方法 | |
CN101937072A (zh) | 基于量子特性的卫星定位系统与方法 | |
CN103472457A (zh) | 稀疏孔径压缩计算关联飞行时间的三维成像系统及方法 | |
CN104142506B (zh) | 一种基于压缩感知的激光雷达成像系统 | |
CN105929399A (zh) | 一种干涉sar数据成像及高程估计方法 | |
CN103472455A (zh) | 稀疏孔径压缩计算关联飞行时间四维光谱成像系统及方法 | |
CN109407094A (zh) | 车载超宽带雷达前视成像系统 | |
CN106802403A (zh) | 声矢量传感器二维阵列music解相干参数估计方法 | |
CN107300697A (zh) | 基于无人机的运动目标ukf滤波方法 | |
CN105403731B (zh) | 一种高精度宽量程的飞秒激光多普勒三维流速传感方法 | |
CN105404128B (zh) | 多帧相移数字全息方法及装置 | |
CN108254760A (zh) | 一种基于三颗量子卫星的定位与导航方法与系统 | |
CN106774901A (zh) | 基于超声定位的远距离pc体感输入方法 | |
CN110161501A (zh) | 一种多时序sar图像的目标区域地表起伏信息提取方法 | |
CN106019274A (zh) | 一种新型多普勒雷达成像装置及方法 | |
CN106370302A (zh) | 基于时空联合调制和mzi光波导阵列的成像光谱仪 | |
Anghel et al. | Compact ground-based interferometric synthetic aperture radar: Short-range structural monitoring | |
CN106197391A (zh) | 基于弱测量技术的光子陀螺仪 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address | ||
CP03 | Change of name, title or address |
Address after: Room 820, 8th floor, Zhongchuang Plaza, No. 385 Aerospace Middle Road, Aerospace Base, Xi'an City, Shaanxi Province, 710100 Patentee after: Xi'an Shenrui Intelligent Technology Co.,Ltd. Country or region after: China Address before: 710199 Shaanxi Xi'an Aerospace base 385 space Road 8, 820 floor, Chong Chuang square. Patentee before: XI'AN SHENQIONG PHOTOELECTRIC TECHNOLOGY CO.,LTD. Country or region before: China |