JP3955899B2 - 電子的スペックル干渉法を用いた変形計測方法および装置 - Google Patents
電子的スペックル干渉法を用いた変形計測方法および装置 Download PDFInfo
- Publication number
- JP3955899B2 JP3955899B2 JP2002275775A JP2002275775A JP3955899B2 JP 3955899 B2 JP3955899 B2 JP 3955899B2 JP 2002275775 A JP2002275775 A JP 2002275775A JP 2002275775 A JP2002275775 A JP 2002275775A JP 3955899 B2 JP3955899 B2 JP 3955899B2
- Authority
- JP
- Japan
- Prior art keywords
- phase
- intensity signal
- component
- speckle
- deformation
- 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 - Lifetime
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02015—Interferometers characterised by the beam path configuration
- G01B9/02022—Interferometers characterised by the beam path configuration contacting one object by grazing incidence
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
- G01B11/161—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge by interferometric means
- G01B11/162—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge by interferometric means by speckle- or shearing interferometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02001—Interferometers characterised by controlling or generating intrinsic radiation properties
- G01B9/02007—Two or more frequencies or sources used for interferometric measurement
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02083—Interferometers characterised by particular signal processing and presentation
- G01B9/02084—Processing in the Fourier or frequency domain when not imaged in the frequency domain
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02094—Speckle interferometers, i.e. for detecting changes in speckle pattern
- G01B9/02095—Speckle interferometers, i.e. for detecting changes in speckle pattern detecting deformation from original shape
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Mathematical Physics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Description
【発明の属する技術分野】
本発明は、電子的スペックル干渉法を用いた変形計測方法および装置に関し、詳しくは、動的物体の時間的変形についてのスペックル干渉画像を解析する際に、位相アンラッピングを良好に行い得る電子的スペックル干渉法を用いた変形計測方法および装置に関するものである。
【0002】
【従来の技術】
干渉法によって鏡面状態にある物体の表面状態を測定するために、位相分布(表面状態に対応)を簡単に求めたいという要求は、近年の技術の高度化に伴い光学分野や電子分野を中心に極めて強いものとなっている。特に、干渉測定において、このような観察物体の位相分布を求める手法については、主として位相シフト法と空間的フーリエ変換法(空間キャリアを導入したフーリエ変換法;以下同様)を用いたものが従来より知られている。
【0003】
位相シフト法は、干渉計の物体光と参照光の間で、一般には2πを整数分の一に分割した位相角ずつ位相をずらした各干渉パターン画像のパターン情報に基づき観察物体の位相分布を求めるもので、2πをN等分したN枚の干渉パターン画像から得られる所定の空間領域における強度信号は下記(1)式で与えられる。
【0004】
【数2】
【0005】
ここで、φ(x,y)は求めるべき位相である。また、I0は平均強度、γは干渉パターンの可視度(モデュレーション)であり、いずれも未知量である。簡単な例として、N=4の場合は、下記(2)式より求められる。
【0006】
【数3】
【0007】
次に、フーリエ変換法について説明する。干渉計の一方の光路をy軸周りに微小角θだけ傾けると、下記(3)式で与えられる空間的な強度分布が得られる。
【0008】
【数4】
【0009】
ここで、f0=sinθ/λは、光路を傾けたことによって生じたキャリア縞の空間周波数である。上記(3)式をx方向にフーリエ変換すると、下記(4)式が求まる。
【0010】
【数5】
【0011】
この(4)式の右辺は、f0が十分大きければ空間周波数軸上で3つの項を互いに分離できることを示している。そこで、f0を十分大きくとって右辺の第2項のみを取り出し、他の項をカットするフィルタに通すと、下記(5)式が得られる。
【0012】
【数6】
【0013】
この(5)式の右辺の実部と虚部の比(逆正接)から、下記(6)式が求められる。
【0014】
【数7】
上述した従来の2方法は、いずれも、求めようとする未知量φ(x,y)が他の未知量I0、γとは無関係に求められる。位相値は[−π、π]の間の主値として求められる。位相シフト法を用いた際に、位相が線形的に変化する様子を図8に示す。上記(6)式の分母が余弦関数、分子が正弦関数であり、その比の逆正接により、2πごとに位相跳びが生じる鋸歯状の位相分布が得られる。位相跳びの位置を判定し、この各位値の右側のデータに2πの位相値を加える(または差し引く)ことによって位相跳びを補正し、物体の形状に比例した位相分布を得ることができる。この処理を位相アンラッピングと称する。
【0015】
一般に、位相シフト法では2次元干渉パターンの1点の位相を、その点の位相シフトした複数のデータから解析することができる。したがって、空間的に別の点の影響を受けることはない。それに対して、フーリエ変換法においては、空間キャリアからなる1直線上のデータ全てに亘ってフーリエ変換の演算を施さなければならないので、各点の位相を他の点と独立に求めることはできない。
【0016】
その一方で、位相シフト法では複数枚の位相シフトした干渉パターン画像を取り込む期間に亘って観察物体が静止していなければならない。それに対して、フーリエ変換法においては、キャリア成分さえ作っておけば1枚の干渉パターン画像から位相解析ができるので、動的な現象に適しているといえる。換言すれば、位相シフト法は観察対象に対して空間的な制限は緩いが時間的には定常であることが必須であり、フーリエ変換法においては観察対象は空間的には空間キャリア周期に比べて十分緩やかな分布あるいは一様とする必要があるが、時間的には動いていてもかまわないといえる。
【0017】
しかしながら、観察対象が時間的に変化し、かつ空間的にも一様でないあるいは急激に変化する場合にはこれらの方法を適用することはできない。例えば、材料の塑性変形から破壊に至る過程は非線形的で、その歪測定においては、時間的かつ空間的な変形分布を測定する必要があるが、このような測定において上記2つの方法を適用することは理論的に難しい。
【0018】
このような時間的、かつ空間的な変動がある場合に有効な干渉法としては、動的スペックル干渉法が知られている。
【0019】
スペックル干渉法は、粗面物体をレーザ光で照射したときに観察面に生じる斑点状の模様(スペックルパターン)を利用する干渉法である。スペックルパターンは一般の結像系では画像ノイズとして好ましくないものとされる。しかし、位相情報を担っておりその変化から変形を見積もることができる。また、スペックル干渉法では光の波長を基準とした高精度な変形測定が可能になる。
【0020】
図9に、2光束照射型のスペックル干渉装置を示す。粗面物体である観察物体100は、xz面内で略対称に配された、レーザ光源101からの2光束102A、102Bによって照明される。観察物体100によって拡散反射した光束はCCDカメラ103の結像面上に干渉スペックルパターンを形成する。この後、得られた干渉スペックルパターン画像を解析して、観察物体100の表面形状に応じた位相解析を行うこととなる。
この位相解析手法として、和差法なる手法が知られている。
【0021】
一般に、スペックルパターン画像の位相解析においては、観察物体100の変形前後のスペックルパターンを撮像し、それらの画像点毎の強度の差を計算する。相関の強い場所、すなわち変形による位相変化が0または2πの整数倍の位置では差強度が0に近くなり、一方、相関の弱い場所では差強度が大きな値をとることから、2画像の強度差Isubの絶対値を計算することによって、変形量に応じた相関パターンが得られる。
【0022】
上述した和差法においては、強度差Isubの情報に加えて強度和Iaddの情報を利用することによって、可視度(モデュレーション)γとは無関係に位相を求める。
【0023】
ここで、図9に示すように、2つの光束102A、102Bの各光路上にそれぞれシャッタ104A、104Bを置き、一方の光束のみで物体を照射したときの時間領域における画像点毎の各強度分布I1(x;t)、I2(x;t)をあらかじめ測定しておく。
スペックル干渉装置の光学系(面内変形あるいは面外変形測定系)において得られる干渉パターンI(x;t)は一般的に下記(7)式のように表される。
【0024】
【数8】
【0025】
ここで、I0(x;t)はI1(x;t)およびI2(x;t)の平均強度、θ(x;t)はランダムなスペックル位相、γ(x;t)は変調度(モデュレーション)、φ(x;t)は物***相である。
【0026】
次に、物体変形前の強度をIbefore(t=t1)とし、変形後の強度をIafter(t=t2)とするとそれぞれ下記(8)式、(9)式のように表される。
【0027】
【数9】
【0028】
これらの強度パターンに対して差Isubおよび和Iaddを計算し同時に局所平均を取ると、下記(10)式、(11)式のように表される。
【0029】
【数10】
【0030】
ここで、< >は局所的空間平均を表す。c,c′は定数である。平均強度I0は局所的な時間平均により求めることができる。また、Δφ=φ2−φ1であり時刻t1と時刻t2の間の物体の変形を表す。
これら2つの成分に対して下記(12)式で表される演算を行うことにより、物***相が求められることになる。
【0031】
【数11】
【0032】
しかしながら、上述した和差法を用いた場合には、上述した(10)式および(11)式に示すように正弦および余弦の各成分に対して絶対値をとっているため、図10に示すような直線的に増加する物***相に対して、実際に得られる位相曲線は、前述したような位相跳びは生じないが、その代わりに図10の実線で示されるような位相の折り返しを有する。これに対しては、折り返しの各区画に番号nを付け、下記(13)式を用いて位相を接続する(位相アンラッピングする)ことで元の物***相を復元することができる。
【0033】
【数12】
【0034】
【発明が解決しようとする課題】
ところで、一般的な干渉法では、物***相を求めるための、上記(12)式に対応する演算(例えば上記(6)式)に用いる正弦成分および余弦成分は絶対値を用いないため、得られる位相の折りたたみは図11(あるいは図8)に示すように鋸歯状の位相跳びが生じる。この場合は2πの明確な不連続点が生じるため物***相復元のための位相接続処理は比較的容易に実行できる。すなわち不連続点を自動的に(ソフトウエアにより)検出しその区間に2πの整数倍のオフセットを加算あるいは減算すればよい。しかしながら、上述した和差法では図10に示すように特異な位相の折りたたみが発生する。よって、不連続点が生じないため、折りたたみの区間をソフトウエアにより自動処理することは容易ではない。一般に、スペックル干渉法では図12に示すようにスペックルノイズが顕著である。特に、スペックル相関パターン(縞)が密になる領域、すなわち変形量が大きい領域では、和差法での位相接続処理が特に困難となり多くの場合自動化処理が不可能となる。
【0035】
また、和差法では上記(10)、(11)式に見られるように、局所的な空間平均を取ることが原理的に不可欠であり、この操作によって空間分解能の低下が生じる。その結果として、位相の折り返し点を検出することがさらに困難となり、位相接続処理は極めて困難なものとなる。
【0036】
本願発明者等は、このような和差法の位相接続処理の困難性に鑑み、その改良手法を提案しており(特開 2001-311613 号公報)、それによって位相の折り返し点の検出精度が向上しているが、位相の折り返し点を自動検出し得るまでには至っていない。
【0037】
なお、空間キャリアを導入する空間的フーリエ変換法とは異なり、時間キャリアを導入する時間的フーリエ変換法が近年注目されているが、相対的にノイズ成分が大きい信号を用いる上述したスペックル干渉法に対してこの時間的フーリエ変換法を適用した場合には、上記信号について最適なバンドパスフィルタを自動的に決定して処理を行なうことが困難となる。
【0038】
本発明は、上記事情に鑑みなされたもので、電子的スペックル干渉法においても、位相アンラッピング処理を行う際に、位相分布曲線における位相の接続点の検出を自動化する程度まで容易なものとすることができ、観察物体の動的な変形、振動、歪等を高精度に解析し得る電子的スペックル干渉法を用いた変形計測方法および装置を提供することを目的とするものである。
【0039】
【課題を解決するための手段】
本発明の電子的スペックル干渉法を用いた変形計測方法は、スペックル干渉法を用いて得られた、動的被観察体の位相情報を担持したスペックルパターン画像に基づき、所定位相範囲に位相ラッピングされた被観察体の位相変化曲線を解析により求め、この後、該位相変化曲線を位相アンラッピングする方法において、
所定時間毎に得られた複数の前記スペックルパターン画像に基づき、画像各点毎の時間領域における強度信号を求め、該強度信号の余弦成分を抽出し、抽出された該余弦成分に、時間領域におけるヒルベルト変換処理を施して前記強度信号の正弦成分を求め、求められた該正弦成分と前記余弦成分の比に基づいて前記画像各点毎の位相変化を求めて、前記被観察体の位相変化曲線を求めることを特徴とするものである。
この場合において、前記ヒルベルト変換処理は下式を用いて行われる。
【0040】
【数13】
【0041】
また、前記画像各点毎の時間領域における強度信号I(x;t)の位相項に、該位相項が単調増加または単調減少となるような時間キャリアを導入する位相成分ωtを加算または減算することが好ましい。
また、本発明の電子的スペックル干渉法を用いた変形計測方法は、前記動的被観察体が引っ張り試験に供された試験片であるような場合に特に有用である。
【0042】
また、本発明の電子的スペックル干渉法を用いた変形計測装置は、前記スペックル干渉法を用いて得られた、動的被観察体の位相情報を担持したスペックルパターン画像に基づき、所定位相範囲に位相ラッピングされた被観察体の位相変化曲線を解析により求めるとともに、該位相変化曲線に位相アンラッピング処理を施す装置であって、
所定時間毎に得られた複数の前記スペックルパターン画像に基づき、画像各点毎の時間領域における強度信号I(x;t)を求める強度信号演算手段と、
前記強度信号演算手段から出力された強度信号I(x;t)から、所定の平均強度信号I0(x;t)を差し引いて該強度信号の余弦成分Ic(x;t)を算出する平均成分除去手段と、
前記平均成分除去手段から出力された該強度信号の余弦成分に対し、時間領域におけるヒルベルト変換処理を施して、前記強度信号の正弦成分を算出するヒルベルト変換演算手段と、
前記平均成分除去手段において算出された該強度信号の余弦成分と、前記ヒルベルト変換演算手段において算出された該強度信号の正弦成分との比の逆正接を求める演算を行い、前記被観察体の物***相を求める物***相決定手段と、
を備えたことを特徴とするものである。
【0043】
また、前記強度信号I(x;t)に対して適正な時間キャリア成分を重畳せしめる時間キャリア重畳手段は、例えば、スペックル干渉を生ぜしめるための2系の照明光束のうちの一方に、これら2系の照明光束の間に所定の光路長差を発生させる光路長差発生手段とする。
【0044】
【発明の実施の形態】
以下、本発明に係る電子的スペックル干渉法を用いた変形計測方法の実施形態について図面を用いて説明する。
【0045】
<基本構成>
図1および図2は、それぞれ観察物体の面内変位および面外変位をスペックル干渉測定する光学系の一例を示すものである。
【0046】
図1に示す面外変位計測用の干渉光学系装置では、半導体レーザ光源(LD)11より出射したレーザ光が観察物体12の前方に配置された拡散板(Diffuser)13により一部が散乱されてランダムなスペックルパターンを形成する。一方、拡散板13を透過したレーザ光は観察物体12を照明し、その表面で散乱されて、上記と同様にスペックルパターンを形成する。これら2つの散乱光はCCDカメラ14の撮像面上に、ランダムな干渉スペックルパターンの像を形成し、この干渉パターンが光電的に読み取られる。また、図2に示す面内変位計測用の干渉光学系装置では、観察物体12が、ハーフミラー17により分割された、物体法線に対して略対称な2つのレーザ光により照明されて各々独立なスペックルパターンが形成され、それら2つのスペックルパターンの干渉によりCCDカメラ14の撮像面上に干渉スペックルパターンの像が形成される。この後、上記と同様に、この干渉スペックルパターンが光電的に読み取られる。なお、PZT15については後述する。
【0047】
上述したいずれの装置においてもCCDカメラ14により検出されるランダムな干渉パターンの強度I(x;t)
は下記(14)式のように表される。
【0048】
【数14】
【0049】
ここで、xは座標、tは時刻、I0(x;t)は平均強度、γ(x;t)は変調度、θ(x;t)はランダムなスペックル位相、φ(x;t)は観察物体の変位に伴う物***相である。
【0050】
なお、本実施形態では、上記物***相φ(x;t)を高精度で求めることが目標となる。
そこで本実施形態においては、まず、上記(14)式の強度に対して局所的な時間平均により平均強度I0(x;t)を求め、干渉パターンの強度I(x;t)から差し引くことにより干渉パターン強度の余弦成分のみを取り出す。これをIc(x;t)とおけば下記(15)式のように表すことができる。
【0051】
【数15】
次に、上記(15)式で表される時間信号に対してHilbert変換を適用する。
関数f(t)に対するHilbert変換は下記(16)式で定義される。
【0052】
【数16】
【0053】
上記(15)式の余弦成分に対してHilbert変換を施すことにより、下記(17)式のような正弦成分が得られる。
【0054】
【数17】
【0055】
上述した余弦成分Ic(x;t)および正弦成分Is(x;t)の比を求め、さらにこの比の逆正接を求めることにより下記(18)式のように物***相が決定される。
【0056】
【数18】
【0057】
通常、ランダムなスペックル位相θ(x;t) の変化は物***相の変化に比べて遅いので下記(19)式のように近似することができる。
【0058】
【数19】
【0059】
よって求める物***相は下記(20)式のように決定される。
【0060】
【数20】
【0061】
また、観測時間to(0<t<to)が長く(19)式が成り立たない場合は観測時間toをスペックル位相が一定と見なせるn個の区間td(to=ntd)に分割し、(20)式においてθ(x;0)をθ(x;ktd),k=0,1,…n−1と置き換えることにより各区間において物***相を求める。すなわち適宜参照位相の更新を行なう。最終的に各区間を接続することにより長時間に亘る動的な現象を解析することが可能となる。
このようにして求められた、時間領域における位相曲線は、前述したように、位相値が−πおよびπを主値として、ラッピングされた(折りたたまれた)状態となっている。そこで、このようなラッピングされた(折りたたまれた)位相曲線に対して、これを本来の曲線に戻す位相アンラッピング処理を施す。
【0062】
以上の演算処理をスペックル干渉画像上の画像各点毎に行うことにより、観察物体の変形を動的に、かつ高精度に定量化して解析することができる。
【0063】
<位相項の符号の確定>
ところで、上述した原理の説明においては、式(15)から明らかなように位相項の符号に不確定性が生じる。この不確定性を取り除くためには、式(15)の位相項に意図的に既知の位相成分ωtを導入する。すなわち、下記(21)式において位相項であるθ(x;t)+φ(x;t)+ωtが単調増加あるいは単調減少となるように決定された位相成分ωtを、いわゆる時間キャリア成分として導入する。
【0064】
【数21】
【0065】
具体的には、例えば図1に示すような、半導体レーザ11の注入電流を変調する方法、あるいは図2に示すような、ミラー16を圧電素子(PZT)15によって移動させることにより光路長を変化させる方法等の種々の周知の手法により行うことができる。
【0066】
<装置構成>
図3は、本発明の実施形態に係る電子的スペックル干渉法を用いた変形計測装置を示す概略構成図である。なお、この装置は観察物体の面内変位計測を行う装置の一例である。
【0067】
図2においても説明したように、観察物体12は、ハーフミラー17により分割された、物体法線に対して略対称な2つのレーザ光により照明され、それら2つのレーザ光の観察物体12上からの散乱により各々独立にランダムなスペックルパターンが形成される。これら2つのスペックルパターンの干渉によりCCDカメラ14の撮像面上に干渉パターンが形成される。ここで圧電素子(PZT)15は前述した時間キャリア成分を導入する目的で、干渉計の片方の系の光路長を変化させるために用いる。圧電素子15、CCDカメラ14および半導体レーザ(LD)11はタイミング制御装置(コントローラ)21により制御され、CCDカメラ14により取り込まれた干渉パターンの画像データはコンピュータ22内のメモリに格納される。コンピュータ22においては、この画像データに対し、ヒルベルト変換法を用いた位相解析処理を施し、観察物体の変形歪量等の各情報を得る。得られた観察物体に関する各情報は表示装置(ディスプレイ)23に表示される。
【0068】
なお、スペックル画像は、図4に示すように、所定時間毎に取り込まれるようになっており、これにより、この画像内の画像各点毎に強度I(x,y;t)の時間変化(図4右端のグラフ参照)が得られることになる。
ここで、上記画像内の1点(x,y)の強度I(x,y;t)は下記(22)式により表される。
【0069】
【数22】
【0070】
ここで、ωtは時間キャリア成分であるがこれを導入しない場合はω=0とする。
【0071】
また、上記コンピュータ22における各演算処理は、図5のブロック図に示されるように、入力信号としての上記強度I(x,y;t)(以下、単にI(x;t)と記す)に対し、平均成分除去手段32、ヒルベルト変換演算手段33、物***相導出手段34および位相接続手段35において順次なされるようになっている。
なお、上記強度I(x;t)は、前述したPZT等の時間キャリア重畳手段により適切な時間キャリア成分が重畳されたものとされている。
【0072】
入力された強度I(x;t)は、まず平均成分除去手段32において上記(15)式を用いて平均強度I0(x;t)を除去する演算処理がなされる。
また、ヒルベルト変換演算手段33においては、上記(17)式を用いて、ヒルベルト変換による正弦成分の算出が行われる。
【0073】
さらに、物***相導出手段34においては、上記(18)式を用いて余弦成分と正弦成分の比の逆正接を求める演算がなされ、次に、上記(20)式を用いて物***相が決定される。
【0074】
この後、位相接続手段35において、物***相曲線のアンラッピング演算処理が行われて、観察物体の動的変形情報が求められる。すなわち、前述したように、位相曲線の位相跳びの位置を判定し、その判定位置毎に、その位置の右側のデータに対して、さらに2πの位相値を加える(または差し引く)ことによって位相跳びを補正し、物体の変形に比例した位相分布を得る。
【0075】
最後に、3次元的変形分布形成手段36により、上記で得られた、画像各点についての、位相アンラッピング処理がなされて得られた位相曲線から、図6に示す如き3次元的変形分布を得、このデータをディスプレイ23上で3次元的に表示し得るような形態に変換して出力する。
なお、上記各手段32〜36はコンピュータ22内に格納されたプログラムにより構成されている。
【0076】
なお、本発明方法は、例えば材料の引っ張り試験(あるいは疲労試験)を行う際等の変形計測に適用すると特に有効であるから、以下、アルミニウム合金の引張り試験において刻々変化するその変形の解析処理に本発明方法を適用した場合の効果について、具体的な実施例を用いて説明する。
【0077】
<実施例>
一般に、アルミニウム合金等の材料を引っ張った場合、弾性変形領域では一様な伸びが生じる。これをスペックル干渉法(ESPI)で観察すると、図7(A)に示すような一様な縞が認識される。位置によって縞の方向や密度が異なるのは、変形以外に予期できない試験片の設置誤差等による回転が生じていることによるものである。変形が降伏点を越え、塑性変形領域になると、変形はすべり帯と称される、一般には引張り方向に略45°傾いた狭い領域に局在する。これをスペックル干渉法(ESPI)で観察すると、図7(B)に示すように、局在部に縞が密集して現れる。
【0078】
さらに、変形速度に対して相対的に長い差分時間をもって観察すると、内部の縞構造が消失したような帯状の部分が現れ、時間的に略一定の速度で上下に移動していく様子が観察される。この部分について差分時間を十分に短くとると、縞構造が現れてくる。この部分を拡大して詳細に観察すると、図7(C)(最左図から最右図まで3.6秒毎に観察)に示すように、局在する変形が複雑に非線型的に伸展している様子がわかる。しかし、このような従来のスペックル干渉法(ESPI)では、変形の時間的変化を定性的に把握することは可能であるが、変形の大きさを定量的に求めることはできない。
【0079】
これに対し、その変形の解析処理に本発明方法を適用すれば、S/N比が比較的小さい場合等においても、上記変形の大きさを定量的に、高精度かつ容易に把握することができる。
【0080】
以下、上記実施形態装置を用い、下記実験条件で、観察物体の面内変形による位相分布を解析した。その結果を図6(A)に示す。
なお、比較例として、同様の実験条件で、前述した和差法を用いて得られた結果を図6(B)に示す。
【0081】
(1) 実験条件
物体:アルミ合金(厚み 5mm)
引っ張り速度:0.5μm/秒
観察領域:80mm×30mm
フレームレート:30フレーム/秒
照明光:半導体レーザ光
照明光の波長:532nm
【0082】
(2) 実験結果
図6(A)と図6(B)の比較から明らかなように、和差法を用いた場合には、位相接続が正しく行われていないために極めて不正確な解析結果となっているのに対し、本発明方法を用いた場合には、変形の大きさを正確に得ることができる。
【0083】
<他の手法を用いた場合との比較>
上述したように、電子的スペックル干渉法を用いて、時間領域における観察物体の変形計測を行う場合に、本発明方法に替えて他の手法を用いると種々の問題が生じる。
【0084】
すなわち、位相シフト法や空間的フーリエ変換法を用いた場合には、いずれも時間的かつ空間的に変化するスペックルパターンを扱うこと自体が困難である。また、時間的フーリエ変換法においては、S/N比が小さいスペックル干渉法における信号に対して、自動的に最適なバンドパスフィルタを決定することが困難である。さらに、和差法を用いた場合には位相接続点の特定が難しく、多くの場合自動化処理が不可能であり、データ取得条件が良好な場合を除いて、上記比較例に示すように、極めて不正確な解析結果しか得られない。
【0085】
これに対して、本発明方法を用いた場合には、上述したような各従来手法のような問題が生じる虞がなく、良好な解析結果を得ることができる。
【0086】
<態様の変更>
なお、本発明の実施形態方法においては種々の態様の変更が可能であり、例えば、位相ラッピングされた位相データに対してローパスフィルタによる平滑化処理や、スケルトニング処理(細線化処理)を施すようにしてもよい。
【0087】
また、本発明の適用対象は動的物体全般に及ぶものであり、例えば、上述した材料の引っ張り試験等の項目を含む、材料の疲労試験全般に適用可能であることは勿論のこと、動植物の各構成部分、特に人体組織の微小な動的変化にも適用可能である
【発明の効果】
以上、詳細に説明したように、本発明に係る電子的スペックル干渉法を用いた変形計測方法および装置によれば、被観察体の干渉スペックルパターン強度の余弦成分を抽出し、抽出された該余弦成分に、時間領域におけるヒルベルト変換処理を施して前記強度の正弦成分を求め、求められた正弦成分と余弦成分に基づいて前記画像各点毎の位相を求めて、前記被観察体の位相分布曲線を求めるようにしており、時間的かつ空間的に変化する被観察体であっても、容易に位相分布曲線を求めることができる上、スペックル干渉法に典型的なノイズを多く含む信号からも位相接続点の特定が容易であるから、自動化処理も可能で、極めて高精度な解析結果を得ることができる。
【図面の簡単な説明】
【図1】本発明方法の概念を説明するための面外変位スペックル干渉測定光学系の一例を示す概略図
【図2】本発明方法の概念を説明するための面内変位スペックル干渉測定光学系の一例を示す概略図
【図3】本発明の実施形態装置の構成を示す概略図
【図4】本発明の実施形態装置に取り込まれる干渉パターン画像を概念的に示す図
【図5】本発明の実施形態装置のコンピュータに含まれる各手段を示すブロック図
【図6】本実施例および比較例による、観察物体の面内変形による位相分布の解析結果を示す図
【図7】アルミニウム合金等の材料を引張った場合における、従来のスペックル干渉法(ESPI)で観察した場合の変形分布を示す図
【図8】位相シフト法における位相ラッピング曲線を示すグラフ
【図9】2光束照射型のスペックル干渉装置(測定光学系)を示す概略図
【図10】和差法における位相の折りたたみ状態を示す図
【図11】一般的な干渉法における位相の折りたたみ状態を示す図
【図12】スペックル干渉法により得られた画像におけるスペックルノイズを示す図
【符号の説明】
11 半導体レーザ光源(LD)
12,100 観察物体
13 拡散板(Diffuser)
14、103 CCDカメラ
15 圧電素子(PZT)
16 ミラー
17 ハーフミラー
32 平均成分除去手段
33 ヒルベルト変換演算手段
34 物***相導出手段
35 位相接続手段
36 3次元的変形分布形成手段
101 レーザ光源
102A、102B 光束
104A、104B シャッタ
Claims (6)
- スペックル干渉法を用いて得られた、動的被観察体の位相情報を担持したスペックルパターン画像に基づき、所定位相範囲に位相ラッピングされた被観察体の位相変化曲線を解析により求め、この後、該位相変化曲線を位相アンラッピングする、電子的スペックル干渉法を用いた変形計測方法において、
所定時間毎に得られた複数の前記スペックルパターン画像に基づき、画像各点毎の時間領域における強度信号を求め、該強度信号の余弦成分を抽出し、抽出された該余弦成分に、時間領域におけるヒルベルト変換処理を施して前記強度信号の正弦成分を求め、求められた該正弦成分と前記余弦成分の比に基づいて前記画像各点毎の位相変化を求めて、前記被観察体の位相変化曲線を得ることを特徴とする電子的スペックル干渉法を用いた変形計測方法。 - 前記画像各点毎の時間領域における強度信号I(x;t)の位相項に、該位相項が単調増加または単調減少となるような時間キャリアを導入する位相成分ωtを加算または減算することを特徴とする請求項1または2記載の電子的スペックル干渉法を用いた変形計測方法。
- 前記動的被観察体が引っ張り試験に供された試験片であることを特徴とする請求項1から3のうちいずれか1項記載の電子的スペックル干渉法を用いた変形計測方法。
- スペックル干渉法を用いて得られた、動的被観察体の位相情報を担持したスペックルパターン画像に基づき、所定位相範囲に位相ラッピングされた被観察体の位相変化曲線を解析により求めるとともに、該位相変化曲線に位相アンラッピング処理を施す、電子的スペックル干渉法を用いた変形計測装置において、
所定時間毎に得られた複数の前記スペックルパターン画像に基づき、画像各点毎の時間領域における強度信号I(x;t)を求める強度信号演算手段と、
前記強度信号演算手段から出力された強度信号I(x;t)から、所定の平均強度信号I0(x;t)を差し引いて該強度信号の余弦成分Ic(x;t)を算出する平均成分除去手段と、
前記平均成分除去手段から出力された該強度信号の余弦成分に対し、時間領域におけるヒルベルト変換処理を施して、前記強度信号の正弦成分を算出するヒルベルト変換演算手段と、
前記平均成分除去手段において算出された該強度信号の余弦成分と、前記ヒルベルト変換演算手段において算出された該強度信号の正弦成分との比の逆正接を求める演算を行い、前記被観察体の物***相を求める物***相決定手段と、
を備えたことを特徴とする電子的スペックル干渉法を用いた変形計測装置。 - 前記強度信号I(x;t)に対して適切な時間キャリア成分を重畳せしめる時間キャリア重畳手段を備え、該時間キャリア重畳手段は、スペックル干渉を生ぜしめるための2系の照明光束のうちの一方に、これら2系の照明光束の間に所定の光路長差を発生させる光路長差発生手段を設けたことを特徴とする請求項5記載の電子的スペックル干渉法を用いた変形計測装置。
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002275775A JP3955899B2 (ja) | 2002-09-20 | 2002-09-20 | 電子的スペックル干渉法を用いた変形計測方法および装置 |
US10/334,784 US6943870B2 (en) | 2002-09-20 | 2003-01-02 | Deformation measuring method and apparatus using electronic speckle pattern interferometry |
EP03002029A EP1400779A1 (en) | 2002-09-20 | 2003-01-28 | Deformation measuring method and apparatus using electronic speckle pattern interferometry |
US10/358,285 US20040059526A1 (en) | 2002-09-20 | 2003-02-05 | Deformation measuring method and apparatus using electronic speckle pattern interferometry |
KR1020030008447A KR100947015B1 (ko) | 2002-09-20 | 2003-02-11 | 전자적 스페클 간섭법을 이용한 변형계측방법 및 장치 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002275775A JP3955899B2 (ja) | 2002-09-20 | 2002-09-20 | 電子的スペックル干渉法を用いた変形計測方法および装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2004109075A JP2004109075A (ja) | 2004-04-08 |
JP3955899B2 true JP3955899B2 (ja) | 2007-08-08 |
Family
ID=31944618
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002275775A Expired - Lifetime JP3955899B2 (ja) | 2002-09-20 | 2002-09-20 | 電子的スペックル干渉法を用いた変形計測方法および装置 |
Country Status (4)
Country | Link |
---|---|
US (2) | US6943870B2 (ja) |
EP (1) | EP1400779A1 (ja) |
JP (1) | JP3955899B2 (ja) |
KR (1) | KR100947015B1 (ja) |
Families Citing this family (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3955899B2 (ja) * | 2002-09-20 | 2007-08-08 | 国立大学法人埼玉大学 | 電子的スペックル干渉法を用いた変形計測方法および装置 |
GB2411001B (en) * | 2004-02-10 | 2007-03-28 | Statoil Asa | Seismic exploration |
JP2007071790A (ja) * | 2005-09-08 | 2007-03-22 | Saitama Univ | Memsの動作診断装置及び動作診断方法 |
US7280187B1 (en) * | 2005-10-12 | 2007-10-09 | Southeastern Louisiana University | Method for resolving phase in electronic speckle interferometry |
US8198604B2 (en) * | 2007-09-28 | 2012-06-12 | Trustees Of Boston University | System and method for providing enhanced background rejection in thick tissue with differential-aberration two-photon microscopy |
KR100944082B1 (ko) * | 2008-01-09 | 2010-02-24 | 한국원자력연구원 | 디지털 전자 스펙클 간섭을 이용한 위상 측정 방법 및 위상측정 시스템 |
US20090298206A1 (en) * | 2008-05-28 | 2009-12-03 | International Business Machines | Method and apparatus to minimize stress during reflow process |
TW201043935A (en) * | 2009-06-02 | 2010-12-16 | Univ Nat Taiwan Science Tech | Measuring system |
MD298Z (ro) * | 2010-03-03 | 2011-07-31 | Институт Прикладной Физики Академии Наук Молдовы | Sistem cu fibră optică a semnalizării de pază |
WO2011135698A1 (ja) * | 2010-04-28 | 2011-11-03 | キヤノン株式会社 | 変形計測方法 |
JP5637738B2 (ja) * | 2010-06-04 | 2014-12-10 | キヤノン株式会社 | 変形計測装置および変形計測方法 |
WO2011153973A1 (de) * | 2010-06-10 | 2011-12-15 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Verfahren zur berührungslosen, zerstörungsfreien bestimmung der härte, porosität und/oder mechanischen spannungen an werkstoffen oder verbundwerkstoffen |
US9464882B2 (en) * | 2012-04-05 | 2016-10-11 | Nanowave, Inc. | Interferometer with continuously varying path length measured in wavelengths to the reference mirror |
TWI461650B (zh) * | 2012-04-17 | 2014-11-21 | Univ Nat Taiwan | 影像處理系統及其方法 |
NO20121458A1 (no) * | 2012-12-05 | 2014-05-19 | Optonor As | Fremgangsmåte for vibrasjonsmåling og interferometer |
DE102014224852B4 (de) | 2013-12-05 | 2016-08-25 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Verfahren zur berührungslosen, zerstörungsfreien Bestimmung von Inhomogenitäten und/oder Defekten an Oberflächen von Bauteilen oder Proben |
CN104457581B (zh) * | 2014-08-28 | 2017-03-22 | 深圳奥比中光科技有限公司 | 一种全场z向位移测量*** |
WO2016145406A1 (en) * | 2015-03-11 | 2016-09-15 | Massachusetts Institute Of Technology | Methods and apparatus for modeling deformations of an object |
CN106247952B (zh) * | 2016-08-31 | 2018-11-27 | 昆明理工大学 | 一种基于傅里叶变换相位的散斑照相面内位移量算法 |
US10380745B2 (en) | 2016-09-01 | 2019-08-13 | Massachusetts Institute Of Technology | Methods and devices for measuring object motion using camera images |
KR102016740B1 (ko) * | 2017-10-11 | 2019-10-21 | 한국철도기술연구원 | 시편에 대한 다중 위치에서의 다중 요인 측정 시스템 및 그 동작 방법 |
JP7306409B2 (ja) * | 2019-01-29 | 2023-07-11 | 株式会社島津製作所 | 変位計測装置及び欠陥検出装置 |
CN110926356B (zh) * | 2019-12-05 | 2021-11-02 | 武汉科技大学 | 一种用于定型耐火材料高温形变测试的散斑制备方法 |
CN111121615B (zh) * | 2020-01-02 | 2021-10-01 | 南京大学 | 一种基于希尔伯特变换的相移干涉条纹图批量选取方法 |
US11391563B2 (en) * | 2020-02-06 | 2022-07-19 | Purdue Research Foundation | Ultra-sensitive speckle analyzing system |
US11525667B2 (en) | 2020-06-02 | 2022-12-13 | United States Of America As Represented By The Secretary Of The Air Force | Electronic speckle pattern interferometer (ESPI) for long-range measurement of displacement of materials within hazardous environments |
CN113160088B (zh) * | 2021-04-30 | 2022-08-12 | 河南大学 | 一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法 |
CN114485448B (zh) * | 2022-01-25 | 2022-11-15 | 上海大学 | 一种自适应散斑干涉测量方法及*** |
CN114964578B (zh) * | 2022-07-27 | 2022-11-15 | 华能(浙江)能源开发有限公司玉环分公司 | 基于数字散斑的水冷壁在线应力监测方法及装置 |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4594003A (en) | 1983-07-20 | 1986-06-10 | Zygo Corporation | Interferometric wavefront measurement |
US5363849A (en) * | 1994-01-26 | 1994-11-15 | Cardiovascular Imaging Systems, Inc. | Enhancing intravascular ultrasonic blood vessel image |
DE4427692A1 (de) | 1994-08-04 | 1996-02-08 | Bayerische Motoren Werke Ag | Verfahren zum Bestimmen des Schwingungsverhaltens eines Körpers |
JP2823116B2 (ja) | 1996-06-05 | 1998-11-11 | ヒーハイスト精工株式会社 | ディジタルスペックル干渉法を利用した物体の変形の動的観察方法 |
DE69713318T2 (de) * | 1997-07-21 | 2003-01-16 | European Atomic Energy Community (Euratom), Luxemburg/Luxembourg | Verfahren und Vorrichtung zum Messen der Deformation eines mechanischen Prüfstücks |
US6549801B1 (en) * | 1998-06-11 | 2003-04-15 | The Regents Of The University Of California | Phase-resolved optical coherence tomography and optical doppler tomography for imaging fluid flow in tissue with fast scanning speed and high velocity sensitivity |
US6524249B2 (en) * | 1998-11-11 | 2003-02-25 | Spentech, Inc. | Doppler ultrasound method and apparatus for monitoring blood flow and detecting emboli |
GB9826751D0 (en) * | 1998-12-05 | 1999-01-27 | British Aerospace | Method for increasing the signal to noise ratio in non-destructive testing |
JP2001311613A (ja) | 2000-02-24 | 2001-11-09 | Satoru Toyooka | スペックル干渉の画像解析における位相アンラッピング方法 |
WO2001082225A2 (en) * | 2000-04-24 | 2001-11-01 | Washington University | Method and apparatus for probabilistic model of ultrasonic images |
US6624893B1 (en) * | 2001-06-06 | 2003-09-23 | Veeco Instruments Inc. | Correction of scanning errors in interferometric profiling |
US6685641B2 (en) * | 2002-02-01 | 2004-02-03 | Siemens Medical Solutions Usa, Inc. | Plane wave scanning reception and receiver |
JP3955899B2 (ja) * | 2002-09-20 | 2007-08-08 | 国立大学法人埼玉大学 | 電子的スペックル干渉法を用いた変形計測方法および装置 |
KR101443607B1 (ko) * | 2007-05-17 | 2014-09-23 | 엘지전자 주식회사 | 세탁물 처리 장치 |
-
2002
- 2002-09-20 JP JP2002275775A patent/JP3955899B2/ja not_active Expired - Lifetime
-
2003
- 2003-01-02 US US10/334,784 patent/US6943870B2/en not_active Expired - Fee Related
- 2003-01-28 EP EP03002029A patent/EP1400779A1/en not_active Withdrawn
- 2003-02-05 US US10/358,285 patent/US20040059526A1/en not_active Abandoned
- 2003-02-11 KR KR1020030008447A patent/KR100947015B1/ko not_active IP Right Cessation
Also Published As
Publication number | Publication date |
---|---|
KR100947015B1 (ko) | 2010-03-10 |
EP1400779A1 (en) | 2004-03-24 |
US20040057054A1 (en) | 2004-03-25 |
JP2004109075A (ja) | 2004-04-08 |
US20040059526A1 (en) | 2004-03-25 |
KR20040025808A (ko) | 2004-03-26 |
US6943870B2 (en) | 2005-09-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3955899B2 (ja) | 電子的スペックル干渉法を用いた変形計測方法および装置 | |
JP6707258B2 (ja) | 応力可視化装置および力学物性値可視化装置 | |
WO2017221324A1 (ja) | 音波伝搬映像化装置及び方法 | |
US7177029B2 (en) | Stroboscopic interferometry with frequency domain analysis | |
JPH04249752A (ja) | 被検体の光学的検査方法、被検体の光学的検査装置、および被検体の光学的検査用干渉計 | |
JP3871309B2 (ja) | 位相シフト縞解析方法およびこれを用いた装置 | |
Dhanasekar et al. | Digital speckle interferometry for assessment of surface roughness | |
US8275573B1 (en) | Large-surface defect detection by single-frame spatial-carrier interferometry | |
JP5130513B2 (ja) | 3次元変位ひずみ計測方法及び装置 | |
EP0921371A2 (en) | Method of using a microscopic digital strain gauge | |
Martı́nez et al. | Fracture detection by grating moiré and in-plane ESPI techniques | |
EP2416112B1 (en) | Shape measuring apparatus | |
Cordero et al. | Uncertainty analysis of displacements measured by phase-shifting moiré interferometry | |
JP2001311613A (ja) | スペックル干渉の画像解析における位相アンラッピング方法 | |
Morimoto et al. | Theory and application of sampling moiré method | |
JP5518187B2 (ja) | 変形計測方法 | |
JP4144828B2 (ja) | 干渉縞測定解析方法 | |
Martínez et al. | Analysis of optical configurations for ESPI | |
JP4298105B2 (ja) | 干渉縞測定解析方法 | |
Quan et al. | Wavelet analysis of digital shearing speckle patterns with a temporal carrier | |
Li et al. | Matched correlation sequence analysis in temporal speckle pattern interferometry | |
Carletti et al. | Measuring velocity of speckle fields with digital speckle pattern interferometry | |
JP2006119099A (ja) | 周期可動物の変位測定装置 | |
Martınez-Celorio et al. | Out-of-plane displacement measurement by electronic speckle pattern interferometry in presence of large in-plane displacement | |
Chien et al. | Nanoscale deformation measurement by using the hybrid method of gray-level and holographic interferometry |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050602 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050620 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20070315 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20070411 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 3955899 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
EXPY | Cancellation because of completion of term |