JPH0629703B2 - 変形の測定方法 - Google Patents

変形の測定方法

Info

Publication number
JPH0629703B2
JPH0629703B2 JP15429888A JP15429888A JPH0629703B2 JP H0629703 B2 JPH0629703 B2 JP H0629703B2 JP 15429888 A JP15429888 A JP 15429888A JP 15429888 A JP15429888 A JP 15429888A JP H0629703 B2 JPH0629703 B2 JP H0629703B2
Authority
JP
Japan
Prior art keywords
deformation
speckle
cross
pattern
speckle 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 - Fee Related
Application number
JP15429888A
Other languages
English (en)
Other versions
JPH01320412A (ja
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.)
Hamamatsu Photonics KK
Original Assignee
Hamamatsu Photonics KK
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 Hamamatsu Photonics KK filed Critical Hamamatsu Photonics KK
Priority to JP15429888A priority Critical patent/JPH0629703B2/ja
Priority to US07/367,659 priority patent/US4967093A/en
Priority to DE8989111389T priority patent/DE68904993T2/de
Priority to EP89111389A priority patent/EP0347912B1/en
Publication of JPH01320412A publication Critical patent/JPH01320412A/ja
Publication of JPH0629703B2 publication Critical patent/JPH0629703B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Description

【発明の詳細な説明】 【産業上の利用分野】
本発明は、変形の測定方法に係り、特に、変形前後の物
体表面の一部をレーザビームで照射してスペツクル模様
を現出し、該変形前後のスペツクル模様をそれぞれ光電
変換して得られる信号間の相互相関関数を求め、該相互
相関関数の極値の位置として求められるスペツクル模様
の移動量から物体の変形量を決定する変形の測定方法の
改良に関するものである。
【従来の技術】
レーザビームを粗面に当てた時、拡散光の干渉によつて
生ずるスペツクル模様は、表面に変位や変形が起きる
と、徐々に変形しながら移動する。そこで、スペツクル
模様を光電的に走査し、得られる信号の相関ピーク位置
からスペツクル移動を求め、このスペツクル移動と表面
の変位や変形との関係を利用して、物体の並進、回転、
歪み等による微小な変形を測定する、いわゆるスペツク
ル相関法が提案されている(特公昭59−52963、
レーザー科学研究No.6(1984)pp152〜15
4、最新精密計測技術(昭和62年7月1日)pp241
〜244)。 このスペツクル相関法で、最も実用に近いのは、第8図
に示す如く、一次元イメージセンサ15とマイクロコン
ピユータ16を用いたものであり、これによつて、1μ
m以上の平行移動と10-5rad程度の回転が測定されて
いる。 この装置においては、物体10上の測定点を、レーザ源
12で発生した、直径1m程度のレーザビーム13で、
必要に応じて拡大レンズ14を介して照射し、その拡散
反射光の中に一次元イメージセンサ15を配置してお
く。その際、ビーム径wとセンサ距離Loを調節して、
およそλLo/w(λはレーザビームの波長)で与えられる
センサ15上でのスペツクルの平均径を、そのピツチ
(10〜20μm)より大きくとつておく。又、一次元
イメージセンサ15の軸は、光学系と物体変位の種類
(平行移動、回転、歪みの方向)で決まるスペツクル移
動の方向に合致させておく。 一次元イメージセンサ15の出力をA/D変換してマイ
クロコンピユータ16に入れ、相関器18で物体の変形
前後の出力の間の相互相関関数を計算すると、そのピー
ク位置としてスペツクル移動がほぼ実時間で求められ
る。この際、相関関数の計算時間を短縮するために、一
次元イメージセンサ15の出力信号をその平均値に関し
て2値化し、いわゆる特性相関を計算することも提案さ
れている。スペツクルのコントラストが高いので、こう
して得られるピーク位置は通常の相関関数のそれと常に
一致する。従つて、相互相関関数の極値の位置からスペ
ツクル移動を検出することができる。 前記相互相関関数を求めるに際して、従来のスペツクル
相関法では、物体の変形前のスペツクル模様を基準デー
タとして固定し、物体の変形が持続中のスペツクル模様
を比較データとして相互相関関数を求めるか、又は、物
体の変形が持続中のスペツクル模様を比較データとし、
この比較データの1つ前の走査で得られたデータを基準
データとして、基準データを毎回更新して相互相関関数
を求めていた。 しかしながら、前者のように基準データを固定してしま
うと、物体の変形に伴つてスペツクル模様も変化し、得
られる相互相関関数の極値が周囲の無関係のピークより
低下してしまい、正しい極値の位置が得られなくなるた
め、測定範囲が限定されるという問題点を有していた。 一方、後者のように基準データを毎回更新する方法で
は、前回と今回の2つの走査時間の間に一次元イメージ
センサ15のピツチ間隔の半分以下の変位が生じたとき
には、基準データと比較データが同一となるため、極値
の位置は移動しない。従つて、物体の変位が無視され
る。この誤差は、毎回の走査で繰返されるため、2つの
走査時間内における変位がイメージセンサのピツチ間隔
の半分以下であるような低速の変位は検出できないとい
う問題点を有していた。
【発明が達成しようとする課題】
本発明は、前記従来の問題点を解消するべくなされたも
ので、変位測定の範囲を誤差を生じることなく拡大し
て、測定精度を向上することが可能な変形の測定方法を
提供することを目的とする。
【課題を達成するめの手段】
本発明は、変形前後の物体表面の一部をレーザビームで
照射してスペツクル模様を現出し、該変形前後のスペツ
クル模様をそれぞれ光電変換して得られる信号間の相互
相関関数を求め、該相互相関関数の極値の位置として求
められるスペツクル模様の移動量から物体の変形量を決
定する変形の測定方法において、第1図にその要旨を示
す如く、前記相互相関関数の極値が予め設定された値よ
り低下するか、又は、前記極値の位置の移動が予め設定
された範囲より逸脱した時に、相互相関の基準パターン
を更新することによつて、前記目的を達成したものであ
る。
【作用及び効果】
第2図に示す如く、物体10の測定領域Oを、レーザ源
12からのレーザビーム13で必要に応じて拡大レンズ
14を介して照射し、得られるスペツクル模様を観察面
30で観察する場合を考える。 ここで、物体面上の座標軸をx、y、z、レーザビーム
13の発散点の距離OS=Ls、発散点の方向をlsx、l
sy、lsz、物体面と観察面30の距離をLo、観察点P
の方向をlx、ly、lz、レーザビーム13で照射し
た領域における物体10の並進、回転、歪みの成分をそ
れぞれ(ax、ay、az)、(Ωx、Ωy、Ωz)、(εx
x、εyx、εyy)とする。 この条件下で、物体10が変形を受ける前後における観
察点Pでのスペツクル模様の強度分布I1(x、y)と
2(x、y)の間の相互相関関数C(、)を計算
する。 C(、)=<I1(x、y) ×I2(x+y+)>…(1) ここで、<>は集合平均を意味する。 この(1)式を計算すると、C(、)が、=A
x、=Ayで最大値をとることがわかる。ここでA
x、Ayは、次式で与えられ、物理的には物体変形によ
るスペツクル模様の移動量に相当する。 Ax=−ax[(Lo/Ls)(lsx2−1)+lx2
−1] −ay[(Lo/Ls)lsxlsy+lxly] −az[(Lo/Ls)lsxlsz+lxlz] −Lo[Ωz(lsy+ly) −Ωy(1sz+1z)+εxx(1sx+1x)+εxy(lsy+
ly)]……(2) Ay=−ax[((Lo/Ls)(lsylsx+lylx] −ay[(Lo/Ls)(lsy2−1)+ly2
−1] −az[(Lo/Ls)lsylsz+lxlz] −Lo[−Ωz(lsx+lx) −Ωx(lsz+lz)+εyy(lsy+ly) +εxy(lsx+lx)]……(3) 従つて、前記観察面30に一次元イメージセンサを配置
してスペツクルの移動量Ax、Ayを観測すれば、該一
次元イメージセンサの出力波形は、物体変位前後で第3
図(A)に示す如く変化し、その自己相関波形は、第3
図(B)に示す如くとなり、相互相関波形は第3図
(C)に示す如くとなる。 このような装置を用いて物体の変形を測定するに際し
て、相互相関関数を求める際に、本発明においては、従
来技術のように、基準データを固定してしまつたり、あ
るいは、基準データを毎回更新するのではなく、相互相
関関数の極値が設定値より低下するか、又は、極値の位
置の移動が設定範囲より逸脱した時に、基準パターンを
更新するようにしている。従つて、物体の変形に伴うス
ペツクル模様の移動が大きい場合には、基準データが更
新されるので、測定範囲が限定されることはない。逆
に、物体の変形が小さく、極値の位置の移動が設定範囲
内にある場合には、基準パターンが固定されるので、基
準データを毎回更新する場合に比べて、低速の変位も検
出できる。従つて、変位測定の範囲を誤差なく拡大し
て、測定精度を向上することができる。
【実施例】
以下、図面を参照して、本発明の実施例を詳細に説明す
る。 本実施例は、第4図に示す如く、物体10の表面にレー
ザビーム13を照射してスペツクル模様を発生させるレ
ーザ源12を備えている。このレーザ源12によつて発
生されたスペツクル模様を光電変換するための受光素子
としては、第5図に示すような、例えば短辺が13μ
m、長辺が2.5mm、配設ピツチが25μmの短冊型の
受光素子(・・・22n-1、22n、22n+1・・・)を
有する受光素子列20、24が、列をなす角度を変えて
2組配設されている。 該2組の受光素子列20、24で、それぞれ光電変換さ
れたスペツクル模様は、それぞれ本発明に係る相関器4
0、42に入力され、スペツクル模様の移動の前後の相
互相関関数の極値の位置の変化が検出される。 該相関器40及び42は、第6図に示す如く、前記受光
素子列20(相関器40)又は24(相関器42)より
得られるスペツクル模様の2値化された電気信号入力A
の1フレーム前のデータを一時的に保持するためのメモ
リ40Aと、該メモリ40Aから出力される基準パター
データBと前記受光素子列20又は24から入力される
比較データCが入力され、この基準パターンデータBと
比較データCの相関値Dを出力する相関用IC(例えば
TRW社のTDC1023)40Bと、基準パターンデ
ータBを更新するためのゲート40Cと、各シフトクロ
ツクCL毎の相関値Dのうち、極値Eとそのタイミング
出力Fを出力する最大値検出回路40Dと、最大値のタ
イミング出力FとシフトクロツクCLが入力され、極値
の位置Gを出力するカウンタ40Eと、前記最大値検出
回路40Dから出力される極値Eとその閾値を比較し、
極値Eが閾値未満である場合に基準パターン交替信号H
を出力する比較回路40Fと、前記カウンタ40Eから
出力される極値の位置Gとその最大範囲の閾値を比較
し、極値の位置が最大範囲を超えている場合に基準パタ
ーン交替信号Iを出力する比較回路40Gと、同じく極
値の位置Gとその最小範囲の閾値を比較し、極値の位置
Gが最小範囲未満である場合に基準パターン交替信号J
を出力する比較回路40Hと、前記基準パターン交替信
号H、I、Jの論理和を前記ゲート40C及び後出FI
FOメモリ40Jに出力するORゲート40Iと、前記
カウンタ40E出力の極値の位置Gと基準パターン交替
信号H+I+Jを一時的に格納して、CPUバス40K
に出力するFIFO(フアーストインフアーストアウ
ト)メモリ40Jとから構成されている。 以上の構成により、極値の値Eがある閾値より低下する
か、極値の位置Gが設定範囲より逸脱した場合に、基準
パターン更新が自動的に行われる。 前記相関器40、42で検出された極値の位置の変化は
コンピユータ44に入力され、データ処理を行つて、ス
ペツクル模様の移動量から検出すべき変形量を算出す
る。 コンピユータ44は、ステツピングモータコントローラ
46に指令を与えて、物体10が載置されているリニア
ステージ48をx軸方向に移動させ、スペツクル模様を
x軸方向に移動させると共に、タイミング回路50を介
して、前記相関器40、42に必要なタイミング信号を
入力する。 第7図に、リニアステージ48をx軸方向に6mm移動し
た場合の、極値Eに対する閾値を100%、90%、8
0%、70%、60%以上とした測定例を比較して示
す。100%の測定例は、毎回基準パターンが更新さ
れ、1つ前のフレームを基準パターンとしている従来例
に等しく、一方、60%の測定例は、閾値が低いため基
準パターンの更新が発生せず、基準パターンを固定した
従来例とほぼ同様である。 これに対して、基準パターンを本発明により設定値70
%〜90%として更新した場合には、リニアステージ4
8の実際の移動量に対して、スペツクル模様の移動量か
ら求めたリニアステージ移動量(測定値)もほぼ対応し
ていることが明らかである。 本実施例においては、受光素子列20、24の各素子を
短冊型とすると共に、受光素子列を2組設けているの
で、スペツクル模様の移動の向きと受光素子列の向きを
一致させる必要がなく、又、各方向の成分をそれぞれ測
定することができる。なお、受光素子の形状や受光素子
列の組数は前記実施例に限定されず、例えば所定のx軸
方向の移動成分のみを抽出すればよい時には、受光素子
列を1組のみとしてもよい。又、受光素子列を3組以上
設けて、更に高精度の測定を行うことも可能である。
【図面の簡単な説明】
第1図は、本発明に係る変形の測定方法の要旨を示す流
れ図、第2図は、スペツクル測定法の測定原理を示す斜
視図、第3図(A)、(B)、(C)は、同じく一次元
イメージセンサの出力波形、自己相関波形、相互相関波
形をそれぞれ示す線図、第4図は、本発明に係る変形の
測定方法を実施するための測定装置の構成を示す、一部
ブロツク線図を含む斜視図、第5図は、前記実施例にお
ける受光素子列を示す正面図、第6図は、同じく相関器
の構成を示すブロツク線図、第7図は、前記実施例にお
ける極値の閾値を変えた場合のリニアステージ移動量と
測定値の関係を示す線図、第8図は、従来のスペツクル
相関法の測定原理を示す斜視図である。 10……物体、 12……レーザ源、 13……レーザビーム、 20、24……受光素子列、 40、42……相関器、 40A……メモリ、 40B……相関用IC、 40C……ゲート、 40D……最大値検出回路、 40E……カウンタ、 40F、40G、40H……比較回路、 B……基準パターンデータ、 C……比較データ、 D……相関値、 44……コンピユータ、 48……リニアステージ、 50……タイミング回路。

Claims (1)

    【特許請求の範囲】
  1. 【請求項1】変形前後の物体表面の一部をレーザビーム
    で照射してスペツクル模様を現出し、該変形前後のスペ
    ツクル模様をそれぞれ光電変換して得られる信号間の相
    互相関関数を求め、該相互相関関数の極値の位置として
    求められるスペツクル模様の移動量から物体の変形量を
    決定する変形の測定方法において、 前記相互相関関数の極値が予め設定された値より低下す
    るか、又は、前記極値の位置の移動が予め設定された範
    囲より逸脱した時に、相互相関の基準パターンを更新す
    ることを特徴とする変形の測定方法。
JP15429888A 1988-06-22 1988-06-22 変形の測定方法 Expired - Fee Related JPH0629703B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP15429888A JPH0629703B2 (ja) 1988-06-22 1988-06-22 変形の測定方法
US07/367,659 US4967093A (en) 1988-06-22 1989-06-19 Deformation measuring method and device using cross-correlation function between speckle patterns with reference data renewal
DE8989111389T DE68904993T2 (de) 1988-06-22 1989-06-22 Deformationsmessverfahren und vorrichtung, die die kreuzkorrelationsfunktion zwischen speckle-bildern benuetzt.
EP89111389A EP0347912B1 (en) 1988-06-22 1989-06-22 Deformation measuring method and device using cross-correlation function between speckle patterns

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP15429888A JPH0629703B2 (ja) 1988-06-22 1988-06-22 変形の測定方法

Publications (2)

Publication Number Publication Date
JPH01320412A JPH01320412A (ja) 1989-12-26
JPH0629703B2 true JPH0629703B2 (ja) 1994-04-20

Family

ID=15581075

Family Applications (1)

Application Number Title Priority Date Filing Date
JP15429888A Expired - Fee Related JPH0629703B2 (ja) 1988-06-22 1988-06-22 変形の測定方法

Country Status (1)

Country Link
JP (1) JPH0629703B2 (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03274402A (ja) * 1990-03-26 1991-12-05 Rikagaku Kenkyusho 2軸レーザースペックル歪計
US7295324B2 (en) * 2004-07-13 2007-11-13 Mitutoyo Corporation System and method for improving accuracy in a speckle-based image correlation displacement sensor

Also Published As

Publication number Publication date
JPH01320412A (ja) 1989-12-26

Similar Documents

Publication Publication Date Title
US4967093A (en) Deformation measuring method and device using cross-correlation function between speckle patterns with reference data renewal
EP0134597B2 (en) Measuring system based on the triangulation principle for the dimensional inspection of an object
US4862894A (en) Apparatus for monitoring bloodstream
JPH0428005Y2 (ja)
US5329458A (en) Device and method for detecting position of edge of cutting tool
JPS61116611A (ja) 距離測定方法
JP2771594B2 (ja) 物体の変位測定方法及び装置
US5061860A (en) Deformation measuring method and device using comb-type photosensitive element array
JPH0629703B2 (ja) 変形の測定方法
JPS6244203B2 (ja)
JPH0914935A (ja) 3次元物体の測定装置
JPH07253304A (ja) 多軸位置決めユニットおよびこれにおける測長方法
JP2793931B2 (ja) 光の二次元測定方法および表面粗さ測定方法と測定装置
JPH0635921B2 (ja) スペックル測長計
JPH0629704B2 (ja) 変形の測定方法
EP0157431A1 (en) Procedure to measure the dimensions of a body in movement in a three-dimensional field, and an optoelectronic device to carry out such procedure
JP3192461B2 (ja) 光学的測定装置
JP2626611B2 (ja) 物体形状測定方法
JP2946381B2 (ja) 表面粗さ測定方法及び装置
US4422737A (en) Device for obtaining topographic picture of surface of rotating object
KR0140663B1 (ko) 인쇄회로기판의 솔더 도포형상 감지장치
SU731278A1 (ru) Устройство дл измерени малых размеров
JPS61151407A (ja) 非接触径測定装置
JPH0471161B2 (ja)
CN117249766A (zh) 六自由度位移测量方法及装置

Legal Events

Date Code Title Description
S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313532

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

LAPS Cancellation because of no payment of annual fees