(57)【要約】【課題】  対象の基準フレームにおける慣性センサの向きにどのような前提条件も課さずに、慣性センサを用いて行動認識を行なう。【解決手段】  対象の行動を検出する行動認識装置が提供される。本装置は、直線運動を検知するよう構成された複数の直線運動センサおよび複数の回転運動センサを有し、該直線運動が互いに直交であり、該回転運動が互いに直交であるセンサユニット１０と、該センサユニットに含まれたセンサから信号を受信して処理し、該対象の行動を検知するよう構成された演算ユニット２０とを具備する。センサユニット１０は、前記対象に対して任意の向きで、該対象によって直接あるいは間接的に保持されている。演算ユニット２０は、直線運動センサおよび回転運動センサの両方からの信号を用いる演算を行ない、センサユニットの向きに依存することなく行動を判定する。【選択図】  図１
2006340903.tif 000002
【特許請求の範囲】
【請求項１】
対象の行動を検知するための行動認識装置において、
互いに直交する直線運動を検知するように構成された複数の直線運動センサおよび複数の回転運動センサを有し、該回転運動が互いに直交している、センサユニットと、
前記対象の行動を検知するために、前記センサユニットに含まれたセンサから信号を受信し処理するように構成された演算ユニットとを具備し、
前記センサユニットは、前記対象に対して任意の向きで、該対象により直接あるいは間接的に支持され、
前記演算ユニットは、前記センサユニットの向きに依存することなく前記対象の行動を判定するために、前記直線運動センサおよび前記回転運動センサの両方からの信号を使用した計算を行なうこと
を特徴とする行動認識装置。
【請求項２】
前記複数の直線運動センサが、３つの直交軸で加速度を測定する３つの加速度計であり、
前記複数の回転運動センサが、前記３つの直交軸のまわりの角速度を測定するための３つのジャイロセンサであり、
前記演算ユニットが、前記対象の現在の行動を判定するために、前記加速度計およびジャイロセンサからの６つの出力を使用すること
を特徴とする請求項１記載の行動認識装置。
【請求項３】
前記演算ユニットは、予め定めた期間複数組のセンサ出力から形成された観察シーケンスを格納するように構成されたバッファと、
異なる行動の動作特徴量に対応する複数の基準シーケンスを格納するように構成された記憶部と、
前記観察シーケンスと前記基準シーケンスとのマッチングを行い、最適マッチングした基準シーケンスを見つけるように構成されたマッチングプロセッサとを具備し、
前記マッチングプロセッサは、最適マッチングした基準シーケンスを得るために、最適の回転、および前記観察シーケンスと前記基準シーケンスの１つとの時間的な対応関係を見つけること
を特徴とする請求項１または２記載の行動認識装置。
【請求項４】
前記基準シーケンスが、隠れマルコフモデルでの状態を表わす複数のベクトルを含み、
前記観察シーケンスと前記基準シーケンスとの間のマッチングは、最適の状態経路を得るためのビタビアルゴリズムを用いて行なわれること
を特徴とする請求項３記載の行動認識装置。
【請求項５】
前記演算ユニットは、さらに、前記直線運動センサからの信号の自己相関性および前記回転運動センサからの信号の自己相関性を計算し、これらの自己相関性を合計して簡略化された自己相関関数を得ることで、前記対象の歩行レートを検知し、
前記簡略化された自己相関関数は、前記対象の歩行レートを決定するために使用されること
を特徴とする請求項１または２記載の行動認識装置。
【請求項６】
前記演算ユニットは、さらに、前記回転運動センサからの信号を積分することで第１の回転を算出し、重力モデルおよび前記直線運動センサからの信号を用いて第２の回転を算出し、該第１および第２の回転の加重補間を行なうことで、前記センサユニットの回転を検知し、
前記回転運動センサからの信号における変化がより小さい場合、前記第１の回転の重みがより大きく、
前記直線運動センサからの信号における変化がより小さい場合、前記第２の回転の重みがより大きいこと
を特徴とする請求項１または２記載の行動認識装置。
【請求項７】
基準行動の基準動作特徴量とそれに対応する位置とを格納するように構成された記憶部をさらに具備し、
前記演算ユニットは、さらに、前記センサからの信号で表わされる現在の行動が前記基準行動であるかどうかを判定し、該現在の行動が前記基準行動である場合には、最適マッチングした基準動作特徴量に対応する前記基準動作特徴量と前記現在の行動に対応する動作特徴量とをマッチングして、前記センサユニットの位置を決定すること
を特徴とする請求項１または２記載の行動認識装置。
【請求項８】
対象の行動を判定する行動認識方法において、
３つの直交軸での加速度および３つの直交軸のまわりの角速度をサンプリングし、
前記サンプリングされた加速度および角速度をすべて用いて、前記対象の行動を判定するものであり、
前記加速度および角速度は、前記対象に対して任意の向きで、該対象によって直接あるいは間接的に支持された単一のユニットに統合されたセンサによって検知され、
前記行動の判定は、前記センサユニットの向きに依存することなく前記行動を判定するために、前記加速度および角速度の両方を用いて行なわれること
を特徴とする行動認識方法。
【請求項９】
対象の行動を判定することにより、コンピュータに行動認識を行なわせるプログラムにおいて、
３つの直交軸での加速度および３つの直交軸のまわりの角速度のデータを受信し、
前記加速度および角速度の受信データをすべて用いて、前記対象の行動を判定するものであり、
前記加速度および角速度は、前記対象に対して任意の向きで、該対象によって直接あるいは間接的に支持されている単一のユニットに統合されたセンサによって検知され、
前記行動の判定は、前記センサユニットの向きに依存することなく前記行動を判定するために、前記加速度および角速度の両方を用いて行なわれること
を特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、行動認識装置、方法、プログラムおよび記憶媒体に関する。より詳細には、本発明は、行動認識装置、行動認識方法、該方法をコンピュータに実行させるプログラム、該プログラムを格納する記憶媒体またはメモリ、および該方法をコンピュータに実行させるためにコード化された信号に関する。
【背景技術】
【０００２】
人間の行動の分類は、カメラ、マイクロホン（非特許文献５）、慣性センサ（非特許文献４）およびゴニオメータにより試みられてきた。各センサあるいはセンサの組合せについては、研究者は、程度の差はあるものの、行動を分類することができるカスタムアルゴリズムも開発している。人間行動の分類の一部は人間の動作の正確な測定を得ることにあるので、使用される技術は、キャラクタアニメーション、運動科学およびバイオメカニクスに使用される方法とオーバラップしている。しかしながら、本願の発明者は、モバイルで拘束の無い状況下で、人間の行動をいかに判定し分類するかの問題に取り組んでいる。少数のセンサモダリティを除くすべてのものが適切である。
【０００３】
コンピュータ・ビジョン
そのように抑制された方法の一群においては、静止カメラまたは着用可能なカメラ（非特許文献１０）による対象者のモニタリングが含まれている。多くの場合、視覚的なトラッキングを単純化するために光学マーカが対象者の身体に付けられる。対象者の動きは、上記マーカの動きから逆運動学モデル（inverse-kinetic model）によって推論することができる。
【０００４】
この方法は、極めて豊富で完全な測定を提供することができるが、それは、カメラの視野内の、明るい空間で使用することができるだけで、移動性のモバイルアプリケーションでの使用は妨げられている。
【０００５】
コンピュータ聴覚
マイクロホン（身体上、および環境中）を使用して、研究者がユーザの行動の自動判定を試みた例（非特許文献２および非特許文献９を参照）もある。しかしながら、モーションセンサの主な利点は、より低いデータ速度（それははるかに低い計算負荷となる）およびユーザの動作だけを測定する能力である。オーディオにおいては、ユーザの行動によって引き起こされた音は、常に変化する環境音の背景と混ざっている。これら前景と背景の分離は、モーションセンサを使用する場合には生じない未解決の問題である。
【０００６】
絶対関節角度検知
この方法の群おいては、対象者の関節の角度を測定するための、対象者の身体へのセンサの直接装着が含まれている。可能なセンサとしては、電位差計および回転エンコーダ（これらは共に、対象者の関節角度の変化をセンサに送信するように外骨格（exoskeleton）に要求する）、曲げセンサ（曲げおよび屈曲を検知する材料片）並びに慣性センサ（ニュートンの第一運動法則によって主に速度変化を測定する加速度計およびジャイロ）である。病院関係では、重力に対する角度を測定するためにゴニオメータ（チルトセンサに相当）がしばしば使用される。
【０００７】
慣性検知
これら中でも、慣性センサおよびチルトセンサだけが、衣類（外骨格または曲げ検知材料片）への広範囲な改変を必要とせず、あるいは特別の付属品を必要とせずに、対象者の動きの測定に使用することができる。更に、慣性センサは、ユーザがポケットまたはバッグの中で単に保持するか運ぶことができる他の装置の内部に完全に収納することができる。このため、モバイル機器の中であるいは自立型機器としての使用においては、これらセンサは魅力的なものとなる。
【０００８】
他の関連研究、例えば、非特許文献７、非特許文献８、特許文献１、特許文献２も開示されている。
【０００９】
【非特許文献１】Dimitri P. Bertsekas. Nonlinear Programming. Athena Scientific, 2nd edition, 1995.
【非特許文献２】Brian Clarkson, Alex Pentland, and Kenji Mase. Recognizing user context via wearable sensors. In Fourth International Symposium on Wearable Computers （ISWC), 2000.
【非特許文献３】Masataka Ishikawa. United states patent 5,297,028: Method and apparatus for correcting drift errors in an angular rate sensor. Filed by Zexel Corporation Daihatsu-Nissan, Tokyo, Japan, August 1991.
【非特許文献４】Holger Junker, Jamie Ward, Paul Lukowicz, and Gerhard Troster. Continuous recognition of arm activities with body-worn inertial sensors. In Proceedings of the 8th International Symposium on Wearable Computers, 2004.
【非特許文献５】Paul Lukowicz, Jamie Ward, Holger Junker, Mathias Stager, Gerhard Troster, Amin Atrash, and Thad Starner. Recognizing workshop activity using body-worn microphones and accelerometers. In Second IEEE International Conference on Pervasive Computing and Communications, 2004.
【非特許文献６】Thomas P. Minka. Old and new matrix algebra useful for statistics. Technical report, Massachusetts Institute of Technology, 2000.
【非特許文献７】Bruno Scherzinger. United states patent 6,853,909: Walking stick navigator for position determination. Filed by Applanix Corporation, Inc., Ontario, Canada, November 2002.
【非特許文献８】Douglas Smith, Aaron Joseph, David Boone, Robert Borchers, and Ernest Burgess. United states patent 5,485,402: Gait activity monitor. Filed by Prosthetics Research Study, Seattle, WA, March 1994.
【非特許文献９】Mathias Staeger, Paul Lukowicz, and Gerhard Troester. Implementation and evaluation of a low-power sound-based user activity recognition system. In Eighth International Symposium on Wearable Computers （ISWC), 2004.
【非特許文献１０】Thad Starner, Josh Weaver, and Alex Pentland. Real-time american sign language recognition using desk and wearable computer-based video. IEEE Transactions on Pattern Analysis and Machine Intelligence, 20（12):1371-175, December 1998.
【非特許文献１１】Tero Tolonen and Matti Karjalainen. A computationally efficient multipitch analysis model. IEEE Transactions on Speech and Audio Processing, 8（6), November 2000.
【特許文献１】Robert Levi and Thomas Judd. United States Patent 5,583,776: Dead reckoning navigational system using accelerometer to measure foot impacts. Filed by Point Research Corporation, Santa Ana, CA, March 1995.
【特許文献２】Christopher Townsend and Steven Arms. United States Patent 6,834,436: Posture and body movement system. Filed by MicroStrain, Inc., Williston, VT, February 2002.
【発明の開示】
【発明が解決しようとする課題】
【００１０】
しかしながら、慣性センサで正確な動作測定をすることは、これらのセンサが持つ、衝撃、温度、電気的なノイズなどに対する感度のため、極めて困難であると歴史的に分かっている。例えば、ジャイロ（角速度センサ）は、温度に依存しているが、これは非特許文献３で説明されている。キャリブレーションも必要である。というのは、これらの装置が、典型的にはその出力に影響する任意のオフセットおよびスケールファクターを有するからである。
【００１１】
さらに、慣性センサは、対象者の基準フレームに関する概念を持っておらず、測定は、常に物理的なセンサに対して行われている。１つの例外としては、適切に向けられた加速度計が地球の重力場を測定することができるので、下方向を判定する方法を提供し得ることである。これが働くために、加速度計は他の未知の加速下にあってはならない。
【００１２】
従って、上述の課題を考慮して、対象者の基準フレームに関する慣性センサの向きについてどのような前提条件も課することなく、慣性センサを使用して行動認識を行ない、行動認識システムを実用に際しての柔軟性を高めることが望ましい。
【課題を解決するための手段】
【００１３】
本発明の一実施形態によれば、対象の行動を判定するための行動認識装置が提供される。前記行動認識装置は、互いに直交する直線運動を検知するように構成された複数の直線運動センサおよび複数の回転運動センサを有し、該回転運動が互いに直交しているセンサユニットと、前記対象の行動を検知するために前記センサユニットに含まれたセンサから信号を受信し処理するように構成された演算ユニットとを具備する。前記センサユニットは、前記対象に対して任意の向きで、該対象により直接あるいは間接的に支持され、前記演算ユニットは、前記センサユニットの向きに依存することなく前記対象の行動を判定するために、前記直線運動センサおよび前記回転運動センサの両方からの信号を使用した計算を行なう。
【００１４】
本発明の他の実施形態によれば、対象の行動を判定する行動認識方法が提供される。本方法は、３つの直交軸での加速度および３つの直交軸のまわりの角速度をサンプリングし、前記サンプリングされた加速度および角速度をすべて用いて、前記対象の行動を判定することを含むものである。本方法においては、前記加速度および角速度は、前記対象に対して任意の向きで、該対象によって直接あるいは間接的に支持された単一のユニットに統合されたセンサによって検知され、前記行動の判定は、前記センサユニットの向きに依存することなく前記行動を判定するために、前記加速度および角速度の両方を用いて行なわれる。
【００１５】
本発明の他の実施形態によれば、対象の行動を判定することにより、コンピュータに行動認識を行なわせるプログラムが提供される。本プログラムは、３つの直交軸での加速度および３つの直交軸のまわりの角速度のデータを受信し、前記加速度および角速度の受信データをすべて用いて、前記対象の行動を判定することを含むものである。前記加速度および角速度は、前記対象に対して任意の向きで、該対象によって直接あるいは間接的に支持されている単一のユニットに統合されたセンサによって検知され、前記行動の判定は、前記センサユニットの向きに依存することなく前記行動を判定するために、前記加速度および角速度の両方を用いて行なわれる。
【００１６】
本発明の実施形態においては、前記センサユニットは、任意の向きで、前記対象により直接あるいは間接的に支持される。言いかえれば、前記センサユニットからの信号は、前記対象の座標における前記センサの向きについての情報を何ら持っていない。しかしながら、前記センサユニット全体の異なる動きは、前記対象の異なる行動に対応する。
【００１７】
従って、前記センサユニット全体の動きを検知し、該検知された動きと、前記対象の対応する行動あるいは該行動に関係した対応する動作特徴量（signature）とを関連させることにより、前記センサユニットの向きに依存することなく前記対象の行動を判定することができる。
【発明の効果】
【００１８】
本発明によれば、対象の基準フレームに関する慣性センサの向きについてどのような前提条件も課することなく、慣性センサを使用して行動認識を行ない、行動認識システムを実用に際しての柔軟性を高めることができる。
【発明を実施するための最良の形態】
【００１９】
図面を参照して、本発明の実施形態について述べる。
【００２０】
図１は、本発明の一実施形態による、対象の行動を判定するための行動認識装置の構成例を示す。本装置は、センサユニット１０および演算ユニット２０を有する。センサユニット１０は、直線運動を検知するよう構成された複数の直線運動センサおよび複数の回転運動センサを有しており、該直線運動は互いに直交であり、該回転運動は互いに直交である。演算ユニット２０は、センサユニットに含まるセンサから信号を受信して処理し、対象の行動を出力するように構成されている。センサユニット１０は、対象に対して任意の向きで、対象により直接あるいは間接的に支持される。演算ユニット２０は、直線運動センサおよび回転運動センサの両方からの信号を用いて計算を行ない、センサユニットの向きに依存することなく行動を判定する。
【００２１】
前記複数の直線運動センサは、３つの直交軸での加速度を測定する３つの加速度計であってもよく、また複数の回転運動センサは、３つの直交軸のまわりの角速度を測定する３つのジャイロセンサであってもよい。演算ユニット２０は、加速度計およびジャイロセンサからの６つの出力を用いて、対象の現在の行動を判定するようにしてもよい。
【００２２】
前記演算ユニットは、予め定めた期間複数組のセンサ出力から形成された観察シーケンスを格納するように構成されたバッファと、異なる行動の動作特徴量に対応する複数の基準シーケンスを格納するように構成された記憶部と、前記観察シーケンスと前記基準シーケンスとのマッチングを行い、最適マッチングした基準シーケンスを見つけるように構成されたマッチングプロセッサとを具備してもよい。さらに、前記マッチングプロセッサは、最適マッチングした基準シーケンスを得るために、最適の回転、および前記観察シーケンスと前記基準シーケンスの１つとの時間的な対応関係を見つけてもよい。
【００２３】
さらに、前記基準シーケンスが、隠れマルコフモデルでの状態を表わす複数のベクトルを含み、前記観察シーケンスと前記基準シーケンスとの間のマッチングは、最適の状態経路を得るためのビタビアルゴリズムを用いて行なわれていてもよい。
【００２４】
また、前記演算ユニットは、さらに、前記直線運動センサからの信号の自己相関性および前記回転運動センサからの信号の自己相関性を計算し、これらの自己相関性を合計して簡略化された自己相関関数を得ることで、前記対象の歩行レートを検知し、前記簡略化された自己相関関数は、前記対象の歩行レートを決定するために使用されていてもよい。
【００２５】
また、前記演算ユニットは、さらに、前記回転運動センサからの信号を積分することで第１の回転を算出し、重力モデルおよび前記直線運動センサからの信号を用いて第２の回転を算出し、該第１および第２の回転の加重補間を行なうことで、前記センサユニットの回転を検知しても良い。さらに、前記回転運動センサからの信号における変化がより小さい場合、前記第１の回転の重みがより大きく、前記直線運動センサからの信号における変化がより小さい場合、前記第２の回転の重みがより大きい。
【００２６】
また、前記行動認識装置は、基準行動の基準動作特徴量とそれに対応する位置とを格納するように構成された記憶部をさらに具備してもよい。さらに、前記演算ユニットはさらに、前記センサからの信号で表わされる現在の行動が前記基準行動であるかどうかを判定し、該現在の行動が前記基準行動である場合には、最適マッチングした基準動作特徴量に対応する前記基準動作特徴量と前記現在の行動に対応する動作特徴量とをマッチングして、前記センサユニットの位置を決定してもよい。
【００２７】
本発明の他の実施形態によれば、対象の行動を判定する行動認識方法が提供される。本方法は、３つの直交軸での加速度および３つの直交軸のまわりの角速度をサンプリングし、前記サンプリングされた加速度および角速度をすべて用いて、前記対象の行動を判定することを含むものである。本方法においては、前記加速度および角速度は、前記対象に対して任意の向きで、該対象によって直接あるいは間接的に支持された単一のユニットに統合されたセンサによって検知され、前記行動の判定は、前記センサユニットの向きに依存することなく前記行動を判定するために、前記加速度および角速度の両方を用いて行なわれる。
【００２８】
本発明の他の実施形態によれば、対象の行動を判定することにより、コンピュータに行動認識を行なわせるプログラムが提供される。本プログラムは、３つの直交軸での加速度および３つの直交軸のまわりの角速度のデータを受信し、前記加速度および角速度の受信データをすべて用いて、前記対象の行動を判定することを含むものである。前記加速度および角速度は、前記対象に対して任意の向きで、該対象によって直接あるいは間接的に支持されている単一のユニットに統合されたセンサによって検知され、前記行動の判定は、前記センサユニットの向きに依存することなく前記行動を判定するために、前記加速度および角速度の両方を用いて行なわれる。
【００２９】
本発明のさらに他の実施形態を以下に説明する。以下の実施形態では、行動認識が可能な、センサユニットおよび演算ユニットを含むシステムが提供される。以下の実施形態では、センサユニットが対象によって「支持され」ており、センサユニットからの出力に基づいて、演算ユニットが対象の現在の行動のタイプを判定するものと想定している。
【００３０】
（１）システム概略
本実施形態においては、図２で示されるように、リアルタイムで様々な共通的人間行動を分類するために、６自由度（ＤＯＦ）慣性センサが、センサユニット１０として使用される。以下の説明では、用語「センサ」は６−ＤＯＦ慣性センサを指し、また、「ユーザ」はセンサを持っているか着用している対象者１を指し。更に、用語「装置」は、慣性計測装置（ＩＭＵ）と、センサおよびユーザ行動の分類方法のステップを実行するための前記演算ユニット２０とを組み合わせた本システムを指す。分類の結果は出力しても良く、あるいはアプリケーションに転送されても良い。
【００３１】
人間の行動を測定するために慣性センサを使用することが、以前に試みられたことはあるが、本システムを定義付ける特徴は、任意のセンサの向きで人体のほとんどいかなる場所でも完全に自動的に作動するという能力である。本実施形態において、本システムは、ユーザの行動および装置の状態の異なる態様のための４つの相異なる方法を統合している。その各々の方法は、センサの向きに関していかなる前提条件をも課すことなく、実行することができる。これら方法としては、歩行周波数測定方法２１、ユーザ上における装置の位置測定方法２２、動作トラッキング方法２３、および回転不変ビタビ（Rotation-invariant Viterbi）マッチング方法２４が含まれる（図２参照）。
【００３２】
センサは、ユーザにより直接または間接的に支持または保持されていても良い。例えば、装着の方法としては以下（完全に網羅しているものではない）が可能である。
・ユーザは、センサを手に持っている。
・ユーザは、センサを身に着けている財布かバックパック中に入れている。
・センサは、ユーザの体または衣類に固くまたはやや固く付けられている。
・センサは、ユーザの首かウエストに懸けられている。
・センサは、ユーザの靴、宝石、帽子、手袋、携帯用電子機器、その他ユーザが着用または持ち歩くことが通常の利用形態に含まれている物に組み込まれている。
・センサは、その時点でユーザを移動させている乗り物の内部に取り付けられているかまたはそこに置かれている。
【００３３】
また、本実施形態は、センサの動作がユーザの動作と関係を持つものであれば、センサとユーザの組み合わせがいかなる形式であっても、その適用が可能である。
【００３４】
（２）物理的装置
本実施形態の簡単で具体的な説明のために、６−ＤＯＦ慣性センサが、ｘ、ｙおよびｚ方向の線形加速度並びにｘ、ｙおよびｚ軸のまわりの角速度を測定する部品からできていると仮定する。ただし、以下に説明する方法は、６−ＤＯＦを測定するための他の方式に簡単に適用することができる。
【００３５】
本実施形態において、図３（ａ）で示されるように、本システムは、微小センサ１０ａおよびホスト装置３０ａに埋め込まれた演算ユニット２０を備えていてもよい。微小センサ１０ａは、６−ＤＯＦ慣性センサ１０ａ−１および演算ユニット２０ａにセンサデータを送信するための無線トランスミッタ１０ａ−２を含んでいる。あるいは、無線トランスミッタ１０ａ−２の代わりに、一定期間センサデータを格納するためのメモリが設けられていてもよい。微小センサ１０ａがホスト装置３０ａに接続された場合、格納されたデータが演算ユニット２０ａに転送されてもよい。
【００３６】
あるいは、６−ＤＯＦ慣性センサ１０ｂおよび演算ユニット２０ｂを含む本システム全体は、図３（ｂ）で示されるような携帯電話またはポータブルゲーム装置のような単一のホスト装置３０ｂに埋め込まれていてもよい。
【００３７】
（３）センサ信号調整およびキャリブレーション
以下、リアルタイムでユーザの行動の分類を行うための、センサからの生のセンサデータの演算処理を説明する。最初、センサの生の出力を表わすいくつかのシンボルが定義される：
【数１】
2006340903.tif 000003
【数２】
2006340903.tif 000004
【００３８】
本装置のセンサからの各サンプルは、数式（１）および（２）で特定される６つの値を含んでいる。ｕ（ｔ）センサ値は、ｉ−、ｊ−およびｋ−軸に沿った本装置の加速に直線的に比例する。これらは、例えばＭＥＭＳ加速度計か同様の機器によって測定されると考えることができる。ｖ（ｔ）センサ値は、ｉ−、ｊ−およびｋ−軸のまわりの本装置の角速度に直線的に比例する。これらは、例えばＭＥＭＳジャイロか同様の装置から得られると考えることができる。図４は、本装置自体のローカル座標システム（ｉ−ｊ−ｋ）での基準方向に関してすべての測定がなされることを示す。
【００３９】
ｐ（ｔ）∈Ｒ3が、任意の世界座標系（例えば図４の中のｘ−ｙ−ｚ）での、本装置の位置であり、ｇ（ｔ）∈Ｒ3が、ローカル座標システムでの地球の重力場の加速度であり、ｑ（ｔ）∈Ｒ3が、そのローカル座標システムでの本装置の向きであると仮定すると、以下のように書くことができる：
【数３】
2006340903.tif 000005
【数４】
2006340903.tif 000006
【００４０】
勾配（gradient）はすべて、本装置のローカル座標システムに関して得られる。ＡとＢは、３×３拡大縮小行列（scaling matrix）である。対角線成分は、センサ出力（例えば電圧）を標準物理単位（例えばｍ／ｓ2およびo／ｓ）に変換する倍率（感度としても知られている）を特定する。同時に、非対角成分は、多少の不整合（例えば、θおよびΦジャイロ間の角度は、正確には９０度ではない）によるセンサ間のクロストークについて記述している。典型的なセンサによって示された任意のオフセット（バイアスとしても知られている）は、ｕ0とｖ0によって示される。量子化雑音を含んでいるセンサノイズは、ランダム変数∈u（ｔ）および∈v（ｔ）のシーケンスとして、各センサについて特定される。
【００４１】
数式（３）および（４）は、本装置の実際の動作がセンサの測定にどのようにして変換されるのか説明している。（３）および（４）の式が加速度と角速度について解かれた場合、次の式が得られる：
【数５】
2006340903.tif 000007
【数６】
2006340903.tif 000008
【００４２】
信号キャリブレーションには、Ａ、ｕ0、Ｂ、ｖ0の値を決定することが含まれている。ＩＭＵ内部で使用される実際の部品に依存して、これらの値は、生産条件、周囲温度、駆動電圧などに依存している可能性がある。しかしながら、動作トラッキング方法だけは、これらのパラメータが既知であることを必要とする。ここで与えられた他の方法２１〜２２および２４のすべては、キャリブレーションを行なっていないセンサ値を使用することができる。
【００４３】
最後に、ＩＭＵからの連続的な信号があらかじめ定められたレートでサンプリングされ演算ユニットへ転送されるか、あるいは、該信号が予めＩＭＵでデジタル化される。本明細書の残りの部分では、ｒimuを、ＩＭＵの出力のデジタルサンプリングレートとする。したがって、ＩＭＵのデジタル化された出力は、次のように書くことができる：
【数７】
2006340903.tif 000009
【数８】
2006340903.tif 000010
【数９】
2006340903.tif 000011
【００４４】
（４）歩行周波数測定方法
このセクションでは、ユーザの行動を判定し、例えばユーザが歩いているか、走っているか、移動していないのかどうか等を判定し、その人の歩行レートをトラッキングする我々の方法について説明する。この方法は、ユーザの身体上の任意の位置に置かれた装置で有効である。歩調検出およびトラッキングのための本方法は、加速度センサの出力において明らかな反復動作の周期のトラッキングを含んでいる。
【００４５】
周期的信号の周期を検知・測定するための一般的方法では、信号の自己相関性（autocorrelation）におけるピークを見つける。加速度センサ出力（ｕ）の自己相関性が、ＩＭＵの実際の加速度の自己相関性にどのように関連しているかを確かめるために、ｉ番目のフレームの▽2ｐ（多変数）自己相関性は次のように書くことができる：
【数１０】
2006340903.tif 000012
【００４６】
ここで、＊が相互相関（cross-correlation）演算子を示し、そして、
【数１０−１】
2006340903.tif 000013
【００４７】
数式（１０）の第１項Φ1は、ＩＭＵおよび重力の加速度の周期性（ＩＭＵの向きのみに依存する）を含む最も重要な項である。第２項Φ2は、加速度と重力の間の相互相関である。人体に関するどんな種類の自然な反復運動についても、ｇ（ｔ）および▽2ｐは、同じ周期を共有することが、ほとんど保証されるので、その自己相関性Φ1および相互相関性Φ2は、建設的に加わる。最後の項Φ3は、様々な項のノイズおよびオフセットの相関性を含んでいる。これらの相関性は周期性を示さず、数式（１０）の最終の自己相関性に単にノイズおよび一定のバイアスを加えるであろう。ジャイロセンサ出力の自己相関性の有効性、ｖは正確に同じ議論に従う。
【００４８】
実用においては、自己相関性分析が行われるべき短い時間間隔が選ばれ、分析が行われるべきレートが選ばれる。Ｔgaitを１つの歩行分析フレーム（例えばセンサデータの間隔）の時間長とし、ｒgaitを歩行分析のレートとする。そうすれば、数式（１０）において、以下は１つの分析フレーム中のベクトルの数である。
【数１１】
2006340903.tif 000014
【００４９】
本方法には、ｕとｖの自己相関性の計算と、ＩＭＵ全体の動作の簡略化自己相関関数（ＳＡＣＦ）を得るために結果集計とが含まれる：
【数１２】
2006340903.tif 000015
【００５０】
時間的分解能の効率および改善のために、ウィーナ・ヒンチン定理（Wiener-Khinchin Theorem）が、ウインドウ付きの自己相関性を計算に使用される。一般に、ｘ［ｔ］が任意のシーケンスで、Ｘ［ω］＝Ｆ｛ｘ｝が離散的フーリエ変換（ＤＦＴ）であり、Ｗ［ω］＝Ｆ｛ｗ｝がウインドウ関数（例えばハミングウインドウのような）のＤＦＴであれば、次式のとおりである：
【数１３】
2006340903.tif 000016
【００５１】
高速フーリエ変換（ＦＦＴ）を用いて数式（１３）の中のＤＦＴを計算することで、数式（１２）の自己相関性の演算が極めて効率的となる。
【００５２】
ＡＣＦのピークは、データ中の任意の周期的な動作の周期に相当する。しかしながら、ＡＣＦが、任意周期における実周期の整数倍でのピークをも含むので、非特許文献１１で導入されたものに類似するピーク除去手順が、増強された簡略化自己相関関数（ＥＳＡＣＦ）を計算するために使用される。このピークを除去する技術は以下のとおりである：
＿＿＿＿＿＿＿＿＿___________________________________________________
ＰＲＵＮＥ−ＨＡＲＭＯＮＩＣＳ（ＳＡＣＦ）はＥＳＡＣＦを返す。
１．ＳＡＣＦを正の値にクリップする。
２．因数２で、アップサンプリングおよび補間をする。
３．クリップしたオリジナルから、アップサンプリングされたものを差し引く。
４．再度、正の値にクリップする。
５．３、５、７などの素数について、２〜４を繰り返す。
＿＿＿＿＿＿＿＿＿___________________________________________________
この手順の目的は、連続して各ピークの倍数をすべて削除し、その結果、基本周期だけを残すことにある。
【００５３】
最も高いピーク（非ゼロ遅延）およびしきい値が選択された場合、周期的な動作（例えば歩行のような）がある時間や、それがどの程度のピッチであるかをか決めることが可能である。τ*が、ＥＳＡＣＦで最も高いピークの遅延であれば、歩調ピッチは以下のとおりである：
【数１４】
2006340903.tif 000017
【００５４】
あるしきい値（分析期間内の信号におけるエネルギの割合として定義される）を超える値を有するピークがない場合、本システムは、ユーザが歩いてもいないし走ってもいない（歩調ピッチ＝０ヘルツ）と報告することができる。歩調ピッチは、経時的に変化するので、上記の分析は、連続する分析間隔のそれぞれについて繰り返し行われ、ｒgaitのレートで、歩調ピッチの評価がなされる。
【００５５】
（５）動作トラッキング方法
このセクションでは、下方向トラッキングを行ない、本装置の向きがいかに変わっても、または動いても、連続的に方向決定（heading）を行なう方法２３が提供される。この方法では、ＩＭＵからの出力はキャリブレーションされたものであると仮定する。
【００５６】
３軸ジャイロ検知データは、ＩＭＵの向きをトラッキングするための十分な情報を含んでいる。しかしながら、ジャイロにはすべて、ドリフトやその他のタイプのエラーの蓄積という課題がある。したがって、向きをトラッキングするために単なるジャイロ情報を用いると、短期では極めて正確であるが、長期では大きなトラッキングエラーとなる。
【００５７】
３軸加速検知データは、地球の重力が働いているのがどの方向であるかについての情報を含んでいる。地球の引力の方向についての知識によれば、最大２つの自由度で、ドリフトなしにＩＭＵの向きをトラッキングすることができる。この方法で加速度計を単に使用することについての問題は、ＩＭＵに作用が働いた（ある外部の力によって移動された）場合、地球の引力の推定が破綻するということである。
【００５８】
人間にとって一定の加速あるいは速度で移動したり、または他の対象者を移動させることは非常に難しく、また、不自然である。同様に、回転に関して、一定角速度で回転したり、または他の対象者を回転させることは人間にとって非常に難しい。
【００５９】
下方向をロバストに推定するためには、加速度計信号が一定の場合、加速度測定にもっと依存することが望ましい。というのは、人間が本装置に対して作用していない可能性が極めて高いからである。反対に、加速度計信号が一定の値ではない場合、ジャイロ測定により大きく依存することが望ましい。
【００６０】
この原理は、さらに、ジャイロ信号のオンライン再ゼロ化にも適用できる。言いかえれば、ジャイロ測定が一定の場合、そのジャイロ測定は再ゼロ化されることが望ましい。これは幸いにも、トラッキングの加速度測定にほとんど依存している場合でも真実である。
【００６１】
図６は、動作トラッキング方法のための演算パイプラインを示す。センサを含むＩＭＵ４０からの線形加速度データ４１０は、重力だけがＩＭＵ４０に作用しているという信頼度４１２を計算するために使用されると共に、垂直モデル４１４（後述する）への入力として使用される。ＩＭＵ４０からの角速度データ４２０は、ＩＭＵ４０が回転していないという信頼度４２２を計算するために使用されると共に、ジャイロ積分部４２４への入力として使用される。
【００６２】
垂直モデル４１４およびジャイロ積分４２４は、ＩＭＵの向きに関する個別の推定、つまりＱplumbおよびＱgyro（後述）を出力する。その後、ＱplumbとＱgyroは、信頼度スコアを用いて組合せられ、一方対他方（４３０）に重みを付けて、該重み付けされた補間を行ない、ＩＭＵの推定向きを計算する。ジャイロ積分は、前のタイムステップにおける向きに対する相対的な向きを計算するため、フィードバックが必要となる。
【００６３】
垂直モデル４１４は、向きの絶対的な推定を与えるので、フィードバックは必要としない。
【００６４】
（５.１）垂直モデル
２つの自由度でのみとはいえ、垂直モデルは、向きを測定するために地球の引力を用いる。要するに、このモジュールは、ＩＭＵ４０の座標系において、現在の下降ベクトル（down vector）を単位基準ベクトル（r）にマッピングするための回転を出力する（図４参照）。この基準ベクトルは、任意にあるいはゼロ化中に選ぶことができる。その結果として生じる回転は、次のような四元数（quaternion）形式で与えられる:
【数１５】
2006340903.tif 000018
【００６５】
ここで、ａ［ｔ］＝（ａx［ｔ］，ａy［ｔ］，ａz［ｔ］）Tは回転軸であり、θ［ｔ］は回転角である。下降ベクトルは、単に、正規化された線形加速度（重力以外の力がＩＭＵ４０に作用していないと仮定して）である：
【数１６】
2006340903.tif 000019
【００６６】
基準ベクトルおよび下降ベクトルの外積は回転軸を与える:
【数１７】
2006340903.tif 000020
【００６７】
ｓｉｎ（θ［ｔ］/２）（またｃｏｓ（θ［ｔ］/２）は、次のように、図７から直接計算される（埋込型アプリケーションで実行するのが難しい逆三角関数を使用しないで）：
【数１８】
2006340903.tif 000021
【数１９】
2006340903.tif 000022
【００６８】
（５.２）ジャイロ積分
ジャイロの角速度測定の積分は単純である。入力として、一つ前のＩＭＵの向き、Ｑ［ｔ−１］が必要である。最初、軸の回転はそれぞれ、ローカル座標（図４中のi−j−k）で四元数回転として表わされる：
【数２０】
2006340903.tif 000023
【数２１】
2006340903.tif 000024
【数２２】
2006340903.tif 000025
【００６９】
その後、これらのローカルの四元数は、一つ前のＩＭＵの向きＱ［ｔ−１］を用いて世界座標系へ回転される：
【数２３】
2006340903.tif 000026
【数２４】
2006340903.tif 000027
【数２５】
2006340903.tif 000028
【００７０】
そして、それらを組み合せることで、次の向きを表わす１つの四元数となる。
【数２６】
2006340903.tif 000029
【００７１】
（６）動作特徴量の比較方法
このセクションで示された方法２４（図２参照）は、未知の向きのＩＭＵにより集められた同期されていない動作特徴量をどのように比較するかという問題を解決する。例えば、典型的なアプリケーションでは、リアルタイムの信号を動作特徴量のデータベースにマッチングさせることを要求するであろう。しかしながら、ＩＭＵの向きは、データベース中の動作特徴量を生成するために使用されるＩＭＵの向きとは任意に異なっているかもしれない。更に、データベース中の動作特徴量に対してテスト信号を時間的に整合させる適切な方法を決定しなければならない。
【００７２】
本実施形態による回転不変ビタビ（Rotation-invariant Viterbi）アルゴリズム（ＲＩＶＩＴ）は、最適の回転および２つの動作特徴量の時間的な対応関係を見つける目的で設計されている。列が加速度および角度のコンポーネントの積み重ねられた値である行列として、考慮中の２つのシーケンスが表わされると仮定すると、
【数２７】
2006340903.tif 000030
【００７３】
シーケンスＸおよびＹの長さは、それぞれＴXとＴYである。
【００７４】
図８（ａ）および図８（ｂ）は、異なる長さおよびそれらの間の任意の回転を持っている可能性のある、２つの多次元シーケンスの例を示す。特に、この例は、座った状態での行動の動作特徴量のシーケンスであるが、ＩＭＵは各場合で異なる向きを持っていた。さらに、図８（ａ）中の座った状態での行動は、図８（ｂ）の場合より速かった。
【００７５】
ある任意の回転およびある任意の時間伸縮（time-warping）内で、これらの２つのシーケンス間の類似を計算する場合、問題は次の基準（measure）の最小化と述べてもよい：
【数２８】
2006340903.tif 000031
【００７６】
上記は、回転行列Ｒおよび対応行列（correspondence matrix）Ｃであるその自由パラメータについてであり、ここでシンボル「’」は転置行列を意味する。対応行列は、Ｘの時間次元（列）をＹの時間次元（列）上にマッピングする、特別の種類の射影行列（projection matrix）である。｛０，１，．．．，ＴX−１｝の各整数を、｛０，１，．．．，ＴY−１｝の整数に割り当てる関数（c（ｔ））を持っていると仮定すると、その対応行列は以下のように構成することができる：
【数２９】
2006340903.tif 000032
【００７７】
ＤR,C（Ｘ，Ｙ）を最小化する場合、Cを対応するｃ（ｔ）が減少しない行列に抑制することが望ましい。図８（ｃ）は、可能なマッピングの例を挙げる。最適Ｃのマップは、Ｙと整合させるためにＸをいかに伸ばし圧縮すべきか（いくつかのセクションを削除し、他のものを繰り返すことも可能）を決定する。最適の回転行列Ｒは、ＸのすべてのベクトルがＹのベクトルとマッチングするにはどのように回転すべきかを決定する。
【００７８】
数式（２８）を最小化にするためにここで用いられる方法は、座標降下法（Coordinate descent method：非特許文献１）と呼ばれる。並行演算に容易に適応される単純な反復法である。その主旨は、各反復中の１つの単一軸に沿って誤差関数を最小化にすることである。数式からわかるように、これは最小化問題を極めて単純化するが、各軸を順次処理していく場合、連続的に各座標を最小化することができる。本数式は極小値に収束することが知られている。そのときＲとＣに制約がなかった場合、ＤR,C（Ｘ，Ｙ）は凸状であるので、この極小はグローバルな最小でもある。しかしながら、本実施形態では制約があるので、ローカルの状態（localness）だけが保証される。次のセクションでは、各軸ＲおよびＣに関する誤差関数を最小化にする方法について述べる。
【００７９】
（６.１）Ｒに対するＤR,C（Ｘ，Ｙ）の最小化
Ｒだけに対してＤR,C（Ｘ，Ｙ）を最小化する簡単な方法がある。以下のように定義する：
【数３０】
2006340903.tif 000033
【００８０】
この最小化問題の解は、極分解（Polar decomposition：詳細な導き出し方については非特許文献６参照）として知られており、その解は次のとおりである：
【数３１】
2006340903.tif 000034
【００８１】
行列ＶおよびＵは、ＸＣＹ’＝ＵＤＶ’の特異値分解（Singular Value Decomposition：ＳＶＤ）から得られる。なお、ＸＣＹ’の大きさは６×６であるので、ＳＶＤのコストはシーケンスの長さに依存しない。
【００８２】
本願の発明者は、同じ回転がジャイロおよび加速度計の部品に適用されるように、Ｒを抑制することが有益であることを知見した。これは、Ｒの構造が、次のように対角行列であることが望ましいことを意味する：
【数３２】
2006340903.tif 000035
【００８３】
ここで、Ｑは３×３回転行列である。この制約は、Ｑ＝ＶＵ’の計算により容易に満たすことができる。ただし、ここでは、ＶとＵは、ＵXＣＵY'＋ＶXＣＶY'＝ＵＤＶ’のＳＶＤから得られる。
【００８４】
（６.２）Ｃに対するＤR,C（Ｘ，Ｙ）の最小化
動的時間伸縮（ＤＴＷ：dynamic time warping）を用いて、Ｃだけに対するＤR,C（Ｘ，Ｙ）を最小化する方法がある。以前と同様に、次のように定義する：
【数３３】
2006340903.tif 000036
【００８５】
Ｃを行列として直接見つける代わりに、基本的なマッピング関数c（ｔ）が計算される。Ｙシーケンスのベクトルが、左右トポロジ（left-right topology）を有する隠れマルコフモデル（ＨＭＭ）の中で状態を表わすとする。シーケンスＸは、最適経路シーケンスを提供するために、ビタビアルゴリズムを使用してＨＭＭとマッチングが行われる。Ｘのｊ番目のベクトルでの、ｉ番目のＨＭＭ状態（Ｙのｉ番目のベクトルを表わす）のスコアは、ユークリッド距離の２乗として計算される：
【数３４】
2006340903.tif 000037
【００８６】
これは、各ＨＭＭ状態、ｓt∈｛０，１，．．．，ＴＹ−１｝の出力分布を表わすために、Ｘのベクトルを中心とした単位ガウス分散（unit variance Gaussians）を使用することに相当する。遷移確率は次のとおりである：
【数３５】
2006340903.tif 000038
【００８７】
ここで、Ｋは実験的に決定した定数（我々の実験では、約（１００/１２８）Ｈｚ）であり、εは極めて小さい数である。あるいは、遷移確率を学習データの使用により決定してもよい。
【００８８】
Ｙを現わすとして上記のように決定されたパラメータを備えたＨＭＭおよびＸとして現された観察シーケンスを用いて、最適の（最大尤度）状態経路｛ｓt*｝を見つけるためにＨＭＭを介した標準ビタビアルゴリズムが使用される。この状態経路は、c（ｔ）＝ｓt*にしたがって、ＹのベクトルにＸの各ベクトルを割り当てる。対応行列は、上に述べられるようなｃ（ｔ）から直接計算することができる。なお、実際の運用においては、ＸからＸＣを構築するために、c（ｔ）を直接使用することで、行列積ＸＣがはるかに効率的に計算されることが注目されるべきである。
【００８９】
（６.３）初期条件の決定
上記座標降下法は、すべての降下法のように、座標の少なくとも１つについて初期値が必要となる。座標降下によって見つかった極小の質は、この初期値に極めて影響を受けやすい。ＤR,C（Ｘ，Ｙ）の最適化は、Ｃに関する初期推定から始められる。Ｃの初期値は、セクション（６.２）と同じ手順を用いて計算されるが、数式（３４）は次式で置き換えられる。
【数３６】
2006340903.tif 000039
【００９０】
したがって、回転行列Ｒに依存しない。このように、ＸとＹとの間での近似マッピングc（ｔ）を見つける場合、回転行列は必要とされない。
【００９１】
（６.４）ＸとＹとの間の歪みの計算
以下、誤差関数を最小化するためのＤＩＳＴＯＲＴＩＯＮと呼ばれる手順が提示される。これによって、２つの動作特徴量間の歪み指標（distortion measure）を計算する。ｉが反復を表し、Ｃ0がセクション（６.３）で計算されるＣについての初期推定を表わすとすると、εは収束を検知するために使用した停止基準である：
【数３７】
2006340903.tif 000040
【００９２】
Ｒ*（Ｃ）およびＣ*（Ｒ）はそれぞれセクション（６.１）および（６.２）で与えられた処理である。要約すれば、ＤＩＳＴＯＲＴＩＯＮ手順は、上に概説された座標降下法を用いて、ＲおよびＣに対する極小ＤR,C（Ｘ、Ｙ）を計算する。返された値Ｄfinalは、任意の回転または時間伸縮に対して不変である、ＸとＹの間の相違の指標である。
【００９３】
（６.５）装置構成
動作特徴量の比較について、このセクションで説明される方法は、演算ユニット２０によって実現され得る。該ユニットは、本方法の処理ステップを実行するためにＣＰＵとメモリを含んでいてもよい。より具体的には、例えば、演算ユニット２０は、バッファ２２、回転不変ビタビ（ＲＩＶＩＴ）処理部２４、隠れマルコフモデル（ＨＭＭ）記憶部２６および出力部２８を備えていてもよい。
【００９４】
バッファ２２は、対象者の現在の行動に対応する動作特徴量を示す観察シーケンスＸを構築するために、サンプリングされたセンサ出力を受信し格納する。本観察シーケンスＸは、最適の回転およびシーケンスＹとの時間的な対応（ＨＭＭ記憶部２６に格納されたＨＭＭによって表わされるシーケンスのうちの１つとの）を見つけるために、ＲＩＶＩＴ処理部２４に送られる。ＲＩＶＩＴ処理部２４では、最適の回転が最初に決定されてもよい。次に、最適の時間的な対応は、ビタビを用いてシーケンスＸをＨＭＭとマッチングすることにより決定され、最適経路シーケンスが提供されてもよい。
【００９５】
出力部２８は、マッチング結果を受信し、行動認識結果として、ユーザの現在の行動を示す信号を出力する。本実施形態では、行動は、前もってラベル付けられている又は分類されているか、ＨＭＭのシーケンスまたは状態に対してマッピングされていると想定している。マッピングテーブルは、出力部２８で格納されていてもよいし、あるいはＨＭＭのシーケンスまたは状態がラベル付けされていてもよい。
【００９６】
（７）身体上のセンサの位置を決定する方法
このセクションでは、本方法２２は、ユーザの身体に対してＩＭＵの位置を推定する。ユーザが行うどんな動きに対しても、ユーザの身体のどこにＩＭＵが着けられているかによって異なる測定結果が得られる。例えば、ポケットの中、左手での保持、ネックレス上などである。これら同じ行動に対する測定結果の差は、ＩＭＵの位置を推定するために使用することができる。もちろんこの目的のためには、あるタイプの動きは他のものよりよい。例えば、単に左手が関与する動作は、ＩＭＵがユーザの左のポケットにある状態と、ＩＭＵがユーザの右のポケットにある状態とを識別するために使用することはできない。しかしながら、歩き（walking）のような動作は、左右の対称性を除いて、ユーザの身体全体に亘っての唯一の測定結果を生ずる。
【００９７】
このセクションでは、ＩＭＵ位置を決定するために用いられる行動または動きは、基準動作と呼ばれる。任意の基準動作Ｘを仮定すると、全体的な処理パイプラインは、図５で示される。
【００９８】
基準動作の位置不変検出器（position-invariant detector）は、ＩＭＵ（５０）の出力をゲート制御するために使用され、このようにして、基準動作を含んでいるテスト特徴量のシーケンスは、次のモジュール２４に至る。例えば、基準動作が歩きの場合は、歩行周波数測定（セクション（４））のための方法２１を位置不変検出器として使用することができる。
【００９９】
モジュール２４では、テスト特徴量は、位置に関する基準動作特徴量の網羅的データベース５４とマッチングされる。最良のマッチングは、ＩＭＵの推定位置を示す。動作特徴量の総括的なマッチングについては、セクション（６）の方法が使用できる。さらに、特にそれが歩きのような周期運動である場合、特定の基準動作用のマッチングアルゴリズムをカスタムデザインすることも可能である。
【０１００】
本実施形態によれば、その最も基礎的な例示では、本装置は、ユーザがたった今行なったまたは行っている動作を算出して継続的に出力している。特に、出力するものは、本装置が利用可能な動作のマッチング用データベースによって決定される。この動作データベースは、アプリケーションの要求に応じて静的または動的なものとなりうる。
【０１０１】
本装置は、どんな人間のユーザにも付けられないが、むしろ機械や任意の動いている物体に付けられる構成で使用することも可能である。この例としては、乗り物（例えば、自転車、自動車、およびボート）、工場施設（例えば、ロボットアーム、およびＣＮＣ機械類）、動植物（例えば、犬および猫）、積荷（例えば、卵のクレート、または揮発性の液体のドラム）が含まれる。
【０１０２】
分類結果は、ユーザが行っていることまたは今行なっていることに基づいて、それ自体を構築する必要があるアプリケーションを動かすために直接使用することができる。これらのタイプのアプリケーションとしては、状況適応型（context-aware）と呼ばれ、それらの状況適応性は、この場合、本装置により可能になる。
【０１０３】
なお、添付の特許請求の範囲またはそれと均等の限りでは、様々な変更、組合せ、副次的組合せおよび改変が、設計必要条件および他の要因に依って生じ得ることは、当業者によって理解され得る。
【図面の簡単な説明】
【０１０４】
【図１】本発明の一実施形態による行動認識装置のブロック図である。
【図２】本発明の一実施形態による行動認識装置の概念図である。
【図３】図３（ａ）および３（ｂ）は異なる利用シナリオ用の行動認識装置の概要図である。
【図４】本発明の一実施形態による物理的モジュール（ＩＭＵ）の概念図である。
【図５】基準動作Ｘを用いる場合に、ユーザの身体上のＩＭＵの位置を決定する方法のための処理パイプラインの概略図である。
【図６】ＩＭＵの動作をトラッキングする方法のための処理パイプラインの概略図である。
【図７】垂直モデル中の四元数表現に必要な半角の正弦および余弦を直接計算するために使用される幾何学的構造を示す概要図である。
【図８】図８（ａ）および８（ｂ）は、動作特徴量のシーケンス例を示し、図８（ｃ）は、本発明の一実施形態による２つのシーケンスのマッチングを示す。
【図９】本発明の一実施形態による動作特徴量を比較する方法のための演算ユニットの概要図である。
【符号の説明】
【０１０５】
１...対象者
１０...センサユニット
２０...演算ユニット
２１...歩行周波数測定方法
２２...対象者上での装置の位置測定方法
２３...動作トラッキング方法
２４...回転不変ビタビマッチング方法
４０...慣性計測装置（ＩＭＵ）
【図１】
2006340903.tif 000041
【図２】
2006340903.tif 000042
【図３】
2006340903.tif 000043
【図４】
2006340903.tif 000044
【図５】
2006340903.tif 000045
【図６】
2006340903.tif 000046
【図７】
2006340903.tif 000047
【図８】
2006340903.tif 000048
【図９】
2006340903.tif 000049
【外国語明細書】
2006340903000001.pdf