(57)【特許請求の範囲】
【請求項１】
複数の測定光を被検査物に照射することで得られる複数の戻り光と参照光とを合成し、当該合成により得られた複数の合成光の信号から、前記被検査物の断層画像を取得する断層撮像装置であって、
光学系を介して前記複数の合成光を検出して、当該複数の合成光の信号を取得するセンサと、
前記複数の合成光の信号間の相違であって、前記光学系と前記センサの構成による前記被検査物の深さ方向の相違を補正した信号を生成する補正手段と、
前記補正された前記複数の合成光の信号から、前記被検査物の断層画像を生成する生成手段とを有することを特徴とする断層撮像装置。
【請求項２】
前記生成手段は、前記補正された前記複数の合成光の信号から、前記被検査物の複数の断層画像を生成し、
前記複数の断層画像の位置合わせを行い、当該位置合わせされた前記複数の断層画像を合成する合成手段を有することを特徴とする請求項１に記載の断層撮像装置。
【請求項３】
前記相違は、前記センサ上における前記複数の合成光それぞれの分解能および減衰特性の少なくともいずれか一つの相違であることを特徴とする請求項１に記載の断層撮像装置。
【請求項４】
前記相違は、前記断層撮像装置内の分散補償用の光学系の構成および前記複数の合成光を回折させる回折格子の構成の少なくとも一つに基づく相違であることを特徴とする請求項１に記載の断層撮像装置。
【請求項５】
前記複数の合成光の信号それぞれは、前記複数の測定光の照射方向に平行な前記被検査物の断面を示すことを特徴とすることを特徴とする請求項１に記載の断層撮像装置。
【請求項６】
複数の測定光を被検査物に照射することで得られる複数の戻り光と参照光とを合成し、当該合成により得られた複数の合成光の信号から、前記被検査物の断層画像を取得する断層撮像装置の制御装置であって、
光学系を介してセンサにより検出された前記複数の合成光の信号を取得する取得手段と、
前記複数の合成光の信号間の相違であって、前記光学系と前記センサの構成による前記被検査物の深さ方向の相違を補正した信号を生成する補正手段とを有することを特徴とすることを特徴とする断層撮像装置の制御装置。
【請求項７】
コンピュータを、
複数の測定光を被検査物に照射することで得られる複数の戻り光と参照光とを合成し、当該合成により得られた複数の合成光の信号から、前記被検査物の断層画像を取得する断層撮像装置の制御装置であって、
光学系を介してセンサにより検出された前記複数の合成光の信号を取得する取得手段と、
前記複数の合成光の信号間の相違であって、前記光学系と前記センサの構成による前記被検査物の深さ方向の相違を補正した信号を生成する補正手段とを有することを特徴とすることを特徴とする断層撮像装置の制御装置として機能させるコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、断層撮像装置及び断層撮像の制御装置に関する。
【背景技術】
【０００２】
現在、光学機器を用いた眼科用機器として、様々なものが使用されている。例えば、眼を観察する光学機器として、前眼部撮影機、眼底カメラ、共焦点レーザー走査検眼鏡等がある。これらの中でも、低コヒーレンス光を利用した光コヒーレンストモグラフィ（ＯＣＴ：Ｏｐｔｉｃａｌ  Ｃｏｈｅｒｅｎｃｅ  Ｔｏｍｏｇｒａｐｈｙ）による光断層撮像装置は、被検査物の断層像を高解像度に得ることができる装置である。よって、眼科用機器として網膜の専門外来では必要不可欠な装置になりつつある。以下、本装置をＯＣＴ装置として説明する。
【０００３】
上記ＯＣＴ装置は、低コヒーレント光を、参照光と測定光とに分け、測定光を被検査物に照射し、その被検査物からの戻り光と参照光を干渉させることによって被検査物の断層を測定することができる。つまり、測定光を被検査物上でスキャンすることで２次元や３次元の断層像を得ることができる。ただし、被検査物が眼のような生体である場合、眼の動きによる画像の歪みが問題となる。そこで、高速で高感度に測定することが求められている。
【０００４】
高速かつ高感度に測定するための方法の一つとして、被検査物の複数点を同時に測定する方法が特許文献１に公開されている。それによると、一つの光源からの光をスリットで分割することによって複数の光源を作り出す。そしてそれらの光を、ビームスプリッタにより、複数の測定光と参照光にそれぞれ分ける。測定光は被検査物に照射され、被検査物からの戻り光と参照光をビームスプリッタで合成する。そして、複数の合成光はグレーティングに入射されて、２次元センサで同時に検出される。このように特許文献１に開示された技術では、複数の測定光による同時測定で高速化を可能としている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特表２００８−５０８０６８
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、複数点の測定を同時に行った断層画像を結合して一枚の画像を作る場合、光学系の構成によって、接続部分が目立つ画像になる。つまり、複数点の測定それぞれに対する光学系が完全に等価なものであれば問題にならないが、そうでない場合には、断層像の深さ方向に依存してコントラストが異なったり、分解能が異なったりすることが起こる。
【０００７】
本発明は上記課題を解決するためになされたものであり、複数の測定光を用いて被検査物の複数の断層画像を生成する場合に、断層撮像装置内の光学系により生じる断層画像間の相違を目立ちにくくすることを目的とする。
【課題を解決するための手段】
【０００８】
本発明の断層撮像装置は、複数の測定光を被検査物に照射することで得られる複数の戻り光と参照光とを合成し、当該合成により得られた複数の合成光の信号から、前記被検査物の断層画像を取得する断層撮像装置であって、光学系を介して前記複数の合成光を検出して、当該複数の合成光の信号を取得するセンサと、前記複数の合成光の信号間の相違であって、前記光学系と前記センサの構成による前記被検査物の深さ方向の相違を補正した信号を生成する補正手段と、前記補正された前記複数の合成光の信号から、前記被検査物の断層画像を生成する生成手段とを有することを特徴とする。
【発明の効果】
【０００９】
本発明によれば、複数の測定光を用いて被検査物の複数の断層画像を生成する場合に、断層撮像装置内の光学系により生じる断層画像間の相違を目立ちにくくすることが出来る。
【図面の簡単な説明】
【００１０】
【図１】第一の実施形態における光断層撮像装置の構成を示す図である。
【図２】第一の実施形態における分光器の構成を示す図である。
【図３】Ｒｏｌｌ−Ｏｆｆの例を説明する図である。
【図４】第一の実施形態における信号処理工程を説明する図である。
【図５】第一の実施形態における、（ａ）眼底、（ｂ）ＡＡ’断面、（ｃ）ＢＢ’断面を説明する図である。
【図６】第二の実施形態における光断層撮像装置の構成を示す図である。
【図７】分散補償を行った後の深さ分解能を説明する図である。
【発明を実施するための形態】
【００１１】
（第一の実施形態）
以下、第一の実施形態について、図面を用いて詳細に説明する。図１は、本実施形態における光断層画像装置の構成を示す図である。ＯＣＴ装置１００は、図１に示されるように、全体としてマイケルソン干渉計を構成している。本実施形態は、分光器の構成によって発生する特性の相違による画像の接続部分を目立ちにくくする。尚、本実施形態および他の実施形態の各機能の処理は、コンピュータが記録媒体からコンピュータプログラムを読み取り、行ってもよい。
【００１２】
図１を用いて、本実施形態における断層撮像装置および断層撮像装置の制御装置の構成について説明する。
【００１３】
（光学系）
光源１０１から照射した照射光である出射光１０４はシングルモードファイバ１１０に導かれて光カプラ１５６に入射し、光カプラ１５６にて第１の光路と第２の光路と第３の光路との３つの光路を通る出射光１０４−１〜３に分割される。
【００１４】
さらに、この３つの出射光１０４−１〜３のそれぞれは、偏光コントローラ１５３−１を通過し、光カプラ１３１−１〜３にて参照光１０５−１〜３と測定光１０６−１〜３とに分割される。このように分割された３つの測定光１０６−１〜３は、観察対象である被検眼１０７における網膜１２７等のそれぞれの測定個所によって反射あるいは散乱された戻り光１０８−１〜３となって戻される。
【００１５】
そして、光カプラ１３１−１〜３によって、参照光路を経由してきた参照光１０５−１〜３と合波され合成光１４２−１〜３となる。合成光１４２−１〜３は、透過型回折格子１４１によって波長毎に分光され、ラインセンサ１３９に入射される。ラインセンサ１３９はセンサ素子毎に各波長の光強度を電圧に変換し、その信号を用いて、被検眼１０７の断層像が構成される。
【００１６】
ここで、光源１０１の周辺について説明する。光源１０１は代表的な低コヒーレント光源であるＳＬＤ（Ｓｕｐｅｒ  Ｌｕｍｉｎｅｓｃｅｎｔ  Ｄｉｏｄｅ）である。波長は眼を測定することを鑑みると、近赤外光が適する。さらに波長は、得られる断層像の横方向の分解能に影響するため、なるべく短波長であることが望ましく、ここでは中心波長８４０ｎｍ、波長幅５０ｎｍとする。観察対象の測定部位によっては、他の波長を選んでも良い。光源の種類は、ここではＳＬＤを選択したが、低コヒーレント光が出射できればよく、ＡＳＥ（Ａｍｐｌｉｆｉｅｄ  Ｓｐｏｎｔａｎｅｏｕｓ  Ｅｍｉｓｓｉｏｎ）等も用いることができる。
【００１７】
次に、参照光１０５の参照光路について説明する。光カプラ１３１−１〜３によって分割された３つの参照光１０５−１〜３のそれぞれは、偏光コントローラ１５３−２を通過し、レンズ１３５−１にて略平行光となって、出射される。次に、参照光１０５−１〜３は分散補償用ガラス１１５を通過し、レンズ１３５−２にて、ミラー１１４に集光される。次に、参照光１０５−１〜３はミラー１１４にて方向を変え、再び光カプラ１３１−１〜３に向かう。そして、参照光１０５−１〜３は光カプラ１３１−１〜３を通過し、ラインセンサ１３９に導かれる。なお、分散補償用ガラス１１５は被検眼１０７および走査光学系を測定光１０６が往復した時の分散を、参照光１０５に対して補償するものである。ここでは、例えば、日本人の平均的な眼球の直径として代表的な値を想定し、Ｌ＝２３ｍｍとする。さらに、１１７は電動ステージであり、矢印で図示している方向に移動することができ、参照光１０５の光路長を、調整・制御することができる。また、電動ステージ１１７はコンピュータ１２５により制御される。なお、本実施形態においては、ミラー１１４や電動ステージ１１７および分散補償用ガラス１１５を、３つの光路で同じのものを使用したが、それぞれ別々の構成としても良い。
【００１８】
次に、測定光１０６の測定光路について説明する。光カプラ１３１−１〜３によって分割された測定光１０６−１〜３のそれぞれは、偏光コントローラ１５３−４を通過し、レンズ１２０−３にて、略平行光となって出射され、走査光学系を構成するＸＹスキャナ１１９のミラーに入射される。ここでは、簡単のため、ＸＹスキャナ１１９は一つのミラーとして記したが、実際にはＸスキャン用ミラーとＹスキャン用ミラーとの２枚のミラーが近接して配置され、網膜１２７上を光軸に垂直な方向にラスタースキャンするものである。また、測定光１０６−１〜３のそれぞれの中心はＸＹスキャナ１１９のミラーの回転中心とほぼ一致するようにレンズ１２０−１、３等が調整されている。レンズ１２０−１、１２０−２は測定光１０６−１〜３が網膜１２７を走査するための光学系であり、測定光１０６を角膜１２６の付近を支点として、網膜１２７をスキャンする役割がある。測定光１０６−１〜３はそれぞれ網膜上の任意の位置に結像するように構成されている。
【００１９】
また、１１７−２は電動ステージであり、矢印で図示している方向に移動することができ、付随するレンズ１２０−２の位置を、調整・制御することができる。レンズ１２０−２の位置を調整することで、被検眼１０７の網膜１２７の所望の層に測定光１０６−１〜３のそれぞれを集光し、観察することができる。測定光１０６−１〜３は被検眼１０７に入射すると、網膜１２７からの反射や散乱により戻り光１０８−１〜３となり、光カプラ１３１−１〜３を通過し、ラインセンサ１３９に導かれる。なお、電動ステージ１１７−２はコンピュータ１２５により制御される。以上の構成をとることにより、３つの測定光を同時にスキャンすることができる。
【００２０】
次に、検出系の構成について説明する。網膜１２７にて反射や散乱された戻り光１０８−１〜３と参照光１０５−１〜３とは光カプラ１３１−１〜３により合波される。そして、合波された合成光１４２−１〜３は分光器に入射し、スペクトルが得られる。これらのスペクトルに対し、コンピュータ１２５が信号処理を行うことで断層像を得ることができる。
【００２１】
（分光器）
ここで、分光器について具体的に説明する。この構成では、複数の合成光を一個のラインセンサで処理するため、２次元センサに比べて低コストを実現することができる。
【００２２】
図２に、図１に示した分光器の部分を詳細に説明するために、３つの合成光（１４２−１〜３）が分光器に入射した場合の構成を示す。ファイバ端１６０−１〜３がお互いに離れて配置してあり、ファイバ端１６０−１〜３からそれぞれ合成光１４２−１〜３が入射する。この際、合成光がレンズ１３５の主面に垂直に入射するように、即ちテレセントリックになるようにファイバ端１６０−１〜３の向きが予め調整されている。入射した合成光は、レンズ１３５に入射する。レンズ１３５で３つの合成光１４２−１〜３はそれぞれ略平行光になり、３つの合成光１４２−１〜３はともに透過型回折格子１４１に入射する。光量の損失を低減するために、透過型回折格子１４１の位置は光学系の瞳近傍に配置し、透過型回折格子１４１の表面に絞りを設ける必要がある。また透過型回折格子１４１はレンズ１３５の主面に対して傾いて配置されるため、透過型回折格子１４１の表面で光束は楕円形となる。そのため透過型回折格子１４１の表面に設けた絞りは楕円形にする必要がある。透過型回折格子１４１で回折された合成光１４２−１〜３は、それぞれレンズ１４３に入射する。ここで図２における回折した合成光は、中心波長のみの光束を示しており、その他の波長の回折した合成光は簡単のため主光線のみを記載している。レンズ１４３に入射したそれぞれの回折した合成光１４２−１〜３は、ラインセンサ１３９上に結像され、矢印１６１−１〜３で示した位置にそれぞれ分光スペクトルが観察される。表１は、本実施形態で使用する測定光の波長の上限と下限および中心波長の８４０ｎｍについてまとめたものである。入射角の違いにより回折角が異なり、その結果、結像位置が合成光によって異なることがわかる。さらに１画素あたり１２μｍのセンサ素子で検出するときの画素数がそれぞれの合成光によって変わることとなる。つまり、断層撮像装置内の光学系などの構成に依存して、ラインセンサ１３９上における合成光それぞれの分布特性が変化することになる。
【００２３】
【表１】
0005637721.tif 000002
【００２４】
ここで、分光器で得られるスペクトルにより、ラインセンサ上においてＯＣＴの信号がどのように歪むかを簡単なモデルで説明する。なお、分光器は波長に対して等間隔になるように設計しているが、信号処理で波数に対して等間隔に変換するため、ここでは波数に対して等間隔になっているものとする。まず、分光された後のスペクトルを波数ｋを用いてｓ（ｋ）と表わす。分光器のラインセンサが有限の大きさであることを考えて、そのウィンドウ関数をｇ（ｋ）で表わすと、分光器で得られるスペクトルは次のようになる。
【００２５】
【数１】
0005637721.tif 000003
【００２６】
ＯＣＴは波数スペクトルのフーリエ変換によって得られる。そのため、式１がコンボリュージョンであることから、式２のようにフーリエ変換後のそれぞれの関数の掛け算と等価である。
【００２７】
【数２】
0005637721.tif 000004
【００２８】
ここで、ウィンドウ関数ｇ（ｋ）が幅Ｗ高さ１の方形波であるとした場合、そのフーリエ変換は式３のように表わされる。
【００２９】
【数３】
0005637721.tif 000005
【００３０】
すなわち、理想的なＯＣＴ像がＦＦＴ（ｓ）であるのに対して、式３のようなＳｉｎ関数が掛け合わされているため原点から最初の節まで減衰することになる。これは一般的にＲｏｌｌ−ＯＦＦ（減衰特性）などと呼ばれている。また、このＲｏｌｌ−Ｏｆｆは幅Ｗによって変化する。つまり、Ｗは分光器の分解能（波数）に相当するため、分解能がよければ減衰が緩くなり、分解能が悪ければ減衰がきつくなる。
【００３１】
図３にＲｏｌｌ−ＯＦＦの例を示す。横軸が距離で、縦軸が強度（デジタル値／１２ビットセンサ）である。この測定で用いた光学系は、図１の参照光路で、分散補償用ガラス１１５がない光路が２本ある状態で、さらにスキャナなどを排した等価なものである。測定光路のミラーの位置は、コヒーレンスゲートに対し−２０００から２０００μｍの間で離散的に移動した。そして、それぞれの位置においてコヒーレンス関数を測定し、それらのデータをつなぎあわせたものである。ここで、コヒーレンスゲートとは、参照光路と測定光路の光路長が一致する位置である。また、原点付近は光源などの自己相関関数であるため、データを排除している。点線は、コヒーレンス関数のピークを模式的につないだ包絡線である。コヒーレンスゲートから離れるにつれ減衰しており、Ｒｏｌｌ−ＯＦＦが起こっていることがわかる。
【００３２】
（信号処理）
図４に第一の実施形態における信号処理の工程を示す。
【００３３】
工程Ａ１で、測定を開始する。この状態は、ＯＣＴ装置が起動されていて、被検眼が配置されている。さらに測定に必要な調整が操作者によって行われ、測定が開始された状態である。
【００３４】
工程Ａ２では、３本の測定光１０６−１〜３をＸＹスキャナ１１９により走査することにより得られる信号をラインセンサ１３９で検出し、検出したデータをコンピュータ１２５が第一の取得手段として機能して取得する。
【００３５】
ここで、眼底５０１と３本の測定光の走査範囲の模式図を図５（ａ）に示す。眼底５０１には、黄斑５０２、視神経乳頭５０３および血管５０４などがある。３本の測定光は、それぞれに対応して第１の走査範囲５０５、第２の走査範囲５０６、第３の走査範囲５０７を走査する。それぞれの領域は、第１の走査範囲−第２の走査範囲間の重なり５０８、第２の走査範囲−第３の走査範囲間の重なり５０９のように２０％程度重なっている。座標軸は図に示したように設定され、ｘ方向がＦａｓｔ−Ｓｃａｎでｙ方向がＳｌｏｗ−Ｓｃａｎで、ｚ方向は、紙面の裏から表に向かっているとする。ここでは、一本の測定光あたりｘ方向には例えば５１２ライン、ｙ方向には例えば２００ライン走査するものとする。ただし、ｙ方向は、重なり部分を除けば３本の測定光で５１２ラインを走査している。ラインセンサ１３９には３本の測定光による合成光１４２−１〜３が入射し、一次元配列のデータ４０９６画素が取得される。そして、ｘ方向の連続する５１２ライン分を、二次元配列のデータを単位として保存する（４０９６×５１２、１２ビット）。走査が終了すればこのデータが１測定光につき２００個保存されることになる。
【００３６】
工程Ａ３では、コンピュータ１２５は、ラインセンサ１３９から取得したデータから各測定光に対応する断層像の生成を行う。生成される断層像は、測定光の照射方向に平行な断面の断層像である。尚、コンピュータ１２５は、断層像補正のために、ラインセンサ１３９上における合成光それぞれの分布特性に関するデータを取得する第二の取得手段としても機能する。
【００３７】
ここで、光学系の構成による眼底の深さ方向の相違として、断層像の物理的な分解能について説明する。この分解能は、一般的に光源の帯域幅によって決まる。ＳＤ−ＯＣＴにおいては、信号処理に使用する画素の最大と最小がそれぞれ光源の波数の最大と最小に一致していれば、分解能は式４で表わされる。
【００３８】
【数４】
0005637721.tif 000006
【００３９】
従って８１５ｎｍ−８６５ｎｍであれば、分解能は空気中で７μｍとなる。また、この値は、一画素あたりの距離と一致する。例えば一画素あたりの距離が７μｍであれば、図３の１０００μｍの位置は、１４２画素ということになる。ただし、表１のように合成光によって画素数が異なる場合、３本の測定光で画像サイズが異なるという不都合が生じる。そのために、画素数を増やして、同じ画像サイズに合わせる。内、高速フーリエ変換をするためにゼロの画素を追加して（ゼロパッティング）全体として２のｎ乗のデータにすると便利である。一方、これは数値的に帯域幅が増加したことになり、一画素あたりの換算距離は短くなる。例えば、合成光１４２−２では、有効な画素が８７０で、これにゼロの要素を１５４個追加することになる。前後にこの要素を等価に分配すれば、概算で８１０ｎｍ−８６９ｎｍの帯域とみなせ、一画素あたりの換算距離（一画素を対応する物理距離に換算）が６μｍとなる。当然、計算に使用する画素が光源の帯域より狭ければ一画素あたりの換算距離が物理的な分解能より悪くなる。
【００４０】
断層像の生成は、このように１ライン当たりの画素数を合わせた後（ここでは例えば１０２４）に行う。断層像の生成は、一般的なＯＣＴ像の生成工程を用いればよい。例えば、固定ノイズ除去、波長波数変換、フーリエ変換などである。次に、図５（ｂ）は、ＡＡ’断面のＢ−Ｓｃａｎ像の模式図を示す。このＢ−Ｓｃａｎ像は、単一の測定光で得られた画像であるため違和感のないものである。しかし、図５（ｃ）は、ＢＢ’断面のＢ−Ｓｃａｎ像の模式図であり、異なる合成光によるＢ−Ｓｃａｎ画像を合成したものである。一画素あたりの分解能の違いによって断面に不連続が発生することになる。これはＣＣ’断面のＣ−Ｓｃａｎ像においては重大で、境界面で血管などの構造がなくなったり、出現したりすることになる。なお、分解能のほかにＲｏｌｌ−Ｏｆｆの違いによるコントラストの違いも発生する。
【００４１】
本工程（工程Ａ３）では、それぞれの合成光に対応してＤｂ（ｐ，ｑ、ｒ）の断層像が得られる。ｐはｚ方向で、１０２４の要素となるが、フーリエ変換の性質によって対称となるため０−５１１だけを抜きだす。ｑはｘ方向（０−５１１）である。ｒはｙ方向であり０−１９９である。またｂは合成光の番号（１−３）である。
【００４２】
なお、データの拡張方法として、スペクトルデータを先に補間して、１０２４画素のスペクトルにしてからフーリエ変換を行ってもよい。また、１ライン当たりの画素数を表１と同じにし、それぞれの断層像を作ったのち補間をしてもよい。
工程Ａ４では、深さ方向の補正を行う。
【００４３】
まず、ｚ方向のリサンプリングを行う。これは、一画素あたりの換算距離を３本の測定光の間で揃えるためのものである。ここでは、一画素あたりの基準となる距離は２番目の測定光（測定領域が中央部の測定光）の換算距離とする。直線補間の場合、各測定光の一画素あたりの換算距離をＬｂとし、ガウス記号を用いて式５のように表わされる。なおガウス記号を用いて［ｘ］は、ｘを超えない最大の整数である。また、ｑ、ｒに関しては同様であるので、ｚ方向ｐのみを使用している。
【００４４】
【数５】
0005637721.tif 000007
【００４５】
補間した結果、各測定光に対する要素数が異なるが、最小のものに合わせれば良く、さらに少なくしてもよい。特に、被検査物が眼の場合、一画素あたりの換算距離が６μｍであるため、４００画素で２．４ｍｍとなり網膜を十分計測できる。ここでは、ｉを０−３９９とする。
【００４６】
次に、ｚ方向の深さ方向のコントラストを合せる規格化を行う。すなわち、すべての測定光のＲｏｌｌ−ＯＦＦ特性をあらかじめ測定またはシュミレーションして求めておく。Ｒｏｌｌ−ＯＦＦ特性をＲｂ（ｐ）とすれば、コントラストは式６のように表わされる。
【００４７】
【数６】
0005637721.tif 000008
【００４８】
なお、Ｒｏｌｌ−ＯＦＦの特性を２番目の測定光に合わせてもよい。
【００４９】
【数７】
0005637721.tif 000009
【００５０】
工程Ａ５では、各測定光間の位置合わせを行う。すなわち、被検査物が眼のように動くものである時、測定の時間差によって画像の位置がずれる場合がある。つまり、図５の第１の領域から第３の領域は、同時に図の左上からｘ方向に走査する。この時、重なり領域５０８、５０９では、同じ位置のデータを取得したいにもかかわらず測定光によって位置がずれる場合がある。このような場合は、重ねあっている領域にある血管などの特徴点を合わせるようにする。
【００５１】
なお、工程Ａ４で、式７を用いて規格化した場合、深さ方向の位置を移動させることは、Ｒｏｌｌ−ＯＦＦの特性で、コントラストが変化することを意味する。そのため、コントラストの調整は工程Ａ５の後であってもよい。
【００５２】
なお、動かないものを観察する場合には、位置ずれがないようにあらかじめ装置の調整がなされているものとする。
【００５３】
工程Ａ６では、断層画像の生成を行う。上記信号処理工程により３Ｄにおけるボリュームデータを得ることができる。そして、ＡＡ’断面、ＢＢ’断面、ＣＣ’断面の任意の位置での切り出しても、特に接続部分に違和感のない画像を生成することができる。
工程Ａ７で、測定の終了をする。
なお、別の被検者があれば上記工程を繰り返す。
当然、測定データがあれば、この測定データに信号処理を加えるだけで断層像を得ることができる。
【００５４】
以上のようにすることで、分光器の特性を主とする画像間の相違を少なくし、接続部分の目立たない断層像を得ることが可能になる。
【００５５】
（第二の実施形態）
第二の実施形態について説明する。ここでは、主に第一の実施形態との差異のみについて説明する。まず、本実施形態のＯＣＴ装置６００は、図６に示されるように、第一の実施形態と同様のマイケルソン干渉系を構成している。しかし、分散補償用ガラス６０１の厚みが各測定光で異なっている点と分光器が測定光に対して等価な３台を使用している点が第一の実施形態と異なる。
【００５６】
ここで、広画角測定を行う場合の課題について説明する。まず、測定光路において、３本の測定光１０６−１〜３がレンズ１２０−１、１２０−２、１２０−３を通る位置が異なる。これは、レンズの収差による問題が発生することを意味している。この問題に対応するために、参照光１０５−１および１０５−３は参照光１０５−２に対して分散補償用ガラスが短くなるように設定されている。つまり、広画角測定で第一の実施形態のように分散補償用ガラスが均一の場合は、周辺部にいくにつれ深さおよび横分解能の低下などが生じる。その理由は、各測定光は走査によって２次元でガラスの厚みの異なる位置を通るが、分散補償用ガラスは平均的な厚みに設定されているためである。広画角の場合は、周辺部のガラスの厚みの差が特に大きくなる。一方、本実施形態のように分散補償用ガラスの厚みを変えた場合には、境界面における画像の接続部分が目立つという課題がある。なお、検出光路においては等価な分光器を利用しているため、これに起因する接続部分の問題は最小限にとどめることができる。
【００５７】
（分散補償）
ここで、分散の影響についてさらに詳しく説明する。先の説明に用いた図３の包絡線は、測定誤差を加味してもプラス側とマイナス側で厳密には対称になっていない。この理由は、干渉計に用いる光カプラやファイバなど個々の部材の機差によって生じるものである。このようにシンプルな光学系においても若干の機差が影響する。そのため、本実施例２のように参照光路において積極的に分散補償用ガラスの厚みを変え、広画角に対応した測定光路を用いる場合には、減衰曲線が異なるだけでなく、深さ分解能も異なる場合がある。
【００５８】
原因が分散の違いによる場合には、分散補償を信号処理で行うことができる。ただし、大きな分散はガラスを用いなければ補正できないが、小さな分散であれば信号処理で補正ができるという現実的な対応である。ここでは、解析関数であるヒルベルト変換を使う。すなわち、式１の分光スペクトルを実部、式１の分光スペクトルをヒルベルト変換（ＨＴ）したものを虚部として、虚数単位ｉを用いて式８のような解析関数を定義する。
【００５９】
【数８】
0005637721.tif 000010
【００６０】
式８の位相成分に対して、式９のような２次（ａ２）および３次（ａ３）の位相成分に対して補正を行う。
【００６１】
【数９】
0005637721.tif 000011
【００６２】
なお、ｋ０は波数の中心で、φ０は初期位相である。補正を行った後の式８の実部を新しい分光スペクトルと置き換えることによって、信号処理によって位相補償が行える。
【００６３】
ここで、先の実験で用いたシンプルな光学系に対して、１７ｍｍと１８ｍｍの厚みの異なるガラスをそれぞれの参照光路と測定光路に挿入した場合について説明する。得られた信号をプラス側の領域で、分解能が良くなるように分散補償のａ２、ａ３パラメータを決めた場合を図７に示す。原点（コヒーレンスゲート）を境に減衰の様子が大きく異なることがわかる。つまり、包絡線は、原点に対して非対称である。また、マイナス側の深さ分解能がプラス側の深さ分解能に対して悪くなっている。すなわち、分散補償を行った場合には、式４で表わされる分解能とは必ずしも一致しないことになる。次に、分散補償を補正する信号処理について説明する。
【００６４】
（信号処理）
本実施形態の信号処理においては、第一の実施形態に対して、工程Ａ３および工程Ａ４の処理が異なる。それぞれＡ３’、Ａ４’と置き換える。
【００６５】
工程Ａ１で、測定を開始する。
工程Ａ２では、３本の測定光と参照光による合成光をラインセンサ１３９で検出し、検出したデータをコンピュータ１２５が取得する。
【００６６】
工程Ａ３’では、コンピュータ１２５が、ラインセンサ１３９のデータから各測定光に対応する断層像の生成を行う。ここでは、境界における分解能が同じになるように分散補償のパラメータを調整する。すなわち、境界領域５０８および５０９を用いて、分解能が略等しくなるようにパラメータを調整する。なお、本実施形態はハードウェアに起因する誤差を処理するために、パラメータは各エリア５０６、５０７、５０８に対応してあらかじめ用意しておく。当然、Ｂ−Ｓｃａｎ像毎に用意してもよいし、さらに細かくライン毎に用意してもよい。
【００６７】
なお、被検査物の影響による分散を補正する場合には、それぞれの画像を比較しながらパラメータを決定することになる。
【００６８】
工程Ａ４’では、深さ方向の補正を行う。分散補償のパラメータに対応した包絡線をあらかじめ用意しておく。この曲線に合わせて式６や式７のような処理を施す。
工程Ａ５では、各測定光間の位置合わせを行う。
工程Ａ６では、断層画像の生成を行う。
工程Ａ７で、測定の終了をする。
【００６９】
以上の処理により、分散の違いを主とする画像間の相違を少なくし、接続部分の目立たない断層像を得ることが可能になる。
【図１】
0005637721.tif 000012
【図２】
0005637721.tif 000013
【図３】
0005637721.tif 000014
【図４】
0005637721.tif 000015
【図５】
0005637721.tif 000016
【図６】
0005637721.tif 000017
【図７】
0005637721.tif 000018