Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7570623B2 - Satellite-based positioning method and positioning device - Google Patents
[go: Go Back, main page]

JP7570623B2 - Satellite-based positioning method and positioning device - Google Patents

Satellite-based positioning method and positioning device Download PDF

Info

Publication number
JP7570623B2
JP7570623B2 JP2020068055A JP2020068055A JP7570623B2 JP 7570623 B2 JP7570623 B2 JP 7570623B2 JP 2020068055 A JP2020068055 A JP 2020068055A JP 2020068055 A JP2020068055 A JP 2020068055A JP 7570623 B2 JP7570623 B2 JP 7570623B2
Authority
JP
Japan
Prior art keywords
satellite
satellites
positioning
periodic
difference
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
JP2020068055A
Other languages
Japanese (ja)
Other versions
JP2021165640A (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.)
Kinkei System Corp
Original Assignee
Kinkei System 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 Kinkei System Corp filed Critical Kinkei System Corp
Priority to JP2020068055A priority Critical patent/JP7570623B2/en
Publication of JP2021165640A publication Critical patent/JP2021165640A/en
Application granted granted Critical
Publication of JP7570623B2 publication Critical patent/JP7570623B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)

Description

本発明は、衛星を用いた測位方法及び測位装置に関する。 The present invention relates to a positioning method and device using satellites.

斜面崩壊の監視においては、前兆となる地表面の動きが小さいため数mmの変位を捉える必要がある。地表面の変位を捉えるために、GNSS(Global Navigation Satellite System)を用いた測位技術を利用することが知られている。測位技術の中で、キネマティック法の一種であるRTK(Real Time Kinematic)法は、リアルタイムに高精度測位が可能である。特許文献1には、RTK法を用いてmmレベルの測定精度を実現する測位方法が開示されている。特許文献1の方法では、測定時と、測定時から恒星日の整数倍異なる時点である基準時において、双方の時点に共通する衛星の組合せを用いてそれぞれ算出された位置情報の差から変位を測定する。 In monitoring slope collapse, since the movements of the ground surface that are precursors are small, it is necessary to capture displacements of a few millimeters. It is known to use positioning technology using the Global Navigation Satellite System (GNSS) to capture displacements of the ground surface. Among positioning technologies, the Real Time Kinematic (RTK) method, which is a type of kinematic method, is capable of high-precision positioning in real time. Patent Document 1 discloses a positioning method that achieves measurement accuracy on the millimeter level using the RTK method. In the method of Patent Document 1, displacement is measured from the difference in position information calculated at the time of measurement and at a reference time, which is a time that differs from the measurement time by an integer multiple of a sidereal day, using a combination of satellites that are common to both times.

RTK法は、測定点と基準点でそれぞれ受信する電波の行路差に対応する整数値バイアスを確定する(Fix解を得る)初期化演算ステップと、その後に測定点の位置を計測する測位演算ステップの2段階の演算工程から構成される。初期化時間は、空が開けている場所では10~30秒であるが、周囲に構造部等の障害物が存在すると、反射波に起因するマルチパスの影響によってミスFixや初期化に数分を要する等の問題を生じる。非特許文献1では、信号の質の悪い衛星を排除することで測位結果を改善している。具体的には、マルチパスノイズが大きい低い仰角の衛星を排除する仰角マスク、受信感度の低い衛星を排除するCN(搬送波電力対雑音電力密度比)マスク、衛星システムの選択、測位開始タイミングの選択等である。測位開始タイミングの選択は、特定の衛星で起きるサイクルスリップの見逃しを改善する。 The RTK method is composed of two calculation steps: an initialization calculation step in which an integer bias corresponding to the path difference of the radio waves received at the measurement point and the reference point is determined (a fix solution is obtained), and then a positioning calculation step in which the position of the measurement point is measured. The initialization time is 10 to 30 seconds in an open sky area, but if there are obstacles such as structural parts in the vicinity, problems such as miss fixes and initialization taking several minutes occur due to the influence of multipath caused by reflected waves. In Non-Patent Document 1, the positioning results are improved by eliminating satellites with poor signal quality. Specifically, an elevation angle mask that eliminates satellites with low elevation angles with large multipath noise, a CN 0 (carrier power to noise power density ratio) mask that eliminates satellites with low reception sensitivity, satellite system selection, and positioning start timing selection. The selection of the positioning start timing improves the overlooking of cycle slips that occur in specific satellites.

特開平11-352210号公報Japanese Patent Application Publication No. 11-352210

「測位計算ソフトウェアの設定変更による測位性改善」土木学会論文集F3(土木情報学),Vol. 72, No. 2, II_47-II_54,2016"Improving positioning performance by changing the settings of positioning calculation software" Journal of the Japan Society of Civil Engineers, Vol. 72, No. 2, II_47-II_54, 2016

斜面崩壊の監視は、24時間連続監視が必要である上、推定される斜面領域に十分な数の測定点を設置する必要があり、また、斜面崩壊の監視場所は、通常、電源確保が困難である。したがって、バッテリや太陽電池等を用いた長期計測を行うことになり、電力消費を抑制できるシステムが好ましい。さらに、斜面崩壊の監視場所は、空が開けていない山間部が多いため、RTK法の初期化に時間を要することになるが、電力消費の観点から初期化を早く終えることが望ましい。 Monitoring for slope collapse requires continuous 24-hour monitoring, and a sufficient number of measurement points must be installed in the estimated slope area. Furthermore, it is usually difficult to secure a power source at the slope collapse monitoring locations. Therefore, long-term measurements are performed using batteries or solar cells, and a system that can reduce power consumption is preferable. Furthermore, since the slope collapse monitoring locations are often mountainous areas with no open sky, initialization of the RTK method takes time, but from the standpoint of power consumption it is desirable to complete the initialization quickly.

上述した衛星の仰角マスク、CNマスク、衛星システムの選択等は、いずれもノイズの大きい特定の衛星を排除することによって初期化時間を短縮している。しかしながら、いずれの衛星排除方法も、初期化演算又は測位演算を実行する前に、当該方法を適用すべきかが否かを判定できないという問題がある。従来は、初期化演算又は測位演算の時間短縮のため、1又は複数の衛星排除方法が適用される。 The above-mentioned satellite elevation mask, CN 0 mask, satellite system selection, etc., all reduce the initialization time by excluding specific satellites with large noise. However, each satellite exclusion method has a problem in that it is not possible to determine whether or not the method should be applied before performing the initialization calculation or positioning calculation. Conventionally, one or more satellite exclusion methods are applied to reduce the time of the initialization calculation or positioning calculation.

上記の問題点に鑑み、本発明は、周期性衛星からの電波を用いてRTK法による演算により測定点の測位を行う測位方法及び測位装置において、演算に利用する衛星の中から排除すべき衛星を効率的に選択でき、演算時間を短縮しかつ精度のよい測位を行えることを目的とする。 In view of the above problems, the present invention aims to provide a positioning method and device that uses radio waves from periodic satellites to perform calculations based on the RTK method to determine the position of a measurement point, and that can efficiently select satellites to be excluded from those used in the calculations, thereby shortening the calculation time and performing highly accurate positioning.

上記の問題を解消するべく、本発明は、以下の構成を提供する。
・ 本発明による衛星を用いた測位方法の態様は、測定点及び基準点において複数の衛星からそれぞれ受信する電波を用いて初期化演算ステップと測位演算ステップとを有するRTK法により測定点の測位を行う測位方法において、
前記初期化演算ステップの前に、又は、前記初期化演算ステップと前記測位演算ステップとの間に、その後の演算に利用する衛星を選択するために実行される衛星選択ステップを有し、前記衛星選択ステップが、
恒星日毎に周期性のある周期性衛星以外の衛星を、その後の演算に利用しない衛星として排除するステップと、
各前記周期性衛星について、測定点における電波の搬送波電力対雑音電力密度比であるCNを検出し記憶するステップと、
各前記周期性衛星について、検出されたCNを、1又は複数の恒星日前の記憶されたCNと比較し、所定値以上の差があるか否かを判定するステップと、
CNに所定値以上の差があると判定された周期性衛星を、その後の演算に利用しない衛星として排除するステップと、を有することを特徴とする。
・ 本発明による衛星を用いた測位方法の別の態様は、測定点及び基準点において複数の衛星からそれぞれ受信する電波を用いて初期化演算ステップと測位演算ステップとを有するRTK法により測定点の測位を行う測位方法において、
前記初期化演算ステップの前に、又は、前記初期化演算ステップと前記測位演算ステップとの間に、その後の演算に利用する衛星を選択するために実行される衛星選択ステップを有し、前記衛星選択ステップが、
恒星日毎に周期性のある周期性衛星以外の衛星を、その後の演算に利用しない衛星として排除するステップと、
各前記周期性衛星について、測定点における電波の搬送波電力対雑音電力密度比であるCNを検出し記憶するステップと、
各前記周期性衛星について、検出されたCNを、1又は複数の恒星日前の記憶されたCNと比較し、所定値以上の差があるか否かを判定するステップと、
CNに所定値以上の差があると判定された周期性衛星について、基準点における同一タイミングのCNと、測定点と同じ1又は複数の恒星日前の基準点におけるCNとの差を取得するステップと、
CNに所定値以上の差があると判定された周期性衛星について、当該CNの差が基準点におけるCNの差と同じであるか否かを判定するステップと、
CNに所定値以上の差があると判定された周期性衛星のうち、当該CNの差が基準点におけるCNの差と同じではないと判定された周期性衛星を、その後の演算に利用しない衛星として排除するステップと、を有することを特徴とする。
・ 上記いずれかの方法の態様において、前記衛星選択ステップが、
測定点の周囲環境情報に基づいて、測定点における排除対象の方位範囲を決定するステップと、
決定された前記方位範囲に該当する衛星を、その後の演算に利用しない衛星として排除するステップと、をさらに有することが、好適である。
・ 上記いずれかの方法の態様において、前記初期化演算ステップと前記測位演算ステップとを、一定時間の間隔で行うことが、好適である。
・ 上記いずれかの方法の態様において、前記測位演算ステップが、
利用する周期性衛星のうち5個以上の周期性衛星を含む全ての組合せについて、所定の時間範囲内に所定の時間間隔毎の電波を用いて測位演算を行うステップと、
各前記組合せにおける前記時間間隔毎の測位演算結果について1恒星日前との差分である変位量を計算し、前記時間範囲内の変位量の平均を計算するステップと、
変位量が閾値以上の組合せが過半数以上であるか否かを判定するステップと、
変位量が閾値以上の組合せが過半数以上でない場合、変位なしと判定するステップと、を有することが、好適である。
・ 上記いずれかの方法の態様において、変位量が閾値以上の組合せが過半数以上である場合、特定の周期性衛星を含む組合せのみが閾値を超えているか否かをさらに判定するステップと、
特定の周期性衛星を含む組合せのみが閾値を超えているのではない場合、変位ありと判定するステップと、を有することが、好適である。
・ 本発明による衛星を用いた測位装置の態様は、測定点及び基準点において複数の衛星からそれぞれ受信する電波を用いて、RTK法により測定点の測位を行うべく初期化演算を行う初期化演算部と測位演算を行う測位演算部とを有する測位装置において、
初期化演算の前に、又は、初期化演算と測位演算との間に、その後の演算に利用する衛星の選択を行う衛星選択部を有し、前記衛星選択部が、
恒星日毎に周期性のある周期性衛星以外の衛星を、その後の演算に利用しない衛星として排除する手段と、
各前記周期性衛星について、測定点における電波の搬送波電力対雑音電力密度比であるCNを検出して記憶する手段と、
各前記周期性衛星について、検出されたCNを、1又は複数の恒星日前の記憶されたCNと比較し、所定値以上の差があるか否かを判定する手段と、
CNに所定値以上の差があると判定された周期性衛星を、その後の演算に利用しない衛星として排除する手段と、を有することを特徴とする。
・ 本発明による衛星を用いた測位装置の別の態様は、測定点及び基準点において複数の衛星からそれぞれ受信する電波を用いて、RTK法により測定点の測位を行うべく初期化演算を行う初期化演算部と測位演算を行う測位演算部とを有する測位装置において、
前記初期化演算の前に、又は、前記初期化演算と前記測位演算との間に、その後の演算に利用する衛星を選択するために実行される衛星選択部を有し、前記衛星選択部が、
恒星日毎に周期性のある周期性衛星以外の衛星を、その後の演算に利用しない衛星として排除する手段と、
各前記周期性衛星について、測定点における電波の搬送波電力対雑音電力密度比であるCNを検出して記憶する手段と、
各前記周期性衛星について、検出されたCNを、1又は複数の恒星日前の記憶されたCNと比較し、所定値以上の差があるか否かを判定する手段と、
CNに所定値以上の差があると判定された周期性衛星について、基準点における同一タイミングのCNと、測定点と同じ1又は複数の恒星日前の基準点におけるCNとの差を取得する手段と、
CNに所定値以上の差があると判定された周期性衛星について、当該CNの差が基準点におけるCNの差と同じであるか否かを判定する手段と、
CNに所定値以上の差があると判定された周期性衛星のうち、当該CNの差が基準点におけるCNの差と同じではないと判定された周期性衛星を、その後の演算に利用しない衛星として排除する手段と、を有することを特徴とする。
・ 上記いずれかの装置の態様において、前記衛星選択部が、
測定点の周囲環境情報に基づいて、測定点における排除対象の方位範囲を決定する手段と、
決定された前記方位範囲に含まれる衛星を、その後の演算に利用しない衛星として排除する手段と、をさらに有することが、好適である。
・ 上記いずれかの装置の態様において、前記初期化演算及び前記測位演算を行うために、一定時間の間隔でウェークアップする手段と、
ウェークアップ後に前記初期化演算及び前記測位演算を含む一連の処理を実行し、完了後にスリープ状態とする手段と、をさらに有することが、好適である。
In order to solve the above problems, the present invention provides the following configuration.
An embodiment of a positioning method using satellites according to the present invention is a positioning method for positioning a measurement point by an RTK method having an initialization calculation step and a positioning calculation step using radio waves received from a plurality of satellites at a measurement point and a reference point,
A satellite selection step is executed before the initialization calculation step or between the initialization calculation step and the positioning calculation step to select a satellite to be used in a subsequent calculation, and the satellite selection step includes:
a step of excluding satellites other than the periodic satellites having a periodicity every sidereal day as satellites not to be used in subsequent calculations;
Detecting and storing a carrier power to noise power density ratio (CN0 ) of radio waves at a measurement point for each of the periodic satellites;
for each said periodic satellite, comparing the detected CN0 with a stored CN0 from one or more sidereal days ago to determine whether there is a difference of at least a predetermined value;
and excluding a periodic satellite determined to have a difference in CN0 of a predetermined value or more as a satellite not to be used in subsequent calculations.
Another aspect of the positioning method using satellites according to the present invention is a positioning method for performing positioning of a measurement point by an RTK method having an initialization calculation step and a positioning calculation step using radio waves received from a plurality of satellites at a measurement point and a reference point,
A satellite selection step is executed before the initialization calculation step or between the initialization calculation step and the positioning calculation step to select a satellite to be used in a subsequent calculation, and the satellite selection step includes:
a step of excluding satellites other than the periodic satellites having a periodicity every sidereal day as satellites not to be used in subsequent calculations;
Detecting and storing a carrier power to noise power density ratio (CN0 ) of radio waves at a measurement point for each of the periodic satellites;
for each said periodic satellite, comparing the detected CN0 with a stored CN0 from one or more sidereal days ago to determine whether there is a difference of at least a predetermined value;
For a periodic satellite determined to have a difference of a predetermined value or more in CN 0 , a step of acquiring a difference between CN 0 at the same timing at a reference point and CN 0 at a reference point one or more sidereal days ago that is the same as the measurement point;
For a periodic satellite determined to have a difference in CN0 of a predetermined value or more, determining whether or not the difference in CN0 is the same as the difference in CN0 at a reference point;
and excluding, from among the periodic satellites determined to have a difference in CN0 of a predetermined value or more, a periodic satellite whose difference in CN0 is determined to be not the same as the difference in CN0 at the reference point, as a satellite not to be used in subsequent calculations.
In any of the above methods, the satellite selection step comprises:
determining a range of directions to be excluded at the measurement point based on ambient environment information of the measurement point;
It is preferable that the method further comprises the step of excluding satellites that fall within the determined azimuth range as satellites that will not be used in subsequent calculations.
In any of the above-described methods, it is preferable that the initialization calculation step and the positioning calculation step are performed at regular time intervals.
In any one of the above methods, the positioning calculation step includes:
performing positioning calculations using radio waves at predetermined time intervals within a predetermined time range for all combinations including five or more periodic satellites among the periodic satellites to be used;
A step of calculating a displacement amount, which is a difference from one sidereal day ago, for the positioning calculation result for each of the combinations for each time interval, and calculating an average of the displacement amounts within the time range;
determining whether or not the majority of combinations have a displacement amount equal to or greater than a threshold;
It is preferable to further include a step of determining that there is no displacement when the number of combinations in which the displacement amount is equal to or greater than the threshold is not the majority or more.
In any of the above methods, when the majority or more of the combinations have a displacement amount equal to or greater than the threshold, further determining whether or not only a combination including a specific periodic satellite exceeds the threshold;
It is preferable to have a step of determining that there is a displacement if only the combination including a specific periodic satellite does not exceed the threshold.
An embodiment of a positioning device using satellites according to the present invention is a positioning device having an initialization calculation unit that performs initialization calculations to perform positioning of a measurement point by the RTK method using radio waves received from a plurality of satellites at a measurement point and a reference point, and a positioning calculation unit that performs positioning calculations,
A satellite selection unit is provided for selecting a satellite to be used in a subsequent calculation before an initialization calculation or between the initialization calculation and the positioning calculation, the satellite selection unit being:
a means for excluding satellites other than a periodic satellite having a periodicity every sidereal day as satellites not to be used in subsequent calculations;
A means for detecting and storing a carrier power to noise power density ratio (CN0 ) of radio waves at a measurement point for each of the periodic satellites;
means for comparing, for each of said periodic satellites, the detected CN0 with a stored CN0 from one or more sidereal days ago to determine whether there is a difference between the detected CN0 and the stored CN0 by more than a predetermined value;
and means for excluding a periodic satellite determined to have a difference in CN0 of a predetermined value or more as a satellite not to be used in subsequent calculations.
Another aspect of the satellite-based positioning device according to the present invention is a positioning device having an initialization calculation unit that performs initialization calculations to perform positioning of a measurement point by the RTK method using radio waves received from a plurality of satellites at a measurement point and a reference point, and a positioning calculation unit that performs positioning calculations,
A satellite selection unit is executed before the initialization operation or between the initialization operation and the positioning operation to select a satellite to be used in a subsequent operation, the satellite selection unit being:
a means for excluding satellites other than a periodic satellite having a periodicity every sidereal day as satellites not to be used in subsequent calculations;
A means for detecting and storing a carrier power to noise power density ratio (CN0 ) of radio waves at a measurement point for each of the periodic satellites;
means for comparing, for each of said periodic satellites, the detected CN0 with a stored CN0 from one or more sidereal days ago to determine whether there is a difference between the detected CN0 and the stored CN0 by more than a predetermined value;
For a periodic satellite determined to have a difference of a predetermined value or more in CN 0 , a means for acquiring a difference between CN 0 at the same timing at a reference point and CN 0 at a reference point one or more sidereal days ago that is the same as the measurement point;
A means for determining whether or not the difference in CN0 of a periodic satellite determined to have a difference of a predetermined value or more is the same as the difference in CN0 at a reference point;
and a means for excluding, from among the periodic satellites determined to have a difference in CN0 of a predetermined value or more, a periodic satellite whose difference in CN0 is determined to be not the same as the difference in CN0 at the reference point, as a satellite not to be used in subsequent calculations.
In any of the above-mentioned device aspects, the satellite selection unit
A means for determining a range of directions to be excluded at the measurement point based on ambient environment information of the measurement point;
It is preferable that the satellite detection apparatus further comprises a means for excluding satellites included in the determined azimuth range as satellites not to be used in subsequent calculations.
In any of the above aspects of the device, a means for waking up at a constant time interval in order to perform the initialization calculation and the positioning calculation;
It is preferable that the device further comprises a means for executing a series of processes including the initialization calculation and the positioning calculation after waking up, and going into a sleep state after completion of the processes.

本発明によれば、GNSSを用いた測位方法及び測位装置において、RTK法による演算に利用する衛星の中から排除すべき衛星を効率的に判別できる。 According to the present invention, in a positioning method and a positioning device using GNSS, it is possible to efficiently identify satellites to be excluded from among the satellites used in calculations using the RTK method.

図1は、本発明の第1の実施形態の測位装置を含む測位システム全体の一例を概略的に示す構成図である。FIG. 1 is a diagram showing an example of a configuration of an overall positioning system including a positioning device according to a first embodiment of the present invention. 図2は、第1の実施形態の測位方法の一例を概略的に示したフロー図である。FIG. 2 is a flow diagram illustrating an example of a positioning method according to the first embodiment. 図3は、図2の続きのフロー図である。FIG. 3 is a continuation of the flow chart of FIG. 図4は、図2及び図3に示した測位方法の変形例を概略的に示すフロー図である。FIG. 4 is a flow diagram that illustrates a schematic modification of the positioning method shown in FIGS. 図5は、本発明の第2の実施形態の測位装置を含む測位システム全体の一例を概略的に示す構成図である。FIG. 5 is a diagram showing an example of a configuration of an overall positioning system including a positioning device according to the second embodiment of the present invention. 図6は、第2の実施形態における測位方法の一例を概略的に示したフロー図である。FIG. 6 is a flow diagram illustrating an example of a positioning method according to the second embodiment. 図7は、本発明の第3の実施形態の測位装置を含む測位システム全体の一例を概略的に示す構成図である。FIG. 7 is a diagram showing an example of a configuration of an overall positioning system including a positioning device according to the third embodiment of the present invention. 図8は、第3の実施形態の測位方法の一例を概略的に示したフロー図である。FIG. 8 is a flow diagram illustrating an example of a positioning method according to the third embodiment. 図9は、本発明の第4の実施形態の測位装置を含む測位システム全体の一例を概略的に示す構成図である。FIG. 9 is a diagram showing an example of a configuration of an overall positioning system including a positioning device according to the fourth embodiment of the present invention.

以下、実施例を示した図面を参照しつつ本発明の実施形態を説明する。図面においては、各実施形態における同じ又は類似の構成要素に同じ符号を付している。 Embodiments of the present invention will be described below with reference to the drawings showing examples. In the drawings, the same or similar components in each embodiment are given the same reference numerals.

測位法の一つであるRTK法では、測定点で受信される電波と基準点(既知点)で受信される電波の行路差に対応する整数値バイアスを確定する初期化演算ステップと、その後、測定点の位置を計測する測位演算ステップの2段階の演算工程から構成される。 The RTK method, which is one of the positioning methods, consists of a two-step calculation process: an initialization calculation step that determines an integer bias corresponding to the path difference between the radio waves received at the measurement point and the radio waves received at a reference point (known point), and then a positioning calculation step that measures the position of the measurement point.

(1)第1の実施形態
図1は、本発明の第1の実施形態の測位装置を含む測位システム全体の一例を概略的に示す構成図である。測定点に設けられる移動局システム1は、アンテナ10、受信装置20及び移動局端末30を有する。ここでは、移動局端末30が本発明の測位装置に相当する。
(1) First embodiment Fig. 1 is a schematic diagram showing an example of an overall positioning system including a positioning device according to a first embodiment of the present invention. A mobile station system 1 provided at a measurement point has an antenna 10, a receiving device 20, and a mobile station terminal 30. Here, the mobile station terminal 30 corresponds to the positioning device of the present invention.

移動局システム1におけるアンテナ10は、測定点に設置される。アンテナ10は、例えば、高周波増幅器を内蔵し、測定点において複数の衛星S1、S2..から電波(搬送波)を受信し、増幅した後、アンテナ10に接続された受信装置20に送信する。受信装置20は、例えば、高周波増幅回路、コード復調回路、デジタル信号処理回路等から構成されており、受信した電波に所定の付加データを追加して、移動局端末30に送信する。一例として、受信装置20を移動体端末30に組み込み、一体の装置とすることもできる。 The antenna 10 in the mobile station system 1 is installed at the measurement point. The antenna 10, for example, has a built-in high-frequency amplifier, receives radio waves (carrier waves) from multiple satellites S1, S2, etc. at the measurement point, amplifies the waves, and then transmits them to the receiving device 20 connected to the antenna 10. The receiving device 20 is composed of, for example, a high-frequency amplifier circuit, a code demodulation circuit, a digital signal processing circuit, etc., and adds predetermined additional data to the received radio waves and transmits them to the mobile station terminal 30. As an example, the receiving device 20 can be incorporated into the mobile terminal 30 to form an integrated device.

移動局端末30は、中央処理装置31を備える。中央処理装置31は、図示しないが、CPUとRAM等のメモリを有し、CPUが所定のアプリケーションをメモリに読み込み実行することによって、本発明の測位装置の各機能を実現する演算処理を行う。中央処理装置31は、主な機能処理部として、衛星選択部1A、初期化演算部1B、測位演算部1C及び変位判定部1Dを有する。 The mobile station terminal 30 is equipped with a central processing unit 31. Although not shown, the central processing unit 31 has a CPU and memory such as RAM, and the CPU loads and executes a specific application into the memory to perform calculation processing that realizes each function of the positioning device of the present invention. The central processing unit 31 has, as its main functional processing units, a satellite selection unit 1A, an initialization calculation unit 1B, a positioning calculation unit 1C, and a displacement determination unit 1D.

衛星選択部1Aは、その後の初期化演算部1B及び測位演算部1Cにおける演算に利用する衛星を適切に選択する機能を有する。初期化演算部1Bは、利用する衛星の電波を用いて初期化演算を行ってFix解を得る機能を有する。測位演算部1Cは、利用する衛星の電波を用いて測定点の位置座標を演算する機能を有する。変位判定部1Dは、測定点の位置の変位の有無を判定する機能を有する。各機能処理部の詳細は、後述する図2及び図3を参照して測位方法と共に説明する。 The satellite selection unit 1A has a function of appropriately selecting a satellite to be used for subsequent calculations in the initialization calculation unit 1B and the positioning calculation unit 1C. The initialization calculation unit 1B has a function of performing initialization calculations using radio waves from the satellite to be used to obtain a fix solution. The positioning calculation unit 1C has a function of calculating the position coordinates of the measurement point using radio waves from the satellite to be used. The displacement determination unit 1D has a function of determining whether or not the position of the measurement point has been displaced. Details of each functional processing unit will be described below with reference to Figures 2 and 3, along with the positioning method.

移動局端末30はさらに、好ましくは無線通信機能を備えた通信装置32、処理状況及び処理結果等を表示する表示装置33、受信した電波及び演算結果データ等の種々のデータを保存するための記憶装置34を有する。これらの各装置32、33、34も中央処理装置31によって制御される。 The mobile station terminal 30 further includes a communication device 32, preferably equipped with wireless communication capabilities, a display device 33 for displaying the processing status and processing results, and a storage device 34 for storing various data such as received radio waves and calculation result data. Each of these devices 32, 33, and 34 is also controlled by the central processing unit 31.

移動局端末30は、本発明の機能のみを組み込んだ専用装置として実施できる。しかしながら、例えばパーソナルコンピュータ、タブレットコンピュータ、スマートフォン等に所定のアプリケーションを導入することによっても実施することができる。移動局端末30は、ユーザによって容易に携帯可能であることが好ましい。 The mobile station terminal 30 can be implemented as a dedicated device incorporating only the functions of the present invention. However, it can also be implemented by installing a specific application on, for example, a personal computer, a tablet computer, a smartphone, etc. It is preferable that the mobile station terminal 30 is easily portable by the user.

一方、基準点に設けられる基準局システム5は、アンテナ51と、受信装置52と、通信装置53とを有する。 On the other hand, the reference station system 5 installed at the reference point has an antenna 51, a receiving device 52, and a communication device 53.

基準点におけるアンテナ51は、位置座標が既に判明している基準点に設置される。アンテナ51も、例えば高周波増幅器を内蔵し、基準点において複数の衛星S1、S2..から電波(搬送波)を受信し、増幅した後、アンテナ51に接続された受信装置52に送信する。受信装置52は、受信装置20と同様の構成を有しており、受信した電波に所定のデータを付加して通信装置53に送信する。基準点の通信装置53と、測定点の移動局端末30の通信装置32は、例えば無線通信ネットワークによってデータ伝送可能に接続されており、基準点の電波はリアルタイムで移動局端末30に送信することができる。 An antenna 51 at the reference point is installed at a reference point whose position coordinates are already known. The antenna 51 also has a built-in high-frequency amplifier, for example, and receives radio waves (carrier waves) from multiple satellites S1, S2, etc. at the reference point, amplifies the waves, and then transmits them to a receiving device 52 connected to the antenna 51. The receiving device 52 has a configuration similar to that of the receiving device 20, and adds predetermined data to the received radio waves and transmits them to a communication device 53. The communication device 53 at the reference point and the communication device 32 of the mobile station terminal 30 at the measurement point are connected to each other, for example, via a wireless communication network, so that data can be transmitted, and the radio waves at the reference point can be transmitted to the mobile station terminal 30 in real time.

図2及び図3は、第1の実施形態の測位方法の一例を概略的に示したフロー図である。
ステップ101において、移動局端末30がウェークアップされる。斜面崩壊の監視等のために使用される場合は電源確保が困難であることが多いため、移動局端末30は、バッテリを電源とすることが好ましい。したがって、電力消費を抑制するために、移動局端末30が、連続的ではなく断続的に測位を行うことが好ましい。移動局端末30は、例えば30分毎にウェークアップされ、測位の一連の手順を完了したならばスリープ状態とされる。スリープ状態では、ウェークアップ用の制御機能のみを稼動させ、それ以外の機能を停止させる。斜面崩壊の監視の場合、ウェークアップの間隔は60分毎が好ましい。RTK法による測位の一連の手順は、初期化を1分以内、その後の測位計測を3分以内で完了することが好ましい。なお、ウェークアップの間隔は60分毎に限られず、必要に応じて、例えば30分毎など、一定時間の間隔とすることができる。電力供給に支障がない場合は、連続的に測位することもできる。図示しないが、移動局端末30と通信ネットワークを介して通信可能な遠隔端末から移動局端末30をウェークアップさせることも可能である。そのような遠隔端末から移動局端末30に制御信号を送信することにより、ウェークアップ間隔を設定及び変更することも可能である。
2 and 3 are flow charts each showing an outline of an example of the positioning method according to the first embodiment.
In step 101, the mobile station terminal 30 is woken up. When the mobile station terminal 30 is used for monitoring slope collapse, etc., it is often difficult to secure a power source, so it is preferable that the mobile station terminal 30 uses a battery as a power source. Therefore, in order to suppress power consumption, it is preferable that the mobile station terminal 30 performs positioning intermittently, not continuously. The mobile station terminal 30 is woken up, for example, every 30 minutes, and is put into a sleep state after completing a series of positioning procedures. In the sleep state, only the control function for waking up is operated, and other functions are stopped. In the case of monitoring slope collapse, it is preferable that the wake-up interval is every 60 minutes. It is preferable that the series of positioning procedures using the RTK method are initialized within 1 minute, and the subsequent positioning measurement is completed within 3 minutes. Note that the wake-up interval is not limited to every 60 minutes, and can be set to a fixed time interval, such as every 30 minutes, as necessary. If there is no problem with the power supply, it is also possible to perform positioning continuously. Although not shown, it is also possible to wake up the mobile station terminal 30 from a remote terminal that can communicate with the mobile station terminal 30 via a communication network. It is also possible to set and change the wake-up interval by transmitting a control signal from such a remote terminal to the mobile station terminal 30.

ウェークアップ後、移動局端末30は、そのタイミングにアンテナ10が受信する全ての衛星からの電波の受信を開始すると共に、基準局システム5から送信される電波の受信を開始する。基準局システム5からの電波は、同一タイミングに同じ全ての衛星からアンテナ51が受信したものである。なお、本明細書における「タイミング」の用語は、衛星の周期中の一つの時点を意味する。 After waking up, the mobile station terminal 30 starts receiving radio waves from all satellites that the antenna 10 receives at that timing, and also starts receiving radio waves transmitted from the reference station system 5. The radio waves from the reference station system 5 are received by the antenna 51 from all the same satellites at the same timing. Note that the term "timing" in this specification means a point in time during a satellite cycle.

ステップ102~105の処理は、衛星選択部1Aによって行われる。衛星選択部1Aは、その後の初期化演算及び測位演算においてその電波を利用すべきでない衛星を排除するために行われる。これは、ステップ106、107の初期化演算ステップの前に行われる。利用すべきでない電波とは、ノイズが多く質の悪い電波である。そのような電波は、初期化完了の遅延又は不能、ミスFix、位置精度の低下などの要因となる。 The processing of steps 102 to 105 is performed by the satellite selection unit 1A. The satellite selection unit 1A performs this in order to eliminate satellites whose radio waves should not be used in the subsequent initialization calculations and positioning calculations. This is performed before the initialization calculation steps of steps 106 and 107. Radio waves that should not be used are radio waves with a lot of noise and poor quality. Such radio waves can cause delays or inability to complete initialization, missed fixes, and reduced position accuracy.

ステップ102において、恒星日毎に周期性のある複数の衛星(以下「周期性衛星」と称する)以外の衛星を、その後の演算に利用しない衛星として排除する。周期性衛星は、毎恒星日の周期中の同一タイミングに同じ位置となる衛星である。本発明の測位方法では、周期性衛星のみを利用してRTK法を適用する。 In step 102, satellites other than a number of satellites that are periodic every sidereal day (hereinafter referred to as "periodic satellites") are excluded as satellites that will not be used in subsequent calculations. A periodic satellite is a satellite that is in the same position at the same time during each sidereal day cycle. In the positioning method of the present invention, the RTK method is applied using only periodic satellites.

GPS、QZSSなどのGNSSの衛星は、1恒星日毎に同じ位置に配置される。また、静止衛星(SBASや準天頂衛星3号機、BeiDouのGEO等)も、常に同じ位置に配置されるので周期性衛星に含まれるものとする。アンテナが1恒星日前の周期中の同一タイミングに同じ位置にあれば、これらの衛星からの電波に対するマルチパスや障害物による影響に関して同じ状況が再現されるはずである。したがって、電波を取得する毎に、その電波を、受信したタイミングの、例えば1恒星日前の電波と比較することによって、条件の変化した衛星を判別することができる。別の例として、1恒星日前に替えて、複数恒星日前とすることもできる。例えば、2恒星日前、7恒星日前などである。あるいは、測定点に移動局を設置した当初に、連続した1恒星日に亘って電波を受信し記憶しておき、その記憶された電波の連続データのうちの同じタイミングの電波を比較基準として用いることもできる。以下、本例では、1恒星日前の電波と比較する場合について説明する。ここで、搬送波電力対雑音電力密度比CNは、衛星からの電波の強度を表しており、障害物の変動によるマルチパスの変化や電波の遮断があると変化する。本発明では、CNを電波に対する影響の評価パラメータとして利用する。 GNSS satellites such as GPS and QZSS are located at the same position every sidereal day. Geostationary satellites (SBAS, Quasi-Zenith Satellite 3, BeiDou GEO, etc.) are also considered to be periodic satellites because they are always located at the same position. If an antenna is located at the same position at the same timing during the period one sidereal day ago, the same situation regarding the influence of multipath and obstacles on radio waves from these satellites should be reproduced. Therefore, each time radio waves are acquired, the radio waves can be compared with the radio waves received at the timing, for example, one sidereal day ago, to identify satellites whose conditions have changed. As another example, multiple sidereal days ago can be used instead of one sidereal day ago. For example, two sidereal days ago, seven sidereal days ago, etc. Alternatively, radio waves can be received and stored for one consecutive sidereal day at the beginning of installing a mobile station at a measurement point, and radio waves at the same timing among the continuous data of the stored radio waves can be used as a comparison standard. In the following, in this example, a comparison with radio waves one sidereal day ago will be described. Here, the carrier power to noise power density ratio CN0 represents the strength of radio waves from a satellite, and changes when there is a change in the multipath due to the movement of obstacles or when radio waves are blocked. In the present invention, CN0 is used as an evaluation parameter for the influence on radio waves.

ステップ103において、各周期性衛星の電波に基づいてCNを検出する。そして検出されたCNを、電波を受信したタイミングと対応付けて記憶装置34に記憶する。CNは、以下のように定義される。
CN=C/N
C:搬送波電力
:雑音電力密度
単位:dBHz
In step 103, CN 0 is detected based on the radio waves of each periodic satellite. Then, the detected CN 0 is stored in the storage device 34 in association with the timing at which the radio waves were received. CN 0 is defined as follows.
CN0 = C/ N0
C: Carrier power
N 0 : Noise power density
Unit: dBHz

ステップ104において、ステップ103で検出された各周期性衛星のCNを、1恒星日前の当該周期性衛星のCNと比較する。1恒星日前のCNは、記憶装置34に記憶されている。比較によって、双方の間に所定値以上の差があるか否かを、例えば1dBHz以上の差があるか否かを判定する。所定値以上の差が無い場合は、そのまま初期化演算に進む。 In step 104, the CN 0 of each periodic satellite detected in step 103 is compared with the CN 0 of the corresponding periodic satellite one sidereal day ago. The CN 0 of one sidereal day ago is stored in the storage device 34. By the comparison, it is determined whether there is a difference between the two that is equal to or greater than a predetermined value, for example, whether there is a difference of 1 dBHz or more. If there is no difference that is equal to or greater than the predetermined value, the process proceeds directly to the initialization calculation.

ステップ104でCNに所定値以上の差がある場合は、その周期性衛星を、その後の演算に利用しない衛星として排除する。そして、残りの周期性衛星の電波を用いた初期化演算に進む。 If there is a difference of a predetermined value or more in CN0 in step 104, that periodic satellite is excluded as a satellite that will not be used in the subsequent calculations, and the process proceeds to initialization calculations using the radio waves of the remaining periodic satellites.

1恒星日の前後におけるCNの差を評価パラメータとして衛星選択を行うことによって、利用可能な周期性衛星を初期化演算の前に選択することができる。これによって、初期化演算におけるミスFixの回避及び初期化時間の短縮が可能となる。さらに、精度悪化の要因となる衛星を排除することによって、その後の測位演算結果のばらつきを予め抑制することができる。1恒星日前からCNに変化のない衛星を利用して測位することによって、アンテナ10の位置及び変位した距離を高精度に計測することができる。 By selecting satellites using the difference in CN 0 before and after one sidereal day as an evaluation parameter, available periodic satellites can be selected before the initialization calculation. This makes it possible to avoid misfixes in the initialization calculation and shorten the initialization time. Furthermore, by eliminating satellites that cause accuracy degradation, it is possible to suppress in advance the variation in the positioning calculation results thereafter. By performing positioning using satellites whose CN 0 has not changed since one sidereal day ago, the position of the antenna 10 and the displaced distance can be measured with high accuracy.

ステップ106、107の処理は、初期化演算部1Bによって行われる。ステップ106において、衛星選択部1Aによって選択された、利用する全ての周期性衛星の電波を用いて初期化演算を行う。ステップ107においてRTK法におけるFix解を得る。移動局端末30が30分毎にウェークアップする場合、初期化演算において1分以内にFix解を得ることが好ましい。初期化後に数cmの精度となったところで次のステップの測位計測を行う。 The processes of steps 106 and 107 are performed by the initialization calculation unit 1B. In step 106, the initialization calculation is performed using radio waves from all periodic satellites to be used that are selected by the satellite selection unit 1A. In step 107, a fix solution is obtained in the RTK method. If the mobile station terminal 30 wakes up every 30 minutes, it is preferable to obtain a fix solution within one minute in the initialization calculation. When an accuracy of several centimeters is achieved after initialization, the next step of positioning measurement is performed.

ステップ108、109は、測位演算部1Cによって行われる。ステップ108、109において測定点における測位計測を行う。ステップ108では、衛星選択部1Aによって選択された(すなわち初期化演算に用いた)、利用する全ての周期性衛星の電波(ここでは「計測データ」と称する)を取得する。例えば、1秒毎に計測データを取得し、3分間の計測を行う。n個の周期性衛星を利用する場合、1秒毎にn個の計測データが得られ、3分間にn×180個の計測データが得られる。計測データは、電波を受信したタイミングと対応付けられて記憶装置34に記憶される。 Steps 108 and 109 are performed by the positioning calculation unit 1C. In steps 108 and 109, positioning measurements are performed at the measurement point. In step 108, radio waves (herein referred to as "measurement data") from all periodic satellites to be used that have been selected by the satellite selection unit 1A (i.e., used in the initialization calculation) are acquired. For example, measurement data is acquired every second, and measurements are performed for three minutes. When n periodic satellites are used, n pieces of measurement data are obtained every second, and n x 180 pieces of measurement data are obtained in three minutes. The measurement data is stored in the storage device 34 in association with the timing at which the radio waves were received.

ステップ109では、ステップ108で得られた計測データを用いた測位演算を行って位置座標を算出する。好ましくは、n個の周期性衛星を利用する場合、5個以上の全ての組合せについて測位演算を行う。5個以上の全ての組合せの数は、nCk(kは5以上n以下の自然数)個である。この結果、各組合せについて、3分間の1秒毎の位置座標が得られる。この段階で得られた位置座標は、測定点の仮の位置座標として記憶装置34に保存しておく。 In step 109, a positioning calculation is performed using the measurement data obtained in step 108 to calculate the position coordinates. Preferably, when n periodic satellites are used, the positioning calculation is performed for all combinations of 5 or more. The total number of combinations of 5 or more is nCk (k is a natural number between 5 and n). As a result, position coordinates are obtained for each combination every second for 3 minutes. The position coordinates obtained at this stage are stored in the storage device 34 as provisional position coordinates of the measurement point.

図3のステップ110~117は、変位判定部1Dによって行われる。ステップ110~117において、測定点における1恒星日前からの変位を判定する。ステップ110では、衛星の各組合せについて、1秒毎の位置座標の、1恒星日前の周期中の同一タイミングの位置座標との差分すなわち変位量を計算する。すなわち、同じ衛星の組合せにおける1恒星日前の位置座標との差分を導出する。同じ衛星の組合せとすることによってマルチパスの影響を除去して変位のみを抽出することができる。ステップ111では、各組合せについて、3分間における1秒毎の変位量の平均を計算する。 Steps 110 to 117 in FIG. 3 are performed by the displacement determination unit 1D. In steps 110 to 117, the displacement at the measurement point from one sidereal day ago is determined. In step 110, for each combination of satellites, the difference between the position coordinates for each second and the position coordinates at the same time during the period one sidereal day ago, i.e., the amount of displacement, is calculated. In other words, the difference from the position coordinates for the same combination of satellites one sidereal day ago is derived. By using the same combination of satellites, it is possible to remove the effects of multipath and extract only the displacement. In step 111, the average amount of displacement for each combination for each second over a three-minute period is calculated.

ステップ112において、全ての組合せの各々の変位量を閾値(例えば3mm)と比較し、変位量が閾値を超えている組合せが過半数以上であるか否かを判定する。閾値を超えている組合せが過半数以上でない場合は、ステップ113において「変位なし」と判定すると共に、上記ステップ109で得た各組合せの仮の位置座標から、3分間の時間的な平均及び全ての組合せの平均を計算するなどによって、測定点の位置座標を確定する。 In step 112, the amount of displacement for each of all combinations is compared with a threshold value (e.g., 3 mm) to determine whether or not the majority of combinations have a displacement amount that exceeds the threshold value. If the majority of combinations do not have a displacement amount that exceeds the threshold value, in step 113, it is determined that there is "no displacement," and the position coordinates of the measurement points are determined by calculating the temporal average over three minutes and the average of all combinations from the provisional position coordinates of each combination obtained in step 109 above.

一方、ステップ112において、変位量が閾値を超えている組合せが過半数以上である場合、さらにステップ114において、特定の衛星を含む組合せのみが閾値を超えているか否かを判定する。特定の衛星を含む組合せのみが閾値を超えている場合は、当該特定の衛星は測位に利用するには不適切な衛星であると判定し、ステップ115において当該特定の衛星を含む組合せを排除する。そして、ステップ116において、再度、変位量が閾値を超えている組合せが過半数以上であるか否かを判定する。 On the other hand, if the majority of combinations have a displacement amount exceeding the threshold in step 112, then in step 114 it is determined whether only combinations including a specific satellite exceed the threshold. If only combinations including a specific satellite exceed the threshold, it is determined that the specific satellite is an unsuitable satellite for use in positioning, and in step 115 the combinations including the specific satellite are eliminated. Then, in step 116 it is determined again whether the majority of combinations have a displacement amount exceeding the threshold.

ステップ116において、変位量が閾値を超えている組合せが過半数以上である場合、ステップ117において「変位あり」と判定する。そして、上記ステップ109で得た各組合せの仮の位置座標から、この段階で残っている組合せのみを用いて測定点の位置座標を確定すると共に、測定点の変位量も確定する。 In step 116, if the majority of combinations have a displacement amount exceeding the threshold, it is determined in step 117 that "displacement exists." Then, from the provisional position coordinates of each combination obtained in step 109, the position coordinates of the measurement points are determined using only the combinations remaining at this stage, and the displacement amount of the measurement points is also determined.

ステップ114おいて、特定の衛星を含む組合せのみが閾値を超えているのではない場合、ステップ117において「変位あり」と判定し、測定点の位置座標及び変位量を確定する。 If in step 114 only the combination including a specific satellite does not exceed the threshold, then in step 117 it is determined that "displacement exists" and the position coordinates and amount of displacement of the measurement point are determined.

ステップ116において、変位量が閾値を超えている組合せが過半数以上でない場合は、ステップ113において「変位なし」と判定すると共に測定点の位置座標を確定する。 If the majority of combinations do not have a displacement amount exceeding the threshold in step 116, then in step 113 it is determined that there is no displacement and the position coordinates of the measurement point are determined.

測定点の位置座標の確定、又は、測定点の位置座標及び変位量の確定後、移動局端末30をスリープ状態とする。 After the position coordinates of the measurement point are determined, or the position coordinates of the measurement point and the amount of displacement are determined, the mobile station terminal 30 is put into a sleep state.

図4は、図2及び図3に示した測位方法の変形例を概略的に示すフロー図である。図2及び図3のフロー図とは異なる点のみを説明する。図4の測位方法では、図1の衛星選択部1Aと初期化演算部1Bの各々による処理を、図2の測位方法とは逆の順序で行う。すなわち、ステップ401、402において、先ず、受信した全ての衛星の電波を用いて初期化演算を行ってRTKのFix解を得る。その後、ステップ403~406において、図2のステップ102~105と同様にCNを評価パラメータとして衛星選択を行う。それ以降の測位演算部1C及び変位判定部1Dによる処理は、図2及び図3と同様である。 FIG. 4 is a flow diagram that shows a schematic modification of the positioning method shown in FIG. 2 and FIG. 3. Only the points different from the flow diagrams of FIG. 2 and FIG. 3 will be described. In the positioning method of FIG. 4, the processes by the satellite selection unit 1A and the initialization calculation unit 1B of FIG. 1 are performed in the reverse order to the positioning method of FIG. 2. That is, in steps 401 and 402, first, an initialization calculation is performed using radio waves of all received satellites to obtain a fix solution of RTK. Then, in steps 403 to 406, satellite selection is performed with CN 0 as an evaluation parameter, as in steps 102 to 105 of FIG. 2. The subsequent processes by the positioning calculation unit 1C and the displacement determination unit 1D are the same as those in FIG. 2 and FIG. 3.

図4の変形例の測位方法では、初期化演算ステップと測位演算ステップの間に衛星選択を行うことによって、その後の測位演算に利用可能な周期性衛星を、測位演算ステップの前に選択することができる。これによって、その後の測位演算結果のばらつきを予め抑制することができる。 In the modified positioning method of FIG. 4, satellite selection is performed between the initialization calculation step and the positioning calculation step, so that a periodic satellite that can be used for the subsequent positioning calculation can be selected before the positioning calculation step. This makes it possible to suppress in advance the variation in the results of the subsequent positioning calculation.

図示しないが、本発明の測位方法は、さらに従来技術である仰角マスク、CNマスク、衛星システムの選択、測位開始タイミングの選択等を用いて不適切な衛星の排除を行うことを組み合わせることができる(以下の実施形態においても同様)。これらの従来技術は、初期化演算ステップ又は測位演算ステップの開始後又は完了後に、演算状況及び演算結果に基づいて実行される。 Although not shown, the positioning method of the present invention can further be combined with the conventional techniques of excluding inappropriate satellites using an elevation mask, a CN0 mask, selection of a satellite system, selection of a positioning start timing, etc. (The same applies to the following embodiments.) These conventional techniques are executed based on the calculation status and calculation results after the start or completion of the initialization calculation step or the positioning calculation step.

(2)第2の実施形態
図5は、本発明の第2の実施形態の測位装置を含む測位システム全体の一例を概略的に示す構成図である。測定点に設けられる移動局システム1は、アンテナ10、受信装置20及び移動局端末30を有する。ここでも、移動局端末30が本発明の測位装置に相当する。移動局システム1の構成は、上述した第1の実施形態とほぼ同じであるので説明を省略する。第2の実施形態では、第1の実施形態とは中央処理装置31の衛星選択部1Aにおける処理が異なるが、これについては、図6を参照して測位方法と共に説明する。
(2) Second embodiment Fig. 5 is a schematic diagram showing an example of an overall positioning system including a positioning device according to a second embodiment of the present invention. A mobile station system 1 provided at a measurement point has an antenna 10, a receiving device 20, and a mobile station terminal 30. Again, the mobile station terminal 30 corresponds to the positioning device of the present invention. The configuration of the mobile station system 1 is almost the same as that of the first embodiment described above, so a description thereof will be omitted. In the second embodiment, the processing in the satellite selection unit 1A of the central processing unit 31 is different from that of the first embodiment, but this will be described together with the positioning method with reference to Fig. 6.

一方、基準点に設けられる基準局システム5は、アンテナ51と、受信装置52と、基準局端末60とを有する。 On the other hand, the reference station system 5 installed at the reference point has an antenna 51, a receiving device 52, and a reference station terminal 60.

基準点におけるアンテナ51及び受信装置52は、第1の実施形態と同じであるので説明を省略する。但し、受信装置52は、受信した電波を基準局端末60に送信する。 The antenna 51 and the receiving device 52 at the reference point are the same as those in the first embodiment, so a description is omitted. However, the receiving device 52 transmits the received radio waves to the reference station terminal 60.

基準局端末60は、中央処理装置61と、記憶装置62と、通信装置63とを備える。中央処理装置61は、図示しないが、CPUとRAM等のメモリを有し、CPUが所定のアプリケーションをメモリに読み込み実行することによって、CN検出部6A、CN記録部6Bの各機能処理部を実現する。CN検出部6Aは、受信装置52から受信した基準点における各衛星の電波のCNを検出する。CN記録部6Bは、各衛星について検出されたCNを、電波を受信したタイミングと対応付けて記憶装置62に記憶する。 The reference station terminal 60 comprises a central processing unit 61, a storage device 62, and a communication device 63. The central processing unit 61 has a CPU and memory such as RAM, not shown, and the CPU loads a predetermined application into the memory and executes it to realize each functional processing unit of a CN0 detection unit 6A and a CN0 recording unit 6B. The CN0 detection unit 6A detects the CN0 of the radio waves of each satellite at the reference point received from the receiving device 52. The CN0 recording unit 6B stores the CN0 detected for each satellite in the storage device 62 in association with the timing of receiving the radio waves.

基準局端末60の通信装置63は、中央処理装置61の制御にしたがって、基準点における電波を移動局端末30の通信装置32に送信する。その際、通信装置63は、電波と共に、当該電波のCNと当該電波の例えば1恒星日前のCNとを送信する。なお、第1の実施形態で説明したように、1恒星日前のCNとすることは一例であり、複数恒星日前のCNでもよい。 The communication device 63 of the reference station terminal 60 transmits radio waves at the reference point to the communication device 32 of the mobile station terminal 30 under the control of the central processing unit 61. At that time, the communication device 63 transmits the CN 0 of the radio waves and, for example, the CN 0 of the radio waves one sidereal day ago together with the radio waves. Note that, as explained in the first embodiment, the CN 0 of one sidereal day ago is only an example, and the CN 0 of several sidereal days ago may be used.

別の例として、中央処理装置61が、電波のCNを検出してその電波の1恒星日前の同一タイミングのCNとの差を計算し、計算されたCNの差を、通信装置63から移動局端末30の通信装置32に送信するように構成してもよい。 As another example, the central processing unit 61 may be configured to detect the CN 0 of a radio wave, calculate the difference between that radio wave and the CN 0 at the same timing one sidereal day earlier, and transmit the calculated difference in CN 0 from the communication device 63 to the communication device 32 of the mobile station terminal 30.

基準局端末60は、上記の機能のみを組み込んだ専用装置として実施できる。しかしながら、例えばパーソナルコンピュータ、タブレットコンピュータ、スマートフォン等に所定のアプリケーションを導入することによっても実施することができる。 The reference station terminal 60 can be implemented as a dedicated device incorporating only the above functions. However, it can also be implemented by installing a specific application on, for example, a personal computer, tablet computer, smartphone, etc.

図6は、第2の実施形態における測位方法の一例を概略的に示したフロー図である。図6では、主として第1の実施形態とは異なる部分を示している。ウェークアップは、図2のステップ101と同様である。 Figure 6 is a flow diagram that shows an example of a positioning method in the second embodiment. Figure 6 mainly shows the parts that are different from the first embodiment. Wake-up is similar to step 101 in Figure 2.

ステップ601~606の処理は、衛星選択部1Aによって行われ、その後の初期化演算及び測位演算においてその電波を利用すべきでない衛星を排除する。 The processing of steps 601 to 606 is performed by the satellite selection unit 1A, which excludes satellites whose radio waves should not be used in subsequent initialization and positioning calculations.

ステップ601において、周期性衛星以外の衛星を、その後の演算に利用しない衛星として排除する。このステップは、図2のステップ102と同様である。ステップ602において、周期性衛星の各々の電波に基づいてCNを検出し、電波を受信したタイミングと対応付けて記憶装置34に記憶する。このステップは、図2のステップ103と同様である。 In step 601, satellites other than periodic satellites are excluded as satellites that will not be used in the subsequent calculations. This step is the same as step 102 in Fig. 2. In step 602, CN0 is detected based on the radio waves of each periodic satellite, and stored in storage device 34 in association with the timing of receiving the radio waves. This step is the same as step 103 in Fig. 2.

ステップ603において、ステップ602で検出された各周期性衛星のCNを、1恒星日前の当該周期性衛星のCNと比較する。1恒星日前のCNは、記憶装置34に記憶されている。比較によって、双方の間に所定値以上の差があるか否かを、例えば1dBHz以上の差があるか否かを判定する。所定値以上の差が無い場合は、そのまま初期化演算に進む。 In step 603, the CN 0 of each periodic satellite detected in step 602 is compared with the CN 0 of the corresponding periodic satellite one sidereal day ago. The CN 0 of one sidereal day ago is stored in the storage device 34. By the comparison, it is determined whether there is a difference between the two of a predetermined value or more, for example, a difference of 1 dBHz or more. If there is no difference of the predetermined value or more, the process proceeds directly to the initialization calculation.

ステップ603でCNに所定値以上の差がある場合、ステップ604において基準点における当該周期性衛星のCNを取得する。ここで取得される基準点におけるCNは、測定点と同一タイミングのCNと、その1恒星日前のCNである。これらのCNのデータは、基準局端末60の通信装置63から移動局端末30の通信装置32に送信される。衛星選択部1Aは、受信した基準点における同一タイミングのCNと1恒星日前のCNの差を算出する。別の例として、基準点における同一タイミングのCNとその1恒星日前のCNとの差を、基準局端末60において算出し、CNの差のデータを移動局端末30に送信してもよい。 If there is a difference in CN 0 of a predetermined value or more in step 603, the CN 0 of the periodic satellite at the reference point is acquired in step 604. The CN 0 at the reference point acquired here is the CN 0 at the same timing as the measurement point and the CN 0 one sidereal day before. These CN 0 data are transmitted from the communication device 63 of the reference station terminal 60 to the communication device 32 of the mobile station terminal 30. The satellite selection unit 1A calculates the difference between the CN 0 at the same timing at the received reference point and the CN 0 one sidereal day before. As another example, the difference between the CN 0 at the same timing at the reference point and the CN 0 one sidereal day before may be calculated in the reference station terminal 60, and the data of the difference in CN 0 may be transmitted to the mobile station terminal 30.

ステップ605において、測定点におけるCNの差と、基準点におけるCNの差とを比較し、同じか否かを判定する。「同じ」とは、完全に一致する場合に加え、許容可能な差はあっても実質的に同じである場合を含む。CNは、障害物等の周囲環境の変化だけでなく、雨や雪などの気象条件の変化によっても影響を受ける。この気象条件の変化は、多数の衛星のCNが影響を受けるため、測定点の周囲環境の変化による影響との区別がつきにくい。この場合、気象条件の変化によるCNの影響は、基準点においても同様に現れる(測定点と基準点は気象条件が共通することが前提)。したがって、ステップ605において、測定点におけるCNの差が基準点におけるCNの差と同じであれば、測定点におけるCNの差は気象条件に起因するもの(正常である)と判定し、初期化演算に進む。 In step 605, the difference in CN 0 at the measurement point is compared with the difference in CN 0 at the reference point to determine whether they are the same. "The same" includes cases where they are completely the same, as well as cases where there is a tolerable difference but they are substantially the same. CN 0 is affected not only by changes in the surrounding environment such as obstacles, but also by changes in meteorological conditions such as rain and snow. This change in meteorological conditions affects CN 0 of many satellites, so it is difficult to distinguish it from the effect of changes in the surrounding environment of the measurement point. In this case, the effect of CN 0 due to changes in meteorological conditions is also seen at the reference point (assuming that the meteorological conditions are the same as those of the measurement point and the reference point). Therefore, in step 605, if the difference in CN 0 at the measurement point is the same as the difference in CN 0 at the reference point, it is determined that the difference in CN 0 at the measurement point is due to meteorological conditions (is normal), and the initialization calculation is started.

ステップ605において、測定点におけるCNの差が基準点におけるCNの差と同じではない場合、気象条件以外の要因(障害物等による測定点の周囲環境の変化)によると判定し、その周期性衛星を、その後の演算に利用しない衛星として排除する。そして、残りの周期性衛星の電波を用いた初期化演算に進む。 In step 605, if the difference in CN0 at the measurement point is not the same as the difference in CN0 at the reference point, it is determined that the difference is due to factors other than meteorological conditions (changes in the surrounding environment of the measurement point due to obstacles, etc.), and the periodic satellite is excluded as a satellite that will not be used in subsequent calculations. Then, the initialization calculation is performed using the radio waves of the remaining periodic satellites.

上述した第1の実施形態では、測定点におけるCNの、1又は複数の恒星日の前後における差のみを評価基準としているので、排除する必要のない衛星までも排除してしまう場合がある。第2の実施形態では、測定点と基準点の双方におけるCNの、1又は複数の恒星日の前後における差を比較することによって、気象条件によるCNへの影響を、測定点における周囲環境の変化による影響から区別することができる。この結果、排除する必要のない衛星までも排除することを回避できる。 In the first embodiment described above, the evaluation criterion is only the difference in CN 0 at the measurement point before and after one or more sidereal days, so that satellites that do not need to be excluded may be excluded. In the second embodiment, by comparing the difference in CN 0 at both the measurement point and the reference point before and after one or more sidereal days, it is possible to distinguish the influence of meteorological conditions on CN 0 from the influence of changes in the surrounding environment at the measurement point. As a result, it is possible to avoid excluding satellites that do not need to be excluded.

その後の初期化演算部1Bによる初期化演算は、第1の実施形態における図2のステップ106、107と同様である。さらにその後の測位演算部1C及び変位判定部1Dよる処理は、第1の実施形態における図2のステップ108以降と同様である。 The subsequent initialization calculation by the initialization calculation unit 1B is similar to steps 106 and 107 in FIG. 2 in the first embodiment. Further, the subsequent processing by the positioning calculation unit 1C and the displacement determination unit 1D is similar to step 108 and subsequent steps in FIG. 2 in the first embodiment.

第2の実施形態の変形例として、上述した第1の実施形態の変形例として示した図4と同様に、衛星選択部1Aによる処理と、初期化演算部1Bによる処理の順序を逆にしてもよい。 As a modification of the second embodiment, the order of the processing by the satellite selection unit 1A and the processing by the initialization calculation unit 1B may be reversed, as in FIG. 4 shown as a modification of the first embodiment described above.

(3)第3の実施形態
図7は、本発明の第3の実施形態の測位装置を含む測位システム全体の一例を概略的に示す構成図である。測定点に設けられる移動局システム1は、アンテナ10、受信装置20、移動局端末30及び周囲環境監視装置40を有する。ここでも、移動局端末30が本発明の測位装置に相当する。基準局システム5は、第1の実施形態と同様である。
(3) Third embodiment Fig. 7 is a schematic diagram showing an example of an overall positioning system including a positioning device according to a third embodiment of the present invention. A mobile station system 1 provided at a measurement point has an antenna 10, a receiving device 20, a mobile station terminal 30, and a surrounding environment monitoring device 40. Again, the mobile station terminal 30 corresponds to the positioning device of the present invention. The reference station system 5 is the same as that of the first embodiment.

第3の実施形態では、移動局システム1が周囲環境監視装置40を有する点で第1の実施形態と異なる。この点に関連して第3の実施形態は、第1の実施形態とは中央処理装置31の衛星選択部1Aにおける処理が異なるが、これについては、図8を参照して測位方法と共に説明する。 The third embodiment differs from the first embodiment in that the mobile station system 1 has an ambient environment monitoring device 40. In relation to this, the third embodiment differs from the first embodiment in the processing in the satellite selection unit 1A of the central processing unit 31, which will be described together with the positioning method with reference to FIG. 8.

周囲環境監視装置40は、例えば、ビデオカメラ、物体を検出するセンサ等である。周囲環境監視装置40は、例えば、測定点の周囲の人や車などの障害物の有無を監視する。そして、周囲環境に変化があったとき、例えば新たな障害物の出現及び既存の障害物の消失を検出したとき、当該障害物の位置情報(アンテナ10から視た方位範囲、高さ等)とタイミングを含む周囲環境情報を生成する。 The surrounding environment monitoring device 40 is, for example, a video camera, a sensor that detects objects, etc. The surrounding environment monitoring device 40 monitors, for example, the presence or absence of obstacles such as people and cars around the measurement point. Then, when there is a change in the surrounding environment, for example when it detects the appearance of a new obstacle or the disappearance of an existing obstacle, it generates surrounding environment information including the position information of the obstacle (azimuth range, height, etc. as seen from the antenna 10) and timing.

周囲環境監視装置40は、移動局端末30と有線又は無線で通信可能に接続されている。周囲環境監視装置40は、周囲環境に変化があったとき周囲環境情報を移動局端末30に通信する。移動局端末30は、受信した周囲環境情報を記憶装置34に記憶する。一例として、移動局端末30が30分毎にウェークアップする場合、移動局端末30はウェークアップしたときに周囲環境監視装置40と通信を行って新たな周囲環境情報を取得し、記憶するようにしてもよい。 The surrounding environment monitoring device 40 is connected to the mobile station terminal 30 so as to be able to communicate with it via wired or wireless communication. The surrounding environment monitoring device 40 communicates surrounding environment information to the mobile station terminal 30 when there is a change in the surrounding environment. The mobile station terminal 30 stores the received surrounding environment information in the storage device 34. As an example, if the mobile station terminal 30 wakes up every 30 minutes, the mobile station terminal 30 may communicate with the surrounding environment monitoring device 40 when it wakes up to obtain and store new surrounding environment information.

図8は、第3の実施形態の測位方法の一例を概略的に示したフロー図である。図8では、主として第1の実施形態とは異なる部分を示している。ウェークアップは、図2のステップ101と同様である。 Figure 8 is a flow diagram that shows an example of a positioning method according to the third embodiment. Figure 8 mainly shows the parts that are different from the first embodiment. Wake-up is similar to step 101 in Figure 2.

ステップ801~807の処理は、衛星選択部1Aによって行われ、その後の初期化演算及び測位演算においてその電波を利用すべきでない衛星を排除する。 The processing of steps 801 to 807 is performed by the satellite selection unit 1A, which excludes satellites whose radio waves should not be used in subsequent initialization and positioning calculations.

ステップ801において、周囲環境情報を受信する。ステップ802において、受信した周囲環境情報から新たな障害物の出願が確認された場合、その障害物によって電波を正常に受信できないと予想される方位範囲を排除対象として決定する。決定された排除対象の方位範囲を記憶装置34に記憶する。そして、ステップ803において、排除対象の方位範囲に該当する衛星を、その後の演算に利用しない衛星として排除する。 In step 801, surrounding environment information is received. In step 802, if the presence of a new obstacle is confirmed from the received surrounding environment information, the direction range in which radio waves are expected to be unable to be received normally due to the obstacle is determined as a direction range to be excluded. The determined direction range to be excluded is stored in the storage device 34. Then, in step 803, satellites that fall within the direction range to be excluded are excluded as satellites that will not be used in subsequent calculations.

ステップ804~807は、第1の実施形態における図2のステップ102~105と同様である。別の例として、ステップ804~807に替えて、第2の実施形態における図6のステップ601~606を行ってもよい。その後の初期化演算部1Bによる初期化演算も、図2のステップ106、107と同様である。さらに、その後の測位演算部1C及び変位判定部1Dによる処理も、図2のステップ108以降と同様である。 Steps 804 to 807 are similar to steps 102 to 105 in FIG. 2 in the first embodiment. As another example, steps 601 to 606 in FIG. 6 in the second embodiment may be performed instead of steps 804 to 807. The subsequent initialization calculation by the initialization calculation unit 1B is also similar to steps 106 and 107 in FIG. 2. Furthermore, the subsequent processing by the positioning calculation unit 1C and the displacement determination unit 1D is also similar to step 108 and subsequent steps in FIG. 2.

第3の実施形態の変形例として、上述した第1の実施形態の変形例として示した図4と同様に、衛星選択部1Aによる処理と、初期化演算部1Bによる処理の順序を逆にしてもよい。 As a modification of the third embodiment, similar to FIG. 4 shown as a modification of the first embodiment described above, the order of the processing by the satellite selection unit 1A and the processing by the initialization calculation unit 1B may be reversed.

(4)第4の実施形態
図9は、本発明の第4の実施形態の測位装置を含む測位システム全体の一例を概略的に示す構成図である。測定点に設けられる移動局システム1は、アンテナ10、受信装置20及び移動局端末30を有する。一方、基準点に設けられる基準局システム5は、アンテナ51、受信装置52及び基準局端末60を有する。第4の実施形態では、移動局端末30及び基準局端末60は少なくとも、通信装置32、63をそれぞれ備え、衛星から受信した電波を適宜の通信ネットワークを介して遠隔演算装置70にそれぞれ送信する。第4の実施形態では、遠隔演算装置70が本発明の測位装置に相当する。遠隔演算装置70は、例えば所定のアプリケーションをインストールされたコンピュータであり、任意の場所に設置することができる。
(4) Fourth embodiment Fig. 9 is a schematic diagram showing an example of the overall positioning system including a positioning device according to a fourth embodiment of the present invention. A mobile station system 1 provided at a measurement point has an antenna 10, a receiving device 20, and a mobile station terminal 30. On the other hand, a reference station system 5 provided at a reference point has an antenna 51, a receiving device 52, and a reference station terminal 60. In the fourth embodiment, the mobile station terminal 30 and the reference station terminal 60 each have at least a communication device 32, 63, and transmit radio waves received from a satellite to a remote computing device 70 via an appropriate communication network. In the fourth embodiment, the remote computing device 70 corresponds to the positioning device of the present invention. The remote computing device 70 is, for example, a computer with a predetermined application installed, and can be installed at any location.

遠隔演算装置70は、中央処理装置71を備える。中央処理装置71は、図示しないが、CPUとRAM等のメモリを有し、CPUが所定のアプリケーションをメモリに読み込み実行することによって、本発明の測位装置の各機能を実現する演算処理を行う。中央処理装置71は、主な機能処理部として、衛星選択部1A、初期化演算部1B、測位演算部1C及び変位判定部1Dを有する。中央処理装置71は、移動局端末30及び基準局端末30から受信した電波に基づいて、上述した第1、第2又は第3の実施形態における移動局端末30の中央処理装置31と同様の処理(図2~図4、図6、図8参照)を実行することができる。すなわち、第4の実施形態では、衛星選択、初期化演算、測位演算及び変位判定の各処理が遠隔演算装置70で行われる。その場合、通信装置72、表示装置73、記憶装置74は、第1、第2又は第3の実施形態における移動局端末30の通信装置32、表示装置33、記憶装置34と同様に機能することができる。 The remote computing device 70 includes a central processing unit 71. Although not shown, the central processing unit 71 has a CPU and a memory such as a RAM, and the CPU loads and executes a predetermined application into the memory to perform computational processing to realize each function of the positioning device of the present invention. The central processing unit 71 has a satellite selection unit 1A, an initialization computation unit 1B, a positioning computation unit 1C, and a displacement determination unit 1D as main functional processing units. The central processing unit 71 can perform the same processing (see Figures 2 to 4, 6, and 8) as the central processing unit 31 of the mobile station terminal 30 in the first, second, or third embodiment described above, based on the radio waves received from the mobile station terminal 30 and the reference station terminal 30. That is, in the fourth embodiment, each of the processes of satellite selection, initialization computation, positioning computation, and displacement determination is performed by the remote computing device 70. In that case, the communication device 72, the display device 73, and the storage device 74 can function in the same way as the communication device 32, the display device 33, and the storage device 34 of the mobile station terminal 30 in the first, second, or third embodiment.

なお、中央処理装置71が、第2の実施形態の中央処理装置31(図5、図6参照)と同様の処理を行う場合、基準点における各衛星の電波のCNの検出は、中央処理装置71が基準局端末60から受信した電波に基づいて行ってもよく、又は、第2の実施形態と同様に基準局端末60において基準点における電波のCNの検出を行ってその結果を遠隔演算装置70に送信してもよい。 When the central processing unit 71 performs processing similar to that of the central processing unit 31 of the second embodiment (see Figures 5 and 6), the detection of CN0 of the radio waves of each satellite at the reference point may be performed by the central processing unit 71 based on the radio waves received from the reference station terminal 60, or, as in the second embodiment, the detection of CN0 of the radio waves at the reference point may be performed in the reference station terminal 60 and the result transmitted to the remote computing device 70.

第4の実施形態の変形例として、遠隔演算装置70が、衛星選択、初期化演算、測位演算及び変位判定の各処理の全てを実行するのではなく、一部のみを実行することも可能である。例えば、衛星選択の処理は移動局端末30で行い、移動局端末30は選択された衛星の電波のみを遠隔演算装置70に送信し、遠隔演算装置70が初期化演算以降の処理を実行するように構成してもよい。 As a variation of the fourth embodiment, the remote computing device 70 may execute only some of the processes of satellite selection, initialization calculation, positioning calculation, and displacement determination, rather than all of them. For example, the satellite selection process may be performed by the mobile station terminal 30, which may transmit only radio waves from the selected satellite to the remote computing device 70, and the remote computing device 70 may execute the processes from the initialization calculation onwards.

以上、本発明の各実施形態を、一つの構成例を参照して説明したが、本発明の実施形態は、上述した例に限られず、本発明の主旨から逸脱しない限り多様な変形形態が可能であり、それらの変形形態も本発明の範囲に含まれる。 Although each embodiment of the present invention has been described above with reference to one configuration example, the embodiments of the present invention are not limited to the above-mentioned examples, and various modifications are possible without departing from the spirit of the present invention, and such modifications are also included in the scope of the present invention.

1 移動局システム
10 アンテナ
20 受信装置
30 移動局端末
31 中央処理装置
32 通信装置
33 表示装置
34 記憶装置
40 周囲環境監視装置
5 基準局システム
51 アンテナ
52 受信装置
53 通信装置
60 基準局端末
61 中央処理装置
62 記憶装置
63 通信装置
REFERENCE SIGNS LIST 1 Mobile station system 10 Antenna 20 Receiving device 30 Mobile station terminal 31 Central processing unit 32 Communication device 33 Display device 34 Storage device 40 Surrounding environment monitoring device 5 Reference station system 51 Antenna 52 Receiving device 53 Communication device 60 Reference station terminal 61 Central processing unit 62 Storage device 63 Communication device

Claims (10)

測定点及び基準点において複数の衛星からそれぞれ受信する電波を用いて初期化演算ステップと測位演算ステップとを有するRTK法により測定点の測位を行う測位方法において、
前記初期化演算ステップの前に、又は、前記初期化演算ステップと前記測位演算ステップとの間に、その後の演算に利用する衛星を選択するために実行される衛星選択ステップを有し、前記衛星選択ステップが、
恒星日毎に周期性のある周期性衛星以外の衛星を、その後の演算に利用しない衛星として排除するステップと、
各前記周期性衛星について、測定点における電波の搬送波電力対雑音電力密度比であるCNを検出し記憶するステップと、
各前記周期性衛星について、検出されたCNを、1又は複数の恒星日前の記憶されたCNと比較し、所定値以上の差があるか否かを判定するステップと、
CNに所定値以上の差があると判定された周期性衛星を、その後の演算に利用しない衛星として排除するステップと、を有することを特徴とする衛星を用いた測位方法。
A positioning method for measuring the position of a measurement point by an RTK method, the method having an initialization calculation step and a positioning calculation step, using radio waves received from a plurality of satellites at a measurement point and a reference point,
A satellite selection step is executed before the initialization calculation step or between the initialization calculation step and the positioning calculation step to select a satellite to be used in a subsequent calculation, and the satellite selection step includes:
a step of excluding satellites other than the periodic satellites having a periodicity every sidereal day as satellites not to be used in subsequent calculations;
Detecting and storing a carrier power to noise power density ratio (CN0 ) of radio waves at a measurement point for each of the periodic satellites;
for each said periodic satellite, comparing the detected CN0 with a stored CN0 from one or more sidereal days ago to determine whether there is a difference of at least a predetermined value;
and excluding a periodic satellite determined to have a difference in CN0 of a predetermined value or more as a satellite not to be used in subsequent calculations.
測定点及び基準点において複数の衛星からそれぞれ受信する電波を用いて初期化演算ステップと測位演算ステップとを有するRTK法により測定点の測位を行う測位方法において、
前記初期化演算ステップの前に、又は、前記初期化演算ステップと前記測位演算ステップとの間に、その後の演算に利用する衛星を選択するために実行される衛星選択ステップを有し、前記衛星選択ステップが、
恒星日毎に周期性のある周期性衛星以外の衛星を、その後の演算に利用しない衛星として排除するステップと、
各前記周期性衛星について、測定点における電波の搬送波電力対雑音電力密度比であるCNを検出し記憶するステップと、
各前記周期性衛星について、検出されたCNを、1又は複数の恒星日前の記憶されたCNと比較し、所定値以上の差があるか否かを判定するステップと、
CNに所定値以上の差があると判定された周期性衛星について、基準点における同一タイミングのCNと、測定点と同じ1又は複数の恒星日前の基準点におけるCNとの差を取得するステップと、
CNに所定値以上の差があると判定された周期性衛星について、当該CNの差が基準点におけるCNの差と同じであるか否かを判定するステップと、
CNに所定値以上の差があると判定された周期性衛星のうち、当該CNの差が基準点におけるCNの差と同じではないと判定された周期性衛星を、その後の演算に利用しない衛星として排除するステップと、を有することを特徴とする衛星を用いた測位方法。
A positioning method for measuring the position of a measurement point by an RTK method, the method having an initialization calculation step and a positioning calculation step, using radio waves received from a plurality of satellites at a measurement point and a reference point,
A satellite selection step is executed before the initialization calculation step or between the initialization calculation step and the positioning calculation step to select a satellite to be used in a subsequent calculation, and the satellite selection step includes:
a step of excluding satellites other than the periodic satellites having a periodicity every sidereal day as satellites not to be used in subsequent calculations;
Detecting and storing a carrier power to noise power density ratio (CN0 ) of radio waves at a measurement point for each of the periodic satellites;
for each said periodic satellite, comparing the detected CN0 with a stored CN0 from one or more sidereal days ago to determine whether there is a difference of at least a predetermined value;
For a periodic satellite determined to have a difference of a predetermined value or more in CN 0 , a step of acquiring a difference between CN 0 at the same timing at a reference point and CN 0 at a reference point one or more sidereal days ago that is the same as the measurement point;
For a periodic satellite determined to have a difference in CN0 of a predetermined value or more, determining whether or not the difference in CN0 is the same as the difference in CN0 at a reference point;
and excluding, from among the periodic satellites determined to have a difference in CN0 of a predetermined value or more, a periodic satellite whose difference in CN0 is determined to be not the same as the difference in CN0 at a reference point, as a satellite not to be used in subsequent calculations.
前記衛星選択ステップが、
測定点の周囲環境情報に基づいて、測定点における排除対象の方位範囲を決定するステップと、
決定された前記方位範囲に該当する衛星を、その後の演算に利用しない衛星として排除するステップと、をさらに有することを特徴とする請求項1又は2に記載の衛星を用いた測位方法。
The satellite selection step includes:
determining a range of directions to be excluded at the measurement point based on ambient environment information of the measurement point;
3. The method of claim 1, further comprising the step of excluding satellites that fall within the determined azimuth range as satellites that will not be used in subsequent calculations.
前記初期化演算ステップと前記測位演算ステップとを、一定時間の間隔で行うことを特徴とする請求項1~3のいずれかに記載の衛星を用いた測位方法。 The satellite-based positioning method according to any one of claims 1 to 3, characterized in that the initialization calculation step and the positioning calculation step are performed at regular time intervals. 前記測位演算ステップが、
利用する周期性衛星のうち5個以上の周期性衛星を含む全ての組合せについて、所定の時間範囲内に所定の時間間隔毎の電波を用いて測位演算を行うステップと、
各前記組合せにおける前記時間間隔毎の測位演算結果について1恒星日前との差分である変位量を計算し、前記時間範囲内の変位量の平均を計算するステップと、
変位量が閾値以上の組合せが過半数以上であるか否かを判定するステップと、
変位量が閾値以上の組合せが過半数以上でない場合、変位なしと判定するステップと、を有することを特徴とする
請求項1~4のいずれかに記載の衛星を用いた測位方法。
The positioning calculation step includes:
performing positioning calculations using radio waves at predetermined time intervals within a predetermined time range for all combinations including five or more periodic satellites among the periodic satellites to be used;
calculating a displacement amount, which is a difference from one sidereal day ago, for the positioning calculation result for each of the combinations for each time interval, and calculating an average of the displacement amounts within the time range;
determining whether or not the majority of combinations have a displacement amount equal to or greater than a threshold;
5. The method of claim 1, further comprising a step of determining that there is no displacement when the majority of the combinations do not have a displacement amount equal to or greater than the threshold value.
変位量が閾値以上の組合せが過半数以上である場合、特定の周期性衛星を含む組合せのみが閾値を超えているか否かをさらに判定するステップと、
特定の周期性衛星を含む組合せのみが閾値を超えているのではない場合、変位ありと判定するステップと、を有することを特徴とする
請求項5に記載の衛星を用いた測位方法。
If the majority or more of the combinations have a displacement amount equal to or greater than the threshold, a further step of determining whether or not only a combination including a specific periodic satellite exceeds the threshold;
The method for positioning using satellites according to claim 5, further comprising a step of determining that there is displacement when only a combination including a specific periodic satellite does not exceed the threshold value.
測定点及び基準点において複数の衛星からそれぞれ受信する電波を用いて、RTK法により測定点の測位を行うべく初期化演算を行う初期化演算部と測位演算を行う測位演算部とを有する測位装置において、
初期化演算の前に、又は、初期化演算と測位演算との間に、その後の演算に利用する衛星の選択を行う衛星選択部を有し、前記衛星選択部が、
恒星日毎に周期性のある周期性衛星以外の衛星を、その後の演算に利用しない衛星として排除する手段と、
各前記周期性衛星について、測定点における電波の搬送波電力対雑音電力密度比であるCNを検出して記憶する手段と、
各前記周期性衛星について、検出されたCNを、1又は複数の恒星日前の記憶されたCNと比較し、所定値以上の差があるか否かを判定する手段と、
CNに所定値以上の差があると判定された周期性衛星を、その後の演算に利用しない衛星として排除する手段と、を有することを特徴とする衛星を用いた測位装置。
A positioning device having an initialization calculation unit that performs initialization calculations to perform positioning of a measurement point by an RTK method using radio waves received from a plurality of satellites at a measurement point and a reference point, and a positioning calculation unit that performs positioning calculations,
A satellite selection unit is provided for selecting a satellite to be used in a subsequent calculation before an initialization calculation or between the initialization calculation and the positioning calculation, the satellite selection unit being:
a means for excluding satellites other than a periodic satellite having a periodicity every sidereal day as satellites not to be used in subsequent calculations;
A means for detecting and storing a carrier power to noise power density ratio (CN0 ) of radio waves at a measurement point for each of the periodic satellites;
means for comparing, for each of said periodic satellites, the detected CN0 with a stored CN0 from one or more sidereal days ago to determine whether there is a difference between the detected CN0 and the stored CN0 by more than a predetermined value;
means for excluding a periodic satellite determined to have a difference in CN0 of a predetermined value or more as a satellite not to be used in subsequent calculations.
測定点及び基準点において複数の衛星からそれぞれ受信する電波を用いて、RTK法により測定点の測位を行うべく初期化演算を行う初期化演算部と測位演算を行う測位演算部とを有する測位装置において、
前記初期化演算の前に、又は、前記初期化演算と前記測位演算との間に、その後の演算に利用する衛星を選択するために実行される衛星選択部を有し、前記衛星選択部が、
恒星日毎に周期性のある周期性衛星以外の衛星を、その後の演算に利用しない衛星として排除する手段と、
各前記周期性衛星について、測定点における電波の搬送波電力対雑音電力密度比であるCNを検出して記憶する手段と、
各前記周期性衛星について、検出されたCNを、1又は複数の恒星日前の記憶されたCNと比較し、所定値以上の差があるか否かを判定する手段と、
CNに所定値以上の差があると判定された周期性衛星について、基準点における同一タイミングのCNと、測定点と同じ1又は複数の恒星日前の基準点におけるCNとの差を取得する手段と、
CNに所定値以上の差があると判定された周期性衛星について、当該CNの差が基準点におけるCNの差と同じであるか否かを判定する手段と、
CNに所定値以上の差があると判定された周期性衛星のうち、当該CNの差が基準点におけるCNの差と同じではないと判定された周期性衛星を、その後の演算に利用しない衛星として排除する手段と、を有することを特徴とする衛星を用いた測位装置。
A positioning device having an initialization calculation unit that performs initialization calculations to perform positioning of a measurement point by an RTK method using radio waves received from a plurality of satellites at a measurement point and a reference point, and a positioning calculation unit that performs positioning calculations,
A satellite selection unit is executed before the initialization operation or between the initialization operation and the positioning operation to select a satellite to be used in a subsequent operation, the satellite selection unit being:
a means for excluding satellites other than a periodic satellite having a periodicity every sidereal day as satellites not to be used in subsequent calculations;
A means for detecting and storing a carrier power to noise power density ratio (CN0 ) of radio waves at a measurement point for each of the periodic satellites;
means for comparing, for each of said periodic satellites, the detected CN0 with a stored CN0 from one or more sidereal days ago to determine whether there is a difference between the detected CN0 and the stored CN0 by more than a predetermined value;
For a periodic satellite determined to have a difference of a predetermined value or more in CN 0 , a means for acquiring a difference between CN 0 at the same timing at a reference point and CN 0 at a reference point one or more sidereal days ago that is the same as the measurement point;
A means for determining whether or not the difference in CN0 of a periodic satellite determined to have a difference of a predetermined value or more is the same as the difference in CN0 at a reference point;
and a means for excluding, from among periodic satellites determined to have a difference in CN0 of a predetermined value or more, a periodic satellite whose difference in CN0 is determined to be not the same as the difference in CN0 at a reference point, as a satellite not to be used in subsequent calculations.
前記衛星選択部が、
測定点の周囲環境情報に基づいて、測定点における排除対象の方位範囲を決定する手段と、
決定された前記方位範囲に含まれる衛星を、その後の演算に利用しない衛星として排除する手段と、をさらに有することを特徴とする請求項7又は8に記載の衛星を用いた測位装置。
The satellite selection unit:
A means for determining a range of directions to be excluded at the measurement point based on ambient environment information of the measurement point;
9. The satellite-based positioning device according to claim 7, further comprising: means for excluding satellites included in the determined azimuth range as satellites not to be used in subsequent calculations.
前記初期化演算及び前記測位演算を行うために、一定時間の間隔でウェークアップする手段と、
ウェークアップ後に前記初期化演算及び前記測位演算を含む一連の処理を実行し、完了後にスリープ状態とする手段と、をさらに有することを特徴とする請求項7~9のいずれかに記載の衛星を用いた測位装置。
a means for waking up at a constant time interval in order to perform the initialization calculation and the positioning calculation;
The satellite-based positioning device according to any one of claims 7 to 9, further comprising: a means for executing a series of processes including the initialization calculation and the positioning calculation after waking up, and going into a sleep state after completion of the processes.
JP2020068055A 2020-04-06 2020-04-06 Satellite-based positioning method and positioning device Active JP7570623B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2020068055A JP7570623B2 (en) 2020-04-06 2020-04-06 Satellite-based positioning method and positioning device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020068055A JP7570623B2 (en) 2020-04-06 2020-04-06 Satellite-based positioning method and positioning device

Publications (2)

Publication Number Publication Date
JP2021165640A JP2021165640A (en) 2021-10-14
JP7570623B2 true JP7570623B2 (en) 2024-10-22

Family

ID=78021933

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020068055A Active JP7570623B2 (en) 2020-04-06 2020-04-06 Satellite-based positioning method and positioning device

Country Status (1)

Country Link
JP (1) JP7570623B2 (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006132003A1 (en) 2005-06-06 2006-12-14 National University Corporation Tokyo University Of Marine Science And Technology Gps reception device and gps positioning correction method
JP2012058185A (en) 2010-09-13 2012-03-22 Nec Corp Satellite navigation reinforcement system using only high-reliability gps positioning signal
WO2012176723A1 (en) 2011-06-20 2012-12-27 古野電気株式会社 Gnss analysis device, gnss analysis system, gnss analysis program, and gnss analysis method
JP2014215167A (en) 2013-04-25 2014-11-17 古野電気株式会社 Gnss state display device and gnss state display method
JP2017182692A (en) 2016-03-31 2017-10-05 セコム株式会社 Autonomous mobile robot
US20180302148A1 (en) 2017-04-18 2018-10-18 Airbus Defence and Space GmbH Near-real time EIRP monitoring device/system
JP2018179734A (en) 2017-04-12 2018-11-15 清水建設株式会社 Displacement measurement method and displacement measurement system
WO2019171633A1 (en) 2018-03-09 2019-09-12 Necソリューションイノベータ株式会社 Moving body positioning system, method, and program
WO2019208592A1 (en) 2018-04-27 2019-10-31 国立大学法人静岡大学 Satellite positioning system

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006132003A1 (en) 2005-06-06 2006-12-14 National University Corporation Tokyo University Of Marine Science And Technology Gps reception device and gps positioning correction method
JP2012058185A (en) 2010-09-13 2012-03-22 Nec Corp Satellite navigation reinforcement system using only high-reliability gps positioning signal
WO2012176723A1 (en) 2011-06-20 2012-12-27 古野電気株式会社 Gnss analysis device, gnss analysis system, gnss analysis program, and gnss analysis method
JP2014215167A (en) 2013-04-25 2014-11-17 古野電気株式会社 Gnss state display device and gnss state display method
JP2017182692A (en) 2016-03-31 2017-10-05 セコム株式会社 Autonomous mobile robot
JP2018179734A (en) 2017-04-12 2018-11-15 清水建設株式会社 Displacement measurement method and displacement measurement system
US20180302148A1 (en) 2017-04-18 2018-10-18 Airbus Defence and Space GmbH Near-real time EIRP monitoring device/system
WO2019171633A1 (en) 2018-03-09 2019-09-12 Necソリューションイノベータ株式会社 Moving body positioning system, method, and program
WO2019208592A1 (en) 2018-04-27 2019-10-31 国立大学法人静岡大学 Satellite positioning system

Also Published As

Publication number Publication date
JP2021165640A (en) 2021-10-14

Similar Documents

Publication Publication Date Title
US7561101B1 (en) Last known position reporting for always-on global positioning system receiver
US9510318B2 (en) Method and apparatus for ascertaining a location of a personal portable wireless communication device
JP5957025B2 (en) User receiver positioning method
US9482739B2 (en) Indoor/outdoor transition determination
JP5641279B2 (en) GPS receiver and method of processing GPS signals
US8775076B2 (en) Position calculating method and position calculating device
US6570533B2 (en) Method for determining the phase of information, and an electronic device
JP2002320254A (en) Mobile communication apparatus and its position detection method
US20120223860A1 (en) Use of Motion or Accelerometer Sensors in Low Power Positioning System
US8462046B2 (en) Managing the measurement of signals
RU2625819C1 (en) Method of auxiliary holding in the joint navigation system radio set
JP2006215027A (en) Navigation satellite receiver system and method for accelerating warm start of high sensitivity navigation satellite receiver in weak signal environment
US20100207810A1 (en) Position calculating method and position calculating device
US11041962B2 (en) Adaptive multi-rate navigation
US20130257651A1 (en) Navigation bit boundary determination apparatus and a method therefor
US6597985B2 (en) Positioning apparatus and navigation apparatus, and positioning method and navigation method
US20120218146A1 (en) Dynamic Sleep Time Calculation for GNSS Receiver
JP4985310B2 (en) POSITIONING METHOD, PROGRAM, POSITIONING DEVICE, AND ELECTRONIC DEVICE
JP2010112834A (en) Displacement measurement system, satellite signal receiving device, and satellite positioning method
US20090156205A1 (en) Method and Apparatus for Establishing a Wireless Network Signal Acquisition Rate
JP2012503764A (en) How to optimize the acquisition of spread spectrum signals from satellites by mobile receivers
CN111781619A (en) Positioning method, device, equipment and storage medium based on near field communication network
JPWO2019208592A1 (en) Satellite positioning system
CA3054590A1 (en) Automatic pressure sensor output calibration for reliable altitude determination
Raghuvanshi et al. Improving smartphone positioning accuracy by adapting measurement covariance with t-test on innovations

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220603

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20220603

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20221207

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20230920

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20231031

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20240229

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240508

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20240627

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20241001

R150 Certificate of patent or registration of utility model

Ref document number: 7570623

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150