JP7109644B2 - 体積ｏｃｔ画像データ処理 - Google Patents

Description

本発明の例示的態様は、一般に眼科用光干渉断層計（ＯＣＴ）撮像システムに関する。より具体的には、光干渉断層計（ＯＣＴ）撮像システムにより取得された撮像ターゲットのＢスキャンのシーケンスから成るＣスキャンデータを処理して、ＯＣＴ撮像システムによるＢスキャン取得時のＯＣＴ撮像システムと撮像ターゲットの相対運動による両者の距離変動に起因するＢスキャンシーケンス内のＢスキャン間の軸方向変位を補償するための補正データ生成の方法及び装置に関する。

患者の眼の一部、例えば網膜の典型的な体積光干渉断層計（ＯＣＴ）スキャン（Ｃスキャンとも称する）には、約１～５秒かかる。この間に、ＯＣＴ撮像システムは、第１走査方向にＯＣＴビームを反復走査し、（２次元の）Ｂスキャンのシーケンスを記録する。そのそれぞれは網膜表面上の第１走査方向に沿う個々の点で記録される一続きの軸方向走査から成る。ＢスキャンシーケンスにおけるＢスキャンは、通常第１走査に垂直な方向に配列される。Ｃスキャンの取得中に、通常、患者の不随意運動によって、眼が軸方向（ＯＣＴビームの方向）に移動する場合がある。Ｃスキャン画像をうまく画像化し、かつそのＣスキャン画像を基に行われる後続の計測の精度を確保するために、Ｃスキャン撮像中の患者の動きによる運動アーチファクトを補償することが必要となる場合がある。患者の動きによる、いくつかのＢスキャン間での軸方向シフトを補正するためのＢスキャンデータへの補償が行われない場合には、網膜層識別の精度及び後続する診断計測の精度に悪影響を与える可能性がある。

本発明者らは、本明細書の第１の例示的態様にしたがって、ＯＣＴ撮像システムにより取得された撮像ターゲットのＢスキャンのシーケンスからなるＣスキャンデータを処理して、ＯＣＴ撮像システムによるＢスキャン取得中にその間の距離を変動させる、ＯＣＴ撮像システムと撮像ターゲットとの相対運動により生じる、Ｂスキャンシーケンス内のＢスキャン同士の間の軸方向変位を補償するための補正データを生成し、またさらにその生成された補正データの信頼性を示す信頼性指標を生成する方法を考案した。この方法は、シーケンス内の複数の隣接Ｂスキャンペアの内の各隣接Ｂスキャンペアに対し、隣接するＢスキャン内の共通の視覚特徴のそれぞれの表示の間の軸方向シフトのそれぞれの指標を決定することによって補正データを生成することを含む。この方法はさらに、Ｂスキャンシーケンス内のＢスキャンの複数ペアを使用して、そのＢスキャンの複数ペアが取得されたときのＯＣＴ撮像システムに対する撮像ターゲットの相対的なそれぞれの速度又はそれぞれの加速度を示すメトリックのそれぞれの値を計算し、そのメトリックの少なくとも所定数の計算値が閾値を超えるか否かを判定し、メトリックの少なくとも所定数の計算値が閾値を超えると判定される場合に、信頼性指標を補正データが信頼できないことを示すように設定し、メトリックの少なくとも所定数の計算値が閾値を超えないと判定される場合に、信頼性指標を補正データが信頼できることを示すように設定する、ことによって、信頼性指標を生成することを含む。

この補正データは、シーケンス内の隣接するＢスキャンの対応するペアの位置による決定された指標の変動から、ＯＣＴ撮像システムによるＢスキャンの取得中におけるＯＣＴ撮像システムと撮像ターゲットの相対運動を示す、変動の第１周波数成分を決定することにより生成され得る。

軸方向シフトのそれぞれの指標は、各隣接Ｂスキャンペアに対して、隣接するＢスキャンペア間の相互相関を計算し、計算された相互相関におけるピークに対応するＢスキャンの間のオフセットを指標として決定するか、あるいは、隣接するＢスキャンペアのＢスキャンにおける共通の視覚特徴のそれぞれの位置を特定し、撮像システムの軸方向を表す、Ｂスキャン軸に沿う特定された位置同士の間の変位を決定することにより、決定され得る。

撮像ターゲットは曲率を有する場合があり、本方法はさらに、撮像ターゲットの曲率を示す、変動の第２周波数成分を決定することを含み得る。第２周波数成分は、決定された指標の変動に対してｍ次多項式をフィッティングさせることにより決定され得る。第１周波数成分は、決定された指標の変動内の指標から、ｍ次多項式の値を引くことにより補正された指標変動を生成し、補正された指標変動にｎ次の多項式をフィッティングさせることにより決定され得る。ここで、ｍとｎは整数であって、ｍがｎより小さい。

本方法は、補正データが信頼できることを示す信頼性指標が設定されている場合には、補正データを使用して、ＯＣＴ撮像システムと撮像ターゲットの相対運動により生じるＢスキャンシーケンス内のＢスキャン間の軸方向変位を補償することを更に含み得る。

第１の例示的態様の実施形態において、第１周波数成分は、（ｉ）各残差値が複数の指標の変動における指標と、ｎ次多項式の対応する値との間の差として計算される、複数の残差値を計算するステップと、（ｉｉ）その複数の残差値に、第１の正の閾値より高いか第１の負の閾値より低い外れ値が含まれるか否かを判定するステップと、（ｉｉｉ）その複数の残差値に外れ値が含まれると判定される場合、外れ値に対応する指標を指標の変動から除去して、更新された指標変動を生成し、また、残差が外れ値を含まないと判定される場合には、ｎ次多項式を第１周波数成分として決定して処理を終了するステップと、（ｉｖ）ｎ次多項式を更新された指標変動にフィッティングさせるステップと、から成る処理の少なくとも２回の反復を更に実行することにより決定され得る。ここで、複数の残差値における各残差値は、処理の第１の反復において、補正された指標変動における指標と、指標の補正された変動にフィッティングされたｎ次多項式の対応する値との間の差として計算され、複数の残差値における各残差値は、処理の残りの１以上の反復の各反復において、処理のそれ以前の反復において生成された、更新された指標変動における指標と、処理のそれ以前の反復において生成された、更新された指標変動にフィッティングされたｎ次多項式の対応する値との間の差として計算される。

第１の例示的態様の実施形態において本方法は、信頼性指標が補正データが信頼できることを示すように設定されている場合において、第２の正の閾値が第１の正の閾値より小さく第２の負の閾値が第１の負の閾値より大きい場合、複数の残差値の中で第２の正の閾値よりも大きいか第２の負の閾値より小さい大きさを有する残差値の数を決定することと、決定された残差値の数が第３の閾値より小さい場合に、Ｂスキャンシーケンス内のＢスキャンに第１周波数成分に基づくオフセットを適用することによってＢスキャンシーケンス内のＢスキャン間の軸方向変位を補償することと、決定された残差値の数が第３の閾値より小さくない場合に、Ｂスキャンシーケンス内のＢスキャン間の軸方向変位を補償しないように決定すること、とをさらに含む。

さらに本発明者らは、本明細書に記載の第２の例示的態様に従って、撮像ターゲットを横断して延在する、平行な複数の走査線に沿って走査するＯＣＴ撮像システムにより取得された撮像ターゲットのＢスキャンシーケンスからなるＣスキャンデータと、走査線を横切るトレース走査線に沿って走査を実行するＯＣＴ撮像システムによって取得されたトレーススキャンデータとを処理して、ＯＣＴ撮像システムによるＢスキャンの取得中にその間の距離を変動させるＯＣＴ撮像システムと撮像ターゲットとの相対運動により生じるＢスキャンシーケンス内のＢスキャン同士の間の軸方向変位を補償するための補正データを生成し、かつ、さらに生成された補正データの信頼性を示す信頼性指標を生成する方法を考案した。本方法は、シーケンスにおける少なくともいくつかのＢスキャンの各Ｂスキャンに対し、Ｂスキャン及びトレーススキャンデータ内の共通の視覚特徴のそれぞれの表示の間の軸方向シフトのそれぞれの指標を決定することにより補正データを生成することを含む。この方法はさらに、少なくともいくつかのＢスキャン内のＢスキャンのペアに対して決定された指標を使用して、Ｂスキャンの複数ペアが取得されたときのＯＣＴ撮像システムに対する撮像ターゲットの相対的なそれぞれの速度又はそれぞれの加速度を表すメトリックのそれぞれの値を計算し、そのメトリックの少なくとも所定数の計算値が閾値を超えるか否かを判定し、メトリックの少なくとも所定数の計算値が閾値を超えると判定される場合に、補正データが信頼できないことを示す信頼性指標を設定し、メトリックの少なくとも所定数の計算値が閾値を超えないと判定される場合に、補正データが信頼できることを示す信頼性指標を設定する、ことによって、信頼性指標を生成することを含む。

第２の例示的態様の方法は、補正データが信頼できることを示す信頼性指標が設定されている場合には、補正データを使用して、ＯＣＴ撮像システムと撮像ターゲットの相対運動により生じる、Ｂスキャンシーケンス内のＢスキャン間の軸方向変位を補償することを更に含み得る。追加又は代替として、第２の例示的態様の実施形態では、補正データが、（ｉ）各残差値が複数の指標の変動における指標と、ｎ次多項式の対応する値との間の差として計算される、複数の残差値を計算するステップと、（ｉｉ）その複数の残差値に、第１の正の閾値より高いか第１の負の閾値より低い外れ値が含まれるか否かを判定するステップと、（ｉｉｉ）複数の残差値に外れ値が含まれると判定される場合、外れ値に対応する指標を指標の変動から除去して更新された指標変動を生成し、また、残差が外れ値を含まないと判定される場合には、ｎ次多項式を補正データとして決定して処理を終了するステップと、（ｉｖ）ｎ次多項式を更新された指標変動にフィッティングさせるステップと、から成る処理の少なくとも２回の反復を更に実行することにより生成され得る。ここで、複数の残差値における各残差値は、処理の第１の反復において補正された指標変動における指標と指標の補正された変動にフィッティングされたｎ次多項式の対応する値との間の差として計算され、複数の残差値における各残差値は、処理の残りの１以上の反復の各反復において処理のそれ以前の反復において生成された更新された指標変動における指標と、処理のそれ以前の反復において生成された更新された指標変動にフィッティングされたｎ次多項式の対応する値との間の差として計算される。

第２の例示的態様の実施形態では、本方法は、信頼性指標が前記補正データが信頼できることを示すように設定されている場合において、第２の正の閾値が第１の正の閾値より小さく、第２の負の閾値が第１の負の閾値より大きい場合、複数の残差値の中で第２の正の閾値よりも大きいか第２の負の閾値より小さい大きさを有する残差値の数を決定することと、決定された残差値の数が第３の閾値より小さい場合に、Ｂスキャンシーケンス内のＢスキャンに補正データに基づくオフセットを適用することによってＢスキャンシーケンス内のＢスキャン間の軸方向変位を補償することと、決定された残差値の数が第３の閾値より小さくない場合に、Ｂスキャンシーケンス内のＢスキャン間の軸方向変位を補償しないように決定することと、をさらに含む。

さらに本発明者らは、本明細書に記載の第３の例示的態様に従って、プロセッサによって実行されるとき、上記の方法の少なくとも１つをプロセッサに実行させる、コンピュータプログラム命令を含む、コンピュータプログラムを考案した。コンピュータプログラムは、非一時的なコンピュータ可読記憶媒体に格納されてもよいし、信号によって搬送されてもよい。

さらに本発明者らは、本明細書に記載の第４の例示的態様に従って、ＯＣＴ撮像システムにより取得された撮像ターゲットのＢスキャンシーケンスからなるＣスキャンデータを、
ＯＣＴ撮像システムによりＢスキャンを取得中にその間の距離を変動させるＯＣＴ撮像システムと撮像ターゲットとの相対運動により生じるＢスキャンシーケンス中のＢスキャン同士の間の軸方向変位を補償するための補正データを生成するために、またさらに生成された補正データの信頼性を示す信頼性指標を生成するために、処理するように構成されたデータ処理装置を考案した。この装置は、シーケンス内の複数の隣接Ｂスキャンペアの内の各隣接Ｂスキャンペアに対し、隣接するＢスキャン内の共通の視覚特徴の各表示間の軸方向シフトのそれぞれの指標を決定することによって補正データを決定するように構成された補正データ生成モジュールを備える。この装置はさらに、Ｂスキャンシーケンス内のＢスキャンの複数ペアを使用して、そのＢスキャンの複数ペアが取得されたときのＯＣＴ撮像システムに対する撮像ターゲットの相対的なそれぞれの速度又はそれぞれの加速度を示すメトリックのそれぞれの値を計算し、そのメトリックの少なくとも所定数の計算値が閾値を超えるか否かを判定し、メトリックの少なくとも所定数の計算値が閾値を超えると判定される場合に、信頼性指標を補正データが信頼できないことを示すように設定し、メトリックの少なくとも所定数の計算値が閾値を超えないと判定される場合に、信頼性指標を補正データが信頼できることを示すように設定することによって、信頼性指標を生成するように構成された信頼性指標生成モジュールを更に備える。

補正データ生成モジュールは、シーケンス内の隣接するＢスキャンの対応するペアの位置による決定された指標の変動から、ＯＣＴ撮像システムによるＢスキャンの取得中におけるＯＣＴ撮像システムと撮像ターゲットの相対運動を示す変動の第１周波数成分をさらに決定することにより、補正データを生成するように構成することができる。

補正データ生成モジュールは、隣接Ｂスキャンのペア間の相互相関を計算し、かつ計算された相互相関のピークに対応するＢスキャン間のオフセットを指標として決定することにより、又は、隣接ＢスキャンのペアのＢスキャンにある共通の視覚特徴のそれぞれの位置を特定し、かつ特定された位置間の、撮像システムの軸方向を表すＢスキャンの軸に沿う変位を決定することにより、隣接Ｂスキャンの各ペアに対する軸方向シフトのそれぞれの指標を決定するように構成することができる。

撮像ターゲットは曲率を有する場合があり、補正データ生成モジュールはさらに撮像ターゲットの曲率を表す、変動の第２周波数成分を決定するように構成することができる。補正データ生成モジュールは、決定された指標の変動にｍ次多項式をフィッティングして第２周波数成分を決定し、その決定された指標の変動の指標からｍ次多項式の値を引き算して、補正された指標変動を生成し、その補正された指標変動に対してｎ次多項式をフィッティングして第１周波数成分を決定するように構成することができる。ここで、ｍとｎは整数であって、ｍがｎより小さい。

上記のように設定されたデータ処理装置はさらに変位補償モジュールを備えることができる。これは、信頼性指標が、補正データが信頼できることを示すように設定された場合に、その補正データを使用して、ＯＣＴ撮像システムと撮像ターゲットとの間の相対運動により生じるＢスキャンシーケンス内のＢスキャン間の軸方向変位を補償するように構成される。

第４の例示的態様の実施形態において補正データ生成モジュールは、（ｉ）各残差値が指標変動における指標とｎ次多項式の対応する値との差として計算される、複数の残差値を計算するステップと、（ｉｉ）その複数の残差値に、第１の正の閾値より高いか第１の負の閾値より低い外れ値が含まれるか否かを判定するステップと、（ｉｉｉ）その複数の残差値に外れ値が含まれると判定される場合、外れ値に対応する指標を指標の変動から除去して、更新された指標変動を生成し、また、残差が外れ値を含まないと判定される場合には、ｎ次多項式を第１周波数成分として決定して処理を終了するステップと、（ｉｖ）ｎ次多項式を更新された指標変動にフィッティングさせるステップと、から成る処理の少なくとも２回の反復を更に実行することにより、第１周波数成分を決定するように構成可能である。ここで、複数の残差値における各残差値は、処理の第１の反復において、補正された指標変動における指標と、その指標の補正された変動にフィッティングされたｎ次多項式の対応する値との間の差として計算され、複数の残差値における各残差値は、処理の残りの１以上の反復の各反復において、処理のそれ以前の反復において生成された更新された指標変動における指標と、処理のそれ以前の反復において生成された更新された指標変動にフィッティングされたｎ次多項式の対応する値との間の差として計算される。

第４の例示的態様の実施形態では、信頼性指標が補正データは信頼できることを示すように設定されている場合、補正データ生成モジュール（２）が、複数の残差値の中で、第２の正の閾値より大きいか第２の負の閾値より小さい大きさの残差値の数を決定するように更に構成可能である。ここで第２の正の閾値は第１の正の閾値より小さく、第２の負の閾値は第１の負の閾値より大きい。決定された残差値の数が第３の閾値より小さい場合には、変位補償モジュールが、Ｂスキャンシーケンス内のＢスキャンに対する補正データに基づいてオフセットを適用することによって、Ｂスキャンシーケンス内のＢスキャン間の軸方向変位に対する補償を行うように構成可能である。そして、決定された残差値の数が、第３の閾値より小さくない場合には、変位補償モジュールは、Ｂスキャンシーケンス内のＢスキャン間の軸方向変位の補償を行わないように構成可能である。

さらに本発明者らは、本明細書に記載の第５の例示的態様に従って、撮像ターゲットを横断して延在する平行な複数の走査線に沿って走査するＯＣＴ撮像システムにより取得された撮像ターゲットのＢスキャンシーケンスからなるＣスキャンデータと、走査線を横切るトレース走査線に沿って走査を実行するＯＣＴ撮像システムによって取得されたトレーススキャンデータとを処理して、ＯＣＴ撮像システムによるＢスキャンの取得中にその間の距離を変動させるＯＣＴ撮像システムと撮像ターゲットとの相対運動により生じる、Ｂスキャンシーケンス内のＢスキャン同士の間の軸方向変位を補償するための補正データを生成し、かつ、さらに生成された補正データの信頼性を示す信頼性指標を生成するように構成されたデータ処理装置を考案した。この装置は、シーケンスにおける少なくともいくつかのＢスキャンの各Ｂスキャンに対し、Ｂスキャンとトレーススキャンデータ内の共通の視覚特徴のそれぞれの表示の間の軸方向シフトのそれぞれの指標を決定することにより補正データを生成するように構成された補正データ生成モジュールを備える。この装置はさらに、少なくともいくつかのＢスキャン内のＢスキャンのペアに対して決定された指標を使用して、Ｂスキャンの複数ペアが取得されたときのＯＣＴ撮像システムに対する撮像ターゲットの相対的なそれぞれの速度又はそれぞれの加速度を表すメトリックのそれぞれの値を計算し、そのメトリックの少なくとも所定数の計算値が、閾値を超えるか否かを判定し、メトリックの少なくとも所定数の計算値が閾値を超えると判定される場合に、信頼性指標を補正データが信頼できないことを示すように設定し、メトリックの少なくとも所定数の計算値が閾値を超えないと判定される場合に、信頼性指標を補正データが信頼できることを示すように設定することによって、信頼性指標を生成するように構成された信頼性指標生成モジュールを更に備える。

第５の例示的態様のデータ処理装置はさらに変位補償モジュールを備えることができる。これは信頼性指標が補正データが信頼できることを示すように設定された場合に、補正データを使用して、ＯＣＴ撮像システムと撮像ターゲットとの間の相対運動により生じるＢスキャンシーケンス内のＢスキャン間の軸方向変位を補償するように構成される。

第５の例示的態様の実施形態において、正データ生成モジュールは、（ｉ）各残差値が指標変動における指標とｎ次多項式の対応する値との差として計算される複数の残差値を計算するステップと、（ｉｉ）その複数の残差値に、第１の正の閾値より高いか第１の負の閾値より低い外れ値が含まれるか否かを判定するステップと、（ｉｉｉ）その複数の残差値に外れ値が含まれると判定される場合、外れ値に対応する指標を指標の変動から除去して更新された指標変動を生成し、また、残差が外れ値を含まないと判定される場合には、ｎ次多項式を補正データとして決定して処理を終了するステップと、（ｉｖ）ｎ次多項式を更新された指標変動にフィッティングさせるステップと、から成る処理の少なくとも２回の反復を実行することにより、補正データを生成するように構成可能である。ここで、複数の残差値における各残差値は、処理の第１の反復において、補正された指標変動における指標と、その指標の補正された変動にフィッティングされたｎ次多項式の対応する値との間の差として計算され、複数の残差値における各残差値は、処理の残りの１以上の反復の各反復において、処理のそれ以前の反復において生成された更新された指標変動における指標と、処理のそれ以前の反復において生成された更新された指標変動にフィッティングされたｎ次多項式の対応する値との間の差として計算される。

第５の例示的態様の実施形態では、信頼性指標が補正データは信頼できることを示すように設定されている場合、補正データ生成モジュールが、複数の残差値の中で、第２の正の閾値より大きいか第２の負の閾値より小さい大きさの残差値の数を決定するように構成可能である。ここで第２の正の閾値は第１の正の閾値より小さく、第２の負の閾値は第１の負の閾値より大きい。決定された残差値の数が第３の閾値より小さい場合には、変位補償モジュールが、Ｂスキャンシーケンス内のＢスキャンに対する補正データに基づいてオフセットを適用することによって、Ｂスキャンシーケンス内のＢスキャン間の軸方向変位に対する補償を行うように構成可能である。そして決定された残差値の数が、第３の閾値より小さくない場合に、変位補償モジュールは、Ｂスキャンシーケンス内のＢスキャン間の軸方向変位の補償を行わないように構成可能である。

次に、本発明の例示的実施形態を以下に記載の添付図面を参照して、単なる非限定的な例として詳細に説明する。異なる図面に現れる同様の参照符号は、特に断らない限り同一又は機能的に類似の要素を表すことができる。

本明細書の第１の例示的実施形態による、Ｃスキャンデータを処理するためのデータ処理装置の概略図である。 ＯＣＴ撮像システムによる網膜の撮像中に、患者の網膜とＯＣＴ撮像システムの間の距離を変化させる患者の網膜の動きの一例を示す図である。 ＯＣＴ撮像システムにより取得されたＢスキャンのシーケンス内の隣接Ｂスキャンペアにおける、共通の視覚特徴間の軸方向シフトの一例を示す図である。 本明細書の例示的実施形態のデータ処理装置の、プログラム可能信号処理ハードウェアにおける例示的実装形態を示すブロック図である。 図１のデータ処理装置が、Ｂスキャンのシーケンスを含むＣスキャンデータを処理して、シーケンス内のＢスキャン間の軸方向変位を補償し、さらに補正データの信頼性を示す信頼性指標を生成するための方法を示すフロー図である。 本明細書に記載の第１の例示的実施形態によるＢスキャン取得中のＯＣＴ撮像システムと網膜の相対的運動を示す変動の第１周波数成分を決定するための、データ処理装置によって実行可能な方法を示す図である。 Ｂスキャンシーケンス内のＢスキャンのペアにおける共通の網膜特徴の間の軸方向シフトを示す各指標値の複数のシーケンスをプロットした図である。このプロットは、シーケンス内における隣接するＢスキャンの対応ペアの位置によって指標値が変化する様子を示す。 本明細書に記載の第２の例示的実装形態による、データ処理装置が変動の第１周波数成分を決定するために実行し得る方法を示す図である。 本明細書に記載の第３の例示的実装形態による、データ処理装置が第１周波数成分を決定するために実行し得る方法を示す図である。 本明細書に記載の例示的実施形態による、Ｂスキャンのシーケンス内のＢスキャン間の軸方向変位への補償をするか否かを決定するためにデータ処理装置が実行可能な方法を示す図である。 本明細書の第２の例示的実施形態による、Ｃスキャンデータを処理するためのデータ処理装置の概略図である。 第２の例示的実施形態のデータ処理装置がＣスキャンデータを処理して、Ｂスキャン間の軸方向変位を補償するための補正データと、生成された補正データの信頼性を示す信頼性指標とを生成するための方法を示すフロー図である。 Ｃスキャン内のＡスキャンの位置とトレーススキャンデータのＡスキャンの位置との間の例示的関係を示す図である。 Ｃスキャン内のＡスキャンの位置とトレーススキャンデータのＡスキャンの位置との間の別の例示的関係を示す図である。

［例示的実施形態１］
図１は、第１の例示的実施形態による、データ処理装置１０の概略図である。データ処理装置１０は、ＯＣＴ撮像システム３０によって取得された撮像ターゲット２０のＢスキャンのシーケンスを含むＣスキャンを処理して、Ｂスキャンのシーケンス内のＢスキャン間の軸方向変位を補償するための補正データを生成し、さらに生成された補正データの信頼性を示す信頼性指標５を生成するように構成される。

撮像ターゲット２０は、本例示的実施形態においては、患者の眼２５の網膜であってよい。ただしこれに代わって、例えば眼２５の前方領域などの眼２５の別の部分であってもよい。軸方向変位は、ＯＣＴ撮像システム３０によるＢスキャンの取得中にＯＣＴ撮像システム３０と撮像ターゲット２０の間の距離を変化させる、ＯＣＴ撮像システム３０と撮像ターゲット２０との相対運動によって生じる。軸方向変位は、網膜撮像のためのＯＣＴ撮像システム３０使用時における、眼２５へ入射するＯＣＴ光ビーム４０の伝播方向に沿う変位として理解される。

ＯＣＴ撮像システム３０は眼科用スキャナを使用してＯＣＴ撮像光ビーム４０を撮像ターゲット２０上で走査し、Ｃスキャンデータを取得し、それをデータ処理装置１０で処理する。データ処理装置１０は、本例示的実施形態のように、ＰＣやラップトップなどの単体型の処理装置として提供され得る。これがＯＣＴ撮像システム３０に通信可能に（直接またはインターネットなどのネットワークを介して）結合されて、そこからのＣスキャンデータを受信可能である。あるいは、データ処理装置３０はＯＣＴ撮像システム３０の一部として提供されてもよい。ＯＣＴ撮像システム３０は当業者に周知のＯＣＴスキャナの任意の種類のものであってよく、患者の眼２５からＯＣＴデータを取得することができる。

図１に示すように、本例示的実施形態のデータ処理装置１０は、補正データ生成モジュール２を備え、これがシーケンス内の複数の隣接Ｂスキャンペアの内の各隣接Ｂスキャンペアに対し、隣接するＢスキャン内の共通の視覚特徴のそれぞれの表示間の軸方向シフトのそれぞれの指標を決定することによって補正データを生成するように構成される。本例示的実施形態における軸方向シフトとは、隣接するＢスキャン内にある共通視覚特徴間の、軸方向（すなわち、Ｂスキャンを構成する複数のＡスキャンの各Ａスキャンのデータが配列される方向であり、このデータは眼２５に入射するＯＣＴ光ビーム４０の伝播方向に沿う様々な深さでの測定により取得される）に対応するＢスキャン軸に沿うシフト（すなわちオフセット）を表す。本例示的実施形態において、補正データ生成モジュール２によって処理されるＢスキャンのシーケンスは、ＯＣＴ撮像システム３０によって取得されたＣスキャンを形成するＢスキャンの全シーケンスを含む。ただし、補正データ生成モジュール２によって処理されるＢスキャンのシーケンスは、上に述べた全シーケンスである必要はなく、代わりに全シーケンスの部分集合、例えばＢスキャンの全シーケンス内の１つ置きのＢスキャンを含む部分集合であってもよい。

データ処理装置１０は、信頼性指標５を生成するように構成された信頼性指標生成モジュール４を更に備える。より具体的には、信頼性指標生成モジュール４は、Ｂスキャンのシーケンス内のＢスキャンの複数ペアを使用して、Ｂスキャンの複数ペアが取得されたときのＯＣＴ撮像システムに対する撮像ターゲットの相対的なそれぞれの速度又はそれぞれの加速度を示すメトリックのそれぞれの値を計算するように構成される。信頼性指標生成モジュール４はさらに、このメトリックの計算値の少なくとも所定数が閾値を超えるか否かを判定するように構成される。メトリックの計算値の少なくとも所定数が閾値を超えると判定される場合に、信頼性指標生成モジュール４はその信頼性指標５を、その補正データが信頼できないことを示す設定とするように構成される。ただし、メトリックの計算値の少なくとも所定数が閾値を超えないと判定される場合には、信頼性指標生成モジュール４はその信頼性指標５を補正データが信頼でることを示す設定とするように構成される。

データ処理装置１０は、図１に示す本例示的実施形態のように、変位補償モジュール６を更に備える。これは、ＯＣＴ撮像システム３０によるＢスキャンの取得中にＯＣＴ撮像システム３０と撮像ターゲット２０との相対運動によって生じる、Ｂスキャンのシーケンス内のＢスキャン間の軸方向変位に対して、補正データを使って補償をするように構成される。いくつかの例示的実施形態では変位補償モジュール６は、信頼性指標５が、補正データが信頼できることを示すように設定された場合にのみ前述の補償をするように構成される。

Ｃスキャンはレンダリングにより眼２５の一部の３次元画像を提供することができ、通常、眼２５の２次元領域をＯＣＴ光ビーム４０をラスターパターンなどで走査することにより取得されるＢスキャンのシーケンスを含む。各Ｂスキャンは、ＯＣＴ光ビーム４０を単一方向（例えば網膜の表面のＸ軸に沿って）走査して取得され、眼２５のその領域の２次元断面図（Ｘ軸とＺ軸に沿った）を記録する。各Ｂスキャンは複数のＡスキャンを含み、各Ａスキャンは眼２５の横方向の単一ポイントに関する、眼２５の軸方向／深さ方向（すなわちＺ軸に沿う）を表す画像データを提供する。

図２Ａは、ＯＣＴ撮像システム３０に対する撮像ターゲット２０（すなわち本例においては網膜）の動きの一例を示し、これによりＣスキャンデータ取得中に網膜とＯＣＴ撮像システム３０の間の距離の変化が生じる。図２Ａの例に示すように、Ｂスキャン取得中の患者の動きは、網膜とＯＣＴ撮像システム３０との間の距離を距離Ｄだけ（軸方向に沿って）減少させる。

図２Ｂは、共通の視覚特徴の位置の軸方向シフトを示し、これは図２Ａに示す眼２５とＯＣＴ撮像システム３０の相対的な動きの前後でキャプチャされた、符号２１０及び２２０で表す隣接Ｂスキャンペアにおいて観察される。本例において、隣接Ｂスキャンペア内の共通視覚特徴は眼２５の網膜層２３０である。図２Ｂに示すように、Ｂスキャンの取得中に発生する網膜の動きにより、Ｂスキャン２２０での網膜層２３０の位置はＢスキャン２１０での位置に対してＢスキャンのＹ軸に沿ってｄピクセルの距離だけオフセットする。図２Ｂでは、Ｂスキャン２１０、２２０のＹ軸は軸方向を表し、Ｘ軸は網膜表面に沿った横方向を表す。もし補正されなければ、隣接するＢスキャン同士における網膜層２３０の軸方向シフトは、Ｃスキャンのレンダリングにおいて、運動アーチファクトを出現させ、網膜に潜んでいる特徴の診断または測定が適切に遂行されることを妨げる可能性がある。

図３は、プログラム可能な信号処理ハードウェア３００の概略図である。これは本明細書に記載の手法を用いてＣスキャンデータを処理するように構成され、第１の例示的実施形態の補正データ生成モジュール２、信頼性指標生成モジュール４、及び変位補償モジュール６として機能することができる。

プログラム可能な信号処理装置３００が、ＯＣＴ撮像システム３０と通信してそこからＣスキャンデータを受信するための通信インタフェース（Ｉ／Ｆ）３１０を備える。信号処理装置３００はさらに、プロセッサ（例えば中央処理ユニットＣＰＵ、及び／又はグラフィック処理ユニットＧＰＵ）３２０、作業メモリ３３０（例えばランダムアクセスメモリ）、及びコンピュータ可読命令から成るコンピュータプログラム３９０を格納する命令記憶部３４０を備え、コンピュータ可読命令はプロセッサ３２０によって実行されると、本明細書に記載の補正データ生成モジュール２、信頼性指標生成モジュール４、及び変位補償モジュール６などを含む様々な機能をプロセッサ３２０に実行させる。作業メモリ３３０は、コンピュータプログラム３９０を実行中にプロセッサ３２０によって使用される情報を格納する。命令記憶部３４０は、コンピュータ可読命令を前もってロードしたＲＯＭ（例えば電気的に消去可能なプログラマブル読み出し専用メモリ（ＥＥＰＲＯＭ）又はフラッシュメモリ）を含んでもよい。あるいは、命令記憶部３４０はＲＡＭ又は類似の種類のメモリを備えてもよく、コンピュータプログラム３９０のコンピュータ可読命令を、ＣＤ－ＲＯＭ、ＤＶＤ－ＲＯＭなどの形態の、非一時的なコンピュータ可読記憶媒体３５０などのコンピュータプログラム製品から、又はコンピュータ可読命令を搬送するコンピュータ可読信号３６０から、そこに入力することができる。いずれにせよ、コンピュータプログラム３９０はプロセッサ３２０で実行されるとき、プロセッサ３２０に本明細書に記載のＣスキャンデータの処理方法を実行させる。ただし、補正データ生成モジュール２、信頼性指標生成モジュール４、変位補償モジュール６は、これに代わって、特定用途向け集積回路（ＡＳＩＣ）などの非プログラマブルハードウェアに実装されてもよいことに留意されたい。

図４は図１のデータ処理装置が、Ｂスキャンのシーケンスを含むＣスキャンデータを処理して、Ｂスキャンのシーケンス内のＢスキャン間の軸方向変位を補償するための補正データと、生成された補正データの信頼性を示す信頼性指標とを生成するための方法を示すフロー図である。

図４のステップＳ１０において、補正データ生成モジュール２が、シーケンス内の複数の隣接Ｂスキャンペアの内の各隣接Ｂスキャンペアに対し、隣接するＢスキャン内の共通の視覚特徴のそれぞれの表示間の軸方向シフトのそれぞれの指標を決定することによって補正データ生成する。軸方向シフトのそれぞれの指標は、本例示的実施形態におけるように、隣接Ｂスキャンペアにおける共通の視覚特徴のそれぞれの座標間の、上で述べた軸方向に対応するＢスキャン軸に沿った距離であってよい。ただし、軸方向シフトの指標はこれに限るものではなく、代わりに距離に基づいて決定される値であってもよい。一例として、Ｎ＋１個のＢスキャンＢ_ｋ、ここでｋ＝１，２，…，Ｎ＋１に関し、Ｂスキャンのシーケンス内のＮ個のＢスキャンのペアから、対応するＮ個の指標Ｓ_ｉ、ここでｉ＝１，２，…，Ｎの組を求めることが可能である。ここで指標Ｓ_ｉはＢスキャンのＢ_ｋ＝ｉとＢ_{ｋ＝ｉ＋１}とから求められ、これがＢスキャンのＢ_ｋ＝ｉとＢ_{ｋ＝ｉ＋１}のペア内にある共通の視覚特徴のそれぞれの表示の間の軸方向シフトを示す。指標のインデックスｉは、したがって、指標Ｓ_ｉをもたらすＢスキャンのシーケンス内の隣接Ｂスキャンペアの位置に対応する。

軸方向シフトの指標は、本例示的実施形態におけるように、隣接Ｂスキャンペア間の相互相関を計算し、計算された相互相関のピークに対応するＢスキャン間のオフセットを指標として決定することにより、各隣接Ｂスキャンペアに対して決定可能である。相互相関は、例えば、（Ｂスキャンペアの内の）片方のＢスキャン内にある視覚特徴の表示の位置の、もう一方のＢスキャン内の同一の視覚特徴の表示の位置に対するオフセット（Ｂスキャンの両方の軸に沿う）を決定する、隣接Ｂスキャンペアの規格化された２次元相互相関であってよい。ただし、軸方向に対応するＢスキャン軸に沿う決定されたオフセットのみが、隣接Ｂスキャンペア内の共通の視覚特徴のそれぞれの表示の間の軸方向シフトの指標として採用され得る。いくつかの状況においては、Ｂスキャン取得中の眼の横方向の動きが無視できる場合には、隣接Ｂスキャンペアの１次元相互相関を実行可能であり、この計算された１次元相互相関におけるピークに対応するＢスキャン間のオフセットが、軸方向シフトの指標として採用され得る。

軸方向シフトの指標は、本例示的実施形態では相互相関を使用して計算されるが、このほかの任意の適切な方法を代替的に使用可能である。例えば、代替的な実施形態では、隣接Ｂスキャンペアの内の複数のＢスキャンにおける共通の視覚的特徴のそれぞれの位置を特定し、Ｂスキャン軸に沿うその特定された位置の間での、前述した軸方向に対応する変位を決定することにより、各隣接Ｂスキャンペアに対する軸方向シフトの指標を決定してもよい。隣接Ｂスキャンペアにおける共通の視覚的特徴のそれぞれの位置は、例えば機械学習アルゴリズム、又は他の任意の適切な処理アルゴリズムを用いて特定され得る。

図４のステップＳ２０において、信頼性指標生成モジュール４はＢスキャンシーケンス内の複数のＢスキャンペアを使用して、複数のＢスキャンペアが取得されたときのＯＣＴ撮像システム３０に対する撮像ターゲット２０の速度及び／又は加速度を表すメトリックのそれぞれの値を計算する。このメトリックは、本例示的実施形態におけるように、ＯＣＴ撮像システム３０に対する撮像ターゲット２０の速度（軸方向への）を表す速度のメトリックであってよい。

速度メトリックは、本例示的実施形態のように、Ｂスキャンシーケンス内の隣接Ｂスキャンペアを用いて図４のステップＳ１０で決定された指標により計算可能である。各指標は隣接Ｂスキャンペアにある共通の視覚的特徴の軸方向シフトを示すので、指標値は、隣接Ｂスキャン間の時間で割ると、撮像ターゲットのＯＣＴ撮像システム３０に対する速度を表す。ただし、速度メトリックは隣接Ｂスキャンを使用して計算する必要はなく、代わりにＢスキャンシーケンス内の、非隣接のＢスキャンペアにある共通の視覚的特徴の軸方向シフトに基づいて計算してもよいことを理解されたい。より一般的には、Ｂスキャンペアにある共通視覚特徴の決定された軸方向シフトをＢスキャンペアのＢスキャンの撮影された時間間隔で割れば、共通視覚特徴の軸方向移動速度の指標が与えられ、したがって、撮像ターゲット２０とＯＣＴ撮像システム３０の相対運動の速度が示される。本例示的実施形態のように、このメトリックはシーケンス内の各隣接Ｂスキャンペアに対して計算可能であるが、代替的に、前に述べたように取得されたすべてのＢスキャンのペアの部分集合に対してのみ計算されてもよい。

図４のステップＳ３０において、信頼性指標生成モジュール４は、メトリックの計算値の少なくとも所定数が閾値を超えるか否かを判定する。メトリックの計算値の少なくとも所定数が閾値を超えると判定される場合に、信頼性指標生成モジュール４は、図４のステップＳ４０Ａにおいて、信頼性指標５をその補正データが信頼できないことを示すように設定する。他方で、ステップＳ３０において所定数より少ないメトリックの計算値しか閾値を超えていないと判定される場合には、信頼性指標生成モジュール４は図４のステップＳ４０Ｂにおいて、信頼性指標５をその補正データが信頼できることを示すように設定する。メトリックが速度メトリックである本例示的実施形態では、図４のステップＳ３０で使用される閾値は、ＯＣＴスキャンデータのキャプチャ中に発生することが予想可能な、撮像ターゲット２０の物理的に可能な最大速度に相当する。この閾値を超えるいかなる計算値も、撮像ターゲット２０とＯＣＴ撮像システム３０の実際の相対運動というよりもＢスキャン間の擬似相関によって生じる異常とみなすことができる。図４のステップＳ３０の所定数はしたがって、補正データが信頼できないとみなされる前に受け入れ可能な最大数に基づいて設定され得る。

本例示的実施形態ではメトリックは速度メトリックであるが、図４のステップＳ２０で計算されるメトリックは、他の形、例えばＯＣＴ撮像システム３０に対する撮像ターゲットの加速度を表す加速度メトリックであってもよい。図４のステップＳ２０において、メトリックとして加速度メトリックが用いられるときは、ステップＳ３０における閾値は、撮像ターゲット２０の加速度の現実的又は物理的に可能な最大値に基づいて設定され得る。信頼性指標生成モジュール４は、Ｂスキャンの第１のペア内の共通視覚特徴の軸方向シフト速度を示す第１の値と、Ｂスキャンの第２のペア内の共通視覚特徴の軸方向シフト速度を示す第２の値とを決定することにより加速度メトリックを評価可能である。ここで、第１のペアは第２のペアとは少なくとも１つのＢスキャンだけ異なっている。信頼性指標生成モジュール４はさらに、第１の値と第２の値の間の差に基づいて加速度メトリックを評価可能である。例えば、時刻Ｔ１と時刻Ｔ２においてキャプチャされ、Ａ１の軸方向シフト値の決定された速度を有するＢスキャンの第１のペアと、時刻Ｔ３と時刻Ｔ４においてキャプチャされ（これはＴ１とＴ２より後で発生）、Ａ２の軸方向シフト値の決定された速度を有するＢスキャンの第２のペアに関しては、加速度メトリックは次のように計算される：

ただし、加速度メトリックは上記の形に限られるものではなく、Ａ２とＡ２の間の差と、Ｂスキャンの２つのペアの間の時間的隔たりの別の尺度、例えばＴ３－Ｔ１あるいはＴ４－Ｔ２とに基づいて計算可能であることを理解されたい。

いくつかの例示的実施形態において、信頼性指標５が補正データが信頼できることを示すように設定されている場合には、図４のステップＳ５０Ｂに示すように、変位補償モジュール６が補正データを使用して、ＯＣＴ撮像システム３０と撮像ターゲット２０との相対運動により生じる、Ｂスキャンシーケンス内のＢスキャン間の軸方向変位を補償することが可能である。他方、信頼性指標５が補正データは信頼できないことを示すように設定されている場合には、ステップＳ５０Ａに示すように、変位補償モジュール６は前述の補正データを使用したＢスキャン間の軸方向変位の補償を実行できない。例えば、信頼性指標が、補正データが信頼できないことを示すように設定されている場合、Ｂスキャンのシーケンスは破棄されるか、ユーザに信頼できないことを示すかの少なくともいずれかが行われ、Ｂスキャンの新規のシーケンスから成る新規のＣスキャンデータが、撮像ターゲット２０の新規のＯＣＴスキャンを行うことによりキャプチャされ得る。Ｂスキャンの新規のシーケンスは、図４の処理を用いて処理され、新規の補正データ及び新規の信頼性指標が生成され得る。

シーケンス内の隣接するＢスキャンの対応するペアの位置に伴って、決定された指標が変動することは、雑音データ及び／又はまばたきや線により生じる画像アーチファクトで汚染されたデータにより生じる高周波成分と、その結果としての隣接するＢスキャン間の擬似相関を含む可能性がある。したがって、隣接するＢスキャンに基づいて直接生成される指標は、外れ値や異常値を何らかの形で前以ってフィルタにかけて除去することなしでは、Ｂスキャン軸方向変位の補償への使用には適さない場合がある。したがっていくつかの例示的実施形態においては、Ｂスキャン取得中のＯＣＴ撮像システム３０と撮像ターゲット２０との相対運動を示す変動の周波数成分が、ステップＳ５０Ｂで判定されて補償の実行に使用され得る。

図５は、本明細書に記載の第１の例示的実施形態による、Ｂスキャンの取得中におけるＯＣＴ撮像システム３０と撮像ターゲット２０の相対的運動を示す変動の第１周波数成分を決定するために補正データ生成モジュール２によって実行可能な方法を示すフロー図である。

図５のステップＳ２１０に示すように、補正データ生成モジュール２は、決定された指標の、シーケンス内の隣接するＢスキャンの対応するペアの位置による変動から、ＯＣＴ撮像システム３０によるＢスキャンの取得中のＯＣＴ撮像システム３０と撮像ターゲット２０の相対運動を示す変動の第１周波数成分を決定することにより補正データを生成し得る。図５のステップＳ２２０では、変位補償モジュール６は、Ｂスキャンシーケンス内のＢスキャンに第１周波数成分に基づくオフセットを適用することにより、ＯＣＴ撮像システム３０と撮像ターゲット２０の相対運動により生じるＢスキャンシーケンス内のＢスキャン間の変位に対する補償をすることが可能である。

例えば、図５のステップＳ２１０において、補正データ生成モジュール２は、指標の変動に対してｎ次多項式をフィッティングさせることにより、変動の第１周波数成分を生成可能である。そうして、図５のステップＳ２２０において、変位補償モジュール６が、ｎ次多項式に基づいて、例えばそのｎ次多項式の１以上の値を用いて、Ｂスキャンシーケンス内の１以上のＢスキャンにオフセットを与えることによって、シーケンス内の隣接Ｂスキャン間の軸方向変位を補償し得る。いくつかの例示的実施形態において、図５のステップＳ２２０において実行される補償は、信頼性指標が図４のステップＳ４０Ａにおいて補正データが信頼できることを示すように設定されていることが前提であり得ることに留意されたい。

また、第１周波数成分は図５のステップＳ２１０での多項式フィッティングによって決定される必要はなく、他の適切な方法が使用可能であることにも留意されたい。例えば、代替実施形態では、補正データ生成モジュール２は代わりに、撮像ターゲット２０とＯＣＴ撮像システム３０の相対運動に起因しない変動の高周波成分を除去するために、ステップＳ２１０において決定された指標の変動にローパスフィルタをかけて、第１周波数成分を生成してもよい。

図６は、眼２５の撮像ターゲット２０の６つの異なるＢスキャンシーケンスから決定された、６つの異なる指標の組に対応する、６つのプロット例を示す。各プロットは、対応するＢスキャンのシーケンスに関する決定された指標の変動を示し、ここで各指標は、シーケンス内の隣接Ｂスキャンペアを相互相関させることで計算され、また隣接Ｂスキャンペア間の軸方向シフトを表している。図６の例では、各隣接Ｂスキャンペアに対する決定された指標は、指標の指数ｉに対してプロットされ、これはまた隣接するＢスキャンペアに対応し、したがってＢスキャンシーケンス内の隣接するＢスキャンペアの位置を表す。

図６に示すように、Ｂスキャンシーケンス内の対応する隣接Ｂスキャンペアの位置に伴う決定された指標の変動は、様々な原因による複数の周波数成分を含み得る。例えば、Ｂスキャンシーケンスの取得中の眼２５とＯＣＴ撮像システム３０との相対運動によって生じる周波数成分に加え、変動には様々なそのほかの周波数成分が含まれ得る。例えば、網膜の曲率によって生じる低周波成分であり、それはプロット５１０、５２０、５３０に見ることができる。具体的には、網膜の曲率に起因する周波数成分は、眼２５の動きにより生じる周波数成分より低い可能性がある。それは、網膜の曲率によるＢスキャンシーケンス内の視覚特徴の軸方向シフトの速度は、患者の動きから来る軸方向シフトよりも小さい可能性があるからである。さらに、シーケンス内の対応する隣接Ｂスキャンペアの位置に伴う指標の変動は、隣接Ｂスキャン間の擬似相関から生じる高周波成分も含み得る。そのような高周波成分がプロット５４０、５５０、５６０に見られ、そのそれぞれは、これらの指標値の変動内の高周波成分の存在に寄与する多数のピークを表す。

いくつかの例示的実施形態では、撮像ターゲット２０は曲率を有し得る。これらの実施形態では、（隣接Ｂスキャンの相互相関により、あるいは隣接Ｂスキャンを使用して特徴の特定をすることにより決定される）指標を使用してＢスキャンを直接補償すると、撮像ターゲット２０がその曲率なしでレンダリングされる可能性がある。したがって、いくつかの実施形態では、視覚特徴のより正確な診断及び／又は測定のために、撮像ターゲット２０の曲率をレンダリングされるＣスキャン（補償されたＢスキャンで形成される）に保持しておくことが望ましい場合がある。

図７は、撮像ターゲット２０が曲率を有する実施形態において、図５のステップＳ２１０において第１の周波数成分を決定するために補正データ生成モジュール２が実行し得る方法を示す。図７のステップＳ３１０において、補正データ生成モジュール２が変動の第２周波数成分を決定する。これは撮像ターゲット２０（本例では網膜）の曲率を表す。補正データ生成モジュール２は、本例示的実施形態のように、決定された指標の変動にｍ次多項式のフィッティングを行って第２周波数成分を決定することが可能である。ここでｍは整数である。例えば、決定された指標のシーケンスＳ_ｉ，ｉ＝１，２，…，Ｎは、Ｎ個のデータ点（ｘ_ｉ，Ｓ_ｉ），ｉ＝１，２，…，Ｎの組として表示され得る。次いで、指標のシーケンスを表す等価なデータ点（ｘ_ｉ，Ｓ_ｉ），ｉ＝１～Ｎ、の組にｍ次多項式をフィッティングさせることによって、指標のシーケンスＳ_ｉ，ｉ＝１，２，…，Ｎにｍ次多項式Ｐ^ｍ（ｘ）をフィッティングさせることができる。ここで、ｘ_ｉは、データ点（ｘ_ｉ，Ｓ_ｉ）での指標Ｓ_ｉに関するＸ座標の値である。

図７のステップＳ３２０において、補正データ生成モジュール２が決定された指標の変動内の指標からｍ次多項式の値を引き算して、補正された指標変動を生成する。例えば、データ点形式（ｘ_ｉ，Ｓ_ｉ），ｉ＝１，２，…，Ｎで表された、決定された指標のシーケンスＳ_ｉ，ｉ＝１，２，…，Ｎ及び指標のシーケンスにフィッティングされたｍ次多項式Ｐ^ｍ（ｘ）に対して、補正された指標変動が、Ｃ_ｉ＝Ｓ_ｉ－Ｐ^ｍ（ｘ_ｉ），ｉ＝１，２，…，Ｎの値の組で与えられる。本実施例では、網膜の曲率を示す、変動の第２周波数成分は、決定された指標の変動に２次多項式をフィッティングさせることにより決定される。ただし、撮像ターゲット２０の予想される曲率によっては、異なる次数の多項式が使用される場合がある。２次多項式は、本例示的実施形態のように、指標の変動と２次多項式の間の残差値の２乗の合計を最小化する、最小二乗法を用いてフィッティング可能である。ただし、この他に任意の適切な多項式回帰法を使用可能である。

図７のステップＳ３３０において、補正データ生成モジュール２が、指標Ｃ_ｉ，ｉ＝１，２，…，Ｎの補正された変動に対してｎ次多項式をフィッティングさせることにより、指標変動の第１周波数成分を決定する。ここで、ｎはｍより大きい整数である。ＯＣＴ撮像システム３０と撮像ターゲット２０の相対運動を表す第１周波数成分は、こうしてｎ次多項式で表される。一例として、指標Ｃ_ｉ，ｉ＝１，２，…，Ｎの補正された変動はデータ点（ｘ_ｉ，Ｃ_ｉ），ｉ＝１，２，…，Ｎで表すことができ、ｎ次多項式Ｐ^ｎ（ｘ）はデータ点（ｘ_ｉ，Ｃ_ｉ），ｉ＝１，２，…，Ｎにフィッティングさせることができる。いくつかの実施形態では、擬似相関から生じるデータ点への過剰フィッティングを回避するために、補正データ生成モジュール２が指標の補正変動へｎ次多項式をフィッティングさせる前に、（例えば移動平均フィルタを使用して）平滑化操作をまず行ってもよい。本例示的実施形態では、ｎ次多項式として５次多項式が指標の補正変動にフィッティングされる。ただし、ｎの値はこれに限るものではなく、任意の適切な方法で選択可能である。いくつかの実施形態では、ｎの値は、ＯＣＴ撮像システム３０がＣスキャンデータを取り込むために使用する走査パラメータ、例えば走査密度（網膜の単位面積当たりでキャプチャされるＢスキャンの数）、及び走査継続時間（すなわちＣスキャンデータ取得にかかる総時間）など、に基づいて選択可能である。原則として、走査密度が高く走査継続時間が長いほど、ｎ次多項式に関してより大きなｎの値が選択され得る。

本実施形態では網膜の曲率は、指標変動に対する低次多項式関数のフィッティングにより決定されるが、それに代わって他の方法も使用可能である。例えば、補正データ生成モジュール２がこれに代わって、指標の変動に対する周波数領域のサンプルを決定するために、指標（ｘ_ｉ，Ｓ_ｉ），ｉ＝１，２，…，Ｎに離散的フーリエ変換を行って、網膜の曲率を表す、第２の周波数成分を決定してもよい。さらに補正データ生成モジュール２が、網膜の曲率に関係する所定の周波数域に対応する周波数領域サンプルの部分集合を抽出してもよい。この周波数域は、例えば網膜の予想される曲率を基に経験的に決定することも可能である。周波数決定モジュール４はさらに、周波数領域サンプルの部分集合に対して逆フーリエ変換を行って、第２周波数成分に対応する値を取得可能である。これを決定された指標の変動内の指標から差し引いて、補正された指標変動を生成可能である。

いくつかの例示的実施形態では、周波数領域サンプルの部分集合に対して逆フーリエ変換を行って、得られた第２周波数成分に対応する値を指標の変動から差し引く代わりに、補正データ生成モジュール２は、周波数領域サンプルの部分集合をゼロに設定してから、指標の補正変動を直接取得するためにその処理された周波数領域のサンプルに逆フーリエ変換を行ってもよい。あるいは、補正データ生成モジュール２は、補正された指標変動を決定するために指標変動に帯域フィルタをかけて網膜の曲率に対応する第２周波数成分を除去し、その後その補正された指標変動に対してｎ次多項式をフィッティングさせてもよい。帯域フィルタ処理の低域カットオフ周波数は、網膜の予想曲率に基づいて選択可能である。

ただし、走査の空間範囲が網膜の曲率に比べて小さい場合には、網膜の曲率は無視し得ると仮定でき、したがって、第２周波数成分は決定する必要がないことに留意されたい。この場合、ｎ次多項式は決定された指標のシーケンスＳ_ｉ，ｉ＝１，２，…，Ｎに直接フィッティングさせて、ＯＣＴ撮像システム３０と撮像ターゲット２０の相対運動を示す第１周波数成分として決定することが可能である。

図７のステップＳ３３０でｎ次多項式を決定すると、図７のステップＳ３４０で、ＯＣＴ撮像システム３０によるＢスキャン取得時のＯＣＴ撮像システム３０と撮像ターゲット２０の相対運動によって生じるＢスキャンシーケンス内のＢスキャン間の変位を、変位補償モジュール６が、ｎ次多項式を補正データとして使用してｎ次多項式によるオフセットをＢスキャンシーケンス内のＢスキャンに適用することで、補償可能である。

本例示的実施形態では、第１周波数成分はｎ次多項式Ｐ^ｎ（ｘ）で与えられる。したがって、変位補償モジュール６は、Ｂスキャンシーケンス内の各Ｂスキャンにこのｎ次多項式の値に基づく相殺を行ってオフセットを適用することができる。例えば、図７のステップＳ３３０において、指標Ｃ_ｉ，ｉ＝１，…，Ｎの補正された変動のシーケンスに対してフィッティングされた（あるいはより具体的にはデータ点（ｘ_ｉ，Ｃ_ｉ），ｉ＝１，２，…，Ｎの対応する組に対してフィッティングされた）ｎ次多項式Ｐ^ｎ（ｘ）に関して、ｉ＝１，２，…，Ｎに対するｎ次多項式の値のＰ^ｎ（ｘ_ｉ）の組を決定することができる。ｉ＝１，２，…，Ｎに関する値Ｐ^ｎ（ｘ_ｉ）の組と、シーケンス内のＢスキャンとの間のマッピングを使用して、値Ｐ^ｎ（ｘ_ｉ），ｉ＝１，２，…，Ｎの組の各値をＢスキャンシーケンス内のＢスキャンのオフセットに使用できる。この例では、指標Ｓ_ｉはＢスキャンＢ_ｋ＝ｉとＢ_ｋ＝ｉ＋１から求められるので、多項式Ｐ^ｎ（ｘ_ｉ）の値は、ＢスキャンのＢ_ｋ＝ｉ＋１のオフセットに使用可能である。ただし、各指標はシーケンス内の隣接するＢスキャン間の軸方向シフトを示すので、Ｂスキャンシーケンス内の各Ｂスキャンに適用されるオフセットは、そのシーケンス内の先行する各Ｂスキャン（すなわち、シーケンス内のｋより小さいインデックスを有するＢスキャン）に適用されたオフセットの累積和も含むはずである。言い換えると、シーケンス内の各ＢスキャンＢ_ｋに対し、軸方向変位を補正するために、Σ_ｉ＝１ ^ｋ－１Ｐ^ｎ（ｘ_ｉ）がＢスキャンに適用され得る。

一例として、Ｂ_ｋ，ｋ＝１～１００で表記される１００個のＢスキャンのシーケンスから成るＣスキャンに関して、ｎ次多項式Ｐ^ｎ（ｘ）が指標Ｓ_ｉ，ｉ＝１，２，…，９９の変動から決定可能である。この指標の変動Ｓ_ｉ，ｉ＝１，２，…，９９は、先ず、指標の変動から網膜の曲率を表すｍ次多項式を差し引いて補正された指標変動Ｃ_ｉ，ｉ＝１，２，…，９９を求め、次にこの補正された指標変動に対してｎ次多項式をフィッティングさせることによって、１００個のＢスキャンから計算される。変位補償モジュール６はさらに、補正された指標変動からｎ次多項式Ｐ^ｎ（ｘ）の、Ｐ^ｎ（ｘ_ｉ），ｉ＝１，２，…，９９の値の組を決定可能である。Ｐ^ｎ（ｘ_ｉ），ｉ＝１，２，…，９９の値の組はその後、ＢスキャンのＢ_ｋ，ｋ＝１～１００のそれぞれのオフセットに使用される。例えば、変位補償モジュール６は、１００個のＢスキャンの間の変位を、ＢスキャンＢ２をＰ^ｎ（ｘ１）の値だけオフセットし、ＢスキャンＢ３をＰ^ｎ（ｘ_１）＋Ｐ^ｎ（ｘ_２）の値だけオフセットすることにより、より一般的にはＢスキャンＢ_ｋをΣ_ｉ＝１ ^ｋ－１Ｐ^ｎ（ｘ_ｉ）だけオフセットすることにより、補償可能である。

図７のステップＳ３４０ではＢスキャンのオフセットは指標Ｃ_ｉ，ｉ＝１，…，Ｎの補正された変動のシーケンスにフィッティングされたｎ次多項式の値に基づいているが、他のいくつかの例示的実施形態では、指標Ｃ_ｉ，ｉ＝１，…，Ｎの補正された変動のシーケンスが第１の周波数成分とされて、Ｂスキャンシーケンス内のＢスキャンを直接オフセットさせるために使用され得ることを理解されたい。例えば、シーケンス内の各ＢスキャンＢ_ｋに対し、眼２５の軸方向変位を補正するために、Σ_ｉ＝１ ^ｋ－１Ｃ_ｉがそのＢスキャンに適用され得る。

いくつかの例示的実施形態では、指標Ｃ_ｉ，ｉ＝１，…，Ｎの補正された変動にフィッティングされたｎ次多項式は、図７のステップＳ３４０でのＢスキャンをオフセットするための第１の周波数成分としては直接には使用されない。その代わり補正データ生成モジュール２が、第２の例示的実装に従って図５のステップＳ２１０で、以下のステップから成る処理の少なくとも２つの反復を更に実行することにより第１の周波数成分を決定可能である。（ｉ）各残差値が複数の指標の変動における指標とｎ次多項式の対応する値との間の差として計算される、複数の残差値を計算するステップ。（ｉｉ）複数の残差値に、第１の正の閾値を超えるか第１の負の閾値を下回る外れ値が含まれるか否かを判定するステップ。（ｉｉｉ）複数の残差値に外れ値が含まれると判定される場合、その外れ値に対応する指標を指標の変動から除去して更新された指標変動を生成し、また、複数の残差値が外れ値を含まないと判定される場合には、その処理を終了してｎ次多項式を第１周波数成分として決定するステップ。（ｉｖ）このｎ次多項式を指標の更新された変動にフィッティングするステップ。

複数の残差値における各残差値は、上記処理の第１の反復において、補正された指標変動における指標と補正された指標変動にフィッティングされたｎ次多項式の対応する値との間の差として計算される。さらに、複数の残差値における各残差値は、上記処理の残りの１以上の反復の各反復において、処理のそれ以前の反復において生成された更新された指標変動における指標と、処理のそれ以前の反復において生成された更新された指標変動にフィッティングされたｎ次多項式の対応する値との間の差として計算される、

図８は、補正データ生成モジュール２が、第２の例示的実装に従って第１の周波数成分の決定のために実行し得る処理を示すフロー図である。図８に示すように、補正データ生成モジュール２は、図７において記述された処理のステップＳ３１０、Ｓ３２０、Ｓ３３０をまず実行してよい。ただし、第２の例示的実装では、ステップＳ３３０においてｎ次多項式が決定された後、ｎ次多項式はステップＳ３４０においてＢスキャンのオフセットには直接使用されない。代わりに、補正データ生成モジュール２がさらに、図８のステップＳ４４０において各残差値が複数の指標の変動における指標とｎ次多項式の対応する値との間の差として計算される、複数の残差値を計算する。この例においては、補正された指標変動Ｃ_ｉ，ｉ＝１，２，…，Ｎは、データ点（ｘ_ｉ，Ｃ_ｉ），ｉ＝１，２，…，Ｎとして表され、ｎ次多項式はＰ^ｎ（ｘ）として表されるが、複数の残差値は、図７のステップＳ４４０において、Ｉ＝１，２，…，Ｎに対してＲ（ｘ_ｉ）＝Ｃ_ｉ－Ｐ^ｎ（ｘ_ｉ）として決定可能である。

図８のステップＳ４５０において、補正データ生成モジュール２は、第１の正の所定の閾値より高いか第１の負の所定の閾値より低い外れ値が複数の残差値に含まれるか否かを判定する。複数の残差値が外れ値を含むという判定に応答して、補正データ生成モジュール２が、図８のステップＳ４６０において外れ値に対応する補正された指標Ｃ_ｉを補正された指標変動Ｃ_ｉ，ｉ＝１，２，…，Ｎから除去して、ｘ＝１～Ｎに対して更新された指標変動Ｕｉを生成する。外れ値に対応する補正された指標Ｃ_ｉを除去することは、例えば、補正された指標Ｃ_ｉの値をゼロに設定することを含む。例えば、補正された指標Ｃ_ｉ＝ｐに対応する残差値が外れ値であると判定されると、データ点（ｘ_ｐ，Ｃ_ｐ）は、補正された指標変動を表すデータ点（ｘ_ｉ，Ｃ_ｉ），ｉ＝１，２，…，Ｎの組においてゼロに設定可能である。ただし、いくつかの実施形態では、外れ値に対応する補正された指標Ｃ_ｐをゼロに設定する代わりに、補正データ生成モジュール２が、外れ値に対応する補正された指標Ｃ_ｉを新しい値に代わりに置き換えることができる。この値はステップＳ４６０において外れ値としては判定されない残差値を有するＣ_ｉの１以上の他の値に基づくものであって、例えばそれらの値の内挿により得られるものであってよい。

図８のステップＳ４７０において、補正データ生成モジュール２が、ｎ次多項式をｘ＝１，２，…，Ｎに対する更新された指標変動Ｕｉにフィッティングする。そうして、ステップＳ４４０の処理を反復して、今度は複数の残差値を計算する。ここで、各残差値は、（ステップＳ４６０で決定された）更新された指標変動内の指標と、ステップＳ４７０で更新された指標変動に対してフィッティングされたｎ次多項式の対応する値との差として計算される。ステップＳ４７０が実行された後、反復するように決定された複数の残差値が外れ値を含まないことがステップＳ４５０において判定されるまで（あるいは代替の実装においては、所定数より少ない外れ値の数が含まれるまで）、図８のステップＳ４４０、Ｓ４５０、Ｓ４６０、Ｓ４７０が反復されてよい。ステップＳ４５０において、複数の残差値において外れ値がないと判定される場合は、処理はステップＳ４５０Ａに進み、そこで、複数の残差値の生成に使用されたｎ次多項式が第１周波数成分として採用される。それが次に図７のステップＳ３４０において、網膜とＯＣＴ撮像システム３０との（軸方向の）相対運動により生じるＢスキャンシーケンス内のＢスキャン間の変位を補償するためにＢスキャンのオフセットに使用される。具体的には、補正データ生成モジュール２は、図７のステップＳ３４０に関して前に説明したようにして、ｎ次多項式の値に基づいてＢスキャンをオフセットさせることが可能である。

いくつかの例示的実施形態では、図８のステップＳ４５０で、反復のための複数の残差値がいかなる外れ値も含まないと判定されたとき、複数の残差成分を取得するために使用されたｎ次多項式を第１周波数成分として採用可能であり、図７のステップＳ３４０でＢスキャンのオフセットへの使用の可否を判定するために、補正データ生成モジュール２が、複数の残差値を更に評価してもよい。より具体的には、図９を参照すると、補正データ生成モジュール２が、図９のステップＳ５１０において、複数の残差値の中で、第２の正の閾値より大きいか第２の負の閾値より小さい大きさを有する残差値の数を決定し得る。ここで、第２の正の閾値は、図８のステップＳ４５０における第１の正の閾値より小さく、また第２の負の閾値は図８のステップＳ４５０における第１の負の閾値より大きい。図９のステップＳ５２０において、補正データ生成モジュール２は、図９のステップＳ５１０において決定された残差値の数が第３の閾値を超えるか否かを判定する。図９のステップＳ５１０において判定された残差値の数が第３の閾値を超えないと判定されることに応じて、補正データ生成モジュール２が、図９のステップＳ５３０Ａにおいて、第１周波数成分として複数の残差値（これは外れ値を含まない）の生成に使用されたｎ次多項式を選択し、第１周波数成分としてｎ次多項式使用することにより図７のステップＳ３４０を実行する。他方、補正データ生成モジュール２が、図９のステップＳ５１０で判定された残差値の数が第３の閾値を超えると判定する場合、図７のステップＳ３４０は実行されない。これは、多数の残差値が比較的大きい場合には、これは補正された指標変動（あるいは更新された指標変動）へのｎ次多項式のフィッティングがよくないことを示すものと受け取ることができるからである。したがって、ｎ次多項式を用いて、Ｂスキャンシーケンスにおける軸方向変位を補正する代わりに、そのＢスキャンのシーケンスを破棄して、撮像ターゲット２０への新規のＯＣＴスキャンを実行することにより、新規のＢスキャンシーケンスを含む新規のＣスキャンをキャプチャしてもよい。Ｂスキャンの新規のシーケンスは既述の任意の方法で処理して、ＯＣＴ撮像システム３０と撮像ターゲット２０との相対運動によって生じる、Ｂスキャンシーケンス内のＢスキャン間の軸方向変位を補償するための補正データを生成し得る。
［例示的実施形態２］

上記の第１の例示的実施形態とその変形例においては、補正データはＢスキャンシーケンス内の隣接Ｂスキャンペアから導くことができる。それに対し、この例示的実施形態のデータ処理装置では、Ｂスキャン間の軸方向変位を補償するために使用される補正データは、補正データ決定のために取得されたＢスキャンに対して比較可能な参照データとして機能するトレーススキャンデータを使用して生成される。つまり、シーケンス内の複数のＢスキャンの内の各隣接Ｂスキャンペアに対し、その隣接するＢスキャン内の共通の視覚特徴のそれぞれの表示間での相対的な軸方向シフトのそれぞれの指標を決定することによって補正データを生成するのではなく、それに代わって、シーケンス内の少なくともいくつかのＢスキャンの各Ｂスキャンに関して、Ｂスキャンとトレーススキャンデータ内の共通の視覚特徴のそれぞれの表示の間での軸方向シフトのそれぞれの指標によって、補正データを決定可能である。

図１０は本明細書の第２の例示的実施形態による、Ｃスキャンデータ処理のためのデータ処理装置１０’の概略図である。データ処理装置１０’は、撮像ターゲット２０を横断して延在する平行な走査線に沿って走査するＯＣＴ撮像システム３０によって取得された、撮像ターゲット２０のＢスキャンシーケンスから成るＣスキャンデータを、トレーススキャンデータと共に処理し、Ｂスキャンシーケンス内のＢスキャン間の軸方向変位を補償するための補正データを生成し、かつ生成された補正データの信頼性を示す信頼性指標５’を生成するように構成される。

トレーススキャンデータは、ＯＣＴ撮像システム３０が上記の走査線と交差するトレース走査線に沿って走査を実行することにより取得される。第１の例示的実施形態のように軸方向変位は、ＯＣＴ撮像システム３０によるＢスキャンの取得中に、その間の距離を変化させるＯＣＴ撮像システム３０と撮像ターゲット２０との相対運動によって生じる。図１０に示すように、データ処理装置１０’は、補正データ生成モジュール２’と信頼性データ生成モジュール４’とを備える。さらに、データ処理装置１０’は本例示的実施形態のように、変位補償モジュール６’を備えてもよい。

図１１は、図１０のデータ処理装置１０’がＢスキャンシーケンスから成るＣスキャンを処理して、補正データと信頼性指標５’を生成する方法を示すフロー図である。図１１のステップＳ１０’において補正データ生成モジュール２’が、シーケンス内の少なくともいくつかのＢスキャンの各Ｂスキャンに対し、そのＢスキャンとトレーススキャンデータ内の共通の視覚特徴のそれぞれの表示の間の軸方向シフトのそれぞれの指標を決定することにより補正データを生成する。

図１２は、第１の例による、Ｃスキャンを形成するＢスキャンシーケンスのＢスキャン内のＡスキャンの位置と、トレーススキャンデータからのＡスキャンとの間の関係を示す図である。図１２において、マトリックス１２００は、Ｃスキャンでカバーされる画像化ターゲットの横方向のすべての位置にわたるＣスキャンを上から見た図を表す。マトリックス１２００内の各セルはＣスキャン内のＡスキャンに対応する。さらに、マトリックス１２００の各行はＣスキャン内のＢスキャンに対応する。マトリックス１２００の複数の行に対応するＢスキャンのシーケンスは、撮像ターゲット２０に沿って延在する複数の走査線に沿ってＯＣＴスキャンを実行することにより得られる。この点を図示するために、図１２ではマトリックス１２００に重ねてトレース走査線１２２０を示す。これに沿ってＯＣＴスキャンが実行されて、対応するトレーススキャンデータが生成される。さらには図１２において、灰色の各セル１２０５がトレーススキャンデータのＡスキャンを表し、マトリックス１２００で表されるＣスキャン上に重ねられる。図１２のトレーススキャンデータは撮像ターゲット２０を横断して延在するトレース走査線に沿ってＯＣＴスキャンを実行することにより生成される。

図１２に示すように、トレース走査線１２２０は複数の走査線１２１０と交差し、マトリックス１２０００内の各Ｂスキャン（すなわちマトリックス１２００の行）が、トレーススキャンデータの対応するＡスキャンと重なるか交差するＡスキャンを有するようになっている。したがって、ＢスキャンのＡスキャンと、トレーススキャンデータの対応するＡスキャンとを使用して、Ｃスキャン内のＡスキャンをトレーススキャンデータ内の対応する場所のＡスキャンに写像するマッピングに基づいて、ステップＳ１０’で説明した指標を決定可能である。補正データ生成モジュール２’はこうして、Ｂスキャンシーケンスの１つのＢスキャン内のＡスキャンと対応する位置にあるトレーススキャンデータのＡスキャンとを相互相関させることにより、ステップＳ１０’で指標を決定することができる。あるいは、補正データ生成モジュール２’は共通特徴のそれぞれの表示間の（軸方向への）オフセットを判定する前に、適切な特徴識別アルゴリズムを使用して、ＢスキャンのＡスキャンと、対応する位置にあるトレーススキャンデータのＡスキャンとにある共通特徴の表示を識別してもよい。

ただし、トレース走査線は図１２に示すような直線に限るものではなく、Ｂスキャンシーケンス生成のためにスキャンされる走査線と交差（又は重複）しさえすれば、画像ターゲット２０を横切る任意の軌跡に沿ってより一般的に延在し得ることに留意されたい。

図１３は、円形のトレーススキャンが使用される別の例における、Ｃスキャン内のＡスキャンの位置と、トレーススキャンデータを形成するＡスキャンの位置との間の関係を示す。マトリックス１２００は、図１２に関する説明と同様にＣスキャンを表す。ただし、図１３の例においては、トレーススキャンデータは、ＯＣＴ撮像光により円形状に走査される、撮像ターゲット２０上の複数の横方向位置に対応する複数のＡスキャン１３１０からなる。言い換えると、ＯＣＴ撮像システム３０がＯＣＴ撮像光ビームを走査してトレーススキャンデータを生成するトレース走査線は、この例では撮像ターゲット２０の表面に沿って円形軌跡を描く。さらには、図１３に示すように、トレーススキャンデータを形成するＡスキャンの密度によっては、トレーススキャンデータの複数のＡスキャンが、Ｃスキャンを形成するＢスキャンシーケンス内のＢスキャンのＡスキャンと重なり得る。そのような場合、複数のＡスキャンのそれぞれに対して指標を計算可能であり、例えば、複数のＡスキャンにより得られる指標の平均を計算してＢスキャンに対する指標として使用可能である。

ここで再び図１１を参照すると、図１１のステップＳ１０’で補正データを決定すると、信頼性データ生成モジュール４’が、ステップＳ２０’とＳ３０’、並びにＳ４０Ａ’とＳ４０Ｂ’の内の１つを実行することにより信頼性指標５’を生成する。より具体的にはステップＳ２０’において補正データ生成モジュール２’が、少なくともいくつかのＢスキャン内のＢスキャンのペア２１０、２２０に対して決定された指標を用いて、このＢスキャンのペア２１０、２２０が取得されたときのＯＣＴ撮像システム３０に対する画像ターゲット２０のそれぞれの速度又はそれぞれの加速度を表すメトリックのそれぞれの値を計算することにより、信頼性指標５’を生成する。

このメトリックは、本例示的実施形態におけるように、ＯＣＴ撮像システム３０に対する撮像ターゲット２０の速度（軸方向への）を表す速度のメトリックであってよい。速度メトリックは、本例示的実施形態のように、図１１のステップＳ１０’で決定された複数の指標から指標の１つのペアを用いて計算可能であり、ここでこのペアの指標はＢスキャンシーケンス内の２つのＢスキャンに対応する。具体的には、２つの指標の値の差を、この２つの指標に対応する２つのＢスキャンを取得する間の時間で割ると、ＯＣＴ撮像システム３０に対する撮像ターゲット２０の速度となる。この２つの指標は、２つの隣接するＢスキャンの指標である必要はない。代わりにＢスキャンシーケンス内の２つの隣接しないＢスキャンに対応する指標であってもよい。

いくつかの実施形態において、速度メトリックの代わりに、図４’のステップＳ２０’で計算されるメトリックは、代わりに加速度メトリックであってもよい。これはＯＣＴ撮像システム３０に対する撮像ターゲットの加速度を表す。信頼性指標生成モジュール４’は、Ｂスキャンの第１のペア内の共通視覚特徴の軸方向シフト速度を示す第１の値と、Ｂスキャンの第２のペア内の共通視覚特徴の軸方向シフト速度を示す第２の値とを決定することにより加速度メトリックを計算可能である。ここで、第１のペアは第２のペアとは少なくとも１つのＢスキャンだけ異なっている。本例においては、第１の値は、ステップＳ１０’で決定された複数の指標の内のＢスキャンの第１のペアに対応する指標の第１のペアを用いて決定可能である。第１の値は、例えば指標の第１のペアの間の差として計算され得る。同様に、第２の値は、ステップＳ１０’で決定された複数の指標の内のＢスキャンの第２のペアに対応する指標の第２のペアを用いて決定可能である。第２の値は、例えば指標の第２のペアの間の差として計算され得る。信頼性指標生成モジュール４’はさらに、第１の値と第２の値の間の差に基づいて加速度メトリックを評価可能である。例えば、時刻Ｔ１と時刻Ｔ２においてキャプチャされ、Ａ１の軸方向シフト値の決定された速度を有するＢスキャンの第１のペアと、時刻Ｔ３と時刻Ｔ４においてキャプチャされ（これはＴ１とＴ２より後で発生）、Ａ２の軸方向シフトの決定された速度を有するＢスキャンの第２のペアに関しては、加速度メトリックは次のように計算可能である：

ただし、加速度メトリックは上記の形に限られるものではなく、Ａ２とＡ２の間の差と、Ｂスキャンの２つのペアの間の時間的隔たりの別の尺度、例えばＴ３－Ｔ１あるいはＴ４－Ｔ２とに基づいて計算可能であることを理解されたい。ステップＳ２０’において、メトリックが加速度メトリックである場合には、ステップＳ３０’で設定される閾値は、現実的又は物理的に可能な、加速度の最大値に基づいて設定され得る。

ステップＳ３０’において、信頼性指標生成モジュール４’は、メトリックの計算値の少なくとも所定数が閾値を超えるか否かを判定する。ステップＳ３０’は図１のステップＳ３０と同じであり、唯一の違いは、Ｓ２０’でのメトリックの計算方法がステップＳ２０でのメトリックに比べて異なるために、ステップＳ３０’で使用される閾値及び所定数の値が図１のステップＳ３０における対応する値とは異なり得ることである。

図１１に戻ると、メトリックの計算値の少なくとも所定数が閾値を超えると判定される場合に、信頼性指標生成モジュール４’は、ステップＳ４０Ａ’において、信頼性指標５’をその補正データが信頼できないことを示すように設定する。ただし、メトリックの計算値の少なくとも所定数が閾値を超えないと判定される場合には、信頼性指標生成モジュール４’は信頼性指標５’をその補正データが信頼できることを示すように設定する。

いくつかの例示的実施形態において、補正データが信頼性指標５’は信頼できることを示すように設定されている場合には、図１１のステップＳ５０Ｂ’に示すように、変位補償モジュール６’は補正データを使用して、ＯＣＴ撮像システム３０と撮像ターゲット２０との相対運動により生じる、Ｂスキャンシーケンス内のＢスキャン間の軸方向変位を補償することが可能である。他方、補正データが信頼性指標５’は信頼できないことを示すように設定されている場合には、変位補償モジュール６’は前述の補正データを使用したＢスキャン間の軸方向変位の補償を実行しない。

図１０のデータ処理装置１０’は、図１の装置とは以下の点で異なっていることを理解されたい。すなわち、シーケンス内の複数の隣接するＢスキャンの内の各隣接Ｂスキャンペアに対し、その隣接するＢスキャン内の共通の視覚特徴のそれぞれの表示間での軸方向シフトのそれぞれの指標を決定することによって補正データを生成するのではなく、それに代わって、シーケンス内の少なくともいくつかのＢスキャンの各Ｂスキャンに対して、Ｂスキャンとトレーススキャンデータ内の共通の視覚特徴のそれぞれの表示の間での軸方向シフトのそれぞれの指標を決定することによって、補正データが生成される。これに関し、上で説明したデータ処理装置１０のすべての実装の詳細は、データ処理装置１０’にもまた適用可能である。例えば、図５を参照して説明した処理は、データ処理装置１０’にもまた適用可能である。ただし、本例示的実施形態においては、第１周波数成分は、決定された指標の、シーケンス内の対応するＢスキャンの位置による変動から決定可能であることだけが異なっている。さらに、多項式が指標にフィッティングされて第１周波数成分を決定する場合、図１の補正データ生成モジュール２に対比して、補正データ生成モジュール２’がその指標を違う方法で生成するとすれば、本例示的実施形態では異なる次数の多項式が選択され得る。

データ処理装置１０’に関しては、図７のステップＳ３１０からステップＳ３３０に関連して説明したような、撮像ターゲット２０の曲率を表す第２周波数成分の決定を実行する必要がない。それは、トレーススキャンデータを用いて生成した補正データに基づく補償の実行は、撮像ターゲット２０のいかなる曲率の保持も許容し得るからである。これは、曲率が、トレーススキャンデータとＢスキャンのシーケンスの両方でキャプチャされるからであって、トレーススキャンデータとＢスキャンシーケンスの間の任意の計算された軸方向シフトは、スキャン取得中の相対運動による結果と前述したような擬似相関の結果のみのものであり得るからである。ただし、図８のステップＳ４４０からＳ４７０、及びこれらのステップに関連して上で説明したすべての例及び変形もまた、データ処理装置１０’に適用可能である。より具体的には、データ処理装置１０’に関しては、図１１のステップＳ１０’で生成された、変動指標に対してｎ次多項式をフィッティングさせることによって第１周波数成分を決定した後、データ処理データ装置１０’が、以下のステップを含む処理を少なくとも２回反復することにより、補正データを生成可能である。（ｉ）各残差値が軸方向シフトの決定された指標の変動における指標と、その指標の変動にフィッティングされたｎ次多項式の対応する値との間の差として計算される、複数の残差値を計算するステップ。（ｉｉ）上記の複数の残差値に、第１の正の閾値より高いか第１の負の閾値より低い外れ値が含まれるか否かを判定するステップ。（ｉｉｉ）上記の複数の残差値に外れ値が含まれると判定される場合、その外れ値に対応する指標を指標の変動から除去して更新された指標変動を生成し、また、残差が外れ値を含まないと判定される場合には、ｎ次多項式を補正データとして決定して処理を終了するステップ。（ｉｖ）このｎ次多項式を指標の更新された変動にフィッティングするステップ。

複数の残差値における各残差値は、上記処理の第１の反復において、補正された指標変動における指標と、補正された指標変動にフィッティングされたｎ次多項式の対応する値との間の差として計算される。さらに、複数の残差値における各残差値は、上記処理の残りの１以上の反復の各反復において、処理のそれ以前の反復において生成された、更新された指標変動における指標と、処理のそれ以前の反復において生成された更新された指標変動にフィッティングされたｎ次多項式の対応する値との間の差として計算される。

さらに、図１１のＳ４０Ｂ’において信頼性指標が補正データは信頼できることを示すように設定されている場合、補正データ生成モジュール２’が、複数の残差値の中で、第２の正の閾値より大きいか第２の負の閾値より小さい大きさの残差値の数を決定し得る。ここで第２の正の閾値は第１の正の閾値より小さく、第２の負の閾値は第１の負の閾値より大きい。この決定された残差値の数が第３の閾値より小さい場合には、変位補償モジュール６’が、Ｂスキャンシーケンス内のＢスキャンに対する補正データに基づいてオフセットを適用することによって、Ｂスキャンシーケンス内のＢスキャン間の軸方向変位に対する補償を行い得る。ただし、決定された残差値の数が、第３の閾値より小さくない場合に、変位補償モジュール６’は、Ｂスキャンシーケンス内のＢスキャン間の軸方向変位の補償を行わない場合がある。

ここで述べた例示的態様は制約、特に従来のＯＣＴデータ処理に関するコンピュータ技術に根差す制約を回避する。この制約は、一般に体積ＯＣＴデータ取得中の患者の不随意的な動きによって生じる運動アーチファクトを表示する、レンダリングされた体積ＯＣＴデータをもたらす可能性がある。これらの運動アーチファクトは、例えば視覚特徴の識別及びそれに続く任意の診断測定の精度に負の影響を与える可能性がある。本明細書に記載の例示的態様によって、Ｃスキャンを構成するＢスキャンシーケンス内のＢスキャン間の軸方向変位、これはＯＣＴ撮像システムと撮像ターゲットとの相対運動によって生じ得るものであるが、これを補償する補正データが生成される。具体的には、Ｂスキャンシーケンス内の隣接Ｂスキャンの各ペアに対して、その隣接Ｂスキャン内の共通視覚特徴のそれぞれの表示の間の軸方向シフトのそれぞれの指標が決定される。さらに、補正データの信頼性を表す信頼性指標もまた生成される。信頼性指標は、相対運動の速度又は加速度を表すメトリックに基づく。少なくともいくつかの実施形態において、補正データは、信頼性指標がその補正データが信頼できることを示す場合にのみ、Ｂスキャン軸方向変位の補償に使用される。したがって、信頼性のある補正データのみが軸方向変位の補償に使用可能であるので、Ｂスキャン間の軸方向変位を補償するＢスキャンの処理は改善することができる。少なくともいくつかの実施形態において、シーケンス内の隣接Ｂスキャンの複数のペアの、各隣接Ｂスキャンペアに対して、隣接Ｂスキャン内の共通の視覚特徴のそれぞれの表示の間の軸方向シフトに対する個々の指標を決定して補正データを生成するのではなく、それに代わって、シーケンス内の少なくともいくつかのＢスキャンの各Ｂスキャンに対して、Ｂスキャンとトレーススキャンデータ内の共通の視覚特徴のそれぞれの表示の間での軸方向シフトの個々の指標を決定することによって、補正データが決定可能である。さらに、少なくともいくつかの例示的実施形態において、決定された指標変動にフィッティングされた多項式の残差を用いて補正データの品質を推定し、そうして、その補正データを用いて補償を行うべきか否かを判定する。また、本明細書に記載の例示的態様の上記の機能は、これはコンピュータ技術に根差すものであるが、それによって、本明細書に記載の例示的態様は、コンピュータ及びコンピュータ処理／機能性を改善し、また少なくとも画像処理、光干渉断層計（ＯＣＴ）とデータ処理、並びにＯＣＴ画像データの処理の分野もまた改善する。

上記の説明においては、例示的態様をいくつかの例示的実施形態を参照して説明した。したがって本明細書は、制限的ではなく、例示的であるとみなされるべきである。同様に、例示的実施形態の機能及び利点を強調する、図面に示す形状は例示のみを目的として提示される。例示的実施形態のアーキテクチャは十分に柔軟かつ設定可能であり、添付図に示すもの以外の方法で利用可能である。

本明細書に提示した例のソフトウェア実施形態は、命令又は命令のシーケンスを有する一つ又は複数のプログラムなどのコンピュータプログラム又はソフトウェアとして提供可能である。これらは一実施形態においては、機械でアクセス可能若しくは機械で読み取り可能な媒体、命令記憶部、又はコンピュータ可読記憶装置などの製品に含まれるか格納され、このそれぞれが非一時的であってよい。非一時的な機械アクセス可能媒体、機械可読媒体、命令記憶部、又はコンピュータ可読記憶装置上のプログラム又は命令は、コンピュータシステム又は他の電子デバイスのプログラムに使用可能である。機械可読媒体又はコンピュータ可読媒体、命令記憶部、及び記憶装置には、これに限らないが、フロッピディスク（登録商標）、光ディスク、及び光磁気ディスク、あるいは、電子的命令を格納若しくは送信するのに適した他の種類の、媒体／機械可読媒体／命令記憶部／記憶装置が含まれ得る。本明細書に記載の技術は、いかなる特定のソフトウェア構成にも限定されない。これらは任意のコンピューティング環境又は処理環境への適用性を見出し得る。本明細書に使用されている「コンピュータ可読」、「機械アクセス可能媒体」、「機械可読媒体」、「命令記憶部」及び「コンピュータ可読記憶装置」という用語は、機械、コンピュータ又はコンピュータプロセッサによって実行するために命令又は命令シーケンスを格納、コード化、又は送信することが可能で、かつ本明細書に記載の任意の方法を、機械／コンピュータ／コンピュータプロセッサに実行させる、任意の媒体を含むものとする。さらに、当分野においては、ソフトウェアは、ある形態又は別の形態であっても（例えばプログラム、手順、処理、アプリケーション、モジュール、ユニット、ロジックなど）、動作を起こし、結果をもたらすものを指すのが一般的である。そのような表現は、処理システムによるソフトウェアの実行が、プロセッサに動作を実行させて結果を生成させることの簡略的な記述に過ぎない。

いくつかの実施形態は、特定用途向け集積回路、フィールドプログラマブルゲートアレイを用意することによって、又は従来のコンポーネント回路の適切なネットワークを相互接続することによっても、実装することができる。

いくつかの実施形態にはコンピュータプログラム製品が含まれる。コンピュータプログラム製品は、命令が格納された記憶媒体、命令記憶部、又は記憶装置であってよい。それを使用して、コンピュータ又はコンピュータプロセッサが本明細書に記載の例示的実施形態の任意の手順を実行することを制御し、又は実行させることができる。記憶媒体／命令記憶部／記憶装置には、これに限るものではないが一例として、光ディスク、ＲＯＭ、ＲＡＭ、ＥＰＲＯＭ、ＥＥＰＲＯＭ、ＤＲＡＭ、ＶＲＡＭ、フラッシュメモリ、フラッシュカード、磁気カード、光カード、ナノシステム、分子メモリ集積回路、ＲＡＩＤ、リモートデータストレージ／アーカイブ／ウェアハウジング、及び／又は命令及び／又はデータの格納に好適な他の任意の種類のデバイスが含まれ得る。

コンピュータ可読媒体、命令記憶部又は記憶装置のいずれかに格納されたいくつかの実装には、システムのハードウェアを制御し、システム又はマイクロプロセッサに人間のユーザ又は本明細書に記載の例示的実施形態の結果を利用する他の機構との対話可能とする、ソフトウェアが含まれる。そのようなソフトウェアには、これに限らないが、デバイスドライバ、オペレーティングシステム、及びユーザアプリケーションが含まれ得る。究極的には、そのようなコンピュータ可読媒体又は記憶装置には、上で説明したような本発明の例示的態様を実行するためのソフトウェアが更に含まれる。

プログラミング及び／又はシステムのソフトウェアには、本明細書に記載の手順を実行するためのソフトウェアモジュールが含まれる。本明細書のいくつかの例示的実施形態において、モジュールにはソフトウェアが含まれる。ただし本明細書の他の例示的実施形態ではモジュールにはハードウェア又はハードウェアとソフトウェアの組み合わせが含まれる。

以上本発明の様々な例示的実施形態を述べたが、これらは例示のために提示したものであり、限定するためのものではないことを理解されたい。形式及び詳細において様々な変更をなし得ることは当業者には明らかであろう。したがって、本発明は上記のいかなる例示的実施形態によっても制限されるものではなく、以下に述べる特許請求の範囲及びその均等物によってのみ規定されるべきである。

さらに、要約は、特許庁及び一般人、特に特許又は法律用語又は文体に精通していない科学者、技術者及び実務者が、大まかな調査で本出願の技術的開示の本質及び要点を迅速に判定できるようにすることを目的とするものである。要約は、いずれにしても本明細書に提示する例示的実施形態の範囲に関して制限的であることを意図するものではない。また、特許請求の範囲に説明される任意の手順は、必ずしも提示された順番で遂行されることを必要としないことも理解されたい。

本明細書には多くの特定の実施形態の詳細が含まれているが、それらは任意の発明又は特許請求されるものの範囲を制限するものとみなすべきではなく、むしろ本明細書に記載の具体的な実施形態に固有の特徴の記述として理解されるべきである。本明細書において別々の実施形態の文脈で記載されている特定の特徴は、組み合わせて単一の実施形態に実装することもできる。その逆に、単一の実施形態の文脈で記載されている様々な特徴は、複数の実施形態に別々に、又は任意の適切な部分的な組み合わせにして実装することもできる。さらに、上記において特徴が特定の組み合わせで作用するものとして記述され、また最初にそのように特許請求されているとしても、特許請求された組合せからの１つ以上の特徴は場合によってはその組み合わせから除外することが可能であり、特許請求された組み合わせは部分的な組み合わせ又は部分的な組み合わせの変形に向けることも可能である。

特定の状況下では、マルチタスク及び並列処理をすることも有利であり得る。さらに、上記の実施形態における様々な構成要素の分離は、すべての実施形態においてそのような分離を必要とすると理解するべきではなく、また、上記のプログラム要素及びシステムは一般に単体のソフトウェア製品に統合し、あるいは複数のソフトウェア製品にパッケージ化され得ることを理解されたい。

ここまでにいくつかの例示的実施形態及び実施形態を説明したが、前述したものは例示であって限定ではなく、一例として提示されたことは明らかである。特に、本明細書に示した実施例の多くは、装置又はソフトウェア要素の特定の組み合わせを含むが、これらの要素は同じ目的を達成するために違う形で組み合わされてもよい。一つの実施形態に関連してのみ議論した作用、要素及び特徴は、他の実施形態における同様の役割から除外されることを意図するものではない。

本明細書に記載の装置は、その特性から逸脱することなく他の特定の形態で具体化されてもよい。前述の実施形態は、記述したシステム及び方法を限定するためではなく、説明するためのものである。本明細書に記載の装置の範囲は、したがって前述の説明よりもむしろ添付の特許請求の範囲に示され、特許請求の範囲の意味するところ及び等価物の範囲内にある変更は、その中に含まれる。

Claims (15)

1. 光干渉断層計（ＯＣＴ）撮像システムにより取得された撮像ターゲットのＢスキャンシーケンスからなるＣスキャンデータを、前記ＯＣＴ撮像システムによりＢスキャンを取得中にその間の距離を変動させる前記ＯＣＴ撮像システムと前記撮像ターゲットとの相対運動により生じる、前記Ｂスキャンシーケンス内のＢスキャン同士間の軸方向変位を補償するための補正データを生成するために、及び、さらに前記生成された補正データの信頼性を示す信頼性指標を生成するために、処理する方法であって、
   前記Ｂスキャンシーケンス内の複数の隣接Ｂスキャンペアの内の各隣接Ｂスキャンペアに対し、前記隣接するＢスキャン内の共通の視覚特徴のそれぞれの表示の間の軸方向シフトのそれぞれの指標を決定することによって前記補正データを生成することと、
   前記信頼性指標を生成することであって、
   前記Ｂスキャンシーケンス内のＢスキャンの複数ペアを使用して、前記Ｂスキャンの複数ペアが取得されたときの前記ＯＣＴ撮像システムに対する前記撮像ターゲットの相対的なそれぞれの速度又はそれぞれの加速度を示すメトリックのそれぞれの値を計算し、
   前記メトリックの少なくとも所定数の前記計算値が、閾値を超えるか否かを判定し、
   前記メトリックの少なくとも前記所定数の計算値が前記閾値を超えると判定される場合に、前記信頼性指標を前記補正データが信頼できないことを示すように設定し、
   メトリックの少なくとも所定数の計算値が閾値を超えないと判定される場合に、信頼性指標を、補正データが信頼できることを示すように設定する、
   ことによって、信頼性指標を生成することを含む方法。
2. 前記補正データは、前記Ｂスキャンシーケンス内の前記隣接するＢスキャンの対応するペアの位置による前記決定された指標の変動から、前記ＯＣＴ撮像システムによる前記Ｂスキャンの取得中における前記ＯＣＴ撮像システムと前記撮像ターゲットの相対運動を示す、前記変動の第１周波数成分を決定することにより生成される、請求項１に記載の方法。
3. 前記軸方向シフトの前記それぞれの指標は、前記各隣接Ｂスキャンペアに対して、
   隣接するＢスキャンペア間の相互相関を計算し、計算された相互相関におけるピークに対応するＢスキャンの間のオフセットを指標として決定するか、あるいは、
   前記隣接するＢスキャンペアの前記Ｂスキャンにおける前記共通の視覚特徴のそれぞれの位置を特定し、前記撮像システムの軸方向を表す、Ｂスキャン軸に沿う前記特定された位置同士の間の変位を決定する、
   ことにより決定される、請求項１又は請求項２に記載の方法。
4. 前記撮像ターゲットは曲率を有し、前記方法は、前記撮像ターゲットの前記曲率を示す、前記変動の第２周波数成分を決定することを更に含む、請求項２又は請求項３に記載の方法。
5. 前記第２周波数成分は、前記決定された指標の前記変動に対してｍ次多項式をフィッティングさせることにより決定され、
   第１周波数成分は、前記決定された指標の前記変動内の前記指標から、前記ｍ次多項式の値を引くことにより補正された指標変動を生成し、前記補正された指標変動にｎ次の多項式をフィッティングさせることにより決定され、
   ここで、ｍとｎは整数であって、ｍがｎより小さい、請求項４に記載の方法。
6. 前記補正データが信頼できることを示す信頼性指標が設定されている場合には、前記補正データを使用して、前記ＯＣＴ撮像システムと前記撮像ターゲットの相対運動により生じる、前記Ｂスキャンシーケンス内のＢスキャン間の前記軸方向変位を補償することを更に含む、請求項１～請求項５のいずれか一項に記載の方法。
7. 前記第１周波数成分は、
   （ｉ）各残差値が指標変動における指標とｎ次多項式の対応する値との差として計算される、複数の残差値を計算するステップと、
   （ｉｉ）前記複数の残差値に、第１の正の閾値より高いか第１の負の閾値より低い外れ値が含まれるか否かを判定するステップと、
   （ｉｉｉ）その複数の残差値に外れ値が含まれると判定される場合、外れ値に対応する指標を指標の変動から除去して、更新された指標変動を生成し、また、残差が外れ値を含まないと判定される場合には、ｎ次多項式を第１周波数成分として決定して処理を終了するステップと、
   （ｉｖ）ｎ次多項式を更新された指標変動にフィッティングさせるステップと、
   から成る処理の少なくとも２回の反復を更に実行することにより決定され、
   複数の残差値における各残差値は、処理の第１の反復において、補正された指標変動における指標と、指標の補正された変動にフィッティングされたｎ次多項式の対応する値との間の差として計算され、
   前記複数の残差値における各残差値は、前記処理の残りの１以上の反復の各反復において、前記処理の以前の反復において生成された前記更新された指標変動における指標と、前記処理の以前の反復において生成された前記更新された指標変動にフィッティングされたｎ次多項式の対応する値との間の差として計算される、請求項５に記載の方法。
8. 信頼性指標が、前記補正データが信頼できることを示すように設定されている場合において、
   第２の正の閾値が前記第１の正の閾値より小さく、第２の負の閾値が前記第１の負の閾値より大きい場合、前記複数の残差値の中で、第２の正の閾値よりも大きいか、第２の負の閾値より小さい大きさを有する残差値の数を決定することと、
   前記決定された残差値の数が第３の閾値より小さい場合に、前記Ｂスキャンシーケンス内のＢスキャンに前記第１周波数成分に基づくオフセットを適用することによって、前記Ｂスキャンシーケンス内のＢスキャン間の前記軸方向変位を補償することと、
   前記決定された残差値の数が、第３の閾値より小さくない場合に、前記Ｂスキャンシーケンス内のＢスキャン間の前記軸方向変位を補償しないように決定することと、
   をさらに含む、請求項７に記載の方法。
9. 撮像ターゲットを横断して延在する平行な複数の走査線に沿って走査する光干渉断層計（ＯＣＴ）撮像システムにより取得された前記撮像ターゲットのＢスキャンシーケンスからなるＣスキャンデータと、前記走査線を横切るトレース走査線に沿って走査を実行する前記ＯＣＴ撮像システムによって取得されたトレーススキャンデータとを、前記ＯＣＴ撮像システムによるＢスキャンの取得中にその間の距離を変動させる前記ＯＣＴ撮像システムと前記撮像ターゲットとの相対運動により生じる、前記Ｂスキャンシーケンス内のＢスキャン同士の間の軸方向変位を補償するための補正データを生成し、かつ、さらに前記生成された補正データの信頼性を示す信頼性指標を生成するために処理する方法であって、
   前記Ｂスキャンシーケンスにおける少なくともいくつかの各Ｂスキャンに対し、前記Ｂスキャン及び前記トレーススキャンデータ内の共通の視覚特徴のそれぞれの表示の間の軸方向シフトのそれぞれの指標を決定することにより前記補正データを生成することと、
   信頼性指標を生成することであって、
   前記少なくともいくつかのＢスキャン内のＢスキャンのペアに対して決定された前記指標を使用して、前記Ｂスキャンの複数ペアが取得されたときの前記ＯＣＴ撮像システムに対する前記撮像ターゲットの相対的なそれぞれの速度又はそれぞれの加速度を示すメトリックのそれぞれの値を計算し、
   前記メトリックの計算値の少なくとも所定数が、閾値を超えるか否かを判定し、
   前記メトリックの計算値の少なくとも所定数が前記閾値を超えると判定される場合に、前記補正データが信頼できないことを示す信頼性指標を設定し、
   前記メトリックの計算値の少なくとも所定数が前記閾値を超えないと判定される場合に、前記補正データが信頼できることを示す信頼性指標を設定する、
   ことによって、信頼性指標を生成することと、
   を含む、処理方法。
10. 前記補正データが信頼できることを示す信頼性指標が設定されている場合には、前記補正データを使用して、前記ＯＣＴ撮像システムと前記撮像ターゲットとの間の相対運動により生じる、前記Ｂスキャンシーケンス内のＢスキャン間の前記軸方向変位を補償することを更に含む、請求項９に記載の方法。
11. 前記補正データは、
    （ｉ）各残差値が複数の指標の変動における指標と、ｎ次多項式の対応する値との間の差として計算される、複数の残差値を計算するステップと、
    （ｉｉ）前記複数の残差値に、第１の正の閾値より高いか第１の負の閾値より低い外れ値が含まれるか否かを判定するステップと、
    （ｉｉｉ）前記複数の残差値に外れ値が含まれると判定される場合、前記外れ値に対応する前記指標を前記指標の変動から除去して更新された指標変動を生成し、また、前記残差値が外れ値を含まないと判定される場合には、前記ｎ次多項式を前記補正データとして決定して処理を終了するステップと、
    （ｉｖ）ｎ次多項式を更新された指標変動にフィッティングさせるステップと、
    から成る処理の少なくとも２回の反復を更に実行することにより生成され、
    ここで、複数の残差値における各残差値は、処理の第１の反復において、補正された指標変動における指標と、指標の補正された変動にフィッティングされたｎ次多項式の対応する値との間の差として計算され、
    前記複数の残差値における各残差値は、前記処理の残りの１以上の反復の各反復において、前記処理の以前の反復において生成された前記更新された指標変動における指標と、前記処理の以前の反復において生成された前記更新された指標変動にフィッティングされたｎ次多項式の対応する値との間の差として計算される、請求項９又は請求項１０に記載の方法。
12. 信頼性指標が、前記補正データが信頼できることを示すように設定されている場合において、
    第２の正の閾値が前記第１の正の閾値より小さく、第２の負の閾値が前記第１の負の閾値より大きい場合、前記複数の残差値の中で、第２の正の閾値よりも大きいか、第２の負の閾値より小さい大きさを有する残差値の数を決定することと、
    前記決定された残差値の数が第３の閾値より小さい場合に、前記Ｂスキャンシーケンス内のＢスキャンに前記補正データに基づくオフセットを適用することによって、前記Ｂスキャンシーケンス内のＢスキャン間の前記軸方向変位を補償することと、
    前記決定された残差値の数が、第３の閾値より小さくない場合に、前記Ｂスキャンシーケンス内のＢスキャン間の前記軸方向変位を補償しないように決定することと、
    をさらに含む、請求項１１に記載の方法。
13. プロセッサによって実行されるとき、請求項１～請求項１２の少なくとも１項に記載の方法を前記プロセッサに実行させるコンピュータプログラム命令を含む、コンピュータプログラム。
14. 光干渉断層計（ＯＣＴ）撮像システムにより取得された撮像ターゲットのＢスキャンシーケンスからなるＣスキャンデータを、前記ＯＣＴ撮像システムによりＢスキャンを取得中にその間の距離を変動させる前記ＯＣＴ撮像システムと前記撮像ターゲットとの相対運動により生じる前記Ｂスキャンシーケンス中のＢスキャン同士の間の軸方向変位を補償するための補正データを生成するために、またさらに前記生成された補正データの信頼性を示す信頼性指標を生成するために、処理するように構成されたデータ処理装置であって、
    前記Ｂスキャンシーケンス内の複数の隣接Ｂスキャンペアの内の各隣接Ｂスキャンペアに対し、隣接するＢスキャン内の共通の視覚特徴の各表示間の軸方向シフトのそれぞれの指標を決定することによって前記補正データを決定するように構成された補正データ生成モジュールと、
    信頼性指標生成モジュールであって、
    前記Ｂスキャンシーケンス内のＢスキャンの複数ペアを使用して、前記Ｂスキャンの複数ペアが取得されたときの前記ＯＣＴ撮像システムに対する前記撮像ターゲットの相対的なそれぞれの速度又はそれぞれの加速度を示すメトリックのそれぞれの値を計算し、
    前記メトリックの少なくとも所定数の前記計算値が、閾値を超えるか否かを判定し、
    前記メトリックの少なくとも前記所定数の計算値が前記閾値を超えると判定される場合に、前記信頼性指標を前記補正データが信頼できないことを示すように設定し、
    メトリックの少なくとも所定数の計算値が閾値を超えないと判定される場合に、信頼性指標を、補正データが信頼できることを示すように設定する、
    ことによって、信頼性指標を生成するように構成された信頼性指標生成モジュールと、
    を含む、データ処理装置。
15. 撮像ターゲットを横断して延在する平行な走査線に沿って走査する光干渉断層計（ＯＣＴ）撮像システムにより取得された撮像ターゲットのＢスキャンシーケンスからなるＣスキャンデータと、前記走査線を横切るトレース走査線に沿って走査を実行する前記ＯＣＴ撮像システムによって取得されたトレーススキャンデータとを、前記ＯＣＴ撮像システムによるＢスキャンの取得中にその間の距離を変動させる前記ＯＣＴ撮像システムと前記撮像ターゲットとの相対運動により生じる前記Ｂスキャンシーケンス内のＢスキャン同士の間の軸方向変位を補償するための補正データを生成し、かつ、さらに前記生成された補正データの信頼性を示す信頼性指標を生成するために、処理するように構成されたデータ処理装置であって、
    前記Ｂスキャンシーケンスにおける少なくともいくつかのＢスキャンの各Ｂスキャンに対し、前記Ｂスキャン及び前記トレーススキャンデータ内の共通の視覚特徴のそれぞれの表示の間の軸方向シフトのそれぞれの指標を決定することにより前記補正データを生成するように構成された補正データ生成モジュールと、
    信頼性生成モジュールであって、
    前記少なくともいくつかのＢスキャン内のＢスキャンのペアに対して決定された前記指標を使用して、前記Ｂスキャンの複数ペアが取得されたときの前記ＯＣＴ撮像システムに対する前記撮像ターゲットの相対的なそれぞれの速度又はそれぞれの加速度を示すメトリックの、それぞれの値を計算し、
    そのメトリックの少なくとも所定数の計算値が、閾値を超えるか否かを判定し、
    そのメトリックの少なくとも所定数の計算値が閾値を超えると判定される場合に、信頼性指標を、補正データが信頼できないことを示すように設定し、
    そのメトリックの少なくとも所定数の計算値が閾値を超えないと判定される場合に、信頼性指標を、補正データが信頼できることを示すように設定する、
    ことによって、信頼性指標を生成するように構成された信頼性生成モジュールと、
    を含む、データ処理装置。

Images