目次第1章

緒……………………………………………………………… 2第2章

歩行解析法

言

………………………………………………………… 32.1

測定システム……………………………………………………

32.2

足の軌跡の算出…………………………………………………

32.3

実験方法

…………………………………………………………

5第3章

実験結果

…………………………………………………………… 63.1

机上実験…………………………………………………………

63.2

歩行実験

…………………………………………………………

6第4章

考察

………………………………………………………………… 74.1

机上実験…………………………………………………………

74.2

歩行実験

…………………………………………………………

7

第5章

結言

………………………………………………………………… 81

第 1 章 緒言
今後，我が国においては，加齢に伴う動作能力の低下により日常生活動作に支障をきたす高齢者の増加が考えられる．そこで，高齢者の日常生活動作の特徴を明らかにすることが重要な課題となっている．動作能力低下のメカニズムを理解し，得られた知見をもとに適切な施策を講ずることで，より効果的な支援を実現することが可能だからである．特に移動に関する能力は高齢者の社会生活の自立にとって重要である．歩行時に転倒する危険性は，年齢とともに大きくなる．歩行時に転倒することで股関節骨折などの重篤な障害に繋がり，生活が困難になる．また，現在得られているデータは医学的，リハビリテーション的なものが多く，日常生活動作に即した動的なデータが少ない．そのため，高齢者が利用する製品や生活環境の改善には活用されにくいというのが実情である．
高齢者の歩行特性を明らかにすることにおいて,日常生活動作の動的なデータが不可欠である．現在は，主にディジタルビデオカメラや高速度カメラから得られる映像情報から三次元位置座標を算出する DLT(Direct Linear Transformation)法が用いられている．しかし，DLT 法を用いた計測は高価かつ大規模である実験室的環境が必要である．一方，慣性センサを用いた方法は，低コストかつ簡易で計測場所が制限されない計測が可能である．
慣性センサは装置の小型化，低価格化，高精度化が進んでいる．また，身体部位に取り付けることにより，姿勢を簡易に計測可能である事から，スポーツ工学，医用生体工学分野等において，慣性センサを活用した運動計測法が提案されている[1][2]．慣性センサから得られる加速度と角速度の計測情報を組み合わせることで三次元空間におけるセンサの位置と向きの時間変化を構成することができる．その結果，歩行評価のための歩幅や歩行速度，姿勢情報等を得ることができる．
そこで本研究では，足に装着した慣性センサによって，歩行中の足の軌跡を測定するシステム開発を目的とし，試用実験を行った．

2

第 2 章 歩行解析法
2. 1 測定システム
本研究で用いた計測システムを第 2.1-1 図に示す．本計測システムは,3 軸加速度センサ(InvenSense，MPU-9250，測定範囲：±2[g]）,3 軸ジャイロセンサ(InvenSense,MPU-9250,測定範囲：±250[deg/sec]), 3 軸地磁気センサ（旭化成エレクトロニクス,AK8963,測定範囲：±4800[μT]）が搭載された 9 軸センサモジュールを用いており,3 軸の加速度,角速度,地磁気を計測可能である．本研究では,9 軸のうち 6 軸（加速度,ジャイロ）を使用している．センサは第 2.1-2 図のように装着する．各センサからの情報を取得し,PC へ送信するために,Arduino UNO R2（第 2.1-3 図）を使用している．Arduino は記録のほかに,センサに必要な 3.3V 電源を供給する役割も担っている．通信方法には I2C 通信という方法を用いている．I2C（Inter-Integrated Circuit）とはフィリップス社が提唱した周辺デバイスとのシリアル通信の方式で,シリアル EEPROM メモリ等との高速通信を実現する方式である．本システ
ムは,身体部位に直接取り付けることも出来るが,本研究ではあらかじめ履物に取り付けておき,3 次元の角速度,加速度を計測した．
2. 2 足の軌跡の計算
本研究では, Kitagawa らの方法[3]を参考に測定されたセンサデータに基づいて足の軌跡をオフラインで計算を行った．最初に足底接地時（足底面がが地面に接触している期間）に角速度の大きさが 35[rad/s]より小さくなると定義した．また,足が動作しているかどうか判断するためのしきい値として使用して,角速度がしきい値を超え続けた場合は動作中とした．
歩行中の軌道は,2 つの足底接地の間の並進加速度を時間積分することで計算した．測定された加速度ベクトルは,時間変化するセンサ固定座標系（𝚺S）（第 2.2-1 図）によって表される．しかし,測定された加速度ベクトルは基点座標系内に配列されなければならない．したがって,足底接地時のセンサ座標系の向き（𝚺F）を基点座標系とした．センサは基本的に各足底接地時に静止するので,足底接地時にドリフト誤差を取り除くために回転行列内の値をすべて初期値になるようにした．
測定された加速度ベクトルと基点座標系に対するセンサの向きを表す行列をそれぞれ𝑆𝐚(𝑡)と[𝐢 𝐣 𝐤] = 𝐹𝑹𝑆 とする．上付き文字は,ベクトルが表される座標系を示している．𝚺F座
標系 F𝐚(𝑡)で表される加速度ベクトルは,式(1)に示す．
𝐚(𝑡) = 𝐹 𝑹𝑆 𝑆𝐚(𝑡)F(1)

本研究では,足底接地時からのジャイロセンサの信号（𝛚,センサの角速度ベクトル）を時間積分することにより姿勢を求めることが出来る．歩行中の基点座標系に対する姿勢の時間変化を推定した 𝐹𝑹𝑆 (𝑡) は,基点座標系に対する回転行列で表すことができる．

3

cos 𝜙𝑹𝑆 = [ 0− sin 𝜙F0 sin 𝜙 110 ] [00 cos 𝜙 0cos 𝜓0− sin 𝜃] [ sin 𝜓cos 𝜃00cos 𝜃sin 𝜃− sin 𝜓cos 𝜓000]
1(2)cos 𝜙 sin 𝜓 + sin 𝜙 sin 𝜃 sin 𝜓cos 𝜃 sin 𝜓=[− sin 𝜙 cos 𝜓 + cos 𝜙 sin 𝜃 sin 𝜓− cos 𝜙 sin 𝜓 + sin 𝜙 sin 𝜃 cos 𝜓cos 𝜃 cos 𝜓sin 𝜙 sin 𝜓 + cos 𝜙 sin 𝜃 cos 𝜓sin 𝜙 cos 𝜃− sin 𝜃 ]cos 𝜙 cos 𝜃

ここで,θ,φ,ψはそれぞれX軸,Y軸,Z軸回りの角度である．まず，Y 軸回りにφ回転させ，回転後の X’軸回りにθ回転，回転後の Z”軸回りにψ回転させた．足底接地時に,式(2)の行列が初期化され,3 つの角度が 0 に設定される．足底接地時からのセンサの向き（3 つの回転角度）の変化は,式(3)を時間積分することによって推定できる．
𝜙̇cos 𝜃 sin 𝜓[ 𝜃̇ ] = [cos 𝜃 cos 𝜓− sin 𝜃𝜓̇cos 𝜓− sin 𝜓00 −10] 𝝎1(3)しかし,測定された加速度ベクトルは足の並進加速度だけでなく重力加速度も含んでいる．積分を行う場合,重力加速度は測定された加速度ベクトルから差し引かれるべきである．
センサは基本的に足底接地時に静止するので,足底接地している間の並進加速度は無視できるほど小さい．したがって,足底接地時の測定加速度は重力加速度ベクトルと等しいとみなすことができ,足の並進加速度ベクトルは式(4)のようになる．
𝐚̅(𝑡) = 𝐹𝑹𝑆 (𝑡) 𝑆𝐚(𝑡) − 𝐹𝐚(𝑡0 )FF(4)𝐚̅(𝑡)は基点座標系で表された足の並進加速度ベクトルである．しかし,センサが足の甲に取り付けられているので,この座標系は地面座標系に平行でも垂直でもない．したがって,Z軸が重力加速度と平行である地面座標系（𝚺G）にしなければならない．
𝐹𝐳=F𝐚(𝑡0 )| 𝐚(𝑡0 )|(5)𝒛×𝒊=𝒚|𝒛 × 𝒊|(6)𝒚× 𝒛=𝒙(7)𝐹

𝑹𝐺 = [𝒙𝒚𝒛]
4(8)

𝐹

𝑹𝐺 は基底座標系によって表される地面固定座標系の方位行列である．したがって,並進加速度ベクトルは,最終的に新しい座標系で式(9)のように表すことができる．
G𝐚̅(𝑡) = [ 𝐹𝑹𝐺 ]𝑇 𝐹𝐚̅ (𝑡)(9)G𝐚̅(𝑡)は地面固定座標系で表された足の並進加速度ベクトルである．
本研究では,この加速度ベクトルは台形公式を用いて足底接地時（𝑡0 ）から次の足底接地時（𝑡𝑓 ）に時間積分した．積分ドリフトを補正するために,足が再び足底接地した時のセンサの垂直速度は 0 である仮定した．𝑡0 と𝑡𝑓 の間の累積垂直速度から減算される速度オフセ
ットcは式(10)を解くことで推定することができる．
𝑡𝑓

𝑡𝑓

0 = ∫ (𝑎𝑧 (𝑡) − 𝑐) 𝑑𝑡 = ∫ 𝑎𝑧 (𝑡)𝑑𝑡 − 𝑐(𝑡𝑓 − 𝑡0 )
𝑡0

(10)

𝑡0

𝑎𝑧 は重力加速度の補正がされた並進加速度ベクトル G𝐚̅(𝑡)の垂直成分である．そして,歩行する地面の傾きが無視できると仮定して,足が再び地面接地時のセンサの垂直変位は 0 であると仮定した．𝑡0 と𝑡𝑓 の間の累積垂直変位から減算された変位オフセットbは,式(11)を解くことで推定することができる．
𝑡𝑓

𝑡𝑓

0 = ∫ (𝜈𝑧 (𝑡) − 𝑏) 𝑑𝑡 = ∫ 𝜈𝑧 (𝑡)𝑑𝑡 − 𝑏(𝑡𝑓 − 𝑡0 )𝑡0(11)𝑡0

𝜈𝑧 はセンサの垂直速度である．速度に関してはX軸,Y軸ともに地面接地時に 0 であると仮定できるので,式(10)と同様に推定する．
2. 3

実験方法

本研究で提唱した方法の精度を示すために,測定システムの試用実験を行った．まず,歩行動作を対象としない机上での実験を試みた．机上実験として,回転を伴わない動作と回転を伴う動作の 2 種類を測定した．回転を伴わない動作は,X軸,Y軸,Z軸で 600mm の平行移動を5 回ずつ行った．回転を伴う動作では, ,X軸回転させながらX軸とY軸に 600mm の平行移動を 3 回ずつ行い,Y軸回転とZ軸回転させた場合も同様に行った．
歩行動作について,通常歩行のほかに,床とのクリアランスが小さいすり足歩行,歩幅の狭い小刻み歩行,および内外側に足を振る分廻し歩行の 4 種類の歩行を試した．歩行は全て左足から踏み出しの 5 歩を 5 回行った．また,歩幅は通常歩行で約 0.8m,その他の歩行が約 0.4mであった．
通常歩行は第 2.3-1 図,すり足歩行は第 2.3-2 図,分廻し歩行は第 2.3-3 図に示す．
5

第 3 章 実験結果
3.1 机上での実験
机上での実験の回転を伴わない各軸平行移動を第 3.1-1 表に示す．各軸の変位の最大誤差は X 軸 65mm,Y 軸 60mm,Z 軸 35mm,平均誤差は X 軸 12.7mm,Y 軸 22.4mm,Z 軸 7.0mm であった．また,変位の誤差は各軸とも 5%以内であった．
回転を伴う平行移動を第 3.1-2 表, 第 3.1-3 表に示す．表の角度は回転させた軸の角度を表している．また,ドリフトの発生を評価するため,1 例として X 軸回転 X 軸移動の角度,変位を第 3.1-1 図と第 3.1-2 図に示す．X 軸移動の各軸回転の変位の最大誤差は X 軸 114mm,Y軸 203mm,Z 軸 148mm,平均誤差は X 軸 24.3mm,Y 軸 60mm,Z 軸 10mm であった．回転を伴わない X 軸移動に比べて最大誤差と平均誤差ともに値が大きくなり,X 軸移動は回転を伴うことで誤差が大きくなることが確認できた．また,X 軸移動の変位に関しては誤差が 15%程度にまで拡大した．
Y 軸移動の場合の最大誤差は X 軸 44mm,Y 軸 67mm,Z 軸 67mm,平均誤差は X 軸 11.2mm,Y軸 10.4mm,Z 軸 32.2mm となった．回転を伴わない Y 軸移動と比べて最大誤差と平均誤差の変化は小さなものであった．
第 3.1-2 図から Y 軸の変位のドリフト誤差が確認できた．

3.2 歩行動作
水平面上の軌跡と矢状面上の軌跡をそれぞれ第 3.2-1 図 と第 3.2-2 図 に示す．通常歩行の歩幅は約 0.8m であり,推定した軌道は長くなった．その他の歩幅は実際の歩幅とほぼ一致した．
小刻み歩行では通常歩行に比べて水平面の変位が小さくなり,すり足歩行では足と地面との間の垂直変位は他の歩行よりも小さくなった．さらに,分廻し歩行では水平面上の軌跡が曲線になった．
歩行動作の各軸角度を第 3.2-3 図 , 第 3.2-4 図 , 第 3.2-5 図 に示す．足底は,分廻し歩行の場合を除いてつま先が離れた時には底屈し,分廻し歩行とすり足歩行の場合を除いて背屈していた．
Z 軸回りの角度は,膝の回転の動きを示すことがあり,膝が曲がると足が内側に回転し,膝が伸びると外側に回転することを示した．X 軸回りの角度は他の軸の角度に対して小さくなった．
各歩行の水平面上ばらつきを第 3.2-6 図 , 第 3.2-8 図, 第 3.2-10 図 , 第 3.2-12 図に示し,鉛直変位の平均と最大,最小を第 3.2-7 図 , 第 3.2-9 図, 第 3.2-11 図 , 第 3.2-12 図に示した．

6

第 4 章 考察
4.1 机上での実験
机上での実験において,動作に回転が伴うことで変位の誤差が大きくなったのは,重力加速度の影響が完全に除去できなかったことが原因であると考えられる．本システムでは
Arduino でのローパスフィルタが用いられていたが,ハイパスフィルタによる重力加速度の除去[4]の検証が必要であると考えられる．
4.2 歩行動作
第 3.2-2 図から通常歩行と小幅歩行では,つま先が地面から離れる前に軌道が下方向に出ていることが確認できる．この誤差は,足が素早く回転した時に見られたため,僅かな姿勢のずれが原因である可能性がある．これは,本研究で用いているシステムはスポーツのような速い動きには使用できないことを示している．しかし,本システムの測定間隔は 0.017 秒であったが,さらに測定間隔を短くすることによって改善されるかもしれない．考慮すべき他の解決策として,式(10)または式(11)を使用して,変位の場合と同様に角度のオフセット誤差を調整することであった．センサ軸が動いているときに角度オフセットをどのように導入するかが考えられる．
水平面上に通常歩行の軌跡を分析したところ,得られたデータの大半が Y 軸方向で誤差が出ていた．原因を特定することは出来なかったが,重力加速度の影響であると仮定した．
実際の歩幅は X 軸に沿っていたので,Y 軸の変位は重力加速度によって引き起こされた可能性がある．それは,角度の誤差が 5°であったとしても,1 秒間で 0.85m もの変位の誤差が発生するからである．
4 種類の歩行は本システムによって判別することができ,軌道と角度は実際の動きとあらかた一致して推定することができた．本研究では,歩容の特徴が定量的ではなく定性的にどのように測定されたかを評価した．今後の課題は,本システムの精度を評価するためにVICON 等の光学システムを使用することである．また,実際のリハビリテーションの状況で本システムを使用して,対応する動作の評価しようと考えている．

7

第 5 章 結言
本研究では,慣性センサを用いた足の軌跡を測定するシステムの制作にあたり,加速度,角速度の値を用いてプログラムを作成し,試用実験を行った．
得られた実験結果は,通常歩行,小刻み歩行,すり足歩行,分廻し歩行の 4 種類の歩行のそれぞれの特徴を測定することができ,歩容の特徴を判定することができていた．しかし,回転を伴う動作では変位に誤差が生じてしまうことが確認できた．
今後の課題は,変位の誤差の原因であると考えられる重力加速度を除去する手法の導入である．また,本研究では歩行実験の測定精度が評価できていないため,測定精度を検証する必要がある．

8

参考文献
[1] 穂苅真樹，渡辺嘉二郎，栗原陽介，瀬川友輔，鳴尾丈司，“スポーツフォームの運動解
析と計測 : ゴルフドライバースイングフォームの計測”，計測自動制御学会論文集, Vol.
38, No. 11 (2002), pp.922-930.
[2] 古瀬則夫，渡邊高志，星宮望，“圧電式ジャイロスコープを用いた下肢関節角度の簡易計測法”, 生体医工学：日本 エム・イー学会誌, Vol.43, No.4 (2005), pp.538-543.
[3] N. Kitagawa and N. Ogihara, “Estimation of foot trajectory duaring human walking by awearable inertial measurement unit mounted to the foot,” Gait and Posture, Vol.45,pp. 110114,2016.
[4] 藤本裕, 豊島健, 南慶一郎, 橋本通, “3 軸加速度センサーによる姿勢検出アルゴリズムの改良と実用性の検討” 人口臓器 28(1), 78-82, (1999).

9

第 2.1-1 図. 計測システム

第 2.1-2 図. センサ装着図10

第 2.1-3 図. Arduino UNO R3

第 2.2-1 図. 足の軌道算出に用いた座標系

第 2.3-1 図. 通常歩行11

第 2.3-2 図. すり足歩行第 2.3-3 図. 分廻し歩行12

第 3.1-1 表. 回転を伴わない各軸平行移動

AverageAv.error

X軸平行移動620.4 ± 28.4

12.7 [mm]Y軸平行移動604.6 ± 50.1

22.4 [mm]Z軸平行移動586.2 ± 15.6

6.9 [mm]

第 3.1-2 表. 回転を伴う X 軸平行移動

Average
X軸回転
Y軸回転
Z軸回転

Av.error544.7 ± 42.0

24.3 [mm]33.6 ± 0.5

0.3 [deg]

518.5 ± 97.0

60.0 [mm]47.5 ± 3.7

2.1 [deg]475.1 ± 16.6

9.6 [mm]43.0 ± 2.3

1.3 [deg]

第 3.1-3 表. 回転を伴う Y 軸平行移動

Average
X軸回転
Y軸回転
Z軸回転

Av.error622.8± 19.5

11.2 [mm]29.7 ± 3.7

2.1 [deg]

650.0 ± 18.0

10.4 [mm]33.5 ± 0.70.4 [deg]597.1 ± 59.032.2 [mm]47.8 ± 2.11.2 [deg]13403530X軸

25Y軸[deg]

20Z軸15105-51815222936435057647178859299106113
120127134
141148155

0[ms]第 3.1-1 図. Ⅹ軸回転 X 軸移動の角度

[ms]

0-0.1

18152229364350576471
78
859299106113120127134141148155[m]

0.1-0.2-0.3-0.4

X軸-0.5

Y軸-0.6

Z軸-0.7

第 3.1-2 図. Ⅹ軸回転 X 軸移動の変位

14Y軸[m]0.30.20.1X軸[m]

0-0.9

-0.8

-0.7-0.6-0.5-0.4

-0.3

-0.2

-0.1

00.1

-0.1
-0.2
-0.3
-0.4

すり足

小刻み

分廻し

通常歩行

第 3.2-1 図. 歩行動作の水平面上の変位

0.04

鉛直軸 [m]

0.03

0.02

0.01

000.2

0.4

0.6

0.8

分廻し

通常歩行

-0.01

水平面 [m]
-0.02

すり足

小刻み

第 3.2-2 図. 歩行動作の矢状面上の変位

15

1

-10

17
51
85
119
153
187
221
255
289
323
357
391
425
459
493
527
561
595
629
663
697
731
765
799
833
867
901
935
969
1003
1037

[deg]

第 3.2-3 図.

第 3.2-4 図.

160-20-30

歩行動作での Y 軸角度
1037

0986

935

884833782731

通常歩行

680

分廻し

629

小刻み

578

すり足

527

476

425

374

323

272

221

170

20

119

68

17[deg]25

15

105[ms]-5-10

歩行動作での X 軸角度

70

60

5040

すり足
小刻み

30分廻し
通常歩行

2010[ms]

105[ms]0[deg]-5-10-15-20-25

すり足

小刻み

分廻し

通常歩行-30-35

第 3.2-5 図.

歩行動作での Z 軸角度

17

0.2

[m]
-1.2

0-1-0.8-0.6-0.4-0.2

00.2-0.2-0.4

Y軸 -327±85.2[mm]-0.6-0.8[m]X軸 -947±117[mm]-1

第 3.2-6 図.

通常歩行の水平面上ばらつき

0.06

平均 6.4±1.8[mm]
最大 41.6±6.8[mm]
最小 -7.4±3.8[mm]

0.04

[m]

0.02

[m]

000.2

0.4

0.6

0.8

11.2-0.02

第 3.2-7 図.

矢状面上の平均と最大・最小

18

1.41.6

X軸 -361±103[mm][m]0.20.15

Y軸 114±31.6[mm]

0.1

0.05

[m]

0-0.6

-0.5

-0.4

-0.3

-0.2

-0.1

0-0.05

第 3.2-8 図.

すり足歩行の水平面上ばらつき

0.01

平均 -0.05±1.8[mm]
最大 2.9±3.4[mm]0.005

最小 -1.9±2.1[mm]

[m]

[m]

000.1

0.2

0.3

0.4

0.5

0.6

0.7

-0.005-0.01

第 3.2-9 図.

すり足歩行の鉛直変位の平均と最大・最小19

[m]

0.2

X軸 -520±67[mm]Y軸 101±45[mm]

0.15

0.1

0.05

[m]
0-0.7

-0.6

-0.5

-0.4

-0.3

-0.2

-0.1

00.1

-0.05

第 3.2-10 図.

小刻み歩行の水平面上ばらつき

0.05

平均 3.9±0.7[mm]0.04

最大 34±5[mm]
最小 -0.9±0.7[mm]

[m]

0.03

0.02

0.01

[m]

000.1

0.2

0.3

0.4

0.5

0.6

0.7

-0.01

第 3.2-11 図.

小刻み歩行の鉛直変位の平均と最大・最小20

0.2

[m]

X軸 -217±19[mm]Y軸 -232±23[mm]

0.1

[m]
-0.6

0-0.5

-0.4

-0.3

-0.2

-0.1

00.1

-0.1

-0.2

-0.3

第 3.2-12 図.

分廻し歩行の水平面上ばらつき

0.06

平均 4.8±1.1 [mm]0.05

最大 39.5±8.6 [mm]最小 -5.6±1.3 [mm]

0.04

[m]

0.03
0.02
0.01

[m]

000.1

0.2

0.3

0.4

0.5

0.6

-0.01

第 3.2-13 図.

分廻し歩行の鉛直変位の平均と最大・最小21