JP2004213098A - Congestion prediction system, congestion prediction method and congestion prediction program - Google Patents
Congestion prediction system, congestion prediction method and congestion prediction program Download PDFInfo
- Publication number
- JP2004213098A JP2004213098A JP2002378589A JP2002378589A JP2004213098A JP 2004213098 A JP2004213098 A JP 2004213098A JP 2002378589 A JP2002378589 A JP 2002378589A JP 2002378589 A JP2002378589 A JP 2002378589A JP 2004213098 A JP2004213098 A JP 2004213098A
- Authority
- JP
- Japan
- Prior art keywords
- transmitter
- event
- information
- congestion prediction
- moving object
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
【課題】予測モデルをリアルタイムで更新し、精度良く、将来の混雑予測を行う混雑予測システムを提供する。
【解決手段】混雑予測システム1は、通信ネットワーク25を介して、発信機21を携帯した歩行者22が移動する観測エリア26内の複数の受信機20a、20b、…、20eから送信される発信機情報を受信する受信機制御部13と、発信機情報に基づき、移動イベントを生成する外部イベント取込部12と、移動体の行動を予測するために用いる移動体モデルから移動イベントを生成する模擬イベント発生部18と、外部イベント取込部12と模擬イベント発生部18から送られる移動イベントを切り替えて出力する入力切替部14と、移動体情報を更新することと平行して、将来の混雑予測を行うイベント処理部16とを備える。
【選択図】 図1A congestion prediction system for updating a prediction model in real time and accurately predicting future congestion is provided.
A congestion prediction system (1) transmits via a communication network (25) a plurality of receivers (20a, 20b, ..., 20e) in an observation area (26) in which a pedestrian (22) carrying a transmitter (21) moves. A receiver control unit 13 for receiving device information, an external event capturing unit 12 for generating a mobile event based on transmitter information, and a mobile event generated from a mobile model used to predict the behavior of the mobile. Simulated event generator 18, input switcher 14 for switching and outputting the moving event sent from external event capturing unit 12 and simulated event generator 18, and congestion in the future in parallel with updating the moving body information And an event processing unit 16 for performing prediction.
[Selection diagram] Fig. 1
Description
【0001】
【発明の属する技術分野】
本発明は、ある観測エリア内の移動体の行動を予測する混雑予測システム、混雑予測方法及び混雑予測プログラムに関する。
【0002】
【従来の技術】
従来、イベント会場等の集客施設内に受信機を複数設置し、施設内を移動する歩行者が発信機を所持することにより、歩行者の位置を推定することが行われている。発信機から送信された電波は、その周辺エリアの受信機が検出し、各受信機は、一つの位置検出システムに、無線あるいは有線で発信機から得た情報を送信する。このように、移動通信システムを利用し、歩行者などの移動体の位置の検出が可能な位置検出システム及び位置検出方法は開示されている(例えば、特許文献1参照。)。
【0003】
一方、集客施設において、歩行者の動線の管理は、重要な課題である。ここで、「動線」とは、移動体の移動履歴を指し、場所と時刻の組み合わせで表現できる。集客施設の構造、運用時のレイアウトや案内表示は、歩行者の動線を考慮して設計される。このため、歩行者の動線を予測することが必要となり、予測シミュレーションを使用する場合もある。シミュレーション装置として、人間の集団行動を実際のものに近い形で表現し、分析し、予測する装置は開示されている(例えば、特許文献2参照。)。
【0004】
【特許文献1】
特開平11−178041号公報(第5−6頁、図1)
【0005】
【特許文献2】
特開2000−20499号公報(第4−6頁、図1)
【0006】
【発明が解決しようとする課題】
しかし、従来のシミュレーションで用いられる予測モデルは、集客施設の生データを用いていないため、その信頼性に問題があった。又、新規到来者を考慮しないか、考慮しても単純なポアソン到着モデルを用いていたため、集客施設の時間帯による到来パターンの変化を考慮できず、精度の高い予想ができなかった。
【0007】
上記の問題に鑑み、本発明は、予測モデルをリアルタイムで更新し、精度良く、将来の混雑予測を行う混雑予測システム、混雑予測方法及び混雑予測プログラムを提供することを目的とする。
【0008】
【課題を解決するための手段】
上記目的を達成するため、本発明の第1の特徴は、(イ)通信ネットワークを介して、発信機を携帯した移動体が移動する観測エリア内の複数の受信機から送信される、発信機のIDを含む発信機情報を受信する受信機制御部と、(ロ)発信機情報に基づき、発信機ID、発信機情報を送信した受信機のID、受信機の検出エリアに対する発信機の到着あるいは離脱を示す検出フラグ及び時刻情報を含む移動イベントを生成する外部イベント取込部と、(ハ)移動体の行動を予測するために用いる移動体モデルに基づいて、移動イベントを生成する模擬イベント発生部と、(ニ)外部イベント取込部と模擬イベント発生部から送られる移動イベントを切り替えて出力する入力切替部と、(ホ)入力切替部を介して外部イベント取込部から送られた移動イベントを順次取り出し、移動体の過去の移動履歴で表現される移動体情報を更新し、移動体の現在位置の推定を行うことと平行して、混雑予測を要求された時には、入力切替部を介して模擬イベント発生部から送られた移動イベントを用い、将来の混雑予測を行うイベント処理部とを備える混雑予測システムであることを要旨とする。
【0009】
本発明の第1の特徴に係る混雑予測システムによると、ある観測エリア内を移動する移動体の動線を観測して得られた情報から、移動体モデルをリアルタイムに更新するので、精度良く、将来の混雑予測を行うことができる。
【0010】
又、第1の特徴に係る混雑予測システムは、移動体情報、移動体の検出エリア間の移動状況を示すOD表及び検出エリア毎の検出エリアに滞留した時間を表す滞留時間表を保持する移動体状態保持部と更に備えても良い。
【0011】
更に、第1の特徴に係る混雑予測システムのイベント処理部は、イベント処理部は、混雑予測を要求された時に、非定常ポアソン過程モデルを用いて観測エリアへ新規に到来した移動体を仮発生させ、移動体モデルを作成しても良い。この混雑予測システムによると、観測エリアへの時間帯による到来パターンの変化を考慮し、精度の高い予想をすることができる。
【0012】
本発明の第2の特徴は、通信ネットワークを介して、発信機を携帯した移動体が移動する観測エリア内の複数の受信機から送信される、発信機のIDを含む発信機情報を受信する混雑予測システムにおいて、(イ)発信機情報に基づき、発信機ID、発信機情報を送信した受信機のID、受信機の検出エリアに対する発信機の到着あるいは離脱を示す検出フラグ及び時刻情報を含む移動イベントを順次取り出すステップと、(ロ)移動イベントに基づいて、移動体の過去の移動履歴で表現される移動体情報を更新するステップと、(ハ)移動体情報を用いて、移動体の現在位置の推定を行うステップと、(ニ)混雑予測を要求された時には、混雑予測要求時点での移動体情報を用いて、移動体の行動を予測するために用いる移動体モデルを作成するステップと、(ホ)移動体モデルに基づいて生成された移動イベントを用い、将来の混雑予測を行うステップとを含む混雑予測方法であることを要旨とする。
【0013】
第2の特徴に係る混雑予測方法によると、ある観測エリア内を移動する移動体の動線を観測して得られた情報から、移動体モデルをリアルタイムに更新するので、精度良く、将来の混雑予測を行うことができる。
【0014】
又、第2の特徴に係る混雑予測方法における移動体モデルを作成するステップは、混雑予測要求時点での移動体の検出エリア間の移動状況を示すOD表、及び、混雑予測要求時点での検出エリア毎の検出エリアに滞留した時間を表す滞留時間表を用いて、移動体モデルを作成しても良い。この混雑予測方法によると、実際の動作により近い移動体モデルが作成でき、精度良く、混雑予測を行うことができる。
【0015】
更に、第2の特徴に係る混雑予測方法における移動体モデルを作成するステップは、非定常ポアソン過程モデルを用いて観測エリアへ新規に到来した移動体を仮発生させ、移動体モデルを作成しても良い。この混雑予測方法によると、観測エリアへの時間帯による到来パターンの変化を考慮し、精度の高い予想をすることができる。
【0016】
本発明の第3の特徴は、通信ネットワークを介して、発信機を携帯した移動体が移動する観測エリア内の複数の受信機から送信される、発信機のIDを含む発信機情報を受信する混雑予測システムに、(イ)発信機情報に基づき、発信機ID、発信機情報を送信した受信機のID、受信機の検出エリアに対する発信機の到着あるいは離脱を示す検出フラグ及び時刻情報を含む移動イベントを順次取り出す手順と、(ロ)移動イベントに基づいて、移動体の過去の移動履歴で表現される移動体情報を更新する手順と、(ハ)移動体情報を用いて、移動体の現在位置の推定を行う手順と、(ニ)混雑予測を要求された時には、混雑予測要求時点での移動体情報を用いて、移動体の行動を予測するために用いる移動体モデルを作成する手順と、(ホ)移動体モデルに基づいて生成された移動イベントを用い、将来の混雑予測を行う手順とを実行させる混雑予測プログラムであることを要旨とする。
【0017】
又、第3の特徴に係る混雑予測プログラムにおける移動体モデルを作成する手順は、混雑予測要求時点での移動体の検出エリア間の移動状況を示すOD表、及び、混雑予測要求時点での検出エリア毎の検出エリアに滞留した時間を表す滞留時間表を用いて、移動体モデルを作成しても良い。
【0018】
更に、第3の特徴に係る混雑予測プログラムにおける移動体モデルを作成する手順は、非定常ポアソン過程モデルを用いて観測エリアへ新規に到来した移動体を仮発生させ、移動体モデルを作成しても良い。
【0019】
本発明の第3の特徴に係る混雑予測プログラムを読み出すことにより、混雑予測システム等に上記の手順を実行させることが可能となる。
【0020】
【発明の実施の形態】
次に、図面を参照して、本発明の実施の形態を説明する。以下の図面の記載において、同一又は類似の部分には同一又は類似の符号を付している。但し、図面は模式的なものであることに留意すべきである。
【0021】
本発明の実施の形態に係る混雑予測システムは、ある観測エリア内を移動する歩行者等の移動体の動線を観測して得られた情報から、予測モデルをリアルタイムに更新し、将来の混雑予測を行う。又、予測開始時点ではまだ存在しない歩行者の到来を、非定常ポアソン過程モデルを用いて実際的に仮発生させ、予測に供する。
【0022】
(観測エリア内の受信機及び発信機)
以下において、移動体として、ある観測エリアを移動する歩行者を例にとり、歩行者の動線を観測し、将来の混雑予測を行うシステムについて説明する。本発明の実施の形態に係る混雑予測システム1は、図1に示すように、ある観測エリア26内の歩行者22の動線を観測する。観測エリア26内には、位置を検出したい複数のエリア毎に受信機20a、20b、…、20eを設置する。一方、被観測対象である歩行者22には発信機21を携帯させる。
【0023】
この発信機21は、具体的には、無線タグ等が使用可能である。無線タグは、例えば、いわゆるRF−ID(Radio Frequency Identification)等からなり、図2に示すように、無線タグ30毎に異なるIDを保持する識別情報保持部31と、識別情報保持部31に保持されているIDに応じた所定の周波数の送信信号を生成する送信部32と、送信アンテナ33とを備えている。送信部32は、識別情報保持部31に保持されているIDに応じた送信信号を生成し、送信アンテナ33を介して送信する。
【0024】
この無線タグ30がIDに応じた信号を送信する間隔は、歩行者22の移動速度と、各受信機20a、20b、…、20eの配置間隔に応じて決定し、例えば通常の建物内を歩行によって移動する場合では数秒に1回程度とする。あるいは、会場内を歩行で移動する場合には受信機20a、20b、…、20eの配置間隔を狭くし、IDに応じた信号を送信する間隔を短くすることによって歩行者22の位置を検出する精度を向上させることができる。又、この送信間隔は、送信部32がIDに応じた信号を送信する間隔を制御することによって調整することができる。
【0025】
一方、受信機20a、20b、…、20eも、それぞれ固有のIDを持つ。受信機20a、20b、…、20eは、発信機21から発信機固有のIDを含む電波を検出する。受信機20a、20b、…、20eが電波を検出する領域は調整可能であり、調整によって検出エリアであるエリアA、エリアB、…、エリアEの範囲が変化する。受信機20a、20b、…、20eは、発信機21の検出に関する情報を、通信ネットワーク25を介して混雑予測システム1内の受信機制御部13に送出する。
【0026】
次に、受信機20a、20b、…、20eが発信機21から受信する情報の処理方法について説明する。受信機20a、20b、…、20eは、発信機に関する情報(以下、「発信機情報」という。)のリストを保持する。発信機情報には、図3に示すように、発信機のID(タグID)、現在検出中か否かを示すフラグ、検出フラグが変化した時刻からの経過時間を示すカウンタが含まれる。
【0027】
受信機20a、20b、…、20eは、発信機21の電波を検出すると、発信機IDをキーに発信機情報リストを検索し、見つかれば発信機情報の更新処理を行い、見つからなければそのIDに関する発信機情報を新規に作成する。例えば、発信機ID「XXX」の発信機21の電波を検出し、その発信機情報がリストに無かった場合は、図3に示すように、発信機ID「XXX」の検出フラグをtrue、カウンタを0として、発信機ID「XXX」のリストを新規に作成する。以降、発信機ID「XXX」の発信機情報が検出されるたびカウンタを0にセットするとともに、検出フラグがfalseの場合には検出フラグをtrueに変更する。受信機20a、20b、…、20eは、所定の時間間隔tで図3のテーブルを調査し、すべての発信機についてカウンタを1加算する。又、受信機20a、20b、…、20eは、図3のテーブルにおいて、検出フラグがtrueの発信機についてそのカウンタを調べ、そのカウンタ値が所定の値cthを超えていれば、検出フラグをfalseに変更し、カウンタを0とする。例えば、図3の発信機ID「YYY」においては、検出フラグがfalseに変更された直後の状態を示し、一方発信機ID「ZZZ」においては検出フラグがfalseに変更されてから15ステップ時間が経過していることを示している。即ち、検出フラグがtrueの場合には、発信機IDに対応する発信機21が該当エリアで検出されていることが分かる。又、検出フラグがfalseの場合には、発信機IDに対応する発信機21が該当エリアにカウンタ値が示す時間検出されていないことが分かる。
【0028】
次に、受信機20a、20b、…、20eが、ネットワーク25を介して混雑予測システム1内の受信機制御部13に情報を送出するタイミングについて、説明する。受信機20a、20b、…、20eは次の2種類のケースにおいて、発信機情報、発信機情報を送信した受信機のIDを「受信機情報」として、受信機制御部13に送出する。
【0029】
1.検出フラグがtrueに変化した時点
2.検出フラグがfalseに変化した時点
検出フラグがtrueに変化した時点とは、あるエリアへ歩行者22が到着した時点を指す。又、検出フラグがfalseに変化した時点とは、あるエリアから歩行者22が離脱した時点を指す。検出フラグがfalseに変化する時点は、trueの状態におけるカウンタ値が所定の値cthを経過した時点であるので、検出エリアの変動や電波の混線ではなく、確実に歩行者22がエリアを離脱したとみなせる時点である。この処理によって、混雑予測システム1には、あるエリアへの歩行者の到着と、同エリアからの離脱が伝えられることになる。全エリアからの情報を総合すれば、観測エリア26内の歩行者22の移動履歴を得ることができる。
【0030】
(混雑予測システムの構成)
本発明の実施の形態に係る混雑予測システム1は、図1に示すように、通信ネットワーク25を介して、発信機21を携帯した歩行者22が移動する観測エリア26内の複数の受信機20a、20b、…、20eから送信される発信機情報を受信する受信機制御部13と、発信機情報に基づき、移動イベントを生成する外部イベント取込部12と、移動体の行動を予測するために用いる移動体モデルに基づいて、移動イベントを生成する模擬イベント発生部18と、外部イベント取込部12と模擬イベント発生部18から送られる移動イベントを切り替えて出力する入力切替部14と、移動体情報を更新し、歩行者の現在位置の推定を行うことと平行して、将来の混雑予測を行うイベント処理部16とを備える。更に、本発明の実施の形態に係る混雑予測システム1は、記録イベント取込部11、移動体状態保持部15、環境情報保持部17とを備える。
【0031】
受信機制御部13は、通信ネットワーク25を介して観測エリア26内にある受信機20a、20b、…、20eの動作制御を行うとともに、受信機20a、20b、…、20eから送られる受信機情報を外部イベント取込部12へ送出する。発信機情報に発信機情報を送信した受信機のIDを付加したものを「受信機情報」という。
【0032】
外部イベント取込部12は、受信機制御部13から送られてくる受信機情報に時刻情報を付加して入力切替部14に配送する。ここで、受信機情報に時刻情報を付加したものを、「移動イベント」と呼ぶ。
【0033】
記録イベント取込部11は、記録ファイル5に記述された移動イベントを読み取り、入力切替部14に送出する。記録ファイル5は、過去の移動イベントを記録したファイルであり、FDなどの記録媒体に保存されていても良い。外部イベント取込部12が観測エリア26内の現在の情報を取り込むことに対し、記録イベント取込部11は、観測エリア26内の過去の情報を取り込む。
【0034】
模擬イベント発生部18は、移動体モデル8に基づいて移動イベントを生成し、入力切替部14に送出する。移動体モデル8は、観測エリア26内の歩行者(移動体)の行動を予測するために用いる予測モデルである。移動体モデル8の構築方法は、後に詳述する。
【0035】
入力切替部14は、外部イベント取込部12、記録イベント取込部11、模擬イベント発生部18から送られる移動イベントを、モードに応じて切り替え、イベント処理部16へ送信する。混雑予測システム1は、以下の3種類の動作モードを備える。
【0036】
1.モニタモード:観測エリア内の現在の様子を観測する。外部イベント取込部12から送られた移動イベントを用いる。
【0037】
2.再現モード:過去に記録した観測エリア内の情報から当時の様子を観測する。記録イベント取込部11から送られた移動イベントを用いる。
【0038】
3.シミュレーションモード:予測モデルを用いて、仮想環境における仮想歩行者の様子を予測する。模擬イベント発生部18から送られた移動イベントを用いる。但し、シミュレーションモードであっても、平行して、現在の観測エリア内の観察は継続される。
【0039】
モードの切替えは、入力部9等のユーザインタフェースを介して行われる。入力部9は、具体的には、キーボード、マウス等の機器を指す。入力部9から入力操作が行われると対応するキー情報が入力切替部14に伝達される。
【0040】
移動体状態保持部15は、全歩行者22の現在の状態と動線を保持する。具体的には、図4に示す移動体情報や移動体の現在位置を示す状態バッファを備える。例えば、図5に示すように、受信機20aによる検出エリアA、受信機20bによる検出エリアB、受信機20cによる検出エリアC、受信機20dによる検出エリアD、受信機20eによる検出エリアE、受信機20fによる検出エリアFを有する観測エリア26を歩行者が移動したとする。このとき、移動体状態保持部15は、各受信機20a、20b、…、20fから得られた情報から、各発信機ID毎の移動体情報を保持する。図4では、対象となる発信機21は、エリアAに15:30に到着し、30分間留まり、その後、エリアBに16:05に到着し、10分留まり、次に、エリアC、エリアE、エリアFに移動したことを示している。
【0041】
更に、移動体状態保持部15は、歩行者のエリア間の移動状況を示すOD表やエリア毎のエリアに滞留した時間を表す滞留時間表等を保持する。OD表は、例えば、滞留エリアの数をNとしたとき、整数を要素とするN×N次元マトリクスとなる。i及びjを1〜Nまでの整数とすると、i行j列要素は滞留エリアiから滞留エリアjへ移動した人の累積数とする。歩行者の移動イベントが起こるたびに対応する要素を1ずつインクリメントする。滞留時間表は、例えば、3つのN次元ベクトル、あるいは3×N次元行列で構成する。3つのN次元ベクトルを母数ベクトル、累積ベクトル、二乗累積ベクトルとし、滞留エリアiにおける新たな滞留時間をxiとしたとき、母数ベクトルの第i要素ni、累積ベクトルの第i要素ai、二乗累積ベクトルの第i要素siは以下のように表される。
【0042】
ni = ni + 1
ai = ai + ln(xi + 1)
si = si + {ln(xi + 1)}2
記録される滞留は、OD表と同一の条件を満たすイベントのみとする。
【0043】
移動体状態保持部15は、RAM等の内部記憶装置を用いても良く、HDやFD等の外部記憶装置を用いても良い。
【0044】
環境情報保持部17は、観測エリア26に関する環境情報を保持する。例えば、検出エリア間の位置関係や距離、通路情報を保持する。図5に示すような観測エリア26であると、エリアAとエリアCの間には仕切り27があるという情報を保持する。又、観測エリア26に推奨移動順路が提示されている場合は、その順路を保持する。更に、多くの発信機の動線から、全体として通ることが多い順路を保持する。これらの環境情報を用いて、イベント処理部16は、移動体状態保持部15の情報を更新する。環境情報保持部17は、移動体状態保持部15と同様に、RAM等の内部記憶装置を用いても良く、HDやFD等の外部記憶装置を用いても良い。
【0045】
イベント処理部16は、移動イベントから移動体状態保持部15に保持された情報を更新し、歩行者の位置の推定を行う。この更新処理には、環境情報保持部17に保持された環境情報を用いる。又、イベント処理部16は、図8に示すように、イベントキュー(待ち行列)とハンドラを備える。入力切替部14から供給される移動イベントは、イベントキューの最後に順次格納される。移動イベントを処理する際、ハンドラは、イベントキューの先頭から移動イベントを取り出し、順次処理を行う。そして、移動イベントを処理する毎に、移動体状態保持部15に保持される移動体情報、OD表、滞留時間表等を作成、更新する。
【0046】
イベント処理部16によって得られた歩行者ごとの動線は、動線記録ファイル6あるいは表示部7へ出力される。ここで、動線記録ファイル6とは、歩行者の動線、即ち、場所と時刻の組み合わせで表現できる歩行者の移動履歴を記録したファイルである。又、表示部7とは、モニタなどの画面を指し、液晶表示装置(LCD)、発光ダイオード(LED)パネル、エレクトロルミネッセンス(EL)パネル等が使用可能である。又、表示部7はプリンターなどでも構わない。
【0047】
更に、イベント処理部16は、混雑予測を要求された時には、要求時点での移動体状態保持部15に保持される移動体情報、OD表、滞留時間表等を初期値として用い、予測モデルとなる移動体モデル8をリアルタイムで更新する。そして、この更新された移動体モデル8によって生成された移動イベントを用い、将来の混雑予測を行う。又、移動体モデル8を生成する際は、非定常ポアソン過程モデルを用いて新規到来者を仮発生させるので、観測エリアへの時間帯による到来パターンの変化を考慮し、精度の高い予想をすることができる。
【0048】
尚、図示してはいないが、本発明の実施の形態に係る混雑予測システムは、予測モデルの更新や混雑予測を実行させるためのプログラムを保存するプログラム保持部を備えていても良い。プログラム保持部は、RAM等の内部記憶装置を用いても良く、HDやFD等の外部記憶装置を用いても良い。
【0049】
本発明の実施の形態に係る混雑予測システムによると、ある観測エリア内を移動する歩行者の動線を観測して得られた情報から、移動体モデルをリアルタイムに更新するので、精度良く、将来の混雑予測を行うことができる。
【0050】
(動線観測方法)
次に、動線観測方法について、図6を用いて説明する。イベント処理部16は、入力切替部14から受け取る移動イベントを順次処理し、発信機21毎の動線を生成する。
【0051】
(イ)まず、ステップS101において、イベント処理部16は、イベントキューの先頭から移動イベントをひとつ取り出す。
【0052】
(ロ)次に、ステップS102において、移動体状態保持部15に保持された移動体情報及び状態バッファから、移動イベントの発信機IDに該当する移動体情報及び状態バッファを検索する。
【0053】
(ハ)次に、ステップS103において、移動体情報及び状態バッファを検索できなかった場合には、ステップS104に進み、新規に発信機ID用の移動体情報及び状態バッファを生成し、新規登録する。
【0054】
(ニ)次に、ステップS105において、移動イベントに基づいて、環境情報保持部17に保持された環境情報を考慮し、状態バッファを更新する。例えば、エリア間の位置関係や距離、通路情報を利用して、2エリア間の可能経路を絞り込んだ上で、判定処理を行えば、現実にはあり得ない移動を推定することはなくなる。又、推奨移動順路の提示される環境においては、順路とは逆方向に移動する可能性は順方向よりも低いと考えられるため、推定時点までの動線と照合することにより、より確からしい存在位置推定が可能である。更に、他の発信機の動線から全体として通ることが多い順路と、推定対象の発信機の動線を照合することにより、より確からしい存在位置推定が可能である。“より確からしい”推定と行うというのは、どちらのエリアが正しいのか判断ができない場合に、確率的にひとつのエリアを選択することで、その際に利用する確率を順路や動線で決めるということである。更に、移動体状態保持部15に保持されるOD表や滞留時間表の作成、更新を行う。その他、検出エリアにCCDカメラ等の他のデバイスを設置し、この検出結果を考慮して、更新処理を行っても構わない。
【0055】
(ホ)次に、ステップS106において、更新された状態バッファにもとづいて、発信機の位置、即ち、歩行者の存在エリアを推定する。そして、推定された存在エリアを移動体情報のエリア部分等に格納する。
【0056】
(へ)次に、ステップS107において、ステップS105及びS106で状態バッファあるいは存在エリアが変化した場合には、ステップS108に進み、動線記録ファイル6を出力する。あるいは、表示部7に動線情報を表示する。変化が無い場合は、ステップS101に戻り、上記の処理を繰り返す。
【0057】
本発明の実施の形態に係る動線観測方法によると、歩行者の位置を精度良く検出し、歩行者の動線を途切れることなく観測できる。
【0058】
(混雑予測方法)
次に、混雑予測方法について、図7及び図8を用いて説明する。混雑の予測要求が行われると、本発明の実施の形態に係る混雑予測システム1は、シミュレーションモードに切り替わる。イベント処理部16は、予測要求時のイベントキュー及びハンドラ、移動体状態保持部15の移動体情報、OD表、滞留時間表等をコピーし、予測用シミュレーションインスタンスを新たに生成する。動線観測を行っている状態は、そのままモニタ用インスタンスとして存在し、引き続き動線観測を継続する。即ち、イベント処理部16及び移動体状態保持部15は、仮想的に2つのインスタンスを有することとなる。ここで、「インスタンス」とは、ある処理を行う一つの状態を指す。
【0059】
(イ)まず、図7のステップS201において、将来の混雑予測を行うための要求がイベント処理部16に伝達される。この予測要求は、図1に示す入力部9等を用いて、入力切替部14を介してイベント処理部16に伝達されても良い。このとき、現在からどのくらいの時間が経過した後の状態を予測するか、入力部9等を通じてイベント処理部16に伝達する。予測要求が行われると、イベント処理部16は、図8に示すように、現在の状態を元に、予測用シミュレーションインスタンスを新たに生成する。シミュレーションモードにおいては、入力切替部14は、外部イベント取込部12からの移動イベントをモニタ用インスタンスへ供給し、模擬イベント発生部18からの移動イベントを予測用シミュレーションインスタンスへ供給する。
【0060】
(ロ)ステップS202において、モニタ用インスタンスは、動線観測処理を継続する。即ち、図6のフローチャートに示す手順で、図8に示すように、入力切替部14から送られた移動イベントをイベントキューに格納し、ハンドラによって、移動イベントの処理を行う。そして、移動イベント毎に移動体状態保持部15の移動体情報、OD表、滞留時間表等の更新を行う。
【0061】
(ハ)一方、ステップS203において、予測用シミュレーションインスタンスは、予測要求時の移動体情報、OD表、滞留時間表をコピーし、これらを初期値として、移動体モデル8のパラメータを計算する。
【0062】
(ニ)次に、ステップS204において、移動体モデル8を構築する。このとき、予測開始時点でまだ存在していない歩行者の到来を、非定常ポアソン過程モデルを用いて実際的に仮発生させる。非定常ポアソン過程モデルを用いることにより、観測エリアへの時間帯による到来パターンの変化を考慮し、精度の高い予想をすることができる。移動体モデル8の構築方法は、後に詳述する。予測用シミュレーションインスタンスは、構築した移動体モデル8を模擬イベント発生部18へ送る。模擬イベント発生部18は、移動体モデル8に基づいて移動イベントを生成し、入力切替部14に送出する。入力切替部14は、この移動イベントを予測用シミュレーションインスタンスへ供給する。
【0063】
(ホ)次に、ステップS205において、予測用シミュレーションインスタンスは、移動イベントに基づき、指定時間経過後の予測シミュレーションを行う。この予測シミュレーションについては、後に詳述する。
【0064】
本発明の実施の形態に係る混雑予測方法によると、ある観測エリア内を移動する歩行者の動線を観測して得られた情報から、移動体モデルをリアルタイムに更新するので、精度良く、将来の混雑予測を行うことができる。
【0065】
(移動体モデルの構築方法)
ここで、移動体情報、OD表、滞留時間表等を元に、移動体モデルを構築する方法について説明する。
【0066】
本発明の実施の形態に係る混雑予測システムで用いる移動体モデルは、以下のサブモデルから構成される。
【0067】
・来場モデル…歩行者の到着時刻を決定
・移動モデル…見学する会場を選択
・滞留モデル…各会場での見学所要時間を決定
移動体モデルの模式図を図9に示す。来場モデルに従って、歩行者は来場し、移動モデルに従って、退場時刻になるまで会場を選択・移動し、滞留モデルに従って、見学する各会場で滞留する。以下にそれぞれのモデルの構築例を説明する。ここでは、図10に示すような展示会場を見学する際の歩行者の動線を予測する。
【0068】
<来場モデル>
見学者の来場密度を一定とみなすのは不適切なので、定式化には非定常ポアソン過程(NHPP: Non-homogeneous Poisson Process)を導入する。NHPPの詳細は、「D. P. Heyman, M. J. Sobel, Eds.: “Stochastic Models”, Handbooks in Operations Research and Management Science, Vol. 2, Elsevier Science Publishers B.V. (North-Holland), 1990.」に記載されている。NHPPによると、時間によって変化する来場密度をもつポアソン到着モデルが作成できる。時刻Tにおける次の来場確率分布は、λを来場密度として、次式で与えられる。
【0069】
【数1】
時刻Tから次の来場時刻までの時間は、本式をtについて解けばよい。式(1)を変形し、一様乱数uを導入して次式を得る。
【0070】
【数2】
来場は午前と午後で分け、それぞれの来場時刻分布を正規分布で近似する。これによると、来場密度λ(t) は、μを平均来場時刻、σを標準偏差、AMを午前、PMを午後として、次式で表される。
【0071】
【数3】
来場モデルの同定としては、実際の午前、午後の来場時刻分布に正規分布をフィッティングさせることにより、例えば、次のような来場モデルが得られる。平均値は9:00を0秒としたときの経過時間を示す。
【0072】
μAM=3855, σAM=2703
μPM=18838, σPM=3573 ………(4)
来場者の来場時刻を推定するには、午前来場総数NAMおよび午後来場総数NPMが必要だが、これについては別途与えられるものとする。
【0073】
<移動モデル>
来場者の移動モデルは、次の2つのモデルを合成して定式化する。
【0074】
・退場モデル
・次会場選択モデル
即ち、来場者がある会場での見学を終了したとき、まず退場モデルを用いて退場するか否かを決定する。もし退場しない結論が出た場合、次会場選択モデルを用いて次に訪問する会場を決定する。
【0075】
退場モデルは、退場時刻をtL、来場時刻をtA、滞在時間をtS、閉場時効をTCとして、次式で定義する。
【0076】
【数4】
DL(t)がtrueの場合、退場を選択する。falseの場合、次会場選択モデルにより次の移動先会場を決定する。滞在時間は次式で定義される確率分布(正規分布)により決める。
【0077】
【数5】
平均値μSと標準偏差σSは来場者の来場時刻tAの一次関数として与える。これは、午前の来場者は午前中に、午後の来場者は閉場時刻までに、それぞれ見学を完了しようとする傾向を表現している。
【0078】
【数6】
退場時刻tLは来場者の状態に依存しない。
【0079】
次会場選択モデルは、マルコフ連鎖モデルを用いて定式化する。従って、過去の見学履歴は考慮しないモデルとなる。マルコフ連鎖モデルの詳細は、「金子秀二著、“展示会における人の流れのシミュレーション”, アドバンス, No. 6, pp. 33-37, 東芝アドバンストシステム, 1998.」に記載されている。
【0080】
各来場者の会場間移動を表すマルコフ連鎖をXkとする。来場のための受付をひとつの会場とみなし、受付を含めた展示会場数がMだとする。受付を会場1、その他の会場を会場2〜会場Mと呼ぶ。来場者の状態数はMである。いま、会場iに滞在する状態をbiで表し、時点kおいて会場iに滞留する確率分布を次式のように書く。
【0081】
【数7】
これはM次元ベクトルとして、以下のように書ける。
【0082】
p[k] = [p1[k], p2[k], …,pM[k]] ………(9)
全来場者の初期分布は必ずp[0]=[1,0,…,0]である。推移確率行列をM×M次元行列Pで表す。推移確率行列Pのi行j列要素は、会場iから会場jへの推移確率である。
【0083】
pij = P(Xk+1 = bj | Xk = bi) ………(10)
各来場者の時点kでの確率分布は次式で表される。
【0084】
p[k] = p[0]Pk ………(11)
本モデルでは推移確率行列Pを、対象となる展示会開催中、不変とみなす。
【0085】
退場モデルの同定としては、閉場時刻TCを17:00とし、会場時刻9:00からの秒数で表すと次式となる。
【0086】
TC = 28800 ………(12)
滞在時間tSの分布については、以下のように求める。
【0087】
1.サンプルをそれぞれの来場時刻によって、午前と午後の2グループに分割する。
【0088】
2.午前のグループから全日見学者を除く。
【0089】
3.両グループについて、一時間ごとに滞在時間の平均値と標準偏差を計算する。
【0090】
4.両グループについて、平均値と標準偏差の時間変化を一次関数で同定する。この際、午後グループには以下のサンプルを追加しておく。
【0091】
μS(TC) = 0, σS(TC) = 0 ………(13)
標準偏差の計算は実際の平均による固定値とする。この結果、例えば、滞在時間tSの確率分布は以下のパラメータ(単位は秒)で表す正規分布として得られる。
【0092】
αAM = -0.3521, αPM = -0.6696
βAM = 9752.4, βPM = 19904.4
ξAM = 0, ξPM = -0.2233
ζAM = 3704.4, ξPM = 6508.8 ………(14)
次会場選択モデルの同定としては、推移確率行列P は、実測したOD 表から、次の手順で導出する。
【0093】
1.S列を削除し13×12 次元行列とする。
【0094】
2.各行について正規化する。
【0095】
3.行列を転置する。
【0096】
4.先頭行に1 行挿入し、13×13 次元正方行列とする。挿入する行ベクトルは[0,…,0]である。
【0097】
処理結果を図11に示す。
【0098】
<滞留モデル>
各会場における滞在時間は、ある時不変な確率分布にしたがうものとして定式化する。滞在時間の対数変換結果が正規分布にしたがうものと考える。即ち、会場iにおける滞在時間TSiは次式で与えられる。
【0099】
【数8】
xiは次式で定義される確率分布により決める。
【0100】
【数9】
各会場の滞在時間記録を対数変換し、それぞれについて正規分布にフィッティングした結果の一例を図12に示す。又、平均値を逆変換すると図13を得る。
【0101】
上記のように作成された来場モデル、移動モデル、滞留モデルを参照し、移動体モデルを構築する。
【0102】
(予測シミュレーション)
上記で得られた移動体モデルに基づき、移動イベントを生成し、会場の混雑予測は行われる。以下に、シミュレーションによる会場混雑の再現を行うアルゴリズムを示す。
【0103】
1. 午前と午後の来場者数を、それぞれN_AM,N_PM として設定する。
【0104】
2. 時刻の初期化t=0(午前9:00);来場者カウンタの初期化
i=0
3. while (i < N_AM) {
4. i=i+1
5. 来場モデルにより来場者i 発生;時刻t_i=t+発生間隔
6. 退場モデルにより来場者i の退場時刻t_L[i]決定
7. while (t_i < t_L[i]) {
8. 次会場選択モデルにより訪問会場決定,移動
9. 滞留モデルにより滞留時間t_S 決定;t_i=t_i+t_S
10. }
11. 退場処理
12. }
13. 2〜12 の処理を午後について同様に実施;t=10800(正午)
14. 全来場者の行動データを時刻順に並べ替え
15. 集計処理
シミュレーションプログラムはJava(登録商標)等で作成可能である。
【0105】
シミュレーションにより得られた歩行者情報、OD表、滞留時間表等は、移動体状態保持部15に保存されても良く、シミュレーション結果保持部(図示せず)等に保持されても良い。あるいは、表示部7等に表示されても良い。これらのシミュレーション結果を用いて得られた来場時刻分布の一例を図14に示す。図14では、NAM=354、NPM=660を実際の来場者数として用いている。又、同様に、滞留時間表等のシミュレーション結果を用いて得られた見学時間分布を図15に示す。又、時間毎の各エリアの混雑推移(滞在者数推移)を図16に示す。R、Q、…、Aは、図10の展示会場エリアに相当する。更に、シミュレーション結果から得られたOD表を図17に示す。各数値は、出発会場から到着会場に移動した歩行者の数を示す。これらのシミュレーション結果から、会場の混雑を予測し、集客施設の構造、運用時のレイアウトや案内表示の適切な設計が可能となる。
【0106】
(その他の実施の形態)
本発明は上記の実施の形態によって記載したが、この開示の一部をなす論述及び図面はこの発明を限定するものであると理解すべきではない。この開示から当業者には様々な代替実施の形態、実施例及び運用技術が明らかとなろう。
【0107】
例えば、本発明の実施の形態において、観測エリアを移動する移動体は、歩行者であると説明したが、これに限らないことは勿論である。
【0108】
又、本発明の実施の形態において、外部イベント取込部12で時刻情報を付加すると説明したが、受信機20a、20b、…、20eが受信機制御部13へ情報を送出する際に、時刻情報を付加しても構わない。但し、受信機20a、20b、…、20eよりも、外部イベント取込部12で時刻情報を付加するほうが、受信機20a、20b、…、20e間の時刻の誤差を考慮する必要がないという利点がある。
【0109】
又、本発明の実施の形態において、移動体状態保持部15、環境情報保持部17を分けて備えると説明したが、一つの保持部で共用しても構わない。又、本発明の実施の形態に係る混雑予測システムは、状態バッファやイベントキューを保持する状態バッファ保持部やイベントキュー保持部を別個に備えても構わない。このとき、状態バッファ保持部やイベントキュー保持部は、移動体状態保持部15等と同様に、RAM等の内部記憶装置を用いても良く、HDやFD等の外部記憶装置を用いても良い。
【0110】
このように、本発明はここでは記載していない様々な実施の形態等を含むことは勿論である。従って、本発明の技術的範囲は上記の説明から妥当な特許請求の範囲に係る発明特定事項によってのみ定められるものである。
【0111】
【発明の効果】
本発明によると、予測モデルをリアルタイムで更新し、精度良く、将来の混雑予測を行う混雑予測システム、混雑予測方法及び混雑予測プログラムを提供することができる。
【図面の簡単な説明】
【図1】本発明の実施の形態に係る混雑予測システムの構成ブロック図である。
【図2】本発明の実施に形態に係る発信機に備えられた無線タグの構成ブロック図である。
【図3】本発明の実施の形態に係る発信機情報の一例である。
【図4】本発明の実施の形態に係る移動体情報の一例である。
【図5】本発明の実施の形態に係る観測エリアの一例である。
【図6】本発明の実施の形態に係る動線観測方法のフローチャートである。
【図7】本発明の実施の形態に係る混雑予測方法のフローチャートである。
【図8】本発明の実施の形態に係る混雑予測方法の2つのシミュレーションインスタンスを示す図である。
【図9】本発明の実施の形態に係る移動体モデルの概念図である。
【図10】展示会場のレイアウトの一例である。
【図11】本発明の実施の形態に係る移動体モデルを構築する際に用いる推移確率行列である。
【図12】本発明の実施の形態に係る滞留モデルの同定に用いる滞在時間分布である。
【図13】本発明の実施の形態に係る滞留モデルの滞留滞在時間平均値を示す図である。
【図14】本発明の実施の形態に係る混雑予測方法によって得られる来場時刻分布の一例である。
【図15】本発明の実施の形態に係る混雑予測方法によって得られる滞在時刻分布の一例である。
【図16】本発明の実施の形態に係る混雑予測方法によって得られる混雑推移の一例である。
【図17】本発明の実施の形態に係る混雑予測方法によって得られるOD表の一例である。
【符号の説明】
1 混雑予測システム
5 記録ファイル
6 動線記録ファイル
7 表示部
8 移動体モデル
9 入力部
11 記録イベント取込部
12 外部イベント取込部
13 受信機制御部
14 入力切替部
15 移動体状態保持部
16 イベント処理部
17 環境情報保持部
18 模擬イベント発生部
20a、20b、…、20e 受信機
21 発信機
22 歩行者
25 通信ネットワーク
26 観測エリア
27 仕切り
30 無線タグ
31 識別情報保持部
32 送信部
33 アンテナ[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to a congestion prediction system, a congestion prediction method, and a congestion prediction program for predicting the behavior of a moving object in a certain observation area.
[0002]
[Prior art]
2. Description of the Related Art Conventionally, a plurality of receivers are installed in a facility for attracting customers, such as an event venue, and a pedestrian moving in the facility has a transmitter to estimate the position of the pedestrian. Radio waves transmitted from the transmitter are detected by receivers in the surrounding area, and each receiver transmits information obtained from the transmitter wirelessly or by wire to one position detection system. As described above, a position detection system and a position detection method capable of detecting the position of a moving object such as a pedestrian using a mobile communication system are disclosed (for example, see Patent Document 1).
[0003]
On the other hand, in a customer attracting facility, management of pedestrian traffic lines is an important issue. Here, the “flow line” indicates a movement history of the moving object, and can be expressed by a combination of a place and a time. The structure of the attracting facilities, the layout during operation, and the guidance display are designed in consideration of pedestrian traffic lines. Therefore, it is necessary to predict a pedestrian's flow line, and a prediction simulation may be used. As a simulation device, a device that expresses, analyzes, and predicts human collective behavior in a form close to an actual one has been disclosed (for example, see Patent Document 2).
[0004]
[Patent Document 1]
JP-A-11-178041 (pages 5-6, FIG. 1)
[0005]
[Patent Document 2]
JP-A-2000-20499 (page 4-6, FIG. 1)
[0006]
[Problems to be solved by the invention]
However, since the prediction model used in the conventional simulation does not use the raw data of the customer attracting facility, there is a problem in its reliability. In addition, since a new arrival is not taken into account or a simple Poisson arrival model is used even if it is taken into account, it is not possible to take into account changes in arrival patterns depending on the time zone of the customer attraction facility, and it is not possible to make highly accurate predictions.
[0007]
In view of the above problems, an object of the present invention is to provide a congestion prediction system, a congestion prediction method, and a congestion prediction program for updating a prediction model in real time and accurately predicting future congestion.
[0008]
[Means for Solving the Problems]
In order to achieve the above object, a first feature of the present invention is that (a) a transmitter, which is transmitted from a plurality of receivers in an observation area where a mobile body carrying the transmitter moves via a communication network; A receiver control unit that receives the transmitter information including the ID of the transmitter, and (b) the transmitter ID, the ID of the receiver that transmitted the transmitter information, and the arrival of the transmitter in the detection area of the receiver based on the transmitter information. Alternatively, an external event capturing unit that generates a moving event including a detection flag indicating time of departure and time information, and (c) a simulated event that generates a moving event based on a moving object model used to predict the behavior of the moving object A generating unit, (d) an input switching unit for switching and outputting a moving event sent from the external event capturing unit and the simulated event generating unit, and (e) transmitting from the external event capturing unit via the input switching unit. In parallel with extracting moving events, updating the moving object information expressed in the past moving history of the moving object, and estimating the current position of the moving object, input switching is performed when congestion prediction is requested. The present invention provides a congestion prediction system including an event processing unit that performs future congestion prediction using a movement event transmitted from a simulation event generation unit via a unit.
[0009]
According to the congestion prediction system according to the first aspect of the present invention, a moving body model is updated in real time from information obtained by observing a flow line of a moving body moving in a certain observation area. Predict future congestion.
[0010]
In addition, the congestion prediction system according to the first feature is a movement that holds moving object information, an OD table indicating a moving state of the moving object between detection areas, and a stay time table indicating a time of staying in the detection area for each detection area. It may further include a body condition holding unit.
[0011]
Further, the event processing unit of the congestion prediction system according to the first feature is configured such that, when the congestion prediction is requested, the event processing unit temporarily generates a moving object newly arriving at the observation area using the unsteady Poisson process model. And a moving body model may be created. According to the congestion prediction system, highly accurate prediction can be performed in consideration of a change in an arrival pattern depending on a time zone to an observation area.
[0012]
A second feature of the present invention is to receive, via a communication network, transmitter information including a transmitter ID, which is transmitted from a plurality of receivers in an observation area where a mobile body carrying the transmitter moves. In the congestion prediction system, (a) based on the transmitter information, includes a transmitter ID, an ID of a receiver that transmitted the transmitter information, a detection flag indicating arrival or departure of the transmitter from a detection area of the receiver, and time information. (B) updating the moving body information represented by the past moving history of the moving body based on the moving event; and (c) using the moving body information, E) estimating the current position; and (d) creating a mobile model used to predict the behavior of the mobile using the mobile information at the time of the congestion prediction request when congestion prediction is requested. Comprising the steps of, and summarized in that a (e) using the moving events generated based on the mobile body model, congestion prediction method and performing a future congestion prediction.
[0013]
According to the congestion prediction method according to the second feature, the mobile object model is updated in real time from the information obtained by observing the flow line of the mobile object moving in a certain observation area, so that the congestion of the future can be accurately and accurately determined. You can make predictions.
[0014]
The step of creating a moving object model in the congestion prediction method according to the second feature includes the steps of: an OD table indicating a movement state of the moving object between detection areas at the time of the congestion prediction request; A moving object model may be created using a stay time table indicating the time spent in the detection area for each area. According to this congestion prediction method, a mobile body model closer to the actual operation can be created, and congestion prediction can be performed with high accuracy.
[0015]
Further, the step of creating a moving object model in the congestion prediction method according to the second feature includes temporarily generating a moving object newly arriving at the observation area using the unsteady Poisson process model, and creating a moving object model. Is also good. According to the congestion prediction method, highly accurate prediction can be performed in consideration of a change in an arrival pattern depending on a time zone to an observation area.
[0016]
A third feature of the present invention is to receive, via a communication network, transmitter information including a transmitter ID, which is transmitted from a plurality of receivers in an observation area where a mobile body carrying a transmitter moves. The congestion prediction system includes (a) a transmitter ID, an ID of a receiver that transmitted the transmitter information, a detection flag indicating arrival or departure of the transmitter with respect to a detection area of the receiver, and time information based on the transmitter information. (B) updating the moving body information represented by the past moving history of the moving body based on the moving event; and (c) using the moving body information to identify the moving body. A procedure for estimating the current position, and (d) when a congestion prediction is requested, a procedure for creating a moving body model used for predicting the behavior of the moving body using the moving body information at the time of the congestion prediction request And (ho Using a moving event generated based on the mobile body model, and summarized in that a congestion prediction program for executing a procedure for future congestion prediction.
[0017]
The procedure for creating a moving object model in the congestion prediction program according to the third feature includes an OD table indicating a moving state of the moving object between detection areas at the time of the congestion prediction request, and a detection at the time of the congestion prediction request. A moving object model may be created using a stay time table indicating the time spent in the detection area for each area.
[0018]
Further, the procedure for creating a moving object model in the congestion prediction program according to the third feature includes the step of temporarily generating a newly arrived moving object to the observation area using the non-stationary Poisson process model, and creating a moving object model. Is also good.
[0019]
By reading the congestion prediction program according to the third feature of the present invention, it is possible to cause a congestion prediction system or the like to execute the above procedure.
[0020]
BEST MODE FOR CARRYING OUT THE INVENTION
Next, an embodiment of the present invention will be described with reference to the drawings. In the following description of the drawings, the same or similar parts are denoted by the same or similar reference numerals. However, it should be noted that the drawings are schematic.
[0021]
The congestion prediction system according to the embodiment of the present invention updates a prediction model in real time from information obtained by observing a flow line of a moving body such as a pedestrian moving in a certain observation area, and updates future congestion. Make predictions. Also, the arrival of a pedestrian that does not yet exist at the start of the prediction is actually temporarily generated using the non-stationary Poisson process model, and is used for the prediction.
[0022]
(Receiver and transmitter in observation area)
Hereinafter, a system will be described in which a pedestrian moving in a certain observation area is taken as an example of a moving body, and a flow line of the pedestrian is observed to predict a future congestion. The
[0023]
Specifically, this
[0024]
The interval at which the
[0025]
On the other hand, each of the
[0026]
Next, a method of processing information received by the
[0027]
When the
[0028]
Next, the timing at which the
[0029]
1. When the detection flag changes to true
2. When the detection flag changes to false
The point in time when the detection flag changes to true indicates a point in time when the
[0030]
(Configuration of congestion prediction system)
As shown in FIG. 1, a
[0031]
The
[0032]
The external
[0033]
The recording
[0034]
The simulated
[0035]
The
[0036]
1. Monitor mode: Observes the current situation in the observation area. The movement event sent from the external
[0037]
2. Reproduction mode: Observes the situation at that time from information recorded in the observation area in the past. The movement event sent from the recording
[0038]
3. Simulation mode: A state of a virtual pedestrian in a virtual environment is predicted using a prediction model. The movement event sent from the
[0039]
Switching of the mode is performed via a user interface such as the
[0040]
The moving body
[0041]
Furthermore, the moving object
[0042]
ni= ni + 1
ai= ai + ln (xi + 1)
si= si + (ln (xi + 1)}Two
Only events that satisfy the same conditions as those in the OD table are recorded.
[0043]
The moving object
[0044]
The environment
[0045]
The
[0046]
The flow line for each pedestrian obtained by the
[0047]
Further, when the congestion prediction is requested, the
[0048]
Although not shown, the congestion prediction system according to the embodiment of the present invention may include a program holding unit for storing a program for updating a prediction model or executing a congestion prediction. The program holding unit may use an internal storage device such as a RAM or an external storage device such as an HD or FD.
[0049]
According to the congestion prediction system according to the embodiment of the present invention, the moving body model is updated in real time from information obtained by observing the trajectory of a pedestrian moving in a certain observation area, so that it is accurate and future Can be predicted.
[0050]
(Traffic line observation method)
Next, a flow line observation method will be described with reference to FIG. The
[0051]
(A) First, in step S101, the
[0052]
(B) Next, in step S102, the mobile body information and state buffer corresponding to the transmitter ID of the mobile event are searched from the mobile body information and state buffer held in the mobile body
[0053]
(C) Next, in step S103, if the mobile unit information and the state buffer cannot be searched, the process proceeds to step S104, where the mobile unit information and the state buffer for the transmitter ID are newly generated and newly registered. .
[0054]
(D) Next, in step S105, based on the movement event, the state buffer is updated in consideration of the environment information held in the environment
[0055]
(E) Next, in step S106, the position of the transmitter, that is, the area where the pedestrian is present is estimated based on the updated state buffer. Then, the estimated presence area is stored in the area portion or the like of the mobile object information.
[0056]
(F) Next, in step S107, if the status buffer or the existing area has changed in steps S105 and S106, the process proceeds to step S108, and the flow
[0057]
According to the flow line observing method according to the embodiment of the present invention, the position of a pedestrian can be accurately detected and the pedestrian's flow line can be observed without interruption.
[0058]
(Congestion prediction method)
Next, a congestion prediction method will be described with reference to FIGS. When a congestion prediction request is made, the
[0059]
(A) First, in step S201 in FIG. 7, a request for predicting future congestion is transmitted to the
[0060]
(B) In step S202, the monitor instance continues the flow line observation process. That is, according to the procedure shown in the flowchart of FIG. 6, as shown in FIG. 8, the movement event sent from the
[0061]
(C) On the other hand, in step S203, the simulation instance for prediction copies the moving body information, the OD table, and the residence time table at the time of the prediction request, and calculates the parameters of the moving
[0062]
(D) Next, in step S204, the
[0063]
(E) Next, in step S205, the simulation instance for prediction performs a prediction simulation after a lapse of a specified time based on the movement event. This prediction simulation will be described later in detail.
[0064]
According to the congestion prediction method according to the embodiment of the present invention, the moving body model is updated in real time from the information obtained by observing the flow line of the pedestrian moving in a certain observation area, so that it is accurate and future Can be predicted.
[0065]
(How to build a mobile model)
Here, a method of constructing a moving object model based on the moving object information, the OD table, the residence time table, and the like will be described.
[0066]
The moving object model used in the congestion prediction system according to the embodiment of the present invention includes the following sub-models.
[0067]
・ Visiting model: Determines the arrival time of pedestrians
・ Movement model: Select the venue to visit
・ Residence model: Determines the tour time at each venue
FIG. 9 shows a schematic diagram of the moving object model. Pedestrians arrive according to the arrival model, select and move to a venue according to the movement model until the time of leaving, and stay at each venue to be observed according to the retention model. Hereinafter, an example of constructing each model will be described. Here, the flow line of the pedestrian when observing the exhibition hall as shown in FIG. 10 is predicted.
[0068]
<Visiting model>
Since it is inappropriate to assume that the visitor's visitor density is constant, a non-stationary Poisson process (NHPP) is introduced into the formulation. Details of NHPP are described in "DP Heyman, MJ Sobel, Eds .:" Stochastic Models ", Handbooks in Operations Research and Management Science, Vol. 2, Elsevier Science Publishers BV (North-Holland), 1990." . According to the NHPP, a Poisson arrival model having a time-varying visit density can be created. The next visit probability distribution at time T is given by the following equation, where λ is the visit density.
[0069]
(Equation 1)
The time from the time T to the next visit time may be obtained by solving this equation for t. Equation (1) is modified to introduce a uniform random number u to obtain the following equation.
[0070]
(Equation 2)
The visits are divided into morning and afternoon, and each visit time distribution is approximated by a normal distribution. According to this, the visit density λ (t) is represented by the following equation, where μ is the average visit time, σ is the standard deviation, AM is AM, and PM is PM.
[0071]
(Equation 3)
For the identification of the visit model, for example, the following visit model is obtained by fitting a normal distribution to the actual visit time distribution in the morning and afternoon. The average value indicates the elapsed time when 9:00 is set to 0 second.
[0072]
μAM= 3855, σAM= 2703
μPM= 18838, σPM= 3573 ……… (4)
To estimate the visitor's arrival time, the total number of visitors in the morning NAMAnd N in the afternoonPMIs required, but will be given separately.
[0073]
<Moving model>
The visitor's movement model is formulated by combining the following two models.
[0074]
・ Exit model
・ Next venue selection model
That is, when a visitor finishes a tour at a certain venue, it is first determined whether or not to leave using a leaving model. If it is decided not to leave, the next venue to be visited is determined using the next venue selection model.
[0075]
The exit model uses the exit time as tL, Visit timeA, Stay time tS, TCIs defined by the following equation.
[0076]
(Equation 4)
DLIf (t) is true, exit is selected. In the case of false, the next destination venue is determined by the next venue selection model. The stay time is determined by a probability distribution (normal distribution) defined by the following equation.
[0077]
(Equation 5)
Average value μSAnd standard deviation σSIs the visitor's arrival time tAIs given as a linear function of. This expresses the tendency of visitors in the morning to complete the tour in the morning and visitors in the afternoon to close by the closing time.
[0078]
(Equation 6)
Exit time tLDoes not depend on the state of the visitors.
[0079]
The next venue selection model is formulated using a Markov chain model. Therefore, the model does not consider the past visit history. The details of the Markov chain model are described in Shuji Kaneko, "Simulation of the flow of people in an exhibition", Advance, No. 6, pp. 33-37, Toshiba Advanced System, 1998.
[0080]
X represents a Markov chain representing each visitor's movement between venueskAnd Assume that the reception for the visit is one venue, and the number of exhibition venues including the reception is M. The reception is called
[0081]
(Equation 7)
This can be written as an M-dimensional vector as follows.
[0082]
p [k] = [p1[k], pTwo[k],…, pM[k]] ……… (9)
The initial distribution of all visitors is always p [0] = [1,0, ..., 0]. The transition probability matrix is represented by an M × M-dimensional matrix P. The i-th row and j-th column element of the transition probability matrix P is the transition probability from the venue i to the venue j.
[0083]
pij = P (Xk + 1= bj | Xk = bi) ……… (10)
The probability distribution of each visitor at time k is represented by the following equation.
[0084]
p [k] = p [0] Pk ……… (11)
In this model, the transition probability matrix P is regarded as invariable during the target exhibition.
[0085]
As the identification of the leaving model, the closing time TCIs 17:00, and the number of seconds from the venue time 9:00 is as follows.
[0086]
TC = 28800 ……… (12)
Stay time tSIs obtained as follows.
[0087]
1. The samples are divided into two groups, morning and afternoon, according to their arrival times.
[0088]
2. Exclude all-day visitors from the morning group.
[0089]
3. For both groups, calculate the average value and standard deviation of the stay time every hour.
[0090]
4. For both groups, the time change of the mean value and the standard deviation is identified by a linear function. At this time, the following samples are added to the afternoon group.
[0091]
μS(TC) = 0, σS(TC) = 0 (13)
The calculation of the standard deviation is a fixed value based on the actual average. As a result, for example, the stay time tSIs obtained as a normal distribution represented by the following parameters (unit is seconds).
[0092]
αAM = -0.3521, αPM = -0.6696
βAM = 9752.4, βPM = 19904.4
ξAM = 0, ξPM = -0.2233
ζAM = 3704.4, ξPM = 6508.8 ……… (14)
To identify the next venue selection model, the transition probability matrix P is derived from the actually measured OD table in the following procedure.
[0093]
1. The S column is deleted to make a 13 × 12 dimensional matrix.
[0094]
2. Normalize each row.
[0095]
3. Transpose the matrix.
[0096]
4. One row is inserted in the first row to make a 13 × 13-dimensional square matrix. The row vector to be inserted is [0, ..., 0].
[0097]
FIG. 11 shows the processing result.
[0098]
<Residence model>
The staying time at each venue is formalized as following a certain probability distribution. It is assumed that the logarithmic conversion result of the stay time follows a normal distribution. That is, the stay time T at the venue iSiIs given by the following equation.
[0099]
(Equation 8)
xiIs determined by the probability distribution defined by the following equation.
[0100]
(Equation 9)
FIG. 12 shows an example of the result of logarithmic conversion of the staying time record at each venue and fitting to a normal distribution for each. FIG. 13 is obtained by inversely converting the average value.
[0101]
With reference to the visit model, the movement model, and the staying model created as described above, a moving object model is constructed.
[0102]
(Predictive simulation)
Based on the moving object model obtained above, a moving event is generated, and the congestion of the venue is predicted. The algorithm for reproducing the congestion of the venue by simulation is shown below.
[0103]
1. Set the number of visitors in the morning and afternoon as N_AM and N_PM, respectively.
[0104]
2. Time initialization t = 0 (9:00 am); visitor counter initialization
i = 0
3.while (i <N_AM) {
4. i = i + 1
5. Visitor i occurs according to the visit model; time t_i = t + occurrence interval
6. Departure time t_L [i] of visitor i is determined by the exit model
7. while (t_i <t_L [i]) {
8. Determine the visiting venue by the next venue selection model and move
9. Determination of residence time t_S by residence model; t_i = t_i + t_S
Ten. }
11. Exit processing
12.}
13. Perform
14. Sort the behavior data of all visitors in chronological order
15. Aggregation processing
The simulation program can be created by Java (registered trademark) or the like.
[0105]
The pedestrian information, the OD table, the residence time table, and the like obtained by the simulation may be stored in the moving object
[0106]
(Other embodiments)
Although the present invention has been described with the above embodiments, it should not be understood that the description and drawings forming part of this disclosure limit the present invention. From this disclosure, various alternative embodiments, examples, and operation techniques will be apparent to those skilled in the art.
[0107]
For example, in the embodiment of the present invention, the moving body moving in the observation area is described as a pedestrian, but is not limited to this.
[0108]
Also, in the embodiment of the present invention, it has been described that the time information is added by the external
[0109]
Further, in the embodiment of the present invention, it has been described that the moving object
[0110]
As described above, the present invention naturally includes various embodiments and the like not described herein. Therefore, the technical scope of the present invention is determined only by the invention specifying matters according to the claims that are appropriate from the above description.
[0111]
【The invention's effect】
According to the present invention, it is possible to provide a congestion prediction system, a congestion prediction method, and a congestion prediction program for updating a prediction model in real time and accurately predicting future congestion.
[Brief description of the drawings]
FIG. 1 is a configuration block diagram of a congestion prediction system according to an embodiment of the present invention.
FIG. 2 is a configuration block diagram of a wireless tag provided in the transmitter according to the embodiment of the present invention.
FIG. 3 is an example of transmitter information according to the embodiment of the present invention.
FIG. 4 is an example of moving object information according to the embodiment of the present invention.
FIG. 5 is an example of an observation area according to the embodiment of the present invention.
FIG. 6 is a flowchart of a flow line observing method according to the embodiment of the present invention.
FIG. 7 is a flowchart of a congestion prediction method according to the embodiment of the present invention.
FIG. 8 is a diagram showing two simulation instances of the congestion prediction method according to the embodiment of the present invention.
FIG. 9 is a conceptual diagram of a moving object model according to the embodiment of the present invention.
FIG. 10 is an example of a layout of an exhibition hall.
FIG. 11 is a transition probability matrix used when constructing the moving object model according to the embodiment of the present invention.
FIG. 12 is a stay time distribution used for identifying a staying model according to the embodiment of the present invention.
FIG. 13 is a diagram showing an average staying time of the staying model according to the embodiment of the present invention.
FIG. 14 is an example of a visit time distribution obtained by the congestion prediction method according to the embodiment of the present invention.
FIG. 15 is an example of a stay time distribution obtained by the congestion prediction method according to the embodiment of the present invention.
FIG. 16 is an example of a congestion transition obtained by the congestion prediction method according to the embodiment of the present invention.
FIG. 17 is an example of an OD table obtained by the congestion prediction method according to the embodiment of the present invention.
[Explanation of symbols]
1 Congestion prediction system
5 Record files
6 Traffic flow record file
7 Display
8 Moving object model
9 Input section
11 Recorded event capture unit
12 External event capture unit
13 Receiver control unit
14 Input switching unit
15 Moving body state holding unit
16 Event processing part
17 Environmental Information Storage Department
18 Simulated event generator
20a, 20b, ..., 20e Receiver
21 transmitter
22 pedestrian
25 Communication Network
26 Observation area
27 Partition
30 wireless tags
31 Identification information holding unit
32 transmission unit
33 Antenna
Claims (9)
前記発信機情報に基づき、前記発信機ID、前記発信機情報を送信した受信機のID、前記受信機の検出エリアに対する前記発信機の到着あるいは離脱を示す検出フラグ及び時刻情報を含む移動イベントを生成する外部イベント取込部と、
前記移動体の行動を予測するために用いる移動体モデルに基づいて、前記移動イベントを生成する模擬イベント発生部と、
前記外部イベント取込部と前記模擬イベント発生部から送られる前記移動イベントを切り替えて出力する入力切替部と、
前記入力切替部を介して前記外部イベント取込部から送られた移動イベントを順次取り出し、前記移動体の過去の移動履歴で表現される移動体情報を更新し、前記移動体の現在位置の推定を行うことと平行して、混雑予測を要求された時には、前記入力切替部を介して前記模擬イベント発生部から送られた移動イベントを用い、将来の混雑予測を行うイベント処理部と
を備えることを特徴とする混雑予測システム。Via a communication network, a receiver control unit that receives transmitter information including an ID of the transmitter, which is transmitted from a plurality of receivers in an observation area in which a moving object carrying the transmitter moves,
Based on the transmitter information, a mobile event including the transmitter ID, the ID of the receiver that transmitted the transmitter information, a detection flag indicating the arrival or departure of the transmitter with respect to the detection area of the receiver, and time information including the time information. An external event capture unit to generate,
A simulated event generator that generates the movement event based on a moving body model used to predict the behavior of the moving body,
An input switching unit that switches and outputs the movement event sent from the external event capture unit and the simulation event generation unit,
The moving events sent from the external event capturing unit are sequentially taken out via the input switching unit, and the moving object information represented by the past moving history of the moving object is updated, and the current position of the moving object is estimated. In parallel with performing a congestion prediction, an event processing unit that performs a future congestion prediction using a movement event sent from the simulation event generating unit via the input switching unit when a congestion prediction is requested. A congestion prediction system characterized by the following.
前記発信機情報に基づき、前記発信機ID、前記発信機情報を送信した受信機のID、前記受信機の検出エリアに対する前記発信機の到着あるいは離脱を示す検出フラグ及び時刻情報を含む移動イベントを順次取り出すステップと、
前記移動イベントに基づいて、前記移動体の過去の移動履歴で表現される移動体情報を更新するステップと、
前記移動体情報を用いて、前記移動体の現在位置の推定を行うステップと、
混雑予測を要求された時には、混雑予測要求時点での前記移動体情報を用いて、前記移動体の行動を予測するために用いる移動体モデルを作成するステップと、
前記移動体モデルに基づいて生成された移動イベントを用い、将来の混雑予測を行うステップと
を含むことを特徴とする混雑予測方法。In a congestion prediction system that receives transmitter information including an ID of the transmitter, which is transmitted from a plurality of receivers in an observation area where a mobile body carrying the transmitter moves via a communication network,
Based on the transmitter information, a mobile event including the transmitter ID, the ID of the receiver that transmitted the transmitter information, a detection flag indicating the arrival or departure of the transmitter with respect to the detection area of the receiver, and time information including the time information. Sequentially taking out;
Based on the movement event, updating the moving body information expressed in the past movement history of the moving body,
Estimating a current position of the moving object using the moving object information;
When a congestion prediction is requested, using the mobile information at the time of the congestion prediction request, creating a mobile model used to predict the behavior of the mobile,
Performing a future congestion prediction using a movement event generated based on the moving object model.
前記発信機情報に基づき、前記発信機ID、前記発信機情報を送信した受信機のID、前記受信機の検出エリアに対する前記発信機の到着あるいは離脱を示す検出フラグ及び時刻情報を含む移動イベントを順次取り出す手順と、
前記移動イベントに基づいて、前記移動体の過去の移動履歴で表現される移動体情報を更新する手順と、
前記移動体情報を用いて、前記移動体の現在位置の推定を行う手順と、
混雑予測を要求された時には、混雑予測要求時点での前記移動体情報を用いて、前記移動体の行動を予測するために用いる移動体モデルを作成する手順と、
前記移動体モデルに基づいて生成された移動イベントを用い、将来の混雑予測を行う手順と
を実行させることを特徴とする混雑予測プログラム。Via a communication network, transmitted from a plurality of receivers in the observation area where the mobile body carrying the transmitter moves, the congestion prediction system to receive the transmitter information including the ID of the transmitter,
Based on the transmitter information, a mobile event including the transmitter ID, the ID of the receiver that transmitted the transmitter information, a detection flag indicating the arrival or departure of the transmitter with respect to the detection area of the receiver, and time information including the time information. Steps to take out sequentially
Based on the movement event, a procedure for updating moving body information expressed in a past movement history of the moving body,
A procedure of estimating a current position of the moving body using the moving body information;
When congestion prediction is requested, using the mobile information at the time of congestion prediction request, a procedure to create a mobile model used to predict the behavior of the mobile,
A congestion prediction program for performing prediction of future congestion using a movement event generated based on the moving object model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002378589A JP2004213098A (en) | 2002-12-26 | 2002-12-26 | Congestion prediction system, congestion prediction method and congestion prediction program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002378589A JP2004213098A (en) | 2002-12-26 | 2002-12-26 | Congestion prediction system, congestion prediction method and congestion prediction program |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2004213098A true JP2004213098A (en) | 2004-07-29 |
Family
ID=32815375
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002378589A Pending JP2004213098A (en) | 2002-12-26 | 2002-12-26 | Congestion prediction system, congestion prediction method and congestion prediction program |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2004213098A (en) |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2006025392A1 (en) * | 2004-08-31 | 2006-03-09 | Osaka University | Mobile node simulator and program for mounting the same |
| JP2006106875A (en) * | 2004-09-30 | 2006-04-20 | Nifty Corp | System for predicting movement of visitor in closed space |
| WO2011024379A1 (en) * | 2009-08-24 | 2011-03-03 | ソフトバンクBb株式会社 | Population mobility estimation system, population mobility estimation method, and population mobility estimation program |
| JP2014112292A (en) * | 2012-12-05 | 2014-06-19 | Hitachi Ltd | Human flow survey support system and method |
| JP2016110360A (en) * | 2014-12-05 | 2016-06-20 | 株式会社Nttドコモ | Congestion prediction device, congestion prediction system and congestion prediction method |
| JP2017021855A (en) * | 2016-10-19 | 2017-01-26 | 株式会社Agoop | Information processing system, population flow rate estimating device, program, information processing method, and population flow rate estimating method |
| JP2017146829A (en) * | 2016-02-18 | 2017-08-24 | 日本電信電話株式会社 | Distribution estimation device, distribution estimation method, and distribution estimation program |
| JP2018073236A (en) * | 2016-11-01 | 2018-05-10 | 日本電信電話株式会社 | Human flow rate prediction device, human flow rate prediction method, and human flow rate prediction program |
| JP2018112774A (en) * | 2017-01-06 | 2018-07-19 | 株式会社アルテ | Space-time data mining apparatus |
| KR101897513B1 (en) * | 2017-12-18 | 2018-09-12 | 가천대학교 산학협력단 | Method, apparatus and computer-readable recording medium for selecting feature using relational matrix |
| JP2018195215A (en) * | 2017-05-19 | 2018-12-06 | 日本電信電話株式会社 | Human flow prediction device, human flow prediction method, and human flow prediction program |
| JP2019082757A (en) * | 2017-10-27 | 2019-05-30 | 日本電信電話株式会社 | Estimation device, estimation method, and program |
| WO2019189207A1 (en) * | 2018-03-27 | 2019-10-03 | 日本電信電話株式会社 | Estimation device, estimation method, and program |
| JPWO2020261450A1 (en) * | 2019-06-26 | 2020-12-30 | ||
| JPWO2020261451A1 (en) * | 2019-06-26 | 2020-12-30 | ||
| JP2021002252A (en) * | 2019-06-24 | 2021-01-07 | 株式会社日立製作所 | Traffic flow prediction support device, traffic flow prediction support method, and traffic flow prediction support program |
| CN113271375A (en) * | 2021-04-28 | 2021-08-17 | 杭州智科飞创信息科技有限公司 | Scientific and technological innovation park operation management method |
-
2002
- 2002-12-26 JP JP2002378589A patent/JP2004213098A/en active Pending
Cited By (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2006025392A1 (en) * | 2004-08-31 | 2006-03-09 | Osaka University | Mobile node simulator and program for mounting the same |
| JP2006106875A (en) * | 2004-09-30 | 2006-04-20 | Nifty Corp | System for predicting movement of visitor in closed space |
| WO2011024379A1 (en) * | 2009-08-24 | 2011-03-03 | ソフトバンクBb株式会社 | Population mobility estimation system, population mobility estimation method, and population mobility estimation program |
| JP2014112292A (en) * | 2012-12-05 | 2014-06-19 | Hitachi Ltd | Human flow survey support system and method |
| JP2016110360A (en) * | 2014-12-05 | 2016-06-20 | 株式会社Nttドコモ | Congestion prediction device, congestion prediction system and congestion prediction method |
| JP2017146829A (en) * | 2016-02-18 | 2017-08-24 | 日本電信電話株式会社 | Distribution estimation device, distribution estimation method, and distribution estimation program |
| JP2017021855A (en) * | 2016-10-19 | 2017-01-26 | 株式会社Agoop | Information processing system, population flow rate estimating device, program, information processing method, and population flow rate estimating method |
| JP2018073236A (en) * | 2016-11-01 | 2018-05-10 | 日本電信電話株式会社 | Human flow rate prediction device, human flow rate prediction method, and human flow rate prediction program |
| JP2018112774A (en) * | 2017-01-06 | 2018-07-19 | 株式会社アルテ | Space-time data mining apparatus |
| JP2018195215A (en) * | 2017-05-19 | 2018-12-06 | 日本電信電話株式会社 | Human flow prediction device, human flow prediction method, and human flow prediction program |
| JP2019082757A (en) * | 2017-10-27 | 2019-05-30 | 日本電信電話株式会社 | Estimation device, estimation method, and program |
| KR101897513B1 (en) * | 2017-12-18 | 2018-09-12 | 가천대학교 산학협력단 | Method, apparatus and computer-readable recording medium for selecting feature using relational matrix |
| WO2019189207A1 (en) * | 2018-03-27 | 2019-10-03 | 日本電信電話株式会社 | Estimation device, estimation method, and program |
| JP2019174997A (en) * | 2018-03-27 | 2019-10-10 | 日本電信電話株式会社 | Estimation device, estimation method, and program |
| US11859992B2 (en) | 2018-03-27 | 2024-01-02 | Nippon Telegraph And Telephone Corporation | Estimation apparatus, estimation method and program |
| JP2021002252A (en) * | 2019-06-24 | 2021-01-07 | 株式会社日立製作所 | Traffic flow prediction support device, traffic flow prediction support method, and traffic flow prediction support program |
| JP7372058B2 (en) | 2019-06-24 | 2023-10-31 | 株式会社日立製作所 | Traffic flow prediction support device, traffic flow prediction support method, and traffic flow prediction support program |
| JPWO2020261450A1 (en) * | 2019-06-26 | 2020-12-30 | ||
| JPWO2020261451A1 (en) * | 2019-06-26 | 2020-12-30 | ||
| CN113271375A (en) * | 2021-04-28 | 2021-08-17 | 杭州智科飞创信息科技有限公司 | Scientific and technological innovation park operation management method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2004213098A (en) | Congestion prediction system, congestion prediction method and congestion prediction program | |
| Salimi et al. | Occupancy prediction model for open-plan offices using real-time location system and inhomogeneous Markov chain | |
| Liebig et al. | Dynamic route planning with real-time traffic predictions | |
| CN112470199A (en) | System and method for point-to-point traffic prediction | |
| US10169993B1 (en) | Forecasting with matrix powers | |
| JP6761789B2 (en) | People flow forecasters, systems and programs | |
| JP6939911B2 (en) | Methods and devices for adaptive vehicle control | |
| JP2011059924A5 (en) | ||
| EP3662423A1 (en) | System and method for predictive cleaning | |
| JP2012003494A (en) | Information processing device, information processing method and program | |
| KR20190062443A (en) | Passive visit detection | |
| JP2009133849A (en) | Method, computer program, and data processing system for deciding estimated position in space time of rfid (radio-frequency identification) tag | |
| US8519850B2 (en) | Method and system for the acquisition, transmission and assessment of remote sensor data for trend analysis, prediction and remediation | |
| CN104869581A (en) | Determining Faulty Nodes via Label Propagation in Wireless Sensor Networks | |
| US12210987B2 (en) | System and method for controlling motion of a bank of elevators | |
| JP2018008758A (en) | Elevator system and car call estimation method | |
| JP2004214737A (en) | Flow line display system, flow line display method, and flow line display program | |
| CN120706577A (en) | A method for predicting community resident behavior based on IoT pre-processing | |
| KR102377148B1 (en) | Smart quarantine grid system and method thereof | |
| JP7380921B1 (en) | Facility guidance interlocking building system, information processing device, facility guidance method, and computer-readable recording medium | |
| Sun et al. | Real-time bus travel time prediction and reliability quantification: A hybrid Markov model | |
| Bawaneh et al. | Machine Learning-Based Anomaly Detection in Smart City Traffic: Performance Comparison and Insights. | |
| JP2007279887A (en) | Singular pattern detection system, model learning device, singular pattern detection device, singular pattern detection method, and computer program | |
| JP7047174B1 (en) | Prediction system, prediction device, prediction method and prediction program | |
| JP2008211781A (en) | System for modeling movement of objects in certain environment and method executed by computer |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20040604 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20061026 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20061031 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20061228 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20080115 |