(57)【要約】【課題】個人毎に正確な覚醒度を判定すること。【解決手段】覚醒度判定装置１００は、心拍検出部１０１が被験者から心拍信号を検出し、心拍間隔算出部１０２が心拍間隔を算出し、スペクトル算出部１０３が、周波数毎のスペクトル密度を算出する。そして、ピーク周波数算出部１０４、ピーク周波数クラスタリング部１０５により、特徴点を特定し、グラフ調整部１０７が、特徴点および被験者の行動に関する情報を基にしてグラフを眠気レベル毎に複数の領域に分割し、システム応答判定部１０８が、グラフ上の特徴点に基づいて被験者が実際に眠くなっているか否かを判定する。【選択図】図２
2010155072.tif 000002
【特許請求の範囲】
【請求項１】
被験者の心拍信号から振幅ピークを検出し、検出した各振幅ピークの間隔を心拍間隔として検出する心拍間隔検出手段と、
時間経過に伴う前記心拍間隔の変動に基づいて当該変動の周波数毎にスペクトル密度を算出し、算出したスペクトル密度が極大となる特徴点を判定する特徴点判定手段と、
前記被験者の行動に関する情報を含んだ被験者情報と前記特徴点とを基にして前記被験者の覚醒度を判定する覚醒度判定手段と
を有する覚醒度判定装置。
【請求項２】
前記覚醒度判定手段は、前記周波数とスペクトル密度とを軸とするグラフを作成し、当該グラフ上の前記特徴点の位置と前記被験者情報とを基にして、前記グラフ上の領域を前記覚醒度に応じて分割するグラフ作成手段と、前記グラフ上に前記特徴点をプロットし、当該特徴点の位置する領域に基づいて前記被験者の覚醒度を判定する判定手段とを有する請求項１に記載の覚醒度判定装置。
【請求項３】
前記グラフ作成手段は、前記特徴点の移動範囲に基づいて前記グラフのスケールを調整する請求項２に記載の覚醒度判定装置。
【請求項４】
前記グラフ上における前記特徴点の移動方向を基にして前記被験者に警告を出力するか否かを判定する警告判定手段と、前記警告判定手段の判定結果に基づいて前記被験者にシステム応答を実行する手段とを更に有する請求項２または３に記載の覚醒度判定装置。
【請求項５】
覚醒度判定装置が、
被験者の心拍信号から振幅ピークを検出し、検出した各振幅ピークの間隔を心拍間隔として検出して記憶装置に記憶する心拍間隔検出ステップと、
時間経過に伴う前記心拍間隔の変動に基づいて当該変動の周波数毎にスペクトル密度を算出し、算出したスペクトル密度が極大となる特徴点を判定する特徴点判定ステップと、
前記被験者の行動に関する情報を含んだ被験者情報と前記特徴点とを基にして前記被験者の覚醒度を判定する覚醒度判定ステップと
を含んだ覚醒度判定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、被験者の覚醒度を判定する覚醒度判定装置および覚醒度判定方法に関するものである。
【背景技術】
【０００２】
近年、被験者の眠気や覚醒度を計測する様々な技術が考案されている。例えば、運転手の撮影画像の顔面領域から、運転手の上瞼、及び下瞼を含む目の周辺領域を抽出し、眼球と瞼との輝度差に基づき、上瞼の最高点と下瞼の最下点との距離（以下、各瞼の距離と表記する）を算出することで、覚醒度を求める技術が存在する（例えば、特許文献１参照）。
【０００３】
また、被験者の脈拍信号を取得し、脈拍間隔の揺らぎ周波数の変化に基づいて、覚醒度を判定する技術も知られている。かかる技術では、脈拍間隔の揺らぎ周波数の分布を解析して低周波成分に比較して高周波成分が優位になる場合に、覚醒度が低くなる（被験者が眠い状態である）と判定している。いわば、高周波成分と低周波成分の比率の変化で判断していた。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００８−１８６２６３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、上述した従来の技術では、個人毎に正確な覚醒度を判定することができないという問題があった。
【０００６】
具体的には、被験者毎の各瞼の距離には個人差があるにもかかわらず、上述した手法では、被験者から算出した各瞼の距離によって求められる覚醒度を、例えば、統計的手法で求められた各瞼の距離と覚醒度との関係から特定しているため、特定した覚醒度は、必ずしも各被験者の覚醒度に合致していない。また、お化粧などの影響で下瞼と上瞼の位置を正確にとらえることができずに、距離計算が誤ったものになるなどの影響がでていた。
【０００７】
また、脈拍間隔の揺らぎ周波数の変化に基づいて、覚醒度を判定する場合であっても、覚醒度が下がるにつれてゆらぎ周波数の比率がどの程度変化するかは個人差があり、個人毎に精度よく覚醒度を算出するためには、検出した揺らぎ周波数の変化の度合が、その個人にとってどの程度の覚醒度の変化度合に相当するのかを求める必要がある。
【０００８】
この発明は、上述した従来技術による問題点を解消するためになされたものであり、個人毎に正確な覚醒度を判定することができる覚醒度算出装置および覚醒度判定方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
上述した課題を解決し、目的を達成するため、この覚醒度判定装置は、被験者の心拍信号から振幅ピークを検出し、検出した各振幅ピークの間隔を心拍間隔として検出する心拍間隔検出手段と、時間経過に伴う前記心拍間隔の変動に基づいて当該変動の周波数毎にスペクトル密度を算出し、算出したスペクトル密度が極大となる特徴点を判定する特徴点判定手段と、前記被験者の行動に関する情報を含んだ被験者情報と前記特徴点とを基にして前記被験者の覚醒度を判定する覚醒度判定手段とを有することを要件とする。
【発明の効果】
【００１０】
この覚醒度判定装置によれば、個人毎に正確な覚醒度を判定することができる。
【図面の簡単な説明】
【００１１】
【図１】図１は、本実施例１にかかる覚醒度判定装置の概要および特徴を説明するための図である。
【図２】図２は、本実施例１にかかる覚醒度判定装置の構成を示す機能ブロック図である。
【図３】図３は、心拍検出部が検出する心拍信号の一例を示す図である。
【図４】図４は、心拍間隔変動データを心拍間隔−時間平面上で示した図である。
【図５】図５は、周波数とスペクトル密度との関係を示す図である。
【図６】図６は、各極大周波数および各極大周波数に対応するスペクトル密度を時系列で表した図である。
【図７】図７は、クラスタリングデータ記憶部が記憶するクラスタリングデータのデータ構造の一例を示す図である。
【図８】図８は、眠気判定テーブルのデータ構造の一例を示す図である。
【図９】図９は、スケールを最適に調整されたグラフの一例を示す図である。
【図１０】図１０は、グラフ調整部の処理を説明するための図である。
【図１１】図１１は、覚醒度と周波数とを軸とするグラフの一例を示す図である。
【図１２】図１２は、本実施例１にかかる覚醒度判定装置の処理手順を示すフローチャートである。
【図１３】図１３は、実施例にかかる覚醒度判定装置を構成するコンピュータのハードウェア構成を示す図である。
【発明を実施するための形態】
【００１２】
以下に添付図面を参照して、この発明に係る覚醒度判定装置および覚醒度判定方法の好適な実施の形態を詳細に説明する。
【実施例１】
【００１３】
まず、本実施例１にかかる覚醒度判定装置の概要について説明する。図１は、本実施例１にかかる覚醒度判定装置の概要を説明するための図である。覚醒度判定装置は、まず、被験者（例えば、車を運転中の運転手）の心拍信号から振幅ピークを検出し、検出した各振幅ピークの間隔を心拍間隔として検出する。
【００１４】
そして、覚醒度判定装置は、時間経過に伴う心拍間隔の変動に基づいて当該変動の周波数毎にスペクトル密度を算出し、算出したスペクトル密度が極大となる特徴点を判定する。覚醒度判定装置は、周波数と覚醒度とを軸とするグラフを作成し（図１参照）、被験者情報を基にして、グラフ上の領域を眠気レベルに応じて分割すると共に、特徴点の移動範囲に応じて、グラフのスケールを調整する。
【００１５】
ここで、被験者情報には、被験者の行動に関する情報が含まれており、例えば、被験者の顔画像から抽出した瞼の閉じ具合や、車線逸脱時間、フラフラ運転の度合等が含まれる。また、図１に示す縦軸の覚醒度は、スペクトル密度の値の上下関係を反転したものである（覚醒度とスペクトル密度との関係は、スペクトル密度が大きくなると、覚醒度が小さくなるという関係にある）。
【００１６】
図１に示す例では、覚醒度判定装置は、眠気レベルが低い順に眠気レベル１〜眠気レベル５にグラフの領域を分割している。ここでは一例として、眠気レベルは、レベル５＞レベル４＞レベル３＞レベル２＞レベル１とする。ここで、眠気レベル１の覚醒度が最も高く、眠気レベル５の覚醒度が最も低くなる。なお、覚醒度が高い場合には、被験者は目が覚めている状態であり、覚醒度が低い場合には、被験者は眠い状態にある。
【００１７】
続いて、覚醒度判定装置は、被験者の心拍信号に応じて、時々刻々と変化する特徴点をグラフ上にプロットし、被験者の覚醒度を判定する。例えば、覚醒度判定装置は、現在の特徴点のグラフ上の位置が眠気レベル２の状態であり、かかる特徴点が眠気レベル２から眠気レベル３方向に移動した場合には、被験者は眠気が強くなっており危険であるため、被験者に何らかのシステム応答（アクション）をおこす（アクションの具体的内容は後述する）。一方、現在の特徴点のグラフ上の位置が眠気レベル３の状態であっても、かかる特徴点が眠気レベル３から眠気レベル２の方向に移動した場合には、被験者は目が覚めつつあるので、被験者の気分を害さないために、被験者に対するアクションを起こさない。
【００１８】
このように、本実施例１にかかる覚醒度判定装置は、被験者にあわせてグラフを眠気レベル毎に複数の領域に分割し、グラフ上の特徴点に基づいて被験者が実際に眠くなっているか否かを判定するので、被験者毎に正確な覚醒度を算出することができる。
【００１９】
また、本実施例１にかかる覚醒度判定装置は、グラフ上の特徴点の移動方向に応じて、被験者に対するアクションを切り替えるので、被験者の気分を害すことなく、被験者に各種のアクションを与え、運転事故等を予防することができる。
【００２０】
次に、本実施例１にかかる覚醒度判定装置の構成について説明する。図２は、本実施例１にかかる覚醒度判定装置１００の構成を示す機能ブロック図である。図２に示すように、この覚醒度判定装置１００は、心拍検出部１０１と、心拍間隔算出部１０２と、スペクトル算出部１０３と、ピーク周波数算出部１０４と、ピーク周波数クラスタリング部１０５と、クラスタリングデータ記憶部１０５ａと、眠気レベル判定部１０６と、グラフ調整部１０７と、システム応答判定部１０８と、出力部１０９とを有する。
【００２１】
このうち、心拍検出部１０１は、被験者の心拍信号を検出する手段である。この心拍検出部１０１は、例えば、被験者に接触している電極に対して電圧を印加し、各電極の電位差から被験者の心拍信号を取得することができる。図３は、心拍検出部１０１が検出する心拍信号の一例を示す図である。心拍検出部１０１は、心拍信号のデータ（以下、心拍信号データ）を心拍間隔算出部１０２に出力する。
【００２２】
心拍間隔算出部１０２は、心拍信号データに基づいて心拍信号の振幅ピークを検出し、検出したタイミングの間隔（以下、心拍間隔と表記する）を検出する手段である。図３を用いて心拍間隔算出部１０２の処理を説明する。同図に示すように、心拍間隔算出部１０２は、心拍信号の振幅が閾値以上となるポイントＲ、すなわち振幅ピークを検出し、検出した各ポイントＲの間隔（Ｒ−Ｒインターバルに対応する時間間隔）を算出する。
【００２３】
心拍間隔算出部１０２は、算出した心拍間隔のデータ（以下、心拍間隔データと表記する）をスペクトル算出部１０３に出力する。なお、振幅ピークの検出方法は上述の方法に限るものではなく、例えば、心拍信号の微分係数が正から負に変わるゼロクロス点を使う方法、振幅波形につきパターンマッチングを行ってピークを検出する方法などを用いても構わない。
【００２４】
スペクトル算出部１０３は、心拍間隔データを基にして心拍間隔の変動に対するスペクトル密度を算出する手段である。ここで、スペクトル算出部１０３の行う処理を具体的に説明すると、まず、スペクトル算出部１０３は、心拍間隔データを基にして、時間変化によって変動する心拍間隔のデータ（以下、心拍間隔変動データと表記する）を生成する。図４は、心拍間隔変動データを心拍間隔−時間平面上で示した図である。図４の下段に示すように、心拍間隔は時間変化に伴って変動している。
【００２５】
続いて、スペクトル算出部１０３は、心拍間隔変動データを基にして各周波数に対応するスペクトル密度を算出する。図５は、周波数とスペクトル密度との関係を示す図である。同図に示すように、スペクトル密度は複数の周波数において極大となっている。スペクトル算出部１０３は、算出したスペクトル密度と周波数との関係を示すデータ（以下、スペクトル密度データと表記する）をピーク周波数算出部１０４に出力する。
【００２６】
なお、スペクトル算出部１０３がスペクトル密度を算出する方法はどのような方法を用いても構わないが、例えば、ＡＲ（Autoregressive）モデルを用いてスペクトル密度を算出することができる。ＡＲモデルは、非特許文献（佐藤俊輔、吉川昭、木竜徹、”生体信号処理の基礎”、コロナ社）などに開示されてあるように、ある時点の状態を過去の時系列データの線形和で表すモデルであり、フーリエ変換と比較して少ないデータ数でも明瞭な極大点が得られるという特徴がある。
【００２７】
時系列ｘ（ｓ）のｐ次のＡＲモデルは、過去の値に対する重みであるＡＲ係数ａ（ｍ）および誤差項ｅ（ｓ）を用いて
【数１】
2010155072.tif 000003
によって表すことができる（理想的にはｅ（ｓ）はホワイトノイズである）。
【００２８】
そして、ｐを同定次数、ｆsをサンプリング周波数、εpを同定誤差とし、
【数２】
2010155072.tif 000004
をｋ次のＡＲ係数とすると、スペクトル密度ＰAR（ｆ）は
【数３】
2010155072.tif 000005
によって表すことができる。スペクトル算出部１０３は、式（２）および心拍間隔変動データを基にして、スペクトル密度データを算出することができる。
【００２９】
図２の説明に戻ると、ピーク周波数算出部１０４は、スペクトル密度データを基にしてスペクトル密度が極大となる周波数（以下、極大周波数と表記する）および極大周波数に対応するスペクトル密度（以下、極大スペクトル密度と表記する）を算出する手段である。ピーク周波数算出部１０４は、
【数４】
2010155072.tif 000006
となる周波数ｆを極大周波数として算出し、この極大周波数を式（２）に代入することによって極大スペクトル密度を得る。
【００３０】
図６は、各極大周波数および各極大周波数に対応するスペクトル密度を時系列で表した図である。例えば、スペクトル算出部１０３が１０秒間隔でスペクトル密度データを算出している場合には、図６に示す時系列方向の点の間隔は１０秒間隔となる。ピーク周波数算出部１０４は、極大周波数および極大スペクトル密度のデータをピーク周波数クラスタリング部１０５に出力する。
【００３１】
ピーク周波数クラスタリング部１０５は、極大周波数および極大スペクトル密度を複数のクラスタに分類し、分類した結果をクラスタリングデータ記憶部１０５ａに記憶する手段である。図７は、クラスタリングデータ記憶部１０５ａが記憶するクラスタリングデータのデータ構造の一例を示す図である。
【００３２】
図７に示すように、このクラスタリングデータは、極大周波数が類似する（所定帯域に含まれる）極大周波数ごとに各クラスタに分類されており、それぞれの極大周波数には極大スペクトル密度が対応付けられている。このクラスタリングデータは、計測開始後一定時間の極大周波数および極大スペクトル密度のデータを取得して初期登録してもよいし、個人ごとに予め作成されたデータベースから個人固有の極大周波数および極大スペクトル密度のデータを基にして初期登録してもよい。
【００３３】
また、各クラスタには、各クラスタに代表周波数および代表スペクトル密度が設定されている。代表周波数は、クラスタに含まれる極大周波数の代表値であり、代表スペクトル密度は、クラスタに含まれる極大スペクトル密度の代表値である。
【００３４】
ピーク周波数クラスタリング部１０５は、代表周波数を、クラスタに含まれる極大周波数の平均によって算出することができ、代表スペクトル密度も、クラスタに含まれる極大スペクトル密度の平均によって算出することができる。ピーク周波数クラスタリング部１０５は、クラスタリングデータに含まれる各クラスタ１〜４のうち、予め設定されたクラスタを選択し、選択したクラスタの代表周波数と代表スペクトル密度との関係を特徴点データとして、グラフ調整部１０７およびシステム応答判定部１０８に出力する。
【００３５】
眠気レベル判定部１０６は、車内に取り付けられたカメラ（図示略）を用いて被験者の顔をモニタリングし、被験者の眠気レベルを判定する手段であり、判定結果をグラフ調整部１０７に出力する。具体的に、眠気レベル判定部１０６は、被験者の瞼の状態（例えば、瞼が閉じている時間、瞼が開閉する周期）と眠気判定テーブルとを比較して、被験者の眠気レベルを判定する。図８は、眠気判定テーブルのデータ構造の一例を示す図である。
【００３６】
図８に示すように、この眠気判定テーブルは、眠気レベルと瞼閉時間（瞼が閉じられている間の時間）と傾向とが対応付けられて記憶されている。眠気レベル判定部１０６は、被験者の瞼閉時間が５秒以上の場合には、被験者の眠気レベルは眠気レベル５であると判定する。
【００３７】
眠気レベル判定部１０６は、被験者の瞼閉時間が２秒以上５秒未満の場合、あるいは、所定時間あたりの、瞼の大きな開閉行動（たとえば、上瞼の頂点と下瞼の頂点の距離が閾値以上となる瞼の開閉行動）が断続的に増加する場合には、被験者の眠気レベルは眠気レベル４であると判定する。
【００３８】
眠気レベル判定部１０６は、被験者の瞼閉時間が１秒以上２秒未満の場合、あるいは、所定時間あたりの、瞼の閉時間が増加傾向で断続的な増加である場合には、被験者の眠気レベルは眠気レベル３であると判定する。
【００３９】
眠気レベル判定部１０６は、被験者の瞼閉時間が０．１秒以上１秒未満の場合、あるいは、所定時間あたりの、瞼の開時間が増加傾向の場合には、被験者の眠気レベルは眠気レベル２であると判定する。眠気レベル判定部１０６は、被験者の瞼の開時間および閉時間が周期的である場合には、眠気レベルは眠気レベル１であると判定する。上記の秒数の閾値レベルは、一例であって、レベル２を０．１〜０．５秒、レベル３を０．５〜１．０秒、レベル４を１．０〜２秒、レベル５を２秒以上などとしてもよく、フィールド試験やユーザからの要望などから、最適値を設定することができる。
【００４０】
ところで、眠気レベル判定部１０６が被験者の眠気レベルを判定する処理は、上記の処理に限られるものではない。例えば、眠気レベル判定部１０６は、被験者のあくびに基づいて眠気レベルを判定しても良い。眠気レベル判定部１０６が被験者のあくびを検出する方法としては、例えば、被験者の顔をモニタリングして口元を監視する方法と、被験者の脈拍を監視する方法がある。
【００４１】
まず、眠気レベル判定部１０６が被験者の顔をモニタリングして口元を監視することで、被験者のあくびを検出する方法について説明する。例えば、眠気レベル判定部１０６は、被験者の上唇と下唇とが閾値以上はなれている場合に、被験者が口を開けていると判定する。そして、眠気レベル判定部１０６は、被験者の口が所定時間以上開いている場合に、被験者があくびをしたと判定してもよい。または、眠気レベル判定部１０６は、被験者の口が所定の時間以上の周期で開閉した場合に、被験者があくびをしたと判定しても良い。
【００４２】
続いて、眠気レベル判定部１０６が被験者の脈拍を監視することで、被験者のあくびを検出する方法について説明する。例えば、眠気レベル判定部１０６は、被験者があくびをする際の心拍変動に現れる特徴的な波形を記憶しておき、かかる波形と被験者から検出される波形とを比較して、被験者があくびをしたか否かを判定しても良い。または、眠気レベル判定部１０６は、特開２００７−３１９３７８号公報に開示された技術を利用して、被験者の脈拍から被験者があくびをしたか否かを判定しても良い。
【００４３】
ここで、眠気レベル判定部１０６が被験者のあくびの状態により眠気レベルを判定する場合の一例について説明する。眠気レベル判定部１０６は、被験者のあくびを初めて検出した場合に、眠気レベルを「レベル３」と判定する。また、眠気レベル判定部１０６は、１分間に被験者のあくびを２回以上検出した場合に、眠気レベルを「レベル４」と判定する。眠気レベル判定部１０６は、眠気レベルの判定結果をグラフ調整部１０７に出力する。
【００４４】
グラフ調整部１０７は、スペクトル密度（覚醒度）および周波数を軸とするグラフを作成し、眠気レベル判定部１０６の判定結果に基づいてグラフ上の領域を眠気レベル毎に分割する手段である。以下において、グラフ調整部１０７の処理を具体的に説明する。
【００４５】
まず、グラフ調整部１０７は、グラフのスケールをデフォルトに設定し、所定時間の間、ピーク周波数クラスタリング部１０５から順次出力される特徴点データを基にして、グラフ上に特徴点をプロットする。グラフ調整部１０７は、グラフ上にプロットした各特徴点の中から、最上点、最下点、最右点、最左点を検出する。
【００４６】
ここで、最上点は、グラフにプロットした特徴点の内、最もスペクトル密度が大きい特徴点であり、最下点は、グラフにプロットした特徴点の内、最もスペクトル密度が小さい特徴点である。最右点は、グラフにプロットした特徴点の内、最も周波数が大きい特徴点であり、最左点は、グラフにプロットした特徴点の内、最も周波数が小さい特徴点となる。
【００４７】
グラフ調整部１０７は、検出した最上点、最下点、最右点、最左点が、グラフ上に最適に含まれるように、グラフのスケールを調整する。例えば、グラフ調整部１０７は、グラフに含める周波数の最小値を、最左点の周波数から所定値減算した周波数とし、グラフに含める周波数の最大値を、最右点の周波数に所定値加算した周波数にする。
【００４８】
また、グラフ調整部１０７は、グラフに含めるスペクトル密度の最小値を、最下点のスペクトル密度から所定値減算したスペクトル密度とし、グラフに含めるスペクトル密度の最大値を、最上点のスペクトル密度に所定値加算したスペクトル密度とする。図９は、スケールを最適に調整されたグラフの一例を示す図である。
【００４９】
続いて、グラフ調整部１０７は、グラフ（スケールを最適に調整されたグラフ）上の領域を眠気レベル毎に分割する。図１０は、グラフ調整部１０７の処理を説明するための図である。まず、グラフ調整部１０７は、デフォルト値に基づいて、グラフを眠気レベル１〜５に分割する。
【００５０】
そして、グラフ調整部１０７は、ピーク周波数クラスタリング部１０５から出力される特徴点データをグラフにプロットし、グラフにプロットした領域と、眠気レベル判定部１０６の判定結果に基づいて、グラフの眠気レベルの領域を調整する。
【００５１】
例えば、グラフ調整部１０７は、特徴点Ａ（図１０参照）をグラフにプロットした時点における、眠気レベル判定部１０６の判定結果が、眠気レベル５の場合には、判定結果と特徴点Ａの眠気レベルとが等しいので、グラフの眠気レベルの領域をそのままの状態にする。
【００５２】
また、例えば、グラフ調整部１０７がグラフに特徴点データをプロットした時点において、眠気レベル判定部１０６が被験者のあくびを初めて検出した場合のグラフの調整について説明する。この場合、眠気レベルは「レベル３」となるので、グラフ調整部１０７は、眠気レベル２と３との境界線が、プロットした特徴点データを通過するように、境界線を平行移動させる。なお、グラフ調整部１０７は、眠気レベル２と３との境界線を平行移動させた場合に、かかる移動の距離に合わせて、眠気レベル１と２との境界線、眠気レベル３と４との境界線、眠気レベル４と５との境界線を平行移動させても良い。
【００５３】
また、グラフ調整部１０７が、グラフ調整部１０７がグラフに特徴点データをプロットした時点において、被験者のあくびを基に眠気レベルを「レベル４」と判定した場合のグラフの調整について説明する。この場合、グラフ調整部１０７は、眠気レベル４と５との境界線が、プロットした特徴点データを通過するように、境界線データを平行移動させる。なお、グラフ調整部１０７は、眠気レベル４と５との境界線を平行移動させた場合に、かかる移動の距離に合わせて、眠気レベル１と２との境界線、眠気レベル２と３との境界線、眠気レベル３と４との境界線を平行移動させても良い。
【００５４】
一方、グラフ調整部１０７は、特徴点Ｂ（図１０参照）をグラフにプロットした時点における、眠気レベル判定部１０６の判定結果が、眠気レベル３の場合には、判定結果と特徴点Ｂの眠気レベルとが異なるので、グラフの眠気レベル領域を調整する。この場合には、グラフ上の眠気レベル４の領域を特徴点Ｂまで押し下げる。
【００５５】
グラフ調整部１０７は、調整したグラフのデータをシステム応答判定部１０８および出力部１０９に出力する。なお、グラフ調整部１０７は、グラフのデータを出力する場合に、スペクトル密度の上下関係を反転したデータを作成し、作成したデータをシステム応答判定部１０８および出力部１０９に出力してもよい。スペクトル密度の上下関係を反転することにより、スペクトル密度と周波数とを軸とするグラフが、覚醒度と周波数とを軸とするグラフとなる。図１１は、覚醒度と周波数とを軸とするグラフの一例を示す図である。
【００５６】
システム応答判定部１０８は、グラフ調整部１０７から出力されるグラフのデータと、ピーク周波数クラスタリング部１０５から出力される特徴点データとを取得し、グラフの各眠気レベルの領域を移動する特徴点に基づいて、被験者に警報アクションを行う手段である。
【００５７】
例えば、システム応答判定部１０８は、グラフ上の特徴点が眠気レベル２の領域から眠気レベル３に移行した場合には、被験者は眠りにつきつつあり危険な状態であるため、音声による注意発起を行う（マイクから「少し休憩をしたほうがよさそうです」等の音声を出力する）。一方、システム応答判定部１０８は、グラフ上の特徴点が眠気レベル３の領域から眠気レベル２の領域に移行した場合には、被験者は目が覚めつつあるので、被験者の気分を害さないために、被験者に対するアクションを起こさない。
【００５８】
出力部１０９は、各種の情報を出力する出力手段であり、例えば、ディスプレイ、マイクなどが含まれる。出力部１０９は、グラフ調整部１０７からグラフのデータを取得した場合には、取得したグラフをディスプレイに表示する。また、出力部１０９は、システム応答判定部１０８からの制御命令に応じて、各種のアクションを被験者に対して起こす（例えば、音声による注意発起を行う）。
【００５９】
次に、本実施例１にかかる覚醒度判定装置１００の処理手順について説明する。図１２は、本実施例１にかかる覚醒度判定装置１００の処理手順を示すフローチャートである。図１２に示すように、覚醒度判定装置１００は、心拍検出部１０１が心拍信号データを取得し（ステップＳ１０１）、心拍間隔算出部１０２が心拍間隔を検出する（ステップＳ１０２）。
【００６０】
そして、スペクトル算出部１０３が、各周波数に対応するスペクトル密度を算出し（ステップＳ１０３）、ピーク周波数算出部１０４が、スペクトル密度が極大となる周波数を算出し（ステップＳ１０４）、ピーク周波数クラスタリング部１０５が、極大周波数をクラスタリングする（ステップＳ１０５）。
【００６１】
続いて、グラフ調整部１０７は、ピーク周波数クラスタリング部１０５から出力される特徴点データの移動範囲に基づいて、グラフのスケールを調整し（ステップＳ１０６）、眠気レベル判定部１０６から出力される判定結果に基づいて、グラフ上に眠気レベルをマッピングする（ステップＳ１０７）。システム応答判定部１０８は、グラフのデータと特徴点とを基にして、眠気レベルを判定し（ステップＳ１０８）、判定結果に応じてユーザにアクションを起こす（ステップＳ１０９）。
【００６２】
上述してきたように、本実施例１にかかる覚醒度判定装置１００は、グラフ調整部１０７が、被験者にあわせてグラフを眠気レベル毎に複数の領域に分割し、システム応答判定部１０８が、グラフ上の特徴点に基づいて被験者が実際に眠くなっているか否かを判定するので、被験者毎に正確な覚醒度を算出することができる。
【００６３】
また、本実施例１にかかる覚醒度判定装置１００は、システム応答判定部１０８が、グラフ上の特徴点の移動方向に応じて、被験者に対するアクションを切り替えるので、被験者の気分を害すことなく、被験者に各種のアクションを与え、運転事故等を予防することができる。
【実施例２】
【００６４】
さて、これまで本発明の実施例について説明したが、本発明は上述した実施例１以外にも種々の異なる形態にて実施されてよいものである。そこで、以下では実施例２として本発明に含まれる他の実施例を説明する。
【００６５】
（１）被験者に対するアクション
上述した実施例１では、システム応答判定部１０８が、グラフ上の各眠気レベルの領域を移動する特徴点に基づいて、被験者に対して警報を出力していたが、システム応答判定部１０８が起こすアクションは、実施例１の内容に限定されるものではない。以下において、システム応答判定部１０８が実行するアクションのバリエーションについて説明する。なお、システム応答判定部１０８が実行するアクションには、音声（音響）や音楽による注意発起の他に、車内に酸素を送り込み酸素濃度を増加させるアクション、ミストを車内に送り込み湿度を上昇させるアクション、被験者に対する帯電を防止・あるいは除去するアクション、ノイズ増強を行うアクション、香りを含んだ空気を車内あるいは被験者に送り込む、ハンドルやシートを振動させる、窓の開閉を行う、軽い電気ショックを起こさせる、シートベルトを軽く引っ張るなどのアクションが含まれる。
【００６６】
眠気レベルが２、１、２、・・・と断続的に変化している場合、すなわち、特徴点が眠気レベル２と１の境界を行ったりきたりしている場合には、システム応答判定部１０８は、被験者は自発的に目覚めはじめていると判定し、強度を下げてアクションを実行する（例えば、音量を絞って、注意発起を行う）、もしくは、アクションを実行しない。
【００６７】
眠気レベルが２、２、・・・、３、２のように変化し、一瞬眠気レベル３になり、眠気レベル２に復帰するように特徴点が移動している場合には、システム応答判定部１０８は、被験者は眠気が増大していることを自覚していないと判定し、強度を上げてアクションを実行する（例えば、音量を上げて、注意発起を行う）と共に、他のアクション（例えば、車内の酸素濃度を上げる）をあわせて実行する。
【００６８】
眠気レベルが２、３、３、・・・、３のように変化し、特徴点が定常的に眠気レベル３に停滞している場合には、更に違うアクションを提示する。例えば、システム応答判定部１０８は、「眠気レベルが上昇しています。車線変更などを行い、気分転換することをお勧めします」、等の音声を出力して被験者に車線変更を促すことで、被験者の覚醒度の上昇を図る。また、システム応答判定部１０８は、準備してあった「アイテム」の取得を促すと共に、カーナビゲーションシステムと連動し、予め登録しておいた被験者の嗜好傾向にあるメニュー情報を音声で提示し、会話による覚醒効果を狙う。ここで、アイテムとしては、ガム、ドリンク、アメ、通信カラオケ、タバコなどの嗜好品が含まれる。
【００６９】
眠気レベルが３、４、３、・・・、３、４のように変化し、特徴点が眠気レベル４の領域まで達するようになった場合には、システム応答判定部１０８は、強い注意発起を行う。例えば、システム応答判定部１０８は、「棚上げ効果」を狙い、周辺車両状況として他車両の接近を強調して通知する。また、準備してあった「アイテム」を強く勧めると共に、新たなアイテムの導入を提案する。
【００７０】
眠気レベル４，５が頻発する場合には、警報判定装置１０８は、被験者自身も強く眠気を自覚するほどの眠気を検知したことを被験者に通知し、安全のため休憩を促すと共に、周辺の車両に眠気注意者であると通知する。また、警報アクションとしてあらゆる手段の強度を上げる。
【００７１】
眠気レベルが３、２、３、２、・・・に変化し、特徴点が徐々に覚醒方向（眠気レベル１方向に近付いている場合）に近づいている場合は、さらに覚醒方向への遷移を助長させるアイテムを提示する。ここでの助長アイテムとしては、オペレータとの会話、カラオケなど、おしゃべりによる覚醒効果を得られるアイテムが含まれる。
【００７２】
このように、眠気レベルの変化の特徴に応じて、システム応答判定部１０８がアクションを切り替え、被験者に注意発起を行うので、被験者の気分を害すことなく、眠気による危険を防止することができる。
【００７３】
なお、警報判定装置１０８は、被験者の嗜好と覚醒効果とを予め登録しておき、高速道路等のサービスエリアで放送されているＣＭを受信すると共に、受信したＣＭの中から被験者の覚醒度を高める内容のＣＭを選択して、車内のスピーカから出力させても良い。
【００７４】
（２）システムの構成など
ところで、本実施例において説明した各処理のうち、自動的に行われるものとして説明した処理の全部または一部を手動的に行うこともでき、あるいは、手動的に行われるものとして説明した処理の全部あるいは一部を公知の方法で自動的に行うこともできる。この他、上記文書中や図面中で示した処理手順、制御手順、具体的名称、各種のデータやパラメータを含む情報については、特記する場合を除いて任意に変更することができる。
【００７５】
また、実施例に示した覚醒度判定装置１００の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部または一部を、各種の負荷や使用状況などに応じて、任意の単位で機能的または物理的に分散・統合して構成することができる。さらに、各装置にて行われる各処理機能は、その全部または任意の一部がＣＰＵおよび当該ＣＰＵにて解析実行されるプログラムにて実現され、あるいは、ワイヤードロジックによるハードウェアとして実現され得る。
【００７６】
図１３は、実施例にかかる覚醒度判定装置を構成するコンピュータ１０のハードウェア構成を示す図である。図１３に示すように、このコンピュータ（覚醒度判定装置）１０は、入力装置１１、モニタ１２、他の装置とデータ通信を行う通信制御装置１３、ＲＡＭ（Random  Access  Memory）１４、ＲＯＭ（Read  Only  Memory）１５、被験者から心拍信号を取得する心拍センサ１６、各種のアクションを実行する警報出力装置１７、ＣＰＵ（Central  Processing  Unit）１８、ＨＤＤ（Hard  Disk  Drive）１９をバス２０で接続している。
【００７７】
そして、ＨＤＤ１９には、上述した覚醒度判定装置１００の機能と同様の機能を発揮する覚醒度判定プログラム１９ｂが記憶されている。ＣＰＵ１８が、覚醒度判定プログラム１９ｂを読み出して実行することにより、覚醒度判定プロセス１８ａが起動される。ここで、覚醒度判定プロセス１８ａは、実施例１に示した心拍検出部１０１、心拍間隔算出部１０２、スペクトル算出部１０３、ピーク周波数算出部１０４、ピーク周波数クラスタリング部１０５、眠気レベル判定部１０６、グラフ調整部１０７、システム応答判定部１０８に対応する。
【００７８】
なお、ＨＤＤ１９は、実施例１のクラスタリングデータ記憶部１０５ａ等に記憶されたデータに対応する各種データ１９ａを記憶している。ＣＰＵ１８は、ＨＤＤ１９に記憶された各種データ１９ａをＲＡＭ１４に読み出し、各種データ１４ａおよび入力装置等から入力される被験者の行動に関する情報に基づいて、被験者の眠気レベルを判定し、警報出力装置１７を制御して各種のアクションを被験者に対して実行する。
【００７９】
ところで、図１３に示した覚醒度判定プログラム１９ｂは、必ずしも最初からＨＤＤ１９に記憶させておく必要はない。たとえば、コンピュータに挿入されるフレキシブルディスク（ＦＤ）、ＣＤ−ＲＯＭ、ＤＶＤディスク、光磁気ディスク、ＩＣカードなどの「可搬用の物理媒体」、または、コンピュータの内外に備えられるハードディスクドライブ（ＨＤＤ）などの「固定用の物理媒体」、さらには、公衆回線、インターネット、ＬＡＮ、ＷＡＮなどを介してコンピュータに接続される「他のコンピュータ（またはサーバ）」などに覚醒度判定プログラム１９ｂを記憶しておき、コンピュータがこれらから覚醒度判定プログラム１９ｂを読み出して実行するようにしてもよい。
【００８０】
以上の実施例を含む実施形態に関し、更に以下の付記を開示する。
【００８１】
（付記１）被験者の心拍信号から振幅ピークを検出し、検出した各振幅ピークの間隔を心拍間隔として検出する心拍間隔検出手段と、
時間経過に伴う前記心拍間隔の変動に基づいて当該変動の周波数毎にスペクトル密度を算出し、算出したスペクトル密度が極大となる特徴点を判定する特徴点判定手段と、
前記被験者の行動に関する情報を含んだ被験者情報と前記特徴点とを基にして前記被験者の覚醒度を判定する覚醒度判定手段と
を有する覚醒度判定装置。
【００８２】
（付記２）前記覚醒度判定手段は、前記周波数とスペクトル密度とを軸とするグラフを作成し、当該グラフ上の前記特徴点の位置と前記被験者情報とを基にして、前記グラフ上の領域を前記覚醒度に応じて分割するグラフ作成手段と、前記グラフ上に前記特徴点をプロットし、当該特徴点の位置する領域に基づいて前記被験者の覚醒度を判定する判定手段とを有する付記１に記載の覚醒度判定装置。
【００８３】
（付記３）前記グラフ作成手段は、前記特徴点の移動範囲に基づいて前記グラフのスケールを調整する付記２に記載の覚醒度判定装置。
【００８４】
（付記４）前記グラフ上における前記特徴点の移動方向を基にして前記被験者に警告を出力するか否かを判定する警告判定手段と、前記警告判定手段の判定結果に基づいて前記被験者にシステム応答を実行する手段とを更に有する付記２または３に記載の覚醒度判定装置。
【００８５】
（付記５）覚醒度判定装置が、
被験者の心拍信号から振幅ピークを検出し、検出した各振幅ピークの間隔を心拍間隔として検出して記憶装置に記憶する心拍間隔検出ステップと
時間経過に伴う前記心拍間隔の変動に基づいて当該変動の周波数毎にスペクトル密度を算出し、算出したスペクトル密度が極大となる特徴点を判定する特徴点判定ステップと、
前記被験者の行動に関する情報を含んだ被験者情報と前記特徴点とを基にして前記被験者の覚醒度を判定する覚醒度判定ステップと
を含んだ覚醒度判定方法。
【００８６】
（付記６）前記覚醒度判定ステップは、前記周波数とスペクトル密度とを軸とするグラフを作成し、当該グラフ上の前記特徴点の位置と前記被験者情報とを基にして、前記グラフ上の領域を前記覚醒度に応じて分割するグラフ作成ステップと、前記グラフ上に前記特徴点をプロットし、当該特徴点の位置する領域に基づいて前記被験者の覚醒度を判定する判定ステップとを含む付記５に記載の覚醒度判定方法。
【００８７】
（付記７）前記グラフ作成ステップは、前記特徴点の移動範囲に基づいて前記グラフのスケールを調整する付記６に記載の覚醒度判定方法。
【００８８】
（付記８）前記グラフ上における前記特徴点の移動方向を基にして前記被験者に警告を出力するか否かを判定する警告判定ステップと、前記警告判定ステップの判定結果に基づいて前記被験者にシステム応答を実行するステップとを更に含む付記６または７に記載の覚醒度判定方法。
【００８９】
（付記９）コンピュータに、
被験者の心拍信号から振幅ピークを検出し、検出した各振幅ピークの間隔を心拍間隔として検出して記憶装置に記憶する心拍間隔検出手順と、
時間経過に伴う前記心拍間隔の変動に基づいて当該変動の周波数毎にスペクトル密度を算出し、算出したスペクトル密度が極大となる特徴点を判定する特徴点判定手順と、
前記被験者の行動に関する情報を含んだ被験者情報と前記特徴点とを基にして前記被験者の覚醒度を判定する覚醒度判定手順と
を実行させる覚醒度判定プログラム。
【００９０】
（付記１０）前記覚醒度判定手順は、前記周波数とスペクトル密度とを軸とするグラフを作成し、当該グラフ上の前記特徴点の位置と前記被験者情報とを基にして、前記グラフ上の領域を前記覚醒度に応じて分割するグラフ作成手順と、前記グラフ上に前記特徴点をプロットし、当該特徴点の位置する領域に基づいて前記被験者の覚醒度を判定する判定手順とを含む付記９に記載の覚醒度判定プログラム。
【符号の説明】
【００９１】
１０    コンピュータ
１１    入力装置
１２    モニタ
１３    通信制御装置
１４    ＲＡＭ
１４ａ，１９ａ  各種データ
１５    ＲＯＭ
１６    心拍センサ
１７    警報出力装置
１８    ＣＰＵ
１８ａ  覚醒度判定プロセス
１９    ＨＤＤ
１９ｂ  覚醒度判定プログラム
１００    覚醒度判定装置
１０１    心拍検出部
１０２    心拍間隔算出部
１０３    スペクトル算出部
１０４    ピーク周波数算出部
１０５    ピーク周波数クラスタリング部
１０６    眠気レベル判定部
１０７    グラフ調整部
１０８    システム応答判定部
１０９    出力部
【図１】
2010155072.tif 000007
【図２】
2010155072.tif 000008
【図３】
2010155072.tif 000009
【図４】
2010155072.tif 000010
【図５】
2010155072.tif 000011
【図６】
2010155072.tif 000012
【図７】
2010155072.tif 000013
【図８】
2010155072.tif 000014
【図９】
2010155072.tif 000015
【図１０】
2010155072.tif 000016
【図１１】
2010155072.tif 000017
【図１２】
2010155072.tif 000018
【図１３】
2010155072.tif 000019