特許 6230023 ３次元光コヒーレンス断層撮影インターフェログラムデータから２次元画像を生成する装置および方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6230023

(24)【登録日】2017年10月27日

(45)【発行日】2017年11月15日

(54)【発明の名称】３次元光コヒーレンス断層撮影インターフェログラムデータから２次元画像を生成する装置および方法

(51)【国際特許分類】

   A61B 3/10 20060101AFI20171106BHJP

   G01N 21/17 20060101ALI20171106BHJP

【ＦＩ】

   A61B3/10 RZDM

   G01N21/17 630

【請求項の数】15

【全頁数】16

(21)【出願番号】特願2014-38896(P2014-38896)

(22)【出願日】2014年2月28日

(65)【公開番号】特開2014-188373(P2014-188373A)

(43)【公開日】2014年10月6日

【審査請求日】2016年12月6日

(31)【優先権主張番号】13/851,612

(32)【優先日】2013年3月27日

(33)【優先権主張国】US

(73)【特許権者】

【識別番号】000220343

【氏名又は名称】株式会社トプコン

(74)【代理人】

【識別番号】100112427

【弁理士】

【氏名又は名称】藤本 芳洋

(72)【発明者】

【氏名】ライスマン，チャールズ エー．

【審査官】 九鬼 一慶

(56)【参考文献】

【文献】 特開２００７−１１７７１４（ＪＰ，Ａ）

【文献】 特開２００８−１７５６９８（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ａ６１Ｂ ３／００−３／１８

Ｇ０１Ｎ ２１／１７

(57)【特許請求の範囲】

【請求項1】

光コヒーレンス断層撮影インターフェログラムデータから強度マップを得るための装置であって、
光ビームを生成するための光源と、
前記光ビームをサンプル経路と参照経路とに沿って分割するビームスプリッタと、
前記サンプル経路に沿って位置する、前記光ビームをＸ−Ｙ平面内の異なる位置に向けるためのスキャナと、
前記サンプル経路および前記参照経路の各々から返される光を受け取り、複数の出力セットを生成するための検出器であって、前記複数の出力セットの各々は、前記光ビームが特定のＸ−Ｙ平面位置に向けられるときに前記光源の異なる波長で受け取られる光強度に対応し、前記光強度は、前記特定のＸ−Ｙ平面位置における深さ方向Ｚの対象物内の光反射率分布に関する情報を含むものである、前記検出器と、
コンピュータプログラム命令を記憶するメモリと、
前記メモリに通信可能に結合された、前記コンピュータプログラム命令を実行するように構成されたプロセッサであって、前記コンピュータプログラム命令は、前記プロセッサ上で実行されると前記プロセッサに、
前記光ビームを特定のＸ−Ｙ平面位置に向けることにより生成された出力セットを高域フィルタリングしてフィルタ出力セットを生成するステップと、
前記フィルタ出力セットを、逆累積分布関数に基づいて、前記特定のＸ−Ｙ平面位置についての前記深さ方向Ｚの単一の推定強度値に変換するステップであって、１つもしくは複数のＸ−Ｙ平面位置に対応する１つもしくは複数の単一の推定強度値が前記対象物の２次元画像強度マップを生成するのに適する、前記変換するステップと
を含む方法を実行させるものである、前記プロセッサと
を備える装置。

【請求項2】

前記方法は、前記フィルタ出力セットを、前記対象物の２次元画像強度マップのリアルタイム表示を生成するための単一の推定強度値に変換するステップをさらに含む、請求項１に記載の装置。

【請求項3】

前記方法は、前記対象物の前記２次元画像強度マップをディスプレイにおいて提示するステップをさらに含む、請求項１に記載の装置。

【請求項4】

前記２次元画像強度マップの少なくとも一部が、画像診断法と関連付けられた走査取り込み整合プロセスに使用される、請求項１に記載の装置。

【請求項5】

前記フィルタ出力セットを単一の推定強度値に変換する前記ステップは、
少なくとも１つの選択パーセンタイルに対応する前記フィルタ出力セットのうちの少なくとも１つの出力を選択するステップと、
前記少なくとも１つの出力を前記単一の推定強度値に変換するステップと
をさらに含む、請求項１に記載の装置。

【請求項6】

前記２次元画像強度マップは、別の画像診断法によって得られた画像を位置合わせするのに使用される、請求項１に記載の装置。

【請求項7】

前記光源は可同調掃引光ビームを生成し、前記光源の異なる波長における光強度が経時的に得られる、請求項１に記載の装置。

【請求項8】

前記方法は、前記フィルタ出力セットをダウンサンプリングするステップをさらに含む、請求項１に記載の装置。

【請求項9】

前記方法は、前記フィルタ出力セットを切り捨てるステップをさらに含む、請求項１に記載の装置。

【請求項10】

前記フィルタ出力セットを変換する前記ステップは、範囲、偏り、標準偏差、分散、またはエントロピーのうちの１つを測定するステップおよび推定するステップのうちの一方を含む、請求項１に記載の装置。

【請求項11】

光コヒーレンス断層撮影インターフェログラムデータから強度マップを得るための装置であって、
光ビームを生成するための光源と、
前記光ビームをサンプル経路と参照経路とに沿って分割するビームスプリッタと、
前記サンプル経路に沿って位置する、前記光ビームをＸ−Ｙ平面内の異なる位置に向けるためのスキャナと、
前記サンプル経路および前記参照経路の各々から返される光を受け取り、複数の出力セットを生成するための検出器であって、前記複数の出力セットの各々は、前記光ビームが特定のＸ−Ｙ平面位置に向けられるときに前記光源の異なる波長で受け取られる光強度に対応し、前記光強度は、前記特定のＸ−Ｙ平面位置における深さ方向Ｚの対象物内の光反射率分布に関する情報を含むものである、前記検出器と、
コンピュータプログラム命令を記憶するメモリと、
前記メモリに通信可能に結合された、前記コンピュータプログラム命令を実行するように構成されたプロセッサであって、前記コンピュータプログラム命令は、前記プロセッサ上で実行されると前記プロセッサに、
出力セットを、逆累積分布関数に基づいて、前記特定のＸ−Ｙ平面位置についての前記深さ方向Ｚの単一の推定強度値に変換するステップであって、１つもしくは複数のＸ−Ｙ平面位置に対応する１つもしくは複数の単一の推定強度値が前記対象物の２次元画像強度マップを生成するのに使用される、前記単一の推定強度値に変換するステップと、
前記出力セットを、前記特定のＸ−Ｙ平面位置についての深さ方向Ｚの強度情報を有する３次元データセットに変換するステップであって、１つもしくは複数の３次元データセットが前記対象物の１つもしくは複数のＯＣＴ断面画像を生成するのに使用される、前記３次元データセットに変換するステップと
を含む方法を実行させるものである、前記プロセッサと
を備える装置。

【請求項12】

前記対象物の前記２次元画像強度マップは、別の画像診断法によって得られた前記対象物の１つもしくは複数のＯＣＴ断面画像を位置合わせするのに使用される、請求項１１に記載の装置。

【請求項13】

前記方法は、前記対象物の前記２次元画像強度マップを前記対象物の１つもしくは複数のＯＣＴ断面画像と並列に表示するステップをさらに含む、請求項１１に記載の装置。

【請求項14】

前記方法は、１つもしくは複数のランドマークに基づいて前記３次元データセットをセグメント化するステップをさらに含む、請求項１１に記載の装置。

【請求項15】

光コヒーレンス断層撮影インターフェログラムデータから強度マップを得るための方法であって、
光ビームを光源により生成するステップと、
前記光ビームをサンプル経路と参照経路とに沿って分割するステップと、
前記光ビームを前記サンプル経路に沿ってＸ−Ｙ平面内の異なる位置に向けるステップと、
前記サンプル経路および前記参照経路の各々から返された光を受け取るステップと、
複数の出力セットを生成するステップであって、前記複数の出力セットの各々は、前記光ビームが特定のＸ−Ｙ平面位置に向けられるときに前記光源の異なる波長で受け取られる光強度に対応し、前記光強度は、前記特定のＸ−Ｙ平面位置における深さ方向Ｚの対象物内の光反射率分布に関する情報を含むものである、前記複数の出力セットを生成するステップと、
前記光ビームを特定のＸ−Ｙ平面位置に向けることにより生成された出力セットを高域フィルタリングしてフィルタ出力セットを生成するステップと、
前記フィルタ出力セットを、逆累積分布関数に基づいて、前記特定のＸ−Ｙ平面位置についての前記深さ方向Ｚの単一の推定強度値に変換するステップであって、１つもしくは複数のＸ−Ｙ平面位置に対応する１つもしくは複数の単一の推定強度値が前記対象物の２次元画像強度マップを生成するのに適する、前記変換するステップと
を含む方法。

【発明の詳細な説明】

【技術分野】

【0001】

本開示は、一般に、コヒーレント波形に基づく撮像を対象とし、より詳細には、光コヒーレンス断層撮影法の撮像システムに関する。

【背景技術】

【0002】

光コヒーレンス断層撮影法（Ｏｐｔｉｃａｌ ｃｏｈｅｒｅｎｃｅ ｔｏｍｏｇｒａｐｈｙ：ＯＣＴ）は、光信号の取得および処理の方法である。ＯＣＴは、表面下の画像が得られるように観察対象の光散乱物体（生物組織など）の表面を越えて貫通させるのに使用することができる干渉に基づく技法である。

【0003】

ＯＣＴシステムは、ある特定の医学的状態を診断し、モニタするのに十分な解像度のものである断面画像または断面画像の集合体（３次元画像など）を提供することができる。従来、ＯＣＴシステムは、観察対象物から受け取られる後方散乱および後方反射光からのエコー時間遅延の測定値を提供することによって動作する。そのようなＯＣＴシステムは通常、干渉計および機械的に走査される参照光路とを含み、それらのシステムを一般に時間領域ＯＣＴと呼ぶ。

【0004】

スペクトル領域ＯＣＴシステムまたは波長掃引光源型フーリエ領域ＯＣＴシステムは、観察対象物から測定される光と固定参照経路からの光との間の干渉スペクトルからの光のエコー時間遅延の測定値を提供することによって動作する。スペクトル領域ＯＣＴシステムは通常、光分散性部品と、観察対象物から受け取られる干渉スペクトルを測定するための、電荷結合デバイス（ＣＣＤ）カメラといった検出器アレイとからなる分光計を含む。一方、波長掃引光源型システムは通常、高速波長同調レーザーと、干渉スペクトルを測定するための検出・高速データ取得装置とを含む。どちらの種類のシステムでも、観察対象物からの後方散乱および後方反射光のエコー時間遅延は、干渉スペクトルのフーリエ変換を計算することによって求められる。

【0005】

フーリエ領域ＯＣＴシステムは時間領域ＯＣＴシステムに優る改善である。というのは、観察対象物の異なる軸位置における後方散乱および後方反射光を、順次にではなく同時に測定することができるからである。したがって、ある特定の医学的状態を診断し、モニタするための撮像速度および感度が改善されている。しかし、より効率的で正確な医学的診断、モニタリング他の機能を提供するには、フーリエ領域ＯＣＴデータを使用した測定法におけるさらなる改善が有益となるはずである。

【0006】

さらに、眼科ＯＣＴシステムは通常、走査レーザー検眼鏡や赤外線眼底カメラといった２次的画像診断法を利用して全般的診断のために眼底を撮像し、正面眼底像を作成することによってＯＣＴ走査整合を補助する。しかし、被験者または患者においては、小さい瞳孔や眼球混濁といった撮像状況が存在し、その場合には、たとえ眼科ＯＣＴシステムが許容できる画像を生成しうるとしても、すべての２次的画像診断法が有効に働くとは限らない。

【発明の概要】

【発明が解決しようとする課題】

【0007】

したがって、２次的画像診断法を置き換え、またはその代用をするように、ＯＣＴ走査データを利用して正面眼底画像の表示を生成する方法を作成することが望ましい。正面眼底画像の表示は、臨床家がより確実に意図される位置を走査することができる能力を提供し、また、他のデータセットもしくは画像診断法との位置合わせにも、特徴セグメンテーション計算を可能にするためにも使用することのできるＯＣＴ走査データの本来的に重ね合わされたマップも提供する。

【課題を解決するための手段】

【0008】

３次元ＯＣＴインターフェログラムデータから２次元眼底画像を生成するための方法および装置が提供される。一実施形態によれば、光源からの光ビームがサンプル経路と参照経路とに沿って分割される。光ビームはサンプル経路に沿ってＸ−Ｙ平面内の異なる位置に向けられる。サンプル経路および参照経路の各々から返された光が受け取られる。複数の出力セットが生成され、複数の出力セットの各々は、光ビームが特定のＸ−Ｙ平面位置に向けられるときに光源の異なる波長で受け取られる光強度に対応し、光強度は、特定のＸ−Ｙ平面位置における深さ方向Ｚの対象物内の光反射率分布に関する情報を含む。光ビームを特定のＸ−Ｙ平面位置に向けることにより生成される出力セットが高域フィルタリングされてフィルタ出力セットが生成される。フィルタ出力セットは、任意選択で、ダウンサンプリングされ、または切り捨てられてもよい。フィルタ出力セットは、逆累積分布関数に基づいて、特定のＸ−Ｙ平面位置についての深さ方向Ｚの単一の推定強度値に変換され、その場合、１つもしくは複数のＸ−Ｙ平面位置に対応する１つもしくは複数の単一の推定強度値が、対象物の２次元画像強度マップを生成するのに適する。フィルタ出力セットは二乗され、またはソートされてよい。光源は広帯域光ビームを生成してよく、検出器は、回折格子と検出器アレイとを含む分光計とすることができる。また光源は、可同調掃引光ビームを生成してもよく、光源の異なる波長での光強度が時間の関数として得られてもよい。

【0009】

一実施形態によれば、フィルタ出力セットはリアルタイムで単一の推定強度値に変換されてよい。

【0010】

一実施形態によれば、対象物の２次元画像強度マップはディスプレイにおいて提示されてよい。

【0011】

一実施形態によれば、対象物の２次元画像強度マップの少なくとも一部が、画像診断法と関連付けられた走査取り込み整合プロセスに使用されてよい。また対象物の２次元画像強度マップは、別の画像診断法によって得られた対象物の１つもしくは複数の画像を位置合わせするのに使用されてもよい。画像診断法は、ＯＣＴまたは眼底撮像カメラ、走査レーザー検眼鏡、眼底自発蛍光または血管造影法のうちの１つとすることができる。

【0012】

一実施形態によれば、フィルタ出力セットを単一の推定強度値に変換することは、少なくとも１つの選択パーセンタイル値に対応するフィルタ出力セットのうちの少なくとも１つの出力を選択すること、少なくとも１つの出力を単一の推定強度値に変換することをさらに含む。少なくとも１つの選択パーセンタイル値は、最小値、中間値もしくは最大値のうちの１つに対応してもよく、フィルタ出力セット内の任意の事前選択パーセンタイル値としてもよい。

【0013】

一実施形態によれば、光源からの光ビームがサンプル経路と参照経路とに沿って分割される。光ビームはサンプル経路に沿ってＸ−Ｙ平面内の異なる位置に向けられる。サンプル経路および参照経路の各々から返された光が受け取られる。複数の出力セットが生成され、複数の出力セットの各々は、光ビームが特定のＸ−Ｙ平面位置に向けられるときに光源の異なる波長で受け取られる光強度に対応し、光強度は、特定のＸ−Ｙ平面位置における深さ方向Ｚの対象物内の光反射率分布に関する情報を含む。出力セットが高域フィルタリングされ、フィルタ出力セットは、逆累積分布関数に基づいて、特定のＸ−Ｙ平面位置についての深さ方向Ｚの単一の推定強度値に変換され、その場合、１つもしくは複数のＸ−Ｙ平面位置に対応する１つもしくは複数の単一の推定強度値が、対象物の２次元画像強度マップを生成するのに適し、出力セットは、特定のＸ−Ｙ平面位置についての深さ方向Ｚの強度情報を有する３次元データセットに変換され、その場合、１つもしくは複数の３次元データセットが対象物の１つもしくは複数のＯＣＴ断面画像を生成するのに使用される。３次元データセットは、１つもしくは複数のランドマークに基づいてセグメント化されてよく、ランドマークは深さ方向Ｚの強度情報から識別されうる。ランドマークは、対象物の１つもしくは複数の物理境界を含みうる。

【0014】

一実施形態によれば、対象物の２次元画像強度マップは、対象物の１つもしくは複数のＯＣＴ断面画像を位置合わせするのに使用されてよく、対象物の２次元画像強度マップは、対象物の１つもしくは複数のＯＣＴ断面画像と並列に表示されうる。

【0015】

一実施形態によれば、対象物の眼球マップが、２次元画像強度マップおよびセグメント化３次元データセットのうちの１つに基づいて生成されうる。加えて、対象物の液体充満空間が、セグメント化３次元データセットに基づいて定量化されてもよい。

【0016】

一実施形態によれば、光源からの光ビームがサンプル経路と参照経路とに沿って分割される。光ビームはサンプル経路に沿ってＸ−Ｙ平面内の異なる位置に向けられる。サンプル経路および参照経路の各々から返された光が受け取られる。複数の出力セットが生成され、複数の出力セットの各々は、光ビームが特定のＸ−Ｙ平面位置に向けられるときに光源の異なる波長で受け取られる光強度に対応し、光強度は、特定のＸ−Ｙ平面位置における深さ方向Ｚの対象物内の光反射率分布に関する情報を含む。光ビームを特定のＸ−Ｙ平面位置に向けることにより生成される出力セットが高域フィルタリングされてフィルタ出力セットが生成される。フィルタ出力セットは、任意選択で、ダウンサンプリングされ、または切り捨てられてもよい。フィルタ出力セットは、特定のＸ−Ｙ平面位置についての深さ方向Ｚの単一の推定強度値に変換され、その場合、１つもしくは複数のＸ−Ｙ平面位置に対応する１つもしくは複数の単一の推定強度値が、対象物の２次元画像強度マップを生成するのに適し、対象物の２次元画像強度マップの少なくとも一部が、画像診断法と関連付けられた走査取り込み整合プロセスに使用されうる。

【0017】

本発明の上記その他の利点は、以下の詳細な説明および添付の図面を参照すれば当業者には明らかになるであろう。

【図面の簡単な説明】

【0018】

【図1】一実施形態による、３次元ＯＣＴインターフェログラムデータから２次元画像を生成するのに使用されうるシステムを示す図である。

【0019】

【図2】一実施形態による、３次元ＯＣＴインターフェログラムデータから２次元画像を生成するための機能ステップを示す図である。

【0020】

【図3】一実施形態による、ＯＣＴインターフェログラムデータを処理するためのワークフローを示すワークフロー図である。

【0021】

【図4】一実施形態による、３次元ＯＣＴインターフェログラムデータから２次元画像を生成するための流れを示す流れ図である。

【0022】

【図5】一実施形態による、２次元画像強度マップが利用される撮像セッションの流れを示す流れ図である。

【0023】

【図6】３次元ＯＣＴインターフェログラムデータから２次元画像を生成するのに使用されうる例示的コンピュータを示すハイレベルブロック図である。

【発明を実施するための形態】

【0024】

様々な実施形態によれば、フーリエ領域（例えば分光計や波長掃引光源型の）光断層撮影収束（ＯＣＴ）データが、ＯＣＴ撮像システムで入力として受け取られ、２次元画像強度マップが、逆累積分布関数に基づいてインターフェログラムデータから生成されうる。図１に、一実施形態による、３次元ＯＣＴインターフェログラムデータから２次元画像を生成するのに使用されうるシステムの図を示す。図１では、システム１００は、光ビームを生成するための光源１０２（広帯域光源など）を含む。ビームスプリッタ１０４が光源１０２からの光ビームをサンプル経路１０６と参照経路１０８とに沿って分割する。参照経路１０８は光源１０２からの参照光ビームを調整するための（例えば、最大干渉を達成するための）偏光コントローラ１１０と、レンズ１１４を通して反射鏡１１５に向けられる参照光ビームを平行にするためのコリメータ１１２とを含んでいてよく、空間的に調整可能でありうる。

【0025】

サンプル経路１０６は、光源１０２からの光ビームを、コリメータ１１７および１つもしくは複数の対物レンズ１１８を介して対象物１２０上のＸ−Ｙ平面内の異なる位置を照明するように向けるための２次元スキャナ１１６を含む。一実施形態では、対象物１２０は（図示のように）患者の眼とすることができ、システム１００は一般に、眼科画像を得ることを対象としうる。しかし、本明細書の様々な実施形態は、（それだけに限らないが）表皮、歯科、および脈管系の臨床撮像を含む様々な分野に適用されるため、対象物１２０は、対象となる様々な臨床分野（人間の組織、歯など）または非臨床分野のいずれを含んでいてもよい。したがって、本明細書のいくつかの例は眼底撮像および患者の眼である対象物１２０に言及しうるが、システム１００は眼科用途だけに限定されるものと解釈すべきではない。むしろ、様々な実施形態は様々な用途に関連したものとすることができ、本明細書に記載する技法は広範に適用することができる。

【0026】

光ビームがスキャナ１１６によって対象物１２０上の特定のＸ−Ｙ平面位置に向けられた結果として、インターフェログラム検出部１２２はサンプル経路１０６から返された光を受け取る。またインターフェログラム検出部１２２は、対象物１２０から測定された光と参照経路１０８から測定された光との間の干渉スペクトルから光のエコー時間遅延を測定するために、参照経路１０８から返される光も受け取る。一実施形態では、光源１０２は広帯域光ビームを生成してよく、インターフェログラム検出部１２２は、検出器アレイと、サンプル経路１０６から返される光および参照経路１０８から返される光を波長の関数として角度分散させるための回折格子とを含む分光計とすることができる。あるいは、光源１０２は可同調掃引光ビームを生成してもよく、光源の異なる波長の光強度が、インターフェログラム検出部１２２によって時間の関数として得られてもよい。

【0027】

一実施形態では、インターフェログラム検出部１２２は、サンプル経路１０６および参照経路１０８から受け取られる光（すなわちインターフェログラムデータ）に基づいて複数の出力セットを生成する。例えば、インターフェログラム検出部１２２によって生成される複数の出力セットの各々は、光源１０２の異なる波長で受け取られる光強度に対応しうる。したがって、光源１０２からの光ビームがスキャナ１１６によって特定のＸ−Ｙ平面位置に向けられるときに、検出される光強度は、特定のＸ−Ｙ平面位置における深さ方向Ｚの対象物１２０内の光反射分布に関する情報を含むことができる。したがって、各出力セットは３次元データセットを含む。

【0028】

処理装置１２４は、インターフェログラム検出部１２２によって生成された複数の出力セットを、コンピュータバスによって、あるいは無線で、あるいは、公衆もしくは専用のネットワークを介して（例えばインターネットやアクセス制限イントラネットワークを介して）など、様々な通信手段のいずれかによって受け取る。一実施形態では、処理装置１２４は、例えば、相互にリモートに位置しうる複数の処理装置を含みうる。

【0029】

図２に、一実施形態による、３次元ＯＣＴインターフェログラムデータから２次元画像を生成するための機能ステップの図を示す。例えば、ＯＣＴインターフェログラムデータ（インターフェログラム検出部１２２によって生成された複数の３次元出力セットなど）を使用して（２００）、処理装置１２４において、医学的診断用の眼底マップといった、２次元強度マップＯＣＴ画像を生成し（２０２）、表示する（２０４）ためのコンピュータプログラム命令が実行されうる。一実施形態では、生成画像は、例えば、生成画像を表示し、位置合わせし、セグメント化するための、対象物１２０の２次元画像強度マップ、またはそのセグメント化もしくは変換の結果を生成するのに適したものとすることができる。あるいは、生成画像は、ＯＣＴや眼底撮像カメラといった画像診断法と関連付けられた走査取り込み整合プロセスを実施するのに適したものであってもよい。加えて、または代替として、３次元断面画像を生成する（２０８）ためにＯＣＴインターフェログラムデータを変換する（２０６）コンピュータプログラム命令が実行されてもよい。

【0030】

フーリエ変換計算を利用して１つもしくは複数のフィルタ出力セットを生成することができる。一般的なフーリエ領域ＯＣＴインターフェログラム方程式が次式で与えられる。

【数1】

式中、ｖは光ビームの周波数であり、Ｒ_nおよびＲ_mは、それぞれ、対象物１２０からの、特定のＸ−Ｙ平面位置ｎおよびｍにおける強度反射であり、Ｇ_s（ｖ）は光源１０２のスペクトル密度であり、参照アームの反射は均一であり、距離は、サンプル経路１０６における伝搬時間τ_nおよびτ_mならびに参照経路１０８における伝搬時間τ_rによって表されている。括弧内の第３項は対象物１２０内で散乱するすべての光についての相互干渉であり、最後の項は、対象物１２０からの散乱光と参照経路１０８との間の干渉を含み、そこからＡ走査（すなわち軸方向長走査）が計算される。

【0031】

インターフェログラム検出部１２２が分光計である実施形態では、パラメータｖは、ライン・スキャン・カメラ上の位置も表す。したがって、最初の２項はカメラ画素の線形アレイの端から端までの緩慢な変動を表し、最後の２項は振動（すなわち干渉じま）を表しうる。同様に、光源１０２が可同調掃引光源であり、ある範囲の周波数ｖを掃引される実施形態では、１次方程式の最初の２項は検出される信号の緩慢な時間的変動を表し、最後の２項は振動を表す。

【0032】

一実施形態では、高域フィルタリングにより式の右辺の第３項および第４項を除くすべてが消去され、低反射率を有する対象物がサンプル経路１０６に沿って存在するときには第４項の大きさが通常は第３項よりもずっと大きいものと仮定して、事実上、第４項が分離されることになる。さらに、高域フィルタリングされたスペクトルの範囲、偏り、標準偏差、分散、またはエントロピーを使用して、対象物１２０の２次元画像強度マップが生成されてもよい。

【0033】

図３に、一実施形態による、ＯＣＴインターフェログラムデータを処理するためのワークフロー図を示す。図３に示す方法は、高域フィルタリングされたインターフェログラムデータの分布を推定し、定量化して、対象物１２０上の１つもしくは複数のＸ−Ｙ平面位置に対応する１つもしくは複数の単一の推定強度値を求めるステップを含む。例えば、３００で、インターフェログラムデータ（複数の出力セットなど）がインターフェログラム検出部１２２から受け取られる。インターフェログラムデータは、フーリエ領域インターフェログラム式の各項を簡約して前述のように第４項を事実上分離するために、３０２で高域フィルタリングされてよい。一実施形態では、高域フィルタリングされたインターフェログラムデータのデータ率またはサイズは、例えば、より高い処理効率やより高速な計算を可能にするために、３０４で、様々な公知のダウンサンプリングまたは切り捨ての技法によってさらに低減されてよい。あるいは、データは、インターフェログラムデータの高域フィルタリングの前に切り捨てられ、またはダウンサンプリングされてもよい。例えば、（高域フィルタリングの前または後の）インターフェログラムデータは、ｎ番目ごとの出力を保持し、他のすべてを廃棄することによってダウンサンプリングされてもよい。別の例では、インターフェログラムデータは、選択された範囲内のデータを保持し、他のすべてのデータを廃棄することによって切り捨てられてもよい。あるいは、１サンプルがｎサンプルごとの中からランダムに、または疑似ランダムに選択される疑似ランダム・ダウンサンプリング・プロセスが用いられてもよい。

【0034】

次いでフィルタ出力セットは、逆累積分布関数（ＣＤＦ）に基づいて、特定のＸ−Ｙ平面位置についての深さ方向Ｚの単一の推定強度値に変換されうる。一般に、ＣＤＦ関数は、所与の確率分布を有する実数値の確率変数ｘがｘ以下の値において見つかる確率を示す。確率変数の逆ＣＤＦは、分位関数とも呼ばれ、所与の確率について、確率変数がその確率でその値以下になる値を特定する。

【0035】

あるいは、逆ＣＤＦは、クイックセレクト（ｑｕｉｃｋ ｓｅｌｅｃｔ）関数またはホーア（Ｈｏａｒｅ）の選択関数とすることもでき、その場合一般に、全フィルタ出力セットをソートする必要はない。そうではなく、フィルタ出力セットの一部分だけ、例えば、所望の分位点選択を含むピボット値の前後のサブアレイなどが実際にソートされる。

【0036】

別の代替形態では、逆ＣＤＦまたは分位関数は、フィルタ出力セットのソートを必要とせずに分位点の値を推定するように動作するヒューリスティックセットを用いてもよい。

【0037】

最初の逆ＣＤＦの手法に戻って、フィルタ出力セットは、３０６で二乗され（すなわち、正の値および負の値が同じ範囲に存在するように）、３１０の逆ＣＤＦ関数のために３０８でソートされてよい。一実施形態では、フィルタ出力セットは、３０８で、逆ＣＤＦ計算のために選択されうる任意のパーセンタイルに対応する値を求めるために、昇順または降順にソートされてよい。例えば、０番目のパーセンタイル値は最小パーセンタイル値および最小フィルタ出力値に対応していてよく、５０番目のパーセンタイルは中央パーセンタイル値および中央フィルタ出力値に対応していてよく、１００番目のパーセンタイル値は最大パーセンタイル値および最大フィルタ出力値に対応していてよい。したがって、２０４８個のフィルタ出力があり、９０番目のパーセンタイルが選択される場合、１８４３番目のソート値が選択値である。一例示的実施形態では、少なくとも１つの選択パーセンタイル値は、５０番目のパーセンタイル値と１００番目のパーセンタイル値との間の値に対応しうる。しかし、少なくとも１つの選択パーセンタイル値は、フィルタ出力セット内の任意の事前選択パーセンタイル値に対応しうる。したがって、任意のパーセンタイルに対応するソート値が選択されると、ソート値は３１２で単一の推定強度値に変換されうる。

【0038】

最初の逆ＣＤＦの手法の一変形では、３０８で、フィルタ出力セットのソートを必要とせずに分位点の値を推定するように動作するヒューリスティックセットが用いられてよい。推定値は次いで３１２で単一の推定強度値に変換されうる。

【0039】

代替の例示的実施形態では、単一の推定強度値を求めるのに、逆ＣＤＦの手法と概念的に類似した、最大値および最小値に基づく手法が用いられうる。最大値および最小値に基づく手法では、フィルタ出力セットは、単一の推定強度値を生成するためにソートされる必要がない。例えば、最大パーセンタイルに基づく手法は一般に、１００番目のパーセンタイルが選択されるときに逆ＣＤＦの手法に対応する。一実施形態では、３１４で負の値の可能性を無くすためにフィルタ出力の二乗または絶対値が求められてよく、３１６でこれらの結果から求められた最大（１００番目の）パーセンタイル値が単一の推定強度値（３１２）を表しうる。最大値および最小値に基づく一手法では、３１８で最大パーセンタイル値および最小パーセンタイル値が求められ、３２０で負の値の可能性を無くすために最大パーセンタイル値および最小パーセンタイル値の二乗または絶対値が求められうる。次いで３２２で、結果として得られる値が組み合わされて単一の推定強度値（３１２）が求められうる。例えば、組み合わせ演算は、２つの結果として得られる値の最小値、最大値、または平均値を取ることを含んでいてよい。また、単一の推定強度値（３１２）の大きさは、例えば、生成値の対数または根（二乗根、三乗根、四乗根など）を取ることによってさらに変換されてもよい。さらに、最小値および最大値に基づく一手法は、ブロック３１８で上位ｎ個の最小値および上位ｎ個の最大値が求められるように一般化されてもよい。この一般化された実施形態では、３２２の組み合わせ演算は、２ｎ個の値の平均値を取ること、または逆ＣＤＦ法を使用して事前選択パーセンタイル値を取ることからなっていてもよいはずである。１つもしくは複数のＸ−Ｙ平面位置（すなわち複数のＡライン）に対応する１つもしくは複数の単一の推定強度値が、対象物１２０の２次元画像強度マップを生成するのに適しうる。

【0040】

図４に、一実施形態による、３次元ＯＣＴインターフェログラムデータから２次元画像を生成するための流れ図を示す。図１のシステムを例として使用すると、４０２で、光源１０２からの光ビームが生成され、サンプル経路１０６および参照経路１０８に沿って分割される。例えば、光源１０２は広帯域光ビームを生成してよく、インターフェログラム検出部１２２は、回折格子と検出器アレイとを含む分光計とすることができる。あるいは、光源１０２は、可同調掃引光ビームを生成してもよく、光源の異なる波長の光強度が、インターフェログラム検出部１２２によって時間の関数として得られてもよい。

【0041】

４０４で、光ビームは、スキャナ１１６により、サンプル経路１０６に沿って、対象物１２０上のＸ−Ｙ平面内の異なる位置に向けられる。４０６で、サンプル経路１０６および参照経路１０８の各々から返された光が、インターフェログラム検出部１２２で受け取られる。４０８でインターフェログラム検出部１２２は複数の出力セットを生成し、複数の出力セットの各々は、光ビームが特定のＸ−Ｙ平面位置に向けられるときに光源の異なる波長で受け取られる光強度に対応する。例えば、光強度は、特定のＸ−Ｙ平面位置における深さ方向Ｚの対象物１２０内の光反射率分布に関する情報を含みうる。

【0042】

４１０で、光ビームを特定のＸ−Ｙ平面位置に向けることにより生成される出力セットは、任意選択で、フィルタ出力セットを生成するために高域フィルタリングされてもよい。４１２で、フィルタリングされた、またはフィルタリングされていない出力のセットは、上記で図３に関連して説明したように、逆累積分布関数に基づいて、特定のＸ−Ｙ平面位置についての深さ方向Ｚの単一の推定強度値に変換される。４１４で、対象物の２次元画像強度マップが、１つもしくは複数のＸ−Ｙ平面位置に対応する１つもしくは複数の単一の推定強度値から生成される。２次元画像強度マップは、ディスプレイにおける提示に適したものとすることができる。一実施形態では、単一の推定強度値は、画像強度マップのリアルタイム表示のために変換されうる。

【0043】

図５に、一実施形態による、２次元画像強度マップが利用される撮像セッションのための流れ図を示す。対象物１２０の２次元画像強度マップは、画像診断法と関連付けられた走査取り込み整合プロセスに使用されうる。例えば、画像診断法は、ＯＣＴまたは別のシステム（眼科のための眼底撮像カメラなど）のうちの１つとすることができる。さらに、走査取り込み整合プロセスは自動的に行われてもよく、画像診断法は、例えば技術者などによって、手作業で整合されてもよい。５０２で、被験者（患者など）がＯＣＴ（眼底など）画像のために位置決めされ、５０４で画像診断法が固定される（例えば、画像診断法のシャシが特定の位置で固定して設定される）。５０６で、必要に応じて、画像診断法の焦点が合わせられる。

【0044】

上記図４で説明したように、５０８で３次元ＯＣＴインターフェログラムが取り込まれてよく、５１０で２次元画像強度マップが自動的に生成され、リアルタイムで表示されてよい。一実施形態では、対象物１２０のＯＣＴ画像が、それだけに限らないが、逆累積分布関数を含む方法を含む様々な手段のいずれかによって生成されうる。５１２でＯＣＴ走査または眼底画像取り込みが開始される。一実施形態では、画像診断法が手作業で整合されてよく、５０６で（例えば、別の画像を取り込むのに必要な場合などに）画像診断法の焦点を再度合わせる前に、５１４で固定状態が調整されてよい。

【0045】

様々な実施形態において、走査取り込み整合プロセスは、非散瞳状態においてさえも、正常な撮像セッションを可能にするように機能しうる。また対象物１２０の２次元画像強度マップは、別の画像診断法によって得られた対象物１２０の画像を位置合わせするのに使用されてもよい。

【0046】

加えて、（フィルタリングされていない）出力のセットは、特定のＸ−Ｙ平面位置についての深さ方向Ｚの強度情報を有する３次元データセットに変換されてもよく、その場合、１つもしくは複数の３次元データセットが対象物１２０のＯＣＴ断面画像を生成するのに使用される。例えば、３次元データセットは、１つもしくは複数のランドマーク、例えば対象物１２０の１つもしくは複数の物理境界などに基づいてセグメント化されてよく、ランドマークは深さ方向Ｚの強度情報から識別されうる。様々な実施形態において、セグメント化３次元データセットは、対象物１２０の部分強度画像もしくは眼球マップを生成し、または対象物１２０の液体充満空間を定量化するのに利用されうる。

【0047】

さらに別の実施形態において、対象物１２０の２次元画像強度マップは、例えば、対象物１２０の２次元画像強度マップを対象物１２０のＯＣＴ断面画像と並列に表示することによって、対象物１２０のＯＣＴ断面画像を位置合わせするのに使用されうる。

【0048】

本明細書に記載するシステム、装置、および方法は、ディジタル回路を使用して、または周知のコンピュータプロセッサ、メモリ装置、記憶装置、コンピュータソフトウェア、および他の構成要素を使用した１台もしくは複数のコンピュータを使用して実施されうる。通常、コンピュータは、命令を実行するためのプロセッサと、命令およびデータを記憶するための１もしくはそれ以上のメモリとを含む。またコンピュータは、１つもしくは複数の磁気ディスク、内蔵ハードディスクおよび取り外し可能ディスク、光磁気ディスク、光ディスクなどといった１つもしくは複数の大容量記憶装置も含んでいてよく、またはそれらに結合されていてもよい。

【0049】

本明細書に記載するシステム、装置、および方法は、クライアント−サーバ関係で動作するコンピュータを使用して実施されうる。通常、そのようなシステムでは、クライアントコンピュータはサーバコンピュータから離れて位置し、ネットワークを介して対話する。クライアント−サーバ関係は、それぞれのクライアントコンピュータおよびサーバコンピュータ上で走るコンピュータプログラムによって定義され、制御されうる。

【0050】

本明細書に記載するシステム、装置、および方法は、ネットワークベースのクラウド・コンピューティング・システム内で使用されうる。そのようなネットワークベースのクラウド・コンピューティング・システムでは、ネットワークに接続されているサーバまたは別のプロセッサが、ネットワークを介して１つもしくは複数のクライアントコンピュータと通信する。クライアントコンピュータは、例えば、クライアントコンピュータ上に常駐し、そこで動作するネットワーク・ブラウザ・アプリケーションによってサーバと通信してよい。クライアントコンピュータは、サーバ上にデータを記憶し、ネットワークを介してデータにアクセスしてよい。クライアントコンピュータは、ネットワークを介してサーバに、データを求める要求、またはオンラインサービスを求める要求を送ってよい。サーバは、要求されたサービスを実行し、（１台もしくは複数の）クライアントコンピュータにデータを提供してよい。またサーバは、クライアントコンピュータに、例えば、計算を行わせる、画面上に指定のデータを表示させるなど、指定の機能を果たさせるように適合されたデータを送ってもよい。例えば、サーバは、クライアントコンピュータに、図４および図５のステップのうちの１つもしくは複数を含む本明細書に記載した方法ステップのうちの１つもしくは複数を実行させるように適合された要求を送ってもよい。図４および図５のステップのうちの１つもしくは複数を含む本明細書に記載した方法ステップのうちのある特定のステップが、ネットワークベースのクラウド・コンピューティング・システム内のサーバまたは別のプロセッサによって実行されてもよい。図４および図５のステップのうちの１つもしくは複数を含む本明細書に記載した方法ステップのうちのある特定のステップが、ネットワークベースのクラウド・コンピューティング・システム内のクライアントコンピュータによって実行されてもよい。図４および図５のステップのうちの１つもしくは複数を含む本明細書に記載した方法ステップのうちのステップは、ネットワークベースのクラウド・コンピューティング・システム内のサーバによって、かつ／またはクライアントコンピュータによって、任意の組み合わせで実行されてよい。

【0051】

本明細書に記載するシステム、装置、および方法は、プログラマブルプロセッサが実行するための、情報担体において、例えば非一時的な機械可読記憶装置において有形的に具現化されたコンピュータプログラム製品を使用して実施されてよく、図４および図５のステップのうちの１つもしくは複数を含む本明細書に記載する方法ステップは、そのようなプロセッサが実行することのできる１つもしくは複数のコンピュータプログラムを使用して実施されてよい。コンピュータプログラムは、コンピュータにおいて、ある特定の活動を行い、またはある特定の結果を生じさせるために、直接、または間接的に使用することができるコンピュータプログラム命令のセットである。コンピュータプログラムは、コンパイラ型言語またはインタープリタ型言語を含む、任意の形のプログラミング言語で書くことができ、独立型プログラム、またはモジュール、コンポーネント、サブルーチン、もしくはコンピューティング環境での使用に適した他のユニットを含む任意の形で配備することができる。

【0052】

本明細書に記載するシステム、装置および方法を実施するのに使用されうる例示的コンピュータのハイレベルブロック図が図６に示されている。コンピュータ６００は、データ記憶装置６２０とメモリ６３０とに動作可能に結合されたプロセッサ６１０を備える。プロセッサ６１０は、コンピュータ６００の全般的動作を定義するコンピュータプログラム命令を実行することによってそのような動作を制御する。コンピュータプログラム命令は、データ記憶装置６２０に、もしくは他のコンピュータ可読媒体に記憶されており、コンピュータプログラム命令の実行が求められるときにメモリ６３０にロードされてよい。例えば、図１を参照すると、処理装置１２４は、コンピュータ６００の１つもしくは複数の構成要素を備えていてよい。よって、図４および図５の方法ステップは、メモリ６３０および／またはデータ記憶装置６２０に記憶されたコンピュータプログラム命令によって定義され、コンピュータプログラム命令を実行するプロセッサ６１０によって実行されうる。例えば、コンピュータプログラム命令は、図４および図５の方法ステップによって定義されるアルゴリズムを実行するように当業者によってプログラムされたコンピュータ実行可能コードとして実装することができる。したがって、コンピュータプログラム命令を実行することによって、プロセッサ６１０は、図４および図５の方法ステップによって定義されるアルゴリズムを実行する。またコンピュータ６００は、ネットワークを介して他の装置と通信するための１つもしくは複数のネットワークインターフェース６４０も含む。またコンピュータ６００は、コンピュータ６００とのユーザ対話を可能にする１つもしくは複数の入力／出力装置６５０（ディスプレイ、キーボード、マウス、スピーカ、ボタンなど）も含む。

【0053】

プロセッサ６１０は、汎用マイクロプロセッサと専用マイクロプロセッサの両方を含んでいてよく、コンピュータ６００の唯一のプロセッサであっても、その複数のプロセッサのうちの１つであってもよい。プロセッサ６１０は、例えば、１つもしくは複数の中央処理装置（ＣＰＵ）を含んでいてよい。プロセッサ６１０、データ記憶装置６２０、および／またはメモリ６３０は、１つもしくは複数の特定用途向け集積回路（ＡＳＩＣ）、１つもしくは複数のフィールド・プログラマブル・ゲート・アレイ（ＦＰＧＡ）、および／または１つもしくは複数のディジタル信号プロセッサ（ＤＳＰ）ユニットを含み、またはそれらによって補足され、またはそれらに組み込まれていてよい。

【0054】

データ記憶装置６２０およびメモリ６３０は、各々、有形の非一時的なコンピュータ可読記憶媒体を備える。データ記憶装置６２０およびメモリ６３０は、各々、ダイナミック・ランダム・アクセス・メモリ（ＤＲＡＭ）、スタティック・ランダム・アクセス・メモリ（ＳＲＡＭ）、ダブル・データ・レート・シンクロナス・ダイナミック・ランダム・アクセス・メモリ（ＤＤＲ ＲＡＭ）、他のランダムアクセス固体メモリデバイスといった高速ランダム・アクセス・メモリを含んでいてよく、内蔵ハードディスクや取り外し可能ディスクといった１つもしくは複数の磁気ディスク記憶装置、光磁気ディスク記憶装置、光ディスク記憶装置、フラッシュ・メモリ・デバイス、消去可能プログラム可能読取り専用メモリ（ＥＰＲＯＭ）や電気的消去可能プログラム可能読取り専用メモリ（ＥＥＰＲＯＭ）といった半導体メモリデバイス、コンパクトディスク読取り専用メモリ（ＣＤ−ＲＯＭ）ディスク、ディジタル多用途ディスク読取り専用メモリ（ＤＶＤ−ＲＯＭ）ディスク、あるいは他の不揮発性固体記憶装置といった不揮発性メモリを含んでいてよい。

【0055】

入力／出力装置６５０は、プリンタ、スキャナ、表示画面といった周辺装置を含んでいてよい。例えば、入力／出力装置６５０は、ユーザに情報を表示するための陰極線管（ＣＲＴ）、プラズマもしくは液晶ディスプレイ（ＬＣＤ）モニタといった表示装置、キーボード、およびユーザがコンピュータ６００に入力を提供するためのマウスやトラックボールといったポインティングデバイスを含んでいてよい。

【0056】

処理装置１２４およびインターフェログラム検出部１２２を含む本明細書で論じたシステムおよび装置のいずれかまたは全部が、コンピュータ６００といったコンピュータを使用して実施されてよい。

【0057】

実際のコンピュータまたはコンピュータシステムの実施態様は、他の構造を有していてよく、他の構成要素も含んでいてよいこと、および図６は例示のためのそのようなコンピュータの構成要素のうちのいくつかのハイレベル表現であることを当分野の技術者は理解するであろう。

【0058】

以上の詳細な説明は、あらゆる点において、限定のためではなく、例証および例示のためのものであり、本願で開示する発明の範囲は、詳細な説明から決定されるべきではなく、特許法によって許容される最大限の広さに従って解釈される特許請求の範囲から決定されるべきものである。本明細書で図示し、説明した実施形態は、本発明の原理を例示するためのものにすぎず、様々な改変が本発明の範囲および趣旨を逸脱することなく当業者によって実施されうることを理解すべきである。当業者は、本発明の範囲および趣旨を逸脱することなく、様々な他の特徴の組み合わせを実施することもできるはずである。

【符号の説明】

【0059】

１００ システム
１０２ 光源
１０４ ビームスプリッタ
１０６ サンプル経路
１０８ 参照経路
１１０ 偏光コントローラ
１１２ コリメータ
１１４ レンズ
１１５ 反射鏡
１１６ ２次元スキャナ
１１７ コリメータ
１１８ 対物レンズ
１２０ 対象物
１２２ インターフェログラム検出部
１２４ 処理装置
２００ ＯＣＴインターフェログラムデータ
６００ コンピュータ
６１０ プロセッサ
６２０ データ記憶装置
６３０ メモリ
６４０ ネットワークインターフェース
６５０ 入力／出力装置

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】