(57)【要約】      （修正有）【課題】瞳孔領域を容易に特定でき、特定した瞳孔領域を用いて白内障レベルを判別できる白内障検査装置および白内障判定プログラムを提供する。【解決手段】白内障検査装置は、リング状投射光２０として眼球の瞳孔領域に略正面から投射する投光手段と、リング状投射光２０の眼球のリング状反射光１０，１１，１２を撮像し画像情報とする撮像手段と、リング状反射光を検出する反射光検出手段と、リング状反射光を用いて白内障を判定する白内障判定手段を備える。また、リング状反射光のうち最大径かつ最大輝度の反射光１０の内輪郭の内側の瞳孔領域の色特徴を解析する色特徴解析手段を備え、瞳孔領域の色特徴から核白内障あるいは後嚢下白内障を判定する。瞳孔領域からまつげを判別し、まつげを除いて模様が存在している領域面積および色から混濁度合いを解析し、瞳孔領域の混濁度合いから皮質白内障を判定する。【選択図】図２
2015216939.tif 000002
【特許請求の範囲】
【請求項１】
リング状投射光として眼球の瞳孔領域に略正面から投射する投光手段と、
前記リング状投射光の眼球のリング状反射光を撮像し画像情報とする撮像手段と、
前記リング状反射光を検出する反射光検出手段と、
前記リング状反射光を用いて白内障を判定する白内障判定手段と、を備えることを特徴とする白内障検査装置。
【請求項２】
前記リング状反射光のうち最大径かつ最大輝度の反射光の内輪郭の内側の瞳孔領域の色特徴を解析する色特徴解析手段を更に備え、
前記白内障判定手段において、前記瞳孔領域の色特徴から核白内障あるいは後嚢下白内障を判定することを特徴とする請求項１に記載の白内障検査装置。
【請求項３】
前記瞳孔領域からまつげを判別して、まつげを除いて模様が存在している領域面積および色から混濁度合いを解析する混濁解析手段を更に備え、
前記白内障判定手段において、前記瞳孔領域の混濁度合いから皮質白内障を判定することを特徴とする請求項２に記載の白内障検査装置。
【請求項４】
前記瞳孔領域におけるＲＧＢ画像をＨＳＶ画像に変換する手段と、
前記ＲＧＢ画像を二値化画像に変換する手段と、
前記ＨＳＶ画像におけるＨ成分，Ｓ成分，Ｖ成分と各成分の最大値，平均値，中央値を算出する色特徴量算出手段と、
或は、
前記二値化画像における中心部からの距離に応じてＮ等分（Ｎは２以上）して分けた領域毎に混濁の疑義がある部分の画素数を算出して特徴ベクトルを算定する混濁度合い算定手段と、
を更に備え、
前記白内障判定手段が、前記特徴ベクトルを、機械学習による識別器を用いて、核白内障、皮質白内障、その他の白内障の有無と程度を判定することを特徴とする請求項２又は３に記載の白内障検査装置。
【請求項５】
前記白内障判定手段において、眼内の各透光体の境界面で反射される、前記リング状反射光の個数により、白内障と判定することを特徴とする請求項１〜４の何れかに記載の白内障検査装置。  
【請求項６】
前記撮像手段および前記投光手段を一端に備えるハウジングであって、該ハウジングは、他端が開放され被験者の眼球の周囲に覆い被せた場合に、前記光源からの投射光以外の入射を遮るよう構成され、
被験者が前記光源方向を見るようにする位置調整指標が設けられ、
前記ハウジングの端部は、被験者のまぶたを開くようにまぶたの上下方向に変形する機構を備えたことを特徴とする請求項１〜５の何れかに記載の白内障検査装置。
【請求項７】
コンピュータを用いて眼球の白内障を判定する白内障判定プログラムであって、
前記コンピュータを、
リング状投射光として眼球の瞳孔領域に略正面から投射したときの前記眼球の角膜および水晶体でのリング状反射光の画像情報を取得する画像情報取得手段、
前記リング状反射光を検出する反射光検出手段、
前記リング状反射光を用いて白内障を判定する白内障判定手段、
として機能させる白内障判定プログラム。
【請求項８】
前記コンピュータを、更に、
前記リング状反射光のうち最大径かつ最大輝度の反射光の内輪郭の内側の瞳孔領域の色特徴を解析する色特徴解析手段として機能させ、
前記白内障判定手段において、前記瞳孔領域の色特徴から核白内障あるいは後嚢下白内障を判定することを特徴とする請求項７に記載の白内障判定プログラム。
【請求項９】
前記コンピュータを、更に、
前記瞳孔領域からまつげを判別するまつげ判別手段、
まつげを除いて模様が存在している領域面積および色から混濁度合いを解析する混濁度合い解析手段、として機能させ、
前記白内障判定手段において、前記瞳孔領域の混濁度合いから皮質白内障を判定することを特徴とする請求項８に記載の白内障判定プログラム。
【請求項１０】
前記白内障判定手段において、
前記コンピュータを、
前記瞳孔領域におけるＲＧＢ画像をＨＳＶ画像に変換する手段、
前記ＲＧＢ画像を二値化画像に変換する手段、
前記ＨＳＶ画像におけるＨ成分，Ｓ成分，Ｖ成分と各成分の最大値，平均値，中央値を算出する色特徴量算出手段と、
或は、
前記二値化画像における中心部からの距離に応じてＮ等分（Ｎは２以上）して分けた領域毎に混濁の疑義がある部分の画素数を算出して特徴ベクトルを算定する混濁度合い算定手段、
前記特徴ベクトルを、機械学習による識別器を用いて、核白内障、皮質白内障、その他の白内障の有無と程度を判定する手段、として機能させることを特徴とする請求項８又は９に記載の白内障判定プログラム。
【請求項１１】
前記白内障判定手段において、眼内の各透光体の境界面で反射される、前記リング状反射光の個数により、白内障と判定することを特徴とする請求項７〜１０の何れかに記載の白内障検査プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、白内障診断を簡易に行うための白内障検査装置および白内障判定プログラムに関するものである。
【背景技術】
【０００２】
従来、白内障診断を行うための機器としては、下記の特許文献１や特許文献２に開示されているものが知られているが、いずれも白内障の進行度合いを正確に診断するために分光光学系や光感知検知器などの特別な機器を必要としていた。そのため、特に発展途上国で使用する場合には、機器導入コストや機器を扱う熟練した医師の不足などが問題となっている。
また、分光光学系や光感知検知器などの特別な機器を必要とせず、簡便な機材を用いる構成で、誰でもどこでも簡便に白内障の検査が行える白内障検査装置が知られている（特許文献３）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００２−２２４０４１号公報
【特許文献２】特表２００４−５３５８８０号公報
【特許文献３】特許第５３０５４０９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
上述した特許文献３に開示された白内障検査装置では、点形状の光源を使用して被験者の眼の画像を撮影し、画像解析によって白内障発病の有無を判定している。そのため、瞳孔領域の特定が困難であり、瞳孔領域画像を解析して白内障の重症度の判定を行うのは困難であった。
かかる状況に鑑みて、本発明は、簡便な機材を用いる構成で、瞳孔領域を容易に特定でき、特定した瞳孔領域を用いて白内障レベルを判別できる白内障検査装置および白内障判定プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００５】
本発明者らは、白内障の簡易検査（スクリーニング）を行うための装置および方法を鋭意検討した結果、本発明を完成した。
本発明の白内障検査装置は、下記１）〜４）の手段を備える。
１）リング状投射光として眼球の瞳孔領域に略正面から投射する投光手段
２）リング状投射光の眼球のリング状反射光を撮像し画像情報とする撮像手段
３）リング状反射光を検出する反射光検出手段
４）リング状反射光を用いて白内障を判定する白内障判定手段
【０００６】
本発明の白内障検査装置では、光源からの光をリング状投射光とするのが特徴である。光源からの光をリング状投射光とすることにより、眼球を撮影した画像の中で瞳孔領域を容易に特定することができ、瞳孔領域の解析による白内障の簡易検査（スクリーニング）を行える。リング状反射光を用いて白内障の症状を判定することにより、リング状反射光を用いない場合に比べて、白内障の判定精度を向上することができる。
【０００７】
ここで、リング状投射光により、なぜ瞳孔領域を容易に特定することができるのかという点について説明する。一般的に、眼球の撮影画像において、虹彩の内輪郭を判別すれば瞳孔領域を特定することが可能である。しかし、白内障眼の場合、虹彩の内輪郭に濁りがあると、虹彩の内輪郭が不明確になる。特に、皮質白内障の患者の場合、水晶体の周りから瞳の中に向かってまだらに濁る症状であり、虹彩の内輪郭が不明確になる。このため、リング状投射光を眼球の正面より投射して、瞳孔中心に同心円状の反射光を得る。角膜表面はほぼ半球形をしていることから、正面から投射されるリング状投射光は、角膜表面から強く反射されてリング状反射光となり、円状の画像が得られる。
【０００８】
リング状反射光は、画像処理上では真白に色が飽和した領域、すなわち、ＲＧＢ値が限りなく２５５に近い領域となるため、リング状反射光の内輪郭を抽出するのは容易となる。
また、リング状反射光は、瞳孔中心から同心円状に観測されるのが本来の形であり、瞳孔中心から反射光がずれている場合、それは撮影時の位置合わせの問題に起因する。正常眼において、リング状反射光は４つの境界面（角膜前面、角膜後面、水晶体前面、水晶体後面）のそれぞれで観測されることから、撮影画像には４つのリング状反射光が存在することになるが、水晶体前面は非常に反射率が低いことから、水晶体前面によるリング状反射光は観測が困難である。従って、正常眼におけるリング状反射光は、主に角膜前面による反射光が最も強く観測され、後は、角膜後面と水晶体後面による反射光の３つのリング状反射光が観測できることになる。一方、白内障眼、特に中心性の核白内障の場合は、角膜前面と角膜後面による反射光の２つのリング状反射光しか観測されない。このように撮影した画像にあるリング状の反射光個数を測定することにより、簡易的に正常眼と白内障眼を分けることができる。
上記の判別精度を向上するには、瞳孔を真正面から中心を合わせて撮影する。瞳孔を真正面から中心を合わせて撮影することにより、リング状反射光が同心円状に並ぶため、リング状反射光の識別および反射光の個数の識別が容易となる。
【０００９】
投光手段は、撮像手段と連動するカメラフラッシュであることが好ましい。撮像手段として汎用的なデジタルカメラを使用し、カメラに内蔵された若しくは外付けのフラッシュで、白内障の検査（スクリーニング）を行うことができる。
【００１０】
本発明の白内障検査装置は、リング状反射光のうち最大径かつ最大輝度の反射光の内輪郭の内側の瞳孔領域の色特徴を解析する色特徴解析手段を更に備え、白内障判定手段において、瞳孔領域の色特徴から核白内障あるいは後嚢下白内障を判定する。
核白内障とは、核と呼ばれる水晶体の真ん中（中心）が濁る症状であり、瞳孔領域の周辺部と中心部の色特徴の違いが明瞭であることから判定することができる。一方、白内障のうち、最も多いのが、皮質白内障であり、水晶体の周りから瞳の中に向かってまだらに濁る症状である。また、後嚢下白内障は、水晶体の裏側に近いところで真ん中からスリガラスのように濁る症状である。
【００１１】
また、本発明の白内障検査装置は、瞳孔領域からまつげを判別して、まつげを除いて模様が存在している領域面積および色から混濁度合いを解析する混濁解析手段を更に備え、白内障判定手段において、瞳孔領域の混濁度合いから皮質白内障を判定する。皮質白内障は、水晶体の周りから瞳の中に向かってまだらに濁る症状であり、混濁の領域が広範囲であり、混濁エリアも水晶体の周りに多く分布する傾向があり、それらの特徴から皮質白内障を判別する。
【００１２】
また、まつげは特徴的な形をしていることから、その形状の特徴を用いて、まつげの画像の取り除き、まつげの画像を皮質白内障の混濁と峻別することにより、虹彩の内輪郭付近の瞳孔領域の混濁から、まつげによる混濁を排除し、皮質白内障の判定精度を向上する。
ここで、模様とは色ムラのことであり、混濁度合いにより色ムラが生じることから、それらを模様として捉えている。
【００１３】
また、本発明の白内障検査装置は、下記ａ）〜ｄ）の手段を更に備え、白内障判定手段が、下記の特徴ベクトルを、機械学習による識別器を用いて、核白内障、皮質白内障、その他の白内障の有無と程度を判定することが好ましい。ここで、ｃ）とｄ）の手段は、少なくとも一方の手段を有することで白内障を判定できる。
ａ）瞳孔領域におけるＲＧＢ画像をＨＳＶ画像に変換する手段
ｂ）ＲＧＢ画像を二値化画像に変換する手段
ｃ）ＨＳＶ画像におけるＨ成分，Ｓ成分，Ｖ成分と各成分の最大値，平均値，中央値を算出する色特徴量算出手段
ｄ）二値化画像における中心部からの距離に応じてＮ等分（Ｎは２以上）して分けた領域毎に混濁の疑義がある部分の画素数を算出して特徴ベクトルを算定する混濁度合い算定手段
【００１４】
また、本発明の白内障検査装置の白内障判定手段において、眼内の各透光体の境界面で反射される、リング状反射光の個数により、白内障と判定することができる。
上述の如く、リング状反射光は、瞳孔中心から同心円状に観測され、正常眼の場合、リング状反射光は、主に角膜前面による反射光と角膜後面と水晶体後面による反射光の３つのリング状反射光が観測されるのに対して、白内障眼、特に中心性の核白内障の場合は、角膜前面と角膜後面による反射光の２つのリング状反射光しか観測されない。従って、リング状反射光の個数を測定することにより、正常眼と白内障眼を判別できる。なお、眼内の各透光体の境界面で反射される反射光より、角膜前面で反射されるリング状反射光を用いて瞳孔領域を特定し、その中の反射光の個数により白内障の有無を判定すると共に、上述のａ）〜ｃ）の手段を用いて、核白内障、皮質白内障、その他の白内障の有無と程度を判定することもできる。
なお、下記表１に、角膜前面（角膜表皮）、角膜後面（角膜内皮）、水晶体前面および水晶体後面のそれぞれの反射面の反射像輝度比を示す。水晶体前面の反射像の輝度が、他の反射面よりも一桁以上、輝度が小さく、反射が少ないことがわかる。
【００１５】
【表１】
2015216939.tif 000003
【００１６】
また、本発明の白内障検査装置は、撮像手段および投光手段を一端に備えるハウジングであって、該ハウジングは、他端が開放され被験者の眼球の周囲に覆い被せた場合に、光源からの投射光以外の入射を遮るよう構成され、被験者が光源方向を見るようにする位置調整指標が設けられ、ハウジングの端部は、被験者のまぶたを開くようにまぶたの上下方向に変形する機構を備えたことが好ましい。
リング状反射光内に虹彩が映り込んでしまうと、それが濁りと判定されるので、これを回避すべく、リング状反射光は必ず虹彩の領域の内側になるよう位置調整することが必要である。そのため、散瞳状態に近い状態、すなわち、瞳孔が大きく開いている状態にすべく、ハウジングにより暗所をつくり、暗順応により瞳孔が開いている状態で撮影する。
【００１７】
ハウジングを用いることで、検査の際に、周囲を暗くしなくてもハウジングの筐体で外部光の入射を遮断し、暗室状態にできる。ここで、ハウジングの形状は、例えば、双眼鏡の形状や、筒状、箱状のものでもよい。なお、ハウジングを暗箱状態にして、所定秒数以上経過後に、撮像手段と投光手段を動作することがより好ましい。瞳孔を大きく開かせるためである。所定秒数とは、例えば３〜５秒、或は、５秒以上の時間経過などである。
また、リング状投射光は眼球の正面より投射して、瞳孔中心に同心円状のリング状反射光を得ることが好ましく、そのため、被験者が光源方向を見るようにする位置調整指標を設ける。ここで、位置調整指標とは、例えば、眼球の正面でなければ見えないＬＥＤ点滅マーカーなどである。
【００１８】
さらに、ハウジングの端部は、被験者のまぶたを開くようにまぶたの上下方向に変形する機構を備える。まぶたの開きが小さいと、まぶたによって瞳孔領域が隠され、また、まつげの影響を受けやすくなる。そのため、被験者のまぶたを強制的に開くように、まぶたの上下の皮膚と接して持ち上げ／持ち下げる機構によって、まぶたが開くようにする。
【００１９】
次に、本発明の白内障判定プログラムについて説明する。
本発明の白内障判定プログラムは、コンピュータを下記（１）〜（３）の手段として機能させ、眼球の白内障を判定するプログラムである。
（１）リング状投射光として眼球の瞳孔領域に略正面から投射したときの眼球の角膜および水晶体でのリング状反射光の画像情報を取得する画像情報取得手段
（２）リング状反射光を検出する反射光検出手段
（３）リング状反射光を用いて白内障を判定する白内障判定手段
【００２０】
本発明の白内障判定プログラムにおいて、コンピュータを、更に、リング状反射光のうち最大径かつ最大輝度の反射光の内輪郭の内側の瞳孔領域の色特徴を解析する色特徴解析手段として機能させ、白内障判定手段において、瞳孔領域の色特徴から核白内障あるいは後嚢下白内障を判定することが好ましい。
【００２１】
本発明の白内障判定プログラムにおいて、コンピュータを、更に、瞳孔領域からまつげを判別するまつげ判別手段と、まつげを除いて模様が存在している領域面積および色から混濁度合いを解析する混濁度合い解析手段として機能させ、白内障判定手段において、瞳孔領域の混濁度合いから皮質白内障を判定することが好ましい。
【００２２】
本発明の白内障判定プログラムにおける白内障判定手段において、コンピュータを、下記（ａ）〜（ｅ）の手段として機能させることが好ましい。ここで、ｃ）とｄ）の手段は、少なくとも一方の手段として機能させることで白内障を判定できる。
（ａ）瞳孔領域におけるＲＧＢ画像をＨＳＶ画像に変換する手段
（ｂ）ＲＧＢ画像を二値化画像に変換する手段
（ｃ）ＨＳＶ画像におけるＨ成分，Ｓ成分，Ｖ成分と各成分の最大値，平均値，中央値を算出する色特徴量算出手段
（ｄ）二値化画像における中心部からの距離に応じてＮ等分（Ｎは２以上）して分けた領域毎に混濁の疑義がある部分の画素数を算出して特徴ベクトルを算定する混濁度合い算定手段
（ｅ）特徴ベクトルを、機械学習による識別器を用いて、核白内障、皮質白内障、その他の白内障の有無と程度を判定する手段
【００２３】
本発明の白内障判定プログラムにおける白内障判定手段において、コンピュータを、眼内の各透光体の境界面で反射される、リング状反射光の個数により、白内障と判定する手段として機能させることが好ましい。コンピュータの画像処理で画像中のリング状反射光の個数を測定して、正常眼と白内障眼を判別する。
すなわち、本発明の白内障判定プログラムにおける白内障判定手段においては、眼内の各透光体の境界面で反射される反射光より、角膜前面で反射されるリング状反射光を用いて瞳孔領域を特定し、その中の反射光の個数により白内障の有無を判定すると共に、前述の（ａ）〜（ｄ）の手段を用いて、核白内障、皮質白内障、その他の白内障の有無と程度を判定する。
【発明の効果】
【００２４】
本発明の白内障検査装置および白内障判定プログラムによれば、リング状の光を用いることにより瞳孔領域を容易に特定でき、瞳孔領域の色特徴および混濁度合いを用いて白内障のレベルを判別できるといった効果がある。
発展途上国など眼科医の不足するところで、診断装置設備が不足する場所においても、白内障発病患者に対して、デジタルカメラとコンピュータの組み合わせといった簡便な機器を用いて、白内障の種類まで自動的に推定することができるようになる。
【図面の簡単な説明】
【００２５】
【図１】実施例１の白内障検査装置の概念図
【図２】実施例１の白内障検査装置の構成断面図
【図３】角膜および水晶体前面のリング状反射光の観測画像の説明図
【図４】実施例１の白内障検査装置の処理フロー図
【図５】色特徴の解析処理フロー図
【図６】瞳孔領域画像と極座標展開画像の説明図
【図７】瞳孔領域画像のＲＧＢ画像とＨＳＶ色空間の説明図
【図８】極座標展開画像の３つの特徴量の説明図
【図９】極座標展開画像の説明図
【図１０】極座標展開画像の一例
【図１１】まつげと混濁の識別についての説明図
【図１２】混濁の解析の説明図
【図１３】特徴ベクトルの説明図
【図１４】識別器を用いた判定の概念図
【図１５】１０分割交差検定による評価方法の説明図
【図１６】決定木による識別フローを可視化した図（１）
【図１７】決定木による識別フローを可視化した図（１）
【発明を実施するための最良の形態】
【００２６】
以下、本発明の実施形態の一例を、図面を参照しながら詳細に説明していく。なお、本発明の範囲は、以下の実施例や図示例に限定されるものではなく、幾多の変更及び変形が可能である。
【実施例１】
【００２７】
図１は実施例１の白内障検査装置の概念図を示している。また、図２は白内障検査装置の構成断面図を示している。図１，２に示すように、実施例１の白内障検査装置は、光源からの光をリング状投射光２０として眼球１に投射する投光手段であるリング状フラッシュ２２と、リング状フラッシュ２２の光の眼球１の水晶体２及び角膜６でのリング状反射光を撮像して画像情報とする撮像手段であるカメラ２１と、コンピュータ２３から構成される。コンピュータ２３は、画像情報を取得する画像情報取得手段２３ａと、リング状反射光１０を検出するリング状反射光検出手段２３ｂと、リング状反射光のうち最大径かつ最大輝度の反射光の内輪郭の内側の瞳孔領域の色特徴を解析する色特徴解析手段２３ｃと、瞳孔領域の混濁度合いを解析する混濁解析手段２３ｄと、リング状反射光１０の個数を解析するリング状反射光個数解析手段２３ｅと、リング状反射光１０を用いて、白内障を判定する白内障判定手段２３ｆから構成される。カメラ２１とリング状フラッシュ２２と被験者の眼球１はハウジング７で覆われている。ハウジング７は、暗箱や黒幕など外光を遮断するものである。
【００２８】
図１，２に示すように、リング状フラッシュ２２は、リング状の開孔を有するリング状スリット２２ａと遮蔽板２２ｂと発光部２２ｄ及び発光部をマウントするマウント部２２ｃで構成される。リング状フラッシュ２２の発光部２２ｄを覆うように遮光板２２ｂ及びマウント部２２ｃがあり、その遮光板２２ｂとフラッシュ本体の隙間がリング状スリット２２ａになっている。そして、カメラレンズ２１ａの光学的中心とリング光フラッシュ２２の幾何学的中心が同一線上に配置されるようにする。すなわち、カメラレンズの光学的中心とリング光投射光の光学的中心が同一線上になるようにする。これにより、眼球の略正面からリング状投射光２０を投射でき、角膜および水晶体の正面反射を得られるようにする。なお、発光部２２ｄの電力は、バッテリーユニット２５から電源ケーブル２６を介して供給している。
【００２９】
カメラ２１は、汎用のデジタルカメラを使用できる。また、リング状フラッシュ２２はカメラ２１の外付けのフラッシュを用いてもよく、光ファイバーケーブルによりリング状フラッシュ２２に導光してもよい。そして、カメラ２１のシャッターと連動してリング状フラッシュ２２の光源からリング状投射光を投射する。
コンピュータ２３は、カメラ２１から画像情報データを取り込めるよう、ＵＳＢ（Universal Serial Bus）ケーブル（データ通信インタフェース２４）を介してカメラ２１と接続されている。
【００３０】
カメラ２１で、リング状反射光を観測した場合、例えば、図３（１）あるいは図３（２）に示すような観測画像が得られる。
最大径かつ最大輝度の反射光であるリング状反射光１０は、眼球１の角膜６の前面６ａの反射光である。正常眼では図２に示すように、リング状反射光は４つの境界面（角膜前面６ａ、角膜後面６ｂ、水晶体前面２ａ、水晶体後面２ｂ）のそれぞれで観測され、撮影画像には４つのリング状反射光が存在する。しかし、水晶体前面２ａは非常に反射率が低いことから、水晶体前面２ａによるリング状反射光は観測が一般的に困難であり、図２および図３（２）に示すように、３つの境界面（角膜前面６ａ、角膜後面６ｂ、水晶体後面２ｂ）の３つのリング状反射光（角膜前面６ａのリング状反射光１０、角膜後面６ｂのリング状反射光１１、水晶体後面２ｂのリング状反射光１２）が観測できる。一方、白内障眼、特に中心性の核白内障の場合は、２つの境界面（角膜前面６ａ、角膜後面６ｂ）の２つのリング状反射光（角膜前面６ａのリング状反射光１０、角膜後面６ｂのリング状反射光１１）しか観測されない。
本実施例では、図３（１）に示すように、最大径かつ最大輝度の角膜前面６ａの反射光であるリング状反射光１０だけを捉え、その内輪郭の内側の瞳孔領域の色特徴や混濁度合いで白内障を判別する。
【００３１】
図４は、実施例１の白内障検査装置の処理フロー図を示している。
先ず、リング状フラッシュを点灯させて、リング状投射光を眼球に投射し、眼球を撮影する（ステップＳ０１）。次に、角膜前面でのリング状反射光の撮影画像から瞳孔領域を検出する（ステップＳ０２）。そして、角膜前面でのリング状反射光を用いて、その内輪郭の内側の領域を瞳孔領域画像として切り出す（ステップＳ０３）。
具体的には、瞳孔領域の検出は、撮影画像を白黒濃淡画像に変換して２値化し、細線化してＨｏｕｇｈ変換によりリング状反射光の円形を検出する。
【００３２】
そして、リング状反射光の個数を解析し（ステップＳ０４）、或は、色特徴を解析し（ステップＳ０５）、或は、混濁度合いを解析する（ステップＳ０６）。リング反射光の個数だけを解析して（ステップＳ０４）、リング反射光の個数を検出するだけでも、簡易的に白内障発病の有無を判定することができる。しかし、画像処理による自動判定を精度よく行うために、色特徴の解析（ステップＳ０５）や混濁度合いの解析（ステップＳ０６）もあわせて利用する。
色特徴の解析（ステップＳ０５）や混濁度合いの解析（ステップＳ０６）の結果と、予め学習データから抽出していた判定基準とを照らし合せて、白内障発病（白内障眼）の判定を行う（ステップＳ０７）。判定の結果、白内障発病（白内障眼）または健常眼を判定する（ステップＳ０８，０９）。
【００３３】
なお、本実施例では、白内障発症の有無のみを識別して、発症している場合は眼科医への受診を促すようなシステムを想定している。すなわち、被験者の眼の撮影画像から、白内障の発症の有無を推定し、医療機関に対して、眼球の画像情報と推定結果情報を統合された情報として、被験者の識別子（ＩＤ）と共に医療機関（医師の下）に送信し、被検者の受診を促して、白内障の治療が行われるようにするシステムを想定している。
【００３４】
次に、色特徴の解析処理フローについて図５を参照して説明する。
切り出した瞳孔領域画像から、ＲＧＢ画像を取得し（ステップＳ２１）、それをＨＳＶ画像に変換してＨＳＶ画像を取得する（ステップＳ２２）。そして、ＨＳＶ画像のそれぞれの平均値、最大値、中央値を測定する（ステップＳ２３）。測定した平均値、最大値、中央値を用いて、分類して色特徴を解析する（ステップＳ２４）。
【００３５】
次に、図６〜９を参照して、切り出した瞳孔領域画像の解析について説明する。
図６（１）は、眼球の撮影画像においてリング状反射光の形状・位置に基づいて、瞳孔領域を検出している様子を示している。
図６（２）は、切り出した瞳孔領域画像を示している。この瞳孔領域画像から色特徴を解析する。また、図６（３）は、切り出した瞳孔領域画像（円形状）を極座標展開した極座標展開画像を示している。極座標展開画像からリング状反射光のリング状フラッシュ模様を解析する。
なお、瞳孔領域画像において、健常者であればリング状反射光は円形状を示す。一方、白内障の患者であれば、リング状反射光は輝度が著しく低下するか、一部が欠落した形状を示す。
【００３６】
先ず、色特徴の検出方法について説明する。
図７（１）は、切り出した瞳孔領域画像（ＲＧＢ画像）を示す。上述の如く、ＲＧＢ画像から、それをＨＳＶ画像に変換してＨＳＶ画像を取得して、ＨＳＶ色空間で色特徴を解析する。
ＨＳＶ色空間は、色相（Ｈ），彩度（Ｓ）および明度（Ｖ）の３つの成分からなる色空間であり、円錐を用いて視覚化できる。円錐を用いた視覚表現では、色相（Ｈ）は色環の三次元円錐状の構造に描かれ、彩度（Ｓ）はその円錐の中心からの距離、明度（Ｖ）は円錐の頂点からの距離で表される。瞳孔領域画像（ＲＧＢ画像）におけるＲ，Ｇ，Ｂは、それぞれ０．０から１．０の範囲にある。(Ｒ，Ｇ，Ｂ)で定義された色が与えられたとすると、それに相当する(Ｈ，Ｓ，Ｖ)のカラーは数式により導出できる。一方、色相（Ｈ）は０．０から３６０．０まで変化し、色相（Ｈ）が示された色環に沿った角度で表現される。彩度（Ｓ）および明度（Ｖ）は０．０から１．０までの範囲で変化する。
【００３７】
ＲＧＢ画像（ＲＧＢ色空間）では、明るさと色味が混ざって表現されているので、例えば照明の微妙な変化によって明るさが変化した場合、ＲＧＢすべてに影響が及ぶことになる。これに対して、ＨＳＢ色空間では、色相（Ｈ），彩度（Ｓ）および明度（Ｖ）に分解されているので、瞳孔の色みの変化を正しく検出することが可能である。
色特徴を検出し、水晶体全体が白くなっているか否かを判定することにより、核白内障や後嚢下白内障を判別する。また、色特徴から模様を検出し、水晶体内部の混濁度合いを定量化することにより、皮質白内障を判別する。なお、模様の検出の際、後述するように、まつげ除去を行ってノイズ除去を図る。
【００３８】
極座標展開画像の解析は、図８に示すように、３つの特徴量に着目して解析する。図８において、横方向のラインの特徴であるｖ、縦方向のラインの特徴であるｈ、そしてｖ−ｈの特徴を３つの特徴量として解析する。
ここで、特徴量ｖ、ｈを導出するために、図８に示す横方向および縦方向の画像フィルタを用いる。極座標展開画像において、ある点を中心に画像フィルタを適用して、画像フィルタで白画素になっている画素の明度の平均値 ｍ１と黒画素になっている画素の明度平均値 ｍ２ を計算して、最後に、ｍ１−ｍ２ を計算して特徴量ｖ，ｈ とする。
【００３９】
図８の左の極座標展開画像において、四角枠で囲んでいる部分の場合、横方向の縞模様が多く検出できるため、特徴量ｖは大きな値をとり、反対に、縦方向の縞模様はあまり検出されないため、特徴量ｈは特徴量ｖほど大きな値にはならないと期待される。
このように、特徴量ｖの大きさと特徴量ｖ−ｈを用いて、候補点を算出する。
そして、候補点に算出された画素の画素数、すなわち、リング状反射光と判定された画素の画素数と、それぞれの画素に対して横方向の画像フィルタで計算した結果の平均値を用いて、反射光の明瞭度合いを数値化し、これをリング状反射光の強さを示す特徴量とする。
【００４０】
次に、リング状反射光の検出方法について説明する。
リング状反射光は、瞳孔の中心と同心円状に現れる。そのため、瞳孔領域画像を極座標展開画像に展開し、横方向の縞模様としてリング状反射光を検出する。すなわち、図９に示すように、極座標展開画像から、それぞれの点が円を構成しているか否かを検証し、各候補点に対応する円の中心を「投票」し、矛盾の無いものを検出する。
具体的には、以下の（処理１）〜（処理４）を行う。
【００４１】
（処理１）先ず、円の中心位置と半径を投票するための投票空間 Ｍ（ｘ，ｙ，ｒ）を準備し、Ｍ（ｘ，ｙ，ｒ）＝０  に初期化する。Ｍの値がスコアとなり、スコアが大きいほど、そこにリング状反射光が存在することになる。
【００４２】
（処理２）予め検出された候補点（図９における曲線）について、下記の（処理２−１）および（処理２−２）を実行する。
（処理２−１) 極座標 Ｐ(r, θ) を平面座標 Ｉ(x, y) に変換する。
（処理２−２) 予め決めておいた範囲 [ｒ１,ｒ２] で半径 r を変化させて、（処理２−２−１)および（処理２−２−２)を実行する。
（処理２−２−１) Ｉから距離ｒにある点を算出する。
（処理２−２−２) その点においてスコア Ｍ（ｘ，ｙ，ｒ）に１を加算する。
【００４３】
（処理３）以上の処理を繰り返した後、Ｍ（ｘ，ｙ，ｒ）が最大になる（ｘ，ｙ，ｒ）＝（ｘｍ, ｙｍ, ｒｍ) を算出する。
（処理４) (ｘｍ, ｙｍ, ｒｍ) に投票した候補点を「矛盾のないもの」とする。
【００４４】
図１０（ａ）〜（ｆ）は、様々な極座標展開画像の例を示している。図１０（ａ），（ｃ），（ｅ），（ｆ）は、縦方向に細い線が多数見られるが、これは「まつげ」の反射像が写りこんでいるからである。
この「まつげ」の反射像は、瞳孔領域の混濁を検出する際のノイズになるため、解析では、これを画像から除外する。
【００４５】
図１１を参照して、まつげと混濁の識別について説明する。図１１（１）は、「まつげ」の反射像が写りこんでいる瞳孔領域画像の極座標展開画像の一例を示している。瞳孔領域画像のブロック単位で、高速フーリエ変換（ＦＦＴ）を用いて、模様が存在している領域を検出する。ここでは、６４×６４の画像ブロックに分割し、横方向成分と縦方向成分に分けて高速フーリエ変換（ＦＦＴ）を行う。その際、直流成分は平均値（平均明度）に相当するので無視することにし、また、極座標系に展開したときの横方向の高周波成分は、まつげの写り込みの影響を強く受けるので無視することにする。そして、残りの成分の平均値が閾値より大きければ混濁が存在する候補領域とする。
図１１（２）では、上記の処理で混濁が存在する候補になったブロック部分を示している。そして、混濁が存在する候補になったブロック部分それぞれに対して二値化処理を施す。
【００４６】
まつげの誤検出を回避すべく、さらに以下の処理を行って、まつげの除去を行う。
１）画像の中でまつげが存在しうる領域を予め設定しておき、二値化結果を抽出する。
２）まつげは点在しているものではなく、連続していることから、連結領域を導出する。
３）連結領域で、「同じモーメントをもつ楕円の偏心率（ａ）」と「その楕円の角度（ｂ）」を求める。
４）ａ＞ ａ０、かつ、ｂ＞ｂ０であれば、まつげとして除去する。ここで、ａ０，ｂ０はまつげの形状の判定閾値である。
このように、瞳孔領域画像のブロック単位で高速フーリエ変換（ＦＦＴ）を用いて、模様が存在している領域を検出し、二値化結果の形状特徴から「まつげ」と「混濁」を識別する。
【００４７】
「混濁」に関しては、瞳孔領域のどの部分に混濁が現れているかによって白内障の疑いの有無や白内障の程度が異なるため、瞳孔領域の中心からの距離に応じて混濁度合いを調べることにする。例えば、図１２（１）に示す瞳孔領域のような模様に関して、図１２（２）のように瞳孔領域の中心からの距離に応じてドーナツ状に５等分（Ｐ１〜Ｐ５）に分けて、「混濁度合い」を評価する。
【００４８】
そして、色特徴のパラメータを、Ｈ，Ｓ，Ｖのそれぞれの平均値、最大値、中央値、すなわち、３×３の９個のパラメータ（ｆ１〜ｆ９）とする。また、混濁度合いのパラメータを、瞳孔領域の中心からの距離に応じてドーナツ状に５等分（Ｐ１〜Ｐ５）した５個のパラメータ（ｆ１０〜ｆ１４）とする。リング状反射光の画素数（ｒ_１）とフィルタ応答値とこれらの色特徴のパラメータおよび混濁度合いのパラメータを一つの特徴ベクトルｖ＝（ｆ１，ｆ２，ｆ３，・・・，ｆ１６）とする（図１３を参照）。
【００４９】
そして、上記の特徴ベクトルｖ＝（ｆ１，ｆ２，ｆ３，・・・，ｆＮ）と眼科専門医による画像診断の結果を組にして、学習データとしてデータベース化していき、サーポートベクターマシン（ＳＶＭ）などの機械学習手法を用いて判定規則を学習し、白内障眼と健常眼を判別する。
【００５０】
交差検定による評価方法は、機械学習における識別器のモデルを評価する手法の一つであり、与えられたデータを分割し、その一部をまず学習データとして利用して、残る部分でその識別器のテストを行い、学習の妥当性の検証・確認を行うもので、分類器（識別器）がどれだけ本当に母集団に対処できるかを良い近似で検証・確認を行う。
交差検定による評価方法のうち、Ｋ−分割交差検定では、標本群をＫ個に分割し、そのうちの１つをテスト事例とし、残る Ｋ−１個を学習事例とし、組み合わせを変更してＫ回評価を繰り返す。図１５は、１０分割交差検定による評価を示している。
【００５１】
下記表２〜４は、実施例１の白内障検査装置を用いて、健常者と白内障患者について判定・評価した結果を示す。ここで表では、白内障の眼のうち正しく判定できた枚数（Ｔｒｕｅ−Ｐｏｓｉｔｉｖｅ）、白内障の眼を健常な眼と誤判定した枚数（Ｆａｌｓｅ−Ｐｏｓｉｔｉｖｅ）を示している。ここで、白内障と判定したものには、核白内障と皮質白内障の双方を含んでいる。また、健常な眼のうち正しく判定できた枚数（Ｔｒｕｅ−Ｎｅｇａｔｉｖｅ）、健常眼を白内障眼と誤判定した枚数（Ｆａｌｓｅ−Ｎｅｇａｔｉｖｅ）を示している。
【００５２】
下表２〜４において、「Ｒｅｃａｌｌ」は再現率であり、実際に正であるもののうち、正であると予測されたものの割合を意味し、見つけなければいけないものをどれだけ見つけたかを数値化したものであり、最小値が０で、１あるいは１に近いほど優れている。また、「Ｐｒｅｃｉｓｉｏｎ」は白内障と判定したものの正確さを表す指標であり、ＴＰ／（ＴＰ＋ＦＰ）から算術する。「Ｐｒｅｃｉｓｉｏｎ」も、最小値が０で、１に近いほど優れている。また、「Ｆ値（Ｆ−ｎｕｍｂｅｒ）」は、機械学習における予測結果の評価尺度の一つであり、精度と再現率の調和平均を示している。具体的には、Ｆ値（Ｆ−ｎｕｍｂｅｒ）は、Ｒｅｃａｌｌの値を“Ａ”、Ｐｒｅｃｉｓｉｏｎの値を“Ｂ”とすると、２×Ｂ×Ａ／（Ａ＋Ｂ）から算術する。「Ｆ値（Ｆ−ｎｕｍｂｅｒ）」も、最小値が０で、１に近いほど優れている。
【００５３】
表２は、サポートベクターマシーン（ＳＶＭ）を用いて、色特徴と混濁の解析を行った結果（リング状フラッシュの反射像を特徴量に含めない場合）を示している。一方、表３は、リング状反射光を観測して、サポートベクターマシーン（ＳＶＭ）を用いて、色特徴と混濁の解析を行った結果を示している。また、表４は、リング状フラッシュを使用しない場合とリング状フラッシュを使用した場合の「Ｒｅｃａｌｌ」，「Ｐｒｅｃｉｓｉｏｎ」および「Ｆ値（Ｆ−ｎｕｍｂｅｒ）」を比較したものである。
【００５４】
ここで、表２〜表４の左列に示しているＮＯ閾値およびＣ閾値は学習における教師データを与える際に用いた白内障発症の有無を判断する境界である。先ず、評価に用いた画像すべてについて、白内障分類に熟練した視能訓練士によって、ＬＯＣＳIII（The Lens Opacities Classification System III）白内障分類に沿って、ＮＯおよびＣの判定を行った。ＮＯは水晶体の濁りの程度を示す指標であり、Ｃは皮質白内障の程度を示す指標である。それぞれ値が大きい方ほど、程度が高いことを表す。
これらの指標は主観的なものであり、特定の閾値より大きければ白内障とするという統一的な基準はないため、概ね妥当と判断できる範囲、例えば、ＮＯについては２から３、Ｃについて２から２．４に変化させて、それぞれの値を境界値として学習データを作成しテストを行った。
【００５５】
また、サポートベクターマシン（ＳＶＭ）の学習について説明する。サポートベクターマシンは非常に精度が高い識別器として知られているが、十分な性能を得るためにはパラメータを慎重に設定する必要がある。ここでは広く用いられているラジアル基底関数（ＲＢＦ；radial basis function）カーネルを使用し、そのパラメータを様々に設定して最良の結果となるものを示した。
【００５６】
表４から、全ての閾値について、リング状フラッシュを用いると「Ｒｅｃａｌｌ」，「Ｐｒｅｃｉｓｉｏｎ」および「Ｆ値（Ｆ−ｎｕｍｂｅｒ）」の結果が向上していることが確認できた。すなわち、リング状フラッシュを使用し、リング状反射光を観測することにより、白内障の判定精度が向上することが示された。
【００５７】
【表２】
2015216939.tif 000004
【００５８】
【表３】
2015216939.tif 000005
【００５９】
【表４】
2015216939.tif 000006
【００６０】
下記表５，６に、決定木を使用した識別を行った結果を示す。決定木を用いることにより、色模様，混濁またはリング状フラッシュの内、どの特徴がどのように効いているかがわかる。表５，６から、リング状フラッシュを使った方がほぼすべての場合で、Ｆ値（Ｆ−ｎｕｍｂｅｒ） が向上していることがわかった。
また、表５，６において、最良のケース（グレーで示した行）を決定木で可視化した。図１６，図１７に示すように、決定木を可視化したところ、色特徴を使った識別が主に動作し、その他の情報を用いた識別がそれを補う働きをしていることがわかった。
【００６１】
なお、可視化の際には交差検定は用いずにすべてのデータを使って学習している。また、図１６，１７において、決定木の節点（ノード）は丸枠が対応し、決定木の葉（識別結果）には四角枠が対応している。入力データは最上部に位置する節点（根ノード）に入力され、それぞれの節点での判定条件にしたがってどの子節点に進むかが決定される。丸枠の模様がどの特徴量を使って判定しているかを示している。
【００６２】
【表５】
2015216939.tif 000007
【００６３】
【表６】
2015216939.tif 000008
【００６４】
上述した実施例１の白内障検査装置では、カメラ２１とリング状フラッシュ２２とコンピュータ２３で簡単に構成できるものであったが、カメラ２１とリング状フラッシュ２２とコンピュータ２３のハードウェアが既に存在する場合に、白内障判定プログラムが提供されることにより、カメラとリング状フラッシュとコンピュータのハードウェアが白内障検査装置として作動可能となる。
【００６５】
白内障判定プログラムは、図１に示すように、画像情報取得手段２３ａと反射光検出手段２３ｂと色特徴解析手段２３ｃと混濁解析手段２３ｄとリング状反射光個数解析手段２３ｅと白内障判定手段２３ｆから構成される。画像情報取得手段３１は、カメラ２１で撮影した画像を取り込むものである。例えば、汎用のデジタルカメラはコンピュータと接続可能なようにＵＳＢなどのデータ通信インタフェース２４を有している。画像情報取得手段２３ａは、カメラとデータ通信を行い、カメラ内の画像情報を取り込み、画像情報をコンピュータ内で解析できるようにするものである。
【産業上の利用可能性】
【００６６】
本発明は、眼球の医療診断・検査（スクリーニング）装置として有用である。
【符号の説明】
【００６７】
１    眼球
２    水晶体
３    虹彩
４    網膜
５    瞳孔
６    角膜
７    ハウジング
１０〜１２  リング状反射光
２０  リング状投射光
２１  カメラ
２１ａ  カメラレンズ
２１ｂ  レンズ本体
２１ｃ  ＣＣＤイメージセンサ
２２  リング状フラッシュ
２２ａ  リング状スリット
２２ｂ  遮蔽板
２２ｃ  マウント部
２２ｄ  発光部
２３  コンピュータ
２４  データ通信インタフェース
２５  バッテリーユニット
２６  電源ケーブル
【図１】
2015216939.tif 000009
【図２】
2015216939.tif 000010
【図３】
2015216939.tif 000011
【図４】
2015216939.tif 000012
【図５】
2015216939.tif 000013
【図６】
2015216939.tif 000014
【図７】
2015216939.tif 000015
【図８】
2015216939.tif 000016
【図９】
2015216939.tif 000017
【図１０】
2015216939.tif 000018
【図１１】
2015216939.tif 000019
【図１２】
2015216939.tif 000020
【図１３】
2015216939.tif 000021
【図１４】
2015216939.tif 000022
【図１５】
2015216939.tif 000023
【図１６】
2015216939.tif 000024
【図１７】
2015216939.tif 000025