知財求人 - 知財ポータルサイト「IP Force」
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2024170260
(43)【公開日】2024-12-06
(54)【発明の名称】ＯＤ決定方法およびＯＤ決定プログラム
(51)【国際特許分類】
G08G 1/00 20060101AFI20241129BHJP
【ＦＩ】
G08G1/00 C
【審査請求】未請求
【請求項の数】5
【出願形態】ＯＬ
(21)【出願番号】P 2023087326
(22)【出願日】2023-05-26
(71)【出願人】
【識別番号】000005223
【氏名又は名称】富士通株式会社
(74)【代理人】
【識別番号】100104190
【弁理士】
【氏名又は名称】酒井 昭徳
(72)【発明者】
【氏名】陳 彬
(72)【発明者】
【氏名】瀬川 英吾
【テーマコード（参考）】
5H181
【Ｆターム（参考）】
5H181AA01
5H181BB04
5H181DD04
5H181EE03
(57)【要約】
【課題】ＯＤ交通量の推定にかかる処理時間の短縮化を図ること。
【解決手段】情報処理装置１０１は、第１のＯＤ交通量ｑ１を用いて、第１の時間フレームごとに、各リンクの交通量を推定する第１のシミュレーション１１０を実行する。情報処理装置１０１は、第１のＯＤ交通量ｑ１を調整しながら第１のシミュレーション１１０を繰り返して、第１の時間フレームごとのＯＤ交通量を決定する。情報処理装置１０１は、決定した第１の時間フレームごとのＯＤ交通量に基づいて、第１の時間フレームを分割した第２の時間フレームごとに、第２のＯＤ交通量ｑ２を算出する。情報処理装置１０１は、算出した第２のＯＤ交通量ｑ２を用いて、各リンクの交通量を推定する第２のシミュレーション１２０を実行する。情報処理装置１０１は、第２のＯＤ交通量ｑ２を調整しながら第２のシミュレーション１２０を繰り返して、第２の時間フレームごとのＯＤ交通量を決定する。
【選択図】図１
【解決手段】情報処理装置１０１は、第１のＯＤ交通量ｑ１を用いて、第１の時間フレームごとに、各リンクの交通量を推定する第１のシミュレーション１１０を実行する。情報処理装置１０１は、第１のＯＤ交通量ｑ１を調整しながら第１のシミュレーション１１０を繰り返して、第１の時間フレームごとのＯＤ交通量を決定する。情報処理装置１０１は、決定した第１の時間フレームごとのＯＤ交通量に基づいて、第１の時間フレームを分割した第２の時間フレームごとに、第２のＯＤ交通量ｑ２を算出する。情報処理装置１０１は、算出した第２のＯＤ交通量ｑ２を用いて、各リンクの交通量を推定する第２のシミュレーション１２０を実行する。情報処理装置１０１は、第２のＯＤ交通量ｑ２を調整しながら第２のシミュレーション１２０を繰り返して、第２の時間フレームごとのＯＤ交通量を決定する。
【選択図】図１
【特許請求の範囲】
【請求項1】
対象期間を分割した第１の時間フレームごとに、当該第１の時間フレームに第１の地点を出発して第２の地点に到達する第１のＯＤ交通量を用いて、前記第１の地点から前記第２の地点までの経路に含まれるリンクごとの交通量を推定する第１の交通シミュレーションを実行し、
推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第１のＯＤ交通量を調整しながら前記第１の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第１の時間フレームごとのＯＤ交通量を決定し、
決定した前記第１の時間フレームごとのＯＤ交通量に基づいて、前記第１の時間フレームを分割した第２の時間フレームごとに、当該第２の時間フレームに前記第１の地点を出発して前記第２の地点に到達する第２のＯＤ交通量を算出し、
算出した前記第２のＯＤ交通量を用いて、前記第２の時間フレームごとに、前記リンクごとの交通量を推定する第２の交通シミュレーションを実行し、
推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第２のＯＤ交通量を調整しながら前記第２の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第２の時間フレームごとのＯＤ交通量を決定する、
処理をコンピュータが実行することを特徴とするＯＤ決定方法。
推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第１のＯＤ交通量を調整しながら前記第１の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第１の時間フレームごとのＯＤ交通量を決定し、
決定した前記第１の時間フレームごとのＯＤ交通量に基づいて、前記第１の時間フレームを分割した第２の時間フレームごとに、当該第２の時間フレームに前記第１の地点を出発して前記第２の地点に到達する第２のＯＤ交通量を算出し、
算出した前記第２のＯＤ交通量を用いて、前記第２の時間フレームごとに、前記リンクごとの交通量を推定する第２の交通シミュレーションを実行し、
推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第２のＯＤ交通量を調整しながら前記第２の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第２の時間フレームごとのＯＤ交通量を決定する、
処理をコンピュータが実行することを特徴とするＯＤ決定方法。
【請求項2】
前記算出する処理は、
決定した前記第１の時間フレームごとのＯＤ交通量を、所定の比率に応じて前記第２の時間フレームごとに配分することにより、前記第２のＯＤ交通量を算出する、ことを特徴とする請求項１に記載のＯＤ決定方法。
決定した前記第１の時間フレームごとのＯＤ交通量を、所定の比率に応じて前記第２の時間フレームごとに配分することにより、前記第２のＯＤ交通量を算出する、ことを特徴とする請求項１に記載のＯＤ決定方法。
【請求項3】
前記第１の地点と前記第２の地点とのペアと対応付けて、決定した前記第２の時間フレームごとのＯＤ交通量を出力する、
処理を前記コンピュータが実行することを特徴とする請求項１に記載のＯＤ決定方法。
処理を前記コンピュータが実行することを特徴とする請求項１に記載のＯＤ決定方法。
【請求項4】
前記比率は、前記第１の時間フレームと前記第２の時間フレームとの時間幅をもとに設定される、ことを特徴とする請求項２に記載のＯＤ決定方法。
【請求項5】
対象期間を分割した第１の時間フレームごとに、当該第１の時間フレームに第１の地点を出発して第２の地点に到達する第１のＯＤ交通量を用いて、前記第１の地点から前記第２の地点までの経路に含まれるリンクごとの交通量を推定する第１の交通シミュレーションを実行し、
推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第１のＯＤ交通量を調整しながら前記第１の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第１の時間フレームごとのＯＤ交通量を決定し、
決定した前記第１の時間フレームごとのＯＤ交通量に基づいて、前記第１の時間フレームを分割した第２の時間フレームごとに、当該第２の時間フレームに前記第１の地点を出発して前記第２の地点に到達する第２のＯＤ交通量を算出し、
算出した前記第２のＯＤ交通量を用いて、前記第２の時間フレームごとに、前記リンクごとの交通量を推定する第２の交通シミュレーションを実行し、
推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第２のＯＤ交通量を調整しながら前記第２の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第２の時間フレームごとのＯＤ交通量を決定する、
処理をコンピュータに実行させることを特徴とするＯＤ決定プログラム。
推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第１のＯＤ交通量を調整しながら前記第１の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第１の時間フレームごとのＯＤ交通量を決定し、
決定した前記第１の時間フレームごとのＯＤ交通量に基づいて、前記第１の時間フレームを分割した第２の時間フレームごとに、当該第２の時間フレームに前記第１の地点を出発して前記第２の地点に到達する第２のＯＤ交通量を算出し、
算出した前記第２のＯＤ交通量を用いて、前記第２の時間フレームごとに、前記リンクごとの交通量を推定する第２の交通シミュレーションを実行し、
推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第２のＯＤ交通量を調整しながら前記第２の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第２の時間フレームごとのＯＤ交通量を決定する、
処理をコンピュータに実行させることを特徴とするＯＤ決定プログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、ＯＤ決定方法およびＯＤ決定プログラムに関する。
【背景技術】
【0002】
近年、自治体や企業などにおいて、都市部の交通渋滞の軽減や交通によるＣＯ２排出量の削減などを実現するために、様々な交通施策が検討されている。例えば、様々な交通状況をデジタルツイン上で再現して、交通移動需要の推定を行って、施策効果を評価することが行われている。交通移動需要は、例えば、Ｏｒｉｇｉｎ（出発地点）を出発してＤｅｓｔｉｎａｔｉｏｎ（到着地）に到達する車両数などのＯＤ交通量によって表される。
【0003】
先行技術としては、プローブ情報や統計情報に基づいてＯＤ交通量を推定し、道路ネットワーク情報、信号情報、推定されたＯＤ交通量に基づいて、交通状態をシミュレーションし、シミュレーションされた交通状態に基づいて算出されるリンクの交通量の計算値と、リンクの交通量の目標値との差が小さくなるように、ＯＤ交通量を修正し、修正されたＯＤ交通量を用いて、交通状態をシミュレーションするものがある。
【先行技術文献】
【特許文献】
【0004】
【特許文献1】特開２００９－２５９１５８号公報
【発明の概要】
【発明が解決しようとする課題】
【0005】
しかしながら、従来技術では、ＯＤ交通量（交通移動需要）の推定に膨大な処理時間がかかる場合がある。
【0006】
一つの側面では、本発明は、ＯＤ交通量の推定にかかる処理時間の短縮化を図ることを目的とする。
【課題を解決するための手段】
【0007】
１つの実施態様では、対象期間を分割した第１の時間フレームごとに、当該第１の時間フレームに第１の地点を出発して第２の地点に到達する第１のＯＤ交通量を用いて、前記第１の地点から前記第２の地点までの経路に含まれるリンクごとの交通量を推定する第１の交通シミュレーションを実行し、推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第１のＯＤ交通量を調整しながら前記第１の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第１の時間フレームごとのＯＤ交通量を決定し、決定した前記第１の時間フレームごとのＯＤ交通量に基づいて、前記第１の時間フレームを分割した第２の時間フレームごとに、当該第２の時間フレームに前記第１の地点を出発して前記第２の地点に到達する第２のＯＤ交通量を算出し、算出した前記第２のＯＤ交通量を用いて、前記第２の時間フレームごとに、前記リンクごとの交通量を推定する第２の交通シミュレーションを実行し、推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第２のＯＤ交通量を調整しながら前記第２の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第２の時間フレームごとのＯＤ交通量を決定する、ＯＤ決定方法が提供される。
【発明の効果】
【0008】
本発明の一側面によれば、ＯＤ交通量の推定にかかる処理時間の短縮化を図ることができるという効果を奏する。
【図面の簡単な説明】
【0009】
【発明を実施するための形態】
【0010】
以下に図面を参照して、本発明にかかるＯＤ決定方法およびＯＤ決定プログラムの実施の形態を詳細に説明する。
【0011】
（実施の形態）
図１は、実施の形態にかかるＯＤ決定方法の一実施例を示す説明図である。図１において、情報処理装置１０１は、出発地点と目的地点とのペアについて、指定された時間フレームごとのＯＤ（Ｏｒｉｇｉｎ Ｄｅｓｔｉｎａｔｉｏｎ）交通量を決定するコンピュータである。
図１は、実施の形態にかかるＯＤ決定方法の一実施例を示す説明図である。図１において、情報処理装置１０１は、出発地点と目的地点とのペアについて、指定された時間フレームごとのＯＤ（Ｏｒｉｇｉｎ Ｄｅｓｔｉｎａｔｉｏｎ）交通量を決定するコンピュータである。
【0012】
ここで、ＯＤ交通量とは、交通移動需要を表す指標であり、ある時間フレームに出発地点を出発して目的地点に到達する交通量を表す。ＯＤ交通量は、例えば、交通シミュレーションを利用して推定される。時間フレームは、ある期間を等間隔に分割した時間区間（タイムスロット）である。交通量は、例えば、車両数や人数によって表される。
【0013】
既存の交通需要推定技術として、例えば、交通シミュレーションを利用して、時間フレームごとのＯＤ交通量を推定する技術（従来技術１）がある。交通シミュレーションは、例えば、道路上で起こる交通現象をコンピュータ上で再現し、交通渋滞や交通の乱れなどの交通状況がどのように変化するのかを模擬することである。
【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】
（ＯＤペア表の具体例）
つぎに、ＯＤ推定装置２０１が用いるＯＤペア表の具体例について説明する。
つぎに、ＯＤ推定装置２０１が用いるＯＤペア表の具体例について説明する。
【0059】
図８は、ＯＤペア表の具体例を示す説明図である。図８において、ＯＤペア表８００は、道路ネットワークＮＷ１におけるＯＤペアを示す。ＯＤペアは、Ｏ（Ｏｒｉｇｉｎ）とＤ（Ｄｅｓｔｉｎａｔｉｏｎ）との組み合わせである。図８中の＃は、ＯＤペアを識別する識別子である。
【0060】
ＯＤペアは、例えば、図６に示した出発地点表６００と、図７に示した目的地点表７００とから特定される。例えば、ＯＤペア１は、Ｏ：ノードＮ１（ｎｏｄｅ１）とＤ：ノードＮ４（ｎｏｄｅ４）との組み合わせである。
【0061】
（経路表の具体例）
つぎに、ＯＤ推定装置２０１が用いる経路表の具体例について説明する。
つぎに、ＯＤ推定装置２０１が用いる経路表の具体例について説明する。
【0062】
図９は、経路表の具体例を示す説明図である。図９において、経路表９００は、道路ネットワークＮＷ１におけるＯＤペアに対応する経路（ｐａｔｈ）を示す。経路（ｐａｔｈ）は、Ｏ（Ｏｒｉｇｉｎ）からＤ（Ｄｅｓｔｉｎａｔｉｏｎ）までの経路である。図９中の＃は、経路（ｐａｔｈ）を識別する識別子である。
【0063】
例えば、ｐａｔｈ１は、Ｏ：ノードＮ１（ｎｏｄｅ１）からＤ：ノードＮ４（ｎｏｄｅ４）までの経路「ｎｏｄｅ１－ｌｉｎｋ１－ｌｉｎｋ３－ｎｏｄｅ４」である。なお、ｐａｔｈでは、経由地点となるノード（例えば、ｎｏｄｅ３）の表示は省略されている。
【0064】
（ＯＤ推定装置２０１の機能的構成例）
つぎに、ＯＤ推定装置２０１の機能的構成例について説明する。
つぎに、ＯＤ推定装置２０１の機能的構成例について説明する。
【0065】
図１０は、ＯＤ推定装置２０１の機能的構成例を示すブロック図である。図１０において、ＯＤ推定装置２０１は、取得部１００１と、設定部１００２と、第１の実行部１００３と、第１の決定部１００４と、算出部１００５と、第２の実行部１００６と、第２の決定部１００７と、出力部１００８と、を含む。取得部１００１～出力部１００８は制御部１０００となる機能であり、具体的には、例えば、図３に示したメモリ３０２、ディスク３０４、可搬型記録媒体３０７などの記憶装置に記憶されたプログラムをＣＰＵ３０１に実行させることにより、または、通信Ｉ／Ｆ３０５により、その機能を実現する。各機能部の処理結果は、例えば、メモリ３０２、ディスク３０４などの記憶装置に記憶される。
【0066】
取得部１００１は、道路ネットワーク表を取得する。ここで、道路ネットワーク表は、交通需要の推定対象となる道路（例えば、図４に示した道路ネットワークＮＷ１）を表す情報である。具体的には、例えば、取得部１００１は、図２に示したクライアント装置２０２から受信することにより、図５に示したような道路ネットワーク表５００を取得する。また、取得部１００１は、不図示の入力装置を用いたユーザの操作入力により、道路ネットワーク表５００を取得してもよい。
【0067】
また、取得部１００１は、出発地点表および目的地点表を取得する。ここで、出発地点表は、出発地点を表す情報であり、例えば、道路ネットワーク内のＯ（Ｏｒｉｇｉｎ）となるノードを表す。目的地点表は、目的地点を表す情報であり、例えば、道路ネットワーク内のＤ（Ｄｅｓｔｉｎａｔｉｏｎ）となるノードを表す。
【0068】
具体的には、例えば、取得部１００１は、クライアント装置２０２から受信することにより、図６に示したような出発地点表６００および図７に示したような目的地点表７００を取得する。また、取得部１００１は、不図示の入力装置を用いたユーザの操作入力により、出発地点表６００および目的地点表７００を取得してもよい。
【0069】
また、取得部１００１は、ＯＤペア表を取得する。ここで、ＯＤペア表は、Ｏ（Ｏｒｉｇｉｎ）とＤ（Ｄｅｓｔｉｎａｔｉｏｎ）との組み合わせ（ＯＤペア）を表す情報である。具体的には、例えば、取得部１００１は、出発地点表６００および目的地点表７００から作成することにより、図８に示したようなＯＤペア表８００を取得する。
【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】
また、所定の比率は、例えば、出発地点（第１の地点）から目的地点（第２の地点）までの経路に含まれる第１のリンクについて観測された第２の時間フレームごとの交通量の比をもとに設定されてもよい。例えば、第１のリンクは、例えば、出発地点（第１の地点）のノードに接続されたリンクである。
【0102】
第２のＯＤフローの算出例については、図１１を用いて後述する。
【0103】
第２の実行部１００６は、算出された第２のＯＤフローを用いて、第２の交通シミュレーションを実行する。第２の交通シミュレーションは、第１の時間フレームを分割した第２の時間フレームごとに、出発地点（第１の地点）から目的地点（第２の地点）までの経路に含まれるリンクごとの交通量を推定する交通シミュレーションである。
【0104】
具体的には、例えば、第２の実行部１００６は、道路ネットワーク表（道路ネットワーク）、第２のＯＤフロー、経路表および経路選択表を、交通シミュレータに与えることにより、第２の交通シミュレーションを実行する。これにより、第２の時間フレームごとに、出発地点（第１の地点）から目的地点（第２の地点）までの経路に含まれるリンクごとの交通量を表すリンクフローが得られる。
【0105】
第２の決定部１００７は、推定されたリンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、第２のＯＤフローを調整しながら、第２の交通シミュレーションを繰り返し実行する。これにより、第１の決定部１００４は、ＯＤペアについて、第２の時間フレームごとのＯＤフローを決定する。
【0106】
具体的には、例えば、第２の決定部１００７は、推定されたリンクの交通量（リンクフロー）と当該リンクについて観測された交通量との誤差を算出する。より詳細に説明すると、例えば、第２の決定部１００７は、観測行列を用いて、道路ネットワークから交通量を観測可能なリンクを抽出する。
【0107】
つぎに、第２の決定部１００７は、第２の時間フレームごとに、推定された各リンクのリンクフローのうち、抽出したリンクのリンクフローを抽出する。そして、第２の決定部１００７は、第２の時間フレームごとに、抽出したリンクのリンクフロー（推定リンクフロー）と当該リンクについて観測されたリンクフロー（実測リンクフロー）との誤差を算出する。
【0108】
つぎに、第２の決定部１００７は、算出した誤差が許容範囲内であるか否かを判断する。ここで、算出した誤差が許容範囲内ではない場合、第２の決定部１００７は、算出した誤差に基づいて、第２のＯＤフローを更新する。より詳細に説明すると、例えば、第２の決定部１００７は、推定されたリンクのリンクフローと、算出した誤差とに基づいて、第２のＯＤフローの補正値を算出する。つぎに、第２の決定部１００７は、算出した補正値を用いて、第２のＯＤフローを更新する。
【0109】
なお、第２のＯＤフローの更新例については、図３０を用いて後述する。
【0110】
そして、第２の決定部１００７は、更新した第２のＯＤフローを用いて、第２の交通シミュレーションを実行する。また、第２の決定部１００７は、算出した誤差が許容範囲内の場合、直前の第２の交通シミュレーションに用いた第２のＯＤフローを、第２の時間フレームごとのＯＤフローに決定する。
【0111】
決定された第２の時間フレームごとのＯＤフローの具体例については、図３１を用いて後述する。
【0112】
出力部１００８は、ＯＤペアと対応付けて、決定した第２の時間フレームごとのＯＤフローを出力する。ＯＤペアは、出発地点（第１の地点）と目的地点（第２の地点）との組み合わせである。出力部１００８の出力形式としては、例えば、メモリ３０２、ディスク３０４などの記憶装置への記憶、通信Ｉ／Ｆ３０５による他のコンピュータへの送信、不図示のディスプレイへの表示、不図示のプリンタへの印刷出力などがある。
【0113】
具体的には、例えば、出力部１００８は、図３１に示すようなＯＤフロー２２０２（ＯＤ推定結果）をクライアント装置２０２に送信する。これにより、クライアント装置２０２において、ユーザは、ＯＤペア（出発地点、目的地点）について、第２の時間フレームごとのＯＤフロー（例えば、車両数）を参照することができる。
【0114】
なお、第２の時間フレームよりもさらに細かい時間幅のＯＤフローを推定したい場合がある。そこで、第２の決定部１００７は、第２の時間フレームが、あらかじめ指定された時間フレームであるか否かを判断してもよい。時間フレームの指定は、例えば、クライアント装置２０２から受け付ける。
【0115】
例えば、第２の時間フレームが３０分間隔の時間フレームであるとする。また、指定された時間フレームが１５分間隔の時間フレームであるとする。この場合、第２の決定部１００７は、第２の時間フレームが、指定された時間フレームではないと判断する。算出部１００５は、指定された時間フレームではないと判断された場合、決定された第２の時間フレームごとのＯＤフローを、第１の時間フレームごとのＯＤフローとする。
【0116】
そして、算出部１００５は、当該第１の時間フレームを分割した第２の時間フレームごとに、第２のＯＤフローを算出してもよい。これにより、ＯＤ推定装置２０１は、例えば、段階的に時間フレームを細かくしながら、ＯＤフローの推定を繰り返して、指定された時間フレームごとのＯＤフローを推定することができる。
【0117】
（第２のＯＤフローの算出例）
つぎに、図１１を用いて、第２のＯＤフローの算出例について説明する。
つぎに、図１１を用いて、第２のＯＤフローの算出例について説明する。
【0118】
図１１は、第２のＯＤフローの算出例を示す説明図である。図１１において、対象期間を「６：００～１１：００」とする。また、第１の時間フレームの時間幅を「６０分」とし、第２の時間フレームの時間幅を「３０分」とする。
【0119】
この場合、第１の時間フレームは、例えば、「６：００～７：００」、「７：００～８：００」などである。また、第１の時間フレームを「６：００～７：００」とすると、第２の時間フレームは、「６：００～６：３０」と「６：３０～７：００」となる。また、第１の時間フレームを「７：００～８：００」とすると、第２の時間フレームは、「７：００～７：３０」と「７：３０～８：００」となる。
【0120】
ここで、あるＯＤペアの「６：００～７：００」の第１のＯＤフローを「ｑ１＿０」とし、「７：００～８：００」の第１のＯＤフローを「ｑ２＿０」とする。この場合、同じＯＤペアの「６：００～６：３０」の第２のＯＤフローを「ｑ１＿１」とすると、ｑ１＿１は、下記式（１）を用いて算出される。また、同じＯＤペアの「６：３０～７：００」の第２のＯＤフローを「ｑ２＿１」とすると、ｑ２＿１は、下記式（２）を用いて算出される。
【0121】
また、同じＯＤペアの「７：００～７：３０」の第２のＯＤフローを「ｑ３＿１」とすると、ｑ３＿１は、下記式（３）を用いて算出される。また、同じＯＤペアの「７：３０～８：００」の第２のＯＤフローを「ｑ４＿１」とすると、ｑ４＿１は、下記式（４）を用いて算出される。ただし、ｗ１は、０より大きく、１未満の定数である。ｗ２は、０より大きく、１未満の重み（定数）である。
【0122】
ｑ１＿１＝ｑ１＿０＊ｗ１ ・・・（１）
ｑ２＿１＝ｑ１＿０＊（１．０－ｗ１）・・・（２）
ｑ３＿１＝ｑ２＿０＊ｗ２ ・・・（３）
ｑ４＿１＝ｑ２＿０＊（１．０－ｗ２）・・・（４）
ｑ２＿１＝ｑ１＿０＊（１．０－ｗ１）・・・（２）
ｑ３＿１＝ｑ２＿０＊ｗ２ ・・・（３）
ｑ４＿１＝ｑ２＿０＊（１．０－ｗ２）・・・（４）
【0123】
重みｗ１，ｗ２は、例えば、第１の時間フレームと第２の時間フレームとの時間幅をもとに設定される。図１１の例では、第１の時間フレームと第２の時間フレームとの時間幅の比は、「２：１」である。この場合、重みｗ１は、「ｗ１＝１÷２＝０．５」となる。同様に、重みｗ２は、「ｗ２＝１÷２＝０．５」となる。
【0124】
また、重みｗ１，ｗ２は、例えば、出発地点（第１の地点）から目的地点（第２の地点）までの経路に含まれるリンクについて観測された第２の時間フレームごとの交通量の比をもとに設定されてもよい。
【0125】
ここで、出発地点を道路ネットワークＮＷ１内のｎｏｄｅ１とし、出発地点を道路ネットワークＮＷ１内のｎｏｄｅ４とする（図４参照）。ｎｏｄｅ１（Ｏ）からｎｏｄｅ４（Ｄ）までの経路は、「ｎｏｄｅ１－ｌｉｎｋ１－ｌｉｎｋ３－ｎｏｄｅ４」である。
【0126】
例えば、第１の時間フレーム「６：００～７：００」の第１のＯＤフローｑ１＿０が「ｑ１＿０＝１０両」であったとする。また、ｌｉｎｋ１に設置されたセンサ（車両カウンタ）により観測された、第２の時間フレーム「６：００～６：３０」の交通量（リンクフロー）が「２両」であり、第２の時間フレーム「６：３０～７：００」の交通量が「８両」であったとする。
【0127】
第２の時間フレーム「６：００～６：３０」の重みｗ１は、例えば、「ｗ１＝２／１０＝０．２」となる。この場合、第２の時間フレーム「６：００～６：３０」の第２のＯＤフローｑ１＿１は、上記式（１）から「１０×０．２＝２両」となる。また、第２の時間フレーム「６：３０～７：００」の第２のＯＤフローｑ２＿１は、上記式（２）から「１０×（１－０．２）＝８両」となる。
【0128】
（ＯＤフローの推定例）
つぎに、ＯＤ推定装置２０１におけるＯＤフローの推定例について説明する。ここでは、図５～図９に示した道路ネットワーク表５００、出発地点表６００、目的地点表７００、ＯＤペア表８００および経路表９００をもとに、図４に示した道路ネットワークＮＷ１における各ＯＤペアに対応する細かい時間フレーム（第２の時間フレーム）でのＯＤフローを推定する場合について説明する。
つぎに、ＯＤ推定装置２０１におけるＯＤフローの推定例について説明する。ここでは、図５～図９に示した道路ネットワーク表５００、出発地点表６００、目的地点表７００、ＯＤペア表８００および経路表９００をもとに、図４に示した道路ネットワークＮＷ１における各ＯＤペアに対応する細かい時間フレーム（第２の時間フレーム）でのＯＤフローを推定する場合について説明する。
【0129】
・粗い時間フレームでのＯＤフロー推定
ＯＤ推定装置２０１は、第１の交通シミュレーションにより、粗い時間フレーム（第１の時間フレーム）でのＯＤフローを推定する。ここでは、粗い時間フレーム（第１の時間フレーム）を６０分ごとの時間フレームとする。また、説明の都合上、対象期間を「ｈ１：６：００～７：００」とし、対象期間内に時間フレーム（第１の時間フレーム）が１つのみ含まれる場合を例に挙げて説明する。また、時間フレーム「ｈ１：６：００～７：００」を、単に時間フレームｈ１と表記する場合がある。
ＯＤ推定装置２０１は、第１の交通シミュレーションにより、粗い時間フレーム（第１の時間フレーム）でのＯＤフローを推定する。ここでは、粗い時間フレーム（第１の時間フレーム）を６０分ごとの時間フレームとする。また、説明の都合上、対象期間を「ｈ１：６：００～７：００」とし、対象期間内に時間フレーム（第１の時間フレーム）が１つのみ含まれる場合を例に挙げて説明する。また、時間フレーム「ｈ１：６：００～７：００」を、単に時間フレームｈ１と表記する場合がある。
【0130】
以下、粗い時間フレーム（第１の時間フレーム）でのＯＤフローを推定する際の具体的な処理内容について説明する。
【0131】
【0132】
まず、ＯＤ推定装置２０１は、事前準備として、初期ＯＤフロー１２０１を設定する。初期ＯＤフロー１２０１は、第１の交通シミュレーションに初期値として与えるＯＤフローであり、各ＯＤペアに対応する６０分ごとのＯＤフロー（例えば、乱数または定数）を示す。
【0133】
図１３は、初期ＯＤフローの具体例を示す説明図である。図１３において、初期ＯＤフロー１２０１は、各ＯＤペア１～３に対応する時間フレームｈ１でのＯＤフローを示す。各ＯＤフローは、車両数を示す。例えば、符号１３０１で示す「２」は、時間フレームｈ１にｎｏｄｅ２を出発してｎｏｄｅ４に到達する車両数（ＯＤフロー）を示す。
【0134】
つぎに、ＯＤ推定装置２０１は、Ｆｏｒｗａｒｄ処理１２１０（図１２参照）において、初期ＯＤフロー１２０１をＯＤフローｑ（第１のＯＤフロー）とする。そして、ＯＤ推定装置２０１は、ＯＤフローｑを用いて、第１の交通シミュレーションを実行する。
【0135】
具体的には、例えば、ＯＤ推定装置２０１は、道路ネットワーク表５００、ＯＤフローｑ（初期ＯＤフロー１２０１）、経路表９００、ＯＤペア表８００および経路選択表（図１４参照）を、交通シミュレータに与えることにより、第１の交通シミュレーションを実行する。
【0136】
ここで、図１４を用いて、経路選択表の具体例について説明する。
【0137】
図１４は、経路選択表の具体例を示す説明図（その１）である。図１４において、経路選択表１４００は、時間フレームｈ１：６：００～７：００の間に、道路ネットワークＮＷ１における各ＯＤペアに対応する経路（ｐａｔｈ）が選択される確率（経路選択確率）を示す。図１４中の＃は、経路（ｐａｔｈ）を識別する識別子である。
【0138】
例えば、ＯＤ＜ｎｏｄｅ１，ｎｏｄｅ４＞に対応するｐａｔｈ１（ｎｏｄｅ１－ｌｉｎｋ１－ｌｉｎｋ３－ｎｏｄｅ４）が選択される確率は、「１．０」である。各経路の経路選択確率は、時間帯に依存する。例えば、同じ経路であっても、朝６時台と昼１２時台での経路選択確率が異なる場合がある。
【0139】
図１４の例では、各ＯＤペアに対応する経路は、１本しかない。このため、各ＯＤペアについて、ＯからＤへ移動するときに他に選択できる経路がなく、各経路の経路選択確率は、「１．０」となる。ただし、同一のＯＤペアについて、複数の経路（ｐａｔｈ）が存在する場合がある。この場合、複数の経路の経路選択確率の合計は、「１．０」となる。
【0140】
第１の交通シミュレーションでは、例えば、「６：００～７：００」の間に５分ごとに仮想車両を出発地点から出発させ、目的地点に向けて走行させる。ただし、同じ時間フレーム内に出発した車両は、同じ時間フレーム内に目的地点に到達するものとする。例えば、「６：００～７：００」にｎｏｄｅ１（出発地点）を出発した車両は、「７：００」までにｎｏｄｅ４（目的地点）に全て到達する。渋滞は、ないものとする。
【0141】
第１の交通シミュレーションが完了すると、時間フレームｈ１における、各ＯＤペアに対応する経路（ｐａｔｈ）に含まれるリンクごとのリンクフロー（交通量）の推定結果が得られる。リンクフローは、例えば、リンクを通過した車両数である。
【0142】
【0143】
図１５は、リンクフロー（割合）の具体例を示す説明図（その１）である。図１５において、リンクフロー１５００は、第１の交通シミュレーションにより得られる情報（ρ）であり、時間フレームｈ１：６：００～７：００における、各ＯＤペアに対応する経路（ｐａｔｈ）に含まれるリンクごとのリンクフロー（割合）を表す。
【0144】
リンクフロー（割合）は、各ＯＤペアに対応する経路（ｐａｔｈ）を通って、出発地点から目的地点に向かう車両のうち、時間フレームｈ１に各リンクを通過した車両の割合を示す。リンクフロー（割合）は、例えば、図１４に示した経路選択表１４００内の各ＯＤペアに対応する経路（ｐａｔｈ）の経路選択確率に応じた値となる。
【0145】
例えば、符号１５０１で示す「１．０」は、時間フレームｈ１にｐａｔｈ１を通って、ｎｏｄｅ１からｎｏｄｅ４に向かう車両のうち、 ｌｉｎｋ１を通過する車両の割合を示す。また、符号１５０２で示す「０．０」は、時間フレームｈ１にｐａｔｈ３を通って、ｎｏｄｅ３からｎｏｄｅ４に向かう車両のうち、 ｌｉｎｋ２を通過する車両の割合を示す。
【0146】
ＯＤ推定装置２０１は、例えば、下記式（５）を用いて、リンクフロー１５００に基づいて、リンクフロー（交通量）を算出する。ただし、ｘは、リンクフロー（交通量）を示す。ρは、リンクフロー（割合）を示す。ｐは、経路選択確率を示す。ｑは、ＯＤフローを示す。
【0147】
ｘ＝ρｐｑ ・・・（５）
【0148】
【0149】
【数1】
【0150】
図１６は、リンクフロー（交通量）の具体例を示す説明図（その１）である。図１６において、リンクフロー１６００は、第１の交通シミュレーションにより推定される情報（ｘ）であり、時間フレームｈ１：６：００～７：００における、各ＯＤペアに対応する経路（ｐａｔｈ）に含まれるリンクごとのリンクフロー（交通量）を表す。
【0151】
リンクフロー（交通量）は、時間フレームｈ１に各リンクを通過した車両数を示す。例えば、符号１６０１で示す「２」は、時間フレームｈ１にｌｉｎｋ２を通った車両数を示す。これにより、ＯＤ推定装置２０１は、各ＯＤペアに対応する経路（ｐａｔｈ）に含まれる全てのリンクでの交通量（車両数）を推定することができる。
【0152】
ここで、シミュレータ上では全てのリンクの交通量を算出可能なものの、現実世界では、道路ネットワーク上の全てのリンク（道路）にセンサ（車両カウンタなど）を設置することは難しい。このため、ＯＤ推定装置２０１は、実際に交通量を観測可能なリンクの交通量を抽出する。
【0153】
ここでは、道路ネットワークＮＷ１内のｌｉｎｋ１，ｌｉｎｋ３のみにセンサ（車両カウンタ）が設置されていると仮定する。このため、ＯＤ推定装置２０１は、図１７に示すような観測行列１７００を用いて、道路ネットワークＮＷ１から、交通量を観測可能なリンクを抽出する。
【0154】
図１７は、観測行列の一例を示す説明図（その１）である。図１７において、観測行列１７００は、道路ネットワークＮＷ１内のｌｉｎｋ１，ｌｉｎｋ２，ｌｉｎｋ３のうち、時間フレームｈ１に交通量を観測可能なｌｉｎｋ１，ｌｉｎｋ３を抽出するための行列である。例えば、符号１７０１で示す「１」は、時間フレームｈ１にｌｉｎｋ１を通る交通量（車両数）を観測可能であることを示すフラグである。
【0155】
なお、図１７中、ｂ１１は、時間フレームｈ１のｌｉｎｋ１を抽出するための情報に対応する。ｂ１２は、時間フレームｈ１のｌｉｎｋ３を抽出するための情報に対応する。
【0156】
ＯＤ推定装置２０１は、例えば、下記式（７）を用いて、道路ネットワークＮＷ１内のｌｉｎｋ１，ｌｉｎｋ２，ｌｉｎｋ３のうち、時間フレームｈ１に交通量を観測可能なｌｉｎｋ１，ｌｉｎｋ３のリンクフロー（交通量）を抽出する。ただし、ｙは、観測可能なリンクフロー（交通量）を示す。Ｌは、観測行列を示す。ｘは、シミュレータ上で算出されたリンクフロー（交通量）を示す。
【0157】
ｙ＝Ｌｘ ・・・（７）
【0158】
【0159】
【数2】
【0160】
図１８は、観測可能なリンクフロー（交通量）の具体例を示す説明図（その１）である。図１８において、リンクフロー１８００は、時間フレームｈ１：６：００～７：００に実際に観測されたｌｉｎｋ１，ｌｉｎｋ３のリンクフロー（交通量）を表す。例えば、符号１８０１で示す「６」は、時間フレームｈ１にｌｉｎｋ２を通った車両数を示す。図１８中の＃は、観測されたリンクフロー（交通量）を識別する識別子である。
【0161】
つぎに、ＯＤ推定装置２０１は、例えば、下記式（９）を用いて、シミュレータ上で算出（推定）されたリンクフロー（交通量）と、実際に観測されたリンクフロー（交通量）との誤差を算出する。ただし、Δｙは、リンクフロー（交通量）の誤差を示す。ｙ’は、実際に観測されたリンクフロー（交通量）を示す。ｙは、シミュレータ上で算出されたリンクフロー（交通量）を示す。
【0162】
Δｙ＝ｙ’－ｙ ・・・（９）
【0163】
【0164】
図１９は、リンクフロー（交通量）の誤差の具体例を示す説明図（その１）である。図１９において、リンクフロー（交通量）の誤差１９００は、時間フレームｈ１：６：００～７：００における、交通量を観測可能なｌｉｎｋ１，ｌｉｎｋ３のリンクフロー（交通量）の誤差を表す。なお、リンクフロー１９０１は、実際に観測されたリンクフロー（交通量）であり、リンクフロー観測データから特定される。
【0165】
ここで、算出されたリンクフロー（交通量）の誤差１９００が、あらかじめ決められた許容範囲内ではなかった場合を想定する。
【0166】
この場合、ＯＤ推定装置２０１は、Ｂａｃｋｗａｒｄ処理１２２０（図１２参照）において、リンクフロー（交通量）の誤差１９００に基づいて、ＯＤフローｑを更新する。この時点の更新前のＯＤフローｑは、図１３に示した初期ＯＤフロー１２０１である。
【0167】
具体的には、例えば、まず、ＯＤ推定装置２０１は、下記式（１０）を用いて、ＯＤフローｑの補正値を算出する。ただし、Δｑは、ＯＤフローｑの補正値を示す。ｗは、重み（定数）を示す。ｗは、例えば、０～１の値である。ｐTは、ｐの転置行列を示す。ρTは、ρの転置行列を示す。ＬTは、Ｌの転置行列を示す。Δｙは、リンクフロー（交通量）の誤差を示す。
【0168】
Δｑ＝ｗｐTρTＬTΔｙ ・・・（１０）
【0169】
ここで、ｗを「ｗ＝０．５」とし、ｐを図１４に示した経路選択表１４００とし、ρを、図１５に示したリンクフロー１５００とし、Ｌを、図１７に示した観測行列１７００とし、Δｙを、図１９に示したリンクフロー（交通量）の誤差１９００とする。この場合、Δｑは、下記式（１１）のように計算される。
【0170】
【数3】
【0171】
つぎに、ＯＤ推定装置２０１は、例えば、下記式（１２）を用いて、ＯＤフローｑを更新する。
【0172】
ｑ＝Δｑ＋ｑ ・・・（１２）
【0173】
ここで、Δｑを、上記式（１１）のように計算された値とすると、ＯＤフローｑは、図２０に示すように更新される。
【0174】
図２０は、ＯＤフローの更新例を示す説明図（その１）である。図２０において、ＯＤフロー２０００は、図１３に示した初期ＯＤフロー１２０１を、上記式（１０）のように計算されたΔｑ（ＯＤフローｑの補正値２００１）を用いて更新した更新後のＯＤフローｑを示す。
【0175】
つぎに、ＯＤ推定装置２０１は、ＯＤフロー２０００をＯＤフローｑ（第１のＯＤフロー）とする。そして、ＯＤ推定装置２０１は、Ｆｏｒｗａｒｄ処理１２１０（図１２参照）に戻って、ＯＤフローｑを用いて、第１の交通シミュレーションを実行する。
【0176】
ＯＤ推定装置２０１は、Ｆｏｒｗａｒｄ処理１２１０およびＢａｃｋｗａｒｄ処理１２２０を、算出されるリンクフロー（交通量）の誤差が許容範囲内となるまで反復する。ここでは、反復を４回行った結果、ＯＤフロー１２０２が推定された場合を想定する。
【0177】
図２１は、粗い時間フレームでのＯＤフローの推定結果を示す説明図である。図２１において、ＯＤフロー１２０２は、粗い時間フレーム（第１の時間フレーム）でのＯＤフローの推定結果を示す。粗い時間フレーム（第１の時間フレーム）は、６０分ごとの時間フレームである。
【0178】
・細かい時間フレームでのＯＤフロー推定
つぎに、ＯＤ推定装置２０１は、第２の交通シミュレーションにより、細かい時間フレーム（第２の時間フレーム）でのＯＤフローを推定する。ここでは、細かい時間フレーム（第２の時間フレーム）を３０分ごとの時間フレームとする。また、粗い時間フレーム「６：００～７：００」を３０分間隔で分割した細かい時間フレームを、「ｈ１：６：００～６：３０」および「ｈ２：６：３０～７：００」と表記する場合がある。また、時間フレーム「ｈ１：６：００～６：３０」を、単に時間フレームｈ１と表記し、時間フレーム「ｈ２：６：３０～７：００」を、単に時間フレームｈ２と表記する場合がある。
つぎに、ＯＤ推定装置２０１は、第２の交通シミュレーションにより、細かい時間フレーム（第２の時間フレーム）でのＯＤフローを推定する。ここでは、細かい時間フレーム（第２の時間フレーム）を３０分ごとの時間フレームとする。また、粗い時間フレーム「６：００～７：００」を３０分間隔で分割した細かい時間フレームを、「ｈ１：６：００～６：３０」および「ｈ２：６：３０～７：００」と表記する場合がある。また、時間フレーム「ｈ１：６：００～６：３０」を、単に時間フレームｈ１と表記し、時間フレーム「ｈ２：６：３０～７：００」を、単に時間フレームｈ２と表記する場合がある。
【0179】
以下、細かい時間フレーム（第２の時間フレーム）でのＯＤフローを推定する際の具体的な処理内容について説明する。
【0180】
図２２は、ＯＤ推定装置２０１の第２の動作例を示す説明図である。図２２において、ＯＤ推定装置２０１が、細かい時間フレーム（第２の時間フレーム）でのＯＤフローを推定する際に反復して行う動作が示されている。
【0181】
まず、ＯＤ推定装置２０１は、図２１に示したＯＤフロー１２０２（粗い時間フレームでのＯＤフロー）に基づいて、第２の交通シミュレーションに初期値として与えるＯＤフロー２２０１を設定する。
【0182】
具体的には、例えば、ＯＤ推定装置２０１は、図２１に示したＯＤフロー１２０２（粗い時間フレームでのＯＤフロー）を、細かい時間フレーム（第２の時間フレーム）ごとに配分することにより、ＯＤフロー２２０１を算出する。ここでは、重みを０．５として、粗い時間フレームでのＯＤフローを「１：１」で、細かい時間フレーム（第２の時間フレーム）に配分する場合を想定する。
【0183】
図２３は、ＯＤフローの配分例を示す説明図である。図２３において、ＯＤフロー２２０１は、粗い時間フレームでのＯＤフロー１２０２（図２１参照）を、細かい時間フレームごとに配分することにより算出されたＯＤフローを示す。
【0184】
例えば、符号２３０１で示す「３１５．５」は、時間フレームｈ１にｎｏｄｅ２を出発してｎｏｄｅ４に到達する車両数（ＯＤフロー）を示す。また、符号２３０２で示す「３１５．５」は、時間フレームｈ２にｎｏｄｅ３を出発してｎｏｄｅ４に到達する車両数（ＯＤフロー）を示す。
【0185】
つぎに、ＯＤ推定装置２０１は、Ｆｏｒｗａｒｄ処理２２１０（図２２参照）において、ＯＤフロー２２０１をＯＤフローｑ（第２のＯＤフロー）とする。そして、ＯＤ推定装置２０１は、ＯＤフローｑを用いて、第２の交通シミュレーションを実行する。
【0186】
具体的には、例えば、ＯＤ推定装置２０１は、道路ネットワーク表５００、ＯＤフローｑ（ＯＤフロー２２０１）、経路表９００、ＯＤペア表８００および経路選択表（図２４参照）を、交通シミュレータに与えることにより、第２の交通シミュレーションを実行する。
【0187】
図２４は、経路選択表の具体例を示す説明図（その２）である。図２４において、経路選択表２４００は、各時間フレームｈ１：６：００～６：３０，ｈ２：６：３０～７：００の間に、道路ネットワークＮＷ１における各ＯＤペアに対応する経路（ｐａｔｈ）が選択される確率（経路選択確率）を示す。図２４中の＃は、経路（ｐａｔｈ）を識別する識別子である。
【0188】
例えば、時間フレームｈ１：６：００～６：３０に、ＯＤ＜ｎｏｄｅ１，ｎｏｄｅ４＞に対応するｐａｔｈ１（ｎｏｄｅ１－ｌｉｎｋ１－ｌｉｎｋ３－ｎｏｄｅ４）が選択される確率は、「１．０」である。また、時間フレームｈ２：６：３０～７：００に、ＯＤ＜ｎｏｄｅ３，ｎｏｄｅ４＞に対応するｐａｔｈ３（ｎｏｄｅ３－ｌｉｎｋ３－ｎｏｄｅ４）が選択される確率は、「１．０」である。
【0189】
第２の交通シミュレーションでは、例えば、「６：００～７：００」の間に５分ごとに仮想車両を出発地点から出発させ、目的地点に向けて走行させる。ただし、同じ時間フレーム内に出発した車両は、同じ時間フレーム内に目的地点に到達するものとする。例えば、「６：００～７：００」にｎｏｄｅ１（出発地点）を出発した車両は、「７：００」までにｎｏｄｅ４（目的地点）に全て到達する。渋滞は、ないものとする。
【0190】
第２の交通シミュレーションが完了すると、各時間フレームｈ１，ｈ２における、各ＯＤペアに対応する経路（ｐａｔｈ）に含まれるリンクごとのリンクフロー（交通量）の推定結果が得られる。
【0191】
【0192】
図２５は、リンクフロー（割合）の具体例を示す説明図（その２）である。図２５において、リンクフロー２５００は、第２の交通シミュレーションにより得られる情報（ρ）であり、各時間フレームｈ１：６：００～６：３０，ｈ２：６：３０～７：００における、各ＯＤペアに対応する経路（ｐａｔｈ）に含まれるリンクごとのリンクフロー（割合）を表す。
【0193】
リンクフロー（割合）は、例えば、図２４に示した経路選択表２４００内の各ＯＤペアに対応する経路（ｐａｔｈ）の経路選択確率に応じた値となる。
【0194】
例えば、符号２５０１で示す「１．０」は、時間フレームｈ１にｐａｔｈ２を通って、ｎｏｄｅ２からｎｏｄｅ４に向かう車両のうち、 時間フレームｈ１にｌｉｎｋ３を通過する車両の割合を示す。また、符号２５０２で示す「０．０」は、時間フレームｈ１にｐａｔｈ１を通って、ｎｏｄｅ１からｎｏｄｅ４に向かう車両のうち、時間フレームｈ２にｌｉｎｋ２を通過する車両の割合を示す。
【0195】
ＯＤ推定装置２０１は、例えば、上記式（５）を用いて、リンクフロー２５００に基づいて、リンクフロー（交通量）を算出する。ここで、ρを、図２５に示したリンクフロー１５００とすると、ｘは、下記式（１３）のように計算され、図２６に示すようなリンクフロー２６００となる。
【0196】
【数4】
【0197】
図２６は、リンクフロー（交通量）の具体例を示す説明図（その２）である。図２６において、リンクフロー２６００は、第２の交通シミュレーションにより推定される情報（ｘ）であり、各時間フレームｈ１：６：００～６：３０，ｈ２：６：３０～７：００における、各ＯＤペアに対応する経路（ｐａｔｈ）に含まれるリンクごとのリンクフロー（交通量）を表す。
【0198】
例えば、符号２６０１で示す「３１５．５」は、時間フレームｈ１にｌｉｎｋ２を通過した車両数を示す。また、符号２６０２で示す「９１６．５」は、時間フレームｈ２にｌｉｎｋ３を通過した車両数を示す。これにより、ＯＤ推定装置２０１は、各ＯＤペアに対応する経路（ｐａｔｈ）に含まれる全てのリンクでの交通量（車両数）を推定することができる。
【0199】
つぎに、ＯＤ推定装置２０１は、実際に交通量を観測可能なリンクの交通量を抽出する。ここでは、道路ネットワークＮＷ１内のｌｉｎｋ１，ｌｉｎｋ３のみにセンサ（車両カウンタ）が設置されていると仮定する。このため、ＯＤ推定装置２０１は、図２７に示すような観測行列２７００を用いて、道路ネットワークＮＷ１から、交通量を観測可能なリンクを抽出する。
【0200】
図２７は、観測行列の一例を示す説明図（その２）である。図２７において、観測行列２７００は、道路ネットワークＮＷ１内のｌｉｎｋ１，ｌｉｎｋ２，ｌｉｎｋ３のうち、各時間フレームｈ１，ｈ２に交通量を観測可能なｌｉｎｋ１，ｌｉｎｋ３を抽出するための行列である。
【0201】
例えば、符号２７０１で示す「１」は、時間フレームｈ１にｌｉｎｋ１を通る交通量（車両数）を観測可能であることを示すフラグである。また、符号２７０２で示す「１」は、時間フレームｈ１にｌｉｎｋ３を通る交通量を観測可能であることを示すフラグである。また、符号２７０３で示す「１」は、時間フレームｈ２にｌｉｎｋ１を通る交通量を観測可能であることを示すフラグである。また、符号２７０４で示す「１」は、時間フレームｈ２にｌｉｎｋ３を通る交通量を観測可能であることを示すフラグである。
【0202】
なお、図２７中、ｂ１１は、時間フレームｈ１のｌｉｎｋ１を抽出するための情報に対応する。ｂ１２は、時間フレームｈ１のｌｉｎｋ３を抽出するための情報に対応する。ｂ２１は、時間フレームｈ２のｌｉｎｋ１を抽出するための情報に対応する。ｂ２２は、時間フレームｈ２のｌｉｎｋ３を抽出するための情報に対応する。
【0203】
ＯＤ推定装置２０１は、例えば、上記式（７）を用いて、道路ネットワークＮＷ１内のｌｉｎｋ１，ｌｉｎｋ２，ｌｉｎｋ３のうち、各時間フレームｈ１，ｈ２に交通量を観測可能なｌｉｎｋ１，ｌｉｎｋ３のリンクフロー（交通量）を抽出する。
【0204】
【0205】
【数5】
【0206】
図２８は、観測可能なリンクフロー（交通量）の具体例を示す説明図（その２）である。図２８において、リンクフロー２８００は、各時間フレームｈ１：６：００～６：３０，ｈ２：６：３０～７：００に実際に観測されたｌｉｎｋ１，ｌｉｎｋ３のリンクフロー（交通量）を表す。
【0207】
例えば、符号２８０１で示す「９１６．５」は、時間フレームｈ１にｌｉｎｋ３を通った車両数を示す。図２８中の＃は、観測されたリンクフロー（交通量）を識別する識別子である。
【0208】
つぎに、ＯＤ推定装置２０１は、例えば、上記式（９）を用いて、シミュレータ上で算出（推定）されたリンクフロー（交通量）と、実際に観測されたリンクフロー（交通量）との誤差を算出する。
【0209】
【0210】
図２９は、リンクフロー（交通量）の誤差の具体例を示す説明図（その２）である。図２９において、リンクフロー（交通量）の誤差２９００は、各時間フレームｈ１：６：００～６：３０，ｈ２：６：３０～７：００における、交通量を観測可能なｌｉｎｋ１，ｌｉｎｋ３のリンクフロー（交通量）の誤差を表す。なお、リンクフロー２９０１は、実際に観測されたリンクフロー（交通量）であり、リンクフロー観測データから特定される。
【0211】
ここで、算出されたリンクフロー（交通量）の誤差２９００が、あらかじめ決められた許容範囲内ではなかった場合を想定する。
【0212】
この場合、ＯＤ推定装置２０１は、Ｂａｃｋｗａｒｄ処理２２２０（図２２参照）において、リンクフロー（交通量）の誤差２９００に基づいて、ＯＤフローｑを更新する。この時点の更新前のＯＤフローｑは、図２３に示したＯＤフロー２２０１である。
【0213】
具体的には、例えば、まず、ＯＤ推定装置２０１は、上記式（１０）を用いて、ＯＤフローｑの補正値を算出する。ここで、ｗを「ｗ＝０．５」とし、ｐを図２４に示した経路選択表２４００とし、ρを、図２５に示したリンクフロー２５００とし、Ｌを、図２７にしめした観測行列２７００とし、Δｙを、図２９に示したリンクフロー（交通量）の誤差２９００とする。この場合、Δｑは、下記式（１５）のように計算される。
【0214】
【数6】
【0215】
つぎに、ＯＤ推定装置２０１は、例えば、上記式（１２）を用いて、ＯＤフローｑを更新する。
【0216】
ここで、Δｑを、上記式（１５）のように計算された値とすると、ＯＤフローｑは、図３０に示すように更新される。
【0217】
図３０は、ＯＤフローの更新例を示す説明図（その２）である。図３０において、ＯＤフロー３０００は、図２３に示したＯＤフロー２２０１を、上記式（１４）のように計算されたΔｑ（ＯＤフローｑの補正値３００１）を用いて更新した更新後のＯＤフローｑを示す。
【0218】
つぎに、ＯＤ推定装置２０１は、ＯＤフロー３０００をＯＤフローｑ（第２のＯＤフロー）とする。そして、ＯＤ推定装置２０１は、Ｆｏｒｗａｒｄ処理２２１０（図２２参照）に戻って、ＯＤフローｑを用いて、第２の交通シミュレーションを実行する。
【0219】
ＯＤ推定装置２０１は、Ｆｏｒｗａｒｄ処理２２１０およびＢａｃｋｗａｒｄ処理２２２０を、算出されるリンクフロー（交通量）の誤差が許容範囲内となるまで反復する。ここでは、反復を１０回行った結果、ＯＤフロー２２０２が推定された場合を想定する。
【0220】
図３１は、細かい時間フレームでのＯＤフローの推定結果を示す説明図である。図３１において、ＯＤフロー２２０２は、細かい時間フレーム（第２の時間フレーム）でのＯＤフローの推定結果を示す。細かい時間フレーム（第２の時間フレーム）は、３０分ごとの時間フレームである。
【0221】
これにより、ＯＤ推定装置２０１は、ＯＤフローの推定にかかる処理時間を抑えて、細かい時間フレーム（第２の時間フレーム）でのＯＤフローを推定することができる。
【0222】
【0223】
図３２および図３３は、ＯＤ推定装置２０１のＯＤ決定処理手順の一例を示すフローチャートである。図３２のフローチャートにおいて、まず、ＯＤ推定装置２０１は、第１の交通シミュレーションを実行するにあたり、初期ＯＤフローを設定する（ステップＳ３２０１）。
【0224】
つぎに、ＯＤ推定装置２０１は、第１の交通シミュレーションにおける第１の時間フレームの幅（時間幅）を設定する（ステップＳ３２０２）。第１の時間フレームの幅には、初期ＯＤフローの時間フレームの幅が設定される。
【0225】
そして、ＯＤ推定装置２０１は、初期ＯＤフローを第１のＯＤフローに設定し（ステップＳ３２０３）、第１のＯＤフローを用いて、第１の交通シミュレーションを実行する（ステップＳ３２０４）。つぎに、ＯＤ推定装置２０１は、第１の交通シミュレーションにより得られる第１の時間フレームごとのリンクフロー（割合）から、第１の時間フレームごとのリンクフロー（交通量）を算出する（ステップＳ３２０５）。
【0226】
そして、ＯＤ推定装置２０１は、算出したリンクフロー（交通量）から、交通量を観測可能なリンクのリンクフロー（交通量）を抽出する（ステップＳ３２０６）。つぎに、ＯＤ推定装置２０１は、抽出したリンクのリンクフロー（交通量）と当該リンクについて実際に観測されたリンクフロー（交通量）との誤差を算出する（ステップＳ３２０７）。実際に観測されたリンクフロー（交通量）は、リンクフロー観測データから特定される。
【0227】
そして、ＯＤ推定装置２０１は、算出したリンクフロー（交通量）の誤差が許容範囲内であるか否かを判断する（ステップＳ３２０８）。ここで、誤差が許容範囲内ではない場合（ステップＳ３２０８：Ｎｏ）、ＯＤ推定装置２０１は、算出した誤差に基づいて、第１のＯＤフローの補正値を算出する（ステップＳ３２０９）。
【0228】
つぎに、ＯＤ推定装置２０１は、算出した補正値を用いて、第１のＯＤフローを更新する（ステップＳ３２１０）。そして、ＯＤ推定装置２０１は、ステップＳ３２０３に戻って、更新後の第１のＯＤフローを、第１のＯＤフローに設定して、以降の処理を実行する。
【0229】
また、ステップＳ３２０８において、算出したリンクフロー（交通量）の誤差が許容範囲内の場合（ステップＳ３２０８：Ｙｅｓ）、ＯＤ推定装置２０１は、直前の第１の交通シミュレーションに用いた第１のＯＤフローを、第１の時間フレームごとのＯＤフローに決定する（ステップＳ３２１１）。
【0230】
そして、ＯＤ推定装置２０１は、決定した第１の時間フレームごとのＯＤフローを、第１の時間フレームを分割した第２の時間フレームごとに配分することにより、第２のＯＤフローを算出して（ステップＳ３２１２）、図３３に示すステップＳ３３０１に移行する。
【0231】
図３３のフローチャートにおいて、まず、ＯＤ推定装置２０１は、第２の交通シミュレーションを実行するにあたり、初期ＯＤフローを設定する（ステップＳ３３０１）。初期ＯＤフローには、図３２に示したステップＳ３２１２において算出された第２のＯＤフローが設定される。
【0232】
つぎに、ＯＤ推定装置２０１は、第２の交通シミュレーションにおける第２の時間フレームの幅を設定する（ステップＳ３３０２）。そして、ＯＤ推定装置２０１は、初期ＯＤフローを第２のＯＤフローに設定し（ステップＳ３３０３）、第２のＯＤフローを用いて、第２の交通シミュレーションを実行する（ステップＳ３３０４）。
【0233】
つぎに、ＯＤ推定装置２０１は、第２の交通シミュレーションにより得られる第２の時間フレームごとのリンクフロー（割合）から、第２の時間フレームごとのリンクフロー（交通量）を算出する（ステップＳ３３０５）。そして、ＯＤ推定装置２０１は、算出したリンクフロー（交通量）から、交通量を観測可能なリンクのリンクフロー（交通量）を抽出する（ステップＳ３３０６）。
【0234】
つぎに、ＯＤ推定装置２０１は、抽出したリンクのリンクフロー（交通量）と当該リンクについて実際に観測されたリンクフロー（交通量）との誤差を算出する（ステップＳ３３０７）。そして、ＯＤ推定装置２０１は、算出したリンクフロー（交通量）の誤差が許容範囲内であるか否かを判断する（ステップＳ３３０８）。
【0235】
ここで、誤差が許容範囲内ではない場合（ステップＳ３３０８：Ｎｏ）、ＯＤ推定装置２０１は、算出した誤差に基づいて、第２のＯＤフローの補正値を算出する（ステップＳ３３０９）。つぎに、ＯＤ推定装置２０１は、算出した補正値を用いて、第２のＯＤフローを更新する（ステップＳ３３１０）。
【0236】
そして、ＯＤ推定装置２０１は、ステップＳ３３０３に戻って、更新後の第２のＯＤフローを、第２のＯＤフローに設定して、以降の処理を実行する。
【0237】
また、ステップＳ３３０８において、算出したリンクフロー（交通量）の誤差が許容範囲内の場合（ステップＳ３３０８：Ｙｅｓ）、ＯＤ推定装置２０１は、直前の第２の交通シミュレーションに用いた第２のＯＤフローを、第２の時間フレームごとのＯＤフローに決定する（ステップＳ３３１１）。
【0238】
そして、ＯＤ推定装置２０１は、決定した第２の時間フレームごとのＯＤフローを出力して（ステップＳ３３１２）、本フローチャートによる一連の処理を終了する。
【0239】
これにより、ＯＤ推定装置２０１は、粗い時間フレーム（第１の時間フレーム）で推定したＯＤフローを利用して、細かい時間フレーム（第２の時間フレーム）でのＯＤフローを推定することができる。
【0240】
以上説明したように、実施の形態にかかるＯＤ推定装置２０１によれば、対象期間を分割した第１の時間フレームごとに、第１のＯＤフローを用いて、第１の地点から第２の地点までの経路に含まれるリンクごとの交通量を推定する第１の交通シミュレーションを実行することができる。第１のＯＤフローは、第１の時間フレームに第１の地点を出発して第２の地点に到達する交通量（例えば、車両数）をあらわす。そして、ＯＤ推定装置２０１によれば、推定したリンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、第１のＯＤフローを調整しながら第１の交通シミュレーションを繰り返し実行することにより、第１の地点と第２の地点とのペア（ＯＤペア）について、第１の時間フレームごとのＯＤフローを決定することができる。
【0241】
これにより、ＯＤ推定装置２０１は、出発地点（Ｏ）と目的地点（Ｄ）とのＯＤペアについて、本来推定したい時間フレームよりも粗い時間フレーム（第１の時間フレーム）でのＯＤフローを推定することができる。粗い時間フレームでのＯＤフロー推定では、最適解を得るまでの反復回数が多くなっても、取り扱い行列サイズ（Ｌ，ρ，ｐ，ｑ）が比較的小さいため、１反復にかかる処理時間は短くなることが期待できる。
【0242】
また、ＯＤ推定装置２０１によれば、決定した第１の時間フレームごとのＯＤフローに基づいて、第１の時間フレームを分割した第２の時間フレームごとに、当該第２の時間フレームに第１の地点を出発して第２の地点に到達する第２のＯＤフローを算出することができる。算出された第２のＯＤフローは、第２の交通シミュレーションの初期ＯＤフローに相当する。また、ＯＤ推定装置２０１によれば、算出した第２のＯＤフローを用いて、第２の時間フレームごとに、リンクごとの交通量を推定する第２の交通シミュレーションを実行することができる。そして、ＯＤ推定装置２０１によれば、推定したリンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、第２のＯＤフローを調整しながら第２の交通シミュレーションを繰り返し実行することにより、第１の地点と第２の地点とのペアについて、第２の時間フレームごとのＯＤフローを決定することができる。
【0243】
これにより、ＯＤ推定装置２０１は、出発地点（Ｏ）と目的地点（Ｄ）とのＯＤペアについて、粗い時間フレーム（第１の時間フレーム）でのＯＤフローを利用して、本来推定したい時間フレーム（第２の時間フレーム）でのＯＤフローを推定することができる。ＯＤ推定装置２０１は、細かい時間フレームでのＯＤフロー推定の初期値として、粗い時間フレームで推定されたＯＤフローを利用することで、最適解を得るまでの反復回数を減らすことができる。
【0244】
また、ＯＤ推定装置２０１によれば、決定した第１の時間フレームごとのＯＤフローを、所定の比率に応じて第２の時間フレームごとに配分することにより、第２のＯＤフローを算出することができる。
【0245】
これにより、ＯＤ推定装置２０１は、第１の時間フレームごとのＯＤフローを、当該第１の時間フレームから分割された第２の時間フレームそれぞれに分配して利用することができる。
【0246】
また、ＯＤ推定装置２０１によれば、第１の地点と第２の地点とのペア（ＯＤペア）と対応付けて、決定した第２の時間フレームごとのＯＤフローを出力することができる。
【0247】
これにより、ＯＤ推定装置２０１は、出発地点（Ｏ）と目的地点（Ｄ）とのＯＤペアについて、細かい時間フレーム（第２の時間フレーム）でのＯＤフローをＯＤ推定結果として出力することができる。
【0248】
また、ＯＤ推定装置２０１によれば、第１の時間フレームと第２の時間フレームとの時間幅をもとに設定される比率に応じて、決定した第１の時間フレームごとのＯＤフローを第２の時間フレームごとに配分することにより、第２のＯＤフローを算出することができる。
【0249】
これにより、ＯＤ推定装置２０１は、第１の時間フレームと第２の時間フレームとの時間幅の比に応じて、第１の時間フレームごとのＯＤフローを、第２の時間フレームそれぞれに分配して利用することができる。例えば、第１の時間フレームと第２の時間フレームとの時間幅の比が「２：１」の場合、第１の時間フレームのＯＤフローを２分割したものが第２の時間フレームそれぞれに分配される。このため、ＯＤ推定装置２０１は、第２の交通シミュレーションの初期値として利用するＯＤフローの算出精度を高めることができる。
【0250】
また、ＯＤ推定装置２０１によれば、第１の地点から第２の地点までの経路に含まれる第１のリンクについて観測された、第２の時間フレームごとの交通量をもとに設定される比率に応じて、決定した第１の時間フレームごとのＯＤフローを第２の時間フレームごとに配分することにより、第２のＯＤフローを算出することができる。第１のリンクは、例えば、第１の地点に接続されたリンクである。
【0251】
これにより、ＯＤ推定装置２０１は、第２の時間フレームごとに実際に観測されたリンクフロー（交通量）の比に応じて、第１の時間フレームごとのＯＤフローを、第２の時間フレームそれぞれに分配して利用することができる。このため、ＯＤ推定装置２０１は、第２の交通シミュレーションの初期値として利用するＯＤフローの算出精度を高めることができる。
【0252】
また、ＯＤ推定装置２０１によれば、第１の交通シミュレーションにより推定したリンクの交通量と、当該リンクについて観測された交通量との誤差を算出し、算出した誤差が許容範囲内であるか否かを判断し、誤差が許容範囲内ではない場合、誤差に基づいて、第１のＯＤフローを更新し、更新した第１のＯＤフローを用いて、第１の交通シミュレーションを実行することができる。そして、ＯＤ推定装置２０１によれば、算出した誤差が許容範囲内の場合、直前の第１の交通シミュレーションに用いた第１のＯＤフローを、第１の時間フレームごとのＯＤフローに決定することができる。
【0253】
これにより、ＯＤ推定装置２０１は、第１の交通シミュレーションにより推定したリンクフローと、実際に観測されたリンクフローとの誤差が許容範囲内となるまで、第１のＯＤフローを更新しながら、第１の交通シミュレーションを繰り返すことができる。
【0254】
また、ＯＤ推定装置２０１によれば、第２の交通シミュレーションにより推定したリンクの交通量と、当該リンクについて観測された交通量との誤差を算出し、算出した誤差が許容範囲内であるか否かを判断し、誤差が許容範囲内ではない場合、誤差に基づいて、第２のＯＤフローを更新し、更新した第２のＯＤフローを用いて、第２の交通シミュレーションを実行することができる。そして、ＯＤ推定装置２０１によれば、算出した誤差が許容範囲内の場合、直前の第２の交通シミュレーションに用いた第２のＯＤフローを、第２の時間フレームごとのＯＤフローに決定することができる。
【0255】
これにより、ＯＤ推定装置２０１は、第２の交通シミュレーションにより推定したリンクフローと、実際に観測されたリンクフローとの誤差が許容範囲内となるまで、第２のＯＤフローを更新しながら、第２の交通シミュレーションを繰り返すことができる。
【0256】
これらのことから、ＯＤ推定装置２０１によれば、粗い時間フレームで推定したＯＤフローを利用して、細かい時間フレームでのＯＤフローを推定することができる。これにより、所望の時間フレームでのＯＤフローの推定にかかる処理時間を大幅に短縮することができる。
【0257】
例えば、ＯＤ推定装置２０１は、時間フレームの幅が大きい層では、時間フレームの幅が小さい層に比べて、取り扱い行列サイズを小さくできるため、全体の処理時間を短縮することができる。また、ＯＤ推定装置２０１は、時間フレームの幅の小さい層では、ＯＤフローが既に真値に近い値に近づいているため、最適化されるまでにかかる反復回数を少なくすることができる。
【0258】
また、本ＯＤ決定方法によれば、細かい時間フレームでのＯＤフロー（交通量）を短時間かつ高精度に推定することが可能となる。このため、例えば、ＯＤフローを推定するサービスをクラウドで提供することが可能となる。また、交通局やタクシーなどの配車サービスを提供する企業において、交通施策やサービス向上にＯＤ推定結果を役立てることができる。
【0259】
なお、本実施の形態で説明したＯＤ決定方法は、あらかじめ用意されたプログラムをパーソナル・コンピュータやワークステーション等のコンピュータで実行することにより実現することができる。本ＯＤ決定プログラムは、ハードディスク、フレキシブルディスク、ＣＤ－ＲＯＭ、ＤＶＤ、ＵＳＢメモリ等のコンピュータで読み取り可能な記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行される。また、本ＯＤ決定プログラムは、インターネット等のネットワークを介して配布してもよい。
【0260】
また、本実施の形態で説明したＯＤ推定装置２０１は、スタンダードセルやストラクチャードＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）などの特定用途向けＩＣやＦＰＧＡなどのＰＬＤ（Ｐｒｏｇｒａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅｖｉｃｅ）によっても実現することができる。
【0261】
上述した実施の形態に関し、さらに以下の付記を開示する。
【0262】
（付記１）対象期間を分割した第１の時間フレームごとに、当該第１の時間フレームに第１の地点を出発して第２の地点に到達する第１のＯＤ交通量を用いて、前記第１の地点から前記第２の地点までの経路に含まれるリンクごとの交通量を推定する第１の交通シミュレーションを実行し、
推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第１のＯＤ交通量を調整しながら前記第１の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第１の時間フレームごとのＯＤ交通量を決定し、
決定した前記第１の時間フレームごとのＯＤ交通量に基づいて、前記第１の時間フレームを分割した第２の時間フレームごとに、当該第２の時間フレームに前記第１の地点を出発して前記第２の地点に到達する第２のＯＤ交通量を算出し、
算出した前記第２のＯＤ交通量を用いて、前記第２の時間フレームごとに、前記リンクごとの交通量を推定する第２の交通シミュレーションを実行し、
推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第２のＯＤ交通量を調整しながら前記第２の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第２の時間フレームごとのＯＤ交通量を決定する、
処理をコンピュータが実行することを特徴とするＯＤ決定方法。
推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第１のＯＤ交通量を調整しながら前記第１の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第１の時間フレームごとのＯＤ交通量を決定し、
決定した前記第１の時間フレームごとのＯＤ交通量に基づいて、前記第１の時間フレームを分割した第２の時間フレームごとに、当該第２の時間フレームに前記第１の地点を出発して前記第２の地点に到達する第２のＯＤ交通量を算出し、
算出した前記第２のＯＤ交通量を用いて、前記第２の時間フレームごとに、前記リンクごとの交通量を推定する第２の交通シミュレーションを実行し、
推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第２のＯＤ交通量を調整しながら前記第２の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第２の時間フレームごとのＯＤ交通量を決定する、
処理をコンピュータが実行することを特徴とするＯＤ決定方法。
【0263】
（付記２）前記算出する処理は、
決定した前記第１の時間フレームごとのＯＤ交通量を、所定の比率に応じて前記第２の時間フレームごとに配分することにより、前記第２のＯＤ交通量を算出する、ことを特徴とする付記１に記載のＯＤ決定方法。
決定した前記第１の時間フレームごとのＯＤ交通量を、所定の比率に応じて前記第２の時間フレームごとに配分することにより、前記第２のＯＤ交通量を算出する、ことを特徴とする付記１に記載のＯＤ決定方法。
【0264】
（付記３）前記第１の地点と前記第２の地点とのペアと対応付けて、決定した前記第２の時間フレームごとのＯＤ交通量を出力する、
処理を前記コンピュータが実行することを特徴とする付記１に記載のＯＤ決定方法。
処理を前記コンピュータが実行することを特徴とする付記１に記載のＯＤ決定方法。
【0265】
（付記４）前記比率は、前記第１の時間フレームと前記第２の時間フレームとの時間幅をもとに設定される、ことを特徴とする付記２に記載のＯＤ決定方法。
【0266】
（付記５）前記比率は、前記経路に含まれる第１のリンクについて観測された、前記第２の時間フレームごとの交通量をもとに設定される、ことを特徴とする付記２に記載のＯＤ決定方法。
【0267】
（付記６）前記第１のリンクは、前記第１の地点に接続されたリンクである、ことを特徴とする付記５に記載のＯＤ決定方法。
【0268】
（付記７）前記第１の時間フレームごとのＯＤ交通量を決定する処理は、
前記第１の交通シミュレーションにより推定した前記リンクの交通量と、当該リンクについて観測された交通量との誤差を算出し、算出した前記誤差が許容範囲内であるか否かを判断し、
前記誤差が許容範囲内ではない場合、前記誤差に基づいて、前記第１のＯＤ交通量を更新し、更新した前記第１のＯＤ交通量を用いて、前記第１の交通シミュレーションを実行し、
算出した前記誤差が許容範囲内の場合、直前の前記第１の交通シミュレーションに用いた前記第１のＯＤ交通量を、前記第１の時間フレームごとのＯＤ交通量に決定する、
ことを特徴とする付記１に記載のＯＤ決定方法。
前記第１の交通シミュレーションにより推定した前記リンクの交通量と、当該リンクについて観測された交通量との誤差を算出し、算出した前記誤差が許容範囲内であるか否かを判断し、
前記誤差が許容範囲内ではない場合、前記誤差に基づいて、前記第１のＯＤ交通量を更新し、更新した前記第１のＯＤ交通量を用いて、前記第１の交通シミュレーションを実行し、
算出した前記誤差が許容範囲内の場合、直前の前記第１の交通シミュレーションに用いた前記第１のＯＤ交通量を、前記第１の時間フレームごとのＯＤ交通量に決定する、
ことを特徴とする付記１に記載のＯＤ決定方法。
【0269】
（付記８）前記第２の時間フレームごとのＯＤ交通量を決定する処理は、
前記第２の交通シミュレーションにより推定した前記リンクの交通量と、当該リンクについて観測された交通量との誤差を算出し、算出した前記誤差が許容範囲内であるか否かを判断し、
前記誤差が許容範囲内ではない場合、前記誤差に基づいて、前記第２のＯＤ交通量を更新し、更新した前記第２のＯＤ交通量を用いて、前記第２の交通シミュレーションを実行し、
算出した前記誤差が許容範囲内の場合、直前の前記第２の交通シミュレーションに用いた前記第２のＯＤ交通量を、前記第２の時間フレームごとのＯＤ交通量に決定する、
ことを特徴とする付記１に記載のＯＤ決定方法。
前記第２の交通シミュレーションにより推定した前記リンクの交通量と、当該リンクについて観測された交通量との誤差を算出し、算出した前記誤差が許容範囲内であるか否かを判断し、
前記誤差が許容範囲内ではない場合、前記誤差に基づいて、前記第２のＯＤ交通量を更新し、更新した前記第２のＯＤ交通量を用いて、前記第２の交通シミュレーションを実行し、
算出した前記誤差が許容範囲内の場合、直前の前記第２の交通シミュレーションに用いた前記第２のＯＤ交通量を、前記第２の時間フレームごとのＯＤ交通量に決定する、
ことを特徴とする付記１に記載のＯＤ決定方法。
【0270】
（付記９）対象期間を分割した第１の時間フレームごとに、当該第１の時間フレームに第１の地点を出発して第２の地点に到達する第１のＯＤ交通量を用いて、前記第１の地点から前記第２の地点までの経路に含まれるリンクごとの交通量を推定する第１の交通シミュレーションを実行し、
推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第１のＯＤ交通量を調整しながら前記第１の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第１の時間フレームごとのＯＤ交通量を決定し、
決定した前記第１の時間フレームごとのＯＤ交通量に基づいて、前記第１の時間フレームを分割した第２の時間フレームごとに、当該第２の時間フレームに前記第１の地点を出発して前記第２の地点に到達する第２のＯＤ交通量を算出し、
算出した前記第２のＯＤ交通量を用いて、前記第２の時間フレームごとに、前記リンクごとの交通量を推定する第２の交通シミュレーションを実行し、
推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第２のＯＤ交通量を調整しながら前記第２の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第２の時間フレームごとのＯＤ交通量を決定する、
処理をコンピュータに実行させることを特徴とするＯＤ決定プログラム。
推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第１のＯＤ交通量を調整しながら前記第１の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第１の時間フレームごとのＯＤ交通量を決定し、
決定した前記第１の時間フレームごとのＯＤ交通量に基づいて、前記第１の時間フレームを分割した第２の時間フレームごとに、当該第２の時間フレームに前記第１の地点を出発して前記第２の地点に到達する第２のＯＤ交通量を算出し、
算出した前記第２のＯＤ交通量を用いて、前記第２の時間フレームごとに、前記リンクごとの交通量を推定する第２の交通シミュレーションを実行し、
推定した前記リンクの交通量と当該リンクについて観測された交通量との誤差が小さくなるように、前記第２のＯＤ交通量を調整しながら前記第２の交通シミュレーションを繰り返し実行することにより、前記第１の地点と前記第２の地点とのペアについて、前記第２の時間フレームごとのＯＤ交通量を決定する、
処理をコンピュータに実行させることを特徴とするＯＤ決定プログラム。
【符号の説明】
【0271】
１０１ 情報処理装置
１１０ 第１のシミュレーション
１２０ 第２のシミュレーション
２００ 情報処理システム
２０１ ＯＤ推定装置
２０２ クライアント装置
２１０ ネットワーク
３００ バス
３０１ ＣＰＵ
３０２ メモリ
３０３ ディスクドライブ
３０４ ディスク
３０５ 通信Ｉ／Ｆ
３０６ 可搬型記録媒体Ｉ／Ｆ
３０７ 可搬型記録媒体
５００ 道路ネットワーク表
６００ 出発地点表
７００ 目的地点表
８００ ＯＤペア表
９００ 経路表
１０００ 制御部
１００１ 取得部
１００２ 設定部
１００３ 第１の実行部
１００４ 第１の決定部
１００５ 算出部
１００６ 第２の実行部
１００７ 第２の決定部
１００８ 出力部
１２０１ 初期ＯＤフロー
１２０２，２０００，２２０１，２２０２，３０００ ＯＤフロー
１２１０，２２１０ Ｆｏｒｗａｒｄ処理
１２２０，２２２０ Ｂａｃｋｗａｒｄ処理
１４００，２４００ 経路選択表
１５００，１６００，１８００，１９０１，２５００，２６００，２８００，２９０１ リンクフロー
１７００，２７００ 観測行列
１９００，２９００ 誤差
２００１，３００１ 補正値
１１０ 第１のシミュレーション
１２０ 第２のシミュレーション
２００ 情報処理システム
２０１ ＯＤ推定装置
２０２ クライアント装置
２１０ ネットワーク
３００ バス
３０１ ＣＰＵ
３０２ メモリ
３０３ ディスクドライブ
３０４ ディスク
３０５ 通信Ｉ／Ｆ
３０６ 可搬型記録媒体Ｉ／Ｆ
３０７ 可搬型記録媒体
５００ 道路ネットワーク表
６００ 出発地点表
７００ 目的地点表
８００ ＯＤペア表
９００ 経路表
１０００ 制御部
１００１ 取得部
１００２ 設定部
１００３ 第１の実行部
１００４ 第１の決定部
１００５ 算出部
１００６ 第２の実行部
１００７ 第２の決定部
１００８ 出力部
１２０１ 初期ＯＤフロー
１２０２，２０００，２２０１，２２０２，３０００ ＯＤフロー
１２１０，２２１０ Ｆｏｒｗａｒｄ処理
１２２０，２２２０ Ｂａｃｋｗａｒｄ処理
１４００，２４００ 経路選択表
１５００，１６００，１８００，１９０１，２５００，２６００，２８００，２９０１ リンクフロー
１７００，２７００ 観測行列
１９００，２９００ 誤差
２００１，３００１ 補正値
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23】
【図24】
【図25】
【図26】
【図27】
【図28】
【図29】
【図30】
【図31】
【図32】
【図33】