(57)【要約】【課題】  動物の歩様把握のために採取された三次元加速度データをわかりやすく視覚的に表現し、「どの方向に、どの程度アンバランスなのか」を含め、歩様を定量的に把握する好適な方法と装置の提供する。【解決手段】  極座標変換で得られた動物動作中の極座標（ｒ、θ、Φ）加速度ベクトルデータのｒ、θ、Φの３成分のうちｒ成分でランキング（順位）表を作成、ランキング上位であるｒ成分を与える注目θ値の前後のｒ成分（加速度ベクトル）データを角度と対応付けて、角度を円周に割り当てた円と該円周と直交するｒ軸によって表す。【選択図】  図６
2008109966.tif 000002
【特許請求の範囲】
【請求項１】
動物に装着された三次元直交座標（ｘ、ｙ、ｚ）加速度センサーから採取される加速度ベクトルデータによって動物動作を極座標（ｒ、θ、Φ）を利用して定量化する方法であって、
極座標（ｒ、θ、Φ）の緯度（仰角θ）90度の方向を重力加速度（上下）方向とし、動物が動作中に加速度センサーから採取される直交座標（ｘ、ｙ、ｚ）加速度ベクトルデータを前記極座標（ｒ、θ、Φ）に変換する工程、前記変換で得られた動物動作中の極座標（ｒ、θ、Φ）加速度ベクトルデータのｒ、θ、Φの３成分のうち、２つの成分を対応付けて表す図を描画する工程を有する動物動作の定量化方法。
【請求項２】
２つの成分を対応付けて表す図が、加速度ベクトルの大きさであるｒ成分のθ分布を、θの角度を円周に割り当てた円と該円周と直交するｒ軸によってｒ成分とθとを対応付けて表す図である、および／または加速度ベクトルの大きさであるｒ成分のΦ分布を、Φの角度を円周に割り当てた円と該円周に直交するｒ軸によってｒ成分とΦとを対応付けて表す図である請求項１の動物動作の定量化方法。
【請求項３】
極座標変換で得られた動物動作中の極座標（ｒ、θ、Φ）加速度ベクトルデータのｒ、θ、Φの３成分のうちｒ成分（加速度ベクトル）の大きさでランキング（順位）表を作成する工程、前記作成されたランキング（順位）表の上位にある加速度ベクトルデータのθ、Φの２成分の値θ0、Φ0に注目し、加速度ベクトルデータのθ成分が、θ0−Δ＜θ0＜θ0＋Δ（Δはπ／１０以下）であるデータ、および／または加速度ベクトルデータのΦ成分が、Φ0−Δ＜Φ0＜Φ0＋Δ（Δはπ／１０以下）であるデータを抽出する工程、注目θ値であるθ0から抽出された加速度ベクトルデータのｒ成分（加速度ベクトル）データをΦと、および／または注目Φ値であるΦ0から抽出された加速度ベクトルデータのｒ成分（加速度ベクトル）データをθと対応付けて表す図を描画する工程をさらに有する請求項１または請求項２の動物動作の定量化方法。
【請求項４】
動物に装着された三次元直交座標（ｘ、ｙ、ｚ）加速度センサーから採取される時系列の加速度ベクトルデータによって動物動作を極座標（ｒ、θ、Φ）を利用して定量化する方法であって、
直交座標（ｘ、ｙ、ｚ）のｚ軸および極座標（ｒ、θ、Φ）の緯度（仰角θ）90度の方向を重力加速度（上下）方向とし、動物が動作中に加速度センサーから採取される時系列の直交座標（ｘ、ｙ、ｚ）加速度ベクトルデータを前記極座標（ｒ、θ、Φ）に変換する工程、前記変換で得られた動物動作中の時系列の極座標（ｒ、θ、Φ）加速度ベクトルデータのｒ成分（加速度ベクトル）の大きさでランキング（順位）表を作成する工程、前記ランキング（順位）表の上位にあるベクトルデータのｒ成分が大きい加速度ベクトルデータのθ、Φの２成分の値θ0、Φ0に注目し、時系列の極座標（ｒ、θ、Φ）加速度ベクトルデータのθ成分が、θ0−Δ＜θ0＜θ0＋Δ（Δはπ／１０以下）である時系列データ、および／または時系列の極座標（ｒ、θ、Φ）加速度ベクトルデータのΦ成分が、Φ0−Δ＜Φ0＜Φ0＋Δ（Δはπ／１０以下）である時系列データをフーリエ変換する工程を有する動物動作の定量化方法。
【請求項５】
移動面に移動方向傾斜角（α）および／または移動方向に直交する方向傾斜角（β）を与える傾斜手段を具備した無限軌道面をもつトレッドミル上を動作する動物に装着された三次元直交座標（ｘ、ｙ、ｚ）加速度センサーから採取される加速度ベクトルデータによって動物動作を極座標（ｒ、θ、Φ）を利用して定量化する方法であって、
極座標（ｒ、θ、Φ）の緯度（仰角θ）90度の方向を重力加速度（上下）方向とし、動作方向傾斜および／または動作方向に直交する方向傾斜がある動作面を動物が動作中に加速度センサーから採取される直交座標（ｘ、ｙ、ｚ）加速度ベクトルデータを前記極座標（ｒ、θ、Φ）に変換する工程、前記変換で得られた動物動作中の極座標（ｒ、θ、Φ）加速度ベクトルデータを移動方向傾斜角（α）および／または移動方向に直交する方向傾斜角（β）をパラメータとして表す図を描画する工程を有する動物動作の定量化方法。
【請求項６】
極座標（ｒ、θ、Φ）の緯度（仰角θ）90度の方向を重力加速度（上下）方向とした次の段階で、かかる緯度（仰角θ）90度の方向を重力加速度（上下）方向から移動方向傾斜角および／または移動方向に直交する方向傾斜角だけシフトする工程、加速度センサーから採取されていた直交座標（ｘ、ｙ、ｚ）加速度ベクトルデータを前記シフト後の極座標（ｒ、θ、Φ）に変換する工程、前記シフト後の変換で得られた動物動作中のシフト後極座標（ｒ、θ、Φ）加速度ベクトルデータを表す図を描画する工程をさらに有する請求項５の動物動作の定量化方法。
【請求項７】
２つ以上の加速度センサーが動物に装着されていて、該２つ以上の加速度センサーから得られる２つ以上の加速度データによる動物動作定量化データの差を出力する工程をさらに有する請求項１から請求項６のいずれかに記載された動物動作の定量化方法。
【請求項８】
動物に装着された三次元直交座標（ｘ、ｙ、ｚ）加速度センサーから採取される加速度ベクトルデータによって動物動作を極座標（ｒ、θ、Φ）を利用して定量化する装置であって、
直交座標（ｘ、ｙ、ｚ）のｚ軸および極座標（ｒ、θ、Φ）の緯度（仰角θ）90度の方向を重力加速度（上下）方向とし、動物が動作中に加速度センサーから採取される直交座標（ｘ、ｙ、ｚ）加速度ベクトルデータを前記極座標（ｒ、θ、Φ）に変換する手段、前記変換で得られた動物動作中の極座標（ｒ、θ、Φ）加速度ベクトルデータのｒ、θ、Φの３成分のうち、２つの成分を対応付けて表す図を描画して表示する手段を有する動物動作の定量化装置。
【請求項９】
極座標変換で得られた動物動作中の極座標（ｒ、θ、Φ）加速度ベクトルデータのｒ、θ、Φの３成分のうちｒ成分（加速度ベクトル）の大きさでランキング（順位）表を作成する手段、前記作成されたランキング（順位）表の上位にある加速度ベクトルデータのθ、Φの２成分の値θ0、Φ0に注目し、加速度ベクトルデータのθ成分が、θ0−Δ＜θ0＜θ0＋Δ（Δはπ／１０以下）であるデータ、および／または加速度ベクトルデータのΦ成分が、Φ0−Δ＜Φ0＜Φ0＋Δ（Δはπ／１０以下）であるデータを抽出する手段、前記抽出手段で注目θ値であるθ0から抽出された加速度ベクトルデータのｒ成分（加速度ベクトル）データをΦと、および／または前記抽出手段で注目Φ値であるΦ0から抽出された加速度ベクトルデータのｒ成分（加速度ベクトル）データをθと対応付けて表す図を描画する手段をさらに有する請求項８の動物動作の定量化装置。
【請求項１０】
動物に装着された三次元直交座標（ｘ、ｙ、ｚ）加速度センサーから採取される時系列の加速度ベクトルデータによって動物動作を極座標（ｒ、θ、Φ）を利用して定量化する装置であって、
直交座標（ｘ、ｙ、ｚ）のｚ軸および極座標（ｒ、θ、Φ）の緯度（仰角θ）90度の方向を重力加速度（上下）方向とし、動物が動作中に加速度センサーから採取される時系列の直交座標（ｘ、ｙ、ｚ）加速度ベクトルデータを前記極座標（ｒ、θ、Φ）に変換する手段、前記変換で得られた動物動作中の時系列の極座標（ｒ、θ、Φ）加速度ベクトルデータのｒ成分（加速度ベクトル）の大きさでランキング（順位）表を作成する手段、前記ランキング（順位）表の上位にあるベクトルデータのｒ成分が大きい加速度ベクトルデータのθ、Φの２成分の値θ0、Φ0に注目し、時系列の極座標（ｒ、θ、Φ）加速度ベクトルデータのθ成分が、θ0−Δ＜θ0＜θ0＋Δ（Δはπ／１０以下）である時系列データ、および／または時系列の極座標（ｒ、θ、Φ）加速度ベクトルデータのΦ成分が、Φ0−Δ＜Φ0＜Φ0＋Δ（Δはπ／１０以下）である時系列データをフーリエ変換する手段を有する動物動作の定量化装置。
【請求項１１】
移動面に移動方向傾斜角（α）および／または移動方向に直交する方向傾斜角（β）を与える傾斜手段を具備した無限軌道面をもつトレッドミル上を動作する動物に装着された三次元直交座標（ｘ、ｙ、ｚ）加速度センサーから採取される加速度ベクトルデータによって動物動作を極座標（ｒ、θ、Φ）を利用して定量化する装置であって、
極座標（ｒ、θ、Φ）の緯度（仰角θ）90度の方向を重力加速度（上下）方向とし、移動方向傾斜および／または移動方向に直交する方向傾斜がある移動面を動物が動作中に加速度センサーから採取される直交座標（ｘ、ｙ、ｚ）加速度ベクトルデータを前記極座標（ｒ、θ、Φ）に変換する手段、前記変換で得られた動物動作中の極座標（ｒ、θ、Φ）加速度ベクトルデータを歩行方向傾斜角（α）および／または歩行方向に直交する方向傾斜角（β）をパラメータとして表す図を描画する手段を有する動物動作の定量化装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、モーションセンサーを装着した動物の動作を分析する一般技術である。限定した用途で効果を具体的に示すのが適切であるので、「動物」を「二足または四足歩行動物」とし、「動作モード」を「歩様」として説明する。他の動作解析に関しては後述する。ここで、「歩様（英文はＧＡＩＴ）」とは、歩行動物の歩き方、歩行の態様を示す動物行動学上の用語である。
【０００２】
歩行動物は種を問わず、二足歩行するヒト、四足歩行の蓄獣を含む。これら動物にモーションセンサーを装着し該センサーから得られるデータから、正常な身体バランスがとれているか、アンバランスで転倒しやすい状態であるか、の弁別、さらに単なるバランス・アンバランスの弁別以上の「どの程度アンバランスなのか」、といったバランスに関する量的指標の抽出は、臨床診断・動物薬効試験などで必要である。
【背景技術】
【０００３】
ヒト歩様を把握することは、転倒防止や歩行障害治療のリハビリで重要視され、力学センサー・画像センサー・加速度センサーでデータを採取して解析されてきた。これらセンサーのうち、フォースプレート（床反力の測定装置）などの力学センサー方式、モーションキャプチャーなどの動画像解析を組み合せた画像センサー方法（特許文献１参照）は、コスト・パフォーマンスで加速度センサー方式に対し、やや優位性を欠いている。特に画像方式は関節部位にマーカ装着する手間と画像記録や解析に高価な装置を要する問題がある。そこで、データ採取の容易性および経済性から、加速度センサーを用いた方式（特許文献２から６、非特許文献１から５参照）に注目する。
【０００４】
特許文献２では、ヒトの脛骨および踵に加速度センサーを装着し歩行中に観測される複数の垂直方向時刻のピークに注目した歩様解析を開示している。さらに同文献にて、加速度の時系列データをフーリエ変換して周波数解析することが開示されている。
【０００５】
特許文献３では、パーキンソン患者の指に加速度センサーを装着し、指同士の接触インターバル時間を解析する動作解析手法が開示されている。
【０００６】
特許文献４では、ヒトの腹部に加速度センサーを装着し、採取されたデータと背屈力との関係で歩様解析する手法が開示されている。
【０００７】
特許文献５では、四足歩行動物の牛に加速度センサーを装着し、加速度ベクトルの終点の移動を示すデータを採取し、該データを牛の前後・左右・上下のうちの２つの座標軸を組み合わせたリサージュ図形で視覚化し、歩様を定量化する方法が開示されている。これは疾病乳牛の治療前後の変化を定量化することで治療法でもある。
【０００８】
こういった加速度データ解析について加速度センサーの原理からくる問題がある。加速度センサーの加速度検知は、加速度による分子構造のズレのエネルギーを電気信号に変換する素子によるものであり、かかる変換素子を三次元直交座標（ｘ、ｙ、ｚ）方向に配設していることから、センサーデータは三次元直交座標（ｘ、ｙ、ｚ）の加速度データである。ゆえに、ベクトル方向は直接的には得られない。
【０００９】
特許文献６に、連続採取された歩行状態の特定期間内の加速度ベクトルデータを用いて、加速度センサ座標系における重力加速度方向および身体軸を自動修正する機能をもつ解析法が開示されている。しかし特許文献６の技術では、歩行状態に異常があると適切な座標方向・座標軸の修正ができない。
【００１０】
なお、２つの生体データの解析にリサージュ図形を利用することは公知である。具体的には、特許文献７に第１ の測定信号をＸ 軸に、第２ の測定信号をＹ 軸にして、生体データのリサージュ波形を描くことでデータ解析する手法が開示されている。さらに本件と同一発明者らによる特許文献８には、歩行動物に装着された三次元直交座標（ｘ、ｙ、ｚ）加速度センサーから採取される加速度ベクトルの３つの成分中の２成分のデータのリサージュ波形を描くことで歩行動物の歩様を解析する手法が具体的に記載されている（図２・図３参照）。
【００１１】
非特許文献１から６は、加速度による歩様バランス定量化の学術的アプローチを開示している。特に、非特許文献６で３方向加速度センサの歩行データをフーリエ変換した周波数解析が開示されている。しかし、三次元直交座標（ｘ、ｙ、ｚ）を用いた分析であって、、ベクトル方向は直接的には得られない。その他の座標系での分析は開示されていない。
【００１２】
以上例示した個々の公知技術およびこれらを組み合わせた技術では、歩行中に正常な身体バランスがとれている状態とアンバランスで転倒しやすい状態の「弁別」、その弁別に用いる具体的かつ視覚的にわかりやすい量的指標、特に、どの方向に、どの程度アンバランスなのか、というバランスに関する具体的かつ視覚的にわかりやすい量的指標については、まだ改善の余地がある。
【００１３】
どの方向に、どの程度アンバランスなのか、というバランスに関する量的指標を三次元直交座標（ｘ、ｙ、ｚ）以外の座標系で具体的に分析した先行文献は見当たらない。本発明では極座標（ｒ、θ、Φ）を用いるが、それ自体は公知である（図１参照）。直交座標（ｘ、ｙ、ｚ）と極座標（ｒ、θ、Φ）とは、「ｘ＝ｒ ｓｉｎθ ｃｏｓΦ、ｙ＝ｒ ｓｉｎθ ｓｉｎΦ、ｚ＝ｒ ｃｏｓθ」の関係式で容易に相互変換される。
【００１４】
【特許文献１】国際公開WO05-96939号公報「リハビリ支援用計測システム、及びリハビリ支援用計測方法」財団法人新産業創造研究機構
【特許文献２】特開2004-261525号公報「変形性膝関節症の判定方法および判定装置 」マイクロストーン株式会社
【特許文献３】特開2005-152053号公報「動作解析装置およびその利用」独立行政法人科学技術振興機構
【特許文献４】特開2006-87735号公報 「歩行解析装置」アイシン精機株式会社
【特許文献５】特開2006-218122号公報「脚状態診断システムとその診断方法」財団法人新産業創造研究機構
【特許文献６】特開2006−175206号公報「身体状態検出装置、その検出方法及び検出プログラム」独立行政法人産業技術総合研究所
【特許文献７】特開2005−245574号公報「信号処理方法及びそれを用いたパルスフォトメータ」日本光電工業株式会社
【特許文献８】特願2005-317524号「診断システム」財団法人新産業創造研究機構
【非特許文献１】神先秀人、外４名、「機器を使用した歩行の評価」、理学療法学、2003年6月20日、第30巻、第4号、p-258、図7-8
【非特許文献２】小林哲平、外3名「加速度センサを用いた運動学的歩行分析システム一般関節疾患の術後リハビリにおけるWalk―Mite有効性評価への適用」計測自動制御学会論文集、vol.42、No.5、p567 (2006)
【非特許文献３】松原正明、外5名「三次元加速度センサー計測装置による片側変形性股関節症術前・術後の歩行の検討」 リハビリテーション医学, vol.42, suppl,pp.S152 (2005) 第42回日本リハビリテーション医学会学術集会
【非特許文献４】和田義明、外5名「三次元加速度センサーを用いた計測装置による歩行の検討」 リハビリテーション医学, vol.42, suppl, pp.S335 (2005)第42回日本リハビリテーション医学会学術集会
【非特許文献５】下釜みずき、三宅美博, "共創型歩行介助システムへの加速度センサ導入による多様な歩行障害への対応," 第4回SICEシステムインテグレーション部門講演会講演論文集(SI2003), pp.175-176 (2003)
【非特許文献６】黄  健、関根正樹、外3名「３方向加速度センサを用いた歩行計測と周波数解析」日本機械学会  第76期全国大会講演会講演論文集(II)，No.98-3(1998-10), pp.285-286.
【発明の開示】
【発明が解決しようとする課題】
【００１５】
加速度センサーをもちいて歩様を定量化する多くの試みがなされているが、歩様バランスに注目し、そのバランス判定の明確な手続きが確立されるに至ってはいない。とりわけ、「どの方向に、どの程度アンバランスなのか」、というバランスに関する具体的かつ視覚的にわかりやすい量的指標については、技術的に未完成でがある。本発明は、三次元加速度データをわかりやすく視覚的に表現し、「どの方向に、どの程度アンバランスなのか」を含め、歩様を定量的に把握する好適な方法と装置を提案するものである。
【課題を解決するための手段】
【００１６】
本発明の本質は、加速度の方向を分析するには直交座標（ｘ、ｙ、ｚ）よりも極座標（ｒ、θ、Φ）を利用したほうが、より適切であることによる。ただし単に座標変換しただけではなく、データをわかりやすく視覚的に表現する図示方法とその手順を提供している。これを以下に説明する。
【００１７】
本発明は（請求項１）、動物に装着された三次元直交座標（ｘ、ｙ、ｚ）加速度センサーから採取される加速度ベクトルデータによって動物動作を極座標（ｒ、θ、Φ）を利用して定量化する方法であって、直交座標（ｘ、ｙ、ｚ）のｚ軸および極座標（ｒ、θ、Φ）の緯度（仰角θ）90度の方向を重力加速度（上下）方向とし、動物が動作中に加速度センサーから採取される直交座標（ｘ、ｙ、ｚ）加速度ベクトルデータを前記極座標（ｒ、θ、Φ）に変換する工程、前記変換で得られた動物動作中の極座標（ｒ、θ、Φ）加速度ベクトルデータのｒ、θ、Φの３成分のうち、２つの成分を対応付けて表す図を描画する工程を有する動物動作の定量化方法である（図１４のフローチャート参照）。
【００１８】
また本発明は（請求項２）、「２つの成分を対応付けて表す図」が、加速度ベクトルの大きさであるｒ成分のθ分布を、θの角度を円周に割り当てた円と該円周と直交するｒ軸によってｒ成分とθとを対応付けて表す図である、および／または加速度ベクトルの大きさであるｒ成分のΦ分布を、Φの角度を円周に割り当てた円と該円周に直交するｒ軸によってｒ成分とΦとを対応付けて表す図とした方法である。
【００１９】
図４が、本発明の極座標（ｒ、θ、Φ）であって、図２・図３同様に「側面」「正面（背面）」「平面」図を表示している。角度そのものの情報が表示されることに特徴がある。真上・真下・真前・真後ろ・真横（左右）を示すための角度設定の代表値を示す。すなわち、真上はθ=90度（π/2）、真下はθ=-90度（-π/2）、真前はΦ=0度、真後ろはΦ=180度（π）、真横（左右）はΦ=90,270度（π/2,3π/2）で表される。
【００２０】
図５が、本発明の加速度データ表示の例図であって、図５上図は、θの角度を円周に割り当てた円１１に１１の円周と直交するｒ軸３を記した図、図５下図は、Φの角度を円周に割り当てた円１２に１２の円周と直交するｒ軸３を記した図である。ｒ軸３は、３つの角度例（θ1θ2θ3、Φ1Φ2Φ3）のみ記し、該軸上に加速度データ例（○、●、△、▲、▽、▲、□、■）を記入した。ｒ軸３は略してもかまわない。また、円１１円１２の円周がｒ軸のゼロライン（ｒ=0）であるのが好適であるがその限定はない。
【００２１】
図５上図を詳しく見ると、上図の左半面はΦ=0、右半面はΦ=πのｒ軸データを抽出し記載したものである（図４側面図参照）。同様に、下図の右図はθ=0のｒ軸データを抽出し記載したものである（図４側面図・正面図参照）。ここで注意を要することは、ｒ軸データを抽出する際の角度成分が厳密な0、π（0.00000・・・・、3.14159265・・・・）の角度成分をもつものを抽出するということではなく、任意の数値幅をもつものを抽出してかまわない、ということである。
【００２２】
たとえばΦ=0の加速度データは、真前方向の加速度成分をもつものであるが、歩様解析において厳密な真前のみを注目することは無意味である。つまり、数度の角度以内（−λ＜Φ＜λ：λは数度）の範囲のデータがすべて真前方向であるとしてかまわない。同様に、θ=0の加速度データは、水平方向の加速度成分をもつものであるが、歩様解析において厳密な水平のみを注目することは無意味である。つまり、数度の角度以内（−Δ＜θ＜Δ：Δはπ／１０以下）の範囲のデータがすべて水平方向であるとしてかまわない。厳密な0、π（0.00000・・・・、3.14159265・・・・）の抽出とするとわずかなデータしか条件にフィットしないが、その分角度の分析分解能はあがる。逆に、Δの設定が歩様解析の分解能の条件であるので、Δは歩様定量化の条件として記録されるべきものである。
【００２３】
同心円状に図５上図（左半面Φ=0、右半面はΦ=180度（π））と下図（θ=0）のｒ軸データが記載された図を重ね描画してもよい（図示略）。また、同心円状かつ半円、３分の１円、４分の1円等々の組み合わせで、注目すべきΦ=Φ1、Φ=Φ2、Φ=Φ3、Φ=Φ4等々のｒ軸データが記載された円図または部分円図を組み合わせ描画してもよい（図示略）。さらにΦ=Φi、とθ=θjの任意の注目伏角仰角の角度のｒ軸データが記載された円図または部分円図を組み合わせ描画してもよい（図示略）。ここでΦ=Φi、とθ=θjは、正しくはΦi−Δ＜Φ＜Φi＋Δ、θj−Δ＜θ＜θj＋Δであって、Δはπ／１０以下が好適である。
【００２４】
注目すべきΦ=Φ1、Φ=Φ2、Φ=Φ3、Φ=Φ4等々を如何にして選ぶか、であるが、以下のようにすればよい。すなわち、本発明では（請求項３）、極座標変換で得られた動物動作中の極座標（ｒ、θ、Φ）加速度ベクトルデータのｒ、θ、Φの３成分のうちｒ成分（加速度ベクトル）の大きさでランキング（順位）表を作成する工程、前記作成されたランキング（順位）表の上位にある加速度ベクトルデータのθ、Φの２成分の値θ0、Φ0に注目し、加速度ベクトルデータのθ成分が、θ0−Δ＜θ0＜θ0＋Δ（Δはπ／１０以下）であるデータ、および／または加速度ベクトルデータのΦ成分が、Φ0−Δ＜Φ0＜Φ0＋Δ（Δはπ／１０以下）であるデータを抽出する工程、注目θ値であるθ0から抽出された加速度ベクトルデータのｒ成分（加速度ベクトル）データをΦと、および／または注目Φ値であるΦ0から抽出された加速度ベクトルデータのｒ成分（加速度ベクトル）データをθと対応付けて表す図を描画する工程をさらに有すると好適である（（図１５のフローチャート参照）。
【００２５】
ここでΔは角度の小さい値であって、分析の精度または分解能にあたる値である。この値が大きいほど大雑把な分析、小さいほど角度について小さな差異にも敏感な分析になる一方、より多くのデータが必要になる。
【００２６】
繰り返しになるが、Δは角度切り出しの角度範囲であって、分析の分解能である。アンバランス角度方向の考察では、プラスマイナス１度といったこまかい角度の差は無意味であるので、角度切り出しは、５度や１０度の角度で切り出して問題ない。すなわち、０≦Φ＜５、５≦Φ＜１０・・・・・８５≦Φ≦９０度といった範囲で切り分ければ十分である。逆に、ある歩様行動個体と他の類似した歩様行動をとる個体との歩様を厳密に弁別したければ、十分に多くのデータを長い時間をかけて採取して、上記の角度区分を、０≦Φ＜１、２≦Φ＜３・・・・・８９≦Φ≦９０度といった細かい範囲で切り分けて解析することも可能ではある。
【００２７】
図６の上図は、左半面Φ=0、右半面Φ=πの側面における加速度の大きさ（ｒ成分）とθとを対応付けて表す図であり、図６の左下図は、左半面Φ=90度（π/2）、右半面Φ=270度（3π/2）の側面における加速度の大きさ（ｒ成分）とθとを対応付けて表す図である。図６のデータの歩様モードは、θ=0の水平面でそれぞれ真前（Φ=0）、真後ろ（Φ=180度（π））への加速度、ならびに真横（左右）への加速度（Φ=90,270度（π/2,3π/2））の加速度が大きいモードであってバランスが良好であることがわかる。
【００２８】
図７の左上図は、左半面Φ=Φ1  (０＜Φ1＜π/2）、右半面Φ=Φ3（π＜Φ3＜3π/2）の側面における加速度の大きさ（ｒ成分）とθとを対応付けて表す図であり、図７左下図は、左半面Φ=Φ2  (π/2＜Φ2＜π）、右半面Φ=Φ4（3π/2＜Φ4＜2π）の側面における加速度の大きさ（ｒ成分）とθとを対応付けて表す図である。図７のデータの歩様モードは、θ=π/2近傍で、Φ=Φ3（左斜め後方）とΦ=Φ2（右斜め後方）の加速度が大きく、ややアンバランスであることがわかる。
【００２９】
図８の左上図は、上半面θ=θ3（0＜θ3＜π/2）、下半面θ=θ1  (-π/2＜θ1＜0）の斜側面における加速度の大きさ（ｒ成分）とΦとを、図８左中図は、θ=0の水平面における加速度の大きさ（ｒ成分）とΦとを、図８左下図は、上半面θ=θ2  （0＜θ3＜π/2）、下半面θ=θ4(-π/2＜θ1＜0）の斜側面における加速度の大きさ（ｒ成分）とΦとを、それぞれ対応付けて表す図である。
【００３０】
図９は、図８に加速度データをプロットした例図であって、本図のデータの歩様モードは左上図よりθ=π/4近傍の加速度が大きく、ややアンバランスであることがわかる。
【発明の効果】
【００３１】
本発明によって、加速度センサーで採取された動物動作データが、従来よりもわかりやすく視覚化（グラフ化）され、歩様を定量的に把握できるようになる。歩様定量化について、より明確な手続きが確立され、様々な動作データの統一的評価ができるようになる。特に「どの方向に」歩様がアンバランスであるか、については、θとΦという角度で明確化され、それらとベクトルの大きさｒ値が対応しているので、バランス良好データのｒ値とバランス不良データのｒ値とを比較することで「どの程度アンバランスなのか」、も定量把握できる。
【発明を実施するための最良の形態】
【００３２】
発明実施の際に好適な他の工夫を説明する。本発明にても、従来からなされているように時系列データからパワスペクトルを得て定量化するのが好適である。すなわち（請求項４）、動物が動作中に加速度センサーから採取される時系列の直交座標（ｘ、ｙ、ｚ）加速度ベクトルデータを前記極座標（ｒ、θ、Φ）に変換する工程、前記変換で得られた動物動作中の時系列の極座標（ｒ、θ、Φ）加速度ベクトルデータのｒ成分（加速度ベクトル）の大きさでランキング（順位）表を作成する工程、前記ランキング（順位）表の上位にあるベクトルデータのｒ成分が大きい加速度ベクトルデータのθ、Φの２成分の値θ0、Φ0に注目し、時系列の極座標（ｒ、θ、Φ）加速度ベクトルデータのθ成分が、θ0−Δ＜θ0＜θ0＋Δ（Δはπ／１０以下）である時系列データ、および／または時系列の極座標（ｒ、θ、Φ）加速度ベクトルデータのΦ成分が、Φ0−Δ＜Φ0＜Φ0＋Δ（Δはπ／１０以下）である時系列データをフーリエ変換する工程を実施することも好適である。
【００３３】
図１０の上図は、従来の直交座標（ｘ、ｙ、ｚ）で表された時系列の加速度ベクトルデータをフーリエ変換した対象例図。図１０の下図は上図と同じ時系列の加速度ベクトルデータを極座標（ｒ、θ、Φ）に変換し、変換で得られた動物動作中の極座標（ｒ、θ、Φ）加速度ベクトルデータのｒ成分（加速度ベクトル）の大きさでランキング（順位）表を作成し、作成されたランキング（順位）表の上位にある加速度ベクトルデータのθ、Φの２成分の値に注目し、かかる注目θ値または注目Φ値の前後、すなわち、（注目θ値−Δ）＜（注目θ値）＜（注目θ値＋Δ）の範囲、または（注目Φ値−Δ）＜（注目Φ値）＜（注目Φ値＋Δ）の範囲にある加速度ベクトルデータを抽出し、さらに抽出された加速度ベクトルデータのｒ成分（加速度ベクトル）時系列をフーリエ変換した図である。同じデータでありながら、本発明による図１０下図の方がパワスペクトルの大きなピークを検知しており、より適切な歩様の把握と定量化がなされているのがわかる。
【００３４】
また、本発明では（請求項５）、傾斜したトレッドミルを利用してもよい。すなわち、移動面に移動方向傾斜角（α）および／または移動方向に直交する方向傾斜角（β）を与える傾斜手段を具備した無限軌道面をもつトレッドミル上を動作する動物に装着された三次元直交座標（ｘ、ｙ、ｚ）加速度センサーから採取される加速度ベクトルデータによって動物動作を極座標（ｒ、θ、Φ）を利用して定量化する方法であって、極座標（ｒ、θ、Φ）の緯度（仰角θ）90度の方向を重力加速度（上下）方向とし、動作方向傾斜および／または動作方向に直交する方向傾斜がある動作面を動物が動作中に加速度センサーから採取される直交座標（ｘ、ｙ、ｚ）加速度ベクトルデータを前記極座標（ｒ、θ、Φ）に変換する工程、前記変換で得られた動物動作中の極座標（ｒ、θ、Φ）加速度ベクトルデータを移動方向傾斜角（α）および／または移動方向に直交する方向傾斜角（β）をパラメータとして表す図を描画する工程を実施することも好適である。
【００３５】
図１１が、歩行面に歩行方向傾斜角（α）を与える傾斜手段を具備した無限軌道回転移動面をもつトレッドミルの説明図、図１２が、歩行面に直交する方向傾斜角（β）を与える傾斜手段を具備した無限軌道回転移動面をもつトレッドミルの説明図である。これらを利用してデータを採取し、歩行面の歩行方向傾斜角（α）であって上り勾配角度αup、歩行面の歩行方向傾斜角（α）であって下り勾配角度αdown、歩行面の歩行方向に直交する方向傾斜角（β）であって右上がり（右下がり）左右傾斜（カント）角度βright、歩行面の歩行方向に直交する方向傾斜角（β）であって左上がり（左下がり）左右傾斜（カント）角度βleftと採取データを関連付けた解析をおこなってもよい。
【００３６】
すなわち（請求項６）、極座標（ｒ、θ、Φ）の緯度（仰角θ）90度の方向を重力加速度（上下）方向とした次の段階で、かかる緯度（仰角θ）90度の方向を重力加速度（上下）方向から移動方向傾斜角および／または移動方向に直交する方向傾斜角だけシフトする工程、加速度センサーから採取されていた直交座標（ｘ、ｙ、ｚ）加速度ベクトルデータを前記シフト後の極座標（ｒ、θ、Φ）に変換する工程、前記シフト後の変換で得られた動物動作中のシフト後極座標（ｒ、θ、Φ）加速度ベクトルデータを表す図を描画する工程を実施してもよい（図１３参照）。
【００３７】
さらに、複数の加速度センサーを動物に装着し、それらの加速度センサーの信号から、センサ装着位置の中間または２つのセンサー装着位置の延長線上にある加速度の極小点（不動点）の推定、も可能であろう。すなわち（請求項７）、２つ以上の加速度センサーが動物に装着されていて、該２つ以上の加速度センサーから得られる２つ以上の加速度データによる動物動作定量化データの差を出力する工程をさらに実施してもよい。ただし加速度の極小点（不動点）の推定のための方法・アルゴリズムは現時点では明確でない。
【００３８】
本発明の実施装置は、基本的にソフトウェアであって、ここまでに説明した工程をなすアルゴリズムを実行するプログラムをパソコンなどにインストールすればよい。それらプログラムが各工程を実施する手段となる。請求項８の装置は請求項１、請求項９の装置は請求項３、請求項１０の装置は請求項４、請求項１１の装置は請求項５の工程を手段と置換したものであるので説明は省く。
【００３９】
本発明は、動物動作を分析する一般技術であり、例として「動物」を「二足または四足歩行動物」とし、「動作モード」を「歩様」として説明した。しかしこれらに限定されることはなく、二足・四足歩行動物以外の動物、あるいはロボットを含む運動メカニズムに加速度センサーを装着して適用できる。分析対象の「動作モード」も無制限であって、たとえばヒトのスポーツ動作では、バッティング・スイング動作、水泳のストローク動作、臨床用途では、四肢の運動障害の動作分析に利用できる。
【図面の簡単な説明】
【００４０】
【図１】極座標（ｒ、θ、Φ）と直交座標（ｘ、ｙ、ｚ）および両座標の座標変換式
【図２】三次元直交座標（ｘ、ｙ、ｚ）を利用して歩行動物の歩様バランスを定量化するための図例であって、加速度ベクトルの３成分中の２成分のデータの対応を描くための歩行動物「側面」「正面（背面）」「平面」図。
【図３】図２の歩行動物「側面」「正面（背面）」「平面」図に加速度センサーから採取される加速度ベクトルの終点を描いた例であって、３成分中の２成分がそれぞれ縦軸と横軸に反映されている。いわゆるリサージュ図形であってこの図自体にベクトルの角度情報はない。
【図４】本発明の極座標（ｒ、θ、Φ）であって、図２・図３同様に「側面」「正面（背面）」「平面」図を表示している。角度そのものの情報が表示されることに特徴点がある。真上・真下・真前・真後ろ・真横（左右）を示すための角度設定の代表値を示す。すなわち、真上はθ=90度（π/2）、真下はθ=-90度（-π/2）、真前はΦ=0度、真後ろはΦ=180度（π）、真横（左右）はΦ=90,270度（π/2,3π/2）で表される。
【図５】本発明の加速度データ表示の例図であって、上図は、θの角度を円周に割り当てた円１１に１１の円周と直交するｒ軸３を記した図、下図は、Φの角度を円周に割り当てた円１２に１２の円周と直交するｒ軸３を記した図である。ｒ軸３は、３つの角度例（θ1θ2θ3、Φ1Φ2Φ3）のみ記し、該軸上に加速度データ例（○、●、△、▲、▽、▲、□、■）を記入した。
【図６】左上図は、左半面Φ=0、右半面Φ=πの側面における加速度の大きさ（ｒ成分）とθとを対応付けて表す図であり、左下図は、左半面Φ=90度（π/2）、右半面Φ=270度（3π/2）の側面における加速度の大きさ（ｒ成分）とθとを対応付けて表す図である。本図のデータの歩様モードは、θ=0の水平面でそれぞれ真前（Φ=0）、真後ろ（Φ=180度（π））への加速度、ならびに真横（左右）への加速度（Φ=90,270度（π/2,3π/2））の加速度が大きいモードであってバランスが良好であることがわかる。
【図７】左上図は、左半面Φ=Φ1  (０＜Φ1＜π/2）、右半面Φ=Φ3（π＜Φ3＜3π/2）の側面における加速度の大きさ（ｒ成分）とθとを対応付けて表す図であり、左下図は、左半面Φ=Φ2  (π/2＜Φ2＜π）、右半面Φ=Φ4（3π/2＜Φ4＜2π）の側面における加速度の大きさ（ｒ成分）とθとを対応付けて表す図である。本図のデータの歩様モードは、θ=π/2近傍で、Φ=Φ3（左斜め後方）とΦ=Φ2（右斜め後方）の加速度が大きく、ややアンバランスであることがわかる。
【図８】左上図は、上半面θ=θ3（0＜θ3＜π/2）、下半面θ=θ1  (-π/2＜θ1＜0）の斜側面における加速度の大きさ（ｒ成分）とΦとを、左中図は、θ=0の水平面における加速度の大きさ（ｒ成分）とΦとを、左下面は、上半面θ=θ2  （0＜θ3＜π/2）、下半面θ=θ4(-π/2＜θ1＜0）の斜側面における加速度の大きさ（ｒ成分）とΦとを、それぞれ対応付けて表す図である。
【図９】図８に加速度データをプロットした例図であって、本図のデータの歩様モードは左上図よりθ=π/4近傍の加速度が大きく、ややアンバランスであることがわかる。
【図１０】上図は従来の直交座標（ｘ、ｙ、ｚ）で表された時系列の加速度ベクトルデータをフーリエ変換した対象例図。下図は上図と同じ時系列の加速度ベクトルデータを極座標（ｒ、θ、Φ）に変換し、変換で得られた動物動作中の極座標（ｒ、θ、Φ）加速度ベクトルデータのｒ成分（加速度ベクトル）の大きさでランキング（順位）表を作成し、作成されたランキング（順位）表の上位にある加速度ベクトルデータのθ、Φの２成分の値に注目し、かかる注目θ値または注目Φ値の前後、すなわち、（注目θ値−Δ）＜（注目θ値）＜（注目θ値＋Δ）の範囲、または（注目Φ値−Δ）＜（注目Φ値）＜（注目Φ値＋Δ）の範囲にある加速度ベクトルデータを抽出し、さらに抽出された加速度ベクトルデータのｒ成分（加速度ベクトル）時系列をフーリエ変換した図である。同じデータでありながら、本発明による下図の方がパワスペクトルの大きなピークを検知しており、適切な歩様定量化がなされているのがわかる。
【図１１】歩行面に歩行方向傾斜角（α）を与える傾斜手段を具備した無限軌道回転移動面をもつトレッドミルの説明図
【図１２】歩行面に直交する方向傾斜角（β）を与える傾斜手段を具備した無限軌道回転移動面をもつトレッドミルの説明図
【図１３】極座標（ｒ、θ、Φ）の緯度（仰角θ）90度の方向を重力加速度（上下）方向とした次の段階で、かかる緯度（仰角θ）90度の方向を重力加速度（上下）方向から移動方向傾斜角（α）および／または移動方向に直交する方向傾斜角（β）だけシフトすることの説明図
【図１４】請求項１のフローチャートであって図中の「サブフローＡ」は請求項３のフローチャートである。
【図１５】請求項３のフローチャート
【符号の説明】
【００４１】
１１  θの角度を円周に割り当てた円
１２  Φの角度を円周に割り当てた円
３  １１または１２の円周と直交するｒ軸
○、●、△、▲、▽、▲、□、■は加速度ベクトルのデータ
Ｈ  歩行中の歩行動物（たとえばヒト）
αup    歩行面の歩行方向傾斜角（α）であって上り勾配角度
αdown  歩行面の歩行方向傾斜角（α）であって下り勾配角度
βright   歩行面の歩行方向に直交する方向傾斜角（β）であって右上がり（右下がり）左右傾斜（カント）角度
βleft     歩行面の歩行方向に直交する方向傾斜角（β）であって左上がり（左下がり）左右傾斜（カント）角度  
Ｔｒ  無限軌道回転歩行面を有するトレッドミル
ｔup    歩行面に歩行方向傾斜角αupがある時間
ｔdown  歩行面に歩行方向傾斜角αdownがある時間
【図１】
2008109966.tif 000003
【図２】
2008109966.tif 000004
【図３】
2008109966.tif 000005
【図４】
2008109966.tif 000006
【図５】
2008109966.tif 000007
【図６】
2008109966.tif 000008
【図７】
2008109966.tif 000009
【図８】
2008109966.tif 000010
【図９】
2008109966.tif 000011
【図１０】
2008109966.tif 000012
【図１１】
2008109966.tif 000013
【図１２】
2008109966.tif 000014
【図１３】
2008109966.tif 000015
【図１４】
2008109966.tif 000016
【図１５】
2008109966.tif 000017