JP4289973B2 - 細隙灯顕微鏡から取得された白内障被検眼のカラー画像データに基づき、白内障の症状を定量化するための色補正装置 - Google Patents

Description

本発明は、撮像手段を備えた細隙灯顕微鏡により取得した白内障被検眼のカラー画像データに基づき、白内障の症状を定量化するための色補正装置に関するものである。

眼科の分野において、同一患者の被検眼のカラー画像データを時系列的に記録し、電子ファイリングやコンピュータによる診断支援等を行うために、撮像手段を用いてカラー画像データとして取り込むようにした細隙灯顕微鏡等の眼科装置が開発されている。しかし、このような眼科装置においては、光源および撮像手段が同一であっても、被検眼像の撮像時の照明条件が異なると、同じ被検眼像を撮像しても異なるカラー画像データが得られてしまう。また、光源や撮像手段が異なると、同一の照明条件下において、同一の被検眼像を撮像する場合であっても、異なるカラー画像データが得られてしまう。このように眼科装置によって取り込まれたカラー画像データにバラツキが生じると、同一患者について記録された一連のカラー画像データに基づいて、症状が進行しているのかあるいは治癒しつつあるのかの判断ができず、正確な診断等が行えない。

このため、例えば、撮像手段と、撮像手段により変換された画像信号をデジタル値に変換する変換手段と、変換手段からの出力を保存する被検眼画像データ記憶手段と、被検眼画像データからアパーチャ領域の画像データを抽出する画像データ抽出手段と、アパーチャ領域の画像データに基づいて補正値を算出する演算手段と、被検眼画像データから補正値を減算する減算手段を備えた眼科撮影装置が提案されている（特許文献１参照）。この装置によれば、被検眼像の撮影時に、光源および撮像手段が同一で、照明条件が異なる場合、並びに、光源および撮像手段が異なり、照明条件が同一の場合のいずれにおいても、色補正されたバラツキのない画像データが得られる。

ところで、白内障眼の水晶体は微妙な色合い（必ずしも白色とは限らない）で濁ってくるので、白内障の症状の進行度や回復状況を的確に判断するためには、同一患者の白内障眼における水晶体の微妙な色の変化が正確かつ的確に把握できなければならない。このため、異なる光源および撮像手段、並びに照明条件で撮像された一連のカラー画像データをできるだけ高い精度で色補正し、同一の条件（光源、撮像手段および照明条件）下に撮像されたカラー画像データに変換しなければならない。

しかしながら、上述の従来の装置では、カラー画像データの黒レベルのみを基準としてそのカラー画像データの全体を色補正するので、時間を追って取得された一連の画像データを同一の条件（光源、撮像手段および照明条件）下に撮像されたカラー画像データに変換する精度はそれほど高くない。したがって、この装置では、同一患者の被検眼の一連の画像データに基づいて、白内障の症状の進行度や回復状況の正確な判断が可能な程度の精度で、白内障の症状を定量化することはできない。
特開２００３−２１０４１２号公報

したがって、本発明の課題は、撮像手段によって撮像された、細隙灯顕微鏡による白内障眼の像のカラー画像データを用いて、高い精度で白内障の症状の定量化を行うことができるようにすることにある。

上記課題を解決するため、本発明は、撮像手段を備えた細隙灯顕微鏡により取得された白内障被検眼のカラー画像データに基づき、白内障の症状を定量化するための色補正装置であって、細隙灯顕微鏡に備えられた、線形な感度特性を有する撮像手段によって撮像され、取得されたカラーチャートのカラー画像データのカラーパッチ毎の生の色度ベクトル値（Ｐ_1ｉ）_ｊを求める第１色度ベクトル値検出部と、既知の色特性を有する所定の光源、および既知の線形な感度特性を有する所定の撮像手段を使用して前記カラーチャートを撮像し、取得した前記カラーチャートのカラー画像データに基づいて求められた、前記カラーチャートのカラーパッチ毎の初期参照色度ベクトル値（Ｐ_0ｉ）_ｊを格納する初期参照色度ベクトル値格納部と、前記第１色度ベクトル値検出部によって求められた前記生の色度ベクトル値（Ｐ_1ｉ）_ｊと、前記初期参照色度ベクトル値格納部に格納された前記初期参照色度ベクトル値（Ｐ_0ｉ）_ｊとから、前記細隙灯顕微鏡の光源の色特性Ｆと、前記細隙灯顕微鏡に備えられた撮像手段の感度特性Ｓの積からなるマトリクスＡ_ｊｋ＝Ｆ×Ｓを、

を用いて求める第１演算部と、前記細隙灯顕微鏡に備えられた撮像手段によって撮像され、取得された白内障被検眼のカラー画像データを構成する各画素の生の色度ベクトル値（Ｑ_0ｉ）_ｊを求める第２色度ベクトル値検出部と、前記第２色度ベクトル値検出部によって求められた前記生の色度ベクトル値（Ｑ_0ｉ）_ｊに対する補正色度ベクトル値（ＣＱ_0ｉ）_ｊを、前記第１演算部によって求められたマトリクスＡ_ｊｋの値、および

を用いて求める第２演算部と、前記第２演算部によって求められた補正色度ベクトル値（ＣＱ_0ｉ）_ｊに基づいて、前記白内障被検眼のカラー画像データを補正する第１画像データ補正部と、前記第１画像データ補正部によって補正された前記白内障被検眼のカラー画像データにおける虹彩部分の補正色度ベクトル値の平均値（Ｗ_ａｖ）_ｊを求める第３色度ベクトル値検出部と、前記白内障被検眼のカラー画像データが第１回目に取得されたとき、前記第３色度ベクトル値検出部によって求められた平均値（Ｗ_ａｖ）_ｊを補正参照色度ベクトル値（Ｗ₀）_ｊとして記録するとともに、前記第１画像データ補正部によって補正された前記白内障被検眼のカラー画像データを記録するが、前記白内障被検眼のカラー画像データが第２回目以降に取得されたときには、前記第３色度ベクトル値検出部によって求められた平均値（Ｗ_ａｖ）_ｊと前記補正参照色度ベクトル値（Ｗ₀）_ｊとの平均値を求めて新たな補正参照色度ベクトル値（Ｗ₀）_ｊとし、さらに、前記平均値（Ｗ_ａｖ）_ｊと前記新たな補正参照色度ベクトル値（Ｗ₀）_ｊとの差を求め、その差に基づいて、前記第１画像データ補正部によって補正された前記白内障被検眼のカラー画像データをさらに補正し、それを記録する第２画像データ補正部とを有していることを特徴とする装置を構成したものである。

上記構成において、色度ベクトル値としては、例えば、ＲＧＢ値、ＸＹＺ値、Ｌａｂ値またはＬ^＊ａ^＊ｂ^＊値が挙げられる。また、撮像手段としては、例えば、ＣＣＤまたはＣＭＯＳが挙げられる。

本発明によれば、光源および撮像手段を備えた細隙灯顕微鏡からなる光学系を、白内障眼のカラー画像データの撮像に使用する前に、予めカラーチャートを用いて色補正し、さらに、この色補正を行った光学系を用いて白内障眼のカラー画像データを取得する際に、、既に取得された同一の患者の白内障眼の画像データの虹彩部分の色度ベクトル値の平均値を基準として、取得された画像データの色補正を行うようにしたので、異なる条件（光源、撮像手段および照明条件）下に取得された同一患者の白内障眼の一連の画像データを、非常に高い精度で、同一の条件下に撮像されたカラー画像データに変換することができ、同一患者の被検眼の一連の画像データに基づいて、白内障の症状を定量化し、それによって白内障の症状の進行度や回復状況を正確に判断することができる。

以下、添付図面を参照して本発明の好ましい実施例について説明する。図１は、本発明の１実施例による方法のフロー図である。図１を参照して、本発明によれば、まず最初、既知の色特性を有する所定の光源および既知の線形な感度特性を有する所定の撮像手段を使用してカラーチャートが撮像され、取得されたカラーチャートのカラー画像データのカラーパッチ毎の初期参照色度ベクトル値（Ｐ_0ｉ）_ｊが求められる（図１のステップ（ａ）参照）。

次に、ステップ（ａ）で使用されたものと同一のカラーチャートが細隙灯顕微鏡によって観察され、その観察像が、細隙灯顕微鏡に備えられた、線形な感度特性を有する撮像手段によって撮像され、取得されたカラーチャートのカラー画像データのカラーパッチ毎の生の色度ベクトル値（Ｐ_1ｉ）_ｊが求められる（図１のステップ（ｂ）参照）。
図５には、本発明の方法において使用するのに適したカラーチャートの１例が示してある。図５を参照して、カラーチャート１０は、中央の細隙灯顕微鏡用カラーチャート部分４１と、細隙灯顕微鏡用カラーチャート部分４１の両側に配置された初期参照色度ベクトル値測定用カラーチャート部分４０ａ、４０ｂとから構成されている。細隙灯顕微鏡用カラーチャート部分４１は、細隙灯（光源）による光照射範囲内に収まるようなサイズ（幅約２ｍｍ）に形成され、また、初期参照色度ベクトル値測定用カラーチャート部分４０ａ、４０ｂには、細隙灯顕微鏡用カラーチャート部分４１に配列されたカラーパッチと同じカラーパッチが配列されている。この場合、細隙灯顕微鏡用カラーチャート部分４１の一部に孔４２を設けることにより、黒色パッチとすることができる。こうして、ステップ（ａ）では、このカラーチャート１０の初期参照色度ベクトル測定用カラーチャート部分４０ａ、４０ｂが使用され、ステップ（ｂ）では、細隙灯顕微鏡用カラーチャート部分４２１が使用される。

次に、ステップ（ａ）で求められたカラーチャートのカラーパッチ毎の初期参照色度ベクトル値（Ｐ_0ｉ）_ｊと、ステップ（ｂ）で求められたカラーチャートのカラーパッチ毎の生の色度ベクトル値（Ｐ_1ｉ）_ｊから、細隙灯顕微鏡の光源の色特性Ｆと、細隙灯顕微鏡に備えられた撮像手段の感度特性Ｓの積からなるマトリクスＡ_ｊｋ＝Ｆ×Ｓが、

を用いて求められる（図１のステップ（ｃ）参照）。そして、細隙灯顕微鏡から取得されたカラー画像データをマトリクスＡによって変換することにより、色補正がなされたカラー画像データを得ることができる。

こうして予めカラーチャートによる色補正がなされた細隙灯顕微鏡を用いて、白内障被検眼が第１回目に観察され、観察像が細隙灯顕微鏡に備えられた撮像手段によって撮像され、取得された白内障被検眼のカラー画像データを構成する各画素の生の色度ベクトル値（Ｑ_0ｉ）_ｊから、その補正色度ベクトル値（ＣＱ_0ｉ）_ｊが、

を用いて求められ、補正色度ベクトル値（ＣＱ_0ｉ）_ｊに基づいて白内障被検眼のカラー画像データが補正され、それが記録されるとともに、補正後の白内障被検眼のカラー画像データにおける虹彩部分の補正色度ベクトル値の平均値が求められ、それが補正参照色度ベクトル値（Ｗ₀）_ｊとされる（図１のステップ（ｄ）参照）。

そして、ステップ（ｄ）で使用されたものと同一の細隙灯顕微鏡を用いて、白内障被検眼が第２回目以後に観察されるたびに、観察像が細隙灯顕微鏡に備えられた撮像手段によって撮像され、取得された白内障被検眼のカラー画像データを構成する各画素の生の色度ベクトル値（Ｑ_1ｉ）_ｊから、その補正色度ベクトル値（ＣＱ_1ｉ）_ｊが、

を用いて求められ、補正色度ベクトル値（ＣＱ_1ｉ）_ｊに基づいて白内障被検眼のカラー画像データが補正され、補正後の白内障被検眼のカラー画像データにおける虹彩部分の補正色度ベクトル値の平均値（Ｗ_ａｖ）_ｊが求められ、その平均値（Ｗ_ａｖ）_ｊと補正参照色度ベクトル値（Ｗ₀）_ｊとの平均値が求められて新たな補正参照色度ベクトル値（Ｗ₀）_ｊとされ、さらに、平均値（Ｗ_ａｖ）_ｊと新たな補正参照色度ベクトル値（Ｗ₀）_ｊとの差が求められ、その差に基づいて、補正後の白内障被検眼のカラー画像データがさらに補正され、それが記録される（図１のステップ（ｅ）参照）。

本発明の方法による、カラーチャートを用いてマトリクスＡ_ｊｋを求める一連の手順（上記ステップ（ａ）〜（ｃ））の理論的根拠について少し説明すると、各波長毎の光の強度を表したものを分光分布といいλを波長として、波長の関数ｓ（λ）で表わされる。今、３８０ｎｍ〜７８０ｎｍの波長領域を１０ｎｍ毎に区分し、ｊ＝０〜４０として、

と表すことにより分光分布ｓを４１次元のベクトルで表すことができる。すなわち、
ｓ＝（ｓ_０，ｓ_１，ｓ_２，・・・，ｓ_４０）

人間の知覚する物体の色は、光源の分光分布、物体の分光反射率、人間の分光感度により決定される。
今、ある物体の分光反射率のベクトル表現をＢ、光源の分光分布のベクトル表現をＦ、人間の各錐体の分光感度のベクトル表現をそれぞれｐ、ｄ、ｔとおくと三刺激値Ｐ、Ｄ、Ｔは

となる。
また物体から反射して目に入ってくる光の分光エネルギーのベクトル表現をＯとし、分光感度をまとめて

とすると

と内積の形に表現できる。
人間の目にとっては、光の分光分布が異なっても三刺激値が同一であれば同じ色として感じられる。

全ての色光は、独立した３種の単色光で等色（人間の目で等しいと感じられる）できることが知られている。ここで、「独立した」とは、ある色が他の色によって等色されないことを意味する。
今、ＣＩＥの標準により
赤の原色光Ｒとしてλ＝７００ｎｍの単色光、
緑の原色光Ｇとしてλ＝５４６．１ｎｍの単色光、
青の原色光Ｂとしてλ＝４３５．８ｎｍの単色光、
を三原色として、３８０ｎｍ〜７８０ｎｍの範囲内の各波長の単色光を三原色で等色した際のＲＧＢの分光分布、等色関数と呼ばれる。
等色関数は三種の錐体の分光感度の線形変換となっている。

ここで、Aはマトリクスをあらわす。

今、２４色のカラーパッチをＤ６５光源のもと、等色関数でセンスされた撮像手段で撮像した際のＲＧＢ値を（Ｐ_0ｉ）_ｊとする。ここで、ｉは各パッチの色を示す添字で０〜２３となり、ｊは色度を示すベクトルの次元を表し、ここではＲＧＢの３次元であるため０〜２となる。
また、未知の光源で未知の感度特性を有する撮像手段によりカラーパッチを撮像した際のＲＧＢ値を同様に（Ｐ_1ｉ）_ｋとする。ここで、双方とも色空間を同じＲＧＢとしているためｊ＝ｋである。
（Ｐ_0ｉ）_ｊと（Ｐ_1ｉ）_ｋとの関係を、（Ｐ_0ｉ）_ｊ＝Ａ_ｉｊｋ×（Ｐ_1ｉ）_ｋとし、未知の感度特性が等色関数でセンスされたものであると仮定すると、（Ｐ_0ｉ）_ｊと（Ｐ_1ｉ）_ｋとは線形の関係にあり、Ａ_０ｊｋ＝Ａ_１ｊｋ＝Ａ_２ｊｋ＝・・・＝Ａ_２３ｊｋとなり、最小二乗法によりパラメータＡ_ｊｋを求めることができる。

上述の方法を具体例を用いて説明すると、ステップ（ａ）では、所定の光源としてＤ６５を用い、また、所定の撮像手段として等色関数によってセンスされたＣＣＤを用いて、２４個の異なるカラーパッチを備えたカラーチャートを撮像して取得されたカラー画像データのカラーパッチ毎の初期参照ＲＧＢ値［（Ｐ_0ｉ）_１＝Ｒ_ｉ，（Ｐ_0ｉ）_２＝Ｇ_ｉ，（Ｐ_0ｉ）_３＝Ｂ_ｉ，ｉ＝１〜２４］、すなわち、
［（Ｒ_１，Ｇ_１，Ｂ_１）（Ｒ_２，Ｇ_２，Ｂ_２）・・・（Ｒ_２４，Ｇ_２４，Ｂ_２４）］
が求められる。

ステップ（ｂ）では、細隙灯顕微鏡の撮像手段（ＣＣＤ）によって取得されたステップ（ａ）と同一のカラーチャートのカラー画像データのカラーパッチ毎の生のＲＧＢ値［（Ｐ_0ｉ）_１＝ｒ_ｉ，（Ｐ_0ｉ）_２＝ｇ_ｉ，（Ｐ_0ｉ）_３＝ｂ_ｉ，ｉ＝１〜２４］、すなわち、
［（ｒ_１，ｇ_１，ｂ_１）（ｒ_２，ｇ_２，ｂ_２）・・・（ｒ_２４，ｇ_２４，ｂ_２４）］
が求められる。

細隙灯顕微鏡のＣＣＤの分光感度が等色関数によってセンスされていると仮定すると、３次元空間において線形に写像されているものと考えられるから、

として、
Ｙ＝ＡＸ
が成立し、これから最小二乗法を用いて、
Ａ（ＸＸ^Ｔ）＝ＹＸ^Ｔ、
Ａ＝ＹＸ^Ｔ（Ｘ^ＴＸ）^−１
から、マトリクスＡが求められる。

こうして、予め色補正がなされた細隙灯顕微鏡を用いて、白内障被検眼が第１回目に観察され、観察像が細隙灯顕微鏡に備えられた撮像手段によって撮像され、取得された白内障被検眼のカラー画像データを構成する各画素の生のＲＧＢ値（ｒ_ｉ，ｇ_ｉ，ｂ_ｉ）から、その補正ＲＢＧ値（Ｒ_ｉ，Ｇ_ｉ，Ｂ_ｉ）が、
Ｒ_ｉ＝ａｒ_ｉ＋ｄｇ_ｉ＋ｇｂ_ｉ
Ｇ_ｉ＝ｂｒ_ｉ＋ｅｇ_ｉ＋ｈｂ_ｉ
Ｂ_ｉ＝ｃｒ_ｉ＋ｆｇ_ｉ＋ｉｂ_ｉ
を用いて求められ、補正ＲＧＢ値（Ｒ_ｉ，Ｇ_ｉ，Ｂ_ｉ）に基づいて白内障被検眼のカラー画像データが補正され、それが記録される。また、補正後の白内障被検眼のカラー画像データにおける虹彩部分の補正ＲＧＢ値（Ｒ_ｉ，Ｇ_ｉ，Ｂ_ｉ）の平均値が求められ、それが補正参照ＲＧＢ値（Ｗ_Ｒ，Ｗ_Ｇ，Ｗ_Ｂ）とされる。

この場合、補正後の白内障被検眼のカラー画像データにおける虹彩部分の補正ＲＧＢ値の平均値は次のようにして求められる。図３は、同一の患者の白内障被検眼のカラー画像の一例を示した平面図であり、（Ａ）は第１回目に撮像された画像を示し、（Ｂ）は第２回目に撮像された画像を示している。図３（Ａ）において、２０は白内障被検眼であり、２１は虹彩であり、２２は瞳孔であり、２３は白内障によって白濁化した水晶体部分である。図３（Ｂ）において、３０は白内障被検眼であり、３１は虹彩であり、３２は瞳孔であり、３３は白内障によって白濁化した水晶体部分である。
今、適当なディスプレイに、図３（Ａ）に示された白内障被検眼のカラー画像２０が表示され、そして、この画像２０における虹彩２１の適当な領域（例えば、図３（Ａ）領域のＭ_１）が、医師により適当なポインティング・デバイスを用いて指定される。そして、指定された領域Ｍ_１を構成する画素のＲＧＢ値の平均値が求められ、求められた平均ＲＧＢ値が補正参照ＲＧＢ値（Ｗ_Ｒ，Ｗ_Ｇ，Ｗ_Ｂ）とされる。なお、この例では、画像２０中の虹彩２１における適当な領域を医師が指定するようにしているが、この画像データに基づき、コンピュータを用いて自動的に虹彩部分を検出し、その画素のＲＧＢ値の平均値を求めるようにすることもできる。

次に、同一の細隙灯顕微鏡を用いて、同一患者の白内障被検眼が第２回目に観察され、第１回目の場合と同様にして、取得された白内障被検眼のカラー画像データを構成する各画素の生のＲＧＢ値からその補正ＲＧＢ値が求められ、補正ＲＧＢ値に基づいて白内障被検眼のカラー画像データが補正される。そして、ディスプレイに、図３（Ｂ）に示された白内障被検眼のカラー画像３０が表示され、この画像３０における虹彩３１の適当な領域（例えば、図３（Ｂ）の領域Ｍ_２）が、医師により適当なポインティング・デバイスを用いて指定される。そして、指定された領域Ｍ_２を構成する画素のＲＧＢ値の平均値が求められ、さらに、その平均値と補正参照ＲＧＢ値（Ｗ_Ｒ，Ｗ_Ｇ，Ｗ_Ｂ）との平均値が求められて新たな補正参照ＲＧＢ値とされ、さらに、前記平均値と新たな補正参照ＲＧＢ値との差が求められ、その差に基づいて、補正後の白内障被検眼のカラー画像データがさらに補正され、それが記録される。同様にして、同一の細隙灯顕微鏡による同一患者の白内障被検眼の第３回目以降の観察がなされ、取得されたカラー画像データが色補正され、色補正後のカラー画像データが記録される。

こうして、本発明の方法によれば、光源および撮像手段を備えた細隙灯顕微鏡からなる光学系を予めカラーチャートを用いて色補正し、この色補正がなされた光学系を用いて白内障眼のカラー画像データを取得する際に、既に取得された同一の患者の白内障眼の画像データの虹彩部分の色度ベクトル値の平均値を基準として、取得された画像データの色補正を行うようにしたので、異なる光源、撮像手段により、または異なる照明条件の下に取得された同一患者の白内障眼の一連の画像データを、非常に高い精度で、同一の光源、撮像手段によって、同一の照明条件の下に撮像されたカラー画像データに変換することができ、同一患者の被検眼の一連の画像データに基づいて白内障の症状を定量化し、白内障の症状の進行度や回復状況を正確に判断することができる。

図２は、本発明の１実施例による色補正装置の構成を示すブロック図である。図２を参照して、番号１は本発明の色補正装置を示しており、第１色度ベクトル値検出部３は、細隙灯顕微鏡１２に備えられた、線形な感度特性を有する撮像手段１４によって撮像され、取得されたカラーチャート１０のカラー画像データのカラーパッチ毎の生の色度ベクトル値（Ｐ_1ｉ）_ｊを求める。初期参照色度ベクトル値格納部２は、既知の色特性を有する所定の光源、および既知の線形な感度特性を有する所定の撮像手段を使用してカラーチャート１０を撮像し、取得したカラーチャート１０のカラー画像データに基づいて求められた、カラーチャートのカラーパッチ毎の初期参照色度ベクトル値（Ｐ_0ｉ）_ｊを格納する。
第１演算部４は、第１色度ベクトル値検出部３によって求められた生の色度ベクトル値（Ｐ_1ｉ）_ｊと、初期参照色度ベクトル値格納部２に格納された初期参照色度ベクトル値（Ｐ_0ｉ）_ｊとから、細隙灯顕微鏡１２の光源１３の色特性Ｆと、細隙灯顕微鏡１２に備えられた撮像手段１４の感度特性Ｓの積からなるマトリクスＡ_ｊｋ＝Ｆ×Ｓを、

を用いて求める。

第２演算部６は、細隙灯顕微鏡１２に備えられた撮像手段１４によって撮像され、取得された白内障被検眼１１のカラー画像データを構成する各画素の生の色度ベクトル値（Ｑ_0ｉ）_ｊを求める第２色度ベクトル値検出部５と、第２色度ベクトル値検出部５によって求められた生の色度ベクトル値（Ｑ_0ｉ）_ｊに対する補正色度ベクトル値（ＣＱ_0ｉ）_ｊを、第１演算部４によって求められたマトリクスＡ_ｊｋの値、および

を用いて求める。また、第１画像データ補正部７は、第２演算部６によって求められた補正色度ベクトル値（ＣＱ_0ｉ）_ｊに基づいて、白内障被検眼１１のカラー画像データを補正する。

第３色度ベクトル値検出部８は、第１画像データ補正部７によって補正された白内障被検眼１１のカラー画像データにおける虹彩部分の補正色度ベクトル値の平均値（Ｗ_ａｖ）_ｊを求める。第２画像データ補正部９は、白内障被検眼１１のカラー画像データが第１回目に取得されたとき、第３色度ベクトル値検出部８によって求められた平均値（Ｗ_ａｖ）_ｊを補正参照色度ベクトル値（Ｗ₀）_ｊとして記録するとともに、第１画像データ補正部７によって補正された白内障被検眼１１のカラー画像データを記録するが、白内障被検眼１１のカラー画像データが第２回目以降に取得されたときには、第３色度ベクトル値検出部８によって求められた平均値（Ｗ_ａｖ）_ｊと補正参照色度ベクトル値（Ｗ₀）_ｊとの平均値を求めて新たな補正参照色度ベクトル値（Ｗ₀）_ｊとし、さらに、前記平均値（Ｗ_ａｖ）_ｊと前記新たな補正参照色度ベクトル値（Ｗ₀）_ｊとの差を求め、その差に基づいて、第１画像データ補正部７によって補正された白内障被検眼１１のカラー画像データをさらに補正し、それを記録する。

図１および図２に示した実施例は、いずれも、細隙灯顕微鏡に備えられた撮像手段が線形な感度特性を有している場合（等色関数でセンスされている場合）に有効である。図４には、細隙灯顕微鏡に備えられた撮像手段が非線形な感度特性を有している場合に有効な方法のフロー図を示したものである。
図４を参照して、まず最初、既知の色特性を有する所定の光源、および既知の線形な感度特性を有する所定の撮像手段を使用してカラーチャートが撮像され、取得されたカラー画像データのカラーパッチ毎の初期参照色度ベクトル値（Ｐ_0ｉ）_ｊが求められる（図４のステップ（ａ’）参照）。

次に、ステップ（ａ’）で使用されたものと同一のカラーチャートが細隙灯顕微鏡によって観察され、その観察像が細隙灯顕微鏡に備えられた、非線形な感度特性を有する撮像手段によって撮像され、取得されたカラー画像データのカラーパッチ毎の生の色度ベクトル値（Ｐ_1ｉ）_ｊが求められる（図４のステップ（ｂ’）参照）。

そして、ステップ（ａ’）で求められた初期参照色度ベクトル値（Ｐ_0ｉ）_ｊと、ステップ（ｂ）で求められた生の色度ベクトル値（Ｐ_1ｉ）_ｊからルックアップ・テーブルが作成される（図４のステップ（ｃ’）参照）。

こうして予め色補正がなされた細隙灯顕微鏡を用いて、白内障被検眼が第１回目に観察され、観察像が細隙灯顕微鏡に備えられた撮像手段によって撮像され、取得されたカラー画像データを構成する各画素の生の色度ベクトル値（Ｑ_0ｉ）_ｊの補正色度ベクトル値（ＣＱ_0ｉ）_ｊが、ルックアップ・テーブルを使用した補間方法を用いて求められる。さらに、補正色度ベクトル値（ＣＱ_0ｉ）_ｊに基づいて白内障被検眼のカラー画像データが補正され、それが記録されるとともに、補正後の白内障被検眼のカラー画像データにおける虹彩部分の補正色度ベクトル値の平均値が求められ、それが補正参照色度ベクトル値（Ｗ₀）_ｊとされる（図４のステップ（ｄ’）参照）。

そして、ステップ（ｄ’）で使用されたものと同一の細隙灯顕微鏡を用いて、白内障被検眼が第２回目以後に観察されるたびに、観察像が細隙灯顕微鏡に備えられた撮像手段によって撮像され、取得されたカラー画像データを構成する各画素の生の色度ベクトル値（Ｑ_1ｉ）_ｊの補正色度ベクトル値（ＣＱ_1ｉ）_ｊが、ルックアップ・テーブルを使用した補間方法を用いて求められる。そして、補正色度ベクトル値（ＣＱ_1ｉ）_ｊに基づいて白内障被検眼のカラー画像データが補正され、補正後の白内障被検眼のカラー画像データにおける虹彩部分の補正色度ベクトル値の平均値（Ｗ_ａｖ）_ｊが求められ、その平均値（Ｗ_ａｖ）_ｊと補正参照色度ベクトル値（Ｗ₀）_ｊとの平均値が求められて新たな補正参照色度ベクトル値（Ｗ₀）_ｊとされる。さらに、平均値（Ｗ_ａｖ）_ｊと新たな補正参照色度ベクトル値（Ｗ₀）_ｊとの差が求められ、その差に基づいて、補正後の白内障被検眼のカラー画像データがさらに補正され、それが記録される（図４のステップ（ｅ’）参照）。

この場合、補間方法としては、線形補間やスプライン補間、あるいは、ガウス関数等を基底関数に用い非線形な関数を近似して求めたい色度を推測する方法等を適用することができる。

本発明の１実施例による色補正方法のフロー図である。 本発明の１実施例による色補正装置のブロック図である。 図１の色補正方法における被検眼の補正カラー画像データ中の虹彩部分のＲＧＢ値の平均値を求めるステップを説明するための図であり、 本発明の別の実施例による色補正方法のフロー図である。 カラーチャートの１例の平面図である。

符号の説明

１ 色補正装置
２ 初期参照色度ベクトル値格納部
３ 第１色度ベクトル値検出部
４ 第１演算部
５ 第２色度ベクトル値検出部
６ 第２演算部
７ 第１画像データ補正部
８ 第３色度ベクトル値検出部
９ 第２画像データ補正部
１０ カラーチャート
１１ 白内障被検眼
１２ 細隙灯顕微鏡
１３ 光源
１４ 撮像手段

Claims (1)

1. 撮像手段を備えた細隙灯顕微鏡により取得された白内障被検眼のカラー画像データに基づき、白内障の症状を定量化するための色補正装置であって、
   細隙灯顕微鏡に備えられた、線形な感度特性を有する撮像手段によって撮像され、取得されたカラーチャートのカラー画像データのカラーパッチ毎の生の色度ベクトル値（Ｐ _1ｉ ） _ｊ を求める第１色度ベクトル値検出部と、
   既知の色特性を有する所定の光源、および既知の線形な感度特性を有する所定の撮像手段を使用して前記カラーチャートを撮像し、取得した前記カラーチャートのカラー画像データに基づいて求められた、前記カラーチャートのカラーパッチ毎の初期参照色度ベクトル値（Ｐ _0ｉ ） _ｊ を格納する初期参照色度ベクトル値格納部と、
   前記第１色度ベクトル値検出部によって求められた前記生の色度ベクトル値（Ｐ _1ｉ ） _ｊ と、前記初期参照色度ベクトル値格納部に格納された前記初期参照色度ベクトル値（Ｐ _0ｉ ） _ｊ とから、前記細隙灯顕微鏡の光源の色特性Ｆと、前記細隙灯顕微鏡に備えられた撮像手段の感度特性Ｓの積からなるマトリクスＡ _ｊｋ ＝Ｆ×Ｓを、
   

   

   を用いて求める第１演算部と、
   前記細隙灯顕微鏡に備えられた撮像手段によって撮像され、取得された白内障被検眼のカラー画像データを構成する各画素の生の色度ベクトル値（Ｑ _0ｉ ） _ｊ を求める第２色度ベクトル値検出部と、
   前記第２色度ベクトル値検出部によって求められた前記生の色度ベクトル値（Ｑ _0ｉ ） _ｊ に対する補正色度ベクトル値（ＣＱ _0ｉ ） _ｊ を、前記第１演算部によって求められたマトリクスＡ _ｊｋ の値、および
   

   

   を用いて求める第２演算部と、
   前記第２演算部によって求められた補正色度ベクトル値（ＣＱ _0ｉ ） _ｊ に基づいて、前記白内障被検眼のカラー画像データを補正する第１画像データ補正部と、
   前記第１画像データ補正部によって補正された前記白内障被検眼のカラー画像データにおける虹彩部分の補正色度ベクトル値の平均値（Ｗ _ａｖ ） _ｊ を求める第３色度ベクトル値検出部と、
   前記白内障被検眼のカラー画像データが第１回目に取得されたとき、前記第３色度ベクトル値検出部によって求められた平均値（Ｗ _ａｖ ） _ｊ を補正参照色度ベクトル値（Ｗ ₀ ） _ｊ として記録するとともに、前記第１画像データ補正部によって補正された前記白内障被検眼のカラー画像データを記録するが、前記白内障被検眼のカラー画像データが第２回目以降に取得されたときには、前記第３色度ベクトル値検出部によって求められた平均値（Ｗ _ａｖ ） _ｊ と前記補正参照色度ベクトル値（Ｗ ₀ ） _ｊ との平均値を求めて新たな補正参照色度ベクトル値（Ｗ ₀ ） _ｊ とし、さらに、前記平均値（Ｗ _ａｖ ） _ｊ と前記新たな補正参照色度ベクトル値（Ｗ ₀ ） _ｊ との差を求め、その差に基づいて、前記第１画像データ補正部によって補正された前記白内障被検眼のカラー画像データをさらに補正し、それを記録する第２画像データ補正部と、を有していることを特徴とする装置。

Images