(57)【要約】      （修正有）【課題】３次元地磁気センサから出力される地磁気データのオフセット誤差を補正する磁気データ処理装置を提供する。【解決手段】３次元磁気センサ２０から出力される磁気データを順次取得しながら前記磁気データを母集団として繰り返し蓄積する蓄積手段９１と、前記母集団が蓄積される度に、球面の方程式から導かれる目的関数ｆ（ｐ）の最小値を前記母集団の信頼指数Ｓとして導出し、前記母集団を十分信頼できるか否かを前記信頼指数Ｓを用いて判定する判定手段９２と、前記母集団を十分信頼できる場合、前記母集団に基づいて前記磁気データのオフセットを導出するオフセット導出手段９４a、９４bと、を備える。【選択図】図１
2011022073.tif 000002
【特許請求の範囲】
【請求項１】
３次元磁気センサから出力される磁気データを順次取得しながら前記磁気データを母集団ｑ１，・・・，ｑＮ（Ｎ≧４）として繰り返し蓄積する蓄積手段と、
前記母集団が蓄積される度にｆ（ｐ）の最小値を前記母集団の信頼指数Ｓとして導出し、前記母集団を十分信頼できるか否かを前記信頼指数Ｓを用いて判定する判定手段と、
前記母集団を十分信頼できる場合、前記母集団に基づいて前記磁気データのオフセットを導出するオフセット導出手段と、
を備え、
【数１８】
2011022073.tif 000020
である磁気データ処理装置。
【請求項２】
前記判定手段は、対称行列Ａの固有値をλ１、λ２、λ３（λ１≧λ２≧λ３）とし、λ１、λ２、λ３にそれぞれ対応する大きさ１の固有ベクトルをｕ１，ｕ２，ｕ３、αを所定の閾値とするとき、
【数１９】
2011022073.tif 000021
として前記信頼指数Ｓを導出し、
【数２０】
2011022073.tif 000022
である請求項１に記載の磁気データ処理装置。
【請求項３】
前記判定手段は、前記母集団が立体的に拡散しているか否かを閾値を用いて判定し、
前記オフセット導出手段は、前記母集団が立体的に拡散している場合、ｆ（ｐ）が最小値をとるｐを前記オフセットとして導出する、
請求項１又は２に記載の磁気データ処理装置。
【請求項４】
前記判定手段は、前記母集団が平面的に拡散しているか否かを閾値を用いて判定し、
前記オフセット導出手段は、前記母集団が平面的に拡散している場合、対称行列Ａの固有値をλ１、λ２、λ３（λ１≧λ２≧λ３）とし、λ１、λ２、λ３にそれぞれ対応する大きさ１の固有ベクトルをｕ１，ｕ２，ｕ３、前回導出した前記オフセットをｐ０とするとき、
ｐ＝ｐ0＋β1ｕ1＋β2ｕ2（β1、β2は実数）
の制約条件の下でｆ（ｐ）が最小値をとるｐを前記オフセットとして導出し、
Ａ＝ＸＴＸ
である請求項１又は２に記載の磁気データ処理装置。
【請求項５】
前記３次元磁気センサをさらに備える、
請求項１から４のいずれか一項に記載の磁気データ処理装置。
【請求項６】
３次元磁気センサから出力される磁気データを順次取得しながら前記磁気データを母集団ｑ１，・・・，ｑＮ（Ｎ≧４）として繰り返し蓄積し、
前記母集団が蓄積される度にｆ（ｐ）の最小値を前記母集団の信頼指数Ｓとして導出し、前記母集団を十分信頼できるか否かを前記信頼指数Ｓを用いて判定し、
前記母集団を十分信頼できる場合、前記母集団に基づいて前記磁気データのオフセットを導出する、
ことを含み、
【数２１】
2011022073.tif 000023
である磁気データ処理方法。
【請求項７】
３次元磁気センサから出力される磁気データを順次取得しながら前記磁気データを母集団ｑ１，・・・，ｑＮ（Ｎ≧４）として繰り返し蓄積する蓄積手段と、
前記母集団が蓄積される度にｆ（ｐ）の最小値を前記母集団の信頼指数Ｓとして導出し、前記母集団を十分信頼できるか否かを前記信頼指数Ｓを用いて判定する判定手段と、
前記母集団を十分信頼できる場合、前記母集団に基づいて前記磁気データのオフセットを導出するオフセット導出手段と、
してコンピューターを機能させ、
【数２２】
2011022073.tif 000024
である磁気データ処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は磁気データ処理装置、磁気データ処理方法及び磁気データ処理プログラムに関し、特に３次元地磁気データのオフセット補正に関する。
【背景技術】
【０００２】
従来、携帯型電話機等の移動体に搭載され、地磁気の方向を検出する３次元地磁気センサが知られている。磁気データ処理装置が磁気データに基づいて方位を導出するとき、移動体の着磁による測定誤差を打ち消すために磁気データを補正する処理が必要である。この補正処理の操作値は、複数の磁気データに基づいて導出され、オフセットと呼ばれる（特許文献１，２，３参照）。
【０００３】
ところで、磁気データには移動体の着磁によるオフセット成分以外にも、移動体に搭載されている電子回路が発生源となる磁気や、ガウス分布に従う磁気センサの出力のゆらぎなどによる誤差成分が含まれている。したがって、移動体の着磁による真のオフセットを正確に導出するためには、広い範囲に立体的に拡散した磁気データを蓄積することが必要となる。しかし、このような磁気データが蓄積されるような移動体の運動は特殊であるから、一般には、移動体に特殊な運動をさせるための操作がユーザーに要求される。そして、真のオフセットは不定期に変動するため、ユーザーが頻繁にそのような操作をしない限り、磁気データ処理装置によって導出される方位と実際の方位とがずれる可能性が高まる。これまでに本件発明者は、ユーザーに特殊な操作を要求することなく、真のオフセットに近いオフセットを導出できる発明を創作している（特許文献１参照）。特許文献１には、母集団としての磁気データの分布が立体的である場合には母集団に基づいて３次元的にオフセットを導出し、母集団としての磁気データの分布が平面的である場合には２次元的にオフセットを導出する方法が開示されている。この方法は、磁気データが立体的に拡散していない場合であっても、磁気データが平面的に拡散している場合であれば、磁気データが拡散している平面に平行な方向について、過去に導出されたオフセットを２次元的に補正するものである。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００７−２４０２７０号公報
【特許文献２】特開２００７−１０７９２１号公報
【特許文献３】特開２００７−１３９７１５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、特許文献１、２，３のいずれに記載の方法においても、母集団がオフセットを導出するために十分信頼できるか否かが判定される。特許文献１，２，３に記載の方法では、その判定に用いる指数として、母集団に基づいて導出したオフセット候補と母集団を構成しているそれぞれの磁気データとの距離の平均ｒの関数を用いている。しかし、特許文献１に記載の方法によると、過去に導出されたオフセットを２次元的に補正したオフセット候補の位置は、過去に導出されたオフセットの位置に依存する。したがって、オフセット候補が導出される母集団が同じであったとしても、過去に導出されたオフセットの位置が異なれば、母集団に基づいて導出したオフセット候補と母集団を構成しているそれぞれの磁気データとの距離の平均ｒも異なることになる。その結果、オフセット候補が導出される母集団が同じであったとしても、その母集団がオフセットを導出するために十分信頼できるかを示す指数が異なることになる。また、特許文献１，２，３に記載の方法では、その母集団がオフセットを導出するために十分信頼できるか否かを判定するには、オフセット候補を導出する必要があるため、オフセット候補が棄却される場合にはオフセット候補を導出したことが無駄になる。
本発明は、磁気データのオフセットを導出するために用いる母集団の信頼性判定に一貫性を持たせるとともに、信頼性判定に費やされる計算量を低減することを目的の１つとする。
【課題を解決するための手段】
【０００６】
（１）上記目的を達成するための磁気データ処理装置は、３次元磁気センサから出力される磁気データを順次取得しながら前記磁気データを母集団ｑ１，・・・，ｑＮ（Ｎ≧４）として繰り返し蓄積する蓄積手段と、前記母集団が蓄積される度にｆ（ｐ）の最小値を前記母集団の信頼指数Ｓとして導出し、前記母集団を十分信頼できるか否かを前記信頼指数Ｓを用いて判定する判定手段と、前記母集団を十分信頼できる場合、前記母集団に基づいて前記磁気データのオフセットを導出するオフセット導出手段と、を備える。ここで、
【数１】
2011022073.tif 000003
である。
関数ｆ（ｐ）の最小値は母集団が近傍に集中している球面を基準にして母集団がばらついているほど大きくなる。したがって関数ｆ（ｐ）の最小値を信頼指数Ｓとして用いて母集団の信頼性を判定することができる。そして本発明によると、過去に導出されたオフセットと無関係に母集団によって決まるＸおよびｊを定数とする関数ｆ（ｐ）の最小値を用いて母集団の信頼性を判定するため、母集団の信頼性判定に一貫性がある。また本発明によると、オフセット候補を導出することなく母集団の信頼性を判定できるため、母集団の信頼性判定に費やされる計算量を低減できる。
【０００７】
（２）上記目的を達成するための磁気データ処理装置において、前記判定手段は、対称行列Ａの固有値をλ１、λ２、λ３（λ１≧λ２≧λ３）とし、λ１、λ２、λ３にそれぞれ対応する大きさ１の固有ベクトルをｕ１，ｕ２，ｕ３、αを所定の閾値とするとき、
【数２】
2011022073.tif 000004
として前記信頼指数Ｓを導出してもよい。ここで、
【数３】
2011022073.tif 000005
である。
【０００８】
（３）上記目的を達成するための磁気データ処理装置において、前記判定手段は、前記母集団が立体的に拡散しているか否かを閾値を用いて判定し、前記オフセット導出手段は、前記母集団が立体的に拡散している場合、ｆ（ｐ）が最小値をとるｐを前記オフセットとして導出してもよい。
本発明によると、母集団の信頼性を判定するために用いる関数に現れる係数と、オフセットを導出するために用いる連立方程式の係数とが一致するため、オフセットを導出するために費やされる計算量を低減できる。
（４）上記目的を達成するための磁気データ処理装置において、前記判定手段は、前記母集団が平面的に拡散しているか否かを閾値を用いて判定し、前記オフセット導出手段は、前記母集団が平面的に拡散している場合、対称行列Ａの固有値をλ１、λ２、λ３（λ１≧λ２≧λ３）とし、λ１、λ２、λ３にそれぞれ対応する大きさ１の固有ベクトルをｕ１，ｕ２，ｕ３、前回導出した前記オフセットをｐ０とするとき、
ｐ＝ｐ0＋β1ｕ1＋β2ｕ2（β1、β2は実数）
の制約条件の下でｆ（ｐ）が最小値をとるｐを前記オフセットとして導出してもよい。ここでＡ＝ＸＴＸである。
本発明によると、母集団の信頼性を判定するために用いる関数に現れる係数と、オフセットを導出するために用いる連立方程式の係数とが一致するため、オフセットを導出するために費やされる計算量を低減できる。また、母集団が三次元的に拡散していない場合であっても、二次元的に拡散していれば真のオフセットにより近いオフセットを導出することができる。
【０００９】
なお、請求項に記載された各手段の機能は、構成自体で機能が特定されるハードウェア資源、プログラムにより機能が特定されるハードウェア資源、又はそれらの組み合わせにより実現される。また、これら各手段の機能は、各々が物理的に互いに独立したハードウェア資源で実現されるものに限定されない。さらに、本発明は磁気データ処理方法としても、磁気データ処理プログラムとしても、磁気データ処理プログラムの記録媒体としても成立する。むろん、そのコンピュータプログラムの記録媒体は、磁気記録媒体であってもよいし光磁気記録媒体であってもよいし、今後開発されるいかなる記録媒体であってもよい。
【図面の簡単な説明】
【００１０】
【図１】本発明の実施形態にかかるブロック図。
【図２】本発明の実施形態にかかるフローチャート。
【図３】本発明の実施形態にかかる模式図。
【発明を実施するための形態】
【００１１】
以下、本発明の実施の形態を添付図面を参照しながら説明する。尚、各図において対応する構成要素には同一の符号が付され、重複する説明は省略される。
１．第一実施形態
（概要）
本発明の第一実施形態では、母集団としての磁気データの分布が立体的である場合には母集団に基づいて３次元的にオフセットが導出される。一方、母集団としての磁気データの分布が平面的である場合には過去に導出されたオフセット（旧オフセット）を母集団に基づいて２次元的に補正することによってオフセットが導出される。そして、母集団としての磁気データが近傍に集中している球面を基準としてばらついている場合には、その母集団を用いてオフセットを導出する処理は実行されない。
【００１２】
（磁気データ処理装置の構成）
図１は本発明の磁気データ処理装置の一実施形態を示すブロック図である。磁気データ処理装置１は携帯電話機、ＰＮＤ（Personal Navigation Device）、電子コンパス、デジタルカメラなどの携帯型情報端末に搭載される。
【００１３】
磁気データ処理装置１は磁気センサ２０とマイクロコンピューターとで構成されている。マイクロコンピューターは、ＣＰＵ４０、ＲＯＭ４２、ＲＡＭ４４、入出力機構（Ｉ／Ｏ）３０などからなり、磁気センサ２０から出力される磁気データを入力し、オフセット補正された磁気データに基づいて進行方位や走行予定経路を画像や音声によってユーザーに報知するための方位データを出力する。表示部６０は方位を示す画像を表示するためのディスプレイとディスプレイドライバーとを含む。
【００１４】
磁気センサ２０は、磁界ベクトルの直交する３軸の成分をそれぞれ検出するｘ軸センサ２１とｙ軸センサ２２とｚ軸センサ２３とインターフェース２４とを備えた３次元磁気センサである。ｘ軸センサ２１、ｙ軸センサ２２、ｚ軸センサ２３は、いずれも磁気抵抗素子、ホール素子等で構成され、指向性のある１次元磁気センサであればどのようなものであってもよい。ｘ軸センサ２１、ｙ軸センサ２２およびｚ軸センサ２３は、それぞれの感度方向が互いに直交するように固定されている。インターフェース２４は、ｘ軸センサ２１、ｙ軸センサ２２およびｚ軸センサ２３の出力を時分割して取得し、ｘ軸センサ２１、ｙ軸センサ２２およびｚ軸センサ２３からの入力を増幅してＡＤ変換し、磁気データｑ＝（ｑｘ、ｑｙ、ｑｚ）を出力する。インターフェース２４から出力されるディジタル信号である磁気データｑはＲＡＭ４４の所定のアドレスに格納される。不揮発性の記憶媒体であるＲＯＭ４２には、ＲＡＭ４４にロードされ、ＣＰＵ４０によって実行される磁気データ処理プログラム９０や、携帯型情報端末の機能を実現するための種々のプログラムが格納されている。
【００１５】
磁気データ処理プログラム９０は、磁気データに基づいて方位データを出力するためのプログラムであって、ＲＯＭ４２に格納されている。方位データは地磁気の方向を示すベクトルデータである。磁気データ処理プログラム９０は、蓄積モジュール９１、判定モジュール９２、３次元オフセット更新モジュール９４ａ、２次元オフセット更新モジュール９４ｂ、方位導出モジュール９５等のモジュール群で構成されている。
【００１６】
蓄積モジュール９１は、３次元磁気センサから出力される磁気データを順次取得しながら４つ以上のＮ個の磁気データｑ１、・・・、ｑＮを母集団として繰り返し蓄積する機能を実現するプログラム部品である。すなわち蓄積モジュール９１は、磁気センサ２０から順次出力される磁気データを一定の時間間隔で順次取得し、取得した磁気データをバッファに蓄積する機能を実現する。バッファは、Ｎ個の磁気データを記憶するためのＲＡＭ４４の記憶領域である。
【００１７】
判定モジュール９２は、母集団が蓄積される度に母集団の信頼性を判定し、母集団が立体的に拡散しているか否かを判定し、母集団が平面的に拡散しているか否かを判定する機能を実現するプログラムモジュールである。
【００１８】
本実施形態において、判定モジュール９２は、母集団ｑ１、・・・、ｑＮの分布の主値に相当する対称行列Ａの固有値をλ１、λ２、λ３（λ１≧λ２≧λ３）、固有値λ１、λ２、λ３にそれぞれ対応する固有ベクトルをｕ１，ｕ２，ｕ３、α１を予め決められた閾値とするとき、次式（１）によって定義されるＳを母集団の信頼性を示す信頼指数とし、信頼指数Ｓと予め決められた閾値α２がＳ＞α２を満たす場合に、母集団を十分信頼できると判定する。
【数４】
2011022073.tif 000006
ただし、
【数５】
2011022073.tif 000007
である。なお、行列Ａは、次式のようにも書けるため、母集団の分散共分散行列のＮ倍に相当する。
【数６】
2011022073.tif 000008
【００１９】
ここで、母集団が近傍に集中している球面を基準にしてどの程度ばらついているかを信頼指数Ｓが表すことについて、その理由を説明する。
母集団が同一平面上にない４つの磁気データｑ１、・・・、ｑ４で構成されている場合、母集団を含む球面は統計的手法を用いることなく一意に特定される。この球面の中心の位置ベクトルｐ＝（ｐx、ｐy、ｐｙ）と４つの磁気データｑ１、・・・、ｑ４との間には次式（１１）が成り立つ。したがってこの球面の中心は式（１１）を変形した連立方程式（１２）を解くことによって得られる。尚、３変数に対して等号制約が４つあるが、等号制約の１つは冗長になっているため方程式（１２）は必ず解を持つ。
【数７】
2011022073.tif 000009
【００２０】
母集団が５つ以上の磁気データで構成されている場合、ｐについての連立一次方程式（１４）が解を持てば、その解ｐ＝（ｐx、ｐy、ｐｙ）は、母集団データ群を含む球面の中心である。
Ｘｐ=ｊ・・・（１４）
【００２１】
しかし、磁気センサ２０自体の測定誤差を考慮すると、現実には、方程式（１４）が解を持つことはまれである。そこで、統計的な手法により尤もらしい解を得るために、次式（１５）で定義されるベクトルｅを導入する。
ｅ＝Ｘｐ−ｊ・・・（１５）
【００２２】
||ｅ||22（すなわちｅTｅ）を最小にするｐが、母集団との距離が最も近い球面、すなわち母集団が近傍に集中している球面の中心であるといえる。||ｅ||22を最小にするｐを求める問題は、行列Ａが正則のときには次式（１６）の目的関数を最小にする最適化問題となる。
【数８】
2011022073.tif 000010
【００２３】
そして、ｆ（ｐ）は、母集団が近傍に集中する特定の球面を基準として母集団がばらついているほど大きくなり、母集団を構成している磁気データが全てその球面に含まれる場合には０となる。したがって、ｆ（ｐ）の最小値は、母集団が近傍に集中する球面を基準として母集団がばらついている程度を示すことになる。すなわち、ｆ（ｐ）の最小値は母集団の信頼指数Ｓとして用いることができる。ここで式（７）のようにｂ、ｃを定義すると、式（１６）は、次式（１７）のように２次形式で表現できる。
【数９】
2011022073.tif 000011
【００２４】
対称行列Ａの固有値が０を含む場合であっても、式（４）のように行列Ａを定義し、式（７）のようにｂを定義する場合には、目的関数ｆ（ｐ）は最小値を必ず持つ。ここで対称行列Ａの固有値λ１、λ２、λ３（λ１≧λ２≧λ３）にそれぞれ対応する大きさ１の固有ベクトルをｕ１，ｕ２，ｕ３とし、Ｌを式（１０）によって定義すると、対称行列Ａは式（１８）のように対角化できる。
【数１０】
2011022073.tif 000012
【００２５】
またｂについてもＬを用いて式（８）の座標変換を施すと、対称行列Ａが正定値のとき（固有値が０を含まないとき）には目的関数ｆ（ｐ）の最小値ｍｉｎは式（１９）と表される。
【数１１】
2011022073.tif 000013
そして式（９）のようにｂ'を表して式（１９）を書き直すと、目的関数ｆ（ｐ）の最小値ｍｉｎは次式（２０）によって表される。
【数１２】
2011022073.tif 000014
【００２６】
コンピューターによる計算誤差を考慮すると、対称行列Ａが正定値であるという条件は、０に極めて近い閾値α１を定めることによって次式（２１）と表される。
λ３＞α１・・・（２１）
【００２７】
次に対称行列Ａが重複していない固有値０を持つ場合、すなわちλ１≧λ２＞０、λ３＝０の場合について考える。この場合、証明は省略するが、式（２０）の括弧中の第３項は０になる。なお、λ３が極めて小さく、かつ式（２１）が満たされる場合には、信頼指数Ｓに与えるλ３の影響は無視できないほど大きい。したがって、閾値α１は０で除する計算において例外処理が発生しない範囲での最小値に設定することが好ましい。
以上述べたように、式（１６）で定義される目的関数ｆ（ｐ）の最小値を母集団の信頼指数Ｓとするならば、信頼指数Ｓは式（１）のように表されることになる。すなわち、判定モジュール９２は、信頼指数Ｓと予め決められた閾値α２との大小関係を比較し、閾値α２よりも信頼指数Ｓが大きい場合には、母集団を十分信頼できると判定する。
【００２８】
また本実施形態においては、母集団が立体的に拡散している程度を示す指数をλ３／λ１とする。また母集団が平面的に拡散している程度を示す指数をλ２／λ１とする。そして判定モジュール９２は、λ３／λ１が予め決められた閾値を越える場合には、母集団が立体的に拡散していると判定する。また判定モジュール９２は、λ３／λ１がその閾値を越えずλ２／λ１が予め決められた別の閾値を越える場合には、母集団が平面的に拡散していると判定する。
【００２９】
三次元オフセット更新モジュール９４ａは、母集団に基づいて３次元的にオフセットを導出する機能を実現するプログラム部品である。具体的には式（１６）の目的関数ｆ（ｐ）を最小にするｐを求めることにより、オフセットとしてのｐが導出される。目的関数ｆ（ｐ）を最小にするｐは、本実施形態で想定している（ＸＴＸ）が正則のときは式（２２）のように書くことができる。
【数１３】
2011022073.tif 000015
【００３０】
二次元オフセット更新モジュール９４ｂは、母集団に基づいて二次元的に過去のオフセット（旧オフセット）を補正することによって最新のオフセットを導出する機能を実現するプログラム部品である。すなわち、二次元オフセット更新モジュール９４ｂは、蓄積モジュール９１によって保持されている母集団と、過去のオフセットとに基づいて最新のオフセットを導出し、過去のオフセットを最新のオフセットに更新する機能を実現する。具体的には次の通りである。
【００３１】
母集団が特定の平面の近傍に集中して拡散している場合（分布が二次元的である場合）、旧オフセットに対する補正方向をその平面に平行で互いに直交する２方向に制限して新オフセットが導出される。母集団がある特定の平面近傍に集中し、その平面に垂直な方向から見て拡散している場合、その平面と平行な方向については母集団を十分信頼できる一方でその平面に垂直な方向についてはオフセット補正のために用いるデータ群として母集団を信頼できないことになる。このような場合には、その平面に垂直な方向については旧オフセットを補正しないことにより、信頼に値しない情報に基づいてオフセットが更新されることを防止できる。
【００３２】
母集団がある特定の平面近傍に集中し、その平面に垂直な方向から見て拡散している場合、その平面に垂直な方向は、最小固有値λ3に対応する固有ベクトルｕ3の方向に一致し、その平面に平行で互いに直交する方向は、最大固有値λ1、中間固有値λ2にそれぞれ対応するｕ1、ｕ2の方向に一致する。したがって、その平面に垂直な方向について旧オフセットｐ０を補正せずに新オフセットｐを導出するため、次式（２３）の制約条件のもとで式（１６）の目的関数が最小になる新オフセットｐを求める。
ｐ＝ｐ0＋β1ｕ1＋β2ｕ2（β1、β2は実数）・・・（２３）
なお、式（２３）は次式（２４）と等価である。
【数１４】
2011022073.tif 000016
【００３３】
式（２１）の制約条件のもとで式（１６）の最適化問題を解く式は、ラグランジュの未定乗数法によって等価な連立方程式に変形できる。未定乗数ρを導入し、
【数１５】
2011022073.tif 000017
とすると、ｘの連立一次方程式（２６）がその方程式となる。
【数１６】
2011022073.tif 000018
【００３４】
方位導出モジュール９５は、磁気センサ２０から順次取得する磁気データを最新のオフセットを用いて補正して方位データを生成するプログラム部品である。具体的には、方位導出モジュール９５は、ベクトルデータである磁気データの各成分からオフセットの各成分を引き算して得られるベクトルデータを方位データとして出力する。
【００３５】
（磁気データ処理方法）
次に磁気データ処理プログラム９０を実行することによって実現される磁気データ処理方法について図２を参照しながら説明する。図２に示す処理は、一定時間間隔（例えば１０ミリ秒間隔）で繰り返し実行される。
【００３６】
はじめに、磁気センサ２０から出力された磁気データが蓄積モジュール９１によって取得され、バッファに蓄積される（Ｓ１００）。
【００３７】
次に、予め決められた数の磁気データが母集団としてバッファに蓄積されているか否かが判定モジュール９２によって判定される（Ｓ１０１）。予め決められた数の磁気データが母集団としてバッファに蓄積されていない場合、所定時間後に次の磁気データがバッファに蓄積される（Ｓ１００）。
【００３８】
予め決められた数の磁気データが母集団としてバッファに蓄積されている場合、母集団を十分信頼できるか否かが判定モジュール９２によって判定される（Ｓ１０２）。すなわち、式（１）で表される信頼指数Ｓが導出され、信頼指数Ｓと予め決められた閾値α２とがＳ＞α２を満たす場合、母集団は十分信頼できると判定され、Ｓ＞α２を満たさない場合、母集団は十分信頼できないと判定される。
【００３９】
母集団が十分信頼できる場合、母集団が立体的に拡散しているかが判定モジュール９２によって判定される（Ｓ１０３）。すなわち、λ３／λ１が予め決められた閾値を越える場合には、母集団が立体的に拡散していると判定され、λ３／λ１が予め決められた閾値を越えない場合には、母集団が立体的に拡散していないと判定される。
【００４０】
母集団が立体的に拡散している場合、母集団のみに基づいてオフセットが３次元的に更新される（Ｓ１０４）。すなわち、バッファに記憶されている磁気データを母集団として式（２２）によって表されるｐが３次元オフセット更新モジュール９４ａによって最新のオフセットとして導出され、過去のオフセットが最新のオフセットに更新される。
【００４１】
母集団が立体的に拡散していない場合、母集団が平面的に拡散しているか否かが判定モジュール９２によって判定される（Ｓ１０５）。すなわち、λ３／λ１が予め決められた閾値を越えずλ２／λ１が予め決められた別の閾値を越える場合には、母集団は平面的に拡散していると判定され、λ３／λ１も予め決められた閾値を越えずλ２／λ１も予め決められた別の閾値を越えない場合には、母集団は平面的に拡散していないと判定される。
【００４２】
母集団が平面的に拡散している場合、旧オフセットを母集団に基づいて二次元的に補正することによってオフセットが更新される（Ｓ１０６）。すなわち、二次元オフセット更新モジュール９２ｂによって、連立一次方程式（２６）の解ｐが最新のオフセットとして導出され、旧オフセットが導出された最新のオフセットに更新される。
【００４３】
母集団が十分信頼できない場合、あるいは、母集団が立体的にも平面的にも拡散していない場合、蓄積モジュール９１によって母集団が破棄される（Ｓ１０７）。
【００４４】
以上の処理では、最新のオフセットｐを導出する前に母集団の信頼性を判定できるため、オフセット候補を導出した後に母集団の信頼性を判定する場合に比べると、母集団の信頼性を判定するために費やされる計算量が低減される。
【００４５】
ここで仮に、母集団に基づいて導出したオフセット候補ｐと、母集団を構成しているそれぞれの磁気データとの距離の平均ｒの関数を信頼指数Ｓ'として次式（２９）のように定めるとする。
【数１７】
2011022073.tif 000019
【００４６】
いま、図３に黒点で示すように母集団が立体的に拡散しておらず平面的に拡散しているとし、上述したように旧オフセットを直交する２方向についてのみ補正することによってオフセット候補を導出するとする。すると、旧オフセットがｐ０ａの場合、オフセット候補はｐａとなり、旧オフセットがｐ０ｂの場合、オフセット候補はｐｂとなる。オフセット候補ｐａ、ｐｂは、いずれも母集団が近傍に集中している円Ｃの中心Ｏ'を通り、円Ｃを含む平面に対して垂直な直線上にあるが、２つの旧オフセットがその直線に対して垂直な面上にない限り、一致しない。したがって、母集団を構成しているそれぞれの磁気データとオフセット候補との距離の平均ｒａ、ｒｂの関数である信頼指数Ｓ'は評価対象が同一であるにも関わらず、旧オフセットの位置に応じて異なる値となる。このため、信頼指数Ｓ'と閾値との大小関係に応じて母集団の信頼性を判定するならば、母集団が近傍に集中する平面と旧オフセットとの距離が取り得る範囲を加味して閾値を設定しなければならない。閾値を緩く（小さく）設定すれば、母集団がノイズ成分を多く含んでばらついている場合であっても、母集団が近傍に集中している平面と旧オフセットの距離が遠ければ、母集団は信頼できると判定されることになる。一方、閾値を厳しく（大きく）設定すれば、母集団がノイズ成分をふくんでいなくても、母集団が近傍に集中している平面と旧オフセットの距離が近ければ、母集団は信頼できないと判定されることになる。
これに対し、本実施形態においては、母集団の信頼性は、式（１）で表されるように母集団のみの関数である信頼指数Ｓを用いて判定されるため、母集団が同じであれば、母集団の信頼指数Ｓは常に一致する。したがって、信頼指数Ｓと比較される閾値α２を、旧オフセットを考慮した信頼指数Ｓ'のばらつきと無関係に厳密に定めることができる。このため、本実施形態によると母集団の信頼性を精度良く判定できる。
【００４７】
２．他の実施形態
尚、本発明の技術的範囲は、上述した実施の形態に限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変更を加え得ることは勿論である。
例えば、母集団に基づいて３次元的にオフセットを導出する手段も、母集団に基づいて２次元的に旧オフセットを補正する手段も、上記実施形態に開示した手段に限定されず、どのような公知の手段を用いてもよいし、今後開発されるいかなる手段を用いても良い。公知の手段としては、統計的手法を用いて２次元的または３次元的にオフセットを導出するものや、３つまたは４つの磁気データを母集団から選抜し、選抜された母集団からオフセットとしての球面の中心または旧オフセットを補正するための仮オフセットとしての円の中心を求めるもの等がある。尚、磁気データのオフセットを導出する方法は、本件発明者によっても複数提案されている（特開２００７−２４０２７０号公報、特開２００７−２０５９４４号公報、特開２００７−１３９７１５号公報、特開２００７−１０７９２１号公報、特願２００７−３３９４７８など）。
【００４８】
また例えば、母集団が立体的に拡散しているか平面的に拡散しているかを判定する手段は、上記実施形態に開示した手段に限定されず、どのような公知の手段を用いてもよいし、今後開発されるいかなる手段を用いても良い。公知の手段としては、母集団の分散の大きさを閾値と比較する手段などがある。
【００４９】
また例えば、母集団に基づいてオフセットを導出した後に、さらに別の指数と閾値とを用いてそのオフセットの信頼性を判定するとともに信頼性が低い場合には、導出したオフセットを棄却してもよい。
【符号の説明】
【００５０】
１：磁気データ処理装置、２０：磁気センサ、２１：ｘ軸センサ、２２：ｙ軸センサ、２３：ｚ軸センサ、２４：インターフェース、３０：Ｉ／Ｏ、４０：ＣＰＵ、４２：ＲＯＭ、４４：ＲＡＭ、６０：表示部、９０：磁気データ処理プログラム、９１：蓄積モジュール、９２：判定モジュール、９４ａ：３次元オフセット更新モジュール、９４ｂ：二次元オフセット更新モジュール、９５：方位導出モジュール
【図１】
2011022073.tif 000025
【図２】
2011022073.tif 000026
【図３】
2011022073.tif 000027