(57)【特許請求の範囲】
【請求項１】
携帯端末を所持したユーザの有意圏を推定することができる、広域無線通信網に接続されたユーザ有意圏推定装置であって、
携帯端末が配下に位置する基地局から、日時刻及び基地局識別子の多数のデータを含む接続履歴を収集し、携帯端末毎に区分する接続履歴収集手段と、
前記接続履歴を所定の時間窓（時間区間）に分割し、時間窓毎に各基地局識別子の出現数を計数する時間窓分割手段と、
前記時間窓毎に、複数の基地局識別子を要素とし、出現数を値とするベクトルを生成し、基地局識別子の出現数に依存するユーザの潜在状態（移動状態／滞留状態）を確率分布によって推定し、各時間窓の潜在状態を各時間窓が所属するクラスタとして、当該ベクトルをクラスタに凝集するクラスタリング手段と、
前記クラスタ毎に、当該クラスタに含まれるベクトルの各基地局識別子の出現数の偏りに基づいてユーザの有意圏を推定するユーザ有意圏推定手段と
を有することを特徴とするユーザ有意圏推定装置。
【請求項２】
前記クラスタ毎に、当該クラスタのエントロピー（クラスタ内の基地局分布のランダムさ）が低いクラスタほど「滞留状態」と判定する滞留判定手段を更に有し、
前記ユーザ有意圏推定手段は、滞留状態と判定されたクラスタに基づく地域圏を、当該ユーザの生活圏として推定する
ことを特徴とする請求項１に記載のユーザ有意圏推定装置。
【請求項３】
前記クラスタ毎に、当該クラスタに含まれるベクトルの分散が低いクラスタほど「滞留状態」と判定する滞留判定手段を更に有し、
前記ユーザ有意圏推定手段は、滞留状態と判定されたクラスタに基づく地域圏を、当該ユーザの生活圏として推定する
ことを特徴とする請求項１に記載のユーザ有意圏推定装置。
【請求項４】
前記クラスタリング手段は、時間窓を「文書」とし且つ基地局識別子を「単語」として、文書及び単語に対する潜在状態（トピック）を推定するＨＤＰ(Hierarchical Dirichlet Process)−ＬＤＡ(Latent Dirichlet Allocation) を用いる
ことを特徴とする請求項１から３のいずれか１項に記載のユーザ有意圏推定装置。
【請求項５】
前記接続履歴収集手段は、前記接続履歴を、カレンダ又はユーザ設定による平日及び休日の日属性に基づいて区分し、
前記ユーザ有意圏推定手段は、
前記滞留状態について、自宅のような「住所」、及び、職場・学校のような「居所」に区分し、
全日滞在率が高い第１のクラスタを前記住所と規定し、第１のクラスタ以外に平日滞在率が高い第２のクラスタを前記居所と規定するか、
全日滞在率が高い２つのクラスタを選択し、平日滞在率が高い一方のクラスタを前記居所と規定し、他方のクラスタを前記住所と規定するか、又は、
全日滞在率が高い２つのクラスタを選択し、休日滞在率が高い一方のクラスタを前記住所と規定し、他方のクラスタを前記居所と規定する
ことによって、滞留状態と判定されたクラスタの地域圏が、住所／居所を推定することを特徴とする請求項１から４のいずれか１項に記載のユーザ有意圏推定装置。
【請求項６】
前記広域無線通信網は、携帯電話網であり、
前記携帯端末は、携帯電話機であり、
前記ユーザ有意圏推定装置は、携帯電話網の通信事業者設備として設置されたものである
ことを特徴とする請求項１から５のいずれか１項に記載のユーザ有意圏推定装置。
【請求項７】
携帯端末を所持したユーザの有意圏を推定することができる、広域無線通信網に接続された通信事業者設備装置におけるユーザ有意圏推定方法であって、
携帯端末が配下に位置する基地局から、日時刻及び基地局識別子の多数のデータを含む接続履歴を収集し、携帯端末毎に区分する第１のステップと、
前記接続履歴を所定の時間窓（時間区間）に分割し、時間窓毎に各基地局識別子の出現数を計数する第２のステップと、
前記時間窓毎に、複数の基地局識別子を要素とし、出現数を値とするベクトルを生成し、基地局識別子の出現数に依存するユーザの潜在状態（移動状態／滞留状態）を確率分布によって推定し、各時間窓の潜在状態を各時間窓が所属するクラスタとして、当該ベクトルをクラスタに凝集する第３のステップと、
前記クラスタ毎に、当該クラスタに含まれるベクトルの各基地局識別子の出現数の偏りに基づいてユーザの有意圏を推定する第４のステップと
を有することを特徴とするユーザ有意圏推定方法。
【請求項８】
携帯端末を所持したユーザの有意圏を推定することができる、広域無線通信網に接続された通信事業者設備装置に搭載されたコンピュータにおけるユーザ有意圏推定プログラムであって、
携帯端末が配下に位置する基地局から、日時刻及び基地局識別子の多数のデータを含む接続履歴を収集し、携帯端末毎に区分する接続履歴収集手段と、
前記接続履歴を所定の時間窓（時間区間）に分割し、時間窓毎に各基地局識別子の出現数を計数する時間窓分割手段と、
前記時間窓毎に、複数の基地局識別子を要素とし、出現数を値とするベクトルを生成し、基地局識別子の出現数に依存するユーザの潜在状態（移動状態／滞留状態）を確率分布によって推定し、各時間窓の潜在状態を各時間窓が所属するクラスタとして、当該ベクトルをクラスタに凝集するクラスタリング手段と、
前記クラスタ毎に、当該クラスタに含まれるベクトルの各基地局識別子の出現数の偏りに基づいてユーザの有意圏を推定するユーザ有意圏推定手段と
してコンピュータを機能させることを特徴とするユーザ有意圏推定プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、携帯端末を所持したユーザの移動に伴う位置を推定する技術に関する。
【背景技術】
【０００２】
近年、携帯電話機に代表される携帯端末には、ＧＰＳ(Global Positioning System)のような測位機能が、一般的に搭載されてきている。そのため、ユーザは、携帯端末を用いて、現在位置を測位できると共に、その位置をネットワークを介してサーバへ送信することによって、位置に応じた様々なサービス情報を受信することができる。
【０００３】
従来、携帯端末のＧＰＳ機能によって取得された位置情報をサーバへ送信し、当該サーバが、そのユーザの行動履歴から行動範囲を算出し、その行動範囲を反映した情報を提供する技術がある（例えば特許文献１参照）。この技術によれば、携帯端末によって計測された位置情報同士の距離に基づいてクラスタリングされる。
【０００４】
このような技術によれば、携帯端末を所持したユーザの行動範囲を算出することができる一方で、その位置がどのような意味を持つのか（有意な位置であるか否か）についてまで認識することはできない。サーバが、携帯端末へその位置に応じたサービス情報を提供する場合であっても、その位置における「意味付け」（有意性）によっては、提供すべきサービス情報も異なる。例えば、滞在時間が長い場所であっても、自宅（住所）と職場（居所）とでは、提供すべきサービス情報の内容も全く異なる。通常、測位位置に対する「意味付け」は、ユーザによって指定されるか、又は、他の登録情報を用いる必要がある。
【０００５】
これに対し、携帯端末のＧＰＳ機能によって取得された位置情報に基づく行動履歴から、ユーザの有意位置を学習する技術がある（例えば非特許文献１参照）。この技術によれば、位置情報は、ＧＰＳ機能によって定期的に取得される必要がある。
【０００６】
また、携帯端末のＧＰＳ機能によって取得された位置情報に基づく行動履歴から、混合ガウス分布を用いて、ユーザの有意位置を学習する技術がある（例えば非特許文献２参照）。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２０１０−４９２９５号公報
【非特許文献】
【０００８】
【非特許文献１】遠山緑生、服部隆志、荻野達也、「携帯電話の測位機能を用いた有意位置の学習」、情報処理学会論文誌、vol.46 No.12、pp.2915-2924、2005
【非特許文献２】Petteri Nurmi、SouravBhattacharya、「Identifying Meaningful Places: The Non-parametric Way」、Pervasive2008、pp.111-127、2008
【発明の概要】
【発明が解決しようとする課題】
【０００９】
前述した従来技術によれば、ユーザの位置情報は、携帯端末のＧＰＳ機能によって取得されることを前提としている。しかしながら、ユーザによって所持された携帯端末について、ＧＰＳ機能及びそのアプリケーションを常に又は定期的に起動させることは、携帯端末の電池の消耗を早めだけでなく、携帯端末からのパケットの送出量を増やすという問題があった。そのために、ユーザの行動圏に基づくサービスを広く普及させることも難しい。
【００１０】
また、非特許文献１に記載の技術によれば、一定時間内での移動距離（速度）に基づいて滞在状態を判定するために、携帯端末における位置情報の計測時間間隔が一定である必要がある。しかしながら、位置情報の計測時間間隔が一定でない場合、直線的な移動であると仮定できる十分に短い時間間隔でない限り、実際の速度との誤差が大きくなる。即ち、滞在位置の判定の精度が低下する。
【００１１】
更に、非特許文献２に記載の技術によれば、空間的に疎な位置情報履歴（位置情報同士の地理的距離が比較的長い）を用いた場合、混合ガウス分布のパラメータ推定の性質によって、離れた位置情報同士を、同一のクラスタに含めてしまうという傾向がある。これによっても、滞在位置の判定精度が低下する。
【００１２】
このような従来技術に対して、通信事業者側としては、携帯端末によって取得された位置情報ではなく、その携帯端末が配下となる基地局の位置情報によって、携帯端末の有意位置を推定できることが好ましい。この場合、携帯端末が常に又は定期的にＧＰＳ機能を起動させる必要もない。
【００１３】
尚、特許文献１に記載の技術によれば、携帯端末が配下となる基地局の遷移履歴から、ユーザの行動履歴を算出することも想定できる。しかしながら、定期的な測位位置が得られない状況では、アルゴリズムを適用することが難しい。
【００１４】
そこで、本発明では、携帯端末の測位機能を起動させることなく、通信事業者設備によって取得可能な、空間的粒度が粗く且つ時間間隔が一定でない基地局位置情報を用いて、ユーザにとって有意な位置の集合であるところのユーザの有意圏を推定することができる装置、プログラム及び方法を提供することを目的とする。
【課題を解決するための手段】
【００１５】
本発明によれば、携帯端末を所持したユーザの有意圏を推定することができる、広域無線通信網に接続されたユーザ有意圏推定装置であって、
携帯端末が配下に位置する基地局から、日時刻及び基地局識別子の多数のデータを含む接続履歴を収集し、携帯端末毎に区分する接続履歴収集手段と、
接続履歴を所定の時間窓（時間区間）に分割し、時間窓毎に各基地局識別子の出現数を計数する時間窓分割手段と、
時間窓毎に、複数の基地局識別子を要素とし、出現数を値とするベクトルを生成し、基地局識別子の出現数に依存するユーザの潜在状態（移動状態／滞留状態）を確率分布によって推定し、各時間窓の潜在状態を各時間窓が所属するクラスタとして、当該ベクトルをクラスタに凝集するクラスタリング手段と、
クラスタ毎に、当該クラスタに含まれるベクトルの各基地局識別子の出現数の偏りに基づいてユーザの有意圏を推定するユーザ有意圏推定手段と
を有することを特徴とする。
【００１６】
本発明のユーザ有意圏推定装置における他の実施形態によれば、
クラスタ毎に、当該クラスタのエントロピー（クラスタ内の基地局分布のランダムさ）が低いクラスタほど「滞留状態」と判定する滞留判定手段を更に有し、
ユーザ有意圏推定手段は、滞留状態と判定されたクラスタに基づく地域圏を、当該ユーザの生活圏として推定することも好ましい。
【００１７】
本発明のユーザ有意圏推定装置における他の実施形態によれば、
クラスタ毎に、当該クラスタに含まれるベクトルの分散が低いクラスタほど「滞留状態」と判定する滞留判定手段を更に有し、
ユーザ有意圏推定手段は、滞留状態と判定されたクラスタに基づく地域圏を、当該ユーザの生活圏として推定することも好ましい。
【００１８】
本発明のユーザ有意圏推定装置における他の実施形態によれば、
クラスタリング手段は、時間窓を「文書」とし且つ基地局識別子を「単語」として、文書及び単語に対する潜在状態（トピック）を推定するＨＤＰ(Hierarchical Dirichlet Process)−ＬＤＡ(Latent Dirichlet Allocation) を用いる
ことも好ましい。
【００１９】
本発明のユーザ有意圏推定装置における他の実施形態によれば、
接続履歴収集手段は、接続履歴を、カレンダ又はユーザ設定による平日及び休日の日属性に基づいて区分し、
ユーザ有意圏推定手段は、
滞留状態について、自宅のような「住所」、及び、職場・学校のような「居所」に区分し、
全日滞在率が高い第１のクラスタを住所と規定し、第１のクラスタ以外に平日滞在率が高い第２のクラスタを居所と規定するか、
全日滞在率が高い２つのクラスタを選択し、平日滞在率が高い一方のクラスタを居所と規定し、他方のクラスタを住所と規定するか、又は、
全日滞在率が高い２つのクラスタを選択し、休日滞在率が高い一方のクラスタを住所と規定し、他方のクラスタを居所と規定する
ことによって、滞留状態と判定されたクラスタの地域圏が、住所／居所を推定することも好ましい。
【００２０】
本発明のユーザ有意圏推定装置における他の実施形態によれば、
広域無線通信網は、携帯電話網であり、
携帯端末は、携帯電話機であり、
ユーザ有意圏推定装置は、携帯電話網の通信事業者設備として設置されたものであることも好ましい。
【００２１】
本発明によれば、携帯端末を所持したユーザの有意圏を推定することができる、広域無線通信網に接続された通信事業者設備装置におけるユーザ有意圏推定方法であって、
携帯端末が配下に位置する基地局から、日時刻及び基地局識別子の多数のデータを含む接続履歴を収集し、携帯端末毎に区分する第１のステップと、
接続履歴を所定の時間窓（時間区間）に分割し、時間窓毎に各基地局識別子の出現数を計数する第２のステップと、
時間窓毎に、複数の基地局識別子を要素とし、出現数を値とするベクトルを生成し、基地局識別子の出現数に依存するユーザの潜在状態（移動状態／滞留状態）を確率分布によって推定し、各時間窓の潜在状態を各時間窓が所属するクラスタとして、当該ベクトルをクラスタに凝集する第３のステップと、
クラスタ毎に、当該クラスタに含まれるベクトルの各基地局識別子の出現数に基づいてユーザの有意圏を推定する第４のステップと
を有することを特徴とする。
【００２２】
本発明によれば、携帯端末を所持したユーザの有意圏を推定することができる、広域無線通信網に接続された通信事業者設備装置に搭載されたコンピュータにおけるユーザ有意圏推定プログラムであって、
携帯端末が配下に位置する基地局から、日時刻及び基地局識別子の多数のデータを含む接続履歴を収集し、携帯端末毎に区分する接続履歴収集手段と、
接続履歴を所定の時間窓（時間区間）に分割し、時間窓毎に各基地局識別子の出現数を計数する時間窓分割手段と、
時間窓毎に、複数の基地局識別子を要素とし、出現数を値とするベクトルを生成し、基地局識別子の出現数に依存するユーザの潜在状態（移動状態／滞留状態）を確率分布によって推定し、各時間窓の潜在状態を各時間窓が所属するクラスタとして、当該ベクトルをクラスタに凝集するクラスタリング手段と、
クラスタ毎に、当該クラスタに含まれるベクトルの各基地局識別子の出現数に基づいてユーザの有意圏を推定するユーザ有意圏推定手段と
してコンピュータを機能させることを特徴とするユーザ有意圏推定プログラム。
【発明の効果】
【００２３】
本発明の装置、プログラム及び方法によれば、携帯端末の測位機能を起動させることなく、通信事業者設備によって取得可能な、空間的粒度が粗く且つ時間間隔が一定でない基地局位置情報を用いて、ユーザの有意位置を推定することができる。
【図面の簡単な説明】
【００２４】
【図１】携帯端末の移動を表す空間的な外観図である。
【図２】本発明におけるユーザ有意圏推定装置の機能構成図である。
【図３】本発明のフローチャートに基づく説明図である。
【図４】ＬＤＡの確率ネットワーク図である。
【図５】ＨＤＰ−ＬＤＡの確率ネットワーク図である。
【発明を実施するための形態】
【００２５】
以下、本発明の実施の形態について、図面を用いて詳細に説明する。
【００２６】
図１は、携帯端末の移動を表す空間的な外観図である。
【００２７】
ユーザに所持された携帯端末（例えば携帯電話機）は、いずれの位置にあっても、常に、基地局の配下にあってその基地局からの電波を受信し続けている。図１によれば、ユーザは、自宅の住所が「埼玉県ふじみ野市」にあり、職場の居所が「東京都港区飯田橋」にある。そのユーザは、自宅と職場との間を、池袋を経由して通勤している。また、そのユーザは、訪問先として「東京都千代田区大手町」へも頻繁に移動している。
【００２８】
多数の基地局を統合する通信事業者設備では、携帯端末毎に、空間的粒度が粗く、且つ、時間間隔が一定でない基地局位置情報を常に収集することができる。「空間的粒度が粗く」とは、位置情報同士の地理的な距離が比較的長いことを意味する。また、「時間間隔が一定でない」とは、位置情報の取得時間間隔が比較的ばらついていることを意味する。
【００２９】
本発明について、「住所」とは、ユーザ各人の生活の本拠をいい、例えば生活の中心となっている自宅の所在地を意味する。また、「居所」とは、継続して居るものの、生活の本拠というほどその場所との結びつきが強くない場所をいい、例えば職場又は学校の所在地をいう。
【００３０】
図２は、本発明におけるユーザ有意圏推定装置の機能構成図である。
図３は、本発明のフローチャートに基づく説明図である。
【００３１】
広域無線通信網（携帯電話網）に接続された基地局３は、無線フレーム信号を交換することによって、その配下に位置する携帯端末（携帯電話機）２を認識し、その日時刻を取得する。そして、基地局３は、携帯電話機２の「端末識別子」（アドレス、電話番号、識別番号等）と、「基地局識別子」と、「日時刻」とを含むデータを、ユーザ有意圏推定装置１へ送信する。
データ（端末識別子、基地局識別子、日時刻）
【００３２】
ユーザ有意圏推定装置１は、広域無線通信網（携帯電話網）に接続されており、携帯端末（携帯電話機）を所持したユーザの有意圏を推定することができる。図２によれば、ユーザ有意圏推定装置１は、通信事業者網に接続する通信インタフェース部１０と、接続履歴収集部１１と、時間窓分割部１２と、クラスタリング部１３と、滞留判定部１４と、ユーザ有意圏推定部１５と、ユーザ有意圏登録部１６と、アプリケーション処理部１７とを有する。アプリケーション処理部１７は、本発明によって推定されたユーザ毎の有意圏に基づいて、様々なサービスを実行する。通信インタフェース部を除くこれら機能構成部は、装置に搭載されたコンピュータを機能させるプログラムを実行することによって実現される。
【００３３】
［接続履歴収集部］
接続履歴収集部１１（図３のＳ１１）は、携帯端末２が配下に位置する基地局３から、日時刻及び基地局識別子の多数のデータを含む接続履歴を収集し、携帯端末毎に区分する。
【００３４】
表１は、１つの携帯端末について収集された接続履歴を表す。
【表１】
0005523274.tif 000002
【００３５】
［時間窓分割部］
時間窓分割部１２（図３のＳ１２）は、接続履歴を所定の時間窓（時間区間）に分割する。時間窓は、時間幅Ｔ及びシフト幅Ｓによって決定される。シフト幅Ｓとは、開始時刻をＳだけ遅らせたものである。即ち、Ｔ＞Ｓの場合、時間窓は、Ｔ−Ｓだけ重畳することとなる。
【００３６】
表２は、Ｔ＝６０分及びＳ＝１５分とした場合における、各時間窓の開始時刻及び終了時刻を表す。
【表２】
0005523274.tif 000003
【００３７】
次に、時間窓毎に、各基地局識別子の出現数を計数する。
18:00:00〜18:59:59の時間窓について、基地局ＩＤの個数を計数すると、「基地１」は４個、「基地２」は３個となる。
18:15:00〜19:14:59の時間窓について、基地局ＩＤの個数を計数すると、「基地１」は２個、「基地２」は３個となる。
以下同様に計数される。
【００３８】
表３は、時間窓毎に、各基地局ＩＤの出現数を表す。
【表３】
0005523274.tif 000004
【００３９】
［クラスタリング部］
クラスタリング部１３（図３のＳ１３）は、時間窓毎に、複数の基地局識別子を要素とし、出現数を値とするベクトルを生成し、当該ベクトルをクラスタに凝集する。
【００４０】
クラスタリング部１３は、時間窓を「文書」とし且つ基地局識別子を「単語」として、文書及び単語に対する潜在状態（トピック）を推定するＨＤＰ(Hierarchical Dirichlet Process)−ＬＤＡ(Latent Dirichlet Allocation)を用いる。これによって、各時間窓について推定された潜在状態を、各時間窓が所属するクラスタとする。
【００４１】
ここで、各時間窓における基地局ＩＤの出現数は、ユーザの潜在状態に依存する確率分布に従って生成されたと想定できる。ユーザの潜在状態は、時間窓毎に、異なる多項分布に従って生成される。
【００４２】
例えば、携帯端末を所持するユーザが、滞留点Ａから滞留点Ｂへ移動した場合、潜在状態として以下の３つが想定できる。
・「滞留点Ａ」
・「滞留点Ｂ」
・「移動中」（適宜さらに分割される場合もあるが、仮にひとつのまとまりとする）
【００４３】
潜在状態が「滞留点」である場合には、携帯端末は、その付近の少数個の基地局との間で通信イベントが発生することが想定される。一方で、潜在状態が「移動中」の場合には、移動経路上の多くの基地局との間で通信イベントが発生することが想定される。
【００４４】
表４は、時間窓毎における各基地局の出現数に基づいて、潜在状態を表す。
【表４】
0005523274.tif 000005
【００４５】
＜ＬＤＡ＞
図４は、ＬＤＡの確率ネットワーク図である。
ｗij：時間窓iで、j番目に観測された基地局ＩＤ
ｚij：潜在状態
θi：時間窓毎の潜在状態分布パラメータ
α、β：超パラメータ
ｚij、θiの事後確率分布と、α、βの最適値（第二種最尤推定値）を求めることによって、携帯端末毎に基地局との接続履歴から、ユーザの潜在状態を推定し、時間窓毎に潜在状態を分類する（各潜在状態への所属確率を求める）ことができる。βについては、更に事前分布を設定することもある。
【００４６】
＜ＨＤＰ−ＬＤＡ＞
図５は、ＨＤＰ−ＬＤＡの確率ネットワーク図である。
ＬＤＡによれば、潜在状態の数を予め与える必要がある。これに対し、ＨＤＰ−ＬＤＡによれば、データの複雑さに合わせて必要な数の潜在状態数を自動的に決めることができる。そのために、ＨＤＰ−ＬＤＡによれば、Ｇ0、θi、ｚijの分布と、α及びβの最適値を推定する。これによって、θiの次元、即ち、潜在クラス数は、Dirichlet過程のパラメータαに従ってデータから決まるため、予め決定しておく必要もない。ＬＤＡ又はＨＤＰ−ＬＤＡのいずれについても、実際の推定計算には、変分ベイズ法又はマルコフ連鎖モンテカルロ法が用いられる。
【００４７】
表５は、ＨＤＰ−ＬＤＡのクラスタリングの計算結果を表す。これは、時間窓毎の各クラスタへの所属確率を表す。
【表５】
0005523274.tif 000006
【００４８】
表６は、ＨＤＰ−ＬＤＡを用いたクラスタリングの過程で計算される超パラメータβの結果を表す。βは、各潜在状態について、各基地局と通信する確率を表すパラメータの集合である。
【表６】
0005523274.tif 000007
【００４９】
［滞留判定部］
滞留判定部１４（図３のＳ１４）は、クラスタ毎に、「滞留状態」か否かを判定する。ここで、「滞留状態」を判定するために、以下の２つの方法がある。
＜１＞エントロピー（クラスタ内の基地局分布のランダムさ）が低いクラスタほど「滞留状態」と判定する。逆に、エントロピーが高いクラスタほど「移動状態」と判定する。
＜２＞当該クラスタに含まれるベクトルの分散が低いクラスタほど「滞留状態」と判定する。逆に、当該クラスタに含まれるベクトルの分散が高いクラスタほど「移動状態」と判定する。
【００５０】
＜エントロピーを用いる方法＞
HDP-LDAを用いたクラスタリングの過程で計算される超パラメータβ（各潜在状態において各基地局と通信をする確率を表すパラメータの集合）の結果を用いて、各クラスタ内の基地局分布のランダムさの程度を表すエントロピーを計算し、その大きさをもとに滞留判定を行う。
【００５１】
「エントロピー」とは、情報理論の概念であって、事象の生起のランダムさを表す尺度である。どのイベントが起こるかはっきりしている（ランダムさが小さい）状態（例えば滞留状態）には何らかの意味がある可能性が高く、どのイベントが起こるかはっきりしない（ランダムさが大きい）状態にはあまり意味がない。
【００５２】
（１）クラスタiについて、各基地局jと通信する確率ｐijを用いてクラスタiのエントロピーentropyiを計算する。エントロピーは、クラスタの分散の程度を表す。
entropyi＝−Σjｐijlogｐij
例えば、クラスタ１の場合、各基地局と通信する確率は、表６によれば以下のようになる。
{0.49, 0.49, 0.01, 0.01, 0.0, 0.0, 0.0}
従って、エントロピーは、以下のように表される。
entropyi＝-0.49log(0.49)-0.49log(0.49)-0.01log(0.01)-0.01log(0.01)
≒0.34
（２）次に、クラスタを、エントロピーの昇順でソートする。
（３）そして、エントロピーが閾値th以上のクラスタを「移動状態クラスタ」、th以下のクラスタを「滞留状態クラスタ」と判定する。
【００５３】
表７は、エントロピーを用いた場合における、クラスタ毎の滞留状態の判定結果を表す。ここでは、例えばth=0.5とした場合であって、クラスタ１及び２は「滞留状態」と判定され、「クラスタ３」は「移動状態」と判定される。
【表７】
0005523274.tif 000008
【００５４】
＜分散を用いる方法＞
表８は、各基地局の位置情報を表す。このように、予め固定されている基地局の位置情報を用いて、クラスタの位置情報の分散共分散行列を計算し、その最大固有値の大きさに基づいて滞留状態を判定する。
【表８】
0005523274.tif 000009
【００５５】
（１）エリアに合わせて、緯度経度から分散共分散行列計算のための基地局位置情報を求める。表９は、補正後の基地局位置を表す。ここでは、前述の緯度経度の平均位置である緯度経度（35.4, 135.8）付近の緯度1度≒111.3km、経度1度≒90.8kmに基づいて、以下の式により、基地局jの緯度latj，経度lonjから基地局jの補正位置ｘj，ｙjを得る。楕円体を用いて、更に正確に補正することも好ましい。
【数１】
0005523274.tif 000010
【表９】
0005523274.tif 000011
【００５６】
（２）次に、時間窓毎の各クラスタへの所属確率（表５参照）と、時間窓毎の基地局の個数（表６参照）とに基づいて、クラスタiで基地局jとの通信回数の期待値ｎijを計算する。
例えば、ｎ11＝0.99×4.0＋0.98×2.0＝5.92となる。
【００５７】
表１０は、クラスタ毎の各基地局との通信回数の期待値を表す。
【表１０】
0005523274.tif 000012
【００５８】
（３）以下の式によって、クラスタiの分散共分散行列Ｓiを計算する。
【数２】
0005523274.tif 000013
例えば、以下のようになる。
【数３】
0005523274.tif 000014
【００５９】
（４）次に、クラスタiの分散共分散行列の最大固有値（第１主成分得点の分散）λiを計算する。例えば、λ1＝178.83となる。
【００６０】
（５）次に、クラスタを、分散共分散行列の最大固有値の昇順でソートする。
【００６１】
（６）そして、分散共分散行列の最大固有値が閾値th以上のクラスタを「移動状態クラスタ」、th以下のクラスタを「滞留状態クラスタ」と判定する。
例えばth=2000とした場合、クラスタ１及び２「滞留状態クラスタ」と判定され、クラスタ３は「移動状態クラスタ」と判定される。
【００６２】
表１１は、分散を用いた場合における滞留状態の判定結果を表す。
【表１１】
0005523274.tif 000015
【００６３】
［ユーザ有意圏推定部］
ユーザ有意圏推定部１５（図３のＳ１５）は、クラスタ毎に、当該クラスタに含まれるベクトルの各基地局識別子の出現数に基づいてユーザの有意圏を推定する。有意圏は、滞留状態と判定されたクラスタに基づく地域圏であるとする。ここで、ユーザにとっての有意圏とは、「住所」又は「居所」とする。また、携帯端末から基地局に対する接続履歴は、カレンダ又はユーザ設定による平日及び休日の日属性に基づいて変換する。
【００６４】
表１２は、ユーザの有意圏の推定を表す。
【表１２】
0005523274.tif 000016
【００６５】
ここで、平日／休日の定義は、カレンダ情報に基づいて土曜日及び日曜日に設定されたものであってもよいし、ユーザに利用態様に応じて設定されたものであってもよい。携帯端末が、携帯電話機である場合、アプリケーションとして搭載されたカレンダ／スケジューラの機能を用いることもできる。
【００６６】
次に、クラスタ毎に、日属性に基づいて滞在率を算出する。滞在率は、例えば以下のように算出される。
Ｄ：総出現回数
Ｄw：平日出現回数
Ｄh：休日出現回数
ｎｄ(c)：クラスタcにおける総出現回数
ｎｄw(c)：クラスタcにおける平日出現回数
ｎｄh(c)：クラスタcにおける休日出現回数
Ｒ(c)：クラスタcにおける全日滞在率
Ｒ(c)＝ｎｄ(c)／Ｄ
Ｒw(c)：クラスタcにおける平日滞在率
Ｒw(c)＝ｎｄw(c)／Ｄw
Ｒh(c)：クラスタcにおける休日滞在率
Ｒh(c)＝ｎｄh(c)／Ｄh
【００６７】
ユーザ有意圏推定部１５は、滞在状態にあると判定されたクラスタについて、自宅のような「住所」、及び、職場・学校のような「居所」に区分する。そして、以下のような３つのいずれかの方法で、滞留状態と判定されたクラスタの地域圏に対して、住所／居所を推定することができる。
【００６８】
（１）平日休日を問わず、全日滞在率が高いクラスタＡを住所と規定し、クラスタＡ以外に平日滞在率が高いクラスタＢを居所と規定する。
【００６９】
（２）全日滞在率が高い２つのクラスタＡ及びクラスタＢを選択し、平日滞在率が高い一方のクラスタを「居所」と規定し、他方のクラスタを「住所」と規定する。
Ｒw(A)＞Ｒw(B)の場合：クラスタＡ「住所」、クラスタＢ「居所」
Ｒw(A)＜Ｒw(B)の場合：クラスタＡ「居所」、クラスタＢ「住所」
Ｒw(A)＝Ｒw(B)の場合：休日滞在率が高い方のクラスタを「住所」
【００７０】
（３）全日滞在率が高い２つのクラスタＡ及びクラスタＢを選択し、休日滞在率が高い一方のクラスタを「住所」と規定し、他方のクラスタを「居所」と規定する。
Ｒh(A)＞Ｒh(B)の場合：クラスタＡ「住所」、クラスタＢ「居所」
Ｒh(A)＜Ｒh(B)の場合：クラスタＡ「居所」、クラスタＢ「住所」
Ｒh(A)＝Ｒh(B)の場合：平日滞在率が高い方のクラスタを「居所」
【００７１】
例えば、表１２によれば、クラスタ２が「住所」として規定され、クラスタ１が「居所」として規定されている。
【００７２】
［ユーザ有意圏登録部］
ユーザ有意圏登録部１６は、クラスタＩＤ、基地局ＩＤ、所属率、緯度（任意）、経度（任意）、及び、有意圏ラベル（任意）から構成されたデータを登録する。ここで、所属率は、クラスタ毎に各基地局と通信をする確率であり、クラスタリング部１３によって算出されたβの値を格納する。また、クラスタについては現回数の降順にソートし、基地局については所属率の降順にソートする。
【００７３】
表１３は、有意圏登録データを表す。ここでは、所属率=0.05の閾値が設定されており、その閾値以上の基地局ＩＤが登録されている。また、所属率の代わりに、クラスタ毎に各基地局と通信する回数の期待値（表１０参照）に対して閾値を設けてもよい。
【表１３】
0005523274.tif 000017
【００７４】
尚、ユーザ有意圏登録部が、滞在状態（住所／居所）と判定されたクラスタについて、地域名称に対応付けて記憶することも好ましい。例えば、そのクラスタに属する基地局周辺の駅名を地域名称として記憶するものであってもよい。これら地域名称は、アプリケーション処理部１１３によって様々に利用することができる。
【００７５】
例えば、図１を参照すると、住所位置に最も近い「ふじみ野」駅を、住所位置の地域名称と規定される。また、居所位置に最も近い「飯田橋」駅を、居所位置の地域名称と規定される。従って、居所位置から離れた「水道橋」駅及び「神保町」駅は、居所位置の地域名称として選択されない。クラスタの重心位置が正しい住所位置からずれている場合であっても、地域名称は正しい名称をつけることができる。
【００７６】
以上、詳細に説明したように、本発明の装置、プログラム及び方法によれば、携帯端末の測位機能を起動させることなく、通信事業者設備によって取得可能であって、空間的粒度が粗く、且つ、時間間隔が一定でない基地局位置情報を用いて、ユーザの有意位置を推定することができる。
【００７７】
本発明によれば、携帯端末がＧＰＳのような測位機能を起動する必要がないので、携帯端末の電池の消耗を考慮する必要がない。即ち、ユーザの有意圏を推定するための全ての情報は、通信事業者側のみで取得できるものである。ユーザの有意圏を推定することによって提供できるサービスとしては、例えば、ユーザ毎に生活場所に応じたクーポン情報等を配信するパーソナライズド情報提供サービスなどがある。
【００７８】
前述した本発明の種々の実施形態について、本発明の技術思想及び見地の範囲の種々の変更、修正及び省略は、当業者によれば容易に行うことができる。前述の説明はあくまで例であって、何ら制約しようとするものではない。本発明は、特許請求の範囲及びその均等物として限定するものにのみ制約される。
【符号の説明】
【００７９】
１  ユーザ有意圏推定装置
１０  通信インタフェース部
１１  接続履歴収集部
１２  時間窓分割部
１３  クラスタリング部
１４  滞留判定部
１５  ユーザ有意圏推定部
１６  ユーザ有意圏登録部
１７  アプリケーション処理部
２  携帯端末、携帯電話機
３  基地局
【図１】
0005523274.tif 000018
【図２】
0005523274.tif 000019
【図３】
0005523274.tif 000020
【図４】
0005523274.tif 000021
【図５】
0005523274.tif 000022