CN102959996A - 终端数量估计装置及终端数量估计方法 - Google Patents
终端数量估计装置及终端数量估计方法 Download PDFInfo
- Publication number
- CN102959996A CN102959996A CN2011800302519A CN201180030251A CN102959996A CN 102959996 A CN102959996 A CN 102959996A CN 2011800302519 A CN2011800302519 A CN 2011800302519A CN 201180030251 A CN201180030251 A CN 201180030251A CN 102959996 A CN102959996 A CN 102959996A
- Authority
- CN
- China
- Prior art keywords
- position data
- observation
- unit
- information
- terminal quantity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W64/00—Locating users or terminals or network equipment for network management purposes, e.g. mobility management
- H04W64/006—Locating users or terminals or network equipment for network management purposes, e.g. mobility management with additional information processing, e.g. for direction or speed determination
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W64/00—Locating users or terminals or network equipment for network management purposes, e.g. mobility management
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W24/00—Supervisory, monitoring or testing arrangements
- H04W24/10—Scheduling measurement reports ; Arrangements for measurement reports
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
- H04W4/021—Services related to particular areas, e.g. point of interest [POI] services, venue services or geofences
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
- H04W4/029—Location-based management or tracking services
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W8/00—Network data management
- H04W8/02—Processing of mobility data, e.g. registration information at HLR [Home Location Register] or VLR [Visitor Location Register]; Transfer of mobility data, e.g. between HLR, VLR or external networks
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L67/00—Network arrangements or protocols for supporting network services or applications
- H04L67/50—Network services
- H04L67/52—Network services specially adapted for the location of the user terminal
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Business, Economics & Management (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Quality & Reliability (AREA)
- Physics & Mathematics (AREA)
- Development Economics (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Educational Administration (AREA)
- Tourism & Hospitality (AREA)
- Game Theory and Decision Science (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- Mobile Radio Communication Systems (AREA)
- Telephonic Communication Services (AREA)
Abstract
一种终端数量估计装置,其具有:位置数据取得部,其取得包含识别信息、位置信息和位置取得时刻信息的位置数据;前后位置数据取得部,其对于某个第1位置数据,取得包括相同的识别信息的位置数据中、该第1位置数据的紧前面的第2位置数据和紧后面的第3位置数据的位置取得时刻信息;特征量计算部,其根据第1~第3位置数据的位置取得时刻信息中两个以上的位置取得时刻信息,计算有关第1位置数据的特征量;观测对象取得部,其取得包括观测开始时刻以后而且是观测结束时刻以前的位置取得时刻信息、而且包含与观测区域信息相对应的位置信息的位置数据;以及终端数量估计部,其根据有关观测对象位置数据的特征量及观测期间长度,估计在观测期间中位于观测区域的终端数量。
Description
技术领域
本发明涉及终端数量估计装置及终端数量估计方法,使用从便携式终端(例如移动电话)的网络设备得到的有关便携式终端的位置信息,估计位于某个区域内的终端数量。
背景技术
在移动电话运营商的网络设备中,生成移动电话的位置数据或利用者的属性数据等的运用数据,以向移动电话的利用者提供电气通信服务。通过对这些运用数据实施累计等统计处理,能够得到“人口分布”、“人口变动”、“人口构成”等有关人口动态的估计值。其中,“人口分布”是分布于各个区域的人口的相关信息,“人口变动”是某个区域中的沿着时间轴的人口变动的相关信息,“人口构成”是例如按照性别或年龄段等区分的人口分布或人口变动的相关信息。
作为上述的位置数据例如有位置登记信号。该位置登记信号是从移动电话向归属基站大致周期性地发送的信号,在某个基站接收到某个移动电话的位置登记信号的情况下,能够估计到该移动电话在接收时刻位于该基站的电波到达范围即基站扇区中。
另外,作为位置数据的另一个示例有GPS信息。该GPS信息是从移动电话向归属基站周期性地发送的信息、或者响应于来自终端的操作或者移动电话网络的请求而发送的GPS测位结果的相关信息。同样地,在接收到GPS信息的时刻,能够估计到该移动电话位于GPS测位结果所表示的位置的周围(例如参照专利文献1)。
如果能够根据诸如上述的位置数据的观测结果来估计位于某个地理区域的移动电话的台数(终端数量),通过再考虑移动电话的签约率等,预计能够得到有关上述的各种人口动态的估计值。
现有技术文献
专利文献
专利文献1:日本特开2003-44969号公报
发明内容
发明要解决的问题
可是,根据诸如上述的位置数据准确估计终端数量并非易事。这是因为位置登记信号或GPS信息等终端的位置数据不是始终发送给移动电话网络,而是间隔某种程度的时间间隔来发送的以及发送的时间间隔不固定。
例如,分析使用位置登记信号作为位置数据来估计终端数量的情况。如果能够假设位置登记信号是完全被周期性地发送的,则某个扇区的在设定的观测时间内接收的位置登记信号的数量与扇区内的终端数量成比例。但是,实际上尽管位置登记信号是以按照例如移动电话内的定时器来周期性地发送为基本,但是无论定时器的状态如何,均借助跨越某个特定的基站扇区之间的契机来进行位置登记信号的发送,有时由于通话或服务区外等的影响而使得发送延迟。另外,对于GPS信息也一样,由于服务区外或终端的操作等各种影响,发送接收的周期不固定。
另外,分析并用位置登记信号和GPS信息作为位置数据的情况。在这种情况下,能够在估计终端数量时利用的信息量增加,因而预计能够实现精度更高的终端数量估计。但是,在这种情况下,并用两种信息的位置数据的发送接收频次与分别独立使用时相比,随机程度进一步增大,如果假设发送接收的周期为固定来估计终端数量,则担忧估计精度反而恶化。
因此,为了高精度地估计终端数量,需要不仅单纯地统计接收信号的数量,而且要考虑接收信号的接收间隔的变动。
本发明正是鉴于上述情况而提出的,其目的在于,在使用位置数据估计终端数量时,对接收间隔的变动所产生的影响进行校正,并且高精度地估计终端数量。
用于解决问题的手段
本发明的一个方面的终端数量估计装置的特征在于,该终端数量估计装置具有:位置数据取得单元,其取得包含识别便携式终端的识别信息、有关便携式终端的位置的位置信息和取得位置信息的位置取得时刻信息的位置数据;前后位置数据取得单元,其对于某个第1位置数据,取得包括与该第1位置数据相同的识别信息的位置数据中、该第1位置数据的紧前面的位置数据即第2位置数据的位置取得时刻信息、以及该第1位置数据的紧后面的位置数据即第3位置数据的位置取得时刻信息;特征量计算单元,其根据第1位置数据的位置取得时刻信息、第2位置数据的位置取得时刻信息及第3位置数据的位置取得时刻信息中两个以上的位置取得时刻信息,计算有关第1位置数据的特征量;观测对象取得单元,其取得包括有关待观测的观测期间的观测开始时刻以后而且是观测结束时刻以前的位置取得时刻信息、而且包括与有关待观测的观测区域的观测区域信息相对应的位置信息在内的一个或多个位置数据,作为观测对象位置数据;以及终端数量估计单元,其根据有关观测对象位置数据的特征量以及观测开始时刻与观测结束时刻之差即观测期间长度,估计在观测期间中位于观测区域内的终端数量。另外,“特征量”是指与有关便携式终端生成的位置数据的估计生成密度对应的信息,其中的“估计生成密度”是指生成该位置数据的便携式终端在该位置数据的生成时刻(相当于上述的位置取得时刻)左右于每单位时间内生成的位置数据的数量的估计值。在上述的终端数量估计装置中也可以是,特征量计算单元计算第2位置数据的位置取得时刻与第3位置数据的位置取得时刻之差,作为有关第1位置数据的特征量,终端数量估计单元将有关观测对象位置数据的特征量的总和除以观测期间长度的2倍,将相除得到的数值估计为终端数量。具体原理将在后面叙述,特征量计算单元计算第2、第3位置数据的位置取得时刻之差作为有关第1位置数据的特征量,终端数量估计单元将有关观测对象位置数据的特征量的总和除以观测期间长度的2倍,将相除得到的数值估计为终端数量,由此能够对接收间隔的变动所产生的影响进行校正,并且高精度地估计终端数量。
另外,也可以是,在第1位置数据的位置取得时刻与第2位置数据的位置取得时刻之差大于预定值的情况下,特征量计算单元将从第1位置数据的位置取得时刻起向过去追溯预定的时间得到的时刻用作第2位置数据的位置取得时刻,来计算有关第1位置数据的特征量。同样,也可以是,在第1位置数据的位置取得时刻与第3位置数据的位置取得时刻之差大于预定值的情况下,特征量计算单元将从第1位置数据的位置取得时刻起向未来前进预定的时间得到的时刻用作第3位置数据的位置取得时刻,来计算有关第1位置数据的特征量。通过使特征量计算单元按照以上所述进行动作,在位置数据的取得时间间隔由于便携式终端位于服务区外或便携式终端的电源切断等而变得异常长时,能够防止由于该异常长的取得时间间隔而造成的影响过大。
另外,也可以是,特征量计算单元判别第1位置数据是否包括由于跨越位置登记区域边界而生成的位置登记信息、以及第3位置数据是否包括由于跨越位置登记区域边界而生成的位置登记信息,根据第1位置数据是否包括由于跨越位置登记区域边界而生成的位置登记信息的判别结果、和第3位置数据是否包括由于跨越位置登记区域边界而生成的位置登记信息的判别结果,使用第1位置数据的位置取得时刻信息、第2位置数据的位置取得时刻信息及第3位置数据的位置取得时刻信息中的两个以上的位置取得时刻信息,计算有关第1位置数据的特征量。在这种情况下,具体原理将在后面叙述,通过考虑与由于跨越位置登记区域边界而生成的位置登记信息相关的生成定时的特性,能够得到精度较高的特征量。另外,所谓“由于跨越位置登记区域边界而生成的位置登记信息”,是指由于便携式终端跨越位置登记区域边界而生成的位置登记信息。
更具体地讲,也可以是,特征量计算单元在第1位置数据包括由于跨越位置登记区域边界而生成的位置登记信息的情况下,将第1位置数据的位置取得时刻设定为第1变量,在第1位置数据不包括由于跨越位置登记区域边界而生成的位置登记信息的情况下,将第1位置数据的位置取得时刻与第2位置数据的位置取得时刻之间的中点时刻设定为第1变量,在第3位置数据包括由于跨越位置登记区域边界而生成的位置登记信息的情况下,将第3位置数据的位置取得时刻设定为第2变量,在第3位置数据不包括由于跨越位置登记区域边界而生成的位置登记信息的情况下,将第1位置数据的位置取得时刻与第3位置数据的位置取得时刻之间的中点时刻设定为第2变量,根据所设定的第1变量与第2变量之差来计算有关第1位置数据的特征量。
另外,也可以是,在第1位置数据的位置取得时刻与第1变量之差大于预定值的情况下,特征量计算单元将从第1位置数据的位置取得时刻起向过去追溯预定的时间得到的时刻用作第1变量,来计算有关第1位置数据的特征量。同样,也可以是,在第1位置数据的位置取得时刻与第2变量之差大于预定值的情况下,特征量计算单元将从第1位置数据的位置取得时刻起向未来前进预定的时间得到的时刻用作第2变量,来计算有关第1位置数据的特征量。通过使特征量计算单元按照以上所述进行动作,在位置数据的取得时间间隔由于便携式终端位于服务区外或便携式终端的电源被切断等而变得异常长时,能够防止由于该异常长的取得时间间隔而造成的影响过大。
另外,在终端数量估计装置中,作为特征量的计算对象,也可以将由观测对象取得单元取得的观测对象位置数据作为对象,还可以将由位置数据取得单元取得的所有位置数据作为对象。其中,在将由观测对象取得单元取得的观测对象位置数据作为对象的情况下,前后位置数据取得单元将由观测对象取得单元取得的各个观测对象位置数据作为第1位置数据,并取得有关该第1位置数据的第2位置数据的位置取得时刻信息、以及有关该第1位置数据的第3位置数据的位置取得时刻信息,特征量计算单元计算有关各个观测对象位置数据的特征量,终端数量估计单元使用计算得到的有关观测对象位置数据的特征量来估计终端数量。
另一方面,在将由位置数据取得单元取得的所有位置数据作为对象的情况下,前后位置数据取得单元将由位置数据取得单元取得的所有位置数据的各个作为第1位置数据,并取得有关该第1位置数据的第2位置数据的位置取得时刻信息、以及有关该第1位置数据的第3位置数据的位置取得时刻信息,特征量计算单元计算有关所有位置数据各自的特征量,终端数量估计单元使用计算得到的有关所有位置数据各自的特征量中有关观测对象位置数据的特征量来估计终端数量。
本发明的一个方面的终端数量估计装置也可以构成为,该终端数量估计装置具有:位置数据取得单元,其取得包含识别便携式终端的识别信息、有关便携式终端的位置的位置信息和取得位置信息的位置取得时刻信息的位置数据;前后位置数据取得单元,其对于某个第1位置数据,取得包括与该第1位置数据相同的识别信息的位置数据中、该第1位置数据的紧前面的位置数据即第2位置数据的位置取得时刻信息、以及该第1位置数据的紧后面的位置数据即第3位置数据的位置取得时刻信息;特征量计算单元,其根据至少第2位置数据的位置取得时刻信息及第3位置数据的位置取得时刻信息,计算有关第1位置数据的特征量;观测对象取得单元,其取得包括有关待观测的观测期间的观测开始时刻以后而且是观测结束时刻以前的位置取得时刻信息、而且包括与有关待观测的观测区域的观测区域信息相对应的位置信息在内的一个或多个位置数据,作为观测对象位置数据;以及终端数量估计单元,其根据有关观测对象位置数据的特征量以及观测开始时刻与观测结束时刻之差即观测期间长度,估计在观测期间中位于观测区域内的终端数量。
另外,也可以是,终端数量估计装置还具有存储用于将终端数量变换为人口的放大系数的放大系数存储单元,终端数量估计单元根据有关观测对象位置数据的特征量、观测期间长度及放大系数,估计在观测期间中位于观测区域的人口、以及每人口估计单位的人口中的至少一方,其中该人口估计单位是估计人口的单位。关于上述的“人口估计单位”,例如可以列举属性、场所、时间段等。另外,放大系数可以使用存储在放大系数存储单元中的放大系数,也可以按照下面所述来导出放大系数。作为一例,放大系数能够使用“在服务区率与终端的普及率之积(即,在服务区率与人口的比率)”的倒数。其中的“在服务区率”是指在服务区数量与签约台数的比率,“普及率”是指签约台数与人口的比率。优选按照上述的每人口估计单位来导出这种放大系数,但不是必须的。
另外,也可以是,例如使用按照下面所述根据特征量及观测期间长度而估计出的终端数量(在服务区数量)来导出放大系数。即,从位置数据求出特征量,根据特征量及观测期间长度来累计每种放大系数计算单位的终端数量,由此得到用户数量金字塔数据,并且作为统计数据(例如居民基本登记册等)来取得预先求出的相同的每种放大系数计算单位的人口金字塔数据。并且,在用户数量金字塔数据和人口金字塔数据中计算每种放大系数计算单位的位置数据的取得率(即在服务区数量/人口)。在此得到的“位置数据的取得率(即在服务区数量/人口)”相当于前述的“在服务区率与终端的普及率之积”。能够导出这样得到的“位置数据的取得率”的倒数作为放大系数。另外,关于计算放大系数的放大系数计算单位,作为一例,也可以采用住所的都道府县、间隔5岁或10岁的年龄段、男女、作为时间段的1小时等,还可以采用将其中两个以上的单位进行组合的单位。例如,在将放大系数计算单位设为“居住在东京都的二十几岁的男性”的情况下,通过提取日本全国的与居住在东京都的(即,用户属性中的住所信息是东京都)二十几岁的男性对应的位置数据并累计终端数量,能够得到用户数量金字塔数据,并且从统计数据取得与居住在东京都的二十几岁的男性相关的人口金字塔数据。另外,在得到上述用户数量金字塔数据时,关于条件“居住在东京都”,不是仅提取在东京都服务区的用户的位置数据,而是提取用户属性中的住所信息是东京都的位置数据。并且,根据用户数量金字塔数据和人口金字塔数据计算放大系数计算单位(此处指居住在东京都的二十几岁的男性)的位置数据的取得率(即,在服务区数量/人口),并导出所得到的“位置数据的取得率”的倒数作为放大系数。另外,在本申请中,将放大系数计算单位和人口估计单位作为相同单位进行说明,但这毕竟仅是一例,并不限于此。
另外,也可以是,终端数量估计装置还具有变换单元,该变换单元根据与观测区域不同的输出单位和观测区域的重合区域在观测区域中所占的面积比,将通过所述终端数量估计单元的估计而得到的每个观测区域的估计值变换为每输出单位的估计值。也可以是,在室内站的通信区域以及利用电波到达范围不同的频带的多个室外站的通信区域中有两个以上的通信区域重叠存在于地理上相同的观测区域的情况下,上述的变换单元对于重叠存在的各个通信区域,基于所述面积比变换为每输出单位的估计值、并针对各个所述通信区域对变换后的估计值进行合计,由此得到每输出单位的估计值。
也可以是,终端数量估计单元根据有关观测对象位置数据的特征量、观测期间长度、用于将终端数量变换为人口的放大系数、以及与观测区域不同的输出单位和观测区域的重合区域在观测区域中所占的面积比,区分为每输出单位及每人口估计单位来估计人口。
作为区分为每输出单位及每种人口估计单位来估计人口的方式,也可以是,在终端数量估计单元估计人口之前,使特征量、放大系数、面积比以及与该面积比相关的输出单位ID的组合与位置数据相关联,终端数量估计单元对于与相同输出单位ID相关联的位置数据计算(特征量×面积比×放大系数),并按人口估计单位来对所得到的每输出单位的(特征量×面积比×放大系数)进行累计,根据所得到的每输出单位且每人口估计单位的累计值和观测期间长度,估计每输出单位且每人口估计单位的人口。
另外,也可以是,终端数量估计装置还具有:观测期间取得单元,其取得包括观测开始时刻与观测结束时刻的组的观测期间信息;以及观测区域取得单元,其取得与一个或多个位置信息相对应的观测区域信息。
另外,也可以是,终端数量估计装置还具有输出所得到的估计值的输出单元。输出单元的输出形式能够采用表示人口分布的图、表示按时序人口变动的图、以及表示人口构成的图中的至少一种,输出单元的输出单位能够根据便携式终端的用户的属性、时间段、场所中的至少一种来设定。
另外,也可以是,终端数量估计装置还具有非识别化单元,其针对由位置数据取得单元取得的位置数据中所包含的识别信息,进行包含如下变换的非识别化处理,该变换是向基于单向性函数的不可逆编码的变换,非识别化单元在进行使用了便携式终端的用户的属性信息的处理的情况下,在该处理之前对属性信息进行非识别化处理。
另外,也可以是,终端数量估计装置还具有隐匿处理单元,在所得到的估计值被输出之前,该隐匿处理单元根据预先设定的基准对估计值进行隐匿处理。在这种情况下,也可以是,隐匿处理单元判定取得源终端数量是否小于判定为需要隐匿处理的基准值,在某个区域的位置数据的取得源终端数量小于基准值的情况下,将有关该区域的估计值隐匿,其中,所述取得源终端数量表示作为估计的基础的每个区域的位置数据是从哪台终端取得的。另外,上述的取得源终端数量表示去除了有关同一终端的重复的唯一的终端数量。
关于其中的隐匿方法,例如可以采用将估计值设为零的方法、利用预定的文字或记号(例如“×”等)表示估计值的方法等。另一方面,在某个区域的位置数据的取得源终端数量不小于基准值的情况下,隐匿处理单元可以不对估计值进行隐匿处理,也可以进行如下所述的圆整。即,隐匿处理单元根据在估计值输出中使用的多个级别中该区域的估计值所属的级别中的上限值、下限值、级幅以及该估计值,将该区域的估计值圆整为分别伴随有概率值的所述上限值及所述下限值,其中所述概率值对应于与所述上限值的差分以及与所述下限值的差分。
上述的终端数量估计装置的发明能够作为由终端数量估计装置执行的终端数量估计方法的发明,并发挥相同的作用及效果。具体地讲,根据是将观测对象位置数据作为对象来计算特征量、还是将所取得的所有位置数据作为对象来计算特征量,能够进行如下所述的记述。
本发明的一个方面的终端数量估计方法是由终端数量估计装置执行的终端数量估计方法,其特征在于,该终端数量估计方法包括:位置数据取得步骤,取得包含识别便携式终端的识别信息、有关便携式终端的位置的位置信息和取得位置信息的位置取得时刻信息的位置数据;观测对象取得步骤,取得包括有关待观测的观测期间的观测开始时刻以后而且是观测结束时刻以前的位置取得时刻信息、而且包括与有关待观测的观测区域的观测区域信息相对应的位置信息在内的一个或多个位置数据,作为观测对象位置数据;前后位置数据取得步骤,将所取得的各个观测对象位置数据作为第1位置数据,对于各个第1位置数据,取得包括与该第1位置数据相同的识别信息的位置数据中、该第1位置数据的紧前面的位置数据即第2位置数据的位置取得时刻信息、以及该第1位置数据的紧后面的位置数据即第3位置数据的位置取得时刻信息;特征量计算步骤,根据第1位置数据的位置取得时刻信息、第2位置数据的位置取得时刻信息及第3位置数据的位置取得时刻信息中两个以上的位置取得时刻信息,计算有关各个观测对象位置数据的特征量;以及终端数量估计步骤,根据有关通过计算得到的观测对象位置数据的特征量以及观测开始时刻与观测结束时刻之差即观测期间长度,估计在观测期间中位于观测区域内的终端数量。
本发明的一个方面的终端数量估计方法是由终端数量估计装置执行的终端数量估计方法,其特征在于,该终端数量估计方法包括:位置数据取得步骤,取得包含识别便携式终端的识别信息、有关便携式终端的位置的位置信息和取得位置信息的位置取得时刻信息的位置数据;前后位置数据取得步骤,将所取得的所有位置数据分别作为第1位置数据,对于各个第1位置数据,取得包括与该第1位置数据相同的识别信息的位置数据中、该第1位置数据的紧前面的位置数据即第2位置数据的位置取得时刻信息、以及该第1位置数据的紧后面的位置数据即第3位置数据的位置取得时刻信息;特征量计算步骤,根据第1位置数据的位置取得时刻信息、第2位置数据的位置取得时刻信息及第3位置数据的位置取得时刻信息中两个以上的位置取得时刻信息,计算有关所有位置数据各自的特征量;观测对象取得步骤,取得包括有关待观测的观测期间的观测开始时刻以后而且是观测结束时刻以前的位置取得时刻信息、而且包括与有关待观测的观测区域的观测区域信息相对应的位置信息在内的一个或多个位置数据,作为观测对象位置数据;以及终端数量估计步骤,根据有关通过计算得到的所有位置数据各自的特征量中有关观测对象位置数据的特征量、以及观测开始时刻与观测结束时刻之差即观测期间长度,估计在观测期间中位于观测区域内的终端数量。
本发明的一个方面的终端数量估计方法是由终端数量估计装置执行的终端数量估计方法,其特征在于,该终端数量估计方法包括:位置数据取得步骤,取得包含识别便携式终端的识别信息、有关便携式终端的位置的位置信息和取得位置信息的位置取得时刻信息的位置数据;观测对象取得步骤,取得包括有关待观测的观测期间的观测开始时刻以后而且是观测结束时刻以前的位置取得时刻信息、而且包括与有关待观测的观测区域的观测区域信息相对应的位置信息在内的一个或多个位置数据,作为观测对象位置数据;前后位置数据取得步骤,将所取得的各个观测对象位置数据作为第1位置数据,对于各个第1位置数据取得包括与该第1位置数据相同的识别信息的位置数据中、该第1位置数据的紧前面的位置数据即第2位置数据的位置取得时刻信息、以及该第1位置数据的紧后面的位置数据即第3位置数据的位置取得时刻信息;特征量计算步骤,根据至少所述第2位置数据的位置取得时刻信息及所述第3位置数据的位置取得时刻信息,计算有关各个观测对象位置数据的特征量;以及终端数量估计步骤,根据有关通过计算得到的观测对象位置数据的特征量以及观测开始时刻与观测结束时刻之差即观测期间长度,估计在观测期间中位于观测区域内的终端数量。
本发明的一个方面的终端数量估计方法是由终端数量估计装置执行的终端数量估计方法,其特征在于,该终端数量估计方法包括:位置数据取得步骤,取得包含识别便携式终端的识别信息、有关便携式终端的位置的位置信息和取得位置信息的位置取得时刻信息的位置数据;前后位置数据取得步骤,将所取得的所有位置数据分别作为第1位置数据,对于各个第1位置数据取得包括与该第1位置数据相同的识别信息的位置数据中、该第1位置数据的紧前面的位置数据即第2位置数据的位置取得时刻信息、以及该第1位置数据的紧后面的位置数据即第3位置数据的位置取得时刻信息;特征量计算步骤,根据至少第2位置数据的位置取得时刻信息及第3位置数据的位置取得时刻信息,计算有关所有位置数据各自的特征量;观测对象取得步骤,取得包括有关待观测的观测期间的观测开始时刻以后而且是观测结束时刻以前的位置取得时刻信息、而且包括与有关待观测的观测区域的观测区域信息相对应的位置信息在内的一个或多个位置数据,作为观测对象位置数据;以及终端数量估计步骤,根据有关通过计算得到的所有位置数据各自的特征量中有关观测对象位置数据的特征量、以及观测开始时刻与观测结束时刻之差即观测期间长度,估计在观测期间中位于观测区域内的终端数量。
发明效果
根据本发明,能够对接收间隔的变动所产生的影响进行校正,并且高精度地估计终端数量。
附图说明
图1是示出第1~第6实施方式的通信***的***结构的图。
图2是示出第1实施方式的终端数量估计装置的结构的图。
图3是用于说明终端数量估计的第1思路的图。
图4是用于说明终端数量估计的第1计算方法的图。
图5是示出第1实施方式的终端数量估计处理的流程图。
图6是示出第1、第2实施方式的特征量的计算处理的流程图。
图7是示出第2实施方式的终端数量估计装置的结构的图。
图8是示出第2实施方式的终端数量估计处理的流程图。
图9是用于说明终端数量估计的第2思路的图。
图10是示出第3实施方式的特征量的计算处理的流程图。
图11是示出变量s、e的调整(適正化)处理的流程图。
图12是示出第4实施方式的终端数量估计装置的结构的图。
图13是示出第4实施方式的人口估计处理的流程图。
图14是示出第5实施方式的终端数量估计装置的结构的图。
图15是示出第5实施方式的人口估计处理的流程图。
图16是示出人口估计处理的输出示例的图。
图17是示出第6实施方式的终端数量估计装置的结构的图。
图18是用于说明网格与区域图的合成的图。
图19是用于说明各个分割区域的面积及面积比的计算的图。
图20是用于说明某个网格内的分割区域的人口的总和计算的图。
图21是用于说明第7实施方式的估计值的变换处理的图。
图22是示出用于变换为作为对象的输出单位的估计人口的矩阵的图。
图23是示出第8实施方式的终端数量估计装置的结构的图。
图24是用于说明非识别化处理的一例的图。
图25是示出第9实施方式的终端数量估计装置的结构的图。
图26是示出隐匿处理的一例的流程图。
图27是示出作为输出方式的示例的人口分布、人口变动及人口构成的图。
图28是用于说明终端数量估计的第2思路的图。
具体实施方式
参照附图说明本发明的实施方式。在尽可能的情况下,对相同的部分标注相同的标号,并省略重复说明。
[第1实施方式]
[通信***的结构]
图1是本实施方式的通信***1的***结构图。如图1所示,该通信***1构成为包括便携式终端100、BTS(基站)200、RNC(无线控制装置)300、交换机400、各种处理节点700及管理中心500。并且,该管理中心500构成为包括社会感测单元501、PETA挖掘单元(peta mining unit)502、流动人口统计(demography)单元503及可视化处理单元504。
交换机400通过BTS 200、RNC 300收集有关便携式终端100的后述的位置信息。RNC 300在与便携式终端100之间进行通信连接时,能够使用RRC连接请求信号中的延迟值来测定便携式终端100的位置。交换机400能够在便携式终端100执行通信连接时接收这样测定到的便携式终端100的位置信息。交换机400存储所接收到的位置信息,并按照预定的定时或者来自管理中心500的请求,将收集到的位置信息输出给管理中心500。
各种处理节点700通过RNC 300和交换机400取得便携式终端100的位置信息,根据情况还进行位置的再计算等,并按照预定的定时或者来自管理中心500的请求,将收集到的位置信息输出给管理中心500。
关于本实施方式中的便携式终端100的位置信息,能够采用表示利用位置登记信号而求出的归属扇区的扇区号、由GPS测位***或在采用PRACHPD的位置信息取得***中得到的位置测位数据等。便携式终端100的位置数据包括:上述那样的位置信息、以及识别该便携式终端的识别信息(例如线路号码等与便携式终端对应的信息)、和取得位置信息的位置信息取得时刻信息。另外,在使用线路号码作为识别信息的情况下,不是直接采用线路号码,而是优选使用与线路号码对应的值(例如线路号码的散列值(Hash值)等)(即,将线路号码非识别化)。在使用这样与线路号码对应的值(例如线路号码的散列值(Hash值)等)来进行用户的每种属性的处理的情况下,即使是包含于属性信息中的能够特定用户的信息,也需要采用与该能够特定用户的信息对应的值(即,将该能够特定用户的信息非识别化)。关于这种非识别化,将在第8实施方式中详细说明。
管理中心500如上所述构成为包括社会感测单元501、PETA挖掘单元502、流动人口统计单元503及可视化处理单元504,在各个单元中进行在便携式终端100的位置信息中使用的统计处理。另外,后述的终端数量估计装置10(图2)例如能够由管理中心500构成。
社会感测单元501是服务器装置,从各个交换机400及各种处理节点700或者离线收集包含便携式终端100的位置信息等的数据。该社会感测单元501构成为接收从各个交换机400及各种处理节点700定期输出的数据,或者按照在社会感测单元501中预先设定的定时从交换机400及各种处理节点700取得数据。
PETA挖掘单元502是服务器装置,将从社会感测单元501接收到的数据变换为预定的数据形式。例如,PETA挖掘单元502以用户ID为密钥进行分类处理,或按照每个区域进行分类处理。
流动人口统计单元503是服务器装置,进行针对在PETA挖掘单元502中处理的数据的累计处理即各个项目的统计处理。例如,流动人口统计单元503计数处于某个区域的用户数量,并且累计在服务区分布。
可视化处理单元504是服务器装置,能够可视地处理在流动人口统计单元503中被累计处理的数据。例如,可视化处理单元504能够在地图上对所累计的数据进行映射处理。由该可视化处理单元504进行处理后的数据被提供给企业、政府机构或者个人等,用于店铺开发、道路交通调查、灾害对策、环境对策等。另外,这样进行统计处理后的信息当然是以不特定个人等的方式进行加工的,以便不侵害隐私。
另外,社会感测单元501、PETA挖掘单元502、流动人口统计单元503及可视化处理单元504如前面所述均由服务器装置构成,虽然省略了图示,但是当然具有通常的信息处理装置的基本结构(即CPU、RAM、ROM、键盘和鼠标等输入装置、与外部进行通信的通信装置、存储信息的存储装置、以及显示器或打印机等输出装置)。
[终端数量估计装置的结构]
下面,说明本实施方式的终端数量估计装置。图2示出终端数量估计装置10的功能结构。如该图2所示,终端数量估计装置10具有位置数据取得部11(位置数据取得单元)、蓄积部12、观测期间取得部13(观测期间取得单元)、观测区域取得部14(观测区域取得单元)、观测对象取得部15(观测对象取得单元)、前后位置数据取得部16(前后位置数据取得单元)、特征量计算部17(特征量计算单元)、终端数量估计部18(终端数量估计单元)、及终端数量输出部19(输出单元)。
下面,说明图2所示的终端数量估计装置10的各部分的功能。位置数据取得部11从外部取得上述的位置数据并保存在蓄积部12中。蓄积部12保存有关多个用户(便携式终端)的多个时刻的位置数据。观测期间取得部13取得包含观测开始时刻与观测结束时刻的组的观测期间信息。观测区域取得部14取得与一个或多个位置信息相对应的观测区域信息。此处的观测区域信息被赋予为例如扇区号、经度纬度、地理范围(例如市区村镇)等,优选观测区域取得部14具有管理对应关联了取得的观测区域信息的表述形式与位置信息的表述形式的信息的管理数据库(例如扇区号与经度纬度的对应关系信息等)。
观测对象取得部15从蓄积部12取得包括有关待观测的观测期间的观测开始时刻以后而且是观测结束时刻以前的位置取得时刻信息、而且包括与有关待观测的观测区域的观测区域信息相对应的位置信息在内的一个或多个位置数据,作为观测对象位置数据。另外,观测对象位置数据也可以利用另外赋予的条件(例如便携式终端的用户的年龄段等)再进行锁定。
前后位置数据取得部16针对求出特征量的对象的位置数据(下面称为第1位置数据),取得包括与该第1位置数据相同的识别信息的位置数据中、该第1位置数据的紧前面的位置数据(下面称为第2位置数据)的位置取得时刻信息、以及该第1位置数据的紧后面的位置数据(下面称为第3位置数据)的位置取得时刻信息。另外,前后位置数据取得部16不是必须全部取得第2或者第3位置数据,只要至少取得位置数据中所包含的位置取得时刻信息即可。
在第1实施方式中,前后位置数据取得部16将由观测对象取得部15取得的观测对象位置数据作为第1位置数据,取得与该第1位置数据相关的第2、第3位置数据的位置取得时刻信息,后述的特征量计算部17计算有关观测对象位置数据的特征量。即,第1实施方式是将计算特征量的对象的位置数据锁定为观测对象位置数据的实施方式。与此相对,在后述的第2实施方式中,说明将所取得的所有位置数据作为对象来计算特征量的示例。
特征量计算部17计算有关各个第1位置数据(在第1实施方式中指观测对象位置数据)的特征量。例如,特征量计算部17计算第2位置数据的位置取得时刻与第3位置数据的位置取得时刻之差,作为有关该第1位置数据的特征量。并且,在第2位置数据的位置取得时刻是异常值的情况下,在此作为一例是在第1位置数据的位置取得时刻与第2位置数据的位置取得时刻之差大于预定的基准值(例如1小时)的情况下,特征量计算部17将从第1位置数据的位置取得时刻起向过去追溯预定的时间(例如1小时)得到的时刻,用作第2位置数据的位置取得时刻,来计算有关第1位置数据的特征量。同样,在第3位置数据的位置取得时刻是异常值的情况下,在此作为一例是在第1位置数据的位置取得时刻与第3位置数据的位置取得时刻之差大于预定的基准值(例如1小时)的情况下,特征量计算部17将从第1位置数据的位置取得时刻起向未来前进预定的时间(例如1小时)得到的时刻,用作第3位置数据的位置取得时刻,来计算有关第1位置数据的特征量。这种在第2、第3位置数据的位置取得时刻是异常值时的处理不是必须的处理,但通过进行上述处理,在由于便携式终端100位于服务区外或便携式终端100的电源被切断等,位置数据的取得时间间隔变得异常地长时,能够防止由于该异常长的取得时间间隔而造成的影响过大。
终端数量估计部18根据有关观测对象位置数据的特征量以及观测开始时刻与观测结束时刻之差即观测期间长度,估计在观测期间中位于观测区域的终端数量。详细情况将在后面叙述,终端数量估计部18将有关观测对象位置数据的特征量的总和除以观测期间长度的2倍,将相除得到的数值估计为终端数量。
终端数量输出部19输出通过估计而得到的终端数量。此处的输出包括显示输出、语音输出、打印输出等各种输出形式。
[终端数量估计的思路及计算方法]
下面,说明终端数量估计的思路及计算方法。如图3所示的模型那样,假设在某个观测期间(长度T)中,n个终端a1、a2、…、an通过扇区S,各个终端ai在观测期间内在扇区S的滞留时间为ti(0<ti≦T)。此时,位于扇区S中的终端数量m(实际是位于扇区S中的终端数量m的观测期间内的平均值)用下式(1)表示。
[数式1]
即,将各个终端ai在观测期间内在扇区S的滞留时间ti的总和除以观测期间的长度T,将相除的结果估计为终端数量m。其中,虽然终端ai在观测期间内在扇区S的滞留时间ti的真正值是不能观测到的,但是各个终端ai发出信号(例如位置登记信号),这些信号是能够观测到的。
将终端ai在观测期间内在扇区S发出的信号按照时刻顺序设为
[数式2]
(xi表示终端ai在观测期间内在扇区S发出的信号的总数),终端数量的估计一定是根据观测到的信号qij(j为1以上xi以下的整数)估计m的值。
另外,根据图4说明终端数量估计的计算方法。将从终端ai发送信号qij的密度(即,每单位时间的信号数)设为pi。此时,如果发送信号的概率相对于扇区是独立的,则终端ai在观测期间内在扇区S发出的信号的总数xi的期望值E(xi)是E(xi)=ti×pi,因而关于终端ai在观测期间内在扇区S的滞留时间ti的期望值E(xi),下面的式(2)成立。
E(ti)=xi/pi (2)
其中,在设信号qij的发送时刻为uij时,信号qij的密度pij通过下式(3)来赋值。
pij=2/(ui(j+1)-ui(j-1)) (3)
其中,在设信号qij为有关第1位置数据的信号时,信号qi(j-1)相当于有关第2位置数据的信号,信号qi(j+1)相当于有关第3位置数据的信号。在本实施方式中,将有关第2位置数据的信号qi(j-1)的发送时刻ui(j-1)与有关第3位置数据的信号qi(j+1)的发送时刻ui(j+1)之差即上式(3)中的(ui(j+1)-ui(j-1)),作为有关第1位置数据的特征量wij。因此,上式(3)表示如下。即,能够与密度pij的倒数相对应地计算出特征量wij。
pij=2/(ui(j+1)-ui(j-1))=2/wij (4)
此时,密度pi用下式赋值,
[数式3]
因而能够用下式(6)计算出终端数量m的估计值E(m)。
[数式4]
如图4的示例所示,在观测期间内而且是终端ai滞留于扇区S的期间内,假设终端ai发送信号qi1、qi2、qi3,在信号qi1的紧前面发送信号qi0,在信号qi3的紧后面发送信号qi4,将信号qi0、qi1、qi2、qi3、qi4的发送时刻分别设为ui0、ui1、ui2、ui3、ui4,上述思路相当于将终端ai在观测期间内在扇区S的滞留时间ti估计为从(ui0与ui1的中点)到(ui3与ui4的中点)的期间。另外,虽然终端ai不在观测期间内,但是在滞留于扇区S时发送信号qi4。另外,为了保持滞留时间ti的估计量的不偏倚性,在此作为一例,说明进行不将滞留时间ti的结束时刻估计为与观测期间T的结束时刻相同的处理。
[终端数量估计处理]
下面,说明本发明的终端数量估计方法的终端数量估计处理。在此,作为一例,假设便携式终端的位置数据中所包含的位置信息被赋予了该便携式终端所处扇区的扇区号。
如图5所示,首先位置数据取得部11从外部取得位置数据并保存在蓄积部12中(图5的步骤S1)。由此,蓄积部12保存有关多个用户(便携式终端)多个时刻的位置数据。另外,在执行步骤S1的处理之后,自步骤S2起的处理也可以间隔时间来执行。即,也可以执行步骤S1作为自步骤S2起的处理的事前准备。
然后,观测期间取得部13取得包括观测开始时刻和观测结束时刻的组的观测期间信息,并且观测区域取得部14取得与一个或多个位置信息相对应的观测区域信息(步骤S2)。在此,取得观测开始时刻T1和观测结束时刻T2的组作为观测期间信息,取得扇区号S作为观测区域信息。
然后,观测对象取得部15从蓄积部12取得包含观测开始时刻T1以后而且是观测结束时刻T2以前的位置取得时刻信息、而且包含与观测区域信息即扇区号S相对应的位置信息在内(例如位置信息是扇区号S)的一个或多个位置数据,作为观测对象位置数据(步骤S3)。即,观测对象取得部15取得符合下面的条件的位置数据作为观测对象位置数据。
条件1:位置取得时刻在观测开始时刻T1以后而且在观测结束时刻T2以前。即包含在观测期间内。
条件2:位置信息是扇区S。
然后,针对所取得的各个观测对象位置数据执行下面的步骤S4、S5的处理。在步骤S4,前后位置数据取得部16对于观测对象位置数据中作为求出特征量的对象的位置数据(第1位置数据),取得包括与该第1位置数据相同的识别信息的位置数据中、从位置取得时刻来看是该第1位置数据的紧前面的位置数据(第2位置数据)的位置取得时刻信息以及是该第1位置数据的紧后面的位置数据(第3位置数据)的位置取得时刻信息。另外,前后位置数据取得部16不是必须取得全部第2、第3位置数据,只要取得第2、第3位置数据中所包含的位置取得时刻信息即可。
另外,在步骤S5,特征量计算部17计算有关第1位置数据的特征量。关于其处理内容将使用图6进行说明。在此,将第1、第2、第3位置数据的位置取得时刻分别设为t1、t2、t3。并且,将作为用于判定第2位置数据的位置取得时刻t2是异常值的基准的预定的基准值(与第1、第2位置数据的位置取得时刻之差相关的基准值)设为基准值A(例如1小时),将作为用于判定第3位置数据的位置取得时刻t3是异常值的基准的预定的基准值(与第1、第3位置数据的位置取得时刻之差相关的基准值)设为基准值B(例如1小时)。
特征量计算部17计算第1、第2位置数据的位置取得时刻之差(即时刻t1与t2之差)Da、以及第1、第3位置数据的位置取得时刻之差(即时刻t1与t3之差)Db(图6的步骤S11)。并且,特征量计算部17判定第1、第2位置数据的位置取得时刻之差Da是否大于预定的基准值A(例如1小时)(步骤S12),如果差Da大于基准值A,则将从第1位置数据的位置取得时刻t1起向过去追溯预定的时间(例如1小时)得到的时刻,作为第2位置数据的位置取得时刻t2(步骤S13)。然后,特征量计算部17判定第1、第3位置数据的位置取得时刻之差Db是否大于预定的基准值B(例如1小时)(步骤S14),如果差Db大于基准值B,则将从第1位置数据的位置取得时刻t1起向未来前进预定的时间(例如1小时)得到的时刻,作为第3位置数据的位置取得时刻t3(步骤S15)。并且,特征量计算部17计算第2位置数据的位置取得时刻t2与第3位置数据的位置取得时刻t3之差,作为有关第1位置数据的特征量(步骤S16)。通过以上处理,结束有关某一个观测对象位置数据(第1位置数据)的步骤S4、S5的处理。
以后,针对各个观测对象位置数据执行上述的步骤S4、S5的处理,在对所有观测对象位置数据执行了处理之后(步骤S6为肯定判定),进入到步骤S7。
在步骤S7,终端数量估计部18按照前述的式(6)所示,将有关观测对象位置数据的特征量wij的总和除以观测期间长度T的2倍,将相除得到的数值估计为终端数量。另外,如根据式(6)可知,也可以是,终端数量估计部18将有关观测对象位置数据的特征量wij分别除以2,求出(特征量wij/2)的总和,将所求出的总和除以观测期间长度T,将相除得到的数值估计为终端数量。但是,如本实施方式这样将有关观测对象位置数据的特征量wij的总和除以观测期间长度T的2倍的计算方法,能够以极少的除法运算的次数完成处理,因而具有能够减轻处理负荷的优点。
另外,终端数量输出部19输出通过估计而得到的终端数量(步骤S8)。
根据以上的第1实施方式,在使用位置数据估计终端数量时,通过实施使用了前后的位置数据的取得时刻信息的校正,能够对接收间隔的变动所产生的影响进行校正,并且高精度地估计终端数量。
另外,在特征量的计算处理中,通过进行前述的第2、第3位置数据的位置取得时刻是异常值时的处理,在由于便携式终端100位于服务区外或便携式终端100的电源切断等,位置数据的取得时间间隔异常地长时,能够防止由于该异常长的取得时间间隔而造成的影响过大。
[第2实施方式]
在上述的第1实施方式中说明了将用于计算特征量的对象的位置数据锁定为观测对象位置数据的示例,而在下面的第2实施方式中说明将所取得的全部位置数据作为对象来计算特征量的示例,即在实施终端数量估计之前预先计算有关所有位置数据的特征量,使用其中有关观测对象位置数据的特征量来估计终端数量的示例。另外,第2实施方式的通信***的***结构与图1所示的第1实施方式的***结构相同,因而省略其***结构的说明。
如图7所示,第2实施方式的终端数量估计装置10具有与第1实施方式的终端数量估计装置(图2)相同的构成要素,各个构成要素的功能基本相同,因而在此以与第1实施方式的终端数量估计装置的不同之处为中心进行说明。
第2实施方式的观测对象取得部15取得包含有关待观测的观测期间的观测开始时刻以后而且是观测结束时刻以前的位置取得时刻信息、而且包含与有关待观测的观测区域的观测区域信息相对应的位置信息在内的一个或多个位置数据,作为观测对象位置数据,然后向终端数量估计部18输出观测对象位置数据。
前后位置数据取得部16将由位置数据取得部11取得的所有位置数据分别作为第1位置数据,并取得有关该第1位置数据的第2位置数据(紧前面的位置数据)以及第3位置数据(紧后面的位置数据)各自的位置取得时刻信息。另外,上述由位置数据取得部11取得的位置数据可以是在由位置数据取得部11取得后被蓄积在蓄积部12中的位置数据,也可以是不蓄积在蓄积部12中而从位置数据取得部11发送给前后位置数据取得部16的位置数据。
特征量计算单元17将由位置数据取得部11取得的所有位置数据分别作为第1位置数据,计算有关该第1位置数据的特征量。由于其计算结果是庞大的量,因而如图7所示,特征量计算部17具有用于保存计算结果的特征量的特征量保存部17A,优选由特征量保存部17A保存计算结果的特征量。另外,特征量计算部17计算第2、第3位置数据的位置取得时刻之差作为有关第1位置数据的特征量、以及进行如图6所示的第2或者第3位置数据的位置取得时刻是异常值时的处理,这些与第1实施方式相同。
终端数量估计部18从预先计算并保存在特征量保存部17A中的有关所有位置数据的特征量中,提取从观测对象取得部15接收到的有关观测对象位置数据的特征量,并根据有关该观测对象位置数据的特征量以及观测开始时刻与观测结束时刻之差(观测期间长度),估计在观测期间中位于观测区域的终端数量。具体地讲,与第1实施方式相同地,终端数量估计部18将有关观测对象位置数据的特征量的总和除以观测期间长度的2倍,将相除得到的数值估计为终端数量。
下面,说明第2实施方式的终端数量估计处理。在此,假设便携式终端的位置数据中所包含的位置信息被赋予了该便携式终端所在扇区的扇区号。
如图8所示,首先位置数据取得部11从外部取得位置数据并保存在蓄积部12中(图8的步骤S21)。但是,位置数据向蓄积部12的保存不是必须的,也可以直接从位置数据取得部11向前后位置数据取得部16转发位置数据,并执行后述的步骤S22的处理。并且,自步骤S22起的处理也可以在执行步骤S21的处理之后隔开时间间隔来执行。
然后,针对所取得的所有位置数据分别执行下面的步骤S22~S24的处理。在步骤S22,前后位置数据取得部16对于作为求出特征量的对象的一个位置数据(第1位置数据),取得包括与该第1位置数据相同的识别信息的位置数据中、从位置取得时刻来看是该第1位置数据的紧前面的位置数据(第2位置数据)的位置取得时刻信息以及是该第1位置数据的紧后面的位置数据(第3位置数据)的位置取得时刻信息。另外,前后位置数据取得部16不是必须取得全部第2、第3位置数据,只要取得第2、第3位置数据中所包含的位置取得时刻信息即可。另外,在步骤S23,特征量计算部17按照与第1实施方式相同的图6所示的步骤,计算有关第1位置数据的特征量。步骤S23的处理与前述第1实施方式的步骤S5的处理相同,因而省略说明。然后,在步骤S23得到的特征量被保存在特征量保存部17A中(步骤S24)。
通过以上处理,完成有关某一个观测对象位置数据(第1位置数据)的步骤S22~S24的处理。
以后,针对所有位置数据分别执行步骤S22~S24的处理。并且,在对所有位置数据完成步骤S22~S24的处理后(步骤S25为肯定判定),计算有关所有位置数据的特征量并保存在特征量保存部17A中。这样,能够在实施终端数量估计之前预先计算有关所有位置数据的特征量并进行保存。
在后面的步骤S26,观测期间取得部13取得包括观测开始时刻和观测结束时刻的组的观测期间信息,并且观测区域取得部14取得与一个或多个位置信息相对应的观测区域信息。在此,取得观测开始时刻T1和观测结束时刻T2的组作为观测期间信息,取得扇区号S作为观测区域信息。
然后,观测对象取得部15从蓄积部12取得包括观测开始时刻T1以后而且是观测结束时刻T2以前的位置取得时刻信息、而且包括与观测区域信息即扇区号S相对应的位置信息在内(例如位置信息是扇区号S)的一个或多个位置数据,作为观测对象位置数据(步骤S27)。即,观测对象取得部15取得符合下面的条件的位置数据作为观测对象位置数据。
条件1:位置取得时刻在观测开始时刻T1以后而且在观测结束时刻T2以前。即包含在观测期间内。
条件2:位置信息是扇区S。
然后,观测对象位置数据的取得也可以按照下面所述来进行。即,特征量计算部17将通过计算得到的特征量与计算对象的位置数据(上述第1位置数据)相关联,并将关联起来的带特征量的位置数据保存在特征量保存部17A中。另外,也可以是,终端数量估计部18经由观测对象取得部15取得观测期间信息和观测区域信息,从特征量保存部17读出符合观测期间及观测区域的条件的带特征量的位置数据,并作为观测对象位置数据进行取得。
返回到图8,在后面的步骤S28中,终端数量估计部18按照前述的式(6)所示,将有关观测对象位置数据的特征量wij的总和除以观测期间长度T的2倍,将相除得到的数值估计为终端数量。另外,如根据式(6)可知,也可以是,终端数量估计部18将有关观测对象位置数据的特征量wij分别除以2,求出(特征量wij/2)的总和,将所求出的总和除以观测期间长度T,将相除得到的数值估计为终端数量。但是,如本实施方式这样将有关观测对象位置数据的特征量wij的总和除以观测期间长度T的2倍的计算方法,能够以极少的除法运算的次数完成处理,因而具有能够减轻处理负荷的优点。
另外,终端数量输出部19输出通过估计而得到的终端数量(步骤S29)。
根据以上的第2实施方式,与第1实施方式相同地,在使用位置数据估计终端数量时,通过实施使用了前后的位置数据的取得时刻信息的校正,能够对接收间隔的变动所产生的影响进行校正,并且高精度地估计终端数量。
尤其是在第2实施方式中,在实施终端数量估计之前预先计算有关所有位置数据的特征量并进行保存,因而具有如下优点:能够缩短从终端数量估计装置取得观测期间信息及观测区域信息并开始终端数量的估计处理、到得到估计结果的终端数量的时间。
另外,在图8的处理中,不是必须在步骤S25之后执行步骤S26~S27的处理,也可以是同时并行地执行步骤S22~S25的处理和步骤S26~S27的处理。
[第3实施方式]
在第3实施方式中说明有关终端数量估计的第2方法以及基于该方法的特征量的计算处理。第3实施方式的通信***及终端数量估计装置的结构与第1、第2实施方式相同,因而省略说明。
图9是表示有关终端数量估计的第2思路的图。在该图9中,qij表示由终端ai生成的位置数据,其中qi1、qi2、qi3表示在观测期间内由处于扇区S的终端ai生成的位置数据。其中,假设qi1包含由于终端ai跨越位置登记区域(Location Area)边界而生成的位置登记信息(下面称为“LA跨越位置登记信息”),下面称为“LA跨越位置登记信息qi1”。在这种情况下下述的思路成立,即由于能够判定终端ai在生成了LA跨越位置登记信息qi1的定时进入到扇区S,因而将有关LA跨越位置登记信息qi1的特征量wi1作为LA跨越位置登记信息qi1的生成时刻与紧后面的位置数据qi2的生成时刻之差分,而不是如前述第1、第2实施方式那样的紧前面的位置数据qi0的生成时刻与紧后面的位置数据qi2的生成时刻之差分。
基于这种思路,终端ai在观测期间内处于扇区S的在服务区期间ti是图9中用粗线示出的期间,相比第1、第2实施方式的在服务区期间(图9中用虚线示出的期间)缩短了(LA跨越位置登记信息qi1的生成时刻与紧前面的位置数据qi0的生成时刻之差分的二分之一)。
并且,如图28所示,假设位置数据qi4包含LA跨越位置登记信息,下面称为“LA跨越位置登记信息qi4”。在这种情况下,能够判定终端ai在生成了LA跨越位置登记信息qi4的定时从扇区S退出。因此,在计算有关LA跨越位置登记信息qi4的紧前面的位置数据qi3的特征量wi3时,有关位置数据qi3的特征量wi3延长了(位置数据qi3的生成时刻与紧后面的LA跨越位置登记信息qi4的生成时刻之差分的二分之一)。即,终端ai在观测期间内位于扇区S服务区的在服务区期间ti是图28中用粗线示出的期间,相比第1、第2实施方式的在服务区期间(图28中用虚线示出的期间)延长了(位置数据qi3的生成时刻与紧后面的LA跨越位置登记信息qi4的生成时刻之差分的二分之一)。
关于如上所述的基于终端数量估计的第2思路的特征量的计算处理,使用图10进行说明。下面,将计算特征量的对象的位置数据称为“计算对象位置数据”。
如图10所示,特征量计算部17首先根据例如计算对象位置数据中所包含的呼叫类别信息,判定计算对象位置数据是否包含LA跨越位置登记信息(步骤S31)。其中,在计算对象位置数据包含LA跨越位置登记信息的情况下,将计算对象位置数据的位置取得时刻设定为用于计算特征量的第1变量s(下面称为“变量s”)(步骤S32),在计算对象位置数据不包含LA跨越位置登记信息的情况下,将计算对象位置数据的位置取得时刻与紧前面的位置数据的位置取得时刻的中点时刻设定为变量s(步骤S33)。
然后,特征量计算部17根据例如紧后面的位置数据中所包含的呼叫类别信息,判定紧后面的位置数据是否包括LA跨越位置登记信息(步骤S34)。其中,在紧后面的位置数据包括LA跨越位置登记信息的情况下,设定紧后面的位置数据的位置取得时刻为用于计算特征量的第2变量e(下面称为“变量e”)(步骤S35),在紧后面的位置数据不包括LA跨越位置登记信息的情况下,设定计算对象位置数据的位置取得时刻与紧后面的位置数据的位置取得时刻的中点时刻为变量e(步骤S36)。另外,不是必须根据呼叫类别信息进行上述的步骤S31、S34的判定处理,也可以根据其它信息来进行。例如,也可以预先保存表示位置登记区域的范围的区域信息,根据计算对象位置数据及紧后面的位置数据的位置信息和区域信息来进行判定处理。
然后,特征量计算部17进行图11所示的变量s、e的调整处理(步骤S37)。其中,将计算对象位置数据的位置取得时刻设为t1,将作为用于判定变量s是异常值的基准的预定的基准值设为基准值C(例如0.5小时),将作为用于判定变量e是异常值的基准的预定的基准值设为基准值D(例如0.5小时)。
特征量计算部17计算变量s与时刻t1之差Dc以及变量e与时刻t1之差Dd(图11的步骤S41)。并且,特征量计算部17判定变量s与时刻t1之差Dc是否大于预定的基准值C(例如0.5小时)(步骤S42),如果差Dc大于基准值C,则对变量s设定从时刻t1起向过去追溯预定的时间(例如0.5小时)得到的时刻(步骤S43)。然后,特征量计算部17判定变量e与时刻t1之差Dd是否大于预定的基准值D(例如0.5小时)(步骤S44),如果差Dd大于基准值D,则对变量e设定从时刻t1起向未来前进预定的时间(例如0.5小时)得到的时刻(步骤S45)。这样通过进行变量s、e的调整处理,在位置数据的取得时间间隔由于便携式终端100位于服务区外或便携式终端100的电源被切断等而变得异常地长时,能够防止由于该异常长的取得时间间隔而造成的影响过大。
然后返回到图10,特征量计算部17计算将(变量e-变量s)乘以2倍的值,并作为有关计算对象位置数据的特征量(步骤S38)。通过以上处理,能够得到有关计算对象位置数据的特征量。
根据以上所述的第3实施方式,在计算对象位置数据及其紧后面的位置数据中至少一方是LA跨越位置登记信息的情况下,通过考虑在生成该LA跨越位置登记信息的定时发生了向扇区S的流入或者从扇区S流出,能够得到高精度的特征量。
另外,在第3实施方式中说明的特征量的计算方法,也能够应用于如第1实施方式那样锁定为观测对象位置数据来计算特征量并估计终端数量的情况,以及如第2实施方式那样预先计算有关所有位置数据的特征量,并使用其中有关观测对象位置数据的特征量来估计终端数量的情况。
[第4实施方式]
在下面的第4、第5实施方式中,对使用用于将通过估计得到的终端数量变换为人口的系数(下面称为“放大系数”)来估计人口的示例进行说明。其中,在第4实施方式中说明将特征量乘以放大系数,并使用该相乘结果的累计值求出人口的实施方式,在第5实施方式中说明按照每人口估计单位(例如每种属性或每个时间段)来累计特征量,将对应于人口估计单位的放大系数与该累计值相乘,使用该相乘结果求出人口的实施方式。
在自第4实施方式起的实施方式中,以预先计算有关所有位置数据的特征量,并使用其中有关观测对象位置数据的特征量来估计人口或终端数量的处理为基础进行说明。但是,也能够应用于锁定为观测对象位置数据来计算特征量并估计人口或终端数量的处理。
如图12所示,第4实施方式的终端数量估计装置10的功能块结构在以下方面与第2实施方式的终端数量估计装置的功能块结构(图7)不同,因而对该不同之处进行说明。终端数量估计装置10具有根据特征量估计人口的人口估计部21来取代终端数量估计部18,并且具有输出人口估计部21的估计值的人口输出部22来取代终端数量输出部19。但是,权利要求书中的“终端数量估计单元”对应于终端数量估计部18及人口估计部21,“输出单元”对应于终端数量输出部19及人口输出部22。
另外,终端数量估计装置10具有属性/放大系数存储部23(放大系数存储单元),其存储了每个便携式终端用户的属性信息及预先求出的每种属性的放大系数。特征量计算部17以位置数据的用户识别信息(例如已散列的电话号码)为密钥,从属性/放大系数存储部23读出用户的属性信息及与该属性信息相关的放大系数,使计算出的特征量和读出的属性信息及放大系数与位置数据相关联,将关联后的位置数据保存在特征量保存部17A中。
下面,使用图13说明第4实施方式的终端数量估计装置10的处理。图13所示的第4实施方式的处理在以下方面与第2实施方式的处理(图8)不同,因而对该不同之处进行说明。
如图13所示,当在步骤S23进行特征量的计算处理之后,在步骤S24A,特征量计算部17以位置数据的用户识别信息(例如已散列的电话号码)为密钥,从属性/放大系数存储部23读出用户的属性信息及与该属性信息相关的放大系数,使计算出的特征量和读出的属性信息及放大系数与位置数据相关联,将关联后的位置数据保存在特征量保存部17A中。通过对所有位置数据依次执行图13的步骤S22~S24A的处理,将关联了各自的特征量、属性信息及放大系数的位置数据保存在特征量保存部17A中。
然后,当在步骤S26取得观测期间信息和观测区域信息后,在步骤S27,人口估计部21经由观测对象取得部15接收观测期间信息和观测区域信息,并从特征量保存部17A取得符合观测期间信息和观测区域信息的条件的位置数据(即观测对象位置数据)。并且,在步骤S28A,人口估计部21对于各个观测对象位置数据将特征量和放大系数相乘,将得到的相乘结果的总和除以(观测长度×2),将相除得到的值估计为在观测期间中位于观测区域内的人口。并且,人口估计部21根据与各个观测对象位置数据相关联的属性信息和各个观测对象位置数据中所包含的取得时刻信息,将有关各个观测对象位置数据的上述相乘结果按照每种属性进行累计,将其累计值除以(观测长度×2),将相除得到的值估计为每种属性的人口。另外,在此叙述了按照每种属性估计人口的示例,但是估计人口的单位不限于属性,也可以采用场所或时间段等。下面,将这些估计单位(属性、场所、时间段等)称为“人口估计单位”。
另外,在后面的步骤S29,人口输出部22输出通过步骤S28A的估计而得到的、在观测期间中位于观测区域内的人口及按照每人口估计单位的人口。在此,例如输出如图16所示的针对观测区域估计出的人口(在图16中表述为“估计人口”)、或按照性别/年龄段/住所等每种属性的人口。并且,也可以采用将多种属性进行组合的条件(例如,将性别和住所进行组合的条件“居住在东京都的女性”)。另外,此处的输出包括显示输出、语音输出、打印输出等各种输出形式。
根据以上所述的第4实施方式,能够得到在观测期间中位于观测区域内的人口及按照每人口估计单位(例如属性或时间段)的人口。
另外,在上述的步骤S28A、S29,叙述了估计观测期间中的观测区域内的人口和按照每人口估计单位的人口两者并进行输出的示例,然而不是必须对这两者进行估计并输出,也可以只对一方进行估计并输出。
并且,在第4实施方式中说明了预先求出放大系数的示例,放大系数也可以按照下面所述来取得。作为一例,放大系数能够采用在服务区率与终端的普及率之积(即,在服务区数量与人口的比率)的倒数。其中,“在服务区率”是指在服务区数量与签约台数的比率,“普及率”是指签约台数与人口的比率。优选按照上述的每种放大系数计算单位来导出这种放大系数,但不是必须的。另外,也可以是,例如按照下面所述使用根据特征量及观测期间长度而估计出的终端数量(在服务区数量)来导出放大系数。利用在第1~第3实施方式中叙述的方法从位置数据求出特征量,根据特征量及观测期间长度来累计每种放大系数计算单位的终端数量,由此得到用户数量金字塔数据,并且取得预先求出的相同的放大系数计算单位中的人口金字塔数据作为统计数据(例如居民基本登记册等)。并且,在用户数量金字塔数据和人口金字塔数据中计算每种放大系数计算单位的位置数据的取得率(即在服务区数量/人口)。在此得到的“位置数据的取得率(即在服务区数量/人口)”相当于前述的“在服务区率与终端的普及率之积”。能够导出这样得到的“位置数据的取得率”的倒数作为放大系数。另外,关于计算放大系数的放大系数计算单位,作为一例,也可以采用住所的都道府县、间隔5岁或10岁的年龄段、男女单位、作为时间段的1小时等,还可以采用将其中两种以上单位进行组合的单位。例如,在将放大系数计算单位设为“居住在东京都的二十几岁的男性”的情况下,通过提取日本全国的与居住在东京都的(即,用户属性中的住所信息是东京都)二十几岁的男性对应的位置数据并累计终端数量,能够得到用户数量金字塔数据,并且从统计数据中取得与居住在东京都的二十几岁的男性相关的人口金字塔数据。另外,在得到上述用户数量金字塔数据时,关于条件“居住在东京都”,不是仅提取处于东京都服务区的用户的位置数据,而是提取用户属性中的住所信息是东京都的位置数据。并且,根据用户数量金字塔数据和人口金字塔数据计算放大系数计算单位(此处指居住在东京都的二十几岁的男性)的位置数据的取得率(即,在服务区数量/人口),并导出所得到的“位置数据的取得率”的倒数作为放大系数。另外,在本申请中,将放大系数计算单位和人口估计单位作为相同单位进行说明,但这毕竟仅是一例,并不限于此。
另外,在第4实施方式中说明了以第2实施方式为基础来求出人口的处理,然而第4实施方式也能够应用于前述的第1、第3实施方式。
[第5实施方式]
在第5实施方式中说明这样的实施方式,即按照每人口估计单位(例如每种属性或每个时间段)来累计特征量,将对应于人口估计单位的放大系数与该累计值相乘,由此求出人口。
如图14所示,第5实施方式的终端数量估计装置10的功能块结构与前述第4实施方式的终端数量估计装置(图12)的功能块结构相似,但是与第4实施方式的不同之处在于,关于放大系数,不是由特征量计算部17,而是由人口估计部21从属性/放大系数存储部23读出,并进行使用了放大系数的人口的估计。即,特征量计算部17及人口估计部21的处理与第4实施方式不同,因而使用图15对该不同之处进行说明。
如图15所示,当在步骤S23中进行特征量的计算处理之后,在步骤S24B,特征量计算部17以用户识别信息(例如已散列的电话号码)为密钥,根据位置数据确定该位置数据的用户的属性信息,从属性/放大系数存储部23读出该属性信息,使计算出的特征量和读出的属性信息与位置数据相关联,将关联后的位置数据保存在特征量保存部17A中。通过对所有位置数据依次执行图15的步骤S22~S24B的处理,将关联了各自的特征量及属性信息的位置数据保存在特征量保存部17A中。
然后,当在步骤S26取得观测期间信息和观测区域信息后,在步骤S27,人口估计部21从特征量保存部17A取得与观测期间信息和观测区域信息的条件相符的位置数据(即观测对象位置数据)。并且,在步骤S28B,人口估计部21按照下面所述进行使用了放大系数的人口的估计。
首先,人口估计部21确定所取得的观测对象位置数据中符合某种人口估计单位(属性或时间段)的条件的观测对象位置数据,对确定出的观测对象位置数据累计特征量(步骤S51),从属性/放大系数存储部23读出有关该人口估计单位的放大系数,将该放大系数和特征量累计值相乘,将相乘结果(放大系数×特征量累计值)除以(观测期间长度×2),将相除得到的值估计为该人口估计单位的人口(步骤S52)。并且,针对各种人口估计单位执行上述步骤S51~S52,在对所有人口估计单位执行完成后,进入到步骤S54。在该时刻,作为估计结果能够得到有关每人口估计单位的人口。另外,将有关每人口估计单位的人口进行总和,将该总和估计为在观测期间中位于观测区域内的人口(步骤S54)。
通过以上的步骤S28B,能够得到在观测期间中位于观测区域内的人口以及每人口估计单位的人口。并且,在后面的步骤S29,人口输出部22输出通过步骤S28A的估计而得到的、在观测期间中位于观测区域内的人口及按照每人口估计单位的人口。
根据以上所述的第5实施方式,例如能够得到如图16所示的针对观测区域而估计出的人口(在图16中表述为“估计人口”)、和按照属性或时间段等每人口估计单位的人口。
与第4实施方式相同地,在步骤S28A、S29叙述了估计观测期间中的观测区域内的人口和每人口估计单位的人口两者并进行输出的示例,然而不是必须对这两者进行估计并输出,也可以只对一方进行估计并输出。并且,在第5实施方式中说明了以第2实施方式为基础来求出人口的处理,然而第5实施方式也能够应用于前述的第1、第3实施方式。
[第6实施方式]
在下面的第6、第7实施方式中,对将作为累计单位的每个区域的估计值(人口或终端数量)变换为每种输出单位(在此作为一例是网格(mesh))的估计值的示例进行说明。其中,在第6实施方式中说明利用单个频带的室外站的通信区域所在的环境中的处理,在第7实施方式中说明室内站的通信区域以及利用电波到达范围不同的频带的多个室外站的通信区域中两个以上的通信区域在地理上重复地存在的环境中的处理。另外,在输出单位和累计单位相同的情况下,不需要在下面的第6、第7实施方式中叙述的变换处理。
在第6实施方式中说明对第4实施方式的估计人口的终端数量估计装置增加了向上述每种网格的估计值进行变换的功能的终端数量估计装置及其处理。
如图17所示,第6实施方式的终端数量估计装置10的功能块结构构成为在第4实施方式的终端数量估计装置(图12)中,在人口估计部21与人口输出部22之间增加了变换部24(变换单元)。变换部24通过后述的处理,将通过人口估计部21的估计而得到的每个区域的人口变换为每个网格的人口。
下面,根据图18具体说明变换部24的处理。图18(a)是表示区域的区域范围的图,图18(b)是表示网格的图,图18(c)是将区域和网格合成的合成图。
变换部24将根据预先存储的区域边界信息而再现的区域图(参照图18(a))、和根据预定的区划规则而再现的二维网格(参照图18(b))进行合成,得到如图18(c)所示的合成图。然后,变换部24在上述合成图中根据网格边界将各个区域进行分割。例如,如图19所示,图18(a)的区域A根据网格边界被分割为4个分割区域A-1、A-2、A-3、A-4。并且,变换部24计算各个分割区域的面积,并计算各个分割区域的面积比。例如,如图19所示,假设计算出分割区域A-1、A-2、A-3、A-4的面积分别是10m2、50m2、100m2、40m2,并计算出分割区域A-1、A-2、A-3、A-4的面积比(例如百分比)分别是5%、25%、50%、20%。
另外,变换部24不是必须计算各个分割区域的面积比。例如,也可以构成为,预先求出各个分割区域的面积比,变换部24能够从终端数量估计装置10内的未图示的表或者外部参照各个分割区域的面积比的信息。
然后,变换部24计算各个分割区域的人口。例如,假设图18(a)中的区域A的人口是800人,如图18所示计算出分割区域A-2的人口为200人(即800人×25%)。同样,假设区域B、C的人口分别是500人、750人,计算出区域B的面积比80%的分割区域B-1的人口是400人(即500人×80%),计算出区域C的面积比80%的分割区域C-4的人口是600人(即750人×80%)。
另外,变换部24计算一个网格中所包含的多个分割区域的人口的总和,由此计算该网格的人口。在图20的示例中,计算出一个网格中所包含的分割区域A-2、B-1、C-4的人口的总和是1200人(即200人+400人+600人),将该1200人作为该网格的人口。
如上所述,能够在利用单个频带的室外站的通信区域所在的环境中,将每累计单位的人口变换为每输出单位的人口。
另外,在第6实施方式中说明了以第4实施方式为基础,将每累计单位的人口变换为每输出单位的人口的处理,然而第6实施方式也能够适用于前述的第5实施方式。并且,在第6实施方式中叙述的变换处理也能够应用于将每累计单位的终端数量变换为每输出单位的终端数量的情况,通过应用于前述的第1~第3实施方式,能够将每累计单位的终端数量变换为每输出单位的终端数量。
[第7实施方式]
在第7实施方式中说明这样的示例,在室内站的通信区域以及利用电波到达范围不同的频带的多个室外站的通信区域中两个以上的通信区域在地理上重复存在的环境中,将作为累计单位的每个区域(扇区)的估计值(终端数量或者人口)变换为每输出单位(在此作为一例是指网格)的估计值。
第7实施方式的终端数量估计装置的功能块结构与第6实施方式相同,但是变换部24的处理不同,因而根据图21、图22来说明变换部24的处理。
如图21所示,在室内站的通信区域以及利用电波到达范围不同的频带(室外2GHz/1.7GHz和室外800MHz)的多个室外站的通信区域中两个以上的通信区域在地理上重复存在的环境中,变换部24针对各个通信区域进行在第6实施方式中叙述的变换处理,由此求出有关各个通信区域的每输出单位(网格)的人口,最后按照每输出单位将有关各个通信区域的人口进行合计,由此得到每输出单位的人口。
在图21的示例中,变换部24首先对利用室外2GHz/1.7GHz的室外站的通信区域、利用室外800MHz的室外站的通信区域、室内站的通信区域,分别进行在第6实施方式中叙述的变换处理。例如,在利用室外2GHz/1.7GHz的室外站的通信区域中,假设输出单位Q与区域A重合的分割区域占区域A整体的面积比是40%,通过将区域A的估计人口100人乘以面积比0.4,能够得到输出单位Q与区域A重合的分割区域的估计人口是40人。同样,能够得到输出单位Q与区域B重合的分割区域的估计人口是3人(区域B的估计人口30人×面积比0.1)、以及输出单位Q与区域C重合的分割区域的估计人口是5人(区域C的估计人口100人×面积比0.05)。同样,对于利用室外800MHz的室外站的通信区域,能够得到输出单位Q与区域D重合的分割区域的估计人口是3人(区域D的估计人口10人×面积比0.3)、以及输出单位Q与区域F重合的分割区域的估计人口是9人(区域F的估计人口30人×面积比0.3)。另一方面,关于室内站,各个室内站的电波到达范围的区域非常小,在图21的示例中,一个室内站的区域L整体与输出单位Q重合,因而能够认为面积比是100%。因此,通过向区域L的估计人口10乘以面积比1.0,能够得到输出单位Q与区域L重合的区域(在该示例中是区域L整体)的估计人口是10人。
最后,变换部24将如上所述得到的输出单位Q与各个区域重合的区域的估计人口进行合计,得到输出单位Q的估计人口是70人。通过以上所述的处理,能够从每累计单位的人口变换为输出单位Q的累计人口。
图21表示向一个输出单位Q的估计人口的变换,但是通过对其它的输出单位执行相同的处理,能够进行向作为对象的所有输出单位的估计人口的变换。
图22表示用于进行向作为对象的n个输出单位的估计人口的变换的矩阵。即,图22中的式子的右边的
[数式5]
(j表示1以上m以下的整数(m表示与作为对象的n个输出单位中任意一个输出单位重合的累计单位的数量))表示通过估计而求出的每累计单位的人口(估计人口),左边的
[数式6]
(i表示1以上n以下的整数)表示每输出单位的人口,右边的矩阵中的
[数式7]
表示用于从累计单位bj的人口变换为输出单位ai的人口的变换系数。其中的变换系数相当于前述的分割区域在原来的区域整体中所占的面积比。
图22中的各个变换系数能够根据估计单位(例如区域)与输出单位(例如网格)的位置关系而预先求出,通过预先求出各个变换系数并存储图22的式子,能够使用图22的式子容易且快速地从通过估计而求出的每累计单位的人口(估计人口)变换为每输出单位的人口。
通过以上所述的处理,即使是在室内站的通信区域以及利用电波到达范围不同的频带的多个室外站的通信区域中两个以上的通信区域在地理上重复存在的环境中,也能够将每累计单位的人口变换为每输出单位的人口。
另外,在第6、第7实施方式中说明了将预先通过估计而得到的每累计单位的人口变换为每输出单位的人口的示例,但也可以考虑以下的变形例。
将特征量和放大系数、以及与该位置数据所属的区域(累计单位)在地理上重合的多个输出单位分别相关的输出单位ID及与该输出单位的重合部分的面积比(即,重合部分在累计单位整体中所占的面积比)的一个以上的组合,与位置数据相关联地保存在特征量保存部17A中。
在此,也可以是,对于与相同的输出单位ID相关联的位置数据计算(特征量×面积比),将所得到的(特征量×面积比)按照每输出单位进行累计,将每输出单位的累计值除以(观测期间长度×2),由此估计每输出单位的终端数量。
并且,也可以是,对于与相同的输出单位ID相关联的位置数据计算(特征量×面积比×放大系数),将所得到的(特征量×面积比×放大系数)按照每输出单位进行累计,将每输出单位的累计值除以(观测期间长度×2),由此估计每输出单位的人口。此时,将每输出单位的(特征量×面积比×放大系数)按照每人口估计单位(例如属性或时间段等)进行累计,将所得到的每输出单位及每人口估计单位的累计值除以(观测期间长度×2),由此估计每输出单位及每人口估计单位的人口。
并且,也可以是,对于与相同的输出单位ID相关联的位置数据,按照每人口估计单位(例如属性或时间段等)计算(特征量×面积比),将所得到的每人口估计单位的(特征量×面积比)按照每输出单位进行累计,将所得到的每输出单位及每人口估计单位的(特征量×面积比)的累计值与有关该人口估计单位的放大系数相乘,将所得到的相乘结果除以(观测期间长度×2),由此估计每输出单位及每人口估计单位的人口。但是,在这种情况下,不是必须将放大系数与位置数据相关联,在上述的乘法运算时,有关各个人口估计单位的放大系数也可以使用例如与位置数据不同的预先存储的放大系数或者从外部取得的放大系数。通过上述的处理,如果能够得到每输出单位及每人口估计单位的人口,通过按照每输出单位将每输出单位及每人口估计单位的人口进行合计,也可以估计每输出单位的人口。
并且,如前面所述,不是必须按照人口估计单位来求出放大系数,例如也可以不使用全体共用的放大系数等每人口估计单位的放大系数。在这种情况下,作为一例,也可以是,将有关作为对象的所有位置数据的(特征量×面积比)进行累计,向所得到的累计值乘以共同的放大系数,将相乘结果除以(观测期间长度×2),由此估计人口。
另外,在上述的各种变形例中使用的特征量可以利用前述的第1~第3实施方式的任意一种方法计算出。
另外,在第7实施方式中说明了以第4实施方式为基础将每累计单位的人口变换为每输出单位的人口的处理,但是第7实施方式也能够适用于前述的第5实施方式。并且,在第7实施方式中叙述的变换处理也能够应用于将每累计单位的终端数量变换为每输出单位的终端数量的情况,通过应用于前述的第1~第3实施方式,能够将每累计单位的终端数量变换为每输出单位的终端数量。
[第8实施方式]
在第8实施方式中说明进行非识别化处理的实施方式,该非识别化处理用于从位置数据或属性信息中去除具有个人识别性质的信息。
如图23所示,第8实施方式的终端数量估计装置10的功能块结构构成为,在第4实施方式的终端数量估计装置的功能块结构(图12中),在位置数据取得部11与蓄积部12之间增加了进行非识别化处理的非识别化部25(非识别化单元)。
非识别化部25针对位置信息中所包含的识别信息(例如电话号码),进行包含如下变换的非识别化处理,该变换是向基于单向性函数的不可逆编码的变换。此处采用单向性函数是为了防止根据变换后的信息进行复原,关于单向性函数例如可以采用如图24所示的、由国内外的评价项目或评价机构推荐的基于散列函数(Hash Function)的带密钥散列函数。
并且,如图24所示,在进行采用便携式终端的用户的属性信息的处理的情况下,非识别化部25在进行该处理之前,针对属性信息中能够确定个人的号码(例如电话号码),进行包含向基于单向性函数的不可逆编码的变换的非识别化处理。另外,也可以是,非识别化部25删除属性信息中的姓名信息,将出生年月日信息置换为年龄信息,将住所信息置换为删除了门牌号信息的街道号级别的住所信息。
通过以上所述的非识别化部25的非识别化处理,能够从位置数据或属性信息中去除具有个人识别性质的信息,因而能够防止从位置数据或属性信息识别出个人。
另外,在第8实施方式中说明了以第4实施方式为基础的非识别化处理,但是第8实施方式也能够应用于前述的第5实施方式。并且,在第8实施方式中叙述的非识别化处理不仅能够应用于人口估计,而且也能够应用于在第1~第3实施方式中说明的终端数量估计。并且,在第8实施方式中叙述的非识别化处理也能够应用于进行在第6、第7实施方式中说明的变换处理的情况。
[第9实施方式]
在第9实施方式中说明这样的实施方式,即在输出某个估计值(终端数量或人口等)之前,根据预先设定的基准对估计值进行隐匿处理。
如图25所示,第9实施方式的终端数量估计装置10的功能块结构构成为,在第8实施方式的终端数量估计装置的功能块结构(图23中),在人口估计部21与人口输出部22之间增加了进行隐匿处理的隐匿处理部26(隐匿处理单元)。
隐匿处理部26在从人口估计部21接收到估计值(估计人口等)时,进行例如图26所示的隐匿处理。即,隐匿处理部26判定取得源终端数量是否小于用于判定为需要隐匿处理的预先设定的基准值(作为一例是10)(图26的步骤S61),所述取得源终端数量表示作为估计的基础的每个区域(小区)的位置数据是从哪台终端取得的。另外,取得源终端数量表示去除了有关同一终端的重复的唯一的终端数量。在进行此处的判定时需要每个区域(小区)的位置数据的取得源终端数量,但作为一例也可以是,人口估计部21在累计与位置数据相对应的特征量时,计数该位置数据中的识别信息(例如,前述的非识别化部25的非识别化处理之后的已散列的电话号码)的数量,将所得到的已散列的电话号码总数的信息转发给隐匿处理部26,隐匿处理部26将已散列的电话号码总数用作每个区域(小区)的位置数据的取得源终端数量。
在步骤S61,如果作为估计的基础的每个区域(小区)的位置数据的取得源终端数量小于基准值,则隐匿处理部26将有关该区域(小区)的估计值设为零,由此将该估计值隐匿化(步骤S62)。另外,此处的隐匿方法不限于将估计值设为零,也可以采用利用预定的文字或记号(例如“×”等)表示估计值的方法等其它方法。
另一方面,在步骤S61,如果作为估计的基础的每个区域(小区)的位置数据的取得源终端数量为基准值以上,隐匿处理部26针对有关该区域(小区)的估计值,进行如下所述的在输出估计值时使用的级幅的概率圆整(步骤S63)。即,将有关该区域(小区)的估计值设为x,将级幅设为k,在kn≦x<k(n+1)(n为整数)的情况下,隐匿处理部26用概率(x-kn)/k将该估计值x圆整为k(n+1),用概率(k(n+1)-x)/k将该估计值x圆整为kn。
例如,在估计值x为23、级幅k为10的情况下,由于k×2≦x<k(2+1),因而n=2,用概率0.3(30%的概率)将估计值[23]圆整为[30],用概率0.7(70%的概率)将估计值[23]圆整为[20]。
通过以上所述的隐匿处理部26进行的隐匿处理,能够防止根据估计结果来确定个人,提高估计结果的有用性。并且,能够防止利用其它的值估计出被隐匿化的值的问题。
另外,关于概率圆整中的级幅,也可以从在累计中使用的位置数据中提取放大系数最大的放大系数,将把该放大系数设为预定倍(例如10倍)后得到的值设定为级幅。并且,还可以预先确定级幅来实施概率圆整。此时,也可以采用将包括放大系数超过预先设定的级幅的预定比率(例如1/10)的位置数据在内的区域(小区)的数据舍去的处理规则。
并且,隐匿处理部26的隐匿处理不限于图26所示的处理,也可以采用其它处理。例如,在图26的步骤S61,也可以取代“作为估计的基础的每个区域(小区)的位置数据的取得源终端数量”,而判定“作为估计的基础的每个区域(小区)的位置数据的个数”是否小于基准值,还可以判定“所输出的估计值(人口或者终端数量)”是否小于基准值。并且,当在图26的步骤S61判定为否定的情况下,也可以省略步骤S63的处理。
另外,在第9实施方式中以第8实施方式(进行非识别化处理的方式)为基础说明了隐匿处理,但也能够应用于不进行如第8实施方式所示的非识别化处理的情况。并且,第9实施方式也能够适用于在第4、第5实施方式中叙述的人口估计,还能够应用于在第1~第3实施方式中叙述的终端数量估计。并且,第9实施方式也能够应用于进行在第6、第7实施方式中叙述的变换处理的情况。
另外,关于前述的各种实施方式的人口估计单位、以及与输出相关的输出单位及输出形式说明如下。
人口估计单位能够根据用户的属性(年龄、性别、住所等)、时间段、场所(小区或网格等)中至少一种要素进行设定。例如,在对于人口估计单位设定了某种属性的情况下,如在第4实施方式中叙述的那样,人口估计部21根据与各个观测对象位置数据相关联的属性信息,按照被设为人口估计单位的每种属性,将与各个观测对象位置数据相关联的特征量和放大系数相乘,将该相乘结果的总和除以(观测长度×2),由此能够估计出被设为人口估计单位的每种属性的人口。例如,如图16所示,能够输出针对观测区域估计出的人口或性别/年龄段/住所这样的每种属性的人口。并且,也可以采用将多种属性进行组合的条件(例如,将性别和住所进行组合得到的条件“居住在东京都的女性”)。并且,在作为人口估计单位设定了例如每天间隔1小时的时间段的情况下,人口估计部21根据各个观测对象位置数据中所包含的取得时刻信息来确定各个观测对象位置数据属于哪个时间段,按照间隔1小时的每个时间段,将与各个观测对象位置数据相关联的特征量和放大系数相乘,将该相乘结果的总和除以(观测长度×2),由此能够估计出每个时间段的人口。另外,在设定时间段作为人口估计单位的情况下,也可以设定按照时间序列而不连续的时间段(例如,每天13:00~14:00的时间段,每周星期六及星期日的10:00~11:00的时间段等)。
另一方面,与输出相关的输出单位也与估计单位相同地,能够根据用户的属性(年龄、性别、住所等)、时间段、场所(扇区或网格等)中至少一种要素进行设定。
并且,关于与输出相关的输出形式,可以选择例如图27(a)所示的表示人口分布的图、图27(b)所示的表示按时间序列的人口变动的图、图27(c)所示的表示人口构成的图等各种输出形式。并且,关于这些各种输出形式,也可以按照用户的属性(年龄、性别、住所等)、时间段、场所(扇区或网格等)中的一个或者多个要素的每种组合进行输出,例如也可以按照图27(c)所示输出间隔5岁的每个年龄段的男女分类的人口构成。
下面,叙述与特征量相关的变形例。在前述的第1、第2实施方式中示出了这样的示例,计算求出特征量的对象的位置数据(第1位置数据)前后的位置数据的时间差(第2位置数据与第3位置数据的时间差),作为第1位置数据的特征量。用式子表示该示例,则特征量可以用下面的式(7)表示。另外,下面的式(7)仅仅是将前述的式(4)进行变形,因而与式(4)等效(即,不是变更了式(4)的思路的方式)。
wij=ui(j+1)-ui(j-1) (7)
本变形例表示在特征量计算部17中计算出的特征量的计算方法的另一个变形。
在本变形例中,特征量计算部17在求出上述的第1位置数据的特征量时,考虑有关第2位置数据及第3位置数据的类别信息(例如后述的位置数据的生成原因(生成定时))。具体地讲,特征量计算部17计算向第3位置数据与第1位置数据的时间差乘以与第3位置数据的类别信息(此处是生成原因)对应的校正系数α得到的值,并且计算向第1位置数据与第2位置数据的时间差乘以与第2位置数据的类别信息(此处是生成原因)对应的校正系数β得到的值。但是,除上述情况之外,特征量计算部17也可以根据第1位置数据的类别信息确定校正系数α或β,并且也可以根据第1及第2位置数据的类别信息确定校正系数β,还可以根据第1及第3位置数据的类别信息确定校正系数α。并且,特征量计算部17将通过这些乘法运算得到的值进行合计,将合计值作为第1位置数据的特征量。如果用式子表示特征量计算部17中特征量的计算处理,则可以用下式(8)表示。
wij=α(ui(j+1)-uij)+β(uij-ui(j-1)) (8)
关于有关第2位置数据及第3位置数据的类别信息,例如在位置数据是位置登记信息的情况下,可以举出与该位置登记信息的生成原因相关的信息,与该生成原因相关的信息包含在所生成的位置登记信息中。关于位置登记信息的生成原因,可以举出终端跨越位置登记区域(Location Area)边界、根据周期性地进行的位置登记而生成、由于终端的电源接通等而执行连接(attach)处理、由于终端的电源切断等而执行断开(detach)处理等,与这些生成原因相对应地,预先确定校正系数α及β的设定值。并且,特征量计算部17可以根据与第3位置数据的生成原因相关的信息,设定有关第3位置数据的校正系数α,并根据与第2位置数据的生成原因相关的信息,设定有关第2位置数据的校正系数β。另外,也可以是,校正系数α及β均预先设定为0以上1以下的值。但是,不一定必须是该数值范围。
例如,在诸如基于周期性进行的位置登记的位置登记信息那样、终端的位置与位置登记信息的生成契机没有关系的位置登记信息的情况下,可以认为在当前的扇区中滞留的时间的期望值在生成该位置登记信息的前后是相同的。另一方面,在由于终端跨越位置登记区域边界而生成的位置登记信息的情况下,能够判定为至少在生成该位置登记信息之前终端没有滞留在当前的扇区中。因此,认为在生成该位置登记信息之前终端滞留在当前的扇区中的时间为0,如果第1位置数据的类别信息(生成原因)是“跨越位置登记区域边界”,则能够将上述式(8)中的校正系数β(即,有关与紧前面的位置数据的时间差的校正系数β)设定为0。由此,能够计算出更切合实际的特征量。另外,这样在第1位置数据的类别信息(生成原因)是“跨越位置登记区域边界”的情况下,如果将校正系数β设为0来计算特征量,则能够得到与前述第3实施方式相同的修改。
这样,特征量计算部17在计算有关对象的位置数据(第1位置数据)的特征量的情况下,根据有关第1位置数据的前后的位置数据即第2和第3位置数据的类别信息(作为一例是位置数据的生成原因),对第2位置数据与第3位置数据的时间差进行校正,并使用校正后的时间差计算特征量。由此,能够根据位置数据的类别信息更加高精度地计算特征量。
标号说明
1通信***;10终端数量估计装置;11位置数据取得部;12蓄积部;13观测期间取得部;14观测区域取得部;15观测对象取得部;16前后位置数据取得部;17特征量计算部;17A特征量保存部;18终端数量估计部;19终端数量输出部;21人口估计部;22人口输出部;23属性/放大系数存储部;24变换部;25非识别化部;26隐匿处理部;100便携式终端;200BTS;300RNC;400交换机;500管理中心;501社会感测单元;502PETA挖掘单元;503流动人口统计单元;504可视化处理单元;700各种处理节点。
Claims (26)
1.一种终端数量估计装置,其特征在于,该终端数量估计装置具有:
位置数据取得单元,其取得包含识别便携式终端的识别信息、有关所述便携式终端的位置的位置信息和取得所述位置信息的位置取得时刻信息的位置数据;
前后位置数据取得单元,其对于某个第1位置数据,取得包括与该第1位置数据相同的识别信息的位置数据中、该第1位置数据的紧前面的位置数据即第2位置数据的位置取得时刻信息、以及该第1位置数据的紧后面的位置数据即第3位置数据的位置取得时刻信息;
特征量计算单元,其根据所述第1位置数据的位置取得时刻信息、所述第2位置数据的位置取得时刻信息及所述第3位置数据的位置取得时刻信息中两个以上的位置取得时刻信息,计算有关所述第1位置数据的特征量;
观测对象取得单元,其取得包括有关待观测的观测期间的观测开始时刻以后而且是观测结束时刻以前的位置取得时刻信息、而且包括与有关待观测的观测区域的观测区域信息相对应的位置信息在内的一个或多个位置数据,作为观测对象位置数据;以及
终端数量估计单元,其根据有关所述观测对象位置数据的特征量以及所述观测开始时刻与所述观测结束时刻之差即观测期间长度,估计在观测期间中位于观测区域内的终端数量。
2.根据权利要求1所述的终端数量估计装置,其特征在于,
所述特征量计算单元计算所述第2位置数据的位置取得时刻与所述第3位置数据的位置取得时刻之差,作为有关所述第1位置数据的特征量,
所述终端数量估计单元将有关所述观测对象位置数据的特征量的总和除以所述观测期间长度的2倍,将相除得到的数值估计为所述终端数量。
3.根据权利要求2所述的终端数量估计装置,其特征在于,
在所述第1位置数据的位置取得时刻与所述第2位置数据的位置取得时刻之差大于预定值的情况下,所述特征量计算单元将从所述第1位置数据的位置取得时刻起向过去追溯预定的时间得到的时刻,用作所述第2位置数据的位置取得时刻,来计算有关所述第1位置数据的特征量,
在所述第1位置数据的位置取得时刻与所述第3位置数据的位置取得时刻之差大于预定值的情况下,所述特征量计算单元将从所述第1位置数据的位置取得时刻起向未来前进预定的时间得到的时刻,用作所述第3位置数据的位置取得时刻,来计算有关所述第1位置数据的特征量。
4.根据权利要求1所述的终端数量估计装置,其特征在于,
所述特征量计算单元判别所述第1位置数据是否包括由于跨越位置登记区域边界而生成的位置登记信息、以及所述第3位置数据是否包括由于跨越位置登记区域边界而生成的位置登记信息,
所述特征量计算单元根据有关所述第1位置数据是否包括由于跨越位置登记区域边界而生成的位置登记信息的判别结果、和有关所述第3位置数据是否包括由于跨越位置登记区域边界而生成的位置登记信息的判别结果,使用所述第1位置数据的位置取得时刻信息、所述第2位置数据的位置取得时刻信息及所述第3位置数据的位置取得时刻信息中的两个以上的位置取得时刻信息,计算有关所述第1位置数据的特征量。
5.根据权利要求4所述的终端数量估计装置,其特征在于,
在所述第1位置数据包括由于跨越位置登记区域边界而生成的位置登记信息的情况下,所述特征量计算单元将所述第1位置数据的位置取得时刻设定为第1变量,
在所述第1位置数据不包括由于跨越位置登记区域边界而生成的位置登记信息的情况下,所述特征量计算单元将所述第1位置数据的位置取得时刻与所述第2位置数据的位置取得时刻之间的中点时刻设定为第1变量,
在所述第3位置数据包括由于跨越位置登记区域边界而生成的位置登记信息的情况下,所述特征量计算单元将所述第3位置数据的位置取得时刻设定为第2变量,
在所述第3位置数据不包括由于跨越位置登记区域边界而生成的位置登记信息的情况下,所述特征量计算单元将所述第1位置数据的位置取得时刻与所述第3位置数据的位置取得时刻之间的中点时刻设定为第2变量,
所述特征量计算单元根据所设定的第1变量与第2变量之差来计算有关所述第1位置数据的特征量。
6.根据权利要求4或5所述的终端数量估计装置,其特征在于,
在所述第1位置数据的位置取得时刻与第1变量之差大于预定值的情况下,所述特征量计算单元将从所述第1位置数据的位置取得时刻起向过去追溯预定的时间得到的时刻用作第1变量,来计算有关所述第1位置数据的特征量,
在所述第1位置数据的位置取得时刻与第2变量之差大于预定值的情况下,所述特征量计算单元将从所述第1位置数据的位置取得时刻起向未来前进预定的时间得到的时刻用作第2变量,来计算有关所述第1位置数据的特征量。
7.根据权利要求1~6中任意一项所述的终端数量估计装置,其特征在于,
所述前后位置数据取得单元将由所述观测对象取得单元取得的各个观测对象位置数据作为所述第1位置数据,并取得有关该第1位置数据的第2位置数据的位置取得时刻信息、以及有关该第1位置数据的第3位置数据的位置取得时刻信息,
所述特征量计算单元计算有关各个所述观测对象位置数据的特征量,
所述终端数量估计单元使用计算得到的有关所述观测对象位置数据的特征量来估计所述终端数量。
8.根据权利要求1~6中任意一项所述的终端数量估计装置,其特征在于,
所述前后位置数据取得单元将由所述位置数据取得单元取得的所有位置数据的各个作为所述第1位置数据,并取得有关该第1位置数据的第2位置数据的位置取得时刻信息、以及有关该第1位置数据的第3位置数据的位置取得时刻信息,
所述特征量计算单元计算有关所述所有位置数据各自的特征量,
所述终端数量估计单元使用计算得到的有关所述所有位置数据各自的特征量中有关所述观测对象位置数据的特征量来估计所述终端数量。
9.一种终端数量估计装置,其特征在于,该终端数量估计装置具有:
位置数据取得单元,其取得包含识别便携式终端的识别信息、有关所述便携式终端的位置的位置信息和取得所述位置信息的位置取得时刻信息的位置数据;
前后位置数据取得单元,其对于某个第1位置数据取得包括与该第1位置数据相同的识别信息的位置数据中、该第1位置数据的紧前面的位置数据即第2位置数据的位置取得时刻信息、以及该第1位置数据的紧后面的位置数据即第3位置数据的位置取得时刻信息;
特征量计算单元,其根据至少所述第2位置数据的位置取得时刻及所述第3位置数据的位置取得时刻,计算有关所述第1位置数据的特征量;
观测对象取得单元,其取得包括有关待观测的观测期间的观测开始时刻以后而且是观测结束时刻以前的位置取得时刻信息、而且包括与有关待观测的观测区域的观测区域信息相对应的位置信息在内的一个或多个位置数据,作为观测对象位置数据;以及
终端数量估计单元,其根据有关所述观测对象位置数据的特征量以及所述观测开始时刻与所述观测结束时刻之差即观测期间长度,估计在观测期间中位于观测区域内的终端数量。
10.根据权利要求1~9中任意一项所述的终端数量估计装置,其特征在于,
所述终端数量估计装置还具有存储用于将终端数量变换为人口的放大系数的放大系数存储单元,
所述终端数量估计单元根据有关所述观测对象位置数据的特征量、所述观测期间长度及所述放大系数,估计在观测期间中位于观测区域的人口、以及每人口估计单位的人口中的至少一方,其中该人口估计单位是估计人口的单位。
11.根据权利要求10所述的终端数量估计装置,其特征在于,所述放大系数是使用根据特征量及观测期间长度估计出的终端数量而导出的。
12.根据权利要求1~11中任意一项所述的终端数量估计装置,其特征在于,所述终端数量估计装置还具有变换单元,该变换单元根据与观测区域不同的输出单位和观测区域的重合区域在观测区域中所占的面积比,将通过所述终端数量估计单元的估计而得到的每个观测区域的估计值变换为每输出单位的估计值。
13.根据权利要求12所述的终端数量估计装置,其特征在于,在室内站的通信区域以及利用电波到达范围不同的频带的多个室外站的通信区域中有两个以上的通信区域重叠存在于地理上相同的观测区域的情况下,所述变换单元对于重叠存在的各个通信区域,基于所述面积比变换为每输出单位的估计值、并针对各个所述通信区域对变换后的估计值进行合计,由此得到每输出单位的估计值。
14.根据权利要求1~13中任意一项所述的终端数量估计装置,其特征在于,所述终端数量估计单元根据有关所述观测对象位置数据的特征量、所述观测期间长度、用于将终端数量变换为人口的放大系数、以及与观测区域不同的输出单位和观测区域的重合区域在观测区域中所占的面积比,区分为每输出单位及每人口估计单位来估计人口。
15.根据权利要求14所述的终端数量估计装置,其特征在于,
在所述终端数量估计单元估计人口之前,使所述特征量、所述放大系数、所述面积比以及与该面积比相关的输出单位ID的组合与位置数据相关联,
所述终端数量估计单元对于与相同输出单位ID相关联的位置数据计算(特征量×面积比×放大系数),并按人口估计单位来对所得到的每输出单位的(特征量×面积比×放大系数)进行累计,根据所得到的每输出单位且每人口估计单位的累计值和观测期间长度,估计每输出单位且每人口估计单位的人口。
16.根据权利要求1~15中任意一项所述的终端数量估计装置,其特征在于,所述终端数量估计装置还具有:
观测期间取得单元,其取得包括观测开始时刻与观测结束时刻的组的观测期间信息;以及
观测区域取得单元,其取得与一个或多个位置信息相对应的观测区域信息。
17.根据权利要求1~16中任意一项所述的终端数量估计装置,其特征在于,所述终端数量估计装置还具有输出所得到的估计值的输出单元。
18.根据权利要求17所述的终端数量估计装置,其特征在于,所述输出单元的输出形式能够采用表示人口分布的图、表示按时序人口变动的图、以及表示人口构成的图中的至少一种,所述输出单元的输出单位能够根据便携式终端的用户的属性、时间段、场所中至少一种来设定。
19.根据权利要求1~18中任意一项所述的终端数量估计装置,其特征在于,
所述终端数量估计装置还具有非识别化单元,该非识别化单元针对由所述位置数据取得单元取得的位置数据中所包含的识别信息,进行包括向基于单向性函数的不可逆编码的变换在内的非识别化处理,
所述非识别化单元在进行使用了便携式终端的用户的属性信息的处理的情况下,在该处理之前对所述属性信息进行所述非识别化处理。
20.根据权利要求1~19中任意一项所述的终端数量估计装置,其特征在于,所述终端数量估计装置还具有隐匿处理单元,在所得到的估计值被输出之前,该隐匿处理单元根据预定的基准对所述估计值进行隐匿处理。
21.根据权利要求20所述的终端数量估计装置,其特征在于,所述隐匿处理单元判定取得源终端数量是否小于用于判定为需要隐匿处理的基准值,在某个区域的位置数据的取得源终端数量小于基准值的情况下,将有关该区域的估计值隐匿,其中所述取得源终端数量表示作为估计的基础的每个区域的位置数据是从哪台终端取得的。
22.根据权利要求21所述的终端数量估计装置,其特征在于,在某个区域的位置数据的取得源终端数量不小于基准值的情况下,所述隐匿处理单元根据在估计值输出中使用的多个级别中该区域的估计值所属的级别中的上限值、下限值、级幅以及该估计值,将该区域的估计值圆整为分别伴随有概率值的所述上限值及所述下限值,其中所述概率值对应于与所述上限值的差分以及与所述下限值的差分。
23.一种由终端数量估计装置执行的终端数量估计方法,其特征在于,该终端数量估计方法包括:
位置数据取得步骤,取得包含识别便携式终端的识别信息、有关所述便携式终端的位置的位置信息和取得所述位置信息的位置取得时刻信息的位置数据;
观测对象取得步骤,取得包括有关待观测的观测期间的观测开始时刻以后而且是观测结束时刻以前的位置取得时刻信息、而且包括与有关待观测的观测区域的观测区域信息相对应的位置信息在内的一个或多个位置数据,作为观测对象位置数据;
前后位置数据取得步骤,将所取得的各个观测对象位置数据作为第1位置数据,对于各个第1位置数据取得包括与该第1位置数据相同的识别信息的位置数据中、该第1位置数据的紧前面的位置数据即第2位置数据的位置取得时刻信息、以及该第1位置数据的紧后面的位置数据即第3位置数据的位置取得时刻信息;
特征量计算步骤,根据所述第1位置数据的位置取得时刻信息、所述第2位置数据的位置取得时刻信息及所述第3位置数据的位置取得时刻信息中两个以上的位置取得时刻信息,计算有关各个观测对象位置数据的特征量;以及
终端数量估计步骤,根据有关计算得到的所述观测对象位置数据的特征量以及所述观测开始时刻与所述观测结束时刻之差即观测期间长度,估计在观测期间中位于观测区域内的终端数量。
24.一种由终端数量估计装置执行的终端数量估计方法,其特征在于,该终端数量估计方法包括:
位置数据取得步骤,取得包含识别便携式终端的识别信息、有关所述便携式终端的位置的位置信息和取得所述位置信息的位置取得时刻信息的位置数据;
前后位置数据取得步骤,将所取得的所有位置数据分别作为第1位置数据,对于各个第1位置数据取得包括与该第1位置数据相同的识别信息的位置数据中、该第1位置数据的紧前面的位置数据即第2位置数据的位置取得时刻信息、以及该第1位置数据的紧后面的位置数据即第3位置数据的位置取得时刻信息;
特征量计算步骤,根据所述第1位置数据的位置取得时刻信息、所述第2位置数据的位置取得时刻信息及所述第3位置数据的位置取得时刻信息中两个以上的位置取得时刻信息,计算有关所述所有位置数据各自的特征量;
观测对象取得步骤,取得包括有关待观测的观测期间的观测开始时刻以后而且是观测结束时刻以前的位置取得时刻信息、而且包括与有关待观测的观测区域的观测区域信息相对应的位置信息在内的一个或多个位置数据,作为观测对象位置数据;以及
终端数量估计步骤,根据有关计算得到的所述所有位置数据各自的特征量中有关所述观测对象位置数据的特征量、以及所述观测开始时刻与所述观测结束时刻之差即观测期间长度,估计在观测期间中位于观测区域内的终端数量。
25.一种由终端数量估计装置执行的终端数量估计方法,其特征在于,该终端数量估计方法包括:
位置数据取得步骤,取得包含识别便携式终端的识别信息、有关所述便携式终端的位置的位置信息和取得所述位置信息的位置取得时刻信息的位置数据;
观测对象取得步骤,取得包括有关待观测的观测期间的观测开始时刻以后而且是观测结束时刻以前的位置取得时刻信息、而且包括与有关待观测的观测区域的观测区域信息相对应的位置信息在内的一个或多个位置数据,作为观测对象位置数据;
前后位置数据取得步骤,将所取得的各个观测对象位置数据作为第1位置数据,对于各个第1位置数据取得包括与该第1位置数据相同的识别信息的位置数据中、该第1位置数据的紧前面的位置数据即第2位置数据的位置取得时刻信息、以及该第1位置数据的紧后面的位置数据即第3位置数据的位置取得时刻信息;
特征量计算步骤,根据至少所述第2位置数据的位置取得时刻及所述第3位置数据的位置取得时刻,计算有关各个观测对象位置数据的特征量;以及
终端数量估计步骤,根据有关计算得到的所述观测对象位置数据的特征量以及所述观测开始时刻与所述观测结束时刻之差即观测期间长度,估计在观测期间中位于观测区域内的终端数量。
26.一种由终端数量估计装置执行的终端数量估计方法,其特征在于,该终端数量估计方法包括:
位置数据取得步骤,取得包含识别便携式终端的识别信息、有关所述便携式终端的位置的位置信息和取得所述位置信息的位置取得时刻信息的位置数据;
前后位置数据取得步骤,将所取得的所有位置数据分别作为第1位置数据,对于各个第1位置数据,取得包括与该第1位置数据相同的识别信息的位置数据中、该第1位置数据的紧前面的位置数据即第2位置数据的位置取得时刻信息、以及该第1位置数据的紧后面的位置数据即第3位置数据的位置取得时刻信息;
特征量计算步骤,根据至少所述第2位置数据的位置取得时刻及所述第3位置数据的位置取得时刻,计算有关所述所有位置数据各自的特征量;
观测对象取得步骤,取得包括有关待观测的观测期间的观测开始时刻以后而且是观测结束时刻以前的位置取得时刻信息、而且包括与有关待观测的观测区域的观测区域信息相对应的位置信息在内的一个或多个位置数据,作为观测对象位置数据;以及
终端数量估计步骤,根据有关计算得到的所述所有位置数据各自的特征量中有关所述观测对象位置数据的特征量、以及所述观测开始时刻与所述观测结束时刻之差即观测期间长度,估计在观测期间中位于观测区域内的终端数量。
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010240551 | 2010-10-27 | ||
JP2010-240551 | 2010-10-27 | ||
JP2011-084796 | 2011-04-06 | ||
JP2011084796 | 2011-04-06 | ||
PCT/JP2011/073584 WO2012056900A1 (ja) | 2010-10-27 | 2011-10-13 | 端末数推計装置および端末数推計方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102959996A true CN102959996A (zh) | 2013-03-06 |
CN102959996B CN102959996B (zh) | 2015-07-22 |
Family
ID=45993625
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201180030251.9A Expired - Fee Related CN102959996B (zh) | 2010-10-27 | 2011-10-13 | 终端数量估计装置及终端数量估计方法 |
Country Status (8)
Country | Link |
---|---|
US (1) | US8942724B2 (zh) |
EP (1) | EP2635053A4 (zh) |
JP (1) | JP5470469B2 (zh) |
KR (1) | KR101456906B1 (zh) |
CN (1) | CN102959996B (zh) |
BR (1) | BR112013000913A2 (zh) |
RU (1) | RU2530349C1 (zh) |
WO (1) | WO2012056900A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107423912A (zh) * | 2017-08-02 | 2017-12-01 | 长威信息科技发展股份有限公司 | 一种基于***网格动态划分的方法和*** |
CN109587629A (zh) * | 2018-12-18 | 2019-04-05 | 浙江每日互动网络科技股份有限公司 | 一种基于移动终端WiFi位置数据的突发灾害判断方法 |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10200844B2 (en) * | 2012-08-23 | 2019-02-05 | Nokia Solutions And Networks Oy | Massive discovery of devices |
JP6048204B2 (ja) * | 2013-02-20 | 2016-12-21 | 株式会社リコー | 情報処理装置、情報処理システム、及び情報処理プログラム |
JP6348019B2 (ja) * | 2014-08-28 | 2018-06-27 | ルネサスエレクトロニクス株式会社 | 通信システム、通信装置、自動車および通信方法 |
JP6110960B2 (ja) | 2014-10-24 | 2017-04-05 | 株式会社Agoop | 人口推定装置、プログラム及び人口推定方法 |
JP6130524B2 (ja) * | 2014-10-24 | 2017-05-17 | 株式会社Agoop | 人口推定装置、プログラム及び人口推定方法 |
US9615242B2 (en) | 2015-07-07 | 2017-04-04 | T-Mobile Usa, Inc. | Determining a service leakage rate within a wireless communication network |
EP3142393A1 (en) * | 2015-09-14 | 2017-03-15 | BASE Company | Method and system for obtaining demographic information |
JP6634511B2 (ja) | 2016-03-30 | 2020-01-22 | 株式会社Agoop | 人口推定装置、プログラム及び人口推定方法 |
AU2018320434B2 (en) * | 2017-08-22 | 2020-12-17 | Nippon Telegraph And Telephone Corporation | Secure computation device, secure computation method, program, and recording medium |
JP7060537B2 (ja) * | 2019-03-18 | 2022-04-26 | ソフトバンク株式会社 | 情報処理システム、移動局数推計方法、人口推計方法及びプログラム |
US20230169414A1 (en) * | 2020-04-23 | 2023-06-01 | Ntt Docomo, Inc. | Population extraction device |
JP7142058B2 (ja) * | 2020-07-28 | 2022-09-26 | ヤフー株式会社 | 情報処理装置、情報処理方法および情報処理プログラム |
GB2623492A (en) * | 2022-10-10 | 2024-04-24 | Nokia Technologies Oy | Estimating position of user device |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003122877A (ja) * | 2001-10-09 | 2003-04-25 | Oki Electric Ind Co Ltd | 地域内人口算出システム |
JP4364936B1 (ja) * | 2009-02-27 | 2009-11-18 | ソフトバンクBb株式会社 | 流動性人口推定システム、流動性人口推定方法及び流動性人口推定プログラム |
WO2010116900A1 (ja) * | 2009-04-06 | 2010-10-14 | 株式会社エヌ・ティ・ティ・ドコモ | 通信システム、情報分析装置、及び情報分析方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4442007B2 (ja) * | 2000-03-30 | 2010-03-31 | Kddi株式会社 | 移動端末にグループ通信を提供するシステム |
US7181225B1 (en) * | 2000-12-19 | 2007-02-20 | Bellsouth Intellectual Property Corporation | System and method for surveying wireless device users by location |
JP2003044969A (ja) | 2001-08-01 | 2003-02-14 | Marubeni Energy Corp | 緊急対策用ネットワークシステム |
US7123918B1 (en) * | 2001-08-20 | 2006-10-17 | Verizon Services Corp. | Methods and apparatus for extrapolating person and device counts |
US6865395B2 (en) * | 2002-08-08 | 2005-03-08 | Qualcomm Inc. | Area based position determination for terminals in a wireless network |
KR100595585B1 (ko) * | 2003-05-13 | 2006-07-03 | 엘지전자 주식회사 | 이동통신시스템에서의 멀티미디어 방송 및 멀티캐스트서비스를 위한 무선자원관리방법 |
US20060221918A1 (en) * | 2005-04-01 | 2006-10-05 | Hitachi, Ltd. | System, method and computer program product for providing content to a remote device |
US7546120B1 (en) * | 2005-10-14 | 2009-06-09 | Sprint Spectrum L.P. | Method and system for managing transmission of media to multiple subscribers |
US8688396B2 (en) * | 2010-02-02 | 2014-04-01 | Microsoft Corporation | Targeted data collection for positioning systems |
US20110285591A1 (en) * | 2010-05-19 | 2011-11-24 | Palm, Inc. | Correlating contextual data and position data to improve location based services |
US8478289B1 (en) * | 2010-06-01 | 2013-07-02 | Google Inc. | Predicting geographic population density |
-
2011
- 2011-10-13 KR KR1020127030467A patent/KR101456906B1/ko not_active IP Right Cessation
- 2011-10-13 BR BR112013000913A patent/BR112013000913A2/pt not_active IP Right Cessation
- 2011-10-13 WO PCT/JP2011/073584 patent/WO2012056900A1/ja active Application Filing
- 2011-10-13 JP JP2012540768A patent/JP5470469B2/ja active Active
- 2011-10-13 RU RU2013109951/08A patent/RU2530349C1/ru not_active IP Right Cessation
- 2011-10-13 EP EP11836039.5A patent/EP2635053A4/en not_active Withdrawn
- 2011-10-13 US US13/704,446 patent/US8942724B2/en not_active Expired - Fee Related
- 2011-10-13 CN CN201180030251.9A patent/CN102959996B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003122877A (ja) * | 2001-10-09 | 2003-04-25 | Oki Electric Ind Co Ltd | 地域内人口算出システム |
JP4364936B1 (ja) * | 2009-02-27 | 2009-11-18 | ソフトバンクBb株式会社 | 流動性人口推定システム、流動性人口推定方法及び流動性人口推定プログラム |
WO2010116900A1 (ja) * | 2009-04-06 | 2010-10-14 | 株式会社エヌ・ティ・ティ・ドコモ | 通信システム、情報分析装置、及び情報分析方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107423912A (zh) * | 2017-08-02 | 2017-12-01 | 长威信息科技发展股份有限公司 | 一种基于***网格动态划分的方法和*** |
CN109587629A (zh) * | 2018-12-18 | 2019-04-05 | 浙江每日互动网络科技股份有限公司 | 一种基于移动终端WiFi位置数据的突发灾害判断方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102959996B (zh) | 2015-07-22 |
KR101456906B1 (ko) | 2014-10-31 |
KR20130010017A (ko) | 2013-01-24 |
JP5470469B2 (ja) | 2014-04-16 |
WO2012056900A1 (ja) | 2012-05-03 |
US20130090132A1 (en) | 2013-04-11 |
EP2635053A1 (en) | 2013-09-04 |
EP2635053A4 (en) | 2017-04-26 |
US8942724B2 (en) | 2015-01-27 |
BR112013000913A2 (pt) | 2016-05-17 |
RU2530349C1 (ru) | 2014-10-10 |
JPWO2012056900A1 (ja) | 2014-03-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102959996A (zh) | 终端数量估计装置及终端数量估计方法 | |
CN103154980A (zh) | 加入终端数量估计装置及加入终端数量估计方法 | |
US20110055216A1 (en) | Tempo spatial data extraction from network connected devices | |
CN103348705A (zh) | 终端数量推算装置和终端数量推算方法 | |
US8903791B2 (en) | Statistical information generation system and statistical information generation method | |
KR20130052014A (ko) | 단말기 수 추계 장치 및 단말기 수 추계 방법 | |
CN106464706A (zh) | 用于通过可从电信网络获得的数据识别显著地点的方法和*** | |
KR101411451B1 (ko) | 관리 서버, 통신 시스템 및 통계 처리 방법 | |
CN109885636B (zh) | 一种用户画像方法和服务器 | |
CN103039116A (zh) | 信息分析装置和信息分析方法 | |
CN103037424A (zh) | 一种3g网络覆盖的评估方法及装置 | |
JP2012054921A (ja) | 移動機分布算出システム及び移動機分布算出方法 | |
CN103120001A (zh) | 终端数量估计装置及终端数量估计方法 | |
CN101631343B (zh) | 一种识别移动用户归属区域的方法及装置 | |
CN102638821A (zh) | 基于手机基站数据挖掘人的移动模式的方法 | |
Yakimov | Features of the Use of Geoanalytical Data in the Development of Transport Planning Documents | |
Lind et al. | Human mobility patterns at the smallest scales | |
AU2014386813A1 (en) | Method and system for generating a user activity grid table | |
Chen et al. | Forecasting Individual Demand in Cellular Networks | |
Xie et al. | Pattern mining, semantic label identification and movement prediction using mobile phone data | |
JP2013138417A (ja) | 情報分析装置および情報分析方法 | |
Hong et al. | A grid-based node split algorithm for managing current location data of moving objects |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150722 Termination date: 20171013 |
|
CF01 | Termination of patent right due to non-payment of annual fee |