(57)【特許請求の範囲】
【請求項１】
眼振データにおいてファジー理論を適用して緩徐相を抽出する方法であって、
ｎ（ｎは２以上の整数）個目の眼振データと（ｎ−１）個目の眼振データとの差分データを算出するステップ、
前記差分データの絶対値に基づいてメンバーシップ値を決定するステップ、及び
前記メンバーシップ値を所定の閾値と比較することによって緩徐相を抽出するステップ、
を有することを特徴とする緩徐相抽出方法。
【請求項２】
前記眼振データが、回転刺激が与えられている場合の眼振データであるか、回転刺激が与えられていない場合の眼振データであるかを選択するステップ、を更に有し、
前記選択するステップにおいて回転刺激が与えられていない場合の眼振データが選択された場合に、前記算出するステップ、前記決定するステップ、及び前記抽出するステップが実行される
ことを特徴とする請求項１に記載の緩徐相抽出方法。
【請求項３】
眼振データにおいてファジー理論を適用して緩徐相を抽出する眼振データ解析システムであって、
ｎ（ｎは２以上の整数）個目の眼振データと（ｎ−１）個目の眼振データとの差分データを算出し、前記差分データの絶対値に基づいてメンバーシップ値を決定し、前記メンバーシップ値を所定の閾値と比較することによって緩徐相を抽出するデータ処理部
を有することを特徴とする眼振データ解析システム。
【請求項４】
前記データ処理部は、前記眼振データが、回転刺激が与えられていない場合の眼振データである場合に、ｎ（ｎは２以上の整数）個目の眼振データと（ｎ−１）個目の眼振データとの差分データを算出し、前記差分データの絶対値に基づいてメンバーシップ値を決定し、前記メンバーシップ値を所定の閾値と比較することによって緩徐相を抽出する
よう構成されていることを特徴とする請求項３に記載の眼振データ解析システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、眼振データにおける緩徐相抽出方法及びそれを用いた眼振データ解析システムに関する。
【背景技術】
【０００２】
医学的な眼球運動解析手法の１つに、前庭動眼反射を用いた解析手法がある。前庭動眼反射とは、頭部を回転した場合に、その回転とは反対方向に眼球が動く反射のことである。一定方向の加速度刺激が持続すると、眼球は逆方向に徐々に変位し、急速に元に戻るという運動が規則的に反復する。
【０００３】
眼球運動において、徐々に変位するなめらかな成分は「緩徐相」、急速に元に戻る成分は「急速相」と呼ばれる。これら２つの相が交互に現れる現象は眼振として知られている。この眼振の様子を観測した画像系列から得られたデータを、以下「眼振データ」と言う。
前庭動眼反射解析のためには、小脳等の中枢の働きによって引き起こされる急速相成分を除去し、前庭動眼反射によって引き起こされる緩徐相成分のみを抽出する必要がある。
【０００４】
従来、眼振データから急速相成分を除去し、緩徐相成分を抽出するための種々の手法が提案されてきた。例えば、非特許文献１に、眼振データにおいてファジー集合理論を適用して緩徐相を自動的に抽出する手法が開示されている。非特許文献１の手法等、従来の手法においては、被験者の頭部へ回転刺激（正弦波刺激）を与えた場合の眼振を観測した眼振データに対して、正弦波近似曲線に基づいてファジー集合理論を適用しメンバーシップ値を決定することによって緩徐相を抽出するものであった。
【０００５】
【非特許文献１】Mohammad Arzi, Michel Magnin: "A Fuzzy Set Theoretical Approach to Automatic Analysis of Nystagmic Eye Movements" IEEE Trans. On Biomedical Engineering, vol.36, NO.9, pp.954-963
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかしながら、眼振は、被験者の頭部へ回転刺激（正弦波刺激）を与えた場合以外に、三半規管や脳の障害等の疾患等によっても引き起こされる。上記従来の手法は、被験者の頭部へ回転刺激が与えられている場合の眼振を観測した眼振データから緩徐相を抽出することはできるが、回転刺激が与えられていない場合の眼振を観測した眼振データからは精度良く緩徐相の抽出を行うことができない、という問題があった。
【０００７】
本発明は、上記問題に鑑み、回転刺激が与えられていない場合の眼振データからも高精度に緩徐相を抽出する、信頼性の高い緩徐相抽出方法及び眼振データ解析システムを提供することを課題とする。
【課題を解決するための手段】
【０００８】
上記問題を解決するために、本発明は以下の構成を有する。
本発明の１つの観点による緩徐相抽出方法は、眼振データにおいてファジー理論を適用して緩徐相を抽出する方法であって、ｎ（ｎは２以上の整数）個目の眼振データと（ｎ−１）個目の眼振データとの差分データを算出するステップ、前記差分データの絶対値に基づいてメンバーシップ値を決定するステップ、及び前記メンバーシップ値を所定の閾値と比較することによって緩徐相を抽出するステップ、を有する。
所定の閾値は、任意の閾値決定方法により決定され得るが、例えば、Kittlerの方法が高精度であるため好ましい。
この発明によれば、回転刺激が与えられている場合の眼振データのみならず、回転刺激が与えられていない場合の眼振データからも高精度に緩徐相を抽出する、信頼性の高い緩徐相抽出方法を実現できる。
【０００９】
好ましくは、前記眼振データが、回転刺激が与えられている場合の眼振データであるか、回転刺激が与えられていない場合の眼振データであるかを選択するステップ、を更に有し、前記選択するステップにおいて回転刺激が与えられていない場合の眼振データが選択された場合に前記算出するステップ、前記決定するステップ、及び前記抽出するステップが実行される。
この発明によれば、回転刺激が与えられている場合の眼振データのみならず、回転刺激が与えられていない場合の眼振データからも高精度に緩徐相を抽出する、信頼性の高い緩徐相抽出方法を実現できる。
【００１０】
本発明の１つの観点による眼振データ解析システムは、眼振データにおいてファジー理論を適用して緩徐相を抽出する眼振データ解析システムであって、ｎ（ｎは２以上の整数）個目の眼振データと（ｎ−１）個目の眼振データとの差分データを算出し、前記差分データの絶対値に基づいてメンバーシップ値を決定し、前記メンバーシップ値を所定の閾値と比較することによって緩徐相を抽出するデータ処理部を有する。
所定の閾値は、任意の閾値決定方法により決定され得るが、例えば、Kittlerの方法が高精度であるため好ましい。
この発明によれば、回転刺激が与えられている場合の眼振データのみならず、回転刺激が与えられていない場合の眼振データからも高精度に緩徐相を抽出する、信頼性の高い眼振データ解析システムを実現できる。
【００１１】
好ましくは、前記データ処理部は、前記眼振データが、回転刺激が与えられていない場合の眼振データである場合に、ｎ（ｎは２以上の整数）個目の眼振データと（ｎ−１）個目の眼振データとの差分データを算出し、前記差分データの絶対値に基づいてメンバーシップ値を決定し、前記メンバーシップ値を所定の閾値と比較することによって緩徐相を抽出するよう構成されている。
この発明によれば、回転刺激が与えられている場合の眼振データのみならず、回転刺激が与えられていない場合の眼振データからも高精度に緩徐相を抽出する、信頼性の高い眼振データ解析システムを実現できる。
【発明の効果】
【００１２】
本発明によれば、回転刺激が与えられている場合の眼振データのみならず、回転刺激が与えられていない場合の眼振データからも高精度に緩徐相を抽出する、信頼性の高い緩徐相抽出方法及び眼振データ解析システムを実現できる、という有利な効果を奏する。
【発明を実施するための最良の形態】
【００１３】
以下本発明の実施をするための最良の形態を具体的に示した実施例について、図面とともに記載する。
【実施例】
【００１４】
図１〜図１９を参照して、本発明の実施例における緩徐相抽出方法を用いた眼振データ解析システムについて説明する。
まず、図１を参照して、本実施例における緩徐相抽出方法を用いた眼振データ解析システムで扱う眼振データについて説明する。
【００１５】
図１は、眼振データにおける緩徐相及び急速相を説明するための図である。眼振データには、通常、図１に示すように、緩やかな変位である「緩徐相」（符号１で示す）、及び、緩徐相とは反対の方向に急激に元に戻る変位である「急速相」（符号２で示す）がある。
【００１６】
本実施例において、被験者の頭部へ回転刺激（正弦波刺激）を与えた場合の眼振を観察した眼振データ、及び、回転刺激が与えられていない場合の眼振を観察した眼振データの２種類の眼振データを取り扱う。両者は、後述するそれぞれの条件の下で被験者の眼球運動をビデオカメラで撮影した動画像から秒間３０フレームで画像系列を獲得し、そこから眼振データを算出することで得られる。（即ち、１秒間にそれぞれ３０個の眼振眼位データが得られる。）
【００１７】
本実施例における眼振データとして、画像系列の各画像における眼球の位置を示す眼振眼位データ（以下、単に「眼位データ」と記す）と、眼振時の眼球の変位の速度を示す眼振速度データ（以下、単に「速度データ」と記す）とを用いる。
【００１８】
眼位データを算出する方法として、例えば、ある瞬間の画像から虹彩紋裡のパターンを抽出し、そのパターンと正面視時の画像のパターンとを比較することで眼位の差を求める方法を用いる。速度データを算出する方法として、例えば、画像系列において連続する２つの画像間での眼球の位置の変位を測定し、それを各データ間の時間（秒間３０フレームの場合、（１／３０）秒）で割ることによって求める方法を用いる。
【００１９】
被験者の頭部へ回転刺激を与えた場合の眼振データについて説明する。
被験者の頭部を、６０度後屈させ、かつ、右又は左の方向に４５度回旋させた姿勢に保ち、時計回り又は反時計回りに被験者を回転させる。このようにして被験者の頭部へ回転刺激を与えることにより、半規管を刺激して眼振を誘発する。回転中に生じる眼球運動をビデオカメラで撮影する。
回転刺激の周波数は、例えば、０．１Ｈｚ（すなわち、１０秒間に１回転）と０．５Ｈｚ（すなわち、２秒間に１回転）である。一般的に、このように被験者の頭部に回転刺激を与えた場合、得られる眼振データ中の緩徐相成分は、回転刺激の周波数に近い周波数の正弦波状になることが知られている。特に速度データに関して、この特徴が顕著に現れる。
【００２０】
回転刺激が与えられていない場合の眼振データについて説明する。
半規管に疾患を有する良性発作性頭位めまい症（BPPV : Benign Paroxysmal Positional Vertigo）患者を対象に、後半規管の面で被験者の頭位を最初に１回だけ変化させて、眼振を誘発する。その後生じる眼球運動をビデオカメラで撮影する。
【００２１】
図２及び図１９を参照して、本実施例における緩徐相抽出方法を用いた眼振データ解析システムの構成について説明する。図１９は、本実施例における緩徐相抽出方法を用いた眼振データ解析システムの構成を示すブロック図である。図２は、本実施例における緩徐相抽出方法を用いた眼振データ解析システムの処理の流れを示すフローチャートである。
【００２２】
図１９に示す、本実施例における緩徐相抽出方法を用いた眼振データ解析システムは、眼振データ抽出部２０００、操作部２００１、データ処理部２００２、記憶部２００３、画像処理部２００４、及び表示部２００５を有する。
【００２３】
眼振データ抽出部２０００は、撮影された画像データを入力して、眼振データを抽出する。
操作部２００１は、操作者が、撮影された画像データの種別、即ち、回転刺激が与えられている場合のデータであるかどうか、あるいは、眼位データであるか速度データであるか等を入力するための操作入力部である。
【００２４】
データ処理部２００２は、操作部２００１から入力した眼振データ種別を用いて、記憶部２００３からそれぞれの種別に対応するフローチャート（処理手順）を読み出し、眼振データ抽出部２０００から入力した眼振データに対して緩徐相を抽出するための処理を行う（詳細は後述）。
【００２５】
画像処理部２００４は、データ処理部２００２が抽出した緩徐相を眼振データ抽出部２０００が抽出した元の眼振データに重畳して、表示画像データを作成する。具体的には、眼振データをグラフ化した画像に、緩徐相の部分を強調した画像を重ね合わせる等により表示画像データを作成する。
表示部２００５は、画像処理部２００４が作成した表示画像データをディスプレイ等に出力する。操作者（例えば、医師等）は、表示部２００５に表示された緩徐相抽出結果である画像データによって、BPPVを診断する。
【００２６】
本実施例において、後半規管の面で被験者の頭位を最初に１回だけ変化させた場合の眼振データの緩徐相を抽出することによって、BPPVを診断しているが、これに限らず、他の疾患によって出現する眼振（注視眼振、自発眼振等）の特徴に基づいて別の方法によって眼振を誘発することによって、別の疾患を診断することもできる。
【００２７】
図２を参照して、本実施例における緩徐相抽出方法を用いた眼振データ解析システムの処理の流れについて説明する。
ステップＳ２０１において、眼振データ抽出部２０００は、眼球運動をビデオカメラで撮影して得られた画像系列から眼振データを抽出する。
ステップＳ２０２において、データ処理部２００２は、操作部２００１から入力された眼振データ種別に基づいて、抽出された眼振データが回転刺激を与えられた場合の眼振データであるか否かを調べる。眼振データが回転刺激を与えられた場合の眼振データである場合、ステップＳ２０３に進む。そうでない場合は、ステップＳ２０６に進む。
【００２８】
ステップＳ２０３において、データ処理部２００２は、操作部２００１からの眼振データ種別に基づいて、抽出された眼振データが速度データであるか眼位データであるかを調べる。眼振データが速度データである場合は、ステップＳ２０４に進む。眼位データである場合は、ステップＳ２０５に進む。
【００２９】
各ステップＳ２０４、Ｓ２０５、及びＳ２０６において、データ処理部２００２は、眼振データ種別に基づいて記憶部２００３から読み出したフローチャート（後に図３、図８、及び図９を用いて説明する。）を実行し、緩徐相を抽出する。ステップＳ２０４、Ｓ２０５、及びＳ２０６における緩徐相抽出方法についての詳細な説明は後述する。
【００３０】
ステップＳ２０７において、画像処理部２００４は、各ステップＳ２０４、Ｓ２０５、及びＳ２０６で抽出された緩徐相を眼振データ抽出部２０００が抽出した元の眼振データに重畳して、表示画像を作成する。
【００３１】
ステップＳ２０８において、表示部２００５は、画像処理部２００４が作成した表示画像データを出力する。
【００３２】
図３〜図９を参照して、本実施例の緩徐相抽出方法について説明する。
まず、図３〜図７を参照して、ステップＳ２０４の回転刺激が与えられている場合の速度データの緩徐相抽出方法について、詳細に説明する。図３は、回転刺激が与えられている場合の速度データについての緩徐相抽出方法を示すフローチャートである。
【００３３】
ステップＳ３０１において、速度データ中の急激なスパイクを除去する。
例えば、ローパスフィルタを用いる方法がある。元データｙｉとローパスフィルタ通過後のデータｙｉ’との差を取り、ノイズ成分ｑｉとする。ノイズ成分ｑｉは、以下の式（１）で示される。
ｑｉ＝ｙｉ−ｙｉ’  （ｉ＝２，３，・・・，Ｎ）  ・・・・・・  （１）
【００３４】
各ノイズ成分ｑｉの絶対値と閾値ｑｒｅｆとを比較し、次式（２）を満たすデータを急激なスパイクと見なして除去する。
｜ｑｉ｜＞ｑｒｅｆ    ・・・・・・  （２）
【００３５】
ステップＳ３０２において、ステップＳ３０１でスパイクの除去を行った後のデータを、複数の断片に周期分割する。
本実施例では、まず、周期分割のために次式（３）及び（４）で表される曲線ｆ１（ｘｉ）を定義する。次式（３）において、Ａ、Ｂ、及びＣはフィッティングのためのパラメータである。ωは、１秒当たり３０個のデータが得られる場合における、回転刺激の周波数Ｆに対する角周波数であり、次式（４）で表される。
ｆ１（ｘｉ）＝Ａ・ｓｉｎ（ω・ｘｉ）＋Ｂ・ｃｏｓ（ω・ｘｉ）＋Ｃ  ・・・・・・  （３）
ω＝２πＦ／３０  ・・・・・・  （４）
【００３６】
速度データ（ｘｉ，ｙｉ）に対して、次式（５）によって示された、定義した曲線ｆ１（ｘｉ）と各データｙｉとの二乗誤差の総和Sが最小となるように、パラメータＡ、Ｂ、及びＣを決定する。Ｎは速度データの全データ数を示す。
【数１】
0004681377.tif 000002
【００３７】
例えば、図４に示される速度データに対して上記の処理を行い、上記式（３）のパラメータＡ、Ｂ、及びＣを決定した場合の曲線を図５に示す。この曲線の周期に従って速度データを分割する。図５の曲線の周期に従って速度データを分割した様子を図６に示す。分割された断片毎に以降のファジー集合理論を用いた解析を行うことで、眼振データが図４のように正弦波形が上下に変位している場合でも、緩徐相の抽出精度を維持することが出来る。
【００３８】
ステップＳ３０３において、ステップＳ３０２で周期分割された最初の断片を選択する。また、パラメータの初期値設定としてｌｏｏｐ＝１とする。
ステップＳ３０４において、選択された断片における速度データを近似する曲線ｆ２（ｘｉ）を決定する。
まず、周波数Ｆ（単位はＨｚ）の回転刺激に対して次式（６）で表される曲線（以下、「慣性曲線」と記す）ｆ２（ｘｉ）を定義する。次式（６）において、Ｐ、Ｑ、及びＲはフィッティングパラメータである。ωは、１秒当たり３０個のデータが得られる場合における、回転刺激の周波数Ｆに対する角周波数であり、上記式（４）で表される。
ｆ２（ｘｉ）＝Ｐ・ｓｉｎ（ω・ｘｉ）＋Ｑ・ｃｏｓ（ω・ｘｉ）＋Ｒ  ・・・・・・  （６）
【００３９】
選択された断片における速度データに対して、次式（７）によって示された、定義した慣性曲線ｆ２（ｘｉ）と各データｙｉとの二乗誤差の総和Ｓが最小となるように、慣性曲線ｆ２（ｘｉ）のパラメータＰ、Ｑ、及びＲが決定される。次式（７）において、ａ及びｂは周期分割された断片の始点及び終点、μｉはメンバーシップ値をそれぞれ示す。初期状態において、選択された断片における速度データの全ての点でμｉ＝１である。
【数２】
0004681377.tif 000003
【００４０】
ステップＳ３０５において、パラメータが決定された上記の慣性曲線ｆ２（ｘｉ）を用いてメンバーシップ値を決定する。
まず、速度データ（ｘｉ，ｙｉ）に対して、次式（８）を満たすｎｉの値を決定する。Δｙは、次式（９）で定義される解像度である。ｍａｘ（ｙｉ）は、データｙｉのうち、最大値であるデータを示す。ｍｉｎ（ｙｉ）は、データｙｉのうち、最小値であるデータを示す。
Δｙとして任意の値が設定されうるが、本実施例では、一例として、眼位データの最大値ｍａｘ（ｙｉ）及び最小値ｍｉｎ（ｙｉ）の差を２０で割った値とする。ｎｉは、速度データの各点における、慣性曲線ｆ２（ｘｉ）に対する、Δｙを単位とした速度データｙｉの分離度を示す値である。
ｎｉ・Δｙ≦（ｙｉ−ｆ２（ｘｉ））＜（ｎｉ＋１）・Δｙ  ・・・・・・  （８）
Δｙ＝（ｍａｘ（ｙｉ）−ｍｉｎ（ｙｉ））／２０  ・・・・・・  （９）
【００４１】
図７を参照して、メンバーシップ値の求め方の一例を説明する。
まず、図７に示すように、所定の速度データ個数（図７では速度データ３個分）に相当する幅ｗを有し、互いの重なり幅が（ｗ／３）である２つの窓Ｗ１及びＷ２を定義する。１つ目の窓Ｗ１の中にある速度データの各点（現在の注目点と以前の２つの点）の分離度の総和をｍ１ｉ（＝ｎｉ−２＋ｎｉ−１＋ｎｉ）、２つ目の窓Ｗ２の中にある眼振データの各点（現在の注目点と以後の２つの点）の分離度の総和をｍ２ｉ（＝ｎｉ＋ｎｉ＋１＋ｎｉ＋２）、１つ目の窓Ｗ１及び２つ目の窓Ｗ２の重なり部分の中にある眼振データの点（現在の注目点）の分離度の総和をｍ３ｉ（＝ｎｉ）とする。ｍ１ｉ、ｍ２ｉ、及びｍ３ｉは、上記式（８）から算出される。
【００４２】
算出されたｍ１ｉ、ｍ２ｉ、及びｍ３ｉを、次式（１０）に示すように全て掛け合わせる。これにより、１つ目の窓Ｗ１及び２つ目の窓Ｗ２の重なり部分に属する点（現在の注目点）における、重み付けを行ったメンバーシップ値μｉ’が求められる。
μｉ’＝ｍ１ｉ・ｍ２ｉ・ｍ３ｉ  （ｉ＝ａ＋２，ａ＋３，・・・，ｂ−３，ｂ−２）  ・・・（１０）
【００４３】
周期分割された断片において、２つの窓Ｗ１及びＷ２を順次ずらしていくことによって、速度データの各点におけるメンバーシップ値μｉ’が求められる。最終的に、求められたメンバーシップ値μｉ’に対して次式（１１）に示すように正規化を行い、速度データの各点におけるメンバーシップ値μｉを算出する。ｍａｘ（μｉ’）は、選択された断片において、上記のようにして求めたメンバーシップ値μｉ’のうち、最大値であるメンバーシップ値を示す。メンバーシップ値μｉは、１に近いほど慣性曲線ｆ２（ｘｉ）に近く、緩徐相である可能性が高いことを示す。
μｉ＝１−（μｉ’／ｍａｘ（μｉ’））  ・・・・・・  （１１）
【００４４】
ステップＳ３０６において、パラメータｌｏｏｐが３であるか否かを調べ、ｌｏｏｐ＝３であればステップＳ３０８に進む。そうでなければ、パラメータｌｏｏｐに１を足して（ステップＳ３０７）、ステップＳ３０４に戻り、メンバーシップ値μｉを求めるステップＳ３０４及びＳ３０５を繰り返す。ステップＳ３０４及びＳ３０５を複数回繰り返すことによって、最適なメンバーシップ値μｉが得られる。本実施例においては、ステップＳ３０４及びＳ３０５を３回繰り返しているが、繰り返す回数は、１回以上の任意の回数であって良い。
【００４５】
ステップＳ３０８において、現在選択されている断片が最後の断片であるか否かを調べる。最後の断片でなければ、次の断片を選択してパラメータｌｏｏｐを１に初期設定して（ステップＳ３０９）、ステップ３０４に戻り、次の断片に対して同様の処理を繰り返す。最後の断片であれば、ステップ３１０へ進む。
【００４６】
ステップＳ３１０において、緩徐相を抽出するための閾値を決定する。本実施例において、閾値決定法としてKittlerの方法を用いる。Kittlerの方法は、ある閾値によってメンバーシップ値の分布を２つのクラスに分割した場合に、あるメンバーシップ値がどちらのクラスに属するのかを示すあいまい度を、最悪の分布（正規分布）に対して最小にする方法である。
【００４７】
閾値ｋ（０＜ｋ＜１）に対して、閾値ｋ以上のメンバーシップ値を持つクラスをＣ０、閾値ｋ未満のメンバーシップ値を持つクラスをＣ１とする。以下の式（１２）によって示される、２クラスのどちらに属するのかのあいまい度Ｅ（ｋ）を最小にする値ｋを閾値として選定する。Nは全メンバーシップ値の数（即ち、全速度データの数）、Ｎ０はクラスＣ０に属するメンバーシップ値の数、Ｎ１はクラスＣ１に属するメンバーシップ値の数、σ０２（ｋ）及びσ１２（ｋ）はクラスＣ０及びＣ１におけるそれぞれのメンバーシップ値の分散を示す。また、ω０＝Ｎ０／Ｎ、ω１＝Ｎ１／Ｎの各式が成り立つ。
Ｅ（ｋ）＝ω０ｌｏｇ｛σ０／ω０｝＋ω１ｌｏｇ｛σ１／ω１｝  ・・・・・・  （１２）
【００４８】
なお、Kittlerの方法に代えて、既知の大津の方法等、他の閾値決定方法を用いても良い。しかし、Kittlerの方法が、緩徐相抽出精度が高いため好ましい。
【００４９】
ステップＳ３１１において、緩徐相の抽出が実行される。ステップＳ３１０で決定された閾値ｋと、各速度データの点におけるメンバーシップ値μｉとを比較し、メンバーシップ値μｉが閾値ｋ以上である場合に、そのデータ点が緩徐相であると見なす。緩徐相として抽出されたデータ点は、元データのグラフ上に黒点として印される等によって元データ上に反映され、表示画像データが作成される。
以上ステップＳ３０１〜Ｓ３１１の処理が終了した後、図２のステップＳ２０７に進む。図２のステップＳ２０７については、既に説明した。
【００５０】
次に、図８を参照して、ステップＳ２０５の回転刺激が与えられている場合の眼位データの緩徐相抽出方法について、詳細に説明する。図８は、回転刺激が与えられている場合の眼位データについての緩徐相抽出方法を示すフローチャートである。
【００５１】
回転刺激が与えられている場合の眼位データは、速度データに比べて正弦波形が保たれていないため、図３に示したフローチャートのステップに、差分データを算出するステップ（ステップＳ８０１）を追加する。
【００５２】
ステップＳ８０１において、眼位データの各点に対して、その眼位データと１つ前の眼位データとの差を取り、差分データを得る。回転刺激が与えられている場合の眼位データは、変位が激しいため、差分データを得ることで、眼位データをより正弦波状に近づけることができる。眼位データ（ｘｉ，ｙｉ）に対して、差分データｄｉは以下の式（１）で定義される。Ｎは、全眼位データ数である。
ｄｉ＝ｙｉ−ｙｉ−１  （ｉ＝２，３，・・・，Ｎ）  ・・・・・・  （１３）
【００５３】
ステップＳ８０２〜Ｓ８１２は、速度データに代えて上記のステップＳ８０１で得られた差分データに対して実行される点以外、上記した図３のステップＳ３０１〜Ｓ３１１と同様であるので、説明を省略する。以上ステップＳ８０１〜Ｓ８１２の処理が終了した後、図２のステップＳ２０７に進む。図２のステップＳ２０７については、既に説明した。ステップＳ８０１は、ステップＳ８０２の直前に実行されているが、ステップＳ８０２よりも後のタイミングで実行されても良い。
【００５４】
なお、図３及び図８に示したフローチャートにおいて、抽出された緩徐相に対する最終フィッティング曲線を求めるステップを追加しても良い。最終フィッティング曲線の振幅と位相が医学的に有用だからである。その場合、抽出された緩徐相に対して次式（１４）の曲線を定義し、この式のパラメータａ、ｂ、ｃ、及びｄを次式（１５）のＧを最小にするように決定することで最終フィッティング曲線を求める。なお、μｉは最終的に得られたメンバーシップ値、ｚｉは抽出された緩徐相成分である。
ｆ（ｘｉ）＝ａ＋ｂ・ｘｉ＋ｃ・ｓｉｎ（ω・ｘｉ＋ｄ）  ・・・・・・  （１４）
【数３】
0004681377.tif 000004
【００５５】
次に、図９を参照して、ステップＳ２０６の回転刺激が与えられていない場合の眼振データの緩徐相抽出方法について、詳細に説明する。図９は、回転刺激が与えられていない場合の眼振データについての緩徐相抽出方法を示すフローチャートである。
【００５６】
回転刺激が与えられていない場合の眼振データについては、正弦波形が得られないため、緩徐相の形状を予測することが困難で、回転刺激が与えられている場合のように慣性曲線を定義できない。そこで本実施例では、眼振データの各点に対して、その眼振データと１つ前の眼振データとの差の絶対値に着目し、慣性曲線の代わりに差分の絶対値を用いる。「差分の絶対値が小さい」ということは、すなわちその眼振データが緩徐相に属する度合いが高く、「差分の絶対値が大きい」ということは、すなわちその眼振データが急速相に属する度合いが高い。
【００５７】
なお、この手法は、眼位データ及び速度データのいずれにも適用可能である。以下、眼位データに適用した場合についてのみ記載し、速度データについての説明は同様であるため省略する。
【００５８】
ステップＳ９０１において、眼位データの各点に対して、その眼位データ及び１つ前の眼位データの差を算出し、差分データを得る。眼位データ（ｘｉ，ｙｉ）に対して、差分ｄｉは以下の式（１６）で定義される。Ｎは全眼位データ数である。
ｄｉ＝ｙｉ−ｙｉ−１  （ｉ＝２，３，・・・，Ｎ）  ・・・・・・  （１６）
【００５９】
ステップＳ９０２において、ステップＳ９０１で算出された差分データｄｉの絶対値を用いて、メンバーシップ値を決定する。
まず、差分データの各点において、次式（１７）のｎｉの値を決定する。Δｙは、次式（１８）で定義される解像度である。ｍａｘ（ｄｉ）は、上記のようにして求めた差分データｄｉのうち、最大値である差分データを示す。ｍｉｎ（ｄｉ）は、上記のようにして求めた差分データｄｉのうち、最小値である差分データを示す。
Δｙは、任意の値を設定しうるが、本実施例では、一例として、差分データの最大値ｍａｘ（ｄｉ）及び最小値ｍｉｎ（ｄｉ）の差を２０で割った値とする。ｎｉは、Δｙを単位とした差分データｄｉの絶対値の分離度を示す値である。
ｎｉ・Δｙ≦｜ｄｉ｜＜（ｎｉ＋１）・Δｙ  ・・・・・・  （１７）
Δｙ＝（ｍａｘ（ｄｉ）−ｍｉｎ（ｄｉ））／２０  ・・・・・・  （１８）
【００６０】
分離度ｎｉの分布の求め方の一例を説明する。
まず、所定の差分データ個数（例えば３個分）に相当する幅ｗを有し、互いの重なり幅が（ｗ／３）である２つの窓Ｗ１及びＷ２を定義する。１つ目の窓Ｗ１の中にある差分データの絶対値の各点（現在の注目点と以前の２つの点）の分離度の総和をｍ１ｉ（＝ｎｉ−２＋ｎｉ−１＋ｎｉ）、２つ目の窓Ｗ２の中にある差分データの絶対値の各点（現在の注目点と以後の２つの点）の分離度の総和をｍ２ｉ（＝ｎｉ＋ｎｉ＋１＋ｎｉ＋２）、１つ目の窓Ｗ１及び２つ目の窓Ｗ２の重なり部分の中にある差分データの絶対値の点（現在の注目点）の分離度の総和をｍ３ｉ（＝ｎｉ）とする。ｍ１ｉ、ｍ２ｉ、及びｍ３ｉは、上記式（１７）から算出される。
【００６１】
算出されたｍ１ｉ、ｍ２ｉ、及びｍ３ｉを、次式（１９）に示すように全て掛け合わせる。これにより、１つ目の窓Ｗ１及び２つ目の窓Ｗ２の重なり部分に属する差分データの絶対値の点（現在の注目点）における、重み付けを行ったメンバーシップ値μｉ’が求められる。
μｉ’＝ｍ１ｉ・ｍ２ｉ・ｍ３ｉ  （ｉ＝３，４，・・・，Ｎ−１，Ｎ−２）・・・・・・（１９）
【００６２】
２つの窓Ｗ１及びＷ２を順次ずらしていくことによって、差分データの絶対値の各点におけるメンバーシップ値μｉ’が求められる。最終的に、求められたメンバーシップ値μｉ’に対して次式（２０）に示されるように正規化を行い、差分データの絶対値の各点におけるメンバーシップ値μｉを算出さする。ｍａｘ（μｉ’）は、上記のようにして求めたメンバーシップ値μｉ’のうち、最大値であるメンバーシップ値を示す。メンバーシップ値μｉは、１に近いほど緩徐相である可能性が高いことを示す。
μｉ＝１−（μｉ’／ｍａｘ（μｉ’））  ・・・・・・  （２０）
【００６３】
ステップＳ９０３及びＳ９０４は、上記した図３のステップＳ３１０及びＳ３１１、上記した図８のステップＳ８１１及びＳ８１２と同様であるので、説明を省略する。以上ステップＳ９０１〜Ｓ９０４の処理が終了した後、図２のステップＳ２０７に進む。図２のステップＳ２０７については、既に説明した。
【００６４】
以上のように、本実施例における緩徐相抽出方法によれば、回転刺激が与えられていない場合の眼振データに対しても緩徐相を抽出することができる。なお、図９に示した回転刺激が与えられていない場合の眼振データについての緩徐相抽出方法は、回転刺激が与えられている場合の眼振データについても適用することができる。その場合、図９のフローチャートによる処理は、図３及び図８のフローチャートによる処理との比較からも明らかなように、処理時間が短くて済むという効果が得られる。
【００６５】
図１０〜図１８を参照して、本実施例の緩徐相抽出方法の精度について評価する。
まず、図１０を用いて、本実施例の緩徐相抽出方法をデータモデルに対して適用した場合の精度について述べる。図１０は、４つのデータモデルを示す図である。
【００６６】
図１０（ａ）は、回転刺激が与えられている場合のデータモデルであって、緩徐相成分の振幅が一定かつ水平方向に変位がない速度データモデルを示す。図１０（ｂ）は、回転刺激が与えられている場合のデータモデルであって、緩徐相成分の振幅が一定であるが、斜めに傾いている速度データモデルを示す。図１０（ｃ）は、回転刺激が与えられている場合のデータモデルであって、緩徐相成分の振幅が一定ではない速度データモデルを示す。図１０（ｄ）は、眼振データの形状がのこぎり状である眼位データモデルを示す。
図１０（ａ）、（ｂ）、（ｃ）、及び（ｄ）に示したデータモデルを、それぞれ以下モデルＡ、Ｂ、Ｃ、及びＤと記す。
【００６７】
モデルＡ、Ｂ、及びＣについては、緩徐相成分を正弦波とし、その正弦波信号に急速相成分に相当するランダムノイズを合成して作成した。モデルＤについては、緩徐相としての傾きがゆるやかな曲線及び急速相としての傾きが急な曲線をつなぎあわせて作成した。
【００６８】
上記の緩徐相抽出方法をモデルＡ、Ｂ、Ｃ、及びＤへ適用して抽出された緩徐相データと、上記のモデル化した眼振データ中の既知の緩徐相データを比較照合することにより、緩徐相抽出精度を定量化する。緩徐相抽出精度ｐを次式（２１）で定義する（ｐの単位は％である）。なお、抽出された緩徐相の数をＵ、全データ数をＮ、人為的に作り出した急速相の数をＶとする。
ｐ＝（Ｕ／（Ｎ−Ｖ））×１００  ・・・・・・  （２１）
【００６９】
図３のフローチャート及び図８のフローチャートを用いる緩徐相抽出方法をアルゴリズム１、図９のフローチャートを用いる緩徐相抽出方法をアルゴリズム２とする。
【００７０】
アルゴリズム１をモデルＡ、Ｂ、及びＣに対して適用して緩徐相抽出精度ｐを計算した結果、及び、アルゴリズム２をモデルＡ、Ｂ、Ｃ、及びＤに対して適用して緩徐相抽出精度ｐを計算した結果を表１に示す。アルゴリズム１をモデルＤに適用しない理由は、アルゴリズム１が緩徐相が正弦波状になるという仮定の下に解析を行うもので、モデルＤに対しては適用不可能であるためである。
表１から明らかなように、いずれも８５％以上の緩徐相抽出精度が得られた。モデルＤに対しては、９０％以上の緩徐相抽出精度が得られた。
【表１】
0004681377.tif 000005
【００７１】
次に、図１１〜図１８を参照して、本実施例の緩徐相抽出方法を実データに適用した場合の精度にについて述べる。
まず、図１１及び図１２を参照して、本実施例の緩徐相抽出方法を回転刺激が与えられている場合の眼振データへ適用した結果について説明する。図１１及び図１２は、０．１Ｈｚ及び０．５Ｈｚの回転刺激が与えられている場合の実測眼位データ及びそれに対する緩徐相抽出結果をそれぞれ示す図である。図１１及び図１２において、（ａ）は元の眼位データ、（ｂ）は本実施例の緩徐相抽出方法を用いて緩徐相を抽出した結果、（ｃ）は（ｂ）の一部拡大図をそれぞれ示している。
図１１及び図１２の（ｂ）及び（ｃ）に示されるように、０．１Ｈｚ及び０．５Ｈｚの両方のデータにおいて、眼位データ中の緩やかな変位部分が抽出されており、良好な緩徐相抽出結果が得られた。
【００７２】
図１３及び図１４を参照して、本実施例の緩徐相抽出方法を回転刺激が与えられている場合の速度データへ適用した結果について説明する。図１３及び図１４は、０．１Ｈｚ及び０．５Ｈｚの回転刺激が与えられている場合の実測速度データ及びそれに対する緩徐相抽出結果をそれぞれ示す図である。図１３及び図１４において、（ａ）は元の速度データ、（ｂ）は本実施例の緩徐相抽出方法を用いて緩徐相を抽出した結果、（ｃ）は（ｂ）の一部拡大図をそれぞれ示している。速度データに対しては、抽出された緩徐相に対して最終フィッティング曲線を求める処理を行った。
図１３及び図１４の（ｂ）及び（ｃ）に示されるように、０．１Ｈｚ及び０．５Ｈｚの両方のデータにおいて、速度データ中の緩やかな変位部分が抽出されており、良好な緩徐相抽出結果が得られた。
【００７３】
図１５、図１６、図１７、及び図１８を参照して、本実施例の緩徐相抽出方法を刺激が与えられていない場合の眼振データへ適用した結果について説明する。図１５、図１６、及び図１７は、後半規管型良性発作性頭位めまい症の症例（症例ｋｔ）、後半規管型良性発作性頭位めまい症のうち特に回旋成分の強い症例（症例ｇｔ）、及び外側半規管型良性発作性頭位めまい症の症例（症例ｒｉｅｄ）の実測眼位データ及びそれに対する緩徐相抽出結果を示す図である。図１８は、上記症例ｒｉｅｄの実測速度データ及びそれに対する緩徐相抽出結果を示す図である。図１５、図１６、図１７及び図１８において、（ａ）は元データ、（ｂ）は本実施例の緩徐相抽出方法を用いた緩徐相抽出の結果、（ｃ）は（ｂ）の一部拡大図をそれぞれ示している。
【００７４】
図１５、図１６、図１７、及び図１８に示されるように、いずれの症例のデータにおいても、一部で急速相が取り除けていない、あるいは、緩徐相が抽出されていない部分はあるが、全体として良好な結果が得られた。特に、症例ｒｉｅｄにおいては、緩徐相成分が指数的に減衰するという特徴が知られているが、このことと図１８に示された緩徐相抽出結果が合致しており、良好な結果が示された。
【００７５】
以上のように、本発明の緩徐相抽出方法は、刺激が与えられている場合の眼振データのみならず、刺激が与えられていない場合の眼振データからも高精度に緩徐相の自動抽出を行う、信頼性の高い緩徐相抽出方法及び眼振データ解析システムを実現することができる。
【産業上の利用可能性】
【００７６】
本発明にかかる緩徐相抽出方法及び眼振データ解析システムは、例えば、医療目的の眼振データ解析システム等に利用することができる。
【図面の簡単な説明】
【００７７】
【図１】眼振データにおける緩徐相及び急速相を説明するための図
【図２】本発明の実施例における、緩徐相抽出方法を示すフローチャート
【図３】本発明の実施例における、回転刺激が与えられている場合の速度データについての緩徐相抽出方法を示すフローチャート
【図４】正弦波形が上下に変位している状態の速度データを示す図
【図５】図４の速度データに対して処理を行い、パラメータを決定した場合の曲線を示す図
【図６】図５の曲線の周期に従って速度データを分割した様子を示す図
【図７】メンバーシップ値を求める方法を説明するための図
【図８】本発明の実施例における、回転刺激が与えられている場合の眼位データについての緩徐相抽出方法を示すフローチャート
【図９】本発明の実施例における、刺激が与えられていない場合の眼位データ及び速度データについての緩徐相抽出方法を示すフローチャート
【図１０】本発明の実施例における、眼振データの４つのデータモデルを示す図
【図１１】被験者に０．１Ｈｚの回転刺激が与えられている場合の実測眼位データ及びそれに対する緩徐相抽出結果を示す図
【図１２】被験者に０．５Ｈｚの回転刺激が与えられている場合の実測眼位データ及びそれに対する緩徐相抽出結果を示す図
【図１３】被験者に０．１Ｈｚの回転刺激が与えられている場合の実測速度データ及びそれに対する緩徐相抽出結果を示す図
【図１４】被験者に０．５Ｈｚの回転刺激が与えられている場合の実測速度データ及びそれに対する緩徐相抽出結果を示す図
【図１５】症例ｋｔの実測眼位データ及びそれに対する緩徐相抽出結果を示す図
【図１６】症例ｇｔの実測眼位データ及びそれに対する緩徐相抽出結果を示す図
【図１７】症例ｒｉｅｄの実測眼位データ及びそれに対する緩徐相抽出結果を示す図
【図１８】症例ｒｉｅｄの実測速度データ及びそれに対する緩徐相抽出結果を示す図
【図１９】本発明の実施例における緩徐相抽出方法を用いた眼振データ解析システムの構成を示すブロック図
【符号の説明】
【００７８】
１  緩徐相
２  急速相
【図１】
0004681377.tif 000006
【図２】
0004681377.tif 000007
【図３】
0004681377.tif 000008
【図４】
0004681377.tif 000009
【図５】
0004681377.tif 000010
【図６】
0004681377.tif 000011
【図７】
0004681377.tif 000012
【図８】
0004681377.tif 000013
【図９】
0004681377.tif 000014
【図１０】
0004681377.tif 000015
【図１１】
0004681377.tif 000016
【図１２】
0004681377.tif 000017
【図１３】
0004681377.tif 000018
【図１４】
0004681377.tif 000019
【図１５】
0004681377.tif 000020
【図１６】
0004681377.tif 000021
【図１７】
0004681377.tif 000022
【図１８】
0004681377.tif 000023
【図１９】
0004681377.tif 000024