(57)【要約】【課題】磁気データ群の分布に関わらず、蓄積された磁気データ群からオフセットを簡素な処理で補正できる磁気データ処理装置、方法及びプログラム、並びに磁気測定装置を提供する。【解決手段】磁気センサから出力される磁気データｑｉを順次入力する入力手段（９２）と、前記入力手段によって入力される前記磁気データを順次蓄積する蓄積手段（９２）と、前記蓄積手段に蓄積された複数の前記磁気データを母集団とする確率論的手法によりオフセットを導出するオフセット導出手段（９４）と、を備え、前記オフセット導出手段（９４）は、前記各磁気データに含まれる誤差Δｐｉを確率分布μに従う確率変数ｖｉとみなし、確率分布μに対応する期待値を用いて前記オフセットを導出する、磁気データ処理装置。【選択図】図４
2011203243.tif 000002
【特許請求の範囲】
【請求項１】
磁気センサから出力される磁気データｑｉを順次入力する入力手段と、
前記入力手段によって入力される前記磁気データを順次蓄積する蓄積手段と、
前記蓄積手段に蓄積された複数の前記磁気データを母集団とする確率論的手法によりオフセットを導出するオフセット導出手段と、を備え、
前記オフセット導出手段は、前記各磁気データに含まれる誤差Δｐｉを確率分布μに従う確率変数ｖｉとみなし、確率分布μに対応する期待値を用いて前記オフセットを導出する、
磁気データ処理装置。
【請求項２】
前記オフセット導出手段は、前記オフセットから前記母集団の誤差を除いた真の磁気データまでの距離のばらつきを示す指標を導出し、当該導出した指標の期待値を最大または最小にする前記オフセットを導出する、
請求項１に記載の磁気データ処理装置。
【請求項３】
前記指標を実数値関数ｆｖ（ｃ）、ｆｖ（ｃ）の期待値をｆ（ｃ）とするとき、
前記オフセット導出手段は、時間について独立なｖｉを含むｆｖ（ｃ）の項をμに対応する期待値によって置換されたｆ（ｃ）の値を最大または最小にするｃを前記オフセットとして導出する、
請求項２に記載の磁気データ処理装置。
【請求項４】
【数４７】
2011203243.tif 000049
ただし、
【数４８】
2011203243.tif 000050
である、
請求項３に記載の磁気データ処理装置。
【請求項５】
前記オフセット導出手段は、Ａ'が正定値のとき、
【数４９】
2011203243.tif 000051
または
【数５０】
2011203243.tif 000052
を前記オフセットとして導出する、ただし、
Ｖはμの分散共分散行列であって、
【数５１】
2011203243.tif 000053
である、
請求項３または４に記載の磁気データ処理装置。
【請求項６】
前記指標を実数値関数ｆｖ（ｃ）、ｆｖ（ｃ）の期待値をｆ（ｃ）とするとき、
前記オフセット導出手段は、時間について独立なｖｉを含むｆｖ（ｃ）の項をμに対応する期待値によって置換されたｆ（ｃ）の値を制約条件付きで最大または最小にするｃを前記オフセットとして導出する、
請求項２に記載の磁気データ処理装置。
【請求項７】
【数５２】
2011203243.tif 000054
ただし、
【数５３】
2011203243.tif 000055
であって、
前記蓄積手段は、前記磁気データの旧オフセットを新オフセットに更新するために複数の前記磁気データを順次蓄積し、
前記オフセット導出手段は、前記旧オフセットを用いずに前記母集団から導出される仮オフセットの前記旧オフセットに対する位置ベクトルであって前記母集団の分布の主軸方向の基本ベクトル群の一次結合である位置ベクトルの各係数を前記母集団の分布の主値に応じて重み付けした値を係数とする前記基本ベクトル群の一次結合を補正ベクトルとするとき、前記旧オフセットと前記補正ベクトルの和となる前記新オフセットを求めることを前記制約条件として、前記旧オフセットと前記母集団とに基づいて前記新オフセットを導出する、
請求項６に記載の磁気データ処理装置。
【請求項８】
磁気センサから出力される磁気データｑｉを順次入力する入力ステップと、
前記入力ステップにおいて入力される前記磁気データを順次蓄積する蓄積ステップと、
前記蓄積ステップにおいて蓄積された複数の前記磁気データを母集団とする確率論的手法によりオフセットを導出するオフセット導出ステップと、を含み、
前記オフセット導出ステップにおいて、前記各磁気データに含まれる誤差Δｐｉを確率分布μに従う確率変数ｖｉとみなし、確率分布μに対応する期待値を用いて前記オフセットを導出する、
磁気データ処理方法。
【請求項９】
磁気センサから出力される磁気データｑｉを順次入力する入力手段と、
前記入力手段によって入力される前記磁気データを順次蓄積する蓄積手段と、
前記蓄積手段に蓄積された複数の前記磁気データを母集団とする確率論的手法によりオフセットを導出するオフセット導出手段と、
してコンピューターを機能させ、
前記オフセット導出手段は、前記各磁気データに含まれる誤差Δｐｉを確率分布μに従う確率変数ｖｉとみなし、確率分布μに対応する期待値を用いて前記オフセットを導出する、
磁気データ処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は磁気データ処理装置、方法及びプログラムに関し、特に磁気センサのオフセットを補正する技術に関する。
【背景技術】
【０００２】
従来、携帯型電話機、車両等の移動体に搭載され、地磁気の方向を検出する磁気センサが知られている。一般に３次元磁気センサは、磁力を互いに直交する３方向の成分に分解するための３つの磁気センサモジュールを備えている。３次元磁気センサの出力である磁気データは、そのような３つの磁気センサモジュールのそれぞれの出力の組み合わせからなり、互いに直交する単位ベクトルの一次結合である３次元ベクトルデータである。磁気データの方向と大きさが、３次元磁気センサが検出している磁力の方向と大きさである。磁気センサの出力に基づいて地磁気の方向又は大きさを特定するとき、移動体の着磁による測定誤差を打ち消すために磁気センサの出力を補正する処理が必要である。この補正処理の操作値はオフセットと呼ばれている。オフセットは磁気センサが検出している移動体の着磁による磁力成分を表すベクトルデータであって、磁気センサの出力である磁気データからオフセットが引き算されることによって測定誤差が打ち消される。３次元磁気センサでは磁気データ群が分布する球面の中心を求めることによって、２次元磁気センサでは磁気データ群が分布する円の中心を求めることによって、オフセットを算出することができる。
【０００３】
特許文献１には、磁気データ群の分布形状に応じてオフセットを補正する方向を制限することによって、磁気データ群の分布が狭い場合でもオフセットを真値に近づけるオフセット補正方法が記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００７−２４０２７０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、仮にオフセットが真値であったとしても、オフセット補正された磁気データ群のそれぞれは誤差を含む。磁気センサの出力自体がガウス分布に従う測定誤差を有しているし、現実には完全に一様な磁界は存在しないためオフセットを算出するために必要な磁気データが蓄積される期間中に磁気センサが測定している磁界が変動するし、磁気センサの出力をディジタル値で取り出すまでの計算誤差があるからである。すなわち、磁気データｑｉ（ｉ＝１，２，・・・）は真値ｐｉにオフセットとは無関係な誤差Δｐｉが重畳された値である。つまり、
ｑｉ＝ｐｉ＋Δｐｉ
が磁気データｑｉとして磁気データから出力される。そして真のオフセットは、真値ｐｉが分布する円または球面の中心に相当する。
しかし、磁気データｑｉを真値ｐｉと誤差Δｐｉとに分離することができないために、従来のオフセット補正方法では誤差Δｐｉが０とはいえないにも関わらず、Δｐｉを０とみなしているため、従来のオフセット補正方法によって導出されるオフセットは真値と異なるものになる。
【０００６】
本発明は、蓄積された磁気データ群からオフセットを正確に導出できる磁気データ処理装置、方法及びプログラム、並びに磁気測定装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
（１）上記目的を達成するための磁気データ処理装置は、磁気センサから出力される磁気データｑｉを順次入力する入力手段と、前記入力手段によって入力される前記磁気データを順次蓄積する蓄積手段と、前記蓄積手段に蓄積された複数の前記磁気データを母集団とする確率論的手法によりオフセットを導出するオフセット導出手段と、を備え、前記オフセット導出手段は、前記各磁気データに含まれる誤差Δｐｉを確率分布μに従う確率変数ｖｉとみなし、確率分布μに対応する期待値を用いて前記オフセットを導出する。
本発明によると、複数の磁気データからなる母集団から確率論的手法によりオフセットを導出するとき、オフセットとは無関係に磁気データに含まれる誤差Δｐｉを確率分布μに従う確率変数ｖｉとみなし、確率分布μに対応する期待値を用いるため、オフセットを正確に導出することができる。具体的には、磁気データｑｉ（ｉ＝１，２，・・・）を真値ｐｉと確率分布μに対応する誤差ｖｉとが重畳された値として扱う確率論的な演算によってオフセットを導出する過程において、確率変数である誤差ｖｉを含む項を確率分布μに対応する期待値に置換する。すると、オフセットとは無関係の誤差ｖｉの傾向（確率分布）を加味した正確なオフセットを導出することができる。
（２）複数の磁気データを母集団とする確率論的手法により磁気センサのオフセットを導出するにあたっては、数理計画問題の解としてオフセットを導出することができる。数理計画問題において目的関数の係数の中に確率的要素が含まれているとき、その問題は確率計画問題といわれる。オフセットを数理計画問題の解として正確に導出するためには誤差Δｐｉを用いて目的関数を定義しなければならない。すなわち磁気センサのオフセットを正確に導出する数理計画問題は確率論的モデルになる。従来、磁気センサのオフセットを導出する数理計画問題は、本来確率論的モデルであるにも関わらず、オフセットと無関係な誤差と真値とを分離できないが故に決定論的モデルとして扱われていた。確率計画問題の解として導出するために定義される目的関数に含まれる確率変数である誤差ｖｉを含む項を確率分布μに対応する期待値によって置換することにより、確率論的モデルを決定論的モデルに置き換えて解を導出することが可能になる。このように、期待値を用いて決定論的モデルに置き換えた数理計画問題の解を磁気センサのオフセットとして導出することによって、確率変数である誤差ｖｉの傾向（確率分布）を加味した正確なオフセットを導出することが可能になる。
そこで上記目的を達成するための磁気データ処理装置において、前記オフセット導出手段は、前記オフセットから前記母集団の誤差を除いた真の磁気データまでの距離のばらつきを示す指標を導出し、当該導出した指標の期待値を最大または最小にする前記オフセットを導出してもよい。
（３）また、前記指標を実数値関数ｆｖ（ｃ）、ｆｖ（ｃ）の期待値をｆ（ｃ）とするとき、前記オフセット導出手段は、時間について独立なｖｉを含むｆｖ（ｃ）の項をμに対応する期待値によって置換されたｆ（ｃ）の値を最大または最小にするｃを前記オフセットとして導出してもよい。
（４）磁気データをｑｉ＝ｐｉ＋ｖｉ（ｉ＝１，２，・・・）とすると、真値ｐｉからオフセットｃまでの距離の２乗であるｄｉ、すなわち母集団が分布する球面または円の半径の２乗は、
【数１】
2011203243.tif 000003
となる。ここでｄｉの平均を、
【数２】
2011203243.tif 000004
と定義すると、次の実数値関数ｆｖ（ｃ）はｄｉのばらつきの指標となる。
【数３】
2011203243.tif 000005
ｄｉのばらつきの指標である関数ｆｖ（ｃ）を最小にするｃがオフセットとして最も確からしいといえる。ここで、
【数４】
2011203243.tif 000006
とすると、
【数５】
2011203243.tif 000007
である。
そこで上記目的を達成するための磁気データ処理装置において、
【数６】
2011203243.tif 000008
としてもよい。
（５）現実の磁気センサのオフセットと無関係な誤差Δｐｉは、平均が０、分散共分散行列がＶで表される正規分布に従い時間について独立な確率変数ｖｉによって近似することができる。
そこで上記目的を達成するための磁気データ処理装置において、誤差ｖｉを、平均が０、分散共分散行列がＶで表される正規分布に従い時間について独立な確率変数とみなし、制約条件無しに目的関数ｆｖ（ｃ）を最小にするｃを求めてもよい。すなわち上記目的を達成するための磁気データ処理装置において、前記導出手段は、Ａ'が正定値のとき、
【数７】
2011203243.tif 000009
または
【数８】
2011203243.tif 000010
を前記オフセットとして導出してもよい。ただし、
Ｖはμの分散共分散行列であって、
【数９】
2011203243.tif 000011
である。
（６、７）既に述べたとおり、次の関数ｆｖ（ｃ）はｄｉのばらつきの指標となる。
【数１０】
2011203243.tif 000012
そしてｄｉのばらつきの指標である関数ｆｖ（ｃ）を最小にするｃを導出することにより正確にオフセットを導出できる。しかし、母集団の分布が特定の平面や特定の直線の近傍に偏って分布している場合には、その平面または直線に対して垂直な方向にはオフセットを補正しない方が正確にオフセットを導出できる可能性が高い。そこで、母集団の分布が特定の平面や特定の直線の近傍に偏って分布している程度に応じて、その平面または直線に対して垂直な方向にオフセットを補正する量を制限してもよい。母集団の分布が特定の平面や特定の直線の近傍に偏って分布している程度を表す指標としては"母集団の分布の主値"を用いることができる。
すなわち、上記目的を達成するための磁気データ処理装置において、前記指標を実数値関数ｆｖ（ｃ）、ｆｖ（ｃ）の期待値をｆ（ｃ）とするとき、
前記オフセット導出手段は、時間について独立なｖｉを含むｆ（ｃ）の項をμに対応する期待値によって置換し、前記期待値によって置換されたｆ（ｃ）の値を制約条件付きで最大または最小にするｃを前記オフセットとして導出してもよい。
さらに、上記目的を達成するための磁気データ処理装置において、
【数１１】
2011203243.tif 000013
ただし、
【数１２】
2011203243.tif 000014
であって、前記蓄積手段は、前記磁気データの旧オフセットを新オフセットに更新するために複数の前記磁気データを順次蓄積し、前記導出手段は、前記旧オフセットを用いずに前記母集団から導出される仮オフセットの前記旧オフセットに対する位置ベクトルであって前記母集団の分布の主軸方向の第二の基本ベクトル群の一次結合である位置ベクトルの各係数を前記母集団の分布の主値に応じて重み付けした値を係数とする前記第二の基本ベクトル群の一次結合を補正ベクトルとするとき、前記旧オフセットと前記補正ベクトルの和となる前記新オフセットを求めることを前記制約条件として、前記旧オフセットと前記母集団とに基づいて前記新オフセットを導出してもよい。
ここで採用しているアルゴリズムを図１を参照しながら詳細に説明する。このアルゴリズムの要点は、分散が大きな主軸方向に分布している磁気データ群ほど、オフセットを更新するために用いる母集団要素として重く評価され、分散が小さな主軸方向に分布している磁気データ群ほど母集団要素として軽く評価されることにある。具体的には次のとおりである。旧オフセットｃ0、新オフセットｃ及び仮オフセット（原点０に対するｇの終点の位置）は、磁気データの基本ベクトル群の一次結合である３次元の位置ベクトルデータである。新オフセットｃは、旧オフセットｃ0を新オフセットｃに更新するために蓄えられた磁気データ群と旧オフセットｃ0とに基づいて導出される。
【０００８】
新オフセットｃに更新するために蓄えられた複数の磁気データからなる母集団は予め決められた期間内に蓄積された磁気データ群で構成されてもよいし、予め決められた個数の磁気データ群で構成されてもよいし、あるタイミング（例えばオフセットの更新要求が発生したタイミング）で蓄積されている任意の個数の磁気データ群で構成されてもよい。
旧オフセットｃ0は、新オフセットｃと同様の方法で導出されるものでもよいし、予め決められているものでもよい。
【０００９】
仮オフセットは、旧オフセットｃ0を用いずに母集団から導出されると定義するが、この定義は新オフセットｃを導出する制約条件を定義するための定義であって、実際に導出されなければならないデータではない。仮オフセットは、旧オフセットｃ0を用いずに母集団から導出されるものであるから、実際に仮オフセットが導出されるとすれば、母集団が近傍に分布する球面の中心の位置ベクトルとなる。しかし、母集団がそれから導出される球面の一部に片寄って分布している場合には、母集団がそれぞれ含んでいる誤差がその球面の導出結果に強く影響するため、真のオフセットから遠い仮オフセットが導出される可能性が高い。例えば、図１に示すように、互いに直交する固有ベクトルがｕ1，ｕ2，ｕ3となるドーナツ状に母集団が分布し、最小の主値に対応する固有ベクトルｕ3の方向について母集団の分散が最小となる場合について考える。この場合、分布の固有ベクトルｕ3の方向については母集団の分散が小さいため、固有ベクトルｕ3の方向については、そのような母集団から導出される仮オフセットと真のオフセットの位置とが遠くなる可能性が高い。一方この場合、分布の固有ベクトルｕ1の方向については母集団の分散が大きいため、固有ベクトルｕ1の方向については、そのような母集団から導出される仮オフセットと真のオフセットの位置とが近くなる可能性が高い。
【００１０】
分布の主軸方向の分散は分布の主値（本明細書では、測定値の分散共分散行列から誤差の分散共分散行列を引いた行列の固有値に相当する。）λ1、λ2、λ3を指標として表すことができるため、この装置では、分布の主値λ1、λ2、λ3に応じて各主値に対応する方向に分布する母集団を評価する。具体的には、まず分布の主軸方向と一致する座標軸α、β、γを持つ座標系において、旧オフセットｃ0に対する新オフセットｃの位置ベクトルである補正ベクトルｆと旧オフセットｃ0に対する仮オフセットの位置ベクトルｇを定義することができる。すなわち、補正ベクトルｆと位置ベクトルｇとは分布の主軸方向の基本ベクトルの一次結合として定義することができる。これは主軸変換に相当する。そして、旧オフセットｃ0に対する仮オフセットの位置ベクトルｇの各成分ｇα、ｇβ、ｇγを分布の対応する主値ｕ1，ｕ2，ｕ3の大きさに応じて重み付けして補正ベクトルｆの各成分ｆα、ｆβ、ｆγを導出すると、分散が大きい方向について母集団の信頼性を高め、分散が小さい方向について母集団の信頼性を低くして補正ベクトルｆを導出することができる。ただし、このような補正ベクトルｆ及び位置ベクトルｇの定義も、新オフセットｃを導出する制約条件を定義するための定義であって、実際に導出する必要があるデータではない。
【００１１】
結局、このように定義される補正ベクトルｆと旧オフセットｃ0との和として新オフセットｃを求めることを制約条件として新オフセットｃを導出することにより、分散が大きな主軸方向に分布している磁気データ群ほど、オフセットを更新するために用いる母集団として重く評価し、分散が小さな主軸方向に分布している磁気データ群ほど母集団として軽く評価しながら、母集団の分布の偏りに関わらずオフセットをより確からしい値に更新することができる。
【００１２】
尚、本発明に備わる複数の手段の各機能は、構成自体で機能が特定されるハードウェア資源、プログラムにより機能が特定されるハードウェア資源、又はそれらの組み合わせにより実現される。また、これら複数の手段の各機能は、各々が物理的に互いに独立したハードウェア資源で実現されるものに限定されない。また、本発明は装置の発明として特定できるだけでなく、プログラムの発明としても、そのプログラムを記録した記録媒体の発明としても、方法の発明としても特定することができる。また、請求項に記載された方法の各動作の順序は、技術上の阻害要因がない限り、記載順に限定されるものではなく、どのような順番で実行されてもよく、また同時に実行されてもよい。
【図面の簡単な説明】
【００１３】
【図１】本発明の実施形態に係る模式図。
【図２】本発明の実施形態に係る模式図。
【図３】本発明の実施形態に係るブロック図。
【図４】本発明の実施形態に係るブロック図。
【図５】本発明の第一実施形態に係るフローチャート。
【図６】本発明の第一実施形態に係る模式図。
【図７】本発明の第一実施形態に係る模式図。
【図８】本発明の第二実施形態に係るフローチャート。
【発明を実施するための形態】
【００１４】
以下、本発明の実施の形態を以下の順に説明する。
Ａ．第一の実施形態
［１．全体説明］
１−１．ハードウェア構成
１−２．ソフトウェア構成
［２．処理の流れ］
２−１．全体の流れ
２−２．バッファ更新
２−３．分布の評価
２−４．数理計画問題による新オフセットの導出
２−５．分布が立体的である場合の新オフセットの導出
２−６．分布が平面的である場合の新オフセットの導出
２−７．分布がやや直線的である場合の新オフセットの導出
２−８．まとめ
Ｂ．第二の実施形態
・概要
・分布の評価
・新オフセットの導出
Ｃ．他の実施形態
【００１５】
［１．全体説明］
１−１．ハードウェア構成
図２は、本発明が適用される移動体の一例である携帯型電話機３の外観を示す模式図である。携帯型電話機３には３次元磁気センサ４が搭載されている。３次元磁気センサは互いに直交するｘ、ｙ、ｚの３方向について磁界の強さを検出することによって地磁気の方向及び強さを検出する。携帯型電話機３のディスプレイ２には、文字や画像の各種情報が表示される。例えば、ディスプレイ２には、地図と、方位を示す矢印や文字が表示される。
【００１６】
図３は、３次元磁気センサ４と磁気データ処理装置１とを備える磁気測定装置を示すブロック図である。３次元磁気センサ４は、地磁気による磁界ベクトルのｘ方向成分、ｙ方向成分、ｚ方向成分をそれぞれ検出するｘ軸センサ３０とｙ軸センサ３２とｚ軸センサ３４とを備えている。ｘ軸センサ３０、ｙ軸センサ３２、ｚ軸センサ３４は、いずれも磁気抵抗素子、ホール素子等で構成され、指向性のある１次元磁気センサであればどのようなものであってもよい。ｘ軸センサ３０、ｙ軸センサ３２及びｚ軸センサ３４は、それぞれの感度方向が互いに直交するように固定されている。ｘ軸センサ３０、ｙ軸センサ３２及びｚ軸センサ３４の出力は、時分割して磁気センサＩ／Ｆ（InterFace）２２に入力される。磁気センサＩ／Ｆ２２では、ｘ軸センサ３０、ｙ軸センサ３２及びｚ軸センサ３４からの入力が増幅された後にＡＤ変換される。磁気センサＩ／Ｆ２２から出力されるディジタルの磁気データは、バス５を介して磁気データ処理装置１に入力される。
【００１７】
磁気データ処理装置１は、ＣＰＵ４０とＲＯＭ４２とＲＡＭ４４とを備えている所謂コンピュータである。ＣＰＵ４０は、例えば携帯型電話機３の全体制御を司るプロセッサである。ＲＯＭ４２は、ＣＰＵ４０によって実行される磁気データ処理プログラムや、移動体の機能を実現するための種々のプログラム、例えばナビゲーションプログラムが格納されている、不揮発性の記憶媒体である。ＲＡＭ４４はＣＰＵ４０の処理対象となるデータを一時的に保持する揮発性の記憶媒体である。尚、磁気データ処理装置１と３次元磁気センサ４とを１チップの磁気測定装置として構成することもできる。
【００１８】
１−２．ソフトウェア構成
図４は、磁気データ処理プログラム９０の構成を示すブロック図である。磁気データ処理プログラム９０は、ナビゲーションプログラム９８に方位データを提供するためのプログラムであって、ＲＯＭ４２に格納されている。方位データは地磁気の方向を示す２次元ベクトルデータである。尚、方位データは例えば移動体の姿勢認識のための３次元ベクトルデータとして他のアプリケーションプログラムに提供されてもよい。磁気データ処理プログラム９０は、バッファ管理モジュール９２、オフセット導出モジュール９４、方位導出モジュール９６等のモジュール群で構成されている。
【００１９】
バッファ管理モジュール９２は、磁気センサ４から順次出力される磁気データを入力し、入力した磁気データをオフセット更新に用いるためにバッファに蓄積するプログラム部品であって、ＣＰＵ４０、ＲＡＭ４４及びＲＯＭ４２を入力手段及び蓄積手段として機能させる。尚、このバッファはハードウェアで構成されてもよいし、ソフトウェアで構成されてもよい。このバッファに蓄積される複数の磁気データを母集団というものとする。
【００２０】
オフセット導出モジュール９４は、バッファ管理モジュール９２によって保持されている母集団と、自らが保持している旧オフセットとに基づいて新オフセットを導出し、旧オフセットを新オフセットに更新するプログラム部品であって、ＣＰＵ４０、ＲＡＭ４４及びＲＯＭ４２をオフセット導出手段として機能させる。尚、旧オフセットが新オフセットに更新された時点でその新オフセットは旧オフセットになるため、誤解のない文脈では旧オフセットのことを単にオフセットというものとする。実際には、方位データの補正に用いられるオフセットは１つの変数に設定され、新オフセットはその変数とは別の変数として導出され、導出された時点で方位データの補正に用いられる変数に設定されるので、方位データの補正に用いられる変数は旧オフセットが格納される変数である。
【００２１】
方位導出モジュール９６は、磁気センサ４から順次出力される磁気データをオフセット導出モジュール９４が保持しているオフセットによって補正して方位データを生成するプログラム部品であって、ＣＰＵ４０、ＲＡＭ４４及びＲＯＭ４２を方位導出手段として機能させる。具体的には、方位導出モジュール９６は、３次元ベクトルデータである磁気データの各成分からオフセットの各成分を引き算して得られる３成分のうちの２成分又は３成分全部を方位データとして出力する。
【００２２】
ナビゲーションプログラム９８は、目的地までの誘導経路を探索し、地図上に誘導経路を表示する周知のプログラムである。地図の認識のし易さから、地図は現実の方位に地図上の方位が一致するように画面表示される。したがって例えば、携帯型電話機３が回転すると、ディスプレイ２に表示される地図は地面に対して回転しないようにディスプレイ２に対して回転する。このような地図の表示処理に方位データが用いられる。もちろん、方位データは、単に東西南北を文字や矢印で表示するためにのみ用いられてもよい。
【００２３】
［２．処理の流れ］
２−１．全体の流れ
図５は、新オフセットを導出する処理の流れを示すフローチャートである。図５に示す処理は、オフセットの更新要求が発生したときにＣＰＵ４０がオフセット導出モジュール９４を実行することによって進行する。
【００２４】
２−２．バッファ更新
ステップＳ１００では、新オフセットの導出に用いられる複数の磁気データ（母集団）を蓄積するためのバッファに蓄積されている全ての磁気データが消去される。すなわちこのとき、旧オフセットの導出に用いられた複数の磁気データ（母集団）が消去される。
ステップＳ１０２では、新オフセットの導出に用いられる磁気データが入力されバッファへ格納される。携帯型電話機３の姿勢がほとんど変化していない状況において、短い時間間隔で順次磁気センサ４から磁気データを入力すると、連続入力される２つの磁気データ間の距離が近くなる。距離が互いに近い複数の磁気データが限られた容量のバッファに格納されることは、メモリ資源の浪費であるし、無駄なバッファの更新処理を発生させる。また、互いの距離が近い磁気データ群に基づいて新オフセットを導出すると、偏った分布を持つ母集団に基づいて精度の低い新オフセットが導出される可能性がある。そこで、バッファの更新必要性が次のように判定されてもよい。例えば、直前にバッファに格納された磁気データと最後に入力された磁気データとの距離があるしきい値より小さければ、バッファの更新必要性がないと判定され、最後に入力された磁気データはバッファに格納されることなく破棄される。
【００２５】
ステップＳ１０４では、精度よく新オフセットを導出するために十分な規定個数の磁気データがバッファに蓄積されたかが判定される。すなわち、母集団の要素数は予め決められている。予め決める母集団の要素数を少なく設定することにより、オフセットの更新要求に対するレスポンスが向上する。規定個数の磁気データがバッファに蓄積されるまで、ステップＳ１０２とステップＳ１０４の処理が繰り返される。
【００２６】
２−３．分布の評価
規定個数の磁気データがバッファに蓄積されると、母集団の分布が評価される（ステップＳ１０６、ステップＳ１０８）。分布は、分布の主値に基づいて評価される。母集団としての磁気データ群を次式（１）で表し、磁気データの真値ｐｉと誤差Δｐｉを導入し、誤差Δｐｉを、確率分布μに従う確率変数ｖｉで近似する。すなわち、ｑｉ（ｉ＝１，２，・・・）は次式（１）となる。
qi=ｐｉ＋ｖｉ＝(qix,qiy,qiz)  (i=1,2,・・・)・・・（１）
分布の主値は、母集団の重心（平均）を始点とし各磁気データを終点とするベクトルを用いて式（２）、（３）、（４）、（５）で定義される対称行列Ａ'の固有値である。
ただし、確率分布μを平均が０、分散共分散行列がＶで表される正規分布とし、
【数１３】
2011203243.tif 000015
である。
行列Ａは、式（６）とも書ける。
【数１４】
2011203243.tif 000016
【００２７】
行列Ａ'の固有値を大きい順にλ1、λ2、λ3とする。またλ1、λ2、λ3に対応し、大きさ１で正規化され互いに直交する固有ベクトルをｕ1、ｕ2、ｕ3とする。２つ以上の固有値が０になるとき、言い換えると行列Ａ'の階数が１以下のときは、母集団の要素数が１であるか、分布が完全な直線であるため、考慮する必要はない。
母集団の分布は、最小固有値であるλ3と、中間固有値であるλ2とに基づいて評価される。
【００２８】
ステップＳ１０６では母集団の分布が十分立体的であるかが判定される。具体的には、次式（７）を満足する場合に肯定判定され、満足しない場合に否定判定される。
λ3＞ｔ1  かつ  λ2＞ｔ2・・・（７）
ｔ1、ｔ2は、予め決められている定数である。ｔ1、ｔ2の値をどのように設定するかは設計事項であり、オフセットの導出特性をどのように定めるかによって任意に設定することができる。式（７）が満たされる場合、母集団の重心から等方的に母集団が分布していることになる。母集団の重心から等方的に母集団が分布していることは、特定の球面近傍に満遍なく母集団が分布していることを意味する。
【００２９】
ステップＳ１０８では母集団の分布が十分平面的であるかが判定される。具体的には、次式（８）を満たす場合に肯定判定され、満足しない場合に否定判定される。
λ3≦ｔ1  かつ  λ2＞ｔ2・・・（８）
式（８）が満たされる場合、特定の平面近傍に限定した範囲内では、母集団の重心から等方的に母集団が分布していることになる。特定の平面近傍に限定した範囲内において、母集団の重心から等方的に母集団が分布していることは、特定の球面の断面円の円周近傍に偏って母集団が分布していることを意味する。
【００３０】
ステップＳ１０８で否定判定される場合、母集団の分布はやや直線的である。母集団の分布がやや直線的であることは、特定の球面の断面円の短い弧に母集団が偏って分布しているか、その断面円の直径の両端に母集団が偏って分布していることを意味する。
【００３１】
２−４．数理計画問題による新オフセットの導出
ここで新オフセットを導出するための数理計画問題について説明する。
母集団が同一平面上にない４つの磁気データで構成されている場合、母集団が分布する球面は統計的手法を用いることなく一意に特定される。この球面の中心の位置ベクトルｃ＝（ｃx、ｃy、ｃz）は連立方程式（９）を解くことによって得られる。尚、３変数に対して等号制約が４つあるが、等号制約の１つは冗長になっているため方程式（９）は必ず解を持つ。
【数１５】
2011203243.tif 000017
ただし、
【数１６】
2011203243.tif 000018
【００３２】
母集団の要素数が５個以上あるときについて、ｊを次式（１１）で定義する。
【数１７】
2011203243.tif 000019
このとき、ｃについての連立一次方程式（１２）が解を持てば、その解は、母集団が分布する球面の中心である。
Ｘｃ=ｊ・・・（１２）
【００３３】
しかし、３次元磁気センサ４の測定誤差を考慮すると、現実には、方程式（１２）が解を持つことはまれである。そこで、ｄｉのばらつきの指標である次式（１３）の関数ｆｖ（ｃ）を最小にするｃが新オフセットとして最も確からしいといえる。
【数１８】
2011203243.tif 000020
ただし、
【数１９】
2011203243.tif 000021
である。
ここで、
【数２０】
2011203243.tif 000022
とすると、
【数２１】
2011203243.tif 000023
である。
【００３４】
確率分布μを平均が０、分散共分散行列がＶで表される正規分布とし、時間について独立な確率変数ｖｉを含む項を確率分布μに対応する期待値で置換すると、式（２１）で表された確率論的モデルは、次式（２２）で表される決定論的モデルに置換される。ただし、既に述べたとおりＡ'＝Ａ−（ｎ−１）Ｖである。式（２２）は式（２１）の期待値であって、式（２１）の時間について独立な確率変数ｖｉを含む項を確率分布μに対応する期待値で置換することにより求められる。
【数２２】
2011203243.tif 000024
【００３５】
ここで、Ｅ［ ］を期待値を意味する記号とすると、式（２１）から式（２２）は以下の式（２３）によって表現される。
【数２３】
2011203243.tif 000025
ここで、確率変数ｖｉを含むｆ（ｃ）の各項は、
【数２４】
2011203243.tif 000026
と置換されている。
【００３６】
２−５．分布が立体的である場合の新オフセットの導出
分布が立体的である場合、旧オフセットに対する補正方向を制限することなく新オフセットが導出される（ステップＳ１１４）。分布が立体的である場合、すなわち、母集団の重心からみてある程度全方位的に母集団が分布している場合、母集団は全方向について十分信頼できる。したがって、このような場合には、新オフセットを導出するために必ずしも旧オフセットを用いる必要が無く、旧オフセットを用いることなく母集団に基づいて新オフセットを導出することもできる。旧オフセットを用いることなく母集団に基づいて新オフセットを導出するアルゴリズムは、種々提案されており、統計的手法を用いるアルゴリズムであってよいし、本件出願人によって既に出願されている特開２００５−３３７４１２号公報及び特開２００６−４４２８９号公報に開示されているように統計的手法を用いないアルゴリズムであってもよい。
【００３７】
本実施形態では、確率論的手法により新オフセットが導出される。すなわち、ステップＳ１１４では、式（２２）の目的関数を制約条件無しで最小にする数理計画問題の解として新オフセットｃが導出される。
Ａ'が正定値のとき、式（２２）で表される目的関数ｆ（ｃ）を最小にするｃは、
【数２５】
2011203243.tif 000027
である。Ａ'が正定値でないとき、式（２２）で表される目的関数ｆ（ｃ）を最小にするｃは存在しないか、あるいは一意に定まらない。
【００３８】
ところで、
【数２６】
2011203243.tif 000028
とすると、式（２２）で表される目的関数ｆ（ｃ）を最小にするｃは、
【数２７】
2011203243.tif 000029
である。
すなわち、母集団の分布が十分立体的である場合、ステップＳ１１４では、式（２７）、（３０）のいずれかを満たすｃを新オフセットとして導出すればよい。式（２７）、（３０）のいずれを用いて新オフセットを導出するかについては設計事項である。
【００３９】
２−６．分布が平面的である場合の新オフセットの導出
図６に示すように母集団の分布が平面的である場合、旧オフセットに対する補正方向を互いに直交する２方向に制限して新オフセットが導出される（ステップＳ１１２）。母集団がある特定の平面近傍に分布し、その平面の法線方向から見て離散的に分布している場合、その平面と平行な方向については母集団を十分信頼できる一方でその平面の法線方向については母集団を信頼できないことになる。このような場合には、その平面の法線方向については旧オフセットを補正しないことにより、信頼に値しない情報に基づいてオフセットが更新されることを防止できる。
【００４０】
母集団がある特定の平面近傍に分布し、その平面の法線方向から見て離散的に分布している場合、その平面の法線方向は、最小固有値λ3に対応する固有ベクトルｕ3の方向に一致し、その平面に平行で互いに直交する方向は、最大固有値λ1、中間固有値λ2にそれぞれ対応するｕ1、ｕ2の方向に一致する。したがって、その平面の法線方向について旧オフセットを補正せずに新オフセットを導出するため、次式（３１）の制約条件のもとで式（２２）の目的関数が最小になる新オフセットｃを求める。
ｃ＝ｃ0＋β1ｕ1＋β2ｕ2（β1、β2は実数）・・・（３１）
式（３１）は次式（３２）と等価である。
【数２８】
2011203243.tif 000030
【００４１】
式（３１）の制約条件のもとで式（２２）の数理計画問題を解く式は、ラグランジュの未定乗数法によって等価な連立方程式に変形できる。未定乗数ρを導入し、
【数２９】
2011203243.tif 000031
とすると、ｘの連立一次方程式（３４）がその方程式となる。
【数３０】
2011203243.tif 000032
ただし、
【数３１】
2011203243.tif 000033
【００４２】
すなわち、母集団の分布が平面的である場合、ステップＳ１１２において新オフセットを導出するには、連立一次方程式（３４）を解けばよい。尚、行列Ｂ4の階数（rank）は必ず４になるため、解ｘは必ず一意に決定される。また最終的に導出するのはｃであるが、これも当然一意に決定される。
なお、
【数３２】
2011203243.tif 000034
とし、ｃ'、未定乗数ρ'についての連立方程式（３８）の解を用いて、
【数３３】
2011203243.tif 000035
として新オフセットｃを導出しても良い。
【００４３】
２−７．分布がやや直線的である場合の新オフセットの導出
図７に示すように母集団の分布がやや直線的である場合、旧オフセットに対する補正方向を分布の主方向に制限して新オフセットが導出される（ステップＳ１１０）。母集団がある特定の直線近傍に分布し、その直線方向には離散的に分布している場合、その直線の方向については母集団を十分信頼でき、それ以外の方向については母集団を信頼できないことになる。このような場合には、その直線の方向以外の方向については旧オフセットを補正しないことにより、信頼に値しない情報に基づいてオフセットが更新されることを防止できる。
【００４４】
母集団がある特定の直線近傍に分布し、その直線方向には離散的に分布している場合、その直線の方向は、最大固有値λ1に対応する固有ベクトルｕ1の方向に一致し、それ以外の方向のベクトルは中間固有値λ2、最小固有値λ3にそれぞれ対応する固有ベクトルｕ2、ｕ3の一次結合となる。したがって、その直線の方向についてのみ旧オフセットを補正して新オフセットを導出するため、次式（４０）の制約条件のもとで式（２２）の目的関数が最小になる新オフセットｃを求める。
【数３４】
2011203243.tif 000036
式（４０）は式（４１）と式（４２）とを同時に満足することと等価である。
【数３５】
2011203243.tif 000037
【００４５】
式（４０）の制約条件のもとで式（２２）の数理計画問題を解く式は、ラグランジュの未定乗数法によって等価な連立方程式に変形できる。未定乗数ρ1、ρ2を導入し、
【数３６】
2011203243.tif 000038
とすると、ｘの連立一次方程式（４４）がその方程式となる。
【数３７】
2011203243.tif 000039
ただし、
【数３８】
2011203243.tif 000040
とする。
【００４６】
すなわち、母集団の分布がやや直線的である場合、ステップＳ１１０において新オフセットを導出するには、連立一次方程式（４４）を解けばよい。尚、行列Ｂ5の階数は必ず５になるため、解ｘは必ず一意に決定される。また最終的に導出するのはｃであるが、これも当然一意に決定される。
なお、
【数３９】
2011203243.tif 000041
とし、ｃ'、未定乗数ρ１'、ρ２'についての連立方程式（４８）の解を用いて、
【数４０】
2011203243.tif 000042
として新オフセットｃを導出しても良い。
【００４７】
２−８．まとめ
ここで図１、６、７に基づいて、ステップＳ１１０、Ｓ１１２、Ｓ１１４の処理を空間的な概念を用いて説明する。母集団を完全に信頼できるものとすれば、旧オフセットｃ0に対する、母集団のみから導出される球面の中心の位置ベクトルｇと、旧オフセットｃ0との和として新オフセットｃを考えることにより、新オフセットｃは次式（５０）によって定義される。
ｃ＝ｃ0＋ｇ・・・（５０）
【００４８】
式（２１）の目的関数を制約条件無しで最小にする数理計画問題の解として得られる位置ベクトルｇは、分布の固有ベクトルｕ1、ｕ2、ｕ3と同じ方向の基本ベクトルの一次結合である。このため、位置ベクトルｇの係数ｇα、ｇβ、ｇγを、対応する主軸方向についての母集団の信頼度に応じて重み付けすることにより、位置ベクトルｇを各成分の信頼度に応じて補正したベクトルに相当する補正ベクトルｆを得ることができる（図１参照）。
【００４９】
図６に示すように母集団の分布が平面的である場合に実行されるステップＳ１１２では、旧オフセットｃ0と母集団とに基づいて新オフセットを導出するにあたり次の制約条件を付している。その制約条件は、分布の最大の主値、すなわち最大固有値λ1に対応する分布の主軸方向の成分ｇαと、分布の中間の主値、すなわち中間固有値λ2に対応する主軸方向の成分ｇβについて、位置ベクトルｇの重み係数をともに１とし、分布の最小の主値、すなわち最小固有値λ3に対応する主軸方向の成分ｇγについて、位置ベクトルｇの重み係数を０として求まる補正ベクトルｆと旧オフセットｃ0の和として新オフセットｃを求めることである。
【００５０】
図７に示すように母集団の分布がやや直線的である場合に実行されるステップＳ１１０では、旧オフセットｃ0と母集団とに基づいて新オフセットを導出するにあたり次の制約条件を付している。その制約条件は、分布の最大の主値、すなわち最大固有値λ1に対応する分布の主軸方向（主方向）の成分ｇαについて、位置ベクトルｇの重み係数を１とし、分布の中間の主値、すなわち中間固有値λ2に対応する主軸方向の成分ｇβと、分布の最小の主値、すなわち最小固有値λ3に対応する主軸方向の成分ｇγについて、位置ベクトルｇの重み係数をともに０として求まる補正ベクトルｆと旧オフセットｃ0の和として新オフセットｃを求めることである。
【００５１】
母集団の分布が立体的である場合に実行されるステップＳ１１４では、特に制約条件が付されていない。すなわちステップＳ１１４では、式（２２）の目的関数を制約条件無しで最小にする数理計画問題の解として得られる位置ベクトルｇと旧オフセットｃ0との和として新オフセットｃが導出されるが、結果的に旧オフセットｃ０を用いずに新オフセットが導出される。
【００５２】
Ｂ．第二の実施形態
・概要
第一の実施形態では、母集団の分布を離散的に評価し、分布が平面的である場合には主値が最小になる主軸方向の成分について補正ベクトルｆを０とし、分布が直線的である場合には主値が中間値と最小値になる２つの主軸方向の成分について補正ベクトルｆを０として新オフセットｃを導出した。第二の実施形態では、第一の実施形態のように分布の評価に応じて異なる処理を実行する必要がなく、また母集団を効率よく用いてより確からしい新オフセットを導出できる、簡素かつ高精度なアルゴリズムが説明される。
【００５３】
図８は、新オフセットを導出する処理の流れを示すフローチャートである。図８に示す処理は、第一実施形態と同様に、オフセットの更新要求が発生したときにＣＰＵ４０がオフセット導出モジュール９４を実行することによって進行する。尚、ステップＳ２００の処理は第一実施形態において説明したステップＳ１００の処理と同一である。ステップＳ２０２の処理は第一実施形態において説明したステップＳ１０２の処理と同一である。ステップＳ２０４の処理は第一実施形態において説明したステップＳ１０４の処理と同一である。
【００５４】
・分布の評価
ステップＳ２０６では、母集団の分布指標が導出される。例えば、次式（５１）、（５２）によって定義するｍ2、ｍ3を分布指標として導出することにより、母集団の分布が評価される。
【数４１】
2011203243.tif 000043
【００５５】
ｋ2、ｋ3は予め決められている正の定数である。ｋ2、ｋ3の値によって主値とそれに対応する母集団の主軸方向の信頼度の対応関係が定まる。ここで、ｍ2、ｍ3は次の条件を全て満たす。また次の条件を全て満たす限り、母集団の分布指標としてのｍ2、ｍ3をどのように定めても良い。
λ２＜０のとき、ｍ２＝１
λ３＜０のとき、ｍ３＝１
ｍ２，ｍ３の値域が［０，１］となるか、その部分集合となること
ｍ２がλ２の増加に対して広義の単調減少をすること
ｍ３がλ３の増加に対して広義の単調減少をすること
【００５６】
対応する主値が大きい分布の主軸方向成分の順に位置ベクトルｇの係数をｇα、ｇβ、ｇγとし、対応する主値が大きい主軸方向成分の順に補正ベクトルｆの係数をｆα、ｆβ、ｆγとするとき、位置ベクトルｇと補正ベクトルｆとｍ2、ｍ3との関係は次式（５３）、（５４）、（５５）の通りである。
【数４２】
2011203243.tif 000044
【００５７】
・新オフセットの導出
数理計画問題の解を特定の制約条件の下で導出することが困難である場合、その制約条件を緩めて数理計画問題を解くことに相当する緩和問題が導入されることがある。本実施形態では、このような緩和問題を応用することにより、母集団の分布の主値に対応する重み係数によって前述した位置ベクトルｇの係数ｇα、ｇβ、ｇγ（図１参照）を重み付けして得られる補正ベクトルｆと、旧オフセットｃ0との和として新オフセットｃを導出する処理が実現される。具体的には次の通りである。
【００５８】
途中計算に必要な変数として未定乗数ρ1、ρ2を導入し、ｃ、ρ1、ρ2を１つのベクトルにまとめ、次式（５６）によりｘを定義する。
【数４３】
2011203243.tif 000045
【００５９】
また、行列Ｂを式（５７）、ベクトルｂを式（５８）によって定義する。
【数４４】
2011203243.tif 000046
【００６０】
ステップＳ２０８において新オフセットを導出する処理は、次の連立方程式（５９）の解を求める処理となる。尚、行列Ｂは必ず正則であるので、ｘは一意に決定される。したがって当然にｃも一意に決定される。
Ｂｘ＝ｂ・・・（５９）
なお、
【数４５】
2011203243.tif 000047
としてｃ'、ρ'１、ρ'２についての連立方程式（６１）の解を用いて、
【数４６】
2011203243.tif 000048
として新オフセットｃを導出しても良い。
連立方程式（５９）の解を求めることは、母集団の分布の主値に対応している重み係数ｆα／ｇα、ｆβ／ｇβ、ｆγ／ｇγによって、各主値に対応する分布の主軸方向の位置ベクトルｇの係数を重み付けした値を各成分とする補正ベクトルｆと旧オフセットｃ0との和として新オフセットを求めることを制約条件として、式（２２）の目的関数を最小にする数理計画問題を解くことと等価である。
【００６１】
このように、第二実施形態では、新オフセットを導出する処理を母集団の分布に応じて分岐させる必要がないため、オフセット導出モジュール９４の開発や改良が容易であり、オフセット導出モジュール９４のデータサイズが小さくなる。また、第二実施形態では、母集団の分布の主値が０以下でない限り、分布の主軸方向に旧オフセットを主値に対応する距離だけ補正することができるため、オフセット導出モジュール９４による母集団の利用効率が高く、方位導出モジュールにおいて最も確からしいオフセットを用いて磁気データを補正できる。
【００６２】
Ｃ．他の実施形態
本発明は、上記実施形態に限定されるものではなく、その要旨を逸脱しない範囲で種々の形態で実施可能である。例えば、本発明がＰＤＡに搭載される磁気センサや車両に搭載される磁気センサにも適用できることは当然である。また本発明は３次元磁気センサのみならず、２次元磁気センサにも当然に適用できる。
また確率論的モデルである確率計画問題としての目的関数を決定論的モデルである目的関数に期待値を用いて置換できる確率分布μは、正規分布に限らず、二項分布、一様分布、ポアソン分布などであってもよい。これらの確率分布のうち磁気データの誤差の分布にもっともよく当てはまるものに対応する期待値によって、確率変数である誤差を含む目的関数の項を置換すればよい。さらに、本願発明は正規分布でなくても正規分布とみなし、その分布に対応する期待値を用いてオフセットを導出するようにしてもよい。
また本発明は、磁気データに含まれる誤差Δｐｉを確率分布μに従う確率変数ｖｉとみなし、確率分布μに対応する期待値を用いてオフセットを導出する限りにおいて、いかなる確率論的手法によりオフセットを導出してもよい。例えば国際公開２００６−００９２４７、特開２００７−２２５３２５、特開２００７−２２５３２９、特開２００７−２２５３３２等の公報に開示された統計的手法に誤差の期待値を導入することによって、本発明を実施することができる。
【符号の説明】
【００６３】
１：磁気データ処理装置、３：携帯型電話機（移動体）、４：３次元磁気センサ、９０：磁気データ処理プログラム、９２：バッファ管理モジュール（入力手段、蓄積手段）、９４：オフセット導出モジュール（オフセット導出手段）、ｃ0：旧オフセット、ｃ：新オフセット、ｇ：位置ベクトル、ｆ：補正ベクトル
【図１】
2011203243.tif 000056
【図２】
2011203243.tif 000057
【図３】
2011203243.tif 000058
【図４】
2011203243.tif 000059
【図５】
2011203243.tif 000060
【図６】
2011203243.tif 000061
【図７】
2011203243.tif 000062
【図８】
2011203243.tif 000063