JP6259370B2

JP6259370B2 - 光干渉断層装置

Info

Publication number: JP6259370B2
Application number: JP2014151157A
Authority: JP
Inventors: 上野　雅浩; 雅浩上野; 坂本　尊; 尊坂本; 雄三佐々木; 豊田　誠治; 誠治豊田; 小林　潤也; 潤也小林
Original assignee: Nippon Telegraph and Telephone Corp
Current assignee: Nippon Telegraph and Telephone Corp
Priority date: 2014-07-24
Filing date: 2014-07-24
Publication date: 2018-01-10
Anticipated expiration: 2034-07-24
Also published as: JP2016024176A

Description

本発明は、光干渉断層装置に関し、詳細には、主に医療で使用されている３次元断層画像を取得するための光干渉断層装置に関する。

時間に対して光の波数が連続的に変化（掃引）する光源を波長掃引光源と呼ぶが、その様な光源において、時間に対して光の波数が線形的に（1次関数的に）掃引する光源を波数リニアな波長掃引光源と呼ぶ。

光干渉断層装置としては、波数リニアではない波長可変光源を用いた場合に得た干渉波形信号を、波数リニアとなる波長可変光源を用いた場合と同等の干渉波形信号に変換（リサンプリング）して３次元断層画像を取得するＳｗｅｐｔＳｏｕｒｃｅＯｐｔｉｃａｌＣｏｈｅｒｅｎｃｅＴｏｍｏｇｒａｐｈｙ（ＳＳ−ＯＣＴ）システムが知られている（非特許文献１）。非特許文献１には図１に示す構成のＳＳ−ＯＣＴシステムが記載されている。図１に示すＳＳ−ＯＣＴシステムは、櫛歯型のスペクトルとなるように波長掃引を行なう光源１０１と、カップラ１０２と、干渉計１０３と、フォトディテクタ１０４と、ピーク時刻取得部１１０６と、リサンプリング部１０７と、フーリエ変換部１０８とを備えている。さらに干渉計１０３は、バランスフォトディテクタ１０３１と、サーキュレータ１０３２と、カップラ１０３３と、ファイバコリメータ１０３８と、ミラー１０３５と、ファイバコリメータ１０３９と、光偏向器１０３７とを有している。図中の各構成は光ファイバ（二重線）または電気ケーブル（一本線）で接続されている。カップラ１０３３からファイバコリメータ１０３８を通ってミラー１０３５までの光路を参照アーム１０３４と呼び、また、カップラ１０３３からファイバコリメータ１０３９と光偏向器１０３７を通って観察対象Ｏの中の反射面Ｂまでの光路をサンプルアーム１０３６呼ぶこととする。参照アーム１０３４にあるファイバコリメータ１０３８からミラー１０３４の間と、サンプルアーム１０３６にあるファイバコリメータ１０３９から光偏向器１０３７を経由した反射面Ｂまでの光路は空間系となっており、空間中（通常は空気中）を光が伝播する。

光源１０１は上述の通り櫛歯型のスペクトルとなるように波長掃引を行う光源であるが、これを実現するため、レーザ共振器内に、所定の自由スペクトル間隔（ｆｒｅｅｓｐｅｃｔｒａｌｒａｎｇｅ：ＦＳＲ）を持つフィルタ（櫛歯型スペクトルの透過率を持つフィルタ：櫛歯フィルタ）を入れている。このようなフィルタをレーザ共振器内に挿入したレーザは、フィルタの透過率に応じて出力光強度が変化する。たとえば、フィルタの透過率が極大／極小となる光の周波数となった場合は、出力光強度は極大／極小となる。このため、レーザ共振器の光の共振周波数を連続的に時間変化させた場合に、共振する光の周波数に応じて出力光強度が変わるが、この様子を分散分光方式（モノクロメータを使用する方式）の光スペクトルアナライザで観測すると、光源１０１のスペクトルは、ＦＳＲの間隔でピークを持つ櫛歯型となる。

ＳＳ−ＯＣＴでは、光源のコヒーレンス長が、ＳＳ−ＯＣＴで得られる断層画像で観測できる奥行（深達長）を決め、コヒーレンス長が長ければ深達長が長くなる。このため、コヒーレンス長の長い光源であることがＳＳ−ＯＣＴにおいて重要である。

上記に示す櫛歯型フィルタを共振器に挿入したレーザは、フィルタの透過率が極大、つまり、レーザ出力が極大となる波長の出力光のコヒーレンス長が長くなることが知られている。したがって、光源光強度の強い光、つまり、コヒーレンス長が長い光で得られた干渉計出力のみを抽出して信号処理すると、深達長の長い観測が可能となる。以下に説明する方法は上記抽出を行う方法であり、光源光１０１の出力光の強度変調のピーク値で干渉計１０３の出力をサンプリングする方法である。

光源１０１からは、周期的に波長が掃引された光が出力される。その光強度は、レーザ共振器内に挿入された櫛歯型フィルタの透過率とレーザ共振器の共振周波数に応じて、時間的に変化している。その光強度の極大値（ピーク値）の間隔は、時間をｔ、光源出力光の波数をｋ（ｔ）、ｈを定数としたときの、ｃｏｓ（２ｈｋ（ｔ））の正のピーク間隔と同じである。ｈの値は観測対象において観測したいもっとも深い位置についての情報によって設計するが、これについては後述する。

光源１０１からの出力光はカップラ１０２で分岐する。カップラ１０２から出た一方の光を干渉計１０３に入れ、サーキュレータ１０３２を経由後、カップラ１０３３にて光を２つに分け、参照アーム１０３４とサンプルアーム１０３６に入れる。参照アーム１０３４ではミラー１０３５にて反射する。サンプルアーム１０３６では観測対象の物体Ｏの中の反射面にて反射される。参照アーム１０３４とサンプルアーム１０３６から返ってきた光は再度カップラ１０３３に入って合波されて干渉光（フリンジ光）となる。この干渉光（フリンジ光）は一方はサーキュレータ１０３２を経由し、もう一方は直接バランスフォトディテクタ１０３１に入れ光電変換し、それによって生じた電気信号であるフリンジ信号をリサンプリング部１０７に入れる。

参照アーム１０３４とサンプルアーム１０３６との光路長差ｚは、参照アーム中のミラーの配置位置である基準面Ａと物体Ｏ内の反射面Ｂとの間の距離によって決まる光路長差である。そのフリンジ信号には、光路長差ｚに比例する周波数を持つ信号成分が含まれている。たとえば、時間をｔとしてｋ（ｔ）という関数で光源１０１の出力光の波数を表すと、フリンジ信号にはｃｏｓ（２ｚｋ（ｔ））の信号が含まれる。つまり、時刻ｔにはｚｋ’（ｔ）／πとなる周波数の信号が含まれる。ここで、ｋ’（ｔ）はｋ（ｔ）を時間微分した関数であり、ｋ’（ｔ）＝ｄｋ／ｄｔである。物体Ｏの中の反射面Ｂは、一般に１つとは限らず、複数ある場合がある。たとえば、ｚ₁、…、ｚ_NのＮ個の反射面がある場合は、フリンジ信号はｃｏｓ（２ｚ₁ｋ（ｔ））、…、ｃｏｓ（２ｚ_Nｋ（ｔ））のＮ個の成分を持つこととなる。

カップラ１０２から出たもう一方の光はフォトディテクタ１０４が光電変換し、ピーク時刻取得部１１０６へ渡す。ピーク時刻取得部１１０６は、入力した電気信号のピークの時刻を測定してリサンプリング部１０７に出力する。

図２は図１に示す構成のＳＳ−ＯＣＴシステムの、フォトディテクタ１０４と干渉計１０３から出力された電気信号の処理手順である。

非特許文献１に記載のピーク時刻取得方法は、ピーク検出対象のデータ（フォトディテクタ１０４から出力された電気信号）をアップサンプリングした上でローパスフィルタをかけるものであり、高周波ノイズを除去してアップサンプリングのサンプリング間隔の時間精度でピーク時刻を取得する（Ｓ−１２０１）。

リサンプリング部１０７は、ピーク時刻取得部１０６から入力したピーク時刻に従って、バランスフォトディテクタ１０３１から入力したフリンジ信号をサンプリングし、フーリエ変換部１０８に出力する。（Ｓ−１２０２）フーリエ変換部１０８は、リサンプリング部１０７から入力した、フリンジ信号のサンプリング信号をフーリエ変換して出力する。（Ｓ−１２０３）
リサンプリング部１０７はフリンジ信号に含まれるｃｏｓ（２ｚｋ（ｔ））の成分がｃｏｓ（２ｚ（Δｋ／ΔＴ）ｔ）となるようにフリンジ信号を時間的に伸縮するものである。ここで、Δｋは光源の掃引波数幅、ΔＴはΔｋだけ掃引する時間である。リサンプリング前のフリンジ信号にはｚｋ’（ｔ）／πで表される周波数を含み、ｋ’（ｔ）がｔに対して線形でなければ（ｔの一次関数でなければ）時間によって周波数が変わる。リサンプリングはフリンジ信号をｚ（Δｋ／ΔＴ）ｔ／πという時間変動しない周波数を持つ信号になるように、フリンジ信号を時間的に伸縮処理するものである。これは、光源１０１の出力光がピークとなる時刻がｃｏｓ（２ｈｋ（ｔ））の正のピークとなる時刻となる性質を使って、そのピーク時刻でフリンジ信号をサンプリングことにより得られる。このとき、リサンプリングによってフリンジ信号は、時間当たりΔｋ／ΔＴとなるような、ｋの時間変動が時間に対して線形となるように変換される。このようにｋが時間に対して線形になることを波数リニアと呼ぶ。

光源１０１の出力光の強度変調の周期を決めるｈは、観測対象において観測したいもっとも深い位置における、参照アーム１０３４とサンプルアーム１０３６の光路長差ｚ_maxによって決める。リサンプリング部１０７は、光源１０１の出力光の強度のピーク時刻でフリンジ信号をサンプリングするが、ｚに関するフリンジ光の強度変調の時間ｔにおける周波数はｚｋ’（ｔ）／πであり、光源出力光の強度変調のそれはｈｋ’（ｔ）／πであるため、ｚが０からｚ_maxとなる全てのフリンジ信号ｃｏｓ（２ｚｋ（ｔ））（ただし、０≦ｚ≦ｚ_max）を取得するために、サンプリング定理に基づいてサンプリング周波数はｚ_maxの周波数の２倍以上、つまり、２ｚ_maxｋ’（ｔ）／π≦ｈｋ’（ｔ）／π、すなわち、２ｚ≦ｈとなるようにｈを決める。光源のｈをどのように実現するかについては、たとえば、非特許文献１では、ファイバレーザ中にファイバファブリーペローフィルタ（ｆｉｂｅｒＦａｂｒｙ−Ｐｅｒｏｔｆｉｌｔｅｒ：ＦＦＰＦ）のＦＳＲがｃ／（２ｈ）となるように設計した例がある。ここでｃは光速である。

フーリエ変換部１０８から出力された信号は、ＳＳ−ＯＣＴの深さ方向のスキャン信号であり、この深さ信号のスキャンのことをＡスキャンという。また、このスキャン信号のことをＡスキャン信号と呼ぶこととする。Ａスキャンとは光偏向器等で偏向する等、物理的な光の操作を意味するのではなく、上記の通り、フリンジ信号をフーリエ変換することにより奥行方向の情報を得ることを言う。

Ａスキャン信号はリサンプリングされたフリンジ信号のフーリエ変換であるため、周波数に対する強度を示す信号であるが、周波数ｆは深さ方向の距離ｚに換算でき、それらの関係はｚ＝πｆ／ｋ’と表される。ここで、ｋ’はリサンプリング後のＡスキャン信号においては、理想的にはｋ’＝Δｋ／ΔＴとなるので、リサンプリング後はｚ＝πｆΔＴ／Δｋとなる。

観測対象の物体Ｏ内の反射面は、Ａスキャン信号では信号強度で表されるが、Ａスキャン信号上では、反射面がある位置に相当する周波数の信号が強くなる。反射面が複数ある場合は、それら反射面に対応する、参照アームとサンプルアームの光路長差に応じた成分の信号がフリンジ信号に含まれる。たとえば、Ｎ個の反射点があり、それらに対応する光路長差がｚ₁、…、ｚ_Nであった場合、フリンジ信号にはｃｏｓ（２ｚ₁ｋ（ｔ））、…、ｃｏｓ（２ｚ_Nｋ（ｔ））のＮ個の成分が含まれる。リサンプリングされたフリンジ信号にも同様にｃｏｓ（２ｚ₁Δｋ／ΔＴｔ）、…、ｃｏｓ（２ｚ_NΔｋ／ΔＴｔ）のＮ個の成分が含まれ、それらの周波数はｆ₁＝ｚ₁Δｋ／（πΔＴ）、…、ｆ_N＝ｚ_NΔｋ／（πΔＴ）となる。このリサンプリングされたフリンジ信号をフーリエ変換すると、観測対象の物体Ｏ内の反射強度に応じて、上記ｆ₁、…、ｆ_Nに相当する周波数の位置の強度が大きい信号が現れる。これらの周波数をｚ＝πｆΔＴ／Δｋで換算すると、奥行距離と反射の強度の対応を取ることができる。

観測対象の物体Ｏ内の反射面が１つのみであり、その反射面がこの時の参照光とサンプル光の光路長差ｚに相当する箇所にある場合は、Ａスキャン信号のｚに相当する箇所、つまり、周波数ｆ＝ｚ（Δｋ／ΔＴ）／πの信号が強くなり、それ以外の周波数の信号は０（デルタ関数）になることが望ましい。このような、１つの反射点を持つ場合のフリンジ信号から得られるＡスキャン信号を点広がり関数（ｐｏｉｎｔｓｐｒｅａｄｆｕｎｃｔｉｏｎ：ＰＳＦ）と呼ぶ。しかし、ΔＴ内での光源光強度の変調や、光学的および電気的なノイズ、リサンプリング等の信号処理の不十分さ等が原因で、実際のＰＳＦはデルタ関数とはならず、ＰＳＦの半値幅は大きくなり、また、ＳＮＲ（ＰＳＦのピーク／ノイズフロアレベル）は低く（悪く）なる。

リサンプリングとは、波数リニアでないフリンジ信号を波数リニアとなるようにする処理であるが、サンプリングする時刻に誤差が生じると（正確なリサンプリングができないと）、正確な波数リニアなフリンジ信号を生成することはできなくなり、ＰＳＦの半値幅は大きくなり、また、ＰＳＦのＳＮＲは悪化する。

Tsung-Han Tsaia, Chao Zhoua, Desmond Adlera, and James G. Fujimotoa, "Frequency Comb Swept Lasers for Optical Coherence Tomography," Proc. of SPIE , Vol. 7554, pp. 7554E-1-10 (2010).

ピーク時刻取得部１１０６において光源光の強度変調信号のピークが現れる時刻を取得する際のピーク時刻の検出精度はサンプリング周波数に依存する。したがってサンプリング周波数を高くすればピーク時刻の検出精度は向上するが、その分ピーク時刻検出処理に多大なリソースを要することになる。一方、低いサンプリング周波数でピーク時刻の検出を行うと、ピーク時刻検出精度が劣化してしまう。

すなわち、ピーク時刻検出精度が低い場合は、その後の処理で正確なリサンプリングが難しくなるため、光源１０１の出力光のコヒーレンス長の長い光で干渉計１０３で干渉させて得たフリンジ信号のみを抽出することが難しくなるので、Ａスキャン信号の深達長が長くならない、つまり、断層画像が深くまで見ることができないという問題がある。また、同じ原因により、フーリエ変換部１０８から出力されるＡスキャン信号のＰＳＦの半値幅が大きくなりやすいという問題と、ＳＮＲが悪くなりやすいという問題がある。

本発明は上記従来の問題に鑑みてなされたものであり、本発明の課題は、簡素な処理で正確なリサンプリングを行うことができる光干渉断層装置を提供することにある。

上記の課題を解決するために、一実施形態に記載された発明は、櫛歯型のスペクトルとなるように波長掃引する光源と、前記光源から出力された光を２分岐する第１のカップラと、前記２分岐された一方の光を観測対象の物体に入射して観測対象の物体の構造に応じた干渉信号であるフリンジ信号を出力する干渉計と、前記２分岐された他方の光の強度信号の位相を算出する位相算出部と、前記位相算出部で算出された位相が２ｎπとなる時刻であるピーク時刻を取得する位相２ｎπ時刻取得部と、前記ピーク時刻におけるフリンジ信号をサンプリングするリサンプリング部とを備えることを特徴する光干渉断層装置である。

従来のＳＳ−ＯＣＴシステムの構成例を示す図である。図１のＳＳ−ＯＣＴシステムのフォトディテクタから出力された信号の処理フローを示す図である。第１の実施形態のＳＳ−ＯＣＴシステムの構成例を示す図である。空間光学系で構成した光源の構成例を示す図である。図３の装置のフォトディテクタから出力された信号の処理フローを示す図である。強度信号と位相の経時変化を示す図である。

以下、本発明の実施の形態について、詳細に説明する。本発明の光干渉断層装置は、櫛歯型のスペクトルとなるように波長掃引する光源と、前記光源から出力された光を２分岐する第１のカップラと、前記２分岐された一方の光を観測対象の物体に入射して観測対象の物体の構造に応じた干渉信号であるフリンジ信号を出力する干渉計と、前記２分岐された他方の光の強度信号の位相を算出する位相算出部と、前記位相算出部で算出された位相が２ｎπとなる時刻であるピーク時刻を取得する位相２ｎπ時刻取得部と、前記ピーク時刻におけるフリンジ信号をサンプリングするリサンプリング部とを備えている。

上記光干渉断層装置によれば、光源出力光のピーク時刻を正確に得られるため、正確なリサンプリングを行うことができる。これにより、光源の出力光強度のピーク時刻に現れるコヒーレンス長の長い光の干渉信号を正確にサンプリングすることが可能となり、断層画像の深達長を長くすること、つまり、観察対象の物体を深くまで観測することが可能となる。

また、上記の通り、正確なリサンプリングが可能となることにより、リサンプリング部から出力されるフリンジ信号は光源が正確に波数リニアな掃引した場合のフリンジ信号と同一となるため、フーリエ変換部から出力されるＡスキャン信号のＰＳＦの半値幅は細くなり、また、ＳＮＲが大きくなるため、鮮明度は高くなる。

（第１の実施形態）
第１の実施形態の光干渉断層装置について説明する。図３は、第１の実施形態のＳＳ−ＯＣＴシステムの構成例を示す図である。本実施形態のＳＳ−ＯＣＴシステムは、図３に示すように、櫛歯型のスペクトルとなるように波長掃引を行なう光源１０１と、カップラ１０２と、干渉計１０３と、フォトディテクタ１０４と、位相算出部１０５と、位相２ｎπ時刻取得部１０６と、リサンプリング部１０７と、フーリエ変換部１０８とを備えている。さらに干渉計１０３は、バランスフォトディテクタ１０３１と、サーキュレータ１０３２と、カップラ１０３３と、ファイバコリメータ１０３８と、ミラー１０３５と、ファイバコリメータ１０３９と、光偏向器１０３７とを有している。カップラ１０３３からファイバコリメータ１０３８を通ってミラー１０３５までの光路が参照アーム１０３４であり、カップラ１０３３からファイバコリメータ１０３９と光偏向器１０３７を通って観察対象Ｏの中の反射面Ｂまでの光路がサンプルアーム１０３６である。図中の各構成は光ファイバ（二重線）または電気ケーブル（一本線）で接続されている。

光源１０１は、時間的に波長（周波数）が推移する（掃引する）レーザであり、その発振スペクトルのピークは周波数に対して等間隔に離散的な値をとる。その発振スペクトルを１掃引分重ね合わせた（１掃引分を時間積分した）スペクトルは、周波数に対して等間隔にピークがある櫛歯型のスペクトルとなる。この光源１０１は、周波数が時間的に推移するので、レーザ発振は時間的に強くなったり弱くなったりを繰り返す。つまりこの光源１０１の出力光は強度変調する。この強度変調はｃｏｓ（２ｈｋ（ｔ））の基本成分とその高調波成分を持つ。ここで、ｃｏｓ（ｘ）はｘの余弦、ｈは定数、ｋ（ｔ）は光の波数であり、時間ｔにより変化する。このｈは光源内のレーザ共振器中にある光学部品のパラメータであるが、測定対象の物体の最大深さに基づいて決定される定数である。

強度変調がｃｏｓ（２ｈｋ（ｔ））となる光源１０１は、その光源１０１で生じる櫛歯型のスペクトルの隣同士の櫛歯の周波数間隔（自由スペクトル間隔：ＦＳＲ）がｃ／（２ｈ）となる光源である。ここでｃは光速である。その光源１０１の構成としては、ファイバレーザ中にファイバファブリーペローフィルタ（ｆｉｂｅｒＦａｂｒｙ−Ｐｅｒｏｔｆｉｌｔｅｒ：ＦＦＰＦ）を挿入したものや、空間光学系で構成したレーザ共振器内にエタロンを挿入したものがある。このときのＦＦＰＦやエタロンはそのＦＳＲがｃ／（２ｈ）となるものを使用することで、ＦＳＲがｃ／（２ｈ）となる光源が実現する。

図４に空間光学系で構成した光源１の構成例を示す。図４に示す構成は空間光学系で構成されたＬｉｔｔｍａｎ―Ｍｅｔｃａｌｆ型共振器であり、半導体光増幅器（ｓｅｍｉｃｏｎｄｕｃｔｏｒｏｐｔｉｃａｌａｍｐｌｉｆｉｌｅｒ：ＳＯＡ）３０４から発した光は、コリメータレンズ３０６でコリメートされた後、エタロン３０７を通して光偏向器３０８で偏向されて回折格子３１０に入射し、回折した光をミラー３１１に送る。ミラー３１１に垂直入射された光のみが、逆の経路をたどってＳＯＡ３０４に戻り、ＳＯＡ反射面３０５で折り返され、最終的にはＳＯＡ反射面３０５とミラー３１１の間でレーザ共振する。ＳＯＡ反射面３０５はハーフミラーとなっており、レーザ発振した光の一部がＳＯＡ反射面を通って結合レンズ３０３、アイソレータ３０２を介して光ファイバ３０１に至り、共振器外部に出力される。なお、エタロンの光路長はｈである物を使う。たとえば、エタロンの屈折率をｎ、エタロンの厚みをｄとすると、ｈ＝ｎｄとなる。レーザ共振器の光の共振周波数は、光偏向器３０８から回折格子３１０への入射角度によって変化するので、光偏向器制御器３０９によって光偏向器３０８から出射する光の角度（偏向角）を時間的に連続的に変えることによって、光源１０１の出力光の周波数を連続的に変える（掃引する）。

カップラ１０２は、光源１０１から入力された光を２分岐し、一方を干渉計１０３に入力し、他方をフォトディテクタ１０４に入力する。

カップラ１０２で分岐された一方の光が入力される干渉計１０３では、入力された光をサーキュレータ１０３２を介してカップラ１０３３に入れ、カップラ１０３３で２分岐した後、参照アーム１０３４とサンプルアーム１０３６に入力する。カップラ１０３３で分岐された一方の光は、参照アーム１０３４の基準位置Ａに設けられたミラー１０３５で反射して再びカップラ１０３３に戻る。カップラ１０３３で分岐されたもう一方の光は、観測対象の物体Ｏ内の反射面Ｂ（反射面Ｂは物体内の一か所とは限らない）で反射されて再びカップラ１０３３に戻る。光源が１掃引するごとに光偏向器１０３７で幅方向ＷにＡスキャンを行う位置を少しずつ移動させる。このように、物理的に光を観測対象にあてる位置（Ａスキャンを行う位置）を連続的に変える操作をＢスキャンという。

カップラ１０３３では参照アーム１０３４およびサンプルアーム１０３６からの２つの戻り光を合波して干渉させ、その干渉光（フリンジ光）を、お互いに位相が反転した状態で２分岐し、一方はサーキュレータ１０３２を介して、他方は直接バランスフォトディテクタ１０３１に入射する。バランスフォトディテクタ１０３１は、それぞれの入射光を電気信号に変換し、それらの差分信号をリサンプリング部１０７に出力する。この出力信号を干渉信号（フリンジ信号）と呼ぶこととする。

干渉計１０３から出力されるフリンジ信号（電気信号）は、参照アームとサンプルアームとの光路長差ｚに応じて変化する信号である。具体的には、フリンジ信号はｃｏｓ（２ｚｋ（ｔ））の成分を持つ。また、光源１０１からの出力光はｃｏｓ（２ｈｋ（ｔ））で強度変調しているので、フリンジ信号はｃｏｓ（２ｚｋ（ｔ））およびｃｏｓ（２ｈｋ（ｔ））で表される２つの基本成分とその高調波成分を持つこととなる。観測対象の物体Ｏ内の反射面の基準面Ａからの最深の位置ｚをｚ_maxとすると、光源１０１の強度変調の周波数を決定し、また、リサンプリングする時間間隔を決定するｈは２ｚ_max以上（ｈ≧２ｚ_max）となるようにする。このようにする理由は次に示すような原理による。ｚに対するフリンジ信号はｃｏｓ（２ｚｋ（ｔ））で変調されており、その時刻ｔにおける周波数ｆｚ（ｔ）はｚｋ’（ｔ）／πとなる。ただし、ｋ’（ｔ）は時刻ｔにおけるｋ（ｔ）の微分した関数でありｋ’（ｔ）＝ｄｋ／ｄｔである。一方、光源１０１の光出力はｃｏｓ（２ｈｋ（ｔ））で変調されており、その時刻ｔにおける周波数ｆｈ（ｔ）はｈｋ’（ｔ）／πとなる。光源１０１の光出力のピークの時刻でフリンジ信号をサンプリングするので、時刻ｔあたりではｈｋ’（ｔ）／πで表されるサンプリング周波数でフリンジ信号をサンプリングすることになる。観測対象の物体Ｏ内の反射点によるフリンジの最高周波数はｚ_maxとなる反射点からの光で生じる干渉光の周波数ｚ_maxｋ’（ｔ）／πであるので、この周波数の信号成分の情報をすべて取得するには、サンプリング定理から、２ｚ_maxｋ’（ｔ）／π≦ｈｋ’（ｔ）／πである必要がある。したがって、ｈは２ｚ_max以上（、ｈ≧２ｚ_max）となるようにする。

カップラ１０２で分岐された他方の光が入力されるフォトディテクタ１０４では、入力された光を光強度に応じた光源強度信号（電気信号）に変換し、位相算出部１０５に送る。

図５は、位相算出部１０５、位相２ｎπ時刻取得部１０６、リサンプリング部１０７によるフォトディテクタ１０４と干渉計１０３とから出力された電気信号に対する処理手順である。

位相算出部１０５は、フォトディテクタ１０４から出力された光源光強度信号を入力し、信号の時刻ごとの位相を算出し、位相２ｎπ時刻取得部１０６に出力する（Ｓ−２０１）。

位相算出の一例を以下に示す。ただし、フォトディテクタ出力信号が光源強度に比例する場合は、位相算出部１０５の入力信号はｃｏｓ（２ｈｋ（ｔ））に比例するが、説明を簡単にするため、ここでは余弦関数の係数を１として、つまり、入力信号ｓ（ｔ）をｃｏｓ（２ｈｋ（ｔ））として位相算出手段とその原理を以下に説明する。

１）入力信号ｓ（ｔ）をフーリエ変換する。その結果Ｓ（ω）は以下のようになる。

ここで、

はフーリエ変換を示し、ｊは虚数単位を表す。上記式から、Ｓ（ω）は時刻ｔにおいて、±ｚｋ’（ｔ）／πの２つの周波数（±２ｚｋ’（ｔ）の２つの角周波数）を持つ信号であることが分かる。

２）Ｓ（ω）のマイナスの周波数成分を削除する。つまり、Ｓ（ω）の

を削除する。マイナスの周波数成分が削除された信号をＳ’（ω）とすると、Ｓ’（ω）は以下のようになる。

３）Ｓ’（ω）を逆フーリエ変換する。その結果をｓ’（ｔ）とするとｓ’（ｔ）は以下のようになる。

４）ｓ’（ｔ）の位相を算出する。具体的な位相の算出方法としては、たとえば、ｓ’（ｔ）の実部と虚部から位相を算出する。算出した位相は、上式中の２ｚｋ（ｔ）と一致する。

位相２ｎπ時刻取得部１０６は、位相算出部１０５から得た位相が２ｎπとなる（ただしｎは整数）となる時刻をピーク時刻としてリサンプリング部１０７に出力する。位相が２ｎπとなる時刻について説明すると、光源強度信号はｃｏｓ（２ｚｋ（ｔ））に比例するため、２ｚｋ（ｔ）＝２ｎπとなる時刻ｔとは光源強度は正のピークとなる時刻のことである（Ｓ−２０２）。

図６は強度信号ｓ（ｔ）と位相２ｚｋ（ｔ）を示したものである。ただし、説明を簡単化するために、位相を２πでラッピングして表記おり、図６では位相は−π〜πの範囲で示されている。したがって、図６では位相２ｎπは０で示されている。図６に示すように、信号の正のピーク位置において位相が０となっていることが分かる。したがって、位相が０、すなわち２ｎπとなる時刻を強度信号ｓ（ｔ）が正のピーク位置となるピーク時刻として検出することができる。

図６では、強度信号ｓ（ｔ）に重畳したノイズのために、位相が０付近で行ったり来たりして、強度信号ｓ（ｔ）のピーク位置の付近で複数回位相が０となる現象がでているが、この場合は、ローパスフィルタ等で強度信号ｓ（ｔ）から高周波ノイズを除去することにより、ピーク位置付近で位相を１回のみ０にすることができる。

リサンプリング部１０７は、位相２ｎπ時刻取得部１０６で取得した時刻で、干渉計１０３から出力されたフリンジ信号をサンプリングする。（Ｓ−２０３）
フーリエ変換部１０８は、リサンプリング部１０７でサンプリングしたフリンジ信号をフーリエ変換する。このフーリエ変換後の信号がＳＳ−ＯＣＴのＡスキャン信号となる。（Ｓ−２０４）

以上のＯＣＴシステムによれば、サンプリング周波数を高くすることなく、また、それによってピーク検出に多くの計算リソースを使用することなく、光源出力光のピーク時刻を正確に得られるため、簡素な処理で正確なリサンプリングを行うことができる。

１０１光源
１０２カップラ
１０３干渉計
１０３１バランスフォトディテクタ
１０３２サーキュレータ
１０３３カップラ
１０３４参照アーム
１０３５ミラー
１０３６サンプルアーム
１０３７光偏向器
１０３８ファイバコリメータ
１０３９ファイバコリメータ
１０４フォトディテクタ
１０５位相算出部
１０６位相２ｎπ時刻取得部
１０７リサンプリング部
１０８フーリエ変換部
１１０６ピーク時刻取得部
３０１光ファイバ
３０２アイソレータ
３０３結合レンズ
３０４ＳＯＡ
３０５ＳＯＡ反射面
３０６コリメータレンズ
３０７エタロン
３０８光偏向器
３０９光偏向器制御器
３１０回折格子
３１１ミラー

Claims

櫛歯型のスペクトルとなるように波長掃引する光源と、
前記光源から出力された光を２分岐する第１のカップラと、
前記２分岐された一方の光を観測対象の物体に入射して観測対象の物体の構造に応じた干渉信号であるフリンジ信号を出力する干渉計と、
前記２分岐された他方の光の強度信号の時刻ごとの位相を算出する位相算出部と、
前記位相算出部で算出された位相が２ｎπとなる時刻であるピーク時刻を取得する位相２ｎπ時刻取得部と、
前記ピーク時刻におけるフリンジ信号をサンプリングするリサンプリング部とを備えることを特徴する光干渉断層装置。
前記干渉計は、入力された光を２分岐し一方の光を参照アームに入力して基準位置Ａに配置されたミラーで反射させ、他方の光をサンプルアームに入力して観測対象の物体の反射面Ｂで反射させた２つの光を合波して干渉した干渉光を出力する第２のカップラと、該干渉光をフリンジ信号に変換して出力するバランスフォトディテクタとを備えることを特徴とする請求項１記載の光干渉断層装置。