CN107113560A - 用于实时计数公共集会处的参与者的数量的方法和*** - Google Patents

用于实时计数公共集会处的参与者的数量的方法和*** Download PDF

Info

Publication number
CN107113560A
CN107113560A CN201480084457.3A CN201480084457A CN107113560A CN 107113560 A CN107113560 A CN 107113560A CN 201480084457 A CN201480084457 A CN 201480084457A CN 107113560 A CN107113560 A CN 107113560A
Authority
CN
China
Prior art keywords
interest
area
radius
user equipment
time
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
Application number
CN201480084457.3A
Other languages
English (en)
Other versions
CN107113560B (zh
Inventor
M·科洛纳
M·麻美
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Telecom Italia SpA
Original Assignee
Telecom Italia SpA
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Telecom Italia SpA filed Critical Telecom Italia SpA
Publication of CN107113560A publication Critical patent/CN107113560A/zh
Application granted granted Critical
Publication of CN107113560B publication Critical patent/CN107113560B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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/021—Services related to particular areas, e.g. point of interest [POI] services, venue services or geofences
    • H04W4/022—Services related to particular areas, e.g. point of interest [POI] services, venue services or geofences with dynamic range variability
    • 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
    • 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
    • G06Q30/00—Commerce
    • G06Q30/02—Marketing; Price estimation or determination; Fundraising
    • G06Q30/0201—Market modelling; Market analysis; Collecting market data
    • G06Q30/0204—Market segmentation
    • G06Q30/0205—Location or geographical consideration
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W84/00—Network topologies
    • H04W84/02—Hierarchically pre-organised networks, e.g. paging networks, cellular networks, WLAN [Wireless Local Area Network] or WLL [Wireless Local Loop]
    • H04W84/04—Large scale networks; Deep hierarchical networks
    • H04W84/042—Public Land Mobile systems, e.g. cellular systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Mobile Radio Communication Systems (AREA)
  • Telephonic Communication Services (AREA)

Abstract

提出了一种估计为了参加在日期(g)的时间间隔([Ts；Te])期间的公共集会(S)而聚集在感兴趣区域(107)处的人数的方法。所述感兴趣区域(107)由感兴趣区域中心(C)和感兴趣区域半径(Ra)来定义，并且由具有多个通信站(105a)的移动电信网络(105)覆盖，每个通信站适于管理其中细分移动电信网络(105)的一个或多个服务区域(105b)中的用户设备的通信。该方法包括以下步骤：a)定义(504、506、524‑528)感兴趣区域半径(Ra)的多个计算出的半径值(Rk)，以及对于每个计算出的半径值(Rk)：b)识别(506)包含在感兴趣区域(107)中的服务区域(105b)中的相关服务区域(405a‑d)的数量；c)识别(510、512)与在感兴趣区域(107)内的公共集会的日期(g)的时间间隔([Ts；Te])期间在用户设备和移动通信网络(105)之间发生的交互的对应事件(ei)的至少一个事件记录(eri)相关联的用户设备的第一数量(Uk)；d)识别(514、516)在感兴趣区域(107)内的公共集会的日期(g)之前预定数量(P)个先前日期(gp)中的每一个的时间间隔([Ts；Te])期间在用户设备和移动通信网络(105)之间发生的交互的对应事件(ei)的至少一个事件记录(eri)相关联的用户设备的第二数量(Upk)；e)将用户设备的第一数量(Uk)和用户设备的第二数量(Upk)组合(518、520)以获得统计量(Zk)；f)通过相对于相关服务区域(405a‑d)的半径的总和(Rtotk)归一化(522)统计量(Zk)来计算归一化统计量(Z'k)；g)将感兴趣区域半径(Ra)的最佳半径值(Ro)计算(530)为由归一化统计量(Z'k)加权的计算出的半径值(Rk)的平均值；h)计数(532‑566)为了参加在具有等于最佳半径值(Ro)的感兴趣区域半径(Ra)的感兴趣区域内的公共集会(S)而聚集的人数。

Description

用于实时计数公共集会处的参与者的数量的方法和***
技术领域
本发明涉及人群计数，即，涉及用于计数或估计人群中的人数的技术。在本说明书中，为了本发明的目的，通过“人群”意指聚集在特定地点(例如，参加具有最不同性质的公共事件或集会)的特定数量的人的聚集，像例如(并且非穷举)直播电视公共集会、艺术/娱乐表演、文化展览、戏剧演出、体育比赛、音乐会、电影、示威等。
特别地，本发明涉及利用由无线或移动电信网络提供的信息的人群计数技术。
背景技术
在城市规划、活动管理(例如，运输***管理和紧急情况管理)以及旅游和本地营销的任务中，了解在特定地点或感兴趣区域(简称AoI，例如，建筑物，诸如例如体育场或剧院或电影院、其周围环境、城市或城镇或乡村的广场或(一条或多条)街道、辖区等)处聚集的人的数量是有用的，例如，因为他们参加了在感兴趣区域内发生的(例如，与文化、娱乐、政治或体育相关的)比如节目的公共集会。事实上，这种知识允许例如更有效地规划后续的相同类型的公共集会。特别地，这种知识允许更有效地规划和管理与未来可能发生的类似的公共集会(诸如，例如，在体育场定期发生的体育比赛)直接或间接相关的资源和活动(诸如基础设施、运输***和安全)。此外，从商业角度来看，这种知识允许更好地管理旨在促进未来可能发生的类似事件的营销活动。
如今，移动通信装置(以下称为移动电话或UE，包括蜂窝电话、智能电话、平板电脑等)已经在许多国家的人口中得到彻底的传播，并且移动电话所有者几乎始终携带其移动电话与他们在一起。由于移动电话与移动电话网络的多个基站通信，并且每个基站覆盖(即，服务)移动通信服务提供商(例如移动电话网络所有者或虚拟移动电话服务提供商)已知的一个或多个预定服务区域或小区，因此移动电话结果成为作为跟踪装置的最佳候选者，用于收集对识别参加一个或多个公共集会的人的数量有用的数据。
在本领域中，为了收集关于个人的用户设备(UE，例如，移动电话、智能电话、平板电脑等)连接到移动电话网络(例如，用于进行语音呼叫或发送文本消息)的时间和地点的信息，并且使用这种收集到的信息以便得到与特定公共集会具有多少参加者相关的信息，已经提出了许多***和方法。
例如，2006年的Francesco Calabrese，Carlo Ratti，“Real Time Rome”，Networks and Communications Studies 20(3-4)，第247-258页公开了在第十届意大利威尼斯国际建筑展(International Architecture Exhibition)上呈现的实时罗马(RealTime Rome)项目。实时罗马项目是收集并处理由电信网络和运输***提供的数据以便了解罗马日常生活的模式的城市范围实时监视***的第一个示例。观察城镇中的实时日常生活成为了解现在和预测未来城市环境的手段。
从http://mox.polimi.it可获得的F.Manfredini，P.Pucci，P.Secchi，P.Tagliolato，S.Vantini，V.Vitelli，“Treelet decomposition of mobile phone datafor deriving city usage and mobility pattern in the Milan urban region”，MOX-Report No.25/2012，MOX，Department of Mathematics“F.Brioschi”，Politecnico diMilano公开了针对识别关于移动电话使用的隐藏模式的有用信息的地理统计学无监督学习技术。这些隐藏的模式涉及城市在时间和空间上与个人移动性相关的不同使用，概述了该技术用于城市规划社区的潜力。该方法允许获得参考基础，该参考基础报告一些活动对所记录的Erlang数据的特定影响以及显示每个活动对本地Erlang信号的贡献的一组映射。已经选择了对于解释特定移动性和城市使用模式(通勤、夜间活动、住宅分布、非***移动性)有重要意义的结果，并且已在米兰城市地区规模从城市分析和规划视角对其意义及其解释进行了测试。
2012年6月的第二届普及城市应用(Pervasive Urban Applications，PURBA)研讨会，Ramon Caceres，James Rowland，Christopher Small和Simon Urbanek，“Exploringthe Use of Urban Greenspace through Cellular Network Activity”公开了使用蜂窝网络活动的匿名记录来研究城市区域中的人口密度的时空模式。该文章呈现了该努力的愿景和一些早期的结果。首先，描述了纽约大都市地区中的六个月的活动的数据集。第二，呈现了用于估计网络覆盖区域的技术。第三，描述了分析在那些区域内的活动量变化所使用的途径。最后，呈现了关于中央公园周围人口密度变化的初步结果。
2009年的F.Girardin，A.Gerber，A.Vaccari，A.Biderman，C.Ratti，“Towardsestimating the presence of visitors from the aggregate mobile phone networkactivity they generate”，International Conference on Computers in UrbanPlanning and Urban Management审查了纽约市在2008年的“Waterfalls”公共展览附近的本地和非本地注册移动电话的使用。研究了与覆盖展览及其邻近地区的网络部分相关的总体统计(即，呼叫的数量)。利用传统调查技术(诸如现场计数)的未来贡献来校准这些移动电话网络测量，其目的是开发技术来估计游客在时间和空间上的总体移动和地点，同时确保其隐私。
F.Calabrese，F.C.Pereira，G.Di Lorenzo，L.Liu，C.Ratti，“The Geography ofTaste：Analyzing Cell-Phone Mobility in Social Events”，Pervasive Computing，LNCS 6030，Springer，2010，pp.22-37公开了在特殊事件期间对人群移动性的分析。已经分析了近一百万部蜂窝电话的踪迹，并且将其目的地与社会事件相关联。已经观察到，参加事件的人的来源与事件的类型密切相关，这在城市管理方面具有启示，因为附加的流动的知识对于做出关于事件管理和拥塞缓解的决策会是关键的信息。
2011年10月9-11日的Traag，V.A.；Browet，A.；Calabrese，F.；Morlot，F.，“SocialEvent Detection in Massive Mobile Phone Data Using Probabilistic LocationInference”，2011 IEEE Third International Conference on Privacy，Security，Riskand Trust(Passat)和2011IEEE Third International Conference on SocialComputing(Socialcom)，pp.625,628，专注于人们的异常大型聚集，即，异常社会事件。基于贝叶斯地点推理框架，介绍了在大规模移动电话数据中检测这种社会事件的方法。更具体地，还开发了用于决定谁参加事件的框架。演示了关于几个示例的方法。最后，讨论了用于事件检测的一些可能的未来途径，以及对检测到的社会事件的一些可能的分析。
发明内容
申请人已经观察到，通常，本领域已知的方法和***提供了不能令人满意的结果，因为它们不能够确定(或者在确定方面具有有限的能力)UE所有者是否已经出于参加其中的原因或出于其他原因(例如，由于UE所有者居住在该感兴趣区域附近或其内或者UE所有者在该感兴趣区域附近或其内有业务)而在其中已经举行一个或多个公共集会的感兴趣区域(AoI)中。另外，由已知解决方案提供的结果受到在一个或多个公共集会中被选择用于分析参加者的数量的感兴趣区域的尺寸的强烈影响。换句话说，如果感兴趣区域具有大的尺寸，那么实际上不是人群的一部分的特定数量的UE所有者在公共集会的参加者数量的评估中将被考虑。相反，如果感兴趣区域具有小的尺寸，那么实际上是人群的一部分的特定数量的UE所有者将被排除在人群中的人数的评估之外。
因此，基于由本领域已知的方法和***获得的结果对(上面提到的类型的)资源和活动进行后续的规划和管理将由于其有限的精确度而实现有限的效率。
因此，申请人已经处理了设计适于克服影响现有技术解决方案的问题的***和方法的问题。
申请人已经发现，可以基于在一个或多个公共集会的过程期间和在一个或多个公共集会之前的特定天数中与UE相关的操作信息来确定最佳感兴趣区域的尺寸。
特别地，本发明的一个方面提出了估计为了参加在日期g的时间间隔[Ts；Te]期间的公共集会S而聚集在感兴趣区域处的人数的方法。所述感兴趣区域由感兴趣区域中心C和感兴趣区域半径Ra来定义，并且由具有多个通信站的移动电信网络覆盖，每个通信站适于管理其中细分移动电信网络的一个或多个服务区域中的用户设备的通信。该方法包括以下步骤：a)定义感兴趣区域半径Ra的多个计算出的半径值Rk，以及对于每个计算出的半径值Rk：b)识别包含在感兴趣区域中的服务区域中的相关服务区域的数量；c)识别与在感兴趣区域内的公共集会S的日期g的时间间隔[Ts；Te]期间在用户设备和移动通信网络之间发生的交互的对应事件ei的至少一个事件记录eri相关联的用户设备的第一数量Uk；d)识别在感兴趣区域内的公共集会的日期g之前预定数量P个先前日期gp中的每一个先前日期的时间间隔[Ts；Te]期间在用户设备和移动通信网络之间发生的交互的对应事件ei’的至少一个事件记录eri’相关联的用户设备的第二数量Upk；e)将用户设备的第一数量Uk和用户设备的第二数量Upk组合以获得统计量Zk；f)通过相对于相关服务区域的半径归一化统计量Zk来计算归一化统计量Z'k；g)将感兴趣区域半径Ra的最佳半径值Ro计算为由归一化统计量Z'k加权的计算出的半径值Rk的平均值；h)计数为了参加在具有等于最佳半径值Ro的感兴趣区域半径Ra的感兴趣区域内的公共集会S而聚集的人数。
本发明的优选特征在从属权利要求中阐述。
在本发明的一个实施例中，通过相对于相关服务区域的半径归一化统计量Zk来计算归一化统计量Z'k的步骤f)包括通过相对于相关服务区域的半径的总和Rtotk归一化统计量Zk来计算归一化统计量Z'k。
在本发明的一个实施例中，如果服务区域验证了以下不等式，那么该服务区域被识别为相关服务区域：
Dist(C,B)≤|Rc+Rk|，
其中C是感兴趣区域的中心，B是服务区域的中心，Dist(C,B)是感兴趣区域的中心C与服务区域的中心B之间的地理距离，Rc是服务区域的半径，并且Rk是计算出的半径值。
在本发明的一个实施例中，组合用户设备的第一数量Uk和用户设备的第二数量Upk以便获得统计量Zk的步骤d)包括：组合先前日期gpn中的每一个先前日期的用户设备的第二数量Upk，以便确定平均用户设备数量μk和用户设备数量标准偏差σk。
在本发明的一个实施例中，组合用户设备的第一数量Uk和用户设备的第二数量Upk以便获得统计量Zk的步骤d)还包括：将统计量Zk计算为：
Zk＝(Uk–μk)/σk，
其中Uk是第一数量，μk是平均用户设备数量并且σk是用户设备数量标准偏差。
在本发明的一个实施例中，多个计算出的半径值Rk的范围从最小半径值Rmin到最大半径值Rmax，每个计算出的半径值Rk与下一个半径值分开迭代宽度Δ。
在本发明的一个实施例中，计数为了参加在具有等于最佳半径值Ro的感兴趣区域半径Ra的感兴趣区域内的公共集会S而聚集的人数的步骤h)包括：i)定义包含在移动通信网络中的服务区域中的相关服务区域的数量，所述相关服务区域是至少部分地叠加在具有等于最佳半径值Ro的感兴趣区域半径Ra的感兴趣区域上的服务区域。
在本发明的一个实施例中，如果服务区域验证了以下不等式，那么该服务区域被识别为相关服务区域：
Dist(C,B)≤|Rc+Ro|，
其中C是感兴趣区域的中心，B是服务区域的中心，Dist(C,B)是感兴趣区域的中心C与服务区域的中心B之间的地理距离，Rc是服务区域的半径，并且Ro是最佳半径值。
在本发明的一个实施例中，计数为了参加在具有等于最佳半径值Ro的感兴趣区域半径Ra的感兴趣区域内的公共集会S而聚集的人数的步骤h)包括：j)构建用户设备列表uL，所述用户设备列表包括与在感兴趣区域内在日期g的时间间隔[Ts；Te]中发生的事件记录eri中的至少一个事件记录eri相关联的每个用户设备UEj的标识符。
在本发明的一个实施例中，计数为了参加在具有等于最佳半径值Ro的感兴趣区域半径Ra的感兴趣区域内的公共集会S而聚集的人数的步骤h)还包括，对于用户设备列表uL中的每个用户设备UEj：k)计算与用户设备UEj相关联的连续事件记录eri、eri'之间的平均iat中间到达时间。
在本发明的一个实施例中，平均中间到达时间iat基于在日期g和在先前日期gp两者发生的事件记录eri、eri'来计算。
在本发明的一个实施例中，计数为了参加在具有等于最佳半径值Ro的感兴趣区域半径Ra的感兴趣区域内的公共集会S而聚集的人数的步骤h)还包括，对于用户设备列表uL中的每个用户设备UEj：l)识别在感兴趣区域内在日期g的时间间隔[Ts；Te]中发生的、分别参考第一事件记录erf和最后事件记录erI的第一时间数据tdf和最后时间数据tdI，以及m)组合第一时间数据tdf、最后时间数据tdI和平均中间到达时间iat，以便确定第一时间片段f1，所述第一时间片段指示用户设备UEj在日期g在时间间隔[Ts；Te]期间在感兴趣区域内花费的时间段。
在本发明的一个实施例中，第一时间片段被确定为：
其中tdf是第一时间数据，tdI是最后时间数据，iat是平均中间到达时间，Ts是公共集会S的开始时间并且Te是公共集会S的结束时间。
在本发明的一个实施例中，计数为了参加在具有等于最佳半径值Ro的感兴趣区域半径Ra的感兴趣区域内的公共集会S而聚集的人数的步骤h)还包括，对于用户设备列表uL中的每个用户设备UEj：n)识别在感兴趣区域内在先前日期gp的时间间隔[Ts；Te]中发生的、分别参考第一事件记录erf’和最后事件记录erI’的第一先前时间数据tdf’和最后先前时间数据tdI’，以及o)组合第一先前时间数据tdf’、最后先前时间数据tdI’和平均中间到达时间iat，以便确定第二时间片段f2，所述第二时间片段f2指示用户设备UEj在先前日期gp期间在感兴趣区域内花费的时间段。
在本发明的一个实施例中，第二时间片段f2被确定为：
其中tdpf是第一先前时间数据，tdpI是最后先前时间数据，iat是平均中间到达时间，并且Tgp是先前日期gp的总持续时间。
在本发明的一个实施例中，计数为了参加在具有等于最佳半径值Ro的感兴趣区域半径Ra的感兴趣区域内的公共集会S而聚集的人数的步骤h)还包括，对于用户设备列表uL中的每个用户设备UEj：p)组合第一时间片段f1和第二时间片段f2，以便确定用户设备UEj的所有者参加公共集会S的概率pj。
在本发明的一个实施例中，用户设备UEj的所有者参加公共集会S的概率pj被计算为：
pj＝f1*(1–f2)，
其中f1是第一时间片段，并且f2是第二时间片段。
在本发明的一个实施例中，计数为了参加在具有等于最佳半径值Ro的感兴趣区域半径Ra的感兴趣区域内的公共集会S而聚集的人数的步骤h)还包括：q)将为了参加公共集会S而聚集的人数A计算为针对用户设备列表uL的每个用户设备UEj确定的概率pj的总和。
本发明的另一个方面提出了耦合到无线电信网络的用于估计在感兴趣区域处聚集的人数的***。该***包括计算引擎、储存库以及管理员界面，所述计算引擎适于处理从移动电话网络检索到的数据，所述储存库适于存储关于用户设备和移动电话网络之间的交互的数据、由计算引擎生成的计算结果以及可能的由***生成的和/或提供给***的任何处理数据，所述管理员界面可操作用于修改由计算引擎使用的参数和/或算法，和/或访问存储在储存库中的数据。根据本发明的***还包括存储被配置为通过***实现以上方法的软件程序产品的存储器元件。
在本发明的一个实施例中，***还包括至少一个用户界面，所述至少一个用户界面适于接收来自***的用户的输入，并向***的用户提供输出，用户包括一个或多个人和/或由***提供的服务的一个或多个外部计算***订户。
根据本发明的解决方案的优点之一在于它是计算上简单的，仅涉及计数和代数运算的操作。
附图说明
通过阅读仅仅通过非限制性示例的方式提供的、要结合附图阅读的本发明的实施例的以下详细描述，将更好地理解根据本发明的解决方案的这些和其他特征及优点，其中：
图1是根据本发明的实施例的人群估计***的示意图；
图2A至图2E是根据本发明的实施例的其中移动通信网络的小区可以被建模的示例性形状；
图3A至图3E是根据本发明的实施例的待确定的AoI可以采用的示例性形状；
图4A至图4D是根据本发明的实施例的移动通信网络105相对于AoI的小区中的相关小区，以及
图5A至图5C是根据本发明的实施例的人群评估算法的示意流程图。
具体实施方式
参考附图，图1是根据本发明的示例性实施例的在下文中简单地表示为***100的人群估计***的示意图。
人群估计***允许进行为了参加具有最不同性质的一个或多个公共集会而聚集的人群中的人数的估计，像例如(并且非穷举)直播电视公共集会、艺术/娱乐表演、文化展览、戏剧演出、体育比赛、音乐会、电影、示威等。
***100耦合到诸如(2G、3G、4G或更高代)移动电话网络的移动通信网络105，并且被配置用于从移动通信网络105接收位于感兴趣地理区域(简称AoI)中的个人的以下(例如，移动电话、智能电话、具有2G-3G-4G连接的平板电脑等)中的每个用户设备UE的定位数据，感兴趣的地理区域在图1中示意性地示为点划线107内的区域(例如，建筑物及其周围环境、城市的一部分、整个城市、直辖市、辖区等)。
(在下面进一步描述的)AoI 107通常可以包括其中已经发生公共集会的核心地点(例如，体育场、剧院、城市广场等)，并且可能地包括核心地点的周围环境(例如，附近的停车场、附近的街道、附近的运输站等)。
移动通信网络105包括地理上跨AoI 107分布的多个(两个或更多个)通信站105a(例如，移动电话网络的无线电基站)。每个通信站105a适于管理在一个或多个服务区域或小区105b(在所述示例中，每个通信站105a服务三个小区)中的UE(未示出，诸如例如移动电话)的通信，如下面将更详细地讨论的。
甚至更一般地，移动通信网络105的每个通信站105a适于(例如，在电源接通/断开时、在地点区域更新时、在来电/去电呼叫时、在发送/接收SMS和/或MMS时、在因特网接入时等)与位于由这种通信站105a服务的小区105b之一内的任何UE交互。在下文中，UE和移动通信网络105之间的这些交互将通常被表示为事件ei(i＝1，...，I；其中I是整数)。
***100包括被配置为适于处理从移动通信网络105检索到的数据的计算引擎110，和被配置为适于(通常以二进制格式)存储关于UE和移动通信网络105之间的交互的数据、由计算引擎110生成的计算结果以及可能的由***100生成和/或提供给***100的任何处理数据的储存库115(诸如数据库、文件***等)。***100设置有管理员接口120(例如，计算机)，所述管理员接口120被配置并且可操作用于修改由计算引擎110使用的参数和/或算法和/或访问存储在储存库115中的数据。
优选地，***100包括一个或多个用户界面125(例如，用户终端、在连接到***100的远程终端上运行的软件)，所述一个或多个用户界面125适于从***100的用户接收输入并向***100的用户提供输出。如本公开中所使用的术语“***的用户”可以指一个或多个人和/或作为由***100提供的服务的订户并且使得能够访问***100的第三方的外部计算***(诸如计算机网络，未示出)--例如，在与***100的服务提供商所有者签订合同的情况下，并且与由通过管理员接口120操作的***100的管理员拥有的访问权相比，通常具有减少的对***100的访问权。
应当认识到，***100可以以任何已知的方式实现；例如，***100可以包括单个计算机、或分布式计算机的网络，或者是物理类型(例如，具有实现计算引擎110和储存库115的一个或多个主机器，连接到实现管理员和用户界面120和125的其他机器)或者是虚拟类型(例如，通过在计算机网络中实现一个或多个虚拟机)。
***100适于针对在AoI 107内在UE和移动通信网络105之间(通过其通信站105a之一)发生的每个事件ei检索(和/或接收)事件记录eri。优选地，由***100从移动通信网络105检索到的每个事件记录eri包括–以非限制性方式–对应事件ei中涉及的UE的标识符(例如，UE标识符可以被选择作为国际移动设备身份-IMEI、国际移动订户身份-IMSI和移动订户ISDN号-MSISDN码之中的一个或多个)、指示对应事件ei发生的时间的时间数据(也表示为时间戳)以及UE地理位置数据，例如基于在对应事件ei发生时UE所位于的小区105b的空间指示。
在本发明的一个实施例中，事件记录eri中涉及的UE的UE标识符可以作为加密信息提供，以便确保UE所有者的隐私。无论如何，如果需要出现，那么可以通过实现合适的解密算法来解密加密的信息(即，对应于UE标识符的UE的所有者的身份)，诸如例如在2012年3月6日，National Institute of Standards and Technology FIPS-180-4，“Secure HashStandard(SHS)”中描述的算法SHA256。
***100可以通过获取在移动通信网络105中生成并使用的数据的记录从移动通信网络105检索(和/或接收)与通用UE相关的事件记录eri。例如，在移动通信网络105是GSM网络的情况下，可以从移动通信网络105检索也被称为呼叫数据记录的计费数据记录(CDR)和/或访问者地点记录(VLR)，并重新用作事件记录eri。CDR是包含特定于在UE和移动通信网络105之间进行的电话呼叫或其他通信事务的单个实例的属性的数据记录(通常用于由通过移动通信网络105操作的移动电话服务提供商的计费目的)。VLR是列出已漫游到移动通信网络105的移动交换中心(MSC，未示出)的管辖范围内的UE的数据库，MSC是移动通信网络105的管理多个通信站105a上的事件的管理元件。移动通信网络105中的每个通信站105a通常与相应的VLR相关联。
相反，如果移动通信网络105是LTE网络，那么与通用UE的事件记录eri相关联的数据的记录由包含在移动通信网络105中的移动性管理实体(或MME)生成，MME负责LTE网络(其中没有实现VLR)中的UE跟踪和寻呼过程。
应当注意的是，本公开中描述的方法可以通过使用可以从中获得事件记录eri的(例如，由一个或多个WiFi网络提供的)任何数据源来实现，事件记录eri包括个人的单义标识符(诸如，上面提到的UE标识符)、这种个人的位置指示以及在其间已发生这种事件的时刻的时间指示。
在操作中，事件记录eri可以由***100从移动通信网络105连续地检索。可替代地，事件记录eri可以由***100例如针对预定的时间段(例如，每特定数量的小时，以每天或每周为基础)定期地收集。例如，事件记录eri可以在它们被生成时以一种“推送”模式从移动通信网络105被传递到***100，或者事件记录eri可以每天在移动通信网络105中被收集，并且然后定期地或根据***100的请求打包并传递到***100。
从移动通信网络105检索到的事件记录eri存储在储存库115中，其中使得它们对于计算引擎110可用以用于处理。优选地，由相同UE生成的事件记录eri在储存库115中被分组到一起，即，如果它们包括公共的UE标识符，那么事件记录eri被分组到一起并且在下文中被表示为事件记录组ergl(例如，l＝0，...，L，L≥0)。
优选地，计算引擎110根据由例如包含在图1的示例中的计算引擎110中的、存储在***100的存储器元件110a中的软件程序产品实现的(在下文中描述的)人群估计算法来处理事件记录eri(如下面详细讨论的)，即使该软件程序产品也可以存储在储存库115中(或者在设置在***100中的任何其他存储器元件中)。
甚至更优选地，事件记录eri根据由***管理员(通过管理员接口120)提供的(例如存储在储存库115中的)指令，并且可能根据由用户(通过用户界面125)提供的指令被处理，如在下面描述的。最后，计算引擎110通过用户界面125向用户提供对事件记录eri进行的处理的结果，并且可选地将这些处理结果存储在储存库115中。
现在转到图2A至图2E，它们是根据本发明的实施例的其中移动通信网络105的小区105b可以被建模的示例性形状。
为了本发明的目的，移动通信网络105的每个小区105b可以被建模为具有相应的小区中心B(不一定对应于通信站105a的地理位置)以及相应的小区半径Rc的区域(如图2A所示)，所述区域包围由对应的通信站105a服务的有效服务区域(未示出)(例如，其中每个点由通信站105a传送的无线电信号可到达的区域)。可替代地，小区半径Rc可以对应于包围有效服务区域的大部分的圆周的半径，诸如有效服务区域的85％或更多，诸如例如有效服务区域的90％。
应当注意的是，小区105b不限于盘状形状，事实上，小区105b可以具有优选地规则的多边形的形状。在这种情况下，小区中心B对应于多边形的质心(或形心(centroid))，而小区半径Rc对应于连接多边形质心(即，小区中心B)与多边形的顶点(如图2B和图2D中所示)或多边形的一侧的中点(如图2C和图2E中所示)的段。
有效服务区域，并且因此小区半径Rc可以借助由移动通信网络105的提供商所使用的众所周知的网络规划软件工具来定义，或者可以基于(全向或定向的，诸如具有120°辐射角的)天线辐射图和简单的辐射模型(诸如，例如在Theodore S.Rappaport，“WirelessCommunications”，Prentice Hall，1996年中所描述的辐射模型)来计算。
可替代地，可以借助Voronoi棋盘(tessellation)图对移动通信网络105进行建模，其中每个Voronoi小区对应于移动通信网络105的小区105b(由于Voronoi棋盘图在本领域中是众所周知的，因此本文不对它们进一步讨论)。
优选地，移动通信网络105的小区105b的列表和数量由管理员通过管理员界面120输入到***100。
在根据本发明的实施例的解决方案中，***100适于基于通过在AoI 107内的UE和服务该UE的移动通信网络105之间的交互生成的事件ei来识别个体是否参加了在AoI 107内进行的一个或多个公共集会。
现在转到图3A至图3E，它们是根据本发明的实施例的待确定的AoI 107可以采取的示例性形状。
通常，公共集会的AoI 107可以被建模为具有AoI中心C和AoI半径Ra的区域。例如，AoI 107可以由以AoI中心C为中心并且具有AoI半径Ra作为圆周半径的圆周界定(如图3A中所示)。
应当注意的是，AoI 107可以具有与圆周不同的形状。例如，AoI 107可以具有优选地规则的多边形的形状。在这种情况下，按照与上述小区105b建模类似的方式，AoI中心C对应于多边形的质心(或形心)，而AoI半径Ra对应于连接多边形的质心与多边形的顶点(如图3B和图3D中所示)或者多边形的一侧的中点(如图3C和图3E中所示)的段。
AoI中心C可以(例如，由用户通过用户界面125或由***管理员通过管理员界面120)被设置为AoI 107的(地理)中心点(例如，核心地点的地理中心点)、设置为公共集会的核心地点的地址、设置为通过地图软件(诸如web地图服务(例如，Google mapsTM，OpenStreetMapTM等)提供的点。
如将在下面更详细描述的，AoI半径Ra可以取零或负值以及正值。在AoI半径Ra取零或负值的情况下，AoI 107限于AoI中心C(即，一个或多个公共集会的核心地点)。AoI半径Ra的零值或负值的含义将通过参考下面描述的实施例中的这些零值或负值进一步阐明。
下面描述的算法被配置为确定AoI 107的AoI半径Ra的最佳半径值Ro。在本发明的一个实施例中，最佳半径值Ro借助从最小半径值Rmin开始到最大半径值Rmax的迭代步骤来确定(如下文中所述)。优选地，最小半径值Rmin和最大半径值Rmax由***100的管理员通过管理员界面120设置。
在本发明的实施例中，基于由申请人进行的关于多个过去的公共集会的经验数据的统计分析，最小半径值Rmin被设置为等于-1500m(Rmin＝-1500m)，而最大半径值Rmax被设置为等于1500m(Rmax＝1500m)。
在已经定义了移动通信网络105的小区105b的形状和AoI 107的形状之后，现在将介绍根据本发明的实施例的相关小区的概念，即，被认为至少部分地属于AoI 107的移动通信网络105的小区105b。
图4A至图4D是根据本发明的实施例的移动通信网络105相对于AoI 107的小区105b中的相关服务区域或小区405a-d。
在本发明的一个实施例中，假定AoI 107具有AoI中心C并且通用小区105b具有小区中心B和小区半径Rc，如果验证了以下不等式，那么通用小区105b可以被认为是AoI 107的相关小区405a-d：
Dist(C,B)≤|Rc+Ra|， (1)
其中Dist(C,B)是AoI中心C和小区中心B之间的地理距离。
根据AoI 107的AoI半径Ra的值，不等式(1)可以采取三个不同的含义。
即，如果AoI 107的AoI半径Ra大于零(即，Ra>0)，那么不等式(1)简化为：
Dist(C,B)≤(Rc+Ra), (2)
并且如果AoI 107和通用小区105b的区域至少部分地重叠(即使AoI中心C落在通用小区105b之外)，那么通用小区105b被认为是具有大于零的AoI半径Ra的AoI 107的相关小区(诸如图4A中的相关小区405a的情况)。
如果AoI 107的AoI半径Ra等于零(即，Ra＝0)，那么不等式(1)简化为：
Dist(C,B)≤Rc, (3)
并且如果AoI 107的AoI中心C被包含在通用小区105b中，那么通用小区105b被认为是具有等于零的AoI半径Ra的AoI 107的相关小区(诸如图4B和图4C中的相关小区405b和405c的情况)。
最后，如果AoI 107的AoI半径Ra小于零(即，Ra<0)，那么如果AoI 107的AoI中心C被包含在距离小区中心B等于或小于Rc-|Ra|的通用小区105b内，则通用小区105b被认为是具有AoI半径Ra小于零的AoI 107的相关小区(诸如，在图4D中的相关小区405d的情况)。
除了被保持在特定地点(即，AoI 107)之外，(通用)公共集会S具有开始时间Ts和结束时间Te。因此，为了本发明的目的，公共集会S具有等于观察时间间隔[Ts，Te]的相关持续时间(即，从开始时间Ts开始并在结束时间Te结束的时间间隔，持续Te-Ts时间单位，例如秒、分钟或小时)。
可以定义开始时间Ts和结束时间Te两者，以便对应于为该公共集会S计划的正式(正式宣布的)开始和结束时间。但是，申请人已经观察到，通过相对于公共集会S的正式开始时间来预期开始时间Ts，可以考虑人们(即，参加公共集会S的UE所有者)在公共集会S的正式开始时间之前到达AoI 107的事实，这对于收集关于参加者到达公共集会S的流动的时间上的趋势的数据可能有用。例如，基于先前的公共集会的经验数据，申请人已经发现将开始时间Ts预期为在公共集会S的正式开始时间之前的60分钟可能是有用的，以便考虑参加者到达公共集会S的趋势。
类似地，申请人已经观察到，结束时间Te可以相对于公共集会S的正式结束时间被推迟，以便考虑人们在公共集会的正式结束时间之后离开AoI 107的事实，这对于收集关于参加者离开公共集会S的流动的时间上的趋势的数据可能有用。例如，基于先前的公共集会的经验数据，申请人已经发现将结束时间Ts延迟到公共集会S的官方结束时间之后30分钟可能是有用的，以便考虑参加者离开公共集会S的趋势。
无论如何，管理员通过管理员界面120和/或用户通过用户界面125可以为公共集会S设置任何定制的开始时间Ts和结束时间Te。例如，开始时间Ts和结束时间Te可以被设置以便定义比公共集会S的有效持续时间短(比整个事件的持续时间短)的观察时间间隔[Ts，Te]，以便分析仅在公共集会S的整个持续时间的子部分期间到公共集会S的参加者的数量或变化。
另外，在公共集会的过程期间，管理员可以实时地改变结束时间Te，以便获得到公共集会S的参加者的即时数量。例如，管理员可以将结束时间Te设置为当前时刻以便确定到当前时刻为止到公共集会的参加者的数量。此外，管理员可以确定直到多个后续结束时间Te为止到公共集会S的参加者的数量，以便识别随着时间到公共集会S的参加者的数量上的趋势。
已经描述了***100，以及公共集会S的时间(即，开始时间Ts和结束时间Te)和空间(即，AoI 107的AoI中心C和AoI半径Ra)特性，现在将通过参考图5A至图5C来描述根据本发明的实施例的人群估计算法(或人群计数算法)，图5A至图5C是其示意性框图。
人群估计算法的第一部分被配置为基于关于所考虑的公共集会S的数据来确定AoI 107的AoI半径Ra的最佳半径值Ro。
初始地(步骤502)，公共集会S的AoI中心C、观察日期g和开始时间Ts和结束时间Te例如由用户通过用户界面125或由管理员通过管理员界面120被输入到***100。
之后(步骤504)，迭代变量k被初始化为零(即，k＝0)，并且计算出的半径值Rk初始地被设置为最小半径值Rmin(即，Rk＝Rmin)。迭代变量k说明算法的第一部分的迭代次数。
接下来(步骤506)，借助如上所述的不等式(1)来识别具有等于计算出的半径值Rk的AoI半径Ra(Ra＝Rk)的AoI 107的相关小区405a-d。
之后(步骤508)，通过组合在先前步骤506处识别出的相关小区405a-b的半径Rc来计算总半径Rtotk。例如，总半径被简单地定义为：
Rtotk＝∑rRcr， (4)
其中Rcr表示在步骤506处识别出的第r个相关小区405a-b(1≤r≤R，R是正整数)的半径。
从储存库115中检索参考观察日期g在观察时间间隔[Ts，Te]期间和在步骤506处确定的相关小区405a-d内已发生的所有事件记录eri(步骤510)。
随后(步骤512)，将第一UE数量Uk计算为对应于(即，与之相关联的)参考已在先前步骤506被检索到的相关小区405a-d的事件记录eri中的至少一个事件记录eri的UE的数量(第一UE数量Uk取决于相关小区，并且因此，取决于计算出的半径值Rk)。
类似地，从储存库115中检索参考在观察日期g之前的先前日期gp在观察时间间隔[Ts，Te]期间并且在步骤506处确定的相关小区405a-d内已发生的所有事件记录eri’(步骤514)。
在本发明的一个实施例中，对于公共集会S，考虑了在观察日期g之前的P个先前日期gp(其中1≤p≤P并且P是整数)的集合。所考虑的先前日期gp的数量P优选地由管理员(通过管理员界面120)设置。在本发明的一个实施例中，管理员根据储存库115的存储能力(即，为了能够存储关于P个先前日期gp的所有数据)和/或基于计算引擎110的计算能力(即，为了能够处理关于P个先前日期gp的所有数据)设置先前日期gp的数量P。优选地，管理员还基于对同类(即，文化、娱乐、政治或体育公共集会)的过去公共集会的统计分析来设置先前日期gp的数量。
申请人已经发现，将先前日期gp的数量P设置为等于6(即，P＝6)为大多种类的公共集会提供了良好的结果(但是这不应当被解释为对于本发明是限制性的)。
然后(步骤516)，对于先前日期gp中的每一个，将第二UE数量Upk计算为对应于参考已在先前步骤506处被检索到的相关小区405a-d的事件记录eri’中的至少一个事件记录eri’的UE的数量(第二UE数量Upk取决于相关小区，并且因此取决于计算出的半径值Rk)。
刚计算出的第二UE数量Upk被组合(步骤518)，以便确定考虑在P个先前日期gp期间在相关小区内的UE数量的平均UE数量μk(其中)以及UE数量标准偏差σk(其中
平均UE数量μk和UE数量标准偏差σk与第一UE数量Uk被组合(步骤520)，以便获得定义为z-分值Zk的(统计)量(其取决于计算出的半径值Rk)：
Zk＝(Uk–μk)/σk, (5)
刚计算出的z-分值Zk相对于在步骤508处计算出的总半径Rtotk被归一化(步骤522)：
Z’k＝Zk/Rtotk, (6)
接下来，变量k被增加单位元素(步骤524；即，k＝k+1)并且计算出的半径值Rk被增加(步骤526)：
Rk＝Rmin+kΔ， (7)
其中Δ是可以由管理员定义的迭代步骤或宽度(例如，Δ＝100m)，因此每个计算出的半径值Rk与下一个计算出的半径值分开迭代宽度Δ。应当注意的是，迭代宽度Δ定义迭代变量k的最大迭代值kmax，-并且因此，用于确定最佳半径值Ro的最大迭代次数-为：
kmax＝(|Rmin|+Rmax)/Δ， (8)
应该注意的是，***管理员可以使用迭代宽度Δ来调整利用其来确定最佳半径值Ro的粒度(即，细度)，即，由管理员设置的迭代宽度Δ越小，由最大迭代值kmax定义的迭代数量越大，并且因此，算法的粒度越细。
在本发明的实施例中，由于最小半径值Rmin被设置为-1500m，因此最大半径值Rmax被设置为1500m，并且迭代宽度Δ被设置为100m，对于迭代变量k的最大迭代值kmax结果等于30，并且因此，用于确定最佳半径值Ro的最大迭代次数被限制为30。
之后，检查(步骤528)计算出的半径值Rk是否低于或等于最大半径值Rmax：
Rk≤Rmax， (9)
在肯定的情况(判定框528的退出分支Y)下，即，计算出的半径值Rk低于或等于最大半径值Rmax(即，Rk≤Rmax)，操作返回到步骤506，以便基于(在步骤526)刚被增加另外的第k个迭代宽度Δ的计算出的半径值Rk开始算法的第一部分的新的迭代。
在否定的情况(判定框528的退出分支N)下，即，计算出的半径值Rk大于最大半径值Rmax(即，Rk>Rmax)，最佳半径值Ro被计算(步骤530)为由归一化z-分值Z’k加权的计算出的半径值Rk(其中1≤k≤kmax)的平均值，所述平均值按以下公式计算，或者：
迭代算法的第一部分的步骤506至528，直到计算出的半径值Rk大于最大半径值Rmax(即，Rk>Rmax)为止，并且(在步骤530)计算最佳半径值Ro。
随着在步骤530处的最佳半径值Ro的计算，算法的第一部分结束，并且算法的第二部分(在下面描述的步骤532处)开始。在算法的第一部分结束时，AoI 107由AoI中心C和由被设置为等于最佳半径值Ro的AoI半径Ra(Ra＝Ro)正确地定义。
根据本发明的实施例的算法的第二部分被配置为确定正在参加公共集会S的人数。
在步骤530处已经计算出最佳半径值Ro之后，定义实际相关小区405a-d的集合(步骤532)。该集合包括当AoI半径Ra被设置为等于最佳半径值Ro时，要对其验证不等式(1)的移动通信网络105的所有小区105b，或者：
Dist(C,B)≤|Rc+Ro| (11)
然后(步骤534)，从储存库115中检索参考观察日期g在观察时间间隔[Ts，Te]期间并且在步骤532处确定的实际相关小区405a-d内已经发生的新的和所有事件记录eri。
随后(步骤536)，构建了UE列表uL。UE列表uL包括与参考在先前步骤534处已经检索到的相关小区的事件记录eri中的至少一个事件记录eri对应的每个UE的标识符。
UE列表uL提供在公共集会S处的可能的参加者(即，包含在UE列表uL中的UE的所有者)的第一估计。因此，在本发明的简化的实施例中，算法可以简单地计算公共集会S的UE列表uL，并且通过用户界面125向用户提供结果产生的UE列表uL。
一旦UE列表uL已被建立，UE变量j就被初始化为单位元素(即，j＝1)，并且人数A被初始化为零(即，A＝0)(步骤538)。UE变量j用于扫描包含在UE列表uL中的所有用户，而人数A说明在公共集会S处的参加者的数量(如下所述)。
从储存库115中检索参考在P个先前日期gp中的每一个中在观察时间间隔[Ts，Te]期间记录的UE UEj并且在移动通信网络105的小区105b中的任何一个小区内已发生的所有事件记录eri’(步骤540)。
然后(步骤542)，为UE UEj计算连续事件记录eri之间的平均中间到达时间iat。在本发明的一个实施例中，UE UEj的中间到达时间被计算为两个连续事件记录eri的时间数据(即，时间戳)之间的差。优选地，平均中间到达时间iat基于在步骤534处检索到的观察日期g期间记录的事件记录eri和eri’以及在步骤540处检索到的P个先前日期gp来计算。
为UE UEj识别参考观察日期g在观察时间间隔[Ts，Te]期间并且在步骤532处确定的实际相关小区405a-d内已发生的第一事件记录er1和最后事件记录erI(步骤544)，并且从其中检索相应的第一观察时间数据tdf和最后观察时间数据tdI(步骤546)。
第一观察时间数据tdf、最后观察时间数据tdI和平均中间到达时间iat被组合(步骤548)，以便确定UE UEj在观察日期g在观察时间间隔[Ts，Te]期间在AoI 107内花费的第一时间片段f1：
随后，为UE UEj识别参考P个先前观察日期gp在观察时间间隔[Ts，Te]期间并且在步骤532处确定的实际相关小区405a-d内已发生的所有事件记录eri’中的第一事件记录erf’和最后事件记录erI’(步骤550)，并且从其中检索相应的第一先前时间数据tdp1和最后先前时间数据tdpI(步骤552)。
第一先前时间数据tdpf、最后先前时间数据tdpI和平均中间到达时间iat被组合(步骤554)，以便确定UE UEj在P个先前日期gp期间在AoI 107内花费的第二时间片段f2：
其中Tgp是P个先前日期gp的总持续时间，其可以根据用于算法中的时间量(诸如，例如，第一先前时间数据tdpf、最后先前时间数据tdpI和平均中间到达时间iat)的时间单位(例如，秒、分钟或小时：例如，如果所考虑的时间单位是一小时，那么Tgp＝P×24)，例如以秒、分钟或小时来计算。
之后(步骤556)，通过组合第一时间片段f1和第二时间片段f2来计算参加公共集会S的UE UEj的所有者的人概率pj：
pj＝f1*(1-f2)， (14)
因此，第一时间片段f1和第二时间片段f2可以被认为是概率。即，第一时间片段f1可以被解释为UE UEj的所有者在公共集会S期间在AoI 107中的概率，而第二时间片段f2可以被解释为UE UEj的所有者在先前日期gp期间在AoI 107中的概率。
然后通过将人概率pj添加到人数A的实际值来更新人数A(步骤558)，或者：
A＝A+pj (15)
应当注意的是，根据本发明的人数A被计算为参考在UE列表uL中列出的每个UEUEj的每个相应所有者的人概率pj的总和。
UE变量j被增加单位元素(步骤560；即，j＝j+1)，并且检查(步骤562)UE变量j是否低于或等于在UE列表uL中列出的所列出的用户J的总数(其中J是整数)：
j≤J (16)
在肯定的情况(判定框562的退出分支Y)下，即，UE变量j小于或等于列出的用户J的数量J(j≤J)，那么操作返回到步骤540，以便分析参考下一个UE UEj的事件记录eri。
在否定的情况(判定框562的退出分支N)下，即，UE变量j大于列出的用户的总数J(j>J)，并且UE列表uL已被完全扫描。因此，参考在观察日期g举行的公共集会S的人数A被存储(步骤564)在储存库115中。
优选地，通过用户终端125向用户供给(步骤566)用于检查和/或进一步处理的结果，即，在步骤536处确定的人数A和可能的UE列表uL以及UE列表uLn中的每个UE的相应的第一和第二时间片段f1和f2，人群估计算法终止。
迭代算法的第二部分的步骤540至562，直到UE列表uL中的所有UE UEj都被分析为止，并且因此(在步骤566处)算法终止，其中通过用户终端125向用户供给结果。
总之，人群估计算法(或人群计数算法)包括第一部分和第二部分。
人群估计算法的第一部分包括扫描(步骤506-528)公共集会的最小半径值Rmin和最大半径值Rmax之间的所有计算出的半径值Rk的第一循环。对于每个计算出的半径值Rk，确定相应的相关小区405a-d和归一化z-分值Z’k。基于这种数据(即，相应的相关小区405a-d和归一化z-分值Z’k)，识别最佳半径值Ro。在人群估计算法的第一部分的结束处，定义具有最佳半径值Ro的AoI 107。
人群估计算法的第二部分包括扫描(步骤540-562)在观察日期g中在观察时间间隔[Ts，Te]期间在至少一个相关小区405a-d(即，AoI 107)中生成事件记录eri的所有UEUEj。对于每个UE UEj，确定在公共集会S的观察日期g在观察时间间隔[Ts，Te]期间在包含在具有等于最佳半径值Ro的AoI半径Ra的AoI 107内的相关小区405a-d中的一个或多个相关小区内花费的时间(即，第一时间片段f1)，以及在P个先前日期gp期间在相同AoI 107内花费的另外的时间(即，第二时间片段f2)。基于在AoI 107中花费的这些时间的知识，确定UE UEj的所有者正在参加公共集会S的人概率pj，并且基于该人概率pj，确定正在参加公共集会S的人数。
根据本发明的实施例的人群估计***100和人群估计算法允许以可靠的方式对在公共集会S处的参加者的数量进行实时计算，并且(通过确定最佳半径值Ro)正确地识别与公共集会S相关联的AoI107的有效扩展。

Claims (19)

1.一种估计为了参加在一个日期的时间间隔期间的公共集会而聚集在感兴趣区域(107)处的人数的方法，其中所述感兴趣区域(107)由感兴趣区域中心和感兴趣区域半径来定义，并且由具有多个通信站(105a)的移动电信网络(105)覆盖，每个通信站适于管理其中细分移动电信网络(105)的一个或多个服务区域(105b)中的用户设备的通信，所述方法包括以下步骤：
a)定义(504、506、524-528)感兴趣区域半径的多个计算出的半径值，以及
对于每个计算出的半径值：
b)识别(506)包含在感兴趣区域(107)中的服务区域(105b)中的相关服务区域(405a-d)的数量；
c)识别(510、512)与在感兴趣区域(107)内的公共集会的所述日期的时间间隔期间在用户设备和移动通信网络(105)之间发生的交互的对应事件的至少一个事件记录相关联的用户设备的第一数量；
d)识别(514、516)在感兴趣区域(107)内的公共集会的所述日期之前的预定数量的先前日期中的每一个先前日期的时间间隔期间在用户设备和移动通信网络(105)之间发生的交互的对应事件的至少一个事件记录相关联的用户设备的第二数量；
e)将用户设备的第一数量和用户设备的第二数量组合(518、520)以获得统计量；
f)通过相对于相关服务区域(405a-d)的半径归一化(522)统计量来计算归一化统计量；
g)将感兴趣区域半径的最佳半径值计算(530)为由归一化统计量加权的计算出的半径值的平均值；
h)计数(532-566)为了参加在具有等于最佳半径值的感兴趣区域半径的感兴趣区域内的公共集会而聚集的人数。
2.如权利要求1所述的方法，其中如果服务区域(105b)验证了以下不等式，那么服务区域(105b)被识别为相关服务区域(405a-d)：
Dist(C,B)≤|Rc+Rk|，
其中C是感兴趣区域(107)的中心，B是服务区域(105b)的中心，Dist(C,B)是感兴趣区域的中心C与服务区域的中心B之间的地理距离，Rc是服务区域的半径，并且Rk是计算出的半径值。
3.如权利要求1或2所述的方法，其中组合(518、520)用户设备的第一数量和用户设备的第二数量以便获得统计量的步骤d)包括：
-组合(518)先前日期中的每一个先前日期的用户设备的第二数量，以便确定平均用户设备数量和用户设备数量标准偏差。
4.如权利要求3所述的方法，其中组合(518、520)用户设备的第一数量和用户设备的第二数量以便获得统计量的步骤d)还包括：
-将统计量计算为：
Zk＝(Uk–μk)/σk，
其中Uk是第一数量，μk是平均用户设备数量并且σk是用户设备数量标准偏差。
5.如前述权利要求中任何一项所述的方法，其中所述多个计算出的半径值的范围从最小半径值到最大半径值，每个计算出的半径值与下一个半径值分开迭代宽度。
6.如前述权利要求中任何一项所述的方法，其中，计数(532-566)为了参加在具有等于最佳半径值的感兴趣区域半径的感兴趣区域内的公共集会而聚集的人数的步骤h)包括：
i)定义(532)包含在移动通信网络(105)中的服务区域(105b)中的相关服务区域(405a-d)的数量，其中所述相关服务区域是至少部分地叠加在具有等于最佳半径值的感兴趣区域半径的感兴趣区域(107)上的服务区域。
7.如权利要求6所述的方法，其中，如果服务区域(105b)验证了以下不等式，那么服务区域(105b)被识别为相关服务区域(405a-d)：
Dist(C,B)≤|Rc+Ro|，
其中C是感兴趣区域(107)的中心，B是服务区域(105b)的中心，Dist(C,B)是感兴趣区域的中心C与服务区域的中心B之间的地理距离，Rc是服务区域的半径，并且Ro是最佳半径值。
8.如前述权利要求1至7中任何一项所述的方法，其中，计数(532-566)为了参加在具有等于最佳半径值的感兴趣区域半径的感兴趣区域内的公共集会而聚集的人数的步骤h)还包括：
j)构建(534、536)用户设备列表，所述用户设备列表包括与在感兴趣区域(107)内在所述日期的时间间隔中发生的事件记录中的至少一个事件记录相关联的每个用户设备的标识符。
9.如权利要求8所述的方法，其中，计数(532-566)为了参加在具有等于最佳半径值的感兴趣区域半径的感兴趣区域内的公共集会而聚集的人数的步骤h)还包括，对于用户设备列表中的每个用户设备：
k)计算(542)与用户设备相关联的连续事件记录之间的平均中间到达时间。
10.如权利要求9所述的方法，其中所述平均中间到达时间基于在所述日期和在所述先前日期两者发生的事件记录来计算。
11.如权利要求9或10所述的方法，其中，计数(532-566)为了参加在具有等于最佳半径值的感兴趣区域半径的感兴趣区域内的公共集会而聚集的人数的步骤h)还包括，对于用户设备列表中的每个用户设备：
l)识别(544、546)在感兴趣区域(107)内在所述日期的所述时间间隔中发生的、分别参考第一事件记录和最后事件记录的第一时间数据和最后时间数据，以及
m)组合(548)第一时间数据、最后时间数据和平均中间到达时间，以便确定第一时间片段，所述第一时间片段指示用户设备在所述日期在所述时间间隔期间在感兴趣区域(107)内花费的时间段。
12.如权利要求11所述的方法，其中所述第一时间片段被确定为：
<mrow> <mi>f</mi> <mn>1</mn> <mo>=</mo> <mfrac> <mrow> <mo>|</mo> <mrow> <msub> <mi>td</mi> <mi>I</mi> </msub> <mo>-</mo> <msub> <mi>td</mi> <mi>f</mi> </msub> <mo>+</mo> <mi>i</mi> <mi>a</mi> <mi>t</mi> </mrow> <mo>|</mo> </mrow> <mrow> <mi>T</mi> <mi>e</mi> <mo>-</mo> <mi>T</mi> <mi>s</mi> </mrow> </mfrac> </mrow>
其中tdf是第一时间数据，tdI是最后时间数据，iat是平均中间到达时间，Ts是公共集会的开始时间并且Te是公共集会的结束时间。
13.如前述权利要求10至12中任何一项所述的方法，其中，计数(532-566)为了参加在具有等于最佳半径值的感兴趣区域半径的感兴趣区域内的公共集会而聚集的人数的步骤h)还包括，对于UE列表中的每个用户设备：
n)识别(550、552)在感兴趣区域(107)内在先前日期的时间间隔中发生的、分别参考第一事件记录和最后事件记录的第一先前时间数据和最后先前时间数据，以及
o)组合(554)第一先前时间数据、最后先前时间数据和平均中间到达时间，以便确定第二时间片段，所述第二时间片段指示用户设备在先前日期期间在感兴趣区域(107)内花费的时间段。
14.如权利要求13所述的方法，其中所述第二时间片段被确定为：
<mrow> <mi>f</mi> <mn>2</mn> <mo>=</mo> <mfrac> <mrow> <mo>|</mo> <mrow> <msub> <mi>tdp</mi> <mi>I</mi> </msub> <mo>-</mo> <msub> <mi>tdp</mi> <mi>f</mi> </msub> <mo>+</mo> <mi>i</mi> <mi>a</mi> <mi>t</mi> </mrow> <mo>|</mo> </mrow> <mrow> <mi>T</mi> <mi>g</mi> <mi>p</mi> </mrow> </mfrac> </mrow>
其中tdpf是第一先前时间数据，tdpI是最后先前时间数据，iat是平均中间到达时间，并且Tgp是先前日期的总持续时间。
15.如权利要求13或14所述的方法，如取决于权利要求11或12，其中，计数(532-566)为了参加在具有等于最佳半径值的感兴趣区域半径的感兴趣区域内的公共集会而聚集的人数的步骤h)还包括，对于用户设备列表中的每个用户设备：
p)组合(554)第一时间片段和第二时间片段，以便确定用户设备的所有者参加公共集会的概率。
16.如权利要求15所述的方法，其中所述用户设备的所有者参加公共集会的概率被计算为：
pj＝f1*(1–f2)，
其中f1是第一时间片段，并且f2是第二时间片段。
17.如权利要求15或16所述的方法，其中，计数(532-566)为了参加在具有等于最佳半径值的感兴趣区域半径的感兴趣区域内的公共集会而聚集的人数的步骤h)还包括：
q)将为了参加公共集会而聚集的人数计算(556)为针对用户设备列表的每个用户设备确定的概率的总和。
18.一种耦合到无线电信网络(105)的用于估计在感兴趣区域(107)处聚集的人数的***(100)，所述***包括：
计算引擎(110)，所述计算引擎(110)适于处理从移动电话网络(105)检索到的数据；
储存库(115)，所述储存库(115)适于存储关于用户设备和移动电话网络之间的交互的数据、由计算引擎生成的计算结果以及可能的由***生成的和/或提供给***的任何处理数据，以及
管理员界面(120)，所述管理员界面(120)能够操作用于修改由计算引擎使用的参数和/或算法，和/或访问存储在储存库中的数据，
其特征在于
还包括存储器元件(110a)，所述存储器元件(110a)存储被配置为通过所述***(100)实现如权利要求1至17中任何一项所述的方法的、用于估计在感兴趣区域处聚集的人数的软件程序产品。
19.如权利要求18所述的***，还包括适于接收来自所述***的用户的输入和向所述***的用户提供输出的至少一个用户界面(125)，所述用户包括一个或多个人和/或由所述***提供的服务的一个或多个外部计算***订户。
CN201480084457.3A 2014-12-30 2014-12-30 用于实时计数公共集会处的参与者的数量的方法和*** Active CN107113560B (zh)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/EP2014/079435 WO2016107647A1 (en) 2014-12-30 2014-12-30 Method and system for a real-time counting of a number of participants at a public happening

Publications (2)

Publication Number Publication Date
CN107113560A true CN107113560A (zh) 2017-08-29
CN107113560B CN107113560B (zh) 2020-03-20

Family

ID=52432769

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201480084457.3A Active CN107113560B (zh) 2014-12-30 2014-12-30 用于实时计数公共集会处的参与者的数量的方法和***

Country Status (5)

Country Link
US (1) US10285002B2 (zh)
EP (1) EP3241368B1 (zh)
CN (1) CN107113560B (zh)
BR (1) BR112017014024B1 (zh)
WO (1) WO2016107647A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107624249A (zh) * 2015-03-30 2018-01-23 意大利电信股份公司 用于通过电信网络的聚合数据实时计数人群中的人数的方法和***
CN110781743A (zh) * 2019-09-23 2020-02-11 腾讯科技（深圳）有限公司 一种人群聚集事件的处理方法、装置及电子设备

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9179185B2 (en) 2013-03-14 2015-11-03 The Nielsen Company (Us), Llc Methods and apparatus to determine a number of people in an area
US9396354B1 (en) 2014-05-28 2016-07-19 Snapchat, Inc. Apparatus and method for automated privacy protection in distributed images
US9113301B1 (en) 2014-06-13 2015-08-18 Snapchat, Inc. Geo-location based event gallery
US10824654B2 (en) 2014-09-18 2020-11-03 Snap Inc. Geolocation-based pictographs
US11216869B2 (en) 2014-09-23 2022-01-04 Snap Inc. User interface to augment an image using geolocation
US9015285B1 (en) 2014-11-12 2015-04-21 Snapchat, Inc. User interface for accessing media at a geographic location
US9385983B1 (en) 2014-12-19 2016-07-05 Snapchat, Inc. Gallery of messages from individuals with a shared interest
US10311916B2 (en) 2014-12-19 2019-06-04 Snap Inc. Gallery of videos set to an audio time line
WO2016149594A1 (en) 2015-03-18 2016-09-22 Allen Nicholas Richard Geo-fence authorization provisioning
US10135949B1 (en) 2015-05-05 2018-11-20 Snap Inc. Systems and methods for story and sub-story navigation
WO2017025134A1 (en) * 2015-08-11 2017-02-16 Telecom Italia S.P.A. Method and system for computing an o-d matrix obtained through radio mobile network data
US10354425B2 (en) 2015-12-18 2019-07-16 Snap Inc. Method and system for providing context relevant media augmentation
US10771924B2 (en) 2016-12-22 2020-09-08 Telecom Italia S.P.A. Method and system for counting people at a transport hub by means of data of a telecommunication network
US10915911B2 (en) * 2017-02-03 2021-02-09 Snap Inc. System to determine a price-schedule to distribute media content
US10582277B2 (en) 2017-03-27 2020-03-03 Snap Inc. Generating a stitched data stream
IT201900015096A1 (it) 2019-08-27 2021-02-27 Telecom Italia Spa Sistema di rete auto-organizzante

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101751553A (zh) * 2008-12-03 2010-06-23 中国科学院自动化研究所 一种大规模人群密度分析和预测方法
CN102184512A (zh) * 2011-05-10 2011-09-14 天津大学 利用手机数据发现城市活动中异常事件的方法
CN102655654A (zh) * 2011-03-04 2012-09-05 中兴通讯股份有限公司 一种话务统计方法和装置、组网方法和装置
WO2013090747A1 (en) * 2011-12-16 2013-06-20 Marvell World Trade Ltd. Current balancing circuits for light-emitting-diode-based illumination systems
US20130226543A1 (en) * 2010-06-01 2013-08-29 At&T Intellectual Property I, L.P. Systems, Methods, and Computer Program Products for Estimating Crowd Sizes Using Information Collected From Mobile Devices in a Wireless Communications Network
CN103634775A (zh) * 2012-08-24 2014-03-12 中兴通讯股份有限公司 一种获得行政区域间用户流动的方法及装置
CN104244307A (zh) * 2014-09-18 2014-12-24 大唐移动通信设备有限公司 异常事件上报、处理方法、装置、基站及管理服务器

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130225196A1 (en) * 2012-02-27 2013-08-29 Bryan J. James Geo-Fence Generation and Updating Based on Device Movement Patterns

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101751553A (zh) * 2008-12-03 2010-06-23 中国科学院自动化研究所 一种大规模人群密度分析和预测方法
US20130226543A1 (en) * 2010-06-01 2013-08-29 At&T Intellectual Property I, L.P. Systems, Methods, and Computer Program Products for Estimating Crowd Sizes Using Information Collected From Mobile Devices in a Wireless Communications Network
CN102655654A (zh) * 2011-03-04 2012-09-05 中兴通讯股份有限公司 一种话务统计方法和装置、组网方法和装置
CN102184512A (zh) * 2011-05-10 2011-09-14 天津大学 利用手机数据发现城市活动中异常事件的方法
WO2013090747A1 (en) * 2011-12-16 2013-06-20 Marvell World Trade Ltd. Current balancing circuits for light-emitting-diode-based illumination systems
CN103634775A (zh) * 2012-08-24 2014-03-12 中兴通讯股份有限公司 一种获得行政区域间用户流动的方法及装置
CN104244307A (zh) * 2014-09-18 2014-12-24 大唐移动通信设备有限公司 异常事件上报、处理方法、装置、基站及管理服务器

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107624249A (zh) * 2015-03-30 2018-01-23 意大利电信股份公司 用于通过电信网络的聚合数据实时计数人群中的人数的方法和***
CN107624249B (zh) * 2015-03-30 2020-07-03 意大利电信股份公司 通过电信网络的聚合数据实时计数聚集人数的方法和***
CN110781743A (zh) * 2019-09-23 2020-02-11 腾讯科技（深圳）有限公司 一种人群聚集事件的处理方法、装置及电子设备

Also Published As

Publication number Publication date
US20170359686A1 (en) 2017-12-14
EP3241368B1 (en) 2018-10-31
WO2016107647A1 (en) 2016-07-07
BR112017014024A2 (pt) 2018-01-02
CN107113560B (zh) 2020-03-20
BR112017014024B1 (pt) 2023-01-17
US10285002B2 (en) 2019-05-07
EP3241368A1 (en) 2017-11-08

Similar Documents

Publication Publication Date Title
CN107113560A (zh) 用于实时计数公共集会处的参与者的数量的方法和***
CN107113561A (zh) 用于估计人群中的人数的方法和***
US10375516B2 (en) Method and system for a real-time counting of a number of persons in a crowd by means of aggregated data of a telecommunication network
Noulas et al. Exploiting foursquare and cellular data to infer user activity in urban environments
US9706363B2 (en) Method and system for identifying significant locations through data obtainable from a telecommunication network
EP3278579B1 (en) Method and system for estimating a posteriori a number of persons in one or more crowds by means of aggregated data of a telecommunication network
US10278015B2 (en) Method and system for real-time computing of origin-destination matrices relating to attendees at a public happening through analysis of mobile communication network data
US9955299B2 (en) Method and system for a posteriori computation of origin-destination matrices relating to gathering of people through analysis of mobile communication network data
US10887728B2 (en) Method and system for determining the mobility status of a user of a wireless communication network
Nguyen et al. Crowdsensing in opportunistic mobile social networks: A context-aware and human-centric approach
Keleher et al. Connecting isolated communities: Quantitative evidence on the adoption of community cellular networks in the Philippines
Bartolacci et al. Promoting Resiliency in Emergency Communication Networks: A Network Interdiction Modeling Approach
BR112017020838B1 (pt) Método para estimar um número de pessoas reunidas em uma área de interesse, e, sistema acoplado com uma rede de telecomunicação sem fio para estimar um número de pessoas reunidas em uma área de interesse

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant