(57)【要約】【課題】乗降施設のような対象空間の出入口で計測した人数の情報から、対象空間における人の移動状況を高い精度で推定する技術に関する。【解決手段】人流動推定システムは、観測の対象となる対象空間の出入口を通過する人の数である流動人数を計測する計測部と、前記流動人数に基づき、所定の流動推定演算手法により、所定の旅客状態の人の数を推定した比較用推定人数と、前記対象空間での人の流動を推定した流動推定情報とを取得する推定部と、前記旅客状態の人の数を他の方法で取得した比較対象人数と前記比較用推定人数を比較し、比較結果に基づいて、前記流動推定情報の推定精度を評価する評価部と、を有している。【選択図】図２
2016045919.tif 000002
【特許請求の範囲】
【請求項１】
観測の対象となる対象空間の出入口を通過する人の数である流動人数を計測する計測部と、
前記流動人数に基づき、所定の流動推定演算手法により、所定の旅客状態の人の数を推定した比較用推定人数と、前記対象空間での人の流動を推定した流動推定情報とを取得する推定部と、
前記旅客状態の人の数を他の方法で取得した比較対象人数と前記比較用推定人数を比較し、比較結果に基づいて、前記流動推定情報の推定精度を評価する評価部と、
を有する人流動推定システム。
【請求項２】
前記計測部は、前記出入口で前記対象空間に流入する人の数である流入人数と、前記出入口で前記対象空間から流出する人の数である流出人数とを計測し、
前記推定部は、前記流入人数と前記流出人数のいずれか一方に基づき他方を推定することにより得られた結果を前記比較用推定人数とし、
前記評価部は、前記比較用推定人数と計測された前記他方とを比較することにより、前記流動推定情報の推定精度を評価する、
請求項１に記載の人流動推定システム。
【請求項３】
前記対象空間は、旅客が乗物に乗降する乗降施設内の空間であり、
前記出入口には、乗物の乗降場所に接続される第１出入口と、前記乗降施設の出入口である第２出入口とが含まれ、
前記人流動推定システムは、前記乗物の到着および出発に関する運行情報を保持する記憶部を更に有し、
前記推定部は、前記対象空間に流入する出入口と前記対象空間から流出する出入口との対に対応づけた、前記対象空間に流入して前記対象空間を移動して前記対象空間から流出する人の数を示すＯＤ情報を算出するとき、前記第１出入口からの流入に関する第１ＯＤ情報は前記乗物の到着に対応付けて算出し、前記第２出入口からの流入に関する第２ＯＤ情報は前記時間単位で算出する、
請求項２に記載の人流動推定システム。
【請求項４】
前記記憶部は、到着した乗物から前記第１出入口を通って前記対象空間に流入する人が他の各出入口に向かう第１分岐率の情報と、前記第２出入口から前記対象空間に流入する人が他の各出入口に向かう第２分岐率の情報とを更新可能に更に保持し、
前記推定部は、前記第１分岐率および前記第２分岐率に基づいて、前記第１ＯＤ情報と、前記第２ＯＤ情報と、前記第１ＯＤ情報および前記第２ＯＤ情報に基づき出入口間の移動時間を考慮して算出された、各出入口から流出する人の数である流出推定人数と、を推定し、
前記評価部は、前記流出推定人数を前記比較用推定人数として前記流出人数と比較することにより、前記第１ＯＤ情報および前記第２ＯＤ情報の推定精度を評価し、評価結果に基づいて前記第１分岐率と前記第２分岐率の少なくとも一方を更新するか否か判定する、
請求項３に記載の人流動推定システム。
【請求項５】
前記評価部は、前記評価結果において前記流出推定人数と前記流出人数の差分が最も大きい出入口の第１分岐率および第２分岐率を更新する、
請求項４に記載の人流動推定システム。
【請求項６】
前記計測部は、前記出入口で前記対象空間に流入する人の数である流入人数と、前記出入口で前記対象空間から流出する人の数である流出人数とを計測し、
前記推定部は、前記流入人数に基づき前記対象空間を格子状の単位空間に分割した格子空間における人の群集密度である第１群集密度を推定し、前記流出人数に基づき格子空間における人の群集密度である第２群集密度を推定し、
前記評価部は、前記第１群衆密度と前記第２群集密度を比較することにより、前記流動推定情報の推定精度を評価する、
請求項１に記載の人流動推定システム。
【請求項７】
前記推定部は、前記対象空間を格子状の単位空間に分割した格子空間を表わす空間情報と前記第１ＯＤ情報および前記第２ＯＤ情報とに基づき前記格子空間の単位空間に人を疑似するエージェントを出現させ、単位空間の間を移動させることにより、前記出入口間の移動時間を算出する、請求項４に記載の人流動推定システム。
【請求項８】
観測の対象となる対象空間の出入口を通る人を観測する観測システムを更に有し、
前記計測部は、前記観測システムの観測結果に基づいて、前記出入口を通過する人の数である流動人数を計測する、
請求項１に記載の人流動推定システム。
【請求項９】
計測手段が、観測の対象となる対象空間の出入口を通過する人の数である流動人数を計測し、
推定手段が、前記流動人数に基づき、所定の流動推定演算手法により、所定の旅客状態の人の数を推定した比較用推定人数と、前記対象空間での人の流動を推定した流動推定情報とを取得し、
評価手段が、前記旅客状態の人の数を他の方法で取得した比較対象人数と前記比較用推定人数を比較し、比較結果に基づいて、前記流動推定情報の推定精度を評価する、
人流動推定方法。
【請求項１０】
観測の対象となる対象空間の出入口を通過する人の数である流動人数を計測する手順と、
前記流動人数に基づき、所定の流動推定演算手法により、所定の旅客状態の人の数を推定した比較用推定人数と、前記対象空間での人の流動を推定した流動推定情報とを取得する手順と、
前記旅客状態の人の数を他の方法で取得した比較対象人数と前記比較用推定人数を比較し、比較結果に基づいて、前記流動推定情報の推定精度を評価する手順と、
をコンピュータに実行させるための人流動推定プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は空間内の人の流動を推定する技術に関する。
【背景技術】
【０００２】
大都市圏の鉄道駅や空港において通勤時や繁忙期の旅客の混雑が問題となっている。特に鉄道では、駅が混雑していると旅客の乗降に時間がかかって列車の運行が遅れたり、あるいは旅客が混雑したホームから線路に転落したりすることが懸念される。
【０００３】
鉄道駅や空港など乗降施設における旅客の混雑を効率的に解決していくためには、施設内での旅客の移動の状況を詳細に把握するのが好ましい。ここで把握すべき旅客の移動状況として、例えば鉄道駅においては、時間帯毎の改札から入場する旅客数および改札から出場する旅客数だけでなく、列車から列車へ乗り換える旅客数を含めた全ての駅利用者数、路線毎の旅客数、乗り換えのためにホーム間を移動する旅客数、列車毎の乗降する旅客数などが挙げられる。
【０００４】
旅客の移動を把握する手段として、レーザセンサを用いて通行人の移動を追跡し、行動解析を行う装置が特許文献１に示されている。
【０００５】
また、特定の空間における旅客数を算出する手段として、その空間の出入口を通過する人数をカメラにより計測し、空間内に存在する旅客の人数を計測する手段を備えた人数計測システムが特許文献２に開示されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００８−２１７２８９号公報
【特許文献２】特開２０１２−１９６９８７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかし、特許文献１の装置で把握できる旅客の移動は、レーザセンサのレーザ照射範囲までであり、駅構内全域を把握するには、階段などの高低差や柱や通路などの遮蔽物によってレーザが遮断され、届かないような場所がないように、多数のレーザセンサを設置する必要がある。しかし、レーザセンサを多数設置するにはコストがかさみ実現が難しい。
【０００８】
また、特許文献２のシステムでは推定したい空間全域をくまなくカメラで観測する必要はないが、対象空間内にいる旅客数を計測できるのみであり、乗り換えでホーム間を移動する旅客数など、その内訳は計測できない。
【０００９】
本発明の目的は、乗降施設のような対象空間の出入口で計測した人数の情報から、対象空間における人の移動状況を高い精度で推定する技術に関する。
【課題を解決するための手段】
【００１０】
本発明の一態様による人流動推定システムは、観測の対象となる対象空間の出入口を通過する人の数である流動人数を計測する計測部と、前記流動人数に基づき、所定の流動推定演算手法により、所定の旅客状態の人の数を推定した比較用推定人数と、前記対象空間での人の流動を推定した流動推定情報とを取得する推定部と、前記旅客状態の人の数を他の方法で取得した比較対象人数と前記比較用推定人数を比較し、比較結果に基づいて、前記流動推定情報の推定精度を評価する評価部と、を有している。
【発明の効果】
【００１１】
本発明によれば、出入口での計測に基づき流動推定演算手法を用いて対象空間での人の流動を推定するとともに、同じ流動推定演算手法で流動推定情報と比較用推定人数を推定し、比較用推定人数を他の手法で取得した比較対象人数と比較し、比較結果で流動推定情報の推定精度を評価するので、高い精度で流動推定情報を取得することができる。
【図面の簡単な説明】
【００１２】
【図１】本実施例における旅客流動推定装置１０００の構成の一例を示したブロック図である。
【図２】旅客流動推定装置における処理とデータの流れの概要を示した図である。
【図３】対象空間を示したレイアウト情報のデータの一例を示した図である。
【図４】出入口と、その出入口における対象空間の接続先との対応を示した出入口接続先対応表のデータの一例を示した図である。
【図５】出入口間の移動に要する所要時間を示す出入口間所要時間対応表のデータの一例を示した図である。
【図６】対象空間の出入口へ流入する流入人数情報１１０１のデータの一例を示す図である。
【図７】対象空間の出入口から流出する流出人数情報１１０２のデータの一例を示す図である。
【図８】公共交通の運行計画または運行実績を示す運行情報のデータの一例を示す図である。
【図９】降車人数比情報のデータの一例を示した図である。
【図１０】イベントにより発生する時刻以降に集中して旅客が流入する特性を有する出入口における流入人数を表すイベントＯＤデータの一例を示した図である。
【図１１】時間的な偏りが少なく、比較的定常的に旅客が流入する出入口における旅客の流入人数を表す定常ＯＤデータの一例を示した図である。
【図１２】各列車からの降車旅客の移動先となる各出入口への分岐の割合を示した降車旅客分岐情報のデータの一例を示した図である。
【図１３】時間的な偏りが少ない出入口から流入する旅客の他の各出入口への分岐率を示した流入分岐情報の一例を示した図である。
【図１４】推定部１２０１で推定され出力される流出人数情報の一例を示した図である。
【図１５】定常ＯＤデータの初期値の一例を示す図である。
【図１６】イベントＯＤデータの初期値の一例を示す図である。
【図１７】流出人数情報の初期値の一例を示す図である。
【図１８】推定部１２０１の全体処理の概要を説明するための図である。
【図１９】ＯＤ推定部１２１１が実行するＯＤ推定処理のフローチャートである。
【図２０】定常ＯＤ情報の推定処理を示すフローチャートである。
【図２１】ホームａの列車発着と、ホームａと接続する出入口Ａでの計測結果との関係の一例を示した図である。
【図２２】列車割当を含むイベントＯＤ推定処理のフローチャートである。
【図２３】移動推定部１２１２の処理のフローチャートである。
【図２４】格子空間上に設定される歩行者エージントの移動方向の一例を示す図である。
【図２５】移動する歩行者エージェント間の相互作用について説明するための図である。
【図２６】評価部１２０３での評価の一例について説明するための図である。
【図２７】評価部１２０３における分岐情報更新処理を示すフローチャートである。
【図２８】群集密度表示端末の画面の一例を示す図である。
【発明を実施するための形態】
【００１３】
本発明の一実施例について図面を用いて説明する。
＜概要＞
【００１４】
まず、本発明における旅客流動推定の一例について概要を説明する。
【００１５】
一般に、任意の空間を対象空間として設定したとき、対象空間の出入口（境界）を通過して対象空間内に流入した人は、対象空間内を移動し、いずれかの出入口（境界）を通過して対象空間外に流出する。前記対象空間内に流入するすべての旅客は基本的に対象空間内を移動し、いずれかの出口から流出することは自明であるが、どの出入口から流出するかを推定する手段を本発明では実現する。
【００１６】
本発明では、旅客流動を推定したい任意の対象空間に対して、出入口の流入人数と流出人数を計測し、対象空間内の移動を推定する。対象空間内の移動を推するにあたっては、まず前記流出データを用いず、流入データをもとに移動を推定する。そして、その推定の結果と、推定の演算に入力に含んでいない流出データとを比較することで、推定精度を評価し、この評価結果を改善するように、パラメータを更新しながら推定演算を繰り返すことで、推定精度を向上させる。
【００１７】
図１は、本実施例における旅客流動推定装置１０００の構成の一例を示したブロック図である。
【００１８】
旅客流動推定装置１０００は、流入人数情報１１０１、流出人数情報１１０２、空間情報１１０３、運行情報１１０４、分岐情報１１０６、および評価結果情報１１０７を少なくとも記録する記録部１１００と、推定部１２０１および評価部１２０３を少なくとも有する演算部１２００と、対象空間の出入口を通る人を観測する観測システムと接続する人数計測部１３００と、出力部１００１とを有している。
【００１９】
観測システムは、カメラシステム１３１０と自動改札機システム１３２０のいずれか一方または両方である。前記記録部１１００、演算部１２００、人数計測部１３００および出力部１００１はそれぞれ相互に接続している。記録部１１００、演算部１２００、および人数計測部１３００は、それぞれネットワークを介して設置される別個の装置に実装されてもよい。
【００２０】
図２は、旅客流動推定装置における処理とデータの流れの概要を示した図である。以下、図２を用いて、処理とデータの流れの概要を説明する。
【００２１】
まず、人数計測部１３００が対象空間の出入口の通過人数を流入人数と流出人数に分けて計測し、流入人数情報１１０１と流出人数情報１１０２を出力する。
【００２２】
次に推定部１２０１は、計測した流入人数情報１１０１、空間情報１１０３、運行情報１１０４、分岐情報１１０６のデータを入力として、所定の流動推定演算手法により、対象空間内の人の移動を推定する。また、推定部１２０１は、比較評価のために、同じ流動推定演算手法で算出した推定結果から出入口の流出人数（流出推定人数）を算出し、評価部１２０３に出力する。
【００２３】
評価部１２０３は、人数計測部１３００で計測した流出人数情報１１０２に示された流出人数と、推定部１２０１で推定した流出人数を入力とし、推定部１２０１で推定した流出人数の精度を評価する。ここで、評価部１２０３は、評価結果が予め設定しておいた一定の条件を満たしていれば、推定部１２０１から出力された推定結果を、評価結果と対応させて評価結果情報１２０７として出力し、評価結果がその条件を満たしていなければ、評価結果に応じて、推定部１２０１の入力パラメータのひとつである分岐情報１１０６を更新し、処理を推定部１２０１に戻し、再度対象空間内の移動を推定し直すという処理を繰り返させる（矢印１４００）。
【００２４】
上記評価に用いる一定の条件の例としては、流出人数同士を比較した評価の結果、推定部１２０１で推定した流出人数情報の精度が任意の精度以上であること、推定演算の繰り返し回数が一定回数以下であること、推定演算に要した計算時間が予め設定しておいた時間以下であること、などが挙げられる。また、これらのうち複数の条件を組み合わせて評価に用いてもよい。
【００２５】
最後に、出力部１００１が評価結果情報１２０７を出力部１００１から出力する。
【００２６】
以上で、図２を用いた処理とデータの流れの概要についての説明を終える。
＜構成＞
【００２７】
以下、図１を用いて旅客流動推定装置１０００の構成について説明する。
【００２８】
旅客流動推定装置１０００では、人数計測部１３００はカメラシステム１３１０に接続されたカメラ１３１１、１３１３、１３１５と、自動改札機システム１３２０に接続された自動改札機１３２１とにより、対象空間１３０１の出入口１３１２、１３１４、１３１６、１３２２から前記出入口の流入人数情報１１０１および流出人数情報１１０２を計測し、演算部１２００が有する推定部１２０１にて旅客の移動を推定し、前記推定の結果を評価部１２０３で評価し、推定結果を記録部１１００へ保存し、出力部１００１は記録部１１００に保存された推定結果の情報を出力する。
【００２９】
なお、評価部１２０３での評価結果で、予め定めた条件を満たさない場合、演算処理に入力されるパラメータを更新して前記条件を満たすまで推定部１２０１での旅客の移動の推定を繰り返してもよい。前記推定の繰り返しにあたり、推定結果を変化させることを目的として、分岐情報１１０６に記録されるパラメータ情報を更新する。各部での処理やデータの詳細は後述する。
【００３０】
以上で、図１を用いた旅客流動推定装置１０００の構成についての説明を終える。
＜データの説明＞
【００３１】
本実施例で用いる各種データについて説明する。
【００３２】
図３は、対象空間を示したレイアウト情報のデータの一例を示した図である。レイアウト情報は、記録部１１００の空間情報１１０３にあらかじめ記録されている。
【００３３】
本実施例では、対象空間を含む空間を格子状に分割した格子空間８０００としてレイアウト情報を表現する。格子空間８０００は通路８０１１、壁８０１２、階段８０１３、改札８０１４、および出入口８０１５の単位空間を組み合わせて構成している。そのうち、壁８０１２以外は通行可能である。格子空間８０００中では、位置８０２１は通路８０１１、位置８０２２が壁８０１２、位置８０２３が階段８０１３、位置８０２４が改札８０１４、位置８０２５から８０２７が対象空間への出入口８０１５の単位空間となっている。
【００３４】
本実施例のデータの例では、位置８０２５とその周囲に隣接する位置の出入口８０１５をまとめて出入口Ａ、位置８０２６とその周囲に隣接する位置の出入口８０１５をまとめて出入口Ｂ、位置８０２７とその周囲に隣接する位置の出入口８０１５をまとめて出入口Ｃ、位置８０２８とその周囲に隣接する位置の出入口８０１５をまとめて出入口Ｄとしている。推定部１２０１はこれらの出入口からの流出人数を演算により推定し、人数計測部１３００はこの出入口に対応する現実の空間において流入人数及び流出人数を計測する。
【００３５】
図４は、出入口と、その出入口における対象空間の接続先との対応を示した出入口接続先対応表のデータの一例を示した図である。出入口接続先対応表は記録部１１００の空間情報１１０３にあらかじめ記録されている。出入口接続先対応表には、出入口９０００、その接続先９００１、出入口から接続先までの平均移動時間９００３が記録されている。出入口とその接続先までの平均移動時間は、例えば、接続先がホームの場合、列車を降車した後、対象空間の出入口に対応する計測位置に到達するまでに要する平均移動時間とする。
【００３６】
例えば、レコード９０１１は、出入口Ａが１番線ホームに接続し、出入口Ａから１番線ホームまでの平均移動時間は20秒であることを示している。レコード９０１２は、出入口Ｃが２番線ホームに接続し、出入口Ｃから２番線ホームまでの平均移動時間は25秒であることを示す。本実施例ではテーブルとしてデータを定義したが、同様の情報を算出できるデータと処理によって代替してもよい。
【００３７】
図５は、出入口間の移動に要する所要時間を示す出入口間所要時間対応表のデータの一例を示した図である。出入口間所要時間対応表は記録部１１００の空間情報１１０３にあらかじめ記録されている。
【００３８】
縦方向１０００１には対象空間に流入する出入口、横方向１０００２は対象空間から流出する出入口の項目が示されている。対象空間に流入する出入口である移動開始地点と、対象空間から流出する出入口である移動終了地点との組み合わせに対する、計測位置間の平均移動時間を示す。例えば、１００１０は出入口Ｄから流入し、出入口Ｂから流出するまでの移動平均時間が50秒であることを示す。出入口間の平均移動時間は、図３に示したレイアウト情報からの算出により決定してもよい。
【００３９】
図６、図７は人数計測部１３００で計測されるデータを示している。図６は対象空間の出入口へ流入する流入人数情報１１０１のデータの一例を示す図である。図７は対象空間の出入口から流出する流出人数情報１１０２のデータの一例を示す図である。
【００４０】
流入人数情報１１０１のデータは、計測期間３００１と、各出入口の計測位置３００２〜３００５とによって構成される。流出人数情報１１０２のデータは、計測期間４００１と、各出入口の計測位置４００２〜４００５とによって構成される。データの項目自体は流入人数情報１１０１と「流出人数情報１１０２は同じである。出入口を通過する人の進行方向が対象空間に対して流入方向であれば流入人数として計測し、流出方向であれば流出人数として計測する。
【００４１】
図８は、公共交通の運行計画または運行実績を示す運行情報のデータの一例を示す図である。運行情報は記録部１１００の運行情報１１０４にあらかじめ記録されている。運行情報は、列車を一意に特定できる列車ＩＤ５００１、列車が到着する番線５００２、列車の停車駅を示す列車種別５００３、駅に到着する列車の到着時刻５００４、および駅を出発する列車の出発時刻５００５により構成される。
【００４２】
図９は、降車人数比情報のデータの一例を示した図である。降車人数比情報は、運行情報に示された各列車に対応する列車ＩＤ６００１と、降車人数比６００２により構成される。降車人数比情報は、記録部１１００の分岐情報１１０６にあらかじめ記録されている。
【００４３】
降車人数比は、各列車から同一のホームに降車する旅客の人数の比を示す。降車人数比は、複数の列車が連続して到着し、同一のホームに旅客が降車した場合に利用される。例えば、ＩＤ＝４とＩＤ＝５の列車が同時に到着した場合、ＩＤ＝４の列車からの降車人数とＩＤ＝５の列車からの降車人数の比は３：２となることを示している。
【００４４】
図１０は、イベントにより発生する時刻以降に集中して旅客が流入する特性を有する出入口における流入人数を表すイベントＯＤデータの一例を示した図である。イベントＯＤデータは、例えば、ホームに接続する階段を対象空間の出入口とする場合など、列車の発着といったイベントに起因して対象空間へまとまって流入する旅客の人数を表すデータである。イベントＯＤデータは、推定部１２０１の推定結果のひとつとして出力され、評価結果情報にも記録される。イベントＯＤデータは、流入する出入口１４００１、流入が開始される時刻（流入時刻）１４００２、移動先出口別出現人数１４００３、１４００４、１４００５、１４００６により構成される。流入時刻１４００２は、例えば列車の到着時刻である。
【００４５】
図１１は、時間的な偏りが少なく、比較的定常的に旅客が流入する出入口における旅客の流入人数を表す定常ＯＤデータの一例を示した図である。定常ＯＤデータは、改札口や通路からの流入旅客など平均的に対象空間へ流入する旅客を表すデータである。定常ＯＤデータは、推定部１２０１の推定結果のひとつとして出力され、評価結果情報１１０７にも記録される。流入する出入口１３００１、流入する期間（流入期間）１３００２、移動先出口別出現人数１３００３、１３００４、１３００５、１３００６により構成される。
【００４６】
図１２は、各列車からの降車旅客の移動先となる各出入口への分岐の割合を示した降車旅客分岐情報のデータの一例を示した図である。降車旅客分岐情報は記録部１１００の分岐情報１１０６に記録され、推定部１２０１では入力の一つとしてイベントＯＤデータの推定に用いられる。降車旅客分岐情報は運行情報に対応する列車ＩＤ１１００１と、該当列車の降車旅客が移動先（流出する出入口）の割合を移動先毎に設定した降車旅客分岐率１１００２〜１１００５により構成される。例えば、データ１１０１１は、列車ＩＤが６の列車（運行情報よりホームａに7:34:53到着、7:35:27出発する列車(図８参照)）から降車し、出入口Ｃへ向かう人数が34%であることを示す。降車旅客分岐情報は演算部１２００での処理の中で更新されるデータである。降車旅客分岐情報の初期値を任意に設定して分岐情報１１０６に記録しておく。初期値の設定は例えば均等な割合を設定するなどでよい。
【００４７】
図１３は、時間的な偏りが少ない出入口から流入する旅客の他の各出入口への分岐率を示した流入分岐情報の一例を示した図である。流入分岐情報は記録部１１００の分岐情報１１０６に記録され、推定部１２０１では入力の一つとして定常ＯＤデータの推定に用いられる。流入分岐情報は、対象期間１２００１、流入位置１２００２、流出位置毎に設定した出入口から流入後の移動先の割合１２００３〜１２００５により構成される。例えば、データ１２０１１は、7:35から7:36の間に出入口Ｄからの流入する旅客のうち、出入口Ｂへ移動する割合は35%であることを示す。
【００４８】
図１４は、推定部１２０１で推定され出力される流出人数情報の一例を示した図である。流出人数情報は、対象期間１５００１毎に各出入口１５００２、１５００３、１５００４、１５００５を通過すると推定部１２０１が推定した流出人数（比較用推定人数）を示す情報である。そのデータの形式は、人数計測部１３００で計測する流出人数情報と同様である。評価部１２０３では、この推定した流出人数情報と、人数計測部１３００で計測した流出人数情報とを比較して、推定精度の評価を行う。
【００４９】
図１５は、定常ＯＤデータの初期値の一例を示す図である。図１６は、イベントＯＤデータの初期値の一例を示す図である。図１７は、流出人数情報の初期値の一例を示す図である。
【００５０】
本実施例で用いるデータについての説明は以上である。
＜演算処理＞
【００５１】
演算部１２００での処理について説明する。以下の説明ではステップをＳと省略する。
【００５２】
図１８は、推定部１２０１の全体処理の概要を説明するための図である。まず、図１８を用いて推定部の概要について説明する。
【００５３】
推定部１２０１は、ＯＤ推定部１２１１と移動推定部１２１２を有し、２段階で移動を推定する構成とする。ＯＤ推定部１２１１は、流入人数情報１１０１、空間情報１１０３、運行情報１１０４、および分岐情報１１０６を入力とし、ＯＤ情報１２２１を推定し、移動推定部１２１２へ出力する。
【００５４】
移動推定部１２１２は、ＯＤ情報１２２１と、分岐情報１１０６を入力とし、旅客の移動を推定し、その推定結果（移動推定結果）を出力する。移動推定結果には少なくとも流出人数情報を含んでおり、更に、入力としたＯＤ情報１２２１を含んでもよい。
【００５５】
以上で、図１８を用いた推定部１２０１の概要の説明を終える。
【００５６】
以下では、まずＯＤ推定部１２１１の処理について説明する。
【００５７】
図１９は、ＯＤ推定部１２１１が実行するＯＤ推定処理のフローチャートである。以下、図１９を用いてＯＤ推定部の処理について説明する。
【００５８】
Ｓ３１０１では、ＯＤ推定部１２１１は、すべての出入口に対してそれぞれＳ３１０５に到達するまでの処理を繰り返し実行する。以下、処理の対象とする出入口をｅとする。
【００５９】
Ｓ３１０２では、出入口ｅの接続先がホームであるか否か判定する。出入口ｅの接続先がホームであるかどうかの判断は、図４で説明した出入口接続先対応表を参照することで行う。
【００６０】
ＯＤ推定部１２１１は、出入口ｅの接続先がホームであるときＳ３１０３へ進み、出入口ｅの接続先がホームではないときＳ３１０４へ進む。
【００６１】
Ｓ３１０３では、ＯＤ推定部１２１１は、流入人数情報のうち、出入口ｅに対応する流入人数を入力とし、イベントＯＤ情報の推定を行う。イベントＯＤ情報の推定については後述する。
【００６２】
Ｓ３１０４では、ＯＤ推定部１２１１は、流入人数情報のうち、出入口ｅに対応する流入人数を入力とし、定常ＯＤ情報の推定を行う。定常ＯＤ情報の推定は、前記流入人数と、分岐情報１１０６の流入分岐情報から算出する。定常ＯＤ情報の推定については後述する。
【００６３】
Ｓ３１０３あるいはＳ３１０４の処理の後、ＯＤ推定部１２１１は、Ｓ３１０５へ進み、すべての出入口に対して上記処理を行うまで、Ｓ３１０１に戻り処理を繰り返す。すべての出入口に対して処理を完了すると、ＯＤ推定部１２１１は、Ｓ３１０５へ進み、各出入口に対して推定したＯＤ情報を出力し、ＯＤ推定処理は終了する。
【００６４】
図１９を用いたＯＤ推定処理についての説明は以上である。
【００６５】
図１７は、定常ＯＤ情報の推定処理を示すフローチャートである。図２０を用いて定常ＯＤ情報の推定処理について以下で説明する。
【００６６】
Ｓ２７０１のループで、ＯＤ推定部１２１１は、流入人数情報を時系列に走査し、処理する。以下、処理の対象とする期間をｔとする。
【００６７】
Ｓ２７０２では、ＯＤ推定部１２１１は、流入人数情報の期間ｔに出入口ｅから流入人数を取得し、ｎとする。
【００６８】
Ｓ２７０３のループで、ＯＤ推定部１２１１は、流出先の出入口に対して走査し、処理する。以下、処理の対象とする流出先の出入口ｅ’とする。
【００６９】
Ｓ２７０４では、ＯＤ推定部１２１１は、流入分岐情報の流入位置が出入口ｅかつ、期間がｔであるレコードの出入口ｅ’の流出割合を取得しｒとする。
【００７０】
Ｓ２７０５では、ＯＤ推定部１２１１は、推定する定常ＯＤ情報において期間ｔにおいて出入口ｅから出入口ｅ’への移動人数をｒ×ｎとして設定する。
【００７１】
Ｓ２７０６では、ＯＤ推定部１２１１は、すべての流出先ｅ’の出入口に対して処理が完了するまで、Ｓ２７０３へ戻って繰返し処理し、処理が完了するとＳ２７０７へ進む。
【００７２】
Ｓ２７０７では、ＯＤ推定部１２１１は、流入人数情報に定義されているすべての期間ｔに対して処理を行うまで、Ｓ２７０１へ戻って繰返し処理し、処理が完了するとＳ２７０８へ進む。
【００７３】
Ｓ２７０８では、ＯＤ推定部１２１１は、算出した定常ＯＤ情報を出力する。
【００７４】
以上で、図２０を用いた定常ＯＤ情報の推定処理の説明を終える。
【００７５】
図２１は、ホームａの列車発着と、ホームａと接続する出入口Ａでの計測結果との関係の一例を示した図である。イベントＯＤ情報の推定では、ある期間に計測された流入人数がどの列車から降車してきた人であるかを推定し、推定した列車に対して計測した流入人数を列車に割り当てる列車割当の処理が必要である。
【００７６】
以下、図２１を用いてイベントＯＤ推定における列車割当の処理について説明する。図中のグラフの横軸２１００１は流入人数の計測時刻、縦軸２１００２は計測された流入人数を示している。２１０２１から２１０２６は、計測期間毎の流入人数を示している。例えば、２１０２１は7:30から7:31の期間に計測された流入人数に対応する。
【００７７】
太線２１０１１から２１０１３は列車の発着を表す線であり、この線では縦軸２１００２は列車の位置を示しており、この線と横軸２１００１との交差点のうち、時刻の早い左側の交点は駅への到着時刻を示し、時刻の遅い右側の交点は出発時刻を表す。
【００７８】
流入人数として計測されるタイミングは、列車の到着時刻以後であるため、流入人数の計測期間以前に到着した列車のうち、到着時刻が計測期間から一定時間以内であるという条件を満たす列車に流入人数を割り当てる。
【００７９】
到着時刻が計測期間から一定時間以内である列車が複数ある場合、図９に示した、運行情報１１０４の降車人数比情報を用いて各列車から降車した旅客の割合を算出し、その割合に従って各列車に流入人数を割り当てる。なお、配分しようとする流入人数に対して、それ以前に計測された流入人数の配分結果を考慮することとする。
【００８０】
例えば、7:35の流入人数２１０２５は、列車２１０１２（列車ＩＤが４）および列車２１０１３（列車ＩＤが６）のどちらから降車してきたかの推定が必要である。列車２１０１２から降車した旅客の人数と、列車２１０１３から降車した旅客の人数の比率は、図９に示された３：５である。
【００８１】
列車２１０１２から降車した旅客のうち２３人は7:34に計測されている。列車２１０１３から降車した旅客のうち２人は7:36に計測されている。列車２１０１２から降車し7:35に計測された旅客の人数をｘとすると、以下の関係式が成り立つ。
【００８２】
（２３＋ｘ）：（８０−ｘ＋２）＝３：５
【００８３】
この式を解くと、ｘ≒１６人となり、7:35の流入人数のうち、列車２１０１２に配分される人数は１６人、列車２１０１２に配分される人数は６４人と求めることができる。
【００８４】
以上で、図２１を用いたイベントＯＤ推定の列車割当の説明を終える。
【００８５】
図２２は、上述した列車割当を含むイベントＯＤ推定処理のフローチャートである。以下、図２２を用いてイベントＯＤ推定の処理について説明する。
【００８６】
Ｓ２２０１のループで、ＯＤ推定部１２１１は、流入人数情報を時系列に走査し処理する。以下、処理の対象とする期間をｔとする。
【００８７】
Ｓ２２０２では、ＯＤ推定部１２１１は、割当対象となる列車を取得し、割当対象となる列車の集合Ｃを取得する。割当対象となる列車は、計測期間以前に到着した列車のうち、あらかじめ設定しておいた時間範囲内に到着する列車であり、その情報は運行情報１１０４を探索して取得することができる。ただし、この取得方法で対象となる列車が一件も取得できない場合、計測期間以前に到着した列車のうち、到着時刻が計測期間にもっとも近い列車を取得することとする。
【００８８】
Ｓ２２０３のループで、ＯＤ推定部１２１１は、割当対象の列車集合Ｃに対して走査し処理する。以下、処理の対象とする列車を列車ｃとする。
【００８９】
Ｓ２２０４では、ＯＤ推定部１２１１は、流入人数情報から期間ｔの流入人数を取得し、前記の列車割当により列車ｃからの流入人数ｎを推定する。
【００９０】
Ｓ２２０５のループで、ＯＤ推定部１２１１は、流出先の各出入口に対して走査し処理する。以下、処理の対象とする流出先の出入口ｅ’とする。
【００９１】
Ｓ２２０６では、ＯＤ推定部１２１１は、降車旅客分岐情報の列車ｃの列車ＩＤに対応するレコードの出入口ｅ’の流出割合を取得し、ｒとする。
【００９２】
Ｓ２２０７では、ＯＤ推定部１２１１は、推定するイベントＯＤ情報において、列車ｃの降車旅客のうち、出入口ｅ’への移動人数をr×nとして設定する。
【００９３】
Ｓ２２０８では、ＯＤ推定部１２１１は、すべての流出先の出入口に対して処理が完了するまで、Ｓ２２０５へ戻って繰返し処理し、処理が完了するとＳ２２０９へ進む。
【００９４】
Ｓ２２０９では、ＯＤ推定部１２１１は、割当対象列車集合Ｃのすべての列車に対して、処理を行うまで、Ｓ２２０３へ戻って繰返し処理し、処理が完了するとＳ２２１０へ進む。
【００９５】
Ｓ２２１０では、流入人数情報に定義されているすべての期間に対して処理を行うまで、Ｓ２２０１へ戻って繰返し処理し、処理が完了するとＳ２２１１へ進む。
【００９６】
Ｓ２２１１では、ＯＤ推定部１２１１は、算出したイベントＯＤ情報を出力する。
【００９７】
以上で、図２２を用いたイベントＯＤ情報の推定処理の説明を終える。
【００９８】
ＯＤ推定部１２１１の処理についての説明も以上である。
【００９９】
次に、移動推定部１２１２の処理の一例について説明する。
【０１００】
移動推定部１２１２では、図３で説明したレイアウト情報により表現される格子空間を対象空間に見立て、計算上の時刻を小刻みにΔtずつ進めていき、時刻毎に対象空間内を移動する人を推定する。前記Δtは本実施例では数ミリ秒程度として設定するが、要求される計算時間や精度によって変更してもよい。
【０１０１】
このような推定手法はエージェントベースシミュレーションとして知られる公知の技術であり、現実世界で対象空間内を移動する歩行者をシミュレーション上の格子空間上で表したものを以下では歩行者エージェントとする。歩行者エージェントは、少なくとも格子空間上の位置、移動先の目的地および一意に特定できるＩＤを有し、それぞれ前記移動先の出入口に向かって移動することとする。格子空間内を移動する人は、レイアウト情報内に定義された出入口からの流入として格子空間に出現し、移動先の出入口まで移動し、移動先の出入口に到達すると流出として格子空間から消失する。
【０１０２】
図２３は、移動推定部１２１２の処理のフローチャートである。図２３を用いて移動推定部１２１２の処理の流れについて以下で説明する。
【０１０３】
Ｓ２５０１では、移動推定部１２１２は、計算上の時刻ｔを移動推定の開始時刻に初期化する。移動推定の開始時刻は、入力とするＯＤデータのもっとも早い流入時刻以前に設定すればよい。
【０１０４】
Ｓ２５０２では、移動推定部１２１２は、歩行者エージェントの流入および流出の処理を行う。流入は、各出入口に対応したイベントＯＤ情報および定常ＯＤ情報に基づき決定する。
【０１０５】
流入処理は、出入口毎に行う。出入口がホームに接続した出入口か、ホームに接続していない出入口かによって流入処理は異なる。以下、それぞれについて説明する。
【０１０６】
対象とする出入口がホームに接続しない出入口の場合、移動推定部１２１２は、定常ＯＤ情報に基づいて流入処理を行う。定常ＯＤ情報の対象となる期間内に、移動先別出現人数に設定されている人数に従い、前記移動先に向かう歩行者エージェントを確率的に時間的な偏りない人数分だけ流入として出現させる。
【０１０７】
対象とする出入口がホームに接続する出入口の場合、移動推定部１２１２は、イベントＯＤ情報に基づいて流入処理を行う。イベントＯＤ情報の流入時刻以降、各出入口に移動先別出現人数に従い、歩行者エージェントを流入として出現させる。イベントＯＤ情報により出現する歩行者エージェントは、流入時刻の直後に集中して出現することになる。
【０１０８】
流入処理において、移動推定部１２１２は、出現した歩行者エージェントを、空間内に存在する歩行者エージェント集合Ａに追加する。
【０１０９】
流入処理についての説明は以上とし、次に流出処理の説明に移る。
【０１１０】
流出処理は、出入口としてレイアウト情報に設定されている格子空間上に存在する歩行者エージェントのうち、前記歩行者エージェントの目的地と、歩行者エージェントが存在している格子空間に設定させている出入口とが一致する場合に、流出として歩行者エージェントを消失させるという処理である。流出として消失させる歩行者エージェントは空間内に存在する歩行者エージェント集合Ａからも削除する。
【０１１１】
Ｓ２５０３のループで、移動推定部１２１２は、格子空間内に存在するすべての歩行者エージェントに対して走査し処理する。以下、処理の対象とする歩行者エージェントをａとする。
【０１１２】
Ｓ２５０４では、移動推定部１２１２は、歩行者エージェントａの移動処理を行う。歩行者エージェントａの移動処理では、一定時間間隔以上、格子空間内で移動していない歩行者エージェントのみを移動させる。歩行者エージェントの移動方向は、格子空間上設定されている各出入口への移動方向に従い、歩行者エージェントの移動先の出入口に対応した移動方向に隣接する格子空間へと移動させる。
【０１１３】
図２４は、格子空間上に設定される歩行者エージントの移動方向の一例を示す図である。図２４を用いて移動方向について説明する。矢印２３００２、２３００３はその向きによって、格子空間上に設定されるエージェントの移動方向を示している。図２４の例では、出入口Ｂを移動先とする歩行者エージェントの移動方向を示したものである。矢印２３００２は下方向に直線的に進めば出入口Ｂに到着できるので矢印は下方向に設定される。矢印２３００３は、その位置からでは壁を回り込まなければ出口Ｂに到着できないため、回り込む方向に矢印が設定されている。この矢印の向きの設定は任意の方法で行ってよい。例えば、人手により予め決定しておくことや、公知の技術であるDijkstraのアルゴリズムを用いて出入口まで移動可能な格子空間のみを移動した場合の最短距離経路をレイアウト情報から自動的に算出し、その最短経路を示すような矢印を設定させてもよい。
【０１１４】
図２４を用いた格子空間上に設定される移動方向についての説明を終える。
【０１１５】
図２５は、移動する歩行者エージェント間の相互作用について説明するための図である。
【０１１６】
空間図２４００１は、時刻ｔの時点での歩行者エージェントの格子空間上の位置を示している。空間図２４００２は時刻ｔ＋Δｔの時点での歩行者エージェントの格子空間上の位置を示している。
【０１１７】
２４０１０、２４０２０、２４０３０、２４０４０の丸は歩行者エージェントを示している。空間図２４００１と空間図２４００２を比べると、歩行者エージェント２４０１０、２４０３０、２４０４０は時刻ｔから時刻ｔ＋Δｔの間では移動しておらず、歩行者エージェント２４０２０のみが移動していることが分かる。
【０１１８】
本実施例の移動推定処理では、１つの格子空間には１人歩行者エージェントのみが存在できることとし、歩行者エージェントは進もうとする方向に他の歩行者エージェントが既に存在する場合、移動させずに待機させる。
【０１１９】
例えば、時刻ｔ時点で、歩行者エージェント２４０３０および歩行者エージェント２４０４０が位置する格子空間上の単位空間の矢印は、歩行者エージェント２４０２０が位置する単位空間を向いている。そのため、時刻ｔ時点の歩行者エージェント２４０３０および２４０４０は、移動させることができない。ただし、時刻ｔ＋Δｔ時点では歩行者エージェント２４０２０が移動しているため、歩行者エージェント２４０３０および２４０４０は移動可能となる。ただし、どちらかの歩行者エージェントが単位空間２４１０１へ移動すると他方の歩行者エージェントは移動できなくなり、単位空間２４１０１に歩行者エージェントが存在しなくなるまで待つことになる。
【０１２０】
以上で図２５を用いた歩行者エージェント間の相互作用についての説明を終え、図２３の移動推定部の処理のフローチャートの説明に戻る。
【０１２１】
Ｓ２５０５では、格子空間内に存在するすべての歩行者エージェントに対して、時刻ｔ時点で移動処理を終えるまでＳ２５０３へ戻り繰返し処理し、処理が完了するとＳ２５０６へ進む。
【０１２２】
Ｓ２５０６では、時刻ｔをΔｔ進める。
【０１２３】
Ｓ２５０７では、時刻ｔが終了時刻以前であれば、Ｓ２５０２へ戻り処理を繰り返し、時刻ｔが終了時刻となると、移動推定部の処理を終了させる。
【０１２４】
Ｓ２５０６において、推定した歩行者エージェントの位置は、逐次記録しておき、終了時にまとめてデータや動画として出力してもよい。また記録せずにＳ２５０６の処理と同時に画面を持つ端末に出力し表示してもよい。
【０１２５】
以上で、図２３を用いた移動推定部１２１２の処理の流れについての説明を終える。
【０１２６】
図２３で説明した移動推定部の処理は一例であり、例えば、他にも出入口間の平均移動時間に乱数を加えた時間で出入口間を移動するとして移動を推定するなど他の方法も考えられる。
【０１２７】
移動推定部１２１２の処理についての説明も以上である。
【０１２８】
次に、評価部１２０３の処理について説明する。
【０１２９】
評価部１２０３は、図２で示した通り、計測部１３００で計測された流出人数情報１１０２と推定部１２０１から出力される推定結果とを入力とし、ＯＤ情報の推定精度を評価し、評価結果に基づき、評価結果情報１２０７へ出力するか、あるいは、あらかじめ定めた繰り返し条件に従い、分岐情報を更新して再度推定処理をやり直すかを判断する。
【０１３０】
降車旅客分岐率情報と流入分岐率情報の初期値は、一例として均等な割合が設定されているとすると、ＯＤ推定部１２１１で算出される定常ＯＤ情報およびイベントＯＤ情報は図１５および図１６にそれぞれ示されるように各移動先に均等に旅客が出現するというデータとなる。その分岐情報の初期値によって算出された定常ＯＤ情報およびイベントＯＤ情報に基づいて算出される流出人数情報が例えば図１７に示されているものであるとする。
【０１３１】
図２６は、評価部１２０３での評価の一例について説明するための図である。以下、図２６を用いて、評価部１２０３での評価の処理について説明する。
【０１３２】
評価部１２０３では、人数計測部１３００で計測した流出人数の計測期間ごとに評価を行う。図２６は、一例として、７：３０から７：３１の期間における流出人数の推定処理の評価結果を示している。縦軸２６００２は７：３０−７：３１の流出人数を示している。凡例にも示すように、図２６には、各出入口での流出人数の計測結果２６０１１と推定結果２６０１２を比較したグラフが示されている。
【０１３３】
計測結果２６０１１は人数計測部１３００で計測した流出人数を示し、推定結果２６０１２は推定部１２０１で推定した流出人数を示す。図２６には、計測された図７の流出人数情報と、推定された図１７の流出人数情報とが対比されている。
【０１３４】
各出入口についての推定結果の精度を誤差率（流出人数に対する誤差の割合）で評価すると、出入口Ａは誤差率＋１３％であり、出入口Ｂは誤差率が±０％であり、出入口Ｃは誤差率が＋２８％であり、出入口Ｄは誤差率が−１８％である。図２６の説明は以上である。
【０１３５】
評価部１２０３は、この評価結果と、推定部１２０１から出力された推定結果とを、評価結果情報として出力する。推定結果として、移動推定部１２１２に入力された定常ＯＤ情報及びイベントＯＤ情報などをそのまま出力してもよい。
【０１３６】
また、評価部１２０３は、あらかじめ定めておいた繰り返し条件に従い、分岐情報１１０６を更新し、推定部１２０１での推定処理を再度実行するかを判断する。
【０１３７】
前記繰り返し条件の例として、計算時間が一定時間内であること、計算回数が一定回数内であること、各出入口の推定結果の誤差率の絶対値が一定割合以上となるものを含むこと、などの条件のひとつあるいは複数の組み合わせが挙げられる。設定した条件が満たされた場合、評価部１２０３は、分岐情報１１０６を更新し、推定部１０２１に推定処理を再度実行させる。
【０１３８】
例えば、各出入口の評価結果の誤差率の絶対値が２０％以上であるものを含む場合に推定を繰り返すという条件が設定されていた場合、評価部１２０３での評価の結果として図２６の評価結果が得られたら、評価部１２０３は、分岐情報１１０６を更新し、推定部１２０１に推定処理を再度実行させることになる。
【０１３９】
次に、分岐情報１１０６の更新について説明する。
【０１４０】
評価部１２０３は、図２６に示した出入口毎の評価結果の誤差率に応じて、分岐情報１１０６に記録されている降車旅客分岐率情報（図１２）と流入分岐情報（図１３）を更新する。
【０１４１】
図２７は、評価部１２０３における分岐情報更新処理を示すフローチャートである。以下、図２７を用いて分岐情報更新処理について説明する。
【０１４２】
Ｓ２８０１では、評価部１２０３は、評価結果のうち、流出人数の誤差（計測結果と推定結果の差分）が最大となっている出入口ｅおよび期間t1-t2を算出し、またその時の誤差率をｒとする。
【０１４３】
Ｓ２８０２のループで、評価部１２０３は、すべての出入口を走査し処理する。以下、処理対象の出入口をｅ’とする。
【０１４４】
Ｓ２８０３では、評価部１２０３は、出入口ｅ’から流入し、出入口ｅへ流出するとき、出入口ｅでの流出時刻が期間t1-t2となりそうな、出入口ｅ’での流入期間t3-t4を算出する。時刻t3は、時刻ｔ1−（出入口ｅ’から出入口ｅの平均移動時間×(１＋余裕率)）として設定し、時刻t4は、時刻ｔ２−（出入口ｅ’から出入口ｅの平均移動時間×(１＋余裕率)）として設定する。計測期間に対してのズレがある場合、時刻ｔ３は切り捨て、時刻ｔ４は切り上げて対応させる。出入口間の平均移動時間は、空間情報１１０３に記録された出入口間所要時間対応表から取得し、余裕率は予め適当な値を設定しておく。
【０１４５】
Ｓ２８０４では、評価部１２０３は、降車旅客列車分岐率において、出入口ｅ’に接続するホームに期間t3-t4に到着する列車の分岐確率のうち、出入口ｅへ向かう旅客についての値を1/(1+r)倍する。
【０１４６】
Ｓ２８０５では、評価部１２０３は、流入分岐情報に対して、流入出入口がｅ’であり、かつ期間t3-t4の出入口ｅへの分岐率を1/(1+r)倍する。
【０１４７】
Ｓ２８０６では、評価部１２０３は、すべての出入口に対して処理を終えるまでＳ３３０２へ戻って繰返し処理し、処理が完了するとＳ３３０７へ進む。
【０１４８】
Ｓ２８０７では、評価部１２０３は、流入分岐率および降車旅客分岐率の出入口毎の分岐率の比を維持したまま、分岐率の合計が１００％となるように調整し、その後、分岐率更新処理を終了する。分岐率の合計を１００％に維持する処理は、それぞれの分岐率を合計した値で、各分岐率を割ればよい。
【０１４９】
以上で、図２７を用いた分岐情報更新処理の説明を終了する。分岐情報更新処理は、ここに例示した方式に限らない。他の例として、評価結果を用いずに、乱数を用いて、流入分岐率および降車旅客分岐率の出入口毎の分岐率を設定する方式で実現してもよい。
【０１５０】
評価部１２０３の処理についての説明は以上である。また、演算部１２００での処理についての説明も以上である。
【０１５１】
出力部１００１は、評価結果情報１１０７に記録された情報を、本旅客流動推定装置と接続する他装置あるいは他システムに対して、データとして出力する。
【０１５２】
評価結果情報１１０７のデータは各種データを表示する端末に出力してもよい。例えば、評価結果情報１１０７を表形式および／またはグラフ形式で表示する端末に出力し、表示させてもよい。また、評価結果情報１１０７を、移動推定部１２１２の推定処理の課程として、図２５に示したような対象空間内の歩行者の移動の様子を動画で表示する端末（移動状況表示端末）に出力し、表示させてもよい。また、鉄道駅の混雑状況を駅員や旅客に情報提供するため対象空間内の群集密度を表示する端末（群集密度表示端末）に評価結果情報１１０７を出力し、表示させてもよい。
【０１５３】
上述した移動状況表示端末の一実現例として、移動推定部１２１２における推定処理において、逐次更新する前記歩行者エージェントの位置情報を元に動画を作成し、評価部１２０３を経由して、評価結果情報１１０７として出力部１００１から端末に出力し、画面表示させてもよい。
【０１５４】
また、上述した群集密度表示端末の一実現例として、移動推定部１２１２における推定処理において、逐次更新する歩行者エージェントの位置情報を元に、格子空間における単位面積あたりの滞在人数（平均滞在人数）について一定時間の平均値を算出し、評価部１２０３を経て、評価結果情報１１０７として出力部１００１から出力することにしてもよい。端末では、その平均滞在人数に応じて、格子空間の各単位空間を異なる色や塗りつぶしパターンで表示することにより、群集密度の大小を視覚的に画面表示してもよい。
【０１５５】
図２８は、群集密度表示端末の画面の一例を示す図である。以下、図２８を用いて群集密度表示端末の画面について説明する。凡例に示す例７００１、７００２のように、算出した群集密度に応じて格子空間を表示する色や塗りつぶしパターンなどの表現を変化させて、レイアウト情報が示す対象空間中の群集密度の分布を表示させている。単位空間７０１１の群集密度は0.5人/ｍ２以上1.5人/ｍ２未満（例７００１）であり、単位空間７０１２の群集密度は1.5人/ｍ２以上（例７００２）である。この表示によって、対象空間中の旅客によって混雑しやすい場所を確認することができる。
＜作用、効果＞
【０１５６】
本実施例の旅客流動推定装置により、対象空間の出入口での流入人数および流出人数を計測することで、対象空間内の旅客の移動を直接追跡するのではなく、推定することができる。
【０１５７】
これによって、従来把握することが困難であった時間帯毎の乗り換えのためのホーム間の旅客の移動状況として、各ホーム間の移動の割合がＯＤ情報として推定でき、このＯＤ情報を用いることで、現状の駅の各空間の状況を定量的に把握することができ、より旅客にとって利便性のある運行計画や、旅客にとって快適な空間の設計や改修を精度よく評価することができる。
【０１５８】
以上の説明では、計測される旅客の流動人数のうち、計測した流入人数情報を入力のひとつとして流出人数情報を推定し、推定した流出人数情報と、計測した流出人数情報とを比較することにより推定精度を評価したが、本発明がこれに限定されることはない。他の例として、流入と流出を逆にして、計測した流出人数情報を入力のひとつとして流出人数情報を推定し、推定した流入人数情報と、計測した流入人数情報とを比較することにしてもよい。流出人数情報から流入情報を推定する場合、移動推定処理において、繰り返し処理の時間進行を逆にすればよい。
【０１５９】
また、計測した流入人数情報を入力のひとつとして格子空間の群集密度を推定し、また計測した流出人数情報を入力のひとつとして格子空間の群集密度を逆推定し、流入人数情報から推定した格子空間の群集密度と、流出人数情報から推定した格子空間の群集密度を比較することにより、推定精度を評価してもよい。
【０１６０】
以上、本発明の実施例について述べてきたが、本発明は、これらの実施例だけに限定されるものではなく、本発明の技術思想の範囲内において、これらの実施例を組み合わせて使用したり、一部の構成を変更したりしてもよい。
【０１６１】
本実施例では、鉄道駅における旅客の流動を推定するシステムを例示したが、本発明がこれに限定されることはない。本発明は乗物の乗降場所における乗客の流動を推定するシステムに広く適用可能である。他の例として、航空旅客機が発着する空港における旅客の流動を推定するシステムに適用することもできる。
【０１６２】
また、本実施例では、図２７に示したように、降車旅客分岐率と流入分岐率の両方を更新する例を示したが、本発明がこの例に限定されることはない。他の例として、降車旅客分岐率または流入分岐率のいずれか一方のみを更新することにしてもよい。また、降車旅客分岐率と流入分岐率の更新演算にそれぞれ別個の式を用いてもよい。
＜付記＞
【０１６３】
上記の実施例は以下の付記のように整理して表現することも可能である。ただし、本発明が以下の付記に限定されるものでないことは言うまでもない。
（付記１）
【０１６４】
観測の対象となる対象空間の出入口を通過する人の数である流動人数を計測する計測部と、
前記流動人数に基づき、所定の流動推定演算手法により、所定の旅客状態の人の数を推定した比較用推定人数と、前記対象空間での人の流動を推定した流動推定情報とを取得する推定部と、
前記旅客状態の人の数を他の方法で取得した比較対象人数と前記比較用推定人数を比較し、比較結果に基づいて、前記流動推定情報の推定精度を評価する評価部と、
を有する人流動推定システム。
【０１６５】
これによれば、出入口での計測に基づき流動推定演算手法を用いて対象空間での人の流動を推定するとともに、同じ流動推定演算手法で比較用推定人数を推定し、比較用推定人数を他の手法で取得した比較対象人数と比較し、比較結果で流動推定情報の推定精度を評価するので、高い精度で流動推定情報を取得することができる。
（付記２）
【０１６６】
前記計測部は、前記出入口で前記対象空間に流入する人の数である流入人数と、前記出入口で前記対象空間から流出する人の数である流出人数とを計測し、
前記推定部は、前記流入人数と前記流出人数のいずれか一方に基づき他方を推定することにより得られた結果を前記比較用推定人数とし、
前記評価部は、前記比較用推定人数と計測された前記他方とを比較することにより、前記流動推定情報の推定精度を評価する、
付記１に記載の人流動推定システム。
【０１６７】
これによれば、出入口で計測された一方の計測値から流動推定演算手法で他方の計測値に相当する推定値を推定し、その他方の計測値と推定値を比較するので、計測値に基づいて推定精度を正確に評価することができる。
（付記３）
【０１６８】
前記対象空間は、旅客が乗物に乗降する乗降施設内の空間であり、
前記出入口には、乗物の乗降場所に接続される第１出入口と、前記乗降施設の出入口である第２出入口とが含まれ、
前記人流動推定システムは、前記乗物の到着および出発に関する運行情報を保持する記憶部を更に有し、
前記推定部は、前記対象空間に流入する出入口と前記対象空間から流出する出入口との対に対応づけた、前記対象空間に流入して前記対象空間を移動して前記対象空間から流出する人の数を示すＯＤ情報を算出するとき、前記第１出入口からの流入に関する第１ＯＤ情報は前記乗物の到着に対応付けて算出し、前記第２出入口からの流入に関する第２ＯＤ情報は前記時間単位で算出する、
付記２に記載の人流動推定システム。
【０１６９】
これによれば、乗物の到着時に対象空間にまとまって人が流入する特性を有する第１出入口については、乗物の到着に対応付けてＯＤ情報を算出するので、実際の人の流れに、より忠実な推定情報を得ることができる。
（付記４）
【０１７０】
前記記憶部は、到着した乗物から前記第１出入口を通って前記対象空間に流入する人が他の各出入口に向かう第１分岐率の情報と、前記第２出入口から前記対象空間に流入する人が他の各出入口に向かう第２分岐率の情報とを更新可能に更に保持し、
前記推定部は、前記第１分岐率および前記第２分岐率に基づいて、前記第１ＯＤ情報と、前記第２ＯＤ情報と、前記第１ＯＤ情報および前記第２ＯＤ情報に基づき出入口間の移動時間を考慮して算出された、各出入口から流出する人の数である流出推定人数と、を推定し、
前記評価部は、前記流出推定人数を前記比較用推定人数として前記流出人数と比較することにより、前記第１ＯＤ情報および前記第２ＯＤ情報の推定精度を評価し、評価結果に基づいて前記第１分岐率と前記第２分岐率の少なくとも一方を更新するか否か判定する、
付記３に記載の人流動推定システム。
【０１７１】
これによれば、誤差が含まれ易い分岐率を更新しながらＯＤ情報の精度を上げていくことができるので、段階的に実際の状況に近づけるように推定精度を収束させることができる。
（付記５）
【０１７２】
前記評価部は、前記評価結果において前記流出推定人数と前記流出人数の差分が最も大きい出入口の第１分岐率および第２分岐率を更新する、
付記４に記載の人流動推定システム。
【０１７３】
これによれば、誤差が最大の出入口の分岐率を更新して推定精度を上げていくので、短時間でＯＤ情報の精度を向上させることができる。
（付記６）
【０１７４】
前記計測部は、前記出入口で前記対象空間に流入する人の数である流入人数と、前記出入口で前記対象空間から流出する人の数である流出人数とを計測し、
前記推定部は、前記流入人数に基づき前記対象空間を格子状の単位空間に分割した格子空間における人の群集密度である第１群集密度を推定し、前記流出人数に基づき格子空間における人の群集密度である第２群集密度を推定し、
前記評価部は、前記第１群衆密度と前記第２群集密度を比較することにより、前記流動推定情報の推定精度を評価する、
付記１に記載の人流動推定システム。
【０１７５】
これによれば、流入人数と流出人数の両方から求めた群集密度を比較して推定精度を評価するので、流入人数からの推定の精度と、流出人数からの推定の精度を同時に評価することができる。
（付記７）
【０１７６】
前記推定部は、前記対象空間を格子状の単位空間に分割した格子空間を表わす空間情報と前記第１ＯＤ情報および前記第２ＯＤ情報とに基づき前記格子空間の単位空間に人を疑似するエージェントを出現させ、単位空間の間を移動させることにより、前記出入口間の移動時間を算出する、付記４に記載の人流動推定システム。
【０１７７】
これによれば、推定部は、格子空間上のエージェントにより移動時間を算出するので、容易に移動時間を算出し、流出推定人数の算出に用いることができる。
（付記８）
【０１７８】
観測の対象となる対象空間の出入口を通る人を観測する観測システムを更に有し、
前記計測部は、前記観測システムの観測結果に基づいて、前記出入口を通過する人の数である流動人数を計測する、
付記１に記載の人流動推定システム。
（付記９）
【０１７９】
計測手段が、観測の対象となる対象空間の出入口を通過する人の数である流動人数を計測し、
推定手段が、前記流動人数に基づき、所定の流動推定演算手法により、所定の旅客状態の人の数を推定した比較用推定人数と、前記対象空間での人の流動を推定した流動推定情報とを取得し、
評価手段が、前記旅客状態の人の数を他の方法で取得した比較対象人数と前記比較用推定人数を比較し、比較結果に基づいて、前記流動推定情報の推定精度を評価する、
人流動推定方法。
（付記１０）
【０１８０】
観測の対象となる対象空間の出入口を通過する人の数である流動人数を計測する手順と、
前記流動人数に基づき、所定の流動推定演算手法により、所定の旅客状態の人の数を推定した比較用推定人数と、前記対象空間での人の流動を推定した流動推定情報とを取得する手順と、
前記旅客状態の人の数を他の方法で取得した比較対象人数と前記比較用推定人数を比較し、比較結果に基づいて、前記流動推定情報の推定精度を評価する手順と、
をコンピュータに実行させるための人流動推定プログラム。
【符号の説明】
【０１８１】
１０００…旅客流動推定装置、１００１…出力部、１０２１…推定部、１１００…記録部、１１００１…列車ＩＤ、１１００２…降車旅客分岐率、１１０１…流入人数情報、１１０２…流出人数情報、１１０３…空間情報、１１０４…運行情報、１１０６…分岐情報、１１０７…評価結果情報、１２００…演算部、１２０１…推定部、１２０３…評価部、１２１１…ＯＤ推定部、１２１２…移動推定部、１３００…計測部、１３００…人数計測部、１３０１…対象空間、１３１０…カメラシステム、１３１１…カメラ、１３１２…出入口、１３２０…自動改札機システム、１３２１…自動改札機、８０００…格子空間、８０１１…通路、８０１２…壁、８０１３…階段、８０１４…改札、８０１５…出入口
【図１】
2016045919.tif 000003
【図２】
2016045919.tif 000004
【図３】
2016045919.tif 000005
【図４】
2016045919.tif 000006
【図５】
2016045919.tif 000007
【図６】
2016045919.tif 000008
【図７】
2016045919.tif 000009
【図８】
2016045919.tif 000010
【図９】
2016045919.tif 000011
【図１０】
2016045919.tif 000012
【図１１】
2016045919.tif 000013
【図１２】
2016045919.tif 000014
【図１３】
2016045919.tif 000015
【図１４】
2016045919.tif 000016
【図１５】
2016045919.tif 000017
【図１６】
2016045919.tif 000018
【図１７】
2016045919.tif 000019
【図１８】
2016045919.tif 000020
【図１９】
2016045919.tif 000021
【図２０】
2016045919.tif 000022
【図２１】
2016045919.tif 000023
【図２２】
2016045919.tif 000024
【図２３】
2016045919.tif 000025
【図２４】
2016045919.tif 000026
【図２５】
2016045919.tif 000027
【図２６】
2016045919.tif 000028
【図２７】
2016045919.tif 000029
【図２８】
2016045919.tif 000030