JP4620443B2 - 眼科装置 - Google Patents

Description

本発明は眼科装置に関する。詳しくは、角膜表面に形成される涙液層の状態を評価することができる眼科装置に関する。

角膜表面や角膜形状を観察・測定する眼科装置としては、ケラトメータや角膜形状解析装置等が知られている。この種の眼科装置では、一般的に、角膜表面に任意のパターンの光を投影し、その投影光の角膜表面からの反射像を観察・撮影・解析することで角膜表面の状態や角膜形状等を測定している。
ところで、角膜表面は通常、涙液等に覆われ、その表面にはなめらかな液状の層が形成されている。この涙液層は、瞬目して開瞼した後、まぶたを開け続けると時間の経過に伴って徐々に厚みが減少し、やがて消滅する。モニタ等を長時間凝視した場合、まばたきの回数が減り、涙液層が消失した状態（いわゆる、ドライアイと呼ばれる乾燥状態）となる。また、モニタ等を凝視しない場合でも、涙液層の消失する速度には個人差があり、短時間で乾燥状態となってしまう人々もいる。このように短時間でドライアイとなる人々は、日常生活において角膜又はまぶたに異物感を感じ、眼科を受診することが多い。
眼科においてドライアイ状態を主訴する人々を診察する場合、涙液層の角膜表面上における厚みを特定し、その時間変化を観察することが重要となる。開瞼後の涙液層の経時的な変化を捉えることで、涙液層の消失状況を知ることができるためである。このため、角膜表面に形成される涙液層の状態を測定するための眼科装置が開発されている（例えば、特許文献１〜３）。

特許文献１に記載の眼科装置は、投光光学系と受光光学系により構成され、被検眼を観察するための細隙灯顕微鏡とＣＣＤカメラ等を備える。この眼科装置では、被検眼に蛍光染色液を点眼し、角膜表面の蛍光反射の濃淡をＣＣＤカメラで撮影する。そして、撮影された蛍光反射の濃淡から涙液および蛍光染色液の混合液によって構成された層の厚みを解析する。この眼科装置は、蛍光反射の濃度値（蛍光強度）を測定し、その蛍光強度から涙液層の厚みを推測するものであり、涙液層の厚みを定量的に解析するものではなかった。また、角膜に蛍光染色液を点眼しなければならないため侵襲的検査となっている。
特許文献２に記載の眼科装置では、被検眼の涙液最表層（脂質層）の所定点に光が投影され、脂質層の表面と裏面で反射された反射光の干渉による干渉模様がＣＣＤカメラによって撮影される。ＣＣＤカメラで撮影された画像はモニタに表示され、検者はモニタに表示される画像から脂質層の状態や様子を観察して簡易的に診断を行う。この眼科装置では、ＣＣＤカメラによって観察された画像と指標となる参照画像とを見比べることで、被検眼の涙液層の状態を検者が主観的に決定するものである。したがって、非侵襲的検査ではあるが、涙液層を定量的に解析するものではない。
特許文献３に記載の眼科装置では、角膜下部の涙液メニスカスに対して所定のパターンの投影を行い、その反射パターンをＣＣＤカメラで撮影する。撮影された反射パターンは、涙液メニスカスによって変位しているため、その変位量から涙液メニスカスでの涙液量を評価する。この眼科装置は、涙液メニスカスの涙液量から角膜表面の涙液層の状態を間接的に評価するものであり、角膜表面上の涙液層を定量的に解析するものではなかった。
特開平６−２７７１７９号 特開平７−１３６１２０号 特開平１１−２６７１０２号

本発明は上述した実情に鑑みてなされたものであり、その目的は、角膜表面に形成される涙液層の時間の経過に伴う変化を非侵襲的に計測し、かつ、定量的に解析を行うことを可能にする眼科装置を提供することである。

本発明の眼科装置は、被検眼の角膜表面に形成される涙液層の時間の経過に伴う変化を測定する。この眼科装置は、角膜表面に所定パターンの光を投影する手段と、その投影された光の角膜表面からの反射像を撮影する手段と、撮影手段を所定の時間間隔で繰り返し実行させる手段と、撮影手段で撮影された複数の画像の中の１つを基準画像として、その基準画像に対する他の画像の相関度を算出する画像解析手段と、を有している。そして、画像解析手段は、（１）基準画像内に所定パターンと交差するように設定された解析線上の各点の座標値とその点における濃度値の関係を表す「座標値−濃度値」の関数を周波数解析して得られる「周波数−強度」の関数と、（２）比較対象となる他の画像内に設定された対応する解析線上の「座標値−濃度値」の関数を周波数解析して得られる「周波数−強度」の関数との相関度を算出することを特徴とする。
また、本明細書で開示する眼科装置の一態様は、被検眼の角膜表面に形成される涙液層の時間の経過に伴う変化を測定する眼科装置であり、角膜表面に所定パターンの光を投影する手段と、その投影された光の角膜表面からの反射像を撮影する手段と、撮影手段を所定の時間間隔で繰り返し実行させる手段と、撮影手段で撮影された複数の画像の中の１つを基準画像として、その基準画像に対する他の画像の相関度を算出する画像解析手段と、を有する。
この眼科装置では、被検眼の角膜表面に所定パターンの光を投影し、角膜表面からの反射像を所定の時間間隔で複数回撮影する。撮影によって得られた複数の画像は、これらの画像の中の１つを基準画像とし、他の画像を基準画像との相関度で評価する。すなわち、角膜表面から反射される反射像のパターンは、涙液層の状態が同一であれば同一パターンとして撮影され、涙液層の状態が変化すれば異なるパターンとして撮影される。したがって、基準画像との相関度を算出することで、その変化した程度を定量的に評価することが可能となる。この眼科装置では、撮影された画像から角膜表面における涙液層の厚み等を直接算出するのではなく、基準画像との相関度を算出することで角膜表面における涙液層の時間の経過に伴う変化を定量的に評価する。

上記の一態様に係る眼科装置の画像解析手段としては、例えば、（１）基準画像内に所定パターンと交差するように設定された解析線上の各点の座標値とその点における濃度値の関係を表す「座標値−濃度値」の関数と、（２）比較対象となる他の画像内に設定された対応する解析線上の「座標値−濃度値」の関数との相関度を算出することが好ましい。
このような構成では、撮影された画像内に所定パターンと交差するように解析線が設定され、その解析線上の各点の濃度値が求められる。撮影された画像内においては、光が反射されている部分と光が反射されていない部分の濃度値は大きく異なるため、解析線上の各点の濃度値は振動する波形となる。涙液層の状態が変化すると反射光が撮影される位置も変化することから、この解析線上の濃度値波形が変化する。したがって、この濃度値波形（すなわち、「座標値−濃度値」の関数）の相関度を算出することで、涙液層の変化を定量的に評価することができる。

上記の構成を採用する場合、画像解析手段は、（ａ）基準画像から得られた「座標値−濃度値」の関数を周波数解析して得られる「周波数−強度」の関数と、（ｂ）比較画像から得られた「座標値−濃度値」の関数を周波数解析して得られる「周波数−強度」の関数との相関度を算出することが好ましい。
このような構成によると、濃度値波形（「座標値−濃度値」の関数）を周波数解析することで、その振動特性がより明確化され、両者の相関度を適切に評価することができる。また、周波数解析を行うことでデータ数も少なくなり、両者の相関度を少ない計算量で算出することができる。

上記の一態様に係る眼科装置においては、所定パターンが角膜中心を中心として同心円状に略等間隔で配置された複数のリングパターンで構成されると共に、解析線が角膜中心から外周方向に直線状に設定されており、画像解析手段は、「周波数−強度」の関数の所定周波数領域について相関度を算出することが好ましい。
このような構成によると、濃度値波形（「座標値−濃度値」の関数）が所定の周期で振動する波形となるため、特定の周波数領域にピークを有することとなる。このため、周波数解析によって得られた「周波数−強度」の関数に対して所定周波数領域についての相関度を算出することで、２つの画像の相関度を評価することができる。また、所定周波数領域についてのみ相関度を算出することで、その他の領域に含まれるノイズの影響等を除去することができる。
所定周波数領域についての相関度を算出する場合は、例えば、画像解析手段は、「周波数−強度」の関数に所定周波数領域に帯域を有する帯域フィルタを乗算して得られる波形値の相関度を算出することができる。帯域フィルタを乗算することで、所定周波数領域の波形のみを取り出すことができる。

なお、上記の一態様に係る眼科装置においては、撮影手段で撮影された画像内に複数の解析線が設定されており、画像解析手段は解析線毎に相関度を算出することが好ましい。画像内に複数の解析線を設定することで、角膜表面全体の涙液層の状態を評価することができる。
このような構成を採用する場合、例えば、解析線が角膜中心から放射状に周方向に連続して複数設定されていることが好ましい。解析線を角膜中心から放射状に周方向に連続して複数設定することで、角膜全面をもれなく評価することができる。

また、本願は、被検眼の角膜表面に形成される涙液層の時間の経過に伴う変化を測定する眼科装置を制御するコンピュータに用いられる新規なプログラムを提供する。すなわち、本願に係るプログラムは、コンピュータを、照明光学系から角膜表面に所定パターンの光を投影させる照明制御手段と、その投影された光の角膜表面からの反射像を撮影装置に撮影させる撮影制御手段と、撮影制御手段を所定の時間間隔で繰り返し実行させる繰り返し実行手段と、撮影手段で撮影された複数の画像の中の１つを基準画像として、その基準画像に対する他の画像の相関度を算出する画像解析手段として機能させるためのプログラムであり、前記画像解析手段は、（１）基準画像内に所定パターンと交差するように設定された解析線上の各点の座標値とその点における濃度値の関係を表す「座標値−濃度値」の関数を周波数解析して得られる「周波数−強度」の関数と、（２）比較対象となる他の画像内に設定された対応する解析線上の「座標値−濃度値」の関数を周波数解析して得られる「周波数−強度」の関数との相関度を算出することを特徴とする。
このプログラムを従来からある角膜形状解析装置等にインストールすることで、涙液層の時間経過に伴う変化を定量的に評価することができる。

以下、本発明を具現化した一実施形態に係る眼科装置（角膜形状解析装置）について説明する。図１は本実施形態に係る眼科装置の光学系の構成を示す模式図であり、図２は図１に示す光学系に配置されるコーンの斜視図であり、図３は本実施形態の眼科装置において角膜に投影される投影パターンを示す図である。
図１に示すように、本実施形態の眼科装置は、被検眼Ｃの角膜に対向して位置決めされるコーン（０１）と、コーン（０１）の背面に配置された照明装置（０２）（例えば、ＬＥＤ等）を備えている。図２に示すようにコーン（０１）は、内部が中空の円錐形状の筒体であり、透明樹脂によって形成されている。コーン（０１）の内壁面には同心円状のパターンが印刷された透明フィルムが貼り付けられる一方で、コーン（０１）の外壁面には光を反射する塗装が施されている。したがって、コーン（０１）の背面に配置された照明装置（０２）による照明光は、コーン（０１）内で散乱し、内壁面に貼り付けられた透明フィルムによって光の一部が遮られ、被検眼Ｃの角膜に投影される。これによって、被検眼の角膜には、図３に示すような同心円状のパターン光が投影される。なお、コーン（０１）によって所望のパターンを投影する方法としては、コーン（０１）の内側に所望のパターンを直接切削し、その部分に遮光性の塗装を施す方法を採ることもできる。
照明装置（０２）の背面にはレンズ（０３）及びハーフミラー（０４）が配置される。ハーフミラー（０４）には、固視灯点光源（０６）の光がレンズ（０５）を介して入射するようになっている。固視灯点光源（０６）の光は、ハーフミラー（０４）で反射され、レンズ（０３）を透過してコーン（０１）に入射する。コーン（０１）に入射した光は、コーン（０１）の中央開口部から被検眼Ｃの角膜に照射される。図３に示すように、固視灯点光源（０６）の光はコーン（０１）によって形成される同心円状パターンの中心に位置するように調整されており、角膜観察時の固視灯となる。また、被検眼Ｃの角膜表面で反射された光（すなわち、角膜表面に投影された同心円状パターンの光と固視灯の光の反射光）は、は、コーン（０１）、レンズ（０３）及びハーフミラー（０４）を透過し、レンズ（０７）によって焦点調整が行われた後にＣＣＤカメラ（０８）で観察される。
なお、上述した眼科装置で角膜形状を正確に計測するためには、被検眼Ｃの角膜頂点を中心とする同心円状にリングパターンを投影する必要がある。このため、角膜形状計測時には固視灯点光源（０６）を点灯し、その光を被検者に固視させる。検者はＣＣＤカメラ（０８）で観察される固視灯点光源（０６）の光が画像中央となるように、被検眼Ｃに対して光学系の位置を調整する。これによって、被検眼Ｃの角膜頂点を中心とする同心円状にリングパターンが投影されることとなる。

上述した眼科装置の制御系の構成について説明する。図４は本実施形態に係る眼科装置の制御系の構成を示している。図４に示すように、眼科装置の制御はコントローラ（１２）によって行われている。
コントローラ（１２）は、ＣＰＵ，ＲＯＭ，ＲＡＭ等を備えたコンピュータによって構成することができる。コントローラ１２には、ＣＣＤカメラ（０８）、撮影開始スイッチ（１１）及びメモリ（１０）が接続されている。撮影開始スイッチ（１１）は検者によって操作される。検者が撮影開始スイッチ（１１）を操作すると、コントローラ（１２）は所定のプログラムを起動して、角膜表面に形成される涙液層の時間経過に伴う変化量を測定する。ＣＣＤカメラ（０８）によって撮影された画像データはコントローラ（１２）に入力され、メモリ（１０）等に記憶される。
また、コントローラ１２には、照明装置（０２）、固視灯点光源（０６）及び表示装置（０９）が接続されている。コントローラ（１２）は、照明装置（０２）及び固視灯点光源（０６）のＯＮ／ＯＦＦを行い、また、ＣＣＤカメラ（０８）で撮影された画像等を表示装置（０９）に表示する。

次に、上述した眼科装置を用いて角膜表面に形成される涙液層の時間的変化を測定する手順について説明する。図５は測定手順を示すフローチャートである。
図５に示すように、まず、検者は固視灯点光源（０６）を点灯し、被検者に固視灯を固視するよう指示する（Ｓ０１）。固視灯点光源（０６）が点灯されると、ＣＣＤカメラ（０８）で観察される画像が表示装置（０９）に表示される。表示装置（０９）には、固視灯点光源（０６）による輝点と、被検眼Ｃの角膜の画像が表示される。なお、固視灯点光源（０６）と照明装置（０２）の点灯は同時に行うようにしてもよい。
次に、検者は表示装置（０９）に表示される被検眼の前眼部が画面上の中央に位置するように（すなわち、固視灯点光源（０６）による反射光が角膜Ｃの中央に位置決めされるように）、被検眼Ｃに対して光学系の位置調整を行う（Ｓ０２）。
光学系の位置調整が終了すると照明装置（０２）を点灯し、被検眼Ｃの角膜表面に同心円状のパターン光を投影する（Ｓ０３）。次に、検者は、被検者にいったん閉瞼させて均質な涙液層を角膜表面に形成させ、被検者が開瞼すると同時に撮影開始スイッチ（１１）を操作する（Ｓ０４）。
検者によって撮影開始スイッチ（１１）が操作されると、コントローラ（１２）は、ＣＣＤカメラ（０８）を操作して、角膜表面からの反射パターンを所定の時間間隔で所定コマ数だけ撮影する（Ｓ０５）。本実施形態では、撮影開始スイッチ（１１）が操作されてから１０秒が経過するまで、１秒毎に角膜表面からの反射パターンを撮影する。したがって、ステップＳ０５の処理によって１１枚の角膜表面画像が撮影される。撮影された一連の画像は、メモリ（１０）に時系列順に格納される。図６にはステップＳ０５の処理によって取得された一連の画像が時系列順に示されている。なお、ステップＳ０５の撮影処理が完了するまでは、検者に瞬目をしないよう指示しておく。

ステップＳ０５の撮影処理が完了すると、次に、コントローラ（１２）はステップＳ０５で撮影された各画像を解析する画像解析処理を開始する（Ｓ０６）。図７は画像解析処理の手順を示すフローチャートである。図７に示すようにコントローラ（１２）は、まず、撮影開始スイッチ（１１）の操作と同時に撮影された画像（経過時間ｔ＝０秒）を基準画像としてメモリ（１０）より読み出す（Ｓ１０）。

次に、読み出した画像データに対して、その画像データ内に設定された解析線毎に、その解析線上の各点の濃度値をそれぞれ取得し、各解析線の濃度値波形を取得する（Ｓ１２）。図８に示すように本実施形態では、撮影された画像の中心（角膜中心）から放射状に２５６本の解析線が設定されている。これらの解析線は周方向に等間隔で設定されており、隣接する解析線がなす中心角は（３６０／２５６）°である。これらの解析線によって角膜表面全体の涙液層が評価できるようになっている。ステップＳ１２では、これら２５６本の解析線のそれぞれについて濃度値波形を取得する。具体的には、解析線毎に、解析線上の各点の座標値ｎ（角膜中心からの距離）を撮影された画像内の座標値（θ，Ｒ）に変換し、その座標値（θ，Ｒ）における濃度値（グレイレベル）を各点の濃度値とする。これによって、各解析線について、解析線上の座標値ｎと濃度値（グレイレベル）の関係（濃度値波形）を取得することができる。
図９にはＣＣＤカメラ（０８）で撮影された画像と、この画像内に設定された解析線が例示されている（ただし、図９では水平方向に引かれた１本の解析線のみが図示されている。）。図９から明らかなように、解析線は同心円状の反射パターンと直交し、解析線上には反射パターンと交差する部分（明るい部分）と反射パターンと交差しない部分（暗い部分）が交互に現れている。また、解析線は、角膜中心から同心円状の反射パターンの８番目のリングと交差する位置まで設定されている。８番目のリングと交差する位置までとしたのは、外周部が瞼にかかる度合いが多くなるためであり、また、８本目と９本目とのリングの間隔が若干広いためである。
図１０には、図９に示す解析線の濃度値波形（位置と濃度値の関係）を示している。図１０から明らかなように、解析線の濃度値波形は振動し、反射パターンと交差する位置で濃度値が高くなり、反射パターンと交差しない位置で濃度値が低くなっている。

ステップＳ１２で各解析線の濃度値波形を取得すると、コントローラ（１２）は、それら取得された濃度値波形をそれぞれフーリエ変換する（Ｓ１４）。フーリエ変換することで、濃度値波形の振動特性が抽出されることとなる。図１１にはフーリエ変換された後の波形が示されている。図中、周波数軸は横軸に、各周波数における強度を縦軸にとっている。図から明らかなように、フーリエ変換で得られた波形は、特定の周波数（周波数−２０と＋２０）にピークを有している。これは、反射パターンのリング間隔が略等間隔であり、このリング間隔に応じた周波数帯で周波数強度が大きくなるためである。
ステップＳ１６に進むと、コントローラ（１２）は、解析線毎に、ステップＳ１４で得られたフーリエ変換波形に帯域フィルタを乗算する（Ｓ１６）。ステップＳ１６で乗算する帯域フィルタは、図１１に示す三角形状の帯域フィルタであり、フーリエ変換波形のピークとなる周波数（周波数−２０と＋２０）を頂点としている。このような帯域フィルタをフーリエ変換波形に乗じることで、フーリエ変換波形の振動部分（反射パターンと交差する部分）の周波数特性のみを抽出することができる。
ステップＳ１８では、解析線毎に、帯域フィルタを乗じることで得られた波形値の絶対値を加算し、各解析線の評価値を算出する（Ｓ１８）。この評価値が、各解析線の相似度を評価するための値として用いられる。なお、波形値の二乗和を算出し、この値を評価地としてもよい。

ステップＳ２０に進むと、コントローラ（１２）は、基準画像と比較する比較画像（経過時間ｔ＝１，・・９）の一つをメモリ（１０）より読み出す。比較画像を読み出すと、その読み出した比較画像に対してステップＳ１２〜Ｓ１８の処理を行う（Ｓ２２）。これによって、比較画像の各解析線について評価値が算出される。
ステップＳ２４では、コントローラ（１２）は比較画像の各解析線について、基準画像の対応する解析線に対する相関度（インデックス値）を算出する。具体的には、比較画像の評価値を基準画像の対応する評価値で除算する。なお、ステップＳ２４で算出されたインデックス値が１．０を超える場合は、その解析線のインデックス値は１．０とする。
ステップＳ２６に進むと、コントローラ（１２）は全ての画像の各解析線についてインデックス値を算出したか否かを判断する。全ての画像の各解析線についてインデックス値を算出している場合（ステップＳ２６でＹＥＳ）はそのまま処理を終了し、全ての画像の各解析線についてインデックス値を算出していない場合（ステップＳ２６でＮＯ）はステップＳ２０に戻ってステップＳ２０からの処理を繰返す。これによって、全ての比較画像の各解析線についてインデックス値が算出される。

画像解析処理が終了すると、図５に戻ってコントローラ（１１）は表示装置（０９）に画像解析処理の結果を表示する（Ｓ０７）。本実施形態では、画像解析処理によって得られた各解析線のインデックス値を所定の角度範囲で平均を算出し、その平均値を画面に表示する。図１２は表示装置（１２）に表示される画像の一例を示している。図１２に示す例では、角膜を表す円を半径方向に２等分すると共に、周方向に８等分している。内側の円内には５秒後の画像のインデックス値が表示され、外側の円内には９秒後の画像のインデックス値が表示されている。また、８等分した角度範囲内では、その角度範囲に所属する解析線についてインデックス値を平均し、その平均値を表示している。インデックス値の表示は、図１３に示すように、その値に応じて適当な色彩を配色して表示することができる。

図１４には涙液層が安定している被検眼に対して測定した結果が示されている。図１４Ｂは開瞼から１０秒後の画像を示している。図１４Ａから明らかなように、角膜表面で反射された反射パターンのリング形状が明確に認識でき、涙液層が非常に安定していることが分かる。また、図１４Ａでは、内側の円内に５秒後のインデックス値が示され、外側の円内に９秒後のインデックス値が示されている。図１４Ｂから明らかなように、５秒後の状態と９秒後の状態は略同一であり、涙液層が安定していることを一目で把握することができる。
図１５には涙液層が安定していない被検眼に対して測定した結果が示されている。図１５Ｂに示す１０秒後の画像では反射パターンのリング形状が滲んでおり、特に、角膜の下半分においてその乱れが顕著である。また、図１５Ａから明らかなように、内側の円内の５秒後の状態に比べて、外側の円内の９秒後の状態ではインデックス値が低くなっている。なお、角膜の上部はまぶたがかかるために、アーチファクトが観察される場合も考えられる。

上述した説明から明らかなように、本実施形態の眼科装置では、ＣＣＤカメラ（０８）で撮影された画像データ全体を情報源として、基準画像（経過時間ｔ＝０秒）と比較画像（経過時間ｔ＝１〜９）との相関度を算出する。すなわち、角膜に投影されたリングパターンの反射像のうち、角膜表面の涙液層の乱れによって滲んだり隣同士で癒着したりしたような部位の情報を活用して、基準画像と比較画像の相関度を算出する。このため、基準画像と比較画像の相関度が角膜表面の涙液層の状態を表すインデックスとなり、そのインデックス値によって涙液層の時間の経過に伴う変化を定量的に評価することができる。
すなわち、従来の角膜形状解析装置では、角膜に投影されたリングパターンのうち、角膜表面の涙液層の乱れによって、そのリングパターンの反射像が滲んだり隣同士で癒着したりしたような部位は、通常、計測欠点として扱われ、周囲の数値からその部位の数値が補間されていた（すなわち、計測欠点の情報は解析に用いられていなかった）。これに対して、本実施形態では、従来は用いられなかった計測欠点を生じる部位の情報をあえて活用し、これによって、角膜をおおう涙液層の時間的変化を定量的に計測している。

以上、本発明の好適ないくつかの実施形態について詳細に説明したが、これらは例示に過ぎず、本発明は当業者の知識に基づいて種々の変更、改良を施した形態で実施することができる。
例えば、角膜表面に投影される投影パターンは同心円状のものに限られず、その形状等は任意に設定することができる。なお、投影パターンは、角膜全面に、均一で、かつ、偏在することなく光を投影できるものであることが好ましい。角膜全面を均一に評価できるためである。
また、画像内に設定される解析線は、角膜中心から放射状に設定するものに限られず、種々に設定することができる。その際は、角膜表面に投影される投影パターンと解析線が交差するように設定することが好ましい。このような構成とすることで、投影パターンの反射像の滲みや癒着等が解析線の濃度値波形に反映され、涙液層の時間的変化を定量的に算出することができる。
また、角膜表面からの反射パターンを撮影する回数や、その時間間隔は種々に設定することができる。例えば、最大開瞼時間内に５〜１０程度の画像を撮影することができるように決定することができる。また、撮影間隔も等間隔とする必要は必ずしもない。
また、上述した実施形態では、被検眼に対する光学系のアライメントをマニュアルで行うようにしたが、光学系のアライメントは自動で行うようにしてもよい。さらに、本実施形態では、濃度値波形をフーリエ変換したが、この他にも種々の周波数解析の方法を用いることができる。

なお、本明細書または図面に説明した技術要素は、単独であるいは各種の組み合わせによって技術的有用性を発揮するものであり、出願時請求項記載の組み合わせに限定されるものではない。また、本明細書または図面に例示した技術は複数目的を同時に達成するものであり、そのうちの一つの目的を達成すること自体で技術的有用性を持つものである。

本実施形態に係る眼科装置の光学系の全体構成を示す模式図。 図１に示す光学系に配置されるコーンの斜視図。 コーンによって投影される投影パターンを示す図。 図１に示す眼科装置の制御系の構成を示すブロック図。 角膜表面に形成される涙液層の時間的変化を測定するためのフローチャート。 図５の測定によって撮影される連続画像を時系列順に表示した図。 画像解析処理の手順を示すフローチャート。 撮影された画像内に設定された解析線を模式的に示す図。 実際に撮影された画像とその画像内に設定された解析線を併せて示す図。 図９に示す解析線の濃度値波形を示す図。 図９に示す濃度値波形のフーリエ変換波形と、そのフーリエ変換波形に乗算する帯域フィルタとを併せて示す図。 表示装置に表示される機能画像を説明するための図。 インデックス値と、そのインデックス値に対応する配色との関係を示す表。 涙液層が安定している症例のインデックス値の状態を示す機能画像。 涙液層が安定している症例を示す図（開瞼後９秒経過したときに撮影された画像）。 涙液層が不安定な症例のインデックス値の状態を示す機能画像。 涙液層が不安定な症例を示す図（開瞼後９秒経過したときに撮影された画像）。

符号の説明

０１：コーン
０２：照明装置
０３：レンズ
０４：ハーフミラー
０６：固視灯点光源
０８：ＣＣＤカメラ
１０：メモリ
１２：コントローラ

Claims (6)

1. 被検眼の角膜表面に形成される涙液層の時間の経過に伴う変化を測定する眼科装置であって、
   角膜表面に所定パターンの光を投影する手段と、
   その投影された光の角膜表面からの反射像を撮影する手段と、
   撮影手段を所定の時間間隔で繰り返し実行させる手段と、
   撮影手段で撮影された複数の画像の中の１つを基準画像として、その基準画像に対する他の画像の相関度を算出する画像解析手段と、を有しており、
   画像解析手段は、（１）基準画像内に所定パターンと交差するように設定された解析線上の各点の座標値とその点における濃度値の関係を表す「座標値−濃度値」の関数を周波数解析して得られる「周波数−強度」の関数と、（２）比較対象となる他の画像内に設定された対応する解析線上の「座標値−濃度値」の関数を周波数解析して得られる「周波数−強度」の関数との相関度を算出することを特徴とする眼科装置。
2. 所定パターンが角膜中心を中心として同心円状に略等間隔で配置された複数のリングパターンによって構成されると共に、解析線が角膜中心から外周方向に直線状に設定されており、画像解析手段は、「周波数−強度」の関数の所定周波数領域について相関度を算出することを特徴とする請求項１に記載の眼科装置。
3. 画像解析手段は、「周波数−強度」の関数に所定周波数領域に帯域を有する帯域フィルタを乗算して得られる波形値の相関度を算出することを特徴とする請求項２に記載の眼科装置。
4. 撮影手段で撮影された画像内には複数の解析線が設定されており、画像解析手段は解析線毎に相関度を算出することを特徴とする請求項１〜３のいずれかに記載の眼科装置。
5. 解析線が角膜中心から放射状に周方向に連続して複数設定されていることを特徴とする請求項４に記載の眼科装置。
6. コンピュータを、
   照明光学系から角膜表面に所定パターンの光を投影させる照明制御手段と、
   その投影された光の角膜表面からの反射像を撮影装置に撮影させる撮影制御手段と、
   撮影制御手段を所定の時間間隔で繰り返し実行させる繰り返し実行手段と、
   撮影手段で撮影された複数の画像の中の１つを基準画像として、その基準画像に対する他の画像の相関度を算出する画像解析手段として機能させるためのプログラムであって、
   前記画像解析手段は、（１）基準画像内に所定パターンと交差するように設定された解析線上の各点の座標値とその点における濃度値の関係を表す「座標値−濃度値」の関数を周波数解析して得られる「周波数−強度」の関数と、（２）比較対象となる他の画像内に設定された対応する解析線上の「座標値−濃度値」の関数を周波数解析して得られる「周波数−強度」の関数との相関度を算出することを特徴とするプログラム。