JP4566685B2 - 光画像計測装置及び光画像計測方法 - Google Patents

光画像計測装置及び光画像計測方法 Download PDF

Info

Publication number
JP4566685B2
JP4566685B2 JP2004299036A JP2004299036A JP4566685B2 JP 4566685 B2 JP4566685 B2 JP 4566685B2 JP 2004299036 A JP2004299036 A JP 2004299036A JP 2004299036 A JP2004299036 A JP 2004299036A JP 4566685 B2 JP4566685 B2 JP 4566685B2
Authority
JP
Japan
Prior art keywords
light
auxiliary
measurement
signal
measurement depth
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
Application number
JP2004299036A
Other languages
English (en)
Other versions
JP2006112864A (ja
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.)
Topcon Corp
Original Assignee
Topcon Corp
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
Application filed by Topcon Corp filed Critical Topcon Corp
Priority to JP2004299036A priority Critical patent/JP4566685B2/ja
Priority to US11/246,394 priority patent/US7492466B2/en
Priority to CN200510112767.4A priority patent/CN100582741C/zh
Priority to EP05022289A priority patent/EP1650528A3/en
Publication of JP2006112864A publication Critical patent/JP2006112864A/ja
Application granted granted Critical
Publication of JP4566685B2 publication Critical patent/JP4566685B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02001Interferometers characterised by controlling or generating intrinsic radiation properties
    • G01B9/02007Two or more frequencies or sources used for interferometric measurement
    • 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
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02001Interferometers characterised by controlling or generating intrinsic radiation properties
    • G01B9/02002Interferometers characterised by controlling or generating intrinsic radiation properties using two or more frequencies
    • G01B9/02003Interferometers characterised by controlling or generating intrinsic radiation properties using two or more frequencies using beat frequencies
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02001Interferometers characterised by controlling or generating intrinsic radiation properties
    • G01B9/02012Interferometers characterised by controlling or generating intrinsic radiation properties using temporal intensity variation
    • G01B9/02014Interferometers characterised by controlling or generating intrinsic radiation properties using temporal intensity variation by using pulsed light
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02015Interferometers characterised by the beam path configuration
    • G01B9/02027Two or more interferometric channels or interferometers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02041Interferometers characterised by particular imaging or detection techniques
    • G01B9/02044Imaging in the frequency domain, e.g. by using a spectrometer
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/02055Reduction or prevention of errors; Testing; Calibration
    • G01B9/02062Active error reduction, i.e. varying with time
    • G01B9/02067Active error reduction, i.e. varying with time by electronic control systems, i.e. using feedback acting on optics or light
    • G01B9/02069Synchronization of light source or manipulator and detector
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00Measuring instruments characterised by the use of optical techniques
    • G01B9/02Interferometers
    • G01B9/0209Low-coherence interferometers
    • G01B9/02091Tomographic interferometers, e.g. based on optical coherence
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B2290/00Aspects of interferometers not specifically covered by any group under G01B9/02
    • G01B2290/45Multiple detectors for detecting interferometer signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B2290/00Aspects of interferometers not specifically covered by any group under G01B9/02
    • G01B2290/65Spatial scanning object beam
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B2290/00Aspects of interferometers not specifically covered by any group under G01B9/02
    • G01B2290/70Using polarization in the interferometer

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Optics & Photonics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Eye Examination Apparatus (AREA)

Description

本発明は、特に光散乱媒質の被測定物体に光ビームを照射し、その反射光もしくは透過光を用いて被測定物体の表面形態や内部形態を計測し、その画像を形成する光画像計測装置及び光画像計測方法に関する。特に、光ヘテロダイン検出法を利用して被測定物体の表面形態や内部形態を計測し画像を形成する光画像計測装置及び光画像計測方法に関するものである。
近年、レーザ光源等を用いて被測定物体の表面や内部の画像を形成する光画像計測技術が注目を集めている。この光画像計測技術は、従来からのX線CTのような人体に対する有害性を持たないことから、特に医療分野における応用の展開が期待されている。
光画像計測技術における代表的な手法の一例として低コヒーレンス干渉法(光コヒーレンス断層画像化法などとも呼ばれる)がある。この手法は、例えばスーパールミネセントダイオード(Super Luminescent Diode;SLD)のように広いスペクトル幅をもつ広帯域光源の低干渉性を利用するもので、被測定物体からの反射光や透過光をμmオーダーの優れた距離分解能で検出可能とするものである(例えば下記の非特許文献1を参照)。
この低コヒーレンス干渉法を利用した装置の一例として、マイケルソン型の干渉計に基づく従来の光画像計測装置の基本構成を図8に示す。この光画像計測装置200は、広帯域光源201、鏡202、ビームスプリッタ203及び光検出器204を含んで構成されている。被測定物体205は、散乱媒質により形成されている。広帯域光源201からの光ビームは、ビームスプリッタ203により、鏡202に向かう参照光Rと被測定物体205に向かう信号光Sの2つに分割される。参照光Rはビームスプリッタ203による反射光であり、信号光Sはビームスプリッタ203の透過光である。
ここで、図8に示すように、信号光Sの進行方向をz軸として定めるとともに、信号光Sの進行方向に対する直交面をx−y面として定義する。鏡202は、同図中の両側矢印方向に変位(z−スキャン)可能とされている。
参照光Rは、鏡202に反射される際にそのz−スキャンによりドップラー周波数シフトを受ける。一方、信号光Sは、被測定物体205に照射されると、その表面及び内部層により反射される。被測定物体は散乱媒質であるので、信号光Sの反射光は多重散乱を含む乱雑な位相をもった拡散波面であると考えられる。被測定物体205を経由した信号光と、鏡202を経由し周波数シフトを受けた参照光は、ビームスプリッタ203によって重畳されて干渉光を生成する。
低コヒーレンス干渉方法を用いた画像計測では、信号光Sと参照光Rの光路長差が光源のμmオーダーのコヒーレント長(可干渉距離)以内でありかつ参照光Rと位相相関のある信号光Sの成分のみが参照光Rと干渉を生じる。すなわち、信号光Sのコヒーレントな信号光成分のみが選択的に参照光Rと干渉し合う。この原理から、鏡202の位置を移動させて参照光Rの光路長を変化させることにより、被測定物体205の様々なz座標(計測深度)における反射光の情報を含み干渉光が生成され、それにより被測定物体205の内部層の光反射プロフィールが測定される。更に、被測定物体205へ照射される信号光Sをx−y面方向にも走査する。このようなz方向及びx−y面方向のスキャンを行いながら干渉光を光検出器204で検出し、その検出結果として出力される電気信号(ヘテロダイン信号)を解析することによって、被検体205の2次元断層画像が取得される(例えば、非特許文献1を参照)。
なお、ビームスプリッタ203によって重畳される参照光R及び信号光Sの強度をそれぞれIとIとし、両光波の間の周波数差及び位相差をそれぞれfif及びΔθとすると、光検出器からは次式に示すようなヘテロダイン信号が出力される(例えば、非特許文献2を参照)。
Figure 0004566685
式(1)の右辺第3項は交流電気信号であり、その周波数fifは参照光Rと信号光Sのうなり(ビート)の周波数に等しい。ヘテロダイン信号の交流成分の周波数fifは、ビート周波数などと呼ばれる。また、式(1)の右辺第1項及び第2項はヘテロダイン信号の直流成分であり、干渉光の背景光の信号強度に対応している。
しかしながら、このような従来の低コヒーレンス干渉法で2次元断層画像を取得するためには、被測定物体205に対して光ビームを走査することにより、被測定物体205の計測深度方向(z方向)及び断層面方向(x−y面上)の各部位からの反射光波を順次に検出する必要がある。したがって、被測定物体205を計測するためには長時間を要し、また、その計測原理を勘案すると計測時間の短縮を図ることは困難である。
このような問題点に鑑み、計測時間の短縮を図るための光画像計測装置が考案されている。図9は、そのような装置の一例の基本構成が示されている。同図に示す光画像計測装置300は、広帯域光源301、鏡302、ビームスプリッタ303、光検出器としての2次元光センサアレイ304、及びレンズ306、307を含んで構成されている。光源301から出射された光ビームは、レンズ306、307により平行光束にされるとともにそのビーム径が拡大され、ビームスプリッタ303によって参照光Rと信号光Sとに二分される。参照光Rには鏡302のz−スキャンによりドップラー周波数シフトが付与される。一方、信号光Sは、ビーム径が広げられているので、x−y面の広い範囲に亘って被測定物体305に入射される。よって、信号光Sは、当該入射範囲における被測定物体305の表面や内部の情報を含んだ反射光となる。参照光Rと信号光Sは、ビームスプリッタ303により重畳され、2次元光センサアレイ304上に並列配置された素子(光センサ)により検出される。したがって、光ビームをx−y方向に走査させることなく、被測定物体305の2次元断層画像をリアルタイムに取得することが可能となる。
このような非走査型の光画像計測装置としては、非特許文献3に記載のものが知られている。同文献に記載の装置では、2次元光センサアレイから出力される複数のヘテロダイン信号を、並列配置された複数の信号処理系に入力して、各ヘテロダイン信号の振幅と位相を検出するように構成されている。
しかし、このような構成において画像の空間分解能を高めるためにはアレイの素子数を増加させる必要があり、更に、その素子数に対応するチャンネル数を備えた信号処理系を用意しなければならない。したがって、高分解能の画像を必要とする医療分野や工業分野等においては実用化する上で難があると思われる。
そこで、本発明者らは、下記の特許文献1において、次のような非走査型の光画像計測装置を提案した。当該文献に記載の光画像計測装置は、光ビームを出射する光源と、該光源から出射された光ビームを、被検体が配置される被検体配置位置を経由する信号光と、前記被検体配置位置を経由する光路とは異なる光路を経由する参照光とに二分するとともに、前記被検体配置位置を経由した後の信号光と、前記異なる光路を経由した参照光とを互いに重畳することにより干渉光を生成する干渉光学系と、該干渉光学系が、前記信号光の周波数と前記参照光の周波数を相対的にシフトさせる周波数シフタと、前記干渉光学系が、前記干渉光を受光するために、前記干渉光を二分割して、さらに、該二分割された干渉光を周期的に遮断することにより、互いの位相差が90度である2列の干渉光パルスを生成する光遮断装置と、前記2列の干渉光パルスをそれぞれ受光する光センサと、該光センサが、空間的に配列され、それぞれが独立に受光信号を得る複数の受光素子を有するものであり、前記光センサで得られた複数の受光信号を統合して前記被検体配置位置に配置された被検体の表面もしくは内部層の、前記信号光の伝搬経路上の各関心点に対応する信号を生成する信号処理部を具備している。
この光画像計測装置は、参照光と信号光の干渉光を二分して2台の光センサ(2次元光センサアレイ)で受光するとともに、両センサアレイの前にそれぞれ光遮断装置を配置して干渉光をサンプリングするように構成されている。そして、分割された2つの干渉光のサンプリング周期にπ/2の位相差を設けることにより、干渉光の背景光を構成する信号光と参照光の強度と、干渉光の位相の直交成分(sin成分とcos成分)とを検出するとともに、両センサアレイからの出力に含まれる背景光の強度を両センサアレイの出力から差し引くことにより、干渉光の2つの位相直交成分を算出し、その算出結果を用いて干渉光の振幅を求めるようになっている。
なお、以上のような光画像計測装置の2次元光センサアレイとしてはCCD(Charge−Coupled Device)カメラなどの市販のイメージセンサが広く用いられている。しかし、現在市販されているCCDカメラは周波数応答特性が低く、数KHzから数MHz程度のヘテロダイン信号のビート周波数に追従できないという問題点が従来から認識されていた。特許文献1記載の光画像計測装置は、当該問題点を十分に認識した上で、その応答特性の低さを利用して計測を行っている点が特徴的であるといえる。
特許文献1記載のような光画像計測装置において、x−y断層像1フレームを取得するには1ミリ秒程度の時間が掛かる。また、被測定物体の3次元画像、計測深度方向の断層画像(x−z断層画像、y−z断層画像)、あるいは或る軸に斜交する方向の断層画像を形成するためには、多くの計測深度(z座標)におけるx−y断層画像を取得する必要があり、また、被測定物体の3次元画像等を好適な精度で取得するためには、例えば5マイクロメートル程度の一定の間隔で計測深度方向(z方向)の走査を行う必要がある。なお、被測定物体の3次元画像等は、各x−y断層画像が当該一定間隔で計測されたものと仮定して形成される。
このように、3次元画像等を形成するためには、被測定物体の種類等にも依るが、1秒ないしはそれ以上の時間を要する。しかし、人眼等の生体組織などの可動な被測定物体を計測する場合、計測中に被測定物体が動いてしまうと、各x−y断層画像の計測位置にずれが生じて3次元画像等の精度が劣化するおそれがある。
特に、計測深度方向におけるx−y断層画像の計測間隔にずれが生じると、3次元画像等の形成に係る上記の仮定が崩れることとなるため、画像の精度が著しく低下してしまう。そこで、計測深度方向(z方向)におけるx−y断層画像を関係付ける何らかのプロファイルが必要となる。
このような被測定物体のz方向の情報を取得する手法としては、例えば図10に表すような光計測装置が公知となっている(例えば、非特許文献4を参照)。同図に示す光計測装置400は、広帯域光源401からの光ビームをハーフミラー402により、被測定物体404に向かう信号光と固定配置された鏡(固定鏡)403に向かう参照光とに分割し、被測定物体404及び固定鏡403によりそれぞれ反射された信号光及び参照光を重畳して干渉光を生成し、この干渉光を回折格子405によって波長λ1〜λnの波長成分に分光して、各波長成分を1次元光センサアレイ406によって検出するようになっている。1次元光センサアレイ406を構成する各光センサは、検出した波長成分の光強度の検出信号をコンピュータ407に出力する。
コンピュータ407は、1次元光センサアレイ406から出力された干渉光の各波長成分の検出信号に基づいて、この干渉光についての波長−光強度の関係、すなわち干渉光の光強度分布(波長スペクトル)を求める。図11(A)は、干渉光の波長スペクトルの一例の概略を表すグラフである。
更に、コンピュータ407は、求めた干渉光の波長スペクトルをフーリエ変換する。それにより、図11(B)に示すような、被測定物体404のz座標(計測深度)を変数とした干渉信号の強度分布が得られる。これは、被測定物体404の計測深度に依存する情報となっている。なお、1次元光センサアレイ406は、通常、MHz以上(すなわちμsec以下)の読み出しスピードを有しており、計測深度を変数とした干渉信号の強度分布の取得も同等の速さで行うことができる。
特開2001−330558号公報(請求項、明細書段落[0068]−[0084]、第1図) 丹野直弘、「光学」、28巻3号、116(1999) 吉沢、瀬田編、「光ヘテロダイン技術(改訂版)」、新技術コミュニケーションズ(2003)、p.2 K.P.Chan、M.Yamada、H.Inaba、「Electronics Letters」、Vol.30 1753、(1994) A.F.Fercher、C.K.Hitzenberger、G.Kamp、S.Y.Elzaiat、「Optics Communication」Vol.117、pp.43−48(1995)
本発明は、以上のような事情に鑑みてなされたもので、計測中に被測定物体が動いたり移動してしまったような場合であっても、被測定物体の様々な深さ(z座標)の領域について取得されたx−y断層画像に基づいて高い精度の3次元画像等を形成することが可能な光画像計測装置及び光画像計測方法を提供することを目的とする。
上記目的を達成するために、請求項1に記載の発明は、低コヒーレント光源から出力された光ビームを被測定物体に向かう信号光と参照物体に向かう参照光とに分割する分割手段と、前記信号光の周波数と前記参照光の周波数とを相対的にシフトさせる周波数シフト手段と、前記参照光の光路長を変更する変更手段と、前記被測定物体を経由した前記信号光と前記参照物体を経由した前記参照光とを重畳させて前記周波数のシフト量に応じたビート周波数の干渉光を生成する重畳手段と、前記生成された干渉光を受光して検出信号を出力する検出手段と、前記出力された検出信号に基づいて前記信号光の進行方向に直交する前記被測定物体の断層画像を形成する画像形成手段とを有する光画像計測装置であって、前記形成される前記被測定物体の前記断層画像について、前記信号光の進行方向における計測深度を求める取得手段と、前記断層画像と、前記断層画像について取得された前記計測深度の情報とを関連付けて記憶する情報記憶手段と、を備え、前記取得手段は、低コヒーレントな補助光ビームを出力する補助光源と、前記出力された補助光ビームを、前記被測定物体を経由する補助信号光と、補助参照物体を経由する補助参照光とに分割する補助分割手段と、前記被測定物体を経由した前記補助信号光と、前記補助参照物体を経由した前記補助参照光とを重畳して補助干渉光を生成する補助重畳手段と、前記生成された前記補助干渉光から複数の波長成分を分離する波長成分分離手段と、前記分離された前記複数の波長成分を受光して補助検出信号を出力する補助検出手段と、前記検出された前記補助検出信号に基づいて前記断層画像の前記計測深度を算出する計測深度算出手段と、を備えることを特徴とする。
また、請求項2に記載の発明は、前記情報記憶手段に記憶された前記断層画像と関連付けられた計測深度の情報に基づいて、前記断層画像を前記計測深度方向に配列させる画像処理手段を更に備えることを特徴とする。
また、請求項に記載の発明は、請求項1又は請求項に記載の光画像計測装置であって、前記取得手段は、前記補助分割手段により分割された前記補助信号光を前記信号光に合成して前記被測定物体に入射させ、前記被測定物体を経由した前記補助信号光を前記信号光から分離して前記補助重畳手段に導く合成分離手段を更に含むことを特徴とする。
また、請求項に記載の発明は、請求項に記載の光画像計測装置であって、前記取得手段は、前記合成分離手段により前記信号光に合成された前記補助信号光を、前記信号光の進行方向に対して直交する方向に走査する走査手段を更に含み、前記計測深度算出手段は、前記走査手段による前記補助信号光の走査に基づいて前記補助検出手段により出力される複数の前記補助検出信号に基づいて前記計測深度を算出する、ことを特徴とする。
また、請求項に記載の発明は、請求項に記載の光画像計測装置であって、前記走査手段は、前記変更手段により前記参照光の光路長が変更される度毎に、前記補助信号光を複数の位置に走査させ、前記計測深度算出手段は、前記複数の位置のそれぞれに走査された前記補助信号光に基づく前記補助検出信号に基づいて、前記複数の位置のそれぞれに対応する複数の前記計測深度を算出し、前記算出された前記複数の計測深度に基づいて前記断層画像の前記計測深度を算出する、ことを特徴とする。
また、請求項に記載の発明は、請求項又は請求項に記載の光画像計測装置であって、前記走査手段は、前記補助信号光を反射する反射鏡と、前記反射鏡の反射面の向きを変更する駆動部とを含むガルバノミラーであることを特徴とする。
また、請求項に記載の発明は、請求項ないし請求項のいずれか一項に記載の光画像計測装置であって、前記合成分離手段は、前記信号光の光路に斜設された波長フィルタあるいはビームスプリッタであることを特徴とする。
また、請求項に記載の発明は、請求項ないし請求項のいずれか一項に記載の光画像計測装置であって、前記波長成分分離手段は、回折格子であることを特徴とする。
また、請求項に記載の発明は、請求項ないし請求項のいずれか一項に記載の光画像計測装置であって、前記補助検出手段は、1次元光センサアレイであることを特徴とする。
また、請求項10に記載の発明は、請求項ないし請求項のいずれか一項に記載の光画像計測装置であって、前記計測深度算出手段は、前記補助検出信号に基づいて前記補助干渉光の波長スペクトルを求める波長スペクトル取得手段と、前記求められた前記波長スペクトルをフーリエ変換して、前記計測深度に応じた前記補助干渉光の強度分布を示す計測深度情報を求める計測深度情報取得手段と、を含み、前記画像形成手段により形成される前記断層画像のそれぞれについて前記求められた前記計測深度情報に基づいて前記断層画像のそれぞれの前記計測深度を算出する、ことを特徴とする。
また、請求項11に記載の発明は、請求項10に記載の光画像計測装置であって、前記計測深度算出手段は、前記画像形成手段により連続して形成される一対の前記断層画像のそれぞれについて前記計測深度情報取得手段により求められた前記計測深度情報を比較して、前記一対の前記断層画像の間の計測間隔を算出することにより、前記計測深度を求めることを特徴とする。
また、請求項12に記載の発明は、低コヒーレント光源から出力された光ビームを被測定物体に向かう信号光と参照物体に向かう参照光とに分割し、前記信号光の周波数と前記参照光の周波数とを相対的にシフトさせ、前記被測定物体を経由した前記信号光と前記参照物体を経由した前記参照光とを重畳させて前記周波数のシフト量に応じたビート周波数の干渉光を生成し、前記生成された干渉光を受光して検出信号を出力し、前記出力された検出信号に基づいて前記信号光の進行方向に直交する前記被測定物体の断層画像を形成する光画像計測方法であって、前記形成される前記被測定物体の前記断層画像について、前記信号光の進行方向における計測深度を求めるステップと、前記断層画像と、前記断層画像について取得された前記計測深度の情報とを関連付けて記憶するステップと、前記断層画像と関連付けられた計測深度の情報に基づいて、前記断層画像を前記計測深度方向に配列させるステップと、を有し、前記計測深度を求めるステップは、低コヒーレントな補助光ビームを出力するステップと、前記出力された補助光ビームを、前記信号光とともに前記被測定物体を経由する補助信号光と、補助参照物体を経由する補助参照光とに分割するステップと、前記被測定物体を経由した前記補助信号光と、前記補助参照物体を経由した前記補助参照光とを重畳して補助干渉光を生成するステップと、前記生成された前記補助干渉光から複数の波長成分を分離するステップと、前記分離された前記複数の波長成分を受光して補助検出信号を出力するステップと、前記検出された前記補助検出信号に基づいて前記断層画像の前記計測深度を算出するステップと、を含むことを特徴とする。
また、請求項13に記載の発明は、請求項12に記載の光画像計測方法であって、前記計測深度を求めるステップは、前記補助信号光を、前記信号光の進行方向に対して直交する方向に走査するステップを更に含み、前記計測深度を算出するステップは、前記補助信号光の走査に基づいて前記出力される複数の前記補助検出信号に基づいて前記計測深度を算出する、ことを特徴とする。
また、請求項14に記載の発明は、請求項12又は請求項13のいずれか一項に記載の光画像計測方法であって、前記計測深度を算出するステップは、前記補助検出信号に基づいて前記補助干渉光の波長スペクトルを求めるステップと、前記求められた前記波長スペクトルをフーリエ変換して、前記計測深度に応じた前記補助干渉光の強度分布を表す計測深度情報を求めるステップと、を含み、前記形成される前記断層画像のそれぞれについて前記求められた前記計測深度情報に基づいて前記断層画像のそれぞれの前記計測深度を算出する、ことを特徴とする。
本発明によれば、被測定物体の断層画像について信号光の進行方向における計測深度を求めることができるので、計測中に被測定物体が動いたり移動してしまったような場合であっても、その求めた計測深度に基づいて高精度の3次元画像等を形成することができる。
特に、請求項2に記載の本発明によれば、参照光の光路長を変更させることにより得られる複数の断層画像を、これら複数の断層画像のそれぞれについて求められた計測深度に基づいて計測深度方向に配列させるように構成されているので、高精度の3次元画像等を容易に取得できる。
また、請求項あるいは請求項13に記載の本発明によれば、信号光の進行方向に対して直交する方向に補助信号光を走査するとともに、この補助信号光の走査に基づく複数の補助検出信号に基づいて計測深度を算出するように構成されているので、より高精度の3次元画像等を形成することが可能となる。
本発明に係る光画像計測装置及び光画像計測方法の好適な実施形態の一例について、図面を参照しながら詳細に説明する。
〈第1の実施形態〉
[装置構成]
本発明に係る光画像計測装置の第1の実施形態の構成を図1、2を参照して説明する。図1は本実施形態の光画像計測装置の光学系の概略構成を表し、図2はその制御系の構成を表している。本実施形態の光画像計測装置は、例えば医療や工業分野における被測定物体の断層画像や表面画像、更にはその3次元画像の形成等に利用される装置である。被測定物体は、例えば医療分野においては人眼等の散乱媒質からなる物体である。
〔光学系の構成〕
まず、図1を参照して本実施形態の光画像計測装置の光学系の構成について説明する。同図に示す光画像計測装置1は、低コヒーレントな光ビームを出力する広帯域光源2と、この光ビームの偏光特性を直線偏光に変換する偏光板3と、光ビームを平行光束とするとともにそのビーム径を拡大するレンズ4、5と、光ビームを信号光Sと参照光Rとに分割するとともに、それらを重畳して干渉光Lを生成するハーフミラー6と、参照光Rの偏光特性を直線偏光から円偏光に変換する波長板7と、参照光Rの周波数をシフトさせる周波数シフタ8と、参照光Rの進行方向に対して直交する反射面により参照光Rを全反射する参照鏡9と、参照鏡9の反射面の背面に設けられたピエゾ素子9Aとを含んでいる。
広帯域光源2は、本発明にいう「低コヒーレント光源」に相当し、SLDやLED(発光ダイオード)等により構成される。なお、市販の近赤外域SLDのコヒーレント長は30μm程度、LEDの場合には10μm程度である。
図1中に示すxyz座標系は、広帯域光源2から出力された光ビームの進行方向をz方向とし、それに直交する光ビームの振動面をx−y平面として定義している。x方向、y方向は、光ビームの電場(電界)成分の振動面、磁場(磁界)成分の振動面に一致するように定義される。また、z方向は、被測定物体Oに向かう信号光Sの進行方向として定義されるとともに、被測定物体Oの計測深度方向としても定義される。
偏光板3は、広帯域光源3からの光ビームの所定方向の振動成分を透過させる偏光変換素子である。本実施形態における偏光板3は、上記xyz座標系のx軸及びy軸に対してそれぞれ45°をなす角度方向の振動成分を透過させるように構成される。それにより、偏光板3を透過した光ビームは、角度45°の直線偏光を有するものとなる。したがって、当該光ビームのx軸方向及びy軸方向における偏光成分は、それぞれ等しい振幅を有している。換言すれば、当該光ビームのP偏光成分とS偏光成分とは、それぞれ等しい振幅を有する。
ハーフミラー6は、平行光束とされた直線偏光の光ビームを、被測定物体Oに向かう信号光Sと参照鏡9に向かう参照光Rとに分割するように作用し、本発明の「分割手段」を構成する。ハーフミラー6は、光ビームの一部(半分)を反射して参照光Rを形成するとともに、その残りを透過させて信号光Sを形成する。
また、ハーフミラー6は、被測定物体Oを経由した信号光Sの一部を反射するとともに参照鏡9を経由した参照光Rの一部を透過させることにより、信号光Sと参照光Rとを重畳させて干渉光Lを生成する、本発明の「重畳手段」を構成している。なお、信号光Sの光路上に斜設された符号45の光学素子は、信号光Sの波長を透過する特性を有する波長フィルタ45である(詳細は後述する)。
本実施形態においては、反射体としての被測定物体O及び参照鏡9と、ハーフミラー6とによって形成されるマイケルソン型の干渉計を用いていることから、分割手段と重畳手段とを同一のハーフミラー6(の異なる反射面)により構成している。一方、マッハツェンダー型などの他の干渉計を採用する場合には、分割手段と重畳手段とは、それぞれ別々の光学素子によって構成されていてもよい。また、分割手段及び重畳手段としては、出力された光ビームや信号光Sや参照光Rの偏光特性に影響を与えない無偏光型の任意のビームスプリッタを適用することが望ましい。
波長板7は、直線偏光の参照光Rの偏光特性を変換する偏光変換素子である。本実施形態における波長板7としては、1/8波長板が用いられる。それにより、参照光Rには、波長板7を通過するときに、そのP偏光成分とS偏光成分との間に位相差π/4が与えられる。参照光Rは、ハーフミラー6から参照鏡9に向かうときと、参照鏡9に反射されてハーフミラー6に再入射するときに、それぞれ当該位相差を与えられるので、結果として位相差π/2が付与される。したがって、45°の直線偏光を有する参照光Rに対して1/4波長板と同様に作用することから、ハーフミラー6に再入射される参照光Rは円偏光に変換されることとなる。なお、上述のようにマッハツェンダー型などの他の干渉計を用いる場合には、1/4波長板を適用することができる。
周波数シフタ8は、本発明にいう「周波数シフト手段」を構成し、参照鏡9に反射される前後の参照光Rにそれぞれ周波数シフトを与える。この周波数シフタ8は、例えば電気光学変調器や音響光学変調器などにより構成される。なお、後述するように、本発明に係る光画像計測装置として周波数シフタ8を含まない構成を採用することも可能である。その場合、参照鏡9を移動させることによって参照光Rの周波数をシフトさせる。また、本実施形態の周波数シフタ8は参照光Rの光路上に配置されているが、信号光Sの光路上に配置させるようにしてもよい。すなわち、周波数シフタは、干渉光Lのビートを形成するためのものであるから、信号光Sと参照光Rとの間に相対的に周波数差を形成できるように構成されていれば十分である。
参照鏡9は、本発明の「参照物体」を構成し、参照光Rの光路方向に移動されることにより、被測定物体Oの様々な深度(z座標)による信号光Sの反射光を抽出する。より具体的に説明すると、広帯域光源2からの光ビームは低コヒーレント光であるから、参照光Rとほぼ等距離を経由した信号光Sのみが干渉光Lの生成に寄与する。つまり、ハーフミラー6に対して参照鏡9とほぼ等距離にある被測定物体Oのz座標の領域からの反射光のみが参照光Rと干渉して干渉光Lを生成する。したがって、参照鏡9の位置を移動(z−スキャン)させることにより、被測定物体Oの様々なz座標の領域からの反射光を逐次抽出するようになっている。
また、参照鏡9は、ピエゾ素子9Aによって参照光Rの光路方向に連続的に移動されて参照光Rの周波数をシフトさせるように作用する。このような参照鏡9の移動によって付与される周波数シフトをドップラー周波数シフトと呼ぶことがある。このとき、参照鏡9とピエゾ素子9Aとは、本発明の「周波数シフト手段」を構成している。また、参照鏡9とピエゾ素子9Aは、参照光Rの光路長を変更する本発明の「変更手段」を構成するものである。
なお、参照鏡9とピエゾ素子9Aを周波数シフト用に使用せず、周波数シフタ8のみによって周波数シフトを付与する構成を採用することもできる。また、周波数シフタ8を設けずに、参照鏡9とピエゾ素子9Aによってのみ周波数シフトを付与するように構成することも可能である。
本実施形態の光画像計測装置1には、更に、重畳手段としてのハーフミラー6により生成された干渉光Lを結像させるための結像用レンズ群8と、干渉光Lから異なる複数の偏光成分を分離する偏光ビームスプリッタ11と、分離された干渉光Lの各偏光成分の光路上に設けられたCCD(カメラ)21、22と、CCD21、22のそれぞれによる検出結果を処理する信号処理部20とを含んでいる。
偏光ビームスプリッタ11は、干渉光Lから異なる複数の偏光成分を分離するように作用する。より具体的には、偏光ビームスプリッタ11は、干渉光LのS偏光成分L1を反射してCCD21に受光させるとともに、P偏光成分L2を透過させてCCD22に受光させるように作用する。ここで、干渉光LのS偏光成分L1及びP偏光成分L2は、互いに等しい振幅(つまり最大強度)を有している。
CCD21、22は、本発明にいう「検出手段」を構成するものであり、干渉光検出用の蓄積型の2次元光センサアレイである。CCD21は、偏光ビームスプリッタ11により反射された干渉光LのS偏光成分L1を受光し、それを光電変換して検出信号を生成して信号処理部20に出力する。同様に、CCD22は、偏光ビームスプリッタ11を透過した干渉光LのP偏光成分L2を受光し、それを光電変換して検出信号を生成して信号処理部20に出力する。各CCD21、22から出力される検出信号はヘテロダイン信号である。
信号処理部20は、CCD21、22からそれぞれ出力される検出信号に基づいて後述する演算処理や、当該演算処理の結果を解析して被測定物体Oの2次元断層画像を形成する処理を行う。信号処理部20により形成されるこの2次元断層画像は、信号光Sの光路が参照光Rのそれとほぼ等しくなる被測定物体Oの深度(z座標)の領域におけるx−y断層画像である。このように、光画像計測装置1によれば、被測定物体Oの任意の深さ領域のx−y断層画像を一度の計測で(つまりx−y方向のスキャンを行うことなく)取得することができる。
また、詳細については後述するが、信号処理部20は、様々な深度のx−y断層画像に基づいて、被測定物体Oの3次元画像や、計測深度方向の断層画像(x−z断層画像、y−z断層画像)、あるいは或る軸に斜交する方向の断層画像などの各種画像を形成する。信号処理部20によって形成された画像は、モニタ装置等の表示装置(図示省略)により表示されるようになっている。
以上のような処理を実行する信号処理部20は、例えば、所定のプログラムを格納したROMやハードディスクドライブ等の記憶装置と、当該プログラムを実行するCPU等の演算制御装置とを有するコンピュータ等によって構成される。信号処理部20の構成や作用の詳細については後述する。
更に、本実施形態の光画像計測装置1は、参照光Rに付与される周波数シフトをモニタして、広帯域光源2からの光ビームを強度変調するための構成として、光源31、ビームスプリッタ32、反射鏡33、フォトディテクタ(PD)34及びパルス駆動器35を備えている。
光源31は、広帯域光源2より長いコヒーレント長を持つレーザ光を発するレーザダイオード等によって構成される。ビームスプリッタ32は、光源31からのレーザ光を、周波数シフタ8及び参照鏡9を経由する第1のレーザ光(反射光)と、固定配置された反射鏡33を経由する第2のレーザ光(透過光)とに分割するとともに、周波数シフタ8等により周波数シフトを受けた第1のレーザ光と反射鏡33にて反射された第2のレーザ光とを重畳して干渉光を生成する。
フォトディテクタ34は、ビームスプリッタ32により生成された干渉光を検出し、当該干渉光と等しい周波数の電気信号を出力する。また、パルス駆動器35は、フォトディテクタ34から出力された電気信号に同期した周波数(例えば等しい周波数)のパルス信号を生成して広帯域光源2に出力する。
広帯域光源2は、パルス駆動器34から出力されたパルス信号によって駆動され、当該パルス信号と同期した周波数のパルス状の光ビームを出力する。このとき、広帯域光源2からの光ビームは、例えば、50%dutyの矩形列のパルス光として出力される。
更に、光画像計測装置1には、被測定物体Oの複数のx−y断層画像の計測深度(z座標)を取得するために用いられる計測深度プロファイルを作成するための構成として、広帯域光源41、ハーフミラー42、固定鏡43、ガルバノミラー44、波長フィルタ45、回折格子46及び1次元光センサアレイ47が設けられている。なお、詳細については後述するが、この計測深度プロファイルは、本発明にいう「計測深度情報」を構成するものである。
広帯域光源41は、計測用の広帯域光源2と異なる波長の光ビーム(「補助光ビーム」と呼ぶ)を出力する本発明の「補助光源」である。広帯域光源41から出力された補助光ビームの一部(「補助信号光」と呼ぶ)は、本発明の「補助分割手段」として作用するハーフミラー42により反射されてガルバノミラー44に向かって進行し、また、補助光ビームの残りの部分(「補助参照光」と呼ぶ)は、ハーフミラー42を透過して固定鏡43に向かう。
ガルバノミラー44は、本発明にいう「走査手段」を構成し、反射鏡と、この反射鏡を駆動して反射面の向きを変更する駆動部とを備えている(図2を参照)。また、波長フィルタ45は、広帯域光源2からの光ビーム(したがって信号光S)の波長を透過させ、広帯域光源41からの補助光ビーム(したがって補助信号光)の波長を反射する特性を有する。
それにより、ハーフミラー42による反射光からなる補助信号光は、ガルバノミラー44によって反射され、波長フィルタ45により更に反射されて信号光Sとともに被測定物体Oに照射される。被測定物体Oに照射された補助信号光は、信号光Sとともに様々な深度(z座標)の領域にて反射されて被測定物体Oから出射し、波長フィルタ45によって信号光Sから分離されるとともに、ガルバノミラー44により反射されてハーフミラー42に入射する。なお、波長フィルタ45は、本発明にいう「合成分離手段」を構成している。
被測定物体Oを経由してハーフミラー42に入射した補助信号光は、固定鏡43により反射された補助参照光と重畳されて干渉光(「補助干渉光」と呼ぶ)を形成する。ハーフミラー42は、本発明の「補助重畳手段」を構成している。
補助干渉光は、回折格子46により各波長成分に分離され、1次元光センサアレイ47によって受光される。1次元光センサアレイ406を構成する各光センサは、検出した波長成分の光強度の検出信号(「補助検出信号」と呼ぶ)を信号処理部20に出力する。ここで、回折格子46は、本発明の「波長成分分離手段」を構成し、1次元光センサアレイ47は、本発明の「補助検出手段」を構成している。
信号処理部20は、1次元光センサアレイ47から出力された各波長成分に対応する補助検出信号に基づいて、補助干渉光に関する波長−光強度の関係、換言すれば補助干渉光の光強度分布(波長スペクトル)を求める(図11(A)を参照)。
更に、信号処理部20は、求めた補助干渉光の波長スペクトルをフーリエ変換することにより、被測定物体Oのz座標(計測深度)を変数とした補助干渉光の強度分布を求める(図11(B)を参照)。この強度分布は、前述の計測深度プロファイル(計測深度情報)に相当する。そして、この強度分布を参照して、各x−y断層画像の計測深度を算出する。以上のような計測深度取得処理は、MHz以上の速さ、すなわちμsec以下の時間で実行される。
なお、ガルバノミラー44により補助信号光を走査しつつ計測を行う場合、各x−y断層画像について、複数の走査位置に対応する複数の補助検出信号が信号処理部20に入力される。その場合、信号処理部20は、この複数の補助検出信号に基づいて当該x−y断層画像の計測深度を算出する。
信号処理部20は、取得した複数のx−y断層画像と、各x−y断層画像について算出される計測深度に基づいて複数のx−y断層画像を配列させるとともに、その配列された複数のx−y断層画像に対して3次元化処理(補完処理)等を施して被測定物体Oの3次元画像等を形成して表示装置に表示させる。以上に説明した信号処理部20による処理の詳細については後述する。
なお、本実施形態において、被測定物体Oのx−y断層画像の計測深度を求めるように作用する本発明の「取得手段」は、広帯域光源41、ハーフミラー42、固定鏡43、ガルバノミラー44、波長フィルタ45、回折格子46、1次元光センサアレイ47及び信号処理部20を含んで構成されている。
また、信号処理部20は、CCD21、22からの検出信号に基づいて被測定物体Oのx−y断層画像を形成する本発明の「画像形成手段」を含むとともに、取得された計測深度に基づいて複数のx−y断層画像を配列させる本発明の「画像処理手段」を含んでいる。
〔制御系の構成〕
次に、図2を参照して、光画像計測装置1の制御系の構成について説明する。光画像計測装置1の制御系は、信号処理部20、CCD21、22、ピエゾ素子9A(及び参照鏡9)、広帯域光源41、ガルバノミラー44、1次元光センサアレイ47、及び表示装置60を含んで構成される。
ガルバノミラー44は、上述のように、広帯域光源41からの光ビームに基づく干渉光の光路上に配置された反射鏡44Aと、この反射鏡44Aを駆動してその反射面の向きを変更する駆動部44Bとを有している。
信号処理部20は、前述のようなコンピュータにより構成されており、制御部50、画像形成部51、走査制御部52、計測深度算出部53、情報記憶部56、画像処理部57を含んでいる。
制御部50は、信号処理部20による各種の演算処理や画像処理や制御処理等を制御するものであり、CPU等によって構成される。
(画像形成部)
画像形成部51は、本発明にいう「画像形成手段」を構成し、CCD21、22から出力された検出信号に基づいて被測定物体Oのx−y断層画像を形成する処理を行う。この画像形成部51による画像形成処理は、従来と同様の手順にて実行される。画像形成部51には、参照鏡9の移動に伴う計測深度に対応する検出信号がCCD21、22から順次入力される。画像形成部51は、各検出信号に対して画像形成処理を実行してx−y断層画像を順次形成する。ここでx−y断層画像の計測枚数をNとする。画像形成部51により逐次形成されるN個のx−y断層画像G1〜GNは、制御部50により情報記憶部56に記憶される。
(走査制御部)
走査制御部52は、ガルバノミラー44の駆動部44Bに制御信号を送信して、反射鏡44Aの反射面の向きを変更して補助信号光を信号光Sの進行方向に対して直交する方向に走査させる。図3は、補助信号光の走査の一例を表す。同図は、図1において波長フィルタ45側から被測定物体Oを見たときの状態を表している。同図には、被測定物体O上に円を描くように補助信号光を走査させたときの補助信号光の軌跡が示されている。すなわち、補助信号光は、例えば垂直上方位置を走査の起点T1とし、時計回りに一定間隔でT2、T3、T4、・・・・・、T(M−1)、TMと走査される。このとき、ガルバノミラー44の反射鏡44Aは、その反射面の法線が一定速度で回転するように駆動され、かつ、広帯域光源41は、反射鏡44Aが一回転する間にM回の補助光ビームを一定間隔で出力するように制御される。それにより、1次元ラインセンサ47からは、一定間隔で補助干渉光の補助検出信号が出力される。
なお、連続的な補助光ビームを出力する広帯域光源の前に液晶シャッタ等を配置して一定間隔で開閉させることにより、補助光ビームを断続的に出力させるようにしてもよい。また、広帯域光源41から連続的な補助光ビームを出力して、補助信号光を連続的に走査させるようにしてもよい。
(計測深度算出部)
計測深度算出部53は、本発明にいう「計測深度算出手段」を構成し、1次元光センサアレイ47から出力される補助検出信号に基づいて各x−y断層画像G1〜GNの計測深度を算出する処理を行う。
より詳しく説明すると、計測深度算出部53は、補助検出信号に基づいて補助干渉光の波長スペクトルを求めるスペクトル取得部54と、この波長スペクトルをフーリエ変換して、計測深度に応じた補助干渉光の強度分布、すなわち計測深度プロファイル(計測深度情報)を求めるプロファイル取得部55とを有している。この計測深度プロファイルは、例えば図3に示した補助信号光の走査態様を適用する場合には、各x−y断層画像G1〜GNについてそれぞれM個ずつ作成される。なお、スペクトル取得部54は、本発明にいう「波長スペクトル取得手段」を構成し、プロファイル取得部55は、本発明の「計測深度情報取得手段」を構成する。
更に、計測深度算出部53は、各x−y断層画像G1〜GNに関するM個の計測深度プロファイルを比較して、連続して形成されたx−y断層画像Gi、G(i+1)の間の計測間隔を算出することにより、各x−y断層画像G1〜GNの計測深度を求める(i=1〜N−1)。求めた計測深度は、制御部50により、計測深度算出結果P1〜PNとして情報記憶部56に記憶される。このとき、各計測深度算出結果Pjは、対応するx−y断層画像Gjと関連付けられて記憶される(j=1〜N)。
この計測深度算出部53による処理は、例えば次のようにして行う。ここで、第i枚目のx−y断層画像Gi(i=1〜N)が計測された深度における補助信号光の各走査位置Tk(k=1〜M)をTk(i)で表すこととする。また、図4は、x−y断層画像G1〜GNを取得するときのz方向への走査態様の一例を表し、連続する画像Gi、G(i+1)の計測間隔をd=一定とする(被測定物体Oは静止状態と仮定する)。なお、画像の計測間隔は、参照鏡9のz方向への移動速度と、広帯域光源2による光ビームの出力周期とによって設定される。この仮定の計測間隔dは、制御部50により把握され、計測深度算出部53に送られる。なお、参照鏡9を連続的に移動させる代わりに、段階的に参照鏡9を移動させて各画像Giの計測を逐次行うようにしてもよい。
さて、図4に示す第i枚目のx−y断層画像Giの計測を行うとき、走査制御部52は、ガルバノミラー44を制御して、補助信号光を図3に示す走査位置T1〜TMに走査させる。各走査位置Tk(k=1〜M)の様々な深さ領域にて反射された補助信号光に基づく補助干渉光は、回折格子46により波長成分にλ1〜λnに分離されて1次元光センサアレイ47により検出される。計測深度算出部53のスペクトル取得部54は、1次元光センサアレイ47からの補助検出信号に基づいて図11(A)に示すような波長スペクトルを求める。
次に、プロファイル取得部55は、この波長スペクトルをフーリエ変換して図5(A)に示すような計測深度プロファイルを取得する。なお、本実施形態の光画像計測装置1では、被測定物体Oが静止していると仮定した場合、各走査位置における計測深度プロファイルは、画像Giに依らずに常に同一となる。
図5(B)は、画像G(i+1)の計測時に、被測定物体Oが移動してしまった場合に取得される計測深度プロファイルの態様を表している。図5(B)に示す計測深度プロファイルは、図5(A)に示すプロファイルと同じ波形が計測深度方向に変位したものとなっている。ここで、各計測位置Tkについて、画像Giと画像G(i+1)におけるプロファイルの変位量をΔdk(i+1)で表す(i=1〜N−1)。なお、画像Gi、G(i+1)の計測の間に被測定物体Oが静止していたときはΔdk(i+1)=0であり、+z方向に移動するとΔdk(i+1)>0、−z方向に移動するとΔdk(i+1)<0となる。
計測深度算出部53は、連続して計測される各画像Gi、G(i+1)の計測深度プロファイルを比較して、その変位量Δdk(i+1)を算出する。当該算出処理は、例えば、各プロファイルをパターン認識したり、平行移動させることにより行う。また、各プロファイルの例えばピーク値の計測深度方向への変位量を算出するようにしてもよい。
更に、計測深度算出部53は、この変位量dk(i+1)を用いて、画像Gi、G(i+1)の計測間隔d(前述のように一定と仮定されている)を変更(補正)する。すなわち、変位量Δdk(i+1)=0である場合、画像Giと画像G(i+1)との「実際の」計測間隔をdとし、変位量dk(i+1)≠0である場合には、「実際の」計測間隔はd−Δdk(i+1)とする。求められた実際の計測間隔は、画像G(i+1)と関連付けられて計測深度算出結果P(i+1)として情報記憶部56に記憶される。それにより、各画像G(i+1)について、走査位置T1〜TMに対応するΔd1(i+1)〜ΔdM(i+1)に基づくM個の計測間隔が算出され、計測深度算出結果P(i+1)として記憶される。
なお、例えば第1枚目の画像G1の計測深度を基準(z=z1=0)するとともに、第2枚目の画像G2の変位量Δdk(2)に基づいて画像G2のz座標z2を算出し、同様にして第i枚目の画像Giのz座標ziを算出することにより、各画像G1〜GNの計測深度(z座標)z1〜zNを求め、各画像G1〜GNと関連付けて計測深度算出結果P1〜PNとして記憶することもできる。
また、各走査位置T1〜TMについて得られる変位量Δd1(i+1)〜ΔdM(i+1)の平均値Δd(i+1)を算出し、この平均値Δd(i+1)を用いて画像Gi、G(i+1)の計測間隔を求めるようにしてもよい。
(情報記憶部)
情報記憶部56は、イメージメモリやハードディスクドライブ等の記憶装置により構成される。情報記憶部56に記憶される情報にはディレクトリが割り振られ、ディレクトリ同士を関連付けることにより、記憶情報を関連付けるようになっている。特に、x−y断層画像Giと計測深度算出結果Piとは、このようにして互いに関連付けられて記憶される(i=1〜N)。なお、情報記憶部56に対する記憶処理及び関連付け処理は、制御部50により実行される。
(画像処理部)
画像処理部57は、本発明にいう「画像処理手段」を構成するもので、計測深度算出部53により求められたx−y断層画像G1〜GNのそれぞれの計測深度算出結果P1〜PNに基づいて、画像G1〜GNを計測深度方向(z方向)に配列させる処理を行う。更に、画像処理部57は、計測深度方向に配列された画像G1〜GNに対して画像処理(3次元化処理等の補完処理など)を施すことにより、被測定物体Oの3次元画像や、計測深度方向の断層画像(x−z断層画像、y−z断層画像)、あるいはx、y、z軸に斜交する方向の断層画像を形成する処理を行う。
より具体的に説明すると、画像処理部57は、x−y断層画像G1〜GNを情報記憶部56から読み出すとともに、計測深度算出結果P1〜PNに含まれる画像Giと画像G(i+1)(i=1〜N−1)との間の計測間隔Δd1(i+1)〜ΔdM(i+1)を参照して、画像G1〜GNをz方向に配列させる。なお、計測深度算出結果Pi(i=1〜N)として各画像Giのz座標ziが記憶されているときには、各x−y断層画像Giをz座標ziに配置させる。更に、画像処理部57は、計測間隔あるいはz座標に応じて配列されたx−y断層画像G1〜GNについて、連続する画像Gi、G(i+1)の間の画像を補完することにより被測定物体Oの3次元画像を形成する。
また、被測定物体Oの(任意のy座標y0における)x−z断層画像を形成する場合、画像処理部57は、各x−y断層画像Gi(i=1〜N)のy=Y0における部分画像を、計測深度算出結果Piを参照してz方向に配列させるとともに、補完処理を施してx−z断層画像を形成する。被測定物体Oの(任意のx座標x0における)y−z断層画像を形成する場合も同様である。
また、x、y、z軸に斜交する方向における被測定物体Oの断層画像を形成する場合については、例えば、その斜交面(断層面)と交差する3次元画像の部分のみを抽出することにより画像化できる。なお、x−z断層画像やy−z断層画像を形成する場合についても、同様に、3次元画像のスライス画像を形成するようにしてもよい。
画像処理部57により形成された3次元画像等は、制御部50によって表示装置60に表示される。また、形成された3次元画像等を情報記憶部56等に格納するようにしてもよい。
〔測定形態〕
続いて、本実施形態の光画像計測装置1により実行される干渉光Lの信号強度や位相の空間分布の測定形態、すなわちヘテロダイン信号の信号強度やその位相情報の測定形態について説明する。以下に詳述する信号処理は、図1に示した信号処理部20によって実行されるものである。
本実施形態の光画像計測装置1は、偏光特性の異なる信号光Sと参照光Rを形成し、それらの干渉光Lをヘテロダイン信号として検出することにより、被測定物体Oの表面画像や断層画像を得ることを特徴としている。
(測定原理)
まず、光画像計測装置1による測定形態の基本原理を説明する。広帯域光源2から出力された光ビームは、偏光板3により上記x軸に対して45°をなす角度方向の直線偏光に変換され、レンズ4、5によってビーム径を拡大され、かつ、平行光束とされてハーフミラー6に入射して信号光Sと参照光Rとに2分される。
信号光Sは、散乱媒質からなる被測定物体Oに入射し、その表面や様々な深さの断層面にて反射される。被測定物体Oからの反射光波の一部はハーフミラー6により反射されて結像用レンズ群10に伝送される。
一方、参照光Rは波長板7を通過して参照鏡9へと伝送される。このとき、参照鏡9は、ピエゾ素子9Aによって参照光Rの光路方向に駆動(z−スキャン)されている。また、参照光Rは、周波数シフタ8によって所定量の周波数シフトを受ける。参照鏡9からの反射光波は、参照鏡9のz-スキャンに伴うドップラー周波数シフト、更には周波数シフタ8による周波数シフトを受け、波長板7を通過する。ここで、参照光Rの偏光特性は角度45°の直線偏光であり、波長板7は1/8波長板であることから、波長板7を2回通過した参照光Rの偏光特性は円偏光に変換されることとなる。円偏光とされた参照光Rの一部はハーフミラー6を透過して結像用レンズ群10に伝送される。
このとき、ハーフミラー6は、被測定物体Oにて反射された直線偏光の信号光Sと、周波数がシフトされ、円偏光に変換された参照光Rとを重畳して干渉光Lを生成し、この干渉光Lが結像用レンズ群10に伝送されることとなる。干渉光Lは、結像用レンズ群10を経由して偏光ビームスプリッタ11に伝搬される。
偏光ビームスプリッタ11は、干渉光LのS偏光成分L1を反射し、P偏光成分L2を透過するように作用する。干渉光LのS偏光成分L1はCCD21により検出され、P偏光成分L2はCCD22によって検出される。ここで、干渉光LのS偏光成分L1は、信号光SのS偏光成分Essと、参照光RのS偏光成分Ersとを含んでおり、干渉光LのP偏光成分L2は、信号光SのP偏光成分Espと、参照光RのP偏光成分Erpとを含んでいる。信号光SのS偏光成分Ess及びP偏光成分Espと、参照光RのS偏光成分Ers及びP偏光成分Erpとは、次式のように表される。
Figure 0004566685
ここで、fは広帯域光源2から出力される光ビームの周波数を表し、fは周波数シフトを表し、φは信号光Sの初期位相を、φ′は参照光Rの初期位相をそれぞれ表す。更に、信号光Sと参照光Rとの初期位相の差をΔφ(=φ−φ′)と表すこととする。[数2]に示す式(2)〜(5)を参照すると、干渉光LのS偏光成分L1とP偏光成分L2とは、CCD21、22により、それぞれ次式のようなヘテロダイン信号i、iとして検出される。
Figure 0004566685
式(6)、(7)を比較すると、各式の第3項の交流信号は、同位相のcos関数とsin関数であることから90°の位相差があることが分かる。本発明は、この特徴を利用するとともに、周期的に強度変調された光ビームを測定光として使用することにより、シャッタによるサンプリング処理を用いない光ヘテロダイン検出を実現可能とし、それにより干渉光Lの信号強度および位相の空間分布を測定するものである。従来の光画像計測技術においては、単一の干渉光を位相の異なる複数の関数でサンプリングすることによりそのcos成分とsin成分とを検出していたが、本発明においては、参照光Rや信号光Sの偏光特性を変換して位相の異なる複数(ここでは2つ)の干渉光を生成し、それらを別々に検出するように構成されている点が特徴的である。以下に、本発明における測定原理を説明する。
本実施形態では、光源31、ビームスプリッタ32、反射鏡33、フォトディテクタ(PD)34及びパルス駆動器35を用いて、広帯域光源2から周期的に強度変調された光ビームを出力する。
光源31から出力されたレーザ光は、ビームスプリッタ32によって参照鏡9方向の光路(反射光)と反射鏡33方向の光路(透過光)とに分割される。参照光9方向の光路のレーザ光は、周波数シフタ8と参照鏡9を経由するようになっており、それらによる周波数シフトを受けてビームスプリッタ32に再入射する。一方、反射鏡33方向の光路のレーザ光は、反射鏡33による反射光として(周波数をシフトされることなく)ビームスプリッタ32に再入射する。両光路を経由したレーザ光は、ビームスプリッタ33によって重畳され干渉光となってフォトディテクタ34により検出される。
フォトディテクタ34により検出される干渉光は、参照光Rと同様に、周波数シフタ8に寄る周波数シフトと参照鏡9によるドップラー周波数シフトとを受けるので、参照光Rと(ほぼ)同量の周波数シフトを受ける。したがって、当該干渉光は、信号光Sと参照光Rとからなる干渉光Lと(ほぼ)同一のビート周波数を有している。
フォトディテクタ34は、検出した干渉光に対応する電気信号をパルス駆動器35に出力する。この電気信号は、式(1)に示すヘテロダイン信号と同様に直流成分と交流成分とを有しており、その交流成分の周波数は上述のように干渉光Lのビート周波数とほぼ同じである。パルス駆動器35は、フォトディテクタ34からの電気信号を受け、それと同じ周波数を有するパルス信号を広帯域光源2に出力する。広帯域光源2は、パルス駆動器35からのパルス信号によって駆動されて、当該パルス信号と同じ周波数のパルス状の光ビームを出力する。
このように、本実施形態では、参照光Rに付与される周波数シフトのシフト量をモニタし、当該シフト量と(ほぼ)等しい周波数のパルス状の光ビームを用いて被測定物体Oの測定を行うようになっている。広帯域光源2からの光ビームは、上述のように、例えば50%dutyの矩形列のパルス光として出力される。なお、当該光ビームのduty比は50%に限定されるものではなく、パルス形状は矩形列以外(例えば三角形列や台形列等)でもよい。また、出力強度0と100とを切り換えて得られるパルス光に代えて、例えば出力強度を50と100との間で変調して得られる光ビームを適用することもできる。すなわち、ここで重要なのは、光ビームの強度変調の度合ではなく、その強度変調の周波数が干渉光Lのビート周波数とほぼ等しくなるように制御されることである。
次に、図6に示すグラフを参照して、本実施形態の光画像計測装置1における干渉光Lの検出態様について説明する。以下、広帯域光源2から出力される光ビームの強度の変調周波数をfとする。また、前述したように、fは参照光Rに付与される周波数シフト(干渉光Lのビート周波数)を表し、光ビームの変調周波数fは周波数シフトfと等しいか、それに近い値とされている。
図6(A)は、周波数fで強度変調されて広帯域光源2から出力される光ビームの時間波形を表す。図6(B)は、光ビームが連続光であり、よって参照光Rと信号光Sがともに連続光である場合における干渉光LのS偏光成分L1(ビート周波数f)の時間波形を表す。図6(C)は、参照光Rと信号光Sがともに連続光である場合における干渉光LのP偏光成分L2の時間波形を表す。ここで、図6(B)、(C)に示すS偏光成分L1とP偏光成分L2との位相差は90°である。
また、図6(D)は、広帯域光源2からの光ビームが図6(A)のように強度変調される場合における干渉光LのS偏光成分L1の時間波形を表す(図6(B)に対応する)。図6(E)は、広帯域光源2からの光ビームが図6(A)のように強度変調される場合における干渉光LのP偏光成分L2の時間波形を表す(図6(C)に対応する)。図(D)、(E)に示すS偏光成分L1とP偏光成分L2とは90°の位相差を有する。
CCD21は、図6(D)に示す時間波形のS偏光成分L1を検出する。広帯域光源2からの光ビームは周波数f、50%dutyの矩形列の光パルスであり、その変調周波数fと干渉光Lのビート周波数fとの差δf=|f−f|が蓄積型光センサであるCCD21の応答周波数に比べて十分小さいときには、CCD21から出力されるS偏光成分L1の検出信号は、検出時間内に蓄積された光電荷量に比例するものとなり、次式のように与えられる。(例えば、M.Akiba、K.P.Chan、N.Tanno、Japanese Journal of Applied Physics、Vol.39、L1194(2000)参照)。
Figure 0004566685
ここで、<・>はCCD21の蓄積効果による時間平均を表し、Kは偏光ビームスプリッタ11の反射率とCCD21の光電変換率を含めた光検出効率、m(t)は広帯域光源2の出力を強度変調する関数(矩形のパルスを表す関数)、またβは測定における初期位相値を表す。式(8)から分かるように、CCD21から出力される検出信号には、信号光Sと参照光Rの強度に関する項(背景光成分)の他に、干渉光LのS偏光成分L1の振幅√(Issrs)及び位相2πδft+βに関する項が含まれている。
同様に、CCD22は、図6(E)に示す時間波形のP偏光成分L2を検出し、次式のような検出信号を出力する。
Figure 0004566685
ここで、Kは偏光ビームスプリッタ11の透過率とCCD22の光電変換率を含めた光検出効率である。
次に、CCD21、22からそれぞれ出力される検出信号(8)、(9)に基づく、干渉光Lの信号強度の算出処理について説明する。
参照光Rは、波長板7により円偏光に変換されているので、そのS偏光成分Ersの強度IrsとP偏光成分Erpの強度Irpとは等しいと考えられる(Irs=Irp=Iと示す)。
一方、信号光Sについて、被測定物体Oからの反射光は入射光の偏光特性に顕著に依存しないと考えられることから、そのS偏光成分Essの強度IssとP偏光成分Espの強度Ispとは等しいか、あるいは近い値であると考えられる(Iss=Isp=Iと示す)。また、信号光Sは被測定物体Oによって散乱、吸収されることから、その強度は一般的に参照光Rより十分に小さい(I<<I)と考えることができる。
また、式(8)及び式(9)の右辺の第1項と第2項は背景光の強度を表し、その値は、事前に若しくは別途に測定することができる。例えば、広帯域光源2により連続光からなる光ビームを出力してCCD21等により検出し、それを1波長分(あるいはその整数倍)だけ積分して第3項(交流成分;位相直交成分)をキャンセルすることにより、背景光の強度(背景光成分)を取得することができる。
取得された背景光成分を各CCD21、22からの検出信号の強度から除算することにより、各検出信号の位相直交成分、すなわち、干渉光LのS偏光成分L1及びP偏光成分L2の位相直交成分S′(t)、S′(t)が算出される(次式を参照)。
Figure 0004566685
これら式(10)、(11)を用いると、干渉光LのS偏光成分L1及びP偏光成分L2の振幅は次式のように表される。
Figure 0004566685
更に、本実施形態の光画像計測装置1は、次のようにして干渉光Lの位相の空間分布を画像化する。
或る測定時間t=tにおいて、干渉光LのS偏光成分L1の位相直交成分S′(t1)がCCD21により検出され、P偏光成分L2の位相直交成分S′(t1)がCCD22により検出されたとすると、これら両位相直交成分の比を取ると次のような信号が得られる。
Figure 0004566685
この式(13)に示す信号Sは、干渉光Lの振幅に依存せず、位相情報のみから構成されていることが分かる。本実施形態では、複数のピクセルが2次元的に配列された受光面を持つCCD21、22によりS偏光成分L1とP偏光成分L2を検出しているので、各ピクセルで検出される信号の位相β(x、y、t)は、次式のように表される(ここで、(x、y)は、各ピクセルの受光面上における座標を表す)。
Figure 0004566685
この式(14)の第2項は、ゼロあるいはほぼゼロの周波数δf(≒0)を有する交流信号の測定時間tにおける瞬時位相値であるから、CCD21、22のピクセルの位置、つまり座標x、yに依らずに均一であると考えられる。したがって、例えばCCD21、22の受光面上の或る特定点(x=x、y=y)に位置するピクセルで検出される位相φ(x、y、t)を基準に、各ピクセルで検出される検出信号の位相差を求めれば、ヘテロダイン信号の位相差の空間分布、すなわち干渉光Lの位相差の空間分布を画像化できる。
他方、干渉光Lの位相情報からその周波数情報を取得することも可能である。2つの測定時間t=t及びt=tにおける干渉光LのS偏光成分L1及びP偏光成分L2の位相をそれぞれβ(x、y、t)及びβ(x、y、t)とすると、干渉光Lのビート周波数fと、広帯域光源2からの光ビームの変調周波数fとの差δfは、次式のように表される。
Figure 0004566685
ここで、光ビームの変調周波数fは既知であるので、式(10)や式(11)からヘテロダイン周波数、つまり干渉光Lのビート周波数fを算出できる。このようなヘテロダイン周波数の測定方法は、例えばヘテロダイン干渉法を用いたドップラー速度計測に有効であると考えられる。
[作用効果]
本実施形態の光画像計測装置1は、以上のような原理によって取得される被測定物体Oの複数のx−y断層画像をz方向(計測深度方向)に配列させるとともに、その配列された複数のx−y断層画像に画像処理を施して被測定物体Oの3次元画像等を形成する。ここで、複数のx−y断層画像は、各画像毎に算出された計測間隔(あるいはz座標)に基づいて配列される。したがって、計測中に被測定物体Oが動いたり移動した場合であっても、各x−y断層画像の計測深度を良好に変更(補正)することができるので、高精度の3次元画像等を形成することができる。
また、補助信号光を走査することにより、各x−y断層画像の複数の位置における計測深度を取得することができるので、被測定物体Oが回転運動を行ったような場合であっても、x−y断層画像をその回転運動による変位を打ち消す方向に回転させて補正を行うことができる。例えば、図3における走査位置T1と走査位置T(M/2+1)(Mは偶数とする)のような、対向位置にある一対の走査位置の計測深度プロファイルの変位を比較することにより、被測定物体Oの回転方向を求めることができ、その求めた回転方向を打ち消す方向にx−y断層画像を回転させて補正することができる。また、x、y、zの各方向において対向する走査位置をそれぞれ比較して、x方向への回転量、y方向への回転量、z方向への回転量を求め、それらを合成することにより被測定物体Oの3次元的な回転を求めるようにしてもよい。このような処理を行えば、被測定物体Oの3次元画像等をより高い精度で形成することが可能となる。
なお、補助信号光の走査は任意の軌跡上を行うことができる。また、本発明に係る光画像計測装置としては、補助信号光を走査させる必要はなく、各x−y断層画像について1点のみを測定すれば十分である。
〔変形例〕
上記の実施形態は、ガルバノミラーを用いて補助信号光を走査するように構成されているが、本発明の走査手段はガルバノミラーに限定されるものではなく、補助信号光を好適に走査できる任意の構成を採用することができる。
また、補助信号光を信号光に合成するとともに、補助信号光を信号光から分離する本発明の合成分離手段として波長フィルタが用いられているが、ビームスプリッタ等を代替的に用いることができる。
更に、補助干渉光を検出する本発明の補助検出手段として1次元光センサアレイが用いられているが、CCD等の2次元光センサアレイ等を代わりに用いてもよい。ただし、形成される3次元画像等の精度を考慮すると、より多数の走査位置における計測深度を取得するために、読み出しスピードの速い1次元光センサアレイを用いることが望ましいと思われる。
また、上記構成では、参照光Rに周波数シフトを付与するための構成として、周波数シフタ8と、参照鏡9及びピエゾ素子9Aとの双方を考慮したが、これらの内の一方のみを備えていてもよい。例えば、周波数シフタ8を有さない光画像計測装置を形成し、参照鏡9のz−スキャンのみによって参照光Rに周波数シフトを付与するように構成しても、同様の測定を実行することができる。また、周波数シフタ8を用いる場合、信号光Sの光路上に設けるようにしてもよい。本発明に係る画像計測においては、重畳時における信号光Sの周波数と参照光Rの周波数とが相対的にシフトされていれば十分だからである。
また、上記構成では、広帯域光源2からの光ビームをまず直線偏光とし、それから信号光Sと参照光Rとに分割するようになっているが、光ビームの分割後に信号光Sと参照光Rとをそれぞれ直線偏光に変換するようにしてもよい。ただし、その場合には、信号光Sと参照光Rの双方の光路上に偏光板を設ける必要があり、上記構成よりも若干複雑な構成となるため、実用上は上記構成の方が好適であると思われる。
また、上記構成では、参照光Rの偏光特性を円偏光に変換するようになっているが、信号光Sの方を円偏光に変換し、直線偏光のままの参照光と重畳させるようにすることも可能である。しかし、上述のように、信号光Sの被測定物体Oによる反射光は参照光Rと比較して微弱であるので、信号光Sの光路上に波長板を配置させると、それを通過するときに信号光Sが弱められてしまう。このように被測定物体Oの情報を含んだ信号光Sの強度を弱めることは測定の感度に悪影響を及ぼす可能性がある。したがって、上記構成のように参照光Rの偏光特性を変換する方が有利といえる。なお、周波数シフタの配置についても同様である。
また、上記構成では、参照光Rの周波数のシフト量をモニタする光源31、ビームスプリッタ32、反射鏡33及びフォトディテクタ34が設けられ、そのモニタ結果を光ビームの強度変調にフィードバックするようになっているが、例えば参照光Rに付与する周波数シフト量が設定されているときなどには、当該シフト量と(ほぼ)等しい周波数のパルス信号を自発的に生成するパルス駆動器35を設けて光ビームの強度変調を制御するようにしてもよい。
また、パルス駆動器35による広帯域光源2のパルス駆動に代えて、連続的な光ビーム(連続光)を発する広帯域光源2と、この連続的な光ビームを周期的に遮断するシャッタとを設けることにより、光ビームの強度を周期的に変調させるようにしてもよい。
〈第2の実施形態〉
本発明に係る光画像計測装置の第2の実施形態について説明する。第1の実施形態は、x−y断層画像の計測深度を取得するために、画像計測用の光源(広帯域光源2)とは別の光源(広帯域光源41)を備える構成とされていたが、本実施形態は、画像計測用の光源に基づく干渉光の一部を検出することにより、x−y断層画像の計測深度を求める構成を有するものである。
図7は、本実施形態の光画像計測装置の構成の一例を表す。同図に示す光画像計測装置100は、第1の実施形態の光画像計測装置1とほぼ同様の構成を備えている。図7において、図1と同じ構成要素は同一の符号で示されている。
本実施形態の光画像計測装置100は、ハーフミラー6により生成された干渉光Lの光路上に斜設され、その干渉光Lの一部を分離するビームスプリッタ71を備えている。このビームスプリッタ71は、本発明にいう「干渉光分離手段」を構成するものである。
ビームスプリッタ71により分離された干渉光Lの一部は、その光路上に斜設された回折格子72(波長成分分離手段)により複数の波長成分に分離されて1次元光センサアレイ73(補助検出手段)によって検出される。1次元光センサアレイ73は、その検出信号(補助検出信号)を信号処理部20に出力する。
信号処理部20は、第1の実施形態と同様の構成を備え(図2を参照)、CCD21、22からの検出信号に基づいて形成される各x−y断層画像の計測深度を計測深度算出部53により算出し、画像処理部57によりx−y断層画像を計測深度に基づいてz方向(計測深度方向)に配列させる。更に、信号処理部20は、その配列された複数のx−y断層画像に対し画像処理を施すことにより被測定物体Oの3次元画像等を形成し、表示装置に表示させる。本実施形態において、本発明の「取得手段」は、ビームスプリッタ71、回折格子72、1次元光センサアレイ73及び信号処理部20(計測深度算出部53)により構成される。
なお、本実施形態においては、計測深度プロファイルの変位は参照光Rの光路長の変化を加味したものとなるので、計測深度プロファイルの変位から参照光Rの光路長の変化分を差し引いた値が被測定物体Oの変位量となる。なお、図2に示す制御系の構成から分かるように、ピエゾ素子9Aによる参照光Rの光路長の変化は、制御部50によって制御されるので、当該光路長の変化分は制御部50によって把握されている。計測深度算出部53(計測深度算出手段)は、制御部50から当該変化分の情報を受けてx−y断層画像の計測深度を算出する。
このような本実施形態の光画像計測装置100によれば、第1の実施形態と同様に、計測中に被測定物体Oが動いたり移動した場合であっても、各x−y断層画像の計測深度を良好に変更(補正)することができるので、高精度の3次元画像等を形成することが可能である。
〈その他の変形例〉
以上に詳述した構成は、本発明に係る光画像計測装置を実施するための構成の一例に過ぎないものである。したがって、本発明の要旨の範囲内において各種の変形を施すことが可能である。
例えば、上記の各実施形態の構成では、断層画像を取得するための計測と断層画像の計測深度を取得するための計測とを同時に実行し、それぞれの計測の結果を対のデータ(図2に示すx−y断層画像G1〜GNと計測深度算出結果P1〜PN)として記憶するように構成されているが、本発明はこれに限定されるものではなく、それぞれの計測を異なるタイミングにて行うように構成することも可能である。また、断層画像の計測中にその計測深度を複数回計測し、その複数回の計測の結果に基づいて被測定物体の動きを推定して補正するように構成することもできる。
また、上記実施形態の光画像計測装置は、断層画像の計測深度を算出するための情報(補助検出信号)を取得する構成と、その取得した情報から計測深度を算出する構成との双方を備えているが、それら構成を別々に設けるようにしてもよい。例えば、第1、2の実施形態の光画像計測装置1、100の変形例として、信号処理部20を有さない構成を採用することができる。そのとき、信号処理部20は、当該変形例の光画像計測装置とは別個に設けられたコンピュータ装置等により構成され、光画像計測装置は、この外部のコンピュータ装置等に対して補助検出信号を出力する。当該変形例を採用する場合、本発明にいう「取得手段」は、第1の実施形態に係る当該変形例においては、広帯域光源41、ハーフミラー42、固定鏡43、ガルバノミラー44、波長フィルタ45、回折格子46及び1次元光センサアレイ47により構成され、第2の実施形態に係る当該変形例においては、ビームスプリッタ71、回折格子72及び1次元光センサアレイ73により構成される。これは、当該変形例の取得手段により得られる補助検出信号は、断層画像の計測深度を反映した情報を含み、また、計測深度はこの補助検出信号に基づいて得られるものであるから、補助検出信号自体を計測深度とみなせるからである。また、当該変形例の構成を採用する場合には、本発明にいう「画像形成手段」についても具備している必要はなく、例えば、検出手段からの検出信号を外部のコンピュータ装置等に出力するとともに、このコンピュータ装置等によって複数の断層画像の配列処理を実行するように構成することができる。
また、図1や図7に示した構成において、信号光Sの光路上、つまりハーフミラー6と被測定物体Oとの間に波長板(1/2波長板)を設けることにより、被測定物体Oを経由するときの位相の変化に起因する信号光Sの偏光方向の傾きを補正するように構成すれば、計測精度の更なる向上を図ることができる。
本発明に係る光画像計測装置の検出手段は、前述のCCD21、22に限定されるものではなく、例えば積算回路を備えたラインセンサなど、干渉光を検出して光電変換する機能と、検出した電荷を蓄積する機能との双方を備えているものであればよく、また1次元のセンサでも2次元のセンサでもよい。
以上の実施形態においては、マイケルソン型の干渉計を備えた光画像計測装置について説明したが、例えばマッハツェンダー型などその他の干渉計を採用することも可能である(例えば、本発明者らによる特許第3245135号を参照)。
また、本発明者らによる上述の特許文献1(特開2001−330558号公報)に記載された、シャッタを用いて干渉光をサンプリングする構成の光画像計測装置に、本発明の構成を適用することも当然に可能である。その場合、例えば、図1に示したような広帯域光源41〜1次元光センサアレイ47からなる光学系を追加的に設けるとともに、図2に示したような制御系を適用すればよい。更に、本発明に特徴的な構成は、特許請求の範囲に記載した本発明の構成要素を満たすあらゆるタイプの光画像計測装置に適用することが可能である。
また、干渉計の一部に光ファイバ(バンドル)を設けて導光部材として用いることにより、装置設計上の自由度を高めたり、装置のコンパクト化を図ったり、あるいは、被測定物体の配置の自由度を高めたりすることができる(例えば、上記の特許第3245135号を参照)。
本発明の光画像計測装置を例えば眼科の分野に応用すると、上述した眼底の血流状態の測定のほか、網膜や角膜の2次元断層像などを得ることができる。それにより、例えば角膜の内皮細胞数などを測定することが可能となる。なお、その他にも各種の応用が可能であることはいうまでもない。
本発明に係る光画像計測方法は、例えば以上に詳述したような本発明の光画像計測装置によって実施することが可能である。
本発明に係る第1の実施形態の光画像計測装置の構成の一例を示す概略図である。 本発明に係る第1の実施形態の光画像計測装置の制御系の構成の一例を示すブロック図である。 本発明に係る第1の実施形態の光画像計測装置による補助信号光の走査態様の一例を示す概略図である。 本発明に係る第1の実施形態の光画像計測装置によるx−y断層画像の計測態様の一例を示す概略図である。 本発明に係る第1の実施形態の光画像計測装置により求められる計測深度プロファイルの一例を示す概略図である。 本発明に係る光画像計測装置における干渉光の検出態様を説明するためのグラフ図である。図6(A)は、周波数が強度変調されて広帯域光源から出力される光ビームの時間波形を表す。図6(B)は、広帯域光源から出力される光ビームが連続光であるときの干渉光のS偏光成分の時間波形を表す。図6(C)は、広帯域光源から出力される光ビームが連続光であるときの干渉光のP偏光成分の時間波形を表す。図6(D)は、広帯域光源から出力される光ビームが強度変調される場合における干渉光のS偏光成分の時間波形を表す図6(E)は、広帯域光源から出力される光ビームが強度変調される場合における干渉光のP偏光成分の時間波形を表す。 本発明に係る第2の実施形態の光画像計測装置の構成の一例を示す概略図である。 従来の光画像計測装置の構成を示す概略図である。 従来の光画像計測装置の構成を示す概略図である。 従来の光計測装置の構成を示す概略図である。 従来の光計測装置により取得される情報の態様を表す概略図である。図11(A)は、干渉光の波長スペクトルの一例の概略を表すグラフである。また、図11(B)は、被測定物体のz座標(計測深度)を変数とした干渉信号の強度分布の一例の概略を表すグラフである。
符号の説明
1、100 光画像計測装置
2 広帯域光源
3 偏光板
4、5 レンズ
6 ハーフミラー
7 波長板
8 周波数シフタ
9 参照鏡
9A ピエゾ素子
10 結像用レンズ群
11 偏光ビームスプリッタ
20 信号処理部
21、22 CCD
31 光源
32 ビームスプリッタ
33 反射鏡
34 フォトディテクタ
35 パルス駆動器
41 広帯域光源
42 ハーフミラー
43 固定鏡
44 ガルバノミラー
44A 反射鏡
44B 駆動部
45 波長フィルタ
46 回折格子
47 1次元光センサアレイ
50 制御部
51 画像形成部
52 走査制御部
53 計測深度算出部
54 スペクトル取得部
55 プロファイル取得部
56 情報記憶部
57 画像処理部
60 表示装置
71 ビームスプリッタ
72 回折格子
73 1次元光センサアレイ
R 参照光
S 信号光
L 干渉光
L1 S偏光成分
L2 P偏光成分
O 被測定物体
G1〜GN x−y断層画像
P1〜PN 計測深度算出結果
T1〜TM 補助信号光の走査位置

Claims (14)

  1. 低コヒーレント光源から出力された光ビームを被測定物体に向かう信号光と参照物体に向かう参照光とに分割する分割手段と、前記信号光の周波数と前記参照光の周波数とを相対的にシフトさせる周波数シフト手段と、前記参照光の光路長を変更する変更手段と、前記被測定物体を経由した前記信号光と前記参照物体を経由した前記参照光とを重畳させて前記周波数のシフト量に応じたビート周波数の干渉光を生成する重畳手段と、前記生成された干渉光を受光して検出信号を出力する検出手段と、前記出力された検出信号に基づいて前記信号光の進行方向に直交する前記被測定物体の断層画像を形成する画像形成手段とを有する光画像計測装置であって、
    前記形成される前記被測定物体の前記断層画像について、前記信号光の進行方向における計測深度を求める取得手段と、
    前記断層画像と、前記断層画像について取得された前記計測深度の情報とを関連付けて記憶する情報記憶手段と、を備え、
    前記取得手段は、
    低コヒーレントな補助光ビームを出力する補助光源と、
    前記出力された補助光ビームを、前記被測定物体を経由する補助信号光と、補助参照物体を経由する補助参照光とに分割する補助分割手段と、
    前記被測定物体を経由した前記補助信号光と、前記補助参照物体を経由した前記補助参照光とを重畳して補助干渉光を生成する補助重畳手段と、
    前記生成された前記補助干渉光から複数の波長成分を分離する波長成分分離手段と、
    前記分離された前記複数の波長成分を受光して補助検出信号を出力する補助検出手段と、
    前記検出された前記補助検出信号に基づいて前記断層画像の前記計測深度を算出する計測深度算出手段と、
    を備えることを特徴とする光画像計測装置。
  2. 前記情報記憶手段に記憶された前記断層画像と関連付けられた計測深度の情報に基づいて、前記断層画像を前記計測深度方向に配列させる画像処理手段を更に備えることを特徴とする請求項1に記載の光画像計測装置。
  3. 前記取得手段は、前記補助分割手段により分割された前記補助信号光を前記信号光に合成して前記被測定物体に入射させ、前記被測定物体を経由した前記補助信号光を前記信号光から分離して前記補助重畳手段に導く合成分離手段を更に含むことを特徴とする請求項1又は2に記載の光画像計測装置。
  4. 前記取得手段は、
    前記合成分離手段により前記信号光に合成された前記補助信号光を、前記信号光の進行方向に対して直交する方向に走査する走査手段を更に含み、
    前記計測深度算出手段は、前記走査手段による前記補助信号光の走査に基づいて前記補助検出手段により出力される複数の前記補助検出信号に基づいて前記計測深度を算出する、
    ことを特徴とする請求項3に記載の光画像計測装置。
  5. 前記走査手段は、前記変更手段により前記参照光の光路長が変更される度毎に、前記補助信号光を複数の位置に走査させ、
    前記計測深度算出手段は、前記複数の位置のそれぞれに走査された前記補助信号光に基づく前記補助検出信号に基づいて、前記複数の位置のそれぞれに対応する複数の前記計測深度を算出し、前記算出された前記複数の計測深度に基づいて前記断層画像の前記計測深度を算出する、
    ことを特徴とする請求項4に記載の光画像計測装置。
  6. 前記走査手段は、前記補助信号光を反射する反射鏡と、前記反射鏡の反射面の向きを変更する駆動部とを含むガルバノミラーであることを特徴とする請求項4又は請求項5に記載の光画像計測装置。
  7. 前記合成分離手段は、前記信号光の光路に斜設された波長フィルタあるいはビームスプリッタであることを特徴とする請求項3ないし請求項6のいずれか一項に記載の光画像計測装置。
  8. 前記波長成分分離手段は、回折格子であることを特徴とする請求項1ないし請求項7のいずれか一項に記載の光画像計測装置。
  9. 前記補助検出手段は、1次元光センサアレイであることを特徴とする請求項1ないし請求項8のいずれか一項に記載の光画像計測装置。
  10. 前記計測深度算出手段は、
    前記補助検出信号に基づいて前記補助干渉光の波長スペクトルを求める波長スペクトル取得手段と、
    前記求められた前記波長スペクトルをフーリエ変換して、前記計測深度に応じた前記補助干渉光の強度分布を表す計測深度情報を求める計測深度情報取得手段と、
    を含み、
    前記画像形成手段により形成される前記断層画像のそれぞれについて前記求められた前記計測深度情報に基づいて前記断層画像のそれぞれの前記計測深度を算出する、
    ことを特徴とする請求項1ないし請求項9のいずれか一項に記載の光画像計測装置。
  11. 前記計測深度算出手段は、前記画像形成手段により連続して形成される一対の前記断層画像のそれぞれについて前記計測深度情報取得手段により求められた前記計測深度情報を比較して、前記一対の前記断層画像の間の計測間隔を算出することにより、前記計測深度を求めることを特徴とする請求項10に記載の光画像計測装置。
  12. 低コヒーレント光源から出力された光ビームを被測定物体に向かう信号光と参照物体に向かう参照光とに分割し、前記信号光の周波数と前記参照光の周波数とを相対的にシフトさせ、前記被測定物体を経由した前記信号光と前記参照物体を経由した前記参照光とを重畳させて前記周波数のシフト量に応じたビート周波数の干渉光を生成し、前記生成された干渉光を受光して検出信号を出力し、前記出力された検出信号に基づいて前記信号光の進行方向に直交する前記被測定物体の断層画像を形成する光画像計測方法であって、
    前記形成される前記被測定物体の前記断層画像について、前記信号光の進行方向における計測深度を求めるステップと、
    前記断層画像と、前記断層画像について取得された前記計測深度の情報とを関連付けて記憶するステップと、
    前記断層画像と関連付けられた計測深度の情報に基づいて、前記断層画像を前記計測深度方向に配列させるステップと、を有し、
    前記計測深度を求めるステップは、
    低コヒーレントな補助光ビームを出力するステップと、
    前記出力された補助光ビームを、前記信号光とともに前記被測定物体を経由する補助信号光と、補助参照物体を経由する補助参照光とに分割するステップと、
    前記被測定物体を経由した前記補助信号光と、前記補助参照物体を経由した前記補助参照光とを重畳して補助干渉光を生成するステップと、
    前記生成された前記補助干渉光から複数の波長成分を分離するステップと、
    前記分離された前記複数の波長成分を受光して補助検出信号を出力するステップと、
    前記検出された前記補助検出信号に基づいて前記断層画像の前記計測深度を算出するステップと、
    を含むことを特徴とする光画像計測方法。
  13. 前記計測深度を求めるステップは、
    前記補助信号光を、前記信号光の進行方向に対して直交する方向に走査するステップを更に含み、
    前記計測深度を算出するステップは、前記補助信号光の走査に基づいて前記出力される複数の前記補助検出信号に基づいて前記計測深度を算出する、
    ことを特徴とする請求項12に記載の光画像計測方法。
  14. 前記計測深度を算出するステップは、
    前記補助検出信号に基づいて前記補助干渉光の波長スペクトルを求めるステップと、
    前記求められた前記波長スペクトルをフーリエ変換して、前記計測深度に応じた前記補助干渉光の強度分布を表す計測深度情報を求めるステップと、
    を含み、
    前記形成される前記断層画像のそれぞれについて前記求められた前記計測深度情報に基づいて前記断層画像のそれぞれの前記計測深度を算出する、
    ことを特徴とする請求項12又は請求項13のいずれか一項に記載の光画像計測方法。
JP2004299036A 2004-10-13 2004-10-13 光画像計測装置及び光画像計測方法 Expired - Fee Related JP4566685B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2004299036A JP4566685B2 (ja) 2004-10-13 2004-10-13 光画像計測装置及び光画像計測方法
US11/246,394 US7492466B2 (en) 2004-10-13 2005-10-11 Optical image measuring apparatus and optical image measuring method
CN200510112767.4A CN100582741C (zh) 2004-10-13 2005-10-12 光图像计测装置及光图像计测方法
EP05022289A EP1650528A3 (en) 2004-10-13 2005-10-12 Apparatus and method of heterodyne interferometry for imaging

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2004299036A JP4566685B2 (ja) 2004-10-13 2004-10-13 光画像計測装置及び光画像計測方法

Publications (2)

Publication Number Publication Date
JP2006112864A JP2006112864A (ja) 2006-04-27
JP4566685B2 true JP4566685B2 (ja) 2010-10-20

Family

ID=35677235

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004299036A Expired - Fee Related JP4566685B2 (ja) 2004-10-13 2004-10-13 光画像計測装置及び光画像計測方法

Country Status (4)

Country Link
US (1) US7492466B2 (ja)
EP (1) EP1650528A3 (ja)
JP (1) JP4566685B2 (ja)
CN (1) CN100582741C (ja)

Families Citing this family (59)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1576337B1 (en) * 2002-12-18 2010-07-14 Koninklijke Philips Electronics N.V. Method and arrangement for optical coherence tomography
JP4409332B2 (ja) * 2004-03-30 2010-02-03 株式会社トプコン 光画像計測装置
JP4916779B2 (ja) * 2005-09-29 2012-04-18 株式会社トプコン 眼底観察装置
US7295293B2 (en) * 2005-10-21 2007-11-13 Hewlett-Packard Development Company, L.P. Apparatus and method for testing a reflector coating
CN101322025B (zh) * 2005-12-07 2011-08-03 株式会社拓普康 光图像测量装置
CA2649065A1 (en) * 2006-05-01 2007-11-15 Physical Sciences, Inc. Hybrid spectral domain optical coherence tomography line scanning laser ophthalmoscope
JP4969925B2 (ja) * 2006-06-28 2012-07-04 株式会社トプコン 眼底観察装置
JP5095167B2 (ja) * 2006-09-19 2012-12-12 株式会社トプコン 眼底観察装置、眼底画像表示装置及び眼底観察プログラム
JP4895277B2 (ja) * 2006-09-29 2012-03-14 富士フイルム株式会社 光断層画像化装置
JP5406427B2 (ja) * 2006-11-17 2014-02-05 株式会社トプコン 断層画像処理方法、装置およびプログラムならびにこれを用いた光断層画像化システム
JP5007114B2 (ja) * 2006-12-22 2012-08-22 株式会社トプコン 眼底観察装置、眼底画像表示装置及びプログラム
JP4996918B2 (ja) * 2006-12-26 2012-08-08 株式会社トプコン 光画像計測装置及び光画像計測装置を制御するプログラム
JP4996917B2 (ja) * 2006-12-26 2012-08-08 株式会社トプコン 光画像計測装置及び光画像計測装置を制御するプログラム
EP1962051A1 (de) * 2007-02-21 2008-08-27 Agfa HealthCare N.V. System und Verfahren zur optischen Kohärenztomographie
EP1962082A1 (de) * 2007-02-21 2008-08-27 Agfa HealthCare N.V. System und Verfahren zur optischen Kohärenztomographie
JP5058627B2 (ja) * 2007-02-26 2012-10-24 株式会社トプコン 眼底観察装置
JP5523658B2 (ja) * 2007-03-23 2014-06-18 株式会社トプコン 光画像計測装置
JP5448353B2 (ja) 2007-05-02 2014-03-19 キヤノン株式会社 光干渉断層計を用いた画像形成方法、及び光干渉断層装置
CN100458358C (zh) * 2007-07-10 2009-02-04 浙江大学 一种基于轴向立体视觉的逆向测量方法与装置
WO2009033111A2 (en) * 2007-09-06 2009-03-12 Lensx Lasers, Inc. Precise targeting of surgical photodisruption
WO2009114852A2 (en) * 2008-03-14 2009-09-17 The Trustees Of Columbia University In The City Of New York Systems and methods for dynamic imaging of tissue using digital optical tomography
JP4236123B1 (ja) * 2008-04-21 2009-03-11 株式会社林創研 三次元画像取得装置
US7884946B2 (en) * 2008-04-28 2011-02-08 Lumetrics, Inc. Apparatus for measurement of the axial length of an eye
JP5179265B2 (ja) * 2008-06-02 2013-04-10 株式会社ニデック 眼科撮影装置
JP5306075B2 (ja) * 2008-07-07 2013-10-02 キヤノン株式会社 光干渉断層法を用いる撮像装置及び撮像方法
US8471895B2 (en) * 2008-11-25 2013-06-25 Paul S. Banks Systems and methods of high resolution three-dimensional imaging
JP5259374B2 (ja) * 2008-12-19 2013-08-07 富士フイルム株式会社 光構造観察装置及びその構造情報処理方法
US9700210B2 (en) * 2009-03-02 2017-07-11 Canon Kabushiki Kaisha Image processing apparatus and method for controlling the same
JP5400481B2 (ja) 2009-06-03 2014-01-29 株式会社トプコン 光画像計測装置
US9492322B2 (en) 2009-11-16 2016-11-15 Alcon Lensx, Inc. Imaging surgical target tissue by nonlinear scanning
US8265364B2 (en) * 2010-02-05 2012-09-11 Alcon Lensx, Inc. Gradient search integrated with local imaging in laser surgical systems
US8414564B2 (en) * 2010-02-18 2013-04-09 Alcon Lensx, Inc. Optical coherence tomographic system for ophthalmic surgery
JP6039156B2 (ja) * 2010-06-08 2016-12-07 キヤノン株式会社 画像処理装置、画像処理方法、及びプログラム
US8398236B2 (en) 2010-06-14 2013-03-19 Alcon Lensx, Inc. Image-guided docking for ophthalmic surgical systems
CN101949689B (zh) * 2010-06-22 2012-05-30 深圳市斯尔顿科技有限公司 一种oct***校正方法
JP5657941B2 (ja) * 2010-07-30 2015-01-21 株式会社トプコン 光断層画像化装置及びその作動方法
US9532708B2 (en) 2010-09-17 2017-01-03 Alcon Lensx, Inc. Electronically controlled fixation light for ophthalmic imaging systems
JP5794664B2 (ja) * 2011-01-20 2015-10-14 キヤノン株式会社 断層画像生成装置及び断層画像生成方法
WO2012102887A2 (en) * 2011-01-24 2012-08-02 The Board Of Trustees Of The University Of Illinois Computational adaptive optics for interferometric synthetic aperture microscopy and other interferometric imaging
US8459794B2 (en) 2011-05-02 2013-06-11 Alcon Lensx, Inc. Image-processor-controlled misalignment-reduction for ophthalmic systems
US9622913B2 (en) 2011-05-18 2017-04-18 Alcon Lensx, Inc. Imaging-controlled laser surgical system
US8398238B1 (en) 2011-08-26 2013-03-19 Alcon Lensx, Inc. Imaging-based guidance system for ophthalmic docking using a location-orientation analysis
US9023016B2 (en) 2011-12-19 2015-05-05 Alcon Lensx, Inc. Image processor for intra-surgical optical coherence tomographic imaging of laser cataract procedures
US9066784B2 (en) 2011-12-19 2015-06-30 Alcon Lensx, Inc. Intra-surgical optical coherence tomographic imaging of cataract procedures
US9677869B2 (en) 2012-12-05 2017-06-13 Perimeter Medical Imaging, Inc. System and method for generating a wide-field OCT image of a portion of a sample
CN103630077B (zh) * 2013-12-12 2016-05-11 哈尔滨工业大学 一种使用双频激光的两轴光栅位移测量***
WO2015164868A1 (en) 2014-04-26 2015-10-29 Tetravue, Inc. Method and system for robust and extended illumination waveforms for depth sensing in 3d imaging
FR3033643B1 (fr) * 2015-03-13 2020-07-17 Unity Semiconductor Dispositif et procede pour detecter des defauts dans des zones de liaison entre des echantillons tels que des wafers
US10337983B2 (en) 2015-04-12 2019-07-02 Taiwan Biophotonic Corporation Module, device and method for optical measurement
CN112969031B (zh) 2015-06-17 2022-11-15 松下知识产权经营株式会社 摄像装置
WO2017068878A1 (ja) * 2015-10-19 2017-04-27 ソニー株式会社 測距装置及び撮像システム
WO2017086442A1 (ja) * 2015-11-18 2017-05-26 株式会社ニコン 撮像素子、計測装置および計測方法
JP2019513976A (ja) 2016-02-29 2019-05-30 テトラビュー, インコーポレイテッド 3d結像システムおよび方法
DE102016218290A1 (de) * 2016-07-15 2018-01-18 Carl Zeiss Meditec Ag Verfahren zur hochsensitiven Messung von Abständen und Winkeln im menschlichen Auge
CN108209867B (zh) 2016-12-15 2022-03-18 松下知识产权经营株式会社 摄像装置
EP3655748B1 (en) 2017-07-18 2023-08-09 Perimeter Medical Imaging, Inc. Sample container for stabilizing and aligning excised biological tissue samples for ex vivo analysis
CN109856108B (zh) * 2017-11-30 2021-05-18 北京华泰诺安探测技术有限公司 一种拉曼光谱检测***和方法
US11212512B2 (en) 2017-12-28 2021-12-28 Nlight, Inc. System and method of imaging using multiple illumination pulses
US10925767B2 (en) 2018-03-06 2021-02-23 Lutronic Vision Inc Laser doppler vibrometry for eye surface vibration measurement to determine cell damage

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001066245A (ja) * 1999-08-26 2001-03-16 Japan Science & Technology Corp 光波反射断層像観測装置
JP2001330558A (ja) * 2000-05-25 2001-11-30 Japan Science & Technology Corp 2次元光ヘテロダイン検出法を用いた光画像計測装置
JP2003035660A (ja) * 2001-07-19 2003-02-07 Japan Science & Technology Corp 分光機能を備えた光干渉断層画像計測装置
JP2003050199A (ja) * 2001-08-03 2003-02-21 Olympus Optical Co Ltd 光イメージング装置

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02173505A (ja) * 1988-12-26 1990-07-05 Brother Ind Ltd 光波干渉型微細表面形状測定装置
JPH0339605A (ja) * 1989-07-05 1991-02-20 Brother Ind Ltd 光学式表面形状測定装置
US5051226A (en) * 1989-09-18 1991-09-24 The Boeing Company Method of curing composite parts
US5838485A (en) * 1996-08-20 1998-11-17 Zygo Corporation Superheterodyne interferometer and method for compensating the refractive index of air using electronic frequency multiplication
US6002480A (en) * 1997-06-02 1999-12-14 Izatt; Joseph A. Depth-resolved spectroscopic optical coherence tomography
JP4138027B2 (ja) * 1997-06-02 2008-08-20 イザット,ジョーゼフ,エイ. 光学コヒーレンス断層撮影法を用いたドップラー流の撮像
JP2000292705A (ja) * 1999-04-05 2000-10-20 Olympus Optical Co Ltd 走査型顕微鏡
EP1149556B1 (en) * 2000-04-26 2007-08-15 FUJIFILM Corporation Optical tomograph
DE10123844A1 (de) * 2001-04-09 2002-10-17 Bosch Gmbh Robert Interferometrische Messvorrichtung
WO2003011764A2 (en) * 2001-08-03 2003-02-13 Volker Westphal Real-time imaging system and method
US7006232B2 (en) * 2002-04-05 2006-02-28 Case Western Reserve University Phase-referenced doppler optical coherence tomography
AU2003247550A1 (en) * 2002-06-17 2003-12-31 Zygo Corporation Interferometry methods and systems having a coupled cavity geometry for use with an extended source
JP4381847B2 (ja) * 2004-02-26 2009-12-09 株式会社トプコン 光画像計測装置
JP4522724B2 (ja) * 2004-03-16 2010-08-11 株式会社トプコン 光画像計測装置
US7126693B2 (en) * 2004-03-29 2006-10-24 Carl Zeiss Meditec, Inc. Simple high efficiency optical coherence domain reflectometer design
JP4409332B2 (ja) * 2004-03-30 2010-02-03 株式会社トプコン 光画像計測装置
JP4409384B2 (ja) * 2004-08-03 2010-02-03 株式会社トプコン 光画像計測装置及び光画像計測方法
JP4633423B2 (ja) * 2004-09-15 2011-02-16 株式会社トプコン 光画像計測装置
JP4563130B2 (ja) * 2004-10-04 2010-10-13 株式会社トプコン 光画像計測装置
JP4494160B2 (ja) * 2004-10-14 2010-06-30 株式会社トプコン 光画像計測装置
JP4597744B2 (ja) * 2004-11-08 2010-12-15 株式会社トプコン 光画像計測装置及び光画像計測方法
US7355719B2 (en) * 2005-08-16 2008-04-08 Agilent Technologies, Inc. Interferometer for measuring perpendicular translations

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001066245A (ja) * 1999-08-26 2001-03-16 Japan Science & Technology Corp 光波反射断層像観測装置
JP2001330558A (ja) * 2000-05-25 2001-11-30 Japan Science & Technology Corp 2次元光ヘテロダイン検出法を用いた光画像計測装置
JP2003035660A (ja) * 2001-07-19 2003-02-07 Japan Science & Technology Corp 分光機能を備えた光干渉断層画像計測装置
JP2003050199A (ja) * 2001-08-03 2003-02-21 Olympus Optical Co Ltd 光イメージング装置

Also Published As

Publication number Publication date
CN100582741C (zh) 2010-01-20
US7492466B2 (en) 2009-02-17
CN1760663A (zh) 2006-04-19
EP1650528A3 (en) 2012-01-25
JP2006112864A (ja) 2006-04-27
EP1650528A2 (en) 2006-04-26
US20060077395A1 (en) 2006-04-13

Similar Documents

Publication Publication Date Title
JP4566685B2 (ja) 光画像計測装置及び光画像計測方法
JP4409384B2 (ja) 光画像計測装置及び光画像計測方法
JP4633423B2 (ja) 光画像計測装置
JP4494160B2 (ja) 光画像計測装置
JP4837300B2 (ja) 光画像計測装置
JP4409331B2 (ja) 光画像計測装置
US7486402B2 (en) Optical image measuring apparatus
EP1568963B1 (en) Interferometric apparatus for measuring shapes
EP1577637B1 (en) Optical image measuring apparatus
JP4602372B2 (ja) 光画像計測装置及び光画像計測方法
JP5400481B2 (ja) 光画像計測装置
US20050219545A1 (en) Optical image measuring apparatus
JP5903903B2 (ja) 光コヒーレンストモグラフィー装置
JP4409334B2 (ja) 光画像計測装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20071012

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20081217

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20091224

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20100105

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100304

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20100525

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100715

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: 20100803

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: 20100804

R150 Certificate of patent or registration of utility model

Ref document number: 4566685

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130813

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees