(57)【特許請求の範囲】
【請求項１】
網膜の異なる位置で複数の断層画像を取得する画像取得手段と、
前記取得された複数の断層画像各々において前記網膜に含まれる層を検出する網膜層検出手段と、
前記検出された層における複数の位置で前記検出された層に対する垂直な方向を演算する垂直方向演算手段と、
前記複数の位置で前記演算された方向における前記検出された層の厚さを計測する層厚計測手段と、
前記複数の位置の前記計測された厚さに基づいて、層厚マップを作成するマップ作成手段と、
を有することを特徴とする画像処理装置。
【請求項２】
前記マップ作成手段は、前記網膜の水平方向の面の各位置に対応する前記複数の位置の前記計測された厚さに基づいて、前記層厚マップを作成することを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記マップ作成手段は、前記複数の位置の前記厚さに基づいて、前記水平方向の面における前記各位置以外の位置に対応する前記検出された層の厚さを補間することにより、前記層厚マップを作成することを特徴とする請求項２に記載の画像処理装置。
【請求項４】
前記画像取得手段は、前記網膜上の任意の点を中心として角度を変えながら断層画像を撮影するラディアル・スキャンにより前記複数の断層画像を取得し、
前記層厚計測手段は、前記任意の点を中心とする同心円と前記ラディアル・スキャンにおけるスキャン方向との複数の交点での層厚を用いて、前記同心円上の他の位置の層厚を求めることを特徴とする請求項１乃至３のいずれか１項に記載の画像処理装置。
【請求項５】
網膜の断層画像において前記網膜に含まれる層を検出する網膜層検出手段と、
前記検出された層における複数の位置で前記網膜の層の傾きと前記断層画像の縦横の解像度とに基づいて前記検出された層の厚さを計測する方向を演算する垂直方向演算手段と、
前記複数の位置で前記演算された方向における前記検出された層の厚さを計測する層厚計測手段と、
を有することを特徴とする画像処理装置。
【請求項６】
前記網膜層検出手段は、前記断層画像が表示手段に表示される際の縦横の解像度が同じである場合の断層画像において前記層を検出し、
前記垂直方向演算手段は、前記複数の位置で前記検出された層に対する垂直な方向を前記計測する方向として演算することを特徴とする請求項５に記載の画像処理装置。
【請求項７】
網膜の断層画像から前記網膜の網膜層境界を特定する網膜層検出手段と、
前記断層画像における垂直な方向の解像度と水平な方向の解像度とが等しくされた所定断層画像における前記網膜層境界に対して垂直な方向を演算する垂直方向演算手段と、
前記所定断層画像における前記網膜層境界に対して垂直な方向の実所定長さを測定し、表示される前記断層画像における演算された前記網膜層境界に対して垂直な方向の画像上の表示所定長さと、前記実所定長さと、の相関を求める手段と、
前記表示所定長さと前記実所定長さとの相関を示すスケールを前記表示される断層画像に重畳させる処理手段と、を有することを特徴とする画像処理装置。
【請求項８】
網膜の異なる位置で複数の断層画像を取得する画像取得工程と、
前記取得された複数の断層画像各々において前記網膜に含まれる層を検出する網膜層検出工程と、
前記検出された層における複数の位置で前記検出された層に対する垂直な方向を演算する垂直方向演算工程と、
前記複数の位置で前記演算された方向における前記検出された層の厚さを計測する層厚計測工程と、
前記複数の位置の前記計測された厚さに基づいて、層厚マップを作成するマップ作成工程と、
を有することを特徴とする画像処理方法。
【請求項９】
前記マップ作成工程では、前記網膜の水平方向の面の各位置に対応する前記複数の位置の前記計測された厚さに基づいて、前記層厚マップを作成することを特徴とする請求項８に記載の画像処理方法。
【請求項１０】
前記マップ作成工程では、前記複数の位置の前記厚さに基づいて、前記水平方向の面における前記各位置以外の位置に対応する前記検出された層の厚さを補間することにより、前記層厚マップを作成することを特徴とする請求項９に記載の画像処理方法。
【請求項１１】
前記画像取得工程では、前記網膜上の任意の点を中心として角度を変えながら断層画像を撮影するラディアル・スキャンにより前記複数の断層画像を取得し、
前記層厚計測工程では、前記任意の点を中心とする同心円と前記ラディアル・スキャンにおけるスキャン方向との複数の交点での層厚を用いて、前記同心円上の他の位置の層厚を求めることを特徴とする請求項８乃至１０のいずれか１項に記載の画像処理方法。
【請求項１２】
網膜の断層画像において前記網膜に含まれる層を検出する網膜層検出工程と、
前記検出された層における複数の位置で前記網膜の層の傾きと前記断層画像の縦横の解像度とに基づいて前記検出された層の厚さを計測する方向を演算する垂直方向演算工程と、
前記複数の位置で前記演算された方向における前記検出された層の厚さを計測する層厚計測工程と、
を有することを特徴とする画像処理方法。
【請求項１３】
前記網膜層検出工程では、前記断層画像が表示手段に表示される際の縦横の解像度が同じである場合の断層画像において前記層を検出し、
前記垂直方向演算工程では、前記複数の位置で前記検出された層に対する垂直な方向を前記計測する方向として演算することを特徴とする請求項１２に記載の画像処理方法。
【請求項１４】
網膜の断層画像から前記網膜の網膜層境界を特定する網膜層検出工程と、
前記断層画像における垂直な方向の解像度と水平な方向の解像度とが等しくされた所定断層画像における前記網膜層境界に対して垂直な方向を演算する垂直方向演算工程と、
前記所定断層画像における前記網膜層境界に対して垂直な方向の実所定長さを測定し、表示される前記断層画像における演算された前記網膜層境界に対して垂直な方向の画像上の表示所定長さと、前記実所定長さと、の相関を求める工程と、
前記表示所定長さと前記実所定長さとの相関を示すスケールを前記表示された断層画像に重畳させる処理工程と、
を有することを特徴とする画像処理方法。
【請求項１５】
請求項８乃至１４に記載の画像処理方法の各工程をコンピュータに実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、眼部の画像診断を支援する画像処理装置に関する。
【背景技術】
【０００２】
生活習慣病や失明原因の上位を占める疾病の早期診断を目的として、眼部の検査が広く行われている。光干渉断層計（ＯＣＴ；Ｏｐｔｉｃａｌ  Ｃｏｈｅｒｅｎｃｅ  Ｔｏｍｏｇｒａｐｈｙ）を用いた撮像装置（以下、ＯＣＴ装置とも呼ぶ。）は、眼底の網膜の内部構造の状態を３次元的に観察することが可能であるため、疾病の診断に有用である。眼底の網膜は複数の層構造を示しており、各層の厚みが疾病の進行を示す指標となることが知られている。ＯＣＴ装置により取得した網膜の断層画像を用いることにより、網膜の層構造を定量的に計測することができる。このため、疾病の進行を正確に把握することができるため、疾病の早期発見への期待が大きくなっている。ここで、網膜の断層画像を用いて、各層同士の深度方向（断層画像を表す図３の縦方向）の距離を網膜の層厚として計測する画像処理法が、特許文献１に開示されている。
【０００３】
ところで、近年、アジア地域で多くみられる、近視眼のＯＣＴによる観察に注目が集まっている。近視眼では通常の眼底よりも、網膜の湾曲が大きくなる場合があり、疾病との相関に関心が集まっている。
【０００４】
ここで、図３は、網膜の湾曲が大きい場合について、近視眼眼底の黄斑付近の断層画像の模式図である。Ｌ１〜Ｌ４は、網膜の層構造の境界を示している。Ｌ１は内境界膜とその上部の組織との境界（以下ＩＬＭとする）、Ｌ２は神経線維層とその下部の層との境界（以下ＮＦＬとする）、Ｌ３は視細胞内節概説接合部との上部の層との境界（以下ＩＳ／ＯＳとする）、Ｌ４は網膜色素上皮とその下部の組織との境界（以下ＲＰＥとする）を示している。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００９−６６０１５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
一般的に、網膜の湾曲が大きい場合、断層画像の上部における各層同士の垂線方向（所定の層に対する垂直方向）の距離が短くなる（図３参照）。しかしながら、特許文献１のように、各層同士の深度方向の距離を網膜の層厚として計測すると、湾曲の大きな網膜の状態を正確に把握することは困難である。これでは、疾病の進行を正確に把握することは難しい。
【課題を解決するための手段】
【０００７】
本発明は上記課題を解決するためになされたものであり、
網膜の異なる位置で複数の断層画像を取得する画像取得手段と、
前記取得された複数の断層画像各々において前記網膜に含まれる層を検出する網膜層検出手段と、
前記検出された層における複数の位置で前記検出された層に対する垂直な方向を演算する垂直方向演算手段と、
前記複数の位置で前記演算された方向における前記検出された層の厚さを計測する層厚計測手段と、
前記複数の位置の前記計測された厚さに基づいて、層厚マップを作成するマップ作成手段と、
を有することを特徴とする。
【発明の効果】
【０００８】
本発明によれば、近視等の理由により網膜が大きく湾曲しているような場合にも、網膜の状態をより正確に計測することができる。
【図面の簡単な説明】
【０００９】
【図１】実施形態１に係る画像処理装置１０の機能構成を示す図。
【図２】実施形態１に係る画像処理装置１０の処理手順を示すフローチャート。
【図３】湾曲の大きな網膜断層画像の模式図。
【図４】断層画像上の角度や層厚と、縦横の解像度を補正した画像上の角度や層厚の関係を示す模式図。
【図５】黄斑を中心としたラディアル・スキャンの模式図。
【図６】ラディアル・スキャンから作成した層厚マップのイメージ図。
【図７】実施形態２に係る画像処理装置１０の機能構成を示す図。
【図８】実施形態２に係る画像処理装置１０の処理手順を示すフローチャート。
【図９】湾曲の大きな網膜断層画像上でカーソルの位置にスケールを提示した場合の例。
【発明を実施するための形態】
【００１０】
［実施例１］
本発明は、近視眼など湾曲が大きくなることが知られている疾病眼の網膜の層厚を定量的に計測する場合に、縦横の解像度が異なるＯＣＴ断層画像から正確な層厚を計測することを目的としたものである。より具体的には、断層画像から網膜層を検出してその傾きを求め、該傾きと、ＯＣＴ断層画像の縦横の解像度の値から、層厚を計測する方向を算出して、その方向に計測を行う。このような計測法を用いることにより、湾曲が大きな網膜の場合でも正確な層厚を計測することが可能となり、該層厚により作成される層厚マップがより診断に有効な情報を提供する。
【００１１】
図１は、本実施形態に係る画像処理装置１０の機能構成を示したものである。図中の１００は画像取得部であり、不図示の断層画像取得装置により撮影された画像、もしくはデータベースに保存された画像を取得する。また取得した画像から画像の縦及び横方向の解像度を取得し、制御部１１０を通じて記憶部１６０に保存する。１１０は制御部、１２０は網膜層検出部、１３０は垂直方向演算部、１４０は層厚計測部、１５０はマップ作成部、１６０は記憶部であり、取得された画像に撮影されている網膜の層厚を計測し、結果を保存する。具体的には、網膜層検出部１２０により網膜層を検出し、記憶部１６０に保存された縦横の解像度比を考慮して層厚計測角度を垂直方向演算部１３０で算出する。算出された角度に基づき層厚計測部１４０で正確な層厚を計測し、マップ作成部１４０で層厚マップを作成する。１７０は出力部であり、作成された層厚マップを不図示のモニタ等に出力する他、記憶部１６０に保存されている処理結果を不図示のデータベースに保存する。
【００１２】
次に、図２のフローチャートを参照して、本実施形態の画像処理装置１０の処理手順を説明する。
【００１３】
＜ステップＳ２１０＞
ステップＳ２１０において、画像取得部１１０は、画像処理装置１０に接続される不図示のＯＣＴ断層画像取得装置、もしくは該ＯＣＴ断層画像取得装置により撮像された断層画像を保存するデータベースから、解析対象となる画像を取得する。そして、取得した断層画像を制御部１１０へと送信するとともに、画像から縦横の解像度を取得して記憶部１３０へと送信する。ここで解像度の情報は、Ｔｉｆｆ画像のように画像内のタグに埋め込まれている情報を取得する場合もあれば、ｍｈｄファイルのような形で別ファイルに保存されている情報を同時に取得する場合や、不図示の操作者による入力値を取得する場合もある。
【００１４】
＜ステップＳ２２０＞
ステップＳ２２０において、網膜層検出部１２０は、記憶部１６０に保存された断層画像から網膜層境界を検出する。網膜層検出部１２０は、本発明において網膜の断層画像から網膜層境界を検出し、特定する網膜層検出手段とし作用し、網膜層検出工程を実行する。層のセグメンテーション手法に関しては様々な手法が知られているが、本実施形態ではエッジ強調フィルタを用いて層境界となるエッジを抽出し、その後網膜層に関する医学知識を用いて検出されたエッジと層境界を対応づける手法を用いた場合について説明する。またここではＩＬＭとＲＰＥの検出について述べるが、他の層境界についても、同様の手法により検出することが可能である。
【００１５】
まず網膜層検出部１２０は、断層画像に対して平滑化フィルタ処理を行い、ノイズ成分を除去する。そして、エッジ検出フィルタ処理を行い断層画像からエッジ成分を検出し、層の境界に相当するエッジを抽出する。さらに、エッジ検出をした二次元断層画像から背景領域を特定し、背景領域の輝度値特徴を二次元断層画像から抽出する。そして、次に、エッジ成分のピーク値と、ピーク間の輝度値特徴を用いることで、各層の境界を判定する。
【００１６】
例えば網膜層検出部１２０は、硝子体側から眼底の深度方向にエッジを探索し、エッジ成分のピークと、その上下の輝度特徴と、背景の輝度特徴から、硝子体と網膜層との境界（ＩＬＭ）を判定する。さらに、眼底の深度方向にエッジを探索し、エッジ成分のピークと、ピーク間の輝度特徴と、背景の輝度特徴を参照して、網膜色素上皮層境界（ＲＰＥ）を判定する。以上の処理によって、層の境界を検出することができる。
【００１７】
このようにして検出した各層の境界（コントロールポイント）を制御部１１０へ送信するとともに、記憶部１３０に保存する。
【００１８】
＜ステップＳ２３０＞
ステップＳ２３０において、垂直方向演算部１３０は、ステップＳ２２０において検出されたＲＰＥとＩＬＭを用いて、網膜層に対して垂直となる方向を算出する。ここで網膜層に対して垂直となる方向とは、縦横の解像度が等しい画像において、網膜に水平な方向と直行する方向とする。ステップＳ２２０において検出されたＲＰＥの各点に対して、ＲＰＥに接する方向と水平方向のなす角θを計算する。さらに、このθから以下の式を用いて網膜層に対して垂直となる方向が、画像の垂直方向となす角度θ”を算出する。
【００１９】
【数１】
0005952564.tif 000002
【００２０】
このように各コントロールポイントに対して演算された角度θ及びθ”を、制御部１１０へ送信するとともに、記憶部１３０に保存する。垂直方向演算部１３０は、本発明において特定された前記網膜層境界に対して垂直な方向の解像度と前記網膜層境界に対して平行の方向の解像度とが等しくされた所定断層画像、即ち縦横比を反映した断層画像における網膜層の境界に対して垂直な方向を演算する垂直方向演算手段として作用し、垂直方向演算工程を実行する。
【００２１】
＜ステップＳ２４０＞
ステップＳ２４０において、層厚計測部１４０は、ステップS２３０において算出された網膜層に対して垂直となる方向の角度θ”に基づき、層厚を計測する。具体的には、ステップＳ２２０で求められたＲＰＥの各点から、ステップＳ２３０で算出された角度の方向に、ＩＬＭまでの距離Ｌを求める。そして求められた距離Ｌから、以下の式を用いて層厚Ｌ’を算出する。
【００２２】
【数２】
0005952564.tif 000003
このように各コントロールポイントに対して計測された層厚Ｌ’を、制御部１１０へ送信するとともに、記憶部１３０に保存する。
即ち、層厚計測部１４０は、本発明における、所定の長さである網膜層の厚さについて、所定断層画像における前記垂直な方向の前記所定長さである実所定長さを測定する、後述する相関を求める手段の一部を構成する。ここで述べる実所定長さとは、網膜層境界に対して垂直な方向（縦）の解像度と、水平の方向（横）の解像度とが等しいところの本発明において定義する所定断層画像における網膜層境界に対して設定された該網膜層境界に垂直な方向での網膜層の厚さに対応する。
【００２３】
ここで、上記ステップＳ２３０、Ｓ２４０において用いた（式１）、(式２)について、このようにして角度や層厚を計算する理由を説明する。
【００２４】
図４においてＬ１、Ｌ４は、補正前の、縦横で解像度が異なる場合のＩＬＭ、ＲＰＥをそれぞれ示していて、Ｌ１’、Ｌ４’は、縦横の解像度が等しくなるように補正した後の補正ＩＬＭ、補正ＲＰＥを示している。点Ａは、ＲＰＥ上のある点を表し、点ＡにおけるＲＰＥの水平方向とのなす角をθとする。点Ａに対応する補正ＲＰＥ上の点を点Ａ’とする。この点Ａ’における補正ＲＰＥの水平方向とのなす角をθ’とする。ここで層厚の計測は、縦横の解像度が等しくなるように補正した画像において、網膜(ＲＰＥ)に対して直角方向に計測した場合のＲＰＥからＩＬＭまでの距離とする。この場合、点Ａ’における補正ＲＰＥに対する直角方向は、点Ａ’から垂直方向に対してθ’の方向となる。この点Ａ’から補正ＲＰＥに対して直角方向にひかれた直線が、補正ＩＬＭと交わる点を点Ｂ’とする。そしてこの点Ｂ’に対応する、ＩＬＭ上の点を点Ｂとする。これより、補正前の画像において点Ａにおける層厚を計測したい場合は、点Ａから点Ｂ、つまり点Ａに対して垂直方向から角度θ”の方向に計測すればよいことがわかる。
【００２５】
ここで、θとθ”の間には以下のような関係がある（式１）。
【数３】
0005952564.tif 000004
ここでｒは、横方向に対する縦方向の解像度の大きさであり，補正前の画像は、縦方向にｒ倍に引き伸ばされて表示されていることを示す。
この関係を用いれば、ステップＳ２２０において検出されたＲＰＥに対して、網膜層に対して垂直となる方向の角度θ”を計測することができる。
【００２６】
さらに、層厚に関しても、点Ａと点Ｂの距離をＬ、点Ａ’と点Ｂ’の距離をＬ’としたとき、ＬとＬ’の間には以下のような関係がある（式２）。
【数４】
0005952564.tif 000005
これにより、補正前の画像を用いて計測した距離Ｌから、補正後の距離Ａ’Ｂ’を求めることができる。
【００２７】
以上の工程は、断層画像における演算された垂直な方向の画像上の所定長さである表示上の所定長さである表示所定長さ、即ち網膜層の表示厚さと、具体的な測定により得られる実際の厚さである実所定長さ（厚さ）と、の相関を求める、上述した相関を求める手段により為される。
【００２８】
＜ステップＳ２５０＞
ステップＳ２５０において、マップ作成部１５０は、ステップＳ２４０において算出された網膜の層厚（ここではＲＰＥからＩＬＭまでの距離）を記憶部１３０より取得し、層厚マップを作成する。ここでは、ある点を中心として角度を変えながら断層画像を撮影するラディアル・スキャンによって撮影を行った場合を例として、層厚マップの作成法を示す。図５は、黄斑を中心として、１５度ごとに１２枚のラディアル・スキャンを行った場合の、撮像位置を示している。ここでは水平方向のスキャンをスキャン１、そこから１５度回転させたスキャンをスキャン２としている。
【００２９】
ステップＳ２２０〜Ｓ２４０によって、１２枚のラディアル・スキャン上で、層厚の算出値が取得されている。ここでは、１０ピクセルおきにコントロールポイントが設定されていて、各コントロールポイント上に層厚の値が得られているとする。図６（ａ）において、Ａ１、Ｂ１はスキャン１上のコントロールポイント、Ａ２，Ｂ２はスキャン２上のコントロールポイントとする。Ａ１とＡ２はラディアル・スキャンの中心を中心とした同心円上に配置されるから、Ａ１とＡ２における層厚の値をもとにして、Ａ１とＡ２を結ぶ同心円上の辺Ａ１Ａ２上の層厚値を補間する。ここで補間の方法としては様々な方法が考えられる。Ａ１とＡ２からの距離に応じた線形補間や、Ａ１とＡ２・・・Ａ１２を考慮したスプライン補間などが考えられる。このようにして、同心円上に値を補間することにより、２次元の層厚マップを作成する。図６（ｂ）に作成した層厚マップの模式図を示す。
【００３０】
以上の工程を実行するマップ作成部１５０は、本発明における表示上の網膜層厚さと実際に測定、計測された網膜層厚さを用いてこれらの相関を反映させた断層画像である層厚マップを作成する処理手段として作用する。
【００３１】
＜ステップＳ２６０＞
ステップＳ２６０において、出力部１７０は、ステップＳ２５０で作成した層厚マップを不図示のモニタに表示するとともに、記憶部１６０に保存された各ラディアル・スキャンの各コントロールポイントにおける、層厚を計測する方向を示す角度及び層厚の値を不図示のデータサーバに保存する。
【００３２】
以上の構成により、湾曲の大きな網膜の断層画像からでも、正確な層厚を反映した層厚マップを提示することができるという効果がある。
【００３３】
ここで、網膜の所定の層に対する垂直方向における各層同士の距離を網膜の層厚として計測することが、本出願人により提案されている（特願2009-199018参照）。このとき、断層画像において網膜が傾いて表示されている場合でも、正確に網膜の状態を把握することができる。
【００３４】
ところで、一般的に、断層画像における縦方向（所定の層に対する垂直方向）の解像度は、横方向（所定の層に平行な方向）よりも高い。このとき、縦方向と横方向との解像度を等しくするため、表示される断層画像は、実際の網膜の形状を縦方向に拡大した状態となる。この場合、全体的に層が潰れた印象を与える画像になる。また、上述した本出願人による提案（特願2009-199018）の手法で層厚を計測しても、網膜の状態を正しく把握しにくい場合がある。これに対し、本実施例では、縦方向の解像度を横方向に合わせない（解像度比の補正を行わない）断層画像の層厚マップを提示することとしている。これにより、医師は、見易い画像を用いて網膜内部の詳細の観察を行うことが可能となり、より正確に診断することができる。
【００３５】
［実施例２］
実施例１では、縦横の解像度が異なるＯＣＴ断層画像の湾曲した網膜において、正確な層厚を計測し、層厚マップを作成して提示する例を示した。しかし湾曲した網膜において、断層画像をより詳細に観察したい場合もある。本実施例では、断層画像上に網膜に対して垂直となる方向や、網膜と水平方向と垂直方向での長さの比を示す例について説明する。
【００３６】
図７は、本実施形態に係る画像処理装置１０の機能構成を示したものである。但し図７の、画像取得部１００、制御部１１０、網膜層検出部１２０、記憶部１６０、出力部１７０は図１と同じであるため説明は省略する。入力部７０５は不図示の操作者により不図示のモニタ上で断層画像上に提示された位置、例えばカーソルの位置を取得する。領域分割部７３０は実施例で算出した層厚を計測する角度に基づき、網膜領域を分割する。長さ比演算部７４０は、層検出部１２０により検出された網膜層上において、網膜層と水平方向及び垂直方向の長さの変化量を算出する。
【００３７】
次に、図８のフローチャートを参照して、本実施形態の画像処理装置１０の処理手順を説明する。但し図８のステップＳ８１０〜Ｓ８２０は、図２のステップＳ２１０〜Ｓ２２０と同じであるため説明は省略する。
【００３８】
なお、網膜層検出部１２０は、本発明において網膜の断層画像から網膜層境界を検出し、特定する網膜層検出手段とし作用し、網膜層検出工程を実行する。
【００３９】
＜ステップＳ８３０＞
ステップＳ８３０において、領域分割部７３０は、実施例１のステップＳ２３０と同じ方法で網膜層に対して垂直となる方向の角度θ”を算出し、該角度に基づき、網膜領域を分割する。具体的には、ステップＳ８２０で検出し保存されているＲＰＥの各コントロールポイントに対して算出した角度θ”の方向にひいた直線がＩＬＭと交差する点を求める(図４のＢに対応)。ＲＰＥ上のコントロールポイントをＡｎとし、対応するＩＬＭ上の交点をＢｎとすれば、網膜内領域（ＲＰＥとＩＬＭに囲まれた領域）を領域Ｒｎ（Ａｎ，Ｂｎ，Ａｎ＋１，Ｂｎ＋１で囲まれた領域）に分割することができる。
【００４０】
なお、本実施例では、実施例１における垂直方向演算部１３０に換えて領域分割部７３０が、本発明において特定された前記網膜層境界に対して垂直な方向の解像度と前記網膜層境界に対して平行の方向の解像度とが等しくされた所定断層画像、即ち縦横比を反映した断層画像における網膜層の境界に対して垂直な方向を演算する垂直方向演算手段として作用し、垂直方向演算工程を実行する。また、該領域分割部７３０は、演算された垂直な方向に基づいて網膜層を複数の領域に分割する工程を実施する領域分割手段として作用し、前述した相関を求める手段に包含される。
【００４１】
＜ステップＳ８４０＞
ステップＳ８４０において、入力部７０５は、不図示の操作者により不図示のモニタ上に指定された位置を取得し、制御部１１０を通じて記憶部１６０に保存する。具体的には、不図示のモニタ上にＯＣＴ断層画像が提示されている場合に、不図示の操作者が断層画像上にカーソルを移動した場合のカーソル位置Ｃを取得することによって得られる。
【００４２】
＜ステップＳ８５０＞
ステップＳ８５０において、制御部１１０は、ステップＳ８４０で取得されたカーソル位置Ｃが、網膜内領域に含まれるか否かを判定する。判定の方法としては、カーソル位置Ｃとして取得された座標が、ステップＳ８２０で検出されたＲＰＥとＩＬＭに挟まれた領域に含まれるかで判定するなどの方法がある。そして網膜内領域に含まれる場合には、ステップＳ８３０で分割された領域のうち、カーソル位置Ｃを含む領域（Ｒｍとする）を求め、ステップＳ８６０に進む。含まれない場合には、ステップＳ８８０に進む。
【００４３】
＜ステップＳ８６０＞
ステップＳ８６０において、長さ比演算部７４０は、ステップＳ８４０で取得された網膜層内のカーソル位置Ｃに対して、網膜層に水平な方向とのなす角θと、網膜層に垂直となる方向の角度θ”それぞれの方向について、単位長さｄの線分の画像上での長さを算出する。この比例演算工程は、前述した相関を求める手段に含まれた、分割された領域内の点に対して単位長さの縦横の画像上の比を演算する比例演算手段の一例である長さ比演算部７４０により実施される。
【００４４】
ここで、もしもカーソル位置ＣがＲＰＥ上のコントロールポイントと一致している場合には、該コントロールポイントに対してすでに取得されているθ、θ”を用いて、θ方向の長さｄ’は、
【数５】
0005952564.tif 000006
【００４５】
θ”方向の長さｄ”は、
【数６】
0005952564.tif 000007
と表わされる（ステップＳ２４０及び図４の説明を参照）。（式３）、（式４）によれば、網膜が水平方向（θ＝０）であるとき、ｄ’＝ｄ、ｄ”＝ｒｄとなり、水平方向は単位長さｄのまま、垂直方向は単位長さがｒ倍に拡大された値となっている。
【００４６】
カーソル位置がコントロールポイントに一致しない場合には、ステップＳ８５０において取得したカーソル位置を含む領域Ｒｍに対して、ＡｍとＡｍ＋１の値より補間した値を算出する。補間の手法は複数の方法が考えられるが、ここでは線形補間の場合を例として説明する。但し補間の方法が線形補間に限定されるものではない。
【００４７】
例えばカーソル位置Ｃから線分ＡｍＡｍ＋１に下ろした垂線が、線分ＡｍＡｍ＋１をｈ：１−ｈに分割する場合に、カーソル位置Ｃにおける角度を以下のように近似する。
【数７】
0005952564.tif 000008
上記近似値を用いて、（式３）、（式４）より画像上での長さｄ’、ｄ”を取得する。
【００４８】
＜ステップＳ８７０＞
ステップＳ８７０において、出力部１７０は、ステップＳ８６０において算出されたθ及びθ”の方向における、単位長さｄの画像上での見た目の長さｄ’、ｄ”に基づき、断層画像上にスケールを提示する。例えば、カーソル位置Ｃを原点として、画像の水平方向から角度θの方向に長さｄ’の線分を引き、垂直方向から角度θ”の方向に長さｄ”の線分を引いて、断層画像上に表示する。図９はスケールが提示された断層画像の一例を示しており、図中矢印で示されるカーソルの先端部を原点として個々の領域でのスケールが表示されている。線分の向きに関しては、カーソル位置ＣがＲＰＥに近い場合にはｄ”の線分の向きを上向きに、ＩＬＭに近い場合には下向きに、画面の右側にある場合には、ｄ’の線分の向きを左向きに、左側にある場合には右向きにするなどが考えられる。
【００４９】
以上の表示する工程は、上述した処理手段に含まれた、演算した長さ比に基づいて作成したスケールを網膜像の上の点に表示する制御手段１１０及び出力手段１７０によって実施される。
【００５０】
＜ステップＳ８８０＞
ステップＳ８８０において、制御部１１０は、断層画像による観察が終了したが否かを判定する。判定の方法は、不図示の操作者により終了を示す入力を受け取ることで行われる。断層画像の観察が終了していない場合には、ステップＳ８４０に戻って再度カーソル位置を取得する。カーソル位置に変更がある場合にはステップＳ８５０〜Ｓ８７０の操作を繰り返して新たなスケールを提示する。観察が終了したという入力があった場合には、処理を終了する。
【００５１】
（その他の実施形態）
本発明の目的は、前述した実施形態の機能を実現するソフトウェアのプログラムによって、或いは該プログラムコードを記憶した記憶媒体を、システムあるいは装置に供給する。そして、そのシステムあるいは装置のコンピュータ（またはＣＰＵやＭＰＵ）が記憶媒体に格納されたプログラムコードを読出し実行することによっても、達成されることは言うまでもない。
【００５２】
（その他の実施例）
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。
【符号の説明】
【００５３】
１０      画像処理装置
１００    画像取得部
１１０    制御部
１２０    網膜層検出部
１３０    垂直方向演算部
１４０    層厚計測部
１５０    マップ作成部
１６０    記憶部
１７０    出力部
７０５    入力部
７３０    領域分割部
７４０    長さ比演算部
【図１】
0005952564.tif 000009
【図２】
0005952564.tif 000010
【図３】
0005952564.tif 000011
【図４】
0005952564.tif 000012
【図５】
0005952564.tif 000013
【図６】
0005952564.tif 000014
【図７】
0005952564.tif 000015
【図８】
0005952564.tif 000016
【図９】
0005952564.tif 000017