(57)【要約】【課題】  計測者の主観や計測箇所によるバラツキの問題を解消しながら、血流幅比率を客観的に計測することが可能な眼底画像処理装置を提供する。【解決手段】  眼底画像と、外形が四角形である計測用図形を表示する表示手段２０、外部指示に基づき、眼底画像に重ねて配置された計測用図形に基づき計測範囲を設定する計測範囲設定手段４０、計測用図形の対向する２辺について、各辺に所定の間隔でサンプルを設定し、各サンプルに対応する座標を算出して１次元の座標配列を作成し、座標配列の各サンプルについて、眼底画像より対応する輝度値を取得して輝度配列を作成し、輝度配列を用いて隣接するサンプルの輝度値との差分値を用いて微分輝度値の配列である微分配列を作成し、微分配列に対して最大値をとる最大サンプルと最小値をとる最小サンプルを探索し、最大サンプルと最小サンプルを基に血流領域端点を特定する血流領域端点特定手段５０、を有する。【選択図】  図２
2018011726.tif 000002
【特許請求の範囲】
【請求項１】
眼底画像に対して処理を行う装置であって、
前記眼底画像と、外形が台形である計測用図形を表示する表示手段と、
外部からの指示に基づき、眼底画像に重ねて配置された前記計測用図形に基づき、当該計測用図形の平行な２辺の少なくとも一方の辺を計測範囲として設定する計測範囲設定手段と、
前記計測範囲として設定された辺について、前記辺をサンプリングし、得られた各サンプルに対応する座標を算出して１次元の座標配列を作成し、
前記座標配列の各座標を用いて前記眼底画像より輝度値を取得し、各サンプルに対応する輝度配列を作成し、
前記輝度配列を用いて隣接するサンプルの輝度値との差分値を用いて微分輝度値の配列である微分配列を作成し、
前記微分配列に対して最大値をとる最大サンプルと最小値をとる最小サンプルを探索し、
前記最大サンプルと最小サンプルを基に血流領域の２端点を特定する血流領域端点特定手段と、
を有することを特徴とする眼底画像処理装置。
【請求項２】
動静脈径比を算出する動静脈径比算出手段を更に有し、
前記計測範囲設定手段は、外部からの指示に基づき、静脈に対応する計測範囲と動脈に対応する計測範囲の各々を設定するものであり、
前記動静脈径比算出手段は、動脈に対応する計測範囲について算出された血流領域の２端点間の距離と、静脈に対応する計測範囲について算出された血流領域の２端点間の距離との比率に基づいて、動静脈径比を算出することを特徴とする請求項１に記載の眼底画像処理装置。
【請求項３】
静脈径比を算出する静脈径比算出手段を更に有し、
前記計測範囲設定手段は、外部からの指示に基づき、静脈に対応して少なくとも２箇所の計測範囲を設定するものであり、
前記計測範囲として設定された計測用図形の平行な２辺により設定された２箇所の計測範囲について算出された２つの血流領域の２端点間の距離どうしの比率に基づいて、静脈径比を算出することを特徴とする請求項１または請求項２に記載の眼底画像処理装置。
【請求項４】
外部からの指定に基づいて、前記眼底画像の指定された領域をトリミングし、当該領域が所定の画素数になるように前記眼底画像の拡大処理を行う画像拡大手段を更に備えることを特徴とする請求項１から請求項３のいずれか一項に記載の眼底画像処理装置。
【請求項５】
前記血流領域端点特定手段は、
輝度配列を作成する際、前記眼底画像がＲＧＢの各成分が所定の階調をもつカラー画像の場合、前記眼底画像に対して、Ｇ成分の階調値とＢ成分の階調値を反転した値との積の平方根によりグレースケール画像に変換し、更にネガポジ反転した階調値を前記輝度値として与えることを特徴とする請求項１から請求項４のいずれか一項に記載の眼底画像処理装置。
【請求項６】
前記計測用図形は、前記台形の平行な２辺の中点を結ぶ基準線を備えることを特徴とする請求項１から請求項５のいずれか一項に記載の眼底画像処理装置。
【請求項７】
前記血流領域端点特定手段は、
輝度配列を作成する際、
前記計測用図形の平行な２辺のいずれかの辺をサンプリングし、前記辺に直交する２方向に各サンプルを所定の距離だけ移動させた点を２箇所特定し、特定された２点における輝度値を用いて前記各サンプルに対応する輝度値を補正することを特徴とする請求項１から請求項６のいずれか一項に記載の眼底画像処理装置。
【請求項８】
前記血流領域端点特定手段は、
前記最大サンプルより微分輝度値の小さい極大点が前記最大サンプルから見て前記最小サンプルの方向に存在する場合、前記極大点を最大サンプルとして置き換え、前記最小サンプルより微分輝度値の大きい極小点が前記最小サンプルから見て前記最大サンプルの方向に存在する場合、前記極小点を最小サンプルとして置き換えることを特徴とする請求項１から請求項７のいずれか一項に記載の眼底画像処理装置。
【請求項９】
前記血流領域端点特定手段は、
前記最大サンプルより微分輝度値が小さく第１所定値より微分輝度値が大きい極大点から見て先頭のサンプル寄りで第１所定値より小さい微分輝度値をとるサンプルを前記血流領域の第１の端点として特定し、前記最小サンプルより微分輝度値が大きく第１所定値より小さい第２所定値より微分輝度値が小さい極小点から見て末尾のサンプル寄りで第２所定値より大きい微分輝度値をとるサンプルを前記血流領域の第２の端点として特定することを特徴とする請求項１から請求項８のいずれか一項に記載の眼底画像処理装置。    
【請求項１０】
前記表示手段は、前記平行な２辺の各々に特定された前記血流領域の４端点を頂点とする外形が台形の図形を、前記眼底画像に重ねて表示することを特徴とする請求項１から請求項９のいずれか一項に記載の眼底画像処理装置。
【請求項１１】
コンピュータを、請求項１から請求項１０のいずれか一項に記載の眼底画像処理装置として機能させるためのプログラム。
【請求項１２】
眼底画像に含まれる血管の血流領域を特定する際に、計測用図形を表示する表示装置であって、
外部からの指示に基づき、外形が台形であり、平行な２辺の中点を結ぶ基準線を備える計測用図形を表示し、
外部からの指示に基づき配置された計測用図形の平行な２辺それぞれを計測範囲として設定し、設定された計測範囲の２辺に基づいて血流領域の２端点を各々特定し、特定された血流領域の４端点を４頂点とする外形が台形形状の図形を前記計測用図形に代えて表示することを特徴とする計測用図形表示装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、眼底画像の処理技術に関し、特に、動脈硬化や高血圧の指標の１つである動静脈径比と動静脈交叉部の静脈径比を、眼底画像を処理することにより解析するための技術に関する。
【背景技術】
【０００２】
生活習慣病として代表的な高血圧・脂質異常症・糖尿病・動脈硬化などいわゆるメタボリックシンドロームにおいて、血圧・脂質・血糖は簡便に計測でき自己検診も可能である。しかし、動脈硬化の計測（通称：血管年齢検査）については、現状ではPWV検査（脈波伝播速度）、ABI検査（足関節上腕血圧比）、頚動脈エコー検査など、検査設備を備えた循環器専門のクリニックでないと計測ができない。これに対して、眼科クリニックが保有する眼底カメラを用いて眼底血管口径を計測することにより、高血圧や動脈硬化を簡便に計測する手法も知られている。眼底は人体の中で唯一血管を直接観察できる箇所であり、眼底動脈の硬化が全身の動脈と同様に進行すると仮定すれば、眼底写真の撮影により動脈硬化の計測が可能になる。眼底カメラは小型化・低価格化が進み、海外では既にスマートフォンでも撮影できる特殊なレンズが販売されており、眼底写真の自己撮影のインフラも整いつつある。
【０００３】
眼底写真を用いた高血圧や動脈硬化の計測手法として、Scheie分類が提案されており、動静脈径比、動静脈交叉部の静脈径比、動脈の血柱反射比の３者が主要な計測項目であり、これらに対応した計測支援手法も提案されている（特許文献１参照）。また、動静脈径比の計測を完全自動で行う手法も提案されている（特許文献２参照）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００７−３１９４０３号公報
【特許文献２】特開平１０−２４３９２４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
特許文献１に記載の技術では、血管幅の計測について、ユーザにより円形で指示された中心から、隣接画素との画素値の差が大きい箇所を境界とみなし、最近傍の境界点までの距離をもとに、血管の半径を算出している。しかし、眼底の病態により血管境界がぼけている場合、境界点を特定できず、血管幅が小さ目に計測される問題がある。特許文献２に記載の技術では、動脈と静脈では輝度に差があり（背景に比べ双方とも輝度が低いが、静脈は動脈より顕著に低い）、同一のしきい値で二値化して血管抽出すると、動脈幅は小さ目に静脈幅は大き目になる。そのため、動静脈径比が正しく計測できない。また、動静脈判別を自動で行う手法を提案し、誤判定箇所を血管の連続性により自動修正する機能も提案されているが、血管の分岐箇所では連続性の判定も誤る可能性がある。
【０００６】
上記のような問題により、現状では眼科医が眼底写真の上にゲージを当てて目視で計測する手法が一般的である。そのため、動静脈径比は１／２、１／３といった大まかな数値になり、進行度を定量的に比較できないという問題がある。一方、読み取り精度を上げても、計測者の主観によるバラツキ、計測箇所によるバラツキが多くなり、多くの計測箇所をサンプリングする必要があるが、眼科医の作業負荷が大きく現実的ではない。また、血流幅比率としては、動脈と静脈の間での動静脈径比だけでなく、動静脈交叉部における静脈の血流狭窄比率を求めたいという要望もある。
【０００７】
そこで、本発明は、計測者の主観や計測箇所によるバラツキの問題を解消しながら、血流幅比率を客観的に計測することが可能な眼底画像処理装置を提供することを課題とする。
【課題を解決するための手段】
【０００８】
上記課題を解決するため、本発明第１の態様では、
眼底画像に対して処理を行う装置であって、
前記眼底画像と、外形が台形である計測用図形を表示する表示手段と、
外部からの指示に基づき、眼底画像に重ねて配置された前記計測用図形に基づき、当該計測用図形の平行な２辺の少なくとも一方の辺を計測範囲として設定する計測範囲設定手段と、
前記計測範囲として設定された辺について、前記辺をサンプリングし、得られた各サンプルに対応する座標を算出して１次元の座標配列を作成し、
前記座標配列の各座標を用いて前記眼底画像より輝度値を取得し、各サンプルに対応する輝度配列を作成し、
前記輝度配列を用いて隣接するサンプルの輝度値との差分値を用いて微分輝度値の配列である微分配列を作成し、
前記微分配列に対して最大値をとる最大サンプルと最小値をとる最小サンプルを探索し、
前記最大サンプルと最小サンプルを基に血流領域の２端点を特定する血流領域端点特定手段と、
を有することを特徴とする眼底画像処理装置を提供する。
ここで、計測用図形の外形は台形であるが、台形とは少なくとも対向する一組の辺が平行であることを意味し、４頂点が直角になる長方形を含む意味である。
【０００９】
本発明第１の態様によれば、外部からの指示に基づき、眼底画像に重ねて配置された計測用図形に基づき、計測用図形の平行な２辺の少なくとも一方の辺を計測範囲として設定し、計測用範囲として設定された辺について、辺をサンプリングし、各サンプルに対応する座標を算出して１次元の座標配列を作成し、座標配列の各座標を用いて眼底画像より輝度値を取得し、各サンプルに対応する輝度配列を作成し、輝度配列を用いて隣接するサンプルの輝度値との差分値を用いて微分輝度値の配列である微分配列を作成し、微分配列に対して最大値をとる最大サンプルと最小値をとる最小サンプルを探索し、最大サンプルと最小サンプルを基に血流領域の２端点を特定するようにしたので、計測者の主観の問題を解消しながら、血流幅を測定するための血流領域を的確に特定することが可能となる。特に、計測用図形の平行な２辺の双方を計測範囲として設定した場合には、計測用図形の平行な２辺の各々に基づいて、２箇所の血流領域の２端点を同時に特定することができるため、２箇所の計測値の平均値を求めることにより、計測箇所によるバラツキを軽減したり、単一血管における２箇所の口径比率を算出できる。また、血流領域を特定するために、２次元方向に処理を行うことなく、１次元の配列を用いて処理を行うため、少ない処理負荷で血流領域の２端点を特定することが可能となる。
【００１０】
また、本発明第２の態様では、動静脈径比を算出する動静脈径比算出手段を更に有し、
前記計測範囲設定手段は、外部からの指示に基づき、静脈に対応する計測範囲と動脈に対応する計測範囲の各々を設定するものであり、
前記動静脈径比算出手段は、動脈に対応する計測範囲について算出された血流領域の２端点間の距離と、静脈に対応する計測範囲について算出された血流領域の２端点間の距離との比率に基づいて、動静脈径比を算出することを特徴とする。
【００１１】
本発明第２の態様によれば、外部からの指示に基づき、静脈に対応する計測範囲と動脈に対応する計測範囲の各々を設定し、動脈に対応する計測範囲について算出された血流領域の２端点間の距離と、静脈に対応する計測範囲について算出された血流領域の２端点間の距離との比率に基づいて、動静脈径比を算出するようにしたので、動静脈径比を的確に算出することが可能となる。
【００１２】
また、本発明第３の態様では、静脈径比を算出する静脈径比算出手段を更に有し、
前記計測範囲設定手段は、外部からの指示に基づき、静脈に対応して少なくとも２箇所の計測範囲を設定するものであり、
前記計測範囲として設定された計測用図形の平行な２辺により設定された２箇所の計測範囲について算出された２つの血流領域の２端点間の距離どうしの比率に基づいて、静脈径比を算出することを特徴とする。
【００１３】
本発明第３の態様によれば、外部からの指示に基づき、静脈に対応して少なくとも２箇所の計測範囲を設定し、前記計測範囲として設定された計測用図形の平行な２辺により設定された２箇所の計測範囲について算出された２つの血流領域の２端点間の距離どうしの比率に基づいて、静脈径比を算出するようにしたので、動静脈交叉部における静脈径比を簡便かつ的確に算出することが可能となる。
【００１４】
また、本発明第４の態様では、外部からの指定に基づいて、前記眼底画像の指定された領域をトリミングし、当該領域が所定の画素数になるように前記眼底画像の拡大処理を行う画像拡大手段を更に備えることを特徴とする。
【００１５】
本発明第４の態様によれば、外部からの指定に基づいて、眼底画像の指定された領域をトリミングし、その領域が所定の画素数になるように拡大処理を行い、拡大された眼底画像に対して計測範囲を設定するようにしたので、画素数の単位で算出される動脈と静脈の２種類の血流の幅が大きな整数値になり、これらの比率である動静脈径比や、動静脈交叉部における静脈径比を算出する際に、整数値どうしの割り算による丸め誤差を小さくすることが可能となる。
【００１６】
また、本発明第５の態様では、前記血流領域端点特定手段は、
輝度配列を作成する際、前記眼底画像がＲＧＢの各成分が所定の階調をもつカラー画像の場合、前記眼底画像に対して、Ｇ成分の階調値とＢ成分の階調値を反転した値との積の平方根によりグレースケール画像に変換し、更にネガポジ反転した階調値を前記輝度値として与えることを特徴とする。
【００１７】
本発明第５の態様によれば、輝度配列を作成する際、眼底画像がＲＧＢの各成分が所定の階調をもつカラー画像の場合、眼底画像に対して、Ｇ成分の階調値とＢ成分の階調値を反転した値との積の平方根によりグレースケール画像に変換し、更にネガポジ反転した階調値を輝度値として与えるようにしたので、血流領域の背景に対するコントラストを向上させ、血流領域の端点を精度よく特定することが可能となる。
【００１８】
また、本発明第６の態様では、前記計測用図形は、前記台形の平行な２辺の中点を結ぶ基準線を備えることを特徴とする。
【００１９】
本発明第６の態様によれば、計測用図形は、台形の平行な２辺の中点を結ぶ基準線を備えるようにしたので、利用者が計測用図形を血管走行に沿って設定する操作を容易にするとともに、設定される計測用図形の平行な２辺が血管走行方向と直交し、２箇所の血流領域の２端点を基に血流幅を高精度に計測することが可能となる。
【００２０】
また、本発明第７の態様では、前記血流領域端点特定手段は、
輝度配列を作成する際、
前記計測用図形の平行な２辺のいずれかの辺をサンプリングし、前記辺に直交する２方向に各サンプルを所定の距離だけ移動させた点を２箇所特定し、特定された２点における輝度値を用いて前記各サンプルに対応する輝度値を補正することを特徴とする。
【００２１】
本発明第７の態様によれば、輝度配列を作成する際、計測用図形の平行な２辺のいずれかの辺をサンプリングし、前記辺に直交する２方向に各サンプルを所定の距離だけ移動させた点を２箇所づつ特定し、特定された２点における輝度値を用いて当該サンプルに対応する輝度値を補正するようにしたので、計測用図形の平行な２辺のいずれかのサンプルにスポット的な輝度ムラがあったとしても、その影響を緩和し、正しい計測を行うことが可能となる。
【００２２】
また、本発明第８の態様では、前記血流領域端点特定手段は、
前記最大サンプルより微分輝度値の小さい極大点が前記最大サンプルから見て前記最小サンプルの方向に存在する場合、前記極大点を最大サンプルとして置き換え、前記最小サンプルより微分輝度値の大きい極小点が前記最小サンプルから見て前記最大サンプルの方向に存在する場合、前記極小点を最小サンプルとして置き換えることを特徴とする。
【００２３】
本発明第８の態様によれば、最大サンプルより微分輝度値の小さい極大点が最小サンプルの方向に存在する場合、前記極大点を最大サンプルとして置き換え、最小サンプルより微分輝度値の大きい極小点が前記最大サンプルの方向に存在する場合、前記極小点を最小サンプルとして置き換えるようにしたので、計測範囲内に非計測対象の血管が含まれるような場合に、中央の計測対象の血管を対象として計測を行うことが可能となる。
【００２４】
また、本発明第９の態様では、前記血流領域端点特定手段は、
前記最大サンプル（ｉｍａｘ）より微分輝度値が小さく第１所定値（Ｈmax）より微分輝度値が大きい極大点（ｉp）から見て先頭（ｓ＋１個のサンプルのうちサンプルｉ＝０）のサンプル寄りで第１所定値（Ｈmax）より小さい微分輝度値をとるサンプルを前記血流領域の第１の端点として特定し、前記最小サンプル（ｉｍｉｎ）より微分輝度値が大きく第１所定値（Ｈmax）より小さい第２所定値（Ｈmin）より微分輝度値が小さい極小点から見て末尾（ｓ＋１個のサンプルのうちサンプルｉ＝ｓ）のサンプル寄りで第２所定値（Ｈmin）より大きい微分輝度値をとるサンプルを前記血流領域の第２の端点として特定することを特徴とする。
【００２５】
本発明第９の態様によれば、最大サンプルより微分輝度値が小さく第１所定値（Ｈmax）より微分輝度値が大きい極大点（ｉp）から見て先頭（ｉ＝０）のサンプル寄りで第１所定値（Ｈmax）より小さい微分輝度値をとるサンプルを血流領域の第１の端点として特定し、最小サンプル（ｉｍｉｎ）より微分輝度値が大きく第１所定値（Ｈmax）より小さい第２所定値（Ｈmin）より微分輝度値が小さい極小点（ｉn）から見て末尾（ｉ＝ｓ）のサンプル寄りで第２所定値（Ｈmin）より大きい微分輝度値をとるサンプルを血流領域の第２の端点として特定するようにしたので、計測範囲に含まれる血管の血流領域を的確に特定することが可能となる。
【００２６】
また、本発明第１０の態様では、前記表示手段は、前記平行な２辺の各々に特定された前記血流領域の４端点を頂点とする外形が台形の図形を、前記眼底画像に重ねて表示することを特徴とする。
【００２７】
本発明第１０の態様によれば、平行な２辺の各々に特定された血流領域の４端点を頂点とする外形が台形の図形を、眼底画像に重ねて表示するようにしたので、計測用図形が血管の血流領域を適切に指示しているかを視覚的に確認することが可能となる。
【００２８】
また、本発明第１１の態様では、コンピュータを、上記眼底画像処理装置として機能させるためのプログラムを提供する。コンピュータとしては、同一の筐体に内蔵される画像入力装置で撮影された眼底画像に対して画像処理をするものも含まれる。
【００２９】
本発明第１１の態様によれば、プログラムを組み込むことにより、眼底画像処理装置を汎用のコンピュータにより実現することが可能となる。タブレットやスマートフォンなどの携帯型コンピュータに、アプリと呼ばれるプログラムを組み込んだ場合には、内蔵されるカメラで撮影された眼底画像に対して画像処理を実現でき、眼底画像入力装置および眼底画像処理装置の双方の機能を単一のコンピュータにより実現することが可能となる。
【００３０】
また、本発明第１２の態様では、
眼底画像に含まれる血管の血流領域を特定する際に、計測用図形を表示する表示装置であって、
外部からの指示に基づき、外形が台形であり、平行な２辺の中点を結ぶ基準線を備える計測用図形を表示し、
外部からの指示に基づき配置された計測用図形の平行な２辺それぞれを計測範囲として設定し、設定された計測範囲の２辺に基づいて血流領域の２端点を各々特定し、特定された血流領域の４端点を４頂点とする外形が台形形状の図形を前記計測用図形に代えて表示することを特徴とする計測用図形表示装置を提供する。
【００３１】
本発明第１２の態様によれば、眼底画像に含まれる血管の血流領域を特定する際に、外部からの指示に基づき、外形が台形であり、平行な２辺の中点を結ぶ基準線を備える計測用図形を表示し、外部からの指示に基づき配置された計測用図形の平行な２辺それぞれを計測範囲として設定し、設定された計測範囲の２辺に基づいて血流領域の２端点を各々特定し、特定された血流領域の４端点を４頂点とする外形が台形形状の図形を計測用図形に代えて表示するようにしたので、計測用図形がいわゆる計測ゲージのように、変形して血流領域を計測したかのような印象を与えるとともに、血流領域における血流幅を直感的に把握し易くなる。
【発明の効果】
【００３２】
本発明によれば、計測者の主観や計測箇所によるバラツキの問題を解消しながら、血流幅比率を客観的に計測することが可能となる。
【図面の簡単な説明】
【００３３】
【図１】本発明の一実施形態に係る眼底画像処理装置１００のハードウェア構成図である。
【図２】本発明の一実施形態に係る眼底画像処理装置の構成を示す機能ブロック図である。
【図３】本発明の一実施形態に係る眼底画像処理装置の処理概要を示すフローチャートである。
【図４】眼底画像と解析範囲を示す図である。
【図５】計測用図形の初期状態を示す図である。
【図６】計測用図形の操作による変形の様子を示す図である。
【図７】解析範囲、計測用図形、計測範囲設定画面を示す図である。
【図８】ステップＳ４００における血流領域端点の特定の処理を示すフローチャートである。
【図９】計測処理前と処理後における計測用図形の変化を示す図である。
【図１０】始点辺ＳＨまたは終点辺ＥＨと１．５画素分移動させた画素の関係を示す図である。
【図１１】計測範囲の眼底画像と輝度値、微分輝度値の関係を示す図である。
【図１２】両側の血管の一部が計測範囲に含まれる場合の輝度値、微分輝度値の関係を示す図である。
【図１３】血管が一本だけ計測範囲に含まれる場合の輝度値、微分輝度値の関係を示す図である。
【図１４】実際の画像に計測用図形を重ねて表示させた例を示す図である。
【図１５】解析範囲内に６箇所（動脈３箇所、静脈３箇所）の計測用図形が配置された例を示す図である。
【図１６】解析範囲内に６箇所（動脈３箇所、静脈３箇所）の計測用図形が配置された例を示す図である。
【図１７】眼底画像と解析範囲（動静脈交叉部）を示す図である。
【図１８】解析範囲、計測用図形（動静脈交叉部）、計測範囲設定画面を示す図である。
【図１９】解析範囲の静脈上に、処理後の計測用図形を重ねて表示させた例を示す図である。
【図２０】解析範囲内に２箇所（静脈２箇所）の計測用図形が配置された例を示す図である。
【図２１】解析範囲内に２箇所（静脈２箇所）の計測用図形が配置された例を示す図である。
【発明を実施するための形態】
【００３４】
以下、本発明の好適な実施形態について図面を参照して詳細に説明する。
＜１．装置構成＞
図１は、本発明の一実施形態に係る眼底画像処理装置のハードウェア構成図である。本実施形態に係る眼底画像処理装置１００は、タブレットやスマートフォンなどカメラや無線通信機能を備えた携帯型の汎用のコンピュータで実現することができる。本実施形態に係る眼底画像処理装置１００は、図１に示すように、ＣＰＵ（Central Processing Unit）１と、コンピュータのメインメモリであるＲＡＭ（Random Access Memory）２と、ＣＰＵ１が実行するプログラムやデータを記憶するためのフラッシュメモリ等の大容量の記憶装置３と、解析対象の血流領域を指示したり操作メニューを指示したりするためのタッチパネル入力Ｉ／Ｆ（インターフェース）４と、無線通信等により外部の眼底カメラから画像入力したり、解析結果を外部のサーバーコンピュータやプリンターに出力するためのデータ入出力Ｉ／Ｆ（インターフェース）５と、画像データ等を表示するための液晶ディスプレイ、有機ＥＬディスプレイ等の表示デバイスである表示部６と本体内蔵のカメラ（眼底部を撮影するための光源とマクロレンズを装着して眼底カメラとしても使用できる）である画像入力部７と、を備え、互いにバスを介して接続されている。もちろん、同等の機能をもつものであれば、据置型の汎用のコンピュータで実現することもできる。
【００３５】
図２は、本実施形態に係る眼底画像処理装置の構成を示す機能ブロック図である。図２において、１０は指示入力手段、２０は表示手段、３０は画像拡大手段、４０は計測範囲設定手段、５０は血流領域端点特定手段、６０は動静脈径比算出手段、７０は静脈径比算出手段、８０は眼底画像記憶手段、９０は情報記憶手段である。
【００３６】
指示入力手段１０は、外部からの指示入力を各手段に伝えるものであり、タッチパネル入力Ｉ／Ｆ（インターフェース）４とＣＰＵ１により実現される。表示手段２０は、記憶手段に記憶された情報や各手段により処理された情報を表示するためのものであり、表示部６とＣＰＵ１により実現される。画像拡大手段３０は、眼底画像記憶手段８０から読み込まれた眼底画像のうち設定された解析範囲をトリミングし、所定の画素数になるように拡大処理を行う。これは、血流幅の計測が画素単位に行われるが、原画の画素数のままでは値が小さ過ぎて、２種の血流幅どうしで割り算を行うと、丸め誤差により高血圧や動脈硬化の度合を判定するのに十分な精度が得られないためである。計測範囲設定手段４０は、眼底画像に対して、動脈または静脈を含む計測用図形に基づき、当該計測用図形の平行な２辺の少なくとも一方の辺を計測範囲として設定する。血流領域端点特定手段５０は、計測範囲から血流に相当する領域である血流領域の２端点を特定する。尚、単一の計測用図形から平行な２辺に基づき１次元の血流領域である計測範囲は同時に最大２箇所特定される。
【００３７】
動静脈径比算出手段６０は、動脈に対応する計測範囲について算出された血流領域の２端点間の距離と、静脈に対応する計測範囲について算出された血流領域の２端点間の距離との比率に基づいて、動静脈径比を算出する。実際には、動脈に対応する計測範囲および静脈に対応する計測範囲からは、計測用図形の平行な２辺の各々が計測範囲に対応して、各々２箇所の血流領域が特定され２つの距離が算出されるので、各々２つの距離の平均値どうしの比率に基づいて、動静脈径比を算出する。ここで、１次元の血流領域の２端点を算出するのに、２次元の計測用図形を用いて１次元の計測範囲を設定する理由は、血流領域が計測対象の血管の走行方向と直交するようにさせるためで、１次元の直線で画面上で指定することにより計測範囲を設定すると、血流領域が計測対象の血管の走行方向と斜めになりがちになるためである。また、後述の静脈径比を正確に計測するためには、同一静脈の血管で２箇所の平行な１次元の血流領域を特定することが求められるが、血管走行方向と直交し平行な２辺を備える計測用図形に基づき容易かつ高精度に計測できる。静脈径比算出手段７０は、単一の計測用図形の一方の辺に基づいてについて算出された血流領域の２点間の距離と、他方の辺に基づいて算出された血流領域の２点間の距離との比率に基づいて、静脈径比を算出する。「動静脈径比」「静脈径比」はいずれも血流幅比率である。なお、本発明において「血流」領域を特定し、「血流」幅比率を算出しているのは、血管の幅が透明であることから、血管壁の特定、血管幅比率の算出が難しいためである。しかし、動脈硬化により一般に血管壁は厚くなり、血流幅は狭くなるが、血管の幅はあまり変化しない。そのため、むしろ血管幅比率を算出するよりも血流領域を特定し、血流幅比率を算出する方法の方が、高血圧や動脈硬化の度合を判定するのには有用である。
【００３８】
画像拡大手段３０、計測範囲設定手段４０、血流領域端点特定手段５０、動静脈径比算出手段６０、静脈径比算出手段７０は、ＣＰＵ１が、記憶装置３に記憶されているプログラムを実行することにより実現される。眼底画像記憶手段８０は、データ入出力Ｉ／Ｆ５または画像入力部７により可視光または赤外線・光源方式の眼底カメラを用いてフルカラー（可視光の場合）・モノクロ（赤外光の場合）で撮影された、血流領域の抽出対象となるフルカラー眼底画像を記憶した記憶手段であり、記憶装置３により実現される。尚、赤外線光源方式の眼底カメラで撮影された眼底画像は、通常モノクロで撮影されるが、ＲＧＢに同一の値をもたせたフルカラー画像として記憶装置３には記憶される。また、赤外線光源方式の眼底カメラは、コンピュータとは外付けになり、データ入力Ｉ／Ｆ５を介して記憶装置３に取り込まれる形態になるが、可視光光源方式の眼底カメラは、コンピュータに内蔵する画像入力部７に眼底撮影用の光源とマクロレンズを装着して眼底カメラとして使用することもできる。眼底画像処理装置にフルカラーの眼底画像を読み込ませて、そのまま処理を行う場合は、ＲＡＭ２が眼底画像記憶手段８０としての役割を果たす。フルカラーの眼底画像とは、ＲＧＢの３成分により記録された画像データであり、被験者の眼底が撮影されたものである。本実施形態では、ＲＧＢ各色８ビット２５６階調で記録されたものを用いている。情報記憶手段９０は、血流領域端点特定手段５０により特定された血流領域端点の情報、動静脈径比算出手段６０、静脈径比算出手段７０により算出された血流幅、血流幅比率等を記憶する記憶手段であり、記憶装置３により実現される。
【００３９】
図２に示した各構成手段は、現実には図１に示したように、コンピュータおよびその周辺機器等のハードウェアに専用のプログラムを搭載することにより実現される。すなわち、コンピュータが、専用のプログラムに従って各手段の内容を実行することになる。コンピュータに内蔵するカメラに眼底撮影用の光源やマクロレンズを装着すれば眼底カメラの機能を代替でき（現状、日本国内では認可されていないが）、眼底撮影から眼底画像処理まで一貫して単一の携帯型コンピュータで実現できる。なお、本明細書において、コンピュータとは、ＣＰＵ等の演算処理部を有し、データ処理が可能な装置を意味し、カメラや無線通信機能を備えたタブレット・スマートフォン等の携帯型端末だけでなく、パーソナルコンピュータ等の汎用のコンピュータも含む。
【００４０】
図１に示した記憶装置３には、ＣＰＵ１を動作させ、コンピュータを、眼底画像処理装置として機能させるための専用のプログラムが実装されている。この専用のプログラムを実行することにより、ＣＰＵ１は、画像拡大手段３０、計測範囲設定手段４０、血流領域端点特定手段５０、動静脈径比算出手段６０、静脈径比算出手段７０としての機能、および指示入力手段１０、表示手段２０の一部としての機能を実現することになる。また、記憶装置３は、眼底画像記憶手段８０、情報記憶手段９０として機能するだけでなく、眼底画像処理装置としての処理に必要な様々なデータを記憶する。
【００４１】
＜２．処理動作＞
＜２．１．前処理＞
まず、処理対象とする眼底画像を撮影または用意する。血管計測を目的とした眼底画像としては、赤外線光源により撮影されるグレースケール（モノクロ）画像（眼底血管造影像）であることが理想であるが、撮影機器が高額なため、本実施形態では普及している可視光光源により撮影されるカラー画像を処理対象とする。グレースケール（モノクロ）画像の場合もフルカラー画像の形式に統一される。可視光光源で撮影されるカラー画像も、後述するように後にグレースケール画像に変換するが、色の特徴を活用した方が、血流領域を高精度に抽出できるためである。カラー画像としては、Ｒ，Ｇ，Ｂ各成分少なくとも８ビット２５６階調以上のフルカラー画像が好ましい。フルカラーの眼底画像としては、デジタル方式の眼底カメラによりフルカラーで撮影した画像ファイルがあれば、そのまま使用できる。アナログ方式の眼底カメラにより写真媒体に記録された過去のものであれば、保管されていたアナログのモノクロまたはカラーのネガ・ポジフィルム、印画紙、インスタント写真等をスキャナによりフルカラーで読み取る等してデジタルの眼底画像ファイルを取得する。この時、たとえ原画がモノクロであっても、フルカラーで読み取る。最近では、可視光・光源方式のデジタル眼底カメラを用いてフルカラーで撮影することにより眼底画像がデジタルファイルの形式で得られる。取得した眼底画像は、眼底画像処理装置の眼底画像記憶手段８０に記憶させる。本実施形態では、眼底画像としてＲ，Ｇ，Ｂ各成分８ビット２５６階調のフルカラー画像を用意する。
【００４２】
＜２．２．処理概要＞
次に、図１、図２に示した眼底画像処理装置の処理動作について説明する。図３は、本発明の一実施形態に係る眼底画像処理装置の処理概要を示すフローチャートである。まず、画像拡大手段３０が、解析範囲の設定を行う（ステップＳ１００）。具体的には、まず、眼底画像を表示部６に表示させ、利用者に対して解析範囲の指定を促す。利用者が、指示入力手段１０であるタッチパネル入力Ｉ／Ｆ４を介して眼底画像内の所定の範囲を指定すると、画像拡大手段３０は、指定された範囲を解析範囲として設定する。所定の範囲の指定は、例えば、利用者がタッチパネル等を用いて、眼底画像上に表示される長方形に対してドラッグ操作等を行って解析範囲を特定することにより行うことができる。このときの眼底画像の表示例を図４に示す。図４の例では、横１２８０画素×縦１０２４画素の眼底画像に、横１０３画素×縦１０３画素の解析範囲が設定された状態を示している。解析範囲は、後に設定される複数の計測範囲全てを含む範囲として設定される。すなわち、計測範囲は、設定された解析範囲の中で設定されることになる。
【００４３】
解析範囲が設定されたら、画像拡大手段３０が、設定された解析範囲をトリミングし、トリミングした部分の拡大処理を行う（ステップＳ２００）。具体的には、解析範囲の眼底画像の画素数を増加させる。この時、見かけ上は解像度が上がるが、拡大処理に特殊な補間を加えて解像度そのものを向上させる必要はない。処理対象とする画像の画素数は、血流幅が数十画素程度になるように設定することが好ましい。横１２８０画素×縦１０２４画素の画像に眼底全体が収まるような場合には、血流幅が数画素程度となり、そのままでは、血流幅どうしで割り算をすると、的確な解析を行うことができない。そのため、元の眼底画像が横１２８０画素×縦１０２４画素である場合、縦横ともに３００％〜５００％拡大する。ここでは、事前に拡大率を５００％に設定している。本実施形態では解像度を上げる必要はないので、拡大処理を行うためには、増加させた画素に適当な画素値を与えればよく、高度な補間手法を用いる必要はない。そのため、本実施形態では、公知のバイリニアー法を用いている。
【００４４】
解析範囲のトリミングと拡大処理が行われたら、計測範囲の設定を行う（ステップＳ３００）。この際、まず、計測範囲設定手段４０は、利用者に対して計測範囲の指定を促す。具体的には、拡大した解析範囲の画像に重ねて、外形が長方形（正方形も含む意味）状の計測用図形を配置した表示を行う。この計測用図形は、利用者から見ると、計測ゲージのような印象を与える図形である。計測用図形が表示されたら、利用者は、指示入力手段１０により計測用図形を縦横に変倍、移動、回転させて、平行な２辺が計測対象の血管幅をまたがるように配置する。
【００４５】
図５に表示部６に表示される計測用図形を示す。計測用図形は、図５に示すように頂点１１、頂点１２、頂点２１、頂点２２の４頂点からなる長方形となっている。頂点１１と頂点１２を結ぶ線分の中点を始点、頂点２１と頂点２２を結ぶ線分の中点を終点と呼ぶことにし、始点と中点を結ぶ線分を基準線ＫＨとする。また、始点を中点とする線分を始点辺ＳＨ、終点を中点とする線分を終点辺ＥＨとする。
【００４６】
計測用図形の６点のいずれかを、指示入力手段１０により選択すると、その点が移動可能な状態になり、指定位置で再度操作指示を行うと固定化される。計測用図形の上を選択すると、計測用図形全体が移動可能な状態になる。計測用図形の始点または終点を選択した場合は、長さの伸縮と回転が可能となる。図６（ａ）の例では、○で囲った終点を選択した場合を示している。終点を選択した場合、図６（ａ）の上図に示すように、右方向への操作指示により、右方向に計測用図形全体の長さの伸縮を行う。また、終点を選択した場合、図６（ａ）の下図に示すように、上下方向への操作指示により、始点を中心として計測用図形全体の回転を行う。図６（ａ）において、×を付した点は固定されて動かない点を示している。図６（ａ）は、終点を選択した場合を示しているが、始点を選択した場合も同様に、伸縮・回転が可能である。
【００４７】
計測用図形の頂点１１または頂点１２を選択した場合、すなわち、始点辺ＳＨの両端のいずれかを選択した場合は、図６（ｂ）に示すように、上下方向への操作指示により、幅の一括伸縮を行う。計測用図形の頂点２１または頂点２２を選択した場合、すなわち、終点辺ＥＨの両端のいずれかを選択した場合は、図６（ｃ）に示すように、上下方向への操作指示により、終点辺ＥＨのみの幅の伸縮を行う。図６（ｂ）（ｃ）においても、×を付した点は固定されて動かない点を示している。この操作は、計測用図形を手動で血管幅に合わせて計測できるようにするために設けた機能で、計測結果に操作者の主観が入るため、本発明では通常使用しない。
【００４８】
図５に示したような初期状態のまま、はじめに、利用者は、計測用図形全体を計測対象の血管の位置に移動させる。そして、基準線ＫＨを計測対象の血管走行の中心線（血管幅の中心を走行方向に結んだ仮想的な曲線）に概ね合わせるように、計測用図形全体の回転を行う。この時、計測用図形に含まれる中心線がほぼ直線になり、基準線ＫＨとほぼ一致するように、計測用図形全体の長さの伸縮も併せて行う。最後に、頂点１１と頂点２１を結ぶ辺と頂点１２と頂点２２を結ぶ辺が、計測対象の血管の外側に十分はみ出るように、幅の一括伸縮を行う。この時、計測用図形の内部には計測対象の血管領域のみが含まれていることが望ましいが、計測用図形の内部に非計測対象の血管領域が多少含まれていても計測に支障ない。そうすると、始点辺ＳＨと終点辺ＥＨが血管とほぼ直角に交差する位置に配置される。眼底画像内の解析範囲の所定の位置に計測用図形を配置すると、計測範囲設定手段４０は、指定された計測用図形に基づき、計測用図形の平行な２辺それぞれを１次元の計測範囲として設定する。このときの解析範囲と計測用図形の表示例を図７（ａ）に示す。図７（ａ）の例では、解析範囲である横５１５画素×縦５１５画素の眼底画像に、外形が長方形状の計測用図形が配置された状態を示している。後述するように、計測は、始点辺ＳＨ、終点辺ＥＨの２辺上およびその近傍における画素を用いて行われる。したがって、２次元の計測用図形を配置することによって、計測用図形の始点辺ＳＨ、終点辺ＥＨが１次元の計測範囲として設定されることになる。１次元の計測範囲を設定するために、画面上で１次元の線分で指定すると、血管と直交する計測範囲を適切に設定できない。これに対して、対向する辺が平行である２次元の計測用図形を画面上で配置して指定することにより、血管と直交する計測範囲を適切に設定することができる。本実施形態では、１つの計測用図形の平行な２辺をそれぞれ計測範囲として設定するが、平行な２辺のうち、動静脈径比の計測の場合には、どちらか１辺のみを計測範囲として設定することも可能である。動静脈径比の計測の場合には、１辺だけであっても計測は可能であるが、血管の幅や傾きが常に一定であるわけでなないことから、平行な２辺を計測範囲とする方がより正確な計測を行うことができる。
【００４９】
計測範囲が設定されたら、計測範囲設定手段４０は、計測範囲設定画面を表示させ、利用者に対して計測範囲が動脈か静脈かのいずれかの選択指示を促す。図７（ｂ）に計測範囲設定画面の一例を示す。計測範囲設定手段４０は、利用者に対して複数の計測用図形に基づく計測範囲の指定、配置された各計測用図形に対して動脈か静脈かのいずれかの選択指示を促し、計測範囲を設定する。したがって、１つの計測用図形に対応する２つの計測範囲については、同一の設定が行われる。動静脈径比を求める場合は、計測範囲の設定として、少なくとも動脈、静脈それぞれに対応する箇所に１つ以上、トータルで２つ以上計測用図形を配置して設定されるように行う。この時、動脈と静脈に配置される計測用図形は互いに近傍であることが望ましいが、個数は同数である必要は無い。動静脈交叉部・静脈径比を求める場合は、計測範囲の設定として、上側に動脈が交叉している動静脈交叉部の静脈に対応する箇所に１つ以上設定されるように行う。この時、始点辺ＳＨは動静脈交叉部近辺の静脈狭窄部に、終点辺ＥＨは狭窄が起こっていない動静脈交叉部より離れた同一の静脈上の箇所に位置するように計測用図形を伸縮させる（始点辺ＳＨと終点辺ＥＨの位置関係は逆でも良い。）この時、計測用図形の内部には計測対象の静脈領域のみが含まれていることが望ましいが、計測用図形の内部に交叉している動脈領域など非計測対象の血管領域が多少含まれていても計測に支障ない。
【００５０】
図７（ｂ）に示した計測範囲設定画面では、計測範囲が動脈か静脈かのいずれかの選択指示だけでなく、計測範囲の解析条件について様々な指定を行うことができるようになっているが、動脈と静脈のいずれであるかの指定は必須となっている。
【００５１】
計測範囲の設定が行われたら、次に、血流領域端点特定手段５０が、計測範囲設定画面を利用して入力された内容に従って、血流領域端点の特定を行う（ステップＳ４００）。具体的には、計測範囲である始点辺ＳＨ、終点辺ＥＨを用いて血流領域端点の特定を行う。
【００５２】
図８は、ステップＳ４００における血流領域端点の特定の処理を示すフローチャートである。まず、血流領域端点特定手段５０は、計測範囲である始点辺ＳＨ、終点辺ＥＨの２辺に対応する計測用図形の情報を取得する（ステップＳ４１０）。具体的には、まず、計測用図形の始点の座標（Ｘｏ１，Ｙｏ１）と終点の座標（Ｘｏ２，Ｙｏ２）、始点辺ＳＨの頂点（Ｘｏ１１，Ｙｏ１１）と（Ｘｏ１２，Ｙｏ１２）、終点辺ＥＨの頂点（Ｘｏ２１，Ｙｏ２１）と（Ｘｏ２２，Ｙｏ２２）を取得するとともに、始点辺ＳＨまたは終点辺ＥＨのいずれか大きい方の長さの１／２を計測用図形の幅Ｗとして取得する。図９（ａ）は、配置された計測用図形と各点の座標値を示す図である。図９（ａ）に示すように、計測用図形は、長方形内部に始点辺ＳＨの中点と終点辺ＥＨの中点を結ぶ線分を備えた形状となっている。長方形内部の線分は基準線ＫＨであり、利用者が、眼底画像に計測用図形を配置する際の目安となるもので、実際の計測処理には使用しない。始点の座標（Ｘｏ１，Ｙｏ１）と終点の座標（Ｘｏ２，Ｙｏ２）は、眼底画像の画素の位置を特定するものである。
【００５３】
上述のように始点側の２頂点（Ｘｏ１１，Ｙｏ１１）、（Ｘｏ１２，Ｙｏ１２）を結んで得られる線分が始点辺ＳＨ、終点側の２頂点（Ｘｏ２１，Ｙｏ２１）、（Ｘｏ２２，Ｙｏ２２）を結んで得られる線分が終点辺ＥＨである。計測用図形が長方形でなく一組の対向する辺のみが平行である台形に変形した状態で配置されて、始点辺ＳＨと終点辺ＥＨの長さが異なる場合は、長い方の辺の長さを用いて伸ばす処理を行う。すなわち、計測範囲でもある短い方の辺の両頂点の位置が中点から同距離ずつ離れるようにして再設定される。この結果、始点側の２頂点（Ｘｏ１１，Ｙｏ１１）、（Ｘｏ１２，Ｙｏ１２）間の距離と、終点側の２頂点（Ｘｏ２１，Ｙｏ２１）、（Ｘｏ２２，Ｙｏ２２）間の距離は同一となる。そして、以下のステップＳ４２０〜ステップＳ４６０の処理を始点辺ＳＨ、終点辺ＥＨそれぞれについて実行する。
【００５４】
はじめに、輝度配列を作成する（ステップＳ４２０）。具体的には、まず、
血管と交差して配置された辺（始点辺ＳＨまたは終点辺ＥＨ）をサンプリングする。すなわち、血管と交差して配置された辺（始点辺ＳＨまたは終点辺ＥＨ）の一方の頂点の座標を（Ｘｓ，Ｙｓ）、辺の他方の頂点の座標を（Ｘｅ，Ｙｅ）とし、各辺に所定の間隔でｉ＝０，・・・，ｓのｓ＋１個のサンプル（Ｘ（ｉ），Ｙ（ｉ））を設定する。具体的には、処理対象の辺が始点辺ＳＨの場合、Ｘｓ＝Ｘｏ１１、Ｙｓ＝Ｙｏ１１、Ｘｅ＝Ｘｏ１２、Ｙｅ＝Ｙｏ１２となり、処理対象の辺が終点辺ＥＨの場合、Ｘｓ＝Ｘｏ２１、Ｙｓ＝Ｙｏ２１、Ｘｅ＝Ｘｏ２２、Ｙｅ＝Ｙｏ２２となる。各サンプルの座標（Ｘ（ｉ），Ｙ（ｉ））は、Ｄｘ＝Ｘｅ−Ｘｓ、Ｄｙ＝Ｙｅ−Ｙｓとして、Ｘ（ｉ）＝Ｄｘ・ｉ／ｓ＋Ｘｓ、Ｙ（ｉ）＝Ｄｙ・ｉ／ｓ＋Ｙｓで定義される。
【００５５】
これにより、１次元の座標配列（Ｘ（ｉ），Ｙ（ｉ））（ｉ＝０，・・・，ｓ）が作成される。続いて、各サンプルに対応する画素の輝度値を取得し、ｓ＋１個の輝度値からなる輝度配列Ｇ（ｉ）（ｉ＝０，・・・，ｓ）を得る。この際、より正確性を高めるため、辺に対して直交方向の両側に各サンプルから１．５画素分移動させた点を設定する。１．５画素分とするのは、辺が最大４５°傾いた場合であっても、１画素以上異なる画素の輝度値を取得するためである。移動させる長さは、１．５画素に限らず、１．５画素以上であれば、適宜設定することができる。図１０は、始点辺ＳＨまたは終点辺ＥＨと１．５画素分移動させた画素の関係を示す図である。
【００５６】
図１０において、中央の線分は、始点辺ＳＨまたは終点辺ＥＨであり、（Ｘｓ，Ｙｓ）から（Ｘｅ，Ｙｅ）までｓ＋１個のサンプル（Ｘ（ｉ），Ｙ（ｉ））が設定されている。また、破線で囲った網掛け部分は、血流領域を示している。そして、各サンプル（Ｘ（ｉ），Ｙ（ｉ））について、始点辺ＳＨまたは終点辺ＥＨに直交する方向に１．５画素分移動させた位置に点（Ｘｌ，Ｙｌ）、（Ｘｒ，Ｙｒ）を設定する。各サンプル（Ｘ（ｉ），Ｙ（ｉ））について、点（Ｘｌ，Ｙｌ）、（Ｘｒ，Ｙｒ）を設定することにより、始点辺ＳＨまたは終点辺ＥＨから１．５画素分移動した２本の線分上の点の輝度値が始点辺ＳＨまたは終点辺ＥＨのサンプルの輝度値に反映される。
【００５７】
このように、始点辺ＳＨまたは終点辺ＥＨから１．５画素分だけ移動させたサンプルの輝度値を反映させるのは、始点辺ＳＨまたは終点辺ＥＨだけであると、始点辺ＳＨまたは終点辺ＥＨのサンプルにスポット的な輝度ムラがあると、その影響により、正しい測定ができないためである。このような点（Ｘｌ，Ｙｌ）、（Ｘｒ，Ｙｒ）は、以下の〔数式１〕に従った処理を実行することにより算出する。
【００５８】
〔数式１〕
Ｘｌ＝−Ｄｙ・１．５＋Ｘ（ｉ）
Ｙｌ＝Ｄｘ・１．５＋Ｙ（ｉ）
Ｘｒ＝Ｄｙ・１．５＋Ｘ（ｉ）
Ｙｒ＝−Ｄｘ・１．５＋Ｙ（ｉ）
【００５９】
そして、（Ｘ（ｉ），Ｙ（ｉ））（Ｘｌ，Ｙｌ）、（Ｘｒ，Ｙｒ）の輝度値を平均化したものを、各サンプル（Ｘ（ｉ），Ｙ（ｉ））の輝度値とする。すなわち、特定された２点における輝度値を用いて各サンプル（Ｘ（ｉ），Ｙ（ｉ））に対応する輝度値を補正する。眼底画像がフルカラー画像である場合、輝度値を取得するためにグレースケールに変換する必要がある。本実施形態では、以下のようにして、対象の画素をグレースケール変換する。
【００６０】
フルカラーの眼底画像は、各色８ビット２５６階調の画像データである。したがって、ｘ方向の画素数Ｘｇ、ｙ方向の画素数Ｙｇの眼底画像は、色成分を示す変数ｃ＝０（Ｒｅｄ），１（Ｇｒｅｅｎ），２（Ｂｌｕｅ）とすると、Ｉｍａｇｅ（ｘ，ｙ，ｃ）＝０〜２５５（ｘ＝０，・・・，Ｘｇ−１；ｙ＝０，・・・，Ｙｇ−１；ｃ＝０，１，２）と定義される。
【００６１】
グレースケールの画像への変換の際、フルカラーの眼底画像Ｉｍａｇｅ（ｘ，ｙ，ｃ）に対してＧｒｅｅｎ成分のみを使用する、即ち、Ｇｒ（ｘ，ｙ）＝Ｉｍａｇｅ（ｘ，ｙ，１）により変換する方法が最良であることが知られている。理由は、眼底カメラで使用されている３原色フィルタではＧｒｅｅｎ成分に比べＲｅｄとＢｌｕｅ成分の解像度が低く、これらの分色画像では血管の輪郭が不鮮明になりがちなためである。しかし、本発明者は、多くの眼底カメラ画像を分析した結果、Ｇｒｅｅｎ成分は照明光が十分にあたっている中央付近では鮮明であるが、照明光が弱い周辺部ではＧｒｅｅｎ成分よりむしろＢｌｕｅ成分の方が鮮明であり、更にＢｌｕｅ成分はＧｒｅｅｎ成分と異なり、血管部の輝度が周囲に比べて高くなるという反転現象を見出した（眼底画像全般に血管部の輝度は周囲に比べて低くなる）。そこで、本実施形態では、以下の〔数式２〕に従った処理を実行することによりグレースケール変換する手法を提案する。ただし、〔数式２〕は全てのフルカラーの眼底画像に適用になるわけではなく、従来から言われているＧｒ（ｘ，ｙ）＝Ｉｍａｇｅ（ｘ，ｙ，１）により変換する方法の方が適切な眼底画像も少なからず存在するため、これらは与えられた眼底画像ごとに適宜選択する必要がある。血流領域端点特定手段５０は、以下の〔数式２〕に従った処理を実行することによりグレースケール変換する。
【００６２】
〔数式２〕
Ｇｒ（ｘ，ｙ）＝[Ｉｍａｇｅ（ｘ，ｙ，１）・（２５５−Ｉｍａｇｅ（ｘ，ｙ，２））]1/2
【００６３】
上記〔数式２〕において、Ｉｍａｇｅ（ｘ，ｙ，１）は、フルカラーの眼底画像のうち、Ｇ（グリーン、緑）の成分を示し、Ｉｍａｇｅ（ｘ，ｙ，２）は、フルカラーの眼底画像のうち、Ｂ（ブルー、青）の成分を示している。したがって、上記〔数式２〕においては、各画素のＧ成分と、Ｂ成分をネガポジ反転したものの相乗平均値（積の平方根）を求め、グレースケール形式の画像Ｇｒ（ｘ，ｙ）を得ている。
【００６４】
さらに、血流領域端点特定手段５０は、グレースケール形式の画像Ｇｒ（ｘ，ｙ）に対して、以下の〔数式３〕に従った処理を実行することによりネガポジ反転してグレースケールの眼底画像Ｇｒａｙ（ｘ，ｙ）を得る。
【００６５】
〔数式３〕
Ｇｒａｙ（ｘ，ｙ）＝２５５−Ｇｒ（ｘ，ｙ）
【００６６】
上記〔数式３〕に従った処理を実行して、ネガポジ反転することにより、周囲に比べ輝度が低い血流領域の画素値が高い状態に変換される。ここで、ネガポジ反転とは、元の画素値（階調値）の大きいもの程小さい値に変換する処理を意味する。
【００６７】
以上のようにして、グレースケール化により（Ｘ（ｉ），Ｙ（ｉ））（Ｘｌ，Ｙｌ）、（Ｘｒ，Ｙｒ）の輝度値が得られたら、以下の〔数式４〕に従った処理を実行することにより輝度配列Ｇ（ｉ）を算出する。
【００６８】
〔数式４〕
Ｇ（ｉ）＝[Ｇｒａｙ（Ｘ（ｉ），Ｙ（ｉ））＋Ｇｒａｙ（Ｘｌ，Ｙｌ）＋Ｇｒａｙ（Ｘｒ，Ｙｒ）]／３
【００６９】
上述のように、各サンプルは、計算の便宜上、始点辺ＳＨまたは終点辺ＥＨの両端の間にｓ＋１個設定したものである。また、始点辺ＳＨまたは終点辺ＥＨは、画素の配列であるｘｙ軸とは傾いて設定されることもあるが、サンプルの座標値は整数値に近似されるため、そのまま画素の輝度値を取得できる。
【００７０】
各サンプルｉについて輝度値を取得して、輝度配列Ｇ（ｉ）が得られたら、次に、輝度配列Ｇ（ｉ）を用いて微分輝度値の配列に変換した微分配列Ｄ（ｉ）を算出する（ステップＳ４３０）。具体的には、以下の〔数式５〕に従った処理を実行することにより微分配列Ｄ（ｉ）を算出する。ただし、ｉ＜３またはｉ＞ｓ−２の場合、Ｄ（ｉ）＝０とする。
【００７１】
〔数式５〕
Ｄ（ｉ）＝Ｇ（ｉ）＋Ｇ（ｉ＋１）＋Ｇ（ｉ＋２）−Ｇ（ｉ−１）−Ｇ（ｉ−２）−Ｇ（ｉ−３）
【００７２】
上記〔数式５〕では、サンプルｉに隣接するサンプルｉ＋１、サンプルｉ−１を含み、サンプルｉ，ｉ＋１，ｉ＋２，ｉ−１，ｉ−２，ｉ−３の計６個のサンプルを用いてサンプルｉにおける微分輝度値Ｄ（ｉ）を算出しているが、より多くのサンプルを用いて微分輝度値Ｄ（ｉ）を算出するようにしてもよい。また、上記〔数式５〕においては、隣接するサンプルとの加算値＋Ｇ（ｉ＋１）や、差分値−Ｇ（ｉ−１）を用いている。
【００７３】
微分配列Ｄ（ｉ）が得られたら、次に、微分配列Ｄ（ｉ）の最大値、最小値を特定する（ステップＳ４４０）。具体的には、以下の〔数式６〕に従った処理を実行することにより微分配列Ｄ（ｉ）の最大値Ｄmaxおよび最小値Ｄminを算出する。
【００７４】
〔数式６〕
Ｄmax＝ＭＡＸi=0,S/2Ｄ（ｉ）＝Ｄ（ｉmax）
Ｄmin＝ＭＩＮi=S/2,SＤ（ｉ）＝Ｄ（ｉmin）
【００７５】
上記〔数式６〕において、第１式の“ＭＡＸi=0,S/2”は、添え字のｉ＝０からｉ＝ｓ／２までの各サンプルｉにおけるＤ（ｉ）の最大値”を求めることを示している。また、上記〔数式６〕において、第２式の“ＭＩＮi=S/2,S”は、添え字のｉ＝ｓ／２からｉ＝ｓまでの各サンプルｉにおけるＤ（ｉ）の最小値”を求めることを示している。最大値Ｄmaxをとるサンプルｉｍａｘが最大サンプル、最小値Ｄminをとるサンプルｉｍｉｎが最小サンプルとして特定される。
【００７６】
ここで、計測範囲における眼底画像と輝度値、微分輝度値の関係を図１１に示す。図１１において、横軸は、サンプルｉを示している。図１１の上段は、計測範囲近傍の眼底画像であり、網掛けした部分は、血流領域を示している。すなわち、図１１の例では、３本の血管を含む範囲に計測範囲が設定された場合を示している。図１１の中段は、輝度値Ｇ（ｉ）であり、血流領域の中央ほど輝度値が高いことを示している。図１１の下段は、輝度値Ｇ（ｉ）の微分値である微分輝度値Ｄ（ｉ）である。ステップＳ４４０においては、微分輝度値の最大値Ｄmaxおよび最小値Ｄminが算出される。図１１のように計測範囲が設定された場合、中央の血流領域の幅を計測したいのが、利用者の意図であると考えられる。ところが、現時点で得られた微分輝度値の最大値Ｄmax、最小値Ｄminは、両側の血流領域に位置している。そのため、本実施形態では、さらなる探索が行われることになる。
【００７７】
ステップＳ４４０において微分配列Ｄ（ｉ）の最大値Ｄmaxおよび最小値Ｄminが得られたら、次に、微分配列Ｄ（ｉ）の最大値Ｄmaxおよび最小値Ｄminよりもサンプル中央値ｓ／２寄りの極大値Ｄmax2、極小値Ｄmin2を特定する（ステップＳ４５０）。具体的には、以下の〔数式７〕に従った処理を実行することにより極大値Ｄmax2、極小値Ｄmin2を探索する。
【００７８】
〔数式７〕
Ｄmax2＝Ｄ（ｉmax2）＞Ｄmax／２、かつ０＜ｉ＜ｉmax2およびｉmax2＜ｉ＜ｓ／２でＤ（ｉ）＜Ｄmax2を満たすｉmax2が存在すれば、Ｄmax＝Ｄmax2とする。
Ｄmin2＝Ｄ（ｉmin2）＜Ｄmin／２、かつ０＜ｉ＜ｉmin2およびｉmin2＜ｉ＜ｓ／２でＤ（ｉ）＞Ｄmin2を満たすｉmin2が存在すれば、Ｄmin＝Ｄmin2とする。
【００７９】
極大値とは、両側のサンプルの値がその極大値より小さくなる値であり、極小値とは、両側のサンプルの値がその極小値より大きくなる値である。サンプル中央値ｓ／２寄りとは、計測範囲における中点側を意味する。ステップＳ４５０の処理により、図１１に示すように、計測範囲に複数の血管が含まれている場合に、対象とする血流領域の微分輝度値の最大値、微分輝度値の最小値以外に、サンプル中央値ｓ／２寄りに極大値Ｄmax2または極小値Ｄmin2が見つかることがある（いずれか一方しか見つからない場合や双方とも見つからない場合もある）。〔数式７〕で指定の条件を満たす極大値Ｄmax2または極小値Ｄmin2が見つかった場合、〔数式７〕では、見つかった極大値または極小値を最大値Ｄmaxまたは最小値Ｄminに置き換えている。すなわち、ステップＳ４５０においては、血流領域端点特定手段５０は、最大サンプルより微分輝度値の小さい極大点が最大サンプルから見て最小サンプルの方向に存在する場合、極大点を最大サンプルとして置き換え、最小サンプルより微分輝度値の大きい極小点が最小サンプルから見て最大サンプルの方向に存在する場合、極小点を最小サンプルとして置き換える処理を行っている。
【００８０】
対象とする血流領域の微分輝度値の最大値および最小値として、最大値Ｄmaxおよび最小値Ｄminが決定したら、次に、血流領域の２端点を特定する（ステップＳ４６０）。具体的には、以下の〔数式８〕に従った処理を実行することにより血流領域の２端点を特定する。
【００８１】
〔数式８〕
Ｈmax＝Ｄmax・α
ｉ＝ｓ／２からｉ＝０の順に、Ｄ（ｉp）≧Ｄ（ｉp＋１）かつＤ（ｉp）≧Ｄ（ｉp−１）かつＤ（ｉp）≧Ｈmaxを満たす極大値Ｄ（ｉp）を探索し、ｉ１＜ｉpでＤ（ｉ１）＜Ｈmaxを満たすサンプルｉ１を血流領域の第１の端点として特定
Ｈmin＝Ｄmin・α
ｉ＝ｓ／２からｉ＝ｓの順に、Ｄ（ｉn）≦Ｄ（ｉn＋１）かつＤ（ｉn）≦Ｄ（ｉn−１）かつＤ（ｉn）≦Ｈminを満たす極小値Ｄ（ｉn）を探索し、ｉ２＞ｉnでＤ（ｉ２）＞Ｈminを満たすサンプルｉ２を血流領域の第２の端点として特定
【００８２】
上記〔数式８〕において、ｉ＝ｓ／２からｉ＝０の順に極大値を探索し、ｉ＝ｓ／２からｉ＝ｓの順に極小値を探索することは、中点付近のサンプルから先頭または末尾のサンプル側に探索して血流領域の２端点を特定することを示している。また、上記〔数式８〕において、αは０以上１以下の実数値である。本実施形態では、α＝０．８に設定してある。係数αを用いることにより、第１所定値Ｈmax、第２所定値Ｈminが定まる。図１１の下段において、Ｄ（ｉ１）、Ｄ（ｉ２）で示す矢印の大きさが微分輝度値であり、これらの微分輝度値Ｄ（ｉ１）、Ｄ（ｉ２）をとるサンプルｉ１、ｉ２が３本ある血管のうちの中央の血管の血流領域の２端点となる。なお、上記ステップＳ４５０において、Ｄmax2をＤmax、Ｄmin2をＤminと置き換えているが、図１１においては、置き換えていない状態でＤmax2、Ｄmin2を図示している。
【００８３】
結局、ステップＳ４６０においては、血流領域端点特定手段５０は、最大サンプルより微分輝度値が小さく第１所定値Ｈmaxより微分輝度値が大きい極大点から見て先頭（ｉ＝０）のサンプル寄りで第１所定値Ｈmaxより小さい微分輝度値をとるサンプルを血流領域の第１の端点として特定し、最小サンプルより微分輝度値が大きく第１所定値Ｈmaxより小さい第２所定値Ｈminより微分輝度値が小さい極小点から見て末尾（ｉ＝ｓ）のサンプル寄りで第２所定値Ｈminより大きい微分輝度値をとるサンプルを血流領域の第２の端点として特定する処理を行っている。始点辺ＳＨについて、ステップＳ４２０〜ステップＳ４６０の処理を終えたら、計測範囲を終点辺ＥＨに変更して、同様に終点辺ＥＨに対してステップＳ４２０〜ステップＳ４６０の処理を実行し、終点辺ＥＨと交わる血流領域の２端点を特定する。
【００８４】
図１２は、両側の血管の一部が計測範囲に含まれる場合の輝度値、微分輝度値の関係を示す図である。図１２の上段のように、３本の血管のうち、中央の血管が計測範囲に完全に含まれ、両側の血管の一部が計測範囲に含まれる場合、ステップＳ４４０において算出された微分配列Ｄ（ｉ）の最大値Ｄmaxおよび最小値Ｄminは、いきなり対象とする中央の血流領域の微分輝度値の最大値および最小値となる。このため、ステップＳ４５０の処理を実行しても、〔数式７〕の条件を満たす極大値Ｄmax2、極小値Ｄmin2は存在せず、最大値Ｄmaxおよび最小値Ｄminは、置き換えられずに、ステップＳ４６０の処理が行われることになる。ステップＳ４６０の処理の結果、図１２に示したような、微分輝度値Ｄ（ｉ１）、Ｄ（ｉ２）をとるサンプルｉ１、ｉ２が中央の血管の血流領域の２端点として特定される。図１２の変形例として、図１２の上段の左端または右端の血管のいずれかが、図１１の上段の左端または右端の血管のいずれかのように、計測範囲に完全に含まれる場合が考えられる。その場合は、ステップＳ４５０の処理を実行すると、〔数式７〕の条件を満たす極大値Ｄmax2または極小値Ｄmin2のいずれか一方が存在し、最大値Ｄmaxまたは最小値Ｄminのいずれかが、置き換えられて、ステップＳ４６０の処理が行われることになる。
【００８５】
図１３は、血管が一本だけ計測範囲に含まれる場合の輝度値、微分輝度値の関係を示す図である。図１３の上段のように、血管が一本だけ計測範囲に含まれる場合、ステップＳ４４０において算出された微分配列Ｄ（ｉ）の最大値Ｄmaxおよび最小値Ｄminは、対象とする中央の血流領域の微分輝度値の最大値および最小値となる。このため、ステップＳ４５０の処理を実行しても、〔数式７〕の条件を満たす極大値Ｄmax2、極小値Ｄmin2は存在せず、最大値Ｄmaxおよび最小値Ｄminは、置き換えられずに、ステップＳ４６０の処理が行われることになる。ステップＳ４６０の処理の結果、図１３に示したような、微分輝度値Ｄ（ｉ１）、Ｄ（ｉ２）をとるサンプルｉ１、ｉ２が一本だけ含まれる血管の血流領域端点として特定される。図１３の第１の変形例として、図１２の上段の左端または右端の血管のように、もう１つの血管が計測範囲に部分的に含まれる場合が考えられる。その場合も、ステップＳ４５０の処理を実行しても、〔数式７〕の条件を満たす極大値Ｄmax2、極小値Ｄmin2は存在せず、最大値Ｄmaxおよび最小値Ｄminは、置き換えられずに、ステップＳ４６０の処理が行われることになる。更に、図１３の第２の変形例として、図１１の上段の左端または右端の血管のように、もう１つの血管が計測範囲に完全に含まれる場合が考えられる。その場合は、ステップＳ４５０の処理を実行すると、〔数式７〕の条件を満たす極大値Ｄmax2または極小値Ｄmin2のいずれか一方が存在し、最大値Ｄmaxまたは最小値Ｄminのいずれかが、置き換えられて、ステップＳ４６０の処理が行われることになる。
【００８６】
始点辺ＳＨ、終点辺ＥＨの双方に対して血流領域の２端点を特定してステップＳ４００の処理を終えたら、次に、動静脈径比算出手段６０が、血流幅比率である動静脈径比の算出を行う（ステップＳ５００）。具体的には、以下の〔数式９〕に従った処理を実行することにより動静脈径比を算出する。
【００８７】
〔数式９〕
動静脈径比＝（動脈に対応する血流領域の２端点間の距離の平均値）／（静脈に対応する血流領域の２端点間の距離の平均値）
【００８８】
上記〔数式９〕に示すように、動静脈径比は、（動脈に対応する血流領域の２端点間の距離の平均値）を（静脈に対応する血流領域の２端点間の距離の平均値）で除算することにより得られる。この平均値を算出することには、２つの意味があり、単一の計測用図形に基づいて、始点辺ＳＨより算出される血流領域の２端点間の距離と、終点辺ＥＨより算出される血流領域の２端点間の距離の２箇所の距離が同時に算出されるため、これらの平均値である平均距離を求めることと、動脈または静脈に対して計測用図形が複数配置されている場合、更に複数の計測用図形ごとに算出される平均距離どうしの平均値を求める必要がある。動脈に対して計測用図形がＮａ個、静脈に対して計測用図形がＮｖ個配置されて計測範囲が設定された場合（Ｎａ，Ｎｖ≧１でＮａ＝Ｎｖでなくても良い。）、動脈に対して配置された各計測用図形を基に算出される平均距離の総和をＮａで除算した値と、静脈に対して配置された各計測用図形を基に算出される平均距離の総和をＮｖで除算した値を除算したものが、動静脈径比となる。動静脈径比算出手段６０は、算出した動静脈径比を情報記憶手段９０に記憶させる。並行して表示手段２０が算出された動静脈径比を表示する。
【００８９】
動静脈径比の算出を終えたら、次に、静脈径比算出手段７０が、血流幅比率である静脈径比の算出を行う（ステップＳ６００）。具体的には、以下の〔数式９〕に従った処理を実行することにより動静脈交叉・静脈径比を算出する。
【００９０】
〔数式１０〕
動静脈交叉・静脈径比＝（各計測用図形に基づいて算出される２つの血流領域の２端点間の距離の小さい方の総和値）／（各計測用図形に基づいて算出される２つの血流領域の２端点間の距離の大きい方の総和値）
【００９１】
上記〔数式１０〕に示すように、動静脈交叉・静脈径比は、（各計測用図形に基づいて算出される２つの血流領域の２端点間の距離の小さい方の総和値）を（各計測用図形に基づいて算出される２つの血流領域の２端点間の距離の大きい方の総和値）で除算することにより得られる。ここで、総和値を算出するのは、複数の計測用図形が同一の動静脈交叉部または同一解析範囲内で異なる動静脈交叉部の静脈の複数個所に跨って配置された場合に対応するためで、その場合は、計測用図形ごとに算出される２つの距離の小さい方の値どうしの総和値および大きい方の値どうしの総和値を求める必要がある。計測用図形が１箇所のみ配置された計測範囲が設定された場合、当該計測用図形に基づいて算出される２つの距離の各々がそのまま総和値となり、２つの距離の小さい方の値を大きい方の値で除算することにより、動静脈交叉・静脈径比を算出する。静脈径比算出手段７０は、算出した動静脈交叉・静脈径比を情報記憶手段９０に記憶させる。並行して表示手段２０が算出された動静脈交叉・静脈径比を表示する。
【００９２】
静脈径比の算出を行ったら、次に、計測用図形の表示を行う（ステップＳ７００）。具体的には、始点辺に基づいて特定された血流領域の２端点、終点辺に基づいて特定された血流領域の２端点を結ぶ台形と、２対の２端点の中点どうしを結ぶ線分とからなる計測用図形を表示手段２０に表示する。具体的には、まず、ステップＳ４６０にて特定された血流領域の２端点ｉ１，ｉ２を用いて、始点辺については、Ｘ１１＝Ｘ（ｉ１），Ｙ１１＝Ｙ（ｉ１），Ｘ１２＝Ｘ（ｉ２），Ｙ１２＝Ｙ（ｉ２）とし、終点辺については、Ｘ２１＝Ｘ（ｉ１），Ｙ２１＝Ｙ（ｉ１），Ｘ２２＝Ｘ（ｉ２），Ｙ２２＝Ｙ（ｉ２）として、血流領域端点に対応する４頂点を特定する。
【００９３】
次に、表示手段２０が、特定された血流領域の４端点（Ｘ１１，Ｙ１１），（Ｘ１２，Ｙ１２），（Ｘ２１，Ｙ２１），（Ｘ２２，Ｙ２２）と、始点辺側の２端点の中点（Ｘ１，Ｙ１）＝（（Ｘ１１＋Ｘ１２）／２，（Ｙ１１＋Ｙ１２）／２）、終点辺側の２端点の中点（Ｘ２，Ｙ２）＝（（Ｘ２１＋Ｘ２２）／２，（Ｙ２１＋Ｙ２２）／２）に対応する画素を特定し、４端点にそれぞれ対応する４頂点を結ぶ台形と、中点を結ぶ線分に対応する画素を所定の画素値に置き換えて表示する処理を行う。この結果、利用者には、配置された計測用図形のどの範囲が血流領域として特定されたか、血流領域の２端点の値が適切に算出されたかを把握することが可能となる。図９（ｂ）に、特定された血流領域に合わせて再表示された計測用図形の一例を示す。網掛けした部分は重ねて表示される眼底画像の血流領域に対応する部分である。ステップＳ４００により血流領域の２端点が特定されると、上述のように、ステップＳ７００において、２つの辺に対応して特定された血流領域の各々２端点を結ぶ２辺を含む台形と始点辺側の２端点の中点と終点辺側の２端点の中点を結ぶ線分からなる図形が表示されることになる。この図形は、実際には、演算により特定された血流領域の４端点を結ぶ台形と、始点辺側の２端点の中点と終点辺側の２端点の中点を結ぶ線分を表示したものである。しかし、演算処理は人間の感覚からすれば一瞬で終わるので、計測用図形を配置した後、計測処理の実行を指示すると、瞬時に長方形の計測用図形が、血流領域に合わせた幅の台形の計測ゲージに変化したように見える。
【００９４】
ステップＳ５００、Ｓ６００、Ｓ７００の処理は、必ずしも行う必要はなく、また実行する場合にいずれを先に行っても良い。ステップＳ４００により得られた血流領域の端点の情報を、血流領域端点特定手段５０が情報記憶手段９０に直接記憶させ、処理を終えるようにすることも可能である。
【００９５】
図１４は、実際の眼底画像の解析範囲の静脈上に、処理後の計測用図形を重ねて表示させた例を示す図である。図７（ａ）に示した状態で配置された計測用図形は、図１４に示したように、血流領域の幅に合わせて表示される。始点辺ＳＨと終点辺ＥＨにおける幅が異なる場合は、計測用図形は台形として表示される。したがって、図７（ａ）に示した長方形の計測用図形が血流領域にフィッティング（位置合わせ）されて、台形に変形したように見える。
【００９６】
図１５は、解析範囲内の６箇所（動脈３箇所、静脈３箇所）に計測用図形が配置された例を示す図である。このうち、図１５（ａ）は、解析範囲の画像を示し、図１５（ｂ）は、複数の計測用図形に基づく動静脈径比（血流幅比率）の算出結果を示している。一般に、動脈の方が静脈よりも細いため（健常者の中には動脈の方が太い場合もある）、動静脈径比は、０．５以上が正常である。図１５の例では、動静脈径比が０．８６であるため正常となっている。
【００９７】
図１６は、解析範囲内の６箇所（動脈３箇所、静脈３箇所）に計測用図形が配置された例を示す図である。このうち、図１６（ａ）は、眼底画像全体を示し、図１６（ｂ）は、解析範囲の画像を示し、図１６（ｃ）は、複数の計測用図形に基づく動静脈径比（血流幅比率）の算出結果を示している。図１６の例では、動静脈径比が０．３７であり、０．５未満であるため動脈硬化進行となっている。
【００９８】
動静脈交叉・静脈径比を算出する場合の例について説明しておく。図１７は、図４と同一の眼底画像であるが、解析範囲が動静脈交叉部に設定されている。ステップＳ１００で図１７に示したような解析範囲が設定された場合における、解析範囲の表示例を図１８（ａ）に示す。図１８（ａ）の例では、解析範囲である横５１５画素×縦５１５画素の眼底画像の動静脈交叉部に、長方形状の計測用図形が配置された状態を示している。図１８（ｂ）は、図１８（ａ）に示した計測用図形に対応する計測範囲設定画面を示す図である。図１８（ｂ）に示した計測範囲設定画面では、図７（ｂ）に示した計測範囲設定画面と同様の設定を行うことができる。
【００９９】
図１９は、図４と同様、実際の眼底画像の解析範囲の静脈上に、処理後の計測用図形を重ねて表示させた例を示す図である。図１８（ａ）に示した状態で配置された計測用図形は、図１９に示したように、血流領域の幅に合わせて表示される。始点辺ＳＨと終点辺ＥＨにおける幅が異なる場合は、計測用図形は台形として表示される。
【０１００】
図２０は、解析範囲内の２箇所（単一の動静脈交叉から対向して出ている同一静脈２箇所）に計測用図形が配置された例を示す図である。このうち、図２０（ａ）は、解析範囲の画像を示し、図２０（ｂ）は、２箇所の計測用図形に基づく動静脈交叉・静脈径比（血流幅比率）の算出結果を示している。動静脈交叉・静脈径比も、０．５以上が正常である。図２０の例では、動静脈交叉・静脈径比が０．６０であるため正常となっている。
【０１０１】
図２１は、解析範囲内の２箇所（単一の動静脈交叉から対向して出ている同一静脈２箇所）に計測用図形が配置された例を示す図である。このうち、図２１（ａ）は、眼底画像全体を示し、図２１（ｂ）は、解析範囲の画像を示し、図２１（ｃ）は、２箇所の計測用図形に基づく動静脈交叉・静脈径比（血流幅比率）の算出結果を示している。図２１の例では、動静脈交叉・静脈径比が０．４９であり、０．５未満であるため動脈硬化進行となっている。
【０１０２】
＜３．計測用図形表示装置としての機能＞
上記実施形態では、眼底画像処理装置について説明したが、眼底画像処理装置は、計測用図形を表示する表示装置である計測用図形表示装置としての機能も果たす。すなわち、眼底画像に含まれる血管の血流領域端点を特定する際に、指示入力手段１０からの指示に基づき、外形が台形（一態様として長方形を含む意味）であり、平行な２辺の中点を結ぶ基準線を備える計測用図形を表示手段２０に表示する。そして、指示入力手段１０からの指示に基づき、計測範囲設定手段４０が、配置された計測用図形の対向する２辺それぞれを計測範囲として設定する。さらに、血流領域端点特定手段５０が、設定された計測範囲内の平行な２辺の各々に基づいて血流領域の２端点を特定し、２辺に対して特定された血流領域の４端点を４頂点とする外形が台形形状の図形を計測用図形に代えて表示手段２０に表示する。すなわち、図１、図２に示した眼底画像処理装置を計測用図形表示装置として機能させることができる。この場合、計測用図形がいわゆる計測ゲージのように、変形して血流領域にフィッティングしたかのような印象を与えるとともに、血流領域における血流幅を直感的に把握し易くなる。
【０１０３】
＜４．変形例等＞
以上、本発明の好適な実施形態について説明したが、本発明は上記実施形態に限定されず、種々の変形が可能である。例えば、上記実施形態では、撮影により得られた眼底画像としてフルカラー画像を用いたが、赤外線光源のデジタル眼底カメラで撮影されるモノクロ画像も利用されつつあり、その場合は、眼底血管が鮮明に造影されるため、血流領域を特定する処理における二値化処理を単純な判別分析法で簡便に高精度に行うことができる。また、近年業務用のデジタルカメラの階調は１０ビット以上に拡大しており、各色１０２４階調等のフルカラー画像を取得することも可能になっており、より階調数の多いカラー画像を用いてもよい。
【符号の説明】
【０１０４】
１・・・ＣＰＵ（Central Processing Unit）
２・・・ＲＡＭ（Random Access Memory）
３・・・記憶装置
４・・・タッチパネル入力Ｉ／Ｆ
５・・・データ入出力Ｉ／Ｆ
６・・・表示部
７・・・画像入力部
１０・・・指示入力手段１０
２０・・・表示手段
３０・・・画像拡大手段
４０・・・計測範囲設定手段
５０・・・血流領域端点特定手段
６０・・・動静脈径比算出手段
７０・・・静脈径比算出手段
８０・・・眼底画像記憶手段
９０・・・情報記憶手段
１００・・・眼底画像処理装置
【図１】
2018011726.tif 000003
【図２】
2018011726.tif 000004
【図３】
2018011726.tif 000005
【図４】
2018011726.tif 000006
【図５】
2018011726.tif 000007
【図６】
2018011726.tif 000008
【図７】
2018011726.tif 000009
【図８】
2018011726.tif 000010
【図９】
2018011726.tif 000011
【図１０】
2018011726.tif 000012
【図１１】
2018011726.tif 000013
【図１２】
2018011726.tif 000014
【図１３】
2018011726.tif 000015
【図１４】
2018011726.tif 000016
【図１５】
2018011726.tif 000017
【図１６】
2018011726.tif 000018
【図１７】
2018011726.tif 000019
【図１８】
2018011726.tif 000020
【図１９】
2018011726.tif 000021
【図２０】
2018011726.tif 000022
【図２１】
2018011726.tif 000023