JP4782958B2 - 表面形状測定装置及びその方法、プログラム並びに記憶媒体 - Google Patents
表面形状測定装置及びその方法、プログラム並びに記憶媒体 Download PDFInfo
- Publication number
- JP4782958B2 JP4782958B2 JP2001289786A JP2001289786A JP4782958B2 JP 4782958 B2 JP4782958 B2 JP 4782958B2 JP 2001289786 A JP2001289786 A JP 2001289786A JP 2001289786 A JP2001289786 A JP 2001289786A JP 4782958 B2 JP4782958 B2 JP 4782958B2
- Authority
- JP
- Japan
- Prior art keywords
- data
- surface shape
- measured
- light
- measurement
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/2441—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures using interferometry
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Instruments For Measurement Of Length By Optical Means (AREA)
Description
【０００１】
【発明の属する技術分野】
この発明は、被測定物の表面形状を測定する表面形状測定装置及びその方法、プログラム並びに記憶媒体に関する。
【０００２】
【従来の技術】
物体の表面形状を光の波長レベルで測定可能な技術として、レーザ干渉計（例えば、Ｚｙｇｏ社のレーザ干渉システム等）がある。このレーザ干渉計は、参照光の位相に変調をかけながら、ＣＣＤなどの撮像素子を用いて複数の干渉縞画像を収録して、それらの干渉縞から被測定物の反射光の位相を求め、その位相データに基づいて光の波長レベルで物体形状を求めるものである。
【０００３】
しかしながら、物体が光源波長の１／２を超えるような段差を含む場合は、その段差の位置において、位相のアンラップができなくなり、これにより位相の関連性がとぎれるため、形状測定値に不確定性が含まれてしまうという不具合がある。
【０００４】
走査型白色干渉計（例えば、Ｚｙｇｏ社の三次元表面構造解析顕微鏡等）は、かかる不具合を解消するもので、複数の波長でそれぞれ求めた位相情報の相互関係から、単一波長での光軸方向の測定範囲λ／２を拡大する方法を用いている（多波長測定）。そして段差がさらに大きくなった場合には、参照光と物体反射光の光路長がほぼ一致したときにのみ白色干渉縞が発生する性質を利用して、干渉計、被測定物体のどちらかを干渉計の光軸方向に走査し、白色干渉縞の位置と、干渉計又は被測定物の移動量を検知することで、絶対測長を実施している。
【０００５】
しかし、このような走査型白色干渉計により大きな段差を測定する場合には、測定のために被測定物又は測定光学系を測定光学系の光軸方向に移動させる必要があるために、測定時間がかかるという不具合がある。また、測定中の移動であり、かつ、干渉計又は被測定物の移動量を検出する必要があるため、測定誤差要因が増え、また移動部の移動範囲内に光軸方向測定範囲が制限されるという不具合がある。
【０００６】
これに対し、特開2000-221013に開示の技術は、波長を超える段差や絶対距離を含んだ形状の測定を高精度かつ高速に測定ミスすることなく行なうため、走査型白色干渉計において、取得すべき干渉縞画像を低減し、測定の高速化を図ろうとしたものである。
【０００７】
また、特開2001-41724に開示の技術では、干渉計を用いた表面形状測定装置において、測定対象物を移動させることなく、測定対象物表面の三次元形状を高精度に測定できるようにするため、複屈折位相板を利用した位相シフト法を実施しており、測定対象物を移動させることなく三次元形状の測定を可能としている。また、位相板の代わりに電気光学素子を利用して電気的に位相シフトをさせれば、高速測定が可能としている。
【０００８】
さらに、物体をＣＣＤなどで観測した強度データをもとに、物体形状を求める技術として、depth from focus理論に基づく方法がある。これはレンズの焦点距離を変化させながら収録した複数の撮像画像をもとに、物体表面の合焦位置を求め、それにより表面形状を求めるものである（例えば、石原満宏、佐々木博美「合焦法による高速三次元形状計測」精密工学会誌，vol.63,No.1,1997を参照）。
【０００９】
【発明が解決しようとする課題】
しかしながら、特開2000-221013に開示の技術では、サブフリンジオーダーで形状を求めるのに、位相変調法（位相シフト法）を用いているため、装置内に機械的な駆動部分を有しており、また、複数の画像を取得する必要があるため、走査型白色干渉計よりは速いとはいえ、高速化に限界がある。そのため、動きを伴う物体の動的形状を測定することは不可能である。また、機械的な駆動を伴う測定であるために、測定誤差要因が増えるという不具合は依然として残る。
また、特開2001-41724に開示の技術では、原理的には位相シフト法を用いたレーザ干渉測定であり、上述したような光源波長の１／２を超える段差を含む面の測定には不確定性がある。また、高速化したといっても、複数の画像データを必要とするために高速化に限界があり、動的形状の測定は不可能である。
【００１０】
さらに、前記のdepth from focus理論に基づく技術では、一般にレーザ干渉に比べると測定分解能が低いため、サブフリンジオーダーの形状が要求される部品の評価には適さないという不具合がある。また、光学系の倍率を上げることで、測定分解能をある程度向上させることができるが、レーザ干渉による分解能と同じレベルまで上げようとすると、横方向の測定範囲（視野）が狭くなってしまうという不具合がある。またレンズの焦点距離を変えながら複数の撮像画像を収録する方法では、収録方法を工夫したとしても複数の画像を取り込むだけの時間がかかるうえ、レンズの焦点距離を正確に変化させるための構成が必要となり、装置構成が複雑化する。また、移動部分が存在するために誤差要因も増え、光軸方向の測定範囲も移動部の可動範囲に限定されるといった不具合がある。
【００１１】
この発明の目的は、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することである。
【００１２】
この発明の目的は、２次元的な振幅データではコントラストが低下して合焦点はずれが検知できない場合にも、被測定物の表面が平滑で、振幅データでは光学的にテクスチャを観測できない場合に、被測定面の形状測定ができるようにすることである。
【００１３】
この発明の目的は、被測定物の位置、姿勢を観察しやすくし、測定の操作性を向上させることである。
【００１４】
この発明の目的は、高速に全焦点画像を生成でき、より自然な被測定物の画像の観測が可能となるようにすることである。
【００１５】
この発明の目的は、より正確な全焦点画像の生成ができ、測定における被測定物の位置、姿勢の調整がより正確に実施可能となり、また、測定の高速化、高精度化を図ることができるようにすることである。
【００１６】
この発明の目的は、測定の高速化を図れるようにすることである。
【００１７】
この発明の目的は、被測定物の動きが高速な場合においても、正確に物体の動的形状を測定することができるようにすることである。
【００１８】
この発明の目的は、デフォーカス収差を除去して、干渉光学系の拡大倍率を高めた測定を正確に実施することである。
【００１９】
この発明の目的は、被測定物表面のより広い範囲を一度に測定できるようにすることである。
【００２０】
この発明の目的は、被測定物表面のより広い範囲を一度に測定することを容易に行なえるようにすることである。
【００２１】
【課題を解決するための手段】
請求項１に記載の発明は、被測定物に照射する光を発する光源と、その光の被測定物における反射光と所定の参照光との間で干渉縞を発生させる干渉光学系と、この干渉縞を撮像して当該干渉縞の画像データである干渉縞データを出力する撮像素子と、前記干渉縞データに基づいて前記反射光の光軸方向の複数位置における前記被測定物の反射光の位相を示す位相データ及び振幅を示す振幅データを求める光学データ検出手段と、前記振幅データの微分値の和である合焦測度を算出し、前記合焦測度が最大となる合焦距離に基づき、前記被測定物の第１の表面形状データを求める第１の表面形状データ作成手段と、前記位相データを長さ単位に変換することで、前記被測定物の第２の表面形状データを求める第２の表面形状データ作成手段と、前記第１の表面形状データと前記第２の表面形状データとを合成又は比較して、前記被測定物の表面形状のデータである第３の表面形状データを求める合成又は比較手段と、を備えている表面形状測定装置である。
【００２２】
したがって、被測定物の光源波長の１／２を超える段差については第１の表面形状データで受け持ち、それ以外の領域の形状測定については第２の表面形状データで受け持って、この両者を合成又は比較することにより、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【００３１】
請求項２に記載の発明は、請求項１に記載の表面形状測定装置において、前記干渉光学系は、その拡大倍率の２乗に反比例する感度で得られる前記振幅データによる前記被測定物の表面形状における前記反射光の光軸方向の測定分解能を、前記光源の波長の１／２より小さくなるように所定の拡大倍率に設定されている。
【００３２】
したがって、測定分解能を光源波長の１／２より小さく設定することにより、光源波長の１／２を超える段差も測定することができる。
【００３３】
請求項３に記載の発明は、請求項１又は２に記載の表面形状測定装置において、前記データ検出手段は、前記反射光の光軸方向の複数位置を、それぞれの位置の間隔が前記光源の波長の１／２より小さいものとしている。
【００３４】
したがって、反射光の光軸方向の複数位置は、それぞれの位置の間隔が光源の波長の１／２より小さくなるように設定しているので、光源波長の１／２を超える段差も測定することができる。
【００３５】
請求項４に記載の発明は、請求項１〜３のいずれかの一に記載の表面形状測定装置において、前記合成又は比較手段は、前記合成として、
ｍ×（λ／２）＋δ
（ｍ：次数，λ：光源波長，δ：前記第２の表面形状データ）における次数ｍを決定する処理を行なうものである。
【００３６】
したがって、第１と第２の表面形状データを的確に合成して、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【００３７】
請求項５に記載の発明は、請求項１〜４のいずれかの一に記載の表面形状測定装置において、前記光学データ検出手段は、前記位相データ及び振幅データをキャリヤ周波数を与えた単一の前記干渉縞の干渉縞データから求める。
【００３８】
したがって、観測する干渉縞にキャリヤ周波数を与えて、単一の干渉縞データから被測定物の反射光の振幅データ、位相データを計算することにより、測定の高速化を図ることができる。
【００３９】
請求項６に記載の発明は、請求項１〜５のいずれかの一に記載の表面形状測定装置において、前記光源は、前記被測定物にパルス光を照射するものである。
【００４０】
したがって、パルスの発光時間を被測定物の動きに対応して十分短くすることで、被測定物が光を反射する瞬間に被測定物がみかけ上静止した状態となり、被測定物からの反射光を安定させて、正確に反射光の振幅データ、位相データを求めることができるので、被測定物の動きが高速な場合においても、正確に物体の動的形状を測定することができる。
【００４１】
請求項７に記載の発明は、請求項１〜６のいずれかの一に記載の表面形状測定装置において、前記変換光学系は、前記反射光及び前記参照光をそれぞれ球面波に変換し、前記撮像素子による干渉縞の撮像位置での前記両球面波の曲率をほぼ一致させて前記干渉光学系の倍率を設定している。
【００４２】
したがって、デフォーカス収差を除去して、干渉光学系の拡大倍率を高めた測定を正確に実施することができる。
【００４３】
請求項８に記載の発明は、請求項１〜７のいずれかの一に記載の表面形状測定装置において、前記干渉光学系の拡大倍率を可変する拡大倍率可変装置を備え、前記第１の表面形状データ作成手段及び前記第２の表面形状データ作成手段のうち前者による前記被測定物の測定分解能又は前記反射光の光軸方向の測定範囲のみを設定できるようにしている。
【００４４】
したがって、光源波長の１／２より小さい測定分解能のもとで、被測定物の大きさに応じて干渉光学系の拡大倍率を変化させるようにすれば、被測定物表面のより広い範囲を一度に測定することができる。
【００４５】
請求項９に記載の発明は、請求項８に記載の表面形状測定装置において、前記拡大倍率可変装置は、前記撮像素子による前記干渉縞の撮像位置における前記反射光の曲率と前記参照光の曲率とがほぼ一致するよう連動させて曲率を変化させることにより、前記干渉光学系の倍率を可変としている。
【００４６】
したがって、被測定物表面のより広い範囲を一度に測定することを容易に行なうことができる。
【００４７】
請求項１０に記載の発明は、被測定物に光を照射したときの反射光と所定の参照光との間で発生させた干渉縞を撮像して当該干渉縞の画像データである干渉縞データを取得する撮像工程と、前記干渉縞データに基づいて前記反射光の光軸方向の複数位置における前記被測定物の反射光の位相を示す位相データ及び振幅を示す振幅データを求める光学データ検出工程と、前記振幅データの微分値の和である合焦測度を算出し、前記合焦測度が最大となる合焦距離に基づき、前記被測定物の第１の表面形状データを求める第１の表面形状データ作成工程と、前記位相データを長さ単位に変換することで、前記被測定物の第２の表面形状データを求める第２の表面形状データ作成工程と、前記第１の表面形状データと前記第２の表面形状データとを合成又は比較して、前記被測定物の表面形状のデータである第３の表面形状データを求める合成又は比較工程と、を含んでなる表面形状測定方法である。
【００４８】
したがって、被測定物の光源波長の１／２を超える段差については第１の表面形状データで受け持ち、それ以外の領域の形状測定については第２の表面形状データで受け持って、この両者を合成又は比較することにより、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【００５７】
なお、前記データ検出工程は、前記反射光の光軸方向の複数位置を、それぞれの位置の間隔が前記光源の波長の１／２より小さいものとすることができる。
【００５８】
これにより、反射光の光軸方向の複数位置は、それぞれの位置の間隔が光源の波長の１／２より小さくなるように設定しているので、光源波長の１／２を超える段差も測定することができる。
【００５９】
前記合成又は比較工程は、前記合成として、
ｍ×（λ／２）＋δ
（ｍ：次数，λ：光源波長，δ：前記第２の表面形状データ）
における次数ｍを決定する処理を行なうものである。
【００６０】
これにより、第１と第２の表面形状データを的確に合成して、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【００６１】
前記光学データ検出工程は、前記位相データ及び振幅データをキャリヤ周波数を与えた単一の前記干渉縞の干渉縞データから求めることができる。
【００６２】
これにより、観測する干渉縞にキャリヤ周波数を与えて、単一の干渉縞データから被測定物の反射光の振幅データ、位相データを計算することにより、測定の高速化を図ることができる。
【００６３】
前記撮像工程は、前記被測定物に照射する光をパルス光とすることができる。
【００６４】
これにより、パルスの発光時間を被測定物の動きに対応して十分短くすることで、被測定物が光を反射する瞬間に被測定物がみかけ上静止した状態となり、被測定物からの反射光を安定させて、正確に反射光の振幅データ、位相データを求めることができるので、被測定物の動きが高速な場合においても、正確に物体の動的形状を測定することができる。
【００６５】
請求項１１に記載の発明は、被測定物に光を照射したときの反射光と所定の参照光との間で発生させた干渉縞を撮像した画像データである干渉縞データに基づいて前記反射光の光軸方向の複数位置における前記被測定物の反射光の位相を示す位相データ及び振幅を示す振幅データを求める光学データ検出処理と、前記振幅データの微分値の和である合焦測度を算出し、前記合焦測度が最大となる合焦距離に基づき、前記被測定物の第１の表面形状データを求める第１の表面形状データ作成処理と、前記位相データを長さ単位に変換することで、前記被測定物の第２の表面形状データを求める第２の表面形状データ作成処理と、前記第１の表面形状データと前記第２の表面形状データとを合成又は比較して、前記被測定物の表面形状のデータである第３の表面形状データを求める合成又は比較処理と、をコンピュータに実行させるコンピュータに読み取り可能なプログラムである。
【００６６】
したがって、被測定物の光源波長の１／２を超える段差については第１の表面形状データで受け持ち、それ以外の領域の形状測定については第２の表面形状データで受け持って、この両者を合成又は比較することにより、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【００７５】
なお、前記データ検出処理は、前記反射光の光軸方向の複数位置を、それぞれの位置の間隔が前記光源の波長の１／２より小さいものとすることができる。
【００７６】
これにより、反射光の光軸方向の複数位置は、それぞれの位置の間隔が光源の波長の１／２より小さくなるように設定しているので、光源波長の１／２を超える段差も測定することができる。
【００７７】
前記合成又は比較処理は、前記合成として、
ｍ×（λ／２）＋δ
（ｍ：次数，λ：光源波長，δ：前記第２の表面形状データ）
における次数ｍを決定する処理を行なうものである。
【００７８】
これにより、第１と第２の表面形状データを的確に合成して、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【００７９】
前記光学データ検出処理は、前記位相データ及び振幅データをキャリヤ周波数を与えた単一の前記干渉縞の干渉縞データから求めることができる。
【００８０】
これにより、観測する干渉縞にキャリヤ周波数を与えて、単一の干渉縞データから被測定物の反射光の振幅データ、位相データを計算することにより、測定の高速化を図ることができる。
【００８１】
前記撮像処理は、前記被測定物に照射する光をパルス光とすることができる。
【００８２】
これにより、パルスの発光時間を被測定物の動きに対応して十分短くすることで、被測定物が光を反射する瞬間に被測定物がみかけ上静止した状態となり、被測定物からの反射光を安定させて、正確に反射光の振幅データ、位相データを求めることができるので、被測定物の動きが高速な場合においても、正確に物体の動的形状を測定することができる。
【００８３】
請求項１２に記載の発明は、請求項１１に記載のプログラムを記憶した記憶媒体である。
【００８４】
したがって、請求項１１に記載の発明と同様の作用、効果を奏する。
【００８５】
【発明の実施の形態】
この発明の一実施の形態について説明する。
【００８６】
図１は、この発明の一実施の形態である表面形状測定装置の概略構成を示すブロック図である。図１に示すように、この表面形状測定装置は、光源となるＨｅ−Ｎｅレーザ１と、その出力するレーザ光の被測定物１４への照射光強度を調整するためのＮＤフィルタ２と、レーザ光を拡大するためのビームエキスパンダ３とを備えている。ビームエキスパンダ３はスペイシャルフィルタとしても作用する。
【００８７】
ビームエキスパンダ３にて拡大された光は、ミラー４にて折り返され、ビームスプリッタ５に入射する。ビームスプリッタ５で反射した光は、光を球面波に変換するためのレンズ６を通り、ミラー７で折り返され、光強度を調整するためのＮＤフィルタ８、ハーフミラー９を介して撮像素子であるＣＣＤ１０に入射する。この光は、被測定物１４の反射光と干渉するための参照光となる。
【００８８】
一方、ビームスプリッタ５を透過した光は、ミラー１１にて折り返され、レンズ１２により球面波に変換され、ハーフミラー９を透過して、対物レンズ１３によりほぼ平行光に変換されて、被測定物１４に照射する。被測定物１４で反射した光は、対物レンズ１３を通過し、ハーフミラー９で反射して物体光としてＣＣＤ１０に到達し、前記の参照光との間に干渉をおこして干渉縞を発生する。したがって、これらの光学素子により干渉光学系を実現している。この干渉縞はＣＣＤ１０にて撮像され、フレームグラバー１５を介してコンピュータ１６に転送され、コンピュータ１６に記録される。
【００８９】
ビームスプリッタ５で分岐された後の参照光路と物体光路の長さは、光源であるＨｅ−Ｎｅレーザ１のコヒーレンス長以下になるように設定されている。また、被測定物１４からの反射光がＣＣＤ１０で撮像するのに適した強度になるように、ＮＤフィルタ２により被測定物１４への照射光強度を調整してあり、さらに、被測定物１４からの反射光と参照光との干渉により発生する干渉縞のコントラストが高くなるように、ＮＤフィルタ８により参照光強度が調整されている。
【００９０】
対物レンズ１３は被測定物１４の反射光を球面波に変換するように作用し、これにより後述のように被測定物１４の反射光の複素振幅（被測定物１４の像）を再生したときに、像が拡大して観測される。レンズ１２と対物レンズ１３の位置関係は、被測定物１４にほぼ平行光が照射されるように調整されている。レンズ６は、参照光を球面波に変換するように作用し、対物レンズ１３により球面波に変換された被測定物１４の反射光の曲率と、参照光の曲率とが、ＣＣＤ１０の撮像位置においてほぼ一致するように、位置が調整されている。
【００９１】
被測定物１４の像を拡大しない場合は、対物レンズ１３は不要であり、それに伴い、レンズ６、レンズ１２も不要である。また後で詳細を説明するが、ここでは干渉縞をホログラムとして扱い、ＣＣＤ１０にて撮像したホログラム干渉縞に仮想的に参照光を照射して、その回折光波により被測定物１４からの反射光（被測定物１４の像）の複素振幅を再生する。単一の干渉縞データをＣＣＤ１０にて撮像し、被測定物１４の反射光を再生する場合、ホログラムを透過する０次回折光と実像と虚像の３つが再生されるため、この３つを分離して再生できるように、被測定物１４からの反射光の光軸と参照光の光軸との間に適当な傾きを与えて、干渉縞にキャリヤ周波数がのるようにしている。前記の両光軸に傾きを与えるためには、例えば、ミラー７、被測定物１４又はハーフミラー９かの傾きを調整すればよい。
【００９２】
次に、被測定物１４からの反射光（被測定物１４の像）の複素振幅（この明細書において、光の振幅と位相を併せて複素振幅と呼ぶ）を再生計算する方法について説明する。
【００９３】
まず、被測定物１４からの反射光と参照光との干渉により発生する干渉縞のＣＣＤ１０の撮像位置での強度は、（１）式で与えられる。
【００９４】
【数１】
【００９５】
ここで、Ｉは干渉縞の強度であり、ｘ，ｙはＣＣＤ１０の撮像位置のｘ，ｙ座標、Ｒは参照光の複素振幅であり、Ａ，φはそれぞれ被測定物１４の反射光の振幅と位相である。
【００９６】
干渉縞をＣＣＤ１０にて撮像し、ホログラム画像データとしてコンピュータ１６に記憶する。乾板を用いたホログラフィでは、干渉縞を記録した乾板（ホログラム）に改めて参照光を照射すると、参照光がホログラムに記録された干渉縞にて回折し、それが記録時の物体反射光として振舞うので、これにより、物体反射光（物体像）が再生されることになった。ここでは、実際に参照光を照射せずに、仮想的に参照光がホログラムに照射されたものとして、ＣＣＤ１０にて記録したホログラム画像からフレネル近似のもとに物体反射光（物体像）を再生する。ホログラムに参照光として平行光が照射されたとすると、ホログラムによる回折光、すなわち物体反射光は（２）式のように表される。
【００９７】
【数２】
【００９８】
ここで、Ｕは再生距離ｄだけ離れた像面位置での物体反射光波の複素振幅、ｘ’，ｙ’は像面位置でのｘ，ｙ座標、ｃは複素定数、λは光源波長を表す。
【００９９】
したがって、再生距離ｄを入力し、収録した干渉縞データを用いて（２）式を計算することにより、再生距離ｄにおける物体反射光の複素振幅が再生されることになる。被測定物１４の像を拡大しない場合は、再生距離ｄにＣＣＤ１０の撮像位置と被測定物１４の表面との距離を代入すればよい。図１に示した装置のように、被測定物１４の像を拡大する場合は、対物レンズ１３による被測定物１４の拡大像面とＣＣＤ１０の撮像位置との距離を、再生距離ｄとして入力する。
【０１００】
図２は、図１における被測定物１４、対物レンズ１３、ホログラム面（ＣＣＤ１０の撮像位置）、対物レンズ１３による被測定物１４の拡大像それぞれの位置関係を示す説明図である。図２において、符号１７は被測定物１４の位置を示し、符号１８は対物レンズ１３の位置、符号１９はホログラム面（ＣＣＤ１０の撮像位置）、符号２０は対物レンズ１３による被測定物１４の拡大像面である。
【０１０１】
図２における距離ｄ’が再生距離となり、それを（２）式に入力して、被測定物１４の像を再生計算する。図２の距離ｓ及び距離ｓ’はレンズ結像の下記（３）式の関係にあるため、対物レンズ１８の焦点距離ｆと距離ｓをあらかじめ測定しておくことにより、距離ｓ’を求めることができる。そして、対物レンズ１３（位置１８）からほぼ距離ｓ’の位置にＣＣＤ１０の撮像位置を設置すれば、距離ｄ’はほぼ０となり、目視ではほぼ焦点のあった被測定物１４の像が得られ、被測定物１４の位置などの調整が容易になる。また、距離ｄ’が大きくなると、再生した像の周辺部に歪が生じる場合があるため、歪の少ない再生した像を得るためにも効果がある。
【０１０２】
【数３】
【０１０３】
再生像には、ホログラムを透過する光（０次回折光）による像、実像、虚像の３つの像が含まれていて、被測定物１４の形状を求めるためには、そのうちの実像に注目する。撮像時において物体反射光と参照光との光軸が平行に近いと、その３つの像が重なって再生されるため、前記したように、物体反射光と参照光の光軸に相対的な傾きを与えて、干渉縞にキャリヤ周波数をのせた状態でホログラム（干渉縞画像）を取り込み、これにより像面位置でそれらの３つの像を離れた位置に再生することができる。
【０１０４】
（２）式で表される複素振幅の実部を“Real｛Ｕ（ｘ’，ｙ’）｝”とし、虚部を“Imaginary｛Ｕ（ｘ’，ｙ’）｝”とすると、被測定物１４の反射光の振幅Ａ、位相φは、それぞれ次の（４）、（５）式にて与えられる。
【０１０５】
【数４】
【０１０６】
（４）式にて計算される振幅Ａと、（５）式にて計算される位相φが、それぞれ、振幅データ、位相データに相当する。この位相データの示す位相φに、“λ／（４π）”を乗じて長さ単位に変換し、位相の飛びがある場合、それをアンラップすることにより、被測定物１４の物体形状を光源の波長λ以下の精度で求めることができる。
【０１０７】
しかしながら、被測定物１４の物体面に、レーザ光の波長λに対し、λ／２を超えるような大きな段差があると、その段差位置で位相φ（形状）の関連性が途切れるため、測定値に不確定性が生じてくる。
【０１０８】
そのため、この形状測定装置では、（４）式による振幅データを活用することによっても、被測定物１４の表面形状を測定し、振幅データによる表面形状の測定値と位相データによる表面形状の測定値とを比較し、又は、合成することによって、λ／２を超える段差を含む被測定面であっても、表面形状の測定を可能なものとしている。かかる点について以下に説明する。
【０１０９】
まず、振幅データを用いた表面形状の測定について説明する。簡単のため、被測定物１４を、光学系を用いて拡大しない場合における被測定物１４、ホログラム面（ＣＣＤ１０の撮像位置）、及び、再生像の位置関係を表す図３を参照して説明する。図３において、符号２１は被測定物１４の位置、符号２２はホログラム面（ＣＣＤ１０の撮像位置）、符号２３は被測定物１４の再生像を示す。
【０１１０】
ホログラム（ＣＣＤ１０の撮像位置）を原点とした場合、再生距離ｄを−ｘ１とすると、ｘ１の位置にある面の実像が焦点のあった状態で再生され、同様に、再生距離ｄを−ｘ２，−ｘ３とすると、それぞれｘ２，ｘ３の位置にある面の実像が焦点のあった状態で再生されることを示している。その場合、“ｄ＝−ｘ１”とした場合には、ｘ２，ｘ３にある面は焦点ずれすることになる。したがって焦点があっているかどうかを検出することにより、ｘ１，ｘ２，ｘ３の相対的な距離の差（すなわち段差）を求めることができる。
【０１１１】
その焦点が合っているかどうかの判定のために、ここでは、（４）式で求めた振幅データを用いて、depth from focus理論に基づいた演算を実施する。depth from focus理論に基づく演算方法としてはさまざまな種類のものがあるが、例えば、“石原満宏、佐々木博美「合焦法による高速三次元形状計測」精密工学会誌，vol.63,No.1,1997”に示された方法を用いることができる。
【０１１２】
すなわち、撮像した干渉縞データを用いて、再生距離ｄを微小量（Δｄ）ずつ変えながら、複数の複素振幅を（２）式を用いて計算し、各々の再生距離ｄにおける振幅データを（４）式により計算する。そして各々の振幅データにおいて、ＣＣＤ１０の撮像面の任意の小領域Ｌにおける振幅値Ａの微分値の和を合焦測度ｖとして求め、これらの合焦測度ｖのうちピークを示す再生距離ｄを、焦点のあった位置と判定する。
【０１１３】
図４に、撮像画像の任意の小領域Ｌにおいて、再生距離ｄを変化させたときの合焦測度ｖの変化を示すが、図４の距離Ｄが焦点のあう再生距離ｄ（合焦距離Ｄ）となる。したがって、任意の小領域Ｌにおいて振幅値Ａの微分値の和がピークになる再生距離ｄを求め、それを画像全体にわたって繰り返すことにより、被測定物１４の表面全体の各小領域Ｌにおける合焦距離Ｄが計算できる。そして基準位置を任意の小領域にとり、基準位置の再生距離との差を求めることによって、被測定物１４全体の表面形状を求めることができる。よって、撮像画像の各画素について、被測定物１４の表面形状データＳijが得られる。ここで、ｉ，ｊは、撮像画像の画素単位で表わされる座標を示す。
【０１１４】
合焦測度ｖを求めるとき、再生距離ｄを少しずつ変化させながら、被測定物１４の反射光の略光軸方向における複数位置での被測定物１４の反射光の超えるを計算により求めるが、その場合、再生距離ｄの変化の間隔（ピッチ）Δｄは、光源波長λの１／２より小さいことが望ましい。ただし、図４に示したように、複数の合焦測度ｖのプロットをガウス関数などの関数に近似し、関数の頂点を求める作業により合焦距離Ｄを求めてもよく、その場合は、前記の再生距離ｄの変化の間隔Δｄは、厳密に光源波長λの１／２より小さくなくてもよい。
【０１１５】
ここで、振幅データを用いた測定の縦方向の感度γは、定数ｋと光学系の倍率βを用いて、
【０１１６】
【数５】
【０１１７】
と表わすことができ、光学系の倍率の２乗に反比例して分解能が上がる。
【０１１８】
一方、光源波長λの１／２の段差を、振幅データを用いる方法により測定しようとすると、光源波長λの１／２より小さい測定分解能が、振幅データを用いる方法に要求される。
【０１１９】
光の波長は、可視光の場合は数百ｎｍであるため、振幅データを用いた測定の測定分解能はｎｍオーダーでなければならない。それを考慮すると、被測定物１４の像を拡大せずに振幅データを用いた測定分解能の要求を満たそうとするのは困難であると考えられる。そこで、図１の装置構成では、対物レンズ１３を用いて被測定物１４の像を拡大して測定する装置構成を示した。ただし、光源であるＨｅ−Ｎｅレーザ１の波長をμｍやｍｍのオーダーまで大きくした場合などにおいては、被測定物１４の像を拡大せずとも、振幅データを用いた測定の測定分解能の要求は満足でき、測定が可能となる。
次に、（５）式により計算した位相データを用いて、被測定物１４の表面におけるλ／２を下回る微小形状を演算する方法について説明する。以下では、図１における、レンズ６、レンズ１２、対物レンズ１３を用いずに被測定物１４の像を拡大しない場合について説明する。
【０１２０】
すなわち、任意の小領域において、振幅データを用いる方法で焦点のあった再生距離ｄで再生した複素振幅を用い、（４）式により位相データを計算する。そして計算した位相φに“λ／４π”を乗じて長さ単位に変換して、被測定物１４の表面形状を求める。この場合も撮像画像の各画素について、被測定物１４の表面形状データが得られる。
図１に示した装置のように、対物レンズ１３を用いて被測定物１４の像を拡大した場合において、被測定物１４の反射光の曲率と参照光の曲率が一致していないと、前記の曲率の差に起因して、位相φの測定値に図５に示すような同心円状のデフォーカス収差が含まれてくる。被測定物１４の表面形状を正確に求めるためには、それを除去しなければならない。
【０１２１】
図１においてレンズ６の位置を調整して、被測定物１４の反射光の曲率と参照光の曲率をほぼ一致させると、両者の曲率の差はなくなるため、それにより位相データにデフォーカス収差は含まれなくなる。したがって任意の小領域Ｌにおいて、振幅データを用いる方法で焦点のあった再生距離ｄで再生した複素振幅を用い、（５）式により位相データを計算する。そして、この計算した位相φに“λ／４π”を乗じて長さ単位に変換して、被測定物１４の表面形状を求める。この場合も、撮像画像の各画素について被測定物１４の表面形状データＴijが得られる。ここで、ｉ，ｊは画像の画素単位で表わされる座標を示す。
【０１２２】
次に、振幅データにより求めた形状データＳijと、位相データにより求めた形状データＴijとを合成する方法について説明する。まず、振幅データより得られた形状データＳijの値をλ／２で除算し、その商ｍを求める。そして、形状データＳijの値にｍとλ／２の積を加えたものを、形状データＺijの値に代入する。この処理を、撮像画像の各画素にわたって実施することにより、被測定物１４のλ／２を超える段差を含む表面であっても、表面形状の正しい測定値を得ることができる。
【０１２３】
以上の処理は、言いかえると、次のように表わすこともできる。すなわち、光源の波長λの１／２を超える段差を含む面の正しい形状データをＺijとし、光源波長をλ、次数をｍij、位相データにより求めた形状をδijとすると、次の（７）式に示す関係が成り立つので、撮像画像の各画素において次数ｍijを求めればよい。
【０１２４】
【数６】
【０１２５】
以上のようにして測定を行なうことにより、被測定物１４や光学系を測定中に移動させることなく、光源波長λの１／２を超えるような段差を含む面をサブフリンジオーダーで高精度に測定可能となる。
【０１２６】
次に、以上説明した方法を用いて、コンピュータ１６が行なう形状データＺijの測定処理について概要を整理して説明する。
【０１２７】
まず、コンピュータ１６の構成について説明する。図６は、コンピュータ１６の電気的な接続を示すブロック図である。図６に示すように、コンピュータ１６は、各種演算を行ないコンピュータ１６の各部を集中的に制御するＣＰＵ３１と、各種のＲＯＭ、ＲＡＭからなるメモリ３２とが、バス３３で接続されている。
【０１２８】
バス３３には、所定のインターフェイスを介して、ハードディスクなどの磁気記憶装置３４と、マウス、キーボード等により構成される入力装置３５と、ＬＣＤ、ＣＲＴ等の表示装置３６と、光ディスクなどの記憶媒体３７を読み取る記憶媒体読取装置３８とが接続され、また、インターネットなどのネットワーク３９と通信を行なう所定の通信インターフェイス４０が接続されている。なお、記憶媒体３７は、この発明の記憶媒体を実施するものであり、ＣＤ，ＤＶＤなどの光ディスク、光磁気ディスク、フロッピーディスクなどの各種方式のメディアを用いることができる。また、記憶媒体読取装置３８は、具体的には記憶媒体３７の種類に応じて光ディスク装置、光磁気ディスク装置、フロッピーディスク装置などが用いられる。
【０１２９】
磁気記憶装置３４には、この発明のプログラムを実現する測定演算プログラムが記憶されている。この測定演算プログラムは、記憶媒体３７に記憶されていたものを記憶媒体読取装置３８により読み取るか、あるいは、インターネットなどのネットワーク３９からダウンロードするなどして、磁気記憶装置３４にインストールしたものである。このインストールによりコンピュータ１６は、形状データＺijの測定処理の実行が可能な状態となる。この測定演算プログラムは、所定のＯＳ上で動作するものであってもよい。
【０１３０】
以下では、測定演算プログラムに基づいてコンピュータ１６が行なう形状データＺijの測定処理の内容について、図７のフローチャートを参照して説明する。
【０１３１】
図７に示すように、まず、ＣＰＵ３１は、ＣＣＤ１０にて撮像された干渉縞の干渉縞データを取り込んで、メモリ３２に記憶し（ステップＳ１）、メモリ３２の所定領域に記憶される所定のカウント値ｋを０にリセットし（ステップＳ２）、干渉縞データを用いて以下の処理を行なう。ステップＳ１により、撮像工程を実現している。
【０１３２】
まず、再生距離ｄの初期値ｄ０に対し、再生距離ｄの微小量Δｄにカウント値ｋを乗算した値を加算し（ステップＳ３）、フレネル回折計算を行って、“再生距離ｄ０＋Δｄ×ｋ”での複素振幅を取得する（ステップＳ４）。ステップＳ３，Ｓ４により光学データ検出手段、光学データ検出工程、光学データ検出処理を実現している。そして、その複素振幅の振幅データにおける微小領域Ｌでの合焦測度を求め（ステップＳ５）、カウント値ｋを＋１だけインクリメントする（ステップＳ６）、という処理を、カウント値ｋが所定値Ｋになるまで繰り返す（ステップＳ７のＮ）。
【０１３３】
カウント値ｋが所定値Ｋになったときは（ステップＳ７のＹ）、各微小領域Ｌでの合焦距離を求めて、全微小領域Ｌでの合焦距離Ｄを求め（ステップＳ８）、合焦距離Ｄから表面形状データＳijを求める（ステップＳ９）。ステップＳ８，Ｓ９により第１の表面形状データ作成手段、第１の表面形状データ作成工程、第１の表面形状データ作成処理を実現している。
【０１３４】
次に、各画素での次数ｍijを計算し（ステップＳ１０）、微小領域Ｌにおける合焦距離Ｄを再生距離ｄとして（ステップＳ１１）、フレネル回折計算を行い、各微小領域Ｌにおける複素振幅を取得する（ステップＳ１２）。ステップＳ１０〜Ｓ１２により、光学データ検出手段、光学データ検出工程、光学データ検出処理を実現している。
【０１３５】
そして、その複素振幅の位相データにより、各画素での位相φの値を表面形状データＴijに変換する（ステップＳ１３）。ステップＳ１３により、第２の表面形状データ作成手段、第２の表面形状データ作成工程、第２の表面形状データ作成処理を実現している。
【０１３６】
最後に、以上のようにして求めた表面形状データＳijと表面形状データＴijとを、前記のように合成又は比較して、表面形状データＺijを求める（ステップＳ１４）。ステップＳ１４により、合成又は比較手段、合成又は比較工程、合成又は比較処理を実現している。
【０１３７】
なお、以上の処理では、振幅データを用いた合焦法を実施した結果、得られた各微小領域Ｌにおける合焦距離Ｄを用いて、フレネル回折計算を繰り返し、各微小領域Ｌでの複素振幅を再計算して（ステップＳ１２）、その位置での位相データを再度求めている。正確に計算するためにはそのほうが望ましいが、フレネル回折計算には時間がかかるため、全画素で位相データを計算するのに多大な時間を要する。そこで、求めた合焦距離Ｄが初期位置に対して微小な場合は、初期位置での位相データを代用しても誤差が小さいため、求めた合焦距離Ｄの値によっては、ステップＳ１２のフレネル回折計算を実施しなくてもよい。これにより演算時間、測定時間の短縮が図れる。
【０１３８】
図１の表面形状測定装置では、キャリヤ周波数をのせた干渉縞を用いて被測定物１４の反射光の複素振幅を求めたが、キャリヤ周波数をのせない同軸の干渉縞においても、参照光、被測定物１４の反射光の位相を変調することによって、被測定物１４の反射光の複素振幅を独立して求めることができる。この場合は、このような複素振幅を用いてもよいし、また、光源であるＨｅ−Ｎｅレーザ１のレーザ光の波長λを大きくした場合は、直接ＣＣＤ１０で被測定物１４の反射光の複素振幅を求められるようになるため、このようにして直接取得した複素振幅を用いてもよい。
【０１３９】
図１の表面形状測定装置によると、ＣＣＤ１０の撮像時間より遅い速度で動く物体に対しては、その動いている瞬間の動的形状の測定が可能になる。ただし、物体の動きが速くなり、ＣＣＤ１０の撮像時間より速くなると、撮像している最中に被測定物１４からの反射光の強度や位相が変化するため、正確に反射光の複素振幅を求めることができなくなる。
【０１４０】
この場合は、被測定物１４の動く速度に対して、十分短い時間で被測定物１４にパルス光を照射し、その反射光を受光するようにすると、光の照射中は、被測定物１４は見かけ上静止している状態とみなせるため、被測定物１４の反射光が安定し、正確にその複素振幅を求めることができて、正確な被測定物１４の動的表面形状測定が可能になる。
【０１４１】
被測定物１４にパルス光を照射する光源には、ルビーレーザやＹＡＧレーザのようにパルス光を出力する固体レーザや、パルス光を出力する半導体レーザを用いることができ、また、ＣＷ光を出力する半導体レーザをパルス変調駆動してパルス光を生成してもよいし、Ｈｅ−ＮｅレーザやアルゴンレーザからのＣＷ光を、回転チョッパを用いてパルス化して用いてもよい。
【０１４２】
また、パルス光を被測定物１４に照射するのに代えて、図１の表面形状測定装置におけるＣＣＤ１０を高速撮像が可能な高速カメラに置き換えることにより、光源をパルス化するのと同様の効果を得ることもできる。
【０１４３】
さらに、（３）式の関係において、光学系倍率は“距離ｓ’／距離ｓ”で決まるため、対物レンズ１３の位置や焦点距離を変えることによって、光学系倍率を変えることができる。光学系倍率を変えることにより、振幅データを用いた表面形状測定方法の測定分解能を変化させることができるため、被測定物１４の大きさに応じて光学系倍率を変化させるようにして、被測定物１４の表面のより広い範囲（視野）を一度に測定できるようにすることができる。
【０１４４】
図１の装置構成に代えて、図８のような装置構成としてもよい。図８において、図１と同一符号の部材は図１の装置と同様であるため、詳細な説明は省略する。
【０１４５】
すなわち、図８の装置構成が図１のものと相違するのは、被測定物１４を光学系の略光軸方向に移動させる拡大倍率可変装置であるステージ２４を設け、これにより被測定物１４と対物レンズ１３との間隔を変化させるようにしている点である。
【０１４６】
被測定物１４と対物レンズ１３との間隔が変化すると（３）式における距離ｓの値が変化し、それに合わせて（３）式が成り立つように距離ｓ’が変化するため、もって光学系の倍率を変化させることができる。
【０１４７】
また、図８の構成は、レンズ６を略光軸方向に移動させるための拡大倍率可変装置であるステージ２５を設けている点も相違する。すなわち、ステージ２４の移動に伴う被測定物１４と対物レンズ１３との間隔の変化に合わせてステージ２５を移動させ、ＣＣＤ１０の撮像位置での被測定物１４の反射光と参照光の曲率をほぼ一致させる。これにより、干渉縞データから計算した位相データにデフォーカス収差（図５参照）が含まれずに、光学系の倍率を変化させることが可能になる。また、距離ｓ’の変化にあわせて、ＣＣＤ１０の位置を変化させるためのステージを設けるようにしてもよい。
【０１４８】
さらに、図１の装置構成に代えて、図９のような装置構成とすることも考えられる。図９において、図１と同一符号の部材は図１の装置と同様であるため、詳細な説明は省略する。
【０１４９】
すなわち、図９の装置構成が図１のものと相違するのは、被測定物１４の反射光を球面波に変換する可変焦点レンズ２６と、参照光を球面波に変換する可変焦点レンズ２７と、レンズ１２をほぼ光軸方向に移動させるための拡大倍率可変装置であるステージ２８とを備えている点である。
【０１５０】
そして、可変焦点レンズ２６とレンズ１２との間にアフォーカル系の関係が成り立つように、可変焦点レンズ２６の焦点距離の変化に応じて、ステージ２８でレンズ１２を移動させることで、可変焦点レンズ２６から被測定物１４に照射される光がほぼ平行化されるようになっている。
【０１５１】
可変焦点レンズ２６又は２７のどちらかの焦点距離を変え、球面波の曲率を変化させることにより、拡大倍率を変化させることができるが、被測定物１４の反射光と参照光のＣＣＤ１０の撮像位置での曲率に差が生じると、それに応じて、干渉縞から計算した位相データにデフォーカス収差が含まれてしまう。そのため、被測定物１４の反射光と参照光との間にＣＣＤ１０の撮像位置での曲率に差が生じないように、可変焦点レンズ２６、２７の焦点距離の変化を連動させるようにする。
【０１５２】
これにより、干渉縞データから計算した位相データにデフォーカス収差が含まれることなく、光学系の倍率を変化させることが可能になる。また、可変焦点レンズ２６又は２７を固定焦点レンズにして、それをほぼ光軸方向に移動可能なステージに搭載して移動させることによって、被測定物１４の反射光と参照光の曲率の差が生じないように調整してもよい。さらに、対物レンズ１３の焦点距離を変えたことにより、（３）式の関係から距離ｓが変化するが、距離ｓ’の変化に応じてＣＣＤ１０の位置を変化させるためのステージを設けてもよい。
【０１５３】
ところで、以上説明した例では、合焦法の原理により求めた表面形状データと、干渉計測の原理により求めた表面形状データＴijとを合成又は比較することによって表面形状データＺijを求めて、光源波長λの１／２を超えるような段差を含む被測定物１４の表面形状をサブフリンジオーダーで測定できるようにしている。
【０１５４】
しかし、合焦法の原理により表面形状データを求める手段としては、前記の例に限定されるものではない。すなわち、前記の例では合焦法の原理により表面形状データＳijを求めるために振幅データを用いたが、これに代えて位相データを２次元画像データとして捉え、これを用いて合焦法を実施してもよい。
【０１５５】
合焦法では、物体像のボケ具合（合焦測度）をコントラストの低下として検知する。しかし、被測定面が鏡面のような平滑面であると、被測定物１４の表面の明確なテクスチャ（模様）を観察できない場合がある。そして、テクスチャが観察できないと、像のボケによるコントラストの低下が検出できなくなるため、合焦法が実施できず、表面形状測定が不可能になる。
【０１５６】
このような場合、振幅データに比べて比較的物体像をコントラストよく観察できる位相データ（位相像）を２次元強度データとして捉えて用いることができる。すなわち、平滑な面であっても、その表面を光学的に拡大していけば、表面粗さが光学系の空間解像力を下回らない限り表面粗さのテクスチャが観察されるようになる。そのため、その像をコントラストの高い位相像として観察することによって、合焦法による形状測定が可能となる。
【０１５７】
そのために、コンピュータ１６で行なう処理としては、図７のフローチャートにおいて、振幅データの微小領域Ｌで合焦法を行なうのに（ステップＳ５）代えて、位相データの微小領域Ｌにおいて合焦法を行なえばよい。
【０１５８】
また、合焦法により表面形状測定を実施する場合には、焦点の合う画像領域を被測定撮像位置で探索していくこととなる。そのため、焦点の合った画像領域をつなぎ合わせることで、全視野で焦点の合った全焦点画像を生成することができる。顕微鏡などのように拡大倍率が大きい光学系では焦点深度が浅くなり、そのため、被測定物１４の位置や姿勢の調整がしづらくなる。このような場合に、全焦点画像を生成することにより、調整作業が容易になる。
【０１５９】
全焦点画像の生成の方法としては、フレネル回折計算を繰り返して複数のデフォーカス画像を取得した後、画像における各微小領域Ｌで合焦距離Ｄを求める。そして、各微小領域Ｌにおいて、求めた合焦距離Ｄと最も近い距離で得られたデフォーカス画像の画素濃度を、表示するための画像における同じ位置の微小領域Ｌに入れる。その作業を画像の全域にわたって繰り返すことにより、画像上のすべての微小領域Ｌで焦点のあった全焦点画像を、高速に取得することができる。
【０１６０】
また、全焦点画像を生成する別の処理方法としては、各微小領域Ｌにおいて求めた合焦点距離Ｄを回折距離として入力し、取得した干渉縞データを用いてフレネル回折計算を実施するようにしてもよい。そして、フレネル回折計算の結果得られた画像における微小領域Ｌの画素濃度を、表示するための画像における同じ位置の微小領域Ｌに入れる。その作業を画像の全域にわたって繰り返すことによって、画像上のすべての微小領域Ｌで、より正確に焦点のあった全焦点画像を取得することができる。
【０１６１】
以上説明した、この実施の形態の表面形状測定装置によれば、次のような作用効果を奏することがわかる。
【０１６２】
すなわち、図１に示す表面形状測定装置によれば、撮像した干渉縞データにより、被測定物１４の反射光のほぼ光軸方向における複数位置での複素振幅を求めて、そのうちの振幅データを用いて、合焦状態の検知により表面形状を求めるdepth from focus理論に基づいた演算処理を施すことによって、被測定物１４の表面形状データＳijを測定する。
【０１６３】
その場合、測定分解能および光軸方向の測定範囲は光学系の拡大倍率により任意に設定することができるため、その測定分解能を光源波長λの１／２より小さく設定することにより、レーザ干渉により測定できない光源波長λの１／２を超える段差を測定することができる。
【０１６４】
一方、干渉縞データから求めた複素振幅のうち、位相データを用いれば、通常のレーザ干渉と同様にしてサブフリンジオーダーの微小形状を測定することができる。したがって、被測定物１４の表面において光源波長λの１／２を超える段差については振幅データを用いた測定で受け持ち（表面形状データＳij）、それ以外の領域の形状測定については位相データを用いた測定で受け持ち（表面形状データＴij）、両測定結果を合成又は比較することにより（表面形状データＺij）、光源波長λの１／２を超える段差をもつ被測定面であっても、サブフリンジオーダーで測定可能となる。
【０１６５】
また、従来の走査型白色干渉計に比べて、光源波長λの１／２を超える段差を含む面を測定するために、干渉光学系や被測定物１４を機械的な駆動で移動させる必要がないため、測定誤差要因が減るうえ、移動範囲の制限によって光軸方向の測定範囲が制限されない。
振幅データを用いた測定については、レンズの焦点距離を変化させる振幅データを用いた従来法に対して、干渉縞データを用いて、光学系の光軸方向における複数位置での被測定物１４の反射光の複素振幅を計算し、そのうちの振幅データを用いてdepth from focus法を実施するため、数値演算でdepth from focus法を実施していることになり、depth from focus法を実施するために、やはり従来のように機械的な駆動を伴わないため、その誤差の影響を受けず、機械的駆動部分の可動範囲で決まる測定範囲の制限も生じない。
【０１６６】
また、測定のために機械的な駆動を要さないため、高速測定が可能であり、動きを伴う物体の動的形状が測定可能である。
【０１６７】
さらに、観測する干渉縞にキャリヤ周波数を与えて単一の干渉縞データから被測定物１４の反射光の複素振幅を計算するようにすると、そのデータが得られるまでに要する時間はＣＣＤ１０による撮像時間のみで決まることになり、前記した複数の画像を撮像する必要がある走査型白色干渉計やレンズ焦点距離を変化させる方法に比べて、高速測定が可能になる。
【０１６８】
このようにして、高速測定が可能になっても、被測定物１４の動きの速度がＣＣＤ１０による撮像時間より速くなると、被測定物１４からの反射光が安定しないため、正確に反射光の複素振幅を求めることができなくなり、正確な測定ができなくなる。そこで、前記のように、被測定物１４にパルス光を照射し、そのパルスの発光時間を被測定物１４の動きに対応して十分短くすることにより、被測定物１４が光を反射する瞬間において、被測定物１４をみかけ上静止した状態にして、被測定物１４からの反射光を安定させることができるため、反射光の複素振幅を正確に求めることができ、正確な測定が可能になる。
【０１６９】
変形、変位、あるいは振動中の物体の形状測定に関し、単一画像から物体反射光の複素振幅を求める方法として、パルス光源を用いたデジタルホログラフィによる物体の動的測定法（例えば、G.Pedrini,H.J.Tiziani,「Quantitative evaluation of two-dimensional dynamic deformations using digital holography」,Optics＆Laser Technology,Vol.29,No.5,pp249-256,1997を参照）があるが、かかる方法においては、粗面物体を測定対象としており、光源波長レベルでの物体形状の測定や、また光源波長λの１／２を超える段差を含む面を光源波長レベルで測定する方法については示されていない。
【０１７０】
かかる方法に比べ、この表面形状測定装置による測定は、被測定物１４が微小であっても、鏡面又はそれに近い面であっても測定可能で、被測定物１４の表面に光源波長λの１／２を超える段差が含まれていても測定可能である。
【０１７１】
また、この表面形状測定装置によれば、光学系倍率を設定することにより、振幅データを用いた表面形状測定と位相データを用いた表面形状測定のうち、振幅データを形状測定についてのみ、その測定分解能を設定することができる。
【０１７２】
一方、振幅データを用いた表面形状測定により光源波長λの１／２の段差を測定するには、その測定分解能は光源波長λの１／２より小さいほうが望ましく、通常の可視光を光源に用いると波長はｎｍオーダーになるため、干渉光学系の拡大倍率を上げたほうがよい。
【０１７３】
その場合、被測定物１４からの反射光又は参照光のどちらか一方を球面波にすることによって干渉光学系の拡大倍率を設定し、前記振幅データを用いた形状測定の測定分解能を光源波長の１／２より小さくすることができるが、被測定物１４からの反射光か参照光かのどちらか一方を球面波にすると、被測定物１４の反射光の波面と参照光の波面との曲率の差に応じたデフォーカス収差が位相データに含まれてくる。デフォーカス収差が含まれると、位相データから正確に被測定物１４の表面形状を求めることができなくなるため、デフォーカス収差を除去する必要がある。
【０１７４】
そのため、この表面形状測定装置では、被測定物１４からの反射光か参照光のどちらか一方を球面波にしたとき、他方の光についても球面波にし、両者の曲率をほぼ一致させることにより、上記のデフォーカス収差を除去することができる。これにより、干渉光学系の拡大倍率をあげた場合でも、測定を正確に実施することができる。
【０１７５】
この表面形状測定装置では、光学系倍率を設定することにより、振幅データを用いた表面形状測定と位相データを用いた表面形状測定のうち、振幅データを用いた表面形状測定についてのみ、その測定分解能を設定する。
【０１７６】
しかし、その測定分解能を過剰に小さくすると光学系倍率が大きくなるため、横方向の測定範囲（視野）は狭くなり、被測定物１４の表面における広い範囲（視野）を一度に測定することができなくなる。
【０１７７】
そこで、この表面形状測定装置では、光源波長の１／２より小さい測定分解能をもつという条件のもとで、被測定物１４の大きさに応じて光学系倍率を変化させるようにし、被測定物１４の表面のより広い範囲（視野）を一度に測定できるようにしている。
【０１７８】
合焦法の原理を用いて、合焦点はずれを検出することにより被測定物１４の表面形状を測定する場合、被測定物１４の表面が平滑で光学的にテクスチャを観測できないと、合焦点はずれの検知が困難になることがある。
【０１７９】
合焦点はずれは被測定物１４の表面のパターンのコントラストの低下により検出されるものである。２次元的な振幅データではコントラストが低下して合焦点はずれが検出できない場合でも、２次元的な位相データではコントラストが得られる場合がある。
【０１８０】
そのような場合は、前記のように位相データを２次元強度データとして用いて、合焦法の原理により被測定物１４の表面形状を測定する。これにより被測定物１４表面が平滑で、振幅データでは光学的にテクスチャを観測できない場合でも、被測定面の形状測定が可能となり、測定の汎用性を向上させることができる。
【０１８１】
物体形状を測定する場合、測定者が被測定物１４を観察しながら測定できたほうが、被測定物１４の設置、姿勢の調整などの点において便利である。被測定物１４が微小になると肉眼で観察することは難しくなるため、光学的に拡大して観察できたほうがよい。しかしながら、拡大光学系は焦点深度が浅く、光軸方向で焦点のあう範囲が狭いため、観察しづらいという欠点がある。
【０１８２】
これに対し、この表面形状測定装置では、振幅データを用いて合焦法の原理により被測定物１４の表面形状を測定するが、その場合、２次元画像内で合焦点の位置を求める処理を実施するため、測定のために取得した合焦点の情報を用いて、２次元的な全領域で焦点のあった全焦点画像を、生成、表示することができ、被測定物１４の位置、姿勢を観察しやすくすることにより、測定の操作性を向上させることができる。
【０１８３】
この場合、高速に画像を生成、表示することが重要な課題となる。この表面形状測定装置では、全焦点画像における微小領域の焦点画像を、合焦状態から求めた合焦距離に対して最も近い距離にある振幅データに置き換えて生成することによって、高速に全焦点画像を生成することができるので、より自然な被測定物１４の観測が可能となり、測定の操作性を向上させることができる。
【０１８４】
また、全焦点画像を生成する場合、全焦点画像における微小領域の焦点画像を、合焦状態から求めた合焦距離を用いてフレネル回折計算することにより、より正確な全焦点画像の生成が可能となり、測定における被測定物１４の位置、姿勢の調整が、より正確に実施可能となるので、測定の高速化、高精度化を図ることができる。
【０１８５】
【発明の効果】
請求項１に記載の発明は、被測定物の光源波長の１／２を超える段差については第１の表面形状データで受け持ち、それ以外の領域の形状測定については第２の表面形状データで受け持って、この両者を合成又は比較することにより、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【０１９０】
請求項２に記載の発明は、請求項１に記載の表面形状測定装置において、測定分解能を光源波長の１／２より小さく設定することにより、光源波長の１／２を超える段差も測定することができる。
【０１９１】
請求項３に記載の発明は、請求項１又は２に記載の表面形状測定装置において、反射光の光軸方向の複数位置は、それぞれの位置の間隔が光源の波長の１／２より小さくなるように設定しているので、光源波長の１／２を超える段差も測定することができる。
【０１９２】
請求項４に記載の発明は、請求項１〜３のいずれかの一に記載の表面形状測定装置において、第１と第２の表面形状データを的確に合成して、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【０１９３】
請求項５に記載の発明は、請求項１〜４のいずれかの一に記載の表面形状測定装置において、観測する干渉縞にキャリヤ周波数を与えて、単一の干渉縞データから被測定物の反射光の振幅データ、位相データを計算することにより、測定の高速化を図ることができる。
【０１９４】
請求項６に記載の発明は、請求項１〜５のいずれかの一に記載の表面形状測定装置において、パルスの発光時間を被測定物の動きに対応して十分短くすることで、被測定物が光を反射する瞬間に被測定物がみかけ上静止した状態となり、被測定物からの反射光を安定させて、正確に反射光の振幅データ、位相データを求めることができるので、被測定物の動きが高速な場合においても、正確に物体の動的形状を測定することができる。
【０１９５】
請求項７に記載の発明は、請求項１〜６のいずれかの一に記載の表面形状測定装置において、デフォーカス収差を除去して、干渉光学系の拡大倍率を高めた測定を正確に実施することができる。
【０１９６】
請求項８に記載の発明は、請求項１〜７のいずれかの一に記載の表面形状測定装置において、光源波長の１／２より小さい測定分解能のもとで、被測定物の大きさに応じて干渉光学系の拡大倍率を変化させるようにすれば、被測定物表面のより広い範囲を一度に測定することができる。
【０１９７】
請求項９に記載の発明は、請求項８に記載の表面形状測定装置において、被測定物表面のより広い範囲を一度に測定することを容易に行なうことができる。
【０１９８】
請求項１０に記載の発明は、被測定物の光源波長の１／２を超える段差については第１の表面形状データで受け持ち、それ以外の領域の形状測定については第２の表面形状データで受け持って、この両者を合成又は比較することにより、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【０２０３】
請求項１１に記載の発明は、被測定物の光源波長の１／２を超える段差については第１の表面形状データで受け持ち、それ以外の領域の形状測定については第２の表面形状データで受け持って、この両者を合成又は比較することにより、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【０２０８】
請求項１２記載の発明は、請求項１１に記載の発明と同様の作用、効果を奏する。
【図面の簡単な説明】
【図１】この発明の一実施の形態である表面形状測定装置の概略構成を示すブロック図である。
【図２】前記表面形状測定装置における光学系の位置関係を示す説明図である。
【図３】同説明図である。
【図４】前記表面形状測定装置が行なう処理で再生距離と合焦測度との関係を示すグラフである。
【図５】デフォーカス収差について説明する説明図である。
【図６】前記表面形状測定装置を構成するコンピュータの電気的な接続を示すブロック図である。
【図７】コンピュータが実行する処理のフローチャートである。
【図８】前記表面形状測定装置の他の例について概略構成を示すブロック図である。
【図９】同ブロック図である。
【符号の説明】
１ 光源
５〜９，１１〜１３ 干渉光学系
１０ 撮像素子
１４ 被測定物
２４，２５，２８ 拡大倍率可変装置
３７ 記憶媒体
【発明の属する技術分野】
この発明は、被測定物の表面形状を測定する表面形状測定装置及びその方法、プログラム並びに記憶媒体に関する。
【０００２】
【従来の技術】
物体の表面形状を光の波長レベルで測定可能な技術として、レーザ干渉計（例えば、Ｚｙｇｏ社のレーザ干渉システム等）がある。このレーザ干渉計は、参照光の位相に変調をかけながら、ＣＣＤなどの撮像素子を用いて複数の干渉縞画像を収録して、それらの干渉縞から被測定物の反射光の位相を求め、その位相データに基づいて光の波長レベルで物体形状を求めるものである。
【０００３】
しかしながら、物体が光源波長の１／２を超えるような段差を含む場合は、その段差の位置において、位相のアンラップができなくなり、これにより位相の関連性がとぎれるため、形状測定値に不確定性が含まれてしまうという不具合がある。
【０００４】
走査型白色干渉計（例えば、Ｚｙｇｏ社の三次元表面構造解析顕微鏡等）は、かかる不具合を解消するもので、複数の波長でそれぞれ求めた位相情報の相互関係から、単一波長での光軸方向の測定範囲λ／２を拡大する方法を用いている（多波長測定）。そして段差がさらに大きくなった場合には、参照光と物体反射光の光路長がほぼ一致したときにのみ白色干渉縞が発生する性質を利用して、干渉計、被測定物体のどちらかを干渉計の光軸方向に走査し、白色干渉縞の位置と、干渉計又は被測定物の移動量を検知することで、絶対測長を実施している。
【０００５】
しかし、このような走査型白色干渉計により大きな段差を測定する場合には、測定のために被測定物又は測定光学系を測定光学系の光軸方向に移動させる必要があるために、測定時間がかかるという不具合がある。また、測定中の移動であり、かつ、干渉計又は被測定物の移動量を検出する必要があるため、測定誤差要因が増え、また移動部の移動範囲内に光軸方向測定範囲が制限されるという不具合がある。
【０００６】
これに対し、特開2000-221013に開示の技術は、波長を超える段差や絶対距離を含んだ形状の測定を高精度かつ高速に測定ミスすることなく行なうため、走査型白色干渉計において、取得すべき干渉縞画像を低減し、測定の高速化を図ろうとしたものである。
【０００７】
また、特開2001-41724に開示の技術では、干渉計を用いた表面形状測定装置において、測定対象物を移動させることなく、測定対象物表面の三次元形状を高精度に測定できるようにするため、複屈折位相板を利用した位相シフト法を実施しており、測定対象物を移動させることなく三次元形状の測定を可能としている。また、位相板の代わりに電気光学素子を利用して電気的に位相シフトをさせれば、高速測定が可能としている。
【０００８】
さらに、物体をＣＣＤなどで観測した強度データをもとに、物体形状を求める技術として、depth from focus理論に基づく方法がある。これはレンズの焦点距離を変化させながら収録した複数の撮像画像をもとに、物体表面の合焦位置を求め、それにより表面形状を求めるものである（例えば、石原満宏、佐々木博美「合焦法による高速三次元形状計測」精密工学会誌，vol.63,No.1,1997を参照）。
【０００９】
【発明が解決しようとする課題】
しかしながら、特開2000-221013に開示の技術では、サブフリンジオーダーで形状を求めるのに、位相変調法（位相シフト法）を用いているため、装置内に機械的な駆動部分を有しており、また、複数の画像を取得する必要があるため、走査型白色干渉計よりは速いとはいえ、高速化に限界がある。そのため、動きを伴う物体の動的形状を測定することは不可能である。また、機械的な駆動を伴う測定であるために、測定誤差要因が増えるという不具合は依然として残る。
また、特開2001-41724に開示の技術では、原理的には位相シフト法を用いたレーザ干渉測定であり、上述したような光源波長の１／２を超える段差を含む面の測定には不確定性がある。また、高速化したといっても、複数の画像データを必要とするために高速化に限界があり、動的形状の測定は不可能である。
【００１０】
さらに、前記のdepth from focus理論に基づく技術では、一般にレーザ干渉に比べると測定分解能が低いため、サブフリンジオーダーの形状が要求される部品の評価には適さないという不具合がある。また、光学系の倍率を上げることで、測定分解能をある程度向上させることができるが、レーザ干渉による分解能と同じレベルまで上げようとすると、横方向の測定範囲（視野）が狭くなってしまうという不具合がある。またレンズの焦点距離を変えながら複数の撮像画像を収録する方法では、収録方法を工夫したとしても複数の画像を取り込むだけの時間がかかるうえ、レンズの焦点距離を正確に変化させるための構成が必要となり、装置構成が複雑化する。また、移動部分が存在するために誤差要因も増え、光軸方向の測定範囲も移動部の可動範囲に限定されるといった不具合がある。
【００１１】
この発明の目的は、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することである。
【００１２】
この発明の目的は、２次元的な振幅データではコントラストが低下して合焦点はずれが検知できない場合にも、被測定物の表面が平滑で、振幅データでは光学的にテクスチャを観測できない場合に、被測定面の形状測定ができるようにすることである。
【００１３】
この発明の目的は、被測定物の位置、姿勢を観察しやすくし、測定の操作性を向上させることである。
【００１４】
この発明の目的は、高速に全焦点画像を生成でき、より自然な被測定物の画像の観測が可能となるようにすることである。
【００１５】
この発明の目的は、より正確な全焦点画像の生成ができ、測定における被測定物の位置、姿勢の調整がより正確に実施可能となり、また、測定の高速化、高精度化を図ることができるようにすることである。
【００１６】
この発明の目的は、測定の高速化を図れるようにすることである。
【００１７】
この発明の目的は、被測定物の動きが高速な場合においても、正確に物体の動的形状を測定することができるようにすることである。
【００１８】
この発明の目的は、デフォーカス収差を除去して、干渉光学系の拡大倍率を高めた測定を正確に実施することである。
【００１９】
この発明の目的は、被測定物表面のより広い範囲を一度に測定できるようにすることである。
【００２０】
この発明の目的は、被測定物表面のより広い範囲を一度に測定することを容易に行なえるようにすることである。
【００２１】
【課題を解決するための手段】
請求項１に記載の発明は、被測定物に照射する光を発する光源と、その光の被測定物における反射光と所定の参照光との間で干渉縞を発生させる干渉光学系と、この干渉縞を撮像して当該干渉縞の画像データである干渉縞データを出力する撮像素子と、前記干渉縞データに基づいて前記反射光の光軸方向の複数位置における前記被測定物の反射光の位相を示す位相データ及び振幅を示す振幅データを求める光学データ検出手段と、前記振幅データの微分値の和である合焦測度を算出し、前記合焦測度が最大となる合焦距離に基づき、前記被測定物の第１の表面形状データを求める第１の表面形状データ作成手段と、前記位相データを長さ単位に変換することで、前記被測定物の第２の表面形状データを求める第２の表面形状データ作成手段と、前記第１の表面形状データと前記第２の表面形状データとを合成又は比較して、前記被測定物の表面形状のデータである第３の表面形状データを求める合成又は比較手段と、を備えている表面形状測定装置である。
【００２２】
したがって、被測定物の光源波長の１／２を超える段差については第１の表面形状データで受け持ち、それ以外の領域の形状測定については第２の表面形状データで受け持って、この両者を合成又は比較することにより、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【００３１】
請求項２に記載の発明は、請求項１に記載の表面形状測定装置において、前記干渉光学系は、その拡大倍率の２乗に反比例する感度で得られる前記振幅データによる前記被測定物の表面形状における前記反射光の光軸方向の測定分解能を、前記光源の波長の１／２より小さくなるように所定の拡大倍率に設定されている。
【００３２】
したがって、測定分解能を光源波長の１／２より小さく設定することにより、光源波長の１／２を超える段差も測定することができる。
【００３３】
請求項３に記載の発明は、請求項１又は２に記載の表面形状測定装置において、前記データ検出手段は、前記反射光の光軸方向の複数位置を、それぞれの位置の間隔が前記光源の波長の１／２より小さいものとしている。
【００３４】
したがって、反射光の光軸方向の複数位置は、それぞれの位置の間隔が光源の波長の１／２より小さくなるように設定しているので、光源波長の１／２を超える段差も測定することができる。
【００３５】
請求項４に記載の発明は、請求項１〜３のいずれかの一に記載の表面形状測定装置において、前記合成又は比較手段は、前記合成として、
ｍ×（λ／２）＋δ
（ｍ：次数，λ：光源波長，δ：前記第２の表面形状データ）における次数ｍを決定する処理を行なうものである。
【００３６】
したがって、第１と第２の表面形状データを的確に合成して、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【００３７】
請求項５に記載の発明は、請求項１〜４のいずれかの一に記載の表面形状測定装置において、前記光学データ検出手段は、前記位相データ及び振幅データをキャリヤ周波数を与えた単一の前記干渉縞の干渉縞データから求める。
【００３８】
したがって、観測する干渉縞にキャリヤ周波数を与えて、単一の干渉縞データから被測定物の反射光の振幅データ、位相データを計算することにより、測定の高速化を図ることができる。
【００３９】
請求項６に記載の発明は、請求項１〜５のいずれかの一に記載の表面形状測定装置において、前記光源は、前記被測定物にパルス光を照射するものである。
【００４０】
したがって、パルスの発光時間を被測定物の動きに対応して十分短くすることで、被測定物が光を反射する瞬間に被測定物がみかけ上静止した状態となり、被測定物からの反射光を安定させて、正確に反射光の振幅データ、位相データを求めることができるので、被測定物の動きが高速な場合においても、正確に物体の動的形状を測定することができる。
【００４１】
請求項７に記載の発明は、請求項１〜６のいずれかの一に記載の表面形状測定装置において、前記変換光学系は、前記反射光及び前記参照光をそれぞれ球面波に変換し、前記撮像素子による干渉縞の撮像位置での前記両球面波の曲率をほぼ一致させて前記干渉光学系の倍率を設定している。
【００４２】
したがって、デフォーカス収差を除去して、干渉光学系の拡大倍率を高めた測定を正確に実施することができる。
【００４３】
請求項８に記載の発明は、請求項１〜７のいずれかの一に記載の表面形状測定装置において、前記干渉光学系の拡大倍率を可変する拡大倍率可変装置を備え、前記第１の表面形状データ作成手段及び前記第２の表面形状データ作成手段のうち前者による前記被測定物の測定分解能又は前記反射光の光軸方向の測定範囲のみを設定できるようにしている。
【００４４】
したがって、光源波長の１／２より小さい測定分解能のもとで、被測定物の大きさに応じて干渉光学系の拡大倍率を変化させるようにすれば、被測定物表面のより広い範囲を一度に測定することができる。
【００４５】
請求項９に記載の発明は、請求項８に記載の表面形状測定装置において、前記拡大倍率可変装置は、前記撮像素子による前記干渉縞の撮像位置における前記反射光の曲率と前記参照光の曲率とがほぼ一致するよう連動させて曲率を変化させることにより、前記干渉光学系の倍率を可変としている。
【００４６】
したがって、被測定物表面のより広い範囲を一度に測定することを容易に行なうことができる。
【００４７】
請求項１０に記載の発明は、被測定物に光を照射したときの反射光と所定の参照光との間で発生させた干渉縞を撮像して当該干渉縞の画像データである干渉縞データを取得する撮像工程と、前記干渉縞データに基づいて前記反射光の光軸方向の複数位置における前記被測定物の反射光の位相を示す位相データ及び振幅を示す振幅データを求める光学データ検出工程と、前記振幅データの微分値の和である合焦測度を算出し、前記合焦測度が最大となる合焦距離に基づき、前記被測定物の第１の表面形状データを求める第１の表面形状データ作成工程と、前記位相データを長さ単位に変換することで、前記被測定物の第２の表面形状データを求める第２の表面形状データ作成工程と、前記第１の表面形状データと前記第２の表面形状データとを合成又は比較して、前記被測定物の表面形状のデータである第３の表面形状データを求める合成又は比較工程と、を含んでなる表面形状測定方法である。
【００４８】
したがって、被測定物の光源波長の１／２を超える段差については第１の表面形状データで受け持ち、それ以外の領域の形状測定については第２の表面形状データで受け持って、この両者を合成又は比較することにより、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【００５７】
なお、前記データ検出工程は、前記反射光の光軸方向の複数位置を、それぞれの位置の間隔が前記光源の波長の１／２より小さいものとすることができる。
【００５８】
これにより、反射光の光軸方向の複数位置は、それぞれの位置の間隔が光源の波長の１／２より小さくなるように設定しているので、光源波長の１／２を超える段差も測定することができる。
【００５９】
前記合成又は比較工程は、前記合成として、
ｍ×（λ／２）＋δ
（ｍ：次数，λ：光源波長，δ：前記第２の表面形状データ）
における次数ｍを決定する処理を行なうものである。
【００６０】
これにより、第１と第２の表面形状データを的確に合成して、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【００６１】
前記光学データ検出工程は、前記位相データ及び振幅データをキャリヤ周波数を与えた単一の前記干渉縞の干渉縞データから求めることができる。
【００６２】
これにより、観測する干渉縞にキャリヤ周波数を与えて、単一の干渉縞データから被測定物の反射光の振幅データ、位相データを計算することにより、測定の高速化を図ることができる。
【００６３】
前記撮像工程は、前記被測定物に照射する光をパルス光とすることができる。
【００６４】
これにより、パルスの発光時間を被測定物の動きに対応して十分短くすることで、被測定物が光を反射する瞬間に被測定物がみかけ上静止した状態となり、被測定物からの反射光を安定させて、正確に反射光の振幅データ、位相データを求めることができるので、被測定物の動きが高速な場合においても、正確に物体の動的形状を測定することができる。
【００６５】
請求項１１に記載の発明は、被測定物に光を照射したときの反射光と所定の参照光との間で発生させた干渉縞を撮像した画像データである干渉縞データに基づいて前記反射光の光軸方向の複数位置における前記被測定物の反射光の位相を示す位相データ及び振幅を示す振幅データを求める光学データ検出処理と、前記振幅データの微分値の和である合焦測度を算出し、前記合焦測度が最大となる合焦距離に基づき、前記被測定物の第１の表面形状データを求める第１の表面形状データ作成処理と、前記位相データを長さ単位に変換することで、前記被測定物の第２の表面形状データを求める第２の表面形状データ作成処理と、前記第１の表面形状データと前記第２の表面形状データとを合成又は比較して、前記被測定物の表面形状のデータである第３の表面形状データを求める合成又は比較処理と、をコンピュータに実行させるコンピュータに読み取り可能なプログラムである。
【００６６】
したがって、被測定物の光源波長の１／２を超える段差については第１の表面形状データで受け持ち、それ以外の領域の形状測定については第２の表面形状データで受け持って、この両者を合成又は比較することにより、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【００７５】
なお、前記データ検出処理は、前記反射光の光軸方向の複数位置を、それぞれの位置の間隔が前記光源の波長の１／２より小さいものとすることができる。
【００７６】
これにより、反射光の光軸方向の複数位置は、それぞれの位置の間隔が光源の波長の１／２より小さくなるように設定しているので、光源波長の１／２を超える段差も測定することができる。
【００７７】
前記合成又は比較処理は、前記合成として、
ｍ×（λ／２）＋δ
（ｍ：次数，λ：光源波長，δ：前記第２の表面形状データ）
における次数ｍを決定する処理を行なうものである。
【００７８】
これにより、第１と第２の表面形状データを的確に合成して、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【００７９】
前記光学データ検出処理は、前記位相データ及び振幅データをキャリヤ周波数を与えた単一の前記干渉縞の干渉縞データから求めることができる。
【００８０】
これにより、観測する干渉縞にキャリヤ周波数を与えて、単一の干渉縞データから被測定物の反射光の振幅データ、位相データを計算することにより、測定の高速化を図ることができる。
【００８１】
前記撮像処理は、前記被測定物に照射する光をパルス光とすることができる。
【００８２】
これにより、パルスの発光時間を被測定物の動きに対応して十分短くすることで、被測定物が光を反射する瞬間に被測定物がみかけ上静止した状態となり、被測定物からの反射光を安定させて、正確に反射光の振幅データ、位相データを求めることができるので、被測定物の動きが高速な場合においても、正確に物体の動的形状を測定することができる。
【００８３】
請求項１２に記載の発明は、請求項１１に記載のプログラムを記憶した記憶媒体である。
【００８４】
したがって、請求項１１に記載の発明と同様の作用、効果を奏する。
【００８５】
【発明の実施の形態】
この発明の一実施の形態について説明する。
【００８６】
図１は、この発明の一実施の形態である表面形状測定装置の概略構成を示すブロック図である。図１に示すように、この表面形状測定装置は、光源となるＨｅ−Ｎｅレーザ１と、その出力するレーザ光の被測定物１４への照射光強度を調整するためのＮＤフィルタ２と、レーザ光を拡大するためのビームエキスパンダ３とを備えている。ビームエキスパンダ３はスペイシャルフィルタとしても作用する。
【００８７】
ビームエキスパンダ３にて拡大された光は、ミラー４にて折り返され、ビームスプリッタ５に入射する。ビームスプリッタ５で反射した光は、光を球面波に変換するためのレンズ６を通り、ミラー７で折り返され、光強度を調整するためのＮＤフィルタ８、ハーフミラー９を介して撮像素子であるＣＣＤ１０に入射する。この光は、被測定物１４の反射光と干渉するための参照光となる。
【００８８】
一方、ビームスプリッタ５を透過した光は、ミラー１１にて折り返され、レンズ１２により球面波に変換され、ハーフミラー９を透過して、対物レンズ１３によりほぼ平行光に変換されて、被測定物１４に照射する。被測定物１４で反射した光は、対物レンズ１３を通過し、ハーフミラー９で反射して物体光としてＣＣＤ１０に到達し、前記の参照光との間に干渉をおこして干渉縞を発生する。したがって、これらの光学素子により干渉光学系を実現している。この干渉縞はＣＣＤ１０にて撮像され、フレームグラバー１５を介してコンピュータ１６に転送され、コンピュータ１６に記録される。
【００８９】
ビームスプリッタ５で分岐された後の参照光路と物体光路の長さは、光源であるＨｅ−Ｎｅレーザ１のコヒーレンス長以下になるように設定されている。また、被測定物１４からの反射光がＣＣＤ１０で撮像するのに適した強度になるように、ＮＤフィルタ２により被測定物１４への照射光強度を調整してあり、さらに、被測定物１４からの反射光と参照光との干渉により発生する干渉縞のコントラストが高くなるように、ＮＤフィルタ８により参照光強度が調整されている。
【００９０】
対物レンズ１３は被測定物１４の反射光を球面波に変換するように作用し、これにより後述のように被測定物１４の反射光の複素振幅（被測定物１４の像）を再生したときに、像が拡大して観測される。レンズ１２と対物レンズ１３の位置関係は、被測定物１４にほぼ平行光が照射されるように調整されている。レンズ６は、参照光を球面波に変換するように作用し、対物レンズ１３により球面波に変換された被測定物１４の反射光の曲率と、参照光の曲率とが、ＣＣＤ１０の撮像位置においてほぼ一致するように、位置が調整されている。
【００９１】
被測定物１４の像を拡大しない場合は、対物レンズ１３は不要であり、それに伴い、レンズ６、レンズ１２も不要である。また後で詳細を説明するが、ここでは干渉縞をホログラムとして扱い、ＣＣＤ１０にて撮像したホログラム干渉縞に仮想的に参照光を照射して、その回折光波により被測定物１４からの反射光（被測定物１４の像）の複素振幅を再生する。単一の干渉縞データをＣＣＤ１０にて撮像し、被測定物１４の反射光を再生する場合、ホログラムを透過する０次回折光と実像と虚像の３つが再生されるため、この３つを分離して再生できるように、被測定物１４からの反射光の光軸と参照光の光軸との間に適当な傾きを与えて、干渉縞にキャリヤ周波数がのるようにしている。前記の両光軸に傾きを与えるためには、例えば、ミラー７、被測定物１４又はハーフミラー９かの傾きを調整すればよい。
【００９２】
次に、被測定物１４からの反射光（被測定物１４の像）の複素振幅（この明細書において、光の振幅と位相を併せて複素振幅と呼ぶ）を再生計算する方法について説明する。
【００９３】
まず、被測定物１４からの反射光と参照光との干渉により発生する干渉縞のＣＣＤ１０の撮像位置での強度は、（１）式で与えられる。
【００９４】
【数１】
ここで、Ｉは干渉縞の強度であり、ｘ，ｙはＣＣＤ１０の撮像位置のｘ，ｙ座標、Ｒは参照光の複素振幅であり、Ａ，φはそれぞれ被測定物１４の反射光の振幅と位相である。
【００９６】
干渉縞をＣＣＤ１０にて撮像し、ホログラム画像データとしてコンピュータ１６に記憶する。乾板を用いたホログラフィでは、干渉縞を記録した乾板（ホログラム）に改めて参照光を照射すると、参照光がホログラムに記録された干渉縞にて回折し、それが記録時の物体反射光として振舞うので、これにより、物体反射光（物体像）が再生されることになった。ここでは、実際に参照光を照射せずに、仮想的に参照光がホログラムに照射されたものとして、ＣＣＤ１０にて記録したホログラム画像からフレネル近似のもとに物体反射光（物体像）を再生する。ホログラムに参照光として平行光が照射されたとすると、ホログラムによる回折光、すなわち物体反射光は（２）式のように表される。
【００９７】
【数２】
ここで、Ｕは再生距離ｄだけ離れた像面位置での物体反射光波の複素振幅、ｘ’，ｙ’は像面位置でのｘ，ｙ座標、ｃは複素定数、λは光源波長を表す。
【００９９】
したがって、再生距離ｄを入力し、収録した干渉縞データを用いて（２）式を計算することにより、再生距離ｄにおける物体反射光の複素振幅が再生されることになる。被測定物１４の像を拡大しない場合は、再生距離ｄにＣＣＤ１０の撮像位置と被測定物１４の表面との距離を代入すればよい。図１に示した装置のように、被測定物１４の像を拡大する場合は、対物レンズ１３による被測定物１４の拡大像面とＣＣＤ１０の撮像位置との距離を、再生距離ｄとして入力する。
【０１００】
図２は、図１における被測定物１４、対物レンズ１３、ホログラム面（ＣＣＤ１０の撮像位置）、対物レンズ１３による被測定物１４の拡大像それぞれの位置関係を示す説明図である。図２において、符号１７は被測定物１４の位置を示し、符号１８は対物レンズ１３の位置、符号１９はホログラム面（ＣＣＤ１０の撮像位置）、符号２０は対物レンズ１３による被測定物１４の拡大像面である。
【０１０１】
図２における距離ｄ’が再生距離となり、それを（２）式に入力して、被測定物１４の像を再生計算する。図２の距離ｓ及び距離ｓ’はレンズ結像の下記（３）式の関係にあるため、対物レンズ１８の焦点距離ｆと距離ｓをあらかじめ測定しておくことにより、距離ｓ’を求めることができる。そして、対物レンズ１３（位置１８）からほぼ距離ｓ’の位置にＣＣＤ１０の撮像位置を設置すれば、距離ｄ’はほぼ０となり、目視ではほぼ焦点のあった被測定物１４の像が得られ、被測定物１４の位置などの調整が容易になる。また、距離ｄ’が大きくなると、再生した像の周辺部に歪が生じる場合があるため、歪の少ない再生した像を得るためにも効果がある。
【０１０２】
【数３】
再生像には、ホログラムを透過する光（０次回折光）による像、実像、虚像の３つの像が含まれていて、被測定物１４の形状を求めるためには、そのうちの実像に注目する。撮像時において物体反射光と参照光との光軸が平行に近いと、その３つの像が重なって再生されるため、前記したように、物体反射光と参照光の光軸に相対的な傾きを与えて、干渉縞にキャリヤ周波数をのせた状態でホログラム（干渉縞画像）を取り込み、これにより像面位置でそれらの３つの像を離れた位置に再生することができる。
【０１０４】
（２）式で表される複素振幅の実部を“Real｛Ｕ（ｘ’，ｙ’）｝”とし、虚部を“Imaginary｛Ｕ（ｘ’，ｙ’）｝”とすると、被測定物１４の反射光の振幅Ａ、位相φは、それぞれ次の（４）、（５）式にて与えられる。
【０１０５】
【数４】
（４）式にて計算される振幅Ａと、（５）式にて計算される位相φが、それぞれ、振幅データ、位相データに相当する。この位相データの示す位相φに、“λ／（４π）”を乗じて長さ単位に変換し、位相の飛びがある場合、それをアンラップすることにより、被測定物１４の物体形状を光源の波長λ以下の精度で求めることができる。
【０１０７】
しかしながら、被測定物１４の物体面に、レーザ光の波長λに対し、λ／２を超えるような大きな段差があると、その段差位置で位相φ（形状）の関連性が途切れるため、測定値に不確定性が生じてくる。
【０１０８】
そのため、この形状測定装置では、（４）式による振幅データを活用することによっても、被測定物１４の表面形状を測定し、振幅データによる表面形状の測定値と位相データによる表面形状の測定値とを比較し、又は、合成することによって、λ／２を超える段差を含む被測定面であっても、表面形状の測定を可能なものとしている。かかる点について以下に説明する。
【０１０９】
まず、振幅データを用いた表面形状の測定について説明する。簡単のため、被測定物１４を、光学系を用いて拡大しない場合における被測定物１４、ホログラム面（ＣＣＤ１０の撮像位置）、及び、再生像の位置関係を表す図３を参照して説明する。図３において、符号２１は被測定物１４の位置、符号２２はホログラム面（ＣＣＤ１０の撮像位置）、符号２３は被測定物１４の再生像を示す。
【０１１０】
ホログラム（ＣＣＤ１０の撮像位置）を原点とした場合、再生距離ｄを−ｘ１とすると、ｘ１の位置にある面の実像が焦点のあった状態で再生され、同様に、再生距離ｄを−ｘ２，−ｘ３とすると、それぞれｘ２，ｘ３の位置にある面の実像が焦点のあった状態で再生されることを示している。その場合、“ｄ＝−ｘ１”とした場合には、ｘ２，ｘ３にある面は焦点ずれすることになる。したがって焦点があっているかどうかを検出することにより、ｘ１，ｘ２，ｘ３の相対的な距離の差（すなわち段差）を求めることができる。
【０１１１】
その焦点が合っているかどうかの判定のために、ここでは、（４）式で求めた振幅データを用いて、depth from focus理論に基づいた演算を実施する。depth from focus理論に基づく演算方法としてはさまざまな種類のものがあるが、例えば、“石原満宏、佐々木博美「合焦法による高速三次元形状計測」精密工学会誌，vol.63,No.1,1997”に示された方法を用いることができる。
【０１１２】
すなわち、撮像した干渉縞データを用いて、再生距離ｄを微小量（Δｄ）ずつ変えながら、複数の複素振幅を（２）式を用いて計算し、各々の再生距離ｄにおける振幅データを（４）式により計算する。そして各々の振幅データにおいて、ＣＣＤ１０の撮像面の任意の小領域Ｌにおける振幅値Ａの微分値の和を合焦測度ｖとして求め、これらの合焦測度ｖのうちピークを示す再生距離ｄを、焦点のあった位置と判定する。
【０１１３】
図４に、撮像画像の任意の小領域Ｌにおいて、再生距離ｄを変化させたときの合焦測度ｖの変化を示すが、図４の距離Ｄが焦点のあう再生距離ｄ（合焦距離Ｄ）となる。したがって、任意の小領域Ｌにおいて振幅値Ａの微分値の和がピークになる再生距離ｄを求め、それを画像全体にわたって繰り返すことにより、被測定物１４の表面全体の各小領域Ｌにおける合焦距離Ｄが計算できる。そして基準位置を任意の小領域にとり、基準位置の再生距離との差を求めることによって、被測定物１４全体の表面形状を求めることができる。よって、撮像画像の各画素について、被測定物１４の表面形状データＳijが得られる。ここで、ｉ，ｊは、撮像画像の画素単位で表わされる座標を示す。
【０１１４】
合焦測度ｖを求めるとき、再生距離ｄを少しずつ変化させながら、被測定物１４の反射光の略光軸方向における複数位置での被測定物１４の反射光の超えるを計算により求めるが、その場合、再生距離ｄの変化の間隔（ピッチ）Δｄは、光源波長λの１／２より小さいことが望ましい。ただし、図４に示したように、複数の合焦測度ｖのプロットをガウス関数などの関数に近似し、関数の頂点を求める作業により合焦距離Ｄを求めてもよく、その場合は、前記の再生距離ｄの変化の間隔Δｄは、厳密に光源波長λの１／２より小さくなくてもよい。
【０１１５】
ここで、振幅データを用いた測定の縦方向の感度γは、定数ｋと光学系の倍率βを用いて、
【０１１６】
【数５】
と表わすことができ、光学系の倍率の２乗に反比例して分解能が上がる。
【０１１８】
一方、光源波長λの１／２の段差を、振幅データを用いる方法により測定しようとすると、光源波長λの１／２より小さい測定分解能が、振幅データを用いる方法に要求される。
【０１１９】
光の波長は、可視光の場合は数百ｎｍであるため、振幅データを用いた測定の測定分解能はｎｍオーダーでなければならない。それを考慮すると、被測定物１４の像を拡大せずに振幅データを用いた測定分解能の要求を満たそうとするのは困難であると考えられる。そこで、図１の装置構成では、対物レンズ１３を用いて被測定物１４の像を拡大して測定する装置構成を示した。ただし、光源であるＨｅ−Ｎｅレーザ１の波長をμｍやｍｍのオーダーまで大きくした場合などにおいては、被測定物１４の像を拡大せずとも、振幅データを用いた測定の測定分解能の要求は満足でき、測定が可能となる。
次に、（５）式により計算した位相データを用いて、被測定物１４の表面におけるλ／２を下回る微小形状を演算する方法について説明する。以下では、図１における、レンズ６、レンズ１２、対物レンズ１３を用いずに被測定物１４の像を拡大しない場合について説明する。
【０１２０】
すなわち、任意の小領域において、振幅データを用いる方法で焦点のあった再生距離ｄで再生した複素振幅を用い、（４）式により位相データを計算する。そして計算した位相φに“λ／４π”を乗じて長さ単位に変換して、被測定物１４の表面形状を求める。この場合も撮像画像の各画素について、被測定物１４の表面形状データが得られる。
図１に示した装置のように、対物レンズ１３を用いて被測定物１４の像を拡大した場合において、被測定物１４の反射光の曲率と参照光の曲率が一致していないと、前記の曲率の差に起因して、位相φの測定値に図５に示すような同心円状のデフォーカス収差が含まれてくる。被測定物１４の表面形状を正確に求めるためには、それを除去しなければならない。
【０１２１】
図１においてレンズ６の位置を調整して、被測定物１４の反射光の曲率と参照光の曲率をほぼ一致させると、両者の曲率の差はなくなるため、それにより位相データにデフォーカス収差は含まれなくなる。したがって任意の小領域Ｌにおいて、振幅データを用いる方法で焦点のあった再生距離ｄで再生した複素振幅を用い、（５）式により位相データを計算する。そして、この計算した位相φに“λ／４π”を乗じて長さ単位に変換して、被測定物１４の表面形状を求める。この場合も、撮像画像の各画素について被測定物１４の表面形状データＴijが得られる。ここで、ｉ，ｊは画像の画素単位で表わされる座標を示す。
【０１２２】
次に、振幅データにより求めた形状データＳijと、位相データにより求めた形状データＴijとを合成する方法について説明する。まず、振幅データより得られた形状データＳijの値をλ／２で除算し、その商ｍを求める。そして、形状データＳijの値にｍとλ／２の積を加えたものを、形状データＺijの値に代入する。この処理を、撮像画像の各画素にわたって実施することにより、被測定物１４のλ／２を超える段差を含む表面であっても、表面形状の正しい測定値を得ることができる。
【０１２３】
以上の処理は、言いかえると、次のように表わすこともできる。すなわち、光源の波長λの１／２を超える段差を含む面の正しい形状データをＺijとし、光源波長をλ、次数をｍij、位相データにより求めた形状をδijとすると、次の（７）式に示す関係が成り立つので、撮像画像の各画素において次数ｍijを求めればよい。
【０１２４】
【数６】
以上のようにして測定を行なうことにより、被測定物１４や光学系を測定中に移動させることなく、光源波長λの１／２を超えるような段差を含む面をサブフリンジオーダーで高精度に測定可能となる。
【０１２６】
次に、以上説明した方法を用いて、コンピュータ１６が行なう形状データＺijの測定処理について概要を整理して説明する。
【０１２７】
まず、コンピュータ１６の構成について説明する。図６は、コンピュータ１６の電気的な接続を示すブロック図である。図６に示すように、コンピュータ１６は、各種演算を行ないコンピュータ１６の各部を集中的に制御するＣＰＵ３１と、各種のＲＯＭ、ＲＡＭからなるメモリ３２とが、バス３３で接続されている。
【０１２８】
バス３３には、所定のインターフェイスを介して、ハードディスクなどの磁気記憶装置３４と、マウス、キーボード等により構成される入力装置３５と、ＬＣＤ、ＣＲＴ等の表示装置３６と、光ディスクなどの記憶媒体３７を読み取る記憶媒体読取装置３８とが接続され、また、インターネットなどのネットワーク３９と通信を行なう所定の通信インターフェイス４０が接続されている。なお、記憶媒体３７は、この発明の記憶媒体を実施するものであり、ＣＤ，ＤＶＤなどの光ディスク、光磁気ディスク、フロッピーディスクなどの各種方式のメディアを用いることができる。また、記憶媒体読取装置３８は、具体的には記憶媒体３７の種類に応じて光ディスク装置、光磁気ディスク装置、フロッピーディスク装置などが用いられる。
【０１２９】
磁気記憶装置３４には、この発明のプログラムを実現する測定演算プログラムが記憶されている。この測定演算プログラムは、記憶媒体３７に記憶されていたものを記憶媒体読取装置３８により読み取るか、あるいは、インターネットなどのネットワーク３９からダウンロードするなどして、磁気記憶装置３４にインストールしたものである。このインストールによりコンピュータ１６は、形状データＺijの測定処理の実行が可能な状態となる。この測定演算プログラムは、所定のＯＳ上で動作するものであってもよい。
【０１３０】
以下では、測定演算プログラムに基づいてコンピュータ１６が行なう形状データＺijの測定処理の内容について、図７のフローチャートを参照して説明する。
【０１３１】
図７に示すように、まず、ＣＰＵ３１は、ＣＣＤ１０にて撮像された干渉縞の干渉縞データを取り込んで、メモリ３２に記憶し（ステップＳ１）、メモリ３２の所定領域に記憶される所定のカウント値ｋを０にリセットし（ステップＳ２）、干渉縞データを用いて以下の処理を行なう。ステップＳ１により、撮像工程を実現している。
【０１３２】
まず、再生距離ｄの初期値ｄ０に対し、再生距離ｄの微小量Δｄにカウント値ｋを乗算した値を加算し（ステップＳ３）、フレネル回折計算を行って、“再生距離ｄ０＋Δｄ×ｋ”での複素振幅を取得する（ステップＳ４）。ステップＳ３，Ｓ４により光学データ検出手段、光学データ検出工程、光学データ検出処理を実現している。そして、その複素振幅の振幅データにおける微小領域Ｌでの合焦測度を求め（ステップＳ５）、カウント値ｋを＋１だけインクリメントする（ステップＳ６）、という処理を、カウント値ｋが所定値Ｋになるまで繰り返す（ステップＳ７のＮ）。
【０１３３】
カウント値ｋが所定値Ｋになったときは（ステップＳ７のＹ）、各微小領域Ｌでの合焦距離を求めて、全微小領域Ｌでの合焦距離Ｄを求め（ステップＳ８）、合焦距離Ｄから表面形状データＳijを求める（ステップＳ９）。ステップＳ８，Ｓ９により第１の表面形状データ作成手段、第１の表面形状データ作成工程、第１の表面形状データ作成処理を実現している。
【０１３４】
次に、各画素での次数ｍijを計算し（ステップＳ１０）、微小領域Ｌにおける合焦距離Ｄを再生距離ｄとして（ステップＳ１１）、フレネル回折計算を行い、各微小領域Ｌにおける複素振幅を取得する（ステップＳ１２）。ステップＳ１０〜Ｓ１２により、光学データ検出手段、光学データ検出工程、光学データ検出処理を実現している。
【０１３５】
そして、その複素振幅の位相データにより、各画素での位相φの値を表面形状データＴijに変換する（ステップＳ１３）。ステップＳ１３により、第２の表面形状データ作成手段、第２の表面形状データ作成工程、第２の表面形状データ作成処理を実現している。
【０１３６】
最後に、以上のようにして求めた表面形状データＳijと表面形状データＴijとを、前記のように合成又は比較して、表面形状データＺijを求める（ステップＳ１４）。ステップＳ１４により、合成又は比較手段、合成又は比較工程、合成又は比較処理を実現している。
【０１３７】
なお、以上の処理では、振幅データを用いた合焦法を実施した結果、得られた各微小領域Ｌにおける合焦距離Ｄを用いて、フレネル回折計算を繰り返し、各微小領域Ｌでの複素振幅を再計算して（ステップＳ１２）、その位置での位相データを再度求めている。正確に計算するためにはそのほうが望ましいが、フレネル回折計算には時間がかかるため、全画素で位相データを計算するのに多大な時間を要する。そこで、求めた合焦距離Ｄが初期位置に対して微小な場合は、初期位置での位相データを代用しても誤差が小さいため、求めた合焦距離Ｄの値によっては、ステップＳ１２のフレネル回折計算を実施しなくてもよい。これにより演算時間、測定時間の短縮が図れる。
【０１３８】
図１の表面形状測定装置では、キャリヤ周波数をのせた干渉縞を用いて被測定物１４の反射光の複素振幅を求めたが、キャリヤ周波数をのせない同軸の干渉縞においても、参照光、被測定物１４の反射光の位相を変調することによって、被測定物１４の反射光の複素振幅を独立して求めることができる。この場合は、このような複素振幅を用いてもよいし、また、光源であるＨｅ−Ｎｅレーザ１のレーザ光の波長λを大きくした場合は、直接ＣＣＤ１０で被測定物１４の反射光の複素振幅を求められるようになるため、このようにして直接取得した複素振幅を用いてもよい。
【０１３９】
図１の表面形状測定装置によると、ＣＣＤ１０の撮像時間より遅い速度で動く物体に対しては、その動いている瞬間の動的形状の測定が可能になる。ただし、物体の動きが速くなり、ＣＣＤ１０の撮像時間より速くなると、撮像している最中に被測定物１４からの反射光の強度や位相が変化するため、正確に反射光の複素振幅を求めることができなくなる。
【０１４０】
この場合は、被測定物１４の動く速度に対して、十分短い時間で被測定物１４にパルス光を照射し、その反射光を受光するようにすると、光の照射中は、被測定物１４は見かけ上静止している状態とみなせるため、被測定物１４の反射光が安定し、正確にその複素振幅を求めることができて、正確な被測定物１４の動的表面形状測定が可能になる。
【０１４１】
被測定物１４にパルス光を照射する光源には、ルビーレーザやＹＡＧレーザのようにパルス光を出力する固体レーザや、パルス光を出力する半導体レーザを用いることができ、また、ＣＷ光を出力する半導体レーザをパルス変調駆動してパルス光を生成してもよいし、Ｈｅ−ＮｅレーザやアルゴンレーザからのＣＷ光を、回転チョッパを用いてパルス化して用いてもよい。
【０１４２】
また、パルス光を被測定物１４に照射するのに代えて、図１の表面形状測定装置におけるＣＣＤ１０を高速撮像が可能な高速カメラに置き換えることにより、光源をパルス化するのと同様の効果を得ることもできる。
【０１４３】
さらに、（３）式の関係において、光学系倍率は“距離ｓ’／距離ｓ”で決まるため、対物レンズ１３の位置や焦点距離を変えることによって、光学系倍率を変えることができる。光学系倍率を変えることにより、振幅データを用いた表面形状測定方法の測定分解能を変化させることができるため、被測定物１４の大きさに応じて光学系倍率を変化させるようにして、被測定物１４の表面のより広い範囲（視野）を一度に測定できるようにすることができる。
【０１４４】
図１の装置構成に代えて、図８のような装置構成としてもよい。図８において、図１と同一符号の部材は図１の装置と同様であるため、詳細な説明は省略する。
【０１４５】
すなわち、図８の装置構成が図１のものと相違するのは、被測定物１４を光学系の略光軸方向に移動させる拡大倍率可変装置であるステージ２４を設け、これにより被測定物１４と対物レンズ１３との間隔を変化させるようにしている点である。
【０１４６】
被測定物１４と対物レンズ１３との間隔が変化すると（３）式における距離ｓの値が変化し、それに合わせて（３）式が成り立つように距離ｓ’が変化するため、もって光学系の倍率を変化させることができる。
【０１４７】
また、図８の構成は、レンズ６を略光軸方向に移動させるための拡大倍率可変装置であるステージ２５を設けている点も相違する。すなわち、ステージ２４の移動に伴う被測定物１４と対物レンズ１３との間隔の変化に合わせてステージ２５を移動させ、ＣＣＤ１０の撮像位置での被測定物１４の反射光と参照光の曲率をほぼ一致させる。これにより、干渉縞データから計算した位相データにデフォーカス収差（図５参照）が含まれずに、光学系の倍率を変化させることが可能になる。また、距離ｓ’の変化にあわせて、ＣＣＤ１０の位置を変化させるためのステージを設けるようにしてもよい。
【０１４８】
さらに、図１の装置構成に代えて、図９のような装置構成とすることも考えられる。図９において、図１と同一符号の部材は図１の装置と同様であるため、詳細な説明は省略する。
【０１４９】
すなわち、図９の装置構成が図１のものと相違するのは、被測定物１４の反射光を球面波に変換する可変焦点レンズ２６と、参照光を球面波に変換する可変焦点レンズ２７と、レンズ１２をほぼ光軸方向に移動させるための拡大倍率可変装置であるステージ２８とを備えている点である。
【０１５０】
そして、可変焦点レンズ２６とレンズ１２との間にアフォーカル系の関係が成り立つように、可変焦点レンズ２６の焦点距離の変化に応じて、ステージ２８でレンズ１２を移動させることで、可変焦点レンズ２６から被測定物１４に照射される光がほぼ平行化されるようになっている。
【０１５１】
可変焦点レンズ２６又は２７のどちらかの焦点距離を変え、球面波の曲率を変化させることにより、拡大倍率を変化させることができるが、被測定物１４の反射光と参照光のＣＣＤ１０の撮像位置での曲率に差が生じると、それに応じて、干渉縞から計算した位相データにデフォーカス収差が含まれてしまう。そのため、被測定物１４の反射光と参照光との間にＣＣＤ１０の撮像位置での曲率に差が生じないように、可変焦点レンズ２６、２７の焦点距離の変化を連動させるようにする。
【０１５２】
これにより、干渉縞データから計算した位相データにデフォーカス収差が含まれることなく、光学系の倍率を変化させることが可能になる。また、可変焦点レンズ２６又は２７を固定焦点レンズにして、それをほぼ光軸方向に移動可能なステージに搭載して移動させることによって、被測定物１４の反射光と参照光の曲率の差が生じないように調整してもよい。さらに、対物レンズ１３の焦点距離を変えたことにより、（３）式の関係から距離ｓが変化するが、距離ｓ’の変化に応じてＣＣＤ１０の位置を変化させるためのステージを設けてもよい。
【０１５３】
ところで、以上説明した例では、合焦法の原理により求めた表面形状データと、干渉計測の原理により求めた表面形状データＴijとを合成又は比較することによって表面形状データＺijを求めて、光源波長λの１／２を超えるような段差を含む被測定物１４の表面形状をサブフリンジオーダーで測定できるようにしている。
【０１５４】
しかし、合焦法の原理により表面形状データを求める手段としては、前記の例に限定されるものではない。すなわち、前記の例では合焦法の原理により表面形状データＳijを求めるために振幅データを用いたが、これに代えて位相データを２次元画像データとして捉え、これを用いて合焦法を実施してもよい。
【０１５５】
合焦法では、物体像のボケ具合（合焦測度）をコントラストの低下として検知する。しかし、被測定面が鏡面のような平滑面であると、被測定物１４の表面の明確なテクスチャ（模様）を観察できない場合がある。そして、テクスチャが観察できないと、像のボケによるコントラストの低下が検出できなくなるため、合焦法が実施できず、表面形状測定が不可能になる。
【０１５６】
このような場合、振幅データに比べて比較的物体像をコントラストよく観察できる位相データ（位相像）を２次元強度データとして捉えて用いることができる。すなわち、平滑な面であっても、その表面を光学的に拡大していけば、表面粗さが光学系の空間解像力を下回らない限り表面粗さのテクスチャが観察されるようになる。そのため、その像をコントラストの高い位相像として観察することによって、合焦法による形状測定が可能となる。
【０１５７】
そのために、コンピュータ１６で行なう処理としては、図７のフローチャートにおいて、振幅データの微小領域Ｌで合焦法を行なうのに（ステップＳ５）代えて、位相データの微小領域Ｌにおいて合焦法を行なえばよい。
【０１５８】
また、合焦法により表面形状測定を実施する場合には、焦点の合う画像領域を被測定撮像位置で探索していくこととなる。そのため、焦点の合った画像領域をつなぎ合わせることで、全視野で焦点の合った全焦点画像を生成することができる。顕微鏡などのように拡大倍率が大きい光学系では焦点深度が浅くなり、そのため、被測定物１４の位置や姿勢の調整がしづらくなる。このような場合に、全焦点画像を生成することにより、調整作業が容易になる。
【０１５９】
全焦点画像の生成の方法としては、フレネル回折計算を繰り返して複数のデフォーカス画像を取得した後、画像における各微小領域Ｌで合焦距離Ｄを求める。そして、各微小領域Ｌにおいて、求めた合焦距離Ｄと最も近い距離で得られたデフォーカス画像の画素濃度を、表示するための画像における同じ位置の微小領域Ｌに入れる。その作業を画像の全域にわたって繰り返すことにより、画像上のすべての微小領域Ｌで焦点のあった全焦点画像を、高速に取得することができる。
【０１６０】
また、全焦点画像を生成する別の処理方法としては、各微小領域Ｌにおいて求めた合焦点距離Ｄを回折距離として入力し、取得した干渉縞データを用いてフレネル回折計算を実施するようにしてもよい。そして、フレネル回折計算の結果得られた画像における微小領域Ｌの画素濃度を、表示するための画像における同じ位置の微小領域Ｌに入れる。その作業を画像の全域にわたって繰り返すことによって、画像上のすべての微小領域Ｌで、より正確に焦点のあった全焦点画像を取得することができる。
【０１６１】
以上説明した、この実施の形態の表面形状測定装置によれば、次のような作用効果を奏することがわかる。
【０１６２】
すなわち、図１に示す表面形状測定装置によれば、撮像した干渉縞データにより、被測定物１４の反射光のほぼ光軸方向における複数位置での複素振幅を求めて、そのうちの振幅データを用いて、合焦状態の検知により表面形状を求めるdepth from focus理論に基づいた演算処理を施すことによって、被測定物１４の表面形状データＳijを測定する。
【０１６３】
その場合、測定分解能および光軸方向の測定範囲は光学系の拡大倍率により任意に設定することができるため、その測定分解能を光源波長λの１／２より小さく設定することにより、レーザ干渉により測定できない光源波長λの１／２を超える段差を測定することができる。
【０１６４】
一方、干渉縞データから求めた複素振幅のうち、位相データを用いれば、通常のレーザ干渉と同様にしてサブフリンジオーダーの微小形状を測定することができる。したがって、被測定物１４の表面において光源波長λの１／２を超える段差については振幅データを用いた測定で受け持ち（表面形状データＳij）、それ以外の領域の形状測定については位相データを用いた測定で受け持ち（表面形状データＴij）、両測定結果を合成又は比較することにより（表面形状データＺij）、光源波長λの１／２を超える段差をもつ被測定面であっても、サブフリンジオーダーで測定可能となる。
【０１６５】
また、従来の走査型白色干渉計に比べて、光源波長λの１／２を超える段差を含む面を測定するために、干渉光学系や被測定物１４を機械的な駆動で移動させる必要がないため、測定誤差要因が減るうえ、移動範囲の制限によって光軸方向の測定範囲が制限されない。
振幅データを用いた測定については、レンズの焦点距離を変化させる振幅データを用いた従来法に対して、干渉縞データを用いて、光学系の光軸方向における複数位置での被測定物１４の反射光の複素振幅を計算し、そのうちの振幅データを用いてdepth from focus法を実施するため、数値演算でdepth from focus法を実施していることになり、depth from focus法を実施するために、やはり従来のように機械的な駆動を伴わないため、その誤差の影響を受けず、機械的駆動部分の可動範囲で決まる測定範囲の制限も生じない。
【０１６６】
また、測定のために機械的な駆動を要さないため、高速測定が可能であり、動きを伴う物体の動的形状が測定可能である。
【０１６７】
さらに、観測する干渉縞にキャリヤ周波数を与えて単一の干渉縞データから被測定物１４の反射光の複素振幅を計算するようにすると、そのデータが得られるまでに要する時間はＣＣＤ１０による撮像時間のみで決まることになり、前記した複数の画像を撮像する必要がある走査型白色干渉計やレンズ焦点距離を変化させる方法に比べて、高速測定が可能になる。
【０１６８】
このようにして、高速測定が可能になっても、被測定物１４の動きの速度がＣＣＤ１０による撮像時間より速くなると、被測定物１４からの反射光が安定しないため、正確に反射光の複素振幅を求めることができなくなり、正確な測定ができなくなる。そこで、前記のように、被測定物１４にパルス光を照射し、そのパルスの発光時間を被測定物１４の動きに対応して十分短くすることにより、被測定物１４が光を反射する瞬間において、被測定物１４をみかけ上静止した状態にして、被測定物１４からの反射光を安定させることができるため、反射光の複素振幅を正確に求めることができ、正確な測定が可能になる。
【０１６９】
変形、変位、あるいは振動中の物体の形状測定に関し、単一画像から物体反射光の複素振幅を求める方法として、パルス光源を用いたデジタルホログラフィによる物体の動的測定法（例えば、G.Pedrini,H.J.Tiziani,「Quantitative evaluation of two-dimensional dynamic deformations using digital holography」,Optics＆Laser Technology,Vol.29,No.5,pp249-256,1997を参照）があるが、かかる方法においては、粗面物体を測定対象としており、光源波長レベルでの物体形状の測定や、また光源波長λの１／２を超える段差を含む面を光源波長レベルで測定する方法については示されていない。
【０１７０】
かかる方法に比べ、この表面形状測定装置による測定は、被測定物１４が微小であっても、鏡面又はそれに近い面であっても測定可能で、被測定物１４の表面に光源波長λの１／２を超える段差が含まれていても測定可能である。
【０１７１】
また、この表面形状測定装置によれば、光学系倍率を設定することにより、振幅データを用いた表面形状測定と位相データを用いた表面形状測定のうち、振幅データを形状測定についてのみ、その測定分解能を設定することができる。
【０１７２】
一方、振幅データを用いた表面形状測定により光源波長λの１／２の段差を測定するには、その測定分解能は光源波長λの１／２より小さいほうが望ましく、通常の可視光を光源に用いると波長はｎｍオーダーになるため、干渉光学系の拡大倍率を上げたほうがよい。
【０１７３】
その場合、被測定物１４からの反射光又は参照光のどちらか一方を球面波にすることによって干渉光学系の拡大倍率を設定し、前記振幅データを用いた形状測定の測定分解能を光源波長の１／２より小さくすることができるが、被測定物１４からの反射光か参照光かのどちらか一方を球面波にすると、被測定物１４の反射光の波面と参照光の波面との曲率の差に応じたデフォーカス収差が位相データに含まれてくる。デフォーカス収差が含まれると、位相データから正確に被測定物１４の表面形状を求めることができなくなるため、デフォーカス収差を除去する必要がある。
【０１７４】
そのため、この表面形状測定装置では、被測定物１４からの反射光か参照光のどちらか一方を球面波にしたとき、他方の光についても球面波にし、両者の曲率をほぼ一致させることにより、上記のデフォーカス収差を除去することができる。これにより、干渉光学系の拡大倍率をあげた場合でも、測定を正確に実施することができる。
【０１７５】
この表面形状測定装置では、光学系倍率を設定することにより、振幅データを用いた表面形状測定と位相データを用いた表面形状測定のうち、振幅データを用いた表面形状測定についてのみ、その測定分解能を設定する。
【０１７６】
しかし、その測定分解能を過剰に小さくすると光学系倍率が大きくなるため、横方向の測定範囲（視野）は狭くなり、被測定物１４の表面における広い範囲（視野）を一度に測定することができなくなる。
【０１７７】
そこで、この表面形状測定装置では、光源波長の１／２より小さい測定分解能をもつという条件のもとで、被測定物１４の大きさに応じて光学系倍率を変化させるようにし、被測定物１４の表面のより広い範囲（視野）を一度に測定できるようにしている。
【０１７８】
合焦法の原理を用いて、合焦点はずれを検出することにより被測定物１４の表面形状を測定する場合、被測定物１４の表面が平滑で光学的にテクスチャを観測できないと、合焦点はずれの検知が困難になることがある。
【０１７９】
合焦点はずれは被測定物１４の表面のパターンのコントラストの低下により検出されるものである。２次元的な振幅データではコントラストが低下して合焦点はずれが検出できない場合でも、２次元的な位相データではコントラストが得られる場合がある。
【０１８０】
そのような場合は、前記のように位相データを２次元強度データとして用いて、合焦法の原理により被測定物１４の表面形状を測定する。これにより被測定物１４表面が平滑で、振幅データでは光学的にテクスチャを観測できない場合でも、被測定面の形状測定が可能となり、測定の汎用性を向上させることができる。
【０１８１】
物体形状を測定する場合、測定者が被測定物１４を観察しながら測定できたほうが、被測定物１４の設置、姿勢の調整などの点において便利である。被測定物１４が微小になると肉眼で観察することは難しくなるため、光学的に拡大して観察できたほうがよい。しかしながら、拡大光学系は焦点深度が浅く、光軸方向で焦点のあう範囲が狭いため、観察しづらいという欠点がある。
【０１８２】
これに対し、この表面形状測定装置では、振幅データを用いて合焦法の原理により被測定物１４の表面形状を測定するが、その場合、２次元画像内で合焦点の位置を求める処理を実施するため、測定のために取得した合焦点の情報を用いて、２次元的な全領域で焦点のあった全焦点画像を、生成、表示することができ、被測定物１４の位置、姿勢を観察しやすくすることにより、測定の操作性を向上させることができる。
【０１８３】
この場合、高速に画像を生成、表示することが重要な課題となる。この表面形状測定装置では、全焦点画像における微小領域の焦点画像を、合焦状態から求めた合焦距離に対して最も近い距離にある振幅データに置き換えて生成することによって、高速に全焦点画像を生成することができるので、より自然な被測定物１４の観測が可能となり、測定の操作性を向上させることができる。
【０１８４】
また、全焦点画像を生成する場合、全焦点画像における微小領域の焦点画像を、合焦状態から求めた合焦距離を用いてフレネル回折計算することにより、より正確な全焦点画像の生成が可能となり、測定における被測定物１４の位置、姿勢の調整が、より正確に実施可能となるので、測定の高速化、高精度化を図ることができる。
【０１８５】
【発明の効果】
請求項１に記載の発明は、被測定物の光源波長の１／２を超える段差については第１の表面形状データで受け持ち、それ以外の領域の形状測定については第２の表面形状データで受け持って、この両者を合成又は比較することにより、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【０１９０】
請求項２に記載の発明は、請求項１に記載の表面形状測定装置において、測定分解能を光源波長の１／２より小さく設定することにより、光源波長の１／２を超える段差も測定することができる。
【０１９１】
請求項３に記載の発明は、請求項１又は２に記載の表面形状測定装置において、反射光の光軸方向の複数位置は、それぞれの位置の間隔が光源の波長の１／２より小さくなるように設定しているので、光源波長の１／２を超える段差も測定することができる。
【０１９２】
請求項４に記載の発明は、請求項１〜３のいずれかの一に記載の表面形状測定装置において、第１と第２の表面形状データを的確に合成して、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【０１９３】
請求項５に記載の発明は、請求項１〜４のいずれかの一に記載の表面形状測定装置において、観測する干渉縞にキャリヤ周波数を与えて、単一の干渉縞データから被測定物の反射光の振幅データ、位相データを計算することにより、測定の高速化を図ることができる。
【０１９４】
請求項６に記載の発明は、請求項１〜５のいずれかの一に記載の表面形状測定装置において、パルスの発光時間を被測定物の動きに対応して十分短くすることで、被測定物が光を反射する瞬間に被測定物がみかけ上静止した状態となり、被測定物からの反射光を安定させて、正確に反射光の振幅データ、位相データを求めることができるので、被測定物の動きが高速な場合においても、正確に物体の動的形状を測定することができる。
【０１９５】
請求項７に記載の発明は、請求項１〜６のいずれかの一に記載の表面形状測定装置において、デフォーカス収差を除去して、干渉光学系の拡大倍率を高めた測定を正確に実施することができる。
【０１９６】
請求項８に記載の発明は、請求項１〜７のいずれかの一に記載の表面形状測定装置において、光源波長の１／２より小さい測定分解能のもとで、被測定物の大きさに応じて干渉光学系の拡大倍率を変化させるようにすれば、被測定物表面のより広い範囲を一度に測定することができる。
【０１９７】
請求項９に記載の発明は、請求項８に記載の表面形状測定装置において、被測定物表面のより広い範囲を一度に測定することを容易に行なうことができる。
【０１９８】
請求項１０に記載の発明は、被測定物の光源波長の１／２を超える段差については第１の表面形状データで受け持ち、それ以外の領域の形状測定については第２の表面形状データで受け持って、この両者を合成又は比較することにより、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【０２０３】
請求項１１に記載の発明は、被測定物の光源波長の１／２を超える段差については第１の表面形状データで受け持ち、それ以外の領域の形状測定については第２の表面形状データで受け持って、この両者を合成又は比較することにより、被測定物や光学系を測定中に移動させることなく、光源波長の１／２を超えるような段差を含む被測定物の表面もサブフリンジオーダーで高精度かつ容易に測定することができる。
【０２０８】
請求項１２記載の発明は、請求項１１に記載の発明と同様の作用、効果を奏する。
【図面の簡単な説明】
【図１】この発明の一実施の形態である表面形状測定装置の概略構成を示すブロック図である。
【図２】前記表面形状測定装置における光学系の位置関係を示す説明図である。
【図３】同説明図である。
【図４】前記表面形状測定装置が行なう処理で再生距離と合焦測度との関係を示すグラフである。
【図５】デフォーカス収差について説明する説明図である。
【図６】前記表面形状測定装置を構成するコンピュータの電気的な接続を示すブロック図である。
【図７】コンピュータが実行する処理のフローチャートである。
【図８】前記表面形状測定装置の他の例について概略構成を示すブロック図である。
【図９】同ブロック図である。
【符号の説明】
１ 光源
５〜９，１１〜１３ 干渉光学系
１０ 撮像素子
１４ 被測定物
２４，２５，２８ 拡大倍率可変装置
３７ 記憶媒体
Claims (12)
- 被測定物に照射する光を発する光源と、
その光の被測定物における反射光と所定の参照光との間で干渉縞を発生させる干渉光学系と、
この干渉縞を撮像して当該干渉縞の画像データである干渉縞データを出力する撮像素子と、
前記干渉縞データに基づいて前記反射光の光軸方向の複数位置における前記被測定物の反射光の位相を示す位相データ及び振幅を示す振幅データを求める光学データ検出手段と、
前記振幅データの微分値の和である合焦測度を算出し、前記合焦測度が最大となる合焦距離に基づき、前記被測定物の第１の表面形状データを求める第１の表面形状データ作成手段と、
前記位相データを長さ単位に変換することで、前記被測定物の第２の表面形状データを求める第２の表面形状データ作成手段と、
前記第１の表面形状データと前記第２の表面形状データとを合成又は比較して、前記被測定物の表面形状のデータである第３の表面形状データを求める合成又は比較手段と、を備えている表面形状測定装置。 - 前記干渉光学系は、その拡大倍率の２乗に反比例する感度で得られる前記振幅データによる前記被測定物の表面形状における前記反射光の光軸方向の測定分解能を、前記光源の波長の１／２より小さくなるように所定の拡大倍率に設定されている請求項１に記載の表面形状測定装置。
- 前記データ検出手段は、前記反射光の光軸方向の複数位置を、それぞれの位置の間隔が前記光源の波長の１／２より小さいものとしている請求項１又は２に記載の表面形状測定装置。
- 前記合成又は比較手段は、前記合成として、
ｍ×（λ／２）＋δ
（ｍ：次数，λ：光源波長，δ：前記第２の表面形状データ）における次数ｍを決定する処理を行うものである請求項１〜３のいずれかの一に記載の表面形状測定装置。 - 前記光学データ検出手段は、前記位相データ及び振幅データをキャリヤ周波数を与えた単一の前記干渉縞の干渉縞データから求める請求項１〜４のいずれかの一に記載の表面形状測定装置。
- 前記光源は、前記被測定物にパルス光を照射するものである請求項１〜５のいずれかの一に記載の表面形状測定装置。
- 前記変換光学系は、前記反射光及び前記参照光をそれぞれ球面波に変換し、前記撮像素子による干渉縞の撮像位置での前記両球面波の曲率をほぼ一致させて前記干渉光学系の倍率を設定している請求項１〜６のいずれかの一に記載の表面形状測定装置。
- 前記干渉光学系の拡大倍率を可変する拡大倍率可変装置を備え、前記第１の表面形状データ作成手段及び前記第２の表面形状データ作成手段のうち前者による前記被測定物の測定分解能又は前記反射光の光軸方向の測定範囲のみを設定できるようにしている請求項１〜７のいずれかの一に記載の表面形状測定装置。
- 前記拡大倍率可変装置は、前記撮像素子による前記干渉縞の撮像位置における前記反射光の曲率と前記参照光の曲率とがほぼ一致するよう連動させて曲率を変化させることにより、前記干渉光学系の倍率を可変としている請求項８に記載の表面形状測定装置。
- 被測定物に光を照射したときの反射光と所定の参照光との間で発生させた干渉縞を撮像して当該干渉縞の画像データである干渉縞データを取得する撮像工程と、
前記干渉縞データに基づいて前記反射光の光軸方向の複数位置における前記被測定物の反射光の位相を示す位相データ及び振幅を示す振幅データを求める光学データ検出工程と、
前記振幅データの微分値の和である合焦測度を算出し、前記合焦測度が最大となる合焦距離に基づき、前記被測定物の第１の表面形状データを求める第１の表面形状データ作成工程と、
前記位相データを長さ単位に変換することで、前記被測定物の第２の表面形状データを求める第２の表面形状データ作成工程と、
前記第１の表面形状データと前記第２の表面形状データとを合成又は比較して、前記被測定物の表面形状のデータである第３の表面形状データを求める合成又は比較工程と、を含んでなる表面形状測定方法。 - 被測定物に光を照射したときの反射光と所定の参照光との間で発生させた干渉縞を撮像した画像データである干渉縞データに基づいて前記反射光の光軸方向の複数位置における前記被測定物の反射光の位相を示す位相データ及び振幅を示す振幅データを求める光学データ検出処理と、
前記振幅データの微分値の和である合焦測度を算出し、前記合焦測度が最大となる合焦距離に基づき、前記被測定物の第１の表面形状データを求める第１の表面形状データ作成処理と、
前記位相データを長さ単位に変換することで、前記被測定物の第２の表面形状データを求める第２の表面形状データ作成処理と、
前記第１の表面形状データと前記第２の表面形状データとを合成又は比較して、前記被測定物の表面形状のデータである第３の表面形状データを求める合成又は比較処理と、をコンピュータに実行させるコンピュータに読み取り可能なプログラム。 - 請求項１１に記載のプログラムを記憶した記憶媒体。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001289786A JP4782958B2 (ja) | 2001-09-21 | 2001-09-21 | 表面形状測定装置及びその方法、プログラム並びに記憶媒体 |
| US10/251,830 US6885460B2 (en) | 2001-09-21 | 2002-09-23 | Apparatus for and method of measuring surface shape of an object |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001289786A JP4782958B2 (ja) | 2001-09-21 | 2001-09-21 | 表面形状測定装置及びその方法、プログラム並びに記憶媒体 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2003097922A JP2003097922A (ja) | 2003-04-03 |
| JP4782958B2 true JP4782958B2 (ja) | 2011-09-28 |
Family
ID=19112231
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001289786A Expired - Fee Related JP4782958B2 (ja) | 2001-09-21 | 2001-09-21 | 表面形状測定装置及びその方法、プログラム並びに記憶媒体 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US6885460B2 (ja) |
| JP (1) | JP4782958B2 (ja) |
Families Citing this family (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005156434A (ja) * | 2003-11-27 | 2005-06-16 | Fujinon Corp | 計算機ホログラムを用いた光波干渉測定方法およびこれを用いた干渉計装置 |
| JP2005315850A (ja) * | 2004-03-29 | 2005-11-10 | Kgt Inc | 三次元位置計測装置 |
| DE102004017232A1 (de) * | 2004-04-05 | 2005-10-20 | Bosch Gmbh Robert | Interferometrische Messvorrichtung |
| JP4514209B2 (ja) * | 2004-10-15 | 2010-07-28 | キヤノン株式会社 | 位置検出装置及び方法 |
| US7436520B1 (en) * | 2005-01-18 | 2008-10-14 | Carl Zeiss Smt Ag | Method of calibrating an interferometer optics and of processing an optical element having an optical surface |
| US7884947B2 (en) * | 2005-01-20 | 2011-02-08 | Zygo Corporation | Interferometry for determining characteristics of an object surface, with spatially coherent illumination |
| CN100464172C (zh) * | 2005-04-08 | 2009-02-25 | 中国科学院上海光学精密机械研究所 | 数字干涉位相差放大方法 |
| WO2007002898A2 (en) * | 2005-06-29 | 2007-01-04 | University Of South Florida | Variable tomographic scanning with wavelength scanning digital interface holography |
| US20070034297A1 (en) * | 2005-07-14 | 2007-02-15 | Greg Zielke | Lumber processing system |
| US8124280B2 (en) * | 2006-12-22 | 2012-02-28 | Uchicago Argonne, Llc | Intermetallic electrodes for lithium batteries |
| TWI366749B (en) * | 2007-12-10 | 2012-06-21 | Ind Tech Res Inst | Multi-color off-axis digital holographic system and the imaging method thereof |
| JP5174697B2 (ja) * | 2008-01-31 | 2013-04-03 | 株式会社日立ハイテクノロジーズ | 欠陥検査装置 |
| JP5019478B2 (ja) * | 2008-09-26 | 2012-09-05 | 独立行政法人日本原子力研究開発機構 | マーカ自動登録方法及びシステム |
| TWI391624B (zh) * | 2008-11-28 | 2013-04-01 | Univ Nat Sun Yat Sen | Striped Projection Profiler |
| EE05614B1 (et) * | 2008-12-16 | 2012-12-17 | Nikolay Voznesenskiy | Interferomeeter optiliste süsteemide ja näidiste testimiseks |
| JP5439218B2 (ja) * | 2010-02-17 | 2014-03-12 | 株式会社ミツトヨ | 画像測定機 |
| GB201205491D0 (en) | 2012-03-28 | 2012-05-09 | Ucl Business Plc | Measuring surface curvature |
| JP6193546B2 (ja) * | 2012-08-24 | 2017-09-06 | ウシオ電機株式会社 | 三次元位置測定方法、速度測定方法、三次元位置測定装置及び速度測定装置 |
| WO2016035108A1 (ja) * | 2014-09-03 | 2016-03-10 | オリンパス株式会社 | 撮像装置 |
| JP6572099B2 (ja) * | 2015-11-06 | 2019-09-04 | キヤノン株式会社 | 撮像装置及びその制御方法、ならびにプログラム |
| JP6246875B1 (ja) * | 2016-08-24 | 2017-12-13 | Ckd株式会社 | 計測装置 |
| JP6805930B2 (ja) * | 2017-03-29 | 2020-12-23 | 株式会社島津製作所 | 振動測定装置 |
| JP6423032B2 (ja) * | 2017-03-29 | 2018-11-14 | 株式会社トプコン | ３次元測量装置 |
| WO2021188277A1 (en) * | 2020-03-17 | 2021-09-23 | Roof Asset Management Usa Ltd. | Method for evaluating artificial lighting of a surface |
| EP4361685A1 (en) * | 2022-10-25 | 2024-05-01 | Nextlens Switzerland AG | An optomechanical assembly and operating procedure of said assembly |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5321497A (en) * | 1992-03-09 | 1994-06-14 | Wyko Corporation | Interferometric integration technique and apparatus to confine 2π discontinuity |
| JP3237309B2 (ja) * | 1992-06-17 | 2001-12-10 | キヤノン株式会社 | システムエラー測定方法及びそれを用いた形状測定装置 |
| TW283201B (ja) * | 1993-08-13 | 1996-08-11 | Ricoh Kk | |
| US5844670A (en) * | 1995-07-28 | 1998-12-01 | Ricoh Co., Ltd. | Method of and systems for measuring eccentricity of an aspherical lens surface |
| CN1272622C (zh) * | 1998-04-22 | 2006-08-30 | 株式会社理光 | 双折射测定方法及其装置 |
| US6344898B1 (en) * | 1998-09-14 | 2002-02-05 | Nikon Corporation | Interferometric apparatus and methods for measuring surface topography of a test surface |
| JP2000097666A (ja) * | 1998-09-22 | 2000-04-07 | Nikon Corp | 面形状計測用干渉計、波面収差測定機、前記干渉計及び前記波面収差測定機を用いた投影光学系の製造方法、及び前記干渉計の校正方法 |
| US6312373B1 (en) * | 1998-09-22 | 2001-11-06 | Nikon Corporation | Method of manufacturing an optical system |
| US6262818B1 (en) * | 1998-10-07 | 2001-07-17 | Institute Of Applied Optics, Swiss Federal Institute Of Technology | Method for simultaneous amplitude and quantitative phase contrast imaging by numerical reconstruction of digital holograms |
| JP2000221013A (ja) * | 1999-02-03 | 2000-08-11 | Fuji Xerox Co Ltd | 干渉計測方法および干渉計測装置 |
| JP4386543B2 (ja) | 1999-05-26 | 2009-12-16 | 株式会社神戸製鋼所 | 表面形状測定装置 |
| JP4183219B2 (ja) * | 1999-12-21 | 2008-11-19 | フジノン株式会社 | フーリエ変換を用いた縞解析方法 |
-
2001
- 2001-09-21 JP JP2001289786A patent/JP4782958B2/ja not_active Expired - Fee Related
-
2002
- 2002-09-23 US US10/251,830 patent/US6885460B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| US20030067609A1 (en) | 2003-04-10 |
| US6885460B2 (en) | 2005-04-26 |
| JP2003097922A (ja) | 2003-04-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4782958B2 (ja) | 表面形状測定装置及びその方法、プログラム並びに記憶媒体 | |
| Ma et al. | Numerical reconstruction of digital holograms for three-dimensional shape measurement | |
| JP5214883B2 (ja) | 三次元分光的符号化撮像のための方法と装置 | |
| US10845760B2 (en) | Digital holographic reconstruction device and method using single generation phase shifting method | |
| JP7231946B2 (ja) | 表面形状計測装置および表面形状計測方法 | |
| CN107121084B (zh) | 测量方法和测量程序 | |
| JP2005539256A (ja) | 複合画像間の差異を検出するためのシステムおよび方法 | |
| JP3816402B2 (ja) | 表面形状測定装置及び表面形状測定方法 | |
| JP4100553B2 (ja) | 動的形状及び動的位置の同時測定装置及び同時測定方法 | |
| KR101505745B1 (ko) | 이중 검출 반사 공초점 현미경 및 이를 사용하는 시편의 높이의 정보를 검출하는 방법 | |
| KR20170023363A (ko) | 디지털 홀로그래피 마이크로스코프를 이용한 고단차 측정 방법 | |
| JP2012145361A (ja) | デジタルホログラフィ装置 | |
| JP4069204B2 (ja) | デジタルホログラフィを利用した変位分布計測方法 | |
| JP3958099B2 (ja) | ホログラフィ装置 | |
| JP4025878B2 (ja) | 物体の再生像を得る装置、位相シフトデジタルホログラフィ変位分布計測装置及びパラメータを同定する方法 | |
| JP4023666B2 (ja) | 動的形状測定方法及び装置 | |
| KR20170079441A (ko) | 진동 환경에서 큰 단차를 갖는 샘플의 높낮이 측정을 위한 디지털 홀로그램 기록재생장치 및 기록재생방법 | |
| JP2022162306A (ja) | 表面形状計測装置および表面形状計測方法 | |
| JP4177188B2 (ja) | 動的形状及び動的位置の同時測定方法、装置、光学素子 | |
| JP2020153992A (ja) | 白色干渉計による形状測定装置 | |
| JP3233723B2 (ja) | 位相パターン差判別装置 | |
| Shin et al. | Frequency domain zero padding for accurate autofocusing based on digital holography | |
| JPH07190737A (ja) | シアリング干渉計測方法及びシアリング干渉計 | |
| KR20240051022A (ko) | 디지털 홀로그래픽 현미경을 제어하는 방법 및 장치 | |
| JP2002228426A (ja) | 曲率半径測定方法及び装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20041005 |
|
| RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7421 Effective date: 20051021 |
|
| RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7421 Effective date: 20080111 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080902 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20110323 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110329 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110524 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20110614 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20110708 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140715 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |