特許 7278445 光コヒーレンストモグラフィ画像を用いた３次元解析

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

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

(11)【特許番号】

(24)【登録日】2023-05-11

(45)【発行日】2023-05-19

(54)【発明の名称】光コヒーレンストモグラフィ画像を用いた３次元解析

(51)【国際特許分類】

   A61B 3/10 20060101AFI20230512BHJP

   G06T 7/10 20170101ALI20230512BHJP

   G06T 7/136 20170101ALI20230512BHJP

【ＦＩ】

A61B3/10 100

G06T7/10

G06T7/136

A61B3/10 ZDM

【請求項の数】 5

(21)【出願番号】P 2022030645

(22)【出願日】2022-03-01

(62)【分割の表示】P 2020134831の分割

【原出願日】2020-08-07

(65)【公開番号】P2022082541

(43)【公開日】2022-06-02

【審査請求日】2022-03-03

(31)【優先権主張番号】16/845,307

(32)【優先日】2020-04-10

(33)【優先権主張国・地域又は機関】US

(73)【特許権者】

【識別番号】000220343

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

(74)【代理人】

【識別番号】100124626

【弁理士】

【氏名又は名称】榎並 智和

(72)【発明者】

【氏名】メイ・ソン

(72)【発明者】

【氏名】マオ・ザイシン

(72)【発明者】

【氏名】スイ・シン

(72)【発明者】

【氏名】ワング・ゼングォ

(72)【発明者】

【氏名】チャン・キンプイ

【審査官】▲高▼木 尚哉

(56)【参考文献】

【文献】特表２０１３－５４２８４０（ＪＰ，Ａ）

【文献】特開平１１－２７２８６５（ＪＰ，Ａ）

【文献】特表２００９－５０２３５４（ＪＰ，Ａ）

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

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

Ｇ０６Ｔ ７／００－７／９０

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

【請求項1】

被検者の対象の３次元光コヒーレンストモグラフィボリューメトリックデータを１回のスキャンによって取得するステップと、
前記３次元光コヒーレンストモグラフィボリューメトリックデータに複数の前処理法を用いた前処理を施すことによって前処理データを生成するステップと、
前記前処理データに複数のセグメンテーション法を用いたセグメンテーションを適用することで前記対象の生理学的要素を前記前処理データからセグメント化し、前記複数のセグメンテーション法で得られた複数のデータを合成することによって３次元セグメントデータを生成するステップと、
前記３次元セグメントデータを解析することによって前記３次元光コヒーレンストモグラフィボリューメトリックデータの２次元メトリックを求めるステップと、
前記２次元メトリックの可視化を行うステップと
を含み、
前記複数のセグメンテーション法は、
前記前処理データから得られた第１の２次元画像であるＢスキャン画像に対する第１の局所閾値処理セグメンテーション法と、
前記前処理データから得られた第２の２次元画像であるＣスキャン画像に対する第２の局所閾値処理セグメンテーション法と
を含み、
前記第１の局所閾値処理セグメンテーション法で前記Ｂスキャン画像から得られたデータと前記第２の局所閾値処理セグメンテーション法で前記Ｃスキャン画像から得られたデータとを含む複数のデータを合成することによって前記３次元セグメントデータを生成する、
３次元定量化方法。

【請求項2】

前記第１の局所閾値処理セグメンテーション法は、前記Ｂスキャン画像から前記生理学的要素のセグメントデータを生成し、
前記第２の局所閾値処理セグメンテーション法は、前記Ｃスキャン画像から前記生理学的要素のセグメントデータを生成する、
請求項１の方法。

【請求項3】

前記複数のセグメンテーション法は、前記前処理データの全体に対する大域閾値処理セグメンテーション法を更に含む、
請求項１又は２の方法。

【請求項4】

前記大域閾値処理セグメンテーション法は、前記前処理データから前記生理学的要素のセグメントデータを生成する、
請求項３の方法。

【請求項5】

前記対象は、網膜であり、
前記生理学的要素は、脈絡膜血管系である、
請求項１～４のいずれかの方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、光コヒーレンストモグラフィ画像を用いた３次元解析に関する。

【背景技術】

【0002】

光コヒーレンストモグラフィー（ＯＣＴ）は、様々な生体組織の（例えば、２次元スライス及び／又は３次元ボリュームとしての）インビボイメージング及び解析のための技術である。３次元（３Ｄ）ボリュームＯＣＴデータから作成された画像は、画像化された組織の更なる要素の異なる見え方／明るさを示す。更なる解析及び／又は可視化のために、それらの要素をこの違いに基づき画像からセグメント化することができる。例えば、脈絡膜血管系は、ＯＣＴ画像において、脈絡膜間質よりも暗く見える。したがって、ＯＣＴ画像における脈絡膜血管系は、強度閾値を適用することによってセグメント化することができる。しかし、ＯＣＴイメージングの固有の性質のため、閾値を画像に直接に適用すると、血管のセグメンテーションにアーチファクトが発生する。このようにして、ＯＣＴデータの要素をセグメント化するための他の技術が開発されてきたが、これらにおいても様々な欠陥や限界が問題となっている。

【0003】

例えば、局所二値化法によって脈絡膜の内腔領域及び間質領域を決定する場合、解析される深さにおいて十分な品質を達成し、且つ、利用されるＯＣＴシステムのタイプに応じて生じるノイズを回避するためには、特殊な画像収集プロトコルと平均化されたラインスキャンとが要求される。更に、最終的な閾値は手作業で適用される。２次元（２Ｄ）投影画像における脈絡膜血管密度測定を用いると、深さ分解能が不足し、シャドウアーチファクトが問題となる可能性がある。同様に、血管境界の自動検出（機械学習を用いる場合においても）は、シャドウアーチファクトの影響を受ける可能性があることに加え、２次元（２Ｄ）Ｂスキャンのみに適用が限定され、また、大きな血管への適用に限定される。更に、ボリューム全体にセグメンテーションが適用されるのではなく、ボリューム内の各Ｂスキャンに対してセグメンテーションが繰り返されるため、セグメント化された血管の連続性が悪化する可能性がある。このため、セグメント化されたボリュームを生成するためには、セグメント化されたＢスキャンのそれぞれをスプライシングするか、又は他の方法でつなぎ合わせなければならないことがある。他のセグメンテーション技術は、正常眼（非疾患眼）にのみ適用可能であり、また、疾患により網膜構造が変化した場合にはエラーを生じる。更に、いくつかのセグメンテーションは、基礎となるデータへのノイズ低減フィルタの適用に関連した不正確性を生じる。

【0004】

要するに、ノイズ除去（ノイズ低減）を行わずに、脈絡膜血管系が適切にセグメント化され得るデータを生成するためには、繰り返し実行されたＢスキャンの平均化又は深さ方向に沿った平均化が必要である。その結果、セグメンテーションは、次元及び位置によって制限され得る。更に、３Ｄデータへの適用においては、計算時間が非常に長くなり、解析可能なデータが限られてしまう。

【0005】

これらの制限のために、多くの臨床的に価値のある脈絡膜血管系の可視化及び定量化を提供することは、実用的ではなく、及び／又は可能でさえなかった。例えば、３Ｄボリューメトリックデータ又は結果として得られる画像に対して定量的解析が実行されても、その結果として得られる計量（メトリック）は３Ｄ情報を単一の値に圧縮するものである。これは、データの価値を大幅に低下させ、データを完全に利用しないものである。他の例では、定量化は、解析を正確に行うにはノイズが残り過ぎているＯＣＴデータから得られ、多くのボリュームから取得された平均を利用しており、それはやはりノイズの悪影響を受けており且つ（平均が取られる反復的なボリューム形態それぞれについて）スキャン時間の増加が必要とされ、又は比較的小さな関心領域に限定されている（例えば、単一のＢスキャンにおいて中心窩の下方１．５ミリメートル）。したがって、医療従事者は、臨床的に適切な情報を取得可能な３ＤボリューメトリックＯＣＴデータを完全に理解することができていない。

【発明の概要】

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

【0006】

本発明の目的の一つは、３次元ボリューメトリックＯＣＴデータの臨床的に価値のある解析及び可視化を提供することにある。

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

【0007】

実施形態の３次元定量化方法は、被検者の対象の３次元光コヒーレンストモグラフィボリューメトリックデータを１回のスキャンによって取得するステップと、前記３次元光コヒーレンストモグラフィボリューメトリックデータに複数の前処理法を用いた前処理を施すことによって前処理データを生成するステップと、前記前処理データに複数のセグメンテーション法を用いたセグメンテーションを適用することで前記対象の生理学的要素を前記前処理データからセグメント化することによって３次元セグメントデータを生成するステップと、前記３次元セグメントデータを解析することによって前記３次元光コヒーレンストモグラフィボリューメトリックデータの２次元メトリックを求めるステップと、前記２次元メトリックの可視化を行うステップとを含み、前記複数のセグメンテーション法は、前記前処理データから得られたＢスキャン画像に対する第１の局所閾値処理セグメンテーション法と、前記前処理データから得られたＣスキャン画像に対する第２の局所閾値処理セグメンテーション法とを含む。

【発明の効果】

【0008】

例示的な実施形態によれば、３次元ボリューメトリックＯＣＴデータの臨床的に価値のある解析及び可視化を提供することが可能である。

【図面の簡単な説明】

【0009】

【図1】本開示による例示的な方法のフローチャートを示す。

【図2】本開示による前処理及びセグメンテーションの適用例を示す。

【図3】本開示にしたがって生成された合成画像の例を示す。

【図4】本開示による可視化の例を示す。

【図5A】本開示による可視化の例としての脈絡血管の２次元ボリュームマップの例を示す。

【図5B】本開示による可視化の例としての脈絡血管の２次元ボリュームマップの例を示す。

【図6】本開示による可視化の例としての脈絡膜ボリュームのトレンドを示す。

【図7】本開示による可視化の例としての血管ボリュームを示す。

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

【0010】

本開示は、既知の技術では以前には実用的でなく及び／又は不可能であった３次元（３Ｄ）ボリューメトリックＯＣＴデータの臨床的に価値のある解析及び可視化に関するものである。そのような解析及び可視化は、疾患の診断、監視、及び治療管理を行う当業者の能力を改善し得る。簡単に言えば、一連の前処理に続いて、３次元における関心要素（例えば、脈絡膜血管系）についてのＯＣＴデータをセグメント化することによって解析が実行され、可視化が行われる。前処理に続いてデータにセグメンテーションを適用することができ、その後、所望の要素の最終的且つ完全な３次元セグメンテーションを生成するために合成することができる。その後、セグメント化された要素に対して、スムージング技術のような後処理を適用してもよい。本明細書ではＯＣＴデータの脈絡膜血管系について特に論じるが、本開示はこれに限定されるものではない。

【0011】

本開示による臨床的に価値のある解析及び可視化を行うための例示的な方法を図１に示す。同図に示すように、３ＤボリューメトリックＯＣＴデータが取得され、且つ、対応する生画像が、被検眼の撮像（１００）によって生成される。画像はグラフィック形成の基礎となるデータ表現であることに鑑み、以下、用語「画像」及び用語「データ」を互換的に使用する。撮像に続いて、個々の２次元画像（又は、３次元ボリュームとしてひとまとめにされる多数の２次元画像）に前処理が適用される（１０２）。この前処理（１０２）は、例えば、ディープラーニングに基づくノイズ低減法を適用することによってデータ及び画像のスペックル及び他のノイズに対処するものであってよい。このノイズ低減法は、例えば、２０２０年２月２１日に出願された「Ｉｍａｇｅ Ｑｕａｌｉｔｙ Ｉｍｐｒｏｖｅｍｅｎｔ Ｍｅｔｈｏｄｓ ｆｏｒ Ｏｐｔｉｃａｌ Ｃｏｈｅｒｅｎｃｅ Ｔｏｍｏｇｒａｐｈｙ（光コヒーレンストモグラフィのための画像品質改善方法）」と題する米国特許出願１６／７９７８４８号に記載されたものであってよい。米国特許出願１６／７９７８４８号の内容は、本開示の一部を構成するものとして、その全体が参照により援用される。更に、シャドウ（影）アーチファクト及びプロジェクション（投影）アーチファクトについては、画像処理及び／又はディープラーニング技術を適用することによって低減することができる。このアーチファクト低減法は、例えば、２０１９年９月２８日に出願された「３Ｄ Ｓｈａｄｏｗ Ｒｅｄｕｃｔｉｏｎ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｍｅｔｈｏｄ ｆｏｒ Ｏｐｔｉｃａｌ Ｃｏｈｅｒｅｎｃｅ Ｔｏｍｏｇｒａｐｈｙ （ＯＣＴ） Ｉｍａｇｅｓ（ＯＣＴ画像のための３Ｄ影低減信号処理方法）」と題する米国特許出願１６／５７４４５３号に記載されたものであってよい。米国特許出願１６／５７４４５３号の内容は、本開示の一部を構成するものとして、その全体が参照により援用される。もちろん、他のノイズ除去技術を適用してもよい。

【0012】

深さの次元（深さ方向）に沿った強度減衰については、強度補償技術及びコントラスト強調技術を適用することによって対処することができる。そのような技術は、例えば、所望の深さ及び関心領域（２Ｄ又は３Ｄ）における局所ラプラシアンフィルターとして、局所的に適用されてもよい。これに加えて、又はこれに代えて、コントラスト制限適応ヒストグラム均等化法（ｃｏｎｔｒａｓｔ－ｌｉｍｉｔｅｄ ａｄａｐｔｉｖｅ ｈｉｓｔｏｇｒａｍ ｅｑｕａｌｉｚａｔｉｏｎ、ＣＬＡＨＥ）を、コントラストを強調するために用いてもよい。もちろん、他のコントラスト強調技術（局所的又は大域的に適用される）及び／又は他の前処理技術を適用してもよい。

【0013】

前処理（１０２）は、画像又はボリュームの全体に適用されてもよいし、選択された関心領域のみに適用されてもよい。その結果、前処理（１０２）に入力された生画像又はボリュームのそれぞれについて、複数の前処理画像を生成することができる。換言すると、生ボリューメトリックＯＣＴデータから得られたＢスキャン又はＣスキャンのそれぞれに対して異なる前処理法を適用することで、複数の前処理画像を生成することができる。前処理（１０２）に続いて、前処理画像（又は画像の基礎となるデータ）は、画像／データの所望の要素（例えば、脈絡膜血管系）についてセグメント化される（１０４）。セグメンテーション（１０４）は、１つ以上の異なる技術を利用してもよい。この場合、適用される各セグメンテーション法は、個別に、比較的単純且つ高速に実行されてよく、異なる長所及び短所を有してよい。

【0014】

例えば、いくつかのセグメンテーション法は、異なる閾値レベルを利用してもよく、及び／又は、異なるビュー（例えば、Ｂスキャン又はＣスキャン）の解析に基づくものであってもよい。より詳細には、Ｃスキャンに対してセグメンテーションを実行する場合、各Ｃスキャン画像がボリュームの視野全体の情報を含むので、Ｂスキャンにセグメンテーションを施す場合と比べて血管の連続性を向上させることが可能である。これは更に、Ｂスキャンでのセグメンテーションに比べてより小さい血管のセグメンテーションを可能にするとともに、ユーザが行うセグメンテーションの手動検証を容易にする。しかしながら、Ｃスキャンのセグメンテーションは、ボリューメトリックデータを平坦化するために使用される先行のブルッフ膜のセグメンテーションの正確度に依存することがある。

【0015】

上記を考慮すると、複数の前処理画像のうちの１つ以上に対して、異なるセグメンテーション技術を選択的に適用することができる。更に、上述したように、ノイズ及び減衰（例えば、アーチファクトを引き起こすもの）の影響により、ＯＣＴボリューム全体に対して大域的なセグメンテーションを施すことは、実用上、不可能であった。しかしながら、上述した前処理の適用に続いて、ボリュームから得られた個々のＢスキャン又はＣスキャンではなく、ＯＣＴボリューム全体に対してセグメンテーションを適用してもよい。いずれにしても、各々のセグメンテーションは、前処理画像／データにおける所望の要素をセグメント化するものである。ボリューム全体に適用されるセグメンテーションは、複数のセグメンテーションをつなぎ合わせる必要がないため、セグメンテーションの接続性を更に向上させることができるが、そのような複数のセグメンテーションは、比較的低いコントラストを有するボリュームの局所領域に対する感度が低いかもしれない。ただし、上述した深度補償技術及びコントラスト強調技術によってこれを軽減し得る。

【0016】

一つの例示的な実施形態では、個別に前処理が施された画像／データに対して各セグメンテーション法を適用してもよい。他の実施形態では、異なる関心領域に対応する画像／データに対して複数のセグメンテーション法を選択的に適用してもよい。例えば、第１のセグメンテーション法にしたがって第１の２つの前処理画像をセグメント化しつつ、第２のセグメンテーション法にしたがって第２の２つの前処理画像をセグメント化することができる。他の実施形態では、３ＤボリューメトリックＯＣＴデータの前処理を任意個数の手法にしたがって実行した後、脈絡膜血管系の第１の決定を得るために、前処理された３ＤボリューメトリックＯＣＴデータから得られたＢスキャン画像に対して局所閾値処理セグメンテーション法が適用され、脈絡膜血管系の第２の決定を得るために、前処理された３ＤボリューメトリックＯＣＴデータから得られたＣスキャン画像に局所的な閾値処理法が適用され、更に、脈絡膜血管系の第３の決定を得るために、前処理された３Ｄボリューメトリックデータの全体に大域的な閾値処理法が適用される。

【0017】

適用される前処理法及びセグメンテーション法の個数にかかわらず、これらのセグメンテーションは、合成セグメント画像又はデータを生成するために組み合わされる。この合成セグメント画像又はデータは、解析（１０８）の一部として異なる定量的メトリックを求めるための処理と、セグメンテーション及び／又はメトリックの可視化（１１０）との両方について十分な品質を有し、且つ、アーチファクトを含まないものである。このようにして、合成画像は、前処理法及びセグメンテーション法の全てを含んでもよく、結合（ｕｎｉｏｎ）、交差（ｉｎｔｅｒｓｅｃｔｉｏｎ）、重み付け（ｗｅｉｇｈｔｉｎｇ）、投票（ｖｏｔｉｎｇ）などの任意の方法にしたがって結合されてもよい。セグメンテーション（１０４）に続いて、例えば平滑化のための後処理をセグメント画像又はデータに対して更に施してもよい。

【0018】

上記した前処理とセグメンテーションとの組み合わせが図２に模式的に示されている。図２に示す例は、それぞれが共通の３ＤボリューメトリックＯＣＴデータセットに由来する、生画像／生データの２つのサブセットを用いる。画像／データのこれらのサブセットは、関心領域、ビュー（例えば、Ｂスキャン及びＣスキャン）などにしたがって分けられてよい。図２の例によれば、第１のサブセット２００は第１の前処理２０２を受け、第２のサブセット２０４は第２の前処理２０６を受ける。図２に破線で示した他の実施形態では、各サブセット２００、２０４は、利用可能な前処理２０２、２０６のいずれかに供されてよい。このようにして、第１のサブセット２００に関連付けられたデータは少なくとも１つの前処理データサブセットをもたらし、第２のサブセット２０４に関連付けられたデータは少なくとも２つの前処理データサブセットをもたらす。前処理に続いて、結果として得られるデータセットのそれぞれが、任意の利用可能なセグメンテーション法（図２には３つのセグメンテーション法が例示されている）によって同様にセグメント化される。図示されているように、それぞれの前処理の結果は、異なるセグメンテーション法２０８、２１０、２１２によって別々にセグメント化されるが、破線で示した他の実施形態では、セグメンテーション法２０８、２１０、２１２のうちの１つ以上が、複数の前処理画像／複数の前処理データのうちのいずれかに適用されてよい。最後に、各セグメンテーション法２０８、２１０、２１２の出力が、上述したように合成されて（２１４）、合成セグメンテーションが生成される。上記を考慮すると、生画像／生データが由来する３ＤボリューメトリックＯＣＴデータの単一の合成セグメンテーションを生成するための方法の一部として、異なる前処理及び／又はセグメンテーション技術を共通の生画像／生データに適用してもよい。

【0019】

上述したように、単一の複雑な前処理及びセグメンテーションを実行するのではなく、複数の前処理法及び複数のセグメンテーション法を用いて合成結果を生成することによって、前処理及びセグメンテーションの総時間及び計算量が低減される。それにもかかわらず、同等の品質を達成することが可能であり、且つ、３Ｄボリュームの全体にセグメンテーションを適用することが可能である。したがって、結果として得られるセグメンテーションは、ノイズ及びシャドウアーチファクトを無くすことができ、可視化及び定量化のために十分な品質を提供することができる（後述）。上記の例にしたがう例示的な合成画像を図３に示す。図３の合成画像には、３ＤボリューメトリックＯＣＴデータからセグメント化された脈絡膜血管が３Ｄビューで描出されている。

【0020】

図１を再び参照する。合成画像又はボリュームは、次に、ボリュームの定量解析のために以前に使用されたＢスキャンの２次元データではなく、ボリューメトリックＯＣＴデータ全体に基づいて、多数の定量可能なメトリック１０８を生成して解析するために処理されてもよい。これらのメトリックは、上述の前処理が施されセグメント化されたＯＣＴデータから生成されるので、メトリックは、従来の技術にしたがってＯＣＴデータから得られたものよりも有意に正確性が高くなるであろう。更に、メトリック（及び図３に示すようなセグメント化された可視化、及びメトリックから生成された任意の可視化）は、ボリューム又はボリューム全体の複数の２Ｄ画像にわたる比較的大きな領域（例えば、単一のＢスキャンの１．５ｍｍよりも大きい領域）に関して決定されてよく、且つ、（複数のスキャンの平均ではなく、単一のスキャンから取得された）単一のＯＣＴボリュームから決定されてもよい。

【0021】

例えば、３次元ボリューム内において、合成セグメント画像内のセグメント化されたデータをセグメント化されていないデータと比較することによって、血管系の空間ボリューム（及び、それに関連して、密度は、血管系又は同様のセグメント化された要素に相当する所定の領域内のボリューム全体の割合である）、直径、長さ、体積比（インデックスとも呼ばれる）などを特定することができる。例えば、セグメント化されたピクセルの個数を数えることで、関心領域内の血管系の量（例えば、体積又は密度）の指標を提供することができる。それらのメトリックのみを深さなどの１次元に投影（例えば、最大値、最小値、平均値、和などを演算）することにより、体積マップ、直径マップ、インデックスマップなどを生成することができる。このようなマップにより、網膜上の各位置についてのメトリックの定量化された値を視覚的に示すことが可能である。更に、それらのメトリックを（例えば、マップ全体にわたって）一つの次元方向に総計することにより、全体の体積、代表的なインデックスなどを特定することが可能である。単一のＯＣＴボリュームからそのようなメトリックを広い範囲にわたって定量化することにより、被検者間におけるボリューメトリックＯＣＴデータの比較や、個々の被検者の経時的な比較といった、従来は不可能であった比較が可能になる。

【0022】

これらのメトリックは、比較によるもの（相対的なもの）であってもよい。例えば、比較メトリックの基になるＯＣＴボリュームのメトリックは、異なる時間に単一の被検者から取得されたもの、異なる眼（例えば、単一の被検者の右眼及び左眼）から取得されたもの、（例えば、個体群を代表する集合的個体のそれぞれの間の）複数の被検者から取得されたもの、又は、同じ眼の異なる関心領域（例えば、異なる層）から得られたものであってよい。これらの比較は、比較の各要素についてのメトリックを決定して任意の統計的比較手法を適用することによって行うことができる。例えば、比較メトリックは、比較データの比、比較データ間の差、比較データの平均、比較データの合計などのいずれかであってよい。比較は、一般的にボリューメトリックデータ全体に対して行われてもよいし、場所ごとに行われてもよい（例えば、比較マップのピクセル位置ごとに行われてよい）。

【0023】

共通の関心領域から得られたメトリックを比較する場合、比較される要素（異なるデータセット、画像、ボリューム、メトリックなど）は、好ましくは、同様の比較を行うことができるように、互いに位置合わせ（レジストレーション）される。言い換えれば、レジストレーションを行うことによって、各要素の対応する部分を比較することが可能になる。いくつかの例では、例えば、脈絡膜血管系の変化を比較する場合、血管系が各要素において必ずしも同一ではないため、血管系自体に基づくレジストレーションが行われなくてよい（例えば、比較される期間にわたる治療に起因する）。より一般的に言えば、レジストレーションは、好ましくは、要素間で異なる可能性がある情報、又は比較されるメトリックにおいて使用される情報に基づいて実行されない。これを考慮して、いくつかの実施形態では、レジストレーションは、比較される各要素の生のＯＣＴボリューム（例えば、前処理が施されていないＯＣＴボリューム）から生成されたアンファス画像（ｅｎ ｆａｃｅ ｉｍａｇｅ、正面画像）に基づいて実行されてもよい。これらのアンファス画像は、レジストレーションのために使用される領域内の各Ａラインに沿った強度を加算、平均化などすることによって生成されてもよい。網膜血管は影（シャドウ）を投影することがあり、したがってＯＣＴアンファス画像上に影を生じることがあるので、アンファス画像はレジストレーションに有用であり、比較的安定した状態のままでいる暗い網膜血管系をランドマークとして用いることができる。更に、本来、任意のメトリック、脈絡膜血管系画像、又はＯＣＴボリュームから生成される同様の画像は、アンファス画像とともにレジストレーションされる。これは、それらが同じボリュームに由来するからである。例えば、第１のボリューム内の表面血管と第２のボリューム内の表面血管とをレジストレーションすることができ、第１のボリューム内の脈絡膜血管（又は、脈絡膜血管のメトリック）と第２のボリューム内の脈絡膜血管とを比較することができる。

【0024】

その後、これらのメトリックの可視化が実行されて表示に表示されてもよく（１１０）、後で参照するために保存されてもよい。すなわち、本開示に係る技術は、ボリューメトリックＯＣＴデータのセグメント化された要素（例えば、脈絡膜血管）の可視化だけでなく、そのセグメント化された要素に関連する定量化されたメトリックの可視化（例えば、マップ、グラフ）を生成することも可能である。これらの定量化されたメトリックの可視化は、上述の比較を更に単純化するものである。そのような可視化は、３Ｄボリューメトリック情報を表すメトリックの２Ｄ表現、並びに／又は、２以上のＯＣＴボリューム間の変化及び／若しくは相違（差分）を表す比較メトリックの表現であってよい。上述のメトリックを考慮すると、可視化は、例えば、脈絡膜血管インデックスマップ、脈絡膜厚さマップ、若しくは血管体積マップ、及び／又は、これらのいずれかの比較であってよい。

【0025】

可視化において情報が様々な形態で符号化されてもよい。例えば、可視化における各ピクセルの強度は、そのピクセルに対応する位置におけるメトリックの値を示してよく、また、色は、その値のトレンドを示してよい（又は、トレンドの表現に強度を使用し、値の表現に色を使用してもよい）。更に他の実施形態では、異なるメトリック情報を識別するために異なるカラーチャネルを使用してもよい（例えば、メトリックごとに異なる色を使用し、そのメトリックのトレンド又は値を強度で示す）。更に他の実施形態では、様々な形態の色相（色調）、彩度、及び値（ＨＳＶ）でのエンコーディング、及び／又は、色相（色調）、彩度、及び光（ＨＳＬ）でのエンコーディングを利用してもよい。更に他の実施形態では、付加的な情報をエンコードするために透明度を用いてもよい。例示的な可視化を図４～図７に示す。

【0026】

図４は、本開示に係る第１の例示的な可視化を示す。図４の可視化は、脈絡膜血管の２Ｄ画像であり、各ピクセルの強度は、メトリックに対応し、且つ、色は、以前のスキャンとの比較としてのメトリックの局所的なトレンドを示す。例えば、各ピクセルの強度は、３次元ボリューメトリックデータの血管体積、血管長、血管厚、又は同様の測定値に対応していてよい。色は、以前に取得された３Ｄボリューメトリックデータに基づく以前のメトリック測定との比較としての各ピクセルの変化を示していてもよい。例えば、赤色は、以前の測定から得られた血管系測定結果の拡大を示すために使用されてよく、紫色は、血管系測定結果の縮小を示すために使用されてよい。青及び緑は、比較的一貫した測定結果（すなわち、変化がほとんどないか、又は全くない）を示してよい。図４の白黒画像には色が明確に示されていないので、縮小（例えば、紫として識別される）及び拡大（例えば、赤として識別される）に対応する例示的な領域が、参考のために明示的に示されている。以前の測定値との比較は、単純な差、複数の測定値の平均に対する相対的な変化、標準偏差、及び／又は、同様の統計演算として求められてよい。もちろん、色と変化との間の相関関係は別途に設定されてもよい。

【0027】

図５Ａ及び５Ｂは、それぞれ、本開示に係る例示的な可視化としての例示的な脈絡膜血管２Ｄボリュームマップを示す。３Ｄボリューメトリックデータセットの脈絡膜血管のボリュームは、３Ｄボリュームデータの各Ａラインについての脈絡膜血管に対応するピクセルの個数に各ピクセルの解像度を乗算したものとして決定されてよい。この積算（累算）が深さにわたって行われる場合、ボリュームマップの各ピクセルは、３次元ボリュームデータセットの１つのＡラインに対応する。図４の例と同様に、ボリュームマップの各ピクセルの強度は、対応する位置での血管体積に相当し、色は、以前のスキャンとの比較としてのこの血管体積の局所的なトレンドに相当する。同様に、セグメント化されたピクセルの個数を脈絡膜（又は他の領域）のピクセルの総数と比較することによって、その領域にわたる血管密度（又は、他の要素の密度）の定量化を提供することができる。一般に、体積及び密度は、ともに増加及び減少するであろう。

【0028】

上述したように、２Ｄ可視化マップを生成するために使用されるメトリックは、更なる解析のために、関心領域にわたって更に累算されてもよい。例えば、中心窩（半径１ｍｍ程度）、傍中心窩（上方、鼻側、下方、耳側）（中心窩からの半径１－３ｍｍ程度）、周中心窩（上方、鼻側、下方、耳側）（中心窩からの半径３－５ｍｍ程度）、及び／又は、同様の領域について、メトリック値及び／又はピクセル強度を累算してもよい。累算は、加算、標準偏差などの任意の統計演算によって求めることができる。累算された数値が異なる複数の時点について収集された場合、トレンド解析を行い、対応するトレンドの可視化を生成することができる。累算された数値を患者間で比較することや、累算された数値を標準値（ｎｏｒｍａｔｉｖｅ ｖａｌｕｅ（ｓ））と比較することも可能である。

【0029】

中心窩及び周中心窩鼻側のそれぞれにおける脈絡膜体積トレンドの例示的な可視化を図６に示す。同図から分かるように、４週間の期間にわたって毎週、中心窩領域及び周中心窩鼻側領域のそれぞれにおいて、脈絡膜体積が集められた。この可視化から明らかなように、この被検者について、第１週目と第２週目の間に周中心窩鼻側において脈絡膜体積が増加を呈し、同じく第１週目と第２週目の間において中心窩における体積が対応する減少を呈していることが分かる。しかしながら、第３週目から中心窩での血管体積が増加し始めると、周中心窩鼻側における体積は元の値より減少している。第３週目から第４週目にかけては、双方の部位の体積が増加している。

【0030】

別の例の可視化を図７に示す。この可視化には、脈絡膜血管の総体積が、脈絡膜の異なるセクターについて示されている：中心窩（中心）、鼻側－上方（ＮＳ）、鼻側（Ｎ）、鼻側－下方（ＮＩ）、耳側－下方（ＴＩ）、耳側（Ｔ）、及び、耳側－上方（ＴＳ）。これらの総体積は、各セクター内の脈絡膜血管系ピクセルの総数を合計することによって決定されてよい。更に、３Ｄデータの解像度に基づいて、ピクセルの総数を物理的な寸法（立方ミリメートルなど）に変換してもよい。図７の可視化によれば、患者の治療前、治療の１ヶ月後、及び治療の１年後の体積が示されている。これから分かるように、治療後の各セクターにおいて血管系の体積が大きく減少している。

【0031】

もちろん、異なるメトリクスについて同様の２Ｄマップ及びトレンド可視化を生成してもよい。例えば、３Ｄボリューメトリックデータセットの各Ａラインについて脈絡膜血管系ピクセルの総数を求めることにより、血管厚マップ及びトレンドの可視化を生成してもよい。また、領域（脈絡膜など）内の非血管ピクセルの総数を求めることにより、非血管インデックスマップ及びトレンドの可視化を生成してもよい。

【0032】

上述した態様は、プロセッサによってハードウェア及び／又はソフトウェアを介して実装されることが想定される。プロセッサは、例えば、抵抗器、トランジスタ、コンデンサ、インダクタなどを含む、任意の数の電気部品で構成される任意の電気回路、又はその一部であってもよい。回路は、例えば、集積回路、集積回路の集合体、マイクロコントローラ、マイクロプロセッサ、プリント回路基板（ＰＣＢ）上のディスクリート電子部品の集合体などを含む任意の形態であってもよい。プロセッサは、ランダムアクセスメモリ、フラッシュメモリ、デジタルハードディスクなどのような揮発性又は不揮発性の形態のメモリに格納されたソフトウェア命令を実行可能であってよい。このプロセッサは、ＯＣＴ等の撮像システムにおけるプロセッサと一体化されていてもよいが、単独で設けられもよいし、画像データの処理以外の動作に使用されるコンピュータの一部であってもよい。

【図1】

【図2】

【図3】

【図4】

【図5A】

【図5B】

【図6】

【図7】