(57)【要約】【課題】より詳細な移動機の分布を精度良く算出する。【解決手段】移動通信システム１０の補正データ算出部６０４が、分割エリア内に存在する移動機１００の移動機数に基づいて、セクタ内における移動機１００の分布を算出するための補正データを算出する。そして、分布算出部６０８が、算出された補正データとセクタ内の移動機数とに基づいて、セクタ内における移動機１００の分布を算出する。これにより、セクタ内において、より詳細な移動機１００の分布を精度良く算出することができる。【選択図】図２
2012054921.tif 000002
【特許請求の範囲】
【請求項１】
所定エリア内の移動機の移動機数を取得する第１取得手段と、
前記所定エリアを構成する複数の分割エリア内の前記移動機の移動機数を前記分割エリア毎に取得する第２取得手段と、
前記第２取得手段が取得する前記移動機数に基づいて、前記第１取得手段によって取得された前記移動機数から前記所定エリア内における前記移動機の分布を算出するための補正データを算出する補正データ算出手段と、
前記補正データと前記第１取得手段によって取得された前記移動機数とに基づいて、前記所定エリア内における前記移動機の分布を算出する分布算出手段と、
を備えることを特徴とする移動機分布算出システム。
【請求項２】
前記補正データ算出手段は、前記補正データとして、前記第２取得手段によって取得された前記移動機数に基づき、前記分割エリア毎に前記所定エリア内における前記移動機の出現確率を算出し、
前記分布算出手段は、前記出現確率と前記第１取得手段によって取得された前記移動機数とに基づき、前記所定エリア内の前記移動機の移動機数を前記分割エリア毎に振り分けることにより、前記所定エリア内における前記移動機の分布を算出する、
ことを特徴とする請求項１に記載の移動機分布算出システム。
【請求項３】
前記第２取得手段は、前記移動機の存在位置を示す座標情報を取得するものであり、取得した前記座標情報に基づいて前記移動機がどの前記分割エリア内に存在しているかを把握し、把握結果に基づいて前記移動機数を前記分割エリア毎に取得することを特徴とする請求項１又は２に記載の移動機分布算出システム。
【請求項４】
前記第１取得手段は、前記移動機が位置登録を行う際に出力された位置登録信号に基づいて前記移動機数を取得するものであり、
前記第２取得手段は、前記座標情報のうち、前記移動機が前記座標情報を取得した時刻と、前記位置登録が行われた時刻とを比較し、前記位置登録が行われた時刻との差が所定時間以内の前記座標情報に基づいて、前記分割エリア毎に前記移動機の移動機数を取得することを特徴とする請求項３に記載の移動機分布算出システム。
【請求項５】
前記第２取得手段は、一の前記移動機から前記座標情報を複数受信した場合、受信した前記各座標情報に対して前記移動機数の計数についての重み付けを行い、該重み付けの結果に基づいて前記分割エリア毎に前記移動機数を取得することを特徴とする請求項３又は４に記載の移動機分布算出システム。
【請求項６】
前記移動機を識別する識別子と、前記移動機の存在位置に関する情報と、前記存在位置に関する情報が取得された位置取得時刻情報と、を含む位置情報を記憶する第１位置情報記憶手段を更に備え、
前記第１取得手段は、
前記第１位置情報記憶手段に記憶された第１の位置情報について、当該第１の位置情報と同一の識別子を含む位置情報のうち、当該第１の位置情報の直前の位置情報である第２の位置情報の位置取得時刻情報、及び当該第１の位置情報の直後の位置情報である第３の位置情報の位置取得時刻情報を取得する第１前後位置情報取得手段と、
前記第１の位置情報の位置取得時刻情報、前記第２の位置情報の位置取得時刻情報及び前記第３の位置情報の位置取得時刻情報のうち２つ以上に基づいて、前記第１の位置情報についての特徴量を計算する特徴量計算手段と、
集計すべき集計時間帯に関する集計開始時刻以降であり且つ集計終了時刻以前である前記位置取得時刻情報を含むと共に前記移動機の存在位置に関する情報が前記所定エリアに対応づけられる１又は複数の位置情報についての特徴量、及び前記集計開始時刻と前記集計終了時刻との差である集計時間帯の時間長に基づいて、集計時間帯における前記所定エリア内の移動機数を推計する移動機数推計手段と、
を備えることを特徴とする請求項１又は２に記載の移動機分布算出システム。
【請求項７】
前記移動機を識別する識別子と、前記移動機の存在位置に関する情報と、前記存在位置に関する情報が取得された位置取得時刻情報と、を含む位置情報を記憶する第２位置情報記憶手段を更に備え、
前記第２取得手段は、
前記第２位置情報記憶手段に記憶された第１の位置情報について、当該第１の位置情報と同一の識別子を含む位置情報のうち、当該第１の位置情報の直前の位置情報である第２の位置情報の位置取得時刻情報、及び当該第１の位置情報の直後の位置情報である第３の位置情報の位置取得時刻情報を取得する第２前後位置情報取得手段と、
前記第１の位置情報の位置取得時刻情報、前記第２の位置情報の位置取得時刻情報及び前記第３の位置情報の位置取得時刻情報のうち２つ以上に基づいて、前記第１の位置情報についての特徴量を計算する第２特徴量計算手段と、
集計すべき集計時間帯に関する集計開始時刻以降であり且つ集計終了時刻以前である前記位置取得時刻情報を含むと共に前記移動機の存在位置に関する情報が前記分割エリアに対応づけられる１又は複数の位置情報についての特徴量、及び前記集計開始時刻と前記集計終了時刻との差である集計時間帯の時間長に基づいて、集計時間帯における前記分割エリア内の移動機数を推計する第２移動機数推計手段と、
を備えることを特徴とする請求項１，２及び６のいずれか一項に記載の移動機分布算出システム。
【請求項８】
前記移動機を識別する識別子と、前記移動機の存在位置に関する情報と、前記存在位置に関する情報が取得された位置取得時刻情報と、を含む位置情報を記憶する第１位置情報記憶手段を更に備え、
前記第１取得手段は、
前記第１位置情報記憶手段に記憶された位置情報のうち、前記位置取得時刻情報が示す時刻が集計すべき集計時間帯内又は集計時間帯を所定幅拡張した拡張期間内にある前記位置情報に基づいて、集計時間帯の少なくとも一部において前記所定エリアに滞在したと推定される移動機、又は当該移動機にて集計時間帯又は拡張期間内に生成された１つの位置情報を抽出する抽出手段と、
抽出された位置情報数又は移動機の数に基づいて、集計時間帯内における前記所定エリア内の移動機数を推計する移動機数推計手段と、
を備えることを特徴とする請求項１又は２に記載の移動機分布算出システム。
【請求項９】
前記第１取得手段及び前記第２取得手段は、それぞれ所定の集計時間帯における前記移動機数を取得することを特徴とする請求項１〜８のいずれか一項に記載の移動機分布算出システム。
【請求項１０】
前記移動機のユーザの属性情報を記憶する属性情報記憶手段を更に備え、
前記第１取得手段及び前記第２取得手段は、前記属性情報記憶手段に記憶された前記属性情報の属性毎の前記移動機数を取得することを特徴とする請求項１〜９のいずれか一項に記載の移動機分布算出システム。
【請求項１１】
前記分布算出手段は、移動機数を人口に変換するための拡大係数を用いて、移動機数の分布を人口の分布に変換することを特徴とする請求項１〜１０のいずれか一項に記載の移動機分布算出システム。
【請求項１２】
所定エリア内の移動機の移動機数を取得する第１取得ステップと、
前記所定エリアを構成する複数の分割エリア内の前記移動機の移動機数を前記分割エリア毎に取得する第２取得ステップと、
前記第２取得ステップで取得された前記移動機数に基づいて、前記第１取得ステップで取得された前記移動機数から前記所定エリア内における前記移動機の分布を算出するための補正データを算出する補正データ算出ステップと、
前記補正データと前記第１取得ステップで取得された前記移動機数とに基づいて、前記所定エリア内における前記移動機の分布を算出する分布算出ステップと、
を備えることを特徴とする移動機分布算出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、移動機の分布を算出する移動機分布算出システム及び移動機分布算出方法に関する。
【背景技術】
【０００２】
従来、移動機から発信された位置登録信号に基づいて人口分布を把握する装置として、例えば特許文献１に記載されたものがある。この装置では、移動機から発信される位置登録信号を各基地局が管理するエリア毎に収集し、該エリア毎に人口を算出することにより人口分布を把握している。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００２−３４２５５７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ところで、広範囲の人口分布メッシュ（例えば１ｋｍ×１ｋｍ程度のメッシュ）を、より詳細な（狭い範囲の）人口分布メッシュ（例えば１００ｍ×１００ｍ程度のメッシュ）に分割する技術がある。この技術は、国勢調査等で得られた世帯データや電話帳データ等を利用し、広範囲の人口分布メッシュの中で、どの部分に人口が集中しているか等を把握し、集中度合いに応じて、広範囲の人口分布メッシュ内の人口をより狭い範囲の人口分布メッシュ毎に振り分けるものである。このような技術を用いることにより、特許文献１に記載されているような移動機から発信された位置登録信号に基づいて把握した人口分布から、より詳細な人口分布を得ることが考えられる。
【０００５】
しかしながら、例えば、昼間等、人が活動をしている時間帯においては、就労のために居住地から移動している等により、世帯データや電話帳データから得られる人口の集中度合いと、実際の人口の集中度合いとの差が大きくなる傾向があるので、より詳細な人口分布を精度良く得ることは困難であるという問題がある。
【０００６】
そこで本発明は、より詳細な移動機の分布を精度良く算出することができる移動機分布算出システム及び移動機分布算出方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明は、所定エリア内の移動機の移動機数を取得する第１取得手段と、所定エリアを構成する複数の分割エリア内の移動機の移動機数を分割エリア毎に取得する第２取得手段と、第２取得手段が取得する移動機数に基づいて、第１取得手段によって取得された移動機数から所定エリア内における移動機の分布を算出するための補正データを算出する補正データ算出手段と、補正データと第１取得手段によって取得された移動機数とに基づいて、所定エリア内における移動機の分布を算出する分布算出手段と、を備えることを特徴とする。
【０００８】
また、本発明は、所定エリア内の移動機の移動機数を取得する第１取得ステップと、所定エリアを構成する複数の分割エリア内の移動機の移動機数を分割エリア毎に取得する第２取得ステップと、第２取得ステップで取得された移動機数に基づいて、第１取得ステップで取得された移動機数から所定エリア内における移動機の分布を算出するための補正データを算出する補正データ算出ステップと、補正データと第１取得ステップで取得された移動機数とに基づいて、所定エリア内における移動機の分布を算出する分布算出ステップと、を備えることを特徴とする。
【０００９】
これらの発明にあっては、分割エリア内に存在する移動機の移動機数に基づいて、所定エリア内における移動機の分布を算出するための補正データが算出される。そして、算出された補正データと、所定エリア内の移動機の移動機数とに基づいて、所定エリア内における移動機の分布が算出される。これにより、所定エリア内において、より詳細な移動機の分布を精度良く算出することができる。例えば、第２取得手段（第２取得ステップ）によって取得される移動機数が実際に分割エリア内に存在する移動機の移動機数より少ない場合であっても、取得された移動機数に基づいて、所定エリア内における移動機の分布を算出するための補正データを算出することができる。従って、算出された補正データを用いることにより、所定エリア内におけるより詳細な移動機の分布を精度良く算出することができる。
【００１０】
また、補正データ算出手段は、補正データとして、第２取得手段によって取得された移動機数に基づき、分割エリア毎に所定エリア内における移動機の出現確率を算出し、分布算出手段は、出現確率と第１取得手段によって取得された移動機数とに基づき、所定エリア内の移動機の移動機数を分割エリア毎に振り分けることにより、所定エリア内における移動機の分布を算出する、ことが好適である。これにより、出現確率と第１取得手段によって取得された移動機数とに基づいて、所定エリア内における移動機の分布を分割エリア毎に精度良く算出することができる。
【００１１】
また、第２取得手段は、移動機の存在位置を示す座標情報を取得するものであり、取得した座標情報に基づいて移動機がどの分割エリア内に存在しているかを把握し、把握結果に基づいて移動機数を分割エリア毎に取得することが好適である。このように第２取得手段は、移動機の存在位置を示す座標情報に基づいて、分割エリア毎に移動機数を容易に取得することができる。
【００１２】
また、第１取得手段は、移動機が位置登録を行う際に出力された位置登録信号に基づいて移動機数を取得するものであり、第２取得手段は、座標情報のうち、移動機が座標情報を取得した時刻と、位置登録が行われた時刻とを比較し、位置登録が行われた時刻との差が所定時間以内の座標情報に基づいて、分割エリア毎に移動機の移動機数を取得することが好適である。例えば、移動機が位置登録を行ったタイミングと、移動機が座標情報を取得したタイミングとが時間的に近い場合、位置登録を行った時と、移動機が座標情報を取得した時とにおいて移動機はほぼ同じ位置にいる（ほとんど移動していない）と考えられる。このため、両タイミングの時間差が所定時間以内の座標情報を用いることにより、位置登録が行われた時に実際に移動機が存在していた分割エリアを把握し、把握結果に基づいて分割エリア毎に移動機数を取得することができる。これにより、座標情報に基づいて算出される補正データが、位置登録信号に基づいて取得された移動機数の情報に対して時間的により正確に対応するものとなり、所定エリア内の移動機の分布をより精度良く算出することができる。
【００１３】
また、第２取得手段は、一の移動機から座標情報を複数受信した場合、受信した各座標情報に対して移動機数の計数についての重み付けを行い、該重み付けの結果に基づいて分割エリア毎に移動機数を取得することが好適である。これにより、重み付けの結果に基づいて分割エリア毎に移動機の移動機数をより現実に則して取得することができる。
【００１４】
移動機を識別する識別子と、移動機の存在位置に関する情報と、存在位置に関する情報が取得された位置取得時刻情報と、を含む位置情報を記憶する第１位置情報記憶手段を更に備え、第１取得手段は、第１位置情報記憶手段に記憶された第１の位置情報について、当該第１の位置情報と同一の識別子を含む位置情報のうち、当該第１の位置情報の直前の位置情報である第２の位置情報の位置取得時刻情報、及び当該第１の位置情報の直後の位置情報である第３の位置情報の位置取得時刻情報を取得する第１前後位置情報取得手段と、第１の位置情報の位置取得時刻情報、第２の位置情報の位置取得時刻情報及び第３の位置情報の位置取得時刻情報のうち２つ以上に基づいて、第１の位置情報についての特徴量を計算する特徴量計算手段と、集計すべき集計時間帯に関する集計開始時刻以降であり且つ集計終了時刻以前である位置取得時刻情報を含むと共に移動機の存在位置に関する情報が所定エリアに対応づけられる１又は複数の位置情報についての特徴量、及び集計開始時刻と集計終了時刻との差である集計時間帯の時間長に基づいて、集計時間帯における所定エリア内の移動機数を推計する移動機数推計手段と、を備えることが好ましい。この場合には、特徴量に基づいて、所定エリア内の移動機数をより正確に推計することができる。
【００１５】
移動機を識別する識別子と、移動機の存在位置に関する情報と、存在位置に関する情報が取得された位置取得時刻情報と、を含む位置情報を記憶する第２位置情報記憶手段を更に備え、第２取得手段は、第２位置情報記憶手段に記憶された第１の位置情報について、当該第１の位置情報と同一の識別子を含む位置情報のうち、当該第１の位置情報の直前の位置情報である第２の位置情報の位置取得時刻情報、及び当該第１の位置情報の直後の位置情報である第３の位置情報の位置取得時刻情報を取得する第２前後位置情報取得手段と、第１の位置情報の位置取得時刻情報、第２の位置情報の位置取得時刻情報及び第３の位置情報の位置取得時刻情報のうち２つ以上に基づいて、第１の位置情報についての特徴量を計算する第２特徴量計算手段と、集計すべき集計時間帯に関する集計開始時刻以降であり且つ集計終了時刻以前である位置取得時刻情報を含むと共に移動機の存在位置に関する情報が分割エリアに対応づけられる１又は複数の位置情報についての特徴量、及び集計開始時刻と集計終了時刻との差である集計時間帯の時間長に基づいて、集計時間帯における分割エリア内の移動機数を推計する第２移動機数推計手段と、を備えることが好ましい。この場合には、特徴量に基づいて、分割エリア内の移動機数をより正確に推計することができ、この移動機数に基づいて補正データの算出が可能となる。
【００１６】
移動機を識別する識別子と、移動機の存在位置に関する情報と、存在位置に関する情報が取得された位置取得時刻情報と、を含む位置情報を記憶する第１位置情報記憶手段を更に備え、第１取得手段は、第１位置情報記憶手段に記憶された位置情報のうち、位置取得時刻情報が示す時刻が集計すべき集計時間帯内又は集計時間帯を所定幅拡張した拡張期間内にある位置情報に基づいて、集計時間帯の少なくとも一部において所定エリアに滞在したと推定される移動機、又は当該移動機にて集計時間帯又は拡張期間内に生成された１つの位置情報を抽出する抽出手段と、抽出された位置情報数又は移動機の数に基づいて、集計時間帯内における所定エリア内の移動機数を推計する移動機数推計手段と、を備えることが好ましい。この場合には、２重カウントを回避しつつ移動機数を精度よく推計することができる。
【００１７】
また、第１取得手段及び第２取得手段は、それぞれ所定の集計時間帯における移動機数を取得することが好適である。これにより、所定の集計時間帯における所定エリア内の移動機の移動機数と、所定の集計時間帯における分割エリア毎の移動機の移動機数に基づいて算出された補正データとを得ることができる。従って、この移動機数及び補正データが時間的に整合しているため、これらのデータに基づいて所定の集計時間帯における所定エリア内のより詳細な移動機の分布をより精度良く算出することができる。
【００１８】
また、移動機のユーザの属性情報を記憶する属性情報記憶手段を更に備え、第１取得手段及び第２取得手段は、属性情報記憶手段に記憶された属性情報の属性毎の移動機数を取得することが好適である。これにより、所望の属性のユーザについての移動機の分布を算出することができ、データを有効利用することが可能となる。
【００１９】
分布算出手段は、移動機数を人口に変換するための拡大係数を用いて、移動機数の分布を人口の分布に変換することが好ましい。この場合には、移動機数から人口を求めることができる。
【発明の効果】
【００２０】
本発明によれば、より詳細な移動機の分布を精度良く算出することができる。
【図面の簡単な説明】
【００２１】
【図１】第１実施形態における移動通信システムの概略構成図である。
【図２】第１実施形態及び変形例における移動通信システムの機能構成を示すブロック図である。
【図３】ＢＴＳとセクタとの関係を示す図である。
【図４】ＧＰＳ位置情報を示すデータテーブルである。
【図５】分割エリア対応表記憶部に記憶された分割エリア対応表を示すデータテーブルである。
【図６】ＧＰＳ位置情報と分割エリア対応表との照合結果を示すデータテーブルである。
【図７】補正データを示すデータテーブルである。
【図８】位置登録情報を示すデータテーブルである。
【図９】セクタ毎の移動機数を示すデータテーブルである。
【図１０】セクタ内における移動機の詳細な分布を算出する手順を示す概念図である。
【図１１】分割エリア毎の移動機数を示すデータテーブルである。
【図１２】集計エリア対応表を示すデータテーブルである。
【図１３】分割エリア毎の移動機数から集計エリア毎の移動機数を集計する手順を示す概念図である。
【図１４】集計エリア毎に対応付けられた移動機数を示すデータテーブルである。
【図１５】集計エリア毎の移動機数を示すデータテーブルである。
【図１６】第１実施形態における補正データを算出する処理の流れを示すフローチャートである。
【図１７】第１実施形態における集計エリア内の移動機数を算出する処理の流れを示すフローチャートである。
【図１８】算出された分割エリア毎の移動機数を示すデータテーブルである。
【図１９】算出された分割エリア毎の移動機数から集計エリア毎の移動機数を集計する手順を示す概念図である。
【図２０】セクタ範囲と集計エリアとの位置関係を示す図である。
【図２１】第２実施形態における移動通信システムの機能構成を示すブロック図である。
【図２２】ＧＰＳ位置情報を示すデータテーブルである。
【図２３】第２実施形態における補正データを算出する処理の流れを示すフローチャートである。
【図２４】第３実施形態における情報分析装置の機能構成を示すブロック図である。
【図２５】第３実施形態に係る第２取得部の機能構成を示すブロック図である。
【図２６】第３実施形態に係る第２取得部の特徴量計算部が作成するデータベースを示す図である。
【図２７】第３実施形態に係る第２取得部の分割エリア照合部が作成するデータベースを示す図である。
【図２８】移動機数推計の考え方を説明するための図である。
【図２９】移動機数推計に係る計算方法を説明するための図である。
【図３０】第３実施形態に係る第１取得部の機能構成を示すブロック図である。
【図３１】第３実施形態に係る第１取得部の特徴量計算部が作成するデータベースを示す図である。
【図３２】第４実施形態に係る情報分析装置の機能構成を示すブロック図である。
【図３３】第４実施形態に係る第１取得部の機能構成を示すブロック図である。
【図３４】推定滞在時間の算出方法を説明するための図である。
【図３５】位置情報の抽出方法を説明するための図である。
【発明を実施するための形態】
【００２２】
以下、図面を参照しつつ本発明の好適な実施形態について詳細に説明する。なお、図面の説明において同一の要素には同一の符号を付し、重複する説明を省略する。
【００２３】
［第１実施形態］
まず、本発明に係る移動機分布算出システム及び移動機分布算出方法を適用した移動通信システムの好適な第１実施形態について説明する。図１は、第１実施形態に係る移動通信システム１０の概略構成図である。図１に示すように、移動通信システム１０は、移動機１００、ＢＴＳ（基地局）２００、ＲＮＣ（無線ネットワーク制御装置）３００、交換機４００、及び管理センタ５００を含んで構成されている。また、この管理センタ５００は、社会センサユニット５０１、ペタマイニングユニット５０２、モバイルデモグラフィユニット５０３、及び可視化ソリューションユニット５０４から構成されている。
【００２４】
交換機４００は、ＢＴＳ２００、ＲＮＣ３００を介して、移動機１００が位置登録要求や通信接続等を行った場合に移動機１００の位置情報（詳しくは後述する位置登録情報及びＧＰＳ位置情報）を取得する。交換機４００は受け取った位置情報を記憶しておき、所定のタイミング、または管理センタ５００からの要求に応じて、収集した位置情報を管理センタ５００に出力する。ここで、一般的に、ＲＮＣ３００は、約千個からなるものであり、日本全国に配置されている。一方で、交換機４００は、３００個程度日本国内に配置されている。
【００２５】
管理センタ５００は、上述したとおり、社会センサユニット５０１、ペタマイニングユニット５０２、モバイルデモグラフィユニット５０３、及び可視化ソリューションユニット５０４を含んで構成されており、各ユニットでは、移動機１００の分布情報を用いた統計処理を行う。
【００２６】
社会センサユニット５０１は、各交換機４００から移動機１００の位置情報等を含んだデータを収集するサーバ装置である。この社会センサユニット５０１は、交換機４００から定期的に出力されたデータを受信したり、または社会センサユニット５０１において予め定められたタイミングに従って交換機４００からデータを取得したりできるように構成されている。
【００２７】
ペタマイニングユニット５０２は、社会センサユニット５０１から受信したデータを所定のデータ形式に変換するサーバ装置である。例えば、ペタマイニングユニット５０２は、ユーザＩＤをキーにソーティング処理を行ったり、データの属するエリア毎にソーティング処理を行ったりする。
【００２８】
モバイルデモグラフィユニット５０３は、ペタマイニングユニット５０２において処理されたデータに対する集計処理、即ち各項目のカウンティング処理を行うサーバ装置である。例えば、モバイルデモグラフィユニット５０３は、あるエリアに在圏する移動機１００のユーザ数をカウントしたり、また在圏分布を集計したりすることができる。
【００２９】
可視化ソリューションユニット５０４は、モバイルデモグラフィユニット５０３において集計処理されたデータを可視可能に処理するサーバ装置である。例えば、可視化ソリューションユニット５０４は、集計されたデータを地図上にマッピング処理することができる。この可視化ソリューションユニット５０４にて処理されたデータは、企業、官公庁または個人等に提供され、店舗開発、道路交通調査、災害対策、環境対策等に利用される。なお、このように統計処理された情報は、当然にプライバシーを侵害しないように個人等は特定されないように加工されている。
【００３０】
なお、社会センサユニット５０１、ペタマイニングユニット５０２、モバイルデモグラフィユニット５０３及び可視化ソリューションユニット５０４はいずれも、前述したようにサーバ装置により構成され、図示は省略するが、通常の情報処理装置の基本構成、即ち、ＣＰＵ、ＲＡＭ、ＲＯＭ、キーボードやマウス等の入力デバイス、外部との通信を行う通信デバイス、情報を記憶する記憶デバイス、及び、ディスプレイやプリンタ等の出力デバイス等を備えることは言うまでもない。
【００３１】
次に、第１実施形態に係る移動通信システム１０の機能ブロック構成について説明する。図２は、移動通信システム１０の機能構成を示すブロック図である。図２に示すように、移動通信システム１０は、移動機１００と、ＢＴＳ２００と、ＲＮＣ３００と、交換機４００と、情報分析装置６００とを含んで構成されている。なお、本実施形態の情報分析装置６００は、図１に示される管理センタ５００の各構成ユニット（社会センサユニット５０１、ペタマイニングユニット５０２、モバイルデモグラフィユニット５０３、及び可視化ソリューションユニット５０４）の全部もしくは一部により実現される。即ち、図２に示される情報分析装置６００の機能ブロックは、管理センタ５００を構成する各構成ユニットにおいて分散して配置されており、各構成ユニットにより一つの情報分析装置６００が構成されている。なお、一の構成ユニットにより情報分析装置６００が構成されてもよい。
【００３２】
移動機１００は、ＧＰＳ機能部１０１と、位置登録部１０２とを含んで構成されている。ＧＰＳ機能部１０１は、ＧＰＳを用いて移動機１００の存在位置を示す詳細な位置情報（座標情報）（以下「ＧＰＳ位置情報」と称する。）を取得するものである。ＧＰＳ機能部１０１によって取得されたＧＰＳ位置情報は、移動機１００に搭載された所定のアプリケーション（ＧＰＳ位置情報を用いるアプリケーション）の実行時に利用される。また、ＧＰＳ機能部１０１は、ＧＰＳ位置情報を取得すると、取得したＧＰＳ位置情報を交換機４００に対して出力する。位置登録部１０２は、位置登録要求を行うための位置登録信号を出力するものである。なお、位置登録信号は、例えば５４分毎に出力される。
【００３３】
ＢＴＳ２００は、通信エリアであるセクタ（所定エリア）を形成する。このセクタ内に、複数の移動機１００が在圏している。ここで、図３に、ＢＴＳ２００とセクタとの関係を示す。図３に示すように、円形で示されている領域の中心にＢＴＳ２００が位置するものであり、それを中心に複数の等分されたものがセクタである。例えば図３では、ＢＴＳ２００の通信エリアは最大６セクタからなるものであり、ＲＮＣ３００は、いずれのセクタに移動機１００が在圏しているかを把握することができる。
【００３４】
ＲＮＣ３００は、通信制御部３０１と、セクタ識別部３０２とを含んで構成されている。通信制御部３０１は、移動機１００とＢＴＳ２００を介して通信接続を行う部分であり、例えば、移動機１００からの発信処理若しくは着信処理に基づいた通信接続処理、及び位置登録要求に基づいた通信接続処理を行う部分である。セクタ識別部３０２は、移動機１００が位置登録のために位置登録信号を出力した際に当該移動機１００が在圏するセクタを識別し、識別結果を通信制御部３０１を通じて交換機４００へ出力する。また、セクタ識別部３０２は、移動機１００がＧＰＳ位置情報を交換機４００に出力したときに、当該移動機１００が在圏するセクタを識別し、識別結果を通信制御部３０１を通じて交換機４００へ出力する。
【００３５】
交換機４００は、通信制御部４０１と、位置情報管理部４０２とを含んで構成されている。通信制御部４０１は、ＲＮＣ３００の通信制御部３０１と同様に、通信接続処理を行う部分である。位置情報管理部４０２は、移動機１００が位置登録を行ったときに、セクタ識別部３０２によって識別されたセクタ（セクタ識別子）と、移動機１００を識別するためのユーザ識別子と、移動機１００が位置登録を行った時刻とを対応付けし、対応付けした結果を位置登録情報として記憶する。また、位置情報管理部４０２は、移動機１００がＧＰＳ位置情報を出力したときに、セクタ識別部３０２によって識別されたセクタ（セクタ識別子）と、移動機１００を識別するためのユーザ識別子と、移動機１００によってＧＰＳ位置情報が取得された時刻と、緯度・経度とを対応付けし、対応付けした結果をＧＰＳ位置情報として記憶する。位置情報管理部４０２は、記憶している位置情報（位置登録情報及びＧＰＳ位置情報）を、所定のタイミング、または情報分析装置６００からの要求に応じて情報分析装置６００に出力する。
【００３６】
情報分析装置６００は、ＧＰＳ位置情報記憶部６０１と、第２取得部（第２取得手段）６０２と、分割エリア対応表記憶部６０３と、補正データ算出部（補正データ算出手段）６０４と、補正データ記憶部６０５と、位置登録情報記憶部６０６と、第１取得部（第１取得手段）６０７と、分布算出部（分布算出手段）６０８と、分布情報記憶部６０９と、集計部６１０と、集計エリア対応表記憶部６１１と、集計情報記憶部６１２とを含んで構成されている。
【００３７】
ＧＰＳ位置情報記憶部６０１は、交換機４００の位置情報管理部４０２から出力されたＧＰＳ位置情報を記憶するものである。具体的には、ＧＰＳ位置情報記憶部６０１は、図４に示すように、ユーザ識別子と、移動機１００によってＧＰＳ位置情報が取得された時刻と、セクタ識別子と、緯度・経度とが対応付けられたＧＰＳ位置情報とを、ＧＰＳ位置情報が取得された時刻毎に記憶している。
【００３８】
第２取得部６０２は、セクタを構成する複数の分割エリア内の移動機数を分割エリア毎に取得するものである。なお、セクタは、集計処理上便宜的に区分けした複数の分割エリアによって構成されている。具体的には、第２取得部６０２は、ＧＰＳ位置情報記憶部６０１に記憶されたＧＰＳ位置情報（図４参照）のうち、所定の抽出条件に対応するＧＰＳ位置情報を抽出する。ここで、抽出条件として、詳しくは後述する集計部６１０において移動機１００の詳細な分布を得ようとする時間帯（以下「集計時間帯」と称する）に対応する時間帯を設定することが好ましい。本実施形態では、集計時間帯を、平日の００：００〜００：５９とし、抽出条件として、過去数ヶ月間の集計時間帯とする。なお、過去数ヶ月間の集計時間帯におけるＧＰＳ位置情報を取得するのは、ＧＰＳ位置情報のデータ数が少ない場合であっても、より精度の良い出現確率（詳しくは後述する）を算出するために十分なデータ数を確保するためである。
【００３９】
抽出の後、第２取得部６０２は、抽出したＧＰＳ位置情報に含まれる緯度・経度と、分割エリア対応表記憶部６０３に記憶された分割エリア対応表（図５参照。詳しくは後述する。）とを用い、抽出したＧＰＳ位置情報のそれぞれがどの分割エリアに対応するかを照合する。第２取得部６０２は、照合の結果として、図６に示すように、ユーザ識別子と、時刻と、セクタ識別子と、分割エリア識別子とを対応付ける。
【００４０】
そして、第２取得部６０２は、対応付けした結果（図６参照）に基づき、セクタ識別子が同一のレコード数Ａと、セクタ識別子が同一かつ分割エリア識別子が同一のレコード数Ｂとを計数する。このレコード数Ａが、セクタ内の移動機１００の移動機数となり、レコード数Ｂが、分割エリア内の移動機１００の移動機数となる。なお、レコード数Ａ，Ｂは、同一の移動機１００が複数のＧＰＳ位置情報を出力した場合の重複を含む。第２取得部６０２は、抽出条件（過去数ヶ月間の平日の００：００〜００：５９）と、レコード数Ａと、当該レコード数Ａに対応するセクタ識別子と、を対応付け、抽出条件と、レコード数Ｂと、当該レコード数Ｂに対応するセクタ識別子及び分割エリア識別子と、を対応付けて、それぞれ補正データ算出部６０４へ出力する。
【００４１】
分割エリア対応表記憶部６０３は、分割エリア対応表を記憶する。分割エリア対応表は、図５に示すように、分割エリア識別子と、分割エリアの位置・形状データを表すポリゴン表現とが対応付けられたものである。ここで、ポリゴン表現とは、緯度・経度の情報と、エリアとを結合するための情報である。
【００４２】
補正データ算出部６０４は、第２取得部６０２から出力された情報に基づいて、第１取得部６０７によって取得されたセクタ内の移動機数（詳しくは後述する）からセクタ内における移動機１００の分布を算出するための補正データを算出する。具体的には、補正データ算出部６０４は、取得したレコード数Ｂを、当該レコード数Ｂに対応付けられた分割エリアを含むセクタのレコード数Ａで除算（例えば、分割エリア識別子Ｙ１のレコード数Ｂ÷分割エリア識別子Ｙ１を含むセクタＸ１のレコード数Ａ）することにより、空間的な出現確率（補正データ）を算出する。この出現確率の算出を分割エリア毎に行うことにより、所定の抽出条件の場合における、セクタに対する移動機ユーザの平均的出現確率分布を算出することができる。
【００４３】
補正データ記憶部６０５は、補正データ算出部６０４が算出した出現確率を、図７に示すように、分割エリア毎に、セクタ識別子と、種別（集計時間帯の平日又は休日の種別）と、集計時間帯と、分割エリア識別子と、算出した出現確率とを対応付けて記憶する。
【００４４】
位置登録情報記憶部６０６は、交換機４００の位置情報管理部４０２から出力された位置登録情報を記憶するものである。具体的には、位置登録情報記憶部６０６は、図８に示すように、ユーザ識別子と、位置登録がされた時刻と、セクタ識別子とが対応付けられた情報を、位置登録がされた時刻毎に記憶している。
【００４５】
第１取得部６０７は、セクタ内の移動機１００の移動機数を取得するものである。具体的には、第１取得部６０７は、集計時間帯（平日の００：００〜００：５９）に位置登録を行った移動機１００の位置登録情報を、位置登録情報（図８参照）の時刻の情報に基づいて位置登録情報記憶部６０６から抽出する。
【００４６】
そして、第１取得部６０７は、抽出した位置登録情報に基づいて、セクタ識別子が同一のレコード数を計数することにより、セクタ毎に、集計時間帯におけるセクタ内の移動機１００の移動機数を計数する。第１取得部６０７は、図９に示すように、種別（集計時間帯の平日又は休日の種別）と、集計時間帯と、セクタ識別子と、セクタ毎の移動機数とを対応付けし、対応付けした結果を分布算出部６０８へ出力する。
【００４７】
分布算出部６０８は、補正データ記憶部６０５に記憶された補正データと、第１取得部６０７から出力されたセクタ毎の移動機数（図９参照）とに基づいて、セクタ内における移動機１００の分布を算出するものである。
【００４８】
ここで、分布算出部６０８が行う処理の概念について説明する。図１０は、セクタ内における移動機１００の詳細な分布を算出する手順を示す概念図である。分布算出部６０８は、図１０（ａ）に示すように、第１取得部６０７からセクタＸについての移動機数（ここでは、１００台とする）を取得する。分布算出部６０８は、補正データ記憶部６０５から、所望の集計時間帯に対応し、かつセクタＸ内の分割エリアＹ毎の出現確率を取得する（図１０（ｂ）参照）。なお、図１０（ｂ）は、最も色が濃い分割エリアの出現確率を「０．１」、次に色が濃い分割エリアの出現確率を「０．０５」、次に色が濃い分割エリアの出現確率を「０．０１」、白色の分割エリアの出現確率を「０」とした場合を示している。
【００４９】
分布算出部６０８は、セクタＸの移動機数（１００台）と、各分割エリアＹの出現確率（「０．１」，「０．０５」，「０．０１」）とを乗算し、図１０（ｃ）に示すように、各分割エリア毎の移動機数（１０台，５台，1台）を算出する。即ち、セクタ内の移動機数（１００台）を、出現確率を用いて分割エリア毎に振り分けることにより、セクタ内における移動機の分布（ユーザの分布）を算出する。
【００５０】
詳細には、分布算出部６０８は、第１取得部６０７から取得したセクタ毎の移動機数（図９参照）と、補正データ記憶部６０５に記憶された補正データ（図７参照）のうち、対応するセクタ内の分割エリアの出現確率とを乗算し、分割エリア毎の移動機数を算出する。
【００５１】
分布情報記憶部６０９は、分布算出部６０８が算出した分割エリア毎の移動機数を、図１１に示すように、分割エリア毎に、セクタ識別子と、種別と、集計時間帯と、分割エリア識別子と、出現確率と、算出した分割エリア毎の移動機数とを対応付けし、対応付けした結果を分布情報として記憶する。
【００５２】
集計エリア対応表記憶部６１１は、図１２に示すように、分割エリア識別子と、集計エリア識別子とが対応付けられた集計エリア対応表を記憶している。
【００５３】
集計部６１０は、分布情報記憶部６０９に記憶された分布情報から、所定の集計エリア内の移動機数を集計するものである。ここで、集計部６１０が行う処理の概念について説明する。図１３は、分割エリア毎の移動機数から集計エリア毎の移動機数を集計する手順を示す概念図である。図１３（ａ）に示すように、集計部６１０は、分布情報記憶部６０９から、所望の集計時間帯かつ所定の集計エリアを含むセクタＸ内の移動機数を分割エリアＹ毎に取得する。なお、図１３（ａ）は、最も色が濃い分割エリアの移動機数を１０台、次に色が濃い分割エリアの移動機数を５台、次に色が濃い分割エリアの移動機数を１台、白色の分割エリアの移動機数を０台とした場合を示している。
【００５４】
図１３（ｂ）に示すように、集計部６１０は、集計エリア対応表記憶部６１１に記憶された集計エリア対応表に基づいて、各分割エリアＹが集計エリアａ〜ｄのうちどの集計エリアに属するかを照合する。そして、集計部６１０は、図１３（ｃ）に示すように、各分割エリアＹが属する集計エリア毎に分割エリア毎に算出された移動機数を合算し、集計エリア内の移動機数を集計する。
【００５５】
より詳細には、集計部６１０は、分布情報記憶部６０９から分布情報（図１１参照）と、集計エリア対応表記憶部６１１に記憶された集計エリア対応表（図１２参照）とを取得する。集計部６１０は、分布情報（図１１参照）の分割エリア識別子を、集計エリア対応表に基づいて対応する集計エリア識別子に置き換えることにより、図１４に示すように、セクタ識別子と、種別と、集計時間帯と、集計エリア識別子と、出現確率と、分割エリア毎の移動機数とが対応付けられた情報を得る。
【００５６】
そして、集計部６１０は、対応付けした結果（図１４参照）のうち集計エリア識別子と分割エリア毎の移動機数とに基づいて、集計エリア識別子毎に分割エリア毎の移動機数を合算することによって集計エリア識別子毎の移動機数を集計情報として算出する。
【００５７】
集計情報記憶部６１２は、集計部６１０が算出した集計情報を、図１５に示すように、集計エリア識別子と、集計時間帯と、集計した集計エリア毎の移動機数とを対応付し、対応付けした結果を記憶する。なお、集計部６１０が算出した集計情報を集計情報記憶部６１２に記憶する以外にも、集計情報を外部のサーバ装置等に出力しても良い。これにより、例えば、集計エリアとして、ある市区町村に対応するエリアが設定されている場合、集計エリア毎の移動機数を市区町村毎の移動機数として求めることができ、また、集計エリア内における分割エリア毎の移動機数に基づいて市区町村内における移動機の分布を求めることができる。
【００５８】
次に、情報分析装置６００で行われる補正データ算出処理及び集計エリア内の移動機数を算出する処理の流れについて説明する。まず、補正データ算出処理について説明する。図１６は、補正データを算出するための処理の流れを示すフローチャートである。なお、本処理は、例えば１週間に１回等、定期的に行われるものである。
【００５９】
まず、第２取得部６０２は、ＧＰＳ位置情報記憶部６０１に記憶されたＧＰＳ位置情報（図４参照）を、抽出条件に従って抽出する（ステップＳ１０１）。そして、第２取得部６０２は、分割エリア対応表記憶部６０３に記憶された分割エリア対応表（図５参照）に基づいて、抽出したＧＰＳ位置情報のそれぞれがどの分割エリアに対応するかを照合する（ステップＳ１０２）。その後、第２取得部６０２は、セクタ毎の移動機数（レコード数Ａ）と、分割エリア毎の移動機数（レコード数Ｂ）とを計数する（ステップＳ１０３：第２取得ステップ）。
【００６０】
次に、補正データ算出部６０４は、補正データとして、レコード数Ｂをレコード数Ａで除算することにより分割エリア毎の出現確率を算出する（ステップＳ１０４：補正データ算出ステップ）。そして、補正データ記憶部６０５は、補正データ算出部６０４が算出した補正データ（図７参照）を記憶する（ステップＳ１０５）。
【００６１】
次に、集計エリア内の移動機数を算出する処理の流れについて説明する。図１７は、集計エリア内の移動機数を算出するための処理の流れを示すフローチャートである。なお、本処理は、１時間に１回等、補正データ算出処理の処理周期とは異なる周期で行われるものである。まず、第１取得部６０７が、位置登録情報記憶部６０６に記憶された位置登録情報（図８参照）から、抽出条件に対応する位置登録情報を抽出する（ステップＳ２０１）。そして、第１取得部６０７は、抽出した位置登録情報に基づいて、セクタ毎に、集計時間帯におけるセクタ内の移動機１００の移動機数を計数（図９参照）する（ステップＳ２０２：第１取得ステップ）。
【００６２】
次に、分布算出部６０８は、補正データ記憶部６０５から、所望の集計時間帯に対応する補正データ（図７参照）を取得する（ステップＳ２０３）。そして、分布算出部６０８は、第１取得部６０７で計数された各セクタ毎の移動機数（図９参照）と、補正データ（図７参照）の出現確率とを乗算することにより、セクタ内の移動機数を分割エリア毎に振り分け、セクタ内における移動機の分布（図１１参照）を算出する（ステップＳ２０４：分布算出ステップ）。分布情報記憶部６０９は、分布算出部６０８が算出したセクタ内における移動機の分布（図１１参照）を記憶する（ステップＳ２０５）。
【００６３】
次に、集計部６１０は、集計エリア対応表記憶部６１１に記憶された集計エリア対応表（図１２参照）に基づいて、分布情報記憶部６０９に記憶された分布情報（図１１参照）において、分割エリアがどの集計エリアに対応するかを識別（図１４参照）する（ステップＳ２０６）。そして、集計部６１０は、集計エリア毎の移動機数を集計情報として算出する（ステップＳ２０７）。集計情報記憶部６１２は、集計部６１０によって集計された集計情報（図１５）を記憶する（ステップＳ２０８）。
【００６４】
次に、第１実施形態に係る移動通信システム１０の作用効果について説明する。移動通信システム１０によれば、分割エリア内に存在する移動機１００の移動機数に基づいて、セクタ内における移動機１００の分布を算出するための補正データが算出される。そして、算出された補正データと、セクタ内の移動機数とに基づいて、セクタ内における移動機１００の分布が算出される。これにより、セクタ内において、より詳細な移動機１００の分布を精度良く算出することができ、さらに、算出された移動機１００の分布に基づき、集計エリア毎の移動機数を精度良く算出することができる。
【００６５】
例えば、詳細な位置情報であるＧＰＳ位置情報は、ＧＰＳ位置情報を用いる所定のアプリケーションが実行された場合でないと取得されないため、全ての移動機１００から詳細な位置情報を取得し、詳細な移動機１００の分布を得ることは困難である。これに対し、位置登録信号を用いることにより、セクタ単位ではあるが移動機１００の位置を把握することができ、移動機数を算出することができる。そこで、本実施形態における移動通信システム１０の補正データ算出部６０４は、詳細な位置情報であるＧＰＳ位置情報に基づいて、セクタ内における移動機１００の分布を算出するための補正データを算出する。算出された補正データと、セクタ毎の移動機数（位置登録により得られる移動機数）とを用いることにより、セクタ内におけるより詳細な移動機１００の分布を精度良く算出することができる。
【００６６】
また、補正データとして、分割エリア毎の出現確率を算出し、算出した出現確率とセクタ内の移動機数とに基づいて、セクタ内の移動機数を分割エリア毎に振り分ける。これにより、セクタ内における移動機１００の分布を分割エリア毎に精度良く算出することができる。
【００６７】
また、第２取得部６０２が、ＧＰＳ位置情報記憶部６０１に記憶されたＧＰＳ位置情報に基づいて移動機数を計数することにより、分割エリア毎の移動機数を容易に取得することができる。
【００６８】
また、第２取得部６０２は、ＧＰＳ位置情報記憶部６０１から所定の集計時間帯のＧＰＳ位置情報を抽出し、分割エリア毎の移動機数を算出する。第１取得部６０７は、位置登録情報記憶部６０６から所定の集計時間帯の位置登録情報を抽出し、セクタ毎の移動機数を算出する。これにより、所定の集計時間帯におけるセクタ内の移動機数と、所定の集計時間帯における分割エリア毎の移動機数に基づいて算出された補正データとが得られる。このように、セクタ内の移動機数と補正データとが時間的に整合しているため、これらのデータを用いることにより所定の集計時間帯におけるセクタ内のより詳細な移動機の分布をより精度良く算出することができる。
【００６９】
なお、第１実施形態では、第２取得部６０２がＧＰＳ位置情報から分割エリア内の移動機数を算出する際に、同一の移動機１００が複数のＧＰＳ位置情報を出力したときの重複を含む場合について説明した。これに替えて、例えば、一の移動機１００からＧＰＳ位置情報を複数受信した場合、受信したＧＰＳ位置情報に対して移動機数の計数についての重み付けを行い、該重み付けの結果に基づいて分割エリア毎に移動機数を計数しても良い。この重み付けは、例えば、ＧＰＳ位置情報の緯度・経度の情報より得られる位置に基づいて行うことができる。これにより、重み付けの結果に基づいて分割エリア毎に移動機数をより現実に則して取得することができる。
【００７０】
また、重み付けの他の例として、一の移動機１００が出力したＧＰＳ位置情報の個数（回数）に基づいて重み付けを行うこともできる。この場合の具体例について説明する。ここでは、補正データ算出部６０４が、ある週の平日の００：００〜０１：００の補正データを算出するものとし、このため第２取得部６０２が、ある週の平日の００：００〜０１：００における分割エリア毎の移動機数を計数する場合について説明する。また、ある一の移動機１００が、月曜日の００：００〜０１：００の間に５個のＧＰＳ位置情報を出力し、火曜日の００：００〜０１：００の間に４個のＧＰＳ位置情報を出力しているものとする。
【００７１】
この場合、第２取得部６０２は、月曜日の００：００〜０１：００の間において、一の移動機１００がＧＰＳ位置情報を出力した分割エリアについては、一の移動機１００が出力したＧＰＳ位置情報の信号１個に付き０．２台（１／ＧＰＳ位置情報の数５個＝０．２）として移動機１００の計数を行う。また、第２取得部６０２は、火曜日の００：００〜０１：００の間において、一の移動機１００がＧＰＳ位置情報を出力した分割エリアについては、一の移動機１００が出力したＧＰＳ位置情報の信号１個に付き０．２５台（１／ＧＰＳ位置情報の数４個＝０．２５）として移動機１００の計数を行う。このようにして、月〜金曜日毎に、各分割エリア毎の移動機１００の数を計数する。
【００７２】
そして、月〜金曜日毎の移動機１００の数を分割エリア毎に加算し、加算した数を５（月〜金曜日までの日数）で除算することにより、ある週の平日の００：００〜０１：００における分割エリア毎の移動機１００の数を算出することができる。セクタ毎の移動機１００の数についても、上記と同様の重み付けを行って算出する。補正データ算出部６０４は、セクタ毎の移動機数と分割エリア毎の移動機数とより補正データを算出する。このように移動機１００が出力したＧＰＳ位置情報の信号の数に基づいて重み付けを行うことで、より精度の高い補正データを得ることができる。
【００７３】
なお、第１実施形態において、複数の移動機１００として、異なる周波数帯（例えば、２ＧＨｚ、１．７ＧＨｚ、８００ＭＨｚ）を利用するものが混在していてもよい。この場合であっても、セクタ毎に補正データを求め、セクタ毎に移動機の分布を求めることができる。また、分割エリア毎の移動機の分布に基づいて集計エリア内の移動機数を合算することで、集計エリア内の全移動機数を求めることができる。
【００７４】
ここで、周波数帯が異なると、周波数帯毎にＢＴＳ２００からの電波の到達範囲（セクタの範囲）が異なる。また、ある分割エリアが複数の周波数帯によってカバーされている場合、ある分割エリアをカバーするセクタの識別子も周波数帯毎にそれぞれ異なる。
【００７５】
第２取得部６０２が補正データを算出する場合、上述したように、所定の集計時間帯に対応するＧＰＳ位置情報を抽出してセクタ毎に補正データを算出する。
【００７６】
第１取得部６０７は、上述したように、セクタ毎に、セクタ内の移動機数を取得する。
【００７７】
分布算出部６０８は、補正データ算出部６０４においてセクタ毎に求められた補正データと、第１取得部６０７においてセクタ毎に求められた移動機数と、に基づいて、分割エリア毎の移動機数（図１８参照）を算出する。図１８は、分布算出部６０８において算出された分割エリア毎の移動機数を示すデータテーブルである。なお、図１８において、セクタ識別子が数字(セクタ１、セクタ２、セクタ３、…)のデータは２ＧＨｚの周波数帯を利用する移動機１００についてのデータであり、セクタ識別子が英文字(セクタＡ、セクタＢ、セクタＣ、…)のデータは８００ＭＨｚの周波数帯を利用する移動機１００についてのデータであるものとする。また、図１８に示す、ある２つのデータにおいて、一方のデータのエリア識別子が数字であり他方のデータのエリア識別子が英文字であっても、即ち、移動機が利用する周波数帯が互いに異なるデータであっても、分割エリアの地理的な位置が同じである場合には、分割エリア識別子は同じとなっている。但し、分割エリアの地理的な位置が同じであっても、周波数帯毎に分割エリア識別子が異なっていてもよい。
【００７８】
分布算出部６０８によって算出された分割エリア毎の移動機数から、集計エリア内の移動機数を集計する手順は、図１９に示す概念図によって表される。図１９は、分割エリア毎の移動機数から集計エリア内の移動機数を集計する手順を示す概念図である。図１９（ａ）は、セクタと集計エリアｅとの位置関係を示しており、セクタ１（セクタ識別子が「セクタ１」のセクタ）の一部とセクタ２（セクタ識別子が「セクタ２」のセクタ）の一部とが集計エリアｅに含まれ、更に、セクタＡ（セクタ識別子が「セクタＡ」のセクタ）の一部とセクタＢ（セクタ識別子が「セクタＢ」のセクタ）の一部とが集計エリアｅに含まれている様子を示している。
【００７９】
集計部６１０は、分布算出部６０８によって算出された分割エリア毎の移動機数と、集計エリア対応表記憶部６１１に記憶された集計エリア対応表と、に基づいて、集計エリアｅに属する分割エリアの移動機数を分割エリア毎に抽出する（図１９（ｂ）参照）。なお、図１９（ａ）及び図１９（ｂ）では、移動機１００が利用する周波数帯毎に移動機数の分布を示しているが、集計部６１０としては、周波数帯を意識する必要はなく、集計エリアｅに含まれる分割エリアの移動機数を抽出すればよい。
【００８０】
集計部６１０は、抽出した分割エリア毎の移動機数のうち、分割エリア識別子が同じ分割エリアの移動機数を合算する（図１９（ｃ）参照）。この図１９（ｃ）に示す移動機数が、集計エリアｅ内における移動機の分布となる。また、集計部６１０が、分割エリア毎に求められている集計エリアｅ内の移動機数を合算することで、集計エリアｅ内における全移動機数を求めることができる（図１９（ｄ）参照）。また、集計部６１０が集計エリアｅ内の分割エリア毎の移動機数を合算する際に、例えば、セクタＡ，Ｂについての分割エリア内の移動機数（図１９（ｂ）の下段に示す分割エリア内の移動機数）のみを合算することで、集計エリアｅ内における８００ＭＨｚを利用する移動機１００の台数を求めることができる。
【００８１】
なお、集計エリアｅは、図１９（ａ）及び図２０に示すようにセクタ間を跨って設定されていてもよい。集計エリアｅ内の移動機数を求める場合、上述のように、分割エリア毎の移動機数が分かっていれば求めることができる。このため、集計エリアがセクタ間を跨っていても集計エリアに対応する分割エリア内の移動機数を用いることで、セクタ間の跨りに関係なく、集計エリア内の移動機数や分布を求めることができる。また、図２０に示すように、集計エリア内に少なくとも一部が含まれることとなるセクタの数が、周波数帯毎に異なっていても良い。例えば、図２０に示す例において、集計エリアｅは、２ＧＨｚを利用する移動機１００をカバーするセクタのうち３つのセクタ１〜３に跨り、８００ＭＨｚを利用する移動機１００をカバーするセクタのうち２つのセクタＢ，Ｃに跨っている。
【００８２】
また、上記では、移動機１００が利用する周波数帯として、２ＧＨｚ、１．７ＧＨｚ、８００ＭＨｚを例に示したが、これ以外にも、例えば屋内に設置される比較的狭いエリアを通信範囲とする屋内局が使用する周波数帯を含めても良い。
【００８３】
次に、第１実施形態の変形例について説明する。本変形例は、第１実施形態における情報分析装置６００の構成を一部変更したものであり、第１実施形態から変更された部分についてのみ説明する。また、第１実施形態と同一構成要素については同一番号を付して説明を省略する。
【００８４】
図２は、移動通信システム１０Ａの機能構成を示すブロック図である。図２に示すように、移動通信システム１０Ａは、移動機１００と、ＢＴＳ２００と、ＲＮＣ３００と、交換機４００と、情報分析装置６００Ａとを含んで構成されている。
【００８５】
情報分析装置６００Ａは、移動機のユーザの属性情報を予め記憶する属性情報記憶部（属性情報記憶手段）６１５Ａを備える。属性情報記憶部６１５Ａは、ユーザの属性情報として、例えば、年齢層、性別、及び住所情報等と、移動機１００のユーザ識別子と、を対応付けて予め記憶している。情報分析装置６００Ａにおいて、属性情報記憶部６１５Ａ以外の構成要素は第１実施形態と同じであり、説明を省略する。
【００８６】
次に、情報分析装置６００Ａで行われる処理について説明する。第２取得部６０２は、ＧＰＳ位置情報記憶部６０１からＧＰＳ位置情報を抽出するとともに属性情報記憶部６１５Ａから属性情報を取得し、属性情報のユーザ識別子とＧＰＳ位置情報のユーザ識別子とを照合することにより、各ＧＰＳ位置情報毎のユーザの属性を把握する。そして、第２取得部６０２は、把握した各ＧＰＳ位置情報毎のユーザの属性に基づいて、属性毎及び分割エリア毎に移動機数を計数する。補正データ算出部６０４は、属性毎及び分割エリア毎の補正データを算出する。
【００８７】
第１取得部６０７は、位置登録情報記憶部６０６から位置登録情報を抽出するとともに属性情報記憶部６１５Ａから属性情報を取得し、属性情報のユーザ識別子と位置登録情報のユーザ識別子とを照合することにより、各位置登録情報毎のユーザの属性を把握する。そして、第１取得部６０７は、把握した各位置登録情報毎のユーザの属性に基づいて、属性毎及びセクタ毎に移動機数を計数する。分布算出部６０８は、セクタ内の移動機１００の分布を属性毎に算出する。集計部６１０は、属性毎に集計エリア毎の移動機１００の分布を算出する。
【００８８】
このように本変形例によれば、第１実施形態における効果に加え、ユーザの属性（例えば、年齢層、性別、住所情報等）毎に移動機１００の詳細な分布を得ることができ、データを有効利用することが可能となる。
【００８９】
［第２実施形態］
次に、本発明に係る移動機分布算出システム及び移動機分布算出方法を適用した移動通信システムの好適な第２実施形態について説明する。第２実施形態は、第１実施形態における情報分析装置６００が補正データを算出する際の処理構成を変更したものである。また、第１実施形態と同一構成要素については同一番号を付して説明を省略する。
【００９０】
第２実施形態における移動通信システム（１０Ｂ）の概略構成は、第１実施形態において図１を用いて説明した移動通信システム１０と同様であり、説明を省略する。
【００９１】
次に、第２実施形態に係る移動通信システムの機能ブロック構成について説明する。図２１は、移動通信システム１０Ｂの機能構成を示すブロック図である。移動通信システム１０Ｂは、移動機１００と、ＢＴＳ２００と、ＲＮＣ３００Ｂと、交換機４００と、情報分析装置６００Ｂとを含んで構成されている。
【００９２】
ＲＮＣ３００Ｂは、通信制御部３０１と、セクタ識別部３０２Ｂとを含んで構成されている。セクタ識別部３０２Ｂは、移動機１００が位置登録を行った際に当該移動機１００が在圏するセクタを識別し、識別結果を交換機４００へ出力する。なお、第２実施形態におけるセクタ識別部３０２Ｂは、第１実施形態におけるセクタ識別部３０２とは異なり、移動機１００がＧＰＳ位置情報を交換機４００に出力したときに当該移動機１００が在圏するセクタの識別処理を行わない。
【００９３】
情報分析装置６００Ｂは、ＧＰＳ位置情報記憶部６０１Ｂと、第２取得部６０２Ｂと、結合位置情報記憶部６１６Ｂと、分割エリア対応表記憶部６０３と、補正データ算出部６０４Ｂと、補正データ記憶部６０５と、位置登録情報記憶部６０６と、第１取得部６０７と、分布算出部６０８と、分布情報記憶部６０９と、集計部６１０と、集計エリア対応表記憶部６１１と、集計情報記憶部６１２とを含んで構成されている。
【００９４】
ＧＰＳ位置情報記憶部６０１Ｂは、交換機４００の位置情報管理部４０２から出力されたＧＰＳ位置情報を記憶するものである。具体的には、ＧＰＳ位置情報記憶部６０１Ｂは、図２２に示すように、ユーザ識別子と、移動機１００によってＧＰＳ位置情報が取得された時刻と、緯度・経度とが対応付けられた情報を、ＧＰＳ位置情報が取得された時刻毎に記憶している。
【００９５】
第２取得部６０２Ｂは、ＧＰＳ位置情報記憶部６０１Ｂに記憶されたＧＰＳ位置情報（図２２参照）と、位置登録情報記憶部６０６に記憶された位置登録情報（図８参照）とを取得する。そして、第２取得部６０２Ｂは、取得したＧＰＳ位置情報及び位置登録情報のうち、ユーザ識別子が同一、かつＧＰＳ位置情報の時刻と位置登録情報の時刻との差が所定時間以内（第２実施形態では、例えば、１分以内とする。）の情報を抽出する。第２取得部６０２Ｂは、時刻の差が所定時間以内のＧＰＳ位置情報及び位置登録情報を結合して結合位置情報を算出する。
【００９６】
結合位置情報記憶部６１６Ｂは、第２取得部６０２Ｂが算出した結合位置情報を、図４に示すように、ユーザ識別子と、時刻と、セクタ識別子と、緯度・経度とを対応付けて記憶する。なお、結合位置情報（図４参照）の時刻の情報は、ＧＰＳ位置情報の時刻及び位置登録情報の時刻に差がある場合、例えば、ＧＰＳ位置情報の時刻又は位置登録情報の時刻等を用いる（第２実施形態では、ＧＰＳ位置情報の時刻を用いるものとする。）。
【００９７】
次に、第２取得部６０２Ｂは、結合位置情報記憶部６１６Ｂに記憶された結合位置情報を用いて、分割エリア毎の移動機数を取得する。具体的には、第２取得部６０２Ｂは、結合位置情報記憶部６１６Ｂに記憶された結合位置情報（図４参照）のうち、所定の抽出条件（例えば、過去数ヶ月間の平日の００：００〜００：５９）に対応する結合位置情報を抽出する。そして、第２取得部６０２Ｂは、抽出した結合位置情報に含まれる緯度・経度と、分割エリア対応表記憶部６０３に記憶された分割エリア対応表（図５参照）とを用い、抽出した結合位置情報のそれぞれがどの分割エリアに対応するかを照合する。この照合の結果、第２取得部６０２Ｂは、図６に示すように、ユーザ識別子と、時刻と、セクタ識別子と、分割エリア識別子とが対応付けられた情報を得る。
【００９８】
そして、第２取得部６０２Ｂは、第１実施形態と同様に、図６に示す情報に基づいて、セクタ識別子が同一のレコード数Ａと、セクタ識別子が同一かつ分割エリア識別子が同一のレコード数Ｂとを計数する。第２取得部６０２Ｂは、抽出条件（過去数ヶ月間の集計時間帯）と、レコード数Ａと、当該レコード数Ａに対応するセクタ識別子と、レコード数Ｂと、当該レコード数Ｂに対応するセクタ識別子及び分割エリア識別子とを、それぞれ対応付けて補正データ算出部６０４Ｂへ出力する。
【００９９】
補正データ算出部６０４Ｂは、第１実施形態と同様に、第２取得部６０２Ｂから出力された情報に基づいて、分割エリア毎の出現確率を含む補正データ（図７参照）を算出する。
【０１００】
なお、第１取得部６０７がセクタ内の移動機数を取得する処理、分布算出部６０８が補正データ記憶部６０５に記憶された補正データを用いて、セクタ内における移動機１００の分布を算出する処理、集計部６１０が集計エリア内の移動機数を集計する処理は、第１実施形態と同様であるため、説明を省略する。
【０１０１】
次に、情報分析装置６００Ｂで行われる補正データ算出処理の流れについて説明する。図２３は、補正データを算出するための処理の流れを示すフローチャートである。なお、本処理は、例えば１週間に１回等、定期的に行われるものである。
【０１０２】
まず、第２取得部６０２Ｂは、ＧＰＳ位置情報記憶部６０１Ｂに記憶されたＧＰＳ位置情報（図２２参照）を取得し（ステップＳ３０１）、さらに、位置登録情報記憶部６０６に記憶された位置登録情報（図８参照）を取得する（ステップＳ３０２）。そして、第２取得部６０２Ｂは、時刻の差が所定時間以内のＧＰＳ位置情報及び位置登録情報を結合し結合位置情報を算出する（ステップＳ３０３）。結合位置情報記憶部６１６Ｂは、第２取得部６０２によって算出された結合位置情報（図４参照）を記憶する（ステップＳ３０４）。
【０１０３】
次に、第２取得部６０２Ｂは、結合位置情報記憶部６１６Ｂに記憶された結合位置情報を、抽出条件に従って抽出する（ステップＳ３０５）。そして、第２取得部６０２Ｂは、分割エリア対応表記憶部６０３に記憶された分割エリア対応表（図５参照）に基づいて、抽出した結合位置情報のそれぞれがどの分割エリアに対応するかを照合（図６参照）する（ステップＳ３０６）。その後、第２取得部６０２Ｂは、各セクタ毎の移動機数（レコード数Ａ）と、各分割エリア毎の移動機数（レコード数Ｂ）とを計数する（ステップＳ３０７：第２取得ステップ）。
【０１０４】
次に、補正データ算出部６０４Ｂは、補正データとして、レコード数Ｂをレコード数Ａで除算することにより各分割エリア毎の出現確率を補正データ（図７参照）として算出する（ステップＳ３０８：補正データ算出ステップ）。そして、補正データ記憶部６０５は、補正データ算出部６０４Ｂによって算出された補正データを記憶する（ステップＳ３０９）。
【０１０５】
なお、情報分析装置６００Ｂで行われる集計エリア内の移動機数を算出する処理の流れは、第１実施形態と同様（図１７のフローチャート参照）であるため、説明を省略する。
【０１０６】
次に、第２実施形態に係る移動通信システム１０Ｂの作用効果について説明する。移動通信システム１０Ｂによれば、分割エリア内に存在する移動機１００の移動機数に基づいて、セクタ内における移動機１００の分布を算出するための補正データが算出される。そして、算出された補正データと、セクタ内の移動機数とに基づいて、セクタ内における移動機１００の分布が算出される。これにより、セクタ内において、より詳細な移動機１００の分布を精度良く算出することができ、さらに、算出された移動機１００の分布に基づき、集計エリア毎の移動機数を精度良く算出することができる。
【０１０７】
また、第２取得部６０２Ｂにおいて、ＧＰＳ位置情報及び位置登録情報のうち、ユーザ識別子が同一、かつＧＰＳ位置情報の時刻と位置登録情報の時刻との差が所定時間以内の情報が抽出される。ここで、例えば、移動機１００が位置登録を行ったタイミングと、移動機１００がＧＰＳ位置情報を取得したタイミングとが時間的に近い（例えば、１分等）場合、位置登録を行った時と、移動機がＧＰＳ位置情報を取得した時とにおいて移動機１００はほぼ同じ位置にいる（ほとんど移動していない）と考えられる。このため、両タイミングの時間差が所定時間以内のＧＰＳ位置情報を用いることにより、位置登録が行われた時に実際に移動機１００が存在していた分割エリアを把握し、把握結果に基づいて分割エリア毎に移動機数を取得することができる。これにより、ＧＰＳ位置情報に基づいて算出される補正データが、位置登録信号に基づいて取得された移動機数の情報に対して時間的により正確に対応するものとなり、セクタ内の移動機１００の分布をより精度良く算出することができる。なお、第２実施形態において、複数の移動機１００として異なる周波数帯を利用するものが混在している場合であっても、第１実施形態のなお書きで説明した処理と同様の処理により、集計エリア内の全移動機数等を算出することができる。
【０１０８】
［第３実施形態］
第３実施形態では、第１取得部及び第２取得部において取得される移動機数を、特徴量を用いて算出する実施形態について説明する。なお、「特徴量」とは、移動機によって生成された位置登録情報又はＧＰＳ位置情報についての推定生成密度に対応する情報である。また、ここでの「推定生成密度」とは、当該位置登録情報又はＧＰＳ位置情報を生成した移動機が、位置登録情報又はＧＰＳ位置情報が生成された時刻周辺で単位時間あたりに生成する位置登録情報又はＧＰＳ位置情報の数の推定値を意味する。
【０１０９】
図２４に示すように、第３実施形態における情報分析装置６００Ｃの機能ブロック構成は、第１実施形態における情報分析装置６００の機能ブロック構成（図２）において、第１取得部６０７及び第２取得部６０２に代えて第１取得部６０７Ｃ及び第２取得部６０２Ｃを備えた構成とされている。第１取得部６０７Ｃ及び第２取得部６０２Ｃ以外の構成は第１実施形態で説明した構成と同様であり、同一番号を付して説明を省略する。
【０１１０】
まず、第２取得部６０２Ｃの詳細について説明する。図２５は、第２取得部の機能構成を示すブロック図である。図２５に示すように第２取得部６０２Ｃは、集計対象取得部６２１、前後位置情報取得部（第２前後位置情報取得手段）６２２、特徴量計算部（第２特徴量計算手段）６２３、分割エリア照合部６２４、及び移動機数推計部（第２移動機数推計手段）６２５を備えている。
【０１１１】
集計対象取得部６２１は、ＧＰＳ位置情報記憶部（第２位置情報記憶手段）６０１から、ＧＰＳ位置情報が取得された時刻（位置取得時刻情報）が集計時間帯の集計開始時刻以降であり且つ集計終了時刻以前である１又は複数のＧＰＳ位置情報を、集計対象ＧＰＳ位置情報として取得する。なお、集計時間帯は、第１実施形態で述べたように集計部６１０において移動機１００の詳細な分布を得ようとする時間帯である。
【０１１２】
前後位置情報取得部６２２は、特徴量を求める対象の集計対象ＧＰＳ位置情報（以下「第１のＧＰＳ位置情報」という）について、当該第１のＧＰＳ位置情報と同一のユーザ識別子を含むＧＰＳ位置情報のうち、当該第１のＧＰＳ位置情報の直前のＧＰＳ位置情報（以下「第２のＧＰＳ位置情報」という）が取得された時刻、及び当該第１のＧＰＳ位置情報の直後のＧＰＳ位置情報（以下「第３のＧＰＳ位置情報」という）が取得された時刻を取得する。なお、前後位置情報取得部６２２は、第２又は第３のＧＰＳ位置情報の全体を取得することは必須ではなく、少なくとも、ＧＰＳ位置情報に含まれる時刻についての情報を取得すればよい。
【０１１３】
特徴量計算部６２３は、第１のＧＰＳ位置情報それぞれについての特徴量を計算する。例えば、特徴量計算部６２３は、第２のＧＰＳ位置情報が取得された時刻と第３のＧＰＳ位置情報が取得された時刻との差を、当該第１のＧＰＳ位置情報についての特徴量として計算する。また、特徴量計算部６２３は、第２のＧＰＳ位置情報が取得された時刻が異常値である場合、ここでは一例として第１のＧＰＳ位置情報が取得された刻と第２のＧＰＳ位置情報が取得された時刻との差が所定の基準値（例えば１時間）より大きい場合に、第１のＧＰＳ位置情報が取得された時刻から予め定められた時間（例えば１時間）だけ過去に遡った時刻を第２のＧＰＳ位置情報が取得された時刻として用いて、第１のＧＰＳ位置情報についての特徴量を計算する。同様に、特徴量計算部６２３は、第３のＧＰＳ位置情報が取得された時刻が異常値である場合、ここでは一例として第１のＧＰＳ位置情報が取得された時刻と第３のＧＰＳ位置情報が取得された時刻との差が所定の基準値（例えば１時間）より大きい場合に、第１のＧＰＳ位置情報が取得された時刻から予め定められた時間（例えば１時間）だけ未来に進めた時刻を第３のＧＰＳ位置情報が取得された時刻として用いて、第１のＧＰＳ位置情報についての特徴量を計算する。このような第２、第３のＧＰＳ位置情報が取得された時刻が異常値である場合の処理は、必須の処理ではないが、上記処理を行うことで、移動機１００が圏外に位置していることや移動機１００の電源がオフされていること等に起因してＧＰＳ位置情報の取得時間間隔が異常に長くなった際に、当該異常に長くなった取得時間間隔による影響が過大に出ることを防ぐことができる。また、特徴量計算部６２３は、図２６に示すように、ＧＰＳ位置情報（第１のＧＰＳ位置情報）と、算出した特徴量と、を対応付けたデータベースを作成する。
【０１１４】
分割エリア照合部６２４は、特徴量が対応付けされたＧＰＳ位置情報（図２６参照）に含まれる緯度・経度と、分割エリア対応表記憶部６０３に記憶された分割エリア対応表（図５参照）とを用い、特徴量が対応付けされたＧＰＳ位置情報のそれぞれがどの分割エリアに対応するかを照合する。分割エリア照合部６２４は、照合の結果として、図２７に示すように、ユーザ識別子と、時刻と、セクタ識別子と、分割エリア識別子と、特徴量と、を対応付けたデータベースを作成する。
【０１１５】
移動機数推計部６２５は、ＧＰＳ位置情報についての特徴量、及び集計開始時刻と集計終了時刻との差である集計時間帯の時間長に基づいて、移動機数を推計する。詳細には、移動機数推計部６２５は、分割エリア照合部６２４によって対応付けされた結果（図２７参照）を用い、セクタ識別子が同一のＧＰＳ位置情報を抽出し、抽出したＧＰＳ位置情報に対応付けられた特徴量に基づいて移動機数Ａを推計し、セクタ識別子が同一かつ分割エリア識別子が同一のＧＰＳ位置情報を抽出し、抽出したＧＰＳ位置情報に対応付けられた特徴量に基づいて移動機数Ｂを推計する。この移動機数Ａが、セクタ内の移動機１００の移動機数となり、移動機数Ｂが、分割エリア内の移動機１００の移動機数となる。
【０１１６】
詳細は後述するが、移動機数推計部６２５は、移動機数Ａを推計する際に、図２７に示すデータベースにおいて、セクタ識別子が同一のＧＰＳ位置情報の特徴量の総和を集計時間帯の時間長の２倍によって除して得られた数値を移動機数Ａとして推計する。同様に、移動機数推計部６２５は、移動機数Ｂを推計する際に、図２７に示すデータベースにおいて、セクタ識別子が同一かつ分割エリア識別子が同一のＧＰＳ位置情報の特徴量の総和を集計時間帯の時間長の２倍によって除して得られた数値を移動機数Ｂとして推計する。
【０１１７】
［移動機数推計の考え方および計算方法］
ここで、移動機数推計の考え方および計算方法を説明する。図２８に示すモデルのように、ある集計時間帯（長さＴ）の間に、ｎ個の移動機ａ１，ａ２，…，ａｎがセクタＳを通過し、各移動機ａｉの集計時間帯内のセクタＳの滞在時間がｔｉ（０＜ｔｉ≦Ｔ）であったとする。このとき、セクタＳに存在する移動機数ｍ（実際にはセクタＳに存在する移動機数ｍの集計時間帯内における平均値）は、以下の式（１）で表わされる。
【数１】
2012054921.tif 000003
即ち、各移動機ａｉの集計時間帯内のセクタＳの滞在時間ｔｉの総和を集計時間帯の長さＴで除した結果を、移動機数ｍとして推計する。ただし、移動機ａｉの集計時間帯内のセクタＳの滞在時間ｔｉの真の値は観測不能であるが、各移動機ａｉはＧＰＳ位置情報を発信し、それらのＧＰＳ位置情報は観測可能である。
【０１１８】
移動機ａｉが集計時間帯内にセクタＳで発信したＧＰＳ位置情報を、時刻順に
【数２】
2012054921.tif 000004
（ｘｉは、移動機ａｉが集計時間帯内にセクタＳで発信したＧＰＳ位置情報の総数）とすると、移動機数の推計とは、観測されたＧＰＳ位置情報ｑｉｊ（ｊは１以上ｘｉ以下の整数）からｍの値を推計することに他ならない。
【０１１９】
さて、図２９に基づき移動機数推計の計算方法を説明する。移動機ａｉからＧＰＳ位置情報ｑｉｊが送信される密度（即ち、単位時間あたりのＧＰＳ位置情報数）をｐｉとする。このとき、ＧＰＳ位置情報が送信される確率がセクタに対して独立であれば、移動機ａｉが集計時間帯内にセクタＳで発信したＧＰＳ位置情報の総数ｘｉの期待値Ｅ（ｘｉ）は、Ｅ（ｘｉ）＝ｔｉ×ｐｉであるため、移動機ａｉの集計時間帯内のセクタＳの滞在時間ｔｉの期待値Ｅ（ｔｉ）について以下の式（２）が成立する。
Ｅ（ｔｉ）＝ｘｉ／ｐｉ  （２）
ここで、ＧＰＳ位置情報ｑｉｊの送信時刻をｕｉｊとしたとき、ＧＰＳ位置情報ｑｉｊの密度ｐｉｊは、以下の式（３）で与えられる。
ｐｉｊ＝２／（ｕｉ（ｊ＋１）−ｕｉ（ｊ−１））  （３）
ここで、ＧＰＳ位置情報ｑｉｊを第１のＧＰＳ位置情報とすると、ＧＰＳ位置情報ｑｉ（ｊ−１）は第２のＧＰＳ位置情報、ＧＰＳ位置情報ｑｉ（ｊ＋１）は第３のＧＰＳ位置情報に相当する。本実施形態では、第２のＧＰＳ位置情報ｑｉ（ｊ−１）の送信時刻ｕｉ（ｊ−１）と第３のＧＰＳ位置情報ｑｉ（ｊ＋１）の送信時刻ｕｉ（ｊ＋１）の差、即ち、上記式（３）の（ｕｉ（ｊ＋１）−ｕｉ（ｊ−１））を、第１のＧＰＳ位置情報についての特徴量ｗｉｊとする。そのため、上記式（３）は、以下となる。即ち、特徴量ｗｉｊは、密度ｐｉｊの逆数に対応づけて算出することができる。
ｐｉｊ＝２／（ｕｉ（ｊ＋１）−ｕｉ（ｊ−１））＝２／ｗｉｊ  （４）
このとき密度ｐｉは、
【数３】
2012054921.tif 000005
で与えられるため、移動機数ｍの推計値Ｅ（ｍ）は以下の式（６）で計算することができる。
【数４】
2012054921.tif 000006
図２９の例に示すように、集計時間帯内であり且つ移動機ａｉがセクタＳに滞在していた期間内に、移動機ａｉはＧＰＳ位置情報ｑｉ１、ｑｉ２、ｑｉ３を送信し、ＧＰＳ位置情報ｑｉ１の直前にＧＰＳ位置情報ｑｉ０を、ＧＰＳ位置情報ｑｉ３の直後にＧＰＳ位置情報ｑｉ４を送信したものとし、ＧＰＳ位置情報ｑｉ０、ｑｉ１、ｑｉ２、ｑｉ３、ｑｉ４の送信時刻をそれぞれｕｉ０、ｕｉ１、ｕｉ２、ｕｉ３、ｕｉ４とすると、上記の考え方は、移動機ａｉの集計時間帯内のセクタＳの滞在時間ｔｉを、（ｕｉ０とｕｉ１の中点）から（ｕｉ３とｕｉ４の中点）までの期間と推計することに相当する。なお、移動機ａｉは、集計時間帯内ではないものの、セクタＳへの滞在中にＧＰＳ位置情報ｑｉ４を送信している。但し、滞在時間ｔｉの推計量の不偏性を維持するために、ここでは一例として、滞在時間ｔｉの終了時刻を集計時間帯Ｔの終了時刻と同じとして推計することは行わない処理を説明する。
【０１２０】
図２５に戻り、移動機数Ａ，Ｂの推計後、移動機数推計部６２５は、抽出条件（集計時間帯の情報）と、移動機数Ａと、当該移動機数Ａに対応するセクタ識別子と、を対応付け、抽出条件と、移動機数Ｂと、当該移動機数Ｂに対応するセクタ識別子及び分割エリア識別子とを、対応付けて、それぞれ補正データ算出部６０４へ出力する。
【０１２１】
補正データ算出部６０４は、第１実施形態においてレコード数Ａとレコード数Ｂとを用いて補正データとしての出現確率を求めた方法と同様に、移動機数Ｂを、当該移動機数Ｂに対応付けられた分割エリアを含むセクタの移動機数Ａで除算することにより、出現確率を算出する。
【０１２２】
次に、第１取得部６０７Ｃの詳細について説明する。図３０は、第１取得部の機能構成を示すブロック図である。図３０に示すように第１取得部６０７Ｃは、集計対象取得部６７１、前後位置情報取得部（第１前後位置情報取得手段）６７２、特徴量計算部（特徴量計算手段）６７３、及び移動機数推計部（移動機数推計手段）６７５を備えている。
【０１２３】
集計対象取得部６７１は、位置登録情報記憶部（第１位置情報記憶手段）６０６から、位置登録信号が取得された時刻（位置取得時刻情報）が集計時間帯の集計開始時刻以降であり且つ集計終了時刻以前である１又は複数の位置登録情報を、集計対象位置登録情報として取得する。なお、集計時間帯は、上述したように集計部６１０において移動機１００の詳細な分布を得ようとする時間帯である。
【０１２４】
前後位置情報取得部６７２は、特徴量を求める対象の集計対象位置登録情報（以下「第１の位置登録情報」という）について、当該第１の位置登録情報と同一のユーザ識別子を含む位置登録情報のうち、当該第１の位置登録情報の直前の位置登録情報（以下「第２の位置登録情報」という）が取得された時刻、及び当該第１の位置登録情報の直後の位置登録情報（以下「第３の位置登録情報」という）が取得された時刻を取得する。なお、前後位置情報取得部６７２は、第２又は第３の位置登録情報の全体を取得することは必須ではなく、少なくとも、位置登録情報に含まれる時刻についての情報を取得すればよい。
【０１２５】
特徴量計算部６７３は、第１の位置登録情報それぞれについての特徴量を計算する。この特徴量の算出は、上述した特徴量計算部６２３で行われる処理と同様の処理によって行うことができる。特徴量計算部６７３は、図３１に示すように、位置登録情報（第１の位置登録情報）と、算出した特徴量と、を対応付けたデータベースを作成する。
【０１２６】
移動機数推計部６７５は、位置登録情報についての特徴量、及び集計開始時刻と集計終了時刻との差である集計時間帯の時間長に基づいて、セクタ毎の移動機数を推計する。詳細には、移動機数推計部６７５は、特徴量計算部６７３による対応付け結果（図３１参照）を用い、セクタ識別子が同一の位置登録情報を抽出し、抽出した位置登録情報の特徴量に基づいて移動機数を推計することで、セクタ毎の移動機数を推計する。この特徴量を用いた移動機数の推計は、上述の移動機数推計部６２５で行われる処理と同様の処理によって行うことができる。移動機数の推計後、移動機数推計部６７５は、図９に示すように、種別（集計時間帯の平日又は休日の種別）と、集計時間帯と、セクタ識別子と、セクタ毎の移動機数とを対応付けし、対応付けした結果を分布算出部６０８へ出力する。
【０１２７】
以上のように、第２取得部６０２ＣがＧＰＳ位置情報についての特徴量に基づいて移動機数Ａ，Ｂを算出し、算出された移動機数Ａ，Ｂに基づいて補正データ算出部６０４が補正データを算出する。また、第１取得部６０７Ｃが、位置登録情報についての特徴量に基づいてセクタ毎の移動機数を算出する。この場合には、特徴量を用いて移動機数をより正確に算出することができ、移動機の分布を精度良く算出することができる。
【０１２８】
また、集計エリア内の移動機数から人口を算出することにより、集計エリア内の人口を精度良く得ることができる。具体的には、集計部６１０が、集計エリア内の移動機数を求めると共に、求めた移動機数に拡大係数を乗算することで、集計エリア内の人口を得る。また、集計部６１０において、分割エリア毎の移動機数に拡大係数を乗算することで、分割エリア毎の人口分布を求めてもよい。この拡大係数は、算出された移動機数を人口に変換するための係数であり、予め定めて記憶しておく以外に、以下のようにして求めてもよい。拡大係数は、一例として、「在圏率」と「移動機の普及率」との積（即ち、人口に対する在圏数の比率）の逆数を用いることができる。ここで「在圏率」とは、移動機の契約台数に対する在圏数の比率を意味し、「移動機の普及率」とは人口に対する移動機の契約台数の比率を意味する。
【０１２９】
また、拡大係数は、例えば、以下のように特徴量及び集計時間帯の時間長に基づいて推計された移動機数（在圏数）を用いて導出してもよい。即ち、位置登録情報から特徴量を求め、特徴量及び集計時間帯の時間長に基づいて、あるエリア内の人口推計単位毎の移動機数を集計することでユーザ数ピラミッドデータを得るともに、統計データ（例えば住民基本台帳など）として予め求められた同じエリアにおける人口ピラミッドデータを取得する。そして、ユーザ数ピラミッドデータ及び人口ピラミッドデータにおいて人口推計単位毎の位置登録情報の取得率（即ち、在圏数／人口）を算出する。ここで得られた「位置登録情報の取得率（即ち、在圏数／人口）」が、前述した「在圏率と移動機の普及率との積」に相当する。このようにして得られた「位置登録情報の取得率」の逆数を拡大係数として導出することができる。なお、拡大係数を求める人口推計単位としては、一例として、住所の都道府県毎、時間帯として１時間毎等を採用してもよい。
【０１３０】
このように拡大係数を用いることで集計エリア毎の人口や人口分布を求めることができ、求めた人口や人口分布を、例えば、移動機１００のアンテナを設置する際の回線容量の基準として用いたりすることができる。なお、本明細書で説明する他の実施形態及び変形例においても、拡大係数を用いることで移動機数から人口を求めることができる。
【０１３１】
また、従来のように、国勢調査や電話帳データを用いて詳細な広範囲の人口分布から詳細な人口分布を得ようとする場合、国勢調査の世帯データや電話帳データは、曜日や時間帯毎に集計されたものではないため、例えば、休日の午後のある時間帯の人口分布を得たいという要求には対応できない。これに対し、本実施形態における移動通信システムでは、所望の時間帯や曜日等に対応する補正データやセクタ毎の移動機数を算出することにより、所望の時間帯や曜日等に対応するより詳細な人口分布を精度良く得ることができる。
【０１３２】
［第４実施形態］
第４実施形態は、第３実施形態の第１取得部６０７Ｃにおいて行われるセクタ毎の移動機数算出処理の内容を変更したものである。
【０１３３】
図３２に示すように、第４実施形態における情報分析装置６００Ｄの機能ブロック構成は、第１実施形態における情報分析装置６００の機能ブロック構成（図２）において、第１取得部６０７に代えて第１取得部６０７Ｄを備えた構成とされている。第１取得部６０７Ｄ以外の構成は第１実施形態で説明した構成と同様であり、同一番号を付して説明を省略する。
【０１３４】
第１取得部６０７Ｄの詳細について説明する。図３３は、第１取得部の機能構成を示すブロック図である。図３３に示すように第１取得部６０７Ｄは、抽出部（抽出手段）６８１、及び入込端末数推計部（移動機数推計手段）６８２を備える。
【０１３５】
抽出部６８１は、位置登録情報記憶部（第１位置情報記憶手段）６０６に記憶された位置登録情報から、位置登録情報に含まれるセクタＩＤ毎に、位置登録が行われた時刻（位置取得時刻情報）が集計時間帯内にある位置登録情報を移動機毎に１つずつ抽出する。又は、抽出部６８１は、位置登録情報に含まれるセクタＩＤ毎に、位置登録が行われた時刻が集計時間帯内にある位置登録情報に対応する移動機を抽出する。その具体的な抽出方法は後述する。抽出部６８１は、詳しくは後述する抽出方法により、集計時間帯の少なくとも一部において、セクタＩＤ毎に当該セクタに滞在したと推定される移動機、又は、当該移動機にて集計時間帯内又は後述する拡張期間内に生成された１つの位置登録情報を抽出する機能を実現する。なお、集計時間帯は、第１実施形態で述べたように集計部６１０において移動機１００の詳細な分布を得ようとする時間帯である。
【０１３６】
入込端末数推計部６８２は、抽出部６８１により抽出された位置登録情報の数又は移動機の数をセクタＩＤ毎にカウントし、得られたカウント数を集計時間帯内における各セクタの入込端末数として推計する。なお、本明細書において、入込端末数とは、集計時間帯の少なくとも一部において集計対象とするエリア（セクタ）に滞在したユニークな移動機数を意味する。また、ここでの「ユニークな」とは、同一移動機の重複カウントを除いた後の数を意味する。
【０１３７】
［入込端末数推計処理］
以下、第１取得部６０７Ｄで行われる入込端末数推計処理の一例を説明する。この一例においては、セクタ内の各移動機１００の推定滞在期間を用いて入込端末数を求める。
【０１３８】
抽出部６８１は、ユーザ識別子が同じ位置登録情報において、位置登録情報が取得された時刻が後述の拡張期間内にあり、セクタＩＤがある所定のセクタを示すセクタ内位置データ、及び位置登録情報が取得された時刻順に時系列に並べた場合に当該セクタ内位置データに隣接しセクタＩＤがある所定のセクタ外を示すセクタ外位置データに基づいて、移動機がある所定のセクタ内に滞在していた推定滞在期間を移動機毎に算出する。なお、上記の「拡張期間」は、ここでは一例として、集計時間帯を前後に所定の時間幅（例えば１時間）だけ拡張した期間、即ち、集計開始時刻ｔ０から所定時間だけ過去に遡った時刻を開始点とし、集計終了時刻ｔ１から所定時間だけ未来に進めた時刻を終了点とする時間帯を意味する。
【０１３９】
図３４に示すように、拡張期間内に取得された同一移動機に関する一連の位置登録情報を当該位置登録情報が取得された時刻順に時系列に並べ、位置登録情報のセクタＩＤがある所定のセクタを示すセクタ内位置データを黒丸で、セクタ内位置データに隣接しセクタＩＤがある所定のセクタ外を示すセクタ外位置データを白丸で、それぞれ示した場合、時系列上で最先のセクタ内位置データの取得時刻ｔａと、該最先のセクタ内位置データに隣接するセクタ外位置データの取得時刻ｔｘとの按分点（一例として中点）に相当する時刻ｔｉｎを推定滞在期間の開始時刻とする。同様に、時系列上で最後のセクタ内位置データの取得時刻ｔｃと、該最後のセクタ内位置データに隣接するセクタ外位置データの取得時刻ｔｙとの按分点（一例として中点）に相当する時刻ｔｏｕｔを推定滞在期間の終了時刻とする。これにより、図３４に長方形で示した期間、即ち、時刻ｔｉｎから時刻ｔｏｕｔまでの期間が、当該移動機の推定滞在期間として算出される。
【０１４０】
そして、抽出部６８１は、算出された推定滞在期間と集計時間帯とが重複する移動機を抽出する。図３５には、拡張期間内に取得された移動機Ａ〜Ｅに関する一連の位置登録情報を移動機毎に、位置登録情報の取得時刻順に時系列に横軸方向に並べた図を示しており、この図３５では、セクタＩＤがある所定のセクタを示すセクタ内位置データを黒丸で、セクタＩＤがある所定のセクタ外を示すセクタ外位置データを白丸でそれぞれ示している。仮に、位置登録情報の取得時刻が集計時間帯（時刻ｔ０〜時刻ｔ１）内にあるセクタ内位置データにのみ着目すると、対応する移動機としては移動機Ｃ〜Ｅの３つが抽出される。
【０１４１】
ところが、図３５で長方形により示した推定滞在期間と集計時間帯（時刻ｔ０〜時刻ｔ１）とが重複する移動機としては、移動機Ａ〜Ｅの５つが挙げられるため、抽出部６８１は、移動機Ａ〜Ｅの５つを抽出する。これにより、黒丸で示すセクタ内位置データにのみ着目した場合に除外される移動機Ａ、Ｂについても、集計時間帯中のどこかである所定のセクタ内に滞在していたと推定されることとなる。この処理を、セクタ毎に順次行うことで、集計時間帯中のどこかで各セクタ内に滞在していた移動機をセクタ毎に抽出することができる。なお、上述した推定滞在期間を用いた方法は一例であり、他の方法を採用しても構わない。他の方法については後述する。
【０１４２】
入込端末数推計部６８２は、抽出部６８１によって抽出された移動機の数をセクタＩＤ毎にカウントし、得られたカウント数を集計時間帯内における各セクタ内の入込端末数として推計する。
【０１４３】
以上のように、第１取得部６０７Ｄにおいてセクタ毎の入込端末数を求める。分布算出部６０８は、第１取得部６０７Ｄで求められた入込端末数と、補正データとに基づいてセクタ内における移動機１００の分布（入込端末数の分布）を求める。また、集計部６１０は、分布算出部６０８によって求められた分割エリア毎の移動機数を用いて、集計エリア内の移動機数（入込端末数）を求める。この場合には、移動機１００の入込端末数の分布や集計エリア内の移動機数（入込端末数）を求めることができる。
【０１４４】
なお、前述した第１取得部６０７Ｄにおける推定滞在期間を用いた入込端末数の推計方法は一例であり、他の方法を採用しても構わない。別の例として、抽出部６８１は、同一移動機に関する位置登録情報において、位置登録情報が取得された時刻が集計時間帯内にあるセクタ内位置データのうち、１つの位置登録情報（一例として取得時刻が最も早い位置登録情報）を抽出してもよい。このようにして移動機毎に１つの位置登録情報を抽出する場合、上述した推定滞在期間を用いた方法よりも抽出の精度がやや劣るものの、移動機毎の推定滞在期間を算出する必要がないため、低い処理負荷で移動機毎に１つの位置登録情報を抽出することができる。この場合、入込端末数推計部６８２は、抽出された位置登録情報の数をカウントすることで入込端末数を推計すればよい。なお、移動機毎に１つの位置登録情報を抽出する場合、取得時刻が最も早い位置登録情報を抽出することは必須ではなく、他の１つの位置登録情報を抽出してもよい。例えば、取得時刻が最も遅い位置登録情報を抽出してもよいし、取得時刻が観測期間のちょうど真ん中の時点に最も近い位置登録情報を抽出してもよい。
【０１４５】
なお、上記各実施形態では、移動機１００の位置情報として、ＧＰＳ位置情報と位置登録情報とを用いるものとしたが、位置を把握できるものであれば、これ以外の情報であってもよい。
【０１４６】
また、上記各実施形態において、ＧＰＳ位置情報や位置登録情報に含まれるユーザ識別子（例えば電話番号）に対し、一方向性関数による不可逆符号への変換を含む非識別化処理を行ってもよい。この一方向性関数として、国内外の評価プロジェクトや評価機関により推奨されているハッシュ関数に基づく鍵付ハッシュ関数を用いることができる。
【０１４７】
また、上記各実施形態において、情報分析装置６００，６００Ａ，６００Ｂ，６００Ｃ，６００Ｄは、位置登録情報やＧＰＳ位置情報を交換機４００から取得するものとしたが、情報分析装置以外の構成は一例であって、位置登録情報やＧＰＳ位置情報を交換機４００以外から取得してもよい。例えば、情報分析装置は、上述したＲＮＣ３００や交換機４００等を備えた通信システムとは異なるシステムから位置登録情報やＧＰＳ位置情報を取得したり、位置登録情報やＧＰＳ位置情報が格納されたデータセンター等からこれらの情報を取得したりすることができる。
【符号の説明】
【０１４８】
１０，１０Ａ，１０Ｂ…移動通信システム、１００…移動機、６０１…ＧＰＳ位置情報記憶部、６０２，６０２Ｂ，６０２Ｃ…第２取得部、６０４，６０４Ｂ…補正データ算出部、６０６…位置登録情報記憶部、６０７，６０７Ｃ，６０７Ｄ…第１取得部、６０８…分布算出部、６１５Ａ…属性情報記憶部、６２２，６７２…前後位置情報取得部、６２３，６７３…特徴量計算部、６２５，６７５…移動機数推計部、６８１…抽出部、６８２…入込端末数推計部。
【図１】
2012054921.tif 000007
【図２】
2012054921.tif 000008
【図３】
2012054921.tif 000009
【図４】
2012054921.tif 000010
【図５】
2012054921.tif 000011
【図６】
2012054921.tif 000012
【図７】
2012054921.tif 000013
【図８】
2012054921.tif 000014
【図９】
2012054921.tif 000015
【図１０】
2012054921.tif 000016
【図１１】
2012054921.tif 000017
【図１２】
2012054921.tif 000018
【図１３】
2012054921.tif 000019
【図１４】
2012054921.tif 000020
【図１５】
2012054921.tif 000021
【図１６】
2012054921.tif 000022
【図１７】
2012054921.tif 000023
【図１８】
2012054921.tif 000024
【図１９】
2012054921.tif 000025
【図２０】
2012054921.tif 000026
【図２１】
2012054921.tif 000027
【図２２】
2012054921.tif 000028
【図２３】
2012054921.tif 000029
【図２４】
2012054921.tif 000030
【図２５】
2012054921.tif 000031
【図２６】
2012054921.tif 000032
【図２７】
2012054921.tif 000033
【図２８】
2012054921.tif 000034
【図２９】
2012054921.tif 000035
【図３０】
2012054921.tif 000036
【図３１】
2012054921.tif 000037
【図３２】
2012054921.tif 000038
【図３３】
2012054921.tif 000039
【図３４】
2012054921.tif 000040
【図３５】
2012054921.tif 000041