(57)【要約】  精密単独測位とリアルタイムキネマティックを組み合わせるＰＰＰ−ＲＴＫ測位方法が開示され、当該方法は、マルチシステム航法衛星及び低軌道コンステレーションが放送する直接信号を検出すると、最初の観測データを決定することと、前記低軌道コンステレーションが放送する航法衛星補強情報、及び低軌道衛星精密軌道と精密時計誤差を受信することと、前記航法衛星補強情報、前記低軌道衛星精密軌道及び前記精密時計誤差、並びに最初の観測データを用いて、精密単独測位を行い、又は、地盤補強総合誤差補正情報を受信すると、前記航法衛星補強情報、前記低軌道衛星精密軌道及び前記精密時計誤差、前記最初の観測データ及び前記地盤補強総合誤差補正情報を用いて、地盤補強の精密単独測位を行うことと、を含む。【選択図】図１
2021526643.tif 000002
【特許請求の範囲】
【請求項１】
マルチシステム航法衛星及び低軌道コンステレーションが放送する直接信号を検出すると、最初の観測データを決定することと、
前記低軌道コンステレーションが放送する航法衛星補強情報、及び低軌道衛星精密軌道と精密時計誤差を受信することと、
前記航法衛星補強情報、前記低軌道衛星精密軌道及び前記精密時計誤差、並びに最初の観測データを用いて、精密単独測位を行い、又は、
地盤補強総合誤差補正情報を受信すると、前記航法衛星補強情報、前記低軌道衛星精密軌道及び前記精密時計誤差、前記最初の観測データ及び前記地盤補強総合誤差補正情報を用いて、地盤補強の精密単独測位を行うことと、
を含む、精密単独測位とリアルタイムキネマティックを組み合わせるＰＰＰ−ＲＴＫの測位方法。
【請求項２】
前記マルチシステム航法衛星及び低軌道コンステレーションが放送する直接信号を検出すると、最初の観測データを決定することは、
マルチシステム航法衛星及び低軌道コンステレーションが放送する直接信号を検出すると、前記直接信号を追跡し、追跡時間を決定することと、
前記追跡時間が予め設定された時間より長い場合、前記直接信号を測定し、最初の観測データを決定することと、を含む、請求項１に記載の方法。
【請求項３】
前記マルチシステム航法衛星は、米国グローバル・ポジショニング・システムＧＰＳ、中国ＢｅｉＤｏｕ、欧州連合ガリレオ、ロシアグローバル・ポジショニング・システムＧＬＯＮＡＳＳの少なくとも一種を含む、請求項１に記載の方法。
【請求項４】
前記航法衛星補強情報は、航法衛星精密軌道及び時計誤差、航法衛星位相小数偏差補正数、低軌道衛星位相小数偏差補正数、及び電離層モデルパラメータ情報の少なくとも一つを含む、請求項１に記載の方法。
【請求項５】
前記最初の観測データは、航法衛星と低軌道衛星擬似距離観測データ、航法衛星と低軌道衛星キャリア位相観測データ、及び航法衛星と低軌道衛星ドップラー観測データの少なくとも一つを含む、請求項１又は２に記載の方法。
【請求項６】
精密単独測位の処理モードは、低軌道衛星補強の不確定性浮動小数点解法モード及び低軌道衛星補強の不確定性決定解法モードの少なくとも一つを含む、請求項１に記載の方法。
【請求項７】
地盤補強総合誤差補正情報は、非差分擬似距離観測総合誤差及び非差分キャリア位相観測総合誤差の少なくとも一つを含む、請求項１に記載の方法。
【請求項８】
少なくとも１つのプロセッサと、前記少なくとも１つのプロセッサと通信可能に接続されたメモリと、を備える、精密単独測位とリアルタイムキネマティックを組み合わせるＰＰＰ−ＲＴＫの測位装置であって、
前記メモリは、前記少なくとも一つのプロセッサが実行可能なコマンドを記憶し、前記コマンドは、前記少なくとも一つのプロセッサに実行される時、請求項１〜６のいずれか一項に記載の方法を実行する、精密単独測位とリアルタイムキネマティックを組み合わせるＰＰＰ−ＲＴＫの測位装置。
【請求項９】
コンピュータ可読記憶媒体であって、
前記記録媒体は、記憶されたプログラムを含み、
前記プログラムが実行する時、請求項１〜６のいずれか一項に記載の方法を実行する、コンピュータ可読記憶媒体。
【請求項１０】
プログラムを実行するように設置され、前記プログラムは請求項１〜６のいずれか一項に記載の方法を実行する、プロセッサ。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、２０１８年０６月０４日に中国特許庁で特許出願された、出願番号が２０１８１０５６４９５２.４である中国特許出願の優先権を主張するものであり、当該出願の全開示内容は援用により本明細書に組み合わせる。
【０００２】
（技術分野）
本発明は衛星航法技術に関し、例えば、精密単独測位とリアルタイムキネマティックを組み合わせる（ＰＰＰ−ＲＴＫ）の測位方法及び装置に関する。
【背景技術】
【０００３】
航法衛星精密単独測位（Ｐｒｅｃｉｓｅ Ｐｏｉｎｔ Ｐｏｓｉｔｉ ｏｎ，ＰＰＰ）の初期化及び中断後の再初期化に長時間（３０分以上）を要することは、高速及びリアルタイムキネマティック高精度分野における該技術の応用を制限する主な要素である。初期化時間の短縮及び測位精度の向上のため、近年、整数不確定性決定技術が提案されて発展され、グローバル監視ネットワークを介して衛星位相の小数偏差補正数をリアルタイムに演算及び発行し、ユーザは偏差補正数を用いて非差分不確定性の整数特性を復元し、さらに関連成熟技術を利用して全周不確定性決定を行う。検討した結果、整数不確定性決定解法技術はＰＰＰの初期化時間を２０分間程度に短縮することができることが判明した。
【０００４】
大気遅延誤差が初期化への影響を低減するために、大気遅延制約を考慮するＰＰＰ測位方法を提案した学習者もあり、電離層モデルで生成した遅延量を観測制約情報として演算性能を改善し、初期化時間をさらに１５分間に短縮することができるが、依然として高精度でリアルタイムな測位要求を満たすことが困難である。収束時間を低減するために、現在主に地盤補強システムが非差分総合補正情報を放送する方式を用いて移動局で対応する誤差を補正し、それにより、不確定性パラメータと位置パラメータの迅速な分離を達成し、いくつかの元期内に不確定性パラメータを決定することができ、しかし、地盤監視ネットワークの配置密度に対する要求が高く、低動的ユーザに適用することが多い。
【発明の概要】
【発明が解決しようとする課題】
【０００５】
本発明は精密単独測位とリアルタイムキネマティックを組み合わせる（Ｐｒｅｃｉｓｅ Ｐｏｉｎｔ Ｐｏｓｉｔｉｏｎ - Ｒｅａｌ-Ｔｉｍｅ Ｋｉｎｅｍａｔｉｃ，ＰＰＰ−ＲＴＫ）の測位方法及び装置を提案しており、低軌道衛星高速移動特性を利用して航法信号を放送し、同時にマルチシステムを利用して航法衛星観測数を増加し、ユーザ観測空間の幾何学モデルを包括的に改善し、地盤補強監視ネットワークがある地域において、領域総合誤差情報を利用してユーザ観測誤差を補正し、ユーザの精密測位の初期化時間を低減し、且つ統一モデルによって低軌道コンステレーション航法補強システムに基づくＰＰＰ及びリアルタイムキネマティック（Ｒｅａｌ-Ｔｉｍｅ Ｋｉｎｅｍａｔｉｃ 、ＲＴＫ）サービスのシームレスな切り替えを実現する。
【課題を解決するための手段】
【０００６】
本発明が提供する低軌道コンステレーション航法補強システムに基づくＰＰＰ−ＲＴＫ測位方法であって、ステップＳ１１、ステップＳ１２及びステップＳ１３（又はステップＳ１３'）を有する。
ステップＳ１１において、マルチシステム航法衛星及び低軌道コンステレーションが放送する直接信号を検出すると、最初の観測データを決定する。
ステップＳ１２において、低軌道コンステレーションが放送する航法衛星補強情報、及び低軌道衛星精密軌道と精密時計誤差を受信する。
ステップＳ１３において、航法衛星補強情報、低軌道衛星精密軌道及び精密時計誤差、並びに最初の観測データを用いて、精密単独測位を行う。
ステップＳ１３’において、地盤補強総合誤差補正情報を受信すると、航法衛星補強情報、低軌道衛星精密軌道及び精密時計誤差、最初の観測データ及び地盤補強総合誤差補正情報を用いて、地盤補強の精密単独測位を行う。
【０００７】
一実施例では、前記マルチシステム航法衛星及び低軌道コンステレーションが放送する直接信号を検出すると、最初の観測データを決定することは、マルチシステム航法衛星及び低軌道コンステレーションが放送する直接信号を検出すると、前記直接信号を追跡し、追跡時間を決定することと、前記追跡時間が予め設定された時間より長い場合、前記直接信号を測定し、最初の観測データを決定することと、を含む。
【０００８】
一実施例では、航法衛星は、米国グローバル・ポジショニング・システム（Ｇｌｏｂａｌ Ｐｏｓｉｔｉｏｎｉｎｇ Ｓｙｓｔｅｍ，ＧＰＳ）、中国ＢｅｉＤｏｕ（北斗）、欧州連合ガリレオ、ロシアグローバル・ポジショニング・システム（ＧＬＯＢＡＬ ＮＡＶＩＧＡＴＩＯＮ ＳＡＴＥＬＬＩＴＥ ＳＹＳＴＥＭ，ＧＬＯＮＡＳＳ）の少なくとも一種を含む。
【０００９】
一実施例では、航法衛星補強情報は、航法衛星精密軌道及び時計誤差、航法衛星位相小数偏差補正数、低軌道衛星位相小数偏差補正数、及び電離層モデルパラメータ情報の少なくとも一つを含む。
【００１０】
一実施例では、最初の観測データは、航法衛星と低軌道衛星擬似距離観測データ、航法衛星と低軌道衛星キャリア位相観測データ、及び航法衛星と低軌道衛星ドップラー観測データの少なくとも一つを含む。
【００１１】
一実施例では、精密単独測位の処理モードは、低軌道衛星補強の不確定性浮動小数点解法モード及び低軌道衛星補強の不確定性決定解法モードの少なくとも一つを含む。
【００１２】
一実施例では、地盤補強総合誤差補正情報は、非差分擬似距離観測総合誤差と非差分キャリア位相観測総合誤差の少なくとも一つを含む。
【００１３】
本発明の実施例は、少なくとも１つのプロセッサと、前記少なくとも１つのプロセッサと通信可能に接続されたメモリと、を備える装置であって、
前記メモリは、前記少なくとも一つのプロセッサが実行可能なコマンドを記憶し、前記コマンドは、前記少なくとも一つのプロセッサに実行され、上記いずれか一項に記載の方法を実行する、装置を提供する。
【００１４】
本発明の実施例は、記憶されたプログラムを含み、前記プログラムが実行する時に、上記いずれか一項の精密単独測位とリアルタイムキネマティックを組み合わせる方法を実行する、コンピュータ可読記憶媒体を提供する。
【００１５】
本発明の実施例は、プログラムを実行させるように設置され、前記プログラムが実行する時に上記いずれか一項の精密単独測位とリアルタイムキネマティックを組み合わせる方法を実行させる、プロセッサを提供する。
【００１６】
本発明が提供する低軌道コンステレーション航法補強システムに基づくＰＰＰ−ＲＴＫ測位方法は、地盤補強監視ネットワークがある地域とグローバルで他の地域との間にシームレスに切り替えることができ、統一された精密単独測位モードを用いて演算する。地盤補強監視ネットワークがある地域でリアルタイム初期化、更に個別元期初期化、センチメートルレベルの測位精度を実現し、グローバルで他の地域でほぼリアルタイムに近い初期化、デシメートルレベル、さらにセンチメートルレベルの測位精度を実現する。
【００１７】
本発明が提供する低軌道コンステレーション航法補強システムに基づくＰＰＰ−ＲＴＫ測位方法は、グローバルにほぼリアルタイムに近い精密測位、速度測定及びタイムスタンプ結果を得ることができ、地盤補強領域にリアルタイムのセンチメートルレベルの測位、速度測定及びタイムスタンプ結果を得ることができ、且つ、地盤補強領域と非地盤補強のグローバル領域との間にシームレスに切り替えることができる。
【図面の簡単な説明】
【００１８】
以下の好ましい実施形態の詳細な説明を読むことにより、様々な他の利点及び長所は当業者にとって明らかになる。図面は好適な実施形態を示すための目的だけであり、本発明を限定すると考えられない。また、図面全体において、同一の構成要素には同一の参照符号を付している。
【００１９】
【図１】本発明の実施例の低軌道コンステレーション航法補強システムに基づくＰＰＰ−ＲＴＫ測位方法を示す方法フローチャートである。
【図２】本発明の実施例の低軌道コンステレーション航法補強システムに基づくＰＰＰ−ＲＴＫ測位方法を示す実装原理の概略図である。
【図３】本発明の他方の実施例の低軌道コンステレーション航法補強システムに基づくＰＰＰ−ＲＴＫ測位を示す具体的な方法のフローチャートである。
【図４】本発明の実施例が提供する装置を示すブロック構成図である。
【発明を実施するための形態】
【００２０】
以下、図面を参照しながら、本発明の例示的な実施例についてより詳細に説明する。図面には本開示の例示的な実施例が示されるが、理解されるように、ここで説明した実施例に限定されず、様々な形態で本発明を実現することができる。逆に、これらの実施例を提供することは本開示をより完全に理解するとともに、本発明の範囲を当業者に完全に伝達することができるためである。
【００２１】
図１は、本発明の実施例の低軌道コンステレーション航法補強システムに基づくＰＰＰ−ＲＴＫ測位方法を示す方法フローチャートである。図１を参照して、本発明の実施例が提供する低軌道コンステレーション航法補強システムに基づくＰＰＰ−ＲＴＫ測位方法は、ステップＳ１１、ステップＳ１２及びステップＳ１３（又はステップＳ１３'）を有する。
【００２２】
ステップＳ１１において、マルチシステム航法衛星及び低軌道コンステレーションが放送する直接信号を検出すると、最初の観測データを決定する。
【００２３】
ステップＳ１２において、低軌道コンステレーションが放送する航法衛星補強情報、及び低軌道衛星精密軌道と精密時計誤差を受信する。
【００２４】
ステップＳ１３において、航法衛星補強情報、低軌道衛星精密軌道及び精密時計誤差、並びに最初の観測データを用いて、精密単独測位を行う。
【００２５】
ステップＳ１３’において、地盤補強総合誤差補正情報を受信すると、航法衛星補強情報、低軌道衛星精密軌道及び精密時計誤差、最初の観測データ及び地盤補強総合誤差補正情報を用いて、地盤補強の精密単独測位を行う。
【００２６】
本発明実施例が提供するＰＰＰ−ＲＴＫ測位方法は、低軌道衛星が高速に移動する特性を利用して航法信号を放送し、同時に、マルチシステムを利用して航法衛星の観測数を増加し、ユーザ観測空間の幾何学モデルを包括的に改善し、グローバル範囲でほぼリアルタイムに近い初期化を実現できる。
【００２７】
一実施例では、地盤補強監視ネットワークがある地域において、該方法は地面補強監視ネットワークの配置密度を効果的に低減させ、現在領域の電離層、対流圏等の総合誤差情報を受信することによりユーザの観測誤差を補正し、統一されたＰＰＰ演算処理モードで、リアルタイム初期化を実現することができる。
【００２８】
図２は、本発明の実施例のＰＰＰ−ＲＴＫ測位方法を示す実装原理の概略図である。その処理過程は、図３に示すように、ステップ２０１〜ステップ２１０を含む。
【００２９】
ステップ２０１において、マルチシステム航法衛星及び低軌道コンステレーションが放送する直接信号を受信し、直接信号に対してキャプチャ、追跡を行う。
【００３０】
ステップ２０２において、元期毎に、航法直接信号を測定し、擬似距離、キャリア位相及びドップラー観測データを生成する。
【００３１】
ステップ２０３において、信号が安定に追跡される前提で、低軌道衛星直接信号メッセージパラメータを復調し、航法衛星補強情報と低軌道衛星精密軌道、時計誤差を取得し、前記航法衛星補強情報は、航法衛星精密軌道、時計誤差、位相小数偏差、及びグローバル電離層モデルパラメータを含む。
【００３２】
即ち、前記マルチシステム航法衛星及び低軌道コンステレーションが放送する直接信号を検出すると、最初の観測データを決定することは、マルチシステム航法衛星及び低軌道コンステレーションが放送する直接信号を検出すると、前記直接信号を追跡し、追跡時間を決定することと、前記追跡時間が予め設定された時間より長い場合、前記直接信号を測定し、最初の観測データを決定することと、を含む。
【００３３】
ステップ２０４において、最初の観測データを利用して観測方程式を確立し、１つの衛星航法システムを基準として、他の衛星航法システム及び低軌道衛星観測データに対して正規化処理を行い、統一された時間基準観測方程式を得る。
【００３４】
ステップ２０６において、非地盤補強領域に位置すれば、航法衛星補強情報及び低軌道衛星精密軌道、時計誤差を用いて観測補正を行う。
【００３５】
ステップ２０７において、地盤補強領域に位置すれば、通信リンクを介して地盤監視ネットワークが放送する非差分総合補正情報を受信する。
【００３６】
ステップ２０８において、受信した非差分総合補正情報に基づいて、ユーザの概略位置が各航法衛星及び低軌道衛星に対する誤差補正パラメータを演算する。
【００３７】
ステップ２０９において、航法衛星補強情報及び低軌道衛星精密軌道、時計誤差、及び上記演算した誤差補正パラメータを用いて観測データの補正を行う。
【００３８】
ステップ２１０において、精密単独測位モードを用いて測位処理を行い、低軌道コンステレーション航法補強システムに基づくＰＰＰ−ＲＴＫ測位、タイムスタンプ及び速度測定の結果並びにキャリア位相不確定性パラメータなどを得る。
【００３９】
本発明が提供するＰＰＰ−ＲＴＫ測位方法は、低軌道衛星が高速に移動する特性を利用して航法信号を放送し、同時にマルチシステムを利用して航法衛星の観測数を増加し、ユーザ観測空間の幾何学モデルを包括的に改善し、グローバル範囲でほぼリアルタイムに近い精密測位、速度測定及びタイムスタンプ結果を得ることができ、地盤補強領域にリアルタイムのセンチメートルレベル測位、速度測定及びタイムスタンプ結果を得ることができ、且つ地盤補強領域と非地盤補強のグローバル領域との間にシームレスに切り替えることができる。
【００４０】
以下、１つの具体的な実施例によって本発明の技術案を詳細に説明する。
【００４１】
グローバル領域の低軌道コンステレーション航法補強システムに基づく高速ＰＰＰ処理の主なプロセスは以下の通りである。
【００４２】
(１)最初の観測データを利用して観測方程式を確立する。
【００４３】
受信機が航法信号を受信して生成した最初の観測データは、マルチコンステレーション多周波数擬似距離、キャリア位相及びドップラー観測データを含み、前記擬似距離及びキャリア位相観測方程式は以下のように示される。
【００４４】
【数１】
2021526643.tif 000003
（１）
【数２】
2021526643.tif 000004
（２）
【００４５】
式中、
Ｇは、衛星航法システム及び低軌道補強システムを表す。
ｉは、信号周波数符号を表し、ｉ＝１，２，３。
ｒ,Ｓは、それぞれ受信機及び衛星符号を表す。
【数３】
2021526643.tif 000005
は、それぞれメートルを単位とする擬似距離とキャリア位相観測値である。
ρｒＧ，Ｓは、衛星から受信機までの幾何学的距離である。
ｄｔｒ，ｄｔＧ，Ｓは、それぞれ、受信機及び時計誤差である。
ｄＴｒｏｐｒＧ，Ｓは、対流圏遅延である。
ｄＩｏｎｒＧ，Ｓは、周波数１以上の電離層遅延である。
ｂｉ，ｒ，ｂｉＧ，Ｓは、それぞれ、受信機及び衛星の擬似距離ハードウェアパス遅延である。
【数４】
2021526643.tif 000006
は、全周不確定性である。
ｈｉ，ｒ，δφｉ，ｒは、それぞれ受信機キャリア位相パス遅延及び初期位相偏差である。
ｈｉＧ，Ｓ，δφｉＧ，Ｓは、それぞれ衛星キャリア位相パス遅延及び初期位相偏差である。
【数５】
2021526643.tif 000007
は、それぞれ擬似距離とキャリア位相観測にモデル化されなかった残留誤差である。
【００４６】
キャリア位相パス遅延と初期位相偏差は分離できないため、通常、２つの量を合わせ、未校正ハードウェア遅延と呼ばれ、受信機及び衛星端末はそれぞれ以下のように表される：
【数６】
2021526643.tif 000008
（３）
【数７】
2021526643.tif 000009
（４）
【００４７】
すると、キャリア位相観測方程式は以下のように表される。
【００４８】
【数８】
2021526643.tif 000010
（５）
【００４９】
(２)電離層がない組み合わせ観測値を確立する。
【００５０】
デュアルバンド観測データを用いて電離層がない組み合わせ観測を確立し、１次電離層の遅延影響を解消し、未知パラメータを低減する。具体的な組み合わせモデルは以下の通りである。
【００５１】
【数９】
2021526643.tif 000011
（６）
【数１０】
2021526643.tif 000012
（７）
【００５２】
ここで、
【数１１】
2021526643.tif 000013
【００５３】
受信機端末の擬似距離ハードウェア遅延は受信機時計誤差に吸収されるため、
【数１２】
2021526643.tif 000014
とされ、上式は以下になる。
【００５４】
【数１３】
2021526643.tif 000015
（８）
【数１４】
2021526643.tif 000016
（９）
【００５５】
マルチシステムの観測データが連携して処理する時、受信機時計誤差パラメータは受信機端末における擬似距離のパス遅延を吸収するが、パス遅延はまた信号に関連するため、異なるシステムが異なる受信機時計誤差
【数１５】
2021526643.tif 000017
に対応し、それにより、低軌道衛星及び他の衛星航法装置に対応する観測方程式は以下のように書き換えてもよい。
【００５６】
【数１６】
2021526643.tif 000018
（１０）
【数１７】
2021526643.tif 000019
（１１）
【数１８】
2021526643.tif 000020
（１２）
【数１９】
2021526643.tif 000021
（１３）
【００５７】
式中、
【数２０】
2021526643.tif 000022
である。ＧＬＯＮＡＳＳは、周波数分割多元接続技術を用いて、異なる周波数の衛星信号による受信機の擬似距離パス遅延が異なり、異なる遅延が受信機時計誤差に完全に吸収されることを引き起こす。しかし、ＧＬＯＮＡＳＳ擬似距離観測値に小さい重みを与えると、これらのパス遅延の差はほぼ残差に含まれると考えられる。従って、これらの変数は観測モデルには現れない。
【００５８】
(３)低軌道衛星が放送する航法衛星補強情報及びモデルを用いて誤差変更を行う。
【００５９】
低軌道衛星が放送する補強情報において、精密衛星軌道製品は、いずれも統一された空間座標参照基準を用い、精密衛星時計誤差製品は、統一された時間基準を用いる。従って、観測モデルにおいて座標基準や時間基準が統一されていないという問題は存在しない。同時に、精密衛星時計誤差製品は、イオンフリーコンビネーション観測値によって生成され、衛星端末擬似距離パス遅延を含む。
【００６０】
また、対流圏遅延は通常、乾燥成分と湿潤成分の二つに分けることができる。乾燥成分は、モデルによって補正され、湿潤成分は、推定対象パラメータとして推定される。評価対象パラメータの数を減らすために、マップ関数を用いてスロープ遅延を天頂方向に射影して、１つの天頂遅延のみを推定してもよい。
【数２１】
2021526643.tif 000023
（１４）
とされる。
【００６１】
提供された航法衛星補強情報、及び相対論効果、地球自転、及びアンテナ位相中心等のモデルを利用して観測方程式を補正し、一部の未知パラメータを除去するとともに、残留した衛星軌道及び時計誤差を無視し、低軌道衛星及び他の衛星航法システムに対応する観測方程式は以下になる。
【００６２】
【数２２】
2021526643.tif 000024
（１５）
【数２３】
2021526643.tif 000025
（１６）
【数２４】
2021526643.tif 000026
（１７）
【数２５】
2021526643.tif 000027
（１８）
【００６３】
式中、
【数２６】
2021526643.tif 000028
はそれぞれ低軌道衛星及び他の衛星航法システムに対応するマップ関数であり、Ｚｒは対流圏天頂湿潤遅延である。
【００６４】
(４)観測方程式の線形化を行う。
【００６５】
受信機の近似位置でテーラー展開を行い、２次項を捨て、以下の線形化観測方程式を得る。
【００６６】
【数２７】
2021526643.tif 000029
（１９）
【数２８】
2021526643.tif 000030
（２０）
【数２９】
2021526643.tif 000031
（２１）
【数３０】
2021526643.tif 000032
（２２）
【００６７】
ここで、
【数３１】
2021526643.tif 000033
（ｘＳ，ｙＳ，ｚＳ）は、低軌道衛星及び航法衛星精密軌道座標であり、（ｘｒ，０，ｙｒ，０，ｚｒ，０）は、受信機の近似位置である。これにより、観測方程式は以下のように簡略化される。
【００６８】
【数３２】
2021526643.tif 000034
（２３）
【００６９】
式中、Ｖは観測残差であり、Ａは係数行列であり、ΔＸは受信機座標補正、受信機時計誤差、対流圏の天頂湿潤遅延、キャリア位相不確定性を含む未知ベクトルであり、Ｌは計算ベクトルである。
【００７０】
(５)パラメータ推定及び不確定性決定処理を行う。
【００７１】
カルマン（Ｋａｌｍａｎ）フィルタを用いて統合ＰＰＰ処理を行う。フィルタリングにおいて、適切な観測値ランダムモデル及び状態ベクトル動的モデルを提供する必要がある。ランダムモデルは、観測値の統計的な特性を記述したものであり、通常、観測値の分散共分散行列で表される。観測方程式から分かるように、イオンフリーコンビネーション観測値は、最初の観測値の線形組合せであり、仮に異なる周波数での観測値が無相関であると、イオンフリーコンビネーション観測値の初期分散は誤差伝播の法則により演算することができる。具体的な分散は初期分散と衛星高度角の関数として定義付けられる。仮に異なる衛星、異なるシステムの観測値は無相関であり、及び、異なるタイプの観測値、すなわち、擬似距離及び位相観測値は無関係であると、観測値の分散共分散行列を得ることができる。
【００７２】
状態ベクトルの動的モデルに対して、静的受信機座標は定数と表すことができ、動的受信機座標及び受信機時計誤差はランダム遊走又は１次ガウスのマルコフ過程と表すことができ、対流圏天頂湿潤遅延はランダム遊走過程と表すことができ、キャリア位相不確定性パラメータは定数と表すことができ、それにより、状態方程式を得る。
【００７３】
【数３３】
2021526643.tif 000035
（２４）
【００７４】
式中、Ｘは推定対象の受信機座標補正、受信機時計誤差等のパラメータであり、Φは状態遷移マトリックスであり、ｗｋ−１は状態遷移騒音である。観測方程式及び状態方程式を包括的にして見れば、標準Ｋａｌｍａｎフィルタプロセスを用いてパラメータ推定を行うことができる。ここでは、衛星位相小数偏差補正を行わないため、キャリア位相不確定性浮動小数点解法の結果のみを得る。更に低軌道衛星補強情報に含まれる衛星位相小数偏差補正を利用して観測方程式の補正を行うと、不確定性の整数特性を復元し、不確定性の決定を実現し、キャリア位相不確定性決定解法の結果を得て、初期化時間を更に短縮し、測位、速度測定及びタイムスタンプ精度を向上させることができる。
【００７５】
低軌道コンステレーション航法の直接信号の観測データを増加するため、低軌道衛星の高速移動特性はユーザが観測する幾何構造を大幅に向上させ、それによりＰＰＰ初期化時間を大幅に低下させる。
【００７６】
地盤補強領域の低軌道コンステレーション航法補強に基づくＰＰＰ−ＲＴＫ処理の主なプロセスは以下の通りである。
【００７７】
地盤補強領域において、全ての参考局は三角分割（Ｄｅｌａｕｎａｙ）方法を利用して若干の三角サブネットワークに分割され、且つ非差分補正数に基づくネットワークＲＴＫ方法に従って各サブネットワークに対して各可視衛星の総合誤差補正情報をそれぞれ構築し、そのうち、各衛星方向の電離層、対流圏、及び衛星に関連するパス遅延、衛星時計誤差、衛星軌道誤差を含み、以下のように示される。
【００７８】
【数３４】
2021526643.tif 000036
（２５）
【数３５】
2021526643.tif 000037
（２６）
【数３６】
2021526643.tif 000038
は、それぞれ擬似距離とキャリア位相総合誤差補正情報を示す。
【００７９】
受信機は概略位置に基づいて周辺少なくとも３つの地盤補強局の総合誤差補正情報に対して平面フィッティングモデリングを行い、補間により得られたローカル誤差補正情報を利用してユーザの擬似距離及びキャリア位相観測値を精密化する。補正された低軌道衛星及び他の衛星航法システムに対応する観測方程式は、以下の通りに記載されてもよい。
【００８０】
【数３７】
2021526643.tif 000039
（２７）
【数３８】
2021526643.tif 000040
（２８）
【数３９】
2021526643.tif 000041
（２９）
【数４０】
2021526643.tif 000042
（３０）
【００８１】
この時、方程式における未知変数は受信機位置座標、受信機時計誤差、受信機パス遅延及びキャリア位相不確定性パラメータを含み、星間の単一差分を利用して受信機時計誤差及びパス遅延をさらに解消することができ、したがって、上記に紹介された線形化処理対策及びパラメータ推定方法を利用して受信機位置の推定及びキャリア位相不確定性の決定を行うことができる。
【００８２】
低軌道コンステレーション航法の直接信号の観測データを増加するため、低軌道衛星の高速移動特性はユーザが観測する幾何構造を大幅に向上させ、それにより、同じ初期化時間及び測位精度の要求で、地盤補強監視ネットワークの配置密度を大幅に低減させ、監視ネットワークの建設コストを低減させる。
【００８３】
グローバル領域及び地盤補強領域の低軌道コンステレーション補強高速速度測定、タイムスタンプ演算処理過程は、測位処理過程と類似し、ここで説明を省略する。
【００８４】
本発明の実施例が提供するＰＰＰ−ＲＴＫ測位方法は、低軌道衛星によって航法直接信号を放送し、低軌道衛星の高速移動特性を利用してユーザ観測空間の幾何学モデルを包括的に改善し、ユーザのＰＰＰ初期化時間を準リアルタイムに短縮することができる。地盤補強領域において、現在領域の総合誤差情報をさらに受信することによってユーザ観測誤差を補正し、ＰＰＰと統一された演算処理モードを採用し、初期化時間をリアルタイムにさらに短縮する。該方法は地面補強監視ネットワークの配置密度を効果的に低減させ、且つ統一されたモデルによって低軌道コンステレーション補強マルチシステムＰＰＰとＲＴＫサービスのシームレスな切り替えを実現する。
【００８５】
図４は、本実施例が提供する装置を示すブロック構成図である。本実施例が提供する装置は、プロセッサ(ｐｒｏｃｅｓｓｏｒ)４０１およびメモリ(ｍｅｍｏｒｙ)４０３を含んでもよく、さらに、通信インターフェース(Ｃｏｍｍｕｎｉｃａｔｉｏｎｓ Ｉｎｔｅｒｆａｃｅ)４０２およびバス４０４を含んでもよい。前記プロセッサ４０１、通信インターフェース４０２、メモリ４０３は、バス４０４を介して相互に通信を行うことができる。通信インターフェース４０２は、情報伝送に用いてもよい。プロセッサ８０１は、上記実施例の精密単独測位とリアルタイムキネマティックを組み合わせるＰＰＰ−ＲＴＫの測位方法を実行させるように、メモリ４０３内の論理コマンドを呼び出すことができる。前記方法は、マルチシステム航法衛星及び低軌道コンステレーションが放送する直接信号を検出すると、最初の観測データを決定することと、前記低軌道コンステレーションが放送する航法衛星補強情報、及び低軌道衛星精密軌道と精密時計誤差を受信することと、前記航法衛星補強情報、前記低軌道衛星精密軌道及び前記精密時計誤差、並びに最初の観測データを用いて、精密単独測位を行い、又は、地盤補強総合誤差補正情報を受信すると、前記航法衛星補強情報、前記低軌道衛星精密軌道及び前記精密時計誤差、前記最初の観測データ及び前記地盤補強総合誤差補正情報を用いて、地盤補強の精密単独測位を行うことと、を含む。
【００８６】
本発明実施例は記録媒体を更に提供し、記録媒体は、記憶したプログラムを含み、前記プログラムが実行するとき、上記実施例が提供する方法を実行する。前記方法は、
マルチシステム航法衛星及び低軌道コンステレーションが放送する直接信号を検出すると、最初の観測データを決定することと、前記低軌道コンステレーションが放送する航法衛星補強情報、及び低軌道衛星精密軌道と精密時計誤差を受信することと、前記航法衛星補強情報、前記低軌道衛星精密軌道及び前記精密時計誤差、並びに最初の観測データを用いて、精密単独測位を行い、又は、地盤補強総合誤差補正情報を受信すると、前記航法衛星補強情報、前記低軌道衛星精密軌道及び前記精密時計誤差、前記最初の観測データ及び前記地盤補強総合誤差補正情報を用いて、地盤補強の精密単独測位を行うことと、を含む。
【００８７】
本発明は、プログラムを実行させるように設置されるプロセッサを更に提供し、前記プログラムは、上記実施例が提供する方法を実行させるように前記プロセッサに設置し、前記方法は、
マルチシステム航法衛星及び低軌道コンステレーションが放送する直接信号を検出すると、最初の観測データを決定することと、前記低軌道コンステレーションが放送する航法衛星補強情報、及び低軌道衛星精密軌道と精密時計誤差を受信することと、前記航法衛星補強情報、前記低軌道衛星精密軌道及び前記精密時計誤差、並びに最初の観測データを用いて、精密単独測位を行い、又は、地盤補強総合誤差補正情報を受信すると、前記航法衛星補強情報、前記低軌道衛星精密軌道及び前記精密時計誤差、前記最初の観測データ及び前記地盤補強総合誤差補正情報を用いて、地盤補強の精密単独測位を行うことと、を含む。
【００８８】
以上の実施形態の説明により、当業者が明らかに分かるように、複数の実施形態はソフトウェア及び必要な汎用ハードウェアプラットフォームの方式によって実現することができ、当然、ハードウェアによって実現することができる。このような理解に基づき、上記技術案の本質的に言うか、関連技術に寄与する部分はソフトウェア製品の形態で実現することができ、該コンピュータソフトウェア製品は、例えば、読み取り専用メモリ（Ｒｅａｄ-Ｏｎｌｙ Ｍｅｍｏｒｙ 、ＲＯＭ）／ランダムアクセスメモリ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ 、ＲＡＭ）、磁気ディスク、光ディスク等のようなコンピュータ可読記憶媒体に記憶することができ、コンピュータ装置（パソコン、サーバ、又はネットワーク装置等であってもよい）に複数の実施例又は実施例の幾つかの部分に記載の方法を実行させるように、複数のコマンドを含む。
【００８９】
また、当業者であれば理解できるように、ここでのいくつかの実施例は他の特徴ではなく、他の実施例に含まれるいくつかの特徴を含むが、異なる実施例の特徴の組み合わせは本発明の範囲内にあり且つ異なる実施例を形成することを意味する。例えば、以下の特許請求の範囲において、特許請求する実施例のいずれかは任意の組み合わせ方式で使用することができる。
【００９０】
（付記）
（付記１）
マルチシステム航法衛星及び低軌道コンステレーションが放送する直接信号を検出すると、最初の観測データを決定することと、
前記低軌道コンステレーションが放送する航法衛星補強情報、及び低軌道衛星精密軌道と精密時計誤差を受信することと、
前記航法衛星補強情報、前記低軌道衛星精密軌道及び前記精密時計誤差、並びに最初の観測データを用いて、精密単独測位を行い、又は、
地盤補強総合誤差補正情報を受信すると、前記航法衛星補強情報、前記低軌道衛星精密軌道及び前記精密時計誤差、前記最初の観測データ及び前記地盤補強総合誤差補正情報を用いて、地盤補強の精密単独測位を行うことと、
を含む、精密単独測位とリアルタイムキネマティックを組み合わせるＰＰＰ−ＲＴＫの測位方法。
【００９１】
（付記２）
前記マルチシステム航法衛星及び低軌道コンステレーションが放送する直接信号を検出すると、最初の観測データを決定することは、
マルチシステム航法衛星及び低軌道コンステレーションが放送する直接信号を検出すると、前記直接信号を追跡し、追跡時間を決定することと、
前記追跡時間が予め設定された時間より長い場合、前記直接信号を測定し、最初の観測データを決定することと、を含む、付記１に記載の方法。
【００９２】
（付記３）
前記マルチシステム航法衛星は、米国グローバル・ポジショニング・システムＧＰＳ、中国ＢｅｉＤｏｕ、欧州連合ガリレオ、ロシアグローバル・ポジショニング・システムＧＬＯＮＡＳＳの少なくとも一種を含む、付記１に記載の方法。
【００９３】
（付記４）
前記航法衛星補強情報は、航法衛星精密軌道及び時計誤差、航法衛星位相小数偏差補正数、低軌道衛星位相小数偏差補正数、及び電離層モデルパラメータ情報の少なくとも一つを含む、付記１に記載の方法。
【００９４】
（付記５）
前記最初の観測データは、航法衛星と低軌道衛星擬似距離観測データ、航法衛星と低軌道衛星キャリア位相観測データ、及び航法衛星と低軌道衛星ドップラー観測データの少なくとも一つを含む、付記１又は２に記載の方法。
【００９５】
（付記６）
精密単独測位の処理モードは、低軌道衛星補強の不確定性浮動小数点解法モード及び低軌道衛星補強の不確定性決定解法モードの少なくとも一つを含む、付記１に記載の方法。
【００９６】
（付記７）
地盤補強総合誤差補正情報は、非差分擬似距離観測総合誤差及び非差分キャリア位相観測総合誤差の少なくとも一つを含む、付記１に記載の方法。
【００９７】
（付記８）
少なくとも１つのプロセッサと、前記少なくとも１つのプロセッサと通信可能に接続されたメモリと、を備える、精密単独測位とリアルタイムキネマティックを組み合わせるＰＰＰ−ＲＴＫの測位装置であって、
前記メモリは、前記少なくとも一つのプロセッサが実行可能なコマンドを記憶し、前記コマンドは、前記少なくとも一つのプロセッサに実行される時、付記１〜６のいずれか一つに記載の方法を実行する、精密単独測位とリアルタイムキネマティックを組み合わせるＰＰＰ−ＲＴＫの測位装置。
【００９８】
（付記９）
コンピュータ可読記憶媒体であって、
前記記録媒体は、記憶されたプログラムを含み、
前記プログラムが実行する時、付記１〜６のいずれか一つに記載の方法を実行する、コンピュータ可読記憶媒体。
【００９９】
（付記１０）
プログラムを実行するように設置され、前記プログラムは付記１〜６のいずれか一つに記載の方法を実行する、プロセッサ。
【符号の説明】
【０１００】
１００  航法コンステレーション
１０１  ＢｅｉＤｏｕ（北斗）航法衛星
１０２  ＧＰＳ航法衛星
１０３  ＧＬＯＮＡＳＳ航法衛星
１０４  他の航法システム衛星
１１０  低軌道補強コンステレーション
１１１  低軌道衛星
１２０  航法信号であり、航法衛星航法信号及び低軌道衛星航法信号を含む
１３０  グローバル地域
１３１  低軌道コンステレーション航法補強システムに基づくＰＰＰ−ＲＴＫ測位方法を用いる航法装置
１４０  地盤補強領域
１４１  地盤補強監視局
１４２  地盤補強情報放送装置
【図１】
2021526643.tif 000043
【図２】
2021526643.tif 000044
【図３】
2021526643.tif 000045
【図４】
2021526643.tif 000046
【国際調査報告】
2021526643.tif 000047
2021526643.tif 000048
2021526643.tif 000049
2021526643.tif 000050
2021526643.tif 000051