(57)【特許請求の範囲】
【請求項１】
空が良く見えるアンテナ（１３）を備えたベースＧＰＳレシーバ（３）と、
１以上のリモートＧＰＳレシーバ（６）と、
複数時間にわたって前記ベースＧＰＳレシーバ（３）、及び前記リモートＧＰＳレシーバ（６）から距離データを収集し、各リモートＧＰＳレシーバの単一の位置、及び関連する位置共分散を計算するデータ処理センター（１２）であって、前記データ処理センターが、前記距離データに含まれる搬送波位相測定値を処理するための浮動不明瞭度フィルタを使用して、各リモートＧＰＳレシーバから受信したデータを、データを処理する複数のパスにおいてさらにバッチ処理し、第１のパスの後、収集された距離データを使用して前回計算された推定位置を使用して前記浮動不明瞭度フィルタが初期化されることを更に特徴とするデータ処理センター（１２）と、
とを含み、
前記複数のパスのうちの前記第1のパスにおいて、前記データ処理センターは、収集されたデータを全て使用し、
前記複数のパスのうちの第２のパスにおいて、前記データ処理センターは、前記搬送波位相測定値のみを使用して、前記計算された位置、及び位置共分散を絞り込み、
前記複数のパスのうちの第３のパスにおいて、前記データ処理センターは、前記第２のパスにおいて計算された最良推定値に位置を固定し、前記搬送波位相測定値を使用して計算された二重差搬送波残差における増大を探し、前記ベースＧＰＳレシーバにより生成された測定値を使用して二重差搬送波不明瞭度を推定し、
前記データ処理センターは、小さな残差増大に関連し、又は残差増大に何も関連しないデータをさらに処理することを選択し、比較的大きな残差増大に関連するデータに無効フラグを付け、
前記データ処理センターは、前記第３のパスの結果に基づいて推定位置を再計算し、再計算された位置と前回再計算された位置との差が所定の閾値未満になるまで、前記第３のパス、及びその結果に基づく推定位置の再計算を繰り返す、ＧＰＳシステム。
【請求項２】
前記データ処理センターは、前記リモートＧＰＳレシーバが視野内に２以上のＧＰＳ衛星を有するコードエポックにおいて、前記浮動不明瞭度フィルタの位置、及び位置共分散を更新し、サイクルスリップが発生したとき、又は少なくとも１つのコードエポックについて視野内に１以下の衛星しか存在しないときに、前記位置、及び位置共分散を再初期化する、請求項１に記載のＧＰＳシステム。
【請求項３】
前記データ処理センターは、データを処理する前記複数のパスのうちの１以上において、各コードエポックについて前記距離データの有効性を判定し、有効であるものと判定されたデータだけを前記推定位置、及び位置共分散の計算において使用する、請求項１に記載のＧＰＳシステム。
【請求項４】
前記データ処理センターは、前記データに重みを割り当てる、請求項３に記載のＧＰＳシステム。
【請求項５】
前記データ処理センターは、前記ベースＧＰＳレシーバ（３）により生成された対応する測定値を使用して搬送波位相測定値の二重差を計算したときの残差の二乗の和に基づいて、前記距離データが有効であるか否かを判定する、請求項３に記載のＧＰＳシステム。
【請求項６】
前記データ処理センターは、二重差搬送波サイクル不明瞭度を推定し、搬送波サイクルスリップが検出された場合、前記搬送波不明瞭度をリセットする、請求項５に記載のＧＰＳシステム。
【請求項７】
カバレッジに制限のある環境においてデータ処理センターによりＧＰＳ位置を計算する方法であって、
１以上のベースＧＰＳレシーバ（３）から空の明瞭な眺めに関連する距離データを受信するステップと、
空の眺めに制限がある１以上のリモートＧＰＳレシーバ（６）から距離データを受信するステップと、
複数時間にわたって受信した前記距離データを、データを処理する複数のパスにおいて、浮動不明瞭度フィルタを使用してバッチ処理することによって所与のリモートＧＰＳレシーバの単一の位置、及び関連する位置共分散を判定するステップであって、前記浮動不明瞭度フィルタが、データを処理する第１のパスの後、収集された距離データを使用して前回計算された推定位置を使用して初期化されるものである、所与のリモートＧＰＳレシーバの単一の位置、及び関連する位置共分散を判定するステップと
を含み、
前記データ処理センターは、収集されたデータを全て使用し、
前記複数のパスのうちの第２のパスにおいて、前記データ処理センターは、前記搬送波位相測定値のみを使用して、前記計算された位置、及び位置共分散を絞り込み、
前記複数のパスのうちの第３のパスにおいて、前記データ処理センターは、前記第２のパスにおいて計算された最良推定値に位置を固定し、二重差搬送波残差における増大を探し、前記１以上のベースＧＰＳレシーバにより生成された測定値を使用して、二重差に基づいて不明瞭度を推定し、
前記データ処理センターは、小さな残差増大に関連する、又はは残差増大に何も関連しないデータをさらに処理することを選択し、比較的大きな残差増大に関連するデータに無効フラグを付け、
前記データ処理センターは、前記第３のパスの結果に基づいて推定位置を再計算し、再計算された位置と前回再計算された位置との差が所定の閾値未満になるまで、前記第３のパス、及びその結果に基づく推定位置の再計算を繰り返すことからなる方法。
からなる方法。
【請求項８】
前記バッチ処理するステップは、前記リモートＧＰＳレシーバが視野内に２以上のＧＰＳ衛星を有するコードエポックにおいて、前記浮動不明瞭度フィルタの位置、及び位置共分散を更新し、サイクルスリップが発生したとき、又は少なくとも１つのコードエポックについて視野内に１以下の衛星しか存在しないときは常に、前記位置、及び位置共分散を再初期化することを更に含む、請求項７に記載の方法。
【請求項９】
前記バッチ処理するステップは、データを処理する前記複数のパスのうちの１以上において、各コードエポックについて前記距離データの有効性を判定し、有効であるものと判定されたデータだけを前記推定位置、及び位置共分散の計算において使用することを更に含む、請求項７に記載の方法。
【請求項１０】
前記バッチ処理するステップは、前記データに重みを割り当てることを更に含む、請求項９に記載の方法。
【請求項１１】
前記バッチ処理するステップは、前記１以上のベースＧＰＳレシーバ（３）により生成された対応する測定値を使用して搬送波位相測定値の二重差を計算したときの残差の二乗の和に基づいて、前記距離データが有効であるか否かを判定することを更に含む，請求項９に記載の方法。
【請求項１２】
前記バッチ処理するステップは、二重差搬送波サイクル不明瞭度を推定し、搬送波サイクルスリップが検出された場合、前記搬送波不明瞭度をリセットすることを更に含む、請求項１１に記載の方法。
【請求項１３】
前記バッチ処理するステップは、ａ）最後のサイクルスリップからの時間が第１の所定の閾値を超え、ｂ）連続して取得される測定値の１つのインターバルにわたって、前記搬送波測定値の残差の二乗の和が正規化された閾値を下回り、かつｃ）連続して取得される測定値の１つのインターバルにわたって、前記残差の二乗の和の増加率が第２の所定の閾値を超えない場合に、前記距離データが有効であると判定することを更に含む、請求項１１に記載の方法。
【請求項１４】
前記空の明瞭な眺めに関連する距離データを受信するステップは、複数のベースＧＰＳレシーバが全体として空の明瞭な眺めを有する複数のベースＧＰＳレシーバから距離データを受信することを含む、請求項７に記載の方法。
【請求項１５】
前記空の明瞭な眺めに関連する距離データを受信するステップは、空が良く見える単一のベースＧＰＳレシーバ（３）から距離データを受信することを含む、請求項７に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
［発明の分野］
本発明は概してＧＰＳレシーバに関し、特に、カバレッジ（受信有効範囲）に制限がある環境で動作するＧＰＳレシーバに関する。
【背景技術】
【０００２】
［背景情報］
固定位置ＧＰＳレシーバを利用して正確な位置情報を求める測量システムその他の測定システムは、「カバレッジに制限のある」環境で動作させることが難しく、場合によっては、動作させることが出来ないこともある。カバレッジに制限のある環境とは、ＧＰＳ衛星からの直接の信号が、間接信号（マルチパス信号）によって乱され、場合によっては、直接の信号が完全に遮断される環境である。カバレッジに制限のある環境の例は、樹木の枝葉、地形、建物などによってＧＰＳレシーバからの空の眺めの一部が不明瞭となる環境や、大きなマルチパス信号がＧＰＳレシーバに供給される環境である。
【０００３】
カバレッジに制限のある環境では、ＧＰＳレシーバが受信する信号群は、乱された信号を含むだけでなく、純粋に間接的な信号、すなわち、視線上で妨害された信号に対応する反射信号を含むことがある。従来の態様で動作するＧＰＳレシーバは、受信信号を全て使用して、わずか数メートルの精度の位置を生成することができる。また、ＧＰＳレシーバは、位置の計算に使用される十分な数のＧＰＳ衛星を視界内に捕らえることが出来ないことがある。したがって、測量システムその他の測定システムは、目的によっては必要とされる誤差範囲内で位置（緯度、経度、及び高度）を判定することが出来ないことがある。
【０００４】
以下で説明するように、我々は、カバレッジに制限のある環境にあるＧＰＳレシーバから提供された距離情報を処理し、ＧＰＳレシーバの位置を必要とされる誤差範囲内で判定する方法を考え出した。
【発明の開示】
【課題を解決するための手段】
【０００５】
本発明は、ＧＰＳレシーバから得られる距離情報をバッチ処理することによって小さな誤差範囲内でＧＰＳ位置を判定するシステムである。このシステムは、空の眺めが障害物によって遮られていないベースＧＰＳレシーバを含む。測定範囲内にあるベースＧＰＳレシーバ、及び他のＧＰＳレシーバは、所与の時刻においてＧＰＳレシーバが追跡する衛星信号に基づいて距離情報を生成し、その距離情報をデータ記録処理センターに送る。この距離情報は、コードと、ＧＰＳレシーバがトラッキングしている各信号に対応する搬送波測定値との両方を含む。データ記録処理センターは、長期間、例えば、数時間、数日、又は数週間にわたって距離情報を収集した後、収集した情報をバッチ処理することによって、ＧＰＳレシーバの正確な緯度、経度、及び高度を計算する。また、後で詳しく説明するように、データ記録処理センターは、バッチ位置計算に使用される情報が十分に信頼出来るものであることを確認するために、収集された距離情報の品質の計算も行う。
【０００６】
長期間の間に、所与のＧＰＳレシーバが、３つ又は４つの比較的短いインターバルの間に異なる空の位置において少なくとも２つの衛星を同時に捕捉していれば、距離情報のバッチ処理は、測量システムや測定システムのような用途に必要とされる比較的小さな許容誤差の範囲内で、ＧＳＰレシーバの位置を生成することができる。バッチ処理により本システムは、ＧＰＳレシーバが複数の衛星からのＧＰＳ衛星信号を常に追跡していなくても、また、ＧＰＳレシーバが同じＧＰＳ衛星の組からの信号を常に追跡していなくても、ＧＰＳレシーバの正確な位置を計算することが出来る。
【０００７】
データ処理センターは、長時間にわたって所与のＧＰＳレシーバから収集された距離データをバッチ処理することにより、単一の位置（すなわち、緯度、経度、高度）、及び関連する位置共分散を計算する。バッチ処理は、収集された距離データを処理する複数のパスを含む。第１のパスは、収集されたデータを全て使用して、すなわち、全ての疑似距離情報、及び搬送波位相情報を使用して、３０メートルから６０メートルの期待精度で、推定地上位置を生成する。上で述べたように、この精度は、距離データの全体的品質の影響を受ける。
【０００８】
データを処理する次のパスにおいて、データ記録制御センターは、マルチパスの影響をあまり受けない搬送波位相測定だけを使用して、計算上の位置、及び位置共分散を絞り込む。ＧＰＳレシーバは、第１のパスから得られた推定位置、及び位置共分散を使用して開始され（集束的態様の計算である）、搬送波サイクル不明瞭度を推定し、最新の推定位置、及び関連する位置共分散を判定する。推定位置、及び関連位置共分散は、２以上の衛星がＧＰＳレシーバの視野内に存在する毎コードエポックにおいて、すなわち、搬送波サイクル不明瞭度を計算するための二重差を計算できるときに、更新される。第２のパスの終了時における位置推定精度は、３メートルから６メートルであることが期待され、大半の誤差は、高度成分に起因する。
【０００９】
データを処理する第３のパスは、第２のパスから得られた最良推定位置に位置を固定し、搬送波位相測定値を使用して計算される残差、及びベースＧＰＳレシーバからの測定値を使用して二重差に基づいて計算される推定不明瞭度の増大を探す。次にプロセスは、残差の増大が少ない、又は残差の増大に全く関与しないデータの更なる処理を選択し、残差の増大が比較的大きいデータに無効フラグを付ける。さらに、プロセスは、後の処理に使用される種々のデータに選択的に重みを付ける。
【００１０】
次にプロセスは、第３のパスの結果に基づいて推定位置を再計算する。システムは、無効フラグが付けられた測定値を計算対象から外し、残りの測定値に適当な重みを割当て、新たな位置推定値、及び関連する位置共分散を生成する。
【００１１】
次に処理システムは、位置、及び位置共分散を新たな推定値に固定し、第３のパスを繰り返す。すなわち、関連する残差に基づいて測定値の有効性、及び重みを判定する処理ステップを繰り返す。このステップにおいてシステムは、前回の推定位置、及び位置共分散に対して移動した思われるが、新たな推定値に対しては動いていないことが判明した特定の測定値を有効なものとして受け入れる場合がある。次にシステムは、最新の重み、及び有効性判定を使用して、次の推定位置、及び関連する位置共分散を判定し、この処理、すなわち、新たな重みと次の推定位置の判定を繰り返し続ける。この処理は、繰り返し行われる処理間で、推定位置の変化が所定の閾値未満になるまで継続される。その結果、所望の測定を行うために必要とされる許容範囲内の精度で位置が得られる。
【００１２】
本発明の説明では、添付の図面を参照する。
【発明を実施するための最良の形態】
【００１３】
［例示的実施形態に関する詳細な説明］
ここで図１を参照すると、位置判定システム１００は、特定地域内の選択された場所にある信号レシーバ７を通して距離情報、すなわち距離測定値、及び搬送波測定値をデータ記録制御センター１２に供給する複数のＧＰＳレシーバ６を含む。信号レシーバ７は、ケーブル１４によってセンター１２に接続される。センター１２は、距離情報を保持、及び処理するための１以上のワークステーション２、及びデータ記憶装置１を有する。種々のケーブル１４（図面には１つしか描かれていない）からの信号は、従来の態様で動作するマルチプレクサ４を通して１以上のワークステーション２に供給される。本明細書において、ＧＰＳレシーバ６は、「リモートレシーバ」と呼ばれることもある。
【００１４】
データ記録制御センター１２は、空が良く見える場所にあるベースＧＰＳアンテナ１３を使用して種々のＧＰＳ衛星９から信号を受信するベースＧＰＳレシーバ３を含む。空が良く見えるＧＰＳアンテナ１３を設置するために、高いタワー５が必要になる場合がある。図に描かれているように、ＧＰＳレシーバ６のうちの特定のもの、又は全ては、様々な樹木の下、又は樹木の付近にあり、任意の所与の時刻において、上空の特定位置にあるＧＰＳ衛星９からの信号１０は、種々のＧＰＳレシーバの位置において得られないか、又は弱められる場合がある。
【００１５】
ベースＧＰＳレシーバ３は、視野内のＧＰＳ衛星９のそれぞれからの信号１０を捕捉し、それらの信号を追跡する。ＧＰＳレシーバ６は、所与の時刻において追跡可能な衛星信号に基づいて距離情報を生成し、その距離情報を信号レシーバ７を介してデータ記録制御センター１２に送る。後で図２を参照して詳しく説明するように、センター１２は、長時間（例えば、数時間、数日、又は数週間）にわたって距離情報を収集した後、収集した情報をバッチ処理することにより、ＧＰＳレシーバの正確な緯度、経度、及び高度を計算する。また、センター１２は、バッチ処理に使用される情報が十分に信頼できるものであることを確認するために、収集された距離情報の品質の計算も行う。
【００１６】
長時間の間に、所与のＧＰＳレシーバ６が、３つ又は４つの比較的短いインターバルの間に上空の異なる位置にある少なくとも２つの衛星を同時に捕捉していて、且つ、データ収集期間中にＧＰＳレシーバが動かなかった場合に限り、距離情報のバッチ処理は、地震の測定や計測のような用途に必要とされる比較的小さな許容誤差の範囲内で、ＧＰＳレシーバの位置を計算することができる。したがって、バッチ処理によって本システムは、ＧＰＳレシーバがＧＰＳ衛星信号を常に追跡していなくても、また、ＧＰＳ衛星の同じ組からの信号を追跡していなくても、ＧＰＳレシーバの正確な位置を計算することが出来る。
【００１７】
データ記録制御センター１２は、浮動不明瞭度フィルタ８を使用してデータを処理する複数のパスを作成し、距離情報をバッチ処理する。後で詳しく説明するように、このフィルタは、搬送波不明瞭度がリセットされたときに、リセットされていない位置、及び位置共分散行列を使用して、前記複数のパスのうちの特定のパスで動作する。その結果、フィルタは、観測情報（ここでは、搬送波情報）の全履歴を使用して、所与のレシーバの静止位置を推定することが出来る。したがって、このフィルタは、連続した搬送波位相測定値の個々のブロックに関連する観測情報の使用に限定されない。
【００１８】
図２を更に参照すると、バッチ処理は、収集された距離情報を利用する次のようなパスを含む。
【００１９】
パス１：  所与のＧＰＳレシーバ６から得られるデータ（疑似距離、及び搬送波位相を含む）を浮動不明瞭度フィルタ８を使用して処理し、第１の推定位置を生成する（ステップ２００）。このパスの出力は、３０メートルから６０メートルの期待精度を持つ位置である。このデータは、長期間（すなわち、数時間、数日など）のうちの種々の時点においてＧＰＳレシーバの視野内にあるＧＰＳ衛星からの信号に基づいて生成される。通常、距離データは、８時間から２４時間の間の時間にわたって収集される。樹木、又は他の不完全な視野が原因で、リモートＧＰＳレシーバの空の眺めには制限があり、ＧＰＳレシーバは、長期間のうちの一部においてしか、視野内に同じ組の衛星を捕捉できないことがあり、また、時刻によっては、視野内に１つしか衛星を捕捉できない場合もある。また、一部の衛星信号は、ＧＰＳレシーバの空の眺めを遮る障害物や、近くの建造物などによる妨害によって反射された信号の影響を受けた多数のマルチパス成分によって、歪められる場合がある。したがって、データによっては信頼できないものもある。
【００２０】
パス２：  浮動不明瞭度フィルタ８は、パス１で計算された推定位置を使用して初期化される。パス２では、疑似距離測定値は使用せず、搬送波位相測定値だけを使用して、正確な位置を計算する（ステップ２０２）。位置、及び位置共分散は、２以上の衛星がレシーバにとって、すなわちレシーバの視野内において利用可能となる毎コードエポックにおいて更新される。サイクルスリップの後、又は少なくとも１つのコードエポックについて１以下の衛星しか利用できないときは常に、このフィルタはリセットされる。フィルタがリセットされるたびに、システムは、最後の推定位置、及び位置共分散を初期値として使用する（ステップ２０４）。このような形で、位置情報は保持され、不明瞭度情報は実質的に破棄される。最終的なパス２位置の精度は３メートルから６メートルであり、大半の誤差は高度成分に関連する。
【００２１】
パス３：  このパスでは、全ての観測搬送波データを使用できるように、位置、及び位置共分散は一定に固定される（ステップ２０６）。この固定位置は、パス２で計算された最後の位置である。位置共分散は、非常に小さな要素（すなわち、０．０００００００１ｍ２）しか持たない対角行列として保持される。位置を固定した状態で、システムは、ベースレシーバからの観測値を使用して、二重差搬送波不明瞭度を推定し、関連する二重差搬送波残差を判定する（ステップ２０８）。残差、及びその増大率の分析に基づいて、データは、無効であるもとして扱われるか、又はそれに従って重み付けされる（ステップ２１０）。
【００２２】
このパスでは、浮動不明瞭度フィルタに変更を加え、搬送波サイクルのスリップ、すなわち、ロック外れが検出された場合にのみ、搬送波不明瞭度をリセットするように浮動不明瞭度フィルタを変更する。搬送波測定値が連続していれば、残差は無制限に増大することが許される。このパスで生成された残差は、次のパスで使用される一連の信号品質インジケータの生成に使用される。信号品質インジケータによれば、各コードエポックにおいて、各衛星からの搬送波測定値に適切な重みを割り当てることが可能となる。下記の条件が満たされたとき、最大の重みが搬送波測定値に割り当てられる。
（ａ）最後のサイクルスリップからの時間が、第１の所定の閾値を超えていて、
（ｂ）連続して取得される測定値の１つのインターバルにわたって、搬送波測定値の残差の二乗の和が、正規化された閾値未満に下回っていて、
（ｃ）連続して取得される測定値の１つのインターバルにわたって、前記残差の二乗の和の増加率が、第２の所定の閾値を超えていないこと。
【００２３】
これらの条件が全て満たされると、サイクルスリップ間、又はロック破壊間のインターバル全体にわたる一連の搬送波観測値に、無効というフラグが付けられる。サイクルスリップ間のあるインターバルにおける一連の搬送波測定値が有効であると思われる場合、すなわち、それらの条件が満たされている場合、そのインターバル内の測定値の一部（インターバルの最後のポイントにおける測定値は含まない）の重みは除去される場合がある。この重みの除去は、連続した測定値における相関のあるマルチパス誤差、すなわち、ホワイトノイズを防止し、推定位置に不適切な影響が加わることを防止するために使用される。この重みの除去は、例えば、特定のコードエポックにわたって測定値の一部だけを使用する形をとる。すなわち、測定値４つにつき１つの測定値を使用するか、又は、関連計算において比較的大きな標準偏差を使用する形をとる。
【００２４】
したがって、システムは、各コードエポックにおいて、各ＰＲＮコードについてフラグを有する重み付けテーブル（図示せず）を作成する。これらのフラグは、データを処理する次のパスにおいて、測定値をどのように使用するかを知らせるためのものである。これらのフラグは、単なる「有効」又は「無効」といったものでもよいし、どの程度の量の測定値から重み付けを除去すべきかを指定するものであってもよい。
【００２５】
パス４：  最後の推定位置、及び位置共分散は、このパスにおける浮動不明瞭度フィルタの初期化に使用される（ステップ２１２）。次にフィルタは、重み付けテーブルを使用して、フィルタが使用する搬送波観測値に重みを割当て、「無効」というフラグが付けられた観測値を前記計算対象から外す。システムは、新たな推定位置、及び位置共分散を生成し（ステップ２１４）、その後、以下で説明する有効性検証を実施する。
【００２６】
パス５及びそれ以降：  パス５はパス３と同じで、最後の推定位置を固定位置として使用する（ステップ２１８）。したがって、二重差搬送波残差は、パス３と同様に計算され、テストされる。そして、対応する観測値の有効、又は無効が再度決定される。次に、新たな重み付けテーブルを作成し、新たな推定位置、及び位置共分散が、パス４と同じ仕方で計算される。次にシステムは、その新たに計算された位置が、最後に計算された位置に対して、所定の閾値を超えて異なるか否かを判定する（ステップ２１６）。異なっている場合、システムはパス５を繰り返す（ステップ２１８）。そうでない場合、システムは、新たな推定値を計算上の位置として使用する（ステップ２２０）。例えば、繰り返される処理間において、高度（大抵の誤差はここで発生する）の変化が０．０５メートル未満である場合、解は集束したものとみなされ、処理は終了する。
【００２７】
テスト説明
上で述べたプロセスをテストするために、制限が非常に厳しいテスト環境を選択した。テスト場所は、ブリティッシュコロンビア州メイプルリッジにあるカナダ・ブリティッシュコロンビア大学（ＵＢＣ）のマルコム・クナップ研究森林（Ｍａｌｃｏｌｍ  Ｋｎａｐｐ Ｒｅｓｅａｒｃｈ  Ｆｏｒｅｓｔ）である。全部で８個のＮｏｖＡｔｅｌ社のＯＥＭ４−Ｇ２Ｌ型レシーバをその森の原生林区画に一列に設置した。この森は、非常に密生していて、高さ５０メートルを超える背丈の高い樹木があり、下草も多い。また、テスト期間の大半を通じて雨が降っていたため、湿度レベルも高く、枝葉はたっぷりと水を蓄えていた。ローバＧＰＳアンテナは、地面の直ぐ近くに設置した。これらのＧＳＰアンテナは、螺旋ロッドに取り付けられ、地面から３インチ以内の位置に配置した。２２時間分のデータセットを収集した。ＮｏｖＡｓｔｅｌ社のＯＥＭ４−Ｇ２Ｌ基地局は、空の眺めが妨げられない状態に設置した。
【００２８】
ＵＢＣの研究森林は、優れたカバレッジ環境を提供するだけでなく、その中に、ＦＥＲＩＣによってテスト範囲として決められた既存の測量ネットワークを有している。従来のＧＰＳの能力は非常に貧弱で、特に、枝葉が多い環境では、高度について貧弱であるため、ＧＰＳレシーバ位置の独立した基準が必要であった。水平位置判定のための総合ステーションと、垂直成分の判定のための水準器とを使用して、ＦＥＲＩＣによるテスト範囲の外のポイントで、従来の測量を実施した。ＦＥＲＩＣ測量ネットワークは、水平方向については１０ｃｍの精度を有していると思われるが、高さについては１メートルの精度しか有していないものと思われる。このテストによって実施された水平方向の測量は、ＦＥＲＩＣネットワークが０．２０メートル未満の精度であることに合意した。レベルループは０．０１７メートル付近であるから、各地点における垂直方向の制御は、０．０１メートル未満となるはずである。本明細書に記載されている位置誤差は全て、この地上測量に関するものである。
【００２９】
テスト結果
８個のローバＧＰＳレシーバは、樹木の生い茂った渓谷の領域に配置され、大きな樹木や切り株等の近くに配置された。したがって、システムは、マルチパス信号が多数存在し、大半の場所で空の眺めに制限がある環境にＧＰＳレシーバを置いた状態でテストされた。表１は、８個のローバＧＰＳレシーバの位置誤差をまとめたものである。
【００３０】
【表１】
0005027658.tif 000002
【００３１】
水平位置の判定結果は、従来の処理技術を使用して達成されるものに比べて良好であるが、垂直位置の判定結果が非常に素晴らしく、平均高度誤差は０．５２メートルであった。
【００３２】
表２は、高度誤差と、それらの高度誤差を実現するために必要とされるパス数をまとめたものである。
【００３３】
【表２】
0005027658.tif 000003
【００３４】
表３に示されているように、高度推定値の集束に関する改善度はパスによって異なることが分かる。
【００３５】
【表３】
0005027658.tif 000004
【００３６】
テスト結果から分かるように、本システムは、レシーバが、たとえカバレッジに極めて厳しい制限のある環境を特徴とする位置に置かれた場合であっても、測量や地震計測といった用途に関係する許容誤差範囲内でＧＰＳ位置を計算する。所与のレシーバによって得られるデータを処理するためのパス数は、レシーバの置かれた位置における条件によって決まり、したがって、データの信頼性に基づいて決まる。
【００３７】
上記の説明は、本発明の特定の実施形態に制限されない。しかしながら、当然ながら、本発明の利点の一部又は全てを維持したまま、本発明に変更や修正を施すことも可能である。例えば、ベースＧＰＳレシーバ観測データは、２以上のＧＰＳレシーバによって得ることもでき、それらのＧＰＳレシーバを合わせて空全体をはっきりと見れるようにしてもよく、位置共分散は、知覚される相対移動等に関する所定の小さな閾値を表わす他の値に固定してもよい。したがって、添付の特許請求の範囲の目的は、そうした変形や修正も全て、本発明の真の思想、及び範囲に含めることにある。
【図面の簡単な説明】
【００３８】
【図１】本発明に従って構成されたシステムの機能ブロック図である。
【図２】図１のシステムによって実施されるバッチ処理を示すフロー図である。
【図１】
0005027658.tif 000005
【図２】
0005027658.tif 000006