(57)【特許請求の範囲】
【請求項１】
移動機が第１周期で登録する第１位置情報に基づき、所定エリア内の移動機数を集計する第１集計手段と、
移動機が第１周期よりも短い第２周期で登録する第２位置情報に基づき、前記所定エリア内の移動機数を集計する第２集計手段と、
前記第１集計手段で集計された移動機のうち、前記第２集計手段で集計された移動機と一致する移動機数を集計する共通部集計手段と、
前記第２集計手段により集計された移動機数および前記共通部集計手段により集計された移動機数を用い、前記第１集計手段により集計された移動機数に対して第１周期よりも短く設定された期間における移動機数の増減分を調整する時間粒度詳細化係数を算出し、前記第１集計手段により集計された移動機数に前記時間粒度詳細化係数を乗算して人口を算出する人口算出手段と、
を備える人口算出装置。
【請求項２】
移動機が第１周期で登録する第１位置情報に基づき、所定エリア内の移動機数を集計する第１集計手段と、
移動機が第１周期よりも短い第２周期で登録する第２位置情報に基づき、前記所定エリア内の移動機数を集計する第２集計手段と、
前記第１集計手段で集計された移動機のうち、前記第２集計手段で集計された移動機と一致する移動機数を集計する共通部集計手段と、
前記第１集計手段により集計された移動機数および前記共通部集計手段により集計された移動機数を用い、前記第１位置情報に対する前記第２位置情報の利用率を算出し、前記第２集計手段により集計された移動機数を前記利用率で除算して人口を算出する人口算出手段と、
を備える人口算出装置。
【請求項３】
前記人口算出手段は、前記第１集計手段により集計された移動機数、前記第２集計手段により集計された移動機数、および前記共通部集計手段により集計された移動機数を用いてユーザ数を算出し、前記ユーザ数を人口に変換するための拡大係数を用いて、前記ユーザ数を人口に変換することで人口を算出することを特徴とする請求項１又は２記載の人口算出装置。
【請求項４】
前記第１集計手段は、所定期間に前記所定エリア内で登録された前記第１位置情報のそれぞれについて、移動機が前記第１位置情報を登録した時刻、同じ移動機が一つ前の前記第１位置情報を登録した時刻、及び同じ移動機が一つ後の前記第１位置情報を登録した時刻、のうち二つ以上を用いて第１特徴量を算出し、前記第１特徴量の総和に基づいて前記所定エリア内の移動機数を集計し、
前記第２集計手段は、前記所定期間に前記所定エリア内で登録された前記第２位置情報のそれぞれについて、移動機が前記第２位置情報を登録した時刻、同じ移動機が一つ前の前記第２位置情報を登録した時刻、及び同じ移動機が一つ後の前記第２位置情報を登録した時刻、のうち二つ以上を用いて第２特徴量を算出し、前記第２特徴量の総和に基づいて前記所定エリア内の移動機数を集計し、
前記共通部集計手段は、前記所定期間に前記所定エリア内で前記第１位置情報を登録した移動機によって、前記所定期間に前記所定エリア内で登録された前記第２位置情報のそれぞれについて、移動機が前記第２位置情報を登録した時刻、同じ移動機が一つ前の前記第２位置情報を登録した時刻、及び同じ移動機が一つ後の前記第２位置情報を登録した時刻、のうち二つ以上を用いて共通特徴量を算出し、共通特徴量の総和に基づいて、前記第１集計手段で集計された移動機のうち、前記第２集計手段で集計された移動機と一致する移動機数を集計することを特徴とする請求項１〜３のいずれか一項記載の人口算出装置。
【請求項５】
人口算出装置が行う人口算出方法において、
移動機が第１周期で登録する第１位置情報に基づき、所定エリア内の移動機数を集計する第１集計ステップと、
移動機が第１周期よりも短い第２周期で登録する第２位置情報に基づき、前記所定エリア内の移動機数を集計する第２集計ステップと、
前記第１集計ステップで集計された移動機のうち、前記第２集計ステップで集計された移動機と一致する移動機数を集計する共通部集計ステップと、
前記第２集計ステップにて集計された移動機数および前記共通部集計ステップにて集計された移動機数を用い、前記第１集計ステップにて集計された移動機数に対して第１周期よりも短く設定された期間における移動機数の増減分を調整する時間粒度詳細化係数を算出し、前記第１集計ステップにて集計された移動機数に前記時間粒度詳細化係数を乗算して人口を算出する人口算出ステップと、
を備える人口算出方法。
【請求項６】
人口算出装置が行う人口算出方法において、
移動機が第１周期で登録する第１位置情報に基づき、所定エリア内の移動機数を集計する第１集計ステップと、
移動機が第１周期よりも短い第２周期で登録する第２位置情報に基づき、前記所定エリア内の移動機数を集計する第２集計ステップと、
前記第１集計ステップで集計された移動機のうち、前記第２集計ステップで集計された移動機と一致する移動機数を集計する共通部集計ステップと、
前記第１集計ステップにて集計された移動機数および前記共通部集計ステップにて集計された移動機数を用い、前記第１位置情報に対する前記第２位置情報の利用率を算出し、前記第２集計ステップにて集計された移動機数を前記利用率で除算して人口を算出する人口算出ステップと、
を備える人口算出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、移動機数に基づいて人口を算出する人口算出装置および人口算出方法に関する。
【背景技術】
【０００２】
移動機の位置登録数に基づいて移動機数を算出し、これを人口数と把握することが知られている。例えば、特許文献１（特開２００３−１２２８７７号公報）には、携帯電話の位置情報ＤＢを利用して携帯電話数を算出し、これを所定領域の人口として算出することが記載されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００３−１２２８７７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、上述従来の技術においては、短時間における急激な人数の増減を把握することができないという問題があった。即ち、従来の技術が利用する携帯電話の位置登録処理は、５４分周期で行われることが通信規格上決まっているため、５４分以内に移動する携帯電話数を把握できなかった。
【０００５】
本発明は、このような問題点を解決するためになされたものであり、短時間における急激な人数の増減を把握できる人口算出装置および人口算出方法の提供を目的とする。
【課題を解決するための手段】
【０００６】
本発明による人口算出装置は、移動機が第１周期で登録する第１位置情報に基づき、所定エリア内の移動機数を集計する第１集計手段と、移動機が第１周期よりも短い第２周期で登録する第２位置情報に基づき、所定エリア内の移動機数を集計する第２集計手段と、第１集計手段で集計された移動機のうち、第２集計手段で集計された移動機と一致する移動機数を集計する共通部集計手段と、第１集計手段により集計された移動機数、第２集計手段により集計された移動機数、および共通部集計手段により集計された移動機数を用い所定エリア内の人口を算出する人口算出手段と、を備える。
【０００７】
また、本発明による人口算出方法は、人口算出装置が行う人口算出方法において、移動機が第１周期で登録する第１位置情報に基づき、所定エリア内の移動機数を集計する第１集計ステップと、移動機が第１周期よりも短い第２周期で登録する第２位置情報に基づき、所定エリア内の移動機数を集計する第２集計ステップと、第１集計ステップで集計された移動機のうち、第２集計ステップで集計された移動機と一致する移動機数を集計する共通部集計ステップと、第１集計ステップにより集計された移動機数、第２集計ステップにより集計された移動機数、および共通部集計ステップにより集計された移動機数を用い所定エリア内の人口を算出する人口算出ステップと、を備える。
【０００８】
この発明によれば、第１集計手段により集計された移動機数、第２集計手段により集計された移動機数、および共通部集計手段により集計された移動機数が用いられ、所定エリア内の人口が算出される。そして、第１周期で登録される第１位置情報と、第２周期で登録される第２位置情報とが互いに補われる。これにより、例えば、第１位置情報のみでは急激な人口の増減をとらえられない一方で、第２位置情報を登録する移動機数が少なく第２位置情報のみでは人口を算出できない場合であっても、第１位置情報と第２位置情報とが互いに補われ、急激な人口の増減を把握すること可能となる。
【０００９】
ここで、上記作用を効果的に奏する構成として、人口算出手段は、第２集計手段により集計された移動機数および共通部集計手段により集計された移動機数を用い、第１集計手段により集計された移動機数に対して第１周期よりも短く設定された期間における移動機数の増減分を調整する時間粒度詳細化係数を算出し、第１集計手段により集計された移動機数に時間粒度詳細化係数を乗算して人口を算出する構成が挙げられる。
【００１０】
また、上記作用を効果的に奏する構成として、人口算出手段は、第１集計手段により集計された移動機数および共通部集計手段により集計された移動機数を用い、第１位置情報に対する第２位置情報の利用率を算出し、第２集計手段により集計された移動機数を利用率で除算して人口を算出する構成が挙げられる。
【００１１】
また、上記作用を効果的に奏する構成として、人口算出手段は、第１集計手段により集計された移動機数、第２集計手段により集計された移動機数、および共通部集計手段により集計された移動機数を用いてユーザ数を算出し、ユーザ数を人口に変換するための拡大係数を用いて、ユーザ数を人口に変換することで人口を算出する構成が挙げられる。
【００１２】
また、上記作用を効果的に奏する構成として、第１集計手段は、所定期間に所定エリア内で登録された第１位置情報のそれぞれについて、移動機が第１位置情報を登録した時刻、同じ移動機が一つ前の第１位置情報を登録した時刻、及び同じ移動機が一つ後の第１位置情報を登録した時刻、のうち二つ以上を用いて第１特徴量を算出し、第１特徴量の総和に基づいて所定エリア内の移動機数を集計し、第２集計手段は、所定期間に所定エリア内で登録された第２位置情報のそれぞれについて、移動機が第２位置情報を登録した時刻、同じ移動機が一つ前の第２位置情報を登録した時刻、及び同じ移動機が一つ後の第２位置情報を登録した時刻、のうち二つ以上を用いて第２特徴量を算出し、第２特徴量の総和に基づいて所定エリア内の移動機数を集計し、共通部集計手段は、所定期間に所定エリア内で第１位置情報を登録した移動機によって、所定期間に所定エリア内で登録された第２位置情報のそれぞれについて、移動機が第２位置情報を登録した時刻、同じ移動機が一つ前の第２位置情報を登録した時刻、及び同じ移動機が一つ後の第２位置情報を登録した時刻、のうち二つ以上を用いて共通特徴量を算出し、共通特徴量の総和に基づいて、第１集計手段で集計された移動機のうち、第２集計手段で集計された移動機と一致する移動機数を集計する構成が挙げられる。
【発明の効果】
【００１３】
本発明によれば、短時間における急激な人数の増減を把握できる。
【図面の簡単な説明】
【００１４】
【図１】本発明に係る人口算出装置の第１実施形態を含む通信システムのシステム構成図である。
【図２】本発明に係る人口算出装置の第１実施形態の機能構成を示すブロック図である。
【図３】図２の人口算出装置の人口算出手順を示すフローチャートである。
【図４】人口算出対象のエリアへの移動機の出入りを示す模式図である。
【図５】長周期データ集計部が有する管理テーブルを示す説明図である。
【図６】短周期データ集計部が有する管理テーブルを示す説明図である。
【図７】本発明に係る人口算出装置の第２実施形態の人口算出手順を示すフローチャートである。
【図８】関連セクタを例示する図である。
【図９】本発明に係る人口算出装置の第３実施形態の人口算出手順を示すフローチャートである。
【図１０】端末数推計の考え方を説明するための図である。
【図１１】端末数推計に係る計算方法を説明するための図である。
【発明を実施するための形態】
【００１５】
以下、本発明に係る人口算出装置の好適な実施形態について添付図面を参照しながら説明する。
【００１６】
図１は、本発明に係る人口算出装置の第１実施形態を含む通信システムのシステム構成図である。図１に示すように、本実施形態の人口算出装置を含む通信システム１０は、移動機１００、ＢＴＳ（基地局）２００、ＲＮＣ（無線制御装置）３００、交換機４００、および管理センタ５００を含んで構成されている。また、この管理センタ５００は、社会センサユニット５０１、ペタマイニングユニット５０２、モバイルデモグラフィユニット５０３、および可視化ソリューションユニット５０４を含んで構成されている。
【００１７】
移動機１００は、第１周期（ここでは５４分周期）で第１位置情報の登録処理を行う。当該第１位置情報は、ＢＴＳ２００、ＲＮＣ３００、交換機４００を介し、後述の社会センサユニット５０１に蓄積される。なお、第１位置情報は、ＢＴＳ２００の通信エリアを複数に等分したセクタや基地局の単位で表され、移動機１００の在圏しているセクタや基地局のＩＤが第１位置情報として記録される。また、第１位置情報は、ＢＴＳ２００などの提供する機能（例えばＰＲＡＣＨ  ＰＤなど）や移動機が有する機能（ＧＰＳなど）によって測位された緯度経度情報でもよい。
【００１８】
また、一部の移動機１００は、上記第１周期よりも短い第２周期（ここでは５分周期）で第２位置情報の登録処理を行う。そして、当該第２位置情報は、ＢＴＳ２００、ＲＮＣ３００、交換機４００を介し、後述の社会センサユニット５０１に蓄積される。なお、第２位置情報はＢＴＳ２００などの提供する機能（例えばＰＲＡＣＨ  ＰＤなど）や移動機が有する機能（ＧＰＳなど）を用いて取得され、移動機１００の位置する緯度および経度が第２位置情報として記録される。また、第２位置情報は、ＢＴＳ２００の通信エリアを複数に等分したセクタや基地局の単位で表され、移動機１００の在圏しているセクタや基地局のＩＤが第２位置情報として記録されてもよい。
【００１９】
社会センサユニット５０１は、移動機１００の位置情報等を含んだデータを蓄積するサーバ装置である。この社会センサユニット５０１は、交換機４００からのデータ受信要求に応じ、または交換機４００にデータ送信要求を出すことでデータを取得し、蓄積する。
【００２０】
ペタマイニングユニット５０２は、社会センサユニット５０１から受信したデータを所定のデータ形式に変換するサーバ装置である。例えば、ペタマイニングユニット５０２は、識別子をキーにソーティング処理を行ったり、エリアごとにソーティング処理を行ったりする。本実施形態による人口算出装置６００は、ペタマイニングユニット５０２内に構成されており、人口算出対象（ここでは人口算出対象エリアおよび時間帯）の人口を移動機１００の数から算出する。なお、人口算出装置６００は、定期的な算出指示により、所定時間帯における所定エリア内の人口を定期的に算出させるように用いることが可能である。また、人口算出装置６００は、外部インターフェースから適宜算出指示を入力し、所望の時間帯における所望のエリア内の人口を算出させるように用いることも可能である。
【００２１】
モバイルデモグラフィユニット５０３は、ペタマイニングユニット５０２において処理されたデータに対する集計・統計解析処理、空間解析処理を行うサーバ装置である。
【００２２】
可視化ソリューションユニット５０４は、モバイルデモグラフィユニット５０３において集計処理されたデータを可視化するサーバ装置である。例えば、可視化ソリューションユニット５０４は、集計されたデータを地図上にマッピング処理することができる。この可視化ソリューションユニット５０４にて処理されたデータは、企業、官公庁または個人等に提供され、店舗開発、道路交通調査、災害対策、環境対策等に利用される。なお、このように統計処理された情報は、当然に個人情報は特定されないように加工されている。
【００２３】
ここで、人口算出装置６００の構成を説明する。図２は、人口算出装置６００の機能構成を示すブロック図である。図２に示すように、人口算出装置６００は、指示受付部６０１、長周期データ集計部６０２、短周期データ集計部６０３、共通データ集計部６０４、係数算出部６０５、人口算出部６０６、および結果出力部６０７を備えている。
【００２４】
指示受付部６０１は、ユーザからの人口算出指示、人口算出対象のエリアおよび時間帯を受け付ける部分である。また、指示受付部６０１は、外部装置から定期的に発信される人口算出指示、人口算出対象のエリアおよび時間帯を受け付けることもできる。
【００２５】
長周期データ集計部６０２は、指示受付部６０１において受け付けられた人口算出指示等に従って、ペタマイニングユニット５０２に蓄積された第１位置情報から、所定エリア内で第１位置情報を登録した移動機数（以下、長周期ユーザ数）を集計する部分である。このように、長周期データ集計部６０２は、移動機１００が第１周期で登録する第１位置情報に基づき、所定エリア内の移動機数を集計する第１集計手段を構成している。このように集計された第１位置情報は、図５に示される管理テーブルに記憶される。図５は、長周期データ集計部６０２が有する管理テーブルを示す説明図である。図５に示されるように、長周期データ集計部６０２の管理テーブルでは、第１位置情報が、識別情報、タイムスタンプ等と対応付けて記憶される。
【００２６】
短周期データ集計部６０３は、指示受付部６０１において受け付けられた人口算出指示等に従って、ペタマイニングユニット５０２に蓄積された第２位置情報から、所定エリア内で第２位置情報を登録した移動機数（以下、短周期ユーザ数）を集計する部分である。このように、短周期データ集計部６０３は、移動機１００が第２周期で登録する第２位置情報に基づき、所定エリア内の移動機数を集計する第２集計手段を構成している。このように集計された第２位置情報は、図６に示される管理テーブルに記憶される。図６は、短周期データ集計部６０３が有する管理テーブルを示す説明図である。図６に示されるように、短周期データ集計部６０３の管理テーブルでは、第２位置情報が、識別情報、タイムスタンプ等と対応付けて記憶される。
【００２７】
なお、第１位置情報及び第２位置情報が対応付けられる識別情報は、例えば、ユーザ識別子（例えば電話番号）に対して非識別化処理を行ったものである。非識別化処理の具体例としては、一方向性関数による不可逆符号への変換を含む処理が挙げられる。一方向性関数の具体例としては、国内外の評価プロジェクトや評価機関により推奨されているハッシュ関数に基づく鍵付ハッシュ関数が挙げられる。
【００２８】
共通データ集計部６０４は、所定エリア内で第１位置情報及び第２位置情報の両方を登録した移動機数（以下、共通ユーザ数）を集計する部分である。即ち、共通データ集計部６０４は、長周期データ集計部６０２で集計された移動機１００のうち、短周期データ集計部６０３で集計された移動機１００と識別情報の一致する移動機数を集計する。このように、共通データ集計部６０４は、第１集計手段で集計された移動機１００のうち、第２集計手段で集計された移動機１００と一致する移動機数を集計する共通部集計手段を構成している。
【００２９】
係数算出部６０５は、時間粒度詳細化係数を算出する部分である。時間粒度詳細化係数は、第１集計手段により集計された移動機数に対して第１周期よりも短く設定された期間における移動機数の増減分を調整する係数であり、下記式（１）により算出される。
【００３０】
時間粒度詳細化係数＝短周期ユーザ数／共通ユーザ数・・・（１）
【００３１】
人口算出部６０６は、係数算出部６０５で算出された時間粒度詳細化係数を用い、下記式（２）により算出されるユーザ数に基づき人口を算出する部分である。
【００３２】
ユーザ数＝長周期ユーザ数×時間粒度詳細化係数・・・（２）
【００３３】
なお、係数算出部６０５は、下記式（３）により、ＧＰＳ利用率、即ち上記第１位置情報に対する上記第２位置情報の利用率を算出する部分であってもよい。
【００３４】
ＧＰＳ利用率＝共通ユーザ数／長周期ユーザ数・・・（３）
【００３５】
この場合、人口算出部６０６は、下記式（４）によりユーザ数を算出する。
【００３６】
ユーザ数＝短周期ユーザ数／ＧＰＳ利用率・・・（４）
【００３７】
このように、係数算出部６０５および人口算出部６０６は、第１集計手段により集計された移動機数、第２集計手段により集計された移動機数、および共通部集計手段により集計された移動機数を用い所定エリア内のユーザ数を算出し、求めたユーザ数にシェア率の逆数を乗算するなどの統計処理を施すことにより、人口を算出する人口算出手段を構成している。
【００３８】
結果出力部６０７は、人口算出部６０６により算出された人口を出力する部分である。
【００３９】
続いて、人口算出装置６００の行う処理について説明する。図３は、人口算出装置６００の人口算出手順を示すフローチャートである。まず、指示受付部６０１により人口算出指示が検出され、人口算出処理が開始される。人口算出処理が開始されると、指示受付部６０１によって人口算出対象のエリアおよび時間帯が取得される（ステップＳ６００）。本実施形態では、人口算出対象のエリアが上記所定エリアとなる。また、人口算出対象のエリアは、上記セクタを一単位として取得されるものとし、人口算出対象エリアの境界線は上記セクタの境界線と一致しているものとする。なお、上記基地局のカバーエリアを一単位として人口算出対象のエリアを取得してもよい。
【００４０】
ステップＳ６００において人口算出対象のエリアおよび時間帯が取得されると、短周期データ集計部６０３により、人口算出対象の時間帯に人口算出対象のエリア内で第２位置情報を登録した移動機１００が抽出される（ステップＳ６０１）。即ち、人口算出対象の時間帯において、上記緯度および経度が人口算出対象のエリア内にある移動機１００が抽出される。そして、抽出された移動機１００の数が短周期ユーザ数として集計される（ステップＳ６０２）。
【００４１】
次に、長周期データ集計部６０２により、人口算出対象の時間帯に人口算出対象のエリア内で第１位置情報を登録した移動機１００が抽出される（ステップＳ６０３）。即ち、人口算出対象の時間帯において、人口算出対象のエリアに含まれるセクタのＩＤを登録した移動機１００が抽出される。そして、抽出された移動機１００の数が長周期ユーザ数として集計される（ステップＳ６０４）。
【００４２】
短周期ユーザ数および長周期ユーザ数が集計されると、共通データ集計部６０４により上記共通ユーザ数が集計される（ステップＳ６０５）。
【００４３】
次に、係数算出部６０５により、上記式（１）により時間粒度詳細化係数が算出され、または上記式（３）によりＧＰＳ利用率が算出される（ステップＳ６０６）。そして、人口算出部６０６により、上記式（２）または式（４）によりユーザ数が算出される（ステップＳ６０７）。
【００４４】
ここで、上記人口算出処理の具体例を示す。図４は、人口算出対象のエリアへの移動機の出入りを示す模式図である。
【００４５】
図４に示すように、人口算出対象のエリアには識別情報ａ，ｂ，ｃ〜ｊの合計１０台の移動機１００が出入している。例えば識別情報ｂの移動機１００は、実線の引かれた期間に人口算出対象のエリアに在圏し、ポイント６１において第１位置情報の登録処理を行い、ポイント６２において第２位置情報の登録処理を行っている。
【００４６】
この状況において、人口算出対象の時間帯Ｔの短周期ユーザとして集計される移動機１００は、識別情報ｂ，ｃ，ｆ，ｇ，ｊの５台の移動機となる（図６参照）。また、長周期ユーザとして集計される移動機１００は識別情報ａ，ｂの２台の移動機となる（図５参照）。そして、共通ユーザとして集計される移動機１００は、長周期ユーザとして集計される識別情報ａ，ｂの移動機うち、短周期ユーザとしても集計される識別情報ｂの１台の移動機となる。
【００４７】
以上のように、短周期ユーザ数が５人、共通ユーザ数が１人であるため、式（１）より時間粒度詳細化係数が５となる。そして、長周期ユーザ数が２人であるため、式（２）よりユーザ数の算出結果は１０人となる。
【００４８】
また、長周期ユーザ数が２人、共通ユーザ数が１人であるため、式（３）よりＧＰＳ利用率が０．５となる。そして、短周期ユーザ数が５人であるため、式（４）よりユーザ数の算出結果は１０人となる。
【００４９】
このようにユーザ数が算出されると、そのユーザ数に拡大係数が乗算され、人口が算出される。拡大係数は、ユーザ数を人口に変換するための係数であり、例えばユーザの属性ごとに予め求められ、記憶されている。拡大係数は、時間帯やエリアに応じて求められていてもよい。
【００５０】
また、予め記憶しておく以外に、以下のようにして拡大係数を求めてもよい。例えば、「在圏率」と「端末の普及率」との積（即ち、人口に対する在圏数の比率）の逆数を拡大係数とすることができる。ここで、「在圏率」とは、契約台数に対する在圏数の比率を意味し、「端末の普及率」とは人口に対する端末の契約台数の比率を意味する。人口に対する在圏数の比率は、例えば、ある集計単位ごとに移動機数を集計することで得られるユーザ数ピラミッドデータと、予め求められた統計データ（例えば住民基本台帳など）から同じ集計単位ごとに取得される人口ピラミッドデータと、を用い、集計単位ごとに求めることができる。集計単位としては、住所の都道府県ごと、５才刻み年齢層ごと、男女ごと、１時間の時間帯ごと、等が挙げられる。なお、拡大係数は、集計単位ごとに求めることが望ましいが、必須ではない。
【００５１】
最後に、結果出力部６０７により算出結果が出力される（ステップＳ６０８）。以上により人口算出処理が終了する。なお、以上の人口算出処理を、性別・年齢層・住所といった属性ごとに行い、属性ごとに人口を出力できるようにしてもよい。
【００５２】
以上説明したように、本発明に係る人口算出装置６００によれば、５４分周期で登録された第１位置情報に基づき長周期ユーザ数が集計され、５分周期で登録された第２位置情報に基づき短周期ユーザ数が集計され、さらに共通ユーザ数が集計される。そして、長周期ユーザ数に対して５４分周期よりも短い５分周期における移動機数の増減分を調整する時間粒度詳細化係数が算出され、長周期ユーザ数に時間粒度詳細化係数を乗算することで、５分間での増減分を含むユーザ数が算出され、当該ユーザ数に基づいて人口が算出される。
【００５３】
なお、時間粒度詳細化係数に代えて、ＧＰＳ利用率即ち第１位置情報に対する第２位置情報の利用率を算出し、短周期ユーザ数からＧＰＳ利用率を除算することで、短周期ユーザ数をＧＰＳ利用率が１００％である場合に換算し、ユーザ数を算出することもできる。
【００５４】
このように、人口算出装置６００によれば、５４分周期で登録される第１位置情報と、５分周期で登録される第２位置情報とが互いに補われる。これにより、第１位置情報のみでは急激な人口の増減をとらえられない一方で、第２位置情報を登録する移動機数が少なく第２位置情報のみでは人口を算出できない状況であっても、第１位置情報と第２位置情報とが互いに補われ、急激な人口の増減を把握することが可能となる。
【００５５】
続いて、本発明に係る人口算出装置の第２実施形態について説明する。第２実施形態は、セクタを一単位としていないエリアの人口を算出することができる点で第１実施形態と異なる。第２実施形態の機能構成は、第１実施形態の機能構成と同様であるため、第２実施形態については人口算出処理手順のみを説明する。
【００５６】
図７に示すように、まず、第１実施形態と同様に、指示受付部６０１によって人口算出対象のエリア・時間帯が取得される（ステップＳ７００）。さらに、人口算出対象のエリアと一部が重複するセクタが関連セクタとして抽出される（ステップＳ７０１）。本実施形態では、抽出された関連セクタの全てを含むエリアが上記所定エリアとなる。図８は、関連セクタを例示する図である。図８では、人口算出対象のエリアＡ１が、例えば周波数２ＧＨｚの基地局のセクタＢ１及びＢ２の一部と重複すると共に、周波数８００ＭＨｚの基地局のセクタＣ１の一部と重複している場合を示している。この場合、セクタＢ１，Ｂ２及びセクタＣ１が関連セクタとして抽出され、これらのセクタの全てを含むエリアＤが上記所定エリアとなる。
【００５７】
次に、短周期データ集計部６０３により、人口算出対象の時間帯に関連セクタ内で第２位置情報を登録した移動機１００が抽出される（ステップＳ７０２）。即ち、人口算出対象の時間帯において、緯度および経度が関連セクタ内にある移動機１００が抽出される。そして、抽出された移動機１００の数が短周期ユーザ数として集計される（ステップＳ７０３）。
【００５８】
次に、長周期データ集計部６０２により、人口算出対象の時間帯に関連セクタ内で第１位置情報を登録した移動機１００が抽出される（ステップＳ７０４）。即ち、人口算出対象の時間帯において、関連セクタのいずれかのＩＤを登録した移動機１００が抽出される。そして、抽出された移動機１００の数が長周期ユーザ数として集計される（ステップＳ７０５）。
【００５９】
次に、共通データ集計部６０４により、関連セクタ内の共通ユーザ数が集計される（ステップＳ７０６）。
【００６０】
次に、係数算出部６０５により、上記式（１）によって関連セクタ全体での時間粒度詳細化係数が算出され、または上記式（３）により関連セクタ全体でのＧＰＳ利用率が算出される（ステップＳ７０７）。
【００６１】
次に、人口算出部６０６により、人口算出対象のエリア内のユーザ数が算出される（ステップＳ７０８）。ここでは、まず上記式（２）または式（４）によって関連セクタごとのユーザ数が算出される。この際に、各関連セクタの時間粒度詳細化係数またはＧＰＳ利用率には、関連セクタ全体での時間粒度詳細化係数またはＧＰＳ利用率が共通して用いられる。更に、各関連セクタと人口算出対象のエリアとの重複領域の面積が、その関連セクタの面積にしめる割合に応じて、各関連セクタのユーザ数が人口算出対象のエリアに按分される。即ち、各関連セクタのユーザ数に上記割合が乗算され、乗算結果が人口算出対象のエリアに割り当てられる。そして、各関連セクタから人口算出対象のエリアに割り当てられたユーザ数が合計され、人口算出対象のエリアのユーザ数が算出される。
【００６２】
例えば、図８のエリアＡ１のユーザ数は次のように算出される。図８では、関連セクタＢ１のユーザ数が１０００人、関連セクタＣ１のユーザ数が５００人、関連セクタＢ２のユーザ数が２０００人であるとする。また、関連セクタＢ１とエリアＡ１との重複領域の面積が関連セクタＢ１の面積にしめる割合が０．２、関連セクタＢ２とエリアＡ１との重複領域の面積が関連セクタＢ２の面積にしめる割合が０．２、関連セクタＣ１とエリアＡ１との重複領域の面積が関連セクタＣ１の面積にしめる割合が０．５であるとする。この場合、
エリアＡ１のユーザ数＝１０００人×０．２＋５００人×０．５＋２０００人×０．２＝８５０人
となる。
【００６３】
次に、第１実施形態と同様に拡大係数を乗算して人口が算出され（ステップＳ７０９）、人口算出結果が出力される（ステップＳ７１０）。以上により第２実施形態の人口算出処理が終了する。
【００６４】
本実施形態によれば、セクタを一単位としていないエリアの人口を算出することができ、例えば市区町村等、様々なエリアの人口を算出することができる。
【００６５】
なお、第２実施形態では、関連セクタ全体で時間粒度詳細化係数またはＧＰＳ利用率を算出している。即ち、第２実施形態では、時間粒度詳細化係数等を関連セクタごとには算出していない。これにより、第２位置情報のみを登録した移動機がいずれのセクタに属しているかを識別することなく、時間粒度詳細化係数等を算出することができる。第２位置情報のみを登録した移動機についても、いずれのセクタに属しているかを識別可能な場合には、長周期ユーザ数、短周期ユーザ数、及び共通ユーザ数を関連セクタごとに集計し、時間粒度詳細化係数等を関連セクタごとに算出してもよい。例えば、第２位置情報を登録した移動機が、登録処理の際にセクタＩＤを併せて登録している場合には、第２位置情報のみを登録した移動機についてもいずれのセクタに属するかを識別できる。
【００６６】
続いて、本発明に係る人口算出装置の第３実施形態について説明する。第３実施形態は、セクタを一単位としていないエリアの人口を算出する処理を、第２実施形態とは異なる手順で行うものである。
【００６７】
図９に示すように、まず、第２実施形態と同様に、人口算出対象のエリア・時間帯が取得され（ステップＳ８００）、関連セクタが抽出される（ステップＳ８０１）。本実施形態では、人口算出対象のエリアが上記所定エリアとなる。
【００６８】
次に、短周期データ集計部６０３により、人口算出対象の時間帯に人口算出対象のエリア内で第２位置情報を登録した移動機１００が抽出される（ステップＳ８０２）。即ち、人口算出対象の時間帯において、上記緯度および経度が人口算出対象のエリア内にある移動機１００が抽出される。そして、抽出された移動機１００の数が短周期ユーザ数として集計される（ステップＳ８０３）。
【００６９】
次に、長周期データ集計部６０２により、人口算出対象の時間帯に関連セクタ内で第１位置情報を登録した移動機１００が抽出される（ステップＳ８０４）。即ち、人口算出対象の時間帯において、関連セクタのいずれかのＩＤを登録した移動機１００が抽出される。そして、抽出された移動機１００の数が、各関連セクタごとに長周期ユーザ数として集計される（ステップＳ８０５）。
【００７０】
次に、各関連セクタと人口算出対象のエリアとの重複領域の面積が、その関連セクタの面積にしめる割合に応じて、各関連セクタの長周期ユーザ数が人口算出対象のエリアに按分される。即ち、各関連セクタの長周期ユーザ数に上記割合が乗算され、乗算結果が人口算出対象のエリアに割り当てられる。そして、各関連セクタから人口算出対象のエリアに割り当てられた長周期ユーザ数が合計され、人口算出対象のエリアの長周期ユーザ数が集計される（ステップＳ８０６）。
【００７１】
次に、共通データ集計部６０４により、人口算出対象のエリア内の共通ユーザ数が集計される（ステップＳ８０７）。
【００７２】
次に、係数算出部６０５により、上記式（１）によって人口算出対象のエリア内での時間粒度詳細化係数が算出され、または上記式（３）により人口算出対象のエリア内でのＧＰＳ利用率が算出される（ステップＳ８０８）。
【００７３】
次に、人口算出部６０６により、上記式（２）または式（４）によって人口算出対象のエリア内のユーザ数が算出される（ステップＳ８０９）。次に、第１実施形態と同様に拡大係数を乗算して人口が算出され（ステップＳ８１０）、人口算出結果が出力される（ステップＳ８１１）。以上により第３実施形態の人口算出処理が終了する。本実施形態も、第２実施形態と同様の効果を奏する。
【００７４】
以上、本発明に係る好適な実施形態について説明してきたが、本発明は必ずしも上述した実施形態に限定されるものではなく、その要旨を変更しない範囲で様々な変更が可能である。例えば、上述した実施形態では、短周期ユーザ数が長周期ユーザ数よりも先に集計されているが、長周期ユーザ数を短周期ユーザ数よりも先に算出してもよい。
【００７５】
また、上述した実施形態では、人口算出装置はペタマイニングユニットに設けられているものとしたがこれに限られない。例えば、人口算出装置はモバイルデモグラフィユニットに設けられていてもよい。また、例えば、人口算出装置のうち結果出力部のみがモバイルデモグラフィユニットに設けられる等、人口算出装置がペタマイニングユニットやモバイルデモグラフィユニットに分散して設けられていてもよい。
【００７６】
また、上述した実施形態では、人口算出対象の時間帯に、所定エリア内で第１位置情報及び第２位置情報の少なくとも一方を登録した移動機が抽出され、長周期ユーザ数、短周期ユーザ数、及び共通ユーザ数が集計されているがこれに限られない。例えば、次に説明する端末数推計により、長周期ユーザ数、短周期ユーザ数、及び共通ユーザ数のそれぞれを集計してもよい。以下に、端末数推計の考え方及び計算方法を説明する。
【００７７】
図１０に示すモデルのように、人口算出対象の時間帯Ｔの間に、ｎ個の端末ａ１，ａ２，…，ａｎがエリアＳを通過し、各端末ａｉの推計対象期間内のエリアＳの滞在時間がｔｉ（０＜ｔｉ≦Ｔ）であったとする。端末数推計によって算出しようとする端末数ｍは、以下の式（５）で表わされる。
【数１】
0005735888.tif 000002
即ち、端末数推計は、各端末ａｉの推計対象期間内のエリアＳの滞在時間ｔｉの総和を推計対象期間の長さＴで除した結果を、端末数ｍとして推計しようとするものである。端末ａｉの滞在時間ｔｉの真の値は観測不能であるが、各端末ａｉが第１位置情報や第２位置情報を登録するために発する信号は観測可能である。
【００７８】
端末ａｉが人口算出対象の時間帯ＴにエリアＳで発信した信号を、時刻順に
【数２】
0005735888.tif 000003
（ｘｉは、端末ａｉが推計対象期間内にエリアＳで発信した信号の総数）とすると、観測された信号ｑｉｊ（ｊは１以上ｘｉ以下の整数）からｍの値を推計することができる。
【００７９】
図１１に基づき端末数推計の計算方法を説明する。端末ａｉから信号ｑｉｊが送信される密度（即ち、単位時間あたりの信号数）をｐｉとする。このとき、信号が送信される確率がエリアに対して独立であれば、端末ａｉが推計対象期間内にエリアＳで発信した信号の総数ｘｉの期待値Ｅ（ｘｉ）は、Ｅ（ｘｉ）＝ｔｉ×ｐｉであるため、端末ａｉの推計対象期間内のエリアＳの滞在時間ｔｉの期待値Ｅ（ｔｉ）について以下の式（６）が成立する。
Ｅ（ｔｉ）＝ｘｉ／ｐｉ・・・（６）
信号ｑｉｊの送信時刻をｕｉｊとしたとき、信号ｑｉｊの密度ｐｉｊは、以下の式（７）で与えられる。
ｐｉｊ＝２／（ｕｉ（ｊ＋１）−ｕｉ（ｊ−１））・・・（７）
ここで、上記式（７）の（ｕｉ（ｊ＋１）−ｕｉ（ｊ−１））を、信号ｑｉｊについての特徴量ｗｉｊとする。すると、上記式（７）は、以下となる。
ｐｉｊ＝２／ｗｉｊ・・・（８）
【００８０】
このとき密度ｐｉは、
【数３】
0005735888.tif 000004
で与えられるため、端末数ｍの推計値Ｅ（ｍ）は以下の式（１０）で計算することができる。
【数４】
0005735888.tif 000005
【００８１】
図１１の例に示すように、集計時間帯内であり且つ移動機ａｉがエリアＳに滞在していた期間内に、移動機ａｉはＧＰＳ位置情報ｑｉ１、ｑｉ２、ｑｉ３を送信し、ＧＰＳ位置情報ｑｉ１の直前にＧＰＳ位置情報ｑｉ０を、ＧＰＳ位置情報ｑｉ３の直後にＧＰＳ位置情報ｑｉ４を送信したものとし、ＧＰＳ位置情報ｑｉ０、ｑｉ１、ｑｉ２、ｑｉ３、ｑｉ４の送信時刻をそれぞれｕｉ０、ｕｉ１、ｕｉ２、ｕｉ３、ｕｉ４とすると、上記の考え方は、移動機ａｉの集計時間帯内のエリアＳの滞在時間ｔｉを、（ｕｉ０とｕｉ１の中点）から（ｕｉ３とｕｉ４の中点）までの期間と推計することに相当する。
【００８２】
このような端末数推計を用いる場合、例えば、長周期データ集計部６０２は、人口算出対象の時間帯Ｔに所定エリア内で登録された第１位置情報のそれぞれについて、第１特徴量を算出する。例えば、移動機が第１位置情報を登録した時刻、同じ移動機が一つ前の第１位値情報を登録した時刻（直前登録時刻）、及び同じ移動機が一つ後の第１位置情報を登録した時刻（直後登録時刻）、を抽出し、
第１特徴量＝直後登録時刻−直前登録時刻
として第１特徴量を算出する。
そして、
長周期ユーザ数＝第１特徴量の総和／２／第１基準時間
として長周期ユーザ数を集計する。第１基準時間の長さは、人口算出の条件に応じて適宜設定できるが、少なくとも第１周期以上に設定することが好ましい。
【００８３】
また、短周期データ集計部６０３は、人口算出対象の時間帯Ｔに所定エリア内で登録された第２位置情報のそれぞれについて、第２特徴量を算出する。例えば、移動機が第２位置情報を登録した時刻、同じ移動機が一つ前の第２位値情報を登録した時刻（直前登録時刻）、及び同じ移動機が一つ後の第２位置情報を登録した時刻（直後登録時刻）、を抽出し、
第２特徴量＝直後登録時刻−直前登録時刻
として第２特徴量を算出する。
そして、
短周期ユーザ数＝第２特徴量の総和／２／第２基準時間
として短周期ユーザ数を集計する。第２基準時間の長さは、人口算出の条件に応じて適宜設定できるが、少なくとも第２周期以上に設定することが好ましい。
【００８４】
また、共通データ集計部６０４は、人口算出対象の時間帯Ｔに所定エリア内で第１位置情報を登録した移動機によって、人口算出対象の時間帯Ｔに所定エリア内で登録された第２位置情報のそれぞれについて、共通特徴量を算出する。例えば、移動機が第２位置情報を登録した時刻、同じ移動機が一つ前の第２位値情報を登録した時刻（直前登録時刻）、及び同じ移動機が一つ後の第２位置情報を登録した時刻（直後登録時刻）、を抽出し、
共通特徴量＝直後登録時刻−直前登録時刻
として共通特徴量を算出する。
そして、
共通ユーザ数＝共通特徴量の総和／２／共通基準時間
として共通ユーザ数を集計する。共通基準時間の長さは、人口算出の条件に応じて適宜設定できるが、少なくとも第２周期以上に設定することが好ましい。このように、共通ユーザ数の集計に、第１位置情報に比べ信号の発信周期が短い第２位置情報を用いることで、共通ユーザ数の精度を高めることができる。
【００８５】
なお、このような端末数推計を、上記拡大係数の算出に適用することもできる。具体的には、端末数推計を用いて、上記集計単位ごとに移動機数を集計し、上記ユーザ数ピラミッドデータを得ることができる。
【符号の説明】
【００８６】
１０…通信システム、１００…移動機、６００…人口算出装置、６０１…指示受付部、６０２…長周期データ集計部、６０３…短周期データ集計部、６０４…共通データ集計部、６０５…係数算出部、６０６…人口算出部、６０７…結果出力部。
【図１】
0005735888.tif 000006
【図２】
0005735888.tif 000007
【図３】
0005735888.tif 000008
【図４】
0005735888.tif 000009
【図５】
0005735888.tif 000010
【図６】
0005735888.tif 000011
【図７】
0005735888.tif 000012
【図８】
0005735888.tif 000013
【図９】
0005735888.tif 000014
【図１０】
0005735888.tif 000015
【図１１】
0005735888.tif 000016