平成 24 年度情報処理学会関西支部

E-06

支部大会

楽 WALK : スマートフォンを用いた歩行時心拍数推定システム
RAKUWALK : heart rate prediction system during walking with smartphone
隅田 麻由1 水本 旭洋1,2 安本 慶一1Mayu Sumida1 Teruhiro Mizumoto1,2 Keiichi Yasumoto1
概要：
ウォーキングは，健康維持・増進に効果的な運動である．しかし，負担の高い無理なウォーキングは逆効果であり，効果的かつ継続的にウォーキングを実施するには，個人の身体能力に応じた歩行ルートの選択が必要である．歩行中の身体的負担度は心拍数から推定可能であるが，心拍計などの特殊デバイスの装着は手間がかかる．本研究では，消費カロリーや歩行時間の制約を満たし，かつ，ユーザにとって負担度が最小となる歩行ルートを推薦するナビゲーションシステム，楽 WALK の実現を目指し，スマートフォンで利用可能な機能のみを用いた身体的負担度の推定法を提案する．提案手法では，個々人の負担度および負担度の変化を推定するために，機械学習を基に，加速度や歩行速度などの歩行データから心拍数を予測する負担度モデルを作成する．実際のウォーキングで計測したデータに本手法を適用した結果，提案手法により平均絶対誤差 11bpm 以内で心拍数の推定ができることが分かった．

Abstract: Walking is eﬀective exercise for maintaing and promoting healthy life, but high perceived exertion
during walking may be counterproductive. For eﬀective and continuous walking, it is necessary to select a
walking route suitable for individual physical ability. Though perceived exertion during walking can be estimated by heart rate, it is costly for a user to equip with a special device such as a heart rate monitor. In
this paper, aiming to realize a health-conscious walking navigation system that recommends a walking routewith minimal perceived exertion satisfying constraints of calorie consumption and walking hours, we proposea system called RAKUWALK which estimates perceived exertion during walking with only available functionsof a smartphone. For this purpose, we build a perceived exertion model which predicts the heart rate fromwalking data including acceleration and walking speed based on machine learning. We applied the proposed
method to actual walking data and conﬁrmed that the method estimates the heart rate within 11bpm.

1. はじめに度は心拍数より算出可能であるが，心拍計の装着に手間がかかり（課題 1），身体能力に応じて同じ運動量でも負担度

生活習慣病の予防には，運動の継続が重要である．特に，には個人差がある（課題 2）．また，様々な歩行ルートに対

誰でも気軽に行うことができるウォーキングが近年注目さする負担度の推定には大量のデータが必要であり，収集のれており，健康づくりのための運動指針 2006[1] では，１コストが高い（課題 3）．そこで，課題 1 および課題 2 の解日１万歩ウォーキングが推奨されている．しかし，自分の決のため，年齢や性別などの属性により作成したカテゴリ

体力にあったルートを選択しなければ，心臓や関節に負担

ごとに，負担度モデルを作成する．負担度モデルとしては，がかかり，ウォーキングに対する意欲の低下や継続性の低

スマートフォンから取得する加速度の振幅や周波数，歩行下につながる．特に高齢者や特定の持病保持者については，

速度や道の勾配などのデータと，心拍計により計測した実負担が高い状態が続くことが危険な場合もある．

際の心拍数データに機械学習を適用することで，心拍数予関連研究として，心拍計および加速度センサによる運動強測モデルを構築する．加えて，課題 3 の解決のため，ユー

度推定 [2][3] や心拍数予測 [4] を行う手法は存在するが，負ザ参加型センシング [12] によりデータを収集および共有し，担度の推定を目的としたものは存在しない．また既存ナビ

歩行ルート検索時に負担度を推定できるようにする．

ゲーションシステムでは，個人の身体能力を考慮したルー
トを検索することはできない．

提案手法の評価では，負担度モデルに与える入力項目の検討および推定精度の評価を行う．ここでは，実際に一般

本研究では，消費カロリーや歩行時間の制約を満たし，か道で測定したデータから作成した負担度モデルを用意し，つ，ユーザにとって負担度が最小となる歩行ルートを推薦

前回の実験 [5] で入力として有用性が確認できた酸素摂取量するナビゲーションシステム，楽 WALK の実現を目指し，について，一般道での実験における有用性を評価する．更様々な歩行ルートに対する負担度推定手法を提案する．負に，実測値とモデルより得られた予測値を比較することで，担度の推定には，次に示す 3 つの課題が挙げられる．負担

予測誤差および負担度の推定精度を算出するとともに，高12奈良先端科学技術大学院大学
Nara Institute Science and Technology
日本学術振興会特別研究員 DC
JSPS Research Fellow DC

い精度が得られる入力パラメタの組合わせを求める．また，求めた入力パラメタの組合せを用いることで，歩行ルートに関わりなく平均絶対誤差 11bpm 以内での心拍数の推定がきると考えられるが，Feng らは心拍数変動が少ない日常生

可能なことが分かった．

活を対象としているため，ここで考案されているモデルは，

2. 関連研究

本研究に必要な運動中の心拍数予測には使用できない．

2.1 節で述べるように，負担度は心拍数と強い関連性があり，また心拍数は運動強度によって変動することが分かっ

2.4 既存のナビゲーションシステムている．したがって本研究では，歩行中の運動強度と心拍既存のナビゲーションシステムとして，ソニーの徒歩ナビ

数を関連付けることで，歩行中に感じる負担度を推定する．

ゲーション NAV-U[10] や，NAVITIME 社の自転車用ルー

そのためには，歩行ルートに対する個々人の運動強度の推ト検索サービス [11] が提供されている．これらのシステム

定および運動強度の変化の推定が必要である．

では，比較的負担の少ないルートを検索できるが，道の勾配を基準としており，個人の身体能力を考慮したルートを

2.1 身体的負担度に関する知見

検索することはできない．そのため，坂道が少ないルート

個人の身体的負担度の測定には，Borg が定義した RPE（主観的運動強度）[6] がよく用いられ，これは表 1 に示すを選択した場合でも，人によっては負担を感じることが考えられ，身体に無理な負担がかかる可能性がある．

15 段階のボルグ・スケール [14] として知られている．RPEと心拍数との間には直線的な相関関係があり，スケールの2.5 本研究の位置付け

数字は心拍数の 1/10 に対応している [7]．よって，これは既存研究では，運動強度の推定や心拍数の予測が行われ主観的な運動強度を数値化する際に有効であり，心拍数とているが，運動中の負担度推定を目的とした研究は存在し負担度の関連付けが可能となる．

ない．更に 2.2 節で示した Tapia ら [3] の研究では，数分間の運動を対象とした運動強度推定を行っているため，運動

表 1 ボルグ・スケール
RPE

主観的運動強度20
19

主観的運動強度

11

楽である

10
かなりきついきつい

14
13

かなり楽である7システムにおいても，個人の身体能力を考慮したルートを検索できるシステムは存在しない．
これより本研究では，スマートフォン搭載センサのみを

8

16
15

9

強度の変化は考慮されておらず，歩行ルートと運動強度の関連付けも行われていない．また，既存のナビゲーション

12
非常にきつい18
17

RPE

非常に楽である

6
ややきつい用いた負担度推定手法を提案し，個々人の運動強度の推定
および運動強度の変化の推定を行う．そしてこの手法に基
づき，ユーザにとって負担が最小となる歩行ルートを推薦

2.2 運動強度推定に関する研究
Oscar ら [2]，Tapia ら [3] は，加速度センサと心拍計を併用し，「ウォーキング」や「ジョギング」などといった活動
種別の推定 [2][3] およびその運動強度の推定 [3] を行っていするナビゲーションシステムを提案する．

3. 負担度を考慮した歩行ナビゲーション
3.1 目的る．心拍数を特徴量とすることで推定精度が向上している

本研究で提案する負担度を考慮したナビゲーションシス

ことから，運動強度の推定には心拍数の利用が有用である

テム（以下，楽 WALK）は，無理なウォーキングによる継と考えられる．

続性の低下，もしくは，怪我などの危険性を回避すること

最近では SUUNTO 社の t6d[8] や POLAR 社の RS100[9]を目的とする．そのため提案システムでは，歩行ルート検のような腕時計型のワイヤレス心拍計を用いることで，手索時に，負担の少ない推薦歩行ルートとそのルートに対す

軽に運動中の心拍数を測定できる．しかしこのような心拍る負担度を提示する．またウォーキング中には，システム

計は，トランスミッターと呼ばれる装置を胸部に直接装着によって負担が高いと判断された場合にのみ警告を表示し，する必要があるため，装着に手間がかかり，煩わしい．加

ユーザに歩行速度の調整を行ってもらう．これにより，ユー

えて，デバイスコストが高いという問題もある．上記で挙ザが推薦歩行ルート以外のルートを歩行する際にも，提案

げた既存研究は，複数の特殊デバイスの装着が必要であるシステムを利用することができる．

ことから，手間やコストの問題が発生すると考えられる．

3.2 要件と基本方針
2.3 心拍数予測に関する研究

楽 WALK には，ユーザに提示する負担度の信頼性が必要

Feng ら [4] は，ニューラルネットワークを用いた心拍数となる．信頼性を向上させるためには，様々な歩行ルート

予測モデルを提案している．ここでは，3 軸加速度センサに対する負担度や歩行データなどの大量なデータが必要でにおける各軸の加速度および 1 ステップ前の心拍数予測値

あり，特殊デバイスを使用しない低コストかつ手軽なデーを入力として与えており，推定平均誤差 5bpm で予測可能タ収集法が求められる．更に，勾配のみを考慮した負担度という結果が報告されている．これよりニューラルネット

ではなく，年齢や性別，運動習慣など個人の身体条件に適し

ワークは運動強度と心拍数の関連付けの手段として利用でた負担度を提示することが望ましい．また，ウォーキング

図 1 楽 WALK の概要

中に警告を表示するためには，負担度の変化を推定する必て運動を行いたいユーザを対象とし，3 軸加速度センサお要がある．歩行中の負担度は，歩行速度や勾配，運動強度

よび GPS センサを搭載したスマートフォン用アプリケーとその継続時間によって変化すると考えられ，例えば，負ションとしての利用を想定する．

担度が高くなるような急な上り坂を歩いたとしても，その図 1 に，楽 WALK の概要図を示す．本システムは，クラ後に平坦な道をしばらく歩けば，負担度は徐々に低くなっ

ウド上に置かれている共有データベースおよびサーバ，そていく．このような連続的な負担度の変化まで推定するた

してユーザが装着するスマートフォンから構成される．

めには，瞬間的な歩行速度や勾配の変化に加えて，運動の継続時間も考慮した負担度の推定を行う必要がある．
以上のことから，提案システムを実現するには，下記の共有データベースには，加速度センサおよび GPS センサ
により測定するユーザの歩行データが蓄積される．ここでの歩行データとは，加速度の振幅および周波数，歩行速度，

要件を満たす必要がある．

勾配，位置情報を指す．加えて，共有データベースには，推( 1 ) 普及デバイスによる負担度の測定

定した心拍数および負担度も蓄積する．年齢や性別などの( 2 ) ウォーキング中の負担度変化の推定

属性により分類されるカテゴリごとにデータを蓄積するこ( 3 ) 様々なルートにおける歩行データの効率良い収集とで，様々なユーザに対する歩行データを共有する．

( 4 ) 年齢，性別，身体能力などの個人差の考慮

本システムの歩行ナビゲーションは，スマートフォン内

これらの要件を満たすための基本方針について述べる．

のアプリケーションを介して行う．まずユーザは，希望す要件（1）については，近年普及しているスマートフォンにる運動時間および消費カロリー，プロフィールの入力，加え着目し，搭載されたセンサのみを用いた負担度推定を行うて，マップ上で出発地点および到着地点の選択を行う（図 1

ことで，特殊デバイスの装着に伴う手間やコストなどの負（1））．そして入力されたプロフィールからユーザが属する担を軽減する．要件（2）については，歩行データ（加速度，カテゴリを決定し，共有データベースに蓄積されている該

歩行速度，勾配，位置などの時系列データ，詳細は後述）か

当カテゴリのデータを基に，推薦歩行ルートの決定，更に，ら負担度と関連性の強い心拍数を予測する負担度モデルを

図 2 に示すように，スマートフォンの画面にルートと負担作成する．これは後述するように，ニューラルネットワー

度の表示を行う（図 1（2）〜（5））．推薦歩行ルートの決定

クなどの機械学習を使用し，歩行速度，勾配などのデータ方法については，3.4 節で述べる．

を入力とするモデルを作成する．更に，歩行速度や勾配のユーザのウォーキング開始後，スマートフォンにより取情報を基に算出する酸素摂取量を入力に加え，運動強度の得された歩行データは，定期的にクラウド上にアップロー

変化と運動の継続時間を考慮する．要件（3）については，

ドされる（図 1（2））
．サーバ上では，負担度モデルにこのユーザ参加型センシング（Participatory Sensing）[12] によ

データを与えることで心拍数および負担度の推定を行い，るデータ収集を行う．ここでは，スマートフォン用のデータ

ユーザに無理な負担がないかどうかを確認する（図 1（3））．

収集アプリをユーザに配布し，広範囲における歩行データ

ユーザの歩行速度が適切でない場合には，「もう少しゆっを効率良く収集する．収集データをサーバ上のデータベー

くり歩いてください」や「やや早めに歩いてください」な

スで共有することで，ユーザが初めて歩行するルートに対

どの歩行速度に関するアドバイスを提示する（図 1（3）〜しても負担度の推定が可能となる．またデータが増えるに（5））．ウォーキング中または終了時に，ユーザの主観的なつれ，システムが提示する負担度の信頼性が向上すると考負担度もしくは心拍数データ（心拍計を装着していた場合えられる．要件（4）については，「20 代男性，運動習慣有のみ）をアップロードしても良い．このデータは，ユーザり」，「30 代女性，運動習慣無し」というように，年齢や性の歩行データと共に共有データベースに蓄積される．ユー

別などの属性によりカテゴリを作成する．同一カテゴリに属するユーザは同程度の負担を感じるという前提の下，カ
テゴリ毎に歩行データの収集や負担度モデルの作成を行う
ことで，負担度に関する個人差を考慮する．

3.3 提案システムの設計
以上の節で述べた目的および要件をふまえ，楽 WALK の概要について述べる．本システムは，ウォーキングを通し

図 2 歩行ルート推薦画面の例

表 2 心拍数と負担度評価の関係

ザビリティを保つという観点から，主観的負担度の評価も
しくは心拍数データの提供は必須ではなく，各ユーザの意志に委ねる．

3.4 推薦歩行ルートの決定方法
以下に示す推薦歩行ルート決定アルゴリズムを用いて，

レベル

主観的評価

心拍数 (bpm)

1

負担を感じない

60〜99

2

やや負担を感じる

100〜119

3

負担は感じるが，運動の継続は可能

120〜1394

歩行速度（運動強度）の低下が必要

140〜1695

休憩が必要

170〜200

ユーザが希望する運動時間および消費カロリーを満たし，かつ，ユーザにとって負担が最小の歩行ルートを推薦する．
手順 1：カテゴリ決定 ユーザが入力したプロフィール
情報により，ユーザのカテゴリを決定する．
手順 2：候補ルートの選択

ユーザがマップ上で選択した出発地点および到着地点の情報を基に，一定距離（2m/

つのルートを選び，表示する．

( 1 ) 手順 4 で求めた消費カロリーが，ユーザの希望する消費カロリー以上である( 2 ) 手順 5 で推定した負担度について，表 2 の 4（歩行速秒×運動時間など）以内の候補となる歩行ルートを全て見度（運動強度）の低下が必要）および 5（休憩が必要）

つけ出し，これらを候補ルート群とする．

の状態になる合計時間が最小である

手順 3：歩行速度の算出 候補ルート群に属する各ルートに対し，歩行速度の算出を行う．ここでは，各ルートの4. スマートフォンによる身体的負担度の推定距離をユーザの希望する運動時間で除算することで歩行速提案手法では，スマートフォンから取得する歩行データ

度を求める．ただしユーザが一定速度で歩行するとは限らから心拍数を予測する負担度モデルを作成する．ここでは，ないため，ルートの勾配を考慮し，勾配に応じて時間を傾

2.1 節で述べたボルグ・スケールに着目し，心拍数から負担斜配分する．勾配は与えられているものとする．

度を定量的に評価できるという仮定の下，歩行データを基

手順 4：消費カロリーの計算 手順 3 で計算した歩行速度および勾配を基に，ACSM（American College of Sports

Medicine）により定義された下記の算出式 [13] を用いて，各ルートに対する消費カロリーを求める．

に心拍数予測を行う．そして，予測した心拍数を表 2 にあてはめることで，負担度を推定する．
加速度の振幅と心拍数の間には，変化の仕方に非線形的な関係が見られる [4]．したがって本研究では，入力である

C = V O2 ∗ 5.01

歩行データと心拍数の非線形な対応関係を再現するニュー

V O2 = R + H + V

ラルネットワークを用い，心拍数の予測を行う．負担度モ

R = 3.5

デルの作成は，（a）スマートフォンによる歩行データの取H = 0.1 ∗ speed

得，ニューラルネットワークに与える入力に関する処理をV = 1.8 ∗ speed ∗ grade

行うための（b）特徴量の抽出，（c）学習データによるモデここで，C[kcal] は消費カロリーを表し，これは酸素摂取量

ルの作成，という 3 つのフェーズに分けられる．以下の節

V O2 [ml/Kg/min] から算出できる．V O2 は，安静時の酸素で（a）〜（c）について述べる．

摂取量 R, 歩行速度に基づく水平方向の酸素摂取量 H ，勾配の増加に伴い必要となる垂直方向の酸素摂取量 V から求められる．また，speed[m/min]，grade[%] はそれぞれ，手順 3 で求めた歩行速度，勾配を表す．
手順 5：心拍数および負担度の推定 各ルートに対し，44.1 （a）スマートフォンによる歩行データの収集
心拍数の予測には，（i）加速度の振幅および周波数，（ii）
歩行速度， （iii）勾配，（iv）酸素摂取量を入力パラメタとして使用する．

章で述べる手法を使用して，歩行時の心拍数および負担度

トレッドミルを用いて歩行中の加速度を測定した実験でを推定する．ここでは，共有データベース内の該当カテゴ

は，歩行速度や勾配によって加速度の振幅，周波数に違い

リにおいて，同一ルートの歩行データを基に心拍数推定をが見られた．また 3.4 節で示したカロリー消費量の算出式行い，負担度とその時間変化について推定する．そして，

より，心拍数の増減と相関が高い酸素摂取量が歩行速度お表 2 のように負担度を設定し，レベル 4（歩行速度（運動よび勾配から算出できることが分かる．これより，心拍数強度）の低下が必要）および 5（休憩が必要）の状態になの変化に関連のあるデータとして，上記（i）〜（iv）の入る場合には，歩行速度を調整した上で，再度，負担度の推力パラメタを考える．

定を行う．なお，表 2 は，表 1 のボルグ・スケール [14] お

本研究では，これらの入力パラメタ（特徴量）を取得す

よび ACSM が定めた運動強度分類表 [13] を基に作成する．

るために，スマートフォンに搭載されている加速度センサ

ASCM より中程度の運動が薦められている [13] ことから，および GPS センサを用いて，歩行中の加速度および位置情

ここでは表 2 のレベル 4 以上を危険レベルとする．

報を収集する．このような時系列のデータがあれば，加速

手順 6：推薦歩行ルートの決定 候補ルート群に属する

度データから加速度の振幅および周波数，歩行速度を，位ルートの中から，下記条件を満たすものを推薦歩行ルート

置情報から勾配を取得することができ，更には，歩行速度とする．推薦ルートが複数ある場合は，ランダムに最大 3および勾配から酸素摂取量が算出できる．

4.2 （b）特徴量の抽出収集した加速度データおよび位置情報から，加速度の振幅および周波数，歩行速度，勾配を算出する．更に，歩行速度および勾配から算出した酸素摂取量を入力として与え，運動強度とその継続時間を考慮した負担度の推定を行う．

4.2.1 加速度の振幅および周波数の抽出
ここでは，3 軸加速度センサにより得られた各軸のデータ
図 3 クロスポイントの検出および 3 軸合成のデータを使用する．予備実験の結果を基に，前処理として周波数 0Hz の直流成分を除去した各デー

る．また AltDif f [m] は，3 秒おきに GPS で取得する 2 点

タを，ハーフ・オーバラップ方式で一定時間のウィンドウ

間の標高差を表す．これは，Google Maps API version3 にに区切り，ウィンドウごとに平均振幅およびパワースペク

よって取得する．

トルにおいてエネルギーが最大となる周波数（以下，周波4.2.4 酸素摂取量の算出

数ピーク）を抽出する．ウィンドウサイズは，使用センサ

酸素摂取量は心拍数と高い相関関係があり，これらは各のサンプリング時間の都合上 600 点の倍数にする必要があ

運動強度に対応した値まで瞬間的に増大するのではなく，るため，ここでは経験的に 1200 点に設定した，運動開始後 2〜3 分で定常状態へと近づく．また歩行状態か

4.2.2 歩行速度の算出

ら安静状態になる場合も同様に，瞬間的に変化するのでは

歩行速度は，加速度データを用いた歩数推定を利用し，なく，徐々に減少する [17]．

歩数推定により得られた歩数およびユーザの歩幅から移動

一 定 強 度 の 運 動 を 行 っ た 場 合 ，酸 素 摂 取 量 の 増 加 量

距離を求めることで算出する．GPS より取得する位置情報

U [ml/kg/min] および減少量 D[ml/kg/min] は，運動に必から移動距離を求める方法も考えられるが，GPS は位置推

要なエネルギーを生み出すために必要な酸素量（以下，酸定における誤差が大きいため，慣性センサ（加速度センサ，

素需要量）K[ml/kg/min]，酸素需要量に収束するまでの時

ジャイロセンサなど）を用いてスタート位置からの相対位間を表すパラメータ T ，時刻 t[s] を用いて次式のように表置を推定するデッド・レコニング [15] を使用する．

される [17]．

U = Ke− t
T
D = K(1 − e− t )T

歩行速度の算出には，次に示す式を用いる．

(4)
(5)

Speed = (Dist/T ime) ∗ 3.6(1)

Dist = Step ∗ Stride length(2)

本研究では，3.4 節の手順 4 で示した ACSM の消費カロ(3)リー算出式を用いて，一定時間 P [s] における平均歩行速度

Stride length = Height − 1

ここで，Speed[km/h] は歩行速度，Dist[m] は移動時間

T ime[s] における移動距離を表す．また，移動距離 Dist[m]は，式（2）に示すように移動時間 T ime[s] における歩数

Step[歩] と歩幅 Stride length[m/歩] から算出する．歩幅はユーザの身長 Height[m] を基に，式（3）より算出する．
歩数推定は，Ying らが提案した推定アルゴリズム [16] を参考にする．このアルゴリズムでは，波形の絶対値が閾値未満になるタイミング（以下，クロスポイント）を検出する
ことで歩数推定を行う．Ying らは閾値を 0 と設定しているが，これでは停止時の検出が難しいため，本研究では，経験的に加速度の分散を 3 倍したものを閾値とする．
図 3 に，加速度データから検出したクロスポイントを赤丸で示す．閾値を通過するポイントを検出しているため，
クロスポイント数の 14 が歩数であると考えることができる．
また，図 3 に示した波形は，Ying らの提案アルゴリズムに従い，Z 軸の加速度データに対し，平滑化フィルタおよび導関数を求めるための微分関数を適応したものである．

4.2.3 勾配の算出
道の勾配は，次に示す式により算出する．

Degree = 100 ∗ AltDif f /Dist
ここで，Degree[％] は勾配，Dist[m] は移動距離を表す．
これは，4.2.2 節で述べた式（1）により算出した値を使用す

および平均勾配から酸素摂取量 Vnow [ml/kg/min] を算出す
る．Vnow が P 秒前の値 Vpre [ml/kg/min] と比較して増加した場合には，K = Vnow として上記の式（4）を用いて増加量を計算する．逆に Vpre と比較して減少した場合には，

K = Vnow として上記の式（5）を用いて減少量を計算する．
Vnow = Vpre の場合には運動強度に変化がないものと考え，
P 秒前の状態が増加傾向であれば K = Vnow として上記の式（4）を，減少傾向であれば K = Vnow として上記の式（5）を用いて変化量を計算する．パラメータ値は，経験的に，T = 53，P = 24 とする．

4.3 （c）学習データによるモデルの作成
負 担 度 モ デ ル の 作 成 に は ，デ ー タ マ イ ニ ン グ ツ ー ル

WEKA3[18] を使用する．このツールでは，入力層，出力層，中間層を持つ階層型ニューラルネットワークが利用できる．ここでの学習は教師あり学習と呼ばれ，学習に必要な教師データ（以下，学習データ）により，予め，入出力の関係を学習させることで負担度モデルを作成する．作成
したモデルに対し，ある時刻 t0 から現在時刻 tn における加速度や歩行速度などの時系列データ（以下，テストデー
タ）を与え，これより算出された時刻 t0 から時刻 tn までの心拍数が出力値（予測値）となる．実際に使用したデータについては，次章で述べる．

5. 負担度モデルによる推定精度の評価
5.1 目的と方法
評価の目的は，4 章で述べた負担度モデルを用いて予測した心拍数および負担度が，提案システムを利用する場合に十分正確であるか調査することである．
これまでに我々は，トレッドミルを用いて測定したデータを用いて，ニューラルネットワークに与える有用な入力についての検討を行い，酸素摂取量の有効性を確認した [5]．
本稿ではこの結果をふまえ，一般道で測定したデータから負担度モデルを作成し，一般道における酸素摂取量の有効性についての評価（評価 1）および負担度の推定精度の評価（評価 2）を行う．
評価 1 では，酸素摂取量を入力に含む場合と含まない場

図 5 測定を行った歩行ルートの標高および距離合の予測結果を比較することで，その有効性について検討する．入力には，各軸の加速度データおよび 3 軸合成の加図 5 にあるように，（A）〜（C）は高低差の大きい道，（D）

速度データから抽出する振幅（AX，AY，AZ，A）および周および（E）は高低差の小さい道であり，歩行時間は（A），

波数ピーク（FX，FY，FZ，F），歩行速度（S），勾配（G），（E）が 20 分，（C），（D）が 40 分，（E）が 7 分である．実

酸素摂取量（VO）の計 11 種類の特徴量を考える．

験では，5 種類のデータのうち１つをテストデータ，残りを評価 2 では，表 2 において心拍数が 20bpm 増加すると学習データとして使用する．

負担度が最大で 2 レベル高くなってしまうことから，平均絶対誤差 10bpm を目指し，負担度の誤検出を 1 レベル以

5.4 評価および考察下に抑えることを目標とする．また，提案システムの実用

5.4.1 酸素摂取量の有効性についての検討の観点から，負担度が低い場合に高いと判定される誤検出

5.1 節で挙げた 11 種類の入力項目のうち，入力数が 6 以（false positive） よりも負担度が高い場合に低いと判定され下となる全組合わせ（1485 通り）について，負担度モデル

る誤検出（false negative）を抑える必要があると考えられを作成し，心拍数予測を行った．テストデータには（A）〜

る．これより本稿では，false positive の出現率（FP（%））（E）をそれぞれ使用し，入力ごとに実測値と予測値の相関および false negative の出現率（FN（%））を算出する．

係数および平均絶対誤差を求めた．
ここで，一般道における酸素摂取量の有用性を評価する

5.2 実験環境ために，酸素摂取量を入力に含む場合（VO あり）と含ま

実験では，腕時計型のワイヤレス心拍計（SUUNTO，t6d）ない場合（VO なし）について，最大相関係数，最小平均を腕に装着し，腰（正面，中央）にスマートフォン（Sony絶対誤差，入力を，表 3，表 4 にそれぞれ示す．このとき，Ericsson，Xperia active）を装着した．スマートフォンは，routeA〜routeE は，それぞれ（A）〜（E）をテストデータ搭載されている加速度センサが図 4 に示す向きになるようとして実験した結果である．

固定した．またサンプリング間隔は，加速度センサを 20 ミ

表 3，表 4 から，相関係数が 0.75–0.90（VO あり），0.55–

リ秒，心拍計を 2 秒，GPS センサを 3 秒とし，被験者を 1

0.79（VO なし），心拍数の平均絶対誤差が 4.75–11.78bpm

人（20 代，女性）とした．

（VO あり），6.83–8.41bpm（VO なし）の推定ができた．これは，目標である誤差 10bpm 以内を概ね満たしている．

5.3 負担度モデル作成に使用するデータ

テストケースごとに「VO あり」および「VO なし」の

一般道におけるデータ測定では，異なる 5 つの歩行ルー

結果を比較すると，相関係数は，routeD を除いて「VO あトを歩いた．図 5 に，各ルートの標高および距離を示す．

り」の方が相関が高い結果となった．一方，平均絶対誤差
は，routeD および routeE では「VO あり」が，残りのテス
トケースについては「VO なし」が，誤差が少ない結果となった．これより，routeA〜routeC では酸素摂取量の有効性が見られたが，routeD〜routeE では有効性が認められなかった.
また，入力に着目すると，「VO あり」の routeA〜routeC
については，歩行速度もしくは勾配を入力に含むことで誤差が小さくなっている．一方，「VO なし」では，振幅およ

図 4 加速度センサの向き

び周波数に加えて，routeA〜routeC は歩行速度および勾配

表 3 最大相関係数および最小平均絶対誤差（VO あり）

表 5 心拍数の予測結果および負担度の推定率（routeA〜routeC の入力：周波数，歩行速度，勾配，酸素摂取量

テスト

相関係数

ケース

最大値

routeArouteB

平均絶対誤差 （bpm）

入力

最小値

0.90

AY, VO

4.75

FZ, S, VO

0.83

F, FX, A, AX, VO

8.78

FZ, A, AX, AY, S, VO

routeD〜routeE の入力：周波数，振幅，歩行速度）

入力

テスト
ケース

相関係数

平均絶対誤差
(bpm)

NE(%)

FP(%)

FN(%)
6.18

routeC

0.77

FZ, VO

8.45

F, FY, D, VO

routeA

0.81

7.23

66.65

27.17

routeD

0.44

FX, FZ, D, S, VO

6.49

F, FY, A, AY, D, VO

routeB

0.69

10.31

45.68

54.32

0.00

routeE

0.75

FX, FZ, AY, AZ, VO

11.78

F, FX, FZ, AX, VO

routeC

0.68

8.83

51.17

30.99

17.84

routeD

0.38

4.55

88.42

5.26

6.32

routeE

0.20

9.91

75.76

0.00

24.24

表 4 最大相関係数および最小平均絶対誤差（VO なし）
テスト

相関係数

平均絶対誤差 （bpm）

最大値

入力

最小値

入力

歩行速度（S），勾配（D），酸素摂取量（VO），高低差の小routeA

0.79

FX, FY, AY, D, S

6.83

F, FX, FY, AX, D, S

さい場合は周波数（FX，FY，FZ），振幅（AZ），歩行速度

routeB

0.75

FX, A, D, S

8.41

FY, FZ, AX, AZ, D, S

routeC

0.60

F, FY, AZ, D

8.93

FX, FY, FZ, D, S（S）となった．

routeD

0.55

FZ, AX, AZ, S

3.02

F, FZ, AZ, S

routeE

0.55

FX, A, AX, S

9.40

FX, FY, FZ, A, AZ, D

ここで，前者の入力を routeA〜routeC に適応した場合，後者の入力を routeD〜routeE に適応した場合の心拍数の予測結果と負担度の推定率を表 5 に示す．負担度の推定率

が，routeD〜routeE は歩行速度が有効に働くことが分かる．

は，表 2 に示す負担度を用いて，負担度を正しく推定した

上記のように，テストケースによって酸素摂取量の有効割合（NE），実際よりも負担度を高く推定した割合（FP），

性に違いが現れたのは，歩行ルートの高低差の違いが原因低く推定した割合（FN）をそれぞれ計算した．表 5 より，であると考えられる．routeA〜routeC のように高低差が大全てのテストケースで心拍数の目標誤差を概ね満たすこと

きい場合は，勾配や歩行速度の変化に伴い心拍数が大きくができた．また，高低差の小さい場合に特に高い正答率が変動し，その特徴を勾配および歩行速度から算出する酸素得られ，全体的に FN の発生率が小さい結果となった．

摂取量によって再現できるため，その有効性が現れた．し更に，高低差が大きい場合の推定例として，routeA の心

かし，routeD〜routeE のように高低差が小さい場合は，勾拍数の予測結果および負担度の分類表を表 6 および図 6 に，配の変化が少ないために酸素摂取量による推定効果が低い．

高低差が小さい場合の推定例として，routeD の心拍数の予また，入力に関しては，routeA〜routeC では歩行速度お測結果および負担度の分類表をそれぞれ表 7 および図 7 によび勾配が，routeD〜routeE では歩行速度が有用であることから，歩行速度や勾配が心拍数変化に大きく影響することが分かる．
正確な推定には歩行速度および勾配が重要であるが，現示す．
分類表は，行方向を実測値の負担度，列方向を予測値の負担度とする．そして表の各欄の値は，対応する負担度（実測レベル，予測レベル）に分類された数の全体に対する割時点では，高低差が小さい所における勾配（ならびに酸素摂合を表わし，表の全ての欄の値を合計すると 100%となる．

取量）の誤差が大きい．また，デッドレコニングによる歩行テストデータのサンプル数は routeA が 81 個，routeD が速度推定には誤差が含まれる．精度を上げるには，勾配お95 個（12 秒おきのサンプル），更に，routeA は負担度のレ

よび歩行速度の正確な計測法，または，入力を歩行ルート

ベル 3 以上，routeD はレベル 5 に該当するデータは含まなによって変える（例えば，勾配を削除）方法が必要である．

い．また，分類表において，NE は対角成分の和，FP は上

5.4.2 評価 2：心拍数および負担度の推定精度

三角成分の和，FN は下三角成分の和として計算する．

評価 1 では，歩行ルートの高低差やテストケースによっ

表 5 および表 6 より，routeA では，負担度の誤推定は最て有用な入力が変わることが分かった．しかし，楽 WALK

大 2 レベル，正答率 67%，負担度 1 レベル以内での正答率

における心拍数予測においては，予測に使用する入力を予は 98.76%という結果が得られた．これより，歩行中に負担め決定しておく必要がある．そこで評価 2 では，高低差のが変化する場合でも，ある程度正しく負担度を推定できる

大きい場合と小さい場合それぞれについて最適な入力を決ことが分かる．また，この場合の FN の発生率は，6.18%で定し，負担度の推定精度を評価する．

あった．

まず，高低差の大きい場合（routeA〜routeC），小さい場

表 5 および表 7 より，routeD では，負担度の誤推定を 1

合（routeD〜routeE）に分け，それぞれに対し最適な入力レベル以内に抑えることができ，推定正答率およそ 88%との決定を行った．前者の場合については，routeA，routeB，いう結果を得ることができた．更に，FN の割合も 6.32%とrouteC のそれぞれに対し，評価 1 で示した 1485 通りの入少なく，歩行ナビゲーションに十分利用できる精度が得ら

力ごとに平均絶対誤差を求めて平均をとり，平均値が最小れたと言える．

となる入力パラメタの組合せを最適な入力とした．後者の場合についても，同様に最適な入力を求めた．その結果，最適な入力は，高低差の大きい場合は，周波数（F，FY），

6. まとめ
本稿では，スマートフォン搭載センサのみを用いた負担

表 6 負担度分類表（高低差が大きい場合）
予

1

測

レ

ベ

ル

2

3

4

5

実

1

7.40

2.47

1.24

0.00

0.00

測

2

1.24

40.74

8.64

0.00

0.00

レ

3

0.00

4.94

17.28

14.82

0.00

ベ

4

0.00

0.00

0.00

1.23

0.00

ル

5

0.00

0.00

0.00

0.00

0.00

[2]

[3]

表 7 負担度分類表（高低差が小さい場合）
予

1

測

レ

2

3

ベ

ル

[4]

4

5

実

1

2.10

4.21

0.00

0.00

0.00

測

2

6.32

86.32

1.05

0.00

0.00

レ

3

0.00

0.00

0.00

0.00

0.00

ベ

4

0.00

0.00

0.00

0.00

0.00

ル

5

0.00

0.00

0.00

0.00

0.00

[5]

[6]

[7]

[8]

図 6 高低差が大きい場合の結果[9]（入力：周波数，歩行速度，勾配，酸素摂取量）

[10]
[11]
[12]

[13]
図 7 高低差が小さい場合の結果（入力：周波数，振幅，歩行速度）度推定手法と，その手法に基づいた，個人の身体能力に応[14]

じた歩行ルートを推薦する楽 WALK について述べた．負担度モデルの評価実験では，歩行ルートの高低差によって有[15]用な入力が変わることが分かり，学習により得られた入力
パラメタの組合せにより，歩行ルートに関わりなく 11bpm
以下の平均絶対誤差で心拍数の推定が可能となった．

[16]

今後は，歩行速度および勾配の算出における誤差を改善すると共に，ルートの高低差に応じた最適な入力についても更に検討を進めたい．また，今回は被験者 1 人のデータ[17]を使用したが，複数の被験者から収集したデータを使用した実験にも取り組みたい．
参考文献
[1]

運動所要量・運動指針の策定検討会：健康づくりの ための運動指針 2006〜生活習慣病予防のために〜, エクササイズガ[18]

イド 2006：厚生労働省, http://www.mhlw.go.jp/bunya/
kenkou/undou01/pdf/data.pdf：（参照 2012-7-10）.
Lara, Ö. D., Përez, A. J., Labrador, M. A. andPosada, J. D.: Centinela: A human activity recognition
system based on acceleration and vital sign data, Journal
on Pervasive and Mobile Computing (2011).
Tapia, E. M., Intille, S. S., Haskell, H. W., Larson, K.,
Wright, J., King, A. and Friedman, R.: Real-Time Recognition of Physical Activities and Their Intensities Using
Wireless Accelerometers and a Heart Rate Monitor, Proc.
of 11th International Symposium on Wearable Computers (ISWC2007), pp. 37–40 (2007).
Xiao, F., Chen, Y., Yuchi, M., Ding, M. and Jo, J.: Heart
Rate Prediction Model Based on Physical Activities Using
Evolutionary Neural Network, Fourth International Conference on Genetic and Evolutionary Computing, pp.
198–201 (2010).
隅田麻由, 今津眞也, 水本旭洋, 安本慶一：健康指向歩行ナビゲーションのためのスマートフォンによる身体的負担度の推定, 情報処理学会研究報告, Vol. 2012-MBL-62,
No. 27, pp. 1–8 (2012).
Borg, G.: Psychophysical scaling with applications inphysical work and the perception of exertion, Scandinavian Journal of Work Environment Health, Vol. 16(SUPPL 1), pp. 55–58 (1990).
Borg, G.: Psychophysical scaling with applications inphysical work and the perception of exertion, Med. Sci.
Sports Exerc, Vol. 5, No. 2, pp. 90–93 (1973).
SUUNTO, : ス ポ ー ツ コ ン ピ ュ ー タ ー TRAINING t6d, http://www.suunto-japan.jp/products/
list/sportscomputer/training/t6d: （参照 2011-0822）.
POLAR, :ハ ー ト レ イ ト モ ニ タ ー RS100,
http://www.polar.co.jp/ja/products/get active/running multisport/RS100: （参照 2011-08-22）.
SONY, : ポ ー タ ブ ル ナ ビ ゲ ー シ ョ ン NAV-U, http://www.sony.jp/nav-u/contents/: （参照 2011-11-29）.
NAVITIME, : 自転車ルート検索, http://www.navitime.
co.jp/bicycle/: （参照 2011-06-12）.
Goldman, J., Shilton, K., Burke, J., Estrin, D.,
Hansen, M., Ramanathan, N., Reddy, S., Samanta, V.
and Srivastiva, M.: Participatory Sensing: A citizenpowered approach to illuminating the patterns that shape
our world, Foresight & Governance Project, White Paper (2009).
Sports Medicine, of A. C.: ACSM’s Guidelines for Exercise Testing and Prescription (7th edition), Lippincott
Williams & Wilkins, Philadelphia (2005).
小野寺孝一, 宮下充正：全身持久性運動における主観的強度と客観的強度の対応性 : Rating of perceived exertion の観点から, 体育學研究, Vol. 21, No. 4, pp. 191–203 (1976).
Montemerli, M., Becker, J., Bhat, S., Dahlkamp, H., Dolgov, D., Ettinger, S. and Haehnel, D.: Junior : thestanford entry in the urban challenge, J Fileld Robot
25(9):569-597 C (2008).
Ying, H., Silex, C., Schnitzer, A. and S. Leonhardt, M. S.:
Automatic Step Detection in the Accelerometer Signal, 4th International Workshop on Wearable and Implantable Body Sensor Networks, Springer Berlin Heidelberg (2007).
佐川貢一, 石原正, 猪岡光, 猪岡英二：歩行形態の違いを考慮した消費カロリーの無拘束推定, 計測自動制御学会, No.
183-7 (1999).
WEKA3, : Data Mining with Open Source Machine
Learning Software in Java, http://www.cs.waikato.ac.
nz/ml/weka: （Accessed 2011-11-29）.