(57)【要約】【課題】機械学習が不要で高精度な瞬き検出を可能とする、ドップラーセンサを用いた瞬き検出システムを提供すること。【解決手段】本瞬き検出システムは、ドップラー信号を得るドップラーセンサと、前記ドップラー信号の振幅信号を求め、前記振幅信号を時間軸方向に連続する複数のセルに分割し、注目するセルに隣接しない前後のセルの振幅値に基づいて第１閾値を決定し、前記注目するセルにおいて前記第１閾値以上の信号成分を検出する第１検出部と、瞬き波形の特徴に基づいて、前記第１閾値以上の信号成分の中から誤検出を除去し、残りの信号成分を瞬きと判定する瞬き判定部と、を有する。【選択図】図１
2017148447.tif 000002
【特許請求の範囲】
【請求項１】
ドップラー信号を得るドップラーセンサと、
前記ドップラー信号の振幅信号を求め、前記振幅信号を時間軸方向に連続する複数のセルに分割し、注目するセルに隣接しない前後のセルの振幅値に基づいて第１閾値を決定し、前記注目するセルにおいて前記第１閾値以上の信号成分を検出する第１検出部と、
瞬き波形の特徴に基づいて、前記第１閾値以上の信号成分の中から誤検出を除去し、残りの信号成分を瞬きと判定する瞬き判定部と、を有する瞬き検出システム。
【請求項２】
前記振幅信号を平滑化した振幅信号を生成し、前記平滑化した振幅信号と第２閾値とを比較し、前記第２閾値以下の信号成分を検出する第２検出部と、
前記第２閾値以下の信号成分のピーク検出を行うピーク検出部と、を有し、
前記瞬き判定部は、瞬き波形の特徴に基づいて、前記第１閾値以上の信号成分の中から誤検出を除去すると共に、前記ピーク検出部で検出されたピークの中から誤検出を除去し、それぞれの残りの信号成分を瞬きと判定する請求項１に記載の瞬き検出システム。
【請求項３】
前記瞬き判定部は、前記第１閾値以上の信号成分に含まれるピークの高さ及び前記ピーク検出部で検出されたピークの高さと基準値とをそれぞれ比較し、前記ピークの高さが前記基準値以下の部分を誤検出と判定するピーク高判定部を備えている請求項２に記載の瞬き検出システム。
【請求項４】
前記瞬き判定部は、前記ドップラー信号を構成するＩ信号及びＱ信号のそれぞれの電圧値の高低差を求め、前記高低差が共に第３閾値以上の場合に誤検出と判定する電圧高低差判定部を備えている請求項１乃至３の何れか一項に記載の瞬き検出システム。
【請求項５】
前記瞬き判定部は、前記ドップラー信号を構成するＩ信号及びＱ信号のそれぞれの電圧値の極大点と極小点の合計数を求め、前記合計数が共に第４閾値以上の場合に誤検出と判定する位相変化判定部を備えている請求項１乃至４の何れか一項に記載の瞬き検出システム。
【請求項６】
ドップラーセンサで被験者からの反射波を受信してドップラー信号を得る第１ステップと、
前記ドップラー信号の振幅信号を求め、前記振幅信号を時間軸方向に連続する複数のセルに分割し、注目するセルに隣接しない前後のセルの振幅値に基づいて第１閾値を決定し、前記注目するセルにおいて前記第１閾値以上の信号成分を検出する第２ステップと、
瞬き波形の特徴に基づいて、前記第１閾値以上の信号成分の中から誤検出を除去し、残りの信号成分を瞬きと判定する第３ステップと、を有する瞬き検出方法。
【請求項７】
前記振幅信号を平滑化した振幅信号を生成し、前記平滑化した振幅信号と第２閾値とを比較し、前記第２閾値以下の信号成分を検出する第４ステップと、
前記第２閾値以下の信号成分のピーク検出を行う第５ステップと、を有し、
前記第３ステップでは、瞬き波形の特徴に基づいて、前記第１閾値以上の信号成分の中から誤検出を除去すると共に、前記第５ステップで検出されたピークの中から誤検出を除去し、それぞれの残りの信号成分を瞬きと判定する請求項６に記載の瞬き検出方法。
【請求項８】
前記第３ステップは、前記第１閾値以上の信号成分に含まれるピークの高さ及び前記第５ステップで検出されたピークの高さと基準値とをそれぞれ比較し、前記ピークの高さが前記基準値以下の部分を誤検出と判定するステップを備えている請求項７に記載の瞬き検出方法。
【請求項９】
前記第３ステップは、前記ドップラー信号を構成するＩ信号及びＱ信号のそれぞれの電圧値の高低差を求め、前記高低差が共に第３閾値以上の場合に誤検出と判定するステップを備えている請求項６乃至８の何れか一項に記載の瞬き検出方法。
【請求項１０】
前記第３ステップは、前記ドップラー信号を構成するＩ信号及びＱ信号のそれぞれの電圧値の極大点と極小点の合計数を求め、前記合計数が共に第４閾値以上の場合に誤検出と判定するステップを備えている請求項６乃至９の何れか一項に記載の瞬き検出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、瞬き検出システム及び瞬き検出方法に関する。
【背景技術】
【０００２】
近年、コンピュータの画像表示端末（ＶＤＴ：Visual Display Terminal）での作業の普及に伴い、多くの作業者から肩や腰、手、特に目の痛みを訴える事例が多く挙げられており、これらの事例はＶＤＴ症候群等と呼ばれている。長時間のＶＤＴ作業は、眼精疲労等の様々な眼症状を引き起こすため、ＶＤＴ作業者の間で問題となっている。
【０００３】
瞬きは涙液膜を広げたり角膜を守ったりする役割があり、瞬き情報の取得は疲労度や集中度、眠気といった身体状態の把握のみならず、ＶＤＴ作業者の眼症状を防ぐという観点でも有効である。
【０００４】
こうした背景から、瞬き情報の取得に注目した様々な取り組みが行われており、例えば、機械学習に基づく手法が提案されている（例えば、非特許文献１参照）。この手法では、瞬き信号について５つの特徴量（電圧最大値、時間幅、分散、生信号の分散、生信号の電圧最大値）を定義し、機械学習の一つであるＳＶＭ（Support Vector Machine）により瞬きと非瞬きを識別する。
【０００５】
これら特徴量は人や環境毎に異なるため、事前に様々な瞬きについて特徴量を学習することで、検出特性が向上する。しかし、機械学習を用いた手法では、人や環境毎に事前学習を行う必要があり、学習データとテストデータが適さない場合には検出精度が劣化してしまう。
【先行技術文献】
【非特許文献】
【０００６】
【非特許文献１】C，Tamba and T．Ohtsuki，"Learning−based Blink Detection Using a Doppler Sensor，"IEICE technical report，vol．114，no．418，ASN2014−123，pp．97−102．Jan．2015．
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本発明は、上記の点に鑑みてなされたものであり、機械学習が不要で高精度な瞬き検出を可能とする、ドップラーセンサを用いた瞬き検出システムを提供することを目的とする。
【課題を解決するための手段】
【０００８】
本瞬き検出システムは、ドップラー信号を得るドップラーセンサと、前記ドップラー信号の振幅信号を求め、前記振幅信号を時間軸方向に連続する複数のセルに分割し、注目するセルに隣接しない前後のセルの振幅値に基づいて第１閾値を決定し、前記注目するセルにおいて前記第１閾値以上の信号成分を検出する第１検出部と、瞬き波形の特徴に基づいて、前記第１閾値以上の信号成分の中から誤検出を除去し、残りの信号成分を瞬きと判定する瞬き判定部と、を有することを要件とする。
【発明の効果】
【０００９】
開示の技術によれば、機械学習が不要で高精度な瞬き検出を可能とする、ドップラーセンサを用いた瞬き検出システムを提供できる。
【図面の簡単な説明】
【００１０】
【図１】本実施の形態に係る瞬き検出システムを例示する機能ブロック図である。
【図２】ドップラーセンサ１０で得たドップラー信号の一例である。
【図３】本実施の形態に係る瞬き検出方法のフローチャートの一例である。
【図４】ＣＡ−ＣＦＡＲ検出について説明する図である。
【図５】振幅信号Ａ、振幅信号Ａを平滑化した振幅信号ＡＭについて説明する図である。
【図６】第２閾値ＴＨ２以下の振幅信号ＡＭのピークの一例である。
【図７】ピーク高判定部の判定について説明する図である。
【図８】Ｉ信号及びＱ信号のそれぞれの高低差ＭＩ及びＭＱについて説明する図である。
【図９】平滑化したＩ信号及びＱ信号から得られた極大点と極小点について説明する図である。
【発明を実施するための形態】
【００１１】
以下、図面を参照して、実施の形態の説明を行う。なお、各図面において、同一構成部分には同一符号を付し、重複した説明を省略する場合がある。
【００１２】
図１は、本実施の形態に係る瞬き検出システムを例示する機能ブロック図である。図１に示す瞬き検出システム１は、主要な構成要素として、ドップラーセンサ１０と、信号処理部２０とを有している。
【００１３】
ドップラーセンサ１０は、ドップラー効果による送信信号と受信信号の周波数シフトを観測することで、観測対象（被験者）の動きを検出するセンサである。本実施の形態では、一例として、送信波として無変調連続波（ＣＷ：Continuous Wave）を用いる。
【００１４】
ドップラーセンサ１０は被験者の近傍に配置され、被験者の瞼又はその近傍で反射された信号（反射波）を受信し、受信信号に送信信号の周波数を掛け合わせてドップラー信号を得る。被験者としては、例えば、画像表示端末（ＶＤＴ）で作業する作業者、車両の運転者等が挙げられる。
【００１５】
図２は、ドップラーセンサ１０で得たドップラー信号の一例である。図２に示す信号は、送信信号と受信信号の間の周波数シフトを表わすドップラー信号を時間の関数として取得したものであり、送信信号と同相（In-phase）成分であるＩ信号、及び直交位相（Quadrature）成分であるＱ信号で構成される。
【００１６】
瞬き、心拍や呼吸、体動（身体による動作）等の身体の表面の動きはドップラーセンサ１０により観測できる。身体の表面積が広い体動の場合、ドップラーセンサ１０の出力信号も大きくなるため、５０Ｈｚ程度の高周波数領域までスペクトルが広がって分布する。
【００１７】
一方、心拍や呼吸等の小さな身体動作は、スペクトルが数Ｈｚ程度の低周波領域に分布する傾向にあり、心拍成分は０．８〜１．５Ｈｚ程度、呼吸成分は０．１５〜０．４Ｈｚ程度の周波数となる。瞬きは、心拍や呼吸よりも高く、体動よりも低い周波数領域となり、約１０〜２０Ｈｚまでスペクトルが広がる。
【００１８】
従って、瞬き検出を行う際には、バンドパスフィルタ等を用いて適宜ノイズ除去を行うことが好ましい。
【００１９】
図１に戻り、信号処理部２０は、ドップラーセンサ１０の出力信号であるドップラー信号に基づいて、被験者の瞬きを検出する。信号処理部２０は、適宜、ドップラーセンサ１０で受信したドップラー信号のＩ信号及びＱ信号をそのまま利用したり、Ｉ信号及びＱ信号に基づいて振幅信号を生成したりすることができる。
【００２０】
信号処理部２０は、ＣＡ−ＣＦＡＲ（Cell Averaging Constant False Alarm Rate）検出部２１と、体動検出部２２と、ピーク検出部２３と、瞬き判定部２４とを有している。又、瞬き判定部２４は、ピーク高判定部２４１と、電圧高低差判定部２４２と、位相変化判定部２４３とを有している。これら各部の機能については、図３のフローチャートを参照しながら後述する。
【００２１】
信号処理部２０は、例えば、ＣＰＵ（Central Processing Unit）、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）、メインメモリ等を含む構成とすることができる。この場合、信号処理部２０の各種機能は、ＲＯＭ等に記録されたプログラムがメインメモリに読み出されてＣＰＵにより実行されることによって実現できる。但し、信号処理部２０の一部又は全部は、ハードウェアのみにより実現されてもよい。又、信号処理部２０は、物理的に複数の装置等により構成されてもよい。
【００２２】
図３は、本実施の形態に係る瞬き検出方法のフローチャートの一例である。図３を中心にして適宜他の図も参照しながら、本実施の形態に係る瞬き検出方法について説明する。
【００２３】
まず、ステップＳ１０１では、ドップラーセンサ１０で被験者からの反射波を受信してドップラー信号を得る。ここで得るドップラー信号は、例えば、図２に示すようなＩ信号及びＱ信号で構成されている。
【００２４】
次に、ステップＳ１０２では、ステップＳ１０１で受信したドップラー信号のノイズ除去を行う。ノイズ除去は、ドップラーセンサ１０と信号処理部２０との間にバンドパスフィルタ等を挿入してハードウェア的に行ってもよいし、ドップラーセンサ１０の出力信号を直接信号処理部２０に入力し、信号処理部２０内においてディジタル信号処理（ディジタルフィルタ等）により行ってもよい。
【００２５】
又、除去するべきノイズの周波数（通過させるべき信号の周波数）は各処理により異なる場合があるから、通過帯域の異なる複数の信号を出力可能に構成しておくことが好ましい。例えば、通過帯域を５〜５０Ｈｚ程度とした信号、通過帯域を５〜２５Ｈｚ程度とした信号、ノイズ除去を行わない信号等を適宜出力させるように構成することができる。但し、ステップＳ１０２は必要に応じて実行すればよい。
【００２６】
次に、ステップＳ１０３では、第１検出部であるＣＡ−ＣＦＡＲ検出部２１は、ステップＳ１０２においてノイズ除去されたドップラー信号に基づいてＣＡ−ＣＦＡＲ検出を行う。まず、ノイズ除去されたドップラー信号（Ｉ信号及びＱ信号）から式（１）により振幅信号Ａを求める。なお、ＣＡ−ＣＦＡＲ検出では、ステップＳ１０２で通過帯域を５〜５０Ｈｚ程度とした信号を用いることが好ましい。前述のように、心拍成分は０．８〜１．５Ｈｚ程度、呼吸成分は０．１５〜０．４Ｈｚ程度であるから、通過帯域を５〜５０Ｈｚ程度とすることで、心拍や呼吸の影響を低減できるからである。
【００２７】
【数１】
2017148447.tif 000003
次に、振幅信号Ａを用いてＣＡ−ＣＦＡＲ検出を行う。ＣＡ−ＣＦＡＲ検出では、所望の誤警報確率に基づいて第１閾値ＴＨ１を決定する。第１閾値ＴＨ１を決定するには、図４に示すように振幅信号Ａを時間軸方向に連続する複数のセルに分割し、注目するセル（Ｃｕｔ：Cell under test）を決定する。
【００２８】
次に、Ｃｕｔに隣接しない前後のセルの振幅値に基づいて第１閾値ＴＨ１を決定する。具体的には、Ｃｕｔの直前の所定数のセルをガードセルＧＦ、Ｃｕｔの直後の所定数のセルをガードセルＧＲとし、更に、ガードセルＧＦの直前の所定数のセルをフロントセルＦ、ガードセルＧＲの直後の所定数のセルをリアセルＲとする。ドップラーセンサ１０で観測された瞬き信号の時間長は、経験的に約０．４秒であることがわかっている。そこで、図４における各セルのセル長は０．２秒程度に設定することが好ましい。
【００２９】
なお、図４では、ガードセルＧＦ、ガードセルＧＲ、フロントセルＦ、及びリアセルＲをそれぞれ３セルで構成しているが、これは一例であり、ガードセルＧＦ、ガードセルＧＲ、フロントセルＦ、及びリアセルＲを構成するセル数は適宜決定することができる。
【００３０】
次に、フロントセルＦの振幅値とリアセルＲの振幅値の平均値ＡＡＶＥを算出し、式（２）に示すように平均値ＡＡＶＥを係数α倍したものを第１閾値ＴＨ１とする。なお、ガードセルＧＦ及びＧＲは第１閾値ＴＨ１の算出には用いない。これは、注目するセルであるＣｕｔ自身が第１閾値ＴＨ１へ影響することを防ぐためである。
【００３１】
【数２】
2017148447.tif 000004
第１閾値ＴＨ１が低すぎると、多くの雑音を瞬きと誤検出しＦＰ（False Positive）が増加する。一方、第１閾値ＴＨ１が高すぎると、瞬きの検出漏れが多くなりＦＮ（False Negative）が増加する。そこで、係数αを適切な値に設定する必要があるが、発明者らの事前実験により、α＝３．４程度とすることでＦＰ及びＦＮを低減できることがわかっている。
【００３２】
以上のようにして第１閾値ＴＨ１が決定されると、ＣＡ−ＣＦＡＲ検出部２１は、Ｃｕｔにおいて第１閾値ＴＨ１以上の信号成分を検出する。引き続きＣｕｔの位置を移動させながら同様の検出を繰り返す。そして、第１閾値ＴＨ１以上の信号成分は、瞬きである可能性があると判断され、ステップＳ１０６（後述）の処理が実行される。
【００３３】
次に、ステップＳ１０４では、第２検出部である体動検出部２２は、体動検出を行う。体動検出は、身体の大きな動きを検出し、瞬き検出処理から除外するために行う処理である。例えば、退屈な直線道路の運転中やＶＤＴ作業集中時等、動作中でない状況での瞬きを検出対象とする。瞬きによる振幅値は小さくなるのに対し、身体の大きな動きによる振幅値は大きくなるため、大きな振幅値を持つ信号は体動によるものであると判別できる。
【００３４】
具体的には、体動検出部２２は、ステップＳ１０２においてノイズ除去されたドップラー信号に基づいて体動検出を行う。まず、ノイズ除去されたドップラー信号から式（１）により振幅信号Ａを求める。ここでは、ステップＳ１０２で通過帯域を５〜２５Ｈｚ程度とした信号（心拍や呼吸、高周波雑音成分を低減した信号）を用いて図５に示す振幅信号Ａを求める。
【００３５】
続いて、振幅信号Ａを移動平均により平滑化した振幅信号ＡＭを生成する。そして、振幅信号ＡＭと第２閾値ＴＨ２とを比較する。そして、第２閾値ＴＨ２より大きな信号成分を体動によるものと判定する。なお、第２閾値ＴＨ２は、事前実験に基づいて適宜決定することができる。
【００３６】
振幅信号ＡＭが第２閾値ＴＨ２以下である信号成分は、体動が検出されていないものと判断され、ステップＳ１０５へ進む。ステップＳ１０５では、ピーク検出部２３は、ステップＳ１０４における第２閾値ＴＨ２以下の振幅信号ＡＭのピーク検出を行う。図６は、ステップＳ１０５で検出された第２閾値ＴＨ２以下の振幅信号ＡＭのピークの一例である。図６に示すように、通常は複数のピークが検出される。
【００３７】
ピーク検出部２３において第２閾値ＴＨ２以下の振幅信号ＡＭでピークが検出された信号成分は、瞬きである可能性があると判断され、ステップＳ１０６（後述）の処理が実行される。
【００３８】
なお、ステップＳ１０３の処理とステップＳ１０４及びＳ１０５の処理は、ステップＳ１０１でドップラーセンサ１０で得た同一のドップラー信号に対して行われる処理であり、並行して行うことができる。
【００３９】
次に、ステップＳ１０６では、瞬き判定部２４は、ステップＳ１０３において第１閾値ＴＨ１以上の信号成分、及びステップＳ１０５においてピークが検出された信号成分について、瞬き波形の特徴に基づいて瞬き判定（誤検出の除去）を行う。
【００４０】
ステップＳ１０３のＣＡ−ＣＦＡＲ検出やステップＳ１０５のピーク検出でステップＳ１０６の処理対象であると判断された信号成分は、瞬きである可能性が高い。しかし、ステップＳ１０３やＳ１０５において心拍や呼吸成分を完全に取り除けず、ステップＳ１０６の処理対象であると判断された信号成分に、心拍や呼吸成分のピークが含まれている場合があるため、いくつかのＦＰが生じる。これらのＦＰを取り除くことがステップＳ１０６の処理の目的である。
【００４１】
まず、ピーク高判定部２４１は、ステップＳ１０３において第１閾値ＴＨ１以上の信号成分、及びステップＳ１０５においてピークが検出された信号成分について、ピーク高判定を行うことでＦＰを取り除く。具体的には、ステップＳ１０３において第１閾値ＴＨ１以上の信号成分に含まれるピークの高さ、及びステップＳ１０５においてピークが検出された信号成分のピークの高さと基準値Ｓｖとをそれぞれ比較し、ピークの高さが基準値Ｓｖ以下の部分を誤検出（非瞬き）と判定する。
【００４２】
例えば、図７に示すように、ピーク高判定部２４１は、ピークの高さが基準値Ｓｖ以下の部分（丸で囲んだ部分）を誤検出（非瞬き）であると判定する。なお、図７は、ステップＳ１０５で検出された第２閾値ＴＨ２以下の振幅信号ＡＭのピークの高さと基準値Ｓｖとを比較した例であるが、ステップＳ１０３において第１閾値ＴＨ１以上の信号成分に含まれるピークの高さについても同様に基準値Ｓｖと比較する。
【００４３】
ここで基準値Ｓｖは、式（３）に示すように、ピークの高さの平均値ＰＡＶＥの定数ｃ倍で定義される。式（３）における定数ｃは経験的に０．７程度に設定することが好ましい。ピークの高さは、瞬きの様子やドップラーセンサ１０からの距離によって変化するが、長時間観測の場合は基準値Ｓｖの値を更新することで対応することができる。
【００４４】
【数３】
2017148447.tif 000005
次に、電圧高低差判定部２４２は、Ｉ信号の電圧値及びＱ信号の電圧値の高低差に基づいてＦＰを取り除く。まず、ハイパスフィルタにより、ステップＳ１０１で得たドップラー信号（Ｉ信号及びＱ信号）から心拍や呼吸の周波数成分を除去する。ここで、ハイパスフィルタのカットオフ周波数は、例えば、５Ｈｚ程度に設定することができる。
【００４５】
そして、図８に示すように、Ｉ信号及びＱ信号のそれぞれにおいて、０．２秒の時間幅内における電圧値の高低差ＭＩ及びＭＱを求める。高低差ＭＩ及びＭＱが共に式（４）により求まる第３閾値ＴＨ３以上の場合、表面積の大きな動きによってドップラーセンサ１０の出力信号（Ｉ信号及びＱ信号）が大きくなったと考えられ、誤検出（体動による信号成分）と判定できる。
【００４６】
【数４】
2017148447.tif 000006
次に、位相変化判定部２４３は、Ｉ信号の電圧値及びＱ信号の電圧値の極大点と極小点の合計数に基づいてＦＰを取り除く。まず、ローパスフィルタにより、ステップＳ１０１で得たドップラー信号（Ｉ信号及びＱ信号）の高周波数成分を除去して平滑化する。ここで、ローパスフィルタのカットオフ周波数は、例えば、２５Ｈｚ程度に設定することができる。
【００４７】
体動による信号波形は極大点と極小点を多く持つ傾向にある。そこで、平滑化したＩ信号及びＱ信号から得られた極大点と極小点の合計数が共に第４閾値ＴＨ４以上の場合は、その信号が誤検出（体動による信号）であると判定する。ここでは、一例として第４閾値ＴＨ４を５点とするが、これには限定されない。
【００４８】
第４閾値ＴＨ４を５点とした場合、例えば、図９（ａ）では極小点が１点であるから瞬きであると判定し、図９（ｂ）では極大点と極小点の合計数が５点であるから体動であると判定する。
【００４９】
以上をまとめると、ステップＳ１０６において、ピーク高判定部２４１の判定でピークの高さが基準値Ｓｖ以下の信号、電圧高低差判定部２４２の判定で高低差ＭＩ及びＭＱが共に第３閾値ＴＨ３以上の信号、及び位相変化判定部２４３の判定で平滑化したＩ信号及びＱ信号から得られた極大点と極小点の合計数が共に第４閾値ＴＨ４以上の信号がＦＰとして除外される。
【００５０】
そして、ステップＳ１０３において第１閾値ＴＨ１以上の信号成分、及びステップＳ１０５においてピークが検出された信号成分のうち、ステップＳ１０６でＦＰとして除外された信号成分を除いた残りの信号成分が瞬きと判定される。ここで、ピーク高判定部２４１の判定、電圧高低差判定部２４２の判定、及び位相変化判定部２４３の判定の処理順
は任意に決定して構わない。
【００５１】
なお、ステップＳ１０１〜Ｓ１０６の全てを実行する必要はなく、最低限ステップＳ１０１、Ｓ１０３、及びＳ１０６（ピーク高判定、電圧高低差判定、位相変化判定の何れか１つ）を実行することで、従来よりも瞬き検出の精度を向上できる。
【００５２】
但し、ステップＳ１０４及びＳ１０５を実行することで瞬き検出の精度を向上できる。又、ステップＳ１０６でピーク高判定、電圧高低差判定、位相変化判定の何れか２つを実行することで更に瞬き検出の精度を向上でき、３つ全てを実行することで更に瞬き検出の精度を向上できる。
【００５３】
このように、本実施の形態に係る瞬き検出システム１では、ＣＡ−ＣＦＡＲ検出部２１によるＣＡ−ＣＦＡＲ検出において第１閾値ＴＨ１を動的に決定することで、身体の揺らぎや雑音に対するロバスト性を高め、瞬きの検出漏れ（ＦＮ）を低減することができる。又、瞬き波形の特徴に基づく瞬き判定部２４の処理により、誤って瞬きと検出しやすい微かな体動を除外し、誤検出（ＦＰ）の数を低減させることができる。そして、これらを組み合わせることで、機械学習が不要で高精度な瞬き検出が可能となる。
【００５４】
［実施例］
ドップラーセンサ１０を用いて３種類の実験を行った。実験の諸元を表１に示す。
【００５５】
【表１】
2017148447.tif 000007
実験１は異なる環境（１）及び（２）、実験２は異なる被験者Ａ〜Ｄでの実験である。実験３は、ＣＡ−ＣＦＡＲ検出で生じるＦＰを低減する瞬き判定部２４による処理の効果を確かめる実験である。実験３において、状況１では被験者は目を瞑りながら着座静止し、身体の揺らぎの影響がある状況とする。状況２では、身体を前後に動かす動作と静止を各２〜３秒間ずつ交互に繰り返す。
【００５６】
又、評価には検出率Ｒｒと適合率Ｐｒを用いる。検出率Ｒｒは実際の瞬きのうち検出できた瞬きの割合を、適合率Ｐｒは瞬きと検出したうちの正解の割合を示し、式（５）及び式（６）で定義される。式（５）及び式（６）において、ＴＰは正しく検出できた瞬きの数、ＦＮは検出漏れした瞬きの数、ＦＰは誤って瞬きとして検出した数を示している。
【００５７】
【数５】
2017148447.tif 000008
【００５８】
【数６】
2017148447.tif 000009
実験１の結果、及び比較として従来法の結果を表２に示す。なお、ここでは、非特許文献１に記載の機械学習を用いた手法を従来法、図３等を参照して説明した本実施の形態に係る瞬き検出方法を提案法と称する。
【００５９】
【表２】
2017148447.tif 000010
従来法ではＳＶＭを用いて様々な瞬き波形を事前に学習することで瞬き検出性能を向上させ、９９．２％と高精度な検出率Ｒｒを達成している。しかし適合率Ｐｒは９０．２％であり、十分な値とはいえない。これは従来法が微かな体動や雑音を誤って瞬きとして検出し、ＦＰが増加したためである。
【００６０】
一方、提案法は検出率Ｒｒ、適合率Ｐｒともに約９９％の高い瞬き検出精度を達成した。提案法では、瞬き信号の検出にＣＡ−ＣＦＡＲ検出を適用することで身体の揺らぎや雑音に対するロバスト性を高めたことでＦＰを低減した。更に、瞬き波形の特徴に基づく瞬き判定部２４の処理により、瞬きとして誤検出しやすい微かな体動を除外し、ＦＰを低減した。このように、提案法が従来法と比較して瞬き検出の精度を向上できることが確認された。
【００６１】
次に、実験２の結果を表３に示す。
【００６２】
【表３】
2017148447.tif 000011
被験者Ｄにおいて適合率Ｐｒが９７．４％と被験者Ａ、Ｂ、及びＤに比べて僅かに低い結果となった。これは、判定過程における各パラメータ（定数ｃ、第１閾値ＴＨ１、第２閾値ＴＨ２、第３閾値ＴＨ３、第４閾値ＴＨ４等）が被験者Ａに基づいて設定されていることによりＦＰ低減の効果が低下したためと考えられる。しかし、この結果から、提案法が平均して９９％以上の高い瞬き検出の精度を達成できることが確認された。
【００６３】
なお、この結果から、判定過程における各パラメータを複数の被験者に基づいて設定することで、適合率Ｐｒが更に向上すると考えられる。又、瞬き判定の前に、被験者毎にキャリブレーションを行い、被験者毎に各パラメータを微調整することで、適合率Ｐｒが更に向上すると考えられる。このキャリブレーションは、従来法のような大掛かりな学習とは異なり、極短時間で各パラメータを微調整できる。
【００６４】
最後に、実験３の結果を表４に示す。
【００６５】
【表４】
2017148447.tif 000012
身体の揺らぎを伴った状況１において、ＣＡ−ＣＦＡＲ検出のみを適用した場合には，身体の揺らぎによる誤検出ＦＰが生じた。これは、身体の揺らぎによる信号がＣＡ−ＣＦＡＲ検出の第１閾値ＴＨ１を超えてしまうためである。又、身体の大きな動作を伴った状況２において、その動作も多くのＦＰの原因となった。しかし、何れの場合にも、瞬き波形の特徴に基づく瞬き判定部２４の処理により、瞬きとして誤検出しやすい体動を除外し、ＦＰを低減できることが確認された。
【００６６】
以上の実験結果から、提案法のようにＣＡ−ＣＦＡＲ検出部２１によるＣＡ−ＣＦＡＲ検出と瞬き判定部２４によるＦＰ除外処理とを組み合わせることで、ＴＰの数を増加させつつＦＰの数を低減させることが可能となり、検出率Ｒｒ及び適合率Ｐｒが何れも約９９％である高精度の瞬き検出を実現できることがわかった。
【００６７】
以上、好ましい実施の形態等について詳説したが、上述した実施の形態等に制限されることはなく、特許請求の範囲に記載された範囲を逸脱することなく、上述した実施の形態等に種々の変形及び置換を加えることができる。
【符号の説明】
【００６８】
１  瞬き検出システム
１０  ドップラーセンサ
２０  信号処理部
２１  ＣＡ−ＣＦＡＲ検出部
２２  体動検出部
２３  ピーク検出部
２４  瞬き判定部
２４１  ピーク高判定部
２４２  電圧高低差判定部
２４３  位相変化判定部
【図１】
2017148447.tif 000013
【図２】
2017148447.tif 000014
【図３】
2017148447.tif 000015
【図４】
2017148447.tif 000016
【図５】
2017148447.tif 000017
【図６】
2017148447.tif 000018
【図７】
2017148447.tif 000019
【図８】
2017148447.tif 000020
【図９】
2017148447.tif 000021