JP3532267B2 - 測位システム - Google Patents
測位システムInfo
- Publication number
- JP3532267B2 JP3532267B2 JP28514794A JP28514794A JP3532267B2 JP 3532267 B2 JP3532267 B2 JP 3532267B2 JP 28514794 A JP28514794 A JP 28514794A JP 28514794 A JP28514794 A JP 28514794A JP 3532267 B2 JP3532267 B2 JP 3532267B2
- Authority
- JP
- Japan
- Prior art keywords
- phase
- center
- observation
- value
- vector
- 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.)
- Expired - Fee Related
Links
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Description
【０００１】
【産業上の利用分野】この発明は、複数地点間の基線ベ
クトルを高精度でリアルタイムに求めることができる測
位システムに関する。
クトルを高精度でリアルタイムに求めることができる測
位システムに関する。
【０００２】
【従来の技術】従来オンラインで遠隔地点の位置をリア
ルタイムで検出するシステムとしては、特開平６−２３
０１０１号公報に示すシステムが提案されている。この
システムは、座標が既知で固定している基準地点および
移動を観測する観測地点にそれぞれ単独測位を行うＧＰ
Ｓ測位装置を設置して互いに無線で接続し、基準地点に
おけるＧＰＳ測位結果とその地点の座標との誤差を補正
ベクトルとし、観測地点から受信した測定結果をこの補
正ベクトルで補正して観測地点の位置を割り出すもので
ある。
ルタイムで検出するシステムとしては、特開平６−２３
０１０１号公報に示すシステムが提案されている。この
システムは、座標が既知で固定している基準地点および
移動を観測する観測地点にそれぞれ単独測位を行うＧＰ
Ｓ測位装置を設置して互いに無線で接続し、基準地点に
おけるＧＰＳ測位結果とその地点の座標との誤差を補正
ベクトルとし、観測地点から受信した測定結果をこの補
正ベクトルで補正して観測地点の位置を割り出すもので
ある。
【０００３】
【発明が解決しようとする課題】しかし、この方式で
は、単独測位方式が基本となっており、基準地点におけ
る単独測位の誤差を用いて観測局における単独測位の誤
差を補正しようとするものであるため、測位精度が１〜
数メートル程度あり、たとえば、地盤の変動や地滑りな
どにおける初期的な微小な変化を速やかに精度よく察知
することが困難である欠点があった。
は、単独測位方式が基本となっており、基準地点におけ
る単独測位の誤差を用いて観測局における単独測位の誤
差を補正しようとするものであるため、測位精度が１〜
数メートル程度あり、たとえば、地盤の変動や地滑りな
どにおける初期的な微小な変化を速やかに精度よく察知
することが困難である欠点があった。
【０００４】この発明は、搬送波位相の積算値を用いた
相対測位を用いることにより、観測地点の位置を高精度
でリアルタイムに求めることができる測位システムを提
供することを目的とする。
相対測位を用いることにより、観測地点の位置を高精度
でリアルタイムに求めることができる測位システムを提
供することを目的とする。
【０００５】
【課題を解決するための手段】この出願の請求項１の発
明は、複数の地点に設けられた観測装置と、各観測装置
とオンラインで結ばれたセンタとからなる測位システム
であって、前記観測装置は、複数のＧＰＳ衛星が送信す
る信号電波を受信してその信号または搬送波の位相の積
算値をカウントする手段と、カウントした位相の積算値
を一定時間毎に出力する手段とを備え、前記センタは、
各観測装置から受信した位相の積算値に基づいて上記複
数の地点間のベクトルを算出するベクトル算出手段と、
該基線ベクトルに基づき、観測装置の最も確からしい位
置を示す最確値および所定確率以上で観測装置の位置が
真となる範囲を示す信頼域を算出し、算出された信頼域
内に観測装置の初期の最確値が存在するかどうかを判定
する判定手段と、を備えたことを特徴とする。
明は、複数の地点に設けられた観測装置と、各観測装置
とオンラインで結ばれたセンタとからなる測位システム
であって、前記観測装置は、複数のＧＰＳ衛星が送信す
る信号電波を受信してその信号または搬送波の位相の積
算値をカウントする手段と、カウントした位相の積算値
を一定時間毎に出力する手段とを備え、前記センタは、
各観測装置から受信した位相の積算値に基づいて上記複
数の地点間のベクトルを算出するベクトル算出手段と、
該基線ベクトルに基づき、観測装置の最も確からしい位
置を示す最確値および所定確率以上で観測装置の位置が
真となる範囲を示す信頼域を算出し、算出された信頼域
内に観測装置の初期の最確値が存在するかどうかを判定
する判定手段と、を備えたことを特徴とする。
【０００６】この出願の請求項２の発明は、前記センタ
に、刻々入力される位相の積算値に基づいて前記ベクト
ル算出手段を繰り返して実行し、判定手段で順次算出さ
れる信頼域内に前記観測装置の初期の最確値が存在しな
いとき警報を発する手段を備えたことを特徴とする。
に、刻々入力される位相の積算値に基づいて前記ベクト
ル算出手段を繰り返して実行し、判定手段で順次算出さ
れる信頼域内に前記観測装置の初期の最確値が存在しな
いとき警報を発する手段を備えたことを特徴とする。
【０００７】
【作用】この発明の測定システムでは、複数地点に設け
られた観測装置が、複数のＧＰＳ衛星が送信する信号電
波を受信してその信号または搬送波の位相の積算値をカ
ウントし、これを一定時間毎に出力する。センタは、複
数の観測装置からそれぞれ複数のＧＰＳ衛星の位相の積
算値を受信する。そしてこれに基づいて前記複数の地点
間の基線ベクトルを算出する。算出の方法は、従来より
周知の２重位相差法を用いればよい。これにより、極め
て高い精度でオンラインでデータを受信してほぼリアル
タイムに複数地点間の基線ベクトルを求めることができ
る。
られた観測装置が、複数のＧＰＳ衛星が送信する信号電
波を受信してその信号または搬送波の位相の積算値をカ
ウントし、これを一定時間毎に出力する。センタは、複
数の観測装置からそれぞれ複数のＧＰＳ衛星の位相の積
算値を受信する。そしてこれに基づいて前記複数の地点
間の基線ベクトルを算出する。算出の方法は、従来より
周知の２重位相差法を用いればよい。これにより、極め
て高い精度でオンラインでデータを受信してほぼリアル
タイムに複数地点間の基線ベクトルを求めることができ
る。
【０００８】また、この発明では、この基線ベクトルを
継続して繰り返し算出し、このベクトルが一定以上変動
した場合には、警報を発するようにした。基線ベクトル
の変動は、観測地点を山岳部に設定した場合には崖崩れ
や雪崩などに対応するものであり、これにより、精度の
高い速やかな地滑り，崖崩れ，雪崩などの観測が可能に
なる。
継続して繰り返し算出し、このベクトルが一定以上変動
した場合には、警報を発するようにした。基線ベクトル
の変動は、観測地点を山岳部に設定した場合には崖崩れ
や雪崩などに対応するものであり、これにより、精度の
高い速やかな地滑り，崖崩れ，雪崩などの観測が可能に
なる。
【０００９】
【実施例】図１はこの発明の実施例である測位システム
の構成を示す図である。この測位システムは、基準地点
に設置される基準局３と観測地点に設置される観測局２
とを有している。基準地点は、正確な位置（緯度，経
度，高度）が分かっており、移動するおそれのない場所
に設定される。この測位システムは、たとえば、堤防や
ダムの決壊の検知，地滑りの検知，崖崩れの検知，埋め
立て地の不等沈下の計測・監視、滑走路・道路の凹凸の
計測・監視などに適用されるが、この場合、基準地点は
山麓や固い地盤の地点に設定され、その地点に基準局３
が設けられる。また、観測地点は決壊や崖崩れを起こし
そうな地点や、不等沈下，凹凸を起こしそうな地点に設
置され、その地点に観測局２が設けられる。
の構成を示す図である。この測位システムは、基準地点
に設置される基準局３と観測地点に設置される観測局２
とを有している。基準地点は、正確な位置（緯度，経
度，高度）が分かっており、移動するおそれのない場所
に設定される。この測位システムは、たとえば、堤防や
ダムの決壊の検知，地滑りの検知，崖崩れの検知，埋め
立て地の不等沈下の計測・監視、滑走路・道路の凹凸の
計測・監視などに適用されるが、この場合、基準地点は
山麓や固い地盤の地点に設定され、その地点に基準局３
が設けられる。また、観測地点は決壊や崖崩れを起こし
そうな地点や、不等沈下，凹凸を起こしそうな地点に設
置され、その地点に観測局２が設けられる。
【００１０】各基準局３，観測局２にはＧＰＳアンテナ
およびＧＰＳ信号の搬送波位相の積算値を出力するＧＰ
Ｓ受信機が設置されている。ＧＰＳ受信機は、複数のＧ
ＰＳ衛星１から送られてくるＧＰＳ信号の２つの搬送波
（Ｌ１波およびＬ２波）の位相積算値を一定時間毎に出
力する装置である。なお、Ｌ１波の搬送波周波数は１５
７５．４２ＭＨｚ、Ｌ２波の搬送波周波数は１２２７．
６ＭＨｚである。観測局２が出力した位相積算値は通信
ケーブル４−中継器５を介して通信制御装置６に送信さ
れる。また、基準局３が出力した位相積算値も通信ケー
ブル４を介して通信制御装置６に入力される。通信制御
装置６は、データ処理装置であるセンタ８に接続されて
おり、通信制御装置６が各観測局２，基準局３から受信
した位相積算値はセンタ８に集約される。センタ８は、
集約した位相積算値に基づいて各観測局２（観測地点）
の位置を割り出す。この位置は基準局３（基準地点）か
らの基線ベクトルとして求められる。基線ベクトルの算
出手法としては、最小２乗法を用いた網平均手法などが
用いられる。センタ８は、一定時間毎に最新の位相積算
値を用いて基線ベクトルの算出を実行しており、新たに
算出された基線ベクトルを過去に求めた基線ベクトルと
比較することにより、観測地点の移動、すなわち、堤防
やダムの決壊，崖崩れ，埋め立て地の不等沈下，滑走路
・道路の凹凸などが判断される。
およびＧＰＳ信号の搬送波位相の積算値を出力するＧＰ
Ｓ受信機が設置されている。ＧＰＳ受信機は、複数のＧ
ＰＳ衛星１から送られてくるＧＰＳ信号の２つの搬送波
（Ｌ１波およびＬ２波）の位相積算値を一定時間毎に出
力する装置である。なお、Ｌ１波の搬送波周波数は１５
７５．４２ＭＨｚ、Ｌ２波の搬送波周波数は１２２７．
６ＭＨｚである。観測局２が出力した位相積算値は通信
ケーブル４−中継器５を介して通信制御装置６に送信さ
れる。また、基準局３が出力した位相積算値も通信ケー
ブル４を介して通信制御装置６に入力される。通信制御
装置６は、データ処理装置であるセンタ８に接続されて
おり、通信制御装置６が各観測局２，基準局３から受信
した位相積算値はセンタ８に集約される。センタ８は、
集約した位相積算値に基づいて各観測局２（観測地点）
の位置を割り出す。この位置は基準局３（基準地点）か
らの基線ベクトルとして求められる。基線ベクトルの算
出手法としては、最小２乗法を用いた網平均手法などが
用いられる。センタ８は、一定時間毎に最新の位相積算
値を用いて基線ベクトルの算出を実行しており、新たに
算出された基線ベクトルを過去に求めた基線ベクトルと
比較することにより、観測地点の移動、すなわち、堤防
やダムの決壊，崖崩れ，埋め立て地の不等沈下，滑走路
・道路の凹凸などが判断される。
【００１１】図２に上記観測局２，基準局３に設置され
るＧＰＳ受信機の構成を示す。ＧＰＳ受信機は、複数個
（一般的には４個）のＧＰＳ衛星のＬ１波，Ｌ２波を受
信し、一定時間（例えば、１０秒）毎の度単位の波数を
位相角単位で積算してその積算値を出力する装置であ
る。この装置は、アンテナ１０，周波数変換部１１，ク
ロック抽出部１２，位相カウンタ１３，通信制御部１４
および時計１５を備えているが、少なくとも、周波数変
換部１１，クロック抽出部１２，位相カウンタ１３は、
４個のＧＰＳ衛星のＬ１波，Ｌ２波を同時受信できるよ
うに８個並列に設けられているものとする。
るＧＰＳ受信機の構成を示す。ＧＰＳ受信機は、複数個
（一般的には４個）のＧＰＳ衛星のＬ１波，Ｌ２波を受
信し、一定時間（例えば、１０秒）毎の度単位の波数を
位相角単位で積算してその積算値を出力する装置であ
る。この装置は、アンテナ１０，周波数変換部１１，ク
ロック抽出部１２，位相カウンタ１３，通信制御部１４
および時計１５を備えているが、少なくとも、周波数変
換部１１，クロック抽出部１２，位相カウンタ１３は、
４個のＧＰＳ衛星のＬ１波，Ｌ２波を同時受信できるよ
うに８個並列に設けられているものとする。
【００１２】ＧＰＳ衛星１から信号電波を受信するＧＰ
Ｓアンテナ１０は、無指向性のアンテナであり、どの方
位にあるＧＰＳ衛星の信号電波をも受信することができ
る。このＧＰＳアンテナ１０には周波数変換部１１が接
続されている。周波数変換部１１は受信したマイクロ波
の信号電波を回路内で処理可能な周波数の高周波信号に
変換する。この信号はクロック抽出部１２に入力され
る。クロック抽出部１２は、中間周波に変換された高周
波信号の搬送波クロックを抽出する回路である。クロッ
ク抽出部１２はＰＬＬ回路で構成されており、入力信号
と位相が同期した発振器信号を抽出搬送波クロック信号
として出力する。
Ｓアンテナ１０は、無指向性のアンテナであり、どの方
位にあるＧＰＳ衛星の信号電波をも受信することができ
る。このＧＰＳアンテナ１０には周波数変換部１１が接
続されている。周波数変換部１１は受信したマイクロ波
の信号電波を回路内で処理可能な周波数の高周波信号に
変換する。この信号はクロック抽出部１２に入力され
る。クロック抽出部１２は、中間周波に変換された高周
波信号の搬送波クロックを抽出する回路である。クロッ
ク抽出部１２はＰＬＬ回路で構成されており、入力信号
と位相が同期した発振器信号を抽出搬送波クロック信号
として出力する。
【００１３】図３にクロック抽出部１２の回路構成を示
す。クロック抽出部１２は、２乗回路２０，フィルタ２
１，位相比較回路２２，ループフィルタ２３およびＶＣ
Ｏ２４からなる位相同期ループ回路で構成されている。
この回路は、所定のＧＰＳ衛星の信号搬送波をそのＧＰ
Ｓ衛星の運行によるドップラ変位を手掛かりに抽出する
回路である。まず、搬送波を抽出したいＧＰＳ衛星の運
行データに基づいてその衛星の軌道を割り出し、その軌
道を移動することによる搬送波周波数のドップラ変位分
を算出する。この算出は図示しない制御部が行う。この
算出値が予測ドップラ値としてＶＣＯ２４に入力され
る。ＶＣＯ２４はこの予測ドップラ値によって規定され
る周波数の高周波信号を発振する。この信号は位相比較
回路２２に入力される。一方、中間周波に変換された受
信信号は２乗回路２０で２乗され２倍の周波数の信号に
変換される。この信号はフィルタ２１で搬送波周波数付
近のみを透過するフィルタ２１で他の周波数と分離さ
れ、搬送波周波数のみが位相比較回路２２に入力され
る。位相比較回路２２では、２乗回路２０で２乗された
入力信号とＶＣＯ２４が発振した信号の位相を比較す
る。位相のずれは電圧値として出力されループフィルタ
２３を介してＶＣＯ２４に入力される。ＶＣＯ２４はこ
の位相差電圧によって発振周波数を修正され、入力され
た搬送波周波数と同一周波数の信号を発振するようにな
る。この信号は位相比較回路２４にフィードバックされ
るとともに抽出搬送波クロックとして次段の位相カウン
タ１３に出力される。
す。クロック抽出部１２は、２乗回路２０，フィルタ２
１，位相比較回路２２，ループフィルタ２３およびＶＣ
Ｏ２４からなる位相同期ループ回路で構成されている。
この回路は、所定のＧＰＳ衛星の信号搬送波をそのＧＰ
Ｓ衛星の運行によるドップラ変位を手掛かりに抽出する
回路である。まず、搬送波を抽出したいＧＰＳ衛星の運
行データに基づいてその衛星の軌道を割り出し、その軌
道を移動することによる搬送波周波数のドップラ変位分
を算出する。この算出は図示しない制御部が行う。この
算出値が予測ドップラ値としてＶＣＯ２４に入力され
る。ＶＣＯ２４はこの予測ドップラ値によって規定され
る周波数の高周波信号を発振する。この信号は位相比較
回路２２に入力される。一方、中間周波に変換された受
信信号は２乗回路２０で２乗され２倍の周波数の信号に
変換される。この信号はフィルタ２１で搬送波周波数付
近のみを透過するフィルタ２１で他の周波数と分離さ
れ、搬送波周波数のみが位相比較回路２２に入力され
る。位相比較回路２２では、２乗回路２０で２乗された
入力信号とＶＣＯ２４が発振した信号の位相を比較す
る。位相のずれは電圧値として出力されループフィルタ
２３を介してＶＣＯ２４に入力される。ＶＣＯ２４はこ
の位相差電圧によって発振周波数を修正され、入力され
た搬送波周波数と同一周波数の信号を発振するようにな
る。この信号は位相比較回路２４にフィードバックされ
るとともに抽出搬送波クロックとして次段の位相カウン
タ１３に出力される。
【００１４】位相カウンタ１３は１０8 程度のカウンタ
であり、搬送波クロックの位相を度単位でカウントす
る。たとえば、最終の中間周波数を１０ｋＨｚに選んだ
とするとωｔに相当する位相回転は毎秒３．６×１０6
度となる。これに種々の変化要素が加味されたとして
も、１０秒間カウントすると１．７〜５．５×１０7 程
度の大きさの数になるため、カウンタがオーバフローし
た場合でもキャリー信号を無視して０からのカウントを
繰り返すようにする。時計１５はこの位相カウンタ１３
に対して一定時間毎にカウント値を出力するようにトリ
ガを入力する。このカウント値は通信制御部１４に対し
て位相積算値として出力される。通信制御部１４は、図
４に示すように、発信元を示す局名，受信したＧＰＳ衛
星の番号およびＬ１波，Ｌ２波の別、観測時刻、および
その時刻におけるカウント数からなる電文を編集してセ
ンタ８の通信制御部６に対して送出する。
であり、搬送波クロックの位相を度単位でカウントす
る。たとえば、最終の中間周波数を１０ｋＨｚに選んだ
とするとωｔに相当する位相回転は毎秒３．６×１０6
度となる。これに種々の変化要素が加味されたとして
も、１０秒間カウントすると１．７〜５．５×１０7 程
度の大きさの数になるため、カウンタがオーバフローし
た場合でもキャリー信号を無視して０からのカウントを
繰り返すようにする。時計１５はこの位相カウンタ１３
に対して一定時間毎にカウント値を出力するようにトリ
ガを入力する。このカウント値は通信制御部１４に対し
て位相積算値として出力される。通信制御部１４は、図
４に示すように、発信元を示す局名，受信したＧＰＳ衛
星の番号およびＬ１波，Ｌ２波の別、観測時刻、および
その時刻におけるカウント数からなる電文を編集してセ
ンタ８の通信制御部６に対して送出する。
【００１５】なお、このＧＰＳ受信機では、Ｌ１，Ｌ２
の２周波を受信することによって電離層の影響による誤
差を除去し、複数の衛星の信号電波を受信することによ
って周波数変換用の局部発振回路の周波数ゆらぎの影響
を除去している。なお、基準局，観測局の設置場所が数
十キロメートルの範囲に納まる場合であれば、全ての局
における電離層の影響は同一であると考えられるため一
方の周波（Ｌ１波）のみを受信しても基線ベクトル算出
時に電離層の影響は相殺される。
の２周波を受信することによって電離層の影響による誤
差を除去し、複数の衛星の信号電波を受信することによ
って周波数変換用の局部発振回路の周波数ゆらぎの影響
を除去している。なお、基準局，観測局の設置場所が数
十キロメートルの範囲に納まる場合であれば、全ての局
における電離層の影響は同一であると考えられるため一
方の周波（Ｌ１波）のみを受信しても基線ベクトル算出
時に電離層の影響は相殺される。
【００１６】センタ８では、各局から受信した積算値に
基づいて相対測位を実行する。相対測位の原理を図５を
参照しつつ説明する。同図において、α，βをＧＰＳ衛
星、Ａ，Ｂをアンテナとする。このうち、Ａが基準局の
アンテナ（観測点），Ｂが観測局のアンテナとする。ま
た、観測点Ａ，Ｂ間の位相差をφαとすると、 φα＝ｒ・Ｓα／λ となる。ここでｒはＡ，Ｂ間の基線ベクトル、Ｓαは観
測点から衛星αへの方向余弦ベクトル，λは搬送波の波
長である。
基づいて相対測位を実行する。相対測位の原理を図５を
参照しつつ説明する。同図において、α，βをＧＰＳ衛
星、Ａ，Ｂをアンテナとする。このうち、Ａが基準局の
アンテナ（観測点），Ｂが観測局のアンテナとする。ま
た、観測点Ａ，Ｂ間の位相差をφαとすると、 φα＝ｒ・Ｓα／λ となる。ここでｒはＡ，Ｂ間の基線ベクトル、Ｓαは観
測点から衛星αへの方向余弦ベクトル，λは搬送波の波
長である。
【００１７】衛星α，βの位相差の差を求めると、
φα−φβ＝ｒ・（Ｓα−Ｓβ）／λ
となる。φα−φβは観測によって求められ、Ｓα−Ｓ
βは軌道計算によって正確に算出することができること
から、いくつかのＧＰＳ衛星の組に対して観測を行い、
同様の式を多数作成することにより、これらを最小二乗
法的に解くことによって正確なｒを求めることができ
る。
βは軌道計算によって正確に算出することができること
から、いくつかのＧＰＳ衛星の組に対して観測を行い、
同様の式を多数作成することにより、これらを最小二乗
法的に解くことによって正確なｒを求めることができ
る。
【００１８】さらに、３以上の地点の位相積算値が入力
されるシステムの場合には、網平均の手法を用いて各地
点の位置（基線ベクトル）を同時に高精度で算出するこ
とができる。網平均は、３以上のＧＰＳ受信機のそれぞ
れを網状に組み合わせ、各ＧＰＳ受信機間の基線ベクト
ルの最確値と誤差をネットワーク全体として最も確から
しい値になるように調整する計算方法である。一般的
に、調整は行列を用いた最小二乗法によって全部の観測
値を同時に調整する。このように、複数の衛星を用いて
網平均によって各点の位置を割り出すようにしたことに
より、各観測点の位置決定精度は数ｍｍ〜数ｃｍに高め
ることができる。
されるシステムの場合には、網平均の手法を用いて各地
点の位置（基線ベクトル）を同時に高精度で算出するこ
とができる。網平均は、３以上のＧＰＳ受信機のそれぞ
れを網状に組み合わせ、各ＧＰＳ受信機間の基線ベクト
ルの最確値と誤差をネットワーク全体として最も確から
しい値になるように調整する計算方法である。一般的
に、調整は行列を用いた最小二乗法によって全部の観測
値を同時に調整する。このように、複数の衛星を用いて
網平均によって各点の位置を割り出すようにしたことに
より、各観測点の位置決定精度は数ｍｍ〜数ｃｍに高め
ることができる。
【００１９】センタ８では、このようにして求められた
観測局２の位置が過去に求められていた位置から大きく
移動していた場合には、たとえば崖崩れがあったと判断
して警報を発する。
観測局２の位置が過去に求められていた位置から大きく
移動していた場合には、たとえば崖崩れがあったと判断
して警報を発する。
【００２０】ここで、図６を参照して位置確定の手法と
移動検出の手法を説明する。最小二乗法による基線ベク
トルの算出によれば、最確値，信頼域という２つの計算
結果を得ることができる。最確値は得られたデータから
最も確からしいと考えられる点である。信頼域は、ある
一定の確率（例えば９９％）で真の値が存在する範囲を
示している。まず、観測の最初に求めた最確値を初期位
置とし、その後周期的にアンテナの最確値と信頼域を求
める。そして、初期位置がそのときの信頼域の範囲内に
存在するかを判定する。その信頼域の範囲内に初期位置
があれば「移動なし」と判断し、無ければ「移動あり」
と判断する。
移動検出の手法を説明する。最小二乗法による基線ベク
トルの算出によれば、最確値，信頼域という２つの計算
結果を得ることができる。最確値は得られたデータから
最も確からしいと考えられる点である。信頼域は、ある
一定の確率（例えば９９％）で真の値が存在する範囲を
示している。まず、観測の最初に求めた最確値を初期位
置とし、その後周期的にアンテナの最確値と信頼域を求
める。そして、初期位置がそのときの信頼域の範囲内に
存在するかを判定する。その信頼域の範囲内に初期位置
があれば「移動なし」と判断し、無ければ「移動あり」
と判断する。
【００２１】図７は同センタ８の動作を示すフローチャ
ートである。
ートである。
【００２２】まず、ｎ１で基線ベクトルを算出する時刻
か否かを判断する。算出時刻であればｎ２以下の動作を
実行する。ｎ２では、観測局２，基準局３から入力され
た同時刻の搬送波位相の積算値を収集して（ｎ２）、２
つの局（ＧＰＳ受信機）の組が作る全ての基線ベクトル
を二重位相差と最小二乗法を用いて求める（ｎ３）。こ
のようにして求められた全ての基線ベクトルを網平均に
よりネットワーク全体として最も確からしい基線ベクト
ルに補正する（ｎ４）。この計算結果を過去の計算結果
を比較して各観測点の移動の有無を判断する（ｎ５）。
複数の観測点２のうちいずれかに移動が検出された場合
には（ｎ６）、警報を発し（ｎ７）、係員に移動を告知
する。移動がない場合には、そのままルーチン動作を継
続する。
か否かを判断する。算出時刻であればｎ２以下の動作を
実行する。ｎ２では、観測局２，基準局３から入力され
た同時刻の搬送波位相の積算値を収集して（ｎ２）、２
つの局（ＧＰＳ受信機）の組が作る全ての基線ベクトル
を二重位相差と最小二乗法を用いて求める（ｎ３）。こ
のようにして求められた全ての基線ベクトルを網平均に
よりネットワーク全体として最も確からしい基線ベクト
ルに補正する（ｎ４）。この計算結果を過去の計算結果
を比較して各観測点の移動の有無を判断する（ｎ５）。
複数の観測点２のうちいずれかに移動が検出された場合
には（ｎ６）、警報を発し（ｎ７）、係員に移動を告知
する。移動がない場合には、そのままルーチン動作を継
続する。
【００２３】なお、同図では、基準局３が１局，観測局
２が５局設置されているが、この設置数はこの例に限定
されず任意でよい。また、同図では基準局３，観測局２
と通信制御装置６とを有線通信で接続しているが、これ
らの間の通信を無線で行うようにしてもよい。
２が５局設置されているが、この設置数はこの例に限定
されず任意でよい。また、同図では基準局３，観測局２
と通信制御装置６とを有線通信で接続しているが、これ
らの間の通信を無線で行うようにしてもよい。
【００２４】さらに、この実施例のＧＰＳ受信機はＬ１
波，Ｌ２波の両波を受信してその搬送波位相の積算値を
出力するものであるが、そのいずれか一方（例えばＬ１
波）のみを受信してその搬送波位相の積算値を出力する
ものを適用することもできる。また、同時受信する衛星
数も４個に限定されるものではなく、多ければ多いほど
測位精度を高くすることができる。また、ＧＰＳ受信機
の構成も搬送波位相の積算値を出力することができるも
のであれば、上記実施例のものに限定されない。また、
このＧＰＳ受信機の動作電源はどのようなものでもよい
が、太陽電池，バッテリおよびバッテリ充電回路を備え
ることによりどのような場所に設置した場合でも継続的
な運用が可能になる。
波，Ｌ２波の両波を受信してその搬送波位相の積算値を
出力するものであるが、そのいずれか一方（例えばＬ１
波）のみを受信してその搬送波位相の積算値を出力する
ものを適用することもできる。また、同時受信する衛星
数も４個に限定されるものではなく、多ければ多いほど
測位精度を高くすることができる。また、ＧＰＳ受信機
の構成も搬送波位相の積算値を出力することができるも
のであれば、上記実施例のものに限定されない。また、
このＧＰＳ受信機の動作電源はどのようなものでもよい
が、太陽電池，バッテリおよびバッテリ充電回路を備え
ることによりどのような場所に設置した場合でも継続的
な運用が可能になる。
【００２５】なお、基準局３は既知の地点に設置される
ことが最も望ましいが、観測局２の変位のみを観測する
場合には、基準局３のＧＰＳ受信機で単独測位を行い、
その測位結果を基準局３の位置としてもよい。この場合
には、単独測位を複数回行ってその平均値を求めるなど
して測位の高精度化を図ればよい。
ことが最も望ましいが、観測局２の変位のみを観測する
場合には、基準局３のＧＰＳ受信機で単独測位を行い、
その測位結果を基準局３の位置としてもよい。この場合
には、単独測位を複数回行ってその平均値を求めるなど
して測位の高精度化を図ればよい。
【００２６】
【発明の効果】以上のようにこの発明によれば、複数の
地点で同時に位相積算値を計測し、それをセンタに集約
して相対測位演算を実行することができるため、ほぼリ
アルタイムに高精度の基線ベクトルを算出することが可
能になる。
地点で同時に位相積算値を計測し、それをセンタに集約
して相対測位演算を実行することができるため、ほぼリ
アルタイムに高精度の基線ベクトルを算出することが可
能になる。
【図１】この発明の実施例である測位システムの構成を
示す図
示す図
【図２】同測位システムに用いられるＧＰＳ受信機の構
成を示す図
成を示す図
【図３】同ＧＰＳ受信機の搬送波クロック抽出部の構成
を示す図
を示す図
【図４】各局がセンタに向けて送信する電文の構成を示
す図
す図
【図５】相対測位の原理を説明する図
【図６】移動検出の手法を説明する図
【図７】センタの動作を示すフローチャート
１−ＧＰＳ衛星、２−観測局、３−基準局、８−センタ
─────────────────────────────────────────────────────
フロントページの続き
(56)参考文献 特開 平６−160509（ＪＰ，Ａ)
特開 平３−68884（ＪＰ，Ａ)
特開 平５−280978（ＪＰ，Ａ)
特開 平５−312934（ＪＰ，Ａ)
(58)調査した分野(Int.Cl.7，ＤＢ名)
G01S 5/00 - 5/14
Claims (2)
- 【請求項１】 複数の地点に設けられた観測装置と、各
観測装置とオンラインで結ばれたセンタとからなる測位
システムであって、 前記観測装置は、複数のＧＰＳ衛星が送信する信号電波
を受信して、その信号または搬送波の位相の積算値をカ
ウントする手段と、カウントした位相の積算値を一定時
間毎に出力する手段とを備え、 前記センタは、各観測装置から受信した位相の積算値に
基づいて上記複数の地点間の基線ベクトルを算出するベ
クトル算出手段と、該基線ベクトルに基づき、前記観測装置の最も確からし
い位置を示す最確値および所定確率以上で前記観測装置
の位置が真となる範囲を示す信頼域を算出し、該信頼域
内に前記観測装置の初期の最確値が存在するかどうかを
判定する判定手段と、 を備えたことを特徴とする測位シ
ステム。 - 【請求項２】 前記センタは、刻々入力される位相の積
算値に基づいて前記ベクトル算出手段を繰り返して実行
し、前記判定手段で順次算出される前記信頼域内に前記
観測装置の初期の最確値が存在しないときに警報を発す
る手段を備えたことを特徴とする請求項１に記載の測位
システム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP28514794A JP3532267B2 (ja) | 1994-11-18 | 1994-11-18 | 測位システム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP28514794A JP3532267B2 (ja) | 1994-11-18 | 1994-11-18 | 測位システム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08146112A JPH08146112A (ja) | 1996-06-07 |
| JP3532267B2 true JP3532267B2 (ja) | 2004-05-31 |
Family
ID=17687710
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP28514794A Expired - Fee Related JP3532267B2 (ja) | 1994-11-18 | 1994-11-18 | 測位システム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3532267B2 (ja) |
Families Citing this family (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH1026660A (ja) * | 1996-07-10 | 1998-01-27 | Sokkia Co Ltd | 位置測定システム |
| JP3838387B2 (ja) * | 1997-03-28 | 2006-10-25 | 株式会社ニコン | 測量システム及び測量方法 |
| JP3758917B2 (ja) * | 1999-11-24 | 2006-03-22 | 日立造船株式会社 | Ｇｐｓによる物体の変位計測方法および変位計測装置 |
| KR100814251B1 (ko) * | 2000-09-20 | 2008-03-17 | 엔엑스피 비 브이 | 이동 유닛, 이동 유닛 위치 결정 방법 및 위치 결정 장치 |
| JP2003090873A (ja) * | 2001-07-10 | 2003-03-28 | Techno Vanguard:Kk | 観測装置及びｇｐｓ位置観測設備及び該観測装置を用いたｇｐｓ位置観測設備 |
| JP2003270318A (ja) * | 2002-03-15 | 2003-09-25 | Japan Radio Co Ltd | 多点観測網を利用した測位システム |
| JP3635270B2 (ja) * | 2002-07-10 | 2005-04-06 | 株式会社パスコ | 地盤変動計測システム |
| JP5043352B2 (ja) * | 2006-03-23 | 2012-10-10 | 前田建設工業株式会社 | 建造物健全性判定システム |
| JP4892304B2 (ja) * | 2006-09-01 | 2012-03-07 | 古野電気株式会社 | 観測装置および土地地盤変位監視システム |
| JP5058594B2 (ja) * | 2006-12-28 | 2012-10-24 | 株式会社日豊 | 測位点推定装置および測位点推定方法ならびにそのプログラムと記録媒体 |
| JP5028096B2 (ja) * | 2007-01-18 | 2012-09-19 | 前田建設工業株式会社 | 長大構造物の健全性判定システム |
| JP5340543B2 (ja) * | 2007-01-18 | 2013-11-13 | 前田建設工業株式会社 | 仮設構造体の健全性判定システム |
| JP5616702B2 (ja) * | 2010-07-02 | 2014-10-29 | 古野電気株式会社 | Ｇｎｓｓ解析システム、ｇｎｓｓ解析装置、及びｇｎｓｓ解析プログラム |
| JP5976413B2 (ja) * | 2012-06-22 | 2016-08-23 | 株式会社東芝 | レーダーシステム、送信機、受信機及びレーダーシステムの送受信方法 |
| JP2017058307A (ja) * | 2015-09-18 | 2017-03-23 | ソニー株式会社 | 情報処理装置、情報処理方法、およびプログラム |
| JP7044459B2 (ja) * | 2016-02-17 | 2022-03-30 | 古野電気株式会社 | 観測データ解析装置、観測データ解析システム、および、観測データ解析方法 |
| JP2020003341A (ja) * | 2018-06-28 | 2020-01-09 | パナソニックＩｐマネジメント株式会社 | データ配信サーバおよびデータ配信システム |
| JP2020003407A (ja) * | 2018-06-29 | 2020-01-09 | パナソニックＩｐマネジメント株式会社 | 構造物監視サーバおよび構造物監視システム |
-
1994
- 1994-11-18 JP JP28514794A patent/JP3532267B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH08146112A (ja) | 1996-06-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3532267B2 (ja) | 測位システム | |
| EP2312331B1 (en) | Method and system of utilizing broadcast fm signal to determine geographical position | |
| US6636744B1 (en) | Obtaining pilot phase offset time delay parameter for a wireless terminal of an integrated wireless-global positioning system | |
| CA1314964C (en) | Navigation and tracking system | |
| US8085201B2 (en) | System for determining position over a network | |
| US5173710A (en) | Navigation and positioning system and method using uncoordinated beacon signals | |
| EP0264440B1 (en) | Method and apparatus for precision dynamic differential positioning | |
| US5774829A (en) | Navigation and positioning system and method using uncoordinated beacon signals in conjunction with an absolute positioning system | |
| RU2565386C2 (ru) | Способ, устройство и система для определения позиции объекта, имеющего приемник глобальной навигационной спутниковой системы, посредством обработки неразностных данных, подобных измерениям фазы несущей, и внешних данных, подобных ионосферным данным | |
| US5347285A (en) | Method and apparatus for tracking the position and velocity of airborne instrumentation | |
| US9798017B2 (en) | Reducing time and increasing reliability of ambiguity resolution in GNSS | |
| US4048563A (en) | Carrier-modulated coherency monitoring system | |
| US20050147196A1 (en) | Time and frequency synchronisations of equipment at different locations | |
| US20070236387A1 (en) | Method of optimization of processing of location data in the presence of a plurality of satellite positioning constellations | |
| US6388613B1 (en) | Portable GPS type distance/speed meter capable of selectively using doppler speed measuring method | |
| US20080204315A1 (en) | System and Method For Monitoring and Surveying Movements of the Terrain, Large Infrastructures and Civil Building Works In General, Based Upon the Signals Transmitted by the Gps Navigation Satellite System | |
| US6720916B2 (en) | Apparatus and method for performing open loop tracking of a signal | |
| US6172638B1 (en) | Satellite signal receiver with detector of incoherence between code phase and carrier frequency measurements | |
| JPH08101265A (ja) | 測位用衛星の異常検知方法および測位システム | |
| JP3557024B2 (ja) | 測位装置 | |
| RU187712U1 (ru) | Устройство для определения ошибки слежения за временем прихода навигационного радиосигнала при его распространении через искусственное ионосферное образование | |
| EP0664008B1 (en) | Method and apparatus for tracking the position and velocity of airborne instrumentation | |
| RU93525U1 (ru) | Устройство определения степени частотно-селективных замираний навигационных радиосигналов спутниковых радионавигационных систем | |
| JP4216419B2 (ja) | 移動速度測定装置 | |
| JPS61226610A (ja) | 自動車用航法装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20031126 |
|
| 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: 20040224 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040303 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |