(57)【要約】【課題】  モーションコントラスト画像上の複数の解析領域に対して実行された解析の結果を示す画像が表示される場合に、複数の解析領域に対する解析の種類として互いに対応する種類の選択が容易に実行されるように構成すること。【解決手段】  画像処理装置は、眼部のモーションコントラスト画像における第１の領域と第１の領域よりも狭い領域を少なくとも含む第２の領域とのいずれか一方に対して選択された解析の種類を示す情報を用いて一方が解析された結果を示す画像を表示手段に表示させ、一方に対する選択の後に操作者からの指示に応じて他方に対して選択された解析の種類と、一方に対して選択された解析の種類とが異なる場合には、一方が解析された結果を示す画像の表示を、他方に対して選択された解析の種類を示す情報を用いて一方及び他方が解析された結果を示す画像の表示に変更する制御を実行する表示制御手段を有する。【選択図】  図９
2019154716.tif 000002
【特許請求の範囲】
【請求項１】
眼部のモーションコントラスト画像における第１の領域と前記第１の領域よりも狭い領域を少なくとも含む第２の領域とのうち少なくとも１つに対する解析を実行する解析手段と、
操作者からの指示に応じて前記第１の領域と前記第２の領域とのいずれか一方に対して選択された解析の種類を示す情報を用いて前記一方が解析された結果を示す画像を表示手段に表示させ、前記一方に対する選択の後に操作者からの指示に応じて他方に対して選択された解析の種類と、前記一方に対して選択された解析の種類とが異なる場合には、前記一方が解析された結果を示す画像の表示領域において、前記一方が解析された結果を示す画像の表示を、前記他方に対して選択された解析の種類を示す情報を用いて前記一方及び前記他方が解析された結果を示す画像の表示に変更する制御を実行する表示制御手段と、
を有することを特徴とする画像処理装置。
【請求項２】
前記表示制御手段は、前記一方に対する選択の後に操作者からの指示に応じて他方に対して選択された解析の種類と、前記一方に対して選択された解析の種類とが異なる場合には、前記一方に対して選択された解析の種類を示す情報の表示を、前記他方に対して選択された解析の種類を示す情報の表示に変更する別の制御を実行することを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記表示制御手段は、複数の検査日に対応する前記眼部の複数のモーションコントラスト画像が時系列に並べて表示される表示領域において、前記複数のモーションコントラスト画像のうちの１つの画像に対して実行される前記制御を、他の画像に対しても実行することを特徴とする請求項１または２に記載の画像処理装置。
【請求項４】
前記解析手段は、操作者からの指示に応じて、複数の検査日に対応する前記眼部の複数のモーションコントラスト画像のうちの１つの画像における前記第１の領域と前記第２の領域との少なくとも１つに対して選択される解析の種類を示す情報を、他の画像に対しても適用することを特徴とする請求項１乃至３のいずれか１項に記載の画像処理装置。
【請求項５】
前記第１の領域と前記第２の領域とのいずれか一方に対して選択された解析の種類が操作者からの指示に応じて非選択に変更された場合には、他方に対して選択された解析の種類は変更されないように構成されることを特徴とする請求項１乃至４のいずれか１項に記載の画像処理装置。
【請求項６】
前記表示制御手段は、前記第１の領域及び前記第２の領域に対する解析の種類が非選択の場合には、解析された結果を示す画像が前記表示領域において非表示の状態で、前記モーションコントラスト画像を前記表示領域に表示させることを特徴とする請求項１乃至５のいずれか１項に記載の画像処理装置。
【請求項７】
操作者からの指示に応じて選択された解析の種類を示す情報を用いて、眼部のモーションコントラスト画像における第１の領域に対する解析を実行する解析手段と、
操作者からの指示に応じて前記モーションコントラスト画像における前記第１の領域よりも狭い領域を少なくとも含む第２の領域に対して解析された結果を示す画像の表示が選択された場合には、前記選択された解析の種類を示す情報を用いて前記第２の領域に対して解析された結果を示す画像を、前記第１の領域に対して解析された結果を示す画像に重畳された状態で表示手段に表示させる表示制御手段と、
を有することを特徴とする画像処理装置。
【請求項８】
前記選択された解析の種類は、前記モーションコントラスト画像において特定される血管領域に関するパラメータを少なくとも含む複数の種類のいずれかであることを特徴とする請求項１乃至７のいずれか１項に記載の画像処理装置。
【請求項９】
前記選択された解析の種類は、前記モーションコントラスト画像において特定される血管領域の面積に関する血管密度と、前記血管領域の長さに関する血管密度とを少なくとも含む複数の種類のいずれかであることを特徴とする請求項１乃至８のいずれか１項に記載の画像処理装置。
【請求項１０】
前記選択された解析の種類は、前記モーションコントラスト画像において特定される無血管領域に関するパラメータを少なくとも含む複数の種類のいずれかであることを特徴とする請求項１乃至９のいずれか１項に記載の画像処理装置。
【請求項１１】
前記第２の領域は、セクタ領域で、
前記第１の領域は、前記セクタ領域よりも広い領域であることを特徴とする請求項１乃至１０のいずれか１項に記載の画像処理装置。
【請求項１２】
前記モーションコントラスト画像は、前記眼部の同一位置を測定光が走査されるように制御して得た複数の３次元のモーションコントラスト画像が合成された画像であることを特徴とする請求項１乃至１１のいずれか１項に記載の画像処理装置。
【請求項１３】
前記モーションコントラスト画像は、前記眼部の３次元のモーションコントラスト画像と、操作者からの指示に応じて設定された深さ範囲に関する情報とを用いて生成された、前記眼部のモーションコントラスト正面画像であり、
前記解析された結果を示す画像は、前記モーションコントラスト正面画像の少なくとも一部の領域に対して解析された結果を示す２次元画像であることを特徴とする請求項１乃至１２のいずれか１項に記載の画像処理装置。
【請求項１４】
眼部のモーションコントラスト画像における第１の領域と前記第１の領域よりも狭い領域を少なくとも含む第２の領域とのうち少なくとも１つに対する解析を実行する工程と、
操作者からの指示に応じて前記第１の領域と前記第２の領域とのいずれか一方に対して選択された解析の種類を示す情報を用いて前記一方が解析された結果を示す画像を表示手段に表示させ、前記一方に対する選択の後に操作者からの指示に応じて他方に対して選択された解析の種類と、前記一方に対して選択された解析の種類とが異なる場合には、前記一方が解析された結果を示す画像の表示領域において、前記一方が解析された結果を示す画像の表示を、前記他方に対して選択された解析の種類を示す情報を用いて前記一方及び前記他方が解析された結果を示す画像の表示に変更する制御を実行する工程と、
を有することを特徴とする画像処理方法。
【請求項１５】
操作者からの指示に応じて選択された解析の種類を示す情報を用いて、眼部のモーションコントラスト画像における第１の領域に対する解析を実行する工程と、
操作者からの指示に応じて前記モーションコントラスト画像における前記第１の領域よりも狭い領域を少なくとも含む第２の領域に対して解析された結果を示す画像の表示が選択された場合には、前記選択された解析の種類を示す情報を用いて前記第２の領域に対して解析された結果を示す画像を、前記第１の領域に対して解析された結果を示す画像に重畳された状態で表示手段に表示させる工程と、
を有することを特徴とする画像処理方法。
【請求項１６】
請求項１４または１５に記載の画像処理方法の各工程をコンピュータに実行させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
開示の技術は、画像処理装置、画像処理方法及びプログラムに関する。
【背景技術】
【０００２】
光干渉断層計（ＯＣＴ；Ｏｐｔｉｃａｌ  Ｃｏｈｅｒｅｎｃｅ  Ｔｏｍｏｇｒａｐｈｙ）などの眼部の断層画像撮影装置を用いると、網膜層内部の状態を３次元的に観察できる。この断層画像撮影装置は、疾病の診断をより的確に行うのに有用であることから眼科診療に広く用いられている。ＯＣＴの形態として、例えば広帯域な光源とマイケルソン干渉計を組み合わせたＴＤ−ＯＣＴ（Ｔｉｍｅ  ｄｏｍａｉｎ  ＯＣＴ）がある。これは、参照ミラーの位置を一定速度で移動させて信号アームで取得した後方散乱光との干渉光を計測し、深さ方向の反射光強度分布を得るように構成されている。しかし、このようなＴＤ−ＯＣＴでは機械的な走査が必要となるため高速な画像取得は難しい。そこで、より高速な画像取得法として広帯域光源を用い、分光器で干渉信号を取得するＳＤ−ＯＣＴ（Ｓｐｅｃｔｒａｌ  ｄｏｍａｉｎ  ＯＣＴ）や高速波長掃引光源を用いることで時間的に分光するＳＳ−ＯＣＴ（Ｓｗｅｐｔ  Ｓｏｕｒｃｅ  ＯＣＴ）が開発され、より広画角な断層画像を取得できるようになっている。
【０００３】
一方、眼科診療では眼底血管の病態を把握するためにこれまで侵襲的な蛍光眼底造影検査が行われてきた。近年は、ＯＣＴを用いて非侵襲に眼底血管を３次元で描出するＯＣＴ  Ａｎｇｉｏｇｒａｐｈｙ（以下、ＯＣＴＡと表記）技術が用いられるようになってきている。ＯＣＴＡでは測定光で同一位置を複数回走査し、赤血球の変位と測定光との相互作用により得られるモーションコントラストを画像化する。図４は主走査方向が水平（ｘ軸）方向で、副走査方向（ｙ軸方向）の各位置（ｙｉ；１≦ｉ≦ｎ）においてｒ回連続でＢスキャンを行うＯＣＴＡ撮影の例を示している。なおＯＣＴＡ撮像において同一位置で複数回走査することをクラスタ走査、同一位置で得られた複数枚の断層画像のことをクラスタと呼ぶ。クラスタ単位でモーションコントラストデータを生成し、１クラスタあたりの断層画像数（略同一位置での走査回数）を増やすと、ＯＣＴＡ画像のコントラストが向上することが知られている。
【０００４】
ここで、モーションコントラストデータを用いて算出された血管解析マップの表示に関する技術が、特許文献に開示されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１７−７７４１４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
ところで、ＯＣＴＡ画像を解析する場合には、解析の種類（例えば、血管領域の面積に関する血管密度（Ｖｅｓｓｅｌ  Ａｒｅａ  Ｄｅｎｓｉｔｙ；ＶＡＤ）、血管長に関する血管密度（Ｖｅｓｓｅｌ  Ｌｅｎｇｔｈ  Ｄｅｎｓｉｔｙ；ＶＬＤ）、無血管領域の面積等）や、解析領域等の条件を設定した上で解析が行われることが望まれる。このとき、解析領域としては、画像全体、ＥＴＤＲＳチャートのような扇形（セクタ）領域、操作者が設定した任意形状の関心領域などが選択されることが望まれる。ここで、例えば、１枚のＯＣＴＡ画像上の各解析領域で異なる種類の解析を選択可能にした場合を考える。このとき、例えば、１枚のＯＣＴＡ画像上に、ＶＡＤマップとＶＬＤセクタマップのような異なる種類の解析による結果を示す画像として解析マップが重畳表示される場合が生じる。このような場合には、操作者がデータの解釈を誤ってしまう可能性がある。
【０００７】
開示の技術は、モーションコントラスト画像上の複数の解析領域に対して実行された解析の結果を示す画像が表示される場合に、複数の解析領域に対する解析の種類として互いに対応する種類の選択が容易に実行されるように構成することを目的の一つとする。
【０００８】
なお、上記目的に限らず、後述する発明を実施するための形態に示す各構成により導かれる作用効果であって、従来の技術によっては得られない作用効果を奏することも本件の他の目的の１つとして位置付けることができる。
【課題を解決するための手段】
【０００９】
上記目的を達成するための、開示の画像処理装置の一つは、
眼部のモーションコントラスト画像における第１の領域と前記第１の領域よりも狭い領域を少なくとも含む第２の領域とのうち少なくとも１つに対する解析を実行する解析手段と、
操作者からの指示に応じて前記第１の領域と前記第２の領域とのいずれか一方に対して選択された解析の種類を示す情報を用いて前記一方が解析された結果を示す画像を表示手段に表示させ、前記一方に対する選択の後に操作者からの指示に応じて他方に対して選択された解析の種類と、前記一方に対して選択された解析の種類とが異なる場合には、前記一方が解析された結果を示す画像の表示領域において、前記一方が解析された結果を示す画像の表示を、前記他方に対して選択された解析の種類を示す情報を用いて前記一方及び前記他方が解析された結果を示す画像の表示に変更する制御を実行する表示制御手段と、を有する。
【発明の効果】
【００１０】
開示の技術の一つによれば、モーションコントラスト画像上の複数の解析領域に対して実行された解析の結果を示す画像が表示される場合に、複数の解析領域に対する解析の種類として互いに対応する種類の選択が容易に実行されるように構成することができる。
【図面の簡単な説明】
【００１１】
【図１】本発明の第一実施形態に係る画像処理装置の構成を示すブロック図である。
【図２】本発明の実施形態に係る画像処理システムや、該画像処理システムを構成する断層画像撮影装置に含まれる測定光学系を説明する図である。
【図３】本発明の第一実施形態に係る画像処理システムが実行可能な処理のフローチャートである。
【図４】本発明の実施形態におけるＯＣＴＡ撮影の走査方法を説明する図である。
【図５】本発明の第一実施形態のＳ３０７で実行される処理を説明する図である。
【図６】本発明の第一実施形態のＳ３０８で実行される処理を説明する図である。
【図７】本発明の第一実施形態において表示手段に表示する基準検査の選択画面及び撮影画面を説明する図である。
【図８】本発明の第一実施形態のＳ３０４における画像処理内容及びＳ３０５において表示手段に表示するレポート画面を説明する図である。
【図９】本発明の第一実施形態において表示手段に表示する計測操作画面とＳ３０８において表示する計測レポート画面を説明する図である。
【図１０】本発明の第一実施形態において特定された血管領域をユーザが修正する場合の操作手順と実行される画像処理内容を説明する図である。
【図１１】本発明の第一実施形態のＳ３１１において表示手段に表示する経時変化計測レポート画面を説明する図である。
【図１２】本発明の第一実施形態の警告メッセージが表示された計測レポート画面を説明する図である。
【発明を実施するための形態】
【００１２】
［第一実施形態］
本実施形態に係る画像処理装置の一つは、同一被検眼に対して異なる日時に略同一撮影条件（経過観察を目的とするフォローアップ撮影の条件）で取得したＯＣＴＡ重ね合わせ画像から生成した網膜表層及び網膜深層の正面モーションコントラスト画像を用いて血管領域特定及び血管密度計測処理を実施する。該特定処理及び計測処理により得られた合成画像及び計測値を複数の深度範囲で時系列並置表示する場合について説明する。なお、本発明において、経過観察を目的とするフォローアップ撮影や、画像の重ね合わせ処理等は必須のものではない。
【００１３】
以下、図面を参照しながら、本発明の第一実施形態に係る画像処理装置を備える画像処理システムについて説明する。
【００１４】
図２は、本実施形態に係る画像処理装置１０１を備える画像処理システム１０の構成を示す図である。図２に示すように、画像処理システム１０は、画像処理装置１０１が、インタフェースを介して断層画像撮影装置１００（ＯＣＴとも言う）、外部記憶部１０２、入力部１０３、表示部１０４と接続されることにより構成されている。
【００１５】
断層画像撮影装置１００は、眼部の断層画像を撮影する装置である。本実施形態においては、断層画像撮影装置１００としてＳＤ−ＯＣＴを用いるものとする。これに限らず、例えばＳＳ−ＯＣＴを用いて構成しても良い。
【００１６】
図２（ａ）において、測定光学系１００−１は前眼部像、被検眼のＳＬＯ眼底像、断層画像を取得するための光学系である。ステージ部１００−２は、測定光学系１００−１を前後左右に移動可能にする。ベース部１００−３は、後述の分光器を内蔵している。
【００１７】
画像処理装置１０１は、ステージ部１００−２の制御、アラインメント動作の制御、断層画像の再構成などを実行するコンピュータである。外部記憶部１０２は、断層撮像用のプログラム、患者情報、撮影データ、過去検査の画像データや計測データなどを記憶する。
【００１８】
入力部１０３はコンピュータへの指示を行い、具体的にはキーボードとマウスから構成される。表示部１０４は、例えばモニタからなる。
【００１９】
（断層画像撮影装置の構成）
本実施形態の断層画像撮影装置１００における測定光学系及び分光器の構成について図２（ｂ）を用いて説明する。
【００２０】
まず、測定光学系１００−１の内部について説明する。被検眼２００に対向して対物レンズ２０１が設置され、その光軸上に第１ダイクロイックミラー２０２及び第２ダイクロイックミラー２０３が配置されている。これらのダイクロイックミラーによってＯＣＴ光学系の光路２５０、ＳＬＯ光学系と固視灯用の光路２５１、及び前眼観察用の光路２５２とに波長帯域ごとに分岐される。
【００２１】
ＳＬＯ光学系と固視灯用の光路２５１は、ＳＬＯ走査手段２０４、レンズ２０５及び２０６、ミラー２０７、第３ダイクロイックミラー２０８、ＡＰＤ（Ａｖａｌａｎｃｈｅ  Ｐｈｏｔｏｄｉｏｄｅ）２０９、ＳＬＯ光源２１０、固視灯２１１を有している。
【００２２】
ミラー２０７は、穴あきミラーや中空のミラーが蒸着されたプリズムであり、ＳＬＯ光源２１０による照明光と、被検眼からの戻り光とを分離する。第３ダイクロイックミラー２０８はＳＬＯ光源２１０の光路と固視灯２１１の光路とに波長帯域ごとに分離する。
【００２３】
ＳＬＯ走査手段２０４は、ＳＬＯ光源２１０から発せられた光を被検眼２００上で走査するものであり、Ｘ方向に走査するＸスキャナ、Ｙ方向に走査するＹスキャナから構成されている。本実施形態では、Ｘスキャナは高速走査を行う必要があるためポリゴンミラーで、Ｙスキャナはガルバノミラーによって構成されている。
【００２４】
レンズ２０５はＳＬＯ光学系及び固視灯２１１の焦点合わせのため、不図示のモータによって駆動される。ＳＬＯ光源２１０は７８０ｎｍ付近の波長の光を発生する。ＡＰＤ２０９は、被検眼からの戻り光を検出する。固視灯２１１は可視光を発生して被検者の固視を促すものである。
【００２５】
ＳＬＯ光源２１０から発せられた光は、第３ダイクロイックミラー２０８で反射され、ミラー２０７を通過し、レンズ２０６及び２０５を通ってＳＬＯ走査手段２０４によって被検眼２００上で走査される。被検眼２００からの戻り光は、照明光と同じ経路を戻った後、ミラー２０７によって反射され、ＡＰＤ２０９へと導かれ、ＳＬＯ眼底像が得られる。
【００２６】
固視灯２１１から発せられた光は、第３ダイクロイックミラー２０８、ミラー２０７を透過し、レンズ２０６及び２０５を通り、ＳＬＯ走査手段２０４によって被検眼２００上の任意の位置に所定の形状を作り、被検者の固視を促す。
【００２７】
前眼観察用の光路２５２には、レンズ２１２及び２１３、スプリットプリズム２１４、赤外光を検知する前眼部観察用のＣＣＤ２１５が配置されている。このＣＣＤ２１５は、不図示の前眼部観察用照射光の波長、具体的には９７０ｎｍ付近に感度を持つものである。スプリットプリズム２１４は、被検眼２００の瞳孔と共役な位置に配置されており、被検眼２００に対する測定光学系１００−１のＺ軸方向（光軸方向）の距離を、前眼部のスプリット像として検出できる。
【００２８】
ＯＣＴ光学系の光路２５０は前述の通りＯＣＴ光学系を構成しており、被検眼２００の断層画像を撮影するためのものである。より具体的には、断層画像を形成するための干渉信号を得るものである。ＸＹスキャナ２１６は光を被検眼２００上で走査するためのものであり、図２（ｂ）では１枚のミラーとして図示されているが、実際はＸＹ２軸方向の走査を行うガルバノミラーである。
【００２９】
レンズ２１７及び２１８のうち、レンズ２１７については光カプラー２１９に接続されているファイバー２２４から出射するＯＣＴ光源２２０からの光を、被検眼２００に焦点合わせするために不図示のモータによって駆動される。この焦点合わせによって、被検眼２００からの戻り光は同時にファイバー２２４の先端に、スポット状に結像されて入射されることとなる。次に、ＯＣＴ光源２２０からの光路と参照光学系、分光器の構成について説明する。２２０はＯＣＴ光源、２２１は参照ミラー、２２２は分散補償硝子、２２３はレンズ、２１９は光カプラー、２２４から２２７は光カプラーに接続されて一体化しているシングルモードの光ファイバー、２３０は分光器である。
【００３０】
これらの構成によってマイケルソン干渉計を構成している。ＯＣＴ光源２２０から出射された光は、光ファイバー２２５を通じ、光カプラー２１９を介して光ファイバー２２４側の測定光と、光ファイバー２２６側の参照光とに分割される。測定光は前述のＯＣＴ光学系光路を通じ、観察対象である被検眼２００に照射され、被検眼２００による反射や散乱により同じ光路を通じて光カプラー２１９に到達する。
【００３１】
一方、参照光は光ファイバー２２６、レンズ２２３、測定光と参照光の波長分散を合わせるために挿入された分散補償ガラス２２２を介して参照ミラー２２１に到達し反射される。そして同じ光路を戻り、光カプラー２１９に到達する。
【００３２】
光カプラー２１９によって、測定光と参照光は合波され干渉光となる。
【００３３】
ここで、測定光の光路長と参照光の光路長がほぼ同一となったときに干渉を生じる。参照ミラー２２１は、不図示のモータおよび駆動機構によって光軸方向に調整可能に保持され、測定光の光路長に参照光の光路長を合わせることが可能である。干渉光は光ファイバー２２７を介して分光器２３０に導かれる。
【００３４】
また、偏光調整部２２８、２２９は、各々光ファイバー２２４、２２６中に設けられ、偏光調整を行う。これらの偏光調整部は光ファイバーをループ状に引きまわした部分を幾つか持っている。このループ状の部分をファイバーの長手方向を中心として回転させることでファイバーに捩じりを加え、測定光と参照光の偏光状態を各々調整して合わせることができる。
【００３５】
分光器２３０はレンズ２３２、２３４、回折格子２３３、ラインセンサ２３１から構成される。光ファイバー２２７から出射された干渉光はレンズ２３４を介して平行光となった後、回折格子２３３で分光され、レンズ２３２によってラインセンサ２３１に結像される。
【００３６】
次に、ＯＣＴ光源２２０の周辺について説明する。ＯＣＴ光源２２０は、代表的な低コヒーレント光源であるＳＬＤ（Ｓｕｐｅｒ  Ｌｕｍｉｎｅｓｃｅｎｔ  Ｄｉｏｄｅ）である。中心波長は８５５ｎｍ、波長バンド幅は約１００ｎｍである。ここで、バンド幅は、得られる断層画像の光軸方向の分解能に影響するため、重要なパラメータである。
【００３７】
光源の種類は、ここではＳＬＤを選択したが、低コヒーレント光が出射できればよく、ＡＳＥ（Ａｍｐｌｉｆｉｅｄ  Ｓｐｏｎｔａｎｅｏｕｓ  Ｅｍｉｓｓｉｏｎ）等を用いることができる。中心波長は眼を測定することを鑑みると近赤外光が適する。また、中心波長は得られる断層画像の横方向の分解能に影響するため、なるべく短波長であることが望ましい。双方の理由から中心波長は８５５ｎｍとした。
【００３８】
本実施形態では干渉計としてマイケルソン干渉計を用いたが、マッハツェンダー干渉計を用いても良い。測定光と参照光との光量差に応じて、光量差が大きい場合にはマッハツェンダー干渉計を、光量差が比較的小さい場合にはマイケルソン干渉計を用いることが望ましい。
【００３９】
（画像処理装置の構成）
本実施形態の画像処理装置１０１の構成について図１を用いて説明する。
【００４０】
画像処理装置１０１は断層画像撮影装置１００に接続されたパーソナルコンピュータ（ＰＣ）であり、画像取得部１０１−０１、記憶部１０１−０２、撮影制御部１０１−０３、画像処理部１０１−０４、表示制御部１０１−０５を備える。また、画像処理装置１０１は演算処理装置ＣＰＵが画像取得部１０１−０１、撮影制御部１０１−０３、画像処理部１０１−０４および表示制御部１０１−０５を実現するソフトウェアモジュールを実行することで機能を実現する。本発明はこれに限定されず、例えば画像処理部１０１−０４をＡＳＩＣ等の専用のハードウェアで実現してもよいし、表示制御部１０１−０５をＣＰＵとは異なるＧＰＵ等の専用プロセッサを用いて実現してもよい。また断層画像撮影装置１００と画像処理装置１０１との接続はネットワークを介した構成であってもよい。
【００４１】
画像取得部１０１−０１は断層画像撮影装置１００により撮影されたＳＬＯ眼底像や断層画像の信号データを取得する。また画像取得部１０１−０１は断層画像生成部１０１―１１及びモーションコントラストデータ生成部１０１−１２を有する。断層画像生成部１０１―１１は断層画像撮影装置１００により撮影された断層画像の信号データ（干渉信号）を取得して信号処理により断層画像を生成し、生成した断層画像を記憶部１０１−０２に格納する。
【００４２】
撮影制御部１０１−０３は、断層画像撮影装置１００に対する撮影制御を行う。撮影制御には、断層画像撮影装置１００に対して撮影パラメータの設定に関して指示することや、撮影の開始もしくは終了に関して指示することも含まれる。
【００４３】
画像処理部１０１−０４は、位置合わせ部１０１−４１、合成部１０１−４２、補正部１０１−４３、画像特徴取得部１０１−４４、投影部１０１−４５、解析部１０１−４６を有する。先に述べた画像取得部１０１−０１及び合成部１０１−４２は、取得手段の一例である。合成部１０１−４２はモーションコントラストデータ生成部１０１−１２により生成された複数のモーションコントラストデータを位置合わせ部１０１−４１により得られた位置合わせパラメータに基づいて合成し、合成モーションコントラスト画像を生成する。補正部１０１−４３はモーションコントラスト画像内に生じるプロジェクションアーチファクトを２次元もしくは３次元的に抑制する処理を行う（プロジェクションアーチファクトについてはＳ３０４で説明する）。画像特徴取得部１０１−４４は断層画像から網膜や脈絡膜の層境界、中心窩や視神経乳頭中心の位置を取得する。投影部１０１−４５は画像特徴取得部１０１−４４が取得した層境界の位置に基づく深度範囲でモーションコントラスト画像を投影し、正面モーションコントラスト画像を生成する。解析部１０１−４６は、強調部１０１−４６１、抽出部１０１−４６２、計測部１０１−４６３、比較部１０１−４６４を有し、正面モーションコントラスト画像から血管領域の抽出や計測処理を行う。ここで、解析部１０１−４６は、眼部のモーションコントラスト画像における第１の領域と第１の領域よりも狭い領域を少なくとも含む第２の領域とのうち少なくとも１つに対する解析を実行する解析手段の一例である。また、第２の領域は、セクタ領域の一例である。また、第１の領域は、セクタ領域よりも広い領域（例えば、画像全体）の一例である。また、解析部１０１−４６は、眼部のモーションコントラスト画像の少なくとも一部の領域に対する解析を実行するものであっても良い。強調部１０１−４６１は血管強調処理を実行する。また、抽出部１０１−４６２は血管強調画像に基づいて血管領域を抽出する。さらに、計測部１０１−４６３は抽出された該血管領域や該血管領域を細線化することで取得した血管中心線データを用いて血管密度等の計測値を算出する。比較部１０１−４６４は記憶部１０１−０２もしくは外部記憶部１０２から異なる検査日に取得された同一被検眼の合成モーションコントラスト画像及び付随する計測データを読み込んで経時比較データを生成する。
【００４４】
また、表示制御部１０１−０５は、操作者からの指示に応じて第１の領域と第２の領域とのいずれか一方に対して選択された解析の種類を示す情報を用いて一方が解析された結果を示す画像を表示部１０４に表示させることが好ましい。ここで、選択された解析の種類は、例えば、血管領域の面積に関する血管密度（Ｖｅｓｓｅｌ  Ａｒｅａ  Ｄｅｎｓｉｔｙ；ＶＡＤ）や血管長に関する血管密度（Ｖｅｓｓｅｌ  Ｌｅｎｇｔｈ  Ｄｅｎｓｉｔｙ；ＶＬＤ）等である。また、解析された結果を示す画像は、例えば、モーションコントラスト正面画像の少なくとも一部の領域に対して解析された結果を示す２次元画像である。また、解析された結果を示す２次元画像は、例えば、ＶＡＤマップ、ＶＬＤマップ、ＶＡＤセクタマップ、ＶＬＤセクタマップ、また、これらの各解析マップがモーションコントラスト正面画像に重畳された画像がある。また、互いに同じ種類の複数の解析マップが重畳された画像や、互いに同じ種類の複数の解析マップがモーションコントラスト正面画像に重畳された画像であっても良い。例えば、ＶＡＤセクタマップがＶＡＤマップに重畳された２次元画像、ＶＡＤセクタマップ及びＶＡＤマップがモーションコントラスト画像に重畳された２次元画像、ＶＬＤセクタマップがＶＬＤマップに重畳された２次元画像、ＶＬＤセクタマップ及びＶＬＤマップがモーションコントラスト画像に重畳された２次元画像がある。なお、解析部１０１−４６による解析を実行するタイミングは、操作者からの指示に応じて解析の種類が選択されたタイミングでも良いし、解析の種類が選択される前に想定される解析の種類に対応する解析を事前に完了させておいても良い。
【００４５】
ここで、一方に対する選択の後に操作者からの指示に応じて他方に対して選択された解析の種類と、一方に対して選択された解析の種類とが異なる場合を考える。このとき、表示制御部１０１−０５は、一方が解析された結果を示す画像の表示領域において、一方が解析された結果を示す画像の表示を、他方に対して選択された解析の種類を示す情報を用いて一方及び他方が解析された結果を示す画像の表示に変更する制御を実行することが好ましい。これにより、モーションコントラスト画像上の複数の解析領域に対して実行された解析の結果を示す画像が表示される場合に、複数の解析領域に対する解析の種類として互いに対応する種類の選択が容易に実行されるように構成することができる。例えば、図９（ａ）右側のＭａｐボタン群９０２においてＶＬＤマップが選択された後に、Ｓｅｃｔｏｒボタン群９０３においてＶＡＤセクタが選択された場合に、ＶＬＤマップがモーションコントラスト正面画像に重畳された２次元画像の表示を、ＶＡＤセクタマップ及びＶＡＤマップがモーションコントラスト正面画像に重畳された２次元画像の表示に変更する制御が実行されることが好ましい。これにより、例えば、ＶＡＤセクタマップ及びＶＬＤマップがモーションコントラスト画像に重畳された２次元画像のように、解析の種類が異なる状態で表示されることがない。すなわち、例えば、重畳される複数の解析マップが同じ解析の種類として確実に選択されるため、解析結果の確認を行い易くすることができる。このとき、当然、解析の種類を示す情報の表示についても、変更する制御が実行されることが好ましい。すなわち、表示制御部１０１−０５は、一方に対する選択の後に操作者からの指示に応じて他方に対して選択された解析の種類と、一方に対して選択された解析の種類とが異なる場合には、一方に対して選択された解析の種類を示す情報の表示を、他方に対して選択された解析の種類を示す情報の表示に変更する別の制御を実行することが好ましい。なお、選択された解析の種類を示す情報の表示は、選択された解析の種類が識別可能なように表示部１０４に表示されていれば何でも良く、例えば、図９（ａ）右側のＭａｐボタン群９０２、Ｓｅｃｔｏｒボタン群９０３がある。また、経過観察を目的とするフォローアップ撮影において、複数の検査日に対応する複数のモーションコントラスト画像が時系列に並べて表示される表示領域において、複数のモーションコントラスト画像のうちの１つの画像に対して実行される上記制御が、他の画像に対しても実行されることが好ましい。また、複数の検査日に対応する複数のモーションコントラスト画像のうちの１つの画像に対して選択される解析の種類を示す情報が、他の画像に対しても適用されることが好ましい。これらにより、例えば、経過観察を目的とするフォローアップ撮影での利便性も向上することができる。なお、セクタ領域には、複数の領域に区分けされていることが好ましく、各領域内には、その領域の解析結果を示す値（例えば、その領域の平均値）をその解析の種類の単位が識別可能な状態で表示させることが好ましい。ここで、Ｍａｐボタン群９０２やＳｅｃｔｏｒボタン群９０３において、操作者からの指示に応じて「Ｎｏｎｅ」が選択されると、解析の種類が非選択の状態となることが好ましい。このとき、解析された結果を示す画像が表示領域において非表示の状態となり、モーションコントラスト画像が表示領域に表示されることが好ましい。また、一方に対して選択された解析の種類が操作者からの指示に応じて非選択に変更された場合には、他方に対して選択された解析の種類は変更されないことが好ましい。
【００４６】
また、操作者からの指示に応じて選択された解析の種類を示す情報を用いて、眼部のモーションコントラスト画像における第１の領域に対する解析が実行された場合を考える。このとき、表示制御部１０１−０５は、操作者からの指示に応じてモーションコントラスト画像における第２の領域に対して解析された結果を示す画像の表示が選択された場合には、選択された解析の種類を示す情報を用いて第２の領域に対して解析された結果を示す画像を、第１の領域に対して解析された結果を示す画像に重畳された状態で表示部１０４に表示させるものであっても良い。これにより、例えば、解析の種類としてＶＡＤが選択され、また、セクタ領域の表示として「Ｏｎ」が選択されると、表示領域においてＶＡＤセクタマップがＶＡＤマップに重畳された２次元画像を表示させることができる。このため、例えば、重畳される複数の解析マップが同じ解析の種類として確実に選択されるため、解析結果の確認を行い易くすることができる。
【００４７】
また、表示制御部１０１−０５は、操作者からの指示に応じて表示部１０４に表示される解析結果を示す画像が、解析に適した複数の条件のうちの少なくとも２つの条件が満足されない状態で得られた画像である場合には、少なくとも２つの条件に関する情報を、複数の条件の優先順位に応じて報知する報知手段の一例であっても良い。これにより、解析に適した複数の条件のうち少なくとも２つの条件が満足されない場合であっても、より精度の高い解析結果が得られるように操作者が対処し易くすることができる。ここで、解析に適した複数の条件には、例えば、モーションコントラスト画像が、眼部の同一位置を測定光が走査されるように制御して得た複数の３次元のモーションコントラスト画像が合成された画像である、という条件が、他の条件よりも優先度の高い条件として含まれることが好ましい。これにより、例えば、高画質な画像を用いた解析結果を確認することを操作者にアドバイスすることができる。また、表示制御部１０１−０５は、少なくとも２つの条件に関する情報を表示部１０４に表示させることが好ましい。このとき、表示制御部１０１−０５は、操作者からの指示に応じて選択された解析の種類を示す情報を用いて解析結果を示す画像を、少なくとも２つの条件に関する情報と並べて表示部１０４に表示させることが好ましい。また、少なくとも２つの条件のうちの優先度の高い条件に関する警告メッセージであることが好ましい。例えば、図１２の右下部に記載されているように、「Ａｖｅｒａｇｅｄ  ＯＣＴＡ  ｉｓ  ｒｅｃｏｍｍｅｎｄｅｄ  ｉｎ  ｃａｌｃｕｌａｔｉｎｇ  ＶＡＤ  ｏｒ  ＶＬＤ．」という警告メッセージを、解析結果を示す画像が表示される表示画面の端等に表示すればよい。また、上記警告メッセージは、例えば、解析結果を示す画像が表示される表示領域上の端等に、例えば、解析結果を示す画像に重畳された状態で表示されても良い。これらにより、例えば、表示部１０４の表示画面において、操作者が最も確認したい解析結果を示す画像を表示させつつ、残りのスペースを効果的に用いて操作者に対して解析に適した条件をアドバイスすることができる。もちろん、報知手段は、少なくとも２つの条件に関する情報として、少なくとも２つの条件それぞれに関する警告メッセージを、優先順位で報知しても良い。
【００４８】
また、外部記憶部１０２は、被検眼の情報（患者の氏名、年齢、性別など）と、撮影した画像（断層画像及びＳＬＯ画像・ＯＣＴＡ画像）や合成画像、撮影パラメータ、血管領域や血管中心線の位置データ、計測値、操作者が設定したパラメータを関連付けて保持している。入力部１０３は、例えば、マウス、キーボード、タッチ操作画面などであり、操作者は、入力部１０３を介して、画像処理装置１０１や断層画像撮影装置１００へ指示を行う。
【００４９】
次に、図３を参照して本実施形態の画像処理装置１０１の処理手順を示す。図３は、本実施形態における本システム全体の動作処理の流れを示すフローチャートである。
【００５０】
＜ステップ３０１＞
操作者は過去検査データが保存されている被検眼に関して基準検査を選択する。また、画像処理装置１０１は選択された基準検査と同一撮影条件になるようにＯＣＴＡ撮影の撮影条件を設定する。
【００５１】
本実施形態においては、図７（ａ）に示すような患者画面７００において、操作者が入力部１０３を操作することで患者リスト（Ｐａｔｉｅｎｔ  Ｌｉｓｔ）から被検者７０１を選択する。さらに操作者が該被検者の検査リスト（Ｅｘａｍｉｎａｔｉｏｎ  Ｌｉｓｔ）からフォローアップ検査における基準検査（Ｂａｓｅｌｉｎｅ）を選択することで基準検査を決定する（図７の７０２）。検査セットやスキャンモードの選択については、操作者が基準検査を選択した状態で撮影画面（ＯＣＴＣａｐｔｕｒｅ）を開くことにより画像処理装置１０１がフォローアップ検査セットを選択し、スキャンモードを基準検査と同一のスキャンモードに設定する。本実施形態においては、図７（ｂ）の撮影画面７１０に示すように、検査セットとして「Ｆｏｌｌｏｗ−ｕｐ」（７１１）、スキャンモードとして「ＯＣＴＡ」モード７１２を選択する。ここで、検査セットとは検査目的別に設定した（スキャンモードを含む）撮像手順や、ＯＣＴ画像やＯＣＴＡ画像の既定の表示法を指す。
【００５２】
画像処理装置は、断層画像撮影装置１００に対して指示するＯＣＴＡ画像の撮影条件を設定する。なお、個々のＯＣＴＡ撮影に関する撮影条件としては以下の１）〜７）に示すような設定項目があり、これらの設定項目を基準検査と同一の値に設定した上で、Ｓ３０２において適宜休憩を挟みながら（同一撮像条件の）ＯＣＴＡ撮影を所定の回数だけ繰り返し実行する。
【００５３】
本実施形態では７）１クラスタあたりのＢスキャン数が４であるようなＯＣＴＡ撮影を３回繰り返すものとする。
１）走査パターン（Ｓｃａｎ  Ｐａｔｔｅｒｎ）
２）走査領域サイズ（Ｓｃａｎ  Ｓｉｚｅ）
３）主走査方向（Ｓｃａｎｎｉｎｇ  Ｄｉｒｅｃｔｉｏｎ）
４）走査間隔（Ｄｉｓｔａｎｃｅ  ｂｅｔｗｅｅｎ  Ｂ−ｓｃａｎｓ）
５）固視灯位置（Ｆｉｘａｔｉｏｎ  Ｐｏｓｉｔｉｏｎ）
６）コヒーレンスゲート位置（Ｃ−Ｇａｔｅ  Ｏｒｉｅｎｔａｔｉｏｎ）
７）１クラスタあたりのＢスキャン数（Ｂ−ｓｃａｎｓ  ｐｅｒ  Ｃｌｓｔｅｒ）
【００５４】
＜ステップ３０２＞
操作者は入力部１０３を操作して図７（ｂ）に示す撮影画面７１０中の撮影開始（Ｃａｐｔｕｒｅ）ボタン７１３を押下することにより、Ｓ３０１で指定した撮影条件による繰り返しＯＣＴＡ撮影を開始する。
【００５５】
撮影制御部１０１−０３は断層画像撮影装置１００に対してＳ３０１で操作者が指示した設定に基づいて繰り返しＯＣＴＡ撮影を実施することを指示し、断層画像撮影装置１００が対応するＯＣＴ断層画像を取得する。
【００５６】
なお、本ステップにおいて断層画像撮影装置１００はＳＬＯ画像の取得も行い、ＳＬＯ動画像に基づく追尾処理を実行する。本実施形態において繰り返しＯＣＴＡ撮影における追尾処理に用いる基準ＳＬＯ画像は１回目の繰り返しＯＣＴＡ撮影において設定した基準ＳＬＯ画像とし、全ての繰り返しＯＣＴＡ撮影において共通の基準ＳＬＯ画像を用いる。
【００５７】
本実施形ではクラスタ走査として、中心窩を撮像中心とする３ｘ３ｍｍの矩形領域内を、水平方向を主走査方向として垂直方向（副走査方向）の各位置において４回連続でＢ−ｓｃａｎ撮影する。また、隣接するクラスタ走査の間隔は０．０１ｍｍであり、硝子体側にコヒーレンスゲートを設定してＯＣＴ断層画像を取得している。なお、本実施形態においてはＢ−ｓｃａｎ１枚が３００Ａスキャンで構成されるものとする。
【００５８】
またＯＣＴＡ繰り返し撮影中は、Ｓ３０１で設定した撮影条件に加えて
・左右眼の選択
・追尾処理の実行有無
についても同じ設定値を用いる（変更しない）ものとする。
【００５９】
＜ステップ３０３＞
画像取得部１０１−０１及び画像処理部１０１−０４は、Ｓ３０２で取得されたＯＣＴ断層画像に基づいてモーションコントラスト画像を生成する。
【００６０】
まず断層画像生成部１０１−１１は画像取得部１０１−０１が取得した干渉信号に対して波数変換及び高速フーリエ変換（ＦＦＴ）、絶対値変換（振幅の取得）を行うことで１クラスタ分の断層画像を生成する。
【００６１】
次に位置合わせ部１０１−４１は同一クラスタに属する断層画像同士を位置合わせし、重ねあわせ処理を行う。画像特徴取得部１０１−４４が該重ね合わせ断層画像から層境界データを取得する。本実施形態では層境界の取得法として可変形状モデルを用いるが、任意の公知の層境界取得手法を用いてよい。なお層境界の取得処理は必須ではなく、例えばモーションコントラスト画像の生成を３次元のみで行い、深度方向に投影した２次元のモーションコントラスト画像を生成しない場合には層境界の取得処理は省略できる。モーションコントラストデータ生成部１０１−１２が同一クラスタ内の隣接する断層画像間でモーションコントラストを算出する。本実施形態では、モーションコントラストとして脱相関値Ｍｘｙを以下の式（１）に基づき求める。
【００６２】
【数１】
2019154716.tif 000003
【００６３】
ここで、Ａｘｙは断層画像データＡの位置（ｘ，ｙ）における（ＦＦＴ処理後の複素数データの）振幅、Ｂｘｙは断層データＢの同一位置（ｘ，ｙ）における振幅を示している。０≦Ｍｘｙ≦１であり、両振幅値の差異が大きいほど１に近い値をとる。式（１）のような脱相関演算処理を（同一クラスタに属する）任意の隣接する断層画像間で行い、得られた（１クラスタあたりの断層画像数−１）個のモーションコントラスト値の平均を画素値として持つ画像を最終的なモーションコントラスト画像として生成する。
【００６４】
なお、ここではＦＦＴ処理後の複素数データの振幅に基づいてモーションコントラストを計算したが、モーションコントラストの計算法は上記に限定されない。例えば複素数データの位相情報に基づいてモーションコントラストを計算してもよいし、振幅と位相の両方の情報に基づいてモーションコントラストを計算してもよい。あるいは、複素数データの実部や虚部に基づいてモーションコントラストを計算してもよい。
【００６５】
また、本実施形態ではモーションコントラストとして脱相関値を計算したが、モーションコントラストの計算法はこれに限定されない。例えば二つの値の差分に基づいてモーションコントラストを計算しても良いし、二つの値の比に基づいてモーションコントラストを計算してもよい。
【００６６】
さらに、上記では取得された複数の脱相関値の平均値を求めることで最終的なモーションコントラスト画像を得ているが、本発明はこれに限定されない。例えば取得された複数の脱相関値の中央値、あるいは最大値を画素値として持つ画像を最終的なモーションコントラスト画像として生成しても良い。
【００６７】
＜ステップ３０４＞
画像処理部１０１−０４は、繰り返しＯＣＴＡ撮影を通して得られたモーションコントラスト画像群（図８（ａ））を３次元的に位置合わせし、加算平均することで図８（ｂ）に示すように高コントラストな合成モーションコントラスト画像を生成する。なお、合成処理は単純加算平均に限定されない。例えば各モーションコントラスト画像の輝度値に対して任意の重みづけをした上で平均した値でもよいし、中央値をはじめとする任意の統計値を算出してもよい。また位置合わせ処理を２次元的に行う場合も本発明に含まれる。
【００６８】
なお、合成部１０１−４２が合成処理に不適なモーションコントラスト画像が含まれているか否かを判定した上で、不適と判定したモーションコントラスト画像を除いて合成処理を行うよう構成してもよい。例えば、各モーションコントラスト画像に対して評価値（例えば脱相関値の平均値や、ｆＳＮＲ）が所定の範囲外である場合に、合成処理に不適と判定すればよい。
【００６９】
本実施形態では合成部１０１−４２がモーションコントラスト画像を３次元的に合成した後、補正部１０１−４３がモーションコントラスト画像内に生じるプロジェクションアーチファクトを３次元的に抑制する処理を行う。
【００７０】
ここで、プロジェクションアーチファクトは網膜表層血管内のモーションコントラストが深層側（網膜深層や網膜外層・脈絡膜）に映り込み、実際には血管の存在しない深層側の領域に高い脱相関値が生じる現象を指す。図８（ｃ）に、３次元ＯＣＴ断層画像上に３次元モーションコントラストデータを重畳表示した例を示す。網膜表層血管領域に対応する高い脱相関値を持つ領域８０１の深層側（視細胞層）に、高い脱相関値を持つ領域８０２が生じている。本来視細胞層に血管は存在しないにもかかわらず、網膜表層で生じている血管影の明滅が視細胞層に映り込み、視細胞層の輝度値が変化することでアーチファクト８０２が生じる。
【００７１】
補正部１０１−４３は、３次元の合成モーションコントラスト画像上に生じたプロジェクションアーチファクト８０２を抑制する処理を実行する。任意の公知のプロジェクションアーチファクト抑制手法を用いてよいが、本実施形態ではＳｔｅｐ−ｄｏｗｎ  Ｅｘｐｏｎｅｎｔｉａｌ  Ｆｉｌｔｅｒｉｎｇを用いる。Ｓｔｅｐ−ｄｏｗｎ  Ｅｘｐｏｎｅｎｔｉａｌ  Ｆｉｌｔｅｒｉｎｇでは、３次元モーションコントラスト画像上の各Ａスキャンデータに対して式（２）で表される処理を実行することにより、プロジェクションアーチファクトを抑制する。
【００７２】
【数２】
2019154716.tif 000004
【００７３】
ここで、γは負の値を持つ減衰係数、Ｄ（ｘ，ｙ，ｚ）はプロジェクションアーチファクト抑制処理前の脱相関値、ＤＥ（ｘ，ｙ，ｚ）は該抑制処理後の脱相関値を表す。
【００７４】
図８（ｄ）にプロジェクションアーチファクト抑制処理後の３次元合成モーションコントラストデータ（灰色）を断層画像上に重畳表示した例を示す。プロジェクションアーチファクト抑制処理前（図８（ｃ））に視細胞層上に見られたアーチファクトが、該抑制処理によって除去されたことがわかる。
【００７５】
次に、投影部１０１−４５はＳ３０３で画像特徴取得部１０１−４４が取得した層境界の位置に基づく深度範囲でモーションコントラスト画像を投影し、正面モーションコントラスト画像を生成する。任意の深度範囲で投影してよいが、本実施形態においては網膜表層及び網膜深層の深度範囲で２種類の２次元合成モーションコントラスト画像を生成する。また、投影法としては最大値投影（ＭＩＰ；  Ｍａｘｉｍｕｍ  Ｉｎｔｅｎｓｉｔｙ  Ｐｒｏｊｅｃｔｉｏｎ）・平均値投影（ＡＩＰ；  Ａｖｅｒａｇｅ  Ｉｎｔｅｎｓｉｔｙ  Ｐｒｏｊｅｃｔｉｏｎ）のいずれかを選択でき、本実施形態では最大値投影で投影するものとする。
【００７６】
最後に、画像処理装置１０１は取得した画像群（ＳＬＯ画像や断層画像）と該画像群の撮影条件データや、生成したモーションコントラスト画像と付随する生成条件データを検査日時、披検眼を同定する情報と関連付けて外部記憶部１０２へ保存する。
【００７７】
＜ステップ３０５＞
表示制御部１０１−０５は、Ｓ３０３で生成した断層画像や、Ｓ３０４で合成したモーションコントラスト画像、撮影条件や合成条件に関する情報を表示部１０４に表示させる。
【００７８】
図８（ｅ）にレポート画面８０３の例を示す。本実施形態ではＳＬＯ画像及び断層画像、Ｓ３０４で合成及び投影することにより生成した異なる深度範囲の正面モーションコントラスト画像、対応する正面ＯＣＴ画像を表示する。
【００７９】
正面モーションコントラスト画像の投影範囲はリストボックスに表示された既定の深度範囲セット（８０５及び８０９）から操作者が選択することで変更できる。また、投影範囲の指定に用いる層境界の種類とオフセット位置を８０６や８１０のようなユーザインターフェースから変更したり、断層像上に重畳された層境界データ（８０７及び８１１）を入力部１０３から操作して移動させたりすることで投影範囲を変更できる。
【００８０】
さらに、画像投影法やプロジェクションアーチファクト抑制処理の有無を例えばコンテキストメニューのようなユーザインターフェースから選択することにより変更してもよい。
【００８１】
＜ステップ３０６＞
操作者が入力部１０３を用いてＯＣＴＡ計測処理の開始を指示する。
【００８２】
本実施形態では図８（ｅ）のレポート画面８０３のモーションコントラスト画像上をダブルクリックすることで、図９（ａ）のようなＯＣＴＡ計測画面に移行する。モーションコントラスト画像が拡大表示され、適宜画像投影法の種類（最大値投影（ＭＩＰ）もしくは平均値投影（ＡＩＰ））や投影深度範囲、プロジェクションアーチファクト除去処理の実施の有無を選択する。次に、図９（ａ）右側のＭａｐボタン群９０２、Ｓｅｃｔｏｒボタン群９０３、Ｍｅａｓｕｒｅｍｅｎｔボタン９０４経由で表示される選択画面９０５、の中から操作者が適切な項目を選択することで計測の種類と対象領域を選択し、解析部１０１−４６が計測処理を開始する。なお、ＯＣＴＡ計測画面を表示させた時点では、いずれの計測対象領域も設定されていない（Ｍａｐボタン群９０２、Ｓｅｃｔｏｒボタン群ともＮｏｎｅが選択され、選択画面９０５が表示されていない状態）。
【００８３】
計測処理の種類として、本実施形態ではＭａｐボタン群９０２もしくはＳｅｃｔｏｒボタン群９０３から
ｉ）Ｎｏｎｅ（計測しない）
ｉｉ）ＶＡＤ（血管が占める面積に基づいて算出される血管密度）
ｉｉｉ）ＶＬＤ（血管の長さの総和に基づいて算出される血管密度）
のいずれかを選択する。これに限らず、例えば血管構造の複雑さを定量化するＦｒａｃｔａｌ  Ｄｉｍｅｎｓｉｏｎや、血管径の分布（血管の瘤や狭窄の分布）を表すＶｅｓｓｅｌ  Ｄｉａｍｅｔｅｒ  Ｉｎｄｅｘを選択可能に構成してもよい。また、Ｍｅａｓｕｒｅｍｅｎｔボタン９０４経由で表示される選択画面９０５からは
ｉ）任意の２点間の距離計測
ｉｉ）無血管領域の面積計測
ｉｉｉ）ＶＡＤ
ｉｖ）ＶＬＤ
のいずれかを選択できる。
【００８４】
計測処理の対象領域として、本実施形態ではＭａｐボタン群９０２の中からＮｏｎｅ以外を選択することにより画像全体を設定でき、Ｓｅｃｔｏｒボタン群９０３の中からＮｏｎｅ以外を選択することによりセクタ領域（固視位置を中心とし、半径の異なる複数の同心円と該固視位置を通る角度の異なる複数の直線で規定される最小円領域及び扇形領域）を設定できる。また、Ｍｅａｓｕｒｅｍｅｎｔボタン９０４経由で表示される選択画面９０５から所望の計測の種類を選択した上で入力部１０３を用いてモーションコントラスト画像上で任意形状の境界位置（図１０（ｂ）の灰色線部１００１）を指定してＯＫボタンを押下することで、任意形状の計測対象領域を設定できる。該領域内に示された数値は該領域内で計測した値（この場合はＶＡＤ値）を示している。なお、関心領域の手動設定時には指定された境界位置（灰色線部１００１）上に該境界位置が編集可能であることを示す丸い制御点が重畳表示され、ＯＫボタンを押下した時点で該丸い制御点は消えて灰色線部１００１のみとなり、該境界位置は編集不可の状態になる。
【００８５】
本実施形態ではＭａｐボタン群９０２、Ｓｅｃｔｏｒボタン群９０３から各々ＶＡＤを選択することで、ＶＡＤマップ（計測の種類がＶＡＤ、計測対象領域が画像全体）とＶＡＤセクタマップ（計測の種類がＶＡＤ、計測対象領域がＥＴＤＲＳグリッドに対応するセクタ領域）を選択する場合について説明する。
【００８６】
ここで、ＶＡＤはＶｅｓｓｅｌ  Ａｒｅａ  Ｄｅｎｓｉｔｙの略であり、計測対象に含まれる血管領域の割合で定義される血管密度（単位：％）である。すなわち、ＶＡＤは、モーションコントラスト画像において特定される血管領域の面積に関する血管密度の一例である。また、ＶＬＤはＶｅｓｓｅｌ  Ｌｅｎｇｔｈ  Ｄｅｎｓｉｔｙの略であり、単位面積あたりに含まれる血管の長さの総和（単位：ｍｍ−１）で定義される血管密度である。すなわち、ＶＬＤは、モーションコントラスト画像において特定される血管領域の長さに関する血管密度の一例である。また、ＶＡＤやＶＬＤは、モーションコントラスト画像において特定される血管領域に関するパラメータの一例である。血管領域に関するパラメータとしては、血管領域の面積、血管長、血管の曲率等が含まれる。ここで、血管密度は、血管の閉塞範囲や血管網の疎密の程度を定量化するための指標であり、ＶＡＤが最もよく用いられている。ただし、ＶＡＤでは計測値に占める大血管領域の寄与分が大きくなるため、糖尿病網膜症のように毛細血管の病態に注目して計測したい場合には（より毛細血管の閉塞に敏感な指標として）ＶＬＤが用いられる。なお、解析の種類としては、血管領域に関するパラメータ以外にも、例えば、モーションコントラスト画像において特定される無血管領域（Ｎｏｎ  Ｐｅｒｆｕｓｉｏｎ  Ａｒｅａ；ＮＰＡ）に関するパラメータ等がある。無血管領域に関するパラメータとしては、無血管領域の面積や形状（長さや円形度）等が含まれる。また、これに限らず、例えば血管構造の複雑さを定量化するＦｒａｃｔａｌ  Ｄｉｍｅｎｓｉｏｎや、血管径の分布（血管の瘤や狭窄の分布）を表すＶｅｓｓｅｌ  Ｄｉａｍｅｔｅｒ  Ｉｎｄｅｘを計測してもよい。
【００８７】
なお、同一のモーションコントラスト画像に対して複数の計測対象領域を設定してもよい。複数の計測対象領域としては、例えば画像全体、セクタ領域、任意形状領域のうちの少なくとも２つや、２以上の深度範囲、あるいはそれらの組み合わせが挙げられる。該複数の計測対象領域に対して異なる種類の計測が選択された場合には、最後に（指定された解析対象領域に対して）選択された解析の種類を他の計測対象領域に対しても連動して適用した上で計測し、該計測結果を表示するように構成してもよい。例えば、ＶＡＤマップとＶＡＤセクタマップが選択された状態でＶＬＤマップに変更する指示を行うと、ＶＬＤセクタマップが自動で選択され、画像全体に対するＶＬＤ計測とＥＴＤＲＳセクタ領域に対するＶＬＤ計測が実行される。このような連動選択操作により、同一画像に対して異なる種類の計測値が重畳されて操作者が表示内容について誤解・混乱することを防止できる。なお、計測の種類のうちＮｏｎｅ（計測しない）に関しては、各計測対象領域で独立に選択する（ある計測対象領域に対して「Ｎｏｎｅ」が選択された場合に、該「Ｎｏｎｅ」の選択は他の計測対象領域に対して連動して適用されない）。また、本発明は最後に選択された解析の種類を全ての解析対象領域に対して連動して適用し、計測・表示することに限定されるものではない。例えば面内方向の複数の計測対象領域に対しては連動させて深度方向の複数の対象領域に対しては連動させない、あるいはその逆（面内方向の複数の計測対象領域に対しては連動させず、深度方向の複数の対象領域に対しては連動させる）になるように適応的に適用して計測し、対応する計測結果を表示するよう構成してもよい。
【００８８】
次に、解析部１０１−４６は、計測処理の前処理として画像拡大及びトップハットフィルタ処理を行う。トップハットフィルタを適用することにより、背景成分の輝度ムラを軽減できる。本実施形態では合成モーションコントラスト画像の画素サイズが約３μｍになるようにＢｉｃｕｂｉｃ補間を用いて画像拡大し、円形の構造要素を用いてトップハットフィルタ処理を行うものとする。
【００８９】
＜ステップ３０７＞
解析部１０１−４６が血管領域の特定処理を行う。本実施形態では、強調部１０１−４６１がヘシアンフィルタ及びエッジ選択鮮鋭化に基づく血管強調処理を行う。次に抽出部１０１−４６２が２種類の血管強調画像を用いて２値化処理を行い、整形処理を行うことで血管領域を特定する。
【００９０】
血管領域特定処理の詳細はＳ５１０〜Ｓ５６０で説明する。
【００９１】
＜ステップ３０８＞
計測部１０１−４６３が、操作者により指定された計測対象領域に関する情報に基づいて単検査の画像に対する血管密度の計測を行う。引き続いて表示制御部１０１−０５が、計測結果を表示部１０４に表示する。
【００９２】
血管密度としてはＶＡＤとＶＬＤの２種類の指標があり、本実施形態ではＶＡＤを計算する場合の手順を例に説明する。なお、ＶＬＤを計算する場合の手順についても後述する。
【００９３】
操作者が入力部１０３から血管領域もしくは血管中心線データを修正するよう指示を入力した場合には、解析部１０１−４６が操作者から入力部１０３を介して指定された位置情報に基づいて血管領域もしくは血管中心線データを修正し、計測値を再計算する。
【００９４】
なお、本ステップで所定の条件を満たさずに計測が実施された場合には、表示制御部１０１−０５が所定の計測条件を満たした状態で計測すべき旨のメッセージ（警告表示）を表示部１０４に出力する。
【００９５】
ＶＡＤ計測処理の詳細はＳ８１０〜Ｓ８３０、ＶＬＤ計測処理の詳細はＳ８４０〜Ｓ８７０で各々説明する。
【００９６】
＜ステップ３０９＞
解析部１０１−４６は、Ｓ３０７で特定した血管領域や血管中心線のデータを修正するか否かの指示を外部から取得する。この指示は例えば入力部１０３を介して操作者により入力される。修正処理が指示された場合はＳ３０８へ、修正処理が指示されなかった場合はＳ３１０へと処理を進める。
【００９７】
＜ステップ３１０＞
比較部１０１−４６４は経時変化計測（Ｐｒｏｇｒｅｓｓｉｏｎ計測）処理を実行する。図１１にＰｒｏｇｒｅｓｓｉｏｎ計測レポートの例を示す。Ｐｒｏｇｒｅｓｓｉｏｎモードタブ１１０１を指定することにより、該レポート画面が表示され、Ｓ３０６で選択した計測の種類及び計測対象領域に基づく経時変化計測処理が開始される。本実施形態ではＰｒｏｇｒｅｓｓｉｏｎ計測対象画像として検査日が新しいものから順に４検査を自動選択するものとする。これに限らず、例えば最も古い検査日の画像と最新検査の画像、両検査日間に撮影された画像で検査間隔が略等しい画像を選択するよう構成してもよい。ここで、計測対象画像の選択条件は優先度が高い順に
ｉ）固視位置が同じ画像であること
ｉｉ）略同一位置における取得断層像数が大きい（例えば４以上）モーションコントラスト画像か、同等のモーションコントラスト画像になるようにＯＣＴＡ重ね合わせ処理を実施して得られる合成モーションコントラスト画像であること
であり、上記選択条件を満たす画像を優先して選択するものとする。
【００９８】
なお、計測対象画像はこれに限らず、例えば図１１の選択ボタン１１０７を選択して選択画面を表示させ、該選択画面上に表示された画像リストから選択してもよい。
【００９９】
次に、比較部１０１−４６４はＳ３０９で実施した単検査の計測内容に対応する過去検査の画像及び計測値に関するデータを外部記憶部１０２から取得する。さらに位置合わせ部１０１−４１はＳ３０８で計測した単検査の画像と過去検査画像間の位置合わせを行い、比較部１０１−４６４が共通した領域に関する計測データ（計測値、計測値マップ、差分マップ、トレンドグラフのうちの少なくとも一つ）を生成する。なお、差分マップは例えば図１１の１１０８に示すような「Ｓｈｏｗ  Ｄｉｆｆｅｒｅｎｃｅ」チェックボックスを指定することにより生成する。
【０１００】
＜ステップ３１１＞
表示制御部１０１−０５は、Ｓ３１０で実施したＰｒｏｇｒｅｓｓｉｏｎ計測に関するレポートを表示部１０４に表示する。
【０１０１】
本実施形態では、図１１に示すＰｒｏｇｒｅｓｓｉｏｎ計測レポートの上段に網膜表層において計測したＶＡＤマップとＶＡＤセクタマップを重畳表示し、下段に網膜深層において計測したＶＡＤマップとＶＡＤセクタマップを重畳表示する。これにより異なる深度位置の血管病態の時系列変化を一覧して把握でき、図１１の時系列並置表示されたＶＡＤ計測結果では、網膜深層に初期病変が生じ、時間の経過とともに網膜表層へ、また中心窩から傍中心窩へ血管閉塞が広がっていく様子が一覧して把握できる。
【０１０２】
また各計測対象画像に関して、略同一位置における断層画像数やＯＣＴＡ重ね合わせ処理の実施条件に関する情報、ＯＣＴ断層画像もしくはモーションコントラスト画像の評価値（画質指標）に関する情報を表示部１０４に表示させてもよい。図１１ではＯＣＴＡ重ね合わせ処理実施済であることを示す印（左上の「Ａｖｅｒａｇｅｄ  ＯＣＴＡ」）が表示されている。なお、図１１の上部に表示されている矢印１１０４は現在選択されている検査であることを示す印であり、基準検査（Ｂａｓｅｌｉｎｅ）はＦｏｌｌｏｗ−ｕｐ撮影の際に選択した検査（図１１の一番左画像の画像）である。もちろん、基準検査を示すマークを表示部１０４に表示させてもよい。またＳ３１０で「Ｓｈｏｗ  Ｄｉｆｆｅｒｅｎｃｅ」チェックボックス１１０８が指定された場合には、基準画像上に基準画像に対する計測値分布（マップもしくはセクタマップ）を表示し、それ以外の検査日に対応する領域に基準画像に対して算出した計測値分布との差分計測値マップを表示する。計測結果としてはレポート画面上にトレンドグラフ（経時変化計測によって得られた各検査日の画像に対する計測値のグラフ）を表示させてもよい。該トレンドグラフの回帰直線（曲線）や対応する数式をレポート画面に表示させてもよい。
【０１０３】
なお、本実施形態では異なる深度範囲として網膜表層及び網膜深層の画像及び計測値を時系列で表示したが、これに限らず例えば網膜表層・網膜深層・網膜外層・脈絡膜の４種類の深度範囲の画像及び計測値を時系列で表示してもよい。
【０１０４】
あるいは、異なる指標の計測結果を並置して時系列で表示してもよい。例えば、上段にＶＡＤマップの時系列表示、下段にＶＬＤマップ（もしくは無血管領域の形状値）の時系列並置表示を行ってもよい。
【０１０５】
なお、時系列並置表示する場合の投影深度範囲は、Ｓ３０５で説明した図８（ｅ）のユーザインターフェース（８０５・８０６や８０９・８１０）の場合と同様に、図１１の１１０２・１１０５・１１０６のユーザインターフェースを用いて変更できる。また、同様に投影方法（ＭＩＰ／ＡＩＰ）やプロジェクションアーチファクト抑制処理についても例えばコンテキストメニューから選択するなどの方法により変更してもよい。さらに、ショートカットメニューから計測の種類及び計測対象領域に関する項目を異なる値に変更することにより、Ｐｒｏｇｒｅｓｓｉｏｎ計測の種類及び計測対象領域を変更して再計測できる。
【０１０６】
例えば、ショートカットメニュー上に図９（ａ）のＭａｐボタン群９０２に相当する項目及びＳｅｃｔｏｒボタン群９０３に相当する項目を表示しておき、各々１項目ずつ（例えば「ＶＬＤ  Ｍａｐ」と「ＶＬＤ  Ｓｅｃｔｏｒ」を）選択する。なお、Ｓ３０６の場合と同様に、複数の計測対象領域が選択され、かつ一方の領域に対する計測の種類が変更された場合は他方の領域に対しても同じ計測の種類が連動して適用され、計測が実行される。また計測対象領域を設定しない（「Ｎｏｎｅ」を選択する）指示に関しては、他方の計測対象領域に対して連動して適用されないものとする。
【０１０７】
また、表示部１０４に表示されたモーションコントラスト画像や、抽出部１０１−４６２や計測部１０１−４６３によって生成された血管領域や血管中心線に関する２値画像、計測値や計測マップを外部記憶部１０２にファイルとして出力し、保存してもよい。なお、比較観察を容易にするため、ファイル出力するモーションコントラスト画像、血管領域や血管中心線に関する２値画像、計測値マップの画像サイズや画素サイズは同一にすることが望ましい。
【０１０８】
さらに、単検査に対する計測の場合と同様の方法で（詳細はＳ８３０で説明する）、所定の条件を満たさない状態で実施された計測結果を計測レポート画面に表示する際には警告メッセージを表示してよい。なお推奨条件としてはＳ８３０で示した条件に限定されるものではなく、例えば「選択した経時変化計測対象画像間で略同一位置における取得断層画像数もしくはモーションコントラスト画像の合成条件もしくは画質指標値の相違が所定値未満であること」を推奨条件として設定し、該条件を満たさない場合に警告表示してもよい。
【０１０９】
＜ステップ３１２＞
画像処理装置１０１はＳ３０１からＳ３１２に至る一連の処理を終了するか否かの指示を外部から取得する。この指示は入力部１０３を介して操作者により入力される。処理終了の指示を取得した場合は処理を終了する。一方、処理継続の指示を取得した場合にはＳ３０２に処理を戻し、次の披検眼に対する処理（または同一披検眼に対する再処理を）行う。
【０１１０】
さらに、図５（ａ）に示すフローチャートを参照しながら、Ｓ３０７で実行される処理の詳細について説明する。
【０１１１】
＜ステップ５１０＞
強調部１０１−４６１は、ステップ３０６の前処理を実施されたモーションコントラスト画像に対してヘッセ行列の固有値に基づく血管強調フィルタ処理を行う。このような強調フィルタはヘシアンフィルタと総称され、例えばＶｅｓｓｅｌｎｅｓｓ  ｆｉｌｔｅｒやＭｕｌｔｉ−ｓｃａｌｅ  ｌｉｎｅ  ｆｉｌｔｅｒが挙げられる。本実施形態ではＭｕｌｔｉ−ｓｃａｌｅ  ｌｉｎｅ  ｆｉｌｔｅｒを用いるが、任意の公知の血管強調フィルタを用いてよい。
【０１１２】
ヘシアンフィルタは強調したい血管の径に適したサイズで画像を平滑化した上で、該平滑化画像の各画素において輝度値の２次微分値を要素として持つヘッセ行列を算出し、該行列の固有値の大小関係に基づいて局所構造を強調する。ヘッセ行列は式（３）で与えられるような正方行列であり、該行列の各要素は例えば式（４）に示すような画像の輝度値Ｉを平滑化した画像の輝度値Ｉｓの２次微分値で表される。ヘシアンフィルタでは、このようなヘッセ行列の「固有値（λ１、λ２）の一方が０に近く、他方が負かつ絶対値が大きい」場合に線状構造とみなして強調する。これはモーションコントラスト画像上の血管領域が持つ特徴、すなわち「走行方向では輝度変化が小さく、走行方向に直交する方向では輝度値が大きく低下する」が成り立つ画素を線状構造とみなして強調することに相当する。
【０１１３】
またモーションコントラスト画像には毛細血管から細動静脈まで様々な径の血管が含まれることから、複数のスケールでガウスフィルタにより平滑化した画像に対してヘッセ行列を用いて線強調画像を生成する。次に式（５）に示すようにガウスフィルタの平滑化パラメータσの二乗を補正係数として乗じた上で最大値演算により合成し、該合成画像Ｉｈｅｓｓｉａｎをヘシアンフィルタの出力とする。
【０１１４】
ヘシアンフィルタはノイズに強く、血管の連続性が向上するという利点がある。一方で実際には事前に画像に含まれる最大血管径が不明の場合が多いため、特に平滑化パラメータが画像中の最大血管径に対して大きすぎる場合に強調された血管領域が太くなりやすいという欠点がある。
【０１１５】
そこで、本実施形態ではＳ５３０で述べる別の血管強調手法で血管領域を強調した画像と演算することにより、血管領域が太くなりすぎることを防止する。
【０１１６】
【数３】
2019154716.tif 000005
【０１１７】
【数４】
2019154716.tif 000006
【０１１８】
【数５】
2019154716.tif 000007
【０１１９】
＜ステップ５２０＞
抽出部１０１−４６２は、Ｓ５１０で生成したヘシアンフィルタによる血管強調画像（以下、ヘシアン強調画像と表記）を２値化する。
【０１２０】
ヘシアン強調画像の輝度統計値（平均値や中央値等）を閾値として２値化する場合、例えば視神経乳頭部では大血管の高輝度領域の影響を受けて閾値が高くなり、乳頭周囲毛細血管（ＲＰＣ；  Ｒａｄｉａｌ  Ｐｅｒｉｐａｐｉｌｌａｒｙ  Ｃａｐｉｌｌａｒｙ）の抽出不足が生じる場合がある。また、網膜深層のような無血管領域の拡大が生じやすい領域においては、閾値が低すぎて無血管領域を血管として誤検出する場合がある。
【０１２１】
そこで、本実施形態では低スケールの強調画像のみで合成したヘシアン強調画像の平均値を閾値とすることで視神経乳頭部において閾値が高くなりすぎることを防ぐ。また該閾値の下限値を設定することにより無血管領域における誤検出を抑制する。
【０１２２】
ここで、視神経乳頭部で閾値が高くなりすぎるのを防止する方法としては、低スケールの強調画像の統計値を閾値として２値化することに限定されない。例えば、ヘシアン強調画像上の輝度値が所定値以上である場合に所定値であるとみなした場合の平均値を閾値として２値化しても同等の効果が期待できる。あるいは、例えばＭ−ｅｓｔｉｍａｔｏｒのようなロバスト推定値を閾値として２値化してもよい。
【０１２３】
また、本実施形態では合成モーションコントラスト画像をヘシアンフィルタで強調しているため、単独のモーションコントラスト画像をヘシアンフィルタで強調する場合に比べて２値化した血管領域の連続性がさらに向上する。
【０１２４】
＜ステップ５３０＞
強調部１０１−４６１は、Ｓ３０６で生成したトップハットフィルタ適用後の合成モーションコントラスト画像に対してエッジ選択鮮鋭化処理を行う。
【０１２５】
ここで、エッジ選択鮮鋭化処理とは画像のエッジ部分に重みを大きく設定した上で重みづけ鮮鋭化処理を行うことを指す。本実施形態では、前記合成モーションコントラスト画像に対してＳｏｂｅｌフィルタを適用した画像を重みとしてアンシャープマスク処理を行うことにより、エッジ選択鮮鋭化処理を実施する。
【０１２６】
小さなフィルタサイズで鮮鋭化処理を実施すると、細い血管のエッジが強調され２値化した際により正確に血管領域を特定できる（血管領域が太くなる現象を防止できる）。一方で特に同一撮影位置での断層像数が少ないモーションコントラスト画像の場合にはノイズが多いため、特に血管内のノイズも一緒に強調してしまう恐れがある。そこで、エッジ選択鮮鋭化を行うことによってノイズの強調を抑制する。
【０１２７】
＜ステップ５４０＞
抽出部１０１−４６２は、Ｓ５３０で生成したエッジ選択鮮鋭化処理を適用した鮮鋭化画像を２値化する。任意の公知の２値化法を用いてよいが、本実施形態では該鮮鋭化画像上の各局所領域内で算出した輝度統計値（平均値もしくは中央値）を閾値として２値化する。
【０１２８】
ただし、視神経乳頭部の大血管領域においては設定される閾値が高すぎて２値画像上の血管領域内に多数の穴が空くため、上記閾値の上限値を設定することにより特に視神経乳頭部において閾値が高くなりすぎるのを防止する。
【０１２９】
また、Ｓ５２０の場合と同様に画像中に無血管領域が占める割合が大きい場合には閾値が低すぎて無血管領域の一部を血管として誤検出する場合が生じる。そこで、上記閾値の下限値を設定することにより誤検出を抑制する。
【０１３０】
なお、Ｓ５２０の場合と同様に合成モーションコントラスト画像をエッジ選択鮮鋭化しているため、単独のモーションコントラスト画像をエッジ選択鮮鋭化する場合に比べて２値化した場合のノイズ状の誤検出領域をより減らすことができる。
【０１３１】
＜ステップ５５０＞
抽出部１０１−４６２は、Ｓ５２０で生成したヘシアン強調画像の２値画像の輝度値と、Ｓ５４０で生成したエッジ選択鮮鋭化画像の２値画像の輝度値の双方が０より大きい場合に血管候補領域として抽出する。当該演算処理により、ヘシアン強調画像に見られる血管径を過大評価している領域と、エッジ選択鮮鋭化画像上に見られるノイズ領域がともに抑制され、血管の境界位置が正確かつ血管の連続性が良好な２値画像を取得できる。
【０１３２】
また双方の２値画像とも合成モーションコントラスト画像に基づく２値画像であることから、単独のモーションコントラスト画像に基づく２値画像に比べて２値化した場合のノイズ状の誤検出領域が減少するとともに、特に毛細血管領域の連続性が向上する。また、合成モーションコントラスト画像であることから検査間の画質や輝度レベルも安定し、血管の抽出能が検査間で安定しやすい。
【０１３３】
＜ステップ５６０＞
抽出部１０１−４６２は、血管領域の整形処理として２値画像のオープニング処理（収縮処理後に膨張処理を行うこと）及びクロージング処理（膨張処理後に収縮処理を行うこと）を実施する。なお、整形処理はこれに限らず例えば２値画像をラベリングした場合の各ラベルの面積に基づく小領域除去を行ってもよい。
【０１３４】
なお、様々な径の血管が含まれるモーションコントラスト画像において血管強調する際のスケールを適応的に決定する方法はＳ５１０〜Ｓ５６０に述べた方法に限定されない。例えば、図５（ｂ）のＳ６１０〜Ｓ６５０に示すように、ヘシアン強調画像の輝度値とエッジ選択鮮鋭化による血管強調画像の輝度値を乗じる演算を適用した画像に対する輝度統計値（例えば平均値）を閾値として２値化することによって血管領域を特定してもよい。該閾値には下限値や上限値を設定できる。
【０１３５】
あるいは、図５（ｃ）のＳ７１０〜Ｓ７４０に示すように、画像の固視位置や深度範囲によってヘシアンフィルタ適用時の平滑化パラメータσの範囲を適応的に変えた上でヘシアンフィルタを適用し、２値化することによって血管強調してもよい。例えば、乳頭部網膜表層ではσ＝１〜１０、黄斑部網膜表層ではσ＝１〜８、黄斑部網膜深層ではσ＝１〜６のように設定できる。
【０１３６】
また、２値化処理は閾値処理に限定されるものではなく任意の公知のセグメンテーション手法を用いてよい。
【０１３７】
さらに、図６（ａ）に示すフローチャートを参照しながら、Ｓ３０８で実行される処理の詳細について説明する。
【０１３８】
＜ステップ８１０＞
操作者は、入力部１０３を介して計測処理における関心領域を設定する。本実施形態ではＳ３０６で計測内容（計測の種類と対象領域）としてＶＡＤマップ（計測の種類がＶＡＤで、計測対象領域が画像全体）とＶＡＤセクタマップ（計測の種類がＶＡＤで、計測対象領域がＥＴＤＲＳグリッドに対応するセクタ領域）が選択されている。従って、関心領域として（ｉ）画像全体（ｉｉ）固視灯位置を中心とするセクタ領域（直径１ｍｍの内円と直径３ｍｍの外円で規定される環状領域内をＳｕｅｒｉｏｒ・Ｉｎｆｅｒｉｏｒ・Ｎａｓａｌ・Ｔｅｍｐｏｒａｌの４つの扇形に分割した領域及び該内円領域）を設定する。
【０１３９】
＜ステップ８２０＞
計測部１０１−４６３は、Ｓ３０７で得られた血管領域の２値画像に基づいて計測処理を行う。本実施形態では、該２値画像の各画素位置において当該画素を中心とした近傍領域内に占める非０画素（白画素）の割合を当該画素における血管密度（ＶＡＤ）として算出する。さらに、各画素で算出した血管密度（ＶＡＤ）の値を持つ画像（ＶＡＤマップ）を生成する。
【０１４０】
また、該２値画像上の（Ｓ８１０で設定した）各セクタ領域における非０画素（白画素）の割合を当該セクタにおける血管密度（ＶＡＤ）として算出する。さらに、各セクタ領域で算出した血管密度（ＶＡＤ）の値を持つマップ（ＶＡＤセクタマップ）を生成する。
【０１４１】
＜ステップ８３０＞
表示制御部１０１−０５は、計測結果としてＳ８２０で生成したＶＡＤマップ及びＶＡＤセクタマップを表示部１０４に表示する。本実施形態では、図９（ｂ）の９０６に網膜表層のＶＡＤマップ、９０８に網膜深層のＶＡＤマップを表示する。さらに９０７に網膜表層のＶＡＤセクタマップ、９０９に網膜深層のＶＡＤセクタマップを重畳表示する。
【０１４２】
なお、本実施形態では図９（ｂ）において実施する計測の推奨条件として
ｉ）選択した計測対象画像の中に略同一位置における取得断層画像数が所定値以上のモーションコントラスト画像もしくは該所定値以上相当の合成モーションコントラスト画像、あるいは画質指標値（Ｑｕａｌｉｔｙ  Ｉｎｄｅｘ）が所定値以上のモーションコントラスト画像に対する計測であること
ｉｉ）最大値投影で生成されたモーションコントラスト画像に対する計測であること
ｉｉｉ）プロジェクションアーチファクト除去処理が実施済であること
ｉｖ）網膜表層、網膜深層、放射状乳頭周囲毛細血管（ＲＰＣ）が含まれる投影深度範囲、のいずれか投影深度範囲で生成されたモーションコントラスト画像に対する計測であること
を設定するものとし、上記ｉ）〜ｉｖ）の少なくとも一つを満たさない状態で実施された計測結果を計測レポート画面に表示する際には（正確に計測することが期待できないような条件で計測が行われたとみなして）警告表示を行う。
【０１４３】
例えば、ｉ）を満たさない状態で計測された結果を表示する際には、図９（ｂ）の右下部に「Ａｖｅｒａｇｅｄ  ＯＣＴＡ  ｉｓ  ｒｅｃｏｍｍｅｎｄｅｄ  ｉｎ  ｃａｌｃｕｌａｔｉｎｇ  ＶＡＤ  ｏｒ  ＶＬＤ．」のように警告メッセージを表示すればよい。（例えば、図１２の右下部）
【０１４４】
また、ｉｉ）を満たさない状態で計測された結果を表示する際には、図９（ｂ）の右下部に「ＭＩＰ  ｉｓ  ｒｅｃｏｍｍｅｎｄｅｄ  ｉｎ  ｃａｌｃｕｌａｔｉｎｇ  ＶＡＤ  ｏｒ  ＶＬＤ．」のように警告メッセージを表示すればよい。
【０１４５】
同様に、ｉｉｉ）を満たさない状態で計測された結果を表示する際には、図９（ｂ）の右下部に「ＰＡＲ  ｉｓ  ｒｅｃｏｍｍｅｎｄｅｄ  ｉｎ  ｃａｌｃｕｌａｔｉｎｇ  ＶＡＤ  ｏｒ  ＶＬＤ．」のように警告メッセージを表示すればよい。
【０１４６】
さらに、ｉｖ）を満たさない状態で計測された結果を表示する際には「Ｓｕｐｅｒｆｉｃｉａｌ  Ｃａｐｉｌｌａｒｙ，  Ｄｅｅｐ  Ｃａｐｉｌｌａｒｙ，  ＲＰＣ  ｃａｎ  ｂｅ  ａｎａｌｙｚｅｄ  ｉｎ  ｃａｌｃｕｌａｔｉｎｇ  ＶＡＤ  ｏｒ  ＶＬＤ．」のように警告メッセージを表示する。該警告メッセージを表示することで、推奨する計測条件を満たさない計測により得られた計測結果は信頼性の低い計測結果となるリスクがあることを周知するとともに、推奨する計測条件を示すことで、より信頼性の高い計測を実施しやすくする。
【０１４７】
なお、レポート画面を多数の警告メッセージが占めるのを避けるために上記推奨条件に優先順位を付け（例えば、ｉ）を最優先、ｉｉ）を２番目、ｉｉｉ）を３番目、ｉｖ）を４番目に重要な条件とし）満たされていない推奨条件のうち最も優先度の高い条件に関する警告表示を行うよう構成してもよい。また図９（ｂ）のように複数の計測結果を表示する場合には、個別の計測に対して各々警告メッセージを表示してもよいし、表示対象の警告メッセージのうち最も優先度の高い警告メッセージのみ表示させるようにしてもよい。あるいは、計測結果の信頼性に影響が大きい条件を理解しやすくしつつ満たされていない条件に関する警告を漏れなく表示するために、満たされていない計測条件に関する警告を優先順位が識別可能な態様で（色や大きさ等を変えて）表示部１０４に表示してもよい。複数の計測結果を表示する場合の例として、レポート画面の上下段に各々計測結果を表示する場合や、同一画像に対し複数の計測対象領域を設定して計測した結果を表示する場合が挙げられる。
【０１４８】
なお、警告メッセージは同一のレポート画面内に表示してもよいし、別画面として表示してもよい。また警告メッセージは文字列に限定されるものではなく静止画像や動画像を表示部１０４に表示したり、音声として出力したりしてもよい。該警告メッセージが表示されたレポート画面をファイル出力したり、印刷出力したりする場合も本発明に含まれる。
【０１４９】
さらに、表示部１０４に表示された警告メッセージのうち消去する警告メッセージを操作者が入力部１０３を用いて選択可能にしたり、警告の優先順位の変更や表示対象外とする警告メッセージの指定を行うためのユーザインターフェースを備えたりしてもよい。
【０１５０】
上記では血管密度としてＶＡＤを計測する場合の手順を例に説明したが、計測値としてＶＬＤマップやＶＬＤセクタマップを生成する場合は、上記Ｓ８１０〜８３０の代わりに図６（ｂ）に示すＳ８４０〜８７０を実行する。
【０１５１】
＜ステップ８４０＞
計測部１０１−４６３は、Ｓ３０７で生成した血管領域の２値画像を細線化処理することにより、血管の中心線に相当する線幅１画素の２値画像（以下、スケルトン画像と表記）を生成する。任意の細線化法もしくはスケルトン処理を用いてよいが、本実施形態では細線化法としてＨｉｌｄｉｔｃｈの細線化法を用いる。
【０１５２】
＜ステップ８５０＞
操作者は、入力部１０３を介してＳ８１０と場合と同様の関心領域を設定する。本実施形態では計測内容（計測の種類と対象領域）としてＶＬＤマップとＶＬＤセクタマップを算出するものとし、Ｓ８１０ではＶＡＤが選択されていたのに対し、本ステップではＶＬＤが選択されていることが異なるだけである。なお、ＶＬＤマップもしくはＶＬＤセクタマップをモーションコントラスト画像上に重畳表示したくない場合は、図９（ａ）のＭａｐもしくはＳｅｃｔｏｒの項目を「Ｎｏｎｅ」に設定すればよい。
【０１５３】
＜ステップ８６０＞
計測部１０１−４６３はＳ８４０で得られたスケルトン画像に基づいて計測処理を行う。本実施形態では、該スケルトン画像の各画素位置において当該画素を中心とした近傍領域における単位面積当たりの非０画素（白画素）の長さの総和［ｍｍ−１］を当該画素における血管密度（ＶＬＤ）として算出する。さらに、各画素で算出した血管密度（ＶＬＤ）の値を持つ画像（ＶＬＤマップ）を生成する。
【０１５４】
また、該スケルトン画像上の（Ｓ８５０で設定した）各セクタ領域における単位面積当たりの非０画素（白画素）の長さの総和［ｍｍ−１］を当該セクタにおける血管密度（ＶＬＤ）として算出する。さらに、各セクタ領域で算出した血管密度（ＶＬＤ）の値を持つマップ（ＶＬＤセクタマップ）を生成する。
【０１５５】
＜ステップ８７０＞
表示制御部１０１−０５は、計測結果としてＳ８６０で生成したＶＬＤマップ及びＶＬＤセクタマップを図９（ｂ）の９０６・９０７もしくは９０８・９０９に表示する。
【０１５６】
なお、Ｓ８３０の場合と同様に、所定の解析に適した条件を満たさない状態で実施された計測結果を計測レポート画面に表示する際には表示部１０４に警告メッセージを表示するものとする。
【０１５７】
また本実施形態では単検査での血管領域特定及び計測結果の表示法として正面モーションコントラスト画像上に計測マップを重畳表示する場合について説明したが、これに限定されるものではない。例えば、特定された血管領域の２値画像やスケルトン画像を図９（ｂ）の９０６や９０８に表示させてもよい。あるいは、９０６や９０８にモーションコントラスト画像を表示しておき、その上に特定された血管領域の２値画像もしくはスケルトン画像を色もしくは透明度を適宜調整した上で重畳表示するよう構成する場合も本発明に含まれる。また２値画像は正面画像として表示することに限定されず、例えばＢスキャン断層画像上に特定された血管領域の２値画像もしくはスケルトン画像を色もしくは透明度を適宜調整した上で重畳表示してもよい。
【０１５８】
なお、Ｓ３０９で操作者が入力部１０３から血管領域もしくは血管中心線データを修正するよう指示を入力した場合には以下のような手順で修正する。
【０１５９】
すなわち、図１０（ａ）に示すような合成モーションコントラスト画像に対して同図（ｃ）に示すような過抽出領域を含む２値画像が得られた場合に、操作者が入力部１０３を経由して指定した位置の白画素を解析部１０１−４６が削除する。追加／削除位置の指定方法の例としては、例えば削除の場合は「ｄ」キーを押しながら、追加の場合は「ａ」キーを押しながらマウスでクリックする方法が挙げられる。あるいは、同図（ｄ）に示すようにモーションコントラスト画像に基づく画像上に修正対象である２値画像（血管領域もしくは血管中心線）の色や透明度を調整して重畳表示し、過抽出もしくは抽出不足の領域が判別しやすい状態にしておく。同図（ｄ）の矩形領域１００２内を拡大した画像を同図（ｅ）に示す。灰色が過抽出した領域で、白色が元のモーションコントラスト画像の脱相関値を示す。該過抽出／抽出不足の領域を操作者が入力部１０３を用いて指定することにより、正確かつ効率的に２値画像上の血管もしく血管中心線領域を修正するよう構成してもよい。なお、２値画像の修正処理は正面画像に限定されない。例えば図９（ａ）の９１０に示すような任意のスライス位置のＢスキャン断層像上にモーションコントラストデータや血管領域の２値データもしくは血管中心線領域を色や透明度の調整後に重畳する。このように過抽出もしくは抽出不足の領域が判別しやすい状態にした上で、操作者が修正（追加／削除）する２値データの３次元位置（ｘ，ｙ，ｚ座標）を入力部１０３により指定して修正してもよい。
【０１６０】
さらに、２値画像（血管領域の２値画像もしくはスケルトン画像）が修正済であることを示す情報もしくは修正位置に関する情報を該２値画像と関連付けて外部記憶部１０２に保存しておき、Ｓ８７０もしくはＳ３１１で血管特定結果及び計測結果を表示部１０４に表示する際に該修正済であることを示す情報もしくは修正位置に関する情報を表示部１０４に表示してもよい。
【０１６１】
また、本実施形態では合成部１０１−４２が繰り返しＯＣＴＡ撮影終了時に合成モーションコントラスト画像を生成する場合について説明したが、合成モーションコントラスト画像の生成手順はこれに限定されない。例えば図８（ｅ）のレポート画面８０３上に合成モーションコントラスト画像生成指示ボタン８１２を配置しておく。ＯＣＴＡ撮影完了後（撮影日より後の日でもよい）に操作者が明示的に該生成指示ボタン８１２を押下した場合に合成部１０１−４２が合成モーションコントラスト画像を生成するよう画像処理装置１０１を構成してもよい。操作者が明示的に合成画像生成指示ボタン８１２を押下して合成画像を生成する場合、図８（ｅ）に示すようなレポート画面８０３上に合成モーションコントラスト画像８０４や合成条件データ、検査画像リスト上に合成画像に関する項目を表示させる。
【０１６２】
また操作者が明示的に該生成指示ボタン８１２を押下する場合は、表示制御部１０１−０５が以下の処理を行う。すなわち、合成対象画像選択画面を表示させ、操作者が入力部１０３を操作して合成対象画像群を指定し、許容ボタンを押下した場合に合成部１０１−４２が合成モーションコントラスト画像を生成し、表示部１０４に表示させる。なお、生成済の合成モーションコントラスト画像を選択して合成する場合も本発明に含まれる。
【０１６３】
また操作者が合成画像生成指示ボタン８１２を押下した場合、３次元モーションコントラスト画像を投影した２次元画像同士を合成することにより２次元合成画像を生成してもよいし、３次元合成画像を生成後に投影することで２次元合成画像を生成してもよい。
【０１６４】
以上述べた構成によれば、画像処理装置１０１は同一被検眼に対して異なる日時に略同一撮影条件で取得したＯＣＴＡ重ね合わせ画像から生成した網膜表層及び網膜深層の正面モーションコントラスト画像を用いて血管領域特定及び血管密度計測処理を実施する。該特定処理及び計測処理により得られた合成画像及び計測値を複数の深度範囲で時系列並置表示する。これにより、例えば、検査ごとのＯＣＴ断層画像の信号強度や画質のばらつきの影響を抑制しながら、血管病態（血管閉塞・新生血管・血管瘤等）の変化を正確に特定・計測できる。また、１枚のＯＣＴＡ画像に対して解析することにより、例えば、血管病変の分布を定量的に把握することができる。
【０１６５】
［第二実施形態］
本実施形態に係る画像処理装置の一つは、第一実施形態における血管領域特定及び計測処理を３次元で実施し、得られた画像及び計測データ（血管領域や血管中心線・計測データ）を時系列で並置表示するよう構成したものである。
【０１６６】
具体的には、脈絡膜新生血管（ＣＮＶ；Ｃｈｏｒｏｉｄａｌ  ＮｅｏＶａｓｕｃｕｌａｒｉｚａｔｉｏｎ）を含む３次元の合成モーションコントラスト画像に対してモーションアーチファクト抑制処理を実施する。次に３次元のモルフォロジーフィルタ及び血管強調フィルタを適用し、２値化することにより３次元でＣＮＶを含む血管領域を特定する。さらに網膜表層及び網膜深層で算出した血管密度と、網膜外層で特定及び計測した脈絡膜新生血管領域の２値画像及び体積値を時系列表示する場合について説明する。
【０１６７】
本実施形態に係る画像処理装置１０１を備える画像処理システム１０の構成及び画像処理フローは第一実施形態の場合と同様であるので省略する。
【０１６８】
また、本実施形態での画像処理フローのうち図３のＳ３０６〜Ｓ３０８、Ｓ３１０〜Ｓ３１１以外は第一実施形態の場合と同様であるので説明は省略する。
【０１６９】
＜ステップ３０６＞
操作者が入力部１０３を用いてＯＣＴＡ計測処理の開始を指示する。
【０１７０】
本実施形態では図８（ｅ）のレポート画面８０３のモーションコントラスト画像上をダブルクリックすることで、図９（ａ）のようなＯＣＴＡ計測画面に移行する。モーションコントラスト画像が拡大表示され、適宜画像投影法の種類（最大値投影（ＭＩＰ）もしくは平均値投影（ＡＩＰ））や投影深度範囲、プロジェクションアーチファクト除去処理の実施の有無を選択する。次に、図９（ａ）右側のＭａｐボタン群９０２、Ｓｅｃｔｏｒボタン群９０３、Ｍｅａｓｕｒｅｍｅｎｔボタン９０４経由で表示される選択画面９０５、の中から操作者が適切な項目を選択することで計測の種類と対象領域を選択し、解析部１０１−４６が計測処理を開始する。
【０１７１】
計測処理の種類として、本実施形態ではＭａｐボタン群もしくはＳｅｃｔｏｒボタン群から
ｉ）Ｎｏｎｅ（計測しない）
ｉｉ）ＶＡＤ（血管が占める面積に基づいて算出される血管密度）
ｉｉｉ）ＶＬＤ（血管の長さの総和に基づいて算出される血管密度）
ｉｖ）Ｖｏｌｕｍｅ（血管領域の体積）
のいずれかを選択する。これに限らず任意の種類の計測を実施してよい。
【０１７２】
例えば、ｉｖ）のＶｏｌｕｍｅの代わりに（２次元モーションコントラスト画像上で特定したり、特定した３次元血管領域を所定の深度範囲（例えば網膜外層）で投影したりして得られる）血管領域（例えば脈絡膜毛細血管）の面積を計測する場合も本発明に含まれる。
【０１７３】
また、Ｍｅａｓｕｒｅｍｅｎｔボタン経由で表示される選択画面からは
ｉ）無血管領域の面積計測
ｉｉ）血管密度（ＶＡＤ）
ｉｉｉ）血管密度（ＶＬＤ）
ｉｖ）血管領域の体積（Ｖｏｌｕｍｅ）
のいずれかを選択する。これに限らず、例えば（２次元モーションコントラスト画像上で特定したり、特定した３次元血管領域を所定の深度範囲（例えば網膜外層）で投影したりして得られる）血管領域（例えば脈絡膜毛細血管）の面積を計測してもよい。
【０１７４】
なお、３次元画像処理による計測は大きく
１）３次元で強調し、２次元投影した該強調画像上で特定した血管領域もしくは血管中心線データに対する２次元計測
２）３次元で強調及び特定した血管領域もしくは血管中心線データを投影した場合の２次元計測
３）３次元で強調及び特定された血管領域もしくは血管中心線データに対する３次元計測
に大別できる。
【０１７５】
１）及び２）の例として、投影画像上で上記無血管領域の面積や血管密度、血管領域の面積や径、長さ、曲率を計測することが挙げられる。計測内容は正面モーションコントラスト画像に対する計測と同様であるものの、正面モーションコントラスト画像を強調・特定して計測する場合よりも血管抽出能が向上するため計測精度が向上する。
【０１７６】
また、３）の例として
３−１）血管の体積計測
３−２）任意方向の断面画像もしくは曲断面画像上の計測
（血管の径もしくは断面積計測も含む）
３−３）血管の長さや曲率計測
が挙げられる。
【０１７７】
本実施形態では、３次元で血管強調及び血管領域の特定処理を実行した後、網膜表層及び網膜深層の深度範囲で投影した各２値画像上でＶＡＤ、網膜外層の深度範囲で血管領域（脈絡膜新生血管領域）の体積を各々計測する。  なお、第一実施形態の場合と同様に、Ｍａｐボタン群から選択した計測の種類、Ｓｅｃｔｏｒボタン群から選択した計測の種類のいずれか一方が変更された場合に、他方も連動して（同じ計測の種類に）変更されるように構成してもよい。次に、解析部１０１−４６は、計測処理の前処理として画像拡大及びトップハットフィルタ処理を行う。本実施形態では、いずれも３次元のＢｉｃｕｂｉｃ補間及びトップハットフィルタ処理を実行する。
【０１７８】
＜ステップ３０７＞
解析部１０１−４６が血管領域の特定処理を行う。本実施形態では、強調部１０１−４６１が３次元ヘシアンフィルタ及び３次元エッジ選択鮮鋭化フィルタ処理に基づく血管強調処理を行う。次に抽出部１０１−４６２が第一実施形態の場合と同様に２種類の血管強調画像を用いて２値化処理を行い、整形処理を行うことで血管領域を特定する。血管領域特定処理の詳細はＳ５１０〜Ｓ５６０で説明する。
【０１７９】
＜ステップ３０８＞
計測部１０１−４６３が、操作者により指定された計測対象領域に関する情報に基づいて単検査の画像に対する計測を行う。引き続いて表示制御部１０１−０５が、計測結果を表示部１０４に表示する。なお、第一実施形態の場合と同様に、所定の解析に適した条件を満たさない状態で実施された計測結果を計測レポート画面に表示する際には表示部１０４に警告メッセージを表示するものとする。また、操作者が入力部１０３から血管領域もしくは血管中心線データを修正するよう指示を入力した場合には、第一実施形態の場合と同様に解析部１０１−４６が操作者から入力部１０３を介して指定された位置情報に基づいて血管領域もしくは血管中心線データを修正し、計測値を再計算する。網膜表層及び網膜深層におけるＶＡＤ計測、網膜外層における脈絡膜新生血管の体積計測についてはＳ８１０〜Ｓ８３０、網膜表層及び網膜深層におけるＶＬＤ計測、網膜外層における脈絡膜新生血管の総血管長計測についてはＳ８４０〜Ｓ８７０で各々説明する。
【０１８０】
＜ステップ３１０＞
比較部１０１−４６４は、第一実施形態の場合と同様の操作で、経時変化計測（Ｐｒｏｇｒｅｓｓｉｏｎ計測）処理を実行する。
【０１８１】
＜ステップ３１１＞
表示制御部１０１−０５は、Ｓ３１０で実施したＰｒｏｇｒｅｓｓｉｏｎ計測に関するレポートを表示部１０４に表示する。本実施形態では、Ｐｒｏｇｒｅｓｓｉｏｎ計測レポートの最上段に網膜表層において計測したＶＡＤマップを表示し、２段目に網膜深層において計測したＶＡＤマップ、３段目に網膜外層において計測した脈絡膜新生血管領域の
ｉ）２値画像（もしくは該２値画像の基準画像との差分画像）
ｉｉ）体積値もしくは総血管長（もしくは該２値画像の基準画像との差分値）
を時系列で並置表示する。これに限らず、例えば４段目に脈絡膜における血管密度（ＶＡＤもしくはＶＬＤ）マップを時系列で並置表示してもよい。これにより、眼底血管の３次元的な病態に関する時系列変化を一覧して把握できる。
【０１８２】
また、第一実施形態と同様に、各計測対象画像に関して略同一位置における断層画像数、ＯＣＴＡ重ね合わせ処理実施の有無、ＯＣＴＡ重ね合わせ処理の実施条件、ＯＣＴ断層画像もしくはモーションコントラスト画像の評価値の情報を表示部１０４に表示させてもよい。さらに、ショートカットメニューから計測の種類及び計測対象領域に関する項目を異なる値に変更することにより、Ｐｒｏｇｒｅｓｓｉｏｎ計測の種類及び計測対象領域を変更して再計測できる。第一実施形態と同様に、複数の計測対象領域が選択され、かつ一方の領域に対する計測の種類が変更された場合は他方の領域に対しても同じ計測の種類が連動して適用され、計測が実行される。さらに、第一実施形態と同様の方法で、所定の条件を満たさない状態で実施された計測結果を計測レポート画面に表示する際には警告メッセージを表示してよい。
【０１８３】
なお、本発明は、異なる深度範囲の正面画像及び該正面画像に対する計測値分布の時系列表示に限られるものではなく、例えば、該正面画像に直交する画像及び該正面画像に直交する画像に対する計測値分布、ボリュームレンダリングした３次元画像及び該３次元画像に対する計測値分布を時系列表示してもよい。
【０１８４】
さらに、図５（ａ）に示すフローチャートを参照しながら、Ｓ３０７で実行される処理の詳細について説明する。
【０１８５】
＜ステップ５１０＞
強調部１０１−４６１は、ステップ３０６の前処理を実施された３次元モーションコントラスト画像に対してヘッセ行列の固有値に基づく３次元血管強調フィルタ処理を行う。本実施形態では３次元のＭｕｌｔｉ−ｓｃａｌｅ  ｌｉｎｅ  ｆｉｌｔｅｒを用いるが、任意の公知の血管強調フィルタを用いてよい。なお、３次元のヘシアンフィルタでは３次元画像上の各画素において算出したヘッセ行列（（６）式）の「固有値（λ１、λ２、λ３）のうち１つが０に近く、残り２つが負かつ絶対値が大きい」場合に線状構造とみなして強調する。
【０１８６】
【数６】
2019154716.tif 000008
【０１８７】
３次元ヘシアンフィルタを用いると、深度方向に屈曲した血管に関しても
「血管走行方向の輝度変化が小さく、血管走行方向に直交する２方向の輝度が大きく低下する」という性質が成り立つため、良好に血管強調できるという利点がある。眼底血管の中には、深度方向に屈曲した血管として例えば
・脈絡膜側から網膜内に侵入する脈絡膜新生血管（ＣＮＶ）
・視神経乳頭部の血管
・網膜表層毛細血管と網膜深層毛細血管との接続部
が挙げられる。上記血管に対して正面モーションコントラスト画像上で２次元のヘシアンフィルタを適用すると、２次元平面内では「該平面内での血管走行方向の輝度変化が小さく、血管走行方向に直交する方向の輝度が大きく低下する」という性質が成り立たないため十分強調されず、血管領域として特定できないという課題がある。３次元ヘシアンフィルタを用いると上記血管についても良好に強調でき、血管検出能が向上する。
【０１８８】
＜ステップ５２０＞
抽出部１０１−４６２は、Ｓ５１０で生成した３次元ヘシアンフィルタによる血管強調画像（以下、３次元ヘシアン強調画像と表記）を２値化する。２値化の手順は第一実施形態の場合と同様であるが、３次元データの２値化であるという点が第一実施形態の場合と異なる。また、合成モーションコントラスト画像をヘシアンフィルタで強調処理した画像であるため、単独のモーションコントラスト画像をヘシアンフィルタで強調した場合に比べて２値化した血管領域の連続性が向上する。
【０１８９】
＜ステップ５３０＞
強調部１０１−４６１は、Ｓ３０６で生成したトップハットフィルタ適用後の合成モーションコントラスト画像に対して３次元エッジ選択鮮鋭化処理を行う。本実施形態では、３次元モーションコントラスト画像に対して３次元のＳｏｂｅｌフィルタを適用した画像を重みとして３次元のアンシャープマスク処理を行うことにより、エッジ選択鮮鋭化処理を実施する。
【０１９０】
＜ステップ５４０＞
抽出部１０１−４６２は、Ｓ５３０で生成したエッジ選択鮮鋭化処理を適用した鮮鋭化画像を２値化する。任意の公知の２値化法を用いてよいが、本実施形態では該３次元鮮鋭化画像上の各３次元局所領域内で算出した輝度統計値（平均値もしくは中央値）を閾値として２値化する。第一実施形態の場合と同様に、該閾値の上限及び下限値を設定することにより、血管領域内における抽出不足や無血管領域における誤抽出を抑制する。また、Ｓ５２０の場合と同様、合成モーションコントラスト画像をエッジ選択鮮鋭化した画像であるため、単独のモーションコントラスト画像をエッジ選択鮮鋭化した場合に比べて２値化した場合のノイズ状の誤検出領域が減少する。
【０１９１】
＜ステップ５５０＞
抽出部１０１−４６２は、Ｓ５２０で生成した３次元ヘシアン強調画像の２値画像の輝度値と、Ｓ５４０で生成した３次元エッジ選択強調画像の２値画像の輝度値の双方が０より大きい場合に血管候補領域として抽出する。当該演算処理により、ヘシアン強調画像に見られる血管径を過大評価している領域と、エッジ選択強調画像上に見られるノイズ領域がともに抑制され、血管の境界位置が正確かつ血管の連続性が良好な２値画像を取得できる。また、双方の２値画像とも合成モーションコントラスト画像に基づく２値画像であることから、単独のモーションコントラスト画像に基づく２値画像に比べて２値化した場合のノイズ状の誤検出領域が減少するとともに、特に毛細血管領域の連続性が向上する。また、合成モーションコントラスト画像であることから検査間の画質や輝度レベルも安定し、血管の抽出能が検査間で安定しやすい。
【０１９２】
＜ステップ５６０＞
抽出部１０１−４６２は、血管領域の整形処理として３次元のオープニング処理（収縮処理後に膨張処理を行うこと）及びクロージング処理（膨張処理後に収縮処理を行うこと）を実施する。なお、整形処理はこれに限らず２値画像をラベリングした場合の各ラベルの面積に基づく小領域除去を行ってもよい。
【０１９３】
なお、第一実施形態の場合と同様に、様々な径の血管が含まれるモーションコントラスト画像において血管強調する際のスケールを適応的に決定する方法はＳ５１０〜Ｓ５６０に述べた方法に限定されない。例えば図５（ｂ）のＳ６１０〜Ｓ６５０に示すように、３次元ヘシアン強調画像の輝度値と３次元エッジ選択鮮鋭化による血管強調画像の輝度値を乗じる演算を適用した画像に対する輝度統計値（例えば平均値）を閾値として２値化することによって血管領域を特定してもよい。該閾値には下限値や上限値を設定できる。あるいは、図５（ｃ）のＳ７１０〜Ｓ７４０に示すように、各画素の３次元位置（固視位置や深度範囲のデータでもよい）に基づいてヘシアンフィルタ適用時の平滑化フィルタのパラメータ（ガウスフィルタの平滑化パラメータσ）を適応的に変えた上でヘシアンフィルタを適用し、２値化することによって血管強調してもよい。また、２値化処理は閾値処理に限定されるものではなく任意の公知のセグメンテーション手法を用いることができる。
【０１９４】
さらに、図６（ａ）に示すフローチャートを参照しながら、Ｓ３０８で実行される処理の詳細について説明する。
【０１９５】
＜ステップ８１０＞
操作者は、入力部１０３を介して計測処理における関心領域を設定する。
【０１９６】
本実施形態では計測内容として
１）網膜表層及び網膜深層におけるＶＡＤマップ及びＶＡＤセクタマップ
２）網膜外層における脈絡膜新生血管の体積
を算出する。従って、関心領域として網膜表層及び網膜深層においては画像全体及び固視灯位置を中心とするセクタ領域を選択する。また網膜外層においては網膜外層に相当する層境界（ＯＰＬ／ＯＮＬ境界と、ブルッフ膜境界を該境界深層側に２０μｍ移動させた位置で囲まれる範囲）を指定する。
【０１９７】
＜ステップ８２０＞
計測部１０１−４６３は、Ｓ３０７で得られた血管領域の２値画像に基づいて計測処理を行う。網膜表層及び網膜深層における計測内容（ＶＡＤマップ及びＶＡＤセクタマップ生成）は基本的に第一実施形態と同様であるが、網膜表層及び網膜深層において特定された３次元血管領域を正面像として投影してから計測する点が異なる。網膜外層では、Ｓ８１０で設定した網膜外層に相当する関心領域内における非０画素（白画素）の体積を算出する。
【０１９８】
＜ステップ８３０＞
表示制御部１０１−０５は、計測結果としてＳ８２０で生成した網膜表層及び網膜深層におけるＶＡＤマップ及びＶＡＤセクタマップ、網膜外層における血管領域の２値画像と該血管領域の体積値を表示部１０４に表示する。なお、第一実施形態と同様の方法で、推奨する計測条件の少なくとも一つを満たさない状態で実施された計測結果を計測レポート画面に表示する際には（正確に計測することが期待できないような条件で計測が行われたとみなして）警告表示を行ってよい。
【０１９９】
上記では特定された３次元血管領域に基づいて体積を計測する場合の手順を例に説明したが、３次元の血管中心線に基づいて計測する場合には、上記Ｓ８１０〜８３０の代わりに図６（ｂ）に示すＳ８４０〜８７０を実行する。
【０２００】
＜ステップ８４０＞
計測部１０１−４６３は、Ｓ８２０で生成した血管領域の２値画像を３次元細線化処理することにより、血管の中心線に相当する線幅１画素のスケルトン画像を生成する。
【０２０１】
＜ステップ８５０＞
操作者は、入力部１０３を介してＳ８１０と場合と同様の関心領域を設定する。なお、本実施形態では計測内容として
１）網膜表層及び網膜深層におけるＶＬＤマップ及びＶＬＤセクタマップ
２）網膜外層における脈絡膜新生血管の総血管長
を算出する。
【０２０２】
＜ステップ８６０＞
計測部１０１−４６３はＳ８４０で得られたスケルトン画像に基づいて計測処理を行う。網膜表層及び網膜深層における計測内容（ＶＬＤマップ及びＶＬＤセクタマップ生成）は基本的に第一実施形態と同様であるが、網膜表層及び網膜深層において特定された３次元スケルトンを正面像として投影してから計測する点が異なる。なお、Ｓ３０７で特定した３次元血管領域を正面像として投影してから２次元細線化処理を実行して計測してもよい。網膜外層では、Ｓ８１０で設定した網膜外層に相当する関心領域内における非０画素（白画素）の長さの総和を算出する。
【０２０３】
＜ステップ８７０＞
表示制御部１０１−０５は、計測結果としてＳ８６０で生成した網膜表層及び網膜深層におけるＶＬＤマップ及びＶＬＤセクタマップ、網膜外層におけるスケルトン画像と該スケルトンの長さの総和を表示部１０４に表示する。また、Ｓ８３０の場合と同様に、所定の解析に適した条件を満たさない状態で実施された計測結果を計測レポート画面に表示する際には警告メッセージを表示するものとする。
【０２０４】
なお、本実施形態では脈絡膜新生血管を３次元抽出して体積や血管長の総和を時系列表示する場合の手順について説明したが、本発明はこれに限定されない。例えば、図１１の「Ｓｈｏｗ  Ｄｉｆｆｅｒｅｎｃｅ」チェックボックスを押下して基準画像における該新生血管の体積とその他の画像における該新生血管の体積との差分画像及び差分値を生成して時系列で表示してもよい。あるいは、視神経乳頭部の深度方向に屈曲した動静脈領域を本実施形態で説明した手順と同様にして特定し、血管径や血管断面積、血管中心線の曲率といった血管形状を計測してもよい。あるいは、網膜表層の毛細血管と網膜深層の毛細血管の連結部を３次元的に抽出して強調表示させたり、連結部の本数を計数したりしてもよい。
【０２０５】
以上述べた構成によれば、画像処理装置１０１は３次元の合成モーションコントラスト画像に対してモーションアーチファクト抑制処理を実施後、３次元のモルフォロジーフィルタ及び血管強調フィルタを適用し、２値化することにより３次元で血管領域を特定する。さらに特定した血管領域の体積を算出して血管領域の２値画像と体積値を時系列表示する。これにより、例えば、検査ごとのＯＣＴ断層画像の信号強度や画質のばらつきの影響を抑制しながら、血管病態の変化を正確に特定・計測できる。
【０２０６】
［その他の実施形態］
上記の各実施形態では、本発明を画像処理装置１０１として実現したが、本発明の実施形態は画像処理装置１０１のみに限定されるものではない。例えば、本発明はシステム、装置、方法、プログラムもしくは記憶媒体等としての実施態様をとることができる。
【図１】
2019154716.tif 000009
【図２】
2019154716.tif 000010
【図３】
2019154716.tif 000011
【図４】
2019154716.tif 000012
【図５】
2019154716.tif 000013
【図６】
2019154716.tif 000014
【図７】
2019154716.tif 000015
【図８】
2019154716.tif 000016
【図９】
2019154716.tif 000017
【図１０】
2019154716.tif 000018
【図１１】
2019154716.tif 000019
【図１２】
2019154716.tif 000020