(57)【特許請求の範囲】
【請求項１】
ユーザの移動に伴って所持される携帯端末であって、ＧＰＳ(Global Positioning System)の測位電波を受信する測位部と、前記測位部を所定時間周期で起動し且つ測位位置を取得する位置取得手段と、時刻に応じた前記測位位置を蓄積する位置蓄積手段とを有する携帯端末であって、
時系列の測位位置の移動に応じて、移動中位置以外の滞留位置のみを抽出する滞留位置抽出手段と、
複数の滞留位置を、空間的なクラスタに区分するクラスタリング手段と、
ＧＰＳの測位容易性の不安定化を吸収するべく、前記クラスタ毎に、前記滞留位置が１回でも存在する滞在日数を用いて、日属性に基づく滞在率を算出する滞在率算出手段と、
日属性に基づく滞在率が高いクラスタについて、当該クラスタの滞在特性を推定する滞在特性推定手段と
を有することを特徴とする携帯端末。
【請求項２】
前記滞在特性推定手段について、
前記日属性は、カレンダ又はユーザ設定による平日及び休日に基づいており、
前記滞在特性は、自宅のような「住所」、及び、職場・学校のような「居所」に基づいており、
全日滞在率が高い第１のクラスタを前記住所と規定し、第１のクラスタ以外に平日滞在率が高い第２のクラスタを前記居所と規定するか、
全日滞在率が高い２つのクラスタを選択し、平日滞在率が高い一方のクラスタを前記居所と規定し、他方のクラスタを前記住所と規定するか、又は、
全日滞在率が高い２つのクラスタを選択し、休日滞在率が高い一方のクラスタを前記住所と規定し、他方のクラスタを前記居所と規定する
ことを特徴とする請求項１に記載の携帯端末。
【請求項３】
前記滞在特性推定手段によって抽出された第１のクラスタ及び第２のクラスタについて、多数の滞留位置から重心位置を算出し、第１のクラスタの重心位置を「住所」位置として規定し、第２のクラスタの重心位置を「居所」位置として規定する重心位置決定手段を更に有する特徴とする請求項２に記載の携帯端末。
【請求項４】
地域名称に対応付けて位置を記憶した地域名称記憶手段と、
前記地域名称記憶手段を用いて、前記住所位置及び居所位置に最も近い「地域名称」を検索し、当該「地域名称」を、前記住所位置の住所及び前記居所位置の居所と規定する名前解決手段と
を更に有することを特徴とする請求項３に記載の携帯端末。
【請求項５】
前記滞留位置抽出手段は、
現時刻tから直近過去所定時間Ｔ以内に取得された複数の測位位置ｐt〜ｐt-nを平滑化し、現時刻tにおける平滑化位置ｐａtを算出し、
現時刻tの平滑化位置ｐａtに対する直近過去の平滑化位置ｐａt-1からの変位距離を算出し、
前記変位距離に対する変位時間から移動速度を算出し、
前記移動速度が所定閾値以下となる平滑化位置ｐａを「滞留位置」として抽出する
ことを特徴とする請求項１から４のいずれか１項に記載の携帯端末。
【請求項６】
前記位置は、緯度・経度のベクトルであり、
前記クラスタリング手段は、ＤＢＳＣＡＮ(Density-Based Spatial Clustering)を用いることを特徴とする請求項１から５のいずれか１項に記載の携帯端末。
【請求項７】
ユーザの移動に伴って所持される携帯端末と通信可能なサーバであり、
前記携帯端末は、ＧＰＳの測位電波を受信する測位部と、前記測位部を所定時間周期で起動し且つ測位位置を取得する位置取得手段と、時刻に応じた前記測位位置を蓄積する位置蓄積手段と、前記位置蓄積手段に蓄積された多数の測位位置を送信する位置送信手段とを有し、当該携帯端末から多数の測位位置を受信する滞在特性推定サーバであって、
時系列の測位位置の移動に応じて、移動中位置以外の滞留位置のみを抽出する滞留位置抽出手段と、
複数の滞留位置を、空間的なクラスタに区分するクラスタリング手段と、
ＧＰＳの測位容易性の不安定化を吸収するべく、前記クラスタ毎に、前記滞留位置が１回でも存在する滞在日数を用いて、日属性に基づく滞在率を算出する滞在率算出手段と、
日属性に基づく滞在率が高いクラスタについて、当該クラスタの滞在特性を推定する滞在特性推定手段と
を有することを特徴とする滞在特性推定サーバ。
【請求項８】
ＧＰＳの測位電波を受信する測位部を有し、ユーザの移動に伴って所持される携帯端末に搭載されたコンピュータを、前記測位部を所定時間周期で起動し且つ測位位置を取得する位置取得手段と、時刻に応じた前記測位位置を蓄積する位置蓄積手段として機能させるプログラムであって、
時系列の測位位置の移動に応じて、移動中位置以外の滞留位置のみを抽出する滞留位置抽出手段と、
複数の滞留位置を、空間的なクラスタに区分するクラスタリング手段と、
ＧＰＳの測位容易性の不安定化を吸収するべく、前記クラスタ毎に、前記滞留位置が１回でも存在する滞在日数を用いて、日属性に基づく滞在率を算出する滞在率算出手段と、
日属性に基づく滞在率が高いクラスタについて、当該クラスタの滞在特性を推定する滞在特性推定手段と
してコンピュータを機能させることを特徴とする携帯端末用のプログラム。
【請求項９】
ＧＰＳの測位電波を受信する測位部を有し、ユーザの移動に伴って所持される携帯端末であって、前記測位部を所定時間周期で起動し、測位位置を取得し、時刻に応じた前記測位位置を蓄積する携帯端末について、ユーザの住所及び居所を推定する方法であって、
時系列の測位位置の移動に応じて、移動中位置以外の滞留位置のみを抽出する第１のステップと、
複数の滞留位置を、空間的なクラスタに区分する第２のステップと、
ＧＰＳの測位容易性の不安定化を吸収するべく、前記クラスタ毎に、前記滞留位置が１回でも存在する滞在日数を用いて、日属性に基づく滞在率を算出する第３のステップと、
日属性に基づく滞在率が高いクラスタについて、当該クラスタの滞在特性を推定する第４のステップと
を有することを特徴とする携帯端末の住所居所推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、携帯端末に搭載された測位機能を用いて、ユーザの移動に伴う位置を取得する技術に関する。
【背景技術】
【０００２】
近年、携帯電話機に代表される携帯端末には、ＧＰＳ(Global Positioning System)のような測位機能が、一般的に搭載されてきている。そのため、ユーザは、携帯端末を用いて、現在位置を測位できると共に、その位置をネットワークを介してサーバへ送信することによって、位置に応じた様々なサービス情報を受信することができる。
【０００３】
従来、携帯端末の測位位置をサーバへ送信し、当該サーバが、その携帯端末を所持するユーザが頻繁に滞在するエリアを検出する技術がある（例えば特許文献１参照）。この技術によれば、エリア毎の滞在頻度を検出するために、地域を事前にエリア単位に区切る必要がある。例えば、東京都のような地域を、一定のメッシュに区切る必要がある。この技術は、エリアの区分の仕方によって、異なる目的・効果に用いられる。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００３−４４５１２号公報
【非特許文献】
【０００５】
【非特許文献１】遠山緑生、服部隆志、萩野達也、「携帯電話の測位機能を用いた有意位置の学習」、情報処理学会論文誌46(12)、pp.2915-2924、2005年
【非特許文献２】Ester M., et.al: “A Density-Based Algorithmfor Discovering Clusters in Large Spatial Databases with Noise”, Proc. 2nd int.Conf. on Knowledge Discovery and Data Mining (KDD‘96), Portland, Oregon, 1996,AAAI Press, 1996.
【非特許文献３】「駅データ」、[online]、［平成２２年２月１４日検索］、インターネット＜URL:http://www.ekidata.jp/index.html＞
【発明の概要】
【発明が解決しようとする課題】
【０００６】
携帯端末の現在位置の緯度経度情報は、ＧＰＳによって取得できる一方で、その現在位置がどのような意味を持つのか（有意な位置であるか否か）についてまで取得することはできない。サーバが、位置に応じたサービス情報を携帯端末へ提供する場合であっても、その位置における「意味付け」（有意性）によっては、提供すべきサービス情報も異なる。例えば、同程度の頻度で訪れる場所であっても、お気に入りの店と、取引先のオフィスとでは、提供すべきサービス情報の内容も全く異なる。有意な位置は、個人によって様々であるが、その中でも、自宅のような「住所」と、職場・学校のような「居所」とは、最も有意な位置と考えられる。特許文献１に記載された技術を用いることによって、「滞在頻度が多い場所」を抽出できたとしても、その有意性についてまで特定することできない。通常、測位位置に対する「意味付け」は、ユーザによって指定されるか、又は、他の登録情報を用いる必要がある。
【０００７】
そこで、本発明は、携帯端末に搭載された測位機能のみを用いて、ユーザの移動に伴う住所／居所を自動的に推定する携帯端末、サーバ、プログラム及び方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
本発明によれば、ユーザの移動に伴って所持される携帯端末であって、ＧＰＳ(Global Positioning System)の測位電波を受信する測位部と、測位部を所定時間周期で起動し且つ測位位置を取得する位置取得手段と、時刻に応じた測位位置を蓄積する位置蓄積手段とを有する携帯端末であって、
時系列の測位位置の移動に応じて、移動中位置以外の滞留位置のみを抽出する滞留位置抽出手段と、
複数の滞留位置を、空間的なクラスタに区分するクラスタリング手段と、
ＧＰＳの測位容易性の不安定化を吸収するべく、クラスタ毎に、滞留位置が１回でも存在する滞在日数を、日属性に基づいて計数する滞在率算出手段と、
滞在日数が多いクラスタから順に、日属性に基づいて、当該クラスタの滞在特性を推定する滞在特性推定手段と
を有することを特徴とする。
【０００９】
本発明の携帯端末における他の実施形態によれば、滞在特性推定手段について、
日属性は、カレンダ又はユーザ設定による平日及び休日に基づいており、
滞在特性は、自宅のような「住所」、及び、職場・学校のような「居所」に基づいており、
全日滞在率が高い第１のクラスタを住所と規定し、第１のクラスタ以外に平日滞在率が高い第２のクラスタを居所と規定するか、
全日滞在率が高い２つのクラスタを選択し、平日滞在率が高い一方のクラスタを居所と規定し、他方のクラスタを住所と規定するか、又は、
全日滞在率が高い２つのクラスタを選択し、休日滞在率が高い一方のクラスタを住所と規定し、他方のクラスタを居所と規定することも好ましい。
【００１０】
本発明の携帯端末における他の実施形態によれば、滞在特性推定手段によって抽出された第１のクラスタ及び第２のクラスタについて、多数の滞留位置から重心位置を算出し、第１のクラスタの重心位置を「住所」位置として規定し、第２のクラスタの重心位置を「居所」位置として規定する重心位置決定手段を更に有することも好ましい。
【００１１】
本発明の携帯端末における他の実施形態によれば、
地域名称に対応付けて位置を記憶した地域名称記憶手段と、
地域名称記憶手段を用いて、住所位置及び居所位置に最も近い「地域名称」を検索し、当該「地域名称」を、住所位置の住所及び居所位置の居所と規定する名前解決手段と
を更に有することも好ましい。
【００１２】
本発明の携帯端末における他の実施形態によれば、滞留位置抽出手段は、
現時刻tから直近過去所定時間Ｔ以内に取得された複数の測位位置ｐt〜ｐt-nを平滑化し、現時刻tにおける平滑化位置ｐａtを算出し、
現時刻tの平滑化位置ｐａtに対する直近過去の平滑化位置ｐａt-1からの変位距離を算出し、
変位距離に対する変位時間から移動速度を算出し、
移動速度が所定閾値以下となる平滑化位置ｐａを「滞留位置」として抽出することも好ましい。
【００１３】
本発明の携帯端末における他の実施形態によれば、位置は、緯度・経度のベクトルであり、クラスタリング手段は、ＤＢＳＣＡＮ(Density-Based Spatial Clustering)を用いることも好ましい。
【００１４】
本発明によれば、ユーザの移動に伴って所持される携帯端末と通信可能なサーバであり、
携帯端末は、ＧＰＳの測位電波を受信する測位部と、測位部を所定時間周期で起動し且つ測位位置を取得する位置取得手段と、時刻に応じた測位位置を蓄積する位置蓄積手段と、位置蓄積手段に蓄積された多数の測位位置を送信する位置送信手段とを有し、当該携帯端末から多数の測位位置を受信する滞在特性推定サーバであって、
時系列の測位位置の移動に応じて、移動中位置以外の滞留位置のみを抽出する滞留位置抽出手段と、
複数の滞留位置を、空間的なクラスタに区分するクラスタリング手段と、
ＧＰＳの測位容易性の不安定化を吸収するべく、クラスタ毎に、滞留位置が１回でも存在する滞在日数を、日属性に基づいて計数する滞在率算出手段と、
滞在日数が多いクラスタから順に、日属性に基づいて、当該クラスタの滞在特性を推定する滞在特性推定手段と
を有することを特徴とする。
【００１５】
本発明によれば、ＧＰＳの測位電波を受信する測位部を有し、ユーザの移動に伴って所持される携帯端末に搭載されたコンピュータを、測位部を所定時間周期で起動し且つ測位位置を取得する位置取得手段と、時刻に応じた測位位置を蓄積する位置蓄積手段として機能させるプログラムであって、
時系列の測位位置の移動に応じて、移動中位置以外の滞留位置のみを抽出する滞留位置抽出手段と、
複数の滞留位置を、空間的なクラスタに区分するクラスタリング手段と、
ＧＰＳの測位容易性の不安定化を吸収するべく、クラスタ毎に、滞留位置が１回でも存在する滞在日数を、日属性に基づいて計数する滞在率算出手段と、
滞在日数が多いクラスタから順に、日属性に基づいて、当該クラスタの滞在特性を推定する滞在特性推定手段と
してコンピュータを機能させることを特徴とする。
【００１６】
本発明によれば、ＧＰＳの測位電波を受信する測位部を有し、ユーザの移動に伴って所持される携帯端末であって、測位部を所定時間周期で起動し、測位位置を取得し、時刻に応じた測位位置を蓄積する携帯端末について、ユーザの住所及び居所を推定する方法であって、
時系列の測位位置の移動に応じて、移動中位置以外の滞留位置のみを抽出する第１のステップと、
複数の滞留位置を、空間的なクラスタに区分する第２のステップと、
ＧＰＳの測位容易性の不安定化を吸収するべく、クラスタ毎に、滞留位置が１回でも存在する滞在日数を、日属性に基づいて計数する第３のステップと、
滞在日数が多いクラスタから順に、日属性に基づいて、当該クラスタの滞在特性を推定する第４のステップと
を有することを特徴とする。
【発明の効果】
【００１７】
本発明の携帯端末、サーバ、プログラム及び方法によれば、携帯端末に搭載された測位機能のみを用いて、ユーザの移動に伴う住所／居所を自動的に推定することができる。
【図面の簡単な説明】
【００１８】
【図１】所定時間周期毎の測位位置における空間的なプロット図である。
【図２】本発明における携帯端末の機能構成図である。
【図３】多数の測位位置から平滑化位置を抽出する説明図である。
【図４】複数の平滑化位置から滞留位置を抽出する説明図である。
【図５】複数の滞留位置を、空間的なクラスタに区分した説明図である。
【図６】クラスタ毎に、滞在日数を計数した説明図である。
【図７】多数の滞留位置から重心位置を算出する説明図である。
【図８】携帯端末及び滞在特性サーバを有するシステム構成図である。
【発明を実施するための形態】
【００１９】
以下、本発明の実施の形態について、図面を用いて詳細に説明する。
【００２０】
図１は、所定時間周期毎の測位位置における空間的なプロット図である。
【００２１】
ユーザの移動に伴って所持される携帯端末（例えば携帯電話機）は、ＧＰＳ衛星からの測位電波を受信することによって、位置を取得することができる。図１によれば、ユーザは、自宅の住所が「埼玉県ふじみ野市」にあり、職場の居所が「東京都港区飯田橋」にある。そのユーザは、自宅と職場との間を、池袋を経由して通勤している。また、そのユーザは、訪問先として「東京都千代田区大手町」へも頻繁に移動している。このような場合、ユーザに所持される携帯端末は、所定時間周期（例えば１５分間隔）で測位することによって、時刻に応じた多数の測位位置が取得される。尚、測位位置は、ＧＰＳを用いた場合、「緯度・経度」として表される。
【００２２】
本発明について、「住所」とは、ユーザ各人の生活の本拠をいい、例えば生活の中心となっている自宅の所在地を意味する。また、「居所」とは、継続して居るものの、生活の本拠というほどその場所との結びつきが強くない場所をいい、例えば職場又は学校の所在地をいう。
【００２３】
図２は、本発明における携帯端末の機能構成図である。
【００２４】
図２によれば、携帯端末１は、測位部１００と、位置取得部１１１と、位置蓄積部１１２と、アプリケーション処理部１１３とを有する。測位部１００は、ＧＰＳ衛星から送信される測位電波を受信し、現在位置を特定する。位置取得部１１１は、測位部１００を所定時間周期で起動し且つ測位位置を取得する。位置蓄積部１１２は、時刻に応じた測位位置を蓄積する。アプリケーション処理部１１３は、本発明によって推定されたユーザに滞在特性に応じて、様々なサービスを実行する。
【００２５】
また、本発明における携帯端末１は、滞留位置抽出部１２１と、クラスタリング部１２２と、滞在率算出部１２３と、滞在特性推定部１２４と、重心位置決定部１２５と、地域名称記憶部１２６と、名前解決部１２７とを更に有する。これら機能構成部は、携帯端末に搭載されたコンピュータを機能させるプログラムを実行することによって実現される。
【００２６】
［滞留位置抽出部］
滞留位置抽出部１２１は、時系列の測位位置の移動に応じて、移動中位置以外の滞留位置のみを抽出する。具体的には、ユーザが自動車や電車等における移動中位置を取り除き、自宅や職場における滞留位置を抽出する。
【００２７】
図３は、多数の測位位置から平滑化位置を抽出する説明図である。
【００２８】
図３（ａ）によれば、所定時間周期Ｉ毎における多数の測位位置ｐ(t)(t=1〜6)が表されている。また、図３（ｂ）によれば、現時刻tにおける平滑化位置ｐａ(t)(t=3〜6)が表されている。平滑化位置ｐａ(t)は、現時刻tから直近過去所定時間Ｔ以内に取得された複数の測位位置ｐ(t)〜ｐ(t-n)を平滑化したものである。ＧＰＳによって得られた測位位置は誤差を含むために、実際に移動していない場合であっても、測位位置が変位している場合がある。そのために、平滑化位置ｐａ(t)を算出することによって、誤差分を吸収することができる。
【００２９】
図３（ｂ）によれば、以下のように表される。
Ｉ：所定時間周期
Ｔ：平滑化のための直近過去所定時間
ｐ(t)：時刻tにおける測位位置    P(t)＝ベクトル(latt, lont)
ｎ(t)：時刻tにおける直近過去所定時間Ｔ以内に測位に成功したデータ数
０≦ｎ(t)≦（Ｔ／Ｉ）
ｐａ(t)：ｎ(t)＞０である場合、時刻tにおける平滑化位置
ｐａ(t)＝(１／ｎ(t))Σt-Ttｐ(t)
【００３０】
図４は、複数の平滑化位置から滞留位置を抽出する説明図である。
【００３１】
図４（ａ）によれば、直近過去の平滑化位置ｐａ(t-1)から、現時刻tの平滑化位置ｐａ(t)への変位距離を算出する。変位距離は、例えば以下のように、ユークリッド距離によって算出される。
ｄ（ｐａ(t)，ｐａ(t-1)）
＝√（(lata(t)-lata(t-1))2＋(lona(t)-lona(t-1))2）
ｄ（a，b）：点aと点bとの間のユークリッド距離
【００３２】
次に、平滑化位置ｐａ間毎の変位距離に対する変位時間から、移動速度を算出する。
移動速度は、例えば以下の式によって算出される。
ｖ(t)＝ｄ（ｐａ(t)，ｐａ(t-1)）／変位時間
変位時間：ｐａ(t)の時刻tと、ｐａ(t-1)の時刻t-1との時間差
【００３３】
図４（ｂ）によれば、最後に、移動速度ｖ(t)が、所定閾値Ｖth以下となる平滑化位置ｐａのみを抽出する。即ち、移動速度ｖ(t)が、所定閾値Ｖthよりも速い平滑化位置を除去することによって、移動中位置を取り除く。そして、移動速度が所定閾値以下となる滞留位置のみを残すことができる。
【００３４】
［クラスタリング部］
クラスタリング部１２２は、複数の滞留位置の群を、空間的な０個以上のクラスタに区分する。
【００３５】
図５は、複数の滞留位置を、空間的なクラスタに区分した説明図である。
【００３６】
例えば、図５によれば、４つのクラスタに区分されている。複数の滞留位置をクラスタリングするために、既知のアルゴリズム、例えばＤＢＳＣＡＮ(Density-Based Spatial Clustering)（例えば非特許文献２参照）を用いることができる。k-meansを用いてもよいが、クラスタ数kを事前に与える必要がある。これに対し、ＤＢＳＣＡＮの場合、クラスタ数を事前に与えことなく、比較的に直感的なクラスタを抽出できる点で適している。
【００３７】
ＤＢＳＣＡＮは、超球状ではない任意形状のクラスタの抽出を目的としたクラスタリング方法であって、２つの点における直接密度到達可能(directly density-rechable)を導出するものである。あるseed点から、直接密度到達可能な関係を推移的に辿って、到達可能な極大集合を１つのクラスタとして抽出する。
【００３８】
ＤＢＳＣＡＮによれば、Eps(epsilon、距離)及びMinPts(minimum points、最低ポイント)の２つのパラメータを用いる。例えばEps=0.0005及びMinPts=10を事前に与えるものであってもよい。また、これらパラメータを、測位の所定時間周期や最終的な住所居所推定の結果によって変更するものであってもよい。
【００３９】
［滞在率算出部］
滞在率算出部１２３は、クラスタ毎に、滞留位置が１回でも存在する滞在日数を用いて、日属性に基づく滞在率を算出する。
【００４０】
図６は、クラスタ毎に、滞在日数を計数した説明図である。
【００４１】
滞在率算出部１２３は、クラスタ毎に、日属性（全日、平日、及び休日）に基づいて、滞留位置が１回でも存在する滞在日数を計数する。
ｎｄ(c)：クラスタcにおける全日滞在日数
ｎｄw(c)：クラスタcにおける平日滞在日数
ｎｄh(c)：クラスタcにおける休日滞在日数
【００４２】
ここで、日数ベースで、滞留位置をカウントする理由について説明する。本発明によれば、測位位置が１日に１回以上存在すれば、滞在日数を１増分する。この場合、測位位置が少なくても、全ての日で１回以上その測位位置が存在すれば、滞在率は１００％となる。逆に、いずれかの特定の１日のみで、測位位置が多く存在しても、その滞在率は低くなる。
【００４３】
滞在率を「測位位置の数」で見ないのは、ＧＰＳの場合、測位位置の数は、滞在時間の長さではなく、測位容易性（測位のし易さ）に依存してしまうためである。
【００４４】
また、滞在率を「滞在時間」で見ないのは、３つの理由がある。
【００４５】
第１の理由は、滞在時間を導出するために、クラスタへの入場時刻と出場時刻とを追跡する必要があり、計算量が多くなるからである。入場時刻とは、クラスタ外からクラスタ内へ測位位置が変位した時刻である。出場時刻とは、クラスタ内からクラスタ外へ測位位置が変位した時刻である。
【００４６】
第２の理由は、ＧＰＳの測位容易性によって、クラスタごとに滞在時間の厳密さに、ばらつきが生じるからである。即ち、測位のし易いクラスタの滞在時間は正確に求められるが、測位のし難いクラスタの滞在率は極めて荒くなる。
【００４７】
第３の理由は、クラスタの面積の大きさが、滞在時間に対して影響を与えやすいからである。即ち、大きい面積のクラスタほど滞在時間が長くなる。従って、クラスタリング精度が、滞在時間の精度に影響を与えてしまう。例えば、同じ測位位置に対するクラスタリングの結果が、パラメータによって、図６のクラスタＡのような場合もあれば、クラスタＢのような場合もあり得る。このとき、滞在時間はクラスタの面積が大きいクラスタＢの方が長くなるが、滞在日数は等しい。
【００４８】
このような理由から、日数ベースで滞在率を算出することによって、ＧＰＳの測位容易性の不安定さや、クラスタリングの精度の不安定さを吸収し、高い精度で、滞在場所を推定することが可能となる。
【００４９】
ここで、平日／休日の定義は、カレンダ情報に基づいて土曜日及び日曜日に設定されたものであってもよいし、ユーザに利用態様に応じて設定されたものであってもよい。携帯端末が、携帯電話機である場合、アプリケーションとして搭載されたカレンダ／スケジューラの機能を用いることもできる。
【００５０】
次に、滞在率算出部１２３は、クラスタ毎に、日属性に基づいて、滞在日数から滞在率を算出する。滞在率は、例えば以下のように算出される。
Ｄ：全日日数
Ｄw：平日日数
Ｄh：休日日数
Ｒ(c)：クラスタcにおける全日滞在率
Ｒ(c)＝ｎｄ(c)／Ｄ
Ｒw(c)：クラスタcにおける平日滞在率
Ｒw(c)＝ｎｄw(c)／Ｄw
Ｒh(c)：クラスタcにおける休日滞在率
Ｒh(c)＝ｎｄh(c)／Ｄh
【００５１】
［滞在特性推定部］
滞在特性推定部１２４は、日属性に基づく滞在率が高いクラスタについて、当該クラスタの滞在特性を推定する。滞在特性は、自宅のような「住所」、及び、職場・学校のような「居所」に基づいている。推定方法としては、例えば以下の３通りがある。
【００５２】
第１に、平日休日を問わず、全日滞在率（滞在日数／記録日数）が高い第１のクラスタを「住所」と規定し、第１のクラスタ以外に平日滞在率が高い第２のクラスタを「居所」と規定するものであってもよい。
【００５３】
第２に、全日滞在率が高い２つのクラスタＡ及びクラスタＢを選択し、平日を基準とした以下の関係に基づいて、「住所／居所」を規定するものであってもよい。
ＲwA＞ＲwBの場合：クラスタＡ「住所」、クラスタＢ「居所」
ＲwA＜ＲwBの場合：クラスタＡ「居所」、クラスタＢ「住所」
ＲwA＝ＲwBの場合：休日滞在率が高い方のクラスタを「住所」
【００５４】
第３の、全日滞在率が高い２つのクラスタＡ及びクラスタＢを選択し、休日を基準とした以下の関係に基づいて、「住所／居所」を規定するものであってもよい。
ＲhA＞ＲhBの場合：クラスタＡ「住所」、クラスタＢ「居所」
ＲhA＜ＲhBの場合：クラスタＡ「居所」、クラスタＢ「住所」
ＲhA＝ＲhBの場合：平日滞在率が高い方のクラスタを「居所」
【００５５】
例えば、図６によれば、クラスタ１が「住所」として規定され、クラスタ３が「居所」として規定されている。
【００５６】
［重心位置決定部］
重心位置決定部１２５は、滞在特性推定部１２４によって抽出された第１のクラスタ及び第２のクラスタについて、多数の滞留位置から重心位置を算出する。
【００５７】
図７は、多数の滞留位置から重心位置を算出する説明図である。
【００５８】
重心位置を算出することによって、クラスタリングの精度の不安定さを吸収し、比較的正確な住所／居所の位置を推定することができる。また、クラスタの面積が広くても、その重心位置はほぼ近い位置とる。例えば、同じ測位位置に対するクラスタリングの結果が、パラメータによって、図７のクラスタＡのような場合もあれば、クラスタＢのような場合もあり得る。このとき、クラスタの示す領域は異なるが、重心位置は大きく変わらない。
【００５９】
あるクラスタcの重心位置ｐｃは、そのクラスタに属する滞留位置の群の平均値を算出すればよい。クラスタcにｎc個の滞留位置（平滑化位置）が含まれている場合、以下の式によって算出される。
ｐｃx＝(１／ｎc)Σpai∈cｐａi
【００６０】
そして、重心位置決定部１２５は、第１のクラスタの重心位置を「住所」位置として規定し、第２のクラスタの重心位置を「居所」位置として規定する。
【００６１】
［地域名称記憶部］
地域名称記憶部１２６は、地域名称に対応付けて位置を記憶する。例えば、駅名を地域名称として、その駅名に対応付けて緯度経度を記憶したデータベースもある（例えば非特許文献３参照）。非特許文献３によれば、例えば以下のようなデータを記憶している。
鉄道概要コード  rr_cd        整数２桁
路線コード      line_cd      整数５桁  鉄道コード＋エリアコード＋連番
駅コード        station_cd   整数７7桁  路線コード＋連番
・・・・
鉄道概要名      rr_name      文字列６４byte
路線名          line_name    文字列１２８byte
駅名            station_name 文字列１２８byte
都道府県コード  pref_cd      整数２桁  北海道〜沖縄＋その他
緯度            lat          数値
経度            lon          数値
【００６２】
［名前解決部］
名前解決部１２７は、地域名称記憶部１２６を用いて、住所位置及び居所位置に最も近い「地域名称」を検索し、当該「地域名称」を、住所位置の住所及び居所位置の居所と規定する。具体的には、駅データを用いた場合、住所位置の緯度・経度に最も近い緯度・経度となる駅名を、住所位置の地域名称とする。また、居所位置の緯度・経度に最も近い緯度・経度となる駅名を、居所位置の地域名称とする。地域名称を規定することによって、重心位置精度の不安定さを吸収することができる。また、これら地域名称は、アプリケーション処理部１１３によって様々に利用することができる。
【００６３】
例えば、図７によれば、住所位置に最も近い「ふじみ野」駅を、住所位置の地域名称と規定される。また、居所位置に最も近い「飯田橋」駅を、居所位置の地域名称と規定される。従って、居所位置から離れた「水道橋」駅及び「神保町」駅は、居所位置の地域名称として選択されない。重心位置が正しい住所位置からずれている場合であっても、地域名称は正しい名称をつけることができる。
【００６４】
図８は、携帯端末及び滞在特性サーバを有するシステム構成図である。
【００６５】
図８によれば、携帯端末１は、ユーザの移動に伴って所持され、所定時間周期で測位位置を取得する。取得された測位位置は、位置情報送信部１１４から、ネットワークを介して、滞在特性推定サーバ２へ送信される。滞在特性推定サーバ２は、前述した滞留点抽出部１２１〜名前解決部１２７の機能を有する。これによって、サーバは、携帯端末の移動における滞在特性を推定し、様々なサービス情報を携帯端末へ送信することができる。
【００６６】
以上、詳細に説明したように、本発明の携帯端末、サーバ、プログラム及び方法によれば、携帯端末に搭載された測位機能のみを用いて、ユーザの移動に伴う住所／居所を自動的に推定することができる。
【００６７】
前述した本発明の種々の実施形態について、本発明の技術思想及び見地の範囲の種々の変更、修正及び省略は、当業者によれば容易に行うことができる。前述の説明はあくまで例であって、何ら制約しようとするものではない。本発明は、特許請求の範囲及びその均等物として限定するものにのみ制約される。
【符号の説明】
【００６８】
１  携帯端末
１００  測位部
１０１  無線通信インタフェース
１０２  通信インタフェース
１１１  位置取得部
１１２  位置蓄積部
１１３  アプリケーション処理部
１１４  位置情報送信部
１１５  位置情報受信部
１２１  滞留位置抽出部
１２２  クラスタリング部
１２３  滞在率算出部
１２４  滞在特性推定部
１２５  重心位置決定部
１２６  地域名称記憶部
１２７  名前解決部
２  滞在特性推定サーバ
【図１】
0005536485.tif 000002
【図２】
0005536485.tif 000003
【図３】
0005536485.tif 000004
【図４】
0005536485.tif 000005
【図５】
0005536485.tif 000006
【図６】
0005536485.tif 000007
【図７】
0005536485.tif 000008
【図８】
0005536485.tif 000009