ホーム > 特許ランキング > オプトス ピーエルシー
知財求人 - 知財ポータルサイト「IP Force」
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2024-08-07
(45)【発行日】2024-08-16
(54)【発明の名称】ＯＣＴＡデータ処理の方法、コンピュータプログラム、及び装置
(51)【国際特許分類】
A61B 3/10 20060101AFI20240808BHJP
A61B 3/12 20060101ALI20240808BHJP
G06T 1/00 20060101ALI20240808BHJP
【ＦＩ】
A61B3/10 100
A61B3/12 300
G06T1/00 290Z
【請求項の数】
15
【外国語出願】
(21)【出願番号】P 2022212215
(22)【出願日】2022-12-28
(65)【公開番号】P2023145332
(43)【公開日】2023-10-11
【審査請求日】2023-05-23
(31)【優先権主張番号】22164712.6
(32)【優先日】2022-03-28
(33)【優先権主張国・地域又は機関】EP
(73)【特許権者】
【識別番号】509012991
【氏名又は名称】オプトス ピーエルシー
(74)【代理人】
【識別番号】110001519
【氏名又は名称】弁理士法人太陽国際特許事務所
(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】
【図1A】網膜の撮像領域における表層のＯＣＴＡ正面投影の例である。
【図1B】網膜の撮像領域におけるより深い高反射層のＯＣＴＡ正面投影の例である。
【図1C】本例示的実施形態の装置により生成された、網膜の撮像領域における上記のより深い高反射層のＯＣＴＡ正面投影の例である。
【図2】網膜の共通撮像領域のＢスキャンを４回繰り返し取得するときの、網膜を通る光線の伝播の模試図である。
【図3A】ＯＣＴＡデータを処理して正面投影の投影アーチファクトの少ない補正ＯＣＴＡデータを生成する、本明細書中の例示的実施形態の装置の模式図である。
【図3B】ＯＣＴＡデータを取得するように構成された撮像部品を有する撮像装置と、ＯＣＴＡデータを処理して補正ＯＣＴＡデータを生成するための例示的実施形態の装置の模式図である。
【図4】ＯＣＴＡデータ生成のためのＯＣＴＡデータ生成アルゴリズムによるＯＣＴの繰り返しＢスキャン処理の模式図である。
【図5】例示的実施形態の装置を、プログラマブル信号処理装置に実装したハードウェア実装例を示す図である。
【図6】例示的実施形態による、正面投影における投影アーチファクトがＯＣＴＡデータより少ない補正ＯＣＴＡデータ生成のためのＯＣＴＡデータ処理方法を示すフロー図である。
【図7】それぞれの繰り返しＢスキャン内の対応する位置にあって、相関値の計算に使用されるＯＣＴデータ要素の模式図である
【図9】投影アーチファクト補正の強度を決定するスケール因子の異なる値を使用して例示的実施形態の装置によって生成された、補正ＯＣＴＡデータのＯＣＴＡ正面投影の第１の組を示す図である。
【図10】異なるスケール因子の値を使用して例示的実施形態の装置によって生成された、補正ＯＣＴＡデータのＯＣＴＡ正面投影の第２の組を示す図である。
【発明を実施するための形態】
【0013】
以下にＯＣＴＡ正面投影から投影アーチファクトを除去するためのボリューム軸基準処理アルゴリズムについて述べる。上で概要を述べた、既に生成された正面投影を画像処理してアーチファクトを除去する従来の投影アーチファクト除去アルゴリズムに比べて、本明細書に記載の投影アーチファクト除去技術は、正面投影を生成する前にＯＣＴＡボリュームのＡスキャン内の偽信号を除去ないし低減し得るものである。記載の技術は、ＯＣＴＡボリューム内にある関連の軸方向情報を使用して、網膜の解剖学的構造に起因する偽信号成分をＯＣＴＡボリュームのボクセル値から低減若しくは除去し、真のＯＣＴＡ信号をより良好に反映する補正ボクセル値をもたらすことが可能である。最終正面画像からその場限りの方式でアーチファクトを除去して生成する正面投影に比べて、この方式ではアーチファクト除去後の投影正面画像はより自然な外観を表すことが可能である。これは補正軸方向信号の平均化によって、生成された正面画像内の急峻な変化を排除し得るからである。さらには、既に投影されたＯＣＴＡ正面画像ではなく、ＯＣＴＡ画像のアーチファクトの源になるＯＣＴＡ信号を補正するので、アーチファクト除去はより正確になり得る。
【0014】
図３Ａは例示的実施形態による装置１００の模式図であり、ＯＣＴＡデータ２００と関連するＯＣＴＡデータ２５０を処理して、ＯＣＴデータ２００よりも正面投影中の投影アーチファクトが低減された補正ＯＣＴＡデータ３００、すなわちＯＣＴＡデータ２００の正面投影よりも投影アーチファクトが少ないか、より目立たない正面投影となる補正ＯＣＴＡデータ３００を生成するように構成される。
【0015】
装置１００は、データ処理専用のスタンドアロン装置（例えばデスクトップＰＣやラップトップコンピュータ）の形態で提供されてもよいし、本例示的実施形態のように、図３Ｂに示すような撮像装置４００の一部を構成してもよい。撮像装置４００は、網膜（または他の身体部分）の領域を撮像し、その撮像領域における脈管構造分布を表すＯＣＴＡデータ２００を取得するように構成された撮像構成部５００を備える。撮像装置４００は、取得したＯＣＴＡデータ２００を処理してＯＣＴＡデータ２００に対して正面投影の投影アーチファクトの少ない補正ＯＣＴＡデータ３００を生成するように構成された、例示的実施形態の装置１００も有する。
【0016】
ＯＣＴＡデータ２００は、ＯＣＴデータ取得装置（図示せず）によって撮像された身体部分の領域にある脈管構造分布を表すデータ要素列の配列から構成される。本例示的実施形態においては、ＯＣＴＡデータ２００は撮像された身体部分の三次元領域における脈管構造分布を表すデータ要素（本明細書では「ボクセル」とも称する）の列の二次元配列から構成される。ただし、本明細書に記載の技術は、身体部分の撮像された二次元領域（スライス）にある脈管構造分布を表すデータ要素（本明細書においては「ピクセル」とも称する）の列の一次元配列の形をしたＯＣＴＡデータの処理にも適用可能である。各データ要素は後方散乱された（ＯＣＴ）信号の変動性を示唆する個々の値を表し、血流領域と静止組織領域との区別に使用可能である。
【0017】
本例示的実施形態においては身体部分は眼の網膜の一部を含む。ただし、本明細書に記載の技術は、ＯＣＴＡ正面投影画像に投影アーチファクトが観察される傾向にある他の任意の身体部分から取得したＯＣＴＡデータの処理に一般的に適用できることを理解されたい。
【0018】
図４は、ｘ軸に沿って配置されたＥ個のスキャン高さの１つにおけるＯＣＴＡデータ２００のスライスを形成する、ＯＣＴＡデータ要素の列２１２－１～２１２－Ｗの１次元配列２１０の模式表示である。ここで、一次元配列２１０を導出する、繰り返しＢスキャン２５１、２５３、２５５、２５７がＯＣＴデータ取得装置により取得されている。これらの繰り返しＢスキャンはそのスキャン高さにおける撮像網膜の共通二次元領域を覆い、その領域は、眼の撮像された三次元領域内の網膜のほぼ平坦部（すなわち眼の硝子体に隣接する網膜の表面）の面内にあるように設定されたデカルト座標系の第１軸（ｙ軸）と、網膜の深さ方向（すなわち、硝子体から脈絡膜と強膜へ至る網膜の深さに跨る方向）である軸方向（デカルト座標系のｚ軸）とに沿って延在する。
【0019】
図４に示すように、ＯＣＴデータ２５０の４つの繰り返しＢスキャンは同一スキャン高さで異なる時間に撮影され（図４の時間軸ｔ方向の間隔で示すように）、ＯＣＴＡデータ生成アルゴリズム２３０によって処理されてＯＣＴＡデータのスライス２１０が生成される。ＯＣＴＡデータ３００の列の配列２１０の各列２１２－１～２１２－Ｗ内のデータ要素は、当業者には周知の任意の種類のＯＣＴＡデータ生成アルゴリズム２３０を用いて、ＯＣＴデータ取得装置によって取得された身体部分の撮像領域に関する繰り返しＢスキャン２５１、２５３、２５５、２５７の各ＢスキャンにおけるＡスキャン２５２－１～２５２－Ｗ、２５４－１～２５４－Ｗ、２５６－１～２５６－Ｗ、２５８－１～２５８－Ｗの対応する位置にあるＯＣＴデータ２５０のデータ要素から生成可能である。すなわち、列の配列２１０の各列２１２－１～２１２－Ｗの各要素は、網膜の撮像領域の繰り返しＢスキャンの組のそれぞれのＢスキャン２５１、２５３、２５５、２５７の対応する位置にあるＡスキャンのＯＣＴデータ２５０のデータ要素から生成される。一例として、配列２１０の列２１２－ｉのデータ要素ｄＯＣＴＡ（ｉ，ｊ）は、身体部分の撮像領域のそれぞれのＢスキャン２５１、２５３、２５５、２５７のＡスキャン２５２－ｉ、２５４－ｉ、２５６－ｉ、２５８－ｉ内の対応する位置にあるＯＣＴデータ２５０のデータ要素ｄＯＣＴ１（ｉ，ｊ）、ｄＯＣＴ２（ｉ，ｊ）、ｄＯＣＴ３（ｉ，ｊ）、ｄＯＣＴ４（ｉ，ｊ）を基にして、ＯＣＴＡデータ生成アルゴリズム２３０によって生成される。
【0020】
図３Ａを再度参照すると、装置１００はハードウェア構成部１１０を備え、これがＯＣＴデータ２５０のデータ要素とＯＣＴＡデータ２００のデータ要素に以下で説明する技術による処理を施して補正ＯＣＴＡデータ３００を生成するように構成されている。装置１００はさらに出力構成部１２０を備え、これが、補正ＯＣＴＡデータ３００を、例えばハードウェア構成部１１０によって生成され得る補正ＯＣＴＡデータ３００の正面投影を、例えばコンピュータモニタなどのディスプレイ装置上に表示するように構成されている。
【0021】
ハードウェア構成部１１０は、以下に記載する例示的実施形態のＯＣＴデータ２５０のデータ要素及びＯＣＴＡデータ２００のデータ要素を処理する方法を実行するように構成された任意の種類のコンピュータハードウェアで構成され得る。ハードウェア構成部１１０は、本例示的実施形態におけるように、プロセッサと、そのプロセッサによって実行されるときにプロセッサにその方法を実行させる命令を格納するメモリを備え得る。あるいは、ハードウェア構成部１１０は、その方法を実行するように構成された特定用途向け集積回路（ＡＳＩＣ）、フィールドプログラマブルゲートアレイ（ＦＰＧＡ）、又はその他の電子回路であってよい。
【0022】
図５は、本明細書に記載の技術を用いてＯＣＴＡデータ２００及び関連するＯＣＴデータ２５０を処理するように構成され、かつ例示的実施形態のハードウェア構成部１１０及び出力構成部１２０として機能し得るプログラマブル信号処理ハードウェア６００を模式的に示す図である。プログラマブル信号処理装置６００は、ＯＣＴＡデータ２００及びＯＣＴデータ２５０を受信し、生成された補正ＯＣＴＡデータ３００及び／又はその補正されたＯＣＴＡデータ３００をコンピュータ画面などのディスプレイに表示するためのグラフィック表示を（例えば補正ＯＣＴＡデータ３００の正面投影の形で）出力する、通信インタフェース（Ｉ／Ｆ）６１０を備える。通信インタフェース（Ｉ／Ｆ）６１０はこのように出力構成部１２０の実装例を提供する。信号処理装置６００はさらに、プロセッサ（例えば中央処理ユニットＣＰＵ、及び／又はグラフィック処理ユニットＧＰＵ）６２０、作業メモリ６３０（例えばランダムアクセスメモリ）、及びコンピュータ可読命令から成るコンピュータプログラム６４５を格納する命令記憶部６４０を備える。コンピュータ可読命令はプロセッサ６２０によって実行されると、本明細書に記載のハードウェア構成部１１０の機能を含む様々な機能をプロセッサ６２０に実行させる。作業メモリ６３０は、コンピュータプログラム６４５の実行時にプロセッサ６２０によって使用される情報を格納する。命令記憶部６４０は、コンピュータ可読命令が予め格納されたＲＯＭ（例えば、電気的消去可能プログラム可能読取り専用メモリ（ＥＥＰＲＯＭ）またはフラッシュメモリの形態）で構成されてもよい。あるいは、命令記憶部６４０はＲＡＭ又は類似の種類のメモリを備えてもよく、コンピュータプログラム６４５のコンピュータ可読命令を、ＣＤ－ＲＯＭ、ＤＶＤＲＯＭなどの形態の、非一時的なコンピュータ可読記憶媒体６５０などのコンピュータプログラム製品から、又はコンピュータ可読命令を伝送するコンピュータ可読信号６６０から、そこに入力することができる。いずれにせよ、コンピュータプログラム６４５はプロセッサ６２０で実行されるとき、プロセッサ６２０に本明細書に記載の補正ＯＣＴＡデータ３００を生成する方法を実行させる。言い換えると、例示的実施形態の装置１００は、コンピュータプロセッサ６２０とコンピュータ可読命令を格納するメモリ６４０とを含むことができ、コンピュータ可読命令はコンピュータプロセッサ６２０によって実行されると、以下で説明するようにＯＣＴＡデータ２００とＯＣＴデータ２５０を処理して、正面投影における投影アーチファクトをＯＣＴＡデータ２００に比べて低減した補正ＯＣＴＡデータ３００を生成する処理方法をコンピュータプロセッサ６２０に実行させる。
【0023】
ただし、ハードウェア構成部１１０と出力構成部１２０の１つ又は両方は、構成部のそれぞれの機能を実行する、ＡＳＩＣ、ＦＰＧＡ、あるいはその他の専用集積回路で代替的に実装されてもよいことに留意されたい。
【0024】
図６は、本例示的実施形態の装置１００が、ＯＣＴＡデータ２００を処理して正面投影中の投影アーチファクトがそのＯＣＴＡデータ２００に比べて低減された補正ＯＣＴＡデータ３００を生成する方法を示すフロー図である。一般に、ハードウェア構成部１１０が、図に示すステップＳ１０、Ｓ２０、Ｓ３０のシーケンスで定義される処理Ｐを繰り返し実行することによって、ＯＣＴデータ２５０のデータ要素とＯＣＴＡデータ２００のデータ要素を処理する。この時、ステップＳ１０、Ｓ２０、Ｓ３０のシーケンスの毎回の実行において、第１のデータ要素としてのＯＣＴＡデータ２００のデータ要素と第２のデータ要素としてのＯＣＴＡデータ２００のデータ要素の組み合わせを変える。この処理のシーケンスは、補正ＯＣＴＡデータ３００を生成するためにいくつかの異なる方法で反復されてもよい。補正ＯＣＴＡデータ３００の生成のための処理シーケンスの反復方法のいくつかの例を以下に示す。
【0025】
図６に示すように、処理シーケンスにはステップＳ１０が含まれる。これは、ハードウェア構成部１１０が、（ｉ）ＯＣＴデータ（２５０）のデータ要素の第１のシーケンスであって、各Ｂスキャン内のこの第１シーケンスのデータ要素の各位置が、（同じ配列の座標を有するという意味で）ＯＣＴＡデータ（２００）の列の配列内の第１データ要素の位置に対応する第１のシーケンスと、（ｉｉ）ＯＣＴデータ（２５０）のデータ要素の第２のシーケンスであって、各Ｂスキャン内のこの第２シーケンスのデータ要素の各位置がＯＣＴＡデータ（２００）の列の配列内の第２データ要素の位置に対応し（同じ配列の座標を有するという意味で）、この第２データ要素は、ＯＣＴＡデータ（２００）の列が延在するＯＣＴＡデータ（２００）の軸（すなわちｚ軸）方向の第１データ要素よりもさらに深いところにある、第２のシーケンスと、の間の相関度合いを表す相関値を計算する。三次元ＯＣＴＡデータ配列の軸（ｚ軸）方向に配列されたＯＣＴＡデータ２００のデータ要素（ボクセル）は、網膜の対応する点から、網膜の深さ方向に、すなわちその光を受ける表面から網膜の内部に向かって、延在するＯＣＴＡデータを含む。このように、第２データ要素は、ＯＣＴＡデータ２００の列の配列において、列が延在するＯＣＴＡデータ配列の（ｚ軸）方向で第１データ要素より下にあり、その結果、第２データ要素が、ＯＣＴＡデータ２００の列の配列内の第２データ要素よりも上のデータ要素の位置に対応する網膜内の位置からの光の散乱によって生じるデータ成分を含む可能性がある。このデータ成分がＯＣＴＡデータ２００の正面投影にアーチファクトを生じさせる。第２ＯＣＴＡデータ要素は、この例示的実施形態では、第１ＯＣＴＡデータ要素と同じＯＣＴＡデータ２００の列内にあってよい。第１シーケンスのデータ要素は、これらのデータ要素を含む各繰り返しＢスキャンの取得順序と同じ順番で現れる。同様に、第２シーケンスのデータ要素は、これらのデータ要素を含む各繰り返しＢスキャンの取得順序と同じ順番で現れる。
【0026】
図７は、それぞれの繰り返しＢスキャン内に対応して位置し、それを基にハードウェア構成部１１０によって相関値が計算される、ＯＣＴデータ要素を模式的に示す図である。
この図では、データ要素の第１シーケンス（図７に示すＢスキャン１～４内のそのそれぞれの位置は、ＯＣＴＡデータ３００の列の配列２１０内のＯＣＴＡデータ要素ｄＯＣＴＡ（ｉ，ｊ）の位置に対応する）は、ＯＣＴデータ要素ｄＯＣＴ１（ｉ，ｊ）、ｄＯＣＴ２（ｉ，ｊ）、ｄＯＣＴ３（ｉ，ｊ）及びｄＯＣＴ４（ｉ，ｊ）で表されている。さらに、データ要素の第２シーケンス（図７に示すＢスキャン１～４内のそのそれぞれの位置は、配列２１０内のＯＣＴＡデータ要素ｄＯＣＴＡ（ｉ，ｋ）の位置に対応する）は、ＯＣＴデータ要素ｄＯＣＴ１（ｉ，ｋ）、ｄＯＣＴ２（ｉ，ｋ）、ｄＯＣＴ３（ｉ，ｋ）、ｄＯＣＴ４（ｉ，ｋ）で表されている。ここで、ｉ、ｊ、ｋは整数であり、かつｋ＞ｊである。図７の右側に示すように、例示的な相関係数ｒｊｋの形で示す相関値は、ステップＳ１０においてＯＣＴデータ要素の第１シーケンスと第２シーケンスを基にハードウェア構成部１１０によって計算される。この値は、ＯＣＴデータ要素の値が２つのシーケンス間で同じように変動する場合（例示値０．９９９７を取る）は、その変動が似ていない場合（例示値０．００２８を取る）よりも高くなる。
この図では、データ要素の第１シーケンス（図７に示すＢスキャン１～４内のそのそれぞれの位置は、ＯＣＴＡデータ３００の列の配列２１０内のＯＣＴＡデータ要素ｄＯＣＴＡ（ｉ，ｊ）の位置に対応する）は、ＯＣＴデータ要素ｄＯＣＴ１（ｉ，ｊ）、ｄＯＣＴ２（ｉ，ｊ）、ｄＯＣＴ３（ｉ，ｊ）及びｄＯＣＴ４（ｉ，ｊ）で表されている。さらに、データ要素の第２シーケンス（図７に示すＢスキャン１～４内のそのそれぞれの位置は、配列２１０内のＯＣＴＡデータ要素ｄＯＣＴＡ（ｉ，ｋ）の位置に対応する）は、ＯＣＴデータ要素ｄＯＣＴ１（ｉ，ｋ）、ｄＯＣＴ２（ｉ，ｋ）、ｄＯＣＴ３（ｉ，ｋ）、ｄＯＣＴ４（ｉ，ｋ）で表されている。ここで、ｉ、ｊ、ｋは整数であり、かつｋ＞ｊである。図７の右側に示すように、例示的な相関係数ｒｊｋの形で示す相関値は、ステップＳ１０においてＯＣＴデータ要素の第１シーケンスと第２シーケンスを基にハードウェア構成部１１０によって計算される。この値は、ＯＣＴデータ要素の値が２つのシーケンス間で同じように変動する場合（例示値０．９９９７を取る）は、その変動が似ていない場合（例示値０．００２８を取る）よりも高くなる。
【0027】
相関係数は、投影アーチファクトのもととなる偽のＯＣＴＡ信号を補正可能としたまま、真のＯＣＴＡ信号の誤った補正を防ぐためにステップＳ１０で計算される。相関係数（その絶対値）は、ＯＣＴＡ信号の真偽を確率的に識別するものと考えることができる。ボクセルにおける偽信号の度合いは上部網膜層内の繰り返しボクセルの反射率のばらつきに相関があるので、この相関を補正重みとして使用可能である。したがって、相関が強いほど信号が偽信号である可能性が高く、またその逆も成り立つ。
【0028】
ハードウェアモジュール１１０で計算される列内のデータ要素ｊとデータ要素ｋとの間の相関値は、本例示的実施形態のように、繰り返しＢスキャンで記録されたこれらのデータ要素の反射率が２つのランダム変数であると仮定して、列内のデータ要素ｊとｋに対するピアソン相関係数ｒｊｋ（ここで－１≦ｒｊｋ≦１）として計算できる。これとは別に、相関値は、例えばスピアマンのρ又はケンドールのτであってもよい。これらは実の値ではなく値のランクを使用することがピアソンの相関係数ｒとの違いである。ランクの使用では状況によっては微妙な相関を検出できないことがあるので、ピアソンの相関係数ｒが有利であろう。
【0029】
再び図６を参照すると、ステップＳ２０においてハードウェア構成部１１０が第１のデータ要素ｄＯＣＴＡ（ｉ，ｊ）と計算した補正値ｒｊｋを用いて、第２のデータ要素ｄＯＣＴＡ（ｉ，ｋ）に対する補正を計算する。ハードウェア構成部１１０は、本例示的実施形態におけるように、第１変数と第２変数の補正関数を評価することによって第２データ要素ｄＯＣＴＡ（ｉ，ｋ）に対する補正値を計算することができる。この補正関数は、第１データ要素ｄＯＣＴＡ（ｉ，ｊ）の値を第１変数の値とし、計算した補正値ｒｊｋを第２変数の値としたときに、第１変数の値の増加、及び第２変数の値の増加と共に増加する。より具体的には、ハードウェア構成部１１０は本例示的実施形態におけるように、第１データ要素ｄＯＣＴＡ（ｉ，ｊ）と計算した補正値ｒｊｋとスケール因子ｓとの積を評価することによって、第２データ要素に対する補正を計算することができる。このように前述の関数は本例示的実施形態においては第１変数と第２変数の積で構成されるが、これは単なる例示であり、関数の形態はこれに限るものではない。
【0030】
スケール因子は固定数であってもよいし、本例示的実施形態のようにユーザが調整可能であって、装置１００のユーザが装置１００（例えば図５に示すプログラマブル信号処理ハードウェアの通信インタフェース（Ｉ／Ｆ）６１０）に接続されたユーザインタフェース（例えばマウス及び／又はキーボード）を介してスケール因子の値を設定可能であってもよい。スケール因子の調整については以下において図１０、１１を参照してさらに議論する。
【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】
前述したように、補正ＯＣＴデータ３００を生成するために、ステップＳ１０～Ｓ３０のシーケンスの各実行において、第１データ要素としてのＯＣＴＡデータ（２００）のデータ要素と第２データ要素としてのＯＣＴＡデータ（２００）のデータ要素との異なる組み合わせを用いてステップＳ１０～Ｓ３０のシーケンスが複数回実行される。言い換えれば、ステップＳ１０～Ｓ３０の各実行において第１データ要素及び第２データ要素として使用されるＯＣＴＡデータのデータ要素の組み合わせは、ステップＳ１０～Ｓ３０の他のすべての実行において第１データ要素と第２データ要素として使用されるＯＣＴＡデータのデータ要素の組み合わせとは異なる。
【0060】
図８Ａと図８Ｂの処理の実行におけるステップＳ２０の各実行において、ハードウェア構成部１１０が、第１データ要素ｄＯＣＴＡ（ｉ，ｊ）と、計算した補正値ｒｊｋとスケール因子ｓとの積として、第２データ要素ｄＯＣＴＡ（ｉ，ｋ）に対する補正を計算する。ここで、第１データ要素ｄＯＣＴＡ（ｉ，ｊ）は、処理中の列内のそれより上のＯＣＴＡデータ要素によってｄＯＣＴＡ（ｉ，ｊ）に加えられた以前の補正を反映する。他方、第２データ要素ｄＯＣＴＡ（ｉ，ｋ）は、ステップＳ２０の他の実行中に、ｄＯＣＴＡ（ｉ，ｊ）のみならず、列中のｄＯＣＴＡ（ｉ，ｊ）より下でｄＯＣＴＡ（ｉ，ｋ）より上のＯＣＴＡデータ要素を基にして補正される。
【0061】
上で述べたように、ステップＳ２０において第２データ要素ｄＯＣＴＡ（ｉ，ｋ）に対する補正の計算に使用されるスケール因子ｓは、ユーザによって調整可能である。図９と図１０は、本例示的実施形態の装置１００によって生成された、補正ＯＣＴＡデータ３００のＯＣＴＡ正面投影の例を示す。図９と図１０のグリッドは、図８Ａ、８ＢのＰＡＲアルゴリズムを２つの各ＯＣＴＡボリュームに適用した結果であり、各グリッドセルは、ＰＡＲアルゴリズムにおいて図に示したスケール因子ｓの値を使用して得られる正面投影を表す。図９に示すように、小さい値のｓ（例えば０．００５と０．０１５）は、投影アーチファクトの十分な抑制を可能とせず、これらは正面投影画像中に「ゴースト」の脈管構造として残る。他方で、比較的高いｓの値は、投影アーチファクトが効果的に抑制された結果をもたらすが、その代償として撮像された（真の）脈管構造の品質が劣化する。図９の例において、ｓの値が０．０７５付近で、一方での効果的な投影アーチファクトの抑制と、他方での脈管構造の画像品質維持との間の良好なバランスが達成される。同様な傾向が図１０において観察される。ここでは、約０．０７５のスケール因子の値が同様に前述の競合する要求間の良好なバランスを達成する。
【0062】
これまでに述べた実施形態のいくつかを、以下の実施例Ｅ１～Ｅ１６にまとめる。
【0063】
Ｅ１． 身体部分の撮像領域の脈管構造分布を表す、データ要素列の配列から成る光干渉断層撮影アンギオグラフィ（ＯＣＴＡ）データ（２００）を処理して、正面投影内の投影アーチファクトがＯＣＴＡデータ（２００）に対して低減された補正ＯＣＴＡデータ（３００）を生成するように構成された装置（１００）
ＯＣＴＡデータ（２００）の列の配列の各列にあるデータ要素は、身体部分の撮像領域の繰り返しＢスキャンの組の中の各ＢスキャンのＡスキャン内の対応位置にあるＯＣＴデータ（２５０）のデータ要素から生成されたものであり、この装置は、
ハードウェア構成部（１１０）であって、
（ｉ）ＯＣＴデータ（２５０）のデータ要素の第１シーケンスであって、各Ｂスキャン内の第１シーケンスのデータ要素の各位置がＯＣＴＡデータ（２００）の列の配列内の第１データ要素の位置に対応する第１のシーケンスと、ＯＣＴデータ（２５０）のデータ要素の第２のシーケンスであって、各Ｂスキャン内の第２シーケンスのデータ要素の各位置がＯＣＴＡデータ（２００）の列の配列内の第２データ要素の位置に対応し、この第２データ要素はＯＣＴＡデータ（２００）の列が延在するＯＣＴＡデータ（２００）の軸方向に第１データ要素よりもさらに深部にある第２のシーケンスと、の間の相関度合いを表す相関値を計算するステップと、
（ｉｉ）第１のデータ要素と計算した相関値とを用いて、第２データ要素に対する補正を計算するステップと、
（ｉｉｉ）計算した補正を第２データ要素に適用するステップと、
によって、ＯＣＴデータ（２５０）のデータ要素とＯＣＴＡデータ（２００）のデータ要素とを処理することを含む方法を実行するように構成され、補正ＯＣＴＡデータ（３００）を生成するために、ステップ（ｉ）～（ｉｉｉ）の各実行ごとに第１データ要素としてのＯＣＴＡデータ（２００）のデータ要素と第２データ要素としてのＯＣＴＡデータ（２００）のデータ要素との異なる組み合わせを用いて、ステップ（ｉ）～（ｉｉｉ）を複数回実行するように構成されたハードウェア構成部（１１０）と、
補正ＯＣＴＡデータ（３００）を出力するように構成された出力構成部（１２０）と、
を備える。
ＯＣＴＡデータ（２００）の列の配列の各列にあるデータ要素は、身体部分の撮像領域の繰り返しＢスキャンの組の中の各ＢスキャンのＡスキャン内の対応位置にあるＯＣＴデータ（２５０）のデータ要素から生成されたものであり、この装置は、
ハードウェア構成部（１１０）であって、
（ｉ）ＯＣＴデータ（２５０）のデータ要素の第１シーケンスであって、各Ｂスキャン内の第１シーケンスのデータ要素の各位置がＯＣＴＡデータ（２００）の列の配列内の第１データ要素の位置に対応する第１のシーケンスと、ＯＣＴデータ（２５０）のデータ要素の第２のシーケンスであって、各Ｂスキャン内の第２シーケンスのデータ要素の各位置がＯＣＴＡデータ（２００）の列の配列内の第２データ要素の位置に対応し、この第２データ要素はＯＣＴＡデータ（２００）の列が延在するＯＣＴＡデータ（２００）の軸方向に第１データ要素よりもさらに深部にある第２のシーケンスと、の間の相関度合いを表す相関値を計算するステップと、
（ｉｉ）第１のデータ要素と計算した相関値とを用いて、第２データ要素に対する補正を計算するステップと、
（ｉｉｉ）計算した補正を第２データ要素に適用するステップと、
によって、ＯＣＴデータ（２５０）のデータ要素とＯＣＴＡデータ（２００）のデータ要素とを処理することを含む方法を実行するように構成され、補正ＯＣＴＡデータ（３００）を生成するために、ステップ（ｉ）～（ｉｉｉ）の各実行ごとに第１データ要素としてのＯＣＴＡデータ（２００）のデータ要素と第２データ要素としてのＯＣＴＡデータ（２００）のデータ要素との異なる組み合わせを用いて、ステップ（ｉ）～（ｉｉｉ）を複数回実行するように構成されたハードウェア構成部（１１０）と、
補正ＯＣＴＡデータ（３００）を出力するように構成された出力構成部（１２０）と、
を備える。
【0064】
Ｅ２． Ｅ１による装置（１００）
上記の方法が、ＯＣＴＡデータ（２００）の各列のデータ要素を処理するステップ（ｉ）～（ｉｉｉ）を複数回実行するその各回において、ステップ（ｉ）の一部としてかつ相関値を計算する前に、
ＯＣＴＡデータ（２００）の列内に開始データ要素から列の軸方向に沿って順番に配置された複数のデータ要素の各データ要素を、そのデータ要素がＯＣＴＡ信号閾値を超えるか否かを判定して、１つのデータ要素がＯＣＴＡ信号閾値を超えると判定されるまで、処理し、
ＯＣＴＡ信号閾値を超えると判定されたデータ要素を第１のデータ要素に設定する、
ことをさらに含む。
上記の方法が、ＯＣＴＡデータ（２００）の各列のデータ要素を処理するステップ（ｉ）～（ｉｉｉ）を複数回実行するその各回において、ステップ（ｉ）の一部としてかつ相関値を計算する前に、
ＯＣＴＡデータ（２００）の列内に開始データ要素から列の軸方向に沿って順番に配置された複数のデータ要素の各データ要素を、そのデータ要素がＯＣＴＡ信号閾値を超えるか否かを判定して、１つのデータ要素がＯＣＴＡ信号閾値を超えると判定されるまで、処理し、
ＯＣＴＡ信号閾値を超えると判定されたデータ要素を第１のデータ要素に設定する、
ことをさらに含む。
【0065】
Ｅ３． Ｅ２による装置（１００）
ステップ（ｉ）～（ｉｉｉ）の複数回の実行が、
最初にステップ（ｉ）を実行し、次いでステップ（ｉ）を繰り返して、ステップ（ｉ）の実行における第１のデータ要素として、ステップ（ｉ）の最初の回と各繰り返しとで、ＯＣＴＡデータ（２００）の列内の異なる各データ要素を設定し、
ＯＣＴＡデータ（２００）の列内の複数のデータ要素の内の１つのデータ要素を第１データ要素として設定するためのステップ（ｉ）の最初の実行においては、ハードウェア構成部（１１０）はＯＣＴＡデータ（２００）の列内の所定のデータ要素を開始データ要素として使用するように構成され、
さらに、ステップ（ｉ）の繰り返し実行のそれぞれにおいてハードウェア構成部（１１０）は、先行するステップ（ｉ）の実行において第１データ要素として設定されたＯＣＴＡデータ２００の各データ要素を、開始データ要素として使用するように構成される。
ステップ（ｉ）～（ｉｉｉ）の複数回の実行が、
最初にステップ（ｉ）を実行し、次いでステップ（ｉ）を繰り返して、ステップ（ｉ）の実行における第１のデータ要素として、ステップ（ｉ）の最初の回と各繰り返しとで、ＯＣＴＡデータ（２００）の列内の異なる各データ要素を設定し、
ＯＣＴＡデータ（２００）の列内の複数のデータ要素の内の１つのデータ要素を第１データ要素として設定するためのステップ（ｉ）の最初の実行においては、ハードウェア構成部（１１０）はＯＣＴＡデータ（２００）の列内の所定のデータ要素を開始データ要素として使用するように構成され、
さらに、ステップ（ｉ）の繰り返し実行のそれぞれにおいてハードウェア構成部（１１０）は、先行するステップ（ｉ）の実行において第１データ要素として設定されたＯＣＴＡデータ２００の各データ要素を、開始データ要素として使用するように構成される。
【0066】
Ｅ４． Ｅ３による装置（１００）
所定のデータ要素は、ＯＣＴＡデータ（２００）の列の端のデータ要素、又は
脈管構造を含む身体部分の第１の部分を表す列のデータ要素の組のデータ要素であって、そのデータ要素の組と、脈管構造を含まない身体部分の第２の部分を表す列内のデータ要素との間の境界に隣接するデータ要素、
のいずれかである。
所定のデータ要素は、ＯＣＴＡデータ（２００）の列の端のデータ要素、又は
脈管構造を含む身体部分の第１の部分を表す列のデータ要素の組のデータ要素であって、そのデータ要素の組と、脈管構造を含まない身体部分の第２の部分を表す列内のデータ要素との間の境界に隣接するデータ要素、
のいずれかである。
【0067】
Ｅ５． Ｅ１～Ｅ４のいずれかによる装置（１００）
上記の方法が、ＯＣＴＡデータ（２００）の各列のデータ要素を処理するステップ（ｉ）～（ｉｉｉ）を複数回実行するその各回において、ステップ（ｉ）の一部としてかつ相関値を計算する前に、
ＯＣＴＡデータ（２００）の列内に初期データ要素から軸方向に順番に配置された複数のデータ要素の各データ要素について、データ要素メトリック、これは繰り返しＢスキャンの各Ｂスキャン内の位置がＯＣＴＡデータ（２００）の列の配列内のデータ要素の位置に対応するＯＣＴデータ要素の少なくともいくつかに基づくものであるが、そのデータ要素メトリックがＯＣＴ値閾値を超えるか否かの判定を、１つのデータ要素におけるデータ要素メトリックがＯＣＴ値閾値を超えると判定されるまで行い、かつデータ要素メトリックがＯＣＴ値閾値を超えると判定されたデータ要素を第２のデータ要素として設定する処理をさらに含む。
上記の方法が、ＯＣＴＡデータ（２００）の各列のデータ要素を処理するステップ（ｉ）～（ｉｉｉ）を複数回実行するその各回において、ステップ（ｉ）の一部としてかつ相関値を計算する前に、
ＯＣＴＡデータ（２００）の列内に初期データ要素から軸方向に順番に配置された複数のデータ要素の各データ要素について、データ要素メトリック、これは繰り返しＢスキャンの各Ｂスキャン内の位置がＯＣＴＡデータ（２００）の列の配列内のデータ要素の位置に対応するＯＣＴデータ要素の少なくともいくつかに基づくものであるが、そのデータ要素メトリックがＯＣＴ値閾値を超えるか否かの判定を、１つのデータ要素におけるデータ要素メトリックがＯＣＴ値閾値を超えると判定されるまで行い、かつデータ要素メトリックがＯＣＴ値閾値を超えると判定されたデータ要素を第２のデータ要素として設定する処理をさらに含む。
【0068】
Ｅ６． Ｅ５に記載の装置（１００）
ステップ（ｉ）～（ｉｉｉ）の複数回の実行が、最初にステップ（ｉ）を実行し、次いでステップ（ｉ）を繰り返して、ＯＣＴＡデータ（２００）の列内の異なる各データ要素をステップ（ｉ）の最初の実行とステップ（ｉ）の各繰り返しとにおける第２データ要素として設定することを含み、
ＯＣＴＡデータ（２００）の列内のデータ要素を第２データ要素として設定するステップＳ１０の最初の実行において、ハードウェア構成部（１１０）は、ＯＣＴＡデータ（２００）の列内のその位置が第１データ要素を含むＯＣＴＡデータ（２００）の列内の第１データ要素の位置に隣接する位置に対応するＯＣＴＡデータ（２００）のデータ要素を初期データ要素として使用するように構成され、
ステップ（ｉ）の各繰り返しにおいてハードウェア構成部（１１０）は、先行するステップ（ｉ）の実行において第２データ要素として設定されていたＯＣＴＡデータ（２００）の各データ要素を初期データ要素として使用するように構成される。
ステップ（ｉ）～（ｉｉｉ）の複数回の実行が、最初にステップ（ｉ）を実行し、次いでステップ（ｉ）を繰り返して、ＯＣＴＡデータ（２００）の列内の異なる各データ要素をステップ（ｉ）の最初の実行とステップ（ｉ）の各繰り返しとにおける第２データ要素として設定することを含み、
ＯＣＴＡデータ（２００）の列内のデータ要素を第２データ要素として設定するステップＳ１０の最初の実行において、ハードウェア構成部（１１０）は、ＯＣＴＡデータ（２００）の列内のその位置が第１データ要素を含むＯＣＴＡデータ（２００）の列内の第１データ要素の位置に隣接する位置に対応するＯＣＴＡデータ（２００）のデータ要素を初期データ要素として使用するように構成され、
ステップ（ｉ）の各繰り返しにおいてハードウェア構成部（１１０）は、先行するステップ（ｉ）の実行において第２データ要素として設定されていたＯＣＴＡデータ（２００）の各データ要素を初期データ要素として使用するように構成される。
【0069】
Ｅ８． Ｅ１による装置（１００）
ハードウェア構成部（１１０）が、データ要素の複数の列の内の各列を、第１データ要素及び第２データ要素として列内の異なるデータ要素の組み合わせを用いてステップ（ｉ）～（ｉｉｉ）を実行することにより個別に処理するように構成される。この処理は、以下のステップによる。
（ａ）ＯＣＴＡデータ（２００）の列内に開始データ要素から列の第１方向に沿って順番に配置された複数のデータ要素の各データ要素について、データ要素がＯＣＴＡ信号閾値を超えるか否かの判定を、１つのデータ要素がＯＣＴＡ信号閾値を超えると判定されるまで行い、ＯＣＴＡ信号閾値を超えると判定されたデータ要素を第１のデータ要素に設定するステップ。
（ｂ）ＯＣＴＡデータ（２００）の列内に初期データ要素から列の第１方向に順番に配置された１以上のデータ要素の各データ要素について、繰り返しＢスキャンの各Ｂスキャン内の位置がＯＣＴＡデータ（２００）の列の配列内のデータ要素の位置に対応するＯＣＴデータ要素の少なくともいくつかに基づくデータ要素メトリックがＯＣＴ値閾値を超えるか否かの判定を、１つのデータ要素のデータ要素メトリックがＯＣＴＡ値閾値を超えると判定されるまで行い、データ要素メトリックがＯＣＴ値閾値を超えると判定されたデータ要素を第２データ要素として設定する、ことによって処理するステップ。
（ｃ）ＯＣＴデータ（２５０）のデータ要素の第１シーケンスであって、
それぞれのＢスキャン内の第１シーケンスのデータ要素のそれぞれの位置がＯＣＴＡデータ（２００）の列の配列内の第１データ要素の位置に対応する、第１シーケンスと、
ＯＣＴデータ（２５０）のデータ要素の第２シーケンスであって、それぞれのＢスキャン内の第２シーケンスのデータ要素のそれぞれの位置がＯＣＴＡデータ（２００）の列の配列内の第２データ要素の位置に対応する、第２のシーケンスと、
の間の相関の度合いを示す相関値を計算するステップ。
（ｄ）第１データと計算した相関値を用いて第２データ要素に対する補正を計算するステップ。
（ｅ）計算した補正を第２データ要素に適用するステップ。
（ｆ）ステップ（ｂ）～（ｅ）の反復ステップであって、
ステップ（ｂ）の最初の実行において、初期データ要素は、列における第１データ要素の位置に隣接する位置にあるＯＣＴＡデータ（２００）のデータ要素であり、
各繰り返しステップ（ｂ）の初期データ要素は、ステップ（ｂ）の先行する実行における第２データ要素として設定された列内のデータ要素であり、
ステップ（ｂ）～（ｅ）は、ＯＣＴＡデータ（２００）の列にステップ（ｂ）で処理すべきデータ要素がそれ以上なくなるまで繰り返される、反復ステップ。
（ｇ）ステップ（ａ）～（ｆ）の反復ステップであって、
ステップ（ａ）の最初の実行において、前記開始データ要素はＯＣＴＡデータ（２００）の列内の所定のデータ要素であり、
各繰り返しステップ（ａ）の開始データ要素は、ステップ（ａ）の先行する実行における第１データ要素として設定された列内のデータ要素であり、
ステップ（ａ）～（ｆ）は、ステップ（ａ）で処理すべきＯＣＴＡデータ（２００）の列内のデータ要素が２つ以下となるまで繰り返される、反復ステップ。
ハードウェア構成部（１１０）が、データ要素の複数の列の内の各列を、第１データ要素及び第２データ要素として列内の異なるデータ要素の組み合わせを用いてステップ（ｉ）～（ｉｉｉ）を実行することにより個別に処理するように構成される。この処理は、以下のステップによる。
（ａ）ＯＣＴＡデータ（２００）の列内に開始データ要素から列の第１方向に沿って順番に配置された複数のデータ要素の各データ要素について、データ要素がＯＣＴＡ信号閾値を超えるか否かの判定を、１つのデータ要素がＯＣＴＡ信号閾値を超えると判定されるまで行い、ＯＣＴＡ信号閾値を超えると判定されたデータ要素を第１のデータ要素に設定するステップ。
（ｂ）ＯＣＴＡデータ（２００）の列内に初期データ要素から列の第１方向に順番に配置された１以上のデータ要素の各データ要素について、繰り返しＢスキャンの各Ｂスキャン内の位置がＯＣＴＡデータ（２００）の列の配列内のデータ要素の位置に対応するＯＣＴデータ要素の少なくともいくつかに基づくデータ要素メトリックがＯＣＴ値閾値を超えるか否かの判定を、１つのデータ要素のデータ要素メトリックがＯＣＴＡ値閾値を超えると判定されるまで行い、データ要素メトリックがＯＣＴ値閾値を超えると判定されたデータ要素を第２データ要素として設定する、ことによって処理するステップ。
（ｃ）ＯＣＴデータ（２５０）のデータ要素の第１シーケンスであって、
それぞれのＢスキャン内の第１シーケンスのデータ要素のそれぞれの位置がＯＣＴＡデータ（２００）の列の配列内の第１データ要素の位置に対応する、第１シーケンスと、
ＯＣＴデータ（２５０）のデータ要素の第２シーケンスであって、それぞれのＢスキャン内の第２シーケンスのデータ要素のそれぞれの位置がＯＣＴＡデータ（２００）の列の配列内の第２データ要素の位置に対応する、第２のシーケンスと、
の間の相関の度合いを示す相関値を計算するステップ。
（ｄ）第１データと計算した相関値を用いて第２データ要素に対する補正を計算するステップ。
（ｅ）計算した補正を第２データ要素に適用するステップ。
（ｆ）ステップ（ｂ）～（ｅ）の反復ステップであって、
ステップ（ｂ）の最初の実行において、初期データ要素は、列における第１データ要素の位置に隣接する位置にあるＯＣＴＡデータ（２００）のデータ要素であり、
各繰り返しステップ（ｂ）の初期データ要素は、ステップ（ｂ）の先行する実行における第２データ要素として設定された列内のデータ要素であり、
ステップ（ｂ）～（ｅ）は、ＯＣＴＡデータ（２００）の列にステップ（ｂ）で処理すべきデータ要素がそれ以上なくなるまで繰り返される、反復ステップ。
（ｇ）ステップ（ａ）～（ｆ）の反復ステップであって、
ステップ（ａ）の最初の実行において、前記開始データ要素はＯＣＴＡデータ（２００）の列内の所定のデータ要素であり、
各繰り返しステップ（ａ）の開始データ要素は、ステップ（ａ）の先行する実行における第１データ要素として設定された列内のデータ要素であり、
ステップ（ａ）～（ｆ）は、ステップ（ａ）で処理すべきＯＣＴＡデータ（２００）の列内のデータ要素が２つ以下となるまで繰り返される、反復ステップ。
【0070】
Ｅ９． Ｅ８による装置（１００）
所定のデータ要素は、ＯＣＴＡデータ（２００）の列の端のデータ要素、又は脈管構造を含む身体部分の第１の部分を表す列のデータ要素の組のデータ要素であって、そのデータ要素の組と、脈管構造を含まない身体部分の第２の部分を表す列内のデータ要素との間の境界に隣接するデータ要素、のいずれかである。
所定のデータ要素は、ＯＣＴＡデータ（２００）の列の端のデータ要素、又は脈管構造を含む身体部分の第１の部分を表す列のデータ要素の組のデータ要素であって、そのデータ要素の組と、脈管構造を含まない身体部分の第２の部分を表す列内のデータ要素との間の境界に隣接するデータ要素、のいずれかである。
【0071】
Ｅ１０． Ｅ８又はＥ９に記載の装置（１００）
ハードウェア構成部（１１０）はデータ要素の複数の列を、ステップ（ａ）～（ｇ）によって並列に処理するように構成される。
ハードウェア構成部（１１０）はデータ要素の複数の列を、ステップ（ａ）～（ｇ）によって並列に処理するように構成される。
【0072】
Ｅ１１． Ｅ８～Ｅ１０のいずれかによる装置（１００）
データ要素メトリックは、繰り返しＢスキャンの各Ｂスキャンにおける位置がＯＣＴＡデータ（２００）の列の配列におけるデータ要素の位置に対応する少なくともいくつかのＯＣＴデータ要素の平均と、繰り返しＢスキャンの各Ｂスキャンにおける位置がＯＣＴＡデータ（２００）の列の配列におけるデータ要素の位置に対応する少なくともいくつかのＯＣＴデータ要素の中央値との内の１つである。
データ要素メトリックは、繰り返しＢスキャンの各Ｂスキャンにおける位置がＯＣＴＡデータ（２００）の列の配列におけるデータ要素の位置に対応する少なくともいくつかのＯＣＴデータ要素の平均と、繰り返しＢスキャンの各Ｂスキャンにおける位置がＯＣＴＡデータ（２００）の列の配列におけるデータ要素の位置に対応する少なくともいくつかのＯＣＴデータ要素の中央値との内の１つである。
【0073】
Ｅ１２． Ｅ１～Ｅ１１のいずれかに記載の装置（１００）
相関値は、ピアソンの相関係数値、スピアマンのランク相関係数値、及びケンドールのランク相関係数値の内の１つである。
相関値は、ピアソンの相関係数値、スピアマンのランク相関係数値、及びケンドールのランク相関係数値の内の１つである。
【0074】
Ｅ１３． Ｅ１～Ｅ１２のいずれかによる装置（１００））
ハードウェア構成部１１０が、第１データ要素と計算した相関値とスケール因子との積を評価することによって、第２データ要素に対する補正を計算するように構成される。
ハードウェア構成部１１０が、第１データ要素と計算した相関値とスケール因子との積を評価することによって、第２データ要素に対する補正を計算するように構成される。
【0075】
Ｅ１４． Ｅ１～Ｅ１３のいずれかに記載の装置（１００）
ハードウェア構成部（１１０）がさらに、補正ＯＣＴＡデータ（３００）の正面投影を生成するように構成され、出力構成部（１２０）が生成した正面投影を出力するように構成される。
ハードウェア構成部（１１０）がさらに、補正ＯＣＴＡデータ（３００）の正面投影を生成するように構成され、出力構成部（１２０）が生成した正面投影を出力するように構成される。
【0076】
Ｅ１５． Ｅ１４による装置（１００）
出力構成部（１２０）が生成した正面投影をディスプレイ装置上に表示するように構成される。
出力構成部（１２０）が生成した正面投影をディスプレイ装置上に表示するように構成される。
【0077】
Ｅ１６． 撮像装置（４００）
身体部分の領域を撮像し、その身体部分の撮像領域の脈管構造分布を表す光干渉断層撮影アンギオグラフィ、ＯＣＴＡ、データ（２００）を取得するように構成された撮像構成部（５００）と、取得したＯＣＴＡデータ（２００）を処理して、正面投影内の投影アーチファクトがＯＣＴＡデータ（２００）に対して減少した補正ＯＣＴＡデータ（３００）を生成するように構成された、Ｅ１～Ｅ１５のいずれかに記載の装置（１００）と、を備える。
身体部分の領域を撮像し、その身体部分の撮像領域の脈管構造分布を表す光干渉断層撮影アンギオグラフィ、ＯＣＴＡ、データ（２００）を取得するように構成された撮像構成部（５００）と、取得したＯＣＴＡデータ（２００）を処理して、正面投影内の投影アーチファクトがＯＣＴＡデータ（２００）に対して減少した補正ＯＣＴＡデータ（３００）を生成するように構成された、Ｅ１～Ｅ１５のいずれかに記載の装置（１００）と、を備える。
【0078】
ＯＣＴＡデータ２００の各列は、上述の例示的実施形態のように別々に処理して、１）マイクロプロセッサ及び／又はグラフィック処理ユニット（ＧＰＵ）で複数のスレッドやコアを持つ計算リソースが利用できる場合、並列化を可能とし、２）処理中の列に対する近接する列内のボクセルによる影響を無視して補正方式を単純化する、ことが可能である。物理的な観点から実際のモデルは、近接する列に広がる網膜上皮層内の血管からの散乱を考慮することが可能である。他の例示的実施形態では、列ａ内で光源散乱特性を満たすすべてのボクセルｖに対して、ボクセルｖの下のａに近接するすべての列内の全標的ボクセルへの補正が適用可能である。ここで、近接度合いは、ｖに頂点又は頂端を有し、ａに平行な軸を有する、ａに関して定義される影響円錐内に列ｂがあるか否かによって定義される。ユーザは、幾何学的パラメータ（直円錐であるとして、高さ及び底面半径など）を用いて円錐形状を定義することができる。本明細書に記載の例示的実施形態では、影響円錐は実質的に同じ列ａの幅だけに限定され、ａ内のｖの下のボクセルのみが補正される。
【0079】
上記の説明においては、例示的態様をいくつかの例示的実施形態を参照して説明した。したがって本明細書は、制限的ではなく、例示的であるとみなされるべきである。同様に、例示的実施形態の機能及び利点を強調する図面に示す形状は、例示のみを目的として提示される。例示的実施形態のアーキテクチャは十分に柔軟かつ設定可能であり、添付図に示すもの以外の方法で利用可能である。
【0080】
例えば、図８Ａを参照して前述したステップＳ１２０において、ハードウェア構成部１１０は、列内の第１データ要素の軸（ｚ軸）方向の隣のＯＣＴＡデータ２００のデータ要素に初期データ要素を設定する。
次に図８ＡのステップＳ１３０において、ハードウェア構成部１１０は、ＯＣＴＡデータ２００の列内に初期データ要素から軸（ｚ）方向に順番に配置された１以上のデータ要素の各データ要素を処理する。ただし、ステップＳ１３０で処理される、列内の第１データ要素の直下のデータ要素の順番は、軸方向に沿うことに限るわけではなく、替わりにその反対方向であってもよい。すなわち、ハードウェア構成部１１０は代替的に、初期データ要素を、第１データ要素から最も遠い、列の下端にあるＯＣＴＡデータ２００のデータ要素に設定してもよい。そうして、図８ＡのステップＳ１３０の変形において、ハードウェア構成部１１０は、ＯＣＴＡデータ２００の列内に初期データ要素から第１データ要素に向かって軸（ｚ軸）の反対方向に順番に配置された１以上のデータ要素の各データ要素を処理してもよい。こうして、列内のデータ要素が、図８Ａ及び８ＢのステップＳ１２０～ステップＳ１５０において、（例示的実施形態のように）列を下るのではなく列を上って処理することができる。
次に図８ＡのステップＳ１３０において、ハードウェア構成部１１０は、ＯＣＴＡデータ２００の列内に初期データ要素から軸（ｚ）方向に順番に配置された１以上のデータ要素の各データ要素を処理する。ただし、ステップＳ１３０で処理される、列内の第１データ要素の直下のデータ要素の順番は、軸方向に沿うことに限るわけではなく、替わりにその反対方向であってもよい。すなわち、ハードウェア構成部１１０は代替的に、初期データ要素を、第１データ要素から最も遠い、列の下端にあるＯＣＴＡデータ２００のデータ要素に設定してもよい。そうして、図８ＡのステップＳ１３０の変形において、ハードウェア構成部１１０は、ＯＣＴＡデータ２００の列内に初期データ要素から第１データ要素に向かって軸（ｚ軸）の反対方向に順番に配置された１以上のデータ要素の各データ要素を処理してもよい。こうして、列内のデータ要素が、図８Ａ及び８ＢのステップＳ１２０～ステップＳ１５０において、（例示的実施形態のように）列を下るのではなく列を上って処理することができる。
【0081】
本明細書に提示した実施例のソフトウェア実施形態は、命令又は命令のシーケンスを有する一つ又は複数のプログラムなどのコンピュータプログラム又はソフトウェアとして提供可能である。これらは一実施形態においては、機械でアクセス可能若しくは機械で読み取り可能な媒体、命令記憶部、又はコンピュータ可読記憶装置などの製品に含まれるか格納され、このそれぞれが非一時的であり得る。非一時的な機械アクセス可能媒体、機械可読媒体、命令記憶部、又はコンピュータ可読記憶装置上のプログラム又は命令は、コンピュータシステム又は他の電子デバイスをプログラムするのに使用可能である。機械可読媒体又はコンピュータ可読媒体、命令記憶部、及び記憶装置は、これに限らないが、フロッピディスク（登録商標）、光ディスク、及び光磁気ディスク、あるいは、電子的命令を格納若しくは送信するのに適した他の種類の媒体／機械可読媒体／命令記憶部／記憶装置が含まれ得る。本明細書に記載の技術はいかなる特定のソフトウェア構成にも限定されない。これらは任意のコンピューティング環境又は処理環境への適用性を見出し得る。本明細書に使用されている「コンピュータ可読」、「機械アクセス可能媒体」、「機械可読媒体」、「命令記憶部」及び「コンピュータ可読記憶装置」という用語は、機械、コンピュータ又はコンピュータプロセッサによって実行するために命令又は命令シーケンスを格納、コード化、又は送信することが可能で、かつ本明細書に記載の任意の方法を、機械／コンピュータ／コンピュータプロセッサに実行させる任意の媒体を含むものとする。さらに、当分野においてソフトウェアについては、ある形態又は別の形態（例えばプログラム、手順、プロセス、アプリケーション、モジュール、ユニット、ロジックなど）で、動作を起こし、結果をもたらすものを指すのが一般的である。そのような表現は、処理システムによるソフトウェアの実行が、プロセッサに動作を実行させて結果を生成させることの簡略的な言い方に過ぎない。
【0082】
いくつかの実施形態は、特定用途向け集積回路、フィールドプログラマブルゲートアレイを用意することによって、又は従来のコンポーネント回路の適切なネットワークを相互接続することによっても、実装することができる。
【0083】
いくつかの実施形態にはコンピュータプログラム製品が含まれる。コンピュータプログラム製品は、命令が格納された記憶媒体、命令記憶部、又は記憶装置であってよい。それを使用して、コンピュータ又はコンピュータプロセッサが本明細書に記載の例示的実施形態の任意の手順を実行することを制御し、又は実行させることができる。記憶媒体／命令記憶部／記憶装置には、これに限るものではないが一例として、光ディスク、ＲＯＭ、ＲＡＭ、ＥＰＲＯＭ、ＥＥＰＲＯＭ、ＤＲＡＭ、ＶＲＡＭ、フラッシュメモリ、フラッシュカード、磁気カード、光カード、ナノシステム、分子メモリ集積回路、ＲＡＩＤ、リモートデータストレージ／アーカイブ／ウェアハウジング、及び／又は命令及び／又はデータの格納に好適な他の任意の種類のデバイスが含まれ得る。
【0084】
コンピュータ可読媒体、命令記憶部又は記憶装置のいずれかに格納されたいくつかの実装には、システムのハードウェアを制御し、システム又はマイクロプロセッサに人間のユーザ又は本明細書に記載の例示的実施形態の結果を利用する他の機構との対話可能とする、ソフトウェアが含まれる。そのようなソフトウェアには、これに限らないが、デバイスドライバ、オペレーティングシステム、及びユーザアプリケーションが含まれ得る。究極的には、そのようなコンピュータ可読媒体又は記憶装置には、上で説明したような本発明の例示的態様を実行するためのソフトウェアが更に含まれる。
【0085】
プログラミング及び／又はシステムのソフトウェアには、本明細書に記載の手順を実行するためのソフトウェアモジュールが含まれる。本明細書のいくつかの例示的実施形態において、モジュールにはソフトウェアが含まれる。ただし本明細書の他の例示的実施形態ではモジュールにはハードウェア又はハードウェアとソフトウェアの組み合わせが含まれる。
【0086】
以上本発明の様々な例示的実施形態を述べたが、これらは例示のために提示したものであり、限定するためのものではないことを理解されたい。形式及び詳細において様々な変更をなし得ることは当業者には明らかであろう。したがって、本発明は上記のいかなる例示的実施形態によっても制限されるものではなく、以下に述べる特許請求の範囲及びその均等物によってのみ規定されるべきである。
【0087】
さらに、要約は、特許庁及び一般人、特に特許又は法律用語又は文体に精通していない科学者、技術者及び実務者が、大まかな調査で本出願の技術的開示の本質及び要点を迅速に判定できるようにすることを目的とするものである。要約は、いずれにしても本明細書に提示する例示的実施形態の範囲に関して制限的であることを意図するものではない。また、特許請求の範囲に説明される任意の手順は、必ずしも提示された順番で遂行されることを必要としないことも理解されたい。
【0088】
本明細書には多くの特定の実施形態の詳細が含まれているが、それらは任意の発明又は特許請求されるものの範囲を制限するものとみなすべきではなく、むしろ本明細書に記載の具体的な実施形態に固有の特徴の記述として理解されるべきである。本明細書において別々の実施形態の文脈で記載されている特定の特徴は、組み合わせて単一の実施形態に実装することもできる。その逆に、単一の実施形態の文脈で記載されている様々な特徴は、複数の実施形態に別々に、又は任意の適切な部分的な組み合わせにして実装することもできる。さらに、上記において特徴が特定の組み合わせで作用するものとして記述され、また最初にそのように特許請求されているとしても、特許請求された組合せからの１つ以上の特徴は場合によってはその組み合わせから除外することが可能であり、特許請求された組み合わせは部分的な組み合わせ又は部分的な組み合わせの変形に向けることも可能である。
【0089】
特定の状況下では、マルチタスク及び並列処理をすることも有利であり得る。さらに、上記の実施形態における様々な構成要素の分離は、すべての実施形態においてそのような分離を必要とすると理解するべきではなく、また、上記のプログラム要素及びシステムは一般に単体のソフトウェア製品に統合し、あるいは複数のソフトウェア製品にパッケージ化され得ることを理解されたい。
【0090】
ここまでにいくつかの例示的実施形態及び実施形態を説明したが、前述したものは例示であって限定ではなく、一例として提示されたことは明らかである。特に、本明細書に示した実施例の多くは、装置又はソフトウェア要素の特定の組み合わせを含むが、これらの要素は同じ目的を達成するために違う形で組み合わされてもよい。一つの実施形態に関連してのみ議論した作用、要素及び特徴は、他の実施形態における同様の役割から除外されることを意図するものではない。
【図1A】
【図1B】
【図1C】
【図2】
【図3A】
【図3B】
【図4】
【図5】
【図6】
【図7】
【図8A】
【図8B】
【図9】
【図10】
