JP4409331B2 - 光画像計測装置 - Google Patents

光画像計測装置 Download PDF

Info

Publication number
JP4409331B2
JP4409331B2 JP2004100741A JP2004100741A JP4409331B2 JP 4409331 B2 JP4409331 B2 JP 4409331B2 JP 2004100741 A JP2004100741 A JP 2004100741A JP 2004100741 A JP2004100741 A JP 2004100741A JP 4409331 B2 JP4409331 B2 JP 4409331B2
Authority
JP
Japan
Prior art keywords
light
interference light
interference
frequency
measured
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 - Lifetime
Application number
JP2004100741A
Other languages
English (en)
Other versions
JP2005283471A (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 JP2004100741A priority Critical patent/JP4409331B2/ja
Priority to US11/091,478 priority patent/US7548320B2/en
Priority to EP05006897A priority patent/EP1582143B1/en
Publication of JP2005283471A publication Critical patent/JP2005283471A/ja
Application granted granted Critical
Publication of JP4409331B2 publication Critical patent/JP4409331B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/12Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes
    • A61B3/1225Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes using coherent radiation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B3/00Apparatus for testing the eyes; Instruments for examining the eyes
    • A61B3/10Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
    • A61B3/102Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for optical coherence tomography [OCT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0062Arrangements for scanning
    • A61B5/0066Optical coherence imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0073Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by tomography, i.e. reconstruction of 3D images from 2D projections
    • 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/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
    • 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
    • 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
    • G01B9/02028Two or more reference or object arms in one interferometer
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Engineering & Computer Science (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Ophthalmology & Optometry (AREA)
  • Instruments For Measurement Of Length By Optical Means (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Description

本発明は、特に光散乱媒質の被測定物体に光ビームを照射し、その反射光もしくは透過光を用いて被測定物体の表面形態や内部形態を計測し、その画像を形成する光画像計測装置に関する。特に、光ヘテロダイン検出法を利用して被測定物体の表面形態や内部形態を計測し画像を形成する光画像計測装置に関するものである。
近年、レーザ光源等を用いて被測定物体の表面や内部の画像を形成する光画像計測技術が注目を集めている。この光画像計測技術は、従来からのX線CTのような人体に対する有害性を持たないことから、特に医療分野における応用の展開が期待されている。
光画像計測技術における代表的な手法の一例として低コヒーレンス干渉法(光コヒーレンス断層画像化法などとも呼ばれる)がある。この手法は、例えばスーパールミネセントダイオード(Super Luminescent Diode;SLD)のように広いスペクトル幅をもつ広帯域光源の低干渉性を利用して、被測定物体からの反射光や透過光をμmオーダーの優れた距離分解能で検出可能とするものである(例えば下記の非特許文献1を参照)。
この低コヒーレンス干渉法を利用した装置の一例として、マイケルソン干渉計に基づく従来の光画像計測装置の基本構成を図6に示す。この光画像計測装置100は、広帯域光源101、鏡102、ビームスプリッタ103及び光検出器104を含んで構成されている。被測定物体105は、散乱媒質により形成されている。広帯域光源101からの光ビームは、ビームスプリッタ103により、鏡102に向かう参照光Rと被測定物体105に向かう信号光Sの2つに分割される。参照光Rはビームスプリッタ103による反射光であり、信号光Sはビームスプリッタ103の透過光である。
ここで、図6に示すように、信号光Sの進行方向をz軸方向として定めるとともに、信号光Sの進行方向に対する直交面をx−y面として定義する。鏡102は、同図中の両側矢印方向(z−スキャン方向)に変位可能とされている。
参照光Rは、鏡102に反射される際にそのz−スキャンによりドップラー周波数シフトを受ける。一方、信号光Sは、被測定物体105に照射されると、その表面及び内部層により反射される。被測定物体は散乱媒質であるので、信号光Sの反射光は多重散乱を含む乱雑な位相をもった拡散波面であると考えられる。被測定物体105を経由した信号光と、鏡102を経由し周波数シフトを受けた参照光は、ビームスプリッタ103によって重畳されて干渉光を生成する。
このような低コヒーレンス干渉方法を用いた画像計測では、信号光Sと参照光Rの光路長差が光源のμmオーダーのコヒーレント長(可干渉距離)以内でありかつ参照光Rと位相相関のある信号光Sの成分のみが参照光Rと干渉を生じる。すなわち、信号光Sのコヒーレントな信号光成分のみが選択的に参照光Rと干渉し合う。この原理から、鏡102の位置をz−スキャンして参照光Rの光路長を変化させることにより、被測定物体105の内部層の光反射プロフィールが測定される。更に、被測定物体105へ照射される信号光Sをx−y面方向にも走査する。このようなz方向及びx−y面方向のスキャンを行いながら干渉光を光検出器104で検出し、その検出結果として出力される電気信号(ヘテロダイン信号)を解析することによって、被検体105の2次元断層画像が取得される(非特許文献1を参照)。
なお、ビームスプリッタ103によって重畳される参照光R及び信号光Sの強度をそれぞれIとIとし、両光波の間の周波数差及び位相差をそれぞれfif及びΔθとすると、光検出器からは次式に示すようなヘテロダイン信号が出力される(例えば、非特許文献2を参照)。
Figure 0004409331
式(1)の右辺第3項は交流電気信号であり、その周波数fifは参照光Rと信号光Sが干渉されて発生するうなり(ビート)の周波数に等しい。ヘテロダイン信号の交流成分の周波数fifは、ビート周波数などと呼ばれる。また、式(1)の右辺第1項及び第2項はヘテロダイン信号の直流成分であり、干渉光の背景光の信号強度に対応している。
しかしながら、このような従来の低コヒーレンス干渉法で2次元断層画像を取得するためには、被測定物体105に対して光ビームを走査することにより、被測定物体105の深度方向(z方向)及び断層面方向(x−y面方向)の各部位からの反射光波を順次に検出する必要がある。したがって、被測定物体105を計測するためには長時間を要し、また、その計測原理を勘案すると計測時間の短縮を図ることは困難である。
このような問題点に鑑み、計測時間の短縮を図るための光画像計測装置が考案されている。図7は、そのような装置の一例の基本構成が示されている。同図に示す光画像計測装置200は、広帯域光源201、鏡202、ビームスプリッタ203、光検出器としての2次元光センサアレイ204、及びレンズ206、207を含んで構成されている。光源201から出射された光ビームは、レンズ206、207により平行光束にされるとともにそのビーム径が拡大され、ビームスプリッタ203によって参照光Rと信号光Sとに二分される。参照光Rには鏡202のz−スキャンによりドップラー周波数シフトが付与される。一方、信号光Sは、ビーム径が広げられているので、x−y面の広い範囲に亘って被測定物体205に入射される。よって、信号光Sは、当該入射範囲における被測定物体205の表面や内部の情報を含んだ反射光となる。参照光Rと信号光Sは、ビームスプリッタ203により重畳され、2次元光センサアレイ204上に並列配置された素子(光センサ)により検出される。したがって、光ビームを走査することなく、被測定物体205の2次元断層画像をリアルタイムに取得することが可能となる。
このような非走査型の光画像計測装置としては、非特許文献3に記載のものが知られている。同文献に記載の装置では、2次元光センサアレイから出力される複数のヘテロダイン信号を並列配置の信号処理系に入力して、各ヘテロダイン信号の振幅と位相を検出するようになっている。
しかしながら、画像の空間分解能を高めるためにはアレイの素子数を増加させなければならず、更に、その素子数に対応するチャンネル数を備えた信号処理系を用意しなければならない。したがって、高分解能の画像を必要とする医療や工業等の分野においては実用化する上で難があると思われる。
そこで、本発明者らは、下記の特許文献1において、次のような非走査型の光画像計測装置を提案した。本提案に係る光画像計測装置は、光ビームを出射する光源と、該光源から出射された光ビームを、被検体が配置される被検体配置位置を経由する信号光と、前記被検体配置位置を経由する光路とは異なる光路を経由する参照光とに二分するとともに、前記被検体配置位置を経由した後の信号光と、前記異なる光路を経由した参照光とを互いに重畳することにより干渉光を生成する干渉光学系と、該干渉光学系が、前記信号光の周波数と前記参照光の周波数を相対的にシフトさせる周波数シフタと、前記干渉光学系が、前記干渉光を受光するために、前記干渉光を二分割して、さらに、該二分割された干渉光を周期的に遮断することにより、互いの位相差が90度である2列の干渉光パルスを生成する光遮断装置と、前記2列の干渉光パルスをそれぞれ受光する光センサと、該光センサが、空間的に配列され、それぞれが独立に受光信号を得る複数の受光素子を有するものであり、前記光センサで得られた複数の受光信号を統合して前記被検体配置位置に配置された被検体の表面もしくは内部層の、前記信号光の伝搬経路上の各関心点に対応する信号を生成する信号処理部を具備している。
この光画像計測装置は、参照光と信号光からなる干渉光を二分して2台の光センサ(2次元光センサアレイ)で受光するとともに、両センサアレイの前にそれぞれ光遮断装置(シャッタ)を配置して干渉光をサンプリングするように構成されている。そして、分割された2つの干渉光のサンプリング周期にπ/2の位相差を設けることにより、干渉光の背景光を構成する信号光と参照光の強度と、干渉光の位相の直交成分(sin成分とcos成分)とを検出し、更に、両センサアレイからの出力に含まれる背景光の強度を両センサアレイの出力から差し引くことにより、干渉光の2つの位相直交成分を算出し、その算出結果を用いて干渉光の振幅を求めるようになっている。
また、本発明者らは、特許文献2において、光ビームを出射する光源と、該光源から出射された光ビームを、被測定検体が配置される被測定検体の配置位置を経由する信号光と、前記被測定検体の配置位置を経由する光路とは異なる光路を経由する参考光とに二分するとともに、前記被測定検体の配置位置を経由した後の信号光と、前記異なる光路を経由した参考光とを互いに重畳することにより前記信号光と前記参照光とが干渉した干渉光を生成する干渉光学系と、該干渉光学系が、前記信号光の周波数と前記参照光の周波数を相対的にシフトさせる周波数シフタと、前記干渉光学系が、前記信号光と前記参照光のうち少なくとも一方の光路上に光を周期的に遮断する光学装置とを備え、前記光学装置の遮断周波数を前記信号光と前記参照光の間の周波数差に等しくなるように構成された光画像計測装置を提案した。当該光画像計測装置によれば、ビート周波数と等しい遮断周波数にて干渉光をサンプリングすることができるため、好適な光ヘテロダイン測定の実現が図られている。
特許文献1、2に記載の光画像計測装置は、被測定物体の単一の深さ領域からの反射光(あるいは透過光)に基づく干渉光を基に、当該深さ領域の画像を取得するものである。したがって、被測定物体の複数の深さ領域の画像を取得する場合には、各深さ領域についてそれぞれスキャンしながら計測しなければならないため、計測作業には長時間を要していた。
特開2001−330558号公報(請求項、明細書段落、[0044]、[0072]−[0077]) 特許第3245135号(請求項、明細書段落[0072]−[0082]) 丹野直弘、「光学」、28巻3号、116(1999) 吉沢、瀬田編、「光ヘテロダイン技術(改訂版)」、新技術コミュニケーションズ(2003)、p.2 K.P.Chan、K.Satori、H.Inaba、「Electronics Letters」、Vol.34、1101(1998)
本発明は、以上のような事情に鑑みてなされたもので、被測定物体の複数の深さ領域に対応する成分を含む干渉光を生成するとともに、当該干渉光から各成分を選択的に検出することにより、各深さ領域の画像を並列的に取得することを可能とする光画像計測装置を提供することを目的としている。
上記目的を達成するために、請求項1に記載の発明は、光ビームを出射する光源と、この光源からの光ビームを平行光束にするとともにそのビーム径を拡大するレンズ系と、このレンズ系により平行光束とされてビーム径が拡大された光ビームを基に、被測定物体の深さ領域毎に異なる周波数を有する複数の干渉光成分を含んだ前記拡大されたビーム径を有する干渉光を生成する干渉系と、前記生成された干渉光から前記各干渉光成分を選択的に検出する検出系と、前記検出された前記各干渉光成分に基づいて、当該各干渉光成分に対応する前記深さ領域における前記被測定物体の2次元断面画像を形成する画像形成手段と、を備えていることを特徴とする光画像計測装置である。
また、請求項2に記載の発明は、請求項1に記載の光画像計測装置であって、前記干渉系は、前記被測定物体からの距離が各々異なるように配置された複数の参照物体と、前記レンズ系により平行光束とされてビーム径が拡大された光ビームを前記被測定物体の複数の前記深さ領域にて反射される前記拡大されたビーム径を有する信号光と、前記複数の参照物体を各々経由する前記拡大されたビーム径を有する複数の参照光とに分割する分割手段と、前記複数の参照光の周波数を各々異なる量だけシフトする周波数シフト手段と、前記被測定物体を経由した前記信号光と、前記周波数シフト手段により周波数が各々シフトされ、前記複数の参照物体を各々経由した前記複数の参照光とを重畳して前記干渉光を生成する重畳手段と、を含んでいることを特徴とする。
また、請求項3に記載の発明は、請求項1又は請求項2に記載の光画像計測装置であって、前記干渉系は、マイケルソン干渉計を用いた構成とされていることを特徴とする。
また、請求項4に記載の発明は、請求項2に記載の光画像計測装置であって、前記検出系は、前記生成された前記干渉光の強度を前記周波数のシフト量に同期した変調周波数で周期的に変調することにより、当該変調周波数に対応する前記干渉光成分を前記干渉光から抽出する強度変調手段と、前記抽出された前記干渉光成分を受光し、電気信号に変換して出力する2次元光センサアレイと、を含んでいることを特徴とする。
また、請求項5に記載の発明は、請求項4に記載の光画像計測装置であって、前記強度変調手段は、前記干渉光を前記変調周波数にて周期的に遮断するシャッタであることを特徴とする。
また、請求項6に記載の発明は、請求項4に記載の光画像計測装置であって、前記検出系は、前記生成された干渉光を複数に分割する干渉光分割手段を更に備え、前記強度変調手段及び前記2次元光センサアレイは、前記分割された前記複数の干渉光の各光路上に各々配置されていることを特徴とする。
また、請求項7に記載の発明は、光ビームを出射する光源と、この光源からの光ビームを平行光束にするとともにそのビーム径を拡大するレンズ系と、このレンズ系により平行光束とされてビーム径が拡大された光ビームを、被測定物体の複数の深さ領域にて反射される前記拡大されたビーム径を有する信号光と所定の参照物体により反射される前記拡大されたビーム径を有する参照光とに分割する第1のビームスプリッタと、前記参照光を複数に分割する第2のビームスプリッタと、前記第2のビームスプリッタにより分割された前記複数の参照光の各光路上に配置され、前記複数の参照光の周波数を各々異なる量だけシフトさせる周波数シフタと、前記複数の参照光の各光路上に配置され、前記被測定物体からの距離が各々異なるように配置された前記参照物体としての反射板と、前記周波数シフタにより前記周波数が各々シフトされ、前記反射板により各々反射された前記複数の前記参照光と、前記被測定物体の前記複数の深さ領域にて反射された前記信号光とを重畳することにより、前記各深さ領域毎に異なる周波数を有する複数の干渉光成分を含んだ前記拡大されたビーム径を有する干渉光を生成する第3のビームスプリッタと、前記生成された前記干渉光を前記周波数のシフト量のいずれかに同期した周波数で周期的に遮断することにより、当該同期した周波数に対応する前記干渉光成分を前記干渉光から抽出するシャッタと、前記抽出された前記干渉光成分を受光し、電気信号に変換して出力する2次元光センサアレイと、前記出力された電気信号に基づいて前記干渉光成分の強度及び位相の分布を演算し、当該強度及び位相の分布に基づいて前記干渉光成分に対応する前記深さ領域における前記被測定物体の反射率分布を演算し、当該反射率分布に基づいて当該深さ領域における2次元断面画像を形成するコンピュータと、を備えていることを特徴とする光画像計測装置である。
なお、本発明において、「周波数のシフト量に同期した(変調)周波数」の「同期」とは、(変調)周波数が当該シフト量と(ほぼ)同一である場合、(変調)周波数が当該シフト量のn倍である場合、及び1/n倍である場合を示すものとする(ここで、nは2以上の整数とする)。
また、「深さ領域毎に異なる周波数」の「異なる」とは、各深さ領域に対応する周波数が、他の深さ領域に対応する周波数と同一でも、n倍でも、1/n倍でもないことを示していることとする。同様に、「周波数を各々異なる量だけシフトする」の「異なる」についても、各周波数のシフト量が、他の周波数のシフト量と同一でも、n倍でも、1/n倍でもないことを示すものとする。
本発明の光画像計測装置によれば、被測定物体の各深さ領域に対応する干渉光成分を含んだ干渉光を生成するとともに、当該干渉光から各干渉光成分を選択的に検出し、各干渉光成分に対応する深さ領域の2次元断面画像を形成することにより、各深さ領域の画像を並列的に取得することができる。
以下、本発明に係る光画像計測装置の実施形態の一例について、図面を参照しながら詳細に説明する。
〔装置の構成〕
図1は、本発明に係る光画像計測装置1の好適な一実施形態の概略構成を示している。光画像計測装置1は、例えば医療用や工業用に利用可能な装置であって、散乱媒質からなる被測定物体Oの複数の深さ領域の2次元断層画像を並列的に取得するための構成を備えている。ここで、被測定物体Oの深さ領域とは、被測定物体Oの或る深さ位置(図1中に示す3次元座標の或るz座標)における領域を示している。具体的には、図1中の符号D1〜Dnで示すように、被測定物体Oの深さ方向(z方向)の或る位置においてxy平面方向に広がる範囲と、その深さ位置前後の所定のz方向の範囲とによって決定される領域である。当該所定のz方向の範囲は、干渉光Lにより検出可能なz方向の幅を示している。
光画像計測装置1は、低コヒーレントな光を出力するSLDや発光ダイオード(LED)等からなる広帯域光源2と、この光源2からの光ビームを平行光束とするとともにそのビーム径を拡大するレンズ系3と、光ビームを信号光Sと参照光Rとに分割するビームスプリッタ4とを含んでいる。光源2は、連続光やパルス光などを適宜発光可能とされている。ちなみに、市販の近赤外域SLDのコヒーレント長は30μm程度、LEDの場合には10μm程度である。信号光Sは、被測定物体Oを経由する光路を伝搬し、被測定物体Oの複数の深さ領域D1〜Dnにおいて反射される光ビームである。参照光Rは、所定の参照物体(後述の反射板7−1〜7−n)を経由する光路を伝搬する光ビームである。
光画像計測装置1には、第1の参照系A1〜第nの参照系Anとして示すn個の参照系が設けられている。第1の参照系A1は、電器光学変調器や音響光学変調器等からなる周波数シフタ6−1と、反射板7−1とを備えている。第2の参照系A2〜第nの参照系Anも同様に構成されている。ここで、nは2以上の任意の整数であるものとする。
ビームスプリッタ4によって分割された参照光Rの光路上には、参照光Rをn個の参照光R1〜Rnに更に分割し、各参照光R1〜Rnを対応する第1の参照系A1〜第nの参照系Anにそれぞれ導く、ビームスプリッタ5−1〜5−(n−1)及び反射鏡5−nが配設されている。なお、参照系Anを参照光Rの光路の延長上に設ける配置構成を採用する場合には、ビームスプリッタ5−(n−1)を透過した参照光を反射する反射鏡5−nは不要である。
図示は省略するが、反射板7−1〜7−nを参照光R1〜Rnの進行方向(図1中の点線矢印方向)に各々独立に移動させる駆動装置が設けられている。各反射板7−1〜7−nは、この駆動装置によって移動されることにより、被測定物体Oに対する距離を変更可能とされている。被測定物体Oの計測を行うとき、各反射板7−1〜7−nは、被測定物体Oに対して各々異なる距離となる各計測位置に配置される。各反射板7−1〜7−nの当該計測位置は、対応する被測定物体Oの深さ領域D1〜Dnの深さ(z座標)に応じて決定される。
また、上記の駆動装置によって反射板7−1等を図1中の点線矢印方向に移動させることにより参照光R1等にドップラー周波数シフトを与える構成としてもよい。その場合には周波数シフタ6−1等は不要である。ただし、参照光Rの周波数のシフト量(つまり干渉光のビート周波数)を大きくする必要がある場合などには、周波数シフタ6−1等による周波数シフトと上記ドップラー周波数シフトとを併用可能に構成することができる。
各参照系A1〜Anを経由した参照光R1〜Rnを重ね合わせて得られる参照光Rと、被測定物体Oの各深さ領域D1〜Dnにおいてそれぞれ反射された信号光Sとは、ビームスプリッタ4により重畳され、干渉光Lが生成される。すなわち、干渉光Lは、ビームスプリッタ4を透過する参照光Rの一部と、ビームスプリッタ4により反射される信号光Sの一部とが重畳されることにより生成される。
ここで、ビームスプリッタ4、ビームスプリッタ5−1〜5−(n−1)、反射鏡5−n、周波数シフタ6−1〜6−n及び反射板7−1〜7−nは、本発明にいう「干渉系」を構成している。また、ビームスプリッタ4、ビームスプリッタ5−1〜5−(n−1)、反射鏡5−nは、本発明の「分割手段」を構成し、周波数シフタ6−1〜6−nは、「周波数シフト手段」を構成している。更に、ビームスプリッタ4は、本発明の「重畳手段」、「第1のビームスプリッタ」及び「第3のビームスプリッタ」を構成している。本実施形態の光画像計測装置1は、マイケルソン干渉計を用いて構成されており、本発明の第1及び第3のビームスプリッタとして共通のビームスプリッタ4を使用しているが、マッハツェンダー干渉計など他の干渉計を用いる場合などには、第1及び第3のビームスプリッタを別々のビームスプリッタで構成してもよい。また、ビームスプリッタ5−1〜5−(n−1)は、本発明の「第2のビームスプリッタ」を構成している。なお、図示は省略するが、光画像計測装置1には、干渉光Lを結像させる結像用レンズ群など他の光学素子が適宜配置されている。
ビームスプリッタ4により生成された干渉光Lは、ビームスプリッタ8−1、8−2により3つの光路に分割される。各光路を伝搬する干渉光Lは、液晶シャッタ等の高速シャッタなどからなるシャッタ9−1、9−2、9−3によりそれぞれ所定の周波数で周期的に遮断されサンプリングされる。サンプリングされた各干渉光Lは、蓄積型の2次元光センサアレイであるCCDカメラ10−1、10−2、10−3により受光される。各CCDカメラ10−1、10−2、10−3は、受光した干渉光Lを光電変換し電気信号として出力する。
シャッタ9−1、9−2、9−3は、それぞれCCDカメラ10−1、10−2、10−3の直前に配置されている必要はなく、ビームスプリッタ8−1、8−2による干渉光Lの分岐点からCCDカメラ10−1、10−2、10−3とを結ぶ各光路上の任意の位置に配設することが可能である。すなわち、シャッタ9−1、9−2、9−3は、各干渉光Lを周期的に遮断することにより、CCDカメラ10−1、10−2、10−3による受光量を周期的に0に変調することが可能な位置に配置されていれば十分である。
ここで、ビームスプリッタ8−1、8−2は、本発明にいう「干渉光分割手段」を構成する。また、シャッタ9−1、9−2、9−3は、本発明の「強度変調手段」を構成し、CCDカメラ10−1、10−2、10−3は、「光検出手段」を構成するものである。更に、シャッタ9−1、9−2、9−3とCCDカメラ10−1、10−2、10−3は、本発明の「検出系」構成している。
更に、光画像計測装置1は、各シャッタ9−1、9−2、9−3の開閉タイミングを制御するタイミング信号を生成するための構成として、パルス信号を発生させるパルス信号発生器と、このパルス信号の位相をシフトさせて当該タイミング信号を各シャッタ9−1、9−2、9−3に出力する位相シフタとを備えている。位相シフタは、例えば各シャッタ9−1、9−2、9−3毎に設けられており、各シャッタ9−1、9−2、9−3に対するタイミング信号の位相をそれぞれ独立にシフトさせる。それにより、各シャッタ9−1、9−2、9−3は、それぞれに独立したタイミングで動作するようになっている。
各シャッタ9−1、9−2、9−3は、位相シフタからのタイミング信号に基づいて所定の周波数(本発明において変調周波数と称することとする)で干渉光Lをそれぞれ周期的に遮断することにより、各干渉光のサンプリングを行う。それにより、各CCDカメラ10−1、10−2、10−3は、対応する干渉光Lを周期的に受光することとなり、各干渉光Lを周期的なパルス光として受光する。このとき、各シャッタ9−1、9−2、9−3はそれぞれ独立して開閉されるため、CCDカメラ10−1、10−2、10−3によりそれぞれ検出される干渉光Lのパルスは、所定の位相差を有するものとなる。CCDカメラ10−1、10−2、10−3は、各ピクセルにて受光された干渉光Lを光電変換し、その変換結果であるヘテロダイン信号をコンピュータ20に出力する。ヘテロダイン信号は、検出された干渉光の強度及び位相を反映した電気信号である。
コンピュータ20は、CCDカメラ10−1、10−2、10−3から出力されるヘテロダイン信号に基づいて後述の演算処理を実行する。更に、コンピュータ20は、当該演算処理の結果を解析して、被測定物体Oの所定の深さ領域における反射率分布を求め、その反射率分布に基づいて当該深さ領域の画像を形成する。形成された画像は、コンピュータ20の図示しないモニタ装置等の表示装置に表示される。このようなコンピュータ20は、上記のような演算を行うための演算プログラムを格納するROMやハードディスクドライブ等の記憶装置と、当該演算プログラムを実行するCPUとを含んで構成される。当該コンピュータ20は、本発明にいう「画像形成手段」を構成している。なお、本発明の画像形成手段としては、このようなコンピュータを光学的構成を格納する筐体の外部に設けたり、同一筐体内に配置させたりすることができる。
〔測定形態〕
続いて、以上のような構成の光画像計測装置1による測定形態について説明する。まず、オペレータは、被測定物体Oの各深さ領域D1〜Dnを計測可能な位置に反射板7−1〜7−nを移動させる。具体的には、反射板7−1〜7−nは、各参照光R1〜Rnと、対応する深さ領域D1〜Dnにて反射された信号光Sの成分が干渉し合う位置に配置される。すなわち、ビームスプリッタ4に対する反射板7−kと深さ領域Dkとの距離がほぼ等しくなるように、反射板7−kは配置される(k=1〜n)。また、各周波数シフタ6−1〜6−nの周波数シフト量が設定される。当該シフト量は、オペレータが手動で設定してもよいし、所定のデフォルト値を使用してもよい。また、各周波数シフタ6−1〜6−nの設定に応じて、各シャッタ9−1〜9−nの開閉動作の周波数が設定される(詳細は後述する)。
光源2から出射された光ビームは、レンズ系3によってそのビーム径が拡大され、ビームスプリッタ4により信号光Sと参照光Rとに分割される。
信号光Sは、被測定物体Oに入射され、その表面及び内部にて反射される。当該反射光は、各深さ領域D1〜Dnにより反射された成分を含んでいる。このとき、被測定物体Oは散乱媒質であるので、被測定物体Oにより反射された信号光Sは多重散乱を含む乱雑な位相をもった拡散波面となる。被測定物体Oを経由した信号光Sは、被測定物体Oのビームスプリッタ4に再び入射される。
一方、参照光Rは、ビームスプリッタ5−1によりその一部(参照光R1)が反射されて第1の参照系A1に入射する。参照光R1は、周波数シフタ6−1を通過して周波数がシフトされた後、反射板7−1に伝送されて反射される。反射された参照光R1は周波数シフタ6−1を再び通過して更なる周波数シフトを受けてビームスプリッタ5−1に入射する。なお、前述のように反射板7−1を移動させることにより参照光R1の周波数をシフトさせてもよい。
ビームスプリッタ5−1を透過した参照光Rは、ビームスプリッタ5−2によりその一部(参照光R2)が反射されて第2の参照系A2に入射する。参照光R2は、周波数シフタ6−2を通過して周波数がシフトされ、反射板7−2により反射され、周波数シフタ6−2を再び通過して更なる周波数シフトを受けてビームスプリッタ5−2に入射する。
同様に、参照光Rは、各ビームスプリッタ5−k(k=3〜n−1)によりその一部(参照光Rk)が反射されて第kの参照系Akに入射され、周波数シフタ6−kにより周波数シフトを受けてビームスプリッタ5−kに再入射する。
ビームスプリッタ5−(n−1)を透過した参照光R(Rn)は、反射鏡5−nにより反射され、第nの参照系Anに入射し、周波数シフタ6−nにより周波数シフトを受けて反射鏡5−nに再入射する。
ここで、周波数シフタ6−1〜6−nは、参照光R1〜Rnに対して各々異なるシフト量だけ周波数を変化させる。各周波数シフタ6−1〜6−nによる周波数のシフト量はあらかじめ設定されている。
参照系A1〜Anを経由した参照光R1〜Rnは、ビームスプリッタ5−1〜5−(n−1)及び反射鏡5−nによって合成され、ビームスプリッタ4に再び入射する。この合成光のことも参照光Rと呼ぶこととする。当該参照光Rは、各々異なる量の周波数シフトを受けた参照光R1〜Rnを周波数成分とする光ビームである。
被測定物体Oの各深さ領域D1〜Dnに対応する成分を含んだ信号光Sと、各参照光R1〜Rnを周波数成分とする参照光Rとは、ビームスプリッタ4によって重畳されて干渉光Lを生成する。このとき、反射板7−1〜7−nは上述のような配置設定であることから、深さ領域Dkに対応する信号光Sの成分と、参照光Rの周波数成分Rkとが干渉することとなる(k=1〜n)。したがって、干渉光Lは、図2に示すように、深さ領域D1に対応する成分と周波数成分R1とからなる干渉光L1、深さ領域D2に対応する成分と周波数成分R2とからなる干渉光L2、・・・・・、深さ領域Dnに対応する成分と周波数成分Rnとからなる干渉光Lnを成分して含んでいる。この干渉光L1〜Lnを「干渉光成分」と称することとする。各干渉光成分L1〜Lnは、上記の式(1)に示したように、干渉光Lの背景光からなる直流成分と、うなりの周波数(ビート周波数)を有する交流成分(ビート信号)とが重畳されたものと考えることができる。
干渉光Lは、ビームスプリッタ8−1によって2つに分割され、その一部がシャッタ9−1を介してCCDカメラ10−1により検出される。
また、ビームスプリッタ8−1を透過した干渉光Lは、次のビームスプリッタ8−2によって更に二分され、その一部がシャッタ9−2を介してCCDカメラ10−2により検出される。
ビームスプリッタ8−2を透過した干渉光Lは、シャッタ9−3を介してCCDカメラ10−3によって検出される。
ここで、ビームスプリッタ8−1による干渉光の分割率、つまり透過される干渉光と反射される干渉光Lとの強度比は2:1であることが望ましい。すなわち、ビームスプリッタ8−1は、入射光の2/3を透過させ、1/3を反射させる特性を有していることが望ましい。また、ビームスプリッタ8−2により透過される干渉光Lと反射される干渉光Lとの強度比は1:1であることが望ましい。すなわち、ビームスプリッタ8−2は、入射光の1/2を透過させ、1/2を反射させる特性を有していることが望ましい。これにより、CCDカメラ10−1、10−2、10−3により検出される干渉光Lは、それぞれ等しい強度レベルとされるので、後述の演算処理を行うのに好適である。ただし、分割される干渉光の強度比はこれに限定されるものではなく、適宜設定することが可能である。
前述の図2と図3〜図5は、シャッタ9−1、9−2、9−3による干渉光Lのサンプリング形態を説明するための図である。図3(A)は、深さ領域D1に対応する成分と周波数成分R1とからなる干渉光成分L1の時間波形を示している。図3(B)、(C)、(D)は、それぞれシャッタ9−1、9−2、9−3の開閉動作を制御するサンプリング関数の一例を示している。これらサンプリング関数の周波数は、干渉光成分L1の周波数に等しいか、もしくはそれに近い値に設定されている。また、周波数シフタ6−1〜6−nは各々異なるシフト量を付与することから、干渉光成分L1〜Lnは各々異なる周波数を有している。したがって、各シャッタ9−1、9−2、9−3は、干渉光Lから干渉光成分L1を選択的に抽出してサンプリングすることができる。
図3(B)に示すサンプリング関数をm(t)とする。シャッタ9−1は、このサンプリング関数m(t)に基づいて周期的に開閉(on−off)させることにより干渉光成分L1のサンプリングを行う。サンプリング関数m(t)は、50%dutyの矩形列からなる波形を有し、その周波数fsmは、上述のように干渉光成分L1の周波数fifに等しいか、もしくはそれに近い値に設定されている(すなわち、fsm=fifもしくはfsm≒fifとされる)。
サンプリング関数m(t)によれば、干渉光成分L1の位相0〜πの範囲がサンプリングされる。ここで、サンプリング関数m(t)の周波数fsmと式(1)に示すような干渉光成分L1のビート周波数fifとの差(δf=|fif−fsm|)は、蓄積型光センサであるCCDカメラ10−1の応答周波数に比べて十分小さく設定されている。それにより、干渉光成分L1の各周期においてほぼ同じ位相範囲がサンプリングされるようになっている。このとき、干渉光成分L1を受光したCCDカメラ10−1からの出力i(t)は、測定時間内にCCDカメラ10−1に蓄積された光電荷量に比例しており、具体的には次式によって与えられる(例えば、M.Akiba、K.P.Chan、N.Tanno、「Optics Letters」、Vol.28、816(2003)を参照)。
Figure 0004409331
ここで、<−>はCCDカメラ10−1の蓄積効果による時間平均を表している。また、φは測定の初期位相値を表し、Kはビームスプリッタ8−1の反射率とCCDカメラ10−1の光電変換率を含めた光検出効率を表している。
同様に、干渉光成分L1は、図3(C)に示すサンプリング関数m(t)に基づいてその開閉タイミングが制御されるシャッタ9−2によってサンプリングが施されてCCDカメラ10−2により検出される。このサンプリング関数m(t)は、サンプリング関数m(t)と同一の周波数fsmを有し、50%dutyの矩形列の波形を有するものである。ここで、サンプリング関数m(t)は、サンプリング関数m(t)に対して位相差Δθ1,2を有する。この位相差Δθ1,2は、位相シフタによる位相のシフト量をあらかじめ設定することによって生成される。以上のような条件の下に、式(2)と同様の原理で、CCDカメラ10−2からは次のような出力i(t)が得られる。
Figure 0004409331
ただし、Kは、ビームスプリッタ8−1の透過率と、ビームスプリッタ8−2の反射率と、CCDカメラ10−2の光電変換率とを含めた光検出効率である。
式(2)と式(3)から分かるように、CCDカメラ10−1、10−2からの出力には、信号光Sと参照光Rの強度I、Iの項とともに、干渉光成分L1の振幅√(I)及び位相(2πδft+φ)、(2πδft+Δθ1,2)に関わる項がそれぞれ含まれている。
更に、干渉光成分L1は、図3(D)に示すサンプリング関数m(t)に基づいてその開閉タイミングが制御されるシャッタ9−3によってサンプリングが施されCCDカメラ10−3により検出される。このサンプリング関数m(t)は、サンプリング関数m(t)と同一の周波数fsmを有し、50%dutyの矩形列の波形を有するものである。ここで、サンプリング関数m(t)は、サンプリング関数m(t)に対して位相差Δθ1,3を有する。この位相差Δθ1,3は、位相シフタによる位相のシフト量をあらかじめ設定することによって生成される。このとき、式(2)と同様の原理で、CCDカメラ10−3からは次のような出力i(t)が得られる。
Figure 0004409331
ただし、Kは、ビームスプリッタ8−1、12のそれぞれの透過率、及びCCDカメラ10−3の光電変換率を含めた光検出効率である。
図4及び図5に示す干渉光成分L2及びLnも、同様の原理によってサンプリングされ、各CCDカメラ10−1、10−2、10−3により検出される。
〔演算処理〕
式(2)、(3)、(4)に示す各CCDカメラ10−1、10−2、10−3から出力された電気信号は、コンピュータ20に送信される。コンピュータ20は、これらの出力結果を用いて以下に説明するような演算を実行することにより、干渉光成分L1の背景光、すなわち干渉光Lの背景光に対応する式(1)に示す直流成分と、干渉光成分L1の信号強度及び位相の空間分布とを算出する。コンピュータ20は、この干渉光成分L1の強度及び位相の分布に基づいて、当該干渉光成分L1に対応する深さ領域D1における被測定物体Oの反射率分布を演算し、この反射率分布に基づいて当該深さ領域D1の2次元断面画像を形成する処理を実行する。なお、干渉光成分L1の強度及び位相情報に基づくコンピュータ20の画像形成処理は、従来の装置と同様の方法で行われるものである。
ここで、図3(B)、(C)、(D)に示すように、サンプリング関数m(t)とサンプリング関数m(t)との位相差Δθ1,2は−π/2に設定され、サンプリング関数m(t)とサンプリング関数m(t)との位相差Δθ1,3はπ/2に設定されているものとする。このとき、干渉光成分L1の背景光からなる直流成分の強度Sと位相直交成分(sin成分とcos成分)S及びSは、それぞれ以下の式により表される。
Figure 0004409331
したがって、3つのCCDカメラ10−1、10−2、10−3のうち2台のCCDカメラ10−2、10−3からの電気信号に基づいて、干渉光Lの背景光に対応する直流成分の強度を算出することができる。
Figure 0004409331
Figure 0004409331
式(6)と式(7)を用いると、干渉光成分L1の振幅は次のように表される。
Figure 0004409331
ここで、右辺に係る比例定数はπ/4である。したがって、式(5)により求めた直流成分を用いることにより、干渉光成分L1の振幅を算出することができる。したがって、上記直流成分とヘテロダイン信号の振幅とを足し合わせることにより、干渉光成分L1の強度を求めることができる。
更に、光画像計測装置1によれば、次のような測定方法により干渉光成分L1の位相の空間分布を求めて画像化することができる。
ある測定時間t=tについて、式(6)及び式(7)に示す干渉光成分L1の干渉成分S(t)及びS(t)が取得されると、これら双方の干渉成分の比を取ることによって次式に示すような信号が得られる。
Figure 0004409331
この式(9)から分かるように、信号Sは干渉光成分L1の振幅には依存せず、その位相情報のみを含んでいる。したがって、2次元光センサアレイであるCCDカメラ10−1、10−2、10−3の各ピクセルにより出力される干渉光成分L1の位相φ(x、y、t)は、次式のように表される。ここで、(x、y)は、これらCCD上に設定された各ピクセルの位置座標を示している。
Figure 0004409331
この式(10)の第2項2πδftは、ゼロあるいはほぼゼロの周波数δfを有する交流信号の測定時間tにおける瞬時位相値であり、CCDカメラ10−1、10−2、10−3のピクセルの位置(すなわち変数x、y)によらず一定であると仮定できる。したがって、CCDカメラ10−1、10−2、10−3上のある座標x=x、y=yに位置するピクセルにおけるヘテロダイン信号の位相φ(x、y、t)に対する、各ピクセルからの信号の位相差を求めることにより、干渉光成分L1の位相の空間分布を画像化することができる。このような干渉光の位相の空間分布測定は、その位相値を基準とする画像計測、例えばヘテロダイン干渉法による鏡状表面の高精度計測に有効であると考えられる。
更に、上記の位相情報を用いると、干渉光成分L1の周波数情報を取得することができる。すなわち、干渉光成分L1に対応するヘテロダイン信号の周波数fifとサンプリング周波数fsmとの周波数差δfは、ある2つの測定時間t=t及びt=tにおける位相の算出結果φ(x、y、t)及びφ(x、y、t)に基づき次式を用いて算出される。
Figure 0004409331
サンプリング周波数fsmは既知であるので、式(11)による算出結果からヘテロダイン信号の周波数fif、つまり干渉光成分L1の周波数を算出することができる。このようなヘテロダイン周波数の測定方法は、例えば被検眼眼底の血流の状態を計測する場合など、ヘテロダイン干渉法を用いたドップラー速度計測に有効に利用可能であると考えられる。
図4及び図5に示す干渉光成分L2及びLnについても、同様の原理によってその強度及び位相の空間分布を求めることができる。
上記実施形態の光画像計測装置1では、本発明の強度変調手段として高速シャッタ等のシャッタ9−1、9−2、9−3を備える構成を採用したが、強度変調手段はこれに限定されるものではない。例えば、干渉光を完全に遮断するこのようなシャッタ手段に代えて、干渉光を透過させる透過率を周期的に増減させることが可能な例えばSLM(Spatial Light Modulator;空間光変調素子)などを設けることにより、CCD等の光検出手段により受光される干渉光の強度を変調させて干渉光をサンプリングすることができる。すなわち、シャッタ手段は、光検出手段に受光される干渉光の強度を0と100(最高強度)との間で切り換えるものであるが、本発明の強度変調手段としては、干渉光の強度を例えば10と80との間で周期的に切り換えるような構成を適用することが可能である。
更に、干渉光の強度変調は、2値間での切り換えだけではなく、3値以上の間を周期的に切り換える方式や、2値間を周期的にかつ連続的に変化させる方式を、サンプリング態様などに応じて採用することが可能である。なお、強度変調の幅は、CCDカメラの感度等を勘案して決定すればよい。本発明の強度変調手段としては、干渉光の強度を周期的に変調させることが可能であれば、いかなる構成を採用してもよい。強度変調手段と光検出手段とを一体的に構成することも可能である。
ビームスプリッタ4、5−1〜5−(n−1)、8−1、8−2としては任意の形態のものを用いることができるが、キューブ型のビームスプリッタを用いると、空気との境界面における反射光がCCDカメラに入射してしまうおそれがあるため、プレート型やウェッジ型などのビームスプリッタを用いることが好ましい。
光画像計測装置1には3つの単体のCCDカメラ10−1、10−2、10−3が設けられているが、例えば3CCDタイプのカラーCCDカメラのような3枚式CCDカメラ(ユニット)を使用し、各CCDチップの前位置に強度変調手段を配設することにより外見上1つのCCDカメラとして構成してもよい。それにより、装置構成の簡略化や装置内の省スペース化などを図ることができる。
また、1つのCCDカメラの受光面を複数の領域に分割し、各領域の前位置に強度変調手段を配置することにより、CCDカメラの各領域を1つのCCDカメラに見立てて干渉光を検出するように構成することもできる。このとき、CCDカメラの当該複数の領域を覆う大きさの液晶SLM等からなる強度変調手段を1つ配置し、CCDカメラの各領域に対応する強度変調手段の領域をそれぞれ制御することにより干渉光の検出を行うようにしてもよい。このような構成によれば、装置構成の簡略化や装置内の省スペース化を図ることができる。また、複数のCCDカメラを互いに同期制御させてサンプリングする必要がなくなるので、制御系の簡略化を図ることも可能となる。
更に、CCDカメラに蓄積された電荷の直流成分のオフセット調整や交流信号のゲイン調整を適宜行うことにより、検出される干渉光が形成する干渉縞のコントラストを向上させるように構成してもよい。
また、参照物体として、2回反射の折り返しミラーやコーナーキューブを適用することにより、参照物体の移動距離を減少させることができる。それにより、装置内の省スペース化を図ることができ、また、参照物体を移動させる上記の駆動装置に掛かる負担や消費電力量を低減することが可能となる。
光画像計測装置1では、サンプリング関数m(t)とサンプリング関数m(t)との位相差Δθ1,2を−π/2に設定し、サンプリング関数m(t)とサンプリング関数m(t)との位相差Δθ1,3をπ/2に設定することにより、各干渉光成分のサンプリングを行うようになっているが、これに限定されるものではない。
また、光画像計測装置1は、ビート周波数と(ほぼ)等しいサンプリング周波数にて干渉光成分のサンプリングを行うように構成されているが、これに限定されるものではない。例えば、サンプリング周波数として、抽出したい干渉光成分のビート周波数の整数倍の周波数を適用すれば、干渉光の複数の位相範囲をそれぞれ周期的にサンプリングすることができる。この方法によれば、干渉光の各周期について複数の位相範囲をサンプリングできるので、干渉光をより詳細に解析することが可能となり、計測精度の向上が期待される。
一方、ビート周波数の整数分の1(1/n)倍のサンプリング周波数を適用することも可能である。この方法によれば、干渉光の所定の位相範囲がn周期おきにサンプリングされるので、強度変調手段による強度の切り換えがビート周波数に追いつかないような場合に有効に利用することができる。
本発明の光画像計測装置に用いられるサンプリング関数は、50%dutyのものを用いることが好ましい。すなわち、50%未満のdutyである場合にはCCDカメラ10−1、10−2、10−3に受光される光量が減少し干渉光の検出効率が低下し、一方、50%を超えるduty比であっても検出効率がダウンしてしまうからである。ただし、50%以外のdutyのサンプリング関数を用いてよいことはいうまでもない。
また、本発明の光画像計測装置に用いられるサンプリング関数としては、シャッタ9−1、9−2、9−3の開閉タイミングを好適に切り換えるために、図3等に示すような矩形の波形を用いることが好ましい。ただし、正弦波や三角波などの矩形以外の波形からなるサンプリング関数を適宜用いることも可能である。
本発明に係る光画像計測装置は、上記実施形態のように干渉光を3つの光路に分割する構成に限定されるものではなく、分割される光路数は任意である。複数の光路に分割する場合には、各光路上にCCD等の光検出手段を設けるとともに、1つの光路を除く全ての光路上に強度変調手段を設ける構成とするか、又は全ての光路上に強度変調手段を設ける構成とすることが好ましい。前者の構成では、強度変調手段が配置されない光路を伝搬する干渉光成分を連続的に受光することにより、その直流成分を求めることが可能である。なお、測定方法や演算方法に応じて、分割された複数の光路のうち任意数の光路上に強度変調手段を設ける構成を適宜採用することができる。
また、本発明に係る光画像計測装置の光検出手段としては、上述したCCDカメラ10−1、10−2、10−3の他にも、積算回路を備えたラインセンサなどを適用することができる。本発明の光検出手段は、干渉光を受光して光電変換する機能と、受光した干渉光に基づく電荷を蓄積する機能との双方を備えている限り、1次元ないし2次元のあらゆる構成のものを用いることが可能である。
以上で説明した本発明の実施形態においては、マイケルソン型の干渉光学系を備えた光画像計測装置について説明したが、例えばマッハツェンダー型などその他の干渉光学系を採用することも当然に可能である(例えば、上記特許文献2参照)。
また、干渉光学系の一部に光ファイバ(バンドル)を設けて導光部材として用いることにより、装置設計上の自由度を高めたり、装置のコンパクト化を図ったり、あるいは、被測定物体の配置の自由度を高めたりすることができる(例えば、上記特許文献2参照)。
本発明の光画像計測装置を例えば眼科の分野に応用すると、被検眼の複数の部分の画像を並列的に取得することができる。例えば、反射板7−1、7−2、7−3を、被検眼の角膜、水晶体、網膜に各々対応する位置に配置させて測定を行うと、干渉光成分L1、L2、L3は、それぞれ角膜、水晶体、網膜の情報を有する信号となる。各干渉光成分に同期した周波数でサンプリングを行うことにより、角膜、水晶体、網膜の2次元断面像をそれぞれ得ることができる。
また、各反射板7−1、7−2、7−3を干渉光成分の可干渉距離分だけ移動させることにより角膜等を深さ方向にスキャンすれば、画像取得の高速化を図ることができる。
また、第1の参照系A1を経由する光路と、第2の参照系A2を経由する光路との光路長差を網膜の厚さ程度に設定するとともに、干渉光成分L1が検出されたら干渉光成分L2の検出に自動的に移行するように制御することもできる。それにより、網膜の前面及び後面の測定を速やかに行うことができる。なお、網膜以外の部分を検出することや、3つ以上の干渉光成分の検出を連動させることも可能である。
本発明の光画像計測装置は、このような使用形態のほかにも、各種の医療的用途や工業的用途に用いることができる。
以上に詳述した構成は、本発明に係る光画像計測装置の実施形態の一例に過ぎないものであり、本発明の要旨の範囲内において各種の変形を施すことができる。
本発明に係る光画像計測装置の構成の一例を示す概略図である。 本発明に係る光画像計測装置により検出される干渉光の時間波形を示すグラフ図である。 本発明に係る光画像計測装置による干渉光の第1の干渉光成分のサンプリング態様の一例を説明するための図である。図3(A)は、第1の干渉光成分の時間波形を示すグラフである。図3(B)は、第1の干渉光成分をサンプリングするための第1のサンプリング関数の波形の一例を示すグラフである。図3(C)は、第1の干渉光成分をサンプリングするための第2のサンプリング関数の波形の一例を示すグラフである。図3(D)は、第1の干渉光成分をサンプリングするための第3のサンプリング関数の波形の一例を示すグラフである。 本発明に係る光画像計測装置による干渉光の第2の干渉光成分のサンプリング態様の一例を説明するための図である。図4(A)は、第2の干渉光成分の時間波形を示すグラフである。図4(B)は、第2の干渉光成分をサンプリングするための第1のサンプリング関数の波形の一例を示すグラフである。図4(C)は、第2の干渉光成分をサンプリングするための第2のサンプリング関数の波形の一例を示すグラフである。図3(D)は、第2の干渉光成分をサンプリングするための第4のサンプリング関数の波形の一例を示すグラフである。 本発明に係る光画像計測装置による干渉光の第nの干渉光成分のサンプリング態様の一例を説明するための図である。図5(A)は、第nの干渉光成分の時間波形を示すグラフである。図5(B)は、第nの干渉光成分をサンプリングするための第1のサンプリング関数の波形の一例を示すグラフである。図5(C)は、第nの干渉光成分をサンプリングするための第2のサンプリング関数の波形の一例を示すグラフである。図5(D)は、第nの干渉光成分をサンプリングするための第3のサンプリング関数の波形の一例を示すグラフである。 従来の光画像計測装置の構成を示す概略図である。 従来の光画像計測装置の構成を示す概略図である。
符号の説明
1 光画像計測装置
2 光源
3 レンズ系
4、5−1〜5−(n−1)、8−1、8−2 ビームスプリッタ
5−n 反射鏡
6−1〜6−n 周波数シフタ
7−1〜7−n 反射板
9−1、9−2、9−3 シャッタ
10−1、10−2、10−3 CCDカメラ
20 コンピュータ
A1〜An 参照系
R、R1〜Rn 参照光
S 信号光
L 干渉光
L1〜Ln 干渉光成分
O 被測定物体
D1、D2、D3 深さ領域

Claims (7)

  1. 光ビームを出射する光源と、
    この光源からの光ビームを平行光束にするとともにそのビーム径を拡大するレンズ系と、
    このレンズ系により平行光束とされてビーム径が拡大された光ビームを基に、被測定物体の深さ領域毎に異なる周波数を有する複数の干渉光成分を含んだ前記拡大されたビーム径を有する干渉光を生成する干渉系と、
    前記生成された干渉光から前記各干渉光成分を選択的に検出する検出系と、
    前記検出された前記各干渉光成分に基づいて、当該各干渉光成分に対応する前記深さ領域における前記被測定物体の2次元断面画像を形成する画像形成手段と、
    を備えていることを特徴とする光画像計測装置。
  2. 前記干渉系は、
    前記被測定物体からの距離が各々異なるように配置された複数の参照物体と、
    前記レンズ系により平行光束とされてビーム径が拡大された光ビームを前記被測定物体の複数の前記深さ領域にて反射される前記拡大されたビーム径を有する信号光と、前記複数の参照物体を各々経由する前記拡大されたビーム径を有する複数の参照光とに分割する分割手段と、
    前記複数の参照光の周波数を各々異なる量だけシフトする周波数シフト手段と、
    前記被測定物体を経由した前記信号光と、前記周波数シフト手段により周波数が各々シフトされ、前記複数の参照物体を各々経由した前記複数の参照光とを重畳して前記干渉光を生成する重畳手段と、
    を含んでいることを特徴とする請求項1に記載の光画像計測装置。
  3. 前記干渉系は、マイケルソン干渉計を用いた構成とされていることを特徴とする請求項1又は請求項2に記載の光画像計測装置。
  4. 前記検出系は、
    前記生成された前記干渉光の強度を前記周波数のシフト量に同期した変調周波数で周期的に変調することにより、当該変調周波数に対応する前記干渉光成分を前記干渉光から抽出する強度変調手段と、
    前記抽出された前記干渉光成分を受光し、電気信号に変換して出力する2次元光センサアレイと、
    を含んでいることを特徴とする請求項2に記載の光画像計測装置。
  5. 前記強度変調手段は、前記干渉光を前記変調周波数にて周期的に遮断するシャッタであることを特徴とする請求項4に記載の光画像計測装置。
  6. 前記検出系は、前記生成された干渉光を複数に分割する干渉光分割手段を更に備え、
    前記強度変調手段及び前記2次元光センサアレイは、前記分割された前記複数の干渉光の各光路上に各々配置されていることを特徴とする請求項4に記載の光画像計測装置。
  7. 光ビームを出射する光源と、
    この光源からの光ビームを平行光束にするとともにそのビーム径を拡大するレンズ系と、
    このレンズ系により平行光束とされてビーム径が拡大された光ビームを、被測定物体の複数の深さ領域にて反射される前記拡大されたビーム径を有する信号光と所定の参照物体により反射される前記拡大されたビーム径を有する参照光とに分割する第1のビームスプリッタと、
    前記参照光を複数に分割する第2のビームスプリッタと、
    前記第2のビームスプリッタにより分割された前記複数の参照光の各光路上に配置され、前記複数の参照光の周波数を各々異なる量だけシフトさせる周波数シフタと、
    前記複数の参照光の各光路上に配置され、前記被測定物体からの距離が各々異なるように配置された前記参照物体としての反射板と、
    前記周波数シフタにより前記周波数が各々シフトされ、前記反射板により各々反射された前記複数の前記参照光と、前記被測定物体の前記複数の深さ領域にて反射された前記信号光とを重畳することにより、前記各深さ領域毎に異なる周波数を有する複数の干渉光成分を含んだ前記拡大されたビーム径を有する干渉光を生成する第3のビームスプリッタと、
    前記生成された前記干渉光を前記周波数のシフト量のいずれかに同期した周波数で周期的に遮断することにより、当該同期した周波数に対応する前記干渉光成分を前記干渉光から抽出するシャッタと、
    前記抽出された前記干渉光成分を受光し、電気信号に変換して出力する2次元光センサアレイと、
    前記出力された電気信号に基づいて前記干渉光成分の強度及び位相の分布を演算し、当該強度及び位相の分布に基づいて前記干渉光成分に対応する前記深さ領域における前記被測定物体の反射率分布を演算し、当該反射率分布に基づいて当該深さ領域における2次元断面画像を形成するコンピュータと、
    を備えていることを特徴とする光画像計測装置。
JP2004100741A 2004-03-30 2004-03-30 光画像計測装置 Expired - Lifetime JP4409331B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2004100741A JP4409331B2 (ja) 2004-03-30 2004-03-30 光画像計測装置
US11/091,478 US7548320B2 (en) 2004-03-30 2005-03-29 Optical image measuring apparatus
EP05006897A EP1582143B1 (en) 2004-03-30 2005-03-30 Optical inteferometer for imaging at several depth regions of an object

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2004100741A JP4409331B2 (ja) 2004-03-30 2004-03-30 光画像計測装置

Publications (2)

Publication Number Publication Date
JP2005283471A JP2005283471A (ja) 2005-10-13
JP4409331B2 true JP4409331B2 (ja) 2010-02-03

Family

ID=34879994

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004100741A Expired - Lifetime JP4409331B2 (ja) 2004-03-30 2004-03-30 光画像計測装置

Country Status (3)

Country Link
US (1) US7548320B2 (ja)
EP (1) EP1582143B1 (ja)
JP (1) JP4409331B2 (ja)

Families Citing this family (53)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4549714B2 (ja) * 2004-03-31 2010-09-22 株式会社トプコン 光画像伝送システム、光画像送信装置、光画像受信装置及び光画像伝送方法
JP4505807B2 (ja) * 2004-08-09 2010-07-21 国立大学法人 筑波大学 多重化スペクトル干渉光コヒーレンストモグラフィー
US7526329B2 (en) * 2004-08-19 2009-04-28 Hogan Josh N Multiple reference non-invasive analysis system
JP2006162366A (ja) * 2004-12-06 2006-06-22 Fujinon Corp 光断層映像装置
US7365856B2 (en) 2005-01-21 2008-04-29 Carl Zeiss Meditec, Inc. Method of motion correction in optical coherence tomography imaging
US7805009B2 (en) 2005-04-06 2010-09-28 Carl Zeiss Meditec, Inc. Method and apparatus for measuring motion of a subject using a series of partial images from an imaging system
ATE516739T1 (de) 2005-12-06 2011-08-15 Zeiss Carl Meditec Ag Interferometrische probenmessung
DE102005058220A1 (de) * 2005-12-06 2007-06-14 Carl Zeiss Meditec Ag Interferometrische Probenmessung
US8009297B2 (en) * 2005-12-07 2011-08-30 Kabushiki Kaisha Topcon Optical image measuring apparatus
JP4890878B2 (ja) * 2006-02-16 2012-03-07 株式会社トプコン 眼底観察装置
JP5406427B2 (ja) * 2006-11-17 2014-02-05 株式会社トプコン 断層画像処理方法、装置およびプログラムならびにこれを用いた光断層画像化システム
EP2102583A2 (en) * 2007-01-19 2009-09-23 The General Hospital Corporation Apparatus and method for controlling ranging depth in optical frequency domain imaging
EP1962080B1 (de) * 2007-02-21 2011-06-01 Agfa HealthCare N.V. System zur optischen Kohärenztomographie
EP2267403A3 (de) * 2007-02-21 2011-04-20 Agfa HealthCare N.V. System und Verfahren zur optischen Kohärenztomographie
EP1962051A1 (de) * 2007-02-21 2008-08-27 Agfa HealthCare N.V. System und Verfahren zur optischen Kohärenztomographie
EP1962049B1 (de) * 2007-02-21 2015-12-23 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
CN101611301B (zh) * 2007-02-28 2012-11-07 日本电信电话株式会社 光反射测定方法以及装置
KR100785802B1 (ko) * 2007-05-29 2007-12-13 (주) 인텍플러스 입체 형상 측정장치
JP5133614B2 (ja) * 2007-06-22 2013-01-30 株式会社ブリヂストン 3次元形状測定システム
JP5117787B2 (ja) * 2007-08-13 2013-01-16 株式会社トプコン 光画像計測装置
JP5154868B2 (ja) * 2007-09-10 2013-02-27 テルモ株式会社 画像診断装置及びその作動方法
JP5339828B2 (ja) * 2007-10-04 2013-11-13 キヤノン株式会社 光干渉断層撮像装置及び光干渉断層撮像方法
US7800759B2 (en) 2007-12-11 2010-09-21 Bausch & Lomb Incorporated Eye length measurement apparatus
KR100939538B1 (ko) * 2007-12-14 2010-02-03 (주) 인텍플러스 입체 형상 측정 장치
DE102008017390A1 (de) * 2008-04-02 2009-10-08 Carl Zeiss Surgical Gmbh Verfahren und Vorrichtung zum Aufnehmen von Farbbildern und LDI-Signalen von einem Objekt
JP5324839B2 (ja) * 2008-06-19 2013-10-23 株式会社トプコン 光画像計測装置
JP2010019630A (ja) * 2008-07-09 2010-01-28 Tokyo Institute Of Technology 顕微分光装置
EP2199734B1 (de) * 2008-12-16 2016-03-09 Agfa HealthCare NV Verfahren und System zur optischen Kohärenztomographie
JP5618533B2 (ja) * 2008-12-26 2014-11-05 キヤノン株式会社 光干渉断層情報取得装置、撮像装置及び撮像方法
JP5339934B2 (ja) * 2009-01-22 2013-11-13 キヤノン株式会社 光断層撮像装置および光断層撮像方法
JP5623028B2 (ja) * 2009-01-23 2014-11-12 キヤノン株式会社 光干渉断層画像を撮る撮像方法及びその装置
JP5725697B2 (ja) * 2009-05-11 2015-05-27 キヤノン株式会社 情報処理装置および情報処理方法
TWI417534B (zh) * 2010-01-27 2013-12-01 私立中原大學 Surface and internal interface of the contrast and measurement device
JP5569963B2 (ja) * 2010-04-15 2014-08-13 国立大学法人宇都宮大学 表面形状測定装置
US9033510B2 (en) 2011-03-30 2015-05-19 Carl Zeiss Meditec, Inc. Systems and methods for efficiently obtaining measurements of the human eye using tracking
US8857988B2 (en) 2011-07-07 2014-10-14 Carl Zeiss Meditec, Inc. Data acquisition methods for reduced motion artifacts and applications in OCT angiography
NL2009273A (en) * 2011-08-31 2013-03-04 Asml Netherlands Bv Level sensor arrangement for lithographic apparatus, lithographic apparatus and device manufacturing method.
DE102011085599B3 (de) * 2011-11-02 2012-12-13 Polytec Gmbh Vorrichtung und Verfahren zur interferometrischen Vermessung eines Objekts
JP5905711B2 (ja) * 2011-11-25 2016-04-20 株式会社トプコン 光画像計測装置
US8934103B2 (en) * 2011-12-22 2015-01-13 General Electric Company Quantitative phase microscopy for label-free high-contrast cell imaging
US9101294B2 (en) 2012-01-19 2015-08-11 Carl Zeiss Meditec, Inc. Systems and methods for enhanced accuracy in OCT imaging of the cornea
JP2013160699A (ja) * 2012-02-08 2013-08-19 Hitachi High-Technologies Corp 光断層画像測定装置
JP6021384B2 (ja) * 2012-03-30 2016-11-09 キヤノン株式会社 光干渉断層撮影装置及び制御方法
WO2014070082A1 (en) * 2012-10-29 2014-05-08 General Electric Company Quantitative phase microscopy for label-free high-contrast cell imaging
KR101357120B1 (ko) * 2013-04-23 2014-02-05 김태민 광 신호를 이용한 거리 측정 방법 및 장치
NL2013701B1 (en) * 2014-10-29 2016-10-04 Cassini B V Interferometer, in particular for optical coherence tomography, comprising a reference arm having optical elements in a fixed positional relationship.
HUP1500616A2 (en) * 2015-12-14 2017-07-28 Pecsi Tudomanyegyetem Data collecting and processing method and arrangement for soft-tomographic examinations
JP6790479B2 (ja) * 2016-06-15 2020-11-25 富士ゼロックス株式会社 眼球の光計測装置
JP7406219B2 (ja) * 2018-09-11 2023-12-27 株式会社トーメーコーポレーション 眼科装置
CN110338756A (zh) * 2019-05-24 2019-10-18 广州永士达医疗科技有限责任公司 一种适用于多光路的oct成像装置以及oct成像系统
IT201900023229A1 (it) * 2019-12-06 2021-06-06 Adige Spa Procedimento e sistema per la determinazione della posizione di un elemento di un sistema ottico in un complesso di lavorazione o misura di un oggetto tramite misure interferometriche parallele
IL273779A (en) * 2020-04-02 2021-10-31 Yehuda Hai Vidal Method and facility for mapping and segmentation based on coherence time comparison

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3044183A1 (de) * 1980-11-24 1982-06-24 Reinhard Dipl.-Phys. Dr. 7250 Leonberg Ulrich Verfahren zur optischen messung von laengen und laengenaenderungen und anordnung zur durchfuehrung des verfahrens
JPH0287005A (ja) * 1988-09-26 1990-03-27 Brother Ind Ltd 光ヘテロダイン干渉表面形状測定装置
US5471303A (en) * 1994-04-29 1995-11-28 Wyko Corporation Combination of white-light scanning and phase-shifting interferometry for surface profile measurements
US5589938A (en) * 1995-07-10 1996-12-31 Zygo Corporation Method and apparatus for optical interferometric measurements with reduced sensitivity to vibration
US5909279A (en) * 1997-03-17 1999-06-01 Hughes Electronics Corporation Ultrasonic sensor using short coherence length optical source, and operating method
US6198540B1 (en) * 1997-03-26 2001-03-06 Kowa Company, Ltd. Optical coherence tomography have plural reference beams of differing modulations
US6088100A (en) * 1997-07-14 2000-07-11 Massachusetts Institute Of Technology Three-dimensional light absorption spectroscopic imaging
FR2783323B1 (fr) * 1998-09-10 2000-10-13 Suisse Electronique Microtech Dispositif interferometrique pour relever les caracteristiques de reflexion et/ou de transmission optiques en profondeur d'un objet
JP3423229B2 (ja) * 1998-11-17 2003-07-07 株式会社ミツトヨ 光波干渉計及び光波干渉計を用いた測長方法
US6195168B1 (en) * 1999-07-22 2001-02-27 Zygo Corporation Infrared scanning interferometry apparatus and method
JP3245135B2 (ja) 1999-08-26 2002-01-07 科学技術振興事業団 光計測装置
JP3594875B2 (ja) 2000-05-25 2004-12-02 独立行政法人 科学技術振興機構 2次元光ヘテロダイン検出法を用いた光画像計測装置
US6924898B2 (en) * 2000-08-08 2005-08-02 Zygo Corporation Phase-shifting interferometry method and system
DE10045535A1 (de) * 2000-09-13 2002-04-04 Gerhard Bonnet Optische Vorrichtung
US20020044713A1 (en) * 2000-10-16 2002-04-18 Henry Hung Multiple wavelength michelson interferometer
US6728571B1 (en) 2001-07-16 2004-04-27 Scimed Life Systems, Inc. Electronically scanned optical coherence tomography with frequency modulated signals
DE10150934A1 (de) * 2001-10-09 2003-04-10 Zeiss Carl Jena Gmbh Verfahren und Anordnung zur tiefenaufgelösten Erfassung von Proben
US6882431B2 (en) * 2002-11-26 2005-04-19 Malvin C. Teich Quantum optical coherence tomography data collection apparatus and method for processing therefor
JP4381847B2 (ja) * 2004-02-26 2009-12-09 株式会社トプコン 光画像計測装置
JP4522724B2 (ja) * 2004-03-16 2010-08-11 株式会社トプコン 光画像計測装置
JP4549714B2 (ja) * 2004-03-31 2010-09-22 株式会社トプコン 光画像伝送システム、光画像送信装置、光画像受信装置及び光画像伝送方法
US7388668B2 (en) * 2004-06-09 2008-06-17 President & Fellows Of Harvard College Phase sensitive heterodyne coherent anti-Stokes Raman scattering micro-spectroscopy and microscopy

Also Published As

Publication number Publication date
JP2005283471A (ja) 2005-10-13
US7548320B2 (en) 2009-06-16
US20050219544A1 (en) 2005-10-06
EP1582143A1 (en) 2005-10-05
EP1582143B1 (en) 2012-09-19

Similar Documents

Publication Publication Date Title
JP4409331B2 (ja) 光画像計測装置
JP4381847B2 (ja) 光画像計測装置
JP4409332B2 (ja) 光画像計測装置
JP4522724B2 (ja) 光画像計測装置
JP4409384B2 (ja) 光画像計測装置及び光画像計測方法
JP4633423B2 (ja) 光画像計測装置
JP4494160B2 (ja) 光画像計測装置
JP4566685B2 (ja) 光画像計測装置及び光画像計測方法
JP4837300B2 (ja) 光画像計測装置
JP4563130B2 (ja) 光画像計測装置
US8009297B2 (en) Optical image measuring apparatus
JP4602372B2 (ja) 光画像計測装置及び光画像計測方法
US8970848B2 (en) Optical image measurement apparatus
JP2010164574A (ja) 多重化スペクトル干渉光コヒーレンストモグラフィー
JP4409334B2 (ja) 光画像計測装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20070328

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

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090818

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20091016

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

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

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

Ref document number: 4409331

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20121120

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20121120

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20131120

Year of fee payment: 4

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

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