JP6527024B2 - 移動運動解析装置及びシステム並びにプログラム - Google Patents

Description

本発明は、加速度センサ（sensor）を用いて人の歩行などの移動運動を評価するための技術に関する。

従来、医療施設や介護施設等において、歩行障害を持つ患者に対するリハビリテーション（rehabilitation）が行われている。

このリハビリテーションでは、一般的に、患者に適した歩行訓練プログラム（program）を作成したり、患者の回復レベル（level）を把握したりするために、理学療法士などの指導員が、患者の歩行運動を繰り返し評価する。

また、近年においては、加速度センサを用いて人の歩行などの移動運動中に生じる加速度を測定し、その測定結果に基づいて移動運動を定量的かつ客観的に評価する試みが成されている（例えば、特許文献１，要約等参照）。

特開２０１３−０５９４８９号公報

上記のように加速度の測定結果に基づいて移動運動を評価する場合、加速度データ（data）における各一歩に対応したデータが左右どちらの足によるものであるかを判別することは非常に重要である。

しかしながら、加速度データには複雑な波形が現れることが多く、一歩に対応したデータが左右どちらの足によるものであるかを判別することは難しい。

このような事情により、移動運動をする人に取り付けられた加速度センサからのデータにおいて、一歩に対応したデータが左右どちらの足によるものであるかを高い精度で判別することができる技術が望まれている。

第１の観点の発明は、
加速度センサを用いて得られた、人の移動運動中における加速度を表すデータに基づいて、前記人の左右方向における加速度の波形を特定する第１の処理と、
前記第１の処理で得られる加速度の波形を積分処理して速度の波形を生成する第２の処理と、
前記第２の処理で得られる速度の波形を積分処理して位置の波形を生成する第３の処理と、
前記第３の処理で得られる位置の波形から、該位置の波形を移動平均処理して得られる成分を除去する第４の処理とを行うことにより、処理済みの位置の波形を生成する第１の生成手段と、
前記処理済みの位置の波形に基づいて、左足及び右足の各々の一歩に対応する区間または時刻を特定する特定手段とを備えた移動運動解析装置を呈要する。

ここで、「第１の処理で得られる加速度の波形」は、第１の処理によって直接得られた加速度の波形に限らず、当該加速度の波形に何らかの処理が施された加速度の波形を含む。また、「第２の処理で得られる速度の波形」は、第２の処理によって直接得られた速度の波形に限らず、当該速度の波形に何らかの処理が施された速度の波形を含む。また、「第３の処理で得られる位置の波形」は、第３の処理によって直接得られた位置の波形に限らず、当該位置の波形に何らかの処理が施された位置の波形を含む。

第２の観点の発明は、
前記第１の生成手段が、前記第２の処理または前記第３の処理の後に、該処理で得られた波形における線形成分を除去する第５の処理を行う、上記第１の観点の移動運動解析装置を提供する。

第３の観点の発明は、
前記第１の生成手段が、前記第１の処理、前記第２の処理または前記第３の処理の前に、該処理で得られる波形における低周波成分を除去する第６の処理を行う上記第１の観点または第２の観点の移動運動解析装置を提供する。

第４の観点の発明は、
前記移動平均処理が、前記人の移動運動における実質的に二歩に相当する時間幅をウィンドウ幅とする処理である上記第１の観点から第３の観点のいずれか一つの観点の移動運動解析装置を提供する。

第５の観点の発明は、
前記二歩に相当する時間幅が、前記位置の波形に自己相関関数を適用して得られる時間である上記第４の観点の移動運動解析装置を提供する。

第６の観点の発明は、
前記二歩に相当する時間幅が、前記データに基づく前記人の上下方向の加速度と前後方向の加速度との線形結合を表す波形の周期に基づいて決定される上記第４の観点の移動運動解析装置を提供する。

第７の観点の発明は、
前記特定手段が、前記処理済みの位置の波形における上に凸となる波形部分と下に凸となる波形部分とを特定することにより、左足及び右足の各々の一歩に対応する区間または時刻を特定する上記第１の観点から第６の観点のいずれか一つの観点の移動運動解析装置を提供する。

第８の観点の発明は、
前記第６の処理が、０．４ヘルツ（Ｈｚ）以下である所定の周波数以下の成分を除去するように設定されたフィルタを用いて、前記低周波成分を除去する処理である上記第１の観点から第７の観点のいずれか一つの観点の移動運動解析装置を提供する。

第９の観点の発明は、
前記データに基づいて、前記人の加速度が表す波形を生成する第２の生成手段をさらに備えており、
前記特定手段が、該生成された波形における左足及び右足の各々の一歩に対応する区間または時刻を特定する上記第１の観点から第８の観点のいずれか一つの観点の移動運動解析装置を提供する。

第１０の観点の発明は、
前記データに基づいて、前記人の加速度に基づく波形を生成する第２の生成手段をさらに備えており、
前記特定手段が、該生成された波形における左足及び右足の各々の一歩に対応する区間または時刻を特定する上記第１の観点から第８の観点のいずれか一つの観点の移動運動解析装置を提供する。

第１１の観点の発明は、
前記特定手段が、前記第２の生成手段により生成された波形における極大値及び極小値のいずれか一方に対応した時刻に対して、左足が着地した時刻と右足が着地した時刻とを対応付ける上記第９の観点または第１０の観点の移動運動解析装置を提供する。

第１２の観点の発明は、
前記加速度センサと、上記第１の観点から第１１の観点のいずれか一つの観点の移動運動解析装置とを備えた移動運動解析システムを提供する。

第１３の観点の発明は、
コンピュータを上記第１の観点から第１１の観点のいずれか一つの観点の移動運動解析装置として機能させるためのプログラムを提供する。

上記観点の発明によれば、移動運動をする人に取り付けられた加速度センサからのデータにおいて、一歩に対応するデータが左右どちらの足によるものであるかを高い精度で認識することができる。

歩行解析システムの構成を概略的に示す図である。 加速度センサモジュール及び歩行解析装置のハードウェアの構成を示す図である。 加速度センサモジュール及び歩行解析装置の機能的な構成を示す機能ブロック図である。 加速度データ解析部の機能的な構成を示す機能ブロック図である。 歩行解析システムにおける処理の流れを示すフロー図である。 サンプル加速度データに基づいて生成された左右方向加速度波形、前後方向加速度波形、及び上下方向加速度波形を示す図である。 サンプル加速度データに基づいて生成された左右方向加速度波形、前後方向加速度波形、及び上下方向加速度波形の拡大図である。 サンプル加速度データに基づいて生成された加速度波形Ｗ_ｘ，Ｗ_y，Ｗ_zに、解析対象として決定された歩行期間Ｒ′を重ねて表した図である。 上下方向加速度波形Ｗ_zと歩行位相との関係を示す図である。 ステップ基準時刻検出部の機能的な構成を示す機能ブロック図である。 ステップ基準時刻検出処理の流れを示すフロー図である。 サンプル加速度データに基づく前後方向加速度波形Ｗ_y′及び上下方向加速度波形Ｗ_z′の拡大図である。 サンプル加速度データに基づく歩行同期波形ＷＢを示す図である。 歩行同期波形ＷＢにおける一歩区間Ｋを示す図である。 歩行同期波形ＷＢと上下加速度反映波形ＷＪ_z′とを示す図である。 上下加速度反映波形ＷＪ_z′における一歩区間Ｋごとに特定された区間内最大値Ｍ_kを示す図である。 ステップ基準時刻として検出された踵着地の時刻を示す図である。 踵着地左右判別部の機能的な構成を示す図である。 踵着地左右判別処理の流れを示すフロー図である。 サンプル加速度データに基づく左右方向加速度波形Ｗ_x′の拡大図である。 左右方向位置波形のトレンド除去前後での周波数分布の比較例を示す図である。 処理済みの左右方向の位置の波形ＷＬ_x′と各種の波形とを示す図である。 踵着地の左右判別の例を示す図である。 ステップ時間グラフ生成部の機能的な構成を示す機能ブロック図である。 ステップ時間グラフ生成処理の流れを示すフロー図である。 左右ステップ時間の特定例を示す図である。 通常歩行時のステップ時間グラフの例を示す図である。 右膝固定時のステップ時間グラフの例を示す図である。

以下、発明の実施形態について説明する。なお、これにより発明は限定されない。

図１は、歩行解析システム（system）１の構成を概略的に示す図である。なお、歩行解析システム１は、発明における移動運動解析システムの一例である。

歩行解析システム１は、図１に示すように、加速度センサモジュール（sensor module）２と、歩行解析装置３とを有している。加速度センサモジュール２は、患者１０の背面の腰部中央等に、粘着パッド（pad）やバンド（band）等により装着される。歩行解析装置３は、操作者１１が携帯したり操作したりして使用される。なお、歩行解析装置３は、発明における移動運動解析装置の一例である。

図２は、加速度センサモジュール２及び歩行解析装置３のハードウェア（hardware）の構成を示す図である。

図２に示すように、加速度センサモジュール２は、プロセッサ（processor）２１と、加速度センサ２２と、メモリ（memory）２３と、通信Ｉ／Ｆ（interface）２４と、バッテリ（battery）２５とを有している。歩行解析装置３は、例えば、スマートフォン（smart phone）、タブレット型コンピュータ（tablet computer）、ノートパソコン（note PC）などのコンピュータ端末であり、プロセッサ３１と、ディスプレイ（display）３２と、操作部３３と、メモリ３４と、通信Ｉ／Ｆ３５と、バッテリ３６とを有している。なお、プロセッサ２１及びプロセッサ３１は、それぞれ、単一のプロセッサに限定されず、複数のプロセッサである場合も考えられる。

図３は、加速度センサモジュール２及び歩行解析装置３の機能的な構成を示す機能ブロック（block）図である。

加速度センサモジュール２は、図３に示すように、加速度センサ部２０１と、サンプリング（sampling）部２０２と、送信部２０３とを有している。なお、サンプリング部２０２及び送信部２０３は、プロセッサ２１がメモリ２３に記憶されている所定のプログラム（program）を読み出して実行することにより実現される。

加速度センサ部２０１は、センサ本体を基準とした３次元直交座標系におけるｘ，ｙ，ｚの各軸方向の加速度成分について、その加速度成分に応じたアナログ（analog）信号をほぼリアルタイム（real time）に出力する。

サンプリング部２０２は、そのアナログ信号を所定のサンプリング周波数でサンプリングしてデジタル（digital）の加速度データに変換する。サンプリング周波数は、例えば１２８Ｈｚである。サンプリング部２０２は、例えば、１ｇ（重力加速度）＝９．８ｍ／ｓ²＝加速度データ値１２８となるスケール（scale）で、加速度データを出力する。

送信部２０３は、サンプリングされた各時刻における加速度成分を表す加速度データをほぼリアルタイムにて無線で送信する。

なお、本例では、加速度センサモジュール２は、センサ本体のｘ軸方向、ｙ軸方向及びｚ軸方向が、それぞれ、患者１０のＲＬ（Right-Left）方向、ＡＰ（Anterior-Posterior）方向及びＳＩ（Superior-Inferior）方向と一致するように取り付けられる。ＲＬ方向、ＡＰ方向及びＳＩ方向は、それぞれサジタル（sagittal）方向、コロナル（coronal）方向及びアキシャル（axial）方向とも言う。また、本例では、加速度センサモジュール２の姿勢（傾き）は、患者１０の歩行中において変化しないものと仮定する。

歩行解析装置３は、図３に示すように、操作部３０１と、ディスプレイ部３０２と、患者情報受付部３０３と、受信部３０４と、加速度データ取得制御部３０５と、加速度データ解析部３０７と、表示制御部３１０と、記憶部３１２とを有している。患者情報受付部３０３、加速度データ取得制御部３０５、加速度データ解析部３０７、及び表示制御部３１０は、プロセッサ３１がメモリ３４に記憶されている所定のプログラムを読み出して実行することにより実現される。

操作部３０１は、操作者１１の操作を受け付ける。操作部３０１は、例えば、タッチパネル（touch panel）、タッチパッド（touch pad）、キーボード（keyboard）、マウス（mouse）などにより構成されている。なお、操作者１１は、例えば、理学療法士などの指導員である。

ディスプレイ部３０２は、画像を表示する。ディスプレイ部３０２は、例えば、液晶パネル、有機ＥＬパネルなどにより構成されている。

患者情報受付部３０３は、患者情報の入力を受け付け、入力された患者情報を記憶部３１２に記憶させる。

受信部３０４は、加速度センサモジュール２の送信部２０３から送信された加速度データを無線で受信する。なお、送信部２０３と受信部３０４との無線通信には、例えば、ブルートゥース（Bluetooth（登録商標））等の規格を用いることができる。

加速度データ取得制御部３０５は、操作者１１による操作に基づいて加速度データを取得するよう受信部３０４及び記憶部３１２を制御する。

加速度データ解析部３０７は、取得された加速度データを解析して、その解析結果を出力する。加速度データ解析部３０７の詳細については後述する。

表示制御部３１０は、ディスプレイ部３０２の画面に、少なくとも加速度データの解析結果を含む種々の画像や文字情報などを表示するようディスプレイ部３０２を制御する。

記憶部３１２は、入力された患者情報、取得された加速度データ、加速度データの解析結果などを記憶する。なお、これらの情報は、必要に応じて、歩行解析装置３に接続されたデータベース（database）４１に転送されたり、外付けのＤＶＤ−ＲＯＭ、メモリカード（memory card）などの媒体や、インターネット（internet）を介して接続された外部の媒体などを含む記憶媒体４２に保存されたりする。

ここで、加速度データ解析部３０７の詳細について説明する。加速度データ解析部３０７は、取得された加速度データに対して解析処理を行い、その解析結果を出力する。解析処理は、複数用意されている。加速度データ解析部３０７は、操作者１１によって指定された解析処理を実行する。本例では、実行する解析処理として、取得された加速度データが担持する加速度成分の時間変化を表す加速度波形を生成し、その加速度波形から患者１０の個々の一歩の前進動作に掛かった時間であるステップ（step）時間を求め、求めたステップ時間時間をグラフ（graph）化する処理を想定する。

なお、一般的に、連続的な左右一歩ずつの前進動作の中には、右足の踵着地、左足のつま先蹴り、左足の踵着地、及び右足のつま先蹴りの各動作（歩行位相）が１つずつ含まれる。また、一歩の前進動作は、一方の足の着地から他方の足の着地までの間の動作として定義される。一歩の前進動作は、ステップ若しくは１ステップともいう。また、連続的な左右一歩ずつから成る二歩の前進動作は、一方の足の着地から他方の足の着地を経て再度の一方の足の着地までの間の動作として定義される。二歩の前進動作は、ストライド（stride）若しくは１ストライドともいう。また、一歩周期は、一歩の前進動作に要する時間として定義される。一歩周期は、一方の足の特定の歩行位相（例えば踵着地やつま先蹴り）から他方の足の同じ歩行位相までの所要時間になる。二歩周期または歩行周期は、一方の足の特定の歩行位相から次の一方の足の同じ歩行位相までの所要時間になる。

図４は、加速度データ解析部３０７の機能的な構成を示す機能ブロック図である。加速度データ解析部３０７は、上記の機能を実現させるため、図４に示すように、加速度成分算出部７１と、加速度波形生成部７２と、歩行期間特定部７３と、ステップ基準時刻検出部７４と、踵着地左右判別部７５と、ステップ時間グラフ生成部７６とを有している。

加速度成分算出部７１は、取得された加速度データに基づいて、データ取得期間の各サンプリング時刻における患者１０の左右方向、前後方向及び上下方向の加速度成分ａ_x，ａ_y，ａ_zをそれぞれ算出する。本例では、これらの加速度成分ａ_x，ａ_y，ａ_zは、重力加速度ｇの成分を除去して、患者１０の純粋な運動により生じた加速度成分として算出することを想定する。ただし、より簡便に、重力加速度ｇの成分を含む形で特定してもよい。また、左右方向、前後方向及び上下方向は、それぞれ、水平左右方向、水平進行方向及び鉛直方向を想定する。ただし、より簡便に、加速度センサモジュール２のセンサ本体を基準としたｘ軸方向、ｙ軸方向及びｚ軸方向としてもよい。なお、ここでは、加速度成分の正負は、左右方向では右側寄り、前後方向では前側寄り、上下方向では上側寄りをそれぞれ正とする。

加速度波形生成部７２は、算出された各方向の各サンプリング時刻における加速度成分ａ_x，ａ_y，ａ_zに基づいて、左右方向の加速度成分ａ_xの時間変化を表す左右加速度波形Ｗ_x、前後方向の加速度成分ａ_yの時間変化を表す前後加速度波形Ｗ_y、上下方向の加速度成分ａ_zの時間変化を表す上下加速度波形Ｗ_zをそれぞれ生成する。

歩行期間特定部７３は、加速度データ取得期間の中で患者１０が実際に歩行を行っている期間（以下、歩行期間ともいう）を特定する。

ステップ基準時刻検出部７４は、解析対象として決定された歩行期間中の加速度波形に基づいて、それぞれが患者１０の一歩の前進動作に対応する複数のステップ基準時刻を検出する。検出されたステップ基準時刻は、加速度波形において一歩分または複数歩分の前進動作に対応する部分波形を抽出する際に用いる。なお、ここでは、ステップ基準時刻として、踵着地のタイミングに対応した時刻を検出する。

踵着地左右判別部７５は、歩行同期波形を生成し、当該波形に基づいて、特定された踵着地が、それぞれ、左右どちらの足によるものであるかを判別する。

ステップ時間グラフ生成部７６は、検出されたステップ基準時刻に基づいて、歩行運動を構成する個々の一歩の前進動作について、その一歩の前進動作に掛かる時間であるステップ時間を求める。そして、このステップ時間の時系列的な変化を表すグラフを生成する。また、ステップ時間グラフ生成部７６は、求められたステップ時間に基づいて、患者１０の歩行評価に有用な種々の特徴量を算出する。

これより、歩行解析システム１における処理の流れについて説明する。

図５は、歩行解析システム１における処理の流れを示すフロー（flow）図である。

ステップ（step）Ｓ１では、患者情報受付部３０３が、患者情報の入力を受け付け、入力された患者情報を記憶部３１２に記憶させる。ここでは、操作者１１が、歩行解析装置３の操作部３０１を操作して、患者１０の患者情報を直接入力する。患者情報受付部３０３は、その直接入力された患者情報を記憶部３１２に記憶させる。患者情報には、例えば、患者のＩＤ番号、氏名、年齢、性別、生年月日などが含まれる。なお、後述する患者１０の加速度データやこの加速度データの解析結果などは、この患者情報と対応付けて記憶部３１２に記憶される。

ステップＳ２では、加速度データ取得制御部３０５が、受信部３０４及び記憶部３１２を制御して、患者１０の各時刻ｔ_iにおける加速度データを取得する。ここでは、まず、操作者１１が、患者１０の腰部に加速度センサモジュール２を取り付ける。そして、操作者１１は、歩行解析装置３の操作部３１により、加速度データの取得開始操作を行う。加速度データ取得制御部３０５は、この操作に応答して、受信部３０４に加速度データの受信を開始させ、記憶部３１２にその受信された加速度データの記憶を開始させる。次に、患者１０に、自身の標準的な歩行速度でしばらく歩行してもらう。歩行は、通常、距離にして５ｍ〜２０ｍ程度、時間にして２０秒〜３分程度、歩数にして１０歩〜４０歩程度である。加速度センサモジュール２のサンプリング部２０２は、加速度センサ部２０１の出力に基づいて、患者１０の歩行中におけるｘ軸方向、ｙ軸方向、ｚ軸方向それぞれの加速度成分Ａ_x，Ａ_y，Ａ_zをサンプリングして計測する。加速度センサモジュール２の送信部２０３は、計測された加速度成分を表す加速度データをほぼリアルタイムで送信する。この間、受信部３０４は、送信部２０３から送信された加速度データを順次受信し、記憶部３１２は、その受信された加速度データを記憶する。患者１０の歩行が終了したら、操作者１１は、操作部３１により加速度データの取得終了操作を行う。加速度データ取得制御部３０５は、この操作に応答して、受信部３０４に加速度データの受信を終了させる。これにより、加速度データの取得開始操作が成されてから取得終了操作が成されるまでの期間が実質的に加速度データ取得期間となり、この期間の各サンプリング時刻における各方向の加速度データが取得される。

ステップＳ３では、加速度データ解析部３０７が、加速度データに対して実行する解析処理を設定する。例えば、操作者１１が、操作部３０１を操作して、取得した加速度データをグラフ化して表示したり、取得した加速度データを解析してその結果を表示したりする複数の機能の中から、実行させたい所望の機能を選択する。加速度データ解析部３０７は、その選択された機能に応じて、実行させる解析処理を設定する。本例では、操作者１１は、患者１０のステップ時間の時系列的な変化を表すグラフと、ステップ時間を用いた特徴量とを表示する機能を選択するものとする。このような機能によれば、患者１０の歩行中における一歩一歩の前進動作に掛かる時間やその時間変化、左右の違い、左右のバランス（balance）などを容易に認識することができる。その結果、操作者１１は、患者１０の歩行運動における左右の足の動きのバランスや安定性、疲れ具合などを理解することができる。

ステップＳ４では、加速度成分算出部７１が、取得された加速度データを記憶部３１２から読み出し、当該加速度データに基づいて、加速度データ取得期間の各サンプリング時刻における患者１０の左右方向、前後方向及び上下方向の加速度成分ａ_x，ａ_y，ａ_zを算出あるいは特定する。なお、ここでは、加速度データが表す加速度から重力加速度ｇの成分を除去する処理を含む所定のアルゴリズム（algorithm）を用いて、各サンプリング時刻及び各方向の加速度成分を算出する。算出された加速度成分は、記憶部３１２に送信され記憶される。

ステップＳ５では、加速度波形生成部７２が、ステップＳ４で算出された患者１０の各サンプリング時刻における左右方向加速度成分ａ_x、前後方向加速度成分ａ_y、及び上下方向加速度成分ａ_zに基づいて、左右方向加速度波形Ｗ_x、前後方向加速度波形Ｗ_y、及び上下方向加速度波形Ｗ_zを生成する。本例では、加速度波形生成部７２は、加速度成分の各方向ごとに、加速度データの取得開始時点からの経過時間（時刻）と加速度成分とを２軸とした２次元座標系において、各時刻ｔ_iでの加速度成分ａ(i)に対応するデータ点[ａ(i), ｔ_i]をそれぞれプロット（plot）することにより加速度波形を生成する。加速度波形は、必要に応じて、平滑化処理やスムージング（smoothing）処理を行って滑らかな曲線にする。

図６は、サンプル加速度データに基づいて生成された左右方向加速度波形Ｗ_x、前後方向加速度波形Ｗ_y、及び上下方向加速度波形Ｗ_zを示す図である。横軸は、加速度データ取得開始から経過した時間ｔ（秒）であり、縦軸は、加速度データ値ａ_x，ａ_y，ａ_z（重力加速度ｇ／１２８）である。このサンプル加速度データは、約４０秒間に渡って取得されたものである。このサンプル加速度データを取得する際に、被検者は、時間ｔ＝７秒あたりで歩行を開始し、途中の時間ｔ＝２０〜２３秒あたりで歩行を一時停止してしゃがみ込み、時間ｔ＝３５秒あたりで歩行を終了している。生成された加速度波形には、被検者のそのような動作による加速度成分の変化が現れている。

図７は、サンプル加速度データに基づいて生成された左右方向加速度波形Ｗ_x、前後方向加速度波形Ｗ_y、及び上下方向加速度波形Ｗ_zの拡大図である。

人の歩行運動では、通常、一方の足の踵着地、他方の足のつま先蹴り、他方の足の踵着地、一方の足のつま先蹴りという４つの動作がこの順番で繰り返し行われる。

上下方向加速度波形Ｗ_zにおいては、図７に示すように、歩行運動を構成する上記４つの動作の各々に対応して、波高値が一定以上となる極大値すなわちピーク波形を取ることが知られている。また、一方（他方）の足の踵着地から他方（一方）の足のつま先蹴りまでの時間は、相対的に短くなり、一方（他方）の足のつま先蹴りから同じ一方（他方）の足の踵着地までの時間は、相対的に長くなる。

前後方向加速度波形Ｗzにおいては、図７に示すように、一方の足の踵着地から他方の足のつま先蹴りまでの一歩の前進動作と、他方の足の踵着地から一方の足のつま先蹴りまでの一歩の前進動作と対応して、波高が一定以上となる極大値すなわちピーク波形を取ることが知られている。

図５に戻り、ステップＳ６では、歩行期間特定部７３が、加速度データ取得期間の中で歩行期間を特定する。一般的に、加速度データ取得期間には、患者１０が歩行を行っている期間と歩行を行っていない期間とが含まれている。歩行を行っていない期間としては、例えば、加速度データの取得を開始してから患者１０が歩行を開始するまでの期間、患者１０が歩行を終了してから加速度データの取得を終了するまでの期間、患者１０が歩行中に一時的に歩行を止めてしまう期間などが挙げられる。一方、解析対象に歩行を行っていない期間の加速度データが含まれていると、正しい解析を行うことができない。そこで、ここでは、解析処理を行う前に、加速度データ取得期間の中で歩行期間を特定し、その歩行期間における加速度データを解析処理の対象として決定する。

一般的に、歩行期間を特定する方法としては、次のような方法が考えられる。

第１の歩行期間特定方法は、操作者１１が加速度波形を見て歩行期間と考える期間を手動で指定し、指定された期間を歩行期間として特定する方法である。

第２の歩行期間特定方法は、サンプリング時刻ごとに患者１０に生じた加速度の大きさを表す特徴量を求め、この特徴量が所定の閾値以上になった時点から当該閾値以下になった時点までを、歩行期間として特定する方法である。加速度の大きさを表す特徴量としては、例えば、重力加速度ｇの成分が除去された各方向の加速度成分ａ_x，ａ_y，ａ_zの平方二乗和が考えられる。

図８は、上記サンプル加速度データに基づいて生成された加速度波形Ｗ_ｘ，Ｗ_y，Ｗ_zに、解析対象として決定された歩行期間Ｒ′を重ねて表した図である。

図５に戻り、ステップＳ７では、ステップ基準時刻検出部７４が、解析対象となる歩行期間における加速度波形に基づいて、それぞれが患者１０の一歩一歩の前進動作に対応する複数のステップ基準時刻ｔｂ_jを検出する。以下に、一般的に考えられる検出方法を説明した後に、本例による検出方法について説明する。

一般的に、ステップ基準時刻を検出する方法としては、例えば、上下方向加速度波形Ｗ_zにおけるピーク（peak）波形や特定の波形パターンなど周期性を有する波形形状を検出することにより、ステップ基準時刻ｔｂ_jを検出する方法が考えられる。

図９は、上下方向加速度波形Ｗ_zと歩行位相との関係を示す図である。ステップ基準時刻ｔｂ_jを検出する方法としては、例えば、上下方向加速度波形Ｗ_zにおいて所定の閾値を超える極大値を特定し、その極大値に対応する時刻を踵着地の時刻として認識し、踵着地の時刻をステップ基準時刻ｔｂ_jとして検出する。

この方法は、アルゴリズムが簡単で分かりやすい。

しかしながら、加速度の閾値判定だけで判断しているので、特定精度に限界がある。例えば、歩行における一歩の前進動作とは異なる動作により患者１０の体が動いていると、何らかの加速度が生じ、誤検出することが考えられる。また例えば、閾値の設定によって検出タイミング（timing）が変化する。また例えば、加速度波形において一歩の前進動作に同期したピーク波形の波高値が単発的に弱くなってしまった場合に、検出漏れが発生したり、逆に一歩の前進動作とは異なる動作に伴って比較的強いピーク波形が現れた場合に、誤検出が発生する。また例えば、患者の年齢、体格、歩行障害の程度などによって動作時に生じる加速度の大きさが異なるため、患者１０によって検出精度にばらつきが生じる。

そこで、本例では、ステップ基準時刻検出部７４は、ステップ基準時刻を精度よく検出することができるように工夫された方法を用いて、ステップ基準時刻を検出する。以下、このようなステップ基準時刻検出部７４の機能的な構成と、そのステップ基準時刻検出処理について説明する。

図１０は、ステップ基準時刻検出部７４の機能的な構成を示す機能ブロック図である。ステップ基準時刻検出部７４は、図１０に示すように、加速度波形読取部７４１と、歩行同期波形生成部７４２と、一歩区間特定部７４３と、上下加速度反映波形生成部７４４と、区間内最大値特定部７４５と、踵着地認識部７４６とを有している。

図１１は、ステップ基準時刻検出処理の流れを示すフロー図である。

ステップＳ７１では、加速度波形読取部７４１が、前後方向加速度波形Ｗ_y及び上下方向加速度波形Ｗ_zを記憶部３１２から読み出す。

加速度波形読取部７４１は、さらに、読み出された前後方向加速度波形Ｗ_y及び上下方向加速度波形Ｗ_zのうち解析対象となる歩行期間の波形部分を、それぞれ前後方向加速度波形Ｗ_y′、上下方向加速度波Ｗ_z′として切り出す。

図１２は、上記のサンプル加速度データに基づく前後方向加速度波形Ｗ_y′及び上下方向加速度波形Ｗ_z′の拡大図である。横軸は、時間ｔ（相対値）を表しており、縦軸は、前後方向加速度成分ａ_y及び上下方向加速度成分ａ_z（相対値）を表している。

ステップＳ７２では、歩行同期波形生成部７４２が、歩行同期波形を生成する。歩行同期波形とは、歩行時の一歩一歩の動作に略同期して所定のパターンが繰り返される波形のことである。本出願人の実験結果によれば、このような波形は、前後方向加速度成分ａ_yと上下方向加速度成分ａ_zとの線形結合を表す波形として得られることが確認されている。そこで、ここでは、ステップＳ７１で得られた前後方向加速度波形Ｗ_y′における加速度成分ａ_yと、ステップＳ７１で得られた上下方向加速度波形Ｗ_z′における加速度成分ａ_zとに基づいて、次式により値Ｂを求め、その時間変化を表す波形ＷＢを歩行同期波形として生成する。
（数１）
Ｂ＝ｄｅｎｏｉｓｅ（α・ａ_y＋β・ａ_z）

ここで、ｄｅｎｏｉｓｅ（）は、ノイズを低減するための関数であり、主に高周波成分を低減する関数である。このような関数としては、ガウシアンフィルタ（Gaussian filter）等に相当する関数を考えることができる。また、α，βは、係数であり、いずれも１〜３程度の値である。

図１３は、サンプル加速度データに基づいて生成された歩行同期波形ＷＢを示す図である。

ステップＳ７３では、一歩区間特定部７４３が、一歩区間Ｋを特定する。一歩区間Ｋとは、歩行における一歩の動作に対応した時間的な区間のことである。一歩区間Ｋの特定方法は、種々考えられるが、基本的な考え方としては、歩行同期波形において一歩の動作に対応した所定の波形パターンあるいは波形条件を見つけ出し、その波形パターンあるいは波形条件に対応する時間的な範囲を一歩区間Ｋとして特定する。なお、歩行における特定の位相に対応した時刻を検出するためには、一歩区間Ｋは、その時刻が明らかに含まれるような区間とする必要がある。本例では、踵着地ＨＣの時刻を検出するので、一歩区間Ｋは、踵着地ＨＣのタイミングが含まれるような区間とする必要がある。

ここでは、図１４に示すように、歩行同期波形ＷＢにおいて、極小値を取る時刻と次の極小値を取る時刻との間の区間を一歩区間Ｋとしてそれぞれ特定する。なお、当然ではあるが、歩行同期波形が−Ｂの時間変化を表す波形である場合には、当該波形における極大値を取る時刻と次の極大値を取る時刻との間の区間を一歩区間Ｋとして特定するようにする。また、歩行同期波形ＷＢを低周波の曲線で近似し、その近似曲線より大きい値を持つ波形部分に対応する区間を、一歩区間Ｋとして特定してもよい。

ステップＳ７４では、上下加速度反映波形生成部７４４が、上下加速度反映波形ＷJ_z′を生成する。上下加速度反映波形ＷJ_z′とは、患者１０の上下方向加速度成分ａ_zが反映された波形のことであり、上下加速度反映成分J_z′の時間変化を表す波形であるともいえる。ここでは、上下加速度反映波形ＷJ′は、上下方向の上側への加速度の増大が正側への変化として現れる波形とする。上下加速度反映波形ＷJ′は、上下方向加速度成分ａ_zを歩行同期波形と同一の時間で表す波形であってもよいし、上下方向加速度成分ａ_zに基づいて演算される量を歩行同期波形と同一の時間で表す波形であってもよい。例えば、演算される量は、上下方向加速度成分ａ_zの項を含む演算式により演算される。また例えば、演算される量は、加速度の変化をより強調するため、上下方向加速度成分ａ_zを二階時間微分した項をｋ乗（ｋ≧１）した項を含む演算式により演算される。ここでは、上下加速度反映波形ＷJ_z′として、上下方向加速度成分ａ_zの時間変化を表す上下方向加速度波形Ｗ_z′に高周波成分を低減するフィルタ（filter）を適用して成る波形を用いる。

図１５は、サンプル加速度データに基づいて生成された、歩行同期波形ＷＢと上下加速度反映波形ＷJ_z′とを同一時間軸上で示した図である。

ステップＳ７５では、区間内最大値特定部７４が、一歩区間Ｋごとに区間内最大値Ｍ_kを特定する。区間内最大値Ｍ_kとは、上下加速度反映波形ＷJ_z′における一歩区間Ｋ内での最大値のことである。踵着地ＨＣは、基本的に、一歩区間Ｋ内に必ず一つ存在すると考えることができる。また、踵着地ＨＣは、上下加速度反映波形ＷJ_z′において、局所的な時間内のピーク波形として現れることが知られている。よって、区間内最大値Ｍ_kに対応する時刻は、踵着地ＨＣの時刻と対応付けることができる。なお、上下加速度反映波形ＷJ_z′が、上下方向における上側への加速度の増大が負側への変化として現れる波形である場合には、区間内最小値を特定し、その区間内最小値に対応する時刻を踵着地ＨＣの時刻と対応づける。

図１６は、上下加速度反映波形ＷJ_z′における一歩区間Ｋごとに特定された区間内最大値Ｍ_kを示す図である。図中、区間内最大値Ｍ_kは、丸印で示されている。

ステップＳ７６では、踵着地認識部７４６が、ステップＳ７５で特定した区間内最大値Ｍ_kに対応する時刻を、踵着地ＨＣの時刻と対応付けて認識する。ここでは、この踵着地ＨＣの時刻を、ステップ基準時刻とする。

図１７は、ステップ基準時刻として検出された踵着地ＨＣの時刻を示す図である。図中、ステップ基準時刻すなわち踵着地の時刻は、四角印で示されている。

なお、歩行同期波形としては、上記の波形ＷＢとは別の波形を用いることもできる。本出願人の実験結果によれば、例えば、次に示すような方法で生成される波形を用いることもできる。

まず、左右方向加速度成分、前後方向加速度成分、または上下方向加速度成分を時間積分処理して速度の波形を求める。次に、当該速度の波形に対して、線形成分の除去と、低周波成分の除去とを行う。これによって得られた波形を再び時間積分処理して位置（変位）の波形を得る。そして、当該位置の波形に対して、一歩周期または二歩周期に相当する時間幅を、平均化する幅すなわちウィンドウ（windows）幅として移動平均処理を行い、当該位置の波形の低周波成分を求め、この低周波成分を元の当該位置の波形から減算することにより、トレンド（trend）を除去する。

ここで、一歩周期は、例えば、一方の足の踵着地から他方の足の踵着地までの間に相当する時間である。また、二歩周期は、例えば、一方の足の踵着地から他方の足の踵着地を経て次の同じ一方の足の踵着地までの間に相当する時間である。

このようにして得られた波形は、歩行時の各一歩の動作に非常によく同期した波形であることが分かっており、歩行同期波形として好適である。

なお、上記の線形成分の除去及び低周波成分の除去は、場合によっては省略可能である。また、上記の線形成分の除去は、速度の波形の生成より後、移動平均を用いたトレンド除去より前であれば、いずれの段階で行ってもよい。また、上記の低周波成分の除去は、加速度の波形の特定より後、移動平均を用いたトレンド除去より前であれば、いずれの段階で行ってもよい。ただし、不要な成分は可能な限り早い段階で除去することが好ましいため、上記の順序で行うことが望ましい。また、一歩周期または二歩周期に相当する時間幅は、例えば、上記数式１によって得られる波形の周期から求めてもよいし、上記位置の波形に自己相関関数を適用して求めてもよい。

このように、歩行同期波形において波形の特徴を基に一歩一歩に対応した時間的な区間を特定し、加速度反映波形において区間内最大値に対応する時刻を特定する方法によれば、加速度反映波形における一歩の動作に対応したピーク波形の波高値が一歩ごとにばらついたり、一歩の動作とは異なる動作に伴うピーク波形が一歩の動作に対応したピーク波形に近接して現れたりしても、検出漏れや誤検出を抑えることができ、踵着地の時刻を精度よく認識することができる。その結果、例えば、患者１０の一歩ごとの加速度成分について、波形を観察したり、波形の解析を行ったりする上で、その精度を向上させることができる。

図５に戻り、ステップＳ８では、踵着地左右判別部７５が、ステップＳ７で認識されたそれぞれの踵着地について、左右どちらの足によるものであるかを判別する。以下に、一般的に考えられる判別方法を説明した後、本例による判別方向について説明する。

一般的に、踵着地の足の左右を判別する方法としては、例えば、左右方向加速度波形Ｗ_xにおいて、左足の一歩に特徴的な波形形状と右足の一歩に特徴的な波形形状とを検出することにより、左足の一歩に対応した時間と右足の一歩に対応した時間とを特定し、これらの時間と踵着地ＨＣの時刻とを対比して判別する方法が考えられる。

この方法は、実質的に既得である波形を用いることができ、考え方もシンプル（simple）なので、利用しやすい。

しかしながら、左右方向加速度波形Ｗ_xは、他の方向の加速度波形よりも複雑な波形を示すことが多く、このような波形から特徴的な波形形状を見つけて、左足の一歩に対応した時間と右足の一歩に対応した時間とを特定することは容易ではない。

一方、本出願人の実験結果によれば、左右方向加速度成分ａ_xの波形を二階時間積分して得られる左右方向の位置の波形Ｌ_xでは、ノイズ（noise）成分や解析に不要な成分が適切に除去されると、右足の一歩と左足の一歩とが、上に凸となる波形部分と下に凸となる波形部分とにかなり正確に対応して現れることが確認されている。

そこで、本例では、踵着地左右判別部７５は、踵着地の左右を精度よく判別することができるように工夫された、次に示すような方法を用いて、踵着地の左右を判別する。

以下、このような踵着地左右判別部７５の機能的な構成と、その踵着地左右判別処理について説明する。

図１８は、踵着地左右判別部７５の機能的な構成を示す図である。踵着地左右判別部７５は、図１８に示すように、加速度波形読取部７５１と、左右位置波形生成部７５２と、左右判別部７５３とを有している。

なお、加速度波形読取部７５１及び左右位置波形生成部７５２は、発明における生成手段の一例であり、左右判別部７５３は、発明における特定手段の一例である。

図１９は、踵着地左右判別処理のフロー図である。

ステップＳ８１では、加速度波形読取部７５１が、左右方向加速度波形Ｗ_xを記憶部３１２から読み出して特定する。

ステップＳ８２では、加速度波形読取部７５１は、さらに、特定された左右方向加速度波形Ｗ_xのうち解析対象となる歩行期間の波形部分を、左右方向加速度波形Ｗ_x′として切り出す。

図２０は、サンプル加速度データに基づく左右方向加速度波形Ｗ_x′の拡大図である。横軸は解析対象として決定された歩行期間における時間ｔを相対値で示しており、縦軸は、左右方向加速度成分ａ_xを相対値で示している。

ステップＳ８３では、左右位置波形生成部７５２が、左右方向加速度波形Ｗ_x′に対して、時間積分処理を行い、左右方向の速度の波形Ｖ_xを求める。

ステップＳ８４では、左右位置波形生成部７５２が、左右方向の速度の波形Ｖ_xに対して、線形成分を除去する処理を行う。この線形成分の除去には、例えば、上記データに基づく波形を１次関数で近似（線形近似）して、元の波形から近似した１次関数の成分を減算する手法を用いる。ここで線形成分を除去する主な目的は、積分前の波形が持つオフセット（off set）成分を除去することや、デジタルサンプリング（digital sampling）によって生じる積分誤差を除去することにある。

ステップＳ８５では、左右位置波形生成部７５２が、ステップＳ８４で得られた左右方向の速度の波形Ｖ_x′に対して、低周波成分を除去する処理を行う。低周波成分の除去には、例えば、ハイパスフィルタ（high pass filter）またはローカットフィルタ（low cut filter）を用いる。これらのフィルタは、例えば、０．４ヘルツ（Ｈｚ）以下である所定の周波数以下の成分を除去するように設定される。所定の周波数は、例えば、０．２ヘルツである。ここで低周波成分を除去する主な目的は、数値積分によって生じる誤差を取り除くことにある。数値積分では、周波数ωの逆数に比例した強さで低周波ノイズが増幅される。このため、周波数ωが０近傍になると、線形成分の除去や後述のトレンド（trend）の除去でも取り切れないほどの大きなノイズが発生する可能性がある。ハイパスフィルタによってこれを強制的に除去することで、計算精度が上がると考えられる。

ステップＳ８６では、左右位置波形生成部７５２が、線形成分及び低周波成分が除去された左右方向の速度の波形Ｖ_x″に対して、さらに時間積分処理を行って、左右方向の位置（変位）の波形Ｌ_xを得る。

ステップＳ８７では、左右位置波形生成部７５２が、ステップＳ８６にて得られた左右方向の位置の波形Ｌ_xに対して、時間軸方向に移動平均処理を行うことにより当該位置の波形Ｌ_xの低周波成分を求め、この低周波成分を元の当該位置の波形Ｌ_xから減算する処理を行う。移動平均処理におけるウィンドウ幅、すなわち平均化する時間幅は、例えば、実質的に歩行周期（２歩周期）に相当する時間幅とする。歩行周期に相当する時間幅は、例えば、一方の足の踵着地から他方の足の踵着地を経て、再び一方の踵着地が行われるまでの間に相当する時間である。この処理により、トレンドが除去された左右方向の位置の波形Ｌ_x′が得られる。ここで移動平均により求めた低周波成分を除去する主な目的は、左右方向の位置の波形Ｌ_xから、患者１０の実際の歩行周期と異なる周期を持つ波形成分を除去することにより、歩行周期と同じ周期を持つ波形だけを極力抽出することにある。実際の歩行では、歩行周期に同期した特定の動作のみが規則正しく繰返し行われるわけではなく、歩行周期に同期しない動作、例えば左右や前後に不規則に振ら付く動作などを含みながら前進する。そのため、位置の波形には、患者１０の歩行に伴って生じる歩行周期に同期した波形と同期しない波形とが含まれることになる。ここで取り扱いたい波形は、歩行周期に同期した波形のみである。そこで、歩行周期に同期しない波形は、トレンドとして極力排除する。

図２１は、トレンド除去前の左右方向の位置Ｌ_xの波形ＷＬ_xにおける周波数分布ＤＬ_xと、トレンド除去後の左右方向の位置Ｌ_x′の波形ＷＬ_x′における周波数分布ＤＬ_x′との比較例を示す図である。横軸は周波数を表しており、縦軸はフーリエ（Fourier）成分ノルム（norm）を表している。

加速度における１ヘルツ以下の振動は、主に患者１０の重心の並進運動に起因するノイズであり、歩行運動を推定する上で本質的な情報ではない。ここでの移動平均処理を用いたトレンド除去によって、１ヘルツ以下の長周期振動を強く減衰させることができる。つまり、これにより、目的の成分以外の解析に不要な成分が除去された処理済みの位置の波形を得ることができる。

図２２は、処理済みの左右方向の位置の波形ＷＬ_x′とステップＳ７で得られた各種の波形とを同一の時間軸上で表す図である。

このようにして得られた左右方向の位置の波形ＷＬ_x′は、上または下に凸となる波形部分が、歩行時の一歩の動作と時間的によく対応した波形であることが分かる。また、当該波形は、波形部分の形状が上に凸であるか下に凸であるかによって、その波形部分に対応した一歩の動作が左右どちらの足によるものであるか判別可能な波形であることも分かる。

なお、ステップＳ８４の線形成分を除去する処理と、ステップＳ８５の低周波成分を除去する処理とは、必須の処理ではない。線形成分や低周波成分は、ステップＳ８７のトレンドを除去する処理により、ある程度の除去が可能であると考えられるからである。しかしながら、トレンドを除去する処理だけで、不要な線形成分及び低周波成分をすべて除去できるわけではない。したがって、歩行周期に同期した波形の抽出精度を向上させるためには、当然ながら、線形成分及び低周波成分の除去は行った方がよい。

また、線形成分を除去する処理を行うタイミングと、低周波成分を除去する処理を行うタイミングとは、本例のタイミングに限定されない。

例えば、線形成分を除去する処理は、速度の波形を生成する処理と、移動平均を用いてトレンドを除去する処理との間であれば、どのタイミングで行ってもよい。つまり、線形成分を除去する処理は、速度の波形を生成する処理または位置の波形を生成する処理の後に、当該処理で得られる波形における線形成分を除去する処理とすることができる。

また例えば、低周波成分を除去する処理は、加速度の波形を特定する処理の後、移動平均を用いてトレンドを除去する処理の前であれば、どのタイミングで行ってもよい。つまり、低周波成分を除去する処理は、加速度の波形を特定する処理、速度の波形を生成する処理または位置の波形を生成する処理の後に、当該処理で得られる波形における低周波成分を除去する処理とすることができる。

ただし、ノイズ成分によるエラー（error）の拡散や積算を抑える、あるいは、取り扱う値を小さく抑えて演算効率を高めるといった観点からは、線形成分を除去する処理及び低周波成分を除去する処理は、できるだけ早い段階で行うことが望ましい。したがって、本例のように、これらの処理を、速度の波形を時間積分処理する前に行うことは、好適な例の一つである。

ステップＳ８８では、左右判別部７５３が、ステップＳ８７で得られた左右方向の位置の波形Ｌ_x′に基づいて、認識されたそれぞれの踵着地の左右の判別を行う。

図２３は、踵着地の左右判別の例を示す図である。

前述のとおり、ここでは、左右方向の加速度成分の正負は、右側寄りを正としている。したがって、左右方向の位置の波形ＷＬ_x′では、正側寄りすなわち上に凸となる波形部分は、右足ステップに対応し、負側寄りすなわち下に凸となる波形部分は、左足ステップに対応する。そこで、左右方向の位置の波形ＷＬ_x′において、上に凸となる波形部分に時間的に対応した踵着地ＨＣを、右足の踵着地と判別する。また、左右方向の位置の波形において、下に凸となる波形部分に時間的に対応した踵着地ＨＣを、左足の踵着地と判別する。

判別アルゴリズムは、種々考えられるが、例えば、次のようなアルゴリズムを考えることができる。
ステップＴ１：判別対象である踵着地ＨＣの時刻が含まれる一歩区間Ｋを、左右方向の位置の波形ＷＬ_x′における凹凸を判定する区間［ｔ₀，ｔ₁］に決定する。
ステップＴ２：座標（ｔ₀，Ｌ_x′(ｔ₀）)と座標（ｔ₁，Ｌ_x′(ｔ₁）)とを結ぶ直線を表す関数ｆ（ｔ）を決める。
ステップＴ３：区間［ｔ₀，ｔ₁］でｆ（ｔ）−Ｌ_x′(ｔ）を積分し、その正負を求める。
ステップＴ４：積分値が正であれば、左右方向の位置の波形ＷＬ_x′の一歩区間Ｋにおける波形部分は下に凸であり、踵着地ＨＣは左足の着地であると判定する。積分値が負であれば、左右方向の位置の波形ＷＬ_x′の一歩区間Ｋにおける波形部分は上に凸であり、踵着地ＨＣは右足の着地であると判定する。

また例えば、次のようなアルゴリズムを考えることもできる。
ステップＵ１：移動平均処理等により、左右方向の位置の波形ＷＬ_x′の低周波成分（直線に近い成分）の波形ＷＣを求める。移動平均処理におけるウィンドウ幅は、一歩周期の数倍程度に設定する。
ステップＵ２：判別対象である踵着地ＨＣの時刻ｔ_HCにおける左右方向の位置の座標Ｌ_x′（ｔ_HC）が、波形ＷＣの正側と負側のどちらにあるかを判定する。
ステップＵ３：座標Ｌ_x′（ｔ_HC）が波形ＷＣの負側にあれば、左右方向の位置の波形ＷＬ_x′の踵着地ＨＣに対する波形部分は下に凸であり、踵着地ＨＣは左足の着地であると判定する。座標Ｌ_x′（ｔ_HC）が波形ＷＣの正側にあれば、左右方向の位置の波形ＷＬ_x′の踵着地ＨＣに対応する波形部分は上に凸であり、踵着地ＨＣは右足の着地であると判定する。

これにより、ステップＳ７で特定された複数のステップ基準時刻を、左足のステップ基準時刻と右足のステップ基準時刻とにそれぞれ分けて特定することができる。

このように、左右方向の位置の波形において波形の特徴を基に左足の一歩に対応した時間的な区間と右足の一歩に対応した時間的な区間とを特定し、検出された踵着地の時刻が左右どちらの区間内にあるかによって踵着地の左右を判別する方法によれば、左足の一歩に対応した波形形状と右足の一歩に対応した波形形状とが安定して現れることが確認されている波形を用いて左足の一歩に対応した区間と右足の一歩に対応した区間とを精度よく求めることができ、踵着地の左右を精度よく判別することができる。その結果、患者の歩行動作における左右の足の動作の違いやバランスを定量的かつ客観的に評価する上で、その精度を向上させることができる。

図５に戻り、ステップＳ９では、ステップ時間グラフ生成部７６が、患者１０の一歩一歩の前進動作に掛かる時間であるステップ時間の時系列的な変化を表すグラフの生成と、そのステップ時間に係る幾つかの特徴量の算出とを行う。なお、ステップ時間は、種々の定義の仕方が考えられるが、本例では、ステップ基準時刻を用いて定義することとし、一歩の前進動作に対応したステップ基準時刻と次の一歩の前進動作に対応したステップ基準時刻との間の時間として求めることにする。

以下、ステップ時間グラフ生成部７６の機能的な構成と、ステップ時間グラフ生成処理について説明する。

図２４は、ステップ時間グラフ生成部７６の機能的な構成を示す機能ブロック図である。ステップ時間グラフ生成部７６は、図２４に示すように、ステップ時間特定部７６１と、グラフ生成部７６２と、特徴量演算部７６３とを有している。

図２５は、ステップ時間グラフ生成処理の流れを示すフロー図である。

ステップＳ９１では、図２６に示すように、ステップ時間特定部７６１が、検出されたステップ基準時刻ｔｂ_jとその足の左右の判別結果とに基づいて、左ステップ時間Ｔ_Liと右ステップ時間Ｔ_Riとを特定する。左ステップ時間Ｔ_Liは、患者１０の左足による一歩の前進動作に対応したステップ時間であり、左の踵着地ＨＣに対応したステップ基準時刻から右足の踵着地ＨＣに対応したステップ基準時刻までの時間である。また、右ステップ時間Ｔ_Riは、患者１０の右足による一歩の前進動作に対応したステップ時間であり、右の踵着地ＨＣに対応したステップ基準時刻から左足の踵着地ＨＣに対応したステップ基準時刻までの時間である。

特定された左ステップ時間Ｔ_Lと右ステップ時間Ｔ_Rとの差異は、左足による一歩の前進動作と右足による一歩の前進動作との違いが反映される。そのため、患者１０の足の前進動作の左右差を、これら左ステップ時間Ｔ_L及び右ステップ時間Ｔ_Rによって定量化して評価することが可能になる。

ステップＳ９２では、グラフ生成部７６２が、求められた複数のステップ時間Ｔ_jの時間変化を表すグラフを生成する。

本例では、横軸をステップ番号ｊ、縦軸をステップ時間（秒）とした棒グラフを生成する。ステップ番号ｊとは、ステップ時間Ｔの時系列的な順番を表す番号である。また、本例では、左ステップ時間Ｔ_Lと右ステップ時間Ｔ_Rとが互いに区別しやすいように、左ステップ時間Ｔ_Lと右ステップ時間Ｔ_Rとで棒の色または模様（柄、パターン）を変えて描くようにする。

図２７及び図２８は、生成されたグラフＧの一例を示す図である。図２７は、通常の歩行によるサンプルであり、図２８は、右膝を金属棒で固定した歩行によるサンプルである。なお、ここで生成するグラフＧは、本例のような棒グラフだけでなく、折れ線グラフやスプライン（spline）補間を用いたスプライン曲線グラフ等であってもよい。また、棒グラフに折れ線グラフやスプライ曲線グラフを重ねるようにしてもよい。

通常の歩行では、図２７から分かるように、左ステップ時間Ｔ_Lと右ステップ時間Ｔ_Rとの間に大きな相違はない。しかし、右膝を固定した歩行では、図２８から分かるように、右ステップ時間Ｔ_Rが左ステップ時間Ｔ_Lよりも長くなっているのが見て取れる。これは、右膝を金属棒で固定されたことによって右足による前進動作が遅くなっていることを意味しており、図２８中の矢印Ｙ１で示すように、左ステップ時間Ｔ_Lと右ステップ時間Ｔ_Rとの間に比較的大きな差が存在することと、その差の程度が定量化されていることが確認できる。

また、図２８では、左ステップ時間Ｔ_Lと右ステップ時間Ｔ_Rとの差が、同じ足によるステップ時間のばらつきの程度よりもはるかに大きいことや、偶然に生じた大きなステップ時間に起因したものでなく再現性があることなどが一目で確認できる。

また、通常の歩行では、図２７から分かるように、ステップ時間Ｔ_jはほとんど変動していない。しかし、右膝を金属棒で固定した歩行では、図２８から分かるように、歩行期間内の前半に比べて後半の方がステップ時間Ｔ_jが長くなっている（矢印Ｙ２）。これは、歩いているうちに前進動作が遅くなってきていることを示しており、右膝が固定されて歩きにくくなっていることによる疲労を表している。

このように、上記のようなステップ時間の時間変化を表すグラフＧでは、時系列的な順番で交互に並ぶ左ステップ時間Ｔ_Lと右ステップ時間Ｔ_Rとが視覚的に表される。そのため、このようなグラフＧによれば、ステップ時間Ｔ_jの左右差を歩行期間の全体に渡っても局所的な一部の期間においても直感的に把握することができ、歩行の仕方、左右の足の前進動作のバランスなどを容易に理解することができる。また、このようなグラフＧによれば、ステップ時間Ｔ_jのばらつきすなわち安定性を、全体としても左右別々としても直感的に把握することができ、全体として安定な歩行が行えているのか、左右どちらかの足に不安定な前進動作が含まれていないかなどを容易に理解することができる。また、このようなグラフＧによれば、ステップ時間Ｔ_jの時間的な変化を、全体としても左右別々としても直感的に把握することができ、歩行開始からの時間の経過にともない、全体的な歩行の仕方、左右の足の前進動作のバランスなどがどのように変化しているか、疲れによってステップ時間Ｔ_jが伸びていないかなどを容易に理解することができる。そして、操作者１１は、患者１０にこのグラフＧを見せながら歩行運動の評価結果を説明することで、患者１０に自身の歩行運動の評価結果を簡単かつ詳細に理解してもらうことができる。

図２５に戻り、ステップＳ９３では、特徴量演算部７６３が、ステップ時間に係る各種の特徴量を、患者１０の歩行運動の評価用として算出する。

本例では、第１の特徴量として、左ステップ時間Ｔ_Lと右ステップ時間Ｔ_Rとの差異の程度が反映される特徴量を算出する。一般的に、左ステップ時間Ｔ_Lと右ステップ時間Ｔ_Rとの差異の程度が大きいほど、患者１０の左足による前進動作と右足による前進動作とのバランスがより崩れていることが多い。逆に、左ステップ時間Ｔ_Lと右ステップ時間Ｔ_Rとの差異の程度が小さいほど、患者１０の左足による前進動作と右足による前進動作とのバランスが取れていることが多い。したがって、このような第１の特徴量によれば、その大きさに基づいて、患者１０の左足による前進動作と右足による前進動作とのバランスの良さを定量的に評価することができる。例えば、第１の特徴量が、左ステップ時間Ｔ_Lと右ステップ時間Ｔ_Rとの差異の程度が大きいほど大きい値を取る場合、第１の特徴量の値が小さいほど患者１０の前進動作の左右のバランスがよいと評価することができる。

また本例では、第２の特徴量として、ステップ時間のばらつきの程度が反映される特徴量を算出する。一般的に、ステップ時間のばらつきの程度が大きいほど、患者１０の前進動作の安定性がより悪いことが多い。逆に、ステップ時間のばらつきの程度が小さいほど、患者１０の前進動作の安定性がより良いことが多い。したがって、このような第２の特徴量によれば、その大きさに基づいて、患者１０の前進動作の安定性の良さを評価することができる。例えば、第２の特徴量が、ステップ時間のばらつきの程度が大きいほど大きい値を取る場合、第２の特徴量が小さいほど患者１０の前進動作の安定性がより高いと評価することができる。

また本例では、第３の特徴量として、全歩行期間のうち時間的に前寄りのステップ時間と時間的に後寄りのステップ時間との差異の程度が反映される特徴量を算出する。一般的に、時間的に前寄りのステップ時間と時間的に後寄りのステップ時間との差異の程度が大きいほど、患者１０が歩行運動中に疲れてステップ時間が長期化していることが多い。逆に、時間的に前寄りのステップ時間と時間的に後寄りのステップ時間との差異の程度が小さいほど、患者１０が歩行運動中に疲れずにステップ時間が長期化しないことが多い。したがって、このような第３の特徴量によれば、その大きさに基づいて、患者１０の疲れやすさを定量的に評価することができる。例えば、第３の特徴量が、時間的に前寄りのステップ時間と時間的に後寄りのステップ時間との差異の程度が大きいほど大きな値を取る場合、第３の特徴量が小さいほど患者１０はより疲れにくいと評価することができる。

図５に戻り、ステップＳ１０では、表示制御部３１０が、ディスプレイ部３０２を制御して、その画面にグラフＧを表示させる。

以上、本実施形態によれば、操作者１１は、加速度データの解析結果に基づいて患者１０の歩行中の動きを詳細に把握し、例えば効果的な歩行訓練プランを作成することができる。

また、本実施形態によれば、特に、ステップ基準時刻検出部７４の構成及びそのステップ基準時刻検出処理により、患者１０の踵着地のタイミングを高い精度で認識し、一歩ごとの前進動作の基準となる時刻を検出することができる。このような基準時刻を用いることで、患者１０の一歩または複数歩ずつの動作が行われているときの加速度データを抽出したり、その動作の所要時間を測定したりすることができ、患者１０の一歩一歩の動作に着目した解析を行うことができる。故に、本実施形態は、患者１０の歩行運動の客観的な評価に有効である。

また、本実施形態によれば、特に、踵着地左右判別部７５の構成及びその踵着地左右判別処理により、認識された患者１０のそれぞれの踵着地が、左右どちらの足によるものかを高い精度で判別することができる。このような判別結果を用いることで、患者１０の歩行動作における左右のバランスを定量的及び客観的に評価することができ、その結果、左右のバランスをよくするための適切なアドバイスを患者１０に提供することも可能になる。

なお、発明は、上記実施形態に限定されず、その趣旨を逸脱しない範囲において種々の変更が可能である。

例えば、本実施形態では、ステップ基準時刻として、踵着地の時刻を検出しているが、これに限定されず、足が着地した時刻や、つま先蹴りの時刻を検出するようにしてもよい。

また例えば、本実施形態は、発明を人の歩行運動に適用した例であるが、発明を人のその他の移動運動、例えば人の走行運動などにも適用することができる。

また例えば、本実施形態は、上述したように人に取り付けられた加速度センサから得られた加速度データを解析する移動運動解析装置であるが、コンピュータをこのような装置として機能させるためのプログラムもまた発明の実施形態の一つである。

１ 歩行解析システム
１０ 患者
１１ 操作者
２ 加速度センサモジュール
２１ プロセッサ
２２ 加速度センサ
２３ メモリ
２４ 通信Ｉ／Ｆ
２５ バッテリ
２０１ 加速度センサ部
２０２ サンプリング部
２０３ 送信部
３ 歩行解析装置
３１ プロセッサ
３２ ディスプレイ
３３ 操作部
３４ メモリ
３５ 通信Ｉ／Ｆ
３６ バッテリ
３０１ 操作部
３０２ ディスプレイ部
３０３ 患者情報受付部
３０４ 受信部
３０５ 加速度データ取得制御部
３０７ 加速度データ解析部
３１０ 表示制御部
３１２ 記憶部
４１ データベース
４２ 記憶媒体
７１ 加速度成分算出部
７２ 加速度波形生成部
７３ 歩行期間特定部
７４ ステップ基準時刻検出部
７４１ 加速度波形読取部
７４２ 歩行同期波形生成部
７４３ 一歩区間特定部
７４４ 上下加速度反映波形生成部
７４５ 区間最大値特定部
７４６ 踵着地認識部
７５ 踵着地左右判別部
７５１ 加速度波形読取部
７５２ 左右位置波形生成部
７５３ 左右判別部
７６ ステップ時間グラフ生成部
７６１ ステップ時間特定部
７６２ グラフ生成部
７６３ 特徴量演算部

Claims (13)

1. 加速度センサを用いて得られた人の移動運動中における加速度を表すデータに基づいて、前記人の左右方向における加速度の波形を特定する第１の処理と、
   前記第１の処理で得られる加速度の波形を積分処理して速度の波形を生成する第２の処理と、
   前記第２の処理で得られる速度の波形を積分処理して位置の波形を生成する第３の処理と、
   前記第３の処理で得られる位置の波形から、該位置の波形を移動平均処理して得られる成分を除去する第４の処理とを行うことにより、処理済みの位置の波形を生成する第１の生成手段と、
   前記処理済みの位置の波形に基づいて、左足及び右足の各々の一歩に対応する区間または時刻を特定する特定手段とを備えた移動運動解析装置。
2. 前記第１の生成手段は、前記第２の処理または前記第３の処理の後に、該処理で得られる波形における線形成分を除去する第５の処理を行う請求項１に記載の移動運動解析装置。
3. 前記第１の生成手段は、前記第１の処理、前記第２の処理または前記第３の処理の後に、該処理で得られる波形における低周波成分を除去する第６の処理を行う請求項１または請求項２に記載の移動運動解析装置。
4. 前記移動平均処理は、前記人の移動運動における二歩に相当する時間幅をウィンドウ幅とする処理である請求項１から請求項３のいずれか一項に記載の移動運動解析装置。
5. 前記二歩に相当する時間幅は、前記位置の波形に自己相関関数を適用して得られる時間である請求項４に記載の移動運動解析装置。
6. 前記二歩に相当する時間幅は、前記データに基づく前記人の上下方向の加速度と前後方向の加速度との線形結合を表す波形の周期に基づいて決定される請求項４に記載の移動運動解析装置。
7. 前記特定手段は、前記処理済みの位置の波形における上に凸となる波形部分と下に凸となる波形部分とを特定することにより、左足及び右足の各々の一歩に対応する区間または時刻を特定する請求項１から請求項６のいずれか一項に記載の移動運動解析装置。
8. 前記第６の処理は、０．４ヘルツ（Ｈｚ）以下である所定の周波数以下の成分を除去するように設定されたフィルタを用いて、前記低周波成分を除去する処理である請求項３に記載の移動運動解析装置。
9. 前記データに基づいて、前記人の加速度を表す波形を生成する第２の生成手段をさらに備えており、
   前記特定手段は、該生成された波形における左足及び右足の各々の一歩に対応する区間または時刻を特定する請求項１から請求項８のいずれか一項に記載の移動運動解析装置。
10. 前記データに基づいて、前記人の加速度に基づく波形を生成する第２の生成手段をさらに備えており、
    前記特定手段は、該生成された波形における左足及び右足の各々の一歩に対応する区間または時刻を特定する請求項１から請求項８のいずれか一項に記載の移動運動解析装置。
11. 前記特定手段は、前記第２の生成手段により生成された波形における極大値及び極小値のいずれか一方に対応した時刻に対して、左足が着地した時刻と右足が着地した時刻とを対応付ける請求項９または請求項１０に記載の移動運動解析装置。
12. 前記加速度センサと、請求項１から請求項１１のいずれか一項に記載の移動運動解析装置とを備えた移動運動解析システム。
13. コンピュータを請求項１から請求項１１のいずれか一項に記載の移動運動解析装置として機能させるためのプログラム。