(57)【要約】  入力画像から層のエッジを検出して検出されたエッジを示す境界線候補画像を作成する。入力画像の輝度値を微分して層の輝度勾配を示す輝度値微分画像を作成する。境界線候補画像と抽出しようとする境界線が存在する存在確率画像から得られる境界線位置確率画像と輝度値微分画像の最適な比率での重み付け計算から得られる評価得点画像を作成し、合計評価得点が最も高くなった経路を境界線として抽出する。このような画像処理装置、画像処理方法によれば、複数の層から構成される対象物体を撮影した画像から層の境界線を高精度で抽出することができる。
2017030056.tif 000002
【特許請求の範囲】
【請求項１】
複数の層から構成される対象物体を撮影した画像を入力画像とし、該入力画像から層のエッジを検出して検出されたエッジを層の境界線候補として示す境界線候補画像を作成する境界線候補画像作成手段と、
前記入力画像の輝度値を微分して層の輝度勾配を示す輝度値微分画像を作成する輝度値微分画像作成手段と、
前記境界線候補画像と抽出しようとする境界線が存在する存在確率画像から作成される境界線位置を示す境界線位置確率画像と、前記輝度値微分画像とに基づき経路探索を行って境界線を抽出する境界線抽出手段と、
を備えたことを特徴とする画像処理装置。
【請求項２】
前記入力画像を鉛直方向に（所望の画素数だけ）シフトした輝度値情報画像を作成する輝度値情報画像作成手段をさらに備えたことを特徴とする、請求項１に記載の画像処理装置。
【請求項３】
前記輝度値微分画像が、抽出しようとする境界線に応じた第１の重み付け係数で重み付けされることを特徴とする請求項１又は２に記載の画像処理装置。
【請求項４】
前記輝度値情報画像が、抽出しようとする境界線に応じた第２の重み付け係数で重み付けされることを特徴とする請求項２に記載の画像処理装置。
【請求項５】
前記輝度値微分画像が、抽出しようとする境界線に応じた第１の重み付け係数で重み付けされ、
前記輝度値情報画像が、抽出しようとする境界線に応じた第２の重み付け係数で重み付けされることを特徴とする請求項２に記載の画像処理装置。
【請求項６】
前記存在確率画像は、層の境界線毎に存在確率画像を格納した存在確率画像格納部から抽出しようとする境界線に応じて読み出されることを特徴とする請求項１から５のいずれか１項に記載の画像処理装置。
【請求項７】
前記第１の重み付け係数は、層の境界線毎に重み付け係数を格納した重み付け係数格納部から抽出しようとする境界線に応じて読み出されることを特徴とする請求項３又は５に記載の画像処理装置。
【請求項８】
前記第２の重み付け係数は、層の境界線毎に重み付け係数を格納した重み付け係数格納部から抽出しようとする境界線に応じて読み出されることを特徴とする請求項４又は５に記載の画像処理装置。
【請求項９】
前記境界線候補画像と抽出しようとする境界線が存在する存在確率画像から作成される境界線位置を示す境界線位置確率画像と、前記輝度値微分画像と、前記第１の重み付け係数とに基づき評価得点画像が作成され、前記境界線抽出手段は、前記作成された評価得点画像から評価得点の合計が最高になる経路を探索して、該合計評価得点が最高になった経路を境界線として抽出することを特徴とする請求項３又は５に記載の画像処理装置。
【請求項１０】
前記境界線候補画像と抽出しようとする境界線が存在する存在確率画像から作成される境界線位置を示す境界線位置確率画像と、前記輝度値微分画像と、前記第１の重み付け係数と、前記輝度値情報画像と、前記第２の重み付け係数とに基づき評価得点画像が作成され、前記境界線抽出手段は、前記作成された評価得点画像から評価得点の合計が最高になる経路を探索して、該合計評価得点が最高になった経路を境界線として抽出することを特徴とする請求項５に記載の画像処理装置。
【請求項１１】
抽出された境界線が修正された場合、前記存在確率画像格納部に格納されている当該境界線用の存在確率画像が修正されることを特徴とする請求項６に記載の画像処理装置。
【請求項１２】
前記第１の重み付け係数が修正された場合、前記重み付け係数格納部に格納されている当該第１の重み付け係数が修正されることを特徴とする請求項３又は５に記載の画像処理装置。
【請求項１３】
前記第２の重み付け係数が修正された場合、前記重み付け係数格納部に格納されている当該第２の重み付け係数が修正されることを特徴とする請求項４又は５に記載の画像処理装置。
【請求項１４】
前記入力画像が網膜の断層画像であることを特徴とする請求項１から１３のいずれか１項に記載の画像処理装置。
【請求項１５】
複数の層から構成される対象物体を撮影した画像を入力画像とし、該入力画像から層のエッジを検出して検出されたエッジを層の境界線候補として示す境界線候補画像を作成するステップと、
前記入力画像の輝度値を微分して層の輝度勾配を示す輝度値微分画像を作成するステップと、
前記境界線候補画像と抽出しようとする境界線が存在する存在確率画像から作成される境界線位置を示す境界線位置確率画像と、前記輝度値微分画像とに基づき経路探索を行って境界線を抽出するステップと、
を備えたことを特徴とする画像処理方法。
【請求項１６】
前記入力画像を鉛直方向に（所望の画素数だけ）シフトした輝度値情報画像を作成するステップをさらに備えたことを特徴とする請求項１５に記載の画像処理方法。
【請求項１７】
前記輝度値微分画像が、抽出しようとする境界線に応じた第１の重み付け係数で重み付けされることを特徴とする請求項１５又は１６に記載の画像処理方法。
【請求項１８】
前記輝度値情報画像が、抽出しようとする境界線に応じた第２の重み付け係数で重み付けされることを特徴とする請求項１６に記載の画像処理方法。
【請求項１９】
前記輝度値微分画像が、抽出しようとする境界線に応じた第１の重み付け係数で重み付けされ、
前記輝度値情報画像が、抽出しようとする境界線に応じた第２の重み付け係数で重み付けされることを特徴とする請求項１６に記載の画像処理方法。
【請求項２０】
前記境界線候補画像と抽出しようとする境界線が存在する存在確率画像から作成される境界線位置を示す境界線位置確率画像と、前記輝度値微分画像と、前記第１の重み付け係数とに基づき評価得点画像が作成され、該作成された評価得点画像から評価得点の合計が最高になる経路が探索され、該合計評価得点が最高になった経路が境界線として抽出されることを特徴とする請求項１７又は１９に記載の画像処理方法。
【請求項２１】
前記境界線候補画像と抽出しようとする境界線が存在する存在確率画像から作成される境界線位置を示す境界線位置確率画像と、前記輝度値微分画像と、前記第１の重み付け係数と、前記輝度値情報画像と、前記第２の重み付け係数とに基づき評価得点画像が作成され、該作成された評価得点画像から評価得点の合計が最高になる経路が探索され、該合計評価得点が最高になった経路が境界線として抽出されることを特徴とする請求項１９に記載の画像処理方法。
【請求項２２】
コンピュータを請求項１から１４のいずれか１項に記載の画像処理装置として機能させるための、あるいはコンピュータに請求項１５から２１のいずれか１項に記載の画像処理方法を実行させるための画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、断層像撮影装置などで撮影した被検眼などの対象物体の断層画像を処理して層の境界線を抽出するための画像処理装置、画像処理方法、及び画像処理プログラムに関するものである。
【背景技術】
【０００２】
眼科診断装置の一つで、眼底の断層像を撮影する光干渉断層計（ＯＣＴ：Optical Coherence Tomography）という光干渉を利用した断層像撮影装置がある。このような断層像撮影装置は広帯域な低コヒーレント光を、測定光として眼底に照射し、眼底からの反射光を参照光と干渉させて眼底の断層像を高感度に撮影することができる。【０００３】
このような断層像撮影装置では、網膜層内部の状態を３次元的に観察することができ、例えば神経線維層などの網膜層の層厚、あるいは網膜色素上皮層における凹凸などの層形状の変化を計測し、緑内障などの疾病の進行度や治療後の回復具合を定量的に診断することが可能になっている。【０００４】
特許文献１には、断層像撮影装置で撮影された断層像から網膜層の境界を検出し、眼底画像から病変の一つである白斑を抽出する構成が記載されている。【０００５】
特許文献２には、眼底の断層画像における偽像（アーティファクト）領域を特定し、偽像領域でない領域における網膜層の境界を検出するとともに、偽像領域における網膜層の境界を異なる方法で輝度値に基づいて検出し、それぞれ検出された境界を示す線を重畳して表示させる構成が記載されている。【０００６】
特許文献３には、被検眼の断層画像から得られる信号光の強度の低い側から高い側へのエッジに基づいて層を検出し、該信号光の強度の高い側から低い側へのエッジに基づいて該層の間に存在する層または層境界を検出する構成が記載されている。【０００７】
特許文献４には、ＭＲＩ画像の３次元空間における脳組織の存在確率をモデル化した存在確率モデルを事前に設定し、信号強度分布モデルと存在確率モデルが共に成立する組織分布モデルを求め、ＭＲＩ画像に含まれるボクセル毎に、ボクセルが白質組織及び灰白質組織に属する度合いを算出する構成が記載されている。【０００８】
非特許文献１には、キャニーエッジ検出法を用いて取得される網膜断層像のエッジ画像と、輝度勾配画像をそれぞれ重み付けして演算し、最短経路を探索して網膜層の境界線を抽出する構成が記載されている。また、同非特許文献１には、複数の網膜層の境界線を抽出するとき、最初に２つの境界線を抽出し、その間に存在する他の境界線を、抽出した両境界線内の狭い範囲内で探索して抽出時間を短縮する構成が記載されている。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】特開２０１０−２７９４３８号公報
【特許文献２】特開２０１２−６１３３７号公報
【特許文献３】特許第５６６５７６８号公報
【特許文献４】特開２０１１−３０９１１号公報
【非特許文献】
【００１０】
【非特許文献１】「Automated layer segmentation of macular OCT images using dual-scale gradient information」27 September 2010/ Vol. 18, No, 20/ OPTICS EXPRESS 21293-21307
【発明の概要】
【発明が解決しようとする課題】
【００１１】
しかしながら、特許文献１〜３の構成では、断層画像のエッジを検出することにより層の境界線を抽出しているので、境界線の抽出精度がエッジ検出に依存し、輝度値の低い領域ではエッジ情報を十分な精度で得られず、信頼性のある境界線抽出を行うのが困難である、という問題がある。【００１２】
特許文献２の構成では、輝度値の低い偽像領域内では、モデル形状に関する評価関数値によって層の位置を求めている。しかし、この評価関数値の精度が不十分で、輝度値の低い領域での境界線抽出精度が劣化するという問題がある。【００１３】
非特許文献１では、網膜断層像のエッジ画像と、輝度勾配画像をそれぞれ重み付けして、境界線を抽出している。しかし、重み付けが一次元的で、特にエッジ画像と、抽出しようとする境界線の存在する確率との関連付けがなく、境界線を精度よく抽出ができない、という欠点がある。【００１４】
特許文献４では、脳組織の存在確率を示す存在確率モデルを用いて、脳組織の領域が、白質組織及び灰白質組織に属する度合いを算出しているが、存在確率モデルは三次元的なモデルで、算出に時間を要するという問題がある。【００１５】
また、非特許文献１では、抽出した２つの境界線間に存在する他の境界線を、抽出した両境界線内の狭い範囲内で探索して抽出時間を短縮している。しかし、非特許文献１では、抽出した結果を基に抽出領域を制限し、更にその抽出結果を基に抽出領域を制限するという処理を繰り返して逐次的に行っていないので、抽出された境界線が、既に抽出された他の境界線と交わってしまう、あるいは不明瞭な境界線ないし消失した境界線を精度よく抽出することができない、という問題がある。【００１６】
本発明は、このような点に鑑みてなされたもので、複数の層から構成される対象物体を撮影した画像から層の境界線を高精度で抽出することが可能な画像処理装置、画像処理方法、及び画像処理プログラムを提供することを課題とする。
【課題を解決するための手段】
【００１７】
本発明は、
複数の層から構成される対象物体を撮影した画像を入力画像とし、該入力画像から層のエッジを検出して検出されたエッジを層の境界線候補として示す境界線候補画像を作成すること、
前記入力画像の輝度値を微分して層の輝度勾配を示す輝度値微分画像を作成すること、
前記境界線候補画像と抽出しようとする境界線が存在する存在確率画像から作成される境界線位置を示す境界線位置確率画像と、前記輝度値微分画像とに基づき経路探索を行って境界線を抽出することを特徴とする。【００１８】
輝度値微分画像は、抽出しようとする境界線に応じた重み付け係数で重み付けされる。
存在確率画像は層の境界線毎に存在確率画像を格納した存在確率画像格納部から抽出しようとする境界線に応じて読み出され、また重み付け係数は層の境界線毎に重み付け係数を格納した重み付け係数格納部から抽出しようとする境界線に応じて読み出される。
前記境界線候補画像と抽出しようとする境界線が存在する存在確率画像から作成される境界線位置確率画像と、前記輝度値微分画像とに基づき評価得点画像が作成され、境界線抽出手段は、前記作成された評価得点画像から評価得点の合計が最高になる経路を探索して、該合計評価得点が最高になった経路を境界線として抽出する。
抽出された境界線が修正された場合、格納部に格納されている当該境界線用の存在確率画像が修正され、また重み付け係数が修正された場合には、格納部に格納されている当該境界線用の重み付け係数が修正される。
【発明の効果】
【００１９】
本発明では、存在確率画像で入力画像の境界線位置情報を求め、輝度値微分画像により入力画像の輝度値情報を求め、この境界線位置情報と輝度値情報を組み合わせることにより境界線を抽出しているので、高精度で境界線を抽出することができる。
また、本発明では、抽出しようとする境界線が存在している部分は、存在確率画像による演算により評価得点が高くなるので、境界線の抽出精度を顕著に向上させることができる。
また、抽出しようとする境界線に応じて輝度情報を重み付けすることにより、特性の異なる境界線を高精度で抽出することが可能となる。
【図面の簡単な説明】
【００２０】
【図１】画像処理装置の全体構成を示したブロック図である。
【図２】眼底をスキャンすることにより眼底網膜の断層画像を取得する状態を示した説明図である。
【図３】取得された断層画像の網膜層ないしその境界線を示した説明図である。
【図４】網膜層の境界線を抽出する工程を示したフローチャートである。
【図５】入力画像に基づいて読み出された存在確率画像、並びに入力画画像から作成された境界線候補画像と輝度値微分画像を示した説明図である。
【図６】入力画像から境界線を抽出して結果画像を取得する過程を示した説明図である。
【図７】存在確率画像格納部の構成を示すブロック図である。
【図８】境界線を抽出するための経路探索を示した説明図である。
【図９】複数の境界線を抽出する過程を示したフローチャートである。
【図１０】（ａ）は探索範囲を固定的に設定して境界線を抽出する方法を示した説明図、（ｂ）は探索範囲を動的に設定して境界線を抽出する方法を示した説明図である。
【図１１】不明瞭ないし消失した境界線を抽出する方法を示した説明図である。
【図１２】抽出した境界線上に画素間隔を変えて制御点を表示させる状態を示した説明図である。
【図１３】抽出した境界線上に画素間隔を変えて制御点を表示させ、制御点を移動させて境界線を修正する状態を示した説明図である。
【図１４】制御点を移動した後に再度経路探索を行って他の境界線を抽出する方法を示した説明図（その１）である。
【図１５】制御点を移動した後に再度経路探索を行って他の境界線を抽出する方法を示した説明図（その２）である。
【発明を実施するための形態】
【００２１】
以下、図面を参照しながら実施例に基づいて本発明を詳細に説明する。ここでは、処理の対象物体の画像は、被検眼眼底の断層画像を例にして説明するが、本発明における処理の対象物体となる画像は、眼底の断層画像に限定されるものでなく、複数の層が撮影される対象物体の画像にも適用されるものである。
【実施例１】
【００２２】
＜全体構成＞
図１は、被検眼眼底の断層画像を取得して画像処理するシステム全体を示すブロック図である。断層像撮影装置１０は、被検眼の眼底の断層像を撮影する装置（ＯＣＴ：Optical Coherence Tomography）であり、例えばフーリエドメイン方式で動作する。断層像撮影装置１０は公知であるので、その詳細な説明は省略するが、断層像撮影装置１０には、低コヒーレンス光源が設けられ、低コヒーレンス光源からの光が参照光と信号光に分割される。信号光は、図２に図示したように、眼底Ｅ上で、例えばＸ、Ｙ方向にラスタースキャンされる。眼底Ｅで走査され反射された信号光は、参照ミラーで反射した参照光と重畳され干渉光を発生し、該干渉光に基づいて眼底の深度方向（Ｚ方向）の情報を示すＯＣＴ信号が発生する。【００２３】
画像処理装置２０は、ＣＰＵ、ＲＡＭ、ＲＯＭなどで構成されたコンピュータで実現される制御部２１を有し、制御部２１は画像処理プログラムを実行することにより全体の画像処理を制御する。また、画像処理装置２０には、断層画像形成部２２が設けられる。【００２４】
断層画像形成部２２は、フーリエドメイン方式などの公知の解析方法を実行する専用の電子回路、または、前述のＣＰＵが実行する画像処理プログラムにより実現され、断層像撮影装置１０で生成されたＯＣＴ信号に基づいて、眼底の断層画像を形成する。【００２５】
例えば、図２に示したように、眼底ＥがＹ方向にｙＮ（Ｎ＝１、２、．．．．．ｎ）の位置で、Ｘ方向にスキャンされた場合、当該スキャン時に複数回（ｍ回）サンプリングが行われる。このＸ方向の各サンプリング時点でそれぞれＺ方向の断層画像（Ａスキャン画像）Ａｈ（ｈ＝１、２、．．．．．ｍ）が取得され、これらのＡスキャン画像Ａｈから断層画像ＢＮ（Ｎ＝１、２、．．．．．ｔ）が形成される。Ａスキャン画像は、例えば、Ｘ方向に１画素幅、Ｚ方向にｎ画素の長さで格納されるので、断層画像ＢＮはｍ×ｎ画素の大きさをもつ画像となり、Ｂスキャン画像とも呼ばれる。【００２６】
断層画像形成部２２で形成されたｔ枚の断層画像ＢＮ、あるいはこれらｔ枚の断層画像ＢＮから構築される３次元のボリューム画像は、例えば半導体メモリ、ハードディスク装置等により構成された記憶部２３に格納される。記憶部２３は、さらに上述した画像処理プログラムなども格納する。【００２７】
画像処理装置２０には、画像処理部３０が設けられる。画像処理部３０は、境界線候補画像作成手段３１、輝度値微分画像作成手段３２、輝度値情報画像作成手段３２ａ、評価得点画像作成手段３３、境界線抽出手段３４、探索範囲設定手段３５、並びに制御点設定手段３６を備えている。【００２８】
後述するように、境界線候補画像作成手段３１は入力画像のエッジを検出して境界線候補画像を作成し、輝度値微分画像作成手段３２は入力画像の輝度値を微分して輝度勾配を示す輝度値微分画像を作成し、輝度値情報画像作成手段３２ａは入力画像を鉛直方向にシフトして輝度情報を示す輝度値情報画像を作成し、評価得点画像作成手段３３は作成された画像並びに読み出された画像に基づいて境界線抽出のための評価得点を示す評価得点画像を作成し、境界線抽出手段３４は評価得点画像から評価得点の合計値が最高になる経路を探索してそれを境界線として抽出する。また、探索範囲設定手段３５は、探索範囲を抽出済みの境界線に合わせて設定し、制御点設定手段３６は抽出された境界線上に所定の画素間隔で制御点を設定する。画像処理部３０における各手段あるいは各画像処理は、専用の電子回路を用いることにより、あるいは画像処理プログラムを実行することに実現される。【００２９】
存在確率画像格納部２６は、後述するように、抽出しようとする境界線毎に当該境界線が存在する存在確率を示す画像を格納する。【００３０】
重み付け係数格納部２７は、抽出しようとする境界線毎に輝度値微分画像に重み付けされる重み付け係数及び輝度値情報画像に重み付けされる重み付け係数を格納する。【００３１】
表示部２４は、例えば、ＬＣＤなどのディスプレイ装置によって構成され、記憶部２３に格納された断層画像、画像処理装置２０で生成あるいは処理された画像、制御点設定手段３６で設定された制御点、被験者に関する情報などの付随する情報などを表示する。【００３２】
操作部２５は、例えば、マウスやキーボード、操作ペン、ポインター、操作パネル等を有し、表示部２４に表示された画像の選択、あるいは操作者が画像処理装置２０などに指示を与えるために用いられる。【００３３】
このような構成で撮影された断層画像のうち、図２に示す眼底Ｅの黄斑部Ｒを通過する走査線ｙｋで取得される断層画像Ｂｋが図３の上部に符号Ｂを付して図示されている。眼底の網膜は、種々の膜あるいは層の組織から構成されており、断層画像Ｂにおいて判別できる膜あるいは層が同図に示されている。【００３４】
図３において、ＩＬＭは内境界膜（internal limiting membrane）を示し、ＮＦＬは神経線維層（nerve fiber layer）を、ＧＣＬは神経節細胞層（ganglion cell layer）を、ＩＰＬは内網状層（inner plexiform layer）を、ＩＮＬは内顆粒層 (inner nuclear layer)を、ＯＰＬは外網状層 (outer plexiform layer)を、ＯＮＬは外顆粒層（outer nuclear layer)を、ＥＬＭは外境界膜 (external limiting membrane)を、ＩＳは視細胞層内節 (inner photoreceptor segments)を、ＯＳは視細胞層外節 (outer photoreceptor segments)を、ＲＰＥは網膜色素上皮 (retinal pigment epithelium)を、それぞれ示す。【００３５】
図３の下方には、これらの膜ないし層間の境界線がそれぞれＬ１〜Ｌ１０でかっこを付して示されている。例えば、Ｌ１は内境界膜ＩＬＭにより形成される境界線であり、Ｌ２は神経線維ＮＦＬと神経節細胞層ＧＣＬ間の境界線である。これらの境界線Ｌ１〜Ｌ１０が本発明で抽出しようとする境界線である。【００３６】
このような断層画像から、神経線維層などの層厚、網膜色素上皮層における凹凸などの層形状の変化を計測できれば、眼科疾患の進行度や治療後の回復具合を定量的に診断することが可能となる。しかし、撮影環境によっては、ＯＣＴ信号が減衰もしくは欠損して各層の境界線が不明瞭になったり、あるいは不連続ないしは消失してしまい正確な層厚あるいは層形状の計測が困難になる。【００３７】
そこで、本発明では、以下に説明する方法で網膜層の境界線を高精度で抽出する。以下に、その方法を図４のフローチャートを用いて説明する。【００３８】
＜境界線抽出処理＞
まず、図４のステップＳ１に示したように、網膜層の境界線を抽出しようとする断層画像ＢＮ（Ｎ＝１、２、．．．．．ｔ）を記憶部２３から読み出して表示部２４に表示させ入力画像Ｂを選択する。この入力画像は、図２に示したｔ枚の断層画像ＢＮのうち一枚、あるいは数枚、あるいは全部でもよく、または任意の画像処理を加えた後の画像であってもよい。複数枚の断層画像の境界線を抽出する場合には、いずれか一枚の断層画像を入力画像として選択する。ここでは、例えば、図２の眼底の黄斑部Ｒをスキャンする走査線ｙｋで取得される断層画像Ｂｋを記憶部２３から読み出し入力画像Ｂとする。【００３９】
入力画像Ｂが選択されたら、この入力画像Ｂから抽出しようとする境界線を定め、その境界線用の存在確率画像を存在確率画像格納部２６から読み出す（ステップＳ２）。【００４０】
この存在確率画像格納部２６が図７に図示されている。存在確率画像格納部２６には、断層画像Ｂｋに対して断層画像Ｂｋの層の境界線Ｌ１〜Ｌ１０毎に、当該境界線が存在する存在確率を示す画像が格納されている。例えば、左上のＩＬＭ（Ｌ１）は内境界膜（ＩＬＭ）の境界線Ｌ１が存在する確率を示す画像であり、その下のＮＦＬ／ＧＣＬ（Ｌ２）は神経線維層（ＮＦＬ）と神経節細胞層（ＧＣＬ）との境界を示す２番目の境界線Ｌ２が存在する確率を示す画像である。以下、その他の存在確率画像も、同様に格納されている。【００４１】
このような存在確率画像は、複数の正常眼に対して、同じ走査線ｙｋで取得される断層画像Ｂｋから境界線が存在する確率を各境界線に対して断層画像の各画素（ｉ、ｊ）［ｉ＝１、２、．．．．．ｍ、ｊ＝１、２、．．．．．ｎ］毎に事前に計算し、断層画像Ｂｋの各画素（ｉ、ｊ）の対応する画素位置にデジタル値として格納したｍ×ｎ画素からなる画像である。【００４２】
例えば、断層画像Ｂｋの境界線Ｌ１の存在確率画像ＩＬＭ（Ｌ１）が図７の下方にＨＩＬＭとして模式化されて図示されており、境界線Ｌ１が存在する確率領域（画素領域）がパーセントで示されている。断層画像Ｂｋの画素（ｉ、ｊ）が決まれば、その位置に当該境界線Ｌ１が存在する確率を、パーセントに応じたデジタル値としてＨＩＬＭから求めることができる。図７には、２５％置きの確率で区分した画素領域が図示されている。何％置きの確率で区分するかは求める精度に応じて決められるが、実際には、図７に示したよりも、更に細かい確率で区分されている。【００４３】
なお、図７には、存在確率画像格納部２６には、断層画像Ｂｋの各層毎の１０個の存在確率画像が格納されていることが図示されている。しかし、ボリューム画像を構成する断層画像Ｂｋ以外の断層画像ＢＮに対しても、それぞれの断層画像毎に、各層用の存在確率画像が格納されても良いし、全画像Ｂ１〜ＢＮに共通の一枚の存在確率画像が格納されるものとしても良い。さらに、存在確率画像は、Ｘ方向は均一でＺ方向は連続的に単調増加あるいは単調減少する画像としてもよい。また、存在確率画像格納部２６に格納される各存在確率画像はマップ化されており、抽出された境界線が修正された場合には、抽出に寄与した存在確率画像も適応的に修正することができ、逐次学習させることができる。【００４４】
また、存在確率画像格納部２６は、画像処理装置２０内に設けるのではなく、外部の格納部とすることができ、例えば、インターネットを介して接続されたサーバーに設けるようにしてもよい。【００４５】
ここでは、断層画像Ｂｋの境界線のうち、最初に内境界膜（ＩＬＭ）の境界線Ｌ１を抽出する例を説明するので、入力画像Ｂ、並びにその境界線Ｌ１で決まる存在確率画像ＩＬＭ（Ｌ１）が存在確率画像格納部２６から読み出され、この読み出された存在確率画像が図５の上部にＨＩＬＭとして図示されている。【００４６】
続いて、画像処理部３０の境界線候補画像作成手段３１により入力画像Ｂのエッジを検出して境界線候補画像を作成する（ステップＳ３）。このエッジ検出には、例えば、公知のキャニーエッジ検出法（Canny Edge Detection）を用いる。このキャニーエッジ検出法で抽出されたエッジは細いエッジとなり、コントラストの高い領域には高いしきい値を設定するなどしきい値を適正に設定すると、入力画像Ｂの境界線候補となる複数の細い線からなる境界線候補画像を作成することができる。この境界線候補画像は、入力画像Ｂのエッジの有無の情報を示す値をデジタル値としてその各画素（ｉ、ｊ）に対応した画素位置に有するｍ×ｎ画素の画像であり、図５にＥＩＬＭとして図示されている。【００４７】
続いて、輝度値微分画像作成手段３２により入力画像Ｂの輝度値をＺ方向に微分して輝度値微分画像を作成する（ステップＳ４）。輝度値を微分することはＺ方向の輝度勾配を演算することになる。この輝度勾配は、例えば、公知のゾーベルフィルタ(Sobel Filter）などの差分型フィルタを用いて演算する。輝度値微分画像は、入力画像Ｂの輝度勾配を示す値をデジタル値としてその各画素（ｉ、ｊ）に対応した画素位置に有するｍ×ｎ画素の画像であり、図５にＧＩＬＭとして図示されている。輝度値を微分すると、網膜層の輝度勾配を検出することができるので、輝度値微分画像は境界線候補画像のエッジ情報が欠損していたり不足している場合に、その情報を補足することができる。【００４８】
続いて、輝度値微分画像用に設定された求める境界線ＩＬＭ用の重み付け係数ＷＩＬＭを重み付け係数格納部２７から読み出す（ステップＳ４’）。【００４９】
さらに、入力画像Ｂを鉛直方向に所望の画素数だけシフトした画像を輝度値情報画像Ｂ’として作成し（ステップＳ５）、当該輝度値情報画像Ｂ’用に設定されたＩＬＭ用の重み付け係数ＱＩＬＭを重み付け係数格納部２７から読み出す（ステップＳ５’）。鉛直方向の移動量や移動方向は求めたい境界線により適宜変更することができるが、例としてＩＬＭの場合は上方向に５画素移動する。すなわち、求めたい境界線の上に輝度情報が重なるように移動方向やシフト量を決定するとよい。【００５０】
続いて、評価得点画像作成手段３３により境界線候補画像ＥＩＬＭ、存在確率画像ＨＩＬＭ、輝度値微分画像ＧＩＬＭ、輝度値微分画像用に設定された重み付け係数ＷＩＬＭ、輝度値情報画像Ｂ’、輝度値情報画像用に設定された重み付け係数ＱＩＬＭに基づいて評価得点画像ＣＩＬＭを、以下の式に基づいて演算し作成する（ステップＳ６）。
ＣＩＬＭ＝ＥＩＬＭ×ＨＩＬＭ＋ＷＩＬＭ×ＧＩＬＭ＋ＱＩＬＭ×Ｂ’【００５１】
この評価得点画像ＣＩＬＭの作成過程が図６に図示されている。評価得点画像ＣＩＬＭは、境界線候補画像ＥＩＬＭと存在確率画像ＨＩＬＭとの演算（掛け算）により得られる境界線位置とその存在確率を含んだ２次元的情報を示す境界線位置確率画像と、重み付け係数ＷＩＬＭで重み付けされた輝度値微分画像ＧＩＬＭと、重み付け係数ＱＩＬＭで重み付けされた輝度値情報画像Ｂ’との演算（加算）により形成される。なお、抽出しようとする境界線によって、エッジ要素、輝度値微分要素、輝度値情報のどれを重視することにより抽出しやすくなるかは異なる。この実態に対応すべく、重み付け係数ＷＩＬＭ、ＱＩＬＭを任意の値に設定することにより、エッジを重視して抽出するか、輝度値微分を重視して抽出するか、輝度値情報を重視して抽出するか、いずれの場合でもより良く抽出できる。評価得点画像ＣＩＬＭは、境界線候補画像ＥＩＬＭ、存在確率画像ＨＩＬＭ、並びに輝度値微分画像ＧＩＬＭの各画素（ｉ、ｊ）毎に演算されるので、評価得点画像ＣＩＬＭは各画素（ｉ、ｊ）に対応した画素位置に上述した式による演算値を有するｍ×ｎ画素からなる画像である。図６における各画像は実際には複雑な形状であるが、図示が困難であるので、その概略が模式的に図示されている。【００５２】
評価得点画像ＣＩＬＭは各画素（ｉ、ｊ）がデジタル値でスコア化されているので、例えば動的計画法（Dynamic Programming）を用いて経路探索し、合計得点が最高になる経路を探索してＩＬＭの境界線を抽出する。【００５３】
この経路探索例が図８に図示されている。例えば、探索範囲を上下に示した点線の範囲に設定する。ｉ＝１の画素位置にあるＺ方向に延びる画素列を１番目の画素列Ｐ１とする。画素列Ｐ２の各画素に対して、Ｐ１の近傍の画素中で最も評価得点が高い画素を探索し、その画素の位置を記憶するとともに、その画素の評価得点とＰ２の各画素の評価得点を加算した値を、Ｐ２の各画素の評価得点として更新する。続いてｉを増分し、各画素列Ｐ３、．．．．．Ｐｍまで順次繰り返す。その後、画素列Ｐｍで評価得点が最も大きい画素を始点とし、ｉを減らし記憶した画素の位置を順に辿る事で、総和が最も高くなるような太い点線で示した経路を境界線として抽出する。【００５４】
図８において、細い湾曲した線は、図４のステップＳ３の処理で検出されたエッジＥＩＬＭと抽出しようとする境界線Ｌの存在確率ＨＩＬＭとの演算により得られる線であり、また幅のある横長の画像は、図４のステップＳ４の処理で得られた輝度値微分画像ＧＩＬＭである。輝度勾配がある領域では、評価得点が高くなり、抽出しようとする境界線Ｌがこの領域を通過することが理解できる。必要な場合には、始点を変化させてすべての経路について総和を求め、そのうち最高得点となった経路を内境界膜（ＩＬＭ）の境界線Ｌ１として抽出する（ステップＳ７）。【００５５】
なお、図８では、経路探索は画素列Ｐ１から開始されているが、最終画素列ｉ＝Ｐｍから開始するようにしてもよい。【００５６】
抽出された境界線Ｌ１は、図６の下段に示したように、入力画像Ｂに重ねられ、結果画像ＲＩＬＭとして表示部２４に表示される（ステップＳ８）。【００５７】
続いて、全ての境界線を抽出したかどうかの判断が行われ（ステップＳ９）、抽出していない境界線がある場合には、ステップＳ２に戻り、全ての境界線を抽出した場合には、処理を終了する。【００５８】
以上説明したように、本実施形態では、境界線候補画像と存在確率画像で入力画像の位置情報を求め、輝度値微分画像により入力画像の輝度値情報を求め、この位置情報と任意の重み付けをした輝度値情報を組み合わせることにより評価得点画像を形成している。この評価得点画像の各画素のうち抽出しようとする境界線が存在している画素は、存在確率画像による演算により評価得点が高くなっており、この評価得点の高い画素を探索して境界線を決定しているので、境界線の抽出精度を顕著に向上させることができる。【００５９】
なお、輝度値微分画像に掛る重み付け係数は、抽出しようとする境界線によっては省略することができる。複数の境界線を抽出する場合には、抽出しようとする境界線に応じて輝度勾配情報を重み付けすることにより、特性の異なる境界線を高精度で抽出することが可能となる。【００６０】
また、上述したステップＳ７において抽出した境界線を入力画像に重ねて表示したとき抽出した境界線が本来の境界線とずれている場合には、後述するように、ユーザーが境界線位置を修正することができる。この場合には、修正に応じて存在確率画像格納部２６に格納されている当該境界線用の存在確率画像を修正して学習させることができる。【００６１】
また、境界線抽出過程で輝度値微分画像に掛る重み付け係数を修正することでより良好な境界線抽出が可能になる。このような場合には、重み付け係数の修正に応じて重み付け係数格納部２７に格納されている当該境界線用の重み付け係数も修正して学習させる。【００６２】
＜複数の境界線抽出＞
図４のステップＳ９で示したように、入力画像から複数の境界線を抽出する場合には、抽出済みの境界線を利用して他の境界線を効果的に抽出することができる。以下に、その実施形態を図９を参照して説明する。【００６３】
境界線のうち最上端の内境界膜ＩＬＭ（Ｌ１）と最下端の網膜色素上皮ＲＰＥ（Ｌ１０）は輝度変化が大きな境界であるので、抽出が容易であることからこれらの境界線を最初に抽出し、その抽出済みの境界線を利用して探索範囲を限定ないし設定して、他の境界線を抽出する。【００６４】
図９において、入力画像Ｂの各層に対して、図４のステップＳ３で説明したようなエッジ検出処理を行い、各層用の境界線候補画像ＥＩＬＭ、ＥＮＦＬ／ＧＣＬ，ＥＧＣＬ／ＩＰＬ、．．．．．を作成する（ステップＴ１）。【００６５】
続いて、入力画像Ｂの各層に対して輝度値微分を行い輝度勾配を抽出して（図４のステップＳ４の処理）、各層用の輝度値微分画像ＧＩＬＭ、ＧＮＦＬ／ＧＣＬ，ＧＧＣＬ／ＩＰＬ、．．．．．を作成する（ステップＴ２）。【００６６】
このような処理を行った後、存在確率画像格納部２６、重み付け係数格納部２７から各層用の存在確率画像ＨＩＬＭ、ＨＮＦＬ／ＧＣＬ，ＨＧＣＬ／ＩＰＬ、．．．．．、輝度値微分画像用に設定された各層用の重み付け係数ＷＩＬＭ、ＷＮＦＬ／ＧＣＬ，ＷＧＣＬ／ＩＰＬ、．．．．．、並びに輝度値情報用に設定された各層用の重み付け係数ＱＩＬＭ，ＱＮＦＬ／ＧＣＬ，ＱＧＣＬ／ＩＰＬ、．．．．．、をそれぞれ読み出し、図４のステップＳ６で説明したような処理を行って各層用の評価得点画像を作成する。【００６７】
まず、最上端のＩＬＭは輝度変化が大きな境界であるので、抽出する順序としてＩＬＭを選択する。探索範囲を入力画像全体として評価得点ＥＩＬＭ×ＨＩＬＭ＋ＷＩＬＭ×ＧＩＬＭ＋ＱＩＬＭ×Ｂ’をパラメータとする合計得点が最高になる経路を探索し、該探索された経路をＩＬＭの境界線Ｌ１として抽出する（ステップＴ３）。このＩＬＭに対する境界線Ｌ１を抽出する処理は、図５、図６、図８に関連して述べた処理に対応する。【００６８】
次に、最下端のＲＰＥを選択し、同様に探索範囲を入力画像全体としてパラメータＥＲＰＥ×ＨＲＰＥ＋ＷＲＰＥ×ＧＲＰＥ＋ＱＲＰＥ×Ｂ’の合計得点が最高になる経路を探索する。探索の結果、合計得点が最高になった経路をＲＰＥの境界線Ｌ１０として抽出する（ステップＴ４）。なお、最初にＲＰＥの境界線Ｌ１０を抽出し、続いてＩＬＭの境界線Ｌ１を抽出するようにしてもよい。【００６９】
続いて、探索範囲を抽出済みの境界線ＩＬＭ（Ｌ１）〜ＲＰＥ（Ｌ１０）の範囲に設定して、パラメータＥＩＳ／ＯＳ×ＨＩＳ／ＯＳ＋ＷＩＳ／ＯＳ×ＧＩＳ／ＯＳ＋ＱＩＳ／ＯＳ×Ｂ’の合計得点が最高になる経路を探索し、該探索された経路をＩＳ／ＯＳの境界線Ｌ８として抽出する（ステップＴ５）。【００７０】
続いて、探索範囲を抽出済みの境界線ＩＬＭ（Ｌ１）〜ＩＳ／ＯＳ（Ｌ８）の範囲に設定して、パラメータＥＯＰＬ／ＯＮＬ×ＨＯＰＬ／ＯＮＬ＋ＷＯＰＬ／ＯＮＬ×ＧＯＰＬ／ＯＮＬ＋ＱＯＰＬ／ＯＮＬ×Ｂ’の合計得点が最高になる経路を探索し、該探索された経路をＯＰＬ／ＯＮＬの境界線Ｌ６として抽出し（ステップＴ６）、また、探索範囲を抽出済みの境界線ＩＳ／ＯＳ（Ｌ８）〜ＲＰＥ（Ｌ１０）の範囲に設定して、パラメータＥＯＳ／ＲＰＥ×ＨＯＳ／ＲＰＥ＋ＷＯＳ／ＲＰＥ×ＧＯＳ／ＲＰＥ＋ＱＯＳ／ＲＰＥ×Ｂ’の合計得点が最高になる経路を探索し、該探索された経路をＯＳ／ＲＰＥの境界線Ｌ９として抽出する（ステップＴ７）。【００７１】
同様に、探索範囲を抽出済みの境界線ＩＬＭ（Ｌ１）〜ＯＰＬ／ＯＮＬ（Ｌ６）の範囲に設定して、パラメータＥＮＦＬ／ＧＣＬ×ＨＮＦＬ／ＧＣＬ＋ＷＮＦＬ／ＧＣＬ×ＧＮＦＬ／ＧＣＬ＋ＱＮＦＬ／ＧＣＬ×Ｂ’の合計得点が最高になる経路を探索し、該探索された経路をＮＦＬ／ＧＣＬの境界線Ｌ２として抽出し（ステップＴ８）、また、探索範囲を抽出済みの境界線ＯＰＬ／ＯＮＬ（Ｌ６）〜ＩＳ／ＯＳ（Ｌ８）の範囲に設定して、パラメータＥＥＬＭ×ＨＥＬＭ＋ＷＥＬＭ×ＧＥＬＭ＋ＱＥＬＭ×Ｂ’の合計得点が最高になる経路を探索し、該探索された経路をＥＬＭの境界線Ｌ７として抽出する（ステップＴ９）。【００７２】
以下同様に、探索範囲を抽出済みの境界線ＮＦＬ／ＧＣＬ（Ｌ２）〜ＯＰＬ／ＯＮＬ（Ｌ６）の範囲に設定して、パラメータＥＩＰＬ／ＩＮＬ×ＨＩＰＬ／ＩＮＬ＋ＷＩＰＬ／ＩＮＬ×ＧＩＰＬ／ＩＮＬ＋ＱＩＰＬ／ＩＮＬ×Ｂ’の合計得点が最高になる経路を探索し、該探索された経路をＩＰＬ／ＩＮＬの境界線Ｌ４として抽出し（ステップＴ１０）、また、探索範囲を抽出済みの境界線ＮＦＬ／ＧＣＬ（Ｌ２）〜ＩＰＬ／ＩＮＬ（Ｌ４）の範囲に設定して、パラメータＥＧＣＬ／ＩＰＬ×ＨＧＣＬ／ＩＰＬ＋ＷＧＣＬ／ＩＰＬ×ＧＧＣＬ／ＩＰＬ＋ＱＧＣＬ／ＩＰＬ×Ｂ’の合計得点が最高になる経路を探索し、該探索された経路をＧＣＬ／ＩＰＬの境界線Ｌ３として抽出し（ステップＴ１１）、最後に、探索範囲を抽出済みの境界線ＩＰＬ／ＩＮＬ（Ｌ４）〜ＯＰＬ／ＯＮＬ（Ｌ６）の範囲に設定して、パラメータＥＩＮＬ／ＯＰＬ×ＨＩＮＬ／ＯＰＬ＋ＷＩＮＬ／ＯＰＬ×ＧＩＮＬ／ＯＰＬ＋ＱＩＮＬ／ＯＰＬ×Ｂ’の合計得点が最高になる経路を探索し、該探索された経路をＩＮＬ／ＯＰＬの境界線Ｌ５として抽出して（ステップＴ１２）、合計１０本の境界線を抽出する。【００７３】
上述の処理から明らかなように、最初に抽出した二つの境界線内境界膜ＩＬＭ（Ｌ１）と網膜色素上皮ＲＰＥ（Ｌ１０）を除けば、抽出結果を基に探索範囲を限定して他の境界線抽出を行い、またその抽出結果を基に探索範囲を限定して他の境界線抽出を行い、更にその抽出結果を基に探索範囲を限定して他の境界線抽出を行うなど同様な処理を逐次的に繰り返して境界線を抽出している。【００７４】
このように逐次的に行われる境界線の抽出は、いずれの抽出も探索範囲が限定されることから抽出処理を高速化できる他に、範囲変更の都度パラメータ（例えば存在確率、重み付け係数）を適切に設定しなおせるので抽出が容易になるという利点がある。また、後述するように、抽出済みの境界線を利用して探索範囲を設定できるので、抽出済みの境界線と交差することを回避したり、あるいは不明瞭ないし消失した境界線を抽出することができる。【００７５】
また、複数の境界線を抽出する場合、先に抽出した境界線を用いて湾曲補正することができる。例えば、図９のステップＴ５においてＩＳ／ＯＳ（Ｌ８）の境界線を抽出した場合、先に抽出した境界線ＩＬＭ（Ｌ１）あるいはＲＰＥ（Ｌ１０）の湾曲に合わせて入力画像を補正（特に、傾斜補正）したのちに別の境界線の抽出を行うことができる。このような湾曲補正を行うと、エッジおよび輝度勾配の方向が揃うので境界線の抽出精度を向上させることができる。【００７６】
なお、上述した処理において、境界線をユーザーが修正したとき、当該修正した境界線以降に抽出した境界線に対しては、再度抽出処理を行う。例えば、図９のステップＴ６においてＯＰＬ／ＯＮＬ（Ｌ６）を抽出し、この抽出した境界線をユーザーが修正したとき、この修正した境界線を用いて探索範囲を設定して境界線を抽出するステップＴ８、Ｔ９、Ｔ１０、Ｔ１２の処理は、再度抽出処理を行うようにする。【００７７】
＜探索範囲の設定＞
画像処理部３０には、探索範囲設定手段３５が設けられており、この探索範囲設定手段３５により、抽出済みの境界線を利用して境界線の探索範囲を動的に設定することができる。【００７８】
その例が図１０、図１１に示されている。抽出済みの境界線を利用しない場合には、例えば、注目画素Ｉ（ｉ、ｊ）から左側に探索を開始する場合、図１０（ａ）に示したように、注目画素Ｉ（ｉ、ｊ）の左隣りの画素（ｉ−１、ｊ）を中心にして±ｓの（２ｓ＋１）画素幅の画素列が探索範囲として設定される。【００７９】
これに対して、抽出済みの境界線を利用して境界線の探索範囲を設定する場合には、抽出済みの境界線Ｌの傾きに応じて探索範囲を動的に設定する。抽出済みの境界線Ｌの傾きとは、図１０（ｂ）に示したように、画素列ｉにおいて抽出済みの境界線Ｌが抽出された画素と、その画素列ｉの画素からｘ方向に１画素幅ずれた画素列ｉ−１において抽出済みの境界線Ｌが抽出された画素とが、ｚ方向にどれくらいずれているかを示すものである。この抽出済みの境界線Ｌの傾きをｄとして、注目画素Ｉ（ｉ、ｊ）の左隣りの画素（ｉ−１、ｊ）から抽出済みの境界線方向へｄだけシフトした、すなわち境界線方向にｓ＋ｄ、逆側にｓ−ｄの（２ｓ＋１）画素幅の画素列が探索範囲に設定される。同様に、他の探索範囲も抽出済みの境界線に合わせて設定する。【００８０】
このように、探索範囲を既に抽出済みの境界線の傾きに合わせて設定することにより、抽出済みの境界線と同様な曲線の境界線を高精度で抽出することが可能となる。【００８１】
また、図１０（ｂ）に図示したように、探索範囲を、探索範囲の上端画素がその画素列にある抽出済みの画素から所定画素ｄ´隔たるように設定することにより、抽出済みの境界線と交わらないような範囲で探索を行うことができる。従って、抽出した境界線同士が交差することを回避することができる。【００８２】
また、図１１の上部に図示したように、抽出済みの境界線をＬ、Ｌ' 、その間に存在する抽出しようとする境界線をＬ"とした場合、境界線Ｌ"が途中で不明瞭、あるいは消失して不連続になる場合がある。【００８３】
探索範囲を抽出済みの境界線に合わせて設定しない場合には、同図の左側に図示したように、探索範囲として、注目画素Ｉ（ｉ、ｊ）の左隣りの画素（ｉ−１、ｊ）を中心にｚ方向に３画素長さの画素列が設定され、以降３画素長さの画素列が左側に順次探索範囲として設定される（この場合はｓ＝１である）。設定された最終的な探索範囲は、左側の中段に図示したような範囲となる。すると、評価得点が最高になるような経路探索を行った場合、抽出した境界線は、下段に示したように、点線で示したような曲線となり、下方の抽出済みの境界線と交わってしまい、連続した境界線を抽出することができなくなってしまう状況が生じ得る。【００８４】
これに対して、探索範囲を抽出済みの境界線に合わせて設定すると、不明瞭あるいは不連続になった境界線を抽出することが可能になる。これが、図１１の右側に図示されている。ここでは、抽出済みの境界線Ｌに合わせて探索範囲を設定している。つまり、右側に示したように、順次設定される探索範囲は、抽出済みの境界線Ｌに合わせて、抽出済みの境界線Ｌ、Ｌ'の画素と交差しないように、また不明瞭ないし消失した部分につながる画素Ｉ'を探索範囲に含むように、設定される。従って、設定された最終的な探索範囲は、右側の中段に図示したような範囲となる。経路探索を行った場合、抽出した境界線は、下段に点線で示したような曲線となり、消失した部分あるいは不明瞭な部分が連結された境界線を抽出することが可能になる。この場合、探索範囲は抽出済みの境界線Ｌ、Ｌ'の画素と交わらないように設定されているので、抽出した境界線が抽出済みの境界線と交わることを回避することができる。【００８５】
なお、図１０、図１１に示した例では、探索は右から左に行うようにしているが、逆方向に行うようにしてもよい。【００８６】
このように、抽出済みの境界線を利用して探索範囲を適切に設定することにより、不明瞭な境界線、あるいは一部消失した境界線を、抽出済みの境界線と交わることなく抽出することが可能になる。【００８７】
＜制御点の設定＞
上述したような処理を介して抽出された各境界線には、図１２、図１３に図示したように、その境界線上に制御点を設定することができる。【００８８】
例えば、ユーザーが操作部２５のマウス、あるいは操作ペンで境界線の一点を指定し、画素間隔Ｄを指定する。制御部２１は、指定された境界線上の画素を特定し、制御点設定手段３６により当該特定された画素を中心にＸ方向左右にそれぞれＤ×ｎ（ｎ＝１、２、．．．．．）番目の画素列と境界線がそれぞれ交わる画素位置に制御点を設定する。または、制御部２１は、指定された境界線上に既定のＸ方向位置で制御点を設定するようにしても良い。【００８９】
このように設定された制御点は表示部２４に表示される。例えば、図１２の下方左図に図示したように、最下端の網膜色素上皮ＲＰＥ（Ｌ１０）上に５画素間隔で黒丸で示した制御点が設定され、また同中央図に図示したように、１０画素間隔で、また右図に図示したように、５０画素間隔でそれぞれ制御点が設定される。【００９０】
図６に図示した方法で抽出した境界線は画素間隔が狭い（１画素間隔）ので、本来は滑らかな境界線でも、細かい凹凸が発生し滑らかな曲線にはならない。しかし、上述したように、境界線上に５画素間隔、１０画素間隔など広い画素間隔で制御点を設定し、各設定された制御点を、例えばスプライン曲線によって結んで、それを境界線とすることにより、境界線を主観的に把握しやすい滑らかな曲線に近づけることができる。【００９１】
制御点間隔を狭くすると、抽出結果を忠実に表現することが可能になるが、制御点の数が多くなるため、境界線の修正に時間を要する。また、制御点間隔を狭くすると、滑らかな境界線でも、細かい凹凸が発生し滑らかな曲線とはならない。従って、制御点をどの程度の画素間隔で設定するかは、抽出する境界線の特徴に応じて、あるいは抽出された境界線の必要とされる修正程度に応じて、ユーザーが設定するようにする。【００９２】
例えば、図１３の上図は、抽出された境界線ＮＦＬ／ＧＣＬ（Ｌ２）上に制御点を狭い画素間隔で設定し、各制御点をスプライン曲線で結んで境界線とした例を示しており、中段は制御点を広い画素間隔で設定して境界線とした例を示している。【００９３】
制御点間隔を広くする場合には、修正する制御点の数が少なくなり、修正時間を短縮することができる。例えば、図１３下段に示したように、制御点ＱをＱ’に移動させて抽出した境界線を修正することができる。しかし、同様な修正を、図１３上図に示したように制御点間隔を狭くして行う場合には、複数の制御点を移動させなければならず、修正に時間を要することになる。つまり抽出結果に対しての正確性、修正に係る手間、および見た目の滑らかさのトレードオフであり、制御点間隔をユーザーが指定することでこのバランスもユーザーの所望に近づけることができる。【００９４】
なお、上述したように、境界線上に所定の画素間隔で制御点を設定した場合、設定された制御点を一つあるいは複数削除して、残る制御点をスプライン曲線で結ぶようにすることができ、逆に設定された制御点の間に更に一つあるいは複数の制御点を追加して、各制御点をスプライン曲線で結んで境界線とすることもできる。このように、制御点を削除、追加、あるいは移動させることにより、より滑らかな、あるいは忠実な境界線を抽出することが可能になる。【００９５】
また、上述したように、図１３下段には、制御点ＱをＱ’に移動して抽出した境界線を修正する状態が図示されている。このような移動は、一つの制御点だけでなく、複数の制御点を移動させることもできる。また、上述したように、制御点は追加ないし削除することができる。このように、制御点を移動、追加、削除した場合、境界線抽出手段３４によって、一つないし全ての境界線に対して再度経路探索を行い、境界線を抽出しなおすようにしてもよい。このとき、経路探索が移動した制御点、あるいは追加した制御点を通って行われるように、移動ないし追加した制御点が位置する画素には、高い評価得点を与えるようにするとともに、削除した制御点を通過しないようにするために、削除した制御点が位置する画素には、低い評価得点を与えるか、もしくは評価得点を与えないようにしておく。【００９６】
また複数の境界線に対して再度経路探索を行う場合、境界線が交差してしまう現象を防ぐため、境界線Ａ上の制御点Ｑ（ｘ、ｚ）がＱ’（ｘ、ｚ’）へ移動した場合は、境界線Ｂの再検出時において次のような措置をとることができる。例えば図１４に示すように境界線Ａ上の制御点Ｑ（ｘ、ｚ）が境界線Ｂを跨いで−ｚ方向へ移動しＱ’（ｘ、ｚ’）となった場合、境界線Ｂの検出時には、Ｑ’（ｘ、ｚ’）を端点として＋ｚ方向に伸びる半直線領域を探索範囲から除外するようにしてもよい。同様に、境界線Ａ上の制御点Ｑ（ｘ、ｚ）が境界線Ｂを跨いで＋ｚ方向へ移動しＱ’（ｘ、ｚ’）となった場合、境界線Ｂの検出時には、Ｑ’（ｘ、ｚ’）を端点として−ｚ方向に伸びる半直線領域を探索範囲から除外するようにしてもよい。また、図１５に示すように境界線Ａ上の制御点Ｑ（ｘ、ｚ）が境界線Ｂを跨がずに＋ｚ方向へ移動しＱ’（ｘ、ｚ’）となった場合、境界線Ｂの検出時には、Ｑ’（ｘ、ｚ’）を端点としてｚ軸に平行で境界線Ｂと交差しない半直線領域を探索範囲から除外するようにしてもよい。境界線Ａ上の制御点Ｑ（ｘ、ｚ）が境界線Ｂを跨がずに−ｚ方向へ移動した場合も同様である。【００９７】
なお、境界線が修正された場合には、上述したように、修正に応じて存在確率画像格納部２６に格納されている当該境界線用の存在確率画像を修正して当該存在確率画像を学習させるようにすることもできる。
【符号の説明】
【００９８】
１０  断層像撮影装置
２０  画像処理装置
２１  制御部
２２  断層画像形成部
２３  記憶部
２４  表示部
２５  操作部
２６  存在確率画像格納部
２７  重み付け係数格納部
３０  画像処理部
３１  境界線候補画像作成手段
３２  輝度値微分画像作成手段
３３  評価得点画像作成手段
３４  境界線抽出手段
３５  探索範囲設定手段
３６  制御点設定手段
【図１】
2017030056.tif 000003
【図２】
2017030056.tif 000004
【図３】
2017030056.tif 000005
【図４】
2017030056.tif 000006
【図５】
2017030056.tif 000007
【図６】
2017030056.tif 000008
【図７】
2017030056.tif 000009
【図８】
2017030056.tif 000010
【図９】
2017030056.tif 000011
【図１０】
2017030056.tif 000012
【図１１】
2017030056.tif 000013
【図１２】
2017030056.tif 000014
【図１３】
2017030056.tif 000015
【図１４】
2017030056.tif 000016
【図１５】
2017030056.tif 000017
【国際調査報告】
2017030056.tif 000018
2017030056.tif 000019
2017030056.tif 000020