JP6477214B2 - 傾斜度測定方法及び装置並びに電子機器及びプログラム - Google Patents

Description

本発明は、傾斜度測定方法及び装置並びに電子機器及びプログラム等に関する。

加速度センサーの出力には、環境変化及び経時変化等によるバイアスのドリフトがある。また、加速度センサーに特有の誤差として、加速度センサーを含む機器が水平面に対して傾斜としている場合には、加速度センサーの出力には重力加速度成分に起因した誤差（オフセット）が生ずる。

特許文献１では、自転車などの移動体に配置した加速度センサー１０を水平方向に１８０°反転させる前後の静止状態（加速度＝０）で加速度検出を行ない、両方の検出値の平均値を取ることにより、傾斜に基づくオフセット量を求めている。

特許文献２は、ターンテーブルの角度位置が０度（Ｓ１）、９０度（Ｓ２）、１８０度（Ｓ３）、２７０度（Ｓ４）での感度軸の出力値に基づいて、感度軸の出力値を補正する較正値を求めている。具体的には、９０度ずつの位相の違う４点（Ｓ１〜Ｓ４）の静止データの統計平均値から２直線を表わす連立一次方程式を解き、多軸センサーの最大・最小出力点と振り幅の中心値を求めることで較正値を算出している。

特許文献３は、検出面が傾いた場合でも角速度を検出する角速度センサーを提供するために、角速度検出部と同じ面内に少なくとも三つの加速度検出部を設置している。これらの加速度検出部で検出した加速度を用いて角速度検出面の傾斜角を算出した後、傾斜角による角速度の変動量を推測して補正している。

特開２００４―３５４２１４号公報（００１３） 特開２０１０―２８１５９８号公報（００４４−００４５） 特開平０６―３２４０６６号公報

特許文献１，２では、加速度センサーを０度と１８０度の各停止位置に回転し（特許文献１）、あるいは加速度センサーを０度、９０度、１８０度、２７０度の各停止位置に回転し（特許文献２）、各停止位置で静止状態にて加速度センサーで加速度を検出し、それらの平均値に基づいてオフセット量や較正値を求めている。このような手法で求めたオフセット量や較正値でもなお誤差が存在している。

加速度センサーを搭載するセンサー機器は水平に設置されるが、例えば地震等の外的要因により設置面が傾斜することがある。あるいは、特許文献１のようにセンサー機器が移動体に搭載される場合、移動面は傾斜している。このような場合、水平面に対する傾斜度（傾斜角や傾斜方位）を測定する必要がある。傾斜度を測定する際、加速度センサーの測定値には傾斜による重力加速度成分が誤差（傾斜誤差とも称する）として重畳される。また、特許文献１，２のようにセンサー機器を回転してオフセット量や較正値を求める場合でも、回転前後の各停止位置にて異なるドリフト（動的誤差とも称する）や、各停止位置にて実質的に一定のバイアス（静的誤差とも称する）を除去する必要がある。

本発明の幾つかの態様は、加速度センサーの出力に重畳する誤差として、回転前後の各停止位置にて異なるドリフト（動的誤差）や、各停止位置にて実質的に一定のバイアス（静的誤差）を除去して、加速度センサーの設置面の傾斜度（傾斜角や傾斜方位）を精度高く測定できる傾斜度測定方法及び装置並びに電子機器及びプログラムを提供することを目的とする。

（１）本発明の一態様は、重力方向と交差する仮想平面に検出軸を有する加速度センサーの前記検出軸を、基準方位に設定した第１測定方位と、前記仮想平面に対する法線の周りの回転位置にて前記第１測定方位とは異なる少なくとも２つの第２測定方位と、前記基準方位に戻した第３測定方位と、の各々に設定した静止状態にて、３つ以上の異なる測定方位を含む計４つ以上の測定方位で測定される前記加速度センサーの出力を用いて、前記第１〜第３測定方位の各々で異なる動的誤差を測定する誤差測定工程と、
前記第１測定方位及び前記少なくとも２つの第２測定方位の各々の測定方位に設定された前記加速度センサーの出力から、前記動的誤差と、前記第１〜第３測定方位の各々で実質的に一定の静的誤差とが除去された、前記仮想平面が水平面に対して傾斜する傾斜角及び傾斜方位角に依存するデータに基づいて、前記傾斜角及び前記傾斜方位角の少なくとも一方を測定する工程と、
を有し、
前記誤差測定工程は、前記第１測定方位及び前記第３測定方位に設定された前記加速度センサーの出力の差分に基づいて、前記動的誤差を測定する工程を含む傾斜度測定方法に関する。

第１測定方位と第３測定方位とは、加速度センサーの検出軸が共に基準方位に設定される同一方位と平行になる。ここで、静止状態で測定される加速度センサーの出力には、その測定方位に固有の動的誤差（ドリフト）と、他の測定方位でも等しい静的誤差（バイアス）と、その測定方位に固有の傾斜誤差が含まれる（後述の式（７）（８）を参照のこと）。また、後述の式（７）（８）に示すように、傾斜誤差は傾斜角φと傾斜方位θに依存するので、未知数は４である。よって、検出軸を計３つ以上の異なる測定方位に設定して計４回以上測定（基準方位で２回測定）している。一軸加速度センサーであれば、一軸検出軸が４つ以上の異なる測定方位に設定され、二軸加速度センサーであれば、交差二軸のうちの各一軸を２つ以上の測定方位に設定すれば、計４つ以上の測定方位に設定したことになる。

ここで、平行な同一方位である第１，第３測定方位での加速度センサーの出力の差分は、第１測定方位と第３測定方位との動的誤差の差分となる。この動的誤差の差分は、測定中の傾斜変化あるいは、温度特性変動、環境振動変化、ランダムドリフトなどが考えられるが、測定自体を比較的短時間に行われ、測定タイミングも環境温度変動が小さい、安定している時間帯を選択し、また測定時に測定端末の消費電流変化も安定させる条件とすることで、測定中の傾斜変化あるいは、環境振動変化、などを殆ど無視できるほどに縮小できる。また測定条件を適宜選択することで、ランダムドリフトを殆ど無視できるほどに縮小できる。そこで、例えば第１測定方位の動的誤差は零とし、第３測定方位では動的誤差が最大または最小となり、第３測定方位と第１測定方位との間の少なくとも２つの第２測定方位の動的誤差は、第１，第３測定方位間の動的誤差の差分に基づいて割り当てることができる。

ランダムドリフトを無視すると、静止状態での加速度センサーの出力は、水平面に対する仮想平面の傾斜に起因した重力加速度成分による傾斜誤差（傾斜角φと傾斜方位θに依存する）、動的誤差及び静的誤差の総和となる。静止状態での加速度センサーの出力から動的誤差及び静的誤差を除去することで、傾斜誤差が求められ、傾斜誤差から傾斜角及び前記傾斜方位角の少なくとも一方を測定することができる。

（２）本発明の一態様では、前記第１〜第３測定方位の各々での前記動的誤差を、前記差分を線形補完して求めることができる。動的誤差が一定のレート（傾き）で変化するものとみなすことで、各測定方位での動的誤差は線形補完により求めることができる。

（３）本発明の一態様では、前記誤差測定工程は、前記第１測定方位と前記少なくとも２つの第２測定方位との各々の測定方位で測定された前記加速度センサーの出力の各々から、前記各々の測定方位について測定された前記動的誤差がそれぞれ減算された４つ以上の一次加工データを生成する工程と、前記４つ以上の一次加工データに基づいて、前記第１〜第３測定方位の各々で実質的に一定の静的誤差を測定する工程と、をさらに含むことができる

ここで、加速度センサーの出力にはランダムドリフトが重畳するが、計測条件を最適化することでランダムドリフトの悪影響を少なくでき、本発明の一態様ではランダムドリフトを無視することができる。そうすると、静止状態での加速度センサーの出力は、傾斜誤差、動的誤差及び静的誤差の総和となる。各測定方位にて静止状態で測定された加速度センサーの出力から動的誤差を減算して一次加工データを生成した上で、静的誤差がそれぞれ等しいことを利用して、４つ以上の一次加工データから傾斜誤差を消去することで、静的誤差を求めることができる。

（４）本発明の一態様では、前記誤差測定工程は、前記４つ以上の一次加工データを平均化して、前記静的誤差を求めることができる。４つ以上の一次加工データを平均化することで、傾斜誤差を消去して静的誤差を求めることができる。

（５）本発明の一態様では、前記加速度センサーの前記検出軸を回転角９０°または１８０°で回転させることできる。傾斜誤差は回転角に依存する正弦波または余弦波の成分を含むが、回転角９０°または１８０°とすることで計算が単純化される。また直交二軸加速度センサーの二軸を異なる測定方位として用いることができる。

（６）本発明の一態様では、前記第１〜第３測定方位での測定を１サイクルとしたとき、サイクル数ｎを２以上とする場合、ｎサイクルに亘って測定し、ｎサイクル目の前記第１測定方位を（ｎ−１）サイクル目の前記第３測定方位として兼用することができる。

このようにｎサイクルに亘って測定値をサンプリングすることで、測定精度を高めることができる上、連続サイクルの境界となる測定方位を兼用することで、隣り合う２サイクルにて同一測定値を利用することができる。このようにＮサイクルに亘って測定値をサンプリングすることで、測定精度を高めることができる上、前後のサイクルで測定方位を兼用することで、前後のサイクルにて同一測定値を利用することができる。

また、本発明の一態様では、本測定は連続して方位を変えながら、静止状態で計測を行うもので、測定は、測定順にほぼ等時間間隔で測定され、回転して方位を変えながら連続して計測した時の測定データ列から1サイクルの区間を抽出し、解析することができる。この抽出した1サイクルの区間は、サイクルの最初と最後の方位が同一であり、その他の方位は、測定に必要なすべての方位を含む区間である。この測定順は、方位の並び順に優先し、測定順に方位が並ぶ必要はない。

（７）本発明の一態様では、前記加速度センサーの前記検出軸は、前記第１測定方位から前記少なくとも２つの第２測定方位に向けて第１方向に回転された後に、前記少なくとも２つの第２測定方位から前記第３測定方位に向けて、前記第１方向とは逆方向の第２方向に回転することができる。こうすると、動的誤差の偏りが平均化される。

本発明の一態様では、抽出した１サイクルの測定順に動的誤差を測定し、測定結果から、動的誤差を補正したデータ列について、１サイクルのデータ中の同一方位のデータについて平均値を求めて、測定方位毎に丸めたデータを生成し、この結果について静的誤差を測定する手順に進むことができる。

本発明の一態様では、回転計測の測定方位の基準方向を、センサーの設置された環境のNorth-East-Down（NED）座標系に関連付けて方位を決定して設置することで、回転計測の結果得られる傾斜方位角を環境の座標系方位で示すことができる。

（８）本発明の他の態様は、
重力方向と交差する仮想平面に検出軸を有する加速度センサーの前記検出軸を、基準方位である第１測定方位と、前記仮想平面に対する法線の周りの回転位置にて前記第１測定方位とは異なる少なくとも２つの第２測定方位と、前記基準方位に戻した第３測定方位と、の各々に設定して、前記検出軸を３つ以上の異なる測定方位を含む計４つ以上の測定方位に設定する指令を出力する手段と、
前記第１〜第３測定方位にそれぞれ設定された前記加速度センサーの出力を受信する手段と、
前記加速度センサーの出力を用いて、前記加速度センサーの出力に重畳する、前記第１〜第３測定方位の各々で異なる動的誤差を測定する誤差測定手段と、
前記前記第１測定方位及び前記少なくとも２つの第２測定方位の各々の測定方位に設定された前記加速度センサーの出力から、前記動的誤差と、前記第１〜第３測定方位の各々で実質的に一定の静的誤差とが除去された、前記仮想平面が水平面に対して傾斜する傾斜角及び傾斜方位角に依存するデータに基づいて、前記傾斜角及び前記傾斜方位角の少なくとも一方を測定する手段と、
を有し、
前記誤差測定手段は、前記第１測定方位及び前記第３測定方位に設定された前記加速度センサーの出力の差分に基づいて、前記動的誤差を測定する傾斜度測定装置に関する。

本発明の他の態様に係る傾斜度測定装置により、本発明の一態様に係る傾斜度測定方法を好適に実施することができる。

（９）本発明のさらに他の態様は、（８）に記載の誤差測定装置を有する電子機器を定義している。

（１０）本発明のさらに他の態様は、
重力方向と交差する仮想平面に検出軸を有する加速度センサーの前記検出軸を、基準方位である第１測定方位に設定し、前記仮想平面に対する法線の周りの回転位置にて前記検出軸を前記第１測定方位とは異なる少なくとも２つの第２測定方位に設定し、前記仮想平面に対する法線の周りの回転位置にて前記検出軸を前記基準方位である第３測定方位に戻す指令を出力する手順と、
前記第１〜第３測定方位にそれぞれ設定された前記加速度センサーからの出力を受信する手順と、
前記第１〜第３測定方位の各々の静止状態で、前記検出軸を３つ以上の異なる測定方位を含む計４つ以上の測定方位に設定して測定される前記加速度センサーの出力を用いて、前記第１〜第３測定方位の各々で異なる動的誤差を測定する手順と、
前記第１測定方位及び前記少なくとも２つの第２測定方位の各々の測定方位に設定された前記加速度センサーの出力から、前記動的誤差と、前記第１〜第３測定方位の各々で実質的に一定の静的誤差とが除去された、前記仮想平面が水平面に対して傾斜する傾斜角及び傾斜方位角に依存するデータに基づいて、前記傾斜角及び前記傾斜方位角の少なくとも一方を測定する手順と、
をコンピューターに実施させ、
前記動的誤差を測定する手順は、前記第１測定方位及び前記第３測定方位での前記加速度センサーの出力の差分に基づいて、前記動的誤差を測定する手順を含むプログラムに関する。

このプログラムをパーソナルコンピューター等の汎用機にインストールすることで、汎用機にて誤差測定を実現できる。

本発明の一実施形態に係る測定システムの全体図である。 図１に示す測定端末装置の外観図である。 図２に示す筐体内に配置される三軸加速度センサーの配置を示す斜視図である。 水平座標系に対して傾斜したセンサー座標系を示す図である。 加速度センサーを回転させた時の回転角度と、各静止位置でのセンサー出力、動的誤差、静的誤差及び傾斜誤差との関係を示す図である。 図１に示す端末装置及び測定端末装置のブロック図である。 ２軸４方位回転測定にてＸ軸加速度センサーで得られた測定値を正規化した値（縦軸）と回転角度との関係を示す特性図である。 ２軸４方位回転測定にて直交二軸（Ｘ軸及びＹ軸）加速度センサーで得られた測定値を正規化した値（縦軸）と回転角度との関係を示す特性図である。 重力加速度に基づく入力Ｉφの傾斜方位角を示す説明である。 ２軸４方位回転測定にて直交二軸（Ｘ軸及びＹ軸）加速度センサーで得られた測定値を示す図である。 ２軸４方位回転測定にて直交二軸（Ｘ軸及びＹ軸）加速度センサーで得られた出力から動的誤差を減算した一次加工データを示す図である。 ２軸４方位回転測定にて直交二軸（Ｘ軸及びＹ軸）加速度センサーで得られた出力から動的誤差及び静的誤差を減算した二次加工データ（傾斜誤差）を示す図である。 図１２に示すデータを傾斜量Ｉφで除算した、規格化された三次加工データを示す図である。 傾斜量Ｉφで規格化された三次加工データを複素平面上に置き換えて示す図である。 傾斜方位θを実数成分と虚数成分とで示すフェザー座標を説明する図である。 ２軸４方位回転測定で計算される三次加工データと傾斜方位角θとの関係を示す図である。 ２軸２方位回転測定で収集される加速度の規格化された値と方位との関係を示す図である。 ２軸２方位回転測定で収集される加速度の測定値を示す図である。 ２軸２方位回転測定で計算される一次加工データを示す図である。 ２軸２方位回転測定で計算される二次加工データを示す図である。 ２軸２方位回転測定で計算される二次加工データの計算例を示す図である。 ２軸２方位回転測定で計算される三次加工データを示す図である。 ２軸２方位回転測定で計算されるノルムＩφ、ｓｉｎθ、ｃｏｓθ、ｔａｎθを示す図である。 ２軸２方位回転測定で計算される三次加工データと傾斜方位角θとの関係を示す図である。 １軸３方位４回回転測定で収集される加速度の測定値を示す図である。 １軸３方位４回回転測定で計算される一次加工データを示す図である。 傾斜量Ｉφと傾斜方位角θとの関係を示す図である。 傾斜方位角θが存在する象限、ｍ及び係数の関係を示す図である。 θ＝±π／２の場合にＩφｃｏｓθ≒０の値をプラス側またはマイナス側に設定する場合の解を示す図である。 North-East-Down（NED）座標系を示す図である。 加速度センサーの出力ノイズ分析を示す図である。 傾斜方位角θの判別を説明するための図である。 回転区間ｔ_Ｓ、測定区間ｔ_Ｒ及び一サイクル区間を説明するための図である。 各サイクルにて最初に設定される方向を同一とした複数サイクルを示す図である。 各サイクルにて最初に設定される方向を異ならせた複数サイクルを示す図である。 図３５の変形例を示す図である。 傾斜方位角θφ、測定方位と傾斜方位角の差ψｎ、測定方位角θｎ及び仮想平面Ｐ_φを説明するための図である。 傾斜角φと傾斜量Ｉφとを説明するための図である。 一次加工データを示す図である。 二次加工データを示す図である。 方位角方向の求め方を示す図である。 静的誤差の温度依存性を示す図である。 更正値の更新動作を示すフローチャートである。 傾斜測定動作を示すフローチャートである。 各サイクルで基準方位を変えながらデータ収集を複数サイクルで実施する例を示す図である。 各サイクルで基準方位を変えながらデータ収集を複数サイクルで実施する他の例を示す図である。 測定開始指令から待ち時間を経てデータ収集する動作フローを示す図である。

以下、本発明の好適な実施の形態について、比較例を参照して詳細に説明する。なお以下に説明する本実施形態は特許請求の範囲に記載された本発明の内容を不当に限定するものではなく、本実施形態で説明される構成の全てが本発明の解決手段として必須であるとは限らない。

１． 測定システム
図１に示す測定システムは、測定データを集中管理するホスト装置１と、ホスト装置１と例えばインターネットや無線、携帯電話で接続された複数の端末装置２と、各端末装置２と例えばWiFi（登録商標）やBluetooth（登録商標）等の無線通信、移動体通信システム、ローカルエリアネットワーク（ＬＡＮ）で接続された一又は複数の測定端末装置３と、を有する。端末装置２は携帯電話やスマートフォン、タブレット、又は端末ＰＣ等の汎用機であっても専用機であっても良い。測定端末装置３は、加速度センサーに接続或いは内蔵している機器、携帯電話やスマートフォン、タブレット、又は端末ＰＣ等の汎用機であっても専用機、加速度計測器であってもよい。

各測定端末３では、設置時での初期較正と、定期的な維持較正が実施される。各測定端末３には、経時的変化や環境変化に起因したドリフトや、設置面の傾斜に起因した誤差が生ずるからである。その較正を実施するために、ホスト装置１または各端末装置２の記憶装置に較正プログラム（誤差測定プログラム、傾斜度測定プログラム、バイアス補正プログラム等を含む）が格納されている。各測定端末３での較正は、較正プログラムに従って実行される。較正プログラムは、ホスト装置１より端末装置２に送信されても良い。

図２に示す測定端末装置２は、底面４Ａを設置面とする筐体４内部に、図３に示す三軸（Ｘ，Ｙ，Ｚ）加速度センサー１０ｘ，１０ｙ，１０ｚを有する。なお、誤差測定や傾斜角・傾斜方位角測定には加速度センサー１０ｚの出力は用いず、後述の通り二軸加速度センサー１０ｘ，１０ｙの一方又は双方を用いることができる。換言すれば、加速度センサー１０ｚの出力を用いずに測定できる、傾斜角φが過度に大きくない範囲で本実施形態は有効である。

三軸加速度センサー１０ｘ，１０ｙ，１０ｚの座標Ｘ，Ｙ，Ｚにおいて、直交Ｘ，Ｙ軸を含むＸ−Ｙ平面（仮想平面ともいう）は底面４Ａと平行であり、Ｚ軸はＸ−Ｙ平面と直交する。筐体４には回転操作部５とインジケーター（ＬＥＤ等）６が設けられている。測定端末２にて、較正プログラムに従った処理が実施される手順を、インジケーター６の点灯または点滅により操作者に指示することができる。指示に従い回転操作部５が回転操作されると、底面（設置面）４Ａと直交する法線（Ｚ軸）回りで、三軸加速度センサー１０ｘ，１０ｙ，１０ｚが回転される。なお、特に定期的な維持較正を無人で実施するために、回転操作部５に代えて、あるいは回転操作部５に加えて、上述の軸を自動回転させるモーター等の回転駆動部８（図６参照）を設けても良い。具体的には、筐体４の外からツマミにより加速度センサーを回転することができ、ツマミが必要な方向で止まるようにクリックがある。また、回転による較正作業の開始をプログラムに通知するトリガ入力部がある。回転させて良い時間を表示したり、較正作業が終了したことを示す表示部がある。また取得した較正値の適用を指示する入力部がある。

２． 加速度センサーの出力に重畳する誤差
図４には、水平座標系（Ｘｒ，Ｙｒ，Ｚｒ）と、センサー座標系（Ｘ，Ｙ，Ｚ）が示されている。水平座標系（Ｘｒ，Ｙｒ，Ｚｒ）の直交二軸Ｘｒ，Ｙｒを含む（Ｘｒ−Ｙｒ）平面は水平面であり、Ｚｒ軸は水平面（Ｘｒ−Ｙｒ）に直交する鉛直軸である。

図４に示すように、水平座標系（Ｘｒ，Ｙｒ，Ｚｒ）のＺｒ軸（鉛直軸）とセンサー座標系（Ｘ，Ｙ，Ｚ）のＺ軸とは傾斜角φであり、センサー座標系（Ｘ，Ｙ，Ｚ）の（Ｘ−Ｙ）平面は水平座標系（Ｘｒ，Ｙｒ，Ｚｒ）の（Ｘｒ−Ｙｒ）平面に対して傾斜しているものと仮定する。回転操作部５により三軸加速度センサー１０ｘ，１０ｙ，１０ｚをＺｒ軸周りで回転させると、Ｚｒ軸に沿って作用する重力加速度Ｇは、図４に示すようにＺ軸を中心に円又は楕円或いは直線Ｃを描く。この円又は楕円或いは直線Ｃの大きさは、回転静止される二軸加速度センサー１０ｘ，１０ｙの測定値Ｄｘ，Ｄｙ値のいずれか一方又は双方から求めることができる。また、（Ｘ−Ｙ）平面と（Ｘｒ−Ｙｒ）平面の交わる直線がＸｒまたはＹｒのどちらか一方と平行の場合、Ｚｒ軸に沿って作用する重力加速度Ｇは、Ｚ軸を中心に直線の軌跡を描く。この直線の大きさは、測定値Ｄｘ，Ｄｙ値のいずれか一方又は双方から求めることができる。なお、加速度センサー１０ｘ，１０ｙ，１０ｚの座標系（Ｘ，Ｙ，Ｚ）では、加速度Ｘと加速度Ｙは右手系座標となっており、−Ｚ方向から見てＹはＸに対して反時計回りに９０度の方向となっている。

加速度センサー１０ｘ，１０ｙ，１０ｚの出力軸を、図３０に示すNorth-East-Down（NED）座標系に方位が関連付けて位置決定された第１測定方位を基準方位とすることができる。図３０に示すN-E平面で示される水平面に投影したセンサーの基準方向を、NED座標系のN-E平面内の方位に関連付けする。加速度センサー１０ｘ，１０ｙ，１０ｚによって傾斜方位角ψを求めるため、この傾斜方位角ψは加速度センサー１０ｘ，１０ｙ，１０ｚが設置された場所で、設置場所のNED座標系と加速度センサー１０ｘ，１０ｙ，１０ｚの出力軸方向が関連付けされていることで、傾斜の方角を示すことができる。つまり、図３０に示すように、加速度センサーの基準方向出力の水平面への投影方向とNED座標系の基準とする方向のなす角をψとする。

ここで、三軸加速度センサー１０ｘ，１０ｙ，１０ｚが静止している場合、入力Iは重力加速度Ｇの傾き成分となる。傾斜角φによる加速度ＧのＸ−Ｙ平面上のバイアスシフトのノルムＩφは、
となる。また、図４の通りＸ軸とＩφとの交差角を、傾斜方位角θとする。加速度センサー１０ｘの測定値Ｄｘと、加速度センサー１０ｙの測定値Ｄｙとは、ノルムＩφと傾斜方位角θを用いて次の式で表すことができる。

加速度センサー１０ｘ，１０ｙの測定値には、実際にはドリフトによる誤差が重畳されている。このドリフトを次のように定義する。

ここで、三軸加速度センサー１０ｘ，１０ｙが静止している場合、入力Iは重力加速度Ｇの傾き成分Ｉφのみとなるので、式（５）は次式（６）となる。
Ｘ，Ｙ軸のバイアスＤＦをＤ_Ｆｘ，Ｄ_Ｆｙとする。ランダムドリフトレートＤ_Ｒは、レート・ホワイト・ノイズＤ_ＲＮ、フリッカー・ノイズＤ_ＲＢ、レート・ランダム・ウォーク・ノイズＤ_ＲＲの総和で示されるが、計測条件を最適化することでこれらの悪影響を少なくすることで、以下の測定では無視することができる。

図３１に示す加速度センサーの出力ノイズ分析から、フリッカー・ノイズＤ_ＲＢが主たる区間Ａ−Ｂをデータ取得時の平均区間として選択する。レート・ホワイト・ノイズＤ_ＲＮは、最小になる測定区間の時間長さ以上に調整します。測定データの平均計算区間をこのように選択することで、レート・ホワイト・ノイズＤ_ＲＮの少ない精度の高い測定データが得られる。レート・ランダム・ウォーク・ノイズＤ_ＲＲのように平均区間が長くなるとドリフトが増大するノイズは、データ測定環境下で最小になる測定区間の時間長さ以下に調整する。測定データの平均計算区間をこのように選択することで、レート・ランダム・ウォーク・ノイズＤ_ＲＲが少ない精度の高い測定データが得られる。フリッカー・ノイズＤ_ＲＢは除去できないので、計測時間が不要に長くならないように適正化します。従って、データの平均値区間と異なり、方位を変えて測定する１サイクル全体がこの適正区間の最長以内に入るように選択することが望ましい。計測条件はランダムドリフトレートＤ_Ｒが最小となるＡ-Ｂ区間の計測時間とするのが望ましい。

環境によるドリフトレートＥは、測定中の温度ドリフトによるオフセットＤｘ_ΔＴΔＴ，Ｄｙ_ΔＴΔＴと表記する。以上のことから、加速度センサー１０ｘの測定値Ｄｘと加速度センサー１０ｙの測定値Ｄｙは、誤差を考慮すると次の式（７）（８）の通りとなる。

式（７）のモデルは、加速度センサー１０ｘを回転させると、回転角度（横軸）と出力（縦軸）の関係は図５の通りとなる。図５に示す加速度センサー１０ｘの測定値Ｄｘの振幅は、傾斜角φに依存する。式（８）のモデルも同様となる。ここで、ＤＦｘ，ＤＦｙは回転前後の各静止位置で実質的に一定であることから、静的誤差とも称する。一方、温度ドリフトによるオフセットＤｘ_ΔＴΔＴ，Ｄｙ_ΔＴΔＴは、回転前後の各静止位置で異なることから、動的誤差とも称する。なお図５では、温度ドリフトＤｘ_ΔＴΔＴを一定レートとして示している。実際の温度ドリフトＤｘ_ΔＴΔＴは一定レートではないと予想されるが、ここでは一定のレートに近似して示している。また、重力加速度Ｇの傾き成分Ｉφ・ｃｏｓθ，Ｉφ・ｓｉｎθは、傾斜による誤差としてセンサー出力に重畳することから、傾斜誤差と称する。

３． 端末装置及び測定端末装置
図６に、端末装置２及び測定端末装置３を示す。測定端末装置３は、三軸加速度センサー１０ｘ，１０ｙ，１０ｚ、回転操作部５、インジケーター６、通信装置７を有する。端末装置２は上述した回転駆動部８を有していても良い。また、測定端末装置３は、回転操作部５と回転駆動部８を外部に有してもよい。

端末ＰＣ等で構成される端末装置２は、上述した較正プログラム（誤差測定プログラム、傾斜度測定プログラム、バイアス補正プログラム等）を記憶する記憶装置２０を有する。端末装置２で実施される誤差測定、傾斜度測定、バイアス補正は、較正プログラムに従って実施される。

例えば端末装置２は、較正プログラムに従って実行制御される回転司令部３０、誤差測定部４０、傾斜度測定部５０、傾斜バイアス補正部６０及び通信装置７０等を有する。回転司令部３０は、較正プログラムに従って、通信装置７０を介して測定端末装置３に送信される回転指令を生成する。

誤差測定部４０は、測定端末装置３からの出力に基づいて、加速度センサー１０ｘ，１０Ｙの出力に重畳している誤差を測定する誤差測定部４０は、ドリフト（動的誤差）を測定する動的誤差測定部４１と、バイアス（静的誤差）を測定する静的誤差測定部４２と、例えば加速度センサー１０Ｘ，１０Ｙの出力から動的誤差及び静的誤差を除去して傾斜誤差を測定する傾斜誤差測定部４３と、を有する。なお、加速度センサー１０Ｘ，１０Ｙの出力を処理する過程で静的誤差を除去できる場合には、静的誤差測定部４２は不要であるか、あるいは不使用とすることができる。

傾斜度測定部５０は、例えば加速度センサー１０Ｘ，１０Ｙの出力から動的誤差及び静的誤差を除去して得られたデータから、傾斜度（傾斜角、傾斜方位）を測定する。傾斜度測定部５０は、傾斜角測定部５１と傾斜方位角測定部５２とを有する。

これら回転司令部３０、誤差測定部４０（動的誤差測定部４１、静的誤差測定部４２、傾斜誤差測定部４３）、傾斜度測定部５０（傾斜角測定部５１、傾斜方位角測定部５２）の機能は、以下にて説明する各種の回転測定方法にて明らかにされる。

傾斜バイアス補正部６０は、傾斜誤差測定部４３で求められる傾斜誤差を傾斜バイアスとして、加速度検出時に、加速度センサーからの出力から傾斜バイアスを取り除くものである。
４． 誤差測定方法
以下に、較正プログラムに従って回転操作部５を操作して設定される回転位置（静止位置）での動的誤差Ｄｘ_ΔＴΔＴ，Ｄｙ_ΔＴΔＴ、静的誤差ＤＦｘ，ＤＦｙ及び傾斜誤差Ｉφ・ｃｏｓθ，Ｉφ・ｓｉｎθの測定方法について説明する。

本実施形態では、測定手法として以下の通り定める。
（ｉ）測定時に設定される加速度センサーの検出軸の方位を測定方位と称する。少なくとも一つの測定方位（基準方位）を重複させて測定し、基準方位とは異なる方位でも測定する。好ましくは最初と最後の測定方位を基準方位とする。

（ii）式（７）（８）中の未知数はそれぞれ４つ（例えば式（７）では、Ｉφ、ｃｏｓθ、ＤＦｘ及びＤｘ_ΔＴΔＴの４つ）であるので、加速度センサーの検出軸を少なくとも４つ以上の測定方位（静止位置）に設定して測定する。つまり、未知数は、傾斜誤差Ｉφｃｏｓθ、バイアス（静的誤差）ＤＦｘ及びドリフト（動的誤差）Ｄｘ_ΔＴΔＴであり、これらの誤差を求めることになる。ただし、一検出軸での測定を１つと数え、二軸加速度センサー１０ｘ，１０ｙの出力を用いれば一静止位置での測定は２つと数えることができる。

（iii）測定間の回転角度は、測定方位の数に依存して任意に定めることができるが、計算の都合上好ましくは測定間の回転角度は９０度または１８０度刻みとする。
測定間の回転角度を９０度または１８０度刻みとする場合には、２軸２方位回転測定（測定方位は０度、１８０度、３６０度）、２軸４方位回転測定（測定方位は０度、９０度、１８０度、２７０度、３６０度）、１軸３方位４回回転測定（（測定方位は０度、９０度、１８０度、３６０度）、１軸４方位回転測定（測定方位は０度、９０度、１８０度、２７０度、３６０度）等の各種測定方法を挙げることができる。

４．１． ２軸４方位回転測定
図７は、２軸４方位で測定した時の測定方位（０度、９０度、１８０度、２７０度、３６０度）毎の加速度センサー１０ｘの測定値を示す。図７に示すように、０度と３６０度で測定値が同じである場合、この測定の間に加速度センサー１０ｘの出力中のドリフトが安定していると仮定する。測定中のドリフトが安定していると仮定した場合、０度方位から＋２７０度方位の間で直交する４方位の測定値は、測定時点の測定系の傾斜量と、加速度センサー１０ｘ（または１０ｙ）のドリフト（静的誤差）とによるオフセットの和であるとする。０度と３６０度で値に変動があるようなら、その差分を動的誤差として各測定に等分に割り当てて引き算する必要がある（動的誤差の除去）。任意の検出軸の直交する４方位での測定値の平均値は、バイアスのドリフトによるオフセット量を示す。任意の検出軸の直交する４方位の出力の回転方位に対する軌跡は、図７の通りサインカーブ上の４点となり、その振幅は水平面に対する傾きとなる。また、図７に示す最大値と最小値を取る方位は、その傾きの方位となる。

なお、図７の縦軸は、測定値そのものではなく、予想される最大振幅と測定結果の比（振幅±１）に正規化されている。もし、センサー座標系（Ｘ，Ｙ，Ｚ）の（Ｘ−Ｙ）平面が水平面に対して平行である場合、図７の縦軸の値は全て０になる。図７と同様にして、直交二軸加速度センサー１０ｘ，１０ｙの出力を示すと、図８の通りとなる。図８に示すように、直交二軸加速度センサー１０ｘ，１０ｙの出力の位相は９０度ずれている。

ここで、センサー座標系の（Ｘ−Ｙ）平面を図９の通り定義する。Ｘ軸方向に対してＹ軸は反時計回りに９０度の方向とする。加速度センサーの回転方向は、例えば反時計回りとし、０度→９０度→１８０度→２７０度→０度（３６０度）の方向に回転させる。リリースする時には、内部配線が伸びきらないように、時計回り方向に戻すことができる。

図８に示すように、０度方向からの傾斜方位角をθとし、傾斜量をＩとし、測定順をｎ（ｎ＝０，１，２，３）とし、測定方位をｎ・π／２とすると、式（７）（８）に示す各測定値は次の通りとなる。
例えば、ｎ＝０の時の０度方向の測定値は、誤差をＸｅ，Ｙｅと表記すると次の式（９）（１０）の通りとなる。

４．１．１． 測定手順１（２軸４方位データ収集）
先ず、図２に示す回転操作部５を操作して直交二軸加速度センサー１０ｘ，１０ｙを回転させ、例えば０度→９０度→１８０度→２７０度→０度（３６０度）の順に各測定方位で停止させる。各測定方位にて、静止状態で２軸４方位データを加速度センサー１０ｘ，１０ｙにて収集する。２軸４方位データを図１０の通り定義する。

４．１．２． 測定手順２（動的誤差の測定）
次に、動的誤差である、測定中のドリフトＤｘ_ΔＴΔＴ，Ｄｙ_ΔＴΔＴを測定する。図７について上述した通り、もし最初（０度）と最後（３６０度＝０度）の測定値が等しい場合には測定中のドリフトは無視できるが、変動があるようなら、最初（０度）と最後（３６０度）の測定値の差分からドリフトＤｘ_ΔＴΔＴ，Ｄｙ_ΔＴΔＴを求めることができる。なぜなら、式（７）（８）から分かる通り、測定方位が０度で共に等しい測定値Ｘ０，Ｘ４は、ドリフトＤｘ_ΔＴΔＴのみが異なり、他の値であるＩφ・ｃｏｓθ及びＤＦｘについては測定値Ｘ０，Ｘ４間で共に等しいからである。本実施形態にて、最初と最後の測定方位を等しくした理由は、全ての測定方位でのドリフト（動的誤差）の測定を可能とするためである。

ここで、測定方位が０度で共に等しい測定値Ｘ０，Ｘ４の差分、測定値Ｘ０，Ｘ４中の動的誤差の差分となる。この動的誤差の差分は、測定中の傾斜変化或いは、温度特性変動、環境振動変化、などが考えられるが、測定自体を比較的短時間に行われ、測定タイミングも環境温度変動が小さい、安定している時間帯を選択し、また測定時に測定端末の消費電流変化も安定させる条件とすることで、殆ど無視できるほどに縮小できる。そこで、最初の０度の静止位置の動的誤差は零とし、最後の３６０度の３静止位置では動的誤差が最大または最小となり、その間の９０度、１８０度、２７０度の各静止位置での動的誤差は、測定値Ｘ０，Ｘ４の差分に基づいて割り当てることができる。

例えば、上述したように測定中のドリフトＤｘ_ΔＴΔＴ，Ｄｙ_ΔＴΔＴをそれぞれ一定レート（傾き）で近似した場合、最初（０度）と最後（３６０度）の測定値の差分（Ｘ４−Ｘ０）を線形補完することで、ｎ回目の測定値Ｘｎに重畳しているドリフト（動的誤差）は、ｎ（Ｘ４−Ｘ０）／４と表すことができる。こうして、全ての測定方位での測定値に重畳しているドリフト（動的誤差）を測定することが可能となる。なお、ｎ番目の測定値Ｘｎからドリフト（動的誤差）が除外されたデータＸｎ’は次の通り示される。また、図１０に示す各測定値からはドリフト（動的誤差）が除外された一次加工データＸｎ’Ｙｎ’は図１１に示すようになる。

４．１．３． 測定手順３（静的誤差の測定）
式（９）（１０）に示す表記を用いると、測定値Ｘ０〜Ｘ３は、次のように示すことができる。

Ｘ０〜Ｘ３の計測値の平均は、
となる。右辺第1項のカッコ内は、測定の前提としての対称性から、
となり、このことから、下記の式が成立する。
以上のことから、次の式が成立する。
このように、直交する測定方位（０度、９０度、１８０度、２７０度）の計測値を加算平均することで、正弦波または余弦波の成分の値は消去されて、誤差Ｘｅのみを取り出すことができる。

誤差Ｘｅ＝ＤＦｘ＋ＤＲｘ＋Ｄｘ_ΔＴΔＴまたは誤差Ｘｅ＝ＤＦｙ＋ＤＲｙ＋Ｄｙ_ΔＴΔＴであるので、Ｘavgに代えて、Ｘ’avgを計算する。
データＸ’は動的誤差Ｄｘ_ΔＴΔＴ，Ｄｙ_ΔＴΔＴが除去されているため、その平均値Ｘ’ａｖｇも動的誤差Ｄｘ_ΔＴΔＴ，Ｄｙ_ΔＴΔＴを含まない誤差（Ｘｅ−Ｄｘ_ΔＴΔＴ）または（Ｘｅ−Ｄｙ_ΔＴΔＴ）を求めることができる。ここで、上述したようにランダムドリフトＤＲｘ，ＤＲｙを無視できるとすれば、平均値Ｘ’ａｖｇ＝ＤＦｘまたは平均値Ｙ’ａｖｇ＝ＤＦｙとなり、平均値Ｘ’ａｖｇまたはＹ’ａｖｇにより静的誤差ＤＦｘまたはＤＦｙが求まることが分かる。

４．１．４． 測定手順４（傾斜誤差の測定）
次に、動的誤差が除去されたデータＸ’からさらに静的誤差ＤＦｘを除去するために次式を計算し、図１２に示す二次加工データを求める。
こうして求められた二次加工データＸ”ｎは、測定値Ｘから動的誤差Ｄｘ_ΔＴΔＴ及び静的誤差ＤＦｘが減算されたものであり、式（７）からＸ”ｎ＝Ｉφ・ｃｏｓθが傾斜誤差として求められる。同様に、式（８）からＹ”ｎ＝Ｉφ・ｓｉｎθが傾斜誤差として求められる。

以上により、２軸４方位データから、動的誤差、静的誤差及び傾斜誤差の全てを求めることができる。

４．１．５． 測定手順５（傾斜角φの測定）
求められた二次加工データＸ”ｎを測定順序ｎを用いて示すと次式の通りである。
この式を用いると、手順４で求めた各二次加工データＸ”ｎの関係は次の通りとなる。
測定誤差等を縮小するために、次の通りの組み合わせで平均を計算する。
このＸ”_０２とＸ”_１３は互いに直交するので、この２値から傾斜量Ｉφを次式の通りに計算する。
この傾斜量Ｉφは傾斜角φの水平面への投影であることから、傾斜角Ｉ_φangleは、次の通りとなる。
測定値がＧ単位である場合、傾斜角φは次の通り表記される。
以上の通り、２軸４方位データから傾斜角φを求めることができる。

４．１．６． 測定手順６（傾斜方位角θの測定）
先ず、手順５で求めた傾斜量Ｉφで、手順４で求められた各二次加工データＸ”ｎを次式の通り除算して、図１３に示す三次加工データを求める。
この三次加工データは、Ｘ'''＝ｃｏｓθ、Ｙ'''＝ｓｉｎθとなり、図８に示す振幅±１に規格化された正弦波または余弦波を意味する。

この正弦波または余弦波を、図１４に示すようにガウス平面に直交座標に並べ替え、その際に測定方向を考慮して極性を付与する。
式（１２）を、式（１１）を用いて変形するとともに、図１４に示すＸ，Ｙの実数毎及びＸ，Ｙの虚数毎に平均する。それにより、Ｘ軸及びＹ軸の複素数表現は次式の通りとなる。

式（１８）に示すＸ軸及びＹ軸の複素数のノルムは、次式の通り１となる。
例えばＸ軸の実数及び虚数とノルムとの関係は図１５に示す通りとなり、図１５から傾斜方位角θについては、実数及び虚数から次式が成立する。

図１５に示す第１〜第４象限の各区間で取り得る値から、次のことが成立する。
第１象限 θｓ＞０，θｃ＞０ 結果θ＝θｃ，θｓ＝θｃ
第２象限 θｓ＞０，θｃ＞０ 結果θ＝θｃ，θｓ≠θｃ
第３象限 θｓ＜０，θｃ＞０ 結果θ＝２π−θｃ，θｓ＝θｃ
第４象限 θｓ＜０，θｃ＞０ 結果θ＝２π−θｃ，θｓ≠θｃ
その結果、θｓ＞０，θｃ＞０の場合には、次式が成立する（単位はradian）。
同様に、θｓ＞０，θｃ＞０の場合には、次式が成立する（単位はradian）。
以上のことから、傾斜方位角θｘｎ，θｙｎは図１６の通りとなる。図１６に示す傾斜方位角θｘｎ，θｙｎの平均値を計算することもできる。
Ｘ平均値θ_{ａｖｇ−ｘ}とＹ平均値θ_{ａｖｇ−ｙ}は前提としては同値になるが、種々の要因で異なる値となる場合があるが概ね等しい。
以上の通り、測定手順１〜６により、式（７）（８）中の未知数４つの全てについて解を求めることができる。

４．２． ２軸２方位回転測定
４．２．１． 測定手順１（２軸２方位データ収集）
先ず、図２に示す回転操作部５を操作して直交二軸加速度センサー１０ｘ，１０ｙを回転させ、０度→１８０度→０度（３６０度）の順に各測定方位で停止させる。各測定方位にて、静止状態で２軸２方位での測定値を加速度センサー１０ｘ，１０ｙにて収集する。図１７は、２軸２方位で測定した時の測定方位（０度、１８０度、３６０度）毎の加速度センサー１０ｘの測定値を示す（Ｙ軸は省略）。図７に示す２軸４方位回転測定と同様に、０度と３６０度で値が同じである場合、この測定の間に加速度センサー１０ｘの出力中のドリフトが安定していると仮定する。もし変動があるようなら、その差分を動的誤差として各測定に等分に割り当てて引き算する。なお、図１７の縦軸も図７と同様に測定値そのものではなく、予想される最大振幅と測定結果の比（振幅±１）に正規化されている。２軸２方位データを図１８の通り定義する。

４．２．２． 測定手順２（動的誤差の測定）
次に、式（７）（８）中の動的誤差である、測定中のドリフトＤｘ_ΔＴΔＴｎ，Ｄｙ_ΔＴΔＴｎを測定する。図１７について上述した通り、もし最初（０度）と最後（３６０度＝０度）の測定値が等しい場合には測定中のドリフトは無視できるが、変動があるようなら、２軸４方位回転測定の測定手順２（４．１．１）と同様にして、最初（０度）と最後（３６０度）の測定値の差分からドリフトＤｘ_ΔＴΔＴｎ，Ｄｙ_ΔＴΔＴｎを求めることができる。

ここで、測定方位が０度で共に等しい測定値Ｘ０，Ｘ２の差分が、測定値Ｘ０，Ｘ２中の動的誤差の差分となる。最初の０度の静止位置の動的誤差は零とし、最後の３６０度の３静止位置では動的誤差が最大または最小となり、その間の１８０度の静止位置での動的誤差は、測定値Ｘ０，Ｘ２の差分に基づいて割り当てることができる。例えば、上述したように測定中のドリフトＤｘ_ΔＴΔＴ，Ｄｙ_ΔＴΔＴをそれぞれ一定レート（傾き）で近似した場合、最初（０度）と最後（３６０度）の測定値の差分（Ｘ２−Ｘ０）を線形補完することで、ｎ回目の測定値Ｘｎに重畳しているドリフト（動的誤差）は、ｎ（Ｘ２−Ｘ０）／２と表すことができる。こうして、全ての測定方位での測定値に重畳しているドリフト（動的誤差）を測定することが可能となる。なお、ｎ番目の測定値Ｘｎからドリフト（動的誤差）が除外されたデータＸｎ’は次の通り示される。また、図１８に示す各測定値からはドリフト（動的誤差）が除外された一次加工データＸｎ’，Ｙｎ’は図１９に示すようになる。

４．２．３． 測定手順３（静的誤差の測定）
２軸４方位回転測定の処理手順３（４．１．３）と同様にして、式（１３）で定義されるＸ’０とＸ’１との平均から静的誤差ＤＦｘを求め、Ｙ’０とＹ’１との平均から静的誤差ＤＦｙを求める。

４．２．４． 測定手順４（傾斜誤差の測定）
次に、動的誤差が除去された一次加工データＸ’からさらに静的誤差ＤＦｘを除去するために次式を計算し、図２０に示す二次加工データを求める。
こうして求められた二次加工データＸ”ｎは、測定値Ｘから動的誤差Ｄｘ_ΔＴΔＴ及び静的誤差ＤＦｘが減算されたものであり、式（７）からＸ”ｎ＝Ｉφ・ｃｏｓθが傾斜誤差として求められる。同様に、式（８）からＹ”ｎ＝Ｉφ・ｓｉｎθが傾斜誤差として求められる。式（１６）を計算することにより、各静止位置での傾斜誤差Ｘ”ｎ，Ｙ”ｎは図２１の通りとなる。
以上により、２軸２方位回転測定から、動的誤差、静的誤差及び傾斜誤差の全てを求めることができる。

４．２．５． 測定手順５（傾斜角φの測定）
処理手順４で求められた二次加工データＸ”ｎ，Ｙ”ｎを次式の通り変形して、Ｘ”０，Ｘ”１，Ｙ”０，Ｙ”１をそれぞれ求める。

以上のことから、次式が成立する。
測定誤差等を縮小するため、Ｘ軸、Ｙ軸についてそれぞれ次式により平均化する。なお、次式では正負の符号を考慮して平均化している。
上述の式を計算して次式を求める。

このＸ”０１とＹ”０１は直交しているので、この２値から次式により傾斜量Ｉφを計算することができる。
この傾斜量Ｉφは傾斜角φの水平面への投影であることから、２軸４方位回転測定の式（１２）（１３）と同様にして傾斜角φを求めることができる。
以上の通り、２軸２方位データから傾斜角φを求めることができる。

４．２．６ 測定手順６（傾斜方位角θの測定）
先ず、手順５で求めた傾斜量Ｉφで、手順４で求められた各二次加工データＸ”ｎを式（１４）の通り除算して規格化して、図２２に示す三次加工データＸ'''ｎを求める。
次に、０度及び１８０度の三次加工データについて、Ｘ軸を実数としＹ軸を虚数として複素数Ｖを表すと次の通りとなる。

０度及び１８０度の三次加工データについての複素数Ｖを用いて、図１５に示す複素数のノルム、ｓｉｎθ、ｃｏｓθ及びｔａｎθを求めると、図２３の通りとなる。また、傾斜方位角θｘｎ，θｙｎは図２４の通りとなる。図２４に示す傾斜方位角θｘｎ，θｙｎの平均値を計算することもできる。
Ｘ平均値θavg-xとＹ平均値θavg-yは前提としては同値になるが、種々の要因で異なる値となる場合があるが概ね等しい。
以上の通り、測定手順１〜６により、２軸２方位回転測定でも式（７）（８）中の未知数４つの全てについて解を求めることができる。なお、２軸４方位回転測定と比べると、２軸２方位回転測定ではデータ数が小さくなるので精度は劣化するが、測定量が約１／２となるので計測上の負荷は減少する。

４．３． 一軸３方位４回回転測定
４．３．１． 測定手順１（２軸２方位データ収集）
先ず、図２に示す回転操作部５を操作して直交二軸加速度センサー１０ｘ，１０ｙを回転させ、０度→９０度→１８０度→０度（３６０度）の順に各測定方位で停止させる。各測定方位にて、静止状態で１軸３方位４回データを加速度センサー１０ｘ，１０ｙのいずれか一方にて収集する。図２５は、１軸３方位で４回測定した時の測定方位（０度、９０度、１８０度、３６０度）毎の加速度センサー１０ｘ，１０ｙの測定値を示す。以降の説明では、Ｘ軸加速度センサー１０の出力を用いて説明する。図７及び図１７に示す回転測定と同様に、０度と３６０度で値が同じである場合、この測定の間に加速度センサー１０ｘの出力中のドリフトが安定していると仮定する。もし変動があるようなら、その差分を動的誤差として各測定に等分に割り当てて引き算する。

４．３．２． 測定手順２（動的誤差の測定）
次に、式（７）（８）中の動的誤差である、測定中のドリフトＤｘ_ΔＴΔＴｎ，Ｄｙ_ΔＴΔＴｎを測定する。図２５について上述した通り、もし最初（０度）と最後（３６０度＝０度）の測定値が等しい場合には測定中のドリフトは無視できるが、変動があるようなら、２軸４方位回転測定の測定手順２（４．１．２）と同様にして、最初（０度）と最後（３６０度）の測定値の差分からドリフトＤｘ_ΔＴΔＴｎ，Ｄｙ_ΔＴΔＴｎを求めることができる。

ここで、測定方位が０度で共に等しい測定値Ｘ０，Ｘ３の差分が、測定値Ｘ０，Ｘ３中の動的誤差の差分となる。最初の０度の静止位置の動的誤差は零とし、最後の３６０度の３静止位置では動的誤差が最大または最小となり、その間の１８０度の静止位置での動的誤差は、測定値Ｘ０，Ｘ３の差分に基づいて割り当てることができる。例えば、上述したように測定中のドリフトＤｘ_ΔＴΔＴｎ，Ｄｙ_ΔＴΔＴｎをそれぞれ一定レート（傾き）で近似した場合、最初（０度）と最後（３６０度）の測定値の差分（Ｘ３−Ｘ０）を線形補完することで、ｎ回目の測定値Ｘｎに重畳しているドリフト（動的誤差）は、ｎ（Ｘ３−Ｘ０）／３と表すことができる。こうして、全ての測定方位での測定値に重畳しているドリフト（動的誤差）を測定することが可能となる。なお、ｎ番目の測定値Ｘｎからドリフト（動的誤差）が除外されたデータＸｎ’は次の通り示される。また、図１８に示す各測定値からはドリフト（動的誤差）が除外された一次加工データＸｎ’，Ｙｎ’は図１９に示すようになる。

４．３．３． 測定手順３（傾斜角及び傾斜方位角の測定）
測定中のドリフト（動的誤差）Ｄｘ_ΔＴΔＴを除いた、回転によって変動しない項を含むモデルを、次の通り定義する。なお、次式中のＸｎは測定値でなく、図２６中の一次加工データＸ’ｎ、つまりデータＸ’０（ｎ＝０）、データＸ’１（ｎ＝２）、またはデータＸ’２（ｎ＝２）を意味する（次式中では一次加工データＸ’ｎのダッシュは省略し、一軸３方位４回回転測定では以降同様とする。）。また、次式中のＸｅは、静的誤差ＤＦｘを意味する。

式（１８）の未知数は３（動的誤差は測定済み）であるから、以降の計算では次式の通りＸ０，Ｘ１，Ｘ２だけを使用する。

一次加工データＸ０，Ｘ１，Ｘ２中のＸｅが共に等しいことを利用して、解を次式の通りも止める。つまり、本実施形態では静的誤差Ｘｅを求めることはせず、静的誤差Ｘｅを式中から消去している。

以上の式より、下記の式が成立する。

式（２０）の関係を図２７に示す。式（２０）から、傾斜量Ｉφ及び傾斜方位角θは次式の通り求められる。
この傾斜量Ｉφは傾斜角φの水平面への投影であることから、２軸４方位回転測定の式（１２）（１３）と同様にして傾斜角φを求めることができる。

ただし、傾斜方位角θに関して、上述の結果は極性に注意が必要である。ここで、式（１９）中のＸｅは同値であるので、次のように定義されるＸｅ１，Ｘｅ２，Ｘｅ３が同値であるか否かを検証する。

ｍ＝０，１として次式を計算し、Ｘｅ１ｍ＝Ｘｅ２ｍ＝Ｘｅ３ｍとなるｍを選択する。
選択されたｍから、算出したθに加算する係数を図２８の通り決めて、θfinal＝θ＋ｍπを決定する。但し、θfinal＜０またはｍ＝０の場合には、反時計回りの角度表現は、θfinal ＝θ＋２πとする。

ここで、θ＝±π／２の場合には、Ｉφｃｏｓθ＝（Ｘ０−Ｘ２）＝０となり、次式によりθを求めることができない。
Ｘ０＝Ｘ２の場合には、Ｉφｃｏｓθ＝−０．０００００００００１などのマイナス側の適当な小さい数を代入し、Ｘｅ１ｍ≒Ｘｅ２ｍ≒Ｘｅ３ｍとなるｍ＝０または１を選択する。こうすると、Ｘ０＝Ｘ２の場合にはｍ＝１となり、図２９の通りθfinal＝θ＋πと一致する。

式（２０）から、θｃ、θｓ、θｔ、Ｉφを計算する。

このθｃ、θｓ、θｔの値域により、図３２を参照して傾斜方位角θの判別を行う。

なお、θｃ＝０の時、例えば次のような小さい値を代入する。

４．３．４． 測定手順４（傾斜誤差の測定）

式（１８）中の傾斜誤差ｃｏｓ（θ−ｎ・π／２）は、既知である測定値Ｘｎ（センサー出力から動的誤差が除去された測定値）と静的誤差Ｘｅとを式（１８）中に代入することで求められる。あるいは、傾斜誤差ｃｏｓ（θ−ｎ・π／２）にｎとθfinalとを代入して求めても良い。

以上の通り、測定手順１〜４により、１軸３方位４回回転測定でも式（７）（８）中の未知数４つのうち静的誤差を除く３つについて解を求めることができる。本実施形態では静的誤差は求める必要がなく、式中から消去することで除去されている。ただし、必要があれば、算出された傾斜角φや傾斜方位角θを式（１８）に代入して静的誤差Ｘｅを算出しても良い。なお、以上の説明ではＸ軸加速度センサー１０ｘの出力を用いたが、Ｙ軸加速度センサー１０ｙの出力を用いてもよい。あるいは、Ｘ軸加速度センサー１０ｘとＹ軸加速度センサー１０ｙの出力からそれぞれ解を求めても良い。

４．４． 一軸４方位回転測定
本実施形態は、４．１．二軸４方位回転測定を、Ｘ軸加速度センサー１０ｘ及びＹ軸加速度センサー１０ｙのいずれか一方を用いて実施する測定である。４．１．二軸４方位回転測定の欄で説明したように、式（７）または式（８）の４つの未知数は、Ｘ軸加速度センサー１０ｘ及びＹ軸加速度センサー１０ｙのいずれか一方を用いて算出することができる。

４．５． 回転角が９０度，１８０度以外の方位回転測定
上述した各実施形態では加速度センサーの回転角度を９０度または１８０度刻みとしたが、これに限定されない。以下、９０度，１８０度以外の回転角を含む方位回転測定について説明する。

４．５．１． 処理手順１（データ収集）
加速度センサーの検出軸を回転させ方向を変えて静止状態で測定する計測を連続して行う。図３３に示すように、加速度センサーの検出軸を回転させて所望の方向に向ける区間ｔ_Ｓと、静止状態の間に加速度を測定する区間ｔ_Ｒとの一対で、一方向の測定を構成する。この一方向の測定が繰り返され、検出軸を回転させて方向を変えながら連続して計測する。測定方位を変えながら測定された連続した測定データから、誤差を測定するデータ処理の区間を切り出し、これを1サイクルとする。計測方向は、円周（２π）を３以上の自然数で割った角度ステップで計測する。円の中心から円に内接する正多角形の頂点への方向を測定方位とする。1サイクルの区間をｍ＝０，１，２，…，ｍ_ｍとする。1サイクル区間は、区間の最初の測定方向（基準方位）ｍ＝０とし、この方向から回転したその他の方向の測定と、最初の測定方向（基準方位）と同じ方向の測定ｍ＝ｍ_ｍまでの区間とする。この1サイクルの区間内に設定される方位は、最初の方位（基準方位）と、基準方位と同じ方向である最後の方位と、その他の方位とは、誤差を測定するために必要なすべての方位を含む。測定に必要な全ての方位が1サイクル中に1回以上含まれる。従って、1サイクルの区間の長さは異なる場合がある。方位を変更する回転の方向は交互に反転してもよい。また、隣り合うサイクル区間の回転は、サイクル区間内の方向順番が対称であってもよいし、あるいは同順であってもよい。

図３４〜図３６に示すように、誤差測定に必要な方向数を例えば５つとし、その方向をａ，ｂ，ｃ，ｄ，ｅ，ｆとし、この５つの方向が［０：２π］の範囲にあるとする。図３４〜図３６に示すように、回転して方位を変えながら連続して計測した時の測定をｓ＝０，１，２，…とする。なお、図３４〜図３６において、ｓ＝０，１，２，…の各回での回転方向は同一方向でも良いし、交互に逆回転させても良い。

図３４及び図３５では、例えばｋ番目のサイクルは、Ｓ＝１４の測定の後、一サイクルの成立が判定され、ｋサイクル区間が成立する。成立したｋサイクル区間について、後述する誤差評価プロセスを行う。一サイクルの最初と最後が同じ方向である。図３４では全てのサイクルの最初の方向は方向ａで同じであるが、図３５での隣り合うサイクルは、サイクルの最初の方向が異なる。図３６でも、サイクルの最初と最後が同じ方向であるが、一サイクル中に５つの方向（ａ，ｂ，ｃ，ｄ，ｅ，ｆ）が各２回ずつ設定されている。

なお、上述した通り、加速度センサーの検出軸を、図３０に示すNorth-East-Down（NED）座標系に方位が関連付けて位置決定された第１測定方位を基準方位とすることができる。重力方向と交差する仮想平面に検出軸を有する加速度センサーの検出軸を、概ね水平面とする。図３０に示すように、加速度センサーの基準方向出力の水平面への投影方向とNED座標系の基準とする方向のなす角をψとする。回転計測の測定方位の基準方向を、加速度センサーの設置された環境のNorth-East-Down（NED）座標系に関連付けて方位を決定する。また、図３１に示すように、計測条件はランダムドリフトレートＤ_Ｒが最小となるＡ-Ｂ区間の計測時間とするのが望ましい。

４．５．２． 処理手順２（動的誤差の測定）
測定結果を次のモデル式でそれぞれのパラメータが表される関係とする。測定データからそれぞれの項が見積もられる。

ここで、図３７には上記式中の傾斜方位角θφ、測定方位と傾斜方位角の差ψｎ、測定方位角θｎが示されている。また、重力方向と交差する仮想平面Ｐ_φ（図３７）に検出軸を有する加速度センサーが、平面Ｐ_φの法線周りに回転した時の、水平面Ｐ_Ｈに対して傾斜している最大角度を図３８の通り傾斜角φとし、この時加速度センサーの検出軸に出力される重力加速度成分を傾斜量Ｉφとする。重力加速度の傾斜角φによる傾斜量Ｉφは、次の通り上述した式（１）と同じとなる。

ここで、傾斜量Ｉφは、傾斜の方向と加速度センサ出力軸を基準方位とした時の仮想平面Ｐ_φ上の方向との関係を、基準方向に対して以下の式の通りとする。

計測した1サイクルのデータについて、測定データの順番ｍを、ｍ＝０，１，２，…，ｍ_ｍとし、ｍ回目の加速度センサー出力をＤｍとし、測定の方位角をθｍとする。１サイクルのデータ区間の最初の測定方向（基準方位）をｍ＝０とし、最後の測定方向（基準方位）をｍ＝ｍ_ｍとする。測定単位区間の基準方向の最初のデータｍ＝０と最後のデータｍ＝ｍ_ｍの差を求め、その差をそれぞれの測定データに測定中に、均等に配分して、ｍ回目の測定中の動的誤差を求める。ｍ回目の加速度センサーの出力Ｄｍから動的誤差を引き算する。

上述の通り動的誤差を減算したデータについて、１サイクルの測定区間の中で測定方位が重複している測定値は、同じ方位で平均し、測定方位角毎のデータすることができる。
測定方位毎のデータは、測定方位番号ｎで区分された測定値Ｄｎと測定方位角θｎの組合せとし、それらを図３９に示す1次加工データＤｎ’とする。

４．５．３． 処理手順３（静的誤差の測定）
次に、一次加工データの平均値を計算して、静的誤差Ｄ_Ｆを算出する。ただし、測定方位番号ｎ＝０，１，２，…，ｎ_ｎであり、測定方位数はｎ_ｎである。

４．５．４． 処理手順４（二次加工データの生成）
一次加工データＤｎ’から静的誤差Ｄ_Ｆを減算して二次加工データＤｎ”を生成する。この二次加工データを図４０に示す。

こうして求められた二次加工データＸ”ｍは、測定値Ｘから動的誤差Ｄｘ_ΔＴΔＴ及び静的誤差ＤＦｘが減算されたものであり、式（７）からＸ”ｍ＝Ｉφ・ｃｏｓθが傾斜誤差として求められる。同様に、式（８）からＹ”ｍ＝Ｉφ・ｓｉｎθが傾斜誤差として求められる。

４．５．５． 処理手順５（二次加工データを用いた計算）
測定方位角θｎと二次加工データＤｎ”の関係は次式で表される。

測定方位θｎは測定条件なので既知数とし、傾斜量Ｉ_φ、傾斜方位角θ_φを未知数Ａ，Ｂとして、任意のＡ，Ｂから二次加工データＣ”ｎを推測する。

推定された二次加工データＣ”ｎと計測された二次加工データＤｎ”との比Ｒｎを求める。

未知数Ｂを任意の範囲で調整して、Ｂが傾斜方位角と等しくなる場合（Ｂ＝θφ）、次の式が成立し、それぞれの測定方位角におけるこの比は一定の値になる。

これを利用して、任意のＢについて二次加工データのそれぞれの測定方位角θｎにおける比Ｒｎの分散σ_Ｒを計算し、この分散σ_Ｒが最小となるＢ_ｍｌを、傾斜方位角θ_φの最尤値とする。

これより傾斜量Ｉ_φを算出する。求めた傾斜方位角θ_φの最尤値Ｂ_ｍｌの時の比Ｒｎの相加平均は、未知数Ａを1とした時、次のように示される。

Ｂ＝Ｂ_ｍｌの時の相加平均を傾斜量Ｉ_φの最尤値とする。図４１に示すように、未知数Ｂを全周の範囲でこの分散σ_Ｒが最小となるＢ_ｍｌを求めると、１８０度方位差の２方向が求まる。従って、未知数Ｂを任意の範囲振って分散σ_Ｒを計算する区間は［０〜π］の区間とし、この区間で最尤値となる方位角にπを加算した方位角の近傍の最尤値を算出する。測定値には誤差が含まれるので、算出した２傾斜方位角の平均を求めて、方位角方向とします。この２つの平行な方向に傾斜方位角があり、傾斜の仰角方向を傾斜方位角方向とすると、この算出された２傾斜方向の相加平均が正の整数区間にある方位を傾斜方位角として選択する。

４．５．６． 処理手順６（平均値の計算）
それぞれの測定サイクルをｋ＝１，２，…，ｋ_ｋ回繰り返して計測したとき、測定サイクルｋで得られた計測値を次の通り定義する。
測定精度を向上させるため、全サイクルｋ＝１〜ｋ_ｋでの計測値の平均を求める。

４．５．７． 処理手順７（温度依存性の評価）
それぞれの測定サイクルをｋ＝１，２，…として、サイクルの測定区間の温度センサー出力Ｔｍの平均を計算し、Ｔｍをサイクル区間の測定温度とする。ｋサイクルの測定温度平均値Ｔｋは次の通りとなる。

それぞれの測定サイクルｋの平均温度Ｔｋと静的誤差の計測値Ｄ_Ｆｋとの相関から、静的誤差Ｄ_Ｆｋの温度依存性は、図４２に示す通りに評価される。図４２に示す相関は、次の通り、各項毎に温度特性補正係数を有する二次関数以上（ｈ≧２）の多項式で近似することができる。

４．５．８． 処理手順８（測定結果に対する補正）
加速度センサーの入出力を次の通り定義して、補正を行うことができる。

バイアス補正は、加速度センサーの入出力を次の通り定義して、補正を行うことができる。なお、温度特性相関を更新するために、較正測定時には温度特性補正を行わない。

４．５．９． 更正値の更新処理
次に、図４３を参照して、上述した処理手順１〜７に従った更正値の更新動作につて説明する。先ず、加速度センサーに関して、図４に示す水平座標系への関連付け、測定条件等が初期設定される（ステップ１）。個別の測定条件として、例えば図３１に示すようにランダムドリフトレートＤ_Ｒが最小となる測定条件を設定し（ステップ２）、図３１に示す区間Ａ−Ｂを初期設定する（ステップ３）。その後、図３３〜図３６に示すようにして加速度センサーを順次回転静止させて測定を行う（ステップ４）。各回転静止毎に、図３３〜図３６に示す一サイクル区間の定義の判定を実施する（ステップ５）。図３３〜図３６に示す一サイクル区間が成立していなければ（ステップ６がＮＯ）、ステップ４に戻る。一サイクル区間が成立していれば（ステップ６がＹＥＳ）、ステップ７に移行する。ステップ７では、加速度センサーの出力Ｄｍから動的誤差を引き算して一次加工データを生成し、さらに一サイクル区間で一次加工データの平均値を生成する。さらに、一次加工データから静的誤差を引き算して二次加工データを求め、二次加工データに基づいて各種の補正値を計算する（ステップ８）。全サイクルについて回転静止位置での測定が終了したか否かが判断され（ステップ９）、ステップ９での判断がＮＯであればステップ４に戻り、ステップ９での判断がＹＥＳであればステップ１０に移行する。ステップ１０では、統計計算として例えば全サイクルｋ＝１〜ｋ_ｋでの計測値の平均が求められ、その統計値が更正値として更新される（ステップ１１）。

４．５．１０． 傾斜角等の測定
図４４は、上述した処理手順１〜８に従った傾斜角、傾斜方位角、傾斜量等の傾斜計測定動作について説明する。図４４において、ステップ１〜ステップ１１は図４３と同じであり、ステップ１２が追加されている。図４４において、ステップ１１にて更正値が更新された後に、傾斜角、傾斜方位角、傾斜量等の傾斜計測定動作が終了したか否かが判断される（ステップ１２）。ステップ１２での判断がＮＯであれば、ステップ４に戻って、ステップ４〜ステップ１１による傾斜角、傾斜方位角、傾斜量等の傾斜計測定動作が行われる。ここで、図４４のステップ１０での統計計算では、測定データ及び更正値を用いて傾斜方位角の平均値θ_φＲ、傾斜量の平均値Ｉ_φＲ、傾斜角の平均値φ_Ｒが求められる。

なお、上記のように本実施形態について詳細に説明したが、本発明の新規事項及び効果から実体的に逸脱しない多くの変形が可能であることは当業者には容易に理解できるであろう。従って、このような変形例は全て本発明の範囲に含まれるものとする。

例えば、基準方位に設定した第１静止位置（例えば０度）から少なくとも２つの第２静止位置（例えば、９０度、１８０度、２７０度）を経て再度基準方位である第３静止位置（３６０度＝０度）にて測定する１サイクルとしたときには、図４５及び図４６に示すようにして複数サイクルで実施される。

図４５に示す例では、０度→９０度→１８０度と反時計回り方向に回転させた後に、１８０度→９０度→０度と時計回転方向に回転させ、逆方向の回転を交互に繰り返している。一方向に連続回転させると内部配線が破綻するからである。また、図４５の例では、１サイクル目（０度→９０度→１８０度→０度）と２サイクル目（１８０度→９０度→０→１８０度）とは、基準方位が０度、１８０度とそれぞれ異なっている。

一方、図４６に示す例では、０度→９０度→１８０→２７０度と反時計回り方向に回転させた後に、２７０度→０度と時計回転方向に戻し回転させている。また、図４６の例でも、１サイクル目（０度→９０度→１８０度→２７０度→０度）と、２サイクル目（９０度→１８０度→２７０度→０度→９０度）と、３サイクル目（１８０度→２７０度→０度→９０度→１８０度）、４サイクル目（２７０度→０度→９０度→１８０度→２７０度）とは、基準方位が０度、９０度、１８０度、２７０度とそれぞれ異なっている。

図４７は、測定開始指令から待ち時間を経てデータ収集する動作フローが示されている。端末ＰＣ２から測定端末装置３へ指令後に、測定端末装置３はデータ出力するが、端末ＰＣ２はそのデータを取り込まない。待ち時間は、上述した通り、環境温度変動が小さい、安定している時間帯を選択するために設けられている。待ち時間経過後に端末ＰＣに取り込まれるデータは、図４７に示すように、所定サイクル毎に平均化することで測定精度を高めることができる。

１ ホスト装置、２ 端末装置、３ 測定端末装置、４ 筐体、４Ａ 設置面、５ 回転操作部、６ インジケーター、７ 通信装置、８ 回転駆動部、１０ｘ，１０ｙ，１０ｚ 加速度センサー、２０ 記憶装置、３０ 回転司令部、４０ 誤差測定部、４１ 動的誤差測定部、４２ 静的誤差測定部、４３ 傾斜誤差測定部、５０ 傾斜度測定部、５１ 傾斜角測定部、５２ 傾斜方位角測定部、６０ 傾斜バイアス補正部、７０ 通信装置、ＸＹＺ センサー座標系、ＸｒＹｒＺｒ 絶対座標系、θ 傾斜方位角、φ 傾斜角、Ｉφ・ｃｏｓθ，Ｉφ・ｓｉｎθ 傾斜誤差、ＤＦｘ，ＤＦｙ 静的誤差、Ｄｘ_ΔＴΔＴ，Ｄｙ_ΔＴΔＴ 動的誤差、Ｘ０，Ｘ１，Ｘ２，…Ｙ０，Ｙ１，Ｙ２… 測定値、Ｘ０，Ｘ’１，Ｘ’２，…Ｙ’０，Ｙ’１，Ｙ’２… 一次加工データ

Claims (10)

1. 重力方向と交差する仮想平面に検出軸を有する加速度センサーの前記検出軸を、基準方位に設定した第１測定方位と、前記仮想平面に対する法線の周りの回転位置にて前記第１測定方位とは異なる少なくとも２つの第２測定方位と、前記基準方位に戻した第３測定方位と、の各々に設定した静止状態にて、３つ以上の異なる測定方位を含む計４つ以上の測定方位で測定される前記加速度センサーの出力を用いて、前記第１〜第３測定方位の各々で異なる動的誤差を測定する誤差測定工程と、
   前記第１測定方位及び前記少なくとも２つの第２測定方位の各々の測定方位に設定された前記加速度センサーの出力から、前記動的誤差と、前記第１〜第３測定方位の各々で実質的に一定の静的誤差とが除去された、前記仮想平面が水平面に対して傾斜する傾斜角及び傾斜方位角に依存するデータに基づいて、前記傾斜角及び前記傾斜方位角の少なくとも一方を測定する工程と、
   を有し、
   前記誤差測定工程は、前記第１測定方位及び前記第３測定方位に設定された前記加速度センサーの出力の差分に基づいて、前記動的誤差を測定する工程を含むことを特徴とする傾斜度測定方法。
2. 請求項１に記載の傾斜度測定方法において、
   前記第１〜第３測定方位の各々での前記動的誤差を、前記差分を線形補完して求めることを特徴とする傾斜度測定方法。
3. 請求項１または２に記載の傾斜度測定方法において、
   前記誤差測定工程は、
   前記第１測定方位と前記少なくとも２つの第２測定方位との各々の測定方位で測定された前記加速度センサーの出力の各々から、前記各々の測定方位について測定された前記動的誤差がそれぞれ減算された４つ以上の一次加工データを生成する工程と、
   前記４つ以上の一次加工データに基づいて、前記第１〜第３測定方位の各々で実質的に一定の静的誤差を測定する工程と、
   をさらに含むことを特徴とする傾斜度測定方法。
4. 請求項３に記載の傾斜度測定方法において、
   前記誤差測定工程は、前記４つ以上の一次加工データを平均化して、前記静的誤差を求めることを特徴とする傾斜度測定方法。
5. 請求項１〜４のいずれか一項に記載の傾斜度測定方法において、
   前記加速度センサーの前記検出軸を回転角９０°または１８０°で回転させることを特徴とする傾斜度測定方法。
6. 請求項１〜５のいずれか一項に記載の傾斜度測定方法において、
   前記第１〜第３測定方位での測定を１サイクルとしたとき、サイクル数ｎを２以上とする場合、ｎサイクルに亘って測定し、ｎサイクル目の前記第１測定方位を（ｎ−１）サイクル目の前記第３測定方位として兼用することを特徴とする傾斜度測定方法。
7. 請求項１〜６のいずれか一項に記載の傾斜度測定方法において、
   前記加速度センサーの前記検出軸は、前記第１測定方位から前記少なくとも２つの第２測定方位に向けて第１方向に回転された後に、前記少なくとも２つの第２測定方位から前記第３測定方位に向けて、前記第１方向とは逆方向の第２方向に回転されることを特徴とする傾斜度測定方法。
8. 重力方向と交差する仮想平面に検出軸を有する加速度センサーの前記検出軸を、基準方位に設定した第１測定方位と、前記仮想平面に対する法線の周りの回転位置にて前記第１測定方位とは異なる少なくとも２つの第２測定方位と、前記基準方位に戻した第３測定方位と、の各々に設定一次して、前記検出軸を３つ以上の異なる測定方位を含む計４つ以上の測定方位に設定する指令を出力する手段と、
   前記第１〜第３測定方位にそれぞれ設定された前記加速度センサーの出力を受信する手段と、
   前記加速度センサーの出力を用いて、前記第１〜第３測定方位の各々で異なる動的誤差を測定する誤差測定手段と、
   前記前記第１測定方位及び前記少なくとも２つの第２測定方位の各々の測定方位に設定された前記加速度センサーの出力から、前記動的誤差と、前記第１〜第３測定方位の各々で実質的に一定の静的誤差とが除去された、前記仮想平面が水平面に対して傾斜する傾斜角及び傾斜方位角に依存するデータに基づいて、前記傾斜角及び前記傾斜方位角の少なくとも一方を測定する手段と、
   を有し、
   前記誤差測定手段は、前記第１測定方位及び前記第３測定方位に設定された前記加速度センサーの出力の差分に基づいて、前記動的誤差を測定することを特徴とする傾斜度測定装置。
9. 請求項８に記載の傾斜度測定装置を有することを特徴とする電子機器。
10. 重力方向と交差する仮想平面に検出軸を有する加速度センサーの前記検出軸を、基準方位である第１測定方位に設定し、前記仮想平面に対する法線の周りの回転位置にて前記検出軸を前記第１測定方位とは異なる少なくとも２つの第２測定方位に設定し、前記仮想平面に対する法線の周りの回転位置にて前記検出軸を前記基準方位である第３測定方位に戻す指令を出力する手順と、
    前記第１〜第３測定方位に設定された前記加速度センサーからの出力を受信する手順と、
    前記第１〜第３測定方位の各々の静止状態で、前記検出軸を３つ以上の異なる測定方位を含む計４つ以上の測定方位に設定して測定される前記加速度センサーの出力を用いて、前記第１〜第３測定方位の各々で異なる動的誤差を測定する手順と、
    前記第１測定方位及び前記少なくとも２つの第２測定方位の各々の測定方位に設定された前記加速度センサーの出力から、前記動的誤差と、前記第１〜第３測定方位の各々で実質的に一定の静的誤差とが除去された、前記仮想平面が水平面に対して傾斜する傾斜角及び傾斜方位角に依存するデータに基づいて、前記傾斜角及び前記傾斜方位角の少なくとも一方を測定する手順と、
    をコンピューターに実施させ、
    前記動的誤差を測定する手順は、前記第１測定方位及び前記第３測定方位での前記加速度センサーの出力の差分に基づいて、前記動的誤差を測定する手順を含むことを特徴とするプログラム。