CN105806240A - 基于光学传递函数的同时测量多个绝对距离的方法 - Google Patents

基于光学传递函数的同时测量多个绝对距离的方法 Download PDF

Info

Publication number
CN105806240A
CN105806240A CN201610353004.7A CN201610353004A CN105806240A CN 105806240 A CN105806240 A CN 105806240A CN 201610353004 A CN201610353004 A CN 201610353004A CN 105806240 A CN105806240 A CN 105806240A
Authority
CN
China
Prior art keywords
diaphragm opening
lambda
diaphragm
mtf
delta
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
CN201610353004.7A
Other languages
English (en)
Other versions
CN105806240B (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201610353004.7A priority Critical patent/CN105806240B/zh
Publication of CN105806240A publication Critical patent/CN105806240A/zh
Application granted granted Critical
Publication of CN105806240B publication Critical patent/CN105806240B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/14Measuring arrangements characterised by the use of optical techniques for measuring distance or clearance between spaced objects or spaced apertures

Landscapes

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

Abstract

本发明涉及基于光学传递函数的同时测量多个绝对距离的方法,是一种高精度、毫米量级测量范围、可实现多个绝对距离同时测量的方法,属于光电检测技术领域。该方法在出瞳设置光阑,光阑上设置多个离散光阑孔,分别采集参考光路和多个测量光路的光波,各光波通过会聚透镜在焦面上发生干涉‑衍射,得到系统焦面光强分布即点扩散函数(Point Spread Function,PSF),继而得到系统的光学调制传递函数(Modulation Transfer Function,MTF)非归一化侧峰值MTFph,再对其归一化得到MTFnph,再由推导得到的MTFnph与绝对距离ΔL的函数关系计算出ΔL,从而实现了高效、高精度的瞬态绝对距离测量。本方法提高了测量效率,测量精度高,可达到纳米量级,测量范围大,可以达到测量所用光源的相干程长。

Description

基于光学传递函数的同时测量多个绝对距离的方法
技术领域
本发明涉及基于光学传递函数的同时测量多个绝对距离的方法,是一种高精度、毫米量级测量范围、可实现多个绝对距离同时测量的方法,属于光电检测技术领域。
技术背景
在天文学、精密测量学、军事技术中,常需要对多个几纳米到几微米的微小台阶高度(即参考镜与被测镜沿光轴方向的绝对距离)进行高精度测量。针对多个微小台阶高度的高精度测量,国内外学者已经提出了很多不同的测量方法,这些方法总体上可以分为两类。
第一类是电学检测方法,常用的有电容测微法和电感测微法。电学检测方法可以实现多个微小台阶的并行测量,但容易受到环境温度、湿度以及电磁场分布的影响,从而影响测量精度。
第二类是光学检测方法,所要测量的实际台阶高度与参考光路和测量光路之间的光程差(OpticalPathDifference,OPD)相对应,对于反射式光路而言OPD是所测台阶高度ΔL的两倍,通过测量OPD即可得到实际微小台阶高度。常用的二维色散条纹法可以实现微小台阶的高精度测量,但这种方法每次只能测量一块被测镜相对参考镜的光程差。
发明内容
本发明的目的是为了解决上述技术问题,提出了一种基于光学传递函数的同时测量多个绝对距离的方法。该方法在出瞳设置光阑,光阑上设置多个离散光阑孔,分别采集参考光路和多个测量光路的光波,各光波通过会聚透镜在焦面上发生干涉-衍射,得到系统焦面光强分布即点扩散函数(PointSpreadFunction,PSF),继而得到系统的光学调制传递函数(ModulationTransferFunction,MTF)非归一化侧峰值MTFph,再对其归一化得到MTFnph,再由MTFnph与ΔL的关系计算出ΔL,从而实现了高效、高精度的瞬态绝对距离测量。为了能够同时测量多个绝对距离,应使与各个被测镜相对应的MTFnph不发生重叠,故需研究非冗余光阑孔的设置,这是本发明的关键所在。
本发明使用的光路,包含平行光光源、分光板、参考镜、多个被测镜、光阑、聚焦透镜,其中光阑放置在光路的出瞳处。光源发出的平行光经过分光板到达参考镜和被测镜上,经参考镜和被测镜反射,携带光程差信息的平面光波分别通过光阑上的离散光阑孔,再通过聚焦透镜在焦平面上,发生干涉和衍射。
本发明的目的通过以下技术方案实现。
基于光学传递函数的同时测量多个绝对距离的方法,包含步骤如下:
步骤一、在系统的出瞳处设置光阑,光阑上设置的N个离散的圆光阑孔分别采集参考光路与N-1路测量光路的光波,携带光程差信息的平面光波通过光阑孔后发生衍射-干涉。定义参考光阑孔与被测光阑孔的中心连线为基线。为了实现光阑孔的非冗余排列,其余的光阑孔设置原则如下:以参考光阑孔作为基准,首先尽量将被测光阑孔设置在不同基线方向上;其次对于在同一基线上的被测光阑孔,应保证任意两个光阑孔中心距不与被测光阑孔和参考光阑孔的中心距相同;当任意两个被测光阑孔中心距与某基线上被测光阑孔和参考光阑孔的中心距相同时,两被测光阑孔中心连线应不与基线平行。
步骤二、光波经过聚焦透镜在焦平面上发生干涉-衍射现象,焦面上的光强分布即为PSF,对PSF进行傅里叶变换、取模,即可得到系统的MTF。MTF包含一个主峰、多个侧峰,标记出每一个测量光波与参考光波形成的MTF侧峰位置,记录MTF侧峰归一化峰值MTFnph的数值。
步骤三、推导MTFnph与ΔL的函数关系。依据傅里叶光学,两个光阑孔形成的PSF为:
PSF c ( x , y ) = 2 ∫ λ 0 - Δ λ 2 λ 0 + Δ λ 2 ( D 2 ) 2 f 2 λ 2 J 1 2 ( π D λ f x 2 + y 2 ) x 2 + y 2 [ 1 + cos k ( Δ L - B f x ) ] d λ - - - ( 1 )
其中,λ0为所用光源的中心波长,Δλ表示带宽,D为光阑孔的直径,J1()表示一阶贝塞尔函数,B为两个光阑孔的距离。
利用微分求和对(1)式进行化简:
PSF c ( x , y ) = = 2 Δ λ n Σ i = 1 n { D 2 2 f 2 λ i 2 J 1 2 ( π D λ i f x 2 + y 2 ) x 2 + y 2 [ 1 + c o s 2 π λ i ( Δ L - B f x ) ] } - - - ( 2 )
其中n为所用波长带宽采样数。
对(2)式进行傅里叶变换、归一化处理,可得侧峰归一化峰值MTFnph与ΔL的函数关系:
MTF n p h = 1 n n + [ Σ j = 1 n - 1 Σ i = j n - 1 cos ( 2 π λ i Δ L + 2 π λ i + 1 Δ L ) ] - - - ( 3 )
步骤四、将步骤二测得的各个侧峰归一化峰值带入公式(3),计算出N-1路绝对距离。
有益效果
本发明对比已有技术具有以下创新点:
通过在出瞳处设置带有N离散个光阑孔的光阑,实现了多个绝对距离的同时测量。
本发明对比已有技术具有以下显著优点:
1)可以同时测量N-1个被测镜与参考镜的绝对距离,提高测量效率。
2)测量精度高,可达到纳米量级;
3)测量范围大,可以达到测量所用光源的相干程长。当使用波长为632.8nm,带宽为1nm时,相干程长为400μm,测量范围0-±400μm。
附图说明
图1为本发明的检测分块镜共相位误差光路原理图(图1A)和主镜图(图1B);
其中:1、望远镜次镜;2、参考镜;3、被测镜1;4、被测镜2;5、准直透镜;6、光阑;7、聚焦透镜;8、探测器
图2为本发明的光阑孔设置图;
图3为本发明的系统光学调制传递函数;
图4为本发明的MTFnph与ΔL的关系曲线;
图5为本发明测量过程流程图。
具体实施方式
下面结合附图与实施例对本发明作进一步的说明。
实施例1
将本发明的方法用于检测分块式主镜望远镜共相位误差,望远镜主镜由18块子镜拼接而成如图1B所示,光路原理如图1A所示。
分块式主镜望远镜的初始共相位误差为亚毫米量级,在此共相位误差对应于被测镜与标准镜间的轴向距离。工作时,望远镜的光源通常是星光,可视为平行光,以其中一块分块子镜作为参考镜,其余子镜作为被测镜,平面波经参考镜和17个被测镜反射,反射光波经过光阑孔,反射光波中携带了参考镜和各被测镜间的绝对距离ΔL即共相位误差的信息,再经望远镜次镜1反射,由准直透镜4变为平行光波,它们之间的光程差OPD为ΔL的2倍。
步骤一、在系统的出瞳处设置光阑如图2所示,光阑上设置的18个离散的圆光阑孔分别采集参考光路与17路测量光路的光波,携带光程差信息的平面光波通过光阑孔后发生衍射-干涉。定义参考光阑孔与被测光阑孔的中心连线为基线。为了实现光阑孔的非冗余排列,其余的光阑孔设置原则如下:以参考光阑孔作为基准,首先尽量将被测光阑孔设置在不同基线方向上;其次对于在同一基线上的被测光阑孔,应保证任意两个光阑孔中心距不与被测光阑孔和参考光阑孔的中心距相同;当任意两个被测光阑孔中心距与某基线上被测光阑孔和参考光阑孔的中心距相同时,两被测光阑孔中心连线应不与基线平行。被测光阑孔的具体排列方式如下:
(1)、在竖直基线方向上设置三个测量光阑孔,参考光阑孔上方设置一个光阑孔,与参考光阑孔的中心距为a;下方设置两个光阑孔,孔中心距分别为b和c;a:b:c=2.27:4:5.73。
(2)、在与竖直基线夹角为60°的基线方向上设置三个测量光阑孔,参考光阑孔的左侧设置一个光阑孔,中心距离为d;参考光阑孔的右侧设置两个光阑孔,与参考光阑孔的中心距分别为e和f;d:e:f=1.11:2.34:4。
(3)、在与竖直基线夹角为-60°的基线方向上光阑孔设置与(2)的原则相同。
(4)、在水平基线方向上,设置两个测量光阑孔,左侧光阑孔与参考光阑孔的中心距为g;右侧光阑孔与参考光阑孔的中心距为h;g:h=2.9:3.7。
(5)、在与竖直基线夹角为19°的方向上,设置两个被测光阑孔,与参考光阑孔的中心距分别为i和j;i:j=2.7:5。
(6)、在在与竖直基线夹角为42°的方向上,设置一个被测光阑孔,距离参考光阑孔的中心距离为k;j:k=5:4.3。
(7)、在与竖直基线夹角为-19°的基线方向上光阑孔设置与(5)的原则相同。
(8)、在与竖直基线夹角为-42°的基线方向上光阑孔设置与(5)的原则相同。
步骤二、光波经过聚焦透镜在焦平面上发生干涉-衍射现象,其焦面上的光强分布即为PSF,对所得到的PSF进行傅里叶变换、取模,即可得到系统的MTF如图3所示。MTF包含一个主峰、多个侧峰,标定出每一个每一束测量光波与参考光波形成的MTF侧峰位置,记录MTF侧峰归一化峰值MTFnph的数值。
步骤三、推导MTFnph与ΔL的函数关系。依据傅里叶光学,两个光阑孔形成的PSF为:
PSF c ( x , y ) = 2 ∫ λ 0 - Δ λ 2 λ 0 + Δ λ 2 ( D 2 ) 2 f 2 λ 2 J 1 2 ( π D λ f x 2 + y 2 ) x 2 + y 2 [ 1 + cos k ( Δ L - B f x ) ] d λ - - - ( 1 )
其中,λ0为所用光源的中心波长,Δλ表示带宽,D为光阑孔的直径,J1()表示一阶贝塞尔函数,B为两个光阑孔的距离。
利用微分求和对(1)式进行化简:
PSF c ( x , y ) = = 2 Δ λ n Σ i = 1 n { D 2 2 f 2 λ i 2 J 1 2 ( π D λ i f x 2 + y 2 ) x 2 + y 2 [ 1 + c o s 2 π λ i ( Δ L - B f x ) ] } - - - ( 2 )
其中n为所用波长带宽采样数。
对(2)式进行傅里叶变换、归一化处理,可得侧峰归一化峰值MTFnph与ΔL的函数关系:
MTF n p h = 1 n n + [ Σ j = 1 n - 1 Σ i = j n - 1 cos ( 2 π λ i Δ L + 2 π λ i + 1 Δ L ) ] - - - ( 3 )
关系曲线如图4所示。
步骤四、将步骤二测得的各个侧峰归一化峰值带入公式(1),计算出17路共相位误差的大小。
本发明的主要特点:
与电学检测方法相比,本发明受环境影响较小,可以保证很高的精度。
与二维色散条纹法相比,本发明在保证大测量量程的前提下,可同时测量N-1块分块镜共相位误差。
本发明是一种非接触式的瞬态绝对距离并行测量方法,测量范围可达毫米量级,可实现N-1个绝对距离的同时测量。特别适用于大口径分块式主镜望远镜的共相位检测。

Claims (1)

1.基于光学传递函数的同时测量多个绝对距离的方法,其特征在于,在出瞳设置光阑,光阑上设置多个离散光阑孔,分别采集参考光路和多个测量光路的光波,各光波通过会聚透镜在焦面上发生干涉-衍射,得到系统焦面光强分布即点扩散函数PSF,PointSpreadFunction,继而得到系统的光学调制传递函数MTF,ModulationTransferFunction的非归一化侧峰值MTFph,再对其归一化得到MTFnph,再由MTFnph与ΔL的关系计算出ΔL;
具体步骤包括如下:
步骤一、在系统的出瞳处设置光阑,光阑上设置的N个离散的圆光阑孔分别采集参考光路与N-1路测量光路的光波,携带光程差信息的平面光波通过光阑孔后发生衍射-干涉;定义参考光阑孔与被测光阑孔的中心连线为基线;为了实现光阑孔的非冗余排列,其余的光阑孔设置原则如下:以参考光阑孔作为基准,首先尽量将被测光阑孔设置在不同基线方向上;其次对于在同一基线上的被测光阑孔,应保证任意两个光阑孔中心距不与被测光阑孔和参考光阑孔的中心距相同;当任意两个被测光阑孔中心距与某基线上被测光阑孔和参考光阑孔的中心距相同时,两被测光阑孔中心连线应不与基线平行;
步骤二、光波经过聚焦透镜在焦平面上发生干涉-衍射现象,焦面上的光强分布即为PSF,对PSF进行傅里叶变换、取模,即可得到系统的MTF;MTF包含一个主峰、多个侧峰,标记出每一个测量光波与参考光波形成的MTF侧峰位置,记录MTF侧峰归一化峰值MTFnph的数值;
步骤三、推导MTFnph与ΔL的函数关系;依据傅里叶光学,两个光阑孔形成的PSF为:
PSF c ( x , y ) = 2 ∫ λ 0 - Δ λ 2 λ 0 + Δ λ 2 ( D 2 ) 2 f 2 λ 2 J 1 2 ( π D λ f x 2 + y 2 ) x 2 + y 2 [ 1 + cos k ( Δ L - B f x ) ] d λ - - - ( 1 )
其中,λ0为所用光源的中心波长,Δλ表示带宽,D为光阑孔的直径,J1()表示一阶贝塞尔函数,B为两个光阑孔的距离;
利用微分求和对(1)式进行化简:
PSF c ( x , y ) = = 2 Δ λ n Σ i = 1 n { D 2 2 f 2 λ i 2 J 1 2 ( π D λ i f x 2 + y 2 ) x 2 + y 2 [ 1 + c o s 2 π λ i ( Δ L - B f x ) ] } - - - ( 2 )
其中n为所用波长带宽采样数;
对(2)式进行傅里叶变换、归一化处理,可得侧峰归一化峰值MTFnph与ΔL的函数关系:
MTF n p h = 1 n n + [ Σ j = 1 n - 1 Σ i = j n - 1 cos ( 2 π λ i Δ L + 2 π λ i + 1 Δ L ) ] - - - ( 3 )
步骤四、将步骤二测得的各个侧峰归一化峰值带入公式(3),计算出N-1路绝对距离。
CN201610353004.7A 2016-05-25 2016-05-25 基于光学传递函数的同时测量多个绝对距离的方法 Active CN105806240B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610353004.7A CN105806240B (zh) 2016-05-25 2016-05-25 基于光学传递函数的同时测量多个绝对距离的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610353004.7A CN105806240B (zh) 2016-05-25 2016-05-25 基于光学传递函数的同时测量多个绝对距离的方法

Publications (2)

Publication Number Publication Date
CN105806240A true CN105806240A (zh) 2016-07-27
CN105806240B CN105806240B (zh) 2018-03-30

Family

ID=56451883

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610353004.7A Active CN105806240B (zh) 2016-05-25 2016-05-25 基于光学传递函数的同时测量多个绝对距离的方法

Country Status (1)

Country Link
CN (1) CN105806240B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109668526A (zh) * 2019-01-29 2019-04-23 北京理工大学 基于光学传递函数的高精度测量倾斜角的方法
CN111220971A (zh) * 2019-10-22 2020-06-02 北京理工大学 不受倾斜角影响的高精度测量绝对距离的方法
CN111811474A (zh) * 2020-07-20 2020-10-23 北京理工大学 基于光学传递函数同时测量倾斜角和绝对距离的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0370770A1 (en) * 1988-11-22 1990-05-30 Hewlett-Packard Company Improved ranging system
CN101221042A (zh) * 2008-01-23 2008-07-16 北京理工大学 用于绝对距离测量的二维色散条纹分析方法
CN105547168A (zh) * 2016-01-27 2016-05-04 北京理工大学 用于绝对距离测量的光学调制传递函数分析方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0370770A1 (en) * 1988-11-22 1990-05-30 Hewlett-Packard Company Improved ranging system
CN101221042A (zh) * 2008-01-23 2008-07-16 北京理工大学 用于绝对距离测量的二维色散条纹分析方法
CN105547168A (zh) * 2016-01-27 2016-05-04 北京理工大学 用于绝对距离测量的光学调制传递函数分析方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
J.F.SIMAR 等: "Single wavelength coarse phasing in segmented telescopes", 《APPLIED OPTICS》 *
曲丹丹 等: "基于衍射理论的分块镜共相位误差的高精度测量", 《红外与激光工程》 *
王姗姗 等: "空间拼接主镜望远镜共相位检测方法", 《光学学报》 *
赵伟瑞 等: "Active cophasing and aligning testbed with segmented mirrors", 《OPTICS EXPRESS》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109668526A (zh) * 2019-01-29 2019-04-23 北京理工大学 基于光学传递函数的高精度测量倾斜角的方法
CN111220971A (zh) * 2019-10-22 2020-06-02 北京理工大学 不受倾斜角影响的高精度测量绝对距离的方法
CN111220971B (zh) * 2019-10-22 2021-10-19 北京理工大学 不受倾斜角影响的高精度测量绝对距离的方法
CN111811474A (zh) * 2020-07-20 2020-10-23 北京理工大学 基于光学传递函数同时测量倾斜角和绝对距离的方法
CN111811474B (zh) * 2020-07-20 2021-07-09 北京理工大学 基于光学传递函数同时测量倾斜角和绝对距离的方法

Also Published As

Publication number Publication date
CN105806240B (zh) 2018-03-30

Similar Documents

Publication Publication Date Title
CN101251484B (zh) 基于调制的微型傅里叶变换光谱仪
CN102944312B (zh) 一种测量部分相干涡旋光束拓扑荷数的方法
CN103197322B (zh) 一种飞秒激光频率梳合成波干涉测距方法及测距系统
CN105181298B (zh) 多次反射式激光共焦长焦距测量方法与装置
CN103884491B (zh) 一种扫描相机摆镜二维动态角测量校准方法与装置
CN102313642B (zh) 一种高精度长焦距透镜的焦距检测装置
CN102564611B (zh) 大功率激光波前测量仪及波前测量方法
CN103471524B (zh) 共焦抛物面顶点曲率半径测量方法
CN104568389A (zh) 双边错位差动共焦元件参数测量方法
CN101852676B (zh) 多焦全息差动共焦超长焦距测量方法与装置
CN102252824B (zh) 基于泰伯效应的复合差分式长焦距测量装置
CN102679895B (zh) 反射式共焦透镜中心厚度测量方法
CN104698468A (zh) 光纤光学相干测距装置和光纤光学测距方法
CN103615971B (zh) 用于检测圆柱体外表面的光学干涉仪
CN106225715A (zh) 一种用于非球面反射镜的五棱镜扫描检测方法
CN105806240A (zh) 基于光学传递函数的同时测量多个绝对距离的方法
CN102901463A (zh) 轴锥镜面形的测量装置和测量方法
CN102818541B (zh) 一种高分辨率的滚转角测量装置及测量方法
CN105547168B (zh) 用于绝对距离测量的光学调制传递函数分析方法
CN103615972A (zh) 用于检测空心圆柱内表面的光学干涉仪
CN106767679A (zh) 一种光电自准直经纬仪
CN103471525A (zh) 差动共焦抛物面顶点曲率半径测量方法
CN103697806A (zh) 用于检测环形导轨外圆弧面的光学干涉仪
CN205003080U (zh) 基于谱域干涉仪的折射率和厚度同步测量系统
CN201203578Y (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