JPH11211474A - 姿勢角検出装置 - Google Patents

姿勢角検出装置

Info

Publication number
JPH11211474A
JPH11211474A JP10023742A JP2374298A JPH11211474A JP H11211474 A JPH11211474 A JP H11211474A JP 10023742 A JP10023742 A JP 10023742A JP 2374298 A JP2374298 A JP 2374298A JP H11211474 A JPH11211474 A JP H11211474A
Authority
JP
Japan
Prior art keywords
angle
axis
attitude
geomagnetic
output
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP10023742A
Other languages
English (en)
Other versions
JP3763435B2 (ja
Inventor
Hiroshi Abe
洋 阿部
Kazutake Mutou
和健 武藤
Naoyuki Sato
尚之 佐藤
Naoharu Yamamoto
直治 山本
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tokin Corp
Original Assignee
Tokin Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tokin Corp filed Critical Tokin Corp
Priority to JP02374298A priority Critical patent/JP3763435B2/ja
Publication of JPH11211474A publication Critical patent/JPH11211474A/ja
Application granted granted Critical
Publication of JP3763435B2 publication Critical patent/JP3763435B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Gyroscopes (AREA)
  • Position Input By Displaying (AREA)

Abstract

(57)【要約】 【課題】 使用する場所や環境の制限を緩和し、幅広く
使用できる移動体の姿勢制御・検出に用いる姿勢角検出
装置を提供すること。 【解決手段】 3個のジャイロスコープ5,6,7と、
2個の加速度センサ10,11と、2個の地磁気センサ
8,9を配置し、前記3個のジャイロスコープが出力す
る角速度に基づいて運動角検出する運動角演算装置31
と、加速度及び地磁気に基づいて姿勢角を検出する静止
角演算装置32と、静止角検出結果の正誤を判別する判
別装置33と、この判別装置の演算結果に応じて運動角
検演算装置と静止角演算装置の演算結果を用いて出力す
べき姿勢角度を演算する姿勢角演算装置30を備え、外
部の信号発生源を使用することなく移動体の姿勢角度を
検出する。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、移動体の姿勢検出
や姿勢制御・バーチャルリアリティー等に使用されるヘ
ッドマウントディスプレーのうち、頭の姿勢角度を検出
するトラッカー、3Dゲームパッド等、3次元空間中で
被測定対象物の姿勢角度(3軸の回転角度)を検出する
のに特に適した姿勢角検出装置に関する。
【0002】
【従来の技術】従来、バーチャルリアリティーに使用さ
れる、頭の動きを検出する姿勢角検出方法として、外部
に設置された交流磁場発生源から微弱交流磁場を発生
し、ヘッドマウントディスプレイ(以下、HMDとい
う)に配置したセンサ部で検知し、制御演算部で頭の動
きを検出する交流磁場を利用する方法や、外部に設置し
た超音波発生源からの超音波信号をHMDに配置したセ
ンサ部で検知し、制御演算部で頭の動きを検出する超音
波を利用する方法がある。
【0003】
【発明が解決しようとする課題】しかしながら、前者の
交流磁場を利用した方法では、信号発生源が微弱交流磁
場のため、ノイズをキャンセルする為にフィルタを多用
することにより応答性が低下し、頭の動きに比べて、H
MDの映像の動きが遅くなり、気分を悪くしたり、酔い
を発生する場合がある。
【0004】また、後者の超音波を利用した方法では、
他の様々な超音波信号の影響を受けての誤動作や、信号
発生源とセンサとの間の遮蔽物により信号を検出できな
くなる可能性があり、センサの前に腕や髪の毛があるだ
けで受信不能になるので、このための誤動作を防止する
ための対策を必要とする等の問題があった。
【0005】本発明は、上記問題点を解消し、使用する
と場所や環境の制限を緩和し、幅広く使用できる移動体
の姿勢制御・検出に用いる姿勢角検出装置を提供するも
のである。
【0006】
【課題を解決するための手段】前記課題を解決するた
め、本発明に係る姿勢角検出装置は、互いに直交する3
軸の回りの角速度を検出する第1、第2及び第3のジャ
イロスコープと、互いに直交する2軸の加速度を検出す
る第1及び第2の加速度センサと、互いに直交する2軸
の地磁気を検出する第1及び第2の地磁気センサを配置
し、前記ジャイロスコープの出力する角速度に基づいて
運動角を検出する運動角検出手段と、加速度及び地磁気
に基づいて姿勢角を検出する静止角検出手段と、静止角
検出手段の検出結果が正しいかどうかを判別する判別装
置と、この判別装置の判別結果と、運動角検出手段と静
止角検出手段の演算結果を用いて、出力すべき姿勢角度
を演算する姿勢角演算装置を備え、外部の信号発生源を
使用することなく移動体の姿勢角度を検出できるように
したことを特徴とする。
【0007】即ち、本発明は、水平面内で直交する軸を
X軸、Y軸とし、そのX軸、Y軸に直交する軸をZ軸と
するとき、これら互いに直交するX軸、Y軸、Z軸が形
成する基準座標系において、互いに直交する3軸の回り
の回転角を検出する姿勢角検出装置において、X軸、Y
軸及びZ軸の回りの角速度を検出する3個のジャイロス
コープと、該ジャイロスコープの角速度に応じた出力に
基づいて単位時間に移動した角度を演算する運動角演算
装置と、X−Y平面で互いに直交する2軸の加速度を検
出するように配置された加速度センサと、X−Y平面で
互いに直交する2軸の地磁気を検出するように配置され
る地磁気センサと、該加速度センサ及び地磁気センサの
出力に基づいて、X軸、Y軸、Z軸回りの回転角を演算
する静止角演算装置と、静止角演算装置による演算結果
の真偽を判別する判別装置と、該判別装置の演算結果に
応じて、運動角演算装置の演算結果と静止角演算装置の
演算結果から出力するべき姿勢角を演算する姿勢角演算
装置とから構成される姿勢角検出装置である。
【0008】また、本発明は、上記姿勢角検出装置のジ
ャイロスコープが、圧電振動ジャイロ用振動子と、その
振動子をジャイロスコープとして動作させるための駆動
検出回路からなり、該振動子が、非磁性材料で構成され
る姿勢角検出装置である。
【0009】また、本発明は、上記の姿勢角検出装置の
ジァイロスコープの出力部に、ハイパスフィルタが接続
され、ハイパスフィルタの出力部は演算装置に入力さ
れ、ハイパスフィルタの出力を用いて運動角を演算する
姿勢角検出装置である。
【0010】また、本発明は、上記ハイパスフィルタ
が、カットオフ周波数可変手段を有する姿勢角検出装置
である。
【0011】また、本発明は、上記姿勢角検出装置にお
いて、2軸の加速度センサと、加速度センサの出力Ax
(n)、Ay(n)から仮のロール角R(n)、仮のピ
ッチ角P(n)を演算するロール角演算手段、ピッチ角
演算手段と、地磁気のx方向成分Mx(n)、y方向成
分My(n)を計測する2軸の地磁気センサと、演算す
る単位時間で1単位時間過去の方位角演算結果Φ(n−
1)を記憶しておくためのバックアップメモリーと、加
速度センサの出力Ax(n)、Ay(n)及びバックア
ップメモリーの内容Φ(n−1)より、傾斜角度の北方
向成分θ2 を演算する北方向成分傾斜角演算手段と、地
磁気のZ方向成分Mzの符号を判定する地磁気Mz符号
判別手段と、地磁気のZ方向成分地磁気Mz絶対値演算
手段と、前記地磁気Mz符号判別手段と、地磁気Mz絶
対値演算手段よりZ軸方向成分の地磁気Mzを計算する
地磁気Mz演算手段と、前記地磁気Mx、My、Mzよ
り、基準座標系での地磁気のX,Y,Z成分HX,H
Y,HZを演算する座変換演算手段と、基準座標系にお
ける前記のHX、HYにより方位角Φを演算する方位演
算手段から構成され、地磁気Mz演算手段により仮想的
に地磁気Mzを求めることにより地磁気センサ2軸によ
り方位角Φ(n)を計測する姿勢角検出装置である。
【0012】また、本発明は、上記姿勢角検出装置にお
いて、前記北方向成分傾斜角演算装置及び地磁気Mz符
号判別手段により、北方向成分傾斜角θ2、原点時の地
磁気の伏角θ0、バックアップメモリの出力をΦ(n−
1)、加速度センサ2軸の出力をそれぞれAx(n)、
Ay(n)とする時、北方向成分傾斜角θ2は、θ2=s
in-1[Ax(n)・cos(−Φ(n−1))+Ay
(n)・sin(−Φ(n−1))]として求め、−9
0度≦θ2≦θ0であれば、地磁気Mz符号はマイナス、
θ0 <θ2 ≦90度であれば、地磁気Mz符号はプラス
と判別する姿勢角検出装置である。
【0013】また、本発明は、上記姿勢角検出装置にお
いて、前記判別装置により、2軸の加速度センサおよ
び、2軸の地磁気センサからなる静止角検出手段から求
められる仮のヨー角Φ、仮のピッチ角P、仮のロール角
Rの検出結果が、正しいかどうかの判別手段として、加
速度センサの現在の出力Ax(n)から演算して得られ
た仮のロール角をR(n)、演算する単位時間で1単位
時間過去の出力Ax(n−1)から得られた仮のロール
角演算結果をR(n−1)とし、同様に加速度センサの
現在の出力Ay(n)から演算して得られた仮のピッチ
角をP(n)、演算する単位時間で1単位時間過去の出
力Ay(n−1)から得られた仮のピッチ角をP(n−
1)とすると、数1あるいは数2のうち、少なくとも一
つの値が、ある値以内にあるとき、静止角検出手段から
の姿勢角情報は正しいと判別する姿勢角検出装置であ
る。
【0014】
【数1】
【0015】
【数2】
【0016】また、本発明は、上記姿勢角検出装置にお
いて、求めるべき出力をロール角α(n)、ピッチ角β
(n)、ヨー角γ(n)、運動角検出手段から求められ
る運動角を、仮のロール角変化△X(n)、仮のピッチ
角変化△Y(n)、仮のヨー角変化△Z(n)、静止角
検出手段から求められる静止角を仮のロール角R
(n)、仮のピッチ角P(n)、仮のヨー角Φ(n)、
誤差より大きくない一定定数をC1,C2,C3とする
と、判別装置の判別結果により、静止角検出手段からの
姿勢角情報は正しいと判別する場合、α(n)=α(n
−1)+△X(n)−C1、β(n)=β(n−1)+
△Y(n)−C2、γ(n)=γ(n−1)+△Z
(n)−C3により、また、静止角検出手段からの姿勢
角情報は誤りと判別する場合は、α(n)=α(n−
1)+△X(n)、β(n)=β(n−1)+△Y
(n)、γ(n)=γ(n−1)+△Z(n)により、
姿勢角を演算する姿勢角検出装置である。
【0017】また、本発明は、上記姿勢角検出装置にお
いて、求めるべき出力をロール角α(n)、ピッチ角β
(n)、ヨー角γ(n)、動的な運動角検出手段から求
められる運動角を、仮のロール角変化△X(n)、仮の
ピッチ角変化△Y(n)、仮のヨー角変化△Z(n)、
静的な静止角検出手段から求められる姿勢角を仮のロー
ル角R(n)、仮のピッチ角P(n)、仮のヨー角Φ
(n)、1以下の比例定数をk1,k2,k3とすると、判
別装置の判別結果により、静止角検出手段からの姿勢角
情報は正しいと判別する場合は、α(n)=α(n−
1)+△X(n)−k1[α(n−1)+△X(n)−
R(n)]、β(n)=β(n−1)+△Y(n)−k
2[β(n−1)+△Y(n)−P(n)]、 γ(n)
=γ(n−1)+△Z(n)−k3[γ(n−1)+△
Z(n)−Φ(n)]により、また、静止角検出手段か
らの姿勢角情報は誤りと判別する場合は、α(n)=α
(n−1)+△X(n)、β(n)=β(n−1)+△
Y(n)、γ(n)=γ(n−1)+△Z(n)によ
り、姿勢角を演算する姿勢角検出装置である。
【0018】
【発明の実施の形態】以下、図面によって本発明の実施
の形態について説明する。
【0019】ここで説明するセンサ座標系から基準座標
系への座標変換演算は、ヨー角γ→ピッチ角α→ロール
角βの順に行うものとして説明する。
【0020】ここで、センサ座標系でのセンサ出力とし
て、次の記号を用いる。 ジャイロスコープの出力:Jx(n)、Jy(n)、J
z(n)。 加速度センサの出力:Ax(n)、Ay(n)。 地磁気センサの出力:Mx(n)、My(n)。 地磁気センサの計算値:Mz(n)。
【0021】また、上記センサ座標系でのセンサ出力を
基準座標系に座標変換して得た値として、次の記号を用
いる。 ジャイロスコープの出力から得た基準座標系X、Y、Z
軸回りの運動角:△X(n)、△Y(n)、△Z
(n)。 加速度センサの出力から得た基準座標系X軸回りの回転
角(仮のロール角):R(n)、同じくY軸回りの回転
角(仮のピッチ角):P(n)。 地磁気センサの出力から得たZ軸回りの回転角(方位角
或いは仮のヨー角):Φ(n)。 上記演算結果より求める最終的に出力するべき基準座標
系X、Y、Z軸回りの回転角度ロール角α(n)、ピッ
チ角β(n)、ヨー角γ(n)とする。
【0022】図1には、基準座標系(X−Y−Z系)と
センサ座標系(x−y−z系)が一致している場合につ
いて示す。図1に示すように、水平面で直交する軸をX
軸、Y軸とし、そのX軸、Y軸それぞれに直交する軸を
Z軸とするとき、これら互いに直交するX軸、Y軸、Z
軸が形成する基準座標系において、X軸回りの回転角度
をロール角α、Y軸回りの回転角度をピッチ角β、Z軸
回りの回転角度をヨー角γと記す。
【0023】図1には、本発明の姿勢角検出装置におけ
るジャイロスコープ、加速度センサ、地磁気センサの配
置を示す。図1で示すセンサ座標系において、互いに直
交する3軸(x軸、y軸及びz軸)の回りの角速度を検
出するための第1のジャイロスコープ5、第2のジャイ
ロスコープ6、及び第3のジャイロスコープ7は、x、
y、z軸各軸に平行に、即ち互いに直交して、配置され
ている。第1の加速度センサ10と第2の加速度センサ
11は、互いに直交する2軸、x軸及びy軸に平行に配
置されている。同様に、ヨー角を検出するための第1の
地磁気センサ8と第2の地磁気センサ9は、互いに直交
する2軸、x軸及びy軸に平行に配置されている。
【0024】図2には、基準座標系(X−Y−Z系)に
対して、センサ座標系(x−y−z系)が、ずれた場合
について示す。センサ座標系に記載されるべき、ジャイ
ロスコープ、加速度センサ、地磁気センサは、省略して
ある。図2では、基準座標系の−Z方向に重力加速度が
働いており、センサ座標系のx軸、y軸方向への分力
が、Ax(n)、Ay(n)であることを示している。
【0025】図3は、本発明の一実施の形態における姿
勢角検出装置の全体構成を示す図である。図3に示すよ
うに、本発明の姿勢角度検出装置は、図1に示すよう
に、直交配置された3個のジャイロスコープと、そのジ
ャイロスコープそれぞれに接続されるハイパスフィルタ
と、ハイパスフィルタを介したジャイロスコープの出力
から単位時間あたりの移動角度(運動角)を演算する運
動角演算装置31と、図1に示すように配置された2個
の加速度センサと、2個の地磁気センサと、前記2個の
加速度センサから仮のロール角R、仮のピッチ角Pを求
め、さらにその仮のロール角R、仮のピッチ角Pから仮
のヨー角Φを求める静止角演算装置32と、静止角演算
装置32の演算結果の真偽を判定するための判別装置3
3と、運動角演算装置31、静止角演算装置32及び判
別装置33からの信号を処理する姿勢角演算装置30と
で構成されている。
【0026】最初に、運動角検出装置31の内容を説明
する。本発明の姿勢角検出装置における運動角検出装置
は、圧電セラミック、圧電単結晶或いはシリコン等の非
磁性材料からなる振動子を用いることが可能である。本
実施例では、(株)トーキン製圧電振動セラミックジャ
イロを用いた場合について説明する。
【0027】図3に示すように、3軸のジャイロスコー
プ5,6,7の出力部には、それぞれ、ハイパスフィル
タ36、37、38が接続され、ハイパスフィルタの出
力部に運動角演算装置31が接続されている。ジャイロ
スコープは、小型で安価な物ほど静止しているときの出
力(オフセット)のばらつき及び変動が大きいが、安定
したものを用いると、姿勢検出装置が高価で大型になっ
てしまう。
【0028】前記ハイパスフィルタは、前記圧電振動ジ
ャイロのオフセットをキャンセルするもので、ハイパス
フィルタを用いた後の出力を用いることで、安価で小型
で高精度の姿勢角検出装置を得ることができるもので、
通常0.1Hz以下の低い周波数に設定される。また、
それぞれのハイパスフィルタは、カットオフ周波数可変
を有している。前記低周波のハイパスフィルタは、電源
投入時から安定するまで時間がかかるので、必要に応じ
てハイパスフィルタのカットオフ周波数を高くすること
により、安定した出力を短時間で得ることができる。カ
ットオフ周波数を高くし安定出力を得た後は、即時にカ
ットオフ周波数は低い状態(測定状態)に戻される。
【0029】運動角演算装置31によって、姿勢角検出
装置(被測定物)の回転角速度に応じた各ジャイロスコ
ープの出力Jx(n)、Jy(n)、Jz(n)は、座
標変換を行って、演算する単位時間で1単位時間過去の
出力β(n−1)とから、次式の如く、基準座標系にお
ける現在の移動角度(運動角)△X(n)に変換され
る。
【0030】 ΔX(n) = cos[β(n-1)]×Jx(n) + sin[β(n-1)]×Jz(n) ΔY(n) = sin[α(n-1)]×sin[β(n-1)]×Jx(n)+cos[α
(n-1)]×Jy(n)-sin[α(n-1)]×cos[β(n-1)]×Jz(n) ΔZ(n) = -cos[α(n-1)]×sin[β(n-1)]×Jx(n)+sin[α
(n-1)]×Jy(n)+cos[α(n-1)]× cos[β(n-1)]×Jz(n)
【0031】次に、図4の静止角演算装置の説明を行
う。静止角演算装置は、先にも述べたように、姿勢角検
出装置の静的動作時の姿勢角である仮のロール角R、仮
のピッチ角P、仮のヨー角Φを出力するものである。
【0032】図1に示す加速度センサ10、11は、図
2に示すように、それぞれ、センサ座標系でのx、y軸
上の重力加速度のx、y成分分力Ax(n)、Ay
(n)を検出するように配置されている。また、図1に
示す地磁気センサ8、9は、それぞれx、y軸上の地磁
気分力Mx(n)、My(n)を検出するように配置さ
れている。
【0033】加速度センサの出力Ax(n)、Ay
(n)は、各々ロール角演算手段、ピッチ角演算手段に
て下式により、基準座標系の水平面(XY平面)との傾
斜角である仮のロール角R(n)及び仮のピッチ角P
(n)を算出する。
【0034】 R(n)=sin-1[Ax(n)/cosP(n)]
P(n)=sin-1Ay(n)
【0035】基準座標系での仮のヨー角Φ(n)は、基
準座標系での地磁気成分を用いて、Φ(n)=tanー1
[HX(n)/HY(n)]で表わされる。姿勢角検出
装置が水平面上で回転する場合は、地磁気センサの出力
Mx(n)及びMy(n)が、それぞれ、HX(n)及びHY
(n)に等しく、 Φ(n)=tanー1[Mx(n)/My(n)] が成り立つ。
【0036】原点時の方位角をΦ(0)とすれば、ヨー
角γ(n)は、 γ(n)=Φ(n)−Φ(0) となる。
【0037】しかし、姿勢角検出装置がロール角Rやピ
ッチ角Pを持ち、水平でないときには、 HX=cos[P(n)]×Mx+sin[P(n)]
×Mz、 HY=sin[R(n)]×sin[P(n)]×Mx
+cos[R(n)]×cos[P(n)]×Mz の関係式から、仮のロール角R(n)、仮のピッチ角P
(n)を用いて座標変換を行い、基準座標系での地磁気
成分を求めて、 Φ(n)=tanー1[HX(n)/HY(n)] に代入し、求める。
【0038】その際、地磁気Z方向成分Mz(n)が必
要となる。
【0039】以下の方法で仮想的にMzを算出する。
【0040】まず、地磁気Htを以下の手順で求める。
最初に、姿勢角検出装置(被測定物)部を90度傾け、
図1の磁気センサ10あるいは磁気センサ11が基準座
標系のZ軸に向くようにする。そして何らかの信号を送
り、例えばキーやボタンを押すなどをし、その時のZ軸
に向いているセンサの値をMz(0)としてメモリに格
納する。
【0041】次に、姿勢角検出装置(被測定物)部を原
点位置に戻し、もう一度、何らかの信号を送り、磁気セ
ンサMx(n),My(n)の値をMx(0),My
(0)としてメモリに格納する。これらMx(0),M
y(0),Mz(0)を用いて、Htは数3から求ま
る。
【0042】
【数3】
【0043】さらに、姿勢角検出装置(被測定物)の地
磁気北方向成分の傾斜角からMzの符号を判別する。こ
のために、バックアップメモリーと北方向成分傾斜角演
算手段を有する。
【0044】演算処理にて算出した、演算する単位時間
で、1単位時間過去の仮のヨー角Φ(n-1)をバックアッ
プメモリーに記憶しておき、この仮のヨー角Φ(n−
1)と加速度センサ出力Ax(n)、Ay(n)の値か
ら、北方向成分傾斜角演算手段で、姿勢角検出装置(被
測定物)の傾斜角度の地磁気北方向成分θ2を次式によ
り算出する。
【0045】θ2=sin-1[Ax(n)×cosΦ(n−
1)+Ay(n)×sinΦ(n−1)]
【0046】により求められ、−90度≦θ2≦θ0であ
れば、地磁気Mz符号はマイナス、θ0<θ2≦90度で
あれば、地磁気Mz符号はプラスと判別する。
【0047】次に、本発明の姿勢角検出装置における、
前記静止角検出装置で演算された仮のロール角R
(n)、仮のピッチ角P(n)、仮のヨー角Φ(n)が
ほぼ正しいかどうかを判別する判別装置を説明をする。
【0048】静止角検出手段における、Ax(n)、A
y(n)は、重力加速度のx,y軸方向成分と運動加速
度のx、y軸方向成分との合成ベクトルとなるために、
運動加速度がある場合には正しいロール角、ピッチ角に
ならない。そこで、R(n)からR(n−1)を引いた
ものが、0の近傍で、ある値以内であり、あるいはP
(n)からP(n−1)を引いたものが、0の近傍で、
ある値以内であるならば、運動加速度がない、即ち静止
角検出手段から演算された仮のロール角R(n)、と仮
のピッチ角P(n)が正しいと判別する。そうでない場
合は、運動加速度があると判別し、静止角検出手段から
の演算結果は正しくないとする。即ち、数1或いは数2
のうち、少なくとも一つの値が、ある値以内で、ほぼ0
であるとき、静止角検出手段からの姿勢角情報は正しい
と判別する判別装置を提供する。
【0049】
【数1】
【0050】
【数2】
【0051】ここでは、加速度センサの出力を演算して
判別する手法について述べたが、磁気センサの出力に誤
差が生じ、大きく変動したとき等の場合も静止角度検出
手段の出力は正しくないと判別する。
【0052】次に、上記判別装置の判別結果に応じて、
運動角検出手段の出力と静止角検出手段の出力を演算し
て最終的な出力を得る姿勢角演算装置で行われる演算方
法について説明する。
【0053】この演算の基本は、判別装置によって静止
角検出手段の出力が誤っていると判別されたときには、
現在の姿勢角検出結果=1単位時間過去の姿勢角検出結
果+運動角となり、判別装置によって静止角検出手段の
出力が正しいと判別されたときには、現在の姿勢角検出
結果=1単位過去の姿勢角検出結果+運動角−補正値と
し、補正値は、誤差×比例定数、或いは一定定数であ
る。ここで、誤差=(1単位時間過去の姿勢角検出結果
+現在の運動角)−(静止角検出手段による現在の検出
結果)である。
【0054】式で表すと、X軸の誤差=α(n−1)+
△X(n)−R(n),Y軸の誤差=β(n−1)+△
Y(n)−P(n),Z軸の誤差=γ(n−1)+△Z
(n)−Φ(n)となる。
【0055】第1の実施例では、本発明の姿勢角検出装
置の求めるべき出力をα(n)、β(n)、γ(n)、
動的な運動角検出手段から得られる運動角を△X
(n)、△Y(n)、△Z(n)、静的な静止角検出手
段から求められる静止角をR(n)、P(n)、Φ
(n)とすると、
【0056】静止角検出手段からの姿勢角情報は正しい
と判別する場合、 α(n)=α(n−1)+△X(n)−C1、 β(n)=β(n−1)+△Y(n)−C2、 γ(n)=γ(n−1)+△Z(n)−C3、 ここで、C1,C2,C3は、任意に選択される誤差より大
きくない定数で、符号は、誤差の符号が、プラスの時、
プラス、マイナスの時、マイナスにより姿勢角を演算す
る。
【0057】また、静止角検出手段からの姿勢角情報が
誤と判別する場合、 α(n)=α(n−1)+△X(n)、 β(n)=β(n−1)+△Y(n)、 γ(n)=γ(n−1)+△Z(n)、 により姿勢角を演算する。
【0058】また、本発明の他の実施の形態として、本
発明の姿勢角検出装置の求めるべき出力をα(n)、β
(n)、γ(n)、動的な運動角検出手段から得られる
運動角を△X(n)、△Y(n)、△Z(n)、静的な
静止角検出手段から求められる静止角をR(n)、P
(n)、Φ(n)とすると、
【0059】静止角検出手段からの姿勢角情報は正しい
と判別する場合、 α(n)=α(n−1)+△X(n)−k1[α(n−1)+△
X(n)−R(n)]、 β(n)=β(n−1)+△Y(n)−k2[β(n−1)+△
Y(n)−P(n)]、 γ(n)=γ(n−1)+△Z(n)−k3[γ(n−1)+△
Z(n)−Φ(n)]、 により姿勢角を演算する。ここで、k1,k2,k3は、
任意に選択される1以下の比例定数である。
【0060】また、静止角検出手段からの姿勢角情報が
誤と判別する場合、 α(n)=α(n−1)+△X(n)、 β(n)=β(n−1)+△Y(n)、 γ(n)=γ(n−1)+△Z(n)、 により姿勢角を演算する。
【0061】このように、静止角検出手段の出力が正し
いときのみ、誤差を補正しながら、運動角の積算で姿勢
角度を求めることにより、高速な応答性と、静止時の安
定性、再現性を同時に満足できるようにした姿勢角検出
装置である。
【0062】図5は、姿勢角検出装置を適用したHMD
の一例示す。このHMDは、目前にあるディスプレイに
表示に表示された映像が頭の動きに連動して変化し、仮
想空間を体験できるように設けられているもので、この
HMDを装着したまま右を向くと、姿勢角検出装置で計
測した頭の姿勢角を映像発生装置に送信し、映像発生装
置はその姿勢角にあった映像をHMDに送信するので右
に展開する映像が映し出され、頭の動きによって全空間
の360度の映像を実感できるように設けられているも
のである。
【0063】図5に示すように、符号1は、HMD本体
を示し、このHMD本体1には、本発明である姿勢角検
出装置(被測定物)2が信号ケーブル3で映像発生装置
4に接続されている。
【0064】
【発明の効果】本発明による姿勢角検出装置によれば、
上記のHMDに適用した場合、外部信号を用いることな
く、また、高速応答で累積誤差の少ない角度情報を得る
ことができ、高性能なHMDを構成することができる。
【0065】本発明による姿勢角検出装置によれば、ジ
ャイロを素子で搭載しているので、小型軽量化が可能で
あり、セラミック振動子を使用しているので、磁気ノイ
ズは受けず、磁気センサと近接させても機能を損なうこ
とがない。また、ハイパスフィルタを使用することによ
り、小型・低価格化、高精度化が可能である。
【0066】さらに、従来、姿勢角検出装置としては、
加速度センサが3個ないと姿勢角の誤差が判断できず、
そのため、複雑な構造になっていたが、本発明によれ
ば、2個の加速度センサと、2個の地磁気センサと、簡
単な計算だけで、誤差成分有無の判別・補正処理を行う
ことができ、高速応答で、高精度な姿勢角検出装置を提
供できる。
【図面の簡単な説明】
【図1】本発明の姿勢角検出装置のジャイロスコープ、
加速度センサ、地磁気センサの配置を示し、基準座標系
(X−Y−Z系)とセンサ座標系(x−y−z系)が一
致している場合を示す図。
【図2】基準座標系(X−Y−Z系)とセンサ座標系
(x−y−z系)が一致していない場合を示す図。
【図3】本発明の姿勢角検出装置の一実施の形態の全体
構成を示す図。
【図4】本発明の姿勢角検出装置のうち、静止角演算装
置の構成を示すブロック図。
【図5】HMDにおける姿勢角検出装置の使用状況の説
明図。
【符号の説明】
1 ディスプレイ本体 2 姿勢角検出装置(被測定物) 3 信号ケーブル 4 映像発生装置 5 (第1の)ジャイロスコープ 6 (第2の)ジャイロスコープ 7 (第3の)ジャイロスコープ 8 第1の地磁気センサ 9 第2の地磁気センサ 10 第1の加速度センサ 11 第2の加速度センサ 31 運動角演算装置 32 静止角演算装置 33 判別装置 30 姿勢角演算装置 36 (第1の)ハイパスフィルタ 37 (第2の)ハイパスフィルタ 38 (第3の)ハイパスフィルタ
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.6 識別記号 FI G06F 3/033 310 G06F 3/033 310Y H04N 5/64 511 H04N 5/64 511A (72)発明者 山本 直治 宮城県仙台市太白区郡山6丁目7番1号 株式会社トーキン内

Claims (9)

    【特許請求の範囲】
  1. 【請求項1】 水平面内で直交する軸をX軸、Y軸と
    し、そのX軸、Y軸に直交する軸をZ軸とするとき、こ
    れら互いに直交するX軸、Y軸、Z軸が形成する基準座
    標系において、互いに直交する3軸の回りの回転角を検
    出する姿勢角検出装置において、X軸、Y軸及びZ軸の
    回りの角速度を検出する3個のジャイロスコープと、該
    ジャイロスコープの角速度に応じた出力に基づいて単位
    時間に移動した角度を演算する運動角演算装置と、X−
    Y平面で互いに直交する2軸の加速度を検出するように
    配置された加速度センサと、X−Y平面で互いに直交す
    る2軸の地磁気を検出するように配置される地磁気セン
    サと、該加速度センサ及び地磁気センサの出力に基づい
    て、X軸、Y軸、Z軸回りの回転角を演算する静止角演
    算装置と、静止角演算装置による演算結果の真偽を判別
    する判別装置と、該判別装置の演算結果に応じて、運動
    角演算装置の演算結果と静止角演算装置の演算結果から
    出力するべき姿勢角を演算する姿勢角演算装置とから構
    成されることを特徴とする姿勢角検出装置。
  2. 【請求項2】 請求項1記載の姿勢角検出装置のジャイ
    ロスコープが、圧電振動ジャイロ用振動子と、その振動
    子をジャイロスコープとして動作させるための駆動検出
    回路からなり、該振動子が、非磁性材料で構成されるこ
    とを特徴とする姿勢角検出装置。
  3. 【請求項3】 請求項1または2記載の姿勢角検出装置
    のジァイロスコープの出力部に、ハイパスフィルタが接
    続され、ハイパスフィルタの出力部は演算装置に入力さ
    れ、ハイパスフィルタの出力を用いて運動角を演算する
    ことを特徴とする姿勢角検出装置。
  4. 【請求項4】 請求項3記載のハイパスフィルタは、カ
    ットオフ周波数可変手段を有することを特徴とする姿勢
    角検出装置。
  5. 【請求項5】 請求項1乃至請求項4のいずれかに記載
    の姿勢角検出装置において、2軸の加速度センサと、加
    速度センサの出力Ax(n)、Ay(n)から仮のロー
    ル角R(n)、仮のピッチ角P(n)を演算するロール
    角演算手段、ピッチ角演算手段と、地磁気のx方向成分
    Mx(n)、y方向成分My(n)を計測する2軸の地
    磁気センサと、演算する単位時間で1単位時間過去の方
    位角演算結果Φ(n−1)を記憶しておくためのバック
    アップメモリーと、加速度センサの出力Ax(n)、A
    y(n)及びバックアップメモリーの内容Φ(n−1)
    より、傾斜角度の北方向成分θ2を演算する北方向成分
    傾斜角演算手段と、地磁気のZ方向成分Mzの符号を判
    定する地磁気Mz符号判別手段と、地磁気のZ方向成分
    地磁気Mz絶対値演算手段と、前記地磁気Mz符号判別
    手段と、地磁気Mz絶対値演算手段よりZ軸方向成分の
    地磁気Mzを計算する地磁気Mz演算手段と、前記地磁
    気Mx、My、Mzより、基準座標系での地磁気のX,
    Y,Z成分HX,HY,HZを演算する座標変換演算手
    段と、基準座標系における前記のHX、HYにより方位
    角Φを演算する方位演算手段から構成され、地磁気Mz
    演算手段により仮想的に地磁気Mzを求めることにより
    地磁気センサ2軸により方位角Φ(n)を計測すること
    を特徴とする姿勢角検出装置。
  6. 【請求項6】 請求項5記載の姿勢角検出装置におい
    て、前記北方向成分傾斜角演算装置及び地磁気Mz符号
    判別手段により、北方向成分傾斜角θ2、原点時の地磁
    気の伏角θ0、バックアップメモリの出力をΦ(n−
    1)、加速度センサ2軸の出力をそれぞれAx(n)、
    Ay(n)とする時、北方向成分傾斜角θ2は、θ2=s
    in-1[Ax(n)・cos(−Φ(n−1))+Ay
    (n)・sin(−Φ(n−1))]として求め、−9
    0度≦θ2≦θ0であれば、地磁気Mz符号はマイナス、
    θ0<θ2≦90度であれば、地磁気Mz符号はプラスと
    判別することを特徴とする姿勢角検出装置。
  7. 【請求項7】 請求項1乃至6の姿勢角検出装置におい
    て、前記判別装置により、2軸の加速度センサおよび、
    2軸の地磁気センサからなる静止角検出手段から求めら
    れる仮のヨー角Φ・仮のピッチ角P・仮のロール角Rの
    検出結果が、正しいかどうかの判別手段として、加速度
    センサの現在の出力Ax(n)から演算して得られた仮
    のロール角をR(n)、演算する単位時間で1単位時間
    過去の出力Ax(n−1)から得られた仮のロール角演
    算結果をR(n−1)とし、同様に加速度センサの現在
    の出力Ay(n)から演算して得られた仮のピッチ角を
    P(n)、演算する単位時間で1単位時間過去の出力A
    y(n−1)から得られた仮のピッチ角をP(n−1)
    とすると、数1あるいは数2のうち、少なくとも一つの
    値が、ある値以内にあるとき、静止角検出手段からの姿
    勢角情報は正しいと判別することを特徴とする姿勢角検
    出装置。 【数1】 【数2】
  8. 【請求項8】 請求項1乃至請求項7の姿勢角検出装置
    において、求めるべき出力をロール角α(n)、ピッチ
    角β(n)、ヨー角γ(n)、運動角検出手段から求め
    られる運動角を、仮のロール角変化△X(n)、仮のピ
    ッチ角変化△Y(n)、仮のヨー角変化△Z(n)、静
    止角検出手段から求められる静止角を、仮のロール角R
    (n)、仮のピッチ角P(n)、仮のヨー角Φ(n)、
    誤差より大きくない一定定数をC1,C2,C3とする
    と、判別装置の判別結果により、静止角検出手段からの
    姿勢角情報は正しいと判別する場合、α(n)=α(n
    −1)+△X(n)−C1、β(n)=β(n−1)+
    △Y(n)−C2、γ(n)=γ(n−1)+△Z
    (n)−C3により、また、静止角検出手段からの姿勢
    角情報は誤りと判別する場合は、α(n)=α(n−
    1)+△X(n)、β(n)=β(n−1)+△Y
    (n)、γ(n)=γ(n−1)+△Z(n)により、
    姿勢角を演算することを特徴とする姿勢角検出装置。
  9. 【請求項9】 請求項1乃至請求項8の姿勢角検出装置
    において、求めるべき出力をロール角α(n)、ピッチ
    角β(n)、ヨー角γ(n)、動的な運動角検出手段か
    ら求められる運動角を、仮のロール角変化△X(n)、
    仮のピッチ角変化△Y(n)、仮のヨー角変化△Z
    (n)、静的な静止角検出手段から求められる姿勢角
    を、仮のロール角R(n)、仮のピッチ角P(n)、仮
    のヨー角Φ(n)、1以下の比例定数をk1,k2,k3
    すると、判別装置の判別結果により、静止角検出手段か
    らの姿勢角情報は正しいと判別する場合は、α(n)=
    α(n−1)+△X(n)−k1[α(n−1)+△X
    (n)−R(n)]、β(n)=β(n−1)+△Y
    (n)−k2[β(n−1)+△Y(n)−P
    (n)]、γ(n)=γ(n−1)+△Z(n)−k3
    [γ(n−1)+△Z(n)−Φ(n)]により、ま
    た、静止角検出手段からの姿勢角情報は誤りと判別する
    場合は、α(n)=α(n−1)+△X(n)、β
    (n)=β(n−1)+△Y(n)、γ(n)=γ(n
    −1)+△Z(n)により、姿勢角を演算することを特
    徴とする姿勢角検出装置。
JP02374298A 1998-01-20 1998-01-20 姿勢角検出装置 Expired - Fee Related JP3763435B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP02374298A JP3763435B2 (ja) 1998-01-20 1998-01-20 姿勢角検出装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP02374298A JP3763435B2 (ja) 1998-01-20 1998-01-20 姿勢角検出装置

Publications (2)

Publication Number Publication Date
JPH11211474A true JPH11211474A (ja) 1999-08-06
JP3763435B2 JP3763435B2 (ja) 2006-04-05

Family

ID=12118773

Family Applications (1)

Application Number Title Priority Date Filing Date
JP02374298A Expired - Fee Related JP3763435B2 (ja) 1998-01-20 1998-01-20 姿勢角検出装置

Country Status (1)

Country Link
JP (1) JP3763435B2 (ja)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002164987A (ja) * 2000-11-29 2002-06-07 Tokin Corp 携帯端末機器用姿勢角度検出装置の取り付け方法
WO2004063667A1 (ja) * 2003-01-15 2004-07-29 National Institute Of Advanced Industrial Science And Technology 姿勢角処理装置および姿勢角処理方法
JP2007257198A (ja) * 2006-03-22 2007-10-04 Nintendo Co Ltd 傾き算出装置および傾き算出プログラムならびにゲーム装置およびゲームプログラム
JP2008517392A (ja) * 2004-10-22 2008-05-22 コミサリヤ・ア・レネルジ・アトミク 少なくとも3つの次元を有する空間においてナビゲートするための独立型のデバイス、システム、および方法
JP2009507272A (ja) * 2005-08-01 2009-02-19 トヨタ自動車株式会社 運動体の姿勢角検出装置
US7616232B2 (en) 2005-12-02 2009-11-10 Fujifilm Corporation Remote shooting system and camera system
KR100943538B1 (ko) * 2008-02-18 2010-02-22 (주)마이크로인피니티 요각 측정 장치 및 방법
JP2011059020A (ja) * 2009-09-11 2011-03-24 Chiba Inst Of Technology 三次元姿勢推定装置、三次元姿勢推定方法及び三次元姿勢推定プログラム
CN102759994A (zh) * 2012-06-07 2012-10-31 歌尔声学股份有限公司 一种控制***和一种控制设备
JP2015169582A (ja) * 2014-03-07 2015-09-28 リコーエレメックス株式会社 姿勢検出装置
JP2017090433A (ja) * 2015-11-10 2017-05-25 愛知製鋼株式会社 ボール回転方向検出システム
CN107085246A (zh) * 2017-05-11 2017-08-22 深圳合优科技有限公司 一种基于mems的人体动作识别方法和装置
CN112987944A (zh) * 2021-03-12 2021-06-18 东莞市攀星智能运动科技有限公司 基于脚部姿势交互的音乐播放控制方法
US11402927B2 (en) 2004-05-28 2022-08-02 UltimatePointer, L.L.C. Pointing device
US11841997B2 (en) 2005-07-13 2023-12-12 UltimatePointer, L.L.C. Apparatus for controlling contents of a computer-generated image using 3D measurements

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102901977B (zh) * 2012-10-24 2014-07-02 北京航天自动控制研究所 一种飞行器的初始姿态角的确定方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0821732A (ja) * 1994-07-05 1996-01-23 Data Tec:Kk 姿勢方位位置計測装置
JPH09106322A (ja) * 1995-10-09 1997-04-22 Data Tec:Kk ヘッドマウントディスプレイにおける姿勢角検出装置
JPH09318382A (ja) * 1996-05-27 1997-12-12 Nec Corp 姿勢角検出装置及び方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0821732A (ja) * 1994-07-05 1996-01-23 Data Tec:Kk 姿勢方位位置計測装置
JPH09106322A (ja) * 1995-10-09 1997-04-22 Data Tec:Kk ヘッドマウントディスプレイにおける姿勢角検出装置
JPH09318382A (ja) * 1996-05-27 1997-12-12 Nec Corp 姿勢角検出装置及び方法

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002164987A (ja) * 2000-11-29 2002-06-07 Tokin Corp 携帯端末機器用姿勢角度検出装置の取り付け方法
WO2004063667A1 (ja) * 2003-01-15 2004-07-29 National Institute Of Advanced Industrial Science And Technology 姿勢角処理装置および姿勢角処理方法
US7222045B2 (en) 2003-01-15 2007-05-22 National Institute Of Advanced Industrial Science And Technology Attitude angle processor and attitude angle processing method
US11755127B2 (en) 2004-05-28 2023-09-12 UltimatePointer, L.L.C. Multi-sensor device with an accelerometer for enabling user interaction through sound or image
US11416084B2 (en) 2004-05-28 2022-08-16 UltimatePointer, L.L.C. Multi-sensor device with an accelerometer for enabling user interaction through sound or image
US11409376B2 (en) 2004-05-28 2022-08-09 UltimatePointer, L.L.C. Multi-sensor device with an accelerometer for enabling user interaction through sound or image
US11402927B2 (en) 2004-05-28 2022-08-02 UltimatePointer, L.L.C. Pointing device
US8072417B2 (en) 2004-10-22 2011-12-06 Commissariat A L'energie Atomique Stand-alone device, system and method for navigating in a space having at least three dimensions
JP2008517392A (ja) * 2004-10-22 2008-05-22 コミサリヤ・ア・レネルジ・アトミク 少なくとも3つの次元を有する空間においてナビゲートするための独立型のデバイス、システム、および方法
US11841997B2 (en) 2005-07-13 2023-12-12 UltimatePointer, L.L.C. Apparatus for controlling contents of a computer-generated image using 3D measurements
JP2009507272A (ja) * 2005-08-01 2009-02-19 トヨタ自動車株式会社 運動体の姿勢角検出装置
JP4751930B2 (ja) * 2005-08-01 2011-08-17 トヨタ自動車株式会社 運動体の姿勢角検出装置
US7616232B2 (en) 2005-12-02 2009-11-10 Fujifilm Corporation Remote shooting system and camera system
US7920985B2 (en) 2006-03-22 2011-04-05 Nintendo Co., Ltd. Inclination calculation apparatus and inclination calculation program, and game apparatus and game program
JP4547346B2 (ja) * 2006-03-22 2010-09-22 任天堂株式会社 傾き算出装置および傾き算出プログラムならびにゲーム装置およびゲームプログラム
JP2007257198A (ja) * 2006-03-22 2007-10-04 Nintendo Co Ltd 傾き算出装置および傾き算出プログラムならびにゲーム装置およびゲームプログラム
KR100943538B1 (ko) * 2008-02-18 2010-02-22 (주)마이크로인피니티 요각 측정 장치 및 방법
JP2011059020A (ja) * 2009-09-11 2011-03-24 Chiba Inst Of Technology 三次元姿勢推定装置、三次元姿勢推定方法及び三次元姿勢推定プログラム
CN102759994A (zh) * 2012-06-07 2012-10-31 歌尔声学股份有限公司 一种控制***和一种控制设备
JP2015169582A (ja) * 2014-03-07 2015-09-28 リコーエレメックス株式会社 姿勢検出装置
JP2017090433A (ja) * 2015-11-10 2017-05-25 愛知製鋼株式会社 ボール回転方向検出システム
US10591637B2 (en) 2015-11-10 2020-03-17 Aichi Steel Corporation Ball rotational direction detecting system
CN107085246A (zh) * 2017-05-11 2017-08-22 深圳合优科技有限公司 一种基于mems的人体动作识别方法和装置
CN112987944A (zh) * 2021-03-12 2021-06-18 东莞市攀星智能运动科技有限公司 基于脚部姿势交互的音乐播放控制方法

Also Published As

Publication number Publication date
JP3763435B2 (ja) 2006-04-05

Similar Documents

Publication Publication Date Title
JP3842316B2 (ja) 位置検出装置および傾斜センサ
KR100533106B1 (ko) 지자계 센서의 자세 오차 보상장치 및 방법
US8552978B2 (en) 3D pointing device and method for compensating rotations of the 3D pointing device thereof
JP3763435B2 (ja) 姿勢角検出装置
US11698687B2 (en) Electronic device for use in motion detection and method for obtaining resultant deviation thereof
JPH09106322A (ja) ヘッドマウントディスプレイにおける姿勢角検出装置
GB1579920A (en) Gyroscope apparatus
JP3960574B2 (ja) 姿勢角検出装置
JP2001208529A (ja) 計測装置及びその制御方法並びに記憶媒体
JP3992206B2 (ja) 姿勢角検出装置
JP4023889B2 (ja) 姿勢角検出装置
JP3797661B2 (ja) 姿勢角度検出装置
JP2006275523A (ja) 電子方位装置および記録媒体
JP4004165B2 (ja) 姿勢角検出装置
JP5457890B2 (ja) 方位検知装置
JP3886005B2 (ja) 姿勢角度検出装置
JP3960575B2 (ja) 姿勢角検出装置
JP2000356520A (ja) 姿勢角検出装置
JP2005331265A (ja) 姿勢角度検出装置
JP3696586B2 (ja) 姿勢角度検出装置
JPH07175583A (ja) 指示入力装置
JP2002164987A (ja) 携帯端末機器用姿勢角度検出装置の取り付け方法
JPH0875442A (ja) 簡易型測長機
RU2130588C1 (ru) Способ измерения магнитного курса подвижного объекта
JP2001033253A (ja) 姿勢角検出装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20040420

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20050930

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20051011

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20051207

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20060110

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20060112

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100127

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110127

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110127

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120127

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130127

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140127

Year of fee payment: 8

LAPS Cancellation because of no payment of annual fees