(57)【要約】【課題】交差点における道路境界を精度よく検出することができるようにする。【解決手段】点群取得部２６によって、車両が交差点を走行したときに撮像装置１２によって取得されたステレオ画像と、測位情報とに基づいて、交差点における３次元点群を取得する。線分検出部２８によって、同一の交差点の異なる走行方向にわたって得られた３次元点群から、線分を各々検出する。線分ペア選択部３２によって、検出された線分から、２つの線分からなる線分ペアを複数選択する。双曲線パラメータ算出部３６によって、線分ペア毎に、線分ペアの２つの線分を漸近線とする片側双曲線のパラメータを算出する。最適化部３８によって、線分ペア毎に算出された片側双曲線のうち、３次元点群に最も適合する片側双曲線を、交差点における道路境界として検出する。【選択図】図１
2014115904.tif 000002
【特許請求の範囲】
【請求項１】
車両が交差点を走行したときに前記車両に搭載されたデータ取得手段によって取得された前記交差点における物体データと、前記車両の測位情報とに基づいて、実座標系の３次元座標を持つ、前記交差点における物体上の複数の点からなる点群データを取得する点群取得手段と、
前記点群取得手段によって取得された前記点群データから、線分を各々検出する線分検出手段と、
同一の交差点の異なる走行方向の前記車両の走行に対して得られた前記点群データから前記線分検出手段によって検出された線分から、２つの線分からなる線分ペアを複数選択する線分選択手段と、
前記線分選択手段によって選択された線分ペア毎に、前記線分ペアの２つの線分を漸近線とする双曲線を算出する双曲線算出手段と、
前記双曲線算出手段によって前記線分ペア毎に算出された双曲線のうち、前記点群取得手段によって前記同一の交差点の異なる走行方向の前記車両の走行に対して取得された前記点群データに最も適合する双曲線の片側曲線を、前記交差点における道路境界として検出する道路境界検出手段と、
を含む交差点道路境界検出装置。
【請求項２】
前記線分検出手段は、前記交差点の走行方向毎に、前記点群取得手段によって取得された前記点群データのうち、前記交差点の前記走行方向の前記車両の走行に対する前記点群データから、線分を各々検出する請求項１記載の交差点道路境界検出装置。
【請求項３】
前記双曲線算出手段は、前記点群取得手段によって取得された前記点群データから、前記線分ペアの各線分の周辺の点を所定個選択し、前記選択された所定個の点と、双曲線との誤差が最小となるように、前記双曲線のパラメータを算出する請求項１又は２記載の交差点道路境界検出装置。
【請求項４】
前記双曲線算出手段は、前記点群取得手段によって取得された前記点群データから、前記線分ペアの２つの線分の交点まで各線分を延長した各延長線の周辺の点を、前記延長線との距離が近い点ほど優先して選択することにより、前記所定個の点を選択し、前記選択された所定個の点と、双曲線との誤差が最小となるように、前記双曲線のパラメータを算出する請求項３記載の交差点道路境界検出装置。
【請求項５】
前記双曲線算出手段は、前記点群取得手段によって取得された前記点群データから、前記線分ペアの各線分の周辺の点を所定個選択することを反復して行い、反復毎に、前記選択された所定個の点と、双曲線との誤差が最小となるように、前記双曲線のパラメータを算出する請求項３又は４記載の交差点道路境界検出装置。
【請求項６】
前記双曲線算出手段は、前記線分選択手段によって選択された線分ペア毎に、前記線分ペアの２つの線分を漸近線とする片側双曲線を算出し、
前記道路境界検出手段は、前記交差点における複数の角に対して、前記双曲線算出手段によって算出された片側双曲線を振り分け、前記交差点の角毎に、前記角に振り分けられた片側双曲線のうち、前記点群取得手段によって取得された前記点群データに最も適合する片側曲線を、前記交差点における前記角の道路境界として検出する請求項１〜請求項５の何れか１項記載の交差点道路境界検出装置。
【請求項７】
前記データ取得手段は、前記車両の周辺を撮像してステレオ画像を生成する撮像手段であって、
前記点群取得手段は、前記車両が交差点を走行したときに撮像された前記ステレオ画像から抽出される相対座標系の３次元座標を持つ複数の点と、前記車両の測位情報とに基づいて、前記点群データを取得する請求項１〜請求項６の何れか１項記載の交差点道路境界検出装置。
【請求項８】
前記データ取得手段は、前記車両の周辺に存在する物体上の点を検出するレーザレーダであって、
前記点群取得手段は、前記車両が交差点を走行したときに検出された前記交差点における物体上の点と、前記車両の測位情報とに基づいて、前記点群データを取得する請求項１〜請求項６の何れか１項記載の交差点道路境界検出装置。
【請求項９】
コンピュータを、
車両が交差点を走行したときに前記車両に搭載されたデータ取得手段によって取得された前記交差点における物体データと、前記車両の測位情報とに基づいて、実座標系の３次元座標を持つ、前記交差点における物体上の複数の点からなる点群データを取得する点群取得手段、
前記点群取得手段によって取得された前記点群データから、線分を各々検出する線分検出手段、
同一の交差点の異なる走行方向の前記車両の走行に対して得られた前記点群データから前記線分検出手段によって検出された線分から、２つの線分からなる線分ペアを複数選択する線分選択手段、
前記線分選択手段によって選択された線分ペア毎に、前記線分ペアの２つの線分を漸近線とする双曲線を算出する双曲線算出手段、及び
前記双曲線算出手段によって前記線分ペア毎に算出された双曲線のうち、前記点群取得手段によって前記同一の交差点の異なる走行方向の前記車両の走行に対して取得された前記点群データに最も適合する双曲線の片側曲線を、前記交差点における道路境界として検出する道路境界検出手段
として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、交差点道路境界検出装置及びプログラムに係り、特に、交差点における道路境界を検出する交差点道路境界検出装置及びプログラムに関する。
【背景技術】
【０００２】
従来より、交差点進入前と進入後の経路が与えられたとき、交差点内の経路を所定の曲率の円弧で補間する道路網データ生成装置が知られている（特許文献１）。
【０００３】
また、スプライン表現により、位置や曲率などを含めた道路形状を表現する方法が知られている(特許文献２)。
【０００４】
また、レーザ等により得た点群データから道路形状を認識する道路形状認識装置が知られている(特許文献３)。この道路形状認識装置では、道路形状を円弧としてモデル化し、最小二乗法を用いたフィッティングにより道路形状を特定している。
【０００５】
また、道路上のカーブの形状を粗い点群データから認識する道路形状認識方法が知られている（特許文献４）。この道路形状認識方法では、カーブの開始点、終了点を求め、その間の点群に円弧をあてはめることでカーブ形状を認識している。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２０１０−２６８７５号公報
【特許文献２】特開２００７−３１６６５４号公報
【特許文献３】特開２０１０−２５０７４３号公報
【特許文献４】特開２００１−１０１５９７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
上記の特許文献２に記載の表現方法では、単路における曲率を表現することができ、自由度が高いが、交差点の道路境界の曲率推定には適さない、という問題がある。
【０００８】
また、上記の特許文献３に記載の技術では、道路が円弧状の場合に道路形状を特定することができるが、交差点では、曲率半径が小さく死角が大きくなること、及び円弧に当てはまる点群が少ないことから、交差点の道路境界を検出する目的で利用することができない。
【０００９】
また、上記の特許文献４に記載の技術では、直線状領域、円弧状領域のつなぎあわせによって形状を記述しているため、カーブの開始点、終了点を求める処理が必要になる。
【００１０】
このように、従来技術では、単路のカーブであって、曲率が大きく穏やかなカーブを対象としたものが多いが、交差点の道路境界の検出では、（１）曲率が小さい、（２）形状が複雑、（３）車両などの物体物が多く検出の妨げになる、など課題が多い。
【００１１】
本発明は、上記の問題点を解決するためになされたもので、交差点における道路境界を精度よく検出することができる交差点道路境界検出装置及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【００１２】
上記の目的を達成するために本発明に係る交差点道路境界検出装置は、車両が交差点を走行したときに前記車両に搭載されたデータ取得手段によって取得された前記交差点における物体データと、前記車両の測位情報とに基づいて、実座標系の３次元座標を持つ、前記交差点における物体上の複数の点からなる点群データを取得する点群取得手段と、前記点群取得手段によって取得された前記点群データから、線分を各々検出する線分検出手段と、同一の交差点の異なる走行方向の前記車両の走行に対して得られた前記点群データから前記線分検出手段によって検出された線分から、２つの線分からなる線分ペアを複数選択する線分選択手段と、前記線分選択手段によって選択された線分ペア毎に、前記線分ペアの２つの線分を漸近線とする双曲線を算出する双曲線算出手段と、前記双曲線算出手段によって前記線分ペア毎に算出された双曲線のうち、前記点群取得手段によって前記同一の交差点の異なる走行方向の前記車両の走行に対して取得された前記点群データに最も適合する双曲線の片側曲線を、前記交差点における道路境界として検出する道路境界検出手段と、を含んで構成されている。
【００１３】
本発明に係るプログラムは、コンピュータを、車両が交差点を走行したときに前記車両に搭載されたデータ取得手段によって取得された前記交差点における物体データと、前記車両の測位情報とに基づいて、実座標系の３次元座標を持つ、前記交差点における物体上の複数の点からなる点群データを取得する点群取得手段、前記点群取得手段によって取得された前記点群データから、線分を各々検出する線分検出手段、同一の交差点の異なる走行方向の前記車両の走行に対して得られた前記点群データから前記線分検出手段によって検出された線分から、２つの線分からなる線分ペアを複数選択する線分選択手段、前記線分選択手段によって選択された線分ペア毎に、前記線分ペアの２つの線分を漸近線とする双曲線を算出する双曲線算出手段、及び前記双曲線算出手段によって前記線分ペア毎に算出された双曲線のうち、前記点群取得手段によって前記同一の交差点の異なる走行方向の前記車両の走行に対して取得された前記点群データに最も適合する双曲線の片側曲線を、前記交差点における道路境界として検出する道路境界検出手段として機能させるためのプログラムである。
【００１４】
本発明によれば、点群取得手段によって、車両が交差点を走行したときに前記車両に搭載されたデータ取得手段によって取得された前記交差点における物体データと、前記車両の測位情報とに基づいて、実座標系の３次元座標を持つ、前記交差点における物体上の複数の点からなる点群データを取得する。線分検出手段によって、前記点群取得手段によって取得された前記点群データから、線分を各々検出する。
【００１５】
そして、線分選択手段によって、同一の交差点の異なる走行方向の前記車両の走行に対して得られた前記点群データから前記線分検出手段によって検出された線分から、２つの線分からなる線分ペアを複数選択する。双曲線算出手段によって、前記線分選択手段によって選択された線分ペア毎に、前記線分ペアの２つの線分を漸近線とする双曲線を算出する。
【００１６】
そして、道路境界検出手段によって、前記双曲線算出手段によって前記線分ペア毎に算出された双曲線のうち、前記点群取得手段によって前記同一の交差点の異なる走行方向の前記車両の走行に対して取得された前記点群データに最も適合する双曲線の片側曲線を、前記交差点における道路境界として検出する。
【００１７】
このように、同一の交差点の異なる走行方向の車両の走行に対して得られた点群データから線分を検出し、線分ペア毎に、線分ペアの２つの線分を漸近線とする双曲線を算出し、点群データに最も適合する双曲線の片側曲線を、交差点における道路境界として検出することにより、交差点における道路境界を精度よく検出することができる。
【００１８】
本発明に係る線分検出手段は、前記交差点の走行方向毎に、前記点群取得手段によって取得された前記点群データのうち、前記交差点の前記走行方向の前記車両の走行に対する前記点群データから、線分を各々検出するようにすることができる。
【００１９】
本発明に係る双曲線算出手段は、前記点群取得手段によって取得された前記点群データから、前記線分ペアの各線分の周辺の点を所定個選択し、前記選択された所定個の点と、双曲線との誤差が最小となるように、前記双曲線のパラメータを算出するようにすることができる。また、上記の双曲線算出手段は、前記点群取得手段によって取得された前記点群データから、前記線分ペアの２つの線分の交点まで各線分を延長した各延長線の周辺の点を、前記延長線との距離が近い点ほど優先して選択することにより、前記所定個の点を選択し、前記選択された所定個の点と、双曲線との誤差が最小となるように、前記双曲線のパラメータを算出するようにすることができる。
【００２０】
また、上記の双曲線算出手段は、前記点群取得手段によって取得された前記点群データから、前記線分ペアの各線分の周辺の点を所定個選択することを反復して行い、反復毎に、前記選択された所定個の点と、双曲線との誤差が最小となるように、前記双曲線のパラメータを算出するようにすることができる。
【００２１】
本発明に係る双曲線算出手段は、前記線分選択手段によって選択された線分ペア毎に、前記線分ペアの２つの線分を漸近線とする片側双曲線を算出し、前記道路境界検出手段は、前記交差点における複数の角に対して、前記双曲線算出手段によって算出された片側双曲線を振り分け、前記交差点の角毎に、前記角に振り分けられた片側双曲線のうち、前記点群取得手段によって取得された前記点群データに最も適合する片側曲線を、前記交差点における前記角の道路境界として検出するようにすることができる。
【００２２】
上記のデータ取得手段は、前記車両の周辺を撮像してステレオ画像を生成する撮像手段であって、前記点群取得手段は、前記車両が交差点を走行したときに撮像された前記ステレオ画像から抽出される相対座標系の３次元座標を持つ複数の点と、前記車両の測位情報とに基づいて、前記点群データを取得するようにすることができる。
【００２３】
上記のデータ取得手段は、前記車両の周辺に存在する物体上の点を検出するレーザレーダであって、前記点群取得手段は、前記車両が交差点を走行したときに検出された前記交差点における物体上の点と、前記車両の測位情報とに基づいて、前記点群データを取得するようにすることができる。
【発明の効果】
【００２４】
以上説明したように、本発明の交差点道路境界検出装置及びプログラムによれば、同一の交差点の異なる走行方向の車両の走行に対して得られた点群データから線分を検出し、線分ペア毎に、線分ペアの２つの線分を漸近線とする双曲線を算出し、点群データに最も適合する双曲線の片側曲線を、交差点における道路境界として検出することにより、交差点における道路境界を精度よく検出することができる、という効果が得られる。
【図面の簡単な説明】
【００２５】
【図１】本発明の実施の形態に係る交差点道路境界検出装置を示すブロック図である。
【図２】同一交差点における３次元点群の鳥瞰図を示すイメージ図である。
【図３Ａ】同一走行方向に対する３次元点群から検出された線分を示す図である。
【図３Ｂ】同一走行方向に対する３次元点群から検出された線分を示す図である。
【図４】部分点群の候補を説明するための図である。
【図５】部分点群を示す図である。
【図６】双曲線パラメータを説明するための図である。
【図７】４方向の片側双曲線を示す図である。
【図８】適合度を説明するための図である。
【図９】適合度を説明するための図である。
【図１０】本発明の実施の形態に係る交差点道路境界検出装置における交差点データ取得処理ルーチンの内容を示すフローチャートである。
【図１１】本発明の実施の形態に係る交差点道路境界検出装置における交差点道路境界検出処理ルーチンの内容を示すフローチャートである。
【図１２】（Ａ）ステレオ画像に含まれる入力画像の例を示す図、及び（Ｂ）交差点の道路境界の検出結果の例を示す図である。
【図１３】３次元点群の例を示す図である。
【図１４】交差点の道路境界が登録された地図の例を示す図である。
【発明を実施するための形態】
【００２６】
以下、図面を参照して本発明の実施の形態を詳細に説明する。
【００２７】
図１に示すように、本発明の実施の形態に係る交差点道路境界検出装置１０は、道路領域を含む自車両前方を撮像してステレオ画像を出力する撮像装置１２と、自車両位置を測定するＧＰＳを含む測位装置１４と、撮像装置１２から得られるステレオ画像及び測位装置１４による測位情報に基づいて、交差点における道路境界を検出する処理を実行するコンピュータ１６と、コンピュータ１６による検出結果を出力する出力部１８とを含んで構成されている。なお、撮像装置１２が、データ取得手段の一例であり、ステレオ画像が、物体データの一例である。
【００２８】
撮像装置１２は、車両に搭載され、道路領域を含む自車両前方を撮像してステレオ画像を出力する。撮像装置１２は、自車両前方を撮像し、画像の画像信号を生成する２つの撮像部（図示省略）と、２つの撮像部で生成された画像信号をＡ／Ｄ変換するＡ／Ｄ変換部（図示省略）と、Ａ／Ｄ変換された画像信号を一時的に格納するための画像メモリ（図示省略）とを備えている。撮像装置１２は、撮像した２枚の画像を、ステレオ画像としてコンピュータ１６に出力する。
【００２９】
測位装置１４は、自車両に搭載されたＧＰＳ装置を用いて、実座標系における自車両の位置を測定し、測定結果をコンピュータ１６に出力する。
【００３０】
コンピュータ１６は、ＣＰＵと、ＲＡＭと、後述する交差点道路境界検出処理ルーチンを実行するためのプログラムを記憶したＲＯＭと、これらを接続するバスとを含んで構成されている。このコンピュータ１６をハードウエアとソフトウエアとに基づいて定まる機能実現手段毎に分割した機能ブロックで説明すると、実座標系の道路平面及び交差点を表わす地図データを記憶する地図データベース２０と、データ取得部２２と、データ記憶部２４と、点群取得部２６と、線分検出部２８と、点群線分記憶部３０と、線分ペア選択部３２と、部分点群選択部３４と、双曲線パラメータ算出部３６と、最適化部３８と、を含んだ構成で表すことができる。なお、最適化部３８は、道路境界検出手段の一例である。
【００３１】
データ取得部２２は、測位装置１４により得られた測位情報と、地図データベース２０の地図データとに基づいて、自車両が交差点を走行しているか否かを判定し、自車両が交差点を走行していると判定している間、撮像装置１２によって撮像された複数フレームのステレオ画像を取得して、複数のフレームの各々に対応して測位装置１４により得られた測位情報と共にデータ記憶部２４に格納すると共に、測位装置１４により得られた測位情報と、地図データベース２０の地図データとに基づいて、自車両が走行している交差点を識別するための情報、及び交差点における走行方向（方角）を示す情報を取得して、上記のステレオ画像と対応付けてデータ記憶部２４に格納する。
【００３２】
データ記憶部２４には、自車両が走行している間に取得したステレオ画像が多数記憶されており、各ステレオ画像には、測位情報、交差点識別情報、及び走行方向情報が対応付けられて記憶されている。
【００３３】
点群取得部２６は、検出対象の交差点についてデータ記憶部２４に記憶されたステレオ画像のうち、走行方向情報が同一である複数のステレオ画像の各々から、エッジ点を抽出し、ステレオマッチングにより得られた視差から、撮像装置１２を基準としたエッジ点の３次元位置を求める。また、点群取得部２６は、検出対象の交差点について走行方向情報が同一となる複数のフレームのステレオ画像について求められた複数の３次元点を、複数のフレームの各々に対応して測位装置１４により得られた測位情報を用いて統合することにより、実空間上の３次元座標を持つ３次元点群を求める（図２参照）。
【００３４】
上記のように３次元点群を求める処理を、検出対象の交差点における走行方向情報毎に行い、検出対象の交差点における走行方向毎に、３次元点群を求め、点群線分記憶部３０に格納する。
【００３５】
なお、画像中のエッジ点は従来既知の手法を用いて抽出すればよく、例えば、Cannyのアルゴリズム等により抽出する。
【００３６】
また、撮像装置１２を基準とした３次元点の位置は、キャリブレーションされたカメラを用い、以下の文献の方法により求めることができる。
【００３７】
参考文献：Ｍ  Ｂｅｒｔｏｚｚ，Ａ  Ｂｒｏｇｇｉ，“Ｓｔｅｒｅｏ  ｉｎｖｅｒｓｅ  ｐｅｒｓｐｅｃｔｉｖｅ  ｍａｐｐｉｎｇ：ｔｈｅｏｒｙ  ａｎｄ  ａｐｐｌｉｃａｔｉｏｎｓ”，ｉｍａｇｅ  Ｖｉｓ．Ｃｏｍｐｕｔ．，ｖｏｌ．８，ｎｏ．１６，ｐｐ．５８５−５９０，１９９８．
【００３８】
線分検出部２８は、検出対象の交差点について走行方向情報が同一となる複数のフレームのステレオ画像から取得した３次元点群から、ガードレール等を表す道路境界の平面を検出し、検出した平面の高さ成分を無視することにより、道路境界候補となる線分を検出する。なお、ガードレール等を表わす道路境界の平面の検出には一般的な直線のハフ変換、または参考文献（石田皓之、目黒淳一、小島祥子、内藤貴志、“Ｃｏｎｆｏｒｍａｌ  Ｇｅｏｍｅｔｒｉｃ  Ａｌｇｅｂｒａを用いた道路物理境界検出”、信学技報ＰＲＭＵ２０１２−３６，ｐｐ．４３−４８，２０１２．）に記載された手法を用いればよい。
【００３９】
線分検出部２８は、検出対象の交差点について走行方向情報が同一となる複数のフレームのステレオ画像から取得した３次元点群から、道路境界候補となる複数の線分を検出して、点群線分記憶部３０に格納する（図３Ａ参照）。
【００４０】
上記のように道路境界候補となる線分を検出する処理を、検出対象の交差点における走行方向情報毎に行い、検出対象の交差点における走行方向毎に、線分を検出し、点群線分記憶部３０に格納する（図３Ａ、図３Ｂ参照）。
【００４１】
なお、本実施の形態では、上記の線分検出処理を、検出対象の交差点について、当該交差点を通過する各回の走行ごとに独立に行い、複数の線分を検出する。
【００４２】
線分ペア選択部３２は、点群線分記憶部３０に記憶された、検出対象の交差点における複数の走行方向の各々に対する線分から、２つの線分を、処理対象の線分ペアとして複数選択する。検出対象の交差点における線分ペアは、検出された線分がｎ本であるとすると、ｎ（ｎ−１）／２組存在する。
【００４３】
選択された全ての線分ペアについて、後述する部分点群選択部３４及び双曲線パラメータ算出部３６の処理で、線分ペアの２つの線分を漸近線とする片側双曲線を算出する処理を行い、各回の走行を統合する。
【００４４】
部分点群選択部３４は、まず、選択された線分ペアについて、当該線分ペアの２つの線分の交点まで、２つの線分の各々を延長した２つの延長線から、所定範囲内にある点群を部分点群の候補として選択する（図４参照）。このときの範囲は２ｍ以内と設定すると有効である。また、交点より遠方にある点群は他車両などを表す可能性が高いので候補としない。
【００４５】
部分点群選択部３４は、上記のように選択された部分点群の候補から、ルーレット選択により所定数の点を選択し部分点群とする。線分ペアの２線分との最小距離の逆数（上限あり：Ｍｍａｘ）をルーレットのマスの大きさとし、乱数を用いて部分点群を選択する。この処理により、道路境界付近に分布する点を選択されやすくすることが可能となる（図５参照）。
【００４６】
ここで、部分点群選択のアルゴリズムの一例について説明する。ただし、点の選択されやすさをＭ、選択される点の数をＮとする。
【００４７】
ステップ１：全候補点に対し、以下の処理を実行する。
（１−１）２線分との最小距離ｄを計算する。
（１−２）Ｍ←ｍｉｎ（１／ｄ、Ｍｍａｘ）
（１−３）ルーレット上にＭ個のマスを追加し、当該候補点のインデクスをＭ個のマスに格納する。
【００４８】
ステップ２：Ｎ回、以下の処理を実行する。
（２−１）ルーレットからランダムでマスを１つ選択し、選択されたマスにインデックスが格納されている候補点を部分点群に追加する。
【００４９】
また、上記のステップ１、２からなる部分点群選択のアルゴリズムの実行を反復し、反復毎に、候補点群から、部分点群が選択される。すなわち、複数の部分点群が選択される。
【００５０】
双曲線パラメータ算出部３６は、部分点群選択部３４によって選択された部分点群に対して、最小二乗法により、部分点群との誤差が最小となる双曲線を当てはめて、双曲線パラメータを算出する。このとき、２線分の交点が双曲線の中心（ｕ、ｖ）であるという制約を与える。ここで、線分（ｘ１、ｙ１）−（ｘ２、ｙ２）と線分（ｘ３、ｙ３）−（ｘ４、ｙ４）の交点（ｕ、ｖ）は、以下の（１）式、（２）式で求められる。
【００５１】
【数１】
2014115904.tif 000003
【００５２】
また、求める双曲線の式は以下の（３）式となる。
【００５３】
【数２】
2014115904.tif 000004
【００５４】
双曲線パラメータ算出部３６は、以下の（４）式により、直接、双曲線のパラメータ（ａ、ｂ、ｄ）を算出する。
【００５５】
【数３】
2014115904.tif 000005
【００５６】
双曲線パラメータ算出部３６は、軸の回転角θおよび漸近線の角度φを、以下の（５）式、（６）式により算出する（図６参照）。
【００５７】
【数４】
2014115904.tif 000006
【００５８】
双曲線パラメータ算出部３６は、以下の（７）式で求められる４方向の双曲線の媒介変数表示ｔに基づいて、４方向の片側双曲線のなかで部分点群と片側双曲線の最短距離の和が最小となる方向１つを選択する（図７参照）。
【００５９】
【数５】
2014115904.tif 000007
【００６０】
双曲線パラメータ算出部３６は、上記のように算出された双曲線パラメータ（ａ、ｂ、ｄ、θ、φ）及び選択された片側双曲線の方向を、線分ペアの２線分を漸近線とし、部分点群に対して当てはめた片側双曲線のパラメータとして出力する。また、双曲線パラメータ算出部３６は、部分点群選択部３４によって選択された部分点群毎に、片側双曲線のパラメータを算出して出力する。
【００６１】
最適化部３８は、双曲線パラメータ算出部３６によって求められた片側双曲線のパラメータと、点群線分記憶部３０に記憶された、検出対象の交差点についての３次元点群とに基づいて、双曲線パラメータ算出部３６によって求められたパラメータが表わす片側双曲線の３次元点群への適合度を算出する。ここで、適合度は、３次元点群の全体集合のなかで、片側双曲線と所定の距離内にある点の数とする（図８、図９参照）。なお、場所ごとに点の数に偏りがある場合、所定範囲内（たとえば１ｍ以内）の点群をまとめて１つとして扱ってもよい。
【００６２】
なお、点と片側双曲線との距離を求めるには、双曲線を所定の小間隔で線分に細分し、点群中の点とサンプリングされた線分との最小距離を求める。この距離が所定の閾値を下回る点の数を数え、適合度とする。閾値はたとえば１０ｃｍとすると有効である。
【００６３】
また、最適化部３８は、双曲線パラメータ算出部３６によって部分点群毎に求められた片側双曲線のパラメータの各々に対して、適合度を算出する。
【００６４】
最適化部３８は、全ての線分ペア、全ての部分点群に対して算出した片側双曲線のパラメータの各々の適合度のうち、適合度最大となる片側双曲線のパラメータが示す片側双曲線を、検出対象の交差点の道路境界として採用し、検出対象の交差点の道路境界の検出結果とする。
【００６５】
ここで、上記の処理は、交差点の一角の境界を求める処理であるが、実際は１つの交差点につき、複数の角が存在し、複数の境界が存在しうる。たとえば十字交差の交差点であれば、最大４つの境界が検出される。そこで、本実施の形態では、自車位置の測位結果を用いて交差点の中心位置を算出し、それを基準として、算出された片側双曲線のパラメータを、片側双曲線のパラメータに基づいて、複数の角（例えば４つの角）に振り分け、交差点の角ごとに、当該角に振り分けられた片側双曲線のうち、適合度最大となる片側双曲線を、検出対象の交差点の当該方角の道路境界として検出する。
【００６６】
次に、本発明の実施の形態に係る交差点道路境界検出装置１０の作用について説明する。まず、撮像装置１２によって、自車両前方の撮像が開始されると共に、測位装置１４によって、自車両の位置が随時測定されると、コンピュータ１６において、図１０に示す交差点データ取得処理ルーチンが繰り返し実行される。
【００６７】
ステップ１００で、測位装置１４で得られた測位情報と、地図データベース２０の地図データとに基づいて、交差点を走行しているか否かを判定する。交差点を走行していると判定されると、ステップ１０２へ進む。
【００６８】
ステップ１０２では、測位装置１４で得られた所定の距離範囲分の測位情報、及び撮像装置１２で撮像された同一の距離範囲分の複数のフレームのステレオ画像を取得する。
【００６９】
次に、ステップ１０４で、上記ステップ１０２で取得した複数のフレームのステレオ画像及び測位情報を、測位装置１４で得られた測位情報と地図データベース２０の地図データとに基づいて得られる交差点識別情報及び交差点方向情報と共に、データ記憶部２４に格納し、交差点データ取得処理ルーチンを終了する。
【００７０】
自車両の走行中に、上記の交差点データ取得処理ルーチンが繰り返し実行されると、各交差点における複数の走行方向に対するデータが、データ記憶部２４に蓄積される。
【００７１】
そして、コンピュータ１６において、図１１に示す道路境界検出処理ルーチンが、各交差点を検出対象として実行される。
【００７２】
まず、ステップ１１０において、検出対象の交差点についてデータ記憶部２４に記憶されたステレオ画像及び測位情報を取得する。ステップ１１２では、上記ステップ１１０で取得した各フレームのステレオ画像から、エッジ点を抽出する。ステップ１１４では、上記ステップ１１２で抽出した各フレームのステレオ画像のエッジ点群の各点について３次元位置を求めると共に、上記ステップ１１０で取得した測位情報を用いて、各フレームのエッジ点群を統合して、３次元点群とする。
【００７３】
そして、ステップ１１６において、上記ステップ１１４で得られた３次元点群のうち、各走行に対する３次元点群から、線分を各々検出し、ステップ１１８において、上記ステップ１１４で得られた３次元点群及び上記ステップ１１６で検出された各線分を、点群線分記憶部３０に格納する。
【００７４】
次のステップ１２０では、上記ステップ１１８で点群線分記憶部３０に格納された線分のうち、２つの線分からなる線分ペアを選択し、処理対象の線分ペアとして設定する。ステップ１２２において、上記ステップ１１８で点群線分記憶部３０に格納された３次元点群から、上記ステップ１２０で設定された処理対象の線分ペアの２つの線分の延長線の周辺の点を、部分点群として選択する。ここで、複数の部分点群が選択される。
【００７５】
そして、ステップ１２４において、上記ステップ１２２で選択された部分点群に基づいて、片側双曲線を示すパラメータを算出し、ステップ１２６で、上記ステップ１２４で算出された片側双曲線パラメータの適合度を算出する。ここで、上記ステップ１２４、１２６の処理は、選択された複数の部分点群の各々について行われる。
【００７６】
ステップ１２８では、上記ステップ１１８で点群線分記憶部３０に格納された線分における全ての線分ペアについて上記ステップ１２０〜ステップ１２６の処理が実行されたか否かを判定し、全ての線分ペアについて上記ステップ１２０〜ステップ１２６の処理が実行されたと判定された場合には、ステップ１３０へ進む。一方、上記ステップ１２０〜ステップ１２６の処理が実行されていない線分ペアが存在する場合には、上記ステップ１２０へ戻り、当該線分ペアを処理対象の線分ペアとして設定する。
【００７７】
ステップ１３０では、上記ステップ１１０で取得した測位情報に基づいて、交差点の中心位置を算出し、算出された中心位置を基準として、上記ステップ１２４で算出された片側双曲線パラメータを、複数の角（例えば４つの角）に分類し、片側双曲線パラメータの各々を、複数の角の何れかに振り分ける。
【００７８】
そして、ステップ１３２では、交差点の角毎に、当該角に振り分けられた片側双曲線パラメータのうち、上記ステップ１２６で算出された適合度が最大となる、片側双曲線パラメータを選択する。次のステップ１３４では、上記ステップ１３２で交差点の角毎に選択された片側双曲線パラメータが示す片側双曲線を、交差点における角毎の道路境界の検出結果として出力し、道路境界検出処理ルーチンを終了する。
【００７９】
道路境界検出処理ルーチンによって検出された交差点における角毎の道路境界は、例えば、地図上に登録される。
【００８０】
たとえば、図１２（Ａ）に示すような入力画像を含むステレオ画像が入力された場合には、図１３に示すような、点群データが取得され、図１２（Ｂ）に示すように、交差点における道路境界であるガードレールが検出される。なお、図１２（Ａ）の画像中の丸印が示すガードレールと、図１２（Ｂ）の丸印が示すガードレールとは対応しており、図１２（Ａ）の画像を撮像したときの車両の位置及び撮像方向が、図１２（Ｂ）中に示されている。
【００８１】
以上説明したように、本発明の実施の形態の交差点道路境界検出装置によれば、同一の交差点の異なる走行方向の車両の走行にわたって得られた３次元点群から線分を検出し、線分ペア毎に、線分ペアの２つの線分を漸近線とする片側双曲線のパラメータを算出し、３次元点群に最も適合する片側双曲線パラメータが示す片側双曲線を、交差点における道路境界として検出することにより、交差点における道路境界を精度よく検出することができる。また、交差点における道路境界を検出して、地図に登録することにより、図１４に示すような高精度な地図を生成することができる。
【００８２】
カメラを搭載した車両で交差点を複数方向に走行し、複数方向に走行して得られた点群を統合利用することで、交差点の道路境界を検出できる。単一方向に走行するだけでは死角が多いため全体を検出できないが、複数方向にわたって取得したデータを統合することで交差点の道路境界の形状を認識できる。また、交差点付近は人、車、バイクなどが停留しやすく、道路境界以外の物体を誤検出しやすいが、本実施の形態に係る交差点道路境界検出装置により、誤検出を抑制でき、静止物の境界のみ検出可能となる。
【００８３】
また、交差点周囲の道路境界を個別の線分として扱うのではなく、交差する道路の線分をそれぞれ漸近線とする単一の片側双曲線として、交差点の道路境界を表現するため、交差する両方の道路境界、および角の形状を併せてモデル化し最適化することができる。
【００８４】
また、本実施の形態では、交差点周囲の道路の直線状領域、円弧状領域をまとめて扱うことのできる双曲線モデルで近似するため、モデルの当てはめが容易になる。また、曲率の小さい直角状の交差点境界も併せて同様の枠組で検出することができる。
【００８５】
また、もともと、双曲線近似はノイズの影響を受けやすく多くの課題がある。本実施の形態では、複数の走行で得られた点群のうち、漸近線となる線分ペアに近い点群を統合して用いることで、車載カメラによって得たノイズを多く含む点群であっても比較的安定して双曲線フィッティング可能である。また、他物体の影響を受けにくいため誤検出を抑制できる。
【００８６】
また、交差点の角ごとに、適合度が最大となる片側双曲線を求めることで、交差点の角ごとの道路境界が得られる。
【００８７】
なお、上記の実施の形態では、車両走行時に、交差点のステレオ画像及び測位情報の取得及び記憶のみを行い、後で、３次元点群を求める場合を例に説明したが、これに限定されるものではない。例えば、車両走行時に、ステレオ画像及び測位情報の取得及び記憶を行うと共に、エッジ点の抽出及び実座標系の３次元座標の算出を行うようにしてもよい。
【００８８】
また、交差点道路境界検出装置が車両に搭載されている場合を例に説明したが、撮像装置、測位装置、データ取得部、地図データベース、及びデータ記憶部を含む車載器とは別の装置で、交差点道路境界検出装置を構成してもよい。この場合には、車載器のデータ記憶部に記憶されたデータを、ネットワーク等を介して、交差点道路境界検出装置に入力するようにすればよい。
【００８９】
また、ステレオカメラである撮像装置の代わりに、レーザレーダを用いて、点群を取得してもよい。この場合には、レーザレーダを用いて、道路領域を含む自車両前方に存在する物体上の３次元点を各々検出し、コンピュータに出力するようにすればよい。また、レーザレーダが、データ取得手段の一例であり、検出される物体上の３次元点が、物体データの一例である。
【００９０】
本発明のプログラムは、記憶媒体に格納して提供するようにしてもよい。
【符号の説明】
【００９１】
１０   交差点道路境界検出装置
１２   撮像装置
１４   測位装置
１６   コンピュータ
２０   地図データベース
２２   データ取得部
２４   データ記憶部
２６   点群取得部
２８   線分検出部
３０   点群線分記憶部
３２   線分ペア選択部
３４   部分点群選択部
３６   双曲線パラメータ算出部
３８   最適化部
【図１】
2014115904.tif 000008
【図２】
2014115904.tif 000009
【図６】
2014115904.tif 000010
【図７】
2014115904.tif 000011
【図１０】
2014115904.tif 000012
【図１１】
2014115904.tif 000013
【図３Ａ】
2014115904.tif 000014
【図３Ｂ】
2014115904.tif 000015
【図４】
2014115904.tif 000016
【図５】
2014115904.tif 000017
【図８】
2014115904.tif 000018
【図９】
2014115904.tif 000019
【図１２】
2014115904.tif 000020
【図１３】
2014115904.tif 000021
【図１４】
2014115904.tif 000022