(57)【要約】  プローブデータから道路の交通量を予測する交通量予測システムでは、道路種別だけで、速度と交通量の関係を定義することは難しいため、求める交通量の精度に限界がある。  プローブデータから自由速度、臨界速度、飽和密度、臨界密度などの道路特性を推定し、この道路特性からプローブデータに対応したクラスタに分類し、各クラスタを構成するプローブデータの速度と加速度を元に、クラスタ毎の走行パターンを作成する。次に、プローブデータを受信した時、そのプローブデータの道路がいずれかのクラスタに属していない場合には、プローブデータから生成した走行パターンと、クラスタ毎の走行パターンを比較して、類似性の高いクラスタの道路特性を取得し、取得した道路特性を使って、プローブデータが収集された道路の交通量を予測する。
2013128486.tif 000002
【特許請求の範囲】
【請求項１】
プローブデータを使って交通量を予測する交通量予測システムであって、
プローブデータから、交通量を表すモデル式で用いられる自由速度、臨界速度、飽和密度、臨界密度のいずれか一つを含む道路特性を推定する道路特性推定装置と、
推定した前記道路特性を互いに類似する道路特性毎に分類したクラスタを生成する分類装置と、
前記クラスタを構成する道路特性の推定に用いた前記プローブデータの速度とその速度変化量により各クラスタについてそのクラスタに属する道路における車両の走行パターンを作成する走行パターン作成装置と、
前記クラスタに属していない道路のプローブデータについて当該プローブデータから作成した車両の走行パターンと、前記走行パターン作成装置により作成したクラスタにおける走行パターンを比較して類似性の高い走行パターンのクラスタの道路特性を取得する道路特性取得装置と、
前記道路特性取得装置で取得した道路特性を元に前記モデル式から前記道路の交通量を予測する交通量予測装置と、
を備えることを特徴とする交通量予測システム。
【請求項２】
請求項１に記載の交通量予測システムにおいて、前記道路特性推定装置は、蓄積した前記プローブデータを統計処理して求めた車両の速度を使って前記モデル式の道路特性を推定することを特徴とする交通量予測システム。
【請求項３】
請求項２に記載の交通量予測システムにおいて、前記走行パターン作成装置は、走行パターンとして、車両の速度と進行方向の速度変化量の関係と、車両の速度と進行方向に対する横方向の速度変化量の関係について前記プローブデータの頻度を求めることを特徴とする交通量予測システム。
【請求項４】
請求項２の道路特性取得装置は、前記クラスタに属していない道路で、かつプローブデータの無い道路の場合、前記クラスタに属している道路の属性情報の内、車線数と道路幅が一致する道路数が多いクラスタの道路特性を、当該道路の道路特性として取得することを特徴とする交通量予測システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、道路の交通量を予測するシステムに関し、特に道路を走行した車両の情報を用いて交通量を予測するシステムに関する。
【背景技術】
【０００２】
日本では、交通量を計測する取り組みとして、５年に１回の間隔で、道路の交通量を人手で測定する道路交通センサスの一般交通量調査が知られている。海外でも同様に、交通量調査は人手で行われている。この様な交通量調査は、多くの作業者と時間を使うため、大きな費用を必要とする。そのため、交通量調査の頻度を上げることは難しい。【０００３】
そこで、特開２００９−１４０００７号公報（特許文献１）には、プローブカーを使って収集した車両の走行履歴を元に、道路区間の交通量を求める発明が開示されている。プローブカーの走行履歴から道路区間の交通量を正確に求めるためには、道路単位に車両の走行速度と交通量の関係を予め作成しておく必要がある。そして車両の走行速度と交通量の関係を作成するためには、実際の交通量と車両の走行速度を同時に測定しなければならない。【０００４】
しかし、全ての道路に対して、実際の交通量と車両の速度を同時に測定するのは、大きな費用を必要とすることから、特許文献１では、道路区間の速度と交通量の関係を、高速道路、国道といった道路種別毎に作成することで、道路種別毎に実際の交通量と車両の走行速度を近似することが記載されている。このように道路種別毎に、交通量と車両の速度を同時に測定すればよいため、交通量調査の費用を抑えることができる。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００９−１４０００７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかし、道路種別だけで、速度と交通量の関係を定義することは難しい。これは、車線数や信号機の数、バス専用レーンの有無、路上駐車台数など、同じ道路種別の道路でも道路の環境が異なるためである。このため、道路種別毎に実際の交通量と車両の走行速度を求めて交通量を予測した場合、その精度に限界がある。
【課題を解決するための手段】
【０００７】
本発明の交通量予測システムは、プローブデータから自由速度、臨界速度、飽和密度、臨界密度などの道路特性を推定し、この道路特性からプローブデータに対応したクラスタに分類し、各クラスタを構成するプローブデータの速度と加速度を元に、クラスタ毎の走行パターンを作成する。次に、プローブデータを受信した時、そのプローブデータの道路がいずれかのクラスタに属していない場合には、プローブデータから生成した走行パターンと、クラスタ毎の走行パターンを比較して、類似性の高いクラスタの道路特性を取得し、取得した道路特性を使って、プローブデータが収集された道路の交通量を予測する。
【発明の効果】
【０００８】
本発明によれば、プローブデータの走行パターンを求めて、予め分類しておいた道路特性のうち類似したグループの道路特性から交通量を予測するため、実際の道路環境を反映して精度良く交通量を予測することができる。
【図面の簡単な説明】
【０００９】
【図１】本発明の実施の形態に係る交通量予測システムの全体構成を示すブロック図である。
【図２】送信装置２００から送信される交通量データのフォーマットを示す図である。
【図３】プローブデータ記憶装置のデータフォーマットを示す図である。
【図４】道路特性推定装置の処理フローを示す図である。
【図５】平均速度と交通密度の関係を示す図である。
【図６】道路特性のクラスタリングの概念図を示す図である。
【図７】幅広道路の走行パターンの例を示す図である。
【図８】幅狭道路の走行パターンの例を示す図である。
【図９】走行パターン作成装置の処理フローを示す図である。
【図１０】走行パターンのテーブルを示す図である。
【図１１】記憶装置のデータフォーマットを示す図である。
【図１２】道路特性検索装置の処理フローを示す図である。
【図１３】本発明の実施の形態の変形に係る交通量予測システムの全体構成を示すブロック図である。
【図１４】受信装置で受信する速度データのデータフォーマットを示す図である。
【図１５】道路情報記憶装置のデータフォーマットを示す図である。
【図１６】道路特性検索装置の処理フローを示す図である。
【発明を実施するための形態】
【００１０】
本発明に係る交通量予測システムの実施の形態について、図面を参照して説明する。
【００１１】
本発明の実施の形態に係る交通量予測システムの全体構成を図１に示す。交通量予測システム１００は、交通量記憶装置１１０、プローブデータ記憶装置１２０、道路特性推定装置１３０、分類装置１４０、走行パターン作成装置１５０、記憶装置１６０、受信装置１７０、道路特性検索装置１８０、交通量予測装置１９０、送信装置２００を備えて構成される。
【００１２】
交通量記憶装置１１０は、ハードディスク、フラッシュメモリ等の記憶装置で、交通量情報を記憶している。交通量情報は、道路特性推定装置１３０による道路特性の推定に利用される。プローブデータ記憶装置１２０は、ハードディスク、フラッシュメモリ等の記憶装置で、プローブデータを記憶している。プローブデータは、車両の識別ＩＤ（以下、車両ＩＤ）、車両の走行する道路の識別ＩＤ（以下、道路ＩＤ）、位置情報、時刻情報、速度情報などで構成される。車両は、携帯電話、無線などの通信機器を介して、車両ＩＤ、位置情報、時刻情報、速度情報を図示されていないプローブセンタにアップリンクする。この様な車両はプローブカーと呼ばれている。
【００１３】
プローブセンタは、各車両から収集した位置情報と時刻情報と地図情報を使って、各車両が走行した道路を特定する。プローブデータ記憶装置１２０に記憶されているプローブデータと、受信装置１７０で受信するプローブデータは、このプローブセンタで処理されたデータである。プローブデータ記憶装置１２０は、道路特性推定装置１３０の道路特性の推定と、走行パターン作成装置１５０の走行パターン作成に利用される。
【００１４】
道路特性推定装置１３０は、交通量記憶装置１１０の交通量情報とプローブデータ記憶装置１２０の速度情報を使って、道路区間単位に、速度と交通密度の関係式のパラメータである道路特性データを推定する。この道路特性データは、交通量を予測する際に使用される。道路特性データは、分類装置１４０へ提供される。
【００１５】
分類装置１４０は、道路特性推定装置１３０からの道路特性データをクラスタリングして、生成されたクラスタの代表道路特性データとクラスタを構成する道路ＩＤ情報を作成する。代表道路特性データは、クラスタリングされた複数の道路特性データから求めた各クラスタの代表値とする道路特性データである。この代表道路特性データは、記憶装置１６０に格納される。また、クラスタを構成する道路特性データと対応付けられている道路ＩＤ情報は、記憶装置１６０に格納されるとともに、走行パターン作成装置１５０へ提供される。
【００１６】
走行パターン作成装置１５０は、分類装置１４０のクラスタを構成する道路ＩＤ情報と、プローブデータ記憶装置１２０のプローブデータを使って、各クラスタの走行パターンデータを生成する。生成した走行パターンデータは、記憶装置１６０に書き込む。
【００１７】
記憶装置１６０は、ハードディスク、フラッシュメモリ等の記憶装置で、分類装置１４０で作成したクラスタを構成する道路ＩＤ情報、クラスタの代表道路特性データ、走行パターン作成装置１５０で作成した走行パターンデータを記憶している。
【００１８】
受信装置１７０は、プローブセンタからのプローブデータを受信する装置である。受信するプローブデータのフォーマットは、プローブデータ記憶装置１２０に格納されるプローブデータのフォーマットと同様である。受信したプローブデータは、プローブデータ記憶装置１２０に蓄積される他、道路特性検索装置１８０や交通量予測装置１９０へ提供される。
【００１９】
道路特性検索装置１８０は、受信装置１７０からのプローブデータにより走行パターンデータを作成し、記憶装置１６０から、この走行パターンデータと類似性の高い走行パターンデータを持つクラスタを検索する。そして該当するクラスタの代表道路特性データが、交通量予測装置１９０へ提供される。
【００２０】
交通量予測装置１９０は、道路特性検索装置１８０から提供された代表道路特性データと、受信装置１７０が受信したプローブデータの走行速度のデータを使用して、交通量を予測する。予測した交通量は送信装置２００へ提供される。この送信装置２００は、交通量予測装置１９０の交通量データを、車載端末や、他のセンターシステムに送信する。
【００２１】
これらの道路特性推定装置１３０、分類装置１４０、走行パターン作成装置１５０、道路特性検索装置１８０、交通量予測装置１９０は、いずれも交通量予測システム１００に搭載されている不図示のマイクロプロセッサや、ＲＡＭ、ＲＯＭなどによって実行されるソフトウェアにより実現される。
【００２２】
図２に送信装置２００から送信される交通量データのフォーマットの例を示す。交通量データは、道路ＩＤ２１、交通量を計測した時間帯の開始時刻２２と終了時刻２３、そしてこの開始時刻から終了時刻までの間の交通量２４で構成される。この交通量は、特定の時間帯に道路を通過する車両の台数である。
【００２３】
次に、交通量記憶装置１１０の詳細を説明する。この交通量記憶装置１１０には、道路の交通量データが格納されている。しかし、全ての道路の交通量データを格納されているわけではなく、道路交通センサスなどの様にして交通量が計測された道路の情報のみ格納されている。この交通量の計測には、観測者が車両の通過台数をカウントする方法や、道路交通インフラのセンサが機械的に車両の通過台数をカウントする方法がある。交通量記憶装置１１０のデータフォーマットは、図２と同じである。
【００２４】
次に、プローブデータ記憶装置１２０の詳細を説明する。プローブデータ記憶装置１２０には、一定期間のプローブデータが格納されている。図３は、プローブデータ記憶装置１２０のデータフォーマットの例である。
【００２５】
プローブデータ記憶装置１２０のプローブデータは、プローブデータが収集された道路を識別するための道路ＩＤ３１、プローブデータを収集した車両を識別するための車両ＩＤ３２、車両がその道路へ進入した時刻である流入時刻３３とその道路から流出した時刻である流出時刻３４、車両が道路を通過した時の平均速度３５、そして詳細データ３６で構成される。
【００２６】
プローブデータ記憶装置１２０の詳細データ３６は、一定の時間間隔で開始されるデータの記録の開始時刻３７、車両の進行方向速度３８、車両の進行方向加速度３９、車両の横方向速度４０、車両の横方向加速度４１の各データで構成される。進行方向速度３８は、プローブカーの車速パルスの情報や、位置情報を時間で微分して求められる情報である。進行方向加速度３９は、プローブカーの加速度センサの情報や、速度データを更に時間で微分して求められる情報である。横方向速度４０は、プローブカーのセンサや位置情報を微分して求められる車両の横方向への移動速度の情報である。横方向加速度４１は、プローブカーの加速度センサの情報や、速度データを微分して求められる、車両の横方向への移動速度の際の加速度情報である。ここで横方向とは、地面と平行で、車両の進行方向と垂直に交わる方向である。
【００２７】
次に、道路特性推定装置１３０の詳細を説明する。道路の交通量Ｑ（単位：台／ｈ）は、ある瞬間の一定距離内に存在する車両の台数である交通密度ｋ（単位：台／km）とその一定距離を通過する車両の平均速度ｖ（単位：km／ｈ）から（式１）で計算される。
【００２８】
Ｑ＝ｋ×ｖ                                     …（式１）
この交通密度ｋは計測が難しいため、交通モデルから平均速度ｖを元に求めることにする。平均速度ｖと交通密度ｋの関係には、さまざまなモデルが提案されており、例えばアンダーウッドのモデル式は次の（式２）で表される。
【００２９】
ｖ＝ｖｆ×ｅｘｐ（−ｋ／ｋ０）                 …（式２）
ｅｘｐ（  ）は指数関数を表し、ｖｆは交通密度が０の時の速度である自由速度、ｋ０は、交通量が最大となる時の交通密度である臨界密度を表す。また、グリーンバーグのモデル式は次の（式３）で表される。
【００３０】
ｖ＝ｃ×ｌｏｇ（ｋｊ／ｋ）                     …（式３）
ここでｌｏｇ（  ）は対数、ｃは定数、ｋｊは、速度が０の時の交通密度である飽和密度である。
【００３１】
また、速度ｖと交通密度ｋの関係を指数関数で表現すると、例えばａ、ｂを定数とした（式４）の様に表される。
【００３２】
ｖ＝ａ×ｅｘｐ（ｂ×ｋ）                       …（式４）
このように、平均速度ｖと交通密度ｋの関係は、モデルとそのパラメータで表現できる。このパラメータは、（式２）のアンダーウッドのモデルでは自由速度ｖｆ、臨界密度ｋ０、（式３）のグリーンバーグのモデルでは定数ｃ、飽和密度ｋｊ、またこれを指数関数で表した（式４）では、定数ａ、ｂである。これらのパラメータを道路特性と呼ぶ。道路特性推定装置１３０は、これらのパラメータを推定する。以下では、アンダーウッドのモデル式を使用するものとして説明を行う。だだし、交通密度ｋと平均速度ｖの関係式は、アンダーウッドのモデル式に限定されるものではない。
【００３３】
道路特性推定装置１３０の処理フローを図４に示す。道路特性推定装置１３０では、プローブデータ記憶装置１２０に蓄積した複数のプローブデータを１ヶ月、１年などの定期的な周期で処理する。以下それぞれのステップについて詳細に説明する。
【００３４】
まずステップＳ１０では、プローブデータ記憶装置１２０に記憶されている全ての道路についてステップＳ２０からＳ４０までの処理を繰り返すループが開始される。ステップＳ２０は、交通量記憶装置１１０に処理対象の道路ＩＤの交通量データがあるかを判定する。交通量データがある場合（Ｓ２０：Ｙｅｓ）は、ステップＳ３０へ進む。交通量データが無い場合（Ｓ２０：Ｎｏ）は、ステップＳ４０へ進む。
【００３５】
ステップＳ３０では、交通量記憶装置１１０から読み出した交通量データとプローブデータを使い、交通量データを教示データとして、処理対象の道路ＩＤの道路特性を推定する。この推定処理では、（式１）と（式２）を使って推定した推定交通量と、教示データである交通量との差が最小になる様に、道路特性である（式２）のパラメータを決定する。そこで道路ＩＤ２１が一致する交通量データの開始時刻２２から終了時刻２３までに収集されたプローブデータをプローブデータ記憶装置１２０から取得する。図５に、横軸にプローブデータの平均速度、縦軸に交通密度を示したグラフを示す。このグラフの縦軸である交通密度は、（式１）により、交通量記憶装置１１０から読み出した交通量の値をプローブデータの速度で割ることで計算できる。図５のグラフにおける１つのサンプル点５１は、プローブデータの速度とその速度の時刻（プローブデータ記憶装置１２０の開始時刻）の交通密度をプロットした点である。１つのプローブデータにつき、図５のグラフ上に１つの点がプロットできる。
【００３６】
次にプロットした点に対し、アンダーウッドのモデル式である（式２）を適用して、複数個のサンプル点５１のデータを使って、自由速度ｖｆと臨界密度ｋ０を誤差が最小になる様に設定した曲線を求める。これが図５の曲線５２である。このようにして、処理対象の道路ＩＤについて、モデル式のパラメータである自由速度ｖｆと臨界密度ｋ０を求める。このパラメータがこの道路の道路特性となる。
【００３７】
図４のステップＳ４０では、対象となる道路について教示データとする交通量データが交通量記憶装置１１０に無いため、プローブデータのみを使って道路特性を推定する。例えばこの場合、自由速度ｖｆは、処理対象の道路ＩＤを走行したプローブデータの速度の中で最大値、又は９０パーセンタイルの値とする。しかし、臨界密度ｋ０は、プローブデータのみでは推定することが難しいため、予め定めた値を使うことにする。例えば、日本の場合、臨界密度は１車線あたり４０〜６０台／kmであることが経験的に知られていることから、その中央値となる５０台／kmの値を用いることとする。
【００３８】
また交通密度と平均速度の関係が（式４）の様な指数関数であるとした場合も、道路特性を表す定数ａ、ｂはプローブデータのみから求めることができない。このため、交通密度と平均速度の関係を、（式４）の様な関数で仮定している場合は、ステップＳ４０では道路特性を作成することなく、ステップＳ５０へ進み、ステップＳ４０の処理を終了する。この場合、プローブデータから道路特性が推定できないため、道路特性はステップＳ３０で推定された値だけが使用されることになる。
【００３９】
図４のステップＳ５０は、全ての道路について処理を完了したか否かを判定し、処理を完了した場合は処理フローを終了する。全ての道路について処理を完了していない場合は、ステップＳ２０へ戻る。
【００４０】
次に、分類装置１４０の詳細を説明する。分類装置１４０は、道路特性推定装置１３０で推定した全ての道路の道路特性データをクラスタリングして、クラスタを作成する。図６にクラスタの概念図を示す。図６（ａ）は、自由速度ｖｆと臨界密度ｋ０の２軸のグラフに、この２つのパラメータを道路特性とする道路をマッピングした結果である。１つのサンプル点６１は、自由速度ｖｆと臨界密度ｋ０のデータを表しており、図５の曲線５２のモデル式から求めたパラメータである。つまり、１つの道路につき、１つの道路特性データが計算でき、この道路特性データを図６（ａ）上に１つの点としてプロットする。図６（ｂ）は、これらの点を既存の方法でクラスタ６２にクラスタリングした結果を示している。クラスタリングは、全ての道路の道路特性データの中で、値が近い道路を同じクラスタに分類し、値が離れている道路は異なるクラスタに分類するように計算する。この計算は、クラスタ数が一定数になったら終了する。各クラスタは、属する道路特性の重心位置（例えば平均値）を計算する。この重心位置をクラスタの代表道路特性とする。
【００４１】
分類装置１４０は、２つの道路特性（自由速度ｖｆと臨界密度ｋ０）に限定されるものではない。アンダーウッドのモデル式では、道路特性が自由速度ｖｆと臨界密度ｋ０の２つであるため、道路特性をクラスタリングする際には、道路特性が図６の様な２次元のグラフにプロットされるが、もし３つの道路特性を扱うモデル式を用いる場合、道路特性は３次元のグラフで表されることになる。
【００４２】
分類装置１４０で計算したクラスタに関する情報は、記憶装置１６０に格納される。各クラスタには識別のためのＩＤ（以下、クラスタＩＤ）が付与される。また、クラスタに属する道路特性に対応するプローブデータが収集された道路の道路ＩＤが、記憶装置１６０に格納される。更に、クラスタに属する道路特性に対応するプローブデータを収集したプローブカーの車両ＩＤを、走行パターン作成装置１５０で使用するため、一時的な記憶領域に格納する。
【００４３】
次に、走行パターン作成装置１５０の詳細を説明する。図７と図８に速度パターンの概念図を示す。走行パターンは、進行方向の速度に対する進行方向の加速度の分布と、進行方向の速度に対する横方向の加速度の分布で定義される。以下では、進行方向の速度に対する進行方向の加速度の分布を第１の走行パターン、進行方向の速度に対する横方向の加速度の分布を第２の走行パターンと定義する。ここで横方向とは、前述のように、地面と平行で、車両の進行方向と垂直に交わる方向のことである。
【００４４】
図７（ａ）は、道路幅が広い道路７１（以下、幅広道路）で、プローブカー７２が、自転車やバイク、あるいは自分より低速で走行する車両などの他の車両を回避しながら走行していく状況を表している。図７（ｂ）は、この幅広道路における第１の走行パターン、図７（ｃ）は第２の走行パターンを表したものであり、横軸が車両の進行方向の速度、縦軸が車両の進行方向の加速度を示している。図７（ｂ）では走行パターンの説明を分かりやすくするために、進行方向の速度に対する加速度の分布を、頻度の高い領域７３（高頻度領域）と、頻度の低い領域７４（低頻度領域）に分けて表している。
【００４５】
図７に示した幅広道路の様な場合には、車両が横方向に自由に移動でき、他の車両を回避しながら走行することができる。このため、図７（ｂ）の様に、進行方向の加速度は、進行方向の速度の高速度域では、幅広く分布するのに対し、進行方向の速度の低速度域では、渋滞しているため、狭く分布する。一方、他の車両を回避するため車両は横方向に比較的大きく移動することから、図７（ｃ）の様に横方向の加速度は、進行方向の速度の速度域に寄らず、幅広く分布する。
【００４６】
図８（ａ）は、道路幅が狭い道路８１（以下、幅狭道路）で、プローブカー８２が他の車両を回避できず、先行する車両と連なって走行していく状況を表している。図８（ｂ）は、この幅狭道路における第１の走行パターン、図８（ｃ）は、第２の走行パターンを表したものであり、図７と同様、進行方向の速度に対する加速度の分布を、頻度の高い領域８３（高頻度領域）と、頻度の低い領域８４（低頻度領域）に分けて表している。図８に示す様な幅狭道路では、車両が横方向に自由に移動できないため、進行方向の加速度と横方向の加速度は、進行方向の速度の速度域に寄らず、狭く分布する。
【００４７】
このように道路によって走行パターンが異なり、走行パターン作成装置１５０では、分類装置１４０で分類したクラスタ毎に、第１の走行パターンと第２の走行パターンを作成する。
【００４８】
走行パターン作成装置１５０の処理フローを図９に示す。走行パターン作成装置１５０では、ステップＳ１００からステップＳ１５０までのループ処理（ステップＳ１１０〜ステップＳ１４０）の中で、分類装置１４０でクラスタリングした全クラスタを処理する。そしてステップＳ１１０からＳ１３０までのループ処理によって、処理対象のクラスタを構成する全てのプローブデータをステップＳ１２０で処理して走行パターンを作成する。この処理対象のクラスタを構成する生成した際に、同じクラスタに属する道路特性に対応するプローブデータを収集したプローブカーの車両ＩＤの情報は、前述の通り、分類装置１４０で作成し、一時記憶領域に保存されている。また処理対象のクラスタに属する道路特性に対応したプローブデータが収集された道路の道路ＩＤの情報は、記憶装置１６０に保存されている。プローブデータ記憶装置１２０から、これらの車両ＩＤと道路ＩＤに該当するプローブデータを抽出する。そして抽出したプローブデータに対してステップＳ１２０の処理を繰り返すことで走行パターンを生成する。
【００４９】
ステップＳ１２０では、処理対象とするプローブデータの詳細データ３６を取得して、各プローブデータの第１の走行パターンと第２の走行パターンの各走行パターンを作成する。各走行パターンは、進行方向速度に対する進行方向加速度／横方向加速度の関係を表すデータであり、図１０に示す様なテーブルを用いて作成する。図１０のテーブルは、速度帯域毎に分割された列１０１と、加速度域毎に区分された行１０２からなる。このテーブルの列１０１では、例えば「０〜１０」の列が、０〜１０km／ｈの速度域を表している。また行１０２では、例えば「−３〜３」の行は、−３〜３ｍ／ｓ2の加速度域を表している。そしてテーブルのセル１０３には、始めは初期値として全て‘０’が入力されているものとする。ステップＳ１２０では、各プローブデータの詳細データ３６から、その進行方向速度３８の値と進行方向加速度３９及び横方向加速度４１の値に該当する速度域、加速度域のセルを見つけ、該当するセルの値に１加える処理を繰り返す。これによりセルの値が大きい程、頻度の高い速度域・加速度域であることを表すことになる。
【００５０】
ステップＳ１３０は、処理対象のクラスタを構成するプローブデータを全て処理したか否かを判定し、１つのクラスタにおける全てのプローブデータについて処理が完了したらステップＳ１４０へ進む。また処理を完了していない場合には、ステップＳ１２０を繰り返す。
【００５１】
次にステップＳ１４０では、ステップＳ１２０で作成した処理対象のクラスタの走行パターンを正規化する。クラスタを構成するプローブデータの数は、クラスタによって異なるため、正規化することで、全てのクラスタを比較できるようにする。図１０に示した走行パターンのテーブルにおけるｉ行目ｊ列目のセルの値をＳｉｊとする。そしてこのセルの正規化後の値をＲｉｊとした時、Ｒｉｊは（式５）の様に表される。
【００５２】
Ｒｉｊ＝Ｓｉｊ／ΣΣＳｉｊ                     …（式５）
この（式５）の分母ΣΣＳｉｊは、全てのセルの値の総和である。このように全てのセルの総和が１となるように、セルの値を正規化する。正規化した第１の走行パターンと第２の走行パターンは、対応するクラスタＩＤと関連付けられて記憶装置１６０に格納される。
【００５３】
そして最後にステップＳ１５０で、分類装置１４０によりクラスタリングした全クラスタを処理したか否かを判定する。全てのクラスタについて処理を完了した場合は、走行パターン作成装置１５０の処理を終了する。しかしまだ処理を完了していないクラスタが残っている場合には、ステップＳ１００へ戻る。
【００５４】
次に、記憶装置１６０の構成を説明する。図１１は、記憶装置１６０のデータフォーマットを示した図である。記憶装置１６０のデータは、分類装置１４０で分類したクラスタＩＤ１１１、クラスタの代表道路特性１１２、クラスタを構成する道路ＩＤ１１３、第１の走行パターン１１４、第２の走行パターン１１５で構成される。代表道路特性１１２は、道路特性推定装置１３０でアンダーウッドのモデル式を使った場合は、自由速度ｖｆと臨界密度ｋ０で構成されているが、他のモデル式を使用する場合には、異なる道路特性が格納される。第１の走行パターン１１４、第２の走行パターン１１５には、それぞれ図１０と同じデータフォーマットで記録され、走行パターン作成装置１５０で作成された第１の走行パターン、第２の走行パターンが格納される。
【００５５】
次に、道路特性検索装置１８０を説明する。道路特性検索装置１８０は、受信装置１７０で受信したプローブデータを使って、そのプローブデータを収集した道路の道路特性に近いクラスタを記憶装置１６０から検索する。
【００５６】
この道路特性検索装置１８０の処理フローを図１２に示す。一定期間内に受信装置１７０で受信したプローブデータが複数個ある場合は、プローブデータの数だけ図１２の処理を実行する。図１２のステップＳ２００では、受信装置１７０からのプローブデータに含まれる走行した道路の道路ＩＤと一致する道路ＩＤが、記憶装置１６０から道路ＩＤ１１３の道路ＩＤリストの中に含まれているデータがあるか否かを判定する。一致するものがある場合、一致する道路ＩＤが道路ＩＤ１１３の道路ＩＤリストの中に含まれているクラスタの代表道路特性１１２の情報を取得して、処理を終了する。一致しない場合、プローブデータの走行パターンからクラスタを検索するために、ステップＳ２１０へ進む。
【００５７】
ステップＳ２１０では、受信装置１７０からのプローブデータの第１の走行パターンと第２の走行パターンを作成する。具体的な処理内容は、図９のステップＳ１２０と同じである。次にステップＳ２２０では、ステップＳ２１０で作成したプローブデータの第１の走行パターンと第２の走行パターンを正規化する。具体的な処理内容は、図９のステップＳ１４０と同じである。正規化することで、記憶装置１６０に格納されている各クラスタの走行パターンと比較することができる。
【００５８】
ステップＳ２３０では、ステップＳ２２０で作成したプローブデータから作成した走行パターンと、記憶装置１６０の各クラスタの走行パターン１１４から類似性の高い走行パターンを抽出し、そのクラスタの道路特性を取得する。第１の走行パターンは、（式６）に示す式で類似性を評価する。プローブデータの正規化した第１の走行パターンでｉ行目ｊ列目のセルの値をＰ１ｉｊ、クラスタＩＤ「ｋ」の第１の走行パターンでｉ行目ｊ列目のセルの値をＲ１ｋｉｊとする。この時クラスタＩＤ「ｋ」の第１の走行パターンの評価値をＷ１ｋとする。
【００５９】
Ｗ１ｋ＝ΣΣ｜Ｐ１ｉｊ−Ｒ１ｋｉｊ｜2          …（式６）
また、第２の走行パターンは、（式７）に示す式で類似性を評価する。プローブデータの正規化した第２の走行パターンでｉ行目ｊ列目のセルの値をＰ２ｉｊ、クラスタＩＤ「ｋ」の第２の走行パターンでｉ行目ｊ列目のセルの値をＲ２ｋｉｊとする。この時クラスタＩＤ「ｋ」の第２の走行パターンの評価値をＷ２ｋとする。
【００６０】
Ｗ２ｋ＝ΣΣ｜Ｐ２ｉｊ−Ｒ２ｋｉｊ｜2          …（式７）
この時評価値の和（Ｗ１ｋ＋Ｗ２ｋ）が最も小さいクラスタを、プローブデータの走行した道路が属するクラスタとして、その代表道路特性１１２の情報を取得する。
【００６１】
次に、交通量予測装置１９０を説明する。交通量予測装置１９０は、受信装置１７０のプローブデータの平均速度３５と、道路特性検索装置１８０から取得した代表道路特性１１２の情報を使って、道路単位に交通量を予測する。この時、道路特性推定装置１３０で設定したモデルを使用する。例えばアンダーウッドのモデル式を使用した場合、プローブデータの平均速度ｖと、代表道路特性１１２の道路特性データ（自由速度ｖｆ、臨界密度ｋ０）を使用して、（式２）から、交通密度ｋを推定する。次に、推定した交通密度ｋとプローブデータの平均速度ｖを使用して、（式１）から交通量を予測する。予測した交通量データと道路ＩＤの情報は送信装置２００へ提供する。
【００６２】
以上で説明した実施の形態によれば、交通量予測システム１００は、プローブデータの走行パターンを使用して、類似したグループの道路特性を使用できるため、大きな費用をかけずに、高精度に交通量を予測することができる様になる。
【００６３】
次に、本発明の実施の形態の変形に係る交通量予測システムについて説明する。この実施形態の変形では、交通インフラのセンサなどから収集した道路単位の速度データを受信し、この速度データを用いて道路特性を検索し、交通量を予測する。
【００６４】
この実施の形態の変形に係る交通量予測システムの構成を図１３に示す。図１３の交通量予測システム３００は、交通量記憶装置１１０、プローブデータ記憶装置１２０、道路特性推定装置１３０、分類装置１４０、走行パターン作成装置１５０、記憶装置１６０、受信装置２１０、道路情報記憶装置２２０、道路特性検索装置２３０、交通量予測装置１９０、送信装置２００で構成される。
【００６５】
受信装置２１０は、外部のサーバーからプローブデータの他に道路単位の速度データを受信する装置である。受信する速度データのフォーマットの例を図１４に示す。速度データは、データを収集した道路の道路ＩＤ１４１、データを収集した時刻１４２、そして収集した速度１４３の各データからなる。受信した速度データは、道路特性検索装置２３０と交通量予測装置１９０へ提供される。
道路情報記憶装置２２０は、ハードディスク、フラッシュメモリ等の記憶装置で、道路の属性情報を記憶しており、道路特性検索装置２３０から参照される。道路情報記憶装置２２０のデータフォーマットの例を図１５に示す。道路情報記憶装置２２０には道路の属性情報として、道路ＩＤ１５１に対応した道路の、車線数１５２、道路幅１５３が格納される。なお、道路の属性情報とは、車線数、道路幅など道路に関する詳細情報であるが、車線数、道路幅に限定されるものではない。
【００６６】
道路特性検索装置２３０は、交通量予測システム３００に搭載されている不図示のマイクロプロセッサや、ＲＡＭ、ＲＯＭなどによって実行されるソフトウェアにより実現される。道路特性検索装置２３０は、受信装置２１０からの速度データを収集した道路ＩＤと、記憶装置１６０に格納している各クラスタを構成する道路ＩＤと、道路情報記憶装置２２０の道路情報を利用して、速度データの道路ＩＤと類似するクラスタを検索する。
【００６７】
道路特性検索装置２３０の処理フローを図１６に示す。まずステップＳ３００では、受信装置２１０からの速度データに含まれている道路ＩＤと一致する道路ＩＤが、道路ＩＤ１１３に記録されている道路ＩＤリストの中に含まれるデータが記憶装置１６０にあるか否かを判定する。一致する道路ＩＤを含む道路ＩＤリストを持ったデータが存在する場合（Ｓ３００：Ｙｅｓ）、そのデータのクラスタＩＤ１１１と代表道路特性１１２の情報を取得して、処理を終了する。一致する道路ＩＤを含む道路ＩＤリストを持ったデータが存在しない場合（Ｓ３００：Ｎｏ）、道路の属性情報から類似する道路に対応したクラスタを検索するために、ステップＳ３１０へ進む。ステップＳ３１０では、受信装置２１０からの速度データの道路ＩＤと一致する道路ＩＤの道路属性情報を道路情報記憶装置２２０から取得する。
【００６８】
ステップＳ３２０では、記憶装置１６０の全クラスタの中から類似性の高いクラスタを抽出し、そのクラスタの道路特性を取得する。具体的には、始めに記憶装置１６０の全てのクラスタについて、各クラスタの道路ＩＤ１１３の道路ＩＤリストに含まれる道路ＩＤを使って、道路情報記憶装置２２０から、一致する道路ＩＤの道路属性情報を取得する。次に、クラスタに含まれる道路の中で、受信装置２１０からの速度データの道路の道路属性と類似性の高い道路を検索して、類似性の最も高いクラスタを見つける。最後に、類似性が高いと判定されたクラスタの道路特性を記憶装置１６０から取得する。
【００６９】
道路属性に関するクラスタの類似性の評価方法を説明する。ここでは、受信装置２１０からの速度データの道路の道路属性情報を、図１５に示す様な車線数と道路幅とする。クラスタを構成するｉ番目の道路の車線数をＡ（ｉ）、道路幅をＢ（ｉ）、評価処理対象の道路属性情報の車線数をＣ、道路幅をＤとした時のクラスタの評価値Ｅを、以下の（式８）により計算する。
【００７０】
Ｅ＝Σ（｜Ａ(ｉ)−Ｃ｜2＋｜Ｂ(ｉ)−Ｄ｜2）     …（式８）
（式８）では、クラスタを構成する全道路について、評価対象道路の車線数と道路幅を比較している。評価対象道路の車線数と道路幅の差分を、クラスタを構成する道路数足しているため、評価値Ｅが小さいクラスタほど、類似性が高くなる。また、道路属性にバスレーンの有無など、数値化できない情報が含まれる場合は、道路属性が同じ場合は１、異なる場合は０として計算する。
【００７１】
これにより、交通量予測システム３００では、プローブカーが走行した道路の属性情報を使用して、類似したグループの道路特性を使用できるため、大きな費用をかけずに、高精度に交通量を予測できる様になる。
【符号の説明】
【００７２】
１００、３００  交通量予測システム
１１０  交通量記憶装置
１２０  プローブデータ記憶装置
１３０  道路特性推定装置
１４０  分類装置
１５０  走行パターン作成装置
１６０  記憶装置
１７０、２１０  受信装置
１８０、２３０  道路特性検索装置
１９０  交通量予測装置
２００  送信装置
２２０  道路情報記憶装置
【図１】
2013128486.tif 000003
【図２】
2013128486.tif 000004
【図３】
2013128486.tif 000005
【図４】
2013128486.tif 000006
【図５】
2013128486.tif 000007
【図６】
2013128486.tif 000008
【図７】
2013128486.tif 000009
【図８】
2013128486.tif 000010
【図９】
2013128486.tif 000011
【図１０】
2013128486.tif 000012
【図１１】
2013128486.tif 000013
【図１２】
2013128486.tif 000014
【図１３】
2013128486.tif 000015
【図１４】
2013128486.tif 000016
【図１５】
2013128486.tif 000017
【図１６】
2013128486.tif 000018
【国際調査報告】
2013128486.tif 000019
2013128486.tif 000020
2013128486.tif 000021
2013128486.tif 000022