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