CN110793440A - Optical deflection transient measurement method - Google Patents

Optical deflection transient measurement method Download PDF

Info

Publication number
CN110793440A
CN110793440A CN201911062078.5A CN201911062078A CN110793440A CN 110793440 A CN110793440 A CN 110793440A CN 201911062078 A CN201911062078 A CN 201911062078A CN 110793440 A CN110793440 A CN 110793440A
Authority
CN
China
Prior art keywords
phase
directions
pattern
shifting
patterns
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
CN201911062078.5A
Other languages
Chinese (zh)
Other versions
CN110793440B (en
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.)
China Jiliang University
Original Assignee
China Jiliang University
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 China Jiliang University filed Critical China Jiliang University
Priority to CN201911062078.5A priority Critical patent/CN110793440B/en
Publication of CN110793440A publication Critical patent/CN110793440A/en
Application granted granted Critical
Publication of CN110793440B publication Critical patent/CN110793440B/en
Expired - Fee Related 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
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/41Refractivity; Phase-affecting properties, e.g. optical path length

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Testing Of Optical Devices Or Fibers (AREA)

Abstract

本发明提供一种基于光学偏折的瞬态测量方法,涉及测量技术领域。本方法利用正交方向的频率载波条纹对x和y方向的四步移相条纹图(一共八幅移相图案)分别进行单独调制,将不同方向的四幅移相条纹的组合图案耦合到不同的颜色通道,最终得到一幅彩色图,即投影屏投影的图案。利用相机采集经过待测元件偏折后的变形图案,采用颜色分离,分离出x、y方向的复合图案,再通过解调操作得到八幅移相图案,即x、y两个方向各四幅移相条纹图,从中经过相位解调相关算法得到相位信息,经过计算得出被测波面的斜率,最后通过积分算法得到重构的波面。本方法解决了传统光学偏折检测中需要在待测元件上连续投影移相条纹而无法进行瞬态测量的技术问题。

Figure 201911062078

The invention provides a transient measurement method based on optical deflection, and relates to the technical field of measurement. This method uses the frequency carrier fringes in orthogonal directions to separately modulate the four-step phase-shifting fringes in the x and y directions (a total of eight phase-shifting patterns), and couples the combined patterns of the four phase-shifting fringes in different directions to different Color channel, and finally get a color image, that is, the pattern projected by the projection screen. A camera is used to collect the deformed pattern after the deflection of the element to be tested, and color separation is used to separate the composite pattern in the x and y directions, and then eight phase-shift patterns are obtained through demodulation operation, that is, four amplitude shift patterns in each of the x and y directions. Phase fringe diagram, from which phase information is obtained through the phase demodulation correlation algorithm, the slope of the measured wave surface is obtained through calculation, and finally the reconstructed wave surface is obtained through the integration algorithm. The method solves the technical problem that in the traditional optical deflection detection, the phase-shifting fringes are continuously projected on the component to be tested, and the transient measurement cannot be performed.

Figure 201911062078

Description

一种光学偏折瞬态测量方法An Optical Deflection Transient Measurement Method

技术领域technical field

本发明涉及测量技术领域,尤其是涉及一种基于频率载波法、相移法和彩色编码法的高精度大动态瞬时的测量方法。The invention relates to the technical field of measurement, in particular to a high-precision large dynamic instantaneous measurement method based on a frequency carrier method, a phase shift method and a color coding method.

背景技术Background technique

随着光学设计和制造的发展,光学元件在测试中的应用越来越广泛。人们利用波前测试方法对元件进行评价。传统的波前检测方法主要是采用干涉法,如Twyman-Green干涉仪、Fizeau干涉仪、点衍射干涉仪等。然而传统干涉仪的动态范围限制了它在实际情况下对大范围传输波前测试的应用。With the development of optical design and manufacturing, the application of optical components in testing is becoming more and more extensive. Components are evaluated using wavefront testing. Traditional wavefront detection methods mainly use interferometry, such as Twyman-Green interferometer, Fizeau interferometer, point diffraction interferometer, etc. However, the dynamic range of conventional interferometers limits its application to large-scale transmission wavefront measurements in practical situations.

近年来,偏折法已被广泛应用,并成功地应用于天文望远镜反射镜的测试,球面、非球面、镜面自由镜、精密X射线镜。与干涉法相比,偏折法对系统配置要求较低,只需要投影屏、CCD相机和待测元件等元件。基于光线追迹方法,确定投影屏通过待测元件后到相机的光线,并通过光线像差得到检测波前。利用二值法、灰度法、相移法、混合方法等方法来创建映射,实现光学偏折测量。多步移相法是常用的检测方法,投影屏连续显示投影相移正弦条纹。但是,该传统方法一次只能显示和采集一个移相图案。在传统的光学偏折测量中,需要在待测元件上连续投影至少三幅正弦移相条纹,因此传统的光学偏折测量无法做到瞬态测量,这也使得该方法在实际的应用中存在局限性,也限制了测量速度。In recent years, the deflection method has been widely used and successfully applied to the test of astronomical telescope reflectors, spherical, aspherical, specular free mirrors, and precision X-ray mirrors. Compared with the interferometric method, the deflection method has lower requirements on the system configuration, and only needs components such as a projection screen, a CCD camera, and the component to be tested. Based on the ray tracing method, the light from the projection screen to the camera after passing through the component to be tested is determined, and the detected wavefront is obtained through the light aberration. Use binary method, gray scale method, phase shift method, hybrid method, etc. to create maps and realize optical deflection measurement. The multi-step phase-shifting method is a commonly used detection method, and the projection screen continuously displays the projected phase-shifted sinusoidal fringes. However, this traditional method can only display and acquire one phase shift pattern at a time. In traditional optical deflection measurement, it is necessary to continuously project at least three sinusoidal phase-shift fringes on the component to be measured, so traditional optical deflection measurement cannot achieve transient measurement, which also makes this method suitable for practical applications. Limitations also limit the measurement speed.

发明内容SUMMARY OF THE INVENTION

针对现有技术的不足,本发明的目的在于解决传统方法中需要在被测对象上连续投影多幅相移条纹而无法进行瞬态测量的技术问题。本发明提供一种基于频率载波法、相移法和彩色编码法的高精度大动态瞬时测量方法。Aiming at the deficiencies of the prior art, the purpose of the present invention is to solve the technical problem that transient measurement cannot be performed because multiple phase shift fringes are continuously projected on the measured object in the traditional method. The invention provides a high-precision large dynamic instantaneous measurement method based on a frequency carrier method, a phase shift method and a color coding method.

本发明通过以下技术方案实现:一种光学偏折瞬态测量方法,包括步骤:The present invention is achieved through the following technical solutions: a method for measuring the transient state of optical deflection, comprising the steps of:

建立基于光学偏折法的光学检测系统,所述光学检测系统包括投影屏、待测元件、相机和计算机;Establish an optical detection system based on the optical deflection method, the optical detection system includes a projection screen, a component to be tested, a camera and a computer;

设计投影屏投影编码图案:利用正交方向的频率载波条纹对四步移相得到的八幅移相图案进行单独调制,将x、y两个不同方向的四幅移相条纹的组合图案耦合到两个不同的颜色通道,得到一幅彩色图;Design the projection coding pattern of the projection screen: use the frequency carrier fringes in orthogonal directions to separately modulate the eight phase-shifting patterns obtained by four-step phase-shifting, and couple the combined pattern of the four phase-shifting fringes in two different directions of x and y to the two. different color channels to get a color image;

图像解调:以所述彩色图投影至所述待测元件,并对经过所述待测元件偏折后相机采集到的变形彩色图进行颜色通道分离,得到x、y方向的复合图案,对所述复合图案进行解调得到x、y两个方向各四幅移相条纹图,一共八幅移相图案;Image demodulation: Project the color map to the component to be tested, and separate the color channels of the deformed color map collected by the camera after the component to be tested is deflected to obtain a composite pattern in the x and y directions. The composite pattern is demodulated to obtain four phase-shift fringe patterns in each of the x and y directions, with a total of eight phase-shift patterns;

通过四步移相算法和相位解包裹算法对上述的八幅移相图案进行处理,得到相位信息,从而得到实际光斑分布,将实际光斑分布与所述光学检测系统中通过光线追迹得到的理想光斑分布进行对比,得出光线像差,所述光线像差为所述待测元件加工误差作用的结果;The above-mentioned eight phase-shifting patterns are processed through the four-step phase-shifting algorithm and the phase unwrapping algorithm to obtain phase information, thereby obtaining the actual spot distribution, and comparing the actual spot distribution with the ideal obtained by ray tracing in the optical detection system The light spot distributions are compared, and the light aberration is obtained, and the light aberration is the result of the processing error of the element to be tested;

根据波像差与所述光线像差的关系得到波前的斜率分布,对波前的斜率进行积分完成波前重构。The slope distribution of the wavefront is obtained according to the relationship between the wave aberration and the light aberration, and the wavefront reconstruction is completed by integrating the slope of the wavefront.

优选的,所述设计投影屏投影编码图案包括步骤:Preferably, the designing the projection screen to project the coding pattern comprises the steps of:

在x,y两个方向各设计四幅强度呈正弦变化的移相图案Ipxn,Ipyn,利用正交方向的频率载波条纹对八幅移相图案进行单独调制,将x方向的四幅移相条纹的组合图案耦合到红色R通道,将y方向的四幅移相条纹的组合图案耦合到蓝色B通道,最终得到一幅彩色图,所述彩色图可以用数学公式表示:Four phase-shifting patterns I pxn and I pyn with sinusoidal intensity are designed in the x and y directions respectively. The combined pattern of is coupled to the red R channel, and the combined pattern of the four phase-shifted fringes in the y direction is coupled to the blue B channel, and finally a color image is obtained, which can be represented by a mathematical formula:

Figure BDA0002258230470000021
Figure BDA0002258230470000021

其中A、B分别为正弦条纹的背景强度和调制强度系数;Ix,Iy分别为x,y方向上的强度;IR,IB,IG分别为R,G,B通道中的强度;n为移相步数,且(n=1,2,3,4);Ifx,n=cos(2πffx, nx),Ify,n=cos(2πffy,ny)分别为x,y方向不同的频率载波;ffx,n,ffy,n分别为x,y方向每个载波的设计频率。Among them, A and B are the background intensity and modulation intensity coefficient of the sinusoidal stripes, respectively; I x , I y are the intensities in the x and y directions, respectively; I R , I B , and I G are the intensities in the R, G, and B channels, respectively ; n is the number of phase shift steps, and (n=1,2,3,4); I fx,n =cos(2πf fx, n x),I fy,n =cos(2πf fy,n y) are respectively Different frequency carriers in the x and y directions; f fx,n , f fy,n are the design frequencies of each carrier in the x and y directions, respectively.

优选的,所述Ipxn和Ipyn分别为:Preferably, the I pxn and I pyn are respectively:

Ipxn=1/2+cos(2πfx+n·π/2)I pxn =1/2+cos(2πfx+n·π/2)

Ipyn=1/2+cos(2πfy+n·π/2)I pyn =1/2+cos(2πfy+n·π/2)

其中f是相移图案的频率,x,y分别是投影屏上的像素坐标。where f is the frequency of the phase shift pattern, and x, y are the pixel coordinates on the projection screen, respectively.

优选的,所述图像解调包括步骤:Preferably, the image demodulation comprises the steps of:

对经过待测元件后相机采集的变形彩色图进行颜色通道分离,得到x、y方向上的复合相移条纹;Perform color channel separation on the deformed color image collected by the camera after passing through the component to be tested, and obtain composite phase-shift fringes in the x and y directions;

利用带通滤波器分离出不同载波频率和相移条纹的复合图案,x、y方向上各有四幅;The composite patterns of different carrier frequencies and phase-shift fringes are separated by band-pass filters, and there are four patterns in each of the x and y directions;

根据巴特沃斯滤波器原理,以fn为中心设计带通滤波器。每个频带的截止频率被设计成fn=(fn-1+fn)/2,其中n=1,2,3,…,N和基带通道f0=0,经过带通滤波后的条纹信息为:According to the principle of Butterworth filter, the band-pass filter is designed with f n as the center. The cutoff frequency of each frequency band is designed as f n =(f n-1 +f n )/2, where n = 1, 2, 3, . . . , N and the baseband channel f 0 =0, the The stripe information is:

其中

Figure BDA0002258230470000032
是沿载波方向以fn频率为中心的带通滤波器;ICP(x,y)是复合相移条纹图的强度;I′n(x,y)是相移条纹图的强度;in
Figure BDA0002258230470000032
is the band-pass filter centered at f n frequency along the carrier direction; I CP (x, y) is the intensity of the composite phase-shift fringe pattern; I′ n (x, y) is the intensity of the phase-shift fringe pattern;

经过低通滤波器进行处理,滤除正交方向的载波条纹信息,从而分离出载波频率和相移条纹图:After processing by a low-pass filter, the carrier fringe information in the orthogonal direction is filtered out, so as to separate the carrier frequency and phase shift fringe pattern:

Figure BDA0002258230470000033
Figure BDA0002258230470000033

其中

Figure BDA0002258230470000034
为最终恢复的相移条纹图的强度,h′LP(x)是低通滤波器。in
Figure BDA0002258230470000034
is the intensity of the final recovered phase-shifted fringe pattern, h' LP (x) is the low-pass filter.

与现有技术相比,本发明的有益效果是:此方法打破了传统的测量速度限制,使用一幅彩色图代替多步移相产生的多幅相移图,实现了瞬态测量,降低了成本。在设计编码图案时,现有方法是三步移相产生三幅相移图分别耦合到R,G,B三个通道中,随之出现的问题是,通道之间会互相串扰,但是此方法中,耦合进通道的条纹图已经是调制完毕的复合相移图,这也解决了通道串扰的问题。由于R,G,B的波长是不同的,三幅相移图在耦合进R,G,B通道后,在后续的移相计算中会引入误差,但在此方法中,同一方向的四幅相移图调制成一幅复合相移图进入同一个通道,也就避开了此类误差。本发明为光学偏折检测提供了一种基于频率载波法、相移法和彩色编码法的高精度大动态瞬时测量方法。Compared with the prior art, the present invention has the beneficial effects that: the method breaks the traditional measurement speed limitation, uses one color image to replace the multiple phase-shift images generated by multi-step phase-shifting, realizes the transient measurement, and reduces the cost. When designing the coding pattern, the existing method is three-step phase-shifting to generate three phase-shifting diagrams, which are respectively coupled to the three channels R, G, and B. The problem is that the channels will crosstalk each other, but this method , the fringe pattern coupled into the channel is already a modulated composite phase shift pattern, which also solves the problem of channel crosstalk. Since the wavelengths of R, G, and B are different, after the three phase shift diagrams are coupled into the R, G, and B channels, errors will be introduced in the subsequent phase shift calculation. However, in this method, the four phase shift diagrams in the same direction will This kind of error is avoided by modulating the shift image into a composite phase shift image into the same channel. The invention provides a high-precision large dynamic instantaneous measurement method based on a frequency carrier method, a phase shift method and a color coding method for optical deflection detection.

附图说明Description of drawings

图1为本发明流程图;Fig. 1 is the flow chart of the present invention;

图2为本发明具体实施方式检测光路系统示意图;2 is a schematic diagram of a detection optical path system according to a specific embodiment of the present invention;

图3为本发明用于瞬时偏折测量投影屏投影的频率载波复合色图编码设计示意图;3 is a schematic diagram of the design of the frequency carrier composite color map coding used for instantaneous deflection measurement projection screen projection according to the present invention;

图4为本发明经过待测元件后相机采集变形彩色图的解调过程图;Fig. 4 is the demodulation process diagram of the camera collecting the deformed color image after the component to be tested of the present invention;

图5为本发明解调后的相位恢复图。FIG. 5 is a phase recovery diagram after demodulation of the present invention.

具体实施方式Detailed ways

以下将结合附图所示的具体实施方式对本发明进行详细描述,但这些实施方式并不限制本发明,本领域的普通技术人员根据这些实施方式所做出的结构、方法、或功能上的变换均包含在本发明的保护范围内。The present invention will be described in detail below with reference to the specific embodiments shown in the accompanying drawings, but these embodiments do not limit the present invention, and those of ordinary skill in the art can make structural, method, or functional transformations according to these embodiments. All are included in the protection scope of the present invention.

请参见图1,本发明提出的一种光学偏折瞬态测量方法,包括步骤:Please refer to FIG. 1 , a method for measuring optical deflection transients proposed by the present invention includes the steps:

S1:建立用于透射元件检测的基于光学偏折法的光学检测系统。如图2所示,该光学检测系统包括投影屏1、待测透射元件2、相机3和计算机。投影屏1、待测元件2、相机3按照逆向哈特曼结构进行摆放。相机3前端带有滤光小孔,克服孔径成像像差。投影屏1显示利用正交方向的频率载波条纹进行调制并将不同方向的四幅移相条纹的组合图案耦合到不同的颜色通道,最终形成的一幅彩色图案,彩色图经过待测元件2,相机3能获得待测元件2的完整成像,采集到成像后的变形彩色图。对变形彩色图进行颜色通道分离,得到x、y方向的复合图案,对复合图案进行解调。相机3、投影屏1分别与计算机连接。S1: Establish an optical detection system based on optical deflection method for transmission element detection. As shown in FIG. 2 , the optical detection system includes a projection screen 1 , a transmission element to be measured 2 , a camera 3 and a computer. The projection screen 1, the component to be tested 2, and the camera 3 are placed according to the reverse Hartmann structure. The front end of the camera 3 has a small filter hole to overcome the aperture imaging aberration. Projection screen 1 displays a color pattern finally formed by modulating the frequency carrier stripes in orthogonal directions and coupling the combined pattern of four phase-shifting stripes in different directions to different color channels. 3. The complete imaging of the component to be tested 2 can be obtained, and the deformed color map after imaging can be collected. The color channel separation is performed on the deformed color image to obtain a composite pattern in the x and y directions, and the composite pattern is demodulated. The camera 3 and the projection screen 1 are respectively connected with the computer.

S2:设计所需投影屏投影编码图案。采用以下所述方法:x、y方向上的四步移相条纹图可表示为:S2: Design the projected coding pattern of the required projection screen. The method described below is used: the four-step phase-shifted fringe pattern in the x and y directions can be expressed as:

Ipxn=1/2+cos(2πfx+n·π/2) (1)I pxn =1/2+cos(2πfx+n·π/2) (1)

Ipyn=1/2+cos(2πfy+n·π/2) (2)I pyn =1/2+cos(2πfy+n·π/2) (2)

其中f是相移图案的频率,x,y分别是投影屏上的像素坐标,n是相移步数(n=1,2,3,4)。Where f is the frequency of the phase shift pattern, x and y are the pixel coordinates on the projection screen respectively, and n is the number of phase shift steps (n=1, 2, 3, 4).

经过四步移相后,得到x、y两个方向各四幅移相图案,利用正交方向的频率载波条纹对八幅移相图案进行单独调制,将不同方向的四幅移相条纹的组合图案耦合到不同的颜色通道:After four-step phase shifting, four phase-shifting patterns are obtained in each of the x and y directions. The eight-amplitude phase-shifting patterns are individually modulated by the frequency carrier fringes in the orthogonal directions, and the combined patterns of the four-amplitude phase-shifting fringes in different directions are coupled. to a different color channel:

Figure BDA0002258230470000051
Figure BDA0002258230470000051

其中Ix,Iy分别为x,y方向上的强度;A、B分别为正弦条纹的背景强度和调制强度系数;Ifx,n=cos(2πffx,nx),Ify,n=cos(2πffy,ny)分别为x,y方向不同的频率载波;ffx,n,ffy,n分别为x,y方向每个载波的设计频率。where I x and I y are the intensities in the x and y directions respectively; A and B are the background intensity and modulation intensity coefficient of the sinusoidal stripes respectively; I fx,n =cos(2πf fx,n x), I fy,n = cos(2πf fy,n y) are the different frequency carriers in the x and y directions respectively; f fx,n , f fy,n are the design frequencies of each carrier in the x and y directions, respectively.

将x方向的四幅移相条纹的组合图案耦合到红色R通道,将y方向的四幅移相条纹的组合图案耦合到蓝色B通道,最终得到一幅彩色图。彩色图可以用数学公式表示:The combined pattern of the four phase-shifted fringes in the x-direction is coupled to the red R channel, and the combined pattern of the four phase-shifted fringes in the y-direction is coupled to the blue B channel, and finally a color image is obtained. A color map can be represented by a mathematical formula:

其中,A、B分别为正弦条纹的背景强度和调制强度系数;Ix,Iy分别为x,y方向上的强度;IR,IB,IG分别为R,G,B通道中的强度;n为移相步数,且(n=1,2,3,4);Ifx,n=cos(2πffx,nx),Ify,n=cos(2πffy,ny)分别为x,y方向不同的频率载波;ffx,n,ffy,n分别为x,y方向每个载波的设计频率。Among them, A and B are the background intensity and modulation intensity coefficient of the sinusoidal stripes, respectively; I x , I y are the intensities in the x and y directions, respectively; I R , I B , and I G are the R, G, and B channels, respectively. Intensity; n is the number of phase shift steps, and (n=1, 2, 3, 4); I fx,n =cos(2πf fx,n x), I fy,n =cos(2πf fy,n y) respectively are different frequency carriers in the x and y directions; f fx,n , f fy,n are the design frequencies of each carrier in the x and y directions, respectively.

S3:图像解调。采用以下方法:以所述彩色图投影至所述待测元件,并对经过所述待测元件偏折后相机采集到的变形彩色图进行颜色通道分离,得到x、y两个方向上的复合条纹,利用基于巴特沃斯的带通滤波器分离出不同方向的复合相移图实现了对相移图的均匀滤波。将低通滤波器应用于频率载波方向来获得变形相移图,最终得到八幅移相图案,即x、y方向各采集到四幅相移图。(以x方向为例,y方向同理可得)根据巴特沃斯滤波器原理,以fn为中心设计带通滤波器。每个频带的截止频率被设计成:S3: Image demodulation. The following method is adopted: project the color map to the component to be tested, and separate the color channels of the deformed color map collected by the camera after the component to be tested is deflected to obtain a composite of x and y directions. Stripes, using a Butterworth-based band-pass filter to separate out the composite phase-shift maps in different directions to achieve uniform filtering of the phase-shift maps. A low-pass filter is applied to the direction of the frequency carrier to obtain the deformed phase shift pattern, and finally eight phase shift patterns are obtained, that is, four phase shift patterns are collected in each of the x and y directions. (Take the x direction as an example, the same can be obtained in the y direction) According to the principle of Butterworth filter, a band-pass filter is designed with f n as the center. The cutoff frequency for each band is designed to be:

Figure BDA0002258230470000062
Figure BDA0002258230470000062

其中n=1,2,3,…,N和基带通道f0=0。经过带通滤波后的条纹信息为:where n=1, 2, 3, ..., N and baseband channel f 0 =0. The fringe information after bandpass filtering is:

Figure BDA0002258230470000063
Figure BDA0002258230470000063

其中

Figure BDA0002258230470000064
是沿正交方向以fn频率为中心的带通滤波器;ICP(x,y)是复合相移条纹图的强度;I′n(x,y)是相移条纹图的强度。经过低通滤波器进行处理,滤除掉正交方向的载波条纹信息,从而分离出载波频率和相移条纹图:in
Figure BDA0002258230470000064
is the bandpass filter centered at the f n frequency along the orthogonal direction; I CP (x,y) is the intensity of the composite phase-shift fringe pattern; I′ n (x,y) is the intensity of the phase-shift fringe pattern. After processing by a low-pass filter, the carrier fringe information in the orthogonal direction is filtered out, so as to separate the carrier frequency and the phase shift fringe pattern:

Figure BDA0002258230470000071
Figure BDA0002258230470000071

其中

Figure BDA0002258230470000072
是最终恢复的相移条纹图的强度;h′LP(x)是低通滤波器。in
Figure BDA0002258230470000072
is the intensity of the final recovered phase-shifted fringe pattern; h′ LP (x) is the low-pass filter.

S4:波前重构:通过四步移相算法和相位解包裹算法得到相位信息,从而得到实际光斑分布,将实际光斑分布与光学偏折检测系统中通过光线追迹得到的理想光斑分布进行对比,得出光线像差,所测得的光线像差为待测元件加工误差作用的结果。从而可利用波像差与光线像差的关系得到波前的斜率分布,对波前的斜率进行积分完成波前重构。S4: Wavefront reconstruction: phase information is obtained through the four-step phase shifting algorithm and phase unwrapping algorithm to obtain the actual spot distribution, and the actual spot distribution is compared with the ideal spot distribution obtained by ray tracing in the optical deflection detection system , the light aberration is obtained, and the measured light aberration is the result of the processing error of the component to be tested. Therefore, the relationship between the wave aberration and the light aberration can be used to obtain the slope distribution of the wavefront, and the slope of the wavefront can be integrated to complete the wavefront reconstruction.

实施例Example

图2是一种用于透射元件检测的光学偏折瞬态测量系统的光路布局。待测元件2为直径25.4mm,折射率1.51的凸透镜。设置测量场尺寸为1080×1080像素,x、y方向的载波频率均设置为40、70、100、130个周期所对应的频率,每幅相移条纹图约有12个周期。用于条纹投影的液晶屏像素数为1920(H)×1080(V),对应的像素大小为0.265mm(H)×0.265mm(V),采用三通道8位彩色CCD相机,分辨率1348(H)×1280(V),像素尺寸为3.63μm(H)×3.63μm(V)。Figure 2 is an optical path layout of an optical deflection transient measurement system for transmission element detection. The component to be tested 2 is a convex lens with a diameter of 25.4 mm and a refractive index of 1.51. The size of the measurement field is set to 1080×1080 pixels, and the carrier frequencies in the x and y directions are set to the frequencies corresponding to 40, 70, 100, and 130 cycles, and each phase shift fringe pattern has about 12 cycles. The number of pixels of the LCD screen used for fringe projection is 1920(H)×1080(V), the corresponding pixel size is 0.265mm(H)×0.265mm(V), a three-channel 8-bit color CCD camera is used, and the resolution is 1348 ( H) × 1280 (V), and the pixel size is 3.63 μm (H) × 3.63 μm (V).

步骤1,如图3所示,设计投影屏投影所需编码图案:x,y两个方向各设计四幅强度呈正弦变化的移相图案Ipxn、Ipyn,利用正交方向的频率载波条纹对八幅移相图案进行单独调制,将x方向的四幅移相条纹的组合图案耦合到红色通道,将y方向的四幅移相条纹的组合图案耦合到蓝色通道,最终得到一幅彩色图。彩色图可以用数学公式表示:Step 1, as shown in Figure 3, design the coding pattern required for projection on the projection screen: design four phase-shifting patterns I pxn and I pyn whose intensity varies sinusoidally in the x and y directions, and use the frequency carrier fringes in the orthogonal directions to pair The eight phase-shifting patterns are individually modulated, and the combined pattern of the four phase-shifting fringes in the x-direction is coupled to the red channel, and the combined pattern of the four-phase-shifting fringes in the y-direction is coupled to the blue channel, and finally a color image is obtained. A color map can be represented by a mathematical formula:

Figure BDA0002258230470000073
Figure BDA0002258230470000073

其中,IR,IB,IG分别为R,G,B通道中的光照强度Among them, IR, IB , IG are the light intensity in R , G , B channels respectively

步骤2,如图4所示对经过待测元件后相机采集的变形彩色图进行颜色通道分离,得到x、y两个方向上的复合条纹,利用基于巴特沃斯的带通滤波器分离出不同方向的复合相移图,实现了对相移图的均匀滤波。将低通滤波器应用于频率载波方向来获得相移图,最终得到八幅相移图,即x、y方向各采集到四幅相移图。Step 2: As shown in Figure 4, the color channel separation is performed on the deformed color image collected by the camera after passing through the component to be tested, to obtain composite stripes in the x and y directions, and use a Butterworth-based bandpass filter to separate out the different color channels. The composite phase shift map of the direction, which realizes uniform filtering of the phase shift map. The low-pass filter is applied to the direction of the frequency carrier to obtain the phase shift diagram, and finally eight phase shift diagrams are obtained, that is, four phase shift diagrams are collected in each of the x and y directions.

步骤3,如图5所示得到八幅移相图案后运用四步移相算法和相位解包裹算法得到相位信息。在x和y方向上采集到的包裹相位图如图5中的(a)、(d)以及相应的解包裹后相位如图5中的(b)、(e)。图5中的(c)、(f)则是针对理想相位的对比数据。由此可以得到实际光斑分布,将实际光斑分布与光学偏折检测系统中通过光线追迹得到的理想光斑分布进行对比,得出光线像差,所测得的光线像差待测元件加工误差作用的结果。从而可利用波像差与光线像差的关系得到波前的斜率分布,对波前的斜率进行积分完成波前重构。Step 3, as shown in Figure 5, after obtaining eight phase-shifting patterns, the four-step phase-shifting algorithm and the phase unwrapping algorithm are used to obtain phase information. The wrapped phase maps collected in the x and y directions are shown in (a) and (d) in Figure 5, and the corresponding unwrapped phases are shown in Figure 5 (b) and (e). (c) and (f) in Figure 5 are the comparison data for the ideal phase. From this, the actual light spot distribution can be obtained, and the actual light spot distribution can be compared with the ideal light spot distribution obtained by ray tracing in the optical deflection detection system, and the light aberration can be obtained. the result of. Therefore, the relationship between the wave aberration and the light aberration can be used to obtain the slope distribution of the wavefront, and the slope of the wavefront can be integrated to complete the wavefront reconstruction.

尽管为示例目的,已经公开了本发明的优选实施方式,但是本领域的普通技术人员将意识到,在不脱离由所附的权利要求书公开的本发明的范围和精神的情况下,各种改进、增加以及取代是可能的。Although the preferred embodiments of the present invention have been disclosed for illustrative purposes, those of ordinary skill in the art will appreciate that various Improvements, additions and substitutions are possible.

Claims (4)

1.一种光学偏折瞬态测量方法,其特征在于:1. an optical deflection transient measurement method is characterized in that: 建立基于光学偏折法的光学检测系统,所述光学检测系统包括投影屏、待测元件、相机和计算机;Establish an optical detection system based on the optical deflection method, the optical detection system includes a projection screen, a component to be tested, a camera and a computer; 设计投影屏投影编码图案:利用正交方向的频率载波条纹对四步移相得到的八幅移相图案进行单独调制,将x、y两个不同方向的四幅移相条纹的组合图案耦合到两个不同的颜色通道,得到一幅彩色图;Design the projection coding pattern of the projection screen: use the frequency carrier fringes in orthogonal directions to separately modulate the eight phase-shifting patterns obtained by four-step phase-shifting, and couple the combined pattern of the four phase-shifting fringes in two different directions of x and y to the two. different color channels to get a color image; 图像解调:以所述彩色图投影至所述待测元件,并对经过所述待测元件偏折后相机采集到的变形彩色图进行颜色通道分离,得到x、y方向的复合图案,对所述复合图案进行解调得到x、y两个方向各四幅移相条纹图,一共八幅移相图案;Image demodulation: Project the color map to the component to be tested, and separate the color channels of the deformed color map collected by the camera after the component to be tested is deflected to obtain a composite pattern in the x and y directions. The composite pattern is demodulated to obtain four phase-shift fringe patterns in each of the x and y directions, with a total of eight phase-shift patterns; 波前重构:通过四步移相算法和相位解包裹算法对上述的八幅移相图案进行处理,得到相位信息,从而得到实际光斑分布,将实际光斑分布与所述光学检测系统中通过光线追迹得到的理想光斑分布进行对比,得出光线像差,所述光线像差为所述待测元件加工误差作用的结果;根据波像差与所述光线像差的关系得到波前的斜率分布,对波前的斜率进行积分完成波前重构。Wavefront reconstruction: The above-mentioned eight phase-shifting patterns are processed through the four-step phase-shifting algorithm and the phase unwrapping algorithm to obtain phase information, thereby obtaining the actual spot distribution, and comparing the actual spot distribution with the light passing through the optical detection system. The ideal spot distributions obtained by tracing are compared, and the light aberration is obtained, which is the result of the processing error of the element to be measured; the slope of the wavefront is obtained according to the relationship between the wave aberration and the light aberration distribution, and the slope of the wavefront is integrated to complete the wavefront reconstruction. 2.根据权利要求1所述的一种光学偏折瞬态测量方法,其特征在于:所述设计投影屏投影编码图案包括步骤:2. A kind of optical deflection transient measurement method according to claim 1, is characterized in that: described designing projection screen projection coding pattern comprises the steps: 在x,y两个方向各设计四幅强度呈正弦变化的移相图案Ipxn,Ipyn,利用正交方向的频率载波条纹对八幅移相图案进行单独调制,将x方向的四幅移相条纹的组合图案耦合到红色R通道,将y方向的四幅移相条纹的组合图案耦合到蓝色B通道,最终得到一幅彩色图,所述彩色图可以用数学公式表示:Four phase-shifting patterns I pxn and I pyn with sinusoidal intensity are designed in the x and y directions respectively. The combined pattern of is coupled to the red R channel, and the combined pattern of the four phase-shifted fringes in the y direction is coupled to the blue B channel, and finally a color image is obtained, which can be represented by a mathematical formula:
Figure FDA0002258230460000011
Figure FDA0002258230460000011
其中A、B分别为正弦条纹的背景强度和调制强度系数;Ix,Iy分别为x,y方向上的强度;IR,IB,IG分别为R,G,B通道中的强度;n为移相步数,且n=(1,2,3,4);Ifx,n=cos(2πffx,nx),Ify,n=cos(2πffy,ny)分别为x,y方向不同的频率载波;ffx,n,ffy,n分别为x,y方向每个载波的设计频率。Among them, A and B are the background intensity and modulation intensity coefficient of the sinusoidal stripes, respectively; I x , I y are the intensities in the x and y directions, respectively; I R , I B , and I G are the intensities in the R, G, and B channels, respectively ; n is the number of phase-shifting steps, and n=(1, 2, 3, 4); I fx,n =cos(2πf fx,n x), I fy,n =cos(2πf fy,n y) are respectively Different frequency carriers in the x and y directions; f fx,n , f fy,n are the design frequencies of each carrier in the x and y directions, respectively.
3.根据权利要求2所述的一种光学偏折瞬态测量方法,其特征在于:所述Ipxn和Ipyn分别为:3. a kind of optical deflection transient measurement method according to claim 2, is characterized in that: described I pxn and I pyn are respectively: Ipxn=1/2+cos(2πfx+n·π/2)I pxn =1/2+cos(2πfx+n·π/2) Ipyn=1/2+cos(2πfy+n·π/2)I pyn =1/2+cos(2πfy+n·π/2) 其中f是相移图案的频率,x,y分别是投影屏上的像素坐标。where f is the frequency of the phase shift pattern, and x, y are the pixel coordinates on the projection screen, respectively. 4.根据权利要求1所述的一种光学偏折瞬态测量方法,其特征在于:所述图像解调包括步骤:4. A kind of optical deflection transient measurement method according to claim 1, is characterized in that: described image demodulation comprises the step: 对经过待测元件后相机采集的变形彩色图进行颜色通道分离,得到x、y方向上的复合相移条纹;Perform color channel separation on the deformed color image collected by the camera after passing through the component to be tested, and obtain composite phase-shift fringes in the x and y directions; 利用带通滤波器分离出不同载波频率和相移条纹的复合图案,x、y方向上各有四幅;The composite patterns of different carrier frequencies and phase-shift fringes are separated by band-pass filters, and there are four patterns in each of the x and y directions; 根据巴特沃斯滤波器原理,以fn为中心设计带通滤波器。每个频带的截止频率被设计成fn=(fn-1+fn)/2,其中n=1,2,3,…,N和基带通道f0=0,经过带通滤波后的条纹信息为:According to the principle of Butterworth filter, the band-pass filter is designed with f n as the center. The cutoff frequency of each frequency band is designed as f n =(f n-1 +f n )/2, where n = 1, 2, 3, . . . , N and the baseband channel f 0 =0, the The stripe information is:
Figure FDA0002258230460000021
Figure FDA0002258230460000021
其中是沿载波方向以fn频率为中心的带通滤波器;ICP(x,y)是复合相移条纹图的强度;I′n(x,y)是相移条纹图的强度;in is the band-pass filter centered at f n frequency along the carrier direction; I CP (x, y) is the intensity of the composite phase-shift fringe pattern; I′ n (x, y) is the intensity of the phase-shift fringe pattern; 经过低通滤波器进行处理,滤除正交方向的载波条纹信息,从而分离出载波频率和相移条纹图:After processing by a low-pass filter, the carrier fringe information in the orthogonal direction is filtered out, so as to separate the carrier frequency and phase shift fringe pattern:
Figure FDA0002258230460000023
Figure FDA0002258230460000023
其中
Figure FDA0002258230460000024
为最终恢复的相移条纹图的强度,h′LP(x)是低通滤波器。
in
Figure FDA0002258230460000024
is the intensity of the final recovered phase-shifted fringe pattern, h' LP (x) is the low-pass filter.
CN201911062078.5A 2019-11-01 2019-11-01 Optical deflection transient measurement method Expired - Fee Related CN110793440B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911062078.5A CN110793440B (en) 2019-11-01 2019-11-01 Optical deflection transient measurement method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911062078.5A CN110793440B (en) 2019-11-01 2019-11-01 Optical deflection transient measurement method

Publications (2)

Publication Number Publication Date
CN110793440A true CN110793440A (en) 2020-02-14
CN110793440B CN110793440B (en) 2021-05-18

Family

ID=69440761

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911062078.5A Expired - Fee Related CN110793440B (en) 2019-11-01 2019-11-01 Optical deflection transient measurement method

Country Status (1)

Country Link
CN (1) CN110793440B (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6304325B1 (en) * 1978-01-12 2001-10-16 Raytheon Company Variable shear A. C. interferometer
WO2007019548A2 (en) * 2005-08-08 2007-02-15 Zetetic Institute Apparatus and methods for reduction and compensation of effects of vibrations and of environmental effects in wavefront interferometry
WO2007100144A1 (en) * 2006-02-28 2007-09-07 Canon Kabushiki Kaisha Measurement apparatus, exposure apparatus having the same, and device manufacturing method
US9310309B1 (en) * 2013-02-21 2016-04-12 The United States Of America, As Represented By The Secretary Of The Army Method of sensing acidic/acid-forming and oxidizable gases for use as a residual filter life indicator
CN107543683A (en) * 2017-07-31 2018-01-05 中国计量大学 A high-precision and large dynamic range measurement system and measurement method for transmission element aberration
CN107560564A (en) * 2017-07-28 2018-01-09 中国计量大学 A kind of free form surface detection method and system
CN108225737A (en) * 2017-12-06 2018-06-29 天津大学 A kind of novel optical lens measurement method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6304325B1 (en) * 1978-01-12 2001-10-16 Raytheon Company Variable shear A. C. interferometer
WO2007019548A2 (en) * 2005-08-08 2007-02-15 Zetetic Institute Apparatus and methods for reduction and compensation of effects of vibrations and of environmental effects in wavefront interferometry
WO2007100144A1 (en) * 2006-02-28 2007-09-07 Canon Kabushiki Kaisha Measurement apparatus, exposure apparatus having the same, and device manufacturing method
US9310309B1 (en) * 2013-02-21 2016-04-12 The United States Of America, As Represented By The Secretary Of The Army Method of sensing acidic/acid-forming and oxidizable gases for use as a residual filter life indicator
CN107560564A (en) * 2017-07-28 2018-01-09 中国计量大学 A kind of free form surface detection method and system
CN107543683A (en) * 2017-07-31 2018-01-05 中国计量大学 A high-precision and large dynamic range measurement system and measurement method for transmission element aberration
CN108225737A (en) * 2017-12-06 2018-06-29 天津大学 A kind of novel optical lens measurement method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
DAODANG WANG: "Transmitted wavefront testing with large dynamic", 《JOURNAL OF OPTICS》 *
XIE ZHONGMING: "Instantaneous wavefront", 《SPIE-INT SOC OPTICAL ENGINEERING》 *
徐平等: "基于相位偏折术的大像差透射波前检测 ", 《仪器仪表学报》 *

Also Published As

Publication number Publication date
CN110793440B (en) 2021-05-18

Similar Documents

Publication Publication Date Title
Su et al. Software configurable optical test system: a computerized reverse Hartmann test
CN104949763B (en) A kind of lens wave-front optical aberration measurement method based on inverse Hartmann's principle
CN107560564B (en) A kind of free form surface detection method and system
CN111537203A (en) Optical lens wavefront measuring method using phase measurement deflection technology
CN109781033B (en) Deep ultraviolet structured light precision detection device for transparent material three-dimensional contour reconstruction
CN106989689A (en) The sub-aperture stitching detection technique and device of heavy-calibre planar optical elements face shape
CN108562240B (en) Digital Moire phase-shifting interference measuring method based on two step carrier waves splicing method
CN102494634A (en) Off-axis aspherical mirror detection method based on fringe reflection
CN108955575B (en) A method of corrugated is accurately restored based on single width interference fringe
CN109709786A (en) A super-resolution digital holographic imaging system and imaging method
CN107543683A (en) A high-precision and large dynamic range measurement system and measurement method for transmission element aberration
CN111307063A (en) Method for eliminating symbol ambiguity problem in single-amplitude interference fringe wave surface recovery
CN102506759A (en) Lonky detection method of aspheric surface with heavy calibre
Quiroga et al. Fourier transform method for automatic processing of moiré deflectograms
CN106767391A (en) The sensitivity intensifier and method of four wavefront lateral shearing interference Wavefront sensors
CN106595529A (en) Measurement method and device for large-curvature-radius non-zero-digit interference based on virtual Newton's ring
Chen et al. A flexible lateral chromatic aberration correction method for a color-encoded fringe projection system
CN110793440B (en) Optical deflection transient measurement method
Islas et al. Development of a dynamic interferometer using recycled components based on polarization phase shifting techniques
CN109186497B (en) Spectrum aliasing noise minimizing technology in the digital Moire phase shift interferometry of Newton iteration
CN104463863B (en) The scaling method and system of movement interference field based on the projection of time heterodyne
CN111121661A (en) A N+1 Amplitude Phase Shift Test Algorithm for Narrow-Band Non-Monochromatic Light for Smooth Surface Topography Measurement
CN112762858B (en) Compensation method for phase error in deflection measurement system
Chen et al. Wavefront measurement of a multilens optical system based on phase measuring deflectometry
JP3146590B2 (en) Shape measuring method and shape measuring system

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210518