JP2006226772A - Correction value parameter generating device, positioning device, and differential satellite positioning method - Google Patents
Correction value parameter generating device, positioning device, and differential satellite positioning method Download PDFInfo
- Publication number
- JP2006226772A JP2006226772A JP2005039237A JP2005039237A JP2006226772A JP 2006226772 A JP2006226772 A JP 2006226772A JP 2005039237 A JP2005039237 A JP 2005039237A JP 2005039237 A JP2005039237 A JP 2005039237A JP 2006226772 A JP2006226772 A JP 2006226772A
- Authority
- JP
- Japan
- Prior art keywords
- correction value
- positioning
- station
- mobile station
- satellite
- 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.)
- Withdrawn
Links
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
Description
本発明は、DGPS(Differential Global Positioning System)を代表とする差動衛星測位方式における補正値パラメータ生成装置,測位装置および差動衛星測位方法に関するものである。 The present invention relates to a correction value parameter generation device, a positioning device, and a differential satellite positioning method in a differential satellite positioning system represented by DGPS (Differential Global Positioning System).
差動衛星測位方式の測位計算は、移動局受信機が行う観測値を用いた衛星測位計算の精度を向上する方式である。基準局受信機の観測値を用いて補正値を生成し、この補正値を移動局受信機の観測値に適用するのが基本的な方法である。補正値の種類は、各測位信号から得られる観測値に対する個別の補正値や、測位計算結果の座標値に対する補正値などいくつかある。 The positioning calculation of the differential satellite positioning method is a method for improving the accuracy of the satellite positioning calculation using the observation values performed by the mobile station receiver. The basic method is to generate a correction value using the observation value of the reference station receiver and apply this correction value to the observation value of the mobile station receiver. There are several types of correction values, such as individual correction values for observation values obtained from each positioning signal and correction values for coordinate values of positioning calculation results.
補正値は、実時間で基準局から移動局に伝送する必要があるため、データ欠落に対応する必要がある。従って、該当時刻の時刻補正値(基準補正値)と、時間経過に対応した補正値時刻変化率を各時刻の補正値パラメータとして伝送する方法が一般的である。ある時刻に補正値が伝送されなかった場合は、過去に入手した時刻補正値に、補正値時刻変化率に経過時間を乗じた値を加算して補正値を生成する。
以下に、時刻tkの補正値パラメータである時刻補正値PRCと補正値時刻変化率RRCから時刻tの補正値PPRCを生成する式(式1)を示す。
PPRC(t) = PRC(tk)+RRC(tk)・(t−tk) (式1)
Since the correction value needs to be transmitted from the reference station to the mobile station in real time, it is necessary to cope with data loss. Therefore, it is common to transmit a time correction value (reference correction value) at a corresponding time and a correction value time change rate corresponding to the passage of time as a correction value parameter at each time. If the correction value is not transmitted at a certain time, a correction value is generated by adding a value obtained by multiplying the correction value time change rate by the elapsed time to the time correction value obtained in the past.
The following formula (Formula 1) for generating the correction value PPRC at time t from the time correction value PRC, which is the correction value parameter at time tk, and the correction value time change rate RRC is shown below.
PPRC (t) = PRC (tk) + RRC (tk) · (t−tk) (Formula 1)
差動衛星測位方式の測位方法は、移動局受信機と基準局受信機が同じ測位信号を同時に観測することで、基準局受信機の観測値に含まれる誤差から移動局受信機の観測値に含まれる誤差を推定して補償するのが原理である。
ここで、基準局受信機は基準局受信機の設置座標と測位衛星の衛星軌道が既知である。従って、衛星軌道から測位衛星の座標を算出し、この既知である座標から求まる測位衛星と基準局受信機との距離と、基準局受信機の観測値から求まる測位衛星と基準局受信機との距離とを比較することで、基準局受信機の観測値に含まれる誤差を検出できる。
上記の原理から補正値を用いて高精度の測位を実現する為には移動局と基準局とが近接していることが必要になる。移動局と基準局との受信機間が離れていると、移動局と基準局とで測位信号の伝播経路が大きく異なり、互いの観測値に含まれる誤差成分の違いが大きく変化する為に、移動局受信機の観測値の誤差と基準局受信機の観測値の誤差とが一致しなくなるからである。
差動衛星測位方式用補正データの算出方式についての文献に特許文献1がある。
Here, the reference station receiver knows the installation coordinates of the reference station receiver and the satellite orbit of the positioning satellite. Therefore, the coordinates of the positioning satellite are calculated from the satellite orbit, and the distance between the positioning satellite obtained from the known coordinates and the reference station receiver and the distance between the positioning satellite and the reference station receiver obtained from the observation value of the reference station receiver are calculated. By comparing, the error included in the observation value of the reference station receiver can be detected.
Based on the above principle, it is necessary that the mobile station and the reference station be close to each other in order to achieve highly accurate positioning using a correction value. If the receivers of the mobile station and the reference station are separated from each other, the propagation path of the positioning signal differs greatly between the mobile station and the reference station, and the difference in the error component included in each observation value changes greatly. This is because the error in the observation value of the receiver does not match the error in the observation value of the reference station receiver.
There is
従って、ある広い区域内で適切な補正値を供給する為には、複数の基準局を区域内に適切に配置する必要がある。そして移動局は、その区域内を移動する際に適切な基準局を順次選択して補正値を入手する必要がある。 Accordingly, in order to supply an appropriate correction value within a certain wide area, it is necessary to appropriately arrange a plurality of reference stations within the area. When the mobile station moves within the area, it is necessary to sequentially select appropriate reference stations to obtain correction values.
しかし、要求精度を達成する為に配置すべき基準局数が、様々な制限により配置できないために、そのままでは要求精度が実現できない場合がある。
また、ある要求仕様の補正値を生成するように配置した基準局を、より高精度な仕様の補正値を生成する目的で使用したい場合がある。
また、移動しながら差動衛星測位方式による測位計算をする場合に、使用する基準局を切替えた際に、補正値の急な変動が発生して測位結果に大きな影響を与える場合がある。
However, since the number of reference stations to be arranged to achieve the required accuracy cannot be arranged due to various restrictions, the required accuracy may not be realized as it is.
In some cases, a reference station arranged to generate a correction value of a certain required specification may be used for the purpose of generating a correction value of a more accurate specification.
In addition, when performing the positioning calculation by the differential satellite positioning method while moving, when the reference station to be used is switched, a sudden change of the correction value may occur, which may greatly affect the positioning result.
本発明は、例えば上記の課題を解決するためになされたもので、差動衛星測位方式による測位の精度を向上することを目的とする。 The present invention has been made, for example, in order to solve the above-described problems, and it is an object of the present invention to improve positioning accuracy by a differential satellite positioning method.
本発明の補正値パラメータ生成装置は、各測位衛星からの測位信号を観測した観測値を基準局から受信する観測値受信部と、前記観測値受信部が入力した観測値に基づいて各測位衛星毎に基準補正値を生成する基準補正値生成部と、単位時間あたりの基準補正値の変化量を示す補正値時刻変化率と、測位対象である移動局と前記基準局との座標差の単位座標あたりの基準補正値の変化量を示す補正値座標変化率との少なくともいずれかの変化率を、前記基準局から各測位衛星に向けた各仰角に基づいて各測位衛星毎に算出する変化率算出部と、前記移動局の座標を測位する際の補正処理に使用する補正値パラメータとして、前記基準補正値生成部が生成した各測位衛星毎の基準補正値と前記変化率算出部が算出した各測位衛星毎の変化率とを前記移動局に送信する補正値パラメータ送信部とを備えたことを特徴とする。 The correction value parameter generation apparatus of the present invention includes an observation value receiving unit that receives an observation value obtained by observing a positioning signal from each positioning satellite from a reference station, and each positioning satellite based on the observation value input by the observation value receiving unit. Per unit coordinate of a reference correction value generator for generating a reference correction value, a correction value time change rate indicating a change amount of the reference correction value per unit time, and a coordinate difference between the mobile station to be positioned and the reference station A rate-of-change calculation unit that calculates a change rate of at least one of a correction value coordinate change rate indicating a change amount of the reference correction value for each positioning satellite based on each elevation angle from the reference station toward each positioning satellite; As a correction value parameter used for correction processing when positioning the coordinates of the mobile station, a reference correction value for each positioning satellite generated by the reference correction value generation unit and each positioning satellite calculated by the change rate calculation unit The rate of change for each Characterized in that a correction value parameter transmission unit which transmits the serial mobile station.
本発明によれば、基準局から測位衛星に向けた仰角に基づく補正値時刻変化率と補正値座標変化率とを算出することで、測位結果に誤差を生じさせる電離層や対流圏の対処をした補正値を算出することができ、精度の高い測位結果を得ることができる。 According to the present invention, by calculating the correction value time change rate and the correction value coordinate change rate based on the elevation angle from the reference station toward the positioning satellite, a correction value for dealing with the ionosphere and troposphere that causes an error in the positioning result Can be calculated, and a highly accurate positioning result can be obtained.
実施の形態1.
実施の形態1において、差動衛星測位方式の測位システムは、測位対象であり測位処理を行う移動局と、1以上の測位衛星の発信した測位信号を観測する複数の基準局と、基準補正値を生成する基準補正値生成局とを備える。
図1は、実施の形態1における差動衛星測位方法を示すフローチャートである。
差動衛星測位方法の1例について図1に基づいて以下に説明する。
In the first embodiment, a differential satellite positioning type positioning system includes a mobile station that is a positioning target and that performs positioning processing, a plurality of reference stations that observe positioning signals transmitted from one or more positioning satellites, and a reference correction value. A reference correction value generation station to be generated.
FIG. 1 is a flowchart showing a differential satellite positioning method in the first embodiment.
An example of the differential satellite positioning method will be described below with reference to FIG.
まず、測位対象である移動局の概算座標を算出する(S101)。
ここで、移動局の概算座標には、例えば、補正値を用いずに移動局が、GPS衛星やガリレオ衛星などの測位衛星を用いた衛星測位や慣性航法装置を用いた慣性測位などにより、測位した結果の座標や、前回測位した時の座標などを利用する。
First, approximate coordinates of a mobile station that is a positioning target are calculated (S101).
Here, in the approximate coordinates of the mobile station, for example, the mobile station can perform positioning without using a correction value by satellite positioning using a positioning satellite such as a GPS satellite or a Galileo satellite or inertial positioning using an inertial navigation device. The coordinates of the result obtained and the coordinates at the previous positioning are used.
また、各基準局で、各測位衛星が発信した測位信号を受信し、受信した測位信号から得られる情報(以下、観測値とする)を取得する。そして、各基準局で取得された観測値を基準補正値生成局に収集する(S102)。
ここで、観測値には、例えば、測位信号を受信した局と測位信号を発信した測位衛星との距離を示す擬似距離や測位衛星の軌道情報などが含まれる。
Each reference station receives a positioning signal transmitted from each positioning satellite, and acquires information (hereinafter referred to as an observation value) obtained from the received positioning signal. Then, the observation values acquired at each reference station are collected at the reference correction value generation station (S102).
Here, the observed value includes, for example, a pseudo distance indicating the distance between the station that has received the positioning signal and the positioning satellite that has transmitted the positioning signal, orbit information of the positioning satellite, and the like.
次に、基準補正値生成局は、各基準局から収集した観測値に基づいて、各基準局毎の補正値の変化率を算出する(S103)。
ここで、算出する補正値の変化率には補正値時刻変化率と補正値座標変化率とがある。
Next, the reference correction value generation station calculates the change rate of the correction value for each reference station based on the observation values collected from each reference station (S103).
Here, the calculated change rate of the correction value includes a correction value time change rate and a correction value coordinate change rate.
実施の形態1における補正値時刻変化率とは、基準局から測位衛星に向けた仰角に基づく単位時間あたりの補正値変化量であり、基準局毎に各測位衛星に対応して算出する。
測位計算で算出される擬似距離は測位信号の搬送時間に光速を乗じることで求められ、擬似距離に誤差を生じさせる大きな要因が電離層・対流圏を通過時の測位信号の速度遅延である。
そこで、基準局から測位衛星に向けた仰角に基づいて電離層・対流圏の状態と測位信号が通過する電離層・対流圏の距離との関係から補正値時刻変化率を算出する。
The correction value time change rate in the first embodiment is a correction value change amount per unit time based on the elevation angle from the reference station toward the positioning satellite, and is calculated for each positioning satellite for each reference station.
The pseudo distance calculated by the positioning calculation is obtained by multiplying the transport time of the positioning signal by the speed of light, and a major factor causing an error in the pseudo distance is the speed delay of the positioning signal when passing through the ionosphere / troposphere.
Therefore, the correction value time change rate is calculated from the relationship between the ionosphere / troposphere state and the ionosphere / troposphere distance through which the positioning signal passes, based on the elevation angle from the reference station toward the positioning satellite.
測位信号を受信する局から衛星への方向を示す視線ベクタの仰角が低い場合は、電離層や対流圏などを通過する距離が長くなる。つまり、低仰角になるほど、電離層の影響が大きくなる。また影響の増加の割合も低仰角になるほど大きくなる。このため、測位衛星が仰角を下げつつ運動している場合、補正値時刻変化率も大きくなるのが適正である。逆も成立する。また厳密には、衛星の移動に伴う仰角の変化を考慮する必要があるため、仰角と仰角変化率を考慮して補正値時刻変化率を算出することが望ましい。さらに、電離層・対流圏の状態は日々変化し、電離層・対流圏の状態により遅延量も変わるため、電離層・対流圏の状態も考慮することが望ましい。
従来の補正値時刻変化率は、前回算出した補正値と今回算出した補正値との差分を変化率で表すだけであるため、長時間の予測計算が可能な信頼性を持たないことが多い。このため、仰角、仰角変化率、電離層・対流圏の状態を考慮することで補正値時刻変化率の品質が上がり、測位精度が改善する。
When the elevation angle of the line-of-sight vector indicating the direction from the station that receives the positioning signal to the satellite is low, the distance passing through the ionosphere, the troposphere, or the like becomes long. That is, the lower the elevation angle, the greater the influence of the ionosphere. Also, the rate of increase in the effect increases as the angle of elevation decreases. For this reason, when the positioning satellite moves while lowering the elevation angle, it is appropriate that the correction value time change rate also increases. The reverse is also true. Strictly speaking, since it is necessary to consider the change in the elevation angle accompanying the movement of the satellite, it is desirable to calculate the correction value time change rate in consideration of the elevation angle and the elevation angle change rate. Furthermore, since the ionosphere / troposphere state changes daily and the delay varies depending on the ionosphere / troposphere state, it is desirable to consider the ionosphere / troposphere state.
Since the conventional correction value time change rate merely represents the difference between the correction value calculated last time and the correction value calculated this time as the change rate, it often does not have the reliability that enables long-term prediction calculation. For this reason, the quality of the correction value time change rate is improved and the positioning accuracy is improved by considering the elevation angle, the elevation angle change rate, and the ionosphere / troposphere state.
図2は、実施の形態1における仰角に対する遅延量の係数を示す図である。
任意の範囲の仰角において、仰角が最大(x=xmax2)の時に係数(W(x))を小さくする。つまり、仰角が最大の時に電離層・対流圏による遅延の影響が一番少ないことを示す。また、仰角が最小(x=xmin1)の時に係数(W(x))を大きくする。つまり、仰角が最小の時に電離層・対流圏による遅延の影響が一番多いことを示す。
FIG. 2 is a diagram showing a coefficient of delay amount with respect to the elevation angle in the first embodiment.
The coefficient (W (x)) is decreased when the elevation angle is the maximum (x = xmax2) in an arbitrary range of elevation angles. In other words, when the elevation angle is maximum, the effect of delay due to the ionosphere and troposphere is the least. Further, the coefficient (W (x)) is increased when the elevation angle is minimum (x = xmin1). In other words, when the elevation angle is at a minimum, the influence of the delay due to the ionosphere / troposphere is the largest.
同様に、電離層・対流圏の状態に対する遅延量の係数と仰角変化率に対する遅延量の係数を求め、各係数に基づいて補正値時刻変化率を算出する。各係数を使用して補正値時刻変化率を算出するためのモデル式は任意である。
ここで、仰角と仰角変化率、電離層・対流圏の状態の少なくともいずれかとに基づいて算出する補正値時刻変化率を、仰角に基づく補正値時刻変化率とする。
Similarly, a delay coefficient for the ionosphere / troposphere state and a delay coefficient for the elevation angle change rate are obtained, and a correction value time change rate is calculated based on each coefficient. The model formula for calculating the correction value time change rate using each coefficient is arbitrary.
Here, the correction value time change rate calculated based on at least one of the elevation angle, the elevation change rate, and the ionosphere / troposphere state is set as a correction value time change rate based on the elevation angle.
次に、図1のS103で算出する変化率の一つである補正値座標変化率について説明する。
実施の形態1における補正値座標変化率とは、基準局と移動局との座標値の単位差あたりの補正値変化量であり、各基準局の補正値に対応して算出する。
つまり、基準局の座標値と移動局の座標値とに対する補正値変化量は、既知である基準局の座標からS101で算出した移動局の概算座標へのベクタの緯度および経度方向の成分に、各成分に対する補正値座標変化率を乗じた値である。
各衛星に対応する補正値座標変化率の算出手順を説明する。既知である基準局の座標からS102で収集した観測値などから算出できる測位衛星の座標へのベクタ(以下、衛星視線ベクタとする)の各成分に対して、仰角、仰角変化率、電離層・対流圏の状態の受信機座標値に依存する変化率に基づいて補正値座標変化率を算出すればよい。また、基準方向から衛星が見える方向への角度を示す方位角(azimuth)に基づいて、衛星視線ベクタの各成分に対する補正値座標変化率を算出してもよい。
図3は、実施の形態1における仰角に対する補正値座標変化率の算出時の係数を示す図である。横軸が仰角で縦軸が係数である。
図4は、実施の形態1における衛星視線ベクタを示す図である。図4は、北を基準方向として時計回りに方位角を表した場合を示している。
Next, the correction value coordinate change rate that is one of the change rates calculated in S103 of FIG. 1 will be described.
The correction value coordinate change rate in the first embodiment is a correction value change amount per unit difference in coordinate values between the reference station and the mobile station, and is calculated corresponding to the correction value of each reference station.
That is, the correction value change amount with respect to the coordinate value of the reference station and the coordinate value of the mobile station is the component of the latitude and longitude directions of the vector from the known coordinate of the reference station to the approximate coordinate of the mobile station calculated in S101. Is a value obtained by multiplying the correction value coordinate change rate for.
A procedure for calculating the correction value coordinate change rate corresponding to each satellite will be described. With respect to each component of the positioning satellite coordinates that can be calculated from the observation values collected in S102 from the coordinates of the known reference station (hereinafter referred to as satellite line-of-sight vectors), the elevation angle, rate of elevation change, ionosphere / troposphere The correction value coordinate change rate may be calculated based on the change rate depending on the receiver coordinate value of the state. Further, the correction value coordinate change rate for each component of the satellite line-of-sight vector may be calculated based on an azimuth indicating an angle from the reference direction to the direction in which the satellite can be seen.
FIG. 3 is a diagram showing coefficients when calculating the correction value coordinate change rate with respect to the elevation angle in the first embodiment. The horizontal axis is the elevation angle and the vertical axis is the coefficient.
FIG. 4 is a diagram showing satellite line-of-sight vectors in the first embodiment. FIG. 4 shows a case where the azimuth is represented clockwise with north as a reference direction.
次に、補正値生成装置は、補正値パラメータの一つである、基準補正値を生成する(S104)。
ここで、基準補正値とは、例えば、既知である基準局と測位衛星との座標から求まる距離と測位信号の搬送時間に光速を乗じて求まる距離との差分や、既知である基準局の座標と基準局が受信した測位信号から求まる基準局の座標との差分である。
Next, the correction value generation device generates a reference correction value that is one of the correction value parameters (S104).
Here, the reference correction value is, for example, the difference between the distance obtained from the coordinates of the known reference station and the positioning satellite and the distance obtained by multiplying the transport time of the positioning signal by the speed of light, the coordinates of the known reference station and the reference station Is the difference from the coordinates of the reference station obtained from the received positioning signal.
次に、基準補正値生成局は、S101で算出した移動局の概算座標と既知である基準局の座標とに基づいて、複合する補正値に対応する基準局を選択する(S105)。 Next, the reference correction value generation station selects a reference station corresponding to the combined correction value based on the approximate coordinates of the mobile station calculated in S101 and the coordinates of the known reference station (S105).
図5、図6は実施の形態1における移動局と基準局との距離の関係を示す図である。
図5および図6は、設定した閾値に対応する、選択する基準局の移動局からの距離を示している。
閾値の設定方法には、例えば、移動局に最も近い基準局までの距離に対応して設定する方法と要求精度に対応して設定する方法とがある。
移動局に最も近い基準局までの距離に対応して設定する方法では、移動局から最も近い基準局までの距離が短い場合に閾値を高く設定するとよい。移動局の近くに基準局がある場合、この基準局の観測値に基づいて生成した補正値を利用することで、移動局の観測値に基づく測位を高精度に補正することができるためである。そこで、閾値を高くして選択する基準局の移動局からの距離を短くする。そして、選択する基準局を減らし、他の基準局の補正値が複合補正値に寄与する率を下げるように制御する。
また、要求精度に対応して設定する方法では、要求精度が高い場合に閾値を低く設定するとよい。閾値を低く設定することで選択する基準局の移動局からの距離を長くする。そして、選択する基準局を増やす。移動しながら測位計算する場合、移動局から各基準局までの距離が変わることにより補正値を複合する基準局の切替えが発生するが、より多くの基準局の補正値を複合することで、基準局の切替えによる複合補正値の変動を抑える効果があるためである。これにより差動衛星測位方式による測位において、連続する測位タイミングで大きく隔たる座標を得るような、不連続な結果を発生しないようにできる。
ただし、基準局を1局以上選択できるように閾値は設定する。
また、図5および図6は閾値と距離との関係を示す非線形関数の一例である。また、線形関数で閾値と距離との関係が示されてもよい。
5 and 6 are diagrams illustrating the relationship between the distance between the mobile station and the reference station in the first embodiment.
5 and 6 show the distance from the mobile station of the reference station to be selected corresponding to the set threshold value.
The threshold setting method includes, for example, a method of setting corresponding to the distance to the reference station closest to the mobile station and a method of setting corresponding to the required accuracy.
In the method of setting corresponding to the distance from the mobile station to the nearest reference station, the threshold value may be set higher when the distance from the mobile station to the nearest reference station is short. This is because, when there is a reference station near the mobile station, the positioning based on the observation value of the mobile station can be corrected with high accuracy by using the correction value generated based on the observation value of the reference station. Therefore, the distance from the mobile station to the selected reference station is shortened by increasing the threshold. Then, control is performed so that the number of reference stations to be selected is reduced and the rate at which the correction values of other reference stations contribute to the composite correction value is decreased.
In the method of setting corresponding to the required accuracy, the threshold value may be set low when the required accuracy is high. The distance from the mobile station of the reference station to be selected is increased by setting the threshold value low. Then, the number of reference stations to be selected is increased. When performing positioning calculation while moving, switching of the reference station that combines the correction values occurs due to the change in the distance from the mobile station to each reference station, but by combining the correction values of more reference stations, This is because there is an effect of suppressing the fluctuation of the composite correction value. Thereby, in the positioning by the differential satellite positioning method, it is possible to prevent the generation of discontinuous results such as obtaining coordinates largely separated at successive positioning timings.
However, the threshold is set so that one or more reference stations can be selected.
FIG. 5 and FIG. 6 are examples of nonlinear functions showing the relationship between the threshold and the distance. Further, the relationship between the threshold value and the distance may be indicated by a linear function.
次に、基準補正値生成局および移動局は、基準局毎に、S103で算出した変化率から補正値変化量を算出し、S104で生成した基準補正値に算出した補正値変化量を加算して、移動局の測位の際の補正処理に使用する補正値(以下、移動局用補正値とする)を生成する(S106)。
例えば、S103で、基準局の座標から移動局の概算座標へのベクタの緯度および経度方向に対する補正値座標変化率を算出した場合の補正値変化量DPRCは以下の式(式2)で算出できる。
Next, the reference correction value generation station and the mobile station calculate the correction value change amount from the change rate calculated in S103 for each reference station, and add the calculated correction value change amount to the reference correction value generated in S104. Then, a correction value (hereinafter, referred to as a mobile station correction value) used for correction processing at the time of positioning of the mobile station is generated (S106).
For example, the correction value change amount DPRC when the correction value coordinate change rate with respect to the latitude and longitude directions of the vector from the coordinates of the reference station to the approximate coordinates of the mobile station is calculated in S103 can be calculated by the following equation (Equation 2).
DPRC = DPRRCφ・dφ+DPRRCλ・dλ (式2)
DPRC:補正値変化量
DPRRCφ:補正値座標変化率の緯度方向成分
DPRRCλ:補正値座標変化率の経度方向成分
dφ:基準局座標と移動局概算座標の緯度差
dλ:基準局座標と移動局概算座標の経度差
DPRC = DPRRCφ · dφ + DPRRCλ · dλ (Formula 2)
DPRC: Correction value change amount DPRRCφ: Latitude direction component of correction value coordinate change rate DPRRCλ: Longitude direction component of correction value coordinate change rate dφ: Latitude difference between reference station coordinates and mobile station approximate coordinates dλ: Reference station coordinates and mobile station approximate coordinates Longitude difference
そして、各基準局の各測位衛星に対する移動局用補正値PPRCqは以下の式(式3)で算出できる。 And the correction value PPRCq for mobile stations for each positioning satellite of each reference station can be calculated by the following equation (Equation 3).
PPRCq=PRCq+RRCq・dt+DPRCq(dφ、dλ、Eq) (式3)
PPRCq:衛星qに対応する補正値
PRCq:衛星qに対応する時刻補正値
RRCq:衛星qに対応する補正値時刻変化率
dt:時刻補正値PRCqに対する時刻差
DPRCq:補正値座標変化率で算出される衛星qに対応する補正値変化量。
dφ:基準局座標と移動局概算座標の緯度差
dλ:基準局座標と移動局概算座標の経度差
Eq:移動局から見た衛星qの仰角
PPRCq = PRCq + RRCq · dt + DPRCq (dφ, dλ, Eq) (Formula 3)
PPRCq: correction value corresponding to satellite q PRCq: time correction value corresponding to satellite q RRCq: correction value time change rate corresponding to satellite q dt: time difference with respect to time correction value PRCq DPRCq: correction value coordinate change rate Correction value change amount corresponding to the satellite q.
dφ: Latitude difference between reference station coordinates and mobile station approximate coordinates dλ: Longitude difference between reference station coordinates and mobile station approximate coordinates Eq: Elevation angle of satellite q as seen from mobile station
上記式3において、RRCqには上記従来の補正値時刻変化率や上記仰角に基づく補正値時刻変化率を用いる。
なお、式2、3に示す補正値は、GPS信号(測位信号)のC/Aコード等を解析して得られる擬似距離の補正値や、GPS信号(測位信号)の搬送波位相擬似距離の補正値を想定しているが、測位結果である位置の補正値であっても良い。
In Equation 3, the conventional correction value time change rate or the correction value time change rate based on the elevation angle is used as RRCq.
Note that the correction values shown in Equations 2 and 3 are correction values for pseudoranges obtained by analyzing C / A codes of GPS signals (positioning signals), and corrections for carrier phase pseudoranges for GPS signals (positioning signals). Although a value is assumed, it may be a position correction value as a positioning result.
また、補正値時刻変化率を使用しない場合の移動局用補正値PCPCqは以下の式(式4)で算出できる。 Further, the correction value PCPCq for the mobile station when the correction value time change rate is not used can be calculated by the following equation (Equation 4).
PCPCq=CPCq+DCPCq(dφ、dλ、Eq) (式4)
PCPCq:衛星qに対応する搬送波位相擬似距離補正値
CPCq:基準局で通常のRTK測位計算用に生成された搬送波位相擬似距離補正値
DCPCq:基準局と移動局の座標を用いて生成される座標依存の補正値変化量
RTK:Real Time Kinematic
PCPCq = CPCq + DCPCq (dφ, dλ, Eq) (Formula 4)
PCPCq: Carrier phase pseudorange correction value corresponding to satellite q CPCq: Carrier phase pseudorange correction value generated for normal RTK positioning calculation at the reference station DCPCq: Coordinate-dependent generated using the coordinates of the reference station and mobile station Correction value change amount RTK: Real Time Kinetic
上記式4は、FKP(面補正パラメータ)方式に対応して搬送波位相擬似距離に対する補正値を算出する場合を示している。 Equation 4 above shows a case where the correction value for the carrier phase pseudorange is calculated in accordance with the FKP (surface correction parameter) method.
次に、基準補正値生成局および移動局は、S105で選択した基準局毎に、S106で生成した移動局用補正値を、複合した複合移動局用補正値を算出する(S107)。
補正値を複合する方法には、例えば、各基準局毎に得られる移動局対応の補正値を平均する方法、移動局からの距離に基づいて各基準局毎の移動局対応の補正値に重み付けを行って重み付けした補正値を所定数(少なくとも2以上)の基準局について合算する方法、測位衛星に対する仰角に基づいて各基準局に対応する補正値に重み付けを行って重み付けした補正値を合算する方法がある。
距離に基づいて補正値を重み付けする方法では、図5や図6で示した非線形関数に基づいて重み(W(x))を設定する。基準局の配置間隔に対応したこのような非線形関数に基づく重みを設定することで、移動局の位置に依存する基準局切替え時の複合補正値の変動を抑える効果がある。また、その他の任意の関数に基づいて重みを設定してもよい。
各基準局に対応する補正値を複合して、補正処理に使用する補正値を生成することができることで、移動局に対して基準局が必ずしも適切に配置されていない場合であっても、精度の高い補正値を適用し、精度の高い測位結果を得ることができる。
特に、移動局と各基準局との距離に基づいて補正値を複合することで、移動局に対して適切に配置された基準局がある場合に、その適切に配置された基準局の補正値を重視した精度の高い補正値を適用することができる。また、補正値を複合する移動局の切替わりによる影響を抑えることができる。
また、各基準局から各測位衛星に向けた各仰角に基づいて補正値を複合することで、電離層・対流圏による影響を考慮した精度の高い補正値を適用することができる。
Next, the reference correction value generation station and the mobile station calculate, for each reference station selected in S105, a composite mobile station correction value obtained by combining the mobile station correction value generated in S106 (S107).
The method of combining correction values includes, for example, a method of averaging correction values corresponding to mobile stations obtained for each reference station, and weighting correction values corresponding to mobile stations for each reference station based on the distance from the mobile station. There are a method of summing the weighted correction values for a predetermined number (at least 2) of reference stations, and a method of weighting the correction values corresponding to the respective reference stations based on the elevation angle with respect to the positioning satellite and summing the weighted correction values.
In the method of weighting the correction value based on the distance, the weight (W (x)) is set based on the nonlinear function shown in FIGS. By setting a weight based on such a nonlinear function corresponding to the arrangement interval of the reference stations, there is an effect of suppressing the fluctuation of the composite correction value at the time of switching the reference station depending on the position of the mobile station. Further, the weight may be set based on any other function.
The correction value used for the correction process can be generated by combining the correction values corresponding to each reference station, so that even when the reference station is not necessarily properly arranged with respect to the mobile station, the accuracy is high. By applying the correction value, a highly accurate positioning result can be obtained.
In particular, by combining the correction values based on the distance between the mobile station and each reference station, when there is a reference station appropriately arranged with respect to the mobile station, the correction value of the appropriately arranged reference station is emphasized A highly accurate correction value can be applied. In addition, it is possible to suppress the influence of switching of mobile stations that combine correction values.
Further, by combining the correction values based on the respective elevation angles from the respective reference stations to the respective positioning satellites, it is possible to apply a high-accuracy correction value in consideration of the influence of the ionosphere / troposphere.
そして、移動局の観測値に基づく測位計算をS107で算出した複合移動局用補正値で補正して移動局の位置を測位する(S108)。 Then, the positioning calculation based on the observation value of the mobile station is corrected with the composite mobile station correction value calculated in S107, and the position of the mobile station is determined (S108).
上記に説明したS101〜S108の処理により移動局の測位処理において精度の高い補正値(複合移動局用補正値)を適用し、精度の高い測位結果を得ることができる。 A highly accurate correction value (compound mobile station correction value) is applied in the positioning process of the mobile station by the processes of S101 to S108 described above, and a highly accurate positioning result can be obtained.
実施の形態2.
実施の形態2では、上記実施の形態1で説明した測位方法を実行するシステムの一例について説明する。
移動局300は携帯端末や車載装置などの形態であり、測位を行う測位装置である。基準補正値生成局200は補正機能、補正情報を提供するサーバ装置である。
Embodiment 2. FIG.
In the second embodiment, an example of a system that executes the positioning method described in the first embodiment will be described.
The
図7は、実施の形態2における測位システムの構成図である。
図7において、測位システムは、基準局100、基準補正値生成局200、移動局300を備え、測位衛星400を利用する。
基準局100は受信機110と通信機120とを備え、基準補正値生成局200は収集局210と生成局220と通信局230とを備え、移動局300は受信機310と複合部320と通信機330とを備える。
FIG. 7 is a configuration diagram of the positioning system in the second embodiment.
In FIG. 7, the positioning system includes a
The
図8は、実施の形態2における測位処理のフローチャートである。 FIG. 8 is a flowchart of the positioning process in the second embodiment.
まず、基準局100の受信機110は、測位衛星400が発信する測位信号を受信する(S201)。
First, the
次に、基準補正値生成局200の収集局210は、各基準局100から観測値を収集する。この時、各基準局100では通信機120が受信機110の受信した測位信号に基づく観測値を基準補正値生成局200に送信する(S202)。この処理は、上記実施の形態1で説明した図1のS102に対応する。
次に、生成局220は、収集局210がS202の処理で収集した各基準局100の観測値に基づいて補正値時刻変化率および補正値座標変化率を算出する(S203)。この処理は、上記実施の形態1で説明した図1のS103に対応する。
次に、生成局220は、収集局210がS202の処理で収集した各基準局100の観測値と既知である各基準局100の座標とに基づいて基準補正値を生成する(S204)。この処理は、上記実施の形態1で説明した図1のS104に対応する。
次に、通信局230は、各基準局100と測位衛星400と毎に補正値パラメータである、基準補正値、補正値時刻変化率、補正値座標変化率を移動局300に配信する(S205)。配信方法は、特定の移動局への配信でもよいし、放送方式による全移動局への配信でもよい。
Next, the
Next, the
Next, the
Next, the
次に、移動局300の通信機330は、基準補正値生成局200がS205の処理で配信した補正値パラメータを受信する(S206)。
次に、複合部320は、移動局300の概算座標を算出する(S207)。この処理は、上記実施の形態1で説明した図1のS101に対応する。
次に、複合部320は、移動局300の概算座標と既知である基準局100の座標とに基づいて複合する補正値に対応する基準局100を選択する(S208)。この処理は、上記実施の形態1で説明した図1のS105に対応する。
次に、複合部320は、受信した補正値パラメータを基に移動局用補正値を生成する(S209)。この処理は、上記実施の携帯1で説明した図1のS106に対応する。
次に、複合部320は、S208の処理で選択した基準局毎に、S209の処理で生成した移動局用補正値を複合し複合移動局用補正値を算出する(S210)。この処理は、上記実施の形態1で説明した図1のS107に対応する。
そして、受信機310は測位衛星400が発信する測位信号を受信し、受信した測位信号に基づく観測値と、S210の処理で算出した複合移動局用補正値とに基づいて測位する(S211)。この処理は、上記実施の形態1で説明した図1のS108に対応する。
Next, the
Next, the
Next, the
Next, the
Next, the
Then, the
上記に説明した測位システムの各構成により移動局300の測位処理において精度の高い補正値を適用し、精度の高い測位結果を得ることができる。
また、移動局300側で補正値の複合処理を行うことで、基準補正値生成局200の負荷を軽減することができる。
With each configuration of the positioning system described above, a highly accurate correction value can be applied in the positioning process of the
In addition, by performing correction value combination processing on the
実施の形態3.
実施の形態3では、上記実施の形態1で説明した測位方法を実行するシステムの一例について説明する。
図9は、実施の形態3における測位処理のフローチャートである。
上記実施の形態2と異なる形態を図9に基づいて以下に説明する。なお、測位システムの構成は上記実施の形態2で説明した図7と同様である。
Embodiment 3 FIG.
In the third embodiment, an example of a system that executes the positioning method described in the first embodiment will be described.
FIG. 9 is a flowchart of the positioning process in the third embodiment.
A different form from the said Embodiment 2 is demonstrated below based on FIG. The configuration of the positioning system is the same as that shown in FIG. 7 described in the second embodiment.
まず、移動局300の複合部320は、移動局300の概算座標を算出する(S301)。この処理は、上記実施の形態1で説明した図1のS101に対応する。
次に、通信機330は、複合部320がS301の処理で算出した概算座標を基準補正値生成局200に送信する(S302)。
First, the
Next, the
また、基準局100の受信機110は、測位衛星400が発信する測位信号を受信する(S303)。
The
そして、基準補正値生成局200の通信局230は、移動局300がS302の処理で送信した移動局300の概算座標を受信する(S304)。
また、収集局210は、各基準局100から観測値を収集する。この時、各基準局100では通信機120が受信機110の受信した測位信号に基づく観測値を基準補正値生成局200に送信する(S305)。この処理は、上記実施の形態1で説明した図1のS102に対応する。
次に、生成局220は、通信局230が受信した移動局300の概算座標と既知である各基準局100の座標とに基づいて、移動局300の周辺に配置されている基準局100を複数選択する(S306)。この処理は、上記実施の形態1で説明した図1のS105に対応する。
次に、生成局220は、通信局230がS304の処理で受信した移動局300の概算座標と収集局210がS305の処理で収集した各基準局100の観測値とに基づいて、S306の処理で選択した各基準局100について補正値時刻変化率および補正値座標変化率を算出する(S307)。この処理は、上記実施の形態1で説明した図1のS103に対応する。
次に、生成局220は、収集局210がS305の処理で収集した各基準局100の観測値と既知である各基準局100の座標とに基づいて、S306の処理で選択した各基準局100について基準補正値を生成する(S308)。この処理は、上記実施の形態1で説明した図1のS104に対応する。
次に、生成局220は、S308の処理で生成した基準補正値に、S307の処理で算出した補正値座標変化率に基づく補正値変化量を加算して、補正値パラメータの一つである、基準補正値を更新する(S309)。この処理は、上記実施の形態1で説明した図1のS106に部分的に対応する。
次に、通信局230は、生成局220がS309の処理で生成した各基準局100と各測位衛星400と毎に補正値パラメータである、基準補正値、補正値時刻変化率を移動局300に送信する(S310)。
Then, the
The
Next, the
Next, the generating
Next, the
Next, the
Next, the
次に、移動局300の通信機330は、基準補正値生成局200がS310の処理で送信した補正値パラメータを受信する(S311)。
次に、複合部320は、S301の処理で算出した移動局300の概算座標と既知である各基準局100の座標とに基づいて複合する補正値に対応する基準局100を選択する(S312)。この処理は、上記実施の形態1で説明した図1のS106に対応する。
次に複合部320は、受信した補正値パラメータを基に、基準局毎に移動体用補正値を生成する(S313)。この処理は、上記実施の形態1で説明した図1のS106に部分的に対応する。
次に、複合部320は、S312の処理で選択した基準局毎に、S313の処理で生成した移動局用補正値を複合し複合移動局用補正値を算出する(S313)。この処理は、上記実施の形態1で説明した図1のS107に対応する。
そして、受信機310は測位衛星400が発信する測位信号を受信し、受信した測位信号に基づく観測値と、S314の処理で算出した複合移動局用補正値とに基づいて測位する(S315)。この処理は、上記実施の形態1で説明した図1のS108に対応する。
Next, the
Next, the
Next, based on the received correction value parameter, the
Next, the
Then, the
上記に説明した測位システムの各構成により移動局300の測位処理において精度の高い補正値を適用し、精度の高い測位結果を得ることができる。また、実施の形態3では移動局300と基準補正値生成局200との双方向通信により、基準補正値生成局200から移動局300に補正値の送信を行った。
With each configuration of the positioning system described above, a highly accurate correction value can be applied in the positioning process of the
上記実施の形態3において、移動局300は、基準補正値を補正値時刻変化率に基づいて移動局用補正値を生成し、生成した移動局用補正値を複合し測位する場合について説明した。
但し、移動局300は、基準補正値と補正値時刻変化率とそれぞれで複合した後に、複合した基準補正値を複合した補正値時刻変化率に基づいて複合移動局用補正値を生成し測位しても構わない。
In the third embodiment, the
However, the
実施の形態4.
実施の形態4では、上記実施の形態1で説明した測位方法を実行するシステムの一例について説明する。
Embodiment 4 FIG.
In the fourth embodiment, an example of a system that executes the positioning method described in the first embodiment will be described.
図10は、実施の形態4における測位システムの構成図である。
図10において、測位システムは、基準局100、基準補正値生成局200、移動局300を備え、測位衛星400を利用する。
基準局100は受信機110と通信機120とを備え、基準補正値生成局200は収集局210と生成局220と通信局230と複合局240とを備え、移動局300は受信機310と通信機330とを備える。
FIG. 10 is a configuration diagram of a positioning system in the fourth embodiment.
In FIG. 10, the positioning system includes a
The
図11は、実施の形態4における測位処理のフローチャートである。 FIG. 11 is a flowchart of the positioning process in the fourth embodiment.
まず、移動局300は、移動局300の概算座標を算出する(S401)。この処理は、上記実施の形態1で説明した図1のS101に対応する。
次に、通信機330は、移動局300がS401の処理で算出した概算座標を基準補正値生成局200に送信する(S402)。
First, the
Next, the
また、基準局100の受信機110は、測位衛星400が発信する測位信号を受信する(S403)。
The
そして、基準補正値生成局200の通信局230は、移動局300がS402の処理で送信した移動局300の概算座標を受信する(S404)。
また、収集局210は、各基準局100から観測値を収集する。この時、各基準局100では通信機120が受信機110の受信した測位信号に基づく観測値を基準補正値生成局200に送信する(S405)。この処理は、上記実施の形態1で説明した図1のS102に対応する。
次に、生成局220は、通信局230がS404の処理で受信した移動局300の概算座標と収集局210がS405の処理で収集した各基準局100の観測値とに基づいて補正値時刻変化率および補正値座標変化率を算出する(S406)。この処理は、上記実施の形態1で説明した図1のS103に対応する。
次に、生成局220は、収集局210がS405の処理で収集した各基準局100の観測値と既知である各基準局100の座標とに基づいて基準補正値を生成する(S407)。この処理は、上記実施の形態1で説明した図1のS104に対応する。
次に、複合局240は、通信局230がS404の処理で受信した移動局300の概算座標と既知である各基準局100の座標とに基づいて複合する補正値に対応する基準局100を選択する(S408)。この処理は、上記実施の形態1で説明した図1のS105に対応する。
次に、生成局220は、S407の処理で生成した基準補正値に、S406の処理で算出した補正値座標変化率に基づく補正値変化量を加算して、補正値パラメータの一つである、基準補正値を更新する(S409)。この処理は、上記実施の形態1で説明した図1のS106に部分的に対応する。
次に、複合局240は、複合局240がS409の処理で選択し生成局220がS409の処理で更新した各基準局100の補正値パラメータである基準補正値と、生成局220がS406の処理で算出した補正値パラメータである補正値時刻変化率をそれぞれ複合し補正値パラメータを算出する(S410)。この処理は、上記実施の形態1で説明した図1のS107に部分的に対応する。
次に、通信局230は、複合局240がS410の処理で算出した補正値パラメータを移動局300に送信する(S411)。
Then, the
The
Next, the
Next, the
Next, the
Next, the
Next, the
Next, the
次に、移動局300の通信機330は、基準補正値生成局200がS411の処理で送信した補正値パラメータを受信する(S412)。
次に受信機310は、通信機330がS412の処理で受信した補正値パラメータを基に移動局用補正値を生成する(S413)。この処理は、上記実施の形態1で説明した図1のS106に部分的に対応する。
そして、受信機310は測位衛星400が発信する測位信号を受信し、受信した測位信号に基づく観測値と、S413の処理で生成した移動局用補正値とに基づいて測位する(S414)。この処理は、上記実施の形態1で説明した図1のS108に対応する。
Next, the
Next, the
Then, the
上記に説明した測位システムの各構成により移動局300の測位処理において精度の高い補正値を適用し、精度の高い測位結果を得ることができる。
With each configuration of the positioning system described above, a highly accurate correction value can be applied in the positioning process of the
上記実施の形態2〜4において、上記実施の形態1で説明した測位方法を実行するシステムの一例について説明した。但し、測位システムの構成は、上記実施の形態2〜4で説明した構成に限らない。また、測位システムの処理の流れは、上記実施の形態2〜4で説明した処理の流れに限らない。
例えば、基準補正値生成局200は複数設置してもよい。その場合は、移動局300が適切な基準補正値生成局200を選択して、選択した基準補正値生成局200から補正値を受信する。
また例えば、複数の基準補正値生成局200と各基準補正値生成局200を管理する主補正値生成局を設置してもよい。その場合は、移動局300は主補正値生成局に問い合わせ、主補正値生成局は移動局300からの問い合わせに対して適切な基準補正値生成局200を選択し、移動局300は主補正値生成局が選択した基準補正値生成局200から補正値を受信する。または、移動局300が主補正値生成局に問い合わせると主補正値生成局は適切な基準補正値生成局を選択して選択した基準補正値生成局の補正値だけを移動局300に配信する。
また例えば、各基準局100は基準補正値を生成し、基準補正値生成局200は各基準局100が生成した補正値を収集し、収集した補正値を元に基準局の選択、補正値の複合、変化率に対応した補正値の生成をしてもよい。
In the said Embodiment 2-4, an example of the system which performs the positioning method demonstrated in the said
For example, a plurality of reference correction
Further, for example, a plurality of reference correction
Also, for example, each
上記の各実施の形態において、複数の基準局が生成する補正値や基準局および移動局の観測値を併用することで、移動局から基準局までの距離が要求より長い場合でも精度を確保する方法について説明した。これは、複数の基準局が生成する補正値や基準局および移動局の観測値の間に、空間的な相関がある事を利用した方法で、複数の基準局の補正値や観測値を適切に組み合わせることで、基準局から離れた場所でも要求精度を達成可能な補正値を得る方法である。
そこで、移動局が当該区域内を移動する際に、当該区域内の任意の地点で要求精度を達成できる補正値を得る複合補正値計算方式について説明した。
また、差動衛星測位方式の計算結果に不連続が発生しないように補正値を得るための複数基準局の選択方式について説明した。
In each of the above embodiments, a method for ensuring accuracy even when the distance from the mobile station to the reference station is longer than required by using correction values generated by a plurality of reference stations and observation values of the reference station and the mobile station together is described. did. This is a method that utilizes the fact that there is a spatial correlation between the correction values generated by multiple reference stations and the observation values of the reference station and mobile station, and appropriately combines the correction values and observation values of multiple reference stations. Thus, it is a method of obtaining a correction value that can achieve the required accuracy even at a location away from the reference station.
Thus, a description has been given of a composite correction value calculation method for obtaining a correction value that can achieve the required accuracy at an arbitrary point in the area when the mobile station moves within the area.
In addition, the selection method of a plurality of reference stations for obtaining correction values so that discontinuity does not occur in the calculation results of the differential satellite positioning method has been described.
図12は、各実施の形態における基準局100、基準補正値生成局200、移動局300のハードウェア構成図である。
図12において、基準局100、基準補正値生成局200、移動局300は、プログラムを実行するCPU(Central Processing Unit)911を備えている。CPU911は、バス912を介してROM913、RAM914、通信ボード915、磁気ディスク装置920と接続されている。
RAM914は、揮発性メモリの一例である。ROM913、磁気ディスク装置920は、不揮発性メモリの一例である。これらは、記憶装置あるいは記憶部の一例である。
通信ボード915はLAN、インターネット等に接続される。また、通信ボード915は情報入力部、出力部の一例である。
FIG. 12 is a hardware configuration diagram of the
In FIG. 12, a
The
The
磁気ディスク装置920には、オペレーティングシステム(OS)921、プログラム群923、ファイル群924が記憶されている。プログラム群923は、CPU911、OS921により実行される。
The
上記プログラム群923には、各実施の形態の説明において「〜部」、「〜局」として説明する機能を実行するプログラムが記憶されている。プログラムは、CPU911により読み出され実行される。
ファイル群924には、各実施の形態の説明において、「〜を判定し」、「〜を判定した結果」、「〜を計算し」、「〜を計算した結果」、「〜を処理し」、「〜を処理した結果」のような表現で説明する結果情報が、「〜ファイル」として記憶されている。
また、各実施の形態の説明において説明するフローチャートの矢印の部分は主としてデータの入出力を示し、そのデータの入出力のためにデータは、磁気ディスク装置920、FD(Flexible Disk cartridge)、光ディスク、CD(コンパクトディスク)、MD(ミニディスク)、DVD(Digital Versatile Disk)、その他の記録媒体に記録される。あるいは、信号線やその他の伝送媒体により伝送される。
The
In the description of each embodiment, the
In addition, the arrows in the flowcharts described in the description of each embodiment mainly indicate input / output of data. For the input / output of the data, the data is a
また、各実施の形態の説明において「〜部」、「〜局」として説明するものは、ROM913に記憶されたファームウェアで実現されていても構わない。或いは、ソフトウェアのみ、或いは、ハードウェアのみ、或いは、ソフトウェアとハードウェアとの組み合わせ、さらには、ファームウェアとの組み合わせで実施されても構わない。
In addition, what is described as “˜unit” and “˜station” in the description of each embodiment may be realized by firmware stored in the
また、各実施の形態を実施するプログラムは、磁気ディスク装置920、FD、光ディスク、CD、MD、DVD、その他の記録媒体による記録装置を用いて記憶されても構わない。
In addition, the program for implementing each embodiment may be stored using a recording device using a
100 基準局、110 受信機、120 通信機、200 基準補正値生成局、210 収集局、220 生成局、230 通信局、240 複合局、300 移動局、310 受信機、320 複合部、330 通信機、400 測位衛星、911 CPU、912 バス、913 ROM、914 RAM、915 通信ボード、920 磁気ディスク装置、921 OS、923 プログラム群、924 ファイル群。 100 reference station, 110 receiver, 120 communication device, 200 reference correction value generation station, 210 collection station, 220 generation station, 230 communication station, 240 composite station, 300 mobile station, 310 receiver, 320 composite unit, 330 communication device, 400 positioning satellite, 911 CPU, 912 bus, 913 ROM, 914 RAM, 915 communication board, 920 magnetic disk unit, 921 OS, 923 program group, 924 file group.
Claims (12)
前記観測値受信部が入力した観測値に基づいて各測位衛星毎に基準補正値を生成する基準補正値生成部と、
単位時間あたりの基準補正値の変化量を示す補正値時刻変化率と、測位対象である移動局と前記基準局との座標差の単位座標あたりの基準補正値の変化量を示す補正値座標変化率との少なくともいずれかの変化率を、前記基準局から各測位衛星に向けた各仰角に基づいて各測位衛星毎に算出する変化率算出部と、
前記移動局の座標を測位する際の補正処理に使用する補正値パラメータとして、前記基準補正値生成部が生成した各測位衛星毎の基準補正値と前記変化率算出部が算出した各測位衛星毎の変化率とを前記移動局に送信する補正値パラメータ送信部と
を備えたことを特徴とする補正値パラメータ生成装置。 An observation value receiving unit for receiving an observation value obtained by observing a positioning signal from each positioning satellite from a reference station;
A reference correction value generating unit that generates a reference correction value for each positioning satellite based on the observation value input by the observation value receiving unit;
Correction value time change rate indicating the change amount of the reference correction value per unit time, and correction value coordinate change rate indicating the change amount of the reference correction value per unit coordinate of the coordinate difference between the mobile station to be positioned and the reference station A rate-of-change calculation unit that calculates the rate of change of at least one of each positioning satellite based on each elevation angle from the reference station toward each positioning satellite;
As a correction value parameter used for correction processing when positioning the coordinates of the mobile station, a reference correction value for each positioning satellite generated by the reference correction value generation unit and each positioning satellite calculated by the change rate calculation unit And a correction value parameter transmitting unit that transmits the rate of change to the mobile station.
前記観測値受信部が入力した観測値に基づいて各基準局と各測位衛星と毎に基準補正値を生成する基準補正値生成部と、
単位時間あたりの基準補正値の変化量を示す補正値時刻変化率を、各基準局から各測位衛星に向けた各仰角に基づいて各基準局と各測位衛星と毎に算出する変化率算出部と、
前記移動局と各基準局との各距離と各基準局から各測位衛星に向けた各仰角との少なくともいずれかに基づいて、前記基準補正値生成部が生成した各基準補正値と前記変化率算出部が算出した各補正値時刻変化率とをそれぞれで重み付けし、重み付けした各基準補正値を各測位衛星毎に合算した値である複合基準補正値と重み付けした各補正値時刻変化率を各測位衛星毎に合算した値である複合補正値時刻変化率とを算出する基準補正値複合部と、
前記移動局の座標を測位する際の補正処理に使用する補正値パラメータとして、前記基準補正値複合部が算出した前記複合基準補正値と前記複合補正値時刻変化率とを前記移動局に送信する補正値パラメータ送信部と
を備えたことを特徴とする補正値パラメータ生成装置。 An observation value receiving unit for receiving an observation value obtained by observing a positioning signal from each positioning satellite from each reference station;
A reference correction value generating unit that generates a reference correction value for each reference station and each positioning satellite based on the observation value input by the observation value receiving unit;
A change rate calculation unit that calculates a correction value time change rate indicating a change amount of the reference correction value per unit time for each reference station and each positioning satellite based on each elevation angle from each reference station toward each positioning satellite;
Each reference correction value generated by the reference correction value generation unit and the change rate calculation unit based on at least one of each distance between the mobile station and each reference station and each elevation angle from each reference station toward each positioning satellite Each correction value time change rate calculated by is weighted individually, and each weighted reference correction value is added to each positioning satellite, and the combined reference correction value and each weighted correction value time change rate are assigned to each positioning satellite. A reference correction value composite unit for calculating a composite correction value time change rate that is a value summed up every time;
The composite reference correction value calculated by the reference correction value combination unit and the composite correction value time change rate are transmitted to the mobile station as correction value parameters used for correction processing when positioning the coordinates of the mobile station. A correction value parameter generation device comprising a correction value parameter transmission unit.
前記基準補正値複合部が重み付けし合算する基準補正値と補正値時刻変化率とを前記移動局と各基準局との距離に基づいて選択する生成装置選択部
を備えたことを特徴とする請求項2記載の補正値パラメータ生成装置。 The correction value parameter generation device further includes:
The generator selection unit for selecting a reference correction value weighted and summed by the reference correction value composite unit and a correction value time change rate based on a distance between the mobile station and each reference station. 2. The correction value parameter generation device according to 2.
各測位衛星毎の特定時刻の補正値である基準補正値と、単位時間あたりの基準補正値の変化量を示す補正値時刻変化率であり各測位衛星毎の補正値時刻変化率とを受信する測位装置受信部と、
前記測位装置受信部が受信した前記補正値時刻変化率に前記特定時刻からの経過時間を乗じた値を各測位衛星毎に算出し、算出した値を前記測位装置受信部が受信した前記基準補正値に加算した値である移動局用補正値を各測位衛星毎に生成する移動局用補正値生成部と、
前記移動局用補正値生成部が生成した各測位衛星毎の移動局用補正値を用いて差動衛星測位方式により自己の位置する座標を測位する測位部と
を備えたことを特徴とする測位装置。 In the positioning device of the mobile station that measures the coordinates where it is located,
A reference correction value that is a correction value at a specific time for each positioning satellite and a correction value time change rate that indicates the amount of change in the reference correction value per unit time and a correction value time change rate for each positioning satellite are received. A positioning device receiver;
A value obtained by multiplying the correction value time change rate received by the positioning device receiving unit by an elapsed time from the specific time is calculated for each positioning satellite, and the reference correction received by the positioning device receiving unit is calculated. A mobile station correction value generating unit that generates a mobile station correction value that is a value added to the value for each positioning satellite;
A positioning unit comprising: a positioning unit that measures the coordinates of the mobile station using a differential satellite positioning method using the mobile station correction value for each positioning satellite generated by the mobile station correction value generation unit. apparatus.
各基準局と各測位衛星と毎の特定時刻の補正値である基準補正値と、単位時間あたりの基準補正値の変化量を示す補正値時刻変化率であり各基準局と各測位衛星と毎の補正値時刻変化率とを受信する測位装置受信部と、
前記測位装置受信部が受信した前記補正値時刻変化率に前記特定時刻からの経過時間を乗じた値を各基準局と各測位衛星と毎に算出し、算出した値を前記測位装置受信部が受信した前記基準補正値に加算した値である移動局用補正値を各基準局と各測位衛星と毎に生成する移動局用補正値生成部と、
前記移動局と各基準局との各距離と各基準局から各測位衛星に向けた各仰角との少なくともいずれかに基づいて、前記移動局用補正値生成部が生成した各移動局用補正値を重み付けし、重み付けした各移動局用補正値を各測位衛星毎に合算した値である複合移動局用補正値を算出する測位装置複合部と、
前記測位装置複合部が算出した各測位衛星毎の複合移動局用補正値を用いて差動衛星測位方式により自己の位置する座標を測位する測位部と
を備えたことを特徴とする測位装置。 In the positioning device of the mobile station that measures the coordinates where it is located,
A reference correction value that is a correction value at a specific time for each reference station and each positioning satellite, and a correction value time change rate that indicates the amount of change of the reference correction value per unit time, and a correction for each reference station and each positioning satellite A positioning device receiving unit for receiving the value time rate of change;
A value obtained by multiplying the correction value time change rate received by the positioning device reception unit by an elapsed time from the specific time is calculated for each reference station and each positioning satellite, and the calculated value is received by the positioning device reception unit. A mobile station correction value generating unit that generates a mobile station correction value that is a value added to the reference correction value for each reference station and each positioning satellite;
Each mobile station correction value generated by the mobile station correction value generation unit is weighted based on at least one of the distance between the mobile station and each reference station and each elevation angle from each reference station toward each positioning satellite. A positioning device composite unit that calculates a composite mobile station correction value that is a sum of the weighted correction values for each mobile station for each positioning satellite;
A positioning device comprising: a positioning unit that measures the coordinates of its own position by a differential satellite positioning method using the composite mobile station correction value for each positioning satellite calculated by the positioning device composite unit.
各基準局と各測位衛星と毎の特定時刻の補正値である基準補正値と、単位時間あたりの基準補正値の変化量を示す補正値時刻変化率であり各基準局と各測位衛星と毎の補正値時刻変化率と共に、前記移動局と各基準局との座標差の単位座標あたりの基準補正値の変化量を示す補正値座標変化率であり各基準局と各測位衛星と毎の補正値座標変化率を受信し、
前記移動局用補正値生成部は、
前記測位装置受信部が受信した前記補正値時刻変化率に前記特定時刻からの経過時間を乗じた値を各基準局と各測位衛星と毎に算出し、前記測位装置受信部が受信した前記基準補正値に算出した値を加算すると共に、前記測位装置受信部が受信した前記補正値座標変化率に前記移動局と各基準局との座標差を乗じた値を加算して移動局用補正値を各基準局と各測位衛星と毎に生成する
ことを特徴とする請求項5記載の測位装置。 The positioning device receiver is
A reference correction value that is a correction value at a specific time for each reference station and each positioning satellite, and a correction value time change rate that indicates the amount of change of the reference correction value per unit time, and a correction for each reference station and each positioning satellite A correction value coordinate change rate indicating a change amount of a reference correction value per unit coordinate of a coordinate difference between the mobile station and each reference station together with a value time change rate, and a correction value coordinate change rate for each reference station and each positioning satellite Receive
The mobile station correction value generator is
A value obtained by multiplying the correction value time change rate received by the positioning device reception unit by an elapsed time from the specific time is calculated for each reference station and each positioning satellite, and the reference correction received by the positioning device reception unit In addition to adding the calculated value to the value, and adding the value obtained by multiplying the correction value coordinate change rate received by the positioning device receiver by the coordinate difference between the mobile station and each reference station, each mobile station correction value is 6. The positioning device according to claim 5, wherein the positioning device is generated for each reference station and each positioning satellite.
前記測位装置複合部が重み付けし合算する移動局用補正値を前記移動局と各基準局との距離に基づいて選択する測位装置選択部
を備えたことを特徴とする請求項5または請求項6いずれかに記載の測位装置。 The positioning device further includes:
7. A positioning device selection unit that selects a mobile station correction value weighted and summed by the positioning device composite unit based on a distance between the mobile station and each reference station. A positioning device according to the above.
各基準局と各測位衛星と毎の特定時刻の補正値である基準補正値と、単位時間あたりの基準補正値の変化量を示す補正値時刻変化率であり各基準局と各測位衛星と毎の補正値時刻変化率とを受信する測位装置受信部と、
前記移動局と各基準局との各距離と各基準局から各測位衛星に向けた各仰角との少なくともいずれかに基づいて、前記測位装置受信部が受信した各基準補正値と各補正値時刻変化率とをそれぞれで重み付けし、重み付けした各基準補正値を各測位衛星毎に合算した値である複合基準補正値と重み付けした各補正値時刻変化率を各測位衛星毎に合算した値である複合補正値時刻変化率とを算出する測位装置複合部と、
前記測位装置複合部が算出した各複合補正値時刻変化率に前記特定時刻からの経過時間を乗じた値を算出し、算出した値を前記測位装置複合部が算出した前記複合基準補正値に加算した値である複合移動局用補正値を各測位衛星毎に生成する移動局用補正値生成部と、
前記移動局用補正値生成部が生成した各測位衛星毎の複合移動局用補正値を用いて差動衛星測位方式により自己の位置する座標を測位する測位部と
を備えたことを特徴とする測位装置。 In the positioning device of the mobile station that measures the coordinates where it is located,
A reference correction value that is a correction value at a specific time for each reference station and each positioning satellite, and a correction value time change rate that indicates the amount of change of the reference correction value per unit time, and a correction for each reference station and each positioning satellite A positioning device receiving unit for receiving the value time rate of change;
Each reference correction value and each correction value time change rate received by the positioning device receiver based on at least one of the distance between the mobile station and each reference station and each elevation angle from each reference station toward each positioning satellite And the weighted reference correction value for each positioning satellite, and the combined reference correction value for each positioning satellite, and the weighted correction value for each positioning satellite. A positioning device composite unit for calculating a value time change rate;
A value obtained by multiplying each composite correction value time change rate calculated by the positioning device composite unit by an elapsed time from the specific time is calculated, and the calculated value is added to the composite reference correction value calculated by the positioning device composite unit. A mobile station correction value generating unit that generates a composite mobile station correction value for each positioning satellite,
And a positioning unit for positioning the coordinates of the mobile station by using a differential satellite positioning method using the composite mobile station correction value for each positioning satellite generated by the mobile station correction value generating unit. Positioning device.
前記測位装置複合部が重み付けし合算する基準補正値と補正値時刻変化率とを前記移動局と各基準局との距離に基づいて選択する測位装置選択部
を備えたことを特徴とする請求項8記載の測位装置。 The positioning device further includes:
9. A positioning device selection unit that selects a reference correction value weighted and summed by the positioning device composite unit and a correction value time change rate based on a distance between the mobile station and each reference station. The described positioning device.
単位時間あたりの補正値の変化量を示す補正値時刻変化率と、前記移動局と前記基準局との座標差の単位座標あたりの補正値の変化量を示す補正値座標変化率との少なくともいずれかの変化率を、前記基準局から各測位衛星に向けた各仰角に基づいて算出し、
算出した前記変化率に基づく変化量を前記補正値に加算して補正値を更新し、
更新した補正値を用いて移動局の座標を測位する
ことを特徴とする差動衛星測位方法。 In the differential satellite positioning method of positioning the coordinates of the mobile station using the correction value based on the observation value for each positioning satellite of the reference station,
At least one of a correction value time change rate indicating the change amount of the correction value per unit time and a correction value coordinate change rate indicating the change amount of the correction value per unit coordinate of the coordinate difference between the mobile station and the reference station. Is calculated based on each elevation angle from the reference station to each positioning satellite,
The amount of change based on the calculated rate of change is added to the correction value to update the correction value,
A differential satellite positioning method characterized by positioning the coordinates of a mobile station using an updated correction value.
前記移動局と各基準局との各距離と各基準局から各測位衛星に向けた各仰角との少なくともいずれかに基づいて各基準局毎の補正値を重み付けし、重み付けした各基準局毎の補正値を合算した値である複合補正値を算出し、
算出した前記複合補正値を用いて移動局の座標を測位する
ことを特徴とする差動衛星測位方法。 In the differential satellite positioning method for positioning the coordinates of the mobile station based on the correction value for each reference station based on the observation value for each positioning satellite of each reference station,
The correction value for each reference station is weighted based on at least one of the distance between the mobile station and each reference station and each elevation angle from each reference station toward each positioning satellite, and the weighted correction values for each reference station are added together. Calculate the composite correction value that is
A differential satellite positioning method characterized by positioning the coordinates of a mobile station using the calculated composite correction value.
前記移動局と各基準局との距離に基づいて各基準局毎の補正値から複数の各基準局毎の補正値を選択し、
前記移動局と各基準局との各距離と各基準局から各測位衛星に向けた各仰角との少なくともいずれかに基づいて、選択した各基準局毎の補正値を重み付けし、重み付けした各基準局毎の補正値を合算した値である複合補正値を算出し、
算出した前記複合補正値を用いて移動局の座標を測位する
ことを特徴とする差動衛星測位方法。 In the differential satellite positioning method for positioning the coordinates of the mobile station based on the correction value for each reference station based on the observation value for each positioning satellite of each reference station,
Based on the distance between the mobile station and each reference station, select a correction value for each reference station from a correction value for each reference station,
Based on at least one of each distance between the mobile station and each reference station and each elevation angle from each reference station toward each positioning satellite, the correction value for each selected reference station is weighted, and the weighted correction for each reference station Calculate the composite correction value that is the sum of the values,
A differential satellite positioning method characterized by positioning the coordinates of a mobile station using the calculated composite correction value.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005039237A JP2006226772A (en) | 2005-02-16 | 2005-02-16 | Correction value parameter generating device, positioning device, and differential satellite positioning method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005039237A JP2006226772A (en) | 2005-02-16 | 2005-02-16 | Correction value parameter generating device, positioning device, and differential satellite positioning method |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2006226772A true JP2006226772A (en) | 2006-08-31 |
Family
ID=36988285
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2005039237A Withdrawn JP2006226772A (en) | 2005-02-16 | 2005-02-16 | Correction value parameter generating device, positioning device, and differential satellite positioning method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2006226772A (en) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7663534B2 (en) | 2004-06-22 | 2010-02-16 | Robert Bosch Gmbh | Radar sensor and method for analyzing objects |
KR101344426B1 (en) | 2012-03-23 | 2013-12-24 | 주식회사 지평스페이스 | An apparatus for processing differential information of differential global navigation satellite system and the method thereof |
JP2014016315A (en) * | 2012-07-11 | 2014-01-30 | Mitsubishi Electric Corp | Transmitter and positioning device |
KR101419339B1 (en) * | 2013-10-07 | 2014-07-16 | 주식회사 지평스페이스 | An apparatus for processing differential information of differential global navigation satellite system and the method thereof |
CN110426723A (en) * | 2019-07-25 | 2019-11-08 | 武汉星源云意科技有限公司 | A kind of method of the acquisition and map publishing of satellite positioning GGA data |
KR20190139616A (en) * | 2018-06-08 | 2019-12-18 | 현대자동차주식회사 | Satellite navigation system and signal processing method thereof |
KR20200063704A (en) * | 2018-11-28 | 2020-06-05 | 한양대학교 산학협력단 | Terminal, base station and location positioning method |
JP2021043175A (en) * | 2019-09-11 | 2021-03-18 | コリア エクスプレスウェイ コーポレーション | Device and method for precise position correction using positioning difference |
CN114185069A (en) * | 2021-12-07 | 2022-03-15 | 西藏金采科技股份有限公司 | Differential positioning method, device and system for improving Beidou positioning accuracy |
CN116840880A (en) * | 2023-06-07 | 2023-10-03 | 武汉大学 | Method for positioning real-time relative position between dynamic targets and computer readable medium |
CN117055087A (en) * | 2023-10-10 | 2023-11-14 | 中国电建集团西北勘测设计研究院有限公司 | GNSS real-time positioning and resolving method for haze influence area |
-
2005
- 2005-02-16 JP JP2005039237A patent/JP2006226772A/en not_active Withdrawn
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7663534B2 (en) | 2004-06-22 | 2010-02-16 | Robert Bosch Gmbh | Radar sensor and method for analyzing objects |
KR101344426B1 (en) | 2012-03-23 | 2013-12-24 | 주식회사 지평스페이스 | An apparatus for processing differential information of differential global navigation satellite system and the method thereof |
JP2014016315A (en) * | 2012-07-11 | 2014-01-30 | Mitsubishi Electric Corp | Transmitter and positioning device |
KR101419339B1 (en) * | 2013-10-07 | 2014-07-16 | 주식회사 지평스페이스 | An apparatus for processing differential information of differential global navigation satellite system and the method thereof |
KR102575716B1 (en) * | 2018-06-08 | 2023-09-07 | 현대자동차주식회사 | Satellite navigation system and signal processing method thereof |
KR20190139616A (en) * | 2018-06-08 | 2019-12-18 | 현대자동차주식회사 | Satellite navigation system and signal processing method thereof |
KR20200063704A (en) * | 2018-11-28 | 2020-06-05 | 한양대학교 산학협력단 | Terminal, base station and location positioning method |
KR102188880B1 (en) * | 2018-11-28 | 2020-12-09 | 한양대학교 산학협력단 | Terminal, base station and location positioning method |
CN110426723A (en) * | 2019-07-25 | 2019-11-08 | 武汉星源云意科技有限公司 | A kind of method of the acquisition and map publishing of satellite positioning GGA data |
CN110426723B (en) * | 2019-07-25 | 2022-12-23 | 武汉星源云意科技有限公司 | Method for acquiring satellite positioning GGA data and publishing map |
US11175408B2 (en) | 2019-09-11 | 2021-11-16 | Korea Expressway Corp. | Apparatus and method for precise position correction using positioning difference |
JP2021043175A (en) * | 2019-09-11 | 2021-03-18 | コリア エクスプレスウェイ コーポレーション | Device and method for precise position correction using positioning difference |
CN114185069A (en) * | 2021-12-07 | 2022-03-15 | 西藏金采科技股份有限公司 | Differential positioning method, device and system for improving Beidou positioning accuracy |
CN116840880A (en) * | 2023-06-07 | 2023-10-03 | 武汉大学 | Method for positioning real-time relative position between dynamic targets and computer readable medium |
CN116840880B (en) * | 2023-06-07 | 2024-02-13 | 武汉大学 | Method for positioning real-time relative position between dynamic targets and computer readable medium |
CN117055087A (en) * | 2023-10-10 | 2023-11-14 | 中国电建集团西北勘测设计研究院有限公司 | GNSS real-time positioning and resolving method for haze influence area |
CN117055087B (en) * | 2023-10-10 | 2024-01-30 | 中国电建集团西北勘测设计研究院有限公司 | GNSS real-time positioning and resolving method for haze influence area |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2006226772A (en) | Correction value parameter generating device, positioning device, and differential satellite positioning method | |
Paziewski et al. | Assessment of GPS+ Galileo and multi-frequency Galileo single-epoch precise positioning with network corrections | |
Wabbena et al. | PPP-RTK: precise point positioning using state-space representation in RTK networks | |
JP4644792B2 (en) | Method and apparatus for correcting ionospheric delay in satellite navigation system. | |
WO2017046914A1 (en) | Positioning satellite selecting device, positioning device, positioning system, positioning information transmitting device and positioning terminal | |
US10393879B2 (en) | Global positioning device | |
WO2018083803A1 (en) | Positioning augmentation device, positioning augmentation system, and positioning augmentation method | |
US20230251388A1 (en) | Navigation satellite system reception device, method for processing navigation satellite signal from same, and program | |
US20090238042A1 (en) | Methods and systems for underwater navigation | |
AU2014398394B2 (en) | Positioning device, positioning method, and program | |
US11460583B2 (en) | Method and apparatus for providing correction data for satellite navigation | |
JP2007171082A (en) | Ionosphere delay correction method, ionosphere delay correction system, earth station, and mobile station | |
US11061143B2 (en) | Global navigation satellite system, navigation terminal, navigation method and program | |
US20150247931A1 (en) | Locally enhanced gnss wide-area augmentation system | |
JP2006258461A (en) | Positioning device, positioning method and positioning program | |
JP2007248177A (en) | Information processor, and correction value compounding method and program | |
JP2016128837A (en) | Gnss positioning device and gnss positioning method | |
JP2011127939A (en) | Moving body position estimating/tracking device, method of estimating/tracking position of moving body, and moving body position estimating/tracking program | |
CN111381253A (en) | System for estimating position of global navigation satellite system receiver and method thereof | |
JP2013181872A (en) | Relative positioning device, relative positioning method, and program | |
JP4413093B2 (en) | Satellite signal receiver | |
KR20170094988A (en) | DGNSS(Differential Global Navigation Satellite System)-BASED POSITIONING METHOD AND GNSS RECEIVER USING THE SAME | |
JP2006162563A (en) | Data filter, data filtering method, data filtering program, positioning device, positioning method, and positioning program | |
JP7081716B2 (en) | Detection method and detection program | |
US20140011516A1 (en) | Method and apparatus for estimating location of terminal using generation of virtual infrastructures |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Effective date: 20071212 Free format text: JAPANESE INTERMEDIATE CODE: A621 |
|
A761 | Written withdrawal of application |
Effective date: 20090616 Free format text: JAPANESE INTERMEDIATE CODE: A761 |