(57)【要約】【課題】  自律航法による歩行ナビゲーションの精度を向上させる。【解決手段】  単位質量に作用する重力と慣性力の合成ベクトルを示すデータを加速度データとして出力する三次元加速度センサから前記加速度データを任意区間内の離散時間に取得し、前記加速度データに基づいて前記三次元加速度センサの加速度の水平成分と鉛直軸成分とを導出し、前記加速度の水平成分の分布の長軸である水平振動軸の前記三次元加速度センサの姿勢を基準とする傾きを統計演算により導出し、前記加速度の前記水平振動軸成分の振動位相と前記加速度の鉛直軸成分の振動位相とのずれに基づいて前記水平振動軸に平行な二方向のうちいずれか一方向を前記任意区間における前記三次元加速度センサの姿勢を基準とする歩行方向として判定する、ことを含む。【選択図】図１
2008116315.tif 000002
【特許請求の範囲】
【請求項１】
単位質量に作用する重力と慣性力の合成ベクトルを示すデータを加速度データとして出力する三次元加速度センサから前記加速度データを任意区間内の離散時間に取得し、
前記加速度データに基づいて前記三次元加速度センサの加速度の水平成分と鉛直軸成分とを導出し、
前記加速度の前記水平成分の分布の長軸である水平振動軸の前記三次元加速度センサの姿勢を基準とする傾きを統計演算により導出し、
前記水平振動軸の傾きに基づいて前記加速度の水平振動軸成分を導出し、
前記加速度の前記水平振動軸成分の振動位相と前記加速度の前記鉛直軸成分の振動位相とのずれに基づいて前記水平振動軸に平行な二方向のうちいずれか一方向を前記任意区間における前記三次元加速度センサの姿勢を基準とする歩行方向として判定する、
ことを含む歩行ナビゲーション方法。
【請求項２】
前記三次元加速度センサに固有の座標系から、静止状態の前記三次元加速度センサが出力する加速度データとみなす静止データの方向がＺ軸と平行になるＸＹＺ直交座標系への写像を導出し、
前記加速度の前記水平成分は、前記ＸＹＺ直交座標系における前記加速度データのＸ、Ｙ成分であって、
前記加速度の前記鉛直軸成分は、前記ＸＹＺ直交座標系における前記加速度データのＺ成分である、
請求項１に記載の歩行ナビゲーション方法。
【請求項３】
複数の前記加速度データの平均を前記静止データとして導出する、
請求項２に記載の歩行ナビゲーション方法。
【請求項４】
前記加速度の前記水平振動軸成分の振動位相と前記加速度の前記鉛直軸成分の振動位相とのずれを、前記加速度の前記水平振動軸成分の振動位相に対する前記加速度の前記鉛直軸成分の振動位相の遅れと前記加速度の前記水平振動軸成分の周期または前記加速度の前記鉛直軸成分の振動周期との比とする、
請求項１から３のいずれか一項に記載の歩行ナビゲーション方法。
【請求項５】
前記三次元加速度センサとともに姿勢変化する磁気センサから磁気データを取得し、
前記磁気データと前記歩行方向とに基づいて前記任意区間における歩行方位を導出する、
請求項１から４のいずれか一項に記載の歩行ナビゲーション方法。
【請求項６】
前記加速度の前記鉛直軸成分に基づいて前記任意区間の歩数を導出し、
基準歩幅と前記任意区間の歩数とに基づいて前記任意区間の歩行距離を導出する、
請求項１から５のいずれか一項に記載の歩行ナビゲーション方法。
【請求項７】
歩数起算地のＧＰＳ位置データと歩数決算地のＧＰＳ位置データとに基づいて前記歩数起算地から前記歩数決算地までの区間距離を導出し、
前記歩数起算値から前記歩数決算値までの歩数である基準歩数を前記加速度の前記鉛直軸成分に基づいて導出し、
前記区間距離と前記基準歩数とに基づいて前記基準歩幅を導出する、
請求項６に記載の歩行ナビゲーション方法。
【請求項８】
前記歩行方位と前記歩行距離とＧＰＳ位置データとに基づいて現在位置を導出する、
請求項６または７に記載の歩行ナビゲーション方法。
【請求項９】
前記任意区間の前記加速度の前記水平成分を（Ｘｐ、Ｙｐ）、ただし（ｐ＝０，１，２...ｓ）とするとき、次式（１）を満たすθ（−９０°＜θ≦９０°）を前記水平振動軸の傾きとして導出する、
請求項１から８のいずれか一項に記載の歩行ナビゲーション方法。
【数１】
2008116315.tif 000018
【請求項１０】
単位質量に作用する重力と慣性力の合成ベクトルを示すデータを加速度データとして出力する三次元加速度センサから前記加速度データを任意区間内の離散時間に取得する手段と、
前記加速度データに基づいて前記三次元加速度センサの加速度の水平成分と鉛直軸成分とを導出する手段と、
前記加速度の前記水平成分の分布の長軸である水平振動軸の前記三次元加速度センサの姿勢を基準とする傾きを統計演算により導出する手段と、
前記加速度の前記水平振動軸成分の振動位相と前記加速度の前記鉛直軸成分の振動位相とのずれに基づいて前記水平振動軸に平行な二方向のうちいずれか一方向を前記任意区間における前記三次元加速度センサの姿勢を基準とする歩行方向として判定する手段と、
を備える歩行ナビゲーションシステム。
【請求項１１】
単位質量に作用する重力と慣性力の合成ベクトルを示すデータを加速度データとして出力する三次元加速度センサから前記加速度データを任意区間内の離散時間に取得する手段と、
前記加速度データに基づいて前記三次元加速度センサの加速度の水平成分と鉛直軸成分とを導出する手段と、
前記加速度の前記水平成分の分布の長軸である水平振動軸の前記三次元加速度センサの姿勢を基準とする傾きを統計演算により導出する手段と、
前記加速度の前記水平振動軸成分の振動位相と前記加速度の前記鉛直軸成分の振動位相とのずれに基づいて前記水平振動軸に平行な二方向のうちいずれか一方向を前記任意区間における前記三次元加速度センサの姿勢を基準とする歩行方向として判定する手段と、
としてコンピュータを機能させる歩行ナビゲーションプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は歩行ナビゲーション方法、システムおよびプログラムに関し、特に自律航法による歩行ナビゲーション方法に関する。
【背景技術】
【０００２】
特許文献１、２、３には、携帯型情報端末において自律航法によって現在地を特定するための技術が開示されている。車両に搭載されるナビゲーションシステムでは移動体である車両の姿勢と磁気センサの姿勢との関係が固定されているのに対し、携帯型情報端末では移動体である歩行者の姿勢に対して磁気センサの姿勢が安定していない。このため携帯型情報端末においては、歩行方位を導出し、歩行方位を用いた自律航法によって現在地を導出することは容易でない。
【０００３】
特許文献１に開示された方法は、三次元加速度センサから出力される加速度データに基づいて、一歩ごとに、三次元加速度センサの加速度の水平成分の微分値が最大となる瞬間における加速度の水平成分の方向を導出し、一歩ごとに、加速度の水平成分が最大となる瞬間における加速度の水平成分の方向を、一歩進む区間の歩行方向として導出する方法である。
【０００４】
【特許文献１】特開２００３−３０２４１９号公報
【特許文献２】特開２００５−２８３３８６号公報
【特許文献３】特開平２−２１６０１１号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかし、特許文献１に開示された方法では、一歩ごとに歩行方向を導出するために誤差が累積しやすいという問題がある。例えば、加速度をサンプリングする時間間隔が広い場合、真の歩行方向と導出される歩行方向とのずれが大きくなる。しかし加速度をサンプリングする時間間隔が狭い場合、演算量が増大するため、サンプリングする時間間隔を闇雲に狭くすることはできない。
【０００６】
しかも特許文献１に開示された方法では、加速度センサが静止した状態において加速度センサの傾きを導出し、そのように導出した加速度センサの傾きを用いることによって、一歩ごとに歩行方向を導出している。しかし、静止状態の加速度センサの姿勢と歩行状態の加速度センサの姿勢が一致する保証がないため、特許文献１に開示された方法では歩行方向の推定誤差が許容範囲を超えて大きくなりうる。
【０００７】
本発明はこのような問題に鑑みて創作されたものであって、自律航法による歩行ナビゲーションの精度を向上させることを目的とする。
【課題を解決するための手段】
【０００８】
（１）上記目的を達成するための歩行ナビゲーション方法は、単位質量に作用する重力と慣性力の合成ベクトルを示すデータとして加速度データを出力する三次元加速度センサから前記加速度データを任意区間内の離散時間に取得し、前記加速度データに基づいて前記三次元加速度センサの加速度の水平成分と鉛直軸成分とを導出し、前記加速度の前記水平成分の分布の長軸である水平振動軸の前記三次元加速度センサの姿勢を基準とする傾きを統計演算により導出し、前記水平振動軸の傾きに基づいて前記加速度の水平振動軸成分を導出し、前記加速度の前記水平振動軸成分の振動位相と前記加速度の前記鉛直軸成分の振動位相とのずれに基づいて前記水平振動軸に平行な二方向のうちいずれか一方向を前記任意区間における前記三次元加速度センサの姿勢を基準とする歩行方向として判定する、ことを含む。
【０００９】
三次元加速度センサが出力する加速度データは三次元ベクトルデータである。三次元加速度センサとともに回転する任意の三軸直交座標系を加速度データについて定めることによって三次元加速度センサの姿勢を基準とする加速度の方向が定まる。任意区間は、時間で区切られても、距離で区切られても、歩数で区切られてもよい任意の区間である。三次元加速度センサの加速度の水平成分とは、三次元加速度センサの加速度の水平成分とみなして処理される二次元ベクトルデータである。三次元加速度センサの加速度の鉛直軸成分とは、三次元加速度センサの加速度の鉛直軸成分とみなして処理される一次元ベクトルデータである。加速度の鉛直軸成分は、加速度の鉛直軸成分と「みなして」処理する一次元ベクトルデータであるから、三次元加速度センサの加速度相当のデータであってもよいし、三次元加速度センサの加速度と重力加速度の合成ベクトル相当のデータであってもよい。加速度データに基づいて三次元加速度センサの加速度の水平成分と鉛直軸成分を導出する方法に制限はなく、例えば特許文献１に記載されている方法を用いてもよい。また本明細書においては、方向と方位を次のように使い分ける。方向というときには、あらかじめ決められた任意の座標系において定義される方向をいい、方位というときには東西南北の四方を基準とする方向をいうものとする。したがって「三次元加速度センサの姿勢を基準とする歩行方向」は、三次元加速度センサとともに座標軸が回転する座標系において定義される方向である。
【００１０】
この歩行ナビゲーション方法では、任意区間において三次元加速度センサから取得する複数の加速度データに基づく統計演算により歩行方向を導出する。このため、母集団を区切る任意区間の長さを適切に設定することにより、一歩毎に歩行方向を導出する方法に比べて精度よく歩行方向を導出することができる。すなわち、特許文献１に開示されているように一歩毎に歩行方向を導出する方法では、一歩ごとに歩行方向を確定するため、歩行方向を確定するたびに確定した誤差が累積する。これに対し、本発明による歩行ナビゲーション方法では、例えば５歩、１０歩と進む任意区間で区切られる加速度データを統計的に処理するため、導出される歩行方向の誤差が小さくなる。
【００１１】
この歩行ナビゲーション方法において統計演算により導出される水平振動軸は本来一意の方向を持たないため、三次元加速度センサの姿勢を基準とする水平振動軸の傾きとして導出されるデータは、歩行方向に平行な方向を一意に特定できる値ではあっても、歩行方向を示している場合もあれば、歩行方向と逆向きを示している場合もある。水平振動軸の傾きとして導出されるデータが歩行方向を示している場合、加速度センサの加速度の水平振動軸成分は加速度センサの加速度の歩行方向軸成分に一致する。尚、ここでいう歩行方向軸とは歩行方向を正にとる向きをもつ座標軸である。後者の場合、加速度センサの加速度の水平振動軸成分は加速度センサの加速度の歩行方向軸成分と大きさが同じで反対向きの関係になる。本発明では歩行者に携行されている加速度センサの加速度の歩行方向軸成分の振動位相と加速度の鉛直軸成分の振動位相がずれる現象に着目し、水平振動軸と平行な２方向のうちいずれか一方向を歩行方向として判定する。もちろん歩行者に携行されている加速度センサの加速度の歩行方向軸成分の振動も鉛直軸成分の振動も完全な周期関数ではないが、歩行方向軸成分の振動においても鉛直軸成分の振動においても一歩ごとに類似した振動パターンが現れ、三次元加速度センサの加速度の歩行方向軸成分の振動位相は鉛直軸成分の振動位相に対して少し進む。したがって例えば、鉛直軸成分の振動位相が水平振動軸成分の振動位相に対して進む場合、水平振動軸の傾きとして導出されたデータが歩行方向と反対の方向を示していると判定し、水平振動軸の傾きを表す角度に１８０°を加えた角度で歩行方向を定めることができる。逆の場合、水平振動軸の傾きとして導出されたデータが歩行方向を示していると判定し、水平振動軸の傾きを表す角度によって歩行方向を定めることができる。
【００１２】
（２）上記目的を達成するための歩行ナビゲーション方法において、前記三次元加速度センサに固有の座標系から、静止状態の前記三次元加速度センサが出力する加速度データとみなす静止データの方向がＺ軸と平行になるＸＹＺ直交座標系への写像を導出してもよい。この場合、前記加速度の前記水平成分は、前記ＸＹＺ直交座標系における前記加速度データのＸ、Ｙ成分であって、前記加速度の前記鉛直軸成分は、前記ＸＹＺ直交座標系における前記加速度データのＺ成分である。
【００１３】
地球の地表において等速直線運動中の質量に作用している力は重力である。加速度を単位質量に作用する重力と慣性力の合成ベクトルを示すデータを加速度データとして出力する三次元加速度センサは、等速直線運動中の加速度を重力加速度と方向が反対の加速度に相当する三次元ベクトルデータとして出力する。すなわち、静止データとは重力加速度と方向が反対の加速度に相当するとみなして処理されるデータである。静止データを導出する方法に特に制限はなく、例えば特許文献１に記載されている方法を用いてもよい。
【００１４】
この歩行ナビゲーション方法では、重力加速度の反対向きの加速度に相当するとみなして処理する静止データの方向がＺ軸と平行になるように加速度データの座標軸を回転させる。すなわち、ＸＹＺ直交座標系において静止データが表現されるとき、加速度データのＸ成分およびＹ成分はともに０になる。このように座標軸を回転させて加速度データを扱うことにより、加速度データの扱いが容易になる。
【００１５】
（３）平地における等速直線運動中に三次元加速度センサが出力する加速度データは重力加速度相当の三次元ベクトルデータである。平地における直線歩行は平地における等速直線運動に近似できる。また平地における曲線歩行は短い時間区間においては平地における直線歩行とみなすことができる。したがって、適切な長さの区間を設定すれば、平地における歩行を平地における等速直線運動に近似することができる。
そこで上記目的を達成するための歩行ナビゲーション方法において、複数の前記加速度データの平均を前記静止データとして導出してもよい。
【００１６】
この方法では、歩行状態で静止データを導出できる。ここで歩行中に取得される加速度データに基づいて導出する静止データの方向と重力加速度の反対方向とのずれと、静止中に取得される加速度データに基づいて導出する静止データの方向と重力加速度の反対方向とのずれについて説明する。静止データの反対方向と重力加速度の方向とのずれは三次元加速度センサの真の姿勢と推定している姿勢とのずれに相当する。静止中に取得される加速度データに基づいて静止データを導出すれば、静止データと重力加速度の反対方向とのずれは、歩行中に取得される加速度データに基づいて静止データを導出する場合に比べて小さくなる。しかし、静止状態での三次元加速度センサの姿勢と歩行中の三次元加速度センサの姿勢とが一致している保証は何もない。静止状態で導出した静止データに基づいて歩行中の加速度センサの姿勢を推定すると、加速度のセンサの静止状態と歩行状態との姿勢差に応じて姿勢の推定誤差が大きくなるため、加速度センサの姿勢の推定誤差が許容できないほど大きくなる可能性がある。その結果、静止状態から推定された加速度センサの姿勢に応じて導出される歩行方向が実際の歩行方向に対する許容できない誤差を含みやすい。これに対し、歩行中に取得される加速度データの平均を静止データとして導出する場合、そのように導出される静止データに対応する三次元加速度センサの姿勢と、三次元加速度センサの真の姿勢との差は、歩行中の加速度センサの運動を等速直線運動に近似することから生ずる誤差でしかないため、許容できないほどの誤差を含みにくい。したがって、本発明によると、水平振動軸の傾きから導出される歩行方向が実際の歩行方向に対する許容できない誤差を含みにくい。
【００１７】
（４）上記目的を達成するための歩行ナビゲーション方法において、前記加速度の前記水平振動軸成分の振動位相と前記加速度の前記鉛直軸成分の振動位相とのずれを、前記加速度の前記水平振動軸成分の振動位相に対する前記加速度の前記鉛直軸成分の振動位相の遅れと前記加速度の前記水平振動軸成分の周期または前記加速度の前記鉛直軸成分の振動周期との比としてもよい。
なお、歩行状態において加速度の水平振動軸成分の振動周期と加速度の鉛直軸成分の振動周期とは一致しているとみなせるため、いずれの振動周期と振動位相の遅れを対比しても実質的な差異はない。
【００１８】
（５）上記目的を達成するための歩行ナビゲーション方法において、前記三次元加速度センサとともに姿勢変化する磁気センサから磁気データを取得し、前記磁気データと前記歩行方向とに基づいて前記任意区間における歩行方位を導出してもよい。
この場合、東西南北の四方を基準とする歩行方位を導出することができる。
【００１９】
（６）前述したとおり歩行状態においては加速度の鉛直軸成分の振動は一歩ごとに類似した振動パターンが現れる周期運動である。したがって、加速度の鉛直軸成分に基づいて歩数を導出可能である。
そこで上記目的を達成するための歩行ナビゲーション方法において、前記加速度の前記鉛直軸成分に基づいて前記任意区間の歩数を導出し、基準歩幅と前記任意区間の歩数とに基づいて前記任意区間の歩行距離を導出してもよい。
【００２０】
この場合、任意区間毎の歩行方向と歩行距離とに基づいて、任意区間の始点を基準とする任意区間の終点の相対位置を導出することが可能になる。すなわち、任意区間毎に始点を基準とする終点の相対位置を導出し、その結果を累積することにより、ＧＰＳ位置データ等で位置が特定されている点を基準として現在地を導出することが可能になる。
尚、基準歩幅は、歩数から歩行距離を推定するために、任意区間の平均歩幅とみなして用いる値であって、演算によって導出してもよいし、ユーザに登録させてもよい。
【００２１】
（７）例えば上記目的を達成するための歩行ナビゲーション方法において、歩数起算地のＧＰＳ位置データと歩数決算地のＧＰＳ位置データとに基づいて前記歩数起算地から前記歩数決算地までの区間距離を導出し、前記歩数起算値から前記歩数決算値までの歩数である基準歩数を前記加速度の前記鉛直軸成分に基づいて導出し、前記区間距離と前記基準歩数とに基づいて前記基準歩幅を導出する、
【００２２】
歩数起算地とは、ＧＰＳ位置データに基づいて歩幅を導出するために暫定的に設定する区間の始点である。歩数決算地とは、ＧＰＳ位置データに基づいて歩幅を導出するために暫定的に設定する区間の終点である。
この場合、基準歩幅や基準歩幅を導出するための身長等の基礎データを登録する手間を省くことができるため、使い勝手が向上する。
【００２３】
（８）上記目的を達成するための歩行ナビゲーション方法において、前記歩行方位と前記歩行距離とＧＰＳ位置データとに基づいて現在位置を導出してもよい。
この場合、衛星航法による現在位置の特定ができない任意区間においても高精度の自律航法によって歩行者の現在位置を導出することができる。
（９）上記目的を達成するための歩行ナビゲーション方法において、前記任意区間の前記加速度の前記水平成分を（Ｘｐ、Ｙｐ）、ただし（ｐ＝０，１，２...ｓ）とするとき、次式（１）を満たすθ（−９０°＜θ≦９０°）を前記水平振動軸の傾きとして導出してもよい。
【数１】
2008116315.tif 000003
【００２４】
ＸＹ座標軸を同一座標平面内でθだけ回転させてｔｕ座標軸に変換する写像は次式（２）によって表される。
【数２】
2008116315.tif 000004
Ｘ軸の回転後の座標軸ｔが水平振動軸となるθを導出するためには、次式（３）で表すｆ（θ）が極大となるθを求めればよい。
【数３】
2008116315.tif 000005
【００２５】
式（１）はｆ（θ）が極値をとる２つのθによって満たされる。式（１）を満たす２つのθのうち一方はｆ（θ）を極大とし、他方はｆ（θ）を極小とする。したがって、式（１）を満たす２つのθのうち、ｆ（θ）に代入するとｆ（θ）がより大きくなる方がｆ（θ）が極大となる値である。このように導出されるθを水平振動軸のＸＹ直交座標系における傾きとする方法は、水平振動軸の傾きを統計演算によって導出する１つの方法である。
【００２６】
（１０）上記目的を達成するための歩行ナビゲーションシステムは、単位質量に作用する重力と慣性力の合成ベクトルを示すデータを加速度データとして出力する三次元加速度センサから前記加速度データを任意区間内の離散時間に取得する手段と、前記加速度データに基づいて前記三次元加速度センサの加速度の水平成分と鉛直軸成分とを導出する手段と、
前記加速度の前記水平成分の分布の長軸である水平振動軸の前記三次元加速度センサの姿勢を基準とする傾きを統計演算により導出する手段と、前記加速度の前記水平振動軸成分の振動位相と前記加速度の前記鉛直軸成分の振動位相とのずれに基づいて前記水平振動軸に平行な二方向のうちいずれか一方向を前記任意区間における前記三次元加速度センサの姿勢を基準とする歩行方向として判定する手段と、を備える。
【００２７】
（１１）上記目的を達成するための歩行ナビゲーションプログラムは、単位質量に作用する重力と慣性力の合成ベクトルを示すデータを加速度データとして出力する三次元加速度センサから前記加速度データを任意区間内の離散時間に取得する手段と、前記加速度データに基づいて前記三次元加速度センサの加速度の水平成分と鉛直軸成分とを導出する手段と、前記加速度の前記水平成分の分布の長軸である水平振動軸の前記三次元加速度センサの姿勢を基準とする傾きを統計演算により導出する手段と、前記加速度の前記水平振動軸成分の振動位相と前記加速度の前記鉛直軸成分の振動位相とのずれに基づいて前記水平振動軸に平行な二方向のうちいずれか一方向を前記任意区間における前記三次元加速度センサの姿勢を基準とする歩行方向として判定する手段と、としてコンピュータを機能させる。
【発明を実施するための最良の形態】
【００２８】
以下、本発明の実施の形態を添付図面を参照しながら次の順序で説明する。
＊＊＊＊＊＊＊＊目次＊＊＊＊＊＊＊＊＊＊＊＊＊
（１）基本原理
（２）歩行ナビゲーションシステムのハードウェア構成
（３）歩行ナビゲーションシステムのデータの流れ
（４）歩行方位の導出方法の詳細
・加速度センサの傾きの推定
・加速度の水平成分および鉛直軸成分の導出
・水平振動軸の傾きの導出
・歩行方向の導出
・歩行方位の導出
（５）効果
（６）他の実施形態
＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊
【００２９】
（１）基本原理
はじめに本発明の実施形態の理解を容易にする基本原理を図１、図２に基づいて説明する。
図１のｘ、ｙ、ｚ軸は加速度センサと磁気センサに固有の座標軸を示している。図１においては、加速度センサと磁気センサとが相対的に固定されており、加速度センサに固有の各座標軸の方向と磁気センサに固有の各座標軸の方向とが一致している（対応する座標軸の方向が一致し、かつ対応する座標軸にとる変数が一致している。）ものと仮定している。
【００３０】
直線歩行中の歩行者に一定の姿勢を保たれて携行されている加速度センサの離散時間の出力をある時間区間についてｘｙｚ座標系にプロットすると、出力は立体的な有限分布を示す。歩行者に携行されている加速度センサは加速度センサが位置する身体部位の運動にともなって三次元運動するが、水平面上の直線歩行中に起きるその運動は二次元的である。すなわち、加速度センサの出力である加速度データは、鉛直軸とそれに垂直な歩行方向軸を含む平面近傍に分布する。この平面を歩行振動面というとき、歩行振動面は水平面に対して垂直な面である。加速度データを、一次元ベクトルである鉛直軸成分と二次元ベクトルである水平成分に分解し、さらに水平成分を歩行振動面と平行な成分と歩行振動面に垂直な成分に分解すると、加速度データの水平成分の歩行振動面と平行な成分は歩行方向軸と平行になる。
【００３１】
水平面（水平面は図１に円で示されている。）に垂直に加速度データを投影したデータの分布を考えると、すなわち加速度データの水平成分の分布を考えると、水平面に含まれる歩行方向軸の方向においてその分布の分散は最大になり、歩行方向軸と垂直な方向においてその分散は最小になる。すなわち、水平面に垂直に加速度データを投影したデータに相当する加速度データのＸＹ成分の分散が最大になる方向は歩行方向軸に平行である。
【００３２】
以上のことから、加速度データを空間情報として扱い、加速度センサの傾きとして加速度センサに固有の座標軸に対する鉛直軸の傾きを推定できれば、推定した水平面に加速度データを投影したデータを母集団とする統計演算によって歩行方向軸を特定することができることがわかる。
【００３３】
一定姿勢の運動状態にある加速度センサの傾きを推定する方法として、静止状態にある加速度センサの傾きを運動状態の傾きとして導出する方法もあるが、本実施形態では運動状態にある加速度センサから取得する加速度データに基づいて運動状態にある加速度センサの傾きを推定する方法について説明する。歩行者に携行されている加速度センサの姿勢は時々刻々変化するが、数秒の短い時間内での姿勢変化はそれほど大きくはない。したがって、歩行中の短い時間区間において加速度センサの姿勢は一定であるとみなす。水平面上での直線歩行は等速直線運動とみなすことができる。平らでない地面上での曲線歩行も、１歩以下ではないがあまり多くはない歩数しか進まない時間区間では等速直線運動と見なすことができる。等速直線運動状態にある加速度センサが検出する加速度は重力加速度と反対向きの加速度である。そこで、ある適切な区間を設定し、その区間における加速度データの平均を重力加速度と反対向きの加速度とみなすこととする。このように重力加速度と反対向きの加速度とみなす加速度を静止データというとき、加速度センサに固有の座標系における静止データの方向が鉛直上向きとして推定されることになる。本実施形態ではこのようにして加速度センサの姿勢を基準として鉛直軸上下の方向を推定し、この方向を加速度センサの傾きとして導出する。このようにして加速度センサの傾きを導出するときに歩行方向軸とみなされる軸が水平振動軸である。
【００３４】
ここまで、歩行方向軸とみなす水平振動軸を導出する原理を説明した。歩行方向として推定すべき方向は水平振動軸に平行な２方向のうちいずれか一方向であるが、歩行データの空間的分布を解析しても歩行方向を推定することはできない。そこで、水平振動軸に平行な２方向のうちいずれか一方向を歩行方向として判定するために加速度データが時系列情報として扱われる。
【００３５】
図２は静止状態から歩行状態にかけて加速度センサが出力する加速度データの水平振動軸成分の大きさと鉛直軸成分の大きさとを時系列にプロットしたグラフである。図２では、静止状態において加速度が０を示すように加速度データから静止データを引き算した実験値をプロットしている。このように静止状態において０を示すように加速度データをオフセットすると、オフセットされた加速度データは、加速度センサが検出している加速度を表すデータになる。
【００３６】
図２に示すように、歩行状態においては、加速度データの水平振動軸成分の振動位相は鉛直軸成分の振動位相に比べて少し進む。したがって、加速度データの水平振動軸成分の振動位相が鉛直軸成分に遅れるのではなく進むように加速度データの水平振動軸の正負の向きを定めることにより、水平振動軸に平行な二方向のいずれが歩行方向であるかを判定することができる。
【００３７】
（２）歩行ナビゲーションシステムのハードウェア構成
図３は本発明の一実施形態としての歩行ナビゲーションシステムのハードウェア構成を示すブロック図である。図３には、携帯電話端末、ＰＤＡ、携帯型ナビゲーション装置などの各種の携帯型情報端末として構成される歩行ナビゲーションシステムの要部を示している。
【００３８】
加速度センサ１０は、直交する３軸の成分に分解して加速度を検出する三次元加速度センサである。加速度センサ１０は、ピエゾ抵抗型、静電容量型、熱検知型などどのような検知方式であってもよい。加速度センサ１０は、重力加速度の反対向きの加速度と加速度センサの運動に固有の加速度を合成した加速度を表す加速度データｇ（ｘ、ｙ、ｚ）を出力する。すなわち例えば、加速度センサ１０のｚ軸の正方向が鉛直上向きになる姿勢で自然落下している状態において加速度センサ１０から出力される加速度データは（０、０、０）であり、同じ姿勢で静止した状態において加速度センサ１０から出力される加速度データは（０、０、ｇ）であり（ｇは正の定数）、同じ姿勢で鉛直上向きに加速度の絶対値が重力加速度の絶対値に等しい等加速度直線運動をしている状態において加速度センサ１０から出力される加速度データは（０、０、２ｇ）である。
【００３９】
磁気センサ１１は、ＭＩ素子、ＭＲ素子などを備え、直交する３軸の成分に分解して磁気を検出するセンサである。磁気センサ１１に固有のｘｙｚ座標系の各軸の方向は、加速度センサ１０に固有のｘｙｚ座標系の各軸の方向と一致している。すなわち、磁気センサ１１から出力される磁気データｈ（ｘ、ｙ、ｚ）＝（１、１、１）と加速度センサ１０から出力される加速度データｇ（ｘ、ｙ、ｚ）＝（１、１、１）とが同じ方向を向くように磁気センサ１１と加速度センサ１０とは例えば共通の基板１７に固定されている。
【００４０】
ＧＰＳ（Global Positioning System）ユニット１２は、衛星航法に用いる３個又は４個の衛星から送られてくる軌道データを受信するアンテナ、現在位置を示す緯度経度データを出力するためのＡＳＩＣ（Application Specific Integrated Circuit）等で構成される。
【００４１】
ＰＵ（Processor Unit）１４は、歩行ナビゲーションプログラム、ＯＳ（Operating System）などを実行することにより、歩行方位を出力したり、歩行ナビゲーションシステムの様々なハードウェアを制御する。
ＲＯＭ（Read Only Memory）１５は、ＰＵ１４が実行する各種のコンピュータプログラムが格納されている不揮発性記憶媒体である。
ＲＡＭ（Random Access Memory）１６は、ＰＵ１４がコンピュータプログラムを実行したり加速度データを格納するためのワークメモリとして機能する記憶媒体である。
ＩＦ（Inter Face）１３は、ＡＤ変換器、ＤＡ変換器などを備え、加速度センサ１０、磁気センサ１１、ＧＰＳユニット１２、ＰＵ１４、ＲＯＭ１５、ＲＡＭ１６などの間で転送されるデータの形式を受取側のデータ形式に変換するモジュールである。
【００４２】
（３）歩行ナビゲーションシステムのデータの流れ
図４は、歩行ナビゲーションシステムのデータの流れの要部を示すデータフロー図である。図４に示す処理は歩行ナビゲーションプログラムを実行するＰＵ１４によって処理される。
【００４３】
傾斜推定プロセス２１は、加速度データｇ（ｘ、ｙ、ｚ）を入力し、加速度センサ１０の傾きの推定値としての（α、β）を加速度データの平均である静止データに基づいて導出する。加速度センサ１０の傾きの推定値は、加速度センサ１０に固有のｘｙｚ座標系において定義する静止データの方向であり、加速度センサ１０に固有のｘ軸と静止データのｘｙ成分とのなす角α（ｙ軸の正方向を正として−１８０°＜α≦１８０°とする。）と、加速度センサ１０に固有のｚ軸と静止データとのなす角β（ｚ軸の正方向を正として−９０°＜β≦９０°とする。）とで表される。
【００４４】
回転プロセス２０は、加速度データｇ（ｘ、ｙ、ｚ）と加速度センサ１０の傾きの推定値（α、β）とを入力し、これらに基づいて加速度データの座標軸を回転することにより、静止状態で加速度センサが検出するとみなす静止データをＺ軸の正方向にとるＸＹＺ直交座標系で表されるデータに加速度データを変換し、その結果として加速度データｇ（Ｘ、Ｙ、Ｚ）を出力する。
【００４５】
水平振動軸導出プロセス２２は、加速度データｇ（Ｘ、Ｙ、Ｚ）のＸＹ成分である加速度の水平成分ｇ（Ｘ、Ｙ）を任意区間について入力し、その分布を解析することにより水平振動軸の傾きを導出する。水平振動軸の傾きはＸＹＺ直交座標系のＸ軸とのなす角θ（−９０°＜θ≦９０°）として導出される。任意区間とはＧＰＳユニット１２が軌道データを受信できない区間である自律航法区間（図５参照）において１つの歩行方位を導出する単位となる区間である。任意区間は、時間で区切ってもよいし、歩数で区切ってもよい。
【００４６】
水平振動軸成分導出プロセス３０は、加速度の水平成分ｇ（Ｘ、Ｙ）と水平振動軸の傾きθとを入力し、これらに基づいて加速度データの座標軸を回転することにより、水平振動軸をｔ軸としｔ軸に垂直で水平面に含まれる軸をｕ軸とするｔｕ直交座標系で表されるデータに加速度データを変換し、その結果として加速度の水平振動軸成分データｇ（ｔ）を導出する。
【００４７】
位相差判定プロセス２３は、加速度の水平振動軸成分であるｇ（ｔ）と加速度データｇ（Ｘ、Ｙ、Ｚ）のＺ成分である加速度の鉛直軸成分ｇ（Ｚ）とを入力し、加速度の水平成分の水平振動軸成分の位相が加速度の鉛直軸成分の位相に対して進んでいるか否かを判定する。
【００４８】
歩行方向導出プロセス２６は、水平振動軸の傾きであるθと位相差判定プロセス２３の判定結果である真／偽とを入力し、歩行方向を導出する。歩行方向はＸＹＺ直交座標系のＸ軸と歩行方向とのなす角であるθ（−９０°＜θ≦９０°）またはθ＋１８０°である。
【００４９】
歩行方位導出プロセス２８は、歩行方向と、磁気データｈ（ｘ、ｙ、ｚ）と、加速度センサ１０の傾き（α、β）とを入力し、これらに基づいて北方向と歩行方向とのなす角である歩行方位を導出する。
【００５０】
以上のプロセスが任意区間の歩行方位を導出するためのプロセスである。次に任意区間の歩行距離を導出するためのプロセスについて説明する。
【００５１】
歩数導出プロセス２４は、歩数起算地から歩数決算地までの歩幅測定区間と任意区間とについて加速度データの鉛直軸成分ｇ（Ｚ）を入力し、これに基づいて任意区間の歩数ｎと、歩数起算地から歩数決算地までの歩数Ｎとを導出する。歩幅測定区間はＧＰＳユニット１２が軌道データを受信できる区間である衛星航法区間（図５参照）において歩行者の平均的な歩幅である基準歩幅が導出される区間である。歩幅測定区間は、時間で区切ってもよいし、歩数で区切ってもよい。任意区間の歩数ｎを導出する方法としては加速度データの大きさが任意区間内で０となる回数を計数し、それに基づいて歩数ｎを導出する方法が知られているが、歩行者に携行されている加速度センサ１０の加速度の歩行方向軸成分と鉛直軸成分の大きさがそれぞれ０になる時点にずれがあるため、成分の大きさを０とみなす適切な閾値を設定することが難しいという問題がある。このため、加速度データの鉛直軸成分のみを使って歩数ｎを導出することが望ましい。
【００５２】
歩幅導出プロセス２５は、歩幅測定区間の歩数Ｎと、歩数起算地と歩数決算地とのＧＰＳ位置データｐ（ｘ、ｙ）とを入力し、これらに基づいて基準歩幅ｗを導出する。
【００５３】
歩行距離導出プロセス２７は、任意区間の歩数ｎと基準歩幅ｗとを入力し、これらを掛け合わせることによって歩行距離ｄを導出する。
【００５４】
以上が任意区間の歩行距離を導出するプロセスである。次に歩行者の現在位置を導出する現在位置導出プロセス２９について説明する。
【００５５】
現在位置導出プロセス２９は、任意区間の歩行方位と任意区間の歩行距離と衛星航法区間のＧＰＳ位置データとを入力し、これらに基づいて現在位置Ｐ（ｘ、ｙ）を導出する。現在位置Ｐ（ｘ、ｙ）は現在地の緯度経度を示すベクトルデータである。現在位置導出プロセス２９は図５に示す衛星航法区間においてはＧＰＳ位置データをそのまま現在位置として出力し、自律航法区間においては次のように現在位置を導出する。すなわち、ＧＰＳユニット１２が軌道データを受信できなくなると、最後に取得したＧＰＳ位置データを最初の任意区間の始点として設定する。その後、任意区間（ｐ３／ｐ４、ｐ４／ｐ５、ｐ５／ｐ６、ｐ６／ｐ７、...）の歩行方位と歩行距離とに基づいて導出する任意区間の歩行ベクトルを、最初の任意区間の始点に積算することによって自律航法区間の現在位置Ｐ（ｘ、ｙ）を導出する。
【００５６】
（４）歩行方位の導出方法の詳細
次に上述のナビゲーションシステムによる歩行方位の導出方法の詳細を説明する。
・加速度センサの傾きの推定
図６は歩行中の任意区間の加速度データを加速度センサ１０に固有のｘｙｚ座標系にプロットした図である。図６に示すように、歩行中の加速度データは原点から中心部が離れた有限範囲に分布する。加速度データの分布の中心部が原点から離れる理由は、重力加速度の反対向きの加速度と歩行運動に固有の加速度を合成した加速度を表すデータが加速度センサ１０から出力されるからである。前述したように、本実施形態では加速度データの平均を重力加速度の反対向きの加速度とみなして静止データとし、加速度センサ１０に固有の座標系で定義される静止データの方向を加速度センサの傾きとして導出する。図７は、静止状態から１０秒間歩行した区間についての加速度データの平均と、歩行開始時の静止状態にある加速度センサ１０から出力される加速度データとの比較結果を示す表である。図７に示す比較結果からは、加速度データの平均を静止データとして導出しても実用上十分信頼できる範囲で加速度センサの傾きを推定できることがわかる。例えば、加速度データの平均を静止データとして導出しても±５°程度の誤差範囲で歩行方位を導出することが可能である。
【００５７】
加速度センサ１０の傾きは具体的には例えば次のようにして推定される。
はじめに、加速度センサ１０の傾きを推定するための区間について加速度データｇ（ｘ、ｙ、ｚ）の平均である静止データｇａｖｅ（ｘ、ｙ、ｚ）を導出する。この区間は、加速度センサ１０の傾き推定値を用いて歩行方位を導出する任意区間と一致していてもよいし、一致していなくてもよいし、連続区間であってもよいし、互いに離れた複数の区間であってもよい。例えば、歩行方位を導出する区間の直前の区間において静止データを導出してもよい。また、ある区間において導出した静止データの方向をその区間の後の複数の区間における歩行方位の導出に利用してもよい。
【００５８】
次に、加速度センサ１０の傾きを、加速度センサ１０に固有の座標系であるｘｙｚ座標系において特定する。具体的には、α、βという２つの値によって静止データの方向を定義することとし、加速度センサ１０の傾きを（α、β）として特定する。αは、静止データのｘｙ成分とｘ軸とのなす角であり、静止データのｘｙ成分であるベクトル（ｘ、ｙ）と（１、０）とのなす角の角度（°）として導出する。αは（０、１）から（１、０）に向かう方向を正とする。βは静止データとｚ軸とのなす角であり、静止データと平行でｚ軸を含む座標平面において表した静止データ（（ｘ２＋ｙ２）１／２、ｚ）と、同一座標平面の単位ベクトル（０、１）とのなす角として導出する。
【００５９】
・加速度の水平成分および鉛直軸成分の導出
図８は、静止データの方向がＺ軸の正の方向となるＸＹＺ直交座標系に加速度データをプロットした図である。このＸＹＺ直交座標系ではＺ軸の正方向が鉛直上向きであると推定され、ＸＹ平面が水平面に平行なものとして推定されている。
【００６０】
加速度センサ１０に固有のｘｙｚ直交座標系からＸＹＺ直交座標系への座標軸変換の写像は次式（４）の行列Ｂと加速度データｇ（ｘ、ｙ、ｚ）との積で表される。
【数４】
2008116315.tif 000006
【００６１】
式（４）の右辺の２つの行列の積を加速度センサ１０の傾き（α、β）に基づいて計算することにより、加速度センサ１０に固有のｘｙｚ座標系からＸＹＺ座標系への写像の行列表現である行列Ｂが導出される。尚、ｘｙｚ直交座標系からＸＹＺ直交座標系への座標軸変換の別の方法として四元数を使用する方法もある。
【００６２】
次に加速度センサ１０に固有のｘｙｚ座標系からＸＹＺ座標系への写像を表す行列Ｂを用いて加速度データの座標軸を変換し、静止データの方向がＺ軸の正の方向となるＸＹＺ座標系で表した加速度データｇ（Ｘ、Ｙ、Ｚ）を導出する。
【００６３】
・水平振動軸の導出
図９は、ＸＹＺ直交座標系のＸＹ平面に加速度データを投影した図である。すなわち、図９に実線で示されるデータは加速度データのＸＹ成分であり、加速度の水平成分として見なしている二次元ベクトルデータである。加速度データの水平成分であるＸＹ成分の分布の長軸である水平振動軸を図１０、図１１に示す。図１０はＸＹ平面が水平面であると推定されているＸＹＺ座標系において水平振動軸を示し、図１１は加速度センサ１０に固有のｘｙｚ座標系において水平振動軸を示している。歩行中の加速度の水平成分は、図１０、図１１に示すように長軸の長さと短軸の長さとが明らかに異なる扁平した分布を示す。前述したように水平振動軸に平行な２方向のいずれかを歩行方向として特定するため、加速度の水平成分の分布を解析し、分布の長軸である水平振動軸とＸ軸とのなす角θ（−９０°＜θ≦９０°）を以下の方法で算出することにより、水平振動軸を導出する。
【００６４】
水平振動軸をｔ軸とするｔｕ直交座標系への加速度データの水平成分ｇ（Ｘ、Ｙ）の座標軸変換の写像は次式（５）によって表される。
【数５】
2008116315.tif 000007
【００６５】
ｔ軸が水平振動軸と平行になるとき、式（５）の写像は加速度データの水平成分の分布の主軸変換に相当する。分布の長軸方向において分散が最大になるため、ｔ軸が水平振動軸と平行であるとき、次式（６）で表すｆ（θ）は極大となる。
【数６】
2008116315.tif 000008
【００６６】
式（６）においては歩行方位を導出する対象となる任意区間の加速度ｇ（Ｘ、Ｙ）の配列Ｑを次式（７）で定義している。
【数７】
2008116315.tif 000009
【００６７】
ｆ（θ）が極値を取るとき次式（８）で表すｆ'（θ）は０である。
【数８】
2008116315.tif 000010
【００６８】
したがって、次式（９）を満たす２つのθ１、θ２を計算してそれぞれｆ(θ)に代入し、θ１、θ２のうちｆ(θ)がより大きくなる角度を特定し、特定した角度を水平振動軸とＸ軸とのなすθとして導出する。
【数９】
2008116315.tif 000011
【００６９】
尚、式（９）の計算が発散することを防ぐため、式（９）の右辺の分母の絶対値と分子の絶対値とを予め計算し、絶対値の大きい方を分母にして計算する必要がある。したがって、式（９）の右辺の分母の絶対値と分子の絶対値とを計算し、分子の絶対値の方が小さい場合には式（９）を用いてθを導出し、分子の絶対値の方が大きい場合には式（９）の代わりに次式（１０）を用いてθを導出する。
【数１０】
2008116315.tif 000012
また、水平加速度の分散が最小となる分布の短軸の傾きを導出し、短軸の傾きから長軸の傾きを導出しても良いことはいうまでもない。
【００７０】
ここまでの処理が終了すると水平振動軸の傾きが導出される。ただし、水平振動軸は正負の向きに意味を持たない軸である。ここまでの処理では、任意区間内の時系列データである加速度データを単なる空間情報として扱うことにより正負の向きに意味を持たない水平振動軸の傾きを導出した。
【００７１】
・歩行方向の導出
図１２は加速度データの水平振動軸成分ｔが加速度データの歩行方向軸成分に一致すると仮定して加速度データの水平振動軸成分ｔをプロットしたグラフである。横軸ｐは、水平加速度ｇ（Ｘ、Ｙ）の配列Ｑの各要素の添え字として用いたｐであるから、加速度センサ１０が加速度データを出力した離散時間を表している。図１３は歩行者に携行されている加速度センサ１０の加速度（この加速度は真の加速度である。）の歩行方向軸成分と鉛直軸成分の振動を表すグラフである。加速度データに計測誤差がなく、加速度センサ１０の推定した傾きにも誤差がない場合、真の加速度の歩行方向軸成分は加速度データの水平振動軸成分ｔまたは−ｔと一致し、真の加速度の鉛直軸成分は加速度データの鉛直成分とみなしているＺ成分と一致する。加速度データの計測誤差は公知のオフセット補正（ここでいうオフセット補正は静止データを差し引く補正ではなく、温度特性等を補正する処理である。）により十分小さくなっている。また前述したように加速度センサ１０の傾きの推定精度は十分高い。したがって、加速度データの水平振動軸成分ｔまたは−ｔの振動位相は、加速度データのＺ成分の振動位相に対して少し進む。そこで加速度データの水平振動軸成分ｔと鉛直軸成分Ｚとを時系列データとして扱う例えば次の方法によって、加速度データの水平振動軸成分ｔの振動位相と加速度データの鉛直軸成分Ｚの振動位相とのずれの大きさを特定し、その結果に基づいて水平振動軸に平行な２方向のいずれか一方向を歩行方向として判定する。
【００７２】
まず加速度データの水平振動軸成分ｔが正から負に向かって０になる時点ｐｔ〔ｉ〕（ｉ＝０，１，２...）と加速度の鉛直軸成分Ｙから静止データを差し引いた値の符号がそれぞれ正から負に向かって０になる時点ｐｚ〔ｊ〕（ｊ＝０，１，２...）とを導出し、各ｉについてｐｚ〔ｊ〕とｐｚ〔ｊ＋１〕との間にあるｐｔ〔ｉ〕を導出する。そして次式（１１）のｆ（ｉ，ｊ）が０．５より大きいか小さいかを判定し、小さければθを歩行方向とし、大きければθ＋１８０°を歩行方向とする。
【数１１】
2008116315.tif 000013
【００７３】
式（１１）を用いた上述の判定方法は、加速度センサ１０の加速度について歩行方向軸成分の振動位相と鉛直軸成分の振動位相との差φ（φは歩行方向軸成分の振動位相が進む方向を正とする値である）が常に加速度センサ１０の加速度の歩行方向軸成分の振動周期Ｔの１／２未満となった実験結果に基づいている。尚、歩行者に携行されている加速度センサ１０の加速度については歩行方向軸成分の振動周期と鉛直軸成分の振動周期とがほぼ一致するため、いずれの周期を式（１１）の左辺の分子に対応させても良い。
【００７４】
ｐｔ〔ｉ〕とｐｚ〔ｊ〕の導出は、例えば次のように行う。
加速度データの水平振動軸成分ｔ、鉛直軸成分Ｚを、加速度センサ１０が加速度データを出力した離散時間を表すｐの関数としてそれぞれ表すと、ｔ〔ｐ〕、Ｚ〔ｐ〕となる。加速度センサ１０が加速度データを出力する時間間隔（すなわち加速度データを読み込む時間間隔）をΔｐとする。
はじめに加速度データの水平振動軸成分が正から負に転ずるタイミングを捕捉するため、次式（１２）を満たすｐを導出する。
【数１２】
2008116315.tif 000014
式（１２）を満たすｐが導出されたら、次にｐとｐ＋１との間においてｔ〔ｐ〕を直線近似して値が０となる時点に対応する水平振動軸成分ｐｔ〔ｉ〕を導出する。具体的には次式（１３）を満たすｐｔ〔ｉ〕を導出する。
【数１３】
2008116315.tif 000015
ｐｚ〔ｊ〕もｐｔ〔ｉ〕と同様にして導出することができる。
【００７５】
ただし、歩行していない状態であっても完全に静止していないと、式（１２）によってノイズとなるタイミングが捕捉されることとなる。そこで式（１２）に加えて次式（１４）をも満たす場合にのみ振動位相の差を導出するための有意なタイミングとして捕捉することが望ましい。
【数１４】
2008116315.tif 000016
ただしＣは、例えば０．１Ｇ相当の正の定数とする。
また加速度データの出力間隔が狭い場合には、１個とばし、２個とばしというようにとびとびの加速度データを式（１２）に代入して加速度データの水平振動軸成分が正から負に転ずるタイミングを捕捉してもよい。
【００７６】
水平振動軸に平行な２方向のいずれか一方向を歩行方向として判定する他の方法としては、式（１１）の左辺の分母分子を入れ替えたｆ（ｉ，ｊ）と２との大小比較によって判定しても良いことはもちろんのこと、加速度データの水平振動軸成分ｔと、加速度データの鉛直方向軸成分Ｚから静止データを差し引いた値との差分の積算値の符号によって判定する方法や、加速度データの水平振動軸成分ｔが０となる時点の加速度データのＹ成分の正負によって判定する方法などが考えられる。
【００７７】
θは静止データをＺ軸の正に取るＸＹＺ直交座標系のＸ軸となす角度であるから、上述した方法によってθまたはθ＋１８０°として導出される歩行方向は、図１４、図１５に示すように、加速度センサ１０の一定とみなしている姿勢を基準とする方向である。
【００７８】
・歩行方位の導出
加速度センサ１０の姿勢を基準とする歩行方向から東西南北の四方を基準とする歩行方位を導出する方法は次の通りである。
磁気センサ１１に固有のｘｙｚ座標系の各軸の方向が加速度センサ１０に固有のｘｙｚ座標系の各軸の方向と一致しているという前提条件の下では、加速度データｇの３軸直交座標系の各軸の方向と磁気データｈの３軸直交座標系の各軸の方向とを一致させる座標軸の回転により、磁気データが示す北方向と加速度センサ１０の姿勢を基準とする歩行方向とのなす角度Ｄ（°）が東西南北を基準とする歩行方位となる（図１参照）。
【００７９】
したがって次式（１５）に示す写像を用いて磁気センサ１１に固有のｘｙｚ直交座標系からＸＹＺ直交座標系に磁気データの座標系を変換し、ＸＹＺ直交座標系の磁気データ（ｈＸ、ｈＹ、ｈＺ）が示す方位角から歩行方向θまたはθ＋１８０°を差し引くことにより、東西南北を基準とする歩行方位Ｄを導出できる。
【数１５】
2008116315.tif 000017
【００８０】
（５）効果
以上説明した歩行ナビゲーション方法では、任意の区間において三次元加速度センサから取得する複数の加速度データに基づいた統計演算により歩行方向を導出した。このため、母集団を区切る任意区間の長さを適切に設定することにより、一歩毎に歩行方向を導出する方法に比べて精度よく歩行方向を導出することができる。すなわち、一歩毎に歩行方向を導出する方法では、一歩ごとに歩行方向を確定するため、歩行方向を確定するたびに確定した誤差が累積する。これに対し、上記実施形態の方法では、例えば５歩、１０歩と進む任意区間で区切られる加速度データを統計的に処理することによって、導出される歩行方向、歩行方位、現在位置の誤差がそれぞれ小さくなる。
また、重力加速度の反対向きの加速度に対応する静止データとして加速度データの平均を導出するため、歩行中において加速度センサ１０の姿勢が変化したとしても、その変化によって増大する歩行方向、歩行方位、現在位置の誤差を歩行中に低減することができる。
（６）他の実施形態
以上、実施形態を用いて本発明を説明したが、本発明の技術的範囲が上述の実施形態に限定されないことはいうまでもなく、本発明の要旨を逸脱しない範囲内において上述の実施形態に対して種々変更を加え得ることは勿論である。
【図面の簡単な説明】
【００８１】
【図１】本発明の実施形態に係る模式図である。
【図２】本発明の実施形態に係るグラフである。
【図３】本発明の実施形態に係るブロック図である。
【図４】本発明の実施形態に係るデータフロー図である。
【図５】本発明の実施形態に係る模式図である。
【図６】本発明の実施形態に係る模式図である。
【図７】本発明の実施形態に係る表である。
【図８】本発明の実施形態に係る模式図である。
【図９】本発明の実施形態に係る模式図である。
【図１０】本発明の実施形態に係る模式図である。
【図１１】本発明の実施形態に係る模式図である。
【図１２】本発明の実施形態に係るグラフである。
【図１３】本発明の実施形態に係るグラフである。
【図１４】本発明の実施形態に係る模式図である。
【図１５】本発明の実施形態に係る模式図である。
【符号の説明】
【００８２】
１０：加速度センサ、１１：磁気センサ、１２：ＧＰＳユニット、１３：ＩＦ、１４：ＰＵ、１５：ＲＯＭ、１６：ＲＡＭ、１７：基板、２０：回転プロセス、２１：傾斜推定プロセス、２２：水平振動軸導出プロセス、２３：位相差判定プロセス、２４：歩数導出プロセス、２５：歩幅導出プロセス、２６：歩行方向導出プロセス、２７：歩行距離導出プロセス、２８：歩行方位導出プロセス、２９：現在位置導出プロセス、Ｄ：角度（歩行方位）、ｄ：歩行距離
【図１】
2008116315.tif 000019
【図２】
2008116315.tif 000020
【図３】
2008116315.tif 000021
【図４】
2008116315.tif 000022
【図５】
2008116315.tif 000023
【図６】
2008116315.tif 000024
【図７】
2008116315.tif 000025
【図８】
2008116315.tif 000026
【図９】
2008116315.tif 000027
【図１０】
2008116315.tif 000028
【図１１】
2008116315.tif 000029
【図１２】
2008116315.tif 000030
【図１３】
2008116315.tif 000031
【図１４】
2008116315.tif 000032
【図１５】
2008116315.tif 000033