(57)【特許請求の範囲】
【請求項１】
被検体上を走査手段によって走査された測定光と、前記測定光に対応する参照光とによるＯＣＴ信号を取得するためのＯＣＴ光学系を備える光コヒーレンストモグラフィ装置であって、
前記測定光を横断方向に走査するＢスキャンを被検体上の同一位置に関して複数回行うことによって、各Ｂスキャン間で時間間隔を空けて取得された複数のＯＣＴ信号のモーションコントラストを取得する取得手段と、
前記モーションコントラストの基礎となる前記複数のＯＣＴ信号間での位相差のプロファイルが連続であるか否かを判定する判定手段と、
前記判定手段によって前記プロファイルが不連続であると判定された場合、前記走査手段の駆動を制御し、各Ｂスキャン間の前記時間間隔を変更する駆動制御手段と、
前記駆動制御手段によって変更された前記時間間隔を空けて取得された前記複数のＯＣＴ信号のモーションコントラストに基づいて前記被検体の血流速度を算出する制御手段と、
を備えることを特徴とする光コヒーレンストモグラフィ装置。
【請求項２】
前記制御手段は、３次元モーションコントラスト画像を生成するために取得した複数のＯＣＴ信号を用いて、ある時点での前記血流速度を求め、前記血流速度の大きさに基づいて前記時間間隔を設定することを特徴とする請求項１の光コヒーレンストモグラフィ装置。
【請求項３】
被検体上を走査手段によって走査された測定光と、前記測定光に対応する参照光とによるＯＣＴ信号を取得するためのＯＣＴ光学系を備える光コヒーレンストモグラフィ装置であって、
前記測定光を横断方向に走査するＢスキャンを被検体上の同一位置に関して複数回行うことによって、各Ｂスキャン間で第１の時間間隔を空けて取得された複数のＯＣＴ信号のモーションコントラストを取得する取得手段と、
前記被検体の血管の３次元構造に基づいて前記血管の径および種類の少なくともいずれかを取得する血管情報取得手段と、を備え、
前記取得手段は、前記血管情報取得手段によって取得された前記血管の径または種類に基づいて、取得された前記血管の径または種類に応じて予め設定された、前記第１の時間間隔とは異なる第２の時間間隔を各Ｂスキャン間で空けて取得された複数のＯＣＴ信号のモーションコントラストを取得することを特徴とする光コヒーレンストモグラフィ装置。
【請求項４】
被検体上を走査手段によって走査された測定光と、前記測定光に対応する参照光とによるＯＣＴ信号を取得するためのＯＣＴ光学系を備える光コヒーレンストモグラフィ装置に用いられる制御プログラムであって、前記光コヒーレンストモグラフィ装置のプロセッサによって実行されることで、
前記測定光を横断方向に走査するＢスキャンを被検体上の同一位置に関して複数回行うことによって、各Ｂスキャン間で時間間隔を空けて取得された複数のＯＣＴ信号のモーションコントラストを取得する取得ステップと、
前記モーションコントラストの基礎となる前記複数のＯＣＴ信号間での位相差のプロファイルが連続であるか否かを判定する判定ステップと、
前記判定ステップにおいて前記プロファイルが不連続であると判定された場合、前記走査手段の駆動を制御し、各Ｂスキャン間の前記時間間隔を変更する駆動制御ステップと、
前記駆動制御ステップにおいて変更された前記時間間隔を空けて取得された前記複数のＯＣＴ信号のモーションコントラストに基づいて前記被検体の血流速度を算出する算出ステップと、
を前記光コヒーレンストモグラフィ装置に実行させることを特徴とする制御プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、被検体のモーションコントラストデータを得る光コヒーレンストモグラフィ装置、及びそれに用いる制御プログラムに関する。
【背景技術】
【０００２】
従来の光コヒーレンストモグラフィ装置（ＯＣＴ装置とも言う）において、網膜の血流を測定する装置が知られている。これは、血管を流れる血球等によって生じるドップラーシフトによって、被検眼に照射された測定光の周波数が変化することを利用したものである。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】国際公開第２０１０／１４３６０１号
【非特許文献】
【０００４】
【非特許文献１】H.C.Hendargo et al. Biomed. Opt. Express, Vol. 4, No.6, p.803/ May 2013
【非特許文献２】Yonghua Zhao et al. OPTICS LETTERS / Vol. 25, No. 2 / January 15, 2000
【非特許文献３】Adrian Mariampillai et al. OPTICS LETTERS / Vol. 33, No. 13 / July 1, 2008
【非特許文献４】Vivek J. Srinivasan et al. OPTICS LETTERS / Vol. 35, No. 1 / January 1, 2010
【非特許文献５】P. Meemon and J. Rolland, "Swept-source based, single-shot, multi-detectable velocity range Doppler optical coherence tomography," Biomed. Opt. Express, vol. 1, no. 3, pp. 3116~3121, 2010.
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ドップラーシフトによる周波数変化は僅かであり、光波の周波数変化はその位相変化として検出される。位相変化の検出範囲は、例えば−π〜πの範囲などに限られることが多い。したがって、検出範囲外の位相変化が生じた場合、実際に生じた位相変化と検出された位相変化とが相違し、位相変化のプロファイルの波形が不連続になることがあった。このように、不連続な波形のプロファイルは、血流の測定に不向きであった。
【０００６】
本開示は、上記問題点を鑑み、被検体の血流速度を好適に取得できる光コヒーレンストモグラフィ装置を提供することを技術課題とする。
【課題を解決するための手段】
【０００７】
上記課題を解決するために、本発明は以下のような構成を備えることを特徴とする。
【０００８】
（１）  被検体上を走査手段によって走査された測定光と、前記測定光に対応する参照光とによるＯＣＴ信号を取得するためのＯＣＴ光学系を備える光コヒーレンストモグラフィ装置であって、前記測定光を横断方向に走査するＢスキャンを被検体上の同一位置に関して複数回行うことによって、各Ｂスキャン間で時間間隔を空けて取得された複数のＯＣＴ信号のモーションコントラストを取得する取得手段と、前記モーションコントラストの基礎となる前記複数のＯＣＴ信号間での位相差のプロファイルが連続であるか否かを判定する判定手段と、前記判定手段によって前記プロファイルが不連続であると判定された場合、前記走査手段の駆動を制御し、各Ｂスキャン間の前記時間間隔を変更する駆動制御手段と、前記駆動制御手段によって変更された前記時間間隔を空けて取得された前記複数のＯＣＴ信号のモーションコントラストに基づいて前記被検体の血流速度を算出する制御手段と、を備えることを特徴とする。
（２）  被検体上を走査手段によって走査された測定光と、前記測定光に対応する参照光とによるＯＣＴ信号を取得するためのＯＣＴ光学系を備える光コヒーレンストモグラフィ装置であって、前記測定光を横断方向に走査するＢスキャンを被検体上の同一位置に関して複数回行うことによって、各Ｂスキャン間で第１の時間間隔を空けて取得された複数のＯＣＴ信号のモーションコントラストを取得する取得手段と、前記被検体の血管の３次元構造に基づいて前記血管の径および種類の少なくともいずれかを取得する血管情報取得手段と、を備え、前記取得手段は、前記血管情報取得手段によって取得された前記血管の径または種類に基づいて、取得された前記血管の径または種類に応じて予め設定された、前記第１の時間間隔とは異なる第２の時間間隔を各Ｂスキャン間で空けて取得された複数のＯＣＴ信号のモーションコントラストを取得することを特徴とする。
（３）  被検体上を走査手段によって走査された測定光と、前記測定光に対応する参照光とによるＯＣＴ信号を取得するためのＯＣＴ光学系を備える光コヒーレンストモグラフィ装置に用いられる制御プログラムであって、前記光コヒーレンストモグラフィ装置のプロセッサによって実行されることで、前記測定光を横断方向に走査するＢスキャンを被検体上の同一位置に関して複数回行うことによって、各Ｂスキャン間で時間間隔を空けて取得された複数のＯＣＴ信号のモーションコントラストを取得する取得ステップと、前記モーションコントラストの基礎となる前記複数のＯＣＴ信号間での位相差のプロファイルが連続であるか否かを判定する判定ステップと、前記判定ステップにおいて前記プロファイルが不連続であると判定された場合、前記走査手段の駆動を制御し、各Ｂスキャン間の前記時間間隔を変更する駆動制御ステップと、前記駆動制御ステップにおいて変更された前記時間間隔を空けて取得された前記複数のＯＣＴ信号のモーションコントラストに基づいて前記被検体の血流速度を算出する算出ステップと、を前記光コヒーレンストモグラフィ装置に実行させることを特徴とする。
【図面の簡単な説明】
【０００９】
【図１】光コヒーレンストモグラフィ装置の構成について説明するブロック図である。
【図２】光学系の概略を示す図である。
【図３】本実施例の処理を説明するためのフローチャートである。
【図４】本実施例の測定について説明するための眼底のイメージ図である。
【図５】断層血管像を示すイメージ図である。
【図６】本実施例の血流測定の制御を説明するフローチャートである。
【図７】モーションコントラスト画像と観察画像の対応付けを説明する図である。
【図８】表示部に表示されたモーションコントラスト画像の一例を示す図である。
【図９】血流速度の求め方を説明する図である。
【図１０】時間と血流速度の関係を示す図である。
【図１１】位相差のプロファイルを説明するための図である。
【発明を実施するための形態】
【００１０】
＜概要＞
以下、本開示に係る光コヒーレンストモグラフィ装置の概要を図１〜１１に基づいて説明する。本開示の光コヒーレンストモグラフィ装置（例えば、光コヒーレンストモグラフィ装置１０）は、被検体のモーションコントラストデータを取得する。モーションコントラストとは、例えば、被検体の動き、時間的な変化等の検出情報である。例えば、フロー画像等もモーションコントラスト画像の一種とする。なお、フロー画像は、例えば、流体等の動きを検出し、画像化したものである。血液の動きを検出して得られた血管位置を造影した血管造影画像等は、モーションコントラスト画像の一種と言える。
【００１１】
本開示の光コヒーレンストモグラフィ装置（以下、本装置と略す場合がある）は、例えば、ＯＣＴ光学系（例えば、ＯＣＴ光学系２００）と、取得部（例えば、制御部７０）と、を主に備える。ＯＣＴ光学系は、例えば、被検体上を走査部（例えば、光スキャナ１０８）によって走査された測定光と、測定光に対応する参照光とによるＯＣＴ信号（干渉信号）を取得してもよい。取得部は、例えば、被検体上の同一位置に関して、第１の時間間隔を空けて取得された複数のＯＣＴ信号のモーションコントラストを取得する。そして、取得部は、モーションコントラストの基礎となる複数のＯＣＴ信号間での位相差のプロファイルが不連続であった場合、第１の時間間隔とは異なる第２の時間間隔を空けて取得された複数のＯＣＴ信号のモーションコントラストを取得してもよい。取得部は、位相差のプロファイルが不連続にならないように時間間隔を変更してもよい。なお、第２の時間間隔は、第１の時間間隔より短い時間間隔であってもよい。このように、連続的なプロファイルを持つモーションコントラストが取得されることによって、検者は、移動物体の分布等を容易に確認できる。
【００１２】
なお、本装置は、判定部（例えば、制御部７０）をさらに備えてもよい。判定部は、例えば、取得部によって取得された位相差のプロファイルの傾きを検出し、検出された傾きに基づいて、位相差のプロファイルが連続であるか否かを判定してもよい。
【００１３】
なお、本装置は、駆動制御部（例えば、制御部７０）をさらに備えてもよい。駆動制御部は、例えば、位相差のプロファイルが不連続であった場合、走査部の駆動を制御し、同一位置に対する照射の時間間隔を変更してもよい。例えば、駆動制御部は、同一位置に対する照射の時間間隔を短く変更してもよい。
【００１４】
なお、取得部は、取得された位相差のプロファイルが不連続であった場合、第１の時間間隔を空けて取得された複数のＯＣＴ信号から、第１の時間間隔とは異なる第２の時間間隔を空けて取得されたＯＣＴ信号の組を少なくとも１つ選出してもよい。そして、取得部は、選出されたＯＣＴ信号の組のモーションコントラストを取得してもよい。
【００１５】
例えば、取得部は、第１の時間間隔で取得された複数のＯＣＴ信号の内、１つ飛ばしで信号を選出し、選出されたＯＣＴ信号に基づいてモーションコントラストを取得してもよい。これによって、取得部は、例えば、第１の時間間隔の２倍である第２の時間間隔を空けて取得された複数のＯＣＴ信号に基づくモーションコントラストを取得する。同様に、取得部は、２つ、３つ、・・・、と飛ばして信号を選出してもよい。これによって、わざわざ走査部の駆動を制御しなくとも演算上のタイムインターバルＴを変更することができる。
【００１６】
なお、取得部は、例えば、モーションコントラストの基礎となる複数のＯＣＴ信号間での位相差のプロファイルが連続であった場合、モーションコントラストの信号雑音比を大きくするために、第１の時間間隔より長い第２の時間間隔を空けて取得された複数のＯＣＴ信号のモーションコントラストを取得してもよい。
【００１７】
なお、本装置は、血管情報取得部（例えば、制御部７０）をさらに備えてもよい。血管情報取得部は、例えば、被検体の血管の３次元構造に基づいて血管の径および種類の少なくともいずれかを取得する。そして、取得部は、モーションコントラストの基礎となる複数のＯＣＴ信号間での位相差のプロファイルが不連続であった場合、血管情報取得部によって取得された血管の径または種類に基づいて、取得された血管の径または種類に応じて予め設定された、第１の時間間隔とは異なる第２の時間間隔を空けて取得された複数のＯＣＴ信号のモーションコントラストを取得してもよい。
【００１８】
なお、本装置は、例えば、プロセッサ（例えば、制御部７０）等を備え、メモリ等に記憶された制御プログラムを実行してもよい。制御プログラムは、例えば、取得ステップを含む。取得ステップは、例えば、被検体上の同一位置に関して、第１の時間間隔を空けて取得された複数のＯＣＴ信号のモーションコントラストを取得し、取得された位相差のプロファイルが不連続であった場合、第１の時間間隔とは異なる第２の時間間隔を空けて取得された複数のＯＣＴ信号の位相差を取得するステップである。
【００１９】
＜実施例＞
以下、典型的な実施例の１つについて、図面を参照して説明する。図１は本実施例に係る光コヒーレンストモグラフィ装置（以下、本装置と呼ぶ場合もある）１０の構成について説明するブロック図である。本装置１０は、一例として、被検眼の眼底の断層像を取得する眼底撮影装置として説明する。
【００２０】
ＯＣＴ制御系１は、ＯＣＴ光学系１００によって取得された検出信号を処理する。ＯＣＴ制御系１は、制御部７０を有する。ＯＣＴ光学系１００は、例えば、被検眼Ｅの眼底Ｅｆの断層像を撮影する。ＯＣＴ光学系１００は、例えば、制御部７０と接続されている。
【００２１】
ＯＣＴ光学系１００を図２に基づいて説明する。ＯＣＴ光学系１００は、眼底に測定光を照射する。ＯＣＴ光学系１００は、眼底から反射された測定光と参照光との干渉状態を受光素子（検出器１２０）によって検出する。ＯＣＴ光学系１００は眼底Ｅｆ上の撮像位置を変更するため、眼底Ｅｆ上における測定光の照射位置を変更する照射位置変更ユニット（例えば、光スキャナ１０８、固視標投影ユニット３００など）を備える。制御部７０は、設定された撮像位置情報に基づいて照射位置変更ユニットの動作を制御し、検出器１２０からの受光信号に基づいて断層像を取得する。
【００２２】
＜ＯＣＴ光学系＞
ＯＣＴ光学系１００は、いわゆる眼科用光断層干渉計（ＯＣＴ：Optical coherence tomography）の装置構成を持ち、眼Ｅの断層像を撮像する。ＯＣＴ光学系１００は、測定光源１０２から出射された光をカップラー（光分割器）１０４によって測定光（試料光）と参照光に分割する。ＯＣＴ光学系１００は、測定光学系１０６によって測定光を眼Ｅの眼底Ｅｆに導き、また、参照光を参照光学系１１０に導く。ＯＣＴ光学系１００は、眼底Ｅｆによって反射された測定光と、参照光との合成による干渉光を検出器（受光素子）１２０に受光させる。
【００２３】
検出器１２０は、測定光と参照光とのＯＣＴ信号を検出する。フーリエドメインＯＣＴの場合では、干渉光のスペクトル強度（スペクトルＯＣＴ信号）が検出器１２０によって検出され、スペクトル強度データに対するフーリエ変換によってＯＣＴ信号が取得される。例えば、複素ＯＣＴ信号における振幅の絶対値を算出することによって、所定範囲における深さ方向（Ａスキャン方向）のプロファイル（Ａスキャン信号）が取得される。光スキャナ１０８によって走査された測定光の各走査位置における深さ方向の輝度プロファイルを並べることによって、ＯＣＴ画像データ（断層画像データ）が取得される。ここで、ＯＣＴ画像データでの輝度プロファイルとは、Ａ−Ｓｃａｎラインで深さ方向に各画素の輝度値を並べた数列であり、深さ方向に対する輝度値のグラフである。
【００２４】
なお、測定光を２次元的に走査することによって、ＯＣＴ３次元データを取得してもよい。また、ＯＣＴ３次元データから、ＯＣＴ正面（Ｅｎｆａｃｅ）画像（例えば、深さ方向に関して積算された積算画像）が取得されてもよい。
【００２５】
また、ＯＣＴ信号の信号処理によって、機能ＯＣＴ信号が取得されてもよい。光スキャナ１０８によって走査された測定光の各走査位置における機能ＯＣＴ信号（モーションコントラストデータ）を並べることによって、機能ＯＣＴ画像データ（モーションコントラスト画像データ）が取得される。さらに、測定光を２次元的に走査することによって、３次元機能ＯＣＴ画像データ（３次元モーションコントラストデータ）を取得してもよい。また、３次元機能ＯＣＴ画像データから、ＯＣＴ機能正面（Ｅｎ−ｆａｃｅ）画像（例えば、ドップラー正面（Ｅｎ−ｆａｃｅ）画像、スペックルバリアンス正面画像）が取得されてもよい。
【００２６】
ＯＣＴ光学系１００は、例えば、Spectral-domain OCT（ＳＤ−ＯＣＴ）、Swept-source OCT（ＳＳ−ＯＣＴ）が挙げられる。また、Time-domain OCT（ＴＤ−ＯＣＴ）であってもよい。
【００２７】
ＳＤ−ＯＣＴの場合、光源１０２として低コヒーレント光源（広帯域光源）が用いられ、検出器１２０には、干渉光を各周波数成分（各波長成分）に分光する分光光学系（スペクトルメータ）が設けられる。スペクトロメータは、例えば、回折格子とラインセンサからなる。
【００２８】
ＳＳ−ＯＣＴの場合、光源１０２として出射波長を時間的に高速で変化させる波長走査型光源（波長可変光源）が用いられ、検出器１２０として、例えば、単一の受光素子が設けられる。光源１０２は、例えば、光源、ファイバーリング共振器、及び波長選択フィルタによって構成される。そして、波長選択フィルタとして、例えば、回折格子とポリゴンミラーの組み合わせ、ファブリー・ペローエタロンを用いたものが挙げられる。
【００２９】
光源１０２から出射された光は、カップラー１０４によって測定光束と参照光束に分割される。そして、測定光束は、光ファイバーを通過した後、空気中へ出射される。その光束は、光スキャナ１０８、及び測定光学系１０６の他の光学部材を介して眼底Ｅｆに集光される。そして、眼底Ｅｆで反射された光は、同様の光路を経て光ファイバーに戻される。
【００３０】
光スキャナ１０８は、眼底上で二次元的に（ＸＹ方向（横断方向））に測定光を走査させる。光スキャナ１０８は、瞳孔と略共役な位置に配置される。光スキャナ１０８は、例えば、２つのガルバノミラーであり、その反射角度が駆動機構５０によって任意に調整される。
【００３１】
これにより、光源１０２から出射された光束はその反射（進行）方向が変化され、眼底Ｅｆ上で任意の方向に走査される。これにより、眼底Ｅｆ上における撮像位置が変更される。光スキャナ１０８としては、光を偏向させる構成であればよい。例えば、反射ミラー（ガルバノミラー、ポリゴンミラー、レゾナントスキャナ）の他、光の進行（偏向）方向を変化させる音響光学素子（ＡＯＭ）等が用いられる。
【００３２】
参照光学系１１０は、眼底Ｅｆでの測定光の反射によって取得される反射光と合成される参照光を生成する。参照光学系１１０は、マイケルソンタイプであってもよいし、マッハツェンダタイプであっても良い。参照光学系１１０は、例えば、反射光学系（例えば、参照ミラー）によって形成され、カップラー１０４からの光を反射光学系により反射することにより再度カップラー１０４に戻し、検出器１２０に導く。他の例としては、参照光学系１１０は、透過光学系（例えば、光ファイバー）によって形成され、カップラー１０４からの光を戻さず透過させることにより検出器１２０へと導く。
【００３３】
参照光学系１１０は、参照光路中の光学部材を移動させることにより、測定光と参照光との光路長差を変更する構成を有する。例えば、参照ミラーが光軸方向に移動される。光路長差を変更するための構成は、測定光学系１０６の測定光路中に配置されてもよい。
【００３４】
＜正面観察光学系＞
正面観察光学系２００は、眼底Ｅｆの正面観察画像を得るために設けられている。観察光学系２００は、例えば、光源から発せられた測定光（例えば、赤外光）を眼底上で二次元的に走査させる光スキャナと、眼底と略共役位置に配置された共焦点開口を介して眼底反射光を受光する第２の受光素子と、を備え、いわゆる眼科用走査型レーザ検眼鏡（ＳＬＯ）の装置構成を持つ。正面観察光学系２００は、例えば、ＯＣＴ光学系１００によるＯＣＴ信号の取得と並行して、正面観察画像を随時取得してもよい。
【００３５】
なお、観察光学系２００の構成としては、いわゆる眼底カメラタイプの構成であってもよい。また、ＯＣＴ光学系１００は、観察光学系２００を兼用してもよい。すなわち、正面観察画像は、二次元的に得られた断層像を形成するデータを用いて取得されるようにしてもよい（例えば、三次元断層像の深さ方向への積算画像、ＸＹ各位置でのスペクトルデータの積算値、ある一定の深さ方向におけるＸＹ各位置での輝度データ、網膜表層画像、等）。
【００３６】
＜固視標投影ユニット＞
固視標投影ユニット３００は、眼Ｅの視線方向を誘導するための光学系を有する。固視標投影ユニット３００は、眼Ｅに呈示する固視標を有し、複数の方向に眼Ｅを誘導できる。
【００３７】
例えば、固視標投影ユニット３００は、可視光を発する可視光源を有し、視標の呈示位置を二次元的に変更させる。これにより、視線方向が変更され、結果的に撮像部位が変更される。例えば、撮影光軸と同方向から固視標が呈示されると、眼底の中心部が撮像部位として設定される。また、撮影光軸に対して固視標が上方に呈示されると、眼底の上部が撮像部位として設定される。すなわち、撮影光軸に対する視標の位置に応じて撮影部位が変更される。
【００３８】
固視標投影ユニット３００としては、例えば、マトリクス状に配列されたＬＥＤの点灯位置により固視位置を調整する構成、光スキャナを用いて光源からの光を走査させ、光源の点灯制御により固視位置を調整する構成等、種々の構成が考えられる。また、投影ユニット３００は、内部固視灯タイプであってもよいし、外部固視灯タイプであってもよい。
【００３９】
＜制御部＞
制御部７０は、ＣＰＵ（プロセッサ）、ＲＡＭ、ＲＯＭ等を備える。制御部７０のＣＰＵは、各構成の各部材など、装置全体（ＯＣＴ制御系１、ＯＣＴ光学系１００）の制御を司る。ＲＡＭは、各種情報を一時的に記憶する。制御部７０のＲＯＭには、装置全体の動作を制御するための各種プログラム、初期値等が記憶されている。なお、制御部７０は、複数の制御部（つまり、複数のプロセッサ）によって構成されてもよい。
【００４０】
制御部７０には、図１に示すように、不揮発性メモリ（記憶手段）７２、操作部（コントロール部）７６、および表示部（モニタ）７５等が電気的に接続されている。不揮発性メモリ（メモリ）７２は、電源の供給が遮断されても記憶内容を保持できる非一過性の記憶媒体である。例えば、ハードディスクドライブ、フラッシュＲＯＭ、着脱可能なＵＳＢメモリ等を不揮発性メモリ７２として使用することができる。メモリ７２には、ＯＣＴ光学系１００による正面画像および断層画像の撮影を制御するための撮影制御プログラムが記憶されている。また、メモリ７２には、ＯＣＴ制御系１によって得られたＯＣＴ信号を信号処理することを可能にする信号処理プログラムが記憶されている。また、メモリ７２には、走査ラインにおける断層像（ＯＣＴデータ）、三次元断層像（三次元ＯＣＴデータ）、眼底正面像、断層像の撮影位置の情報等、撮影に関する各種情報が記憶される。操作部７６には、検者による各種操作指示が入力される。
【００４１】
操作部７６は、入力された操作指示に応じた信号を制御部７０に出力する。操作部７６には、例えば、マウス、ジョイスティック、キーボード、タッチパネル等の少なくともいずれかを用いればよい。
【００４２】
表示部７５は、装置本体に搭載されたディスプレイであってもよいし、本体に接続されたディスプレイであってもよい。パーソナルコンピュータ（以下、「ＰＣ」という。）のディスプレイを用いてもよい。複数のディスプレイが併用されてもよい。また、表示部７５は、タッチパネルであってもよい。表示部７５がタッチパネルである場合、表示部７５が操作部として機能する。表示部７５には、ＯＣＴ光学系１００によって撮影された断層画像および正面画像を含む各種画像が表示される。
【００４３】
＜モーションコントラスト画像の撮影＞
以下、モーションコントラスト画像を撮影するときの本装置１０の操作方法および制御動作について図３を用いて説明する。制御部７０は、例えば、各種制御処理を司るプロセッサ（例えば、ＣＰＵ）と、プログラムを記憶する記憶媒体とを備える。プロセッサは、プログラムに従って、以下に説明する処理を実行する。なお、以下の説明において、制御の各ステップを識別するための番号を付与するが、付与した番号の順番と実際の制御の順番は必ずしも一致しない。
【００４４】
まず、検者は固視標投影ユニット３００の固視標を注視するように被検者に指示した後、図示無き前眼部観察用カメラで撮影される前眼部観察像を表示部７５で見ながら、被検眼の瞳孔中心に測定光軸がくるように、操作部７６（例えば、図示無きジョイスティック）を用いて、アライメント操作を行う。
【００４５】
（ステップ１：ＯＣＴ撮影）
制御部７０は、同じ位置において、時間の異なる少なくとも２フレームのＯＣＴ信号を取得する。例えば、制御部７０は、光スキャナ１０８の駆動を制御し、眼底上で測定光を走査させる。例えば、図４に示す第１の走査ラインＳ１に沿ってｘ方向に測定光を走査させる。なお、横断方向（例えば、ｘ方向）に測定光を走査させることを「Ｂスキャン」と呼ぶ。以下、１フレームのＯＣＴ信号とは、１回のＢスキャンによって得られたＯＣＴ信号として説明する。制御部７０は、走査中に検出器１２０によって検出されたＯＣＴ信号を取得する。なお、図４において、ｚ軸の方向は、測定光の光軸の方向とする。ｘ軸の方向は、ｚ軸に垂直かつ左右の方向とする。ｙ軸の方向は、ｚ軸に垂直かつ上下の方向とする。
【００４６】
１回目の走査が完了すると、制御部７０は、１回目と同じ位置で２回目の走査を行う。例えば、制御部７０は、図４に示す走査ラインＳ１に沿って測定光を走査させた後、再び走査ラインＳ１に沿って測定光を走査させる。制御部７０は、２回目の走査中に検出器１２０によって検出されたＯＣＴ信号を取得する。これによって、制御部７０は、同じ位置における時間の異なる２フレームのＯＣＴ信号を取得することができる。なお、本実施例においては、同じ位置での走査を４回繰り返し、時間の異なる連続する４フレームのＯＣＴ信号を取得する。制御部７０は、例えば、走査ラインＳ１での走査を４回繰り返し、４フレームのＯＣＴ信号を取得する。ただし、４フレームに限らず、時間の異なる少なくも２フレーム以上のＯＣＴ信号を取得すればよい。
【００４７】
なお、１回の走査で、時間の異なる同じ位置の信号を取得することができる場合は、２回目の走査を行わなくてもよい。例えば、所定間隔だけ光軸のずれた２つの測定光を１度に走査させる場合、複数回走査する必要はない。被検体内の同じ位置における時間の異なるＯＣＴ信号を取得することができればよい。
【００４８】
なお、制御部７０は、別の位置においても、同様に、時間の異なる少なくとも２フレームの信号を取得する。図４に示すように、第１の走査ラインＳ１は、例えば、ｙ＝ｙ１であるとする。また、第２の走査ラインＳ２は、例えば、ｙ＝ｙ２であるとする。制御部７０は、第１の走査ラインＳ１において、時間の異なる信号の取得が完了すると、引き続き、第２の走査ラインＳ２において、時間の異なる少なくとも２フレームの信号を取得する。
【００４９】
制御部７０は、図４に示すように、測定光をラスター走査（ラスタースキャン）し、各走査ラインにおいて、時間の異なる少なくとも２フレーム以上のＯＣＴ信号を得る。これによって、眼底内部の３次元的な情報を取得することができる。
【００５０】
なお、ラスタースキャンは、眼底上を測定光が矩形状に走査するパターンである。ラスタースキャンは、例えば、正面（Ｅｎ−ｆａｃｅ）画像スキャンとして用いられる。
【００５１】
（ステップ２：複素ＯＣＴ信号群の取得）
続いて、制御部７０は、ＯＣＴ光学系１００によって取得されたＯＣＴ信号を処理し、複素ＯＣＴ信号を取得する。例えば、制御部７０は、ステップ１において取得されたＯＣＴ信号（干渉信号）をフーリエ変換する。ここで、Ｎフレーム中ｎ枚目の（ｘ，ｚ）の位置の信号をＡｎ（ｘ，ｚ）で表す。制御部７０は、フーリエ変換によって、複素ＯＣＴ信号Ａｎ（ｘ，ｚ）を得る。複素ＯＣＴ信号Ａｎ（ｘ，ｚ）は、実数成分と虚数成分とを含む。
【００５２】
（ステップ３：血管断層像群の取得）
次いで、制御部７０は、ステップ２によって取得された複素ＯＣＴ信号を処理し、血管断層群（モーションコントラスト画像群）を取得する。複素ＯＣＴ信号を処理する方法としては、例えば、複素ＯＣＴ信号の位相差を算出する方法、複素ＯＣＴ信号のベクトル差分を算出する方法、複素ＯＣＴ信号の位相差及びベクトル差分を掛け合わせる方法などが考えられる。本実施例では、位相差を算出する方法を例に説明する。
【００５３】
まず、制御部７０は、同じ位置の少なくとも２つ以上の異なる時間に取得された複素ＯＣＴ信号Ａ（ｘ，ｚ）に対して位相差を算出する。制御部７０は、例えば、式（１）を用いて、位相の変化を算出する。本実施例では、例えば、４つの異なる時間に測定を行うため、Ｔ１とＴ２、Ｔ２とＴ３、Ｔ３とＴ４、の計３回の計算が行われ、３つのデータが算出される。なお、数式中のＡｎは時間Ｔｎに取得された信号を示し、＊は複素共役を示している。
【００５４】
【数１】
0006798095.tif 000002
【００５５】
制御部７０は、例えば、３フレームの信号を加算平均し、ノイズを除去してもよい。ノイズ成分は各フレームにランダムに存在するため、加算平均することによってシグナル成分に比べて小さくなる。
【００５６】
以上のように、制御部７０は複素ＯＣＴ信号の位相差に関する深さ方向のプロファイルを取得し、このプロファイルの大きさに応じて濃淡をつけることによって、被検者の機能ＯＣＴ画像を取得する。制御部７０は、上記の処理を走査ラインごとに繰り返し、図５に示すように、走査ラインごとに断層血管像を生成する。これによって、制御部７０は、３次元のモーションコントラスト画像を取得する。
【００５７】
＜血流の測定＞
以上のようにして撮影されたモーションコントラスト画像を利用して、血流速度を求める方法を図６に基づいて説明する。本実施例においては、３次元のモーションコントラスト画像から得られる血管の３次元構造を用いて、血流の絶対速度を求める。血流の絶対速度を求めるにあたって、制御部７０は、まず、モーションコントラスト画像と、観察光学系２００によって撮影された観察画像（例えば、ＳＬＯ画像、赤外眼底画像など）の対応付けを行う（ステップ２−１）。なお、モーションコントラスト画像と対応づける観察画像は、ＯＣＴ光学系１００によってＯＣＴ信号を取得しているときに撮影された観察画像でもよいし、モーショントラスト画像が生成された後に撮影された観察画像でもよい。
【００５８】
制御部７０は、例えば、図７に示すように、２次元の観察画像Ｄ２におけるＸＹ方向の座標と、３次元のモーションコントラスト画像Ｄ１におけるＸＹ方向の座標を対応付ける。対応付けの方法としては、例えば、固視灯の点灯位置を利用する方法、血管の位置を利用する方法等が挙げられる。固視灯の点灯位置を利用する場合、例えば、モーションコントラスト画像Ｄ１を取得したときの固視灯の点灯位置と、観察画像Ｄ２を取得したときの固視灯の点灯位置との関係から、モーションコントラスト画像Ｄ１と観察画像Ｄ２を対応付けてもよい。血管の位置を利用する場合、例えば、観察画像Ｄ２に映った血管の位置と、モーションコントラスト画像Ｄ１に映った血管の位置を画像処理によってそれぞれ検出し、血管の位置が一致するように両者を対応付けてもよい。
【００５９】
なお、制御部７０は、例えば、モーションコントラスト画像Ｄ１と観察画像Ｄ２の画素位置を対応付けてもよい。例えば、モーションコントラスト画像Ｄ１のある画素の座標が（ｘ１，ｙ１）であるとすると、それに対応する観察画像Ｄ２の画素の座標（ｘ２，ｙ２）が対応付けられてもよい。
【００６０】
このように、制御部７０は、観察画像Ｄ２とモーションコントラスト画像Ｄ１の対応付けを行っておき、これによって取得された対応情報を血流の絶対速度を求める際に用いてもよい。なお、制御部７０は、対応情報をメモリ７２に記憶させてもよい。
【００６１】
続いて、制御部７０は、図８に示すように、前述の方法で取得されたモーションコントラスト画像Ｄ１を表示部７５の画面上に表示させる（ステップ２−２）。表示部７５に表示されるモーションコントラスト画像Ｄ１は、例えば、３次元画像でもよいし、２次元画像でもよい。図８の例では、３次元データをある深度領域において加算処理して得られた２次元画像が表示されている。この場合、加算処理する深度領域を任意に設定できるようにしてもよい。
【００６２】
検者は、モーションコントラスト画像Ｄ１を確認し、血流を測定する血管を指定する。例えば、検者は、操作部７６のポインティングデバイス（例えば、マウス、タッチパネルなど）のクリックまたはタッチ等によってモーションコントラスト画像Ｄ１の位置を入力することによって血管を指定してもよい。図８の例では、制御部７０は、操作部７６の操作に応じて移動可能なポインタ７６ａを表示している。検者は、例えば、画面上に表示されたポインタ７６ａを所望の血管の上に移動させ、クリック等の操作を行うことによって、血管を指定する。この場合、ポインタ７６ａは、表示部７５上における任意の位置を指定するために用いられる。なお、以下の説明において、検者によって指定された血管を指定血管と呼ぶこととする。
【００６３】
制御部７０は、検者の指示を受け付ける（ステップ２−３）。例えば、制御部７０は、操作部７６から出力された信号に基づいて検者の指示を受け付ける。例えば、制御部７０は、操作部７６のポインティングデバイスによって検者が指定した位置（指定位置）を取得してもよい。
【００６４】
制御部７０は、モーションコントラスト画像上に指定された指定位置が、観察画像Ｄ２のどの位置に対応するのかを取得する（ステップ２−４）。例えば、制御部７０は、メモリ７２に記憶された観察画像Ｄ２とモーションコントラスト画像Ｄ１との対応情報に基づいて、モーションコントラスト画像上の指定位置のｘｙ座標に対応する観察画像上の対応位置を取得する。例えば、制御部７０は、指定血管のｘｙ座標に対応する観察画像上のｘｙ座標をメモリ７２に記憶した対応情報から取得する。
【００６５】
観察画像上での指定血管の対応位置が求まると、制御部７０は、対応位置をトラッキング位置として設定し（ステップ２−５）、スキャンを開始する（ステップ２−６）。制御部７０は、常に対応位置をスキャンできるようにトラッキングしながら、時間間隔を空けて複数回スキャンを行う。例えば制御部７０は、測定光のスキャンと並行して観察画像を随時撮影し、更新された観察画像に映る眼底の位置ずれに応じてスキャン位置を補正してもよい。
【００６６】
より詳細には、制御部７０は、モーションコントラスト画像Ｄ１との対応付けに用いた観察画像Ｄ２の静止画と、現在の観察画像Ｄ２と、を比較して、位置ずれ方向及び位置ずれ量を画像処理により検出（演算）する。例えば、制御部７０は、モーションコントラスト画像Ｄ１との対応付け時における観察画像Ｄ２の静止画データを基準画像とし、その基準画像とリアルタイムで取得される観察画像との位置ずれ方向及び位置ずれ量を算出する。これにより、静止画像に対する位置ずれ情報が得られる。
【００６７】
上記のようにして、位置ずれ方向及び位置ずれ量が検出されると、制御部７０は、走査位置のずれが補正されるように、光スキャナ１０８の２つのガルバノミラーを適宜駆動制御する。これによって、走査位置が補正される。以上のようにして、被検眼がずれた場合であっても、走査位置が補正され、常時、制御部７０は、同じスキャン位置において時間の異なる複数のＯＣＴ信号を取得することができる。
【００６８】
制御部７０は、ステップ２−６において取得された複数回のＯＣＴ信号と、３次元のモーションコントラスト画像から得られた血管の３次元構造を用いて、血流の絶対速度を算出する（ステップ２−７）。例えば、制御部７０は、複数回のＯＣＴ信号からドップラー位相シフトを求める。そして、制御部７０は、求められた位相差と血管の３次元構造から得られた血流方向から血流の絶対速度を算出してもよい。
【００６９】
モーションコントラスト画像Ｄ１から得られた血管の３次元構造から血流の方向を求める方法としては、例えば、血管の３次元構造を細線化する方法が挙げられる。例えば、図９に示すように、細線化の処理によって、３次元の血管構造Ｋを１本の線Ｌ０で表すことができる。制御部７０は、例えば、スキャンラインＳＬ上を走査される測定光に対する線Ｌ０の傾きを血流の方向とみなし、前述の位相差から血流の絶対速度を算出してもよい。
【００７０】
＜血流速度の算出＞
血流の絶対速度について説明する。血流の絶対速度vの光軸方向の成分vｚは、測定光と血管のなす角度αを用いて次式（２）で表せる。
【００７１】
【数２】
0006798095.tif 000003
【００７２】
ここで、成分vzは、位相差ΔΦDoppler（−π〜π）、血管組織の屈折率ｎ、中心波数ｋ、タイムインターバルＴを用いて次式（３）で表せる。
【００７３】
【数３】
0006798095.tif 000004
【００７４】
以上の２つの式を用いて、次式（４）のように位相差から血流の絶対速度vが求められる。
【００７５】
【数４】
0006798095.tif 000005
【００７６】
以上のように、本実施例において、制御部７０は、３次元のモーションコントラスト画像Ｄ１と、同じ走査位置に関する時間的に異なる複数のＯＣＴ信号とを用いて被検眼の血流を測定する。制御部７０は、３次元のモーションコントラスト画像Ｄ１によって得られる血管の３次元構造から血流方向を取得できるため、ドップラーシフトから被検眼の絶対的な血流速度を測定することができる。これによって、制御部７０は、例えば、複数の被検者に関して血流速度を比較することによって、被検眼の診断等に役立つ情報を取得できる。
【００７７】
また、本実施例において、制御部７０は、予め取得された３次元モーションコントラスト画像Ｄ１から得られた血流方向を血流の絶対速度の算出に用いる。このため、ＯＣＴ光学系によって被検眼の同一位置において時間的に異なる複数のＯＣＴ信号を継続して取得すれば、制御部７０は、被検者の拍動の影響を考慮した動的な血流の絶対速度を取得することができる（図１０参照）。このように、制御部７０は、時間的な血流速度の変化を取得することによって、動脈硬化等の疾患を発見するのに利用可能な情報を提供できる。
【００７８】
さらに、本実施例において、制御部７０は、モーションコントラスト画像Ｄ１で確認した血管の位置を観察画像Ｄ２に対応付け、その対応情報に基づいてトラッキングを行う。例えば、モーションコントラスト画像Ｄ１は、画像を生成するまでに時間を要するため、被検眼の微動に対して追従するためのずれ情報を取得することが難しい。一方、観察画像Ｄ２は、画像を取得するまでの時間がモーションコントラスト画像に比べて短いため、被検眼の微動に対して追従するためのずれ情報を取得するのに適している。本実施例の制御部７０のように、観察画像Ｄ２では確認しづらい細かい血管の位置をモーションコントラスト画像から対応付けることによって、観察画像Ｄ２では確認しづらい細かい血管であってもトラッキングを行うことができる。これによって、制御部７０は、同一位置における時間の異なる複数のＯＣＴ信号を継続的に取得することが容易となり、被検者の動的な血流の絶対速度を適切に求めることができる。
【００７９】
なお、制御部７０は、検者によって指定された血管をスキャンするときのスキャン幅を、モーションコントラスト画像から求めた血管の径に応じて自動で設定してもよい。もちろん、検者が手動で設定してもよい。
【００８０】
なお、制御部７０は、検者によって指定された血管をスキャンするときのタイムインターバルＴを、モーションコントラスト画像から求めた血管の径に応じて設定してもよい。もちろん、検者が手動で設定してもよい。
【００８１】
なお、検者によって指定された血管をスキャンするときのスキャン方向は、血流方向に対して垂直であってもよい。例えば、制御部７０は、Ｅｎ−ｆａｃｅ（測定光に対して垂直な方向）において、指定された血管の血流方向に対して垂直な方向を算出し、この方向にスキャンを行ってもよい。血管の血流方向は、前述のように血管の３次元構造を細線化することで求めてもよい。細線化には、例えば、モルフォロジ処理、距離変換処理、またはHilditch、Deutschなどの既存アルゴリズムを用いてもよい。
【００８２】
このように、制御部７０は、血流方向に対して垂直にスキャンを行うことによって、血流を算出する際に血流方向に対するスキャン角度を考慮する必要がなくなり、計算処理を簡単にすることができる。制御部７０は、例えば、指定血管の半径方向を自動的にスキャン方向に設定してもよい。もちろん、制御部７０は、検者の操作に応じて操作部７６から出力された操作信号に基づいてスキャンラインを設定してもよい。
【００８３】
なお、図８に示すようなモーションコントラスト画像において、検者が新たな血管を指定する度に、制御部７０は前述のように指定された血管をトラッキングしながらスキャンを行い、血流速度を求めてもよい。
【００８４】
＜位置ずれ検出手法＞
なお、トラッキングの際に２つの観察画像間の位置ずれを検出する手法としては、種々の画像処理手法（各種相関関数を用いる方法、フーリエ変換を利用する方法、特徴点のマッチングに基づく方法）を用いることが可能である。
【００８５】
例えば、所定の基準画像（例えば、過去の正面画像）又は対象画像（現在の正面画像）を１画素ずつ位置ずれさせ、基準画像と対象画像を比較し、両データが最も一致したとき（相関が最も高くなるとき）の両データ間の位置ずれを検出する手法が考えられる。また、所定の基準画像及び対象画像から共通する特徴点を抽出し、抽出された特徴点の位置ずれを検出する手法が考えられる。
【００８６】
また、２つの画像間の位置ずれを求めるための関数として、位相限定相関関数を用いるようにしてもよい。この場合、まず、各画像をフーリエ変換し、各周波数成分の位相と振幅を得る。なお、得られた振幅成分は、各周波数成分に関して大きさ１に正規化しておく。次に、２つの画像間で周波数毎の位相差を算出した後、これらに逆フーリエ変換をかける。
【００８７】
ここで、２つの画像間の位置ずれがなければ、余弦波のみの加算となり、原点位置（０，０）にピークが出現する。また、位置ずれがある場合、位置ずれに対応する位置にピークが出る。そこで、ピークの検出位置を求めることにより２つの画像間の位置ずれが得られる。この手法によれば、正面画像の位置ずれを高精度かつ短時間で検出できる。
【００８８】
なお、前述の式（２）で説明したように、位相差ΔΦDopplerは、逆正接演算により、−πからπの範囲内である。したがって、各走査間の時間間隔（タイムインターバルＴ）が長い場合、血流が速いと位相差ΔΦDopplerが−πからπの範囲を超えてしまい、例えば、図１１（ａ）のように、スキャン方向の位相差ΔΦDopplerのプロファイルが不連続になる。このような場合、検出されたプロファイルから速度測定を含めた被検眼の診断を行うことは難しい。そこで、本実施例のように、例えば、制御部７０は、適切なタイムインターバルＴを探索してスキャンを行う。
【００８９】
以下、ステップ２−６のスキャン信号を取得する処理において、位相差ΔΦDopplerが−π〜πの範囲を超えないように適切なタイムインターバルＴでスキャンを行う方法について説明する。なお、以下の説明において、位相差ΔΦDopplerが−π〜πの範囲を超えて不連続点が生じることを「ラップされる」と表現する。
【００９０】
適切なインターバルでスキャンを行う方法としては、例えば、制御部７０は、タイムインターバルＴを走査毎に順次変更しながらスキャンを行う方法が挙げられる。以下の説明において、制御部７０は、タイムインターバルＴを長い間隔から短い間隔に順次変更しながらスキャンを行う。この場合、制御部７０は、取得された複数のＯＣＴ信号を用いてタイムインターバルＴの異なる位相差ΔΦDopplerを順次算出し、位相差ΔΦDopplerの算出結果がラップラップされるか否かを判定する。そして、制御部７０は、例えば、位相差ΔΦDopplerがラップされなくなったと判定したときのタイムインターバルＴ、またはそれよりも短いタイムインターバルＴに固定してスキャンを続ける。
【００９１】
位相差ΔΦDopplerがラップされるか否か判定する方法としては、例えば、隣り合う画素間における位相差ΔΦDopplerの変化の大きさに基づいて判定する方法等が挙げられる。この場合、制御部７０は、複数のＯＣＴ信号から取得された位相差ΔΦDopplerに対してエッジ検出フィルタ（例えば、ラプラシアンフィルタなどの微分処理）を適用してもよい。例えば、図１１（ａ）に示すように、ラップされている位相差ΔΦDopplerのプロファイルに対してエッジ検出フィルタを適用した場合、図１１（ｂ）に示すように、ラップされている位置の強度が局所的に大きくなる。一方、図１１（ｃ）に示すように、ラップされていない位相差ΔΦDopplerのプロファイルに対してエッジ検出フィルタを適用した場合、図１１（ｄ）に示すように、強度の局所的に大きくなることはない。従って、制御部７０は、例えば、エッジ検出フィルタ適用後のプロファイルの平均値が所定以下になったときに、ラップされなくなったと判定してもよい。もちろんエッジ検出フィルタ適用後の最大値が所定以下になったときにラップされなくなったと判定してもよい。
【００９２】
以上のように、制御部７０は、タイムインターバルＴを適切に設定することによって、血流の速さの違いに関わらずラップされていない位相差ΔΦDopplerを取得するこができる。ラップされていない位相差ΔΦDopplerを取得することによって血流の層流状態を良好に取得でき、被検眼の診断等に利用できる。また、検者が自らタイムインターバルＴの設定を行う手間が省ける。
【００９３】
なお、制御部７０は、できるだけ位相差ΔΦDopplerの信号強度が大きくなるようなタイムインターバルＴを探索してもよい。例えば、制御部７０は、位相差ΔΦDopplerがラップされなくなったと判定した場合、−π〜πの範囲内で位相差ΔΦDopplerのプロファイルの最大値ができるだけ大きくなるように、タイムインターバルＴを徐々に長く変更し、適切なタイムインターバルＴを探索してもよい。
【００９４】
なお、上記の説明において、制御部７０は、タイムインターバルＴを長い間隔から短い間隔に順次変更するものとした。タイムインターバルＴが大きい場合、血流によるドップラー信号の位相差ΔΦDopplerが大きくなり、ラップされる可能性が大きい。このように、わざとラップ状態から徐々にタイムインターバルＴを短くし、ラップされない状態にした場合、−π〜πの範囲で強度の大きい信号を取得できる可能性が大きい。
【００９５】
もちろん、制御部７０は、タイムインターバルＴを短い間隔から長い間隔に変更させてもよい。タイムインターバルＴが短いと、血流によるドップラー信号の位相差ΔΦDopplerが小さく、ラップされない位相差ΔΦDopplerを素早く取得することができる。ただし、タイムインターバルＴが短いと、図１１（ｅ），（ｆ）に示すように、信号強度が小さくなり、信号がノイズフロアに埋もれる可能性があるため、制御部７０は、ラップされない範囲でできるだけ長いタイムインターバルＴを探索してもよい。
【００９６】
なお、前述の例では、タイムインターバルＴを徐々に長く変更する、または徐々に短く変更するものとしたが、これに限らない。例えば、制御部７０は、タイムインターバルＴを長い間隔と短い間隔とで交互に変更してもよいし、ランダムにタイムインターバルＴの間隔を変更してもよい。なお、タイムインターバルＴを短くするか長くするかは、検者の操作によって操作部から出力された操作信号によって設定できるようにしてもよい。
【００９７】
なお、制御部７０は、３次元モーションコントラスト画像Ｄ１を生成するために取得した複数のＯＣＴ信号を用いて、ある時点での血流の絶対速度を求め、その大きさに基づいて、最初のタイムインターバルＴを設定してもよい。その後、前述のように、タイムインターバルＴを順次変更することで適切なタイムインターバルＴを求めてもよい。このように、３次元モーションコントラスト画像Ｄ１を取得した時点での血流の絶対速度からタイムインターバルＴを設定することで、適切なタイムインターバルＴを設定するまでに要する時間を短縮できる。
【００９８】
なお、血流速度の大きさは、血管の太さ、種類（例えば、動脈または静脈など）等によってある程度推定することができる。したがって、制御部７０は、指定血管の太さ、種類等から推定される血流速度に基づいて、タイムインターバルＴを設定してもよい。この場合、例えば、制御部７０は、血管の太さおよび種類等と、血流速度の大きさとの対応関係を過去の測定データなどから取得し、対応表としてメモリ７２に記憶させておく。そして制御部７０は、指定された血管の大きさ、および血管の存在する網膜の深度等から推定される血管の種類の少なくともいずれかを、メモリ７２に記憶させた対応表と照合して指定血管の血流速度を推定してもよい。
【００９９】
なお、制御部７０は、例えば、複数の異なるタイムインターバルＴで複数のＯＣＴ信号の取得を行い、得られた複数のＯＣＴ信号から算出されたΔΦDopplerの内、ラップされていないものを血流速度の算出に用いてもよい。
【０１００】
なお、以上の説明において、例えば、制御部７０はタイムインターバルＴを変更するものとしたが、これに限らない。例えば制御部７０は、一定のタイムインターバルＴでスキャンを行い、得られた複数のＯＣＴ信号の内、算出される位相差ΔΦDopplerがラップされていないＯＣＴ信号の組み合わせを選択して血流の絶対速度を算出してもよい。例えば、制御部７０は、短いタイムインターバルＴでスキャンを行い、取得された複数のＯＣＴ信号からデータを間引いてΔΦDopplerを算出してもよい。例えば、０．５ｍ秒のタイムインターバルＴでスキャンを行った場合、取得された複数のＯＣＴ信号からデータを間引いて、０．５の倍数となるタイムインターバルＴ（例えば、１．０ｍ秒、１．５ｍ秒、２．０ｍ秒など）で取得されたＯＣＴ信号からΔΦDopplerを算出してもよい。この場合、制御部７０は、例えば各タイムインターバルＴにおいて取得されたＯＣＴ信号の内、ΔΦDopplerがラップされていないものを血流速度の算出に用いてもよいし、最大値が大きいものを血流速度の算出に用いてもよい。さらに、各タイムインターバルＴにおいて取得されたＯＣＴ信号の内、ΔΦDopplerの傾き（例えば、微分値）の最小値が小さいものを血流速度の算出に用いてもよい。
【０１０１】
このように、実際に測定光を照射するタイムインターバルＴは変更せずに、演算上のタイムインターバルＴを変更することによって適切なタイムインターバルＴを探索してもよい。
【符号の説明】
【０１０２】
１  光コヒーレンストモグラフィデバイス
７０  制御部
７２  メモリ
７５  モニタ
７６  操作部
１００  干渉光学系（ＯＣＴ光学系）
１０８  光スキャナ
１２０  検出器
２００  正面観察光学系
３００  固視標投影ユニット
【図１】
0006798095.tif 000006
【図２】
0006798095.tif 000007
【図３】
0006798095.tif 000008
【図４】
0006798095.tif 000009
【図５】
0006798095.tif 000010
【図６】
0006798095.tif 000011
【図７】
0006798095.tif 000012
【図８】
0006798095.tif 000013
【図９】
0006798095.tif 000014
【図１０】
0006798095.tif 000015
【図１１】
0006798095.tif 000016