Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JPH0437362B2 - - Google Patents
[go: Go Back, main page]

JPH0437362B2 - - Google Patents

Info

Publication number
JPH0437362B2
JPH0437362B2 JP57168313A JP16831382A JPH0437362B2 JP H0437362 B2 JPH0437362 B2 JP H0437362B2 JP 57168313 A JP57168313 A JP 57168313A JP 16831382 A JP16831382 A JP 16831382A JP H0437362 B2 JPH0437362 B2 JP H0437362B2
Authority
JP
Japan
Prior art keywords
light
discrete
fourier transform
measured
interference pattern
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
Application number
JP57168313A
Other languages
Japanese (ja)
Other versions
JPS5958305A (en
Inventor
Yoshitada Oshida
Tetsuya Kamioka
Tsutomu Kuze
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.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
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 Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP57168313A priority Critical patent/JPS5958305A/en
Publication of JPS5958305A publication Critical patent/JPS5958305A/en
Publication of JPH0437362B2 publication Critical patent/JPH0437362B2/ja
Granted legal-status Critical Current

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/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/255Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures for measuring radius of curvature

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Instruments For Measurement Of Length By Optical Means (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Description

【発明の詳細な説明】 〔発明の利用分野〕 本発明は、被測定物の表面形状又は表面位置を
測定する方法および装置に関するものである。
DETAILED DESCRIPTION OF THE INVENTION [Field of Application of the Invention] The present invention relates to a method and apparatus for measuring the surface shape or surface position of an object to be measured.

〔従来技術〕[Prior art]

被測定物の表面形状を測定する方法として、従
来、トワイマングリーン型干渉計やフイゾウ型干
渉計により、参照光と被測定物により発生する干
渉縞パターンの形状から被測定物の表面形状を求
める方法が広く行われている。
Conventionally, as a method of measuring the surface shape of a workpiece, a Twyman-Green interferometer or a Huizou interferometer is used to determine the surface shape of the workpiece from the shape of the interference fringe pattern generated by the reference light and the workpiece. The method is widely used.

また、測定精度を向上する方法として、参照光
または物体光の光波位相を変化させた時に生ずる
干渉パターンの強度変化を干渉パターン上のサン
プル点で求め、その強度が最小になる時の物体光
の光波位相の変化量や干渉パターンの強度変化の
正弦波的変化の位相差を各サンプル点で求める方
法が行われている。
In addition, as a method to improve measurement accuracy, the change in the intensity of the interference pattern that occurs when the light wave phase of the reference light or object light is changed is determined at sample points on the interference pattern, and the change in the intensity of the object light when the intensity is the minimum is determined. A method is used in which the amount of change in the phase of a light wave or the phase difference in a sinusoidal change in the intensity of an interference pattern is determined at each sample point.

後者の方法として、複素数高速フーリエ変換
(以下、FFTと略称する。)を用いることが可能
であるが、FFTにより得られる値は離散的周波
数サンプル点における周波数成分(複素数)であ
る。
As the latter method, it is possible to use complex fast Fourier transform (hereinafter abbreviated as FFT), but the values obtained by FFT are frequency components (complex numbers) at discrete frequency sample points.

従つて、入力データが非常に多くならなけれ
ば、離散的周波数サンプル点の間隔は大きくな
り、FFTにより得られた離散的周波数サンプル
点での値の絶対値が最大となるデータのみを用い
て位相値を求めると充分な測定精度が得られな
い。
Therefore, unless the amount of input data becomes very large, the interval between discrete frequency sample points will be large, and the phase will be calculated using only the data for which the absolute value of the value at the discrete frequency sample points obtained by FFT is maximum. When calculating the value, sufficient measurement accuracy cannot be obtained.

これを避けるために入力データ数を多数にする
と、すなわち、光波位相を多数の段階に亘り変化
させ、多数の測定データを用いようとすると、莫
大な測定時間および計算時間を要するのみなら
ず、長時間に亘る測定時間における外乱の影響も
無視することができないという問題がある。
In order to avoid this, if we increase the number of input data, that is, if we try to change the optical wave phase in many stages and use a large number of measurement data, it not only takes a huge amount of measurement time and calculation time, but also takes a long time. There is also a problem in that the influence of disturbances on the measurement time over time cannot be ignored.

〔発明の目的〕[Purpose of the invention]

本発明の目的は、被測定物の表面形状又は表面
位置を短時間で、かつ、非常に高精度に測定する
ことができ、また、外乱条件に比較的強い測定が
できる測定方法および装置を提供することにあ
る。
An object of the present invention is to provide a measuring method and apparatus that can measure the surface shape or surface position of an object to be measured in a short time and with very high precision, and that can perform measurements that are relatively resistant to disturbance conditions. It's about doing.

〔発明の概要〕[Summary of the invention]

本発明は、上記目的を達成するために、可干渉
光源から得られる光を2分し、該2分された光路
の内一方の光路の光路長を変化させて光波位相を
変調せしめ、一方の光を被測定物に照射し、該被
測定物から反射又は透過した光と、上記2分した
他方の光について基準面に照射して得られる参照
光とを両光の波面が測定する範囲でほぼ一致する
ようにして干渉させて干渉パターンを発生させ、
該干渉パターンの空間内の所望の複数サンプル点
におけるN個の強度分布を、上記光路長の時間的
変化に応じて正弦波的に変化するN個の離散的デ
ータとして求めて記憶し、該記憶された各点にお
けるN個の離散的データを離散的フーリエ変換
し、該フーリエ変換で得られた離散的フーリエ変
換データの絶対値中の上記正弦波の周期に相当す
る周波数を与えるフーリエ変換データの最大とな
るサンプル点とその隣りのサンプル点における上
記離散的フーリエ変換データの値を補間して干渉
パターンの上記正弦波的成分の情報を求めること
を特徴とする測定方法である。本発明では、被測
定物からの光および参照光のいずれか一方の光波
位相をN回変化した時に得られるN個の干渉パタ
ーンの強度分布を記録し、この干渉パターン上の
各サンプル点からN個取り出した離散的データを
離散的複素数FFTし、FFTして得られたN個の
離散的複素数フーリエ変換データの絶対値が最大
となる離散的周波数と、その周波数の前後の離散
的周波数に対応する複素数フーリエ変換データを
補間し、それぞれのサンプル点での干渉パターン
の位相を求め、それにより被測定物の表面形状を
測定するようにしたことに特徴がある。
In order to achieve the above object, the present invention divides light obtained from a coherent light source into two, changes the optical path length of one of the two divided optical paths, modulates the optical wave phase, and modulates the optical phase of one of the two divided optical paths. Light is irradiated onto the object to be measured, and the light reflected or transmitted from the object to be measured and the reference light obtained by irradiating the reference surface with the other light divided into two are measured within the range in which the wavefronts of both lights are measured. Interfering so that they almost match and generating an interference pattern,
N intensity distributions at a plurality of desired sample points in the space of the interference pattern are determined and stored as N discrete data that changes sinusoidally in accordance with the temporal change in the optical path length; The N pieces of discrete data at each point are subjected to discrete Fourier transform, and the Fourier transform data gives a frequency corresponding to the period of the sine wave in the absolute value of the discrete Fourier transform data obtained by the Fourier transform. This measurement method is characterized in that information on the sinusoidal component of the interference pattern is obtained by interpolating the values of the discrete Fourier transform data at the maximum sample point and the sample points adjacent thereto. In the present invention, the intensity distribution of N interference patterns obtained when the optical wave phase of either the light from the object to be measured or the reference light is changed N times is recorded, and from each sample point on this interference pattern Perform discrete complex FFT on the extracted discrete data, and correspond to the discrete frequency where the absolute value of the N pieces of discrete complex Fourier transform data obtained by FFT is the maximum, and the discrete frequencies before and after that frequency. The feature is that the complex number Fourier transform data obtained by the method is interpolated, the phase of the interference pattern at each sample point is determined, and the surface shape of the object to be measured is thereby measured.

〔発明の実施例〕[Embodiments of the invention]

以下、本発明の実施例を図面により詳細に説明
する。
Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.

第1図は本発明による表面形状測定装置の一実
施例を示すものである。
FIG. 1 shows an embodiment of a surface profile measuring device according to the present invention.

図において、レーザ等の可干渉性光源1を出射
した光を高速シヤツタ2によりオン−オフする。
高速シヤツタ2を通した光をビーム拡大光学系3
により所望のビーム径に拡大し、ビームスプリツ
タ4により2つの光路5および6に分離する。光
路5には、被測定物7が設置されており、被測定
物7が球面からなる場合にはフオーカスレンズ8
が挿入される。このフオーカスレンズ8は入射す
る拡大平行ビームを一点に集束せしめるもので、
この集束点と球面鏡としての被測定物7の曲率中
心が一致するように、被測定物7を設置すれば、
被測定物7で反射した光は再びフオーカスレンズ
8を通過し平行ビームとなり、ビームスプリツタ
4に戻つて来る。
In the figure, light emitted from a coherent light source 1 such as a laser is turned on and off by a high-speed shutter 2.
Beam expanding optical system 3 for light passing through high-speed shutter 2
The beam is expanded to a desired beam diameter, and is separated into two optical paths 5 and 6 by a beam splitter 4. An object to be measured 7 is installed in the optical path 5, and when the object to be measured 7 has a spherical surface, a focus lens 8 is installed.
is inserted. This focus lens 8 focuses the incident expanded parallel beam to one point.
If the object to be measured 7 is installed so that this focal point and the center of curvature of the object to be measured 7 as a spherical mirror coincide,
The light reflected by the object to be measured 7 passes through the focus lens 8 again, becomes a parallel beam, and returns to the beam splitter 4.

ビームスプリツタ4で分離された他方のビーム
の光路6は参照光路となるが、この参照光路には
楔ガラス9が挿入されており、この楔ガラス9は
駆動機構10により一定ピツチずつ移動される。
楔ガラス9を通過し、折返しミラー11で反射さ
れ、楔ガラス9を再び通過したビームは、ビーム
スプリツタ4に戻る。
The optical path 6 of the other beam separated by the beam splitter 4 serves as a reference optical path, and a wedge glass 9 is inserted into this reference optical path, and this wedge glass 9 is moved by a constant pitch by a drive mechanism 10. .
The beam that passes through the wedge glass 9, is reflected by the reflection mirror 11, passes through the wedge glass 9 again, and returns to the beam splitter 4.

被測定物7から戻つたビームおよび参照光路か
らのビーム変換は干渉パターン発生光路12で干
渉する。この干渉パターンは結像レンズ13によ
り、撮像装置14の撮像面15上に投影される。
この結像レンズ13は、ほゞ被測定物7の像を撮
像面15に結像する関係になつている。
The beam returned from the object to be measured 7 and the beam converted from the reference optical path interfere in the interference pattern generation optical path 12. This interference pattern is projected onto the imaging surface 15 of the imaging device 14 by the imaging lens 13.
The imaging lens 13 is arranged to substantially form an image of the object to be measured 7 on an imaging surface 15 .

撮像面15は、第2図に示すように、実効的に
L×L個に分割された撮像サンプル点を有し、干
渉パターン30は一定ピツチでサンプルされ、各
撮像サンプル点における干渉パターン強度が検出
され、制御回路16に送られる。
As shown in FIG. 2, the imaging surface 15 has imaging sample points that are effectively divided into L×L pieces, and the interference pattern 30 is sampled at a constant pitch, and the interference pattern intensity at each imaging sample point is It is detected and sent to the control circuit 16.

制御回路16では、楔ガラス9を一定ピツチで
間欠的に移動しては高速シヤツタ2を開放し、撮
像装置14で得られた干渉パターン強度信号を取
り込む。このようにして、N回の楔ガラス9の位
置移動により、N個の干渉パターン強度信号を取
り込み、それをアナログ/デイジタル変換し、コ
ンピユータ17を通して外部メモリ18に蓄積す
る。
The control circuit 16 moves the wedge glass 9 intermittently at a constant pitch, opens the high-speed shutter 2, and captures the interference pattern intensity signal obtained by the imaging device 14. In this way, by moving the wedge glass 9 N times, N interference pattern intensity signals are captured, converted into analog/digital data, and stored in the external memory 18 via the computer 17.

第ti回目の楔ガラス9の位置に対して、撮像面
15の(l、m)番地の干渉パターン強度を
Iln(ti)とすると、楔ガラス9は一定ピツチで
移動するため、例えば、第2図の撮像面上のAお
よびB点の強度信号はそれぞれ第3図のIAおよび
IBのように変化する。撮像面15上の干渉パター
ン発生部の全サンプル点について、同様の強度変
化が得られ、この正弦波状強度変化の相対的な位
相差が、被測定物の表面形状を表わしている。従
つて、正弦波状強度変化の位相を求めればよい。
For the t i- th position of the wedge glass 9, the interference pattern intensity at address (l, m) on the imaging surface 15 is calculated.
If I l , n (t i ), the wedge glass 9 moves at a constant pitch, so for example, the intensity signals at points A and B on the imaging plane in FIG. 2 are I A and I in FIG. 3, respectively.
Change like I B. A similar intensity change is obtained for all sample points of the interference pattern generating section on the imaging surface 15, and the relative phase difference of this sinusoidal intensity change represents the surface shape of the object to be measured. Therefore, it is only necessary to find the phase of the sinusoidal intensity change.

外部メモリ18に格納された多数のデータか
ら、(l、m)番地のデータのみをN個コンピユ
ータ17に取り出し、これをFFT演算回路19、
例えば、アレイプロセツサに入力し、複素数
FFT演算を行う。
From a large amount of data stored in the external memory 18, only N pieces of data at addresses (l, m) are taken out to the computer 17, and are sent to the FFT calculation circuit 19,
For example, if you input a complex number into an array processor,
Perform FFT calculation.

すなわち、次式に示す離散的フーリエ変換を行
う。
That is, the discrete Fourier transform shown in the following equation is performed.

Ol,n(fk)=N-1ti=0 Il,n(ti)expj2πtifk/N ………(1) (j:虚数、ti,fk=0,1,2,……,N−
1) 離散的フーリエ変換で得られた結果は、第4図
aあるいは第4図bのようになる。
O l,n (f k )= N-1ti=0 I l,n (ti) exp j2πtifk/N ………(1) (j: imaginary number, ti, f k =0,1,2,… ..., N-
1) The results obtained by discrete Fourier transform are as shown in Figure 4a or Figure 4b.

第4図aは、楔ガラス9が一定回数間欠送りさ
れた時、元の干渉パターン強度を再現している場
合である。すなわち、楔ガラスによる1回の間欠
送りによる位相変調量△φが2π/L(L:整数)
で表わされる場合である。この場合には、fp,fk
fN-1-kを除き、全てのサンプル点で、離散的デー
タO(f)の絶対値|O(f)|は0に近い値となり、複
素数Ofkの位相、すなわち、データOfkの実数部
をデータO(fk)の虚数部で割つたもののアーク
タンジエント(arctangent)は求めんとする表面
形状位相となる。
FIG. 4a shows a case where the original interference pattern intensity is reproduced when the wedge glass 9 is intermittently fed a certain number of times. In other words, the amount of phase modulation △φ due to one time of intermittent feeding by the wedge glass is 2π/L (L: integer)
This is the case expressed as. In this case, f p , f k ,
At all sample points except f N-1-k , the absolute value |O(f)| of the discrete data O(f) is close to 0, and the phase of the complex number Of k , that is, the phase of the data Of k The arctangent of the real part divided by the imaginary part of the data O( fk ) becomes the surface shape phase to be sought.

しかるに、一般には、上記のLは整数にならな
い。この場合は、fkおよびfp以外の点でも|O(f)
|は0とならず、干渉パターンの正弦波的変化の
周波数に相当する値は、第4図bに示すように、
FFTの離散的データO(f)の絶対値が最大となる
点fkからずれた点fk′(fk′は整数でない。)とな
る。
However, in general, the above L is not an integer. In this case, |O(f) at points other than f k and f p
| is not 0, and the value corresponding to the frequency of the sinusoidal change in the interference pattern is as shown in Figure 4b.
The point f k ' (f k ' is not an integer) is shifted from the point f k where the absolute value of the FFT discrete data O(f ) is maximum.

従つて、この点fk′を求め、その点におけるデ
ータO(f)の位相を求める必要がある。楔ガラス9
の1回の移動により生ずる光波位相変化△φを
k0/Nとする。こゝで、Nは撮像する干渉パター
ン数であり、k0は実数であり、次の式を満足する
ものとする。
Therefore, it is necessary to find this point f k ' and find the phase of data O(f) at that point. wedge glass 9
The light wave phase change △φ caused by one movement of
Let 2π k0/N . Here, N is the number of interference patterns to be imaged, k 0 is a real number, and it is assumed that the following equation is satisfied.

k0=k1+△k ………(2) (k1:整数、|△k|≦0.5) また、(l,m)番地の位相をφ(l,m)とす
れば、干渉パターンの強度Il,n(ti)は次式で与え
られる。
k 0 = k 1 +△k ………(2) (k 1 : integer, |△k|≦0.5) Also, if the phase of address (l, m) is φ(l, m), the interference pattern is The intensity I l,n (t i ) is given by the following equation.

Il,n(ti) =a0+a1cos{2πk0/Nti+φ(l,m)}………
(3) ti=1,2,3,……,N,a0≧a1 被測定物からの光と参照光とが等しい時のみ、
a0=a1となるが、一般には、a0>a1となるように
設計される。
I l,n (t i ) =a 0 +a 1 cos {2πk 0 /Nt i +φ(l, m)}...
(3) t i = 1, 2, 3, ..., N, a 0 ≧a 1Only when the light from the object to be measured and the reference light are equal,
Although a 0 = a 1 , it is generally designed so that a 0 > a 1 .

FFTを実行した結果は、(3)式を(1)式に代入し
た結果と一致するから、次の式が得られる。
The result of executing FFT matches the result of substituting equation (3) into equation (1), so the following equation is obtained.

Ol,n(fk)=a0exp(iπfkN−1/N)sinπfk/sinπ
fk/N+a1〔expi{φ(l,m) +π(fk+k0)N−1/N}〕sinπ(fk+k0)/sin
πfk+k0/N+expi{−φ(l,m)+π(fk−k0)N
−1/N}sinπ(fk−k0)/sinπfk−k0/N ………(4) この(4)式の第1項は、Nが大きい時、デルタ関
数δ(fk)に比例し、Nの大小に拘わらず、fk
0でN、他の整数fk(0≦fk≦N−1)に対して
は0となる。第2項および第3項はNが大きい
時、デルタ関数δ(fk+k0)およびδ(fk−k0)に
それぞれ比例し、Nの大小に拘わらず、fk=±k0
でNとなる。従つて、k0が整数の時は、第4図a
のように、k0が整数でない時は第4図bのように
なる。第4図bに示すように、FFTにより得ら
れる結果は、fkが整数における値であるから、fk
=±k0(k0≠整数)におけるデータOl,nf(k)の値は
直接出力されない。
O l,n (f k )=a 0 exp (iπf k N-1/N) sinπf k /sinπ
f k /N+a 1 [expi{φ(l,m) +π(f k +k 0 )N-1/N}] sinπ(f k +k 0 )/sin
πf k +k 0 /N+expi{−φ(l,m)+π(f k −k 0 )N
−1/N}sinπ(f k −k 0 )/sinπf k −k 0 /N ………(4) When N is large, the first term of this equation (4) is the delta function δ(f k ) , regardless of the size of N, f k =
0 for N, and 0 for other integers f k (0≦f k ≦N-1). When N is large, the second and third terms are proportional to the delta functions δ(f k +k 0 ) and δ(f k −k 0 ), respectively, and regardless of the size of N, f k =±k 0
becomes N. Therefore, when k 0 is an integer, Fig. 4a
When k 0 is not an integer, the result will be as shown in Figure 4b. As shown in Figure 4b, the result obtained by FFT is that f k is an integer value, so f k
The value of data O l,nf(k) at =±k 0 (k 0 ≠ integer) is not directly output.

しかるに、(l,m)番地におけるfk=k0での
データOl,n(fk)の位相(位相差)が表面形状を表
わしているため、コンピユータ17では、以下の
演算処理を実行する。
However, since the phase (phase difference) of the data O l,n (f k ) at f k = k 0 at address (l, m) represents the surface shape, the computer 17 executes the following calculation process. do.

(1) 楔ガラスの移動量は、予め、ある程度正確に
分つているため、(2)式の整数k1は既知である。
そこで、このk1に対し、△kを−0.5から0.5ま
で一定ステツプで変化させた時のデータOl,n(fo)
の位相Ψ(△k,φ)を次式により求める。
(1) Since the amount of movement of the wedge glass is known in advance with some degree of accuracy, the integer k 1 in equation (2) is known.
Therefore, for this k 1 , the data O l,n(fo) when △k is changed from -0.5 to 0.5 in constant steps is
The phase Ψ(△k, φ) of is determined by the following equation.

Ψ(△k,φ) =tan-1〔In{Ol,n(k1+△k)}/Re{Ol,n(k1
+△k)}〕………(5) (2) 同様にして、△kとφ(l,m)を上述した
ように変化させた時の、次の(6)式で与えられる
Rをテーブルにする。
Ψ(△k,φ) =tan -1 [I n {O l,n (k 1 +△k)}/Re{O l,n (k 1
+△k)}]……(5) (2) Similarly, when △k and φ(l, m) are changed as described above, R given by the following equation (6) is Make it a table.

R(△k,φ) =|Ol,n(k0+S)|/|Ol,n(k0)| ………(6) 但し、こゝで、Sは次の式の条件を満たすも
のとする。
R(△k,φ) =|O l,n (k 0 +S)|/|O l,n (k 0 )| ………(6) However, here, S satisfies the condition of the following equation. shall be met.

|Ol,n(k0+S)|≧|Ol,n(k0−S)| |S|=1 ………(7) (3) 測定されたN個の干渉パターンデータの内の
(l,m)番地のN個のデータをFFTした結果
Ol,n(fk)の絶対値が最大となるfkをfk≠0、fk
≦N/2の範囲で見つけ、その時のfkをk0とする。
|O l,n (k 0 +S)|≧|O l,n (k 0 −S)||S|=1 ………(7) (3) Among the N measured interference pattern data Result of FFT of N pieces of data at address (l, m)
The absolute value of O l,n (f k ) is maximum, f k ≠ 0, f k
Find it in the range ≦N/2, and let f k at that time be k 0 .

このFFTの結果から、fk=k0におけるOl,n
(fk)の位相Ψ(k0)と次式を求める。
From this FFT result, O l,n at f k = k 0
Find the phase Ψ(k 0 ) of (f k ) and the following equation.

R′=|Ol,n(k0+S)|/|Ol,n(k0)|………(
8) 但し、Sは(7)式で与えられる。
R'=|O l,n (k 0 +S)|/|O l,n (k 0 )|......(
8) However, S is given by equation (7).

すなわち、k0+Sはk0の隣りで、2番目に絶
対値が大きなFFT出力のサンプル点である。
That is, k 0 +S is the sample point of the FFT output that is next to k 0 and has the second largest absolute value.

予め作成された、第5図のテーブルを用い、
φ(k0)とR′に一致する△kとφ(l0)をテーブ
ル内から求める処理を行えば、その時のφ(l0
が求める位相となる。
Using the previously created table shown in Figure 5,
If we calculate △k and φ(l 0 ) that match φ(k 0 ) and R′ from the table, then φ(l 0 )
becomes the desired phase.

なお、上述した(3)項の処理を行うに際して、
テーブル内のΨ(△k,φ)とR(△k,φ)の
値が内挿すれば、更に高精度に値を求めること
が可能となる。
In addition, when carrying out the processing in paragraph (3) above,
If the values of Ψ(Δk, φ) and R(Δk, φ) in the table are interpolated, it becomes possible to obtain the values with even higher precision.

この(1)〜(3)の処理を被測定物の干渉パターンが
表われている全番地に亘り行うことにより、被測
定物の形状を求めることができる。
By performing the processes (1) to (3) over all addresses where the interference pattern of the object to be measured appears, the shape of the object to be measured can be determined.

このようにして得られた被測定物の表面形状
は、第1図に示す測定光学系そのものが完全な平
面鏡や安全なビームスプリツタ、フオーカスレン
ズを具備していないことにより、必らずしも正確
な表面形状測定結果を表わしていない。
The surface shape of the object to be measured obtained in this way is not always accurate due to the fact that the measurement optical system shown in Figure 1 itself is not equipped with a perfect plane mirror, a safe beam splitter, or a focus lens. However, it does not represent accurate surface profile measurement results.

このような場合には、被測定物の測定に先が
け、非常に精度の高い原器(第1図では球面原
器)を被測定物として用いて、上述したと同様な
方法により測定する。その測定結果は測定光学系
そのものの不完全性(歪)を表わしているため、
この測定結果を以後の任意の球面測定の補正値と
して用いれば、球面原器の精度を測定絶対精度と
して任意の球面を高精度に測定することができ
る。
In such a case, prior to measuring the object to be measured, a highly accurate prototype (a spherical prototype in FIG. 1) is used as the object to be measured, and the measurement is carried out in the same manner as described above. The measurement results represent imperfections (distortions) in the measurement optical system itself, so
If this measurement result is used as a correction value for any subsequent spherical surface measurement, any spherical surface can be measured with high precision using the accuracy of the spherical prototype as the absolute measurement accuracy.

なお、第1図の例では、球面を被測定物として
いるが、フオーカスレンズを除けば、平面の被測
定物を測定することができる。
In the example shown in FIG. 1, the object to be measured is a spherical surface, but if the focus lens is removed, it is possible to measure a flat object.

また、第1図の例では、反射型の被測定物を示
したが、マツハ・ゼンダー(Mach−Zender)型
干渉計を用いれば、レンズや透過型光学部品の特
性も同様に測定することが可能となる。
In addition, although the example in Figure 1 shows a reflective object to be measured, if a Mach-Zender interferometer is used, the characteristics of lenses and transmission optical components can also be measured in the same way. It becomes possible.

上述した実施例から解るように、比較的少ない
N個の干渉パターンの強度データを離散的複素数
高速フーリエ変換して得られるデータから精度の
高い測定結果が得られるため、比較的短時間の測
定および処理時間で、非常に精度の高い表面形状
測定が可能となるばかりか、外乱の影響も無視で
きる。
As can be seen from the above embodiments, highly accurate measurement results can be obtained from data obtained by performing discrete complex fast Fourier transform on the intensity data of a relatively small number of N interference patterns. Not only is it possible to measure the surface shape with extremely high accuracy within a short processing time, but the effects of external disturbances can also be ignored.

また、楔ガラス等の位相変調器の光波変調の光
波変調ステツプが特定の値でなくても、高精度測
定を可能にすることである。すなわち、従来方法
では数十分のλ程度の精度であつたが、本発明で
は、数百分のλの精度を実現することができる。
Another object of the present invention is to enable highly accurate measurement even if the light wave modulation step of the light wave modulation of a phase modulator such as a wedge glass is not a specific value. That is, while the conventional method had an accuracy of several tens of λ, the present invention can achieve an accuracy of several hundred λ.

さらに、このような高精度測定のためには、従
来方法は、測定光学系や被測定物の固定の安定性
や、被測定物の光路および参照光路の空気のゆら
ぎが非常に安定していることが不可欠であつた
が、上述した実施例では、測定時間内に線型に変
化するものに対しては安定性をそれぞれ要求され
ないため、測定光学系に要求される条件は緩くて
も高精度測定が可能となる。
Furthermore, in order to achieve such high-precision measurements, conventional methods require very stable fixation of the measurement optical system and object to be measured, as well as air fluctuations in the optical path of the object to be measured and the reference optical path. However, in the example described above, stability is not required for things that change linearly within the measurement time, so even if the conditions required for the measurement optical system are loose, high precision measurement is possible. becomes possible.

〔発明の効果〕〔Effect of the invention〕

以上述べたように、本発明によれば、被測定物
の表面形状又は表面位置を短時間に、かつ、高精
度に測定することができ、また、外乱条件に比較
的強い表面形状又は表面位置測定ができる。
As described above, according to the present invention, the surface shape or surface position of the object to be measured can be measured in a short time and with high precision, and the surface shape or surface position can be relatively resistant to disturbance conditions. Can be measured.

【図面の簡単な説明】[Brief explanation of drawings]

図はいずれも本発明の実施例に係るもので、第
1図は測定装置の一実施例の構成図、第2図は第
1図の撮像装置の撮像面上の干渉パターンと撮像
サンプル点の配置を示す図、第3図は干渉パター
ンのサンプル点上の強度変化を示す図、第4図
a,bは干渉パターンのサンプル点上の強度の
FFT結果を示す図、第5図は位相値を求めるた
めのテーブルを示す図である。 1…可干渉性光源、8…被測定物、9…楔ガラ
ス、14…撮像装置、15…撮像面、16…制御
回路、17…コンピユータ。
The figures are all related to embodiments of the present invention. Fig. 1 is a configuration diagram of an embodiment of the measuring device, and Fig. 2 shows the interference pattern on the imaging surface of the imaging device shown in Fig. 1 and the imaging sample point. Figure 3 is a diagram showing the arrangement, Figure 3 is a diagram showing intensity changes on sample points of the interference pattern, and Figures 4 a and b are diagrams showing intensity changes on sample points of the interference pattern.
A diagram showing the FFT results, and FIG. 5 is a diagram showing a table for determining the phase value. DESCRIPTION OF SYMBOLS 1... Coherence light source, 8... Measured object, 9... Wedge glass, 14... Imaging device, 15... Imaging surface, 16... Control circuit, 17... Computer.

Claims (1)

【特許請求の範囲】 1 可干渉光源から得られる光を2分し、該2分
された光路の内一方の光路の光路長を変化させて
光波位相を変調せしめ、一方の光を被測定物に照
射し、該被測定物から反射又は透過した光と、上
記2分した他方の光について基準面に照射して得
られる参照光とを両光の波面が測定する範囲でほ
ぼ一致するようにして干渉させて干渉パターンを
発生させ、該干渉パターンの空間内の所望の複数
サンプル点におけるN個の強度分布を、上記光路
長の時間的変化に応じて正弦波的に変化するN個
の離散的データとして求めて記憶し、該記憶され
た各点におけるN個の離散的データを離散的フー
リエ変換し、該フーリエ変換で得られた離散的フ
ーリエ変換データの絶対値中の上記正弦波の周期
に相当する周波数を与えるフーリエ変換データの
最大となるサンプル点とその隣りのサンプル点に
おける上記離散的フーリエ変換データの値を補間
して干渉パターンの上記正弦波的成分の情報を求
めることを特徴とする測定方法。 2 可干渉光源から得られる光を2分し、一方の
光を被測定物に照射し、該被測定物から反射又は
透過した光と、上記2分した他方の光について該
反射又は透過した光と干渉させる参照光として得
られる他方の光との干渉パターンのN個の強度分
布を、上記被測定物からの光及び参照光のいずれ
か一方の光波位相をN回変化させて求めて記憶
し、干渉パターン上の1つのサンプル点からN個
取り出した離散的データを離散的複素数フーリエ
変換し、該離散的複素数フーリエ変換で得られた
N個の離散的複素数フーリエ変換データの絶対値
が最大となる離散的周波数と該周波数に隣接する
離散的周波数に対応する離散的複素数フーリエ変
換データを補間して上記サンプル点での干渉パタ
ーンの位相を求めることを特徴とする測定方法。 3 可干渉性光源と、該光源からの光を2つの光
路に分離する分離手段と、該2つの光路の内一方
の光路の光路長を変化させて光波位相を変調せし
める変調手段と、上記2つの光路の一方を被測定
物に照射させる照射手段と、上記2分した他方の
光について基準面に照射して参照光を得る参照光
手段と、上記照射手段によつて照射された被測定
物から反射又は透過した光と、上記参照光手段か
ら得られる参照光とを1つの光路に導いて両光の
波面が測定する範囲でほぼ一致するようにして干
渉させて干渉パターンを発生させる発生手段と、
該発生手段で発生された干渉パターンの空間内の
所望の複数サンプル点におけるN個の強度分布
を、上記変調手段で変調された光路長の時間的変
化に応じて正弦波的に変化するN個の離散的デー
タとして検出する検出手段と、該検出手段でN個
の離散的データとして検出されたN個の干渉パタ
ーンの強度分布を記憶する記憶手段と、該記憶手
段に記憶された各点におけるN個の離散的データ
を離散的フーリエ変換するフーリエ変換手段と、
該フーリエ変換手段で得られた離散的フーリエ変
換データの絶対値中の上記正弦波の周期に相当す
る周波数を与えるフーリエ変換データの最大とな
るサンプル点とその隣りのサンプル点における上
記離散的フーリエ変換データの値を補間して干渉
パターンの上記正弦波的成分の情報を求める測定
手段とを備えたことを特徴とする測定装置。 4 可干渉性光源と、該光源からの光を2つの光
路に分離する分離手段と、該2つの光路の一方を
被測定物に照射させる照射手段と、上記2つの光
路の他方の光の光波位相を変調する変調手段と、
上記2つの光路を1つの光路に導いて干渉パター
ンを発生する発生手段と、該発生手段で発生され
た干渉パターンの強度をサンプル点で検出する検
出手段と、該検出手段で検出された、上記光波位
相をN回変調したときのN個の干渉パターンの強
度分布を記憶する記憶手段と、該記憶手段に記憶
されたパターンから、所定サンプル点の強度信号
をN個取り出し、離散的複素数フーリエ変換する
フーリエ変換手段と、該フーリエ変換手段で得ら
れた離散的複素数フーリエ変換データの絶対値が
最大となる離散的周波数と該周波数に隣接する離
散的周波数に対応するフーリエ変換データを補間
して上記サンプル点での干渉パターンの位相を求
める測定手段とを備えたことを特徴とする測定装
置。 5 上記変調手段が、楔形ガラスと、該ガラスを
一定微少量ずつ移動させる手段を有することを特
徴とする特許請求の範囲第4項記載の測定装置。
[Claims] 1. Light obtained from a coherent light source is divided into two, the optical path length of one of the two divided optical paths is changed to modulate the optical wave phase, and one of the light is directed to the object to be measured. The light reflected or transmitted from the object to be measured and the reference light obtained by irradiating the reference surface with the other divided light are made so that the wavefronts of both lights almost match within the range to be measured. to generate an interference pattern, and N intensity distributions at a plurality of desired sample points in the space of the interference pattern are determined by N discrete intensity distributions that vary sinusoidally in accordance with temporal changes in the optical path length. N pieces of discrete data at each of the stored points are subjected to discrete Fourier transform, and the period of the sine wave in the absolute value of the discrete Fourier transform data obtained by the Fourier transform is Information on the sinusoidal component of the interference pattern is obtained by interpolating the values of the discrete Fourier transform data at the maximum sample point of the Fourier transform data giving a frequency corresponding to the frequency and the sample point adjacent thereto. measurement method. 2 The light obtained from the coherent light source is divided into two, one of the lights is irradiated onto the object to be measured, and the light that is reflected or transmitted from the object to be measured and the light that is reflected or transmitted from the other of the above-mentioned two halves. N intensity distributions of interference patterns with the other light obtained as a reference light to be interfered with are obtained by changing the light wave phase of either the light from the object to be measured or the reference light N times and stored. , discrete complex Fourier transform is performed on N pieces of discrete data taken from one sample point on the interference pattern, and the absolute value of the N pieces of discrete complex Fourier transform data obtained by the discrete complex Fourier transform is the maximum. A measurement method characterized by interpolating discrete complex Fourier transform data corresponding to a discrete frequency and a discrete frequency adjacent to the frequency to determine the phase of the interference pattern at the sample point. 3. A coherent light source, a separation means for separating the light from the light source into two optical paths, a modulation means for modulating the optical wave phase by changing the optical path length of one of the two optical paths, and 2. irradiation means for irradiating the object to be measured with one of the two optical paths; reference light means for irradiating the other divided light onto a reference surface to obtain a reference light; and the object to be measured irradiated by the irradiation means. generation means for generating an interference pattern by guiding the light reflected or transmitted from the reference light and the reference light obtained from the reference light means into one optical path so that the wavefronts of both lights almost match within the measurement range, causing interference; and,
N intensity distributions at a plurality of desired sample points in the space of the interference pattern generated by the generation means are changed sinusoidally in accordance with temporal changes in the optical path length modulated by the modulation means. a detection means for detecting as discrete data, a storage means for storing intensity distributions of N interference patterns detected as N discrete data by the detection means; Fourier transform means for performing discrete Fourier transform on N pieces of discrete data;
The discrete Fourier transform at the sample point at which the Fourier transform data has a maximum and the sample point adjacent thereto that gives a frequency corresponding to the period of the sine wave in the absolute value of the discrete Fourier transform data obtained by the Fourier transform means. A measuring device comprising: measuring means for interpolating data values to obtain information on the sinusoidal component of the interference pattern. 4 A coherent light source, a separation means for separating the light from the light source into two optical paths, an irradiation means for irradiating one of the two optical paths onto the object to be measured, and a light wave of the other light of the two optical paths. a modulating means for modulating the phase;
a generating means for generating an interference pattern by guiding the two optical paths into one optical path; a detecting means for detecting the intensity of the interference pattern generated by the generating means at a sample point; A storage means for storing the intensity distribution of N interference patterns when the light wave phase is modulated N times, and N intensity signals at predetermined sample points are taken out from the patterns stored in the storage means and subjected to discrete complex Fourier transformation. and a Fourier transform means for interpolating the Fourier transform data corresponding to a discrete frequency at which the absolute value of the discrete complex Fourier transform data obtained by the Fourier transform means is maximum and a discrete frequency adjacent to the frequency. A measuring device comprising: measuring means for determining the phase of an interference pattern at a sample point. 5. The measuring device according to claim 4, wherein the modulating means includes a wedge-shaped glass and a means for moving the glass by a constant minute amount.
JP57168313A 1982-09-29 1982-09-29 Measuring method and device Granted JPS5958305A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP57168313A JPS5958305A (en) 1982-09-29 1982-09-29 Measuring method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP57168313A JPS5958305A (en) 1982-09-29 1982-09-29 Measuring method and device

Publications (2)

Publication Number Publication Date
JPS5958305A JPS5958305A (en) 1984-04-04
JPH0437362B2 true JPH0437362B2 (en) 1992-06-19

Family

ID=15865710

Family Applications (1)

Application Number Title Priority Date Filing Date
JP57168313A Granted JPS5958305A (en) 1982-09-29 1982-09-29 Measuring method and device

Country Status (1)

Country Link
JP (1) JPS5958305A (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5990009A (en) * 1982-11-15 1984-05-24 Rikagaku Kenkyusho Method and device for measuring shape of body
JPS61294327A (en) * 1985-06-21 1986-12-25 Agency Of Ind Science & Technol Interference method and interferrometer for measuring surface shape of parabolic mirror
JPS62214309A (en) * 1986-03-17 1987-09-21 Tokyo Seimitsu Co Ltd Measuring instrument for surface roughness and shape
JPH02238306A (en) * 1989-03-13 1990-09-20 Ricoh Co Ltd Apparatus for measuring fine displacement
JP2002048506A (en) 2000-08-04 2002-02-15 Matsushita Electric Ind Co Ltd Position sensor for electromagnetic actuator
JP2006242341A (en) 2005-03-04 2006-09-14 Smc Corp Actuator with position detection mechanism
JP7043555B2 (en) 2020-09-04 2022-03-29 Ckd株式会社 3D measuring device

Also Published As

Publication number Publication date
JPS5958305A (en) 1984-04-04

Similar Documents

Publication Publication Date Title
US3694088A (en) Wavefront measurement
US3829219A (en) Shearing interferometer
Koliopoulos Interferometric optical phase measurement techniques
CN106767489B (en) Small dynamic deformation measuring system and measurement method in digital speckle interference face
JP7082137B2 (en) Radius of curvature measurement by spectral control interferometry
CN114459619B (en) A real-time online measurement device and method for phase shift
CN105091781B (en) Method and device for measuring optical surface by single-frame interference fringe pattern
CN116719147A (en) An automatic focusing method, electronic device and medium for digital holographic imaging
JPH0437362B2 (en)
CN116804536B (en) Three-dimensional measurement method and device based on novel dynamic structured light
US4347000A (en) Interferometric system
CN114295203A (en) Vortex intensity measuring device and method for vortex light beam
Quan et al. Application of the holographic carrier fringe and FFT technique for deformation measurement
JP3233723B2 (en) Phase pattern difference discriminator
KR102600341B1 (en) Profilometry Device and Method
Zhang et al. Spatiotemporal phase unwrapping and its application in fringe projection fiber optic phase-shifting profilometry
JP7503972B2 (en) Shape measuring device
CN119618083B (en) Displacement measuring device and displacement measuring method
US3415587A (en) Local reference beam generation for holography
CN111780688A (en) A line structured light projection device and profile measurement method
JPH046882B2 (en)
CN119395962B (en) Microscopic imaging device and phase reconstruction method thereof
JPH11237209A (en) High sensitivity space positioning method
JPS6353481B2 (en)
KR960034979A (en) Optical phase interferometer and surface measurement method using the same