(57)【要約】【課題】分散の影響を減少させ、リフレクトメトリーとOCT測定の分解能を向上させる。【解決手段】高干渉光の波長を連続的または離散的に変化させた光を参照光路と試料光路とに分波し、参照光路を経由した光と試料光路を経由した光とを合波し、合波された光の強度を検出し、波長ごとの光強度信号データを解析して試料の光軸方向の位置と後方散乱強度との関係を求め、試料の光軸方向に１から１０００ミクロン程度に弁別した微小領域ごとに後方散乱光の強度を算出する際、使用する光学部品や試料の屈折率の波長依存性の試料光路と参照光路における差を光学系及び試料の光学的光路長の波長依存性データとして記憶するデータ記憶部３４と、検出された波長ごとの光強度信号データを、位置を変数とする関数から位置を変数とする関数に信号変換する際に光学系及び試料の光学的光路長の波長依存性データに基づいて屈折率の波長依存性を補正する変換処理部３６を備える。【選択図】図６
2005283155.tif 000002
【特許請求の範囲】
【請求項１】
試料の光軸方向の長さに比べて長いコヒーレント長を有する高干渉光を発生し、当該高干渉光の波長を連続的または離散的に変化させる光発生手段と、
前記光発生手段から発生された光を、参照光路と前記試料に導く試料光路とに分波する光分波手段と、前記参照光路を経由した光と前記試料光路を経由した光とを合波する光合波手段と、
前記光合波手段によって合波された光の強度を検出する光強度検出手段と、
前記光強度検出手段により得られた波長ごとの光強度信号データを解析して前記試料の光軸方向の位置と後方散乱強度との関係を求め、前記試料の光軸方向に１から１０００ミクロン程度に弁別した微小領域ごとに後方散乱光の強度を算出する後方散乱光強度解析手段とを有する装置において、
使用する光学部品や試料の屈折率の波長依存性の試料光路と参照光路における差を、光学系及び試料の光学的光路長の波長依存性データとして記憶する記憶手段と、
前記検出された波長ごとの光強度信号データを、位置を変数とする関数に信号変換する際に、前記光学系及び試料の光学的光路長の波長依存性データに基づいて屈折率の波長依存性を補正する信号処理手段と、
を備えたことを特徴とする光干渉断層像撮像法における分散補正装置。
【請求項２】
波長の広がりを持った短コヒーレント光を発生する光発生手段と、
前記光発生手段から発生された光を、参照光路と試料に導く試料光路とに分波する光分波手段と、前記参照光路を経由した光と前記試料光路とを経由した光を合波する光合波手段と、
合波された光を波長ごとの光強度信号に分光する分光手段と、
前記分光手段により分光された波長ごとの光強度信号データを解析して前記試料の光軸方向の位置と後方散乱強度との関係を求め、光軸方向に１から１０００ミクロン程度に弁別した微小領域毎に後方散乱光強度を算出する後方散乱光解析手段において、
使用する光学部品や試料の屈折率の波長依存性の試料光路と参照光路における差を、光学系及び試料の光学的光路長の波長依存性データとして記憶する記憶手段と、
前記検出された波長ごとの光強度信号データを、位置を変数とする関数に信号変換する際に、前記光学系及び試料の光学的光路長の波長依存性データに基づいて屈折率の波長依存性を補正する信号処理手段と、
を備えたことを特徴とする光干渉断層像撮像法における分散補正装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、光干渉計を用いて生体などの組織や物質などの試料の、光軸に沿った奥行き方向の後方反射（試料については、反射には散乱も含まれるものとする。以下同様）位置を測定するリフレクトメーターや、試料に当てる光線の位置を走査しリフレクトメトリーを行うことによって試料の断層像を撮像する装置において、光干渉計に用いる光学部品や試料の屈折率の波長依存性によって生じる装置の分解能の低下を、計算処理によって補正し、分解能を改善した光干渉断層像撮像法における分散補正装置に関する。
【背景技術】
【０００２】
オプティカル・コヒーレンス・トモグラフィー（Optical Coherence Tomography 略してOCTと呼ばれる）技術を用いて生体などの光干渉断層像を撮像する方法としては、現在、主としてオプティカル・コヒーレンス・ドメイン・リフレクトメトリー（Optical Coherence Domain Reflectometry略してOCDR）法が採用されている。以下この用法をOCDR−OCTと呼ぶ。
【０００３】
図１１は、従来より使用されているOCDR−OCTの概略構成を示している。同図を参照してオプティカル・コヒーレンス・ドメイン・リフレクトメトリーについて説明する。同図では、干渉計としてマイケルソン干渉計を使用している。光源としてスペクトル幅の広い、従ってコヒーレンス（干渉性）の悪い低コヒーレンス光源１を用い、ハーフミラー２で光路を参照光路４と、試料光路１２の2つの光路に分割する。参照光路４を進んだ光は、レンズ１１でコリメートされ、参照ミラー３で反射され、同じ光路をハーフミラー２まで戻る。試料光路１２を進んだ光は、試料８に入射し、試料８から後方反射（後方散乱も含む、以下同様）された光は再び試料光路を通ってはハーフミラー２に戻り、参照光路４から戻った光とハーフミラー２で重ね合わされ、干渉する。干渉光の強度は、検出器１３で検出され、アンプ１４で増幅され、A/Dコンバータ（アナログ・デジタル・コンバータ）１５でデジタル信号とし、コンピュータ１６に記憶される。コンピュータ１６に記憶されたデータは、断層像を構成し表示するための計算処理が行われる。結果は、出力／表示手段２８によって、断層画像として表示され、プリンターなどによって出力される。
【０００４】
低コヒーレンス光源１の場合、図１１に示したようなマイケルソン干渉計においては、参照光路４と試料光路１２の光路長が一致するときだけ干渉光が観測される。そこで、参照ミラー３を光軸方向に変位させると、参照光路４と試料光路１２の光路長が一致した試料の奥行き（光軸）方向（縦方向とする）の場所のみ選択して観測できる。試料光路１２の光軸上に、反射面が複数あれば、それらの位置と反射強度が、信号が観測されるときの参照ミラー３の位置から決定できる。これが、オプティカル・コヒーレンス・ドメイン・リフレクトメトリー（OCDR）の原理である。
【０００５】
そして、レンズ９と１０及びガルバノミラー７とを用いて試料の光軸方向と直角方向（横方向とする）の測定位置を選び、走査させながら、オプティカル・コヒーレンス・ドメイン・リフレクトメトリーをそれぞれの測定位置で行えば、試料の奥行き方向とそれに垂直な方向の２次元的な試料の反射面の断層像（トモグラフィー）が得られる。これが、OCDR−OCTである。
【０００６】
OCDR−OCTは、図１２のように、ハーフミラーの替わりにカプラ１７を用い、ファイバー光学系で構成することもできる。
【０００７】
しかし、分散があるとリフレクトメトリーやOCTの分解能が低下する。分散による分解能の低下について次に説明する。図１１の場合も図１２の場合も、光学部品の屈折率が波長とともに変化する（これを、分散があるという）。この分散の大きさは、参照光路４と試料光路１２とでは同じではない。また、試料の中でも分散がある。これは、試料光路１２のみで起こる分散である。これらの分散があると、物理的には波長によって変わらないミラーや反射面の位置が、光学的には波長によって異なって見えてしまう。分散がなければOCDR−OCTは、低コヒーレンス光源１のスペクトル幅が広いほど分解能が上がることが知られている。しかし、分散があると、物理的には同じ参照光路４と試料光路１２との光路差が、光学的には波長によって差が生じるので、波長の変化とともに、試料８と参照ミラー３の相対位置が光学的には波長とともにずれていくかのように見えてしまい、低コヒーレンス光源１のスペクトル幅がある程度以上広がると、分解能がかえって悪くなるという現象が見られる（非特許文献１）。したがって、OCDR−OCTの分解能が、分散で制限されてしまう。
【０００８】
従来、このような分散を補償する方法としては、参照光路４または試料光路１２のどちらかに、分散を補償する光学部品を入れていた。この方法を高分解のOCDR−OCTを実現した文献２の例を図１２に示す。この場合、光ファイバー干渉光学系を用いており、溶融石英５と光学ガラス６（材質はBK7）を参照光路４に入れて、分散を補償している。それぞれの役目は、溶融石英５は光学的な光路長の差を主として補償し、光学ガラス６は両光路に用いられている光学部品の屈折率の２次の波長依存性の差(光の群速度の差)を主として補償する。２つのくさび型のプリズム形状にした部品を向かい合わせて互いにスライドし、重ね合う部分の厚さを増減して、補償の程度を増減する。この光学的な補償によって、文献２では、高分解のOCDR−OCTが実現でき、生体の細胞の細胞核の観測に成功している。
【０００９】
図１１の光学系においても、参照光路４に溶融石英５と光学ガラス６（材質はBK7）を挿入すれば、図１２の場合と同様に、分散を補正することが出来る。
【非特許文献１】C. K. Hitzenberger, A. Baumgartner, W. Drexler, and A.F.Fercher “Dispersion effects in partial coherence interferometry :  Inplications for intraocular ranging ” , Journal of Biomedical Optics vol.4 (1999) 144-151.
【非特許文献２】W. Drexler, U. Morgner, F. X. Kartner, C. Pitris, S. A. Boppart, X. D. Li, E. P. Ippen, and J. G. Fujimoto, “In vivo ultrahigh-resolution optical coherence tomography” Optics Letters vol.24 (1999) 1221-1223.
【非特許文献３】C. K. Hitzenberger and A. F. Fercher, “Alternative OCT Technique”, in Handbook of optical coherence tomography, edts. B. E. Bouma and G. J. Tearney, 359-383 (Marcel Dekker, New York,2002).
【発明の開示】
【発明が解決しようとする課題】
【００１０】
しかしながら、上述した従来方法では、分散の補償のための光学部品を用いているが、この場合、装置の部品数が増して構造が複雑となり、製造コストも上昇する。装置や試料を変更するたびに、分散補償の光学部品を調整する必要もある。また、この補償方法は、試料全体に対して連続的で、試料が層状構造をしていたり、不均一な構造をしていて複数の光学的に異なる部位で構成されていて、異なる部位間で光の分散特性が異なる場合、その光学特性の補償は困難であった。
【００１１】
OCTにおいて、干渉計を構成する光学部品や試料の分散を補償するために、光学部品を付加して行うことは、システムが複雑化し、製造コストが増し、システムに変更を加える度に光学部品の調整が必要となる。
【００１２】
本発明は、上記課題に鑑みてなされたものであり、システムに分散補償用の光学部品を付加することなく、または、システムに分散補償用の光学部品を取り付けても補償しきれない分散を、分散を補正するデータ解析によって分散の影響を減少させ、リフレクトメトリーとOCT測定の分解能を向上させることを可能にした光干渉断層像撮像法における分散補正装置を提供することを目的としている。
【課題を解決するための手段】
【００１３】
上記の目的を達成するために本発明は、請求項１では、試料の光軸方向の長さに比べて長いコヒーレント長を有する高干渉光を発生し、当該高干渉光の波長を連続的または離散的に変化させる光発生手段と、前記光発生手段から発生された光を、参照光路と前記試料に導く試料光路とに分波する光分波手段と、前記参照光路を経由した光と前記試料光路を経由した光とを合波する光合波手段と、前記光合波手段によって合波された光の強度を検出する光強度検出手段と、前記光強度検出手段により得られた波長ごとの光強度信号データを解析して前記試料の光軸方向の位置と後方散乱強度との関係を求め、前記試料の光軸方向に１から１０００ミクロン程度に弁別した微小領域ごとに後方散乱光の強度を算出する後方散乱光強度解析手段とを有する装置において、使用する光学部品や試料の屈折率の波長依存性の試料光路と参照光路における差を、光学系及び試料の光学的光路長の波長依存性データとして記憶する記憶手段と、前記検出された波長ごとの光強度信号データを、位置を変数とする関数に信号変換する際に前記光学系及び試料の光学的光路長の波長依存性データに基づいて屈折率の波長依存性を補正する信号処理手段とを備えたことを特徴としている。
【００１４】
また、請求項２では、波長の広がりを持った短コヒーレント光を発生する光発生手段と、前記光発生手段から発生された光を、参照光路と試料に導く試料光路とに分波する光分波手段と、前記参照光路を経由した光と前記試料光路とを経由した光を合波する光合波手段と、合波された光を波長ごとの光強度信号に分光する分光手段と、前記分光手段により分光された波長ごとの光強度信号データを解析して前記試料の光軸方向の位置と後方散乱強度との関係を求め、光軸方向に１から１０００ミクロン程度に弁別した微小領域毎に後方散乱光強度を算出する後方散乱光解析手段において、使用する光学部品や試料の屈折率の波長依存性の試料光路と参照光路における差を、光学系及び試料の光学的光路長の波長依存性データとして記憶する記憶手段と、前記検出された波長ごとの光強度信号データを、位置を変数とする関数に信号変換する際に、前記光学系及び試料の光学的光路長の波長依存性データに基づいて屈折率の波長依存性を補正する信号処理手段とを備えたことを特徴としている。
【発明の効果】
【００１５】
本発明によれば、光干渉計を用いた生体のOFDR法またFD法に基づくリフレクトメーターとOCT装置において、測定に用いる光学部品や試料の分散のために分解能が低下する現象を、計算処理で分散を補正する解析を行い、システムの単純化、製造コストの削減、システム変更の容易さ、試料の各部位に対して異なる補正が実現できる。
【発明を実施するための最良の形態】
【００１６】
＜発明の原理的説明＞
先ず、OCTについて説明する。現在主流として研究が行われ、また、製品が販売されているOCTは、前述のOCDR−OCTである。OCTを実現する方法としては、この他に、オプティカル・フリーケンシー・ドメイン・リフレクトメトリー（Optical Frequency Domain Reflectometry以下略してOFDR）法とフリーケンシー・ドメイン（Frequency Domain 以下略して FD）法がある（非特許文献３）。フリーケンシー・ドメイン法については、まだ、学会で異なった名称で呼ばれているが、後に説明する図４に示す装置と同等または類似なものを意味するものとする。
【００１７】
これらのOFDR−OCTやFD−OCTを、ここでは、まとめて周波数領域法OCTと呼ぶことにする。周波数領域法OCTでは、干渉信号を波長ごとに別々に標本化（サンプリング）し、得られたデータにフーリエ変換などの計算処理を行って、試料の奥行き方向のそれぞれの位置の信号を得て、断層像を構成する。信号が、波長ごとにサンプリングされているため、波長に依存した屈折率の影響をそれぞれのサンプリングした波長毎に、信号処理の計算過程で補正することができる。
【００１８】
OFDR−OCTの原理の概略を説明する。OFDR−OCTの概略を図１に示す。この干渉計は、ファイバー光学系で構成されるマイケルソン干渉計である。光源に波長可変レーザ１８を用いる。波長可変レーザ１８からの出力光は、カプラ１７で参照光路４と試料光路12に分波される。参照光路４を進んだ光は、参照ミラー３で反射され、カプラ17に戻る。参照光路４の参照ミラー３は、OCDR−OCTの場合と異なり、固定されている。試料光路１２に進んだ光も、試料８の各部位で後方反射されカプラ１７に戻る。試料光路１２と参照光路４から戻り、カプラ１７で合波される光は干渉し、干渉の条件に従って、強めあったり弱めあったりする。
【００１９】
干渉信号について説明する。図１において、試料８が単一のミラーで、ガルバノミラー７を動かさない単純な場合を例として説明する。カプラ１７で合波され干渉した光の強度は、次式で与えられる。
【数１】
2005283155.tif 000003
【００２０】
ここで、ki=2π/λiは波数を表し、添え字ｉはレーザから出力されるｉ番目の光の波長λiに対応することを示す。Ｉr0は参照信号の強度を、Ｉs0は試料８からの信号強度を、ΔLは試料光路１２と参照光路４との光学的光路差を示す。この式の右辺第３項が干渉信号（インターフェログラムとも呼ばれる）である。
【００２１】
干渉した光の強度を検出器１３で検出し、アンプ１４で増幅してA/Dコンバータ１５でデジタル化し、信号処理部２９に取り込み、記憶させる。（１）式の右辺の第１項と第２項は、直流的な信号で、電気信号に変換したのち、電子回路のフィルターで取り除いたり、オフセットをかけた差動アンプで取り除くことが出来る。また、後に説明する図５のように、干渉させるカプラ２６に、参照光と信号光を差動的に入力し、差動アンプ２７を用いて検出すると、差動アンプの出力は、（１）式の第３項の干渉信号に比例した電気信号になる。従って、以下では、この干渉信号の処理について述べる。
【００２２】
〈干渉信号の処理〉
干渉信号の波数依存性と光路差ΔLとの関係について説明する。干渉信号は、波長可変レーザ１８の波数kiを連続的または不連続的に変化させると、干渉条件もそれに伴って変化し、（１）式右辺第３項のように、余弦関数（cos関数）にしたがって強度が増減する。（１）式のように試料が単一のミラーの場合は、波数kiの関数として、参照光路4と試料光路１２の光路差ΔLに逆比例する周期で変動する正弦波の干渉信号が、図３の点線のように得られる。光路差ΔLが大きいほど、点線の正弦波は早く振動する。この振動の周期から、光路差ΔLを求めることができる。
【００２３】
図３に点線で示す波数を変数とする干渉信号の周期を、試料の奥行き方向の位置を変数とする信号に変換する方法について説明する。干渉信号を図１のA/Dコンバーター１５でサンプリングすると、サンプリング間隔で離散化された干渉信号データが得られる。信号処理部２９を図６に示す。周波数領域方式のリフレクトメーターまたはOCT装置において、サンプリングされ離散化された干渉信号は、元データ３０として記憶装置３９に記憶する。この元データ記憶部３０に記憶された元データを、波数変数信号／位置変数信号変換処理部３６で処理し、位置変数信号を求める。
【００２４】
《光学部品と試料の分散が無視できる場合》
光学部品と試料の分散が無視できる場合の、波数変数信号／位置変数信号変換処理部３６の処理方法の一例について説明する。
【００２５】
（１）式の右辺第３項の干渉信号をサンプリングし、元データ記憶部３０に記憶された元データに対して、以下に示すように、離散的なフーリエ解析を行うと、光路差ΔLが求められ、試料光路におけるミラーの位置が求められる。
【００２６】
波数kiにおける元データに、cos(kiz)を乗じてi=lからＮまでのすべてのkiについて和をとったFc(z)と、sin(kiz)を乗じてすべてのkiについて和をとったFs(z)とを次式のように求める。
【数２】
2005283155.tif 000004
【００２７】
ここでＮは、サンプリングされたデータの総数である。これらの計算は、離散的フーリエ変換で、高速に変換する計算アルゴリズムでこれを行う場合は、高速フーリエ変換(FFT)という。（２）式、（３）式は、コンピューターで数値計算すればよいわけであるが、本発明の原理を分かりやすくするために、三角関数の公式を用いて、変形を次のようにする。
【数３】
2005283155.tif 000005
【００２８】
これらの式の２乗の和Ft(z)2を次式のように求める。
【数４】
2005283155.tif 000006
【００２９】
これらの計算は、図６の信号処理部29において、波数変数信号／位置変数信号変換処理部３６で行う。波数を等間隔Δkで走査し、ki=k0＋iΔkとすると、（６）式の計算結果は、次式となる。
【数５】
2005283155.tif 000007
【００３０】
ここで、関数B(z)は、右辺の第１項や第２項に比べて小さく、測定にとっては小さなノイズとなる。（７）式右辺第１項は（４）、（５）式の右辺第１項の和に対応し、（７）式右辺第２項は（４）、（５）式の右辺第２項の和に対応する。
【００３１】
ここで、Ｗkは、走査する波数の全範囲で、次式で与えられる。
【数６】
2005283155.tif 000008
【００３２】
ｚの関数としてFt(z)2を図示すると、図７のようになる。z=2πm/Δk±ΔL（mは整数）で周期的に信号のピークが見られる。周期は2π/Δkである。使用する信号は、この中の１つのピークでよい。例えば、通常は、z=ΔLで観測される信号を用いる。このピークが観測されるｚの値から、光路差ΔLが求まり、参照光路4の参照ミラー3の位置に対する、試料光路12における試料の位置8が求められる。
【００３３】
周期的に観測されるピークが重ならない距離の範囲が望ましい観測可能な光路差ΔLの範囲で、＋方向と−方向にそれぞれπ/Δkで、全体として2π/Δkである。
【００３４】
各ピークの分解能は、ピークの半値全幅（full width at half maximumを略してFWHMと表す）で決まる。全ての波数で強度が一定の場合、この分解能ΔZrは、試料の厚さ（観測される光路差の半分）で見て、次式で与えられる。
【数７】
2005283155.tif 000009
【００３５】
ただし、この分解能は 距離を変数とする関数として、Ft(z)2の替りにFt(z)を図示することがある。この図示の場合のピークのFWHMは、ΔZr=3.79/Ｗkである。
【００３６】
次に、リフレクトメーターについて説明する。
【００３７】
これまでの式は、光学的光路差がΔLの試料の一点からの反射光の場合を式で示した。試料８が、例えば図８（a）に示すように、３枚の透明な板状の物体が重なったもので、それぞれの厚さが、d1，d2，d3とすると、反射面はA，B，C，Dの４面あり、それぞれの面からの反射光が、（５）式のように観測される。ただし、光路差は、最初の面AをΔLとすると、光は往復するので、面B，C，DではそれぞれΔL＋2d1，ΔL＋2(d1＋d2)，ΔL＋2(d1＋d2＋d3)となる。この場合、図７でz=＋ΔLの位置に観測されるピークは４つに分かれ、図８（ｂ）のような信号になる。これらのピークの位置から試料の反射面の位置が、ピークの強さから、反射の強度が測定できる。これがOFDRリフレクトメトリーであり、リフレクトメトリーを行う装置がリフレクトメーターである。
【００３８】
OCTについて説明する。図１のレンズ９、レンズ１０、ガルバノミラー７の光学系の部分において、ガルバノミラー７を回転すると、試料光路１２の光線が試料８に当たる位置が変わる。これを順次走査して、その都度、図８（ｂ）に示すような反射強度の距離依存性を測定し、試料の奥行き方向（縦方向）と横方向の関数として強度を図示すると、図８（a）のような配置構成で、層の間の境界の強度が反射強度に対応するような、断層像が得られる。このような断層像を測定する方法がOCTである。
【００３９】
生体などの組織の場合、後方反射面（散乱部位）は連続的に変化し、それぞれの部分からの反射光がある。この場合、（５）式の信号がｚの関数として連続的に重なって観測される。これが、奥行き方向の、試料の反射光を示す。この測定を、ガルバノミラー７で方向を変えながら測定すれば、生体のOCT像が得られる。
【００４０】
《光学部品や試料に分散がある場合》
次に、光学部品や試料に分散がある場合の、リフレクトメトリーとOCTへの影響について説明する。光学的光路長は、物理的な光路長に屈折率を乗じたものになる。光路差ΔLは、参照光路４と試料光路１２の光学的な光路の差である。波数kiを等間隔で変化させても、波数kiの変化とともに光学的光路差ΔLも変化してしまうので、干渉信号は、図３の点線のような単純な正弦波でなく、実線のように周期が波数とともに変わることになる。この実線は、波数の増加とともに、光路差が小さくなる場合を示している。
【００４１】
波数がk0の時の光学的光路差をΔL0とし、波数がki=k0＋iΔkの時の、分散による光路差の変化をΔLiとし、光路差をΔL0＋ΔLiとすると、分散がある場合の（１）式右辺第３項の干渉信号は次式となる、
【数８】
2005283155.tif 000010
【００４２】
この信号に（２）−（６）に示したような変換を行っても、（７）式のような鋭いピークは得られない。（７）式の各ピークがΔLiの程度広がったピークとなる。
【００４３】
次に、分散の大きさの実験による評価について説明する。図３において、分散がある場合の実線で示す干渉信号は、分散がない場合は点線のような周期的関数になるはずである。波数kiは分かっているので、観測された実線の干渉信号が（１０）式に最適に適合するようにΔL0とΔLiを決定すれば、分散による光学的光路差の変化ΔLiが求められる。これは、図６の波長依存性データ記憶部３４内の光学系及び試料の光学的光路長の波長依存性のデータの一例を示している。光学系及び試料の光学的光路長の波長依存性のデータの求め方については、後述する。
【００４４】
周波数領域方式リフレクトメーターやＯＣＴ装置の場合、光学系及び試料の光学的光路長の波長依存性のデータが定量的に決定できると、光学部品などハードの部品を使わなくても、図６に示す信号処理によって、分散の影響を取り除いた反射強度の距離依存性の信号を求めることができる。
【００４５】
図６の信号処理部について説明する。周波数領域方式リフレクトメーターまたはOCT装置４６の測定データを、元データ３０として記憶装置３９に取り込む。この元データは、各試料位置での図３に示す干渉信号である。干渉信号は波数を変数とする関数である。この信号を、波数変数信号／位置変数信号変換処理部３６で、例えば図８（ｂ）に示すような、位置を変数とした関数に変換する。この変換処理において、（２）、（３）式のような単純なフーリエ変換ではなくて、屈折率の波数依存性を考慮し、その影響を補正する変換を行う。このとき、波長依存性データ記憶部３４内の光学系及び試料の光学的光路長の波長依存性のデータを用いる。
【００４６】
装置４６がリフレクトメーターの場合は、波数変数信号／位置変数信号変換処理されたデータは、リフレクトメーター・データとしてデータ記憶部３１に記憶され、出力／表示手段３７で表示し出力される。装置４６がOCTの場合は、試料の異なる位置での多数のリフレクトメトリー信号に断層画像処理部３３で断層画像処理を行い、OCT画像を構成して、OCTデータとして記憶装置３９のＯＣＴデータ記憶部３２に記憶させ、出力／表示手段３７で表示し出力する。
【００４７】
光学的光路長の波長依存性をパラメータを用いた式で表現し、測定結果を見てその信号の分解能が最高に近づくように、パラメータを光学的光路長の波長依存性のパラメータ設定部３５によって設定し、分散の効果を逐次的に決定することもできる。
【実施例】
【００４８】
先ず、光学的光路差について説明する。図１に示したOFDR-OCTの光学的光路差ΔLの構成を、分かりやすく図２に示す。この図では、光路の向きを変えるためのガルバノミラー７は省略してある。図２でeaの部分はカプラ１７から出た試料光路１２の光ファイバー、ecの部分はカプラ１７から出た参照光路４の光ファイバーを示す。参照光路４の光ファイバーにおいて、試料光路１２の光ファイバーの長さにあたるea’の部分は、試料光路４の光ファイバーeaと相殺して、光学的光路差ΔLには寄与しない。
【００４９】
試料光路１２において、空気中を進む物理的な光路長ｌsfと試料の中やレンズなどの光学部品の中を進む物理的な光路長の和ｌsmが試料光路１２の物理的な光路長である。図２では、簡単のため、レンズ中の光路長がｌsfに含まれるかのように描いてあるが、これらは、ｌsmに含まれるものとする。光学的な光路長は、物理的な光路長に、光の屈折率を乗じて求められる。空気の屈折率は１とみなせる。物質の屈折率はそれぞれの物質で異なるが、ここでは、実効的な屈折率ｎsmを考え、ｎsmｌsmが物理的な光路長ｌsmに対応した光学的な光路長とする。以上より、光路差を求めるための、試料光路１２の光学的光路長は、
【数９】
2005283155.tif 000011
【００５０】
である。
【００５１】
参照光路４においては、空気中を進む物理的な光路長をｌrf、a’からcまでの光ファイバーとレンズなどの光学部品の物理的光路長をｌrm、実効的な屈折率をｎrmとすると、参照光路４の光学的光路長は、
【数１０】
2005283155.tif 000012
【００５２】
である。
【００５３】
参照光路４と試料光路１２との光学的光路差ΔLは、空気の屈折率を１として、（１）式と（１２）式の差をとって、
【数１１】
2005283155.tif 000013
【００５４】
となる。
【００５５】
以上は、参照光路４と試料光路１２に分散がない場合を考えたが、次に、参照光路４と試料光路１２に分散がある場合について次に説明する。
【００５６】
《光学部品や試料に分散がある場合》
光学部品や試料に分散がある場合、屈折率ｎrmやｎsmが波長に依存し、したがって波数に依存する。波数k=k0におけるこれらの屈折率をｎr0およびｎs0とする。ｉ番目の波数ki=k0＋iΔkにおけるこれらの屈折率を波数の関数としてｎrm(k0＋iΔk)，ｎsm(k0＋iΔk)とし、
【数１２】
2005283155.tif 000014
【００５７】
と表す。これらを用いると、（１３）式の光学的光路長の差ΔLは、波数に依存することになり、ｉ番目の波数では、次式となる。
【数１３】
2005283155.tif 000015
【００５８】
右辺第３項が、波数に依存する項である  ここで、
【数１４】
2005283155.tif 000016
【００５９】
とおくと干渉信号は、（１０）式で与えられる次式となる。
【数１５】
2005283155.tif 000017
【００６０】
になる。KiΔLiの項のΔLiは、ｉに依存して値が変わるので、干渉信号はkiの等間隔の変化に対して一定の周期で変化する正弦波でなく、図３の実線のように変化する。このように、正弦波の周期が変化するとき、信号はチャープしていると呼ばれる。図３の実線の場合、kiの増加に対し、正弦波の周期が長くなり、ΔLiが負の場合を示している。信号がチャープしているとき、（２）式、（３）式と同様な処理を行うと、計算は次式となる。
【数１６】
2005283155.tif 000018
【００６１】
これらの式を用いて、（６）式の計算を行うと、結果は（７）式のようにはならず、（７）式の右辺の第１項と第２項の幅が広がった関数になる。分散がないときは、光学的光路差はΔL0であるが、分散があると、ｉが１からＮまで変わる場合、光学的光路差はΔL0＋ΔL1からΔL0＋ΔLNまで変化するので、これらのピークのFWHMで定義される分解能は、分散があるときは分散がないときに比べておよそΔLNの程度だけ悪くなる。
【数１７】
2005283155.tif 000019
【００６２】
を乗じてｉについて和をとると、（１９）と（２０）式に対応する式はそれぞれ次式のようになる。
【数１８】
2005283155.tif 000020
【００６３】
の式の記号±は、どちらかを使うことを意味する。以下、＋の場合について説明する。三角関数の公式を用い、これらの式は次のように変形できる。
【数１９】
2005283155.tif 000021
【００６４】
これらの式の右辺第１項は、分散による光路差の変化の影響ΔLiが取り除かれた式になっており、分散がない場合の（４）、（５）式の右辺第１項と同じように計算できる。
【数２０】
2005283155.tif 000022
【００６５】
を計算すると、次式が得られる。
【数２１】
2005283155.tif 000023
【００６６】
となる。ここで、
【数２２】
2005283155.tif 000024
【００６７】
式（２７）の第一項は、z=ΔL0の位置にあるピークで、分散の影響が取り除かれて信号が測定されており、上記の計算手法で、分散の影響が補正されていることになる。このピークを測定に用いる。
【００６８】
（２７）式右辺第２項のFc2(z)2＋Fs2(z)2の項は、z=-ΔL0にピークがあり、幅がおよそ2ΔLNの程度広がった関数になっている。こちらのピークは、測定には使用しない。
【００６９】
（２７）の右辺第３項は、小さい値で、測定しようとするピークに対してはノイズになる。
【００７０】
以上のような計算処理の結果を、図９に示す。この図は、図７のz=0の両側の2つのピークに対応している。図９に点線で示す曲線は、分散のために広がった2つのピークである。これに、分散を考慮した数値処理を行うと、z=ΔL0のピークは分散が補正されて、（２７）式の右辺第１項で表される鋭いピークになる。これに対し、（２７）式右辺第２項で表されるz=-ΔL0のピークは、逆に幅が広がる。補正後は、z=ΔL0のピークを信号として用いる。
【００７１】
なお、上述した例では、式（２３）、（２４）の符号±の、＋の場合について説明したが、−の場合について計算すると、逆に、z=-ΔL0におけるピークは分散が取り除かれて測定され、z=ΔL0ピーク分散の影響を受けて幅が広がって観測されることになる。
【００７２】
三角関数には、位相が２πの整数倍ずれている式とか、位相がπ／２ずれているがsin関数をcos関数に、cos関数をsin関数にした式、多項式で展開した式、特殊関数で表した式など、表現は異なるが、実質的には同等な数学的表現が多数ある。この発明は、上記の数学的表現のみに限定されるものではなく、実質的には同じ数学的取り扱いを全て含むものである。
【００７３】
また、以上では、三角関数を用いて、フーリエ変換を変形した、（２３）、（２４）の変換式を例として説明したが、フーリエ変換は、より一般的なウェーブレット変換に含まれるものであり、本発明は、ウェーブレット変換を用いた場合も含む。本発明は、計算ソフトを用いて分散の補正をする装置を全て含むものである。
【００７４】
以上が、本発明の計算の原理である。
【００７５】
〔分散による光学的光路長の変化〕
図６の信号処理部２９における、光学系及び試料の光学的光路長の波長依存性データを求める方法について説明する。図１０に光学的光路長を求める手順を示す。
【００７６】
先ず、光学材料の特性値からの計算処理４０について説明する。光学材料物質の屈折率ｎ(λ)の波長λ依存性は、近似的に次式のセルマイヤーの公式で表される。
【数２３】
2005283155.tif 000025
【００７７】
ここで、B1，B2，B3，C1，C2，C3はセルマイヤーの定数とも呼ばれるパラメータで、各光学材料メーカーのカタログに記載されていることが多い。これらのパラメータの厳密な値は、材料の製造工程にも依存するが、平均的な値は、一般的な光学材料の文献値を用いることも出来る。波数kにおける屈折率の値は、λ＝２π／ｋの関係から上式を使い求められる。したがって、光学的光路長の差の波長依存性ΔLiへの寄与も、それぞれの光学部品について、物理的光路長にその波長での屈折率を（３０）式で計算して求められる。
【００７８】
光学系の実測による決定処理４１について説明する。光学材料の屈折率の特性表がないとか、あっても厳密でないとか、光学系が複雑で、各光学部品の計算値を加えただけでは不正確などの場合、実測によって、光学的光路長波長の依存性を決定できる。図１において、試料8として全反射ミラーを用いて干渉信号を測定すると、分散がなければ図３の点線のようになるが、分散があると実線のようになる。この実線は、波数の増加とともに周期が長くなる場合が描いてあるが、分散の符号が逆であれば、波数の増加とともに周期が短くなる。光学的光路差からΔLiを引けば、実線が点線になるようにΔLiを決定すればよい。
【００７９】
図３のように、単一のミラーの測定値から決定したΔLiの値であっても、試料を除く光学系の光路長差は試料のどこの場所でも同じなので、層状で複数のピークが観測される場合でも、連続的に後方散乱強度が観測される場合でも、同じ値が適用できる。
【００８０】
光学系の実測による決定処理４１のより具体的な一例を示す。屈折率の波長依存性が（３０）式の場合、波数の変化量が大きくなければ、波数ki=k0＋iΔk屈折における屈折率の波数依存性を次のように２次式の展開式で近似できる。
【数２４】
2005283155.tif 000026
【００８１】
分散がある場合の干渉信号（１０）は、波数変化iΔkの関数として、パラメータD1，D2を用いて、次式のように表すことができる。
【数２５】
2005283155.tif 000027
【００８２】
この式において、ki=k0＋iΔkは正確に分かっている。cos の[ ]内の引数は、2πの整数倍だけ不定であるが、変化量としての[ ]内の第２項と第３項のパラメータD1，D2は、（32）式を、図３の実線で示す曲線に、非線形最小二乗法による最適化によって求めることができる。余弦関数（32）の振幅が変わる場合は、振幅も最適化のパラメータとして選んでも良い。このように求めたパラメータD1，D2より、光学的光路長の波長依存性は、次式で求められる。
【数２６】
2005283155.tif 000028
【００８３】
この値は、ミラーのある一つの位置で決定すれば、すべての試料に位置に対して、同じ値が使用できる。
【００８４】
光学系のパラメータ設定処理４２について説明する。（33）式のパラメータD1，D2は、実際の測定値に対して、図６の光学的光路長の波長依存性のパラメータ設定部３５によって、少しずつ値を変化させて波数変数信号／位置変数信号変換処理部３６で繰り返して変換処理を行い、出力／表示手段３７から出力される結果が最適になるように決めることもできる。この処理プロセスが、図６において、出力／表示手段３７と光学的光路長の波長依存性のパラメータ設定部３５で設定されたパラメータを結んだ点線で示してある。
【００８５】
波長依存性データ記憶部３４に記憶された試料の光学的光路長の波長依存性データを生成する方法について説明する。
【００８６】
試料の屈折率の波長依存性が分かっている場合は、試料の特性値から計算処理４３によって、試料の光学的光路長の波長依存性のデータを求めることができる。試料が均一でなく、いくつかの光学的に異なる領域に分かれている場合でも、各縦方向と横方向に対して、異なる設定ができる。
【００８７】
試料の実測による決定処理４４について説明する。試料全体の分散を求めるには、試料の後にミラーを置き、試料を通過し反射してもどった光の干渉信号を図３のように測定する。このとき、光学系による光学的光路差は求めておく必要がある。試料の分散によって実線のような干渉信号が得られたとき、これが点線のような周期的な関数になるような最適化によって、試料の光学的光路長の波長依存性を決定できる。
【００８８】
試料が複数の光学的に異なる領域に分かれている場合は、それぞれの領域のみ通過するようにミラーを試料内に挿入し、それぞれの部位の分散を決定できる。
【００８９】
試料のパラメータ設定部４５のパラメータ設定処理について説明する。試料について、分散特性の実測が困難であったり、試料の特性値から計算処理４３や試料の実測による測定処理４４で設定した試料の光学的光路長の波長依存性データが不十分であるときは、試料の分散の効果を、例えば（32）式や（33）式のD1，D2のようにパラメータで表し、解析結果を図６の出力／表示手段３７で出力させ、最適化するように光学的光路長の波長依存性のパラメータ設定部３５によりパラメータを設定し、これを繰り返し行って最適化させる。このように、計算ソフトウェアで試料の分散を設定することの利点は、（21）、（22）式で、試料の測定位置ｚと、光路長の補正ΔLiを組で設定できるので、試料の縦（奥行き）方向と横方向のそれぞれの測定位置でその位置での分散の補正が行えることである。
【００９０】
〔実施の形態１〕
図１はOFDR-OCT装置に対して本発明を実施する、実施の形態１を示す。実施の形態１では、光源に波長可変レーザ１８を用い、干渉計には光ファイバー・カプラ１７を用いている。波長可変レーザ光源１８は、連続的または離散的に出力光の波長を走査する。波長可変レーザからカプラ１７に入った光は、分波され、試料光路１２と参照光路４に進む。参照光路４を進む光は参照ミラー３で反射され、カプラ１７に戻る。試料光路を進んだ波は、試料で後方反射（散乱の場合も含む）され、カプラ１７に戻る。参照光路４及び試料光路１２よりカプラ１７に戻った光は合波され、互いに干渉する。干渉光は検出器１３で検出され、アンプ１４で増幅され、波長可変レーザ１８の各出力波長に対して、A/Dコンバータ１５によってデジタル信号に変換され、信号処理部２９に記憶され、処理される。処理された信号は出力／表示手段で表示・出力される。
【００９１】
この実施例では、光源に波長可変レーザを用いているが、波長可変光源としては、この他に、スペクトル幅の広い光源から、光学的フィルタによってスペクトル幅の狭い、従って、干渉長の長い光を得、光学的フィルタの選択波長を変化させる装置など、レーザ以外の様々な形態が可能である。
【００９２】
光軸に対して横方向の試料の測定部位を変える場合は、光走査光学系を用いる。図１の場合は。ガルバノミラー７とレンズ９，１０を用いているが、音響光学効果など他の方法を用いてもよい。また、参照ミラー３は様々な種類のミラーが可能である。ファイバーの端面を直接ミラーにしてもよいし、ファイバーに周期的な屈折率の変化を導入してブラッグ反射をさせてもよいし、ファイバーから空気中に光を出し、コリメータ用のレンズ１１を用いてミラーに当ててもよい。
【００９３】
信号処理部２９において、図６に示す信号処理を行えば、本発明が実施できる。この実施例は、請求項１の発明の実施例である。
【００９４】
〔実施の形態２〕
図４は、FD-OCTに対して本発明を実施する、実施の形態２を図示したものである。光源として出射光のスペクトル幅が広く、干渉性の低い低コヒーレンス光源１を用いる。光源を出射した光は、ハーフミラー２で参照光路４と試料光路１２に分波される。参照光路に分波された光は、レンズ１１と参照ミラー３によって反射され、参照光路４を戻って、再びハーフミラーに入射する。ハーフミラー２によって分波され試料光路１２に進んだ光は、試料８に入射し、試料内の各部位で後方反射され、試料光路に戻る。試料の断層像を得るために、試料光路１２の光は、試料上を横方向に、レンズ９、レンズ１０及びガルバノミラー７を用いて走査される。
【００９５】
参照光路４及び試料光路１２から戻ってハーフミラー２に入射した光は、ハーフミラー２で合波される。合波された光の各波数の成分は、（１）式で表される干渉信号となる。ハーフミラー２で合波されて生成された干渉信号は、全ての波数（波長の逆数）の成分が全て重ね合わされている。レンズ２３とレンズ２１、及び回折格子２２を用いると、重ね合わされていた各波数の成分は分離される。波数毎に分離した成分を、１次元アレイ検出器１９上に、レンズ２０を用いて結像投影する。１次元アレイ検出器を構成するそれぞれ個々の検出素子はチャンネルと呼ばれる。各チャンネルは単一の検出素子である必要はなく、複数の素子の集まりをまとめたものでもかまわない。各チャンネル上に波数毎に結合された干渉信号は、アンプ１４、Ａ／Ｄコンバータ１５を通して、信号処理部２９に記憶され計算処理される。計算処理された結果は、出力／表示手段３７で出力・表示される。
【００９６】
信号処理部２９において、図６に示す信号処理を行えば、本発明が実施できる。この実施例は、請求項２の発明の実施例である。
【００９７】
干渉信号を、波数毎に分解するために、図４には、反射型の回折格子が示されているが、スペクトル分解する光学部品としては、他に透過型の回折格子、プリズムなどを用いることができるが、本発明はそれらに限定されるものではない。
【００９８】
１次元アレイ検出器１９としては、ＣＣＤ（電荷結合デバイス）やダイオード・アレイを用いることができるが、本発明はそれらに限定されるものでない。
【００９９】
〔実施の形態３〕
図５は、OFDR-OCT装置において本発明を実施する実施の形態3を示したものである。実施の形態１と実施の形態２は、光学干渉計としてマイケルソン干渉計を用いていたが、図５に示す実施の形態３では、マッハツェンダー干渉計を用いる。波長可変レーザ１８を出射した光は、第１カプラ２５で試料光路１２と参照光路４に分波される。参照光路４は、 光ファイバーを通ってそのまま第２カプラ２６に入る。試料光路１２はサーキュレータ２４を用いて、試料８に入り、試料から後方散乱された光は、サーキュレータ２４により、光ファイバーを介して第２カプラ２６に入る。
【０１００】
試料光路１２におけるサーキュレータ２４の役割は、ポートＩに入った光をポートIIに出し、ポートIIに入った光をポートIIIに出す。レンズ９、レンズ１０およびガルバノミラー７の役割は、試料８に照射する光のビームを、横方向に走査し、断層像の測定を可能にする。
【０１０１】
第２カプラ２６に入った参照光路４の光と試料光路１２の光は互いに干渉する。図５に示すように、第２カプラ２６の2つの出口の出力を差動アンプ２７に入力すると、干渉信号の式（１）の、右辺第１項と第２項の光の強度の直流部分は引き算で相殺され、第３項の干渉信号のみが、差動アンプ２７から出力される。出力された干渉信号は、A/Dコンバータ１５を通して信号処理部２９に記憶され、信号処理がなされる。信号処理の結果は、出力／表示手段３７で出力・表示される。
【０１０２】
信号処理部２９に、図６に示す信号処理を行えば、本発明が実施できる。この実施例は、請求項１の発明の実施例である。
【０１０３】
〔他の実施例〕
以上、３つの実施例を述べたが、光干渉計として図１、図２、図４、図１１、図１２のようにマイケルソン干渉計を用いたり、図５のようにマッハツェンダー干渉を用いることが出来、光学系として図４、図１１のように、バルクな光学部品を用いたバルク光学系を用いたり、図１、図２、図５、図１２のように、光ファイバー光学系を用いることが出来、それらの様々な組み合わせが可能である。3つの実施例は、それらの様々な組み合わせの中から一部を例示したものであり、本発明はそれらを含むが、それらに限定されるものでない。
【０１０４】
実施例は、各光路は、１つの光線の場合を示したが、試料の像を２次元検出器で受光し、２次元検出器の各ピクセルの信号に対して、図６に示す信号処理を行っても、本発明を実施できる。
【０１０５】
本発明は、OFDR-OCT装置やFD−OCT装置のような周波数領域OCT装置において実施できる。それぞれの方法と、２種類の干渉計、バルクかファイバーかの光学系の違い、検出器が点か１次元か２次元かの違いの組み合わせは多様であるが、それらの多様な装置信号に本発明の信号処理が適用でき、本発明が実施できる。
【図面の簡単な説明】
【０１０６】
【図１】ファイバー光学系、マイケルソン干渉計を用い、分散補正を計算処理で実行するOFDR−OCT装置の一例を示す構成図。
【図２】図１に示したOFDR−OCT装置における光学距離の差の考え方を示す説明図。
【図３】周波数領域リフレクトメーター又はOCT装置の干渉信号の波数依存性を示す説明図。
【図４】バルクな光学部品、マイケルソン干渉計を用い、分散補正を計算処理で実行するFD−OCT装置の一例を示す構成図。
【図５】ファイバー光学系、マッハツェンダー干渉計を用い、分散補正を計算処理で実行するOFDR−OCT装置の一例を示す構成図。
【図６】図１、図４、図５で示した装置で使用される信号処理部の機能を示す構成図。
【図７】OFDR法のリフレクトメーターの波数変数の干渉信号を変換処理した、位置変数の信号を示す説明図。
【図８】リフレクトメトリーの信号例を示す説明図。
【図９】分散の補正がなされた信号例を示す説明図。
【図１０】光学的光路長の波数依存性のデータの決定法を示す説明図。
【図１１】バルクな光学部品、マイケルソン干渉計を用い、分散の補償を溶融石英と光学ガラスで実行するOCDR−OCT装置の従来例を示す構成図。
【図１２】ファイバー光学系、マイケルソン干渉計を用い、分散の補償を溶融石英と光学ガラスで実行するOCDR−OCT装置の従来例を示す説明図。
【符号の説明】
【０１０７】
１  低コヒーレンス光源
２  ハーフミラー
３  参照ミラー
４  参照光路
５  溶融石英
６  光学ガラス
７  ガルバノミラー
８  試料
９、１０、１１、２０、２１、２３  レンズ
１２  試料光路
１３  検出器
１４  アンプ
１５  A／D コンバータ
１６  コンピュータ
１７  カプラ
１８  波長可変レーザ
１９  １次元アレイ検出器
２２  回折格子
２４  サーキュレータ
２５  第１カプラ
２６  第２カプラ
２７  差動アンプ
２８  出力／表示手段
２９  信号処理部
３０  元データ記憶部
３１  リフレクトメーター・データ記憶部
３２  OCTデータ記憶部
３３  断層像構成処理部
３４  波長依存性データ記憶部
３５  パラメータ設定部
３６  波数変数信号／位置変数信号変換処理部
３７  出力／表示手段
３８  光学部品の光学的光路長の波長依存性決定処理
３９  試料の光学的光路長の波長依存性決定処理
４０  光学材料の特性値から計算処理
４１  光学系の実測による決定処理
４２  光学系のパラメータ設定部
４３  試料の特性値から計算処理
４４  試料の実測による決定処理
４５  試料のパラメータ設定処理
４６  周波数領域方式リフレクトメーターまたはＯＣＴ装置
【図１】
2005283155.tif 000029
【図２】
2005283155.tif 000030
【図３】
2005283155.tif 000031
【図４】
2005283155.tif 000032
【図５】
2005283155.tif 000033
【図６】
2005283155.tif 000034
【図７】
2005283155.tif 000035
【図８】
2005283155.tif 000036
【図９】
2005283155.tif 000037
【図１０】
2005283155.tif 000038
【図１１】
2005283155.tif 000039
【図１２】
2005283155.tif 000040