(57)【特許請求の範囲】
【請求項１】
位置登録が行われたセクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を外部から受信する受信部と、
各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出部と、
各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出部と、
を備え、
前記人口分布算出部は、
各ユーザについて、直前の位置登録時刻と対象時刻との第１の時間差、および、直後の位置登録時刻と対象時刻との第２の時間差を算出し、
各ユーザについて、第１の時間差と第２の時間差に基づいて、直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトを算出し、
集計対象全ユーザについての直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトをセクタごとに集計する、
ことにより、対象時刻におけるセクタ単位の人口分布を算出する、
ことを特徴とする情報分析装置。
【請求項２】
前記人口分布算出部は、
第１の時間差と第２の時間差の和により第２の時間差を除算して得られた値を、直前の位置登録信号のウエイトとして算出し、
第１の時間差と第２の時間差の和により第１の時間差を除算して得られた値を、直後の位置登録信号のウエイトとして算出する、
ことを特徴とする請求項１に記載の情報分析装置。
【請求項３】
前記人口分布算出部は、
前記第１の時間差および前記第２の時間差のうち少なくとも一方が、予め定められた基準値より長い場合、当該基準値より長い時間差に対応する位置登録信号のウエイトを０とする、
ことを特徴とする請求項１に記載の情報分析装置。
【請求項４】
前記人口分布算出部は、
前記第１の時間差および前記第２の時間差のうち一方のみが、前記基準値より長い場合、前記基準値以下の時間差に対応する位置登録信号のウエイトを１とする、
ことを特徴とする請求項３に記載の情報分析装置。
【請求項５】
前記人口分布算出部は、
前記第１の時間差および前記第２の時間差のうち一方のみが、前記基準値より長い場合、前記基準値と当該基準値以下の時間差との差分を当該基準値により除算して得られた値を、当該基準値以下の時間差に対応する位置登録信号のウエイトとして算出する、
ことを特徴とする請求項３に記載の情報分析装置。
【請求項６】
位置登録が行われたセクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を外部から受信する受信部と、
各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出部と、
各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出部と、
を備え、
前記抽出部は、
各ユーザについての位置登録信号のうち、端末が位置登録エリアを跨いだことに起因して生成された位置登録信号を除去する前の位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号である直前の除去前位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号である直後の除去前位置登録信号とを抽出するとともに、各ユーザについての位置登録信号のうち、端末が位置登録エリアを跨いだことに起因して生成された位置登録信号を除去した後の位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号である直前の除去後位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号である直後の除去後位置登録信号とを抽出し、
前記人口分布算出部は、
各ユーザについて、直前の除去前位置登録信号の位置登録時刻と対象時刻との第１の時間差、直後の除去前位置登録信号の位置登録時刻と対象時刻との第２の時間差、直前の除去後位置登録信号の位置登録時刻と対象時刻との第３の時間差、および、直後の除去後位置登録信号の位置登録時刻と対象時刻との第４の時間差を算出し、
各ユーザについて、第１の時間差と第２の時間差に基づいて、直前の除去前位置登録信号のウエイトおよび直後の除去前位置登録信号のウエイトを算出するとともに、第３の時間差と第４の時間差に基づいて、直前の除去後位置登録信号のウエイトおよび直後の除去後位置登録信号のウエイトを算出し、
予め定められた按分パラメータに基づいて、各ユーザについて、直前および直後の除去前位置登録信号のウエイトと直前および直後の除去後位置登録信号のウエイトとを按分し、
集計対象全ユーザについての、按分後の直前および直後の除去前位置登録信号のウエイト、並びに按分後の直前および直後の除去後位置登録信号のウエイトを、セクタごとに集計する、
ことを特徴とする情報分析装置。
【請求項７】
位置登録が行われたセクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を外部から受信する受信部と、
各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出部と、
各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出部と、
を備え、
前記人口分布算出部は、
各ユーザについて、対象時刻の直前の位置登録信号に含まれるセクタＩＤに対応する位置と直後の位置登録信号に含まれるセクタＩＤに対応する位置とを結ぶ線上に、直前の位置登録時刻から直後の位置登録時刻まで等速移動した場合の所定時間おきの通過点に相当する複数の点を生成し、生成された点が含まれるセクタに対し仮想の位置登録信号を生成し、
集計対象全ユーザについての仮想の位置登録信号から対象時刻の位置登録信号を抽出し、
抽出された対象時刻の位置登録信号数をセクタごとに集計する、
ことにより、対象時刻におけるセクタ単位の人口分布を算出する、
ことを特徴とする情報分析装置。
【請求項８】
位置登録が行われたセクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を外部から受信する受信部と、
各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出部と、
各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出部と、
を備え、
前記人口分布算出部は、
各ユーザについて、対象時刻の直前の位置登録信号に含まれるセクタＩＤに対応する位置と直後の位置登録信号に含まれるセクタＩＤに対応する位置とを結ぶ線上に、直前の位置登録時刻から直後の位置登録時刻まで等速移動した場合の所定時間おきの通過点に相当する複数の点を生成し、
対象時刻に対応する点が含まれるセクタに対し仮想の位置登録信号を生成し、
生成された対象時刻の位置登録信号数をセクタごとに集計する、
ことにより、対象時刻におけるセクタ単位の人口分布を算出する、
ことを特徴とする情報分析装置。
【請求項９】
位置登録が行われたセクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を外部から受信する受信部と、
各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出部と、
各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出部と、
を備え、
前記人口分布算出部は、
境界がセクタ境界と対応付け可能な領域における領域内移動速度の領域間相対比率に相関する比率に基づいて、各ユーザについて、対象時刻の直前の位置登録信号に含まれるセクタＩＤに対応する位置と直後の位置登録信号に含まれるセクタＩＤに対応する位置とを結ぶ線上に、直前の位置登録時刻から直後の位置登録時刻まで移動した場合の所定時間おきの通過点に相当する複数の点を生成し、生成された点が含まれるセクタに対し仮想の位置登録信号を生成し、
集計対象全ユーザについての仮想の位置登録信号から対象時刻の位置登録信号を抽出し、
抽出された対象時刻の位置登録信号数をセクタごとに集計する、
ことにより、対象時刻におけるセクタ単位の人口分布を算出する、
ことを特徴とする情報分析装置。
【請求項１０】
位置登録が行われたセクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を外部から受信する受信部と、
各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出部と、
各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出部と、
を備え、
前記人口分布算出部は、
境界がセクタ境界と対応付け可能な領域における領域内移動速度の領域間相対比率に相関する比率に基づいて、各ユーザについて、対象時刻の直前の位置登録信号に含まれるセクタＩＤに対応する位置と直後の位置登録信号に含まれるセクタＩＤに対応する位置とを結ぶ線上に、直前の位置登録時刻から直後の位置登録時刻まで移動した場合の所定時間おきの通過点に相当する複数の点を生成し、
対象時刻に対応する点が含まれるセクタに対し仮想の位置登録信号を生成し、
生成された対象時刻の位置登録信号数をセクタごとに集計する、
ことにより、対象時刻におけるセクタ単位の人口分布を算出する、
ことを特徴とする情報分析装置。
【請求項１１】
前記情報分析装置は、ユーザごとの属性情報を保持した属性情報保持部、をさらに備え、
前記抽出部は、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前および直後の位置登録信号の抽出を行う、
ことを特徴とする請求項１〜１０の何れか１項に記載の情報分析装置。
【請求項１２】
前記情報分析装置は、ユーザごとの属性情報を保持した属性情報保持部、をさらに備え、
前記人口分布算出部は、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として対象時刻におけるセクタ単位の人口分布の算出を行う、
ことを特徴とする請求項１〜１０の何れか１項に記載の情報分析装置。
【請求項１３】
ユーザの位置を示す位置情報、当該位置情報が得られた測位時刻情報およびユーザＩＤを含むポイントデータを外部から受信する受信部と、
各ユーザについてのポイントデータから、測位時刻が対象時刻の直前のポイントデータと、測位時刻が対象時刻の直後のポイントデータとを抽出する抽出部と、
各ユーザについて、対象時刻の直前のポイントデータの示す位置と対象時刻の直後のポイントデータの示す位置との間を補完することにより、対象時刻におけるユーザの位置を推定し、推定した各ユーザの位置に基づいて対象時刻における所定の算出対象エリア単位の人口分布を算出する人口分布算出部と、
を備え、
前記人口分布算出部は、
各ユーザについて、対象時刻の直前のポイントデータに含まれる位置情報の示す位置と直後のポイントデータに含まれる位置情報の示す位置とを結ぶ線上に、直前の測位時刻から直後の測位時刻まで移動した場合の所定時間おきの通過点に相当する複数の点を生成して、対象時刻に対応する点を対象時刻におけるユーザの位置として推定し、推定した各ユーザの位置に基づいて対象時刻における所定の算出対象エリア単位の人口分布を算出する、
ことを特徴とする情報分析装置。
【請求項１４】
前記人口分布算出部は、
前記所定時間おきの通過点に相当する複数の点を生成するにあたって、前記複数の点を生成しようとする生成エリアにおけるユーザの移動特性に基づいて、該生成エリア内に生成する点の間隔を決定する、
ことを特徴とする請求項１３に記載の情報分析装置。
【請求項１５】
前記生成エリアにおけるユーザの移動特性は、該生成エリアにおけるユーザの平均移動速度である、
ことを特徴とする請求項１４に記載の情報分析装置。
【請求項１６】
前記生成エリアにおけるユーザの移動特性は、該生成エリアにおけるユーザの平均滞留時間である、
ことを特徴とする請求項１４に記載の情報分析装置。
【請求項１７】
前記情報分析装置は、ユーザごとの属性情報を保持した属性情報保持部、をさらに備え、
前記抽出部は、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前および直後のポイントデータを抽出する、
ことを特徴とする請求項１３〜１６の何れか１項に記載の情報分析装置。
【請求項１８】
前記情報分析装置は、ユーザごとの属性情報を保持した属性情報保持部、をさらに備え、
前記人口分布算出部は、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として対象時刻における所定の算出対象エリア単位の人口分布を算出する、
ことを特徴とする請求項１３〜１６の何れか１項に記載の情報分析装置。
【請求項１９】
前記人口分布算出部は、さらに、ユーザの端末数を人口に変換するための拡大係数をも基礎として、前記人口分布を算出する、
ことを特徴とする請求項１〜１８の何れか１項に記載の情報分析装置。
【請求項２０】
対象時刻を入力するための対象時刻入力部、
をさらに備える請求項１〜１９の何れか１項に記載の情報分析装置。
【請求項２１】
算出された前記人口分布情報を出力する出力部、
をさらに備える請求項１〜２０の何れか１項に記載の情報分析装置。
【請求項２２】
情報分析装置により実行される情報分析方法であって、
位置登録が行われたセクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を外部から受信する受信ステップと、
各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出ステップと、
各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出ステップと、
を備え、
前記人口分布算出ステップにて、前記情報分析装置は、
各ユーザについて、直前の位置登録時刻と対象時刻との第１の時間差、および、直後の位置登録時刻と対象時刻との第２の時間差を算出し、
各ユーザについて、第１の時間差と第２の時間差に基づいて、直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトを算出し、
集計対象全ユーザについての直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトをセクタごとに集計する、
ことにより、対象時刻におけるセクタ単位の人口分布を算出する、
ことを特徴とする情報分析方法。
【請求項２３】
情報分析装置により実行される情報分析方法であって、
ユーザの位置を示す位置情報、当該位置情報が得られた測位時刻情報およびユーザＩＤを含むポイントデータを外部から受信する受信ステップと、
各ユーザについてのポイントデータから、測位時刻が対象時刻の直前のポイントデータと、測位時刻が対象時刻の直後のポイントデータとを抽出する抽出ステップと、
各ユーザについて、対象時刻の直前のポイントデータの示す位置と対象時刻の直後のポイントデータの示す位置との間を補完することにより、対象時刻におけるユーザの位置を推定し、推定した各ユーザの位置に基づいて対象時刻における所定の算出対象エリア単位の人口分布を算出する人口分布算出ステップと、
を備え、
前記人口分布算出ステップにて、前記情報分析装置は、
各ユーザについて、対象時刻の直前のポイントデータに含まれる位置情報の示す位置と直後のポイントデータに含まれる位置情報の示す位置とを結ぶ線上に、直前の測位時刻から直後の測位時刻まで移動した場合の所定時間おきの通過点に相当する複数の点を生成して、対象時刻に対応する点を対象時刻におけるユーザの位置として推定し、推定した各ユーザの位置に基づいて対象時刻における所定の算出対象エリア単位の人口分布を算出する、
ことを特徴とする情報分析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、人口分布を求める情報分析装置および情報分析方法に関する。
【背景技術】
【０００２】
従来、巨視的な人口分布に関するデータを収集する手法として、５年周期で全国的に実施される国勢調査があった。この調査は、対象者に対しアンケートを配布して回収し、多くの人手を使って集計するなどの非常に手間のかかる一連の作業が必要であり、調査結果を得るまでに時間がかかっていた。また、回答する対象者にとってもアンケートに記入して返送するといった負担を強いるものであった。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００３−４４９６９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
上記のように、従来は巨視的な人口分布に関するデータを収集するだけでも非常に面倒であり、かかるデータを簡易に且つ速やかに収集することは困難であった。【０００５】
一方、携帯端末を利用して人口分布を得る試みとして例えば特許文献１には、ＧＰＳ機能付き携帯端末を利用して人口の分布を求める点について記載されているが、この技術でも、対象とするユーザの全員に対してユーザ全員を一意に識別可能な識別子が付加されたＧＰＳ内蔵機器を配布する、及びＧＰＳを用いて測位した情報を一人一人収集する等の、人口分布を求めるための多大な処理負荷と時間を要する。そのため、より簡易且つ速やかに人口分布に関するデータを収集する技術が待望されていた。また、収集されるデータには、一定レベル以上の精度が求められるのは言うまでもない。なお、上記のようなＧＰＳ機能を用いた場合でも、人口分布を求めたい対象時刻きっかりの測位情報が必ずしも得られるとは限らない。【０００６】
本発明は、上記課題に鑑み、巨視的な人口分布に関する精度の良いデータを簡易に且つ速やかに収集することを目的とする。
【課題を解決するための手段】
【０００７】
出願人は、上記の目的を達成するために、位置登録信号を用いた人口分布の導出に係る発明について、日本出願（特願２００９−９２２２５号）および国際特許出願（ＰＣＴ／ＪＰ２０１０／０５５４２４）を既に提出した。その後も、出願人は、位置登録信号を用いた人口分布の導出に係る発明のさらなる改良を進め、今回、同発明を新たな視点から改良した本件発明をするに至った。【０００８】
本件発明の一側面に係る情報分析装置は、位置登録が行われたセクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を外部から受信する受信部と、各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出部と、各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出部と、を備え、人口分布算出部は、各ユーザについて、直前の位置登録時刻と対象時刻との第１の時間差、および、直後の位置登録時刻と対象時刻との第２の時間差を算出し、各ユーザについて、第１の時間差と第２の時間差に基づいて、直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトを算出し、集計対象全ユーザについての直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトをセクタごとに集計することにより、対象時刻におけるセクタ単位の人口分布を算出することを特徴とする。なお、本件において「対象時刻」とは、人口分布を求めるべき時刻を意味する。
【０００９】
上記情報分析装置では、受信部が、セクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を外部から受信すると、抽出部は、各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出し、そして、人口分布算出部は、各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する。具体的には、人口分布算出部は、各ユーザについて、直前の位置登録時刻と対象時刻との第１の時間差、および、直後の位置登録時刻と対象時刻との第２の時間差を算出し、各ユーザについて、第１の時間差と第２の時間差に基づいて、直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトを算出し、集計対象全ユーザについての直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトをセクタごとに集計することにより、対象時刻におけるセクタ単位の人口分布を算出する。なお、対象とするユーザは、位置登録信号を受信した全ユーザでもよいし、その中から属性情報（例えば年齢、性別、住所など）に応じて限定したユーザでもよい。
【００１０】
以上のように、各ユーザについての対象時刻の直前の位置登録信号と直後の位置登録信号を基礎として、対象時刻におけるセクタ単位の人口分布を算出することにより、巨視的な人口分布に関する精度の良いデータを簡易に且つ速やかに収集することが可能となる。また、このとき、位置登録信号の発生数を基礎とすることなく、各ユーザについての対象時刻の直前の位置登録信号および直後の位置登録信号を基礎として、対象時刻におけるセクタ単位の人口分布を算出するため、生成される時間間隔が変動する位置登録信号から、対象時刻におけるセクタ単位の人口分布を精度良く算出することができる。また、位置登録信号の発生数を基礎とした場合に想定される不都合（例えば、位置登録エリア境界付近で実際の人口数よりもかなり多くの位置登録信号が発生することによる不都合）を防止することができる。【００１１】
なお、本件において「位置登録信号のウエイト」とは、あるユーザが対象時刻において当該位置登録信号に係るセクタに在圏している確率を意味する。実際に、対象時刻ちょうどの位置登録信号があるとは限らないため、位置登録信号の位置登録時刻と対象時刻との時間差が短いほど、ユーザが当該位置登録信号に係るセクタに在圏している確率（即ち、当該位置登録信号のウエイト）が大きくなるように、位置登録信号のウエイトが算出される。例えば、人口分布算出部は、第１の時間差と第２の時間差の和により第２の時間差を除算して得られた値を、直前の位置登録信号のウエイトとして算出し、第１の時間差と第２の時間差の和により第１の時間差を除算して得られた値を、直後の位置登録信号のウエイトとして算出してもよい。例えば、あるユーザについて、直前の位置登録信号のウエイトとして0.4、直後の位置登録信号のウエイトとして0.6が算出された場合、当該ユーザは対象時刻において、直前の位置登録信号に係るセクタに確率「0.4」で在圏し、直後の位置登録信号に係るセクタに確率「0.6」で在圏していることを意味し、換言すれば、当該ユーザについて、直前の位置登録信号に係るセクタに「0.4人」、直後の位置登録信号に係るセクタに「0.6人」が在圏していると捉えることができる。このような直前・直後の位置登録信号のウエイトの算出を各ユーザについて行い、集計対象全ユーザについての直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトをセクタごとに集計することにより、集計対象全ユーザに関する対象時刻におけるセクタ単位の人口分布を算出することができる。
【００１２】
上記の人口分布算出部は、第１の時間差および第２の時間差のうち少なくとも一方が、予め定められた基準値より長い場合、当該基準値より長い時間差に対応する位置登録信号のウエイトを０としてもよい。なお、上記「基準値」としては、例えば、周期的位置登録の周期（５４分）を採用してもよいし、また、通信中に周期位置登録が発生する場合を考慮し、周期的位置登録の周期（５４分）より長めに設定された時間を採用してもよいし、また、集計対象全ユーザの平均位置登録間隔を採用してもよい。このとき、第１の時間差および第２の時間差のうち一方のみが、基準値より長い場合、人口分布算出部は、基準値以下の時間差に対応する位置登録信号のウエイトを１としてもよいし、別の態様として、人口分布算出部は、基準値と当該基準値以下の時間差との差分を当該基準値により除算して得られた値を、当該基準値以下の時間差に対応する位置登録信号のウエイトとして算出してもよい。【００１３】
ある態様の情報分析装置は、位置登録が行われたセクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を外部から受信する受信部と、各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出部と、各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出部と、を備え、抽出部は、各ユーザについての位置登録信号のうち、端末が位置登録エリアを跨いだことに起因して生成された位置登録信号を除去する前の位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号である直前の除去前位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号である直後の除去前位置登録信号とを抽出するとともに、各ユーザについての位置登録信号のうち、端末が位置登録エリアを跨いだことに起因して生成された位置登録信号を除去した後の位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号である直前の除去後位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号である直後の除去後位置登録信号とを抽出し、人口分布算出部は、各ユーザについて、直前の除去前位置登録信号の位置登録時刻と対象時刻との第１の時間差、直後の除去前位置登録信号の位置登録時刻と対象時刻との第２の時間差、直前の除去後位置登録信号の位置登録時刻と対象時刻との第３の時間差、および、直後の除去後位置登録信号の位置登録時刻と対象時刻との第４の時間差を算出し、各ユーザについて、第１の時間差と第２の時間差に基づいて、直前の除去前位置登録信号のウエイトおよび直後の除去前位置登録信号のウエイトを算出するとともに、第３の時間差と第４の時間差に基づいて、直前の除去後位置登録信号のウエイトおよび直後の除去後位置登録信号のウエイトを算出し、予め定められた按分パラメータに基づいて、各ユーザについて、直前および直後の除去前位置登録信号のウエイトと直前および直後の除去後位置登録信号のウエイトとを按分し、集計対象全ユーザについての、按分後の直前および直後の除去前位置登録信号のウエイト、並びに按分後の直前および直後の除去後位置登録信号のウエイトを、セクタごとに集計することを特徴とする。
【００１４】
なお、抽出部により抽出された直前の除去前位置登録信号と直後の除去前位置登録信号の両方とも、端末が位置登録エリアを跨いだことに起因して生成された位置登録信号でない場合は、直前の除去前位置登録信号と直前の除去後位置登録信号とが同一となり、直後の除去前位置登録信号と直後の除去後位置登録信号とが同一となる。従って、この場合は、直前および直後の除去後位置登録信号の抽出、これら直前および直後の除去後位置登録信号に係る第３、第４の時間差の算出、これら直前および直後の除去後位置登録信号それぞれのウエイト算出、並びに、按分パラメータに基づく按分処理は、不要となる。この場合、人口分布算出部は、各ユーザについて、直前の除去前位置登録信号の位置登録時刻と対象時刻との第１の時間差、および、直後の除去前位置登録信号の位置登録時刻と対象時刻との第２の時間差を算出し、各ユーザについて、第１の時間差と第２の時間差に基づいて、直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトを算出し、集計対象全ユーザについての直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトをセクタごとに集計することによって、対象時刻におけるセクタ単位の人口分布を算出すればよい。【００１５】
別の態様の情報分析装置は、位置登録が行われたセクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を外部から受信する受信部と、各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出部と、各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出部と、を備え、人口分布算出部は、各ユーザについて、対象時刻の直前の位置登録信号に含まれるセクタＩＤに対応する位置と直後の位置登録信号に含まれるセクタＩＤに対応する位置とを結ぶ線上に、直前の位置登録時刻から直後の位置登録時刻まで等速移動した場合の所定時間おきの通過点に相当する複数の点を生成し、生成された点が含まれるセクタに対し仮想の位置登録信号を生成し、対象とする全ユーザ（以下「集計対象全ユーザ」という）についての仮想の位置登録信号から対象時刻の位置登録信号を抽出し、抽出された対象時刻の位置登録信号数をセクタごとに集計することによって、対象時刻におけるセクタ単位の人口分布を算出することを特徴とする。
【００１６】
また、別の態様の情報分析装置は、位置登録が行われたセクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を外部から受信する受信部と、各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出部と、各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出部と、を備え、人口分布算出部は、各ユーザについて、対象時刻の直前の位置登録信号に含まれるセクタＩＤに対応する位置と直後の位置登録信号に含まれるセクタＩＤに対応する位置とを結ぶ線上に、直前の位置登録時刻から直後の位置登録時刻まで等速移動した場合の所定時間おきの通過点に相当する複数の点を生成し、対象時刻に対応する点が含まれるセクタに対し仮想の位置登録信号を生成し、生成された対象時刻の位置登録信号数をセクタごとに集計することにより、対象時刻におけるセクタ単位の人口分布を算出することを特徴とする。
【００１７】
また、別の態様の情報分析装置は、位置登録が行われたセクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を外部から受信する受信部と、各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出部と、各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出部と、を備え、人口分布算出部は、境界がセクタ境界と対応付け可能な領域における領域内移動速度の領域間相対比率に相関する比率に基づいて、各ユーザについて、対象時刻の直前の位置登録信号に含まれるセクタＩＤに対応する位置と直後の位置登録信号に含まれるセクタＩＤに対応する位置とを結ぶ線上に、直前の位置登録時刻から直後の位置登録時刻まで移動した場合の所定時間おきの通過点に相当する複数の点を生成し、生成された点が含まれるセクタに対し仮想の位置登録信号を生成し、集計対象全ユーザについての仮想の位置登録信号から対象時刻の位置登録信号を抽出し、抽出された対象時刻の位置登録信号数をセクタごとに集計することにより、対象時刻におけるセクタ単位の人口分布を算出することを特徴とする。
【００１８】
また、別の態様の情報分析装置は、位置登録が行われたセクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を外部から受信する受信部と、各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出部と、各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出部と、を備え、人口分布算出部は、境界がセクタ境界と対応付け可能な領域における領域内移動速度の領域間相対比率に相関する比率に基づいて、各ユーザについて、対象時刻の直前の位置登録信号に含まれるセクタＩＤに対応する位置と直後の位置登録信号に含まれるセクタＩＤに対応する位置とを結ぶ線上に、直前の位置登録時刻から直後の位置登録時刻まで移動した場合の所定時間おきの通過点に相当する複数の点を生成し、対象時刻に対応する点が含まれるセクタに対し仮想の位置登録信号を生成し、生成された対象時刻の位置登録信号数をセクタごとに集計することにより、対象時刻におけるセクタ単位の人口分布を算出することを特徴とする。
【００１９】
さらに別の態様として、人口分布算出部は、各ユーザについて、直前の位置登録時刻と対象時刻との第１の時間差、および、直後の位置登録時刻と対象時刻との第２の時間差を算出し、各ユーザについて、第１の時間差と第２の時間差に基づいて、直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトを算出し、集計対象全ユーザについての直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトをセクタごとに集計することによって、対象時刻におけるセクタ単位の人口分布を算出することができる。なお、上記別の態様における人口分布算出部は、第１の時間差と第２の時間差の和により第２の時間差を除算して得られた値を、直前の位置登録信号のウエイトとして算出し、第１の時間差と第２の時間差の和により第１の時間差を除算して得られた値を、直後の位置登録信号のウエイトとして算出してもよい。【００２０】
なお、上記情報分析装置は、ユーザごとの属性情報を保持した属性情報保持部をさらに備え、抽出部は、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前および直後の位置登録信号の抽出を行ってもよい。【００２１】
別の態様として、上記情報分析装置は、ユーザごとの属性情報を保持した属性情報保持部をさらに備え、人口分布算出部は、属性情報に応じて対象とするユーザを限定し、当該限定したユーザを対象として対象時刻におけるセクタ単位の人口分布の算出を行ってもよい。【００２２】
上述した情報分析装置に係る発明は、情報分析装置により実行される情報分析方法に係る発明として捉え、以下のように記述することができる。【００２３】
本件発明の一側面に係る情報分析方法は、情報分析装置により実行される情報分析方法であって、位置登録が行われたセクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を外部から受信する受信ステップと、各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出ステップと、各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタＩＤに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出ステップと、を備え、人口分布算出ステップにて、情報分析装置は、各ユーザについて、直前の位置登録時刻と対象時刻との第１の時間差、および、直後の位置登録時刻と対象時刻との第２の時間差を算出し、各ユーザについて、第１の時間差と第２の時間差に基づいて、直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトを算出し、集計対象全ユーザについての直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトをセクタごとに集計することにより、対象時刻におけるセクタ単位の人口分布を算出することを特徴とする。
【００２４】
さらに出願人は、位置登録信号を用いた人口分布の導出に係る発明を拡張し、位置登録信号に限定されることなく、ユーザの位置を示す位置情報、当該位置情報が得られた測位時刻情報およびユーザＩＤを含むポイントデータを用いた人口分布の導出に対して広く適用できる以下のような発明をするに至った。【００２５】
本件発明の一側面に係る情報分析装置は、ユーザの位置を示す位置情報、当該位置情報が得られた測位時刻情報およびユーザＩＤを含むポイントデータを外部から受信する受信部と、各ユーザについてのポイントデータから、測位時刻が対象時刻の直前のポイントデータと、測位時刻が対象時刻の直後のポイントデータとを抽出する抽出部と、各ユーザについて、対象時刻の直前のポイントデータの示す位置と対象時刻の直後のポイントデータの示す位置との間を補完することにより、対象時刻におけるユーザの位置を推定し、推定した各ユーザの位置に基づいて対象時刻における所定の算出対象エリア単位の人口分布を算出する人口分布算出部と、を備え、人口分布算出部は、各ユーザについて、対象時刻の直前のポイントデータに含まれる位置情報の示す位置と直後のポイントデータに含まれる位置情報の示す位置とを結ぶ線上に、直前の測位時刻から直後の測位時刻まで移動した場合の所定時間おきの通過点に相当する複数の点を生成して、対象時刻に対応する点を対象時刻におけるユーザの位置として推定し、推定した各ユーザの位置に基づいて対象時刻における所定の算出対象エリア単位の人口分布を算出することを特徴とする。
【００２６】
上記情報分析装置では、受信部が、ユーザの位置を示す位置情報、当該位置情報が得られた測位時刻情報およびユーザＩＤを含むポイントデータを外部から受信すると、抽出部は、各ユーザについてのポイントデータから、測位時刻が対象時刻の直前のポイントデータと、測位時刻が対象時刻の直後のポイントデータとを抽出し、そして、人口分布算出部は、各ユーザについて、対象時刻の直前のポイントデータの示す位置と対象時刻の直後のポイントデータの示す位置との間を補完することにより、対象時刻におけるユーザの位置を推定し、推定した各ユーザの位置に基づいて対象時刻における所定の算出対象エリア単位の人口分布を算出する。さらに、人口分布算出部は、各ユーザについて、対象時刻の直前のポイントデータに含まれる位置情報の示す位置と直後のポイントデータに含まれる位置情報の示す位置とを結ぶ線上に、直前の測位時刻から直後の測位時刻まで移動した場合の所定時間おきの通過点に相当する複数の点を生成して、対象時刻に対応する点を対象時刻におけるユーザの位置として推定し、推定した各ユーザの位置に基づいて対象時刻における所定の算出対象エリア単位の人口分布を算出する。
【００２７】
以上のように、各ユーザについての対象時刻の直前のポイントデータと直後のポイントデータを基礎として、対象時刻における所定の算出対象エリア単位の人口分布を算出することにより、巨視的な人口分布に関する精度の良いデータを簡易に且つ速やかに収集することが可能となる。また、生成される時間間隔が変動するポイントデータから、対象時刻における所定の算出対象エリア単位の人口分布を精度良く算出することができる。なお、所定の算出対象エリアとしては、セクタ、メッシュ、町丁目、市区町村など、さまざまなエリアを採用することができる。【００２９】
所定時間おきの通過点に相当する複数の点を生成するにあたって、人口分布算出部は、上記複数の点を生成しようとする生成エリアにおけるユーザの移動特性に基づいて、該生成エリア内に生成する点の間隔を決定してもよい。上記の生成エリアにおけるユーザの移動特性は、該生成エリアにおけるユーザの平均移動速度であってもよいし、該生成エリアにおけるユーザの平均滞留時間であってもよい。なお、生成エリアは、上記の算出対象エリアと同じものであってもよいし、異なるものであってもよい。【００３０】
なお、上記情報分析装置は、ユーザごとの属性情報を保持した属性情報保持部をさらに備え、抽出部は、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前および直後のポイントデータを抽出してもよい。【００３１】
別の態様として、上記情報分析装置は、ユーザごとの属性情報を保持した属性情報保持部をさらに備え、人口分布算出部は、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として対象時刻における所定の算出対象エリア単位の人口分布を算出してもよい。【００３２】
なお、前述した人口分布算出部は、さらに、ユーザの端末数を人口に変換するための拡大係数をも基礎として、人口分布を算出してもよい。また、前述した情報分析装置は、対象時刻を入力するための対象時刻入力部をさらに備えてもよいし、また、算出された人口分布情報を出力する出力部をさらに備えてもよい。【００３３】
上記の情報分析装置に係る発明は、情報分析装置により実行される情報分析方法に係る発明として捉え、以下のように記述することができる。即ち、本件発明の一側面に係る情報分析方法は、情報分析装置により実行される情報分析方法であって、ユーザの位置を示す位置情報、当該位置情報が得られた測位時刻情報およびユーザＩＤを含むポイントデータを外部から受信する受信ステップと、各ユーザについてのポイントデータから、測位時刻が対象時刻の直前のポイントデータと、測位時刻が対象時刻の直後のポイントデータとを抽出する抽出ステップと、各ユーザについて、対象時刻の直前のポイントデータの示す位置と対象時刻の直後のポイントデータの示す位置との間を補完することにより、対象時刻におけるユーザの位置を推定し、推定した各ユーザの位置に基づいて対象時刻における所定の算出対象エリア単位の人口分布を算出する人口分布算出ステップと、を備え、人口分布算出ステップにて、情報分析装置は、各ユーザについて、対象時刻の直前のポイントデータに含まれる位置情報の示す位置と直後のポイントデータに含まれる位置情報の示す位置とを結ぶ線上に、直前の測位時刻から直後の測位時刻まで移動した場合の所定時間おきの通過点に相当する複数の点を生成して、対象時刻に対応する点を対象時刻におけるユーザの位置として推定し、推定した各ユーザの位置に基づいて対象時刻における所定の算出対象エリア単位の人口分布を算出することを特徴とする。
【発明の効果】
【００３４】
本発明によれば、各ユーザについての対象時刻の直前のポイントデータと直後のポイントデータを基礎として、対象時刻における所定の算出対象エリア単位の人口分布を算出することにより、巨視的な人口分布に関する精度の良いデータを簡易に且つ速やかに収集することが可能となる。【００３５】
特に、ポイントデータとして位置登録信号を用いた場合には、各ユーザについての対象時刻の直前の位置登録信号と直後の位置登録信号を基礎として、対象時刻におけるセクタ単位の人口分布を算出することにより、巨視的な人口分布に関する精度の良いデータを簡易に且つ速やかに収集することが可能となる。また、このとき、位置登録信号の発生数を基礎とすることなく、各ユーザについての対象時刻の直前の位置登録信号および直後の位置登録信号を基礎として、対象時刻におけるセクタ単位の人口分布を算出するため、生成される時間間隔が変動する位置登録信号から、対象時刻におけるセクタ単位の人口分布を精度良く算出することができる。また、位置登録信号の発生数を基礎とした場合に想定される不都合（例えば、位置登録エリア境界付近で実際の人口数よりもかなり多くの位置登録信号が発生することによる不都合）を防止することができる。
【図面の簡単な説明】
【００３６】
【図１】第１、第２実施形態の通信システムのシステム構成を示す図である。
【図２】図１に示す通信システムの機能構成を示す図である。
【図３】ＢＴＳとセクタとの関係を示す図である。
【図４】第１実施形態の処理を示すフローチャートである。
【図５】第１実施形態におけるユーザごとの位置登録信号の補完処理を示すフローチャートである。
【図６】第１実施形態における点の生成方法を説明するための図である。
【図７】第１実施形態で出力される人口分布図の一例を示す図である。
【図８】第２実施形態の処理を示すフローチャートである。
【図９】ユーザごとの位置登録信号のウエイト算出処理を示すフローチャートである。
【図１０】ウエイト算出方法を説明するための図である。
【図１１】第２実施形態で出力される人口分布図の一例を示す図である。
【図１２】セクタ単位の人口分布とともに位置登録エリア単位の人口分布を表す人口分布図の一例を示す図である。
【図１３】第７実施形態の通信システムのシステム構成を示す図である。
【図１４】図５の処理の変形例を示すフローチャートである。
【図１５】図２の機能構成の変形例を示す図である。
【図１６】第３実施形態におけるユーザごとの位置登録信号の補完処理を示すフローチャートである。
【図１７】第３実施形態における点の生成方法の第１の例を説明するための図である。
【図１８】図１６の処理の変形例を示すフローチャートである。
【図１９】平均移動速度の求め方を説明するための図である。
【図２０】第３実施形態における点の生成方法の第２の例を説明するための図である。
【図２１】平均滞留時間の求め方を説明するための図である。
【図２２】第４実施形態におけるウエイト算出処理を示すフローチャートである。
【図２３】第４、第５実施形態におけるウエイト算出処理を説明するための図である。
【図２４】第５実施形態におけるウエイト算出処理を示すフローチャートである。
【図２５】第６実施形態におけるウエイト算出処理およびウエイト集計処理を示すフローチャートである。
【図２６】第６実施形態におけるウエイト算出処理およびウエイト集計処理を説明するための図である。
【図２７】メッシュとエリア図の合成を説明するための図である。
【図２８】各分割エリアの面積および面積比の算出を説明するための図である。
【図２９】あるメッシュ内の分割エリアの人口の総和算出を説明するための図である。
【図３０】屋内局の通信領域および電波到達範囲が異なる周波数帯を利用する複数の屋外局の通信領域が地理的に重複して存在する環境における変換処理を説明するための図である。
【図３１】各セクタにおける推定人口および性別・年齢層・住所といった属性ごとの人口を出力した出力例を示す図である。
【発明を実施するための形態】
【００３７】
添付図面を参照しながら本発明の実施形態を説明する。可能な場合には、同一の部分には同一の符号を付して、重複する説明を省略する。
【００３８】
［第１実施形態］
［通信システムの構成］
図１は、本実施形態の通信システム１０のシステム構成を示す図である。図１に示すように、この通信システム１０は、移動機１００、ＢＴＳ（基地局）２００、ＲＮＣ（無線ネットワーク制御装置）３００、交換機４００、および管理センタ５００を含んで構成されている。また、この管理センタ５００は、社会センサユニット５０１、ペタマイニングユニット５０２、モバイルデモグラフィーユニット５０３、および可視化ソリューションユニット５０４から構成されている。
【００３９】
ＲＮＣ３００は、ＢＴＳ２００を介して移動機１００が送信した位置登録信号を受信し、位置登録信号の信号数をカウントする。なお、標準仕様書「Radio  Resource  Control  (RRC)  Protocol  Specification：  3GPP  TS  25.331」によれば、ＲＮＣ３００においてセクタ単位での信号数カウントについて規定されており、本方法はこれに準拠するものである。
【００４０】
具体的には、例えば、RRC  Connection  Request信号でパラメタRegistrationが設定されているものを計測してもよい。或いは、更に上位の信号内容を確認してもよい。
【００４１】
交換機４００は、ＢＴＳ２００、ＲＮＣ３００を介して、移動機１００が送信した位置登録信号を収集する。ＲＮＣ３００が、セクタ単位で位置登録信号をカウントしているのに対して、交換機４００は、位置登録エリア単位で移動機１００等を管理しており、移動機１００が送信した位置登録信号を収集することにより位置登録エリアに登録している移動機１００の数を位置登録エリア毎に把握及び記憶している。交換機４００は記憶している位置登録エリアに登録している移動機１００の数（登録数）を、所定のタイミング、または管理センタ５００からの要求に応じて管理センタ５００に出力する。
【００４２】
管理センタ５００は、上述したとおり、社会センサユニット５０１、ペタマイニングユニット５０２、モバイルデモグラフィーユニット５０３、および可視化ソリューションユニット５０４を含んで構成されており、各ユニットでは、移動機１００が送信した位置登録信号を用いた統計処理を行う。
【００４３】
社会センサユニット５０１は、各交換機４００から交換機４００が記憶している位置登録エリアに登録している移動機１００の数等を含んだデータを収集するサーバ装置である。この社会センサユニット５０１は、交換機４００から定期的に出力されたデータを受信し、及び社会センサユニット５０１において予め定められたタイミングに従って交換機４００からデータを取得するように構成されている。
【００４４】
ペタマイニングユニット５０２は、社会センサユニット５０１から受信したデータを所定のデータ形式に変換するサーバ装置である。例えば、ペタマイニングユニット５０２は、ユーザＩＤをキーに、或いはエリアごとにソーティング処理を行う。
【００４５】
モバイルデモグラフィーユニット５０３は、ペタマイニングユニット５０２において処理されたデータに対する集計処理、すなわち各項目のカウンティング処理を行うサーバ装置である。例えば、モバイルデモグラフィーユニット５０３は、後述するような位置登録信号を用いた人口分布の導出等の処理を行うことができる。
【００４６】
可視化ソリューションユニット５０４は、モバイルデモグラフィーユニット５０３において集計処理されたデータを可視可能に処理するサーバ装置である。例えば、可視化ソリューションユニット５０４は、集計されたデータを地図上にマッピング処理することができる。この可視化ソリューションユニット５０４にて処理されたデータは、企業、官公庁または個人等に提供され、店舗開発、道路交通調査、災害対策、環境対策などに利用される。なお、このように統計処理された情報は、当然にプライバシーを侵害しないように個人等は特定されないように加工されている。
【００４７】
なお、社会センサユニット５０１、ペタマイニングユニット５０２、モバイルデモグラフィーユニット５０３および可視化ソリューションユニット５０４はいずれも、前述したようにサーバ装置により構成され、図示は省略するが、通常の情報処理装置の基本構成（即ち、ＣＰＵ、ＲＡＭ、ＲＯＭ、キーボードやマウス等の入力デバイス、外部との通信を行う通信デバイス、情報を記憶する記憶デバイス、および、ディスプレイやプリンタ等の出力デバイス）を備えることは言うまでもない。
【００４８】
図２に通信システム１０の機能構成を示す。図２に示すように通信システム１０は、複数のＢＴＳ２００それぞれが制御するセクタに在圏する複数の移動機１００と、ＢＴＳ２００を制御するＲＮＣ３００と、交換機４００と、情報分析装置６００と、を含んで構成される。情報分析装置６００は、前述した図１に示したモバイルデモグラフィーユニット５０３および可視化ソリューションユニット５０４に相当する。図１の社会センサユニット５０１およびペタマイニングユニット５０２に対応する機能に関しては、図２ではこれらの表記を省略している。
【００４９】
まず、ＲＮＣ３００について説明する。ＲＮＣ３００は、ＲＮＣ通信制御部３０２、位置登録信号受信部３０３、及び信号数計測部３０４を含んで構成されている。このうちＲＮＣ通信制御部３０２は、ＢＴＳ２００を介して移動機１００と通信接続を行う部分であり、例えば、移動機１００からの発信処理に基づいた通信接続処理および位置登録要求に基づいた通信接続処理を行う。本実施形態においては、移動機１００は、例えば、移動機１００が位置の登録を要求する場合の単位領域である位置登録エリアを跨いで移動した場合、移動機１００の電源オン等により移動機１００をネットワークに登録するアタッチ処理を行う場合、移動機１００の電源オフ等により移動機１００の登録状態を削除するデタッチ処理を行う場合、一定の周期（例えば５４分）毎、などのタイミングで、位置登録エリアへの位置登録を求める信号である位置登録信号を送信している。更に、本実施形態では、ＲＮＣ通信制御部３０２は、通信接続処理に用いられるInitial  UE  Messageを交換機４００に送信することができる。なお、このInitial  UE  Messageは、発信または位置登録要求を示す指示情報（位置登録信号）、及び移動機１００を一意に特定するテンポラリＩＤ等のＩＤ、を含んでいる。また、Initial  UE  Messageは、移動機１００の位置情報を付加することも可能である。なお、テンポラリＩＤとは、移動機１００がネットワークに接続した際に、交換機４００により払い出されたＩＤ情報である。
【００５０】
図３は、ＢＴＳ２００とセクタとの関係を示す図である。円形で示されている領域の中心にＢＴＳ２００が位置するものであり、それを中心に複数に等分されたものがセクタである。例えば図３では、ＢＴＳ２００の通信エリアは最大６セクタからなるものであり、それぞれのセクタごとに、セクタを一意に識別可能であるセクタＩＤ（セクタ識別子）が割り当てられており、ＲＮＣ３００は、該セクタＩＤを用いることにより移動機１００が何れのセクタにいるかを、ＢＴＳ２００経由で把握可能となっている。
【００５１】
なお、ＲＮＣ３００は、ＲＲＣコネクション要求の処理を行った際に得られる信号の遅延に基づいてセクタ内のどの位置に移動機１００が位置するか、そのＧＡＩ（Geographical  Area  ID）を算出することもできる。そのセクタＩＤとセクタ内の位置に基づいて移動機１００の位置を特定することもできる。
【００５２】
位置登録信号受信部３０３は、移動機１００が送信した位置登録信号をＲＮＣ通信制御部３０２経由で受信する部分である。
【００５３】
信号数計測部３０４は、位置登録信号受信部３０３が受信した位置登録信号から、上記方法によってセクタ毎に位置登録信号があった回数の合計であるセクタ単位信号数（信号数：すなわち該当するセクタにおいて移動機１００が単位時間に送信した位置登録信号の数）を計測する部分である。
【００５４】
信号数計測部３０４は、ＲＮＣ通信制御部３０２を介してセクタ単位信号数を交換機４００へ送信する。
【００５５】
次に、交換機４００について説明する。交換機４００は、交換機通信制御部４０１、変換部４０２、位置登録信号処理部４０４、および記憶部４０３を含んで構成されている。このうち交換機通信制御部４０１は、ＲＮＣ３００から送信されるInitial  UE  Messageを受信し、このInitial  UE  Messageを用いて通信接続処理を行う部分である。
【００５６】
変換部４０２は、交換機通信制御部４０１により受信されたInitial  UE  Messageに含まれているテンポラリＩＤ等のＩＤを電話番号に変換する部分である。変換部４０２は、変換処理に際して、加入者プロファイル情報を記憶する加入者プロファイル情報記憶部（図示せず）から、テンポラリＩＤ等のＩＤに対応付けられている電話番号を抽出し、当該抽出した電話番号に変換する。なお、この加入者プロファイル情報記憶部は、例えば図示しないＨＬＲ（Home  Location  Register）に備えられているものであり、ここではテンポラリＩＤ等のＩＤと電話番号とを対応付けて管理及び記憶している。
【００５７】
位置登録信号処理部４０４は、位置登録信号を移動機１００からＢＴＳ２００を介して受信し、受信した位置登録信号に基づいて位置登録エリアに登録している移動機１００の数である登録数を計測する部分である。
【００５８】
上述の通り本実施形態では移動機１００は、移動機１００が位置登録エリアを跨いで移動した場合に位置登録信号を送信している。それにより、位置登録信号処理部４０４は、位置登録エリア内に存在する移動機１００の実数を把握できることとなる。なお、標準仕様書「Mobile  Application  Part(MAP)  specification:  3GPP  TS  29.002」には、交換機４００が位置登録を管理する方法について記載されている。本実施形態の交換機４００における処理はこの方法に準拠したものである。
【００５９】
記憶部４０３は、位置登録信号処理部４０４が計測した登録数、及び交換機通信制御部４０１を介してＲＮＣ３００から受信したセクタ単位信号数を入力し、記憶する部分である。なお、変換部４０２で変換された電話番号とInitial  UE  Messageに含まれている移動機１００の位置情報と当該位置情報が測位された時刻とを対応付けて記憶することもできる。なお、情報分析装置６００は、例えば上記のようにして取得された移動機１００の位置情報に基づいて、後述の人口分布の導出処理を行う。ただし、位置情報の取得方法は、前述した処理に限定されるものではなく、また、位置情報として、位置登録信号以外に、ユーザのＧＰＳ位置情報を対象とすることもできる。
【００６０】
次に、情報分析装置６００について説明する。情報分析装置６００は、本発明に関連する構成要素として、情報分析装置通信制御部６０１、保持部６０２、抽出部６０３、人口分布算出部６０４、出力部６０５、及び対象時刻入力部６０６を含んで構成されている。このうち情報分析装置通信制御部６０１は、情報分析装置６００と交換機４００との間の通信を制御する部分であり、後述する人口分布の導出処理のために、位置登録が行われたセクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を、交換機通信制御部４０１およびＲＮＣ通信制御部３０２経由で受信する。
【００６１】
保持部６０２は、情報分析装置通信制御部６０１により受信された位置登録信号を保持する部分である。また、保持部６０２は、ユーザごとの属性情報（例えば年齢、性別、住所など）を保持してもよい。以下では、保持部６０２がユーザごとの属性情報を保持しているものとして説明する。
【００６２】
対象時刻入力部６０６は、いつの時点の人口分布を求めるべきか、即ち、求めるべき人口分布の対象時刻をユーザが入力するための部分である。
【００６３】
抽出部６０３は、後述する人口分布算出の対象となる各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する部分である。なお、後述する人口分布算出の対象となるユーザは、受信した位置登録信号の全ユーザであってもよいし、属性情報に応じて限定したユーザであってもよい。即ち、抽出部６０３は、保持部６０２により保持されたユーザごとの属性情報を参照することで、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前・直後の位置登録信号の抽出を行ってもよい。
【００６４】
人口分布算出部６０４は、後述する処理によって、対象時刻におけるセクタ単位の人口分布を算出する部分である。なお、本実施形態における人口分布算出部６０４は、ＲＮＣ３００によりカウントされたセクタ単位の位置登録信号数や、交換機４００により求められた各位置登録エリアの登録数を用いることなく、後述する処理によって位置登録信号に基づいて、対象時刻におけるセクタ単位の人口分布を算出する。
【００６５】
出力部６０５は、人口分布算出部６０４により算出された人口分布情報を出力する部分である。なお、ここでの「出力」は、表示出力および印刷出力を広く含むことは言うまでもない。即ち、人口分布情報は、ディスプレイ等に表示出力されてもよいし、プリンタ等から印刷出力されてもよいし、表示と印刷の両方で出力されてもよい。
【００６６】
［情報分析装置にて実行される処理について］
次に、上記のような通信システム１０の中で情報分析装置６００により実行される処理について、図４〜図７を参照して説明する。
【００６７】
情報分析装置６００においてユーザが対象時刻入力部６０６から人口分布の対象時刻を入力し、所定の処理開始の操作を行うと、図４の処理が情報分析装置６００において実行開始される。ただし、対象時刻の入力は必須条件ではなく、例えば対象時刻として既定の時刻を予め設定しておいてもよい。また、処理開始については、ユーザが所定の処理開始の操作を行うことは必須ではなく、定時動作などで自動的に図４の処理を実行開始してもよい。
【００６８】
図４のステップＳ１０１では、情報分析装置通信制御部６０１が、位置登録が行われたセクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を、交換機通信制御部４０１およびＲＮＣ通信制御部３０２経由で受信し、保持部６０２が、受信された位置登録信号を保持する。
【００６９】
次のステップＳ１０２では、以下に述べる図５のユーザごとの位置登録信号の補完処理が実行される。抽出部６０３は、保持部６０２に保持された、あるユーザの位置登録信号から、対象時刻直前の位置登録信号（以下「直前の位置登録信号」という）および対象時刻直後の位置登録信号（以下「直後の位置登録信号」という）を抽出し（ステップＳ１０２Ａ）、人口分布算出部６０４は、直前に位置登録が行われたセクタに対応する点（即ち、直前の位置登録信号に含まれたセクタＩＤに対応する位置）と、直後に位置登録が行われたセクタに対応する点（即ち、直後の位置登録信号に含まれたセクタＩＤに対応する位置）との間に等間隔に仮想的に点を生成する（ステップＳ１０２Ｂ）。なお、位置登録が行われたセクタに対応する点（即ち、位置登録信号に含まれたセクタＩＤに対応する位置）としては、セクタ中心位置を採用してもよいし、いわゆるＰＲＡＣＨ−ＰＤ測位計算により取得されるポイント（ＧＡＩ（Geographical  Area  ID））を採用してもよいし、セクタの幾何学的な重心位置を採用してもよい。また、ステップＳ１０２Ｂでは、仮想的な点を、直線上に生成してもよいし、仮想的な地図上での道路や線路などに合せた線上に生成してもよい。ここでは、一例として、位置登録が行われたセクタに対応する点としてセクタ中心位置を採用し、仮想的な点を直線上に生成する例を説明する。
【００７０】
図６に示すように、直前・直後の位置登録信号の時間間隔がｍ分で、ｎ分単位の人口（即ちｎ分間隔の人口）を求める場合、直前の位置登録信号に含まれたセクタＩＤに対応する位置（ここではセクタ中心位置Ａ）と直後の位置登録信号に含まれたセクタＩＤに対応する位置（ここではセクタ中心位置Ｂ）とを結ぶ直線を（ｍ／ｎ）等分し、ｎ分間隔で等間隔に点を仮想的に生成する。例えば、対象時刻が0:30で、あるユーザについての直前の位置登録信号に係る位置登録時刻が0:00、直後の位置登録信号に係る位置登録時刻が0:33で、１分単位の人口（即ち１分間隔の人口）を求める場合、直前の位置登録信号に係るセクタ中心位置Ａと、直後の位置登録信号に係るセクタ中心位置Ｂとを結ぶ直線を３３等分し、１分間隔で等間隔に点を仮想的に生成する。そして、人口分布算出部６０４は、生成した点が含まれるセクタに対し、仮想の位置登録信号を生成する（ステップＳ１０２Ｃ）。以上のような図５の処理は、各ユーザについて実行される。
【００７１】
なお、図５のステップＳ１０２Ａにおいて抽出部６０３は、保持部６０２により保持されたユーザごとの属性情報を参照することで、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前・直後の位置登録信号の抽出を行ってもよい。この場合、以降の処理対象となる位置登録信号を絞り込むことができるため、処理負荷の軽減および処理の高速化を図ることができる。
【００７２】
集計対象全ユーザについて図５の処理が実行完了すると、図４へ戻り、次のステップＳ１０３では、人口分布算出部６０４は、集計対象全ユーザを対象にして、図５の処理で生成された位置登録信号から、対象時刻の位置登録信号を抽出する。なお、上述した例では、位置登録時刻が対象時刻0:30の点が含まれるセクタに対し生成された位置登録信号が抽出される。
【００７３】
そして、次のステップＳ１０４で人口分布算出部６０４は、セクタＩＤごとに集計対象全ユーザの位置登録信号数を集計する。これにより、対象時刻におけるセクタＩＤごとの位置登録信号数が得られ、以下に述べる処理により、対象時刻におけるセクタ単位の人口分布が求められる。
【００７４】
ここで、対象時刻におけるセクタＩＤごとの集計対象全ユーザの位置登録信号数（即ち、セクタＩＤごとのユーザ数）から、対象時刻におけるセクタ単位の人口を求める処理（以下「拡大推計処理」という）について説明する。以下、拡大推計処理の一態様として、ユーザ数を人口に変換するための係数（以下「拡大係数」という）を用いて人口を求める処理を説明する。
【００７５】
上記の拡大係数は、例えばユーザの属性ごとに予め求められており、人口分布算出部６０４は、ユーザの属性（例えば性別、年齢層、住所など）ごとの拡大係数を予め記憶している。ステップＳ１０４で、人口分布算出部６０４は、例えば、集計対象ユーザを属性別に分類し、分類後の属性ごとに、対象時刻におけるセクタＩＤごとの位置登録信号数を集計し、その集計値に、当該属性に対応する拡大係数を乗算することで、対象時刻におけるセクタＩＤごとの人口を属性ごとに求める。このとき、セクタ単位の属性ごとの人口が得られる。さらに、人口分布算出部６０４は、属性ごとの人口をセクタＩＤごとに集計することで、対象時刻におけるセクタ単位の人口を求めてもよい。
【００７６】
また、別の導出方法として、ステップＳ１０４で人口分布算出部６０４は、保持部６０２に保持されたユーザごとの属性情報を参照し、集計対象ユーザの位置登録信号それぞれに対応する拡大係数、即ち、集計対象ユーザそれぞれの属性に対応する拡大係数を読み出し、前述した例のように集計対象ユーザを属性別に分類することなく、集計対象全ユーザについて、セクタＩＤごとに拡大係数を合計することで、対象時刻におけるセクタ単位の人口を求めてもよい。
【００７７】
なお、拡大係数は、ユーザの属性に応じて設定する以外にも、例えば、時間帯や場所に応じて設定してもよい。
【００７８】
また、拡大係数は、予め記憶しておく以外に、以下のようにして求めてもよい。拡大係数は、一例として、「在圏率」と「端末の普及率」との積（即ち、人口に対する在圏数の比率）の逆数を用いることができる。ここで「在圏率」とは、契約台数に対する在圏数の比率を意味し、「端末の普及率」とは人口に対する端末の契約台数の比率を意味する。このような拡大係数は、上記の属性、時間帯、場所といった単位で求めることが望ましいが、必須ではない。拡大係数を求める単位としては、一例として、住所の都道府県ごと、５才刻み年齢層ごと、男女ごと、時間帯として１時間ごとなどを採用してもよい。
【００７９】
次のステップＳ１０５では、出力部６０５が、上記のようにして求められたセクタ単位の人口分布を二次元的に表すセクタ単位の人口分布図を出力する。図７にセクタ単位の人口分布図の一例を示す。この図７では、Ａ〜ＬはセクタＩＤを表し、各セクタＩＤの下の数字は当該セクタの人口を表す。また、別の出力例として、出力部６０５は、図３１に示すように、各セクタについて推定された人口（図３１では「推定人口」と表記）や、各セクタにおける性別・年齢層・住所といった属性ごとの人口を、表形式で出力することができる。また、複数の属性を組み合わせた条件（例えば性別と住所を組み合わせた「東京都在住の女性」といった条件）を採用してもよい。なお、ここでの出力には、表示出力、音声出力、印刷出力等さまざまな出力態様が含まれる。
【００８０】
以上説明した第１実施形態により、各ユーザについての対象時刻の直前の位置登録信号と直後の位置登録信号を基礎として、対象時刻におけるセクタ単位の人口分布を算出することにより、巨視的な人口分布に関する精度の良いデータを簡易に且つ速やかに収集することが可能となる。また、このとき、位置登録信号の発生数を基礎とすることなく、各ユーザについての対象時刻の直前の位置登録信号および直後の位置登録信号を基礎として、対象時刻におけるセクタ単位の人口分布を算出するため、生成される時間間隔が変動する位置登録信号から、対象時刻におけるセクタ単位の人口分布を精度良く算出することができる。また、位置登録信号の発生数を基礎とした場合に想定される不都合（例えば、位置登録エリア境界付近で実際の人口数よりもかなり多くの位置登録信号が発生することによる不都合）を防止することができる。
【００８１】
なお、第１実施形態における図５の補完処理のステップＳ１０２Ｃでは、等間隔に生成した複数の点それぞれが含まれるセクタに対して仮想の位置登録信号を生成した例を説明したが、この変形例として、等間隔に生成した複数の点のうち、対象時刻に対応する点が含まれるセクタに対してだけ、仮想の位置登録信号を生成してもよい。つまり、図１４に示すように、図５のステップＳ１０２Ｃに代わるステップＳ１０２Ｄにおいて、対象時刻に対応する点が含まれるセクタに対して仮想の位置登録信号を生成すればよい。この場合、対象時刻に対応する点が含まれるセクタに対する１つの仮想の位置登録信号のみが生成されるため、図４のステップＳ１０３の処理（対象時刻の位置登録信号の抽出）が不要となり、さらなる処理負荷の軽減および処理の高速化を図ることができる。
【００８２】
［第２実施形態］
前述した第１実施形態の処理は、対象時刻の直前・直後の位置登録信号に対応する点の間に仮想的に生成した複数の点を利用してセクタ単位の人口分布を求めるのに対し、以下に述べる第２実施形態の処理は、対象時刻の直前・直後の位置登録信号のウエイトを利用してセクタ単位の人口分布を求める。
【００８３】
なお、第２実施形態では、図１の通信システム１０のシステム構成、および図２の通信システム１０の機能ブロック構成は第１実施形態と同様であるので、重複した説明は省略する。人口分布算出部６０４により実行される処理内容が異なるので、かかる処理内容を含む第２実施形態の処理について、図８〜図１１を参照して、以下に説明する。
【００８４】
情報分析装置６００においてユーザが対象時刻入力部６０６から人口分布の対象時刻を入力し、所定の処理開始の操作を行うと、図８の処理が情報分析装置６００において実行開始される。ただし、対象時刻の入力は必須条件ではなく、例えば対象時刻として既定の時刻を予め設定しておいてもよい。また、処理開始については、ユーザが所定の処理開始の操作を行うことは必須ではなく、定時動作などで自動的に図８の処理を実行開始してもよい。
【００８５】
図８のステップＳ２０１では、情報分析装置通信制御部６０１が、位置登録が行われたセクタＩＤ、位置登録時刻情報およびユーザＩＤを含む位置登録信号を、交換機通信制御部４０１およびＲＮＣ通信制御部３０２経由で受信し、保持部６０２が、受信された位置登録信号を保持する。
【００８６】
次のステップＳ２０２では、以下に述べる図９のユーザごとの位置登録信号のウエイト算出処理が実行される。抽出部６０３は、保持部６０２に保持された、あるユーザの位置登録信号から、対象時刻直前の位置登録信号（以下「直前の位置登録信号」という）および対象時刻直後の位置登録信号（以下「直後の位置登録信号」という）を抽出し（ステップＳ２０２Ａ）、人口分布算出部６０４は、直前の位置登録時刻と対象時刻との第１の時間差、および、直後の位置登録時刻と対象時刻との第２の時間差を算出する（ステップＳ２０２Ｂ）。
【００８７】
そして、人口分布算出部６０４は、第１の時間差と第２の時間差の和により第２の時間差を除算して得られた値を、直前の位置登録信号のウエイトとして算出し、第１の時間差と第２の時間差の和により第１の時間差を除算して得られた値を、直後の位置登録信号のウエイトとして算出する（ステップＳ２０２Ｃ）。このとき、図１０に示すように、直前の位置登録時刻と対象時刻との第１の時間差がｍ分、直後の位置登録時刻と対象時刻との第２の時間差がｎ分である場合、直前の位置登録信号Ａのウエイトは（ｎ／（ｍ＋ｎ））、直後の位置登録信号Ｂのウエイトは（ｍ／（ｍ＋ｎ））となる。例えば、対象時刻が0:30で、あるユーザについての直前の位置登録信号に係る位置登録時刻が0:00、直後の位置登録信号に係る位置登録時刻が0:40の場合、第１の時間差ｍが３０分、第２の時間差ｎが１０分となるため、直前の位置登録信号Ａのウエイトは（１０／（３０＋１０））＝0.25、直後の位置登録信号Ｂのウエイトは（３０／（３０＋１０））＝0.75となる。なお、直前・直後の位置登録信号のウエイトの算出方法は、上記の方法に限定されるものではなく、第１、第２の時間差を用いた別の数式により位置登録信号のウエイトを算出してもよい。
【００８８】
なお、ウエイトにおける小数点以下の数値については、例えば、小数点以下３けた目以降が存在する場合は同３けた目を四捨五入して小数点以下２けた目までの数値とするなど、予め定めた手法に従ってウエイトを求めればよい。
【００８９】
以上のような図９の処理は、集計対象全ユーザについて実行される。なお、図９のステップＳ２０２Ａにおいて抽出部６０３は、保持部６０２により保持されたユーザごとの属性情報を参照することで、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前・直後の位置登録信号の抽出を行ってもよい。この場合、以降の処理対象となる位置登録信号を絞り込むことができるため、処理負荷の軽減および処理の高速化を図ることができる。
【００９０】
集計対象全ユーザについて図９の処理が実行完了すると、図８へ戻り、次のステップＳ２０３では、人口分布算出部６０４は、セクタＩＤごとに集計対象全ユーザの位置登録信号のウエイトを集計する。これにより、対象時刻におけるセクタＩＤごとの位置登録信号のウエイトの合計値が得られ、第１実施形態で述べた拡大推計処理によって、ウエイトの合計値から対象時刻におけるセクタ単位の人口分布が求められる。
【００９１】
次のステップＳ２０４では、出力部６０５が、上記のようにして求められたセクタ単位の人口分布を二次元的に表すセクタ単位の人口分布図を出力する。図１１にセクタ単位の人口分布図の一例を示す。この図１１では、Ａ〜ＬはセクタＩＤを表し、各セクタＩＤの下の数字は当該セクタの人口を表す。図１１の例は、位置登録信号のウエイトを小数点以下１けた目までの数値とし、その後、セクタＩＤごとに集計対象全ユーザの位置登録信号のウエイトを集計した例を示している。
【００９２】
以上説明した第２実施形態により、各ユーザについての対象時刻の直前の位置登録信号と直後の位置登録信号を基礎として、対象時刻におけるセクタ単位の人口分布を算出することにより、巨視的な人口分布に関する精度の良いデータを簡易に且つ速やかに収集することが可能となる。また、このとき、位置登録信号の発生数を基礎とすることなく、各ユーザについての対象時刻の直前の位置登録信号および直後の位置登録信号を基礎として、対象時刻におけるセクタ単位の人口分布を算出するため、位置登録信号の発生数を基礎とした場合に想定される不都合（例えば、位置登録エリア境界付近で実際の人口数よりもかなり多くの位置登録信号が発生することによる不都合）を防止することができる。
【００９３】
なお、上述した第１、第２実施形態では、セクタ単位の人口分布を算出する例を示したが、さらに、同じ位置登録エリアに属するセクタの人口を位置登録エリアごとに集計する処理を行うことで、位置登録エリア単位の人口分布を算出することができる。例えば、図１１におけるセクタＡ〜Ｅが属する位置登録エリアの人口と、セクタＦ〜Ｌが属する位置登録エリアの人口とを算出し、図１２に示すように、各位置登録エリアに、人口数に応じた表記（色分け、斜線、網掛けなど）を施すことにより、セクタ単位の人口分布とともに位置登録エリア単位の人口分布を表す図を出力することができる。
【００９４】
また、上述した第１、第２実施形態では、ユーザごとの属性情報（例えば年齢、性別、住所など）が保持部６０２により保持される例を示したが、これに限定されるものではなく、図１５に示すように、情報分析装置６００が保持部６０２とは別に、ユーザごとの属性情報を保持する属性情報保持部６０７を備える構成を採用してもよい。
【００９５】
また、上述した第１、第２実施形態では抽出部６０３が、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前・直後の位置登録信号の抽出を行う例を示したが、これに限定されるものではなく、人口分布算出部６０５が、例えば図４のステップＳ１０３、Ｓ１０４や図８のステップＳ２０３において属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として人口分布算出を行ってもよい。このように位置登録信号の抽出時点でなく人口分布算出の際に、属性情報に応じて、対象とするユーザを限定する場合、例えば年齢別や男女別など、さまざまな条件で人口分布算出を行うことができる。もちろん、位置登録信号の抽出と人口分布算出の両方で、属性情報に応じた対象ユーザの限定を行ってもよい。
【００９６】
［第３実施形態］
以下の第３実施形態では、第１実施形態の変形態様として、人口分布算出部６０４が、対象時刻の直前の位置登録信号に含まれたセクタＩＤに対応する位置と、直後の位置登録信号に含まれたセクタＩＤに対応する位置との間に、所定時間おきの通過点に相当する複数の仮想的な点を非等間隔に生成して、仮想の位置登録信号を生成する例を説明する。具体的には、セクタ内移動速度のセクタ間相対比率に相関する比率を予め求めておき、その比率に基づいて非等間隔に仮想的な点を生成して仮想の位置登録信号を生成する実施形態である。
【００９７】
上記の「セクタ内移動速度のセクタ間相対比率に相関する比率」としては、後述するセクタ内移動速度のセクタ間相対比率そのものを用いてもよいし、他の比率（例えば、後述するセクタ内滞留時間のセクタ間相対比率）を用いてもよい。なお、上記の相対比率の単位は、セクタ単位であることは必須ではなく、境界がセクタ境界と対応付け可能な領域の単位（例えばメッシュ、町丁目等）であれば用いることができる。ただし、以下では説明を容易にするために、相対比率の単位がセクタ単位である例を説明する。
【００９８】
（(1)セクタ内移動速度のセクタ間相対比率を用いた例）
まず、「セクタ内移動速度のセクタ間相対比率に相関する比率」として、セクタ内移動速度のセクタ間相対比率を用いた例を説明する。以下では、重複した説明を避けるべく、第１実施形態との差異について重点的に説明する。第３実施形態では、第１実施形態における図５の位置登録信号の補完処理が図１６に示す処理となる。即ち、図１６のステップＳ１０２Ａで対象時刻の直前の位置登録信号および直後の位置登録信号が抽出された後、人口分布算出部６０４は、ステップＳ１０２ＢＸにおいて、直前の位置登録信号に含まれたセクタＩＤに対応する位置と、直後の位置登録信号に含まれたセクタＩＤに対応する位置との間に、予め求められたセクタ内移動速度のセクタ間相対比率に基づいて、所定時間おきの通過点に相当する複数の仮想的な点を生成する。
【００９９】
図１７に示すように、対象時刻の直前の位置登録信号に含まれたセクタＩＤに対応する位置Ａがセクタ１にあり、直後の位置登録信号に含まれたセクタＩＤに対応する位置Ｂがセクタ３にあり、セクタ１と３の間にセクタ２が存在する状況で、予め求められたセクタ内移動速度のセクタ間相対比率（セクタ１：セクタ２：セクタ３）が１：２：３であったとする。ここで、位置Ａ、Ｂを結ぶ線とセクタ１、２の境界線との交差点を点Ｘとし、位置Ａ、Ｂを結ぶ線とセクタ２、３の境界線との交差点を点Ｙとした場合、位置Ａ、Ｂを結ぶ線上に、所定時間おき（例えば１分おき）の通過点に相当する複数の仮想的な点を生成する際に、点の間隔の相対比率（ＡＸ間：ＸＹ間：ＹＢ間）が、上記のセクタ内移動速度のセクタ間相対比率、即ち、１：２：３となるようにする。つまり、セクタ内移動速度がセクタ２、３に比べ相対的に低いセクタ１（ここではＡＸ間）では、点の間隔がセクタ２、３に比べ相対的に狭くなる。一方、セクタ内移動速度がセクタ１、２に比べ相対的に高いセクタ３（ここではＹＢ間）では、点の間隔がセクタ１、２に比べ相対的に広くなる。
【０１００】
このようにしてステップＳ１０２ＢＸで位置Ａ、Ｂを結ぶ線上に仮想的に点を生成した後、人口分布算出部６０４は、生成した点が含まれるセクタに対し、仮想の位置登録信号を生成する（ステップＳ１０２Ｃ）。
【０１０１】
以上のような図１６の処理により、セクタ内移動速度のセクタ間相対比率に基づいて、位置Ａ、Ｂを結ぶ線上に非等間隔に仮想的な点を生成して仮想の位置登録信号を生成することができる。
【０１０２】
なお、図１６の補完処理のステップＳ１０２Ｃの変形例として、非等間隔に生成した複数の点のうち、対象時刻に対応する点が含まれるセクタに対してだけ、仮想の位置登録信号を生成してもよい。つまり、図１８に示すように、図１６のステップＳ１０２Ｃに代わるステップＳ１０２Ｄにおいて、対象時刻に対応する点が含まれるセクタに対して仮想の位置登録信号を生成すればよい。この場合、対象時刻に対応する点が含まれるセクタに対する１つの仮想の位置登録信号のみが生成されるため、さらなる処理負荷の軽減および処理の高速化を図ることができる。
【０１０３】
ここで、セクタ内移動速度のセクタ間相対比率の求め方について概説する。一例としては、セクタごとに、まずユーザごとの移動速度を後述の方法により求め、次に全ユーザの移動速度を平均することで当該セクタのセクタ内平均移動速度を求める。そして、セクタごとのセクタ内平均移動速度同士の比率から、セクタ内移動速度のセクタ間相対比率を求めることができる。なお、セクタ内移動速度のセクタ間相対比率を予め求める際には、位置登録信号より粒度の細かい位置情報（例えばユーザのＧＰＳ位置情報等）を用いることが望ましい。また、セクタ内移動速度のセクタ間相対比率は、特定の時間帯ごとに（例えば朝・昼・晩・深夜などの所定時間帯ごとに）求めてもよいし、曜日ごとに求めてもよいし、平日の平均や休日の平均を求めてもよい。
【０１０４】
ユーザごとの移動速度を求める方法としては、例えば以下のような第１、第２の方法を採用できる。なお、ユーザごとの移動速度を求める対象のセクタを、ここでは「対象セクタ」と称する。
【０１０５】
第１の方法では、対象セクタ内で最初に観測された位置情報および最後に観測された位置情報より、それら２つの位置の間の距離および観測時刻の差分を求め、求めた距離を観測時刻の差分により割り算することで、移動速度を求める。例えば図１９（ａ）に示すように、対象セクタ内で最初に観測された位置情報（位置Ｄおよび観測時刻12:00）と、最後に観測された位置情報（位置Ｅおよび観測時刻12:05）より、ＤＥ間の距離（１５０ｍ）および観測時刻の差分（５分）を求め、距離（１５０ｍ）を観測時刻の差分（５分）により割り算することで、移動速度（1.8ｋｍ／時）を求めることができる。
【０１０６】
第２の方法では、対象セクタ内で最初に観測された位置情報および最後に観測された位置情報に加え、対象セクタに入る直前の位置情報および対象セクタから出た直後の位置情報を用いる。例えば図１９（ｂ）に示すように、対象セクタに入る直前の位置情報（位置Ｃおよび観測時刻11:55）および対象セクタ内で最初に観測された位置情報（位置Ｄおよび観測時刻12:00）より、ＣＤ間を等速直線移動したと仮定した場合の対象セクタに入った時刻（即ち、境界上の点Ｐを通過した時刻）（ここでは仮に11:58とする）と、境界上の点Ｐの位置を求める。同様に、対象セクタ内で最後に観測された位置情報（位置Ｅおよび観測時刻12:05）および対象セクタから出た直後の位置情報（位置Ｆおよび観測時刻12:10）より、ＥＦ間を等速直線移動したと仮定した場合の対象セクタから出た時刻（即ち、境界上の点Ｑを通過した時刻）（ここでは仮に12:06とする）と、境界上の点Ｑの位置を求める。そして、点Ｐ、Ｑの位置よりＰＱ間の距離（１８０ｍ）を求めるとともに、対象セクタに入った時刻と対象セクタから出た時刻との差分（８分）を求め、距離（１８０ｍ）を観測時刻の差分（８分）により割り算することで、移動速度（1.35ｋｍ／時）を求めることができる。
【０１０７】
以上のようにしてユーザごとの移動速度を求めることができ、さらに、全ユーザの移動速度を平均することでセクタ内平均移動速度を求め、セクタごとのセクタ内平均移動速度同士の比率から、セクタ内移動速度のセクタ間相対比率を求めることができる。
【０１０８】
（(2)セクタ内滞留時間のセクタ間相対比率を用いた例）
次に、「セクタ内移動速度のセクタ間相対比率に相関する比率」として、セクタ内滞留時間のセクタ間相対比率を用いた例を説明する。この場合、位置登録信号の補完処理は、前述した図１６、図１８の処理とほぼ同様であるため、重複した説明は省略する。そこで、図１６、図１８のステップＳ１０２ＢＸにおける点の生成方法の詳細について、以下、説明する。
【０１０９】
図２０に示すように、対象時刻の直前の位置登録信号に含まれたセクタＩＤに対応する位置Ａがセクタ１にあり、直後の位置登録信号に含まれたセクタＩＤに対応する位置Ｂがセクタ３にあり、セクタ１と３の間にセクタ２が存在する状況で、予め求められたセクタ内滞留時間のセクタ間相対比率（セクタ１：セクタ２：セクタ３）が２：３：４であったとする。ここで、位置Ａ、Ｂを結ぶ線とセクタ１、２の境界線との交差点を点Ｘとし、位置Ａ、Ｂを結ぶ線とセクタ２、３の境界線との交差点を点Ｙとした場合、位置Ａ、Ｂを結ぶ線上に、所定時間おき（例えば１分おき）の通過点に相当する複数の仮想的な点を生成する際に、点の間隔の相対比率（ＡＸ間：ＸＹ間：ＹＢ間）が、上記のセクタ内滞留時間のセクタ間相対比率の逆数の比率、即ち、(1/2)：(1/3)：(1/4)＝６：４：３となるようにする。つまり、セクタ内滞留時間がセクタ２、３に比べ相対的に短いセクタ１（ここではＡＸ間）では、点の間隔がセクタ２、３に比べ相対的に広くなる。一方、セクタ内滞留時間がセクタ１、２に比べ相対的に長いセクタ３（ここではＹＢ間）では、点の間隔がセクタ１、２に比べ相対的に狭くなる。
【０１１０】
このようにして図１６のステップＳ１０２ＢＸで位置Ａ、Ｂを結ぶ線上に複数の仮想的な点を生成した後、人口分布算出部６０４は、生成した点が含まれるセクタに対し、仮想の位置登録信号を生成する（図１６のステップＳ１０２Ｃ）。又は、図１８のステップＳ１０２ＢＸで位置Ａ、Ｂを結ぶ線上に複数の仮想的な点を生成した後、人口分布算出部６０４は、対象時刻に対応する点が含まれるセクタに対して仮想の位置登録信号を生成する（図１８のステップＳ１０２Ｄ）。
【０１１１】
これにより、セクタ内滞留時間のセクタ間相対比率に基づいて、位置Ａ、Ｂを結ぶ線上に非等間隔に仮想的な点を生成して仮想の位置登録信号を生成することができる。
【０１１２】
ここで、セクタ内滞留時間のセクタ間相対比率の求め方について概説する。一例としては、セクタごとに、まずユーザごとの滞留時間を後述の方法により求め、次に全ユーザの滞留時間を平均することで当該セクタのセクタ内平均滞留時間を求める。そして、セクタごとのセクタ内平均滞留時間同士の比率から、セクタ内滞留時間のセクタ間相対比率を求めることができる。なお、セクタ内滞留時間のセクタ間相対比率を予め求める際には、位置登録信号より粒度の細かい位置情報（例えばユーザのＧＰＳ位置情報等）を用いることが望ましい。また、セクタ内滞留時間のセクタ間相対比率は、特定の時間帯ごとに（例えば朝・昼・晩・深夜などの所定時間帯ごとに）求めてもよいし、曜日ごとに求めてもよいし、平日の平均や休日の平均を求めてもよい。
【０１１３】
ユーザごとの滞留時間を求める方法としては、例えば以下のような第１、第２の方法を採用できる。なお、ユーザごとの滞留時間を求める対象のセクタを、ここでは「対象セクタ」と称する。
【０１１４】
第１の方法では、対象セクタ内で最初に観測された位置情報および最後に観測された位置情報より、これら２つの観測時刻の差分を対象セクタにおける滞留時間として求める。例えば図２１（ａ）に示すように、対象セクタ内で最初に観測された位置情報（位置Ｄおよび観測時刻12:00）と、最後に観測された位置情報（位置Ｅおよび観測時刻12:05）より、これら２つの観測時刻の差分（５分）を対象セクタにおける滞留時間として求めることができる。
【０１１５】
第２の方法では、対象セクタ内で最初に観測された位置情報および最後に観測された位置情報に加え、対象セクタに入る直前の位置情報および対象セクタから出た直後の位置情報を用いる。例えば図２１（ｂ）に示すように、対象セクタに入る直前の位置情報（位置Ｃおよび観測時刻11:55）および対象セクタ内で最初に観測された位置情報（位置Ｄおよび観測時刻12:00）より、ＣＤ間を等速直線移動したと仮定した場合の対象セクタに入った時刻（即ち、境界上の点Ｐを通過した時刻）（ここでは仮に11:58とする）を求める。同様に、対象セクタ内で最後に観測された位置情報（位置Ｅおよび観測時刻12:05）および対象セクタから出た直後の位置情報（位置Ｆおよび観測時刻12:10）より、ＥＦ間を等速直線移動したと仮定した場合の対象セクタから出た時刻（即ち、境界上の点Ｑを通過した時刻）（ここでは仮に12:06とする）を求める。そして、対象セクタに入った時刻と対象セクタから出た時刻との差分（８分）を、対象セクタにおける滞留時間として求めることができる。
【０１１６】
以上のようにしてユーザごとの滞留時間を求めることができ、さらに、全ユーザの滞留時間を平均することでセクタ内平均滞留時間を求め、セクタごとのセクタ内平均滞留時間同士の比率から、セクタ内滞留時間のセクタ間相対比率を求めることができる。
【０１１７】
以上説明した第３実施形態によれば、予め求められたセクタ内移動速度のセクタ間相対比率に相関する比率に基づいて、所定時間おきの通過点に相当する複数の仮想的な点を非等間隔に生成して、仮想の位置登録信号を生成することができる。なお、「セクタ内移動速度のセクタ間相対比率に相関する比率」としてセクタ内移動速度のセクタ間相対比率そのものを用いる場合はセクタの大きさによる影響を受けずに精度良く処理を行うことができるという利点がある。一方、「セクタ内移動速度のセクタ間相対比率に相関する比率」としてセクタ内滞留時間のセクタ間相対比率を用いる場合は、セクタ内移動速度よりも小さい処理負荷で簡易にセクタ内滞留時間を求めることができるため、全体的な処理負荷の削減を図れるという利点がある。
【０１１８】
なお、第１〜第３実施形態では、位置登録情報を用いた人口分布の導出について説明したが、本発明は、位置登録情報以外の位置情報（例えばＧＰＳ位置情報やＰＲＡＣＨ−ＰＤ座標情報等）を用いた人口分布の導出にも適用可能である。具体的には、ユーザの位置を示す位置情報、当該位置情報が得られた測位時刻情報およびユーザＩＤを含むポイントデータが得られる環境であれば、本発明は以下のように適用可能である。
【０１１９】
例えば、図１５の情報分析装置６００では、情報分析装置通信制御部６０１が、ユーザの位置を示す位置情報、当該位置情報が得られた測位時刻情報およびユーザＩＤを含むポイントデータを外部から受信し、抽出部６０３が、各ユーザについてのポイントデータから、測位時刻が対象時刻の直前のポイントデータと測位時刻が対象時刻の直後のポイントデータとを抽出し、そして、人口分布算出部６０４が、各ユーザについて、対象時刻の直前のポイントデータの示す位置と対象時刻の直後のポイントデータの示す位置との間を補完することにより、対象時刻におけるユーザの位置を推定し、推定した各ユーザの位置に基づいて、第１実施形態で述べた拡大推計処理により、対象時刻における所定の算出対象エリア単位の人口分布を算出することができる。このように、各ユーザについての対象時刻の直前のポイントデータと直後のポイントデータを基礎として、対象時刻における所定の算出対象エリア単位の人口分布を算出することにより、巨視的な人口分布に関する精度の良いデータを簡易に且つ速やかに収集することが可能となる。なお、所定の算出対象エリアとしては、セクタに限定されるものではなく、メッシュ、町丁目、市区町村など、さまざまなエリアを採用することができる。
【０１２０】
このとき、人口分布算出部６０４は、例えば、各ユーザについて、対象時刻の直前のポイントデータに含まれる位置情報の示す位置と直後のポイントデータに含まれる位置情報の示す位置とを結ぶ線上に、直前の測位時刻から直後の測位時刻まで移動した場合の所定時間おきの通過点に相当する複数の点を生成して、対象時刻に対応する点を対象時刻におけるユーザの位置として推定し、推定した各ユーザの位置に基づいて、第１実施形態で述べた拡大推計処理により、対象時刻における所定の算出対象エリア単位の人口分布を算出してもよい。
【０１２１】
なお、所定時間おきの通過点に相当する複数の点を生成するにあたって、人口分布算出部６０４は、上記複数の点を生成しようとする生成エリアにおけるユーザの移動特性（例えば該生成エリアにおけるユーザの平均移動速度や平均滞留時間）に基づいて、該生成エリア内に生成する点の間隔を決定してもよい。なお、生成エリアは、上記の算出対象エリアと同じものであってもよいし、異なるものであってもよい。
【０１２２】
また、図１５の情報分析装置６００における抽出部６０３は、属性情報保持部６０７に保持されたユーザごとの属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前および直後のポイントデータを抽出してもよい。もしくは、図１５の情報分析装置６００における人口分布算出部６０４は、属性情報保持部６０７に保持されたユーザごとの属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として対象時刻における所定の算出対象エリア単位の人口分布を算出してもよい。
【０１２３】
［第４実施形態］
以下の第４、第５実施形態は、前述した第２実施形態の変形例に相当し、対象時刻の直前・直後の位置登録信号のウエイトを利用して人口分布を求める処理に関する別の実施形態をそれぞれ説明する。
【０１２４】
既に述べたように、移動機１００は、例えば、移動機１００が位置の登録を要求する場合の単位領域である位置登録エリアを跨いで移動したとき、移動機１００の電源オン等により移動機１００をネットワークに登録するアタッチ処理を行うとき、移動機１００の電源オフ等により移動機１００の登録状態を削除するデタッチ処理を行うとき、周期的位置登録の一定の周期（例えば５４分）毎、などのタイミングで、位置登録エリアへの位置登録を求める位置登録信号を生成する。そのため、移動機１００が電源オフや圏外にいる等に起因した通信不能状態でない限り、移動機１００は、少なくとも上記周期的位置登録の一定の周期（以下「周期ｔ」という）毎に位置登録信号を送信する。この点に着目し、第４、第５実施形態では、対象時刻の直前・直後の位置登録信号のうち、位置登録時刻（即ち、位置登録信号の生成時刻）と対象時刻との時間差が予め定められた基準値（例えばここでは「周期ｔ」）よりも長い位置登録信号については、位置登録時刻と対象時刻との時間差が異常値とみなし、時間差が異常値とみなされた位置登録信号（以下「時間差が異常値である位置登録信号」という）についてのウエイトを０とする。なお、上記「基準値」としては、周期的位置登録の周期（５４分）を採用する以外に、通信中に周期位置登録が発生する場合を考慮し、周期的位置登録の周期（５４分）より長めに設定された時間を採用してもよいし、また、集計対象全ユーザの平均位置登録間隔を採用してもよい。以下では一例として、基準値を周期ｔとした例を説明する。
【０１２５】
以下の第４実施形態では、位置登録時刻と対象時刻との時間差が上記基準値（周期ｔ）以下である位置登録信号（以下、「時間差が正常値である位置登録信号」という）については、ウエイトを１とし、第５実施形態では、時間差が正常値である位置登録信号については、周期ｔと当該時間差との差分に応じてウエイトを算出する（詳細は後述する）。
【０１２６】
まず、第４実施形態について説明する。第４実施形態では、図１の通信システム１０のシステム構成、および図２の通信システム１０の機能ブロック構成は第１実施形態と同様であるので、重複した説明は省略する。人口分布算出部６０４により実行される処理の流れは、第２実施形態で説明した図８、図９の処理の流れと同様であるが、図９のステップＳ２０２Ｃにおけるウエイト算出処理が、第２実施形態とは異なるため、第４実施形態に係るウエイト算出処理について、図２２、図２３を参照して、以下に説明する。なお、図２２はウエイト算出処理のフロー図の一例を示すが、ウエイト算出処理のフローは図２２に示すフローに限定されるものではない。
【０１２７】
図２２に示すように、人口分布算出部６０４は、図９のステップＳ２０２Ｂで算出された「直前の位置登録時刻と対象時刻との第１の時間差ｍ」および「直後の位置登録時刻と対象時刻との第２の時間差ｎ」の少なくとも一方が周期ｔより長いか否かを判断する（図２２のステップＳ３０１）。ここで、第１の時間差ｍ、第２の時間差ｎともに周期ｔ以下であれば、直前・直後の位置登録信号はともに時間差が正常値である位置登録信号と判断できるので、人口分布算出部６０４は、第２実施形態と同様に、第１の時間差ｍと第２の時間差ｎの和により第２の時間差ｎを除算して得られた値を、直前の位置登録信号のウエイトとして算出し、第１の時間差ｍと第２の時間差ｎの和により第１の時間差ｍを除算して得られた値を、直後の位置登録信号のウエイトとして算出する（ステップＳ３０２）。
【０１２８】
一方、ステップＳ３０１にて第１の時間差ｍおよび第２の時間差ｎの少なくとも一方が周期ｔより長い場合は、周期ｔより長い時間差に対応する位置登録信号のウエイトを０とする（ステップＳ３０３）。
【０１２９】
次に、第１の時間差ｍ、第２の時間差ｎがともに周期ｔより長いか否かを判断し（ステップＳ３０４）、第１の時間差ｍ、第２の時間差ｎがともに周期ｔより長ければ、図２２の処理を終了する。この場合、ステップＳ３０３の処理により、直前・直後の位置登録信号ともにウエイトは０とされている。
【０１３０】
一方、ステップＳ３０４にて一方の時間差が周期ｔ以下であれば、周期ｔ以下の時間差に対応する位置登録信号（即ち、時間差が正常値である位置登録信号）のウエイトを１とする（ステップＳ３０５Ａ）。
【０１３１】
以上のようにして図２２の処理により、直前・直後の位置登録信号それぞれについてウエイトが設定される。
【０１３２】
ここで、図２３を用いて具体例を説明する。図２３に示すように、ユーザαに関する位置登録信号として、対象時刻（14:00）に対し、直前の位置登録信号（位置登録時刻：13:30、在圏セクタ：エリアＡ）と直後の位置登録信号（位置登録時刻：15:00、在圏セクタ：エリアＢ）とが抽出されたケースを想定する。周期的位置登録の周期ｔは５４分とする。
【０１３３】
このケースでは、直前の位置登録時刻と対象時刻との第１の時間差ｍは３０分、直後の位置登録時刻と対象時刻との第２の時間差ｎは６０分となり、第２の時間差ｎが周期ｔ（５４分）より長いので、直後の位置登録信号についてのウエイトは０となる。一方、第１の時間差ｍは周期ｔ（５４分）以下なので、直前の位置登録信号についてのウエイトは１となる。そのため、ユーザαは対象時刻14:00に、直前の位置登録信号に係る在圏セクタであるエリアＡに１人分のウエイトで存在するという算出結果が得られる。
【０１３４】
以上の第４実施形態により、対象時刻の直前・直後の位置登録信号のうち、位置登録時刻と対象時刻との時間差が周期的位置登録の周期ｔよりも長い位置登録信号について、時間差が異常値である位置登録信号とみなし、ウエイトを０とするため、時間差が異常値である位置登録信号を除外した上で、より精度の良いウエイト算出処理を行うことができる。
【０１３５】
［第５実施形態］
第５実施形態も、第２実施形態における図９のステップＳ２０２Ｃ（ウエイト算出処理）の別の実施形態であるが、第５実施形態では、位置登録時刻と対象時刻との時間差が予め定められた基準値（ここでは第４実施形態と同様に「周期ｔ」とする）以下の位置登録信号（即ち、時間差が正常値である位置登録信号）については、周期ｔと上記時間差との差分に応じてウエイトを算出する。
【０１３６】
第５実施形態では、図１の通信システム１０のシステム構成、および図２の通信システム１０の機能ブロック構成は第１実施形態と同様であるので、重複した説明は省略する。人口分布算出部６０４により実行される処理の流れは、第２実施形態で説明した図８、図９の処理の流れと同様であるが、図９のステップＳ２０２Ｃにおけるウエイト算出処理が、第２実施形態とは異なるため、第５実施形態に係るウエイト算出処理について、図２３、図２４を参照して、以下に説明する。なお、図２４はウエイト算出処理のフロー図の一例を示すが、ウエイト算出処理のフローは図２４に示すフローに限定されるものではない。
【０１３７】
図２４に示すように、人口分布算出部６０４は、図９のステップＳ２０２Ｂで算出された「直前の位置登録時刻と対象時刻との第１の時間差ｍ」および「直後の位置登録時刻と対象時刻との第２の時間差ｎ」の少なくとも一方が周期ｔより長いか否かを判断する（図２４のステップＳ３０１）。ここで、第１の時間差ｍ、第２の時間差ｎともに周期ｔ以下であれば、直前・直後の位置登録信号はともに時間差が正常値である位置登録信号と判断できるので、人口分布算出部６０４は、第２実施形態と同様に、第１の時間差ｍと第２の時間差ｎの和により第２の時間差ｎを除算して得られた値を、直前の位置登録信号のウエイトとして算出し、第１の時間差ｍと第２の時間差ｎの和により第１の時間差ｍを除算して得られた値を、直後の位置登録信号のウエイトとして算出する（ステップＳ３０２）。
【０１３８】
一方、ステップＳ３０１にて第１の時間差ｍおよび第２の時間差ｎの少なくとも一方が周期ｔより長い場合は、周期ｔより長い時間差に対応する位置登録信号のウエイトを０とする（ステップＳ３０３）。
【０１３９】
次に、第１の時間差ｍ、第２の時間差ｎがともに周期ｔより長いか否かを判断し（ステップＳ３０４）、第１の時間差ｍ、第２の時間差ｎがともに周期ｔより長ければ、図２４の処理を終了する。このとき、ステップＳ３０３の処理により、直前・直後の位置登録信号ともにウエイトは０とされている。
【０１４０】
一方、ステップＳ３０４にて一方の時間差が周期ｔ以下であれば、周期ｔ以下の時間差に対応する位置登録信号（即ち、時間差が正常値である位置登録信号）のウエイトを、以下のように、周期ｔと当該時間差との差分に応じて算出する（ステップＳ３０５Ｂ）。ここでは例えば、以下の式（１）を用いて、周期ｔと当該時間差との差分に応じてウエイトを算出する。
ウエイト＝（周期ｔ−当該時間差）／周期ｔ    （１）
以上のようにして図２４の処理により、直前・直後の位置登録信号それぞれについてウエイトが設定される。
【０１４１】
ここで、図２３を用いて具体例を説明する。図２３に示すように、ユーザαに関する位置登録信号として、対象時刻（14:00）に対し、直前の位置登録信号（位置登録時刻：13:30、在圏セクタ：エリアＡ）と直後の位置登録信号（位置登録時刻：15:00、在圏セクタ：エリアＢ）とが抽出されたケースを想定する。周期的位置登録の周期ｔは５４分とする。
【０１４２】
このケースでは、直前の位置登録時刻と対象時刻との第１の時間差ｍは３０分、直後の位置登録時刻と対象時刻との第２の時間差ｎは６０分となり、第２の時間差ｎが周期ｔ（５４分）より長いので、直後の位置登録信号についてのウエイトは０となる。
【０１４３】
一方、第１の時間差ｍは周期ｔ（５４分）以下なので、直前の位置登録信号についてのウエイトは、上記の式（１）により、（５４−３０）／５４＝２４／５４となる。そのため、ユーザαは対象時刻14:00に、直前の位置登録信号に係る在圏セクタであるエリアＡに「２４／５４」人分のウエイトで存在するという算出結果が得られる。
【０１４４】
ここで、図２３の例を用いて、上記のようにしてウエイトを算出する考え方を説明する。例えば図２３のように、対象時刻の直後の位置登録信号が、時間差が異常値である位置登録信号とみなされ、直前のエリアＡでの位置登録信号のみに基づいてウエイトを算出するケースを想定すると、本来は、当該エリアＡでの位置登録信号の生成時（13:30）の５４分後に周期的位置登録の位置登録信号が生成されるはずである。しかし、５４分後に周期的位置登録の位置登録信号が生成されないとすると、５４分後の時点（14:24）では移動機の電源がオフもしくは圏外とみなすことができる。
【０１４５】
このとき、移動機の電源がオフもしくは圏外になった時刻は、上記エリアＡでの位置登録信号の生成時（13:30）から５４分間の時間帯（即ち、13:30〜14:24）におけるいずれかの時刻である。ここで、移動機の電源がオフもしくは圏外になるというイベントの発生する確率が、時系列的に一様であるとすると、例えば、上記時間帯（13:30〜14:24）の途中である対象時刻14:00では、「対象時刻とエリアＡでの位置登録信号の生成時刻との時間差（３０分）／５４分」の確率で、移動機の電源がオフもしくは圏外になっており、一方、「（５４−３０）／５４＝２４／５４」の確率で、移動機がエリアＡに存在する。
【０１４６】
したがって、このような考え方に基づき、上記の例では、対象時刻14:00でユーザαはエリアＡに「２４／５４」人分の重みで存在するものとする。換言すれば、対象時刻14:00でユーザαについて「３０／５４」人分は「移動機の電源がオフもしくは圏外」の扱いにする。
【０１４７】
以上の第５実施形態により、位置登録時刻と対象時刻との時間差が周期的位置登録の周期ｔよりも長い位置登録信号（即ち、時間差が異常値である位置登録信号）を除外した上で、より精度の良いウエイト算出処理を行うことができる。さらに、時間差が正常値である位置登録信号については、周期ｔ（５４分）と「対象時刻〜信号生成時刻の時間差」との差分に応じてウエイトを算出することで、時間差が正常値である位置登録信号のウエイトをより精度良く求めることができる。
【０１４８】
［第６実施形態］
以下の第６実施形態は、前述した第２実施形態の変形例に相当し、位置登録信号のうち、ＬＡ跨りによる位置登録信号を除去する前の位置登録信号から求めたウエイトと、ＬＡ跨りによる位置登録信号を除去した後の位置登録信号から求めたウエイトとを按分し、按分後のウエイトから人口分布を求める処理について説明する。つまり、ＬＡ境界に面したセクタでは、ＬＡ境界に面していないセクタに比べ、ＬＡ跨りによる位置登録信号が生成される確率が相対的に高いため、このようなセクタ間の不均衡に起因した人口分布算出の精度低下を抑えるべく、上記のようにＬＡ跨りによる位置登録信号を除去する前の位置登録信号から求めたウエイトと、ＬＡ跨りによる位置登録信号を除去した後の位置登録信号から求めたウエイトとを按分するものである。
【０１４９】
第６実施形態では、図１の通信システム１０のシステム構成、および図２の通信システム１０の機能ブロック構成は第１実施形態と同様であるので、重複した説明は省略する。人口分布算出部６０４により実行される処理の流れは、第２実施形態で説明した図８の処理の流れと同様であるが、図８のステップＳ２０２のウエイト算出処理およびステップＳ２０３のウエイト集計処理が、第２実施形態とは異なるため、これらの処理について、図２５、図２６を参照して、以下に説明する。
【０１５０】
抽出部６０３は、対象時刻入力部６０６から対象時刻情報を受け取った後、図２５に示すように、各ユーザについて、ＬＡ跨りによる位置登録信号を除去する前の位置登録信号から、対象時刻の直前の位置登録信号（以下「直前の除去前位置登録信号」という）および直後の位置登録信号（以下「直後の除去前位置登録信号」という）を抽出するとともに、ＬＡ跨りによる位置登録信号を除去した後の位置登録信号から、対象時刻の直前の位置登録信号（以下「直前の除去後位置登録信号」という）および直後の位置登録信号（以下「直後の除去後位置登録信号」という）を抽出する（図２５のステップＳ４０１）。
【０１５１】
なお、位置登録信号には、当該信号が生成された要因（ＬＡ跨り、移動機１００の電源オン等によるアタッチ処理の実行、移動機１００の電源オフ等によるデタッチ処理の実行、周期的位置登録など）に応じた情報（例えば識別フラグ）が含まれており、抽出部６０３は、ある位置登録信号がＬＡ跨りによる位置登録信号であるか否かを、例えば当該位置登録信号に含まれる識別フラグによって判別することで、ＬＡ跨りによる位置登録信号を除去することができる。
【０１５２】
抽出された各ユーザについての４つの位置登録信号は人口分布算出部６０４へ渡され、人口分布算出部６０４は、各ユーザについて、直前の除去前位置登録信号の位置登録時刻と対象時刻との第１の時間差ｍ、直後の除去前位置登録信号の位置登録時刻と対象時刻との第２の時間差ｎ、直前の除去後位置登録信号の位置登録時刻と対象時刻との第３の時間差ｑ、および、直後の除去後位置登録信号の位置登録時刻と対象時刻との第４の時間差ｒを算出する（ステップＳ４０２）。
各ユーザについて、第２実施形態と同様の手順で、上記の第１〜第４の時間差から、直前の除去前位置登録信号のウエイトｗ１、直後の除去前位置登録信号のウエイトｗ２、直前の除去後位置登録信号のウエイトｗ３、および、直後の除去後位置登録信号のウエイトｗ４を算出する（ステップＳ４０３）。例えば、以下の式（２）〜（５）によって各ウエイトｗ１〜ｗ４を算出することができる。
直前の除去前位置登録信号のウエイトｗ１＝ｎ／（ｍ＋ｎ）    （２）
直後の除去前位置登録信号のウエイトｗ２＝ｍ／（ｍ＋ｎ）    （３）
直前の除去後位置登録信号のウエイトｗ３＝ｒ／（ｑ＋ｒ）    （４）
直後の除去後位置登録信号のウエイトｗ４＝ｑ／（ｑ＋ｒ）    （５）
【０１５３】
そして、人口分布算出部６０４は、予め定められた按分パラメータｐに基づいて、各ユーザについて、直前および直後の除去前位置登録信号のウエイトと、直前および直後の除去後位置登録信号のウエイトとを按分する（ステップＳ４０４）。ここでの按分パラメータｐおよび按分については、図２６を用いて後述する。
【０１５４】
さらに、人口分布算出部６０４は、集計対象全ユーザについての、按分後の直前および直後の除去前位置登録信号のウエイト、並びに、按分後の直前および直後の除去後位置登録信号のウエイトを、セクタＩＤごとに集計する（ステップＳ４０５）。
【０１５５】
以下、図２６を用いて具体例を説明する。図２６に例示したように、対象時刻14:00の前後の位置登録信号として、エリアＡで13:20に生成された周期位置登録による位置登録信号、エリアＢで13:30に生成されたＬＡ跨りによる位置登録信号、および、エリアＣで14:30に生成された周期位置登録による位置登録信号が存在するケースでは、直前の除去前位置登録信号は、エリアＡで13:20に生成された位置登録信号となり、直前の除去後位置登録信号は、エリアＢで13:30に生成された位置登録信号となる。一方、直後の除去前位置登録信号、直後の除去後位置登録信号はともに、エリアＣで14:30に生成された位置登録信号となる。このとき、上記の第１の時間差ｍ＝３０分、第２の時間差ｎ＝３０分、第３の時間差ｑ＝４０分、第４の時間差ｒ＝３０分となるため、各ウエイトｗ１〜ｗ４は以下のように算出される。
直前の除去前位置登録信号のウエイトｗ１＝３０／（３０＋３０）＝１／２
直後の除去前位置登録信号のウエイトｗ２＝３０／（３０＋３０）＝１／２
直前の除去後位置登録信号のウエイトｗ３＝３０／（４０＋３０）＝３／７
直後の除去後位置登録信号のウエイトｗ４＝４０／（４０＋３０）＝４／７
【０１５６】
ここで、ＬＡ跨りによる位置登録信号を除去する前の位置登録信号から求めたウエイト（直前および直後の除去前位置登録信号のウエイト）ｗ１、ｗ２の集合をｘとし、ＬＡ跨りによる位置登録信号を除去した後の位置登録信号から求めたウエイト（直前および直後の除去後位置登録信号のウエイト）ｗ３、ｗ４の集合をｙとすると、ユーザαの対象時刻14:00の位置を求めるために、上記ｘ、ｙおよび按分パラメータｐ（０＜ｐ＜１）を以下の式（６）に適用することで、上記ｘとｙを按分する。
ｐｘ＋（１−ｐ）ｙ    （６）
【０１５７】
按分パラメータｐは、直前および直後の除去前位置登録信号のウエイトと直前および直後の除去後位置登録信号のウエイトとを按分することで人口分布算出を適正化するためのパラメータであり、その値は、例えば、ユーザの移動速度や、ＬＡ境界におけるＬＡ跨りによる位置登録と周期的位置登録との発生比率などから、予め定められる。
【０１５８】
例えば、按分パラメータｐ＝0.5の場合、図２６の最下段に示すように、ｗ１、ｗ２の集合ｘと、ウエイトｗ３、ｗ４の集合ｙとを、0.5ずつに按分することになる。按分後は、エリアＢで生成された直前の除去前位置登録信号のウエイトが1/4、エリアＣで生成された直後の除去前位置登録信号のウエイトが1/4、エリアＡで生成された直前の除去後位置登録信号のウエイトが3/14、エリアＣで生成された直後の除去後位置登録信号のウエイトが2/7となる。そして、上記の按分後の各ウエイトがセクタＩＤごとに集計される。上記の例では、ユーザαの対象時刻14:00の位置に関するウエイトとして、エリアＡについて3/14、エリアＢについて1/4、エリアＣについて15/28（即ち、1/4と2/7の和）が求められる。
【０１５９】
以上の第６実施形態のように、ＬＡ跨りによる位置登録信号を除去する前の位置登録信号から求めたウエイトと、ＬＡ跨りによる位置登録信号を除去した後の位置登録信号から求めたウエイトとを按分することで、ＬＡ跨りによる位置登録信号の生成される確率がセクタ間（例えばＬＡ境界に面したセクタと面していないセクタとの間）で異なるというセクタ間の不均衡に起因した人口分布算出の精度低下を抑えることができる。
【０１６０】
なお、前述した第４又は第５実施形態のウエイト算出処理は、上記第６実施形態に適用することができる。即ち、図２５のステップＳ４０３のウエイト算出処理に対し、第４又は第５実施形態のウエイト算出処理を適用することができる。その場合、位置登録時刻と対象時刻との時間差が周期的位置登録の周期ｔよりも長い位置登録信号（即ち、時間差が異常値である位置登録信号）を除外した上で、適正なウエイト算出処理を実行できるとともに、上記のようなセクタ間の不均衡に起因した人口分布算出の精度低下を抑えることができる。
【０１６１】
（エリア変換処理について）
また、上記の第１〜第６実施形態により、対象時刻におけるセクタ単位の人口分布が得られるが、集計単位であるセクタ単位の人口分布を、以下のようにして、集計単位とは異なる出力単位（ここでは一例としてメッシュ）ごとの推計値に変換してもよい。以下の変換処理は、人口分布算出部６０４にて実行してもよいし、出力部６０５にて実行してもよいが、以下では、出力部６０５にて実行する例を述べる。
【０１６２】
以下では、(1)単一の周波数帯を利用する屋外局の通信領域が存在する環境における処理と、(2)屋内局の通信領域および電波到達範囲が異なる周波数帯を利用する複数の屋外局の通信領域のうち２つ以上が地理的に重複して存在する環境における処理とを順に説明する。
【０１６３】
まず、(1)単一の周波数帯を利用する屋外局の通信領域が存在する環境における処理について説明する。出力部６０５は、予め記憶したエリア境界情報に基づいて再現されるエリア図（図２７（ａ）参照）と、所定の区画分け規則に基づき再現される２次元的なメッシュ（図２７（ｂ）参照）とを合成し、図２７（ｃ）に示すような合成図を得る。次に、出力部６０５は、上記合成図においてメッシュ境界により各エリアを分割する。例えば、図２８に示すように、図２７（ａ）のエリアＡは、メッシュ境界により４つの分割エリアＡ−１、Ａ−２、Ａ−３、Ａ−４に分割される。そして、出力部６０５は、各分割エリアの面積を算出し、各分割エリアの面積比を算出する。例えば図２８に示すように、分割エリアＡ−１、Ａ−２、Ａ−３、Ａ−４の面積として、それぞれ１０ｍ２、５０ｍ２、１００ｍ２、４０ｍ２が算出されたとすると、分割エリアＡ−１、Ａ−２、Ａ−３、Ａ−４の面積比（例えば百分率）として、５％、２５％、５０％、２０％が算出される。
【０１６４】
なお、出力部６０５が各分割エリアの面積比を算出することは必須ではない。例えば、各分割エリアの面積比は予め求められており、出力部６０５は、各分割エリアの面積比の情報を、情報分析装置６００内の図示しないテーブル又は外部から参照可能とされた構成であってもよい。
【０１６５】
次に、出力部６０５は、各分割エリアの人口を算出する。例えば図２７（ａ）のエリアＡの人口が８００人であったとすると、図２７に示すように分割エリアＡ−２の人口として、２００人（即ち、８００人×２５％）が算出される。同様に、エリアＢ、Ｃの人口がそれぞれ５００人、７５０人であったとすると、エリアＢにおける面積比８０％の分割エリアＢ−１の人口として、４００人（即ち、５００人×８０％）が算出され、エリアＣにおける面積比８０％の分割エリアＣ−４の人口として、６００人（即ち、７５０人×８０％）が算出される。
【０１６６】
さらに、出力部６０５は、１つのメッシュに内包された複数の分割エリアの人口の総和を算出することで、当該メッシュの人口を算出する。図２９の例では、１つのメッシュに内包された分割エリアＡ−２、Ｂ−１、Ｃ−４の人口の総和１２００人（即ち、２００人＋４００人＋６００人）を算出し、この１２００人を当該メッシュの人口とする。
【０１６７】
以上のようにして、単一の周波数帯を利用する屋外局の通信領域が存在する環境において、集計単位ごとの人口を出力単位ごとの人口に変換することができる。
【０１６８】
次に、(2)屋内局の通信領域および電波到達範囲が異なる周波数帯を利用する複数の屋外局の通信領域のうち２つ以上が地理的に重複して存在する環境における処理について説明する。図３０に示すように、屋内局の通信領域および電波到達範囲が異なる周波数帯（屋外2GHz/1.7GHzと屋外800MHz）を利用する複数の屋外局の通信領域が地理的に重複して存在する環境では、出力部６０５は、それぞれの通信領域について、上述した(1)単一の周波数帯を利用する屋外局の通信領域が存在する環境における処理を行うことで、それぞれの通信領域に関する出力単位（メッシュ）ごとの人口を求め、最後に、出力単位ごとに各通信領域に関する人口を合算することで、出力単位ごとの人口を得る。
【０１６９】
図３０の例では、出力部６０５は、まず、屋外2GHz/1.7GHzを利用する屋外局の通信領域、屋外800MHzを利用する屋外局の通信領域、屋内局の通信領域それぞれについて、上述した(1)単一の周波数帯を利用する屋外局の通信領域が存在する環境における処理を行う。例えば、屋外2GHz/1.7GHzを利用する屋外局の通信領域において、出力単位ＱとエリアＡとが重なった分割エリアがエリアＡ全体に対し４０％の面積比であるとすると、エリアＡの推計人口１００人に面積比0.4を乗算することで、出力単位ＱとエリアＡとが重なった分割エリアの推計人口４０人が得られる。同様にして、出力単位ＱとエリアＢとが重なった分割エリアについての推計人口３人（エリアＢの推計人口３０人×面積比0.1）、および、出力単位ＱとエリアＣとが重なった分割エリアの推計人口５人（エリアＣの推計人口１００人×面積比0.05）が得られる。屋外800MHzを利用する屋外局の通信領域についても同様に、出力単位ＱとエリアＤとが重なった分割エリアの推計人口３人（エリアＤの推計人口１０人×面積比0.3）、および、出力単位ＱとエリアＦとが重なった分割エリアの推計人口９人（エリアＦの推計人口３０人×面積比0.3）が得られる。一方、屋内局については、個々の屋内局の電波到達範囲であるエリアは非常に小さくて、図３０の例では、１つの屋内局のエリアＬの全体が出力単位Ｑと重なっているため、１００％の面積比と考えることができる。そこで、エリアＬの推計人口１０人に面積比1.0を乗算することで、出力単位ＱとエリアＬとが重なったエリア（この例ではエリアＬ全体）の推計人口１０人が得られる。
【０１７０】
最後に、出力部６０５は、上記のようにして得られた、出力単位Ｑと各エリアとが重なったエリアの推計人口を合算することで、出力単位Ｑの推計人口７０人を得る。以上のようにして、集計単位ごとの人口から出力単位Ｑの推計人口へ変換することができる。図３０は、１つの出力単位Ｑの推計人口への変換を示しているが、同様の処理を他の出力単位についても実行することで、対象とする全ての出力単位の推計人口への変換を行うことができる。
【０１７１】
以上のようにして、屋内局の通信領域および電波到達範囲が異なる周波数帯を利用する複数の屋外局の通信領域のうち２つ以上が地理的に重複して存在する環境において、集計単位ごとの人口を出力単位ごとの人口に変換することができる。
【０１７２】
なお、上記のエリア変換処理は、セクタ単位の人口を、セクタ単位とは異なる単位（例えばメッシュ単位）の人口に変換する場合だけでなく、第１実施形態で述べた拡大推計処理（端末数を人口に変換する処理）の前に、実行してもよい。つまり、セクタ単位の端末数を、上記のエリア変換処理によって、セクタ単位とは異なる単位（例えばメッシュ単位）の端末数に変換し、その後、変換後の単位（例えばメッシュ単位）の端末数を、拡大推計処理によって人口に変換することで、変換後の単位（例えばメッシュ単位）の人口分布を求めても良い。
【０１７３】
第１〜第６実施形態における人口分布算出では、ユーザの属性ごと、時間帯ごと、場所ごとに、人口分布算出の条件や区分を設定することができる。例えば、人口分布算出に用いられる位置データ（例えば位置登録信号）を、特定の属性を持つユーザの位置データ（例えば位置登録信号）に絞り込むことにより、ユーザの属性ごとに人口分布算出の条件を設定できる。また、対象時刻を特定の時間帯の範囲内に絞り込むことにより、時間帯の条件を設定できる。その際、対象時刻は所定周期に基づく複数の日（例えば平日（月曜日〜金曜日））における特定時刻に設定してもよい。さらに、上記のエリア変換処理によって、セクタ単位の人口（又は端末数）を、セクタ単位とは異なる単位（例えばメッシュ単位）の人口（又は端末数）に変換することで、人口分布算出におけるエリア区分を自在に設定できる。
【０１７４】
［第７実施形態］
図１３は、第７実施形態の通信システム１０ａのシステム構成図である。図１３に示すように、この通信システム１０ａは、新しい通信方式であるＬＴＥ（Long  Term  Evolution）に適用した場合のシステム構成であり、移動機１００、ｅＮＢ（Evolution  Node  B）２５０、交換機４００、および管理センタ５００を含んで構成されている。また、この管理センタ５００は、社会センサユニット５０１、ペタマイニングユニット５０２、モバイルデモグラフィーユニット５０３、および可視化ソリューションユニット５０４から構成されている。なお、ｅＮＢ２５０は、ＢＴＳ２００およびＲＮＣ３００の両方の機能を包含したものである。
【０１７５】
この第７実施形態は、ＬＴＥに適用した場合のシステム構成であって、その具体的な処理内容は、上述の第１〜第６実施形態と同様であるため、その詳細説明は割愛する。なお、第１〜第６実施形態において、そのプロトコルは、ＲＡＮＡＰ（Radio  Access  Network  Application  Part）であるが、第７実施形態においては、ＬＴＥで使用するＳ１ＡＰ（S1  Application  protocol）が用いられており、Initial  UE  Messageについては、Ｓ１ＡＰにおいても同信号が用いられている。
【０１７６】
また、第１〜第７実施形態においては、第三世代携帯電話（３Ｇ）システムを前提として説明していたが、ＧＳＭ（Global  System  for  Mobile  Communications）においても適用可能である。
【符号の説明】
【０１７７】
１０、１０ａ…通信システム、１００…移動機、２００…ＢＴＳ、２５０…ｅＮＢ、３００…ＲＮＣ、３０２…ＲＮＣ通信制御部、３０３…位置登録信号受信部、３０４…信号数計測部、４００…交換機、４０１…交換機通信制御部、４０２…変換部、４０３…記憶部、４０４…位置登録信号処理部、５００…管理センタ、５０１…社会センサユニット、５０２…ペタマイニングユニット、５０３…モバイルデモグラフィーユニット、５０４…可視化ソリューションユニット、６００…情報分析装置、６０１…情報分析装置通信制御部、６０２…保持部、６０３…抽出部、６０４…人口分布算出部、６０５…出力部、６０６…対象時刻入力部、６０７…属性情報保持部。
【図１】
0005497899.tif 000002
【図２】
0005497899.tif 000003
【図３】
0005497899.tif 000004
【図４】
0005497899.tif 000005
【図５】
0005497899.tif 000006
【図６】
0005497899.tif 000007
【図７】
0005497899.tif 000008
【図８】
0005497899.tif 000009
【図９】
0005497899.tif 000010
【図１０】
0005497899.tif 000011
【図１１】
0005497899.tif 000012
【図１２】
0005497899.tif 000013
【図１３】
0005497899.tif 000014
【図１４】
0005497899.tif 000015
【図１５】
0005497899.tif 000016
【図１６】
0005497899.tif 000017
【図１７】
0005497899.tif 000018
【図１８】
0005497899.tif 000019
【図１９】
0005497899.tif 000020
【図２０】
0005497899.tif 000021
【図２１】
0005497899.tif 000022
【図２２】
0005497899.tif 000023
【図２３】
0005497899.tif 000024
【図２４】
0005497899.tif 000025
【図２５】
0005497899.tif 000026
【図２６】
0005497899.tif 000027
【図２７】
0005497899.tif 000028
【図２８】
0005497899.tif 000029
【図２９】
0005497899.tif 000030
【図３０】
0005497899.tif 000031
【図３１】
0005497899.tif 000032