[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

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
Application number
JP2001552052A
Other languages
English (en)
Inventor
ライル シャーリー,
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Massachusetts Institute of Technology
Original Assignee
Massachusetts Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from US09/480,043 external-priority patent/US6690474B1/en
Application filed by Massachusetts Institute of Technology filed Critical Massachusetts Institute of Technology
Publication of JP2003519786A publication Critical patent/JP2003519786A/ja
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/25Measuring 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/2509Color coding
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/2441Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures using interferometry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/25Measuring 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/2518Projection by scanning of the object
    • G01B11/2527Projection by scanning of the object with phase change by in-plane movement of the patern
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/25Measuring 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/2531Measuring 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

(57)【要約】 オブジェクト表面上のポイントの三次元位置情報を測定する装置および方法。一実施形態において、上記方法は、スペクトル分布を有する2つの放射ソースを提供する工程と、上記ソースの各々で上記表面を照射して第1の縞パターンを生成する工程と、上記第1の縞パターンを第2の位置に移動させる工程と、第1のラップされたサイクルマップを生成する工程と、該第1の縞パターン中の縞数を推定する工程と、該第1の縞パターンを変化させる工程と、該第2の縞パターンを第2の位置に移動させる工程と、第2のラップされたサイクルマップを生成する工程と、該第2の縞パターン中の縞数を推定する工程と、上記第2の縞パターンおよび上記第2のラップされたサイクルマップ中の推定された縞数に応答して位置情報を判定する工程とを含む。

Description

【発明の詳細な説明】
【0001】 (関連出願との相互参照) 本願は、1996年2月12日に出願された米国特許第5,870,191号
の一部継続出願である、1999年2月2日に出願された米国特許出願第09/
241,354号の一部継続出願であり、そして1998年6月4日に出願され
た米国仮特許出願第60/087,960号の優先権を主張する。
【0002】 (政府援助) 本明細書に記載する研究は、米国空軍によって授与された、連邦契約第F19
628−95−L−002によって援助された。政府は本発明に特定の権利を有
し得る。
【0003】 (発明の分野) 本発明は、表面測定の分野に関し、具体的には、非接触表面測定の分野に関す
る。
【0004】 (発明の背景) 寸法計測学、すなわち、オブジェクトのサイズおよび形状の測定は、マシーン
が、ほとんどの製造および多くのサブアセンブリから構成される複雑なオブジェ
クトのアセンブリを実行する、今日の製造環境において非常に重要である。構成
要素が適切に嵌合し合うことを保証するには、自動車などの複雑なアセンブリの
各構成要素の形状およびサイズの公差を小さく保持する必要がある。
【0005】 理想的には、形状およびサイズのこのような測定は、物理的接触無しに達成さ
れて、測定を行う際の時間を節減する。多くの非接触測定方法は、利用可能なマ
シーンビジョンシステムを利用する。表面輪郭情報の測定は、深度情報が損失し
たり、解釈が難しい場合が多いため、マシーンビジョンシステムの特に難しい問
題である。深度情報の損失および利用可能な情報を解釈する際の難しさを補償す
るために、多くのマシーンビジョンシステムは、光を用いてオブジェクトの表面
上にモアレ縞を生成し、輪郭情報を得る。モアレ技術の1つの欠点は、異なるサ
イズのオブジェクトに対するその柔軟性の無さである。異なるサイズのオブジェ
クトは、新しい対応する物理的設備を必要とし得る。この欠点に起因して、大規
模なオブジェクトのモアレ技術を用いることが難しい。モアレ技術の別の欠点は
、この技術を用いて得られる解像度が多くの用途にとって十分に高くない場合が
あることである。
【0006】 干渉方法も、表面の詳細な測定が必要な場合に用いられてきた。干渉システム
が表面輪郭情報を提供するが、干渉システムは測定されるオブジェクトおよび用
いられる照射源の両方の変化の影響を受けやすい。
【0007】 必要なのは、縞模様の道理によるあいまい性を解決し、スペックル効果に起因
する劣化を緩和し、そして高い解像度を達成することによって、これらの問題を
回避する技術である。
【0008】 (発明の要旨) 本発明は、表面を有するオブジェクト上において、上記オブジェクトの表面上
のポイントの三次元位置情報を判定する方法に関する。上記方法は、スペクトル
領域の空間分布を有する2つの放射ソースを提供する工程;上記ソースそれぞれ
からの放射で上記表面を照射して、上記表面上の第1の位置に第1の縞パターン
を生成する工程;上記第1の縞パターンを第2の位置に移動させる工程;上記第
1の縞パターンの上記第1の位置および上記第2の位置に応じて、第1のラップ
されたサイクルマップを生成する工程;上記第1の縞パターン中の縞数を推定す
る工程;上記第1の縞パターンを変更して、第1の位置において第2の縞パター
ンを生成する工程;上記第2の縞パターンを第2の位置に移動させる工程;上記
第2の縞パターンの上記第1の位置および上記第2の位置に応じて、第2のラッ
プされたサイクルマップを生成する工程;上記第1の縞パターン中の推定された
縞数に応じて、上記第2の縞パターン中の縞数を推定する工程;および上記第2
の縞パターンおよび上記第2のラップされたサイクルマップ中の推定された縞数
に応じて、上記ポイントの上記表面上での三次元位置情報を判定する工程を含む
【0009】 本発明は、表面上の所定の間隔の縞パターンに対応するラップされたサイクル
マップを合成する方法にさらに関する。上記方法は、スペクトル領域の空間分布
を有する2つの放射ソースを提供する工程;上記ソースそれぞれからの放射で上
記表面を照射して、上記表面上の第1の位置に第1の縞パターンを生成する工程
;上記第1の縞パターンを上記表面上の第2の位置に移動させる工程;上記第1
の縞パターンの上記第1の位置および上記第2の位置に応答して、第1のラップ
されたサイクルマップを生成する工程;上記第1の縞パターンを変更して、第1
の位置において第2の縞パターンを生成する工程;上記第2の縞パターンを第2
の位置に移動させる工程;上記第2の縞パターンの上記第1の位置および上記第
2の位置に応答して、第2のラップされたサイクルマップを生成する工程;上記
第1のラップされたサイクルマップから上記第2のラップされたサイクルマップ
を減算する工程;および上記第2のラップされたサイクルマップと上記第1のラ
ップされたサイクルマップと間の差をラップして、上記所定の間隔の縞パターン
に対応するラップされたサイクルマップを生成する工程を含む。
【0010】 本発明は、表面を有するオブジェクト上において、上記オブジェクトの表面上
のポイントの三次元位置情報を判定する方法にさらに関する。上記方法は、第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のラップされたサ
イクルマップ中の推定された強度サイクルに応じて、三次元位置情報を計算する
工程を含む。
【0011】 本発明は、表面を有するオブジェクト上において、上記オブジェクトの表面上
のポイントの三次元位置情報を判定する方法にさらに関する。上記方法は、互い
にコヒーレントな第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のラップされたサ
イクルマップとの間の差をラップして、ラップされたサイクルマップを生成する
工程を含む。
【0012】 本発明は、オブジェクトの表面上に縞を投射する装置にさらに関する。上記装
置は、スペクトル分布を有する2つの放射ソース;上記2つのソースと光学的に
連絡したコリメータであって、2本の実質的に平行な広域幅放射ビームを生成す
る、コリメータ;上記コリメータと光学的に連絡した回折格子;および上記回折
格子と光学的に連絡したレンズであって、スペクトル領域の空間分布を有する2
つの放射イメージを生成する、レンズを有する。
【0013】 本発明は、オブジェクトの表面上に縞を投射する方法にさらに関する。上記方
法は、ある距離の分だけ隔離された2つの放射ソースを提供する工程であって、
上記放射ソースはスペクトル分布を有する放射を生成する、工程;上記放射を平
行にして、2本の実質的に平行な放射ビームを生成する工程;上記平行な放射ビ
ームのスペクトル成分を描写する工程;およびスペクトル成分の空間分布を有す
る2つの放射イメージを生成する工程を含む。
【0014】 本発明は、表面を有するオブジェクト上において、上記オブジェクトの表面上
のポイントの三次元位置情報を判定する方法にさらに関する。上記方法は、第1
の周期的パターンを投射して、上記表面上の第1の位置に第1の投射パターンを
生成する工程;上記第1の投射パターンを第2の位置に移動させる工程;上記第
1の投射パターンの上記第1の位置および上記第2の位置に応答して、第1のラ
ップされたサイクルマップを生成する工程;上記第1の投射パターン中のサイク
ル数を推定する工程;上記第1の投射パターンを変更して、第1の位置において
第2の投射パターンを生成する工程;上記第2の投射パターンを第2の位置に移
動させる工程;上記第2の投射パターンの上記第1の位置および上記第2の位置
に応じて、第2のラップされたサイクルマップを生成する工程;上記第1の投射
パターン中の推定されたサイクル数に応じて、上記第2の投射パターン中のサイ
クル数を推定する工程;および上記第2の投射パターン中の推定されたサイクル
数および上記第2のラップされたサイクルマップに応じて、上記表面を判定する
工程を含む。
【0015】 本発明は、表面上の所定の間隔の投射された周期的パターンに対応するラップ
されたサイクルマップを合成する方法にさらに関する。上記方法は、第1の周期
的パターンを投射して、上記表面上の第1の位置に第1の投射パターンを生成す
る工程;上記第1の投射パターンを上記表面上の第2の位置に移動させる工程;
上記第1の投射パターンの上記第1の位置および上記第2の位置に応じて、第1
のラップされたサイクルマップを生成する工程;上記第1の投射パターンを変更
して、第1の位置において第2の投射パターンを生成する工程;上記第2の投射
パターンを第2の位置に移動させる工程;上記第2の投射パターンの上記第1の
位置および上記第2の位置に応じて、第2のラップされたサイクルマップを生成
する工程;上記第1のラップされたサイクルマップから上記第2のラップされた
サイクルマップを減算する工程;および上記第2のラップされたサイクルマップ
と上記第1のラップされたサイクルマップと間の差をラップして、上記所定の間
隔の投射された周期的パターンに対応するラップされたサイクルマップを生成す
る工程を含む。
【0016】 本発明は、スペックルがオブジェクトの表面上の測定ポイントに与える影響を
軽減する方法にさらに関する。上記方法は、コヒーレントな縞パターンを生成す
る工程;上記縞パターンが上記オブジェクトの表面を実質的にかするように、上
記コヒーレントな縞パターンを上記オブジェクトの表面上に光路に沿って投射す
る工程;および上記オブジェクトの表面のイメージ中に上記縞パターンおよび上
記スペックルを検出する工程であって、上記オブジェクトの表面に対する法線が
上記光路に実質的に直交する、工程を含む。
【0017】 本発明は、スペックルがオブジェクトの表面上の測定ポイントに与える影響を
軽減する方法にさらに関する。上記方法は、ある距離の分だけ隔離された2つの
放射ソースからコヒーレントな縞パターンを生成する工程;上記コヒーレントな
縞パターンを上記オブジェクトの表面上に投射する工程;上記オブジェクトの表
面のイメージ中に上記縞パターンおよび上記スペックルを検出する工程;上記縞
パターンが実質的に静止したままで、そして上記スペックルが変化するように上
記2つのソースを平行移動させる工程;上記オブジェクトの表面のイメージ中に
上記新しい縞パターンおよび上記変化したスペックルを検出する工程;ならびに
上記検出された縞パターンおよび上記スペックルの変化に応じて、実質的にスペ
ックルの無い上記縞パターンを判定する工程を含む。
【0018】 本発明は、スペックルがオブジェクトの表面上の測定ポイントに与える影響を
軽減する方法にさらに関する。上記方法は、ある距離の分だけ隔離された2つの
放射ソースからコヒーレントな縞パターンを生成する工程;上記コヒーレントな
縞パターンを上記オブジェクトの表面上に投射する工程;上記オブジェクトの表
面のイメージ中に上記縞パターンおよびスペックルを検出する工程;上記縞パタ
ーンが上記オブジェクトの表面に対して実質的に静止したままで、そして上記ス
ペックルが変化するように、上記オブジェクトを上記縞パターンの等位相平面に
平行な経路に沿って上記2つの放射ソースに対して平行移動させる工程;上記オ
ブジェクトの表面のイメージ中に上記新しい縞パターンおよび上記変化したスペ
ックルを検出する工程;ならびに上記検出された縞パターンおよび上記スペック
ルの変化に応じて、上記実質的にスペックルの無い縞パターンを判定する工程を
含む。
【0019】 本発明は、スペックルがオブジェクトの表面上の測定ポイントに与える影響を
軽減する方法にさらに関する。上記方法は、コヒーレントな縞パターンを上記オ
ブジェクトの表面上に投射する工程;透過作用を有するレンズを提供する工程で
あって、上記レンズの透過は上記レンズ端部において徐々に低下する、工程;お
よびオブジェクトの表面のイメージ中に上記縞パターンを検出する工程であって
、上記透過作用は上記スペックルが上記測定に与える影響を実質的に減少させる
、工程を含む。
【0020】 本発明は、スペックルがオブジェクトの表面上の測定ポイントに与える影響を
軽減する方法にさらに関する。上記方法は、第1の周波数を有し、そしてある距
離の分だけ隔離された2つの放射ソースから第1の縞パターンを生成する工程;
上記第1の周波数を第2の周波数に変化させることによって、上記第1の縞パタ
ーンを変化させて、第2の縞パターンを生成する工程;および上記第1の周波数
と上記第2の周波数との差に応答して上記距離を変化させる工程であって、上記
距離と、上記第1の周波数と上記第2の周波数との差との比は実質的に一定であ
る、工程を含む。
【0021】 本発明は、オブジェクトの表面上に縞を投射する方法にさらに関する。上記方
法は、ある距離の分だけ隔離された2つの放射ソースを提供する工程であって、
上記放射ソースは、スペクトル分布を有し、そしてもう一方のソースに対してコ
ヒーレントである、工程;上記ソースそれぞれからの放射で上記オブジェクトの
表面上のポイントを照射する工程;上記ソースのうちの一つをもう一方のソース
に対して移動させる工程;および上記オブジェクトの表面上のポイントによって
散乱された放射を検出する工程を含む。本発明の別の実施形態において、2つの
放射ソースを提供する工程は、スペクトル幅を有する初期放射ビームを提供する
工程;上記初期放射ビームから、第1の放射ビームを第1のビーム角度で、第2
の放射ビームを第2のビーム角度で生成する工程;ならびに上記第1の放射ビー
ムおよび上記第2の放射ビームをイメージングして上記2つの放射ソースを形成
する工程を含む。
【0022】 本発明は、オブジェクトの表面上に縞を投射する装置にさらに関する。上記装
置は、ある距離の分だけ隔離された2つの放射ソースであって、各ソースは、ス
ペクトル分布を有し、そしてもう一方のソースに対してコヒーレントである、放
射ソース;上記ソースそれぞれをもう一方のソースに対して移動させる制御シス
テム;および上記オブジェクトの表面上のポイントから散乱された放射を受け取
るように配置された検出器を含む。本発明の別の実施形態において、2つの放射
ソースは、スペクトル幅を有する放射ビームの初期ソース;上記放射ビームの初
期ソースと光学的に連絡したビーム分離器であって、第1の光学ビームおよび第
2の光学ビームを生成する、ビーム分離器;および上記ビーム分離器に光学的に
連絡したイメージングシステムであって、上記イメージングシステムは、上記2
つの放射ソースを生成し、そしてそれぞれの放射ソースは上記第1の光学ビーム
および上記第2の光学ビームに対応する、イメージングシステムを含む。
【0023】 本発明は、表面を有するオブジェクト上において、上記オブジェクトの表面上
のポイントの三次元位置情報を判定する方法にさらに関する。上記方法は、2つ
の放射ソースを提供する工程であって、上記放射ソースは、スペクトル分布を有
し、そして上記2つの放射ソースのもう一方に対してコヒーレントである、工程
;検出器を上記表面上のポイントに提供する工程;上記オブジェクトの上記表面
のポイントを上記ソースのそれぞれからの上記放射で照射する工程;上記ソース
のそれぞれをもう一方に対して移動させる工程;上記オブジェクトの表面上のポ
イントに上記放射を検出する工程;ならびに上記ソースの移動および上記オブジ
ェクトの表面上のポイントに検出された放射に対して位置情報を計算する工程を
含む。
【0024】 本発明は、表面を有するオブジェクト上において、上記オブジェクトの表面上
のポイントの三次元位置情報を判定する装置にさらに関する。上記装置は、2つ
の放射ソースであって、スペクトル分布を有し、そして上記2つの放射ソースの
もう一方に対してコヒーレントである、放射ソース;上記ソースのそれぞれをも
う一方に対して移動させる制御システム;上記オブジェクトの表面上のポイント
に配置された検出器であって、上記オブジェクトの表面上のポイントを照射する
放射を受け取る、検出器;上記検出器から信号を受け取るプロセッサ;ならびに
上記ソースの移動および上記オブジェクトの表面上のポイントで受け取られた放
射に対して上記オブジェクトの表面上のポイントの位置情報を計算するプロセッ
サを含む。
【0025】 (好適な実施形態の詳細) 本発明は上掲の特許請求の範囲に特に示される。本発明の上述およびさらなる
利点は、添付の図面に関する以下の説明を参照することによってよりよく理解さ
れ得る。
【0026】 本発明の実施形態を説明しながら、「ソース」および「放射ソース」を示す。
これらの用語は、非常に局所化された放射ソースを含む、任意の放射ソースを指
すよう意図される。
【0027】 図1および簡単な概略を参照して、二つの放射ソース、PおよびPは、固
定距離Dの分だけ隔離されており、そしてそれぞれ(x、y、z)および
(x、y、z)の空間座標をそれぞれ有する。それぞれのソース、P
よびPからの放射は、もう一方のソースからの放射に対してコヒーレントであ
る。各ソース、PおよびPは、各発散放射ビーム12および14をオブジェ
クト10の表面上のポイントPに方向付ける。それぞれの放射ソースPおよ
びPから表面上のポイントPまでの距離は、RおよびRによってそれぞ
れ示される。Ψは起点からポイントPまで延びる線と、ソース、PとP
の間を延びる線との間の角度であり、θはz軸とソース、PとPとの間を
延びる線との間の角度であり、αはPから見た場合の2つのソースポイントか
ら延びた辺で囲まれた角度の半分の角度である。各ビーム12、14はもう一方
のビーム14、12と同じ方向に実質的に偏光され、そして、オブジェクト10
上の異なる領域を同時に照射するように独立して走査可能であり得る。あるいは
、オブジェクト10全体が同時に照射され得る。
【0028】 ポイントPによって散乱された光20は光検出器22によって検出される。
一実施形態において、光検出器22は1アレイの光検出器の素子を含んで、測定
されるオブジェクト10の二次元画像を提供する。さらなる実施形態において、
1アレイの光検出器の素子は電荷結合素子(CCD)である。検出器22は、1
つ以上の個々の信号を含み、各信号は検出器22の光検出器の素子のうちの対応
する1つに関連付けられた出力信号26を提供する。
【0029】 好適な実施形態において、焦点調節要素24はオブジェクト10の表面上のポ
イントPと光検出器22との間に配置されて、検出器22上にポイントP
含むオブジェクトの照射された部分をイメージングする。オブジェクトの表面が
粗いため、そして照射放射がコヒーレントであるため、焦点調節されたイメージ
はスペックルを有する。光検出器22からの出力信号26は処理装置28への入
力信号である。
【0030】 一実施形態において、偏光子30は焦点調節要素24と検出器22との間に配
置される。偏光子30の配置方向は、散乱光20の主な偏光成分との合致が最大
化される方向に向けられる。これにより、スペックルのコントラストまたは縞パ
ターンのコントラストが向上する。この構成により、オブジェクト10の表面か
ら散乱された光に関連付けられた信号対雑音比が最大化される。
【0031】 一実施形態において、プロセッサ28は一つのプロセッサであり、検出器アレ
イ22の光検出器の素子のそれぞれに関連付けられた検出器の出力信号26を受
けて動作する。別の実施形態において、プロセッサ28は、複数の個々のプロセ
ッサを有するマルチプロセッサであり、そして各光検出器の素子は入力信号をプ
ロセッサのそれぞれ一つに提供する。検出器22がCCDアレイである、また別
の実施形態において、複数のCCD素子は入力信号をマルチプロセッサの各プロ
セッサに提供する。マルチプロセッサの構成を用いると、複数の個々の光素子か
らの信号の計算が実質的に同時に行われ、これにより、信号処理速度が向上する
【0032】 制御装置32は、放射ソースPおよびPのうちの一つの放射の位相を、オ
ブジェクト10の表面上のポイントPで測定されるもう一方のソースからの放
射の位相に対して変更するように、放射ソースPおよびPの動作を制御する
。プロセッサ28は、信号線またはバス34を介して制御装置32と通信し得る
。例えば、特定の用途において、プロセッサ28が検出器22からの信号を処理
する場合には、オブジェクト10の表面上を介するソースPおよびPの走査
に対応して、またはソースからの放射が掃引される周期に対応して、特定の回数
行うことが望ましい場合がある。このような走査および周波数掃引動作が制御装
置32によって制御されるため、制御装置32とプロセッサ28との間の通信が
これらの環境において望ましい。制御装置32およびプロセッサ28は物理的に
別々の装置であっても、単一の処理システムによって実現されてもよいことが理
解される。
【0033】 ここで図2を参照して、一実施形態において、放射ソースPおよびPは、
チューナブルレーザ40から放出される放射から形成される。チューナブルレー
ザ40によって放出された放射ビーム44はビームスプリッタ48によって分割
される。放射ビーム50はビームスプリッタ48によって反射された後、レンズ
52によって発散するようになる。次いで、発散ビームは、移動可能な照準ミラ
ー54によって反射される。照準ミラー54によって反射される放射ビームは、
コヒーレントな放射ソースのうちの一つ、Pを提供する。同様に、ビームスプ
リッタ48を透過した放射ビーム46は、レンズ58によって発散するようにな
り、発散ビームを第2の移動可能な照準ミラー60に方向付ける。ミラー60に
よって反射された放射ビームは、第2の放射ソースPを提供する。照準ミラー
54および62は、回転可能であり、オブジェクト10の表面を選択的に照射し
得る。照準ミラー54および62はさらに、移動可能であり、ソースPおよび
の位置を変更し得る。
【0034】 図2aを参照して、放射ソースPおよびPの別の実施形態が、放射ビーム
44を提供する回転可能なレーザソース40を含むように図示される。放射ビー
ム44はレンズ62を透過し、これによりビームが発散して、発散ビーム64を
提供する。次いで、発散ビーム64はビームスプリッタ48によって反射されて
、第1のビーム66を提供する。第2のビーム68は、図示するように、ビーム
スプリッタ48を透過する。移動可能な照準ミラー54および60はそれぞれ、
ビーム66および68を反射して、ソースPおよびPを提供する。
【0035】 図2bを参照して、別の実施形態において、放射ソースPおよびPは、光
ファイバスプリッタ56を用いて、回転可能なレーザ40から放出された放射か
ら分割される。ファイバは、端部においてビーム形成要素を有して、2本のビー
ムの発散角度を制御または設定し得る。一実施形態において、ビーム形成要素は
レンズであり得る。ソースPおよびPは、あるいは、周波数が固定された1
対の回転可能なレーザから形成され得る。放射ソースの他の適切な実施形態は、
マイクロ波および音波などの制御可能な位相を有する波を生成する任意のソース
を含む。
【0036】 1実施形態において、小さい膨張係数を有する材料を含むバーの一方の端部に
放射PおよびPの各ソースを取り付けることによって、ソースは互いから一
定の距離Dで維持される。別の実施形態において、放射PおよびPのソース
は一定の距離で保持されないが、その代わりにこれらのソース間の距離Dは非常
に正確であることが公知である。
【0037】 互いに対して一定の距離Dで放射ソースPおよびPを支持するための例示
的な1つのバー70が図3に示される。バー70の両端部にはソケット74が提
供される。玉継手76は、示されているように、ソケット74の各々の内部で旋
回するように(pivotally)配置される。玉継手76の各々は(図2b
に示される)その中に配置される光ファイバースプリッタ56からのファイバー
の端部および発散放射が通過するアパーチャ80を有する。ファイバーはその端
部においてビームを形成する素子を有し得、2つのビームの発散角を制御または
設定し、1実施形態において、このビーム形成素子はレンズである。動作中、玉
継手76は、それぞれのソケット74の中で、矢印78によって示されるように
旋回可能であり、(図1に示される)制御ユニット32によって制御され得る。
この構成を用いて、ファイバーの端部におけるソースPおよびPによって提
供される発散ビーム12および14は、所望の方向に向けられ得、一定の隔離距
離(separation distance)Dが維持される一方で、処理さ
れるべき点Pを含む物体10の全体または部分を照射する。
【0038】 再び図1を参照して、物体10の表面上の点Pの座標は(x、y、z)であ
る。点Pのx座標およびy座標は、通常、検出器22および物体10の幾何学
的形状から直接的に決定されるが、焦点調節素子24をこれらの検出器22と物
体10との間に挿入して任意の倍率を考慮している場合には、深さ座標zを直接
求めることはできない。但し、ここでz軸は結像系の光軸と並行であると定義さ
れる。しかしながら、深さ座標zは、放射ソースPから物体10の表面上の点
および放射ソースPまでの光路差 s=R−R+S (1) をまず考察することによって測定され得る。量Sは、ビームが点PおよびP に到達する前に生じ得る、ビームにおける任意の光路長の差を考慮に入れるた
めに含まれる。
【0039】 sが非ゼロである場合、ソースPおよびPから発せられた放射の周波数を
変更することによって点P0で測定される1方のソースからの放射の位相が、他
方のソースに対して変化する。この位相変化によって、点Pにおける放射の強
度が変調する。強度の変化の1サイクルを完了するために必要とされる周波数Δ
νの変化は、以下の式によって与えられる。 Δν=c/s (2) ここで、cは光の速度である。従って、強度の1振動を引き起こすために必要と
されるレーザの周波数Δνの変化を測定することによって、光路差sが決定され
得る。zの測定は、その後、後述されるように、xおよびyの各値に対してsの
値を決定することに基づく。
【0040】 sのより正確な決定は、複数の振動サイクルにわたってΔνを測定することに
よって得られる。実際には、周波数の変化全体Bによって誘導される(必ずしも
自然数ではない)振動サイクル数Nに関して測定することが好都合である。
【0041】 NはΔνおよびBに関して、以下のように与えられる。 N=B/Δν (3) 式(2)を用いて、式(3)からΔνを消去すると、Nに関するsの以下の数式
が得られる。 s=(c/B)・N (4) 長さBの周波数走査により誘導された振動サイクル数Nは、周波数の走査中、
測定点Pにわたって移動する干渉縞の数に対応することに留意されたい。干渉
縞は、ソースPおよびPから発せられる放射の干渉によって生成される。こ
れらの振動は、結像レンズ24によって干渉縞が解消されてもされなくても生じ
る。
【0042】 Nの測定における不確定度ΔNは、 Δs=(c/B)・ΔN=s・(ΔN/N) (5) のsにおける不確定度Δsに対応する。
【0043】 式(5)は、単一の振動サイクルが決定され得る不確定度ΔNが一定の状態で
ある場合、sにおける不確定度Δsは、測定されたサイクルN数と等しいファク
タだけ減算されることを示す。当業者に公知の種々のレベルの分解能ΔNに対す
るNを決定する複数の方法が存在する。約1の振動サイクルカウント(ΔN=1
)の分解能を得る方法の例は、データシーケンスに高速フーリエ変換(FFT)
を行なうか、または高域フィルタリングされた信号のゼロ交差をカウントするこ
とである。1未満の振動サイクルカウント(ΔN<1)の分解能は、例えば、離
散フーリエ変換(DFT)の大きさが最大化されるところでDFTの引数(ar
gument)を見出すことによってか、または周波数走査の終了時において振
動サイクルの位相を検査することによって改善することができる。位相の精確な
検査に関する当業者に公知の1つの技術は、ビーム経路の1つのレッグ、すなわ
ちビームスプリッタまたは光ファイバースプリッタと、図2、図2(a)および
図2(b)におけるソースPまたはPのうちの1つとの間に位相変調器を挿
入することである。
【0044】 I、IおよびIがそれぞれ−90°、0°および90°の位相変調器に
よって生じる位相シフトに対応する信号強度である場合、振動サイクルの位相φ
は、以下の数式によって与えられる。 φ=tan−1((I−I)/(2I−I−I)) (
6) 位相における不確定度Δφは、2πで除算することによってサイクルカウント
における不確定度ΔNに変換される。
【0045】 チューナブルダイオードレーザのB=15THzの典型的な周波数走査の場合
、ΔN=1サイクルの不確定度では、Δs=20μmの不確定度が提供される。
横方向の分解能に対するsの広がりが、その量よりも小さい場合、ΔN=0.1
サイクルの不確定度は、sにおける不確定度をΔs=2.0μmまで向上する。
物体の表面上の横方向の分解能に対するsにおける広がりがΔsよりも大きい場
合、sの測定における分解能は改善されているので、この場合もやはり、その横
方向の分解能に対するsの平均値または代表値の推定を改善することができる。
【0046】 座標系に置き換えるとsは以下の式によって与えられる。
【0047】
【数1】 計算をより簡略化するために、2つのソースPおよびPが(x、y、z )および(−x、−y、−z)における原点のまわりに対称的に配置さ
れていると仮定する。その後、式(7)は、(x、y、z)に関して、以
下のようになる。
【0048】
【数2】 zについて解くと、式(8)は以下のようになる。
【0049】
【数3】 ここで、Dは2つのソースPとPとの間の距離である。従って、zは、式(
9)に正の根および負の根が存在するために、アンビギティの範囲内であること
が判定される。このアンビギティを回避する1つの方法は、物体10を照射する
ことにより、s=0線(s=0の場合である図1において16で示される)が
結像されるべき物体の領域を2等分しないようにすることである。s=0線を移
動する1つの方法は式(1)におけるsを変更することである。
【0050】 sの変更に対するシステムの感度は、Δs/Δzの比によって示され、ここで
、Δzはsの値における不確定度Δsによって導入されるzにおける不確定度で
ある。この比は、任意の実用的な感度の範囲が十分でないシステムの0と、理論
的には最大システムの2との間にわたる。2の値を達成するのは非現実的である
。なぜなら、物体10の表面は、2つの点ソースPとPとの間にある必要が
あり、各ビームによって表面の片面のみが照射されとるからである。比Δs/Δ
zは、zに対するsの偏導関数を得ることによって計算され、この結果より範囲
分解能に関する以下の数式が取得される。
【0051】
【数4】 式(10)において、図1に示されるように、Rは原点からPまでの距離で
あり、ψは原点から点Pに延びる線と、点Pから点Pに延びる線との間の
角である。良好な範囲の分解能を提供する有用な構成は、ψ=90°を設定する
ことであり、この場合Δzの数式は以下のように簡略化される。
【0052】
【数5】 ここで、θおよびαは図1に示されるとおりである。RおよびDに関しては
tanα=D/(2R)である。式(11)は、角αが増加し、角θが減少
するにつれて範囲分解能は向上することを示す。Δs=5μmの値、α=10°
およびθ=45°の場合、範囲分解能はΔz=20μmである。
【0053】 点Pの横方向の位置(x、y)における不確定度(Δx、Δy)も、範囲分
解能Δzに影響を与える。2つのソース点がx−z平面にある場合、zの測定は
、不確定度Δyの影響を受けない。ψ=90°に関して、xにおける不確定度Δ
xは、zの測定において以下の不確定度を生じさせる。
【0054】 Δz=Δxtanθ (12) 従って、θ=0°近傍の角であれば、点Pの横方向の位置における不確定度
への影響は最も少なくなる。
【0055】 焦点の深さは、最適なシステムの横方向の分解能が向上すると減少するので、
横方向の分解能と物体の最大深さとの間にトレードオフが存在する。物体の深さ
におけるこの制限を低減するための1つの方法は、異なった範囲の平面に連続し
て焦点を合わせ、焦点の深さの中にあるピクセルのみを用いることである。例え
ば、100μmの横方向の分解能は、フィールドの深さを1cmのオーダーに制
限し、10cmの範囲を有する物体は10個の異なった範囲で連続的に焦点を合
わせることによって全分解能で結像され得る。フィールドの深さの影響を最小化
するために、z軸は、範囲の広がりを最小化する、すなわち、物体の表面の平均
的平面に対して垂直の方向に規定され得る。横方向の分解能を損なわずに、結像
され得る横方向の面積を増加させるために、複数のカメラ(すなわち、検出器ア
レイ22)を用いて物体10の対象となる領域全体を覆てもよいし、または対象
となる個別のカメラを用いて検査してもよい。あるいは、単一のレンズの焦点面
に複数の検出器アレイを配置してもよい。これらのアレイは、独立して平行移動
され得、物体の種々の領域を高分解能で検査する。z軸に沿う個別の検出器アレ
イの平行移動または検出器アレイの傾動(tilting)は、異なった深さの
物体の領域に対して焦点を同時に合わせ、物体深さの受容可能な範囲を大きくす
ることができる。
【0056】 図1における光学結像システムの潜在的問題は、ソースと検出器との間のバイ
スタティック角によってシャドウイング効果(shadowing effec
ts)が発生し得ることである。これらの効果は、図4におけるように、レンズ
をソースにより近づけて配置し、かつ検出器が結像面において横方向にオフセッ
トされる、レンズを軸から外した構成で用いることによって低減され得る。この
目的のためにレンズが設計される場合、またはレンズが十分に大きい視野を有す
る場合、軸から外して結像することによって生じる収差は最小化され得る。
【0057】 図5を参照して、本発明の代替的実施形態は可動の放射ソースPおよび静止
放射ソースPを含み、各ソースは発散ビーム150および154を提供し、そ
れらの放射ソースと、物体10の表面上の点Pとの間にRおよびRで表示
された光路長をそれぞれ有する。ソースPおよびPは、単色レーザ等のコヒ
ーレント光の任意の適切なソースによって生成され得、このコヒーレント光は、
2つの点ソースPおよびPを提供するために分割され得る。さらに、図2お
よび図2aのビームスプリッタの実施形態、ならびに図2bの光ファイバースプ
リッタの実施形態等の種々の技術が、コヒーレント放射ソースからの放射を分割
するために適切である。
【0058】 発散ビーム150および154は、位置情報を有する点Pが配置され、測定
されるべき物体10の表面の方向に向けられる。物体10の表面によって散乱さ
れる照射は、焦点調整素子またはレンズ158によって焦点合わせされ、検出器
アレイ22に当たる。レンズは、図4に図示されるように、軸から外された構成
で用いられ得、バイスタティック角に起因するでシャドウイング効果を低減する
。図1と関連付けられた上述のタイプの任意の偏光子(図示せず)は、焦点調整
素子158と検出器アレイ22との間に配置され得、検出器アレイ22に入射す
るスペックル画像または縞パターンのコントラストを向上させる。
【0059】 検出器アレイ22は、後述されるように、検出器に入射する画像を処理するた
めのプロセッサユニット28と通信する。制御ユニット32は、ソースPを軸
160に沿って移動させるために、少なくとも可動ソースPと通信する。上述
のように、制御ユニット32およびプロセッサユニット28は別個のデバイスに
よって実現され得るか、または代替的には、単一システムの部分であり得る。さ
らに、制御ユニット32およびプロセッサユニット28は、特定の用途において
所望され得るように、互いに通信し得る。
【0060】 図1と関連付けて上述されたように、物体10の表面上の点Pと関連付けら
れた深さ座標zは、ソースPおよびPのそれぞれから点Pへのビーム15
0の光路長Rと、ビーム154の光路長Rとの間の差の関数R−Rとし
て決定され得る。図5の実施形態において、可動ソースPからの放射の位相は
、制御ユニット32の制御のもとで、ソースPを軸160に沿って移動させる
ことによって変更される。この構成を用いて、点Pで強度における振動が生成
される。
【0061】 可動点ソースPの瞬間座標は以下のとおりである。
【0062】 x1=al,y1=amおよびz1=an (13
) ここで、aは点ソースPの平行移動の大きさを表し、l、mおよびn
x、yおよびz軸それぞれに関する平行移動の方向を表す方向余弦である。 これらソースPおよびPからの放射が点Pへ伝播した後に測定された、の
位相差は、以下の数式によって与えられる。
【0063】
【数6】 ここで、φは、2つのコヒーレントソースPとPとの間に存在し得る一定
の位相オフセットを表す。Pが軸160に沿って平行移動すると、Rの値は
変化し、φをaの関数として変化させる。
【0064】 ソースPが原点から移動して離れると、点Pにおいて生じる振動の強度(
またはPを横切る干渉縞)の数は、以下の数式によって与えられる。
【0065】
【数7】 ここで、Rは、点Pと座標系の原点との間の距離であり、φ(a)は、ソー
ス間の距離がaの場合の式(14)における光学位相差であり、φ(0)は、a
=0の場合の式(14)における光学位相差である。式(15)を考察すること
によって、ソースPからの移動の結果である振動の強度の数Nは静止ソースP の位置には依存しないことが明らかとなった。この非依存性は、ソースP
よびPが互いに近接して配置されることを可能にする。この構成を用いて、ソ
ースPおよびPそれぞれからの発散ビーム150および154は、大気乱流
および振動等の一般的な妨害を受ける。このようにして、このような妨害の影響
は最小化される。さらに、ビーム150および154は、実質的に同じ偏光を有
した状態で物体10の表面に達する。
【0066】 ソースPとPの間の距離が小さいと、物体10の表面上に縞の間隔が大き
い縞パターンを生成するが、これは結像レンズ24によって容易に解消され得る
ことにも留意されたい。
【0067】 点ソースPの平行移動aの大きさは、Rの値と比較して、比較的小さいの
で、式(15)は、以下のようにa/Rにおける二次まで近似され得る。
【0068】
【数8】 ここで、ψは、原点から点Pへと延びる線と、Pの平行移動の方向によって
規定される線との間の角である。
【0069】 式(16)は、a/Rにおける最下位まであれば、Nが分かるので角ψを決
定することができるということを示す。3つ以上の位置からψが分かれば、P の(x、y、z)座標は三角形分割によって決定され得る。次に、図1に対応す
る実施形態と類似の実施形態が説明される。ここで、x座標およびy座標は、検
出器アレイにおける結像点の位置から決定される。
【0070】 所与の(x、y)の位置に対するzは、Pのある距離の移動に対応する強度
の振動サイクルNを決定することによってか、またはそのような強度振動が生じ
る速度を測定することによって測定され得る。まず、サイクル数Nを決定するこ
とに基づくzの測定を考察されたい。既知のNを用いれば、zを除く式(15)
におけるすべての変数は分かる。zについて式(15)を解くことで、以下の数
式が得られる。
【0071】
【数9】 式(19)は、0と1との間で変化する大きさを有する無次元パラメータを定義
する。これは、Pによって移動される波長単位ごとの振動サイクルNに関する
スペックル強度の平均変調率を表す。aの値が0に近づくと、式(17)は以下
のように近似され得る。
【0072】
【数10】 式17および式20におけるzの数式は、n=0を設定することによって簡
略化され得、その結果、ソースPの平行移動はx−y平面に限定される。この
構成は、後述されるように、ソースPの平行移動に関する良好な実用的選択を
表す。結果として生じるzの数式は、以下のように表され得る。
【0073】
【数11】 ここで、散乱点Pからx、y座標系の原点への距離Rは、完全な数式によっ
て与えられる。
【0074】
【数12】 aが小さいとき、Rは以下のように近似され得る。
【0075】
【数13】 次に、強度の振動が生じる瞬間比を知ることによってzの測定を考察されたい
。瞬間振動率ρは、式(19)における平均振動率と類似の方法で以下のように
表され得る。
【0076】
【数14】 式(15)から得られる強度の振動の数Nの数式を式(24)に代入すること
で以下の式が得られる。
【0077】
【数15】 ここで、関係を分子を簡略化するために用いる。
【0078】 l +m +n =1 (2
6) aの値が小さい場合、ρは以下のように近似される。
【0079】
【数16】 式(25)をzについて解くことで以下の式が得られる。
【0080】
【数17】 ここでは、 O=[xl+ym+a(ρ−1)]n (29)
および
【0081】
【数18】 である。n=0のとき、式(29)は、式(21)の形式で表され得る。
【0082】
【数19】 aの値が小さいとき、式(28)および(31)は、式(20)および式(23
)によってそれぞれ近似され得る。ただし、ρの代わりにρ\を用いる。本明細
書中において「\」とは、「\」の前の文字に上線を付していることと等価であ
り、例えば、ρ\は、
【0083】
【数20】 を示す。
【0084】 範囲分解能を推定するために、測定される量の不確定度(ΔNまたはΔρ)に
よって生じるzの測定における不確定度Δzを考察されたい。簡略化するために
、この計算は、式(16)によって与えられるNに関する近似式に基づく。Δz
を求めるために、zについてN(またはρ)の偏導関数をとり、この偏導関数を
比ΔN/Δz(またはΔρ/Δz)と等価すると、以下の式が得られる。
【0085】
【数21】 は、平行移動の方向および散乱点への方向を考慮に入れたジオメトリカルファク
タである。Gに関する第1の式において、 l=x/R,m=y/Rおよびn=x/R (34
) は、点Pの方向余弦である。Gに関する第2の式において、θは極角およびφ
は方位角であり、球面座標系における原点からPへの方向を表す。同様に、ソ
ース点の平行移動の方向はθおよびφによって与えられる。
【0086】 式(32)の考察によって、物体の距離Rが増加するにつれ、範囲分解能が
低下し、ソースPの平行移動aの大きさが増加するにつれて、範囲分解能が向
上することが明らかとなる。式(33)の考察は、ジオメトリカルファクタGが
1と無限との間で変化することが分かる。ここで、1は達成可能な最良範囲分解
能に対応する。
【0087】 所与の散乱点の方向に対するソースPの平行移動の最適な方向は、Gがl、
mおよびnの所与の値に対して最小化されるようにl、mおよびnを選択
することによって式(33)から得られる。この制約を適用すると以下の式が得
られる。
【0088】
【数22】 これは、最適な平行移動の方向が、原点から散乱点P(ψ=90°)へと延び
る線と直交し、その線およびz軸(φ=φ)によって形成される入射面にある
ことを意味する。式(35)の値を式(33)に代入すると、以下のようになる
【0089】
【数23】 式(36)から、1の達成可能な最良のG値は、n=0(θ=90°)のときに
生じることが分かる。これは、散乱点がx−y平面にあることを意味する。さら
に、z軸上にある分散点に対してGが無限大に近づくにつれ分解能が低下するこ
とも分かる。例えば、θ=30°の場合はG=2、およびθ=10°の場合はG
=5.76である。各点に関する平行移動の方向を変更することなく、画像内の
すべての点について式(35)を満たすことは可能ではないが、最適分解能の条
件は代表的な画像点について式(35)を満たすことによって近似され得る。
【0090】 式(25)および式(27)によって、瞬間変調率ρは、平行移動する点のオ
フセットの大きさaに依存する。測定ρに基づく技術の場合、走査中にρができ
る限り変化せず、その結果、値ρとzとは、ほぼ1対1で対応することが所望さ
れる。従って、標準的スペクトル解析技術は、ρの値を推定しzを決定するため
に適用され得る。走査時に生じるρの不均一性の程度を定量化するために、以下
のように定義される。
【0091】
【数24】 式(27)からのρに関する近似式を式(37)に代入し、aを含む最低次数項
のみを残すと、以下の式が得られる。
【0092】
【数25】 方程式(38)は、変調の不均一性が走査の長さ対物体の距離の比a/Rが線
形に増加することを示す。さらに、不均一性は、ψ=0°であるときに消失し、
ψ=90°であるときに有界を有することなく増加する。しかしながら、ψ=0
°の線におけるすべての点は、範囲、すなわち式(33)におけるG=∞に関係
なく、同じ変調率を有するので、ψ=0°のときに範囲分解能は存在しないこと
が分かる。従って、不均一性を最小化することと、最適範囲分解能との間にトレ
ードオフが存在する。
【0093】 良好な範囲分解能および低減された変調の不均一性を同時に提供する好都合な
測定構成は、n=0を設定し、φ方向におけるオフセット、すなわちφ=φ を有する軸から外れた光学系を用いることである。その後、Gに関する式(3
3)は以下の式に通分される。
【0094】 G=−2/sin(2θ) (39) 測定技術の例示的例として、x−y平面における200mm×200mmであ
る物体を、R=1mの距離から波長λ=0.7μmのレーザを用いて撮像する
ことが所望されることを仮定されたい。n=0および物体の中心がθ=30°
およびφ=φにおいて配置される場合、式(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%よりも少ない。この不均一度は
、走査中に、走査速度を少しだけ変化させ、測定中の周波数の任意の変化を補償
することによってさらに低減され得る。
【0095】 図6は、図1および図5のプロセッサ28によって行われる例示的な一連の工
程を示す。これらの工程を行って、オブジェクト上の各ポイント(x、y)にお
ける深さ座標zを判定する。プロセッサは、工程100から開始し、オブジェク
ト表面上の複数の照射ポイントによって散乱される放射の強度のパラメータを測
定する(工程108)。この情報から、各測定ポイントのz座標を計算する(工
程112)。
【0096】 工程116において、光学フィルタリングプロセスを行うことができる。当業
者に公知の適切なフィルタを挙げると、平滑化フィルタ、メジアンフィルタおよ
びカーブフィッティングフィルタがある(ただし、これらに限定されない)。そ
の後、マッピングされたポイントの表示または出力を当業者に公知の何らかの方
法で行うことができ、その後、本プロセスは図示のように工程124において終
了する。一実施形態において、工程120において、マッピングされたポイント
を、メッシュプロット上の演算されたz情報の関数としてプロットする。
【0097】 図6aも参照して、工程108および工程112の一実施形態が図示されてお
り、これらの工程は、図1の実施形態と共に用いると適切である。工程108’
において、散乱される照射の強度をレーザ周波数オフセットの関数として測定し
、当業者に公知の方法の1つを用いてNを測定する。その後、工程110’にお
いて各ロケーション(x、y)について式(4)を用いてsを計算し、工程11
2’において各ロケーション(x、y)について式(9)を用いてzを計算する
【0098】 図5の実施形態と共に用いられる処理工程108および処理工程112の別の
実施形態を図6bに示す。この場合、工程108’’において測定される強度の
パラメータは、回数N(これは必ずしも整数ではない)であり、この回数の間、
可動ソースP(図5)が平行移動するにつれ、強度は周期的に変化する(cy
cle)。当業者に公知の方法を通じて工程108’’においてNが判定される
と、工程110’’において、このNは式(19)によってρ\に変換される。
その後、工程112’’において式(17)および式(18)を用いてzを計算
する。図5の実施形態と共に用いられる処理工程108および112の別の実施
形態を図6Cに示す。この図6Cの実施形態では、工程108’’において測定
された強度のパラメータは、ソースポイントPが平行移動するときに振動が発
生したときの瞬間振動速度ρである。工程112’’’において、式(28)〜
(30)を通じてρをzに変換する。
【0099】 検出器22およびプロセッサ28は様々に構成することが可能である。図7に
示す一実施形態において、検出器アレイ22の光検出器素子221、1〜22 、m を逐次読み出す。検出器アレイ22のシリアル出力36は、プロセッサ28
への入力を提供する。プロセッサ28は、単一のプロセッサを含み得るか、ある
いは、複数のプロセッサを含むマルチプロセッサであり得る。
【0100】 図7aも参照して、検出器およびプロセッサの別の構成を示す。この実施形態
において、プロセッサ28は、複数のプロセッサ281、1〜28n、mを含む
マルチプロセッサである。検出器アレイ22の光検出器素子221、1〜22 、m はそれぞれ、プロセッサ281、1〜28n、mのうち対応する1つに各出
力信号38を提供する。この構成を用いると、プロセッサ281、1〜28n、 はそれぞれ実質的に同時に動作することができ、有利な性能を実質的に提供す
る。より詳細には、マルチプロセッサユニット28内の各プロセッサ281、1 〜28n、mは、光検出器アレイ22の対応する素子221、1〜22n、m
ら受信したデータに基づいてz座標を計算する機能をする。従って、オブジェク
ト10の表面の各ロケーションのz座標を敏速に判定することができる。
【0101】 図7bは、検出器およびプロセッサコンポーネントのさらに別の実施形態を示
し、これらの検出器およびプロセッサコンポーネントは、図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、 は、検出器221、1〜22n、mそれぞれからの情報を実質的に同時に処理
して、深さ座標を判定する。
【0102】 図8を参照して、本発明の別の実施形態は、検出器22’のアレイを含み、こ
のアレイはオブジェクト10の表面に対向して配置されており、オブジェクト1
0の表面輪郭が測定される。この構成では、オブジェクト10の表面上のポイン
トPからの散乱光を観察してzを判定するのではなく、光の位相シフトの測定
をオブジェクトの表面において直接的に行う。図示はしていないものの、図8の
システムは、ソースPおよびPを制御する制御ユニット28と、検出器22
’に入射する放射を処理するプロセッサ28とを含む。検出器22’については
、図8中に図示しており、図1、5および6と関連して説明している。
【0103】 オブジェクト10の表面上に光検出器素子23を配置する構成およびメカニズ
ムには様々なものがあり得る。図9に示す一実施形態において、アレイ22’の
複数の個々の光検出器素子23を、対象エリア内のオブジェクト10表面上に配
置する。
【0104】 図9aに示す別の実施形態において、支持/制御ユニット88の片側から出て
いるスプリングアーム84上にアレイ22’の個々の光検出器素子23を取り付
ける。これらのスプリングアーム84を制御ユニット88によってオブジェクト
10の表面上で移動させて、特定のポイントまたは対象領域と接触させる。この
ようにスプリングアーム84を片持ち型で支持すると、アーム84がオブジェク
ト10上を移動している間、個々の検出器23はそれぞれ、オブジェクト10の
表面上のロケーションとの接触状態を保つ。すなわち、オブジェクト表面の輪郭
が変化すると、それに応じてスプリングアーム84も上下に移動する。
【0105】 本発明の装置および方法において、放射ソースをさらに2つ以上用いてもよい
ことが理解される。例えば、さらなるソース(単数または複数)を用いて、オブ
ジェクトまたはその一部分に関するx、y座標情報を判定することができる。さ
らに、さらなる放射ソースを用いて、対象領域のシャドーイングに寄与するいか
なる処理精度の不足または処理のあいまい性をも低減することができる。
【0106】 ソースポイントの移動を行う実施形態は他にも変更が可能であることが理解さ
れる。例えば、これらの2つのポイントはどちらとも、反対方向に移動しながら
移動することもできるし、一定の間隔を保ちつつ同じ方向に移動することもでき
るし、共通の中心点の周囲を回転することもできるし、あるいは、制御システム
による切換えが可能なソースポイントのアレイを用いることによって動きをシミ
ュレートすることも可能である。
【0107】 双方のソースを移動させると、単一の可動ソースを用いる場合と比較して式(
38)における変調の不均一性を低減できるなどの利点を得ることができる。こ
の不均一性を低減させることによって企図しているのは、Nがaの奇関数となる
ような様式でポイントを移動させることにより、式(27)中のrに関する線形
項を無くすことである。平行移動を2つのソース間で均等に分割して、これらの
2つのソース間の中心が固定されるようにする(これは、対向移動(oppos
ing motion)と呼ばれる)。その結果、以下の数式が得られる。
【0108】
【数26】 a中に4次までの項を含む式(16)に対応するNを展開すると、以下のよう
になる。
【0109】
【数27】 従って、a中の2次項およびそれよりも高次のa中の偶数項を無くし、3次項
の大きさを1/4だけ低減する。
【0110】 これにより、式(17)および式(22)に基づいてzを正確に数式として表
すと、以下のようになる。
【0111】
【数28】 =0の場合、式(42)は式(21)まで低減し、以下のようになる。
【0112】
【数29】 aの値が小さい場合、式(42)および式(43)はそれぞれ、式(20)お
よび式(23)となる。しかし、双方のソースを移動させることによってaへの
1次依存(first−order dependence)を無くすと、移動
ソースを1つだけ用いる場合よりもこの近似結果の精度が上がる。
【0113】 2つのソースを対向移動させると、式(25)によって表される変調速度の数
式は以下のようになる。
【0114】
【数30】 小規模のaの近似(これは、aまで有効であり、3次のaを含む)を用いると
、式(44)は以下のようになる。
【0115】
【数31】 式(45)中に線形項が存在しないと、スキャン期間中の平均振動速度ρが式
(27)と比較してより安定する。ここで、式(37)において既に規定された
対応する変調不均一性χを以下のように表すことができる。
【0116】
【数32】 式(46)の場合、a/Rの比への2次依存および共通因子中のcosyの
不在のため、得られるχの値は、式(38)のχ値よりもずっと小さい。双方の
ソースポイントを移動させることによって得られる利点の1つとして、レンジ分
解能を最適にした(ここでy=90°)場合、不均一性が小さいままであるとい
う点がある。
【0117】 単一のソースの動きについて既に与えられたレンジ不確定性Δzの結果は、N
およびρの小a近似(small−a approximation)における
第1の項に基づく。これらの近似における第1の項は2つのソースの動きと同一
であるため、先に述べた結果は、2つのソースが移動する場合にも当てはまる。
【0118】 図10を参照して、本発明の別の実施形態は、2つの対向する可動ソースP およびPを含む。レーザソース40からのビーム44は、ビームスプリッタ4
8によって送信ビーム46および反射ビーム50に分割される。送信ビーム46
は、鏡94からレンズ58へと反射される。同様に、反射ビーム50も鏡96か
らレンズ52へと反射される。鏡付き表面91および93を備えた直角プリズム
92を用いて、ビーム46および50をそれぞれレンズ58および52から再度
方向付ける。ビーム46および50はそれぞれ、ポイントPおよびPに集束
し、発散ビーム82および86として延びる。位相シフトは、別個のビーム46
および50がたどるビーム経路間の光路長の差を変更することによって達成する
ことができる。これは、例えば、素子48、94または96のうち少なくとも1
つを圧電変換器(図示せず)によって面外平行移動させるようにすることにより
、達成可能である。あるいは、位相シフトデバイス(図示せず)を、ビーム経路
46または50のうち1つに挿入してもよい。このようなデバイスの一例として
は回転するガラス製ディスクがあり、このガラス製ディスクは、別個の工程にお
いて回転角度を用いて自身の光学的厚さを変化できるようにコーティングされて
いる。プリズム92は軸99に沿って平行移動することが可能であり、これによ
り、ソースPおよびPの対称な動きに影響を与える。あるいは、レーザ40
、ビームスプリッタ48、鏡94および96ならびにレンズ52および58が1
つの群として軸99に沿って平行移動する間にプリズム92を固定状態にして、
ソースPおよびPを対称に平行移動させてもよい。ソースPとソースP との間の距離の変更は、レーザ40とビームスプリッタ48との間のビーム44
中に配置された角度変調デバイス(図示せず)を用いても達成することが可能で
ある。このような角度変調デバイスは、ビーム角度を変化させる。角度変調デバ
イスの例を挙げると、検流計鏡および音響光学変調器がある。
【0119】 別の実施形態(図示せず)において、ソースPおよびソースPは仮想ソー
スである。すなわち、ソースPおよびソースPは、ビーム46および50が
集束状態で通過する場所である光または領域の小ソースではない。その代わりに
、ビーム46および50を発散させる光コンポーネント(すなわち、負レンズ)
を用いて発散ビーム82および86を生成する。その結果、オブジェクトへの照
射を、実際のソースPおよびPを用いた他の実施形態から区別不可能にする
ことが可能となる。
【0120】 図11は、反対方向に移動するソースPおよびPを生成するさらに別の構
成を示す。レーザソース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)放物線型鏡セグメントにして、PおよびPにおいて小さなフォーカルス
ポットを生成させてもよい。あるいは、曲面鏡95、97を球状の素子にしても
よい。球状の素子によって生じる光学収差は、縞接触には影響を与えず、縞位置
のみに影響を与える。これらの光学収差は、解析に取り入れることが可能である
。位相シフトは、ビームスプリッタ48を平行移動させるか、または、上述した
ようなビーム経路46もしくは50のうち1つの中に位相シフトデバイス(図示
せず)を配置することにより、達成可能である。
【0121】 ここで再度図1を少し参照して、上記における導出は、テレセントリックイメ
ージング構成に基づく。焦点面中のイメージポイントPの側方位置は、レンズ
24の倍率によってオブジェクトポイントPの側方位置に関連付けられ、レン
ズ24とオブジェクト10との間の距離から独立する。焦点はずれによってイメ
ージPは不鮮明となるが、これは、検出器22においてイメージPがはっき
りと側方位置にある状態には影響を与えない。そのため、オブジェクトポイント
の側方位置をそのイメージポイントPの位置によって推定する作業をP とPとの間の距離に関係無くおこなうことができる。残念なことに、オブジェ
クトが大型である場合、テレセントリックイメージャ24は用いることはできな
いかもしれない。その理由は、レンズ24の直径は、少なくとも測定対象オブジ
ェクト10の最大寸法と同じくらいの大きさでなければならないからである。従
って、ほとんどのレンズはテレセントリックではない。
【0122】 ここで図12を参照して、この問題を回避するためのポイントPの(x、y
、z)座標を判定するための一般的なアプローチとして、検出器22におけるイ
メージポイントPの位置は、オブジェクトポイントPとフォーカシング素子
24の光学軸98との間の角度αを表すものであると仮定するアプローチがあ
る。座標は、縞数Nと、検出器22におけるイメージポイントのロケーションと
から導出することができる。以下の方程式は正確である(小aは仮定していない
)。 側方座標は、以下によって得られる。
【0123】
【数33】 ここで、αおよびαはそれぞれ、x−z面およびy−z面において測定され
た角度の成分である。 z座標は、以下によって得られる。
【0124】
【数34】 座標(x、y、z)は、ポイントPとPとの間の中間位置Pを表
し、(x、y、z)が座標系の起点以外の場所にある場合に一般化された
解法を提供する。多用性を広げるために、Pおよび他のパラメータを経時変化
させてもよい。
【0125】 上記にて述べた一般化されたオブジェクト座標の判定は、複数の検出器22(
例えば、CCDアレイ)と連絡するプロセッサユニットからのデータを受信する
コンピュータ上のソフトウェアにおいて実施することが可能である。あるいは、
プロセッサユニット28上にアルゴリズムを直接インプリメントしてもよい。
【0126】 式(47)〜(49)を用いてポイントPの(x、y、z)座標を判定する
ためには、角度αおよびαだけではなく、正規化された縞数ρ\(例えば、
式(19)によって規定された縞数ρ\)もわかっていなければならない。上記
の実施例において、これらの角度は、焦点面22内のイメージポイントPのロ
ケーションによって判定している。αおよびαを正確に判定するためには、
レンズ24内に発生する可能性のある任意のイメージひずみを修正することが必
要である。ひずみの修正は、例えば、角度αおよびαをイメージポイントP 上にマッピングするひずみ機能を生成することにより、行うことが可能である
。ひずみマップは、レンズが変わると変動し、また、同じレンズ24を異なるイ
メージング構成において用いた場合にもわずかに変動する可能性がある。
【0127】 本発明の別の実施形態において、さらなるソースヘッドを用いて三角測量に必
要な残りの2つの方向成分を判定することにより、レンズひずみを特徴付けると
いう不便な作業を無くす。例えば、図12aに示すように、さらなるソースヘッ
ドを位置Pに配置し、以前Pにあった(図12を参照)レンズ24を、オブ
ジェクト10上の対象領域をビューイングする際に便利の良い任意の位置に配置
する。第2のソースヘッド内のソースP1xおよびP2xがx軸に平行に方向付
けられている場合、このソースヘッドのポイントPにおける縞数Nを測定す
ることは、ポイントP1xおよびポイントPならびにポイントPおよびP (式(41)を参照)によって規定されたライン間の角度yを測定することに
相当する。その後、式(47)〜(49)に必要な角度αを、α=π/2−
yの関係を用いてyから判定する。同様に、(図12aでは図示していないが
)ソースP1yおよびP2yをy軸に平行に方向付けた状態で位置Pにおいて
ソースヘッドを用いることにより、角度αを判定することも可能である。ソー
スヘッドをこのように方向付けた場合、縞数Nから、α=π/2−yを通
じて角度αを測定することができる。角度αおよびαを提供するソースヘ
ッドは、例えば、z軸周囲を90°回転する同一のソースヘッドであるか、また
は、ビームスプリッタを用いることによって同一のロケーションに現れるように
された別のソースヘッドであり得る。
【0128】 三角測量を複数のソースヘッドを用いて行う場合、レンズ24および検出器2
2のロケーション(またはオブジェクト10上に配置された検出器23のロケー
ション成分)を知る必要は無い。イメージングシステム22および24(または
検出器23)の目的は、Pにおける照射強度を観測することである。検出器ア
レイ22中の各ピクセルは、観測ポイントPを表す。そのため、各ピクセルに
ついて(x、y、z)座標を判定することができる。様々なアプローチを用いて
、異なるソースヘッドに対応するNの測定を別個に行うことが可能である。例え
ば、動作波長が異なる複数のソースヘッドを同時に照射することによって多重化
された異なるソースヘッドまたは波長を逐次照射することにより、測定値を時分
割多重化(time multiplexed)することが可能である。偏波ダ
イバーシチを用いることにより、測定値をさらに区別することも可能である。
【0129】 複数のソースヘッドを用いた三角測量プロセスをさらに説明するため、ここで
式(41)を再度参照する。典型的な用途の場合、αはRよりも小さく、式(
41)は以下のように簡約することができる。
【0130】
【数35】 正規化された縞数ρ\について式(19)を用いてNを書くことにより、ソー
スの間隔aおよび波長λまたは周波数νへの依存を式(57)から無くす。
【0131】
【数36】 従って、特定のソースヘッドに対応する正規化縞数は、ソースヘッドの中間点P から観測ポイントPへの方向の成分yを直接的に示す。Pの三次元座標は
、3つ以上のソースヘッドからの三角測量によって得られた測定値を通じて判定
することが可能である。
【0132】 ソースヘッドが3つの位置Pm1、Pm2およびPm3に配置され、これらの
ロケーションそれぞれに対応するソース対の方向は、方向コサイン(ls1、m s1 、ns1)、(ls2、ms2、ns2)および(ls3、ms3、ns3 )によって得られると仮定する。デカルト座標および方向コサインについて見る
と、式(58)は以下のようになる。
【0133】
【数37】 ここで、下付き文字は、ソースヘッドを区別する。P座標(x、y、z)の
判定は、各ソースヘッドについて正規化縞数ρ\、ρ\およびρ\を測定
し、式(59)によって表される3元連立方程式を解くことにより、達成される
。式(59)の形式主義は、一般化されたソースヘッドを配置および方向付けを
コンパクトな形態で取り扱うため、便利である。
【0134】 特定のピクセル22に対応するオブジェクトポイントPが3つ以上のソース
ヘッドによって照射されない場合、周囲のピクセル22について計算された方向
間で補間を行うことにより、当該ピクセル22の方向情報の欠損を補うことがで
きる。さらに、複数のソースヘッドから得られた方向情報をピクセルロケーショ
ンから得られた方向情報と組み合わせて、測定品質を向上させることもできる。
【0135】 複数のソースヘッドを用いれば、受信器パラメータ(例えば、位置、方向付け
、倍率およびひずみ)が分からなくてもPを判定することが可能であるため、
複数のイメージャ22および24(分かりやすくするため、イメージャは1つし
か図示していない)からの情報を組み合わせる作業が容易になる。これが可能な
のは、各カメラについて生成された点群データ(point cloud)は既
に同じ座標系にあるため、較正プロシージャを行ってデータを登録しなくてもよ
いからである。複数のイメージャ22および24を配置して、特定の対象領域を
任意の所望のビューイング角度および方位分解能でビューイングすることが可能
である。そのため、1つのイメージャ22および24によってオブジェクト10
全体の概要を得て、他のイメージャ22および24によって対象領域を高分解能
で検査することが可能である。
【0136】 上記の主な種類の実施形態の1つにおいて、ソースPおよびPの1つまた
は両方の動きによって生じる強度−振動をカウントすることにより、オブジェク
ト10の表面上のポイントPにおける縞数Nを判定する。サイクルのカウント
が最も基本的な様態で行われる場合、カウントを追跡するためには、1サイクル
あたり少なくとも2回強度をサンプリングしなければならない。また、コストを
最小限にするためには、標準的な市販のコンポーネント(例えば、フレームグラ
ッバーに接続された電荷結合素子(CCD))を用いてサイクルをカウントする
と好ましい場合がある。カウント数が延びるほど、イメージフレーム中のポイン
トPはN=0の位置から来るため、このアプローチを用いると、イメージフレ
ームを何百単位で取得および処理しなければならなくなる可能性がある。
【0137】 Nを判定するために多くのイメージフレームを取得する必要は、より高度なカ
ウント技術を用いることによって解消することが可能である。この技術は累進縞
分割(progressive fringe division)と呼ばれ、
これは、強度変調の周期的性質および予測可能な性質を利用したものである。そ
のため、ソース間隔aの複数の別個の値における強度変調をサンプリングするこ
とにより、Nを判定することが可能である。
【0138】 図13に示すように、累進縞分割は、別個の工程において縞サイズを初期粗大
間隔(Col.1)から最終微細間隔(Col.3)まで変化させる工程からな
る。各工程において、ラップされた(wrapped)位相マップまたはサイク
ルマップ(ロウ(ROW)4)を位相シフトした縞パターン(ロウ1−3)から
生成する。ブートストラッピングオペレーションにおける先行するアンラップの
位相マップに基づいて、連続するラップされた位相マップをそれぞれアンラップ
状態にする。工程間の縞間隔の比は典型的には、10以上のオーダーであるため
、少数のデータフレームを用いて極めて高密度の縞パターンをアンラップ状態に
することができる。高密度の縞パターンから位置座標を計算すると、精度が向上
する。累進縞分割を行うことによって得られる1つの利点は、オブジェクト10
の表面の複雑性に関係無く、アルゴリズムによって縞の順序の番号を各縞に正確
に割り当てるという点である。例えば、オブジェクト10は、表面の不連続部分
、大型の断崖形状部分および穴部を含み得、または、複数のばらばらのオブジェ
クトからなり得る。累進縞分割は表面の連続性に依存しないため、累進縞分割を
用いて、空間中に分散した未接続の複数の別個のポイントの位置をロケートする
ことが可能である。
【0139】 累進縞分割およびその実施形態についてさらに説明するにあたり、以下の3つ
の事実を考えると有用である(簡潔にするため、以下の議論において、ソースP およびソースPは対向移動すると仮定するが、以下の結果は、他の状況(例
えば、1つのソースPのみが移動し、静止状態の周波数がチューニングされた
ソースを用いる)にも当てはまる)。第1に、ソース間隔が初期値であるゼロか
ら最終値aまで増加するときにポイントPにおいて発生する強度振動Nの数は
、最終値aのPにおける縞数Nに等しい。例えば、10が単純なオブジェクト
(例えば、面)であった場合、N=0の位置(これは、R=Rにおいて発生
する)とポイントPとの間のオブジェクトの表面上のN個の縞をカウントする
ことができる。Nの値は、部分的サイクルまたは部分的縞を含むため、必ずしも
整数ではない。
【0140】 この第1の事実については、ソース間隔aがゼロからその最終値まで増加する
ときに縞がポイントPを通過してPとN=0との間の領域内部に集中する縞
の数としてNを解釈することも可能である点に着目することにより、この第1の
事実をさらに説明することができる。言い換えると、a=0である場合、この領
域内に縞は無く、この領域内部に含まれる縞の数は、縞がPを通過するたびに
1つずつ増加する。このように縞数Nと振動サイクルとが対応している状態は、
ロバストな三次元イメージング技術への鍵であり、これにより、ポイントP
囲の表面の連続性について仮定を行う必要なく、任意のポイントPにおけるN
を判定することが可能となる。これにより、表面中の不連続部分(例えば、穴部
および断崖形状部分(すなわち、エッジ))をあいまいさ無く処理する。
【0141】 第2の事実とは、式57によれば、Nは間隔aに直接比例するということであ
る。その結果、aの実行(running)をゼロからその最終値まで行ってN
(a)を判定して全体的スキャンを完了する必要が無くなる。間隔がaである
場合のNの値は、以下の関係を用いて、間隔aにおける値から推定することが
可能である。
【0142】
【数38】 上記の式が示すのは、2つのソース間隔aおよびaにおけるNの値が分かっ
ていると、それらの値を用いて、より小さな差間隔a−aにおけるNの値を
推定することが可能であるということである。
【0143】 第3の事実とは、いかなる中間値のソース間隔aにおいても位相シフト測定を
行うことを可能にすると、視野内の各測定ポイントPにおいて、部分的サイク
ルまたは「ラップされた」コンポーネントの縞数Nの高分解能マップが得られる
点である。こうすると、以下の式により、ラップされたコンポーネントをその位
相から計算する。
【0144】
【数39】 例えば、式(6)によって位相を判定する。これにより、φの範囲が±πである
場合、Nwrappedの範囲は±0.5となる。位相シフト測定ではNwra pped を高精度で判定することが可能であるが、Nが判定されないと、整数値
のオフセットが残る。
【0145】 微細な縞ほど、縞あたりの変動ρ\をより微細に表すため、微細な縞によって
、ρ\を高精度に判定することが可能である(ただし、縞がうまくアンラップ状
態になっていると仮定した場合)。例えば、微細な縞ほど、位相ステップエラー
が表面プロファイルの測定に与える影響を低減する。しかし、微細な縞ほど縞密
度も高いため、アンラップを行うのが困難である。累進縞分割によって得られる
利点は、極めて微細な縞を正確にアンラップして、可能な分解能のうち最も高い
分解能を達成することができる点である。
【0146】 累進縞分割を実施する際、ラップされたサイクルマップ(例えば、図13のロ
ウ4に示すラップされたサイクルマップ)を、複数の別個のソース間隔aそれぞ
れにおけるイメージング対象表面について生成する。これらのサイクルマップを
3つまたはそれ以下の縞サイズ(すなわち、粗大、中程度および微細)について
生成すれば十分である場合が多い。ブートストラッピングオペレーションでは、
アンラップのサイクルマップを連続的に用いて、整数オフセットにエラーを導入
することなく、次のサイクルマップをアンラップする。最も微細な縞パターンに
対応するアンラップのサイクルマップによって、最も高い分解能での測定が得ら
れる。
【0147】 図14は、累進縞分割を実行する際の例示的な一連の工程を示す。このプロシ
ージャは、(工程150)から開始し、ソース間隔を初期値aで設定する(工
程154)。この値は通常は、一連のソース間隔の中で最小であり、オブジェク
ト10上に粗大縞パターンを生成するような値が選択される。その後、式(6)
に示すような位相ステッピング(phase−stepping)測定を行うこ
とにより、ラップされたサイクルマップを生成する(工程158)。式(62)
によって位相をサイクルに変換することにより、工程158を終了する。
【0148】 工程158から独立して、aに対応する粗大縞パターンの初期のアンラップ
された推定値Nを入手する(工程162)。工程162は、工程158と同時に
行ってもよいし、工程158の前に行ってもよいし、または工程158の後に行
ってもよい。しかし、工程158および工程162はどちらとも、次の工程に続
くような様式で行わなければならない。この初期推定値を生成するアプローチに
は多くのアプローチがある。例えば、aの大きさを、オブジェクト全体が1つ
の縞内におさまる(例えば、−0.5<N≦0.5)くらいに小さくすることが
でき、これによりアンラップ作業を不要にする。あるいは、オブジェクト10に
わたって少数の縞が生成されるようにaの大きさを十分に大きくすることができ
、これにより、Nの推定値がアンラップされる全ての測定ポイントにおけるサイ
クルの半分の範囲において正確である限り、視野にわたってNが変動する際の単
純なモデル(例えば、線形傾斜(ramp))が、このNの推定値として十分と
なる。初期推定値を得るための方法の別の例は、一般的な形状または予測される
形状を仮定して、式(40)に基づいてNを計算する方法である。他の例は、先
行測定結果を同一のまたは類似するオブジェクト10に用いる方法、または、サ
イクル計数を用いて視野をサイクルに区分けする方法である。aに対応する縞
の数が小さいため、サイクル計数に必要なフレームの数も少数になる。
【0149】 次の工程(工程166)の目的は、工程158における測定から得られたラッ
プされたサイクルマップをアンラップすることである。工程162からのアンラ
ップされた推定値が1サイクルの半分以内の範囲で正確である場合、以下の数式
を用いてアンラップを達成する。
【0150】
【数40】 式(63)の右辺の2つの項のうち2番目の項は、アンラップに必要な整数オ
フセットに過ぎない点に留意されたい。この項におけるアンラップされた推定値
からラップされた測定値を減算すると、階段状の関数が得られる。この階段状の
関数の値は、整数オフセットにほぼ等しい。最も近い整数まで丸めを行うと、1
/2サイクルの制約内にエラーが収まっている限り、数式は正確になる。1/2
サイクルの制約内にエラーが収まらなくなった場合、式(63)は、Nに整数オ
フセットエラーを導入し、このエラーは、視野にわたって変動し得る。
【0151】 工程166の結果、ソース間隔aが最小である場合のNの推定値が向上する
。測定の品質をさらに上げるために、工程170において新規のソース間隔a を設定し、その結果、より微細な縞パターンがオブジェクト10上に生成される
。その後、新規の位相ステッピング測定を工程174において終了して、より微
細な縞パターンに対応する新規のラップされたサイクルマップを得る。この測定
値をアンラップするためには、aに対応する新規のアンラップされた推定値が
必要となる。工程178において、式(60)を用いて先行するアンラップされ
た推定値(これは、工程166からソース間隔aについて得られる)を適切な
レンジにスケーリングすることにより、この推定値を得る。スケーリングが行わ
れると、工程174からの新規のラップされたサイクルマップが工程182にお
いてアンラップされる。このアンラップ工程は、このサイクルマップと、工程1
78からのスケーリングされた推定値とを式(63)に代入することにより、行
われる。
【0152】 工程170〜工程182を一巡すると、Nの精度の高い推定値が得られ、これ
は、より大きなソース間隔aに対応する。ここで、最適なρ\の分解能が得ら
れるまで、徐々に微細レベルが高くなる縞(すなわち、ソース間隔が大きくなる
)を用いて工程170〜工程182を繰り返すことが可能となる(このループが
終了するたびに、ソース−間隔の下付き文字が1だけインクリメントされる)。
通常は、最適な分解能はこのような繰り返しを1〜3回行うと得られる。典型的
には、新規のソース間隔と古いソース間隔との間の比は10以上のオーダーであ
る。この比がとることのできる最大値を実際に判定する場合、Nの推定値と実際
の値との間のエラーが1/2サイクルの限定内におさまることを要求することに
より、判定を行う。所望の数の反復が完了した後、図14中の工程170〜工程
182によって表されるループを終了する。工程186において、式(19)を
通じてNをρ\に変換する。測定の目的が全ての三次元座標を判定することであ
る場合、例えば、式(47)〜(49)または式(59)を用いて三角測量を行
って、各測定ポイントについて座標を入手する(工程190)。工程194は、
累進縞分割の終了を示す。図14は累進縞分割の前に行われるプロシージャを例
示したものに過ぎず、コンポーネント工程は、複数の異なる順序によって実行す
ることが可能であることが理解されるべきである。例えば、全てのデータを収集
した後に計算を始めることが望ましい。
【0153】 累進縞分割の際に用いられる縞を広範囲のサイズにわたって生成することので
きる能力を得るためには、ソース間隔を縞間隔と同じ比率で変化させなければな
らない。ここで、他の2つの測定値から所望の間隔の縞パターンを合成する技術
について説明する。このアプローチによって得られる1つの利点として、縞パタ
ーンを合成する際に用いられる2つの測定値をより容易に入手することが可能で
ある点がある。
【0154】 縞合成は、式(61)を改変した以下の式に基づく。
【0155】
【数41】 この結果は、式(61)の両辺をラップし、そして、式(61)の右辺のN(a )−N(a)の差をラップした値は、個々にラップされたコンポーネントの
wrapped(a)−Nwrapped(a)の差をラップすることに
相当することを観察することにより、得られる。言い換えると、式(64)中の
角括弧内のラップされた数量間の差は±1の範囲にあり、この差を正しい範囲で
ある±0.5以内に収めるためには、この差をラップしなければならない。
【0156】 方程式(64)が提供するのは、ソース間隔aおよびaにおいて得られた
個々の位相ステッピング測定からの差a−aに等しいソース間隔におけるラ
ップされたNの位相ステッピング測定を合成する手段である。この方程式は、a −aの差が小さい場合に顕著な特性を有する。合成された測定値を構成する
際の個々の測定は不規則になり得、多くの不連続部分を含み得るが、視野にわた
って式(64)の左辺を平滑化して、ラップする不連続部分のうち差a−a に対応する部分の数のみを得ることができる。そのため、縞合成は、累進縞分割
を開始するための初期位相マップを生成する効果的な技術である。また、縞合成
を用いて、累進縞分割に用いられるシーケンス中の残りのアンラップの位相マッ
プのうち任意の部分を合成することも可能である。合成されたサイクルマップか
らρ\を計算する場合、式(19)を改変して以下のような式を得る。
【0157】
【数42】 ここで、Δaは、2つの測定値の間の間隔差を表す。
【0158】 図14aは、図14中の累進縞分割フローチャートに縞合成を取り入れること
を可能にする1つの方法を示す。工程154’は、工程154を改変して、間隔
に加えてベース間隔aを選択するようにした工程である。工程158’は
、工程158を改変して、間隔値aおよびaの両方におけるララップされた
サイクルマップを測定するようにした工程である。工程160’(これは、工程
158の後に挿入される)および工程176’(これは、工程174の後に挿入
される)はさらなる工程であり、式(64)を利用して、差a−aに対応す
るラップされたサイクルマップを合成する工程である。工程186’は、工程1
86を改変した工程であり、この工程186’を用いると、式(65)を用いた
合成されたサイクルマップからρ\を計算することができる。
【0159】 縞合成によって得られる1つの利点は、ソースPとソースPとの間の動き
を生成するシステムへの要求を低減できる点である。例えば、最も広範囲なソー
ス間隔を除々に少しずづ変化させると、極めて小さなソース間隔に対応する縞を
生成することができる。このアプローチは、移動要件を低減するだけではなく、
図10および図11において説明したような、ソースヘッドのいくつかの制約の
可能性も解消する。例えば、これらのソースヘッドを用いると、極めて小さなソ
ース間隔を生成するのが困難となり、ソース間隔を生成しようとすると、振動に
よって縞が不安定となり得る。縞合成によって得られる別の利点は、ソース−間
隔較正の絶対値を維持するよりも2つのソース間隔間の一定の差を繰り返す方が
通常容易であるため、測定の反復性および精度を向上させることができる点であ
る。
【0160】 縞合成のさらに別の利点は、ソース−ヘッドのアライメント不良に起因する動
きエラーを最小限にする点である。これらの動きエラーが発生し得るのは、例え
ば、図10および図11中のフォーカルスポットPおよびPがa=0のノミ
ナル位置においてオーバーラップしない場合である。考慮すべきアライメント不
良は2種類ある。第1に、フォーカルスポットPの1つが他のP2に対してd
の量だけ図面上から移動したと仮定する。これらの2つのフォーカルスポットP およびPはそれぞれ直線状の軌跡を追随し、これらの軌跡が平行であっても
、これらの2つのスポットは、その最近接アプローチ(すなわち、名目上のa=
0の位置)において距離dだけ欠損する。スポットPおよびPが最近接アプ
ローチを通過すると、オブジェクト10上の縞パターンの方向付けは、ノミナル
または所望の方向付けに対して90°だけフリップする。さらに、縞間隔は、ソ
ース間隔dに対応する縞間隔よりも大きくならない(ソースヘッドのアライメン
トは、縞がa=0近隣に適切に方向付けられた様態で保持されるまで変位dを調
節することにより、容易にされる)。
【0161】 第2の種類のアライメント不良において、フォーカルスポットPのうち1つ
が、他のフォーカルスポットPに対して動き方向99に沿って変位する。この
種類のアライメント不良が発生するのは、プリズム92が横方向に変位するため
である(すなわち、この変位は軸99に対して垂直であり、図10および11内
の面にある)。最近接アプローチにおいて、1つのソースPが、他のソースP の後側に直接配置されており、オブジェクト上に同心円状の縞パターンを生成
する。aが増加すると、これらの縞の曲率は小さくなり、自身の通常の外見に近
づく(ソース−ヘッドのアライメントをとっている間のこのエラーソースの最小
化は、小aについてプリズム92を横方向に調節して縞曲率を最小化することに
より、容易化される)。
【0162】 累進縞分割の際に動きエラーが生成する問題は、Nのスケーリングされた推定
値を生成するために式(60)が用いられた際に累進縞分割の結果得られる位相
エラーが大きくなることに起因する。動きエラーがあると、Nの推定値に半サイ
クルエラーを導入することなく用いることが可能なaの比の大きさが限定される
。縞合成を用いると、動きによって誘発される位相エラーはどちらの測定におい
ても共通し、減算において相殺されるため、この制約が解消される。そのため、
所望の形状から変形する縞をこの測定により生成し得る場合にも、合成された縞
は通常のものである。
【0163】 縞合成において2つのビーム間の位相エラーを相殺すると、これらの2つのビ
ーム間のグローバルな位相オフセットエラーを自動的に補償するというさらなる
利点が得られる。これらの位相−オフセットエラーは、例えば、2つの別個のビ
ーム間の経路−長さの差Sによって生じ、波長λの複数倍の整数ではない。位
相オフセットエラーがあると、0°以外の縞位相(すなわち、強度が最大になる
場所以外の場所)において、N=0の位置が発生する。縞合成を行うと、これら
の位相エラーは自動的に相殺され、その結果、縞パターンの位相を較正する必要
が無くなる。あるいは、縞合成によってN=0の正確な位置と判定し、そのポイ
ントにおける強めあう干渉について位相シフト要素を調節することにより、位相
較正を達成することも可能である。
【0164】 合成された縞において位相エラーを相殺することは、個々のソースPおよび
からのオブジェクト10上に降下する照射パターン中の波面エラーを相殺す
る傾向となるという別の理由においても重要である。波面エラーは、ソースの不
完全な使用または長期にわたる使用に起因し、縞パターン中に摂動またはリップ
ルを生成する。これらの摂動は、測定時に小さな高さエラーに変化する。波面エ
ラーの相殺が縞合成において発生する理由は、これらのエラーはソース間の距離
からほとんど独立しているためである。そのため、エラーは、各ソース間隔にお
いて行われる測定に共通し、式(64)における減算によって相殺される。縞合
成は、照射欠陥によるエラーを低減する有効な手段であり、ソース−ヘッド光学
部品の品質への要件を低くすることにより、システムコストを低減する機能をし
得る。
【0165】 縞合成を用いると、累進縞分割を実行することができ、ソースポイントの動く
範囲を小さくすることもでき、そのため、ソースヘッド設計をより簡単にするこ
とができる。一般的には、ソースヘッドは2種類の縞変調(すなわち、縞サイズ
の変動および位相シフトから生じる縞パターンの横方向の動き)を提供する点に
留意されたい。位相シフトのみがある場合、ソースの位置を変更することなく位
相を変更する。しかし、位相を変化させるための別の方法は、1つのソースP を他のソースPに対してオブジェクトに向かってまたはオブジェクトから離れ
て移動させる方法である。このような方法を行うとソースPのロケーションが
変化するが、ソース間の間隔が不十分である場合、この変化は極めて小さく(1
波長未満)、その結果得られるソース対PおよびPの方向付けの変化も小さ
い。縞合成の間ソース間隔は大きいままにすることができるため、ソースポイン
トPおよびPの1つまたは両方の小さな動きを通じて、どちらの種類の変調
も達成可能である。このように動き要件が少ないため、単純なソースヘッドを設
計することが可能となる。例えば、ソースヘッドは、図2bに示すファイバ−光
学部品−スプリッタ構成からなり得る。縞サイズの位相シフトおよび変動は、フ
ァイバの一方の端部を直交方向に移動させることにより、達成される。
【0166】 縞合成は、レーザ−周波数チューニングに基づいた測定に累進縞分割を使用可
能にする重要なさらなる利点を有する。縞合成が無いと、累進縞分割を測定に適
用する際、累進縞分割に必要な縞サイズの大きな変動を達成するために極めて細
密なチューニング帯域が必要となる。縞合成は、従来のチューニングレンジから
ずっと離れた縞間隔変動を生成することを可能にする。このチューニングアプロ
ーチを用いると、ソース動きがなくなり、周波数の測定制御を高い精度で行うこ
とができるため、反復性および精度が極めて高い測定値が得られる。
【0167】 周波数チューニングへの累進縞分割および縞合成の用途について、式(57)
を参照して説明する。Nは、ソース間隔aに依存するのと同様に、レーザ周波数
νに機能面において依存する(式(1)においてSによって表される経路−長
さ差はゼロであり、そうでない場合、方程式においてこのオフセットを考慮する
必要があると仮定する)。従って、Nが周波数νの関数として書かれている場合
、式(60)、(63)および(64)は、累進縞分割および縞合成の際の基本
となり、周波数チューニングにも適用される。周波数チューニングを支配する数
式は以下のようになる。
【0168】
【数43】 方程式(19)および(65)はそれぞれ、以下の改変された形態もとる。
【0169】
【数44】 ここで、Δνは、Nが2つの測定値から合成されている場合の周波数差を表す。
【0170】 図15は、縞合成および周波数チューニングを用いて累進縞分割を実行する際
に行われる工程を示す(νの代わりにaが用いられる場合、図15は縞合成およ
びソース動きを用いた累進縞分割にも適用される点に留意されたい)。この図に
おいて、プロセスの開始時に全てのデータを収集する。プロシージャが開始する
と(工程200)、ラップされたサイクルマップNwrapped(ν)、N wrapped (ν)...Nwrapped(ν)の測定(工程208)
において用いられるレーザ周波数ν、ν,...νを選択する(工程20
4)。その後、方程式(67)を用いて、周波数差ν−ν、ν−ν,.
..ν−νについて、ラップされたサイクルマップを合成する(工程212
)。その後、多くの可能なアプローチのうちの1つ(例えば、図14の工程16
2)を用いて、第1の差ν−νについてアンラップのサイクルマップの推定
値を得る(工程216)。好適なアプローチは、N(ν−ν)を±0.5−
サイクルの制限内におさまらせるくらいに十分に小さな差ν−νを用いて開
始する。工程220において、式(68)を用いて、シリーズ中の次に大きな差
ν−νについてアンラップのサイクルマップを計算する。工程224におい
て式(66)を用いてこのアンラップのサイクルマップをスケーリングし、その
後工程220にフィードバックして、シリーズ中の次の合成されたサイクルマッ
プをアンラップする。最終周波数差ν−νに対応するサイクルマップがアン
ラップされた後、ループは終了する。工程228において、式(70)を用いて
、最終アンラップのサイクルマップNを縞数ρ\に変換する。最終サイクルマッ
プNが合成されなかった場合、方程式(69)を用いる。工程232において、
例えば、式(47)〜(49)または式(59)を用いて、ρ\に対応する方向
情報を各対象ポイントPについて三次元座標に変換する。
【0171】 累進縞分割ならびに累進縞分割および縞合成の組み合わせは、計数のあいまい
さを解消し計数推定値を細密化する強力な一般的な方法であり、これらの技術は
上記にて述べた用途以外の用途にも適用可能であることが理解されるべきである
。詳細には、これらの技術は、データ取得のための上記の方法および装置または
縞解析にも限定されない。
【0172】 図16を参照して、未分解状態の(unresolved)縞またはイメージ
スペックルに対して累進縞分割法および縞合成法を用いた改変例について説明す
る。図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の差による経路長さの差により、スペックル強度
が変調し、これはレーザ周波数と共に変化する。特定のレーザ周波数における変
調の位相を、位相シフト測定を通じて正確に判定することが可能である。
【0173】 図16中のオブジェクト10の表面上に縞が生成されていない理由は、オブジ
ェクト10を照射しているのは1つのレーザビーム50だけであるためである。
2つのビーム46および50からの光のコヒーレントの重ね合わせにより、検出
器22において干渉が発生する。式(1)〜(5)によれば、レーザ周波数が変
化すると、オブジェクト10のイメージ中の個々のスペックルローブの強度が周
期的に変動し、その際の変動速度は、経路長さの差sに比例する。簡潔にするた
め、基準表面102は平面であり、x−y面にあると仮定する。zはオブジェク
ト10上のポイントPの仮想基準表面102’上の高さを表すものとし、この
zを測定する。名目上はs=2z(これは、伝播の往復による)であり、以下の
式によってNからzを計算することができる。
【0174】
【数45】 Nが縞合成を通じて得られる場合、以下のようになる。
【0175】
【数46】 他の場合 図16を用いて三次元イメージを生成するプロシージャは、図15に概要を示
すものと同じであるが、ただし、工程232において、Nをzに変換する際に式
(71)または式(72)を用いている点においてのみ異なる。図16に示すア
プローチは三角測量に依存していないため、ソースPおよびPと受信器22
および24との間のシャドーイング効果を無くすことができるという利点が得ら
れる。このアプローチは特に短距離の用途に有用であり、そのような用途として
は、経路長さの間隔が短い三次元顕微鏡法がある。
【0176】 上記の実施形態は、コヒーレント照射を用いた干渉に基づく。検出器22上に
オブジェクト10のイメージを形成する際にレンズ24を用いる実施形態の場合
、得られるイメージにスペックルがあり、スペックルによってレンジ分解能が制
限され得るという不利点が生じ得る。オブジェクトの表面上に検出器23を直接
配置すると、この制約から逃れることができる点に留意されたい。スペックルが
発生するのは、コヒーレント照射を用いてイメージングを行った場合であり、そ
の理由としては、レンズ24の有限ポイントスプレッド(finite poi
nt−spread)機能により、オブジェクト10の表面上のP周囲の領域
から光が散乱し、検出器22におけるポイントPにおいて干渉するからである
。検出器面におけるスペックルのサイズの平均は、レンズのポイントスプレッド
機能のサイズによって決まる。スペックルにより、P周囲の散乱領域による寄
与は、イメージポイントPにおいて徐々に大きくなる構成にすることができ、
これらの散乱領域により、ρ\の測定値は、一方向または反対方向に引き寄せら
れる。
【0177】 ここで、測定時のスペックルによる影響を軽減するための複数の技術を紹介す
る。実際は、検出器アレイ22の単一のピクセルに対応するエリア内におさまる
スペックルの数が多数である場合、スペックルによる効果は平均化する傾向とな
る。この条件を満たすための1つの方法として、レンズ24によって干渉縞を分
解し、レンズ24のポイントスプレッド機能の大きさを、検出器アレイ22の各
素子上におさまるスペックルの数を多くするくらいに十分に小さくすることを要
求する方法がある。このアプローチは、方位分解能を極限まで使用しなくてもよ
い状況では良好に機能する。一般的には、方位分解能とスペックルの平均化とは
トレードオフする。
【0178】 高い方位分解能が必要な状況の場合、スペックルサイズとピクセルサイズとを
より近密に整合させる必要が出てくる場合がある。一般的には、レンズのポイン
トスプレッド機能のサイドローブの大きさを小さくすることが望まれる(すなわ
ち、アポダイズする)。そうすると、ポイントスプレッド機能の強力なサイドロ
ーブ内に含まれる散乱領域の能力が低減する。このような処理を行わないと、こ
のような散乱領域の能力により、ρ\の測定に偏りが出てくる。アポダイゼーシ
ョンを行うと、方位分解能はわずかに低下し得る。場合によっては、この方位分
解能とレンジ分解能の向上とをトレードオフすると有利である場合がある。
【0179】 スペックルによる効果を低減させるさらなる技術として、個々のピクセル22
に対応するオブジェクト10の表面領域上でのρ\の変動を最小限にする技術が
ある。例えば、ピクセル密度を増加させつつピクセル22あたりのスペックルの
一定比(constant ratio)を維持すると、ピクセル22に対応す
る領域の片側から他方の側へかけてのρ\の変動が低減し、レンジ分解能および
方位分解能の両方が向上する。
【0180】 所与の光学コンフィギュレーションおよびピクセル密度において、ソースヘッ
ドの位置および方向付けを適切に選択することにより、個々のピクセル22あた
りのρ\の変動を最小化することができる。ρ\の変動が最小化されるのは、例
えば、ピクセル22に対応する表面素子が三次元縞パターンの等位相面に沿って
延びる場合である。言い換えると、スペックルによる効果が最小化されるのは、
縞パターンを照射する際に用いられる定位相表面に対して測定対象表面が平行で
あるか(またはほぼ平行である)場合である。この条件は、表面垂線に対して平
行な2つのソースポイントの方向付けを示すベクトルを備えるかすめ入射線を用
いることにより、達成可能である。特別な場合を除き、表面上のどのポイントP についてもこの条件を満たすことは不可能である。しかし、多くの実用上重要
なオブジェクト(例えば、空気力学面(例えば、車両ボディパネルおよび機体パ
ネル))については、この条件に近い条件を得ることが可能である。
【0181】 多くのオブジェクトには広範囲にわたって表面に傾斜および方向があり、その
ため、単一のソースヘッドを用いて表面全体上のρ\を最小化するための上記条
件を満たすことは不可能である。このような状況において、複数のソースヘッド
を用いることにより、この条件をより満足させることが可能である。これらの複
数のソースヘッドの配置は、複数のソースヘッドの1つがオブジェクト10の各
主要領域または最重要対象領域に対して応答性を持つように行われる。各ソース
ヘッドからの情報の出どころを区別するために特定の形態の多重化が行われると
仮定する。同じ領域を照射するソースヘッドが1つ以上ある場合、その結果得ら
れる複数の測定値に対し、当該領域における各測定において予測されるスペック
ルノイズのレベルに従って重み付けを行うことができる。複数のソースヘッドを
用いると、シャドーイングによる制約を解消する際にも有用である場合がある。
複数のソースヘッドからの情報は、オーバーラップ領域をビューイングする所与
の検出器22について容易に組み合わせることが可能である。各ピクセルに対応
する角度情報は各ソースヘッドについて同じであるため、登録は自動的に行われ
る。
【0182】 いくつかの場合においては、オブジェクト10の特定の領域を照射する工程を
ピクセル上のρ\の変動を最小化するための条件を満たすような様式で行うこと
ができない。スペックルによる効果の軽減をこのような制約による影響を受けず
に行うさらなる技術は、スペックル平均化に基づく。スペックル平均化は、異な
るパラメータにわずかな差異を設けた(例えば、光学コンフィギュレーションの
変動)さらなる測定を行うことにより、達成することが可能である。あるいは、
露光の間にパラメータを変化させることも可能である。平均化を行うと、変動す
る一連のパラメータが異なるスペックルパターンの異なる様式の統計的実現を可
能にするくらいに十分に変動した場合、スペックルによる効果は実際に低減する
【0183】 一実施形態において、ソース間の中間位置Pは、調節されたパラメータであ
る。Pの変位方向が図1および12の面から外れると、オブジェクト照射の変
動は最小になる。Pの変動については、式(47)〜(56)および(59)
もおいて明示的に説明している。スペックルによる効果の低減は、例えば、異な
る値のPについて(x、y、z)座標の個々の測定値を平均化することにより
、達成可能である。
【0184】 別の実施形態において、露光期間中、ソース間隔aおよびレーザ周波数νは同
時に変化するため、縞パターンは静止状態のままである。縞を静止状態のままに
しておくための条件は、式(57)においてa/λの比(またはaνの積)が一
定のままであるときに見られる。静止状態の縞を達成するために、フィードバッ
クメカニズム(例えば、縞モニタ)により、ポイントソースPおよびソースP の動きをレーザ40のチューニングに連動させる(またはその反対を行う)。
この縞モニタは、別個の検出器または検出器の線形アレイからなり得、ソースヘ
ッドに組み込まれて、ビームスプリッタによって分離される主要ビームから間隔
が空けられた縞パターンを測定する。
【0185】 スペックル平均化を効果的にするために、用いられるチューニング偏位を、イ
メージスペックルにおいて複数の脱相関(decorrelation)サイク
ルが発生するくらいに十分に大きくする。ソースヘッド中の別個の経路間の経路
長さのオフセットSを最小化するよう、注意が必要である。そうしないと、縞
パターンの位相は、以下の式による周波数シフトΔνと共に変動する。 Δφ=2π(ΔνS)/c あるいは、経路−長さオフセットを故意にゼロ以外の値に設定して、2πの増
分の位相Δφにおいて縞パターンを「ストロボ」して、当該縞の位相をフリーズ
することも可能である。こうすると、ストロボのタイミングを変更することによ
り、異なる位相シフトを達成することができる。
【0186】 サイクル計数、累進縞分割および累進縞分割を縞合成を用いて行うことは、計
数アンビギティーの分解および測定精度の向上を得るたるための一般的な技術で
あり、このような技術の有用性は、上記の用途をはるかに超える点が理解される
。さらに、これらの技術は、一般化された信号(例えば、一次元信号)にも適用
可能であり、シヌソイド変動に限定されない。
【0187】 これらの技術は、任意の手段によって形成された一般的な構造の光パターンに
も適用されることが理解されるべきである。構築された光パターンを生成する手
段の例を挙げると、コヒーレント放射の干渉もしくは回折、異なる波長で形成さ
れた干渉パターンもしくは回折パターンの「白色光」の重ね合わせ、およびコヒ
ーレント放射もしくは非コヒーレント放射のいずれかを用いたレンズによるオブ
ジェクト上へのパターンのイメージングがある。このようなイメージングアプロ
ーチの一例として、デジタルプロジェクタによって生成された白色光パターンの
作製を、当該パターンが拡張および接触してアコーディオン状の動き(acco
rdion motion)を生成し、側方に平行移動して位相変調を生成する
ように行うことが可能である。このような様式で白色光パターンを表面上に投射
するとスペックルによる効果を無くすことはできるが、このアプローチには特定
の不利点がある。イメージングされた光パターンには、縞周期(period)
の小ささに制約があり得る。投射されたパターンのフィールド深さも考慮する必
要がある。オブジェクト全体にわたって焦点が合わされている斜角から微細縞パ
ターンを投射するのも不可能になり得る。一方、干渉または回折によって生成さ
れた縞は、いつもオブジェクト表面に「焦点が合って」おり、この状態は、当該
縞の観測ポイントへの距離またはサイズに関係なく保たれる。
【0188】 本発明のさらなる実施形態では、干渉または回折によって縞を形成することに
よって得られる利点が維持され、広帯域の照射または白色光においてこれらの縞
を形成することにより、スペックルによる効果も無くなる。この実施形態は、照
射スペクトル中の異なるコンポーネント波長において形成される干渉縞強度の非
コヒーレントの重ね合わせに基づく。これらの縞のサイズおよび形状が同じであ
る場合、異なる波長に対応する干渉−縞強度は、パターン中の各ポイントにおい
て互いに強化し合う。
【0189】 式(57)を参照して、a/λの比を固定することにより、サイズおよび形状
が同じ縞を生成する。レンズによって追随される回折格子を用いることにより、
各コンポーネント波長についてもこの目的を同時に達成することができる。回折
角度が小さい場合、格子によって回折される1次ビームの角度偏差は、当該ビー
ムの波長に比例する。レンズは、この角度偏差をソースPまたはPの変位に
変換する。このソースPまたはPの変位もλに比例するため、a/λ比を一
定に維持することが可能となる。この実施形態において、ソースPまたはP のロケーションは、波長において広く分布する(smeared out)。白
色光の位相シフトは、格子を側方に平行移動させるかまたは例えばPancha
ratnam位相に基づいた白色光位相シフタを用いることにより、達成可能で
ある。
【0190】 これらの原理に基づいてソースヘッドをインプリメントするための光学構成に
ついては、多くの種類が可能である。1つの構成では、図10中のビーム46お
よび50中にブラッグ格子を配置する。これらの格子(図示せず)およびプリズ
ム92の配置を、非回折ビームに対応する2つのフォーカルスポットが重複し、
1次ビームがプリズム92のベースに向かって偏向するような配置にする。この
構成により、回折ビームに対応するソース間隔aは、光学波長λに比例する。白
色光位相シフトは、ブラッグ格子の1つを側方に平行移動させることにより、達
成される。縞サイズの変動は、格子の周期を変更することにより、達成される。
一実施形態において、ブラッグセルまたは音響光学変調器を用いて、格子周期を
電子的に変更する。この実施形態において、双方のブラッグセルを同じ周波数で
駆動し、1つの駆動信号の位相を他方の信号に対してシフトさせることにより、
位相シフトを達成する(各ビームは同じドップラーシフトを受けるため、これら
の縞は静止状態である)。別の実施形態において、異なる周期のブラッグ格子を
連続的に配置する。さらに別の実施形態において、単一の回折性素子をビーム4
4中に配置する。別の白色光位相シフタ(図示せず)を、ビーム46または50
のうち1つに配置する。当業者に公知の白色光位相シフタのいくつかの実施形態
において、ビームスプリッタ48は偏光ビームスプリッタである。
【0191】 ここで図17を参照して、図17は、広帯域干渉縞プロジェクタまたは白色光
干渉縞プロジェクタの一実施形態を示す。ソース250は、放射252の平行ビ
ームを実質的に生成する。この平行ビームは、回折格子254に対して実質的に
垂直入射する角度で方向付けられる。回折格子254は、格子周期Dを有する。
例示目的のため、入力ビーム252を、波長コンポーネントλ、λおよびλ から構成されるものとして表す。実際は、ビーム252は任意のスペクトルの
構成を有し得る。回折格子254は、ビーム252を複数の回折ビームに分割す
る。これらの複数の回折ビームの回折する順序は、整数mによって表すことが可
能である。例示目的のため、ビーム252周囲に沿った光線のみを図示している
。これらの回折ビームは、垂直入射に関する以下の格子方程式に従って光学軸2
58に対して角度θで伝播する。
【0192】
【数47】 一実施形態において、回折格子254の設計を、回折ビーム260の回折オー
ダm=+1であり回折ビーム262の回折オーダm=−1であるときの回折効率
を最大かつ均等にするようにな設計にする。他の実施形態において、回折格子2
54の設計を、同じオーダ|m|を有する任意の一連の正ビームおよび負ビーム
の回折効率を最大かつ均等にし、なおかつ、他の全てのオーダーに回折するエネ
ルギーを最小化するするようにな設計にするような設計にする。非回折(m=0
)ビーム264が残留した場合、そのようなビームは全て偏向しない状態で回折
格子254を通過し、レンズ266によって集束され、フォーカルスポット26
8上に来る。
【0193】 フォーカルスポット268のスペクトル成分λ、λおよびλは、実質的
に重複する。一実施形態において、フォーカルスポット268は、光学2重スリ
ット272の中央障害物270によって実質的にブロックすることが可能である
。回折ビーム260および262の異なるスペクトル成分λ、λおよびλ は、レンズ266によってスペクトル領域274および276上に集束される。
スペクトル領域274内のフォーカルスポットと、所与の波長λに対応するスペ
クトル領域276内のフォーカルスポットとの間の距離a(λ)は、波長λに実
質的に比例する。一実施形態において、レンズ266のアパチャストップ278
を用いて、望ましくない高オーダの回折ビームをブロックすることが可能である
。別の実施形態において。光学2重スリット272の不透明領域を用いて、望ま
しくないレンズ266を通過する残留回折オーダのうち望ましくないもの全てを
ブロックすることが可能である。2つのスペクトル領域274および276から
の放射は、伝播して広帯域干渉縞パターン280を形成する際、拡張および重複
する。縞パターン280は、2重スリット272からの代表的距離Rにおいて、
代表的縞周期dを有する。
【0194】 上記の式(57)などの方程式、特定のソース間隔aおよび特定の波長λに基
づき、以下の代入により、これらの方程式を、図17に示す広帯域干渉縞プロジ
ェクタに対する第1のオーダに適用することができる。
【0195】
【数48】 ここでfはレンズ266の焦点距離であり、Dは回折格子254の周期である。
【0196】 図17に示す縞生成方式を用いると、狭帯域またはレーザによる照射を用いて
縞を生成することも可能である点に留意されたい。狭帯域照射用のレンズ266
によって追随される回折格子254を用いることによって得られる利点の1つと
して、縞周期dは波長に対して感度を持たないため、ソースの周波数ずれは測定
を実質的に劣化させない。例えば、レーザダイオードは、比較的低コストであり
入手も容易なソースであるが、その動作波長は温度依存性である。しかし、この
技術は温度依存性の波長シフトに対して感度を持たないため、測定劣化を招くこ
となくレーザダイオードを用いることが可能である。
【0197】 一実施形態において、回折格子254は、肉薄の位相格子であり、方形波の位
相プロファイルを持ち、代表的波長λの場合、その相対位相遅延は、0°と1
80°との間で交互に変化し、その際のデューティサイクルは50%である。格
子254は比較的高効率であり、利用可能なエネルギーのおよそ40.5%をm
=−1およびm=+1の回折オーダそれぞれに回折し、名目上は、0%をm=0
および他の均等な回折オーダに回折する。格子254の相対的な位相遅延は波長
の関数であるため、オーダm=0のときの非回折ビーム264中のエネルギーは
、代表的波長λと異なる波長の場合、増加する。一実施形態において、格子2
54を比較的肉薄に作製すると、この格子254の波長依存レベルが低下するた
め、この格子254は、広範囲の周波数スペクトルにわたって良好に機能する。
【0198】 得られた広帯域の(または狭帯域の)干渉−縞パターン280を位相シフトす
る工程は、回折格子254を図17に示す方向282に平行移動させるだけで達
成される。白色光または広帯域の位相シフトが実現するのは、回折格子254を
格子周期Dの所与の一部分の分だけ平行移動させると、縞パターン280の各ス
ペクトル成分も、当該縞周期dの同じ一部分の分だけシフトするからである。例
えば、格子252をD/4または1/4サイクルだけ平行移動させると、干渉縞
パターン280も1/4サイクル(または90°)だけシフトする。
【0199】 干渉−縞パターン280のアコーディオン状の動き(または縞サイズの変動)
を達成できる方法には様々な種類がある。一実施形態において、回折角度θ
小さい場合、格子254の周期Dを2倍にすると、ビーム260および262の
θの大きさ(これは、式(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)を以下のように改変することが可能
である。
【0200】
【数49】 ここで、Nwrappedはf/D比の関数である。
【0201】 fおよび/またはDについて所定の別個の一連の値を設け、上述したような累
進縞分割を通じて縞数のアンビギティーを分解すると適切である場合が多い。一
実施形態において、累進縞分割は、別個の一連の回折格子254またはレンズ2
66を図17中のそれぞれの位置に連続的に移動させることにより、達成される
。別の実施形態において、格子254またはレンズ266は切換え可能であり、
別個の値のDおよびfをとる。切換え可能な格子のいくつかの例を挙げると、液
晶の空間光変調器、パターン化されたアライメントの液晶空間光変調器、電子的
に切換え可能な回折性光学素子、カスケード式の電子的に切換え可能な回折性光
学素子、および超小型電子システム(MEMS)がある。
【0202】 図18に示す一実施形態において、回折格子254は、長尺かつ平行な位相遅
延器292またはピクセルの線形アレイからなる液晶空間光変調器290であり
、その位相遅延は、周期的な位相−遅延パターン294が生成されるように制御
することが可能である。このようにして、肉薄の位相格子を生成し、電気的に作
動させることができる。例えば、ノミナルの相対位相遅延が50%のデューティ
サイクルで0°および180°に交互に変化する方形波の位相プロファイルを生
成することが可能である。その結果得られた回折格子254の周期Dおよびオフ
セット282は、別個の工程において可変性であり、高い反復性を有する。した
がって、得られた縞パターン280の縞間隔dおよび位相シフトはどちらとも、
別個の反復可能な工程において高速で電気的に制御可能である。
【0203】 図18において、wはピクセル292の幅であり、nblockは、50%
のデューティ−サイクルの肉薄の位相格子を半分生成する際に同じ位相遅延レベ
ルにおいて反復されるピクセル292の数である。一例として、図18において
block=2である。格子周期Dは、以下の式によって得られる。
【0204】
【数50】 格子周期Dは、nblockの値を変更することによって変化させることが可能
である。所与の値のnblockが以下のようになっている場合、回折格子25
4の最小位相シフトΔφおよびそれに起因する干渉−縞パターン280の最小位
相シフトΔφを生成することが可能である。 Δφ=180°/nblock したがって、例えば、nblock値が2である場合、±90°、±180°
および±270°の位相シフトを生成することができ、nblock値が3であ
る場合、±60°、±120°、±180°、±240°および±300°の位
相シフトを生成することができる。3つの任意の位相シフトφ、φおよびφ に関する一般化された形式の式(6)は、以下のように表すことが可能である
【0205】
【数51】 ここで、I、IおよびIは、位相シフトφ、φおよびφに対応する
信号強度である。そのため、位相シフト値を、位相シフト値の利用可能性にした
がって、nblockの各値に応じて個別に選択することができる。
【0206】 縞合成は、Dのさらなる値を得る手段であり、式(78)中の整数nbloc を変化させることによって達成することが可能な手段よりも効果的である。実
際の格子周期Deffの利用可能な値のうち、異なる値のnblockを用いて
2つの測定値から合成することが可能な値は、以下の式によって得られる。
【0207】
【数52】 ここで、nblock1<nblock2である。大きな値のDeffは、式(
81)を用いてnblock1およびnblock2を連続する整数にすること
により、得られる。可能な測定戦略の1つの例示として、ピクセル幅wは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であり、この値は、累進縞
分割において適切なシーケンスである。
【0208】 図17に戻って、nblockが広範囲にわたって変化すると、それに応じて
、スペクトル領域274および276の広がりおよびロケーションも広範囲に変
化する。その場合、一実施形態において、より高度なバージョンの2重スリット
272(例えば、図19に示すビームブロック296)を用いて、不要な残留回
折オーダを無くすことができる。ビームブロック296は、2つのネスト状の2
重スリットを含み、中央障害物268を含む。上記の例における測定シーケンス
において用いられる値nblock0が最小である場合、外側の2重スリット2
75はスペクトル領域274および276を通過し、一方、値nblock1
よびnblock2の値が大きい場合、内側2重スリット277は所望のスペク
トル領域274’および276’を通過する。
【0209】 図20に示すさらに別の実施形態において、レンズ300を光学軸258上に
配置する際、スペクトル領域274および276のイメージ302および304
がイメージ面306に形成されるような配置にすることにより、アコーディオン
状の動きを達成する。別の実施形態において、レンズ300の焦点距離が負であ
る場合、イメージ302および304は仮想イメージ(図示せず)であり得る点
に留意されたい。一実施形態において、レンズ300を光学軸258に沿って移
動させることによってイメージ302および304の倍率を変更することにより
、アコーディオン状の動きを得る。または、別の実施形態において、レンズ30
0の焦点距離を変更することにより、アコーディオン状の動きを得る。さらに別
の実施形態において、レンズ300はズームレンズであり、イメージ面306の
ロケーションを変化させることなく焦点距離を変化させる。さらに別の実施形態
において、焦点距離が異なる別個の一連のレンズを連続的に適切な位置に配置す
ることが可能である。レンズ300を用いると、図17内のコンポーネントのう
ち任意のコンポーネント(例えば、2重スリット272)を変化させることなく
、アコーディオン状の動きが生成される。
【0210】 広帯域干渉−縞プロジェクタのさらに別の実施形態を図21に示す。この実施
形態は、広帯域ソース250と、音響光学変調器(AOM)306と、レンズ2
66とを含む。AOM306は、乗算器310によって生成された合成信号30
8によって駆動される。この乗算器310は、可変周波数fを有するシヌソイ
ド312を、一定周波数fを有するシヌソイド314で乗算する。合成信号3
08は、周期が異なる2つの進行する音波をAOM306中に生成し、すると、
AOM306は角度が離れた2つのビーム260および262を生成する。異な
るスペクトル成分の回折ビーム260および262が、レンズ266によってス
ペクトル領域274および276に集束し、ここで、スペクトル領域274内の
フォーカルスポットと、所与の波長λに対応するスペクトル領域276内のフォ
ーカルスポットとの間の間隔a(λ)は、波長λに実質的に比例する。さらに、
各スペクトル成分について、スペクトル領域274および276のスペクトル成
分間にあるような中間点を、実質的に位置268に局所化させる。fを調節す
ることにより、各波長コンポーネントλの間隔a(λ)を、対称点268の周囲
において、ゼロの間隔からAOM306の帯域が許す最大間隔まで調節する。2
つのスペクトル領域274および276からの放射は、伝播して代表的周期dの
広帯域干渉−縞パターン280を代表的距離Rにおいて形成すると、拡張および
重複する。fを変化させると、縞パターン280は拡張および収縮する。
【0211】 AOM306にかかる圧力変動は進行するため、ビーム274および276は
、互いにドップラーシフトする。その結果、縞パターン280は、合成信号30
8の2つの周波数成分間の周波数差に比例した速度で重複領域を移動する。干渉
パターン280は、振幅変調ソース250により、特定の空間位相において2つ
のビーム260および262の周波数差2fにおいて効果的に凍結することが
可能である。縞280の空間位相は正確に制御することが可能であり、この制御
は、レーザ駆動信号316のAOM制御信号308に対する位相を変化させる(
例えば、遅延を変化させる)ことにより、行われる。あるいは、駆動信号316
を用いて、光路(図示せず)内に配置された外部の振幅変調器を制御するか、ま
たは、縞パターンを観測する際に用いられる検出器22の感度を変調することも
可能である。
【0212】 広帯域干渉−縞プロジェクタの上記の実施形態において、ビーム252は実質
的に平行であるため、回折性素子254(図17)または306(図21)は、
所与の回折オーダmの各スペクトル成分を同じ回折角度θに回折させる。この
コンフィギュレーションを他の様態に改変することも可能であり、このような改
変は、本発明の趣旨および範囲の内におさまる。例えば、別の実施形態において
、ビーム252は発散型または収束型であり得る。さらに別の実施形態において
、ビーム252は、回折性素子254に斜角で入射し得る。さらに別の実施形態
において、屈折素子および回折性素子のオーダを逆にすることもできる。
【0213】 所望の時間コヒーレンスおよび空間コヒーレンスの程度に応じて、他の多くの
光ソース250を用いてもよい。時間コヒーレンスが低い場合(広範囲のスペク
トルの内容の場合)、スペックルによる効果は低減するかまたは無くなる。一方
、空間コヒーレンスが高いと、スペクトル領域274および276中の各スペク
トル成分について、より厳密なフォーカルスポットを生成することが可能になる
。一般的には、フォーカルスポットが厳密であると、縞の可視性が向上し、縞パ
ターン280のフィールド深さが大きくなる。空間的にコヒーレントなソースの
場合、そのフィールド深さは無限である点に留意されたい。さらに、フォーカル
スポットが厳密であると、異なる回折オーダmに対応するスペクトル領域を分離
するのも容易になるため、望ましくないスペクトル領域を代表的マスク272ま
たは296によってブロックすることが可能である。
【0214】 高い空間コヒーレンスを達成するための1つの方法は、光ソース250をレー
ザベースにすることである。一実施形態において、ソース250は、平行な光学
部品を備えるレーザダイオードである。一般的には、スペクトルの幅がおよそ3
nmであるレーザダイオードを利用することが可能である。別の実施形態におい
て、異なる波長で動作するレーザダイオードのアレイを単一のビームに結合する
ことにより、このスペクトル幅を広げることが可能である。さらに別の実施形態
において、高速スキャンのためのチューナブルレーザを用いることも可能である
。レーザのスペクトルを広げるための別のアプローチでは、標準的な単一の−モ
ードシリカファイバを受動型でQスイッチ型のミクロチップレーザによってポン
ピングする。ファイバ中の様々な非線形効果により、ファイバの出力端部におい
て極めて広帯域のスペクトルが生成される。さらに別のアプローチにおいて、モ
ードロックチタン、サファイアレーザによって広帯域照射を達成する。
【0215】 ビーム252を生成するアプローチで従来のソースへの依存度がより高いもの
は他にも多くある。例えばアーク灯を光学ファイバに結合させることにより、低
い(が適切な)空間コヒーレンスを達成することが可能である。ファイバの出力
端部は、長大な(extended)ソースのように機能する。レンズを用いて
ファイバの出力を平行にすると、当該レンズの焦点距離が長くなるほど、ビーム
252の伝播角度の変動は小さくなり、フォーカルスポットも厳密になる。長大
なソースを用いると、レンズ266を光学軸258に沿った距離に配置する際、
回折格子254の表面またはAOM306の活性エリアが照射されているオブジ
ェクト10の表面または表面近隣上にイメージングされるように配置を行うこと
により、空間の非コヒーレンスによる効果が最小化される。このコンフィギュレ
ーションを用いると、縞パターン280の所与のエリアに到達した光線が、入力
ビームを異なる回折オーダに分割させる回折性素子の同一領域から発生するもの
全てと干渉するため、空間の非コヒーレンスによる効果が最小化される。
【0216】 縞パターンを照射する上記の実施形態では、位相シフトφ、φおよびφ の設定可能な精度および反復性は(狭帯域照射に基づくかまたは広帯域照射に基
づくかに関係無く)変化し得る。例えば、圧電変換器によって光コンポーネント
(例えば、図10中の鏡94もしくは96または図17中の格子254)を移動
させることにより位相シフトを行うと、ヒステリシスによって矛盾が生じる。こ
れらの位相シフトがその予想値から変化すると、式(6)または式(80)にお
いてφを計算する際にエラーが発生する。これらのφ値のエラーは周期的な性質
を持ち、再構築された表面プロファイルがリップルが生じた外観になる原因とな
る。また、リップルが生じた外観の原因としては、不完全な照射シーケンスもあ
る。これらの不完全な照射シーケンスは、例えば信号強度I、IおよびI の取得期間の間の露光レベルまたはバイアスレベルの変化によって生じ得る。露
光レベルの変換の原因としては、例えば、ソース強度の変動または完全なシャッ
ターによる露光時間差がある。
【0217】 リップルパターンは、再構築面上に発生し、物理的表面10に投射された縞パ
ターンと同じ方向に発生する。リップルパターンの主要周期は、エラーソースに
依存する。位相シフトエラーがあると、例えば、縞周期について2つのリップル
振動が発生する。リップルは、表面プロファイル測定において主要なエラーソー
スとなり得、特に、広帯域照射を通じてスペックルが無くなった場合または狭帯
域照射からのスペックル効果が上述した軽減技術の1つを通じて抑制された場合
、その傾向は強くなる。
【0218】 ここで、リップルエラーを修正するためのヒストグラムベースのリップル消去
(derippling)と呼ばれる汎用性がありかつロバストな技術を紹介す
る。この技術において、データを解析して、実際の位相オフセットおよび露光レ
ベルの変動を判定し、相関値を用いてφを計算する。この技術を説明するため、
まず、式(80)を書き換えて、推測値および推測値から計算される数量を実際
の値と区別するような様式にする。
【0219】
【数53】 式(82)において、*φ(本明細書中、以下、*φは
【0220】
【数54】 を表すものとする)、*φおよび*φを位相シフトの推定値と仮定し、*φ
を、これらの推測値を用いて計算した位相値とする。
【0221】 強度を以下のように書くことにより、不完全な照射シーケンスによって発生す
る変化を式(82)に取り入れることができる。
【0222】
【数55】 ここで、下付き文字iは1、2または3の値をとり、照射シーケンスの特定の要
素を表す。パラメータAおよびBにより、照射シーケンスの露光レベルおよ
びバイアスレベルをそれぞれ変化させることができる。式(83)において、ノ
ミナルの山−谷振幅の幅(swing)が評価対象ユニットであると一般性を失
うことなく仮定することができる。
【0223】 AおよびBの値が既知であるかまたは*Aおよび*Bとして推定され
る場合、不完全な照射シーケンスを、計算された強度により修正して、
【0224】
【数56】 式(82)中の強度測定値Iの代わりに用いることができる。さらに、推測値
*φ、*φおよび*φではなく実際の値φ、φおよびφを判定して
式(82)において用いることを可能にすると、位相ステップエラーを無くすこ
とが可能となる。リップルを最小化するための1つのアプローチとして、*φに
おけるリップルの程度を示す測定基準(metric)が最小化されるようにパ
ラメータA、Bおよびφを変化させるアプローチがある。測定基準を最小
化する場合、演算要求が多大になり、多次元のサーチアルゴリズムが必要となる
だけではなく、代表的な表面部分にわたって位相関数*φを最初にアンラップし
なければ、最小化のための適切な測定基準を得るのは困難である。アンラップが
行われた後に利用することが可能な測定基準の一例としては、*φと平滑化され
た*φの値との間の差の標準偏差がある。リップルの低減による寄与のため、こ
の測定基準は低減する傾向となる。
【0225】 ヒストグラムベースのリップル消去は、位相アンラップが不要であり、パラメ
ータの修正値を直接生成するため、上記のリップル消去アプローチよりも優れて
いる。ここで、図22を参照して、この技術を引き続き紹介する。図22は、エ
ラー(φ−*φ)/(2π)のプロットを、異なる値のパラメータA、B
よびφについて計算された数量*φ/(2π)の関数として含む。このエラー
は、縞数エラーΔNのリップルによって生じる(ripple−induced
)成分に相当し、式(32)に示すようなレンジ分解能を限定する。これらの曲
線において、推定される位相ステップ値は、*φ=−120°、*φ=0°
および*φ=120°である。φおよび*φの両方を、−π〜+πのラップさ
れた数量であるとみなし、これにより、φ/(2π)および*φ/(2π)レン
ジが−1/2〜+1/2となるようにする。曲線(a)は、振幅変化A=0.
9、Α=1、A=1.2の効果を示し、曲線(b)は、バイアス変化B
0.2、B=0、B3=0.1に相当し、曲線(c)は、推測値と異なる実際
の位相シフトφ=−105°、φ=0°およびφ=100°による効果を
示す。最後に、図22aは、図22中の曲線の傾斜をプロットしたものであり、
これは、*φに関するφの導関数およびによって得られ、φ’(*φ)として示
される。この傾斜は、1の周囲で変動し、推定パラメータがその真なる値に近づ
くと、全ての点において1に収束する。ΔNにおけるエラーが比較的小さい場合
、比較的大きな傾斜の変化が生成されるため、傾斜により、A、Bおよびφ の推測値中のエラーが感度良く示される点に留意されたい。
【0226】 図22および図22aは、A、Bおよびφに関する情報を多量に含むこ
とは明らかである。例えば、これらのプロットは、バイアスが変化すると、縞あ
たり1つのリップルが発生し、位相ステップエラーが発生すると、縞あたり2つ
のリップルが発生することを示す。振幅が変化すると、縞あたり1つおよび2つ
のリップルの混合物が発生する。図22または図22aにおいてグラフ化されて
いる曲線をデータ解析から生成することが可能であれば、A、Bおよびφ 上のこれらの曲線を支配する方程式への依存を用いて、これらのパラメータに関
する情報を回復させることができる。その後、これらの修正されたパラメータを
用いてφを計算することにより、データをリップル消去することができる。この
アプローチの難点は、これらのグラフでは、φの真の値が分かっていることを前
提としている点である。
【0227】 ヒストグラムによるアプローチを用いれば、この難点を解消することができる
。このアプローチでは、先ず、測定(または特定の測定部分)における*φの各
値の発生を相対的に表すヒストグラムを作成する。このヒストグラムは、*φの
値の可能な範囲を有限数のビンに分割し、各ビン中に収められた*φが発生する
回数をカウントすることにより、生成される。*φのヒストグラムを作成すると
、*φの確率密度関数の近似値が得られ、これはP(*φ)と書かれる。以下の
周知の関係を通じて、2つの確率密度関数P(*φ)およびP(φ)を関連付け
ることが可能である。
【0228】
【数57】 同式において、φは*φの一価の関数であると仮定する。式(85)中の導関数
を、図22a中にグラフとして示す。確率分布関数P(*φ)は、*φのヒスト
グラムの計算を通じたデータから利用することが可能である。式(85)中のP
(φ)は不明であるが、多くの場合において、φはその可能な値の範囲にわたっ
て均等に分布していると仮定すると妥当である。サンプリング対象となる表面エ
リアを覆う縞の数が多いほど、この仮定の精度も向上する。実際にも、P(φ)
が均等に分布しているというこの仮定は、縞が比較的少数である場合にも適切で
あることが分かっている。しかし、P(φ)の関数の形態が表面に関するさらな
る情報に基づいて正確になっている場合、ヒストグラムベースのリップル消去に
おいて他の形態のP(φ)を用いることも可能である。
【0229】 P(Φ)が分かると、式(85)を用いて、数量P(*φ)(これは、測定値
*φのヒストグラムを計算することによって得られる)と、導関数φ’(*φ)
の絶対値とを関連付けることができる。その後、この導関数の関数の形態に基づ
いて、A、Bおよびφの実際の値を判定することができる。これを行う方
法を例示するため、ここで、φ’(*φ)の関数の形態をより詳細に解析する。
【0230】 先ず、振幅またはバイアスは変動しない(A=1およびB=0)が、位相
ステップエラーのみがある場合を想定する。この場合、φ’(*φ)は、以下の
ような関数の形態をとることが分かる。
【0231】
【数58】 ここで、K、Kc2およびKs2は定数であり、その値は位相ステップ値*φ および実際の位相ステップ値φに依存する。Kの値(これは、1に近い数であ
る)は他の2つの定数KC2およびKS2から判定することが可能であり、その
結果、式(86)における自由度は2だけとなる。これらの定数KC2およびK s2 は、例えば、このヒストグラムの逆数をとり、フーリエ級数展開のcos(
2*φ)項およびsin(2*φ)項に対応するフーリエ係数を計算することに
より、*φのヒストグラムの形状から得ることができる。その後、フーリエ級数
のDC項を用いてcos(2*φ)項およびsin(2*φ)項のフーリエ係数
を正規化して、Kc2およびKs2をそれぞれ得る。
【0232】 式(86)の自由度は2しかないため、この方程式から3つの位相シフトφ 、φおよびφの値を完全に回復することは不可能である。その代わり、相対
的位相シフトφ−φ、φ−φおよびφ−φを回復することが可能で
あり、これらの相対的位相シフトのうち2つのにより、第3の値を判定する。こ
れらの相対的位相シフトがわかっていれば、リップルを完全に除去するには十分
である。しかし、得られた*φの向上した計算値には、わずかなオフセットエラ
ーが存在する。このオフセットエラーは、φの真の値が分からないことに起因
する。
【0233】 簡潔にするため、推定された位相シフトが等しく、その大きさはφstep
あり、ゼロの周囲に集中していると仮定する。そうすると、これらの位相シフト
を*φ=−φstep、*φ=0、および*φ=φstepと表すことが
できる。デフォルトでφ=0であるとも仮定する。次いで、以下の関係により
、Kc2およびKs2からφおよびφの値を判定することができる。
【0234】
【数59】 上記の式(6)における仮定の場合(φstep=90°)、式(87)および
(88)は以下のように整理される。
【0235】
【数60】 ヒストグラムベースのリップル消去の用途のさらなる例として、位相ステップ
エラーおよび振幅エラーの組み合わせにこの技術を適用する。この場合の数式φ
’(*φ)は式(86)よりもずっと複雑であるが、この数式は、A−1およ
びφ−*φの値が小さい場合、以下のように近似化することができる。
【0236】
【数61】 ここでも、係数Kc1、Ks1、Kc2およびKs2をフーリエ解析を通じて
得ることが可能であり、これらの係数は、φ’(*φ)を近似化する際に用いら
れるヒストグラムの逆数のフーリエ級数展開の対応するフーリエ係数である。
【0237】 *φ=φ=0およびΑ=1と設定し、*φおよび*φの推測値のス
テップサイズを不均等にすると、以下の関係を通じて、Kc1、Ks1、Kc2 、およびKs2からφ、φ、AおよびAの新規近似値を判定することが
できる。
【0238】
【数62】 式(92)〜(95)は、式(91)によって得られた微小エラーの近似に基
づくため、φ、φ、AおよびAの計算値は向上するが、正確ではない。
所望であれば、このプロシージャを繰り返すことにより、精度を高くすることが
可能である。そうするためには、式(92)および(93)から計算されたφ およびφの値は、式(82)において新規の推測値*φおよび*φとなり
、式(94)および(95)から得られたAおよびAの値は、式(84)に
おいて新規の推測値*Aおよび*Aとなる。このプロシージャは高速で収束
するため、数回の反復だけで小数位のレベルでの高い精度が得られる。
【0239】 さらなる利点として、この反復によるアプローチの文脈において、バイアス変
化Bを自動処理する。反復回数が多いほど、バイアス変化Bは、他のパラメ
ータAおよびφに折り畳まれて、その効果は無くなっていく。このアプロー
チを用いると、リップルが無くなり、φが正確には分かっていない場合に生じ
るエラーに類似する微小オフセットエラーが発生する。
【0240】 本発明の好適な実施形態について説明および図示してきたが、本発明のコンセ
プトを取り入れた他の実施形態を用いることも可能であり、また、本発明の範囲
および趣旨内において多くの改変例が可能であることが当業者にとって明らかで
ある。したがって、これらの実施形態は、開示された実施形態に限定されるので
はなく、本明細書内の特許請求の範囲の趣旨および範囲のみによって限定される
べきであることが認識される。
【図面の簡単な説明】
【図1】 図1は、表面輪郭測定を行う本発明の一実施形態のブロック図である。
【図2】 図2は、図1に示す2つの放射ソースを生成するシステムの一実施形態のブロ
ック図である。
【図2a】 図2aは、図1に示す2つの放射ソースを生成するシステムの別の実施形態の
ブロック図である。
【図2b】 図2bは、図1に示す2つの放射ソースを生成するシステムのさらに別の実施
形態のブロック図である。
【図3】 図3は、図1の2つの放射ソースを互いに対して固定した距離に支持する装置
の一実施形態のブロック図である。
【図4】 図4は、図1のイメージングシステムの別の実施形態である。
【図5】 図5は、表面輪郭測定を行う本発明の別の実施形態のブロック図である。
【図6】 図6は、表面輪郭測定を行う際に図1および図5のプロセッサによって用いら
れる工程の一実施形態のフローチャートである。
【図6a】 図6aは、図6のフローチャートの一部の一実施形態である。
【図6b】 図6bは、図6のフローチャートの一部の別の実施形態である。
【図6c】 図6cは、図6のフローチャートの一部のさらに別の実施形態である。
【図7】 図7は、図1および図5のシステムによって用いられる検出器およびプロセッ
サ装置の一実施形態のブロック図である。
【図7a】 図7aは、図1および図5のシステムによって用いられるマルチプロセッサを
含む検出器およびプロセッサ装置の別の実施形態のブロック図である。
【図7b】 図7bは、図1および図5のシステムにおいて用いられる検出器およびプロセ
ッサ装置の別の実施形態のブロック図である。
【図8】 図8は、表面輪郭測定を行う本発明の別の実施形態のブロック図である。
【図9】 図9は、図8の実施形態による、オブジェクトの表面上に配置された光検出器
の素子の一実施形態を示す。
【図9a】 図9aは、図8の実施形態に用いられるスプリングアーム上に配置された光検
出器の素子の別の実施形態を示す。
【図10】 図10は、図1に示す2つの放射ソースを生成するシステムのさらに別の実施
形態のブロック図である。
【図11】 図11は、図1に示す2つの放射ソースを生成するシステムのさらに別の実施
形態のブロック図である。
【図12】 図12は、図1の実施形態のブロック図である。
【図12A】 図12Aは、図1の実施形態のブロック図である。
【図13】 図13は、累進縞分割を実行するために用いられる一連の工程の一実施形態を
示す。
【図14】 図14は、累進縞分割を実行するために用いられる一連の工程の一実施形態の
フローチャートである。
【図14a】 図14aは、縞合成を用いる累進縞分割を実行するために用いられる工程の一
実施形態のフローチャートである。
【図15】 図15は、縞合成および周波数チューニングを用いた累進縞分割の方法の一実
施形態のフローチャートである。
【図16】 図16は、図15の方法で用いられる本発明の1実施形態のブロック図を示す
【図17】 図17は、本発明の広域幅干渉縞プロジェクタの一実施形態を示す。
【図18】 図18は、本発明の一実施形態で用いられる液晶空間光モジュレータの一実施
形態を示す。
【図19】 図19は、本発明の一実施形態で用いられる2つのネスト状の2重スリットの
一実施形態を示す。
【図20】 図20は、アコーディオン状の縞運動を生成する本発明の一実施形態を示す。
【図21】 図21は、本発明の広域幅干渉縞プロジェクタの一実施形態を示す。
【図22】 図22は、本発明の一実施形態の、ヒストグラムベースのリップル消去として
公知の*φ/(2π)の関数として縞数誤差ΔNをグラフィカルに示す表示を含
む。
【図22a】 図22aは、図22に示す曲線の傾斜+1のグラフである。
───────────────────────────────────────────────────── フロントページの続き (81)指定国 EP(AT,BE,CH,CY, DE,DK,ES,FI,FR,GB,GR,IE,I T,LU,MC,NL,PT,SE,TR),OA(BF ,BJ,CF,CG,CI,CM,GA,GN,GW, ML,MR,NE,SN,TD,TG),AP(GH,G M,KE,LS,MW,MZ,SD,SL,SZ,TZ ,UG,ZW),EA(AM,AZ,BY,KG,KZ, MD,RU,TJ,TM),AE,AG,AL,AM, AT,AU,AZ,BA,BB,BG,BR,BY,B Z,CA,CH,CN,CR,CU,CZ,DE,DK ,DM,DZ,EE,ES,FI,GB,GD,GE, GH,GM,HR,HU,ID,IL,IN,IS,J P,KE,KG,KP,KR,KZ,LC,LK,LR ,LS,LT,LU,LV,MA,MD,MG,MK, MN,MW,MX,MZ,NO,NZ,PL,PT,R O,RU,SD,SE,SG,SI,SK,SL,TJ ,TM,TR,TT,TZ,UA,UG,UZ,VN, YU,ZA,ZW (72)発明者 シャーリー, ライル アメリカ合衆国 マサチューセッツ 01719, ボックスボロ, フラッグ ヒ ル ロード 105 Fターム(参考) 2F065 AA04 AA53 BB05 DD03 FF07 GG04 HH06 HH12 HH18 JJ03 JJ26 LL02 LL04 LL13 LL32 LL62 MM16 QQ31 QQ41 UU05 UU07

Claims (120)

    【特許請求の範囲】
  1. 【請求項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】 前記2つの放射ソースは互いにコヒーレントである、請求項
    1に記載の方法。
  3. 【請求項3】 前記2つの放射ソースのスペクトル分布は狭帯域である、請
    求項1に記載の方法。
  4. 【請求項4】 前記2つの放射ソースを提供する工程は、単一の放射ソース
    からの放射を分割する工程を包含する、請求項1に記載の方法。
  5. 【請求項5】 前記2つの放射ソースを提供する工程は、2つのレーザ放射
    ソースを提供する工程を包含する、請求項1に記載の方法。
  6. 【請求項6】 前記工程第1の縞パターンおよび第2の縞パターン中の縞数
    を推定する工程は、数学的モデルを適用する工程を包含する、請求項1に記載の
    方法。
  7. 【請求項7】 前記第1の縞パターンおよび第2の縞パターン中の縞数を推
    定する工程は、線形の推定モデルを適用する工程を包含する、請求項1に記載の
    方法。
  8. 【請求項8】 前記第1の縞パターンおよび第2の縞パターン中の縞数を推
    定する工程は、異なる表面について判定された三次元位置情報を適用する工程を
    包含する、請求項1に記載の方法。
  9. 【請求項9】 前記ソースの1つのスペクトル領域はそれぞれ、該ソースの
    残りの各スペクトル領域から、該スペクトル領域の各波長に比例する各距離の分
    だけ隔離される、請求項1に記載の方法。
  10. 【請求項10】 前記距離は、前記スペクトル領域の波長に線形比例する、
    請求項9に記載の方法。
  11. 【請求項11】 前記第1の縞パターンを変更する工程は、前記スペクトル
    領域それぞれを、該スペクトル領域の残りに対してある距離の分だけ隔離する工
    程を包含する、請求項1に記載の方法。
  12. 【請求項12】 前記距離は、2つのスペクトル領域の各々から等距離であ
    る中間点を含み、該スペクトルの領域それぞれを隔離する工程は、該距離を変更
    して該中間点を固定状態にする工程を包含する、請求項11に記載の方法。
  13. 【請求項13】 前記第2の縞パターンを移動させる工程は、前記ソースの
    1つからの前記空間分布中のスペクトル領域の位相を、該ソースの残りからの空
    間分布中の各スペクトル領域の位相に対して変化させる工程を包含する、請求項
    1に記載の方法。
  14. 【請求項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. 【請求項15】 前記2つの放射ソースは互いにコヒーレントである、請求
    項14に記載の方法。
  16. 【請求項16】 前記2つの放射ソースのスペクトル分布は狭帯域である、
    請求項14に記載の方法。
  17. 【請求項17】 前記2つの放射ソースを提供する工程は、単一の放射ソー
    スからの放射を分割する工程を包含する、請求項14に記載の方法。
  18. 【請求項18】 前記2つの放射ソースを提供する工程は、2つのレーザ放
    射ソースを提供する工程を包含する、請求項14に記載の方法。
  19. 【請求項19】 前記ソースの1つのスペクトル領域はそれぞれ、該ソース
    の残りの各スペクトル領域から、該スペクトル領域の各波長に比例する距離の分
    だけ隔離される、請求項14に記載の方法。
  20. 【請求項20】 前記距離は、前記スペクトル領域の波長に線形比例する、
    請求項19に記載の方法。
  21. 【請求項21】 前記第1の縞パターンを変更する工程は、前記スペクトル
    領域それぞれを、該スペクトル領域の残りに対してある距離の分だけ隔離する工
    程を包含する、請求項14に記載の方法。
  22. 【請求項22】 前記距離は、2つのスペクトル領域の各々から等距離であ
    る中間点を含み、該スペクトルの領域それぞれを隔離する工程は、該距離を変更
    して該中間点を固定状態にする工程を包含する、請求項21に記載の方法。
  23. 【請求項23】 前記第2の縞パターンを移動させる工程は、前記ソースの
    1つからの前記空間分布中のスペクトル領域の位相を、該ソースの残りからの空
    間分布中の各スペクトル領域の位相に対して変化させる工程を包含する、請求項
    14に記載の方法。
  24. 【請求項24】 前記2つの放射ソースは第1の周波数を有し、前記第1の
    縞パターンを変更して第2の縞パターンを生成する工程は、該第1の周波数を第
    2の周波数に変更する工程を包含する、請求項14に記載の方法。
  25. 【請求項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. 【請求項26】 前記第1の干渉パターンを変更して第2の干渉パターンを
    生成する工程は、前記第1の周波数を第2の周波数に変更する工程を包含する、
    請求項25に記載の方法。
  27. 【請求項27】 前記表面をイメージ面上にイメージングする工程をさらに
    包含する、請求項25に記載の方法。
  28. 【請求項28】 前記イメージ面は検出器を含む、請求項27に記載の方法
  29. 【請求項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. 【請求項30】 前記第1の放射ビームおよび第2の放射ビームは第1の周
    波数を有し、前記第1の干渉パターンを変更して第2の干渉パターンを生成する
    工程は、該第1の周波数を第2の周波数に変更する工程を包含する、請求項29
    に記載の方法。
  31. 【請求項31】 前記ラップされたサイクルマップを用いることによって三
    次元位置情報を判定する工程をさらに包含する、請求項29に記載の方法。
  32. 【請求項32】 前記表面をイメージ面上にイメージングする工程をさらに
    包含する、請求項29に記載の方法。
  33. 【請求項33】 前記イメージ面は検出器を含む、請求項32に記載の方法
  34. 【請求項34】 縞をオブジェクト表面上に投射する装置であって、 a)スペクトル分布を有する2つの放射ソースと、 b)該2つのソースと光学的に連絡するコリメータであって、2つの広帯域放
    射の実質的に平行なビームを生成する、コリメータと、 c)該コリメータと光学的に連絡する回折性格子と、 d)該回折性格子と光学的に連絡するレンズであって、スペクトル領域の空間
    分布を有する2つの放射イメージを生成するレンズと、 を備える、装置。
  35. 【請求項35】 前記ソースの1つのスペクトル領域はそれぞれ、該ソース
    の残りの各スペクトル領域から、該スペクトル領域の各波長に比例する各距離の
    分だけ隔離される、請求項34に記載の装置。
  36. 【請求項36】 前記距離は、前記スペクトル領域の波長に線形比例する、
    請求項35に記載の装置。
  37. 【請求項37】 前記距離は、2つのスペクトル領域の各々から等距離であ
    る中間点を含み、該中間点は固定される、請求項35に記載の装置。
  38. 【請求項38】 前記2つの放射ソースは互いにコヒーレントである、請求
    項34に記載の装置。
  39. 【請求項39】 前記2つの放射ソースのスペクトル分布は狭帯域である、
    請求項34に記載の装置。
  40. 【請求項40】 前記オブジェクトの表面上のポイントの三次元位置情報を
    判定する検出器をさらに備える、請求項34に記載の装置。
  41. 【請求項41】 前記2つの放射ソースは単一の放射ソースから生成される
    、請求項34に記載の装置。
  42. 【請求項42】 前記回折性格子に結合されたトランスレータをさらに備え
    、該トランスレータは、前記スペクトル領域の1つの該スペクトル領域の残りに
    対する相対的位相をシフトさせる、請求項34に記載の装置。
  43. 【請求項43】 オーダをブロックするマスクをさらに備える、請求項34
    に記載の縞をオブジェクト表面上に投射する装置。
  44. 【請求項44】 オーダをブロックする2重スリットのマスクをさらに備え
    る、請求項34に記載の縞をオブジェクト表面上に投射する装置。
  45. 【請求項45】 縞をオブジェクト表面上に投射する方法であって、 a)ある距離だけ隔離された2つの放射ソースを提供する工程であって、該ソ
    ースはそれぞれ、スペクトル分布を有する放射を生成する、工程と、 b)該放射を平行にして、2つの実質的に平行な放射ビームを生成する工程と、 c)該平行な放射ビームのスペクトル成分を描写する工程と、 d)スペクトル成分の空間分布を有する2つのイメージ放射を生成する工程と、 を包含する、方法。
  46. 【請求項46】 前記空間分布の1つのスペクトル成分はそれぞれ、該空間
    分布の残りの各スペクトル成分から、該スペクトル成分の各波長に比例する距離
    だけ隔離される、 請求項45に記載の方法。
  47. 【請求項47】 前記距離は、前記スペクトル成分の波長に線形比例する、
    請求項46に記載の方法。
  48. 【請求項48】 前記距離は、2つのスペクトル成分の各々から等距離であ
    る中間点を含み、該中間点は固定される、請求項46に記載の方法。
  49. 【請求項49】 前記2つの放射ソースは互いにコヒーレントである、請求
    項45に記載の方法。
  50. 【請求項50】 前記2つの放射ソースのスペクトル分布は狭帯域である、
    請求項45に記載の方法。
  51. 【請求項51】 前記オブジェクトの表面上のポイントの三次元位置情報を
    判定する工程をさらに包含する、請求項45に記載の方法。
  52. 【請求項52】 前記2つの放射ソースを提供する工程は、単一の放射ソー
    スからの放射を分割する工程を包含する、請求項45に記載の方法。
  53. 【請求項53】 前記スペクトル領域の1つの該スペクトル領域の残りに
    対する相対的位相をシフトさせる工程をさらに包含する、請求項45に記載の方
    法。
  54. 【請求項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. 【請求項55】 前記第1の投射パターンを変更する工程は、第2の周期的
    パターンを投射して、第2の投射パターンを第1の位置において生成する工程を
    包含する、請求項54に記載の方法。
  56. 【請求項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. 【請求項57】 前記第1の投射パターンを変更する工程は、第2の周期的
    パターンを投射して、第2の投射パターンを第1の位置において生成する工程を
    包含する、請求項56に記載の方法。
  58. 【請求項58】 スペックルがオブジェクト表面上の測定ポイントに与える
    影響を軽減する方法であって、 a)コヒーレントな縞パターンを生成する工程と、 b)該コヒーレントな縞パターンが該オブジェクト表面を実質的にかすめるよ
    うに、該縞パターンを光路に沿って該オブジェクト表面に投射する工程と、 c)該オブジェクト表面のイメージ中の該縞パターンおよび該スペックルを検
    出する工程であって、該オブジェクト表面に対する法線は該光学経路に実質的に
    直交する、工程と、 を包含する、方法。
  59. 【請求項59】 前記コヒーレントな縞パターンは、前記オブジェクト表面
    に対して0〜45°の角度で該オブジェクト表面を実質的にかすめる、請求項5
    8に記載の方法。
  60. 【請求項60】 前記2つの放射ソースを提供する工程をさらに包含する、
    請求項58に記載の方法。
  61. 【請求項61】 前記2つのソースは、前記オブジェクト表面の上側に配置
    される、請求項60に記載の方法。
  62. 【請求項62】 前記2つのソースは、前記オブジェクト表面に対する法線
    に沿って実質的に垂直である、請求項60に記載の方法。
  63. 【請求項63】 前記縞パターンはソースによって生成される、請求項58
    に記載の方法。
  64. 【請求項64】 前記2つのソースは互いにコヒーレントである、請求項6
    3に記載の方法。
  65. 【請求項65】 前記2つのソースはレーザソースである、請求項63に記
    載の方法。
  66. 【請求項66】 前記2つのソースは、単一のソースを分割することによっ
    て生成される、請求項63に記載の方法。
  67. 【請求項67】 スペックルがオブジェクト表面上の測定ポイントに与える
    影響を軽減する方法であって、 a)ある距離だけ隔離された2つの放射ソースからコヒーレントな縞パターン
    を生成する工程と、 b)該コヒーレントな縞パターンを該オブジェクト表面上に投射する工程と、 c)該オブジェクト表面のイメージ内において該縞パターンおよび該スペック
    ルを検出する工程と、 d)該縞パターンが実質的に静止状態でありかつ該スペックルが変化するよう
    に、該2つのソースを平行移動させる工程と、 e)該オブジェクト表面のイメージにおいて、該新規の縞パターンと、該変化
    したスペックルとを検出する工程と、 f)該検出された縞パターンおよび該スペックルの変化に応答して、実質的に
    スペックルの発生が無い状態で該縞パターンを判定する工程と、 を包含する、方法。
  68. 【請求項68】 工程d)および工程e)を繰り返し反復する工程をさらに
    包含する、請求項67に記載の方法。
  69. 【請求項69】 前記2つの放射ソースはレーザ放射である、請求項67に
    記載の方法。
  70. 【請求項70】 スペックルがオブジェクト表面上の測定ポイントに与える
    影響を軽減する方法であって、 a)ある距離だけ隔離された2つの放射ソースからコヒーレントな縞パターン
    を生成する工程と、 b)該コヒーレントな縞パターンを該オブジェクト表面上に投射する工程と、 c)該オブジェクト表面のイメージ内において該縞パターンおよび該スペック
    ルを検出する工程と、 d)該オブジェクトを該2つの放射ソースに対して側方に該縞パターンの等位
    相面に平行な経路に沿って平行移動させて、該縞パターンが該オブジェクト表面
    に対して実質的に静止状態でありかつ該スペックルが変化するようにする工程と
    、 e)該オブジェクト表面のイメージにおいて、該新規の縞パターンと、該変化
    したスペックルとを検出する工程と、 f)該検出された縞パターンおよび該スペックルの変化に応答して、実質的に
    スペックルの発生が無い状態で該縞パターンを判定する工程と、 を包含する、方法。
  71. 【請求項71】 工程d)および工程e)を繰り返し反復する工程をさらに
    包含する、請求項70に記載の方法。
  72. 【請求項72】 前記2つの放射ソースはレーザ放射である、請求項70に
    記載の方法。
  73. 【請求項73】 スペックルがオブジェクト表面上の測定ポイントに与える
    影響を軽減する方法であって、 a)該オブジェクト表面上にコヒーレントな縞パターンを投射する工程と、 b)透過関数を有するレンズを提供する工程であって、該レンズの透過は該レ
    ンズの端部において徐々に低下する、工程と、 c)該オブジェクト表面のイメージ内に該縞パターンを検出する工程であって
    、該透過関数は、該スペックルが該測定に与える影響を実質的に低減する、工程
    と、 を包含する、方法。
  74. 【請求項74】 スペックルがオブジェクト表面上の測定ポイントに与える
    影響を軽減する方法であって、 a)第1の周波数を持ちかつある距離だけ隔離された2つの放射ソースから第1
    の縞パターンを生成する工程と、 b)該第1の周波数を第2の周波数に変更することにより、該第1の縞パターン
    を変更して第2の縞パターンを生成する工程と、 c)該第1の周波数と該第2の周波数との間の差に応答して該距離を変更する工
    程であって、該第1の周波数と該第2の周波数との間の差に対する該距離の比は
    実質的に一定である、工程と、 を包含する、方法。
  75. 【請求項75】 縞をオブジェクト表面上に投射する方法であって、 a)ある距離だけ隔離された2つの放射ソースを提供する工程であって、該ソ
    ースはそれぞれ、スペクトル分布を有し、他方のソースに対してコヒーレントで
    ある、工程と、 b)該オブジェクト表面上のポイントを、該ソースそれぞれからの放射によっ
    て照射する工程と、 c)該ソースの1つを他方のソースに対して移動させる工程と、 d)該オブジェクト表面上のポイントによって散乱された放射を検出する工程
    と、 を包含する、方法。
  76. 【請求項76】 前記ソースの移動と、前記オブジェクト表面上のポイント
    によって散乱された放射の検出とに応答して位置情報を計算する工程をさらに包
    含する、請求項75に記載の方法。
  77. 【請求項77】 前記ソースの1つからのスペクトル分布中のスペクトル成
    分の位相を、該ソースのうちの他方のソースからのスペクトル分布中の各スペク
    トル成分の位相を前記オブジェクト表面上のポイントにおいて測定した値に対し
    て変化させる工程をさらに包含する、請求項75に記載の方法。
  78. 【請求項78】 前記ソースの1つを該ソースのうちの他方のソースに対し
    て移動させる工程は、該ソースのうちの1つのスペクトル成分と該ソースのうち
    の他方のソースの各スペクトル成分と間の距離を変化させて、該ソースの各スペ
    クトル成分間の距離を変化させる工程を包含する、請求項75に記載の方法。
  79. 【請求項79】 前記距離は、前記ソースのスペクトル成分の各々から等距
    離である中間点を含み、該ソースのうちの1つを該ソースのうちの他方のソース
    に対して移動させる工程は、該ソースのうちの1つのスペクトル成分と、該ソー
    スのうちの他方のソースの各スペクトル成分と間の距離を変化させて、該ソース
    の各スペクトル成分間の距離を変化させ、該中間点を固定状態に保持する工程を
    包含する、請求項78に記載の方法。
  80. 【請求項80】 前記オブジェクト表面上のポイントの三次元位置情報を判
    定する工程をさらに包含する、請求項75に記載の方法。
  81. 【請求項81】 前記2つの放射ソースを提供する工程は、 a)あるスペクトル値を有する初期放射ビームを提供する工程と、 b)該初期ビーム放射から、第1の放射ビームを第1のビーム角度においてかつ
    第2の放射ビームを第2のビーム角度において生成する工程と、 c)該第1の放射ビームおよび該第2の放射ビームをイメージングして、該2つ
    の放射ソースを形成する工程と、 を包含する、請求項75に記載の方法。
  82. 【請求項82】 前記初期ビーム放射振幅変調する工程をさらに包含する、
    請求項81に記載の方法。
  83. 【請求項83】 前記第1の放射ビームおよび第2の放射ビームを生成する
    工程は、音響フィールドによる音響光学変調を前記初期ビームに行う工程を包含
    する、請求項81に記載の方法。
  84. 【請求項84】 前記音響光学変調を行う工程は、前記音響フィールドに位
    相変調を行う工程を包含する、請求項83に記載の方法。
  85. 【請求項85】 前記音響光学変調を行う工程は、前記音響フィールドを振
    幅変調工程を包含する、請求項83に記載の方法。
  86. 【請求項86】 前記ソースのうちの1つを該ソースのうちの他方のソース
    に対して移動させる工程は、音響フィールドに周波数変調を行うことによって前
    記初期ビームに音響光学変調を行う工程を包含する、請求項75に記載の方法。
  87. 【請求項87】 a)ある距離だけ隔離された2つの放射ソースであって、
    該ソースはそれぞれ、スペクトル分布を有し、該ソースのうちの他方のソースに
    対してコヒーレントである、放射ソースと、 b)該ソースの各々を該ソースのうちの他方のソースに対して移動させる制御シ
    ステムと、 c)オブジェクト表面上のポイントから散乱した放射を受信するような位置に配
    置された検出器と、 を備える、縞をオブジェクト表面上に投射する装置。
  88. 【請求項88】 前記検出器からの信号を受信するプロセッサをさらに備え
    、該プロセッサは、前記ソースの移動と、前記オブジェクト表面上のポイントか
    ら散乱した放射とに応答して位置情報を計算する、請求項87に記載の装置。
  89. 【請求項89】 前記制御システムは、前記放射ソースの1つからのスペク
    トル分布中のスペクトル成分の位相を、該放射ソースのうちの他方の放射ソース
    からのスペクトル分布中の各スペクトル成分の位相を前記オブジェクト表面上の
    ポイントにおいて測定した値に対して変化させる、請求項87に記載の装置。
  90. 【請求項90】 前記制御システムは、前記放射ソースのうちの1つの各ス
    ペクトル成分を、該放射ソースのうちの他方のソースの各スペクトル成分に対し
    て移動させて、該2つの放射ソースの各スペクトル成分間の距離を変化させる、
    請求項87に記載の装置。
  91. 【請求項91】 前記距離は、前記ソースのスペクトル領域の各々から等距
    離である中間点を含み、前記制御システムは、該ソースの各スペクトル成分を該
    ソースのスペクトル成分のうちの他方のスペクトル成分に対して移動させて、該
    距離を変化させ、該中間点を固定状態に保持する、請求項90に記載の装置。
  92. 【請求項92】 前記検出器と、前記オブジェクト表面上のポイントとの間
    に配置されたイメージングシステムをさらに備え、該イメージングシステムは、
    該検出器上のポイントをイメージングする、請求項87に記載の装置。
  93. 【請求項93】 前記検出器は複数の光検出器を備える、請求項87に記載
    の装置。
  94. 【請求項94】 前記プロセッサは、複数の処理ユニットを有するマルチプ
    ロセッサシステムを備え、前記複数の光検出器はそれぞれ、該複数の処理ユニッ
    トの各々と電気的に通信する、請求項93に記載の装置。
  95. 【請求項95】 前記2つの放射ソースは、 a)あるスペクトル幅のビーム放射の初期ソースと、 b)該ビーム放射の初期ソースと光学的に連絡し、第1の光学ビームおよび第
    2の光学ビームを生成するビーム分離器と、 c)該ビーム分離器と光学的に連絡するイメージングシステムであって、該2
    つの放射ソースがそれぞれ該第1の光学ビームおよび該第2の光学ビームに対応
    するように該2つの放射ソースを生成する、イメージングシステムと、 を備える、請求項87に記載の装置。
  96. 【請求項96】 前記ビーム分離器は回折格子である、請求項95に記載の
    装置。
  97. 【請求項97】 1つのソースのスペクトル分布中の波長におけるポイント
    は、他方ソースのスペクトル分布中の波長における各ポイントにおいて前記距離
    だけ隔離され、該距離は、該ポイントの波長に実質的に比例する、請求項87に
    記載の装置。
  98. 【請求項98】 前記制御システムは位相シフタをさらに備える、請求項8
    7に記載の装置。
  99. 【請求項99】 前記位相シフタは、前記ビーム分離器の位置を前記ビーム
    放射に対して変化させるトランスレータである、請求項98に記載の装置。
  100. 【請求項100】 前記制御システムは、焦点距離が可変であるズームレン
    ズを備え、該ズームレンズは、該可変な焦点距離の変化に応答して該隔離距離を
    変化させる、請求項87に記載の装置。
  101. 【請求項101】 前記制御システムは、焦点距離が異なる複数のレンズを
    備え、該隔離距離は、該複数のレンズのうちの1つの各焦点距離に対応する、請
    求項87に記載の装置。
  102. 【請求項102】 前記回折格子は空間光変調器である、請求項96に記載
    の装置。
  103. 【請求項103】 前記空間光変調器は、調節可能な格子周期を有し、該格
    子周期は、該空間光変調器に与えられる電気信号に応答する、請求項102に記
    載の装置。
  104. 【請求項104】 前記ビーム分離器は音響光学変調器である、請求項95
    に記載の装置。
  105. 【請求項105】 前記ビーム放射の初期ソースと光学的に連絡する振幅変
    調器をさらに備える、請求項95に記載の装置。
  106. 【請求項106】 オーダ−ブロックマスクをさらに備える、請求項95に
    記載の装置。
  107. 【請求項107】 前記オーダ−ブロックマスクは2重スリットのマスクで
    ある、請求項106に記載の装置。
  108. 【請求項108】 前記ビーム分離器と光学的に連絡する振幅変調器をさら
    に備え、該振幅変調器は、前景初期ソースの振幅を変調する、請求項95に記載
    の装置。
  109. 【請求項109】 表面を有するオブジェクト上において、該オブジェクト
    の表面上のポイントの三次元位置情報を判定する方法であって、 a)2つの放射ソースを提供する工程であって、該2つの放射ソースはそれぞ
    れ、空間スペクトル分布を有し、該2つのソースのうちの他方のソースに対して
    コヒーレントである、工程と、 b)該表面上のポイントに検出器を設ける工程と、 c)該オブジェクト表面上のポイントを該ソースの各々からの放射で照射する
    工程と、 d)該ソースの各々を互いに移動させる工程と、 e)該オブジェクト表面上のポイントにおいて該放射を検出する工程と、 f)該ソースの移動と、該オブジェクト表面上のポイントにおいて検出された
    放射とに応答して位置情報を計算する工程と、 を包含する、方法。
  110. 【請求項110】 前記2つの放射ソースは、狭帯域である空間スペクトル
    分布を有する、請求項109に記載の方法。
  111. 【請求項111】 前記ソースの1つからの放射の空間スペクトル分布中の
    スペクトル成分の位相を、該ソースのうちの他方のソースからの放射の空間スペ
    クトル分布中の各スペクトル成分の位相を前記オブジェクト表面上のポイントに
    おいて測定した値に対して変化させる工程をさらに包含する、請求項109に記
    載の方法。
  112. 【請求項112】 前記放射は仮想放射ソースである、請求項109に記載
    の方法。
  113. 【請求項113】 前記ソースの各々を移動させる工程は、該ソースのうち
    の1つのスペクトル成分と該ソースのうちの他方のソースの各スペクトル成分と
    間の距離を変化させて、該ソースの各スペクトル成分間の距離を変化させる工程
    を包含する、請求項109に記載の方法。
  114. 【請求項114】 前記距離は、前記ソースのスペクトル成分の各々から等
    距離である中間点を含み、該ソースの各々を移動させる工程は、該ソースのうち
    の1つのスペクトル成分と、該ソースのうちの他方のソースの各スペクトル成分
    と間の距離を変化させて、該ソースの各スペクトル成分間の距離を変化させ、該
    中間点を固定状態に保持する工程を包含する、請求項113に記載の方法。
  115. 【請求項115】 表面を有するオブジェクト上において、該オブジェクト
    の表面上のポイントの三次元位置情報を判定する装置であって、 a)空間スペクトル分布を有し、互いにコヒーレントな2つの放射ソースと、 b)該ソースの各々を相互に移動させる制御システムと、 c)該オブジェクト表面上のポイントを照射する放射を受信するように、該オ
    ブジェクト表面上の該ポイントに配置された検出器と、 d)該検出器からの信号を受信するプロセッサであって、該ソースの移動と、
    該オブジェクト表面上の該ポイントにおいて受信された放射とに応答して該オブ
    ジェクト表面上の該ポイントの位置情報を計算する、プロセッサと、 を備える、装置。
  116. 【請求項116】 前記2つの放射ソースは、狭帯域である空間スペクトル
    分布を有する、請求項115に記載の方法。
  117. 【請求項117】 前記制御システムは、前記ソースの1つからの放射の空
    間スペクトル分布中のスペクトル成分の位相を、該ソースのうちの他方のソース
    からの放射の空間スペクトル分布中の各スペクトル成分の位相を前記オブジェク
    ト表面上のポイントにおいて測定した値に対して変化させる、請求項115に記
    載の方法。
  118. 【請求項118】 前記放射ソースは仮想放射ソースである、請求項115
    に記載の装置。
  119. 【請求項119】 前記制御システムは、前記放射ソースのうちの1つの各
    スペクトル成分を、該放射ソースのうちの他方のソースの各スペクトル成分に対
    して移動させて、該2つの放射ソースの各スペクトル成分間の距離を変化させる
    、請求項115に記載の方法。
  120. 【請求項120】 前記距離は、前記ソースのスペクトル領域の各々から等
    距離である中間点を含み、前記制御システムは、該ソースの各スペクトル成分を
    該ソースのスペクトル成分のうちの他方のスペクトル成分に対して移動させて、
    該距離を変化させ、該中間点を固定状態に保持する、請求項119に記載の装置
JP2001552052A 2000-01-10 2001-01-03 表面輪郭測定のための装置および方法 Pending JP2003519786A (ja)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (1)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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