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