CN101128717B - 用于干涉仪绝对距离测量的相位噪声补偿 - Google Patents

用于干涉仪绝对距离测量的相位噪声补偿 Download PDF

Info

Publication number
CN101128717B
CN101128717B CN200680005923XA CN200680005923A CN101128717B CN 101128717 B CN101128717 B CN 101128717B CN 200680005923X A CN200680005923X A CN 200680005923XA CN 200680005923 A CN200680005923 A CN 200680005923A CN 101128717 B CN101128717 B CN 101128717B
Authority
CN
China
Prior art keywords
interferogram
centerdot
distance measurement
absolute distance
expression
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
CN200680005923XA
Other languages
English (en)
Other versions
CN101128717A (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.)
Leica Geosystems AG
Original Assignee
Leica Geosystems AG
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 Leica Geosystems AG filed Critical Leica Geosystems AG
Publication of CN101128717A publication Critical patent/CN101128717A/zh
Application granted granted Critical
Publication of CN101128717B publication Critical patent/CN101128717B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02001Interferometers characterised by controlling or generating intrinsic radiation properties
    • G01B9/02002Interferometers characterised by controlling or generating intrinsic radiation properties using two or more frequencies
    • G01B9/02004Interferometers characterised by controlling or generating intrinsic radiation properties using two or more frequencies using frequency scans
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02055Reduction or prevention of errors; Testing; Calibration
    • G01B9/0207Error reduction by correction of the measurement signal based on independently determined error sources, e.g. using a reference interferometer
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B2290/00Aspects of interferometers not specifically covered by any group under G01B9/02
    • G01B2290/30Grating as beam-splitter
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B2290/00Aspects of interferometers not specifically covered by any group under G01B9/02
    • G01B2290/60Reference interferometer, i.e. additional interferometer not interacting with object

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Instruments For Measurement Of Length By Optical Means (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Optical Radar Systems And Details Thereof (AREA)

Abstract

本发明涉及用于干涉仪绝对距离测量的相位噪声补偿。在向至少一个待测目标发射调频电磁射线并随后在外差混频下接收由目标散射回的射线的测距方法中,使射线平行经过干涉测量基准长度。此外,在接收时产生由目标散射回的射线的数字化第一干涉图和经过基准长度的射线的数字化第二干涉图。从第二干涉图的相位曲线数据中合成出虚拟干涉图或者虚拟干涉图的相位曲线,在这里,通过比较第一干涉图的相位曲线数据与虚拟干涉图的相位曲线数据来确定距离。

Description

用于干涉仪绝对距离测量的相位噪声补偿
技术领域
本发明涉及根据权利要求1的前序部分的绝对距离测量方法、根据权利要求13的前序部分的绝对距离测量装置、以及计算机程序产品。
背景技术
在电子测距领域,公开了各种不同的原理和方法。一种做法是将经调频的电磁射线(如光线)射向待测目标,随后接收来自反向散射型物体(back-scattering object)的一个或多个反射信号,在理想情况下接收只来自待测目标的反射信号,在这里,待测目标不仅可能具有镜面反散射特性如反射镜特性,而且具有反向散射特性。在外差式干涉仪结构中,可调谐的激光源被用于绝对距离测量。在原理上最简单的实施例中,激光源的光频率是可线性调谐的。接收的信号与从所发出的光信号中导出的另一个信号叠加。所得到的外差混合积(mixed product)(即干涉图)的差频(beat frequency)是到目标物的距离的指标。
所发出的光信号的波长是测量的尺度。它一般是未知的,因而必须通过附加测量来确定。为此,例如使发出的光的一部分经过具有确定的基准长度(reference length)的基准干涉仪。根据所产生的差拍(beaproduct),可以基于已知的基准长度推断出所发出的光信号的随时间的波长变化。
用于实施这些方法的装置通常以激光器为光源,该光源能够被对外共振(如布拉格光栅)或内共振(如分布式反馈DFB或分布式布拉格反射器DBR)的调制所啁啾。在光学领域中,发射光学系统和接收光学系统被用于发射和接收,它们接在用于外差混频的检测器或正交检测器(quadrature detector)、A/D转换器和数字信号处理器的后面。
在实际应用中,主要存在以下多个基本难点:
·激光源的足以用于所需精度的线性调谐只能以高昂成本实现;
·激光源尤其是半导体激光二极管具有大的相位噪声,这意味着尺度长度(scale length)的变化不受控制,这一方面极大限制了测量范围,即将测量限制为测量小于激光相干长度的距离,另一方面,这造成大的测量波动,这种测量波动只随着长的测量时间减小。
·在测量过程中发生的待测目标运动以及测量距离波动在干涉仪混合积中造成额外的频率分量,其可能导致大的测量误差。
现有技术公开了各种做法,这些做法都带有与之相关的问题。
在US4,830,486中描述了用于绝对距离测量的通用型系统。在这里,频率响应曲线的线性化借助控制措施实现,相位噪声抑制通过将干涉仪臂长选定为与目标距离近似一致来实现。测量过程中的待测目标运动可以通过适当选择频率响应曲线来加以区分。该系统允许测量反射性很低的目标,但为了达到所需的高测量精度而需要长时间的测量。此外,要求基准干涉仪的臂长处于待测距离的数量级内,这导致更大的复杂性或者更高的技术成本以及高昂的生产成本。
DE 19522262描述了一种系统,该系统具有两个激光源,以产生合成波长。为此,尽管对干涉仪的结构稳定性以及目标距离的稳定性的要求不高,但该系统由于有两个激光源而更复杂。
在US2003/0107743A1中描述了用于消除信号源的相位噪声的方法。该方法用于测量具有大的光学长度的光学元件的组穿越时间(grouptransit time)。它不适于以所描述的形式测量绝对距离。
发明内容
本发明的目的是提供这样的解决方案,针对上述通用型的绝对距离测量仪,该方案允许更好地识别和掌握频率响应曲线和/或允许减小或避免误差及其影响。具体地,在基准干涉仪内,发出的激光的调谐特性和相位噪声应具有以下特征,即具体地可以通过计算来考虑非线性和噪声,并可以消除其对测量精度的影响。
本发明的另一个目的是也可以实现超过相干长度(coherence length)的干涉仪测量。
通过权利要求1或13以及从属权项的保护主题来完成这些任务或对解决方案的进一步改进。
根据本发明,在基准干涉仪中确定信号源的相位曲线。所确定出的相位曲线数据随后被用于实际消除激光源非理想特性的影响。
对于发出的激光的给定频率响应曲线v(t),基准干涉仪的外差信号的相位曲线在长度(臂长差)为Rr=c·τr/2的情况下如此得到:
Δ φ r ( t k ) = 2 π · ∫ t k - τ r t k v ( t ) dt - - - ( 1 )
在理想的实施方式中,该相位曲线信号在采样点tk被采样和数字化,其中,
tk=α·k·τr         (2)
修正系数α≈1建模了采样间隔和基准长度的滞后时间之间的偏离。就是说,采样间隔在理想情况下精确对应于基准长度的已知的滞后时间τr。根据本发明,可以有利地将采样频率选择为基准长度的滞后时间倒数的整数倍。
由以下公式得到该长度(臂长差)R=c·τm/2的测量距离的外差信号的相位曲线:
Δ φ M ( t k ) = 2 π · ∫ t k - τ m t k v ( t ) dt - - - ( 3 )
根据现有技术的方法,对于静立的目标物,距目标物的距离R通过解以下方程式得到:
R R r ≈ Δ φ M ( t ) Δ φ r ( t ) - - - ( 4 )
方程式(4)的计算可以在给定时刻T(测量间隔结束)进行,或作为优化问题被程式化。
从严格意义上讲,方程式(4)只适用于以下条件之一得到满足时:
·光学频率响应曲线是线性的,该条件不仅因为频率调谐而无法得到确切满足,而且由于相位噪声而无法得到确切满足。
·目标距离R与基准长度一致,R=Rr
根据现有技术,后者意味着,对于各希望的目标距离,必须可以得到具有相应长度的基准干涉仪。这显然实际上只能在特定情况下出现或者只能付出高昂成本来做到。
与方程式(4)不同,根据本发明,依据基准干涉仪(具有比所发出的激光的相干长度短的长度)的相位曲线数据的测量结果合成出等效于测量干涉仪(具有任意目标距离)的相位曲线。在最简单的实施例中,为此目的可应用:
Δ φ m · r ( t k ) = α · Σ j = k - m + 1 k Δ φ r ( t j ) - - - ( 5 )
其中,R≈m·α·Rr,就像对方程式(4)的有效性要求的那样。
在这里,术语“合成”原则上与在多波长干涉测量领域中常用的术语不同。在该领域中,距离信息根据各色光信号的频率差来确定,所述频率差产生“合成”的差频。在那里的术语“合成”表示只用(可测的)频差来工作。而在本发明的上下文中,“合成”表示纯理论上的、非可测的量,即例如计算确定虚拟干涉图(virtual interferogram)或者虚拟干涉图的相位曲线。
现在,方程式(4)被表述为:
R m · α · R r ≈ Δ φ M ( t ) Δ φ m · r ( t ) - - - ( 6 )
并且该方程式可实现对于激光频率调谐和相位噪声的任何理想曲线的较高精度。此外,将方程式(5)带入方程式(6)表明,可以取消修正系数α。在测量数据分析中确定未知参数m,该未知参数表示基准干涉仪和测量干涉仪之间的长度比。
该结构的测量精度能进一步提高,其做法是,方程式(5)由以下内插公式来代替,
Δ φ m · r ( t k ) = α · [ Σ j = k - m + 1 k Δ φ r ( t j ) + R - m · α · R r α · R r Δ φ r ( t k - m ) ] - - - ( 7 )
在这里,m被选择为(m-1)·α·Rr≤R≤m·α·Rr
在本发明中,可以在激光的相干长度之外进行测量。在这种情况下,相位曲线在测量路程中退变为噪声过程。实际测量路程的合成相位曲线的形成产生相关的噪声处理。目标距离通过相关最大值的数字搜索来确定。
为了评估相位曲线数据,优选地通过以下类型的优化问题来替换方程式(6)的解:
min R ~ = [ R ~ m · α · R r Δ φ m · r ( t ) - Δ φ M ( t ) ] 2 - - - ( 8 )
根据本发明,为了简化系统,也可以以逼近方式来实现优化。
方程式(6)和(8)假定可以确定基准干涉图的相位曲线以及测量干涉仪的相位曲线。通常,这适用于基准测量。不过,测量路程上的测量结果通常被噪声严重污染,以至无法从外差信号s(t)≈A·cos(ΔφM(t))的测量中确定相位。在这种情况下,为了确定目标距离,直接依据干涉图的测定的信号的振幅(而不是依据推导出的相位),来解优化问题,
min A ~ , R ~ [ A ~ · cos ( R ~ m · α · R r Δ φ m · r ( t ) ) - s ( t ) ] 2 - - - ( 9 )
根据本发明,也可以以逼近方式实现用于系统简化的优化。
在测量目标以速度V在测量方向上移动时,具有长度R的测量干涉仪的外差信号的相位曲线由下式给出:
Δ φ M ( t k ) = 2 π · [ ∫ t k - τ m t k v ( t ) dt + V · t λ ] - - - ( 10 )
当在对测量数据的分析中未考虑运动时,这导致不正确的距离确定(多普勒效应)。按照本发明,在以下类型的相位曲线数据的优化问题中通过引入运动模型V(t;θ)来考虑测量目标的运动,其中θ为待确定的速度参数:
min R ~ , θ ~ [ R ~ m · α · R r Δ φ m · r ( t ) - Δ φ M ( t ) + 2 π V ( t ; θ ~ ) · t λ ] 2 - - - ( 11 )
当信号受到噪声的严重污染时,为了确定目标距离,直接依据测得的干涉图信号的振幅来解优化问题,而不是依据推导出的相位来解优化问题:
min A ~ , R ~ , θ ~ [ A ~ · cos ( R ~ m · α · R r Δ φ m · r ( t ) + 2 π V ( t ; θ ~ ) λ ) - s ( t ) ] 2 - - - ( 12 )
按照本发明,也可以采用其它运动模型。根据本发明,为了简化系统,优化也可以以逼近方式完成。在已知干涉仪稳定的情况下,可以省去对速度参数的确定。
逼近可以以各种方法实现。尤其是,可以采用速度恒定θ=v的简化运动模型,其中V(t;θ)=v。也可以顺序确定整数关系
Figure G200680005923XD00062
和待测距离
Figure G200680005923XD00063
或者可以省去对发出的激光随时间的变化γ(t)所进行的确定。为了在确定干涉图相位时避免2π含糊性,代替依照方程式(8)和(11)对相位曲线数据的优化,也可以进行对相位差(difference phase)ΔφM(tk)-ΔφM(tk-1)和Δφm·r(tk)-Δφm·r(tk-1)的优化。在干涉图被噪声严重污染的情况下,可以通过适当滤波将按照方程式(12)的干涉图信号的振幅优化减化为按照方程式(11)的相位曲线数据的优化。
基准干涉仪的长度以及数字化的采样间隔确定了测距仪的相对测量精度。虽然可以获得ppm精度的电子振荡器,但基准长度的稳定化是一个大问题。
根据本发明,可以如此解决该问题,即基准长度通过结构措施来稳定,或者利用一个本身恒温的或者温度被稳定的元件(如气室或恒温基准元件)来校准。为此目的,与对基准干涉仪的相位曲线的测量同时地由上述元件在激光频率调谐期间来确定透射极值。该极值的频率间隔Δv通过元件设计来非常精确地确定。在透射极值的周期t1、t2期间基准路程相位曲线的总相位变化Δφ=Δφr(t2)-Δφr(t1)现在是基准长度的指标。
R r = ΔΦ · c 4 π · Δv - - - ( 13 )
该测量与目标物的距离的测量同时或者先后进行。对于同时测量,可以通过滤波来改善测量精度(假定基准长度波动缓慢)。由此一来,相位噪声对基准长度校准的影响可以得到明显抑制,无需其它措施。
现在,通过本发明的做法,对频率调谐线性以及所用激光源的线宽的要求比现有技术的低。此外,测量范围将不会受激光相干长度的限制,测量精度即使在高速测量下也能得以保持。因此,可使用容易实现的、低成本的激光源(相应地,整个系统也具有成本优势),同时可得到高测量精度、宽测量范围和高测量速率。
附图说明
以下,结合附图示意表示的实施例以示例的方式更详细地说明或描述本发明的绝对距离测量方法或者本发明的绝对距离测量装置。
图1是基础测量结构的示意图。
图2是在发出的光信号和由被目标反射的光信号之间形成外差混合积的示意图。
图3表示通过正交接收器(quadrature receiver)记录的干涉图的相位。
图4表示单频信号的频率响应曲线以及时间积分窗口。
图5表示虚拟干涉图的频率响应曲线。
图6示意表示基准元件的第一实施方式。
图7示意表示基准元件的第二实施方式。
图8表示乙炔的吸收光谱。
图9表示乙炔的吸收光谱片断。
图10示意表示测量干涉仪和基准干涉仪的完整的光纤设计,包括用于校准基准长度的基准元件。
具体实施方式
图1表示基础测量构造。信号源SQ通过光学连接器OV1-OV4被耦合输入至测量干涉仪MI、基准干涉仪RI和可选部件校准单元KE。待测目标ZI在测量中是测量干涉仪的组成部分。通过光程OS2、OS3和OS4,信号被传输给检测器单元D2、D3和D4。检测器信号通过电缆L2、L3和L4被送往信号处理器SP。信号处理器SP负责处理和分析所述信号并且还通过控制线SL控制信号源SQ。
图2示出了用于形成外差混合积的干涉仪的示意图。发出的光线的一部分由分光器1耦合输出,在反射器2上反射并经过分光器后,与从反射器3散射回的光信号叠加。对于干涉仪绝对距离测量来说,采用可调谐的激光源。在原理上最简单的实施方式中,取决于时间的发射满足方程式v(t)=v0+γt。由于所经过的路程而被延迟的两个电场(e-field)的叠加产生了外差频率f。
f=v(t2)-v(t1)=(v0+γt2)-(V0+γt1)
t1=2L1/c,t2=2L2/c
f = γ 2 R c
R=L2-L1(路程差)
该基本关系式将待测距离与外差频率f的载波信号联系起来。
图3示出了由正交接收器记录的干涉图的相位。对于γ=40GHz/ms和L=40m的情况来说,作为测量的外差频率,f=10.7MHz。检测器上出现的接收信号强度为:
I=(E1+E2)·(E1+E2)*
I=E1 2+E2 2+2E1E2·cos(2πft)
接收的强度是信号处理的基础。如果信号未受干扰,则通过傅立叶分析能直接得到载频,从而可以直接计算待确定的距离。但实际上信号受到干扰,因此直接由光学的正交接收器测量波动强度的相位φ是有利的。在准确测量时适用线性曲线:
Δφ Δt = f
因此,如图3所示,相位随时间以线性增加。
图4表示单频信号的频率响应曲线。此外,示出了按照式(1)和式(3)的时间相关积分窗口。干涉仪检测发出的信号e(t)和在经过一段路程后又接收的信号r(t)之间的相位改变。为了限定对相位变化的测量,采用单频信号或由多个单频分量组成的信号。这意味着,可以基于信号的振幅E(t)和瞬间频率v(t)实现特性化,它们随着测量过程比较缓慢地变化。
e(t)=E(t)·cosφ(t)
φ ( t ) 2 π = ∫ 0 t v ( t ′ ) dt ′ + φ 0 2 π
其中,e(t)表示发出的信号,E(t)表示信号的振幅,φ(t)表示发出的信号的相位(rad),φ0表示相移(rad),v(t)表示发出的信号的瞬间频率(Hz)。
图4示出了瞬间频率v(t)的随时间变化曲线的例子。理想情况下,接收的信号r(t)是发出的信号e(t)的延迟衰减的复制,通过以下项τ(t)和G(t)来考虑延迟和衰减:
r ( t ) = G ( t ) · e ( t - τ ( t ) ) ≈ G ( t ) · E ( t ) · cos φ ( t - τ ( t ) )
目标距离的确定通过测量从测距仪向目标、并从目标返回测距仪这样的每个循环的相位变化、并进而测量r(t)和e(t)之间相位变化来实现。
Δφ(t)=φ(t)-φ(t-τ(t))
在如上述假定的仅缓慢变化并且忽略可能有的速度分量的情况下,遵循
&Delta;&phi; ( t ) 2 &pi; = &Integral; t - &tau; ( t ) &tau; ( t ) < < v ( t ) / v &CenterDot; ( t ) t v ( t &prime; ) dt &prime;
&ap; v ( t ) &CenterDot; &tau; ( t )
不过,该假定
Figure G200680005923XD00095
仅适用于频率缓慢变化,像例如可能由频率响应曲线的非线性引起的那样。这种假定不足以应对高频噪声分量。
图5示出了虚拟干涉图的频率响应曲线。根据本发明,根据长度短的基准干涉仪测量的相位曲线数据,合成产生对应于测量距离的长度的等同测量数据。在该最简单的实施方式中,计算
&Delta; &phi; m &CenterDot; r ( t k ) = &alpha; &CenterDot; &Sigma; j = k - m + 1 k &Delta; &phi; r ( t j )
作为基准测量的单独的相位间隔之和,其中R≈m·α·Rr。图5中的频率响应曲线示出了这样的虚拟干涉图的相位组成的一个示例。
在图6中示意示出了用于确定基准长度的基准元件的第一实施方式。干涉仪完全整合在微晶玻璃体(Zerodur-Block)4中。在入口窗上装有第一光栅5,可以在出口侧给它补充另一个光栅8。入射光束被分成两个光路,它们对应于零阶和一阶的极值。零阶的第一光路7直接经过微晶玻璃体4传播,而一阶的第二光路6在微晶玻璃体4的界面上多次(全)反射。通过精确确定尺寸,可以在出口窗获得这两个彼此相对相移的光路的叠加。光路的光程差在输出后叠加时提供了具有从微晶玻璃体的材料测量(material measure)中知道的频率的外差信号。射线可以由检测器记录和处理。通过微晶玻璃材料,可以相当程度地摆脱温度的影响。由于例如可以以光刻技术施加的光栅,光传导只能在结构紧凑的恒温的微晶玻璃体4中实现。除了两维分解外,也可以实现传播方向和光路在三条轴线上的三维分解。
图7示意示出了具有光检测器的基准原件的第二实施方式,经过多次折射的光路射向光检测器。基准元件尤其适用于窄的波长范围。实现的基础还是微晶玻璃体4’,第一光栅9和任选的其他光栅10、10’、10”位于其表面上或其表面内。在该实施方式中,光束引导(beam guidance)的角度是灵敏的信息。在所示结构中,在微晶玻璃中,1微米波长折射角度为22.5度,这对应于约4.2微米的光栅周期(grating period)。所需分辨率(1ppm)确定了被照射的光栅部的宽度。按照在此以示例方式示出的四次折射后,图像被送入检测器装置。该检测器约以100兆赫进行读出。各测量点代表调谐激光时的波长的瞬时记录值。
图8和图9示出了作为基准元件的第三实施方式的气室。该实施方式利用自然的波长基准,其由介质中的吸收率限定并与HeNe激光中的稳定化相似。例如,为此目的可以采用乙炔气室,它可以作为元件买到。此外,使用范围为1540纳米、线宽0.1pm。如果调谐超过1纳米(120GHz),可以利用高精度的2-3过渡(2-3transition),因此能很精确地确定激光波长并校准基准干涉仪。由于过渡的原子属性,所以可以忽略与温度相关的影响。
具体来说,图8表示乙炔的吸收光谱,图9表示乙炔的吸收光谱片断。
图10示出了利用在完整的光纤结构中的基准元件的第三实施方式的气室。激光源12的光穿过光绝缘体13,用于避免外部反馈。在绝缘体13后,光在第一分光器14中非对称地被耦合(10dB耦合器)输入两个干涉仪,一个输出在终接件15结束。两个干涉仪又各有一个3×3耦合器16,它造成两个检测器之间的相移,这两个检测器分别构成两个对应的正交接收器17之一。在两个干涉仪之一中设有气室18。两个基准干涉仪中的法拉第镜19-19”此时在叠加中产生正交极化的反射波,调制深度100%。通过借助光纤的实现,获得了射线的良好叠加,而不需要方向校准或调整。单模特性允许通过发射光学系统21发出平的波阵面(wavefront),用于测量目标22。两个光二极管20’用于监视。

Claims (25)

1.一种绝对距离测量方法,该绝对距离测量方法包括以下步骤:
●将经过啁啾的电磁射线,发射
○到至少一个待测目标,并
○经过至少一个基准长度,
●接收从目标所散射回的射线和经过该基准长度的射线,
●将接收的射线转换为信号,
●根据所述信号确定距所述至少一个目标的至少一个距离,
其中,在接收时,记录从目标散射回的射线的数字化的第一干涉图和经过该基准长度的射线的数字化的第二干涉图,
其特征在于,根据所述第二干涉图的相位曲线数据合成出虚拟干涉图或者虚拟干涉图的相位曲线Δφm·r(tk),确定所述至少一个距离的步骤通过将所述第一干涉图的相位曲线数据与所述虚拟干涉图的相位曲线数据相比较来实现。
2.根据权利要求1所述的绝对距离测量方法,其特征在于,以将第一干涉图与第二干涉图相关联的方式、根据所述第二干涉图的相位曲线数据合成出虚拟干涉图,由此减小相位曲线波动对距离确定的干扰影响。
3.根据权利要求1所述的绝对距离测量方法,其特征在于,基准长度与所述数字化的第二干涉图的采样间隔成已知比例。
4.根据权利要求1所述的绝对距离测量方法,其特征在于,所述虚拟干涉图的相位曲线按下式形成,
&Delta; &phi; m &CenterDot; r ( t k ) = &alpha; &CenterDot; &Sigma; j = k - m + 1 k &Delta;&phi; r ( t j ) 或者
&Delta; &phi; m &CenterDot; r ( t k ) = &alpha; &CenterDot; [ &Sigma; j = k - m + 1 k &Delta;&phi; r ( t j ) + R - m &CenterDot; &alpha; &CenterDot; R r &alpha; &CenterDot; R r &Delta;&phi; r ( t k - m ) ]
其中,Rr表示基准长度,m∈N表示基准长度和按照R≈m·α·Rr的测量距离之间比例的整数部分,其中适用(m-1)·α·Rr≤R≤m·α·Rr
Figure F200680005923XC00021
表示第二干涉图的相位曲线,α≈1表示对采样间隔和基准长度的滞后时间之差的偏离进行建模的修正系数,其中τr表示基准长度的滞后时间,v(t)表示发出的射线的给定频率响应曲线。
5.根据权利要求1所述的绝对距离测量方法,其特征在于,所述虚拟干涉图或者所述虚拟干涉图的相位曲线Δφm·r(tk)以多项式方式内插。
6.根据权利要求1所述的绝对距离测量方法,其特征在于,通过对所述第一干涉图的相位曲线数据和所述虚拟干涉图的相位曲线数据进行优化来实现距离确定,采用以下关系式,
min R ~ , m ~ [ R ~ m ~ &CenterDot; &alpha; &CenterDot; R r &Delta;&phi; m &CenterDot; r ( t ) - &Delta;&phi; M ( t ) ] 2
其中,
Figure F200680005923XC00023
表示所述第一干涉图的相位曲线,
Figure F200680005923XC00024
表示所述虚拟干涉图的相位曲线,
Figure F200680005923XC00025
为所述第二干涉图的相位曲线,τr表示基准长度的滞后时间,τm表示测量距离的滞后时间,v(t)表示发出的射线的给定频率响应曲线,Rr表示所述基准长度,
Figure F200680005923XC00026
表示待确定的距离,m∈N表示基准长度和按照R≈m·α·Rr的测量距离之间比例的整数部分,其中适用(m-1)·α·Rr≤R≤m·α·Rr,用对采样间隔和基准长度滞后时间之间的偏离进行建模的修正系数α≈1来修正。
7.根据权利要求1所述的绝对距离测量方法,其特征在于,距离确定步骤通过对第一干涉图的测量信号的振幅进行优化来实现,采用下式
min A ~ , R ~ , m ~ [ A ~ &CenterDot; &gamma; ( t ) &CenterDot; cos ( R ~ m ~ &CenterDot; &alpha; &CenterDot; R r &Delta;&phi; m &CenterDot; r ( t ) ) - s ( t ) ] 2
其中,s(t)表示所述第一干涉图的信号的振幅,
Figure F200680005923XC00028
表示所述第一干涉图的估算信号强度,γ(t)表示发出的作为时间的函数的激光强度曲线的任选测量,在已知输出功率恒定的情况下,可以设定γ(t)=1,表示虚拟干涉图的相位曲线,其中是所述第二干涉图的相位曲线,τr表示基准长度的滞后时间,v(t)表示发出的射线的给定频率响应曲线,Rr表示所述基准长度,表示待确定的距离,m∈N表示基准长度和按照R≈m·α·Rr的测量距离之间比例的整数部分,其中适用(m-1)·α·Rr≤R≤m·α·Rr,用对采样间隔和基准长度滞后时间之间的偏离进行建模的修正系数α≈1来修正。
8.根据权利要求1所述的绝对距离测量方法,其特征在于,在测量目标运动时,对所述距离的确定通过对所述第一干涉图的相位曲线数据和所述虚拟干涉图的相位曲线数据的优化来实现,采用下式
min R ~ , &theta; ~ , m ~ [ R ~ m ~ &CenterDot; &alpha; &CenterDot; R r &Delta;&phi; m &CenterDot; r ( t ) - &Delta;&phi; M ( t ) + 2 &pi; V ( t ; &theta; ~ ) &CenterDot; t &lambda; ] 2 或者
min A ~ , R ~ , &theta; ~ , m ~ [ A ~ &CenterDot; cos ( R ~ m ~ &CenterDot; &alpha; &CenterDot; R r &Delta;&phi; m &CenterDot; r ( t ) + 2 &pi; V ( t ; &theta; ~ ) &CenterDot; t &lambda; ) - s ( t ) ] 2 ,
表示所述第一干涉图的相位曲线,
Figure F200680005923XC00037
表示所述虚拟干涉图的相位曲线,其中
Figure F200680005923XC00038
是所述第二干涉图的相位曲线,τr表示基准长度的滞后时间,τm表示测量距离的滞后时间,V表示在测量方向上测量目标运动的速度,s(t)表示所述第一干涉图的信号的振幅,表示所述第一干涉图的估算信号强度,v(t)表示发出的射线的给定频率响应曲线,Rr表示基准长度,表示待确定的距离,表示待测目标的运动模型,其中
Figure F200680005923XC000312
是速度参数,
Figure F200680005923XC000313
表示基准长度和按照R≈m·α·Rr的测量距离之间比例的整数部分,其中适用(m-1)·α·Rr≤R≤m·α·Rr,用对采样间隔和基准长度的滞后时间之间的偏离进行建模的修正系数α≈1来修正。
9.根据权利要求6所述的绝对距离测量方法,其特征在于,所述优化以逼近方式完成,通过
●顺序地确定整数关系
Figure F200680005923XC00041
和待确定的距离
Figure F200680005923XC00042
●针对相位差ΔφM(tk)-ΔφM(tk-1)和Δφm·r(tk)-Δφm·r(tk-1)进行优化,而不是针对相位曲线数据进行优化。
10.根据权利要求1所述的绝对距离测量方法,其特征在于,对基准长度的校准借助于恒温的或者温度被稳定的基准元件来实现。
11.根据权利要求7所述的绝对距离测量方法,其特征在于,所述优化以逼近方式完成,通过
●顺序地确定整数关系
Figure F200680005923XC00043
和待确定的距离
Figure F200680005923XC00044
●省去对作为时间函数的发出的激光的变化γ(t)进行确定。
12.根据权利要求8所述的绝对距离测量方法,其特征在于,所述优化以逼近方式完成,通过
●速度恒定θ=v的简化运动模型,其中V(t;θ)=v,
●顺序地确定整数关系
Figure F200680005923XC00045
和待确定的距离
Figure F200680005923XC00046
●针对相位差ΔφM(tk)-ΔφM(tk-1)和Δφm·r(tk)-Δφm·r(tk-1)进行优化,而不是针对相位曲线数据进行优化。
13.根据权利要求1所述的绝对距离测量方法,其中,将经过啁啾的激光,发射
○到至少一个待测目标,并
○经过至少一个基准长度。
14.根据权利要求1所述的绝对距离测量方法,其中,利用干涉混合将接收的射线转换为信号。
15.根据权利要求3所述的绝对距离测量方法,其中,滞后时间τr的采样间隔按照Rr=c·τr/2与基准长度Rr相对应,其中c表示光速。
16.一种用于实施根据权利要求1至15之一所述方法的绝对距离测量装置,它具有:
可调制的发射源(SQ),用于产生射线、并向待测目标发射射线,
基准长度,
用于调制所述发射源的信号发生器,
检测器(D2,D3,D4),用于接收反向散射的射线的干涉图,并将其转换成信号,
用于处理信号的信号处理器(SP),
包括混频器,用于执行外差混频方法,
所述检测器(D2,D3,D4)和所述信号处理器(SP)被设置和设计为能记录从目标散射回的射线的数字化第一干涉图和经过该基准长度的射线的第二干涉图,
其特征在于,所述检测器(D2,D3,D4)和所述信号处理器(SP)被设置和设计为能根据所述第二干涉图的相位曲线数据合成出虚拟干涉图或者虚拟干涉图的相位曲线Δφm·r(tk)。
17.根据权利要求16所述的绝对距离测量装置,其特征在于,设有用于校准基准长度的恒温或者温度被稳定的基准元件。
18.根据权利要求17所述的绝对距离测量装置,其特征在于,所述基准元件是气室(18)。
19.根据权利要求17所述的绝对距离测量装置,其特征在于,所述基准元件是恒温的标准器。
20.根据权利要求17所述的绝对距离测量装置,其特征在于,所述基准元件是Zerodur块(4),其在入射面具有用于折射分离入射光束的第一光栅(5)、并在出射面具有用于汇集分离的光束的路径的第二光栅(8),所述分离的光束的光束路径差导致具有限定的周期的外差信号。
21.根据权利要求17所述的绝对距离测量装置,其特征在于,所述基准元件是Zerodur块,其具有用于根据波长将射线折射引导至检测器装置的至少一个光栅(9),该检测器装置使得可以确定信号源的波长。
22.根据权利要求16所述的绝对距离测量装置,其中,发射源(SQ)是可调谐的激光源。
23.根据权利要求16所述的绝对距离测量装置,其中,用于处理信号的信号处理器(SP)是数字信号处理器FPGA或ASIC。
24.根据权利要求16所述的绝对距离测量装置,其中,所述混频器是正交接收器。
25.根据权利要求18所述的绝对距离测量装置,其中所述气室(18)是填充有乙炔。
CN200680005923XA 2005-02-23 2006-02-09 用于干涉仪绝对距离测量的相位噪声补偿 Active CN101128717B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP05101373A EP1696201A1 (de) 2005-02-23 2005-02-23 Phasenrauschkompensation für interferometrische Absolutdistanzmesser
EP05101373.8 2005-02-23
PCT/EP2006/050797 WO2006089845A1 (de) 2005-02-23 2006-02-09 Phasenrauschkompensation für interferometrische absolutdistanzmesser

Publications (2)

Publication Number Publication Date
CN101128717A CN101128717A (zh) 2008-02-20
CN101128717B true CN101128717B (zh) 2010-04-14

Family

ID=34938800

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200680005923XA Active CN101128717B (zh) 2005-02-23 2006-02-09 用于干涉仪绝对距离测量的相位噪声补偿

Country Status (7)

Country Link
US (1) US7619719B2 (zh)
EP (2) EP1696201A1 (zh)
JP (1) JP5108534B2 (zh)
CN (1) CN101128717B (zh)
AU (1) AU2006217996B2 (zh)
CA (1) CA2599128C (zh)
WO (1) WO2006089845A1 (zh)

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2128560B1 (de) * 2008-05-28 2015-07-01 Leica Geosystems AG Interferometrisches Distanzmessverfahren mit spektral trennbarem Doppelchirp und ebensolche Vorrichtung
EP2128561B1 (de) 2008-05-28 2014-10-29 Leica Geosystems AG Interferometrisches Distanzmessverfahren mit verzögertem Chirp-Signal und ebensolche Vorrichtung
CA2783337C (en) * 2009-12-14 2015-01-27 Leica Geosystems Ag Method for speckle mitigation in an interferometric distance meter and corresponding distance meter
DE102010041634A1 (de) * 2010-09-29 2012-03-29 Carl Zeiss Ag Vorrichtung und Verfahren zur Messung des Abstandes eines Objekts von einem Bezugspunkt
EP2541193A1 (de) * 2011-06-27 2013-01-02 Hexagon Technology Center GmbH Interferometrisches Entfernungsmessverfahren zum Vermessen von Oberflächen und ebensolche Messanordnung
JP5868132B2 (ja) * 2011-11-14 2016-02-24 キヤノン株式会社 撮像装置および画像処理方法
DE102012217282B4 (de) * 2012-09-25 2023-03-02 Trimble Jena Gmbh Verfahren und Vorrichtung zur Zuordnung von Messpunkten zu einem Satz von Festpunkten
EP2816315B1 (en) * 2013-06-18 2015-09-23 Hexagon Technology Center GmbH Interferometric determination of distance change with laser diode, high bandwidth detection and fast signal processing
CN103487403B (zh) * 2013-10-14 2015-09-02 北京信息科技大学 带有参考腔补偿的双波长组合光纤激光器气体检测系统
EP3234500A4 (en) 2014-12-19 2018-07-04 University of Utah Research Foundation Interferometry system and associated methods
DE102015110362B4 (de) * 2015-06-26 2018-01-18 n-hands GmbH & Co. KG Verfahren und Vorrichtung zur interferometrischen Absolutmessung einer Entfernung
CN108027312B (zh) * 2015-07-30 2020-08-04 通快光电器件有限公司 用于颗粒尺寸检测的激光传感器
EP3475649A4 (en) 2016-06-23 2020-04-22 University of Utah Research Foundation INTERFEROMETRY SYSTEM AND RELATED METHODS
US11162781B2 (en) 2016-06-23 2021-11-02 University Of Utah Research Foundation Interferometry systems and methods
KR102477195B1 (ko) * 2016-11-30 2022-12-12 블랙모어 센서스 앤드 애널리틱스, 엘엘씨 광 처프 거리 검출의 도플러 검출 및 도플러 보정을 위한 방법 및 장치
CN107741591A (zh) * 2017-09-21 2018-02-27 深圳市盛路物联通讯技术有限公司 一种干涉测距方法及系统
CN111048978B (zh) * 2018-10-12 2021-04-20 华为技术有限公司 一种多波长激光器
CA3091463A1 (en) 2019-08-29 2021-02-28 Institut National De La Recherche Scientifique System and method for stabilization of multi-path optical interferometers
CN111445400B (zh) * 2020-03-13 2023-03-17 五邑大学 存在投影盲区去除绝对相位噪声方法、装置及存储介质
CN116819431B (zh) * 2023-08-31 2023-12-08 杭州岸达科技有限公司 基于反常相模激励的相位干涉仪测向方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0260894A1 (en) * 1986-09-12 1988-03-23 Cogent Limited Optical fibre measuring system
US5371587A (en) * 1992-05-06 1994-12-06 The Boeing Company Chirped synthetic wavelength laser radar
CN1099128A (zh) * 1994-03-04 1995-02-22 清华大学 用双波长激光进行外差干涉测量绝对距离系统
EP0646767A2 (en) * 1993-10-05 1995-04-05 Renishaw plc Interferometric distance measuring apparatus
CN1123573A (zh) * 1993-05-15 1996-05-29 莱卡公开股份有限公司 测距设备
CN1172532A (zh) * 1995-11-15 1998-02-04 莱卡公开股份有限公司 用于绝对距离的电光测量装置
CN1346965A (zh) * 2001-10-15 2002-05-01 天津大学 一种光纤干涉测量距离的方法及测量设备

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4830486A (en) 1984-03-16 1989-05-16 Goodwin Frank E Frequency modulated lasar radar
US4697926A (en) * 1985-02-08 1987-10-06 The Board Of Trustees Of The Leland Stanford Junior University Coherent distributed sensor and method using short coherence length sources
JP2764630B2 (ja) * 1990-03-20 1998-06-11 横河電機株式会社 アブソリュート測長器
JP2725434B2 (ja) * 1990-03-30 1998-03-11 横河電機株式会社 Fmヘテロダイン法を用いたアブソリュート測長方法およびアブソリュート測長器
DE19522262C2 (de) * 1995-06-20 1997-05-22 Zeiss Carl Jena Gmbh Heterodyn-Interferometer-Anordnung
GB9703928D0 (en) * 1997-02-25 1997-04-16 Aberlink Tech Ltd Apparatus for precise length measurement calibration
US6360609B1 (en) * 2000-02-23 2002-03-26 Massachusetts Institute Of Technology Method and system for interpreting and utilizing multimode dispersive acoustic guided waves
US6606182B2 (en) * 2001-11-16 2003-08-12 Oplink Communications, Inc. Tuning and temperature compensation of the air-gap etalon for dense wavelength-division multiplexing application
US6900895B2 (en) * 2001-12-06 2005-05-31 Agilent Technologies, Inc. Phase noise compensation in an interferometric system
US7116430B2 (en) * 2002-03-29 2006-10-03 Georgia Technology Research Corporation Highly-sensitive displacement-measuring optical device
US7009691B2 (en) * 2002-05-29 2006-03-07 Agilent Technologies, Inc. System and method for removing the relative phase uncertainty in device characterizations performed with a polarimeter
EP1707982A1 (en) * 2005-03-31 2006-10-04 AGELLIS Group AB Method for analysing a substance in a container

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0260894A1 (en) * 1986-09-12 1988-03-23 Cogent Limited Optical fibre measuring system
US5371587A (en) * 1992-05-06 1994-12-06 The Boeing Company Chirped synthetic wavelength laser radar
CN1123573A (zh) * 1993-05-15 1996-05-29 莱卡公开股份有限公司 测距设备
EP0646767A2 (en) * 1993-10-05 1995-04-05 Renishaw plc Interferometric distance measuring apparatus
CN1099128A (zh) * 1994-03-04 1995-02-22 清华大学 用双波长激光进行外差干涉测量绝对距离系统
CN1172532A (zh) * 1995-11-15 1998-02-04 莱卡公开股份有限公司 用于绝对距离的电光测量装置
CN1346965A (zh) * 2001-10-15 2002-05-01 天津大学 一种光纤干涉测量距离的方法及测量设备

Also Published As

Publication number Publication date
CN101128717A (zh) 2008-02-20
EP1851504B1 (de) 2016-11-09
AU2006217996B2 (en) 2010-08-05
US20080204698A1 (en) 2008-08-28
AU2006217996A1 (en) 2006-08-31
CA2599128A1 (en) 2006-08-31
CA2599128C (en) 2013-07-09
US7619719B2 (en) 2009-11-17
JP2008532001A (ja) 2008-08-14
EP1696201A1 (de) 2006-08-30
EP1851504A1 (de) 2007-11-07
WO2006089845A1 (de) 2006-08-31
JP5108534B2 (ja) 2012-12-26

Similar Documents

Publication Publication Date Title
CN101128717B (zh) 用于干涉仪绝对距离测量的相位噪声补偿
US5541730A (en) Interferometric measuring apparatus for making absolute measurements of distance or refractive index
EP1853952B1 (en) Compact fiber optic geometry for a counter-chirp fmcw coherent laser radar
EP0506297B1 (en) Three wavelength optical measurement apparatus
EP2606311B1 (en) Apparatus and method for measuring distance
CN101868688B (zh) 干涉仪设备及其运行方法
KR100328007B1 (ko) 다중패스 간섭법을 사용하여 공기의 굴절률을 측정하는 슈퍼헤테 로다인 방법 및 장치
US8363226B2 (en) Optical interference measuring apparatus
JP2015111160A (ja) 対チャープfmcwコヒーレントレーザレーダー用の小型の光ファイバ配置
KR20000068090A (ko) 슈퍼헤테로다인 간섭계 및 전자 주파수 증배를 사용하여 공기의굴절률을 보상하기 위한 방법
JPH0650842A (ja) 光反射率計の較正方法及び光コヒーレンス領域反射率計
US20220390496A1 (en) Photonic rydberg atom radio frequency receiver and measuring a radio frequency electric field
WO2023014740A1 (en) Photonic rydberg atom radio frequency receiver and measuring a radio frequency electric field
JPH03180704A (ja) レーザ干渉計
US20060126073A1 (en) Displacement measuring interferometer system and method using tunable lasers
US7280216B2 (en) Method and apparatus for determining the wavelength of an input light beam
CN100451581C (zh) 利用外差干涉法对激光波长进行测量的方法及装置
JPH09166414A (ja) 光計測装置
Matsumoto Recent interferometric measurements using stabilized lasers
US6847455B1 (en) Heterodyne interferometer with angstrom-level periodic nonlinearity
EP0920599B1 (en) Measuring the effects of the refractive index of a gas using different multiple path interferometry ( superhetrodyne )
JPH05203410A (ja) 光周波数領域反射点測定方法及びその装置
US20240151643A1 (en) Laser Interferometer
JPH01126582A (ja) レーザ測距装置
EP1058810B1 (en) Apparatus and methods for measuring intrinsic optical properties of a gas

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