長期予測軌道データの信頼性判定方法、長期予測軌道データの提供方法及び情報提供装置
【課題】長期予測エフェメリスの信頼性を判定するための手法を提案すること。
【解決手段】測位システム1において、サーバシステム3は、外部システム2から受信したGPS衛星SVの位置を時系列に予測した衛星予測暦に含まれる各予測位置を、測位用衛星の実際の位置を時系列に記憶した衛星精密暦に含まれる実績位置と比較することで、衛星予測暦に含まれる予測位置の正確性の変化を分析する。そして、連続する複数の予測対象期間それぞれにおける予測衛星軌道の衛星軌道パラメータの値を纏めてなる長期予測エフェメリスの予測対象期間それぞれの信頼性を、分析の結果に基づいて判定する。そして、信頼性の判定結果を長期予測エフェメリスに含めて、携帯型電話機4に提供する。
【解決手段】測位システム1において、サーバシステム3は、外部システム2から受信したGPS衛星SVの位置を時系列に予測した衛星予測暦に含まれる各予測位置を、測位用衛星の実際の位置を時系列に記憶した衛星精密暦に含まれる実績位置と比較することで、衛星予測暦に含まれる予測位置の正確性の変化を分析する。そして、連続する複数の予測対象期間それぞれにおける予測衛星軌道の衛星軌道パラメータの値を纏めてなる長期予測エフェメリスの予測対象期間それぞれの信頼性を、分析の結果に基づいて判定する。そして、信頼性の判定結果を長期予測エフェメリスに含めて、携帯型電話機4に提供する。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、長期予測軌道データの信頼性判定方法、長期予測軌道データの提供方法及び情報提供装置に関する。
【背景技術】
【0002】
測位用信号を利用した測位システムとしては、GPS(Global Positioning System)が広く知られており、携帯型電話機やカーナビゲーション装置等に内蔵された測位装置に利用されている。GPSでは、複数のGPS衛星の位置や各GPS衛星から自機までの擬似距離等の情報に基づいて自機の位置を示す3次元の座標値と時計誤差とを求める測位演算を行う。
【0003】
GPSによる測位では、先ず、GPS衛星から発信されるGPS衛星信号に重畳されているアルマナックやエフェメリスといった航法データに基づいて、GPS衛星の位置、速度、移動方向等の衛星情報を算出する。アルマナックは衛星を捕捉する際の有力な手掛かりとはなるが、衛星情報の精度が悪いため測位演算に使用しないのが一般的である。一方、エフェメリスは、衛星情報の精度が良いため、衛星を捕捉する際の有力な手掛かりとなるだけでなく、測位演算にも使用することができる。従って、例えばエフェメリスを保持していない状態で測位を開始した場合には、エフェメリスをGPS衛星信号から取得しなければならず、初回測位時間(TTFF:Time To First Fix)が増大する。
【0004】
そこで、サーバクライアントシステムであって、1週間といった長期間分のエフェメリスである長期予測エフェメリス(長期予測軌道データ)をサーバが予測し、クライアントである測位装置に提供する技術が開発されており、例えば特許文献1や特許文献2に開示されている。
【特許文献1】米国特許出願公開第2002/0188403号明細書
【特許文献2】米国特許出願公開第2005/0212700号明細書
【発明の開示】
【発明が解決しようとする課題】
【0005】
長期予測エフェメリスの定義方法としては、通常のエフェメリスと同様のデータフォーマットで定義する方法が考えられる。すなわち、衛星軌道の近似モデルの1つであるケプラーの楕円軌道モデルを用いて衛星軌道を近似し、その時のモデル式のパラメータ(以下、「衛星軌道パラメータ」と称す。)の値によって、長期予測エフェメリスを定義する方法である。測位用衛星の将来の位置を所定時間おきに時系列に予測した予測位置でなる衛星予測暦(予測位置データ)は、所定の商用システムから提供されている。ケプラーの楕円軌道モデルによる近似計算は、この衛星予測暦を用いて行うことが可能である。
【0006】
しかし、衛星予測暦に含まれる測位用衛星の予測位置は、将来になるほど、測位用衛星の実際位置からずれる傾向がある。そのため、ケプラーの楕円軌道モデルによる近似計算を行って長期予測エフェメリスを生成する場合に、近似計算により求めた衛星軌道は、生成日時から将来のものであるほど、実際の衛星軌道からずれたものとなる可能性がある。従来は、長期予測エフェメリスを一旦生成した後は、生成した長期予測エフェメリスの信頼性(測位に使用するのに適しているか否か)を判定する手法が存在しなかったため、サーバから長期予測エフェメリスを取得した測位装置は、信頼性の低い長期予測エフェメリスを使用して測位を行う場合があり得た。
【0007】
本発明は、上述した課題に鑑みて為されたものであり、長期予測エフェメリスの信頼性を判定するための手法を提案することを目的としている。
【課題を解決するための手段】
【0008】
以上の課題を解決するための第1の発明は、測位用衛星の位置を時系列に予測した予測位置データに含まれる各予測位置を、当該予測位置に対応する前記測位用衛星の実績位置と比較することで、前記予測位置データに含まれる時系列の予測位置の正確性の変化を分析することと、連続する複数の予測対象期間それぞれにおける予測衛星軌道を纏めてなる長期予測軌道データの前記予測対象期間それぞれの信頼性を、前記分析の結果に基づいて判定することと、を含む長期予測軌道データの信頼性判定方法である。
【0009】
また、他の発明として、連続する複数の予測対象期間それぞれにおける予測衛星軌道を纏めてなる長期予測軌道データを生成する生成部と、測位用衛星の位置を時系列に予測した予測位置データに含まれる各予測位置を、当該予測位置に対応する前記測位用衛星の実績位置と比較することで、前記予測位置データに含まれる時系列の予測位置の正確性の変化を分析する分析部と、前記生成部により生成された長期予測軌道データの前記予測対象期間それぞれの信頼性を、前記分析部の分析結果に基づいて判定する判定部と、前記生成部により生成された長期予測軌道データと前記判定部による判定結果とを測位装置に提供する提供部と、を備えた情報提供装置を構成してもよい。
【0010】
この第1の発明等によれば、測位用衛星の位置を時系列に予測した予測位置データに含まれる各予測位置を、当該予測位置に対応する測位用衛星の実績位置と比較することで、予測位置データに含まれる時系列の予測位置の正確性の変化を分析する。そして、連続する複数の予測対象期間それぞれにおける予測衛星軌道を纏めてなる長期予測軌道データの予測対象期間それぞれの信頼性を、分析結果に基づいて判定する。
【0011】
測位用衛星の実績位置は、測位用衛星の実際にいた位置である。そのため、測位用衛星の予測位置と実績位置との差が大きい場合は、予測位置の正確性が低く、この予測位置に基づいて予測した予測衛星軌道の信頼性は低いと言える。従って、予測位置の正確性の変化の分析結果に基づくことで、長期予測軌道データの信頼性を適切に判定することが可能となる。
【0012】
また、第2の発明として、第1の発明の長期予測軌道データの信頼性判定方法であって、前記分析することは、予測した日時が異なる複数の前記予測位置データそれぞれについて、当該予測位置データに含まれる各予測位置を対応する実績位置と比較することと、前記比較の結果を統計処理して、前記測位用衛星の位置を時系列に予測した場合の時間経過に対する予測位置の正確性の変化パターンを算出することと、を含み、前記信頼性を判定することは、前記変化パターンに基づいて、前記長期予測軌道データの前記予測対象期間それぞれの信頼性を判定することである長期予測軌道データの信頼性判定方法を構成してもよい。
【0013】
この第2の発明によれば、予測した日時が異なる複数の予測位置データそれぞれについて、各予測位置を対応する実績位置と比較する。そして、比較の結果を統計処理して予測位置の正確性の変化パターンを算出する。この変化パターンは、いわば予測位置の正確性の変化の様子を表す標準パターンとも言える。そして、算出した変化パターンに基づいて、長期予測軌道データの予測対象期間それぞれの信頼性を判定する。
【0014】
また、第3の発明として、第2の発明の長期予測軌道データの信頼性判定方法であって、前記予測位置データには、複数の前記測位用衛星それぞれの予測位置が含まれ、前記比較することは、前記測位用衛星それぞれの各予測位置を実績位置と比較することであり、前記変化パターンを算出することは、前記測位用衛星それぞれの前記比較の結果を統計処理して、前記測位用衛星間で標準的な前記変化パターンを算出することである長期予測軌道データの信頼性判定方法を構成してもよい。
【0015】
この第3の発明によれば、測位用衛星それぞれの各予測位置を実績位置と比較し、測位用衛星それぞれの比較の結果を統計処理して、測位用衛星間で標準的な前記変化パターンを算出する。
【0016】
また、第4の発明として、第3の発明の長期予測軌道データの信頼性判定方法であって、前記長期予測軌道データは、前記予測対象期間それぞれについての各測位用衛星毎のデータでなり、前記測位用衛星それぞれの前記比較の結果を前記測位用衛星間で相対評価することを更に含み、前記信頼性を判定することは、前記相対評価の結果に基づいて前記変化パターンを前記測位用衛星毎に補正して、前記長期予測軌道データの前記予測対象期間それぞれの信頼性を各測位用衛星毎に判定することを含む長期予測軌道データの信頼性判定方法を構成してもよい。
【0017】
この第4の発明によれば、測位用衛星それぞれの比較の結果を測位用衛星間で相対評価し、相対評価の結果に基づいて変化パターンを測位用衛星毎に補正して、長期予測軌道データの予測対象期間それぞれの信頼性を各測位用衛星毎に判定する。
【0018】
予測位置と実績位置との比較結果を測位用衛星間で相対評価した結果、相対的に比較結果が良いと判定した測位用衛星については、長期予測軌道データの信頼性が高くなるように変化パターンを補正する。逆に、相対的に比較結果が悪いと判定した測位用衛星については、長期予測軌道データの信頼性が低くなるように変化パターンを補正する。これにより、測位用衛星間の比較結果を考慮に入れた適切な信頼性の判定を実現し得る。
【0019】
また、第5の発明として、第3又は第4の発明の長期予測軌道データの信頼性判定方法であって、前記分析することは、前記測位用衛星それぞれについて、前記予測位置データの予測日時毎に、当該予測位置データに関する前記比較の結果を纏めることを含み、前記信頼性を判定することは、前記予測位置データの予測日時毎に纏めた前記比較の結果の前記予測日時の時系列変化に基づいて前記変化パターンを前記測位用衛星毎に補正して、前記長期予測軌道データの前記予測対象期間それぞれの信頼性を各測位用衛星毎に判定することを含む長期予測軌道データの信頼性判定方法を構成してもよい。
【0020】
この第5の発明によれば、測位用衛星それぞれについて、予測位置データの予測日時毎に、当該予測位置データに関する比較の結果を纏める。そして、予測位置データの予測日時毎に纏めた比較の結果の予測日時の時系列変化に基づいて変化パターンを測位用衛星毎に補正して、長期予測軌道データの予測対象期間それぞれの信頼性を各測位用衛星毎に判定する。
【0021】
予測位置と実績位置との比較結果が時間経過に伴って良くなる測位用衛星については、長期予測軌道データの信頼性が高くなるように変化パターンを補正する。逆に、比較結果が時間経過に伴って悪くなる測位用衛星については、長期予測軌道データの信頼性が低くなるように変化パターンを補正する。これにより、比較結果の時間変化を考慮に入れた適切な信頼性の判定を実現し得る。
【0022】
また、第6の発明として、第5の発明の長期予測軌道データの信頼性判定方法であって、前記信頼性を判定することは、前記複数の測位用衛星のうち、当該測位用衛星の予測位置を実績位置と比較した前記比較の結果が所定の不良条件を満たさない測位用衛星についてのみ前記変化パターンの前記補正を行った上で信頼性の判定を行うことを含む長期予測軌道データの信頼性判定方法を構成してもよい。
【0023】
この第6の発明によれば、複数の測位用衛星のうち、当該測位用衛星の予測位置を実績位置と比較した比較の結果が所定の不良条件を満たさない測位用衛星についてのみ変化パターンの補正を行った上で信頼性の判定を行う。予測位置と実績位置との比較結果が所定の不良条件を満たす測位用衛星については、予測位置が実績位置から大きくずれる何らかの原因(例えば衛星の軌道修正)が発生したものと推定して、変化パターンを補正しないことにしたものである。
【0024】
また、第7の発明として、長期予測軌道データを生成することと、第1〜第6の何れかの発明の長期予測軌道データの信頼性判定方法を用いて、前記生成した長期予測軌道データの信頼性を判定することと、前記生成した長期予測軌道データと前記判定した判定結果とを測位装置に提供することと、を含む長期予測軌道データの提供方法を構成してもよい。
【0025】
この第7の発明によれば、長期予測軌道データを生成し、長期予測軌道データの信頼性判定方法を用いて、生成した長期予測軌道データの信頼性を判定する。そして、生成した長期予測軌道データと信頼性の判定結果とを測位装置に提供する。これにより、測位装置は、提供された長期予測軌道データを用いて測位を行うが、この際、信頼性の低いデータを測位に使用しないようにすることができる。
【発明を実施するための最良の形態】
【0026】
以下、図面を参照して、本発明に好適な実施形態の一例を説明する。但し、本発明を適用可能な実施形態がこれに限定されるわけではない。
【0027】
1.システム構成
図1は、本実施形態における測位システム1の概略構成を示す図である。測位システム1は、外部システム2と、情報提供装置の一種であるサーバシステム3と、測位装置を備えた電子機器の一種である携帯型電話機4と、測位用衛星の一種である複数のGPS衛星SV(SV1,SV2,SV3,SV4,・・・)とを備えて構成される。尚、携帯型電話機4が必要なデータをサーバシステム3から取得した後は、携帯型電話機4とGPS衛星SVとで測位が可能であるため、携帯型電話機4とGPS衛星SVとで1つの測位システムが構成されるということもできる。また、地上側のシステムとして、サーバシステム3と携帯型電話機4とで測位システムと呼ぶこともできる。
【0028】
外部システム2は、GPS衛星SVから定期的に衛星信号を受信し、当該衛星信号に含まれる航法データ等に基づいて衛星予測暦を生成してサーバシステム3に提供する公知のシステムである。外部システム2が提供する衛星予測暦は、各GPS衛星SVそれぞれについて、将来の位置を予測した予測位置及びGPS衛星SVに搭載された原子時計の誤差を予測した時計予測誤差を所定時間おき(例えば15分おき)に時系列に並べた位置のデータである。
【0029】
また、外部システム2は、将来のデータとしての衛星予測暦を提供する他に、過去の事実のデータも提供する。すなわち、外部システム2は、GPS衛星SVの実際の位置である実績位置及びGPS衛星SVに搭載された原子時計の実際の誤差である時計実績誤差を含む衛星精密暦を過去の事実のデータとして生成して、サーバシステム3に提供する。実績位置及び時計実績誤差の算出方法については公知であるため、詳細な説明を省略する。外部システム2は、例えば衛星予測暦や衛星精密暦の提供を業務とする民間や公営の団体のコンピュータシステムに相当する。
【0030】
サーバシステム3は、衛星予測暦及び衛星精密暦を外部システム2から取得し、当該衛星予測暦及び衛星精密暦を用いて、全てのGPS衛星SVの予測されるエフェメリスであって、少なくとも1日以上の例えば1週間といった長期間有効なエフェメリス(以下、本実施形態において「長期予測エフェメリス」と称す。長期間有効な軌道でもあるため、長期予測軌道データとも言える。)を生成・提供するサーバを備えたシステムである。
【0031】
携帯型電話機4は、ユーザが通話やメールの送受信等を行うための電子機器であり、通話やメールの送受信といった携帯型電話機としての本来の機能の他、位置を計測する機能(測位機能)をなす測位装置を具備している。携帯型電話機4は、ユーザ操作に従って、サーバシステム3に対して長期予測エフェメリスの要求信号を送信し、サーバシステム3から長期予測エフェメリスを受信する。そして、受信した長期予測エフェメリスを用いてGPS衛星SVを捕捉し、衛星信号に基づく測位演算を実行する。
【0032】
2.原理
サーバシステム3は、外部システム2から取得した衛星予測暦を用いて、長期予測エフェメリスを生成する処理を行う。具体的には、長期予測エフェメリスの生成日時を基準として1週間後までの期間を「生成対象期間」とし、当該生成対象期間を衛星軌道の近似・モデル化を行う複数の期間(以下、「予測対象期間」と称す。)に区切る。本実施形態では、予測対象期間の長さを一律に6時間とする。すなわち、1週間の生成対象期間を6時間毎に28個の予測対象期間(第1予測対象期間〜第28予測対象期間)に区切る。
【0033】
そして、サーバシステム3は、外部システム2から取得した衛星予測暦に含まれる予測位置のうち、各予測対象期間に対応する期間の予測位置を抽出する。そして、抽出した全ての予測位置からの距離の二乗和が最小となるようなケプラーの衛星軌道モデル式(以下、「近似モデル」ともいう。)を各予測対象期間それぞれについて求める。このとき求めた衛星軌道の近似モデル式のパラメータを「衛星軌道パラメータ」と称し、近似モデルを算出する計算を「近似計算」ともいう。また、近似計算により求めた衛星軌道のことを「予測軌道」と称する。長期予測エフェメリスは、全てのGPS衛星の全ての予測対象期間の衛星軌道パラメータの値が格納されたデータである。
【0034】
衛星予測暦に含まれるGPS衛星SVの予測位置は、将来になるほど、GPS衛星SVの実際の位置からずれる傾向がある。そのため、近似計算を行って長期予測エフェメリスを生成した場合に、近似計算により求めた予測軌道は、生成日時から将来のものであるほど、実際の衛星軌道からずれたものとなる可能性がある。
【0035】
本実施形態では、サーバシステム3は、各GPS衛星の各予測対象期間について、予測軌道の信頼性を示す指標値である「予測軌道信頼度」を決定し、決定した予測軌道信頼度を、信頼性パラメータとして長期予測エフェメリスに含めて携帯型電話機4に提供する。本実施形態では、予測軌道信頼度は「0」〜「12」の13段階で表され、「0」が予測軌道の信頼性が最も高く、「12」が予測軌道の信頼性が最も低いことを示している。尚、予測軌道信頼度の数値範囲は適宜設定変更可能であり、例えば「0」〜「15」の16段階で表すこととしてもよい。予測軌道信頼度は、エフェメリスに含まれている「URA index」に相当する値である。
【0036】
予測軌道信頼度の具体的な決定方法について、図面を参照して説明する。本実施形態では、開始日時の異なる複数の衛星予測暦と、それらの衛星予測暦に対応する衛星精密暦(開始日時が同一である衛星精密暦)との組合せ(以下、「衛星暦組合せ」と称す。)を用いて、長期予測エフェメリスに含める予測軌道信頼度を決定する。開始日時とは、衛星予測暦及び衛星精密暦にそれぞれ含まれる衛星位置のうち最も古い衛星位置の日時である。
【0037】
本実施形態では、開始日時が4時間ずつずれた衛星予測暦と、それらに対応する衛星精密暦とを用いて予測軌道信頼度を決定する。また、衛星予測暦及び衛星精密暦には、それぞれ1週間分の衛星位置及び時計誤差のデータが含まれており、便宜的に、この1週間を6時間毎にグループ分けした28個の期間を「第1期間〜第28期間」と表現する。但し、実際には、衛星予測暦及び衛星精密暦のデータは、衛星位置及び時計誤差が羅列的に配列されたデータ構造となっており、期間毎にグループ分けされているわけではない。
【0038】
より具体的には、例えば図2に示すように、現在日時から見て、8日と0時間前を開始日時とする衛星暦組合せ、7日と20時間前を開始日時とする衛星暦組合せ、7日と16時間前を開始日時とする衛星暦組合せ、・・・、1日と0時間前を開始日時とする衛星暦組合せといった具合に、開始日時が異なるN個の衛星予測暦組合せを抽出する。図2において、1つの帯が1つの衛星予測暦及び衛星精密暦を示している。そして、抽出した衛星予測暦に含まれる予測位置と、衛星精密暦に含まれる実績位置とを比較し、予測位置がどれだけ正確であるかを分析することで、予測軌道の信頼性を判定する。
【0039】
(1)観測誤差の算出
予測軌道信頼度を決定するに当たり、先ず、抽出した各衛星暦組合せそれぞれについて、当該衛星予測暦に含まれる各日時(衛星予測暦に含まれる15分おきの日時)における予測位置と、当該衛星精密暦に含まれる各日時(衛星精密暦に含まれる15分おきの日時)における実績位置とを用いて「観測誤差」を算出する。
【0040】
具体的には、先ず、ある日時の予測位置と、同じ日時の実績位置とを用いて観測位置を算出する。観測位置は、予測位置を地球上に投影した位置と、実績位置を地球上に投影した位置との中間位置として算出する。すなわち、予測位置と地球中心を結ぶ線分が地表面(標高0m)と交差する交点の座標と、実績位置と地球中心を結ぶ線分が地表面(標高0m)と交差する交点の座標を求める。そして、この2つの交点座標の中点の座標で表される位置を観測位置とする。いわゆる標高投影法を用いて観測位置を算出することになる。以上の観測位置の算出を、衛星予測暦に含まれている各日時について行う。
【0041】
観測位置を算出したら、衛星予測暦に含まれる各日時について、予測位置と観測位置間の距離(以下、「第1の観測距離」と称す。)、及び、実績位置と観測位置間の距離(以下、「第2の観測距離」と称す。)を算出し、第1の観測距離と第2の観測距離間の距離を「観測距離誤差」として求める。
【0042】
次いで、衛星予測暦に含まれる各日時について、時計予測誤差に光速を乗算することで、時計予測誤差の距離相当値(以下、「第1の時計誤差距離」と称す。)を算出する。同様に、衛星精密暦に記憶されている各日時について、時計実績誤差に光速を乗算することで、時計実績誤差の距離相当値(以下、「第2の時計誤差距離」と称す。)を算出する。そして、第1の時計誤差距離と第2の時計誤差距離との差を「時計距離誤差」として算出する。
【0043】
観測誤差は、観測距離誤差と時計距離誤差との合算値として算出することができる。このようにして、全ての衛星暦組合せについて、全てのGPS衛星の全ての期間における観測誤差を算出する。そして、図11に示すような観測誤差データ356を生成して、観測誤差データベース355に蓄積・記憶する。
【0044】
観測誤差データベース355は、衛星予測暦(衛星精密暦)の開始日時別に観測誤差データ356が蓄積記憶されたデータベースである。観測誤差データ356には、開始日時と対応付けて、第1期間〜第28期間までの各期間、及び、SV1〜SV32までの各GPS衛星における観測誤差が記憶されている。また、観測誤差データ356には、GPS衛星別に観測誤差の平均値を算出した衛星別平均観測誤差31(31−S1,31−S2,・・・,31−S32)と、期間別に観測誤差の平均値を算出した期間別平均観測誤差41(41−P1,41−P2,・・・,41−P28)とが記憶されている。
【0045】
(2)標準パターンの設定
次に、観測誤差を統計処理することで、予測軌道信頼度の標準的な変化パターンである「標準パターン」を設定する。例えば、N個の衛星暦組合せの中から、ランダムに1個の衛星暦組合せを選択して標準パターンを算出してもよいし、N個それぞれについて仮の標準パターンを作成して平均化することで最終的な標準パターンを算出することとしてもよい。以下では、ランダムに1個の衛星暦組合せを選択することとする。
【0046】
具体的には、観測誤差データベース355の中から、選択した衛星暦組合せの開始日時に対応する観測誤差データ356を抽出し、抽出した観測誤差データ356に記憶されている期間別平均観測誤差41を読み出す。
【0047】
図3は、期間別平均観測誤差をプロットしたグラフの一例である。図3において、横軸は日数(期間)、縦軸は期間別平均観測誤差をそれぞれ示している。期間は第1期間〜第28期間まであるが、本実施形態では各期間の長さを6時間としているため、第1期間〜第4期間までが1日目、第5期間〜第8期間までが2日目、・・・、第25期間〜第28期間までが7日目に相当する。
【0048】
衛星予測暦に記憶されている予測位置は、開始日時から将来のものであるほど、正確性が低くなる傾向がある。予測位置の正確性が低いと、予測位置及び観測位置から求めた観測距離(第1の観測距離)と、実績位置及び観測位置から求めた観測距離(第2の観測距離)との差(観測距離誤差)が大きくなるため、観測誤差は大きくなる。そのため、図3に示すように、将来の期間ほど、期間別平均観測誤差が大きくなる傾向がある。
【0049】
観測誤差が大きいほど、予測位置の正確性が低く、予測軌道の信頼性は低いと考えることができる。また、予測軌道信頼度は、値が大きいほど予測軌道の信頼性が低いことを示す指標値である。従って、期間別平均観測誤差が大きな期間ほど大きな標準値が設定されるように、標準パターンを設定する。
【0050】
具体的には、図8に示すような、期間別平均観測誤差と基礎値とが対応付けられたテーブルを予め用意しておく。そして、例えば日単位で、期間別平均観測誤差が含まれる数値範囲を特定し、特定した数値範囲に対応する基礎値を当該日にちに対応付けることで、標準パターンを設定する。図3では、標準パターンとして、1日目及び2日目(第1期間〜第8期間)に「3」、3日目及び4日目(第9期間〜第16期間)に「5」、5日目及び6日目(第17期間〜第24期間)に「7」、7日目(第25期間〜第28期間)に「8」がそれぞれ設定されている。
【0051】
(3)標準パターンの補正
次に、設定した標準パターンを補正することで、予測軌道信頼度を決定する。詳細には、図11の観測誤差データベース355に記憶されている観測誤差データ356を参照し、全てのGPS衛星及び全ての開始日時について算出した衛星別平均観測誤差31をGPS衛星間で相対評価することで、標準パターンの補正値を決定する。
【0052】
標準パターンの補正値は、例えば、相対的に衛星別平均観測誤差が小さいGPS衛星については「−1」、相対的に衛星別平均観測誤差が大きいGPS衛星については「+1」、それ以外のGPS衛星については「±0」とすることができる。勿論、より大きな値(例えば「+2」や「+3」)を補正値としてもよいし、より小さな値(例えば「−2」や「−3」)を補正値としてもよい。
【0053】
具体的に説明すると、先ず、全ての衛星別平均観測誤差の平均値と、衛星別平均観測誤差の標準偏差とを求める。そして、例えば、算出した平均値から標準偏差を減算した値を第1の閾値、算出した平均値に標準偏差を加算した値を第2の閾値とする。衛星別平均観測誤差が第1の閾値よりも小さくなる傾向があるGPS衛星は補正値を「−1」とし、第1の閾値以上第2の閾値以下となる傾向があるGPS衛星は補正値を「±0」とし、第2の閾値よりも大きくなる傾向があるGPS衛星は補正値を「+1」とする。尚、第1の閾値及び第2の閾値は、上述した手法以外にも適宜設定することが可能である。
【0054】
例えば、図4において、GPS衛星「SV1」は相対的に衛星別平均観測誤差が小さいと判定されたため補正値は「−1」とされ、また、GPS衛星「SV4」は相対的に衛星別平均観測誤差が大きいと判定されたため、補正値は「+1」とされている。また、GPS衛星「SV2」は衛星別平均観測誤差が大きくも小さくもないと判定されたため、補正値は「±0」とされている。
【0055】
尚、GPS衛星は、定期或いは不定期に軌道修正を行って、地球の周回軌道を補正している。軌道修正が行われると、衛星予測暦に含まれる予測位置のうち軌道修正が行われた後の予測位置が実績位置から大きくずれるため、軌道修正後の予測位置は、観測誤差が極端に大きな値となる。ここでは、この軌道修正等に起因して極端に大きくなった観測誤差のことを「異常観測誤差」と呼ぶ。
【0056】
本実施形態では、異常観測誤差が多数存在するGPS衛星(観測誤差が不良条件を満たすGPS衛星)については、補正値を「±0」とする。すなわち、異常観測誤差が多数存在するGPS衛星については、例外的に標準パターンの補正を行わないようにする。例えば、図4において、GPS衛星「SV3」は異常観測誤差が多数観測されたため、補正値が「±0」とされている。
【0057】
標準パターンの補正値が決まったら、当該補正値を用いて標準パターンを補正することで、予測軌道信頼度を決定する。具体的には、標準パターンを第1〜第28期間それぞれについて補正値を用いて補正し、補正後の値を、長期予測エフェメリスの第1予測対象期間〜第28予測対象期間の予測軌道信頼度とする。
【0058】
3.機能構成
図5は、サーバシステム3の機能構成を示すブロック図である。サーバシステム3は、CPU(Central Processing Unit)310と、操作部320と、通信部330と、ROM(Read Only Memory)340と、ハードディスク350と、RAM(Random Access Memory)360とを備え、各部がバス370で接続されたコンピュータシステムである。
【0059】
CPU310は、ROM340に記憶されているシステムプログラム等に従ってサーバシステム3の各部を統括的に制御するプロセッサである。本実施形態では、CPU310は、ROM340に記憶されている長期予測エフェメリス提供プログラム341に従って、携帯型電話機4に長期予測エフェメリスを提供する処理を行う。
【0060】
操作部320は、サーバシステム3の管理者による操作指示を受け付け、操作に応じた信号をCPU310に出力する入力装置である。この機能は、例えばキーボードやボタン、マウス等により実現される。
【0061】
通信部330は、システム内部で利用される各種データをインターネット等の通信ネットワークを介して外部システム2や携帯型電話機4とやり取りするための通信装置である。
【0062】
ROM340は、読み取り専用の不揮発性の記憶装置であり、CPU310がサーバシステム3を制御するためのシステムプログラムや、長期予測エフェメリスを携帯型電話機4に提供するためのプログラム、長期予測エフェメリスを生成するためのプログラム等の各種プログラムや各種データ等を記憶している。
【0063】
ハードディスク350は、磁気ヘッド等を用いてデータの読み書きを行う記憶装置であり、ROM340と同様、サーバシステム3が備える各種機能を実現するためのプログラムやデータ等を記憶している。
【0064】
RAM360は、読み書き可能な揮発性の記憶装置であり、CPU310により実行されるシステムプログラム、長期予測エフェメリス提供プログラム、各種処理プログラム、各種処理の処理中データ、処理結果などを一時的に記憶するワークエリアを形成している。
【0065】
4.データ構成
図6は、ROM340に格納されたデータの一例を示す図である。ROM340には、CPU310により読み出され、長期予測エフェメリス提供処理(図15参照)として実行される長期予測エフェメリス提供プログラム341と、観測誤差基礎値対応テーブル343とが記憶されている。また、長期予測エフェメリス提供プログラム341には、長期予測エフェメリス生成処理(図16及び図17参照)として実行される長期予測エフェメリス生成プログラム3411がサブルーチンとして含まれている。
【0066】
長期予測エフェメリス提供処理とは、CPU310が、長期予測エフェメリスデータ359を生成する処理を定期的に行い、携帯型電話機4から長期予測エフェメリスデータ359の要求信号を受信した場合に、生成しておいた長期予測エフェメリスデータ359を要求元の携帯型電話機4に送信する処理である。長期予測エフェメリス提供処理については、フローチャートを用いて詳細に後述する。
【0067】
長期予測エフェメリス生成処理とは、CPU310が、長期予測エフェメリスデータ359を生成する処理である。本実施形態では、CPU310は、長期予測エフェメリスデータ359を4時間に1回生成するものとする。長期予測エフェメリス生成処理についても、フローチャートを用いて詳細に後述する。
【0068】
図8は、観測誤差基礎値対応テーブル343のテーブル構成の一例を示す図である。観測誤差基礎値対応テーブル343には、期間別平均観測誤差が含まれる範囲である期間別平均観測誤差範囲3431と、期間別平均観測誤差が当該期間別平均観測誤差範囲3431に含まれる期間に設定する基礎値3433とが対応付けて記憶されている。例えば、期間別平均観測誤差が「20m〜40m」に含まれる期間は「5」が設定される。長期予測エフェメリス生成処理において、CPU310は、観測誤差基礎値対応テーブル343を用いて標準パターンを設定する。
【0069】
図7は、ハードディスク350に格納されたデータの一例を示す図である。ハードディスク350には、衛星予測暦データベース351と、衛星精密暦データベース353と、観測誤差データベース355と、標準パターンデータ357と、長期予測エフェメリスデータ359とが記憶されている。
【0070】
図9は、衛星予測暦データベース351のデータ構成の一例を示す図である。衛星予測暦データベース351は、複数の衛星予測暦352(352−1,352−2,352−3,・・・)が時系列に記憶されている。衛星予測暦352は、各GPS衛星SVの1週間後までの予測位置及び時計予測誤差が15分おきに記憶された離散的なデータであり、開始日時毎に纏められたデータである。便宜的に、1週間の期間を6時間毎にグループ分けすることで、第1期間〜第28期間の28個の期間が構成される。
【0071】
例えば、衛星予測暦352−1は、開始日時が「2008年8月1日0時00分」のデータである。そして、「2008年8月1日5時45分」におけるGPS衛星「SV2」の予測位置は「(Xp32,Yp32,Zp32)」であり、原子時計の予測誤差は「Δtp32」である。
【0072】
CPU310は、外部システム2から定期的に(例えば4時間に1回)衛星予測暦を受信する。そして、受信した衛星予測暦を衛星予測暦データベース351に蓄積・記憶するためにデータフォーマットの加工処理等を行う。具体的には、開始日時が異なり、長期予測エフェメリスの生成対象期間(例えば1週間)と同じ期間分のデータが格納された複数の衛星予測暦352を生成して、衛星予測暦データベース351に蓄積・記憶させる。
【0073】
図10は、衛星精密暦データベース353のデータ構成の一例を示す図である。衛星精密暦データベース353には、複数の衛星精密暦354(354−1,354−2,354−3,・・・)が時系列に記憶されている。衛星精密暦354は、各GPS衛星SVの1週間分の実績位置及び時計実績誤差が15分おきに記憶された離散的なデータであり、開始日時毎に纏められたデータである。1週間の期間を6時間毎にグループ分けすることで、第1期間〜第28期間の28個の期間が構成される。
【0074】
例えば、衛星精密暦354−1は、開始日時が「2008年8月1日0時00分」のデータである。そして、「2008年8月1日5時45分」におけるGPS衛星「SV2」の実績位置は「(Xm32,Ym32,Zm32)」であり、原子時計の実績誤差は「Δtm32」であった。
【0075】
CPU310は、外部システム2から定期的に(例えば4時間に1回)衛星精密暦を受信する。そして、受信した衛星精密暦を衛星精密暦データベース353に蓄積・記憶するためにデータフォーマットの加工処理等を行う。具体的には、衛星予測暦データベース351に格納されている複数の衛星予測暦352と開始日時が対応する衛星精密暦354を生成して、衛星精密暦データベース353に蓄積・記憶させる。
【0076】
図11は、観測誤差データベース355のデータ構成の一例を示す図である。観測誤差データベース355は、複数の観測誤差データ356(356−1,356−2,356−3,・・・)が開始日時別に記憶されたデータベースである。各観測誤差データ356には、各期間及び各GPS衛星における観測誤差が記憶されている。また、各衛星について全ての期間の観測誤差を平均した衛星別平均観測誤差31(31−S1,31−S2,・・・,31−S32)と、各期間について全てのGPS衛星の観測誤差を平均した期間別平均観測誤差41(41−P1,41−P2,・・・,41−P28)とが記憶されている。
【0077】
例えば、観測誤差データ356−1は、開始日時が「2008年8月1日0時00分」のデータである。また、GPS衛星「SV2」の「第2期間」における観測誤差は「E22」である。また、GPS衛星「SV2」の衛星別平均観測誤差は「ES2」であり、「第2期間」の期間別平均観測誤差は「EP2」である。
【0078】
CPU310は、長期予測エフェメリス提供処理において、上述した原理に従って各GPS衛星の各期間について観測誤差を算出する。また、算出した観測誤差を用いて、衛星別平均観測誤差31及び期間別平均観測誤差41を算出する。そして、これらを開始日時と対応付けた観測誤差データ356を生成して、観測誤差データベース355に蓄積・記憶させる。
【0079】
図12は、標準パターンデータ357のデータ構成の一例を示す図である。標準パターンデータ357には、期間3571と、標準パターン3573とが対応付けて記憶されている。例えば、標準パターンとして、第28期間には「8」が設定されている。
【0080】
CPU310は、長期予測エフェメリス生成処理において、上述した原理に従って標準パターンを設定する。そして、設定した標準パターン3573を期間3571と対応付けて、標準パターンデータ357に記憶させる。
【0081】
図13は、長期予測エフェメリスデータ359のデータ構成の一例を示す図である。長期予測エフェメリスデータ359には、長期予測エフェメリスデータの生成日時3591と、GPS衛星SV1〜SV32の予測エフェメリス3593(3593−1〜3593−32)とが対応付けて記憶されている。
【0082】
図14は、予測エフェメリス3593のデータ構成の一例を示す図である。予測エフェメリス3593(3593−1,3593−2,・・・,3593−32)には、各予測対象期間それぞれについて、軌道長半径や離心率、軌道傾斜角といったケプラーの衛星軌道パラメータの値と、衛星時計の基準時刻、衛星時計のオフセット、衛星時計のドリフト及び衛星時計周波数のドリフトでなるクロック補正パラメータの値と、信頼性パラメータである予測軌道信頼度とが対応付けて記憶されている。
【0083】
CPU310は、長期予測エフェメリス生成処理において、各GPS衛星SVについて、予測対象期間毎に衛星軌道パラメータ、クロック補正パラメータ及び信頼性パラメータの値を算出して予測エフェメリス3593を生成する。そして、全てのGPS衛星SVについて生成した予測エフェメリス3593を纏めて、生成日時3591と対応付けて長期予測エフェメリス359を生成して、ハードディスク350に記憶させる。
【0084】
5.処理の流れ
図15は、ROM340に記憶されている長期予測エフェメリス提供プログラム341がCPU310により読み出されて実行されることで、サーバシステム3において実行される長期予測エフェメリス提供処理の流れを示すフローチャートである。
【0085】
先ず、CPU310は、外部システム2から衛星予測暦及び衛星精密暦を受信したか否かを判定し(ステップA1)、受信しなかったと判定した場合は(ステップA1;No)、ステップA9へと処理を移行する。
【0086】
また、受信したと判定した場合は(ステップA1;Yes)、CPU310は、受信した衛星予測暦及び衛星精密暦を加工処理し、開始日時及び期間が同一である複数の衛星予測暦352及び衛星精密暦354を生成する(ステップA3)。そして、生成した衛星予測暦352及び衛星精密暦354を、ハードディスク350の衛星予測暦データベース351及び衛星精密暦データベース353にそれぞれ蓄積・記憶させる(ステップA5)。
【0087】
次いで、CPU310は、観測誤差算出処理を行う(ステップA7)。具体的には、上述した原理に従って観測誤差を算出し、観測誤差データ356を生成する。そして、生成した観測誤差データ356を観測誤差データベース355に蓄積・記憶させる。
【0088】
次いで、CPU310は、長期予測エフェメリスの生成日時であるか否かを判定する(ステップA9)。本実施形態では、4時間に1回長期予測エフェメリスを生成するものとする。そして、まだ生成日時ではないと判定した場合は(ステップA9;No)、CPU310は、ステップA13へと処理を移行する。
【0089】
また、長期予測エフェメリスの生成日時であると判定した場合は(ステップA9;Yes)、CPU310は、ROM340に記憶されている長期予測エフェメリス生成プログラム3411を読み出して実行することで、長期予測エフェメリス生成処理を行う(ステップA11)。
【0090】
図16及び図17は、長期予測エフェメリス生成処理の流れを示すフローチャートである。
先ず、CPU310は、標準パターン設定処理を行う(ステップB1)。具体的には、観測誤差データベース355に記憶されている観測誤差データ356を用いて、上述した原理に従って、各期間3571について標準パターン3573を設定する。そして、期間3571と標準パターン3573とを対応付けた標準パターンデータ357を、ハードディスク350に記憶させる。
【0091】
次いで、CPU310は、衛星別平均観測誤差をGPS衛星間で相対評価して、各GPS衛星の標準パターンの補正値を決定する(ステップB3)。そして、決定した補正値で標準パターンを補正することで、各期間の信頼性パラメータの値である予測軌道信頼度を決定する(ステップB5)。
【0092】
その後、CPU310は、長期予測エフェメリスの今回の生成日時(現在日時)に基づいて各予測対象期間を判定する(ステップB7)。すなわち、今回の生成日時から1週間後までの期間を生成対象期間とし、当該生成対象期間を6時間毎に区切った各期間を予測対象期間と判定する。
【0093】
次いで、CPU310は、各GPS衛星SVそれぞれについて、ループAの処理を実行する(ステップB9〜B25)。ループAの処理では、CPU310は、ステップB7で判定した各予測対象期間それぞれについて、ループBの処理を実行する(ステップB11〜B21)。
【0094】
ループBの処理では、CPU310は、ハードディスク350の衛星予測暦データベース351に記憶されている最新の衛星予測暦352から、当該GPS衛星SVの当該予測対象期間の各日時における予測位置を読み出す(ステップB13)。
【0095】
そして、CPU310は、読み出した予測位置を用いて、ケプラーの楕円軌道モデルに従って当該GPS衛星SVの当該予測対象期間における予測軌道を算出し、ケプラーの衛星軌道パラメータの値を求める(ステップB15)。尚、予測軌道の具体的な算出方法については公知であるため、詳細な説明を省略する。
【0096】
その後、CPU310は、最新の衛星予測暦352から、当該GPS衛星SVの当該予測対象期間の各日時における時計予測誤差を読み出す(ステップB17)。そして、CPU310は、読み出した時計予測誤差を用いて、当該GPS衛星の当該予測対象期間のクロック補正パラメータの値を求める(ステップB19)。
【0097】
時刻「t」における時計予測誤差「Δt」は、クロック補正パラメータである衛星時計の基準時刻「tc」、衛星時計のオフセット「a0」、衛星時計のドリフト「a1」及び衛星時計周波数のドリフト「a2」を用いて、次式(1)によって近似することができる。
Δt=a0+a1(t−tc)+a2(t−tc)2 ・・・(1)
【0098】
式(1)は、時計予測誤差の経時変化を近似するための時計誤差モデル式である。衛星予測暦351に含まれる各日時の時計予測誤差「Δt」をサンプリングデータとして、例えば最小二乗法を用いた近似計算を行うことで、クロック補正パラメータの値を算出することができる。その後、CPU310は、次の予測対象期間へと処理を移行する。
【0099】
全ての予測対象期間についてステップB13〜B19の処理を行った後、CPU310は、ループBの処理を終了する(ステップB21)。その後、CPU310は、全ての予測対象期間についてステップB15で求めた衛星軌道パラメータの値と、ステップB19で求めたクロック補正パラメータの値と、ステップB5で決定した信頼性パラメータの値(予測軌道信頼度)とを纏めて当該GPS衛星SVの予測エフェメリス3593を生成する(ステップB23)。そして、CPU310は、次のGPS衛星SVへと処理を移行する。
【0100】
全てのGPS衛星SVについてステップB11〜B23の処理を行った後、CPU310は、ループAの処理を終了する(ステップB25)。その後、CPU310は、全てのGPS衛星SVについてステップB23で生成した予測エフェメリス3593を纏めて、生成日時3591と対応付けて長期予測エフェメリスデータ359を生成し、ハードディスク350に記憶させる(ステップB27)。そして、CPU310は、長期予測エフェメリス生成処理を終了する。
【0101】
図15の長期予測エフェメリス提供処理に戻って、長期予測エフェメリス生成処理を行った後、CPU310は、携帯型電話機4から長期予測エフェメリスデータ359の要求信号を受信したか否かを判定する(ステップA13)。そして、受信しなかったと判定した場合は(ステップA13;No)、ステップA1に戻る。
【0102】
また、要求信号を受信したと判定した場合は(ステップA13;Yes)、CPU310は、ハードディスク350に記憶されている長期予測エフェメリスデータ359を要求元の携帯型電話機4に送信する(ステップA15)。そして、CPU310は、ステップA1に戻る。
【0103】
6.作用効果
本実施形態によれば、測位システム1において、サーバシステム3は、外部システム2から受信したGPS衛星SVの位置を時系列に予測した衛星予測暦に含まれる各予測位置を、測位用衛星の実際の位置を時系列に記憶した衛星精密暦に含まれる実績位置と比較することで、衛星予測暦に含まれる予測位置の正確性の変化を分析する。そして、連続する複数の予測対象期間それぞれにおける予測衛星軌道の衛星軌道パラメータの値を纏めてなる長期予測エフェメリスの予測対象期間それぞれの信頼性を、先の分析結果に基づいて判定する。そして、信頼性の判定結果を長期予測エフェメリスに含めて、携帯型電話機4に提供する。
【0104】
より具体的には、予測位置と観測位置間の距離と、実績位置と観測位置間の距離との差で表される観測距離誤差を算出する。また、時計予測誤差の距離相当値と時計実績誤差の距離相当値との差で表される時計距離誤差を算出する。次いで、観測距離誤差と時計距離誤差との合算値を観測誤差として算出する。そして、当該観測誤差を期間別に平均した期間別平均観測誤差を算出し、当該期間別平均観測誤差の時間変化に基づいて、標準パターンを設定する。
【0105】
その後、観測誤差をGPS衛星別に平均した衛星別平均観測誤差を算出し、この衛星別平均観測誤差をGPS衛星間で相対評価することで、標準パターンの補正値を決定する。そして、決定した補正値を用いて標準パターンを補正することで、各予測対象期間それぞれの予測軌道の信頼性の指標値である予測軌道信頼度を決定する。
【0106】
GPS衛星SVの実績位置は、GPS衛星SVのその日時における実際の位置である。そのため、予測位置及び実績位置と用いて算出した観測誤差が大きい場合は、予測位置の正確性が低く、この予測位置に基づいて予測した予測衛星軌道の信頼性は低いと言える。そのため、期間別平均観測誤差が大きい期間ほど信頼性が低くなるように標準パターンを設定した上で、GPS衛星間で衛星別平均観測誤差が相対的に小さい衛星は信頼性が高くなるように標準パターンを補正し、逆に、衛星別平均観測誤差が相対的に大きい衛星は信頼性が低くなるように標準パターンを補正する。これにより、長期予測エフェメリスの信頼性を適確に判定し、適切な信頼性パラメータの値を携帯型電話機4に提供することが可能となる。
【0107】
7.変形例
7−1.測位システム
上述した実施形態では、サーバシステム3と携帯型電話機4を備えた測位システム1を例に挙げて説明したが、本発明を適用可能な測位システムはこれに限られるわけではない。例えば、携帯型電話機4の代わりに、測位装置を備えたノート型パソコンやPDA(Personal Digital Assistant)、カーナビゲーション装置等の電子機器に適用することも可能である。
【0108】
また、上述した実施形態では、情報提供装置の一種としてサーバシステム3を例に挙げて説明したが、情報提供装置はサーバシステム3に限られるわけではない。例えば、汎用のパソコン等であってもよい。
【0109】
7−2.衛星測位システム
また、上述した実施形態では、衛星測位システムとしてGPSを例に挙げて説明したが、WAAS(Wide Area Augmentation System)、QZSS(Quasi Zenith Satellite System)、GLONASS(GLObal NAvigation Satellite System)、GALILEO等の他の衛星測位システムであってもよい。
【0110】
7−3.観測誤差の時系列変化に基づく標準パターンの補正
衛星別観測誤差を相対評価することに代えて又は加えて、衛星別観測誤差の開始日時の時系列変化に基づいて標準パターンを補正することも可能である。
【0111】
図18は、この場合にCPU310が行う処理である第2の長期予測エフェメリス生成処理のうち、図16の長期予測エフェメリス生成処理に対応する部分を示すフローチャートである。尚、長期予測エフェメリス生成処理と同一のステップについては同一の符号を付して説明を省略し、長期予測エフェメリス生成処理とは異なる部分について説明する。
【0112】
第2長期予測エフェメリス生成処理では、CPU310は、ステップB3において標準パターンの補正値を決定した後、各GPS衛星について、開始日時の時系列変化に基づく衛星別平均観測誤差の増減を判定し、標準パターンの補正値を修正する(ステップC4)。
【0113】
具体的には、開始日時が生成日時(現在日時)に近づくにつれて衛星別平均観測誤差が増加する傾向がある場合は、補正値に「1」を加算し、衛星別平均観測誤差は減少する傾向がある場合は、補正値から「1」を減算する。また、どちらとも判断がつかない場合は、補正値をそのままとする。勿論、「1」よりも大きな値(例えば「2」や「3」)を補正値に加算/減算することとしてもよい。衛星別平均観測誤差の増減は、微分値を算出するなどして判断することが可能である。
【0114】
その後、CPU310は、修正した補正値で標準パターンを補正することで各期間の予測軌道信頼度を決定する(ステップC5)。そして、CPU310は、ステップB7以降へと処理を移行する。
【0115】
7−4.衛星位置誤差
上述した実施形態では、観測誤差を用いて予測軌道信頼度を決定するものとして説明したが、観測誤差の代わりに、衛星位置の誤差(以下、「衛星位置誤差」と称す。)を用いて予測軌道信頼度を決定することとしてもよい。衛星位置誤差は、衛星予測暦に含まれる予測位置と、衛星精密暦に含まれる実績位置間の距離として求めることができる。観測誤差を衛星位置誤差に置き換えることで、上述した原理に従って、予測軌道信頼度を決定することができる。
【0116】
7−5.標準パターンの設定
上述した実施形態では、図3に示したように、標準パターンを日単位で設定する場合を一例として説明したが、標準パターンを期間単位で設定してもよいことは勿論である。すなわち、図8の観測誤差基礎値対応テーブル343に従って、各期間それぞれについて、期間別平均観測誤差が何れの期間別平均観測誤差範囲3431に含まれるかを判定する。そして、判定した期間別平均観測誤差範囲3431に対応する基礎値3433を読み出して各期間における基礎値を設定することで、標準パターンを設定する。
【0117】
7−6.平均観測誤差
上述した実施形態では、期間別平均観測誤差に基づいて標準パターンを設定するものとして説明したが、期間別に算出した観測誤差の最大値(期間別最大観測誤差)に基づいて標準パターンを設定することにしてもよい。また、衛星別平均観測誤差を相対評価して標準パターンの補正値を決定するのではなく、衛星別に算出した観測誤差の最大値(衛星別最大観測誤差)を相対評価して標準パターンの補正値を決定することにしてもよい。
【0118】
7−7.長期予測エフェメリスの生成
上述した実施形態では、サーバシステム3が長期予測エフェメリスデータを生成して携帯型電話機4に提供するものとして説明したが、携帯型電話機4が自ら長期予測エフェメリスデータを生成することとしてもよい。すなわち、携帯型電話機4が外部システム2から定期的に衛星予測暦及び衛星精密暦を取得し、取得した衛星予測暦及び衛星精密暦を用いて長期予測エフェメリス生成処理を実行することで、長期予測エフェメリスデータを生成する。携帯型電話機4の代わりに、測位装置を備えたノート型パソコンやPDA、カーナビゲーション装置等の電子機器に適用する場合も同様である。
【0119】
また、上述した実施形態では、サーバシステム3が予め所定の時間間隔(例えば4時間に1回)で長期予測エフェメリスデータを生成しておき、携帯型電話機4から長期予測エフェメリスデータの要求を受けた場合に、生成しておいた長期予測エフェメリスデータを送信するものとして説明した。このような構成ではなく、サーバシステム3が、携帯型電話機4から長期予測エフェメリスデータの要求を受けたときに長期予測エフェメリスデータを生成して、携帯型電話機4に送信する構成としてもよい。
【0120】
7−8.生成対象期間
上述した実施形態では、長期予測エフェメリスの生成日時を基準として1週間後までの期間を生成対象期間として長期予測エフェメリスを生成するものとして説明したが、生成対象期間は1週間よりも長い期間(例えば2週間)としてもよいし、1週間よりも短い期間(例えば3日)としてもよい。GPS衛星SVから送信される航法データとしてのエフェメリスは有効期間が一般に4時間程度であるが、長期予測エフェメリスは少なくともGPS衛星SVから送信される航法データとしてのエフェメリスよりも有効期間が長ければよい。1日以上であれば好適である。
【0121】
7−9.予測対象期間
また、上述した実施形態では、予測対象期間の長さを6時間とするものとして説明したが、これに限定されるわけではなく、4時間や8時間等としてもよく、適宜設定可能であることは勿論である。
【図面の簡単な説明】
【0122】
【図1】測位システムの概略構成を示す図。
【図2】衛星予測暦及び衛星精密暦の説明図。
【図3】期間別平均観測誤差をプロットしたグラフ。
【図4】標準パターンの補正値決定の説明図。
【図5】サーバシステムの機能構成を示すブロック図。
【図6】サーバシステムのROMに格納されたデータの一例を示す図。
【図7】サーバシステムのハードディスクに格納されたデータの一例を示す図。
【図8】観測誤差基礎値対応テーブルのテーブル構成の一例を示す図。
【図9】衛星予測暦データベースのデータ構成の一例を示す図。
【図10】衛星精密暦データベースのデータ構成の一例を示す図。
【図11】観測誤差データベースのデータ構成の一例を示す図。
【図12】標準パターンデータのデータ構成の一例を示す図。
【図13】長期予測エフェメリスデータのデータ構成の一例を示す図。
【図14】予測エフェメリスのデータ構成の一例を示す図。
【図15】長期予測エフェメリス提供処理の流れを示すフローチャート。
【図16】長期予測エフェメリス生成処理の流れを示すフローチャート。
【図17】長期予測エフェメリス生成処理の流れを示すフローチャート。
【図18】第2の長期予測エフェメリス生成処理の流れを示すフローチャート。
【符号の説明】
【0123】
1 測位システム 、 2 外部システム、 3 サーバシステム、
4 携帯型電話機、 310 CPU、 320 操作部、 330 通信部、
340 ROM、 350 ハードディスク、 360 RAM、 370 バス、
SV GPS衛星
【技術分野】
【0001】
本発明は、長期予測軌道データの信頼性判定方法、長期予測軌道データの提供方法及び情報提供装置に関する。
【背景技術】
【0002】
測位用信号を利用した測位システムとしては、GPS(Global Positioning System)が広く知られており、携帯型電話機やカーナビゲーション装置等に内蔵された測位装置に利用されている。GPSでは、複数のGPS衛星の位置や各GPS衛星から自機までの擬似距離等の情報に基づいて自機の位置を示す3次元の座標値と時計誤差とを求める測位演算を行う。
【0003】
GPSによる測位では、先ず、GPS衛星から発信されるGPS衛星信号に重畳されているアルマナックやエフェメリスといった航法データに基づいて、GPS衛星の位置、速度、移動方向等の衛星情報を算出する。アルマナックは衛星を捕捉する際の有力な手掛かりとはなるが、衛星情報の精度が悪いため測位演算に使用しないのが一般的である。一方、エフェメリスは、衛星情報の精度が良いため、衛星を捕捉する際の有力な手掛かりとなるだけでなく、測位演算にも使用することができる。従って、例えばエフェメリスを保持していない状態で測位を開始した場合には、エフェメリスをGPS衛星信号から取得しなければならず、初回測位時間(TTFF:Time To First Fix)が増大する。
【0004】
そこで、サーバクライアントシステムであって、1週間といった長期間分のエフェメリスである長期予測エフェメリス(長期予測軌道データ)をサーバが予測し、クライアントである測位装置に提供する技術が開発されており、例えば特許文献1や特許文献2に開示されている。
【特許文献1】米国特許出願公開第2002/0188403号明細書
【特許文献2】米国特許出願公開第2005/0212700号明細書
【発明の開示】
【発明が解決しようとする課題】
【0005】
長期予測エフェメリスの定義方法としては、通常のエフェメリスと同様のデータフォーマットで定義する方法が考えられる。すなわち、衛星軌道の近似モデルの1つであるケプラーの楕円軌道モデルを用いて衛星軌道を近似し、その時のモデル式のパラメータ(以下、「衛星軌道パラメータ」と称す。)の値によって、長期予測エフェメリスを定義する方法である。測位用衛星の将来の位置を所定時間おきに時系列に予測した予測位置でなる衛星予測暦(予測位置データ)は、所定の商用システムから提供されている。ケプラーの楕円軌道モデルによる近似計算は、この衛星予測暦を用いて行うことが可能である。
【0006】
しかし、衛星予測暦に含まれる測位用衛星の予測位置は、将来になるほど、測位用衛星の実際位置からずれる傾向がある。そのため、ケプラーの楕円軌道モデルによる近似計算を行って長期予測エフェメリスを生成する場合に、近似計算により求めた衛星軌道は、生成日時から将来のものであるほど、実際の衛星軌道からずれたものとなる可能性がある。従来は、長期予測エフェメリスを一旦生成した後は、生成した長期予測エフェメリスの信頼性(測位に使用するのに適しているか否か)を判定する手法が存在しなかったため、サーバから長期予測エフェメリスを取得した測位装置は、信頼性の低い長期予測エフェメリスを使用して測位を行う場合があり得た。
【0007】
本発明は、上述した課題に鑑みて為されたものであり、長期予測エフェメリスの信頼性を判定するための手法を提案することを目的としている。
【課題を解決するための手段】
【0008】
以上の課題を解決するための第1の発明は、測位用衛星の位置を時系列に予測した予測位置データに含まれる各予測位置を、当該予測位置に対応する前記測位用衛星の実績位置と比較することで、前記予測位置データに含まれる時系列の予測位置の正確性の変化を分析することと、連続する複数の予測対象期間それぞれにおける予測衛星軌道を纏めてなる長期予測軌道データの前記予測対象期間それぞれの信頼性を、前記分析の結果に基づいて判定することと、を含む長期予測軌道データの信頼性判定方法である。
【0009】
また、他の発明として、連続する複数の予測対象期間それぞれにおける予測衛星軌道を纏めてなる長期予測軌道データを生成する生成部と、測位用衛星の位置を時系列に予測した予測位置データに含まれる各予測位置を、当該予測位置に対応する前記測位用衛星の実績位置と比較することで、前記予測位置データに含まれる時系列の予測位置の正確性の変化を分析する分析部と、前記生成部により生成された長期予測軌道データの前記予測対象期間それぞれの信頼性を、前記分析部の分析結果に基づいて判定する判定部と、前記生成部により生成された長期予測軌道データと前記判定部による判定結果とを測位装置に提供する提供部と、を備えた情報提供装置を構成してもよい。
【0010】
この第1の発明等によれば、測位用衛星の位置を時系列に予測した予測位置データに含まれる各予測位置を、当該予測位置に対応する測位用衛星の実績位置と比較することで、予測位置データに含まれる時系列の予測位置の正確性の変化を分析する。そして、連続する複数の予測対象期間それぞれにおける予測衛星軌道を纏めてなる長期予測軌道データの予測対象期間それぞれの信頼性を、分析結果に基づいて判定する。
【0011】
測位用衛星の実績位置は、測位用衛星の実際にいた位置である。そのため、測位用衛星の予測位置と実績位置との差が大きい場合は、予測位置の正確性が低く、この予測位置に基づいて予測した予測衛星軌道の信頼性は低いと言える。従って、予測位置の正確性の変化の分析結果に基づくことで、長期予測軌道データの信頼性を適切に判定することが可能となる。
【0012】
また、第2の発明として、第1の発明の長期予測軌道データの信頼性判定方法であって、前記分析することは、予測した日時が異なる複数の前記予測位置データそれぞれについて、当該予測位置データに含まれる各予測位置を対応する実績位置と比較することと、前記比較の結果を統計処理して、前記測位用衛星の位置を時系列に予測した場合の時間経過に対する予測位置の正確性の変化パターンを算出することと、を含み、前記信頼性を判定することは、前記変化パターンに基づいて、前記長期予測軌道データの前記予測対象期間それぞれの信頼性を判定することである長期予測軌道データの信頼性判定方法を構成してもよい。
【0013】
この第2の発明によれば、予測した日時が異なる複数の予測位置データそれぞれについて、各予測位置を対応する実績位置と比較する。そして、比較の結果を統計処理して予測位置の正確性の変化パターンを算出する。この変化パターンは、いわば予測位置の正確性の変化の様子を表す標準パターンとも言える。そして、算出した変化パターンに基づいて、長期予測軌道データの予測対象期間それぞれの信頼性を判定する。
【0014】
また、第3の発明として、第2の発明の長期予測軌道データの信頼性判定方法であって、前記予測位置データには、複数の前記測位用衛星それぞれの予測位置が含まれ、前記比較することは、前記測位用衛星それぞれの各予測位置を実績位置と比較することであり、前記変化パターンを算出することは、前記測位用衛星それぞれの前記比較の結果を統計処理して、前記測位用衛星間で標準的な前記変化パターンを算出することである長期予測軌道データの信頼性判定方法を構成してもよい。
【0015】
この第3の発明によれば、測位用衛星それぞれの各予測位置を実績位置と比較し、測位用衛星それぞれの比較の結果を統計処理して、測位用衛星間で標準的な前記変化パターンを算出する。
【0016】
また、第4の発明として、第3の発明の長期予測軌道データの信頼性判定方法であって、前記長期予測軌道データは、前記予測対象期間それぞれについての各測位用衛星毎のデータでなり、前記測位用衛星それぞれの前記比較の結果を前記測位用衛星間で相対評価することを更に含み、前記信頼性を判定することは、前記相対評価の結果に基づいて前記変化パターンを前記測位用衛星毎に補正して、前記長期予測軌道データの前記予測対象期間それぞれの信頼性を各測位用衛星毎に判定することを含む長期予測軌道データの信頼性判定方法を構成してもよい。
【0017】
この第4の発明によれば、測位用衛星それぞれの比較の結果を測位用衛星間で相対評価し、相対評価の結果に基づいて変化パターンを測位用衛星毎に補正して、長期予測軌道データの予測対象期間それぞれの信頼性を各測位用衛星毎に判定する。
【0018】
予測位置と実績位置との比較結果を測位用衛星間で相対評価した結果、相対的に比較結果が良いと判定した測位用衛星については、長期予測軌道データの信頼性が高くなるように変化パターンを補正する。逆に、相対的に比較結果が悪いと判定した測位用衛星については、長期予測軌道データの信頼性が低くなるように変化パターンを補正する。これにより、測位用衛星間の比較結果を考慮に入れた適切な信頼性の判定を実現し得る。
【0019】
また、第5の発明として、第3又は第4の発明の長期予測軌道データの信頼性判定方法であって、前記分析することは、前記測位用衛星それぞれについて、前記予測位置データの予測日時毎に、当該予測位置データに関する前記比較の結果を纏めることを含み、前記信頼性を判定することは、前記予測位置データの予測日時毎に纏めた前記比較の結果の前記予測日時の時系列変化に基づいて前記変化パターンを前記測位用衛星毎に補正して、前記長期予測軌道データの前記予測対象期間それぞれの信頼性を各測位用衛星毎に判定することを含む長期予測軌道データの信頼性判定方法を構成してもよい。
【0020】
この第5の発明によれば、測位用衛星それぞれについて、予測位置データの予測日時毎に、当該予測位置データに関する比較の結果を纏める。そして、予測位置データの予測日時毎に纏めた比較の結果の予測日時の時系列変化に基づいて変化パターンを測位用衛星毎に補正して、長期予測軌道データの予測対象期間それぞれの信頼性を各測位用衛星毎に判定する。
【0021】
予測位置と実績位置との比較結果が時間経過に伴って良くなる測位用衛星については、長期予測軌道データの信頼性が高くなるように変化パターンを補正する。逆に、比較結果が時間経過に伴って悪くなる測位用衛星については、長期予測軌道データの信頼性が低くなるように変化パターンを補正する。これにより、比較結果の時間変化を考慮に入れた適切な信頼性の判定を実現し得る。
【0022】
また、第6の発明として、第5の発明の長期予測軌道データの信頼性判定方法であって、前記信頼性を判定することは、前記複数の測位用衛星のうち、当該測位用衛星の予測位置を実績位置と比較した前記比較の結果が所定の不良条件を満たさない測位用衛星についてのみ前記変化パターンの前記補正を行った上で信頼性の判定を行うことを含む長期予測軌道データの信頼性判定方法を構成してもよい。
【0023】
この第6の発明によれば、複数の測位用衛星のうち、当該測位用衛星の予測位置を実績位置と比較した比較の結果が所定の不良条件を満たさない測位用衛星についてのみ変化パターンの補正を行った上で信頼性の判定を行う。予測位置と実績位置との比較結果が所定の不良条件を満たす測位用衛星については、予測位置が実績位置から大きくずれる何らかの原因(例えば衛星の軌道修正)が発生したものと推定して、変化パターンを補正しないことにしたものである。
【0024】
また、第7の発明として、長期予測軌道データを生成することと、第1〜第6の何れかの発明の長期予測軌道データの信頼性判定方法を用いて、前記生成した長期予測軌道データの信頼性を判定することと、前記生成した長期予測軌道データと前記判定した判定結果とを測位装置に提供することと、を含む長期予測軌道データの提供方法を構成してもよい。
【0025】
この第7の発明によれば、長期予測軌道データを生成し、長期予測軌道データの信頼性判定方法を用いて、生成した長期予測軌道データの信頼性を判定する。そして、生成した長期予測軌道データと信頼性の判定結果とを測位装置に提供する。これにより、測位装置は、提供された長期予測軌道データを用いて測位を行うが、この際、信頼性の低いデータを測位に使用しないようにすることができる。
【発明を実施するための最良の形態】
【0026】
以下、図面を参照して、本発明に好適な実施形態の一例を説明する。但し、本発明を適用可能な実施形態がこれに限定されるわけではない。
【0027】
1.システム構成
図1は、本実施形態における測位システム1の概略構成を示す図である。測位システム1は、外部システム2と、情報提供装置の一種であるサーバシステム3と、測位装置を備えた電子機器の一種である携帯型電話機4と、測位用衛星の一種である複数のGPS衛星SV(SV1,SV2,SV3,SV4,・・・)とを備えて構成される。尚、携帯型電話機4が必要なデータをサーバシステム3から取得した後は、携帯型電話機4とGPS衛星SVとで測位が可能であるため、携帯型電話機4とGPS衛星SVとで1つの測位システムが構成されるということもできる。また、地上側のシステムとして、サーバシステム3と携帯型電話機4とで測位システムと呼ぶこともできる。
【0028】
外部システム2は、GPS衛星SVから定期的に衛星信号を受信し、当該衛星信号に含まれる航法データ等に基づいて衛星予測暦を生成してサーバシステム3に提供する公知のシステムである。外部システム2が提供する衛星予測暦は、各GPS衛星SVそれぞれについて、将来の位置を予測した予測位置及びGPS衛星SVに搭載された原子時計の誤差を予測した時計予測誤差を所定時間おき(例えば15分おき)に時系列に並べた位置のデータである。
【0029】
また、外部システム2は、将来のデータとしての衛星予測暦を提供する他に、過去の事実のデータも提供する。すなわち、外部システム2は、GPS衛星SVの実際の位置である実績位置及びGPS衛星SVに搭載された原子時計の実際の誤差である時計実績誤差を含む衛星精密暦を過去の事実のデータとして生成して、サーバシステム3に提供する。実績位置及び時計実績誤差の算出方法については公知であるため、詳細な説明を省略する。外部システム2は、例えば衛星予測暦や衛星精密暦の提供を業務とする民間や公営の団体のコンピュータシステムに相当する。
【0030】
サーバシステム3は、衛星予測暦及び衛星精密暦を外部システム2から取得し、当該衛星予測暦及び衛星精密暦を用いて、全てのGPS衛星SVの予測されるエフェメリスであって、少なくとも1日以上の例えば1週間といった長期間有効なエフェメリス(以下、本実施形態において「長期予測エフェメリス」と称す。長期間有効な軌道でもあるため、長期予測軌道データとも言える。)を生成・提供するサーバを備えたシステムである。
【0031】
携帯型電話機4は、ユーザが通話やメールの送受信等を行うための電子機器であり、通話やメールの送受信といった携帯型電話機としての本来の機能の他、位置を計測する機能(測位機能)をなす測位装置を具備している。携帯型電話機4は、ユーザ操作に従って、サーバシステム3に対して長期予測エフェメリスの要求信号を送信し、サーバシステム3から長期予測エフェメリスを受信する。そして、受信した長期予測エフェメリスを用いてGPS衛星SVを捕捉し、衛星信号に基づく測位演算を実行する。
【0032】
2.原理
サーバシステム3は、外部システム2から取得した衛星予測暦を用いて、長期予測エフェメリスを生成する処理を行う。具体的には、長期予測エフェメリスの生成日時を基準として1週間後までの期間を「生成対象期間」とし、当該生成対象期間を衛星軌道の近似・モデル化を行う複数の期間(以下、「予測対象期間」と称す。)に区切る。本実施形態では、予測対象期間の長さを一律に6時間とする。すなわち、1週間の生成対象期間を6時間毎に28個の予測対象期間(第1予測対象期間〜第28予測対象期間)に区切る。
【0033】
そして、サーバシステム3は、外部システム2から取得した衛星予測暦に含まれる予測位置のうち、各予測対象期間に対応する期間の予測位置を抽出する。そして、抽出した全ての予測位置からの距離の二乗和が最小となるようなケプラーの衛星軌道モデル式(以下、「近似モデル」ともいう。)を各予測対象期間それぞれについて求める。このとき求めた衛星軌道の近似モデル式のパラメータを「衛星軌道パラメータ」と称し、近似モデルを算出する計算を「近似計算」ともいう。また、近似計算により求めた衛星軌道のことを「予測軌道」と称する。長期予測エフェメリスは、全てのGPS衛星の全ての予測対象期間の衛星軌道パラメータの値が格納されたデータである。
【0034】
衛星予測暦に含まれるGPS衛星SVの予測位置は、将来になるほど、GPS衛星SVの実際の位置からずれる傾向がある。そのため、近似計算を行って長期予測エフェメリスを生成した場合に、近似計算により求めた予測軌道は、生成日時から将来のものであるほど、実際の衛星軌道からずれたものとなる可能性がある。
【0035】
本実施形態では、サーバシステム3は、各GPS衛星の各予測対象期間について、予測軌道の信頼性を示す指標値である「予測軌道信頼度」を決定し、決定した予測軌道信頼度を、信頼性パラメータとして長期予測エフェメリスに含めて携帯型電話機4に提供する。本実施形態では、予測軌道信頼度は「0」〜「12」の13段階で表され、「0」が予測軌道の信頼性が最も高く、「12」が予測軌道の信頼性が最も低いことを示している。尚、予測軌道信頼度の数値範囲は適宜設定変更可能であり、例えば「0」〜「15」の16段階で表すこととしてもよい。予測軌道信頼度は、エフェメリスに含まれている「URA index」に相当する値である。
【0036】
予測軌道信頼度の具体的な決定方法について、図面を参照して説明する。本実施形態では、開始日時の異なる複数の衛星予測暦と、それらの衛星予測暦に対応する衛星精密暦(開始日時が同一である衛星精密暦)との組合せ(以下、「衛星暦組合せ」と称す。)を用いて、長期予測エフェメリスに含める予測軌道信頼度を決定する。開始日時とは、衛星予測暦及び衛星精密暦にそれぞれ含まれる衛星位置のうち最も古い衛星位置の日時である。
【0037】
本実施形態では、開始日時が4時間ずつずれた衛星予測暦と、それらに対応する衛星精密暦とを用いて予測軌道信頼度を決定する。また、衛星予測暦及び衛星精密暦には、それぞれ1週間分の衛星位置及び時計誤差のデータが含まれており、便宜的に、この1週間を6時間毎にグループ分けした28個の期間を「第1期間〜第28期間」と表現する。但し、実際には、衛星予測暦及び衛星精密暦のデータは、衛星位置及び時計誤差が羅列的に配列されたデータ構造となっており、期間毎にグループ分けされているわけではない。
【0038】
より具体的には、例えば図2に示すように、現在日時から見て、8日と0時間前を開始日時とする衛星暦組合せ、7日と20時間前を開始日時とする衛星暦組合せ、7日と16時間前を開始日時とする衛星暦組合せ、・・・、1日と0時間前を開始日時とする衛星暦組合せといった具合に、開始日時が異なるN個の衛星予測暦組合せを抽出する。図2において、1つの帯が1つの衛星予測暦及び衛星精密暦を示している。そして、抽出した衛星予測暦に含まれる予測位置と、衛星精密暦に含まれる実績位置とを比較し、予測位置がどれだけ正確であるかを分析することで、予測軌道の信頼性を判定する。
【0039】
(1)観測誤差の算出
予測軌道信頼度を決定するに当たり、先ず、抽出した各衛星暦組合せそれぞれについて、当該衛星予測暦に含まれる各日時(衛星予測暦に含まれる15分おきの日時)における予測位置と、当該衛星精密暦に含まれる各日時(衛星精密暦に含まれる15分おきの日時)における実績位置とを用いて「観測誤差」を算出する。
【0040】
具体的には、先ず、ある日時の予測位置と、同じ日時の実績位置とを用いて観測位置を算出する。観測位置は、予測位置を地球上に投影した位置と、実績位置を地球上に投影した位置との中間位置として算出する。すなわち、予測位置と地球中心を結ぶ線分が地表面(標高0m)と交差する交点の座標と、実績位置と地球中心を結ぶ線分が地表面(標高0m)と交差する交点の座標を求める。そして、この2つの交点座標の中点の座標で表される位置を観測位置とする。いわゆる標高投影法を用いて観測位置を算出することになる。以上の観測位置の算出を、衛星予測暦に含まれている各日時について行う。
【0041】
観測位置を算出したら、衛星予測暦に含まれる各日時について、予測位置と観測位置間の距離(以下、「第1の観測距離」と称す。)、及び、実績位置と観測位置間の距離(以下、「第2の観測距離」と称す。)を算出し、第1の観測距離と第2の観測距離間の距離を「観測距離誤差」として求める。
【0042】
次いで、衛星予測暦に含まれる各日時について、時計予測誤差に光速を乗算することで、時計予測誤差の距離相当値(以下、「第1の時計誤差距離」と称す。)を算出する。同様に、衛星精密暦に記憶されている各日時について、時計実績誤差に光速を乗算することで、時計実績誤差の距離相当値(以下、「第2の時計誤差距離」と称す。)を算出する。そして、第1の時計誤差距離と第2の時計誤差距離との差を「時計距離誤差」として算出する。
【0043】
観測誤差は、観測距離誤差と時計距離誤差との合算値として算出することができる。このようにして、全ての衛星暦組合せについて、全てのGPS衛星の全ての期間における観測誤差を算出する。そして、図11に示すような観測誤差データ356を生成して、観測誤差データベース355に蓄積・記憶する。
【0044】
観測誤差データベース355は、衛星予測暦(衛星精密暦)の開始日時別に観測誤差データ356が蓄積記憶されたデータベースである。観測誤差データ356には、開始日時と対応付けて、第1期間〜第28期間までの各期間、及び、SV1〜SV32までの各GPS衛星における観測誤差が記憶されている。また、観測誤差データ356には、GPS衛星別に観測誤差の平均値を算出した衛星別平均観測誤差31(31−S1,31−S2,・・・,31−S32)と、期間別に観測誤差の平均値を算出した期間別平均観測誤差41(41−P1,41−P2,・・・,41−P28)とが記憶されている。
【0045】
(2)標準パターンの設定
次に、観測誤差を統計処理することで、予測軌道信頼度の標準的な変化パターンである「標準パターン」を設定する。例えば、N個の衛星暦組合せの中から、ランダムに1個の衛星暦組合せを選択して標準パターンを算出してもよいし、N個それぞれについて仮の標準パターンを作成して平均化することで最終的な標準パターンを算出することとしてもよい。以下では、ランダムに1個の衛星暦組合せを選択することとする。
【0046】
具体的には、観測誤差データベース355の中から、選択した衛星暦組合せの開始日時に対応する観測誤差データ356を抽出し、抽出した観測誤差データ356に記憶されている期間別平均観測誤差41を読み出す。
【0047】
図3は、期間別平均観測誤差をプロットしたグラフの一例である。図3において、横軸は日数(期間)、縦軸は期間別平均観測誤差をそれぞれ示している。期間は第1期間〜第28期間まであるが、本実施形態では各期間の長さを6時間としているため、第1期間〜第4期間までが1日目、第5期間〜第8期間までが2日目、・・・、第25期間〜第28期間までが7日目に相当する。
【0048】
衛星予測暦に記憶されている予測位置は、開始日時から将来のものであるほど、正確性が低くなる傾向がある。予測位置の正確性が低いと、予測位置及び観測位置から求めた観測距離(第1の観測距離)と、実績位置及び観測位置から求めた観測距離(第2の観測距離)との差(観測距離誤差)が大きくなるため、観測誤差は大きくなる。そのため、図3に示すように、将来の期間ほど、期間別平均観測誤差が大きくなる傾向がある。
【0049】
観測誤差が大きいほど、予測位置の正確性が低く、予測軌道の信頼性は低いと考えることができる。また、予測軌道信頼度は、値が大きいほど予測軌道の信頼性が低いことを示す指標値である。従って、期間別平均観測誤差が大きな期間ほど大きな標準値が設定されるように、標準パターンを設定する。
【0050】
具体的には、図8に示すような、期間別平均観測誤差と基礎値とが対応付けられたテーブルを予め用意しておく。そして、例えば日単位で、期間別平均観測誤差が含まれる数値範囲を特定し、特定した数値範囲に対応する基礎値を当該日にちに対応付けることで、標準パターンを設定する。図3では、標準パターンとして、1日目及び2日目(第1期間〜第8期間)に「3」、3日目及び4日目(第9期間〜第16期間)に「5」、5日目及び6日目(第17期間〜第24期間)に「7」、7日目(第25期間〜第28期間)に「8」がそれぞれ設定されている。
【0051】
(3)標準パターンの補正
次に、設定した標準パターンを補正することで、予測軌道信頼度を決定する。詳細には、図11の観測誤差データベース355に記憶されている観測誤差データ356を参照し、全てのGPS衛星及び全ての開始日時について算出した衛星別平均観測誤差31をGPS衛星間で相対評価することで、標準パターンの補正値を決定する。
【0052】
標準パターンの補正値は、例えば、相対的に衛星別平均観測誤差が小さいGPS衛星については「−1」、相対的に衛星別平均観測誤差が大きいGPS衛星については「+1」、それ以外のGPS衛星については「±0」とすることができる。勿論、より大きな値(例えば「+2」や「+3」)を補正値としてもよいし、より小さな値(例えば「−2」や「−3」)を補正値としてもよい。
【0053】
具体的に説明すると、先ず、全ての衛星別平均観測誤差の平均値と、衛星別平均観測誤差の標準偏差とを求める。そして、例えば、算出した平均値から標準偏差を減算した値を第1の閾値、算出した平均値に標準偏差を加算した値を第2の閾値とする。衛星別平均観測誤差が第1の閾値よりも小さくなる傾向があるGPS衛星は補正値を「−1」とし、第1の閾値以上第2の閾値以下となる傾向があるGPS衛星は補正値を「±0」とし、第2の閾値よりも大きくなる傾向があるGPS衛星は補正値を「+1」とする。尚、第1の閾値及び第2の閾値は、上述した手法以外にも適宜設定することが可能である。
【0054】
例えば、図4において、GPS衛星「SV1」は相対的に衛星別平均観測誤差が小さいと判定されたため補正値は「−1」とされ、また、GPS衛星「SV4」は相対的に衛星別平均観測誤差が大きいと判定されたため、補正値は「+1」とされている。また、GPS衛星「SV2」は衛星別平均観測誤差が大きくも小さくもないと判定されたため、補正値は「±0」とされている。
【0055】
尚、GPS衛星は、定期或いは不定期に軌道修正を行って、地球の周回軌道を補正している。軌道修正が行われると、衛星予測暦に含まれる予測位置のうち軌道修正が行われた後の予測位置が実績位置から大きくずれるため、軌道修正後の予測位置は、観測誤差が極端に大きな値となる。ここでは、この軌道修正等に起因して極端に大きくなった観測誤差のことを「異常観測誤差」と呼ぶ。
【0056】
本実施形態では、異常観測誤差が多数存在するGPS衛星(観測誤差が不良条件を満たすGPS衛星)については、補正値を「±0」とする。すなわち、異常観測誤差が多数存在するGPS衛星については、例外的に標準パターンの補正を行わないようにする。例えば、図4において、GPS衛星「SV3」は異常観測誤差が多数観測されたため、補正値が「±0」とされている。
【0057】
標準パターンの補正値が決まったら、当該補正値を用いて標準パターンを補正することで、予測軌道信頼度を決定する。具体的には、標準パターンを第1〜第28期間それぞれについて補正値を用いて補正し、補正後の値を、長期予測エフェメリスの第1予測対象期間〜第28予測対象期間の予測軌道信頼度とする。
【0058】
3.機能構成
図5は、サーバシステム3の機能構成を示すブロック図である。サーバシステム3は、CPU(Central Processing Unit)310と、操作部320と、通信部330と、ROM(Read Only Memory)340と、ハードディスク350と、RAM(Random Access Memory)360とを備え、各部がバス370で接続されたコンピュータシステムである。
【0059】
CPU310は、ROM340に記憶されているシステムプログラム等に従ってサーバシステム3の各部を統括的に制御するプロセッサである。本実施形態では、CPU310は、ROM340に記憶されている長期予測エフェメリス提供プログラム341に従って、携帯型電話機4に長期予測エフェメリスを提供する処理を行う。
【0060】
操作部320は、サーバシステム3の管理者による操作指示を受け付け、操作に応じた信号をCPU310に出力する入力装置である。この機能は、例えばキーボードやボタン、マウス等により実現される。
【0061】
通信部330は、システム内部で利用される各種データをインターネット等の通信ネットワークを介して外部システム2や携帯型電話機4とやり取りするための通信装置である。
【0062】
ROM340は、読み取り専用の不揮発性の記憶装置であり、CPU310がサーバシステム3を制御するためのシステムプログラムや、長期予測エフェメリスを携帯型電話機4に提供するためのプログラム、長期予測エフェメリスを生成するためのプログラム等の各種プログラムや各種データ等を記憶している。
【0063】
ハードディスク350は、磁気ヘッド等を用いてデータの読み書きを行う記憶装置であり、ROM340と同様、サーバシステム3が備える各種機能を実現するためのプログラムやデータ等を記憶している。
【0064】
RAM360は、読み書き可能な揮発性の記憶装置であり、CPU310により実行されるシステムプログラム、長期予測エフェメリス提供プログラム、各種処理プログラム、各種処理の処理中データ、処理結果などを一時的に記憶するワークエリアを形成している。
【0065】
4.データ構成
図6は、ROM340に格納されたデータの一例を示す図である。ROM340には、CPU310により読み出され、長期予測エフェメリス提供処理(図15参照)として実行される長期予測エフェメリス提供プログラム341と、観測誤差基礎値対応テーブル343とが記憶されている。また、長期予測エフェメリス提供プログラム341には、長期予測エフェメリス生成処理(図16及び図17参照)として実行される長期予測エフェメリス生成プログラム3411がサブルーチンとして含まれている。
【0066】
長期予測エフェメリス提供処理とは、CPU310が、長期予測エフェメリスデータ359を生成する処理を定期的に行い、携帯型電話機4から長期予測エフェメリスデータ359の要求信号を受信した場合に、生成しておいた長期予測エフェメリスデータ359を要求元の携帯型電話機4に送信する処理である。長期予測エフェメリス提供処理については、フローチャートを用いて詳細に後述する。
【0067】
長期予測エフェメリス生成処理とは、CPU310が、長期予測エフェメリスデータ359を生成する処理である。本実施形態では、CPU310は、長期予測エフェメリスデータ359を4時間に1回生成するものとする。長期予測エフェメリス生成処理についても、フローチャートを用いて詳細に後述する。
【0068】
図8は、観測誤差基礎値対応テーブル343のテーブル構成の一例を示す図である。観測誤差基礎値対応テーブル343には、期間別平均観測誤差が含まれる範囲である期間別平均観測誤差範囲3431と、期間別平均観測誤差が当該期間別平均観測誤差範囲3431に含まれる期間に設定する基礎値3433とが対応付けて記憶されている。例えば、期間別平均観測誤差が「20m〜40m」に含まれる期間は「5」が設定される。長期予測エフェメリス生成処理において、CPU310は、観測誤差基礎値対応テーブル343を用いて標準パターンを設定する。
【0069】
図7は、ハードディスク350に格納されたデータの一例を示す図である。ハードディスク350には、衛星予測暦データベース351と、衛星精密暦データベース353と、観測誤差データベース355と、標準パターンデータ357と、長期予測エフェメリスデータ359とが記憶されている。
【0070】
図9は、衛星予測暦データベース351のデータ構成の一例を示す図である。衛星予測暦データベース351は、複数の衛星予測暦352(352−1,352−2,352−3,・・・)が時系列に記憶されている。衛星予測暦352は、各GPS衛星SVの1週間後までの予測位置及び時計予測誤差が15分おきに記憶された離散的なデータであり、開始日時毎に纏められたデータである。便宜的に、1週間の期間を6時間毎にグループ分けすることで、第1期間〜第28期間の28個の期間が構成される。
【0071】
例えば、衛星予測暦352−1は、開始日時が「2008年8月1日0時00分」のデータである。そして、「2008年8月1日5時45分」におけるGPS衛星「SV2」の予測位置は「(Xp32,Yp32,Zp32)」であり、原子時計の予測誤差は「Δtp32」である。
【0072】
CPU310は、外部システム2から定期的に(例えば4時間に1回)衛星予測暦を受信する。そして、受信した衛星予測暦を衛星予測暦データベース351に蓄積・記憶するためにデータフォーマットの加工処理等を行う。具体的には、開始日時が異なり、長期予測エフェメリスの生成対象期間(例えば1週間)と同じ期間分のデータが格納された複数の衛星予測暦352を生成して、衛星予測暦データベース351に蓄積・記憶させる。
【0073】
図10は、衛星精密暦データベース353のデータ構成の一例を示す図である。衛星精密暦データベース353には、複数の衛星精密暦354(354−1,354−2,354−3,・・・)が時系列に記憶されている。衛星精密暦354は、各GPS衛星SVの1週間分の実績位置及び時計実績誤差が15分おきに記憶された離散的なデータであり、開始日時毎に纏められたデータである。1週間の期間を6時間毎にグループ分けすることで、第1期間〜第28期間の28個の期間が構成される。
【0074】
例えば、衛星精密暦354−1は、開始日時が「2008年8月1日0時00分」のデータである。そして、「2008年8月1日5時45分」におけるGPS衛星「SV2」の実績位置は「(Xm32,Ym32,Zm32)」であり、原子時計の実績誤差は「Δtm32」であった。
【0075】
CPU310は、外部システム2から定期的に(例えば4時間に1回)衛星精密暦を受信する。そして、受信した衛星精密暦を衛星精密暦データベース353に蓄積・記憶するためにデータフォーマットの加工処理等を行う。具体的には、衛星予測暦データベース351に格納されている複数の衛星予測暦352と開始日時が対応する衛星精密暦354を生成して、衛星精密暦データベース353に蓄積・記憶させる。
【0076】
図11は、観測誤差データベース355のデータ構成の一例を示す図である。観測誤差データベース355は、複数の観測誤差データ356(356−1,356−2,356−3,・・・)が開始日時別に記憶されたデータベースである。各観測誤差データ356には、各期間及び各GPS衛星における観測誤差が記憶されている。また、各衛星について全ての期間の観測誤差を平均した衛星別平均観測誤差31(31−S1,31−S2,・・・,31−S32)と、各期間について全てのGPS衛星の観測誤差を平均した期間別平均観測誤差41(41−P1,41−P2,・・・,41−P28)とが記憶されている。
【0077】
例えば、観測誤差データ356−1は、開始日時が「2008年8月1日0時00分」のデータである。また、GPS衛星「SV2」の「第2期間」における観測誤差は「E22」である。また、GPS衛星「SV2」の衛星別平均観測誤差は「ES2」であり、「第2期間」の期間別平均観測誤差は「EP2」である。
【0078】
CPU310は、長期予測エフェメリス提供処理において、上述した原理に従って各GPS衛星の各期間について観測誤差を算出する。また、算出した観測誤差を用いて、衛星別平均観測誤差31及び期間別平均観測誤差41を算出する。そして、これらを開始日時と対応付けた観測誤差データ356を生成して、観測誤差データベース355に蓄積・記憶させる。
【0079】
図12は、標準パターンデータ357のデータ構成の一例を示す図である。標準パターンデータ357には、期間3571と、標準パターン3573とが対応付けて記憶されている。例えば、標準パターンとして、第28期間には「8」が設定されている。
【0080】
CPU310は、長期予測エフェメリス生成処理において、上述した原理に従って標準パターンを設定する。そして、設定した標準パターン3573を期間3571と対応付けて、標準パターンデータ357に記憶させる。
【0081】
図13は、長期予測エフェメリスデータ359のデータ構成の一例を示す図である。長期予測エフェメリスデータ359には、長期予測エフェメリスデータの生成日時3591と、GPS衛星SV1〜SV32の予測エフェメリス3593(3593−1〜3593−32)とが対応付けて記憶されている。
【0082】
図14は、予測エフェメリス3593のデータ構成の一例を示す図である。予測エフェメリス3593(3593−1,3593−2,・・・,3593−32)には、各予測対象期間それぞれについて、軌道長半径や離心率、軌道傾斜角といったケプラーの衛星軌道パラメータの値と、衛星時計の基準時刻、衛星時計のオフセット、衛星時計のドリフト及び衛星時計周波数のドリフトでなるクロック補正パラメータの値と、信頼性パラメータである予測軌道信頼度とが対応付けて記憶されている。
【0083】
CPU310は、長期予測エフェメリス生成処理において、各GPS衛星SVについて、予測対象期間毎に衛星軌道パラメータ、クロック補正パラメータ及び信頼性パラメータの値を算出して予測エフェメリス3593を生成する。そして、全てのGPS衛星SVについて生成した予測エフェメリス3593を纏めて、生成日時3591と対応付けて長期予測エフェメリス359を生成して、ハードディスク350に記憶させる。
【0084】
5.処理の流れ
図15は、ROM340に記憶されている長期予測エフェメリス提供プログラム341がCPU310により読み出されて実行されることで、サーバシステム3において実行される長期予測エフェメリス提供処理の流れを示すフローチャートである。
【0085】
先ず、CPU310は、外部システム2から衛星予測暦及び衛星精密暦を受信したか否かを判定し(ステップA1)、受信しなかったと判定した場合は(ステップA1;No)、ステップA9へと処理を移行する。
【0086】
また、受信したと判定した場合は(ステップA1;Yes)、CPU310は、受信した衛星予測暦及び衛星精密暦を加工処理し、開始日時及び期間が同一である複数の衛星予測暦352及び衛星精密暦354を生成する(ステップA3)。そして、生成した衛星予測暦352及び衛星精密暦354を、ハードディスク350の衛星予測暦データベース351及び衛星精密暦データベース353にそれぞれ蓄積・記憶させる(ステップA5)。
【0087】
次いで、CPU310は、観測誤差算出処理を行う(ステップA7)。具体的には、上述した原理に従って観測誤差を算出し、観測誤差データ356を生成する。そして、生成した観測誤差データ356を観測誤差データベース355に蓄積・記憶させる。
【0088】
次いで、CPU310は、長期予測エフェメリスの生成日時であるか否かを判定する(ステップA9)。本実施形態では、4時間に1回長期予測エフェメリスを生成するものとする。そして、まだ生成日時ではないと判定した場合は(ステップA9;No)、CPU310は、ステップA13へと処理を移行する。
【0089】
また、長期予測エフェメリスの生成日時であると判定した場合は(ステップA9;Yes)、CPU310は、ROM340に記憶されている長期予測エフェメリス生成プログラム3411を読み出して実行することで、長期予測エフェメリス生成処理を行う(ステップA11)。
【0090】
図16及び図17は、長期予測エフェメリス生成処理の流れを示すフローチャートである。
先ず、CPU310は、標準パターン設定処理を行う(ステップB1)。具体的には、観測誤差データベース355に記憶されている観測誤差データ356を用いて、上述した原理に従って、各期間3571について標準パターン3573を設定する。そして、期間3571と標準パターン3573とを対応付けた標準パターンデータ357を、ハードディスク350に記憶させる。
【0091】
次いで、CPU310は、衛星別平均観測誤差をGPS衛星間で相対評価して、各GPS衛星の標準パターンの補正値を決定する(ステップB3)。そして、決定した補正値で標準パターンを補正することで、各期間の信頼性パラメータの値である予測軌道信頼度を決定する(ステップB5)。
【0092】
その後、CPU310は、長期予測エフェメリスの今回の生成日時(現在日時)に基づいて各予測対象期間を判定する(ステップB7)。すなわち、今回の生成日時から1週間後までの期間を生成対象期間とし、当該生成対象期間を6時間毎に区切った各期間を予測対象期間と判定する。
【0093】
次いで、CPU310は、各GPS衛星SVそれぞれについて、ループAの処理を実行する(ステップB9〜B25)。ループAの処理では、CPU310は、ステップB7で判定した各予測対象期間それぞれについて、ループBの処理を実行する(ステップB11〜B21)。
【0094】
ループBの処理では、CPU310は、ハードディスク350の衛星予測暦データベース351に記憶されている最新の衛星予測暦352から、当該GPS衛星SVの当該予測対象期間の各日時における予測位置を読み出す(ステップB13)。
【0095】
そして、CPU310は、読み出した予測位置を用いて、ケプラーの楕円軌道モデルに従って当該GPS衛星SVの当該予測対象期間における予測軌道を算出し、ケプラーの衛星軌道パラメータの値を求める(ステップB15)。尚、予測軌道の具体的な算出方法については公知であるため、詳細な説明を省略する。
【0096】
その後、CPU310は、最新の衛星予測暦352から、当該GPS衛星SVの当該予測対象期間の各日時における時計予測誤差を読み出す(ステップB17)。そして、CPU310は、読み出した時計予測誤差を用いて、当該GPS衛星の当該予測対象期間のクロック補正パラメータの値を求める(ステップB19)。
【0097】
時刻「t」における時計予測誤差「Δt」は、クロック補正パラメータである衛星時計の基準時刻「tc」、衛星時計のオフセット「a0」、衛星時計のドリフト「a1」及び衛星時計周波数のドリフト「a2」を用いて、次式(1)によって近似することができる。
Δt=a0+a1(t−tc)+a2(t−tc)2 ・・・(1)
【0098】
式(1)は、時計予測誤差の経時変化を近似するための時計誤差モデル式である。衛星予測暦351に含まれる各日時の時計予測誤差「Δt」をサンプリングデータとして、例えば最小二乗法を用いた近似計算を行うことで、クロック補正パラメータの値を算出することができる。その後、CPU310は、次の予測対象期間へと処理を移行する。
【0099】
全ての予測対象期間についてステップB13〜B19の処理を行った後、CPU310は、ループBの処理を終了する(ステップB21)。その後、CPU310は、全ての予測対象期間についてステップB15で求めた衛星軌道パラメータの値と、ステップB19で求めたクロック補正パラメータの値と、ステップB5で決定した信頼性パラメータの値(予測軌道信頼度)とを纏めて当該GPS衛星SVの予測エフェメリス3593を生成する(ステップB23)。そして、CPU310は、次のGPS衛星SVへと処理を移行する。
【0100】
全てのGPS衛星SVについてステップB11〜B23の処理を行った後、CPU310は、ループAの処理を終了する(ステップB25)。その後、CPU310は、全てのGPS衛星SVについてステップB23で生成した予測エフェメリス3593を纏めて、生成日時3591と対応付けて長期予測エフェメリスデータ359を生成し、ハードディスク350に記憶させる(ステップB27)。そして、CPU310は、長期予測エフェメリス生成処理を終了する。
【0101】
図15の長期予測エフェメリス提供処理に戻って、長期予測エフェメリス生成処理を行った後、CPU310は、携帯型電話機4から長期予測エフェメリスデータ359の要求信号を受信したか否かを判定する(ステップA13)。そして、受信しなかったと判定した場合は(ステップA13;No)、ステップA1に戻る。
【0102】
また、要求信号を受信したと判定した場合は(ステップA13;Yes)、CPU310は、ハードディスク350に記憶されている長期予測エフェメリスデータ359を要求元の携帯型電話機4に送信する(ステップA15)。そして、CPU310は、ステップA1に戻る。
【0103】
6.作用効果
本実施形態によれば、測位システム1において、サーバシステム3は、外部システム2から受信したGPS衛星SVの位置を時系列に予測した衛星予測暦に含まれる各予測位置を、測位用衛星の実際の位置を時系列に記憶した衛星精密暦に含まれる実績位置と比較することで、衛星予測暦に含まれる予測位置の正確性の変化を分析する。そして、連続する複数の予測対象期間それぞれにおける予測衛星軌道の衛星軌道パラメータの値を纏めてなる長期予測エフェメリスの予測対象期間それぞれの信頼性を、先の分析結果に基づいて判定する。そして、信頼性の判定結果を長期予測エフェメリスに含めて、携帯型電話機4に提供する。
【0104】
より具体的には、予測位置と観測位置間の距離と、実績位置と観測位置間の距離との差で表される観測距離誤差を算出する。また、時計予測誤差の距離相当値と時計実績誤差の距離相当値との差で表される時計距離誤差を算出する。次いで、観測距離誤差と時計距離誤差との合算値を観測誤差として算出する。そして、当該観測誤差を期間別に平均した期間別平均観測誤差を算出し、当該期間別平均観測誤差の時間変化に基づいて、標準パターンを設定する。
【0105】
その後、観測誤差をGPS衛星別に平均した衛星別平均観測誤差を算出し、この衛星別平均観測誤差をGPS衛星間で相対評価することで、標準パターンの補正値を決定する。そして、決定した補正値を用いて標準パターンを補正することで、各予測対象期間それぞれの予測軌道の信頼性の指標値である予測軌道信頼度を決定する。
【0106】
GPS衛星SVの実績位置は、GPS衛星SVのその日時における実際の位置である。そのため、予測位置及び実績位置と用いて算出した観測誤差が大きい場合は、予測位置の正確性が低く、この予測位置に基づいて予測した予測衛星軌道の信頼性は低いと言える。そのため、期間別平均観測誤差が大きい期間ほど信頼性が低くなるように標準パターンを設定した上で、GPS衛星間で衛星別平均観測誤差が相対的に小さい衛星は信頼性が高くなるように標準パターンを補正し、逆に、衛星別平均観測誤差が相対的に大きい衛星は信頼性が低くなるように標準パターンを補正する。これにより、長期予測エフェメリスの信頼性を適確に判定し、適切な信頼性パラメータの値を携帯型電話機4に提供することが可能となる。
【0107】
7.変形例
7−1.測位システム
上述した実施形態では、サーバシステム3と携帯型電話機4を備えた測位システム1を例に挙げて説明したが、本発明を適用可能な測位システムはこれに限られるわけではない。例えば、携帯型電話機4の代わりに、測位装置を備えたノート型パソコンやPDA(Personal Digital Assistant)、カーナビゲーション装置等の電子機器に適用することも可能である。
【0108】
また、上述した実施形態では、情報提供装置の一種としてサーバシステム3を例に挙げて説明したが、情報提供装置はサーバシステム3に限られるわけではない。例えば、汎用のパソコン等であってもよい。
【0109】
7−2.衛星測位システム
また、上述した実施形態では、衛星測位システムとしてGPSを例に挙げて説明したが、WAAS(Wide Area Augmentation System)、QZSS(Quasi Zenith Satellite System)、GLONASS(GLObal NAvigation Satellite System)、GALILEO等の他の衛星測位システムであってもよい。
【0110】
7−3.観測誤差の時系列変化に基づく標準パターンの補正
衛星別観測誤差を相対評価することに代えて又は加えて、衛星別観測誤差の開始日時の時系列変化に基づいて標準パターンを補正することも可能である。
【0111】
図18は、この場合にCPU310が行う処理である第2の長期予測エフェメリス生成処理のうち、図16の長期予測エフェメリス生成処理に対応する部分を示すフローチャートである。尚、長期予測エフェメリス生成処理と同一のステップについては同一の符号を付して説明を省略し、長期予測エフェメリス生成処理とは異なる部分について説明する。
【0112】
第2長期予測エフェメリス生成処理では、CPU310は、ステップB3において標準パターンの補正値を決定した後、各GPS衛星について、開始日時の時系列変化に基づく衛星別平均観測誤差の増減を判定し、標準パターンの補正値を修正する(ステップC4)。
【0113】
具体的には、開始日時が生成日時(現在日時)に近づくにつれて衛星別平均観測誤差が増加する傾向がある場合は、補正値に「1」を加算し、衛星別平均観測誤差は減少する傾向がある場合は、補正値から「1」を減算する。また、どちらとも判断がつかない場合は、補正値をそのままとする。勿論、「1」よりも大きな値(例えば「2」や「3」)を補正値に加算/減算することとしてもよい。衛星別平均観測誤差の増減は、微分値を算出するなどして判断することが可能である。
【0114】
その後、CPU310は、修正した補正値で標準パターンを補正することで各期間の予測軌道信頼度を決定する(ステップC5)。そして、CPU310は、ステップB7以降へと処理を移行する。
【0115】
7−4.衛星位置誤差
上述した実施形態では、観測誤差を用いて予測軌道信頼度を決定するものとして説明したが、観測誤差の代わりに、衛星位置の誤差(以下、「衛星位置誤差」と称す。)を用いて予測軌道信頼度を決定することとしてもよい。衛星位置誤差は、衛星予測暦に含まれる予測位置と、衛星精密暦に含まれる実績位置間の距離として求めることができる。観測誤差を衛星位置誤差に置き換えることで、上述した原理に従って、予測軌道信頼度を決定することができる。
【0116】
7−5.標準パターンの設定
上述した実施形態では、図3に示したように、標準パターンを日単位で設定する場合を一例として説明したが、標準パターンを期間単位で設定してもよいことは勿論である。すなわち、図8の観測誤差基礎値対応テーブル343に従って、各期間それぞれについて、期間別平均観測誤差が何れの期間別平均観測誤差範囲3431に含まれるかを判定する。そして、判定した期間別平均観測誤差範囲3431に対応する基礎値3433を読み出して各期間における基礎値を設定することで、標準パターンを設定する。
【0117】
7−6.平均観測誤差
上述した実施形態では、期間別平均観測誤差に基づいて標準パターンを設定するものとして説明したが、期間別に算出した観測誤差の最大値(期間別最大観測誤差)に基づいて標準パターンを設定することにしてもよい。また、衛星別平均観測誤差を相対評価して標準パターンの補正値を決定するのではなく、衛星別に算出した観測誤差の最大値(衛星別最大観測誤差)を相対評価して標準パターンの補正値を決定することにしてもよい。
【0118】
7−7.長期予測エフェメリスの生成
上述した実施形態では、サーバシステム3が長期予測エフェメリスデータを生成して携帯型電話機4に提供するものとして説明したが、携帯型電話機4が自ら長期予測エフェメリスデータを生成することとしてもよい。すなわち、携帯型電話機4が外部システム2から定期的に衛星予測暦及び衛星精密暦を取得し、取得した衛星予測暦及び衛星精密暦を用いて長期予測エフェメリス生成処理を実行することで、長期予測エフェメリスデータを生成する。携帯型電話機4の代わりに、測位装置を備えたノート型パソコンやPDA、カーナビゲーション装置等の電子機器に適用する場合も同様である。
【0119】
また、上述した実施形態では、サーバシステム3が予め所定の時間間隔(例えば4時間に1回)で長期予測エフェメリスデータを生成しておき、携帯型電話機4から長期予測エフェメリスデータの要求を受けた場合に、生成しておいた長期予測エフェメリスデータを送信するものとして説明した。このような構成ではなく、サーバシステム3が、携帯型電話機4から長期予測エフェメリスデータの要求を受けたときに長期予測エフェメリスデータを生成して、携帯型電話機4に送信する構成としてもよい。
【0120】
7−8.生成対象期間
上述した実施形態では、長期予測エフェメリスの生成日時を基準として1週間後までの期間を生成対象期間として長期予測エフェメリスを生成するものとして説明したが、生成対象期間は1週間よりも長い期間(例えば2週間)としてもよいし、1週間よりも短い期間(例えば3日)としてもよい。GPS衛星SVから送信される航法データとしてのエフェメリスは有効期間が一般に4時間程度であるが、長期予測エフェメリスは少なくともGPS衛星SVから送信される航法データとしてのエフェメリスよりも有効期間が長ければよい。1日以上であれば好適である。
【0121】
7−9.予測対象期間
また、上述した実施形態では、予測対象期間の長さを6時間とするものとして説明したが、これに限定されるわけではなく、4時間や8時間等としてもよく、適宜設定可能であることは勿論である。
【図面の簡単な説明】
【0122】
【図1】測位システムの概略構成を示す図。
【図2】衛星予測暦及び衛星精密暦の説明図。
【図3】期間別平均観測誤差をプロットしたグラフ。
【図4】標準パターンの補正値決定の説明図。
【図5】サーバシステムの機能構成を示すブロック図。
【図6】サーバシステムのROMに格納されたデータの一例を示す図。
【図7】サーバシステムのハードディスクに格納されたデータの一例を示す図。
【図8】観測誤差基礎値対応テーブルのテーブル構成の一例を示す図。
【図9】衛星予測暦データベースのデータ構成の一例を示す図。
【図10】衛星精密暦データベースのデータ構成の一例を示す図。
【図11】観測誤差データベースのデータ構成の一例を示す図。
【図12】標準パターンデータのデータ構成の一例を示す図。
【図13】長期予測エフェメリスデータのデータ構成の一例を示す図。
【図14】予測エフェメリスのデータ構成の一例を示す図。
【図15】長期予測エフェメリス提供処理の流れを示すフローチャート。
【図16】長期予測エフェメリス生成処理の流れを示すフローチャート。
【図17】長期予測エフェメリス生成処理の流れを示すフローチャート。
【図18】第2の長期予測エフェメリス生成処理の流れを示すフローチャート。
【符号の説明】
【0123】
1 測位システム 、 2 外部システム、 3 サーバシステム、
4 携帯型電話機、 310 CPU、 320 操作部、 330 通信部、
340 ROM、 350 ハードディスク、 360 RAM、 370 バス、
SV GPS衛星
【特許請求の範囲】
【請求項1】
測位用衛星の位置を時系列に予測した予測位置データに含まれる各予測位置を、当該予測位置に対応する前記測位用衛星の実績位置と比較することで、前記予測位置データに含まれる時系列の予測位置の正確性の変化を分析することと、
連続する複数の予測対象期間それぞれにおける予測衛星軌道を纏めてなる長期予測軌道データの前記予測対象期間それぞれの信頼性を、前記分析の結果に基づいて判定することと、
を含む長期予測軌道データの信頼性判定方法。
【請求項2】
前記分析することは、
予測した日時が異なる複数の前記予測位置データそれぞれについて、当該予測位置データに含まれる各予測位置を対応する実績位置と比較することと、
前記比較の結果を統計処理して、前記測位用衛星の位置を時系列に予測した場合の時間経過に対する予測位置の正確性の変化パターンを算出することと、
を含み、
前記信頼性を判定することは、前記変化パターンに基づいて、前記長期予測軌道データの前記予測対象期間それぞれの信頼性を判定することである、
請求項1に記載の長期予測軌道データの信頼性判定方法。
【請求項3】
前記予測位置データには、複数の前記測位用衛星それぞれの予測位置が含まれ、
前記比較することは、前記測位用衛星それぞれの各予測位置を実績位置と比較することであり、
前記変化パターンを算出することは、前記測位用衛星それぞれの前記比較の結果を統計処理して、前記測位用衛星間で標準的な前記変化パターンを算出することである、
請求項2に記載の長期予測軌道データの信頼性判定方法。
【請求項4】
前記長期予測軌道データは、前記予測対象期間それぞれについての各測位用衛星毎のデータでなり、
前記測位用衛星それぞれの前記比較の結果を前記測位用衛星間で相対評価することを更に含み、
前記信頼性を判定することは、前記相対評価の結果に基づいて前記変化パターンを前記測位用衛星毎に補正して、前記長期予測軌道データの前記予測対象期間それぞれの信頼性を各測位用衛星毎に判定することを含む、
請求項3に記載の長期予測軌道データの信頼性判定方法。
【請求項5】
前記分析することは、前記測位用衛星それぞれについて、前記予測位置データの予測日時毎に、当該予測位置データに関する前記比較の結果を纏めることを含み、
前記信頼性を判定することは、前記予測位置データの予測日時毎に纏めた前記比較の結果の前記予測日時の時系列変化に基づいて前記変化パターンを前記測位用衛星毎に補正して、前記長期予測軌道データの前記予測対象期間それぞれの信頼性を各測位用衛星毎に判定することを含む、
請求項3又は4に記載の長期予測軌道データの信頼性判定方法。
【請求項6】
前記信頼性を判定することは、前記複数の測位用衛星のうち、当該測位用衛星の予測位置を実績位置と比較した前記比較の結果が所定の不良条件を満たさない測位用衛星についてのみ前記変化パターンの前記補正を行った上で信頼性の判定を行うことを含む請求項5に記載の長期予測軌道データの信頼性判定方法。
【請求項7】
長期予測軌道データを生成することと、
請求項1〜6の何れか一項に記載された長期予測軌道データの信頼性判定方法を用いて、前記生成した長期予測軌道データの信頼性を判定することと、
前記生成した長期予測軌道データと前記判定した判定結果とを測位装置に提供することと、
を含む長期予測軌道データの提供方法。
【請求項8】
連続する複数の予測対象期間それぞれにおける予測衛星軌道を纏めてなる長期予測軌道データを生成する生成部と、
測位用衛星の位置を時系列に予測した予測位置データに含まれる各予測位置を、当該予測位置に対応する前記測位用衛星の実績位置と比較することで、前記予測位置データに含まれる時系列の予測位置の正確性の変化を分析する分析部と、
前記生成部により生成された長期予測軌道データの前記予測対象期間それぞれの信頼性を、前記分析部の分析結果に基づいて判定する判定部と、
前記生成部により生成された長期予測軌道データと前記判定部による判定結果とを測位装置に提供する提供部と、
を備えた情報提供装置。
【請求項1】
測位用衛星の位置を時系列に予測した予測位置データに含まれる各予測位置を、当該予測位置に対応する前記測位用衛星の実績位置と比較することで、前記予測位置データに含まれる時系列の予測位置の正確性の変化を分析することと、
連続する複数の予測対象期間それぞれにおける予測衛星軌道を纏めてなる長期予測軌道データの前記予測対象期間それぞれの信頼性を、前記分析の結果に基づいて判定することと、
を含む長期予測軌道データの信頼性判定方法。
【請求項2】
前記分析することは、
予測した日時が異なる複数の前記予測位置データそれぞれについて、当該予測位置データに含まれる各予測位置を対応する実績位置と比較することと、
前記比較の結果を統計処理して、前記測位用衛星の位置を時系列に予測した場合の時間経過に対する予測位置の正確性の変化パターンを算出することと、
を含み、
前記信頼性を判定することは、前記変化パターンに基づいて、前記長期予測軌道データの前記予測対象期間それぞれの信頼性を判定することである、
請求項1に記載の長期予測軌道データの信頼性判定方法。
【請求項3】
前記予測位置データには、複数の前記測位用衛星それぞれの予測位置が含まれ、
前記比較することは、前記測位用衛星それぞれの各予測位置を実績位置と比較することであり、
前記変化パターンを算出することは、前記測位用衛星それぞれの前記比較の結果を統計処理して、前記測位用衛星間で標準的な前記変化パターンを算出することである、
請求項2に記載の長期予測軌道データの信頼性判定方法。
【請求項4】
前記長期予測軌道データは、前記予測対象期間それぞれについての各測位用衛星毎のデータでなり、
前記測位用衛星それぞれの前記比較の結果を前記測位用衛星間で相対評価することを更に含み、
前記信頼性を判定することは、前記相対評価の結果に基づいて前記変化パターンを前記測位用衛星毎に補正して、前記長期予測軌道データの前記予測対象期間それぞれの信頼性を各測位用衛星毎に判定することを含む、
請求項3に記載の長期予測軌道データの信頼性判定方法。
【請求項5】
前記分析することは、前記測位用衛星それぞれについて、前記予測位置データの予測日時毎に、当該予測位置データに関する前記比較の結果を纏めることを含み、
前記信頼性を判定することは、前記予測位置データの予測日時毎に纏めた前記比較の結果の前記予測日時の時系列変化に基づいて前記変化パターンを前記測位用衛星毎に補正して、前記長期予測軌道データの前記予測対象期間それぞれの信頼性を各測位用衛星毎に判定することを含む、
請求項3又は4に記載の長期予測軌道データの信頼性判定方法。
【請求項6】
前記信頼性を判定することは、前記複数の測位用衛星のうち、当該測位用衛星の予測位置を実績位置と比較した前記比較の結果が所定の不良条件を満たさない測位用衛星についてのみ前記変化パターンの前記補正を行った上で信頼性の判定を行うことを含む請求項5に記載の長期予測軌道データの信頼性判定方法。
【請求項7】
長期予測軌道データを生成することと、
請求項1〜6の何れか一項に記載された長期予測軌道データの信頼性判定方法を用いて、前記生成した長期予測軌道データの信頼性を判定することと、
前記生成した長期予測軌道データと前記判定した判定結果とを測位装置に提供することと、
を含む長期予測軌道データの提供方法。
【請求項8】
連続する複数の予測対象期間それぞれにおける予測衛星軌道を纏めてなる長期予測軌道データを生成する生成部と、
測位用衛星の位置を時系列に予測した予測位置データに含まれる各予測位置を、当該予測位置に対応する前記測位用衛星の実績位置と比較することで、前記予測位置データに含まれる時系列の予測位置の正確性の変化を分析する分析部と、
前記生成部により生成された長期予測軌道データの前記予測対象期間それぞれの信頼性を、前記分析部の分析結果に基づいて判定する判定部と、
前記生成部により生成された長期予測軌道データと前記判定部による判定結果とを測位装置に提供する提供部と、
を備えた情報提供装置。
【図1】
【図2】
【図3】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図4】
【図2】
【図3】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図4】
【公開番号】特開2010−101804(P2010−101804A)
【公開日】平成22年5月6日(2010.5.6)
【国際特許分類】
【出願番号】特願2008−274560(P2008−274560)
【出願日】平成20年10月24日(2008.10.24)
【出願人】(000002369)セイコーエプソン株式会社 (51,324)
【Fターム(参考)】
【公開日】平成22年5月6日(2010.5.6)
【国際特許分類】
【出願日】平成20年10月24日(2008.10.24)
【出願人】(000002369)セイコーエプソン株式会社 (51,324)
【Fターム(参考)】
[ Back to top ]