JP2004537732A - 干渉縞を投影し絶対位相写像を評価することによる三次元画像法 - Google Patents
干渉縞を投影し絶対位相写像を評価することによる三次元画像法 Download PDFInfo
- Publication number
- JP2004537732A JP2004537732A JP2003519349A JP2003519349A JP2004537732A JP 2004537732 A JP2004537732 A JP 2004537732A JP 2003519349 A JP2003519349 A JP 2003519349A JP 2003519349 A JP2003519349 A JP 2003519349A JP 2004537732 A JP2004537732 A JP 2004537732A
- Authority
- JP
- Japan
- Prior art keywords
- phase
- frames
- fringe
- fringes
- equation
- 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.)
- Withdrawn
Links
- 238000003384 imaging method Methods 0.000 title description 4
- 238000013507 mapping Methods 0.000 title description 4
- 238000000034 method Methods 0.000 claims abstract description 39
- 238000012545 processing Methods 0.000 claims abstract description 7
- 230000010363 phase shift Effects 0.000 claims description 6
- 238000010408 sweeping Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 14
- 238000005305 interferometry Methods 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 description 18
- 238000013459 approach Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 10
- 230000003287 optical effect Effects 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 238000011105 stabilization Methods 0.000 description 2
- XKRFYHLGVUSROY-UHFFFAOYSA-N argon Substances [Ar] XKRFYHLGVUSROY-UHFFFAOYSA-N 0.000 description 1
- 229910052786 argon Inorganic materials 0.000 description 1
- -1 argon ion Chemical class 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 229910052736 halogen Inorganic materials 0.000 description 1
- 150000002367 halogens Chemical class 0.000 description 1
- CPBQJMYROZQQJC-UHFFFAOYSA-N helium neon Chemical compound [He].[Ne] CPBQJMYROZQQJC-UHFFFAOYSA-N 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000001795 light effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000004441 surface measurement Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000001931 thermography Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/2441—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures using interferometry
-
- 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/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/25—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object
-
- 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/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/25—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object
- G01B11/254—Projection of a pattern, viewing through a pattern, e.g. moiré
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
本発明は、物体の表面上の点の集合に関する三次元表面座標を計算する方法に関する。本方法は以下のステップを備える: 縞の集合を使って前記物体を照明するためのプロジェクタを使用し、前記縞を調整し、異なる縞位相設定でカメラを使用して前記表面の複数の画像を取り込み、前記画像を処理して前記プロジェクタによって照明されると共にカメラに写った前記表面の部分の絶対縞位相写像を生成し、前記縞位相写像を処理して前記物体の表面における点の座標の集合を与える。本方法は、修正されたCarre手法に基づいた暫時位相計測干渉法(TPMI)を使用することによって、表面プロファイリング・計測法を提供する。
Description
【技術分野】
【0001】
本発明は、汎用の3D表面プロファイリング・測定システム(3D surface profiling and ranging system)に関する。
【背景技術】
【0002】
3D表面プロファイリング・測定システムは、たとえば、オンライン生産管理、製品検査、ロボット製造部門、患者を長く静止させておくことが不可能な医療用用途、および大規模3D表面の計測において有用である。
【0003】
3D計測技術には大きく分けて、ステレオ画像法、レーザースキャニングおよび縞投影の三種類がある。ステレオ画像法システムの多くは、コンピュータアルゴリズムを設計して対処する、非常に特殊な状況でのみ首尾よく動作するが、大量の予備知識を必要とする。レーザースキャニングは精度は高いが、費用が高くスピードが遅い。その代替品としての縞投影は、今まで測定機能、汎用性、および縞の次数を一意に識別するための能力に欠け、それ故、接続されていないばらばらの表面投影図の複雑な形状のターゲットにおける絶対縞位相、あるいは明確な表面の特徴のないなめらかに湾曲したプロファイルを持つ大きな3D表面上の絶対縞位相が欠けていた。
【0004】
周辺光の影響または/および光ファイバからの低い光パワー(100mWまで)のために、既存の縞プロジェクタの動作距離は5メートルを超えることはない。能力の高いレーザダイオードは、パルスモードにおいて最大100Wおよびそれ以上のCWパワーを出力することができる。通常、このようなダイオードはマルチモードシングルストライプまたは空間コヒーレンスのよくないキャビティ配列で構成されているため、干渉計に応用することはできない。
【発明の開示】
【0005】
われわれは改良型3D画像法システムを発明した。
本発明によれば、物体の表面における点の集合に関する三次元表面座標を計算する方法が提供され、本方法は、縞の集合を使って前記物体を照明し、前記縞を調整し、異なる縞位相設定でカメラを使って前記表面の複数の画像を取り込み、前記画像を処理してプロジェクタによって照明されると共に前記カメラに写った前記表面の部分の絶対縞位相写像を生成し、前記縞位相写像を処理して前記物体の表面上の点に関する座標の集合を与える。
縞は、例えば、以下に述べる鏡(ロイド鏡)の手法、フレネルの複プリズム、マイケルソン干渉計等、またはマスクから投影された縞等の技術によって生成することができる。
計測処理の間、プロジェクタ、物体およびカメラは動かないように固定される。プロジェクタは照明錐体(illumination cone)内の縞の集合を使って物体を照明する。
【0006】
縞は干渉縞であってよく、生成方法によって、平行あるいは既知の等高線であってよい。これらの縞は、略同等かまたは既知の角距離分離(angular separation)を持ち、システムによって位相および空間周波数を調整することができる。カメラはプロファイリング対象の表面の縞画像を取り込む。縞を調整することによって、異なる縞位相設定で物体の画像をいくつか取り込むことができる。これらの画像は処理され、プロジェクタによって照明されると共にカメラに写った表面の部分の完全な縞位相写像を生成する。位相写像はシステムの形状の詳細と共に処理されることができ、物体の表面の各点に関する座標の集合が与えられる。データの処理手順は、以下に述べる位相解析で述べることとする。
【0007】
縞の画像は、一つはレーザ光源から(あるいはレンズまたは屈折系を介してその像から)直接到達し、他方は光源の像を生成する鏡または他の光学素子(例えば光学くさび)からの反射あるいは屈折を介して到達する二つの波の干渉によって生成することができる。この結果、遠距離フィールドにおける任意の点の強度は、二つの波の物理的経路の差、およびシステムにおける付加的な位相シフトに依存する二つの波の位相差に依存する。
【0008】
照射は、電磁スペクトル(ガンマ線から電波まで)のいずれの部分であってもよく、レーザの任意のタイプ、例えば、気体(エキシマ、アルゴンイオン、ヘリウムネオン、二酸化炭素等)レーザ、またはレーザダイオード、YAG等の固体レーザ、またはLED,ハロゲンランプ等他のレーザ光源を使用することができる。
【0009】
本発明は音響波または波動(例:水の波)など任意の他のタイプでも使うことができる。
カメラは、使用される光源の波長に適するものであればどのようなタイプ(アナログまたはデジタル)でもよい(例:ビジコン、CCD、焦電形、熱探知カメラ等)。
【0010】
縞の調整は位相ステップ単位で行われ、縞解析方法の多くは、それぞれがπ/2の領域内にある位相ステップの集合を必要とする。これらのステップはπ/2周辺で変化し、恐らくπ/6の全体領域をカバーしている。この手段は度で表され、位相ステップを75°から105°の範囲内に限定する。
【0011】
これを実現するためには、干渉計ミラーを単に動かして、光源から直接到達する光および反射されて到達する光とのパスの長さの差を変化させるという最も基本的なステッピング構成は適していないと思われる。これはステップがフィールドにおいて各点で縞の次数に比例するためである。例えば、フィールドの高次の縞の端部が105°のステップを持つ場合、ステップ角が75°より小さければフィールドのほとんど(約70%)が使用できないであろう。たとえπ/6の範囲が例えばπ/3に拡張されても、フィールドの50%が依然として無駄になる。したがって、フィールドの部分の多くを使用しないか、あるいは縞をステッピングしてフィールドにおけるステップサイズの範囲を狭くするより効率的な方法を発見する必要がある。フィールドをフルに使用する一つの方法は、以下のように混合されたステッピング方法となるであろう。
(i) 投影されたフィールドをスイ−プすることにより、全ての縞を一緒にステッピングする。これは、干渉計の後段に配置される第二のミラーを小さい角で回転させることにより実現できるであろう。代わりに、プロジェクタ全体を縞の平面に垂直な平面で回転させることも可能であろう。どちらにしても、縞の間隔はフィールドにおいてほぼ一定であるので、結果として得られる位相シフトは全ての縞に関してほぼ同じになるであろう。
(ii) さらに、干渉計ミラーをステッピングして、フィールドにおいて縞の次数に比例して変化する小さい位相シフトを付け加える。位相シフトは、画像ノイズによって、低い次数の縞においてほぼゼロから恐らくは最大約7π/6またはπ/3に変化するであろう。画像ノイズが小さい場合には、位相シフトを大きくすればよい。
(iii) したがって、全体の位相ステップサイズは、選択された最小値および最大値の間の範囲となるであろう。
【0012】
これは明らかに複雑な方法であり、われわれは、はるかに広い範囲の位相ステップサイズを実現する、改良型好適方法を発明した。
われわれは以下の数式を計算した。
【0013】
【数1】
数式1
ここで、Δdはミラーの動き、λは光の波長、δpはミラーの平面、および光源とミラー間の中間点からの点Pまでの距離との間の角である。
したがって、最大縞次数が例えば30°で存在し、レーザ波長が例えば670nmである場合、この角度における5π/6の各位相ステップに対するミラーの動きは、ワンステップにつき0.28μmになるであろう。
【0014】
6つのフレームが必要とされる場合、五つのステップが必要になり、この場合、例えば約1.4μmの全体の動きが必要となる。
システム設定の間、システムを較正し計測に使用されるdの値を評価するために、d=0までミラーをスキャンする必要がある。したがって、ミラーのステージは、少なくともdの動作値よりも大きい動きが可能となるはずである。実際の所、システムは、より長い(例えば830nm)波長における動作を可能とするはずであり、より数の多い縞(恐らくは30)を使って動作する能力を持つ。必要とされる全体の動きは、適宜な余裕も含めて、恐らくは30から40nmになるであろう。この動きは実用化されている変換ステージの能力内で十分である。
【0015】
位相をシフトしたフレームの集合の点ごとの解析は、Carreの手法と呼ばれる暫時位相計測干渉法(TPMI)(temporal phase measurement inteferometry)によって行うことができ、好ましくは修正されたCarre手法が使用される。
【0016】
基本的なCarre手法において、計測される点における位相ステップ(α)およびラップされる位相値(Ф)は以下のように求められる。
【0017】
【数2】
数式2
および
位相ステップ(a)の最適値はπ/2に近いが、位相ステップの測定は、(ノイズのないシステムでは)α>0からα>πの範囲全体にわたって安定している。
【0018】
位相値(ψ)の計算にπを加えると、位相値は−πから+πではなく0から2πの範囲内になる。
しかし、位相が上記の数式だけを使って計算された場合、結果として得られる位相値はπ/2から+3π/2の範囲内でのみ正しいということに留意されたい。これは、この範囲外の位相値に関しては、計算された逆正接値が再現するからである。2πの範囲全体にわたって正しい位相値を得るためには、計算された位相値がどの四分円に存在するかを識別する必要がある。これは、数式8の分子と分母の符号を求めることによって行うことができる。この情報を使って、連続する位相値の2πの範囲を与えるために訂正を加えることができる。
【0019】
実際には、処理ソフトウェアは四分円情報を抽出するために組み込まれた機能を持つように設計され、これにより位相値の範囲は2π全体となる。
位相および位相ステップの計算に対する基本的なCarre手法のアプローチによって、ノイズのないシステムのための正しい解析が与えられるが、データがノイズを含むという深刻な問題がある。この問題は、縞の位相(0およびπ)のある値において、位相ステップが不定になるということに起因する。ノイズのある実際のシステムでは、計算された位相ステップを評価することのできないこれらの中間点に近い縞の位相値の、問題となる範囲がある。これは、縞の位相に関する情報だけを必要とする従来のアプリケーションでは問題ではない。それらのアプリケーションにとって、位相ステップの値におけるエラーは、位相値0およびπに近い点においては、計算された縞の位相にほとんど影響を与えないという事は幸いである。
【0020】
したがって、中間領域を避けるアプローチが必要となる。一つの方法は、少なくとも一つが縞位相値0あるいはπに近くないように縞位相値が異なるような4フレームの二つ(あるいはそれ以上)の集合を使用する方法である。これは、ここで使用されたアプローチの基本であるが、その実現は、利用可能な集合から一番よい集合を単に選択するよりはるかに効率的である。
【0021】
(例えば、1、2、3、4、5と番号のついた)5つのフレームの集合を得られた場合、これらのフレームは、フレーム1、2、3、4および2、3、4、5からなる4つのフレームの二つの集合として考えることができる。使用することのできる位相ステップの範囲では、少なくとも4フレームの集合の一つが、縞位相の任意の値に対し位相ステップの不確かな値を与えるはずであるということが分かった。
【0022】
同様に、六フレームの集合を四フレームの3つの集合として考えることができる。このアプローチは、勿論任意の数のフレームに応用することができるが、ここでは6を選択することとする。これは、計算の必要条件と光学系の安定化必要条件に過分になることなく良い結果をもたらす。
【0023】
この計算の最も重要な部分は、位相ステップの評価から情報を組み合わせる方法である。この計算では、六つのフレームが4フレームの3つの集合として扱われたが、数式2におけるようにステップの値が4フレームの各集合に対して完全に評価される場合、どの結果が許容できるかを評価するためにその後結果を評価する必要があり、全体の処理はかなり面倒になる。
【0024】
【数3】
数式2
この問題は、不確かさの原因が、数式2内の下記の比に存在することに起因する。
【0025】
【数4】
式2
縞の位相値0あるいはπの近傍では、この比の分子と分母は双方とも0に近づくので、非常にノイズを受けやすくなる。したがって、不確かな領域を使わずに、4フレームの3つの集合からの比を組み合わせることのできる方法が求められた。4フレームの3つの集合にa、b、cと名前をつけた場合、関連する比を以下のように表すことができる。
【0026】
【数5】
位相ステップサイズは、4フレームの各集合に関し、与えられた点において同じであるので、それらのフレームも、組み合わされた以下の比に等しいことが分かる。
【0027】
【数6】
組み合わされた比が示されたように計算された場合、不確かな領域は縞の位相の異なる値に単純にシフトする。これは、AとBの値が正あるいは負であり、それらの和がゼロであるからである。しかし、この応用においては、(位相ステップが各集合で同じであるため)集合a,b,cの比は全て正か負であることに留意されたい。したがって、正の位相ステップ(>0および<π)では、AとBは両方正または負であり、これにより組み合わされた比は、任意のAおよびBの組の符号が両方変化しても変わらない。この結果、AおよびBの各値の係数(?)が組み合わされた比において以下のように使用される場合、この結果が影響されることはない。
【0028】
【数7】
数式3
縞のコントラストまたは位相ステップがゼロに近づく場合、数式3における比の分母は単に0に近づく。これらの条件は両方とも避けられる。図4または5に示されるように、この事が位相ステップ計算に対して与える影響は、実際の画像データの集合の一つのラインから計算された以下(図4と5)のグラフにおいて分かる。
【0029】
図4および5に関する計算では、約4.5の縞にわたり、かつ画像において大きくなる位相ステップを持つフレームの同じ集合からのデータが使われている。計算の結果は、図4および5のグラフにおける9つの不確かな領域となる。ターゲットの物体における物理的ステップに対応する不連続性は図5において明らかに分かり、この場合、縞次数の約三分の一の絶対位相における変化を表す。
【0030】
これによって、位相ステップ計算のための6フレームの数式は、以下のようになる。
【0031】
【数8】
数式4
ここで、
Aa=3(I2−I3)−(I1−I4)
Ba=(I2−I3)+(I1−I4)
Ab=3(I2−I4)−(I2−I5)
Bb=(I2−I4)+(I2−I5)
Ac=3(I4−I5)−(I3−I4)
Bc=(I4−I5)+(I3−I4)
であり、I1からI6は6フレームにおける画像点の強度pfである。
【0032】
位相ステップサイズが不確かさなしに決定できることになったので、縞の位相値に注目してみる。幸い、縞の位相計算はステップ計算よりも問題が少ない。実際、数式3における元の公式は、数式4からの位相ステップのための6フレームの値と共に、強度I2,I3,I4,I5を使って使用することができる。これは有用な結果をもたらすが、利用できるデータの6つのフレームの全てを使用している訳ではない。以下の6フレームの数式は、縞の位相値に関して、多少良い結果をもたらす。
【0033】
【数9】
数式5
絶対位相値計算
原則は、プロジェクタを好適に設計し、確実に位相ステップサイズを絶対縞位相に比例させるということである。各画像点に関して、われわれは二つの値を計算する。これらは位相ステップサイズと、ラップした縞の位相(ψの0から2πまでの再現範囲)である。この画像点における縞の絶対位相を計算するために、これら二つの値に関する比例関係Sの定数を知る必要がある。
【0034】
【数10】
数式6
もっとも単純なアプローチでは、位相ステップは以下のように絶対縞位相の概算値を計算するために使うことができるであろう。Ф=S*a。このように求められた値は別の縞位相計算を全く必要としないが、ラップされた位相値より(ノイズに起因する)はるかに大きい局所的に不確かな値となる。ラップされた位相値は、不明の縞次数n(絶対縞位相値Ф=2nπ+ψであるので、縞位相エラー=2nπ)の不確かさを当然持つこととなる。
【0035】
Фの概算値はラップされた位相値ψと比較することができ、調整されて最終的な絶対位相が与えられる。これは単にФの概算値=sαおよびより正確な値Ф=2nπ+ψとの間で一番適するnの値を選択することである。このアプローチは解析ソフトウェアで使用され、計算された位相ステップ値のエラーが、πより小さい絶対位相エラーに対応するという条件下で機能し、その結果縞の次数が正しく識別される。
【0036】
勿論、データを解析する前にパラメータSを求める必要がある。これは、別の較正処理における既知の縞次数に関する位相ステップを計測することによって直接計測される。較正の精度は、光システムの継続する安定度に依存するであろう。しかし、物体の六つの画像のどの集合も、各計測に関するパラメータSを較正するために必要な全ての情報を既に含んでいる。
【0037】
Sは次のようにして求めることができる。画像におけるいくつかまたは全てのデータポイントを較正に使うことができる。この手順では、計算された各位相ステップの値と対応する計算されたラップ位相値とが比較される。これらの値の組がそれぞれ互いにグラフ化されるとその結果は図6のようになる。
【0038】
このグラフは移動されたラインセグメントの集合で構成される。与えられたいずれの画像のグラフも、物体の表面のプロファイルによって、特徴の一部が失われているが、これは解析に影響を与えず、いくつかの異なる縞の次数による利点がもたらされる。Sを求めることは、隣接するラインセグメント間の中間位相ステップの差を決定することである。これは、同じようなラップされた位相値を持つ二つの点の計測された位相ステップの間の差を計算することによって行われる。この一組の点は、(ラップされた位相の値が非常に近いので)、2πに近い絶対位相差を使って二つの絶対位相値を表す。ここでn=0,1,2,3とする。これらの差の値の集合が(例えば昇順で)サンプリングされグラフ化された場合、図7のようなグラフが得られる。
【0039】
そして、帯域の上限および下限の閾値の組を選択することによって、(例えば)n=1を表す位相ステップ差の値を求めることができる。選択された帯域内の点の平均値が次に計算される。与えられたシステムでは、閾値は固定され、適宜なステップ差の値が閾値の限界内に収まっているという条件下でSを求めることができる。n=1の値がα1であり、隣接する次数間の絶対位相差が2πであることが与えられている場合、数式6からのSの値は下記の式によって与えられる。
【0040】
【数11】
動作パラメータの選択
上記したように、絶対位相(Ф)の計算では、位相ステップaを、計測される点を含む縞の次数を識別するために十分な精度で決定する必要がある。Фの評価を司るパラメータの選択は、縞の次数を誤って識別する確率を許容できるレベルにまで低減しながら、位相計測の精度を最適化するよう注意深く行わなければならない。任意の点における位相ステップサイズの不確かさの範囲は、ある縞の次数と次の次数の間の位相ステップの差を表す閾値よりも小さくなければならない。この不確かさは、下記に示すように、画像における縞の数、および位相ステップサイズの範囲に大きく依存する。
【0041】
位相計測の精度を最適化するための必要条件は以下のとおりである。
― 投影されたフィールドにおける縞の数を最大にする。
― 位相ステップサイズの範囲を最小ノイズ感度値に近い値まで小さくする。
― 画像ノイズを最小にする。
縞次数を誤って識別する確率を最小にするための必要条件は以下のとおりである。
― 投影されたフィールドにおける縞の数を最小にする。
― 位相ステップサイズの範囲を最大にする。
― 画像ノイズを最小にする。
【0042】
したがって、われわれは画像ノイズを最小にし、他の二つのパラメータの最適と思われる値を決定する必要がある。
表面形状計算
表面の位相写像により、画像平面の各ピクセルに関する絶対縞位相値が与えられる。これは、表面座標の写像ではない。計算された位相値と光システムの形状を使って、各ピクセルによって画像化された点の3D座標を求めることができる。
【発明を実施するための最良の形態】
【0043】
本発明による実施形態の光システムの形状の説明、および表面座標の写像を生成するために必要な計算は、添付の図面に図示されている。
図1および3において、物体(3)上の点(P)の座標を求めるための基本的構成では、位置Sにおけるプロジェクタ(1)は縞を物体(3)に投影し、また位置(C)(図3)におけるカメラ(2)は物体の画像を取り込む。図3において、x軸は位置Sにおけるプロジェクタと位置(C)におけるカメラとの間の線であり、z軸はカメラの軸であり、y軸はこれらの軸に直交している。
【0044】
プロジェクタが示されている図2では、照明されている物体が(10)に示されており。レーザ(5)からの光がレンズ(6)および(7)を通り、レーザストライプ(8)をY−Zステージ(9)上のミラーに投影している。その結果、一つはレーザから直接到達し、他方はミラーを介して到達する二本の光線が画像上に投影され、干渉縞を形成する。ステージ(9)を動かしステップ単位にミラーを動かすことにより、二つの光線間の位相差を調整し、これによって干渉縞を調整することができる。
【0045】
表面座標の計算がかなり単純になれば、座標系の選択は注意深く考えられなければならない。基本的なシステムの構成と向きは、上記したように図1および3に示されている。図8は、選択したX軸とZ軸、および関連する系の次元と角を示している。Y軸は図の平面に垂直である。
【0046】
図8は、タイル平面Y=0の構成を示す、しかし、縞は“垂直”、すなわちYにおいて一定であるので、角θA、θS、θXは点PのY座標から独立している。したがって、点PのXおよびZ座標は、Y座標の計算により、任意のYに関して計算することができる。
【0047】
上記の図において、SS’は光源(プロジェクタ)の軸であり、ゼロ次数の縞に対応する方向である。投影された全ての縞はプロジェクタ軸の一つの側にあるので、θSの全ての値は正である。角θSは、以下に述べるように絶対位相Фおよび縞の角距離から求められる。
【0048】
CC’はカメラ軸であり、座標系のz軸に一致する。θXは正および負両方の値を取る。
カメラおよびプロジェクタは、X軸に沿った距離Dの分だけ離れている。
座標の計算
このセクションは、点Pのx、y、z座標を、入力されたパラメータから求める方法について述べる。
【0049】
図8において、カメラ(C)から見ると、
x=z.Tan(θX)
数式7
プロジェクタ(S)から見ると、
【0050】
【数12】
数式8
【0051】
【数13】
数式9
xおよびz座標は、タイル入力パラメータでは以下のように表すことができる。
【0052】
【数14】
数式10
ここで、xは数式7から直接求めることができる。
【0053】
x=z.Tan(θx)
最後に、y座標は次のように表すことができる。
【0054】
【数15】
ここで、θYは、図9に示すように、x−z平面とカメラから物体の点Pまでのタイル方向との間の角として定義される、角θXおよびθYは、以下のように画像のピクセル位置に関連している。
【0055】
画像平面の中央ピクセルがz軸に沿ったビューに対応するようにシステムが位置決めされたと仮定する。カメラが各行においてNx個のピクセルを持つ場合、列は以下のように番号がつけられる。
【0056】
【数16】
ピクセルの位置は、z軸からの角θXにおけるビューに対応する。ピクセル距離がwであり、カメラレンズが焦点距離vを持つ場合、ピクセル番号nxに対するθXは、
【0057】
【数17】
同様に、カメラが各列においてNy個のピクセルを持つ場合、ピクセル番号nyに対するθYは、
【0058】
【数18】
点Pに関する角θSは、絶対縞位相Фに関係する。数式1から
Ф=C*Sin(θS)
ここでCは定数である。
【0059】
Cの値は、縞の次数(またこれによりФ)を一つの角θSに対して計測する別の較正手順から得られる。勿論、Cの実際の値は縞がスキャンされるにつれて変化する。したがって、フレームの集合の中間位置に設定された縞を使ってシステムを較正することが重要である。この中間位置は、位相計算から得られる位置と同じである。
【0060】
角θAはシステムを位置決めすることによって直接設定される。Sにおけるプロジェクタが、カメラの向きによって定義されたようにx軸上に存在することが、有効な解析にとって不可欠である。
【0061】
図10においては、10aが正面図、10bが側面図、10cが平面図であり、縞プロジェクタは、レーザダイオード(26)、ミラー(25)、および圧電作動装置で構成され、レーザキャビティとミラーの表面との間の距離dを調整する。投影された縞パターンは、レーザキャビティSによって形成されたYoungの干渉縞としてみなすことができ、その画像はミラーS’に存在する。ここで、PSおよびPS’間の光路の差は2dsinδpに等しく、δpはミラーの表面の方向およびPとプロジェクタを結ぶ線の方向との間の角である。SからS’まで出力された光の間のPにおける”大域的“位相差Фpは、以下の数式によって与えられる。
【0062】
【数19】
数式11
ここで、λは投影された光の波長である。Pにおける局所位相、すなわち0から2πの不確かな範囲内のФpの値はこれによりφpとして表される。
【0063】
Фp=φp+2Nπ
数式12
ここで、NはPにおける干渉縞の次数である。φpを計測することのできない暫時位相計測干渉法(TPMI)(19)と呼ばれる確立された手法が数多く存在し、これらの干渉法では、制御された位相ステップであるαpjがФpに加えられ、各ステップに対して点Pにおける対応するインターフェログラムの強度が以下のように記録される。
【0064】
lpj=Ap+Bpcos(Фp+αpj)
数式13
ここで、添え字jは位相ステップのシーケンス番号を表す。Carreの手法は、特に、導入されるべき4つの等しいステップαp、1、2、3、4=−3/2αp,−1/2αp、1/2αp、3/2αpを必要とする。これにより下記の数式が与えられる。
【0065】
【数20】
数式14
従来のシステムでは、位相ステップαpは、投影された光全体において同じである。しかし、アルゴリズム自体は統一されている必要があまりないので、ここで提案されたシステムのPにおける縞の干渉次数を識別するためにαpが使われる。ここで、位相ステッピングはΔdの等しいステップ単位でdを変更することによって実現され、数式11にしたがって下記の位相ステップを生成する。
【0066】
【数21】
数式15
数式15と14を組み合わせると、以下のようにФPの概算が与えられ、この概算から縞の干渉次数を以下のように計算することができる。
【0067】
【数22】
数式16
ここでtrunc()は整数に切り捨てることを表す。数式15および14を12に導入すると、Фpのより正確な計測を得ることができる。3D空間の正確な位置Pを絶対グローバルФpを使って不確かさなしに識別することができる。さらに、数式14から|αp|<πであることが分かり、sinδP<sinθであることが分かる。これらの条件を数式15に応用する。各ステップにおけるミラーの変位は次のように推定される。
【0068】
【数23】
数式16
これは、圧電作動装置の範囲内に収まる。
【0069】
図11において、縞プロジェクタのための別の装置が示されており、図11aはその正面図、11bは側面図、11cは平面図である。ここで、二つのレーザダイオード(15)はミラー(16)上に並べられ恒久的に固定されている。これら二つの商d/λは、二つのレーザに関連する大域的位相が物体の表面上の同じ点Pでわずかに異なるように設計され、以下のように表すことができる。
【0070】
【数24】
数式17
ここで、一度にレーザAとBを切り替えることによってより成熟した3フレーム手法を使用し、関連する局所位相ФpAおよびФpBを得ることができる。これは以下のように表される。
【0071】
【数25】
数式18
ここで、制御された一連の小さい傾きをミラーレーザに導入することによって位相ステッピングを実現することができ、縞次数を以下のように計算することができる。
【0072】
【数26】
数式19
同様に、数式(12)と(9)を(2)に導入することによって、大域的位相Фp(この場合ФpA)のより正確な計測を得ることができる。
【0073】
この別の装置では画像を得る時間が少し長くかかるが(合計6フレーム。前の装置では4フレーム)、データ処理がはるかに早いと言う利点がある。これは、Carreの手法と比較すると、位相検索アルゴリズムがより簡単ではるかに成熟しているからであり、これによってより高い精度と恐らくは短い3Dプロファイルサイクルを実現する。さらに、ダイオード−ミラー組立体の全体を大量生産でひとつの基板上に作成することができ、これによってプロジェクタの堅牢さを高め、縞の質を改善して形状と大きさを小さくする。
【0074】
上記した二つの装置の共通の特徴は、レーザダイオード、特にパワーの高いレーザダイオードが光源として使用されることである。これらのダイオードは、パルスモードにおいて最大100Wの光パワーを生成することができ、計測対象の物体の距離および大きさによって、非常に広い範囲にわたって簡単に調整することができる。このようなダイオードはキャビティ配列で構成されるので、通常干渉計に応用することができず、空間コヒーレンスが劣化する。しかし、キャビティ配列がミラーの表面に平行に位置決めされる限り、これは提案されたプロジェクタにおいて問題にはならない。典型的なスペクトルラインの幅である12nmでは、これらのダイオードも非常に短い一時的なコヒーレンスの長さ(約70nm)を持っている。提案されたプロジェクタの固有の構成では、このように短いコヒーレンスの長さは十分に質の高い干渉縞を生成するのに十分なだけでなく、システムに対する特定の抑制(suppression)という利点をさらにもたらす。両方の構成の構造の単純さは、システムの計測精度にとって重要である縞パターンの安定化の改善にも役立つ。
【0075】
さらに、レーザダイオードは非常に効率的な素子である。図1に示されたように、レーザダイオードのライン幅に適するバンドパス干渉フィルタがカメラの絞りに設けられる。カメラのシャッタに対して同期されたレーザパルスと共に、フィッタがシステムの周囲光に対する耐性を少なくとも600倍高めることができる。したがって100Wパルスダイオードは、通常の白色光スライドプロジェクタの60KW電球に相当する。プロジェクタの出力光は非常によく拡散するので、このようなパワーレベルは、設計された動作距離における人間の視覚の安全な限界内に依然として収まる。波長およびパルス幅を注意深く選択することによって、安全性をさらに改善することができる。
【図面の簡単な説明】
【0076】
【図1】図1は、プロファイリングシステムを示す図である。
【図2】図2は、プロジェクタの構成を示す図である。
【図3】図3は、プロファイリングシステムの平面図である。
【図4】図4は、上記で言及された図である。
【図5】図5は、上記で言及された図である。
【図6】図6は、上記で言及された図である。
【図7】図7は、上記で言及された図である。
【図8】図8は、平面Y=0での位相の構成を表す図である。
【図9】図9は、座標の計算を示す図である。
【図10】図10は、装置の一例を示す図である。
【図11】図11は、装置の他の例を示す図である。
【0001】
本発明は、汎用の3D表面プロファイリング・測定システム(3D surface profiling and ranging system)に関する。
【背景技術】
【0002】
3D表面プロファイリング・測定システムは、たとえば、オンライン生産管理、製品検査、ロボット製造部門、患者を長く静止させておくことが不可能な医療用用途、および大規模3D表面の計測において有用である。
【0003】
3D計測技術には大きく分けて、ステレオ画像法、レーザースキャニングおよび縞投影の三種類がある。ステレオ画像法システムの多くは、コンピュータアルゴリズムを設計して対処する、非常に特殊な状況でのみ首尾よく動作するが、大量の予備知識を必要とする。レーザースキャニングは精度は高いが、費用が高くスピードが遅い。その代替品としての縞投影は、今まで測定機能、汎用性、および縞の次数を一意に識別するための能力に欠け、それ故、接続されていないばらばらの表面投影図の複雑な形状のターゲットにおける絶対縞位相、あるいは明確な表面の特徴のないなめらかに湾曲したプロファイルを持つ大きな3D表面上の絶対縞位相が欠けていた。
【0004】
周辺光の影響または/および光ファイバからの低い光パワー(100mWまで)のために、既存の縞プロジェクタの動作距離は5メートルを超えることはない。能力の高いレーザダイオードは、パルスモードにおいて最大100Wおよびそれ以上のCWパワーを出力することができる。通常、このようなダイオードはマルチモードシングルストライプまたは空間コヒーレンスのよくないキャビティ配列で構成されているため、干渉計に応用することはできない。
【発明の開示】
【0005】
われわれは改良型3D画像法システムを発明した。
本発明によれば、物体の表面における点の集合に関する三次元表面座標を計算する方法が提供され、本方法は、縞の集合を使って前記物体を照明し、前記縞を調整し、異なる縞位相設定でカメラを使って前記表面の複数の画像を取り込み、前記画像を処理してプロジェクタによって照明されると共に前記カメラに写った前記表面の部分の絶対縞位相写像を生成し、前記縞位相写像を処理して前記物体の表面上の点に関する座標の集合を与える。
縞は、例えば、以下に述べる鏡(ロイド鏡)の手法、フレネルの複プリズム、マイケルソン干渉計等、またはマスクから投影された縞等の技術によって生成することができる。
計測処理の間、プロジェクタ、物体およびカメラは動かないように固定される。プロジェクタは照明錐体(illumination cone)内の縞の集合を使って物体を照明する。
【0006】
縞は干渉縞であってよく、生成方法によって、平行あるいは既知の等高線であってよい。これらの縞は、略同等かまたは既知の角距離分離(angular separation)を持ち、システムによって位相および空間周波数を調整することができる。カメラはプロファイリング対象の表面の縞画像を取り込む。縞を調整することによって、異なる縞位相設定で物体の画像をいくつか取り込むことができる。これらの画像は処理され、プロジェクタによって照明されると共にカメラに写った表面の部分の完全な縞位相写像を生成する。位相写像はシステムの形状の詳細と共に処理されることができ、物体の表面の各点に関する座標の集合が与えられる。データの処理手順は、以下に述べる位相解析で述べることとする。
【0007】
縞の画像は、一つはレーザ光源から(あるいはレンズまたは屈折系を介してその像から)直接到達し、他方は光源の像を生成する鏡または他の光学素子(例えば光学くさび)からの反射あるいは屈折を介して到達する二つの波の干渉によって生成することができる。この結果、遠距離フィールドにおける任意の点の強度は、二つの波の物理的経路の差、およびシステムにおける付加的な位相シフトに依存する二つの波の位相差に依存する。
【0008】
照射は、電磁スペクトル(ガンマ線から電波まで)のいずれの部分であってもよく、レーザの任意のタイプ、例えば、気体(エキシマ、アルゴンイオン、ヘリウムネオン、二酸化炭素等)レーザ、またはレーザダイオード、YAG等の固体レーザ、またはLED,ハロゲンランプ等他のレーザ光源を使用することができる。
【0009】
本発明は音響波または波動(例:水の波)など任意の他のタイプでも使うことができる。
カメラは、使用される光源の波長に適するものであればどのようなタイプ(アナログまたはデジタル)でもよい(例:ビジコン、CCD、焦電形、熱探知カメラ等)。
【0010】
縞の調整は位相ステップ単位で行われ、縞解析方法の多くは、それぞれがπ/2の領域内にある位相ステップの集合を必要とする。これらのステップはπ/2周辺で変化し、恐らくπ/6の全体領域をカバーしている。この手段は度で表され、位相ステップを75°から105°の範囲内に限定する。
【0011】
これを実現するためには、干渉計ミラーを単に動かして、光源から直接到達する光および反射されて到達する光とのパスの長さの差を変化させるという最も基本的なステッピング構成は適していないと思われる。これはステップがフィールドにおいて各点で縞の次数に比例するためである。例えば、フィールドの高次の縞の端部が105°のステップを持つ場合、ステップ角が75°より小さければフィールドのほとんど(約70%)が使用できないであろう。たとえπ/6の範囲が例えばπ/3に拡張されても、フィールドの50%が依然として無駄になる。したがって、フィールドの部分の多くを使用しないか、あるいは縞をステッピングしてフィールドにおけるステップサイズの範囲を狭くするより効率的な方法を発見する必要がある。フィールドをフルに使用する一つの方法は、以下のように混合されたステッピング方法となるであろう。
(i) 投影されたフィールドをスイ−プすることにより、全ての縞を一緒にステッピングする。これは、干渉計の後段に配置される第二のミラーを小さい角で回転させることにより実現できるであろう。代わりに、プロジェクタ全体を縞の平面に垂直な平面で回転させることも可能であろう。どちらにしても、縞の間隔はフィールドにおいてほぼ一定であるので、結果として得られる位相シフトは全ての縞に関してほぼ同じになるであろう。
(ii) さらに、干渉計ミラーをステッピングして、フィールドにおいて縞の次数に比例して変化する小さい位相シフトを付け加える。位相シフトは、画像ノイズによって、低い次数の縞においてほぼゼロから恐らくは最大約7π/6またはπ/3に変化するであろう。画像ノイズが小さい場合には、位相シフトを大きくすればよい。
(iii) したがって、全体の位相ステップサイズは、選択された最小値および最大値の間の範囲となるであろう。
【0012】
これは明らかに複雑な方法であり、われわれは、はるかに広い範囲の位相ステップサイズを実現する、改良型好適方法を発明した。
われわれは以下の数式を計算した。
【0013】
【数1】
数式1
ここで、Δdはミラーの動き、λは光の波長、δpはミラーの平面、および光源とミラー間の中間点からの点Pまでの距離との間の角である。
したがって、最大縞次数が例えば30°で存在し、レーザ波長が例えば670nmである場合、この角度における5π/6の各位相ステップに対するミラーの動きは、ワンステップにつき0.28μmになるであろう。
【0014】
6つのフレームが必要とされる場合、五つのステップが必要になり、この場合、例えば約1.4μmの全体の動きが必要となる。
システム設定の間、システムを較正し計測に使用されるdの値を評価するために、d=0までミラーをスキャンする必要がある。したがって、ミラーのステージは、少なくともdの動作値よりも大きい動きが可能となるはずである。実際の所、システムは、より長い(例えば830nm)波長における動作を可能とするはずであり、より数の多い縞(恐らくは30)を使って動作する能力を持つ。必要とされる全体の動きは、適宜な余裕も含めて、恐らくは30から40nmになるであろう。この動きは実用化されている変換ステージの能力内で十分である。
【0015】
位相をシフトしたフレームの集合の点ごとの解析は、Carreの手法と呼ばれる暫時位相計測干渉法(TPMI)(temporal phase measurement inteferometry)によって行うことができ、好ましくは修正されたCarre手法が使用される。
【0016】
基本的なCarre手法において、計測される点における位相ステップ(α)およびラップされる位相値(Ф)は以下のように求められる。
【0017】
【数2】
数式2
および
位相ステップ(a)の最適値はπ/2に近いが、位相ステップの測定は、(ノイズのないシステムでは)α>0からα>πの範囲全体にわたって安定している。
【0018】
位相値(ψ)の計算にπを加えると、位相値は−πから+πではなく0から2πの範囲内になる。
しかし、位相が上記の数式だけを使って計算された場合、結果として得られる位相値はπ/2から+3π/2の範囲内でのみ正しいということに留意されたい。これは、この範囲外の位相値に関しては、計算された逆正接値が再現するからである。2πの範囲全体にわたって正しい位相値を得るためには、計算された位相値がどの四分円に存在するかを識別する必要がある。これは、数式8の分子と分母の符号を求めることによって行うことができる。この情報を使って、連続する位相値の2πの範囲を与えるために訂正を加えることができる。
【0019】
実際には、処理ソフトウェアは四分円情報を抽出するために組み込まれた機能を持つように設計され、これにより位相値の範囲は2π全体となる。
位相および位相ステップの計算に対する基本的なCarre手法のアプローチによって、ノイズのないシステムのための正しい解析が与えられるが、データがノイズを含むという深刻な問題がある。この問題は、縞の位相(0およびπ)のある値において、位相ステップが不定になるということに起因する。ノイズのある実際のシステムでは、計算された位相ステップを評価することのできないこれらの中間点に近い縞の位相値の、問題となる範囲がある。これは、縞の位相に関する情報だけを必要とする従来のアプリケーションでは問題ではない。それらのアプリケーションにとって、位相ステップの値におけるエラーは、位相値0およびπに近い点においては、計算された縞の位相にほとんど影響を与えないという事は幸いである。
【0020】
したがって、中間領域を避けるアプローチが必要となる。一つの方法は、少なくとも一つが縞位相値0あるいはπに近くないように縞位相値が異なるような4フレームの二つ(あるいはそれ以上)の集合を使用する方法である。これは、ここで使用されたアプローチの基本であるが、その実現は、利用可能な集合から一番よい集合を単に選択するよりはるかに効率的である。
【0021】
(例えば、1、2、3、4、5と番号のついた)5つのフレームの集合を得られた場合、これらのフレームは、フレーム1、2、3、4および2、3、4、5からなる4つのフレームの二つの集合として考えることができる。使用することのできる位相ステップの範囲では、少なくとも4フレームの集合の一つが、縞位相の任意の値に対し位相ステップの不確かな値を与えるはずであるということが分かった。
【0022】
同様に、六フレームの集合を四フレームの3つの集合として考えることができる。このアプローチは、勿論任意の数のフレームに応用することができるが、ここでは6を選択することとする。これは、計算の必要条件と光学系の安定化必要条件に過分になることなく良い結果をもたらす。
【0023】
この計算の最も重要な部分は、位相ステップの評価から情報を組み合わせる方法である。この計算では、六つのフレームが4フレームの3つの集合として扱われたが、数式2におけるようにステップの値が4フレームの各集合に対して完全に評価される場合、どの結果が許容できるかを評価するためにその後結果を評価する必要があり、全体の処理はかなり面倒になる。
【0024】
【数3】
数式2
この問題は、不確かさの原因が、数式2内の下記の比に存在することに起因する。
【0025】
【数4】
式2
縞の位相値0あるいはπの近傍では、この比の分子と分母は双方とも0に近づくので、非常にノイズを受けやすくなる。したがって、不確かな領域を使わずに、4フレームの3つの集合からの比を組み合わせることのできる方法が求められた。4フレームの3つの集合にa、b、cと名前をつけた場合、関連する比を以下のように表すことができる。
【0026】
【数5】
位相ステップサイズは、4フレームの各集合に関し、与えられた点において同じであるので、それらのフレームも、組み合わされた以下の比に等しいことが分かる。
【0027】
【数6】
組み合わされた比が示されたように計算された場合、不確かな領域は縞の位相の異なる値に単純にシフトする。これは、AとBの値が正あるいは負であり、それらの和がゼロであるからである。しかし、この応用においては、(位相ステップが各集合で同じであるため)集合a,b,cの比は全て正か負であることに留意されたい。したがって、正の位相ステップ(>0および<π)では、AとBは両方正または負であり、これにより組み合わされた比は、任意のAおよびBの組の符号が両方変化しても変わらない。この結果、AおよびBの各値の係数(?)が組み合わされた比において以下のように使用される場合、この結果が影響されることはない。
【0028】
【数7】
数式3
縞のコントラストまたは位相ステップがゼロに近づく場合、数式3における比の分母は単に0に近づく。これらの条件は両方とも避けられる。図4または5に示されるように、この事が位相ステップ計算に対して与える影響は、実際の画像データの集合の一つのラインから計算された以下(図4と5)のグラフにおいて分かる。
【0029】
図4および5に関する計算では、約4.5の縞にわたり、かつ画像において大きくなる位相ステップを持つフレームの同じ集合からのデータが使われている。計算の結果は、図4および5のグラフにおける9つの不確かな領域となる。ターゲットの物体における物理的ステップに対応する不連続性は図5において明らかに分かり、この場合、縞次数の約三分の一の絶対位相における変化を表す。
【0030】
これによって、位相ステップ計算のための6フレームの数式は、以下のようになる。
【0031】
【数8】
数式4
ここで、
Aa=3(I2−I3)−(I1−I4)
Ba=(I2−I3)+(I1−I4)
Ab=3(I2−I4)−(I2−I5)
Bb=(I2−I4)+(I2−I5)
Ac=3(I4−I5)−(I3−I4)
Bc=(I4−I5)+(I3−I4)
であり、I1からI6は6フレームにおける画像点の強度pfである。
【0032】
位相ステップサイズが不確かさなしに決定できることになったので、縞の位相値に注目してみる。幸い、縞の位相計算はステップ計算よりも問題が少ない。実際、数式3における元の公式は、数式4からの位相ステップのための6フレームの値と共に、強度I2,I3,I4,I5を使って使用することができる。これは有用な結果をもたらすが、利用できるデータの6つのフレームの全てを使用している訳ではない。以下の6フレームの数式は、縞の位相値に関して、多少良い結果をもたらす。
【0033】
【数9】
数式5
絶対位相値計算
原則は、プロジェクタを好適に設計し、確実に位相ステップサイズを絶対縞位相に比例させるということである。各画像点に関して、われわれは二つの値を計算する。これらは位相ステップサイズと、ラップした縞の位相(ψの0から2πまでの再現範囲)である。この画像点における縞の絶対位相を計算するために、これら二つの値に関する比例関係Sの定数を知る必要がある。
【0034】
【数10】
数式6
もっとも単純なアプローチでは、位相ステップは以下のように絶対縞位相の概算値を計算するために使うことができるであろう。Ф=S*a。このように求められた値は別の縞位相計算を全く必要としないが、ラップされた位相値より(ノイズに起因する)はるかに大きい局所的に不確かな値となる。ラップされた位相値は、不明の縞次数n(絶対縞位相値Ф=2nπ+ψであるので、縞位相エラー=2nπ)の不確かさを当然持つこととなる。
【0035】
Фの概算値はラップされた位相値ψと比較することができ、調整されて最終的な絶対位相が与えられる。これは単にФの概算値=sαおよびより正確な値Ф=2nπ+ψとの間で一番適するnの値を選択することである。このアプローチは解析ソフトウェアで使用され、計算された位相ステップ値のエラーが、πより小さい絶対位相エラーに対応するという条件下で機能し、その結果縞の次数が正しく識別される。
【0036】
勿論、データを解析する前にパラメータSを求める必要がある。これは、別の較正処理における既知の縞次数に関する位相ステップを計測することによって直接計測される。較正の精度は、光システムの継続する安定度に依存するであろう。しかし、物体の六つの画像のどの集合も、各計測に関するパラメータSを較正するために必要な全ての情報を既に含んでいる。
【0037】
Sは次のようにして求めることができる。画像におけるいくつかまたは全てのデータポイントを較正に使うことができる。この手順では、計算された各位相ステップの値と対応する計算されたラップ位相値とが比較される。これらの値の組がそれぞれ互いにグラフ化されるとその結果は図6のようになる。
【0038】
このグラフは移動されたラインセグメントの集合で構成される。与えられたいずれの画像のグラフも、物体の表面のプロファイルによって、特徴の一部が失われているが、これは解析に影響を与えず、いくつかの異なる縞の次数による利点がもたらされる。Sを求めることは、隣接するラインセグメント間の中間位相ステップの差を決定することである。これは、同じようなラップされた位相値を持つ二つの点の計測された位相ステップの間の差を計算することによって行われる。この一組の点は、(ラップされた位相の値が非常に近いので)、2πに近い絶対位相差を使って二つの絶対位相値を表す。ここでn=0,1,2,3とする。これらの差の値の集合が(例えば昇順で)サンプリングされグラフ化された場合、図7のようなグラフが得られる。
【0039】
そして、帯域の上限および下限の閾値の組を選択することによって、(例えば)n=1を表す位相ステップ差の値を求めることができる。選択された帯域内の点の平均値が次に計算される。与えられたシステムでは、閾値は固定され、適宜なステップ差の値が閾値の限界内に収まっているという条件下でSを求めることができる。n=1の値がα1であり、隣接する次数間の絶対位相差が2πであることが与えられている場合、数式6からのSの値は下記の式によって与えられる。
【0040】
【数11】
動作パラメータの選択
上記したように、絶対位相(Ф)の計算では、位相ステップaを、計測される点を含む縞の次数を識別するために十分な精度で決定する必要がある。Фの評価を司るパラメータの選択は、縞の次数を誤って識別する確率を許容できるレベルにまで低減しながら、位相計測の精度を最適化するよう注意深く行わなければならない。任意の点における位相ステップサイズの不確かさの範囲は、ある縞の次数と次の次数の間の位相ステップの差を表す閾値よりも小さくなければならない。この不確かさは、下記に示すように、画像における縞の数、および位相ステップサイズの範囲に大きく依存する。
【0041】
位相計測の精度を最適化するための必要条件は以下のとおりである。
― 投影されたフィールドにおける縞の数を最大にする。
― 位相ステップサイズの範囲を最小ノイズ感度値に近い値まで小さくする。
― 画像ノイズを最小にする。
縞次数を誤って識別する確率を最小にするための必要条件は以下のとおりである。
― 投影されたフィールドにおける縞の数を最小にする。
― 位相ステップサイズの範囲を最大にする。
― 画像ノイズを最小にする。
【0042】
したがって、われわれは画像ノイズを最小にし、他の二つのパラメータの最適と思われる値を決定する必要がある。
表面形状計算
表面の位相写像により、画像平面の各ピクセルに関する絶対縞位相値が与えられる。これは、表面座標の写像ではない。計算された位相値と光システムの形状を使って、各ピクセルによって画像化された点の3D座標を求めることができる。
【発明を実施するための最良の形態】
【0043】
本発明による実施形態の光システムの形状の説明、および表面座標の写像を生成するために必要な計算は、添付の図面に図示されている。
図1および3において、物体(3)上の点(P)の座標を求めるための基本的構成では、位置Sにおけるプロジェクタ(1)は縞を物体(3)に投影し、また位置(C)(図3)におけるカメラ(2)は物体の画像を取り込む。図3において、x軸は位置Sにおけるプロジェクタと位置(C)におけるカメラとの間の線であり、z軸はカメラの軸であり、y軸はこれらの軸に直交している。
【0044】
プロジェクタが示されている図2では、照明されている物体が(10)に示されており。レーザ(5)からの光がレンズ(6)および(7)を通り、レーザストライプ(8)をY−Zステージ(9)上のミラーに投影している。その結果、一つはレーザから直接到達し、他方はミラーを介して到達する二本の光線が画像上に投影され、干渉縞を形成する。ステージ(9)を動かしステップ単位にミラーを動かすことにより、二つの光線間の位相差を調整し、これによって干渉縞を調整することができる。
【0045】
表面座標の計算がかなり単純になれば、座標系の選択は注意深く考えられなければならない。基本的なシステムの構成と向きは、上記したように図1および3に示されている。図8は、選択したX軸とZ軸、および関連する系の次元と角を示している。Y軸は図の平面に垂直である。
【0046】
図8は、タイル平面Y=0の構成を示す、しかし、縞は“垂直”、すなわちYにおいて一定であるので、角θA、θS、θXは点PのY座標から独立している。したがって、点PのXおよびZ座標は、Y座標の計算により、任意のYに関して計算することができる。
【0047】
上記の図において、SS’は光源(プロジェクタ)の軸であり、ゼロ次数の縞に対応する方向である。投影された全ての縞はプロジェクタ軸の一つの側にあるので、θSの全ての値は正である。角θSは、以下に述べるように絶対位相Фおよび縞の角距離から求められる。
【0048】
CC’はカメラ軸であり、座標系のz軸に一致する。θXは正および負両方の値を取る。
カメラおよびプロジェクタは、X軸に沿った距離Dの分だけ離れている。
座標の計算
このセクションは、点Pのx、y、z座標を、入力されたパラメータから求める方法について述べる。
【0049】
図8において、カメラ(C)から見ると、
x=z.Tan(θX)
数式7
プロジェクタ(S)から見ると、
【0050】
【数12】
数式8
【0051】
【数13】
数式9
xおよびz座標は、タイル入力パラメータでは以下のように表すことができる。
【0052】
【数14】
数式10
ここで、xは数式7から直接求めることができる。
【0053】
x=z.Tan(θx)
最後に、y座標は次のように表すことができる。
【0054】
【数15】
ここで、θYは、図9に示すように、x−z平面とカメラから物体の点Pまでのタイル方向との間の角として定義される、角θXおよびθYは、以下のように画像のピクセル位置に関連している。
【0055】
画像平面の中央ピクセルがz軸に沿ったビューに対応するようにシステムが位置決めされたと仮定する。カメラが各行においてNx個のピクセルを持つ場合、列は以下のように番号がつけられる。
【0056】
【数16】
ピクセルの位置は、z軸からの角θXにおけるビューに対応する。ピクセル距離がwであり、カメラレンズが焦点距離vを持つ場合、ピクセル番号nxに対するθXは、
【0057】
【数17】
同様に、カメラが各列においてNy個のピクセルを持つ場合、ピクセル番号nyに対するθYは、
【0058】
【数18】
点Pに関する角θSは、絶対縞位相Фに関係する。数式1から
Ф=C*Sin(θS)
ここでCは定数である。
【0059】
Cの値は、縞の次数(またこれによりФ)を一つの角θSに対して計測する別の較正手順から得られる。勿論、Cの実際の値は縞がスキャンされるにつれて変化する。したがって、フレームの集合の中間位置に設定された縞を使ってシステムを較正することが重要である。この中間位置は、位相計算から得られる位置と同じである。
【0060】
角θAはシステムを位置決めすることによって直接設定される。Sにおけるプロジェクタが、カメラの向きによって定義されたようにx軸上に存在することが、有効な解析にとって不可欠である。
【0061】
図10においては、10aが正面図、10bが側面図、10cが平面図であり、縞プロジェクタは、レーザダイオード(26)、ミラー(25)、および圧電作動装置で構成され、レーザキャビティとミラーの表面との間の距離dを調整する。投影された縞パターンは、レーザキャビティSによって形成されたYoungの干渉縞としてみなすことができ、その画像はミラーS’に存在する。ここで、PSおよびPS’間の光路の差は2dsinδpに等しく、δpはミラーの表面の方向およびPとプロジェクタを結ぶ線の方向との間の角である。SからS’まで出力された光の間のPにおける”大域的“位相差Фpは、以下の数式によって与えられる。
【0062】
【数19】
数式11
ここで、λは投影された光の波長である。Pにおける局所位相、すなわち0から2πの不確かな範囲内のФpの値はこれによりφpとして表される。
【0063】
Фp=φp+2Nπ
数式12
ここで、NはPにおける干渉縞の次数である。φpを計測することのできない暫時位相計測干渉法(TPMI)(19)と呼ばれる確立された手法が数多く存在し、これらの干渉法では、制御された位相ステップであるαpjがФpに加えられ、各ステップに対して点Pにおける対応するインターフェログラムの強度が以下のように記録される。
【0064】
lpj=Ap+Bpcos(Фp+αpj)
数式13
ここで、添え字jは位相ステップのシーケンス番号を表す。Carreの手法は、特に、導入されるべき4つの等しいステップαp、1、2、3、4=−3/2αp,−1/2αp、1/2αp、3/2αpを必要とする。これにより下記の数式が与えられる。
【0065】
【数20】
数式14
従来のシステムでは、位相ステップαpは、投影された光全体において同じである。しかし、アルゴリズム自体は統一されている必要があまりないので、ここで提案されたシステムのPにおける縞の干渉次数を識別するためにαpが使われる。ここで、位相ステッピングはΔdの等しいステップ単位でdを変更することによって実現され、数式11にしたがって下記の位相ステップを生成する。
【0066】
【数21】
数式15
数式15と14を組み合わせると、以下のようにФPの概算が与えられ、この概算から縞の干渉次数を以下のように計算することができる。
【0067】
【数22】
数式16
ここでtrunc()は整数に切り捨てることを表す。数式15および14を12に導入すると、Фpのより正確な計測を得ることができる。3D空間の正確な位置Pを絶対グローバルФpを使って不確かさなしに識別することができる。さらに、数式14から|αp|<πであることが分かり、sinδP<sinθであることが分かる。これらの条件を数式15に応用する。各ステップにおけるミラーの変位は次のように推定される。
【0068】
【数23】
数式16
これは、圧電作動装置の範囲内に収まる。
【0069】
図11において、縞プロジェクタのための別の装置が示されており、図11aはその正面図、11bは側面図、11cは平面図である。ここで、二つのレーザダイオード(15)はミラー(16)上に並べられ恒久的に固定されている。これら二つの商d/λは、二つのレーザに関連する大域的位相が物体の表面上の同じ点Pでわずかに異なるように設計され、以下のように表すことができる。
【0070】
【数24】
数式17
ここで、一度にレーザAとBを切り替えることによってより成熟した3フレーム手法を使用し、関連する局所位相ФpAおよびФpBを得ることができる。これは以下のように表される。
【0071】
【数25】
数式18
ここで、制御された一連の小さい傾きをミラーレーザに導入することによって位相ステッピングを実現することができ、縞次数を以下のように計算することができる。
【0072】
【数26】
数式19
同様に、数式(12)と(9)を(2)に導入することによって、大域的位相Фp(この場合ФpA)のより正確な計測を得ることができる。
【0073】
この別の装置では画像を得る時間が少し長くかかるが(合計6フレーム。前の装置では4フレーム)、データ処理がはるかに早いと言う利点がある。これは、Carreの手法と比較すると、位相検索アルゴリズムがより簡単ではるかに成熟しているからであり、これによってより高い精度と恐らくは短い3Dプロファイルサイクルを実現する。さらに、ダイオード−ミラー組立体の全体を大量生産でひとつの基板上に作成することができ、これによってプロジェクタの堅牢さを高め、縞の質を改善して形状と大きさを小さくする。
【0074】
上記した二つの装置の共通の特徴は、レーザダイオード、特にパワーの高いレーザダイオードが光源として使用されることである。これらのダイオードは、パルスモードにおいて最大100Wの光パワーを生成することができ、計測対象の物体の距離および大きさによって、非常に広い範囲にわたって簡単に調整することができる。このようなダイオードはキャビティ配列で構成されるので、通常干渉計に応用することができず、空間コヒーレンスが劣化する。しかし、キャビティ配列がミラーの表面に平行に位置決めされる限り、これは提案されたプロジェクタにおいて問題にはならない。典型的なスペクトルラインの幅である12nmでは、これらのダイオードも非常に短い一時的なコヒーレンスの長さ(約70nm)を持っている。提案されたプロジェクタの固有の構成では、このように短いコヒーレンスの長さは十分に質の高い干渉縞を生成するのに十分なだけでなく、システムに対する特定の抑制(suppression)という利点をさらにもたらす。両方の構成の構造の単純さは、システムの計測精度にとって重要である縞パターンの安定化の改善にも役立つ。
【0075】
さらに、レーザダイオードは非常に効率的な素子である。図1に示されたように、レーザダイオードのライン幅に適するバンドパス干渉フィルタがカメラの絞りに設けられる。カメラのシャッタに対して同期されたレーザパルスと共に、フィッタがシステムの周囲光に対する耐性を少なくとも600倍高めることができる。したがって100Wパルスダイオードは、通常の白色光スライドプロジェクタの60KW電球に相当する。プロジェクタの出力光は非常によく拡散するので、このようなパワーレベルは、設計された動作距離における人間の視覚の安全な限界内に依然として収まる。波長およびパルス幅を注意深く選択することによって、安全性をさらに改善することができる。
【図面の簡単な説明】
【0076】
【図1】図1は、プロファイリングシステムを示す図である。
【図2】図2は、プロジェクタの構成を示す図である。
【図3】図3は、プロファイリングシステムの平面図である。
【図4】図4は、上記で言及された図である。
【図5】図5は、上記で言及された図である。
【図6】図6は、上記で言及された図である。
【図7】図7は、上記で言及された図である。
【図8】図8は、平面Y=0での位相の構成を表す図である。
【図9】図9は、座標の計算を示す図である。
【図10】図10は、装置の一例を示す図である。
【図11】図11は、装置の他の例を示す図である。
Claims (13)
- 物体の表面における点の集合に関する3次元の表面座標を計算する方法であって、
縞の集合を使って前記物体を照明し、
前記縞を調整し、
異なる縞の位相設定でカメラを使って前記表面の複数の画像を取り込み、
前記画像を処理して、プロジェクタによって照明されると共にカメラに写った前記表面の部分の絶対縞位相写像を生成し、
前記縞位相写像を処理して、前記物体の表面上における点の座標の集合を与える、
ことを特徴とする方法。 - 前記縞は干渉縞である、
ことを特徴とする請求項1に記載の方法。 - 前記縞は、略同等の角あるいは他の既知の距離を持ち、位相および空間周波数において調整可能であり、
前記カメラは、異なる縞位相設定で前記物体のいくつかの画像を取り込む、
ことを特徴とする請求項1または2に記載の方法。 - 前記縞は、一つはレーザ光源から直接あるいはレンズまたは反射系を介してその像から到達し、他方は前記光源の像を生成する干渉計ミラーからの反射または屈折を介して到達する二つの波の干渉によって生成される、
ことを特徴とする請求項2または3のいずれか一つに記載の方法。 - 前記縞の調整が、前記投影されたフィールドをスイープし、前記干渉ミラーをステッピングして前記フィールド全体の縞の次数の関数として変化する小さい位相シフトを加えることによって位相ステップ単位で実行される、
ことを特徴とする請求項4に記載の方法。 - 前記位相ステップサイズは前記絶対縞位相に関連する、
ことを特徴とする請求項5に記載の方法。 - 前記ミラーのステッピングにおける前記ステップの距離は、明細書中の数式1に従って計算される、
ことを特徴とする請求項5あるいは6に記載の方法。 - 前記ミラーは5つのステップごとにステッピングされ、6つのフレームを生成する、
ことを特徴とする請求項6あるいは7に記載の方法。 - 少なくとも一つが縞位相値0あるいはπに近くないように前記位相値が異なる4フレームの2つ(またはそれ以上)の集合が使用される、
ことを特徴とする請求項7または8のいずれか一つに記載の方法。 - 少なくとも5フレームの集合が、4フレームの複数の集合として獲得され処理される、
ことを特徴とする請求項8または9のいずれか一つに記載の方法。 - 5フレームの集合が獲得され、1、2、3、4、5と番号が付けられ、フレーム1、2、3、4、および2、3、4、5のフレームからなる4フレームの二つの集合としてみなされる、
ことを特徴とする請求項4から9のいずれか一つに記載の方法。 - 少なくとも6フレームの集合が、4フレームの集合として獲得され処理される、
ことを特徴とする請求項4から9のいずれか一つに記載の方法。 - 少なくとも6フレームの集合が、4フレームの3つの集合として獲得され処理される、
ことを特徴とする請求項4から9のいずれか一つに記載の方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GBGB0119036.2A GB0119036D0 (en) | 2001-08-06 | 2001-08-06 | Three dimensional imaging |
PCT/GB2002/003625 WO2003014663A1 (en) | 2001-08-06 | 2002-08-06 | Three dimensional imaging by projecting interference fringes and evaluating absolute phase mapping |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2004537732A true JP2004537732A (ja) | 2004-12-16 |
Family
ID=9919806
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003519349A Withdrawn JP2004537732A (ja) | 2001-08-06 | 2002-08-06 | 干渉縞を投影し絶対位相写像を評価することによる三次元画像法 |
Country Status (6)
Country | Link |
---|---|
US (1) | US20040190001A1 (ja) |
EP (1) | EP1417451A1 (ja) |
JP (1) | JP2004537732A (ja) |
CN (1) | CN1551974A (ja) |
GB (1) | GB0119036D0 (ja) |
WO (1) | WO2003014663A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014126385A (ja) * | 2012-12-25 | 2014-07-07 | Azbil Corp | 光学的測定装置 |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0220228D0 (en) * | 2002-09-02 | 2002-10-09 | Univ Heriot Watt | Absolute phase measuring sensor |
EP1531318A1 (en) * | 2003-11-17 | 2005-05-18 | University of Liege | Process and apparatus for measuring a 3D shape of an object |
WO2005049840A2 (en) * | 2003-11-17 | 2005-06-02 | Universite De Liege | Process and apparatus for measuring the three-dimensional shape of an object |
CN100442140C (zh) * | 2005-12-01 | 2008-12-10 | 上海交通大学 | 投影仪生成平移面阵条纹实现累积三维成像的方法 |
CN100417231C (zh) * | 2006-05-31 | 2008-09-03 | 北京航空航天大学 | 立体视觉半实物仿真***及方法 |
CN103729092A (zh) * | 2012-10-12 | 2014-04-16 | 原相科技股份有限公司 | 手持式指向装置及其控制方法 |
KR101613829B1 (ko) * | 2014-05-08 | 2016-04-20 | 조춘식 | 미분 모아레를 이용한 3차원 형상 측정방법 및 장치 |
JP6658625B2 (ja) * | 2017-03-08 | 2020-03-04 | オムロン株式会社 | 3次元形状計測装置及び3次元形状計測方法 |
US20190011380A1 (en) * | 2017-07-10 | 2019-01-10 | Brigham Young University | Wave interference systems and methods for measuring objects and waves |
CN108050955B (zh) * | 2017-12-14 | 2019-10-18 | 合肥工业大学 | 基于结构光投影与数字图像相关的高温空气扰动滤除方法 |
CN110631510B (zh) * | 2019-09-12 | 2020-07-31 | 中国科学院西安光学精密机械研究所 | 一种基于迈克尔逊结构的高精度测角装置及测角方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4794550A (en) * | 1986-10-15 | 1988-12-27 | Eastman Kodak Company | Extended-range moire contouring |
US4818108A (en) * | 1987-12-14 | 1989-04-04 | Hughes Optical Products, Inc. | Phase modulated ronchi testing of aspheric surfaces |
US5069548A (en) * | 1990-08-08 | 1991-12-03 | Industrial Technology Institute | Field shift moire system |
US5636025A (en) * | 1992-04-23 | 1997-06-03 | Medar, Inc. | System for optically measuring the surface contour of a part using more fringe techniques |
US5557410A (en) * | 1994-05-26 | 1996-09-17 | Lockheed Missiles & Space Company, Inc. | Method of calibrating a three-dimensional optical measurement system |
DE4432313C2 (de) * | 1994-09-12 | 2002-11-14 | Fraunhofer Ges Forschung | Vorrichtung zur Untersuchung von Oberflächentopographien mittels Streifen-Triangulation |
KR19990029064A (ko) * | 1995-07-18 | 1999-04-15 | 낸시 엘. 후체슨 | 확장된 영상 깊이를 갖는 모아레 간섭 시스템 및 방법 |
PT888522E (pt) * | 1996-03-22 | 2001-10-30 | Univ Loughborough | Processo e aparelho para a medicao da forma de objectos |
US6438272B1 (en) * | 1997-12-31 | 2002-08-20 | The Research Foundation Of State University Of Ny | Method and apparatus for three dimensional surface contouring using a digital video projection system |
US6788210B1 (en) * | 1999-09-16 | 2004-09-07 | The Research Foundation Of State University Of New York | Method and apparatus for three dimensional surface contouring and ranging using a digital video projection system |
-
2001
- 2001-08-06 GB GBGB0119036.2A patent/GB0119036D0/en not_active Ceased
-
2002
- 2002-08-06 EP EP02751389A patent/EP1417451A1/en not_active Withdrawn
- 2002-08-06 JP JP2003519349A patent/JP2004537732A/ja not_active Withdrawn
- 2002-08-06 US US10/486,397 patent/US20040190001A1/en not_active Abandoned
- 2002-08-06 WO PCT/GB2002/003625 patent/WO2003014663A1/en not_active Application Discontinuation
- 2002-08-06 CN CNA028173856A patent/CN1551974A/zh active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014126385A (ja) * | 2012-12-25 | 2014-07-07 | Azbil Corp | 光学的測定装置 |
Also Published As
Publication number | Publication date |
---|---|
GB0119036D0 (en) | 2001-09-26 |
US20040190001A1 (en) | 2004-09-30 |
EP1417451A1 (en) | 2004-05-12 |
WO2003014663A1 (en) | 2003-02-20 |
CN1551974A (zh) | 2004-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20230392920A1 (en) | Multiple channel locating | |
US7403650B2 (en) | System for simultaneous projections of multiple phase-shifted patterns for the three-dimensional inspection of an object | |
US8462357B2 (en) | Device and method for obtaining three-dimensional object surface data | |
JP6161714B2 (ja) | 3次元の物体の直線寸法を制御する方法 | |
JP7413372B2 (ja) | 対向配置チャネルを有する三次元センサ | |
JPH04220509A (ja) | 物体表面を無接触で測定する方法および装置 | |
US9131219B2 (en) | Method and apparatus for triangulation-based 3D optical profilometry | |
TW201825864A (zh) | 用於圖案化半導體特徵之特徵化的掃描白光干涉測量系統 | |
JP2004537732A (ja) | 干渉縞を投影し絶対位相写像を評価することによる三次元画像法 | |
CN105091781B (zh) | 一种单桢干涉条纹图测量光学表面的方法与装置 | |
Zhang | Digital multiple wavelength phase shifting algorithm | |
CN116518869A (zh) | 基于光度立体和双目结构光的金属表面测量方法及*** | |
RU125335U1 (ru) | Устройство контроля линейных размеров трехмерных объектов | |
JP2006308452A (ja) | 3次元形状計測方法および装置 | |
CN114076579A (zh) | 一种基于偏振成像的三维粗糙度检测装置及方法 | |
JP4988577B2 (ja) | 鏡面域を有する基準面を備えた干渉系 | |
KR20190136592A (ko) | 다수 조명을 이용한 3차원 데이터 획득 방법 및 장치 | |
Schaffer et al. | Accurate outdoor 3D shape measurements with structured illumination using narrow spectral filtering | |
CA2402849A1 (en) | System for simultaneous projections of multiple phase-shifted patterns for the three-dimensional inspection of an object | |
JPH11271034A (ja) | 三次元形状測定装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A300 | Application deemed to be withdrawn because no request for examination was validly filed |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20051101 |