CN103234479A - 基于压缩感知脉冲编码调制的三维信息获取装置和方法 - Google Patents

基于压缩感知脉冲编码调制的三维信息获取装置和方法 Download PDF

Info

Publication number
CN103234479A
CN103234479A CN201310125563.9A CN201310125563A CN103234479A CN 103234479 A CN103234479 A CN 103234479A CN 201310125563 A CN201310125563 A CN 201310125563A CN 103234479 A CN103234479 A CN 103234479A
Authority
CN
China
Prior art keywords
centerdot
modulation
phi
dimensional
image
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
CN201310125563.9A
Other languages
English (en)
Other versions
CN103234479B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201310125563.9A priority Critical patent/CN103234479B/zh
Publication of CN103234479A publication Critical patent/CN103234479A/zh
Application granted granted Critical
Publication of CN103234479B publication Critical patent/CN103234479B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)
  • Optical Radar Systems And Details Thereof (AREA)
  • Measurement Of Optical Distance (AREA)

Abstract

本发明公开了一种基于压缩感知脉冲编码调制的三维信息获取方法,包括:脉冲激光光源发出脉冲激光照射目标景物并成像到空间光调制器上,然后根据加载其上的二维随机调制图像对目标景物图像进行编码调制;并将调制后图像的光强会聚后得到调制后的不同距离回波的总光强,并由单点探测器接收;利用高速数据采集模块采集单点探测器探测的不同距离回波的总光强;再用一维随机调制矩阵调制不同距离回波的总光强并求和,得到两次调制后的所有距离回波的总光强,重复上述步骤对目标景物进行多次测量,得到多次测量的回波总光强,最后利用压缩感知算法进行一次计算完成目标景物三维信息的获取。本发明还公开了一种基于压缩感知脉冲编码调制的三维信息获取装置。

Description

基于压缩感知脉冲编码调制的三维信息获取装置和方法
技术领域
本发明涉及三维信息获取领域,尤其涉及一种基于压缩感知脉冲编码调制的三维信息获取装置和方法。
背景技术
三维信息的获取如今已经广泛用于三维建模,目标探测等领域,并且已经发展出来各种各样的方法,各有其优缺点及使用范围。如今比较流行的有:
基于计算机视觉的方法,可以从二维的图形获取三维信息,可以实现对大视场三维信息的获取;利用飞行时间法,在一定范围内对目标物体进行扫描,在获取二维信息的同时得到探测脉冲的飞行时间,从而获取目标物体的深度信息,其图像获取速度和分辨率受限于扫描速度和扫描点数;使用结构光对于目标物体进行主动式的测量,扫描速度快,测量精度高,比较适用于室内物体表面反射情况比较好的场合;在远距离三维信息获取上主要采用切片法,对回波信号切片后分别计算不同距离的二维信息,通过多次计算完成三维重构。
基于计算机视觉的方法、结构光法都不适用于远距离三维信息的获取,飞行时间法受限于扫描速度和扫描点数,并且容易受到传播过程中环境的影响,远距离三维信息的获取结果不理想,而切片法需要根据切片次数多次进行压缩感知计算,计算速度慢且繁琐。
压缩感知理论是关于信号采样和数据处理的一门新兴的理论,它指出在一定的条件下,用远低于Nyquist采样定理要求的采样次数对信号进行采样时,也能很好的恢复出原始信号,它可以将图像的压缩采集和压缩过程合二为一。通过M次的测量可以还原N维的信号(M<N)。压缩感知理论的核心内容有两点:目标信号的稀疏表达和测量的非相干性。对于有长度为N的实数信号x(n),可以进行稀疏变换
Figure BDA00003031784600021
或者x=Ψθ。Ψ为相应的稀疏基矩阵。压缩感知并不直接对信号x(n)进行测量,而是通过一个随机投影矩阵Φ进行测量y=Φx。Φ是一个M×N维的矩阵,每一行是一个基向量
Figure BDA00003031784600022
表示对信号x(n)进行一次线性的测量。M表示测量次数,并且满足M<N。由于x可以在Ψ域进行稀疏表示,所以上式也可以表示为y=Φx=ΦΨTθ。求解此方程的问题可以表示为求最小1范数的优化问题:
Figure BDA00003031784600023
约束条件为
Figure BDA00003031784600024
可用的算法有基追踪算法,贪婪追踪算法,凸松弛法,组合算法,TV算法等。
目前还没有将压缩感知理论应用于空间三维编码调制的三维信息获取领域的文献报道。
发明内容
本发明基于压缩感知使用一个空间光调制器对图像二维编码调制并利用一维随机调制向量对不同距离总光强进行编码调制,最后通过一次压缩感知算法实现三维信息获取。
一种基于压缩感知脉冲编码调制的三维信息获取装置,包括:脉冲激光模块、图像生成加载模块、成像系统、空间光调制器、会聚透镜、单点探测器、高速数据采集模块和计算存储模块;
所述脉冲激光模块发出脉冲激光并投射到目标景物上;
所述图像生成加载模块生成多次测量加载的二维随机调制图像,并将二维随机调制图像传递给所述空间光调制器;
所述成像系统用于收集所述目标景物的漫反射光,并将目标景物成像到空间光调制器上;
所述空间光调制器根据加载其上的二维随机调制图像对得到的目标景物的图像进行编码调制;
所述会聚透镜将目标景物图像调制后的光强会聚得到调制后的不同距离回波的总光强;
所述单点探测器用于探测所述会聚透镜会聚得到不同距离回波的总光强;
所述高速数据采集模块用于采集单点探测器多次探测得到的不同距离回波的总光强;
所述计算存储模块根据多次探测得到的不同距离回波的总光强,利用压缩感知算法完成三维信息获取。
所述空间光调制器为透射型空间光调制器或反射型空间光调制器。反射型空间光调制器可以为反射型的硅上液晶器件(LCOS)或者反射型的数字微镜器件(DMD)等,优选为反射型的数字微镜器件(DMD)。
所述成像系统为望远系统、显微系统或照相系统。成像系统用于采集目标景物的漫反射光,并将目标景物成像在空间光调制器上。
所述图像生成加载模块内设有同步模块,实现每切换一张图像脉冲激光模块发一个脉冲激光,由于需要多次采样,同步模块实现图像生成加载模块和脉冲激光的同步。
所述高速数据采集模块和脉冲激光模块之间连接有同步探测模块。
高速数据采集模块为两路信号同时采集,同步探测模块实现脉冲激光模块每发一个脉冲光后两路同步采集光强数据,所述同步探测模块包括探测器和收集透镜或透镜组。
所述脉冲激光模块包括脉冲激光光源、扩束镜和分光棱镜。
本发明还提供了一种基于压缩感知脉冲编码调制的三维信息获取方法,包括以下步骤:
(1)脉冲激光光源发出脉冲激光照射目标景物;
(2)目标景物经成像系统成像到空间光调制器上,空间光调制器根据加载其上的二维随机调制图像对得到的目标景物的图像进行编码调制;
(3)将目标景物调制后图像的光强会聚后得到调制后的不同距离回波的总光强,并由单点探测器接收;
(4)利用高速数据采集模块采集单点探测器探测的不同距离回波的总光强;
(5)重复步骤(1)~(4)对目标景物进行M次测量,得到M次测量的回波总光强,最后利用压缩感知算法完成目标景物三维信息的获取。
在步骤(4)中,需要对高速数据采集模块每次采集的不同距离回波的总光强进行一维随机编码调制,其中,一维随机调制向量为B=1×r,r为一维随机调制向量元素数。
所述M为小于二维随机调制图像总像素数和一维随机调制向量元素数的乘积的自然数。
在步骤(5)中,利用压缩感知算法完成目标景物三维信息获取的具体方法为:
(1)设定二维随机调制图像总像素数A=p×q,p为二维随机调制图像的横向分辨率,q为二维随机调制图像的纵向分辨率;一维随机调制向量B=1×r,目标景物的测量次数为M,生成一个M×A维的矩阵Φ1和一个M×B维的矩阵Φ2,利用矩阵Φ1生成M幅二维随机调制图像Φ1(i),利用矩阵Φ2生成M个一维随机调制向量Φ2(j),其中i,j=1,2,3,4…M;
&Phi; 1 = A 1 ( 1,1 ) A 1 ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; A 1 ( p , q ) A 2 ( 1,1 ) A 2 ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; A 2 ( p , q ) &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; A M ( 1,1 ) A M ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; A M ( p , q ) = &Phi; 1 ( 1 ) &Phi; 1 ( 2 ) &CenterDot; &CenterDot; &CenterDot; &Phi; 1 ( M )
&Phi; 2 = B 1 ( 1,1 ) B 1 ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; B 1 ( 1 , r ) B 2 ( 1,1 ) B 2 ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; B 2 ( 1 , r ) &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; B M ( 1,1 ) B M ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; B M ( 1 , r ) = &Phi; 2 ( 1 ) &Phi; 2 ( 2 ) &CenterDot; &CenterDot; &CenterDot; &Phi; 2 ( M )
其中,Ai(m,n)为第i幅二维随机调制图像中的图像元素,m=1,2,3…p,n=1,2,3…q;Bi(1,s)为第i个一维随机调制向量中的元素,s=1,2,3…r;
(2)空间光调制器依次在脉冲信号光中加载第i幅二维随机调制图像,得到带有第i幅二维随机调制图像信息的调制脉冲信号光;调制后的脉冲信号光被单点探测器探测并由高速数据采集模块第i次采集得到不同距离回波的总光强,第i次采集到的不同距离回波的总光强乘以第i个一维随机调制向量并求和,得到第i次测量两次调制后的所有距离回波的总光强yi,其中i=1,2,3,4…M;
(3)通过矩阵Φ1和Φ2得到测量矩阵Φ为M×(A×B)维的二维矩阵;
&Phi; = B 1 ( 1,1 ) &CenterDot; &Phi; 1 ( 1 ) B 1 ( 1,2 ) &CenterDot; &Phi; 1 ( 1 ) &CenterDot; &CenterDot; &CenterDot; B 1 ( 1 , r ) &CenterDot; &CenterDot; &Phi; 1 ( 1 ) B 2 ( 1,1 ) &CenterDot; &Phi; 1 ( 2 ) B 2 ( 1,2 ) &CenterDot; &Phi; 1 ( 2 ) &CenterDot; &CenterDot; &CenterDot; B 2 ( 1 , r ) &CenterDot; &Phi; 1 ( 2 ) &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; B M ( 1,1 ) &CenterDot; &CenterDot; &Phi; 1 ( M ) B M ( 1,2 ) &CenterDot; &Phi; 1 ( M ) &CenterDot; &CenterDot; &CenterDot; B M ( 1 , r ) &CenterDot; &Phi; 1 ( M )
(4)经过M次测量,M次测量总光强记为Y;
Y = y 1 y 2 &CenterDot; &CenterDot; &CenterDot; y M
其中,yi为第i次测量两次调制后的所有距离回波的总光强;
(5)最后根据Y=ΦX反算出目标物体的三维信息X。
最后对M次测量回波的总光强数据列Y进行压缩感知计算时满足关系如下:
Y = y 1 y 2 &CenterDot; &CenterDot; &CenterDot; y M = B 1 ( 1,1 ) &CenterDot; &Phi; 1 ( 1 ) B 1 ( 1,2 ) &CenterDot; &Phi; 1 ( 1 ) &CenterDot; &CenterDot; &CenterDot; B 1 ( 1 , r ) &CenterDot; &Phi; 1 ( 1 ) B 2 ( 1,1 ) &CenterDot; &Phi; 1 ( 2 ) B 2 ( 1,2 ) &CenterDot; &Phi; 1 ( 2 ) &CenterDot; &CenterDot; &CenterDot; B 2 ( 1 , r ) &CenterDot; &Phi; 1 ( 2 ) &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; B M ( 1,1 ) &CenterDot; &Phi; 1 ( M ) B M ( 1,2 ) &CenterDot; &Phi; 1 ( M ) &CenterDot; &CenterDot; &CenterDot; B M ( 1 , r ) &CenterDot; &Phi; 1 ( M ) x 1,1 x 1,2 &CenterDot; &CenterDot; &CenterDot; x 1 , A x 2,1 &CenterDot; &CenterDot; &CenterDot; x 2 , A &CenterDot; &CenterDot; &CenterDot; x B , 1 &CenterDot; &CenterDot; &CenterDot; x B , A = &Phi;X
其中,xj,i为第j个距离上二维图像的第i个像素;利用压缩感知算法得到每个距离上的二维信息,即得到目标景物的三维信息X。
与现有技术相比,本发明的有益效果体现在:
本发明的基于压缩感知脉冲编码调制的三维信息获取装置不需要对回波信号切片进行多次压缩感知计算,而是利用一维随机编码向量对采集模块采集到的回波信号进行编码调制,然后将调制后的不同距离的光强相加得到所有距离回波的总光强,利用压缩感知算法进行一次计算获得目标景物的三维信息,该系统采用单点探测器,三维信息获取过程测量时间较短。利用压缩感知算法进一步缩短了测量时间,因此利用本发明的方法可快速获得目标景物的三维信息。
附图说明
图1是本发明的空间光调制器采用反射型空间光调制器的基于压缩感知脉冲编码调制的三维信息获取装置的组成图。
图2是本发明的基于压缩感知脉冲编码调制的三维信息获取装置的同步时序图。
图3是空间光调制器加载的二维满足伯努利分布的0,1二值图像。
图4是本发明的空间光调制器采用透射型空间光调制器的基于压缩感知脉冲编码调制的三维信息获取装置的组成图。
图5是本发明基于压缩感知脉冲编码调制的三维信息获取方法的具体实现流程图。
具体实施方式
下面结合实施例和附图来详细说明本发明。
实施例1
如图1所示,一种基于压缩感知脉冲编码调制的三维信息获取装置,该装置中空间光调制器采用反射型空间光调制器,这里的反射型空间光调制器可以为反射型的硅上液晶器件(LCOS)或者反射型的数字微镜器件(DMD)等,在本实施例中优选为反射型的数字微镜器件(DMD)。
上述的三维信息获取装置包括:脉冲激光光源1、成像系统2、DMD空间光调制器3、会聚透镜4、单点探测器5、图像生成加载模块6、高速数据采集模块7、计算存储模块8、同步探测模块9和分光棱镜10。
脉冲激光光源1发出脉冲激光照射目标景物,成像系统2收集目标景物的漫反射光,DMD空间光调制器3位于成像系统的像面上,所以目标景物的像将成在DMD空间光调制器3上,DMD空间光调制器3前面一般还加有TIR棱镜以调整光路空间布置。
图像生成加载模块6用于生成一幅随机二维编码图像Φ1(i),并将其加载到DMD空间光调制器3上。由于DMD空间光调制器3上微镜阵元并不是平行于像面,而是与像面形成一定的夹角,所以由DMD空间光调制器3反射的光束并不是沿着入射光路射出。
随机二维编码图像Φ1(i)和目标景物的像相互作用后由DMD空间光调制器3反射,反射出的光束经会聚透镜4并由单点探测器5接收不同距离回波的总光强,高速数据采集模块6采集单点探测器探测的不同距离回波的总光强。
高速数据采集模块6对采集得到的不同距离回波的总光强进行一维随机编码调制,一维随机调制向量为B=1×r,r为一维随机调制向量元素数。高速数据采集模块6第i次采集到的数据乘以第i个一维随机调制向量Φ2(i)并求和,得到两次调制后的所有距离回波的总光强信息yi
重复上述的步骤,多次测量之后,利用图像生成加载模块6里存储的二维随机图像Φ1(i)和已知的一维随机调制向量Φ2(i)和高速数据采集模块7里得到两次调制后的所有距离回波的总光强信息yi,利用压缩感知算法进行一次计算,就可以得到目标景物的三维信息。
如图2所示,为实现数据的同步更新,图像生成加载模块6和脉冲激光光源1设有同步模块,实现加载一幅二维随机图像脉冲激光光源发出一个脉冲光,高速数据采集模块7和脉冲激光光源1间同样设有同步探测模块9,实现脉冲激光光源每发出一个脉冲光高速数据采集模块7对两路光进行同步采集。
如图3所示,二维随机调制图像Φ1(i)是满足伯努利分布的二维0,1随机图像,一维随机调制向量Φ2(j)是满足伯努利分布的一维0,1随机向量。
如图5所示,利用压缩感知计算目标景物的三维信息的过程为:
(1)设定二维随机调制图像总像素数A=p×q,p为二维随机调制图像的横向分辨率,q为二维随机调制图像的纵向分辨率;
一维随机编码向量B=1×r,目标景物的测量次数为M,生成一个M×A维的矩阵Φ1和一个M×B维的矩阵Φ2,利用矩阵Φ1生成M幅二维随机调制图像Φ1(i),利用矩阵Φ2生成M个一维随机调制向量Φ2(j),其中i,j=1,2,3,4…M;
&Phi; 1 = A 1 ( 1,1 ) A 1 ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; A 1 ( p , q ) A 2 ( 1,1 ) A 2 ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; A 2 ( p , q ) &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; A M ( 1,1 ) A M ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; A M ( p , q ) = &Phi; 1 ( 1 ) &Phi; 1 ( 2 ) &CenterDot; &CenterDot; &CenterDot; &Phi; 1 ( M )
&Phi; 2 = B 1 ( 1,1 ) B 1 ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; B 1 ( 1 , r ) B 2 ( 1,1 ) B 2 ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; B 2 ( 1 , r ) &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; B M ( 1,1 ) B M ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; B M ( 1 , r ) = &Phi; 2 ( 1 ) &Phi; 2 ( 2 ) &CenterDot; &CenterDot; &CenterDot; &Phi; 2 ( M )
其中,Ai(m,n)为第i幅二维随机调制图像中的图像元素,m=1,2,3…p,n=1,2,3…q;Bi(1,s)为第i个一维随机调制向量中的元素,s=1,2,3…r;
(2)空间光调制器依次在脉冲信号光中加载第i幅二维随机调制图像,得到带有第i幅二维随机调制图像信息的调制脉冲信号光;调制后的脉冲信号光被单点探测器探测并由高速数据采集模块第i次采集得到不同距离回波的总光强,第i次采集到的不同距离回波的总光强乘以第i个一维随机调制向量并求和,得到第i次测量两次调制后的所有距离回波的总光强yi,其中i=1,2,3,4…M;
(3)通过矩阵Φ1和Φ2得到测量矩阵Φ为M×(A×B)维的二维矩阵;
&Phi; = B 1 ( 1,1 ) &CenterDot; &Phi; 1 ( 1 ) B 1 ( 1,2 ) &CenterDot; &Phi; 1 ( 1 ) &CenterDot; &CenterDot; &CenterDot; B 1 ( 1 , r ) &CenterDot; &Phi; 1 ( 1 ) B 2 ( 1,1 ) &CenterDot; &Phi; 1 ( 2 ) B 2 ( 1,2 ) &CenterDot; &Phi; 1 ( 2 ) &CenterDot; &CenterDot; &CenterDot; B 2 ( 1 , r ) &CenterDot; &Phi; 1 ( 2 ) &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; B M ( 1,1 ) &CenterDot; &Phi; 1 ( M ) B M ( 1,2 ) &CenterDot; &Phi; 1 ( M ) &CenterDot; &CenterDot; &CenterDot; B M ( 1 , r ) &CenterDot; &Phi; 1 ( M )
(4)经过M次测量,M次测量总光强记为Y;
Y = y 1 y 2 &CenterDot; &CenterDot; &CenterDot; y M
其中,yi为两次调制后的所有距离回波的总光强;
(5)最后对M次测量总光强数据列Y进行压缩感知计算时满足关系如下:
Y = y 1 y 2 &CenterDot; &CenterDot; &CenterDot; y M = B 1 ( 1,1 ) &CenterDot; &Phi; 1 ( 1 ) B 1 ( 1,2 ) &CenterDot; &Phi; 1 ( 1 ) &CenterDot; &CenterDot; &CenterDot; B 1 ( 1 , r ) &CenterDot; &Phi; 1 ( 1 ) B 2 ( 1,1 ) &CenterDot; &Phi; 1 ( 2 ) B 2 ( 1,2 ) &CenterDot; &Phi; 1 ( 2 ) &CenterDot; &CenterDot; &CenterDot; B 2 ( 1 , r ) &CenterDot; &Phi; 1 ( 2 ) &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; B M ( 1,1 ) &CenterDot; &Phi; 1 ( M ) B M ( 1,2 ) &CenterDot; &Phi; 1 ( M ) &CenterDot; &CenterDot; &CenterDot; B M ( 1 , r ) &CenterDot; &Phi; 1 ( M ) x 1,1 x 1,2 &CenterDot; &CenterDot; &CenterDot; x 1 , A x 2,1 &CenterDot; &CenterDot; &CenterDot; x 2 , A &CenterDot; &CenterDot; &CenterDot; x B , 1 &CenterDot; &CenterDot; &CenterDot; x B , A = &Phi;X
xj,i为第j个距离上二维图像的第i个像素;利用压缩感知算法进行一次计算得到每个距离上的二维信息,即得到目标景物的三维信息X。
实施例2
如图4所示,是空间光调制器采用透射型空间光调制器的基于压缩感知的三维信息获取系统。这里透射型空间光调制器可以是透射型的硅上液晶器件(LCOS)等。图4中以透射型液晶器件为例,包括脉冲激光光源1、成像系统2、透射型空间光调制器11、会聚透镜4、单点探测器5、图像生成加载模块6、高速数据采集模块7、计算存储模块8、同步探测模块9和分光棱镜10。
本系统中,采用透射型空间光调制器11作为空间光调制器,脉冲激光光源1发出脉冲激光照射目标景物,成像系统2收集目标景物的漫反射光,透射型空间光调制器11位于成像系统的像面上,所以目标景物的像将成在透射型空间光调制器11上。由图像生成加载模块6生成一幅随机二维编码图像Φ1(i)并将其加载到透射型空间光调制器11上。随机二维编码图像Φ1(i)和目标景物的像相互作用后由透射型空间光调制器11透射,由会聚透镜4会聚并由单点探测器5接收不同距离回波的总光强,高速数据采集模块采集单点探测器探测的不同距离回波的总光强,采集卡第i次采集到的数据乘以第i个一维随机调制向量并求和,得到两次调制后的所有距离回波的总光强信息yi多次测量之后,利用图像生成加载模块6里存储的二维随机图像Φ1(i)和已知的一维随机调制向量Φ2(i)和高速数据采集模块7里得到两次调制后的所有距离回波的总光强信息yi,利用压缩感知算法进行一次计算,就可以得到目标景物的三维信息。

Claims (10)

1.一种基于压缩感知脉冲编码调制的三维信息获取装置,其特征在于,包括:脉冲激光模块、图像生成加载模块(6)、成像系统(2)、空间光调制器(3)、会聚透镜(4)、单点探测器(5)、高速数据采集模块(7)和计算存储模块(8);
所述脉冲激光模块发出脉冲激光并投射到目标景物上;
所述图像生成加载模块(6)生成多次测量加载的二维随机调制图像,并将二维随机调制图像传递给所述空间光调制器(3);
所述成像系统(2)用于收集所述目标景物的漫反射光,并将目标景物成像到空间光调制器(3)上;
所述空间光调制器(3)根据加载其上的二维随机调制图像对得到的目标景物的图像进行编码调制;
所述会聚透镜(4)将目标景物图像调制后的光强会聚得到调制后的不同距离回波的总光强;
所述单点探测器(5)用于探测所述会聚透镜(4)会聚得到不同距离回波的总光强;
所述高速数据采集模块(7)用于采集单点探测器(5)多次探测得到的不同距离回波的总光强;
所述计算存储模块(8)通过对多次探测得到的不同距离回波的总光强进行一维编码调制,利用压缩感知算法完成三维信息获取。
2.如权利要求1所述的基于压缩感知脉冲编码调制的三维信息获取装置,其特征在于,所述空间光调制器(3)为透射型空间光调制器或反射型空间光调制器。
3.如权利要求2所述的基于压缩感知脉冲编码调制的三维信息获取装置,其特征在于,所述的成像系统(2)为望远系统、显微系统或照相系统。
4.如权利要求3所述的基于压缩感知脉冲编码调制的三维信息获取装置,其特征在于,所述图像生成加载模块(6)内设有同步模块,实现每切换一张图像所述脉冲激光模块发一个脉冲激光。
5.如权利要求4所述的基于压缩感知脉冲编码调制的三维信息获取装置,其特征在于,所述高速数据采集模块(7)和脉冲激光模块之间连接有同步探测模块(9)。
6.如权利要求1所述的基于压缩感知脉冲编码调制的三维信息获取装置,其特征在于,所述脉冲激光模块包括脉冲激光光源(1)、扩束镜和分光棱镜(10)。
7.一种基于压缩感知脉冲编码调制的三维信息获取方法,其特征在于,包括以下步骤:
(1)脉冲激光光源发出脉冲激光照射目标景物;
(2)目标景物经成像系统成像到空间光调制器上,空间光调制器根据加载其上的二维随机调制图像对得到的目标景物的图像进行编码调制;
(3)将目标景物调制后图像的光强会聚后得到调制后的不同距离回波的总光强,并由单点探测器接收;
(4)利用高速数据采集模块采集单点探测器探测的不同距离回波的总光强;
(5)重复步骤(1)~(4)对目标景物进行M次测量,得到M次测量的回波总光强,最后利用压缩感知算法完成目标景物三维信息的获取。
8.如权利要求7所述的基于压缩感知脉冲编码调制的三维信息获取方法,其特征在于,在步骤(4)中,需要对高速数据采集模块每次采集的不同距离回波的总光强进行一维随机编码调制,其中,一维随机调制向量为B=1×r,r为一维随机调制向量元素数。
9.如权利要求8所述的基于压缩感知脉冲编码调制的三维信息获取方法,其特征在于,所述M为小于二维随机调制图像总像素数和一维随机调制向量元素数的乘积的自然数。
10.如权利要求9所述的基于压缩感知脉冲编码调制的三维信息获取方法,其特征在于,在步骤(5)中,利用压缩感知算法完成目标景物三维信息获取的具体方法为:
(1)设定二维随机调制图像总像素数A=p×q,p为二维随机调制图像的横向分辨率,q为二维随机调制图像的纵向分辨率;一维随机调制向量B=1×r,目标景物的测量次数为M,生成一个M×A维的矩阵Φ1和一个M×B维的矩阵Φ2,利用矩阵Φ1生成M幅二维随机调制图像Φ1(i),利用矩阵Φ2生成M个一维随机调制向量Φ2(j),其中i,j=1,2,3,4…M;
&Phi; 1 = A 1 ( 1,1 ) A 1 ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; A 1 ( p , q ) A 2 ( 1,1 ) A 2 ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; A 2 ( p , q ) &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; A M ( 1,1 ) A M ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; A M ( p , q ) = &Phi; 1 ( 1 ) &Phi; 1 ( 2 ) &CenterDot; &CenterDot; &CenterDot; &Phi; 1 ( M )
&Phi; 2 = B 1 ( 1,1 ) B 1 ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; B 1 ( 1 , r ) B 2 ( 1,1 ) B 2 ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; B 2 ( 1 , r ) &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; B M ( 1,1 ) B M ( 1,2 ) &CenterDot; &CenterDot; &CenterDot; B M ( 1 , r ) = &Phi; 2 ( 1 ) &Phi; 2 ( 2 ) &CenterDot; &CenterDot; &CenterDot; &Phi; 2 ( M )
其中,Ai(m,n)为第i幅二维随机调制图像中的图像元素,m=1,2,3…p,n=1,2,3…q;Bi(1,s)为第i个一维随机调制向量中的元素,s=1,2,3…r;
(2)空间光调制器依次在脉冲信号光中加载第i幅二维随机调制图像,得到带有第i幅二维随机调制图像信息的调制脉冲信号光;调制后的脉冲信号光被单点探测器探测并由高速数据采集模块第i次采集得到不同距离回波的总光强,第i次采集到的不同距离回波的总光强乘以第i个一维随机调制向量并求和,得到第i次测量两次调制后的所有距离回波的总光强yi,其中i=1,2,3,4…M;
(3)通过矩阵Φ1和Φ2得到测量矩阵Φ为M×(A×B)维的二维矩阵;
&Phi; = B 1 ( 1,1 ) &CenterDot; &Phi; 1 ( 1 ) B 1 ( 1,2 ) &CenterDot; &Phi; 1 ( 1 ) &CenterDot; &CenterDot; &CenterDot; B 1 ( 1 , r ) &CenterDot; &Phi; 1 ( 1 ) B 2 ( 1,1 ) &CenterDot; &Phi; 1 ( 2 ) B 2 ( 1,2 ) &CenterDot; &Phi; 1 ( 2 ) &CenterDot; &CenterDot; &CenterDot; B 2 ( 1 , r ) &CenterDot; &Phi; 1 ( 2 ) &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; B M ( 1,1 ) &CenterDot; &Phi; 1 ( M ) B M ( 1,2 ) &CenterDot; &Phi; 1 ( M ) &CenterDot; &CenterDot; &CenterDot; B M ( 1 , r ) &CenterDot; &Phi; 1 ( M )
(4)经过M次测量,M次测量总光强记为Y;
Y = y 1 y 2 &CenterDot; &CenterDot; &CenterDot; y M
其中,yi为第i次测量两次调制后的所有距离回波的总光强;
(5)最后根据Y=ΦX反算出目标物体的三维信息X。
CN201310125563.9A 2013-04-10 2013-04-10 基于压缩感知脉冲编码调制的三维信息获取装置和方法 Expired - Fee Related CN103234479B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310125563.9A CN103234479B (zh) 2013-04-10 2013-04-10 基于压缩感知脉冲编码调制的三维信息获取装置和方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310125563.9A CN103234479B (zh) 2013-04-10 2013-04-10 基于压缩感知脉冲编码调制的三维信息获取装置和方法

Publications (2)

Publication Number Publication Date
CN103234479A true CN103234479A (zh) 2013-08-07
CN103234479B CN103234479B (zh) 2015-07-08

Family

ID=48882531

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310125563.9A Expired - Fee Related CN103234479B (zh) 2013-04-10 2013-04-10 基于压缩感知脉冲编码调制的三维信息获取装置和方法

Country Status (1)

Country Link
CN (1) CN103234479B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103558604A (zh) * 2013-10-12 2014-02-05 浙江大学 飞行时间原理的调制型漫反射表面反射成像方法与系统
CN103884486A (zh) * 2014-02-27 2014-06-25 中国科学院力学研究所 纹影测量成像系统及方法
CN103913288A (zh) * 2014-03-14 2014-07-09 中国科学院力学研究所 彩虹纹影测量成像系统及方法
CN106019274A (zh) * 2016-05-24 2016-10-12 西安深穹光电科技有限公司 一种新型多普勒雷达成像装置及方法
CN106066481A (zh) * 2016-05-24 2016-11-02 西安深穹光电科技有限公司 一种基于dmd的四维激光雷达成像装置及方法
CN107944422A (zh) * 2017-12-08 2018-04-20 业成科技(成都)有限公司 三维摄像装置、三维摄像方法及人脸识别方法
CN108007385A (zh) * 2017-12-31 2018-05-08 中国人民解放军陆军工程大学 一种大视场弹坑表面形貌成像系统及方法
CN109253997A (zh) * 2018-10-27 2019-01-22 西安电子科技大学 一种基于频率调制和空间编码的拉曼断层成像系统
CN110425986A (zh) * 2019-07-17 2019-11-08 北京理工大学 基于单像素传感器的三维计算成像方法及装置
CN111123284A (zh) * 2019-12-26 2020-05-08 宁波飞芯电子科技有限公司 探测方法及探测设备

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2773729Y (zh) * 2005-03-25 2006-04-19 浙江大学 无扫描器脉冲调制式三维成像系统
CN102914367A (zh) * 2012-10-25 2013-02-06 浙江大学 基于压缩感知的多光谱成像装置和方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2773729Y (zh) * 2005-03-25 2006-04-19 浙江大学 无扫描器脉冲调制式三维成像系统
CN102914367A (zh) * 2012-10-25 2013-02-06 浙江大学 基于压缩感知的多光谱成像装置和方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张硕等: "基于压缩感知的三维物体成像的简单计算方法", 《光学学报》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103558604B (zh) * 2013-10-12 2015-07-15 浙江大学 飞行时间原理的调制型漫反射表面反射成像方法与系统
CN103558604A (zh) * 2013-10-12 2014-02-05 浙江大学 飞行时间原理的调制型漫反射表面反射成像方法与系统
CN103884486A (zh) * 2014-02-27 2014-06-25 中国科学院力学研究所 纹影测量成像系统及方法
CN103884486B (zh) * 2014-02-27 2017-01-11 中国科学院力学研究所 纹影测量成像系统及方法
CN103913288A (zh) * 2014-03-14 2014-07-09 中国科学院力学研究所 彩虹纹影测量成像系统及方法
CN106019274B (zh) * 2016-05-24 2019-01-22 西安深穹光电科技有限公司 一种新型多普勒雷达成像装置及方法
CN106019274A (zh) * 2016-05-24 2016-10-12 西安深穹光电科技有限公司 一种新型多普勒雷达成像装置及方法
CN106066481A (zh) * 2016-05-24 2016-11-02 西安深穹光电科技有限公司 一种基于dmd的四维激光雷达成像装置及方法
CN106066481B (zh) * 2016-05-24 2019-01-18 西安深穹光电科技有限公司 一种基于dmd的四维激光雷达成像装置及方法
CN107944422A (zh) * 2017-12-08 2018-04-20 业成科技(成都)有限公司 三维摄像装置、三维摄像方法及人脸识别方法
CN107944422B (zh) * 2017-12-08 2020-05-12 业成科技(成都)有限公司 三维摄像装置、三维摄像方法及人脸识别方法
CN108007385A (zh) * 2017-12-31 2018-05-08 中国人民解放军陆军工程大学 一种大视场弹坑表面形貌成像系统及方法
CN108007385B (zh) * 2017-12-31 2023-11-14 中国人民解放军陆军工程大学 一种大视场弹坑表面形貌成像系统及方法
CN109253997A (zh) * 2018-10-27 2019-01-22 西安电子科技大学 一种基于频率调制和空间编码的拉曼断层成像系统
CN110425986A (zh) * 2019-07-17 2019-11-08 北京理工大学 基于单像素传感器的三维计算成像方法及装置
CN110425986B (zh) * 2019-07-17 2020-10-16 北京理工大学 基于单像素传感器的三维计算成像方法及装置
CN111123284A (zh) * 2019-12-26 2020-05-08 宁波飞芯电子科技有限公司 探测方法及探测设备
WO2021129413A1 (zh) * 2019-12-26 2021-07-01 宁波飞芯电子科技有限公司 探测方法及探测设备
CN111123284B (zh) * 2019-12-26 2022-02-11 宁波飞芯电子科技有限公司 探测方法及探测设备

Also Published As

Publication number Publication date
CN103234479B (zh) 2015-07-08

Similar Documents

Publication Publication Date Title
CN103234479B (zh) 基于压缩感知脉冲编码调制的三维信息获取装置和方法
JP6771624B2 (ja) 線形モード計算センシングレーザーレーダー
US11166004B2 (en) Three-dimensional computational imaging method and apparatus based on single-pixel sensor, and non-transitory computer-readable storage medium
CN106772430B (zh) 基于多分辨率小波逼近的单像素光子计数三维成像系统及方法
CN103472457A (zh) 稀疏孔径压缩计算关联飞行时间的三维成像系统及方法
CN103472456B (zh) 一种基于稀疏孔径压缩计算关联的主动成像系统及方法
CN102914367B (zh) 基于压缩感知的多光谱成像装置和方法
Georgopoulos et al. Assessing the performance of a structured light scanner
CN103453993B (zh) 基于稀疏孔径压缩计算关联的主动高光谱成像系统及方法
CN102608619B (zh) 光子计数压缩采样相控阵激光三维成像方法
Heide et al. Diffuse mirrors: 3D reconstruction from diffuse indirect illumination using inexpensive time-of-flight sensors
CN102375144A (zh) 单光子计数压缩采样激光三维成像方法
CN102062861B (zh) 基于单探测器关联成像原理的三维成像方法
CN103471718A (zh) 一种基于稀疏孔径压缩计算关联的高光谱成像系统及方法
CN102788979B (zh) 一种基于后向投影InSAR成像配准的GPU实现方法
CN102865833B (zh) 基于等高信息稀疏测量的三维成像装置及方法
CN105676613A (zh) 一种采用单像素桶探测器的数字全息幻影成像系统及其工作方法
CN103616682B (zh) 一种基于曲面投影的多基线InSAR处理方法
Howland et al. Compressive sensing LIDAR for 3D imaging
CN104992424A (zh) 一种基于离散余弦变换的单像素快速主动成像系统
Breuckmann 25 Years of High Definition 3D Scanning: History, State of the Art, Outlook
CN101545976B (zh) 反射层析激光雷达双折线反投影成像方法
CN106791781A (zh) 一种连续波相位测量式单像素三维成像系统及方法
Mei et al. Phase unwrapping and denoising for time-of-flight imaging using generalized approximate message passing
Salin et al. Combined method for measuring 3D wave spectra. I. algorithms to transform the optical-brightness field into the wave-height distribution

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150708

Termination date: 20200410