CN107015466B - 基于tr-music算法的光学扫描全息单点定位方法 - Google Patents
基于tr-music算法的光学扫描全息单点定位方法 Download PDFInfo
- Publication number
- CN107015466B CN107015466B CN201710255761.5A CN201710255761A CN107015466B CN 107015466 B CN107015466 B CN 107015466B CN 201710255761 A CN201710255761 A CN 201710255761A CN 107015466 B CN107015466 B CN 107015466B
- Authority
- CN
- China
- Prior art keywords
- hologram
- indicate
- test target
- target
- matrix
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000003287 optical effect Effects 0.000 title claims abstract description 28
- 238000001093 holography Methods 0.000 title claims abstract description 19
- 238000003384 imaging method Methods 0.000 claims abstract description 35
- 239000011159 matrix material Substances 0.000 claims abstract description 31
- 238000001514 detection method Methods 0.000 claims abstract description 27
- 239000013598 vector Substances 0.000 claims abstract description 19
- 230000009466 transformation Effects 0.000 claims abstract description 8
- 238000001228 spectrum Methods 0.000 claims description 29
- 210000001747 pupil Anatomy 0.000 claims description 26
- 230000010287 polarization Effects 0.000 claims description 20
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 230000017105 transposition Effects 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 5
- 230000021615 conjugation Effects 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 230000035939 shock Effects 0.000 claims description 3
- 239000013604 expression vector Substances 0.000 claims description 2
- 230000010354 integration Effects 0.000 claims 1
- 238000004064 recycling Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 3
- 240000005589 Calophyllum inophyllum Species 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000005693 optoelectronics Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000005622 photoelectricity Effects 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 238000006116 polymerization reaction Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/22—Processes or apparatus for obtaining an optical image from holograms
- G03H1/2249—Holobject properties
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/22—Processes or apparatus for obtaining an optical image from holograms
- G03H1/2249—Holobject properties
- G03H2001/2252—Location of the holobject
- G03H2001/2255—Holobject out of Fourier or hologram planes
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Holo Graphy (AREA)
Abstract
本发明公开了基于TR‑MUSIC算法的光学扫描全息单点定位方法,解决了现有技术定位精准度不够,以及实现定位方式复杂的问题。本发明首先将激光转化成菲涅尔波带板后扫描物体,得到物体的全息图并将该全息图采用傅里叶变换得到矩阵K;然后通过矩阵K得到物体的时间反演矩阵,并求得其特征值和特征向量,从而将全息图分解为信号子空间和噪声子空间;最后利用有限元法将物体等间距离散化,并将各个单元作为测试目标后求得各个单元作为测试目标时的全息图,再利用全息图噪声子空间与信号子空间的正交性得到物体的成像伪谱,最终得到探测目标的位置信息。本发明不仅实现方式简单、便于操作,同时具有很强的实用性,适合推广使用。
Description
技术领域
本发明涉及基于TR-MUSIC算法的光学扫描全息单点定位方法。
背景技术
光学扫描全息技术,简称OSH,是数字全息技术的一个重要分支。它利用光学扫描技术将物体的3维信息切片储存为2维信息,从而得到物体的全息图。该技术是1979年,Poon和Korpel在研究声光外差图像处理器的时候提出。自该技术提出以来,已经在扫描全息显微镜、3D图像识别以及3D光学遥感等领域得到了广泛的应用。
目前光学扫描全息技术实现起来比较复杂,并且在定位的精准上也有所欠缺。
发明内容
本发明要解决的技术问题是:提供基于TR-MUSIC算法的光学扫描全息单点定位方法,实现对单个物体的精准定位,并同时提出了将TR-MUSIC应用到OSH***中以实现对物体定位的相关理论。
为实现上述目的,本发明采用的技术方案如下:
基于TR-MUSIC算法的光学扫描全息单点定位方法,包括以下步骤:
步骤1、将激光采用第一偏振分束器分成两束,之后第一束光依次通过第一光瞳和第一凸透镜后投射至第二偏振分束器,第二束光依次通过第二光瞳和第二凸透镜后投射至第二偏振分束器,第二偏振分束器将投射来的两束光聚光干涉形成菲涅尔波带板;
步骤2、首先采用步骤1中获得的菲涅尔波带板对物体进行扫描,然后利用光电探测器接受扫描后的透射光,经过解调后得到物体的切片全息图,最后将全息图采用傅里叶变换后得到矩阵K;
步骤3、通过步骤2中K矩阵得到物体的时间反演矩阵,并求得其特征值和特征向量,从而将全息图分解为信号子空间和噪声子空间;
步骤4、首先利用有限元法将物体等间距离散化,并将各个单元作为测试目标,然后求得各个单元作为测试目标时的全息图,最后利用全息图噪声子空间与信号子空间的正交性得到物体的成像伪谱,最终得到探测目标的位置信息。
进一步地,步骤1中第一光瞳的函数为矩形1函数,第二光瞳的函数为狄拉克δ函数,第一凸透镜和第二凸透镜的焦距相同。
具体地说,步骤1中激光的光学传递函数如下:
将p1(x,y)=1和p2(x,y)=δ(x,y)代入式(1)中,则式(1)表示为下式:
则式(2)相应的空间冲击响应为:
其中,j表示虚数单位*表示卷积运算,x'和y'分别表示横向和纵向的积分变量,x表示物体的横向坐标,y表示物体的纵向坐标,z表示2D扫描镜到待测物体的距离,表示波数,λ表示光波波长,f表示凸透镜的焦距,kx和ky表示频域坐标,p1(x,y)和p2(x,y)分别表示第一光瞳和第二光瞳函数。
具体地说,步骤2中所述得到物体切片全息图的关系式如下:
其中,|Γ(x,y;z)|2表示物体的复振幅函数,h(x,y;z0)表示扫描位置在z0的点扩散函数,F和F-1分别表示傅里叶变换和傅里叶反变换,*表示卷积运算。
具体地说,步骤2中所述全息图采用傅里叶变换得到矩阵K的公式如下:
K=F{Hc(x,y;z0)}=F{|Γ(x,y;z0)|2*h(x,y;z0)} (5)
其中,其中x表示物体的横向坐标,y表示物体的纵向坐标,z0表示扫描镜到物体的距离并且此处z0为单值,Hc(x,y;z0)表示物体的全息图。
具体地说,步骤3的实现方法如下:
首先将式(5)中求得的矩阵K做奇异值分解即可得时间反演矩阵KHK和KKH,具体如下:
KHK=F-1{F{(|Γ(x,y;z0)|2)H*h(x,y;-z0)}·F{h(x,y;z0)*(|Γ(x,y;z0)|2)}}=F-1{F{(|Γ(x,y;z0)|2)H*h(x,y;-z0)*h(x,y;z0)*(|Γ(x,y;z0)|2)}}
其中,H表示矩阵的共轭转置运算,[h(x,y;z0)]H=h(x,y;-z0),并且OSH中全息图的重建过程为:
|Γ(x,y;z0)|2=(|Γ(x,y;z0)|2)*h(x,y;z0)*h(x,y;-z0) (7)
因此,得到时间反演矩阵KHK和KKH如下:
KHK=(|Γ(x,y;z0)|2)H*|Γ(x,y;z0)|2
KKH=|Γ(x,y;z0)|2*(|Γ(x,y;z0)|2)H (8)
然后,求得时间反演阵KHK和KKH的特征值λ和特征向量v1、v2,非零特征值对应的特征向量对应信号子空间,其余的对应噪声子空间;由式(8)可知,v1携带了物体y方向的特征信息,v2携带了物体x方向的特征信息,其中信号子空间与噪声子空间相互正交,具体如下:
<v1(i=1,...M),v1(j=M+1,...N)>=0
<v2(i=1,...M),v2(j=M+1,...N)>=0
其中,M≤N且M=1,M表示非零特征值的个数,N表示特征值的总个数。
具体地说,步骤4中利用全息图噪声子空间与信号子空间的正交性得到物体成像伪谱的方法如下:
用v2中噪声子空间对应的特征向量的共轭转置左乘测试目标的全息图Ki(Xp)并求和得Qx(Xp),用v1中噪声子空间对应的特征向量右乘测试目标的全息图Ki(Xp)并求和Qy(Xp),即:
其中,Xp为各个单元测试目标,并且Xp为N×N个点,Ki(Xp)为各个单元测试目标Xp通过式(4)求得的全息图,i=1,2,3...N2,m表示向量|Ki(Xp)v1(j)|2或|v2(j)Ki(Xp)|2第m个元素,j表示v1(j)或v2(j)的第j个特征向量,M≤N且M=1,M表示非零特征值的个数,N表示特征值的总个数,Qx(Xp)和Qy(Xp)表示x和y方向测试目标Xp伪谱的一个中间值。
当测试目标在探测目标的位置时,Qx(Xp)和Qy(Xp)等于0,而测试目标不在探测目标位置时,Qx(Xp)和Qy(Xp)的值有限,在求得测试目标Xp的Qx(Xp)和Qy(Xp)之后,即可求得该点的x方向和y方向的成像伪谱Px(Xp)和Py(Xp),具体为:
Px(Xp)=||Ki(Xp)||2/Qx(Xp)
Py(Xp)=||Ki(Xp)||2/Qy(Xp) (10)
其中,Px(Xp)和Py(Xp)分别表示探测目标在x和y方向上的成像伪谱,式(10)中的Px(Xp)和Py(Xp)相乘得到测试目标Xp的成像伪谱P(Xp),具体为:
P(Xp)=Px(Xp)Py(Xp) (11)
当逐点求得P(Xp)之后,得到包含探测目标位置信息的成像伪谱,最终通过成像伪谱得到探测目标的精准位置。
具体地说,步骤4中利用全息图噪声子空间与信号子空间的正交性得到物体成像伪谱的方法如下:
首先,对测试目标Xp的全息图Ki(Xp)做奇异值分解,即求得Ki(Xp)时间反演阵[Ki(Xp)]HKi(Xp)及Ki(Xp)[Ki(Xp)]H的非零特征值对应的特征向量vy和vx;
然后,通过下式求得测试目标Xp的Qx(Xp)和Qy(Xp),具体为:
其中,*表示共轭,T表示矩阵的转置运算,当测试目标在探测目标的位置时,Qx(Xp)和Qy(Xp)等于0,而测试目标不在探测目标位置时,Qx(Xp)和Qy(Xp)的值有限;
最后,通过下式就可求得测试目标Xp的成像伪谱P(Xp),最终通过成像伪谱得到探测目标的精准位置,具体为:
Px(Xp)=||vx||2/Qx(Xp)
Py(Xp)=||vy||2/Qy(Xp)
P(Xp)=Px(Xp)Py(Xp) (13)
其中Px(Xp)和Py(Xp)分别表示探测目标在x和y方向上的成像伪谱,P(Xp)表示包含探测目标位置信息的成像伪谱。
与现有技术相比,本发明具有以下有益效果:
(1)本发明首创性的使用TR-MUSIC算法对全息图进行处理,从而获得物体中探测目标的准确位置。
(2)本发明使用TR-MUSIC算法实现对探测目标的定位,具有较高的精度与定位的准确性。
(3)本发明针对TR-MUSIC算法,利用物体噪声子空间与信号子空间的正交性得到物体的成像伪谱,给出了两种可行方案,并说明了两种方法的关系。给出了在全息扫描技术下使用TR-MUSIC算法实现定位的相关公式,即式(1)-式(13),建立了本发明的理论基础。
(4)本发明不仅实现方式简单、便于操作,同时具有很强的实用性,适合推广使用。
附图说明
图1为本发明流程示意图。
图2为本发明实施例采用的基本结构图。
图3为本发明实施例采用的探测目标位置示意图。
图4为本发明实施例的奇异值分解后,其特征值分布图。
图5为本发明实施例最终得到的成像伪谱。
具体实施方式
下面结合附图说明和实施例对本发明作进一步说明,本发明的方式包括但不仅限于以下实施例。
实施例
如图1,本发明提供的基于TR-MUSIC算法的光学扫描全息单点定位方法,采用TR-MUSIC算法对全息图进行处理,从而实现对探测目标的定位,具有较高的精度并且定位相当准确;本基于TR-MUSIC算法的光学扫描全息单点定位方法不仅实现方式简单、便于操作,同时具有很强的实用性,适合推广使用。
本基于TR-MUSIC算法的光学扫描全息单点定位方法采用如图2所示的结构来实现,该结构具体包括第一偏振分束器BS1、第二偏振分束器BS2、第一光瞳p1(x,y)、第二光瞳p2(x,y)、第一凸透镜L1、第二凸透镜L2、第三凸透镜L3、声光调制器AOFS、第一反射镜M1、第二反射镜M2、扫描器X-Y Scanner、光电二极管PD、光电振荡器、两个LPF、以及PC机。激光光源发射的激光通过第一偏振分束器BS1后分成两束光,其中一束光依次通过第一反射镜、第一光瞳和第一凸透镜射入第二偏振分束器BS2,另一束光依次通过声光调制器、第二反射镜、第二光瞳和第二凸透镜射入第二偏振分束器BS2,两束光经过第二偏振分束器BS2聚合后干涉形成菲涅尔波带板并投入扫描器,扫描器将上述菲涅尔波带板对被测物体(Object)进行扫描,扫描后的透射光通过第三透镜后经过光电二极管转化成电信号并经过光电振荡器后分别经过两个LPF传输至PC机。
本实施例所采用的结构中He-Ne激光器Laser的波长λ=632.8nm,两个凸透镜(L1、L2)的焦距都为400mm,扫描镜到物体的距离为z=400mm,物体如图3所示放置,物体采用的切片尺寸为8mm×8mm,物体经有限元法离散后的阵列为16×16,探测目标放在坐标为(10,5)处。
本基于TR-MUSIC算法的光学扫描全息单点定位方法的具体步骤为:
步骤1、得到菲涅尔波带板。
如图2所示,由同一激光光源发出的角频率为ω的光被第一偏振分束器BS1分成两束,其中一束通过第一光瞳p1(x,y)形成平面波;另一束经过声光调制器产生Ω的频移后再通过第二光瞳p2(x,y)形成球面波;两束光经第二偏振分束器BS2聚合后在待测物体前干涉形成菲涅尔波带板,其可以由式(3)得到,具体为:
首先,激光光源发出的激光经过第一偏振分束器BS1分束,再分别经过第一光瞳p1(x,y)和第二光瞳p2(x,y)后通过第二偏振分束器聚合后在待测物体前干涉形成菲涅尔波带板的光学传递函数为:
本实施例中第一光瞳的函数为矩形1函数,第二光瞳的函数为狄拉克δ函数,并且在第一光瞳和第二偏振分束器之间设有第一凸透镜,在第二光瞳和第二偏振分束器之间设有第二凸透镜,且第一凸透镜和第二凸透镜的焦距相同,于是,将p1(x,y)=1和p2(x,y)=δ(x,y)代入式(1)中,则式(1)表示如下:
则式(2)相应的空间冲击响应即为下式:
其中,x表示物体的横向坐标,y表示物体的纵向坐标,z表示2D扫描镜到待测物体的距离,表示波数,λ表示光波波长,f表示两个凸透镜的焦距,kx和ky表示频域坐标,p1(x,y)和p2(x,y)分别表示第一光瞳和第二光瞳函数。
步骤2、用菲涅尔波带板扫描物体,得到物体的全息图。
首先采用步骤1中获得的菲涅尔波带板对物体进行扫描,然后利用光电探测器接受扫描后的透射光,经过解调后得到物体的切片全息图,最后将全息图采用傅里叶变换后得到矩阵K
为得到物体的全息图,可使用:
又根据物体的全息图Hc(x,y;z0)得到矩阵K,
K=F{Hc(x,y;z0)}=F{|Γ(x,y;z0)|2*h(x,y;z0)} (5)
其中,x表示物体的横向坐标,y表示物体的纵向坐标,z0表示扫描镜到物体的距离(此处z0为单值),Hc(x,y;z0)表示物体的全息图。
步骤3、通过矩阵K得到物体的时间反演矩阵,并求得其特征值和特征向量,从而将全息图分解为信号子空间和噪声子空间。
首先,对矩阵K做奇异值分解,根据式(5)可以求得矩阵K,进而求得时间反演矩阵KHK和KKH,具体为:
KHK=F-1{F{(|Γ(x,y;z0)|2)H*h(x,y;-z0)}·F{h(x,y;z0)*(|Γ(x,y;z0)|2)}}=F-1{F{(|Γ(x,y;z0)|2)H*h(x,y;-z0)*h(x,y;z0)*(|Γ(x,y;z0)|2)}}
其中,H表示矩阵的共轭转置运算,[h(x,y;z0)]H=h(x,y;-z0),并且OSH中全息图的重建过程为:
|Γ(x,y;z0)|2=(|Γ(x,y;z0)|2)*h(x,y;z0)*h(x,y;-z0) (7)
因此,可以得到时间反演矩阵KHK和KKH如下:
KHK=(|Γ(x,y;z0)|2)H*|Γ(x,y;z0)|2
KKH=|Γ(x,y;z0)|2*(|Γ(x,y;z0)|2)H (8)
然后,求得时间反演阵KHK和KKH的特征值λ和特征向量v1、v2,非零特征值对应的特征向量对应信号子空间,其余的对应噪声子空间,如图4所示;且由式(8)可知,v1携带了物体y方向的特征信息,v2携带了物体x方向的特征信息。其中信号子空间与噪声子空间相互正交:
<v1(i=1,...M),v1(j=M+1,...N)>=0
<v2(i=1,...M),v2(j=M+1,...N)>=0
其中,M(M≤N且此处M=1)表示非零特征值的个数,也即是探测目标的个数,N表示特征值的总个数。
步骤4、利用有限元法将物体等间距离散化,并将各个单元作为测试目标;然后求得各个单元作为测试目标Xp(Xp为16×16个点),并通过式(4)可以求得测试目标Xp的全息图Ki(Xp)(i=1,2,3...N2)。
最后,利用物体噪声子空间与信号子空间的正交性得到物体的成像伪谱,如图5所示,并与图3放置的物品位置相吻合。为了得到物体的成像伪谱,本实施例提供了两种相关联的方法,具体如下:
方法1:用v2中噪声子空间对应的特征向量的共轭转置左乘Ki(Xp)并求和得Qx(Xp),用v1中噪声子空间对应的特征向量右乘Ki(Xp)并求和Qy(Xp),即如下式所示:
当测试目标在探测目标的位置时,Qx(Xp)和Qy(Xp)等于0,而测试目标不在探测目标位置时,Qx(Xp)和Qy(Xp)的值有限。求得测试目标Xp的Qx(Xp)和Qy(Xp)之后,即可求得该点的x方向和y方向的成像伪谱Px(Xp)和Py(Xp)为:
Px(Xp)=||Ki(Xp)||2/Qx(Xp)
Py(Xp)=||Ki(Xp)||2/Qy(Xp) (10)
其中Px(Xp)和Py(Xp)分别表示探测目标在x和y方向上的成像伪谱,式(10)中的Px(Xp)和Py(Xp)相乘可以得到测试目标Xp的成像伪谱:
P(Xp)=Px(Xp)Py(Xp) (11)
当逐点求得P(Xp)之后,可以得到包含探测目标位置信息的成像伪谱,最终,通过成像伪谱得到探测目标的位置。
方法2:首先,对测试目标Xp的全息图Ki(Xp)做奇异值分解,即求得Ki(Xp)时间反演阵[Ki(Xp)]HKi(Xp)及Ki(Xp)[Ki(Xp)]H的非零特征值对应的特征向量vy和vx。然后通过式(12)求得测试目标Xp的Qx和Qy,即:
其中,*表示共轭。当测试目标在探测目标的位置时,Qx(Xp)和Qy(Xp)等于0,而测试目标不在探测目标位置时,Qx(Xp)和Qy(Xp)的值有限。最后,通过下式即可求得测试目标Xp的成像伪谱P(Xp),即:
Px(Xp)=||vx||2/Qx(Xp)
Py(Xp)=||vy||2/Qy(Xp)
P(Xp)=Px(Xp)Py(Xp) (13)
其中Px(Xp)和Py(Xp)分别表示探测目标在x和y方向上的成像伪谱,P(Xp)表示包含探测目标位置信息的成像伪谱。
值得注意的是,在步骤4中,利用物体噪声子空间与信号子空间的正交性得到物体的成像伪谱的两种方法是有关联性的,具体分析如下:
对于步骤3和步骤4的过程可以看作是对矩阵K做奇异值分解的过程,因此,有从而有而对应的而对于<v2(j=M+1,...,N),vx>≈0,<v1(j=M+1,...,N),vy>≈0。因此
其中C表示常数,y方向同理。从式(14)知道方法1和方法2之间相差一个常数的倍数关系,而这并不影响对探测目标位置的确定。
本发明首创性的使用TR-MUSIC算法对全息图进行处理,从而实现对探测目标的定位,具有较高的精度并且定位相当准确;同时,本发明针对TR-MUSIC算法,利用物体噪声子空间与信号子空间的正交性得到物体的成像伪谱,给出了两种可行方案,并说明了两种方法的关系。并且首创性的给出了在全息扫描技术下使用TR-MUSIC算法实现定位的相关公式,即式(1)-式(13),建立了本发明的理论基础,具备凸出的实质性特点和显著进步。本发明不仅实现方式简单、操作便捷,同时具有很强的实用性,适合推广使用。
上述实施例仅为本发明的优选实施方式之一,不应当用于限制本发明的保护范围,但凡在本发明的主体设计思想和精神上作出的毫无实质意义的改动或润色,其所解决的技术问题仍然与本发明一致的,均应当包含在本发明的保护范围之内。
Claims (8)
1.基于TR-MUSIC算法的光学扫描全息单点定位方法,其特征在于,包括以下步骤:
步骤1、将激光采用第一偏振分束器分成两束,之后第一束光依次通过第一光瞳和第一凸透镜后投射至第二偏振分束器,第二束光依次通过第二光瞳和第二凸透镜后投射至第二偏振分束器,第二偏振分束器将投射来的两束光聚光干涉形成菲涅尔波带板;
步骤2、首先采用步骤1中获得的菲涅尔波带板对物体进行扫描,然后利用光电探测器接受扫描后的透射光,经过解调后得到物体的切片全息图,最后将全息图采用傅里叶变换后得到矩阵K;
步骤3、通过步骤2中K矩阵得到物体的时间反演矩阵,并求得其特征值和特征向量,从而将全息图分解为信号子空间和噪声子空间;
步骤4、首先利用有限元法将物体等间距离散化,并将各个单元作为测试目标,然后求得各个单元作为测试目标时的全息图,最后利用全息图噪声子空间与信号子空间的正交性得到物体的成像伪谱,最终得到探测目标的位置信息。
2.根据权利要求1所述基于TR-MUSIC算法的光学扫描全息单点定位方法,其特征在于,步骤1中第一光瞳的函数为矩形1函数,第二光瞳的函数为狄拉克δ函数,第一凸透镜和第二凸透镜的焦距相同。
3.根据权利要求2所述基于TR-MUSIC算法的光学扫描全息单点定位方法,其特征在于,步骤1中激光的光学传递函数如下:
将p1(x,y)=1和p2(x,y)=δ(x,y)代入式(1)中,则式(1)表示为下式:
则式(2)相应的空间冲击响应为:
其中,j表示虚数单位表示卷积运算,x'和y'分别表示横向和纵向的积分变量,x表示物体的横向坐标,y表示物体的纵向坐标,z表示2D扫描镜到待测物体的距离,表示波数,λ表示光波波长,f表示凸透镜的焦距,kx和ky表示频域坐标,p1(x,y)和p2(x,y)分别表示第一光瞳和第二光瞳函数。
4.根据权利要求3所述基于TR-MUSIC算法的光学扫描全息单点定位方法,其特征在于,步骤2中所述得到物体切片全息图的关系式如下:
其中,|Γ(x,y;z)|2表示物体的复振幅函数,h(x,y;z0)表示扫描位置在z0的点扩散函数,F和F-1分别表示傅里叶变换和傅里叶反变换,*表示卷积运算。
5.根据权利要求4所述基于TR-MUSIC算法的光学扫描全息单点定位方法,其特征在于,步骤2中所述全息图采用傅里叶变换得到矩阵K的公式如下:
K=F{Hc(x,y;z0)}=F{|Γ(x,y;z0)|2*h(x,y;z0)} (5)
其中,其中x表示物体的横向坐标,y表示物体的纵向坐标,z0表示扫描镜到物体的距离并且此处z0为单值,Hc(x,y;z0)表示物体的全息图。
6.根据权利要求5所述基于TR-MUSIC算法的光学扫描全息单点定位方法,其特征在于,步骤3的实现方法如下:
首先将式(5)中求得的矩阵K做奇异值分解即可得时间反演矩阵KHK和KKH,具体如下:
KHK=F-1{F{(|Γ(x,y;z0)|2)H*h(x,y;-z0)}·F{h(x,y;z0)*(|Γ(x,y;z0)|2)}}
=F-1{F{(|Γ(x,y;z0)|2)H*h(x,y;-z0)*h(x,y;z0)*(|Γ(x,y;z0)|2)}}
其中,H表示矩阵的共轭转置运算,[h(x,y;z0)]H=h(x,y;-z0),并且OSH中全息图的重建过程为:
|Γ(x,y;z0)|2=(|Γ(x,y;z0)|2)*h(x,y;z0)*h(x,y;-z0) (7)
因此,得到时间反演矩阵KHK和KKH如下:
KHK=(|Γ(x,y;z0)|2)H*|Γ(x,y;z0)|2
KKH=|Γ(x,y;z0)|2*(|Γ(x,y;z0)|2)H (8)
然后,求得时间反演阵KHK和KKH的特征值λ和特征向量v1、v2,非零特征值对应的特征向量对应信号子空间,其余的对应噪声子空间;由式(8)可知,v1携带了物体y方向的特征信息,v2携带了物体x方向的特征信息,其中信号子空间与噪声子空间相互正交,具体如下:
<v1(i=1,...M),v1(j=M+1,...N)>=0
<v2(i=1,...M),v2(j=M+1,...N)>=0
其中,M≤N且M=1,M表示非零特征值的个数,N表示特征值的总个数。
7.根据权利要求6所述基于TR-MUSIC算法的光学扫描全息单点定位方法,其特征在于,步骤4中利用全息图噪声子空间与信号子空间的正交性得到物体成像伪谱的方法如下:
用v2中噪声子空间对应的特征向量的共轭转置左乘测试目标的全息图Ki(Xp)并求和得Qx(Xp),用v1中噪声子空间对应的特征向量右乘测试目标的全息图Ki(Xp)并求和Qy(Xp),即:
其中,Xp为各个单元测试目标,并且Xp为N×N个点,Ki(Xp)为各个单元测试目标Xp通过式(4)求得的全息图,i=1,2,3...N2,m表示向量|Ki(Xp)v1(j)|2或|v2(j)Ki(Xp)|2第m个元素,j表示v1(j)或v2(j)的第j个特征向量,M≤N且M=1,M表示非零特征值的个数,N表示特征值的总个数,Qx(Xp)和Qy(Xp)表示x和y方向测试目标Xp伪谱的一个中间值;
当测试目标在探测目标的位置时,Qx(Xp)和Qy(Xp)等于0,而测试目标不在探测目标位置时,Qx(Xp)和Qy(Xp)的值有限,在求得测试目标Xp的Qx(Xp)和Qy(Xp)之后,即可求得该点的x方向和y方向的成像伪谱Px(Xp)和Py(Xp),具体为:
Px(Xp)=||Ki(Xp)||2/Qx(Xp)
Py(Xp)=||Ki(Xp)||2/Qy(Xp) (10)
其中,Px(Xp)和Py(Xp)分别表示探测目标在x和y方向上的成像伪谱,式(10)中的Px(Xp)和Py(Xp)相乘得到测试目标Xp的成像伪谱P(Xp),具体为:
P(Xp)=Px(Xp)Py(Xp) (11)
当逐点求得P(Xp)之后,得到包含探测目标位置信息的成像伪谱,最终通过成像伪谱得到探测目标的精准位置。
8.根据权利要求6所述的基于TR-MUSIC算法的光学扫描全息单点定位方法,其特征在于,步骤4中利用全息图噪声子空间与信号子空间的正交性得到物体成像伪谱的方法如下:
首先,对测试目标Xp的全息图Ki(Xp)做奇异值分解,即求得Ki(Xp)时间反演阵[Ki(Xp)]HKi(Xp)及Ki(Xp)[Ki(Xp)]H的非零特征值对应的特征向量vy和vx;
然后,通过下式求得测试目标Xp的Qx(Xp)和Qy(Xp),具体为:
其中,*表示共轭,T表示矩阵的转置运算,当测试目标在探测目标的位置时,Qx(Xp)和Qy(Xp)等于0,而测试目标不在探测目标位置时,Qx(Xp)和Qy(Xp)的值有限;
最后,通过下式就可求得测试目标Xp的成像伪谱P(Xp),最终通过成像伪谱得到探测目标的精准位置,具体为:
Px(Xp)=||vx||2/Qx(Xp)
Py(Xp)=||vy||2/Qy(Xp)
P(Xp)=Px(Xp)Py(Xp) (13)
其中Px(Xp)和Py(Xp)分别表示探测目标在x和y方向上的成像伪谱,P(Xp)表示包含探测目标位置信息的成像伪谱。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710255761.5A CN107015466B (zh) | 2017-04-19 | 2017-04-19 | 基于tr-music算法的光学扫描全息单点定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710255761.5A CN107015466B (zh) | 2017-04-19 | 2017-04-19 | 基于tr-music算法的光学扫描全息单点定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107015466A CN107015466A (zh) | 2017-08-04 |
CN107015466B true CN107015466B (zh) | 2019-07-19 |
Family
ID=59447321
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710255761.5A Expired - Fee Related CN107015466B (zh) | 2017-04-19 | 2017-04-19 | 基于tr-music算法的光学扫描全息单点定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107015466B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5559603A (en) * | 1992-12-18 | 1996-09-24 | Kabushiki Kaisha Komatsu Seisakusho | three-dimensional image measuring device |
CN101743519A (zh) * | 2007-05-16 | 2010-06-16 | 视瑞尔技术公司 | 全息显示装置 |
JP2011033759A (ja) * | 2009-07-31 | 2011-02-17 | Sony Corp | 三次元画像撮像装置、三次元画像撮像方法、三次元画像情報取得システム、三次元画像情報取得方法 |
CN102087503A (zh) * | 2011-01-11 | 2011-06-08 | 浙江师范大学 | 一种双随机相位光学彩色图像加密装置及方法 |
CN102805613A (zh) * | 2012-08-13 | 2012-12-05 | 电子科技大学 | 一种基于两次扫描的高分辨率光学扫描全息切片成像方法 |
CN102920438A (zh) * | 2012-10-30 | 2013-02-13 | 电子科技大学 | 一种基于可变光瞳的高分辨率光学扫描全息切片成像方法 |
CN104614970A (zh) * | 2015-02-15 | 2015-05-13 | 电子科技大学 | 一种基于双孔光瞳的光学扫描全息图像边缘提取方法 |
CN105204311A (zh) * | 2015-07-06 | 2015-12-30 | 电子科技大学 | 一种基于高斯切趾的光学扫描全息边缘检测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8896840B2 (en) * | 2012-04-25 | 2014-11-25 | Canon Kabushiki Kaisha | Interferometric method and digital holographic microscope |
-
2017
- 2017-04-19 CN CN201710255761.5A patent/CN107015466B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5559603A (en) * | 1992-12-18 | 1996-09-24 | Kabushiki Kaisha Komatsu Seisakusho | three-dimensional image measuring device |
CN101743519A (zh) * | 2007-05-16 | 2010-06-16 | 视瑞尔技术公司 | 全息显示装置 |
JP2011033759A (ja) * | 2009-07-31 | 2011-02-17 | Sony Corp | 三次元画像撮像装置、三次元画像撮像方法、三次元画像情報取得システム、三次元画像情報取得方法 |
CN102087503A (zh) * | 2011-01-11 | 2011-06-08 | 浙江师范大学 | 一种双随机相位光学彩色图像加密装置及方法 |
CN102805613A (zh) * | 2012-08-13 | 2012-12-05 | 电子科技大学 | 一种基于两次扫描的高分辨率光学扫描全息切片成像方法 |
CN102920438A (zh) * | 2012-10-30 | 2013-02-13 | 电子科技大学 | 一种基于可变光瞳的高分辨率光学扫描全息切片成像方法 |
CN104614970A (zh) * | 2015-02-15 | 2015-05-13 | 电子科技大学 | 一种基于双孔光瞳的光学扫描全息图像边缘提取方法 |
CN105204311A (zh) * | 2015-07-06 | 2015-12-30 | 电子科技大学 | 一种基于高斯切趾的光学扫描全息边缘检测方法 |
Non-Patent Citations (1)
Title |
---|
光学全息扫描***中的自聚焦技术研究;潘慧宇;《中国优秀硕士学位论文全文数据库 信息科技辑》;20170228;I138-3026-正文第7-46页 |
Also Published As
Publication number | Publication date |
---|---|
CN107015466A (zh) | 2017-08-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102087411A (zh) | 量子成像方法及量子成像*** | |
CN108140248A (zh) | 用于测向和测频的分布式阵列 | |
CN102054272B (zh) | 一种基于联合分数傅里叶变换的位移矢量探测方法及装置 | |
CN105467806B (zh) | 单像素全息相机 | |
CN105116705B (zh) | 多波长自适应数字全息成像***及方法 | |
Neri et al. | Low-speed cameras system for 3D-DIC vibration measurements in the kHz range | |
Tiziani et al. | From speckle pattern photography to digital holographic interferometry | |
Kadambi et al. | Rethinking machine vision time of flight with ghz heterodyning | |
CN109000781A (zh) | 一种结构微振动线域测量装置及方法 | |
Jiao et al. | Enhanced autofocusing in optical scanning holography based on hologram decomposition | |
EP2313871B1 (en) | An electromagnetic body scanning system, a method and a computer program product | |
Li et al. | Exploiting wavelength diversity for high resolution time-of-flight 3D imaging | |
Pan et al. | Micro-Doppler effect based vibrating object imaging of coherent detection GISC lidar | |
Fong et al. | Phase estimation via riesz transform in laser speckle interferometry for large-area damage imaging | |
Chen et al. | Enhancing speckle statistics for imaging inside scattering media | |
EP3499201B1 (en) | Electromagnetic wave phase/amplitude generation device, electromagnetic wave phase/amplitude generation method, and electromagnetic wave phase/amplitude generation program | |
US11867625B2 (en) | System and method for imaging via scattering medium | |
CN109283821A (zh) | 基于涡旋透镜的相移数字全息单次曝光成像装置及方法 | |
CN107015466B (zh) | 基于tr-music算法的光学扫描全息单点定位方法 | |
CN108931298B (zh) | 一种紧凑型高通量高稳定性干涉成像光谱仪 | |
US10607360B2 (en) | Image processing apparatus and non-transitory computer-readable storage medium | |
EP3582183B1 (en) | Deflectometric techniques | |
CN107145053B (zh) | 基于tr-music算法的光学扫描全息轴向定位方法 | |
CN106886143B (zh) | 基于相位全息相关性的全息扫描空间距离提取方法 | |
CN104949940A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190719 Termination date: 20200419 |
|
CF01 | Termination of patent right due to non-payment of annual fee |