(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2022182658
(43)【公開日】2022-12-08
(54)【発明の名称】軸ずれ推定装置
(51)【国際特許分類】
G01S 7/40 20060101AFI20221201BHJP
G01S 13/931 20200101ALI20221201BHJP
G01S 7/497 20060101ALI20221201BHJP
G01S 17/931 20200101ALI20221201BHJP
【ＦＩ】
G01S7/40 126
G01S13/931
G01S7/497
G01S17/931
【審査請求】未請求
【請求項の数】7
【出願形態】ＯＬ
(21)【出願番号】P 2021090345
(22)【出願日】2021-05-28
(71)【出願人】
【識別番号】000004260
【氏名又は名称】株式会社デンソー
(74)【代理人】
【識別番号】110000578
【氏名又は名称】名古屋国際弁理士法人
(72)【発明者】
【氏名】鈴木 康弘
【テーマコード（参考）】
5J070
5J084
【Ｆターム（参考）】
5J070AB24
5J070AC01
5J070AC02
5J070AC12
5J070AC13
5J070AE01
5J070AF03
5J070AG03
5J070AG07
5J070AH19
5J070AK40
5J084AA04
5J084AA05
5J084AA10
5J084AB01
5J084AC02
5J084BA03
5J084BA11
5J084BA50
5J084CA31
5J084CA32
5J084EA20
(57)【要約】
【課題】軸ずれ量の推定精度を向上させる。
【解決手段】電子制御装置１０は、自車両に搭載されたレーダ装置２０による検出結果に基づいて、車両の前方で走行する先行車両のリフレクタの距離および観測レイヤを取得する。電子制御装置１０は、リフレクタの距離および垂直角度とリフレクタの存在確率との対応関係が予め設定された尤度モデルと、複数のリフレクタの距離および観測レイヤとを用いて、最尤推定法により、レーダ装置２０の軸ずれ量を算出する。
【選択図】図１
【解決手段】電子制御装置１０は、自車両に搭載されたレーダ装置２０による検出結果に基づいて、車両の前方で走行する先行車両のリフレクタの距離および観測レイヤを取得する。電子制御装置１０は、リフレクタの距離および垂直角度とリフレクタの存在確率との対応関係が予め設定された尤度モデルと、複数のリフレクタの距離および観測レイヤとを用いて、最尤推定法により、レーダ装置２０の軸ずれ量を算出する。
【選択図】図１
【特許請求の範囲】
【請求項1】
車両（ＶＨ）に搭載されたレーダ装置（２０）による検出結果に基づいて、前記車両の前方で走行する先行車両のリフレクタの位置を示すリフレクタ位置情報を取得するように構成されたリフレクタ情報取得部（Ｓ１１０，Ｓ１２０）と、
前記リフレクタの位置と前記リフレクタの存在確率との対応関係が予め設定された尤度モデルと、複数の前記リフレクタ位置情報とを用いて、最尤推定法により、前記レーダ装置の軸ずれ量を算出するように構成された軸ずれ量算出部（Ｓ４０～Ｓ６０）と
を備える軸ずれ推定装置（１０）。
前記リフレクタの位置と前記リフレクタの存在確率との対応関係が予め設定された尤度モデルと、複数の前記リフレクタ位置情報とを用いて、最尤推定法により、前記レーダ装置の軸ずれ量を算出するように構成された軸ずれ量算出部（Ｓ４０～Ｓ６０）と
を備える軸ずれ推定装置（１０）。
【請求項2】
請求項１に記載の軸ずれ推定装置であって、
前記尤度モデルにおける前記リフレクタの位置は、前記レーダ装置と前記リフレクタとの間の距離と、前記レーダ装置がレーダ波を送受信する方向となる光軸（ＬＡ）に対して前記リフレクタが存在する垂直方向の角度である垂直角度とを含む軸ずれ推定装置。
前記尤度モデルにおける前記リフレクタの位置は、前記レーダ装置と前記リフレクタとの間の距離と、前記レーダ装置がレーダ波を送受信する方向となる光軸（ＬＡ）に対して前記リフレクタが存在する垂直方向の角度である垂直角度とを含む軸ずれ推定装置。
【請求項3】
請求項１または請求項２に記載の軸ずれ推定装置であって、
前記リフレクタの高さが予め設定された除外判定高を超えているか否かを判断し、前記リフレクタの高さが前記除外判定高を超えている場合には、前記リフレクタの高さが前記除外判定高を超えている前記リフレクタ位置情報を、前記軸ずれ量算出部による前記軸ずれ量の算出から除外するように構成されたリフレクタ高除外部（Ｓ１３０）を備える軸ずれ推定装置。
前記リフレクタの高さが予め設定された除外判定高を超えているか否かを判断し、前記リフレクタの高さが前記除外判定高を超えている場合には、前記リフレクタの高さが前記除外判定高を超えている前記リフレクタ位置情報を、前記軸ずれ量算出部による前記軸ずれ量の算出から除外するように構成されたリフレクタ高除外部（Ｓ１３０）を備える軸ずれ推定装置。
【請求項4】
請求項１～請求項３の何れか１項に記載の軸ずれ推定装置であって、
複数の前記先行車両のうち、一つの前記先行車両を第１先行車両とし、前記第１先行車両とは異なる前記先行車両を第２先行車両として、前記レーダ装置に対して前記第２先行車両の一部が前記第１先行車両に隠されているオクルージョンが発生しているか否かを判断し、前記オクルージョンが発生していると判断した場合には、前記第２先行車両の前記リフレクタ位置情報を、前記軸ずれ量算出部による前記軸ずれ量の算出から除外するように構成されたオクルージョン除外部（Ｓ１４０）を備える軸ずれ推定装置。
複数の前記先行車両のうち、一つの前記先行車両を第１先行車両とし、前記第１先行車両とは異なる前記先行車両を第２先行車両として、前記レーダ装置に対して前記第２先行車両の一部が前記第１先行車両に隠されているオクルージョンが発生しているか否かを判断し、前記オクルージョンが発生していると判断した場合には、前記第２先行車両の前記リフレクタ位置情報を、前記軸ずれ量算出部による前記軸ずれ量の算出から除外するように構成されたオクルージョン除外部（Ｓ１４０）を備える軸ずれ推定装置。
【請求項5】
請求項１～請求項４の何れか１項に記載の軸ずれ推定装置であって、
前記先行車両が虚像であるか否かを判断し、前記先行車両が前記虚像であると判断した場合には、前記虚像であると判断された前記先行車両の前記リフレクタ位置情報を、前記軸ずれ量算出部による前記軸ずれ量の算出から除外するように構成された虚像除外部（Ｓ１５０）を備える軸ずれ推定装置。
前記先行車両が虚像であるか否かを判断し、前記先行車両が前記虚像であると判断した場合には、前記虚像であると判断された前記先行車両の前記リフレクタ位置情報を、前記軸ずれ量算出部による前記軸ずれ量の算出から除外するように構成された虚像除外部（Ｓ１５０）を備える軸ずれ推定装置。
【請求項6】
請求項１～請求項５の何れか１項に記載の軸ずれ推定装置であって、
前記リフレクタの前記存在確率は、販売されている前記車両の前記リフレクタの設置高さの分布に基づいて設定される軸ずれ推定装置。
前記リフレクタの前記存在確率は、販売されている前記車両の前記リフレクタの設置高さの分布に基づいて設定される軸ずれ推定装置。
【請求項7】
請求項１～請求項６の何れか１項に記載の軸ずれ推定装置であって、
前記リフレクタの前記存在確率は、前記レーダ装置が設置される高さに基づいて設定される軸ずれ推定装置。
前記リフレクタの前記存在確率は、前記レーダ装置が設置される高さに基づいて設定される軸ずれ推定装置。
【発明の詳細な説明】
【技術分野】
【0001】
本開示は、レーダ装置の軸ずれ量を推定する軸ずれ推定装置に関する。
【背景技術】
【0002】
特許文献１には、反射波の強度が最大となる反射点における垂直方向の角度と、メインスキャン角度とのずれ角度を演算し、このずれ角度に基づいてメインスキャン角度を修正するレーダ装置が記載されている。
【先行技術文献】
【特許文献】
【0003】
【特許文献1】特許第４８９０９２８号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
アダプティブクルーズコントロール機能の高性能化および高機能化により、レーダ装置は、より遠方の車両を検出することが求められている。このため、レーダ装置が車両に搭載された後におけるレーダ装置の軸ずれ量を高精度に推定する必要がある。
【0005】
レーダ装置によって追跡している車両において反射強度が最も高い反射点（すなわち、リフレクタ）の垂直角度に基づいて軸ずれ量を推定すると、リフレクタの搭載高さのばらつきのために、軸ずれ量の推定精度が低下してしまうという問題があった。
【0006】
本開示は、軸ずれ量の推定精度を向上させることを目的とする。
【課題を解決するための手段】
【0007】
本開示の一態様は、リフレクタ情報取得部（Ｓ１１０，Ｓ１２０）と、軸ずれ量算出部（Ｓ４０～Ｓ６０）とを備える軸ずれ推定装置（１０）である。
リフレクタ情報取得部は、車両（ＶＨ）に搭載されたレーダ装置（２０）による検出結果に基づいて、車両の前方で走行する先行車両のリフレクタの位置を示すリフレクタ位置情報を取得するように構成される。
リフレクタ情報取得部は、車両（ＶＨ）に搭載されたレーダ装置（２０）による検出結果に基づいて、車両の前方で走行する先行車両のリフレクタの位置を示すリフレクタ位置情報を取得するように構成される。
【0008】
軸ずれ量算出部は、リフレクタの位置とリフレクタの存在確率との対応関係が予め設定された尤度モデルと、複数のリフレクタ位置情報とを用いて、最尤推定法により、レーダ装置の軸ずれ量を算出するように構成される。
【0009】
このように構成された本開示の軸ずれ推定装置は、尤度モデルを用いることによって、先行車両のリフレクタの設置位置のばらつきを考慮して軸ずれ量を算出することができるため、軸ずれ量の推定精度を向上させることができる。
【図面の簡単な説明】
【0010】
【図1】車両用障害物認識装置の構成を示すブロック図である。
【図2】レーダ波の照射範囲を示す図である。
【図3】選択視野の構成を示す図である。
【図4】メインルーチンを示すフローチャートである。
【図5】特徴量算出処理を示すフローチャートである。
【図6】自車両の前方に存在する走行車両の強度分布画像を示す図である。
【図7】リフレクタ高判定の方法を説明する図である。
【図8】オクルージョン判定の方法を説明する図である。
【図9】実像および虚像を示す図である。
【図10】リフレクタ搭載高分布を示す図である。
【図11】尤度モデルの算出方法を説明する図である。
【図12】尤度モデルの一例を示す図である。
【図13】軸ずれ量が０°であると仮定したときの確率を算出する方法を説明する図である。
【図14】軸ずれ量が０．６°であると仮定したときの確率を算出する方法を説明する図である。
【図15】平均値および標準偏差の算出方法を説明する図である。
【発明を実施するための形態】
【0011】
以下に本開示の実施形態を図面とともに説明する。
本実施形態の車両用障害物認識装置１は、図１に示すように、電子制御装置１０と、レーダ装置２０と、センサ部３０とを備える。以下、車両用障害物認識装置１を搭載した車両を自車両という。
本実施形態の車両用障害物認識装置１は、図１に示すように、電子制御装置１０と、レーダ装置２０と、センサ部３０とを備える。以下、車両用障害物認識装置１を搭載した車両を自車両という。
【0012】
電子制御装置１０は、ＣＰＵ１１、ＲＯＭおよびＲＡＭ等を備えたマイクロコンピュータを中心に構成された電子制御装置である。マイクロコンピュータの各種機能は、ＣＰＵ１１が非遷移的実体的記録媒体に格納されたプログラムを実行することにより実現される。この例では、ＲＯＭが、プログラムを格納した非遷移的実体的記録媒体に該当する。また、このプログラムの実行により、プログラムに対応する方法が実行される。なお、ＣＰＵ１１が実行する機能の一部または全部を、一つあるいは複数のＩＣ等によりハードウェア的に構成してもよい。また、電子制御装置１０を構成するマイクロコンピュータの数は１つでも複数でもよい。
【0013】
電子制御装置１０は、ＣＰＵ１１と、ＲＡＭ、ＲＯＭおよびフラッシュメモリ等の半導体メモリ１２（以下、メモリ１２）と、通信部１３とを備える。
通信部１３は、通信線を介して他の車載装置との間でデータ通信可能に接続されており、例えばＣＡＮ通信プロトコルに従ってデータの送受信を行う。ＣＡＮは、Controller Area Networkの略である。ＣＡＮは登録商標である。
通信部１３は、通信線を介して他の車載装置との間でデータ通信可能に接続されており、例えばＣＡＮ通信プロトコルに従ってデータの送受信を行う。ＣＡＮは、Controller Area Networkの略である。ＣＡＮは登録商標である。
【0014】
レーダ装置２０は、図２に示すように、自車両ＶＨの前側に搭載され、予め設定された測定周期が経過する毎に、自車両ＶＨの車幅方向に水平および垂直方向の所定角度範囲内でレーダ波を自車両ＶＨの前方へ向けて走査しながら照射し、反射したレーダ波を検出することによって、物体においてレーダ波を反射した箇所（以下、測距点）までの距離と測距点の角度とを測定する。
【0015】
レーダ装置２０は、上記の測定周期を１フレームとして、１フレーム毎に、測距点の距離および角度を測定する。
なお、レーダ装置２０は、例えば、レーダ波としてミリ波帯の電磁波を使用するいわゆるミリ波レーダであってもよいし、レーダ波としてレーザ光を用いるレーザレーダであってもよいし、レーダ波として音波を用いるソナーであってもよい。
なお、レーダ装置２０は、例えば、レーダ波としてミリ波帯の電磁波を使用するいわゆるミリ波レーダであってもよいし、レーダ波としてレーザ光を用いるレーザレーダであってもよいし、レーダ波として音波を用いるソナーであってもよい。
【0016】
レーダ装置２０は、反射したレーダ波をパルス波形で受信する。そしてレーダ装置２０は、受信強度が検出閾値を超えたパルス波形を測距点として検出し、このパルス波形のピーク時の受信強度を反射強度とする。
【0017】
またレーダ装置２０は、パルス波形のピーク時の時間ｔｐに基づいて、測距点の距離を算出する。さらにレーダ装置２０は、レーダ波の走査方向に基づいて、測距点が存在する水平角度および垂直角度を算出する。
【0018】
そしてレーダ装置２０は、測距点の距離、水平角度および垂直角度と、反射強度とを示す測距点情報を電子制御装置１０へ出力する。
また電子制御装置１０は、図１に示すように、通信部１３を介して、レーダ装置２０による測定結果を、例えば、運転支援を行う運転支援装置４０へ送信する。
また電子制御装置１０は、図１に示すように、通信部１３を介して、レーダ装置２０による測定結果を、例えば、運転支援を行う運転支援装置４０へ送信する。
【0019】
センサ部３０は、自車両の挙動を検出する少なくとも１つのセンサを備える。本実施形態では、一例として、センサ部３０は、車速に応じた車速検出信号を電子制御装置１０へ出力する車速センサ３１と、ヨーレートに応じたヨーレート検出信号を電子制御装置１０へ出力するヨーレートセンサ３２とを備える。
【0020】
レーダ装置２０は、水平角度および垂直角度に基づいて、測距点を取得する。そのため、図３の画像Ｇ１で示すように、測距点の位置は、複数の画素を二次元行列状に配列することで表現することができる。１画素は、測距点の距離の情報と、強度情報とを含む。なお、１画素は、反射した複数のレーダ波を受信する場合には、複数の測距点を含む。
【0021】
図３の画像Ｇ１における水平方向軸と垂直方向軸との交点Ｏは、自車両に搭載されたレーダ装置２０がレーダ波を送受信する方向となる光軸ＬＡの延長線上に位置する測距点に対応する。図３は、光軸ＬＡが水平方向ＨＤに対して上方向にずれている状態を示している。光軸ＬＡと水平方向ＨＤとが一致している状態が、軸ずれが無い状態である。
【0022】
電子制御装置１０は、レーダ装置２０により検出可能な視野の一部分を切り出した領域を選択視野とし、選択視野に映っている物体を検出する。
図３に示す画像Ｇ１は、水平方向に沿って２６画素が配列され、且つ、垂直方向に沿って１６画素が配列されることによって形成されている。図３に示す選択視野ＦＶ１，ＦＶ２は、水平方向に沿って２２画素が配列され、且つ、垂直方向に沿って１０画素が配列されることによって形成されている。
図３に示す画像Ｇ１は、水平方向に沿って２６画素が配列され、且つ、垂直方向に沿って１６画素が配列されることによって形成されている。図３に示す選択視野ＦＶ１，ＦＶ２は、水平方向に沿って２２画素が配列され、且つ、垂直方向に沿って１０画素が配列されることによって形成されている。
【0023】
選択視野ＦＶ１は、中心線ＬＣ１を中心として選択された状態で、その中心線ＬＣ１は光軸ＬＡと重なる。一方、選択視野ＦＶ２は、中心線ＬＣ２を中心として選択された状態で、選択視野ＦＶ１よりも下側の領域が選択される。
【0024】
図３は、選択視野ＦＶ１が水平方向ＨＤに対して上方向にずれている状態を示す。この場合に、選択視野ＦＶ２で示すように、選択視野ＦＶ１より下を選択することによって、水平方向ＨＤに対する光軸ＬＡのずれを補正することができる。
【0025】
すなわち、電子制御装置１０は、軸ずれ量に応じて選択視野の位置を変更することによって、軸ずれを補正する。
以下、選択視野において、上から１，２，３，４，５，６，７，８，９，１０番目の行をそれぞれ、第１，２，３，４，５，６，７，８，９，１０観測レイヤＬＹ１，ＬＹ２，ＬＹ３，ＬＹ４，ＬＹ５，ＬＹ６，ＬＹ７，ＬＹ８，ＬＹ９，ＬＹ１０という。図３に記載されている第１～１０観測レイヤＬＹ１～ＬＹ１０は、選択視野ＦＶ１における観測レイヤを示している。
以下、選択視野において、上から１，２，３，４，５，６，７，８，９，１０番目の行をそれぞれ、第１，２，３，４，５，６，７，８，９，１０観測レイヤＬＹ１，ＬＹ２，ＬＹ３，ＬＹ４，ＬＹ５，ＬＹ６，ＬＹ７，ＬＹ８，ＬＹ９，ＬＹ１０という。図３に記載されている第１～１０観測レイヤＬＹ１～ＬＹ１０は、選択視野ＦＶ１における観測レイヤを示している。
【0026】
次に、電子制御装置１０が実行するメインルーチンの手順を説明する。メインルーチンは、電子制御装置１０の動作中において測定周期毎に繰り返し実行される処理である。
メインルーチンが実行されると、電子制御装置１０のＣＰＵ１１は、図４に示すように、まずＳ１０にて、物体追跡を実行する。
メインルーチンが実行されると、電子制御装置１０のＣＰＵ１１は、図４に示すように、まずＳ１０にて、物体追跡を実行する。
【0027】
具体的には、ＣＰＵ１１は、まず、直近のフレーム（すなわち、今回のフレーム）でレーダ装置２０により検出された１または複数の測距点について、測距点の距離および角度に基づいて、測距点の横位置および縦位置を算出する。横位置は、自車両を起点として自車両の車幅方向に沿った位置である。縦位置は、自車両を起点として車幅方向に対して垂直な方向に沿った位置である。
【0028】
そしてＣＰＵ１１は、今回のフレームの測距点（以下、今回測距点）毎に、前回のフレームの測距点（以下、前回測距点）と同一の物標を表すものであるか否かを判定する履歴追尾処理を実行する。
【0029】
具体的には、ＣＰＵ１１は、前回測距点の情報に基づいて、前回測距点に対応する今回測距点の予測位置を算出し、その予測位置と、今回測距点の位置との差分が予め設定された上限値より小さい場合には、履歴接続があるものと判断し、複数のフレーム（例えば５フレーム）に渡って履歴接続があると判断された測距点を物標であると認識する。さらにＣＰＵ１１は、今回測距点の位置と前回測距点の位置との差分と、測定周期（すなわち、１フレームの時間）とに基づいて、今回測距点の相対速度を算出する。
【0030】
次にＣＰＵ１１は、物標であると認識された測距点のうち、同一の物体に起因した測距点を選択するために予め設定された同一物体選択条件を満たす測距点（以下、同一物体測距点）を選択する。なお、本実施形態における同一物体選択条件は、例えば「自車両から最も近い位置に存在している物標の測距点」を代表測距点として、この代表測距点との間で、距離の差が予め設定された距離選択判定値以下あり、且つ、角度の差が予め設定された角度選択判定値以下あり、且つ、相対速度の差が予め設定された相対速度選択判定値以下であることである。
【0031】
さらにＣＰＵ１１は、選択された同一物体測距点の中から、算出された横位置に基づいて、最も右側に位置する測距点（以下、最右端測距点）と、最も左側に位置する測距点（以下、最左端測距点）とを抽出する。そしてＣＰＵ１１は、最右端測距点の横位置と最左端測距点の横位置との中心値を、物体の中心横位置ｘとして決定する。また電子制御装置１０は、最右端測距点の横位置と最左端測距点の横位置との差を、物体の横幅として決定する。
【0032】
さらにＣＰＵ１１は、算出された縦位置に基づいて、最も前側に位置する測距点（以下、最前端測距点）と、最も後側に位置する測距点（以下、最後端測距点）とを抽出する。そしてＣＰＵ１１は、最前端測距点の縦位置と最後端測距点の縦位置との中心値を、物体の中心縦位置ｙとして決定する。またＣＰＵ１１は、最前端測距点の縦位置と最後端測距点の縦位置との差を、物体の奥行として決定する。
【0033】
従って、ＣＰＵ１１は、最右端測距点、最左端測距点、最前端測距点および最後端測距点を囲む矩形を物体として認識する。
さらにＣＰＵ１１は、履歴追尾処理による判定結果に基づいて、前回のフレームにおける物体と今回のフレームにおける物体とで履歴接続がある場合に、前回および今回のフレームにおける物体の中心横位置ｘおよび中心縦位置ｙと、測定周期（すなわち、１フレームの時間）とを用いて、今回のフレームにおける物体の横相対速度ｖｘおよび縦相対速度ｖｙを算出する。
さらにＣＰＵ１１は、履歴追尾処理による判定結果に基づいて、前回のフレームにおける物体と今回のフレームにおける物体とで履歴接続がある場合に、前回および今回のフレームにおける物体の中心横位置ｘおよび中心縦位置ｙと、測定周期（すなわち、１フレームの時間）とを用いて、今回のフレームにおける物体の横相対速度ｖｘおよび縦相対速度ｖｙを算出する。
【0034】
またＣＰＵ１１は、前回のフレームで認識された物体について、今回のフレームで認識された物体と履歴接続できなかった場合には、前回のフレームで認識された物体の中心横位置ｘ、中心縦位置ｙ、横相対速度ｖｘおよび縦相対速度ｖｙに基づいて外挿することにより、今回のフレームで認識された物体の中心横位置ｘ、中心縦位置ｙ、横相対速度ｖｘおよび縦相対速度ｖｙとする。
【0035】
そしてＳ１０の処理が終了すると、ＣＰＵ１１は、Ｓ２０にて、特徴量算出処理を実行する。
ここで、特徴量算出処理の手順を説明する。
ここで、特徴量算出処理の手順を説明する。
【0036】
特徴量算出処理が実行されると、ＣＰＵ１１は、図５に示すように、まずＳ１１０にて、Ｓ１０の処理で認識された物体の中から、予め設定された車両抽出条件に基づいて、車両を抽出する。車両抽出条件は、中心横位置ｘが－２．５ｍ～２．５ｍの範囲内であり、且つ、中心縦位置ｙが３０ｍ～１５０ｍの範囲内であり、且つ、縦絶対速度が４０ｋｍ／ｈ以下となる物体である。縦絶対速度は、縦相対速度ｖｙに自車速を加算した速度である。
【0037】
次にＣＰＵ１１は、Ｓ１２０にて、Ｓ１１０で抽出された車両におけるリフレクタの距離および観測レイヤを特定する。
図６に示す画像Ｇ２は、自車両の前方に存在する走行車両をカメラで撮影することにより得られた画像である。画像Ｇ２において、円ＣＬ１，ＣＬ２で囲まれた領域には、リフレクタが映っている。
図６に示す画像Ｇ２は、自車両の前方に存在する走行車両をカメラで撮影することにより得られた画像である。画像Ｇ２において、円ＣＬ１，ＣＬ２で囲まれた領域には、リフレクタが映っている。
【0038】
図６に示す画像Ｇ３は、自車両の前方に存在する走行車両をレーダ装置２０で検出することによって得られた強度分布画像である。円ＣＬ３，ＣＬ４で囲まれた領域は、リフレクタに対応する測距点である。
【0039】
ＣＰＵ１１は、Ｓ１２０にて、まず、Ｓ１１０で抽出された車両のそれぞれについて、車両を構成する複数の測距点の中から、リフレクタを抽出するために予め設定されたリフレクタ判定値以上の反射強度を有する測距点を抽出する。そしてＣＰＵ１１は、抽出した測距点の中から、最も反射強度が大きい測距点を、リフレクタに対応する測距点として特定する。さらにＣＰＵ１１は、特定した測距点の距離を、リフレクタの距離として設定する。またＣＰＵ１１は、特定した測距点の距離および角度に基づいて、観測レイヤを算出し、この観測レイヤを、リフレクタの観測レイヤとして設定する。
【0040】
Ｓ１２０の処理が終了すると、図５に示すように、ＣＰＵ１１は、Ｓ１３０にて、リフレクタ高判定を実行する。具体的には、ＣＰＵ１１は、まず、選択視野に映っている車両のそれぞれついて、車両を構成する複数の画素の領域（以下、車両画素領域）と、リフレクタが存在する画素（以下、リフレクタ画素）とを特定する。
【0041】
例えば、図７に示す選択視野画像Ｇ４には、２台の車両が映っている。選択視野画像Ｇ４における矩形Ｒ１は第１車両の車両画素領域であり、選択視野画像Ｇ４における矩形Ｒ２は第２車両の車両画素領域である。
【0042】
選択視野画像Ｇ４における矩形Ｒ３は、第１車両において反射強度がリフレクタ判定値以上の画素によって構成される領域である。選択視野画像Ｇ４における矩形Ｒ４は、第４車両において反射強度がリフレクタ判定値以上の画素によって構成される領域である。
【0043】
選択視野画像Ｇ４における矩形Ｒ５は、第１車両のリフレクタ画素である。選択視野画像Ｇ４における矩形Ｒ６は、第２車両のリフレクタ画素である。
そしてＣＰＵ１１は、選択視野に映っている車両のそれぞれついて、車両画素領域における最下端の観測レイヤ（以下、最下端レイヤ）と、リフレクタ画素が位置する観測レイヤ（以下、リフレクタレイヤ）とを特定する。さらにＣＰＵ１１は、最下端観測レイヤとリフレクタレイヤとの間に存在する観測レイヤの数（以下、高さレイヤ数）を算出する。
そしてＣＰＵ１１は、選択視野に映っている車両のそれぞれついて、車両画素領域における最下端の観測レイヤ（以下、最下端レイヤ）と、リフレクタ画素が位置する観測レイヤ（以下、リフレクタレイヤ）とを特定する。さらにＣＰＵ１１は、最下端観測レイヤとリフレクタレイヤとの間に存在する観測レイヤの数（以下、高さレイヤ数）を算出する。
【0044】
図７に示す選択視野画像Ｇ４では、第１車両の最下端レイヤは第９観測レイヤＬＹ９であり、第１車両のリフレクタレイヤは第８観測レイヤＬＹ８である。また、第２車両の最下端レイヤは第９観測レイヤＬＹ９であり、第２車両のリフレクタレイヤは第５観測レイヤＬＹ５である。このため、第１車両の高さレイヤ数は、矢印ＡＬ１で示すように「１」であり、第２車両の高さレイヤ数は、矢印ＡＬ２で示すように「４」である。
【0045】
さらにＣＰＵ１１は、選択視野に映っている車両のそれぞれついて、最下端レイヤに位置する測距点の距離と、高さレイヤ数とを用いて、リフレクタの高さ（以下、リフレクタ高）を算出する。そしてＣＰＵ１１は、選択視野に映っている車両のそれぞれついて、算出されたリフレクタ高が予め設定された除外判定高（例えば、１．５ｍ）を超えるか否かを判断する。そしてＣＰＵ１１は、リフレクタ高が除外判定高を超えている車両を除外し、リフレクタ高が除外判定高以下である車両を推定対象車両に設定する。
【0046】
Ｓ１３０の処理が終了すると、図５に示すように、ＣＰＵ１１は、Ｓ１４０にて、オクルージョン判定を実行する。オクルージョンは、例えば図８に示すように、自車両のレーダ装置２０の前方に第１先行車両ＰＶ１が存在し、さらに第１先行車両ＰＶ１の前方に第２先行車両ＰＶ２が存在している場合において、レーダ装置２０から見て第２先行車両ＰＶ２の一部を第１先行車両ＰＶ１が隠す状態をいう。
【0047】
ＣＰＵ１１は、オクルージョン判定によって、隠されている車両を抽出する。そしてＣＰＵ１１は、推定対象車両に設定されている車両の中から、オクルージョン判定によって抽出した車両を推定対象車両から除外する。
【0048】
以下に、オクルージョン判定の方法の一例を説明する。
図８に示すように、ＣＰＵ１１は、第１先行車両ＰＶ１の存在に起因して矩形状の物体ＲＣ１を認識し、第２先行車両ＰＶ２の存在に起因して矩形状の物体ＲＣ２を認識しているとする。レーダ装置２０から見て物体ＲＣ２が存在する方向と、レーダ装置２０から見て物体ＲＣ１が存在する方向とが一致しており、且つ、物体ＲＣ２の距離は物体ＲＣ１の距離より長い。このため、ＣＰＵ１１は、第２先行車両ＰＶ２の一部を第１先行車両ＰＶ１が隠す状態であると判断することができる。
図８に示すように、ＣＰＵ１１は、第１先行車両ＰＶ１の存在に起因して矩形状の物体ＲＣ１を認識し、第２先行車両ＰＶ２の存在に起因して矩形状の物体ＲＣ２を認識しているとする。レーダ装置２０から見て物体ＲＣ２が存在する方向と、レーダ装置２０から見て物体ＲＣ１が存在する方向とが一致しており、且つ、物体ＲＣ２の距離は物体ＲＣ１の距離より長い。このため、ＣＰＵ１１は、第２先行車両ＰＶ２の一部を第１先行車両ＰＶ１が隠す状態であると判断することができる。
【0049】
隠された車両である第２先行車両ＰＶ２についてはリフクタＲＦ２も隠されている可能性がある。そして、第２先行車両ＰＶ２のリフクタＲＦ２であると誤認識された測距点を用いて軸ずれ量の推定を行うと、軸ずれ量の推定精度が低下する恐れがある。このため、隠されている車両は推定対象車両から除外される。一方、隠されていない車両である第１先行車両ＰＶ１については、リフクタＲＦ１も隠されていないと判断され、推定対象車両から除外されない。
【0050】
Ｓ１４０の処理が終了すると、ＣＰＵ１１は、図５に示すように、Ｓ１５０にて、虚像判定を実行する。虚像は、図９に示すように、レーダ装置２０から送信されたレーダ波が静止物（図９では、トンネル壁や防音壁）で反射した後に物体で反射し、その後、上記の静止物で反射したレーダ波をレーダ装置２０が受信することで、物体と認識される像である。虚像は、壁や測定対象との位置関係によって、高さ方向にずれて見えることがある。このため、虚像を車両として学習してしまうと、軸ずれ量の推定精度が低下する。なお、図９の上側の図は、レーダ装置２０、実像および虚像をレーダ装置２０の上方から見た図である。図９の下側の図は、レーダ装置２０、実像および虚像をレーダ装置２０の後方から見た図である。
【0051】
具体的には、ＣＰＵ１１は、追跡された車両の中から、予め設定されたペア判定条件に基づいて、虚像および実像のペアを特定し、特定したペアのうち、反射強度が小さい方の車両を、虚像として抽出する。虚像は、実像に対して距離および相対速度のそれぞれの差が予め設定された閾値以下（例えば、距離差が５ｍ以下、且つ、速度差が５ｋｍ／ｈ以下）となる特徴があるため、この特徴に基づいてペア判定条件が設定される。
【0052】
ＣＰＵ１１は、虚像判定によって、虚像であると判定された車両を抽出する。そしてＣＰＵ１１は、推定対象車両に設定されている車両の中から、虚像判定によって抽出した車両を推定対象車両から除外する。
【0053】
Ｓ１５０の処理が終了すると、ＣＰＵ１１は、図５に示すように、Ｓ１６０にて、Ｓ１３０～Ｓ１５０の処理によって除外されずに残っている推定対象車両が存在しているか否かを判断する。ここで、推定対象車両が存在していない場合には、ＣＰＵ１１は、特徴量算出処理を終了する。
【0054】
一方、推定対象車両が存在している場合には、ＣＰＵ１１は、Ｓ１７０にて、推定対象車両におけるリフレクタの距離および観測レイヤを、特徴量として、メモリ１２に設けられている特徴量リスト内に記憶して、特徴量算出処理を終了する。
【0055】
特徴量算出処理が終了すると、ＣＰＵ１１は、図４に示すように、Ｓ３０にて、予め設定された推定可能判定条件が成立したか否かを判断する。本実施形態の推定可能判定条件は、後述する第１判定条件、第２判定条件および第３判定条件の全てが成立することである。
【0056】
第１判定条件は、自車両の車速が第１判定速度（例えば、４０ｋｍ／ｈ）以上であることである。なお、自車両の車速は、車速センサ３１によって検出される。
第２判定条件は、自車両が直進していることである。なお、自車両が直進しているか否かは、自車両が走行している道路の曲率半径が予め設定された第２判定半径（例えば、１５００ｍ）以上であるか否かによって判断される。すなわち、曲率半径が第２判定半径以上である場合に、自車両が直進していると判断される。なお、道路の曲率半径は、車速センサ３１によって検出される車速と、ヨーレートセンサ３２によって検出されるヨーレートとに基づいて算出される。
第２判定条件は、自車両が直進していることである。なお、自車両が直進しているか否かは、自車両が走行している道路の曲率半径が予め設定された第２判定半径（例えば、１５００ｍ）以上であるか否かによって判断される。すなわち、曲率半径が第２判定半径以上である場合に、自車両が直進していると判断される。なお、道路の曲率半径は、車速センサ３１によって検出される車速と、ヨーレートセンサ３２によって検出されるヨーレートとに基づいて算出される。
【0057】
第３判定条件は、特徴量リストに、１台以上の推定対象車両におけるリフレクタの距離および観測レイヤが記憶されていることである。
Ｓ３０にて、推定可能判定条件が成立していない場合には、ＣＰＵ１１は、Ｓ７０に移行する。一方、推定可能判定条件が成立している場合には、ＣＰＵ１１は、Ｓ４０にて、軸ずれ量確率分布を算出する。
Ｓ３０にて、推定可能判定条件が成立していない場合には、ＣＰＵ１１は、Ｓ７０に移行する。一方、推定可能判定条件が成立している場合には、ＣＰＵ１１は、Ｓ４０にて、軸ずれ量確率分布を算出する。
【0058】
ここで、軸ずれ量確率分布について説明する。
車両においてリフレクタが搭載されている高さ（以下、リフレクタ搭載高）は、車両毎にばらついている。しかし、リフレクタ搭載高の範囲は保安基準により規定されている。例えば、日本では、「道路運送車両の保安基準の細目を定める告示」における第２１０条で規定され、米国では、Federal Motor Vehicle Safety Standard 108で規定されている。
車両においてリフレクタが搭載されている高さ（以下、リフレクタ搭載高）は、車両毎にばらついている。しかし、リフレクタ搭載高の範囲は保安基準により規定されている。例えば、日本では、「道路運送車両の保安基準の細目を定める告示」における第２１０条で規定され、米国では、Federal Motor Vehicle Safety Standard 108で規定されている。
【0059】
また、販売されている車両の種類（以下、車種）毎のリフレクタ搭載高と、車種毎の販売台数とに基づいて、例えば図１０に示すように、リフレクタ搭載高の分布を算出することが可能である。図１０に示すリフレクタ搭載高分布ＨＤでは、横軸にリフレクタ搭載高が設定され、縦軸に頻度が設定されている。そして、図１０に示すリフレクタ搭載高分布ＨＤでは、リフレクタ搭載高は２００ｍｍ～１５００ｍｍの範囲で分布している。
【0060】
そして、リフレクタ搭載高分布と、レーダ装置２０が自車両に設置されている高さ（以下、レーダ設置高）と、レーダ装置２０の選択視野における観測レイヤとに基づいて、尤度モデルが事前に作成される。尤度モデルは、軸ずれが無い状態において、先行車両のリフレクタの距離および垂直角度に応じたリフレクタの存在確率を定義する。
【0061】
尤度モデルは、レーダ装置２０からの距離を複数の距離区間に分割し、距離区間毎リフレクタ搭載高分布を配置することにより算出される。
説明を簡略化するために、図１１に示すように、選択視野が、上から順に第１，２，３，４，５，６観測レイヤＬＹ１，ＬＹ２，ＬＹ３，ＬＹ４，ＬＹ５，ＬＹ６を備えているとする。第１～６観測レイヤＬＹ１～ＬＹ６における垂直角度範囲は互いに等しい。
説明を簡略化するために、図１１に示すように、選択視野が、上から順に第１，２，３，４，５，６観測レイヤＬＹ１，ＬＹ２，ＬＹ３，ＬＹ４，ＬＹ５，ＬＹ６を備えているとする。第１～６観測レイヤＬＹ１～ＬＹ６における垂直角度範囲は互いに等しい。
【0062】
例えば、レーダ装置２０から距離Ｄ１離れた位置では、リフレクタ搭載高分布ＨＤ１で示すように、第１～４観測レイヤＬＹ１～ＬＹ４での存在確率が高く、第５観測レイヤＬＹ５での存在確率が２番目に低く、第６観測レイヤＬＹ６での存在確率が最も低い０である。
【0063】
また、レーダ装置２０から距離Ｄ２離れた位置では、リフレクタ搭載高分布ＨＤ２で示すように、第５，６観測レイヤＬＹ５，ＬＹ６での存在確率が最も低い０であり、第４観測レイヤＬＹ４での存在確率が３番目に低く、第１観測レイヤＬＹ１での存在確率が４番目に低く、第３観測レイヤＬＹ３での存在確率が最も高い。
【0064】
また、レーダ装置２０から距離Ｄ３離れた位置では、リフレクタ搭載高分布ＨＤ３で示すように、第１，５，６観測レイヤＬＹ１，ＬＹ５，ＬＹ６での存在確率が最も低い０であり、第２観測レイヤＬＹ２での存在確率が４番目に低く、第４観測レイヤＬＹ４での存在確率が５番目に低く、第３観測レイヤＬＹ３での存在確率が最も高い。
【0065】
図１２に示す尤度モデルＬＭでは、横軸にレーダ装置２０からの距離が設定され、縦軸に垂直角度が設定され、距離と垂直角度との組み合わせに応じて尤度が設定される。また尤度モデルＬＭでは、距離区間が５ｍに設定され、垂直角度区間が０．２°に設定されている。
【0066】
例えば、矩形Ｒ１１で示すように、距離が１０ｍであり、垂直角度が２．３°である場合の尤度は「１０」である。また、矩形Ｒ１２で示すように、距離が３５ｍであり、垂直角度が０．３°である場合の尤度は「３３」である。また、矩形Ｒ１３で示すように、距離が７０ｍであり、垂直角度が１．９°である場合の尤度は「０」である。
【0067】
そしてＣＰＵ１１は、Ｓ４０にて、特徴量リスト内に記憶されている複数の特徴量（すなわち、リフレクタの距離および観測レイヤ）と、尤度モデルとを用いて、軸ずれ量確率分布を算出する。
【0068】
軸ずれ量をｘとし、特徴量をｚとし、特徴量リスト内に記憶されている特徴量の数をｍとして、軸ずれ量確率分布Ｐ（ｚ｜ｘ）は、下式（１）で表される。
【0069】
【数1】
【0070】
ここで、軸ずれ量確率分布Ｐ（ｚ｜ｘ）の算出方法を図１３および図１４を用いて説明する。
例えば、特徴量リスト内に第１特徴量、第２特徴量および第３特徴量が記憶されているとする。そして第１特徴量において、リフレクタの距離が３０ｍであり、観測レイヤが第１観測レイヤＬＹ１であるとする。また第２特徴量において、リフレクタの距離が６０ｍであり、観測レイヤが第２観測レイヤＬＹ２であるとする。また第３特徴量において、リフレクタの距離が８０ｍであり、観測レイヤが第３観測レイヤＬＹ３であるとする。
例えば、特徴量リスト内に第１特徴量、第２特徴量および第３特徴量が記憶されているとする。そして第１特徴量において、リフレクタの距離が３０ｍであり、観測レイヤが第１観測レイヤＬＹ１であるとする。また第２特徴量において、リフレクタの距離が６０ｍであり、観測レイヤが第２観測レイヤＬＹ２であるとする。また第３特徴量において、リフレクタの距離が８０ｍであり、観測レイヤが第３観測レイヤＬＹ３であるとする。
【0071】
図１３に示すように、軸ずれ量が０°であると仮定した場合には、第１観測レイヤＬＹ１は垂直角度が１．２°～１．８°の範囲を含むとする。第２観測レイヤＬＹ２は垂直角度が０．６°～１．２°の範囲を含むとする。第３観測レイヤＬＹ３は垂直角度が０．０°～０．６°の範囲を含むとする。第４観測レイヤＬＹ４は垂直角度が－０．６°～０．０°の範囲を含むとする。第５観測レイヤＬＹ５は垂直角度が－１．２°～－０．６°の範囲を含むとする。第６観測レイヤＬＹ６は垂直角度が－１．８°～－１．２°の範囲を含むとする。
【0072】
第１特徴量では、リフレクタの距離が３０ｍであり、観測レイヤが第１観測レイヤＬＹ１である。このため、尤度モデルＬＭにおいて、リフレクタの距離が３０ｍであり、角度が１．３°，１．５°，１．７°であるときの尤度、すなわち、「７」と「８」と「１０」とが抽出される。そして、「７」と「８」と「１０」との総和である「２５」が、軸ずれ量が０°であると仮定したときの第１特徴量に対する尤度となる。この尤度は、Ｌ（ｚ１｜ｘ＝０°）に相当する。
【0073】
また第２特徴量では、リフレクタの距離が６０ｍであり、観測レイヤが第２観測レイヤＬＹ２である。このため、尤度モデルＬＭにおいて、リフレクタの距離が６０ｍであり、角度が０．７°，０．９°，１．１°であるときの尤度、すなわち、「５」と「８」と「１５」とが抽出される。そして、「５」と「８」と「１５」との総和である「２８」が、軸ずれ量が０°であると仮定したときの第２特徴量に対する尤度となる。この尤度は、Ｌ（ｚ２｜ｘ＝０°）に相当する。
【0074】
また第３特徴量では、リフレクタの距離が８０ｍであり、観測レイヤが第３観測レイヤＬＹ３である。このため、尤度モデルＬＭにおいて、リフレクタの距離が８０ｍであり、角度が０．１°，０．３°，０．５°であるときの尤度、すなわち、「２１」と「３５」と「５０」とが抽出される。そして、「２１」と「３５」と「５０」との総和である「１０６」が、軸ずれ量が０°であると仮定したときの第３特徴量に対する尤度となる。この尤度は、Ｌ（ｚ３｜ｘ＝０°）に相当する。
【0075】
このため、ずれ量が０°であると仮定したときの第１，２，３特徴量に対応する尤度の総和、すなわち、「２５」と「２８」と「１０６」との総和である「１５９」が、軸ずれ量が０°であると仮定したときの確率Ｐ（ｚ｜ｘ＝０°）に相当する値となる。
【0076】
また、図１４に示すように、軸ずれ量が０．６°であると仮定した場合には、第１観測レイヤＬＹ１は垂直角度が１．８°～２．４°の範囲を含む。第２観測レイヤＬＹ２は垂直角度が１．２°～１．８°の範囲を含む。第３観測レイヤＬＹ３は垂直角度が０．６°～１．２°の範囲を含む。第４観測レイヤＬＹ４は垂直角度が０．０°～０．６°の範囲を含む。第５観測レイヤＬＹ５は垂直角度が－０．６°～０．０°の範囲を含む。第６観測レイヤＬＹ６は垂直角度が－１．２°～－０．６°の範囲を含む。
【0077】
第１特徴量では、リフレクタの距離が３０ｍであり、観測レイヤが第１観測レイヤＬＹ１である。このため、尤度モデルＬＭにおいて、リフレクタの距離が３０ｍであり、角度が１．９°，２．１°，２．３°であるときの尤度、すなわち、「５」と「４」と「２」とが抽出される。そして、「５」と「４」と「２」との総和である「１１」が、軸ずれ量が０．６°であると仮定したときの第１特徴量に対する尤度となる。この尤度は、Ｌ（ｚ１｜ｘ＝０．６°）に相当する。
【0078】
また第２特徴量では、リフレクタの距離が６０ｍであり、観測レイヤが第２観測レイヤＬＹ２である。このため、尤度モデルＬＭにおいて、リフレクタの距離が６０ｍであり、角度が１．３°，１．５°，１．７°であるときの尤度、すなわち、「０」と「０」と「０」とが抽出される。そして、「０」と「０」と「０」との総和である「０」が、軸ずれ量が０．６°であると仮定したときの第２特徴量に対する尤度となる。この尤度は、Ｌ（ｚ２｜ｘ＝０．６°）に相当する。
【0079】
また第３特徴量では、リフレクタの距離が８０ｍであり、観測レイヤが第３観測レイヤＬＹ３である。このため、尤度モデルＬＭにおいて、リフレクタの距離が８０ｍであり、角度が０．７°，０．９°，１．１°であるときの尤度、すなわち、「８」と「４」と「０」とが抽出される。そして、「８」と「４」と「０」との総和である「１２」が、軸ずれ量が０．６°であると仮定したときの第３特徴量に対する尤度となる。この尤度は、Ｌ（ｚ３｜ｘ＝０．６°）に相当する。
【0080】
このため、軸ずれ量が０．６°であると仮定したときの第１，２，３特徴量に対応する尤度の総和、すなわち、「１１」と「０」と「１２」との総和である「２３」が、軸ずれ量が０．６°であると仮定したときの確率Ｐ（ｚ｜ｘ＝０．６°）に相当する値となる。
【0081】
そして、軸ずれ量がφであると仮定したときの確率Ｐ（ｚ｜ｘ＝φ）の算出を、例えば、尤度モデルＬＭで設定されている角度－１．７°，－１．５°，－１．３°，－１．１°，・・・，－０．１°，０．１°，０．３°，・・・，３．３°，３．５°，３．７°のそれぞれで行うことによって、軸ずれ量確率分布Ｐ（ｚ｜ｘ）が得られる。
【0082】
但し、軸ずれ量確率分布Ｐ（ｚ｜ｘ）は、ｘの全範囲で積分した値が１となるように正規化される。
Ｓ４０の処理が終了すると、ＣＰＵ１１は、図４に示すように、Ｓ５０にて、軸ずれ量確率分布を更新する。具体的には、ＣＰＵ１１は、前回のフレームにおけるＳ５０で更新された軸ずれ量確率分布をＰｔ－１とし、今回のフレームにおけるＳ４０で算出された軸ずれ量確率分布をＰＯとし、今回のフレームにおけるＳ５０で更新される軸ずれ量確率分布をＰｔとし、重み係数をａとして、下式（２）により軸ずれ量確率分布を更新する。
Ｓ４０の処理が終了すると、ＣＰＵ１１は、図４に示すように、Ｓ５０にて、軸ずれ量確率分布を更新する。具体的には、ＣＰＵ１１は、前回のフレームにおけるＳ５０で更新された軸ずれ量確率分布をＰｔ－１とし、今回のフレームにおけるＳ４０で算出された軸ずれ量確率分布をＰＯとし、今回のフレームにおけるＳ５０で更新される軸ずれ量確率分布をＰｔとし、重み係数をａとして、下式（２）により軸ずれ量確率分布を更新する。
【0083】
Ｐｔ ＝ ａ×Ｐｔ－１＋（１－ａ）×ＰＯ ・・・（２）
さらにＣＰＵ１１は、Ｓ６０にて、Ｓ５０で更新された軸ずれ量確率分布に基づいて軸ずれ量を算出し、Ｓ７０に移行する。具体的には、ＣＰＵ１１は、図１５に示すように、Ｓ５０で更新された軸ずれ量確率分布Ｐｔを正規分布ＮＤと仮定して、平均値および標準偏差を算出する。
さらにＣＰＵ１１は、Ｓ６０にて、Ｓ５０で更新された軸ずれ量確率分布に基づいて軸ずれ量を算出し、Ｓ７０に移行する。具体的には、ＣＰＵ１１は、図１５に示すように、Ｓ５０で更新された軸ずれ量確率分布Ｐｔを正規分布ＮＤと仮定して、平均値および標準偏差を算出する。
【0084】
すなわち、ＣＰＵ１１は、軸ずれ量確率分布Ｐｔのピーク位置ｗｐｅａｋを軸ずれ量δとする。またＣＰＵ１１は、下式（３）に示すように、ピーク値Ｐｅａｋに対して６０％の値となる第１位置ｗｍと第２位置ｗｐとの差の半分を標準偏差σとする。なお、第１位置ｗｐの値は第２位置ｗｍの値より大きい。
【0085】
σ ＝ （ｗｐ－ｗｍ）／２ ・・・（３）
Ｓ７０に移行すると、ＣＰＵ１１は、図４に示すように、予め設定された補正開始条件が成立したか否かを判断する。本実施形態の補正条件は、特徴量リスト内に記憶されている特徴量の数が予め設定された第１補正判定値以上であり、且つ、Ｓ６０で算出された標準偏差σが予め設定された第２補正判定値以下であることである。
Ｓ７０に移行すると、ＣＰＵ１１は、図４に示すように、予め設定された補正開始条件が成立したか否かを判断する。本実施形態の補正条件は、特徴量リスト内に記憶されている特徴量の数が予め設定された第１補正判定値以上であり、且つ、Ｓ６０で算出された標準偏差σが予め設定された第２補正判定値以下であることである。
【0086】
ここで、補正開始条件が成立していない場合には、ＣＰＵ１１は、メインルーチンを終了する。一方、補正開始条件が成立している場合には、ＣＰＵ１１は、Ｓ８０にて、軸ずれ補正量を算出する。軸ずれ補正量は、今回のフレームで算出された軸ずれ量から、前回のフレームまでに算出された軸ずれ量の中で最新の軸ずれ量を減算することによって算出される。
【0087】
そしてＣＰＵ１１は、Ｓ９０にて、Ｓ８０で算出された軸ずれ補正量だけ選択視野をＺ軸方向に沿って移動させることによって、軸ずれを補正する。なお、選択視野におけるＺ軸方向に沿った移動は、１観測レイヤ分の垂直角度毎にしか行えない。このため、Ｓ８０にて、軸ずれ補正量は、１観測レイヤ分の垂直角度を単位として算出される。
【0088】
Ｓ９０の処理が終了すると、ＣＰＵ１１は、メインルーチンを終了する。
このように構成された電子制御装置１０は、自車両ＶＨに搭載されたレーダ装置２０による検出結果に基づいて、車両の前方で走行する先行車両のリフレクタの距離および観測レイヤを取得する。
このように構成された電子制御装置１０は、自車両ＶＨに搭載されたレーダ装置２０による検出結果に基づいて、車両の前方で走行する先行車両のリフレクタの距離および観測レイヤを取得する。
【0089】
電子制御装置１０は、リフレクタの距離および垂直角度とリフレクタの存在確率との対応関係が予め設定された尤度モデルＬＭと、複数のリフレクタの距離および観測レイヤとを用いて、最尤推定法により、レーダ装置２０の軸ずれ量を算出する。
【0090】
このような電子制御装置１０は、尤度モデルを用いることによって、先行車両のリフレクタの設置位置のばらつきを考慮して軸ずれ量を算出することができるため、軸ずれ量の推定精度を向上させることができる。
【0091】
また電子制御装置１０は、リフレクタ高が予め設定された除外判定高を超えているか否かを判断し、リフレクタ高が除外判定高を超えている場合には、リフレクタ高が除外判定高を超えているリフレクタの距離および観測レイヤを、軸ずれ量の算出から除外する。
【0092】
また電子制御装置１０は、レーダ装置２０に対して第２先行車両ＰＶ２の一部が第１先行車両ＰＶ１に隠されているオクルージョンが発生しているか否かを判断し、オクルージョンが発生していると判断した場合には、第２先行車両ＰＶ２のリフレクタの距離および観測レイヤを、軸ずれ量の算出から除外する。
【0093】
また電子制御装置１０は、先行車両が虚像であるか否かを判断し、先行車両が虚像であると判断した場合には、虚像であると判断された先行車両のリフレクタの距離および観測レイヤを、軸ずれ量の算出から除外する。
【0094】
これにより、電子制御装置１０は、軸ずれ量の算出に不適切なリフレクタを除外して、軸ずれ量の算出を行うことができるため、軸ずれ量の推定精度を更に向上させることができる。
【0095】
また、リフレクタの存在確率は、販売されている車両のリフレクタの設置高さの分布と、レーダ設置高とに基づいて設定される。これにより、電子制御装置１０は、軸ずれ量の推定精度を更に向上させることができる。
【0096】
以上説明した実施形態において、電子制御装置１０は軸ずれ推定装置に相当し、Ｓ１１０，Ｓ１２０はリフレクタ情報取得部としての処理に相当し、Ｓ４０～Ｓ６０は軸ずれ量算出部としての処理に相当し、リフレクタの距離および観測レイヤはリフレクタ位置情報に相当する。
【0097】
また、Ｓ１３０はリフレクタ高除外部としての処理に相当し、Ｓ１４０は＊オクルージョン除外部としての処理に相当し、Ｓ１５０は虚像除外部としての処理に相当する。
以上、本開示の一実施形態について説明したが、本開示は上記実施形態に限定されるものではなく、種々変形して実施することができる。
以上、本開示の一実施形態について説明したが、本開示は上記実施形態に限定されるものではなく、種々変形して実施することができる。
【0098】
本開示に記載の電子制御装置１０およびその手法は、コンピュータプログラムにより具体化された一つ乃至は複数の機能を実行するようにプログラムされたプロセッサおよびメモリを構成することによって提供された専用コンピュータにより、実現されてもよい。あるいは、本開示に記載の電子制御装置１０およびその手法は、一つ以上の専用ハードウェア論理回路によってプロセッサを構成することによって提供された専用コンピュータにより、実現されてもよい。もしくは、本開示に記載の電子制御装置１０およびその手法は、一つ乃至は複数の機能を実行するようにプログラムされたプロセッサおよびメモリと一つ以上のハードウェア論理回路によって構成されたプロセッサとの組み合わせにより構成された一つ以上の専用コンピュータにより、実現されてもよい。また、コンピュータプログラムは、コンピュータにより実行されるインストラクションとして、コンピュータ読み取り可能な非遷移有形記録媒体に記憶されてもよい。電子制御装置１０に含まれる各部の機能を実現する手法には、必ずしもソフトウェアが含まれている必要はなく、その全部の機能が、一つあるいは複数のハードウェアを用いて実現されてもよい。
【0099】
上記実施形態における１つの構成要素が有する複数の機能を、複数の構成要素によって実現したり、１つの構成要素が有する１つの機能を、複数の構成要素によって実現したりしてもよい。また、複数の構成要素が有する複数の機能を、１つの構成要素によって実現したり、複数の構成要素によって実現される１つの機能を、１つの構成要素によって実現したりしてもよい。また、上記実施形態の構成の一部を省略してもよい。また、上記実施形態の構成の少なくとも一部を、他の上記実施形態の構成に対して付加または置換してもよい。
【0100】
上述した電子制御装置１０の他、当該電子制御装置１０を構成要素とするシステム、当該電子制御装置１０としてコンピュータを機能させるためのプログラム、このプログラムを記録した半導体メモリ等の非遷移的実体的記録媒体、軸ずれ推定方法など、種々の形態で本開示を実現することもできる。
【符号の説明】
【0101】
１０…電子制御装置、２０…レーダ装置
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】