JP2964132B2 - 眼球運動解析装置 - Google Patents
眼球運動解析装置Info
- Publication number
- JP2964132B2 JP2964132B2 JP8309095A JP30909596A JP2964132B2 JP 2964132 B2 JP2964132 B2 JP 2964132B2 JP 8309095 A JP8309095 A JP 8309095A JP 30909596 A JP30909596 A JP 30909596A JP 2964132 B2 JP2964132 B2 JP 2964132B2
- Authority
- JP
- Japan
- Prior art keywords
- vector
- area
- movement
- comparison
- eyeball
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Eye Examination Apparatus (AREA)
Description
【0001】
【発明の属する技術分野】本発明は、眼球運動解析装置
に関し、より詳細には、眼振を含む眼球の運動を正確に
定量的に測定する眼球運動解析装置に関するものであ
る。
に関し、より詳細には、眼振を含む眼球の運動を正確に
定量的に測定する眼球運動解析装置に関するものであ
る。
【0002】
【従来の技術】眼振など眼球の運動を検査するために、
たとえば、Frezel眼鏡を用いて外部から眼球の運動を観
察する手法、或いは、電気眼振検査(ENG)が広く用
いられている。前者は、簡便で手軽に行うことができる
一方、固視の影響を除去することが困難であり、また、
眼振の有無やその性状の判定が、検者の主観的な判断に
依存すしているという問題点がある。後者は、固視の影
響をほぼ完全に除去でき、眼振の定量的測定が可能であ
るが、その原理上、上下方向、左右方向の眼球運動の解
析に限られるという問題点がある。
たとえば、Frezel眼鏡を用いて外部から眼球の運動を観
察する手法、或いは、電気眼振検査(ENG)が広く用
いられている。前者は、簡便で手軽に行うことができる
一方、固視の影響を除去することが困難であり、また、
眼振の有無やその性状の判定が、検者の主観的な判断に
依存すしているという問題点がある。後者は、固視の影
響をほぼ完全に除去でき、眼振の定量的測定が可能であ
るが、その原理上、上下方向、左右方向の眼球運動の解
析に限られるという問題点がある。
【0003】これに対して、赤外線CCDカメラを用い
て眼球を撮影して、これをビデオなどに記録し、記録さ
れたビデオ画像に基づき、眼球運動を解析する装置が提
案されている。
て眼球を撮影して、これをビデオなどに記録し、記録さ
れたビデオ画像に基づき、眼球運動を解析する装置が提
案されている。
【0004】たとえば、特開平4−15032号公報に
は、眼球の撮影画像を記録し、記録された画像データか
ら、追跡対象となる眼球中心位置と眼球中心外位置とを
特定し、これら位置の座標を求め、前述した座標の時間
変化から回旋成分を含む眼球運動を算出するように構成
された眼球運動測定装置が開示されている。このよう
に、所定の時間ごとの、眼球の画像に基づき、画像中の
所定の位置の時間変化を測定して、眼球の移動を決定す
ることにより、比較的短い測定時間にて眼球の移動を測
定することが可能となる。
は、眼球の撮影画像を記録し、記録された画像データか
ら、追跡対象となる眼球中心位置と眼球中心外位置とを
特定し、これら位置の座標を求め、前述した座標の時間
変化から回旋成分を含む眼球運動を算出するように構成
された眼球運動測定装置が開示されている。このよう
に、所定の時間ごとの、眼球の画像に基づき、画像中の
所定の位置の時間変化を測定して、眼球の移動を決定す
ることにより、比較的短い測定時間にて眼球の移動を測
定することが可能となる。
【0005】
【発明が解決しようとする課題】しかしながら、上述し
た眼球運動測定装置においては、追跡対象として、虹彩
紋理のような特定の場所を眼球中心外位置に対応させて
いる、画像中の虹彩紋理の位置に突発的なノイズが生じ
ていた場合などに、眼球の移動を適切に測定することが
できないという問題点があった。また、光の反射により
虹彩紋理の撮影が妨害されることを防止するため、カメ
ラに偏光フィルタを設けるなど、複雑な光学系を要する
という問題点があった。
た眼球運動測定装置においては、追跡対象として、虹彩
紋理のような特定の場所を眼球中心外位置に対応させて
いる、画像中の虹彩紋理の位置に突発的なノイズが生じ
ていた場合などに、眼球の移動を適切に測定することが
できないという問題点があった。また、光の反射により
虹彩紋理の撮影が妨害されることを防止するため、カメ
ラに偏光フィルタを設けるなど、複雑な光学系を要する
という問題点があった。
【0006】本発明の目的は、画像に突発的なノイズが
生じていた場合であっても、正確に、眼球運動を測定可
能であり、かつ、その構成も比較的単純な眼球運動解析
装置を提供することを目的とする。
生じていた場合であっても、正確に、眼球運動を測定可
能であり、かつ、その構成も比較的単純な眼球運動解析
装置を提供することを目的とする。
【0007】
【課題を解決するための手段】本発明の目的は、カメラ
にて眼球を撮影したフレーム画像に基づき、眼球の運動
を測定する眼球運動解析装置であって、前記フレーム画
像から複数の部分領域を得て、各部分領域中の比較対象
とすべき比較領域を決定する比較領域決定手段と、前記
フレーム画像の比較領域と、時間的に隣接する隣接フレ
ーム画像の対応する部分領域中の、前記比較領域と同じ
大きさの領域との間でパターンマッチングを実行して、
前記比較領域の、隣接フレーム画像中の位置を決定する
マッチング実行手段と、マッチング実行手段により得ら
れた、前記比較領域の、隣接フレーム画像中の位置に基
づき、前記比較領域の移動方向および移動距離からなる
ベクトルを生成するベクトル生成手段と、ベクトル生成
手段により得られた、比較領域の各々に関するベクトル
を統計的に処理して、隣接するフレーム画像の間の眼球
の移動を示す眼球移動ベクトルを生成する統計処理手段
と、統計処理手段により得られた眼球移動ベクトルに基
づき、眼球の移動を示すグラフを示すグラフ生成手段と
を備え、当該グラフが出力されるように構成され、か
つ、前記統計処理手段が、前記複数のベクトルの移動方
向および/または移動距離に関する分布を生成し、前記
分布が所定の分布にしたがっているか否かを判断し、分
布にしたがっていないベクトルがある場合には、当該ベ
クトルに関連する部分領域について、前記比較領域決定
手段が、再度比較領域を決定し、当該比較領域に関する
パターンマッチングおよびベクトルの生成が実行される
ように構成されたことを特徴とする眼球運動解析装置に
より達成される。
にて眼球を撮影したフレーム画像に基づき、眼球の運動
を測定する眼球運動解析装置であって、前記フレーム画
像から複数の部分領域を得て、各部分領域中の比較対象
とすべき比較領域を決定する比較領域決定手段と、前記
フレーム画像の比較領域と、時間的に隣接する隣接フレ
ーム画像の対応する部分領域中の、前記比較領域と同じ
大きさの領域との間でパターンマッチングを実行して、
前記比較領域の、隣接フレーム画像中の位置を決定する
マッチング実行手段と、マッチング実行手段により得ら
れた、前記比較領域の、隣接フレーム画像中の位置に基
づき、前記比較領域の移動方向および移動距離からなる
ベクトルを生成するベクトル生成手段と、ベクトル生成
手段により得られた、比較領域の各々に関するベクトル
を統計的に処理して、隣接するフレーム画像の間の眼球
の移動を示す眼球移動ベクトルを生成する統計処理手段
と、統計処理手段により得られた眼球移動ベクトルに基
づき、眼球の移動を示すグラフを示すグラフ生成手段と
を備え、当該グラフが出力されるように構成され、か
つ、前記統計処理手段が、前記複数のベクトルの移動方
向および/または移動距離に関する分布を生成し、前記
分布が所定の分布にしたがっているか否かを判断し、分
布にしたがっていないベクトルがある場合には、当該ベ
クトルに関連する部分領域について、前記比較領域決定
手段が、再度比較領域を決定し、当該比較領域に関する
パターンマッチングおよびベクトルの生成が実行される
ように構成されたことを特徴とする眼球運動解析装置に
より達成される。
【0008】本発明によれば、フレーム画像中の部分領
域中の比較領域の各々が、隣接するフレーム画像間で、
どれだけ移動しているかを示すベクトルを得て、比較領
域に関して得られたベクトルを統計的に処理して、眼球
の移動を示す眼球移動ベクトルを生成している。このた
め、画像中に突発的なノイズが生じている場合であって
も、そのようなノイズが生じた領域に関するベクトルが
フィルタリングにより排除され、適切に眼球の移動を測
定することが可能となる。また、眼球の画像中の特徴点
を切り出すなどの操作或いは処理を必要とすることな
く、眼球の移動を測定することが可能となる。
域中の比較領域の各々が、隣接するフレーム画像間で、
どれだけ移動しているかを示すベクトルを得て、比較領
域に関して得られたベクトルを統計的に処理して、眼球
の移動を示す眼球移動ベクトルを生成している。このた
め、画像中に突発的なノイズが生じている場合であって
も、そのようなノイズが生じた領域に関するベクトルが
フィルタリングにより排除され、適切に眼球の移動を測
定することが可能となる。また、眼球の画像中の特徴点
を切り出すなどの操作或いは処理を必要とすることな
く、眼球の移動を測定することが可能となる。
【0009】また、本発明によれば、ノイズなどが生じ
ており、所定の分布に従わないベクトルに関する部分領
域にて、再度比較領域を決定して、パターンマッチング
によりベクトルを生成するため、得られる眼球移動ベク
トルをより適切化することができる。本発明のさらに好
ましい実施態様においては、前記分布にしたがっていな
いベクトルがある場合に、前記比較領域決定手段が、関
連する部分領域をさらに小領域に分割して、各小領域に
含まれる第2の比較領域を生成するように構成されてい
る。
ており、所定の分布に従わないベクトルに関する部分領
域にて、再度比較領域を決定して、パターンマッチング
によりベクトルを生成するため、得られる眼球移動ベク
トルをより適切化することができる。本発明のさらに好
ましい実施態様においては、前記分布にしたがっていな
いベクトルがある場合に、前記比較領域決定手段が、関
連する部分領域をさらに小領域に分割して、各小領域に
含まれる第2の比較領域を生成するように構成されてい
る。
【0010】本発明のさらに好ましい実施態様において
は、さらに、前記部分領域を分割することにより得られ
た第2の比較領域のそれぞれの移動を示すベクトルを統
計的に処理し、前記部分領域に関連するもとの比較領域
の移動を示すベクトルを得る小領域ベクトル統計処理手
段を備えている。
は、さらに、前記部分領域を分割することにより得られ
た第2の比較領域のそれぞれの移動を示すベクトルを統
計的に処理し、前記部分領域に関連するもとの比較領域
の移動を示すベクトルを得る小領域ベクトル統計処理手
段を備えている。
【0011】或いは、前記分布にしたがっていないベク
トルがある場合に、前記比較領域決定手段が、前記部分
領域において、少なくとも相互に一画素以上離間して位
置する画素をサンプリングすることにより、比較領域を
得るように構成されていても良い。
トルがある場合に、前記比較領域決定手段が、前記部分
領域において、少なくとも相互に一画素以上離間して位
置する画素をサンプリングすることにより、比較領域を
得るように構成されていても良い。
【0012】本発明の好ましい実施態様においては、さ
らに、フレーム画像中の瞳孔の位置を決定する瞳孔位置
決定手段と、隣接するフレーム間の瞳孔の移動を決定
し、瞳孔の移動を示す瞳孔移動ベクトルを生成する瞳孔
移動測定手段と、前記ベクトル生成手段により生成され
た、各比較領域に関するベクトルから、瞳孔移動ベクト
ルを減算し、眼球の回旋を示す比較領域回旋ベクトルを
生成する回旋ベクトル生成手段とを備え、前記統計処理
手段が、回旋ベクトル生成手段により得られた、比較領
域の各々の回旋に関するベクトルを統計的に処理して、
隣接するフレーム画像間の眼球の回旋を示す眼球回旋ベ
クトルを生成するように構成されている。
らに、フレーム画像中の瞳孔の位置を決定する瞳孔位置
決定手段と、隣接するフレーム間の瞳孔の移動を決定
し、瞳孔の移動を示す瞳孔移動ベクトルを生成する瞳孔
移動測定手段と、前記ベクトル生成手段により生成され
た、各比較領域に関するベクトルから、瞳孔移動ベクト
ルを減算し、眼球の回旋を示す比較領域回旋ベクトルを
生成する回旋ベクトル生成手段とを備え、前記統計処理
手段が、回旋ベクトル生成手段により得られた、比較領
域の各々の回旋に関するベクトルを統計的に処理して、
隣接するフレーム画像間の眼球の回旋を示す眼球回旋ベ
クトルを生成するように構成されている。
【0013】この実施態様においては、各比較領域の、
平行移動の成分を取り除き回旋のみを示すベクトルを生
成して、これに基づき眼球の回旋を示す眼球回旋ベクト
ルを得るため、比較的簡単な処理により、適切に眼球の
回旋を測定することが可能となる。
平行移動の成分を取り除き回旋のみを示すベクトルを生
成して、これに基づき眼球の回旋を示す眼球回旋ベクト
ルを得るため、比較的簡単な処理により、適切に眼球の
回旋を測定することが可能となる。
【0014】或いは、画像解析装置は、上述した眼球の
回旋のみを測定するように構成されていても良い。
回旋のみを測定するように構成されていても良い。
【0015】また、本発明の目的は、上記画像解析装置
を作動させるためのプログラムが記憶された記憶媒体に
よっても実現される。
を作動させるためのプログラムが記憶された記憶媒体に
よっても実現される。
【0016】
【発明の実施の形態】以下、添付図面を参照して、本発
明の実施の形態につき詳細に説明を加える。図1は、本
発明の実施の形態にかかる眼球運動解析装置の構成を示
すブロックダイヤグラムである。図1に示すように、眼
球運動解析装置10は、赤外線CCD(Charge Coupled
Device)カメラ12と、赤外線CCDカメラ12からの
画像を記録するビデオテープレコーダ(VTR)14
と、VTR14から出力された画像信号をディジタル信
号に変換するアナログ/ディジタル(A/D)変換器1
6と、A/D変換器16から出力されたディジタルデー
タからフレームごとの画像データを一時的に記憶するフ
レームバッファ18と、画像データをフレームごとに記
憶する画像メモリ20と、フレームバッファ18に記憶
された所定のフレームの画像データを解析するフレーム
画像解析部22と、フレーム画像解析部22により得ら
れた解析結果に基づくグラフなどを作成する解析結果作
成部24と、プロッタ、プリンタなどの出力装置26
と、CRTなどの表示装置28と、これら構成部分を制
御する制御部30とを備えている。制御部30は、後述
するフレームバッファ18、フレーム画像解析部22、
解析結果作成部24などを制御して、これらに所定の動
作をさせるためのプログラムを記憶した記憶媒体(図示
せず)が設けられている。この記憶媒体は、ハードディ
スクであってもよいし、CD−ROM、フロッピーディ
スク、読み書き可能な光ディスクなどであっても良い。
明の実施の形態につき詳細に説明を加える。図1は、本
発明の実施の形態にかかる眼球運動解析装置の構成を示
すブロックダイヤグラムである。図1に示すように、眼
球運動解析装置10は、赤外線CCD(Charge Coupled
Device)カメラ12と、赤外線CCDカメラ12からの
画像を記録するビデオテープレコーダ(VTR)14
と、VTR14から出力された画像信号をディジタル信
号に変換するアナログ/ディジタル(A/D)変換器1
6と、A/D変換器16から出力されたディジタルデー
タからフレームごとの画像データを一時的に記憶するフ
レームバッファ18と、画像データをフレームごとに記
憶する画像メモリ20と、フレームバッファ18に記憶
された所定のフレームの画像データを解析するフレーム
画像解析部22と、フレーム画像解析部22により得ら
れた解析結果に基づくグラフなどを作成する解析結果作
成部24と、プロッタ、プリンタなどの出力装置26
と、CRTなどの表示装置28と、これら構成部分を制
御する制御部30とを備えている。制御部30は、後述
するフレームバッファ18、フレーム画像解析部22、
解析結果作成部24などを制御して、これらに所定の動
作をさせるためのプログラムを記憶した記憶媒体(図示
せず)が設けられている。この記憶媒体は、ハードディ
スクであってもよいし、CD−ROM、フロッピーディ
スク、読み書き可能な光ディスクなどであっても良い。
【0017】CCDカメラ12は、被験者の眼球周辺を
覆い、外部からの光が遮断されるようなゴーグル32に
取り付けられている。このゴーグルは、たとえば、「医
用電子と生体光学」第29巻特別号(1991)に開示されてい
るようなもので良い。このゴーグル内の所定の位置に
は、赤外線光源(図示せず)が設けられ、眼球に向けて
赤外光を照射するようになっている。
覆い、外部からの光が遮断されるようなゴーグル32に
取り付けられている。このゴーグルは、たとえば、「医
用電子と生体光学」第29巻特別号(1991)に開示されてい
るようなもので良い。このゴーグル内の所定の位置に
は、赤外線光源(図示せず)が設けられ、眼球に向けて
赤外光を照射するようになっている。
【0018】画像メモリ20は、たとえば、ハードディ
スクから構成され、解析に必要な複数のフレーム分の画
像データを記憶することができる。
スクから構成され、解析に必要な複数のフレーム分の画
像データを記憶することができる。
【0019】また、フレーム画像解析部22は、後述す
るアルゴリズムにしたがって、比較すべき二つのフレー
ムの画像データに基づき、画像データ中の所定の画素か
らなる領域が、どのくらい移動しているかを判断する。
るアルゴリズムにしたがって、比較すべき二つのフレー
ムの画像データに基づき、画像データ中の所定の画素か
らなる領域が、どのくらい移動しているかを判断する。
【0020】このように構成された眼球運動解析装置の
作動につき、以下に説明する。まず、操作者は、被験者
の眼球が覆われるようにゴーグルを装着し、次いで、赤
外線光源(図示せず)を被験者の眼球に照射する。赤外
線CCDカメラ12は、眼球の表面の反射像を撮像する
ため、操作者は、VTR14を作動させて、赤外線CC
Dカメラ12による眼球の画像を記録する。このように
して、所定の時間の撮影が終了すると、VTR14を再
生して、画像信号をA/Dコンバータ16に与える。
作動につき、以下に説明する。まず、操作者は、被験者
の眼球が覆われるようにゴーグルを装着し、次いで、赤
外線光源(図示せず)を被験者の眼球に照射する。赤外
線CCDカメラ12は、眼球の表面の反射像を撮像する
ため、操作者は、VTR14を作動させて、赤外線CC
Dカメラ12による眼球の画像を記録する。このように
して、所定の時間の撮影が終了すると、VTR14を再
生して、画像信号をA/Dコンバータ16に与える。
【0021】A/Dコンバータ16は、与えられた画像
信号をディジタル信号に変換し、これをフレームバッフ
ァ18に出力し、フレームバッファ18は、1フレーム
当たりのディジタル画像データが、その内部に形成され
ると、これを画像メモリ20に記憶する。このフレーム
ごとの画像の記憶に関する処理は、制御部30により制
御される。この処理を繰り返すことにより、画像メモリ
20に、解析に必要な、所定のフレーム数の画像データ
が記憶される。この実施の形態においては、1フレーム
あたり256×256ドット分の画素データからなるデ
ィジタルデータが得られる。また、毎秒30のフレーム
が形成される。また、上記256×256ドットのデー
タからなる画像には、まぶたなど、眼球運動を解析する
ために不要な部分が含まれるため、A/Dコンバータ1
6から出力されたデータに、マスキング処理を施し、所
定の120×120ドットの画素データからなる画像デ
ータが、画像メモリ20に記憶される。
信号をディジタル信号に変換し、これをフレームバッフ
ァ18に出力し、フレームバッファ18は、1フレーム
当たりのディジタル画像データが、その内部に形成され
ると、これを画像メモリ20に記憶する。このフレーム
ごとの画像の記憶に関する処理は、制御部30により制
御される。この処理を繰り返すことにより、画像メモリ
20に、解析に必要な、所定のフレーム数の画像データ
が記憶される。この実施の形態においては、1フレーム
あたり256×256ドット分の画素データからなるデ
ィジタルデータが得られる。また、毎秒30のフレーム
が形成される。また、上記256×256ドットのデー
タからなる画像には、まぶたなど、眼球運動を解析する
ために不要な部分が含まれるため、A/Dコンバータ1
6から出力されたデータに、マスキング処理を施し、所
定の120×120ドットの画素データからなる画像デ
ータが、画像メモリ20に記憶される。
【0022】画像メモリに必要な画像データが記憶され
ると、制御部30により、フレームバッファに、所定の
フレームの画像データが読み出され、この画像データに
基づき、フレーム画像解析部22による処理が実行され
る。
ると、制御部30により、フレームバッファに、所定の
フレームの画像データが読み出され、この画像データに
基づき、フレーム画像解析部22による処理が実行され
る。
【0023】ここに、フレーム画像解析部22は、以下
の仮定条件に基づき、処理を実行している。すなわち、
(1)時間の変化にかかわらず、眼球中の同じ場所に対応
する画素の画素値(濃淡値)は保存される。(2)眼球運
動は剛体運動とみなすことができる。(3)非直線運動
は、微小な直線運動の集合に近似される。
の仮定条件に基づき、処理を実行している。すなわち、
(1)時間の変化にかかわらず、眼球中の同じ場所に対応
する画素の画素値(濃淡値)は保存される。(2)眼球運
動は剛体運動とみなすことができる。(3)非直線運動
は、微小な直線運動の集合に近似される。
【0024】上述した仮定の下で、この実施の形態にか
かるフレーム画像解析部22は、以下の原理により、フ
レームごとの領域の移動を特定している。画像メモリ2
0に記憶されたフレームごとの画像データ(フレーム画
像)は、ある時間tにおける静止画であり、先に述べた
仮定にしたがえば、静止画を平行移動することにより、
隣接するフレームのフレーム画像を略一致させることが
可能となる。たとえば、図2に示すように、時間的に隣
接する第1のフレームのフレーム画像FR1と第2のフ
レームのフレーム画像FR2との間では、上述した仮定
を満たすならば、一方を平行移動することにより、他方
の画像と一致させることが可能となる。この図2におい
て、右斜め下の方向に、2画素だけ、フレーム画像FR
1を移動させることにより、フレーム画像FR2と一致
させることができる。
かるフレーム画像解析部22は、以下の原理により、フ
レームごとの領域の移動を特定している。画像メモリ2
0に記憶されたフレームごとの画像データ(フレーム画
像)は、ある時間tにおける静止画であり、先に述べた
仮定にしたがえば、静止画を平行移動することにより、
隣接するフレームのフレーム画像を略一致させることが
可能となる。たとえば、図2に示すように、時間的に隣
接する第1のフレームのフレーム画像FR1と第2のフ
レームのフレーム画像FR2との間では、上述した仮定
を満たすならば、一方を平行移動することにより、他方
の画像と一致させることが可能となる。この図2におい
て、右斜め下の方向に、2画素だけ、フレーム画像FR
1を移動させることにより、フレーム画像FR2と一致
させることができる。
【0025】次に、図3に示すように、時間的に隣接す
る二つのフレーム画像FR3、FR4を考える。これら
フレーム画像のうち、時間的に先行するフレーム画像F
R3のある特定の部分領域301(図3において、灰色
にて示した部分)をサンプリングして、この領域と、時
間的に後続するフレーム画像FR4との間で、パターン
マッチングを実行する。すなわち、フレーム画像FR4
において、部分領域301と等しいサイズの領域(たと
えば、311、312)を取り出し、領域を構成する画
素の画素値ごとの差の二乗和を得て、最も二乗和の小さ
い部分が、フレーム画像FR4において部分領域301
と対応する領域であると決定する。すなわち、最も二乗
和が小さくなるような位置に、部分領域301が移動し
ていると決定する。このようにして、時間的に隣接する
フレーム画像の移動を測定することができる。
る二つのフレーム画像FR3、FR4を考える。これら
フレーム画像のうち、時間的に先行するフレーム画像F
R3のある特定の部分領域301(図3において、灰色
にて示した部分)をサンプリングして、この領域と、時
間的に後続するフレーム画像FR4との間で、パターン
マッチングを実行する。すなわち、フレーム画像FR4
において、部分領域301と等しいサイズの領域(たと
えば、311、312)を取り出し、領域を構成する画
素の画素値ごとの差の二乗和を得て、最も二乗和の小さ
い部分が、フレーム画像FR4において部分領域301
と対応する領域であると決定する。すなわち、最も二乗
和が小さくなるような位置に、部分領域301が移動し
ていると決定する。このようにして、時間的に隣接する
フレーム画像の移動を測定することができる。
【0026】次に、より詳細に、フレーム画像解析部2
2の処理につき説明を加える。以下の処理を、場合によ
っては、ブロックマッチング法と称する。図4は、本実
施の形態にかかるフレーム画像解析部22の処理を示す
フローチャートである。図4に示すように、フレーム画
像解析部22においては、まず、フレーム画像の番号を
示すnを初期化する(ステップ401)。次いで、制御
部30に、時間的に隣接するn番目のフレーム画像FR
nと、(n+1)番目のフレーム画像FRn+1を、フレー
ムバッファ18に読み出すことを依頼する。制御部30
は、これに応答して、所定の、時間的に隣接するフレー
ム画像を画像メモリ20から読み出して、フレームバッ
ファ18上に展開する(ステップ402)。
2の処理につき説明を加える。以下の処理を、場合によ
っては、ブロックマッチング法と称する。図4は、本実
施の形態にかかるフレーム画像解析部22の処理を示す
フローチャートである。図4に示すように、フレーム画
像解析部22においては、まず、フレーム画像の番号を
示すnを初期化する(ステップ401)。次いで、制御
部30に、時間的に隣接するn番目のフレーム画像FR
nと、(n+1)番目のフレーム画像FRn+1を、フレー
ムバッファ18に読み出すことを依頼する。制御部30
は、これに応答して、所定の、時間的に隣接するフレー
ム画像を画像メモリ20から読み出して、フレームバッ
ファ18上に展開する(ステップ402)。
【0027】次いで、フレーム画像解析部22は、フレ
ーム画像FRn中の所定の部分領域を決定する(ステッ
プ403)。より詳細には、120×120ドットの画
素を、5×5=25分割し、24×24ドットの画素か
らなる領域を生成する。本発明において、解析すべき対
象は、眼球の緩徐相(ゆっくりとした解析運動)であ
り、フレーム速度を30フレーム/秒としたときには、
時間的に隣接するフレーム間の眼球の移動は、1ないし
2ドットとなる。このため、本実施の形態においては、
24×24ドットの画素の領域うち、上下左右に5ドッ
トずつ小さいサイズの14×14ドットの画素からなる
領域を、パターンマッチングをすべき比較領域とする
(図5(a)参照)。次いで、フレームFRn+1も同様
の手法で、24×24ドットの画素からなる複数の部分
領域を生成する(ステップ404)。
ーム画像FRn中の所定の部分領域を決定する(ステッ
プ403)。より詳細には、120×120ドットの画
素を、5×5=25分割し、24×24ドットの画素か
らなる領域を生成する。本発明において、解析すべき対
象は、眼球の緩徐相(ゆっくりとした解析運動)であ
り、フレーム速度を30フレーム/秒としたときには、
時間的に隣接するフレーム間の眼球の移動は、1ないし
2ドットとなる。このため、本実施の形態においては、
24×24ドットの画素の領域うち、上下左右に5ドッ
トずつ小さいサイズの14×14ドットの画素からなる
領域を、パターンマッチングをすべき比較領域とする
(図5(a)参照)。次いで、フレームFRn+1も同様
の手法で、24×24ドットの画素からなる複数の部分
領域を生成する(ステップ404)。
【0028】この後に、フレームFRnの比較領域と、
フレームFRn+1の部分領域中の同じサイズの領域とを
比較する(ステップ405)。より詳細には、フレーム
画像解析部22は、フレーム画像FRnの部分領域(図
5(a)の領域500)中の14×14ドットの比較領
域(図5(a)の領域501)と、フレーム画像FRn+
1の24×24ドットの対応する部分領域中の、比較領
域と同じサイズを有する14×14ドットの領域とのパ
ターンマッチングを、順次実行する。すなわち、図5
(b)に示すように、FRnの比較領域に含まれる画素
値の各々と、FRn+1(領域502)の対応する領域
(CR1、CR2、・・・、CRj、CRj+1、・・
・)に含まれる画素値の各々との差の二乗和を得て、最
も二乗和の小さい比較領域CRkの位置に、フレーム画
像FRn+1の比較領域が移動したと決定する。フレーム
画像解析部22は、上述したような比較領域の移動を、
25分割されたフレーム画像FRnの部分領域の全てに
関して実行する。
フレームFRn+1の部分領域中の同じサイズの領域とを
比較する(ステップ405)。より詳細には、フレーム
画像解析部22は、フレーム画像FRnの部分領域(図
5(a)の領域500)中の14×14ドットの比較領
域(図5(a)の領域501)と、フレーム画像FRn+
1の24×24ドットの対応する部分領域中の、比較領
域と同じサイズを有する14×14ドットの領域とのパ
ターンマッチングを、順次実行する。すなわち、図5
(b)に示すように、FRnの比較領域に含まれる画素
値の各々と、FRn+1(領域502)の対応する領域
(CR1、CR2、・・・、CRj、CRj+1、・・
・)に含まれる画素値の各々との差の二乗和を得て、最
も二乗和の小さい比較領域CRkの位置に、フレーム画
像FRn+1の比較領域が移動したと決定する。フレーム
画像解析部22は、上述したような比較領域の移動を、
25分割されたフレーム画像FRnの部分領域の全てに
関して実行する。
【0029】このようにして、分割されたフレーム画像
FRnの全てに関して、比較領域が時間的に隣接するフ
レームFRn+1の何れの位置に移動したかを決定する
と、得られた結果、すなわち比較領域の移動を示すベク
トル(移動方向および移動距離)が、前述した仮定条件
を満足しているか否かを判定する(ステップ406)。
これは、時間的に隣接するフレーム画像の間の、ある比
較領域の移動方向および移動距離を示すベクトルが、他
の比較領域を示すベクトルと逆方向を向いている場合
や、このベクトルの大きさが、他のベクトルの大きさと
著しく異なる場合に、当該ベクトルを不適切として排除
するためのフィルタリング処理に対応する。
FRnの全てに関して、比較領域が時間的に隣接するフ
レームFRn+1の何れの位置に移動したかを決定する
と、得られた結果、すなわち比較領域の移動を示すベク
トル(移動方向および移動距離)が、前述した仮定条件
を満足しているか否かを判定する(ステップ406)。
これは、時間的に隣接するフレーム画像の間の、ある比
較領域の移動方向および移動距離を示すベクトルが、他
の比較領域を示すベクトルと逆方向を向いている場合
や、このベクトルの大きさが、他のベクトルの大きさと
著しく異なる場合に、当該ベクトルを不適切として排除
するためのフィルタリング処理に対応する。
【0030】より具体的には、横軸にベクトルの方向或
いは大きさ、縦軸にその頻度を設定して、得られたベク
トルの値をプロットすると、ベクトルの個数が十分に大
きければガウス分布を得ることができることを利用して
いる。したがって、横軸上に二つの閾値を設定し、その
第1の閾値よりも小さい位置に表れたベクトル、或い
は、第2の閾値よりも大きい位置に表れたベクトルを見
つけ出す。すなわち、分布にしたがって起こり得ないで
あろう方向或いは大きさのベクトルを見つけ出す。
いは大きさ、縦軸にその頻度を設定して、得られたベク
トルの値をプロットすると、ベクトルの個数が十分に大
きければガウス分布を得ることができることを利用して
いる。したがって、横軸上に二つの閾値を設定し、その
第1の閾値よりも小さい位置に表れたベクトル、或い
は、第2の閾値よりも大きい位置に表れたベクトルを見
つけ出す。すなわち、分布にしたがって起こり得ないで
あろう方向或いは大きさのベクトルを見つけ出す。
【0031】このようなベクトルを得た部分領域が、仮
定条件を満足していない領域となる。
定条件を満足していない領域となる。
【0032】ステップ406において、ノー(No)と判
断された場合、すなわち、仮定条件を満足していないと
判断された場合には、ステップ403に戻り、フレーム
画像FRnに基づき、比較領域を決定する。この場合に
は、フレーム画像FRnの部分領域中、比較領域の移動
を示すベクトルが仮定条件を満たしていない部分領域の
みに関して、さらに、比較領域を決定する。より詳細に
は、図6に示すように、比較領域の移動を示す結果が仮
定条件を満たしていないような部分領域を、さらに分割
して、分割された小領域601ないし604を生成す
る。図6において、中央の破線で示す領域610が、前
回の処理にて得られた比較領域である。さらに、小領域
601のうち、上下左右に所定のドットずつ小さいサイ
ズの領域を比較領域611ないし614と決定する。こ
の実施の形態においては、24×24ドットのもとの部
分領域を2×2=4分割して、12×12ドットの画素
からなる部分領域を作成し、各部分領域中の、上下左右
に3ドットずつ小さいサイズの6×6ドットの画素から
なる領域を比較領域としている。
断された場合、すなわち、仮定条件を満足していないと
判断された場合には、ステップ403に戻り、フレーム
画像FRnに基づき、比較領域を決定する。この場合に
は、フレーム画像FRnの部分領域中、比較領域の移動
を示すベクトルが仮定条件を満たしていない部分領域の
みに関して、さらに、比較領域を決定する。より詳細に
は、図6に示すように、比較領域の移動を示す結果が仮
定条件を満たしていないような部分領域を、さらに分割
して、分割された小領域601ないし604を生成す
る。図6において、中央の破線で示す領域610が、前
回の処理にて得られた比較領域である。さらに、小領域
601のうち、上下左右に所定のドットずつ小さいサイ
ズの領域を比較領域611ないし614と決定する。こ
の実施の形態においては、24×24ドットのもとの部
分領域を2×2=4分割して、12×12ドットの画素
からなる部分領域を作成し、各部分領域中の、上下左右
に3ドットずつ小さいサイズの6×6ドットの画素から
なる領域を比較領域としている。
【0033】その後、分割された小領域ごとに、フレー
ムFRnの比較領域と、フレームFRn+1の対応するサ
イズの領域との間で、パターンマッチングを実行する。
この処理は、先に述べたステップ404およびステップ
405の処理と同様である。
ムFRnの比較領域と、フレームFRn+1の対応するサ
イズの領域との間で、パターンマッチングを実行する。
この処理は、先に述べたステップ404およびステップ
405の処理と同様である。
【0034】このようにして、分割された小領域のすべ
てに関して、比較領域が、時間的に隣接するフレーム中
の何れの位置に移動したかが決定されると、分割された
小領域のベクトルを統計的に処理して、さきの処理に用
いた、より大きな比較領域に関するベクトルを得る。こ
れは、上述したガウス分布を用いて、最も頻度の高い位
置での値を用いれば良い。
てに関して、比較領域が、時間的に隣接するフレーム中
の何れの位置に移動したかが決定されると、分割された
小領域のベクトルを統計的に処理して、さきの処理に用
いた、より大きな比較領域に関するベクトルを得る。こ
れは、上述したガウス分布を用いて、最も頻度の高い位
置での値を用いれば良い。
【0035】その後に、新たに得られたベクトルを含
む、全ての比較領域に関するベクトルが、先の仮定条件
を満足しているか否かを判定する(ステップ406)。
む、全ての比較領域に関するベクトルが、先の仮定条件
を満足しているか否かを判定する(ステップ406)。
【0036】上述したように、二度のパターンマッチン
グによっても、先の仮定条件を満足していないと判断さ
れた場合には、その領域に関する結果(移動方向および
移動距離)は、「空」となる。
グによっても、先の仮定条件を満足していないと判断さ
れた場合には、その領域に関する結果(移動方向および
移動距離)は、「空」となる。
【0037】前述したガウス分布中の最も頻度の高い位
置での、ベクトルの方向および大きさが、眼球移動ベク
トルと決定され、これにより、フレーム画像FRnとフ
レーム画像FRn+1との間の移動方向および移動距離が
決定される。このようにして得られた眼球移動ベクトル
(移動方向および移動距離)は、メモリ(図示せず)に
記憶される(ステップ407)。
置での、ベクトルの方向および大きさが、眼球移動ベク
トルと決定され、これにより、フレーム画像FRnとフ
レーム画像FRn+1との間の移動方向および移動距離が
決定される。このようにして得られた眼球移動ベクトル
(移動方向および移動距離)は、メモリ(図示せず)に
記憶される(ステップ407)。
【0038】画像メモリ120に記憶されたフレーム画
像の全てについて、ステップ402ないし407の処理
が終了すること(ステップ408でイエス(Yes))によ
り、フレーム画像解析部22による処理が終了する。
像の全てについて、ステップ402ないし407の処理
が終了すること(ステップ408でイエス(Yes))によ
り、フレーム画像解析部22による処理が終了する。
【0039】解析結果作成部24は、メモリに記憶され
た、時間的に隣接する各フレーム間に関する眼球移動ベ
クトル(移動方向および移動距離)のデータに基づき、
眼球の運動を示すグラフを作成する。移動距離および移
動方向に基づき、時間軸を横軸にとり、水平方向の移
動、垂直方向の移動、或いは、角度を縦軸にとることが
できる(図7参照)。たとえば、時間軸を横軸にとり、
かつ、水平方向の移動(右方向の移動を+(プラス)にす
る)を縦軸にとると、図8(a)に示すグラフを得るこ
とができ、その一方、時間軸を横軸にとり、かつ、垂直
方向の移動(上方向の移動を+(プラス)にする)を縦軸
にとると、図8(b)に示すグラフを得ることができ
る。このようにして解析結果作成部24により得られた
グラフは、表示装置28の画面上に表示され、或いは、
出力装置26から印刷物として出力される。
た、時間的に隣接する各フレーム間に関する眼球移動ベ
クトル(移動方向および移動距離)のデータに基づき、
眼球の運動を示すグラフを作成する。移動距離および移
動方向に基づき、時間軸を横軸にとり、水平方向の移
動、垂直方向の移動、或いは、角度を縦軸にとることが
できる(図7参照)。たとえば、時間軸を横軸にとり、
かつ、水平方向の移動(右方向の移動を+(プラス)にす
る)を縦軸にとると、図8(a)に示すグラフを得るこ
とができ、その一方、時間軸を横軸にとり、かつ、垂直
方向の移動(上方向の移動を+(プラス)にする)を縦軸
にとると、図8(b)に示すグラフを得ることができ
る。このようにして解析結果作成部24により得られた
グラフは、表示装置28の画面上に表示され、或いは、
出力装置26から印刷物として出力される。
【0040】また、時間的に隣接するフレーム間に関す
るベクトル表示(図9(a)参照)により、眼球の運動
の大きさ(速さ)およびその方向を表わすことが可能で
ある。これをグラフ化した場合に、時間と略等価である
フレーム番号を横軸にとり、かつ、ベクトルの方向を縦
軸にとると、図9(b)に示すグラフが得られ、その一
方、フレーム番号を横軸にとり、ベクトルの大きさ(移
動距離)を横軸にとると、図9(c)に示すグラフが得
られる。図9(b)においては、右回り方向をプラスに
するように角度が示されている。
るベクトル表示(図9(a)参照)により、眼球の運動
の大きさ(速さ)およびその方向を表わすことが可能で
ある。これをグラフ化した場合に、時間と略等価である
フレーム番号を横軸にとり、かつ、ベクトルの方向を縦
軸にとると、図9(b)に示すグラフが得られ、その一
方、フレーム番号を横軸にとり、ベクトルの大きさ(移
動距離)を横軸にとると、図9(c)に示すグラフが得
られる。図9(b)においては、右回り方向をプラスに
するように角度が示されている。
【0041】本実施の形態によれば、眼球のフレーム画
像を所定の領域に分割し、分割された部分領域の各々に
おいて、一方のフレームの比較領域と、他方のフレーム
の対応する大きさのフレームとのパターンマッチングを
実行し、得られた全ての比較領域の移動方向および移動
距離を統計的処理を施すことにより、時間的に隣接する
フレーム間での眼球の移動(移動方向および移動距離)
を決定している。すなわち、ブロックごとのマッチング
(ブロックマッチング法)を実行し、各ブロックに関す
る移動に基づき、眼球自体の移動を決定している。した
がって、眼球中の特徴点などを抽出する処理や操作など
を経ることなしに、眼球の移動を測定することができ
る。また、全ての比較領域の移動方向および移動距離に
統計的処理を施すため、ノイズや光の反射などが生じた
場合でも、適切に眼球の移動を測定することが可能とな
る。
像を所定の領域に分割し、分割された部分領域の各々に
おいて、一方のフレームの比較領域と、他方のフレーム
の対応する大きさのフレームとのパターンマッチングを
実行し、得られた全ての比較領域の移動方向および移動
距離を統計的処理を施すことにより、時間的に隣接する
フレーム間での眼球の移動(移動方向および移動距離)
を決定している。すなわち、ブロックごとのマッチング
(ブロックマッチング法)を実行し、各ブロックに関す
る移動に基づき、眼球自体の移動を決定している。した
がって、眼球中の特徴点などを抽出する処理や操作など
を経ることなしに、眼球の移動を測定することができ
る。また、全ての比較領域の移動方向および移動距離に
統計的処理を施すため、ノイズや光の反射などが生じた
場合でも、適切に眼球の移動を測定することが可能とな
る。
【0042】次に、本発明の第2の実施の形態にかかる
眼球運動解析装置につき説明を加える。この実施の形態
においては、回旋運動を含む画像の移動を確実に測定す
るようになっている。また、この実施の形態において、
画像運動解析装置は、フレーム画像解析部22での処理
の一部が相違することを除き、図1に示す第1の実施の
形態と同様である。このようの画像運動解析装置10の
作動を、より詳細に説明する。第1の実施の形態と同様
に、この実施の形態においても、VTR14に記録され
た眼球の画像信号が、A/D変換器16を介して、フレ
ームバッファ18に与えられ、さらに、フレームごとの
画像データが、画像メモリ20に記憶される。
眼球運動解析装置につき説明を加える。この実施の形態
においては、回旋運動を含む画像の移動を確実に測定す
るようになっている。また、この実施の形態において、
画像運動解析装置は、フレーム画像解析部22での処理
の一部が相違することを除き、図1に示す第1の実施の
形態と同様である。このようの画像運動解析装置10の
作動を、より詳細に説明する。第1の実施の形態と同様
に、この実施の形態においても、VTR14に記録され
た眼球の画像信号が、A/D変換器16を介して、フレ
ームバッファ18に与えられ、さらに、フレームごとの
画像データが、画像メモリ20に記憶される。
【0043】画像メモリに必要な画像データが記憶され
ると、制御部30により、フレームバッファに、所定の
フレームの画像データが読み出され、この画像データに
基づき、フレーム画像解析部22による処理が実行され
る。図10は、第2の実施の形態にかかるフレーム画像
解析部22による処理を示すフローチャートである。図
10において、この処理は、ステップ1006が付加さ
れたことを除き、図4に示すものと同様である。すなわ
ち、ステップ1001ないし1005は、図4のステッ
プ401ないし405に対応し、ステップ1007ない
しステップ1010は、図4のステップ406ないし4
09に対応する。
ると、制御部30により、フレームバッファに、所定の
フレームの画像データが読み出され、この画像データに
基づき、フレーム画像解析部22による処理が実行され
る。図10は、第2の実施の形態にかかるフレーム画像
解析部22による処理を示すフローチャートである。図
10において、この処理は、ステップ1006が付加さ
れたことを除き、図4に示すものと同様である。すなわ
ち、ステップ1001ないし1005は、図4のステッ
プ401ないし405に対応し、ステップ1007ない
しステップ1010は、図4のステップ406ないし4
09に対応する。
【0044】ステップ1001ないし1005の処理に
より、フレーム画像FRnを分割した部分領域および比
較領域が得られ、かつ、時間的に隣接するフレーム画像
FRn+1を分割した部分領域が得られた後に、フレーム
画像解析部22は、フレーム画像FRnの部分領域(図
5(a)の領域500)中の14×14ドットの比較領
域(図5(a)の領域501)と、フレーム画像FRn+
1の24×24ドットの対応する部分領域中の、比較領
域と同じサイズを有する14×14ドットの領域とのパ
ターンマッチングを、順次実行する。これにより、本実
施の形態においてはフレーム画像を25分割することに
より得られた25の部分領域中の、比較領域の各々の、
時間的に連続するフレーム画像間の移動方向および移動
距離を示すベクトルが得られる。 次いで、画像中の瞳
孔の位置を判定し、その移動を測定する。より詳細に
は、まず、フレーム画像FRnを構成する画素の画素値
(濃淡値)のヒストグラムを作成する。瞳孔は、眼球の
表面において最も黒色の強い部分であるため、黒色の近
傍のを示す画素値(濃淡値)、すなわち、ほぼ最大値と
等しい画素値(濃淡値)を閾値として、それより大きな
値を有する領域を、フレーム画像FRnから取り出すこ
とにより、瞳孔に略対応する領域を特定することが可能
となる。
より、フレーム画像FRnを分割した部分領域および比
較領域が得られ、かつ、時間的に隣接するフレーム画像
FRn+1を分割した部分領域が得られた後に、フレーム
画像解析部22は、フレーム画像FRnの部分領域(図
5(a)の領域500)中の14×14ドットの比較領
域(図5(a)の領域501)と、フレーム画像FRn+
1の24×24ドットの対応する部分領域中の、比較領
域と同じサイズを有する14×14ドットの領域とのパ
ターンマッチングを、順次実行する。これにより、本実
施の形態においてはフレーム画像を25分割することに
より得られた25の部分領域中の、比較領域の各々の、
時間的に連続するフレーム画像間の移動方向および移動
距離を示すベクトルが得られる。 次いで、画像中の瞳
孔の位置を判定し、その移動を測定する。より詳細に
は、まず、フレーム画像FRnを構成する画素の画素値
(濃淡値)のヒストグラムを作成する。瞳孔は、眼球の
表面において最も黒色の強い部分であるため、黒色の近
傍のを示す画素値(濃淡値)、すなわち、ほぼ最大値と
等しい画素値(濃淡値)を閾値として、それより大きな
値を有する領域を、フレーム画像FRnから取り出すこ
とにより、瞳孔に略対応する領域を特定することが可能
となる。
【0045】その後に、得られた瞳孔に略対応する領域
を含むような瞳孔測定用比較領域を決定する。この実施
の形態においては、瞳孔に略対応する領域を取り囲む、
これより僅かに大きな矩形が瞳孔測定用比較領域として
決定される。
を含むような瞳孔測定用比較領域を決定する。この実施
の形態においては、瞳孔に略対応する領域を取り囲む、
これより僅かに大きな矩形が瞳孔測定用比較領域として
決定される。
【0046】次いで、フレーム画像FRn+1において、
フレーム画像FRn+1中の瞳孔測定用比較領域に対応す
る領域、すなわち、等しい位置の領域を特定し、この領
域から上下左右に10ドットづつ大きいサイズの部分領
域を決定する。
フレーム画像FRn+1中の瞳孔測定用比較領域に対応す
る領域、すなわち、等しい位置の領域を特定し、この領
域から上下左右に10ドットづつ大きいサイズの部分領
域を決定する。
【0047】その後に、フレーム画像FRnの瞳孔測定
用比較領域と、フレーム画像FRn+1の部分領域中の、
瞳孔測定用比較領域と等しいサイズの領域とのパターン
マッチングが実行され、画素間の差の二乗和が最も小さ
いような位置に、瞳孔測定用比較領域が移動したと決定
される。
用比較領域と、フレーム画像FRn+1の部分領域中の、
瞳孔測定用比較領域と等しいサイズの領域とのパターン
マッチングが実行され、画素間の差の二乗和が最も小さ
いような位置に、瞳孔測定用比較領域が移動したと決定
される。
【0048】このようにして、瞳孔の移動距離および移
動方向、すなわち、瞳孔の移動を示す瞳孔移動ベクトル
が決定されると、ステップ1005にて得られた、各比
較領域の移動に関するベクトルの各々から、瞳孔移動ベ
クトルが減算され、得られたベクトルの各々が、略純粋
な回旋を示すベクトルとなる。これは、眼球の回旋運動
は、瞳孔を中心になされることに基づいている。このよ
うにして得られたベクトルは、図11に示すように、瞳
孔の位置Pを中心として、ほぼ円周方向の向きとなる。
動方向、すなわち、瞳孔の移動を示す瞳孔移動ベクトル
が決定されると、ステップ1005にて得られた、各比
較領域の移動に関するベクトルの各々から、瞳孔移動ベ
クトルが減算され、得られたベクトルの各々が、略純粋
な回旋を示すベクトルとなる。これは、眼球の回旋運動
は、瞳孔を中心になされることに基づいている。このよ
うにして得られたベクトルは、図11に示すように、瞳
孔の位置Pを中心として、ほぼ円周方向の向きとなる。
【0049】次いで、横軸にベクトルの方向、縦軸にそ
の頻度を設定して、得られた回旋を示すベクトル(移動
方向)の値をプロットする。この分布は、ガウス分布に
ほぼ従う。このため、横軸上に二つの閾値を設定し、第
1の閾値よりも小さい位置に表れたベクトル、或いは、
第2の閾値よりも大きい位置に表れたベクトル、すなわ
ち、分布にしたがって起こり得ないであろう方向のベク
トルを見つけ出す。このようなベクトルを得た部分領域
が仮定条件を満足していない領域となる。ステップ10
07においては、仮定条件を満足していないベクトルが
存在する場合には、ステップ1003に戻る。
の頻度を設定して、得られた回旋を示すベクトル(移動
方向)の値をプロットする。この分布は、ガウス分布に
ほぼ従う。このため、横軸上に二つの閾値を設定し、第
1の閾値よりも小さい位置に表れたベクトル、或いは、
第2の閾値よりも大きい位置に表れたベクトル、すなわ
ち、分布にしたがって起こり得ないであろう方向のベク
トルを見つけ出す。このようなベクトルを得た部分領域
が仮定条件を満足していない領域となる。ステップ10
07においては、仮定条件を満足していないベクトルが
存在する場合には、ステップ1003に戻る。
【0050】二度めに実行されるステップ1003にお
いては、図4のステップ403と同様に、比較領域の移
動を示す結果が、仮定条件を満たしいない部分領域をさ
らに分割して小領域を生成し、続くステップ1004お
よび1005において、小領域中に新たに決定された比
較領域の移動距離および移動方向を決定する。先程の処
理にて瞳孔の移動は測定されているため、これに基づ
き、新たな比較領域の回旋成分のみを抽出し(ステップ
1006)、新たに得られたベクトルを含む、すべての
比較領域に関するベクトルが、先の仮定条件を満たして
いるか否かを判定する(ステップ1007)。
いては、図4のステップ403と同様に、比較領域の移
動を示す結果が、仮定条件を満たしいない部分領域をさ
らに分割して小領域を生成し、続くステップ1004お
よび1005において、小領域中に新たに決定された比
較領域の移動距離および移動方向を決定する。先程の処
理にて瞳孔の移動は測定されているため、これに基づ
き、新たな比較領域の回旋成分のみを抽出し(ステップ
1006)、新たに得られたベクトルを含む、すべての
比較領域に関するベクトルが、先の仮定条件を満たして
いるか否かを判定する(ステップ1007)。
【0051】前述したガウス分布中の最も頻度の高い位
置でのベクトルが、眼球の回旋を示す眼球回旋ベクトル
と決定され、この眼球回旋ベクトルの方向(角度)によ
り、フレーム画像FRnとフレーム画像FRn+1との間
の移動方向(回旋角度)が決まる。このようにして得ら
れた移動方向(回旋角度)は、メモリ(図示せず)に記
憶される(ステップ1008)。
置でのベクトルが、眼球の回旋を示す眼球回旋ベクトル
と決定され、この眼球回旋ベクトルの方向(角度)によ
り、フレーム画像FRnとフレーム画像FRn+1との間
の移動方向(回旋角度)が決まる。このようにして得ら
れた移動方向(回旋角度)は、メモリ(図示せず)に記
憶される(ステップ1008)。
【0052】第1の実施の形態と同様に、画像メモリ1
20に記憶されたフレーム画像の全てについて、この処
理が実行される。
20に記憶されたフレーム画像の全てについて、この処
理が実行される。
【0053】また、解析結果作成部24は、メモリに記
憶された時間的に隣接する各フレーム間に関する眼球移
動ベクトルに対応する移動方向(回旋角度)のデータに
基づき、眼球の回旋を示すグラフを作成する。これは、
たとえば、時間軸を横軸にとり、回旋角度を縦軸にとれ
ば良い。
憶された時間的に隣接する各フレーム間に関する眼球移
動ベクトルに対応する移動方向(回旋角度)のデータに
基づき、眼球の回旋を示すグラフを作成する。これは、
たとえば、時間軸を横軸にとり、回旋角度を縦軸にとれ
ば良い。
【0054】この実施形態によれば、眼球の回旋の中心
をなす瞳孔の移動を示すベクトルを算出し、各比較領域
の移動を示すベクトルから、上記ベクトルを引くことに
より、略純粋な眼球の回旋運動を算出することが可能と
なる。さらに、移動方向(回旋角度)を、統計的処理に
より得るため、ノイズや光の反射などが生じた場合であ
っても、適切に眼球の回旋を測定することが可能とな
る。
をなす瞳孔の移動を示すベクトルを算出し、各比較領域
の移動を示すベクトルから、上記ベクトルを引くことに
より、略純粋な眼球の回旋運動を算出することが可能と
なる。さらに、移動方向(回旋角度)を、統計的処理に
より得るため、ノイズや光の反射などが生じた場合であ
っても、適切に眼球の回旋を測定することが可能とな
る。
【0055】本発明は、以上の実施の形態に限定される
ことなく、特許請求の範囲に記載された発明の範囲内
で、種々の変更が可能であり、それらも本発明の範囲内
に包含されるものであることは言うまでもない。
ことなく、特許請求の範囲に記載された発明の範囲内
で、種々の変更が可能であり、それらも本発明の範囲内
に包含されるものであることは言うまでもない。
【0056】たとえば、前記実施の形態においては、1
フレームあたり256×256ドットの画素データから
なる画像データを得て、これに基づき、120×120
ドットのフレーム画像を得ているが、これに限定される
ものではなく、512×512ドットの画像データか
ら、240×240ドットのフレーム画像を得るなど、
異なる解像度の画像データを用いて、異なる解像度のフ
レーム画像を得てもよいことは明らかである。
フレームあたり256×256ドットの画素データから
なる画像データを得て、これに基づき、120×120
ドットのフレーム画像を得ているが、これに限定される
ものではなく、512×512ドットの画像データか
ら、240×240ドットのフレーム画像を得るなど、
異なる解像度の画像データを用いて、異なる解像度のフ
レーム画像を得てもよいことは明らかである。
【0057】また、前記実施の形態においては、フレー
ム画像を25分割して、部分領域を得ているが、これに
限定するものではなく、分割数は、適宜決定すれば良
い。さらに、必ずしも、フレーム画像を分割して部分領
域を得ることはない。たとえば、その上下左右の端部が
互いにオーバーラップするような部分領域を得て、この
部分領域から比較領域を決定してもよい。
ム画像を25分割して、部分領域を得ているが、これに
限定するものではなく、分割数は、適宜決定すれば良
い。さらに、必ずしも、フレーム画像を分割して部分領
域を得ることはない。たとえば、その上下左右の端部が
互いにオーバーラップするような部分領域を得て、この
部分領域から比較領域を決定してもよい。
【0058】また、前記実施の形態においては、仮定条
件を満たしていないような部分領域において、当該部分
領域を4分割した小領域を生成し、この小領域内に得た
比較領域を用いてパターンマッチングを実行している
が、分割数および比較領域の大きさは、これに限定され
るものではない。
件を満たしていないような部分領域において、当該部分
領域を4分割した小領域を生成し、この小領域内に得た
比較領域を用いてパターンマッチングを実行している
が、分割数および比較領域の大きさは、これに限定され
るものではない。
【0059】さらに、仮定条件を満たしていない部分領
域がある場合に、当該部分領域において、相互に少なく
とも一画素以上離間している位置の画素をサンプリング
することにより比較領域を生成しても良い。たとえば、
一つおきに画素をサンプリングすれば、部分領域を小領
域に分割することなく、画素数の略半分の比較領域を得
ることになる。このように構成することにより、処理時
間を増大させることなく、画像中のノイズを除去して、
適切なベクトルを得ることが可能となる。
域がある場合に、当該部分領域において、相互に少なく
とも一画素以上離間している位置の画素をサンプリング
することにより比較領域を生成しても良い。たとえば、
一つおきに画素をサンプリングすれば、部分領域を小領
域に分割することなく、画素数の略半分の比較領域を得
ることになる。このように構成することにより、処理時
間を増大させることなく、画像中のノイズを除去して、
適切なベクトルを得ることが可能となる。
【0060】また、前記第1の実施の形態および第2の
実施の形態を組合わせて画像解析装置を構成しても良い
ことは明らかである。
実施の形態を組合わせて画像解析装置を構成しても良い
ことは明らかである。
【0061】さらに、前記第2の実施の形態において
は、比較領域ごとの回旋を示すベクトルの方向の分布の
みを考慮して、ガウス分布に従わないようなベクトルを
決定しているが、これに限定されるものではない。たと
えば、各比較領域のベクトルの長さ(移動距離)は、瞳
孔からベクトルまでの距離に比例しているため、この距
離を加味して、ベクトルの長さ(移動距離)の分布を考
え、ガウス分布に従わないベクトルを見つけ出すように
しても良い。
は、比較領域ごとの回旋を示すベクトルの方向の分布の
みを考慮して、ガウス分布に従わないようなベクトルを
決定しているが、これに限定されるものではない。たと
えば、各比較領域のベクトルの長さ(移動距離)は、瞳
孔からベクトルまでの距離に比例しているため、この距
離を加味して、ベクトルの長さ(移動距離)の分布を考
え、ガウス分布に従わないベクトルを見つけ出すように
しても良い。
【0062】また、前記実施の形態においては、赤外線
CCDカメラからの画像をアナログVTRに記録し、記
録された画像信号を、A/D変換器によりディジタル化
して、ディジタル画像データを得ているがこれに限定さ
れるものではない。たとえば、ディジタルスチルカメラ
を用いて、所定の時間ごとにシャッタを作動させること
により、ディジタル画像データを得てもよいし、ディジ
タルビデオカメラを用いても良い。
CCDカメラからの画像をアナログVTRに記録し、記
録された画像信号を、A/D変換器によりディジタル化
して、ディジタル画像データを得ているがこれに限定さ
れるものではない。たとえば、ディジタルスチルカメラ
を用いて、所定の時間ごとにシャッタを作動させること
により、ディジタル画像データを得てもよいし、ディジ
タルビデオカメラを用いても良い。
【0063】また、本明細書において、本明細書におい
て、手段とは必ずしも物理的手段を意味するものではな
く、各手段の機能が、ソフトウェアによって実現される
場合も包含する。また、一つの手段或いは部材の機能
が、二つ以上の物理的手段或いは部材により実現されて
も、若しくは、二つ以上の手段或いは部材の機能が、一
つの物理的手段或いは部材により実現されてもよい。
て、手段とは必ずしも物理的手段を意味するものではな
く、各手段の機能が、ソフトウェアによって実現される
場合も包含する。また、一つの手段或いは部材の機能
が、二つ以上の物理的手段或いは部材により実現されて
も、若しくは、二つ以上の手段或いは部材の機能が、一
つの物理的手段或いは部材により実現されてもよい。
【0064】
【発明の効果】本発明によれば、画像に突発的なノイズ
が生じていた場合であっても、正確に、眼球運動を測定
可能であり、かつ、その構成も比較的単純な眼球運動解
析装置を提供することが可能となる。
が生じていた場合であっても、正確に、眼球運動を測定
可能であり、かつ、その構成も比較的単純な眼球運動解
析装置を提供することが可能となる。
【図1】 図1は、本発明の実施の形態にかかる眼球運
動解析装置の構成を示すブロックダイヤグラムである。
動解析装置の構成を示すブロックダイヤグラムである。
【図2】 図2は、本発明にかかる眼球運動の測定の原
理を説明するための図である。
理を説明するための図である。
【図3】 図3は、本発明にかかる眼球運動の測定の原
理を説明するための図である。
理を説明するための図である。
【図4】 図4は、第1の実施の形態にかかるフレーム
画像解析部の処理を示すフローチャートである。
画像解析部の処理を示すフローチャートである。
【図5】 図5は、第1の実施の形態にかかる、部分領
域中の比較領域の決定およびパターンマッチングを説明
するための図である。
域中の比較領域の決定およびパターンマッチングを説明
するための図である。
【図6】 図6は、本実施の形態にかかる部分領域中の
小領域の決定を説明するための図である。
小領域の決定を説明するための図である。
【図7】 図7は、本実施の形態にかかる処理により得
られた結果を表示するためのグラフを説明するための図
である。
られた結果を表示するためのグラフを説明するための図
である。
【図8】 図8は、本実施の形態にかかる表示装置の画
面上にCRTに表示され、或いは、出力装置に印刷物と
して出力されたグラフの一例を示す図である。
面上にCRTに表示され、或いは、出力装置に印刷物と
して出力されたグラフの一例を示す図である。
【図9】 図9は、本実施の形態にかかる処理に得られ
た結果を表示するグラフの一例を示す図である。
た結果を表示するグラフの一例を示す図である。
【図10】 図10は、第2の実施の形態にかかるフレ
ーム画像解析部の処理を示すフローチャートである。
ーム画像解析部の処理を示すフローチャートである。
【図11】 図11は、第2の実施の形態にかかる処理
により得られた眼球の回旋運動を説明するための図であ
る。
により得られた眼球の回旋運動を説明するための図であ
る。
10 眼球運動解析装置 12 赤外線CCDカメラ 14 VTR 16 A/D変換器 18 フレームバッファ 20 画像メモリ 22 フレーム画像解析部 24 解析結果表示部 26 出力装置 28 表示装置 30 制御部
フロントページの続き (56)参考文献 特開 平4−101126(JP,A) 特開 平7−275206(JP,A) 特開 平5−314262(JP,A) 特開 平5−312819(JP,A) (58)調査した分野(Int.Cl.6,DB名) A61B 3/113
Claims (8)
- 【請求項1】 カメラにて眼球を撮影したフレーム画像
に基づき、眼球の運動を測定する眼球運動解析装置であ
って、 前記フレーム画像から複数の部分領域を得て、各部分領
域中の比較対象とすべき比較領域を決定する比較領域決
定手段と、 前記フレーム画像の比較領域と、時間的に隣接する隣接
フレーム画像の対応する部分領域中の、前記比較領域と
同じ大きさの領域との間でパターンマッチングを実行し
て、前記比較領域の、隣接フレーム画像中の位置を決定
するマッチング実行手段と、 マッチング実行手段により得られた、前記比較領域の、
隣接フレーム画像中の位置に基づき、前記比較領域の移
動方向および移動距離からなるベクトルを生成するベク
トル生成手段と、 ベクトル生成手段により得られた、比較領域の各々に関
するベクトルを統計的に処理して、隣接するフレーム画
像の間の眼球の移動を示す眼球移動ベクトルを生成する
統計処理手段と、 統計処理手段により得られた眼球移動ベクトルに基づ
き、眼球の移動を示すグラフを示すグラフ生成手段とを
備え、当該グラフが出力されるように構成され、かつ、 前記統計処理手段が、前記複数のベクトルの移動方向お
よび/または移動距離に関する分布を生成し、前記分布
が所定の分布にしたがっているか否かを判断し、分布に
したがっていないベクトルがある場合には、当該ベクト
ルに関連する部分領域について、前記比較領域決定手段
が、再度比較領域を決定し、当該比較領域に関するパタ
ーンマッチングおよびベクトルの生成が実行されるよう
に構成されたことを特徴とする眼球運動解析装置。 - 【請求項2】 前記分布にしたがっていないベクトルが
ある場合に、前記比較領域決定手段が、関連する部分領
域をさらに小領域に分割して、各小領域に含まれる第2
の比較領域を生成するように構成されたことを特徴とす
る請求項1に記載の眼球運動解析装置。 - 【請求項3】 さらに、前記部分領域を分割することに
より得られた第2の比較領域のそれぞれの移動を示すベ
クトルを統計的に処理し、前記部分領域に関連するもの
との比較領域を示すベクトルを得る小領域ベクトル統制
処理手段を有することを特徴とする請求項2に記載の画
像運動解析装置。 - 【請求項4】 前記分布にしたがっていないベクトルが
ある場合に、前記比較領域決定手段が、前記部分領域に
おいて、少なくとも相互に一画素以上離間して位置する
画素をサンプリングすることにより、比較領域を得るよ
うに構成されたことを特徴とする請求項1に記載の画像
運動解析装置。 - 【請求項5】 さらに、フレーム画像中の瞳孔の位置を
決定する瞳孔位置決定手段と、隣接するフレーム間の瞳
孔の移動を決定し、瞳孔の移動を示す瞳孔移動ベクトル
を生成する瞳孔移動測定手段と、 前記ベクトル生成手段により生成された、各比較領域に
関するベクトルから、瞳孔移動ベクトルを減算し、眼球
の回旋を示す比較領域回旋ベクトルを生成する回旋ベク
トル生成手段とを備え、 前記統計処理手段が、回旋ベクトル生成手段により得ら
れた、比較領域の各々の回旋に関するベクトルを統計的
に処理して、隣接するフレーム画像間の眼球の回旋を示
す眼球回旋ベクトルを生成するように構成されたことを
特徴とする請求項1ないし4の何れか一項に記載の眼球
運動解析装置。 - 【請求項6】 カメラにて眼球を撮影したフレーム画像
に基づき、眼球の運動を測定する眼球運動解析装置であ
って、 前記フレーム画像から複数の部分領域を得て、各部分領
域中の比較対象とすべき比較領域を決定する比較領域決
定手段と、 前記フレーム画像の比較領域と、時間的に隣接する隣接
フレーム画像の対応する部分領域中の、前記比較領域と
同じ大きさの領域との間でパターンマッチングを実行し
て、前記比較領域の、隣接フレーム画像中の位置を決定
するマッチング実行手段と、 マッチング実行手段により得られた、前記比較領域の、
隣接フレーム画像中の位置に基づき、前記比較領域の移
動方向および移動距離からなるベクトルを生成するベク
トル生成手段と、 フレーム画像中の瞳孔の位置を決定する瞳孔位置決定手
段と、 隣接するフレーム間の瞳孔の移動を決定し、瞳孔の移動
を示す瞳孔移動ベクトルを生成する瞳孔移動測定手段
と、 前記ベクトル生成手段により生成された、各比較領域に
関するベクトルから、瞳孔移動ベクトルを減算し、眼球
の回旋を示す比較領域回旋ベクトルを生成する回旋ベク
トル生成手段と、 回旋ベクトル生成手段により得られた、比較領域の各々
の回旋に関するベクトルを統計的に処理して、隣接する
フレーム画像間の眼球の回旋を示す眼球回旋ベクトルを
生成する統計処理手段と、 統計処理手段により得られた眼球回旋ベクトルに基づ
き、眼球の回旋を示すグラフを示すグラフ生成手段とを
備え、当該グラフが出力されるように構成されたことを
特徴とする画像運動解析装置。 - 【請求項7】 カメラにて眼球を撮影したフレーム画像
に基づき、眼球の運動を測定する眼球運動解析方法であ
って、 前記フレーム画像から複数の部分領域を得て、各部分領
域中の比較対象とすべき比較領域を決定するステップ
と、 前記フレーム画像の比較領域と、時間的に隣接する隣接
フレーム画像の対応する部分領域中の、前記比較領域と
同じ大きさの領域との間でパターンマッチングを実行し
て、前記比較領域の、隣接フレーム画像中の位置を決定
するステップと、 パターンマッチングにより得られた、前記比較領域の、
隣接フレーム画像中の位置に基づき、前記比較領域の移
動方向および移動距離からなるベクトルを生成するステ
ップと、 前記ベクトル生成ステップにより得られた、比較領域の
各々に関するベクトルを統計的に処理して、隣接するフ
レーム画像の間の眼球の移動を示す眼球移動ベクトルを
生成するステップと、 前記眼球移動ベクトルに基づき、眼球の移動を示すグラ
フを示すステップとを備え、 前記統計的な処理処理が、前記複数のベクトルの移動方
向および/または移動距離に関する分布を生成し、前記
分布が所定の分布にしたがっているか否かを判断し、 分布にしたがっていないベクトルがある場合には、当該
ベクトルに関連する部分領域について、再度比較領域を
決定し、当該比較領域に関するパターンマッチングおよ
びベクトルの生成が実行されるように構成されたことを
特徴とする眼球運動解析方法。 - 【請求項8】 カメラにて眼球を撮影したフレーム画像
に基づき、眼球の運動を測定するためのプログラムを格
納した、コンピュータにより読み出し可能な記憶媒体で
あって、 前記フレーム画像から複数の部分領域を得て、各部分領
域中の比較対象とすべき比較領域を決定するステップ
と、 前記フレーム画像の比較領域と、時間的に隣接する隣接
フレーム画像の対応する部分領域中の、前記比較領域と
同じ大きさの領域との間でパターンマッチングを実行し
て、前記比較領域の、隣接フレーム画像中の位置を決定
するステップと、 パターンマッチングにより得られた、前記比較領域の、
隣接フレーム画像中の位置に基づき、前記比較領域の移
動方向および移動距離からなるベクトルを生成するステ
ップと、 前記ベクトル生成ステップにより得られた、比較領域の
各々に関するベクトルを統計的に処理して、隣接するフ
レーム画像の間の眼球の移動を示す眼球移動ベクトルを
生成するステップと、 前記眼球移動ベクトルに基づき、眼球の移動を示すグラ
フを示すステップとを備え、 前記統計的な処理処理が、前記複数のベクトルの移動方
向および/または移動距離に関する分布を生成し、前記
分布が所定の分布にしたがっているか否かを判断し、 分布にしたがっていないベクトルがある場合には、当該
ベクトルに関連する部分領域について、再度比較領域を
決定し、当該比較領域に関するパターンマッチングおよ
びベクトルの生成が実行されるように構成されたプログ
ラムを格納した、コンピュータにより読み出し可能な記
憶媒体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP8309095A JP2964132B2 (ja) | 1996-11-20 | 1996-11-20 | 眼球運動解析装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP8309095A JP2964132B2 (ja) | 1996-11-20 | 1996-11-20 | 眼球運動解析装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH10146319A JPH10146319A (ja) | 1998-06-02 |
JP2964132B2 true JP2964132B2 (ja) | 1999-10-18 |
Family
ID=17988833
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP8309095A Expired - Fee Related JP2964132B2 (ja) | 1996-11-20 | 1996-11-20 | 眼球運動解析装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2964132B2 (ja) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4687247B2 (ja) * | 2005-05-31 | 2011-05-25 | パナソニック株式会社 | 異物検出方法及び眼球観察装置 |
JP5858603B2 (ja) * | 2010-03-12 | 2016-02-10 | キヤノン株式会社 | 眼科装置及びその制御方法 |
US8998412B2 (en) | 2010-03-12 | 2015-04-07 | Canon Kabushiki Kaisha | Ophthalmologic apparatus and control method for the same |
JP5635898B2 (ja) * | 2010-12-17 | 2014-12-03 | キヤノン株式会社 | 眼底撮像装置及びその制御方法 |
CN105324794A (zh) * | 2013-04-10 | 2016-02-10 | 奥克兰联合服务有限公司 | 头部和眼睛跟踪 |
EP3371973B1 (en) * | 2015-11-06 | 2023-08-09 | Facebook Technologies, LLC | Eye tracking using a diffraction pattern of coherent light on the surface of the eye |
-
1996
- 1996-11-20 JP JP8309095A patent/JP2964132B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JPH10146319A (ja) | 1998-06-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3679512B2 (ja) | 画像抽出装置および方法 | |
JP4130820B2 (ja) | 顔中心位置検出装置及び方法並びにプログラム | |
JP6048844B2 (ja) | 眼球回旋測定装置、眼球回旋測定方法、及び、眼球回旋測定プログラム | |
JP2005210288A (ja) | 撮像装置および方法並びにプログラム | |
JPH1031745A (ja) | 自動画像解析方法 | |
JP4513699B2 (ja) | 動画像編集装置、動画像編集方法及びプログラム | |
JP2007241579A (ja) | 特徴点検出装置及びその方法 | |
JP2000298300A (ja) | 手ぶれ画像補正方法、記録媒体及び撮像装置 | |
JP2018025970A (ja) | 画像処理装置、画像処理方法及びプログラム | |
KR100684301B1 (ko) | 화상처리장치 및 방법 | |
JP3993029B2 (ja) | 化粧シミュレーション装置、化粧シミュレーション方法、化粧シミュレーションプログラム及びそのプログラムを記録した記録媒体 | |
CN114286064A (zh) | 一种实时对焦方法、装置、***和计算机可读存储介质 | |
JP2964132B2 (ja) | 眼球運動解析装置 | |
JP2009244944A (ja) | 画像回復装置および撮影装置 | |
EP2693397B1 (en) | Method and apparatus for noise reduction in an imaging system | |
JP2004157778A (ja) | 鼻位置の抽出方法、およびコンピュータに当該鼻位置の抽出方法を実行させるためのプログラムならびに鼻位置抽出装置 | |
JP3963789B2 (ja) | 眼検出装置、眼検出プログラム、そのプログラムを記録する記録媒体及び眼検出方法 | |
JPH11225967A (ja) | 眼球運動解析方法、眼球運動表示診断装置及び個人診断ファイルを記録した記録媒体 | |
JP3789911B2 (ja) | 画像処理装置および顔画像処理装置 | |
JPH09327037A (ja) | 輪郭抽出方法及び画像抽出方法及び装置 | |
JP3571405B2 (ja) | 角膜内皮細胞密度算出方法 | |
JPH04276232A (ja) | 画像処理方法及びこれを用いたシステム | |
KR102669621B1 (ko) | 지능형 고해상도 내시경 영상을 획득하여 저장하는 영상 제공 방법 및 장치 | |
Corkidi et al. | Automated in vivo and online morphometry of human corneal endothelium | |
JP2001092978A (ja) | 人物像の姿勢推定装置および人物像の姿勢推定プログラムを記録した記録媒体 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20070813 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080813 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090813 Year of fee payment: 10 |
|
LAPS | Cancellation because of no payment of annual fees |