知財求人 - 知財ポータルサイト「IP Force」
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2022-12-16
(45)【発行日】2022-12-26
(54)【発明の名称】画像処理装置、画像処理方法及びプログラム
(51)【国際特許分類】
A61B 3/10 20060101AFI20221219BHJP
【ＦＩ】
A61B3/10 100
A61B3/10 300
A61B3/10 ZDM
【請求項の数】
10
(21)【出願番号】P 2018044563
(22)【出願日】2018-03-12
(65)【公開番号】P2019154718
(43)【公開日】2019-09-19
【審査請求日】2021-03-10
(73)【特許権者】
【識別番号】000001007
【氏名又は名称】キヤノン株式会社
(74)【代理人】
【識別番号】100126240
【弁理士】
【氏名又は名称】阿部 琢磨
(74)【代理人】
【識別番号】100124442
【弁理士】
【氏名又は名称】黒岩 創吾
(72)【発明者】
【氏名】今村 裕之
(72)【発明者】
【氏名】内田 弘樹
(72)【発明者】
【氏名】▲高▼橋 理宇眞
【審査官】▲高▼原 悠佑
(56)【参考文献】
【文献】特開２０１７－０７７４１４（ＪＰ，Ａ）
【文献】特開２０１７－１７０１９３（ＪＰ，Ａ）
【文献】特開２０１７－２２１５２５（ＪＰ，Ａ）
【文献】特開２０１７－１９６３０６（ＪＰ，Ａ）
【文献】特開２０１７－０４６９７６（ＪＰ，Ａ）
【文献】特開２０１０－２７９５３６（ＪＰ，Ａ）
【文献】特開２０１４－１４０４９０（ＪＰ，Ａ）
【文献】国際公開第２０１７／１５０５８３（ＷＯ，Ａ１）
【文献】国際公開第２０１５／０６４５４５（ＷＯ，Ａ１）
【文献】特開２０１９－１５４７１６（ＪＰ，Ａ）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ａ６１Ｂ ３／００－３／１８
(57)【特許請求の範囲】
【請求項1】
眼部のモーションコントラスト画像の異なる深度範囲の投影画像であって、複数の投影画像を取得する取得手段と、前記投影画像の第１の領域と、前記第１の領域の部分領域である第２の領域に対する解析を実行する解析手段と、
前記第１の領域に対する解析結果を示す画像に、前記第２の領域に対する解析結果を示す画像を重ねた合成画像を表示手段に表示させる表示制御手段と、を備え、
前記表示制御手段は、第１の深度範囲に対応する合成画像と、前記第１の深度範囲とは異なる第２の深度範囲に対応する合成画像とを並べて表示手段に表示させることを特徴とする画像処理装置。
【請求項2】
前記眼部のモーションコントラスト画像は、前記眼部の同一位置を測定光が走査されるように制御して得た複数の３次元のモーションコントラスト画像を合成した画像であることを特徴とする請求項１に記載の画像処理装置。
【請求項3】
前記表示制御手段は、前記解析手段により解析された結果を示す画像が、前記解析に適した複数の条件のうちの少なくとも２つの条件が満足されない状態で得られた画像である場合には、前記少なくとも２つの条件に関する情報を、前記複数の条件の優先順位に応じて前記表示手段に表示させることを特徴とする請求項２に記載の画像処理装置。
【請求項4】
前記表示制御手段は、前記少なくとも２つの条件に関する情報として、前記少なくとも２つの条件のうちの優先度の高い条件に関する警告メッセージを、操作者からの指示に応じて選択された解析の種類を示す情報を用いて解析された結果を示す画像と並べて前記表示手段に表示させることを特徴とする請求項３に記載の画像処理装置。
【請求項5】
前記複数の条件には、前記モーションコントラスト画像が前記合成した画像であるという条件が、他の条件よりも優先度の高い条件として含まれることを特徴とする請求項３に記載の画像処理装置。
【請求項6】
前記情報は、血管領域、無血管領域および血管中心線の位置のいずれかに基づいて算出された計測値に関する情報と、血管密度の情報とのうち少なくとも１つの情報であることを特徴とする請求項３乃至５のいずれか１項に記載の画像処理装置。
【請求項7】
前記モーションコントラスト画像は、検査日が異なる複数のモーションコントラスト画像を含み、前記表示制御手段は、前記投影画像それぞれから得られた複数の合成画像を、前記検査日に基づく時系列に並べて表示手段に表示させることを特徴とする請求項１乃至６のいずれか１項に記載の画像処理装置。
【請求項8】
操作者からの指示に応じて選択された解析の種類を示す情報を用いて、前記眼部のモーションコントラスト画像における第１の領域に対する解析を実行する解析手段を更に備え、前記表示制御手段は、操作者からの指示に応じて前記モーションコントラスト画像における前記第１の領域よりも狭い領域を少なくとも含む第２の領域に対して解析された結果を示す画像の表示が選択された場合には、前記選択された解析の種類を示す情報を用いて前記第２の領域に対して解析された結果を示す画像を、前記第１の領域に対して解析された結果を示す画像に重畳された状態で表示手段に表示させることを特徴とする請求項１乃至７のいずれか１項に記載の画像処理装置。
【請求項9】
眼部のモーションコントラスト画像の異なる深度範囲の投影画像であって、複数の投影画像を取得する取得工程と、前記投影画像の第１の領域と、前記第１の領域の部分領域である第２の領域に対する解析を実行する解析工程と、
前記第１の領域に対する解析結果を示す画像に、前記第２の領域に対する解析結果を示す画像を重ねた合成画像を表示手段に表示させる表示制御工程と、を含み、
前記表示制御工程において、第１の深度範囲に対応する合成画像と、前記第１の深度範囲とは異なる第２の深度範囲に対応する合成画像とを並べて表示させることを特徴とする画像処理方法。
【請求項10】
請求項９に記載の画像処理方法の各工程をコンピュータに実行させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本明細書の開示は、画像処理装置、画像処理方法及びプログラムに関する。
【背景技術】
【0002】
光干渉断層計（ＯＣＴ；Ｏｐｔｉｃａｌ Ｃｏｈｅｒｅｎｃｅ Ｔｏｍｏｇｒａｐｈｙ）を用いて非侵襲に眼底血管を描出するＯＣＴ Ａｎｇｉｏｇｒａｐｈｙ（以下、ＯＣＴＡと表記）が知られている。ＯＣＴＡでは測定光で同一位置を複数回走査し、複数のＯＣＴ断層画像を取得する。この複数のＯＣＴ断層画像に基づいて赤血球の変位と測定光との相互作用により得られるモーションコントラストデータがＯＣＴＡ画像として画像化される。
【0003】
特許文献１は、取得時期（検査日時）が異なる複数のモーションコントラストデータそれぞれに対して算出した血管解析マップを時系列で並べて表示する技術を開示している。
【先行技術文献】
【特許文献】
【0004】
【文献】特開２０１７－７７４１４号公報
【発明の概要】
【発明が解決しようとする課題】
【0005】
眼部血管を時系列で評価することにより眼底血管の変化を定量的に把握できることが期待される。
【0006】
しかし、検査ごとのＯＣＴ断層画像の信号強度または画質のばらつきに応じて、ＯＣＴＡ画像にも検査ごとに変化が生じてしまう。すなわち、例え眼部血管に経時的な変化が起きていない場合でも、ＯＣＴＡ画像における眼底血管には経時的な変化があらわれてしまう場合があった。すなわち、眼底血管の経時変化を適切に評価できない場合があった。
【0007】
本明細書の開示は、眼底血管に関する経時変化の適切な評価を支援することを目的の１つとする。
【0008】
なお、前記目的に限らず、後述する発明を実施するための形態に示す各構成により導かれる作用効果であって、従来の技術によっては得られない作用効果を奏することも本明細書の開示の他の目的の１つとして位置付けることができる。
【課題を解決するための手段】
【0009】
本発明の目的を達成するために、例えば本明細書に開示の画像処理装置は、
眼部のモーションコントラスト画像の異なる深度範囲の投影画像であって、複数の投影画像を取得する取得手段と、
前記投影画像の第１の領域と、前記第１の領域の部分領域である第２の領域に対する解析を実行する解析手段と、
前記第１の領域に対する解析結果を示す画像に、前記第２の領域に対する解析結果を示す画像を重ねた合成画像を表示手段に表示させる表示制御手段と、を備え、
前記表示制御手段は、第１の深度範囲に対応する合成画像と、前記第１の深度範囲とは異なる第２の深度範囲に対応する合成画像とを並べて表示手段に表示させる。
眼部のモーションコントラスト画像の異なる深度範囲の投影画像であって、複数の投影画像を取得する取得手段と、
前記投影画像の第１の領域と、前記第１の領域の部分領域である第２の領域に対する解析を実行する解析手段と、
前記第１の領域に対する解析結果を示す画像に、前記第２の領域に対する解析結果を示す画像を重ねた合成画像を表示手段に表示させる表示制御手段と、を備え、
前記表示制御手段は、第１の深度範囲に対応する合成画像と、前記第１の深度範囲とは異なる第２の深度範囲に対応する合成画像とを並べて表示手段に表示させる。
【発明の効果】
【0010】
本明細書の開示によれば、眼底血管に関する経時変化の適切な評価を支援することができる。
【図面の簡単な説明】
【0011】
【図1】第１実施形態に係る画像処理装置の構成の一例を示すブロック図である。
【図2】実施形態に係る画像処理システムの一例および該画像処理システムを構成する断層画像撮影装置に含まれる測定光学系の一例を説明する図である。
【図3】第１実施形態に係る画像処理システムが実行可能な処理の一例を示すフローチャートである。
【図4】実施形態におけるＯＣＴＡ撮影の走査方法の一例を説明する図である。
【図5】第１実施形態のＳ３０７で実行される処理の一例を説明する図である。
【図6】第１実施形態のＳ３０８で実行される処理の一例を説明する図である。
【図7】第１実施形態において表示手段に表示する基準検査の選択画面の一例及び撮影画面の一例を説明する図である。
【図8】第１実施形態のＳ３０４における画像処理内容の一例及びＳ３０５において表示手段に表示するレポート画面の一例を説明する図である。
【図9】第１実施形態において表示手段に表示する計測操作画面の一例とＳ３０８において表示する計測レポート画面の一例を説明する図である。
【図10】第１実施形態において特定された血管領域をユーザが修正する場合の操作手順の一例と実行される画像処理内容の一例を説明する図である。
【図11】第１実施形態のＳ３１１において表示手段に表示する経時変化計測レポート画面の一例を説明する図である。
【発明を実施するための形態】
【0012】
［第１の実施形態］
本実施形態に係る画像処理装置は、同一被検眼に対して異なる日時に略同一撮影条件で取得した重ね合わせＯＣＴＡ画像（複数のＯＣＴＡ画像の合成画像）から生成した網膜表層及び網膜深層の正面モーションコントラスト画像を用いて血管領域特定及び血管密度計測処理を実施する。該特定処理及び計測処理により得られた合成画像及び計測値を複数の深度範囲で時系列並置表示する場合について説明する。
本実施形態に係る画像処理装置は、同一被検眼に対して異なる日時に略同一撮影条件で取得した重ね合わせＯＣＴＡ画像（複数のＯＣＴＡ画像の合成画像）から生成した網膜表層及び網膜深層の正面モーションコントラスト画像を用いて血管領域特定及び血管密度計測処理を実施する。該特定処理及び計測処理により得られた合成画像及び計測値を複数の深度範囲で時系列並置表示する場合について説明する。
【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】
【数1】
【0060】
ここで、Ａｘｙは断層画像データＡの位置（ｘ，ｙ）における（ＦＦＴ処理後の複素数データの）振幅、Ｂｘｙは断層データＢの同一位置（ｘ，ｙ）における振幅を示している。０≦Ｍｘｙ≦１であり、両振幅値の差異が大きいほど１に近い値をとる。モーションコントラストデータ生成部１０１－１２は、式（１）のような脱相関演算処理を（同一クラスタに属する）任意の時間的に隣接する断層画像間で行う。モーションコントラストデータ生成部１０１－１２は、得られた（１クラスタあたりの断層画像数－１）個のモーションコントラスト値の平均を画素値として持つ画像を最終的なモーションコントラスト画像として生成する。
【0061】
なお、ここではＦＦＴ処理後の複素数データの振幅に基づいてモーションコントラストを計算したが、モーションコントラストの計算法は上記に限定されない。例えば、モーションコントラストデータ生成部１０１－１２は、複素数データの位相情報に基づいてモーションコントラストを計算してもよいし、振幅と位相の両方の情報に基づいてモーションコントラストを計算してもよい。あるいは、複素数データの実部や虚部に基づいてモーションコントラストを計算してもよい。
【0062】
また、本実施形態ではモーションコントラストとして脱相関値を計算したが、モーションコントラストの計算法はこれに限定されない。例えば二つの値の差分に基づいてモーションコントラストを計算しても良いし、二つの値の比に基づいてモーションコントラストを計算してもよい。
【0063】
さらに、上記では取得された複数の脱相関値の平均値を求めることで最終的なモーションコントラスト画像を得ているが、本発明はこれに限定されない。例えば取得された複数の脱相関値の中央値、あるいは最大値を画素値として持つ画像を最終的なモーションコントラスト画像として生成しても良い。
【0064】
＜ステップ３０４＞
画像処理部１０１－０４は、繰り返しＯＣＴＡ撮影を通して得られたモーションコントラスト画像群を３次元的に位置合わせし、加算平均する。複数のクラスタから得られた複数のモーションコントラスト画像の加算平均を行うことで図８（ｂ）に示すように高コントラストな合成モーションコントラスト画像を生成する。図８（ａ）は、図８（ｂ）との対比のため１つのクラスタから得られるモーションコントラスト画像を示している。なお、合成処理は単純加算平均に限定されない。例えば各モーションコントラスト画像の輝度値に対して任意の重みづけをした上で平均した値でもよいし、中央値をはじめとする任意の統計値を算出してもよい。また位置合わせ処理をＥｎＦａｃｅ画像の状態で行う場合、すなわち位置合わせ処理を２次元的に行う場合も本発明に含まれる。
画像処理部１０１－０４は、繰り返しＯＣＴＡ撮影を通して得られたモーションコントラスト画像群を３次元的に位置合わせし、加算平均する。複数のクラスタから得られた複数のモーションコントラスト画像の加算平均を行うことで図８（ｂ）に示すように高コントラストな合成モーションコントラスト画像を生成する。図８（ａ）は、図８（ｂ）との対比のため１つのクラスタから得られるモーションコントラスト画像を示している。なお、合成処理は単純加算平均に限定されない。例えば各モーションコントラスト画像の輝度値に対して任意の重みづけをした上で平均した値でもよいし、中央値をはじめとする任意の統計値を算出してもよい。また位置合わせ処理をＥｎＦａｃｅ画像の状態で行う場合、すなわち位置合わせ処理を２次元的に行う場合も本発明に含まれる。
【0065】
なお、合成部１０１－４２が合成処理に不適なモーションコントラスト画像が含まれているか否かを判定した上で、不適と判定したモーションコントラスト画像を除いて合成処理を行うよう構成してもよい。例えば、各モーションコントラスト画像に対して評価値（例えば脱相関値の平均値や、ｆＳＮＲ）が所定の範囲外である場合に、合成処理に不適と判定すればよい。
【0066】
本実施形態では合成部１０１－４２がモーションコントラスト画像を３次元的に合成した後、補正部１０１－４３がモーションコントラスト画像内に生じるプロジェクションアーチファクトを３次元的に抑制する処理を行う。
【0067】
ここで、プロジェクションアーチファクトは網膜表層血管内のモーションコントラストが深層側（網膜深層や網膜外層・脈絡膜）に映り込み、実際には血管の存在しない深層側の領域に高い脱相関値が生じる現象を指す。図８（ｃ）に、３次元ＯＣＴ断層画像上に３次元モーションコントラストデータを重畳表示した例を示す。網膜表層血管領域に対応する高い脱相関値を持つ領域８０１の深層側（視細胞層）に、高い脱相関値を持つ領域８０２が生じている。本来視細胞層に血管は存在しないにもかかわらず、網膜表層で生じている血管影の明滅が視細胞層に映り込み、視細胞層の輝度値が変化することでアーチファクト８０９が生じる。
【0068】
補正部１０１－４３は、３次元の合成モーションコントラスト画像上に生じたプロジェクションアーチファクト８０２を抑制する処理を実行する。任意の公知のプロジェクションアーチファクト抑制手法を用いてよいが、本実施形態ではＳｔｅｐ－ｄｏｗｎ Ｅｘｐｏｎｅｎｔｉａｌ Ｆｉｌｔｅｒｉｎｇを用いる。Ｓｔｅｐ－ｄｏｗｎ Ｅｘｐｏｎｅｎｔｉａｌ Ｆｉｌｔｅｒｉｎｇでは、３次元モーションコントラスト画像上の各Ａスキャンデータに対して式（２）で表される処理を実行することにより、プロジェクションアーチファクトを抑制する。
【0069】
【数2】
【0070】
ここで、γは負の値を持つ減衰係数、Ｄ（ｘ，ｙ，ｚ）はプロジェクションアーチファクト抑制処理前の脱相関値、ＤＥ（ｘ，ｙ，ｚ）は該抑制処理後の脱相関値を表す。
【0071】
図８（ｄ）にプロジェクションアーチファクト抑制処理後の３次元合成モーションコントラストデータ（灰色）を断層画像上に重畳表示した例を示す。プロジェクションアーチファクト抑制処理前（図８（ｃ））に視細胞層上に見られたアーチファクトが、該抑制処理によって除去されたことがわかる。
【0072】
次に、投影部１０１－４５はＳ３０３で画像特徴取得部１０１－４４が取得した層境界の位置に基づく深度範囲でモーションコントラスト画像を投影し、正面モーションコントラスト画像を生成する。任意の深度範囲で投影してよいが、本実施形態においては網膜表層及び網膜深層の深度範囲で２種類の２次元合成モーションコントラスト画像を生成する。また、投影部１０１－４５は、投影法としては最大値投影（ＭＩＰ； Ｍａｘｉｍｕｍ Ｉｎｔｅｎｓｉｔｙ Ｐｒｏｊｅｃｔｉｏｎ）・平均値投影（ＡＩＰ； Ａｖｅｒａｇｅ Ｉｎｔｅｎｓｉｔｙ Ｐｒｏｊｅｃｔｉｏｎ）のいずれかを選択でき、本実施形態では最大値投影で投影するものとする。
【0073】
最後に、画像処理装置１０１は取得した画像群（ＳＬＯ画像や断層画像）と該画像群の撮影条件データや、生成したモーションコントラスト画像と付随する生成条件データを検査日時、披検眼を同定する情報と関連付けて外部記憶部１０２へ保存する。
【0074】
＜ステップ３０５＞
表示制御部１０１－０５は、Ｓ３０３で生成した断層画像や、Ｓ３０４で合成したモーションコントラスト画像、撮影条件や合成条件に関する情報を表示部１０４に表示させる。
表示制御部１０１－０５は、Ｓ３０３で生成した断層画像や、Ｓ３０４で合成したモーションコントラスト画像、撮影条件や合成条件に関する情報を表示部１０４に表示させる。
【0075】
図８（ｅ）にレポート画面８０３の例を示す。本実施形態ではＳＬＯ画像及び断層画像、Ｓ３０４で合成及び投影することにより生成した異なる深度範囲の正面モーションコントラスト画像、対応する正面ＯＣＴ画像を表示する。
【0076】
正面モーションコントラスト画像の投影範囲はリストボックスに表示された既定の深度範囲セット（８０５及び８０９）から操作者が選択することで変更できる。図８（ｅ）の例では、リストボックス８０５では網膜表層が選択され、リストボックス８０９では網膜深層が選択されている。８０４は網膜表層のＥｎＦａｃｅ画像、８０８は網膜深層のＥｎＦａｃｅ画像である。また、投影範囲の指定に用いる層境界の種類とオフセット位置を８０６や８１０のようなユーザインターフェースから変更したり、断層像上に重畳された層境界データ（８０７及び８１１）を入力部１０３から操作して移動させたりすることで投影範囲を変更できる。
【0077】
さらに、画像投影法やプロジェクションアーチファクト抑制処理の有無を例えばコンテキストメニューのようなユーザインターフェースから選択することにより変更してもよい。
【0078】
＜ステップ３０６＞
操作者が入力部１０３を用いてＯＣＴＡ計測処理の開始を指示する。
操作者が入力部１０３を用いてＯＣＴＡ計測処理の開始を指示する。
【0079】
本実施形態では図８（ｅ）のレポート画面８０３のモーションコントラスト画像上をダブルクリックすることで、図９（ａ）のようなＯＣＴＡ計測画面に移行する。モーションコントラスト画像が拡大表示され、適宜画像投影法の種類（最大値投影（ＭＩＰ）もしくは平均値投影（ＡＩＰ））や投影深度範囲、プロジェクションアーチファクト除去処理の実施の有無を選択する。次に、図９（ａ）右側のＭａｐボタン群９０２、Ｓｅｃｔｏｒボタン群９０３、Ｍｅａｓｕｒｅｍｅｎｔボタン９０４経由で表示される選択画面９０５、の中から操作者が適切な項目を選択することで計測の種類と対象領域を選択し、解析部１０１－４６が計測処理を開始する。なお、ＯＣＴＡ計測画面を表示させた時点では、いずれの計測対象領域も設定されていない（Ｍａｐボタン群９０２、Ｓｅｃｔｏｒボタン群ともＮｏｎｅが選択され、選択画面９０５が表示されていない状態）。解析部１０１－４６が計測処理を開始する。なお、ＯＣＴＡ計測画面を表示させた時点では、いずれの計測対象領域も設定されていない（Ｍａｐボタン群９０２、Ｓｅｃｔｏｒボタン群ともＮｏｎｅが選択され、選択画面９０５が表示されていない状態）。
【0080】
計測処理の種類として、本実施形態ではＭａｐボタン群９０２もしくはＳｅｃｔｏｒボタン群９０３から以下の（ｉ）～（ｉｉｉ）のいずれかを選択する。
（ｉ） Ｎｏｎｅ（計測しない）
（ｉｉ） ＶＡＤ（血管が占める面積に基づいて算出される血管密度）
（ｉｉｉ）ＶＬＤ（血管の長さの総和に基づいて算出される血管密度）
（ｉ） Ｎｏｎｅ（計測しない）
（ｉｉ） ＶＡＤ（血管が占める面積に基づいて算出される血管密度）
（ｉｉｉ）ＶＬＤ（血管の長さの総和に基づいて算出される血管密度）
【0081】
これに限らず、例えば血管構造の複雑さを定量化するＦｒａｃｔａｌ Ｄｉｍｅｎｓｉｏｎや、血管径の分布（血管の瘤や狭窄の分布）を表すＶｅｓｓｅｌ Ｄｉａｍｅｔｅｒ Ｉｎｄｅｘを選択可能に構成してもよい。また、Ｍｅａｓｕｒｅｍｅｎｔボタン９０４経由で表示される選択画面９０５からは以下の（ｉ）～（ｉｖ）のいずれかを選択できる。
（ｉ） 任意の２点間の距離計測
（ｉｉ）無血管領域の面積計測
（ｉｉｉ）ＶＡＤ
（ｉｖ）ＶＬＤ
（ｉ） 任意の２点間の距離計測
（ｉｉ）無血管領域の面積計測
（ｉｉｉ）ＶＡＤ
（ｉｖ）ＶＬＤ
【0082】
計測処理の対象領域として、本実施形態ではＭａｐボタン群９０２の中からＮｏｎｅ以外を選択することにより画像全体を設定でき、Ｓｅｃｔｏｒボタン群９０３の中からＮｏｎｅ以外を選択することによりセクタ領域（固視位置を中心とし、半径の異なる複数の同心円と該固視位置を通る角度の異なる複数の直線で規定される最小円領域及び扇形領域）を設定できる。またＭｅａｓｕｒｅｍｅｎｔボタン９０４経由で表示される選択画面９０５から所望の計測の種類を選択した上で入力部１０３を用いてモーションコントラスト画像上で任意形状の境界位置（図１０（ｂ）の灰色線部１００１）を指定してＯＫボタンを押下することで、任意形状の計測対象領域を設定できる。該領域内に示された数値は該領域内で計測した値（この場合はＶＡＤ値）を示している。なお、関心領域の手動設定時には指定された境界位置（灰色線部１００１）上に該境界位置が編集可能であることを示す丸い制御点が重畳表示され、ＯＫボタンを押下した時点で該丸い制御点は消えて灰色線部１００１のみとなり、該境界位置は編集不可の状態になる。
【0083】
本実施形態ではＭａｐボタン群９０２、Ｓｅｃｔｏｒボタン群９０３から各々ＶＡＤを選択することで、ＶＡＤマップ（計測の種類がＶＡＤ、計測対象領域が画像全体）とＶＡＤセクタマップ（計測の種類がＶＡＤ、計測対象領域がＥＴＤＲＳグリッドに対応するセクタ領域）を選択する場合について説明する。
【0084】
ここで、ＶＡＤはＶｅｓｓｅｌ Ａｒｅａ Ｄｅｎｓｉｔｙの略であり、計測対象に含まれる血管領域の割合で定義される血管密度（単位：％）である。また、ＶＬＤはＶｅｓｓｅｌ Ｌｅｎｇｔｈ Ｄｅｎｓｉｔｙ の略であり、単位面積あたりに含まれる血管の長さの総和（単位：ｍｍ－１）で定義される血管密度である。
【0085】
血管密度は血管の閉塞範囲や血管網の疎密の程度を定量化するための指標であり、ＶＡＤが最もよく用いられている。ただし、ＶＡＤでは計測値に占める大血管領域の寄与分が大きくなるため、糖尿病網膜症のように毛細血管の病態に注目して計測したい場合には（より毛細血管の閉塞に敏感な指標として）ＶＬＤが用いられる。
【0086】
なお、同一のモーションコントラスト画像に対して複数の計測対象領域を設定してもよい。複数の計測対象領域としては、例えば画像全体、セクタ領域、任意形状領域のうちの少なくとも２つや、２以上の深度範囲、あるいはそれらの組み合わせが挙げられる。該複数の計測対象領域に対して異なる種類の計測が選択された場合には、最後に（指定された解析対象領域に対して）選択された解析の種類を他の計測対象領域に対しても連動して適用した上で計測し、該計測結果を表示するように構成してもよい。例えば、ＶＡＤマップとＶＡＤセクタマップが選択された状態でＶＬＤマップに変更する指示を行うと、ＶＬＤセクタマップが自動で選択され、画像全体に対するＶＬＤ計測とＥＴＤＲＳセクタ領域に対するＶＬＤ計測が実行される。このような連動選択操作により、同一画像に対して異なる種類の計測値が重畳されて操作者が表示内容について誤解・混乱することを防止できる。
【0087】
なお、計測の種類のうちＮｏｎｅ（計測しない）に関しては、各計測対象領域で独立に選択する（ある計測対象領域に対して「Ｎｏｎｅ」が選択された場合に、該「Ｎｏｎｅ」の選択は他の計測対象領域に対して連動して適用されない）。また、本発明は最後に選択された解析の種類を全ての解析対象領域に対して連動して適用し、計測・表示することに限定されるものではない。例えば面内方向の複数の計測対象領域に対しては連動させて深度方向の複数の対象領域（網膜表層と網膜深層）に対しては連動させない、あるいはその逆（面内方向の複数の計測対象領域に対しては連動させず、深度方向の複数の対象領域に対しては連動させる）になるように適応的に適用して計測し、対応する計測結果を表示するよう構成してもよい。
【0088】
なお、ＶＡＤセクタマップおよびＶＬＤセクタマップは入力部１０３の指示に基づいて移動可能であり、移動に伴い値が計測部１０１－４６３により再計算される。
【0089】
次に、解析部１０１－４６は計測処理の前処理として画像拡大及びトップハットフィルタ処理を行う。トップハットフィルタを適用することにより、背景成分の輝度ムラを軽減できる。本実施形態では合成モーションコントラスト画像の画素サイズが約３μｍになるようにＢｉｃｕｂｉｃ補間を用いて画像拡大し、円形の構造要素を用いてトップハットフィルタ処理を行うものとする。
【0090】
＜ステップ３０７＞
解析部１０１－４６が血管領域の特定処理を行う。本実施形態では、強調部１０１－４６１がヘシアンフィルタ及びエッジ選択鮮鋭化に基づく血管強調処理を行う。次に抽出部１０１－４６２が２種類の血管強調画像を用いて２値化処理を行い、整形処理を行うことで血管領域を特定する。
解析部１０１－４６が血管領域の特定処理を行う。本実施形態では、強調部１０１－４６１がヘシアンフィルタ及びエッジ選択鮮鋭化に基づく血管強調処理を行う。次に抽出部１０１－４６２が２種類の血管強調画像を用いて２値化処理を行い、整形処理を行うことで血管領域を特定する。
【0091】
血管領域特定処理の詳細はＳ５１０～Ｓ５６０で説明する。
【0092】
＜ステップ３０８＞
計測部１０１－４６３が、操作者により指定された計測対象領域に関する情報に基づいて単検査の画像に対する血管密度の計測を行う。引き続いて表示制御部１０１－０５が、計測結果を表示部１０４に表示する。
計測部１０１－４６３が、操作者により指定された計測対象領域に関する情報に基づいて単検査の画像に対する血管密度の計測を行う。引き続いて表示制御部１０１－０５が、計測結果を表示部１０４に表示する。
【0093】
血管密度としてはＶＡＤとＶＬＤの２種類の指標があり、本実施形態ではＶＡＤを計算する場合の手順を例に説明する。なお、ＶＬＤを計算する場合の手順についても後述する。
【0094】
操作者が入力部１０３から血管領域もしくは血管中心線データを修正するよう指示を入力した場合には、解析部１０１－４６が操作者から入力部１０３を介して指定された位置情報に基づいて血管領域もしくは血管中心線データを修正し、計測値を再計算する。
【0095】
なお、本ステップで所定の条件を満たさずに計測が実施された場合には、表示制御部１０１－０５が所定の計測条件を満たした状態で計測すべき旨のメッセージ（警告表示）を表示部１０４に出力する。ここで、所定の条件とは、例えばＯＣＴＡ画像の重ね合わせが行われていることなどである。
【0096】
ＶＡＤ計測処理の詳細はＳ８１０～Ｓ８３０、ＶＬＤ計測処理の詳細はＳ８４０～Ｓ８７０で各々説明する。
【0097】
＜ステップ３０９＞
解析部１０１－４６は、Ｓ３０７で特定した血管領域や血管中心線のデータを修正するか否かの指示を外部から取得する。この指示は例えば入力部１０３を介して操作者により入力される。修正処理が指示された場合はＳ３０８へ、修正処理が指示されなかった場合はＳ３１０へと処理を進める。
解析部１０１－４６は、Ｓ３０７で特定した血管領域や血管中心線のデータを修正するか否かの指示を外部から取得する。この指示は例えば入力部１０３を介して操作者により入力される。修正処理が指示された場合はＳ３０８へ、修正処理が指示されなかった場合はＳ３１０へと処理を進める。
【0098】
＜ステップ３１０＞
比較部１０１－４６４は経時変化計測（Ｐｒｏｇｒｅｓｓｉｏｎ計測）処理を実行する。図１１にＰｒｏｇｒｅｓｓｉｏｎ計測レポートの例を示す。Ｐｒｏｇｒｅｓｓｉｏｎモードタブ１１０１を指定することにより、該レポート画面が表示され、Ｓ３０６で選択した計測の種類及び計測対象領域に基づく経時変化計測の処理が開始される。本実施形態ではＰｒｏｇｒｅｓｓｉｏｎ計測対象画像として、比較部１０１－４６４は、検査日が新しいものから順に４検査を自動選択するものとする。これに限らず、例えば最も古い検査日の画像と最新検査の画像、両検査日間に撮影された画像で検査間隔が略等しい画像を選択するよう構成してもよい。なお、最新検査とは例えばＳ３０２における撮影に関する検査である。
比較部１０１－４６４は経時変化計測（Ｐｒｏｇｒｅｓｓｉｏｎ計測）処理を実行する。図１１にＰｒｏｇｒｅｓｓｉｏｎ計測レポートの例を示す。Ｐｒｏｇｒｅｓｓｉｏｎモードタブ１１０１を指定することにより、該レポート画面が表示され、Ｓ３０６で選択した計測の種類及び計測対象領域に基づく経時変化計測の処理が開始される。本実施形態ではＰｒｏｇｒｅｓｓｉｏｎ計測対象画像として、比較部１０１－４６４は、検査日が新しいものから順に４検査を自動選択するものとする。これに限らず、例えば最も古い検査日の画像と最新検査の画像、両検査日間に撮影された画像で検査間隔が略等しい画像を選択するよう構成してもよい。なお、最新検査とは例えばＳ３０２における撮影に関する検査である。
【0099】
ここで、計測対象画像の選択条件は優先度が高い順に以下の（ｉ）、（ｉｉ）がある。なお、選択条件および優先度は下記の例に限定されるものではない。
（ｉ）固視位置が同じ画像であること。
（ｉｉ）略同一位置における取得断層像数が大きい（例えば４以上）モーションコントラスト画像か、または、同等のモーションコントラスト画像になるようにＯＣＴＡ重ね合わせ処理を実施して得られる合成モーションコントラスト画像であること。
（ｉ）固視位置が同じ画像であること。
（ｉｉ）略同一位置における取得断層像数が大きい（例えば４以上）モーションコントラスト画像か、または、同等のモーションコントラスト画像になるようにＯＣＴＡ重ね合わせ処理を実施して得られる合成モーションコントラスト画像であること。
【0100】
比較部１０１－４６４は、上記選択条件を満たす画像を優先して選択するものとする。例えば、固視位置が同一な最新の５検査のうち２番目に新しい検査の画像は重ね合わせがされていないＯＣＴＡ画像であり、他の検査の画像は重ね合わせがされたＯＣＴＡ画像である場合、比較部１０１－４６４は最新の検査および３番目～５番目に新しい検査を選択する。そして、表示制御部１０１－０５は、選択された検査のＯＣＴＡ画像または選択されたＯＣＴＡ画像から得た情報を表示部１０４に表示させる。すなわち、表示制御部１０１－０５は、合成されていないモーションコントラスト画像または合成されていないモーションコントラスト画像から得られた情報は時系列で表示させない。
【0101】
なお、表示制御部１０１－０５は、４つの画像を表示すべきところ選択条件を満たす画像が４未満の場合には、不足する分に関しては表示すべき画像がない旨を表示部１０４に表示させてもよいし、条件を満たしていないことを示す表示を付して画像を表示させることとしてもよい。例えば、２つの画像を表示しなければならないところ１つしか重ね合わせがされたＯＣＴＡ画像がない場合には、重ね合わせがされていないＯＣＴＡ画像を表示することとしてもよい。この場合、重ね合わせがされたＯＣＴＡ画像と重ね合わせがされていないＯＣＴＡ画像とを識別するための情報が表示されていてもよい。例えば、図１１において、ＯＣＴＡ画像が時系列で表示される場合には、重ね合わせＯＣＴＡ画像であることを示す情報（例えば「ＡＶＧ．」という表示）を画像毎（例えば画像の上部）に表示することとし、重ね合わせがされていないＯＣＴＡ画像には「ＡＶＧ．」を表示しないこととしてもよい。
【0102】
また、比較部１０１－４６４は、表示する候補として上記（ｉ）、（ｉｉ）を満たす画像が複数ある場合には、最新の画像を時系列で表示すべき画像として選択する。最新の画像にするメリットは最新の画像が医師等が好ましいと考えた画像であると推測されるからである。
【0103】
なお、計測対象画像はこれに限らず、例えば図１１の選択ボタン１１０７を選択して選択画面を表示させ、該選択画面上に表示された画像リストから選択してもよい。
【0104】
次に、比較部１０１－４６４はＳ３０９で実施した単検査の計測内容に対応する過去検査の画像及び計測値に関するデータを外部記憶部１０２から取得する。なお、計測値はあらかじめ計算されたものを外部記憶部１０２から取得することとしてもよいし、画像の取得後に改めて計算されることとしてもよい。さらに位置合わせ部１０１－４１はＳ３０８で計測した単検査の画像と過去検査画像間の位置合わせを行い、比較部１０１－４６４が共通した領域に関する計測データ（計測値、計測値マップ、差分マップ、トレンドグラフのうち少なくとも一つ）を生成する。なお、位置合わせはＯＣＴＡ画像を用いることとしてもよいし、ＳＬＯ画像を用いることとしてもよい。
【0105】
＜ステップ３１１＞
表示制御部１０１－０５は、Ｓ３１０で実施したＰｒｏｇｒｅｓｓｉｏｎ計測に関するレポートを表示部１０４に表示する。
表示制御部１０１－０５は、Ｓ３１０で実施したＰｒｏｇｒｅｓｓｉｏｎ計測に関するレポートを表示部１０４に表示する。
【0106】
本実施形態では、図１１に示すＰｒｏｇｒｅｓｓｉｏｎ計測レポートの上段に網膜表層において計測したＶＡＤマップとＶＡＤセクタマップを重畳表示し、下段に網膜深層において計測したＶＡＤマップとＶＡＤセクタマップを重畳表示する。これにより異なる深度位置の血管病態の時系列変化を一覧して把握でき、図１１の時系列並置表示されたＶＡＤ計測結果では、網膜深層に初期病変が生じ、時間の経過とともに網膜表層へ、また中心窩から傍中心窩へ血管閉塞が広がっていく様子が一覧して把握できる。なお、表示部１０４は、ＶＡＤマップとＶＡＤセクタマップに替えてＶＬＤマップとＶＬＤセクタマップとを表示することとしてもよい。また、表示部１０４はマップとセクタマップの何れか
一方のみを表示することとしてもよい。また、網膜表層と網膜深層との情報を上下に並べて表示するのではなく、網膜表層のＶＡＤに関する情報と網膜表層のＶＬＤに関する情報とを上下に並べて表示することとしてもよい。すなわち、 表示制御部１０１－０５は、合成画像から異なる計測手法で得られた複数の検査日に関する複数の情報を並置して表示部１０４に表示させてもよい。
一方のみを表示することとしてもよい。また、網膜表層と網膜深層との情報を上下に並べて表示するのではなく、網膜表層のＶＡＤに関する情報と網膜表層のＶＬＤに関する情報とを上下に並べて表示することとしてもよい。すなわち、 表示制御部１０１－０５は、合成画像から異なる計測手法で得られた複数の検査日に関する複数の情報を並置して表示部１０４に表示させてもよい。
【0107】
また、図１１において、表示制御部１０１－０５は、表示されている数値等がＶＡＤであることを示す表示を表示させることとしてもよい。例えば、表示制御部１０１－０５は、表示されている数値等がＶＡＤであることを示す表示としてＶＡＤの単位を表示部１０４に表示させることとしてもよいし、「ＶＡＤ」という文字を表示させることとしてもよい。なお、ＶＬＤについても表示されている数値等がＶＬＤであることを示す表示を表示させることとしてもよい。なお、ＶＡＤセクタマップを重ねる画像はＶＡＤマップではなく、重ね合わせされたＯＣＴＡ画像であってもよい。また、表示制御部１０１－０５は、ＶＡＤマップおよびＶＡＤセクタマップを表示することなく重ね合わせされたＯＣＴＡ画像のみを時系列に表示部１０４に表示させることとしてもよい。すなわち、表示制御部１０１－０５は、複数の検査日に関する複数の合成画像または複数の合成画像から得られた複数の検査日に関する複数の情報を時系列で表示手段に表示させる表示制御手段の一例に相当する。
【0108】
図１１に示す例では、表示制御部１０１－０５は、表示領域の上段に網膜表層のＶＡＤマップおよびＶＡＤセクタマップ、表示領域の下段に網膜深層のＶＡＤマップおよびＶＡＤセクタマップを表示部１０４に表示させている。網膜表層は第１の深度範囲の一例であり、網膜深層は第２の深度範囲の一例である。すなわち、表示制御部１０１－０５は、第１の深度範囲で投影した前記合成画像から得られた前記複数の検査日に関する複数の情報と、前記第１の深度範囲とは異なる第２の深度範囲で投影した前記合成画像から得られた前記複数の検査日に関する複数の情報とを並置して時系列で表示手段に表示させる。
【0109】
また各計測対象画像に関して、略同一位置における断層画像数やＯＣＴＡ重ね合わせ処理の実施条件に関する情報、ＯＣＴ断層画像もしくはモーションコントラスト画像の評価値（画質指標）に関する情報を表示部１０４に表示させてもよい。図１１ではＯＣＴＡ重ね合わせ処理実施済であることを示す印（左上の「Ａｖｅｒａｇｅｄ ＯＣＴＡ」）が表示されている。なお、図１１の上部に表示されている矢印１１０４は現在選択されている検査であることを示す印であり、基準検査（Ｂａｓｅｌｉｎｅ）はＦｏｌｌｏｗ－ｕｐ撮影の際に選択した検査（図１１の一番左画像の画像）である。もちろん、基準検査を示すマークを表示部１０４に表示させてもよい。またＳ３１０で「Ｓｈｏｗ Ｄｉｆｆｅｒｅｎｃｅ」チェックボックス１１０８が指定された場合には、基準画像上に基準画像に対する計測値分布（マップもしくはセクタマップ）を表示し、それ以外の検査日に対応する領域に基準画像に対して算出した計測値分布との差分計測値マップを表示する。計測結果としてはレポート画面上にトレンドグラフ（経時変化計測によって得られた各検査日の画像に対する計測値のグラフ）を表示させてもよい。該トレンドグラフの回帰直線（曲線）や対応する数式をレポート画面に表示させてもよい。ＶＡＤ、ＶＬＤ、無血管領域の大きさ等をトレンドグラフとして表示することが可能である。無血管領域の大きさとしてＦＡＺ（Ｆｏｖｅａｌ Ａｖａｓｃｕｌａｒ Ｚｏｎｅ）を表示してもよい。なお、トレンドグラフは、セクタマップの任意の領域の値を示したグラフであってもよく、トレンドグラフとして表示する領域は入力部１０３により切り替え可能としてもよい。また、トレンドグラフは、セクタマップの各領域のグラフを領域毎に識別可能な状態で同時に表示するものであってもよい。また、トレンドグラフとしてＶＡＤ，ＶＬＤを別々の座標系または同一の座標系に表示させることとしてもよい。このようにすればＶＡＤとＶＬＤとの関係をトレンドグラフから容易に把握可能となる。ここで、トレンドグラフは図１１に示したマップまたは重ね合わせされたＯＣＴＡ画像と同時に表示されることとしてもよいし、トレンドグラフ単独で表示されることとしてもよい。トレンドグラフに含まれる複数の計測値のそれぞれは、比較部１０１－４６４により選択された所定の基準を満たしている画像に関する値であるため、経時変化を正確に把握できる。
【0110】
なお、本実施形態では異なる深度範囲として網膜表層及び網膜深層の画像及び計測値を時系列で表示したが、これに限らず例えば網膜表層、網膜深層、網膜外層、脈絡膜の４種類の深度範囲の画像及び計測値を時系列で表示してもよい。また、表示部１０４は、任意の層の画像および計測値を時系列で表示することとしてもよい。
【0111】
あるいは、異なる指標の計測結果を並置して時系列で表示してもよい。例えば、上段にＶＡＤマップの時系列表示、下段にＶＬＤマップ（もしくは無血管領域の形状値）の時系列並置表示を行ってもよい。
【0112】
なお、時系列並置表示する場合の投影深度範囲は、Ｓ３０５で説明した図８（ｅ）のユーザインターフェース（８０５・８０６や８０９・８１０）の場合と同様に、図１１の１１０２、１１０３、１１０５、１１０６のユーザインターフェースを用いて変更できる。また、同様に投影方法（ＭＩＰ／ＡＩＰ）やプロジェクションアーチファクト抑制処理についても例えばコンテキストメニューから選択するなどの方法により変更してもよい。さらに、ショートカットメニューから計測の種類及び計測対象領域に関する項目を異なる値に変更することにより、Ｐｒｏｇｒｅｓｓｉｏｎ計測の種類及び計測対象領域を変更して再計測できる。
【0113】
例えば、ショートカットメニュー上に図９（ａ）のＭａｐボタン群９０２に相当する項目及びＳｅｃｔｏｒボタン群９０３に相当する項目を表示しておき、各々１項目ずつ（例えば「ＶＬＤ Ｍａｐ」と「ＶＬＤ Ｓｅｃｔｏｒ」を）選択する。なお、Ｓ３０６の場合と同様に、複数の計測対象領域が選択され、かつ一方の領域に対する計測の種類が変更された場合は他方の領域に対しても同じ計測の種類が連動して適用され、計測が実行される。また計測対象領域を設定しない（「Ｎｏｎｅ」を選択する）指示に関しては、他方の計測対象領域に対して連動して適用されないものとする。
【0114】
また、表示部１０４に表示されたモーションコントラスト画像や、抽出部１０１－４６２や計測部１０１－４６３によって生成された血管領域や血管中心線に関する２値画像、計測値や計測マップを外部記憶部１０２にファイルとして出力し、保存してもよい。なお、比較観察を容易にするため、ファイル出力するモーションコントラスト画像、血管領域や血管中心線に関する２値画像、計測値マップの画像サイズや画素サイズは同一にすることが望ましい。
【0115】
さらに、単検査に対する計測の場合と同様の方法で（詳細はＳ８３０で説明する）、推奨条件を満たさない状態で実施された計測結果を計測レポート画面に表示する際には警告メッセージを表示してよい。例えば、図１１における表示にＳ８３０で説明する警告メッセージを表示することとしてもよい。なお推奨条件としてはＳ８３０で示す条件に限定されるものではなく、例えば「選択した経時変化計測対象画像間で略同一位置における取得断層画像数もしくはモーションコントラスト画像の合成条件もしくは画質指標値の相違が所定値未満であること」を推奨条件として設定し、該条件を満たさない場合に警告表示してもよい。
【0116】
＜ステップ３１２＞
画像処理装置１０１はＳ３０１からＳ３１２に至る一連の処理を終了するか否かの指示を外部から取得する。この指示は入力部１０３を介して操作者により入力される。処理終了の指示を取得した場合は処理を終了する。一方、処理継続の指示を取得した場合にはＳ３０２に処理を戻し、次の披検眼に対する処理（または同一披検眼に対する再処理）を行う。
画像処理装置１０１はＳ３０１からＳ３１２に至る一連の処理を終了するか否かの指示を外部から取得する。この指示は入力部１０３を介して操作者により入力される。処理終了の指示を取得した場合は処理を終了する。一方、処理継続の指示を取得した場合にはＳ３０２に処理を戻し、次の披検眼に対する処理（または同一披検眼に対する再処理）を行う。
【0117】
さらに、図５（ａ）に示すフローチャートを参照しながら、Ｓ３０７で実行される処理の詳細について説明する。
【0118】
＜ステップ５１０＞
強調部１０１－４６１は、ステップＳ３０６の前処理を実施されたモーションコントラスト画像（ＯＣＴＡ画像）に対してヘッセ行列の固有値に基づく血管強調フィルタ処理を行う。このような強調フィルタはヘシアンフィルタと総称され、例えばＶｅｓｓｅｌｎｅｓｓ ｆｉｌｔｅｒやＭｕｌｔｉ－ｓｃａｌｅ ｌｉｎｅ ｆｉｌｔｅｒが挙げられる。本実施形態ではＭｕｌｔｉ－ｓｃａｌｅ ｌｉｎｅ ｆｉｌｔｅｒを用いるが、任意の公知の血管強調フィルタを用いてよい。
強調部１０１－４６１は、ステップＳ３０６の前処理を実施されたモーションコントラスト画像（ＯＣＴＡ画像）に対してヘッセ行列の固有値に基づく血管強調フィルタ処理を行う。このような強調フィルタはヘシアンフィルタと総称され、例えばＶｅｓｓｅｌｎｅｓｓ ｆｉｌｔｅｒやＭｕｌｔｉ－ｓｃａｌｅ ｌｉｎｅ ｆｉｌｔｅｒが挙げられる。本実施形態ではＭｕｌｔｉ－ｓｃａｌｅ ｌｉｎｅ ｆｉｌｔｅｒを用いるが、任意の公知の血管強調フィルタを用いてよい。
【0119】
ヘシアンフィルタは強調したい血管の径に適したサイズで画像を平滑化した上で、該平滑化画像の各画素において輝度値の２次微分値を要素として持つヘッセ行列を算出し、該行列の固有値の大小関係に基づいて局所構造を強調する。ヘッセ行列は式（３）で与えられるような正方行列であり、該行列の各要素は例えば式（４）に示すような画像の輝度値Ｉを平滑化した画像の輝度値Ｉｓの２次微分値で表される。ヘシアンフィルタでは、このようなヘッセ行列の「固有値（λ１、λ２）の一方が０に近く、他方が負かつ絶対値が大きい」場合に線状構造とみなして強調する。これはモーションコントラスト画像上の血管領域が持つ特徴、すなわち「走行方向では輝度変化が小さく、走行方向に直交する方向では輝度値が大きく低下する」が成り立つ画素を線状構造とみなして強調することに相当する。すなわち、ヘシアンフィルタは線状構造強調フィルタの一例に相当する。
【0120】
またモーションコントラスト画像には毛細血管から細動静脈まで様々な径の血管が含まれることから、複数のスケールでガウスフィルタにより平滑化した画像に対してヘッセ行列を用いて線強調画像を生成する。例えば、毛細血管の血管径に対応するスケールおよび視神経乳頭付近の血管径に対応するスケールを用いることとしてもよい。次に式（５）に示すようにガウスフィルタの平滑化パラメータσの二乗を補正係数として乗じた上で最大値演算により合成し、該合成画像Ｉｈｅｓｓｉａｎをヘシアンフィルタの出力とする。
【0121】
ヘシアンフィルタはノイズに強く、血管の連続性が向上するという利点がある。一方で実際には事前に画像に含まれる最大血管径が不明の場合が多いため、特に平滑化パラメータが画像中の最大血管径に対して大きすぎる場合に強調された血管領域が太くなりやすいという欠点がある。
【0122】
そこで、本実施形態ではＳ５３０で述べる別の血管強調手法で血管領域を強調した画像と演算することにより、血管領域が太くなりすぎることを防止する。
【0123】
【数3】
【0124】
【数4】
【0125】
【数5】
【0126】
＜ステップ５２０＞
抽出部１０１－４６２は、Ｓ５１０で生成したヘシアンフィルタによる血管強調画像（以下、ヘシアン強調画像と表記）を２値化する。
抽出部１０１－４６２は、Ｓ５１０で生成したヘシアンフィルタによる血管強調画像（以下、ヘシアン強調画像と表記）を２値化する。
【0127】
ヘシアン強調画像の輝度統計値（平均値や中央値等）を閾値として２値化する場合、例えば視神経乳頭部では大血管の高輝度領域の影響を受けて閾値が高くなり、乳頭周囲毛細血管（ＲＰＣ； Ｒａｄｉａｌ Ｐｅｒｉｐａｐｉｌｌａｒｙ Ｃａｐｉｌｌａｒｙ）の抽出不足が生じる場合がある。また、網膜深層のような無血管領域の拡大が生じやすい領域においては、閾値が低すぎて無血管領域を血管として誤検出する場合がある。
【0128】
そこで、本実施形態では低スケール（複数のフィルタのスケールのうち閾値以下の一部のスケール）の強調画像のみで合成したヘシアン強調画像の平均値を閾値とすることで視神経乳頭部において閾値が高くなりすぎることを防ぐ。また該閾値の下限値を設定することにより無血管領域における誤検出を抑制する。
【0129】
ここで、視神経乳頭部で閾値が高くなりすぎるのを防止する方法としては、低スケールの強調画像の統計値を閾値として２値化することに限定されない。例えば、ヘシアン強調画像上の輝度値が所定値以上である場合に所定値であるとみなした場合の平均値を閾値として２値化しても同等の効果が期待できる。あるいは、例えばＭ－ｅｓｔｉｍａｔｏｒのようなロバスト推定値を閾値として２値化してもよい。
【0130】
また、合成モーションコントラスト画像をヘシアンフィルタで強調しているため、単独のモーションコントラスト画像をヘシアンフィルタで強調する場合に比べて２値化した血管領域の連続性がさらに向上する。
【0131】
＜ステップ５３０＞
強調部１０１－４６１は、Ｓ３０６で生成したトップハットフィルタ適用後の合成モーションコントラスト画像に対してエッジ選択鮮鋭化処理を行う。
強調部１０１－４６１は、Ｓ３０６で生成したトップハットフィルタ適用後の合成モーションコントラスト画像に対してエッジ選択鮮鋭化処理を行う。
【0132】
ここで、エッジ選択鮮鋭化処理とは画像のエッジ部分（輝度差が大きい部分）に重みを大きく設定した上で重みづけ鮮鋭化処理を行うことを指す。本実施形態では、前記合成モーションコントラスト画像に対してＳｏｂｅｌフィルタを適用した画像を重みとしてアンシャープマスク処理を行うことにより、エッジ選択鮮鋭化処理を実施する。
【0133】
小さなフィルタサイズで鮮鋭化処理を実施すると、細い血管のエッジが強調され２値化した際により正確に血管領域を特定できる（血管領域が太くなる現象を防止できる）。一方で特に同一撮影位置での断層像数が少ないモーションコントラスト画像の場合にはノイズが多いため、特に血管内のノイズも一緒に強調してしまう恐れがある。そこで、エッジ選択鮮鋭化を行うことによってノイズの強調を抑制する。
【0134】
＜ステップ５４０＞
抽出部１０１－４６２は、Ｓ５３０で生成したエッジ選択鮮鋭化処理を適用した鮮鋭化画像を２値化する。任意の公知の２値化法を用いてよいが、本実施形態では該鮮鋭化画像上の各局所領域内で算出した輝度統計値（平均値もしくは中央値）を閾値として２値化する。
抽出部１０１－４６２は、Ｓ５３０で生成したエッジ選択鮮鋭化処理を適用した鮮鋭化画像を２値化する。任意の公知の２値化法を用いてよいが、本実施形態では該鮮鋭化画像上の各局所領域内で算出した輝度統計値（平均値もしくは中央値）を閾値として２値化する。
【0135】
ただし、視神経乳頭部の大血管領域においては設定される閾値が高すぎて２値画像上の血管領域内に多数の穴が空くため、上記閾値の上限値を設定することにより特に視神経乳頭部において閾値が高くなりすぎるのを防止する。
【0136】
また、Ｓ５２０の場合と同様に画像中に無血管領域が占める割合が大きい場合には閾値が低すぎて無血管領域の一部を血管として誤検出する場合が生じる。そこで、上記閾値の下限値を設定することにより誤検出を抑制する。
【0137】
なお、Ｓ５２０の場合と同様に合成モーションコントラスト画像をエッジ選択鮮鋭化しているため、単独のモーションコントラスト画像をエッジ選択鮮鋭化する場合に比べて２値化した場合のノイズ状の誤検出領域をより減らすことができる。
【0138】
＜ステップ５５０＞
抽出部１０１－４６２は、Ｓ５２０で生成したヘシアン強調画像の２値画像の輝度値と、Ｓ５４０で生成したエッジ選択鮮鋭化画像の２値画像の輝度値の双方が０より大きい場合に血管候補領域として抽出（セグメンテーション）する。当該演算処理により、ヘシアン強調画像に見られる血管径を過大評価している領域と、エッジ選択鮮鋭化画像上に見られるノイズ領域がともに抑制され、血管の境界位置が正確かつ血管の連続性が良好な２値画像を取得できる。
抽出部１０１－４６２は、Ｓ５２０で生成したヘシアン強調画像の２値画像の輝度値と、Ｓ５４０で生成したエッジ選択鮮鋭化画像の２値画像の輝度値の双方が０より大きい場合に血管候補領域として抽出（セグメンテーション）する。当該演算処理により、ヘシアン強調画像に見られる血管径を過大評価している領域と、エッジ選択鮮鋭化画像上に見られるノイズ領域がともに抑制され、血管の境界位置が正確かつ血管の連続性が良好な２値画像を取得できる。
【0139】
また双方の２値画像とも合成モーションコントラスト画像に基づく２値画像であることから、単独のモーションコントラスト画像に基づく２値画像に比べて２値化した場合のノイズ状の誤検出領域が減少するとともに、特に毛細血管領域の連続性が向上する。また、合成モーションコントラスト画像であることから検査間の画質や輝度レベルも安定し、血管の抽出能が検査間で安定しやすい。
【0140】
＜ステップ５６０＞
抽出部１０１－４６２は、血管領域の整形処理として２値画像のオープニング処理（収縮処理後に膨張処理を行うこと）及びクロージング処理（膨張処理後に収縮処理を行うこと）を実施する。なお、整形処理はこれに限らず例えば２値画像をラベリングした場合の各ラベルの面積に基づく小領域除去を行ってもよい。なお、本処理は必須の処理ではない。
抽出部１０１－４６２は、血管領域の整形処理として２値画像のオープニング処理（収縮処理後に膨張処理を行うこと）及びクロージング処理（膨張処理後に収縮処理を行うこと）を実施する。なお、整形処理はこれに限らず例えば２値画像をラベリングした場合の各ラベルの面積に基づく小領域除去を行ってもよい。なお、本処理は必須の処理ではない。
【0141】
上記のＳ５１０～Ｓ５６０により血管領域の２値画像を得ることができる。２値画像は血管と血管以外の部位に異なるラベルが付けられた画像でありセグメンテーションの結果ということができる。
【0142】
なお、様々な径の血管が含まれるモーションコントラスト画像において血管抽出する方法はＳ５１０～Ｓ５６０に述べた方法に限定されない。例えば、図５（ｂ）のＳ６１０～Ｓ６５０に示すように、ヘシアン強調画像の輝度値とエッジ選択鮮鋭化による血管強調画像の輝度値を乗じる演算（Ｓ６３０）を適用した画像に対する輝度統計値（例えば平均値）を閾値として２値化（Ｓ６４０）することによって血管領域を特定してもよい。該閾値には下限値や上限値を設定できる。なお、Ｓ６１０，６２０は、Ｓ５１０，５３０と同様の処理であり、Ｓ６５０はＳ５６０と同様の処理である。
【0143】
あるいは、図５（ｃ）のＳ７１０～Ｓ７４０に示すように、画像の固視位置や深度範囲によってヘシアンフィルタ適用時の平滑化パラメータσの範囲を適応的に変えた上で（Ｓ７１０）、ヘシアンフィルタを適用し（Ｓ７２０）、２値化（Ｓ７３０）することによって血管強調してもよい。なお、Ｓ７４０はＳ５６０と同様の処理である。例えば、乳頭部網膜表層ではσ＝１～１０、黄斑部網膜表層ではσ＝１～８、黄斑部網膜深層ではσ＝１～６のように撮影部位に応じてスケールを設定できる。
【0144】
また、２値化処理は閾値処理に限定されるものではなく任意の公知のセグメンテーション手法を用いてよい。
【0145】
さらに、図６（ａ）に示すフローチャートを参照しながら、Ｓ３０８で実行される処理の詳細について説明する。
【0146】
＜ステップ８１０＞
操作者は、入力部１０３を介して計測処理における関心領域を設定する。本実施形態ではＳ３０６で計測内容（計測の種類と対象領域）としてＶＡＤマップ（計測の種類がＶＡＤで、計測対象領域が画像全体）とＶＡＤセクタマップ（計測の種類がＶＡＤで、計測対象領域がＥＴＤＲＳグリッドに対応するセクタ領域）が選択されている。従って、関心領域として（ｉ）画像全体（ｉｉ）固視灯位置を中心とするセクタ領域（直径１ｍｍの内円と直径３ｍｍの外円で規定される環状領域内をＳｕｅｒｉｏｒ・Ｉｎｆｅｒｉｏｒ・Ｎａｓａｌ・Ｔｅｍｐｏｒａｌの４つの扇形に分割した領域及び該内円領域）を設定する。
操作者は、入力部１０３を介して計測処理における関心領域を設定する。本実施形態ではＳ３０６で計測内容（計測の種類と対象領域）としてＶＡＤマップ（計測の種類がＶＡＤで、計測対象領域が画像全体）とＶＡＤセクタマップ（計測の種類がＶＡＤで、計測対象領域がＥＴＤＲＳグリッドに対応するセクタ領域）が選択されている。従って、関心領域として（ｉ）画像全体（ｉｉ）固視灯位置を中心とするセクタ領域（直径１ｍｍの内円と直径３ｍｍの外円で規定される環状領域内をＳｕｅｒｉｏｒ・Ｉｎｆｅｒｉｏｒ・Ｎａｓａｌ・Ｔｅｍｐｏｒａｌの４つの扇形に分割した領域及び該内円領域）を設定する。
【0147】
＜ステップ８２０＞
計測部１０１－４６３は、Ｓ３０７で得られた血管領域の２値画像に基づいて計測処理を行う。本実施形態では、該２値画像の各画素位置において当該画素を中心とした近傍領域内に占める非０画素（白画素）の割合を当該画素における血管密度（ＶＡＤ）として算出する。さらに、各画素で算出した血管密度（ＶＡＤ）の値を持つ画像（ＶＡＤマップ）を生成する。
計測部１０１－４６３は、Ｓ３０７で得られた血管領域の２値画像に基づいて計測処理を行う。本実施形態では、該２値画像の各画素位置において当該画素を中心とした近傍領域内に占める非０画素（白画素）の割合を当該画素における血管密度（ＶＡＤ）として算出する。さらに、各画素で算出した血管密度（ＶＡＤ）の値を持つ画像（ＶＡＤマップ）を生成する。
【0148】
また、該２値画像上の（Ｓ８１０で設定した）各セクタ領域における非０画素（白画素）の割合を当該セクタにおける血管密度（ＶＡＤ）として算出する。さらに、各セクタ領域で算出した血管密度（ＶＡＤ）の値を持つマップ（ＶＡＤセクタマップ）を生成する。
【0149】
＜ステップ８３０＞
表示制御部１０１－０５は、計測結果としてＳ８２０で生成したＶＡＤマップ及びＶＡＤセクタマップを表示部１０４に表示する。本実施形態では、図９（ｂ）の９０６に網膜表層のＶＡＤマップ、９０８に網膜深層のＶＡＤマップを表示する。さらに９０７に網膜表層のＶＡＤセクタマップ、９０９に網膜深層のＶＡＤセクタマップを重畳表示する。
表示制御部１０１－０５は、計測結果としてＳ８２０で生成したＶＡＤマップ及びＶＡＤセクタマップを表示部１０４に表示する。本実施形態では、図９（ｂ）の９０６に網膜表層のＶＡＤマップ、９０８に網膜深層のＶＡＤマップを表示する。さらに９０７に網膜表層のＶＡＤセクタマップ、９０９に網膜深層のＶＡＤセクタマップを重畳表示する。
【0150】
なお、本実施形態では図９（ｂ）において実施する計測の推奨条件として（ｉ）～（ｉｖ）を設定するものとする。なお、推奨条件は（ｉ）～（ｉｖ）のうち少なくとも１つを用いればよく、下記の条件に限定されるものではない。
（ｉ） 選択した計測対象画像の中に略同一位置における取得断層画像数が所定値以上のモーションコントラスト画像であること。もしくは該所定値以上相当の合成モーションコントラスト画像であること。あるいは画質指標値（Ｑｕａｌｉｔｙ Ｉｎｄｅｘ）が所定値以上のモーションコントラスト画像に対する計測であること。
（ｉｉ） 最大値投影で生成されたモーションコントラスト画像に対する計測であること。
（ｉｉｉ） プロジェクションアーチファクト除去処理（ＰＡＲ）が実施済であること。
（ｉｖ） 網膜表層、網膜深層、放射状乳頭周囲毛細血管（ＲＰＣ）が含まれる投影深度範囲、のいずれか投影深度範囲で生成されたモーションコントラスト画像に対する計測であること。
表示制御部１０１－０５は、上記（ｉ）～（ｉｖ）の少なくとも一つを満たさない状態で実施された計測結果を計測レポート画面に表示する際には正確に計測することが期待できないような条件で計測が行われたとみなして警告表示を行う。
（ｉ） 選択した計測対象画像の中に略同一位置における取得断層画像数が所定値以上のモーションコントラスト画像であること。もしくは該所定値以上相当の合成モーションコントラスト画像であること。あるいは画質指標値（Ｑｕａｌｉｔｙ Ｉｎｄｅｘ）が所定値以上のモーションコントラスト画像に対する計測であること。
（ｉｉ） 最大値投影で生成されたモーションコントラスト画像に対する計測であること。
（ｉｉｉ） プロジェクションアーチファクト除去処理（ＰＡＲ）が実施済であること。
（ｉｖ） 網膜表層、網膜深層、放射状乳頭周囲毛細血管（ＲＰＣ）が含まれる投影深度範囲、のいずれか投影深度範囲で生成されたモーションコントラスト画像に対する計測であること。
表示制御部１０１－０５は、上記（ｉ）～（ｉｖ）の少なくとも一つを満たさない状態で実施された計測結果を計測レポート画面に表示する際には正確に計測することが期待できないような条件で計測が行われたとみなして警告表示を行う。
【0151】
例えば（ｉ）を満たさない状態で計測された結果を表示する際には、表示制御部１０１－０５は、図９（ｂ）の例えば右下部に「Ａｖｅｒａｇｅｄ ＯＣＴＡ ｉｓ ｒｅｃｏｍｍｅｎｄｅｄ ｉｎ ｃａｌｃｕｌａｔｉｎｇ ＶＡＤ ｏｒ ＶＬＤ．」のように警告メッセージを表示部１０４に表示させればよい。
【0152】
また（ｉｉ）を満たさない状態で計測された結果を表示する際には、表示制御部１０１－０５は、図９（ｂ）の例えば右下部に「ＭＩＰ ｉｓ ｒｅｃｏｍｍｅｎｄｅｄ ｉｎ ｃａｌｃｕｌａｔｉｎｇ ＶＡＤ ｏｒ ＶＬＤ．」のように警告メッセージを表示部１０４に表示させればよい。
【0153】
同様に、（ｉｉｉ）を満たさない状態で計測された結果を表示する際には、表示制御部１０１－０５は、図９（ｂ）の例えば右下部に「ＰＡＲ ｉｓ ｒｅｃｏｍｍｅｎｄｅｄ ｉｎ ｃａｌｃｕｌａｔｉｎｇ ＶＡＤ ｏｒ ＶＬＤ．」のように警告メッセージを表示部１０４に表示させればよい。
【0154】
さらに、（ｉｖ）を満たさない状態で計測された結果を表示する際には表示制御部１０１－０５は、「Ｓｕｐｅｒｆｉｃｉａｌ Ｃａｐｉｌｌａｒｙ， Ｄｅｅｐ Ｃａｐｉｌｌａｒｙ， ＲＰＣ ｃａｎ ｂｅ ａｎａｌｙｚｅｄ ｉｎ ｃａｌｃｕｌａｔｉｎｇ ＶＡＤ ｏｒ ＶＬＤ．」のように警告メッセージを表示部１０４に表示させる。該警告メッセージを表示することで、推奨する計測条件を満たさない計測により得られた計測結果は信頼性の低い計測結果となるリスクがあることを周知するとともに、推奨する計測条件を示すことで、より信頼性の高い計測を実施しやすくする。
【0155】
なお、レポート画面を多数の警告メッセージが占めるのを避けるために上記推奨条件に優先順位を付け（例えば（ｉ）を最優先、（ｉｉ）を２番目、（ｉｉｉ）を３番目、（ｉｖ）を４番目に重要な条件とし）満たされていない計測条件のうち最も優先度の高い条件に関する警告表示を行うよう構成してもよい。また図９（ｂ）のように複数の計測結果を表示する場合には、個別の計測に対して各々警告メッセージを表示してもよいし、表示対象の警告メッセージのうち最も優先度の高い警告メッセージのみ表示させるようにしてもよい。あるいは、表示制御部１０１－０５は、計測結果の信頼性に影響が大きい条件を理解しやすくしつつ満たされていない条件に関する警告を漏れなく表示するために、満たされていない計測条件に関する警告を優先順位が識別可能な態様で（色や大きさ等を変えて）表示部１０４に表示させてもよい。複数の計測結果を表示する場合の例として、レポート画面の上下段に各々計測結果を表示する場合や、同一画像に対し複数の計測対象領域を設定して計測した結果を表示する場合が挙げられる。
【0156】
なお、警告メッセージは同一のレポート画面内に表示してもよいし、別画面として表示してもよい。また警告メッセージは文字列に限定されるものではなく静止画像や動画像を表示部１０４に表示したり、音声として出力したりしてもよい。該警告メッセージが表示されたレポート画面をファイル出力したり、印刷出力したりする場合も本発明に含まれる。
【0157】
さらに、表示部１０４に表示された警告メッセージのうち消去する警告メッセージを操作者が入力部１０３を用いて選択可能にしたり、警告の優先順位の変更や表示対象外とする警告メッセージの指定を行うためのユーザインターフェースを備えたりしてもよい。
【0158】
上記では血管密度としてＶＡＤを計測する場合の手順を例に説明したが、計測値としてＶＬＤマップやＶＬＤセクタマップを生成する場合は、上記Ｓ８１０～８３０の代わりに図６（ｂ）に示すＳ８４０～８７０を実行する。
【0159】
＜ステップ８４０＞
計測部１０１－４６３は、Ｓ３０７で生成した血管領域の２値画像を細線化処理することにより、血管の中心線に相当する線幅１画素の２値画像（以下、スケルトン画像と表記）を生成する。任意の細線化法もしくはスケルトン処理を用いてよいが、本実施形態では細線化法としてＨｉｌｄｉｔｃｈの細線化法を用いる。
計測部１０１－４６３は、Ｓ３０７で生成した血管領域の２値画像を細線化処理することにより、血管の中心線に相当する線幅１画素の２値画像（以下、スケルトン画像と表記）を生成する。任意の細線化法もしくはスケルトン処理を用いてよいが、本実施形態では細線化法としてＨｉｌｄｉｔｃｈの細線化法を用いる。
【0160】
＜ステップ８５０＞
操作者は、入力部１０３を介してＳ８１０と場合と同様の関心領域を設定する。本実施形態では計測内容（計測の種類と対象領域）としてＶＬＤマップとＶＬＤセクタマップを算出するものとし、Ｓ８１０ではＶＡＤが選択されていたのに対し、本ステップではＶＬＤが選択されていることが異なるだけである。なお、ＶＬＤマップもしくはＶＬＤセクタマップをモーションコントラスト画像上に重畳表示したくない場合は、図９（ａ）のＭａｐもしくはＳｅｃｔｏｒの項目を「Ｎｏｎｅ」に設定すればよい。
操作者は、入力部１０３を介してＳ８１０と場合と同様の関心領域を設定する。本実施形態では計測内容（計測の種類と対象領域）としてＶＬＤマップとＶＬＤセクタマップを算出するものとし、Ｓ８１０ではＶＡＤが選択されていたのに対し、本ステップではＶＬＤが選択されていることが異なるだけである。なお、ＶＬＤマップもしくはＶＬＤセクタマップをモーションコントラスト画像上に重畳表示したくない場合は、図９（ａ）のＭａｐもしくはＳｅｃｔｏｒの項目を「Ｎｏｎｅ」に設定すればよい。
【0161】
＜ステップ８６０＞
計測部１０１－４６３はＳ８４０で得られたスケルトン画像に基づいて計測処理を行う。本実施形態では、該スケルトン画像の各画素位置において当該画素を中心とした近傍領域における単位面積当たりの非０画素（白画素）の長さの総和［ｍｍ－１］を当該画素における血管密度（ＶＬＤ）として算出する。さらに、各画素で算出した血管密度（ＶＬＤ）の値を持つ画像（ＶＬＤマップ）を生成する。
計測部１０１－４６３はＳ８４０で得られたスケルトン画像に基づいて計測処理を行う。本実施形態では、該スケルトン画像の各画素位置において当該画素を中心とした近傍領域における単位面積当たりの非０画素（白画素）の長さの総和［ｍｍ－１］を当該画素における血管密度（ＶＬＤ）として算出する。さらに、各画素で算出した血管密度（ＶＬＤ）の値を持つ画像（ＶＬＤマップ）を生成する。
【0162】
また、該スケルトン画像上の（Ｓ８５０で設定した）各セクタ領域における単位面積当たりの非０画素（白画素）の長さの総和［ｍｍ－１］を当該セクタにおける血管密度（ＶＬＤ）として算出する。さらに、各セクタ領域で算出した血管密度（ＶＬＤ）の値を持つマップ（ＶＬＤセクタマップ）を生成する。
【0163】
＜ステップ８７０＞
表示制御部１０１－０５は、計測結果としてＳ８６０で生成したＶＬＤマップ及びＶＬＤセクタマップを図９（ｂ）の９０６・９０７もしくは９０８・９０９に表示する。
表示制御部１０１－０５は、計測結果としてＳ８６０で生成したＶＬＤマップ及びＶＬＤセクタマップを図９（ｂ）の９０６・９０７もしくは９０８・９０９に表示する。
【0164】
なおＳ８３０の場合と同様に、所定の解析に適した条件を満たさない状態で実施された計測結果を計測レポート画面に表示する際には表示部１０４に警告メッセージを表示するものとする。
【0165】
また本実施形態では単検査での血管領域特定及び計測結果の表示法として正面モーションコントラスト画像上に計測マップを重畳表示する場合について説明したが、これに限定されるものではない。例えば、特定された血管領域の２値画像やスケルトン画像を図９（ｂ）の９０６や９０８に表示させてもよい。あるいは、９０６や９０８にモーションコントラスト画像を表示しておき、その上に特定された血管領域の２値画像もしくはスケルトン画像を色もしくは透明度を適宜調整した上で重畳表示するよう構成する場合も本発明に含まれる。また２値画像は正面画像として表示することに限定されず、例えばＢスキャン断層画像上に特定された血管領域の２値画像もしくはスケルトン画像を色もしくは透明度を適宜調整した上で重畳表示してもよい。
【0166】
なお、Ｓ３０９で操作者が入力部１０３から血管領域もしくは血管中心線データを修正するよう指示を入力した場合には以下のような手順で修正する。
【0167】
すなわち、図１０（ａ）に示すような合成モーションコントラスト画像に対して同図（ｃ）に示すような過抽出領域を含む２値画像が得られた場合に、操作者が入力部１０３を経由して指定した位置の白画素を解析部１０１－４６が削除する。追加／削除位置の指定方法の例としては、例えば削除の場合は「ｄ」キーを押しながら、追加の場合は「ａ」キーを押しながらマウスでクリックする方法が挙げられる。あるいは、同図（ｄ）に示すようにモーションコントラスト画像に基づく画像上に修正対象である２値画像（血管領域もしくは血管中心線）の色や透明度を調整して重畳表示し、過抽出もしくは抽出不足の領域が判別しやすい状態にしておく。同図（ｄ）の矩形領域１００２内を拡大した画像を同図（ｅ）に示す。灰色が過抽出した領域で、白色が元のモーションコントラスト画像の脱相関値を示す。該過抽出／抽出不足の領域を操作者が入力部１０３を用いて指定することにより、正確かつ効率的に２値画像上の血管もしく血管中心線領域を修正するよう構成してもよい。なお、２値画像の修正処理は正面画像に限定されない。例えば図９（ａ）の９１０に示すような任意のスライス位置のＢスキャン断層像上にモーションコントラストデータや血管領域の２値データもしくは血管中心線領域を色や透明度の調整後に重畳する。このように過抽出もしくは抽出不足の領域が判別しやすい状態にした上で、操作者が修正（追加／削除）する２値データの３次元位置（ｘ，ｙ，ｚ座標）を入力部１０３により指定して修正してもよい。
【0168】
さらに、２値画像（血管領域の２値画像もしくはスケルトン画像）が修正済であることを示す情報もしくは修正位置に関する情報を該２値画像と関連付けて外部記憶部１０２に保存しておき、Ｓ８７０もしくはＳ３１１で血管特定結果及び計測結果を表示部１０４に表示する際に該修正済であることを示す情報もしくは修正位置に関する情報を表示部１０４に表示してもよい。
【0169】
また、本実施形態では合成部１０１－４２が繰り返しＯＣＴＡ撮影終了時に合成モーションコントラスト画像を生成する場合について説明したが、合成モーションコントラスト画像の生成手順はこれに限定されない。例えば図８（ｅ）のレポート画面８０３上に合成モーションコントラスト画像生成指示ボタン８１２を配置しておく。ＯＣＴＡ撮影完了後（撮影日より後の日でもよい）に操作者が明示的に該生成指示ボタン８１２を押下した場合に合成部１０１－４２が合成モーションコントラスト画像を生成するよう画像処理装置１０１を構成してもよい。操作者が明示的に合成画像生成指示ボタン８１２を押下して合成画像を生成する場合、図８（ｅ）に示すようなレポート画面８０３上に合成モーションコントラスト画像８０４や合成条件データ、検査画像リスト上に合成画像に関する項目を表示させる。
【0170】
また操作者が明示的に該生成指示ボタン８１２を押下する場合は、表示制御部１０１－０５が以下の処理を行う。すなわち、合成対象画像選択画面を表示させ、操作者が入力部１０３を操作して合成対象画像群を指定し、許容ボタンを押下した場合に合成部１０１－４２が合成モーションコントラスト画像を生成し、表示部１０４に表示させる。なお、生成済の合成モーションコントラスト画像を選択して合成する場合も本発明に含まれる。
【0171】
また操作者が合成画像生成指示ボタン８１２を押下した場合、３次元モーションコントラスト画像を投影した２次元画像同士を合成することにより２次元合成画像を生成してもよいし、３次元合成画像を生成後に投影することで２次元合成画像を生成してもよい。
【0172】
以上述べた構成によれば、画像処理装置１０１は同一被検眼に対して異なる日時に略同一撮影条件で取得したＯＣＴＡ重ね合わせ画像から生成した網膜表層及び網膜深層の正面モーションコントラスト画像を用いて血管領域特定及び血管密度計測処理を実施する。画像処理装置１０１は、該特定処理及び計測処理により得られた合成画像及び計測値を複数の深度範囲で時系列に並べて表示する。
【0173】
ＯＣＴＡ重ね合わせ画像を用いているため、検査ごとのＯＣＴ断層画像の信号強度や画質のばらつきの影響を抑制することが可能となる。この結果、本実施形態によれば、眼底血管に関する経時変化の適切な評価を支援することができる。具体的には、本実施形態によれば、血管病態の変化を正確に特定・計測できる。
【0174】
［第２の実施形態］
本実施形態に係る画像処理装置は、第１実施形態における血管領域特定及び計測処理を３次元で実施し、得られた画像及び計測データ（血管領域や血管中心線・計測データ）を時系列で並置表示するよう構成したものである。
本実施形態に係る画像処理装置は、第１実施形態における血管領域特定及び計測処理を３次元で実施し、得られた画像及び計測データ（血管領域や血管中心線・計測データ）を時系列で並置表示するよう構成したものである。
【0175】
具体的には、脈絡膜新生血管（ＣＮＶ；Ｃｈｏｒｏｉｄａｌ ＮｅｏＶａｓｕｃｕｌａｒｉｚａｔｉｏｎ）を含む３次元の合成モーションコントラスト画像に対してモーションアーチファクト抑制処理を実施する。次に３次元のモルフォロジーフィルタ及び血管強調フィルタを適用し、２値化することにより３次元でＣＮＶを含む血管領域を特定する。さらに網膜表層及び網膜深層で算出した血管密度と、網膜外層で特定及び計測した脈絡膜新生血管領域の２値画像と体積値を時系列表示する場合について説明する。
【0176】
本実施形態に係る画像処理装置１０１を備える画像処理システム１０の構成及び画像処理フローは第１実施形態の場合と同様であるので省略する。
【0177】
また、本実施形態での画像処理フローのうち図３のＳ３０６～Ｓ３０８、Ｓ３１０～Ｓ３１１以外は第１実施形態の場合と同様であるので説明は省略する。
【0178】
＜ステップ３０６＞
操作者が入力部１０３を用いてＯＣＴＡ計測処理の開始を指示する。
操作者が入力部１０３を用いてＯＣＴＡ計測処理の開始を指示する。
【0179】
本実施形態では図８（ｅ）のレポート画面８０３のモーションコントラスト画像上をダブルクリックすることで、図９（ａ）のようなＯＣＴＡ計測画面に移行する。モーションコントラスト画像が拡大表示され、解析部１０１－４６が計測処理を開始する。
【0180】
計測処理の種類として、本実施形態ではＭａｐボタン群もしくはＳｅｃｔｏｒボタン群から（ｉ）～（ｉｖ）のいずれかを選択する。
（ｉ） Ｎｏｎｅ（計測しない）
（ｉｉ） ＶＡＤ（血管が占める面積に基づいて算出される血管密度）
（ｉｉｉ） ＶＬＤ（血管の長さの総和に基づいて算出される血管密度）
（ｉｖ） Ｖｏｌｕｍｅ（血管領域の体積）
（ｉ） Ｎｏｎｅ（計測しない）
（ｉｉ） ＶＡＤ（血管が占める面積に基づいて算出される血管密度）
（ｉｉｉ） ＶＬＤ（血管の長さの総和に基づいて算出される血管密度）
（ｉｖ） Ｖｏｌｕｍｅ（血管領域の体積）
【0181】
これに限らず任意の種類の計測を実施してよい。
【0182】
例えば（ｉｖ）のＶｏｌｕｍｅの代わりに（２次元モーションコントラスト画像上で特定したり、特定した３次元血管領域を所定の深度範囲（例えば網膜外層）で投影したりして得られる）血管領域（例えば脈絡膜毛細血管）の面積を計測する場合も本発明に含まれる。
【0183】
また、Ｍｅａｓｕｒｅｍｅｎｔボタン経由で表示される選択画面からは（ｉ）～（ｉｖ）のいずれかを選択する。
（ｉ） 無血管領域の面積計測
（ｉｉ） 血管密度（ＶＡＤ）
（ｉｉｉ） 血管密度（ＶＬＤ）
（ｉｖ） 血管領域の体積（Ｖｏｌｕｍｅ）
（ｉ） 無血管領域の面積計測
（ｉｉ） 血管密度（ＶＡＤ）
（ｉｉｉ） 血管密度（ＶＬＤ）
（ｉｖ） 血管領域の体積（Ｖｏｌｕｍｅ）
【0184】
これに限らず、例えば（２次元モーションコントラスト画像上で特定したり、特定した３次元血管領域を所定の深度範囲（例えば網膜外層）で投影したりして得られる）血管領域（例えば脈絡膜毛細血管）の面積を計測してもよい。
【0185】
なお、３次元画像処理による計測は大きく（１）～（３）に大別できる。
（１）３次元で強調し、２次元投影した該強調画像上で特定した血管領域もしくは血管中心線データに対する２次元計測
（２）３次元で強調及び特定した血管領域もしくは血管中心線データを投影した場合の２次元計測
（３）３次元で強調及び特定された血管領域もしくは血管中心線データに対する３次元計測
（１）３次元で強調し、２次元投影した該強調画像上で特定した血管領域もしくは血管中心線データに対する２次元計測
（２）３次元で強調及び特定した血管領域もしくは血管中心線データを投影した場合の２次元計測
（３）３次元で強調及び特定された血管領域もしくは血管中心線データに対する３次元計測
【0186】
（１）及び（２）の例として、投影画像上で上記無血管領域の面積や血管密度、血管領域の面積や径、長さ、曲率を計測することが挙げられる。計測内容は正面モーションコントラスト画像に対する計測と同様であるものの、正面モーションコントラスト画像を強調・特定して計測する場合よりも血管抽出能が向上するため計測精度が向上する。
また、（３）の例として以下の例が挙げられる。
（３－１）血管の体積計測
（３－２）任意方向の断面画像もしくは曲断面画像上の計測
（血管の径もしくは断面積計測も含む）
（３－３）血管の長さや曲率計測
また、（３）の例として以下の例が挙げられる。
（３－１）血管の体積計測
（３－２）任意方向の断面画像もしくは曲断面画像上の計測
（血管の径もしくは断面積計測も含む）
（３－３）血管の長さや曲率計測
【0187】
本実施形態では、３次元で血管強調及び血管領域の特定処理を実行した後、網膜表層及び網膜深層の深度範囲で投影した各２値画像上でＶＡＤ、網膜外層の深度範囲で血管領域（脈絡膜新生血管領域）の体積を各々計測する。３次元モーションコントラスト画像に対してヘシアンフィルタを用いて線状構造を強調するため２次元モーションコントラスト画像で線として描出されていた本来的には線ではない構造を不要に強調することを避けることが可能となる。結果として精度の高い血管領域のセグメンテーション（特定）が可能となる。
【0188】
なお、第１実施形態の場合と同様に、Ｍａｐボタン群から選択した計測の種類、Ｓｅｃｔｏｒボタン群から選択した計測の種類のいずれか一方が変更された場合に、他方も連動して（同じ計測の種類に）変更されるように構成してもよい。
【0189】
次に、解析部１０１－４６は計測処理の前処理として画像拡大及びトップハットフィルタ処理を行う。本実施形態では、いずれも３次元のＢｉｃｕｂｉｃ補間及びトップハットフィルタ処理を実行する。
【0190】
＜ステップ３０７＞
解析部１０１－４６が血管領域の特定処理を行う。本実施形態では、強調部１０１－４６１が３次元ヘシアンフィルタ及び３次元エッジ選択鮮鋭化フィルタ処理に基づく血管強調処理を行う。次に抽出部１０１－４６２が第１実施形態の場合と同様に２種類の血管強調画像を用いて２値化処理を行い、整形処理を行うことで血管領域を特定する。
解析部１０１－４６が血管領域の特定処理を行う。本実施形態では、強調部１０１－４６１が３次元ヘシアンフィルタ及び３次元エッジ選択鮮鋭化フィルタ処理に基づく血管強調処理を行う。次に抽出部１０１－４６２が第１実施形態の場合と同様に２種類の血管強調画像を用いて２値化処理を行い、整形処理を行うことで血管領域を特定する。
【0191】
血管領域特定処理の詳細はＳ５１０～Ｓ５６０で説明する。
【0192】
＜ステップ３０８＞
計測部１０１－４６３が、操作者により指定された計測対象領域に関する情報に基づいて単検査の画像に対する計測を行う。引き続いて表示制御部１０１－０５が、計測結果を表示部１０４に表示する。
計測部１０１－４６３が、操作者により指定された計測対象領域に関する情報に基づいて単検査の画像に対する計測を行う。引き続いて表示制御部１０１－０５が、計測結果を表示部１０４に表示する。
【0193】
なお第１実施形態の場合と同様に、所定の解析に適した条件を満たさない状態で実施された計測結果を計測レポート画面に表示する際には表示部１０４に警告メッセージを表示するものとする。
【0194】
また、操作者が入力部１０３から血管領域もしくは血管中心線データを修正するよう指示を入力した場合には、第１実施形態の場合と同様に解析部１０１－４６が操作者から入力部１０３を介して指定された位置情報に基づいて血管領域もしくは血管中心線データを修正し、計測値を再計算する。
【0195】
網膜表層及び網膜深層におけるＶＡＤ計測、網膜外層における脈絡膜新生血管の体積計測についてはＳ８１０～Ｓ８３０、網膜表層及び網膜深層におけるＶＬＤ計測、網膜外層における脈絡膜新生血管の総血管長計測についてはＳ８４０～Ｓ８７０で各々説明する。
【0196】
＜ステップ３１０＞
比較部１０１－４６４は第１実施形態の場合と同様の操作で、経時変化計測（Ｐｒｏｇｒｅｓｓｉｏｎ計測）処理を実行する。
比較部１０１－４６４は第１実施形態の場合と同様の操作で、経時変化計測（Ｐｒｏｇｒｅｓｓｉｏｎ計測）処理を実行する。
【0197】
＜ステップ３１１＞
表示制御部１０１－０５は、Ｓ３１０で実施したＰｒｏｇｒｅｓｓｉｏｎ計測に関するレポートを表示部１０４に表示する。
表示制御部１０１－０５は、Ｓ３１０で実施したＰｒｏｇｒｅｓｓｉｏｎ計測に関するレポートを表示部１０４に表示する。
【0198】
本実施形態では、Ｐｒｏｇｒｅｓｓｉｏｎ計測レポートの最上段に網膜表層において計測したＶＡＤマップを表示し、２段目に網膜深層において計測したＶＡＤマップ、３段目に網膜外層において計測した脈絡膜新生血管領域の（ｉ）、（ｉｉ）を時系列で並置表示する。
（ｉ）２値画像（もしくは該２値画像の基準画像との差分画像）
（ｉｉ）体積値もしくは総血管長（もしくは該２値画像の基準画像との差分値）
（ｉ）２値画像（もしくは該２値画像の基準画像との差分画像）
（ｉｉ）体積値もしくは総血管長（もしくは該２値画像の基準画像との差分値）
【0199】
これに限らず、例えば４段目に脈絡膜における血管密度（ＶＡＤもしくはＶＬＤ）マップを時系列で並置表示してもよい。
【0200】
これにより、眼底血管の３次元的な病態に関する時系列変化を一覧して把握できる。
【0201】
また第１実施形態と同様に、各計測対象画像に関して略同一位置における断層画像数、ＯＣＴＡ重ね合わせ処理実施の有無、ＯＣＴＡ重ね合わせ処理の実施条件、ＯＣＴ断層画像もしくはモーションコントラスト画像の評価値の情報を表示部１０４に表示させてもよい。さらに、ショートカットメニューから計測の種類及び計測対象領域に関する項目を異なる値に変更することにより、Ｐｒｏｇｒｅｓｓｉｏｎ計測の種類及び計測対象領域を変更して再計測できる。第１実施形態と同様に、複数の計測対象領域が選択され、かつ一方の領域に対する計測の種類が変更された場合は他方の領域に対しても同じ計測の種類が連動して適用され、計測が実行される。さらに、第１実施形態と同様の方法で、所定の条件を満たさない状態で実施された計測結果を計測レポート画面に表示する際には警告メッセージを表示してよい。
【0202】
なお、本発明は異なる深度範囲の正面画像及び該正面画像に対する計測値分布の時系列表示に限られるものではなく、例えば該正面画像に直交する画像及び該正面画像に直交する画像に対する計測値分布、ボリュームレンダリングした３次元画像及び該３次元画像に対する計測値分布を時系列表示してもよい。
【0203】
さらに、図５（ａ）に示すフローチャートを参照しながら、Ｓ３０７で実行される処理の詳細について説明する。
【0204】
＜ステップ５１０＞
強調部１０１－４６１は、ステップ３０６の前処理を実施された３次元モーションコントラスト画像に対してヘッセ行列の固有値に基づく３次元血管強調フィルタ処理を行う。本実施形態では３次元のＭｕｌｔｉ－ｓｃａｌｅ ｌｉｎｅ ｆｉｌｔｅｒを用いるが、任意の公知の血管強調フィルタを用いてよい。
強調部１０１－４６１は、ステップ３０６の前処理を実施された３次元モーションコントラスト画像に対してヘッセ行列の固有値に基づく３次元血管強調フィルタ処理を行う。本実施形態では３次元のＭｕｌｔｉ－ｓｃａｌｅ ｌｉｎｅ ｆｉｌｔｅｒを用いるが、任意の公知の血管強調フィルタを用いてよい。
【0205】
なお、３次元のヘシアンフィルタでは３次元画像上の各画素において算出したヘッセ行列（（６）式）の「固有値（λ１、λ２、λ３）のうち１つが０に近く、残り２つが負かつ絶対値が大きい」場合に線状構造とみなして強調する。
【0206】
【数6】
【0207】
３次元ヘシアンフィルタを用いると、深度方向に屈曲した血管に関しても「血管走行方向の輝度変化が小さく、血管走行方向に直交する２方向の輝度が大きく低下する」という性質が成り立つため、良好に血管強調できるという利点がある。眼底血管の中には、深度方向に屈曲した血管として例えば以下の３点が挙げられる。
・脈絡膜側から網膜内に侵入する脈絡膜新生血管（ＣＮＶ）
・視神経乳頭部の血管
・網膜表層毛細血管と網膜深層毛細血管との接続部
・脈絡膜側から網膜内に侵入する脈絡膜新生血管（ＣＮＶ）
・視神経乳頭部の血管
・網膜表層毛細血管と網膜深層毛細血管との接続部
【0208】
上記血管に対して正面モーションコントラスト画像上で２次元のヘシアンフィルタを適用すると、２次元平面内では「該平面内での血管走行方向の輝度変化が小さく、血管走行方向に直交する方向の輝度が大きく低下する」という性質が成り立たないため十分強調されず、血管領域として特定できないという課題がある。３次元ヘシアンフィルタを用いると上記血管についても良好に強調でき、血管検出能が向上する。
【0209】
＜ステップ５２０＞
抽出部１０１－４６２は、Ｓ５１０で生成した３次元ヘシアンフィルタによる血管強調画像（以下、３次元ヘシアン強調画像と表記）を２値化する。
抽出部１０１－４６２は、Ｓ５１０で生成した３次元ヘシアンフィルタによる血管強調画像（以下、３次元ヘシアン強調画像と表記）を２値化する。
【0210】
２値化の手順は第１実施形態の場合と同様であるが、３次元データの２値化であるという点が第１実施形態の場合と異なる。また、合成モーションコントラスト画像をヘシアンフィルタで強調処理した画像であるため、単独のモーションコントラスト画像をヘシアンフィルタで強調した場合に比べて２値化した血管領域の連続性が向上する。
【0211】
＜ステップ５３０＞
強調部１０１－４６１は、Ｓ３０６で生成したトップハットフィルタ適用後の合成モーションコントラスト画像に対して３次元エッジ選択鮮鋭化処理を行う。本実施形態では、３次元モーションコントラスト画像に対して３次元のＳｏｂｅｌフィルタを適用した画像を重みとして３次元のアンシャープマスク処理を行うことにより、エッジ選択鮮鋭化処理を実施する。
強調部１０１－４６１は、Ｓ３０６で生成したトップハットフィルタ適用後の合成モーションコントラスト画像に対して３次元エッジ選択鮮鋭化処理を行う。本実施形態では、３次元モーションコントラスト画像に対して３次元のＳｏｂｅｌフィルタを適用した画像を重みとして３次元のアンシャープマスク処理を行うことにより、エッジ選択鮮鋭化処理を実施する。
【0212】
＜ステップ５４０＞
抽出部１０１－４６２は、Ｓ５３０で生成したエッジ選択鮮鋭化処理を適用した鮮鋭化画像を２値化する。任意の公知の２値化法を用いてよいが、本実施形態では該３次元鮮鋭化画像上の各３次元局所領域内で算出した輝度統計値（平均値もしくは中央値）を閾値として２値化する。第１実施形態の場合と同様に、該閾値の上限及び下限値を設定することにより、血管領域内における抽出不足や無血管領域における誤抽出を抑制する。
抽出部１０１－４６２は、Ｓ５３０で生成したエッジ選択鮮鋭化処理を適用した鮮鋭化画像を２値化する。任意の公知の２値化法を用いてよいが、本実施形態では該３次元鮮鋭化画像上の各３次元局所領域内で算出した輝度統計値（平均値もしくは中央値）を閾値として２値化する。第１実施形態の場合と同様に、該閾値の上限及び下限値を設定することにより、血管領域内における抽出不足や無血管領域における誤抽出を抑制する。
【0213】
またＳ５２０の場合と同様、合成モーションコントラスト画像をエッジ選択鮮鋭化した画像であるため、単独のモーションコントラスト画像をエッジ選択鮮鋭化した場合に比べて２値化した場合のノイズ状の誤検出領域が減少する。
【0214】
＜ステップ５５０＞
抽出部１０１－４６２は、Ｓ５２０で生成した３次元ヘシアン強調画像の２値画像の輝度値と、Ｓ５４０で生成した３次元エッジ選択強調画像の２値画像の輝度値の双方が０より大きい場合に血管候補領域として抽出する。当該演算処理により、ヘシアン強調画像に見られる血管径を過大評価している領域と、エッジ選択強調画像上に見られるノイズ領域がともに抑制され、血管の境界位置が正確かつ血管の連続性が良好な２値画像を取得できる。
抽出部１０１－４６２は、Ｓ５２０で生成した３次元ヘシアン強調画像の２値画像の輝度値と、Ｓ５４０で生成した３次元エッジ選択強調画像の２値画像の輝度値の双方が０より大きい場合に血管候補領域として抽出する。当該演算処理により、ヘシアン強調画像に見られる血管径を過大評価している領域と、エッジ選択強調画像上に見られるノイズ領域がともに抑制され、血管の境界位置が正確かつ血管の連続性が良好な２値画像を取得できる。
【0215】
また双方の２値画像とも合成モーションコントラスト画像に基づく２値画像であることから、単独のモーションコントラスト画像に基づく２値画像に比べて２値化した場合のノイズ状の誤検出領域が減少するとともに、特に毛細血管領域の連続性が向上する。また、合成モーションコントラスト画像であることから検査間の画質や輝度レベルも安定し、血管の抽出能が検査間で安定しやすい。
【0216】
＜ステップ５６０＞
抽出部１０１－４６２は、血管領域の整形処理として３次元のオープニング処理（収縮処理後に膨張処理を行うこと）及びクロージング処理（膨張処理後に収縮処理を行うこと）を実施する。なお、整形処理はこれに限らず２値画像をラベリングした場合の各ラベルの面積に基づく小領域除去を行ってもよい。
抽出部１０１－４６２は、血管領域の整形処理として３次元のオープニング処理（収縮処理後に膨張処理を行うこと）及びクロージング処理（膨張処理後に収縮処理を行うこと）を実施する。なお、整形処理はこれに限らず２値画像をラベリングした場合の各ラベルの面積に基づく小領域除去を行ってもよい。
【0217】
なお第１実施形態の場合と同様に、様々な径の血管が含まれるモーションコントラスト画像において血管強調する際のスケールを適応的に決定する方法はＳ５１０～Ｓ５６０に述べた方法に限定されない。例えば図５（ｂ）のＳ６１０～Ｓ６５０に示すように、３次元ヘシアン強調画像の輝度値と３次元エッジ選択鮮鋭化による血管強調画像の輝度値を乗じる演算を適用した画像に対する輝度統計値（例えば平均値）を閾値として２値化することによって血管領域を特定してもよい。該閾値には下限値や上限値を設定できる。
【0218】
あるいは、図５（ｃ）のＳ７１０～Ｓ７４０に示すように、各画素の３次元位置（固視位置や深度範囲のデータでもよい）に基づいてヘシアンフィルタ適用時の平滑化フィルタのパラメータ（ガウスフィルタの平滑化パラメータσ）を適応的に変えた上でヘシアンフィルタを適用し、２値化することによって血管強調してもよい。
【0219】
また、２値化処理は閾値処理に限定されるものではなく任意の公知のセグメンテーション手法を用いることができる。
【0220】
さらに、図６（ａ）に示すフローチャートを参照しながら、Ｓ３０８で実行される処理の詳細について説明する。
【0221】
＜ステップ８１０＞
操作者は、入力部１０３を介して計測処理における関心領域を設定する。
操作者は、入力部１０３を介して計測処理における関心領域を設定する。
【0222】
本実施形態では計測内容として
１）網膜表層及び網膜深層におけるＶＡＤマップ及びＶＡＤセクタマップ
２）網膜外層における脈絡膜新生血管の体積
を算出する。従って、関心領域として網膜表層及び網膜深層においては画像全体及び固視灯位置を中心とするセクタ領域を選択する。また網膜外層においては網膜外層に相当する層境界（ＯＰＬ／ＯＮＬ境界と、ブルッフ膜境界を該境界深層側に２０μｍ移動させた位置で囲まれる範囲）を指定する。
１）網膜表層及び網膜深層におけるＶＡＤマップ及びＶＡＤセクタマップ
２）網膜外層における脈絡膜新生血管の体積
を算出する。従って、関心領域として網膜表層及び網膜深層においては画像全体及び固視灯位置を中心とするセクタ領域を選択する。また網膜外層においては網膜外層に相当する層境界（ＯＰＬ／ＯＮＬ境界と、ブルッフ膜境界を該境界深層側に２０μｍ移動させた位置で囲まれる範囲）を指定する。
【0223】
＜ステップ８２０＞
計測部１０１－４６３は、Ｓ３０７で得られた血管領域の２値画像に基づいて計測処理を行う。網膜表層及び網膜深層における計測内容（ＶＡＤマップ及びＶＡＤセクタマップ生成）は基本的に第１実施形態と同様であるが、網膜表層及び網膜深層において特定された３次元血管領域を正面像として投影してから計測する点が異なる。網膜外層では、Ｓ８１０で設定した網膜外層に相当する関心領域内における非０画素（白画素）の体積を算出する。
計測部１０１－４６３は、Ｓ３０７で得られた血管領域の２値画像に基づいて計測処理を行う。網膜表層及び網膜深層における計測内容（ＶＡＤマップ及びＶＡＤセクタマップ生成）は基本的に第１実施形態と同様であるが、網膜表層及び網膜深層において特定された３次元血管領域を正面像として投影してから計測する点が異なる。網膜外層では、Ｓ８１０で設定した網膜外層に相当する関心領域内における非０画素（白画素）の体積を算出する。
【0224】
＜ステップ８３０＞
表示制御部１０１－０５は、計測結果としてＳ８２０で生成した網膜表層及び網膜深層におけるＶＡＤマップ及びＶＡＤセクタマップ、網膜外層における血管領域の２値画像と該血管領域の体積値を表示部１０４に表示する。
表示制御部１０１－０５は、計測結果としてＳ８２０で生成した網膜表層及び網膜深層におけるＶＡＤマップ及びＶＡＤセクタマップ、網膜外層における血管領域の２値画像と該血管領域の体積値を表示部１０４に表示する。
【0225】
なお、第１実施形態と同様の方法で、推奨する計測条件の少なくとも一つを満たさない状態で実施された計測結果を計測レポート画面に表示する際には（正確に計測することが期待できないような条件で計測が行われたとみなして）警告表示を行ってよい。
【0226】
上記では特定された３次元血管領域に基づいて体積を計測する場合の手順を例に説明したが、３次元の血管中心線に基づいて計測する場合には、上記Ｓ８１０～８３０の代わりに図６（ｂ）に示すＳ８４０～８７０を実行する。
【0227】
＜ステップ８４０＞
計測部１０１－４６３は、Ｓ８２０で生成した血管領域の２値画像を３次元細線化処理することにより、血管の中心線に相当する線幅１画素のスケルトン画像を生成する。
計測部１０１－４６３は、Ｓ８２０で生成した血管領域の２値画像を３次元細線化処理することにより、血管の中心線に相当する線幅１画素のスケルトン画像を生成する。
【0228】
＜ステップ８５０＞
操作者は、入力部１０３を介してＳ８１０と場合と同様の関心領域を設定する。なお、本実施形態では計測内容として（１）、（２）を算出する。
（１）網膜表層及び網膜深層におけるＶＬＤマップ及びＶＬＤセクタマップ
（２）網膜外層における脈絡膜新生血管の総血管長
操作者は、入力部１０３を介してＳ８１０と場合と同様の関心領域を設定する。なお、本実施形態では計測内容として（１）、（２）を算出する。
（１）網膜表層及び網膜深層におけるＶＬＤマップ及びＶＬＤセクタマップ
（２）網膜外層における脈絡膜新生血管の総血管長
【0229】
＜ステップ８６０＞
計測部１０１－４６３はＳ８４０で得られたスケルトン画像に基づいて計測処理を行う。網膜表層及び網膜深層における計測内容（ＶＬＤマップ及びＶＬＤセクタマップ生成）は基本的に第１実施形態と同様であるが、網膜表層及び網膜深層において特定された３次元スケルトンを正面像として投影してから計測する点が異なる。なお、Ｓ３０７で特定した３次元血管領域を正面像として投影してから２次元細線化処理を実行して計測してもよい。網膜外層では、Ｓ８１０で設定した網膜外層に相当する関心領域内における非０画素（白画素）の長さの総和を算出する。
計測部１０１－４６３はＳ８４０で得られたスケルトン画像に基づいて計測処理を行う。網膜表層及び網膜深層における計測内容（ＶＬＤマップ及びＶＬＤセクタマップ生成）は基本的に第１実施形態と同様であるが、網膜表層及び網膜深層において特定された３次元スケルトンを正面像として投影してから計測する点が異なる。なお、Ｓ３０７で特定した３次元血管領域を正面像として投影してから２次元細線化処理を実行して計測してもよい。網膜外層では、Ｓ８１０で設定した網膜外層に相当する関心領域内における非０画素（白画素）の長さの総和を算出する。
【0230】
＜ステップ８７０＞
表示制御部１０１－０５は、計測結果としてＳ８６０で生成した網膜表層及び網膜深層におけるＶＬＤマップ及びＶＬＤセクタマップ、網膜外層におけるスケルトン画像と該スケルトンの長さの総和を表示部１０４に表示する。
表示制御部１０１－０５は、計測結果としてＳ８６０で生成した網膜表層及び網膜深層におけるＶＬＤマップ及びＶＬＤセクタマップ、網膜外層におけるスケルトン画像と該スケルトンの長さの総和を表示部１０４に表示する。
【0231】
またＳ８３０の場合と同様に、所定の解析に適した条件を満たさない状態で実施された計測結果を計測レポート画面に表示する際には警告メッセージを表示するものとする。
【0232】
なお、本実施形態では脈絡膜新生血管を３次元抽出して体積や血管長の総和を時系列表示する場合の手順について説明したが、本発明はこれに限定されない。例えば、図１１の「Ｓｈｏｗ Ｄｉｆｆｅｒｅｎｃｅ」チェックボックスを押下して基準画像における該新生血管の体積とその他の画像における該新生血管の体積との差分画像及び差分値を生成して時系列で表示してもよい。
【0233】
あるいは、視神経乳頭部の深度方向に屈曲した動静脈領域を本実施形態で説明した手順と同様にして特定し、血管径や血管断面積、血管中心線の曲率といった血管形状を計測してもよい。あるいは、網膜表層の毛細血管と網膜深層の毛細血管の連結部を３次元的に抽出して強調表示させたり、連結部の本数を計数したりしてもよい。
【0234】
以上述べた構成によれば、画像処理装置１０１は３次元の合成モーションコントラスト画像に対してモーションアーチファクト抑制処理を実施後、３次元のモルフォロジーフィルタ及び血管強調フィルタを適用し、２値化することにより３次元で血管領域を特定する。さらに特定した血管領域の体積を算出して血管領域の２値画像と体積値を時系列表示する。
【0235】
これにより、検査ごとのＯＣＴ断層画像の信号強度や画質のばらつきの影響を抑制しながら、血管病態の変化を正確に特定・計測できる。
【0236】
［第３の実施形態］
上記の実施形態では複数のクラスタそれぞれで得られたＯＣＴＡ画像を重ね合わせて（加算平均）する場合を主として述べたが、本発明はこれに限定されるものではない。例えば、重ね合わせ後のＯＣＴＡ画像に替えて１つのクラスタにおける断層画像の枚数を９以上として得られたＯＣＴＡ画像を用いてもよい。このようにすればクラスタ数を４、各クラスタの断層画像の枚数を３とした場合と同等のＯＣＴＡ画像を得ることが可能となる。
上記の実施形態では複数のクラスタそれぞれで得られたＯＣＴＡ画像を重ね合わせて（加算平均）する場合を主として述べたが、本発明はこれに限定されるものではない。例えば、重ね合わせ後のＯＣＴＡ画像に替えて１つのクラスタにおける断層画像の枚数を９以上として得られたＯＣＴＡ画像を用いてもよい。このようにすればクラスタ数を４、各クラスタの断層画像の枚数を３とした場合と同等のＯＣＴＡ画像を得ることが可能となる。
【0237】
［第４の実施形態］
上記の実施形態では図１１の時系列表示において２次元のＯＣＴＡ画像（モーションコントラスト画像）が表示されていたが、３次元のＯＣＴＡ画像を時系列に表示することとしてもよい。
上記の実施形態では図１１の時系列表示において２次元のＯＣＴＡ画像（モーションコントラスト画像）が表示されていたが、３次元のＯＣＴＡ画像を時系列に表示することとしてもよい。
【0238】
［その他の実施形態］
以上、実施形態例を詳述したが、開示の技術は例えば、システム、装置、方法、プログラム若しくは記録媒体（記憶媒体）等としての実施態様をとることが可能である。具体的には、複数の機器（例えば、ホストコンピュータ、インタフェース機器、撮像装置、ｗｅｂアプリケーション等）から構成されるシステムに適用しても良いし、また、一つの機器からなる装置に適用しても良い。
以上、実施形態例を詳述したが、開示の技術は例えば、システム、装置、方法、プログラム若しくは記録媒体（記憶媒体）等としての実施態様をとることが可能である。具体的には、複数の機器（例えば、ホストコンピュータ、インタフェース機器、撮像装置、ｗｅｂアプリケーション等）から構成されるシステムに適用しても良いし、また、一つの機器からなる装置に適用しても良い。
【0239】
また、本発明の目的は、以下のようにすることによって達成されることはいうまでもない。即ち、前述した実施形態の機能を実現するソフトウェアのプログラムコード（コンピュータプログラム）を記録した記録媒体（または記憶媒体）を、システムあるいは装置に供給する。係る記憶媒体は言うまでもなく、コンピュータ読み取り可能な記憶媒体である。そして、そのシステムあるいは装置のコンピュータ（またはＣＰＵやＭＰＵ）が記録媒体に格納されたプログラムコードを読み出し実行する。この場合、記録媒体から読み出されたプログラムコード自体が前述した実施形態の機能を実現することになり、そのプログラムコードを記録した記録媒体は本発明を構成することになる。
【0240】
なお、上述した実施例および変形例は適宜組み合わせて実施してもよい。
【符号の説明】
【0241】
１０１－１２ モーションコントラストデータ生成部
１０１－４２ 合成部
１０１－０５ 表示制御部
１０１－４２ 合成部
１０１－０５ 表示制御部
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
