移動量計測装置及び位置計測装置
【課題】簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができるようにする。
【解決手段】GPSから、第1時刻及び第2時刻の各々で測定された測位位置と、第1時刻及び第2時刻の各々における受信信号の衛星番号群とを取得し(140)、第1時刻における衛星番号群と、第2時刻における衛星番号群とが同一である場合には(144)、第1時刻及び第2時刻で測定された測位位置の変化に基づいて、自車両の車速を算出する(146)。第1時刻と第2時刻とにおいて受信信号の衛星番号群が同一でない場合には(144)、過去の衛星番号群から、第2時刻における衛星番号群と同一の衛星番号群である時刻を検索し(148)、検索された時刻に測定された測位位置と、第2時刻における測位位置とに基づいて、自車両の車速を算出する(152)。
【解決手段】GPSから、第1時刻及び第2時刻の各々で測定された測位位置と、第1時刻及び第2時刻の各々における受信信号の衛星番号群とを取得し(140)、第1時刻における衛星番号群と、第2時刻における衛星番号群とが同一である場合には(144)、第1時刻及び第2時刻で測定された測位位置の変化に基づいて、自車両の車速を算出する(146)。第1時刻と第2時刻とにおいて受信信号の衛星番号群が同一でない場合には(144)、過去の衛星番号群から、第2時刻における衛星番号群と同一の衛星番号群である時刻を検索し(148)、検索された時刻に測定された測位位置と、第2時刻における測位位置とに基づいて、自車両の車速を算出する(152)。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、移動量計測装置及び位置計測装置に係り、特に、移動体の移動量を計測する移動量計測装置、及び移動体の移動量に基づいて移動体の位置を計測する位置計測装置に関する。
【背景技術】
【0002】
従来より、画像のオプティカルフローを検出して、移動体の運動及び移動量を推定する自己(車)運動認識システムが知られている(例えば、特許文献1)。このシステムでは、ステレオカメラを用いて、移動体の運動と移動量とを推定する。遠方用のステレオカメラと近傍用のステレオカメラとを用意し、遠方用のステレオカメラを用いて運動の回転を求め、近傍用のステレオカメラを用いて移動量としての並進を求める。
【特許文献1】特開2003−247824
【発明の開示】
【発明が解決しようとする課題】
【0003】
しかしながら、上記の特許文献1に記載の技術では、4台のカメラが必要になるため、装置にコストがかかると共に、画像処理に時間がかかる、という問題がある。また、通常のカメラを用いている限りは、周辺の照明環境の影響を受けやすいため、ロバスト性が得られない、という問題がある。
【0004】
本発明は、上記の問題点を解決するためになされたもので、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる移動量計測装置及び位置計測装置を提供することを目的とする。
【課題を解決するための手段】
【0005】
上記の目的を達成するために第1の発明に係る移動量計測装置は、複数の衛星からなる衛星群から複数の信号を受信し、受信した複数の信号に基づいて、移動体の位置を測定する測位手段と、前記測位手段によって同一の衛星群から受信した複数の信号に基づいて測位された前記移動体の複数の位置に基づいて、前記移動体の移動量を算出する移動量算出手段とを含んで構成されている。
【0006】
第1の発明に係る移動量計測装置によれば、測位手段によって、複数の衛星からなる衛星群から複数の信号を受信し、受信した複数の信号に基づいて、移動体の位置を測定する。そして、移動量算出手段によって、測位手段によって同一の衛星群から受信した複数の信号に基づいて測位された移動体の複数の位置に基づいて、移動体の移動量を算出する。
【0007】
このように、同一の衛星群からの信号に基づいて測位された移動体の複数の位置に基づいて、移動体の移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0008】
上記の第1の発明に係る移動量計測装置は、移動体の外部を撮像する単眼の撮像手段を更に含み、移動量算出手段は、第1時刻と該第1時刻より後の第2時刻とにおいて同一の衛星群から複数の信号を受信した場合、測位手段によって第1時刻に測位された移動体の位置と、測位手段によって第2時刻に測位された移動体の位置とに基づいて、第1時刻から第2時刻までの期間における移動体の移動量を算出し、第1時刻と前記第2時刻とにおいて異なる衛星群から複数の信号を受信した場合、撮像手段によって第1時刻に撮像された画像と第2時刻に撮像された画像とに基づいて、第1時刻から第2時刻までの期間における移動体の移動量を算出することができる。これによって、異なる衛星群から信号を受信した場合には、単眼の撮像手段によって撮像された画像に基づいて移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0009】
上記の第1の発明に係る移動量計測装置は、移動体の速度に応じてパルス信号を出力する速度センサを更に含み、移動量算出手段は、第1時刻と該第1時刻より後の第2時刻とにおいて同一の衛星群から複数の信号を受信した場合、測位手段によって第1時刻に測位された移動体の位置と、測位手段によって第2時刻に測位された移動体の位置とに基づいて、第1時刻から第2時刻までの期間における移動体の移動量を算出し、第1時刻と前記第2時刻とにおいて異なる衛星群から複数の信号を受信した場合、第1時刻から第2時刻までの期間に出力されたパルス信号に基づいて、第1時刻から第2時刻までの期間における移動体の移動量を算出することができる。これによって、異なる衛星群から信号を受信した場合には、速度センサのパルス信号に基づいて移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0010】
上記の第1の発明に係る位置計測装置は、移動体の外部を撮像する単眼の撮像手段と、移動体の速度に応じてパルス信号を出力する速度センサとを更に含み、移動量算出手段は、第1時刻と該第1時刻より後の第2時刻とにおいて同一の衛星群から複数の信号を受信した場合、測位手段によって第1時刻に測位された移動体の位置と、測位手段によって第2時刻に測位された移動体の位置とに基づいて、第1時刻から第2時刻までの期間における移動体の移動量を算出し、第1時刻と前記第2時刻とにおいて異なる衛星群から複数の信号を受信した場合であって、かつ、第1時刻から第2時刻までの期間に出力されたパルス信号のパルス数が、所定数以上である場合、第1時刻から第2時刻までの期間に出力されたパルス信号に基づいて、移動体の移動量を算出し、第1時刻と第2時刻とにおいて異なる衛星群から複数の信号を受信した場合であって、かつ、第1時刻から第2時刻までの期間に出力されたパルス信号のパルス数が、所定数未満である場合、撮像手段によって第1時刻に撮像された画像と第2時刻に撮像された画像とに基づいて、第1時刻から第2時刻までの期間における移動体の移動量を算出することができる。これによって、異なる衛星群から信号を受信した場合であって、かつ、パルス信号のパルス数が大きい場合には、速度センサのパルス信号に基づいて移動量を算出し、異なる衛星群から信号を受信した場合であって、かつ、パルス信号のパルス数が小さい場合には、単眼の撮像手段によって撮像された画像に基づいて移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0011】
上記の第1の発明に係る測位手段は、衛星群の各衛星から衛星を識別するための識別情報を含む信号を複数受信し、移動量算出手段は、測位手段によって識別情報の組み合わせが同一の衛星群から受信した複数の信号に基づいて測位された移動体の複数の位置に基づいて、移動体の移動量を算出することができる。このように、衛星からの信号に含まれる識別情報によって、同一の衛星群から信号を受信したか否かを判断することができる。
【0012】
第2の発明に係る移動量計測装置は、移動体の外部を撮像する単眼の撮像手段と、前記移動体の速度に応じてパルス信号を出力する速度センサと、第1時刻から該第1時刻より後の第2時刻までの期間に出力されたパルス信号のパルス数が、所定数以上である場合、前記第1時刻から前記第2時刻までの期間に出力されたパルス信号に基づいて、前記移動体の移動量を算出し、前記第1時刻から前記第2時刻までの期間に出力されたパルス信号のパルス数が、所定数未満である場合、前記撮像手段によって前記第1時刻に撮像された画像と前記第2時刻に撮像された画像とに基づいて、前記第1時刻から前記第2時刻までの期間における前記移動体の移動量を算出する移動量算出手段とを含んで構成されている。
【0013】
第2の発明に係る移動量計測装置によれば、単眼の撮像手段によって、移動体の外部を撮像し、速度センサによって、移動体の速度に応じてパルス信号を出力する。そして、移動量算出手段によって、第1時刻から該第1時刻より後の第2時刻までの期間に出力されたパルス信号のパルス数が、所定数以上である場合、第1時刻から第2時刻までの期間に出力されたパルス信号に基づいて、移動体の移動量を算出する。また、移動量算出手段によって、第1時刻から第2時刻までの期間に出力されたパルス信号のパルス数が、所定数未満である場合、撮像手段によって第1時刻に撮像された画像と第2時刻に撮像された画像とに基づいて、第1時刻から第2時刻までの期間における移動体の移動量を算出する。
【0014】
このように、パルス信号のパルス数が大きい場合には、速度センサのパルス信号に基づいて移動量を算出し、パルス信号のパルス数が小さい場合には、単眼の撮像手段によって撮像された画像に基づいて移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0015】
第3の発明に係る移動量計測装置は、複数の衛星からなる衛星群から複数の信号を受信し、受信した複数の信号に基づいて、移動体の位置を測定する測位手段と、前記移動体の外部を撮像する単眼の撮像手段と、前記撮像手段によって撮像された複数の画像に基づいて、前記画像の特徴点までの距離を算出する距離算出手段と、前記移動体の速度に応じてパルス信号を出力する速度センサと、前記測位手段によって測位された前記移動体の位置、前記距離算出手段によって算出された前記距離、及び前記速度センサから出力されたパルス信号に基づいて、前記移動体の移動量を算出する移動量算出手段とを含んで構成されている。
【0016】
第3の発明に係る移動量計測装置によれば、測位手段によって、複数の衛星からなる衛星群から複数の信号を受信し、受信した複数の信号に基づいて、移動体の位置を測定する。単眼の撮像手段によって、移動体の外部を撮像し、距離算出手段によって、撮像手段によって撮像された複数の画像に基づいて、画像の特徴点までの距離を算出する。速度センサによって、移動体の速度に応じてパルス信号を出力する。
【0017】
そして、移動量算出手段によって、測位手段によって測位された移動体の位置、距離算出手段によって算出された距離、及び速度センサから出力されたパルス信号に基づいて、移動体の移動量を算出する。
【0018】
このように、衛星群から受信した複数の信号に基づく移動体の位置、複数の画像に基づく距離、及び速度センサのパルス信号に基づいて、移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0019】
第4の発明に係る位置計測装置は、上記の移動量計測装置と、前記移動体の運動の3軸角速度の少なくとも1つを推定する運動推定手段と、前記移動量計測装置によって算出された前記移動体の移動量と、前記運動推定手段によって推定された前記移動体の運動の前記3軸角速度の少なくとも1つとに基づいて、前記移動体の位置を計測する位置計測手段とを含んで構成されている。
【0020】
本発明に係る位置計測装置によれば、移動量計測装置によって、移動体の移動量を算出し、また、運動推定手段によって、移動体の運動の3軸角速度の少なくとも1つを推定する。
【0021】
そして、位置計測手段によって、移動量計測装置によって算出された移動体の移動量と、運動推定手段によって推定された移動体の運動の3軸角速度の少なくとも1つとに基づいて、移動体の位置を算出する。
【0022】
このように、ロバスト性を確保して精度よく計測された移動体の移動量と、推定された移動体の運動とに基づいて、移動量の位置を算出することにより、簡易な構成で、ロバスト性を確保して、移動体の位置を精度よく計測することができる。
【0023】
また、上記の運動推定手段は、移動体の運動の3軸角速度の少なくとも1つ及び並進方向を推定し、位置計測手段は、移動量計測装置によって算出された移動体の移動量と、運動推定手段によって推定された移動体の運動の3軸角速度の少なくとも1つ及び並進方向とに基づいて、移動体の位置を計測することができる。
【発明の効果】
【0024】
以上説明したように、本発明の移動量計測装置によれば、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる、という効果が得られる。
【0025】
本発明の位置計測装置によれば、ロバスト性を確保して精度よく計測された移動体の移動量と、推定された移動体の運動とに基づいて、移動量の位置を算出することにより、簡易な構成で、ロバスト性を確保して、移動体の位置を精度よく計測することができる、という効果が得られる。
【発明を実施するための最良の形態】
【0026】
以下、図面を参照して本発明の実施の形態を詳細に説明する。なお、車両に搭載された位置計測装置に本発明を適用した場合を例に説明する。
【0027】
図1に示すように、第1の実施の形態に係る位置計測装置10は、自車両の運動を計測する運動計測部12と、自車両の移動量を計測する移動量計測部14と、計測された自車両の運動及び移動量に基づいて、自車両の位置を推定する位置推定部16とを備えている。
【0028】
運動計測部12は、自車両の前方の画像を撮像する単眼のカメラで構成される画像撮像部20と、撮像画像に基づいて、自車両の運動を推定する運動推定部24とを備えている。
【0029】
運動推定部24では、自車両の運動として、自車両の3軸角速度、あるいは3軸角速度及び並進方向を示す成分を推定する。なお、本実施の形態では、運動推定部24によって、並進量(並進成分スカラ量)は推定しない。
【0030】
本実施の形態では、単眼のカメラを自車両の前方に設置して、自車両の前方の画像を取得する場合を例に説明するが、自車両の外部の画像を撮像すればよく、例えば、単眼のカメラを自車両の後方に設置し、自車両の後方を撮像してもよい。また、単眼のカメラは、通常の画角40度程度のカメラでもよいし、広角のカメラでもよいし、全方位カメラでもよい。また、撮像する画像の波長について、種類は問わない。運動量を適切に推定できる場所であれば、設置場所、画角、波長、及び設置個数について問わない。
【0031】
運動推定部24は、画像撮像部20によって撮像された複数の撮像画像に基づいて、3軸角速度、あるいは3軸角速度及び並進方向を示す成分を算出する。
【0032】
移動量計測部14は、複数の衛星からなる衛星群から複数の信号を受信し、受信した複数の信号に基づいて、自車両の位置を測定するGPS26と、GPS26によって測定された自車両の位置から移動量を推定する移動量推定部28とを備えている。移動量推定部28は、GPS26による測位位置から、自車両の移動量として並進成分スカラ量又は車速を推定する。例えば、第1時刻と第1時刻から計測単位期間Tms経過した第2時刻との各々の測位位置に基づいて、計測単位期間Tmsの車速を推定する。
【0033】
ここで、本発明の原理について説明する。GPS26によって測定される測位位置は、ランダム誤差とオフセット誤差とを有している。ランダム誤差は、時々刻々と変化するものであるが、現在では、SA(Selective Availability)という規制が解除され、かつ、受信機の性能が向上しているため、安定してGPS信号が得られる環境下にいる場合には、ランダム誤差の変動分が比較的小さい。
【0034】
しかし、GPS26によって測定される測位位置は、長周期の誤差も有しているため、短時間で見るとオフセットがあるものと考えられる。この長周期の誤差は、GPS信号が伝播する際に通過する、大気圏や電離層における遅延モデルで生じる誤差によるところが大きい。この誤差は、大気圏や電離層の状態変化や衛星の配置によって変化するため、長周期の誤差となる。
【0035】
同一の衛星群から信号を受信し続けている場合は、時系列間において同一のオフセットが生じていると考えられるが、受信している複数の信号を送信している衛星群に変化があった場合には、オフセットの値に変化が生じる。
【0036】
受信している複数の信号を送信している衛星群が変化しないときに、その測位位置の変化から移動量を算出した場合は、実際の移動量とよく一致した値が得られる。しかし、受信信号の衛星群が変化した場合には、衛星群が変化した前後で、オフセット量やその正負が変化するため、移動量を算出すると誤差が大きくなる。
【0037】
図2に示すように、受信信号の衛星群に変化がない場合には、GPSの測位位置から算出した移動体の速度の時間変化と、移動体の実際の速度の時間変化とは、よく一致している。一方、受信信号の衛星群に変化があった場合には、衛星群が変化した前後で、GPSの測位位置から算出した移動体の速度の誤差が大きくなっている。
【0038】
上記図2に示すようなGPSの性質を考慮して、移動量推定部28では、以下に説明するように、GPS26の測位位置を用いて自車両の移動量を推定する。
【0039】
まず、GPS26によって測定された位置座標と、位置の測定に利用した複数の受信信号の各々に含まれる衛星の識別情報としての衛星番号を保持し、移動量を算出したい時刻間において、受信信号の衛星番号群を比較する。同一の衛星番号群である場合には、同一のオフセットが生じているものとして、移動量を算出したい時刻間の測位位置に基づいて、移動量を算出する。
【0040】
一方、移動量を算出したい時刻間において、衛星番号群に変化がある場合には、過去の測位結果を保持したデータの中から、同一の衛星番号群の信号を受信した時刻を探索する。該当時刻が存在する場合には、現時刻の測位位置と過去の該当時刻の測位位置とに基づいて、移動量を算出し、また、該当時刻から現在の一時刻前までの移動量が既知であれば、一時刻前から現時刻までの間の移動量を算出することができる。
【0041】
ただし、該当時刻と現時刻との間に、数10秒以上の時刻差がある場合、測位位置に大きな隔たりがある場合、又は、該当時刻がない場合については、受信信号の衛星番号群に差異があったとしても、測位位置の変化から移動量を算出すると共に、「低信頼」あるいは「推定誤差大」の可能性があるものとして扱う。
【0042】
位置推定部16は、計測された自車両の運動及び移動量に基づいて、自車両の相対位置の変化量を算出し、過去にGPS26で計測された絶対位置と相対位置の変化量とに基づいて、自車両の絶対位置を推定する。以下に、相対位置の変化量を算出する方法について図3を用いて説明する。
【0043】
まず、図3に示すような座標系において、初期位置を(X(0)、Y(0)、Z(0))とし、時刻tにおいて、移動量(速度V(t))、運動(3軸角速度(ヨー角速度φ(t)、ピッチ角速度θ(t)、ロール角速度ψ(t))、並進方向を示す成分(ax、ay、az))が得られると、下記の(1)式のように、時刻t+Δtにおける相対位置(X(t+Δt)、Y(t+Δt)、Z(t+Δt))を更新することができる。
【0044】
【数1】
ただし、ax/|a|、ay/|a|、及びaz/|a|は、正規化された並進方向を示す成分の各々を表わす。
【0045】
そして、自車両に搭載されたGPS26から、経度及び緯度などで表わされる絶対位置情報が過去に取得されていた場合には、絶対位置情報が示す絶対位置と、算出された相対位置とに基づいて、自車両の現在の絶対位置を計測する。なお、絶対位置情報を、必ずしもGPS26の測位位置から取得する必要はなく、例えば路側の位置が既知である設備から、通信などによって移動体の絶対位置情報を間欠的に取得してもよいし、予め位置が既知であるランドマーク情報が登録された地図などから、自車両の位置を取得してもよい。また、GPS26と同じ頻度で取得する必要もなく、移動量と運動とで補完できる間隔で取得すればよく、取得手段も問わない。
【0046】
次に、第1の実施の形態に係る位置計測装置10の作用について説明する。なお、位置計測装置10を搭載した自車両が走行しているときに、自車両の位置を計測する場合を例に説明する。
【0047】
位置計測装置10において、図4に示す位置計測処理ルーチンが実行される。まず、ステップ100において、GPS26によって自車両の位置を測位し、測位した位置を絶対位置の初期位置として設定する。
【0048】
そして、ステップ102において、運動推定処理を行い、計測単位期間Tms(例えば、第1時刻から第2時刻までの期間)の自車両の運動の3軸角速度の推定値、あるいは3軸角速度及び並進方向を示す成分の推定値を取得する。上記ステップ102の運動推定処理は、図5に示す運動推定処理ルーチンによって実現される。以下に運動推定処理ルーチンについて説明する。
【0049】
まず、ステップ120において、第1時刻及び第2時刻の各々のタイミングで画像撮像部20によって自車両前方を撮像した第1の画像及び第2の画像を取得する。そして、ステップ122において、上記ステップ120で取得した第1の画像を複数の領域に分割する。このステップ122では、例えば、予め与えられた定数y1、y2を各領域の境界線(y軸の切片座標)として、第1の画像を、画像上段の上方領域と、中段の中間領域と、下段の下方領域との3つの領域に分割する。
【0050】
次のステップ124において、上記ステップ122で得られた各分割画像中から、それぞれ所定の数だけ特徴点(不動点であることが期待される点)を抽出する。ここで、特徴点とは、周囲の点とは区別でき、追跡が容易な点のことである。このステップ124では、例えば、二次元的に濃淡変化の勾配値が大きくなる画素を検出する、一般的な周知の方法(例えば、Harrisオペレータなど)を用いて、自動的に特徴点を抽出する。
【0051】
そして、ステップ126において、上記ステップ124で抽出された第1の画像における特徴点の各々を、上記ステップ120で取得した第2の画像において追跡し、第2の画像から、対応する各特徴点をそれぞれ検出する。このとき、第1時刻と第2時刻との間で対応する特徴点は、同一の不動点であることが期待されているので、第1の画像と第2の画像とで対応する点とその周囲の点との輝度は、殆ど変化しないという仮定を用い、この仮定に基づいて各特徴点の追跡を行う。この様な追跡処理を実施する際には、例えば、周知のルーカスカナデ法(Lucas−Kanade法)などを用いると良い。また、各特徴点の追跡処理の後、第1の画像と第2の画像との間で対応する特徴点の各組の各画像内座標を、1つの適当な形式に整理する。
【0052】
そして、ステップ128において、上記ステップ126で得られた第1時刻と第2時刻との間で互いに対応する特徴点の組から、第1時刻と第2時刻との間における画像撮像部20の相対的な位置関係、即ち、画像撮像部20を搭載した自車両の第1時刻と第2時刻との間における運動に関連する情報を推定する。
【0053】
以下、図6を用いて、抽出された特徴点(不動点)と移動体の運動との関係を説明する。図6に示すように、移動体の運動は、第1の画像から第2の画像への回転行列Rと並進ベクトルtとから構成される。例えば、第1の画像上の特徴点x1に対応する第2の画像上の特徴点をx2とすると、これらの特徴点は、同一の不動点を示す点である。
【0054】
ここで、カメラの撮像特性に依る画像の歪みを補正するキャリブレーション行列K、及び基礎行列Fを用いて、対応する2つの特徴点の関係を、以下の(2)式〜(4)式で表すことができる。
【0055】
【数2】
【0056】
上記の(2)式〜(4)式のように、上記の2時刻間における8組以上の互いに対応する特徴点の組(即ち、同一の不動点)が求まっていれば、それに基づいて、上記の基礎行列Fを算出することができる。そして、カメラのキャリブレーション行列Kが既知であれば、上記の式からも分かる様に、その2時刻間における移動体の運動、即ち、その2時刻間における自車両の回転行列Rと並進ベクトルtとを求めることができる。
【0057】
そして、ステップ130において、上記ステップ128で求めた自車両の回転行列Rと並進ベクトルtとに基づいて、自車両の3軸角速度、あるいは3軸角速度及び並進方向を示す成分を算出し、算出した3軸角速度、あるいは3軸角速度及び並進方向を示す成分を、自車両の運動の推定値として出力して、運動推定処理ルーチンを終了する。
【0058】
上述したように、運動推定処理によって、自車両の運動の3軸角速度、あるいは3軸角速度及び並進方向を示す成分が取得される。
【0059】
次に、位置計測処理ルーチンのステップ104で、移動量推定処理を行い、第1時刻と第2時刻との間の自車両の移動量としての車速の推定値を取得する。上記の移動量推定処理は、図7に示す移動量推定処理ルーチンによって実現される。以下に、移動量推定処理ルーチンについて説明する。
【0060】
まず、ステップ140において、GPS26から、第1時刻及び第2時刻の各々で測定された測位位置と、第1時刻及び第2時刻の各々において測位に利用した複数の受信信号に含まれる複数の衛星番号からなる衛星番号群とを取得する。そして、ステップ142で、上記ステップ140で取得した測位位置及び衛星番号群をメモリ(図示省略)に記憶する。
【0061】
次のステップ144では、上記ステップ140で取得した第1時刻における衛星番号群と、第2時刻における衛星番号群とが同一であるか否かを判定する。上記ステップ144で、第1時刻と第2時刻とにおいて受信信号の衛星番号群が同一であると判定された場合には、第1時刻と第2時刻とにおいて同一の衛星群から信号を受信したと判断し、ステップ146において、第1時刻及び第2時刻で測定された測位位置の変化に基づいて、自車両の移動量としての車速を算出して、移動量推定処理ルーチンを終了する。
【0062】
一方、上記ステップ144で、第1時刻と第2時刻とにおいて受信信号の衛星番号群が同一でないと判定された場合には、第1時刻と第2時刻とにおいて異なる衛星群から信号を受信したと判断し、ステップ148において、メモリに記憶された過去の衛星番号群から、第2時刻における衛星番号群と同一の衛星番号群である時刻を検索する。
【0063】
そして、ステップ150において、上記ステップ148で第2時刻と同一の衛星番号群であった時刻が検索されたか否かを判定する。上記ステップ150で、第2時刻と同一の衛星番号群であった時刻が検索されたと判定された場合には、ステップ152において、メモリに記憶された、検索された時刻に測定された測位位置と、第2時刻における測位位置とに基づいて、自車両の車速を算出し、移動量推定処理ルーチンを終了する。
【0064】
一方、上記ステップ150で、第2時刻と同一の衛星番号群であった時刻が検索されなかったと判定された場合には、ステップ154において、第1時刻及び第2時刻の各々で測定された測位位置に基づいて、信頼度の低い推定値として、自車両の車速を算出して、移動量推定処理ルーチンを終了する。
【0065】
上述したように、移動量推定処理によって、自車両の移動量としての車速が取得される。
【0066】
次に、位置計測処理ルーチンのステップ106において、上記ステップ102で得られた自車両の運動の3軸角速度、あるいは3軸角速度及び並進方向を示す成分と、上記ステップ104で得られた自車両の車速とに基づいて、自車両の相対位置の変化量を算出する。次のステップ108では、上記ステップ100で設定された初期位置と、上記ステップ106で算出された相対位置の変化量とに基づいて、自車両の絶対位置を算出する。
【0067】
そして、ステップ112において、上記ステップ102からステップ110によって絶対位置を計測してから、計測単位期間Tmsが経過したか否かを判定し、計測単位期間Tmsが経過すると、ステップ102へ戻り、再び自車両の絶対位置を計測して、自車両の絶対位置を更新する。
【0068】
以上説明したように、第1の実施の形態に係る位置計測装置によれば、同一の衛星群からの信号に基づいて測位された自車両の複数の位置に基づいて、自車両の移動量を算出するため、簡易な構成で、ロバスト性を確保して、自車両の移動量を精度よく計測することができる。
【0069】
また、測位に利用した衛星群に変化がない場合の測位位置の変化を用いて、移動量を計測するため、移動量の計測精度を確保することができる。
【0070】
また、ロバスト性を確保して精度よく計測された自車両の車速と、推定された自車両の運動の3軸角速度及び並進方向を示す成分とに基づいて、自車両の位置を計測することにより、簡易な構成で、ロバスト性を確保して、自車両の位置を精度よく計測することができる。
【0071】
次に、第2の実施の形態について説明する。なお、第1の実施の形態と同様の構成となっている部分については、同一符号を付して説明を省略する。
【0072】
第2の実施の形態では、GPSによる測位位置と、車速センサからの車速パルス信号のパルス数とを組み合わせて、自車両の移動量を推定している点が第1の実施の形態と異なっている。
【0073】
図8に示すように、第2の実施の形態に係る位置計測装置210の移動量計測部214は、GPS26と、車速に応じて車速パルス信号を出力する車速パルスセンサで構成され、自車両の対地速に応じた車速パルス信号を出力する車速センサ226と、GPS26によって測定された自車両の位置、及び車速センサ226の車速パルス信号に基づいて、移動量を推定する移動量推定部228とを備えている。
【0074】
車速センサ226は、例えば、車速が60km/h(回転数が637rpm)のときに、1秒間におけるパルス数が10.6となる車速パルス信号を出力する。
【0075】
移動量推定部228は、計測単位期間Tmsにおいて、GPS26で位置を測定した時に利用した受信信号に含まれる衛星番号群に変化が無い場合、測位位置の変化から車速を算出して移動量の推定値とし、衛星番号群が変化した場合、車速センサ226からの車速パルス信号のパルス数をカウントし、パルス数から車速を算出して移動量の推定値とする。
【0076】
次に、第2の実施の形態に係る移動量推定処理ルーチンについて、図9を用いて説明する。なお、第1の実施の形態と同様の処理については、同一符号を付して説明を省略する。
【0077】
まず、ステップ140において、GPS26から、第1時刻及び第2時刻の各々で測定された測位位置と、第1時刻及び第2時刻の各々において受信した信号に含まれる衛星番号群とを取得する。
【0078】
次のステップ144では、上記ステップ140で取得した第1時刻における衛星番号群と、第2時刻における衛星番号群とが同一であるか否かを判定する。上記ステップ144で、第1時刻と第2時刻とにおける衛星番号群が同一であると判定された場合には、ステップ146において、第1時刻及び第2時刻で測定された測位位置の変化に基づいて、自車両の移動量としての車速を算出して、移動量推定処理ルーチンを終了する。
【0079】
一方、上記ステップ144で、第1時刻と第2時刻とにおける衛星番号群が同一でないと判定された場合には、ステップ250において、第1時刻から第2時刻までの間に車速センサ226から出力された車速パルス信号のパルス数を取得する。そして、ステップ252で、上記ステップ250で取得されたパルス数に基づいて、自車両の車速を算出して、移動量推定処理ルーチンを終了する。
【0080】
なお、位置計測処理ルーチン及び運動推定処理ルーチンについては、第1の実施の形態と同様であるため、説明を省略する。
【0081】
以上説明したように、第2の実施の形態に係る位置計測装置によれば、同一の衛星群からの信号に基づいて測位された自車両の複数の位置に基づいて、自車両の移動量を算出し、異なる衛星群から信号を受信した場合には、速度センサからの車速パルス信号のパルス数に基づいて移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0082】
また、速度センサの車速パルス信号を用いることで、短期的な移動量の計測精度を確保することができる。
【0083】
次に、第3の実施の形態について説明する。なお、第1の実施の形態と同様の構成となっている部分については、同一符号を付して説明を省略する。
【0084】
第3の実施の形態では、GPSによる測位位置と、撮像画像とを組み合わせて、自車両の移動量を推定している点が第1の実施の形態と異なっている。
【0085】
図10に示すように、第3の実施の形態に係る位置計測装置310は、自車両の運動を計測すると共に、自車両の移動量を計測する運動移動量計測部312と、位置推定部16とを備えている。
【0086】
運動移動量計測部312は、画像撮像部20と、運動推定部24と、GPS26と、GPS26によって測定された自車両の位置、及び画像撮像部20によって撮像された撮像画像に基づいて、移動量を推定する移動量推定部328とを備えている。
【0087】
移動量推定部328は、計測単位期間Tmsにおいて、GPS26で位置を測定した時に利用した受信信号に含まれる衛星番号群に変化が無い場合、測位位置の変化から車速を算出して移動量の推定値とし、衛星番号群が変化した場合、画像撮像部20によって撮像された複数の撮像画像に基づいて、車速を算出して移動量の推定値とする。
【0088】
移動量推定部328は、撮像された複数の撮像画像に基づいて、車速を算出する場合、撮像画像から復元した3次元空間の路面平面を推定する。そして、予め計測された、画像撮像部20の単眼カメラが設置されている路面上からの高さや単眼カメラの光軸の角度に基づいて、第1時刻及び第2時刻に撮像された撮像画像の各々から、路面上の近傍の特徴点までの距離を計測し、特徴点までの距離の変化に基づいて、車速を算出する。特に、低速時は移動量が少ないため、自車両周辺の路面を路面平面として推定してもよい。
【0089】
次に、第3の実施の形態に係る移動量推定処理ルーチンについて、図11を用いて説明する。なお、第1の実施の形態と同様の処理については、同一符号を付して、詳細な説明を省略する。
【0090】
まず、ステップ140において、GPS26から、第1時刻及び第2時刻の各々で測定された測位位置と、第1時刻及び第2時刻の各々における衛星番号群とを取得する。
【0091】
次のステップ144では、上記ステップ140で取得した第1時刻における衛星番号群と、第2時刻における衛星番号群とが同一であるか否かを判定する。上記ステップ144で、第1時刻と第2時刻とにおける衛星番号群が同一であると判定された場合には、ステップ146において、第1時刻及び第2時刻で測定された測位位置の変化に基づいて、自車両の移動量としての車速を算出して、移動量推定処理ルーチンを終了する。
【0092】
一方、上記ステップ144で、第1時刻と第2時刻とにおける衛星番号群が同一でないと判定された場合には、ステップ350において、第1時刻より前に撮像された撮像画像を画像撮像部20から取得して、取得した撮像画像に基づいて、3次元空間の路面平面を推定する。
【0093】
そして、ステップ352で、画像撮像部20から第1時刻及び第2時刻の各々で撮像された撮像画像を取得し、ステップ354において、第1時刻に撮像された撮像画像から、上記ステップ350で推定された路面平面上の特徴点までの距離を算出すると共に、第2時刻に撮像された撮像画像から、上記と同じ特徴点までの距離を算出する。
【0094】
次のステップ356では、上記ステップ354で算出された第1時刻に対する距離と第2時刻に対する距離との変化に基づいて、自車両の車速を算出して、移動量推定処理ルーチンを終了する。
【0095】
なお、他の処理については、第1の実施の形態と同様であるため、説明を省略する。
【0096】
このように、同一の衛星群からの信号に基づいて測位された自車両の複数の位置に基づいて、自車両の移動量を算出し、異なる衛星群から信号を受信した場合には、単眼のカメラによって撮像された撮像画像に基づいて移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0097】
なお、上記の実施の形態では、撮像画像から推定された路面平面上の特徴点と、画像撮像部の高さ及び角度とに基づいて、特徴点までの距離を算出する場合を例に説明したが、これに限定されるものではなく、例えば、撮像画像から、大きさが既知の物体(例えば、白線や標識)を抽出し、抽出された物体の撮像画像上の大きさに基づいて、抽出された物体までの距離を算出するようにしてもよい。この場合には、算出された大きさが既知の物体までの距離の変化から、車速を算出するようにすればよい。
【0098】
次に、第4の実施の形態について説明する。なお、第1の実施の形態及び第2の実施の形態と同様の構成となっている部分については、同一符号を付して説明を省略する。
【0099】
第4の実施の形態では、速度センサからの車速パルス信号と、撮像画像とを組み合わせて、自車両の移動量を推定している点が第1の実施の形態と異なっている。
【0100】
図12に示すように、第4の実施の形態に係る位置計測装置410は、自車両の運動を計測すると共に、自車両の移動量を計測する運動移動量計測部412と、位置推定部16と、GPS26とを備えている。
【0101】
運動移動量計測部412は、画像撮像部20と、運動推定部24と、車速センサ226と、車速センサ226から出力された車速パルス信号、及び画像撮像部20によって撮像された撮像画像に基づいて、自車両の移動量を推定する移動量推定部428とを備えている。
【0102】
ここで、本発明の原理について説明する。車速センサ226から出力される車速パルス信号は、連続値ではなく、一定の移動が行われた場合に出力される仕組みである。車速パルス信号から算出される車速の時系列変化と、実際の車速の時系列変化とを比較すると、車速パルス信号から算出される車速の値は、真値に対して大きなズレはないが、連続値ではなく離散値となる。信号処理を工夫して、推定精度を上げることは可能であるが、一般的な車速パルス信号の移動量分解能は40cm程度であるため、低速時に、移動量が得られないことがある。
【0103】
本実施の形態では、上記の車速パルス信号を出力する車速センサ226の性質を考慮し、移動量推定部428は、車速センサ226から出力される車速パルス信号を用いて、以下に説明するように、自車両の移動量としての車速を推定する。
【0104】
まず、車速センサ226から計測単位期間Tmsに出力された車速パルス信号のパルス数をカウントする。このパルス数が、たとえば1以上である場合は、パルス数から換算される車速を移動量として出力する。計測単位期間Tms内にパルス数がカウントされなかった場合(パルス数が0の場合)には、移動の有無の計測が容易ではない。この場合、撮像画像に基づいて計測される計測単位期間Tmsの距離の変化から、移動量としての車速を算出する。なお、パルス数に関する閾値を暫定的に「1」としたが、特に1に限定されるものではなく、他の適切な値を閾値として採用してもよい。
【0105】
次に、第4の実施の形態における移動量推定処理ルーチンについて、図13を用いて説明する。なお、第2の実施の形態及び第3の実施の形態と同様の処理については、同一符号を付して詳細な説明を省略する。
【0106】
まず、ステップ250において、第1時刻から第2時刻までの間に車速センサ226から出力された車速パルス信号のパルス数Npを取得する。そして、ステップ450において、上記ステップ250で取得されたパルス数Npが、1以上であるか否かを判定する。上記ステップ450で、パルス数Npが1以上であると判定された場合には、ステップ252で、上記ステップ250で取得されたパルス数Npに基づいて、自車両の車速を算出して、移動量推定処理ルーチンを終了する。
【0107】
一方、上記ステップ450で、パルス数Npが0であると判定された場合には、ステップ350において、第1時刻より前に撮像された画像に基づいて、3次元空間の路面平面を推定する。そして、ステップ352で、画像撮像部20から第1時刻及び第2時刻の各々で撮像された撮像画像を取得し、ステップ354において、第1時刻及び第2時刻の各々に撮像された撮像画像から、上記ステップ350で推定された路面平面上の特徴点までの距離を各々算出する。
【0108】
次のステップ356では、上記ステップ354で算出された第1時刻に対する距離と第2時刻に対する距離との変化に基づいて、自車両の車速を算出して、移動量推定処理ルーチンを終了する。
【0109】
なお、他の処理については、第1の実施の形態と同様であるため、説明を省略する。
【0110】
以上説明したように、第4の実施の形態に係る位置計測装置によれば、計測単位期間における車速パルス信号のパルス数が1以上である場合には、車速パルス信号のパルス数に基づいて移動量を算出し、計測単位期間における車速パルス信号のパルス数が0である場合には、単眼のカメラによって撮像された撮像画像に基づいて移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0111】
また、自車両が低速時には、低速時の計測精度が高くない車速パルス信号の代わりに、撮像画像から推定した近傍の路面平面上の特徴点までの距離の変化に基づいて、移動量を算出することにより、停止判定も可能になる。
【0112】
次に、第5の実施の形態について説明する。なお、第1の実施の形態及び第4の実施の形態と同様の構成となっている部分については、同一符号を付して説明を省略する。
【0113】
第5の実施の形態では、GPSによる測位位置と、撮像画像と、車速センサからの車速パルス信号とを組み合わせて、自車両の移動量を推定している点が、第4の実施の形態と異なっている。
【0114】
図14に示すように、第5の実施の形態に係る位置計測装置510の運動移動量計測部512は、画像撮像部20と、運動推定部24と、GPS26と、車速センサ226と、GPS26によって測定された測位位置、車速センサ226から出力された車速パルス信号、及び画像撮像部20によって撮像された撮像画像に基づいて、移動量を推定する移動量推定部528とを備えている。
【0115】
移動量推定部528は、計測単位期間Tmsにおいて、GPS26で位置を測定した時に利用した受信信号に含まれる衛星番号群に変化が無い場合、測位位置の変化から車速を算出して移動量の推定値とする。一方、衛星番号群が変化した場合、車速センサ226から計測単位期間Tmsに出力された車速パルス信号のパルス数をカウントし、このパルス数が、たとえば1以上である場合は、車速センサ226の車速パルス信号から車速を算出して移動量の推定値とする。計測単位期間Tms内にパルス数がカウントされなかった場合には、撮像画像に基づいて計測される計測単位期間Tmsの距離の変化から、移動量としての車速を算出する。
【0116】
次に、第5の実施の形態に係る移動量推定処理ルーチンについて、図15を用いて説明する。なお、第1の実施の形態〜第4の実施の形態と同様の処理については、同一符号を付して詳細な説明を省略する。
【0117】
まず、ステップ140において、GPS26から、第1時刻及び第2時刻の各々で測定された測位位置と、第1時刻及び第2時刻の各々における衛星番号群とを取得する。
【0118】
次のステップ144では、上記ステップ140で取得した第1時刻における衛星番号群と、第2時刻における衛星番号群とが同一であるか否かを判定する。上記ステップ144で、第1時刻と第2時刻とにおける衛星番号群が同一であると判定された場合には、ステップ146において、第1時刻及び第2時刻で測定された測位位置の変化に基づいて、自車両の移動量としての車速を算出して、移動量推定処理ルーチンを終了する。
【0119】
一方、上記ステップ144で、第1時刻と第2時刻とにおける衛星番号群が同一でないと判定された場合には、ステップ250において、第1時刻から第2時刻までの間に車速センサ226から出力された車速パルス信号のパルス数を取得する。そして、ステップ450において、上記ステップ250で取得されたパルス数Npが、1以上であるか否かを判定する。上記ステップ450で、パルス数Npが1以上であると判定された場合には、ステップ252で、上記ステップ250で取得されたパルス数Npに基づいて、自車両の車速を算出して、移動量推定処理ルーチンを終了する。
【0120】
一方、上記ステップ450で、パルス数Npが0であると判定された場合には、ステップ350において、第1時刻より前に撮像された撮像画像に基づいて、3次元空間の路面平面を推定する。そして、ステップ352で、画像撮像部20から第1時刻及び第2時刻の各々で撮像された撮像画像を取得し、ステップ354において、第1時刻及び第2時刻の各々に撮像された撮像画像から、上記ステップ350で推定された路面平面上の特徴点までの距離を各々算出する。
【0121】
次のステップ356では、上記ステップ354で算出された第1時刻に対する距離と第2時刻に対する距離との変化に基づいて、自車両の車速を算出して、移動量推定処理ルーチンを終了する。
【0122】
このように、計測単位期間の第1時刻と第2時刻とにおいて、異なる衛星群から信号を受信した場合であって、かつ、車速パルス信号のパルス数が大きい場合には、速度センサからの車速パルス信号に基づいて移動量を算出し、異なる衛星群から信号を受信した場合であって、かつ、車速パルス信号のパルス数が小さい場合には、単眼のカメラによって撮像された撮像画像に基づいて移動量を算出するため、簡易な構成で、ロバスト性を確保して、自車両の移動量を精度よく計測することができる。
【0123】
次に、第6の実施の形態について説明する。なお、第6の実施の形態に係る位置計測装置は、第5の実施の形態と同様の構成となっているため、同一符号を付して説明を省略する。
【0124】
第6の実施の形態では、カルマンフィルタ等を用いた最適処理手法によって、自車両の移動量としての車速を推定している点が第5の実施の形態と異なっている。
【0125】
第6の実施の形態に係る位置計測装置510の移動量推定部528では、画像撮像部20の撮像画像に基づいて計測される距離と、GPS26によって測定される測位位置と、車速センサ226からの車速パルス信号のパルス数とを観測値とし、自車両の車速を推定値として、カルマンフィルタによる最適処理手法を用いた処理を行う。
【0126】
カルマンフィルタによる最適処理手法では、画像撮像部20の単眼カメラの撮像画像に基づいて計測される距離の観測値yimageが得られた場合に、観測値yimageのバラつきを考慮して、観測ノイズwimageを設定し、観測行列Himageについて、観測値yimageと推定値xとの関係式から、観測行列Himageを設定する。
【0127】
また、GPS26によって測位された測位位置の観測値ygpsが得られた場合に、観測値ygpsのバラつきを考慮して、観測ノイズwgpsを設定し、観測行列Hgpsについて、観測値ygpsと推定値xとの関係式から、観測行列Hgpsを設定する。
【0128】
また、車速センサ226から出力された車速パルス信号のパルス数の観測値ypulsが得られた場合に、観測値ypulsのバラつきを考慮して、観測ノイズwpulsを設定し、観測行列Hpulsについて、観測値ypulsと推定値xとの関係式から、観測行列Hpulsを設定する。
【0129】
そして、設定された観測行列、観測値、及び観測ノイズに基づいて、カルマンフィルタで、推定値x及び推定誤差共分散行列を算出して出力する。観測値が得られない場合には、推定値の予測値と推定誤差共分散行列の予測値とを算出して出力する。
【0130】
次に、第6の実施の形態に係る位置計測装置510の作用について説明する。なお、第1の実施の形態と同様の部分については、同一符号を付して詳細な説明を省略する。
【0131】
第6の実施の形態に係る位置計測装置510では、上述した移動量推定処理ルーチンにおけるステップ350からステップ354までの処理と同様の処理が繰り返し実行され、算出された路面平面上の特徴点までの距離が繰り返し出力される。また、GPS26によって自車両の位置が繰り返し測定される。また、車速センサ226から、車速パルス信号が連続して出力され、パルス数が繰り返しカウントされる。
【0132】
また、位置計測装置510において、図16に示す移動量推定処理ルーチンが実行される。まず、ステップ600において、撮像画像に基づいて計測される路面平面上の特徴点までの距離を取得し、ステップ602において、GPS26から測定された測位位置を取得する。そして、ステップ604において、車速センサ226からの車速パルス信号に基づいてカウントされたパルス数を取得する。
【0133】
そして、ステップ606において、上記ステップ600で取得した距離と、上記ステップ602で取得した測位位置と、上記ステップ604で取得したパルス数とに基づいて、移動量としての車速の推定値xに対して、初期値を設定する。次のステップ608では、時間ti−1にΔtを加算して、時間tiに更新する。そして、ステップ610において、1時刻前の推定値xに基づいて、推定値xの予測値を算出すると共に、推定誤差共分散の予測値を算出する。
【0134】
次のステップ612において、現タイミングに、撮像画像に基づいて計測される路面平面上の特徴点までの距離、GPS26から測定された測位位置、及び車速センサ226からの車速パルス信号のパルス数の何れかが、観測値として入力されたか否かを判定し、現タイミングに観測値が入力されなかった場合には、ステップ614において、上記ステップ610で算出された推定値xの予測値の車速を、自車両の移動量の推定値として出力し、ステップ608へ戻る。
【0135】
上記ステップ612において、現タイミングに観測値が入力されたと判定された場合には、ステップ616において、どの観測値が入力されたのかを判定する。
【0136】
上記ステップ616で、観測値としてGPS26から測定された測位位置が入力されたと判定された場合には、ステップ618において、入力された測位位置を観測値として、観測行列Hgps、観測値ygps、及び観測ノイズwgpsを設定する。
【0137】
一方、上記ステップ616で、観測値として、撮像画像に基づいて計測される路面平面上の特徴点までの距離が入力されたと判定された場合には、ステップ620において、入力された特徴点までの距離を観測値として、観測行列Himage、観測値yimage、及び観測ノイズwimageを設定する。
【0138】
また、上記ステップ616で、観測値として、車速センサ226からの車速パルス信号のパルス数が入力されたと判定された場合には、ステップ622において、入力されたパルス数を観測値として、観測行列Hpuls、観測値ypuls、及び観測ノイズwpulsを設定する。
【0139】
そして、ステップ624において、上記ステップ618、上記ステップ620、又は上記ステップ622で設定された観測行列、観測値、及び観測ノイズに基づいて、カルマンフィルタで、カルマンゲインK、推定値x、及び推定誤差共分散行列を算出する。次のステップ626では、上記ステップ624で算出された推定値xの車速を、自車両の移動量の推定値として出力し、ステップ608へ戻る。
【0140】
このように、カルマンフィルタを用いた最適処理手法によって、自車両の外部を撮像した複数の撮像画像に基づいて計測された路面平面上の特徴点までの距離と、GPSによって測定された測位位置と、車速センサからの車速パルス信号のパルス数とに基づいて、自車両の移動量の推定値を算出することにより、簡易な構成で、ロバスト性を確保して、自車両の移動量を精度よく計測することができる。
【0141】
また、ロバスト性を確保して精度よく計測された自車両の移動量と、推定された自車両の運動の3軸角速度及び並進方向を示す成分とに基づいて、自車両の位置を計測することにより、簡易な構成で、ロバスト性を確保して、自車両の位置を精度よく計測することができる。
【0142】
なお、上記の実施の形態では、カルマンフィルタで推定する例を記述したが、周知のパーティクルフィルタなど、他の最適推定手法を用いてもよく、適切に解を推定できる手法であれば、その方法はこの限りではない。
【0143】
また、誤差共分散を、推定された自車両の移動量の推定値に対する信頼度として用いてもよい。例えば、推定値と共に算出された誤差共分散によって、推定値の信頼度が低く、推定値の誤差が大きいと判断される場合に、この推定値を、自車両の移動量の推定値として採用しないようにしてもよい。
【0144】
次に、第7の実施の形態について説明する。なお、第7の実施の形態に係る位置計測装置は、第1の実施の形態と同様の構成となっているため、同一符号を付して説明を省略する。
【0145】
第7の実施の形態では、カルマンフィルタを用いた最適処理手法によって、自車両の位置を計測している点が第1の実施の形態と異なっている。
【0146】
第7の実施の形態に係る位置計測装置10の位置推定部16では、運動推定部24によって推定された運動(3軸角速度、並進方向を示す成分)、移動量推定部28によって推定された移動量(車速)、及びGPS26によって測定された絶対位置(Xg,Yg,Zg)を観測値とし、自車の絶対位置、方位角、車速、及び3軸角速度を推定値として、カルマンフィルタを用いた最適処理手法による位置計測処理を行う。
【0147】
カルマンフィルタを用いた最適処理手法では、運動推定部24によって推定された運動の観測値ymoveが得られた場合に、観測値ymoveのバラつきを考慮して、観測ノイズwmoveを設定し、観測行列Hmoveについて、観測値ymoveと推定値xとの関係式から、観測行列Hmoveを設定する。
【0148】
また、移動量推定部28によって推定された移動量の観測値ytravが得られた場合に、観測値ytravのバラつきを考慮して、観測ノイズwtravを設定し、観測行列Htravについて、観測値ytravと推定値xとの関係式から、観測行列Htravを設定する。
【0149】
また、GPS26によって測定された絶対位置の観測値ygpsが得られた場合に、観測値ygpsのバラつきを考慮して、観測ノイズwgpsを設定し、観測行列Hgpsについて、観測値ygpsと推定値xとの関係式から、観測行列Hgpsを設定する。
【0150】
そして、設定された観測行列、観測値、及び観測ノイズに基づいて、カルマンフィルタで、推定値x及び推定誤差共分散行列を算出して出力する。観測値が得られない場合には、推定値の予測値と推定誤差共分散行列の予測値とを算出して出力する。
【0151】
次に、第7の実施の形態に係る位置計測装置10の作用について説明する。なお、第1の実施の形態と同様の部分については、同一符号を付して詳細な説明を省略する。
【0152】
第7の実施の形態に係る位置計測装置10では、第1の実施の形態で説明した運動推定処理ルーチンが繰り返し実行され、また、上述した移動量推定処理ルーチンが繰り返し実行される。また、GPS26によって、所定時間ごと(例えば1秒ごと)に自車両の絶対位置が測定される。
【0153】
また、位置計測装置10において、図17に示す位置計測処理ルーチンが実行される。まず、ステップ700において、運動推定処理ルーチンによって推定された3軸角速度及び並進方向を示す成分を取得し、ステップ702において、移動量推定処理ルーチンによって推定された車速を取得する。次のステップ704では、GPS26によって測定された自車両の絶対位置を取得する。
【0154】
そして、ステップ706において、上記ステップ700で取得した3軸角速度及び並進方向を示す成分と、上記ステップ702で取得した車速と、上記ステップ704で取得した絶対位置とに基づいて、絶対位置、方位角、車速、及び3軸角速度の推定値xに対して、初期値を設定する。次のステップ708では、時間ti−1にΔtを加算して、時間tiに更新する。そして、ステップ710において、1時刻前の推定値xに基づいて、推定値xの予測値を算出すると共に、推定誤差共分散の予測値を算出する。
【0155】
次のステップ712において、現タイミングに、観測値として、運動推定処理ルーチンによって推定された3軸角速度、移動推定処理ルーチンによって推定された車速、又はGPS26によって測定された絶対位置が入力されたか否かを判定し、現タイミングに観測値が入力されなかった場合には、ステップ714において、上記ステップ710で算出された推定値xの予測値の絶対位置を、自車両の絶対位置の推定値として出力し、ステップ708へ戻る。
【0156】
上記ステップ712において、現タイミングに観測値が入力されたと判定された場合には、ステップ716において、上記ステップ712でどの観測値が入力されたかを判定する。
【0157】
上記ステップ716で、観測値として、運動推定処理ルーチンで推定された3軸角速度及び並進方向を示す成分が入力されたと判定された場合には、ステップ718において、入力された3軸角速度及び並進方向を示す成分を観測値として、観測行列Hmove、観測値ymove、及び観測ノイズwmoveを設定する。
【0158】
一方、上記ステップ716で、観測値として、移動量推定処理ルーチンによって推定された車速が入力されたと判定された場合には、ステップ720において、入力された車速を観測値として、観測行列Htrav、観測値ytrav、及び観測ノイズwtravを設定する。
【0159】
また、上記ステップ716で、観測値として、GPS26によって測定された絶対位置が入力されたと判定された場合には、ステップ722において、入力された絶対位置を観測値として、観測行列Hgps、観測値ygps、及び観測ノイズwgpsを設定する。
【0160】
そして、ステップ724において、上記ステップ718、上記ステップ720、又は上記ステップ722で設定された観測行列、観測値、及び観測ノイズに基づいて、カルマンフィルタで、カルマンゲインK、推定値x、及び推定誤差共分散行列を算出する。次のステップ726では、上記ステップ724で算出された推定値xの絶対位置を、自車両の絶対位置の推定値として出力し、ステップ708へ戻る。
【0161】
このように、カルマンフィルタを用いた最適処理手法によって、推定された自車両の運動の3軸角速度及び並進方向を示す成分と、推定された自車両の車速と、GPSによって測位された位置とに基づいて、自車両の絶対位置の推定値を算出することにより、簡易な構成で、ロバスト性を確保して、自車両の絶対位置を精度よく計測することができる。
【0162】
なお、上記の実施の形態では、カルマンフィルタを用いた位置計測処理において、推定された自車両の運動の3軸角速度及び並進方向を示す成分、推定された自車両の車速、及びGPSによって測位された絶対位置が、観測値として入力される場合を例に説明したが、これに限定されるものではない。例えば、推定された自車両の運動の3軸角速度及び並進方向を示す成分と、及び推定された自車両の車速とに基づいて推定された自車両の相対位置の変化量が、観測値として、カルマンフィルタを用いた位置計測処理に入力されるようになっていてもよい。この場合には、推定された自車両の相対位置の変化量、及びGPSによって測位された絶対位置を観測値とし、自車の絶対位置、方位角、車速、及び3軸角速度を推定値として、カルマンフィルタを用いた最適処理手法による位置計測処理を行うようにすればよい。
【0163】
また、カルマンフィルタに限らず、パーティクルフィルタなどの最適推定手法を用いてもよい。
【0164】
次に、第8の実施の形態について説明する。第8の実施の形態では、ジャイロセンサによって、自車両の運動の3軸角速度を検出している。
【0165】
第8の実施の形態に係る位置計測装置は、ジャイロセンサによって検出された自車両の運動の3軸角速度と、移動量推定部によって推定された自車両の移動量とに基づいて、自車両の位置を推定する。
【0166】
次に、第9の実施の形態について説明する。第9の実施の形態では、ジャイロセンサ及び単眼カメラによって撮像された撮像画像に基づいて、自車両の運動の3軸角速度、あるいは3軸角速度及び並進方向を示す成分を推定している。
【0167】
第9の実施の形態に係る位置計測装置は、ジャイロセンサ及び撮像画像によって推定された自車両の運動の3軸角速度、あるいは3軸角速度及び並進方向を示す成分と、移動量推定部によって推定された自車両の移動量とに基づいて、自車両の位置を推定する。
【0168】
なお、上記の第1の実施の形態〜第9の実施の形態では、撮像画像又はジャイロセンサに基づいて3軸角速度を推定する場合を例に説明したが、これに限定されるものではなく、撮像画像又はジャイロセンサに基づいて、ヨー角速度、ピッチ角速度、及びロール角速度の少なくとも1つを推定するように構成してもよい。
【0169】
また、自車両の移動量として、車速を推定する場合を例に説明したが、これに限定されるものではなく、移動量として、移動距離又は移動ベクトルを推定するようにしてもよい。
【0170】
また、位置計測処理をオンラインで逐次処理する場合を例に説明したが、これに限定されるものではなく、自車両の運動の推定値及び移動量の推定値を取得した結果を用いて、オフラインで後処理として位置計測処理を行うようにしてもよい。
【図面の簡単な説明】
【0171】
【図1】本発明の第1の実施の形態に係る位置計測装置を示すブロック図である。
【図2】GPSの測位位置に基づいて算出される車速の時間変化と実際の車速の時間変化とを示すグラフである。
【図3】相対位置の変化量の算出方法を説明するための座標系を示すイメージ図である。
【図4】本発明の第1の実施の形態に係る位置計測装置における位置計測処理ルーチンの内容を示すフローチャートである。
【図5】本発明の第1の実施の形態に係る位置計測装置における運動推定処理ルーチンの内容を示すフローチャートである。
【図6】抽出された特徴点と移動体の運動との関係を説明するための図である。
【図7】本発明の第1の実施の形態に係る位置計測装置における移動量推定処理ルーチンの内容を示すフローチャートである。
【図8】本発明の第2の実施の形態に係る位置計測装置を示すブロック図である。
【図9】本発明の第2の実施の形態に係る位置計測装置における移動量推定処理ルーチンの内容を示すフローチャートである。
【図10】本発明の第3の実施の形態に係る位置計測装置を示すブロック図である。
【図11】本発明の第3の実施の形態に係る位置計測装置における移動量推定処理ルーチンの内容を示すフローチャートである。
【図12】本発明の第4の実施の形態に係る位置計測装置を示すブロック図である。
【図13】本発明の第4の実施の形態に係る位置計測装置における移動量推定処理ルーチンの内容を示すフローチャートである。
【図14】本発明の第5の実施の形態に係る位置計測装置を示すブロック図である。
【図15】本発明の第5の実施の形態に係る位置計測装置における移動量推定処理ルーチンの内容を示すフローチャートである。
【図16】本発明の第6の実施の形態に係る位置計測装置における移動量推定処理ルーチンの内容を示すフローチャートである。
【図17】本発明の第7の実施の形態に係る位置計測装置における位置計測処理ルーチンの内容を示すフローチャートである。
【符号の説明】
【0172】
10、210、310、410、510 位置計測装置
12、312 運動計測部
14、214 移動量計測部
16 位置推定部
20 画像撮像部
24 運動推定部
26 GPS
28、228、328、428、528 移動量推定部
226 車速センサ
312、412、512 運動移動量計測部
528 移動量推定部
Np パルス数
【技術分野】
【0001】
本発明は、移動量計測装置及び位置計測装置に係り、特に、移動体の移動量を計測する移動量計測装置、及び移動体の移動量に基づいて移動体の位置を計測する位置計測装置に関する。
【背景技術】
【0002】
従来より、画像のオプティカルフローを検出して、移動体の運動及び移動量を推定する自己(車)運動認識システムが知られている(例えば、特許文献1)。このシステムでは、ステレオカメラを用いて、移動体の運動と移動量とを推定する。遠方用のステレオカメラと近傍用のステレオカメラとを用意し、遠方用のステレオカメラを用いて運動の回転を求め、近傍用のステレオカメラを用いて移動量としての並進を求める。
【特許文献1】特開2003−247824
【発明の開示】
【発明が解決しようとする課題】
【0003】
しかしながら、上記の特許文献1に記載の技術では、4台のカメラが必要になるため、装置にコストがかかると共に、画像処理に時間がかかる、という問題がある。また、通常のカメラを用いている限りは、周辺の照明環境の影響を受けやすいため、ロバスト性が得られない、という問題がある。
【0004】
本発明は、上記の問題点を解決するためになされたもので、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる移動量計測装置及び位置計測装置を提供することを目的とする。
【課題を解決するための手段】
【0005】
上記の目的を達成するために第1の発明に係る移動量計測装置は、複数の衛星からなる衛星群から複数の信号を受信し、受信した複数の信号に基づいて、移動体の位置を測定する測位手段と、前記測位手段によって同一の衛星群から受信した複数の信号に基づいて測位された前記移動体の複数の位置に基づいて、前記移動体の移動量を算出する移動量算出手段とを含んで構成されている。
【0006】
第1の発明に係る移動量計測装置によれば、測位手段によって、複数の衛星からなる衛星群から複数の信号を受信し、受信した複数の信号に基づいて、移動体の位置を測定する。そして、移動量算出手段によって、測位手段によって同一の衛星群から受信した複数の信号に基づいて測位された移動体の複数の位置に基づいて、移動体の移動量を算出する。
【0007】
このように、同一の衛星群からの信号に基づいて測位された移動体の複数の位置に基づいて、移動体の移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0008】
上記の第1の発明に係る移動量計測装置は、移動体の外部を撮像する単眼の撮像手段を更に含み、移動量算出手段は、第1時刻と該第1時刻より後の第2時刻とにおいて同一の衛星群から複数の信号を受信した場合、測位手段によって第1時刻に測位された移動体の位置と、測位手段によって第2時刻に測位された移動体の位置とに基づいて、第1時刻から第2時刻までの期間における移動体の移動量を算出し、第1時刻と前記第2時刻とにおいて異なる衛星群から複数の信号を受信した場合、撮像手段によって第1時刻に撮像された画像と第2時刻に撮像された画像とに基づいて、第1時刻から第2時刻までの期間における移動体の移動量を算出することができる。これによって、異なる衛星群から信号を受信した場合には、単眼の撮像手段によって撮像された画像に基づいて移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0009】
上記の第1の発明に係る移動量計測装置は、移動体の速度に応じてパルス信号を出力する速度センサを更に含み、移動量算出手段は、第1時刻と該第1時刻より後の第2時刻とにおいて同一の衛星群から複数の信号を受信した場合、測位手段によって第1時刻に測位された移動体の位置と、測位手段によって第2時刻に測位された移動体の位置とに基づいて、第1時刻から第2時刻までの期間における移動体の移動量を算出し、第1時刻と前記第2時刻とにおいて異なる衛星群から複数の信号を受信した場合、第1時刻から第2時刻までの期間に出力されたパルス信号に基づいて、第1時刻から第2時刻までの期間における移動体の移動量を算出することができる。これによって、異なる衛星群から信号を受信した場合には、速度センサのパルス信号に基づいて移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0010】
上記の第1の発明に係る位置計測装置は、移動体の外部を撮像する単眼の撮像手段と、移動体の速度に応じてパルス信号を出力する速度センサとを更に含み、移動量算出手段は、第1時刻と該第1時刻より後の第2時刻とにおいて同一の衛星群から複数の信号を受信した場合、測位手段によって第1時刻に測位された移動体の位置と、測位手段によって第2時刻に測位された移動体の位置とに基づいて、第1時刻から第2時刻までの期間における移動体の移動量を算出し、第1時刻と前記第2時刻とにおいて異なる衛星群から複数の信号を受信した場合であって、かつ、第1時刻から第2時刻までの期間に出力されたパルス信号のパルス数が、所定数以上である場合、第1時刻から第2時刻までの期間に出力されたパルス信号に基づいて、移動体の移動量を算出し、第1時刻と第2時刻とにおいて異なる衛星群から複数の信号を受信した場合であって、かつ、第1時刻から第2時刻までの期間に出力されたパルス信号のパルス数が、所定数未満である場合、撮像手段によって第1時刻に撮像された画像と第2時刻に撮像された画像とに基づいて、第1時刻から第2時刻までの期間における移動体の移動量を算出することができる。これによって、異なる衛星群から信号を受信した場合であって、かつ、パルス信号のパルス数が大きい場合には、速度センサのパルス信号に基づいて移動量を算出し、異なる衛星群から信号を受信した場合であって、かつ、パルス信号のパルス数が小さい場合には、単眼の撮像手段によって撮像された画像に基づいて移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0011】
上記の第1の発明に係る測位手段は、衛星群の各衛星から衛星を識別するための識別情報を含む信号を複数受信し、移動量算出手段は、測位手段によって識別情報の組み合わせが同一の衛星群から受信した複数の信号に基づいて測位された移動体の複数の位置に基づいて、移動体の移動量を算出することができる。このように、衛星からの信号に含まれる識別情報によって、同一の衛星群から信号を受信したか否かを判断することができる。
【0012】
第2の発明に係る移動量計測装置は、移動体の外部を撮像する単眼の撮像手段と、前記移動体の速度に応じてパルス信号を出力する速度センサと、第1時刻から該第1時刻より後の第2時刻までの期間に出力されたパルス信号のパルス数が、所定数以上である場合、前記第1時刻から前記第2時刻までの期間に出力されたパルス信号に基づいて、前記移動体の移動量を算出し、前記第1時刻から前記第2時刻までの期間に出力されたパルス信号のパルス数が、所定数未満である場合、前記撮像手段によって前記第1時刻に撮像された画像と前記第2時刻に撮像された画像とに基づいて、前記第1時刻から前記第2時刻までの期間における前記移動体の移動量を算出する移動量算出手段とを含んで構成されている。
【0013】
第2の発明に係る移動量計測装置によれば、単眼の撮像手段によって、移動体の外部を撮像し、速度センサによって、移動体の速度に応じてパルス信号を出力する。そして、移動量算出手段によって、第1時刻から該第1時刻より後の第2時刻までの期間に出力されたパルス信号のパルス数が、所定数以上である場合、第1時刻から第2時刻までの期間に出力されたパルス信号に基づいて、移動体の移動量を算出する。また、移動量算出手段によって、第1時刻から第2時刻までの期間に出力されたパルス信号のパルス数が、所定数未満である場合、撮像手段によって第1時刻に撮像された画像と第2時刻に撮像された画像とに基づいて、第1時刻から第2時刻までの期間における移動体の移動量を算出する。
【0014】
このように、パルス信号のパルス数が大きい場合には、速度センサのパルス信号に基づいて移動量を算出し、パルス信号のパルス数が小さい場合には、単眼の撮像手段によって撮像された画像に基づいて移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0015】
第3の発明に係る移動量計測装置は、複数の衛星からなる衛星群から複数の信号を受信し、受信した複数の信号に基づいて、移動体の位置を測定する測位手段と、前記移動体の外部を撮像する単眼の撮像手段と、前記撮像手段によって撮像された複数の画像に基づいて、前記画像の特徴点までの距離を算出する距離算出手段と、前記移動体の速度に応じてパルス信号を出力する速度センサと、前記測位手段によって測位された前記移動体の位置、前記距離算出手段によって算出された前記距離、及び前記速度センサから出力されたパルス信号に基づいて、前記移動体の移動量を算出する移動量算出手段とを含んで構成されている。
【0016】
第3の発明に係る移動量計測装置によれば、測位手段によって、複数の衛星からなる衛星群から複数の信号を受信し、受信した複数の信号に基づいて、移動体の位置を測定する。単眼の撮像手段によって、移動体の外部を撮像し、距離算出手段によって、撮像手段によって撮像された複数の画像に基づいて、画像の特徴点までの距離を算出する。速度センサによって、移動体の速度に応じてパルス信号を出力する。
【0017】
そして、移動量算出手段によって、測位手段によって測位された移動体の位置、距離算出手段によって算出された距離、及び速度センサから出力されたパルス信号に基づいて、移動体の移動量を算出する。
【0018】
このように、衛星群から受信した複数の信号に基づく移動体の位置、複数の画像に基づく距離、及び速度センサのパルス信号に基づいて、移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0019】
第4の発明に係る位置計測装置は、上記の移動量計測装置と、前記移動体の運動の3軸角速度の少なくとも1つを推定する運動推定手段と、前記移動量計測装置によって算出された前記移動体の移動量と、前記運動推定手段によって推定された前記移動体の運動の前記3軸角速度の少なくとも1つとに基づいて、前記移動体の位置を計測する位置計測手段とを含んで構成されている。
【0020】
本発明に係る位置計測装置によれば、移動量計測装置によって、移動体の移動量を算出し、また、運動推定手段によって、移動体の運動の3軸角速度の少なくとも1つを推定する。
【0021】
そして、位置計測手段によって、移動量計測装置によって算出された移動体の移動量と、運動推定手段によって推定された移動体の運動の3軸角速度の少なくとも1つとに基づいて、移動体の位置を算出する。
【0022】
このように、ロバスト性を確保して精度よく計測された移動体の移動量と、推定された移動体の運動とに基づいて、移動量の位置を算出することにより、簡易な構成で、ロバスト性を確保して、移動体の位置を精度よく計測することができる。
【0023】
また、上記の運動推定手段は、移動体の運動の3軸角速度の少なくとも1つ及び並進方向を推定し、位置計測手段は、移動量計測装置によって算出された移動体の移動量と、運動推定手段によって推定された移動体の運動の3軸角速度の少なくとも1つ及び並進方向とに基づいて、移動体の位置を計測することができる。
【発明の効果】
【0024】
以上説明したように、本発明の移動量計測装置によれば、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる、という効果が得られる。
【0025】
本発明の位置計測装置によれば、ロバスト性を確保して精度よく計測された移動体の移動量と、推定された移動体の運動とに基づいて、移動量の位置を算出することにより、簡易な構成で、ロバスト性を確保して、移動体の位置を精度よく計測することができる、という効果が得られる。
【発明を実施するための最良の形態】
【0026】
以下、図面を参照して本発明の実施の形態を詳細に説明する。なお、車両に搭載された位置計測装置に本発明を適用した場合を例に説明する。
【0027】
図1に示すように、第1の実施の形態に係る位置計測装置10は、自車両の運動を計測する運動計測部12と、自車両の移動量を計測する移動量計測部14と、計測された自車両の運動及び移動量に基づいて、自車両の位置を推定する位置推定部16とを備えている。
【0028】
運動計測部12は、自車両の前方の画像を撮像する単眼のカメラで構成される画像撮像部20と、撮像画像に基づいて、自車両の運動を推定する運動推定部24とを備えている。
【0029】
運動推定部24では、自車両の運動として、自車両の3軸角速度、あるいは3軸角速度及び並進方向を示す成分を推定する。なお、本実施の形態では、運動推定部24によって、並進量(並進成分スカラ量)は推定しない。
【0030】
本実施の形態では、単眼のカメラを自車両の前方に設置して、自車両の前方の画像を取得する場合を例に説明するが、自車両の外部の画像を撮像すればよく、例えば、単眼のカメラを自車両の後方に設置し、自車両の後方を撮像してもよい。また、単眼のカメラは、通常の画角40度程度のカメラでもよいし、広角のカメラでもよいし、全方位カメラでもよい。また、撮像する画像の波長について、種類は問わない。運動量を適切に推定できる場所であれば、設置場所、画角、波長、及び設置個数について問わない。
【0031】
運動推定部24は、画像撮像部20によって撮像された複数の撮像画像に基づいて、3軸角速度、あるいは3軸角速度及び並進方向を示す成分を算出する。
【0032】
移動量計測部14は、複数の衛星からなる衛星群から複数の信号を受信し、受信した複数の信号に基づいて、自車両の位置を測定するGPS26と、GPS26によって測定された自車両の位置から移動量を推定する移動量推定部28とを備えている。移動量推定部28は、GPS26による測位位置から、自車両の移動量として並進成分スカラ量又は車速を推定する。例えば、第1時刻と第1時刻から計測単位期間Tms経過した第2時刻との各々の測位位置に基づいて、計測単位期間Tmsの車速を推定する。
【0033】
ここで、本発明の原理について説明する。GPS26によって測定される測位位置は、ランダム誤差とオフセット誤差とを有している。ランダム誤差は、時々刻々と変化するものであるが、現在では、SA(Selective Availability)という規制が解除され、かつ、受信機の性能が向上しているため、安定してGPS信号が得られる環境下にいる場合には、ランダム誤差の変動分が比較的小さい。
【0034】
しかし、GPS26によって測定される測位位置は、長周期の誤差も有しているため、短時間で見るとオフセットがあるものと考えられる。この長周期の誤差は、GPS信号が伝播する際に通過する、大気圏や電離層における遅延モデルで生じる誤差によるところが大きい。この誤差は、大気圏や電離層の状態変化や衛星の配置によって変化するため、長周期の誤差となる。
【0035】
同一の衛星群から信号を受信し続けている場合は、時系列間において同一のオフセットが生じていると考えられるが、受信している複数の信号を送信している衛星群に変化があった場合には、オフセットの値に変化が生じる。
【0036】
受信している複数の信号を送信している衛星群が変化しないときに、その測位位置の変化から移動量を算出した場合は、実際の移動量とよく一致した値が得られる。しかし、受信信号の衛星群が変化した場合には、衛星群が変化した前後で、オフセット量やその正負が変化するため、移動量を算出すると誤差が大きくなる。
【0037】
図2に示すように、受信信号の衛星群に変化がない場合には、GPSの測位位置から算出した移動体の速度の時間変化と、移動体の実際の速度の時間変化とは、よく一致している。一方、受信信号の衛星群に変化があった場合には、衛星群が変化した前後で、GPSの測位位置から算出した移動体の速度の誤差が大きくなっている。
【0038】
上記図2に示すようなGPSの性質を考慮して、移動量推定部28では、以下に説明するように、GPS26の測位位置を用いて自車両の移動量を推定する。
【0039】
まず、GPS26によって測定された位置座標と、位置の測定に利用した複数の受信信号の各々に含まれる衛星の識別情報としての衛星番号を保持し、移動量を算出したい時刻間において、受信信号の衛星番号群を比較する。同一の衛星番号群である場合には、同一のオフセットが生じているものとして、移動量を算出したい時刻間の測位位置に基づいて、移動量を算出する。
【0040】
一方、移動量を算出したい時刻間において、衛星番号群に変化がある場合には、過去の測位結果を保持したデータの中から、同一の衛星番号群の信号を受信した時刻を探索する。該当時刻が存在する場合には、現時刻の測位位置と過去の該当時刻の測位位置とに基づいて、移動量を算出し、また、該当時刻から現在の一時刻前までの移動量が既知であれば、一時刻前から現時刻までの間の移動量を算出することができる。
【0041】
ただし、該当時刻と現時刻との間に、数10秒以上の時刻差がある場合、測位位置に大きな隔たりがある場合、又は、該当時刻がない場合については、受信信号の衛星番号群に差異があったとしても、測位位置の変化から移動量を算出すると共に、「低信頼」あるいは「推定誤差大」の可能性があるものとして扱う。
【0042】
位置推定部16は、計測された自車両の運動及び移動量に基づいて、自車両の相対位置の変化量を算出し、過去にGPS26で計測された絶対位置と相対位置の変化量とに基づいて、自車両の絶対位置を推定する。以下に、相対位置の変化量を算出する方法について図3を用いて説明する。
【0043】
まず、図3に示すような座標系において、初期位置を(X(0)、Y(0)、Z(0))とし、時刻tにおいて、移動量(速度V(t))、運動(3軸角速度(ヨー角速度φ(t)、ピッチ角速度θ(t)、ロール角速度ψ(t))、並進方向を示す成分(ax、ay、az))が得られると、下記の(1)式のように、時刻t+Δtにおける相対位置(X(t+Δt)、Y(t+Δt)、Z(t+Δt))を更新することができる。
【0044】
【数1】
ただし、ax/|a|、ay/|a|、及びaz/|a|は、正規化された並進方向を示す成分の各々を表わす。
【0045】
そして、自車両に搭載されたGPS26から、経度及び緯度などで表わされる絶対位置情報が過去に取得されていた場合には、絶対位置情報が示す絶対位置と、算出された相対位置とに基づいて、自車両の現在の絶対位置を計測する。なお、絶対位置情報を、必ずしもGPS26の測位位置から取得する必要はなく、例えば路側の位置が既知である設備から、通信などによって移動体の絶対位置情報を間欠的に取得してもよいし、予め位置が既知であるランドマーク情報が登録された地図などから、自車両の位置を取得してもよい。また、GPS26と同じ頻度で取得する必要もなく、移動量と運動とで補完できる間隔で取得すればよく、取得手段も問わない。
【0046】
次に、第1の実施の形態に係る位置計測装置10の作用について説明する。なお、位置計測装置10を搭載した自車両が走行しているときに、自車両の位置を計測する場合を例に説明する。
【0047】
位置計測装置10において、図4に示す位置計測処理ルーチンが実行される。まず、ステップ100において、GPS26によって自車両の位置を測位し、測位した位置を絶対位置の初期位置として設定する。
【0048】
そして、ステップ102において、運動推定処理を行い、計測単位期間Tms(例えば、第1時刻から第2時刻までの期間)の自車両の運動の3軸角速度の推定値、あるいは3軸角速度及び並進方向を示す成分の推定値を取得する。上記ステップ102の運動推定処理は、図5に示す運動推定処理ルーチンによって実現される。以下に運動推定処理ルーチンについて説明する。
【0049】
まず、ステップ120において、第1時刻及び第2時刻の各々のタイミングで画像撮像部20によって自車両前方を撮像した第1の画像及び第2の画像を取得する。そして、ステップ122において、上記ステップ120で取得した第1の画像を複数の領域に分割する。このステップ122では、例えば、予め与えられた定数y1、y2を各領域の境界線(y軸の切片座標)として、第1の画像を、画像上段の上方領域と、中段の中間領域と、下段の下方領域との3つの領域に分割する。
【0050】
次のステップ124において、上記ステップ122で得られた各分割画像中から、それぞれ所定の数だけ特徴点(不動点であることが期待される点)を抽出する。ここで、特徴点とは、周囲の点とは区別でき、追跡が容易な点のことである。このステップ124では、例えば、二次元的に濃淡変化の勾配値が大きくなる画素を検出する、一般的な周知の方法(例えば、Harrisオペレータなど)を用いて、自動的に特徴点を抽出する。
【0051】
そして、ステップ126において、上記ステップ124で抽出された第1の画像における特徴点の各々を、上記ステップ120で取得した第2の画像において追跡し、第2の画像から、対応する各特徴点をそれぞれ検出する。このとき、第1時刻と第2時刻との間で対応する特徴点は、同一の不動点であることが期待されているので、第1の画像と第2の画像とで対応する点とその周囲の点との輝度は、殆ど変化しないという仮定を用い、この仮定に基づいて各特徴点の追跡を行う。この様な追跡処理を実施する際には、例えば、周知のルーカスカナデ法(Lucas−Kanade法)などを用いると良い。また、各特徴点の追跡処理の後、第1の画像と第2の画像との間で対応する特徴点の各組の各画像内座標を、1つの適当な形式に整理する。
【0052】
そして、ステップ128において、上記ステップ126で得られた第1時刻と第2時刻との間で互いに対応する特徴点の組から、第1時刻と第2時刻との間における画像撮像部20の相対的な位置関係、即ち、画像撮像部20を搭載した自車両の第1時刻と第2時刻との間における運動に関連する情報を推定する。
【0053】
以下、図6を用いて、抽出された特徴点(不動点)と移動体の運動との関係を説明する。図6に示すように、移動体の運動は、第1の画像から第2の画像への回転行列Rと並進ベクトルtとから構成される。例えば、第1の画像上の特徴点x1に対応する第2の画像上の特徴点をx2とすると、これらの特徴点は、同一の不動点を示す点である。
【0054】
ここで、カメラの撮像特性に依る画像の歪みを補正するキャリブレーション行列K、及び基礎行列Fを用いて、対応する2つの特徴点の関係を、以下の(2)式〜(4)式で表すことができる。
【0055】
【数2】
【0056】
上記の(2)式〜(4)式のように、上記の2時刻間における8組以上の互いに対応する特徴点の組(即ち、同一の不動点)が求まっていれば、それに基づいて、上記の基礎行列Fを算出することができる。そして、カメラのキャリブレーション行列Kが既知であれば、上記の式からも分かる様に、その2時刻間における移動体の運動、即ち、その2時刻間における自車両の回転行列Rと並進ベクトルtとを求めることができる。
【0057】
そして、ステップ130において、上記ステップ128で求めた自車両の回転行列Rと並進ベクトルtとに基づいて、自車両の3軸角速度、あるいは3軸角速度及び並進方向を示す成分を算出し、算出した3軸角速度、あるいは3軸角速度及び並進方向を示す成分を、自車両の運動の推定値として出力して、運動推定処理ルーチンを終了する。
【0058】
上述したように、運動推定処理によって、自車両の運動の3軸角速度、あるいは3軸角速度及び並進方向を示す成分が取得される。
【0059】
次に、位置計測処理ルーチンのステップ104で、移動量推定処理を行い、第1時刻と第2時刻との間の自車両の移動量としての車速の推定値を取得する。上記の移動量推定処理は、図7に示す移動量推定処理ルーチンによって実現される。以下に、移動量推定処理ルーチンについて説明する。
【0060】
まず、ステップ140において、GPS26から、第1時刻及び第2時刻の各々で測定された測位位置と、第1時刻及び第2時刻の各々において測位に利用した複数の受信信号に含まれる複数の衛星番号からなる衛星番号群とを取得する。そして、ステップ142で、上記ステップ140で取得した測位位置及び衛星番号群をメモリ(図示省略)に記憶する。
【0061】
次のステップ144では、上記ステップ140で取得した第1時刻における衛星番号群と、第2時刻における衛星番号群とが同一であるか否かを判定する。上記ステップ144で、第1時刻と第2時刻とにおいて受信信号の衛星番号群が同一であると判定された場合には、第1時刻と第2時刻とにおいて同一の衛星群から信号を受信したと判断し、ステップ146において、第1時刻及び第2時刻で測定された測位位置の変化に基づいて、自車両の移動量としての車速を算出して、移動量推定処理ルーチンを終了する。
【0062】
一方、上記ステップ144で、第1時刻と第2時刻とにおいて受信信号の衛星番号群が同一でないと判定された場合には、第1時刻と第2時刻とにおいて異なる衛星群から信号を受信したと判断し、ステップ148において、メモリに記憶された過去の衛星番号群から、第2時刻における衛星番号群と同一の衛星番号群である時刻を検索する。
【0063】
そして、ステップ150において、上記ステップ148で第2時刻と同一の衛星番号群であった時刻が検索されたか否かを判定する。上記ステップ150で、第2時刻と同一の衛星番号群であった時刻が検索されたと判定された場合には、ステップ152において、メモリに記憶された、検索された時刻に測定された測位位置と、第2時刻における測位位置とに基づいて、自車両の車速を算出し、移動量推定処理ルーチンを終了する。
【0064】
一方、上記ステップ150で、第2時刻と同一の衛星番号群であった時刻が検索されなかったと判定された場合には、ステップ154において、第1時刻及び第2時刻の各々で測定された測位位置に基づいて、信頼度の低い推定値として、自車両の車速を算出して、移動量推定処理ルーチンを終了する。
【0065】
上述したように、移動量推定処理によって、自車両の移動量としての車速が取得される。
【0066】
次に、位置計測処理ルーチンのステップ106において、上記ステップ102で得られた自車両の運動の3軸角速度、あるいは3軸角速度及び並進方向を示す成分と、上記ステップ104で得られた自車両の車速とに基づいて、自車両の相対位置の変化量を算出する。次のステップ108では、上記ステップ100で設定された初期位置と、上記ステップ106で算出された相対位置の変化量とに基づいて、自車両の絶対位置を算出する。
【0067】
そして、ステップ112において、上記ステップ102からステップ110によって絶対位置を計測してから、計測単位期間Tmsが経過したか否かを判定し、計測単位期間Tmsが経過すると、ステップ102へ戻り、再び自車両の絶対位置を計測して、自車両の絶対位置を更新する。
【0068】
以上説明したように、第1の実施の形態に係る位置計測装置によれば、同一の衛星群からの信号に基づいて測位された自車両の複数の位置に基づいて、自車両の移動量を算出するため、簡易な構成で、ロバスト性を確保して、自車両の移動量を精度よく計測することができる。
【0069】
また、測位に利用した衛星群に変化がない場合の測位位置の変化を用いて、移動量を計測するため、移動量の計測精度を確保することができる。
【0070】
また、ロバスト性を確保して精度よく計測された自車両の車速と、推定された自車両の運動の3軸角速度及び並進方向を示す成分とに基づいて、自車両の位置を計測することにより、簡易な構成で、ロバスト性を確保して、自車両の位置を精度よく計測することができる。
【0071】
次に、第2の実施の形態について説明する。なお、第1の実施の形態と同様の構成となっている部分については、同一符号を付して説明を省略する。
【0072】
第2の実施の形態では、GPSによる測位位置と、車速センサからの車速パルス信号のパルス数とを組み合わせて、自車両の移動量を推定している点が第1の実施の形態と異なっている。
【0073】
図8に示すように、第2の実施の形態に係る位置計測装置210の移動量計測部214は、GPS26と、車速に応じて車速パルス信号を出力する車速パルスセンサで構成され、自車両の対地速に応じた車速パルス信号を出力する車速センサ226と、GPS26によって測定された自車両の位置、及び車速センサ226の車速パルス信号に基づいて、移動量を推定する移動量推定部228とを備えている。
【0074】
車速センサ226は、例えば、車速が60km/h(回転数が637rpm)のときに、1秒間におけるパルス数が10.6となる車速パルス信号を出力する。
【0075】
移動量推定部228は、計測単位期間Tmsにおいて、GPS26で位置を測定した時に利用した受信信号に含まれる衛星番号群に変化が無い場合、測位位置の変化から車速を算出して移動量の推定値とし、衛星番号群が変化した場合、車速センサ226からの車速パルス信号のパルス数をカウントし、パルス数から車速を算出して移動量の推定値とする。
【0076】
次に、第2の実施の形態に係る移動量推定処理ルーチンについて、図9を用いて説明する。なお、第1の実施の形態と同様の処理については、同一符号を付して説明を省略する。
【0077】
まず、ステップ140において、GPS26から、第1時刻及び第2時刻の各々で測定された測位位置と、第1時刻及び第2時刻の各々において受信した信号に含まれる衛星番号群とを取得する。
【0078】
次のステップ144では、上記ステップ140で取得した第1時刻における衛星番号群と、第2時刻における衛星番号群とが同一であるか否かを判定する。上記ステップ144で、第1時刻と第2時刻とにおける衛星番号群が同一であると判定された場合には、ステップ146において、第1時刻及び第2時刻で測定された測位位置の変化に基づいて、自車両の移動量としての車速を算出して、移動量推定処理ルーチンを終了する。
【0079】
一方、上記ステップ144で、第1時刻と第2時刻とにおける衛星番号群が同一でないと判定された場合には、ステップ250において、第1時刻から第2時刻までの間に車速センサ226から出力された車速パルス信号のパルス数を取得する。そして、ステップ252で、上記ステップ250で取得されたパルス数に基づいて、自車両の車速を算出して、移動量推定処理ルーチンを終了する。
【0080】
なお、位置計測処理ルーチン及び運動推定処理ルーチンについては、第1の実施の形態と同様であるため、説明を省略する。
【0081】
以上説明したように、第2の実施の形態に係る位置計測装置によれば、同一の衛星群からの信号に基づいて測位された自車両の複数の位置に基づいて、自車両の移動量を算出し、異なる衛星群から信号を受信した場合には、速度センサからの車速パルス信号のパルス数に基づいて移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0082】
また、速度センサの車速パルス信号を用いることで、短期的な移動量の計測精度を確保することができる。
【0083】
次に、第3の実施の形態について説明する。なお、第1の実施の形態と同様の構成となっている部分については、同一符号を付して説明を省略する。
【0084】
第3の実施の形態では、GPSによる測位位置と、撮像画像とを組み合わせて、自車両の移動量を推定している点が第1の実施の形態と異なっている。
【0085】
図10に示すように、第3の実施の形態に係る位置計測装置310は、自車両の運動を計測すると共に、自車両の移動量を計測する運動移動量計測部312と、位置推定部16とを備えている。
【0086】
運動移動量計測部312は、画像撮像部20と、運動推定部24と、GPS26と、GPS26によって測定された自車両の位置、及び画像撮像部20によって撮像された撮像画像に基づいて、移動量を推定する移動量推定部328とを備えている。
【0087】
移動量推定部328は、計測単位期間Tmsにおいて、GPS26で位置を測定した時に利用した受信信号に含まれる衛星番号群に変化が無い場合、測位位置の変化から車速を算出して移動量の推定値とし、衛星番号群が変化した場合、画像撮像部20によって撮像された複数の撮像画像に基づいて、車速を算出して移動量の推定値とする。
【0088】
移動量推定部328は、撮像された複数の撮像画像に基づいて、車速を算出する場合、撮像画像から復元した3次元空間の路面平面を推定する。そして、予め計測された、画像撮像部20の単眼カメラが設置されている路面上からの高さや単眼カメラの光軸の角度に基づいて、第1時刻及び第2時刻に撮像された撮像画像の各々から、路面上の近傍の特徴点までの距離を計測し、特徴点までの距離の変化に基づいて、車速を算出する。特に、低速時は移動量が少ないため、自車両周辺の路面を路面平面として推定してもよい。
【0089】
次に、第3の実施の形態に係る移動量推定処理ルーチンについて、図11を用いて説明する。なお、第1の実施の形態と同様の処理については、同一符号を付して、詳細な説明を省略する。
【0090】
まず、ステップ140において、GPS26から、第1時刻及び第2時刻の各々で測定された測位位置と、第1時刻及び第2時刻の各々における衛星番号群とを取得する。
【0091】
次のステップ144では、上記ステップ140で取得した第1時刻における衛星番号群と、第2時刻における衛星番号群とが同一であるか否かを判定する。上記ステップ144で、第1時刻と第2時刻とにおける衛星番号群が同一であると判定された場合には、ステップ146において、第1時刻及び第2時刻で測定された測位位置の変化に基づいて、自車両の移動量としての車速を算出して、移動量推定処理ルーチンを終了する。
【0092】
一方、上記ステップ144で、第1時刻と第2時刻とにおける衛星番号群が同一でないと判定された場合には、ステップ350において、第1時刻より前に撮像された撮像画像を画像撮像部20から取得して、取得した撮像画像に基づいて、3次元空間の路面平面を推定する。
【0093】
そして、ステップ352で、画像撮像部20から第1時刻及び第2時刻の各々で撮像された撮像画像を取得し、ステップ354において、第1時刻に撮像された撮像画像から、上記ステップ350で推定された路面平面上の特徴点までの距離を算出すると共に、第2時刻に撮像された撮像画像から、上記と同じ特徴点までの距離を算出する。
【0094】
次のステップ356では、上記ステップ354で算出された第1時刻に対する距離と第2時刻に対する距離との変化に基づいて、自車両の車速を算出して、移動量推定処理ルーチンを終了する。
【0095】
なお、他の処理については、第1の実施の形態と同様であるため、説明を省略する。
【0096】
このように、同一の衛星群からの信号に基づいて測位された自車両の複数の位置に基づいて、自車両の移動量を算出し、異なる衛星群から信号を受信した場合には、単眼のカメラによって撮像された撮像画像に基づいて移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0097】
なお、上記の実施の形態では、撮像画像から推定された路面平面上の特徴点と、画像撮像部の高さ及び角度とに基づいて、特徴点までの距離を算出する場合を例に説明したが、これに限定されるものではなく、例えば、撮像画像から、大きさが既知の物体(例えば、白線や標識)を抽出し、抽出された物体の撮像画像上の大きさに基づいて、抽出された物体までの距離を算出するようにしてもよい。この場合には、算出された大きさが既知の物体までの距離の変化から、車速を算出するようにすればよい。
【0098】
次に、第4の実施の形態について説明する。なお、第1の実施の形態及び第2の実施の形態と同様の構成となっている部分については、同一符号を付して説明を省略する。
【0099】
第4の実施の形態では、速度センサからの車速パルス信号と、撮像画像とを組み合わせて、自車両の移動量を推定している点が第1の実施の形態と異なっている。
【0100】
図12に示すように、第4の実施の形態に係る位置計測装置410は、自車両の運動を計測すると共に、自車両の移動量を計測する運動移動量計測部412と、位置推定部16と、GPS26とを備えている。
【0101】
運動移動量計測部412は、画像撮像部20と、運動推定部24と、車速センサ226と、車速センサ226から出力された車速パルス信号、及び画像撮像部20によって撮像された撮像画像に基づいて、自車両の移動量を推定する移動量推定部428とを備えている。
【0102】
ここで、本発明の原理について説明する。車速センサ226から出力される車速パルス信号は、連続値ではなく、一定の移動が行われた場合に出力される仕組みである。車速パルス信号から算出される車速の時系列変化と、実際の車速の時系列変化とを比較すると、車速パルス信号から算出される車速の値は、真値に対して大きなズレはないが、連続値ではなく離散値となる。信号処理を工夫して、推定精度を上げることは可能であるが、一般的な車速パルス信号の移動量分解能は40cm程度であるため、低速時に、移動量が得られないことがある。
【0103】
本実施の形態では、上記の車速パルス信号を出力する車速センサ226の性質を考慮し、移動量推定部428は、車速センサ226から出力される車速パルス信号を用いて、以下に説明するように、自車両の移動量としての車速を推定する。
【0104】
まず、車速センサ226から計測単位期間Tmsに出力された車速パルス信号のパルス数をカウントする。このパルス数が、たとえば1以上である場合は、パルス数から換算される車速を移動量として出力する。計測単位期間Tms内にパルス数がカウントされなかった場合(パルス数が0の場合)には、移動の有無の計測が容易ではない。この場合、撮像画像に基づいて計測される計測単位期間Tmsの距離の変化から、移動量としての車速を算出する。なお、パルス数に関する閾値を暫定的に「1」としたが、特に1に限定されるものではなく、他の適切な値を閾値として採用してもよい。
【0105】
次に、第4の実施の形態における移動量推定処理ルーチンについて、図13を用いて説明する。なお、第2の実施の形態及び第3の実施の形態と同様の処理については、同一符号を付して詳細な説明を省略する。
【0106】
まず、ステップ250において、第1時刻から第2時刻までの間に車速センサ226から出力された車速パルス信号のパルス数Npを取得する。そして、ステップ450において、上記ステップ250で取得されたパルス数Npが、1以上であるか否かを判定する。上記ステップ450で、パルス数Npが1以上であると判定された場合には、ステップ252で、上記ステップ250で取得されたパルス数Npに基づいて、自車両の車速を算出して、移動量推定処理ルーチンを終了する。
【0107】
一方、上記ステップ450で、パルス数Npが0であると判定された場合には、ステップ350において、第1時刻より前に撮像された画像に基づいて、3次元空間の路面平面を推定する。そして、ステップ352で、画像撮像部20から第1時刻及び第2時刻の各々で撮像された撮像画像を取得し、ステップ354において、第1時刻及び第2時刻の各々に撮像された撮像画像から、上記ステップ350で推定された路面平面上の特徴点までの距離を各々算出する。
【0108】
次のステップ356では、上記ステップ354で算出された第1時刻に対する距離と第2時刻に対する距離との変化に基づいて、自車両の車速を算出して、移動量推定処理ルーチンを終了する。
【0109】
なお、他の処理については、第1の実施の形態と同様であるため、説明を省略する。
【0110】
以上説明したように、第4の実施の形態に係る位置計測装置によれば、計測単位期間における車速パルス信号のパルス数が1以上である場合には、車速パルス信号のパルス数に基づいて移動量を算出し、計測単位期間における車速パルス信号のパルス数が0である場合には、単眼のカメラによって撮像された撮像画像に基づいて移動量を算出するため、簡易な構成で、ロバスト性を確保して、移動体の移動量を精度よく計測することができる。
【0111】
また、自車両が低速時には、低速時の計測精度が高くない車速パルス信号の代わりに、撮像画像から推定した近傍の路面平面上の特徴点までの距離の変化に基づいて、移動量を算出することにより、停止判定も可能になる。
【0112】
次に、第5の実施の形態について説明する。なお、第1の実施の形態及び第4の実施の形態と同様の構成となっている部分については、同一符号を付して説明を省略する。
【0113】
第5の実施の形態では、GPSによる測位位置と、撮像画像と、車速センサからの車速パルス信号とを組み合わせて、自車両の移動量を推定している点が、第4の実施の形態と異なっている。
【0114】
図14に示すように、第5の実施の形態に係る位置計測装置510の運動移動量計測部512は、画像撮像部20と、運動推定部24と、GPS26と、車速センサ226と、GPS26によって測定された測位位置、車速センサ226から出力された車速パルス信号、及び画像撮像部20によって撮像された撮像画像に基づいて、移動量を推定する移動量推定部528とを備えている。
【0115】
移動量推定部528は、計測単位期間Tmsにおいて、GPS26で位置を測定した時に利用した受信信号に含まれる衛星番号群に変化が無い場合、測位位置の変化から車速を算出して移動量の推定値とする。一方、衛星番号群が変化した場合、車速センサ226から計測単位期間Tmsに出力された車速パルス信号のパルス数をカウントし、このパルス数が、たとえば1以上である場合は、車速センサ226の車速パルス信号から車速を算出して移動量の推定値とする。計測単位期間Tms内にパルス数がカウントされなかった場合には、撮像画像に基づいて計測される計測単位期間Tmsの距離の変化から、移動量としての車速を算出する。
【0116】
次に、第5の実施の形態に係る移動量推定処理ルーチンについて、図15を用いて説明する。なお、第1の実施の形態〜第4の実施の形態と同様の処理については、同一符号を付して詳細な説明を省略する。
【0117】
まず、ステップ140において、GPS26から、第1時刻及び第2時刻の各々で測定された測位位置と、第1時刻及び第2時刻の各々における衛星番号群とを取得する。
【0118】
次のステップ144では、上記ステップ140で取得した第1時刻における衛星番号群と、第2時刻における衛星番号群とが同一であるか否かを判定する。上記ステップ144で、第1時刻と第2時刻とにおける衛星番号群が同一であると判定された場合には、ステップ146において、第1時刻及び第2時刻で測定された測位位置の変化に基づいて、自車両の移動量としての車速を算出して、移動量推定処理ルーチンを終了する。
【0119】
一方、上記ステップ144で、第1時刻と第2時刻とにおける衛星番号群が同一でないと判定された場合には、ステップ250において、第1時刻から第2時刻までの間に車速センサ226から出力された車速パルス信号のパルス数を取得する。そして、ステップ450において、上記ステップ250で取得されたパルス数Npが、1以上であるか否かを判定する。上記ステップ450で、パルス数Npが1以上であると判定された場合には、ステップ252で、上記ステップ250で取得されたパルス数Npに基づいて、自車両の車速を算出して、移動量推定処理ルーチンを終了する。
【0120】
一方、上記ステップ450で、パルス数Npが0であると判定された場合には、ステップ350において、第1時刻より前に撮像された撮像画像に基づいて、3次元空間の路面平面を推定する。そして、ステップ352で、画像撮像部20から第1時刻及び第2時刻の各々で撮像された撮像画像を取得し、ステップ354において、第1時刻及び第2時刻の各々に撮像された撮像画像から、上記ステップ350で推定された路面平面上の特徴点までの距離を各々算出する。
【0121】
次のステップ356では、上記ステップ354で算出された第1時刻に対する距離と第2時刻に対する距離との変化に基づいて、自車両の車速を算出して、移動量推定処理ルーチンを終了する。
【0122】
このように、計測単位期間の第1時刻と第2時刻とにおいて、異なる衛星群から信号を受信した場合であって、かつ、車速パルス信号のパルス数が大きい場合には、速度センサからの車速パルス信号に基づいて移動量を算出し、異なる衛星群から信号を受信した場合であって、かつ、車速パルス信号のパルス数が小さい場合には、単眼のカメラによって撮像された撮像画像に基づいて移動量を算出するため、簡易な構成で、ロバスト性を確保して、自車両の移動量を精度よく計測することができる。
【0123】
次に、第6の実施の形態について説明する。なお、第6の実施の形態に係る位置計測装置は、第5の実施の形態と同様の構成となっているため、同一符号を付して説明を省略する。
【0124】
第6の実施の形態では、カルマンフィルタ等を用いた最適処理手法によって、自車両の移動量としての車速を推定している点が第5の実施の形態と異なっている。
【0125】
第6の実施の形態に係る位置計測装置510の移動量推定部528では、画像撮像部20の撮像画像に基づいて計測される距離と、GPS26によって測定される測位位置と、車速センサ226からの車速パルス信号のパルス数とを観測値とし、自車両の車速を推定値として、カルマンフィルタによる最適処理手法を用いた処理を行う。
【0126】
カルマンフィルタによる最適処理手法では、画像撮像部20の単眼カメラの撮像画像に基づいて計測される距離の観測値yimageが得られた場合に、観測値yimageのバラつきを考慮して、観測ノイズwimageを設定し、観測行列Himageについて、観測値yimageと推定値xとの関係式から、観測行列Himageを設定する。
【0127】
また、GPS26によって測位された測位位置の観測値ygpsが得られた場合に、観測値ygpsのバラつきを考慮して、観測ノイズwgpsを設定し、観測行列Hgpsについて、観測値ygpsと推定値xとの関係式から、観測行列Hgpsを設定する。
【0128】
また、車速センサ226から出力された車速パルス信号のパルス数の観測値ypulsが得られた場合に、観測値ypulsのバラつきを考慮して、観測ノイズwpulsを設定し、観測行列Hpulsについて、観測値ypulsと推定値xとの関係式から、観測行列Hpulsを設定する。
【0129】
そして、設定された観測行列、観測値、及び観測ノイズに基づいて、カルマンフィルタで、推定値x及び推定誤差共分散行列を算出して出力する。観測値が得られない場合には、推定値の予測値と推定誤差共分散行列の予測値とを算出して出力する。
【0130】
次に、第6の実施の形態に係る位置計測装置510の作用について説明する。なお、第1の実施の形態と同様の部分については、同一符号を付して詳細な説明を省略する。
【0131】
第6の実施の形態に係る位置計測装置510では、上述した移動量推定処理ルーチンにおけるステップ350からステップ354までの処理と同様の処理が繰り返し実行され、算出された路面平面上の特徴点までの距離が繰り返し出力される。また、GPS26によって自車両の位置が繰り返し測定される。また、車速センサ226から、車速パルス信号が連続して出力され、パルス数が繰り返しカウントされる。
【0132】
また、位置計測装置510において、図16に示す移動量推定処理ルーチンが実行される。まず、ステップ600において、撮像画像に基づいて計測される路面平面上の特徴点までの距離を取得し、ステップ602において、GPS26から測定された測位位置を取得する。そして、ステップ604において、車速センサ226からの車速パルス信号に基づいてカウントされたパルス数を取得する。
【0133】
そして、ステップ606において、上記ステップ600で取得した距離と、上記ステップ602で取得した測位位置と、上記ステップ604で取得したパルス数とに基づいて、移動量としての車速の推定値xに対して、初期値を設定する。次のステップ608では、時間ti−1にΔtを加算して、時間tiに更新する。そして、ステップ610において、1時刻前の推定値xに基づいて、推定値xの予測値を算出すると共に、推定誤差共分散の予測値を算出する。
【0134】
次のステップ612において、現タイミングに、撮像画像に基づいて計測される路面平面上の特徴点までの距離、GPS26から測定された測位位置、及び車速センサ226からの車速パルス信号のパルス数の何れかが、観測値として入力されたか否かを判定し、現タイミングに観測値が入力されなかった場合には、ステップ614において、上記ステップ610で算出された推定値xの予測値の車速を、自車両の移動量の推定値として出力し、ステップ608へ戻る。
【0135】
上記ステップ612において、現タイミングに観測値が入力されたと判定された場合には、ステップ616において、どの観測値が入力されたのかを判定する。
【0136】
上記ステップ616で、観測値としてGPS26から測定された測位位置が入力されたと判定された場合には、ステップ618において、入力された測位位置を観測値として、観測行列Hgps、観測値ygps、及び観測ノイズwgpsを設定する。
【0137】
一方、上記ステップ616で、観測値として、撮像画像に基づいて計測される路面平面上の特徴点までの距離が入力されたと判定された場合には、ステップ620において、入力された特徴点までの距離を観測値として、観測行列Himage、観測値yimage、及び観測ノイズwimageを設定する。
【0138】
また、上記ステップ616で、観測値として、車速センサ226からの車速パルス信号のパルス数が入力されたと判定された場合には、ステップ622において、入力されたパルス数を観測値として、観測行列Hpuls、観測値ypuls、及び観測ノイズwpulsを設定する。
【0139】
そして、ステップ624において、上記ステップ618、上記ステップ620、又は上記ステップ622で設定された観測行列、観測値、及び観測ノイズに基づいて、カルマンフィルタで、カルマンゲインK、推定値x、及び推定誤差共分散行列を算出する。次のステップ626では、上記ステップ624で算出された推定値xの車速を、自車両の移動量の推定値として出力し、ステップ608へ戻る。
【0140】
このように、カルマンフィルタを用いた最適処理手法によって、自車両の外部を撮像した複数の撮像画像に基づいて計測された路面平面上の特徴点までの距離と、GPSによって測定された測位位置と、車速センサからの車速パルス信号のパルス数とに基づいて、自車両の移動量の推定値を算出することにより、簡易な構成で、ロバスト性を確保して、自車両の移動量を精度よく計測することができる。
【0141】
また、ロバスト性を確保して精度よく計測された自車両の移動量と、推定された自車両の運動の3軸角速度及び並進方向を示す成分とに基づいて、自車両の位置を計測することにより、簡易な構成で、ロバスト性を確保して、自車両の位置を精度よく計測することができる。
【0142】
なお、上記の実施の形態では、カルマンフィルタで推定する例を記述したが、周知のパーティクルフィルタなど、他の最適推定手法を用いてもよく、適切に解を推定できる手法であれば、その方法はこの限りではない。
【0143】
また、誤差共分散を、推定された自車両の移動量の推定値に対する信頼度として用いてもよい。例えば、推定値と共に算出された誤差共分散によって、推定値の信頼度が低く、推定値の誤差が大きいと判断される場合に、この推定値を、自車両の移動量の推定値として採用しないようにしてもよい。
【0144】
次に、第7の実施の形態について説明する。なお、第7の実施の形態に係る位置計測装置は、第1の実施の形態と同様の構成となっているため、同一符号を付して説明を省略する。
【0145】
第7の実施の形態では、カルマンフィルタを用いた最適処理手法によって、自車両の位置を計測している点が第1の実施の形態と異なっている。
【0146】
第7の実施の形態に係る位置計測装置10の位置推定部16では、運動推定部24によって推定された運動(3軸角速度、並進方向を示す成分)、移動量推定部28によって推定された移動量(車速)、及びGPS26によって測定された絶対位置(Xg,Yg,Zg)を観測値とし、自車の絶対位置、方位角、車速、及び3軸角速度を推定値として、カルマンフィルタを用いた最適処理手法による位置計測処理を行う。
【0147】
カルマンフィルタを用いた最適処理手法では、運動推定部24によって推定された運動の観測値ymoveが得られた場合に、観測値ymoveのバラつきを考慮して、観測ノイズwmoveを設定し、観測行列Hmoveについて、観測値ymoveと推定値xとの関係式から、観測行列Hmoveを設定する。
【0148】
また、移動量推定部28によって推定された移動量の観測値ytravが得られた場合に、観測値ytravのバラつきを考慮して、観測ノイズwtravを設定し、観測行列Htravについて、観測値ytravと推定値xとの関係式から、観測行列Htravを設定する。
【0149】
また、GPS26によって測定された絶対位置の観測値ygpsが得られた場合に、観測値ygpsのバラつきを考慮して、観測ノイズwgpsを設定し、観測行列Hgpsについて、観測値ygpsと推定値xとの関係式から、観測行列Hgpsを設定する。
【0150】
そして、設定された観測行列、観測値、及び観測ノイズに基づいて、カルマンフィルタで、推定値x及び推定誤差共分散行列を算出して出力する。観測値が得られない場合には、推定値の予測値と推定誤差共分散行列の予測値とを算出して出力する。
【0151】
次に、第7の実施の形態に係る位置計測装置10の作用について説明する。なお、第1の実施の形態と同様の部分については、同一符号を付して詳細な説明を省略する。
【0152】
第7の実施の形態に係る位置計測装置10では、第1の実施の形態で説明した運動推定処理ルーチンが繰り返し実行され、また、上述した移動量推定処理ルーチンが繰り返し実行される。また、GPS26によって、所定時間ごと(例えば1秒ごと)に自車両の絶対位置が測定される。
【0153】
また、位置計測装置10において、図17に示す位置計測処理ルーチンが実行される。まず、ステップ700において、運動推定処理ルーチンによって推定された3軸角速度及び並進方向を示す成分を取得し、ステップ702において、移動量推定処理ルーチンによって推定された車速を取得する。次のステップ704では、GPS26によって測定された自車両の絶対位置を取得する。
【0154】
そして、ステップ706において、上記ステップ700で取得した3軸角速度及び並進方向を示す成分と、上記ステップ702で取得した車速と、上記ステップ704で取得した絶対位置とに基づいて、絶対位置、方位角、車速、及び3軸角速度の推定値xに対して、初期値を設定する。次のステップ708では、時間ti−1にΔtを加算して、時間tiに更新する。そして、ステップ710において、1時刻前の推定値xに基づいて、推定値xの予測値を算出すると共に、推定誤差共分散の予測値を算出する。
【0155】
次のステップ712において、現タイミングに、観測値として、運動推定処理ルーチンによって推定された3軸角速度、移動推定処理ルーチンによって推定された車速、又はGPS26によって測定された絶対位置が入力されたか否かを判定し、現タイミングに観測値が入力されなかった場合には、ステップ714において、上記ステップ710で算出された推定値xの予測値の絶対位置を、自車両の絶対位置の推定値として出力し、ステップ708へ戻る。
【0156】
上記ステップ712において、現タイミングに観測値が入力されたと判定された場合には、ステップ716において、上記ステップ712でどの観測値が入力されたかを判定する。
【0157】
上記ステップ716で、観測値として、運動推定処理ルーチンで推定された3軸角速度及び並進方向を示す成分が入力されたと判定された場合には、ステップ718において、入力された3軸角速度及び並進方向を示す成分を観測値として、観測行列Hmove、観測値ymove、及び観測ノイズwmoveを設定する。
【0158】
一方、上記ステップ716で、観測値として、移動量推定処理ルーチンによって推定された車速が入力されたと判定された場合には、ステップ720において、入力された車速を観測値として、観測行列Htrav、観測値ytrav、及び観測ノイズwtravを設定する。
【0159】
また、上記ステップ716で、観測値として、GPS26によって測定された絶対位置が入力されたと判定された場合には、ステップ722において、入力された絶対位置を観測値として、観測行列Hgps、観測値ygps、及び観測ノイズwgpsを設定する。
【0160】
そして、ステップ724において、上記ステップ718、上記ステップ720、又は上記ステップ722で設定された観測行列、観測値、及び観測ノイズに基づいて、カルマンフィルタで、カルマンゲインK、推定値x、及び推定誤差共分散行列を算出する。次のステップ726では、上記ステップ724で算出された推定値xの絶対位置を、自車両の絶対位置の推定値として出力し、ステップ708へ戻る。
【0161】
このように、カルマンフィルタを用いた最適処理手法によって、推定された自車両の運動の3軸角速度及び並進方向を示す成分と、推定された自車両の車速と、GPSによって測位された位置とに基づいて、自車両の絶対位置の推定値を算出することにより、簡易な構成で、ロバスト性を確保して、自車両の絶対位置を精度よく計測することができる。
【0162】
なお、上記の実施の形態では、カルマンフィルタを用いた位置計測処理において、推定された自車両の運動の3軸角速度及び並進方向を示す成分、推定された自車両の車速、及びGPSによって測位された絶対位置が、観測値として入力される場合を例に説明したが、これに限定されるものではない。例えば、推定された自車両の運動の3軸角速度及び並進方向を示す成分と、及び推定された自車両の車速とに基づいて推定された自車両の相対位置の変化量が、観測値として、カルマンフィルタを用いた位置計測処理に入力されるようになっていてもよい。この場合には、推定された自車両の相対位置の変化量、及びGPSによって測位された絶対位置を観測値とし、自車の絶対位置、方位角、車速、及び3軸角速度を推定値として、カルマンフィルタを用いた最適処理手法による位置計測処理を行うようにすればよい。
【0163】
また、カルマンフィルタに限らず、パーティクルフィルタなどの最適推定手法を用いてもよい。
【0164】
次に、第8の実施の形態について説明する。第8の実施の形態では、ジャイロセンサによって、自車両の運動の3軸角速度を検出している。
【0165】
第8の実施の形態に係る位置計測装置は、ジャイロセンサによって検出された自車両の運動の3軸角速度と、移動量推定部によって推定された自車両の移動量とに基づいて、自車両の位置を推定する。
【0166】
次に、第9の実施の形態について説明する。第9の実施の形態では、ジャイロセンサ及び単眼カメラによって撮像された撮像画像に基づいて、自車両の運動の3軸角速度、あるいは3軸角速度及び並進方向を示す成分を推定している。
【0167】
第9の実施の形態に係る位置計測装置は、ジャイロセンサ及び撮像画像によって推定された自車両の運動の3軸角速度、あるいは3軸角速度及び並進方向を示す成分と、移動量推定部によって推定された自車両の移動量とに基づいて、自車両の位置を推定する。
【0168】
なお、上記の第1の実施の形態〜第9の実施の形態では、撮像画像又はジャイロセンサに基づいて3軸角速度を推定する場合を例に説明したが、これに限定されるものではなく、撮像画像又はジャイロセンサに基づいて、ヨー角速度、ピッチ角速度、及びロール角速度の少なくとも1つを推定するように構成してもよい。
【0169】
また、自車両の移動量として、車速を推定する場合を例に説明したが、これに限定されるものではなく、移動量として、移動距離又は移動ベクトルを推定するようにしてもよい。
【0170】
また、位置計測処理をオンラインで逐次処理する場合を例に説明したが、これに限定されるものではなく、自車両の運動の推定値及び移動量の推定値を取得した結果を用いて、オフラインで後処理として位置計測処理を行うようにしてもよい。
【図面の簡単な説明】
【0171】
【図1】本発明の第1の実施の形態に係る位置計測装置を示すブロック図である。
【図2】GPSの測位位置に基づいて算出される車速の時間変化と実際の車速の時間変化とを示すグラフである。
【図3】相対位置の変化量の算出方法を説明するための座標系を示すイメージ図である。
【図4】本発明の第1の実施の形態に係る位置計測装置における位置計測処理ルーチンの内容を示すフローチャートである。
【図5】本発明の第1の実施の形態に係る位置計測装置における運動推定処理ルーチンの内容を示すフローチャートである。
【図6】抽出された特徴点と移動体の運動との関係を説明するための図である。
【図7】本発明の第1の実施の形態に係る位置計測装置における移動量推定処理ルーチンの内容を示すフローチャートである。
【図8】本発明の第2の実施の形態に係る位置計測装置を示すブロック図である。
【図9】本発明の第2の実施の形態に係る位置計測装置における移動量推定処理ルーチンの内容を示すフローチャートである。
【図10】本発明の第3の実施の形態に係る位置計測装置を示すブロック図である。
【図11】本発明の第3の実施の形態に係る位置計測装置における移動量推定処理ルーチンの内容を示すフローチャートである。
【図12】本発明の第4の実施の形態に係る位置計測装置を示すブロック図である。
【図13】本発明の第4の実施の形態に係る位置計測装置における移動量推定処理ルーチンの内容を示すフローチャートである。
【図14】本発明の第5の実施の形態に係る位置計測装置を示すブロック図である。
【図15】本発明の第5の実施の形態に係る位置計測装置における移動量推定処理ルーチンの内容を示すフローチャートである。
【図16】本発明の第6の実施の形態に係る位置計測装置における移動量推定処理ルーチンの内容を示すフローチャートである。
【図17】本発明の第7の実施の形態に係る位置計測装置における位置計測処理ルーチンの内容を示すフローチャートである。
【符号の説明】
【0172】
10、210、310、410、510 位置計測装置
12、312 運動計測部
14、214 移動量計測部
16 位置推定部
20 画像撮像部
24 運動推定部
26 GPS
28、228、328、428、528 移動量推定部
226 車速センサ
312、412、512 運動移動量計測部
528 移動量推定部
Np パルス数
【特許請求の範囲】
【請求項1】
複数の衛星からなる衛星群から複数の信号を受信し、受信した複数の信号に基づいて、移動体の位置を測定する測位手段と、
前記測位手段によって同一の衛星群から受信した複数の信号に基づいて測位された前記移動体の複数の位置に基づいて、前記移動体の移動量を算出する移動量算出手段と、
を含む移動量計測装置。
【請求項2】
前記移動体の外部を撮像する単眼の撮像手段を更に含み、
前記移動量算出手段は、第1時刻と該第1時刻より後の第2時刻とにおいて同一の衛星群から複数の信号を受信した場合、前記測位手段によって前記第1時刻に測位された前記移動体の位置と、前記測位手段によって前記第2時刻に測位された前記移動体の位置とに基づいて、前記第1時刻から前記第2時刻までの期間における前記移動体の移動量を算出し、前記第1時刻と前記第2時刻とにおいて異なる衛星群から複数の信号を受信した場合、前記撮像手段によって前記第1時刻に撮像された画像と前記第2時刻に撮像された画像とに基づいて、前記第1時刻から前記第2時刻までの期間における前記移動体の移動量を算出する請求項1記載の移動量計測装置。
【請求項3】
前記移動体の速度に応じてパルス信号を出力する速度センサを更に含み、
前記移動量算出手段は、第1時刻と該第1時刻より後の第2時刻とにおいて同一の衛星群から複数の信号を受信した場合、前記測位手段によって前記第1時刻に測位された前記移動体の位置と、前記測位手段によって前記第2時刻に測位された前記移動体の位置とに基づいて、前記第1時刻から前記第2時刻までの期間における前記移動体の移動量を算出し、前記第1時刻と前記第2時刻とにおいて異なる衛星群から複数の信号を受信した場合、前記第1時刻から前記第2時刻までの期間に出力されたパルス信号に基づいて、前記第1時刻から前記第2時刻までの期間における前記移動体の移動量を算出する請求項1記載の移動量計測装置。
【請求項4】
前記移動体の外部を撮像する単眼の撮像手段と、
前記移動体の速度に応じてパルス信号を出力する速度センサとを更に含み、
前記移動量算出手段は、第1時刻と該第1時刻より後の第2時刻とにおいて同一の衛星群から複数の信号を受信した場合、前記測位手段によって前記第1時刻に測位された前記移動体の位置と、前記測位手段によって前記第2時刻に測位された前記移動体の位置とに基づいて、前記第1時刻から前記第2時刻までの期間における前記移動体の移動量を算出し、前記第1時刻と前記第2時刻とにおいて異なる衛星群から複数の信号を受信した場合であって、かつ、前記第1時刻から前記第2時刻までの期間に出力されたパルス信号のパルス数が、所定数以上である場合、前記第1時刻から前記第2時刻までの期間に出力されたパルス信号に基づいて、前記移動体の移動量を算出し、前記第1時刻と前記第2時刻とにおいて異なる衛星群から複数の信号を受信した場合であって、かつ、前記第1時刻から前記第2時刻までの期間に出力されたパルス信号のパルス数が、所定数未満である場合、前記撮像手段によって前記第1時刻に撮像された画像と前記第2時刻に撮像された画像とに基づいて、前記第1時刻から前記第2時刻までの期間における前記移動体の移動量を算出する請求項1記載の移動量計測装置。
【請求項5】
移動体の外部を撮像する単眼の撮像手段と、
前記移動体の速度に応じてパルス信号を出力する速度センサと、
第1時刻から該第1時刻より後の第2時刻までの期間に出力されたパルス信号のパルス数が、所定数以上である場合、前記第1時刻から前記第2時刻までの期間に出力されたパルス信号に基づいて、前記移動体の移動量を算出し、前記第1時刻から前記第2時刻までの期間に出力されたパルス信号のパルス数が、所定数未満である場合、前記撮像手段によって前記第1時刻に撮像された画像と前記第2時刻に撮像された画像とに基づいて、前記第1時刻から前記第2時刻までの期間における前記移動体の移動量を算出する移動量算出手段と、
を含む移動量計測装置。
【請求項6】
複数の衛星からなる衛星群から複数の信号を受信し、受信した複数の信号に基づいて、移動体の位置を測定する測位手段と、
前記移動体の外部を撮像する単眼の撮像手段と、
前記撮像手段によって撮像された複数の画像に基づいて、前記画像の特徴点までの距離を算出する距離算出手段と、
前記移動体の速度に応じてパルス信号を出力する速度センサと、
前記測位手段によって測位された前記移動体の位置、前記距離算出手段によって算出された前記距離、及び前記速度センサから出力されたパルス信号に基づいて、前記移動体の移動量を算出する移動量算出手段と、
を含む移動量計測装置。
【請求項7】
請求項1〜請求項6の何れか1項記載の移動量計測装置と、
前記移動体の運動の3軸角速度の少なくとも1つを推定する運動推定手段と、
前記移動量計測装置によって算出された前記移動体の移動量と、前記運動推定手段によって推定された前記移動体の運動の前記3軸角速度の少なくとも1つとに基づいて、前記移動体の位置を計測する位置計測手段と、
を含む位置計測装置。
【請求項8】
前記運動推定手段は、前記移動体の運動の3軸角速度の少なくとも1つ及び並進方向を推定し、
前記位置計測手段は、前記移動量計測装置によって算出された前記移動体の移動量と、前記運動推定手段によって推定された前記移動体の運動の前記3軸角速度の少なくとも1つ及び前記並進方向とに基づいて、前記移動体の位置を計測する請求項7記載の位置計測装置。
【請求項1】
複数の衛星からなる衛星群から複数の信号を受信し、受信した複数の信号に基づいて、移動体の位置を測定する測位手段と、
前記測位手段によって同一の衛星群から受信した複数の信号に基づいて測位された前記移動体の複数の位置に基づいて、前記移動体の移動量を算出する移動量算出手段と、
を含む移動量計測装置。
【請求項2】
前記移動体の外部を撮像する単眼の撮像手段を更に含み、
前記移動量算出手段は、第1時刻と該第1時刻より後の第2時刻とにおいて同一の衛星群から複数の信号を受信した場合、前記測位手段によって前記第1時刻に測位された前記移動体の位置と、前記測位手段によって前記第2時刻に測位された前記移動体の位置とに基づいて、前記第1時刻から前記第2時刻までの期間における前記移動体の移動量を算出し、前記第1時刻と前記第2時刻とにおいて異なる衛星群から複数の信号を受信した場合、前記撮像手段によって前記第1時刻に撮像された画像と前記第2時刻に撮像された画像とに基づいて、前記第1時刻から前記第2時刻までの期間における前記移動体の移動量を算出する請求項1記載の移動量計測装置。
【請求項3】
前記移動体の速度に応じてパルス信号を出力する速度センサを更に含み、
前記移動量算出手段は、第1時刻と該第1時刻より後の第2時刻とにおいて同一の衛星群から複数の信号を受信した場合、前記測位手段によって前記第1時刻に測位された前記移動体の位置と、前記測位手段によって前記第2時刻に測位された前記移動体の位置とに基づいて、前記第1時刻から前記第2時刻までの期間における前記移動体の移動量を算出し、前記第1時刻と前記第2時刻とにおいて異なる衛星群から複数の信号を受信した場合、前記第1時刻から前記第2時刻までの期間に出力されたパルス信号に基づいて、前記第1時刻から前記第2時刻までの期間における前記移動体の移動量を算出する請求項1記載の移動量計測装置。
【請求項4】
前記移動体の外部を撮像する単眼の撮像手段と、
前記移動体の速度に応じてパルス信号を出力する速度センサとを更に含み、
前記移動量算出手段は、第1時刻と該第1時刻より後の第2時刻とにおいて同一の衛星群から複数の信号を受信した場合、前記測位手段によって前記第1時刻に測位された前記移動体の位置と、前記測位手段によって前記第2時刻に測位された前記移動体の位置とに基づいて、前記第1時刻から前記第2時刻までの期間における前記移動体の移動量を算出し、前記第1時刻と前記第2時刻とにおいて異なる衛星群から複数の信号を受信した場合であって、かつ、前記第1時刻から前記第2時刻までの期間に出力されたパルス信号のパルス数が、所定数以上である場合、前記第1時刻から前記第2時刻までの期間に出力されたパルス信号に基づいて、前記移動体の移動量を算出し、前記第1時刻と前記第2時刻とにおいて異なる衛星群から複数の信号を受信した場合であって、かつ、前記第1時刻から前記第2時刻までの期間に出力されたパルス信号のパルス数が、所定数未満である場合、前記撮像手段によって前記第1時刻に撮像された画像と前記第2時刻に撮像された画像とに基づいて、前記第1時刻から前記第2時刻までの期間における前記移動体の移動量を算出する請求項1記載の移動量計測装置。
【請求項5】
移動体の外部を撮像する単眼の撮像手段と、
前記移動体の速度に応じてパルス信号を出力する速度センサと、
第1時刻から該第1時刻より後の第2時刻までの期間に出力されたパルス信号のパルス数が、所定数以上である場合、前記第1時刻から前記第2時刻までの期間に出力されたパルス信号に基づいて、前記移動体の移動量を算出し、前記第1時刻から前記第2時刻までの期間に出力されたパルス信号のパルス数が、所定数未満である場合、前記撮像手段によって前記第1時刻に撮像された画像と前記第2時刻に撮像された画像とに基づいて、前記第1時刻から前記第2時刻までの期間における前記移動体の移動量を算出する移動量算出手段と、
を含む移動量計測装置。
【請求項6】
複数の衛星からなる衛星群から複数の信号を受信し、受信した複数の信号に基づいて、移動体の位置を測定する測位手段と、
前記移動体の外部を撮像する単眼の撮像手段と、
前記撮像手段によって撮像された複数の画像に基づいて、前記画像の特徴点までの距離を算出する距離算出手段と、
前記移動体の速度に応じてパルス信号を出力する速度センサと、
前記測位手段によって測位された前記移動体の位置、前記距離算出手段によって算出された前記距離、及び前記速度センサから出力されたパルス信号に基づいて、前記移動体の移動量を算出する移動量算出手段と、
を含む移動量計測装置。
【請求項7】
請求項1〜請求項6の何れか1項記載の移動量計測装置と、
前記移動体の運動の3軸角速度の少なくとも1つを推定する運動推定手段と、
前記移動量計測装置によって算出された前記移動体の移動量と、前記運動推定手段によって推定された前記移動体の運動の前記3軸角速度の少なくとも1つとに基づいて、前記移動体の位置を計測する位置計測手段と、
を含む位置計測装置。
【請求項8】
前記運動推定手段は、前記移動体の運動の3軸角速度の少なくとも1つ及び並進方向を推定し、
前記位置計測手段は、前記移動量計測装置によって算出された前記移動体の移動量と、前記運動推定手段によって推定された前記移動体の運動の前記3軸角速度の少なくとも1つ及び前記並進方向とに基づいて、前記移動体の位置を計測する請求項7記載の位置計測装置。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【公開番号】特開2009−74861(P2009−74861A)
【公開日】平成21年4月9日(2009.4.9)
【国際特許分類】
【出願番号】特願2007−242705(P2007−242705)
【出願日】平成19年9月19日(2007.9.19)
【出願人】(000003609)株式会社豊田中央研究所 (4,200)
【Fターム(参考)】
【公開日】平成21年4月9日(2009.4.9)
【国際特許分類】
【出願日】平成19年9月19日(2007.9.19)
【出願人】(000003609)株式会社豊田中央研究所 (4,200)
【Fターム(参考)】
[ Back to top ]