(57)【特許請求の範囲】
【請求項１】
複軸センサから得られる時系列の複軸センサデータにおいて検出したい複数の特定信号パターンを辞書データとして予め獲得するとともに、その獲得された特定信号パターンを予め用意された時間窓に基づいて第１の時間分割センサデータとして分割する第１のデータ分割手段と、
前記第１の時間分割センサデータのノイズ除去ならびにデータを正規化することで第１の正規化センサデータを生成する第１の正規化センサデータ生成手段と、
前記第１の正規化センサデータを基底関数展開することで時空間解像度レベルごとに第１の時空間分解センサデータを生成する第１の時空間分解センサデータ生成手段と、
前記第１の時空間分解センサデータの振幅情報に基づいて時空間解像度レベルごとに第１の時空間分解センサデータ情報エントロピーを算出する第１のエントロピー算出手段と、
前記第１の時空間分解センサデータ情報エントロピーの総和が最小になるように予め用意された評価基準に基づいて第１の基底時空間分解センサデータを選択する第１の基底時空間分解センサデータ選択手段と、
前記第１の基底時空間分解センサデータの情報エントロピーに基づいて、前記辞書データの特徴ベクトルである第１のセンサデータ特徴ベクトルを生成する第１のセンサデータ特徴ベクトル生成手段と、
前記第１のセンサデータ特徴ベクトルを自己組織化アルゴリズムにより学習させることで自己組織化された第１の特定信号パターン特徴マップを生成する第１の特定信号パターン特徴マップ生成手段と、
前記第１の特定信号パターン特徴マップの各ノードにおける特徴ノードベクトルに基づいて、前記辞書データの特徴ベクトルを入力した時に発火する特徴ノードベクトルとの対応関係を第１の蓄積結果として蓄積する蓄積手段と、
前記第１の蓄積結果に基づいて、前記第１の特定信号パターン特徴マップの各ノードごとに特定信号パターンが発火する確率を算出し第１の特定信号パターン発火確率マップを生成する第１の特定信号パターン発火確率マップ生成手段と、
前記第１の特定信号パターン発火確率マップのうち予め設定された閾値以上の確率を有するノードを基底特徴ノードベクトルとして抽出する基底特徴ノードベクトル抽出手段と、
前記基底特徴ノードベクトルを発火させる基となった複軸センサデータ特徴ベクトルのうち、基底特徴ノードベクトルとの類似度が予め設定された閾値以上のものを基底複軸センサ特徴ベクトルとして抽出する基底複軸センサ特徴ベクトル抽出手段と、
前記基底複軸センサ特徴ベクトルの基となった複軸の基底センサデータを抽出する基底センサデータ抽出手段と、
前記複軸の基底センサデータを予め設定された刻み幅の角度で各軸回りに回転させることで擬似複軸センサデータを生成する擬似複軸センサデータ生成手段と、
を含み、
さらに
前記第１の時空間分解センサデータ生成手段は、前記疑似複軸センサデータを基底関数展開することで時空間解像度レベルごとに第２の時空間分解センサデータを生成し、
前記第１のエントロピー算出手段は、前記第２の時空間分解センサデータの振幅情報に基づいて時空間解像度レベルごとに第２の時空間分解センサデータ情報エントロピーを算出し、
前記第１の基底時空間分解センサデータ選択手段は、前記第２の時空間分解センサデータ情報エントロピーの総和が最小になるように予め用意された評価基準に基づいて第２の基底時空間分解センサデータを選択し、
前記第１のセンサデータ特徴ベクトル生成手段は、前記第２の基底時空間分解センサデータの情報エントロピーに基づいて第２のセンサデータ特徴ベクトルを生成し、
前記第２のセンサデータ特徴ベクトルを合成して第１の複軸センサデータ特徴ベクトルを生成する第１の複軸センサデータ特徴ベクトル生成手段を備え、
前記第１の特定信号パターン特徴マップ生成手段は、前記第１の複軸センサデータ特徴ベクトルを自己組織化アルゴリズムにより学習させることで自己組織化された第２の特定信号パターン特徴マップを生成し、
前記蓄積手段は、前記第２の特定信号パターン特徴マップの各ノードにおける特徴ノードベクトルに基づいて、前記擬似複軸センサデータの特徴ベクトルを入力した時に発火する特徴ノードベクトルとの対応関係を第２の蓄積結果として蓄積し、
前記第１の特定信号パターン発火確率マップ生成手段は、前記第２の蓄積結果に基づいて、前記第２の特定信号パターン特徴マップの各ノードごとに特定信号パターンが発火する確率を算出し第２の特定信号パターン発火確率マップを生成することを特徴とする特定信号パターン検出システム。
【請求項２】
複軸センサから得られる時系列の複軸センサデータを予め用意された時間窓に基づいて時間分割複軸センサデータとして分割する第２のデータ分割手段と、
前記時間分割複軸センサデータのノイズ除去ならびにデータを正規化することで第２の正規化センサデータを生成する第２の正規化センサデータ生成手段と、
前記第２の正規化センサデータを基底関数展開することで時空間解像度レベルごとに第３の時空間分解センサデータを生成する第２の時空間分解センサデータ生成手段と、
前記第３の時空間分解センサデータの振幅情報に基づいて時空間解像度レベルごとに第３の時空間分解センサデータ情報エントロピーを算出する第２のエントロピー算出手段と、
前記第３の時空間分解センサデータ情報エントロピーの総和が最小になるように予め用意された評価基準に基づいて第３の基底時空間分解センサデータを選択する第２の基底時空間分解センサデータ選択手段と、
前記第３の基底時空間分解センサデータの情報エントロピーに基づいて第３のセンサデータ特徴ベクトルを生成する第２のセンサデータ特徴ベクトル生成手段と、
前記第３のセンサデータ特徴ベクトルを合成して第２の複軸センサデータ特徴ベクトルを生成する第２の複軸センサデータ特徴ベクトル生成手段と、
前記第２の複軸センサデータ特徴ベクトルを事前に得られた複数の特定信号パターン特徴マップに入力して得られる類似度と、前記第２の複軸センサデータ特徴ベクトルを事前に得られた複数の特定信号パターン発火確率マップに入力して得られる発火した特徴マップノードベクトルの発火確率とに基づき、前記発火した特徴マップノードベクトルから特定信号パターン特徴ベクトルを選択する特定信号パターン特徴ベクトル選択手段と、
を含み、
前記事前に得られた複数の特定信号パターン特徴マップは、請求項１に記載の特定信号パターン検出システムにおける前記第１の特定信号パターン特徴マップ及び前記第２の特定信号パターン特徴マップを含み、
前記事前に得られた複数の特定信号パターン発火確率マップは、請求項１に記載の特定信号パターン検出システムにおける前記第１の特定信号パターン発火確率マップ及び前記第２の特定信号パターン発火確率マップを含むことを特徴とする特定信号パターン検出システム。
【請求項３】
複軸センサから得られる時系列の複軸センサデータにおいて検出したい複数の特定信号パターンを辞書データとして予め獲得するとともに、その獲得された特定信号パターンを予め用意された時間窓に基づいて第１の時間分割センサデータとして分割する第１のデータ分割ステップと、
前記第１の時間分割センサデータのノイズ除去ならびにデータを正規化することで第１の正規化センサデータを生成する第１の正規化センサデータ生成ステップと、
前記第１の正規化センサデータを基底関数展開することで時空間解像度レベルごとに第１の時空間分解センサデータを生成する第１の時空間分解センサデータ生成ステップと、
前記第１の時空間分解センサデータの振幅情報に基づいて時空間解像度レベルごとに第１の時空間分解センサデータ情報エントロピーを算出する第１のエントロピー算出ステップと、
前記第１の時空間分解センサデータ情報エントロピーの総和が最小になるように予め用意された評価基準に基づいて第１の基底時空間分解センサデータを選択する第１の基底時空間分解センサデータ選択ステップと、
前記第１の基底時空間分解センサデータの情報エントロピーに基づいて、前記辞書データの特徴ベクトルである第１のセンサデータ特徴ベクトルを生成する第１のセンサデータ特徴ベクトル生成ステップと、
前記第１のセンサデータ特徴ベクトルを自己組織化アルゴリズムにより学習させることで自己組織化された第１の特定信号パターン特徴マップを生成する第１の特定信号パターン特徴マップ生成ステップと、
前記第１の特定信号パターン特徴マップの各ノードにおける特徴ノードベクトルに基づいて、前記辞書データの特徴ベクトルを入力した時に発火する特徴ノードベクトルとの対応関係を第１の蓄積結果として蓄積する第１の蓄積ステップと、
前記第１の蓄積結果に基づいて、前記第１の特定信号パターン特徴マップの各ノードごとに特定信号パターンが発火する確率を算出し第１の特定信号パターン発火確率マップを生成する第１の特定信号パターン発火確率マップ生成ステップと、
前記第１の特定信号パターン発火確率マップのうち予め設定された閾値以上の確率を有するノードを基底特徴ノードベクトルとして抽出する基底特徴ノードベクトル抽出ステップと、
前記基底特徴ノードベクトルを発火させる基となった複軸センサデータ特徴ベクトルのうち、基底特徴ノードベクトルとの類似度が予め設定された閾値以上のものを基底複軸センサ特徴ベクトルとして抽出する基底複軸センサ特徴ベクトル抽出ステップと、
前記基底複軸センサ特徴ベクトルの基となった複軸の基底センサデータを抽出する基底センサデータ抽出ステップと、
前記複軸の基底センサデータを予め設定された刻み幅の角度で各軸回りに回転させることで擬似複軸センサデータを生成する擬似複軸センサデータ生成ステップと、
前記疑似複軸センサデータを基底関数展開することで時空間解像度レベルごとに第２の時空間分解センサデータを生成する第２の時空間分解センサデータ生成ステップと、
前記第２の時空間分解センサデータの振幅情報に基づいて時空間解像度レベルごとに第２の時空間分解センサデータ情報エントロピーを算出する第２のエントロピー算出ステップと、
前記第２の時空間分解センサデータ情報エントロピーの総和が最小になるように予め用意された評価基準に基づいて第２の基底時空間分解センサデータを選択する第２の基底時空間分解センサデータ選択ステップと、
前記第２の基底時空間分解センサデータの情報エントロピーに基づいて第２のセンサデータ特徴ベクトルを生成する第２のセンサデータ特徴ベクトル生成ステップと、
前記第２のセンサデータ特徴ベクトルを合成して第１の複軸センサデータ特徴ベクトルを生成する第１の複軸センサデータ特徴ベクトル生成ステップと、
前記第１の複軸センサデータ特徴ベクトルを自己組織化アルゴリズムにより学習させることで自己組織化された第２の特定信号パターン特徴マップを生成する第２の特定信号パターン特徴マップ生成ステップと、
前記第２の特定信号パターン特徴マップの各ノードにおける特徴ノードベクトルに基づいて、前記擬似複軸センサデータの特徴ベクトルを入力した時に発火する特徴ノードベクトルとの対応関係を第２の蓄積結果として蓄積する第２の蓄積ステップと、
前記第２の蓄積結果に基づいて、前記第２の特定信号パターン特徴マップの各ノードごとに特定信号パターンが発火する確率を算出し第２の特定信号パターン発火確率マップを生成する第２の特定信号パターン発火確率マップ生成ステップと、
を含むことを特徴とする特定信号パターン検出方法。
【請求項４】
複軸センサから得られる時系列の複軸センサデータを予め用意された時間窓に基づいて時間分割複軸センサデータとして分割する第２のデータ分割ステップと、
前記時間分割複軸センサデータのノイズ除去ならびにデータを正規化することで第２の正規化センサデータを生成する第２の正規化センサデータ生成ステップと、
前記第２の正規化センサデータを基底関数展開することで時空間解像度レベルごとに第３の時空間分解センサデータを生成する第２の時空間分解センサデータ生成ステップと、
前記第３の時空間分解センサデータの振幅情報に基づいて時空間解像度レベルごとに第３の時空間分解センサデータ情報エントロピーを算出する第２のエントロピー算出ステップと、
前記第３の時空間分解センサデータ情報エントロピーの総和が最小になるように予め用意された評価基準に基づいて第３の基底時空間分解センサデータを選択する第２の基底時空間分解センサデータ選択ステップと、
前記第３の基底時空間分解センサデータの情報エントロピーに基づいて第３のセンサデータ特徴ベクトルを生成する第２のセンサデータ特徴ベクトル生成ステップと、
前記第３のセンサデータ特徴ベクトルを合成して第２の複軸センサデータ特徴ベクトルを生成する第２の複軸センサデータ特徴ベクトル生成ステップと、
前記第２の複軸センサデータ特徴ベクトルを事前に得られた複数の特定信号パターン特徴マップに入力して得られる類似度と、前記第２の複軸センサデータ特徴ベクトルを事前に得られた複数の特定信号パターン発火確率マップに入力して得られる発火した特徴マップノードベクトルの発火確率とに基づき、前記発火した特徴マップノードベクトルから特定信号パターン特徴ベクトルを選択する特定信号パターン特徴ベクトル選択ステップと、
を含み、
前記事前に得られた複数の特定信号パターン特徴マップは、請求項３に記載の特定信号パターン検出方法における前記第１の特定信号パターン特徴マップ及び前記第２の特定信号パターン特徴マップを含み、
前記事前に得られた複数の特定信号パターン発火確率マップは、請求項３に記載の特定信号パターン検出方法における前記第１の特定信号パターン発火確率マップ及び前記第２の特定信号パターン発火確率マップを含むことを特徴とする特定信号パターン検出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は複軸センサからの特定信号パターン検出システム、複軸センサからの特定信号パターン検出方法に関し、特にユーザのコンテキストを獲得するための特定信号パターン検出システム、特定信号パターン検出方法に関する。
【背景技術】
【０００２】
実世界指向型サービスを実現する上でモバイル環境下での個人の状況（コンテキスト）を理解することは重要である。従来、携帯端末を持ち歩くようなモバイル環境下でユーザのコンテキストを獲得に関する研究は、例えば、非特許文献１に記載されている。
非特許文献１においては、自己組織化特徴マップのノードについて区分けを行い（すなわちクラスタリングを行い）、その区分けを固定している。そして、それに従って処理を行うことで、コンテキスト情報を取得している。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】A. Krause et. al, "Unsupervised, Dynamic Identification of Physiological and Activity Context in Wearable Computing",[online]、[平成２４年５月３日検索]、インターネット＜ＵＲＬ：http://wearables.cc.gatech.edu/paper_of_week/krause-iswc2003-unsup-activity.pdf＞
【発明の概要】
【発明が解決しようとする課題】
【０００４】
上述した非特許文献１に記載されている技術には、以下のような制約があった。
（１）獲得したいコンテキストに応じてセンサを特定の位置に固定することが必要
（２）特定の軸のセンサ情報のみ、もしくは、各軸のセンサデータからセンサデータ全体の大きさにしたものを利用しているため、抽出可能な特定信号パターンが限定
（３）辞書データを部分空間法などの線形手法により特定信号パターンの基底データを選択していたため、ゆらぎのあるデータへ適用した場合は基底データの抽出精度に影響を受け易い。
本発明は上述した従来技術の課題を解決するためになされたものであり、その目的は複軸センサの設置位置に依らず、また、センサデータのゆらぎにロバストであるユーザのコンテキストに相当する特定信号パターンを獲得する特定信号パターン検出システム、特定信号パターン検出方法を提供することである。
【課題を解決するための手段】
【０００５】
本発明の一態様による特定信号パターン検出システムは、
複軸センサから得られる時系列の複軸センサデータにおいて検出したい複数の特定信号パターンを辞書データとして予め獲得するとともに、その獲得された特定信号パターンを予め用意された時間窓に基づいて第１の時間分割センサデータとして分割する第１のデータ分割手段と、
前記第１の時間分割センサデータのノイズ除去ならびにデータを正規化することで第１の正規化センサデータを生成する第１の正規化センサデータ生成手段と、
前記第１の正規化センサデータを基底関数展開することで時空間解像度レベルごとに第１の時空間分解センサデータを生成する第１の時空間分解センサデータ生成手段と、
前記第１の時空間分解センサデータの振幅情報に基づいて時空間解像度レベルごとに第１の時空間分解センサデータ情報エントロピーを算出する第１のエントロピー算出手段と、
前記第１の時空間分解センサデータ情報エントロピーの総和が最小になるように予め用意された評価基準に基づいて第１の基底時空間分解センサデータを選択する第１の基底時空間分解センサデータ選択手段と、
前記第１の基底時空間分解センサデータの情報エントロピーに基づいて、前記辞書データの特徴ベクトルである第１のセンサデータ特徴ベクトルを生成する第１のセンサデータ特徴ベクトル生成手段と、
前記第１のセンサデータ特徴ベクトルを自己組織化アルゴリズムにより学習させることで自己組織化された第１の特定信号パターン特徴マップを生成する第１の特定信号パターン特徴マップ生成手段と、
前記第１の特定信号パターン特徴マップの各ノードにおける特徴ノードベクトルに基づいて、前記辞書データの特徴ベクトルを入力した時に発火する特徴ノードベクトルとの対応関係を第１の蓄積結果として蓄積する蓄積手段と、
前記第１の蓄積結果に基づいて、前記第１の特定信号パターン特徴マップの各ノードごとに特定信号パターンが発火する確率を算出し第１の特定信号パターン発火確率マップを生成する第１の特定信号パターン発火確率マップ生成手段と、
前記第１の特定信号パターン発火確率マップのうち予め設定された閾値以上の確率を有するノードを基底特徴ノードベクトルとして抽出する基底特徴ノードベクトル抽出手段と、
前記基底特徴ノードベクトルを発火させる基となった複軸センサデータ特徴ベクトルのうち、基底特徴ノードベクトルとの類似度が予め設定された閾値以上のものを基底複軸センサ特徴ベクトルとして抽出する基底複軸センサ特徴ベクトル抽出手段と、
前記基底複軸センサ特徴ベクトルの基となった複軸の基底センサデータを抽出する基底センサデータ抽出手段と、
前記複軸の基底センサデータを予め設定された刻み幅の角度で各軸回りに回転させることで擬似複軸センサデータを生成する擬似複軸センサデータ生成手段と、
を含み、
さらに
前記第１の時空間分解センサデータ生成手段は、前記疑似複軸センサデータを基底関数展開することで時空間解像度レベルごとに第２の時空間分解センサデータを生成し、
前記第１のエントロピー算出手段は、前記第２の時空間分解センサデータの振幅情報に基づいて時空間解像度レベルごとに第２の時空間分解センサデータ情報エントロピーを算出し、
前記第１の基底時空間分解センサデータ選択手段は、前記第２の時空間分解センサデータ情報エントロピーの総和が最小になるように予め用意された評価基準に基づいて第２の基底時空間分解センサデータを選択し、
前記第１のセンサデータ特徴ベクトル生成手段は、前記第２の基底時空間分解センサデータの情報エントロピーに基づいて第２のセンサデータ特徴ベクトルを生成し、
前記第２のセンサデータ特徴ベクトルを合成して第１の複軸センサデータ特徴ベクトルを生成する第１の複軸センサデータ特徴ベクトル生成手段を備え、
前記第１の特定信号パターン特徴マップ生成手段は、前記第１の複軸センサデータ特徴ベクトルを自己組織化アルゴリズムにより学習させることで自己組織化された第２の特定信号パターン特徴マップを生成し、
前記蓄積手段は、前記第２の特定信号パターン特徴マップの各ノードにおける特徴ノードベクトルに基づいて、前記擬似複軸センサデータの特徴ベクトルを入力した時に発火する特徴ノードベクトルとの対応関係を第２の蓄積結果として蓄積し、
前記第１の特定信号パターン発火確率マップ生成手段は、前記第２の蓄積結果に基づいて、前記第２の特定信号パターン特徴マップの各ノードごとに特定信号パターンが発火する確率を算出し第２の特定信号パターン発火確率マップを生成することを特徴とする。
このように構成すれば、より質の良いデータを抽出することができる。
【０００６】
本発明の他の態様による特定信号パターン検出システムは、
複軸センサから得られる時系列の複軸センサデータを予め用意された時間窓に基づいて時間分割複軸センサデータとして分割する第２のデータ分割手段と、
前記時間分割複軸センサデータのノイズ除去ならびにデータを正規化することで第２の正規化センサデータを生成する第２の正規化センサデータ生成手段と、
前記第２の正規化センサデータを基底関数展開することで時空間解像度レベルごとに第３の時空間分解センサデータを生成する第２の時空間分解センサデータ生成手段と、
前記第３の時空間分解センサデータの振幅情報に基づいて時空間解像度レベルごとに第３の時空間分解センサデータ情報エントロピーを算出する第２のエントロピー算出手段と、
前記第３の時空間分解センサデータ情報エントロピーの総和が最小になるように予め用意された評価基準に基づいて第３の基底時空間分解センサデータを選択する第２の基底時空間分解センサデータ選択手段と、
前記第３の基底時空間分解センサデータの情報エントロピーに基づいて第３のセンサデータ特徴ベクトルを生成する第２のセンサデータ特徴ベクトル生成手段と、
前記第３のセンサデータ特徴ベクトルを合成して第２の複軸センサデータ特徴ベクトルを生成する第２の複軸センサデータ特徴ベクトル生成手段と、
前記第２の複軸センサデータ特徴ベクトルを事前に得られた複数の特定信号パターン特徴マップに入力して得られる類似度と、前記第２の複軸センサデータ特徴ベクトルを事前に得られた複数の特定信号パターン発火確率マップに入力して得られる発火した特徴マップノードベクトルの発火確率とに基づき、前記発火した特徴マップノードベクトルから特定信号パターン特徴ベクトルを選択する特定信号パターン特徴ベクトル選択手段と、
を含み、
前記事前に得られた複数の特定信号パターン特徴マップは、前記一態様による特定信号パターン検出システムにおける前記第１の特定信号パターン特徴マップ及び前記第２の特定信号パターン特徴マップを含み、
前記事前に得られた複数の特定信号パターン発火確率マップは、前記一態様による特定信号パターン検出システムにおける前記第１の特定信号パターン発火確率マップ及び前記第２の特定信号パターン発火確率マップを含むことを特徴とする。
このように構成すれば、センサが任意の姿勢で設置されても正確なコンテキスト情報を取得することができる。
【０００９】
本発明の一態様による特定信号パターン検出方法は、
複軸センサから得られる時系列の複軸センサデータにおいて検出したい複数の特定信号パターンを辞書データとして予め獲得するとともに、その獲得された特定信号パターンを予め用意された時間窓に基づいて第１の時間分割センサデータとして分割する第１のデータ分割ステップと、
前記第１の時間分割センサデータのノイズ除去ならびにデータを正規化することで第１の正規化センサデータを生成する第１の正規化センサデータ生成ステップと、
前記第１の正規化センサデータを基底関数展開することで時空間解像度レベルごとに第１の時空間分解センサデータを生成する第１の時空間分解センサデータ生成ステップと、
前記第１の時空間分解センサデータの振幅情報に基づいて時空間解像度レベルごとに第１の時空間分解センサデータ情報エントロピーを算出する第１のエントロピー算出ステップと、
前記第１の時空間分解センサデータ情報エントロピーの総和が最小になるように予め用意された評価基準に基づいて第１の基底時空間分解センサデータを選択する第１の基底時空間分解センサデータ選択ステップと、
前記第１の基底時空間分解センサデータの情報エントロピーに基づいて、前記辞書データの特徴ベクトルである第１のセンサデータ特徴ベクトルを生成する第１のセンサデータ特徴ベクトル生成ステップと、
前記第１のセンサデータ特徴ベクトルを自己組織化アルゴリズムにより学習させることで自己組織化された第１の特定信号パターン特徴マップを生成する第１の特定信号パターン特徴マップ生成ステップと、
前記第１の特定信号パターン特徴マップの各ノードにおける特徴ノードベクトルに基づいて、前記辞書データの特徴ベクトルを入力した時に発火する特徴ノードベクトルとの対応関係を第１の蓄積結果として蓄積する第１の蓄積ステップと、
前記第１の蓄積結果に基づいて、前記第１の特定信号パターン特徴マップの各ノードごとに特定信号パターンが発火する確率を算出し第１の特定信号パターン発火確率マップを生成する第１の特定信号パターン発火確率マップ生成ステップと、
前記第１の特定信号パターン発火確率マップのうち予め設定された閾値以上の確率を有するノードを基底特徴ノードベクトルとして抽出する基底特徴ノードベクトル抽出ステップと、
前記基底特徴ノードベクトルを発火させる基となった複軸センサデータ特徴ベクトルのうち、基底特徴ノードベクトルとの類似度が予め設定された閾値以上のものを基底複軸センサ特徴ベクトルとして抽出する基底複軸センサ特徴ベクトル抽出ステップと、
前記基底複軸センサ特徴ベクトルの基となった複軸の基底センサデータを抽出する基底センサデータ抽出ステップと、
前記複軸の基底センサデータを予め設定された刻み幅の角度で各軸回りに回転させることで擬似複軸センサデータを生成する擬似複軸センサデータ生成ステップと、
前記疑似複軸センサデータを基底関数展開することで時空間解像度レベルごとに第２の時空間分解センサデータを生成する第２の時空間分解センサデータ生成ステップと、
前記第２の時空間分解センサデータの振幅情報に基づいて時空間解像度レベルごとに第２の時空間分解センサデータ情報エントロピーを算出する第２のエントロピー算出ステップと、
前記第２の時空間分解センサデータ情報エントロピーの総和が最小になるように予め用意された評価基準に基づいて第２の基底時空間分解センサデータを選択する第２の基底時空間分解センサデータ選択ステップと、
前記第２の基底時空間分解センサデータの情報エントロピーに基づいて第２のセンサデータ特徴ベクトルを生成する第２のセンサデータ特徴ベクトル生成ステップと、
前記第２のセンサデータ特徴ベクトルを合成して第１の複軸センサデータ特徴ベクトルを生成する第１の複軸センサデータ特徴ベクトル生成ステップと、
前記第１の複軸センサデータ特徴ベクトルを自己組織化アルゴリズムにより学習させることで自己組織化された第２の特定信号パターン特徴マップを生成する第２の特定信号パターン特徴マップ生成ステップと、
前記第２の特定信号パターン特徴マップの各ノードにおける特徴ノードベクトルに基づいて、前記擬似複軸センサデータの特徴ベクトルを入力した時に発火する特徴ノードベクトルとの対応関係を第２の蓄積結果として蓄積する第２の蓄積ステップと、
前記第２の蓄積結果に基づいて、前記第２の特定信号パターン特徴マップの各ノードごとに特定信号パターンが発火する確率を算出し第２の特定信号パターン発火確率マップを生成する第２の特定信号パターン発火確率マップ生成ステップと、を含むことを特徴とする。
このような方法によれば、より質の良いデータを抽出することができる。
【００１０】
本発明の他の態様による特定信号パターン検出方法は、
複軸センサから得られる時系列の複軸センサデータを予め用意された時間窓に基づいて時間分割複軸センサデータとして分割する第２のデータ分割ステップと、
前記時間分割複軸センサデータのノイズ除去ならびにデータを正規化することで第２の正規化センサデータを生成する第２の正規化センサデータ生成ステップと、
前記第２の正規化センサデータを基底関数展開することで時空間解像度レベルごとに第３の時空間分解センサデータを生成する第２の時空間分解センサデータ生成ステップと、
前記第３の時空間分解センサデータの振幅情報に基づいて時空間解像度レベルごとに第３の時空間分解センサデータ情報エントロピーを算出する第２のエントロピー算出ステップと、
前記第３の時空間分解センサデータ情報エントロピーの総和が最小になるように予め用意された評価基準に基づいて第３の基底時空間分解センサデータを選択する第２の基底時空間分解センサデータ選択ステップと、
前記第３の基底時空間分解センサデータの情報エントロピーに基づいて第３のセンサデータ特徴ベクトルを生成する第２のセンサデータ特徴ベクトル生成ステップと、
前記第３のセンサデータ特徴ベクトルを合成して第２の複軸センサデータ特徴ベクトルを生成する第２の複軸センサデータ特徴ベクトル生成ステップと、
前記第２の複軸センサデータ特徴ベクトルを事前に得られた複数の特定信号パターン特徴マップに入力して得られる類似度と、前記第２の複軸センサデータ特徴ベクトルを事前に得られた複数の特定信号パターン発火確率マップに入力して得られる発火した特徴マップノードベクトルの発火確率とに基づき、前記発火した特徴マップノードベクトルから特定信号パターン特徴ベクトルを選択する特定信号パターン特徴ベクトル選択ステップと、
を含み、
前記事前に得られた複数の特定信号パターン特徴マップは、前記一態様による特定信号パターン検出方法における前記第１の特定信号パターン特徴マップ及び前記第２の特定信号パターン特徴マップを含み、
前記事前に得られた複数の特定信号パターン発火確率マップは、前記一態様による特定信号パターン検出方法における前記第１の特定信号パターン発火確率マップ及び前記第２の特定信号パターン発火確率マップを含むことを特徴とする。
このような方法によれば、センサが任意の姿勢で設置されても正確なコンテキスト情報を取得することができる。
【００１３】
要するに、本発明では、複軸センサから出力されるデータを取得し、取得データに基づいて、前記複軸センサを含む機器を携行している者のコンテキストを推定するための特定信号パターン特徴マップを作成する場合において、複軸ごとにセンシング可能なセンサデバイス(複軸センサ)の設置位置に依存せずに、その複軸センサから得られる時系列の複軸センサデータにおける特徴的な信号パターン（特定信号パターン）を検出するシステムにおいて、複軸センサから得られる複軸センサデータの中で検出したい特定信号パターンを辞書データとして獲得し、その辞書データを自己組織化アルゴリズムで学習させて得られる特定信号パターン特徴マップならびに特定信号パターン発火確率マップから辞書データの特定信号パターン基底センサデータを抽出している。そして、抽出された特定信号パターン基底センサデータを各軸回りに回転させることで生成した擬似データを自己組織化アルゴリズムで学習させることで特定信号パターン特徴マップならびに特定信号パターン発火確率マップを得ている。さらに、上記で得られた、特定信号パターン特徴マップならびに特定信号パターン発火確率マップを用いて、複軸センサから得られる時系列のセンサデータから特定信号パターンを抽出している。
【発明の効果】
【００１４】
本発明によれば、獲得したいコンテキストごとにモジュールを用意して学習させる等のチューニングをすることなくユーザのコンテキストに相当する特定信号パターンが獲得できるという効果がある。
また、複軸センサを特定の位置に固定する必要がないので、専用の複軸センサを用いる必要がなく、例えば携帯電話機に搭載してある複軸センサを用いてプレゼンス情報を取得することができるという効果がある。
【図面の簡単な説明】
【００１５】
【図１】本発明の実施形態による、複軸センサからの特定信号パターン検出システムの構成を示す図である。
【図２】図１中のユーザＢが使用する携帯端末の内部構成例を示す図である。
【図３】図１中のプレゼンス提供サーバの内部構成例を示す図である。
【図４】図１中のユーザＡが使用する携帯端末３００の内部構成例を示す図である。
【図５】本発明の実施の形態に係る複軸センサからの特定信号パターン検出方法の概略を示す図である。
【図６】図５中の段階（Ｃ−１）の詳細な処理内容を示す図である。
【図７】階段降り動作を行った場合の自己組織化特徴マップの例を示す図である。
【図８】階段昇り動作を行った場合の自己組織化特徴マップの例を示す図である。
【図９】早歩き動作を行った場合の自己組織化特徴マップの例を示す図である。
【図１０】電車内歩行動作を行った場合の自己組織化特徴マップの例を示す図である。
【図１１】電車着座動作を行った場合の自己組織化特徴マップの例を示す図である。
【図１２】走行動作を行った場合の自己組織化特徴マップの例を示す図である。
【図１３】平地歩行動作を行った場合の自己組織化特徴マップの例を示す図である。
【図１４】図５中の段階（Ｃ−２）の詳細な処理内容を示す図である。
【図１５】図５中の段階（Ｃ−３）の詳細な処理内容を示す図である。
【図１６】特定信号パターン検出システムの実施例の構成例を示す図である。
【発明を実施するための形態】
【００１６】
以下、本発明の実施の形態を、図面を参照して説明する。なお、以下の説明において参照する各図では、他の図と同等部分は同一符号によって示されている。
以下、本発明の実施の形態を、図面を参照して説明する。以下の説明において参照する各図では、他の図と同等部分は同一符号によって示されている。
なお、以下の説明において、（→）が付加されている符号はそれがベクトルであることを示している。また、→が上に付加されている符号もそれがベクトルであることを示している。
【００１７】
（基本構成）
本発明の実施形態による、複軸センサからの特定信号パターン検出システムの構成について図１を参照して説明する。同図は、特定信号パターン検出システムの基本的な構成例を示すものである。
同図を参照すると、本システムは、ユーザＢの状況(コンテキスト)を獲得するために必要な複軸センサデータを検出できる加速度センサを搭載した携帯端末２００と、この携帯端末２００に搭載された加速度センサから検出された複軸センサデータからユーザＢの状況をプレゼンス情報として提供する装置であるプレゼンス提供サーバ１００と、このプレゼンス提供サーバ１００で提供されたプレゼンス情報をユーザＡへ提示するための携帯端末３００とを含んで構成されている。なお、同図中のモバイルネットワークＮＷは、これらの間で信号授受を実現するためのネットワークであり、図示せぬ無線基地局装置、無線基地局制御装置等によって実現される。
【００１８】
図１中のユーザＢの携帯端末２００の内部構成例が図２に示されている。同図を参照すると、携帯端末２００は、加速度センサから複軸センサデータを検出するセンサデータ獲得機能部２１と、検出された複軸センサデータをプレゼンス提供サーバへ送信するセンサデータ送信機能部２２とが搭載されている。この携帯端末２００は、図示せぬＣＰＵ、通信部、メモリ、入出力操作部、表示部等、一般的な携帯電話端末としてのハードウェア構成を有しており、これらの構成によって上記の機能を実現している。
【００１９】
図１中のプレゼンス提供サーバ１００の内部構成例が図３に示されている。同図を参照すると、プレゼンス提供サーバ１００は、上記携帯端末２００より送信された複軸センサデータを受信するセンサデータ受信機能部１１と、このセンサデータ受信機能部１１によって受信された複軸センサデータから信号パターン特徴ベクトルを抽出するセンサデータ解析機能部１２と、事前に検出対象となる特定信号パターン特徴ベクトルを抽出し、そのデータを特定信号パターンデータベース１０１へ格納する特定信号パターン獲得機能部１３と、センサデータ解析機能部１２で抽出された信号パターン特徴ベクトルとデータベース１０１に格納されている特定信号パターン特徴ベクトルと照合するセンサデータ照合機能部１４と、センサデータ照合機能部１４で照合された特定信号パターン特徴ベクトルをプレゼンスデータベース１０２に基づいてプレゼンス情報に変換するプレゼンス情報変換機能部１５と、変換されたプレゼンス情報をユーザＡのリクエストに応じて提供するプレゼンス情報送信機能部１６とを含んで構成されている。このプレゼンス提供サーバ１００は、図示せぬＣＰＵ、通信部、メモリ、入出力部、表示部等、一般的なサーバ装置としてのハードウェア構成を有しており、これらの構成によって上記の機能を実現している。
【００２０】
図１中のユーザＡの携帯端末３００の内部構成例が図４に示されている。同図を参照すると、携帯端末３００は、プレゼンス提供サーバ１００から送信されたユーザＢのプレゼンス情報を受信し表示するプレゼンス情報受信機能部２３を含んで構成されている。この携帯端末３００は、図示せぬＣＰＵ、通信部、メモリ、入出力操作部、表示部等、一般的な携帯電話端末としてのハードウェア構成を有しており、これらの構成によって上記の機能を実現している。
【００２１】
（複軸センサからの特定信号パターン検出方法）
図１〜図４に示されている構成からなるシステムによって実現される、複軸センサからの特定信号パターン検出方法の具体的な処理の流れについて、以下説明する。
図５に示されているように、本検出方法による処理は、オフライン処理で学習を行う段階（Ｃ−１）及び段階（Ｃ−２）と、オンライン処理で実測を行う段階（Ｃ−３）の３つの段階に大別することができる。なお、経年変化に対応するため、追加学習を行う処理を
追加してもよい。この処理を追加すれば、ユーザＢの年齢増加にも対応できる。
【００２２】
（Ｃ−１）
この段階は、質の良いデータを抽出する段階である。この段階では、辞書データＤ１を入力とし、前処理（Ｓ１１）、特徴抽出処理（Ｓ１２）、自己組織化学習（Ｓ１３）、学習結果解析処理（Ｓ１４）、基底センサデータ抽出処理（Ｓ１５）を経て、特定信号パターン特徴マップ及び特定信号パターン発火確率マップからなるデータＤ２と、特定信号パターン基底センサデータＤ３とを導出する。
【００２３】
（Ｃ−２）
この段階は、上記の段階（Ｃ−１）で得たデータを基に、より質の良い辞書データを得る段階である。この段階では、上記段階（Ｃ−１）で得られた特定信号パターン基底センサデータＤ３を入力とし、擬似データ生成処理（Ｓ２１）、特徴抽出処理（Ｓ２２）、自己組織化学習（Ｓ２３）、学習結果解析処理（Ｓ２４）を経て、特定信号パターン特徴マップ及び特定信号パターン発火確率マップからなるデータＤ４を導出する。
【００２４】
（Ｃ−３）
この段階は、上記の段階（Ｃ−２）で得たデータを基に、実際の測定を行う段階である。この段階では、実測データＤ５を入力とし、前処理（Ｓ３１）、特徴抽出処理（Ｓ３２）、特定信号パターン検出処理（Ｓ３３）を経て、特定信号パターン分析結果データＤ６を導出する。
【００２５】
以下、各段階について、図６〜図１５を参照して詳細に説明する。
（質の良いデータの抽出）
図６は、図５中の段階（Ｃ−１）の詳細な処理内容を示す図である。
図５中の段階（Ｃ−１）中の前処理（Ｓ１１）は図６中の処理（Ｃ−１−２）及び（Ｃ−１−３）に対応する。図５中の段階（Ｃ−１）中の特徴抽出処理（Ｓ１２）は図６中の処理（Ｃ−１−４）から（Ｃ−１−７）までに対応する。図５中の段階（Ｃ−１）中の自己組織化学習処理（Ｓ１３）は図６中の処理（Ｃ−１−８）に対応する。図５中の段階（
Ｃ−１）中の学習結果解析処理（Ｓ１４）は図６中の処理（Ｃ−１−９）及び（Ｃ−１−１０）に対応する。図５中の段階（Ｃ−１）中の基底センサデータ抽出処理（Ｓ１５）は図６中の処理（Ｃ−１−１１）に対応する。以下、これらの処理の内容について説明する。
【００２６】
処理（Ｃ−１−１）：ユーザＢは複数軸センサデータが獲得できる加速度センサを搭載してある携帯端末を携行している。逐次、ユーザＢの状況（コンテキスト）に応じて変動する複数軸センサデータ（加速度ベクトル）が獲得され、パケット通信でプレゼンス提供サーバへ送信される。
処理（Ｃ−１−２）：プレゼンス提供サーバでは、ユーザＢの携帯端末から送信された加速度ベクトルの時系列データを逐次受信する。時系列データとして受信される加速度ベクトル（→）ａ0（ｔ）をデータ分割する。手法は限定されないが、ここでは、複数の時間窓
【００２７】
【数１】
0005297547.tif 000002
を用意してデータ分割する。
【００２８】
処理（Ｃ−１−３）：その後、データ補正処理を実施する。データ補正処理方法は加速度センサの用途によって異なるが、ここでは、Gaussianフィルターによるノイズ除去処理と、時間窓｛ｔk｝での加速度の正規化を行う。
続いて、センサデータ解析から特定信号パターンとの照合について述べる。
任意の姿勢で設置されたセンサから観測される時刻ｔでの実測データ加速度ベクトルを
（→）ａ0（ｔ）＝（ａ0x（ｔ）、ａ0y（ｔ）、ａ0z（ｔ））  …（２）
理想的な姿勢で設置されたセンサから観測される時刻ｔにおける辞書データとなる加速度ベクトル（重力加速度も含む）を
（→）ａD（ｔ）＝（ａDX（ｔ）、ａDY（ｔ）、ａDZ（ｔ））  …（３）
とおく。
【００２９】
はじめに、辞書データとなる加速度ベクトル（→）ａD（ｔ）から特定信号パターン特徴マップを作成する方法について説明する。
処理（Ｃ−１−４）：まず、加速度ベクトル（→）ａD（ｔ）の特徴を解析するために基底関数展開することで時空間解像度レベルごとに時空間分解センサデータに分解する。ここでは局在波による関数展開であるMother WaveletをHarr関数とするWavelet Packet分解を利用する。上記の式（３）の両辺にWavelet Packet分解を施すと解像度レベルｐでの加速度ベクトル（→）ａD（ｔ）は以下のように表せる。
【００３０】
【数２】
0005297547.tif 000003
【００３１】
ここで、（→）ｕ(p、q)D（ｔ）は、解像度レベルｐのｑ番目の分解波形（基底関数）を表す。
処理（Ｃ−１−５）：次に、これら分解波形である基底関数
（→）ｕ(p、q)D（ｔ）＝
（（→）ｕ(p、q)x（ｔ）、（→）ｕ(p、q)y（ｔ）、（→）ｕ(p、q)z（ｔ））
のうち、加速度ベクトル（→）ａD（ｔ）を表現するのに最良基底関数となる基底関数を選択する。そこで、以下の情報エントロピーコスト関数を定義する。
【００３２】
【数３】
0005297547.tif 000004
【００３３】
式（５）の左辺を最小にするｐとｑとの組み合わせを選択する。これにより、情報エントロピーの増加を防ぎ（冗長な表現ではなく）、波形を正確に分解することができる。さらに、選択された基底関数（→）ｕ(p、q)D_BN（ｔ）の情報エントロピー
【００３４】
【数４】
0005297547.tif 000005
【００３５】
がもっとも大きいものから選び、それらを順次加えていった和が全情報エントロピー（情報エントロピーの総和）に対して閾値ＴSoEに達するまでの基底関数を選択する。
処理（Ｃ−１−６）：そして、最終的に得られた基底関数（→）ｕ(p、q)D_F（ｔ）に基づいて辞書データの特徴ベクトル（→）ＸDを以下のように定義する。
【００３６】
【数５】
0005297547.tif 000006
ここで、
【００３７】
【数６】
0005297547.tif 000007
であり、直流成分Ｘ(i)D_BN（０）の項は用いない。
処理（Ｃ−１−７）：また、このとき特徴ベクトルの分布統計量として、下記のモーメントを算出する。ここでは、分布統計量ベクトルを下記のベクトル（→）Ｍ（（→）Ｘ0）として表す。
（→）Ｍ（（→）Ｘ0）＝（Ｍ(1)（（→）Ｘ0）、Ｍ(2)（（→）Ｘ0）、…、Ｍ(n)（（→
）Ｘ0））…（９）
【００３８】
【数７】
0005297547.tif 000008
ただし、
【００３９】
【数８】
0005297547.tif 000009
【００４０】
そして、辞書データとなる全ての加速度ベクトル（→）ａD（ｔ）について辞書データの特徴ベクトル（→）ＸDを求める。
処理（Ｃ−１−８）：次に、特徴ベクトル（→）ＸDを入力ベクトルとして自己組織化アルゴリズムによって学習させる。ここでは、既に公知であるKohonenの自己組織化アルゴリズムを用いる。
処理（Ｃ−１−９）：学習終了後に得られる自己組織化特徴マップのノード（ｉ，ｊ）における特徴ノードベクトル（→）ＶD(i,j)に基づいて、辞書データの特徴ベクトルを入力した時に発火する特徴ノードベクトルとの対応関係を蓄積する。
【００４１】
ここで、自己組織化特徴マップの例が、図７〜図１３に示されている。図７は、階段降り動作を行った場合の自己組織化特徴マップの例、図８は、階段昇り動作を行った場合の自己組織化特徴マップの例、図９は、早歩き動作を行った場合の自己組織化特徴マップの例、図１０は、電車内歩行動作を行った場合の自己組織化特徴マップの例、図１１は、電車着座動作を行った場合の自己組織化特徴マップの例、図１２は、走行動作を行った場合の自己組織化特徴マップの例、図１３は、平地歩行動作を行った場合の自己組織化特徴マップの例、である。これらは、その動作が行われている最中にいずれかのノードを発火させ、そのノード内に記載されている数値を上昇させたマップである。そして、１つのノードの数値を上昇させる場合、そのノードに隣接する近傍の６つのノードについても数値を上昇させる。このマッピング処理を所定周期（例えば、数秒周期）で繰返し行う。このようにマッピングすることにより、数値の大きなノードほど発火の頻度が高いことになる。すなわち、このマッピング処理によって、発火度数を把握することができる。
【００４２】
非特許文献１においては、自己組織化特徴マップのノードについて区分けを行い（すなわちクラスタリングを行い）、その区分けを固定しそれに従って処理を行っている。このため、図８〜図１３のように発火ノードが比較的集中している場合には、非特許文献１に記載の技術を利用して処理を行っても問題無いが、図７のように発火ノードが分散している傾向がある場合には、非特許文献１に記載の技術を利用して処理を行うと正しいコンテキスト情報が得られない。
【００４３】
これに対し、本発明では、ノードの区分けを行わない。そして、ゆらぎのある質の悪いデータは学習対象とせず、質の良いデータのみを自己組織化特徴マップ中から抽出して学習対象とする。このため、図８〜図１３のように発火ノードが比較的集中している場合はもちろん、図７のように発火ノードが分散している傾向がある場合でも、正しいコンテキスト情報が得られる。
処理（Ｃ−１−１０）：上記対応関係の蓄積結果に基づいて、特徴ノードベクトル（→）ＶD(i,j)が発火した時に辞書データの特徴ベクトル（→）ＸDが特定信号パターンＣｎである確率Ｐ(i,j)（Ｃｎ）を算出する。
【００４４】
【数９】
0005297547.tif 000010
ここで、
【００４５】
【数１０】
0005297547.tif 000011
は特定信号パターンＣｎのうち特徴ノードベクトル（→）ＶD(i,j)で発火する確率を、Ｐ（Ｃn）は特定信号パターンＣnが発生する確率を表し、下記のように定義される。
【００４６】
【数１１】
0005297547.tif 000012
【００４７】
【数１２】
0005297547.tif 000013
ただし、
【００４８】
【数１３】
0005297547.tif 000014
はＺが生じる個数を、（→）ＶD(i,j)は特徴ノード（ｉ，ｊ）で特定信号パターンＣnが生じる事象を表す。
【００４９】
そして、全ての特徴ノード（ｉ，ｊ）についての｛（→）ＶD(i,j)｝を特定信号パターン特徴マップとして、｛Ｐ(i,j)（Ｃn）｝を特定信号パターン発火確率マップとして獲得する。
処理（Ｃ−１−１１）：この時、Ｐ(i,j)（Ｃn）が予め設定された閾値ＴPBを越えている特徴ノードを基底特徴ノードベクトル（→）ＶD(i,j)として抽出する。さらに基底特徴ノードベクトル（→）ＶD(i,j)を発火させた特徴ベクトル（→）ＸDを求め、それを算出する基となった加速度データを基底信号パターンセンサ信号（→）ａ(B)D（ｔ）として抽出する。
【００５０】
（より質の良い辞書データの構築）
図１４は、図５中の段階（Ｃ−２）の詳細な処理内容を示す図である。
図５中の段階（Ｃ−２）中の擬似データ生成処理（Ｓ２１）は図１４中の処理（Ｃ−２−１）に対応する。図５中の段階（Ｃ−２）中の特徴抽出処理（Ｓ２２）は図１４中の処理（Ｃ−２−２）から（Ｃ−２−５）までに対応する。図５中の段階（Ｃ−２）中の自己組織化学習処理（Ｓ２３）は図１４中の処理（Ｃ−２−６）に対応する。図５中の段階（Ｃ−２）中の学習結果解析処理（Ｓ２４）は図１４中の処理（Ｃ−２−７）及び（Ｃ−２−８）に対応する。以下、これらの処理の内容について説明する。
処理（Ｃ−２−１）：次に基底信号パターン（→）ａ(B)D（ｔ）から生成された擬似データである加速度ベクトル（→）ａ(D)D（ｔ）から特定信号パターン特徴マップを作成する方法について説明する。
各軸のまわりに予め設定された角度α、β、γをそれぞれ独立にΔθ刻みで変えながら以下の式に従って基底信号パターン（→）ａ(B)D（ｔ）を回転させることで擬似データ（→）ａ(D)D（ｔ）を生成する。
【００５１】
【数１４】
0005297547.tif 000015
【００５２】
そして、回転角の組み合わせ（α、β、γ）ごとに得られた擬似データ（→）ａ(D)D（ｔ）を、前述の辞書データから特定信号パターン特徴マップと特定信号パターン発火確率マップを獲得したのと同様の手順で擬似データに関する特定信号パターン特徴マップ｛（→）Ｖ(D)D(i,j)｝と特定信号パターン発火確率マップＰ(i,j)(D)（Ｃn）を得る。
【００５４】
処理（Ｃ−２−２）：まず、加速度ベクトル（→）ａD（ｔ）の特徴を解析するために基底関数展開することで時空間解像度レベルごとに時空間分解センサデータに分解する。ここでは局在波による関数展開であるMother WaveletをHarr関数とするWavelet Packet分解を利用する。上記の式（３）の両辺にWavelet Packet分解を施すと解像度レベルｐでの加速度ベクトル（→）ａD（ｔ）は上述した式（４）のように表せる。
処理（Ｃ−２−３）：次に、これら分解波形である基底関数
（→）ｕ(p、q)D（ｔ）＝
（（→）ｕ(p、q)x（ｔ）、（→）ｕ(p、q)y（ｔ）、（→）ｕ(p、q)z（ｔ））
のうち、加速度ベクトル（→）ａD（ｔ）を表現するのに最良基底関数となる基底関数を選択する。そこで、上記式（５）により情報エントロピーコスト関数を定義し、式（５）左辺を最小にするｐとｑとの組み合わせを選択する。
【００５５】
さらに、選択された基底関数（→）ｕ(p、q)D_BN（ｔ）の情報エントロピー（上述した式（６））がもっとも大きいものから選び、それらを順次加えていった和が全情報エントロピー（情報エントロピーの総和）に対して閾値ＴSoEに達するまでの基底関数を選択する。
処理（Ｃ−２−４）：そして、最終的に得られた基底関数（→）ｕ(p、q)D_F（ｔ）に基づいて辞書データの特徴ベクトル加速度ベクトル（→）ＸDを上記式（７）のように定義する。また、上記式（８）と同様であり、直流成分Ｘ(i)D_BN（０）の項は用いない。
処理（Ｃ−２−５）：また、このとき特徴ベクトルの分布統計量として、下記のモーメントを算出する、ここでは、分布統計量ベクトルを上記式（９）、式（１０）と同様にベクトル（→）Ｍ（（→）Ｘ0）として表す。
【００５６】
そして、辞書データとなる全ての加速度ベクトル（→）ａD（ｔ）について辞書データの特徴ベクトル（→）ＸDを求める。
処理（Ｃ−２−６）：次に、特徴ベクトル（→）ＸDを入力ベクトルとして自己組織化アルゴリズムによって学習させる。ここでは、既に公知であるKohonenの自己組織化アルゴリズムを用いる。
【００５７】
処理（Ｃ−２−７）：学習終了後に得られる自己組織化特徴マップのノード（ｉ，ｊ）における特徴ノードベクトル（→）ＶD(i,j)に基づいて、辞書データの特徴ベクトルを入力した時に発火する特徴ノードベクトルとの対応関係を蓄積する。自己組織化特徴マップの例は、すでに述べたので、その説明を省略する。
処理（Ｃ−２−８）：上記対応関係の蓄積結果に基づいて、特徴ノードベクトル（→）ＶD(i,j)が発火した時に辞書データの特徴ベクトル（→）ＸDが特定信号パターンＣｎである確率Ｐ(i,j)（Ｃn）を上記式（１２）によって算出する。
そして、全ての特徴ノード（ｉ，ｊ）についての｛（→）ＶD(i,j)｝を特定信号パターン特徴マップとして、｛Ｐ(i,j)（Ｃn）｝を特定信号パターン発火確率マップとして獲得する。
【００５８】
（実際の測定）
次に、実測データとなる加速度ベクトル（→）ａ0（ｔ）から特定信号パターンを抽出する方法について説明する。
前述で辞書データから特定信号パターン特徴マップと特定信号パターン発火確率マップを獲得したのと同様の手順で、実測データとなる加速度ベクトル（→）ａ0（ｔ）から特徴ベクトル（→）Ｘ0ならびに特徴ベクトルの分布統計量ベクトル（→）Ｍ（（→）Ｘ0）を作成する。
図１５は、図５中の段階（Ｃ−３）の詳細な処理内容を示す図である。
図５中の段階（Ｃ−３）中の前処理（Ｓ３１）は図１５中の処理（Ｃ−３−２）及び（Ｃ−３−３）に対応する。図５中の段階（Ｃ−２）中の特徴抽出処理（Ｓ３２）は図１５中の処理（Ｃ−３−４）から（Ｃ−３−７）までに対応する。図５中の段階（Ｃ−３）中の特定信号パターン検出処理（Ｓ３３）は図１５中の処理（Ｃ−３−８）及び（Ｃ−３−９）に対応する。以下、これらの処理の内容について説明する。
【００５９】
処理（Ｃ−３−１）：ユーザＢは複数軸センサデータが獲得できる加速度センサを搭載してある携帯端末を携行している。逐次、ユーザＢの状況（コンテキスト）に応じて変動する複数軸センサデータ（加速度ベクトル）が獲得され、パケット通信でプレゼンス提供サーバへ送信される。
処理（Ｃ−３−２）：プレゼンス提供サーバでは、ユーザＢの携帯端末から送信された加速度ベクトルの時系列データを逐次受信する。時系列データとして受信される加速度ベクトル（→）ａ0（ｔ）をデータ分割する。手法は限定されないが、ここでは、複数の時間窓を上記式（１）にしたがって用意してデータ分割する。
【００６０】
処理（Ｃ−３−３）：その後、データ補正処理を実施する。データ補正処理方法は加速度センサの用途によって異なるが、ここでは、Gaussianフィルターによるノイズ除去処理と、時間窓｛ｔk｝での加速度の正規化を行う。
続いて、センサデータ解析から特定信号パターンとの照合について述べる。
任意の姿勢で設置されたセンサから観測される時刻ｔでの実測データ加速度ベクトルを上記式（２）とし、理想的な姿勢で設置されたセンサから観測される時刻ｔにおける辞書データとなる加速度ベクトル（重力加速度も含む）を上記式（３）とおく。
【００６１】
はじめに、辞書データとなる加速度ベクトル（→）ａD（ｔ）から特定信号パターン特徴マップを作成する方法について説明する。
処理（Ｃ−３−４）：まず、加速度ベクトル（→）ａD（ｔ）の特徴を解析するために基底関数展開することで時空間解像度レベルごとに時空間分解センサデータに分解する。ここでは局在波による関数展開であるMother WaveletをHarr関数とするWavelet Packet分解を利用する。上記の式（３）の両辺にWavelet Packet分解を施すと解像度レベルｐでの加速度ベクトル（→）ａD（ｔ）は上記（４）のように表せる。
【００６２】
ここで、（→）ｕ(p、q)D（ｔ）は、解像度レベルｐのｑ番目の分解波形（基底関数）を表す。
処理（Ｃ−３−５）：次に、これら分解波形である基底関数
（→）ｕ(p、q)D（ｔ）＝
（（→）ｕ(p、q)x（ｔ）、（→）ｕ(p、q)y（ｔ）、（→）ｕ(p、q)z（ｔ））
のうち、加速度ベクトル（→）ａD（ｔ）を表現するのに最良基底関数となる基底関数を選択する。そこで、上記式（５）により情報エントロピーコスト関数を定義し、式（５）左辺を最小にするｐとｑとの組み合わせを選択する。
【００６３】
さらに、選択された基底関数（→）ｕ(p、q)D_BN（ｔ）の情報エントロピー（上述した式（６））がもっとも大きいものから選び、それらを順次加えていった和が全情報エントロピー（情報エントロピーの総和）に対して閾値ＴSoEに達するまでの基底関数を選択する。
処理（Ｃ−３−６）：そして、最終的に得られた基底関数（→）ｕ(p、q)D_F（ｔ）に基づいて辞書データの特徴ベクトル（→）ＸDを上記式（７）のように定義する。また、上記式（８）と同様であり、直流成分Ｘ(i)D_BN（０）の項は用いない。
処理（Ｃ−３−７）：また、このとき特徴ベクトルの分布統計量として、下記のモーメントを算出する、ここでは、分布統計量ベクトルを上記式（９）、式（１０）と同様にベクトル（→）Ｍ（（→）Ｘ0）として表す。
【００６４】
そして、辞書データとなる全ての加速度ベクトル（→）ａD（ｔ）について辞書データの特徴ベクトル（→）ＸDを求める。
処理（Ｃ−３−８）：そして、この（→）Ｘ0を入力ベクトルとして、前述で生成された辞書データならびに回転角の組み合わせ（α、β、γ）ごとの擬似データに基づいた特定信号パターン特徴マップ｛（→）Ｖ(Dk)D(i,j)｝と特定信号パターン発火確率マップ｛Ｐ(i,j)(Dk)（Ｃn）｝に基づいて特定信号パターン発火確率Ｐ(i,j)（Ｃn、（→）Ｘ0）を求める。
処理（Ｃ−３−９）：そして、以下の条件を満たすＣnを特定信号パターンとして検出する。すなわち、
【００６５】
【数１５】
0005297547.tif 000016
が特定信号パターンＣnとして検出される。
上記式（１６）で検出された特定信号パターンＣnと関連付けられているプレゼンス情報をプレゼンスデータベースから検索する。なお。プレゼンスデータベースは以下のデータ構造を有している。
【００６６】
【表１】
0005297547.tif 000017
【００６７】
ユーザＡの要求に対しプレゼンス提供サーバは、ユーザＢのセキュリティポリシーに従ってプレゼンス情報をユーザＡに対して送信することができる。
（実施例）
以上説明した複軸センサからの特定信号パターン検出方法を実現するための、特定信号パターン検出システムの実施例の構成について、図１６を参照して説明する。同図において、本例の特定信号パターン検出装置は、センサデータ抽出部（Sensory Data Extraction）１００ａと、センサデータ前処理部（Sensory Data Preprocessing）１００ｂと、ウェブレットパケット分解部（Wavelet Packet Decomposition）１００ｃと、最良基底関数選択部（Best Basis Function Selection Based on AIC）１００ｄと、特徴パターン算出部（XYZ−Feature Pattern by composition of Entropy Distribution on each axis）１００ｅと、特徴パターン算出部（XYZ−Feature Pattern by composition of Statistic Values of Entropy Distribution on each axis）１００ｆと、自己組織化特徴マップ（Self−Organizing Feature Map(Training Stage)）１００ｇと、特定信号パターン発生確率マップ（Specific Patterns Occurrence Probabilities Map based on Best Matching Units）１００ｈと、基底特徴パターン抽出部（Basis Feature Pattern Extraction）１００ｉと、基底センサデータ抽出部（Basis Sensory Data Extraction）１００ｊと、擬似センサデータ生成部（Pseudo−Sensory Data Generation）１００ｋと、センサデータ抽出部（Sensory Data Extraction）１００ｍと、自己組織化特徴マップ部（Self−Organizing Feature Map(Clustering Stage)）１００ｎと、特定信号パターン抽出部（Specific Patterns Extraction）１００ｐと、辞書特徴パターン抽出部（Dictionary Feature Map Extraction）１００ｑとを含んで構成されている。
【００６８】
これらの構成において、センサデータ抽出部１００ａは、上述した処理（Ｃ−１−１）を行う。
センサデータ前処理部１００ｂは、上述した処理（Ｃ−１−２）及び（Ｃ−１−３）、（Ｃ−３−２）及び（Ｃ−３−３）を行う。
ウェブレットパケット分解部１００ｃは、上述した処理（Ｃ−１−４）、（Ｃ−２−２）及び（Ｃ−３−４）を行う。
最良基底関数選択部１００ｄは、上述した処理（Ｃ−１−５）及び（Ｃ−２−３）、（Ｃ−３−５）を行う。
特徴パターン算出部１００ｅは、上述した処理（Ｃ−１−６）、（Ｃ−２−４）及び（Ｃ−３−６）を行う。
特徴パターン算出部１００ｆは、上述した処理（Ｃ−１−７）、（Ｃ−２−５）及び（Ｃ−３−７）を行う。
自己組織化特徴マップ１００ｇは、上述した処理（Ｃ−１−８）及び（Ｃ−２−６）を行う。
特定信号パターン発生確率マップ１００ｈは、上述した処理（Ｃ−１−９）及び（Ｃ−１−１０）、（Ｃ−２−７）及び（Ｃ−２−８）を行う。
基底特徴パターン抽出部１００ｉ及び基底センサデータ抽出部１００ｊは、上述した処理（Ｃ−１−１１）を行う。
【００６９】
擬似センサデータ生成部１００ｋは、上述した処理（Ｃ−２−１）を行う。
センサデータ抽出部１００ｍは、上述した処理（Ｃ−３−１）を行う。
自己組織化特徴マップ部１００ｎは、上述した処理（Ｃ−３−８）を行う。
特定信号パターン抽出部１００ｐは、上述した処理（Ｃ−３−９）を行う。
辞書特徴パターン抽出部１００ｑは、辞書特徴マップを出力する処理を行う。
【００７０】
このように、上述した３つの段階（Ｃ−１）、（Ｃ−２）及び（Ｃ−３）を行うため、重複する処理については、その処理を行う構成要素をそれぞれ設けるのではなく、重複する処理に共通して１つの構成要素を設けている。例えば、上記ウェブレットパケット分解部１００ｃは、上述した３回の処理（Ｃ−１−４）、（Ｃ−２−２）及び（Ｃ−３−４）に共通に設けられている。このように重複する処理に共通して１つの構成要素を設ければ、システムのハードウェア規模が増大することを防止できる。
【００７１】
なお、センサデータ抽出部１００ｍは、図２中のセンサデータ獲得機能部２１及びセンサデータ送信機能部２２、並びに、図３中のセンサデータ受信機能部１１に対応する。
また、センサデータ前処理部１００ｂ、ウェブレットパケット分解部１００ｃ、最良基底関数選択部１００ｄ、特徴パターン算出部１００ｅ、及び特徴パターン算出部１００ｆは、図３中のセンサデータ解析機能部１２に対応する。
【００７２】
さらに、センサデータ抽出部１００ａ、センサデータ前処理部１００ｂ、ウェブレットパケット分解部１００ｃ、最良基底関数選択部１００ｄ、特徴パターン算出部１００ｅ、特徴パターン算出部１００ｆ、自己組織化特徴マップ１００ｇ、特定信号パターン発生確率マップ１００ｈ、及び辞書特徴パターン抽出部１００ｑは、図３中のセンサデータ照合機能部１４に対応する。
なお、自己組織化特徴マップ部１００ｎは、図３中の特定信号パターン獲得機能部１３に対応する。
【００７３】
（まとめ）
以上説明したように、本発明によれば、以下の効果が得られる。
（１）センサを特定の場所に固定する必要がないので、携帯端末に搭載するセンサなどモバイル環境での携行方法に制約を与えないことが可能である。
（２）複軸センサデータを構成する各軸のセンサデータを一組のデータとして扱うことで、複軸センサデータの軌跡が非平面である特定パターンにも対応可能である。
（３）局所波を基底関数とする直交変換と情報エントロピーの情報量基準に基づいて特定信号パターンを分解することで、波形形状に応じた特徴抽出が可能である。
【００７４】
（４）自己組織化アルゴリズムで学習された結果を確率的手法に基づいてクラスタリングしているので、クラスターを代表する特徴ノードを学習データサンプルに対して忠実に決定することができる。
（５）任意の姿勢で設置されたセンサを前提に自己組織化アルゴリズムで学習させているので、極端に高速な捻りの動きでなければ、全ての並進運動の動きを検出することが可能である。
なお、上述した図３又は図１６中の構成要素の一部（すなわち機能の一部）を、携帯端末２００内に設けても良い。これにより、プレゼンス提供サーバ１００の処理負荷を軽減することができる。
【産業上の利用可能性】
【００７５】
本発明は、複数の軸方向ごとにセンサデータが検出できるセンサ搭載型携帯端末を用いて個人に適応したコンテキストを獲得する状況認識システムに利用されるものであり、プレゼンスサービス、コンテキスト通信、コンテキストに基づくサービスナビゲーションならびに健康アドバイスサービスなどのヘルスケアの分野に利用することができる。
【符号の説明】
【００７６】
１１  センサデータ受信機能部
１２  センサデータ解析機能部
１３  特定信号パターン獲得機能部
１４  センサデータ照合機能部
１５、１６  プレゼンス情報送信機能部
２１  センサデータ獲得機能部
２２  センサデータ送信機能部
２３  プレゼンス情報受信機能部
１００  プレゼンス提供サーバ
１００ａ  センサデータ抽出部
１００ｂ  センサデータ前処理部
１００ｃ  ウェブレットパケット分解部
１００ｄ  最良基底関数選択部
１００ｅ、１００ｆ  特徴パターン算出部
１００ｇ  自己組織化特徴マップ
１００ｈ  特定信号パターン発生確率マップ
１００ｉ  基底特徴パターン抽出部
１００ｊ  基底センサデータ抽出部
１００ｋ  擬似センサデータ生成部
１００ｍ  センサデータ抽出部
１００ｎ  自己組織化特徴マップ部
１００ｐ  特定信号パターン抽出部
１００ｑ  辞書特徴パターン抽出部
１０１  特定信号パターンデータベース
１０２  プレゼンスデータベース
２００、３００  携帯端末
【図１】
0005297547.tif 000018
【図２】
0005297547.tif 000019
【図３】
0005297547.tif 000020
【図４】
0005297547.tif 000021
【図５】
0005297547.tif 000022
【図６】
0005297547.tif 000023
【図７】
0005297547.tif 000024
【図８】
0005297547.tif 000025
【図９】
0005297547.tif 000026
【図１０】
0005297547.tif 000027
【図１１】
0005297547.tif 000028
【図１２】
0005297547.tif 000029
【図１３】
0005297547.tif 000030
【図１４】
0005297547.tif 000031
【図１５】
0005297547.tif 000032
【図１６】
0005297547.tif 000033