JP2003519786A - 表面輪郭測定のための装置および方法 - Google Patents
表面輪郭測定のための装置および方法Info
- Publication number
- JP2003519786A JP2003519786A JP2001552052A JP2001552052A JP2003519786A JP 2003519786 A JP2003519786 A JP 2003519786A JP 2001552052 A JP2001552052 A JP 2001552052A JP 2001552052 A JP2001552052 A JP 2001552052A JP 2003519786 A JP2003519786 A JP 2003519786A
- Authority
- JP
- Japan
- Prior art keywords
- radiation
- pattern
- sources
- spectral
- source
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 186
- 238000005259 measurement Methods 0.000 title claims description 88
- 230000005855 radiation Effects 0.000 claims abstract description 246
- 230000003595 spectral effect Effects 0.000 claims abstract description 140
- 238000009826 distribution Methods 0.000 claims abstract description 43
- 230000004044 response Effects 0.000 claims abstract description 37
- 230000003287 optical effect Effects 0.000 claims description 54
- 230000001427 coherent effect Effects 0.000 claims description 45
- 230000008859 change Effects 0.000 claims description 40
- 230000000694 effects Effects 0.000 claims description 37
- 230000033001 locomotion Effects 0.000 claims description 37
- 238000003384 imaging method Methods 0.000 claims description 22
- 238000004891 communication Methods 0.000 claims description 17
- 230000000737 periodic effect Effects 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 9
- 230000002194 synthesizing effect Effects 0.000 claims description 8
- 238000001228 spectrum Methods 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 4
- 230000000116 mitigating effect Effects 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 2
- 230000000903 blocking effect Effects 0.000 claims 1
- 238000002955 isolation Methods 0.000 claims 1
- 238000013178 mathematical model Methods 0.000 claims 1
- 230000010363 phase shift Effects 0.000 description 38
- 230000000750 progressive effect Effects 0.000 description 36
- 238000013459 approach Methods 0.000 description 30
- 230000015572 biosynthetic process Effects 0.000 description 26
- 238000003786 synthesis reaction Methods 0.000 description 25
- 230000006870 function Effects 0.000 description 24
- 230000008901 benefit Effects 0.000 description 16
- 230000014509 gene expression Effects 0.000 description 16
- 238000010586 diagram Methods 0.000 description 15
- 238000013519 translation Methods 0.000 description 15
- 239000000835 fiber Substances 0.000 description 14
- 238000005286 illumination Methods 0.000 description 13
- 230000011218 segmentation Effects 0.000 description 11
- 230000008569 process Effects 0.000 description 8
- 230000010355 oscillation Effects 0.000 description 7
- 238000000638 solvent extraction Methods 0.000 description 7
- 238000012935 Averaging Methods 0.000 description 6
- 230000007423 decrease Effects 0.000 description 6
- 238000003491 array Methods 0.000 description 5
- 238000006073 displacement reaction Methods 0.000 description 5
- 239000004973 liquid crystal related substance Substances 0.000 description 5
- 239000000203 mixture Substances 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 230000001419 dependent effect Effects 0.000 description 4
- 238000013461 design Methods 0.000 description 4
- 230000001678 irradiating effect Effects 0.000 description 4
- 230000035945 sensitivity Effects 0.000 description 4
- 230000004075 alteration Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000010287 polarization Effects 0.000 description 3
- 238000004441 surface measurement Methods 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 239000011521 glass Substances 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000006335 response to radiation Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- OVSKIKFHRZPJSS-UHFFFAOYSA-N 2,4-D Chemical compound OC(=O)COC1=CC=C(Cl)C=C1Cl OVSKIKFHRZPJSS-UHFFFAOYSA-N 0.000 description 1
- 241000231739 Rutilus rutilus Species 0.000 description 1
- 102100027340 Slit homolog 2 protein Human genes 0.000 description 1
- 101710133576 Slit homolog 2 protein Proteins 0.000 description 1
- RTAQQCXQSZGOHL-UHFFFAOYSA-N Titanium Chemical compound [Ti] RTAQQCXQSZGOHL-UHFFFAOYSA-N 0.000 description 1
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000004377 microelectronic Methods 0.000 description 1
- 238000000386 microscopy Methods 0.000 description 1
- 230000009022 nonlinear effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 238000009304 pastoral farming Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 229910052594 sapphire Inorganic materials 0.000 description 1
- 239000010980 sapphire Substances 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 239000000377 silicon dioxide Substances 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000009987 spinning Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 238000001308 synthesis method Methods 0.000 description 1
- 229910052719 titanium Inorganic materials 0.000 description 1
- 239000010936 titanium Substances 0.000 description 1
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/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/2509—Color coding
-
- 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
- G01B11/2518—Projection by scanning of the object
- G01B11/2527—Projection by scanning of the object with phase change by in-plane movement of the patern
-
- 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/2531—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 using several gratings, projected with variable angle of incidence on the object, and one detection device
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
Description
の一部継続出願である、1999年2月2日に出願された米国特許出願第09/
241,354号の一部継続出願であり、そして1998年6月4日に出願され
た米国仮特許出願第60/087,960号の優先権を主張する。
628−95−L−002によって援助された。政府は本発明に特定の権利を有
し得る。
る。
が、ほとんどの製造および多くのサブアセンブリから構成される複雑なオブジェ
クトのアセンブリを実行する、今日の製造環境において非常に重要である。構成
要素が適切に嵌合し合うことを保証するには、自動車などの複雑なアセンブリの
各構成要素の形状およびサイズの公差を小さく保持する必要がある。
れて、測定を行う際の時間を節減する。多くの非接触測定方法は、利用可能なマ
シーンビジョンシステムを利用する。表面輪郭情報の測定は、深度情報が損失し
たり、解釈が難しい場合が多いため、マシーンビジョンシステムの特に難しい問
題である。深度情報の損失および利用可能な情報を解釈する際の難しさを補償す
るために、多くのマシーンビジョンシステムは、光を用いてオブジェクトの表面
上にモアレ縞を生成し、輪郭情報を得る。モアレ技術の1つの欠点は、異なるサ
イズのオブジェクトに対するその柔軟性の無さである。異なるサイズのオブジェ
クトは、新しい対応する物理的設備を必要とし得る。この欠点に起因して、大規
模なオブジェクトのモアレ技術を用いることが難しい。モアレ技術の別の欠点は
、この技術を用いて得られる解像度が多くの用途にとって十分に高くない場合が
あることである。
が表面輪郭情報を提供するが、干渉システムは測定されるオブジェクトおよび用
いられる照射源の両方の変化の影響を受けやすい。
する劣化を緩和し、そして高い解像度を達成することによって、これらの問題を
回避する技術である。
のポイントの三次元位置情報を判定する方法に関する。上記方法は、スペクトル
領域の空間分布を有する2つの放射ソースを提供する工程;上記ソースそれぞれ
からの放射で上記表面を照射して、上記表面上の第1の位置に第1の縞パターン
を生成する工程;上記第1の縞パターンを第2の位置に移動させる工程;上記第
1の縞パターンの上記第1の位置および上記第2の位置に応じて、第1のラップ
されたサイクルマップを生成する工程;上記第1の縞パターン中の縞数を推定す
る工程;上記第1の縞パターンを変更して、第1の位置において第2の縞パター
ンを生成する工程;上記第2の縞パターンを第2の位置に移動させる工程;上記
第2の縞パターンの上記第1の位置および上記第2の位置に応じて、第2のラッ
プされたサイクルマップを生成する工程;上記第1の縞パターン中の推定された
縞数に応じて、上記第2の縞パターン中の縞数を推定する工程;および上記第2
の縞パターンおよび上記第2のラップされたサイクルマップ中の推定された縞数
に応じて、上記ポイントの上記表面上での三次元位置情報を判定する工程を含む
。
マップを合成する方法にさらに関する。上記方法は、スペクトル領域の空間分布
を有する2つの放射ソースを提供する工程;上記ソースそれぞれからの放射で上
記表面を照射して、上記表面上の第1の位置に第1の縞パターンを生成する工程
;上記第1の縞パターンを上記表面上の第2の位置に移動させる工程;上記第1
の縞パターンの上記第1の位置および上記第2の位置に応答して、第1のラップ
されたサイクルマップを生成する工程;上記第1の縞パターンを変更して、第1
の位置において第2の縞パターンを生成する工程;上記第2の縞パターンを第2
の位置に移動させる工程;上記第2の縞パターンの上記第1の位置および上記第
2の位置に応答して、第2のラップされたサイクルマップを生成する工程;上記
第1のラップされたサイクルマップから上記第2のラップされたサイクルマップ
を減算する工程;および上記第2のラップされたサイクルマップと上記第1のラ
ップされたサイクルマップと間の差をラップして、上記所定の間隔の縞パターン
に対応するラップされたサイクルマップを生成する工程を含む。
のポイントの三次元位置情報を判定する方法にさらに関する。上記方法は、第1
の周波数を有し、そして互いにコヒーレントな第1の放射ビームおよび第2の放
射ビームを生成する工程;上記表面を上記第2の放射ビームで照射する工程;上
記第1の放射ビームと、上記第2の放射ビームからの放射とが上記表面によって
散乱するのに応じて、第1の位置において第1の干渉パターンを生成する工程;
上記第1の干渉パターンを第2の位置に移動させる工程;上記第1の干渉パター
ンの上記第1の位置および上記第2の位置に応じて、第1のラップされたサイク
ルマップを生成する工程;上記第1の干渉パターン中の強度サイクルを推定する
工程;上記第1の干渉パターンを変更して、第2の干渉パターンを第1の位置に
おいて生成する工程;上記第2の干渉パターンを第2の位置に移動させる工程;
上記第2の干渉パターンの上記第1の位置および上記第2の位置に応じて、第2
のラップされたサイクルマップを生成する工程;上記第1の干渉パターン中の推
定された強度サイクルに応答して、上記第2の干渉パターン中の強度サイクルを
推定する工程;および上記第2の干渉パターンおよび上記第2のラップされたサ
イクルマップ中の推定された強度サイクルに応じて、三次元位置情報を計算する
工程を含む。
のポイントの三次元位置情報を判定する方法にさらに関する。上記方法は、互い
にコヒーレントな第1の放射ビームおよび第2の放射ビームを生成する工程;上
記表面を上記第2の放射ビームで照射する工程;上記第1の放射ビームと、上記
第2の放射ビームからの放射とが上記表面によって散乱するのに応じて、第1の
干渉パターンを生成する工程;上記第1の放射ビームおよび第2の放射ビームの
うち1つを、上記第1の放射ビームおよび第2の放射ビームのもう一方に対して
位相シフトさせて、第1の位相シフトした干渉パターンを生成する工程;上記第
1の干渉パターンおよび上記第1の位相シフトした干渉パターンに応じて第1の
ラップされたサイクルマップを生成する工程;上記第1の干渉パターンを変更し
て、第2の干渉パターンを生成する工程;上記第1の放射ビームおよび第2の放
射ビームのうち1つを、上記第1の放射ビームおよび第2の放射ビームのもう一
方に対して位相シフトさせて、第2の位相シフトした干渉パターンを生成する工
程;上記第2の干渉パターンおよび上記第2の位相シフトした干渉パターンに応
じて、第2のラップされたサイクルマップを生成する工程;上記第2のラップさ
れたサイクルマップを上記第1のラップされたサイクルマップから減算する工程
;および上記第2のラップされたサイクルマップと、上記第1のラップされたサ
イクルマップとの間の差をラップして、ラップされたサイクルマップを生成する
工程を含む。
置は、スペクトル分布を有する2つの放射ソース;上記2つのソースと光学的に
連絡したコリメータであって、2本の実質的に平行な広域幅放射ビームを生成す
る、コリメータ;上記コリメータと光学的に連絡した回折格子;および上記回折
格子と光学的に連絡したレンズであって、スペクトル領域の空間分布を有する2
つの放射イメージを生成する、レンズを有する。
法は、ある距離の分だけ隔離された2つの放射ソースを提供する工程であって、
上記放射ソースはスペクトル分布を有する放射を生成する、工程;上記放射を平
行にして、2本の実質的に平行な放射ビームを生成する工程;上記平行な放射ビ
ームのスペクトル成分を描写する工程;およびスペクトル成分の空間分布を有す
る2つの放射イメージを生成する工程を含む。
のポイントの三次元位置情報を判定する方法にさらに関する。上記方法は、第1
の周期的パターンを投射して、上記表面上の第1の位置に第1の投射パターンを
生成する工程;上記第1の投射パターンを第2の位置に移動させる工程;上記第
1の投射パターンの上記第1の位置および上記第2の位置に応答して、第1のラ
ップされたサイクルマップを生成する工程;上記第1の投射パターン中のサイク
ル数を推定する工程;上記第1の投射パターンを変更して、第1の位置において
第2の投射パターンを生成する工程;上記第2の投射パターンを第2の位置に移
動させる工程;上記第2の投射パターンの上記第1の位置および上記第2の位置
に応じて、第2のラップされたサイクルマップを生成する工程;上記第1の投射
パターン中の推定されたサイクル数に応じて、上記第2の投射パターン中のサイ
クル数を推定する工程;および上記第2の投射パターン中の推定されたサイクル
数および上記第2のラップされたサイクルマップに応じて、上記表面を判定する
工程を含む。
されたサイクルマップを合成する方法にさらに関する。上記方法は、第1の周期
的パターンを投射して、上記表面上の第1の位置に第1の投射パターンを生成す
る工程;上記第1の投射パターンを上記表面上の第2の位置に移動させる工程;
上記第1の投射パターンの上記第1の位置および上記第2の位置に応じて、第1
のラップされたサイクルマップを生成する工程;上記第1の投射パターンを変更
して、第1の位置において第2の投射パターンを生成する工程;上記第2の投射
パターンを第2の位置に移動させる工程;上記第2の投射パターンの上記第1の
位置および上記第2の位置に応じて、第2のラップされたサイクルマップを生成
する工程;上記第1のラップされたサイクルマップから上記第2のラップされた
サイクルマップを減算する工程;および上記第2のラップされたサイクルマップ
と上記第1のラップされたサイクルマップと間の差をラップして、上記所定の間
隔の投射された周期的パターンに対応するラップされたサイクルマップを生成す
る工程を含む。
軽減する方法にさらに関する。上記方法は、コヒーレントな縞パターンを生成す
る工程;上記縞パターンが上記オブジェクトの表面を実質的にかするように、上
記コヒーレントな縞パターンを上記オブジェクトの表面上に光路に沿って投射す
る工程;および上記オブジェクトの表面のイメージ中に上記縞パターンおよび上
記スペックルを検出する工程であって、上記オブジェクトの表面に対する法線が
上記光路に実質的に直交する、工程を含む。
軽減する方法にさらに関する。上記方法は、ある距離の分だけ隔離された2つの
放射ソースからコヒーレントな縞パターンを生成する工程;上記コヒーレントな
縞パターンを上記オブジェクトの表面上に投射する工程;上記オブジェクトの表
面のイメージ中に上記縞パターンおよび上記スペックルを検出する工程;上記縞
パターンが実質的に静止したままで、そして上記スペックルが変化するように上
記2つのソースを平行移動させる工程;上記オブジェクトの表面のイメージ中に
上記新しい縞パターンおよび上記変化したスペックルを検出する工程;ならびに
上記検出された縞パターンおよび上記スペックルの変化に応じて、実質的にスペ
ックルの無い上記縞パターンを判定する工程を含む。
軽減する方法にさらに関する。上記方法は、ある距離の分だけ隔離された2つの
放射ソースからコヒーレントな縞パターンを生成する工程;上記コヒーレントな
縞パターンを上記オブジェクトの表面上に投射する工程;上記オブジェクトの表
面のイメージ中に上記縞パターンおよびスペックルを検出する工程;上記縞パタ
ーンが上記オブジェクトの表面に対して実質的に静止したままで、そして上記ス
ペックルが変化するように、上記オブジェクトを上記縞パターンの等位相平面に
平行な経路に沿って上記2つの放射ソースに対して平行移動させる工程;上記オ
ブジェクトの表面のイメージ中に上記新しい縞パターンおよび上記変化したスペ
ックルを検出する工程;ならびに上記検出された縞パターンおよび上記スペック
ルの変化に応じて、上記実質的にスペックルの無い縞パターンを判定する工程を
含む。
軽減する方法にさらに関する。上記方法は、コヒーレントな縞パターンを上記オ
ブジェクトの表面上に投射する工程;透過作用を有するレンズを提供する工程で
あって、上記レンズの透過は上記レンズ端部において徐々に低下する、工程;お
よびオブジェクトの表面のイメージ中に上記縞パターンを検出する工程であって
、上記透過作用は上記スペックルが上記測定に与える影響を実質的に減少させる
、工程を含む。
軽減する方法にさらに関する。上記方法は、第1の周波数を有し、そしてある距
離の分だけ隔離された2つの放射ソースから第1の縞パターンを生成する工程;
上記第1の周波数を第2の周波数に変化させることによって、上記第1の縞パタ
ーンを変化させて、第2の縞パターンを生成する工程;および上記第1の周波数
と上記第2の周波数との差に応答して上記距離を変化させる工程であって、上記
距離と、上記第1の周波数と上記第2の周波数との差との比は実質的に一定であ
る、工程を含む。
法は、ある距離の分だけ隔離された2つの放射ソースを提供する工程であって、
上記放射ソースは、スペクトル分布を有し、そしてもう一方のソースに対してコ
ヒーレントである、工程;上記ソースそれぞれからの放射で上記オブジェクトの
表面上のポイントを照射する工程;上記ソースのうちの一つをもう一方のソース
に対して移動させる工程;および上記オブジェクトの表面上のポイントによって
散乱された放射を検出する工程を含む。本発明の別の実施形態において、2つの
放射ソースを提供する工程は、スペクトル幅を有する初期放射ビームを提供する
工程;上記初期放射ビームから、第1の放射ビームを第1のビーム角度で、第2
の放射ビームを第2のビーム角度で生成する工程;ならびに上記第1の放射ビー
ムおよび上記第2の放射ビームをイメージングして上記2つの放射ソースを形成
する工程を含む。
置は、ある距離の分だけ隔離された2つの放射ソースであって、各ソースは、ス
ペクトル分布を有し、そしてもう一方のソースに対してコヒーレントである、放
射ソース;上記ソースそれぞれをもう一方のソースに対して移動させる制御シス
テム;および上記オブジェクトの表面上のポイントから散乱された放射を受け取
るように配置された検出器を含む。本発明の別の実施形態において、2つの放射
ソースは、スペクトル幅を有する放射ビームの初期ソース;上記放射ビームの初
期ソースと光学的に連絡したビーム分離器であって、第1の光学ビームおよび第
2の光学ビームを生成する、ビーム分離器;および上記ビーム分離器に光学的に
連絡したイメージングシステムであって、上記イメージングシステムは、上記2
つの放射ソースを生成し、そしてそれぞれの放射ソースは上記第1の光学ビーム
および上記第2の光学ビームに対応する、イメージングシステムを含む。
のポイントの三次元位置情報を判定する方法にさらに関する。上記方法は、2つ
の放射ソースを提供する工程であって、上記放射ソースは、スペクトル分布を有
し、そして上記2つの放射ソースのもう一方に対してコヒーレントである、工程
;検出器を上記表面上のポイントに提供する工程;上記オブジェクトの上記表面
のポイントを上記ソースのそれぞれからの上記放射で照射する工程;上記ソース
のそれぞれをもう一方に対して移動させる工程;上記オブジェクトの表面上のポ
イントに上記放射を検出する工程;ならびに上記ソースの移動および上記オブジ
ェクトの表面上のポイントに検出された放射に対して位置情報を計算する工程を
含む。
のポイントの三次元位置情報を判定する装置にさらに関する。上記装置は、2つ
の放射ソースであって、スペクトル分布を有し、そして上記2つの放射ソースの
もう一方に対してコヒーレントである、放射ソース;上記ソースのそれぞれをも
う一方に対して移動させる制御システム;上記オブジェクトの表面上のポイント
に配置された検出器であって、上記オブジェクトの表面上のポイントを照射する
放射を受け取る、検出器;上記検出器から信号を受け取るプロセッサ;ならびに
上記ソースの移動および上記オブジェクトの表面上のポイントで受け取られた放
射に対して上記オブジェクトの表面上のポイントの位置情報を計算するプロセッ
サを含む。
利点は、添付の図面に関する以下の説明を参照することによってよりよく理解さ
れ得る。
これらの用語は、非常に局所化された放射ソースを含む、任意の放射ソースを指
すよう意図される。
定距離Dの分だけ隔離されており、そしてそれぞれ(x1、y1、z1)および
(x2、y2、z2)の空間座標をそれぞれ有する。それぞれのソース、P1お
よびP2からの放射は、もう一方のソースからの放射に対してコヒーレントであ
る。各ソース、P1およびP2は、各発散放射ビーム12および14をオブジェ
クト10の表面上のポイントP0に方向付ける。それぞれの放射ソースP1およ
びP2から表面上のポイントP0までの距離は、R1およびR2によってそれぞ
れ示される。Ψは起点からポイントP0まで延びる線と、ソース、P1とP2と
の間を延びる線との間の角度であり、θsはz軸とソース、P1とP2との間を
延びる線との間の角度であり、αはP0から見た場合の2つのソースポイントか
ら延びた辺で囲まれた角度の半分の角度である。各ビーム12、14はもう一方
のビーム14、12と同じ方向に実質的に偏光され、そして、オブジェクト10
上の異なる領域を同時に照射するように独立して走査可能であり得る。あるいは
、オブジェクト10全体が同時に照射され得る。
一実施形態において、光検出器22は1アレイの光検出器の素子を含んで、測定
されるオブジェクト10の二次元画像を提供する。さらなる実施形態において、
1アレイの光検出器の素子は電荷結合素子(CCD)である。検出器22は、1
つ以上の個々の信号を含み、各信号は検出器22の光検出器の素子のうちの対応
する1つに関連付けられた出力信号26を提供する。
イントP0と光検出器22との間に配置されて、検出器22上にポイントP0を
含むオブジェクトの照射された部分をイメージングする。オブジェクトの表面が
粗いため、そして照射放射がコヒーレントであるため、焦点調節されたイメージ
はスペックルを有する。光検出器22からの出力信号26は処理装置28への入
力信号である。
置される。偏光子30の配置方向は、散乱光20の主な偏光成分との合致が最大
化される方向に向けられる。これにより、スペックルのコントラストまたは縞パ
ターンのコントラストが向上する。この構成により、オブジェクト10の表面か
ら散乱された光に関連付けられた信号対雑音比が最大化される。
イ22の光検出器の素子のそれぞれに関連付けられた検出器の出力信号26を受
けて動作する。別の実施形態において、プロセッサ28は、複数の個々のプロセ
ッサを有するマルチプロセッサであり、そして各光検出器の素子は入力信号をプ
ロセッサのそれぞれ一つに提供する。検出器22がCCDアレイである、また別
の実施形態において、複数のCCD素子は入力信号をマルチプロセッサの各プロ
セッサに提供する。マルチプロセッサの構成を用いると、複数の個々の光素子か
らの信号の計算が実質的に同時に行われ、これにより、信号処理速度が向上する
。
ブジェクト10の表面上のポイントP0で測定されるもう一方のソースからの放
射の位相に対して変更するように、放射ソースP1およびP2の動作を制御する
。プロセッサ28は、信号線またはバス34を介して制御装置32と通信し得る
。例えば、特定の用途において、プロセッサ28が検出器22からの信号を処理
する場合には、オブジェクト10の表面上を介するソースP1およびP2の走査
に対応して、またはソースからの放射が掃引される周期に対応して、特定の回数
行うことが望ましい場合がある。このような走査および周波数掃引動作が制御装
置32によって制御されるため、制御装置32とプロセッサ28との間の通信が
これらの環境において望ましい。制御装置32およびプロセッサ28は物理的に
別々の装置であっても、単一の処理システムによって実現されてもよいことが理
解される。
チューナブルレーザ40から放出される放射から形成される。チューナブルレー
ザ40によって放出された放射ビーム44はビームスプリッタ48によって分割
される。放射ビーム50はビームスプリッタ48によって反射された後、レンズ
52によって発散するようになる。次いで、発散ビームは、移動可能な照準ミラ
ー54によって反射される。照準ミラー54によって反射される放射ビームは、
コヒーレントな放射ソースのうちの一つ、P1を提供する。同様に、ビームスプ
リッタ48を透過した放射ビーム46は、レンズ58によって発散するようにな
り、発散ビームを第2の移動可能な照準ミラー60に方向付ける。ミラー60に
よって反射された放射ビームは、第2の放射ソースP2を提供する。照準ミラー
54および62は、回転可能であり、オブジェクト10の表面を選択的に照射し
得る。照準ミラー54および62はさらに、移動可能であり、ソースP1および
P2の位置を変更し得る。
44を提供する回転可能なレーザソース40を含むように図示される。放射ビー
ム44はレンズ62を透過し、これによりビームが発散して、発散ビーム64を
提供する。次いで、発散ビーム64はビームスプリッタ48によって反射されて
、第1のビーム66を提供する。第2のビーム68は、図示するように、ビーム
スプリッタ48を透過する。移動可能な照準ミラー54および60はそれぞれ、
ビーム66および68を反射して、ソースP1およびP2を提供する。
ファイバスプリッタ56を用いて、回転可能なレーザ40から放出された放射か
ら分割される。ファイバは、端部においてビーム形成要素を有して、2本のビー
ムの発散角度を制御または設定し得る。一実施形態において、ビーム形成要素は
レンズであり得る。ソースP1およびP2は、あるいは、周波数が固定された1
対の回転可能なレーザから形成され得る。放射ソースの他の適切な実施形態は、
マイクロ波および音波などの制御可能な位相を有する波を生成する任意のソース
を含む。
放射P1およびP2の各ソースを取り付けることによって、ソースは互いから一
定の距離Dで維持される。別の実施形態において、放射P1およびP2のソース
は一定の距離で保持されないが、その代わりにこれらのソース間の距離Dは非常
に正確であることが公知である。
的な1つのバー70が図3に示される。バー70の両端部にはソケット74が提
供される。玉継手76は、示されているように、ソケット74の各々の内部で旋
回するように(pivotally)配置される。玉継手76の各々は(図2b
に示される)その中に配置される光ファイバースプリッタ56からのファイバー
の端部および発散放射が通過するアパーチャ80を有する。ファイバーはその端
部においてビームを形成する素子を有し得、2つのビームの発散角を制御または
設定し、1実施形態において、このビーム形成素子はレンズである。動作中、玉
継手76は、それぞれのソケット74の中で、矢印78によって示されるように
旋回可能であり、(図1に示される)制御ユニット32によって制御され得る。
この構成を用いて、ファイバーの端部におけるソースP1およびP2によって提
供される発散ビーム12および14は、所望の方向に向けられ得、一定の隔離距
離(separation distance)Dが維持される一方で、処理さ
れるべき点P0を含む物体10の全体または部分を照射する。
る。点P0のx座標およびy座標は、通常、検出器22および物体10の幾何学
的形状から直接的に決定されるが、焦点調節素子24をこれらの検出器22と物
体10との間に挿入して任意の倍率を考慮している場合には、深さ座標zを直接
求めることはできない。但し、ここでz軸は結像系の光軸と並行であると定義さ
れる。しかしながら、深さ座標zは、放射ソースP1から物体10の表面上の点
P0および放射ソースP2までの光路差 s=R2−R1+S0 (1) をまず考察することによって測定され得る。量S0は、ビームが点P1およびP 2 に到達する前に生じ得る、ビームにおける任意の光路長の差を考慮に入れるた
めに含まれる。
変更することによって点P0で測定される1方のソースからの放射の位相が、他
方のソースに対して変化する。この位相変化によって、点P0における放射の強
度が変調する。強度の変化の1サイクルを完了するために必要とされる周波数Δ
νの変化は、以下の式によって与えられる。 Δν=c/s (2) ここで、cは光の速度である。従って、強度の1振動を引き起こすために必要と
されるレーザの周波数Δνの変化を測定することによって、光路差sが決定され
得る。zの測定は、その後、後述されるように、xおよびyの各値に対してsの
値を決定することに基づく。
よって得られる。実際には、周波数の変化全体Bによって誘導される(必ずしも
自然数ではない)振動サイクル数Nに関して測定することが好都合である。
が得られる。 s=(c/B)・N (4) 長さBの周波数走査により誘導された振動サイクル数Nは、周波数の走査中、
測定点P0にわたって移動する干渉縞の数に対応することに留意されたい。干渉
縞は、ソースP1およびP2から発せられる放射の干渉によって生成される。こ
れらの振動は、結像レンズ24によって干渉縞が解消されてもされなくても生じ
る。
ある場合、sにおける不確定度Δsは、測定されたサイクルN数と等しいファク
タだけ減算されることを示す。当業者に公知の種々のレベルの分解能ΔNに対す
るNを決定する複数の方法が存在する。約1の振動サイクルカウント(ΔN=1
)の分解能を得る方法の例は、データシーケンスに高速フーリエ変換(FFT)
を行なうか、または高域フィルタリングされた信号のゼロ交差をカウントするこ
とである。1未満の振動サイクルカウント(ΔN<1)の分解能は、例えば、離
散フーリエ変換(DFT)の大きさが最大化されるところでDFTの引数(ar
gument)を見出すことによってか、または周波数走査の終了時において振
動サイクルの位相を検査することによって改善することができる。位相の精確な
検査に関する当業者に公知の1つの技術は、ビーム経路の1つのレッグ、すなわ
ちビームスプリッタまたは光ファイバースプリッタと、図2、図2(a)および
図2(b)におけるソースP1またはP2のうちの1つとの間に位相変調器を挿
入することである。
よって生じる位相シフトに対応する信号強度である場合、振動サイクルの位相φ
は、以下の数式によって与えられる。 φ=tan−1((I1−I3)/(2I2−I1−I3)) (
6) 位相における不確定度Δφは、2πで除算することによってサイクルカウント
における不確定度ΔNに変換される。
、ΔN=1サイクルの不確定度では、Δs=20μmの不確定度が提供される。
横方向の分解能に対するsの広がりが、その量よりも小さい場合、ΔN=0.1
サイクルの不確定度は、sにおける不確定度をΔs=2.0μmまで向上する。
物体の表面上の横方向の分解能に対するsにおける広がりがΔsよりも大きい場
合、sの測定における分解能は改善されているので、この場合もやはり、その横
方向の分解能に対するsの平均値または代表値の推定を改善することができる。
れていると仮定する。その後、式(7)は、(x1、y1、z1)に関して、以
下のようになる。
9)に正の根および負の根が存在するために、アンビギティの範囲内であること
が判定される。このアンビギティを回避する1つの方法は、物体10を照射する
ことにより、s=0線(s0=0の場合である図1において16で示される)が
結像されるべき物体の領域を2等分しないようにすることである。s=0線を移
動する1つの方法は式(1)におけるs0を変更することである。
、Δzはsの値における不確定度Δsによって導入されるzにおける不確定度で
ある。この比は、任意の実用的な感度の範囲が十分でないシステムの0と、理論
的には最大システムの2との間にわたる。2の値を達成するのは非現実的である
。なぜなら、物体10の表面は、2つの点ソースP1とP2との間にある必要が
あり、各ビームによって表面の片面のみが照射されとるからである。比Δs/Δ
zは、zに対するsの偏導関数を得ることによって計算され、この結果より範囲
分解能に関する以下の数式が取得される。
あり、ψは原点から点P0に延びる線と、点P1から点P2に延びる線との間の
角である。良好な範囲の分解能を提供する有用な構成は、ψ=90°を設定する
ことであり、この場合Δzの数式は以下のように簡略化される。
tanα=D/(2R0)である。式(11)は、角αが増加し、角θsが減少
するにつれて範囲分解能は向上することを示す。Δs=5μmの値、α=10°
およびθs=45°の場合、範囲分解能はΔz=20μmである。
解能Δzに影響を与える。2つのソース点がx−z平面にある場合、zの測定は
、不確定度Δyの影響を受けない。ψ=90°に関して、xにおける不確定度Δ
xは、zの測定において以下の不確定度を生じさせる。
への影響は最も少なくなる。
横方向の分解能と物体の最大深さとの間にトレードオフが存在する。物体の深さ
におけるこの制限を低減するための1つの方法は、異なった範囲の平面に連続し
て焦点を合わせ、焦点の深さの中にあるピクセルのみを用いることである。例え
ば、100μmの横方向の分解能は、フィールドの深さを1cmのオーダーに制
限し、10cmの範囲を有する物体は10個の異なった範囲で連続的に焦点を合
わせることによって全分解能で結像され得る。フィールドの深さの影響を最小化
するために、z軸は、範囲の広がりを最小化する、すなわち、物体の表面の平均
的平面に対して垂直の方向に規定され得る。横方向の分解能を損なわずに、結像
され得る横方向の面積を増加させるために、複数のカメラ(すなわち、検出器ア
レイ22)を用いて物体10の対象となる領域全体を覆てもよいし、または対象
となる個別のカメラを用いて検査してもよい。あるいは、単一のレンズの焦点面
に複数の検出器アレイを配置してもよい。これらのアレイは、独立して平行移動
され得、物体の種々の領域を高分解能で検査する。z軸に沿う個別の検出器アレ
イの平行移動または検出器アレイの傾動(tilting)は、異なった深さの
物体の領域に対して焦点を同時に合わせ、物体深さの受容可能な範囲を大きくす
ることができる。
スタティック角によってシャドウイング効果(shadowing effec
ts)が発生し得ることである。これらの効果は、図4におけるように、レンズ
をソースにより近づけて配置し、かつ検出器が結像面において横方向にオフセッ
トされる、レンズを軸から外した構成で用いることによって低減され得る。この
目的のためにレンズが設計される場合、またはレンズが十分に大きい視野を有す
る場合、軸から外して結像することによって生じる収差は最小化され得る。
放射ソースP2を含み、各ソースは発散ビーム150および154を提供し、そ
れらの放射ソースと、物体10の表面上の点P0との間にR1およびR2で表示
された光路長をそれぞれ有する。ソースP1およびP2は、単色レーザ等のコヒ
ーレント光の任意の適切なソースによって生成され得、このコヒーレント光は、
2つの点ソースP1およびP2を提供するために分割され得る。さらに、図2お
よび図2aのビームスプリッタの実施形態、ならびに図2bの光ファイバースプ
リッタの実施形態等の種々の技術が、コヒーレント放射ソースからの放射を分割
するために適切である。
されるべき物体10の表面の方向に向けられる。物体10の表面によって散乱さ
れる照射は、焦点調整素子またはレンズ158によって焦点合わせされ、検出器
アレイ22に当たる。レンズは、図4に図示されるように、軸から外された構成
で用いられ得、バイスタティック角に起因するでシャドウイング効果を低減する
。図1と関連付けられた上述のタイプの任意の偏光子(図示せず)は、焦点調整
素子158と検出器アレイ22との間に配置され得、検出器アレイ22に入射す
るスペックル画像または縞パターンのコントラストを向上させる。
めのプロセッサユニット28と通信する。制御ユニット32は、ソースP1を軸
160に沿って移動させるために、少なくとも可動ソースP1と通信する。上述
のように、制御ユニット32およびプロセッサユニット28は別個のデバイスに
よって実現され得るか、または代替的には、単一システムの部分であり得る。さ
らに、制御ユニット32およびプロセッサユニット28は、特定の用途において
所望され得るように、互いに通信し得る。
れた深さ座標zは、ソースP1およびP2のそれぞれから点P0へのビーム15
0の光路長R1と、ビーム154の光路長R2との間の差の関数R1−R2とし
て決定され得る。図5の実施形態において、可動ソースP1からの放射の位相は
、制御ユニット32の制御のもとで、ソースP1を軸160に沿って移動させる
ことによって変更される。この構成を用いて、点P0で強度における振動が生成
される。
) ここで、aは点ソースP1の平行移動の大きさを表し、ls、msおよびnsは
x、yおよびz軸それぞれに関する平行移動の方向を表す方向余弦である。 これらソースP1およびP2からの放射が点P0へ伝播した後に測定された、の
位相差は、以下の数式によって与えられる。
の位相オフセットを表す。P1が軸160に沿って平行移動すると、R1の値は
変化し、φをaの関数として変化させる。
またはP0を横切る干渉縞)の数は、以下の数式によって与えられる。
ス間の距離がaの場合の式(14)における光学位相差であり、φ(0)は、a
=0の場合の式(14)における光学位相差である。式(15)を考察すること
によって、ソースP1からの移動の結果である振動の強度の数Nは静止ソースP 2 の位置には依存しないことが明らかとなった。この非依存性は、ソースP1お
よびP2が互いに近接して配置されることを可能にする。この構成を用いて、ソ
ースP1およびP2それぞれからの発散ビーム150および154は、大気乱流
および振動等の一般的な妨害を受ける。このようにして、このような妨害の影響
は最小化される。さらに、ビーム150および154は、実質的に同じ偏光を有
した状態で物体10の表面に達する。
い縞パターンを生成するが、これは結像レンズ24によって容易に解消され得る
ことにも留意されたい。
で、式(15)は、以下のようにa/R0における二次まで近似され得る。
規定される線との間の角である。
定することができるということを示す。3つ以上の位置からψが分かれば、P0 の(x、y、z)座標は三角形分割によって決定され得る。次に、図1に対応す
る実施形態と類似の実施形態が説明される。ここで、x座標およびy座標は、検
出器アレイにおける結像点の位置から決定される。
の振動サイクルNを決定することによってか、またはそのような強度振動が生じ
る速度を測定することによって測定され得る。まず、サイクル数Nを決定するこ
とに基づくzの測定を考察されたい。既知のNを用いれば、zを除く式(15)
におけるすべての変数は分かる。zについて式(15)を解くことで、以下の数
式が得られる。
する。これは、P1によって移動される波長単位ごとの振動サイクルNに関する
スペックル強度の平均変調率を表す。aの値が0に近づくと、式(17)は以下
のように近似され得る。
略化され得、その結果、ソースP1の平行移動はx−y平面に限定される。この
構成は、後述されるように、ソースP1の平行移動に関する良好な実用的選択を
表す。結果として生じるzの数式は、以下のように表され得る。
て与えられる。
。瞬間振動率ρは、式(19)における平均振動率と類似の方法で以下のように
表され得る。
で以下の式が得られる。
6) aの値が小さい場合、ρは以下のように近似される。
および
)によってそれぞれ近似され得る。ただし、ρの代わりにρ\を用いる。本明細
書中において「\」とは、「\」の前の文字に上線を付していることと等価であ
り、例えば、ρ\は、
よって生じるzの測定における不確定度Δzを考察されたい。簡略化するために
、この計算は、式(16)によって与えられるNに関する近似式に基づく。Δz
を求めるために、zについてN(またはρ)の偏導関数をとり、この偏導関数を
比ΔN/Δz(またはΔρ/Δz)と等価すると、以下の式が得られる。
タである。Gに関する第1の式において、 l=x/R0,m=y/R0およびn=x/R0 (34
) は、点P0の方向余弦である。Gに関する第2の式において、θは極角およびφ
は方位角であり、球面座標系における原点からP0への方向を表す。同様に、ソ
ース点の平行移動の方向はθsおよびφsによって与えられる。
低下し、ソースP1の平行移動aの大きさが増加するにつれて、範囲分解能が向
上することが明らかとなる。式(33)の考察は、ジオメトリカルファクタGが
1と無限との間で変化することが分かる。ここで、1は達成可能な最良範囲分解
能に対応する。
mおよびnの所与の値に対して最小化されるようにls、msおよびnsを選択
することによって式(33)から得られる。この制約を適用すると以下の式が得
られる。
る線と直交し、その線およびz軸(φs=φ)によって形成される入射面にある
ことを意味する。式(35)の値を式(33)に代入すると、以下のようになる
。
生じることが分かる。これは、散乱点がx−y平面にあることを意味する。さら
に、z軸上にある分散点に対してGが無限大に近づくにつれ分解能が低下するこ
とも分かる。例えば、θ=30°の場合はG=2、およびθ=10°の場合はG
=5.76である。各点に関する平行移動の方向を変更することなく、画像内の
すべての点について式(35)を満たすことは可能ではないが、最適分解能の条
件は代表的な画像点について式(35)を満たすことによって近似され得る。
フセットの大きさaに依存する。測定ρに基づく技術の場合、走査中にρができ
る限り変化せず、その結果、値ρとzとは、ほぼ1対1で対応することが所望さ
れる。従って、標準的スペクトル解析技術は、ρの値を推定しzを決定するため
に適用され得る。走査時に生じるρの不均一性の程度を定量化するために、以下
のように定義される。
のみを残すと、以下の式が得られる。
形に増加することを示す。さらに、不均一性は、ψ=0°であるときに消失し、
ψ=90°であるときに有界を有することなく増加する。しかしながら、ψ=0
°の線におけるすべての点は、範囲、すなわち式(33)におけるG=∞に関係
なく、同じ変調率を有するので、ψ=0°のときに範囲分解能は存在しないこと
が分かる。従って、不均一性を最小化することと、最適範囲分解能との間にトレ
ードオフが存在する。
測定構成は、ns=0を設定し、φs方向におけるオフセット、すなわちφ=φ s を有する軸から外れた光学系を用いることである。その後、Gに関する式(3
3)は以下の式に通分される。
る物体を、R0=1mの距離から波長λ=0.7μmのレーザを用いて撮像する
ことが所望されることを仮定されたい。ns=0および物体の中心がθ=30°
およびφ=φsにおいて配置される場合、式(39)によって、ジオメトリック
ファクタGは、視野にわたって2.1と2.6との間で変化する。式(32)に
よって、a=5mmの平行移動は、カウントの20分の1、すなわちΔN=0.
05の振動の数における不確定度に対して、Δz=16μmの範囲不確定度を(
画像の中心において)生成する。走査全体に対する振動カウントの総数は、式(
16)によりN=3600である。画像の中心における変調の不均一度を推定す
るために、式(38)においてψ=60°を設定し、χ=0.0075を取得す
る。その結果、不均一度は、走査にわたって1%よりも少ない。この不均一度は
、走査中に、走査速度を少しだけ変化させ、測定中の周波数の任意の変化を補償
することによってさらに低減され得る。
程を示す。これらの工程を行って、オブジェクト上の各ポイント(x、y)にお
ける深さ座標zを判定する。プロセッサは、工程100から開始し、オブジェク
ト表面上の複数の照射ポイントによって散乱される放射の強度のパラメータを測
定する(工程108)。この情報から、各測定ポイントのz座標を計算する(工
程112)。
者に公知の適切なフィルタを挙げると、平滑化フィルタ、メジアンフィルタおよ
びカーブフィッティングフィルタがある(ただし、これらに限定されない)。そ
の後、マッピングされたポイントの表示または出力を当業者に公知の何らかの方
法で行うことができ、その後、本プロセスは図示のように工程124において終
了する。一実施形態において、工程120において、マッピングされたポイント
を、メッシュプロット上の演算されたz情報の関数としてプロットする。
り、これらの工程は、図1の実施形態と共に用いると適切である。工程108’
において、散乱される照射の強度をレーザ周波数オフセットの関数として測定し
、当業者に公知の方法の1つを用いてNを測定する。その後、工程110’にお
いて各ロケーション(x、y)について式(4)を用いてsを計算し、工程11
2’において各ロケーション(x、y)について式(9)を用いてzを計算する
。
実施形態を図6bに示す。この場合、工程108’’において測定される強度の
パラメータは、回数N(これは必ずしも整数ではない)であり、この回数の間、
可動ソースP1(図5)が平行移動するにつれ、強度は周期的に変化する(cy
cle)。当業者に公知の方法を通じて工程108’’においてNが判定される
と、工程110’’において、このNは式(19)によってρ\に変換される。
その後、工程112’’において式(17)および式(18)を用いてzを計算
する。図5の実施形態と共に用いられる処理工程108および112の別の実施
形態を図6Cに示す。この図6Cの実施形態では、工程108’’において測定
された強度のパラメータは、ソースポイントP1が平行移動するときに振動が発
生したときの瞬間振動速度ρである。工程112’’’において、式(28)〜
(30)を通じてρをzに変換する。
示す一実施形態において、検出器アレイ22の光検出器素子221、1〜22n 、m を逐次読み出す。検出器アレイ22のシリアル出力36は、プロセッサ28
への入力を提供する。プロセッサ28は、単一のプロセッサを含み得るか、ある
いは、複数のプロセッサを含むマルチプロセッサであり得る。
において、プロセッサ28は、複数のプロセッサ281、1〜28n、mを含む
マルチプロセッサである。検出器アレイ22の光検出器素子221、1〜22n 、m はそれぞれ、プロセッサ281、1〜28n、mのうち対応する1つに各出
力信号38を提供する。この構成を用いると、プロセッサ281、1〜28n、 m はそれぞれ実質的に同時に動作することができ、有利な性能を実質的に提供す
る。より詳細には、マルチプロセッサユニット28内の各プロセッサ281、1 〜28n、mは、光検出器アレイ22の対応する素子221、1〜22n、mか
ら受信したデータに基づいてz座標を計算する機能をする。従って、オブジェク
ト10の表面の各ロケーションのz座標を敏速に判定することができる。
し、これらの検出器およびプロセッサコンポーネントは、図1および図5のシス
テムにおいて、ユニット式の検出器およびプロセッサアレイ25の形態で用いら
れる。アレイ25は、共通基板上に作製されるかもしくは共通基板によって支持
されるか、または、マルチチップモジュール(MCM)として作製されるかもし
くは表面取付け技術(SMT)によって作製される。アレイ25の検出器部分は
光検出器素子221、1〜22n、mを含み、アレイのマルチプロセッサ部分は
、プロセッサ281、1〜28n、mを含む。より詳細には、図示のように、検
出器221、1〜22n、mはそれぞれ、プロセッサ281、1〜28n、mそ
れぞれと関連付けられるかまたはプロセッサ281、1〜28n、mの近隣に配
置され、入力信号を各プロセッサに提供する。プロセッサ281、1〜28n、 m は、検出器221、1〜22n、mそれぞれからの情報を実質的に同時に処理
して、深さ座標を判定する。
のアレイはオブジェクト10の表面に対向して配置されており、オブジェクト1
0の表面輪郭が測定される。この構成では、オブジェクト10の表面上のポイン
トP0からの散乱光を観察してzを判定するのではなく、光の位相シフトの測定
をオブジェクトの表面において直接的に行う。図示はしていないものの、図8の
システムは、ソースP1およびP2を制御する制御ユニット28と、検出器22
’に入射する放射を処理するプロセッサ28とを含む。検出器22’については
、図8中に図示しており、図1、5および6と関連して説明している。
ムには様々なものがあり得る。図9に示す一実施形態において、アレイ22’の
複数の個々の光検出器素子23を、対象エリア内のオブジェクト10表面上に配
置する。
いるスプリングアーム84上にアレイ22’の個々の光検出器素子23を取り付
ける。これらのスプリングアーム84を制御ユニット88によってオブジェクト
10の表面上で移動させて、特定のポイントまたは対象領域と接触させる。この
ようにスプリングアーム84を片持ち型で支持すると、アーム84がオブジェク
ト10上を移動している間、個々の検出器23はそれぞれ、オブジェクト10の
表面上のロケーションとの接触状態を保つ。すなわち、オブジェクト表面の輪郭
が変化すると、それに応じてスプリングアーム84も上下に移動する。
ことが理解される。例えば、さらなるソース(単数または複数)を用いて、オブ
ジェクトまたはその一部分に関するx、y座標情報を判定することができる。さ
らに、さらなる放射ソースを用いて、対象領域のシャドーイングに寄与するいか
なる処理精度の不足または処理のあいまい性をも低減することができる。
れる。例えば、これらの2つのポイントはどちらとも、反対方向に移動しながら
移動することもできるし、一定の間隔を保ちつつ同じ方向に移動することもでき
るし、共通の中心点の周囲を回転することもできるし、あるいは、制御システム
による切換えが可能なソースポイントのアレイを用いることによって動きをシミ
ュレートすることも可能である。
38)における変調の不均一性を低減できるなどの利点を得ることができる。こ
の不均一性を低減させることによって企図しているのは、Nがaの奇関数となる
ような様式でポイントを移動させることにより、式(27)中のrに関する線形
項を無くすことである。平行移動を2つのソース間で均等に分割して、これらの
2つのソース間の中心が固定されるようにする(これは、対向移動(oppos
ing motion)と呼ばれる)。その結果、以下の数式が得られる。
になる。
の大きさを1/4だけ低減する。
すと、以下のようになる。
よび式(23)となる。しかし、双方のソースを移動させることによってaへの
1次依存(first−order dependence)を無くすと、移動
ソースを1つだけ用いる場合よりもこの近似結果の精度が上がる。
式は以下のようになる。
、式(44)は以下のようになる。
(27)と比較してより安定する。ここで、式(37)において既に規定された
対応する変調不均一性χを以下のように表すことができる。
不在のため、得られるχの値は、式(38)のχ値よりもずっと小さい。双方の
ソースポイントを移動させることによって得られる利点の1つとして、レンジ分
解能を最適にした(ここでy=90°)場合、不均一性が小さいままであるとい
う点がある。
およびρの小a近似(small−a approximation)における
第1の項に基づく。これらの近似における第1の項は2つのソースの動きと同一
であるため、先に述べた結果は、2つのソースが移動する場合にも当てはまる。
8によって送信ビーム46および反射ビーム50に分割される。送信ビーム46
は、鏡94からレンズ58へと反射される。同様に、反射ビーム50も鏡96か
らレンズ52へと反射される。鏡付き表面91および93を備えた直角プリズム
92を用いて、ビーム46および50をそれぞれレンズ58および52から再度
方向付ける。ビーム46および50はそれぞれ、ポイントP1およびP2に集束
し、発散ビーム82および86として延びる。位相シフトは、別個のビーム46
および50がたどるビーム経路間の光路長の差を変更することによって達成する
ことができる。これは、例えば、素子48、94または96のうち少なくとも1
つを圧電変換器(図示せず)によって面外平行移動させるようにすることにより
、達成可能である。あるいは、位相シフトデバイス(図示せず)を、ビーム経路
46または50のうち1つに挿入してもよい。このようなデバイスの一例として
は回転するガラス製ディスクがあり、このガラス製ディスクは、別個の工程にお
いて回転角度を用いて自身の光学的厚さを変化できるようにコーティングされて
いる。プリズム92は軸99に沿って平行移動することが可能であり、これによ
り、ソースP1およびP2の対称な動きに影響を与える。あるいは、レーザ40
、ビームスプリッタ48、鏡94および96ならびにレンズ52および58が1
つの群として軸99に沿って平行移動する間にプリズム92を固定状態にして、
ソースP1およびP2を対称に平行移動させてもよい。ソースP1とソースP2 との間の距離の変更は、レーザ40とビームスプリッタ48との間のビーム44
中に配置された角度変調デバイス(図示せず)を用いても達成することが可能で
ある。このような角度変調デバイスは、ビーム角度を変化させる。角度変調デバ
イスの例を挙げると、検流計鏡および音響光学変調器がある。
スである。すなわち、ソースP1およびソースP2は、ビーム46および50が
集束状態で通過する場所である光または領域の小ソースではない。その代わりに
、ビーム46および50を発散させる光コンポーネント(すなわち、負レンズ)
を用いて発散ビーム82および86を生成する。その結果、オブジェクトへの照
射を、実際のソースP1およびP2を用いた他の実施形態から区別不可能にする
ことが可能となる。
成を示す。レーザソース40からのビーム44は、ビームスプリッタキューブ4
8によって送信ビーム46および反射ビーム50に分割される。送信ビーム46
は、曲面鏡95から反射され、直角プリズム92の鏡付き表面91から反射され
る。同様に、反射ビーム50も、曲面鏡97および直角プリズム92の鏡付き表
面93から反射される。ビームスプリッタキューブ48を改変(例えば、ビーム
経路外部のガラスを除去)して、光コンポーネントを密接に取り付けることが可
能である。曲面鏡95、97は、図10中の平面鏡94、96およびレンズ58
、52の機能を行う。これらの曲面鏡95、97を軸のずれた(off−axi
s)放物線型鏡セグメントにして、P1およびP2において小さなフォーカルス
ポットを生成させてもよい。あるいは、曲面鏡95、97を球状の素子にしても
よい。球状の素子によって生じる光学収差は、縞接触には影響を与えず、縞位置
のみに影響を与える。これらの光学収差は、解析に取り入れることが可能である
。位相シフトは、ビームスプリッタ48を平行移動させるか、または、上述した
ようなビーム経路46もしくは50のうち1つの中に位相シフトデバイス(図示
せず)を配置することにより、達成可能である。
ージング構成に基づく。焦点面中のイメージポイントPiの側方位置は、レンズ
24の倍率によってオブジェクトポイントP0の側方位置に関連付けられ、レン
ズ24とオブジェクト10との間の距離から独立する。焦点はずれによってイメ
ージPiは不鮮明となるが、これは、検出器22においてイメージPiがはっき
りと側方位置にある状態には影響を与えない。そのため、オブジェクトポイント
P0の側方位置をそのイメージポイントPiの位置によって推定する作業をP0 とPiとの間の距離に関係無くおこなうことができる。残念なことに、オブジェ
クトが大型である場合、テレセントリックイメージャ24は用いることはできな
いかもしれない。その理由は、レンズ24の直径は、少なくとも測定対象オブジ
ェクト10の最大寸法と同じくらいの大きさでなければならないからである。従
って、ほとんどのレンズはテレセントリックではない。
、z)座標を判定するための一般的なアプローチとして、検出器22におけるイ
メージポイントPiの位置は、オブジェクトポイントP0とフォーカシング素子
24の光学軸98との間の角度αxを表すものであると仮定するアプローチがあ
る。座標は、縞数Nと、検出器22におけるイメージポイントのロケーションと
から導出することができる。以下の方程式は正確である(小aは仮定していない
)。 側方座標は、以下によって得られる。
た角度の成分である。 z座標は、以下によって得られる。
し、(xm、ym、zm)が座標系の起点以外の場所にある場合に一般化された
解法を提供する。多用性を広げるために、Pmおよび他のパラメータを経時変化
させてもよい。
例えば、CCDアレイ)と連絡するプロセッサユニットからのデータを受信する
コンピュータ上のソフトウェアにおいて実施することが可能である。あるいは、
プロセッサユニット28上にアルゴリズムを直接インプリメントしてもよい。
ためには、角度αxおよびαyだけではなく、正規化された縞数ρ\(例えば、
式(19)によって規定された縞数ρ\)もわかっていなければならない。上記
の実施例において、これらの角度は、焦点面22内のイメージポイントPiのロ
ケーションによって判定している。αxおよびαyを正確に判定するためには、
レンズ24内に発生する可能性のある任意のイメージひずみを修正することが必
要である。ひずみの修正は、例えば、角度αxおよびαyをイメージポイントP i 上にマッピングするひずみ機能を生成することにより、行うことが可能である
。ひずみマップは、レンズが変わると変動し、また、同じレンズ24を異なるイ
メージング構成において用いた場合にもわずかに変動する可能性がある。
要な残りの2つの方向成分を判定することにより、レンズひずみを特徴付けると
いう不便な作業を無くす。例えば、図12aに示すように、さらなるソースヘッ
ドを位置P1に配置し、以前P1にあった(図12を参照)レンズ24を、オブ
ジェクト10上の対象領域をビューイングする際に便利の良い任意の位置に配置
する。第2のソースヘッド内のソースP1xおよびP2xがx軸に平行に方向付
けられている場合、このソースヘッドのポイントP0における縞数Nxを測定す
ることは、ポイントP1xおよびポイントP1ならびにポイントP1およびP0 (式(41)を参照)によって規定されたライン間の角度yxを測定することに
相当する。その後、式(47)〜(49)に必要な角度αxを、αx=π/2−
yxの関係を用いてyxから判定する。同様に、(図12aでは図示していないが
)ソースP1yおよびP2yをy軸に平行に方向付けた状態で位置P1において
ソースヘッドを用いることにより、角度αyを判定することも可能である。ソー
スヘッドをこのように方向付けた場合、縞数Nyから、αx=π/2−yxを通
じて角度αyを測定することができる。角度αxおよびαyを提供するソースヘ
ッドは、例えば、z軸周囲を90°回転する同一のソースヘッドであるか、また
は、ビームスプリッタを用いることによって同一のロケーションに現れるように
された別のソースヘッドであり得る。
2のロケーション(またはオブジェクト10上に配置された検出器23のロケー
ション成分)を知る必要は無い。イメージングシステム22および24(または
検出器23)の目的は、P0における照射強度を観測することである。検出器ア
レイ22中の各ピクセルは、観測ポイントP0を表す。そのため、各ピクセルに
ついて(x、y、z)座標を判定することができる。様々なアプローチを用いて
、異なるソースヘッドに対応するNの測定を別個に行うことが可能である。例え
ば、動作波長が異なる複数のソースヘッドを同時に照射することによって多重化
された異なるソースヘッドまたは波長を逐次照射することにより、測定値を時分
割多重化(time multiplexed)することが可能である。偏波ダ
イバーシチを用いることにより、測定値をさらに区別することも可能である。
式(41)を再度参照する。典型的な用途の場合、αはR0よりも小さく、式(
41)は以下のように簡約することができる。
スの間隔aおよび波長λまたは周波数νへの依存を式(57)から無くす。
、3つ以上のソースヘッドからの三角測量によって得られた測定値を通じて判定
することが可能である。
ロケーションそれぞれに対応するソース対の方向は、方向コサイン(ls1、m s1 、ns1)、(ls2、ms2、ns2)および(ls3、ms3、ns3 )によって得られると仮定する。デカルト座標および方向コサインについて見る
と、式(58)は以下のようになる。
判定は、各ソースヘッドについて正規化縞数ρ\1、ρ\2およびρ\3を測定
し、式(59)によって表される3元連立方程式を解くことにより、達成される
。式(59)の形式主義は、一般化されたソースヘッドを配置および方向付けを
コンパクトな形態で取り扱うため、便利である。
ヘッドによって照射されない場合、周囲のピクセル22について計算された方向
間で補間を行うことにより、当該ピクセル22の方向情報の欠損を補うことがで
きる。さらに、複数のソースヘッドから得られた方向情報をピクセルロケーショ
ンから得られた方向情報と組み合わせて、測定品質を向上させることもできる。
、倍率およびひずみ)が分からなくてもP0を判定することが可能であるため、
複数のイメージャ22および24(分かりやすくするため、イメージャは1つし
か図示していない)からの情報を組み合わせる作業が容易になる。これが可能な
のは、各カメラについて生成された点群データ(point cloud)は既
に同じ座標系にあるため、較正プロシージャを行ってデータを登録しなくてもよ
いからである。複数のイメージャ22および24を配置して、特定の対象領域を
任意の所望のビューイング角度および方位分解能でビューイングすることが可能
である。そのため、1つのイメージャ22および24によってオブジェクト10
全体の概要を得て、他のイメージャ22および24によって対象領域を高分解能
で検査することが可能である。
は両方の動きによって生じる強度−振動をカウントすることにより、オブジェク
ト10の表面上のポイントP0における縞数Nを判定する。サイクルのカウント
が最も基本的な様態で行われる場合、カウントを追跡するためには、1サイクル
あたり少なくとも2回強度をサンプリングしなければならない。また、コストを
最小限にするためには、標準的な市販のコンポーネント(例えば、フレームグラ
ッバーに接続された電荷結合素子(CCD))を用いてサイクルをカウントする
と好ましい場合がある。カウント数が延びるほど、イメージフレーム中のポイン
トP0はN=0の位置から来るため、このアプローチを用いると、イメージフレ
ームを何百単位で取得および処理しなければならなくなる可能性がある。
ウント技術を用いることによって解消することが可能である。この技術は累進縞
分割(progressive fringe division)と呼ばれ、
これは、強度変調の周期的性質および予測可能な性質を利用したものである。そ
のため、ソース間隔aの複数の別個の値における強度変調をサンプリングするこ
とにより、Nを判定することが可能である。
間隔(Col.1)から最終微細間隔(Col.3)まで変化させる工程からな
る。各工程において、ラップされた(wrapped)位相マップまたはサイク
ルマップ(ロウ(ROW)4)を位相シフトした縞パターン(ロウ1−3)から
生成する。ブートストラッピングオペレーションにおける先行するアンラップの
位相マップに基づいて、連続するラップされた位相マップをそれぞれアンラップ
状態にする。工程間の縞間隔の比は典型的には、10以上のオーダーであるため
、少数のデータフレームを用いて極めて高密度の縞パターンをアンラップ状態に
することができる。高密度の縞パターンから位置座標を計算すると、精度が向上
する。累進縞分割を行うことによって得られる1つの利点は、オブジェクト10
の表面の複雑性に関係無く、アルゴリズムによって縞の順序の番号を各縞に正確
に割り当てるという点である。例えば、オブジェクト10は、表面の不連続部分
、大型の断崖形状部分および穴部を含み得、または、複数のばらばらのオブジェ
クトからなり得る。累進縞分割は表面の連続性に依存しないため、累進縞分割を
用いて、空間中に分散した未接続の複数の別個のポイントの位置をロケートする
ことが可能である。
の事実を考えると有用である(簡潔にするため、以下の議論において、ソースP 1 およびソースP2は対向移動すると仮定するが、以下の結果は、他の状況(例
えば、1つのソースP1のみが移動し、静止状態の周波数がチューニングされた
ソースを用いる)にも当てはまる)。第1に、ソース間隔が初期値であるゼロか
ら最終値aまで増加するときにポイントP0において発生する強度振動Nの数は
、最終値aのP0における縞数Nに等しい。例えば、10が単純なオブジェクト
(例えば、面)であった場合、N=0の位置(これは、R1=R2において発生
する)とポイントP0との間のオブジェクトの表面上のN個の縞をカウントする
ことができる。Nの値は、部分的サイクルまたは部分的縞を含むため、必ずしも
整数ではない。
ときに縞がポイントP0を通過してP0とN=0との間の領域内部に集中する縞
の数としてNを解釈することも可能である点に着目することにより、この第1の
事実をさらに説明することができる。言い換えると、a=0である場合、この領
域内に縞は無く、この領域内部に含まれる縞の数は、縞がP0を通過するたびに
1つずつ増加する。このように縞数Nと振動サイクルとが対応している状態は、
ロバストな三次元イメージング技術への鍵であり、これにより、ポイントP0周
囲の表面の連続性について仮定を行う必要なく、任意のポイントP0におけるN
を判定することが可能となる。これにより、表面中の不連続部分(例えば、穴部
および断崖形状部分(すなわち、エッジ))をあいまいさ無く処理する。
る。その結果、aの実行(running)をゼロからその最終値まで行ってN
(a)を判定して全体的スキャンを完了する必要が無くなる。間隔がa2である
場合のNの値は、以下の関係を用いて、間隔a1における値から推定することが
可能である。
ていると、それらの値を用いて、より小さな差間隔a2−a1におけるNの値を
推定することが可能であるということである。
行うことを可能にすると、視野内の各測定ポイントP0において、部分的サイク
ルまたは「ラップされた」コンポーネントの縞数Nの高分解能マップが得られる
点である。こうすると、以下の式により、ラップされたコンポーネントをその位
相から計算する。
場合、Nwrappedの範囲は±0.5となる。位相シフト測定ではNwra pped を高精度で判定することが可能であるが、Nが判定されないと、整数値
のオフセットが残る。
、ρ\を高精度に判定することが可能である(ただし、縞がうまくアンラップ状
態になっていると仮定した場合)。例えば、微細な縞ほど、位相ステップエラー
が表面プロファイルの測定に与える影響を低減する。しかし、微細な縞ほど縞密
度も高いため、アンラップを行うのが困難である。累進縞分割によって得られる
利点は、極めて微細な縞を正確にアンラップして、可能な分解能のうち最も高い
分解能を達成することができる点である。
ウ4に示すラップされたサイクルマップ)を、複数の別個のソース間隔aそれぞ
れにおけるイメージング対象表面について生成する。これらのサイクルマップを
3つまたはそれ以下の縞サイズ(すなわち、粗大、中程度および微細)について
生成すれば十分である場合が多い。ブートストラッピングオペレーションでは、
アンラップのサイクルマップを連続的に用いて、整数オフセットにエラーを導入
することなく、次のサイクルマップをアンラップする。最も微細な縞パターンに
対応するアンラップのサイクルマップによって、最も高い分解能での測定が得ら
れる。
ージャは、(工程150)から開始し、ソース間隔を初期値a1で設定する(工
程154)。この値は通常は、一連のソース間隔の中で最小であり、オブジェク
ト10上に粗大縞パターンを生成するような値が選択される。その後、式(6)
に示すような位相ステッピング(phase−stepping)測定を行うこ
とにより、ラップされたサイクルマップを生成する(工程158)。式(62)
によって位相をサイクルに変換することにより、工程158を終了する。
された推定値Nを入手する(工程162)。工程162は、工程158と同時に
行ってもよいし、工程158の前に行ってもよいし、または工程158の後に行
ってもよい。しかし、工程158および工程162はどちらとも、次の工程に続
くような様式で行わなければならない。この初期推定値を生成するアプローチに
は多くのアプローチがある。例えば、a1の大きさを、オブジェクト全体が1つ
の縞内におさまる(例えば、−0.5<N≦0.5)くらいに小さくすることが
でき、これによりアンラップ作業を不要にする。あるいは、オブジェクト10に
わたって少数の縞が生成されるようにaの大きさを十分に大きくすることができ
、これにより、Nの推定値がアンラップされる全ての測定ポイントにおけるサイ
クルの半分の範囲において正確である限り、視野にわたってNが変動する際の単
純なモデル(例えば、線形傾斜(ramp))が、このNの推定値として十分と
なる。初期推定値を得るための方法の別の例は、一般的な形状または予測される
形状を仮定して、式(40)に基づいてNを計算する方法である。他の例は、先
行測定結果を同一のまたは類似するオブジェクト10に用いる方法、または、サ
イクル計数を用いて視野をサイクルに区分けする方法である。a1に対応する縞
の数が小さいため、サイクル計数に必要なフレームの数も少数になる。
プされたサイクルマップをアンラップすることである。工程162からのアンラ
ップされた推定値が1サイクルの半分以内の範囲で正確である場合、以下の数式
を用いてアンラップを達成する。
フセットに過ぎない点に留意されたい。この項におけるアンラップされた推定値
からラップされた測定値を減算すると、階段状の関数が得られる。この階段状の
関数の値は、整数オフセットにほぼ等しい。最も近い整数まで丸めを行うと、1
/2サイクルの制約内にエラーが収まっている限り、数式は正確になる。1/2
サイクルの制約内にエラーが収まらなくなった場合、式(63)は、Nに整数オ
フセットエラーを導入し、このエラーは、視野にわたって変動し得る。
。測定の品質をさらに上げるために、工程170において新規のソース間隔a2 を設定し、その結果、より微細な縞パターンがオブジェクト10上に生成される
。その後、新規の位相ステッピング測定を工程174において終了して、より微
細な縞パターンに対応する新規のラップされたサイクルマップを得る。この測定
値をアンラップするためには、a2に対応する新規のアンラップされた推定値が
必要となる。工程178において、式(60)を用いて先行するアンラップされ
た推定値(これは、工程166からソース間隔a1について得られる)を適切な
レンジにスケーリングすることにより、この推定値を得る。スケーリングが行わ
れると、工程174からの新規のラップされたサイクルマップが工程182にお
いてアンラップされる。このアンラップ工程は、このサイクルマップと、工程1
78からのスケーリングされた推定値とを式(63)に代入することにより、行
われる。
は、より大きなソース間隔a2に対応する。ここで、最適なρ\の分解能が得ら
れるまで、徐々に微細レベルが高くなる縞(すなわち、ソース間隔が大きくなる
)を用いて工程170〜工程182を繰り返すことが可能となる(このループが
終了するたびに、ソース−間隔の下付き文字が1だけインクリメントされる)。
通常は、最適な分解能はこのような繰り返しを1〜3回行うと得られる。典型的
には、新規のソース間隔と古いソース間隔との間の比は10以上のオーダーであ
る。この比がとることのできる最大値を実際に判定する場合、Nの推定値と実際
の値との間のエラーが1/2サイクルの限定内におさまることを要求することに
より、判定を行う。所望の数の反復が完了した後、図14中の工程170〜工程
182によって表されるループを終了する。工程186において、式(19)を
通じてNをρ\に変換する。測定の目的が全ての三次元座標を判定することであ
る場合、例えば、式(47)〜(49)または式(59)を用いて三角測量を行
って、各測定ポイントについて座標を入手する(工程190)。工程194は、
累進縞分割の終了を示す。図14は累進縞分割の前に行われるプロシージャを例
示したものに過ぎず、コンポーネント工程は、複数の異なる順序によって実行す
ることが可能であることが理解されるべきである。例えば、全てのデータを収集
した後に計算を始めることが望ましい。
きる能力を得るためには、ソース間隔を縞間隔と同じ比率で変化させなければな
らない。ここで、他の2つの測定値から所望の間隔の縞パターンを合成する技術
について説明する。このアプローチによって得られる1つの利点として、縞パタ
ーンを合成する際に用いられる2つの測定値をより容易に入手することが可能で
ある点がある。
Nwrapped(a2)−Nwrapped(a1)の差をラップすることに
相当することを観察することにより、得られる。言い換えると、式(64)中の
角括弧内のラップされた数量間の差は±1の範囲にあり、この差を正しい範囲で
ある±0.5以内に収めるためには、この差をラップしなければならない。
個々の位相ステッピング測定からの差a2−a1に等しいソース間隔におけるラ
ップされたNの位相ステッピング測定を合成する手段である。この方程式は、a 2 −a1の差が小さい場合に顕著な特性を有する。合成された測定値を構成する
際の個々の測定は不規則になり得、多くの不連続部分を含み得るが、視野にわた
って式(64)の左辺を平滑化して、ラップする不連続部分のうち差a2−a1 に対応する部分の数のみを得ることができる。そのため、縞合成は、累進縞分割
を開始するための初期位相マップを生成する効果的な技術である。また、縞合成
を用いて、累進縞分割に用いられるシーケンス中の残りのアンラップの位相マッ
プのうち任意の部分を合成することも可能である。合成されたサイクルマップか
らρ\を計算する場合、式(19)を改変して以下のような式を得る。
を可能にする1つの方法を示す。工程154’は、工程154を改変して、間隔
a1に加えてベース間隔a0を選択するようにした工程である。工程158’は
、工程158を改変して、間隔値a0およびa1の両方におけるララップされた
サイクルマップを測定するようにした工程である。工程160’(これは、工程
158の後に挿入される)および工程176’(これは、工程174の後に挿入
される)はさらなる工程であり、式(64)を利用して、差an−a0に対応す
るラップされたサイクルマップを合成する工程である。工程186’は、工程1
86を改変した工程であり、この工程186’を用いると、式(65)を用いた
合成されたサイクルマップからρ\を計算することができる。
を生成するシステムへの要求を低減できる点である。例えば、最も広範囲なソー
ス間隔を除々に少しずづ変化させると、極めて小さなソース間隔に対応する縞を
生成することができる。このアプローチは、移動要件を低減するだけではなく、
図10および図11において説明したような、ソースヘッドのいくつかの制約の
可能性も解消する。例えば、これらのソースヘッドを用いると、極めて小さなソ
ース間隔を生成するのが困難となり、ソース間隔を生成しようとすると、振動に
よって縞が不安定となり得る。縞合成によって得られる別の利点は、ソース−間
隔較正の絶対値を維持するよりも2つのソース間隔間の一定の差を繰り返す方が
通常容易であるため、測定の反復性および精度を向上させることができる点であ
る。
きエラーを最小限にする点である。これらの動きエラーが発生し得るのは、例え
ば、図10および図11中のフォーカルスポットP1およびP2がa=0のノミ
ナル位置においてオーバーラップしない場合である。考慮すべきアライメント不
良は2種類ある。第1に、フォーカルスポットP1の1つが他のP2に対してd
の量だけ図面上から移動したと仮定する。これらの2つのフォーカルスポットP 1 およびP2はそれぞれ直線状の軌跡を追随し、これらの軌跡が平行であっても
、これらの2つのスポットは、その最近接アプローチ(すなわち、名目上のa=
0の位置)において距離dだけ欠損する。スポットP1およびP2が最近接アプ
ローチを通過すると、オブジェクト10上の縞パターンの方向付けは、ノミナル
または所望の方向付けに対して90°だけフリップする。さらに、縞間隔は、ソ
ース間隔dに対応する縞間隔よりも大きくならない(ソースヘッドのアライメン
トは、縞がa=0近隣に適切に方向付けられた様態で保持されるまで変位dを調
節することにより、容易にされる)。
が、他のフォーカルスポットP2に対して動き方向99に沿って変位する。この
種類のアライメント不良が発生するのは、プリズム92が横方向に変位するため
である(すなわち、この変位は軸99に対して垂直であり、図10および11内
の面にある)。最近接アプローチにおいて、1つのソースP1が、他のソースP 2 の後側に直接配置されており、オブジェクト上に同心円状の縞パターンを生成
する。aが増加すると、これらの縞の曲率は小さくなり、自身の通常の外見に近
づく(ソース−ヘッドのアライメントをとっている間のこのエラーソースの最小
化は、小aについてプリズム92を横方向に調節して縞曲率を最小化することに
より、容易化される)。
値を生成するために式(60)が用いられた際に累進縞分割の結果得られる位相
エラーが大きくなることに起因する。動きエラーがあると、Nの推定値に半サイ
クルエラーを導入することなく用いることが可能なaの比の大きさが限定される
。縞合成を用いると、動きによって誘発される位相エラーはどちらの測定におい
ても共通し、減算において相殺されるため、この制約が解消される。そのため、
所望の形状から変形する縞をこの測定により生成し得る場合にも、合成された縞
は通常のものである。
ーム間のグローバルな位相オフセットエラーを自動的に補償するというさらなる
利点が得られる。これらの位相−オフセットエラーは、例えば、2つの別個のビ
ーム間の経路−長さの差S0によって生じ、波長λの複数倍の整数ではない。位
相オフセットエラーがあると、0°以外の縞位相(すなわち、強度が最大になる
場所以外の場所)において、N=0の位置が発生する。縞合成を行うと、これら
の位相エラーは自動的に相殺され、その結果、縞パターンの位相を較正する必要
が無くなる。あるいは、縞合成によってN=0の正確な位置と判定し、そのポイ
ントにおける強めあう干渉について位相シフト要素を調節することにより、位相
較正を達成することも可能である。
P2からのオブジェクト10上に降下する照射パターン中の波面エラーを相殺す
る傾向となるという別の理由においても重要である。波面エラーは、ソースの不
完全な使用または長期にわたる使用に起因し、縞パターン中に摂動またはリップ
ルを生成する。これらの摂動は、測定時に小さな高さエラーに変化する。波面エ
ラーの相殺が縞合成において発生する理由は、これらのエラーはソース間の距離
からほとんど独立しているためである。そのため、エラーは、各ソース間隔にお
いて行われる測定に共通し、式(64)における減算によって相殺される。縞合
成は、照射欠陥によるエラーを低減する有効な手段であり、ソース−ヘッド光学
部品の品質への要件を低くすることにより、システムコストを低減する機能をし
得る。
範囲を小さくすることもでき、そのため、ソースヘッド設計をより簡単にするこ
とができる。一般的には、ソースヘッドは2種類の縞変調(すなわち、縞サイズ
の変動および位相シフトから生じる縞パターンの横方向の動き)を提供する点に
留意されたい。位相シフトのみがある場合、ソースの位置を変更することなく位
相を変更する。しかし、位相を変化させるための別の方法は、1つのソースP1 を他のソースP2に対してオブジェクトに向かってまたはオブジェクトから離れ
て移動させる方法である。このような方法を行うとソースP1のロケーションが
変化するが、ソース間の間隔が不十分である場合、この変化は極めて小さく(1
波長未満)、その結果得られるソース対P1およびP2の方向付けの変化も小さ
い。縞合成の間ソース間隔は大きいままにすることができるため、ソースポイン
トP1およびP2の1つまたは両方の小さな動きを通じて、どちらの種類の変調
も達成可能である。このように動き要件が少ないため、単純なソースヘッドを設
計することが可能となる。例えば、ソースヘッドは、図2bに示すファイバ−光
学部品−スプリッタ構成からなり得る。縞サイズの位相シフトおよび変動は、フ
ァイバの一方の端部を直交方向に移動させることにより、達成される。
能にする重要なさらなる利点を有する。縞合成が無いと、累進縞分割を測定に適
用する際、累進縞分割に必要な縞サイズの大きな変動を達成するために極めて細
密なチューニング帯域が必要となる。縞合成は、従来のチューニングレンジから
ずっと離れた縞間隔変動を生成することを可能にする。このチューニングアプロ
ーチを用いると、ソース動きがなくなり、周波数の測定制御を高い精度で行うこ
とができるため、反復性および精度が極めて高い測定値が得られる。
を参照して説明する。Nは、ソース間隔aに依存するのと同様に、レーザ周波数
νに機能面において依存する(式(1)においてS0によって表される経路−長
さ差はゼロであり、そうでない場合、方程式においてこのオフセットを考慮する
必要があると仮定する)。従って、Nが周波数νの関数として書かれている場合
、式(60)、(63)および(64)は、累進縞分割および縞合成の際の基本
となり、周波数チューニングにも適用される。周波数チューニングを支配する数
式は以下のようになる。
に行われる工程を示す(νの代わりにaが用いられる場合、図15は縞合成およ
びソース動きを用いた累進縞分割にも適用される点に留意されたい)。この図に
おいて、プロセスの開始時に全てのデータを収集する。プロシージャが開始する
と(工程200)、ラップされたサイクルマップNwrapped(νo)、N wrapped (ν1)...Nwrapped(νn)の測定(工程208)
において用いられるレーザ周波数ν0、ν1,...νnを選択する(工程20
4)。その後、方程式(67)を用いて、周波数差ν1−ν0、ν2−ν0,.
..νn−ν0について、ラップされたサイクルマップを合成する(工程212
)。その後、多くの可能なアプローチのうちの1つ(例えば、図14の工程16
2)を用いて、第1の差ν1−ν0についてアンラップのサイクルマップの推定
値を得る(工程216)。好適なアプローチは、N(ν1−ν0)を±0.5−
サイクルの制限内におさまらせるくらいに十分に小さな差ν1−ν0を用いて開
始する。工程220において、式(68)を用いて、シリーズ中の次に大きな差
νi−ν0についてアンラップのサイクルマップを計算する。工程224におい
て式(66)を用いてこのアンラップのサイクルマップをスケーリングし、その
後工程220にフィードバックして、シリーズ中の次の合成されたサイクルマッ
プをアンラップする。最終周波数差νn−ν0に対応するサイクルマップがアン
ラップされた後、ループは終了する。工程228において、式(70)を用いて
、最終アンラップのサイクルマップNを縞数ρ\に変換する。最終サイクルマッ
プNが合成されなかった場合、方程式(69)を用いる。工程232において、
例えば、式(47)〜(49)または式(59)を用いて、ρ\に対応する方向
情報を各対象ポイントP0について三次元座標に変換する。
さを解消し計数推定値を細密化する強力な一般的な方法であり、これらの技術は
上記にて述べた用途以外の用途にも適用可能であることが理解されるべきである
。詳細には、これらの技術は、データ取得のための上記の方法および装置または
縞解析にも限定されない。
スペックルに対して累進縞分割法および縞合成法を用いた改変例について説明す
る。図16において、レーザビーム44はビームスプリッタ48に入射し、ビー
ムスプリッタ48は、ビーム44をビーム46および50に分割する。ビーム5
0は、ビームスプリッタ48を反射し、オブジェクト10上の対象領域を照射す
る。ビーム46は、ビームスプリッタ48を通過して、基準表面102を照射す
る。オブジェクト10および基準表面102から反射された光は、ビームスプリ
ッタ48によって再度組み合わせられて、レンズ24を通過する。レンズ24は
、検出器アレイ22上にスペックルイメージを生成する。別個のビーム経路のう
ちの1つに位相シフト要素(図示せず)を挿入することにより、位相シフトを達
成する。レーザビーム44の周波数は可変であるか、または、2つ以上の別個の
値の間で切り換わることができる。オブジェクト10と基準表面102の仮想イ
メージ102’との間の高さzの差による経路長さの差により、スペックル強度
が変調し、これはレーザ周波数と共に変化する。特定のレーザ周波数における変
調の位相を、位相シフト測定を通じて正確に判定することが可能である。
ェクト10を照射しているのは1つのレーザビーム50だけであるためである。
2つのビーム46および50からの光のコヒーレントの重ね合わせにより、検出
器22において干渉が発生する。式(1)〜(5)によれば、レーザ周波数が変
化すると、オブジェクト10のイメージ中の個々のスペックルローブの強度が周
期的に変動し、その際の変動速度は、経路長さの差sに比例する。簡潔にするた
め、基準表面102は平面であり、x−y面にあると仮定する。zはオブジェク
ト10上のポイントP0の仮想基準表面102’上の高さを表すものとし、この
zを測定する。名目上はs=2z(これは、伝播の往復による)であり、以下の
式によってNからzを計算することができる。
すものと同じであるが、ただし、工程232において、Nをzに変換する際に式
(71)または式(72)を用いている点においてのみ異なる。図16に示すア
プローチは三角測量に依存していないため、ソースP1およびP2と受信器22
および24との間のシャドーイング効果を無くすことができるという利点が得ら
れる。このアプローチは特に短距離の用途に有用であり、そのような用途として
は、経路長さの間隔が短い三次元顕微鏡法がある。
オブジェクト10のイメージを形成する際にレンズ24を用いる実施形態の場合
、得られるイメージにスペックルがあり、スペックルによってレンジ分解能が制
限され得るという不利点が生じ得る。オブジェクトの表面上に検出器23を直接
配置すると、この制約から逃れることができる点に留意されたい。スペックルが
発生するのは、コヒーレント照射を用いてイメージングを行った場合であり、そ
の理由としては、レンズ24の有限ポイントスプレッド(finite poi
nt−spread)機能により、オブジェクト10の表面上のP0周囲の領域
から光が散乱し、検出器22におけるポイントPiにおいて干渉するからである
。検出器面におけるスペックルのサイズの平均は、レンズのポイントスプレッド
機能のサイズによって決まる。スペックルにより、P0周囲の散乱領域による寄
与は、イメージポイントPiにおいて徐々に大きくなる構成にすることができ、
これらの散乱領域により、ρ\の測定値は、一方向または反対方向に引き寄せら
れる。
る。実際は、検出器アレイ22の単一のピクセルに対応するエリア内におさまる
スペックルの数が多数である場合、スペックルによる効果は平均化する傾向とな
る。この条件を満たすための1つの方法として、レンズ24によって干渉縞を分
解し、レンズ24のポイントスプレッド機能の大きさを、検出器アレイ22の各
素子上におさまるスペックルの数を多くするくらいに十分に小さくすることを要
求する方法がある。このアプローチは、方位分解能を極限まで使用しなくてもよ
い状況では良好に機能する。一般的には、方位分解能とスペックルの平均化とは
トレードオフする。
より近密に整合させる必要が出てくる場合がある。一般的には、レンズのポイン
トスプレッド機能のサイドローブの大きさを小さくすることが望まれる(すなわ
ち、アポダイズする)。そうすると、ポイントスプレッド機能の強力なサイドロ
ーブ内に含まれる散乱領域の能力が低減する。このような処理を行わないと、こ
のような散乱領域の能力により、ρ\の測定に偏りが出てくる。アポダイゼーシ
ョンを行うと、方位分解能はわずかに低下し得る。場合によっては、この方位分
解能とレンジ分解能の向上とをトレードオフすると有利である場合がある。
に対応するオブジェクト10の表面領域上でのρ\の変動を最小限にする技術が
ある。例えば、ピクセル密度を増加させつつピクセル22あたりのスペックルの
一定比(constant ratio)を維持すると、ピクセル22に対応す
る領域の片側から他方の側へかけてのρ\の変動が低減し、レンジ分解能および
方位分解能の両方が向上する。
ドの位置および方向付けを適切に選択することにより、個々のピクセル22あた
りのρ\の変動を最小化することができる。ρ\の変動が最小化されるのは、例
えば、ピクセル22に対応する表面素子が三次元縞パターンの等位相面に沿って
延びる場合である。言い換えると、スペックルによる効果が最小化されるのは、
縞パターンを照射する際に用いられる定位相表面に対して測定対象表面が平行で
あるか(またはほぼ平行である)場合である。この条件は、表面垂線に対して平
行な2つのソースポイントの方向付けを示すベクトルを備えるかすめ入射線を用
いることにより、達成可能である。特別な場合を除き、表面上のどのポイントP 0 についてもこの条件を満たすことは不可能である。しかし、多くの実用上重要
なオブジェクト(例えば、空気力学面(例えば、車両ボディパネルおよび機体パ
ネル))については、この条件に近い条件を得ることが可能である。
ため、単一のソースヘッドを用いて表面全体上のρ\を最小化するための上記条
件を満たすことは不可能である。このような状況において、複数のソースヘッド
を用いることにより、この条件をより満足させることが可能である。これらの複
数のソースヘッドの配置は、複数のソースヘッドの1つがオブジェクト10の各
主要領域または最重要対象領域に対して応答性を持つように行われる。各ソース
ヘッドからの情報の出どころを区別するために特定の形態の多重化が行われると
仮定する。同じ領域を照射するソースヘッドが1つ以上ある場合、その結果得ら
れる複数の測定値に対し、当該領域における各測定において予測されるスペック
ルノイズのレベルに従って重み付けを行うことができる。複数のソースヘッドを
用いると、シャドーイングによる制約を解消する際にも有用である場合がある。
複数のソースヘッドからの情報は、オーバーラップ領域をビューイングする所与
の検出器22について容易に組み合わせることが可能である。各ピクセルに対応
する角度情報は各ソースヘッドについて同じであるため、登録は自動的に行われ
る。
ピクセル上のρ\の変動を最小化するための条件を満たすような様式で行うこと
ができない。スペックルによる効果の軽減をこのような制約による影響を受けず
に行うさらなる技術は、スペックル平均化に基づく。スペックル平均化は、異な
るパラメータにわずかな差異を設けた(例えば、光学コンフィギュレーションの
変動)さらなる測定を行うことにより、達成することが可能である。あるいは、
露光の間にパラメータを変化させることも可能である。平均化を行うと、変動す
る一連のパラメータが異なるスペックルパターンの異なる様式の統計的実現を可
能にするくらいに十分に変動した場合、スペックルによる効果は実際に低減する
。
る。Pmの変位方向が図1および12の面から外れると、オブジェクト照射の変
動は最小になる。Pmの変動については、式(47)〜(56)および(59)
もおいて明示的に説明している。スペックルによる効果の低減は、例えば、異な
る値のPmについて(x、y、z)座標の個々の測定値を平均化することにより
、達成可能である。
時に変化するため、縞パターンは静止状態のままである。縞を静止状態のままに
しておくための条件は、式(57)においてa/λの比(またはaνの積)が一
定のままであるときに見られる。静止状態の縞を達成するために、フィードバッ
クメカニズム(例えば、縞モニタ)により、ポイントソースP1およびソースP 2 の動きをレーザ40のチューニングに連動させる(またはその反対を行う)。
この縞モニタは、別個の検出器または検出器の線形アレイからなり得、ソースヘ
ッドに組み込まれて、ビームスプリッタによって分離される主要ビームから間隔
が空けられた縞パターンを測定する。
メージスペックルにおいて複数の脱相関(decorrelation)サイク
ルが発生するくらいに十分に大きくする。ソースヘッド中の別個の経路間の経路
長さのオフセットS0を最小化するよう、注意が必要である。そうしないと、縞
パターンの位相は、以下の式による周波数シフトΔνと共に変動する。 Δφ=2π(ΔνS0)/c あるいは、経路−長さオフセットを故意にゼロ以外の値に設定して、2πの増
分の位相Δφにおいて縞パターンを「ストロボ」して、当該縞の位相をフリーズ
することも可能である。こうすると、ストロボのタイミングを変更することによ
り、異なる位相シフトを達成することができる。
数アンビギティーの分解および測定精度の向上を得るたるための一般的な技術で
あり、このような技術の有用性は、上記の用途をはるかに超える点が理解される
。さらに、これらの技術は、一般化された信号(例えば、一次元信号)にも適用
可能であり、シヌソイド変動に限定されない。
も適用されることが理解されるべきである。構築された光パターンを生成する手
段の例を挙げると、コヒーレント放射の干渉もしくは回折、異なる波長で形成さ
れた干渉パターンもしくは回折パターンの「白色光」の重ね合わせ、およびコヒ
ーレント放射もしくは非コヒーレント放射のいずれかを用いたレンズによるオブ
ジェクト上へのパターンのイメージングがある。このようなイメージングアプロ
ーチの一例として、デジタルプロジェクタによって生成された白色光パターンの
作製を、当該パターンが拡張および接触してアコーディオン状の動き(acco
rdion motion)を生成し、側方に平行移動して位相変調を生成する
ように行うことが可能である。このような様式で白色光パターンを表面上に投射
するとスペックルによる効果を無くすことはできるが、このアプローチには特定
の不利点がある。イメージングされた光パターンには、縞周期(period)
の小ささに制約があり得る。投射されたパターンのフィールド深さも考慮する必
要がある。オブジェクト全体にわたって焦点が合わされている斜角から微細縞パ
ターンを投射するのも不可能になり得る。一方、干渉または回折によって生成さ
れた縞は、いつもオブジェクト表面に「焦点が合って」おり、この状態は、当該
縞の観測ポイントへの距離またはサイズに関係なく保たれる。
よって得られる利点が維持され、広帯域の照射または白色光においてこれらの縞
を形成することにより、スペックルによる効果も無くなる。この実施形態は、照
射スペクトル中の異なるコンポーネント波長において形成される干渉縞強度の非
コヒーレントの重ね合わせに基づく。これらの縞のサイズおよび形状が同じであ
る場合、異なる波長に対応する干渉−縞強度は、パターン中の各ポイントにおい
て互いに強化し合う。
が同じ縞を生成する。レンズによって追随される回折格子を用いることにより、
各コンポーネント波長についてもこの目的を同時に達成することができる。回折
角度が小さい場合、格子によって回折される1次ビームの角度偏差は、当該ビー
ムの波長に比例する。レンズは、この角度偏差をソースP1またはP2の変位に
変換する。このソースP1またはP2の変位もλに比例するため、a/λ比を一
定に維持することが可能となる。この実施形態において、ソースP1またはP2 のロケーションは、波長において広く分布する(smeared out)。白
色光の位相シフトは、格子を側方に平行移動させるかまたは例えばPancha
ratnam位相に基づいた白色光位相シフタを用いることにより、達成可能で
ある。
ついては、多くの種類が可能である。1つの構成では、図10中のビーム46お
よび50中にブラッグ格子を配置する。これらの格子(図示せず)およびプリズ
ム92の配置を、非回折ビームに対応する2つのフォーカルスポットが重複し、
1次ビームがプリズム92のベースに向かって偏向するような配置にする。この
構成により、回折ビームに対応するソース間隔aは、光学波長λに比例する。白
色光位相シフトは、ブラッグ格子の1つを側方に平行移動させることにより、達
成される。縞サイズの変動は、格子の周期を変更することにより、達成される。
一実施形態において、ブラッグセルまたは音響光学変調器を用いて、格子周期を
電子的に変更する。この実施形態において、双方のブラッグセルを同じ周波数で
駆動し、1つの駆動信号の位相を他方の信号に対してシフトさせることにより、
位相シフトを達成する(各ビームは同じドップラーシフトを受けるため、これら
の縞は静止状態である)。別の実施形態において、異なる周期のブラッグ格子を
連続的に配置する。さらに別の実施形態において、単一の回折性素子をビーム4
4中に配置する。別の白色光位相シフタ(図示せず)を、ビーム46または50
のうち1つに配置する。当業者に公知の白色光位相シフタのいくつかの実施形態
において、ビームスプリッタ48は偏光ビームスプリッタである。
干渉縞プロジェクタの一実施形態を示す。ソース250は、放射252の平行ビ
ームを実質的に生成する。この平行ビームは、回折格子254に対して実質的に
垂直入射する角度で方向付けられる。回折格子254は、格子周期Dを有する。
例示目的のため、入力ビーム252を、波長コンポーネントλ1、λ2およびλ 3 から構成されるものとして表す。実際は、ビーム252は任意のスペクトルの
構成を有し得る。回折格子254は、ビーム252を複数の回折ビームに分割す
る。これらの複数の回折ビームの回折する順序は、整数mによって表すことが可
能である。例示目的のため、ビーム252周囲に沿った光線のみを図示している
。これらの回折ビームは、垂直入射に関する以下の格子方程式に従って光学軸2
58に対して角度θmで伝播する。
ダm=+1であり回折ビーム262の回折オーダm=−1であるときの回折効率
を最大かつ均等にするようにな設計にする。他の実施形態において、回折格子2
54の設計を、同じオーダ|m|を有する任意の一連の正ビームおよび負ビーム
の回折効率を最大かつ均等にし、なおかつ、他の全てのオーダーに回折するエネ
ルギーを最小化するするようにな設計にするような設計にする。非回折(m=0
)ビーム264が残留した場合、そのようなビームは全て偏向しない状態で回折
格子254を通過し、レンズ266によって集束され、フォーカルスポット26
8上に来る。
に重複する。一実施形態において、フォーカルスポット268は、光学2重スリ
ット272の中央障害物270によって実質的にブロックすることが可能である
。回折ビーム260および262の異なるスペクトル成分λ1、λ2およびλ3 は、レンズ266によってスペクトル領域274および276上に集束される。
スペクトル領域274内のフォーカルスポットと、所与の波長λに対応するスペ
クトル領域276内のフォーカルスポットとの間の距離a(λ)は、波長λに実
質的に比例する。一実施形態において、レンズ266のアパチャストップ278
を用いて、望ましくない高オーダの回折ビームをブロックすることが可能である
。別の実施形態において。光学2重スリット272の不透明領域を用いて、望ま
しくないレンズ266を通過する残留回折オーダのうち望ましくないもの全てを
ブロックすることが可能である。2つのスペクトル領域274および276から
の放射は、伝播して広帯域干渉縞パターン280を形成する際、拡張および重複
する。縞パターン280は、2重スリット272からの代表的距離Rにおいて、
代表的縞周期dを有する。
づき、以下の代入により、これらの方程式を、図17に示す広帯域干渉縞プロジ
ェクタに対する第1のオーダに適用することができる。
縞を生成することも可能である点に留意されたい。狭帯域照射用のレンズ266
によって追随される回折格子254を用いることによって得られる利点の1つと
して、縞周期dは波長に対して感度を持たないため、ソースの周波数ずれは測定
を実質的に劣化させない。例えば、レーザダイオードは、比較的低コストであり
入手も容易なソースであるが、その動作波長は温度依存性である。しかし、この
技術は温度依存性の波長シフトに対して感度を持たないため、測定劣化を招くこ
となくレーザダイオードを用いることが可能である。
相プロファイルを持ち、代表的波長λ2の場合、その相対位相遅延は、0°と1
80°との間で交互に変化し、その際のデューティサイクルは50%である。格
子254は比較的高効率であり、利用可能なエネルギーのおよそ40.5%をm
=−1およびm=+1の回折オーダそれぞれに回折し、名目上は、0%をm=0
および他の均等な回折オーダに回折する。格子254の相対的な位相遅延は波長
の関数であるため、オーダm=0のときの非回折ビーム264中のエネルギーは
、代表的波長λ2と異なる波長の場合、増加する。一実施形態において、格子2
54を比較的肉薄に作製すると、この格子254の波長依存レベルが低下するた
め、この格子254は、広範囲の周波数スペクトルにわたって良好に機能する。
る工程は、回折格子254を図17に示す方向282に平行移動させるだけで達
成される。白色光または広帯域の位相シフトが実現するのは、回折格子254を
格子周期Dの所与の一部分の分だけ平行移動させると、縞パターン280の各ス
ペクトル成分も、当該縞周期dの同じ一部分の分だけシフトするからである。例
えば、格子252をD/4または1/4サイクルだけ平行移動させると、干渉縞
パターン280も1/4サイクル(または90°)だけシフトする。
を達成できる方法には様々な種類がある。一実施形態において、回折角度θmが
小さい場合、格子254の周期Dを2倍にすると、ビーム260および262の
θmの大きさ(これは、式(74)の小さい角度の近似である)が半減し、その
結果、縞パターン280の周期dは2倍になる。別の実施形態において、レンズ
266の焦点距離fを低下させると、縞パターン280の周期dを増加させるこ
とができる。縞パターン280の周期dは、回折格子254の周期Dと、レンズ
266の焦点距離fと、第1のオーダへの伝播距離Rとに関連する。これを式で
表すと、以下のようになる。 d=(RD)/(2f) 縞パターン280のdを連続的に変化させるかまたは連続する任意の一連のd
値を生成するためには、(76)中のパラメータのうち1つ以上を連続的に変化
させることが可能な能力が必要となる。一実施形態において、レンズ266は、
焦点距離fが変化するズームレンズであり得る。別の実施形態において、格子に
力を与えて格子を伸張させることにより、格子254の周期Dを変化させること
が可能である。さらに別の実施形態において、一例として液晶デバイス中の可変
格子モード(VGM)効果を通じて、格子254の周期Dを電子的に変化させる
ことが可能である。さらに別の実施形態において、格子ラインに対して直交する
軸に沿って格子254の周期を変化させることが可能であり、これにより、格子
254を282に対して垂直方向に側方にスライドさせると、格子254の周期
Dを制御することができる。上述した縞合成を用いることにより、実際の周期d
を(f、DまたはRの変動が小さい様態で)大きく変化させることが可能となる
。例えば、縞合成式(64)および(67)を以下のように改変することが可能
である。
進縞分割を通じて縞数のアンビギティーを分解すると適切である場合が多い。一
実施形態において、累進縞分割は、別個の一連の回折格子254またはレンズ2
66を図17中のそれぞれの位置に連続的に移動させることにより、達成される
。別の実施形態において、格子254またはレンズ266は切換え可能であり、
別個の値のDおよびfをとる。切換え可能な格子のいくつかの例を挙げると、液
晶の空間光変調器、パターン化されたアライメントの液晶空間光変調器、電子的
に切換え可能な回折性光学素子、カスケード式の電子的に切換え可能な回折性光
学素子、および超小型電子システム(MEMS)がある。
延器292またはピクセルの線形アレイからなる液晶空間光変調器290であり
、その位相遅延は、周期的な位相−遅延パターン294が生成されるように制御
することが可能である。このようにして、肉薄の位相格子を生成し、電気的に作
動させることができる。例えば、ノミナルの相対位相遅延が50%のデューティ
サイクルで0°および180°に交互に変化する方形波の位相プロファイルを生
成することが可能である。その結果得られた回折格子254の周期Dおよびオフ
セット282は、別個の工程において可変性であり、高い反復性を有する。した
がって、得られた縞パターン280の縞間隔dおよび位相シフトはどちらとも、
別個の反復可能な工程において高速で電気的に制御可能である。
のデューティ−サイクルの肉薄の位相格子を半分生成する際に同じ位相遅延レベ
ルにおいて反復されるピクセル292の数である。一例として、図18において
nblock=2である。格子周期Dは、以下の式によって得られる。
である。所与の値のnblockが以下のようになっている場合、回折格子25
4の最小位相シフトΔφおよびそれに起因する干渉−縞パターン280の最小位
相シフトΔφを生成することが可能である。 Δφ=180°/nblock したがって、例えば、nblock値が2である場合、±90°、±180°
および±270°の位相シフトを生成することができ、nblock値が3であ
る場合、±60°、±120°、±180°、±240°および±300°の位
相シフトを生成することができる。3つの任意の位相シフトφ1、φ2およびφ 3 に関する一般化された形式の式(6)は、以下のように表すことが可能である
。
信号強度である。そのため、位相シフト値を、位相シフト値の利用可能性にした
がって、nblockの各値に応じて個別に選択することができる。
際の格子周期Deffの利用可能な値のうち、異なる値のnblockを用いて
2つの測定値から合成することが可能な値は、以下の式によって得られる。
81)を用いてnblock1およびnblock2を連続する整数にすること
により、得られる。可能な測定戦略の1つの例示として、ピクセル幅wpは10
μmであり、nblock0=2、nblock1=19およびnblock2 =20の場合の3つの測定のシーケンスを行うと仮定する。これらの3つの測定
のシーケンスを用いて、式(78)においてnblock=2を用い、式(78
)においてnblock=20を用いて400μmの第2の格子周期を用いて、
40μmの第1の格子周期Dを生成することが可能である。ずっと大きな第3の
格子周期は4.2mmであり、これは、式(81)中のnblock1=19お
よびnblock2=20を用いることにより合成することが可能である。この
シーケンス中の連続する格子周期間の比はおよそ10であり、この値は、累進縞
分割において適切なシーケンスである。
、スペクトル領域274および276の広がりおよびロケーションも広範囲に変
化する。その場合、一実施形態において、より高度なバージョンの2重スリット
272(例えば、図19に示すビームブロック296)を用いて、不要な残留回
折オーダを無くすことができる。ビームブロック296は、2つのネスト状の2
重スリットを含み、中央障害物268を含む。上記の例における測定シーケンス
において用いられる値nblock0が最小である場合、外側の2重スリット2
75はスペクトル領域274および276を通過し、一方、値nblock1お
よびnblock2の値が大きい場合、内側2重スリット277は所望のスペク
トル領域274’および276’を通過する。
配置する際、スペクトル領域274および276のイメージ302および304
がイメージ面306に形成されるような配置にすることにより、アコーディオン
状の動きを達成する。別の実施形態において、レンズ300の焦点距離が負であ
る場合、イメージ302および304は仮想イメージ(図示せず)であり得る点
に留意されたい。一実施形態において、レンズ300を光学軸258に沿って移
動させることによってイメージ302および304の倍率を変更することにより
、アコーディオン状の動きを得る。または、別の実施形態において、レンズ30
0の焦点距離を変更することにより、アコーディオン状の動きを得る。さらに別
の実施形態において、レンズ300はズームレンズであり、イメージ面306の
ロケーションを変化させることなく焦点距離を変化させる。さらに別の実施形態
において、焦点距離が異なる別個の一連のレンズを連続的に適切な位置に配置す
ることが可能である。レンズ300を用いると、図17内のコンポーネントのう
ち任意のコンポーネント(例えば、2重スリット272)を変化させることなく
、アコーディオン状の動きが生成される。
形態は、広帯域ソース250と、音響光学変調器(AOM)306と、レンズ2
66とを含む。AOM306は、乗算器310によって生成された合成信号30
8によって駆動される。この乗算器310は、可変周波数fmを有するシヌソイ
ド312を、一定周波数fcを有するシヌソイド314で乗算する。合成信号3
08は、周期が異なる2つの進行する音波をAOM306中に生成し、すると、
AOM306は角度が離れた2つのビーム260および262を生成する。異な
るスペクトル成分の回折ビーム260および262が、レンズ266によってス
ペクトル領域274および276に集束し、ここで、スペクトル領域274内の
フォーカルスポットと、所与の波長λに対応するスペクトル領域276内のフォ
ーカルスポットとの間の間隔a(λ)は、波長λに実質的に比例する。さらに、
各スペクトル成分について、スペクトル領域274および276のスペクトル成
分間にあるような中間点を、実質的に位置268に局所化させる。fmを調節す
ることにより、各波長コンポーネントλの間隔a(λ)を、対称点268の周囲
において、ゼロの間隔からAOM306の帯域が許す最大間隔まで調節する。2
つのスペクトル領域274および276からの放射は、伝播して代表的周期dの
広帯域干渉−縞パターン280を代表的距離Rにおいて形成すると、拡張および
重複する。fmを変化させると、縞パターン280は拡張および収縮する。
、互いにドップラーシフトする。その結果、縞パターン280は、合成信号30
8の2つの周波数成分間の周波数差に比例した速度で重複領域を移動する。干渉
パターン280は、振幅変調ソース250により、特定の空間位相において2つ
のビーム260および262の周波数差2fmにおいて効果的に凍結することが
可能である。縞280の空間位相は正確に制御することが可能であり、この制御
は、レーザ駆動信号316のAOM制御信号308に対する位相を変化させる(
例えば、遅延を変化させる)ことにより、行われる。あるいは、駆動信号316
を用いて、光路(図示せず)内に配置された外部の振幅変調器を制御するか、ま
たは、縞パターンを観測する際に用いられる検出器22の感度を変調することも
可能である。
的に平行であるため、回折性素子254(図17)または306(図21)は、
所与の回折オーダmの各スペクトル成分を同じ回折角度θmに回折させる。この
コンフィギュレーションを他の様態に改変することも可能であり、このような改
変は、本発明の趣旨および範囲の内におさまる。例えば、別の実施形態において
、ビーム252は発散型または収束型であり得る。さらに別の実施形態において
、ビーム252は、回折性素子254に斜角で入射し得る。さらに別の実施形態
において、屈折素子および回折性素子のオーダを逆にすることもできる。
光ソース250を用いてもよい。時間コヒーレンスが低い場合(広範囲のスペク
トルの内容の場合)、スペックルによる効果は低減するかまたは無くなる。一方
、空間コヒーレンスが高いと、スペクトル領域274および276中の各スペク
トル成分について、より厳密なフォーカルスポットを生成することが可能になる
。一般的には、フォーカルスポットが厳密であると、縞の可視性が向上し、縞パ
ターン280のフィールド深さが大きくなる。空間的にコヒーレントなソースの
場合、そのフィールド深さは無限である点に留意されたい。さらに、フォーカル
スポットが厳密であると、異なる回折オーダmに対応するスペクトル領域を分離
するのも容易になるため、望ましくないスペクトル領域を代表的マスク272ま
たは296によってブロックすることが可能である。
ザベースにすることである。一実施形態において、ソース250は、平行な光学
部品を備えるレーザダイオードである。一般的には、スペクトルの幅がおよそ3
nmであるレーザダイオードを利用することが可能である。別の実施形態におい
て、異なる波長で動作するレーザダイオードのアレイを単一のビームに結合する
ことにより、このスペクトル幅を広げることが可能である。さらに別の実施形態
において、高速スキャンのためのチューナブルレーザを用いることも可能である
。レーザのスペクトルを広げるための別のアプローチでは、標準的な単一の−モ
ードシリカファイバを受動型でQスイッチ型のミクロチップレーザによってポン
ピングする。ファイバ中の様々な非線形効果により、ファイバの出力端部におい
て極めて広帯域のスペクトルが生成される。さらに別のアプローチにおいて、モ
ードロックチタン、サファイアレーザによって広帯域照射を達成する。
は他にも多くある。例えばアーク灯を光学ファイバに結合させることにより、低
い(が適切な)空間コヒーレンスを達成することが可能である。ファイバの出力
端部は、長大な(extended)ソースのように機能する。レンズを用いて
ファイバの出力を平行にすると、当該レンズの焦点距離が長くなるほど、ビーム
252の伝播角度の変動は小さくなり、フォーカルスポットも厳密になる。長大
なソースを用いると、レンズ266を光学軸258に沿った距離に配置する際、
回折格子254の表面またはAOM306の活性エリアが照射されているオブジ
ェクト10の表面または表面近隣上にイメージングされるように配置を行うこと
により、空間の非コヒーレンスによる効果が最小化される。このコンフィギュレ
ーションを用いると、縞パターン280の所与のエリアに到達した光線が、入力
ビームを異なる回折オーダに分割させる回折性素子の同一領域から発生するもの
全てと干渉するため、空間の非コヒーレンスによる効果が最小化される。
づくかに関係無く)変化し得る。例えば、圧電変換器によって光コンポーネント
(例えば、図10中の鏡94もしくは96または図17中の格子254)を移動
させることにより位相シフトを行うと、ヒステリシスによって矛盾が生じる。こ
れらの位相シフトがその予想値から変化すると、式(6)または式(80)にお
いてφを計算する際にエラーが発生する。これらのφ値のエラーは周期的な性質
を持ち、再構築された表面プロファイルがリップルが生じた外観になる原因とな
る。また、リップルが生じた外観の原因としては、不完全な照射シーケンスもあ
る。これらの不完全な照射シーケンスは、例えば信号強度I1、I2およびI3 の取得期間の間の露光レベルまたはバイアスレベルの変化によって生じ得る。露
光レベルの変換の原因としては、例えば、ソース強度の変動または完全なシャッ
ターによる露光時間差がある。
ターンと同じ方向に発生する。リップルパターンの主要周期は、エラーソースに
依存する。位相シフトエラーがあると、例えば、縞周期について2つのリップル
振動が発生する。リップルは、表面プロファイル測定において主要なエラーソー
スとなり得、特に、広帯域照射を通じてスペックルが無くなった場合または狭帯
域照射からのスペックル効果が上述した軽減技術の1つを通じて抑制された場合
、その傾向は強くなる。
(derippling)と呼ばれる汎用性がありかつロバストな技術を紹介す
る。この技術において、データを解析して、実際の位相オフセットおよび露光レ
ベルの変動を判定し、相関値を用いてφを計算する。この技術を説明するため、
まず、式(80)を書き換えて、推測値および推測値から計算される数量を実際
の値と区別するような様式にする。
を、これらの推測値を用いて計算した位相値とする。
る変化を式(82)に取り入れることができる。
素を表す。パラメータAiおよびBiにより、照射シーケンスの露光レベルおよ
びバイアスレベルをそれぞれ変化させることができる。式(83)において、ノ
ミナルの山−谷振幅の幅(swing)が評価対象ユニットであると一般性を失
うことなく仮定することができる。
る場合、不完全な照射シーケンスを、計算された強度により修正して、
*φ2、*φ2および*φ3ではなく実際の値φ1、φ2およびφ3を判定して
式(82)において用いることを可能にすると、位相ステップエラーを無くすこ
とが可能となる。リップルを最小化するための1つのアプローチとして、*φに
おけるリップルの程度を示す測定基準(metric)が最小化されるようにパ
ラメータAi、Biおよびφiを変化させるアプローチがある。測定基準を最小
化する場合、演算要求が多大になり、多次元のサーチアルゴリズムが必要となる
だけではなく、代表的な表面部分にわたって位相関数*φを最初にアンラップし
なければ、最小化のための適切な測定基準を得るのは困難である。アンラップが
行われた後に利用することが可能な測定基準の一例としては、*φと平滑化され
た*φの値との間の差の標準偏差がある。リップルの低減による寄与のため、こ
の測定基準は低減する傾向となる。
ータの修正値を直接生成するため、上記のリップル消去アプローチよりも優れて
いる。ここで、図22を参照して、この技術を引き続き紹介する。図22は、エ
ラー(φ−*φ)/(2π)のプロットを、異なる値のパラメータAi、Biお
よびφiについて計算された数量*φ/(2π)の関数として含む。このエラー
は、縞数エラーΔNのリップルによって生じる(ripple−induced
)成分に相当し、式(32)に示すようなレンジ分解能を限定する。これらの曲
線において、推定される位相ステップ値は、*φ1=−120°、*φ2=0°
および*φ3=120°である。φおよび*φの両方を、−π〜+πのラップさ
れた数量であるとみなし、これにより、φ/(2π)および*φ/(2π)レン
ジが−1/2〜+1/2となるようにする。曲線(a)は、振幅変化A1=0.
9、Α2=1、A3=1.2の効果を示し、曲線(b)は、バイアス変化B1=
0.2、B2=0、B3=0.1に相当し、曲線(c)は、推測値と異なる実際
の位相シフトφ1=−105°、φ2=0°およびφ3=100°による効果を
示す。最後に、図22aは、図22中の曲線の傾斜をプロットしたものであり、
これは、*φに関するφの導関数およびによって得られ、φ’(*φ)として示
される。この傾斜は、1の周囲で変動し、推定パラメータがその真なる値に近づ
くと、全ての点において1に収束する。ΔNにおけるエラーが比較的小さい場合
、比較的大きな傾斜の変化が生成されるため、傾斜により、Ai、Biおよびφ i の推測値中のエラーが感度良く示される点に留意されたい。
とは明らかである。例えば、これらのプロットは、バイアスが変化すると、縞あ
たり1つのリップルが発生し、位相ステップエラーが発生すると、縞あたり2つ
のリップルが発生することを示す。振幅が変化すると、縞あたり1つおよび2つ
のリップルの混合物が発生する。図22または図22aにおいてグラフ化されて
いる曲線をデータ解析から生成することが可能であれば、Ai、Biおよびφi 上のこれらの曲線を支配する方程式への依存を用いて、これらのパラメータに関
する情報を回復させることができる。その後、これらの修正されたパラメータを
用いてφを計算することにより、データをリップル消去することができる。この
アプローチの難点は、これらのグラフでは、φの真の値が分かっていることを前
提としている点である。
。このアプローチでは、先ず、測定(または特定の測定部分)における*φの各
値の発生を相対的に表すヒストグラムを作成する。このヒストグラムは、*φの
値の可能な範囲を有限数のビンに分割し、各ビン中に収められた*φが発生する
回数をカウントすることにより、生成される。*φのヒストグラムを作成すると
、*φの確率密度関数の近似値が得られ、これはP(*φ)と書かれる。以下の
周知の関係を通じて、2つの確率密度関数P(*φ)およびP(φ)を関連付け
ることが可能である。
を、図22a中にグラフとして示す。確率分布関数P(*φ)は、*φのヒスト
グラムの計算を通じたデータから利用することが可能である。式(85)中のP
(φ)は不明であるが、多くの場合において、φはその可能な値の範囲にわたっ
て均等に分布していると仮定すると妥当である。サンプリング対象となる表面エ
リアを覆う縞の数が多いほど、この仮定の精度も向上する。実際にも、P(φ)
が均等に分布しているというこの仮定は、縞が比較的少数である場合にも適切で
あることが分かっている。しかし、P(φ)の関数の形態が表面に関するさらな
る情報に基づいて正確になっている場合、ヒストグラムベースのリップル消去に
おいて他の形態のP(φ)を用いることも可能である。
*φのヒストグラムを計算することによって得られる)と、導関数φ’(*φ)
の絶対値とを関連付けることができる。その後、この導関数の関数の形態に基づ
いて、Ai、Biおよびφiの実際の値を判定することができる。これを行う方
法を例示するため、ここで、φ’(*φ)の関数の形態をより詳細に解析する。
ステップエラーのみがある場合を想定する。この場合、φ’(*φ)は、以下の
ような関数の形態をとることが分かる。
る)は他の2つの定数KC2およびKS2から判定することが可能であり、その
結果、式(86)における自由度は2だけとなる。これらの定数KC2およびK s2 は、例えば、このヒストグラムの逆数をとり、フーリエ級数展開のcos(
2*φ)項およびsin(2*φ)項に対応するフーリエ係数を計算することに
より、*φのヒストグラムの形状から得ることができる。その後、フーリエ級数
のDC項を用いてcos(2*φ)項およびsin(2*φ)項のフーリエ係数
を正規化して、Kc2およびKs2をそれぞれ得る。
的位相シフトφ1−φ2、φ2−φ3およびφ3−φ1を回復することが可能で
あり、これらの相対的位相シフトのうち2つのにより、第3の値を判定する。こ
れらの相対的位相シフトがわかっていれば、リップルを完全に除去するには十分
である。しかし、得られた*φの向上した計算値には、わずかなオフセットエラ
ーが存在する。このオフセットエラーは、φ2の真の値が分からないことに起因
する。
あり、ゼロの周囲に集中していると仮定する。そうすると、これらの位相シフト
を*φ1=−φstep、*φ2=0、および*φ3=φstepと表すことが
できる。デフォルトでφ2=0であるとも仮定する。次いで、以下の関係により
、Kc2およびKs2からφ1およびφ3の値を判定することができる。
(88)は以下のように整理される。
エラーおよび振幅エラーの組み合わせにこの技術を適用する。この場合の数式φ
’(*φ)は式(86)よりもずっと複雑であるが、この数式は、Ai−1およ
びφi−*φiの値が小さい場合、以下のように近似化することができる。
得ることが可能であり、これらの係数は、φ’(*φ)を近似化する際に用いら
れるヒストグラムの逆数のフーリエ級数展開の対応するフーリエ係数である。
テップサイズを不均等にすると、以下の関係を通じて、Kc1、Ks1、Kc2 、およびKs2からφ1、φ3、A1およびA3の新規近似値を判定することが
できる。
づくため、φ1、φ3、A1およびA3の計算値は向上するが、正確ではない。
所望であれば、このプロシージャを繰り返すことにより、精度を高くすることが
可能である。そうするためには、式(92)および(93)から計算されたφ1 およびφ3の値は、式(82)において新規の推測値*φ1および*φ3となり
、式(94)および(95)から得られたA1およびA3の値は、式(84)に
おいて新規の推測値*A1および*A3となる。このプロシージャは高速で収束
するため、数回の反復だけで小数位のレベルでの高い精度が得られる。
化Biを自動処理する。反復回数が多いほど、バイアス変化Biは、他のパラメ
ータAiおよびφiに折り畳まれて、その効果は無くなっていく。このアプロー
チを用いると、リップルが無くなり、φ2が正確には分かっていない場合に生じ
るエラーに類似する微小オフセットエラーが発生する。
プトを取り入れた他の実施形態を用いることも可能であり、また、本発明の範囲
および趣旨内において多くの改変例が可能であることが当業者にとって明らかで
ある。したがって、これらの実施形態は、開示された実施形態に限定されるので
はなく、本明細書内の特許請求の範囲の趣旨および範囲のみによって限定される
べきであることが認識される。
ック図である。
ブロック図である。
形態のブロック図である。
の一実施形態のブロック図である。
れる工程の一実施形態のフローチャートである。
サ装置の一実施形態のブロック図である。
含む検出器およびプロセッサ装置の別の実施形態のブロック図である。
ッサ装置の別の実施形態のブロック図である。
の素子の一実施形態を示す。
出器の素子の別の実施形態を示す。
形態のブロック図である。
形態のブロック図である。
示す。
フローチャートである。
実施形態のフローチャートである。
施形態のフローチャートである。
。
形態を示す。
一実施形態を示す。
公知の*φ/(2π)の関数として縞数誤差ΔNをグラフィカルに示す表示を含
む。
Claims (120)
- 【請求項1】 表面を有するオブジェクト上において、該オブジェクトの表
面上のポイントの三次元位置情報を判定する方法であって、 a)スペクトル領域の空間分布を有する2つの放射ソースを提供する工程と、 b)該ソースそれぞれからの放射で該表面を照射して、該表面上の第1の位置
に第1の縞パターンを生成する工程と、 c)該第1の縞パターンを第2の位置に移動させる工程と、 d)該第1の縞パターンの第1の位置および第2の位置に応答して、第1のラ
ップされたサイクルマップを生成する工程と、 e)該第1の縞パターン中の縞数を推定する工程と、 f)該第1の縞パターンを変更して、第1の位置において第2の縞パターンを
生成する工程と、 g)該第2の縞パターンを第2の位置に移動させる工程と、 h)該第2の縞パターンの第1の位置および第2の位置に応答して、第2のラ
ップされたサイクルマップを生成する工程と、 i)該第1の縞パターン中の推定された縞数に応答して、該第2の縞パターン
中の縞数を推定する工程と、 j)該第2の縞パターンおよび該第2のラップされたサイクルマップ中の推定
された縞数に応答して、該ポイントの該表面上での三次元位置情報を判定する工
程と、 を包含する、方法。 - 【請求項2】 前記2つの放射ソースは互いにコヒーレントである、請求項
1に記載の方法。 - 【請求項3】 前記2つの放射ソースのスペクトル分布は狭帯域である、請
求項1に記載の方法。 - 【請求項4】 前記2つの放射ソースを提供する工程は、単一の放射ソース
からの放射を分割する工程を包含する、請求項1に記載の方法。 - 【請求項5】 前記2つの放射ソースを提供する工程は、2つのレーザ放射
ソースを提供する工程を包含する、請求項1に記載の方法。 - 【請求項6】 前記工程第1の縞パターンおよび第2の縞パターン中の縞数
を推定する工程は、数学的モデルを適用する工程を包含する、請求項1に記載の
方法。 - 【請求項7】 前記第1の縞パターンおよび第2の縞パターン中の縞数を推
定する工程は、線形の推定モデルを適用する工程を包含する、請求項1に記載の
方法。 - 【請求項8】 前記第1の縞パターンおよび第2の縞パターン中の縞数を推
定する工程は、異なる表面について判定された三次元位置情報を適用する工程を
包含する、請求項1に記載の方法。 - 【請求項9】 前記ソースの1つのスペクトル領域はそれぞれ、該ソースの
残りの各スペクトル領域から、該スペクトル領域の各波長に比例する各距離の分
だけ隔離される、請求項1に記載の方法。 - 【請求項10】 前記距離は、前記スペクトル領域の波長に線形比例する、
請求項9に記載の方法。 - 【請求項11】 前記第1の縞パターンを変更する工程は、前記スペクトル
領域それぞれを、該スペクトル領域の残りに対してある距離の分だけ隔離する工
程を包含する、請求項1に記載の方法。 - 【請求項12】 前記距離は、2つのスペクトル領域の各々から等距離であ
る中間点を含み、該スペクトルの領域それぞれを隔離する工程は、該距離を変更
して該中間点を固定状態にする工程を包含する、請求項11に記載の方法。 - 【請求項13】 前記第2の縞パターンを移動させる工程は、前記ソースの
1つからの前記空間分布中のスペクトル領域の位相を、該ソースの残りからの空
間分布中の各スペクトル領域の位相に対して変化させる工程を包含する、請求項
1に記載の方法。 - 【請求項14】 表面上の所定の間隔の縞パターンに対応するラップされた
サイクルマップを合成する方法であって、 a)スペクトル領域の空間分布を有する2つの放射ソースを提供する工程と、 b)該ソースそれぞれからの放射で該表面を照射して、該表面上の第1の位置
に第1の縞パターンを生成する工程と、 c)該第1の縞パターンを第2の位置に移動させる工程と、 d)該第1の縞パターンの第1の位置および第2の位置に応答して、第1のラ
ップされたサイクルマップを生成する工程と、 e)該第1の縞パターンを変更して、第1の位置において第2の縞パターンを
生成する工程と、 f)該第2の縞パターンを第2の位置に移動させる工程と、 g)該第2の縞パターンの第1の位置および第2の位置に応答して、第2のラ
ップされたサイクルマップを生成する工程と、 h)該第1のラップされたサイクルマップから該第2のラップされたサイクル
マップを減算する工程と、 i)該第2のラップされたサイクルマップと該第1のラップされたサイクルマ
ップとの間の差をラップして、該所定の間隔の縞パターンに対応するラップされ
たサイクルマップを生成する工程と、 を包含する、方法。 - 【請求項15】 前記2つの放射ソースは互いにコヒーレントである、請求
項14に記載の方法。 - 【請求項16】 前記2つの放射ソースのスペクトル分布は狭帯域である、
請求項14に記載の方法。 - 【請求項17】 前記2つの放射ソースを提供する工程は、単一の放射ソー
スからの放射を分割する工程を包含する、請求項14に記載の方法。 - 【請求項18】 前記2つの放射ソースを提供する工程は、2つのレーザ放
射ソースを提供する工程を包含する、請求項14に記載の方法。 - 【請求項19】 前記ソースの1つのスペクトル領域はそれぞれ、該ソース
の残りの各スペクトル領域から、該スペクトル領域の各波長に比例する距離の分
だけ隔離される、請求項14に記載の方法。 - 【請求項20】 前記距離は、前記スペクトル領域の波長に線形比例する、
請求項19に記載の方法。 - 【請求項21】 前記第1の縞パターンを変更する工程は、前記スペクトル
領域それぞれを、該スペクトル領域の残りに対してある距離の分だけ隔離する工
程を包含する、請求項14に記載の方法。 - 【請求項22】 前記距離は、2つのスペクトル領域の各々から等距離であ
る中間点を含み、該スペクトルの領域それぞれを隔離する工程は、該距離を変更
して該中間点を固定状態にする工程を包含する、請求項21に記載の方法。 - 【請求項23】 前記第2の縞パターンを移動させる工程は、前記ソースの
1つからの前記空間分布中のスペクトル領域の位相を、該ソースの残りからの空
間分布中の各スペクトル領域の位相に対して変化させる工程を包含する、請求項
14に記載の方法。 - 【請求項24】 前記2つの放射ソースは第1の周波数を有し、前記第1の
縞パターンを変更して第2の縞パターンを生成する工程は、該第1の周波数を第
2の周波数に変更する工程を包含する、請求項14に記載の方法。 - 【請求項25】 表面を有するオブジェクト上において、該オブジェクトの
表面上のポイントの三次元位置情報を判定する方法であって、 a)第1の周波数を有しかつ互いにコヒーレントな第1の放射ビームおよび第
2の放射ビームを生成する工程と、 b)該表面を該第2の放射ビームで照射する工程と、 c)第1の放射ビームと、該第2の放射ビームからの放射とが該表面によって
散乱するのに応答して、第1の位置において第1の干渉パターンを生成する工程
と、 d)該第1の干渉パターンを第2の位置に移動させる工程と、 e)該第1の干渉パターンの第1の位置および第2の位置に応答して、第1の
ラップされたサイクルマップを生成する工程と、 f)該第1の干渉パターン中の強度サイクルを推定する工程と、 g)該第1の干渉パターンを変更して、第2の干渉パターンを第1の位置にお
いて生成する工程と、 h)該第2の干渉パターンを第2の位置に移動させる工程と、 i)該第2の干渉パターンの第1の位置および第2の位置に応答して、第2の
ラップされたサイクルマップを生成する工程と、 j)該第1の干渉パターン中の推定された強度サイクルに応答して、該第2の
干渉パターン中の強度サイクルを推定する工程と、 k)該第2の干渉パターンおよび該第2のラップされたサイクルマップ中の推
定された強度サイクルに応答して、該三次元位置情報を計算する工程と、を包含
する、方法。 - 【請求項26】 前記第1の干渉パターンを変更して第2の干渉パターンを
生成する工程は、前記第1の周波数を第2の周波数に変更する工程を包含する、
請求項25に記載の方法。 - 【請求項27】 前記表面をイメージ面上にイメージングする工程をさらに
包含する、請求項25に記載の方法。 - 【請求項28】 前記イメージ面は検出器を含む、請求項27に記載の方法
。 - 【請求項29】 表面を有するオブジェクト上において、該オブジェクトの
表面上のポイントの三次元位置情報を判定する方法であって、 a)互いにコヒーレントな第1の放射ビームおよび第2の放射ビームを生成す
る工程と、 b)該表面を該第2の放射ビームで照射する工程と、 c)該第1の放射ビームと、該第2の放射ビームからの放射とが該表面によっ
て散乱するのに応答して、第1の位置において第1の干渉パターンを生成する工
程と、 d)該第1の放射ビームおよび第2の放射ビームのうちの1つを、該第1の放
射ビームおよび第2の放射ビームのもう一方に対して位相シフトさせて、第1の
位相シフトした干渉パターンを生成する工程と、 e)該第1の干渉パターンおよび該第1の位相シフトした干渉パターンに応答
して、第1のラップされたサイクルマップを生成する工程と、 f)該第1の干渉パターンを変更して、第2の干渉パターンを生成する工程と
、 g)該第1の放射ビームおよび第2の放射ビームのうちの1つを、該第1の放
射ビームおよび第2の放射ビームのもう一方に対して位相シフトさせて、第2の
位相シフトした干渉パターンを生成する工程と、 h)該第2の干渉パターンおよび該第2の位相シフトした干渉パターンに応答
して、第2のラップされたサイクルマップを生成する工程と、 i)該第2のラップされたサイクルマップを該第1のラップされたサイクルマ
ップから減算する工程と、 j)該第2のラップされたサイクルマップと、該第1のラップされたサイクルマ
ップとの間の差をラップして、ラップされたサイクルマップを生成する工程と、 を包含する、方法。 - 【請求項30】 前記第1の放射ビームおよび第2の放射ビームは第1の周
波数を有し、前記第1の干渉パターンを変更して第2の干渉パターンを生成する
工程は、該第1の周波数を第2の周波数に変更する工程を包含する、請求項29
に記載の方法。 - 【請求項31】 前記ラップされたサイクルマップを用いることによって三
次元位置情報を判定する工程をさらに包含する、請求項29に記載の方法。 - 【請求項32】 前記表面をイメージ面上にイメージングする工程をさらに
包含する、請求項29に記載の方法。 - 【請求項33】 前記イメージ面は検出器を含む、請求項32に記載の方法
。 - 【請求項34】 縞をオブジェクト表面上に投射する装置であって、 a)スペクトル分布を有する2つの放射ソースと、 b)該2つのソースと光学的に連絡するコリメータであって、2つの広帯域放
射の実質的に平行なビームを生成する、コリメータと、 c)該コリメータと光学的に連絡する回折性格子と、 d)該回折性格子と光学的に連絡するレンズであって、スペクトル領域の空間
分布を有する2つの放射イメージを生成するレンズと、 を備える、装置。 - 【請求項35】 前記ソースの1つのスペクトル領域はそれぞれ、該ソース
の残りの各スペクトル領域から、該スペクトル領域の各波長に比例する各距離の
分だけ隔離される、請求項34に記載の装置。 - 【請求項36】 前記距離は、前記スペクトル領域の波長に線形比例する、
請求項35に記載の装置。 - 【請求項37】 前記距離は、2つのスペクトル領域の各々から等距離であ
る中間点を含み、該中間点は固定される、請求項35に記載の装置。 - 【請求項38】 前記2つの放射ソースは互いにコヒーレントである、請求
項34に記載の装置。 - 【請求項39】 前記2つの放射ソースのスペクトル分布は狭帯域である、
請求項34に記載の装置。 - 【請求項40】 前記オブジェクトの表面上のポイントの三次元位置情報を
判定する検出器をさらに備える、請求項34に記載の装置。 - 【請求項41】 前記2つの放射ソースは単一の放射ソースから生成される
、請求項34に記載の装置。 - 【請求項42】 前記回折性格子に結合されたトランスレータをさらに備え
、該トランスレータは、前記スペクトル領域の1つの該スペクトル領域の残りに
対する相対的位相をシフトさせる、請求項34に記載の装置。 - 【請求項43】 オーダをブロックするマスクをさらに備える、請求項34
に記載の縞をオブジェクト表面上に投射する装置。 - 【請求項44】 オーダをブロックする2重スリットのマスクをさらに備え
る、請求項34に記載の縞をオブジェクト表面上に投射する装置。 - 【請求項45】 縞をオブジェクト表面上に投射する方法であって、 a)ある距離だけ隔離された2つの放射ソースを提供する工程であって、該ソ
ースはそれぞれ、スペクトル分布を有する放射を生成する、工程と、 b)該放射を平行にして、2つの実質的に平行な放射ビームを生成する工程と、 c)該平行な放射ビームのスペクトル成分を描写する工程と、 d)スペクトル成分の空間分布を有する2つのイメージ放射を生成する工程と、 を包含する、方法。 - 【請求項46】 前記空間分布の1つのスペクトル成分はそれぞれ、該空間
分布の残りの各スペクトル成分から、該スペクトル成分の各波長に比例する距離
だけ隔離される、 請求項45に記載の方法。 - 【請求項47】 前記距離は、前記スペクトル成分の波長に線形比例する、
請求項46に記載の方法。 - 【請求項48】 前記距離は、2つのスペクトル成分の各々から等距離であ
る中間点を含み、該中間点は固定される、請求項46に記載の方法。 - 【請求項49】 前記2つの放射ソースは互いにコヒーレントである、請求
項45に記載の方法。 - 【請求項50】 前記2つの放射ソースのスペクトル分布は狭帯域である、
請求項45に記載の方法。 - 【請求項51】 前記オブジェクトの表面上のポイントの三次元位置情報を
判定する工程をさらに包含する、請求項45に記載の方法。 - 【請求項52】 前記2つの放射ソースを提供する工程は、単一の放射ソー
スからの放射を分割する工程を包含する、請求項45に記載の方法。 - 【請求項53】 前記スペクトル領域の1つの該スペクトル領域の残りに
対する相対的位相をシフトさせる工程をさらに包含する、請求項45に記載の方
法。 - 【請求項54】 表面を有するオブジェクト上において、該オブジェクトの
表面上のポイントの三次元位置情報を判定する方法であって、 a)第1の周期的パターンを投射して、該表面上の第1の位置において第1の
投射パターンを生成する工程と、 b)該第1の投射パターンを第2の位置に移動させる工程と、 c)該第1の投射パターンの第1の位置および第2の位置に応答して、第1の
ラップされたサイクルマップを生成する工程と、 d)該第1の投射パターン中のサイクル数を推定する工程と、 e)該第1の投射パターンを変更して、該第1の位置において第2の投射パタ
ーンを生成する工程と、 f)該第2の投射パターンを該第2の位置に移動させる工程と、 g)該第2の投射パターンの第1の位置および第2の位置に応答して、第2の
ラップされたサイクルマップを生成する工程と、 h)該第1の投射パターン中の推定されたサイクル数に応答して、該第2の投
射パターン中のサイクル数を推定する工程と、 i)該第2の投射パターン中の推定されたサイクル数および該第2のラップさ
れたサイクルマップに応答して、該表面を判定する工程と、 を包含する、方法。 - 【請求項55】 前記第1の投射パターンを変更する工程は、第2の周期的
パターンを投射して、第2の投射パターンを第1の位置において生成する工程を
包含する、請求項54に記載の方法。 - 【請求項56】 表面上に投射された所定の間隔の周期的パターンに対応す
るラップされたサイクルマップを合成する方法であって、 a)第1の周期的パターンを投射して、該表面上の第1の位置において第1の
投射パターンを生成する工程と、 b)該第1の投射パターンを、該表面上の第2の位置に移動させる工程と、 c)該第1の投射パターンの第1の位置および第2の位置に応答して、第1の
ラップされたサイクルマップを生成する工程と、 d)該第1の投射パターンを変更して、第2の投射パターンを第1の位置にお
いて生成する工程と、 e)該第2の投射パターンを該第2の位置に移動させる工程と、 f)該第2の投射パターンの第1の位置および第2の位置に応答して、第2の
ラップされたサイクルマップを生成する工程と、 g)該第1のラップされたサイクルマップから該第2のラップされたサイクル
マップを減算する工程と、 h)該第2のラップされたサイクルマップと該第1のラップされたサイクルマ
ップとの間の差をラップして、該所定の間隔の縞パターンに対応するラップされ
たサイクルマップを生成する工程と、 を包含する、方法。 - 【請求項57】 前記第1の投射パターンを変更する工程は、第2の周期的
パターンを投射して、第2の投射パターンを第1の位置において生成する工程を
包含する、請求項56に記載の方法。 - 【請求項58】 スペックルがオブジェクト表面上の測定ポイントに与える
影響を軽減する方法であって、 a)コヒーレントな縞パターンを生成する工程と、 b)該コヒーレントな縞パターンが該オブジェクト表面を実質的にかすめるよ
うに、該縞パターンを光路に沿って該オブジェクト表面に投射する工程と、 c)該オブジェクト表面のイメージ中の該縞パターンおよび該スペックルを検
出する工程であって、該オブジェクト表面に対する法線は該光学経路に実質的に
直交する、工程と、 を包含する、方法。 - 【請求項59】 前記コヒーレントな縞パターンは、前記オブジェクト表面
に対して0〜45°の角度で該オブジェクト表面を実質的にかすめる、請求項5
8に記載の方法。 - 【請求項60】 前記2つの放射ソースを提供する工程をさらに包含する、
請求項58に記載の方法。 - 【請求項61】 前記2つのソースは、前記オブジェクト表面の上側に配置
される、請求項60に記載の方法。 - 【請求項62】 前記2つのソースは、前記オブジェクト表面に対する法線
に沿って実質的に垂直である、請求項60に記載の方法。 - 【請求項63】 前記縞パターンはソースによって生成される、請求項58
に記載の方法。 - 【請求項64】 前記2つのソースは互いにコヒーレントである、請求項6
3に記載の方法。 - 【請求項65】 前記2つのソースはレーザソースである、請求項63に記
載の方法。 - 【請求項66】 前記2つのソースは、単一のソースを分割することによっ
て生成される、請求項63に記載の方法。 - 【請求項67】 スペックルがオブジェクト表面上の測定ポイントに与える
影響を軽減する方法であって、 a)ある距離だけ隔離された2つの放射ソースからコヒーレントな縞パターン
を生成する工程と、 b)該コヒーレントな縞パターンを該オブジェクト表面上に投射する工程と、 c)該オブジェクト表面のイメージ内において該縞パターンおよび該スペック
ルを検出する工程と、 d)該縞パターンが実質的に静止状態でありかつ該スペックルが変化するよう
に、該2つのソースを平行移動させる工程と、 e)該オブジェクト表面のイメージにおいて、該新規の縞パターンと、該変化
したスペックルとを検出する工程と、 f)該検出された縞パターンおよび該スペックルの変化に応答して、実質的に
スペックルの発生が無い状態で該縞パターンを判定する工程と、 を包含する、方法。 - 【請求項68】 工程d)および工程e)を繰り返し反復する工程をさらに
包含する、請求項67に記載の方法。 - 【請求項69】 前記2つの放射ソースはレーザ放射である、請求項67に
記載の方法。 - 【請求項70】 スペックルがオブジェクト表面上の測定ポイントに与える
影響を軽減する方法であって、 a)ある距離だけ隔離された2つの放射ソースからコヒーレントな縞パターン
を生成する工程と、 b)該コヒーレントな縞パターンを該オブジェクト表面上に投射する工程と、 c)該オブジェクト表面のイメージ内において該縞パターンおよび該スペック
ルを検出する工程と、 d)該オブジェクトを該2つの放射ソースに対して側方に該縞パターンの等位
相面に平行な経路に沿って平行移動させて、該縞パターンが該オブジェクト表面
に対して実質的に静止状態でありかつ該スペックルが変化するようにする工程と
、 e)該オブジェクト表面のイメージにおいて、該新規の縞パターンと、該変化
したスペックルとを検出する工程と、 f)該検出された縞パターンおよび該スペックルの変化に応答して、実質的に
スペックルの発生が無い状態で該縞パターンを判定する工程と、 を包含する、方法。 - 【請求項71】 工程d)および工程e)を繰り返し反復する工程をさらに
包含する、請求項70に記載の方法。 - 【請求項72】 前記2つの放射ソースはレーザ放射である、請求項70に
記載の方法。 - 【請求項73】 スペックルがオブジェクト表面上の測定ポイントに与える
影響を軽減する方法であって、 a)該オブジェクト表面上にコヒーレントな縞パターンを投射する工程と、 b)透過関数を有するレンズを提供する工程であって、該レンズの透過は該レ
ンズの端部において徐々に低下する、工程と、 c)該オブジェクト表面のイメージ内に該縞パターンを検出する工程であって
、該透過関数は、該スペックルが該測定に与える影響を実質的に低減する、工程
と、 を包含する、方法。 - 【請求項74】 スペックルがオブジェクト表面上の測定ポイントに与える
影響を軽減する方法であって、 a)第1の周波数を持ちかつある距離だけ隔離された2つの放射ソースから第1
の縞パターンを生成する工程と、 b)該第1の周波数を第2の周波数に変更することにより、該第1の縞パターン
を変更して第2の縞パターンを生成する工程と、 c)該第1の周波数と該第2の周波数との間の差に応答して該距離を変更する工
程であって、該第1の周波数と該第2の周波数との間の差に対する該距離の比は
実質的に一定である、工程と、 を包含する、方法。 - 【請求項75】 縞をオブジェクト表面上に投射する方法であって、 a)ある距離だけ隔離された2つの放射ソースを提供する工程であって、該ソ
ースはそれぞれ、スペクトル分布を有し、他方のソースに対してコヒーレントで
ある、工程と、 b)該オブジェクト表面上のポイントを、該ソースそれぞれからの放射によっ
て照射する工程と、 c)該ソースの1つを他方のソースに対して移動させる工程と、 d)該オブジェクト表面上のポイントによって散乱された放射を検出する工程
と、 を包含する、方法。 - 【請求項76】 前記ソースの移動と、前記オブジェクト表面上のポイント
によって散乱された放射の検出とに応答して位置情報を計算する工程をさらに包
含する、請求項75に記載の方法。 - 【請求項77】 前記ソースの1つからのスペクトル分布中のスペクトル成
分の位相を、該ソースのうちの他方のソースからのスペクトル分布中の各スペク
トル成分の位相を前記オブジェクト表面上のポイントにおいて測定した値に対し
て変化させる工程をさらに包含する、請求項75に記載の方法。 - 【請求項78】 前記ソースの1つを該ソースのうちの他方のソースに対し
て移動させる工程は、該ソースのうちの1つのスペクトル成分と該ソースのうち
の他方のソースの各スペクトル成分と間の距離を変化させて、該ソースの各スペ
クトル成分間の距離を変化させる工程を包含する、請求項75に記載の方法。 - 【請求項79】 前記距離は、前記ソースのスペクトル成分の各々から等距
離である中間点を含み、該ソースのうちの1つを該ソースのうちの他方のソース
に対して移動させる工程は、該ソースのうちの1つのスペクトル成分と、該ソー
スのうちの他方のソースの各スペクトル成分と間の距離を変化させて、該ソース
の各スペクトル成分間の距離を変化させ、該中間点を固定状態に保持する工程を
包含する、請求項78に記載の方法。 - 【請求項80】 前記オブジェクト表面上のポイントの三次元位置情報を判
定する工程をさらに包含する、請求項75に記載の方法。 - 【請求項81】 前記2つの放射ソースを提供する工程は、 a)あるスペクトル値を有する初期放射ビームを提供する工程と、 b)該初期ビーム放射から、第1の放射ビームを第1のビーム角度においてかつ
第2の放射ビームを第2のビーム角度において生成する工程と、 c)該第1の放射ビームおよび該第2の放射ビームをイメージングして、該2つ
の放射ソースを形成する工程と、 を包含する、請求項75に記載の方法。 - 【請求項82】 前記初期ビーム放射振幅変調する工程をさらに包含する、
請求項81に記載の方法。 - 【請求項83】 前記第1の放射ビームおよび第2の放射ビームを生成する
工程は、音響フィールドによる音響光学変調を前記初期ビームに行う工程を包含
する、請求項81に記載の方法。 - 【請求項84】 前記音響光学変調を行う工程は、前記音響フィールドに位
相変調を行う工程を包含する、請求項83に記載の方法。 - 【請求項85】 前記音響光学変調を行う工程は、前記音響フィールドを振
幅変調工程を包含する、請求項83に記載の方法。 - 【請求項86】 前記ソースのうちの1つを該ソースのうちの他方のソース
に対して移動させる工程は、音響フィールドに周波数変調を行うことによって前
記初期ビームに音響光学変調を行う工程を包含する、請求項75に記載の方法。 - 【請求項87】 a)ある距離だけ隔離された2つの放射ソースであって、
該ソースはそれぞれ、スペクトル分布を有し、該ソースのうちの他方のソースに
対してコヒーレントである、放射ソースと、 b)該ソースの各々を該ソースのうちの他方のソースに対して移動させる制御シ
ステムと、 c)オブジェクト表面上のポイントから散乱した放射を受信するような位置に配
置された検出器と、 を備える、縞をオブジェクト表面上に投射する装置。 - 【請求項88】 前記検出器からの信号を受信するプロセッサをさらに備え
、該プロセッサは、前記ソースの移動と、前記オブジェクト表面上のポイントか
ら散乱した放射とに応答して位置情報を計算する、請求項87に記載の装置。 - 【請求項89】 前記制御システムは、前記放射ソースの1つからのスペク
トル分布中のスペクトル成分の位相を、該放射ソースのうちの他方の放射ソース
からのスペクトル分布中の各スペクトル成分の位相を前記オブジェクト表面上の
ポイントにおいて測定した値に対して変化させる、請求項87に記載の装置。 - 【請求項90】 前記制御システムは、前記放射ソースのうちの1つの各ス
ペクトル成分を、該放射ソースのうちの他方のソースの各スペクトル成分に対し
て移動させて、該2つの放射ソースの各スペクトル成分間の距離を変化させる、
請求項87に記載の装置。 - 【請求項91】 前記距離は、前記ソースのスペクトル領域の各々から等距
離である中間点を含み、前記制御システムは、該ソースの各スペクトル成分を該
ソースのスペクトル成分のうちの他方のスペクトル成分に対して移動させて、該
距離を変化させ、該中間点を固定状態に保持する、請求項90に記載の装置。 - 【請求項92】 前記検出器と、前記オブジェクト表面上のポイントとの間
に配置されたイメージングシステムをさらに備え、該イメージングシステムは、
該検出器上のポイントをイメージングする、請求項87に記載の装置。 - 【請求項93】 前記検出器は複数の光検出器を備える、請求項87に記載
の装置。 - 【請求項94】 前記プロセッサは、複数の処理ユニットを有するマルチプ
ロセッサシステムを備え、前記複数の光検出器はそれぞれ、該複数の処理ユニッ
トの各々と電気的に通信する、請求項93に記載の装置。 - 【請求項95】 前記2つの放射ソースは、 a)あるスペクトル幅のビーム放射の初期ソースと、 b)該ビーム放射の初期ソースと光学的に連絡し、第1の光学ビームおよび第
2の光学ビームを生成するビーム分離器と、 c)該ビーム分離器と光学的に連絡するイメージングシステムであって、該2
つの放射ソースがそれぞれ該第1の光学ビームおよび該第2の光学ビームに対応
するように該2つの放射ソースを生成する、イメージングシステムと、 を備える、請求項87に記載の装置。 - 【請求項96】 前記ビーム分離器は回折格子である、請求項95に記載の
装置。 - 【請求項97】 1つのソースのスペクトル分布中の波長におけるポイント
は、他方ソースのスペクトル分布中の波長における各ポイントにおいて前記距離
だけ隔離され、該距離は、該ポイントの波長に実質的に比例する、請求項87に
記載の装置。 - 【請求項98】 前記制御システムは位相シフタをさらに備える、請求項8
7に記載の装置。 - 【請求項99】 前記位相シフタは、前記ビーム分離器の位置を前記ビーム
放射に対して変化させるトランスレータである、請求項98に記載の装置。 - 【請求項100】 前記制御システムは、焦点距離が可変であるズームレン
ズを備え、該ズームレンズは、該可変な焦点距離の変化に応答して該隔離距離を
変化させる、請求項87に記載の装置。 - 【請求項101】 前記制御システムは、焦点距離が異なる複数のレンズを
備え、該隔離距離は、該複数のレンズのうちの1つの各焦点距離に対応する、請
求項87に記載の装置。 - 【請求項102】 前記回折格子は空間光変調器である、請求項96に記載
の装置。 - 【請求項103】 前記空間光変調器は、調節可能な格子周期を有し、該格
子周期は、該空間光変調器に与えられる電気信号に応答する、請求項102に記
載の装置。 - 【請求項104】 前記ビーム分離器は音響光学変調器である、請求項95
に記載の装置。 - 【請求項105】 前記ビーム放射の初期ソースと光学的に連絡する振幅変
調器をさらに備える、請求項95に記載の装置。 - 【請求項106】 オーダ−ブロックマスクをさらに備える、請求項95に
記載の装置。 - 【請求項107】 前記オーダ−ブロックマスクは2重スリットのマスクで
ある、請求項106に記載の装置。 - 【請求項108】 前記ビーム分離器と光学的に連絡する振幅変調器をさら
に備え、該振幅変調器は、前景初期ソースの振幅を変調する、請求項95に記載
の装置。 - 【請求項109】 表面を有するオブジェクト上において、該オブジェクト
の表面上のポイントの三次元位置情報を判定する方法であって、 a)2つの放射ソースを提供する工程であって、該2つの放射ソースはそれぞ
れ、空間スペクトル分布を有し、該2つのソースのうちの他方のソースに対して
コヒーレントである、工程と、 b)該表面上のポイントに検出器を設ける工程と、 c)該オブジェクト表面上のポイントを該ソースの各々からの放射で照射する
工程と、 d)該ソースの各々を互いに移動させる工程と、 e)該オブジェクト表面上のポイントにおいて該放射を検出する工程と、 f)該ソースの移動と、該オブジェクト表面上のポイントにおいて検出された
放射とに応答して位置情報を計算する工程と、 を包含する、方法。 - 【請求項110】 前記2つの放射ソースは、狭帯域である空間スペクトル
分布を有する、請求項109に記載の方法。 - 【請求項111】 前記ソースの1つからの放射の空間スペクトル分布中の
スペクトル成分の位相を、該ソースのうちの他方のソースからの放射の空間スペ
クトル分布中の各スペクトル成分の位相を前記オブジェクト表面上のポイントに
おいて測定した値に対して変化させる工程をさらに包含する、請求項109に記
載の方法。 - 【請求項112】 前記放射は仮想放射ソースである、請求項109に記載
の方法。 - 【請求項113】 前記ソースの各々を移動させる工程は、該ソースのうち
の1つのスペクトル成分と該ソースのうちの他方のソースの各スペクトル成分と
間の距離を変化させて、該ソースの各スペクトル成分間の距離を変化させる工程
を包含する、請求項109に記載の方法。 - 【請求項114】 前記距離は、前記ソースのスペクトル成分の各々から等
距離である中間点を含み、該ソースの各々を移動させる工程は、該ソースのうち
の1つのスペクトル成分と、該ソースのうちの他方のソースの各スペクトル成分
と間の距離を変化させて、該ソースの各スペクトル成分間の距離を変化させ、該
中間点を固定状態に保持する工程を包含する、請求項113に記載の方法。 - 【請求項115】 表面を有するオブジェクト上において、該オブジェクト
の表面上のポイントの三次元位置情報を判定する装置であって、 a)空間スペクトル分布を有し、互いにコヒーレントな2つの放射ソースと、 b)該ソースの各々を相互に移動させる制御システムと、 c)該オブジェクト表面上のポイントを照射する放射を受信するように、該オ
ブジェクト表面上の該ポイントに配置された検出器と、 d)該検出器からの信号を受信するプロセッサであって、該ソースの移動と、
該オブジェクト表面上の該ポイントにおいて受信された放射とに応答して該オブ
ジェクト表面上の該ポイントの位置情報を計算する、プロセッサと、 を備える、装置。 - 【請求項116】 前記2つの放射ソースは、狭帯域である空間スペクトル
分布を有する、請求項115に記載の方法。 - 【請求項117】 前記制御システムは、前記ソースの1つからの放射の空
間スペクトル分布中のスペクトル成分の位相を、該ソースのうちの他方のソース
からの放射の空間スペクトル分布中の各スペクトル成分の位相を前記オブジェク
ト表面上のポイントにおいて測定した値に対して変化させる、請求項115に記
載の方法。 - 【請求項118】 前記放射ソースは仮想放射ソースである、請求項115
に記載の装置。 - 【請求項119】 前記制御システムは、前記放射ソースのうちの1つの各
スペクトル成分を、該放射ソースのうちの他方のソースの各スペクトル成分に対
して移動させて、該2つの放射ソースの各スペクトル成分間の距離を変化させる
、請求項115に記載の方法。 - 【請求項120】 前記距離は、前記ソースのスペクトル領域の各々から等
距離である中間点を含み、前記制御システムは、該ソースの各スペクトル成分を
該ソースのスペクトル成分のうちの他方のスペクトル成分に対して移動させて、
該距離を変化させ、該中間点を固定状態に保持する、請求項119に記載の装置
。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/480,043 | 2000-01-10 | ||
US09/480,043 US6690474B1 (en) | 1996-02-12 | 2000-01-10 | Apparatus and methods for surface contour measurement |
PCT/US2001/000161 WO2001051886A1 (en) | 2000-01-10 | 2001-01-03 | Apparatus and methods for surface contour measurement |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010267905A Division JP5241806B2 (ja) | 2000-01-10 | 2010-11-30 | 表面輪郭測定のための装置および方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2003519786A true JP2003519786A (ja) | 2003-06-24 |
Family
ID=23906443
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2001552052A Pending JP2003519786A (ja) | 2000-01-10 | 2001-01-03 | 表面輪郭測定のための装置および方法 |
JP2010267905A Expired - Fee Related JP5241806B2 (ja) | 2000-01-10 | 2010-11-30 | 表面輪郭測定のための装置および方法 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010267905A Expired - Fee Related JP5241806B2 (ja) | 2000-01-10 | 2010-11-30 | 表面輪郭測定のための装置および方法 |
Country Status (7)
Country | Link |
---|---|
EP (1) | EP1247070B1 (ja) |
JP (2) | JP2003519786A (ja) |
AT (1) | ATE376165T1 (ja) |
AU (1) | AU2001227574A1 (ja) |
CA (1) | CA2397095C (ja) |
DE (1) | DE60130968T2 (ja) |
WO (1) | WO2001051886A1 (ja) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007113958A (ja) * | 2005-10-18 | 2007-05-10 | Yamatake Corp | 3次元計測装置、3次元計測方法、及び3次元計測プログラム |
JP2008257238A (ja) * | 2007-03-31 | 2008-10-23 | Sony Deutsche Gmbh | 画像生成マイクロディスプレイを均一に照射するための照明装置および方法 |
JP2008281543A (ja) * | 2007-05-08 | 2008-11-20 | Koh Young Technology Inc | 多方向映写式モアレ干渉計及びこれを用いた検査方法 |
JP2016537612A (ja) * | 2013-09-25 | 2016-12-01 | アールト コルケアコウルスエーティ | モデリング構成及び三次元表面のトポグラフィーをモデリングする方法及びシステム |
WO2017183181A1 (ja) * | 2016-04-22 | 2017-10-26 | オリンパス株式会社 | 三次元形状測定装置 |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7496004B2 (en) | 2003-05-02 | 2009-02-24 | Sony Corporation | Data reproducing apparatus, data reproducing method, data recording and reproducing apparatus, and data recording and reproducing method |
EP1548400A1 (de) * | 2003-12-23 | 2005-06-29 | Vogel und Ploetscher GmbH | Messgerät zum Inspizieren von Schienen |
DE102007002880B4 (de) * | 2007-01-15 | 2011-08-25 | DMG Microset GmbH, 33689 | Verfahren und Vorrichtung zum optischen Ausmessen eines Objekts, insbesondere eines Werkstücks oder Werkzeugs |
US8531650B2 (en) | 2008-07-08 | 2013-09-10 | Chiaro Technologies LLC | Multiple channel locating |
EP2166305B1 (de) * | 2008-09-23 | 2012-05-16 | Sick Ag | Beleuchtungseinheit und Verfahren zur Projektion eines Beleuchtungsmusters |
AU2011265572A1 (en) * | 2011-12-23 | 2013-07-11 | Canon Kabushiki Kaisha | Structured light system for robust geometry acquisition |
DE102013201469B4 (de) | 2013-01-30 | 2023-01-12 | Rohde & Schwarz GmbH & Co. Kommanditgesellschaft | Verfahren und Vorrichtung zur Oberflächenbestimmung |
LU92173B1 (en) * | 2013-03-20 | 2014-09-22 | Iee Sarl | Distance determination method |
US9562760B2 (en) | 2014-03-10 | 2017-02-07 | Cognex Corporation | Spatially self-similar patterned illumination for depth imaging |
US10571668B2 (en) | 2015-05-09 | 2020-02-25 | Cognex Corporation | Catadioptric projector systems, devices, and methods |
RU172076U9 (ru) * | 2016-06-14 | 2017-08-22 | Александр Евгеньевич Соколов | Инфракрасное лазерное сканирующее устройство |
US11300402B2 (en) * | 2017-05-31 | 2022-04-12 | Hewlett-Packard Development Company, L.P. | Deriving topology information of a scene |
EP3444782B1 (en) | 2017-08-19 | 2022-03-30 | Cognex Corporation | Coding distance topologies for structured light patterns for 3d reconstruction |
CN109635619B (zh) | 2017-08-19 | 2021-08-31 | 康耐视公司 | 用于三维重建的结构化光图案的编码距离拓扑 |
FR3081592B1 (fr) * | 2018-05-25 | 2021-05-14 | Vit | Systeme de determination d'images tridimensionnelles |
US11143503B2 (en) | 2018-08-07 | 2021-10-12 | Kimball Electronics Indiana, Inc. | Interferometric waviness detection systems |
WO2020076336A1 (en) * | 2018-10-12 | 2020-04-16 | Electric Power Research Institute, Inc. | Method for measuring surface characteristics in optically distorting media |
CN110779461B (zh) * | 2019-10-16 | 2021-09-21 | 中国航空工业集团公司洛阳电光设备研究所 | 一种折反汇聚光路中平面反射镜面型的测试设备及方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999063304A1 (en) * | 1998-06-04 | 1999-12-09 | Massachusetts Institute Of Technology | Apparatus and methods for surface contour measurement |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2641945B2 (ja) * | 1989-08-31 | 1997-08-20 | 日本電気株式会社 | 距離画像取得方法及び装置 |
JPH06241727A (ja) * | 1993-02-19 | 1994-09-02 | Nikon Corp | 位置検出装置 |
US5870191A (en) * | 1996-02-12 | 1999-02-09 | Massachusetts Institute Of Technology | Apparatus and methods for surface contour measurement |
-
2001
- 2001-01-03 AU AU2001227574A patent/AU2001227574A1/en not_active Abandoned
- 2001-01-03 DE DE60130968T patent/DE60130968T2/de not_active Expired - Lifetime
- 2001-01-03 CA CA2397095A patent/CA2397095C/en not_active Expired - Lifetime
- 2001-01-03 WO PCT/US2001/000161 patent/WO2001051886A1/en active IP Right Grant
- 2001-01-03 AT AT01901698T patent/ATE376165T1/de not_active IP Right Cessation
- 2001-01-03 JP JP2001552052A patent/JP2003519786A/ja active Pending
- 2001-01-03 EP EP01901698A patent/EP1247070B1/en not_active Expired - Lifetime
-
2010
- 2010-11-30 JP JP2010267905A patent/JP5241806B2/ja not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999063304A1 (en) * | 1998-06-04 | 1999-12-09 | Massachusetts Institute Of Technology | Apparatus and methods for surface contour measurement |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007113958A (ja) * | 2005-10-18 | 2007-05-10 | Yamatake Corp | 3次元計測装置、3次元計測方法、及び3次元計測プログラム |
JP4688625B2 (ja) * | 2005-10-18 | 2011-05-25 | 株式会社山武 | 3次元計測装置、3次元計測方法、及び3次元計測プログラム |
JP2008257238A (ja) * | 2007-03-31 | 2008-10-23 | Sony Deutsche Gmbh | 画像生成マイクロディスプレイを均一に照射するための照明装置および方法 |
JP2008281543A (ja) * | 2007-05-08 | 2008-11-20 | Koh Young Technology Inc | 多方向映写式モアレ干渉計及びこれを用いた検査方法 |
US7894077B2 (en) | 2007-05-08 | 2011-02-22 | Koh Young Technology Inc. | Multi-directional projection type moire interferometer and inspection method of using the same |
JP2016537612A (ja) * | 2013-09-25 | 2016-12-01 | アールト コルケアコウルスエーティ | モデリング構成及び三次元表面のトポグラフィーをモデリングする方法及びシステム |
WO2017183181A1 (ja) * | 2016-04-22 | 2017-10-26 | オリンパス株式会社 | 三次元形状測定装置 |
JPWO2017183181A1 (ja) * | 2016-04-22 | 2019-02-28 | オリンパス株式会社 | 三次元形状測定装置 |
US10753731B2 (en) | 2016-04-22 | 2020-08-25 | Olympus Corporation | Three-dimensional form measurement device |
Also Published As
Publication number | Publication date |
---|---|
EP1247070B1 (en) | 2007-10-17 |
WO2001051886A1 (en) | 2001-07-19 |
EP1247070A1 (en) | 2002-10-09 |
DE60130968T2 (de) | 2008-07-10 |
CA2397095C (en) | 2010-04-27 |
AU2001227574A1 (en) | 2001-07-24 |
JP5241806B2 (ja) | 2013-07-17 |
WO2001051886A9 (en) | 2002-10-31 |
JP2011053224A (ja) | 2011-03-17 |
ATE376165T1 (de) | 2007-11-15 |
CA2397095A1 (en) | 2001-07-19 |
DE60130968D1 (de) | 2007-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5241806B2 (ja) | 表面輪郭測定のための装置および方法 | |
US6690474B1 (en) | Apparatus and methods for surface contour measurement | |
US6031612A (en) | Apparatus and methods for contour measurement using movable sources | |
US10794689B2 (en) | Autofocus system and method | |
US7777895B2 (en) | Linear-carrier phase-mask interferometer | |
KR101590241B1 (ko) | 광학특성 측정장치 및 광학특성 측정방법 | |
EP1724550B1 (en) | Interferometer and shape measuring method | |
US9494411B2 (en) | Three-dimensional shape measuring device, method for acquiring hologram image, and method for measuring three-dimensional shape | |
US8736847B2 (en) | Method and apparatus for imaging | |
EP1717546B1 (en) | Interferometer and method of calibrating the interferometer | |
US20050046865A1 (en) | Pixelated phase-mask interferometer | |
US20110298896A1 (en) | Speckle noise reduction for a coherent illumination imaging system | |
KR102426103B1 (ko) | 개선된 홀로그래픽 복원 장치 및 방법 | |
US10444004B2 (en) | Phase shift interferometer | |
US7561279B2 (en) | Scanning simultaneous phase-shifting interferometer | |
JP7082137B2 (ja) | スペクトル制御干渉法による曲率半径測定 | |
TW202020400A (zh) | 表面形狀量測裝置以及表面形狀量測方法 | |
KR101554205B1 (ko) | 파장주사 층밀리기 간섭계를 이용한 자유곡면의 삼차원 형상 측정 장치 | |
JP3934131B2 (ja) | 同軸型空間光干渉断層画像計測装置 | |
US20110299090A1 (en) | Real-time interferometer | |
JP2003294418A (ja) | 微小周期構造評価装置及び微小周期構造評価方法 | |
JPH05296879A (ja) | 光学性能測定方法及び装置 | |
Jackson et al. | Demonstration of Interference Pattern Structured Illumination Imaging | |
KR20220107144A (ko) | 개선된 홀로그래픽 복원 장치 및 방법 | |
JP2004116997A (ja) | 被検物面の測定方法及び測定装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20071220 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20100827 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100903 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20101130 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110307 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20110805 |