特許第6804806号(P6804806)IP Force 特許公報掲載プロジェクト 2022.1.31 β版

知財求人 - 知財ポータルサイト「IP Force」

▶ 日本電信電話株式会社の特許一覧

特許6804806測位精度情報算出装置及び測位精度情報算出方法
<>
  • 特許6804806-測位精度情報算出装置及び測位精度情報算出方法 図000002
  • 特許6804806-測位精度情報算出装置及び測位精度情報算出方法 図000003
  • 特許6804806-測位精度情報算出装置及び測位精度情報算出方法 図000004
  • 特許6804806-測位精度情報算出装置及び測位精度情報算出方法 図000005
  • 特許6804806-測位精度情報算出装置及び測位精度情報算出方法 図000006
  • 特許6804806-測位精度情報算出装置及び測位精度情報算出方法 図000007
  • 特許6804806-測位精度情報算出装置及び測位精度情報算出方法 図000008
  • 特許6804806-測位精度情報算出装置及び測位精度情報算出方法 図000009
  • 特許6804806-測位精度情報算出装置及び測位精度情報算出方法 図000010
  • 特許6804806-測位精度情報算出装置及び測位精度情報算出方法 図000011
  • 特許6804806-測位精度情報算出装置及び測位精度情報算出方法 図000012
  • 特許6804806-測位精度情報算出装置及び測位精度情報算出方法 図000013
  • 特許6804806-測位精度情報算出装置及び測位精度情報算出方法 図000014
< >
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】6804806
(24)【登録日】2020年12月7日
(45)【発行日】2020年12月23日
(54)【発明の名称】測位精度情報算出装置及び測位精度情報算出方法
(51)【国際特許分類】
   G01S 19/24 20100101AFI20201214BHJP
【FI】
   G01S19/24
【請求項の数】8
【全頁数】12
(21)【出願番号】特願2017-143842(P2017-143842)
(22)【出願日】2017年7月25日
(65)【公開番号】特開2019-27799(P2019-27799A)
(43)【公開日】2019年2月21日
【審査請求日】2019年6月21日
(73)【特許権者】
【識別番号】000004226
【氏名又は名称】日本電信電話株式会社
(74)【代理人】
【識別番号】110001863
【氏名又は名称】特許業務法人アテンダ国際特許事務所
(72)【発明者】
【氏名】塩田 圭
(72)【発明者】
【氏名】川上 健太
(72)【発明者】
【氏名】森田 章弘
(72)【発明者】
【氏名】大津 茂雄
(72)【発明者】
【氏名】則武 克誌
【審査官】 梶田 真也
(56)【参考文献】
【文献】 特開2011−064523(JP,A)
【文献】 特開2009−042106(JP,A)
【文献】 特開2004−184121(JP,A)
【文献】 特開2016−149142(JP,A)
【文献】 特開昭64−079679(JP,A)
【文献】 特表2010−534849(JP,A)
【文献】 米国特許出願公開第2013/0285849(US,A1)
(58)【調査した分野】(Int.Cl.,DB名)
G01S 5/00 − 5/14
G01S 19/00 − 19/55
(57)【特許請求の範囲】
【請求項1】
複数の領域のそれぞれの位置を表す複数のメッシュ点のうち、一のメッシュ点を代表点とし、前記代表点を含み、前記代表点の近傍の前記複数のメッシュ点を複数の近傍点として定める代表点決定部と、
前記代表点の位置情報を取得すると共に、衛星の軌道情報を取得し、前記代表点における対象時刻の衛星位置の仰角を算出する衛星座標算出部と、
前記複数の近傍点のそれぞれについて、天空画像を取得し、前記天空画像から開空間を識別し、前記開空間に基づいて天頂を中心とする円形の近似開空間を算出する近似開空間算出部と、
前記代表点の前記仰角、及び前記複数の近傍点それぞれの前記近似開空間に基づいて、前記複数の近傍点それぞれの可視衛星数を算出する可視衛星数算出部と、を備える
ことを特徴とする測位精度情報算出装置。
【請求項2】
前記代表点は、前記複数の近傍点に係る前記複数の領域において、重心を含む前記領域の前記メッシュ点である
ことを特徴とする請求項1記載の測位精度情報算出装置。
【請求項3】
前記代表点決定部は、対象地域をメッシュ状に分割し、前記近傍点のグループが異なる場合には、異なる代表点を定める
ことを特徴とする請求項1又は2記載の測位精度情報算出装置。
【請求項4】
前記近似開空間算出部は、三次元地図情報から生成された画像、車載カメラの画像から生成された画像、及び天空を撮影した画像のうちいずれか利用可能なものを取得する
ことを特徴とする請求項1乃至3何れか1項に記載の測位精度情報算出装置。
【請求項5】
前記近似開空間算出部は、複数の天空画像が取得可能である場合に、最新の前記天空画像を取得する
ことを特徴とする請求項1乃至4何れか1項に記載の測位精度情報算出装置。
【請求項6】
前記複数の領域それぞれにおける前記可視衛星数を表した画像を表示する可視衛星数表示部を更に備える
ことを特徴とする請求項1乃至5何れか1項に記載の測位精度情報算出装置。
【請求項7】
前記近似開空間は、閉空間を一部に含む
ことを特徴とする請求項1乃至6何れか1項に記載の測位精度情報算出装置。
【請求項8】
複数の領域のそれぞれの位置を表す複数のメッシュ点のうち、一のメッシュ点を代表点とし、前記代表点を含み、前記代表点の近傍の前記複数のメッシュ点を複数の近傍点として定め、
前記代表点の位置情報を取得すると共に、衛星の軌道情報を取得し、前記代表点における対象時刻の衛星位置の仰角を算出し、
前記複数の近傍点のそれぞれについて、天空画像を取得し、前記天空画像から開空間を識別し、前記開空間に基づいて天頂を中心とする円形の近似開空間を算出し、
前記代表点の前記仰角、及び前記複数の近傍点それぞれの前記近似開空間に基づいて、前記複数の近傍点それぞれの可視衛星数を算出する
ことを特徴とする測位精度情報算出方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、測位精度情報算出装置及び測位精度情報算出方法に関する。
【背景技術】
【0002】
近年GPS(Global Positioning System)に加えてGRONASS(Global Navigation Satellite System)やQZSS(Quasi-Zenith Satellite System)等様々な航法衛星が打ち上げられており、位置情報の利用が広まっている。航法衛星を用いた測位では4基以上の衛星を確保することで測位可能となるが、建物等の構造物の遮蔽による可視衛星数の減少、並びにビル等構造物によるマルチパス及び回折等の影響により、数10mオーダーの測位精度に劣化してしまう場合もある。
【0003】
非特許文献1は、衛星の正確な軌道情報及び3次元デジタル地図を利用して、可視衛星の数、精度低下率(DOP:Dilution Of Precision)及び誤差分布をシミュレートすることについて開示している。
【0004】
非特許文献2は、魚眼レンズを用いて撮影したアンテナ設置位置の天空画像データから衛星信号の受信の障害となる構造物を識別し、航法衛星の軌道情報と合わせて衛星信号の受信特性を推定して、GPSアンテナの設置作業を効率化する方法を検討している。
【先行技術文献】
【非特許文献】
【0005】
【非特許文献1】Tomohiro HAKAMATA, "Development of a Simulation System to Delineate Availability of GNSS with3-D Digital Map", Proceedings of the 23rd Asian Conference on Remote Sensing, 2002, 2002
【非特許文献2】西 啓介、吉田 誠史、広瀬 貴史,“GPSアンテナの設置作業の効率化を実現する衛星信号受信特性推定方式の提案”,信学技報 115(155), 23-28, 2015-07-24
【発明の概要】
【発明が解決しようとする課題】
【0006】
航法衛星は、将来的に自動運転や社会インフラ管理等の様々な分野でより利用の機会が広がると考えられている。このように航法衛星がより頻繁に身近に利用されるようになると、航法衛星を利用したサービスの内容によっては、衛星測位に有利/不利なエリアや、各時間帯で予測される測位精度等の情報を必要とすることが考えられる。また、このような測位精度に関する情報は、衛星電波遮蔽環境改善への基礎データとしても用いることができる。
【0007】
しかしながら、地理的に細かく分けられた領域のそれぞれに対して、時間帯毎の測位精度情報を提供するのは計算の処理負荷が極めて大きかった。
【0008】
本発明は上記事情に鑑みてなされたものであり、その目的とするところは、より処理負荷を軽減して測位精度情報を出力することのできる測位精度情報算出装置及び測位精度情報算出方法を提供することにある。
【課題を解決するための手段】
【0009】
上記目的を達成するために、本願発明に係る測位精度情報算出装置は、複数の領域のそれぞれの位置を表す複数のメッシュ点のうち、一のメッシュ点を代表点とし、前記代表点を含み、前記代表点の近傍の前記複数のメッシュ点を複数の近傍点として定める代表点決定部と、前記代表点の位置情報を取得すると共に、衛星の軌道情報を取得し、前記代表点における対象時刻の衛星位置の仰角を算出する衛星座標算出部と、前記複数の近傍点のそれぞれについて、天空画像を取得し、前記天空画像から開空間を識別し、前記開空間に基づいて天頂を中心とする円形の近似開空間を算出する近似開空間算出部と、前記代表点の前記仰角、及び前記複数の近傍点それぞれの前記近似開空間に基づいて、前記複数の近傍点それぞれの可視衛星数を算出する可視衛星数算出部と、を備えることを特徴とする。
【0010】
上記目的を達成するために、本願発明に係る測位精度情報算出方法は、複数の領域のそれぞれの位置を表す複数のメッシュ点のうち、一のメッシュ点を代表点とし、前記代表点を含み、前記代表点の近傍の前記複数のメッシュ点を複数の近傍点として定め、前記代表点の位置情報を取得すると共に、衛星の軌道情報を取得し、前記代表点における対象時刻の衛星位置の仰角を算出し、前記複数の近傍点のそれぞれについて、天空画像を取得し、前記天空画像から開空間を識別し、前記開空間に基づいて天頂を中心とする円形の近似開空間を算出し、前記代表点の前記仰角、及び前記複数の近傍点それぞれの前記近似開空間に基づいて、前記複数の近傍点それぞれの可視衛星数を算出することを特徴とする。
【発明の効果】
【0011】
本発明によれば、より処理負荷を軽減して航法衛星の測位精度情報を出力することができる。
【図面の簡単な説明】
【0012】
図1】本実施形態に係る測位精度情報算出システムの構成を示す概略図
図2】代表点及び近傍点について説明するための図
図3】衛星位置の仰角について説明する図
図4】天空画像の取得の例について説明するための図
図5】測位精度情報算出処理の例について示すフローチャート
図6】衛星座標算出処理を詳細に示すフローチャート
図7】近似開空間算出処理を詳細に示すフローチャート
図8】開空間及び閉空間を識別しやすくした画像の例を示す図
図9】近似開空間の例を説明する図
図10】近似開空間の例を説明する図
図11】可視衛星数算出処理を詳細に示すフローチャート
図12】可視衛星かどうかの判定方法の例について説明する図
図13】可視衛星数を示した対象地域の画像の例について示す図
【発明を実施するための形態】
【0013】
本発明の一実施の形態に係る測位精度情報算出システム1について図面を参照して説明する。図1は、本実施形態に係る測位精度情報算出システム1の構成を示す概略図である。測位精度情報算出システム1は、測位精度情報算出装置10、衛星情報記憶装置21及び天空画像取得装置22により構成される。これらの測位精度情報算出装置10、衛星情報記憶装置21及び天空画像取得装置22は、それぞれ、主に半導体装置で構成され、CPU(Central Processing Unit)、RAM(Random Access Memory)等の揮発性記憶装置、ハードディスクやフラッシュメモリ等の不揮発性記憶装置、及び外部との通信のための接続を行う通信インターフェースを有する、いわゆる情報処理機器として構成されていてもよい。
【0014】
本実施形態においては、これらの装置はそれぞれ通信ネットワーク31に接続されていることとする。しかしながら、衛星情報記憶装置21及び天空画像取得装置22のいずれか又は両方は、測位精度情報算出装置10内の機能により実現されるものであってもよい。通信ネットワーク31は、LAN(Local Area Network)やインターネット等のTCP(Transmission Control Protocol)/IP(Internet Protocol)プロトコルにより通信接続を行う通信ネットワーク31とすることができる。通信ネットワーク31は、有線及び無線のいずれか又は両方を用いた通信ネットワーク31とすることができる。
【0015】
衛星情報記憶装置21は、GPS、GRONASS及びQZSS等の様々な航法衛星のうち、複数の衛星61の軌道情報を記憶することができる。衛星情報記憶装置21は、通信ネットワーク31を介してアクセスできるサーバ装置として構成されていてもよい。航法衛星の衛星信号により軌道情報を取得する場合には、航法衛星を衛星情報記憶装置21とすることとしてもよい。
【0016】
天空画像取得装置22は、ある位置情報が示す位置における天空画像が要求された際に、その位置での天空画像を送信する。ここで天空画像は、例えば天頂を中心とする平面視による画像、又は魚眼レンズ等を用いて天空を撮影した画像等を含む。ここで天空画像取得装置22は、例えば、三次元地図情報や車載カメラ等が撮影した風景画像に基づいて天空画像を生成する機能を有していてもよい。天空画像取得装置22は、通信ネットワーク31を介してアクセスできるサーバ装置として構成されていてもよい。
【0017】
測位精度情報算出装置10は、代表点決定部11と、衛星座標算出部12と、近似開空間算出部13と、可視衛星数算出部14と、可視衛星数表示部15とを有している。代表点決定部11は、複数の領域54のそれぞれの位置を表す複数のメッシュ点52のうち、一のメッシュ点52を代表点51とし、代表点51を含み、代表点51の近傍の複数のメッシュ点52を複数の近傍点55として定める。ここで、代表点決定部11は、更に、対象地域をメッシュ状に分割し、近傍点55のグループが異なる場合には、異なる代表点51を定めることとしてもよい。
【0018】
図2は、代表点51及び近傍点55について説明するための図である。図2には2つの近傍点領域50が示されており、それぞれメッシュ53により複数の領域54に分割されている。各領域54はその位置を表すメッシュ点52を有している。代表点決定部11は、近傍点領域50に含まれる複数のメッシュ点52のうち一つのメッシュ点52を代表点51として決定し、代表点51を含む、近傍点領域50に含まれるメッシュ点52を近傍点55として決定する。ここで、メッシュ点52は、個々の領域54を2次元平面としたときの重心とすることができる。また、代表点51は、複数の近傍点55に係る複数の領域54(近傍点領域50)において、重心を含む領域54のメッシュ点52とすることができる。なお、図2に示される代表点51及び近傍点55は一例であり、この他の方法で、代表点51、及び代表点51を含む複数の近傍点55として定めることとしてもよい。
【0019】
衛星座標算出部12は、代表点51の位置情報を取得すると共に、例えば衛星情報記憶装置21から衛星61の軌道情報を取得し、代表点51における対象時刻の衛星位置の仰角θを算出する。図3は、衛星位置の仰角θについて説明する図である。この図に示されるように、仰角θは、代表点51において、衛星61を見上げたときの水平面からの角度である。測位衛星は高度約2万kmに軌道を取っているため、代表点51と他の近傍点55の距離は無視できる程度である。このため、各衛星61の仰角θは、代表点51及び代表点51近傍の近傍点55においてほぼ同じであり、仰角θは、近傍点55においても、代表点51の仰角θを用いて計算することとし、近傍点55での仰角θの計算は行わないこととしている。これにより計算の処理負荷をより軽減することができる。
【0020】
具体的には、例えば、与えられた衛星軌道要素等の衛星軌道を同定する情報からECEF(Earth Centered Earth Fixed)座標系での衛星軌道を表す方程式を導出する。そして、この方程式の解を計算し、観測地点中心の座標系に変換する処理を各時刻(例、1分間隔で24時間)ごとに繰り返す。その結果、時間単位で各衛星座標が出力される。ただし、観測地点中心の座標系において、以降では、仰角θのみ使用することから、仰角θのみ計算する。通常の衛星座標の算出であれば、仰角θの他、方位角を計算するところ、方位角の計算を行わないこととしているため、計算の処理負荷をより軽減することができる。
【0021】
近似開空間算出部13は、複数の近傍点55のそれぞれについて、天空画像を取得し、天空画像から開空間65を識別し、開空間65に基づいて天頂を中心とする円形の近似開空間67を算出する。天空画像、開空間65及び近似開空間67については図5のフローチャートの説明の際に詳しく説明する。可視衛星数算出部14は、代表点51の仰角θ、及び複数の近傍点55それぞれの近似開空間67に基づいて、複数の近傍点55それぞれの可視衛星数を算出する。可視衛星数の算出については図5のフローチャートの説明の際に詳しく説明する。可視衛星数表示部15は、複数の領域54それぞれにおける可視衛星数を示した画像を表示する。ここで測位精度情報算出装置10は、可視衛星数表示部15を有しない構成であってもよい。
【0022】
図4は、天空画像取得装置22における天空画像の取得の例について説明するための図である。この図に示されるように、天空画像取得装置22、並びに三次元地図情報記憶装置23、車載カメラ画像記憶装置24及び天空撮影画像記憶装置25は通信ネットワーク32に接続されている。ここで、通信ネットワーク32は、通信ネットワーク31と同じネットワークであってもよく、また、三次元地図情報記憶装置23、車載カメラ画像記憶装置24及び天空撮影画像記憶装置25のいずれか又はすべては、通信ネットワーク32を介してアクセスされるサーバ装置として構成されていてもよい。
【0023】
ここで、三次元地図情報記憶装置23は、建物等の構造物の情報を含んだ三次元地図情報を記憶している。天空画像取得装置22は、三次元地図情報記憶装置23から三次元地図情報を取得し、取得した三次元地図情報に基づいて、位置情報が示す位置における天空画像を生成することができる。
【0024】
車載カメラ画像記憶装置24は、車載カメラ等により撮影された画像を蓄積して記憶している。天空画像取得装置22は、車載カメラ画像記憶装置24から車載カメラ等により撮影された画像(例えば中央に道路があり、両側に建物が写っている画像)を取得し、位置情報が示す位置における天空画像を生成することができる。車載カメラ等により撮影された画像は、すべての位置における画像ではないため、位置情報によっては三次元画像が生成できない場合がある。しかしながら、天空画像そのもの又は天空画像に近い画像が併せて記憶されている場合があるため、車載カメラ等により撮影された画像を利用できる場合には、三次元地図情報から天空画像を生成する場合よりも、処理負荷をも軽減することができることがある。
【0025】
天空撮影画像記憶装置25は、直接撮影された天空画像を記憶している。天空画像取得装置22は、天空撮影画像記憶装置25から天空画像を取得することができる。ある位置情報が示す位置において、直接撮影された天空画像が存在する場合は、車載カメラ等により撮影された画像が利用できる場合より少ないと考えられるが、位置情報に係る位置に天空画像が存在する場合には、ほぼそのまま天空画像として使用することができるため、処理負荷を軽減することができる。
【0026】
上述の三次元地図情報、車載カメラ等により撮影された画像の情報、及び直接撮影された天空画像の情報は、更新頻度、コスト、精度、カバーエリアに違いがあり、測位環境を把握、評価するための情報を作成する際の要件により適切なものは異なる。近似開空間算出部13は、三次元地図情報から生成された画像、車載カメラの画像から生成された画像、及び天空を撮影した画像のうちいずれか利用可能なものを、天空画像取得装置22から取得することとしてもよい。また、近似開空間算出部13は、複数の天空画像が取得可能である場合に、最新の天空画像を天空画像取得装置22から取得することとしてもよい。
【0027】
なお、図4では、三次元地図情報記憶装置23、車載カメラ画像記憶装置24及び天空撮影画像記憶装置25は、通信ネットワーク32を介して天空画像取得装置22に接続されるものとしたが、三次元地図情報記憶装置23、車載カメラ画像記憶装置24及び天空撮影画像記憶装置25のいずれか又はすべては、天空画像取得装置22内の機能により実現されるものであってもよい。
【0028】
図5は、測位精度情報算出装置10による測位精度情報算出処理S1の例について示すフローチャートである。このフローチャートに示すように、測位精度情報算出装置10の代表点決定部11は、ステップS10において、複数の領域54のそれぞれの位置を表す複数のメッシュ点52のうち、一のメッシュ点52を代表点51とし、代表点51を含み、代表点51の近傍の複数のメッシュ点52を複数の近傍点55として定める。ここで、ステップS10において、代表点決定部11は、更に、対象地域をメッシュ状に分割し、近傍点55のグループが異なる場合には、異なる代表点51を定めることとしてもよい。代表点51及び近傍点55の例については、図2で説明したものと同様である。
【0029】
次にステップS20において、ステップS10で複数の代表点51を定めた場合には、未選択の代表点51を選択し、ステップS30において、衛星座標算出部12が、選択された代表点51に関して衛星座標算出処理S30を行う。図6は、衛星座標算出処理S30を詳細に示すフローチャートである。このフローチャートに示されるように、衛星座標算出処理S30では、ステップS31において、選択された代表点51の位置情報を取得する。位置情報は代表点51における緯度経度の情報とすることができる。次にステップS32において、例えば衛星情報記憶装置21から、衛星61の軌道情報を取得する。ステップS33において、代表点51における対象時刻の衛星位置の仰角θを算出する。仰角θについては図3を用いて説明した通りである。ここで複数の衛星61についてそれぞれ仰角θを算出することができる。
【0030】
図5に戻り、衛星座標算出処理S30が終了すると、ステップS40において未選択の近傍点55を選択する。ここでの近傍点55には代表点51も含まれる。近傍点55が選択されると、引き続き、近似開空間算出部13が、近似開空間算出処理S50を行う。図7は、近似開空間算出処理S50を詳細に示すフローチャートである。このフローチャートに示されるように、近似開空間算出処理S50では、ステップS51において、近傍点55のそれぞれについて、例えば、天空画像取得装置22から天空画像を取得する。図4を用いて説明した通り、天空画像取得装置22は、三次元地図情報記憶装置23、車載カメラ画像記憶装置24及び天空撮影画像記憶装置25の情報、並びにその他の天空画像を生成できる情報に基づいて天空画像を取得又は生成することができる。
【0031】
次にステップS52において、近似開空間算出部13は開空間65を識別する。ここで開空間65は、天空画像において衛星61が可視である領域であり、閉空間66は、天空画像において、建物等の構造物の遮蔽により、衛星61が不可視である領域である。図8は、天空画像を画像処理し、開空間65及び閉空間66を識別しやすくした画像の例である。画像処理には、例えば二値化処理、エッジ抽出処理等を用いることができる。このような処理により開空間65と閉空間66とを識別することができる。
【0032】
ステップS53において、近似開空間算出部13は開空間65に基づいて天頂を中心とする円形の近似開空間67を算出する。近似開空間67の算出は、例えば、取得した天空画像上の開空間65内で、天頂を中心とする同心円の面積が最大となる仰角αを算出するものとすることができる。近似開空間67は、算出された仰角αの同心円状領域とすることができる。ここで、近似開空間67を表す値としてその観測地点における仰角αが出力されることとしてもよい。仰角αの算出方法の例としては、まず、天空画像上に矩形を想定し矩形上の周上を走査して、画像中心と距離が最小となる構造物の点を探す。ここで、周上に構造物の点があれば、より内側の矩形で同様の走査をおこない、周上に構造物の点がなければ、より外側の矩形で同様の走査をおこなう。これを繰り返すことにより、仰角αを求めることとしてもよい。
【0033】
図9は近似開空間67の例を説明するための図である。この図に示されるように、近似開空間67は、天頂を中心とする円形状であり、例えば閉空間66を含まないものとすることができる。図10は近似開空間67の別の例を説明するための図である。この図に示されるように、近似開空間67は、天頂を中心とする円形状であり、例えば閉空間66を5%以内等として一部に含むものとすることができる。
【0034】
図5に戻り、近似開空間算出処理S50が終了すると、可視衛星数算出部14が可視衛星数算出処理S60を行う。図11は、可視衛星数算出処理S60を詳細に示すフローチャートである。このフローチャートに示されるように、可視衛星数算出処理S60では、ステップS61において、仰角θを計算した衛星61のうち、未選択の衛星61を選択する。次にステップS62において選択された衛星61が可視衛星かどうかを判定する。図12は、可視衛星かどうかの判定方法の例について説明するための図である。この図に示されるように、例えば選択された近傍点55を頂点として近似開空間67が形成する円錐形状と水平面が形成する仰角αと、選択された衛星61の仰角θを比較することで可視衛星かどうかを判定することとしてもよい。図12の例では、衛星61aの仰角θaは角度αよりも大きいため、衛星61aは可視と判定される。一方、衛星61bの仰角θbは角度αよりも小さいため、衛星61bは不可視と判定される。このように、可視/不可視の判定は、従来の二次元の情報で判定していたものを、一次元の仰角の情報のみで判定しているため、計算の処理負荷をより軽減することができる。
【0035】
図11に戻り、ステップS62において可視衛星として判定された場合には、ステップS63において可視衛星数に1を加算して、ステップS64に移行する。一方ステップS62において可視衛星として判定されなかった場合には、ステップS64に移行する。ステップS64では、すべての衛星61について選択したかどうかを判定し、未だ選択されていない衛星61がある場合にはステップS61に戻り、処理を繰り返す。すべての衛星61について選択された場合には、可視衛星数算出処理S60を終了する。
【0036】
図5に戻り、可視衛星数算出処理S60を終了すると、ステップS70において、すべての近傍点55について選択したかどうかを判定する。ここで未だ選択されていない近傍点55がある場合にはステップS40に戻り、処理を繰り返す。ここで各衛星61の衛星座標(仰角θ)については、近傍点55の数の処理を繰り返さないこととしているため、計算の処理負荷をより軽減することができる。すべての近傍点55について選択された場合には、ステップS80に移行し、すべての代表点51について選択したかどうかを判定する。ここで未だ選択されていない代表点51がある場合にはステップS20に戻り、処理を繰り返す。すべての代表点51について選択されている場合には、可視衛星数表示処理S90に移行する。
【0037】
可視衛星数表示処理S90では、例えば複数の近傍点領域50からなる対象地域について、複数の領域54それぞれにおける可視衛星数を表した対象地域の画像を表示することとしてもよい。これにより、対象地域における測位精度情報を視覚的に認識することができる。図13は、可視衛星数を示した対象地域の画像の例について示す図である。複数の領域54をそれぞれ可視衛星数によって色分け等をして表示することにより、更に測位精度情報を視覚的に認識することができる。ここで、可視衛星数表示処理S70は行われなくてもよいし、リスト等として表示される別の態様であってもよい。
【0038】
以上説明したように、本実施形態の測位精度情報算出装置10では、各衛星61の仰角θのみを用い、近似開空間67内に含まれるかどうかのみで、可視衛星であるかどうかを判定しているため、方位角についての計算を必要とせず、処理負荷を軽減することができる。
【0039】
また、各衛星61の衛星座標(例えば仰角θ)は、代表点51及び近傍点55において、ほとんど変化がないため、本実施形態の測位精度情報算出装置10では、代表点51を除く近傍点55の衛星座標の計算を繰り返さないこととし、代表点51の衛星座標を用いることとしているため、計算の処理負荷をより軽減することができる。
【0040】
また、衛星数の可視/不可視の判定に関し、本実施形態の測位精度情報算出装置10では、従来の二次元の情報で可視/不可視の判定をしていたものを、一次元の仰角の情報のみで判定するようにしているため、計算の処理負荷をより軽減することができる。
【0041】
また、本実施形態の測位精度情報算出装置10では、天空画像について、実際に撮影された天空画像、三次元地図情報から生成された天空画像、車載カメラから撮影された画像から生成された天空画像その他の情報から生成された天空画像を利用することができるため、より多くの地点において、より精度の高い測位精度情報を提供することができる。
【0042】
したがって、本実施形態によれば、より処理負荷を軽減して航法衛星の測位精度情報を出力することができる。
【符号の説明】
【0043】
1…測位精度情報算出システム
10…測位精度情報算出装置
11…代表点決定部
12…衛星座標算出部
13…近似開空間算出部
14…可視衛星数算出部
15…可視衛星数表示部
21…衛星情報記憶装置
22…天空画像取得装置
23…三次元地図情報記憶装置
24…車載カメラ画像記憶装置
25…天空撮影画像記憶装置
図1
図2
図3
図4
図5
図6
図7
図8
図9
図10
図11
図12
図13