位置登録が行われたセクタID、位置登録時刻情報およびユーザIDを含む位置登録信号を外部から受信する受信部と、各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出部と、各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタIDに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタIDに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出部と、を備える情報分析装置。前記人口分布算出部は、各ユーザについて、直前の位置登録時刻と対象時刻との第1の時間差、および、直後の位置登録時刻と対象時刻との第2の時間差を算出し、各ユーザについて、第1の時間差と第2の時間差に基づいて、直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトを算出し、集計対象全ユーザについての直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトをセクタごとに集計する、ことにより、対象時刻におけるセクタ単位の人口分布を算出する、ことを特徴とする請求項1に記載の情報分析装置。前記人口分布算出部は、第1の時間差と第2の時間差の和により第2の時間差を除算して得られた値を、直前の位置登録信号のウエイトとして算出し、第1の時間差と第2の時間差の和により第1の時間差を除算して得られた値を、直後の位置登録信号のウエイトとして算出する、ことを特徴とする請求項2に記載の情報分析装置。前記人口分布算出部は、前記第1の時間差および前記第2の時間差のうち少なくとも一方が、予め定められた基準値より長い場合、当該基準値より長い時間差に対応する位置登録信号のウエイトを0とする、ことを特徴とする請求項2に記載の情報分析装置。前記人口分布算出部は、前記第1の時間差および前記第2の時間差のうち一方のみが、前記基準値より長い場合、前記基準値以下の時間差に対応する位置登録信号のウエイトを1とする、ことを特徴とする請求項4に記載の情報分析装置。前記人口分布算出部は、前記第1の時間差および前記第2の時間差のうち一方のみが、前記基準値より長い場合、前記基準値と当該基準値以下の時間差との差分を当該基準値により除算して得られた値を、当該基準値以下の時間差に対応する位置登録信号のウエイトとして算出する、ことを特徴とする請求項4に記載の情報分析装置。前記抽出部は、各ユーザについての位置登録信号のうち、端末が位置登録エリアを跨いだことに起因して生成された位置登録信号を除去する前の位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号である直前の除去前位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号である直後の除去前位置登録信号とを抽出するとともに、各ユーザについての位置登録信号のうち、端末が位置登録エリアを跨いだことに起因して生成された位置登録信号を除去した後の位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号である直前の除去後位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号である直後の除去後位置登録信号とを抽出し、前記人口分布算出部は、各ユーザについて、直前の除去前位置登録信号の位置登録時刻と対象時刻との第1の時間差、直後の除去前位置登録信号の位置登録時刻と対象時刻との第2の時間差、直前の除去後位置登録信号の位置登録時刻と対象時刻との第3の時間差、および、直後の除去後位置登録信号の位置登録時刻と対象時刻との第4の時間差を算出し、各ユーザについて、第1の時間差と第2の時間差に基づいて、直前の除去前位置登録信号のウエイトおよび直後の除去前位置登録信号のウエイトを算出するとともに、第3の時間差と第4の時間差に基づいて、直前の除去後位置登録信号のウエイトおよび直後の除去後位置登録信号のウエイトを算出し、予め定められた按分パラメータに基づいて、各ユーザについて、直前および直後の除去前位置登録信号のウエイトと直前および直後の除去後位置登録信号のウエイトとを按分し、集計対象全ユーザについての、按分後の直前および直後の除去前位置登録信号のウエイト、並びに按分後の直前および直後の除去後位置登録信号のウエイトを、セクタごとに集計する、ことを特徴とする請求項1に記載の情報分析装置。前記人口分布算出部は、各ユーザについて、対象時刻の直前の位置登録信号に含まれるセクタIDに対応する位置と直後の位置登録信号に含まれるセクタIDに対応する位置とを結ぶ線上に、直前の位置登録時刻から直後の位置登録時刻まで等速移動した場合の所定時間おきの通過点に相当する複数の点を生成し、生成された点が含まれるセクタに対し仮想の位置登録信号を生成し、集計対象全ユーザについての仮想の位置登録信号から対象時刻の位置登録信号を抽出し、抽出された対象時刻の位置登録信号数をセクタごとに集計する、ことにより、対象時刻におけるセクタ単位の人口分布を算出する、ことを特徴とする請求項1に記載の情報分析装置。前記人口分布算出部は、各ユーザについて、対象時刻の直前の位置登録信号に含まれるセクタIDに対応する位置と直後の位置登録信号に含まれるセクタIDに対応する位置とを結ぶ線上に、直前の位置登録時刻から直後の位置登録時刻まで等速移動した場合の所定時間おきの通過点に相当する複数の点を生成し、対象時刻に対応する点が含まれるセクタに対し仮想の位置登録信号を生成し、生成された対象時刻の位置登録信号数をセクタごとに集計する、ことにより、対象時刻におけるセクタ単位の人口分布を算出する、ことを特徴とする請求項1に記載の情報分析装置。前記人口分布算出部は、境界がセクタ境界と対応付け可能な領域における領域内移動速度の領域間相対比率に相関する比率に基づいて、各ユーザについて、対象時刻の直前の位置登録信号に含まれるセクタIDに対応する位置と直後の位置登録信号に含まれるセクタIDに対応する位置とを結ぶ線上に、直前の位置登録時刻から直後の位置登録時刻まで移動した場合の所定時間おきの通過点に相当する複数の点を生成し、生成された点が含まれるセクタに対し仮想の位置登録信号を生成し、集計対象全ユーザについての仮想の位置登録信号から対象時刻の位置登録信号を抽出し、抽出された対象時刻の位置登録信号数をセクタごとに集計する、ことにより、対象時刻におけるセクタ単位の人口分布を算出する、ことを特徴とする請求項1に記載の情報分析装置。前記人口分布算出部は、境界がセクタ境界と対応付け可能な領域における領域内移動速度の領域間相対比率に相関する比率に基づいて、各ユーザについて、対象時刻の直前の位置登録信号に含まれるセクタIDに対応する位置と直後の位置登録信号に含まれるセクタIDに対応する位置とを結ぶ線上に、直前の位置登録時刻から直後の位置登録時刻まで移動した場合の所定時間おきの通過点に相当する複数の点を生成し、対象時刻に対応する点が含まれるセクタに対し仮想の位置登録信号を生成し、生成された対象時刻の位置登録信号数をセクタごとに集計する、ことにより、対象時刻におけるセクタ単位の人口分布を算出する、ことを特徴とする請求項1に記載の情報分析装置。前記情報分析装置は、ユーザごとの属性情報を保持した属性情報保持部、をさらに備え、前記抽出部は、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前および直後の位置登録信号の抽出を行う、ことを特徴とする請求項1~11の何れか1項に記載の情報分析装置。前記情報分析装置は、ユーザごとの属性情報を保持した属性情報保持部、をさらに備え、前記人口分布算出部は、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として対象時刻におけるセクタ単位の人口分布の算出を行う、ことを特徴とする請求項1~11の何れか1項に記載の情報分析装置。ユーザの位置を示す位置情報、当該位置情報が得られた測位時刻情報およびユーザIDを含むポイントデータを外部から受信する受信部と、各ユーザについてのポイントデータから、測位時刻が対象時刻の直前のポイントデータと、測位時刻が対象時刻の直後のポイントデータとを抽出する抽出部と、各ユーザについて、対象時刻の直前のポイントデータの示す位置と対象時刻の直後のポイントデータの示す位置との間を補完することにより、対象時刻におけるユーザの位置を推定し、推定した各ユーザの位置に基づいて対象時刻における所定の算出対象エリア単位の人口分布を算出する人口分布算出部と、を備える情報分析装置。前記人口分布算出部は、各ユーザについて、対象時刻の直前のポイントデータに含まれる位置情報の示す位置と直後のポイントデータに含まれる位置情報の示す位置とを結ぶ線上に、直前の測位時刻から直後の測位時刻まで移動した場合の所定時間おきの通過点に相当する複数の点を生成して、対象時刻に対応する点を対象時刻におけるユーザの位置として推定し、推定した各ユーザの位置に基づいて対象時刻における所定の算出対象エリア単位の人口分布を算出する、ことを特徴とする請求項14に記載の情報分析装置。前記人口分布算出部は、前記所定時間おきの通過点に相当する複数の点を生成するにあたって、前記複数の点を生成しようとする生成エリアにおけるユーザの移動特性に基づいて、該生成エリア内に生成する点の間隔を決定する、ことを特徴とする請求項15に記載の情報分析装置。前記生成エリアにおけるユーザの移動特性は、該生成エリアにおけるユーザの平均移動速度である、ことを特徴とする請求項16に記載の情報分析装置。前記生成エリアにおけるユーザの移動特性は、該生成エリアにおけるユーザの平均滞留時間である、ことを特徴とする請求項16に記載の情報分析装置。前記情報分析装置は、ユーザごとの属性情報を保持した属性情報保持部、をさらに備え、前記抽出部は、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前および直後のポイントデータを抽出する、ことを特徴とする請求項14~18の何れか1項に記載の情報分析装置。前記情報分析装置は、ユーザごとの属性情報を保持した属性情報保持部、をさらに備え、前記人口分布算出部は、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として対象時刻における所定の算出対象エリア単位の人口分布を算出する、ことを特徴とする請求項14~18の何れか1項に記載の情報分析装置。前記人口分布算出部は、さらに、ユーザの端末数を人口に変換するための拡大係数をも基礎として、前記人口分布を算出する、ことを特徴とする請求項1~20の何れか1項に記載の情報分析装置。対象時刻を入力するための対象時刻入力部、をさらに備える請求項1~21の何れか1項に記載の情報分析装置。算出された前記人口分布情報を出力する出力部、をさらに備える請求項1~22の何れか1項に記載の情報分析装置。情報分析装置により実行される情報分析方法であって、位置登録が行われたセクタID、位置登録時刻情報およびユーザIDを含む位置登録信号を外部から受信する受信ステップと、各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出ステップと、各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタIDに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタIDに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出ステップと、を備える情報分析方法。情報分析装置により実行される情報分析方法であって、ユーザの位置を示す位置情報、当該位置情報が得られた測位時刻情報およびユーザIDを含むポイントデータを外部から受信する受信ステップと、各ユーザについてのポイントデータから、測位時刻が対象時刻の直前のポイントデータと、測位時刻が対象時刻の直後のポイントデータとを抽出する抽出ステップと、各ユーザについて、対象時刻の直前のポイントデータの示す位置と対象時刻の直後のポイントデータの示す位置との間を補完することにより、対象時刻におけるユーザの位置を推定し、推定した各ユーザの位置に基づいて対象時刻における所定の算出対象エリア単位の人口分布を算出する人口分布算出ステップと、を備える情報分析方法。上記のように、従来は巨視的な人口分布に関するデータを収集するだけでも非常に面倒であり、かかるデータを簡易に且つ速やかに収集することは困難であった。一方、携帯端末を利用して人口分布を得る試みとして例えば特許文献1には、GPS機能付き携帯端末を利用して人口の分布を求める点について記載されているが、この技術でも、対象とするユーザの全員に対してユーザ全員を一意に識別可能な識別子が付加されたGPS内蔵機器を配布する、及びGPSを用いて測位した情報を一人一人収集する等の、人口分布を求めるための多大な処理負荷と時間を要する。そのため、より簡易且つ速やかに人口分布に関するデータを収集する技術が待望されていた。また、収集されるデータには、一定レベル以上の精度が求められるのは言うまでもない。なお、上記のようなGPS機能を用いた場合でも、人口分布を求めたい対象時刻きっかりの測位情報が必ずしも得られるとは限らない。本発明は、上記課題に鑑み、巨視的な人口分布に関する精度の良いデータを簡易に且つ速やかに収集することを目的とする。出願人は、上記の目的を達成するために、位置登録信号を用いた人口分布の導出に係る発明について、日本出願特願2009-92225号および国際特許出願PCT/JP2010/055424を既に提出した。その後も、出願人は、位置登録信号を用いた人口分布の導出に係る発明のさらなる改良を進め、今回、同発明を新たな視点から改良した本件発明をするに至った。本件発明の一側面に係る情報分析装置は、位置登録が行われたセクタID、位置登録時刻情報およびユーザIDを含む位置登録信号を外部から受信する受信部と、各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出部と、各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタIDに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタIDに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出部と、を備えることを特徴とする。なお、本件において対象時刻とは、人口分布を求めるべき時刻を意味する。上記情報分析装置では、受信部が、セクタID、位置登録時刻情報およびユーザIDを含む位置登録信号を外部から受信すると、抽出部は、各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出し、そして、人口分布算出部は、各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタIDに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタIDに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する。なお、対象とするユーザは、位置登録信号を受信した全ユーザでもよいし、その中から属性情報例えば年齢、性別、住所などに応じて限定したユーザでもよい。以上のように、各ユーザについての対象時刻の直前の位置登録信号と直後の位置登録信号を基礎として、対象時刻におけるセクタ単位の人口分布を算出することにより、巨視的な人口分布に関する精度の良いデータを簡易に且つ速やかに収集することが可能となる。また、このとき、位置登録信号の発生数を基礎とすることなく、各ユーザについての対象時刻の直前の位置登録信号および直後の位置登録信号を基礎として、対象時刻におけるセクタ単位の人口分布を算出するため、生成される時間間隔が変動する位置登録信号から、対象時刻におけるセクタ単位の人口分布を精度良く算出することができる。また、位置登録信号の発生数を基礎とした場合に想定される不都合例えば、位置登録エリア境界付近で実際の人口数よりもかなり多くの位置登録信号が発生することによる不都合を防止することができる。人口分布算出部は、より具体的には以下のようにして、対象時刻におけるセクタ単位の人口分布を算出することができる。例えば、人口分布算出部は、各ユーザについて、直前の位置登録時刻と対象時刻との第1の時間差、および、直後の位置登録時刻と対象時刻との第2の時間差を算出し、各ユーザについて、第1の時間差と第2の時間差に基づいて、直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトを算出し、集計対象全ユーザについての直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトをセクタごとに集計する。なお、本件において位置登録信号のウエイトとは、あるユーザが対象時刻において当該位置登録信号に係るセクタに在圏している確率を意味する。実際に、対象時刻ちょうどの位置登録信号があるとは限らないため、位置登録信号の位置登録時刻と対象時刻との時間差が短いほど、ユーザが当該位置登録信号に係るセクタに在圏している確率即ち、当該位置登録信号のウエイトが大きくなるように、位置登録信号のウエイトが算出される。例えば、人口分布算出部は、第1の時間差と第2の時間差の和により第2の時間差を除算して得られた値を、直前の位置登録信号のウエイトとして算出し、第1の時間差と第2の時間差の和により第1の時間差を除算して得られた値を、直後の位置登録信号のウエイトとして算出してもよい。例えば、あるユーザについて、直前の位置登録信号のウエイトとして04、直後の位置登録信号のウエイトとして06が算出された場合、当該ユーザは対象時刻において、直前の位置登録信号に係るセクタに確率04で在圏し、直後の位置登録信号に係るセクタに確率06で在圏していることを意味し、換言すれば、当該ユーザについて、直前の位置登録信号に係るセクタに04人、直後の位置登録信号に係るセクタに06人が在圏していると捉えることができる。このような直前直後の位置登録信号のウエイトの算出を各ユーザについて行い、集計対象全ユーザについての直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトをセクタごとに集計することにより、集計対象全ユーザに関する対象時刻におけるセクタ単位の人口分布を算出することができる。上記の人口分布算出部は、第1の時間差および第2の時間差のうち少なくとも一方が、予め定められた基準値より長い場合、当該基準値より長い時間差に対応する位置登録信号のウエイトを0としてもよい。なお、上記基準値としては、例えば、周期的位置登録の周期54分を採用してもよいし、また、通信中に周期位置登録が発生する場合を考慮し、周期的位置登録の周期54分より長めに設定された時間を採用してもよいし、また、集計対象全ユーザの平均位置登録間隔を採用してもよい。このとき、第1の時間差および第2の時間差のうち一方のみが、基準値より長い場合、人口分布算出部は、基準値以下の時間差に対応する位置登録信号のウエイトを1としてもよいし、別の態様として、人口分布算出部は、基準値と当該基準値以下の時間差との差分を当該基準値により除算して得られた値を、当該基準値以下の時間差に対応する位置登録信号のウエイトとして算出してもよい。上記の抽出部は、各ユーザについての位置登録信号のうち、端末が位置登録エリアを跨いだことに起因して生成された位置登録信号を除去する前の位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号である直前の除去前位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号である直後の除去前位置登録信号とを抽出するとともに、各ユーザについての位置登録信号のうち、端末が位置登録エリアを跨いだことに起因して生成された位置登録信号を除去した後の位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号である直前の除去後位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号である直後の除去後位置登録信号とを抽出し、人口分布算出部は、各ユーザについて、直前の除去前位置登録信号の位置登録時刻と対象時刻との第1の時間差、直後の除去前位置登録信号の位置登録時刻と対象時刻との第2の時間差、直前の除去後位置登録信号の位置登録時刻と対象時刻との第3の時間差、および、直後の除去後位置登録信号の位置登録時刻と対象時刻との第4の時間差を算出し、各ユーザについて、第1の時間差と第2の時間差に基づいて、直前の除去前位置登録信号のウエイトおよび直後の除去前位置登録信号のウエイトを算出するとともに、第3の時間差と第4の時間差に基づいて、直前の除去後位置登録信号のウエイトおよび直後の除去後位置登録信号のウエイトを算出し、予め定められた按分パラメータに基づいて、各ユーザについて、直前および直後の除去前位置登録信号のウエイトと直前および直後の除去後位置登録信号のウエイトとを按分し、集計対象全ユーザについての、按分後の直前および直後の除去前位置登録信号のウエイト、並びに按分後の直前および直後の除去後位置登録信号のウエイトを、セクタごとに集計してもよい。なお、抽出部により抽出された直前の除去前位置登録信号と直後の除去前位置登録信号の両方とも、端末が位置登録エリアを跨いだことに起因して生成された位置登録信号でない場合は、直前の除去前位置登録信号と直前の除去後位置登録信号とが同一となり、直後の除去前位置登録信号と直後の除去後位置登録信号とが同一となる。従って、この場合は、直前および直後の除去後位置登録信号の抽出、これら直前および直後の除去後位置登録信号に係る第3、第4の時間差の算出、これら直前および直後の除去後位置登録信号それぞれのウエイト算出、並びに、按分パラメータに基づく按分処理は、不要となる。この場合、人口分布算出部は、各ユーザについて、直前の除去前位置登録信号の位置登録時刻と対象時刻との第1の時間差、および、直後の除去前位置登録信号の位置登録時刻と対象時刻との第2の時間差を算出し、各ユーザについて、第1の時間差と第2の時間差に基づいて、直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトを算出し、集計対象全ユーザについての直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトをセクタごとに集計することによって、対象時刻におけるセクタ単位の人口分布を算出すればよい。別の態様として、人口分布算出部は、各ユーザについて、対象時刻の直前の位置登録信号に含まれるセクタIDに対応する位置と直後の位置登録信号に含まれるセクタIDに対応する位置とを結ぶ線上に、直前の位置登録時刻から直後の位置登録時刻まで等速移動した場合の所定時間おきの通過点に相当する複数の点を生成し、生成された点が含まれるセクタに対し仮想の位置登録信号を生成し、対象とする全ユーザ以下集計対象全ユーザというについての仮想の位置登録信号から対象時刻の位置登録信号を抽出し、抽出された対象時刻の位置登録信号数をセクタごとに集計することによって、対象時刻におけるセクタ単位の人口分布を算出することができる。また、別の態様として、人口分布算出部は、各ユーザについて、対象時刻の直前の位置登録信号に含まれるセクタIDに対応する位置と直後の位置登録信号に含まれるセクタIDに対応する位置とを結ぶ線上に、直前の位置登録時刻から直後の位置登録時刻まで等速移動した場合の所定時間おきの通過点に相当する複数の点を生成し、対象時刻に対応する点が含まれるセクタに対し仮想の位置登録信号を生成し、生成された対象時刻の位置登録信号数をセクタごとに集計することにより、対象時刻におけるセクタ単位の人口分布を算出することができる。また、別の態様として、人口分布算出部は、境界がセクタ境界と対応付け可能な領域における領域内移動速度の領域間相対比率に相関する比率に基づいて、各ユーザについて、対象時刻の直前の位置登録信号に含まれるセクタIDに対応する位置と直後の位置登録信号に含まれるセクタIDに対応する位置とを結ぶ線上に、直前の位置登録時刻から直後の位置登録時刻まで移動した場合の所定時間おきの通過点に相当する複数の点を生成し、生成された点が含まれるセクタに対し仮想の位置登録信号を生成し、集計対象全ユーザについての仮想の位置登録信号から対象時刻の位置登録信号を抽出し、抽出された対象時刻の位置登録信号数をセクタごとに集計することにより、対象時刻におけるセクタ単位の人口分布を算出することができる。また、別の態様として、人口分布算出部は、境界がセクタ境界と対応付け可能な領域における領域内移動速度の領域間相対比率に相関する比率に基づいて、各ユーザについて、対象時刻の直前の位置登録信号に含まれるセクタIDに対応する位置と直後の位置登録信号に含まれるセクタIDに対応する位置とを結ぶ線上に、直前の位置登録時刻から直後の位置登録時刻まで移動した場合の所定時間おきの通過点に相当する複数の点を生成し、対象時刻に対応する点が含まれるセクタに対し仮想の位置登録信号を生成し、生成された対象時刻の位置登録信号数をセクタごとに集計することにより、対象時刻におけるセクタ単位の人口分布を算出することができる。さらに別の態様として、人口分布算出部は、各ユーザについて、直前の位置登録時刻と対象時刻との第1の時間差、および、直後の位置登録時刻と対象時刻との第2の時間差を算出し、各ユーザについて、第1の時間差と第2の時間差に基づいて、直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトを算出し、集計対象全ユーザについての直前の位置登録信号のウエイトおよび直後の位置登録信号のウエイトをセクタごとに集計することによって、対象時刻におけるセクタ単位の人口分布を算出することができる。なお、上記別の態様における人口分布算出部は、第1の時間差と第2の時間差の和により第2の時間差を除算して得られた値を、直前の位置登録信号のウエイトとして算出し、第1の時間差と第2の時間差の和により第1の時間差を除算して得られた値を、直後の位置登録信号のウエイトとして算出してもよい。なお、上記情報分析装置は、ユーザごとの属性情報を保持した属性情報保持部をさらに備え、抽出部は、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前および直後の位置登録信号の抽出を行ってもよい。別の態様として、上記情報分析装置は、ユーザごとの属性情報を保持した属性情報保持部をさらに備え、人口分布算出部は、属性情報に応じて対象とするユーザを限定し、当該限定したユーザを対象として対象時刻におけるセクタ単位の人口分布の算出を行ってもよい。上述した情報分析装置に係る発明は、情報分析装置により実行される情報分析方法に係る発明として捉え、以下のように記述することができる。本件発明の一側面に係る情報分析方法は、情報分析装置により実行される情報分析方法であって、位置登録が行われたセクタID、位置登録時刻情報およびユーザIDを含む位置登録信号を外部から受信する受信ステップと、各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する抽出ステップと、各ユーザについての、対象時刻の直前の位置登録信号に含まれる位置登録時刻情報およびセクタIDに対応する位置、対象時刻の直後の位置登録信号に含まれる位置登録時刻情報およびセクタIDに対応する位置、並びに、対象時刻に基づいて、対象時刻におけるセクタ単位の人口分布を算出する人口分布算出ステップと、を備えることを特徴とする。さらに出願人は、位置登録信号を用いた人口分布の導出に係る発明を拡張し、位置登録信号に限定されることなく、ユーザの位置を示す位置情報、当該位置情報が得られた測位時刻情報およびユーザIDを含むポイントデータを用いた人口分布の導出に対して広く適用できる以下のような発明をするに至った。本件発明の一側面に係る情報分析装置は、ユーザの位置を示す位置情報、当該位置情報が得られた測位時刻情報およびユーザIDを含むポイントデータを外部から受信する受信部と、各ユーザについてのポイントデータから、測位時刻が対象時刻の直前のポイントデータと、測位時刻が対象時刻の直後のポイントデータとを抽出する抽出部と、各ユーザについて、対象時刻の直前のポイントデータの示す位置と対象時刻の直後のポイントデータの示す位置との間を補完することにより、対象時刻におけるユーザの位置を推定し、推定した各ユーザの位置に基づいて対象時刻における所定の算出対象エリア単位の人口分布を算出する人口分布算出部と、を備えることを特徴とする。上記情報分析装置では、受信部が、ユーザの位置を示す位置情報、当該位置情報が得られた測位時刻情報およびユーザIDを含むポイントデータを外部から受信すると、抽出部は、各ユーザについてのポイントデータから、測位時刻が対象時刻の直前のポイントデータと、測位時刻が対象時刻の直後のポイントデータとを抽出し、そして、人口分布算出部は、各ユーザについて、対象時刻の直前のポイントデータの示す位置と対象時刻の直後のポイントデータの示す位置との間を補完することにより、対象時刻におけるユーザの位置を推定し、推定した各ユーザの位置に基づいて対象時刻における所定の算出対象エリア単位の人口分布を算出する。以上のように、各ユーザについての対象時刻の直前のポイントデータと直後のポイントデータを基礎として、対象時刻における所定の算出対象エリア単位の人口分布を算出することにより、巨視的な人口分布に関する精度の良いデータを簡易に且つ速やかに収集することが可能となる。また、生成される時間間隔が変動するポイントデータから、対象時刻における所定の算出対象エリア単位の人口分布を精度良く算出することができる。なお、所定の算出対象エリアとしては、セクタ、メッシュ、町丁目、市区町村など、さまざまなエリアを採用することができる。人口分布算出部は、より具体的には以下のようにして人口分布を算出することができる。例えば、人口分布算出部は、各ユーザについて、対象時刻の直前のポイントデータに含まれる位置情報の示す位置と直後のポイントデータに含まれる位置情報の示す位置とを結ぶ線上に、直前の測位時刻から直後の測位時刻まで移動した場合の所定時間おきの通過点に相当する複数の点を生成して、対象時刻に対応する点を対象時刻におけるユーザの位置として推定し、推定した各ユーザの位置に基づいて対象時刻における所定の算出対象エリア単位の人口分布を算出してもよい。所定時間おきの通過点に相当する複数の点を生成するにあたって、人口分布算出部は、上記複数の点を生成しようとする生成エリアにおけるユーザの移動特性に基づいて、該生成エリア内に生成する点の間隔を決定してもよい。上記の生成エリアにおけるユーザの移動特性は、該生成エリアにおけるユーザの平均移動速度であってもよいし、該生成エリアにおけるユーザの平均滞留時間であってもよい。なお、生成エリアは、上記の算出対象エリアと同じものであってもよいし、異なるものであってもよい。なお、上記情報分析装置は、ユーザごとの属性情報を保持した属性情報保持部をさらに備え、抽出部は、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前および直後のポイントデータを抽出してもよい。別の態様として、上記情報分析装置は、ユーザごとの属性情報を保持した属性情報保持部をさらに備え、人口分布算出部は、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として対象時刻における所定の算出対象エリア単位の人口分布を算出してもよい。なお、前述した人口分布算出部は、さらに、ユーザの端末数を人口に変換するための拡大係数をも基礎として、人口分布を算出してもよい。また、前述した情報分析装置は、対象時刻を入力するための対象時刻入力部をさらに備えてもよいし、また、算出された人口分布情報を出力する出力部をさらに備えてもよい。上記の情報分析装置に係る発明は、情報分析装置により実行される情報分析方法に係る発明として捉え、以下のように記述することができる。即ち、本件発明の一側面に係る情報分析方法は、情報分析装置により実行される情報分析方法であって、ユーザの位置を示す位置情報、当該位置情報が得られた測位時刻情報およびユーザIDを含むポイントデータを外部から受信する受信ステップと、各ユーザについてのポイントデータから、測位時刻が対象時刻の直前のポイントデータと、測位時刻が対象時刻の直後のポイントデータとを抽出する抽出ステップと、各ユーザについて、対象時刻の直前のポイントデータの示す位置と対象時刻の直後のポイントデータの示す位置との間を補完することにより、対象時刻におけるユーザの位置を推定し、推定した各ユーザの位置に基づいて対象時刻における所定の算出対象エリア単位の人口分布を算出する人口分布算出ステップと、を備えることを特徴とする。本発明によれば、各ユーザについての対象時刻の直前のポイントデータと直後のポイントデータを基礎として、対象時刻における所定の算出対象エリア単位の人口分布を算出することにより、巨視的な人口分布に関する精度の良いデータを簡易に且つ速やかに収集することが可能となる。特に、ポイントデータとして位置登録信号を用いた場合には、各ユーザについての対象時刻の直前の位置登録信号と直後の位置登録信号を基礎として、対象時刻におけるセクタ単位の人口分布を算出することにより、巨視的な人口分布に関する精度の良いデータを簡易に且つ速やかに収集することが可能となる。また、このとき、位置登録信号の発生数を基礎とすることなく、各ユーザについての対象時刻の直前の位置登録信号および直後の位置登録信号を基礎として、対象時刻におけるセクタ単位の人口分布を算出するため、生成される時間間隔が変動する位置登録信号から、対象時刻におけるセクタ単位の人口分布を精度良く算出することができる。また、位置登録信号の発生数を基礎とした場合に想定される不都合例えば、位置登録エリア境界付近で実際の人口数よりもかなり多くの位置登録信号が発生することによる不都合を防止することができる。添付図面を参照しながら本発明の実施形態を説明する。可能な場合には、同一の部分には同一の符号を付して、重複する説明を省略する。[第1実施形態][通信システムの構成]図1は、本実施形態の通信システム10のシステム構成を示す図である。図1に示すように、この通信システム10は、移動機100、BTS基地局200、RNC無線ネットワーク制御装置300、交換機400、および管理センタ500を含んで構成されている。また、この管理センタ500は、社会センサユニット501、ペタマイニングユニット502、モバイルデモグラフィーユニット503、および可視化ソリューションユニット504から構成されている。RNC300は、BTS200を介して移動機100が送信した位置登録信号を受信し、位置登録信号の信号数をカウントする。なお、標準仕様書RadioResourceControlProtocolSpecification:3GPPTS25331によれば、RNC300においてセクタ単位での信号数カウントについて規定されており、本方法はこれに準拠するものである。具体的には、例えば、RRCConnectionRequest信号でパラメタRegistrationが設定されているものを計測してもよい。或いは、更に上位の信号内容を確認してもよい。交換機400は、BTS200、RNC300を介して、移動機100が送信した位置登録信号を収集する。RNC300が、セクタ単位で位置登録信号をカウントしているのに対して、交換機400は、位置登録エリア単位で移動機100等を管理しており、移動機100が送信した位置登録信号を収集することにより位置登録エリアに登録している移動機100の数を位置登録エリア毎に把握及び記憶している。交換機400は記憶している位置登録エリアに登録している移動機100の数登録数を、所定のタイミング、または管理センタ500からの要求に応じて管理センタ500に出力する。管理センタ500は、上述したとおり、社会センサユニット501、ペタマイニングユニット502、モバイルデモグラフィーユニット503、および可視化ソリューションユニット504を含んで構成されており、各ユニットでは、移動機100が送信した位置登録信号を用いた統計処理を行う。社会センサユニット501は、各交換機400から交換機400が記憶している位置登録エリアに登録している移動機100の数等を含んだデータを収集するサーバ装置である。この社会センサユニット501は、交換機400から定期的に出力されたデータを受信し、及び社会センサユニット501において予め定められたタイミングに従って交換機400からデータを取得するように構成されている。ペタマイニングユニット502は、社会センサユニット501から受信したデータを所定のデータ形式に変換するサーバ装置である。例えば、ペタマイニングユニット502は、ユーザIDをキーに、或いはエリアごとにソーティング処理を行う。モバイルデモグラフィーユニット503は、ペタマイニングユニット502において処理されたデータに対する集計処理、すなわち各項目のカウンティング処理を行うサーバ装置である。例えば、モバイルデモグラフィーユニット503は、後述するような位置登録信号を用いた人口分布の導出等の処理を行うことができる。可視化ソリューションユニット504は、モバイルデモグラフィーユニット503において集計処理されたデータを可視可能に処理するサーバ装置である。例えば、可視化ソリューションユニット504は、集計されたデータを地図上にマッピング処理することができる。この可視化ソリューションユニット504にて処理されたデータは、企業、官公庁または個人等に提供され、店舗開発、道路交通調査、災害対策、環境対策などに利用される。なお、このように統計処理された情報は、当然にプライバシーを侵害しないように個人等は特定されないように加工されている。なお、社会センサユニット501、ペタマイニングユニット502、モバイルデモグラフィーユニット503および可視化ソリューションユニット504はいずれも、前述したようにサーバ装置により構成され、図示は省略するが、通常の情報処理装置の基本構成即ち、CPU、RAM、ROM、キーボードやマウス等の入力デバイス、外部との通信を行う通信デバイス、情報を記憶する記憶デバイス、および、ディスプレイやプリンタ等の出力デバイスを備えることは言うまでもない。図2に通信システム10の機能構成を示す。図2に示すように通信システム10は、複数のBTS200それぞれが制御するセクタに在圏する複数の移動機100と、BTS200を制御するRNC300と、交換機400と、情報分析装置600と、を含んで構成される。情報分析装置600は、前述した図1に示したモバイルデモグラフィーユニット503および可視化ソリューションユニット504に相当する。図1の社会センサユニット501およびペタマイニングユニット502に対応する機能に関しては、図2ではこれらの表記を省略している。まず、RNC300について説明する。RNC300は、RNC通信制御部302、位置登録信号受信部303、及び信号数計測部304を含んで構成されている。このうちRNC通信制御部302は、BTS200を介して移動機100と通信接続を行う部分であり、例えば、移動機100からの発信処理に基づいた通信接続処理および位置登録要求に基づいた通信接続処理を行う。本実施形態においては、移動機100は、例えば、移動機100が位置の登録を要求する場合の単位領域である位置登録エリアを跨いで移動した場合、移動機100の電源オン等により移動機100をネットワークに登録するアタッチ処理を行う場合、移動機100の電源オフ等により移動機100の登録状態を削除するデタッチ処理を行う場合、一定の周期例えば54分毎、などのタイミングで、位置登録エリアへの位置登録を求める信号である位置登録信号を送信している。更に、本実施形態では、RNC通信制御部302は、通信接続処理に用いられるInitialUEMessageを交換機400に送信することができる。なお、このInitialUEMessageは、発信または位置登録要求を示す指示情報位置登録信号、及び移動機100を一意に特定するテンポラリID等のID、を含んでいる。また、InitialUEMessageは、移動機100の位置情報を付加することも可能である。なお、テンポラリIDとは、移動機100がネットワークに接続した際に、交換機400により払い出されたID情報である。図3は、BTS200とセクタとの関係を示す図である。円形で示されている領域の中心にBTS200が位置するものであり、それを中心に複数に等分されたものがセクタである。例えば図3では、BTS200の通信エリアは最大6セクタからなるものであり、それぞれのセクタごとに、セクタを一意に識別可能であるセクタIDセクタ識別子が割り当てられており、RNC300は、該セクタIDを用いることにより移動機100が何れのセクタにいるかを、BTS200経由で把握可能となっている。なお、RNC300は、RRCコネクション要求の処理を行った際に得られる信号の遅延に基づいてセクタ内のどの位置に移動機100が位置するか、そのGAIGeographicalAreaIDを算出することもできる。そのセクタIDとセクタ内の位置に基づいて移動機100の位置を特定することもできる。位置登録信号受信部303は、移動機100が送信した位置登録信号をRNC通信制御部302経由で受信する部分である。信号数計測部304は、位置登録信号受信部303が受信した位置登録信号から、上記方法によってセクタ毎に位置登録信号があった回数の合計であるセクタ単位信号数信号数:すなわち該当するセクタにおいて移動機100が単位時間に送信した位置登録信号の数を計測する部分である。信号数計測部304は、RNC通信制御部302を介してセクタ単位信号数を交換機400へ送信する。次に、交換機400について説明する。交換機400は、交換機通信制御部401、変換部402、位置登録信号処理部404、および記憶部403を含んで構成されている。このうち交換機通信制御部401は、RNC300から送信されるInitialUEMessageを受信し、このInitialUEMessageを用いて通信接続処理を行う部分である。変換部402は、交換機通信制御部401により受信されたInitialUEMessageに含まれているテンポラリID等のIDを電話番号に変換する部分である。変換部402は、変換処理に際して、加入者プロファイル情報を記憶する加入者プロファイル情報記憶部図示せずから、テンポラリID等のIDに対応付けられている電話番号を抽出し、当該抽出した電話番号に変換する。なお、この加入者プロファイル情報記憶部は、例えば図示しないHLRHomeLocationRegisterに備えられているものであり、ここではテンポラリID等のIDと電話番号とを対応付けて管理及び記憶している。位置登録信号処理部404は、位置登録信号を移動機100からBTS200を介して受信し、受信した位置登録信号に基づいて位置登録エリアに登録している移動機100の数である登録数を計測する部分である。上述の通り本実施形態では移動機100は、移動機100が位置登録エリアを跨いで移動した場合に位置登録信号を送信している。それにより、位置登録信号処理部404は、位置登録エリア内に存在する移動機100の実数を把握できることとなる。なお、標準仕様書MobileApplicationPartspecification3GPPTS29002には、交換機400が位置登録を管理する方法について記載されている。本実施形態の交換機400における処理はこの方法に準拠したものである。記憶部403は、位置登録信号処理部404が計測した登録数、及び交換機通信制御部401を介してRNC300から受信したセクタ単位信号数を入力し、記憶する部分である。なお、変換部402で変換された電話番号とInitialUEMessageに含まれている移動機100の位置情報と当該位置情報が測位された時刻とを対応付けて記憶することもできる。なお、情報分析装置600は、例えば上記のようにして取得された移動機100の位置情報に基づいて、後述の人口分布の導出処理を行う。ただし、位置情報の取得方法は、前述した処理に限定されるものではなく、また、位置情報として、位置登録信号以外に、ユーザのGPS位置情報を対象とすることもできる。次に、情報分析装置600について説明する。情報分析装置600は、本発明に関連する構成要素として、情報分析装置通信制御部601、保持部602、抽出部603、人口分布算出部604、出力部605、及び対象時刻入力部606を含んで構成されている。このうち情報分析装置通信制御部601は、情報分析装置600と交換機400との間の通信を制御する部分であり、後述する人口分布の導出処理のために、位置登録が行われたセクタID、位置登録時刻情報およびユーザIDを含む位置登録信号を、交換機通信制御部401およびRNC通信制御部302経由で受信する。保持部602は、情報分析装置通信制御部601により受信された位置登録信号を保持する部分である。また、保持部602は、ユーザごとの属性情報例えば年齢、性別、住所などを保持してもよい。以下では、保持部602がユーザごとの属性情報を保持しているものとして説明する。対象時刻入力部606は、いつの時点の人口分布を求めるべきか、即ち、求めるべき人口分布の対象時刻をユーザが入力するための部分である。抽出部603は、後述する人口分布算出の対象となる各ユーザについての位置登録信号から、位置登録時刻が対象時刻の直前の位置登録信号と、位置登録時刻が対象時刻の直後の位置登録信号とを抽出する部分である。なお、後述する人口分布算出の対象となるユーザは、受信した位置登録信号の全ユーザであってもよいし、属性情報に応じて限定したユーザであってもよい。即ち、抽出部603は、保持部602により保持されたユーザごとの属性情報を参照することで、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前直後の位置登録信号の抽出を行ってもよい。人口分布算出部604は、後述する処理によって、対象時刻におけるセクタ単位の人口分布を算出する部分である。なお、本実施形態における人口分布算出部604は、RNC300によりカウントされたセクタ単位の位置登録信号数や、交換機400により求められた各位置登録エリアの登録数を用いることなく、後述する処理によって位置登録信号に基づいて、対象時刻におけるセクタ単位の人口分布を算出する。出力部605は、人口分布算出部604により算出された人口分布情報を出力する部分である。なお、ここでの出力は、表示出力および印刷出力を広く含むことは言うまでもない。即ち、人口分布情報は、ディスプレイ等に表示出力されてもよいし、プリンタ等から印刷出力されてもよいし、表示と印刷の両方で出力されてもよい。[情報分析装置にて実行される処理について]次に、上記のような通信システム10の中で情報分析装置600により実行される処理について、図4~図7を参照して説明する。情報分析装置600においてユーザが対象時刻入力部606から人口分布の対象時刻を入力し、所定の処理開始の操作を行うと、図4の処理が情報分析装置600において実行開始される。ただし、対象時刻の入力は必須条件ではなく、例えば対象時刻として既定の時刻を予め設定しておいてもよい。また、処理開始については、ユーザが所定の処理開始の操作を行うことは必須ではなく、定時動作などで自動的に図4の処理を実行開始してもよい。図4のステップS101では、情報分析装置通信制御部601が、位置登録が行われたセクタID、位置登録時刻情報およびユーザIDを含む位置登録信号を、交換機通信制御部401およびRNC通信制御部302経由で受信し、保持部602が、受信された位置登録信号を保持する。次のステップS102では、以下に述べる図5のユーザごとの位置登録信号の補完処理が実行される。抽出部603は、保持部602に保持された、あるユーザの位置登録信号から、対象時刻直前の位置登録信号以下直前の位置登録信号というおよび対象時刻直後の位置登録信号以下直後の位置登録信号というを抽出しステップS102A、人口分布算出部604は、直前に位置登録が行われたセクタに対応する点即ち、直前の位置登録信号に含まれたセクタIDに対応する位置と、直後に位置登録が行われたセクタに対応する点即ち、直後の位置登録信号に含まれたセクタIDに対応する位置との間に等間隔に仮想的に点を生成するステップS102B。なお、位置登録が行われたセクタに対応する点即ち、位置登録信号に含まれたセクタIDに対応する位置としては、セクタ中心位置を採用してもよいし、いわゆるPRACH-PD測位計算により取得されるポイントGAIGeographicalAreaIDを採用してもよいし、セクタの幾何学的な重心位置を採用してもよい。また、ステップS102Bでは、仮想的な点を、直線上に生成してもよいし、仮想的な地図上での道路や線路などに合せた線上に生成してもよい。ここでは、一例として、位置登録が行われたセクタに対応する点としてセクタ中心位置を採用し、仮想的な点を直線上に生成する例を説明する。図6に示すように、直前直後の位置登録信号の時間間隔がm分で、n分単位の人口即ちn分間隔の人口を求める場合、直前の位置登録信号に含まれたセクタIDに対応する位置ここではセクタ中心位置Aと直後の位置登録信号に含まれたセクタIDに対応する位置ここではセクタ中心位置Bとを結ぶ直線をm/n等分し、n分間隔で等間隔に点を仮想的に生成する。例えば、対象時刻が030で、あるユーザについての直前の位置登録信号に係る位置登録時刻が000、直後の位置登録信号に係る位置登録時刻が033で、1分単位の人口即ち1分間隔の人口を求める場合、直前の位置登録信号に係るセクタ中心位置Aと、直後の位置登録信号に係るセクタ中心位置Bとを結ぶ直線を33等分し、1分間隔で等間隔に点を仮想的に生成する。そして、人口分布算出部604は、生成した点が含まれるセクタに対し、仮想の位置登録信号を生成するステップS102C。以上のような図5の処理は、各ユーザについて実行される。なお、図5のステップS102Aにおいて抽出部603は、保持部602により保持されたユーザごとの属性情報を参照することで、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前直後の位置登録信号の抽出を行ってもよい。この場合、以降の処理対象となる位置登録信号を絞り込むことができるため、処理負荷の軽減および処理の高速化を図ることができる。集計対象全ユーザについて図5の処理が実行完了すると、図4へ戻り、次のステップS103では、人口分布算出部604は、集計対象全ユーザを対象にして、図5の処理で生成された位置登録信号から、対象時刻の位置登録信号を抽出する。なお、上述した例では、位置登録時刻が対象時刻030の点が含まれるセクタに対し生成された位置登録信号が抽出される。そして、次のステップS104で人口分布算出部604は、セクタIDごとに集計対象全ユーザの位置登録信号数を集計する。これにより、対象時刻におけるセクタIDごとの位置登録信号数が得られ、以下に述べる処理により、対象時刻におけるセクタ単位の人口分布が求められる。ここで、対象時刻におけるセクタIDごとの集計対象全ユーザの位置登録信号数即ち、セクタIDごとのユーザ数から、対象時刻におけるセクタ単位の人口を求める処理以下拡大推計処理というについて説明する。以下、拡大推計処理の一態様として、ユーザ数を人口に変換するための係数以下拡大係数というを用いて人口を求める処理を説明する。上記の拡大係数は、例えばユーザの属性ごとに予め求められており、人口分布算出部604は、ユーザの属性例えば性別、年齢層、住所などごとの拡大係数を予め記憶している。ステップS104で、人口分布算出部604は、例えば、集計対象ユーザを属性別に分類し、分類後の属性ごとに、対象時刻におけるセクタIDごとの位置登録信号数を集計し、その集計値に、当該属性に対応する拡大係数を乗算することで、対象時刻におけるセクタIDごとの人口を属性ごとに求める。このとき、セクタ単位の属性ごとの人口が得られる。さらに、人口分布算出部604は、属性ごとの人口をセクタIDごとに集計することで、対象時刻におけるセクタ単位の人口を求めてもよい。また、別の導出方法として、ステップS104で人口分布算出部604は、保持部602に保持されたユーザごとの属性情報を参照し、集計対象ユーザの位置登録信号それぞれに対応する拡大係数、即ち、集計対象ユーザそれぞれの属性に対応する拡大係数を読み出し、前述した例のように集計対象ユーザを属性別に分類することなく、集計対象全ユーザについて、セクタIDごとに拡大係数を合計することで、対象時刻におけるセクタ単位の人口を求めてもよい。なお、拡大係数は、ユーザの属性に応じて設定する以外にも、例えば、時間帯や場所に応じて設定してもよい。また、拡大係数は、予め記憶しておく以外に、以下のようにして求めてもよい。拡大係数は、一例として、在圏率と端末の普及率との積即ち、人口に対する在圏数の比率の逆数を用いることができる。ここで在圏率とは、契約台数に対する在圏数の比率を意味し、端末の普及率とは人口に対する端末の契約台数の比率を意味する。このような拡大係数は、上記の属性、時間帯、場所といった単位で求めることが望ましいが、必須ではない。拡大係数を求める単位としては、一例として、住所の都道府県ごと、5才刻み年齢層ごと、男女ごと、時間帯として1時間ごとなどを採用してもよい。次のステップS105では、出力部605が、上記のようにして求められたセクタ単位の人口分布を二次元的に表すセクタ単位の人口分布図を出力する。図7にセクタ単位の人口分布図の一例を示す。この図7では、A~LはセクタIDを表し、各セクタIDの下の数字は当該セクタの人口を表す。また、別の出力例として、出力部605は、図31に示すように、各セクタについて推定された人口図31では推定人口と表記や、各セクタにおける性別年齢層住所といった属性ごとの人口を、表形式で出力することができる。また、複数の属性を組み合わせた条件例えば性別と住所を組み合わせた東京都在住の女性といった条件を採用してもよい。なお、ここでの出力には、表示出力、音声出力、印刷出力等さまざまな出力態様が含まれる。以上説明した第1実施形態により、各ユーザについての対象時刻の直前の位置登録信号と直後の位置登録信号を基礎として、対象時刻におけるセクタ単位の人口分布を算出することにより、巨視的な人口分布に関する精度の良いデータを簡易に且つ速やかに収集することが可能となる。また、このとき、位置登録信号の発生数を基礎とすることなく、各ユーザについての対象時刻の直前の位置登録信号および直後の位置登録信号を基礎として、対象時刻におけるセクタ単位の人口分布を算出するため、生成される時間間隔が変動する位置登録信号から、対象時刻におけるセクタ単位の人口分布を精度良く算出することができる。また、位置登録信号の発生数を基礎とした場合に想定される不都合例えば、位置登録エリア境界付近で実際の人口数よりもかなり多くの位置登録信号が発生することによる不都合を防止することができる。なお、第1実施形態における図5の補完処理のステップS102Cでは、等間隔に生成した複数の点それぞれが含まれるセクタに対して仮想の位置登録信号を生成した例を説明したが、この変形例として、等間隔に生成した複数の点のうち、対象時刻に対応する点が含まれるセクタに対してだけ、仮想の位置登録信号を生成してもよい。つまり、図14に示すように、図5のステップS102Cに代わるステップS102Dにおいて、対象時刻に対応する点が含まれるセクタに対して仮想の位置登録信号を生成すればよい。この場合、対象時刻に対応する点が含まれるセクタに対する1つの仮想の位置登録信号のみが生成されるため、図4のステップS103の処理対象時刻の位置登録信号の抽出が不要となり、さらなる処理負荷の軽減および処理の高速化を図ることができる。[第2実施形態]前述した第1実施形態の処理は、対象時刻の直前直後の位置登録信号に対応する点の間に仮想的に生成した複数の点を利用してセクタ単位の人口分布を求めるのに対し、以下に述べる第2実施形態の処理は、対象時刻の直前直後の位置登録信号のウエイトを利用してセクタ単位の人口分布を求める。なお、第2実施形態では、図1の通信システム10のシステム構成、および図2の通信システム10の機能ブロック構成は第1実施形態と同様であるので、重複した説明は省略する。人口分布算出部604により実行される処理内容が異なるので、かかる処理内容を含む第2実施形態の処理について、図8~図11を参照して、以下に説明する。情報分析装置600においてユーザが対象時刻入力部606から人口分布の対象時刻を入力し、所定の処理開始の操作を行うと、図8の処理が情報分析装置600において実行開始される。ただし、対象時刻の入力は必須条件ではなく、例えば対象時刻として既定の時刻を予め設定しておいてもよい。また、処理開始については、ユーザが所定の処理開始の操作を行うことは必須ではなく、定時動作などで自動的に図8の処理を実行開始してもよい。図8のステップS201では、情報分析装置通信制御部601が、位置登録が行われたセクタID、位置登録時刻情報およびユーザIDを含む位置登録信号を、交換機通信制御部401およびRNC通信制御部302経由で受信し、保持部602が、受信された位置登録信号を保持する。次のステップS202では、以下に述べる図9のユーザごとの位置登録信号のウエイト算出処理が実行される。抽出部603は、保持部602に保持された、あるユーザの位置登録信号から、対象時刻直前の位置登録信号以下直前の位置登録信号というおよび対象時刻直後の位置登録信号以下直後の位置登録信号というを抽出しステップS202A、人口分布算出部604は、直前の位置登録時刻と対象時刻との第1の時間差、および、直後の位置登録時刻と対象時刻との第2の時間差を算出するステップS202B。そして、人口分布算出部604は、第1の時間差と第2の時間差の和により第2の時間差を除算して得られた値を、直前の位置登録信号のウエイトとして算出し、第1の時間差と第2の時間差の和により第1の時間差を除算して得られた値を、直後の位置登録信号のウエイトとして算出するステップS202C。このとき、図10に示すように、直前の位置登録時刻と対象時刻との第1の時間差がm分、直後の位置登録時刻と対象時刻との第2の時間差がn分である場合、直前の位置登録信号Aのウエイトはn/mn、直後の位置登録信号Bのウエイトはm/mnとなる。例えば、対象時刻が030で、あるユーザについての直前の位置登録信号に係る位置登録時刻が000、直後の位置登録信号に係る位置登録時刻が040の場合、第1の時間差mが30分、第2の時間差nが10分となるため、直前の位置登録信号Aのウエイトは10/3010=025、直後の位置登録信号Bのウエイトは30/3010=075となる。なお、直前直後の位置登録信号のウエイトの算出方法は、上記の方法に限定されるものではなく、第1、第2の時間差を用いた別の数式により位置登録信号のウエイトを算出してもよい。なお、ウエイトにおける小数点以下の数値については、例えば、小数点以下3けた目以降が存在する場合は同3けた目を四捨五入して小数点以下2けた目までの数値とするなど、予め定めた手法に従ってウエイトを求めればよい。以上のような図9の処理は、集計対象全ユーザについて実行される。なお、図9のステップS202Aにおいて抽出部603は、保持部602により保持されたユーザごとの属性情報を参照することで、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前直後の位置登録信号の抽出を行ってもよい。この場合、以降の処理対象となる位置登録信号を絞り込むことができるため、処理負荷の軽減および処理の高速化を図ることができる。集計対象全ユーザについて図9の処理が実行完了すると、図8へ戻り、次のステップS203では、人口分布算出部604は、セクタIDごとに集計対象全ユーザの位置登録信号のウエイトを集計する。これにより、対象時刻におけるセクタIDごとの位置登録信号のウエイトの合計値が得られ、第1実施形態で述べた拡大推計処理によって、ウエイトの合計値から対象時刻におけるセクタ単位の人口分布が求められる。次のステップS204では、出力部605が、上記のようにして求められたセクタ単位の人口分布を二次元的に表すセクタ単位の人口分布図を出力する。図11にセクタ単位の人口分布図の一例を示す。この図11では、A~LはセクタIDを表し、各セクタIDの下の数字は当該セクタの人口を表す。図11の例は、位置登録信号のウエイトを小数点以下1けた目までの数値とし、その後、セクタIDごとに集計対象全ユーザの位置登録信号のウエイトを集計した例を示している。以上説明した第2実施形態により、各ユーザについての対象時刻の直前の位置登録信号と直後の位置登録信号を基礎として、対象時刻におけるセクタ単位の人口分布を算出することにより、巨視的な人口分布に関する精度の良いデータを簡易に且つ速やかに収集することが可能となる。また、このとき、位置登録信号の発生数を基礎とすることなく、各ユーザについての対象時刻の直前の位置登録信号および直後の位置登録信号を基礎として、対象時刻におけるセクタ単位の人口分布を算出するため、位置登録信号の発生数を基礎とした場合に想定される不都合例えば、位置登録エリア境界付近で実際の人口数よりもかなり多くの位置登録信号が発生することによる不都合を防止することができる。なお、上述した第1、第2実施形態では、セクタ単位の人口分布を算出する例を示したが、さらに、同じ位置登録エリアに属するセクタの人口を位置登録エリアごとに集計する処理を行うことで、位置登録エリア単位の人口分布を算出することができる。例えば、図11におけるセクタA~Eが属する位置登録エリアの人口と、セクタF~Lが属する位置登録エリアの人口とを算出し、図12に示すように、各位置登録エリアに、人口数に応じた表記色分け、斜線、網掛けなどを施すことにより、セクタ単位の人口分布とともに位置登録エリア単位の人口分布を表す図を出力することができる。また、上述した第1、第2実施形態では、ユーザごとの属性情報例えば年齢、性別、住所などが保持部602により保持される例を示したが、これに限定されるものではなく、図15に示すように、情報分析装置600が保持部602とは別に、ユーザごとの属性情報を保持する属性情報保持部607を備える構成を採用してもよい。また、上述した第1、第2実施形態では抽出部603が、属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前直後の位置登録信号の抽出を行う例を示したが、これに限定されるものではなく、人口分布算出部605が、例えば図4のステップS103、S104や図8のステップS203において属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として人口分布算出を行ってもよい。このように位置登録信号の抽出時点でなく人口分布算出の際に、属性情報に応じて、対象とするユーザを限定する場合、例えば年齢別や男女別など、さまざまな条件で人口分布算出を行うことができる。もちろん、位置登録信号の抽出と人口分布算出の両方で、属性情報に応じた対象ユーザの限定を行ってもよい。[第3実施形態]以下の第3実施形態では、第1実施形態の変形態様として、人口分布算出部604が、対象時刻の直前の位置登録信号に含まれたセクタIDに対応する位置と、直後の位置登録信号に含まれたセクタIDに対応する位置との間に、所定時間おきの通過点に相当する複数の仮想的な点を非等間隔に生成して、仮想の位置登録信号を生成する例を説明する。具体的には、セクタ内移動速度のセクタ間相対比率に相関する比率を予め求めておき、その比率に基づいて非等間隔に仮想的な点を生成して仮想の位置登録信号を生成する実施形態である。上記のセクタ内移動速度のセクタ間相対比率に相関する比率としては、後述するセクタ内移動速度のセクタ間相対比率そのものを用いてもよいし、他の比率例えば、後述するセクタ内滞留時間のセクタ間相対比率を用いてもよい。なお、上記の相対比率の単位は、セクタ単位であることは必須ではなく、境界がセクタ境界と対応付け可能な領域の単位例えばメッシュ、町丁目等であれば用いることができる。ただし、以下では説明を容易にするために、相対比率の単位がセクタ単位である例を説明する。セクタ内移動速度のセクタ間相対比率を用いた例まず、セクタ内移動速度のセクタ間相対比率に相関する比率として、セクタ内移動速度のセクタ間相対比率を用いた例を説明する。以下では、重複した説明を避けるべく、第1実施形態との差異について重点的に説明する。第3実施形態では、第1実施形態における図5の位置登録信号の補完処理が図16に示す処理となる。即ち、図16のステップS102Aで対象時刻の直前の位置登録信号および直後の位置登録信号が抽出された後、人口分布算出部604は、ステップS102BXにおいて、直前の位置登録信号に含まれたセクタIDに対応する位置と、直後の位置登録信号に含まれたセクタIDに対応する位置との間に、予め求められたセクタ内移動速度のセクタ間相対比率に基づいて、所定時間おきの通過点に相当する複数の仮想的な点を生成する。図17に示すように、対象時刻の直前の位置登録信号に含まれたセクタIDに対応する位置Aがセクタ1にあり、直後の位置登録信号に含まれたセクタIDに対応する位置Bがセクタ3にあり、セクタ1と3の間にセクタ2が存在する状況で、予め求められたセクタ内移動速度のセクタ間相対比率セクタ1:セクタ2:セクタ3が1:2:3であったとする。ここで、位置A、Bを結ぶ線とセクタ1、2の境界線との交差点を点Xとし、位置A、Bを結ぶ線とセクタ2、3の境界線との交差点を点Yとした場合、位置A、Bを結ぶ線上に、所定時間おき例えば1分おきの通過点に相当する複数の仮想的な点を生成する際に、点の間隔の相対比率AX間:XY間:YB間が、上記のセクタ内移動速度のセクタ間相対比率、即ち、1:2:3となるようにする。つまり、セクタ内移動速度がセクタ2、3に比べ相対的に低いセクタ1ここではAX間では、点の間隔がセクタ2、3に比べ相対的に狭くなる。一方、セクタ内移動速度がセクタ1、2に比べ相対的に高いセクタ3ここではYB間では、点の間隔がセクタ1、2に比べ相対的に広くなる。このようにしてステップS102BXで位置A、Bを結ぶ線上に仮想的に点を生成した後、人口分布算出部604は、生成した点が含まれるセクタに対し、仮想の位置登録信号を生成するステップS102C。以上のような図16の処理により、セクタ内移動速度のセクタ間相対比率に基づいて、位置A、Bを結ぶ線上に非等間隔に仮想的な点を生成して仮想の位置登録信号を生成することができる。なお、図16の補完処理のステップS102Cの変形例として、非等間隔に生成した複数の点のうち、対象時刻に対応する点が含まれるセクタに対してだけ、仮想の位置登録信号を生成してもよい。つまり、図18に示すように、図16のステップS102Cに代わるステップS102Dにおいて、対象時刻に対応する点が含まれるセクタに対して仮想の位置登録信号を生成すればよい。この場合、対象時刻に対応する点が含まれるセクタに対する1つの仮想の位置登録信号のみが生成されるため、さらなる処理負荷の軽減および処理の高速化を図ることができる。ここで、セクタ内移動速度のセクタ間相対比率の求め方について概説する。一例としては、セクタごとに、まずユーザごとの移動速度を後述の方法により求め、次に全ユーザの移動速度を平均することで当該セクタのセクタ内平均移動速度を求める。そして、セクタごとのセクタ内平均移動速度同士の比率から、セクタ内移動速度のセクタ間相対比率を求めることができる。なお、セクタ内移動速度のセクタ間相対比率を予め求める際には、位置登録信号より粒度の細かい位置情報例えばユーザのGPS位置情報等を用いることが望ましい。また、セクタ内移動速度のセクタ間相対比率は、特定の時間帯ごとに例えば朝昼晩深夜などの所定時間帯ごとに求めてもよいし、曜日ごとに求めてもよいし、平日の平均や休日の平均を求めてもよい。ユーザごとの移動速度を求める方法としては、例えば以下のような第1、第2の方法を採用できる。なお、ユーザごとの移動速度を求める対象のセクタを、ここでは対象セクタと称する。第1の方法では、対象セクタ内で最初に観測された位置情報および最後に観測された位置情報より、それら2つの位置の間の距離および観測時刻の差分を求め、求めた距離を観測時刻の差分により割り算することで、移動速度を求める。例えば図19aに示すように、対象セクタ内で最初に観測された位置情報位置Dおよび観測時刻1200と、最後に観測された位置情報位置Eおよび観測時刻1205より、DE間の距離150mおよび観測時刻の差分5分を求め、距離150mを観測時刻の差分5分により割り算することで、移動速度18km/時を求めることができる。第2の方法では、対象セクタ内で最初に観測された位置情報および最後に観測された位置情報に加え、対象セクタに入る直前の位置情報および対象セクタから出た直後の位置情報を用いる。例えば図19bに示すように、対象セクタに入る直前の位置情報位置Cおよび観測時刻1155および対象セクタ内で最初に観測された位置情報位置Dおよび観測時刻1200より、CD間を等速直線移動したと仮定した場合の対象セクタに入った時刻即ち、境界上の点Pを通過した時刻ここでは仮に1158とすると、境界上の点Pの位置を求める。同様に、対象セクタ内で最後に観測された位置情報位置Eおよび観測時刻1205および対象セクタから出た直後の位置情報位置Fおよび観測時刻1210より、EF間を等速直線移動したと仮定した場合の対象セクタから出た時刻即ち、境界上の点Qを通過した時刻ここでは仮に1206とすると、境界上の点Qの位置を求める。そして、点P、Qの位置よりPQ間の距離180mを求めるとともに、対象セクタに入った時刻と対象セクタから出た時刻との差分8分を求め、距離180mを観測時刻の差分8分により割り算することで、移動速度135km/時を求めることができる。以上のようにしてユーザごとの移動速度を求めることができ、さらに、全ユーザの移動速度を平均することでセクタ内平均移動速度を求め、セクタごとのセクタ内平均移動速度同士の比率から、セクタ内移動速度のセクタ間相対比率を求めることができる。セクタ内滞留時間のセクタ間相対比率を用いた例次に、セクタ内移動速度のセクタ間相対比率に相関する比率として、セクタ内滞留時間のセクタ間相対比率を用いた例を説明する。この場合、位置登録信号の補完処理は、前述した図16、図18の処理とほぼ同様であるため、重複した説明は省略する。そこで、図16、図18のステップS102BXにおける点の生成方法の詳細について、以下、説明する。図20に示すように、対象時刻の直前の位置登録信号に含まれたセクタIDに対応する位置Aがセクタ1にあり、直後の位置登録信号に含まれたセクタIDに対応する位置Bがセクタ3にあり、セクタ1と3の間にセクタ2が存在する状況で、予め求められたセクタ内滞留時間のセクタ間相対比率セクタ1:セクタ2:セクタ3が2:3:4であったとする。ここで、位置A、Bを結ぶ線とセクタ1、2の境界線との交差点を点Xとし、位置A、Bを結ぶ線とセクタ2、3の境界線との交差点を点Yとした場合、位置A、Bを結ぶ線上に、所定時間おき例えば1分おきの通過点に相当する複数の仮想的な点を生成する際に、点の間隔の相対比率AX間:XY間:YB間が、上記のセクタ内滞留時間のセクタ間相対比率の逆数の比率、即ち、::=6:4:3となるようにする。つまり、セクタ内滞留時間がセクタ2、3に比べ相対的に短いセクタ1ここではAX間では、点の間隔がセクタ2、3に比べ相対的に広くなる。一方、セクタ内滞留時間がセクタ1、2に比べ相対的に長いセクタ3ここではYB間では、点の間隔がセクタ1、2に比べ相対的に狭くなる。このようにして図16のステップS102BXで位置A、Bを結ぶ線上に複数の仮想的な点を生成した後、人口分布算出部604は、生成した点が含まれるセクタに対し、仮想の位置登録信号を生成する図16のステップS102C。又は、図18のステップS102BXで位置A、Bを結ぶ線上に複数の仮想的な点を生成した後、人口分布算出部604は、対象時刻に対応する点が含まれるセクタに対して仮想の位置登録信号を生成する図18のステップS102D。これにより、セクタ内滞留時間のセクタ間相対比率に基づいて、位置A、Bを結ぶ線上に非等間隔に仮想的な点を生成して仮想の位置登録信号を生成することができる。ここで、セクタ内滞留時間のセクタ間相対比率の求め方について概説する。一例としては、セクタごとに、まずユーザごとの滞留時間を後述の方法により求め、次に全ユーザの滞留時間を平均することで当該セクタのセクタ内平均滞留時間を求める。そして、セクタごとのセクタ内平均滞留時間同士の比率から、セクタ内滞留時間のセクタ間相対比率を求めることができる。なお、セクタ内滞留時間のセクタ間相対比率を予め求める際には、位置登録信号より粒度の細かい位置情報例えばユーザのGPS位置情報等を用いることが望ましい。また、セクタ内滞留時間のセクタ間相対比率は、特定の時間帯ごとに例えば朝昼晩深夜などの所定時間帯ごとに求めてもよいし、曜日ごとに求めてもよいし、平日の平均や休日の平均を求めてもよい。ユーザごとの滞留時間を求める方法としては、例えば以下のような第1、第2の方法を採用できる。なお、ユーザごとの滞留時間を求める対象のセクタを、ここでは対象セクタと称する。第1の方法では、対象セクタ内で最初に観測された位置情報および最後に観測された位置情報より、これら2つの観測時刻の差分を対象セクタにおける滞留時間として求める。例えば図21aに示すように、対象セクタ内で最初に観測された位置情報位置Dおよび観測時刻1200と、最後に観測された位置情報位置Eおよび観測時刻1205より、これら2つの観測時刻の差分5分を対象セクタにおける滞留時間として求めることができる。第2の方法では、対象セクタ内で最初に観測された位置情報および最後に観測された位置情報に加え、対象セクタに入る直前の位置情報および対象セクタから出た直後の位置情報を用いる。例えば図21bに示すように、対象セクタに入る直前の位置情報位置Cおよび観測時刻1155および対象セクタ内で最初に観測された位置情報位置Dおよび観測時刻1200より、CD間を等速直線移動したと仮定した場合の対象セクタに入った時刻即ち、境界上の点Pを通過した時刻ここでは仮に1158とするを求める。同様に、対象セクタ内で最後に観測された位置情報位置Eおよび観測時刻1205および対象セクタから出た直後の位置情報位置Fおよび観測時刻1210より、EF間を等速直線移動したと仮定した場合の対象セクタから出た時刻即ち、境界上の点Qを通過した時刻ここでは仮に1206とするを求める。そして、対象セクタに入った時刻と対象セクタから出た時刻との差分8分を、対象セクタにおける滞留時間として求めることができる。以上のようにしてユーザごとの滞留時間を求めることができ、さらに、全ユーザの滞留時間を平均することでセクタ内平均滞留時間を求め、セクタごとのセクタ内平均滞留時間同士の比率から、セクタ内滞留時間のセクタ間相対比率を求めることができる。以上説明した第3実施形態によれば、予め求められたセクタ内移動速度のセクタ間相対比率に相関する比率に基づいて、所定時間おきの通過点に相当する複数の仮想的な点を非等間隔に生成して、仮想の位置登録信号を生成することができる。なお、セクタ内移動速度のセクタ間相対比率に相関する比率としてセクタ内移動速度のセクタ間相対比率そのものを用いる場合はセクタの大きさによる影響を受けずに精度良く処理を行うことができるという利点がある。一方、セクタ内移動速度のセクタ間相対比率に相関する比率としてセクタ内滞留時間のセクタ間相対比率を用いる場合は、セクタ内移動速度よりも小さい処理負荷で簡易にセクタ内滞留時間を求めることができるため、全体的な処理負荷の削減を図れるという利点がある。なお、第1~第3実施形態では、位置登録情報を用いた人口分布の導出について説明したが、本発明は、位置登録情報以外の位置情報例えばGPS位置情報やPRACH-PD座標情報等を用いた人口分布の導出にも適用可能である。具体的には、ユーザの位置を示す位置情報、当該位置情報が得られた測位時刻情報およびユーザIDを含むポイントデータが得られる環境であれば、本発明は以下のように適用可能である。例えば、図15の情報分析装置600では、情報分析装置通信制御部601が、ユーザの位置を示す位置情報、当該位置情報が得られた測位時刻情報およびユーザIDを含むポイントデータを外部から受信し、抽出部603が、各ユーザについてのポイントデータから、測位時刻が対象時刻の直前のポイントデータと測位時刻が対象時刻の直後のポイントデータとを抽出し、そして、人口分布算出部604が、各ユーザについて、対象時刻の直前のポイントデータの示す位置と対象時刻の直後のポイントデータの示す位置との間を補完することにより、対象時刻におけるユーザの位置を推定し、推定した各ユーザの位置に基づいて、第1実施形態で述べた拡大推計処理により、対象時刻における所定の算出対象エリア単位の人口分布を算出することができる。このように、各ユーザについての対象時刻の直前のポイントデータと直後のポイントデータを基礎として、対象時刻における所定の算出対象エリア単位の人口分布を算出することにより、巨視的な人口分布に関する精度の良いデータを簡易に且つ速やかに収集することが可能となる。なお、所定の算出対象エリアとしては、セクタに限定されるものではなく、メッシュ、町丁目、市区町村など、さまざまなエリアを採用することができる。このとき、人口分布算出部604は、例えば、各ユーザについて、対象時刻の直前のポイントデータに含まれる位置情報の示す位置と直後のポイントデータに含まれる位置情報の示す位置とを結ぶ線上に、直前の測位時刻から直後の測位時刻まで移動した場合の所定時間おきの通過点に相当する複数の点を生成して、対象時刻に対応する点を対象時刻におけるユーザの位置として推定し、推定した各ユーザの位置に基づいて、第1実施形態で述べた拡大推計処理により、対象時刻における所定の算出対象エリア単位の人口分布を算出してもよい。なお、所定時間おきの通過点に相当する複数の点を生成するにあたって、人口分布算出部604は、上記複数の点を生成しようとする生成エリアにおけるユーザの移動特性例えば該生成エリアにおけるユーザの平均移動速度や平均滞留時間に基づいて、該生成エリア内に生成する点の間隔を決定してもよい。なお、生成エリアは、上記の算出対象エリアと同じものであってもよいし、異なるものであってもよい。また、図15の情報分析装置600における抽出部603は、属性情報保持部607に保持されたユーザごとの属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として直前および直後のポイントデータを抽出してもよい。もしくは、図15の情報分析装置600における人口分布算出部604は、属性情報保持部607に保持されたユーザごとの属性情報に応じて、対象とするユーザを限定し、当該限定したユーザを対象として対象時刻における所定の算出対象エリア単位の人口分布を算出してもよい。[第4実施形態]以下の第4、第5実施形態は、前述した第2実施形態の変形例に相当し、対象時刻の直前直後の位置登録信号のウエイトを利用して人口分布を求める処理に関する別の実施形態をそれぞれ説明する。既に述べたように、移動機100は、例えば、移動機100が位置の登録を要求する場合の単位領域である位置登録エリアを跨いで移動したとき、移動機100の電源オン等により移動機100をネットワークに登録するアタッチ処理を行うとき、移動機100の電源オフ等により移動機100の登録状態を削除するデタッチ処理を行うとき、周期的位置登録の一定の周期例えば54分毎、などのタイミングで、位置登録エリアへの位置登録を求める位置登録信号を生成する。そのため、移動機100が電源オフや圏外にいる等に起因した通信不能状態でない限り、移動機100は、少なくとも上記周期的位置登録の一定の周期以下周期tという毎に位置登録信号を送信する。この点に着目し、第4、第5実施形態では、対象時刻の直前直後の位置登録信号のうち、位置登録時刻即ち、位置登録信号の生成時刻と対象時刻との時間差が予め定められた基準値例えばここでは周期tよりも長い位置登録信号については、位置登録時刻と対象時刻との時間差が異常値とみなし、時間差が異常値とみなされた位置登録信号以下時間差が異常値である位置登録信号というについてのウエイトを0とする。なお、上記基準値としては、周期的位置登録の周期54分を採用する以外に、通信中に周期位置登録が発生する場合を考慮し、周期的位置登録の周期54分より長めに設定された時間を採用してもよいし、また、集計対象全ユーザの平均位置登録間隔を採用してもよい。以下では一例として、基準値を周期tとした例を説明する。以下の第4実施形態では、位置登録時刻と対象時刻との時間差が上記基準値周期t以下である位置登録信号以下、時間差が正常値である位置登録信号というについては、ウエイトを1とし、第5実施形態では、時間差が正常値である位置登録信号については、周期tと当該時間差との差分に応じてウエイトを算出する詳細は後述する。まず、第4実施形態について説明する。第4実施形態では、図1の通信システム10のシステム構成、および図2の通信システム10の機能ブロック構成は第1実施形態と同様であるので、重複した説明は省略する。人口分布算出部604により実行される処理の流れは、第2実施形態で説明した図8、図9の処理の流れと同様であるが、図9のステップS202Cにおけるウエイト算出処理が、第2実施形態とは異なるため、第4実施形態に係るウエイト算出処理について、図22、図23を参照して、以下に説明する。なお、図22はウエイト算出処理のフロー図の一例を示すが、ウエイト算出処理のフローは図22に示すフローに限定されるものではない。図22に示すように、人口分布算出部604は、図9のステップS202Bで算出された直前の位置登録時刻と対象時刻との第1の時間差mおよび直後の位置登録時刻と対象時刻との第2の時間差nの少なくとも一方が周期tより長いか否かを判断する図22のステップS301。ここで、第1の時間差m、第2の時間差nともに周期t以下であれば、直前直後の位置登録信号はともに時間差が正常値である位置登録信号と判断できるので、人口分布算出部604は、第2実施形態と同様に、第1の時間差mと第2の時間差nの和により第2の時間差nを除算して得られた値を、直前の位置登録信号のウエイトとして算出し、第1の時間差mと第2の時間差nの和により第1の時間差mを除算して得られた値を、直後の位置登録信号のウエイトとして算出するステップS302。一方、ステップS301にて第1の時間差mおよび第2の時間差nの少なくとも一方が周期tより長い場合は、周期tより長い時間差に対応する位置登録信号のウエイトを0とするステップS303。次に、第1の時間差m、第2の時間差nがともに周期tより長いか否かを判断しステップS304、第1の時間差m、第2の時間差nがともに周期tより長ければ、図22の処理を終了する。この場合、ステップS303の処理により、直前直後の位置登録信号ともにウエイトは0とされている。一方、ステップS304にて一方の時間差が周期t以下であれば、周期t以下の時間差に対応する位置登録信号即ち、時間差が正常値である位置登録信号のウエイトを1とするステップS305A。以上のようにして図22の処理により、直前直後の位置登録信号それぞれについてウエイトが設定される。ここで、図23を用いて具体例を説明する。図23に示すように、ユーザαに関する位置登録信号として、対象時刻1400に対し、直前の位置登録信号位置登録時刻:1330、在圏セクタ:エリアAと直後の位置登録信号位置登録時刻:1500、在圏セクタ:エリアBとが抽出されたケースを想定する。周期的位置登録の周期tは54分とする。このケースでは、直前の位置登録時刻と対象時刻との第1の時間差mは30分、直後の位置登録時刻と対象時刻との第2の時間差nは60分となり、第2の時間差nが周期t54分より長いので、直後の位置登録信号についてのウエイトは0となる。一方、第1の時間差mは周期t54分以下なので、直前の位置登録信号についてのウエイトは1となる。そのため、ユーザαは対象時刻1400に、直前の位置登録信号に係る在圏セクタであるエリアAに1人分のウエイトで存在するという算出結果が得られる。以上の第4実施形態により、対象時刻の直前直後の位置登録信号のうち、位置登録時刻と対象時刻との時間差が周期的位置登録の周期tよりも長い位置登録信号について、時間差が異常値である位置登録信号とみなし、ウエイトを0とするため、時間差が異常値である位置登録信号を除外した上で、より精度の良いウエイト算出処理を行うことができる。[第5実施形態]第5実施形態も、第2実施形態における図9のステップS202Cウエイト算出処理の別の実施形態であるが、第5実施形態では、位置登録時刻と対象時刻との時間差が予め定められた基準値ここでは第4実施形態と同様に周期tとする以下の位置登録信号即ち、時間差が正常値である位置登録信号については、周期tと上記時間差との差分に応じてウエイトを算出する。第5実施形態では、図1の通信システム10のシステム構成、および図2の通信システム10の機能ブロック構成は第1実施形態と同様であるので、重複した説明は省略する。人口分布算出部604により実行される処理の流れは、第2実施形態で説明した図8、図9の処理の流れと同様であるが、図9のステップS202Cにおけるウエイト算出処理が、第2実施形態とは異なるため、第5実施形態に係るウエイト算出処理について、図23、図24を参照して、以下に説明する。なお、図24はウエイト算出処理のフロー図の一例を示すが、ウエイト算出処理のフローは図24に示すフローに限定されるものではない。図24に示すように、人口分布算出部604は、図9のステップS202Bで算出された直前の位置登録時刻と対象時刻との第1の時間差mおよび直後の位置登録時刻と対象時刻との第2の時間差nの少なくとも一方が周期tより長いか否かを判断する図24のステップS301。ここで、第1の時間差m、第2の時間差nともに周期t以下であれば、直前直後の位置登録信号はともに時間差が正常値である位置登録信号と判断できるので、人口分布算出部604は、第2実施形態と同様に、第1の時間差mと第2の時間差nの和により第2の時間差nを除算して得られた値を、直前の位置登録信号のウエイトとして算出し、第1の時間差mと第2の時間差nの和により第1の時間差mを除算して得られた値を、直後の位置登録信号のウエイトとして算出するステップS302。一方、ステップS301にて第1の時間差mおよび第2の時間差nの少なくとも一方が周期tより長い場合は、周期tより長い時間差に対応する位置登録信号のウエイトを0とするステップS303。次に、第1の時間差m、第2の時間差nがともに周期tより長いか否かを判断しステップS304、第1の時間差m、第2の時間差nがともに周期tより長ければ、図24の処理を終了する。このとき、ステップS303の処理により、直前直後の位置登録信号ともにウエイトは0とされている。一方、ステップS304にて一方の時間差が周期t以下であれば、周期t以下の時間差に対応する位置登録信号即ち、時間差が正常値である位置登録信号のウエイトを、以下のように、周期tと当該時間差との差分に応じて算出するステップS305B。ここでは例えば、以下の式1を用いて、周期tと当該時間差との差分に応じてウエイトを算出する。ウエイト=周期t-当該時間差/周期t1以上のようにして図24の処理により、直前直後の位置登録信号それぞれについてウエイトが設定される。ここで、図23を用いて具体例を説明する。図23に示すように、ユーザαに関する位置登録信号として、対象時刻1400に対し、直前の位置登録信号位置登録時刻:1330、在圏セクタ:エリアAと直後の位置登録信号位置登録時刻:1500、在圏セクタ:エリアBとが抽出されたケースを想定する。周期的位置登録の周期tは54分とする。このケースでは、直前の位置登録時刻と対象時刻との第1の時間差mは30分、直後の位置登録時刻と対象時刻との第2の時間差nは60分となり、第2の時間差nが周期t54分より長いので、直後の位置登録信号についてのウエイトは0となる。一方、第1の時間差mは周期t54分以下なので、直前の位置登録信号についてのウエイトは、上記の式1により、54-30/54=24/54となる。そのため、ユーザαは対象時刻1400に、直前の位置登録信号に係る在圏セクタであるエリアAに24/54人分のウエイトで存在するという算出結果が得られる。ここで、図23の例を用いて、上記のようにしてウエイトを算出する考え方を説明する。例えば図23のように、対象時刻の直後の位置登録信号が、時間差が異常値である位置登録信号とみなされ、直前のエリアAでの位置登録信号のみに基づいてウエイトを算出するケースを想定すると、本来は、当該エリアAでの位置登録信号の生成時1330の54分後に周期的位置登録の位置登録信号が生成されるはずである。しかし、54分後に周期的位置登録の位置登録信号が生成されないとすると、54分後の時点1424では移動機の電源がオフもしくは圏外とみなすことができる。このとき、移動機の電源がオフもしくは圏外になった時刻は、上記エリアAでの位置登録信号の生成時1330から54分間の時間帯即ち、1330~1424におけるいずれかの時刻である。ここで、移動機の電源がオフもしくは圏外になるというイベントの発生する確率が、時系列的に一様であるとすると、例えば、上記時間帯1330~1424の途中である対象時刻1400では、対象時刻とエリアAでの位置登録信号の生成時刻との時間差30分/54分の確率で、移動機の電源がオフもしくは圏外になっており、一方、54-30/54=24/54の確率で、移動機がエリアAに存在する。したがって、このような考え方に基づき、上記の例では、対象時刻1400でユーザαはエリアAに24/54人分の重みで存在するものとする。換言すれば、対象時刻1400でユーザαについて30/54人分は移動機の電源がオフもしくは圏外の扱いにする。以上の第5実施形態により、位置登録時刻と対象時刻との時間差が周期的位置登録の周期tよりも長い位置登録信号即ち、時間差が異常値である位置登録信号を除外した上で、より精度の良いウエイト算出処理を行うことができる。さらに、時間差が正常値である位置登録信号については、周期t54分と対象時刻~信号生成時刻の時間差との差分に応じてウエイトを算出することで、時間差が正常値である位置登録信号のウエイトをより精度良く求めることができる。[第6実施形態]以下の第6実施形態は、前述した第2実施形態の変形例に相当し、位置登録信号のうち、LA跨りによる位置登録信号を除去する前の位置登録信号から求めたウエイトと、LA跨りによる位置登録信号を除去した後の位置登録信号から求めたウエイトとを按分し、按分後のウエイトから人口分布を求める処理について説明する。つまり、LA境界に面したセクタでは、LA境界に面していないセクタに比べ、LA跨りによる位置登録信号が生成される確率が相対的に高いため、このようなセクタ間の不均衡に起因した人口分布算出の精度低下を抑えるべく、上記のようにLA跨りによる位置登録信号を除去する前の位置登録信号から求めたウエイトと、LA跨りによる位置登録信号を除去した後の位置登録信号から求めたウエイトとを按分するものである。第6実施形態では、図1の通信システム10のシステム構成、および図2の通信システム10の機能ブロック構成は第1実施形態と同様であるので、重複した説明は省略する。人口分布算出部604により実行される処理の流れは、第2実施形態で説明した図8の処理の流れと同様であるが、図8のステップS202のウエイト算出処理およびステップS203のウエイト集計処理が、第2実施形態とは異なるため、これらの処理について、図25、図26を参照して、以下に説明する。抽出部603は、対象時刻入力部606から対象時刻情報を受け取った後、図25に示すように、各ユーザについて、LA跨りによる位置登録信号を除去する前の位置登録信号から、対象時刻の直前の位置登録信号以下直前の除去前位置登録信号というおよび直後の位置登録信号以下直後の除去前位置登録信号というを抽出するとともに、LA跨りによる位置登録信号を除去した後の位置登録信号から、対象時刻の直前の位置登録信号以下直前の除去後位置登録信号というおよび直後の位置登録信号以下直後の除去後位置登録信号というを抽出する図25のステップS401。なお、位置登録信号には、当該信号が生成された要因LA跨り、移動機100の電源オン等によるアタッチ処理の実行、移動機100の電源オフ等によるデタッチ処理の実行、周期的位置登録などに応じた情報例えば識別フラグが含まれており、抽出部603は、ある位置登録信号がLA跨りによる位置登録信号であるか否かを、例えば当該位置登録信号に含まれる識別フラグによって判別することで、LA跨りによる位置登録信号を除去することができる。抽出された各ユーザについての4つの位置登録信号は人口分布算出部604へ渡され、人口分布算出部604は、各ユーザについて、直前の除去前位置登録信号の位置登録時刻と対象時刻との第1の時間差m、直後の除去前位置登録信号の位置登録時刻と対象時刻との第2の時間差n、直前の除去後位置登録信号の位置登録時刻と対象時刻との第3の時間差q、および、直後の除去後位置登録信号の位置登録時刻と対象時刻との第4の時間差rを算出するステップS402。各ユーザについて、第2実施形態と同様の手順で、上記の第1~第4の時間差から、直前の除去前位置登録信号のウエイトw1、直後の除去前位置登録信号のウエイトw2、直前の除去後位置登録信号のウエイトw3、および、直後の除去後位置登録信号のウエイトw4を算出するステップS403。例えば、以下の式2~5によって各ウエイトw1~w4を算出することができる。直前の除去前位置登録信号のウエイトw1=n/mn2直後の除去前位置登録信号のウエイトw2=m/mn3直前の除去後位置登録信号のウエイトw3=r/qr4直後の除去後位置登録信号のウエイトw4=q/qr5そして、人口分布算出部604は、予め定められた按分パラメータpに基づいて、各ユーザについて、直前および直後の除去前位置登録信号のウエイトと、直前および直後の除去後位置登録信号のウエイトとを按分するステップS404。ここでの按分パラメータpおよび按分については、図26を用いて後述する。さらに、人口分布算出部604は、集計対象全ユーザについての、按分後の直前および直後の除去前位置登録信号のウエイト、並びに、按分後の直前および直後の除去後位置登録信号のウエイトを、セクタIDごとに集計するステップS405。以下、図26を用いて具体例を説明する。図26に例示したように、対象時刻1400の前後の位置登録信号として、エリアAで1320に生成された周期位置登録による位置登録信号、エリアBで1330に生成されたLA跨りによる位置登録信号、および、エリアCで1430に生成された周期位置登録による位置登録信号が存在するケースでは、直前の除去前位置登録信号は、エリアAで1320に生成された位置登録信号となり、直前の除去後位置登録信号は、エリアBで1330に生成された位置登録信号となる。一方、直後の除去前位置登録信号、直後の除去後位置登録信号はともに、エリアCで1430に生成された位置登録信号となる。このとき、上記の第1の時間差m=30分、第2の時間差n=30分、第3の時間差q=40分、第4の時間差r=30分となるため、各ウエイトw1~w4は以下のように算出される。直前の除去前位置登録信号のウエイトw1=30/3030=1/2直後の除去前位置登録信号のウエイトw2=30/3030=1/2直前の除去後位置登録信号のウエイトw3=30/4030=3/7直後の除去後位置登録信号のウエイトw4=40/4030=4/7ここで、LA跨りによる位置登録信号を除去する前の位置登録信号から求めたウエイト直前および直後の除去前位置登録信号のウエイトw1、w2の集合をxとし、LA跨りによる位置登録信号を除去した後の位置登録信号から求めたウエイト直前および直後の除去後位置登録信号のウエイトw3、w4の集合をyとすると、ユーザαの対象時刻1400の位置を求めるために、上記x、yおよび按分パラメータp0<p<1を以下の式6に適用することで、上記xとyを按分する。px1-py6按分パラメータpは、直前および直後の除去前位置登録信号のウエイトと直前および直後の除去後位置登録信号のウエイトとを按分することで人口分布算出を適正化するためのパラメータであり、その値は、例えば、ユーザの移動速度や、LA境界におけるLA跨りによる位置登録と周期的位置登録との発生比率などから、予め定められる。例えば、按分パラメータp=05の場合、図26の最下段に示すように、w1、w2の集合xと、ウエイトw3、w4の集合yとを、05ずつに按分することになる。按分後は、エリアBで生成された直前の除去前位置登録信号のウエイトが14、エリアCで生成された直後の除去前位置登録信号のウエイトが14、エリアAで生成された直前の除去後位置登録信号のウエイトが314、エリアCで生成された直後の除去後位置登録信号のウエイトが27となる。そして、上記の按分後の各ウエイトがセクタIDごとに集計される。上記の例では、ユーザαの対象時刻1400の位置に関するウエイトとして、エリアAについて314、エリアBについて14、エリアCについて1528即ち、14と27の和が求められる。以上の第6実施形態のように、LA跨りによる位置登録信号を除去する前の位置登録信号から求めたウエイトと、LA跨りによる位置登録信号を除去した後の位置登録信号から求めたウエイトとを按分することで、LA跨りによる位置登録信号の生成される確率がセクタ間例えばLA境界に面したセクタと面していないセクタとの間で異なるというセクタ間の不均衡に起因した人口分布算出の精度低下を抑えることができる。なお、前述した第4又は第5実施形態のウエイト算出処理は、上記第6実施形態に適用することができる。即ち、図25のステップS403のウエイト算出処理に対し、第4又は第5実施形態のウエイト算出処理を適用することができる。その場合、位置登録時刻と対象時刻との時間差が周期的位置登録の周期tよりも長い位置登録信号即ち、時間差が異常値である位置登録信号を除外した上で、適正なウエイト算出処理を実行できるとともに、上記のようなセクタ間の不均衡に起因した人口分布算出の精度低下を抑えることができる。エリア変換処理についてまた、上記の第1~第6実施形態により、対象時刻におけるセクタ単位の人口分布が得られるが、集計単位であるセクタ単位の人口分布を、以下のようにして、集計単位とは異なる出力単位ここでは一例としてメッシュごとの推計値に変換してもよい。以下の変換処理は、人口分布算出部604にて実行してもよいし、出力部605にて実行してもよいが、以下では、出力部605にて実行する例を述べる。以下では、単一の周波数帯を利用する屋外局の通信領域が存在する環境における処理と、屋内局の通信領域および電波到達範囲が異なる周波数帯を利用する複数の屋外局の通信領域のうち2つ以上が地理的に重複して存在する環境における処理とを順に説明する。まず、単一の周波数帯を利用する屋外局の通信領域が存在する環境における処理について説明する。出力部605は、予め記憶したエリア境界情報に基づいて再現されるエリア図図27a参照と、所定の区画分け規則に基づき再現される2次元的なメッシュ図27b参照とを合成し、図27cに示すような合成図を得る。次に、出力部605は、上記合成図においてメッシュ境界により各エリアを分割する。例えば、図28に示すように、図27aのエリアAは、メッシュ境界により4つの分割エリアA-1、A-2、A-3、A-4に分割される。そして、出力部605は、各分割エリアの面積を算出し、各分割エリアの面積比を算出する。例えば図28に示すように、分割エリアA-1、A-2、A-3、A-4の面積として、それぞれ10m2、50m2、100m2、40m2が算出されたとすると、分割エリアA-1、A-2、A-3、A-4の面積比例えば百分率として、5、25、50、20が算出される。なお、出力部605が各分割エリアの面積比を算出することは必須ではない。例えば、各分割エリアの面積比は予め求められており、出力部605は、各分割エリアの面積比の情報を、情報分析装置600内の図示しないテーブル又は外部から参照可能とされた構成であってもよい。次に、出力部605は、各分割エリアの人口を算出する。例えば図27aのエリアAの人口が800人であったとすると、図27に示すように分割エリアA-2の人口として、200人即ち、800人×25が算出される。同様に、エリアB、Cの人口がそれぞれ500人、750人であったとすると、エリアBにおける面積比80の分割エリアB-1の人口として、400人即ち、500人×80が算出され、エリアCにおける面積比80の分割エリアC-4の人口として、600人即ち、750人×80が算出される。さらに、出力部605は、1つのメッシュに内包された複数の分割エリアの人口の総和を算出することで、当該メッシュの人口を算出する。図29の例では、1つのメッシュに内包された分割エリアA-2、B-1、C-4の人口の総和1200人即ち、200人400人600人を算出し、この1200人を当該メッシュの人口とする。以上のようにして、単一の周波数帯を利用する屋外局の通信領域が存在する環境において、集計単位ごとの人口を出力単位ごとの人口に変換することができる。次に、屋内局の通信領域および電波到達範囲が異なる周波数帯を利用する複数の屋外局の通信領域のうち2つ以上が地理的に重複して存在する環境における処理について説明する。図30に示すように、屋内局の通信領域および電波到達範囲が異なる周波数帯屋外2GHz17GHzと屋外800MHzを利用する複数の屋外局の通信領域が地理的に重複して存在する環境では、出力部605は、それぞれの通信領域について、上述した単一の周波数帯を利用する屋外局の通信領域が存在する環境における処理を行うことで、それぞれの通信領域に関する出力単位メッシュごとの人口を求め、最後に、出力単位ごとに各通信領域に関する人口を合算することで、出力単位ごとの人口を得る。図30の例では、出力部605は、まず、屋外2GHz17GHzを利用する屋外局の通信領域、屋外800MHzを利用する屋外局の通信領域、屋内局の通信領域それぞれについて、上述した単一の周波数帯を利用する屋外局の通信領域が存在する環境における処理を行う。例えば、屋外2GHz17GHzを利用する屋外局の通信領域において、出力単位QとエリアAとが重なった分割エリアがエリアA全体に対し40の面積比であるとすると、エリアAの推計人口100人に面積比04を乗算することで、出力単位QとエリアAとが重なった分割エリアの推計人口40人が得られる。同様にして、出力単位QとエリアBとが重なった分割エリアについての推計人口3人エリアBの推計人口30人×面積比01、および、出力単位QとエリアCとが重なった分割エリアの推計人口5人エリアCの推計人口100人×面積比005が得られる。屋外800MHzを利用する屋外局の通信領域についても同様に、出力単位QとエリアDとが重なった分割エリアの推計人口3人エリアDの推計人口10人×面積比03、および、出力単位QとエリアFとが重なった分割エリアの推計人口9人エリアFの推計人口30人×面積比03が得られる。一方、屋内局については、個々の屋内局の電波到達範囲であるエリアは非常に小さくて、図30の例では、1つの屋内局のエリアLの全体が出力単位Qと重なっているため、100の面積比と考えることができる。そこで、エリアLの推計人口10人に面積比10を乗算することで、出力単位QとエリアLとが重なったエリアこの例ではエリアL全体の推計人口10人が得られる。最後に、出力部605は、上記のようにして得られた、出力単位Qと各エリアとが重なったエリアの推計人口を合算することで、出力単位Qの推計人口70人を得る。以上のようにして、集計単位ごとの人口から出力単位Qの推計人口へ変換することができる。図30は、1つの出力単位Qの推計人口への変換を示しているが、同様の処理を他の出力単位についても実行することで、対象とする全ての出力単位の推計人口への変換を行うことができる。以上のようにして、屋内局の通信領域および電波到達範囲が異なる周波数帯を利用する複数の屋外局の通信領域のうち2つ以上が地理的に重複して存在する環境において、集計単位ごとの人口を出力単位ごとの人口に変換することができる。なお、上記のエリア変換処理は、セクタ単位の人口を、セクタ単位とは異なる単位例えばメッシュ単位の人口に変換する場合だけでなく、第1実施形態で述べた拡大推計処理端末数を人口に変換する処理の前に、実行してもよい。つまり、セクタ単位の端末数を、上記のエリア変換処理によって、セクタ単位とは異なる単位例えばメッシュ単位の端末数に変換し、その後、変換後の単位例えばメッシュ単位の端末数を、拡大推計処理によって人口に変換することで、変換後の単位例えばメッシュ単位の人口分布を求めても良い。第1~第6実施形態における人口分布算出では、ユーザの属性ごと、時間帯ごと、場所ごとに、人口分布算出の条件や区分を設定することができる。例えば、人口分布算出に用いられる位置データ例えば位置登録信号を、特定の属性を持つユーザの位置データ例えば位置登録信号に絞り込むことにより、ユーザの属性ごとに人口分布算出の条件を設定できる。また、対象時刻を特定の時間帯の範囲内に絞り込むことにより、時間帯の条件を設定できる。その際、対象時刻は所定周期に基づく複数の日例えば平日月曜日~金曜日における特定時刻に設定してもよい。さらに、上記のエリア変換処理によって、セクタ単位の人口又は端末数を、セクタ単位とは異なる単位例えばメッシュ単位の人口又は端末数に変換することで、人口分布算出におけるエリア区分を自在に設定できる。[第7実施形態]図13は、第7実施形態の通信システム10aのシステム構成図である。図13に示すように、この通信システム10aは、新しい通信方式であるLTELongTermEvolutionに適用した場合のシステム構成であり、移動機100、eNBEvolutionNodeB250、交換機400、および管理センタ500を含んで構成されている。また、この管理センタ500は、社会センサユニット501、ペタマイニングユニット502、モバイルデモグラフィーユニット503、および可視化ソリューションユニット504から構成されている。なお、eNB250は、BTS200およびRNC300の両方の機能を包含したものである。この第7実施形態は、LTEに適用した場合のシステム構成であって、その具体的な処理内容は、上述の第1~第6実施形態と同様であるため、その詳細説明は割愛する。なお、第1~第6実施形態において、そのプロトコルは、RANAPRadioAccessNetworkApplicationPartであるが、第7実施形態においては、LTEで使用するS1APS1Applicationprotocolが用いられており、InitialUEMessageについては、S1APにおいても同信号が用いられている。また、第1~第7実施形態においては、第三世代携帯電話3Gシステムを前提として説明していたが、GSMGlobalSystemforMobileCommunicationsにおいても適用可能である。