(57)【要約】【課題】  携帯機器に搭載された３軸電子コンパスを特に自動車に設置して使用する場合においても、オフセット磁界の推定精度、および、方位角の計算精度を向上させること。【解決手段】  測定された地磁気データを、オフセット用フィルタ、方位角用フィルタにそれぞれ導いて補正し、該補正されたオフセット用補正データ、方位角用補正データに基づいて、オフセット磁界の推定と方位角の算出を行う。【選択図】  図１
2006337057.tif 000002
【特許請求の範囲】
【請求項１】
地磁気を検出する地磁気検出手段と、
前記地磁気検出手段により検出された地磁気データを所定の基準に基づいて補正し、該補正により得られた補正データを出力するフィルタ手段と、
前記フィルタ手段から出力された前記補正データに基づき、前記オフセットを算出するオフセット算出手段と、
前記オフセット算出手段により算出された前記オフセットと、前記フィルタ手段から出力された前記補正データとに基づき、前記方位角を算出する方位角算出手段と
を具えたことを特徴とする方位角計測装置。
【請求項２】
前記フィルタ手段は、
前記地磁気検出手段により検出された地磁気データを所定の基準に基づいて補正し、該補正により得られたオフセット用補正データを出力するオフセット用フィルタ手段と、
前記地磁気検出手段により検出された地磁気データを所定の基準に基づいて補正し、該補正により得られた方位角用補正データを出力する方位角用フィルタ手段と
を含むことを特徴とする請求項１記載の方位角計測装置。
【請求項３】
前記フィルタ手段は、
変化量制限手段を備え、前記変化量制限手段は、該変化量制限手段に新たに入力されたデータと、該変化量制限手段が直前に出力したデータとの距離が所定の閾値内に存在するか否かを判断する手段と、
前記距離が前記所定の閾値以下の場合は、前記新たに入力されたデータを次の出力とする手段と、
前記距離が前記所定の閾値を越える場合は、前記直前のデータと前記最新のデータとの中間点で、かつ、該直前のデータからの距離が所定の値である点を次の出力する手段と
を有することを特徴とする請求項１又は２記載の方位角計測装置。
【請求項４】
前記フィルタ手段は、
領域内平均化手段を備え、前記領域内平均化手段は、該領域内平均化手段に新たに入力されたデータが、特定領域内に存在するか否かを判断する手段と、
前記新たに入力されたデータが、前記特定領域外に存在すると判断したとき、該特定領域外となるまでの連続して特定領域内に存在するときのデータの平均値又は中央値を計算し、該計算により得られた値を補正データとして出力すると共に、前記連続して特定領域内に存在するときのデータを破棄し、前記特定領域を前記新たに入力されたデータに基づいて再設定する手段と
を具えたことを特徴とする請求項１ないし３のいずれかに記載の方位角計測装置。
【請求項５】
前記フィルタ手段は、
領域内平均化手段を備え、前記領域内平均化手段は、該領域内平均化手段に新たに入力されたデータが、特定領域内に存在するか否かを判断する手段と、
前記新たに入力されたデータが、前記特定領域内に存在すると判断したとき、該特定領域内のデータの平均値又は中央値を計算し、該計算により得られた値を補正データとして出力する手段と、
前記新たに入力されたデータが前記特定領域内に存在しないと判断したとき、該新たに入力されたデータを補正データとして出力し、前記特定領域内に存在する地磁気データを破棄し、前記特定領域を該新たに入力されたデータに基づいて再設定する手段と
を具えたことを特徴とする請求項１ないし４のいずれかに記載の方位角計測装置。
【請求項６】
前記特定領域は、所定の中心点を中心とする閉じられた領域であり、該中心点は前記フィルタ手段に新たに入力されたデータにより変化することを特徴とする請求項４又は５記載の方位角計測装置。
【請求項７】
前記特定領域は、立方体又は１８面体によって囲まれる領域であることを特徴とする請求項４ないし６のいずれかに記載の方位角計測装置。
【請求項８】
前記フィルタ手段は、
前記フィルタ手段に新たに連続して入力された所定の個数Ｎ（ただし、Ｎは正の整数）のデータからなるデータ群の、所定の平面からの距離のばらつきが、所定の範囲内に存在するか否かを判断する手段と、
前記距離のばらつきが所定の範囲内に存在すると判断したとき、前記所定の個数Ｎのデータからなるデータ群のうち、Ｍ番目（ただし、Ｍは正の整数）に古いデータを次の出力とする手段と
を有する平面データ選択手段
を具えたことを特徴とする請求項１ないし７のいずれかに記載の方位角計測装置。
【請求項９】
前記平面データ選択手段は、
前記所定の個数Ｎのデータからなるデータ群のうち、１番古いデータを次の出力とすることを特徴とする請求項８記載の方位角計測装置。
【請求項１０】
地磁気を検出する工程と、
前記検出された地磁気データを所定の基準に基づいて補正し、該補正により得られた補正データを出力するフィルタ工程と、
前記出力された前記補正データに基づき、前記オフセットを算出するオフセット算出工程と、
前記算出された前記オフセットと、前記出力された前記補正データとに基づき、前記方位角を算出する方位角算出工程と
を具えたことを特徴とする方位角計測方法。
【請求項１１】
前記フィルタ工程は、
前記検出された地磁気データを所定の基準に基づいて補正し、該補正により得られたオフセット用補正データを出力するオフセット用フィルタ工程と、
前記検出された地磁気データを所定の基準に基づいて補正し、該補正により得られた方位角用補正データを出力する方位角用フィルタ工程と
を含むことを特徴とする請求項１０記載の方位角計測方法。
【請求項１２】
前記フィルタ工程は、
フィルタに新たに入力されたデータと、該フィルタが直前に出力したデータとの距離が所定の閾値内に存在するか否かを判断する手段と、
前記距離が前記所定の閾値以下の場合は、前記新たに入力されたデータを次の出力とする手段と、
前記距離が前記所定の閾値を越える場合は、前記直前のデータと前記最新のデータとの中間点で、かつ、該直前のデータからの距離が所定の値である点を次の出力する手段と
を有する変化量制限手段
を具えたことを特徴とする請求項１０又は１１記載の方位角計測方法。
【請求項１３】
前記フィルタ工程は、
フィルタに新たに入力されたデータが、特定領域内に存在するか否かを判断する工程と、
前記新たに入力されたデータが、前記特定領域外に存在すると判断したとき、該特定領域外となるまでの連続して特定領域内に存在するときのデータの平均値又は中央値を計算し、該計算により得られた値を補正データとして出力すると共に、前記連続して特定領域内に存在するときのデータを破棄し、前記特定領域を前記新たに入力されたデータに基づいて再設定する工程と
を有する領域内平均化工程
を具えたことを特徴とする請求項１０ないし１２のいずれかに記載の方位角計測方法。
【請求項１４】
前記フィルタ工程は、
フィルタに新たに入力されたデータが、特定領域内に存在するか否かを判断する工程と、
前記新たに入力されたデータが、前記特定領域内に存在すると判断したとき、該特定領域内のデータの平均値又は中央値を計算し、該計算により得られた値を補正データとして出力する工程と、
前記新たに入力されたデータが前記特定領域内に存在しないと判断したとき、該新たに入力されたデータを出力し、前記特定領域内に存在する地磁気データを破棄し、前記特定領域を該新たに入力されたデータに基づいて再設定する工程と
を有する領域内平均化工程
を具えたことを特徴とする請求項１０ないし１３のいずれかに記載の方位角計測方法。
【請求項１５】
前記特定領域は、所定の中心点を中心とする閉じられた領域であり、該中心点は前記フィルタに新たに入力されたデータにより変化することを特徴とする請求項１３又は１４記載の方位角計測方法。
【請求項１６】
前記特定領域は、立方体又は１８面体によって囲まれる領域であることを特徴とする請求項１３ないし１５のいずれかに記載の方位角計測方法。
【請求項１７】
前記フィルタ工程は、
前記フィルタに新たに連続して入力された所定の個数Ｎ（ただし、Ｎは正の整数）のデータからなるデータ群の、所定の平面からの距離のばらつきが、所定の範囲内に存在するか否かを判断する工程と、
前記距離のばらつきが所定の範囲内に存在すると判断したとき、前記所定の個数Ｎのデータからなるデータ群のうち、Ｍ番目（ただし、Ｍは正の整数）に古いデータを次の出力とする工程と
を有する平面データ選択工程
を具えたことを特徴とする請求項１０ないし１６のいずれかに記載の方位角計測方法。
【請求項１８】
前記平面データ選択工程は、
前記所定の個数Ｎのデータからなるデータ群のうち、１番古いデータを次の出力とすることを特徴とする請求項１７記載の方位角計測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、方位角計測装置、および、方位角計測方法に関する。特に、自動車に設置された状態で地磁気の測定を行い、該測定により取得したデータから算出されるオフセットおよび方位角の計算を行う自動車用の方位角計測装置、および、方位角計測方法に関する。
【背景技術】
【０００２】
近年、地磁気を測定してその方位角を演算する方位角計測装置、いわゆる電子コンパスの開発が盛んになり、小型のものが開発されて、携帯電話やＰＤＡ（Ｐｅｒｓｏｎａｌ  Ｄｉｇｉｔａｌ  Ａｓｓｉｓｔａｎｃｅ）をはじめとする携帯機器に搭載される例がある。携帯機器でも使えるＧＰＳ（Ｇｌｏｂａｌ  Ｐｏｓｉｔｉｏｎｉｎｇ  Ｓｙｓｔｅｍ）装置、およびローカルまたは携帯電話網やワイヤレスＷＡＮ（Ｗｉｄｅ  Ａｒｅａ  Ｎｅｔｗｏｒｋ）網による地図配信・表示アプリケーションと連動することで、端末所持者の現在位置や進行方向を表示して、歩行者ナビゲーションを行う応用例がある。
【０００３】
携帯機器に用いる電子コンパスとしては、Ｘ，Ｙ，Ｚの３軸方向の磁界を検出し、方位角を演算することのできる３軸電子コンパスが有利である。
【０００４】
その第一の理由は、３軸電子コンパスならば、携帯機器の手持ち角度を設定するもしくは傾斜角センサなどから現在の手持ち傾斜角情報を得て、その値を用いて傾斜角補正演算を行うことで、携帯機器を傾けても正しい方位を演算することができ、ユーザにとっての利便性が上がることである。
【０００５】
その第二の理由は、例えば出願人の発明による特許文献１に開示されているような手法を適用することにより、ユーザの自然な動きから得た３軸の磁気センサ出力から携帯機器内部や磁気センサ自身のオフセットを算出し、補正することができるために、ユーザが意識をしてキャリブレーション操作をしなくてもすむからである。
【０００６】
特許文献１に開示されているのは、携帯機器を自由に動かしたときの互いに直交した３軸の磁気センサの出力の軌跡が球の表面に分布することを利用して、その中心座標を、磁気センサの出力点からの距離が等距離になる点を導出することで推定するという方法である。
【０００７】
一方、自動車に設置されたカーナビゲーション装置はよく知られている。これはＧＰＳ装置や、車載の角速度センサ（ジャイロ）、また車軸の回転角センサなどの情報をもとに、自動車の現在位置や進行方向を表示するものである。而して、自動車用のカーナビゲーション装置は得てして高額であり、その普及率は伸びてはいるものの、まだ非搭載車の率も高い。したがって、歩行者ナビゲーションを行う携帯機器を自動車に設置した場合に自動車用のナビゲーションが可能ならば、カーナビゲーション装置をわざわざ購入、設置しなくても簡易的に類似の機能を得ることが可能となり、ユーザにとっては利点となる。
【０００８】
また、携帯機器を自動車内に固定設置した場合のオフセットを求める手法として、特許文献２に開示されたものがある。この特許文献２に開示されているのは、自動車は略水平面を進行するため、地磁気の測定値は球面上ではなくて、球面と平面の交線である空間内の円周上に分布することを用いてオフセットを推定するという方法である。測定された地磁気データ群と自動車内の携帯機器の設置姿勢によって決定される平面との距離を計算し、距離のばらつきが特定の範囲内であればオフセットを計算し、範囲外であれば、ノイズを多く含んだデータとしてオフセットの計算を行わない。携帯機器の設置姿勢によって決定される平面は、携帯機器に搭載されている加速度センサから求めても良いし、測定された地磁気データ群からの距離が和が最小になるように平面を推定し、推定された平面を用いてもよい。
【０００９】
歩行者が携帯機器を所持している際のオフセットは、特許文献１に開示されている方法に加え、特許文献３で出願されている方法で更にオフセットの推定精度を向上させることができる。この特許文献３では、測定された地磁気データ群のなかからよりオフセット推定の為の計算に適したデータを意図的に選択したり、推定されたオフセットの信頼性を算出し、信頼性の如何によりオフセットを破棄することによって、状況にあわせたオフセットの計算を可能にしている。この場合、オフセットが求まっていない場合や、携帯機器の着磁量が事前にわかっている原因によって変化した場合など、精度より時間を優先してオフセットを求めたい場合は、データの選択基準や信頼性の判定基準を甘く設定し、逆に高  精度のオフセットを求めたいときは、基準を厳しくする。
【００１０】
自動車の内部に携帯機器を設定して、３軸電子コンパスを使用する場合でも基本的には上記のような手法を用いることができるが、自動車特有の環境を考慮することによって、自動車内部で使用した場合の精度を向上させることができる。
【００１１】
【特許文献１】国際公開第ＷＯ２００４／００３４７６号パンフレット
【特許文献２】特願２００４−２３５６０８号公報
【特許文献３】ＰＣＴ／２００４／００９３２４
【特許文献４】ＰＣＴ／２００４／０１８８８８
【発明の開示】
【発明が解決しようとする課題】
【００１２】
携帯機器用に設計された３軸電子コンパスによる方位角計測装置を、自動車に設置した状態で使用するためには、下記のような問題点が知られている（例えば、特許文献２参照）。
【００１３】
すなわち、一般に自動車の内部は車体などが磁性体で、着磁をしているために、人間が携帯機器を所持している際の磁界環境と自動車の内部の磁界環境が大きく異なり、そのために方位角センサが感知するオフセット磁界が異なり、方位角の表示の誤差が大きくなる。正確な方位角を表示するためには、自動車内部において、オフセットの補正を行う必要がある。
【００１４】
踏切など大電流が流れる場所や、鉄構造物近くなどの磁気強度が強い場所を通過することによって車体の着磁量即ちオフセット磁界が変化する。
【００１５】
また、携帯機器でナビゲーションを行う場合、フロントパネルの上などの温度変化の激しい場所に設置されることが多い。このため、温度により、携帯機器内部やセンサ自身のオフセットは時々刻々と変化する。
【００１６】
さらに、着磁量が変わらなくとも、自動車の内部で発生する様々な磁気ノイズが原因となったり、走行中、踏切や鉄構造物近くなどの大電流が流れる場所や着磁している構造物の側を通過することが原因となり、磁気センサの出力の軌跡が球面上から大きく離れることがしばしばある。
【００１７】
従って、携帯機器用に設計された３軸電子コンパスによる方位角計測装置を自動車に設置して使用する場合においては、オフセット磁界の推定精度、および、方位角の計算精度を向上させることができず、その結果、正確な方位角を得ることができないという問題がある。
【００１８】
そこで、本発明の目的は、携帯機器に搭載された電子コンパスを特に自動車に設置して使用する場合においても、オフセット磁界の推定精度、および、方位角の計算精度を向上させ、方位角を高精度に算出することが可能な、方位角計測装置、および、方位角計測方法に関する。
【課題を解決するための手段】
【００１９】
本発明は、地磁気を検出する地磁気検出手段と、前記地磁気検出手段により検出された地磁気データを所定の基準に基づいて補正し、該補正により得られた補正データを出力するフィルタ手段と、前記フィルタ手段から出力された前記補正データに基づき、前記オフセットを算出するオフセット算出手段と、前記オフセット算出手段により算出された前記オフセットと、前記フィルタ手段から出力された前記補正データとに基づき、前記方位角を算出する方位角算出手段とを具えることによって、方位角計測装置を構成する。
【００２０】
前記フィルタ手段は、前記地磁気検出手段により検出された地磁気データを所定の基準に基づいて補正し、該補正により得られたオフセット用補正データを出力するオフセット用フィルタ手段と、前記地磁気検出手段により検出された地磁気データを所定の基準に基づいて補正し、該補正により得られた方位角用補正データを出力する方位角用フィルタ手段とを含んでもよい。
【００２１】
前記フィルタ手段は、該フィルタ手段に新たに入力されたデータと、該フィルタ手段が直前に出力したデータとの距離が所定の閾値内に存在するか否かを判断する手段と、前記距離が前記所定の閾値以下の場合は、前記新たに入力されたデータを次の出力とする手段と、前記距離が前記所定の閾値を越える場合は、前記直前のデータと前記最新のデータとの中間点で、かつ、該直前のデータからの距離が所定の値である点を次の出力する手段とを有する変化量制限手段を具えてもよい。
【００２２】
前記フィルタ手段は、該フィルタ手段に新たに入力されたデータが、特定領域内に存在するか否かを判断する手段と、前記新たに入力されたデータが、前記特定領域外に存在すると判断したとき、該特定領域外となるまでの連続して特定領域内に存在するときのデータの平均値又は中央値を計算し、該計算により得られた値を補正データとして出力すると共に、前記連続して特定領域内に存在するときのデータを破棄し、前記特定領域を前記新たに入力されたデータに基づいて再設定する手段とを有する領域内平均化手段を具えてもよい。
【００２３】
前記フィルタ手段は、該フィルタ手段に新たに入力されたデータが、特定領域内に存在するか否かを判断する手段と、前記新たに入力されたデータが、前記特定領域内に存在すると判断したとき、該特定領域内のデータの平均値又は中央値を計算し、該計算により得られた値を補正データとして出力する手段と、前記新たに入力されたデータが前記特定領域内に存在しないと判断したとき、該新たに入力されたデータを出力し、前記特定領域内に存在する地磁気データを破棄し、前記特定領域を該新たに入力されたデータに基づいて再設定する手段とを有する領域内平均化手段を具えてもよい。
【００２４】
前記特定領域は、所定の中心点を中心とする閉じられた領域であり、該中心点は前記フィルタ手段に新たに入力されたデータにより変化させてもよい。
【００２５】
前記特定領域は、立方体又は１８面体によって囲まれる領域としてもよい。
【００２６】
前記フィルタ手段は、前記フィルタ手段に新たに連続して入力された所定の個数Ｎ（ただし、Ｎは正の整数）のデータからなるデータ群の、所定の平面からの距離のばらつきが、所定の範囲内に存在するか否かを判断する手段と、前記距離のばらつきが所定の範囲内に存在すると判断したとき、前記所定の個数Ｎのデータからなるデータ群のうち、Ｍ番目（ただし、Ｍは正の整数）に古いデータを次の出力とする手段とを有してもよい。
【００２７】
前記平面データ選択手段は、前記所定の個数Ｎのデータからなるデータ群のうち、１番古いデータを次の出力としてもよい。
【００２８】
本発明は、地磁気を検出する工程と、前記検出された地磁気データを所定の基準に基づいて補正し、該補正により得られた補正データを出力するフィルタ工程と、前記出力された前記補正データに基づき、前記オフセットを算出するオフセット算出工程と、前記算出された前記オフセットと、前記出力された前記補正データとに基づき、前記方位角を算出する方位角算出工程とを具えることによって、方位角計測方法を提供する。
【発明の効果】
【００２９】
本発明によれば、検出した地磁気データを、オフセット用フィルタ、方位角用フィルタにそれぞれ導いて補正し、該補正されたオフセット用補正データ、方位角用補正データに基づいて、オフセット磁界の推定並びに方位角の算出を行うようにしたので、携帯機器に搭載された３軸電子コンパスを特に自動車に設置して使用する場合においても、オフセット磁界の推定精度、および、方位角の計算精度を向上させることが可能となり、これにより、方位角を高精度に算出することができる。
【発明を実施するための最良の形態】
【００３０】
以下、図面を参照して、本発明の実施の形態を詳細に説明する。
【００３１】
［第１の例］
本発明の第１の実施の形態を、図１〜図５に基づいて説明する。
【００３２】
＜構成＞
図１は、方位角計測装置１００の構成を示す。
【００３３】
本装置１００は、地磁気データを検出する地磁気センサ等からなる３軸地磁気検出部１と、オフセット計算用フィルタ部２と、データ取捨選択部３と、オフセット算出部４と、方位角計算用フィルタ部５と、方位角計算部６と、地磁気オフセット格納部７と、携帯端末設置姿勢格納部８とを備えている。携帯端末設置姿勢格納部８には、傾斜角センサ９により測定された設置姿勢のデータ、あるいは予め決められた設置姿勢が格納される。
【００３４】
＜動作＞
本装置の動作について説明する。
【００３５】
＜歩行時におけるオフセット推定＞
図２は、歩行時におけるオフセット推定のフローチャートを示す。
【００３６】
ステップＳ１では、３軸地磁気検出部１により地磁気を測定し、３軸の地磁気データとして検出する。
【００３７】
ステップＳ２では、その検出された３軸の磁気データは、オフセット計算用データバッファとしてのデータ取捨選択部３に格納する。このデータバッファには、測定された全てのデータが格納されるわけではなく、例えば、新たに測定されたデータと格納されているデータバッファ中のデータの距離が比較され、所定値以上なら新たに測定されたデータを格納する。
【００３８】
ステップＳ３では、そのデータバッファに所定数のデータが蓄積されると、それら蓄積されたデータ群を用いて、オフセット候補が推定される。
ステップＳ４では、その推定されたオフセット候補の信頼性を評価する。オフセット候補は、例えば推定された直近の所定数のオフセット候補群のばらつきで評価される。
オフセット候補群のばらつきが所定数以下であれば、ステップＳ５に進み、所定数を超えていれば、オフセットは更新されない。ばらつきは、データ分布を表現できる指標を用いればよく、例えば、オフセット候補群の各軸の最大値と最小値の差、あるいは、各軸の標準偏差、等を用いることができる。
ステップＳ５では、直近のオフセット候補が新たなオフセットとして採用され、オフセットが更新される。
【００３９】
ステップＳ６では、オフセットが推定されると、測定された磁気データからオフセットを差し引くことによって、地磁気データを得る。地磁気データと、携帯機器の設置姿勢とから、方位角が計算される。設置姿勢は、携帯機器に搭載された傾斜角センサ９で測定された値を用いてもよいし、或いは、水平等の、既定の設置姿勢を用いてもよい。
【００４０】
＜自動車内設置時におけるオフセット推定＞
図３は、自動車内設置時におけるオフセット推定のフローチャートを示す。  
この処理では、前述した図２の歩行時オフセット推定のフローチャートに比べて、オフセット計算用フィルタ部２および方位角計算用フィルタ部５による処理が付加されている。
【００４１】
ステップＳ１０では、図１のオフセット計算用フィルタ部２が、３軸地磁気検出部１により検出された３軸の地磁気データを所定の基準に基づいて補正し、オフセットに関するオフセット用補正データとし、このオフセット用補正データをデータ取捨選択部３に出力する。
【００４２】
ステップＳ１１では、方位角計算用フィルタ部５が、３軸地磁気検出部１により検出された３軸の地磁気データを所定の基準に基づいて補正し、方位角に関する方位角用補正データとし、この方位角用補正データを方位角計算部６に出力する。
【００４３】
ステップ６では、方位角計算部６が、地磁気オフセット格納部７に格納されたオフセットと、携帯端末設置姿勢格納部８に格納された携帯端末設置姿勢のデータと、方位角計算用フィルタ部５から出力された方位角用補正データとに基づき、方位角を算出する。
【００４４】
なお、その他のステップについては、図２と同様と同様に考えられるので、ここでの説明は省略する。
【００４５】
以下、オフセット計算用フィルタによる処理、方位角計算用フィルタによる処理、方位角の計算についてさらに説明する。
【００４６】
（オフセット計算用フィルタによる処理）
自動車内部は、歩行時の磁気環境に比べて、遥かにノイズが多く、また大きい（携帯端末を設置する場所により異なるが、通常カーナビが設置されるようなインスツルメンタルパネル付近では、自動車自身が発生しているノイズ〜１０ｕＴ（参考値）、踏み切りや構造物による外乱ノイズ〜１００ｕＴ（参考値））。特許文献１で開示されている方法では、携帯端末が大きな立体角を動いたときに取得されたデータ群を用いてオフセットを推定した方が精度がよい。
【００４７】
このため、歩行時にオフセットを求める際には、新たに測定されたデータをオフセット推定用のデータバッファに加える前に、新たに測定されたデータとオフセット推定用データバッファ内のデータとの距離を調べ、この距離が所定の基準値以上でない限り新たなデータをデータバッファに加えない。逆に、基準値を大きくすると、データバッファにデータがたまりにくくなり、オフセット推定に要する時間がかかり、携帯機器の着磁量や温度変化に伴うオフセット変化に追随できなくなってしまう。
【００４８】
ところが、自動車の場合、大きなノイズが頻繁にのっているため、自動車が向きを変えていなくても、前述の距離が見かけ上大きくなり、データバッファにノイズを含んだデータばかりが蓄えられ、結果的に誤ったオフセットが推定されてしまう。
【００４９】
上記の理由により、自動車の場合、オフセットを求める際には、測定値をそのまま用いるのではなく、何らかの加工修正（すなわち、本発明に係るオフセット計算用フィルタによる処理）を行い、ノイズを取り除いたデータを用いてオフセット推定するのが好ましい。
【００５０】
（方位角計算用フィルタによる処理）
オフセットが計算されると、測定値からオフセットを差し引くことよって、地磁気を算出することができる。その算出された地磁気から方位角が計算されるが、オフセット計算時と同様の理由により、測定値をそのまま用いると、ノイズにより不必要に方位角算出値が変動してしまう。
【００５１】
そこで、オフセット計算時用のフィルタと同様に、何らかの加工修正（すなわち、本発明に係る方位角計算用フィルタによる処理）フィルタを用いることにより、ノイズを軽減することができる。
【００５２】
（方位角の計算）
方位角の計算の１例について説明する。
【００５３】
図４は、地磁気センサ１、傾斜角センサ（加速度センサ）９を搭載した方位角計測装置１００としての携帯機器を示す。
【００５４】
この携帯機器の測定軸は、携帯機器を水平、北向き（測定Ｘ軸を北向き、測定Ｚ軸を鉛直上向きに一致）に設置時の地磁気測定値（Ｍｘ，Ｍｙ，Ｍｚ）（磁気センサの出力からオフセットを差し引くことによって得られる）は地磁気の大きさをＭ、伏角をＩとすると、（ＭｃｏｓＩ，０，−ＭｓｉｎＩ）となる。
【００５５】
この姿勢から、図５に示すように、測定Ｚ軸を回転軸としてθ、その後、測定Ｙ軸を回転軸としてα回転させた（θとαは回転の向きが逆なので注意）携帯機器で地磁気Ｍｍを測定すると、
【００５６】
【数１】
2006337057.tif 000003
【００５７】
傾斜角センサ９によりαを測定すれば、以下によりθを求めることができる。
慣例により、方位角は、−θによって与えられる。
【００５８】
【数２】
2006337057.tif 000004
【００５９】
（自動車内設置時の地磁気測定値の分布）
携帯機器を水平に設置（測定Ｚ軸を鉛直上向きに一致）したときの測定値Ｍｍは、
【００６０】
【数３】
2006337057.tif 000005
【００６１】
すなわち、平面内の円状に分布する。
携帯機器が、ピッチ角φ、ロール角ηで設置され、自動車が水平面を進行したときの地磁気測定値Ｍｍの分布は、
【００６２】
【数４】
2006337057.tif 000006
【００６３】
となり、その法線を（−ｓｉｎφ，ｃｏｓφｓｉｎη，ｃｏｓφｃｏｓη）とする平面上に分布する。
【００６４】
傾斜角センサ９によって、φ、ηを測定し、平面を求めてもよいし、自動車内の設置姿勢が予め決められているのであれば、その値を用いて平面を決定してもよい。地磁気測定値から平面を決定するには、最小自乗法により測定値を平面に当てはめればよい。
【００６５】
［第２の例］
本発明の第２の実施の形態を、図６〜図１１に基づいて説明する。なお、前述した第１の例と同一部分については、その説明を省略し、同一符号を付す。
【００６６】
本例は、前述した図１のオフセット計算用フィルタ部２の応用例を示す。ここでは、自動車内に方位角計測装置１００としての携帯機器を設置して地磁気を測定する場合、オフセットを求めるのに適したフィルタとして、３つの例を示す。なお、方位角計測装置１００の基本的な構成は、図１と同様であり、ここでの説明は省略する。
【００６７】
以下に示す３つのフィルタの例は、単独で構成してもよいし、任意の順番に並べて構成してもよい。以下の説明では、フィルタの入力は全て地磁気センサ１で測定された値としているが、フィルタを複数個使用する場合は、次のフィルタの入力は、一つ前のフィルタの出力とする。
【００６８】
（オフセット用変化量制限フィルタ）
図６は、オフセット計算用フィルタ部２をオフセット用変化量制限フィルタにより構成した場合におけるフィルタ処理のフローチャートを示す。
【００６９】
2006337057.tif 000007
【００７０】
2006337057.tif 000008
【００７１】
2006337057.tif 000009
【００７２】
ここで、オフセット用変化量制限フィルタについて補足説明する。
【００７３】
歩行時と異なり、車は急激に進路を変更することが難しいので、自動車内に設置した場合の地磁気測定値の変化量は本来小さいはずである。
【００７４】
例えば、自動車内に携帯端末を水平に設置し、地磁気の水平成分を３０ｕＴとする。このとき、自動車は９０度回転するのに３秒かかるとすると、一秒間に１０ｕＴ程度しか変化しない。地磁気の測定周期を１００ｍｓとすれば、連続する測定間隔では１ｕＴ程度しか変化しないことになる。問題となっている自動車内の磁気ノイズは１０ｕＴあるいはそれ以上であり、測定時間間隔が大きくなければ、測定磁気変化量はノイズに比べはるかに小さい。これを利用して、時間的に連続する測定値間で所定値以上測定値が変動した場合には、新しい測定値を下式によって定め、所定値を超えない場合は、測定値をそのまま使用する。すなわち、
【００７５】
【数５】
2006337057.tif 000010
【００７６】
2006337057.tif 000011
【００７７】
上記の方法を用いると、突発的な大きなノイズを排除することができ、また、通常の線形デジタルフィルタのような出力の遅れもない。この場合、自動車が進路を変更すると同時にフィルタ出力も変化する。出力の遅れが無いことは、方位角の演算において重要である。
【００７８】
図７は、地磁気測定値およびフィルタ処理後の地磁気データの変化量の様子を示す。
実際の地磁気測定値は、Ａ１、Ａ２、Ａ３、Ａ４の各点をつなぐ実線で示す軌跡として表される。これに対して、フィルタ後の地磁気データは、Ｂ１、Ｂ２、Ｂ３の各点をつなぐ破線で示す制限された軌跡として表される。
【００７９】
（オフセット用領域内平均フィルタ）
図８は、オフセット計算用フィルタ部２をオフセット用領域内平均フィルタにより構成した場合におけるフィルタ処理のフローチャートを示す。  
【００８０】
測定が開始される前に(一番最初に図８のフローが実行されるより前に)フィルタに入力されたデータ数を表すＮは０に初期化されているとする。
【００８１】
2006337057.tif 000012
【００８２】
2006337057.tif 000013
【００８３】
2006337057.tif 000014
【００８４】
2006337057.tif 000015
【００８５】
2006337057.tif 000016
【００８６】
2006337057.tif 000017
【００８７】
2006337057.tif 000018
【００８８】
2006337057.tif 000019
【００８９】
2006337057.tif 000020
【００９０】
2006337057.tif 000021
【００９１】
また、オフセット用領域内平均フィルタ閾値ｄｐで規定される領域は、立方体又は１８面体によって囲まれる領域としてもよい。
【００９２】
ここで、オフセット用領域内平均フィルタについて補足説明する。
【００９３】
自動車内では、たとえ自動車が直進或いは停止していたとしても、数μＴ以上のランダムノイズが測定される。このようなノイズを排除するため、測定値空間に特定の大きさの領域（以後、特定領域）を設定し、特定領域に入っているデータを平均し、一つの測定値として使用するようにする。自動車が直進中で、ノイズが特定領域に入っている限り、フィルタはデータを出力せず、自動車が曲がり測定値が特定領域を出ると同時に、特定領域に入っていたデータの平均値を計算し出力する。
【００９４】
特定領域を任意の領域に設定すると（例えば球）、測定値が特定領域に入っているか否かを調べ、また、平均値を計算する為には測定したデータを全て蓄えておかねばならないが、この方法は計算量、記憶領域とも負荷が大きく、携帯電話のような小型のシステムには向かない。
【００９５】
計算量、記憶領域の問題を回避できる効果的で単純な方法は、例えば、特定領域を測定値空間上の立方体とすることである（図８のフロー）。すなわち、測定値各軸（ｘ，ｙ，ｚ）の最大値（Ｘｍａｘ，Ｙｍａｘ，Ｚｍａｘ）、最小値（Ｘｍｉｎ，Ｙｍｉｎ，Ｚｍｉｎ）及び和（Ｘｓｕｍ，Ｙｓｕｍ，Ｚｓｕｍ）と、特定領域に存在するデータ数（Ｎ）を格納しておく。新たな測定値と最大、最小値を比べ、結果如何により最大、最小値を更新する。
【００９６】
また、更新された最大値と最小値の差が所定の値以下であれば、和に測定値を加算し、測定値数に１を足す。更新された最大値と最小値の差が所定の値を超えていれば、即ち特定領域外の測定値が取得されたので、それまでの和を測定値数で割った値をフィルタの出力とする。
【００９７】
特定領域を測定値空間上の立方体とするだけでも大抵の場合充分な効果が得られるが、空間内の対称性を良くするために、例えば以下のような１８面体を採用してもよい。
【００９８】
立方体はＸ，Ｙ，Ｚ軸上での測定値の幅を調べたが、さらに、ＸＹ座標系をＺ軸を回転軸として４５度回転させた座標系の座標軸、ＹＺ座標系をＸ軸を回転軸として４５度回転させた座標系の座標軸、ＺＸ座標系をＹ軸を回転軸として４５度回転させた座標系の座標軸、上の測定値の幅を調べる。図１０はこの領域を示したものである。具体的には、
【００９９】
【数６】
2006337057.tif 000022
【０１００】
のそれぞれの最大値、最小値、和、および特定領域に存在するデータ数を格納しておき、それぞれの値が所定の値の範囲を超えるまでは和を蓄えておき、所定の値をどれか一つでも超えた場合、和を測定値数で割り平均値を算出しフィルタの出力とすればよい。
【０１０１】
或いは、測定値Ｘ，Ｙ，Ｚの幅と、√３／３（Ｘ＋Ｙ＋Ｚ），  √３／３（Ｘ＋Ｙ−Ｚ），  √３／３（−Ｘ＋Ｙ＋Ｚ），  √３／３（−Ｘ＋Ｙ−Ｚ）の幅を調べてもよい。図９はこの領域を表したものである。
【０１０２】
フィルタの出力は平均値でなくとも、特定領域を代表する値なら何でもよい。例えば、各軸ごとに測定値を大きさの順に並べ、ちょうど真中の順位の値を出力値としても良い（すなわち、各軸ごとに中央値を出力値とする）。
【０１０３】
特定領域２００は、ＸＹＺ座標系のＸ軸を、まずＺ軸を回転軸として４５度あるいは１３５度回転し、回転後に写ったＹ軸を中心にａｒｃｔａｎ（１／√２）及びａｒｃｔａｎ（−１／√２）回転させた、４つの軸と、Ｘ，Ｙ，Ｚの３軸、合計７軸上での測定値の幅を調べることに相当する。
【０１０４】
図９の例では、特定領域２００は、３軸地磁気測定値を、Ｘ、Ｙ、Ｚとしたとき、
−Ｄ≦Ｘ，Ｙ，Ｚ，√３／３（Ｘ＋Ｙ＋Ｚ），√３／３（Ｘ＋Ｙ−Ｚ），√３／３（−Ｘ＋Ｙ＋Ｚ），√３／３（−Ｘ＋Ｙ−Ｚ）≦Ｄ，Ｄ＞０
を満たす領域である。
【０１０５】
図１０の例では、特定領域２０１は、３軸地磁気測定値を、Ｘ、Ｙ、Ｚとしたとき、
−Ｄ≦√２／２（Ｘ＋Ｙ），√２／２（−Ｘ＋Ｙ），√２／２（Ｚ＋Ｙ），√２／２（−Ｚ＋Ｙ），√２／２（Ｙ＋Ｚ），√２／２（−Ｙ＋Ｚ）≦Ｄ，Ｄ＞０
を満たす領域である。
【０１０６】
実際の応用例では、上式の関係で十分であるが、さらに複雑な領域を設定したければ、測定値の幅を調べる軸を増やせばよい。
【０１０７】
（オフセット用平面データ選択フィルタ）
図１１は、オフセット計算用フィルタ部２をオフセット用平面データ選択フィルタにより構成した場合におけるフィルタ処理のフローチャートを示す。
【０１０８】
ステップＳ４１では、３軸地磁気検出部１によって連続して検出された所定の個数Ｎ（ただし、Ｎは正の整数）の地磁気データからなるデータ群の、所定の平面からの距離のばらつきが、所定の範囲すなわち平面データ選択フィルタ閾値ｄｄ内に存在するか否かを判断する。
【０１０９】
ステップＳ４２では、その所定の平面からの距離のばらつきが、所定の平面データ選択フィルタ閾値ｄｄ内に存在すると判断したとき、所定の個数Ｎの地磁気データからなるデータ群のうち、Ｍ番目（ただし、Ｍは正の整数）に古いデータを次の出力とする。
この場合、所定の個数Ｎの地磁気データからなるデータ群のうち、１番古いデータを次の出力としてもよい。
【０１１０】
ここで、オフセット用平面データ選択フィルタについて補足説明する。
【０１１１】
自動車は、略水平面を進行するので、自動車内に設置した携帯機器で地磁気を測定した場合、測定されるデータ群は、設置姿勢によって決定される平面上に分布する。ただし、この平面は、地磁気の大きさや伏角が変動すると平行移動する（平面の法線方向は変わらない）。地磁気の大きさ及び伏角（厳密に言えば更に偏角）が一定の地域で測定した地磁気は本来、同じ平面上に分布することになり、平面から飛び出したデータはノイズを含んだデータとみなせる。
【０１１２】
上記のフィルタによる処理を利用して、所定の個数（Ｎ）の測定値の、設置姿勢によって決定される平面との距離を計算し、この距離のばらつきが一定値以内であれば、所定の個数（Ｎ）の測定値は信頼できるとする。信頼できるとした所定の個数（Ｎ）の測定値全てをフィルタの出力としてもよいし、Ｎより少ない個数の測定値を出力してもよい。
【０１１３】
例えば、信頼できるＮ個の測定値のうち、一番古いデータ一つを、あるいは、時間軸上の中央値に当たるデータの一つを出力する。設置姿勢によって決定される平面は、携帯機器に搭載されている傾斜角センサ９の測定値から決定してもよいし、地磁気測定値を平面に当てはめ、当てはめられた平面を用いてもよい。また、自動車内の設置姿勢が使用上予め定められているのであれば、定められた設置姿勢から平面を決定してもよい。
【０１１４】
［第３の例］
本発明の第３の実施の形態を、図１２〜図１３に基づいて説明する。なお、前述した各例と同一部分については、その説明を省略し、同一符号を付す。
【０１１５】
本例は、前述した図１の方位角計算用フィルタ部５の応用例を示す。ここでは、自動車内に方位角計測装置１００としての携帯機器を設置して地磁気を測定する場合、方位角を求めるのに適したフィルタとして、２つの例を示す。なお、方位角計測装置１００の基本的な構成は、図１と同様であり、ここでの説明は省略する。
【０１１６】
以下に示す２つのフィルタの例は、単独で構成してもよいし、任意の順番に並べて構成してもよい。以下の説明では、フィルタの入力は全て地磁気センサ１で測定された値としているが、フィルタを複数個使用する場合は、次のフィルタの入力は、一つ前のフィルタの出力とする。
【０１１７】
（方位角用変化量制限フィルタ）
図１２は、方位角計算用フィルタ部５を方位角用変化量制限フィルタにより構成した場合におけるフィルタ処理のフローチャートを示す。
【０１１８】
2006337057.tif 000023
【０１１９】
2006337057.tif 000024
【０１２０】
2006337057.tif 000025
【０１２１】
（方位角用領域内平均フィルタ）
図１３は、方位角計算用フィルタ部５を方位角用領域内平均フィルタにより構成した場合におけるフィルタ処理のフローチャートを示す。
【０１２２】
測定が開始される前に(一番最初に図８のフローが実行されるより前に)フィルタに入力されたデータ数を表すＮは０に初期化されているとする。
【０１２３】
2006337057.tif 000026
【０１２４】
2006337057.tif 000027
【０１２５】
2006337057.tif 000028
【０１２６】
2006337057.tif 000029
【０１２７】
2006337057.tif 000030
【０１２８】
2006337057.tif 000031
【０１２９】
2006337057.tif 000032
【０１３０】
2006337057.tif 000033
【０１３１】
2006337057.tif 000034
【０１３２】
ここで、方位角用領域内平均フィルタについて補足説明する。
【０１３３】
方位角は、通常定期的に（普通は地磁気の測定周期毎に、あるいは、その整数倍）算出される。オフセット計算のための特定領域平均化フィルタを用いると、測定値が特定の範囲を飛び出さない限り、出力が出てこないので、若干アルゴリズムを修正する必要がある。方位角計算用特定領域平均化フィルタは、地磁気測定値が特定の範囲の領域に入っている場合は、その時点での特定領域内のデータの平均を出力し、測定値が飛び出した場合は、測定値をそのまま出力する。上述のフィルタを適用することによって、地磁気の測定毎に、ノイズが削減された出力が得られる。また、出力は測定値の変化に追随する。
【図面の簡単な説明】
【０１３４】
【図１】本発明の第１の実施の形態である、方位角計測装置の構成を示すブロック図である。
【図２】歩行時のオフセット推定処理を示すフローチャートである。
【図３】自動車内設置時でのオフセット推定処理を示すフローチャートである。
【図４】方位角計測装置を、地磁気センサおよび傾斜角センサを搭載した携帯機器として構成した場合の斜視図である。
【図５】携帯機器の回転軸に対する姿勢を示す説明図である。
【図６】本発明の第２の実施の形態である、オフセット用変化量制限フィルタにより構成した場合のフィルタ処理を示すフローチャートである。
【図７】地磁気測定値およびフィルタ処理後の地磁気データの変化量の様子を示す説明図である。
【図８】オフセット用領域内平均フィルタにより構成した場合のフィルタ処理を示すフローチャートである。
【図９】特定領域の領域設定例を示す説明図である
【図１０】特定領域の領域設定例を示す説明図である
【図１１】オフセット用平面データ選択フィルタにより構成した場合のフィルタ処理を示すフローチャートである。
【図１２】本発明の第３の実施の形態である、方位角用変化量制限フィルタにより構成した場合のフィルタ処理を示すフローチャートである。
【図１３】方位角用領域内平均フィルタにより構成した場合のフィルタ処理を示すフローチャートである。
【符号の説明】
【０１３５】
１  ３軸地磁気検出部
２  オフセット計算用フィルタ部
３  データ取捨選択部
４  オフセット算出部
５  方位角計算用フィルタ部
６  方角計算部
７  地磁気オフセット格納部
８  携帯端末設置姿勢格納部
９  傾斜角センサ
１００  方位角計測装置
２００，２０１  特定領域
【図１】
2006337057.tif 000035
【図２】
2006337057.tif 000036
【図３】
2006337057.tif 000037
【図４】
2006337057.tif 000038
【図５】
2006337057.tif 000039
【図６】
2006337057.tif 000040
【図７】
2006337057.tif 000041
【図８】
2006337057.tif 000042
【図９】
2006337057.tif 000043
【図１０】
2006337057.tif 000044
【図１１】
2006337057.tif 000045
【図１２】
2006337057.tif 000046
【図１３】
2006337057.tif 000047