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