(57)【要約】車載機１の自車位置推定部１７は、ライダ２によるランドマークの計測結果である点群データを取得すると共に、地図ＤＢ１０に含まれるランドマークの地物情報を取得する。そして、自車位置推定部１７は、取得した点群データに基づき特定したランドマークのサイズと地物情報が示すランドマークのサイズとの差分に基づいて、ライダ２によるランドマークの計測精度を示す信頼度情報を算出する。
2020085062.tif 000002
【特許請求の範囲】
【請求項１】
計測部による地物の計測結果を取得する第１取得部と、
地図データに含まれる前記地物の地物情報を取得する第２取得部と、
計測した地物のサイズと前記地物情報が示す地物のサイズとの差分に基づいて、前記計測部による前記地物の計測精度を示す精度情報を算出する算出部と、を有する計測精度算出装置。
【請求項２】
前記算出部は、前記差分が大きいほど、前記計測精度が低い精度情報を生成する請求項１に記載の計測精度算出装置。
【請求項３】
前記算出部は、前記差分と、移動体の進行方向と、前記地物情報が示す前記地物の法線方向とに基づき、前記移動体の進行方向と横方向のそれぞれに対する前記地物の計測精度を示す精度情報を生成する請求項１または２に記載の計測精度算出装置。
【請求項４】
前記算出部は、前記進行方向と前記地物の法線方向との角度差が大きいほど、前記横方向に対する前記計測精度への前記差分の影響度合いを低くし、前記角度差が小さいほど、前記進行方向に対する前記計測精度への前記差分の影響度合いを低くする請求項３に記載の計測精度算出装置。
【請求項５】
予測された自己位置を示す予測位置情報を取得する第３取得部と、
前記精度情報に基づき、前記予測された自己位置を補正する補正部と、
をさらに備える請求項１〜４のいずれか一項に記載の計測精度算出装置。
【請求項６】
前記補正部は、移動体から前記地物までの前記計測部による計測距離と、前記地物情報に含まれる前記地物の位置情報に基づき予測された前記移動体から前記地物までの予測距離と、の差分値により前記予測された自己位置を補正する場合の前記差分値に対する利得を、前記精度情報に基づき決定する請求項５に記載の計測精度算出装置。
【請求項７】
計測部による地物の計測結果を取得する第１取得部と、
地図データに含まれる前記地物の地物情報を取得する第２取得部と、
予測された自己位置を示す予測位置情報を取得する第３取得部と、
計測した地物のサイズと前記地物情報が示す地物のサイズとの差分に基づいて、前記予測された自己位置を補正する補正部と、
を備える自己位置推定装置。
【請求項８】
計測精度算出装置が実行する制御方法であって、
計測部による地物の計測結果を取得する第１取得工程と、
地図データに含まれる前記地物の地物情報を取得する第２取得工程と、
計測した地物のサイズと前記地物情報が示す地物のサイズとの差分に基づいて、前記計測部による前記地物の計測精度を示す精度情報を算出する算出工程と、
を有する制御方法。
【請求項９】
コンピュータが実行するプログラムであって、
計測部による地物の計測結果を取得する第１取得部と、
地図データに含まれる前記地物の地物情報を取得する第２取得部と、
計測した地物のサイズと前記地物情報が示す地物のサイズとの差分に基づいて、前記計測部による前記地物の計測精度を示す精度情報を算出する算出部
として前記コンピュータを機能させるプログラム。
【請求項１０】
請求項９に記載のプログラムを記憶した記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、計測部による計測の精度を算出する技術に関する。
【背景技術】
【０００２】
従来から、車両の進行先に設置される地物をレーダやカメラを用いて検出し、その検出結果に基づいて自車位置を校正する技術が知られている。例えば、特許文献１には、計測センサの出力と、予め地図上に登録された地物の位置情報とを照合させることで自己位置を推定する技術が開示されている。また、特許文献２には、カルマンフィルタを用いた自車位置推定技術が開示されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１３−２５７７４２号公報
【特許文献２】特開２０１７−７２４２２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
自己位置推定処理において計測対象となる地物が他の地物と隣接している場合には、計測対象の地物と共に他の地物が計測データに含まれてしまい、計測結果の精度（正確性）が低下してしまうことがある。また、計測対象となる地物の一部にオクルージョンが発生した場合にも同様に、本来計測されるべき地物の一部が計測できなくなるため、計測結果の精度が低下してしまうことがある。【０００５】
本発明は、上記のような課題を解決するためになされたものであり、計測部による地物の計測精度を好適に算出可能な計測精度算出装置及び自己位置推定装置を提供することを主な目的とする。
【課題を解決するための手段】
【０００６】
請求項に記載の発明は、計測精度算出装置であって、計測部による地物の計測結果を取得する第１取得部と、地図データに含まれる前記地物の地物情報を取得する第２取得部と、計測した地物のサイズと前記地物情報が示す地物のサイズとの差分に基づいて、前記計測部による前記地物の計測精度を示す精度情報を算出する算出部と、を有する。【０００７】
また、請求項に記載の発明は、自己位置推定装置であって、計測部による地物の計測結果を取得する第１取得部と、地図データに含まれる前記地物の地物情報を取得する第２取得部と、予測された自己位置を示す予測位置情報を取得する第３取得部と、計測した地物のサイズと前記地物情報が示す地物のサイズとの差分に基づいて、前記予測された自己位置を補正する補正部と、を備える。【０００８】
また、請求項に記載の発明は、計測精度算出装置が実行する制御方法であって、計測部による地物の計測結果を取得する第１取得工程と、地図データに含まれる前記地物の地物情報を取得する第２取得工程と、計測した地物のサイズと前記地物情報が示す地物のサイズとの差分に基づいて、前記計測部による前記地物の計測精度を示す精度情報を算出する算出工程と、を有する。【０００９】
また、請求項に記載の発明は、コンピュータが実行するプログラムであって、計測部による地物の計測結果を取得する第１取得部と、地図データに含まれる前記地物の地物情報を取得する第２取得部と、計測した地物のサイズと前記地物情報が示す地物のサイズとの差分に基づいて、前記計測部による前記地物の計測精度を示す精度情報を算出する算出部として前記コンピュータを機能させる。
【図面の簡単な説明】
【００１０】
【図１】運転支援システムの概略構成図である。
【図２】車載機の機能的構成を示すブロック図である。
【図３】地図データベースのデータ構造の一例を示す。
【図４】状態変数ベクトルを２次元直交座標で表した図である。
【図５】予測ステップと計測更新ステップとの概略的な関係を示す図である。
【図６】自車位置推定部の機能ブロックを示す。
【図７】信頼度情報のグラフを示す。
【図８】オクルージョンが発生した場合のライダのレーザ光の照射範囲を概略的に示した図である。
【図９】ランドマークとなる地物に隣接して他の地物が存在する場合のライダのレーザ光の照射範囲を概略的に示した図である。
【図１０】２次元のワールド座標系及び車両座標系により表された車両とランドマークとの位置関係を示した図である。
【図１１】車両から見てランドマークがほぼ正面を向いているときの信頼度情報の２次元グラフである。
【図１２】車両から見てランドマークがやや横向きに近いときの信頼度情報の２次元グラフである。
【図１３】自車位置推定処理を示すフローチャートである。
【図１４】変形例１に係る信頼度情報のグラフを示す。
【図１５】変形例２に係る信頼度情報のグラフを示す。
【発明を実施するための形態】
【００１１】
本発明の好適な実施形態によれば、計測精度算出装置は、計測部による地物の計測結果を取得する第１取得部と、地図データに含まれる前記地物の地物情報を取得する第２取得部と、計測した地物のサイズと前記地物情報が示す地物のサイズとの差分に基づいて、前記計測部による前記地物の計測精度を示す精度情報を算出する算出部と、を有する。この態様により、計測精度算出装置は、計測対象の地物にオクルージョンが発生して一部の計測データが取得できなかった場合又は計測対象の地物に隣接する他の地物の計測データを計測対象の地物の計測データの一部として取得した場合のいずれにおいても、計測部による地物の計測精度を示す精度情報を的確に算出することができる。
【００１２】
上記計測精度算出装置の一態様では、前記算出部は、前記差分が大きいほど、前記計測精度が低い精度情報を生成する。このようにすることで、計測対象の地物にオクルージョンが発生して一部の計測データが取得できなかった場合又は隣接地物の計測データを計測対象の地物の計測データの一部として取得した場合のいずれにおいても、計測精度が低いことを示す精度情報を好適に生成することができる。
【００１３】
上記計測精度算出装置の一態様では、前記算出部は、前記差分と、移動体の進行方向と、前記地物情報が示す前記地物の法線方向とに基づき、前記移動体の進行方向と横方向のそれぞれに対する前記地物の計測精度を示す精度情報を生成する。一般に、移動体の進行方向に対する地物の向きによって、計測誤差が生じにくい方向と計測誤差が生じやすい方向とがそれぞれ生じる場合がある。従って、この態様によれば、計測精度算出装置は、移動体の進行方向と横方向のそれぞれに対する地物の計測精度を示す精度情報を好適に生成することができる。
【００１４】
上記計測精度算出装置の他の一態様では、前記算出部は、前記進行方向と前記地物の法線方向との角度差が大きいほど、前記横方向に対する前記計測精度への前記差分の影響度合いを低くし、前記角度差が小さいほど、前記進行方向に対する前記計測精度への前記差分の影響度合いを低くする。ここでの地物の法線方向は、地物が形成する平面に対し垂直となる地物の正面方向及び背面方向のうち、移動体の進行方向と９０度差以内となる方向を指すものとする。この場合、移動体の進行方向と地物の法線方向との角度差が大きいほど、地物は移動体に対して横向きとなるため、移動体の横方向に対する計測誤差が生じにくくなる。一方、移動体の進行方向と地物の法線方向との角度差が小さいほど、地物は移動体に対して正面向きとなるため、移動体の進行方向に対する計測誤差が生じにくくなる。よって、計測精度算出装置は、この態様により、移動体の進行方向と横方向のそれぞれに対する地物の計測精度を的確に表す精度情報を生成することができる。
【００１５】
上記計測精度算出装置の他の一態様では、計測精度算出装置は、予測された自己位置を示す予測位置情報を取得する第３取得部と、前記精度情報に基づき、前記予測された自己位置を補正する補正部と、をさらに備える。この態様により、計測精度算出装置は、生成した精度情報を反映した自己位置推定を好適に実行することができる。
【００１６】
上記計測精度算出装置の他の一態様では、前記補正部は、移動体から前記地物までの前記計測部による計測距離と、前記地物情報に含まれる前記地物の位置情報に基づき予測された前記移動体から前記地物までの予測距離と、の差分値により前記予測された自己位置を補正する場合の前記差分値に対する利得を、前記精度情報に基づき決定する。この態様により、計測精度算出装置は、予測された自己位置を計測距離と予測距離との差分値により補正する場合に精度情報に基づきその利得を調整することで、不正確な補正を防止して自己位置推定精度を好適に向上させることができる。
【００１７】
本発明の他の好適な実施形態によれば、自己位置推定装置であって、計測部による地物の計測結果を取得する第１取得部と、地図データに含まれる前記地物の地物情報を取得する第２取得部と、予測された自己位置を示す予測位置情報を取得する第３取得部と、計測した地物のサイズと前記地物情報が示す地物のサイズとの差分に基づいて、前記予測された自己位置を補正する補正部と、を備える。上述の差分は、計測対象の地物にオクルージョンが発生して一部の計測データが取得できなかった場合又は計測対象の地物に隣接する他の地物の計測データを計測対象の地物の計測データの一部として取得した場合のいずれにおいても発生する。従って、自己位置推定装置は、この態様により、予測された自己位置を上述の差分を考慮して好適に補正して最終的な自己位置を決定することができる。
【００１８】
本発明の他の好適な実施形態によれば、計測精度算出装置が実行する制御方法であって、計測部による地物の計測結果を取得する第１取得工程と、地図データに含まれる前記地物の地物情報を取得する第２取得工程と、計測した地物のサイズと前記地物情報が示す地物のサイズとの差分に基づいて、前記計測部による前記地物の計測精度を示す精度情報を算出する算出工程と、を有する。計測精度算出装置は、この制御方法を実行することで、計測対象の地物にオクルージョンが発生して一部の計測データが取得できなかった場合又は計測対象の地物に隣接する他の地物の計測データを計測対象の地物の計測データの一部として取得した場合のいずれにおいても、計測部による地物の計測精度を示す精度情報を的確に算出することができる。
【００１９】
本発明の他の好適な実施形態によれば、コンピュータが実行するプログラムであって、計測部による地物の計測結果を取得する第１取得部と、地図データに含まれる前記地物の地物情報を取得する第２取得部と、計測した地物のサイズと前記地物情報が示す地物のサイズとの差分に基づいて、前記計測部による前記地物の計測精度を示す精度情報を算出する算出部として前記コンピュータを機能させる。コンピュータは、このプログラムを実行することで、計測対象の地物にオクルージョンが発生して一部の計測データが取得できなかった場合又は計測対象の地物に隣接する他の地物の計測データを計測対象の地物の計測データの一部として取得した場合のいずれにおいても、計測部による地物の計測精度を示す精度情報を的確に算出することができる。好適には、上記プログラムは、記憶媒体に記憶される。
【実施例】
【００２０】
以下、図面を参照して本発明の好適な実施例について説明する。なお、任意の記号の上に「＾」または「−」が付された文字を、本明細書では便宜上、「Ａ＾」または「Ａ−」（「Ａ」は任意の文字）と表す。【００２１】
［運転支援システムの概要］
図１は、本実施例に係る運転支援システムの概略構成である。図１に示す運転支援システムは、車両に搭載され、車両の運転支援に関する制御を行う車載機１と、ライダ（Ｌｉｄａｒ：Ｌｉｇｈｔ  Ｄｅｔｅｃｔｉｏｎ  ａｎｄ  Ｒａｎｇｉｎｇ、または、Ｌａｓｅｒ  Ｉｌｌｕｍｉｎａｔｅｄ  Ｄｅｔｅｃｔｉｏｎ  Ａｎｄ  Ｒａｎｇｉｎｇ）２と、ジャイロセンサ３と、車速センサ４と、ＧＰＳ受信機５とを有する。なお、以後において、「地図」とは、従来の経路案内用の車載機が参照するデータに加えて、ＡＤＡＳ（Ａｄｖａｎｃｅｄ  Ｄｒｉｖｅｒ  Ａｓｓｉｓｔａｎｃｅ  Ｓｙｓｔｅｍ）や自動運転に用いられるデータも含むものとする。【００２２】
車載機１は、ライダ２、ジャイロセンサ３、車速センサ４、及びＧＰＳ受信機５と電気的に接続し、これらの出力に基づき、車載機１が搭載される車両の位置（「自車位置」とも呼ぶ。）の推定を行う。そして、車載機１は、自車位置の推定結果に基づき、設定された目的地への経路に沿って走行するように、車両の自動運転制御などを行う。車載機１は、道路データ及び道路付近に設けられた目印となる地物に関する情報である地物情報を記憶した地図データベース（ＤＢ：ＤａｔａＢａｓｅ）１０を記憶する。上述の目印となる地物は、例えば、道路脇に周期的に並んでいるキロポスト、１００ｍポスト、デリニエータ、交通インフラ設備（例えば標識、方面看板、信号）、電柱、街灯などの地物である。そして、車載機１は、この地物情報に基づき、ライダ２等の出力と照合させて自車位置の推定を行う。車載機１は、本発明における「計測精度算出装置」及び「自己位置推定装置」の一例である。尚、地図ＤＢ１０は車載装置１に記憶される代わりに、外部サーバ等の外部の記憶装置に記憶されてもよい。その場合、車載装置１は無線通信等を介して外部の記憶装置から地図ＤＢ１０の少なくとも一部を取得する。【００２３】
ライダ２は、水平方向および垂直方向の所定の角度範囲に対してパルスレーザを出射することで、外界に存在する物体までの距離を離散的に測定し、当該物体の位置を示す３次元の点群情報を生成する。この場合、ライダ２は、照射方向を変えながらレーザ光を照射する照射部と、照射したレーザ光が物体で反射した反射光（散乱光）を受光する受光部と、受光部が出力する受光信号に基づくスキャンデータを出力する出力部とを有する。スキャンデータは、点群データであり、受光部が受光したレーザ光に対応する照射方向と、上述の受光信号に基づき特定される、その照射方向での物体までの距離とに基づき生成される。ライダ２、ジャイロセンサ３、車速センサ４、ＧＰＳ受信機５は、それぞれ、出力データを車載機１へ供給する。ライダ２は、「計測部」の一例である。【００２４】
図２は、車載機１の機能的構成を示すブロック図である。車載機１は、主に、インターフェース１１と、記憶部１２と、入力部１４と、制御部１５と、情報出力部１６と、を有する。これらの各要素は、バスラインを介して相互に接続されている。【００２５】
インターフェース１１は、ライダ２、ジャイロセンサ３、車速センサ４、及びＧＰＳ受信機５などのセンサから出力データを取得し、制御部１５へ供給する。【００２６】
記憶部１２は、制御部１５が実行するプログラムや、制御部１５が所定の処理を実行するのに必要な情報を記憶する。本実施例では、記憶部１２は、地物情報を含む地図ＤＢ１０を記憶する。図３は、地図ＤＢ１０のデータ構造の一例を示す。図３に示すように、地図ＤＢ１０は、施設情報、道路データ、及び地物情報を含む。【００２７】
地物情報は、地物ごとに当該地物に関する情報が関連付けられた情報であり、ここでは、地物のインデックスに相当する地物ＩＤと、位置情報と、向き情報（法線情報）と、サイズ情報とを含む。位置情報は、緯度及び経度（及び標高）等により表わされた地物の絶対的な位置を示す。向き情報及びサイズ情報は、例えば、看板などの平面的な形状を有している地物に対して設けられる情報である。向き情報は、地物の向きを表す情報であり、例えば地物に形成された面に対する法線ベクトル等を示す。サイズ情報は、地物の大きさを表す情報であり、本実施例では、地物に形成された被計測面の横幅及び縦幅の情報を含んでいる。【００２８】
なお、地図ＤＢ１０は、定期的に更新されてもよい。この場合、例えば、制御部１５は、図示しない通信部を介し、地図情報を管理するサーバ装置から、自車位置が属するエリアに関する部分地図情報を受信し、地図ＤＢ１０に反映させる。【００２９】
入力部１４は、ユーザが操作するためのボタン、タッチパネル、リモートコントローラ、音声入力装置等である。情報出力部１６は、例えば、制御部１５の制御に基づき出力を行うディスプレイやスピーカ等である。【００３０】
制御部１５は、プログラムを実行するＣＰＵなどを含み、車載機１の全体を制御する。本実施例では、制御部１５は、インターフェース１１から供給される各センサの出力信号及び地図ＤＢ１０に基づき、自車位置の推定を行う自車位置推定部１７を有する。そして、制御部１５は、自車位置の推定結果に基づいて自動運転制御を含む車両の運転支援に関する制御などを行う。制御部１５は、本発明における「第１取得部」、「第２取得部」、「第３取得部」、「算出部」、「補正部」、及びプログラムを実行する「コンピュータ」の一例である。【００３１】
［自車位置推定処理の概要］
まず、自車位置推定部１７による自車位置の推定処理の概要について説明する。【００３２】
自車位置推定部１７は、計測対象となる地物（「ランドマーク」とも呼ぶ。）に対するライダ２による距離及び角度の計測値と、地図ＤＢ１０から抽出したランドマークの位置情報とに基づき、ジャイロセンサ３、車速センサ４、及び／又はＧＰＳ受信機５の出力データから推定した自車位置を補正する。本実施例では、一例として、自車位置推定部１７は、ベイズ推定に基づく状態推定手法に基づき、ジャイロセンサ３、車速センサ４等の出力データから自車位置を予測する予測ステップと、直前の予測ステップで算出した自車位置の予測値を補正する計測更新ステップとを交互に実行する。これらのステップで用いる状態推定フィルタは、ベイズ推定を行うように開発された様々のフィルタが利用可能であり、例えば、拡張カルマンフィルタ、アンセンテッドカルマンフィルタ、パーティクルフィルタなどが該当する。このように、ベイズ推定に基づく位置推定は、種々の方法が提案されている。【００３３】
以下では、拡張カルマンフィルタを用いた自車位置推定について簡略的に説明する。【００３４】
図４は、状態変数ベクトルを２次元直交座標で表した図である。図４に示すように、ｘｙの２次元直交座標上で定義された平面での自車位置は、座標「（ｘ、ｙ）」、自車の方位（ヨー角）「ψ」により表される。ここでは、ヨー角ψは、車の進行方向とｘ軸とのなす角として定義されている。座標（ｘ、ｙ）は、例えば緯度及び経度の組合せに相当する絶対位置、あるいは所定地点を原点とした位置を示すワールド座標である。【００３５】
図５は、予測ステップと計測更新ステップとの概略的な関係を示す図である。また、図６は、自車位置推定部１７の機能ブロックの一例を示す。図５に示すように、予測ステップと計測更新ステップとを繰り返すことで、自車位置を示す状態変数ベクトル「Ｘ」の推定値の算出及び更新を逐次的に実行する。また、図６に示すように、自車位置推定部１７は、予測ステップを実行する位置予測部２１と、計測更新ステップを実行する位置推定部２２とを有する。位置予測部２１は、デッドレコニングブロック２３及び位置予測ブロック２４を含み、位置推定部２２は、ランドマーク探索・抽出ブロック２５及び位置補正ブロック２６を含む。なお、図５では、計算対象となる基準時刻（即ち現在時刻）「ｔ」の状態変数ベクトルを、「Ｘ−（ｔ）」または「Ｘ＾（ｔ）」と表記している（「状態変数ベクトルＸ（ｔ）＝（ｘ（ｔ）、ｙ（ｔ）、ψ（ｔ））Ｔ」と表記する）。ここで、予測ステップで推定された暫定的な推定値（予測値）には当該予測値を表す文字の上に「−」を付し、計測更新ステップで更新された，より精度の高い推定値には当該値を表す文字の上に「＾」を付す。【００３６】
予測ステップでは、自車位置推定部１７のデッドレコニングブロック２３は、車両の移動速度「ｖ」と角速度「ω」（これらをまとめて「制御値ｕ（ｔ）＝（ｖ（ｔ）、ω（ｔ））Ｔ」と表記する。）を用い、前回時刻からの移動距離と方位変化を求める。自車位置推定部１７の位置予測ブロック２４は、直前の計測更新ステップで算出された時刻ｔ−１の状態変数ベクトルＸ＾（ｔ−１）に対し、求めた移動距離と方位変化を加えて、時刻ｔの自車位置の予測値（「予測自車位置」とも呼ぶ。）Ｘ−（ｔ）を算出する。また、これと同時に、予測自車位置Ｘ−（ｔ）の誤差分布に相当する共分散行列「Ｐ−（ｔ）」を、直前の計測更新ステップで算出された時刻ｔ−１での共分散行列「Ｐ＾（ｔ−１）」から算出する。【００３７】
計測更新ステップでは、自車位置推定部１７のランドマーク探索・抽出ブロック２５は、地図ＤＢ１０に登録された計測対象の地物（ランドマーク）の位置ベクトルとライダ２のスキャンデータとの対応付けを行う。そして、自車位置推定部１７のランドマーク探索・抽出ブロック２５は、この対応付けができた場合に、対応付けができたランドマークのライダ２による計測値「Ｚ（ｔ）」と、予測自車位置Ｘ−（ｔ）及び地図ＤＢ１０に登録されたランドマークの位置ベクトルを用いてライダ２による計測処理をモデル化して求めたランドマークの計測値（「計測予測値」と呼ぶ。）「Ｚ−（ｔ）」と、をそれぞれ取得する。計測値Ｚ（ｔ）は、時刻ｔにライダ２が計測したランドマークの距離及びスキャン角度から、車両の進行方向と横方向を軸とした成分に変換した車両座標系におけるベクトル値である。そして、自車位置推定部１７の位置補正ブロック２６は、計測値Ｚ（ｔ）と計測予測値Ｚ−（ｔ）との差分値を算出する。【００３８】
なお、実際には、ランドマークは所定の大きさを有し、ライダ２のレーザ光が照射された被照射面上の複数の計測点に対応するデータが点群データとして取得される。よって、例えば、ランドマーク探索・抽出ブロック２５は、対象のランドマークを計測した各計測点における計測値を平均化することで、対象のランドマークの中心位置に相当する計測値Ｚ（ｔ）を算出する。【００３９】
そして、自車位置推定部１７の位置補正ブロック２６は、以下の式（１）に示すように、計測値Ｚ（ｔ）と計測予測値Ｚ−（ｔ）との差分値にカルマンゲイン「Ｋ（ｔ）」を乗算し、これを予測自車位置Ｘ−（ｔ）に加えることで、更新された状態変数ベクトル（「推定自車位置」とも呼ぶ。）Ｘ＾（ｔ）を算出する。【００４０】
【数１】
2020085062.tif 000003
また、計測更新ステップでは、自車位置推定部１７の位置補正ブロック２６は、予測ステップと同様、推定自車位置Ｘ＾（ｔ）の誤差分布に相当する共分散行列Ｐ＾（ｔ）（単にＰ（ｔ）とも表記する）を共分散行列Ｐ−（ｔ）から求める。カルマンゲインＫ（ｔ）等のパラメータについては、例えば拡張カルマンフィルタを用いた公知の自己位置推定技術と同様に算出することが可能である。【００４１】
本実施例では、自車位置推定部１７は、カルマンゲインＫ（ｔ）を、計測値Ｚ（ｔ）（即ちランドマークの中心点までの計測距離）の信頼度を表す指標である後述の信頼度情報に応じて補正する。これにより、オクルージョンの発生やランドマークに隣接する地物の存在等に起因して計測値Ｚ（ｔ）の信頼度が低い場合にはカルマンゲインＫ（ｔ）を低くし、予測自車位置Ｘ−（ｔ）に対する補正量を少なくする。この場合、予測自車位置Ｘ−（ｔ）に対する不正確な補正が好適に抑制され、自己位置推定精度が向上する。【００４２】
なお、自車位置推定部１７は、複数の地物に対し、地図ＤＢ１０に登録された地物の位置ベクトルとライダ２のスキャンデータとの対応付けができた場合、選定した任意の一組の計測予測値及び計測値等に基づき計測更新ステップを行ってもよく、対応付けができた全ての計測予測値及び計測値等に基づき計測更新ステップを複数回行ってもよい。なお、複数の計測予測値及び計測値等を用いる場合には、自車位置推定部１７は、ライダ２から遠い地物ほどライダ計測精度が悪化することを勘案し、ライダ２と地物との距離が長いほど、当該地物に関する重み付けを小さくするとよい。【００４３】
このように、予測ステップと計測更新ステップが繰り返し実施され、予測自車位置Ｘ−（ｔ）と推定自車位置Ｘ＾（ｔ）が逐次的に計算されることにより、もっとも確からしい自車位置が計算される。【００４４】
なお、上記の説明において、計測値Ｚ（ｔ）は本発明の「計測距離」の一例であり、計測予測値Ｚ−（ｔ）は本発明の「予測距離」の一例である。【００４５】
［信頼度情報の算出］
次に、計測値Ｚ（ｔ）の信頼度を表す指標である信頼度情報の算出方法について説明する。概略的には、自車位置推定部１７は、地物情報に含まれるサイズ情報が示すランドマークの地図上のサイズと、ライダ２により計測されたランドマークの計測点が分布する範囲（計測範囲）のサイズとの差分を算出し、当該差分を正規化することで、信頼度情報を算出する。これにより、自車位置推定部１７は、計測値Ｚ（ｔ）の信頼度を的確に表した信頼度情報を算出する。【００４６】
（１）方向によらない信頼度情報の算出
まず、計測値Ｚ（ｔ）に対して方向によらない信頼度情報を生成する場合について説明する。【００４７】
この場合、自車位置推定部１７は、地物情報に含まれるランドマークの地図上の横幅「ＷＭ」及び縦幅（高さ方向の幅）「ＨＭ」と、時刻ｔにおけるライダ２の点群データに基づき特定されたランドマークの計測範囲の横幅「ＷＬ（ｔ）」及び縦幅「ＨＬ（ｔ）」と、感度係数「ｋ」とに基づき、以下の式（２）に示すように信頼度情報ａ（ｔ）を定める。【００４８】
【数２】
2020085062.tif 000004
上述の式（２）では、信頼度情報ａ（ｔ）は、０から１までの値域となるように正規化され、計測された横幅ＷＬ（ｔ）及び縦幅ＨＬ（ｔ）が地図上の横幅ＷＭ及び縦幅ＨＭと近いほど最大値の１に近づき、横幅ＷＬ（ｔ）及び縦幅ＨＬ（ｔ）が横幅ＷＭ及び縦幅ＨＭから離れるほど０に近づく。なお、信頼度情報ａ（ｔ）は、１に近いほど計測値Ｚ（ｔ）に対する信頼度が高い（即ち計測値Ｚ（ｔ）が正確）であることを示す。また、感度係数ｋは、計測された横幅ＷＬ（ｔ）及び縦幅ＨＬ（ｔ）と地図上の横幅ＷＭ及び縦幅ＨＭとの差分に対する信頼度情報ａ（ｔ）への感度を調整する係数であり、感度係数ｋが大きいほど算出される信頼度情報ａ（ｔ）が小さくなる。なお、感度係数ｋは、例えば予め定められた値に設定される。【００４９】
図７（Ａ）は、横幅ＷＭが６０ｃｍ、縦幅ＨＭが６０ｃｍとなるランドマークに対して縦幅ＨＬ（ｔ）を６０ｃｍに固定した場合の横幅ＷＬ（ｔ）と式（２）に基づく信頼度情報ａ（ｔ）との関係を示す２次元グラフである。また、図７（Ｂ）は、図７（Ａ）と同一のランドマークに対して横幅ＷＬ（ｔ）を６０ｃｍに固定した場合の縦幅ＨＬ（ｔ）と信頼度情報ａ（ｔ）との関係を示す２次元グラフである。さらに、図７（Ｃ）は、図７（Ａ）、（Ｂ）と同一のランドマークに対する横幅ＷＬ（ｔ）と縦幅ＨＬ（ｔ）と信頼度情報ａ（ｔ）との関係を示す３次元グラフである。また、図７（Ｄ）は、図７（Ａ）における信頼度情報ａ（ｔ）の逆数「１／ａ（ｔ）」を示したグラフである。【００５０】
図７（Ａ）及び図７（Ｃ）に示すように、横幅ＷＬ（ｔ）が横幅ＷＭである６０ｃｍに近づくほど、信頼度情報ａ（ｔ）は最大値である１に近づく。一方、図７（Ｂ）及び図７（Ｃ）に示すように、縦幅ＨＬ（ｔ）が縦幅ＨＭである６０ｃｍに近づくほど、信頼度情報ａ（ｔ）は最大値である１に近づく。なお、図７（Ｄ）に示すように、信頼度情報ａ（ｔ）の逆数は、横幅ＷＬ（ｔ）が横幅ＷＭに近づくほど最小値である１に近づき、横幅ＷＬ（ｔ）が横幅ＷＭから遠ざかるほど大きい値となる。同様に、信頼度情報ａ（ｔ）の逆数は、縦幅ＨＬ（ｔ）が縦幅ＨＭに近づくほど最小値である１に近づき、縦幅ＨＬ（ｔ）が縦幅ＨＭから遠ざかるほど大きい値となる。信頼度情報ａ（ｔ）の逆数は、後述するように、信頼度情報に応じてカルマンゲインを設定する際に用いられる。【００５１】
ここで、式（２）の妥当性について、図８及び図９に示す具体例を用いて補足説明する。【００５２】
図８（Ａ）は、ランドマークとなる地物２０のライダ２による計測時においてオクルージョンが発生した場合のライダ２のレーザ光の照射範囲を概略的に示した図であり、図８（Ｂ）は、ライダ２のレーザ光が照射された被照射面上の計測点を明示した地物２０の平面図である。なお、ここではレーザ光の反射強度が高い道路標識を地物２０の例としている。【００５３】
自車位置推定部１７は、ライダ２の計測範囲内にランドマークとなる地物２０が存在することを認識し、当該地物２０が存在することが予測される範囲を示すウィンドウ２１を設定する。例えば、自車位置推定部１７は、予測自車位置Ｘ−（ｔ）から所定距離以内の位置を示す位置情報が地物情報に含まれる地物２０をランドマークとして特定し、当該地物２０の位置情報が示す位置を中心とする所定サイズの矩形領域をウィンドウ２１として設定する。なお、この場合、自車位置推定部１７は、地物２０のサイズ情報をさらに参照することで、ウィンドウ２１の大きさを決定してもよい。そして、自車位置推定部１７は、ウィンドウ２１内の位置を示す計測点であって、反射強度が所定の閾値以上となる計測点の点群データを、地物２０の点群データとみなす。【００５４】
この場合、図８（Ｂ）に示すように、ウィンドウ２１内のライダ２の点群データに基づき特定する地物２０の横幅ＷＬ（ｔ）は、地図上の地物２０の横幅ＷＭより短くなる。よって、この場合、自車位置推定部１７は、式（２）に基づき、信頼度情報ａ（ｔ）を１より小さい値に設定する（図７（Ａ）及び図７（Ｃ）参照）。【００５５】
一方、自車位置推定部１７は、ウィンドウ２１内の各計測点の計測値を平均化することで、推定自車位置Ｘ＾（ｔ）の算出（式（１）参照）に用いる地物２０の計測値Ｚ（ｔ）を算出する。ここで、地物２０の左側の一部にオクルージョンが発生しているため、上述の平均化により算出される計測値Ｚ（ｔ）は、本来算出されるべき値（即ち地物２０の中心位置を示す値）よりも右側にずれた値となる。【００５６】
このように、式（２）によれば、自車位置推定部１７は、ランドマークにオクルージョンが発生し、計測値Ｚ（ｔ）にずれが生じると推測される場合に、信頼度情報ａ（ｔ）を低い値に設定することができる。これにより、後述するように、自車位置推定部１７は、計測値Ｚ（ｔ）の信頼度が低いときの予測自車位置Ｘ−（ｔ）に対する補正量を少なくし、不正確な補正防止によって自己位置推定精度を好適に向上させる。【００５７】
図９（Ａ）は、ランドマークとなる地物２０に隣接して他の地物２５が存在する場合のライダ２のレーザ光の照射範囲を概略的に示した図であり、図９（Ｂ）は、図９（Ａ）においてライダ２のレーザ光が照射された被照射面上の計測点を明示した地物２０の平面図である。【００５８】
図９（Ａ）の例では、自車位置推定部１７は、ウィンドウ２１を設定後、ウィンドウ２１内において所定閾値以上の反射強度となる計測点の点群データを、地物２０の点群データとみなす。ここで、ウィンドウ２１には、地物２０の計測点に加えて、地物２０に隣接する地物２５の一部の計測点も含まれている。【００５９】
この場合、図９（Ｂ）に示すように、ウィンドウ２１内のライダ２の点群データに基づき特定される横幅ＷＬ（ｔ）は、地図上の地物２０の横幅ＷＭより長くなる。よって、この場合、自車位置推定部１７は、横幅ＷＬ（ｔ）と地図上の横幅ＷＭとの差異が大きくなることから、式（２）に基づき、信頼度情報ａ（ｔ）を１より小さい値に設定する（図７（Ａ）及び図７（Ｃ）参照）。【００６０】
一方、自車位置推定部１７は、ウィンドウ２１内の各計測点の計測値を平均化することで、推定自車位置Ｘ＾（ｔ）の算出（式（１）参照）に用いる地物２０の計測値Ｚ（ｔ）を算出する。ここで、ウィンドウ２１内の点群データには、地物２０を計測した点群データに加えて、地物２０に隣接する地物２５を計測した点群データも含まれているため、上述の平均化により算出される計測値Ｚ（ｔ）は、本来算出されるべき値（即ち地物２０の中心位置を示す値）よりも左側にずれた値となる。【００６１】
このように、式（２）によれば、計測対象の地物に他の地物が隣接することに起因して計測値Ｚ（ｔ）にずれが生じた場合においても、自車位置推定部１７は、信頼度情報ａ（ｔ）を低い値に設定することができる。これにより、後述するように、自車位置推定部１７は、計測値Ｚ（ｔ）の信頼度が低いときの予測自車位置Ｘ−（ｔ）に対する補正量を少なくし、不正確な補正防止によって自己位置推定精度を好適に向上させる。【００６２】
なお、ランドマークとなる地物は、再帰性反射材が使用されているものが多く、ライダ２に戻ってくる反射光の強度が高いため、反射強度に対して閾値を用いたフィルタリングにより、点群データを抽出しやすい対象と言える。一方、走行環境において、視線誘導標や他車両の反射板等の反射強度の高い物体が多数存在し、それらの点群データをランドマークの点群データと誤って検出してしまう可能性がある。このような場合、横幅ＷＬ（ｔ）が横幅ＷＭよりも長くなる、又は／及び、縦幅ＨＬ（ｔ）が縦幅ＨＭよりも長くなる。このような場合において、本実施例では、自車位置推定部１７は、式（２）に基づき、計測値の信頼度が低いことを示す信頼度情報ａ（ｔ）を適切に設定することができる。【００６３】
（２）方向毎の信頼度情報の算出
次に、車両の進行方向と横方向のそれぞれに対して信頼度情報を算出する例について説明する。この場合、概略的には、自車位置推定部１７は、地物情報に含まれる向き情報を参照し、車両の進行方向に対するランドマークの相対的な向きに基づき、車両の進行方向及び横方向の信頼度情報をそれぞれ算出する。以後では、時刻ｔにおける車両の進行方向の信頼度情報を「ａＸ（ｔ）」、横方向の信頼度情報を「ａＹ（ｔ）」とする。【００６４】
ランドマークが車両の進行方向と同一の方向を向いている場合、ランドマークの被照射面が車両の進行方向に対して垂直となる。この場合、計測された横幅ＷＬ（ｔ）及び縦幅ＨＬ（ｔ）と地図上の横幅ＷＭ及び縦幅ＨＭとの差異が大きいと、計測値Ｚ（ｔ）の中心点はずれるものの、車両の進行方向成分の計測値ＬＸ（ｔ）には大きな誤差は生じない。よって、信頼度情報ａＸ（ｔ）の値をあまり小さくする必要はない。一方、車両の横方向成分の計測値ＬＹ（ｔ）には誤差が生じることから、信頼度情報ａＹ（ｔ）の値を小さくする必要がある。【００６５】
以上を勘案し、自車位置推定部１７は、信頼度情報ａＸ（ｔ）と信頼度情報ａＹ（ｔ）とを算出する場合とで異なる感度係数「ｋＸ」、「ｋＹ」をそれぞれ用い、式（２）の感度係数ｋと置き換えて計算を行う。具体的には、自車位置推定部１７は、感度係数ｋＸを用いた以下の式（３）により信頼度情報ａＸ（ｔ）を算出し、感度係数ｋＹを用いた以下の式（４）によりａＹ（ｔ）を算出する。【００６６】
【数３】
2020085062.tif 000005
【００６７】
【数４】
2020085062.tif 000006
そして、自車位置推定部１７は、これらの感度係数ｋＸ、ｋＹを、車両の進行方向に対するランドマークの相対的な向きに基づいて決定する。具体的には、自車位置推定部１７は、ランドマークの向きと車両の進行方向とが平行に近づくほど、信頼度情報ａＸ（ｔ）に対する感度係数ｋＸを小さくし、信頼度情報ａＹ（ｔ）に対する感度係数ｋＹを大きくする。一方、自車位置推定部１７は、ランドマークの向きと車両の進行方向とが垂直に近づくほど、信頼度情報ａＸ（ｔ）に対する感度係数ｋＸを大きくし、信頼度情報ａＹ（ｔ）に対する感度係数ｋＹを小さくする。【００６８】
図１０は、２次元のワールド座標系（絶対座標系）及び車両座標系により表された車両とランドマークとの位置関係を示した図である。ここで、ワールド座標系は、所定地点を原点とし、互いに垂直な座標軸「ｘｗ」及び座標軸「ｙｗ」を有し、車両座標系は、車両の中心を原点とし、車両の進行方向に沿った座標軸「ｘｂ」と車両の側面方向に沿った座標軸「ｙｂ」を有する。そして、図１０では、ワールド座標系における車両のヨー角を「ψ−」、車両の位置を［ｘ−、ｙ−］Ｔとしている。また、ワールド座標系におけるランドマークのヨー角を「Ｍψ」、位置を［Ｍｘ、Ｍｙ］Ｔとし、車両座標系におけるランドマークのヨー角を「Ｌψ」、位置を［Ｌｘ、Ｌｙ］Ｔとしている。【００６９】
この場合、車両から見たランドマークの向きは、車両座標系におけるランドマークのヨー角Ｌψに等しく、ヨー角Ｌψは、図１０に示すように、「Ｍψ−ψ−」に等しい。従って、自車位置推定部１７は、感度係数ｋＸ、ｋＹを、「Ｍψ−ψ−」と定数「ｃ」とを用いてそれぞれ以下の式（５）及び式（６）により定める。【００７０】
【数５】
2020085062.tif 000007
【００７１】
【数６】
2020085062.tif 000008
上記の式（５）及び式（６）によれば、車両から見てランドマークが正面を向いているほど、感度係数ｋＸは最小値ｃに近づき、感度係数ｋＹは大きな値となる。言い換えると、車両から見てランドマークが横を向いているほど、感度係数ｋＸは大きくなり、感度係数ｋＹは最小値ｃに近づく。【００７２】
図１１（Ａ）は、「Ｍψ−ψ−」が１０°の場合において、横幅ＷＭが６０ｃｍ、縦幅ＨＭが６０ｃｍとなるランドマークに対して縦幅ＨＬ（ｔ）を６０ｃｍに固定したときの横幅ＷＬ（ｔ）と式（３）に基づく信頼度情報ａＸ（ｔ）との関係を示す２次元グラフである。また、図１１（Ｂ）は、図１１（Ａ）と同一条件下での横幅ＷＬ（ｔ）と式（４）に基づく信頼度情報ａＹ（ｔ）との関係を示す２次元グラフである。なお、図１１及び後述の図１２の各グラフでは、定数ｃは２０に設定されている。【００７３】
図１１（Ａ）の例では、「Ｍψ−ψ−」が１０°であることから、感度係数ｋＸは、式（５）に基づき「１．０１５ｃ」となる。一方、図１１（Ｂ）の例では、感度係数ｋＹは、式（６）に基づき「５．７５９ｃ」となる。このように、車両から見てランドマークがほぼ正面を向いているときには、感度係数ｋＸはｃに近い値となり、感度係数ｋＹは感度係数ｋＸよりも遥かに大きい値となる。そして、この場合、信頼度情報ａＹ（ｔ）は、計測された横幅ＷＬ（ｔ）及び縦幅ＨＬ（ｔ）と地図上の横幅ＷＭ及び縦幅ＨＭとのずれに対する感度が高くなる。【００７４】
図１２（Ａ）は、「Ｍψ−ψ−」が６０°の場合において、横幅ＷＭが６０ｃｍ、縦幅ＨＭが６０ｃｍとなるランドマークに対して縦幅ＨＬ（ｔ）を６０ｃｍに固定したときの横幅ＷＬ（ｔ）と式（３）に基づく信頼度情報ａＸ（ｔ）との関係を示す２次元グラフである。また、図１２（Ｂ）は、図１２（Ａ）と同一条件下での横幅ＷＬ（ｔ）と式（４）に基づく信頼度情報ａＹ（ｔ）との関係を示す２次元グラフである。【００７５】
図１２（Ａ）の例では、「Ｍψ−ψ−」が６０°であることから、感度係数ｋＸは、式（５）に基づき「２．０ｃ」となる。一方、図１２（Ｂ）の例では、感度係数ｋＹは、式（６）に基づき「１．１５５ｃ」となる。このように、車両から見てランドマークがやや横向きに近いときには、感度係数ｋＸは感度係数ｋＹよりも大きい値となる。そして、この場合、信頼度情報ａＸ（ｔ）は、計測された横幅ＷＬ（ｔ）及び縦幅ＨＬ（ｔ）と地図上の横幅ＷＭ及び縦幅ＨＭとのずれに対する感度が信頼度情報ａＹ（ｔ）と比べて相対的に高くなる。【００７６】
［信頼度情報に応じたカルマンゲイン］
次に、信頼度情報に応じてカルマンゲインを設定する方法について説明する。【００７７】
自車位置推定部１７は、以下の一般式（７）によりカルマンゲインＫ（ｔ）を算出する際に用いる観測雑音行列「Ｒ（ｔ）」の対角成分に、信頼度情報ａＸ（ｔ）、ａＹ（ｔ）の逆数をそれぞれ乗じる。【００７８】
【数７】
2020085062.tif 000009
ここで、対角成分に信頼度情報ａＸ（ｔ）、ａＹ（ｔ）の逆数が乗じられた後の観測雑音行列Ｒ（ｔ）（「修正観測雑音行列Ｒ（ｔ）’」とも呼ぶ。）は、以下の式（８）により表される。【００７９】
【数８】
2020085062.tif 000010
このように、式（３）及び式（４）に基づき信頼度情報ａＸ（ｔ）、ａＹ（ｔ）をそれぞれ生成した場合、観測雑音行列の各成分に乗じる係数は個別となる。なお、自車位置推定部１７は、式（２）に基づき信頼度情報ａ（ｔ）を算出した場合には、式（８）において「ａＸ（ｔ）＝ａＹ（ｔ）＝ａ（ｔ）」とみなして修正観測雑音行列Ｒ（ｔ）’を算出すればよい。【００８０】
式（７）によれば、信頼度情報が最高値の「１」であった場合には、その逆数も１であるため、修正観測雑音行列Ｒ（ｔ）’が元の観測雑音行列Ｒ（ｔ）のままとなる。一方、信頼度情報が「０」に近い値であった場合には、その逆数は１より大きい値となり、修正観測雑音行列Ｒ（ｔ）’の対角成分は増倍されることになる。【００８１】
そして、自車位置推定部１７は、式（７）に示す修正観測雑音行列Ｒ（ｔ）’を用いて、以下の式（９）から適応的なカルマンゲイン「Ｋ（ｔ）’」を算出する。【００８２】
【数９】
2020085062.tif 000011
そして、自車位置推定部１７は、式（９）のカルマンゲインＫ（ｔ）’を用いて、式（１）に基づき予測自車位置Ｘ−（ｔ）を補正した推定自車位置Ｘ＾（ｔ）を算出する。従って、信頼度情報が低い場合は、カルマンゲインＫ（ｔ）’が小さくなるため、予測自車位置Ｘ−（ｔ）に対する補正量が少なくなる。この場合、不正確な補正を防止できるため、推定自車位置Ｘ＾（ｔ）の精度が好適に向上する。【００８３】
なお、自車位置推定部１７は、観測雑音行列Ｒ（ｔ）の対角成分に、信頼度情報ａＸ（ｔ）、ａＹ（ｔ）の逆数をそれぞれ乗じる代わりに、カルマンゲインＫ（ｔ）に対して信頼度情報ａＸ（ｔ）、ａＹ（ｔ）の逆数をそれぞれ乗じてもよい。この場合、自車位置推定部１７は、カルマンゲインＫ（ｔ）’を以下の式（１０）に基づき算出する。【００８４】
【数１０】
2020085062.tif 000012
この式（１０）を用いた場合であっても、信頼度情報が低い場合は、カルマンゲインＫ（ｔ）’が小さくなるため、予測自車位置Ｘ−（ｔ）に対する補正量が少なくなり、不正確な補正を好適に防止することができる。【００８５】
［処理フロー］
図１３は、車載機１の自車位置推定部１７により行われる自車位置推定処理のフローチャートである。車載機１は、図１３のフローチャートの処理を繰り返し実行する。ここでは、一例として、車両の進行方向及び横方向のそれぞれに対する信頼度情報ａＸ（ｔ）、ａＹ（ｔ）を算出する例について説明する。【００８６】
まず、自車位置推定部１７は、ＧＰＳ受信機５等の出力に基づき、自車位置の初期値を設定する（ステップＳ１０１）。次に、自車位置推定部１７は、車速センサ４から車体速度を取得すると共に、ジャイロセンサ３からヨー方向の角速度を取得する（ステップＳ１０２）。そして、自車位置推定部１７は、ステップＳ１０２の取得結果に基づき、車両の移動距離と車両の方位変化を計算する（ステップＳ１０３）。【００８７】
その後、自車位置推定部１７は、１時刻前の推定自車位置Ｘ＾（ｔ−１）に、ステップＳ１０３で計算した移動距離と方位変化を加算し、予測自車位置Ｘ−（ｔ）を算出する（ステップＳ１０４）。さらに、自車位置推定部１７は、予測自車位置Ｘ−（ｔ）に基づき、地図ＤＢ１０の地物情報を参照し、ライダ２の計測範囲内となるランドマークを探索する（ステップＳ１０５）。【００８８】
そして、自車位置推定部１７は、予測自車位置Ｘ−（ｔ）及びステップＳ１０５で探索したランドマークの地物情報が示す位置座標から、計測予測値Ｚ−（ｔ）を算出する（ステップＳ１０６）。さらに、ステップＳ１０６では、自車位置推定部１７は、ステップＳ１０５で探索したランドマークに対するライダ２の計測データから計測値Ｚ（ｔ）を算出する。【００８９】
そして、自車位置推定部１７は、予測自車位置Ｘ−（ｔ）が示す自車の方位（ヨー角）ψ−に相当する車両の進行方向と、地図ＤＢ１０に記録されたランドマークの向き情報とに基づき、車両の進行方向に対する感度係数ｋＸと、車両の横方向に対する感度係数ｋＹとをそれぞれ決定する（ステップＳ１０７）。この場合、例えば、自車位置推定部１７は、式（５）に基づき感度係数ｋＸを算出し、式（６）に基づき感度係数ｋＹを算出する。【００９０】
次に、自車位置推定部１７は、ランドマークの地図上の横幅ＷＭ及び縦幅ＨＭと計測した横幅ＷＬ（ｔ）及び縦幅ＨＬ（ｔ）との各差分と、感度係数ｋＸ、ｋＹとに基づき、式（３）及び式（４）を参照して、信頼度情報ａＸ（ｔ）、ａＹ（ｔ）をそれぞれ算出する（ステップＳ１０８）。そして、自車位置推定部１７は、信頼度情報ａＸ（ｔ）、ａＹ（ｔ）に応じてカルマンゲインＫ（ｔ）’を生成する（ステップＳ１０９）。例えば、自車位置推定部１７は、信頼度情報ａＸ（ｔ）、ａＹ（ｔ）を用いて式（８）に基づき修正観測雑音行列Ｒ（ｔ）’を算出後、式（９）に基づきカルマンゲインＫ（ｔ）’を算出する。その後、自車位置推定部１７は、カルマンゲインＫ（ｔ）´を式（１）のＫ（ｔ）の代わりに用いることで予測自車位置Ｘ−（ｔ）を補正し、推定自車位置Ｘ＾（ｔ）を算出する（ステップＳ１１０）。【００９１】
以上説明したように、本実施例に係る車載機１の自車位置推定部１７は、ライダ２によるランドマークの計測結果である点群データを取得すると共に、地図ＤＢ１０に含まれるランドマークの地物情報を取得する。そして、自車位置推定部１７は、取得した点群データに基づき特定したランドマークのサイズと地物情報が示すランドマークのサイズとの差分に基づいて、ライダ２によるランドマークの計測精度を示す信頼度情報を算出する。これにより、自車位置推定部１７は、推定自車位置Ｘ＾（ｔ）を算出する際の予測自車位置Ｘ−（ｔ）に対する補正量を、ライダ２の計測結果の正確さに応じて的確に定めることができる。【００９２】
［変形例］
以下、実施例に好適な変形例について説明する。以下の変形例は、組み合わせてこれらの実施例に適用してもよい。【００９３】
（変形例１）
自車位置推定部１７は、式（２）に代えて、計測された横幅及び縦幅と地図上の横幅及び縦幅とのそれぞれの差の絶対値｜ＷＬ（ｔ）−ＷＭ｜、｜ＨＬ（ｔ）−ＨＭ｜を用いた以下の式（１１）に基づき信頼度情報ａ（ｔ）を算出してもよい。【００９４】
【数１１】
2020085062.tif 000013
式（１１）に基づく信頼度情報ａ（ｔ）は、式（２）に基づく信頼度情報ａ（ｔ）と同様、０から１までの値域となるように正規化され、計測された横幅ＷＬ（ｔ）及び縦幅ＨＬ（ｔ）が地図上の横幅ＷＭ及び縦幅ＨＭと近いほど最大値の１に近づき、横幅ＷＬ（ｔ）及び縦幅ＨＬ（ｔ）が横幅ＷＭ及び縦幅ＨＭから離れるほど０に近づく。【００９５】
図１４（Ａ）は、横幅ＷＭが６０ｃｍ、縦幅ＨＭが６０ｃｍとなるランドマークに対して計測した縦幅ＨＬ（ｔ）を６０ｃｍに固定した場合の横幅ＷＬ（ｔ）と式（１１）に基づく信頼度情報ａ（ｔ）との関係を示す２次元グラフである。また、図１４（Ｂ）は、図１４（Ａ）と同一のランドマークに対して計測した横幅ＷＬ（ｔ）を６０ｃｍに固定した場合の縦幅ＨＬ（ｔ）と式（１１）に基づく信頼度情報ａ（ｔ）との関係を示す２次元グラフである。さらに、図１４（Ｃ）は、図１４（Ａ）、（Ｂ）と同一のランドマークに対する横幅ＷＬ（ｔ）と縦幅ＨＬ（ｔ）と式（１１）に基づく信頼度情報ａ（ｔ）との関係を示す３次元グラフである。また、図１４（Ｄ）は、図１４（Ａ）における信頼度情報ａ（ｔ）の逆数を示したグラフである。【００９６】
図１４（Ａ）及び図１４（Ｃ）に示すように、横幅ＷＬ（ｔ）は、横幅ＷＭである６０ｃｍに近づくほど信頼度情報ａ（ｔ）が１に近づく。一方、図１４（Ｂ）及び図１４（Ｃ）に示すように、縦幅ＨＬ（ｔ）は、縦幅ＨＭである６０ｃｍに近づくほど信頼度情報ａ（ｔ）が１に近づく。また、図１４（Ｄ）に示すように、信頼度情報ａ（ｔ）の逆数は、横幅ＷＬ（ｔ）が横幅ＷＭに近づくほど１に近づき、横幅ＷＬ（ｔ）が横幅ＷＭから遠ざかるほど大きい値となる。同様に、信頼度情報ａ（ｔ）の逆数は、縦幅ＨＬ（ｔ）が縦幅ＨＭに近づくほど１に近づき、縦幅ＨＬ（ｔ）が縦幅ＨＭから遠ざかるほど大きい値となる。【００９７】
このように、式（１１）に基づく信頼度情報ａ（ｔ）は、式（２）に基づく信頼度情報ａ（ｔ）と同様、最大値が１となり、計測された横幅及び縦幅と地図上の横幅及び縦幅とのそれぞれの差が大きいほど小さい値となる。よって、自車位置推定部１７は、式（１１）に基づく信頼度情報ａ（ｔ）の逆数を用いた式（８）に基づき修正観測雑音行列Ｒ（ｔ）’を算出することで、推定自車位置Ｘ＾（ｔ）を算出する際の予測自車位置Ｘ−（ｔ）に対する補正量を好適に定めることができる。【００９８】
なお、自車位置推定部１７は、実施例と同様、車両の進行方向及び横方向のそれぞれに対する信頼度情報ａＸ（ｔ）及び信頼度情報ａＹ（ｔ）を算出してもよい。この場合、自車位置推定部１７は、感度係数ｋＸ、ｋＹを式（５）、（６）に基づき算出することで、信頼度情報ａＸ（ｔ）及び信頼度情報ａＹ（ｔ）をそれぞれ算出する。【００９９】
（変形例２）
自車位置推定部１７は、式（２）に代えて、信頼度が高いほど最小値の１に近づき、信頼度が低いほど大きな値になるように信頼度情報ａ（ｔ）を定めてもよい。【０１００】
例えば、自車位置推定部１７は、以下の式（１２）に基づき、信頼度情報ａ（ｔ）を算出する。【０１０１】
【数１２】
2020085062.tif 000014
式（１２）に基づく信頼度情報ａ（ｔ）は、１以上の値域となるように正規化され、計測された横幅ＷＬ（ｔ）及び縦幅ＨＬ（ｔ）が地図上の横幅ＷＭ及び縦幅ＨＭと近いほど最小値の１に近づき、横幅ＷＬ（ｔ）及び縦幅ＨＬ（ｔ）が横幅ＷＭ及び縦幅ＨＭから離れるほど大きい値となる。【０１０２】
図１５（Ａ）は、横幅ＷＭが６０ｃｍ、縦幅ＨＭが６０ｃｍとなるランドマークに対して計測した縦幅ＨＬ（ｔ）を６０ｃｍに固定した場合の横幅ＷＬ（ｔ）と式（１２）に基づく信頼度情報ａ（ｔ）との関係を示す２次元グラフである。また、図１５（Ｂ）は、図１５（Ａ）と同一のランドマークに対して計測した横幅ＷＬ（ｔ）を６０ｃｍに固定した場合の縦幅ＨＬ（ｔ）と式（１２）に基づく信頼度情報ａ（ｔ）との関係を示す２次元グラフである。さらに、図１５（Ｃ）は、図１５（Ａ）、（Ｂ）と同一のランドマークに対する横幅ＷＬ（ｔ）と縦幅ＨＬ（ｔ）と式（１２）に基づく信頼度情報ａ（ｔ）との関係を示す３次元グラフである。【０１０３】
図１５（Ａ）及び図１５（Ｃ）に示すように、信頼度情報ａ（ｔ）は、横幅ＷＭである６０ｃｍに横幅ＷＬ（ｔ）が近づくほど最小値である１に近づき、横幅ＷＭから横幅ＷＬ（ｔ）が遠ざかるほど大きい値となる。同様に、図１５（Ｂ）及び図１５（Ｃ）に示すように、信頼度情報ａ（ｔ）は、縦幅ＨＭである６０ｃｍに縦幅ＨＬ（ｔ）が近づくほど最小値である１に近づき、縦幅ＨＭから縦幅ＨＬ（ｔ）が遠ざかるほど大きい値となる。【０１０４】
なお、自車位置推定部１７は、実施例と同様、車両の進行方向及び横方向のそれぞれに対する信頼度情報ａＸ（ｔ）及び信頼度情報ａＹ（ｔ）を算出してもよい。この場合、自車位置推定部１７は、感度係数ｋＸ、ｋＹを式（５）、（６）に基づき算出することで、信頼度情報ａＸ（ｔ）及び信頼度情報ａＹ（ｔ）をそれぞれ算出する。【０１０５】
また、式（１２）に基づく信頼度情報ａ（ｔ）は、式（２）に基づく信頼度情報ａ（ｔ）の逆数（図７（Ｄ）参照）と同様に、最小値が１となり、かつ、横幅ＷＭから横幅ＷＬ（ｔ）が遠ざかる又は縦幅ＨＭから縦幅ＨＬ（ｔ）が遠ざかるほど大きい値となる。よって、自車位置推定部１７は、修正観測雑音行列Ｒ（ｔ）’を算出する際に、信頼度情報ａ（ｔ）を逆数にして観測雑音行列Ｒ（ｔ）の対角成分に乗じる必要がない。具体的には、自車位置推定部１７は、信頼度情報ａＸ（ｔ）及び信頼度情報ａＹ（ｔ）をそれぞれ算出した場合には、以下の式（１３）に基づき、修正観測雑音行列Ｒ（ｔ）’を算出すればよい。【０１０６】
【数１３】
2020085062.tif 000015
これにより、実施例と同様、自車位置推定部１７は、推定自車位置Ｘ＾（ｔ）を算出する際の予測自車位置Ｘ−（ｔ）に対する補正量を好適に定めることができる。【０１０７】
（変形例３）
図１に示す運転支援システムの構成は一例であり、本発明が適用可能な運転支援システムの構成は図１に示す構成に限定されない。例えば、運転支援システムは、車載機１を有する代わりに、車両の電子制御装置が車載機１の自車位置推定部１７の処理を実行してもよい。この場合、地図ＤＢ１０は、例えば車両内の記憶部に記憶され、車両の電子制御装置は、図１３のフローチャートの処理を実行することで推定自車位置を算出する。【０１０８】
（変形例４）
信頼度情報は、推定自車位置Ｘ＾（ｔ）を算出する際の予測自車位置Ｘ−（ｔ）の補正量を決定するパラメータとして用いられることに限定されない。例えば、車載機１は、算出した信頼度情報を、ライダ２が出力する点群データに基づく障害物検知などの他の用途に用いてもよい。
【符号の説明】
【０１０９】
１  車載機
２  ライダ
３  ジャイロセンサ
４  車速センサ
５  ＧＰＳ受信機
１０  地図ＤＢ
【図１】
2020085062.tif 000016
【図２】
2020085062.tif 000017
【図３】
2020085062.tif 000018
【図４】
2020085062.tif 000019
【図５】
2020085062.tif 000020
【図６】
2020085062.tif 000021
【図７】
2020085062.tif 000022
【図８】
2020085062.tif 000023
【図９】
2020085062.tif 000024
【図１０】
2020085062.tif 000025
【図１１】
2020085062.tif 000026
【図１２】
2020085062.tif 000027
【図１３】
2020085062.tif 000028
【図１４】
2020085062.tif 000029
【図１５】
2020085062.tif 000030
【国際調査報告】
2020085062.tif 000031
2020085062.tif 000032