[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

JP5602070B2 - POSITIONING DEVICE, POSITIONING METHOD OF POSITIONING DEVICE, AND POSITIONING PROGRAM - Google Patents

POSITIONING DEVICE, POSITIONING METHOD OF POSITIONING DEVICE, AND POSITIONING PROGRAM Download PDF

Info

Publication number
JP5602070B2
JP5602070B2 JP2011056105A JP2011056105A JP5602070B2 JP 5602070 B2 JP5602070 B2 JP 5602070B2 JP 2011056105 A JP2011056105 A JP 2011056105A JP 2011056105 A JP2011056105 A JP 2011056105A JP 5602070 B2 JP5602070 B2 JP 5602070B2
Authority
JP
Japan
Prior art keywords
calculated
navigation
azimuth
value
dimensional
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.)
Active
Application number
JP2011056105A
Other languages
Japanese (ja)
Other versions
JP2012193965A (en
Inventor
尚幸 梶原
真和 小崎
嘉宏 島
純一 瀧口
隆二郎 黒崎
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.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric 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 Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP2011056105A priority Critical patent/JP5602070B2/en
Publication of JP2012193965A publication Critical patent/JP2012193965A/en
Application granted granted Critical
Publication of JP5602070B2 publication Critical patent/JP5602070B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Navigation (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Description

本発明は、例えば、車両の位置を標定する位置標定装置、位置標定装置の位置標定方法および位置標定プログラムに関するものである。   The present invention relates to, for example, a position locating device for locating the position of a vehicle, a position locating method for the position locating device, and a position locating program.

従来、6自由度の3次元航法は、高価で高精度なIMU(Inertial Measurement Unit)を利用する必要があったため、航空・宇宙関連分野など、限られた分野でしか使用できなかった。   Conventionally, 6-degree-of-freedom three-dimensional navigation has been required to use an expensive and highly accurate IMU (Internal Measurement Unit), and can only be used in limited fields such as aerospace.

そこで、低価格化のために複数台のGPS(Global Positioning System)受信機と車速パルスとを利用して車両の位置・姿勢・方位を標定することも考えられる。
しかし、GPSアンテナ同士をある程度離して相対位置が変化しないように設置する必要がある。この制約は、一般の乗用車に適用するには大きな制約である。
In view of this, it is conceivable to determine the position / orientation / orientation of the vehicle using a plurality of GPS (Global Positioning System) receivers and vehicle speed pulses in order to reduce the price.
However, it is necessary to install the GPS antennas so that the relative position does not change with some distance. This restriction is a great restriction when applied to general passenger cars.

これらの制約が発生する具体的な要因は以下のとおりである。
(1)車両を停止した状態で方位を求める場合、中精度または低精度なIMUでは、地球の自転を計測して方位を求めるといった計算が出来ない。
また、2台以上のGPS受信機を使用して方位を求める場合、GPSアンテナ間の基線解析を行うため、精度を出すためにはGPSアンテナを可能な限り離して相対位置が変化しないよう設置する必要がある。
(2)3次元航法計算では、ある程度正確な初期値(位置・姿勢・方位)を設定しなければ計算結果が発散してしまう可能性が高い。
初期値のうち、姿勢角(ピッチ角、ロール角)は重力観測により得られ、位置はGPS観測により得られるが、方位に関しては上記(1)の理由により適切な値を得ることができない。
(3)正確な初期値が何らかの方法で与えられたとしても、その後の停車時にはジャイロの誤差等によって方位の値がドリフトしてしまう。
The specific factors that cause these restrictions are as follows.
(1) When obtaining the azimuth while the vehicle is stopped, the medium accuracy or low accuracy IMU cannot calculate the azimuth by measuring the rotation of the earth.
Also, when finding the direction using two or more GPS receivers, in order to perform baseline analysis between the GPS antennas, install the GPS antennas as far apart as possible so that the relative position does not change in order to obtain accuracy. There is a need.
(2) In the three-dimensional navigation calculation, the calculation result is likely to diverge unless initial values (position / posture / orientation) that are accurate to some extent are set.
Of the initial values, the posture angle (pitch angle, roll angle) is obtained by gravity observation, and the position is obtained by GPS observation, but an appropriate value cannot be obtained for the direction for the reason (1) above.
(3) Even if an accurate initial value is given by any method, the direction value drifts due to a gyro error or the like when the vehicle stops thereafter.

今後、3次元地図のインフラ整備が急速に進むほか、低価格なMEMS(Micro Electro Mechanical Systems)センサの性能向上によって一般のカーナビゲーションシステムでも車両の位置・姿勢・方位を標定する3次元航法のニーズが高まることが予想されている。   In the future, in addition to the rapid development of 3D map infrastructure, there is a need for 3D navigation to determine the position, orientation, and orientation of vehicles in general car navigation systems by improving the performance of low-cost MEMS (Micro Electro Mechanical Systems) sensors. Is expected to increase.

特開2010−287156号公報JP 2010-287156 A

本発明は、例えば、高精度でないIMUと1台のGPS受信機とを用いて車両の位置を高精度に標定できるようにすることを目的とする。   An object of the present invention is, for example, to be able to determine the position of a vehicle with high accuracy using an IMU that is not highly accurate and a single GPS receiver.

本発明の位置標定装置は、GPS(Global Positioning System)を利用して測位を行うGPS受信機と、速度を計測する速度センサと、加速度を計測する加速度センサと、方位角の角速度を計測する角速度センサとが取り付けられた移動体の位置を標定する。
前記位置標定装置は、
前記角速度センサにより計測された方位角の角速度に基づいて前記移動体の方位角を第一航法方位角として算出し、算出した第一航法方位角と前記速度センサにより計測された速度とに基づいて前記移動体の座標値を第一航法座標値として算出し、算出した第一航法座標値と前記GPS受信機により測位された座標値とに基づいて第1のカルマンフィルタを用いて前記第一航法方位角を補正する補正値を算出し、算出した補正値を用いて前記第一航法方位角を補正し、補正後の第一航法方位角を方位角初期値として出力する初期値算出部と、
前記角速度センサにより計測された方位角の角速度に基づいて方位角の変化量を算出し、算出した方位角の変化量と前記初期値算出部から出力された方位角初期値とに基づいて前記移動体の方位角を第二航法方位角として算出し、算出した第二航法方位角と前記加速度センサにより計測された加速度とに基づいて前記移動体の座標値を第二航法座標値として算出し、算出した第二航法座標値と前記GPS受信機により測位された座標値とに基づいて第2のカルマンフィルタを用いて前記第二航法座標値を補正する補正値を算出し、算出した補正値を用いて前記第二航法座標値を補正し、前記移動体の位置を標定した位置標定値として補正後の第二航法座標値を出力する標定値算出部とを備える。
The position locating device of the present invention includes a GPS receiver that performs positioning using GPS (Global Positioning System), a speed sensor that measures speed, an acceleration sensor that measures acceleration, and an angular speed that measures the angular velocity of an azimuth angle. The position of the moving body to which the sensor is attached is determined.
The position locator is
Based on the angular velocity of the azimuth angle measured by the angular velocity sensor, the azimuth angle of the moving body is calculated as a first navigation azimuth angle, and based on the calculated first navigation azimuth angle and the velocity measured by the speed sensor. The coordinate value of the moving body is calculated as a first navigation coordinate value, and the first navigation direction is calculated using a first Kalman filter based on the calculated first navigation coordinate value and the coordinate value measured by the GPS receiver. An initial value calculation unit that calculates a correction value for correcting an angle, corrects the first navigation azimuth using the calculated correction value, and outputs the corrected first navigation azimuth as an azimuth initial value;
The amount of change in azimuth is calculated based on the angular velocity of the azimuth measured by the angular velocity sensor, and the movement is performed based on the calculated amount of change in azimuth and the initial value of azimuth output from the initial value calculator. Calculating the azimuth of the body as the second navigation azimuth, and calculating the coordinate value of the moving body as the second navigation coordinate value based on the calculated second navigation azimuth and the acceleration measured by the acceleration sensor, A correction value for correcting the second navigation coordinate value is calculated using a second Kalman filter based on the calculated second navigation coordinate value and the coordinate value measured by the GPS receiver, and the calculated correction value is used. An orientation value calculation unit that corrects the second navigation coordinate value and outputs the corrected second navigation coordinate value as a location orientation value obtained by locating the position of the moving body.

本発明によれば、例えば、高精度でないIMU(加速度センサ、角速度センサ)と1台のGPS受信機とを用いて車両(移動体の一例)の位置を高精度に標定することができる。   According to the present invention, for example, the position of a vehicle (an example of a moving body) can be determined with high accuracy using an IMU (acceleration sensor, angular velocity sensor) that is not highly accurate and a single GPS receiver.

実施の形態1における車両200の概要図。1 is a schematic diagram of a vehicle 200 according to Embodiment 1. FIG. 実施の形態1における自己位置姿勢標定装置100の概要図。1 is a schematic diagram of a self-position / posture locating apparatus 100 according to Embodiment 1. FIG. 実施の形態1における自己位置姿勢標定装置100の機能構成図。FIG. 3 is a functional configuration diagram of the self-position / posture locating apparatus 100 according to the first embodiment. 実施の形態1におけるストラップダウン演算部140の機能構成図。FIG. 3 is a functional configuration diagram of a strapdown calculation unit 140 in the first embodiment. 実施の形態1における自己位置姿勢標定方法を示すフローチャート。5 is a flowchart illustrating a self-position / posture locating method according to the first embodiment. 実施の形態1における3自由度2次元航法のアルゴリズムを示す図。FIG. 3 shows an algorithm for three-degree-of-freedom two-dimensional navigation in the first embodiment. 実施の形態1における3自由度2次元航法による方位角の推定結果を示すグラフ。3 is a graph showing an azimuth angle estimation result by three-degree-of-freedom two-dimensional navigation in the first embodiment. 実施の形態1における6自由度3次元航法のアルゴリズムを示す図。FIG. 5 shows an algorithm for six-degree-of-freedom three-dimensional navigation in the first embodiment. 実施の形態1における車両200の停車時の方位角のドリフトを示すグラフ。3 is a graph showing azimuth drift when vehicle 200 is stopped in the first embodiment. 実施の形態1における拡張ZUPTアルゴリズムの概要図。FIG. 3 is a schematic diagram of an extended ZUPT algorithm in the first embodiment. 実施の形態1における自己位置姿勢標定装置100のハードウェア資源の一例を示す図。FIG. 3 is a diagram illustrating an example of hardware resources of the self-position / posture locating apparatus 100 according to the first embodiment.

実施の形態1.
車両の位置・姿勢・方位を標定する自己位置姿勢標定装置について説明する。
Embodiment 1 FIG.
A self-position / posture locating apparatus for locating the position / posture / orientation of a vehicle will be described.

図1は、実施の形態1における車両200の概要図である。
図2は、実施の形態1における自己位置姿勢標定装置100の概要図である。
実施の形態1における車両200の概要と自己位置姿勢標定装置100の概要とについて、図1と図2とに基づいて説明する。
FIG. 1 is a schematic diagram of a vehicle 200 in the first embodiment.
FIG. 2 is a schematic diagram of the self-position posture locating device 100 according to the first embodiment.
The outline of the vehicle 200 and the outline of the self-position / posture locating apparatus 100 in the first embodiment will be described based on FIG. 1 and FIG.

図1において、車両200(移動体の一例)は、GPS受信機210(およびGPSアンテナ211)と慣性センサ220、車速センサ230(速度センサの一例)と自己位置姿勢標定装置100(位置標定装置の一例)を備える。
「GPS」は「Global Positioning System」の略である。
In FIG. 1, a vehicle 200 (an example of a moving body) includes a GPS receiver 210 (and a GPS antenna 211), an inertial sensor 220, a vehicle speed sensor 230 (an example of a speed sensor), and a self-position / posture locator 100 (a position locator). An example).
“GPS” is an abbreviation for “Global Positioning System”.

例えば、GPS受信機210と慣性センサ220とを設置した天板201を車両200の屋根部に設置する。また、自己位置姿勢標定装置100と車速センサ230とを車両200内に設置する。
但し、設置方法はこれに限らない。例えば、GPS受信機210を屋根部以外に設置してもよいし、慣性センサ220を車両200内に設置してもよく、天板201を設けることなく、GPS受信機210と慣性センサ220を車両200の内部または外部に取り付けてもよい。
また、自己位置姿勢標定装置100は車両200とは別個に設けた外部装置であっても構わない。
さらに、車両200は他の移動体(例えば、飛行機、船舶、自動二輪車、自転車、列車)であっても構わない。
For example, the top plate 201 on which the GPS receiver 210 and the inertial sensor 220 are installed is installed on the roof portion of the vehicle 200. In addition, the self-position / posture locator 100 and the vehicle speed sensor 230 are installed in the vehicle 200.
However, the installation method is not limited to this. For example, the GPS receiver 210 may be installed outside the roof, or the inertial sensor 220 may be installed in the vehicle 200. The GPS receiver 210 and the inertial sensor 220 may be installed in the vehicle without providing the top plate 201. It may be attached inside or outside the 200.
Further, the self-position / posture locating device 100 may be an external device provided separately from the vehicle 200.
Furthermore, the vehicle 200 may be another moving body (for example, an airplane, a ship, a motorcycle, a bicycle, or a train).

GPS受信機210は、GPS衛星から発信される測位信号をGPSアンテナ211を用いて観測し、観測結果に基づいて位置を測位する。   The GPS receiver 210 observes a positioning signal transmitted from a GPS satellite using the GPS antenna 211, and measures the position based on the observation result.

慣性センサ220は、3軸方向(車両200の長さ方向と幅方向と高さ方向)の加速度と角速度とを計測する。
慣性センサ220は、IMU(Inertial Measurement Unit)と呼ばれる。
以下、3軸方向の加速度を「三次元加速度」といい、3軸方向の角速度を「三次元角速度」という。
Inertial sensor 220 measures acceleration and angular velocity in three axis directions (the length direction, the width direction, and the height direction of vehicle 200).
The inertial sensor 220 is called an IMU (Internal Measurement Unit).
Hereinafter, the acceleration in the three-axis direction is referred to as “three-dimensional acceleration”, and the angular velocity in the three-axis direction is referred to as “three-dimensional angular velocity”.

車速センサ230は、車軸の回転を車速パルスとして検出し、検出した車速パルスを用いて車両200の速度を計測する。   The vehicle speed sensor 230 detects the rotation of the axle as a vehicle speed pulse, and measures the speed of the vehicle 200 using the detected vehicle speed pulse.

図2に示すように、自己位置姿勢標定装置100は、GPS受信機210のデータと、慣性センサ220のデータと、車速センサ230のデータとに基づいて車両200の位置と姿勢と方位とを標定する。
「位置」は、三次元の座標値(緯度、経度、高度)を意味する。
「姿勢」は、仰角(ピッチ角)と回転角(ロール角)とを意味する。
「方位」は、方位角を意味する。
As shown in FIG. 2, the self-position / posture locating device 100 determines the position, posture, and direction of the vehicle 200 based on the data of the GPS receiver 210, the data of the inertial sensor 220, and the data of the vehicle speed sensor 230. To do.
“Position” means a three-dimensional coordinate value (latitude, longitude, altitude).
“Position” means an elevation angle (pitch angle) and a rotation angle (roll angle).
“Direction” means an azimuth angle.

図3は、実施の形態1における自己位置姿勢標定装置100の機能構成図である。
実施の形態1における自己位置姿勢標定装置100の機能構成について、図3に基づいて説明する。
FIG. 3 is a functional configuration diagram of the self-position posture locating device 100 according to the first embodiment.
A functional configuration of self-position / posture locating apparatus 100 according to the first embodiment will be described with reference to FIG.

自己位置姿勢標定装置100(位置標定装置の一例)は、初期値算出部110、自己位置姿勢標定部120、カルマンフィルタ部130、ストラップダウン演算部140および装置記憶部190を備える。   The self-position / posture locating apparatus 100 (an example of the position locating apparatus) includes an initial value calculating unit 110, a self-position / posture locating unit 120, a Kalman filter unit 130, a strap-down calculating unit 140, and a device storage unit 190.

装置記憶部190は、自己位置姿勢標定装置100で使用するデータを記憶する。
GPSデータ191、速度データ192、加速度データ193または角速度データ194は、装置記憶部190に記憶されるデータの一例である。
The device storage unit 190 stores data used by the self-position / posture orientation device 100.
The GPS data 191, speed data 192, acceleration data 193, or angular velocity data 194 is an example of data stored in the device storage unit 190.

GPSデータ191は、GPS受信機210により計測されたデータであり、時刻に対応付けて「位置」や「速度」を含んでいる。速度は、過去の位置と現在の位置と過去から現在までの経過時間とに基づいて計測される。
速度データ192は、車速センサ230により計測されたデータであり、時刻に対応付けて「速度」を含んでいる。
加速度データ193は、慣性センサ220を構成する加速度センサ221により計測されたデータであり、時刻に対応付けて三次元加速度を含んでいる。
角速度データ194は、慣性センサ220を構成する角速度センサ222により計測されたデータであり、時刻に対応付けて三次元角速度を含んでいる。つまり、角速度データ194は、「仰角」の角速度と「回転角」の角速度と「方位角」の角速度とを含んでいる。角速度センサ222は「ジャイロ」と呼ばれる。
The GPS data 191 is data measured by the GPS receiver 210 and includes “position” and “speed” in association with time. The speed is measured based on the past position, the current position, and the elapsed time from the past to the present.
The speed data 192 is data measured by the vehicle speed sensor 230 and includes “speed” in association with the time.
The acceleration data 193 is data measured by the acceleration sensor 221 constituting the inertial sensor 220, and includes three-dimensional acceleration in association with time.
The angular velocity data 194 is data measured by the angular velocity sensor 222 that constitutes the inertial sensor 220, and includes a three-dimensional angular velocity in association with time. That is, the angular velocity data 194 includes an angular velocity of “elevation angle”, an angular velocity of “rotation angle”, and an angular velocity of “azimuth angle”. The angular velocity sensor 222 is called a “gyro”.

初期値算出部110は、カルマンフィルタ部130を制御して方位角の初期値を算出する(初期値算出処理)。   The initial value calculation unit 110 controls the Kalman filter unit 130 to calculate the initial value of the azimuth (initial value calculation process).

例えば、初期値算出部110は、以下のように方位角の初期値を算出する。
初期値算出部110は、角速度データ194に含まれる方位角の角速度に基づいて、車両200の方位角を第一航法方位角として算出する(後述する方位更新部112の処理)。
初期値算出部110は、算出した第一航法方位角と速度データ192に含まれる速度とに基づいて、車両200の座標値を第一航法座標値として算出する(後述する緯度経度更新部114の処理)。
初期値算出部110は、算出した第一航法座標値とGPSデータ191に含まれる座標値とに基づいて、第1のカルマンフィルタ131を用いて、第一航法方位角を補正する補正値を算出する。
初期値算出部110は、算出した補正値を用いて第一航法方位角を補正し、補正後の第一航法方位角を方位角初期値として出力する(後述する出力補正部115の処理)。
For example, the initial value calculation unit 110 calculates the initial value of the azimuth as follows.
The initial value calculation unit 110 calculates the azimuth angle of the vehicle 200 as the first navigation azimuth angle based on the angular velocity of the azimuth angle included in the angular velocity data 194 (processing of the azimuth update unit 112 described later).
The initial value calculation unit 110 calculates the coordinate value of the vehicle 200 as the first navigation coordinate value based on the calculated first navigation azimuth angle and the speed included in the speed data 192 (the latitude and longitude update unit 114 described later). processing).
The initial value calculation unit 110 calculates a correction value for correcting the first navigation azimuth using the first Kalman filter 131 based on the calculated first navigation coordinate value and the coordinate value included in the GPS data 191. .
The initial value calculation unit 110 corrects the first navigation azimuth using the calculated correction value, and outputs the corrected first navigation azimuth as an azimuth initial value (processing of the output correction unit 115 described later).

自己位置姿勢標定部120は、走行時標定部121(標定値算出部の一例)と停車時標定部122(方位角補正部の一例)と停車判定部123(停止判定部の一例)とを備える。
走行時標定部121は、車両200が走行しているときに、カルマンフィルタ部130とストラップダウン演算部140とを制御して車両200の位置と姿勢と方位とを標定(算出)する。
停車時標定部122は、車両200が停車しているときに、カルマンフィルタ部130とストラップダウン演算部140とを制御して車両200の位置と姿勢と方位とを標定する。
The self-position / posture orientation unit 120 includes a travel time location unit 121 (an example of a location value calculation unit), a stop time location unit 122 (an example of an azimuth angle correction unit), and a stop determination unit 123 (an example of a stop determination unit). .
The running time orientation unit 121 controls (calculates) the position, orientation, and orientation of the vehicle 200 by controlling the Kalman filter unit 130 and the strapdown calculation unit 140 when the vehicle 200 is running.
When the vehicle 200 is stopped, the stop time orientation unit 122 controls the Kalman filter unit 130 and the strap-down calculation unit 140 to determine the position, orientation, and orientation of the vehicle 200.

例えば、走行時標定部121は、以下のように車両200の位置を標定する。
走行時標定部121は、角速度データ194に含まれる方位角の角速度に基づいて、方位角の変化量を算出する(後述する姿勢方位更新部142の処理)。
走行時標定部121は、算出した方位角の変化量と初期値算出部110から出力された方位角初期値とに基づいて、車両200の方位角を第二航法方位角として算出する(後述する姿勢方位更新部142の処理)。
走行時標定部121は、算出した第二航法方位角と加速度データ193に含まれる加速度とに基づいて、車両200の座標値を第二航法座標値として算出する(後述する速度位置更新部144の処理)。
走行時標定部121は、算出した第二航法座標値とGPSデータ191に含まれる座標値とに基づいて、第2のカルマンフィルタ132を用いて、第二航法座標値を補正する補正値を算出する。
走行時標定部121は、算出した補正値を用いて第二航法座標値を補正し、車両200の位置を標定した位置標定値として補正後の第二航法座標値を出力する(後述する出力補正部125の処理)。
For example, the travel time orientation unit 121 locates the position of the vehicle 200 as follows.
The travel time orientation unit 121 calculates the amount of change in the azimuth based on the angular velocity of the azimuth included in the angular velocity data 194 (processing of the orientation / azimuth updating unit 142 described later).
The travel time orientation unit 121 calculates the azimuth angle of the vehicle 200 as the second navigation azimuth angle based on the calculated azimuth angle change amount and the initial azimuth angle value output from the initial value calculation unit 110 (described later). Processing of posture orientation update unit 142).
Based on the calculated second navigation azimuth angle and the acceleration included in the acceleration data 193, the traveling time locating unit 121 calculates the coordinate value of the vehicle 200 as a second navigation coordinate value (of a speed position update unit 144 described later). processing).
Based on the calculated second navigation coordinate value and the coordinate value included in the GPS data 191, the running time orientation unit 121 calculates a correction value for correcting the second navigation coordinate value using the second Kalman filter 132. .
The travel time location unit 121 corrects the second navigation coordinate value using the calculated correction value, and outputs the corrected second navigation coordinate value as a location value obtained by locating the position of the vehicle 200 (output correction described later). Processing of the unit 125).

例えば、走行時標定部121と停車時標定部122とは、以下のように車両200の位置を標定する。
停車時標定部122は、移動していた車両200が移動を停止した場合、第3のカルマンフィルタ133を用いて第二航法方位角を補正する方位角補正値を算出する。
停車時標定部122は、算出した方位角補正値を用いて第二航法方位角を補正して補正後の第二航法方位角を算出する(後述する出力補正部125の処理)。
走行時標定部121は、停止した車両200が移動を再開した場合、補正後の第二航法方位角と新たな方位角の変化量とに基づいて新たな第二航法方位角を算出する(後述する姿勢方位更新部142の処理)。
走行時標定部121は、算出した新たな第二航法方位角に基づいて新たな第二航法座標値を算出する(後述する速度位置更新部144の処理)。
走行時標定部121は、算出した新たな第二航法座標値を補正する位置補正値を第2のカルマンフィルタ132を用いて算出する。
走行時標定部121は、算出した位置補正値を用いて新たな第二航法座標値を補正し、補正後の新たな第二航法座標値を新たな位置標定値として出力する(後述する出力補正部125の処理)。
For example, the travel time location unit 121 and the stop time location unit 122 determine the position of the vehicle 200 as follows.
When the moving vehicle 200 stops moving, the stop time locating unit 122 calculates an azimuth correction value for correcting the second navigation azimuth using the third Kalman filter 133.
The stopping time locating unit 122 corrects the second navigation azimuth by using the calculated azimuth correction value and calculates a corrected second navigation azimuth (processing of the output correction unit 125 described later).
When the stopped vehicle 200 resumes movement, the running time orientation unit 121 calculates a new second navigation azimuth based on the corrected second navigation azimuth and the amount of change in the new azimuth (described later). Processing of posture orientation updating unit 142).
The traveling time locating unit 121 calculates a new second navigation coordinate value based on the calculated new second navigation azimuth angle (processing of a speed position update unit 144 described later).
The running time orientation unit 121 calculates a position correction value for correcting the calculated new second navigation coordinate value using the second Kalman filter 132.
The travel time location unit 121 corrects the new second navigation coordinate value using the calculated position correction value, and outputs the new second navigation coordinate value after correction as a new position location value (output correction described later). Processing of the unit 125).

停車判定部123は、車両200が停車したか又は発車したかを判定する。   The stop determination unit 123 determines whether the vehicle 200 has stopped or started.

例えば、停車判定部123は、車両200が停車したか又は発車したかを以下のように判定する。
停車判定部123は、車速センサ230により計測された速度が所定の速度閾値未満である場合、車両200が移動を停止したと判定する。
停車判定部123は、車両200が移動を停止したと判定した後、車速センサ230により計測された速度が速度閾値以上である場合、車両200が移動を再開したと判定する。
For example, the stop determination unit 123 determines whether the vehicle 200 has stopped or started as follows.
The stop determination unit 123 determines that the vehicle 200 has stopped moving when the speed measured by the vehicle speed sensor 230 is less than a predetermined speed threshold.
After determining that the vehicle 200 has stopped moving, the stop determination unit 123 determines that the vehicle 200 has resumed movement when the speed measured by the vehicle speed sensor 230 is equal to or greater than the speed threshold.

カルマンフィルタ部130は、観測方程式や状態方程式を生成し、生成した観測方程式や状態方程式を用いて各種データの補正値を算出する。
カルマンフィルタ部130は、観測方程式が異なる複数のカルマンフィルタ(符号131−133)を備える。
The Kalman filter unit 130 generates observation equations and state equations, and calculates correction values for various data using the generated observation equations and state equations.
The Kalman filter unit 130 includes a plurality of Kalman filters (reference numerals 131 to 133) having different observation equations.

例えば、第3のカルマンフィルタ133は、前回算出された第二航法方位角と今回算出された第二航法方位角とを一致させる補正値を方位角補正値として算出する。   For example, the third Kalman filter 133 calculates a correction value that matches the previously calculated second navigation azimuth with the currently calculated second navigation azimuth as the azimuth correction value.

例えば、第3のカルマンフィルタ133は、前回算出された第二航法方位角と今回算出された第二航法方位角との差分を表す観測方程式を生成し、生成した観測方程式を用いて方位角補正値を算出する。   For example, the third Kalman filter 133 generates an observation equation representing a difference between the previously calculated second navigation azimuth angle and the currently calculated second navigation azimuth angle, and uses the generated observation equation to calculate an azimuth correction value. Is calculated.

図4は、実施の形態1におけるストラップダウン演算部140の機能構成図である。
実施の形態1におけるストラップダウン演算部140の機能構成について、図4に基づいて説明する。
FIG. 4 is a functional configuration diagram of the strapdown calculation unit 140 in the first embodiment.
The functional configuration of the strapdown calculation unit 140 in the first embodiment will be described with reference to FIG.

ストラップダウン演算部140は、加速度データ193と角速度データ194とに基づいて速度と位置と姿勢と方位とを算出する。   The strapdown calculation unit 140 calculates a speed, a position, an attitude, and an orientation based on the acceleration data 193 and the angular velocity data 194.

ストラップダウン演算部140は、図4に示す構成(符号141−147、148A−C)を備える。   The strapdown calculation unit 140 has the configuration (reference numerals 141 to 147 and 148A-C) shown in FIG.

角速度座標変換部141は、センサ座標系(車体座標系)の三次元角速度を入力する。
角速度座標変換部141は、姿勢方位更新部142により前回算出された姿勢と方位とに基づいて所定の座標変換式(例えば、方向余弦行列)を用いて三次元角速度をセンサ座標系の値から航法座標系(局所水平座標系)の値に変換する。
「センサ座標系」は車両200を基準とする座標系であり、例えば、長さ方向を「x」、幅方向を「y」、高さ方向を「z」として表される。
「航法座標系」は地面を基準とする座標系であり、例えば、「North(北)」「East(東)」「Down(下)」で表される。
The angular velocity coordinate conversion unit 141 inputs a three-dimensional angular velocity of the sensor coordinate system (vehicle body coordinate system).
The angular velocity coordinate conversion unit 141 navigates the three-dimensional angular velocity from the value of the sensor coordinate system using a predetermined coordinate conversion formula (for example, a direction cosine matrix) based on the posture and the azimuth previously calculated by the posture azimuth updating unit 142. Convert to coordinate system (local horizontal coordinate system) value.
The “sensor coordinate system” is a coordinate system based on the vehicle 200, and is represented, for example, as “x” in the length direction, “y” in the width direction, and “z” in the height direction.
The “navigation coordinate system” is a coordinate system based on the ground, and is represented by, for example, “North (north)”, “East (east)”, and “Down (down)”.

第一補正部148Aは、角速度座標変換部141により変換された三次元角速度に航法系回転補正計算部147により算出される補正値を加算して三次元角速度を補正する。   The first correction unit 148A adds the correction value calculated by the navigation system rotation correction calculation unit 147 to the three-dimensional angular velocity converted by the angular velocity coordinate conversion unit 141 to correct the three-dimensional angular velocity.

姿勢方位更新部142は、第一補正部148Aにより補正された方位角の角速度を1回積分して方位角の変化量を算出し、算出した方位角の変化量を方位角の初期値または前回標定された方位角に加算して新たな方位角を算出する。
同様に、姿勢方位更新部142は、新たな姿勢角を算出する。
ここで、姿勢方位更新部142は、初期値算出部110により算出された方位角を方位角の初期値として用いる。
また、姿勢方位更新部142は、加速度センサ221により停車時に計測された加速度を用いて所定の姿勢角算出式を計算し、姿勢角の初期値を算出する。姿勢角算出式は、計測された加速度と重力加速度と姿勢角との関係を表した式である。
The attitude azimuth update unit 142 calculates the azimuth angle change amount by integrating the angular velocity of the azimuth angle corrected by the first correction unit 148A once, and calculates the calculated azimuth angle change amount as the initial value of the azimuth angle or the previous time. A new azimuth angle is calculated by adding to the standardized azimuth angle.
Similarly, the posture orientation update unit 142 calculates a new posture angle.
Here, the posture / azimuth update unit 142 uses the azimuth angle calculated by the initial value calculation unit 110 as the initial value of the azimuth angle.
In addition, the posture orientation updating unit 142 calculates a predetermined posture angle calculation formula using the acceleration measured by the acceleration sensor 221 when the vehicle is stopped, and calculates an initial value of the posture angle. The posture angle calculation formula is a formula that represents the relationship among measured acceleration, gravitational acceleration, and posture angle.

加速度座標変換部143は、センサ座標系の三次元加速度を入力する。
加速度座標変換部143は、姿勢方位更新部142により算出された姿勢角と方位角とに基づいて、所定の座標変換式を用いて、三次元加速度をセンサ座標系の値から航法座標系の値に変換する。
The acceleration coordinate conversion unit 143 inputs the three-dimensional acceleration of the sensor coordinate system.
The acceleration coordinate conversion unit 143 converts the three-dimensional acceleration from the value of the sensor coordinate system to the value of the navigation coordinate system using a predetermined coordinate conversion formula based on the attitude angle and the azimuth calculated by the attitude / azimuth update unit 142. Convert to

第二補正部148Bは、加速度座標変換部143により変換された三次元加速度に重力補正計算部145により算出される補正値を加算して三次元加速度を補正する。   The second correction unit 148B corrects the three-dimensional acceleration by adding the correction value calculated by the gravity correction calculation unit 145 to the three-dimensional acceleration converted by the acceleration coordinate conversion unit 143.

第三補正部148Cは、第二補正部148Bにより補正された三次元加速度にコリオリ補正計算部146により算出される補正値を加算して三次元加速度を補正する。   The third correction unit 148C corrects the three-dimensional acceleration by adding the correction value calculated by the Coriolis correction calculation unit 146 to the three-dimensional acceleration corrected by the second correction unit 148B.

速度位置更新部144は、第三補正部148Cにより補正された三次元加速度を1回積分して速度(速度ベクトル)を算出する。
また、速度位置更新部144は、第三補正部148Cにより補正された三次元加速度を2回積分して位置の変化量を算出し、算出した位置の変化量を位置の初期値または前回標定された位置に加算して新たな位置を算出する。
例えば、速度位置更新部144は、GPS受信機210により停車時に計測された位置を位置の初期値として用いる。
The speed position update unit 144 calculates the speed (speed vector) by integrating the three-dimensional acceleration corrected by the third correction unit 148C once.
Further, the speed position updating unit 144 calculates the amount of change in position by integrating the three-dimensional acceleration corrected by the third correction unit 148C twice, and the calculated amount of change in position is determined as the initial position value or the previous time. A new position is calculated by adding to the new position.
For example, the speed position update unit 144 uses the position measured by the GPS receiver 210 when the vehicle is stopped as the initial position value.

重力補正計算部145は、速度位置更新部144により算出された位置を入力値として所定の重力補正計算式を計算し、三次元加速度の補正値を算出する。この補正値は計測された三次元加速度に含まれる重力加速度を取り除くための値である。   The gravity correction calculation unit 145 calculates a predetermined gravity correction calculation formula using the position calculated by the velocity position update unit 144 as an input value, and calculates a correction value for three-dimensional acceleration. This correction value is a value for removing the gravitational acceleration included in the measured three-dimensional acceleration.

コリオリ補正計算部146は、速度位置更新部144により算出された速度と姿勢方位更新部142により算出された姿勢とを入力値として所定のコリオリ補正計算式を計算し、三次元加速度の補正値を算出する。この補正値はコリオリ力による計測誤差を三次元加速度から取り除くための値である。   The Coriolis correction calculation unit 146 calculates a predetermined Coriolis correction calculation formula using the velocity calculated by the velocity position update unit 144 and the posture calculated by the posture direction update unit 142 as input values, and calculates a correction value of the three-dimensional acceleration. calculate. This correction value is a value for removing the measurement error due to the Coriolis force from the three-dimensional acceleration.

航法系回転補正計算部147は、速度位置更新部144により算出された速度を入力値として所定の航法系回転補正計算式を計算し、三次元角速度の補正値を算出する。この補正値は地球の自転による計測誤差を三次元角速度から取り除くための値である。   The navigation system rotation correction calculation unit 147 calculates a predetermined navigation system rotation correction calculation formula using the speed calculated by the speed position update unit 144 as an input value, and calculates a three-dimensional angular velocity correction value. This correction value is a value for removing a measurement error due to the rotation of the earth from the three-dimensional angular velocity.

図5は、実施の形態1における自己位置姿勢標定方法を示すフローチャートである。
実施の形態1における自己位置姿勢標定方法について、図5に基づいて説明する。
FIG. 5 is a flowchart showing the self-position / posture determination method according to the first embodiment.
The self-position / posture locating method according to the first embodiment will be described with reference to FIG.

まず、自己位置姿勢標定方法の概要について説明する。   First, an outline of the self-position / posture locating method will be described.

初期値算出部110は、車両200の出発時に3自由度2次元航法により方位角の初期値を算出する(S110)。
走行時標定部121は、車両200が停車するまで、方位角の初期値または前回の方位角を用いて6自由度3次元航法により位置と姿勢と方位とを標定する(S120)。
停車時標定部122は、車両200が停車している間、拡張ZUPTアルゴリズムに基づいて6自由度3次元航法により位置と姿勢と方位とを標定する(S140)。
The initial value calculation unit 110 calculates the initial value of the azimuth by three-degree-of-freedom two-dimensional navigation when the vehicle 200 departs (S110).
The travel time orientation unit 121 locates the position, orientation, and orientation by 6-degree-of-freedom three-dimensional navigation using the initial value of the azimuth angle or the previous azimuth angle until the vehicle 200 stops (S120).
While the vehicle 200 is stopped, the stop time location unit 122 locates the position, orientation, and direction by 6-degree-of-freedom three-dimensional navigation based on the extended ZUPT algorithm (S140).

次に、自己位置姿勢標定方法の詳細について説明する。   Next, the details of the self-position orientation determination method will be described.

S110において、初期値算出部110は、車両200の出発時、つまり、自己位置姿勢標定装置100の起動時に3自由度2次元航法(第一航法)により方位角の初期値を算出する。
「3自由度」は「緯度」と「経度」と「方位」とを意味し、「2次元」は「緯度、経度」または「方位」を意味する。
In S110, the initial value calculation unit 110 calculates the initial value of the azimuth by three-degree-of-freedom two-dimensional navigation (first navigation) when the vehicle 200 departs, that is, when the self-position / posture locating device 100 is activated.
“3 degrees of freedom” means “latitude”, “longitude”, and “azimuth”, and “two-dimensional” means “latitude, longitude” or “azimuth”.

図6は、実施の形態1における3自由度2次元航法のアルゴリズムを示す図である。
実施の形態1における3自由度2次元航法について、図6に基づいて説明する。
FIG. 6 is a diagram showing an algorithm for three-degree-of-freedom two-dimensional navigation in the first embodiment.
The three-degree-of-freedom two-dimensional navigation in Embodiment 1 will be described with reference to FIG.

入力補正部111、方位更新部112、座標変換部113、緯度経度更新部114および出力補正部115は、初期値算出部110が備える機能構成である。   The input correction unit 111, the azimuth update unit 112, the coordinate conversion unit 113, the latitude / longitude update unit 114, and the output correction unit 115 are functional configurations included in the initial value calculation unit 110.

入力補正部111は、速度データ192に含まれる速度にカルマンフィルタ部130により算出される補正値を加算して速度(速度ベクトル)を補正する。
さらに、入力補正部111は、角速度データ194に含まれる「方位角」の角速度にカルマンフィルタ部130により算出される補正値を加算して方位角の角速度を補正する。
The input correction unit 111 corrects the speed (speed vector) by adding the correction value calculated by the Kalman filter unit 130 to the speed included in the speed data 192.
Further, the input correction unit 111 corrects the angular velocity of the azimuth by adding the correction value calculated by the Kalman filter unit 130 to the angular velocity of “azimuth” included in the angular velocity data 194.

方位更新部112は、入力補正部111により補正された方位角の角速度を1回積分して方位角の変化量を算出する。
方位更新部112は、算出した方位角の変化量を方位角の初期値または出力補正部115により前回算出された方位角に加算して新たな方位角を算出する。
ここで、方位更新部112は方位角の初期値として所定の値(例えば、0度)を用いる。ここで用いる方位角の初期値は適当な値で構わない。つまり、方位角の初期値は、正確な方位角を示す適切な値であっても、不正確な方位角を示す不適切な値であっても構わない。
The azimuth update unit 112 calculates the amount of change in azimuth by integrating the angular velocity of the azimuth corrected by the input correction unit 111 once.
The azimuth updating unit 112 calculates a new azimuth angle by adding the calculated azimuth angle change amount to the initial value of the azimuth angle or the azimuth angle previously calculated by the output correction unit 115.
Here, the azimuth updating unit 112 uses a predetermined value (for example, 0 degrees) as the initial value of the azimuth angle. The initial value of the azimuth angle used here may be an appropriate value. That is, the initial value of the azimuth angle may be an appropriate value indicating an accurate azimuth angle or an inappropriate value indicating an inaccurate azimuth angle.

座標変換部113は、方位更新部112により算出された姿勢と方位とに基づいて、所定の座標変換式を用いて、入力補正部111により補正された速度をセンサ座標系の値から航法座標系の値に変換する。
これにより、航法座標系の二次元速度(緯度方向の速度、経度方向の速度)が得られる。
The coordinate conversion unit 113 calculates the speed corrected by the input correction unit 111 from the value of the sensor coordinate system using a predetermined coordinate conversion formula based on the attitude and the direction calculated by the direction update unit 112 from the navigation coordinate system. Convert to the value of.
As a result, a two-dimensional velocity (velocity in the latitude direction, velocity in the longitude direction) of the navigation coordinate system is obtained.

緯度経度更新部114は、座標変換部113により得られた二次元速度を1回積分して二次元座標(緯度、経度)の変化量を算出する。
緯度経度更新部114は、算出した二次元座標の変化量を二次元座標の初期値または出力補正部115により前回算出された二次元座標に加算して新たな二次元座標を算出する。
例えば、緯度経度更新部114は、車両200が出発する前、つまり、車両200の停車時にGPS受信機210により測位された二次元座標値を二次元座標の初期値として用いる。
The latitude / longitude update unit 114 integrates the two-dimensional velocity obtained by the coordinate conversion unit 113 once to calculate a change amount of the two-dimensional coordinates (latitude, longitude).
The latitude / longitude update unit 114 calculates a new two-dimensional coordinate by adding the calculated change amount of the two-dimensional coordinate to the initial value of the two-dimensional coordinate or the two-dimensional coordinate previously calculated by the output correction unit 115.
For example, the latitude / longitude update unit 114 uses the two-dimensional coordinate value measured by the GPS receiver 210 before the vehicle 200 departs, that is, when the vehicle 200 stops, as the initial value of the two-dimensional coordinate.

カルマンフィルタ部130は、GPSデータ191に含まれる二次元座標と緯度経度更新部114により算出された二次元座標との差を「座標残差」として算出する。
カルマンフィルタ部130は、GPSデータ191に含まれる二次元速度と座標変換部113により得られた二次元速度とを差を「速度残差」として算出する。
カルマンフィルタ部130は、座標残差と速度残差とを入力値として第1のカルマンフィルタ131を実行し、速度の補正値と角速度の補正値と二次元座標の補正値と方位角の補正値とを算出する。
The Kalman filter unit 130 calculates a difference between the two-dimensional coordinates included in the GPS data 191 and the two-dimensional coordinates calculated by the latitude / longitude update unit 114 as a “coordinate residual”.
The Kalman filter unit 130 calculates a difference between the two-dimensional velocity included in the GPS data 191 and the two-dimensional velocity obtained by the coordinate conversion unit 113 as a “velocity residual”.
The Kalman filter unit 130 executes the first Kalman filter 131 using the coordinate residual and the velocity residual as input values, and obtains the velocity correction value, the angular velocity correction value, the two-dimensional coordinate correction value, and the azimuth angle correction value. calculate.

出力補正部115は、緯度経度更新部114により算出された二次元座標値にカルマンフィルタ部130により算出された補正値を加算して二次元座標値を補正し、補正した二次元座標値を出力する。
さらに、出力補正部115は、方位更新部112により算出された方位角にカルマンフィルタ部130により算出された補正値を加算して方位角を補正し、補正した方位角を出力する。
The output correction unit 115 adds the correction value calculated by the Kalman filter unit 130 to the two-dimensional coordinate value calculated by the latitude / longitude update unit 114 to correct the two-dimensional coordinate value, and outputs the corrected two-dimensional coordinate value. .
Further, the output correction unit 115 corrects the azimuth by adding the correction value calculated by the Kalman filter unit 130 to the azimuth calculated by the azimuth updating unit 112, and outputs the corrected azimuth.

図7は、実施の形態1における3自由度2次元航法による方位角の推定結果を示すグラフである。
実施の形態1における3自由度2次元航法による方位角の推定結果について、図7に基づいて説明する。
FIG. 7 is a graph showing the estimation result of the azimuth angle by the three-degree-of-freedom two-dimensional navigation in the first embodiment.
An estimation result of the azimuth angle by the three-degree-of-freedom two-dimensional navigation in the first embodiment will be described with reference to FIG.

図7において、3自由度2次元航法(図6参照)により算出した方位角を実線で示し、正しい方位角を鎖線で示す。   In FIG. 7, the azimuth angle calculated by three-degree-of-freedom two-dimensional navigation (see FIG. 6) is indicated by a solid line, and the correct azimuth angle is indicated by a chain line.

(1)車両200が「−180度」の方位を向いて停車しているときに方位角の初期値「0度」を初期値算出部110の方位更新部112に設定し、車両200の走行を開始した。
(2)車両200を左回りに回転させると、3自由度2次元航法により算出した方位角(実線)は数秒程度で正しい方位角(鎖線)に近い適切な値になった。
(3)車両200を停車させると、角速度センサ222のバイアス(計測誤差)により方位角の値がドリフトを始めた。
(1) The initial value “0 degree” of the azimuth angle is set in the azimuth update unit 112 of the initial value calculation unit 110 when the vehicle 200 is stopped in the direction of “−180 degrees”. Started.
(2) When the vehicle 200 is rotated counterclockwise, the azimuth angle (solid line) calculated by the three-degree-of-freedom two-dimensional navigation becomes an appropriate value close to the correct azimuth angle (dashed line) in about several seconds.
(3) When the vehicle 200 was stopped, the azimuth value started to drift due to the bias (measurement error) of the angular velocity sensor 222.

図7(1)(2)に示すように、3自由度2次元航法は、方位角の初期値が不適切な値であってもすぐに適切な値を求めることができる。
これは、3自由度2次元航法において、加速度センサ221が計測した加速度ではなく、車速センサ230が計測した速度を用いているためである。
これにより、二次元座標を求める積分計算が1回で済み、速度の座標変換に用いる方位角(の初期値)の誤差の影響が比較的少なく、カルマンフィルタを用いて二次元座標と方位角とを適切に補正できる。
As shown in FIGS. 7A and 7B, the three-degree-of-freedom two-dimensional navigation can immediately obtain an appropriate value even if the initial value of the azimuth is an inappropriate value.
This is because in three-degree-of-freedom two-dimensional navigation, not the acceleration measured by the acceleration sensor 221 but the speed measured by the vehicle speed sensor 230 is used.
As a result, the integration calculation for obtaining the two-dimensional coordinates is only required once, and the influence of the error of the azimuth angle (initial value) used for the coordinate conversion of the speed is relatively small. Can be corrected appropriately.

図5に戻り、S110の説明を続ける。   Returning to FIG. 5, the description of S110 will be continued.

カルマンフィルタ部130は、上記の3自由度2次元航法(図6参照)を所定の初期化時間(例えば、数秒)または所定の航法回数だけ実行し、実行を終了したときに算出した方位角を「方位角初期値」としてストラップダウン演算部140に出力する。
S110の後、S120に進む。
The Kalman filter unit 130 executes the above-described three-degree-of-freedom two-dimensional navigation (see FIG. 6) for a predetermined initialization time (for example, several seconds) or a predetermined number of times of navigation, and the azimuth angle calculated when the execution is completed is It is output to the strapdown calculator 140 as the “azimuth initial value”.
It progresses to S120 after S110.

S120において、走行時標定部121は、方位角の初期値または前回標定した方位角を用いて6自由度3次元航法(第二航法)により位置と姿勢と方位とを標定する。
「6自由度」は「緯度」「経度」「高度」「仰角」「回転角」「方位角」を意味し、「3次元」は「緯度、経度、高度」または「仰角、回転角、方位角」を意味する。
In S120, the travel time location unit 121 locates the position, orientation, and orientation by the 6-degree-of-freedom three-dimensional navigation (second navigation) using the initial value of the azimuth angle or the previously determined azimuth angle.
“6 degrees of freedom” means “latitude” “longitude” “altitude” “elevation angle” “rotation angle” “azimuth angle”, and “three-dimensional” means “latitude, longitude, altitude” or “elevation angle, rotation angle, direction” Means "horn".

図8は、実施の形態1における6自由度3次元航法のアルゴリズムを示す図である。
実施の形態1における6自由度3次元航法のアルゴリズムについて、図8に基づいて説明する。
FIG. 8 is a diagram showing an algorithm for six-degree-of-freedom three-dimensional navigation in the first embodiment.
The 6-degree-of-freedom three-dimensional navigation algorithm in Embodiment 1 will be described with reference to FIG.

停車判定部123、入力補正部124、出力補正部125は自己位置姿勢標定部120が備える機能構成である。   The stop determination unit 123, the input correction unit 124, and the output correction unit 125 are functional configurations provided in the self-position / posture determination unit 120.

入力補正部124は、加速度データ193に含まれる三次元加速度にカルマンフィルタ部130により算出される補正値を加算して三次元加速度を補正する。
さらに、入力補正部124は、角速度データ194に含まれる三次元角速度にカルマンフィルタ部130により算出される補正値を加算して三次元角速度を補正する。
The input correction unit 124 adds the correction value calculated by the Kalman filter unit 130 to the three-dimensional acceleration included in the acceleration data 193 to correct the three-dimensional acceleration.
Further, the input correction unit 124 adds the correction value calculated by the Kalman filter unit 130 to the three-dimensional angular velocity included in the angular velocity data 194 to correct the three-dimensional angular velocity.

ストラップダウン演算部140は、入力補正部124により補正された三次元加速度と入力補正部124により補正された三次元角速度とを用いて位置と姿勢と方位とを算出する(図4参照)。
ここで、ストラップダウン演算部140の姿勢方位更新部142は、初期値算出部110により算出された方位角を方位角の初期値として用いる。
図7で説明したように、初期値算出部110により算出される方位角は正確な方位角に近い値であり、ストラップダウン演算に用いる方位角の初期値として適切な値である。例えば、初期値算出部110の代わりに方位磁石を用いた場合、ストラップダウン演算で用いる初期値として十分な精度の方位角を得ることができない。
ストラップダウン演算部140は、適切な方位角の初期値を用いることにより精度が高い位置と姿勢と方位とを算出することができる。
The strapdown calculation unit 140 calculates the position, orientation, and orientation using the three-dimensional acceleration corrected by the input correction unit 124 and the three-dimensional angular velocity corrected by the input correction unit 124 (see FIG. 4).
Here, the posture / azimuth update unit 142 of the strapdown calculation unit 140 uses the azimuth angle calculated by the initial value calculation unit 110 as the initial value of the azimuth angle.
As described with reference to FIG. 7, the azimuth angle calculated by the initial value calculation unit 110 is a value close to an accurate azimuth angle, and is an appropriate value as the initial value of the azimuth angle used for the strapdown calculation. For example, when an azimuth magnet is used instead of the initial value calculation unit 110, an azimuth angle with sufficient accuracy cannot be obtained as an initial value used in strapdown calculation.
The strapdown computing unit 140 can calculate a position, orientation, and orientation with high accuracy by using an initial value of an appropriate azimuth angle.

出力補正部125は、ストラップダウン演算部140により算出された位置にカルマンフィルタ部130により算出される補正値を加算して位置を補正し、補正した位置を標定結果の位置として出力する。
また、出力補正部125は、ストラップダウン演算部140により算出された姿勢にカルマンフィルタ部130により算出される補正値を加算して姿勢を補正し、補正した姿勢を標定結果の姿勢角として出力する。
さらに、出力補正部125は、ストラップダウン演算部140により算出された方位にカルマンフィルタ部130により算出される補正値を加算して方位を補正し、補正した方位を標定結果の方位として出力する。
The output correction unit 125 corrects the position by adding the correction value calculated by the Kalman filter unit 130 to the position calculated by the strapdown calculation unit 140, and outputs the corrected position as the position of the orientation result.
The output correction unit 125 corrects the posture by adding the correction value calculated by the Kalman filter unit 130 to the posture calculated by the strapdown calculation unit 140, and outputs the corrected posture as the posture angle of the orientation result.
Further, the output correction unit 125 corrects the azimuth by adding the correction value calculated by the Kalman filter unit 130 to the azimuth calculated by the strapdown calculation unit 140, and outputs the corrected azimuth as the azimuth of the orientation result.

停車判定部123は、車両200が停車しているか否かを判定する。   The stop determination unit 123 determines whether or not the vehicle 200 is stopped.

カルマンフィルタ部130は、車両200が停車していない場合、GPSデータ191に含まれる位置とストラップダウン演算部140により算出された位置との差を「位置残差」として算出する。
また、カルマンフィルタ部130は、速度データ192に含まれる速度とストラップダウン演算部140により算出された速度との差を「速度残差」として算出する。
そして、カルマンフィルタ部130は、位置残差と速度残差とを入力値として第2のカルマンフィルタ132を実行し、三次元加速度の補正値と三次元角速度の補正値と位置の補正値と姿勢の補正値と方位の補正値とを算出する。
When the vehicle 200 is not stopped, the Kalman filter unit 130 calculates a difference between the position included in the GPS data 191 and the position calculated by the strapdown calculation unit 140 as a “position residual”.
Further, the Kalman filter unit 130 calculates a difference between the speed included in the speed data 192 and the speed calculated by the strapdown calculation unit 140 as a “speed residual”.
The Kalman filter unit 130 executes the second Kalman filter 132 using the position residual and the velocity residual as input values, and corrects the three-dimensional acceleration correction value, the three-dimensional angular velocity correction value, the position correction value, and the posture correction. The value and the correction value of the direction are calculated.

カルマンフィルタ部130は、車両200が停車している場合、ストラップダウン演算部140により算出された方位に基づいて方位の変化量を「方位残差」として算出する。
そして、カルマンフィルタ部130は、方位残差を入力値として第3のカルマンフィルタ133を実行し、三次元加速度の補正値と三次元角速度の補正値と位置の補正値と姿勢の補正値と方位の補正値とを算出する。
When the vehicle 200 is stopped, the Kalman filter unit 130 calculates an azimuth change amount as an “azimuth residual” based on the azimuth calculated by the strapdown calculation unit 140.
The Kalman filter unit 130 executes the third Kalman filter 133 using the azimuth residual as an input value, and corrects the three-dimensional acceleration correction value, the three-dimensional angular velocity correction value, the position correction value, the posture correction value, and the heading correction. Value.

図5に戻り、自己位置姿勢標定方法の説明を続ける。   Returning to FIG. 5, the description of the self-position / posture locating method will be continued.

S120の後、S130に進む。   It progresses to S130 after S120.

S130において、停車判定部123は、速度データ192に含まれる速度を所定の速度閾値と比較し、速度が速度閾値未満である場合、車両200が停車していると判定し、速度が速度閾値以上である場合、車両200が停車していないと判定する。
但し、停車判定部123は、車両200が停車しているか否かをその他のデータに基づいて判定してもよい。例えば、停車判定部123は、GPSデータ191に含まれる速度に基づいて判定を行ってもよい。また、停車判定部123は、GPSデータ191に含まれる前回の位置と今回の位置とに基づいて位置の変化量が所定の位置閾値未満である場合に車両200が停車していると判定してもよい。また、停車判定部123は、加速度データ193(または角速度データ194)に基づいて加速度(または角速度)の変化量が所定の加速度閾値(または角速度閾値)未満である場合に車両200が停車していると判定してもよい。
車両200が停車している場合(YES)、S140に進む。
車両200が停車していない場合(NO)、S120に戻る。
In S130, the stop determination unit 123 compares the speed included in the speed data 192 with a predetermined speed threshold, and determines that the vehicle 200 is stopped if the speed is less than the speed threshold, and the speed is equal to or greater than the speed threshold. If it is, it is determined that the vehicle 200 is not stopped.
However, the stop determination unit 123 may determine whether or not the vehicle 200 is stopped based on other data. For example, the stop determination unit 123 may perform the determination based on the speed included in the GPS data 191. Further, the stop determination unit 123 determines that the vehicle 200 is stopped based on the previous position and the current position included in the GPS data 191 when the change amount of the position is less than a predetermined position threshold. Also good. Further, the stop determination unit 123 stops the vehicle 200 when the change amount of the acceleration (or angular velocity) is less than a predetermined acceleration threshold (or angular velocity threshold) based on the acceleration data 193 (or angular velocity data 194). May be determined.
When the vehicle 200 is stopped (YES), the process proceeds to S140.
When the vehicle 200 is not stopped (NO), the process returns to S120.

S140において、停車時標定部122は、前回標定された方位角を用いて6自由度3次元航法(図8参照)により位置と姿勢と方位とを標定する。
但し、停車時標定部122は、6自由度3次元航法に「拡張ZUPTアルゴリズム」を適用する。
In S140, the stopping time locating unit 122 locates the position, orientation, and azimuth by the 6-degree-of-freedom three-dimensional navigation (see FIG. 8) using the previously determined azimuth angle.
However, the stopping time locating unit 122 applies the “extended ZUPT algorithm” to the 6-degree-of-freedom three-dimensional navigation.

拡張ZUPTアルゴリズムは、ZUPT(Zero velocity UP daTe)という従来技術を拡張した新規技術である。ZUPTとは、停止中の位置の計測値を固定することによりセンサ誤差による位置(計測値)の変化を防ぐ技術である。
拡張ZUPTアルゴリズムは、停止中の方位角の計測値を固定することによりセンサ誤差による方位角(計測値)の変化を防ぐ。
但し、拡張ZUPTアルゴリズムは、停止中の姿勢角の計測値を固定しない。車両200が停車中であっても人の乗り降りや車内での人の移動によって、車両200の姿勢角は変化するからである。もし、姿勢角の計測値を固定してしまうと姿勢角が実際に変化しているにも関わらず、カルマンフィルタが姿勢角の変化量をセンサ誤差として判断してしまい、間違った補正値を出力してしまう。また、姿勢角はレベリングによって正確な値を求めることができる。
The extended ZUPT algorithm is a new technology that extends the conventional technology called ZUPT (Zero Velocity UP daTe). ZUPT is a technique for preventing a change in position (measured value) due to a sensor error by fixing a measured value at a stopped position.
The extended ZUPT algorithm prevents a change in azimuth angle (measurement value) due to a sensor error by fixing the measurement value of the azimuth angle during a stop.
However, the extended ZUPT algorithm does not fix the measured value of the posture angle while stopped. This is because, even when the vehicle 200 is stopped, the posture angle of the vehicle 200 changes depending on the person getting on and off and the movement of the person inside the vehicle. If the measured value of the posture angle is fixed, the Kalman filter will determine the amount of change in the posture angle as a sensor error even though the posture angle has actually changed, and output an incorrect correction value. End up. Further, the posture angle can be accurately determined by leveling.

図9は、実施の形態1における車両200の停車時の方位角のドリフトを示すグラフである。
車両200の停車時(但し、人の乗り降り及び車内での人の移動は無し)の回転角、仰角、方位角の変化を図9に示す。
図9に示すように、方位角の値は角速度センサ222のバイアス誤差により、時間の経過に伴ってドリフトする。
一方、回転角、仰角の値(姿勢角の値)は、重力加速度と加速度センサ221により計測された加速度との比較によって安定化することができる(レベリング)。
FIG. 9 is a graph showing azimuth drift when the vehicle 200 stops in the first embodiment.
FIG. 9 shows changes in the rotation angle, the elevation angle, and the azimuth angle when the vehicle 200 is stopped (however, no people get on and off and no person moves in the vehicle).
As shown in FIG. 9, the azimuth value drifts with time due to the bias error of the angular velocity sensor 222.
On the other hand, the rotation angle and elevation angle values (posture angle values) can be stabilized by comparing the gravitational acceleration with the acceleration measured by the acceleration sensor 221 (leveling).

図10は、実施の形態1における拡張ZUPTアルゴリズムの概要図である。
実施の形態1における拡張ZUPTアルゴリズムの概要について、図10に基づいて説明する。
FIG. 10 is a schematic diagram of the extended ZUPT algorithm in the first embodiment.
An overview of the extended ZUPT algorithm in the first embodiment will be described with reference to FIG.

拡張ZUPTアルゴリズムでは、車両200の左右に「仮想アンカポイント」を設定し、車両200の長さ方向の2箇所に「仮想センシングポイント」を設定し、仮想アンカポイントと仮想センシングポイントとをワイヤーで結んだ状態を仮定する。
この場合、車両200の姿勢角(x軸またはy軸回りの回転)を変化させることは可能であるが、車両200の方位角(z軸回りの回転)を変化させることはできない。つまり、方位角について車両200が拘束される。
In the extended ZUPT algorithm, “virtual anchor points” are set on the left and right sides of the vehicle 200, “virtual sensing points” are set at two locations in the longitudinal direction of the vehicle 200, and the virtual anchor points and the virtual sensing points are connected by wires. Assume that
In this case, it is possible to change the attitude angle (rotation around the x-axis or y-axis) of the vehicle 200, but it is not possible to change the azimuth angle (rotation around the z-axis) of the vehicle 200. That is, the vehicle 200 is restrained about the azimuth angle.

この状態において、仮想アンカポイントと仮想センシングポイントとの距離の二重差を観測値とする。
観測値を表す観測式(1)−(3)を以下に示す。式(3)に示すように観測値は「0」である。
In this state, the double difference in the distance between the virtual anchor point and the virtual sensing point is taken as the observation value.
Observation formulas (1) to (3) representing observation values are shown below. As shown in Expression (3), the observed value is “0”.

Figure 0005602070
Figure 0005602070

上記式(2)の行列Aを表す式(6)および行列Aの算出式(4)(5)を以下に示す。   Expression (6) representing the matrix A of the above expression (2) and calculation expressions (4) and (5) of the matrix A are shown below.

Figure 0005602070
Figure 0005602070

上記式(4)の右辺の行列は、方位角を車体座標系(xyz)から局所水平座標系(NED)に変換する方向余弦行列である。
また、上記式(4)の右辺の1つ目のベクトルは車両200の右方向の単位行列であり、2つ目のベクトルは車両200の左方向の単位行列である。
The matrix on the right side of the above equation (4) is a direction cosine matrix for converting the azimuth angle from the vehicle body coordinate system (xyz) to the local horizontal coordinate system (NED).
Further, the first vector on the right side of Equation (4) is a unit matrix in the right direction of the vehicle 200, and the second vector is a unit matrix in the left direction of the vehicle 200.

上記式(2)に含まれる行列Aと行列Cとベクトルbとのうち、行列Cはセンサ誤差の影響による方位角の計測値の変化に応じて変化する。
そこで、予測値を以下の観測式(7)で表す。
Of the matrix A, the matrix C, and the vector b included in the above equation (2), the matrix C changes according to the change in the measured value of the azimuth angle due to the influence of the sensor error.
Therefore, the predicted value is expressed by the following observation formula (7).

Figure 0005602070
Figure 0005602070

式(2)に表した観測値と式(7)に表した予測値との差を方位角の観測誤差(方位残差)として以下の式(8)で表す。   The difference between the observed value represented by Expression (2) and the predicted value represented by Expression (7) is represented by the following Expression (8) as an azimuth observation error (azimuth residual).

Figure 0005602070
Figure 0005602070

図5に戻り、S140の説明を続ける。   Returning to FIG. 5, the description of S140 will be continued.

6自由度3次元航法(図8)において、拡張ZUPTアルゴリズムを適用した第3のカルマンフィルタ133が実行される。
第3のカルマンフィルタ133は、上記式(8)を観測方程式として用いてカルマンフィルタ処理により各種の補正値を算出する。
S140の後、S150に進む。
In the six-degree-of-freedom three-dimensional navigation (FIG. 8), the third Kalman filter 133 to which the extended ZUPT algorithm is applied is executed.
The third Kalman filter 133 calculates various correction values by Kalman filter processing using the above equation (8) as an observation equation.
After S140, the process proceeds to S150.

S150において、停車判定部123は、車両200が停車しているか否かをS130と同様に判定する。
車両200が停車している場合(YES)、S140に戻る。
車両200が停車していない場合(NO)、S120に戻る。
In S150, the stop determination unit 123 determines whether or not the vehicle 200 is stopped as in S130.
If the vehicle 200 is stopped (YES), the process returns to S140.
When the vehicle 200 is not stopped (NO), the process returns to S120.

図11は、実施の形態1における自己位置姿勢標定装置100のハードウェア資源の一例を示す図である。
図11において、自己位置姿勢標定装置100は、CPU901(Central Processing Unit)を備えている。CPU901は、バス902を介してROM903、RAM904、通信ボード905、磁気ディスク装置920と接続され、これらのハードウェアデバイスを制御する。
FIG. 11 is a diagram illustrating an example of hardware resources of the self-position / posture locating apparatus 100 according to the first embodiment.
In FIG. 11, the self-position / posture locating apparatus 100 includes a CPU 901 (Central Processing Unit). The CPU 901 is connected to the ROM 903, the RAM 904, the communication board 905, and the magnetic disk device 920 via the bus 902, and controls these hardware devices.

通信ボード905は、有線または無線で、LAN(Local Area Network)、インターネット、電話回線などの通信網に接続している。   The communication board 905 is wired or wirelessly connected to a communication network such as a LAN (Local Area Network), the Internet, or a telephone line.

磁気ディスク装置920には、OS921(オペレーティングシステム)、プログラム群922、ファイル群923が記憶されている。   The magnetic disk device 920 stores an OS 921 (operating system), a program group 922, and a file group 923.

プログラム群922には、実施の形態において「〜部」として説明する機能を実行するプログラムが含まれる。プログラムは、CPU901により読み出され実行される。すなわち、プログラムは、「〜部」としてコンピュータを機能させるものであり、また「〜部」の手順や方法をコンピュータに実行させるものである。   The program group 922 includes programs that execute the functions described as “units” in the embodiments. The program is read and executed by the CPU 901. That is, the program causes the computer to function as “to part”, and causes the computer to execute the procedures and methods of “to part”.

ファイル群923には、実施の形態において説明する「〜部」で使用される各種データ(入力、出力、判定結果、計算結果、処理結果など)が含まれる。   The file group 923 includes various data (input, output, determination result, calculation result, processing result, etc.) used in “˜part” described in the embodiment.

実施の形態において構成図およびフローチャートに含まれている矢印は主としてデータや信号の入出力を示す。   In the embodiment, arrows included in the configuration diagrams and flowcharts mainly indicate input and output of data and signals.

実施の形態において「〜部」として説明するものは「〜回路」、「〜装置」、「〜機器」であってもよく、また「〜ステップ」、「〜手順」、「〜処理」であってもよい。すなわち、「〜部」として説明するものは、ファームウェア、ソフトウェア、ハードウェアまたはこれらの組み合わせのいずれで実装されても構わない。   In the embodiment, what is described as “to part” may be “to circuit”, “to apparatus”, and “to device”, and “to step”, “to procedure”, and “to processing”. May be. That is, what is described as “to part” may be implemented by any of firmware, software, hardware, or a combination thereof.

実施の形態1において、例えば、以下のような自己位置姿勢標定装置100について説明した。   In the first embodiment, for example, the following self-position / posture locating apparatus 100 has been described.

自己位置姿勢標定装置100は、GPSとIMUと車速パルスとを用いて(1)方位を求めるための初期演算と、(2)停車中に実施する方位に関するカルマンフィルタの観測更新アルゴリズム(拡張ZUPT)とを実行する。
(1)(2)により、低価格で中精度なIMUと設置場所に制約がない1台のGPSとを用いて6自由度3次元航法を行うことが可能となる。
The self-position / posture locator 100 uses GPS, IMU, and vehicle speed pulse to (1) an initial calculation for obtaining an azimuth; Execute.
(1) By (2), it is possible to perform six-degree-of-freedom three-dimensional navigation using a low-priced, medium-precision IMU and a single GPS with no restrictions on the installation location.

100 自己位置姿勢標定装置、110 初期値算出部、111 入力補正部、112 方位更新部、113 座標変換部、114 緯度経度更新部、115 出力補正部、120 自己位置姿勢標定部、121 走行時標定部、122 停車時標定部、123 停車判定部、124 入力補正部、125 出力補正部、130 カルマンフィルタ部、131 第1のカルマンフィルタ、132 第2のカルマンフィルタ、133 第3のカルマンフィルタ、140 ストラップダウン演算部、141 角速度座標変換部、142 姿勢方位更新部、143 加速度座標変換部、144 速度位置更新部、145 重力補正計算部、146 コリオリ補正計算部、147 航法系回転補正計算部、148A 第一補正部、148B 第二補正部、148C 第三補正部、190 装置記憶部、191 GPSデータ、192 速度データ、193 加速度データ、194 角速度データ、200 車両、201 天板、210 GPS受信機、211 GPSアンテナ、220 慣性センサ、221 加速度センサ、222 角速度センサ、230 車速センサ、901 CPU、902 バス、903 ROM、904 RAM、905 通信ボード、920 磁気ディスク装置、921 OS、922 プログラム群、923 ファイル群。   DESCRIPTION OF SYMBOLS 100 Self-position posture determination apparatus, 110 Initial value calculation part, 111 Input correction part, 112 Direction update part, 113 Coordinate conversion part, 114 Latitude / longitude update part, 115 Output correction part, 120 Self-position posture determination part, 121 Travel time standardization , 122 Stop time locating unit, 123 Stop determining unit, 124 Input correction unit, 125 Output correction unit, 130 Kalman filter unit, 131 First Kalman filter, 132 Second Kalman filter, 133 Third Kalman filter, 140 Strapdown calculation unit , 141 angular velocity coordinate conversion unit, 142 posture orientation update unit, 143 acceleration coordinate conversion unit, 144 speed position update unit, 145 gravity correction calculation unit, 146 Coriolis correction calculation unit, 147 navigation system rotation correction calculation unit, 148A first correction unit 148B Second correction unit, 148C Third Normal part, 190 device storage part, 191 GPS data, 192 speed data, 193 acceleration data, 194 angular speed data, 200 vehicle, 201 top plate, 210 GPS receiver, 211 GPS antenna, 220 inertial sensor, 221 acceleration sensor, 222 angular speed Sensor, 230 Vehicle speed sensor, 901 CPU, 902 bus, 903 ROM, 904 RAM, 905 communication board, 920 magnetic disk device, 921 OS, 922 program group, 923 file group.

Claims (7)

GPS(Global Positioning System)を利用して測位を行うGPS受信機と、速度を計測する速度センサと、加速度を計測する加速度センサと、方位角の角速度を計測する角速度センサとが取り付けられた移動体の位置を標定する位置標定装置において、
前記角速度センサにより計測された方位角の角速度に基づいて前記移動体の方位角を第一航法方位角として算出し、算出した第一航法方位角に基づいて、前記速度センサにより計測された速さを二次元の速度に変換し、変換後の二次元の速度を積分して二次元座標における変化量を算出し、算出した前記変化量を二次元の位置初期値に加えることにより前記移動体の二次元の座標値を第一航法座標値として算出し、算出した第一航法座標値と前記GPS受信機により測位された座標値に含まれる二次元の座標値とに基づいて第1のカルマンフィルタを用いて前記第一航法方位角を補正する補正値を算出し、算出した補正値を用いて前記第一航法方位角を補正し、補正後の第一航法方位角を方位角初期値として出力する初期値算出部と、
前記角速度センサにより計測された方位角の角速度に基づいて方位角の変化量を算出し、算出した方位角の変化量と前記初期値算出部から出力された方位角初期値とに基づいて前記移動体の方位角を第二航法方位角として算出し、算出した第二航法方位角に基づいて、前記加速度センサにより計測された三次元の加速度を航法座標系の三次元の加速度に変換し、変換後の三次元の加速度を積分して三次元座標における変化量を算出し、算出した前記変化量を三次元の位置初期値に加えることにより前記移動体の三次元の座標値を第二航法座標値として算出し、算出した第二航法座標値と前記GPS受信機により測位された三次元の座標値とに基づいて第2のカルマンフィルタを用いて前記第二航法座標値を補正する補正値を算出し、算出した補正値を用いて前記第二航法座標値を補正し、前記移動体の位置を標定した位置標定値として補正後の第二航法座標値を出力する標定値算出部とを備えたことを特徴とする位置標定装置。
A mobile body equipped with a GPS receiver that performs positioning using GPS (Global Positioning System), a speed sensor that measures speed, an acceleration sensor that measures acceleration, and an angular speed sensor that measures angular velocity of an azimuth angle In a position locator that locates the position of
Wherein calculating the azimuth angle of the mobile object as the first navigation orientation angle based on the angular velocity of azimuth measured by said angular velocity sensor, based on the first navigation orientation angles calculated speed measured by said speed sensor Is converted into a two-dimensional velocity, and the converted two-dimensional velocity is integrated to calculate a change amount in two-dimensional coordinates, and the calculated change amount is added to a two-dimensional position initial value to A two-dimensional coordinate value is calculated as a first navigation coordinate value, and the first Kalman filter is calculated based on the calculated first navigation coordinate value and the two-dimensional coordinate value included in the coordinate value measured by the GPS receiver. A correction value for correcting the first navigation azimuth is calculated, the first navigation azimuth is corrected using the calculated correction value, and the corrected first navigation azimuth is output as an initial azimuth value. An initial value calculation unit;
The amount of change in azimuth is calculated based on the angular velocity of the azimuth measured by the angular velocity sensor, and the movement is performed based on the calculated amount of change in azimuth and the initial value of azimuth output from the initial value calculator. The azimuth angle of the body is calculated as the second navigation azimuth angle, and based on the calculated second navigation azimuth angle , the three-dimensional acceleration measured by the acceleration sensor is converted into the three-dimensional acceleration of the navigation coordinate system, and converted. after integrating the three-dimensional acceleration to calculate the amount of change in the three-dimensional coordinates, calculated second navigation coordinate the coordinate values of the three-dimensional of the movable body by adding the amount of change in position initial value of the three-dimensional And a correction value for correcting the second navigation coordinate value using the second Kalman filter based on the calculated second navigation coordinate value and the three-dimensional coordinate value measured by the GPS receiver. And calculated A ground value calculation unit that corrects the second navigation coordinate value using a positive value and outputs the corrected second navigation coordinate value as a position position value obtained by locating the position of the moving body. Positioning device to do.
前記標定値算出部は、
移動していた前記移動体が移動を停止した場合、前回算出された第二航法方位角と今回算出された第二航法方位角とに基づいて第3のカルマンフィルタを用いて今回算出された前記第二航法方位角を補正する方位角補正値を算出し、算出した方位角補正値を用いて今回算出された前記第二航法方位角を補正して補正後の第二航法方位角を算出し、
停止した前記移動体が移動を再開した場合、前記補正後の第二航法方位角と新たな方位角の変化量とに基づいて新たな第二航法方位角を算出し、算出した新たな第二航法方位角に基づいて、前記加速度センサにより計測された新たな三次元の加速度を航法座標系の三次元の加速度に変換し、変換後の三次元の加速度を積分して三次元座標における変化量を算出し、算出した変化量を前回算出した第二航法座標値に加えることにより、新たな第二航法座標値を算出し、算出した新たな第二航法座標値と前記GPS受信機により測位された新たな三次元の座標値とに基づいて新たな第二航法座標値を補正する位置補正値を前記第2のカルマンフィルタを用いて算出し、算出した位置補正値を用いて前記新たな第二航法座標値を補正し、補正後の新たな第二航法座標値を新たな位置標定値として出力する
ことを特徴とする請求項1記載の位置標定装置。
The orientation value calculation unit
When the moving object that has been moving stops moving, the first calculated by using a third Kalman filter based on the previously calculated second navigation azimuth and the second navigation azimuth calculated this time . Calculate the azimuth correction value for correcting the two navigation azimuths, calculate the second navigation azimuth after correction by correcting the second navigation azimuth calculated this time using the calculated azimuth correction value,
When the stopped moving body resumes moving, a new second navigation azimuth is calculated based on the corrected second navigation azimuth and the amount of change in the new azimuth, and the calculated new second azimuth is calculated. Based on the navigation azimuth, the new three-dimensional acceleration measured by the acceleration sensor is converted into the three-dimensional acceleration of the navigation coordinate system, and the converted three-dimensional acceleration is integrated to change the three-dimensional coordinate. And the calculated change amount is added to the previously calculated second navigation coordinate value to calculate a new second navigation coordinate value, which is measured by the calculated second navigation coordinate value and the GPS receiver. A position correction value for correcting a new second navigation coordinate value based on the new three-dimensional coordinate value is calculated using the second Kalman filter, and the new second value is calculated using the calculated position correction value. The navigation coordinate value is corrected, and the corrected new value Position locating system according to claim 1, wherein the outputting a second navigation coordinate value as a new position location value.
前記標定値算出部は、
前記速度センサにより計測された速度が所定の速度閾値未満である場合、前記移動体が移動を停止したと判定し、
前記移動体が移動を停止したと判定した後、前記速度センサにより計測された速度が前記速度閾値以上である場合、前記移動体が移動を再開したと判定する
ことを特徴とする請求項2記載の位置標定装置。
The orientation value calculation unit
When the speed measured by the speed sensor is less than a predetermined speed threshold, it is determined that the moving body has stopped moving,
3. The method according to claim 2, wherein after determining that the moving body has stopped moving, if the speed measured by the speed sensor is equal to or greater than the speed threshold value, it is determined that the moving body has resumed moving. Position locator.
前記第3のカルマンフィルタは、前回算出された第二航法方位角と今回算出された第二航法方位角とを一致させる補正値を前記方位角補正値として算出する
ことを特徴とする請求項2または請求項3記載の位置標定装置。
The third Kalman filter calculates, as the azimuth angle correction value, a correction value that matches the previously calculated second navigation azimuth angle with the second navigation azimuth angle calculated this time. The position locating device according to claim 3.
前記第3のカルマンフィルタは、前回算出された第二航法方位角と今回算出された第二航法方位角との差分を表す観測方程式を生成し、生成した観測方程式を用いて前記方位角補正値を算出する
ことを特徴とする請求項2から請求項4いずれかに記載の位置標定装置。
The third Kalman filter generates an observation equation representing a difference between the previously calculated second navigation azimuth and the currently calculated second navigation azimuth, and uses the generated observation equation to calculate the azimuth correction value. The position locating device according to any one of claims 2 to 4, wherein the position locating device is calculated.
GPS(Global Positioning System)を利用して測位を行うGPS受信機と、速度を計測する速度センサと、加速度を計測する加速度センサと、方位角の角速度を計測する角速度センサとが取り付けられた移動体の位置を標定する位置標定装置の位置標定方法において、
初期値算出部は、前記角速度センサにより計測された方位角の角速度に基づいて前記移動体の方位角を第一航法方位角として算出し、算出した第一航法方位角に基づいて、前記速度センサにより計測された速さを二次元の速度に変換し、変換後の二次元の速度を積分して二次元座標における変化量を算出し、算出した前記変化量を二次元の位置初期値に加えることにより前記移動体の二次元の座標値を第一航法座標値として算出し、算出した第一航法座標値と前記GPS受信機により測位された座標値に含まれる二次元の座標値とに基づいて第1のカルマンフィルタを用いて前記第一航法方位角を補正する補正値を算出し、算出した補正値を用いて前記第一航法方位角を補正し、補正後の第一航法方位角を方位角初期値として出力し、標定値算出部は、前記角速度センサにより計測された方位角の角速度に基づいて方位角の変化量を算出し、算出した方位角の変化量と前記初期値算出部から出力された方位角初期値とに基づいて前記移動体の方位角を第二航法方位角として算出し、算出した第二航法方位角に基づいて、前記加速度センサにより計測された三次元の加速度を航法座標系の三次元の加速度に変換し、変換後の三次元の加速度を積分して三次元座標における変化量を算出し、算出した前記変化量を三次元の位置初期値に加えることにより前記移動体の三次元の座標値を第二航法座標値として算出し、算出した第二航法座標値と前記GPS受信機により測位された三次元の座標値とに基づいて第2のカルマンフィルタを用いて前記第二航法座標値を補正する補正値を算出し、算出した補正値を用いて前記第二航法座標値を補正し、前記移動体の位置を標定した位置標定値として補正後の第二航法座標値を出力する
ことを特徴とする位置標定装置の位置標定方法。
A mobile body equipped with a GPS receiver that performs positioning using GPS (Global Positioning System), a speed sensor that measures speed, an acceleration sensor that measures acceleration, and an angular speed sensor that measures angular velocity of an azimuth angle In the position locating method of the position locating device for locating the position of
The initial value calculation unit calculates the azimuth angle of the moving body as a first navigation azimuth angle based on the angular velocity of the azimuth angle measured by the angular velocity sensor, and based on the calculated first navigation azimuth angle , the speed sensor Convert the speed measured in step 2 into a two-dimensional speed, integrate the converted two-dimensional speed to calculate the amount of change in two-dimensional coordinates, and add the calculated amount of change to the two-dimensional initial position value Thus, the two-dimensional coordinate value of the moving body is calculated as the first navigation coordinate value, and based on the calculated first navigation coordinate value and the two-dimensional coordinate value included in the coordinate value measured by the GPS receiver. Then, a correction value for correcting the first navigation azimuth is calculated using the first Kalman filter, the first navigation azimuth is corrected using the calculated correction value, and the corrected first navigation azimuth is determined as the azimuth. Output as angle initial value, orientation value The output unit calculates an azimuth angle change amount based on the angular velocity of the azimuth angle measured by the angular velocity sensor, and calculates the calculated azimuth angle change amount and the azimuth angle initial value output from the initial value calculation unit. Based on the calculated second navigation azimuth , the three-dimensional acceleration measured by the acceleration sensor is converted into the three-dimensional acceleration of the navigation coordinate system based on the calculated second navigation azimuth. Converting, integrating the converted three-dimensional acceleration to calculate the amount of change in three-dimensional coordinates, and adding the calculated amount of change to the initial three-dimensional position value to obtain the three-dimensional coordinate value of the moving object. Calculate as the second navigation coordinate value, and correct the second navigation coordinate value using the second Kalman filter based on the calculated second navigation coordinate value and the three-dimensional coordinate value measured by the GPS receiver. Calculate the correction value The position of the position locating device, wherein the second navigation coordinate value is corrected using the calculated correction value, and the corrected second navigation coordinate value is output as a position locating value obtained by locating the position of the moving body. Orientation method.
GPS(Global Positioning System)を利用して測位を行うGPS受信機と、速度を計測する速度センサと、加速度を計測する加速度センサと、方位角の角速度を計測する角速度センサとが取り付けられた移動体の位置を標定する位置標定プログラムにおいて、
前記角速度センサにより計測された方位角の角速度に基づいて前記移動体の方位角を第一航法方位角として算出し、算出した第一航法方位角に基づいて、前記速度センサにより計測された速さを二次元の速度に変換し、変換後の二次元の速度を積分して二次元座標における変化量を算出し、算出した前記変化量を二次元の位置初期値に加えることにより前記移動体の二次元の座標値を第一航法座標値として算出し、算出した第一航法座標値と前記GPS受信機により測位された座標値に含まれる二次元の座表示とに基づいて第1のカルマンフィルタを用いて前記第一航法方位角を補正する補正値を算出し、算出した補正値を用いて前記第一航法方位角を補正し、補正後の第一航法方位角を方位角初期値として出力する初期値算出部と、
前記角速度センサにより計測された方位角の角速度に基づいて方位角の変化量を算出し、算出した方位角の変化量と前記初期値算出部から出力された方位角初期値とに基づいて前記移動体の方位角を第二航法方位角として算出し、算出した第二航法方位角に基づいて、前記加速度センサにより計測された三次元の加速度を航法座標系の三次元の加速度に変換し、変換後の三次元の加速度を積分して三次元座標における変化量を算出し、算出した前記変化量を三次元の位置初期値に加えることにより前記移動体の三次元の座標値を第二航法座標値として算出し、算出した第二航法座標値と前記GPS受信機により測位された三次元の座標値とに基づいて第2のカルマンフィルタを用いて前記第二航法座標値を補正する補正値を算出し、算出した補正値を用いて前記第二航法座標値を補正し、前記移動体の位置を標定した位置標定値として補正後の第二航法座標値を出力する標定値算出部として
コンピュータを機能させることを特徴とする位置標定プログラム。
A mobile body equipped with a GPS receiver that performs positioning using GPS (Global Positioning System), a speed sensor that measures speed, an acceleration sensor that measures acceleration, and an angular speed sensor that measures angular velocity of an azimuth angle In the location program that locates the position of
Wherein calculating the azimuth angle of the mobile object as the first navigation orientation angle based on the angular velocity of azimuth measured by said angular velocity sensor, based on the first navigation orientation angles calculated speed measured by said speed sensor Is converted into a two-dimensional velocity, and the converted two-dimensional velocity is integrated to calculate a change amount in two-dimensional coordinates, and the calculated change amount is added to a two-dimensional position initial value to A two-dimensional coordinate value is calculated as a first navigation coordinate value, and a first Kalman filter is calculated based on the calculated first navigation coordinate value and a two-dimensional seat display included in the coordinate value measured by the GPS receiver. A correction value for correcting the first navigation azimuth is calculated, the first navigation azimuth is corrected using the calculated correction value, and the corrected first navigation azimuth is output as an initial azimuth value. An initial value calculation unit;
The amount of change in azimuth is calculated based on the angular velocity of the azimuth measured by the angular velocity sensor, and the movement is performed based on the calculated amount of change in azimuth and the initial value of azimuth output from the initial value calculator. The azimuth angle of the body is calculated as the second navigation azimuth angle, and based on the calculated second navigation azimuth angle , the three-dimensional acceleration measured by the acceleration sensor is converted into the three-dimensional acceleration of the navigation coordinate system, and converted. after integrating the three-dimensional acceleration to calculate the amount of change in the three-dimensional coordinates, calculated second navigation coordinate the coordinate values of the three-dimensional of the movable body by adding the amount of change in position initial value of the three-dimensional And a correction value for correcting the second navigation coordinate value using the second Kalman filter based on the calculated second navigation coordinate value and the three-dimensional coordinate value measured by the GPS receiver. And calculated The second navigation coordinate value is corrected using a positive value, and the computer is caused to function as an orientation value calculation unit that outputs the corrected second navigation coordinate value as a position orientation value obtained by locating the position of the moving body. Position location program.
JP2011056105A 2011-03-15 2011-03-15 POSITIONING DEVICE, POSITIONING METHOD OF POSITIONING DEVICE, AND POSITIONING PROGRAM Active JP5602070B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2011056105A JP5602070B2 (en) 2011-03-15 2011-03-15 POSITIONING DEVICE, POSITIONING METHOD OF POSITIONING DEVICE, AND POSITIONING PROGRAM

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2011056105A JP5602070B2 (en) 2011-03-15 2011-03-15 POSITIONING DEVICE, POSITIONING METHOD OF POSITIONING DEVICE, AND POSITIONING PROGRAM

Publications (2)

Publication Number Publication Date
JP2012193965A JP2012193965A (en) 2012-10-11
JP5602070B2 true JP5602070B2 (en) 2014-10-08

Family

ID=47086012

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011056105A Active JP5602070B2 (en) 2011-03-15 2011-03-15 POSITIONING DEVICE, POSITIONING METHOD OF POSITIONING DEVICE, AND POSITIONING PROGRAM

Country Status (1)

Country Link
JP (1) JP5602070B2 (en)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5896962B2 (en) * 2013-08-09 2016-03-30 本田技研工業株式会社 Sign information output device
EP3217321B1 (en) 2013-07-31 2019-08-28 Honda Motor Co., Ltd. Sign information output apparatus
CN105823481B (en) * 2015-12-21 2019-02-12 上海华测导航技术股份有限公司 A kind of GNSS-INS vehicle method for determining posture based on single antenna
CN105865450A (en) * 2016-04-19 2016-08-17 武汉理工大学 Zero-speed update method and system based on gait
CN111344590B (en) * 2018-01-30 2024-05-24 古野电气株式会社 Radar antenna device and azimuth measuring method
CN109855617A (en) * 2019-02-28 2019-06-07 深圳市元征科技股份有限公司 A kind of vehicle positioning method, vehicle locating device and terminal device
DE102019203332A1 (en) * 2019-03-12 2020-09-17 Robert Bosch Gmbh Position determination arrangement for a vehicle, vehicle
CN110440805B (en) * 2019-08-09 2021-09-21 深圳市道通智能航空技术股份有限公司 Method and device for fusing yaw angles and aircraft
CN111272170A (en) * 2020-03-17 2020-06-12 电子科技大学 Real-time pedestrian positioning system and method based on SoC
WO2022018964A1 (en) * 2020-07-20 2022-01-27 ソニーグループ株式会社 Information processing device, information processing method, and program
CN112363196B (en) * 2020-10-20 2023-10-31 北京嘀嘀无限科技发展有限公司 Vehicle attribute determining method, device, storage medium and electronic equipment
CN113758502B (en) * 2021-09-24 2024-02-20 广东汇天航空航天科技有限公司 Combined navigation processing method and device
WO2024180582A1 (en) * 2023-02-27 2024-09-06 三菱電機株式会社 Position/attitude calibration device, position/attitude calibration method, and position/attitude calibration program

Family Cites Families (37)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63151813A (en) * 1986-12-16 1988-06-24 Matsushita Electric Ind Co Ltd Apparatus for detecting position of moving body
JPS63302317A (en) * 1987-06-03 1988-12-09 Mitsubishi Heavy Ind Ltd Positional speed measuring apparatus of moving object
JPH04198787A (en) * 1990-11-28 1992-07-20 Matsushita Electric Ind Co Ltd Bearing measuring method for navigation system
JPH04213083A (en) * 1990-12-06 1992-08-04 Mazda Motor Corp Bearing detecting apparatus for vehicle
JPH0526681A (en) * 1991-07-25 1993-02-02 Nissan Motor Co Ltd Navigation device
JPH05157828A (en) * 1991-12-04 1993-06-25 Maspro Denkoh Corp Moving body position detecting device
JPH05215564A (en) * 1992-02-05 1993-08-24 Japan Aviation Electron Ind Ltd Automobile position measuring apparatus
ZA957639B (en) * 1994-10-24 1996-05-24 Caterpillar Inc System and method for precisely determining an operating point for an autonomous vehicle
JP3533745B2 (en) * 1995-03-30 2004-05-31 アイシン精機株式会社 Mobile positioning device
JP3354353B2 (en) * 1995-06-22 2002-12-09 防衛庁技術研究本部長 Adjustment calculation method during movement of inertial navigation system provided on flying vehicle
JPH095106A (en) * 1995-06-22 1997-01-10 Tamagawa Seiki Co Ltd True azimuth detection method of inertial apparatus
JPH09113293A (en) * 1995-10-17 1997-05-02 Japan Radio Co Ltd Vehicle bearing presuming system
US5991692A (en) * 1995-12-28 1999-11-23 Magellan Dis, Inc. Zero motion detection system for improved vehicle navigation system
JPH09189564A (en) * 1996-01-11 1997-07-22 Matsushita Electric Ind Co Ltd Traveling body position speed calculating device
JPH10176932A (en) * 1996-10-14 1998-06-30 Yokogawa Denshi Kiki Kk Inertial navigation apparatus
JP3753833B2 (en) * 1997-03-27 2006-03-08 アジア航測株式会社 Road linear automatic surveying equipment
JPH10307032A (en) * 1997-05-02 1998-11-17 Pioneer Electron Corp Navigator
JPH1194573A (en) * 1997-09-16 1999-04-09 Bio Oriented Technol Res Advancement Inst Position attitude measuring device for mobile body
JP3013309B1 (en) * 1999-02-19 2000-02-28 株式会社ゼンリン Hybrid traveling locus acquisition method and hybrid traveling locus acquisition system
ATE345487T1 (en) * 1999-09-16 2006-12-15 Sirf Tech Inc NAVIGATION SYSTEM AND METHOD FOR TRACKING THE POSITION OF AN OBJECT
US6577952B2 (en) * 2001-01-08 2003-06-10 Motorola, Inc. Position and heading error-correction method and apparatus for vehicle navigation systems
JP3914413B2 (en) * 2001-10-31 2007-05-16 多摩川精機株式会社 Initialization method for flying objects
US6826478B2 (en) * 2002-04-12 2004-11-30 Ensco, Inc. Inertial navigation system for mobile objects with constraints
US7346452B2 (en) * 2003-09-05 2008-03-18 Novatel, Inc. Inertial GPS navigation system using injected alignment data for the inertial system
JP2007156890A (en) * 2005-12-06 2007-06-21 Tamagawa Seiki Co Ltd Initial azimuth setting method for unmanned vehicle
JP4600357B2 (en) * 2006-06-21 2010-12-15 トヨタ自動車株式会社 Positioning device
US20080071476A1 (en) * 2006-09-19 2008-03-20 Takayuki Hoshizaki Vehicle dynamics conditioning method on MEMS based integrated INS/GPS vehicle navigation system
JP2008116370A (en) * 2006-11-06 2008-05-22 Toyota Motor Corp Mobile location positioning device
JP5398120B2 (en) * 2007-03-22 2014-01-29 古野電気株式会社 GPS combined navigation system
JP5022747B2 (en) * 2007-03-22 2012-09-12 古野電気株式会社 Mobile body posture and orientation detection device
JP2009236532A (en) * 2008-03-26 2009-10-15 Seiko Epson Corp Method for geolocation, program, and apparatus for geolocation
JP5322789B2 (en) * 2009-06-15 2013-10-23 三菱電機株式会社 Model generation apparatus, model generation method, model generation program, point cloud image generation method, and point cloud image generation program
JP5521416B2 (en) * 2009-07-08 2014-06-11 富士通株式会社 Autonomous positioning program, autonomous positioning device, and autonomous positioning method
JP5419665B2 (en) * 2009-12-10 2014-02-19 三菱電機株式会社 POSITION LOCATION DEVICE, POSITION LOCATION METHOD, POSITION LOCATION PROGRAM, Velocity Vector Calculation Device, Velocity Vector Calculation Method, and Velocity Vector Calculation Program
JP5586994B2 (en) * 2010-03-11 2014-09-10 三菱電機株式会社 POSITIONING DEVICE, POSITIONING METHOD OF POSITIONING DEVICE, AND POSITIONING PROGRAM
JP2011220844A (en) * 2010-04-09 2011-11-04 Seiko Epson Corp Position calculation method and position calculation device
JP5405417B2 (en) * 2010-08-24 2014-02-05 株式会社小野測器 Attitude angle stabilization device and method

Also Published As

Publication number Publication date
JP2012193965A (en) 2012-10-11

Similar Documents

Publication Publication Date Title
JP5602070B2 (en) POSITIONING DEVICE, POSITIONING METHOD OF POSITIONING DEVICE, AND POSITIONING PROGRAM
EP3109589B1 (en) A unit and method for improving positioning accuracy
JP5419665B2 (en) POSITION LOCATION DEVICE, POSITION LOCATION METHOD, POSITION LOCATION PROGRAM, Velocity Vector Calculation Device, Velocity Vector Calculation Method, and Velocity Vector Calculation Program
JP5328252B2 (en) Position detection apparatus and position detection method for navigation system
US9494428B2 (en) Attitude determination method, position calculation method, and attitude determination device
JP4964047B2 (en) Position detection apparatus and position detection method
JP4989035B2 (en) Error correction of inertial navigation system
US9026263B2 (en) Automotive navigation system and method to utilize internal geometry of sensor position with respect to rear wheel axis
CN104198765B (en) The coordinate system conversion method of vehicle acceleration of motion detection
JP5586994B2 (en) POSITIONING DEVICE, POSITIONING METHOD OF POSITIONING DEVICE, AND POSITIONING PROGRAM
CN111854740B (en) Inertial navigation system capable of dead reckoning in a vehicle
KR101394546B1 (en) Apparatus and method for position information acquisition using data fusion of gps and imu
JP6083279B2 (en) Movement status information calculation method and movement status information calculation device
JP2008058187A (en) Navigation apparatus, and method and program for position detection
JP2018047888A (en) System and method for measuring angular position of vehicle
US11408735B2 (en) Positioning system and positioning method
CN109186597A (en) A kind of localization method of the indoor wheeled robot based on double MEMS-IMU
CN106403952A (en) Method for measuring combined attitudes of Satcom on the move with low cost
KR100525517B1 (en) Car navigation system and control method thereof
JP2016033473A (en) Position calculation method and position calculation device
JP5164645B2 (en) Method and apparatus for repetitive calculation control in Kalman filter processing
CN105928515A (en) Navigation system for unmanned plane
KR20190116500A (en) In-vehicle, Computing Device and Program
JP2013228318A (en) Calibration quality determination apparatus and method
WO2010030565A1 (en) Magnetic sensing device for navigation and detecting inclination

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20130927

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20140425

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20140507

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140703

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: 20140722

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140819

R150 Certificate of patent or registration of utility model

Ref document number: 5602070

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250