CN111093156B - 伪基站的位置定位方法、设备和存储介质 - Google Patents
伪基站的位置定位方法、设备和存储介质 Download PDFInfo
- Publication number
- CN111093156B CN111093156B CN201911197361.9A CN201911197361A CN111093156B CN 111093156 B CN111093156 B CN 111093156B CN 201911197361 A CN201911197361 A CN 201911197361A CN 111093156 B CN111093156 B CN 111093156B
- Authority
- CN
- China
- Prior art keywords
- track
- sub
- points
- point
- terminal
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000005259 measurement Methods 0.000 claims abstract description 48
- 230000001133 acceleration Effects 0.000 claims description 22
- 238000007781 pre-processing Methods 0.000 claims description 16
- 238000012545 processing Methods 0.000 claims description 9
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 5
- 230000002159 abnormal effect Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 238000004891 communication Methods 0.000 description 5
- 238000007726 management method Methods 0.000 description 5
- 238000010295 mobile communication Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000012216 screening Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- HEEACTTWORLLPM-UHFFFAOYSA-N 2-(1h-imidazol-5-yl)ethanol Chemical compound OCCC1=CNC=N1 HEEACTTWORLLPM-UHFFFAOYSA-N 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Images
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/029—Location-based management or tracking services
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W12/00—Security arrangements; Authentication; Protecting privacy or anonymity
- H04W12/12—Detection or prevention of fraud
- H04W12/121—Wireless intrusion detection systems [WIDS]; Wireless intrusion prevention systems [WIPS]
- H04W12/122—Counter-measures against attacks; Protection against rogue devices
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Computer Security & Cryptography (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本申请提供一种伪基站的位置定位方法、设备和存储介质,该方法包括:获取至少一个终端周期性上报的测量报告;所述测量报告包括:所述终端的标识和位置;根据各个所述终端上报的测量报告,生成各个所述终端的至少两条轨迹记录;针对任一所述终端,从至少两条所述轨迹记录中提取满足预设的连续性条件的第一子轨迹,并从任意相邻的两条所述轨迹记录中的所述第一子轨迹中提取作为断点的轨迹点;对所述作为断点的轨迹点进行聚类处理,得到至少一个聚类簇;所述聚类簇包括至少两个轨迹点;将各个所述聚类簇中的轨迹点的位置进行加权处理,得到各个所述聚类簇中伪基站的位置。本申请实施例的方法伪基站的定位准确性和效率较高。
Description
技术领域
本申请涉及通信技术领域,尤其涉及一种伪基站的位置定位方法、设备和存储介质。
背景技术
随着移动通信技术的发展,蜂窝网络制式的升级换代,使移动通信的用户数量也越来越多。相应的,吸引部分不法分子利用伪基站吸纳移动互联网中的正常移动用户,强行向移动用户发送大量诈骗、推销等垃圾短信,给人们生活造成严重危害。
伪基站是一种利用全球移动通信(Global System for Mobile communication,GSM)单向认证缺陷的非法无线电通信设备,主要由主机和笔记本电脑组成,能够搜取以伪基站为中心、一定半径范围内的GSM移动电话信息。只要伪基站运行中,移动用户的信号就会被强制连接到该设备上,无法连接到公用电信网络,影响了移动用户的正常使用。目前伪基站的发现主要是基于人工路测,利用扫频仪、测试终端、移动终端等,在可能存在伪基站的区域进行探测,通过扫频仪测量接收的频率信息,进而定位伪基站的位置。
由于伪基站位置的隐蔽性,且运行时间存在随机性,人工路测的方式需要测量车辆大范围长时间的进行探测,效率较低,且部分区域因为测试车辆不可到达,会导致伪基站的位置无法确定。
发明内容
本申请提供一种伪基站的位置定位方法、设备和存储介质,以提高伪基站的位置定位的效率和准确性。
第一方面,本申请提供一种伪基站的位置定位方法,包括:
获取至少一个终端周期性上报的测量报告;所述测量报告包括:所述终端的标识和位置;
根据各个所述终端上报的测量报告,生成各个所述终端的至少两条轨迹记录;
针对任一所述终端,从至少两条所述轨迹记录中提取满足预设的连续性条件的第一子轨迹,并从任意相邻的两条所述轨迹记录中的所述第一子轨迹中提取作为断点的轨迹点;
对所述作为断点的轨迹点进行聚类处理,得到至少一个聚类簇;所述聚类簇包括至少两个轨迹点;
将各个所述聚类簇中的轨迹点的位置进行加权处理,得到各个所述聚类簇中伪基站的位置。
第二方面,本申请提供一种伪基站的位置定位装置,包括:
获取模块,用于获取至少一个终端周期性上报的测量报告;所述测量报告包括:所述终端的标识和位置;
预处理模块,用于根据各个所述终端上报的测量报告,生成各个所述终端的至少两条轨迹记录;
所述预处理模块,还用于针对任一所述终端,从至少两条所述轨迹记录中提取满足预设的连续性条件的第一子轨迹,并从任意相邻的两条所述轨迹记录中的所述第一子轨迹中提取作为断点的轨迹点;
所述预处理模块,还用于对所述作为断点的轨迹点进行聚类处理,得到至少一个聚类簇;所述聚类簇包括至少两个轨迹点;
处理模块,用于将各个所述聚类簇中的轨迹点的位置进行加权处理,得到各个所述聚类簇中伪基站的位置。
第三方面,本申请实施例提供一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现第一方面中任一项所述的方法。
第四方面,本申请实施例提供一种电子设备,包括:
处理器;以及
存储器,用于存储所述处理器的可执行指令;
其中,所述处理器配置为经由执行所述可执行指令来执行第一方面中任一项所述的方法。
本申请实施例提供的伪基站的位置定位方法、设备和存储介质,获取至少一个终端周期性上报的测量报告;所述测量报告包括:所述终端的标识和位置;根据各个所述终端上报的测量报告,生成各个所述终端的至少两条轨迹记录;针对任一所述终端,从至少两条所述轨迹记录中提取满足预设的连续性条件的第一子轨迹,并从任意相邻的两条所述轨迹记录中的所述第一子轨迹中提取作为断点的轨迹点;对所述作为断点的轨迹点进行聚类处理,得到至少一个聚类簇;所述聚类簇包括至少两个轨迹点;将各个所述聚类簇中的轨迹点的位置进行加权处理,得到各个所述聚类簇中伪基站的位置,上述通过终端上报的测量报告生成轨迹记录,进而根据轨迹记录中出现的断点确定伪基站的位置,准确性和效率较高。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本公开的实施例,并与说明书一起用于解释本公开的原理。
图1是本申请提供的伪基站的位置定位方法一实施例的流程示意图;
图2是本申请提供的方法一实施例的子轨迹示意图;
图3是本申请提供的方法一实施例的异常轨迹点修正示意图;
图4是本申请提供的方法一实施例的第一子轨迹提取原理示意图;
图5是本申请提供的方法一实施例的断点轨迹对生成原理示意图;
图6是本申请提供的伪基站的位置定位装置一实施例的结构示意图;
图7是本申请提供的电子设备实施例的结构示意图。
通过上述附图,已示出本公开明确的实施例,后文中将有更详细的描述。这些附图和文字描述并不是为了通过任何方式限制本公开构思的范围,而是通过参考特定实施例为本领域技术人员说明本公开的概念。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本公开相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本公开的一些方面相一致的装置和方法的例子。
本申请的说明书和权利要求书及所述附图中的术语“包括”和“具有”以及它们任何变形,意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、***、产品或设备没有限定于已列出的步骤或单元,而是可选地还包括没有列出的步骤或单元,或可选地还包括对于这些过程、方法、产品或设备固有的其它步骤或单元。
首先对本申请所涉及的应用场景进行介绍:
随着移动通信网络的快速发展,围绕着伪基站所展开的违法犯罪活动数量不断增加,覆盖范围越来越广,不法分子通过伪基站对移动用户发送大量非法广告和诈骗相关的短信,不仅侵犯了群众隐私、损害群众合法利益,还导致社会公共秩序乃至国家通讯秩序收到严重干扰。因此,对于伪基站的快速精准的定位已经成为网络质量提升的重要方向之一。
用户的终端在连接态时,会根据事件或者周期性上报测量报告(Measure Report,MR),MR中包含采集时间、终端唯一标识、(Assisted GPS,A-GPS)经纬度,可以实时上报用户位置。当接入伪基站时,终端将进行自动重选,接入相对应的伪基站小区,并且在接入伪基站后无法正常上报MR,导致用户连续上报的测量报告出现中断。在离开伪基站覆盖范围后,终端恢复正常上报MR。所以,可以聚类多用户的移动轨迹中出现的断点定位伪基站位置,准确性和效率较高。
本申请提供的方法可由电子设备如处理器执行相应的软件代码实现,也可由电子设备在执行相应的软件代码的同时,通过和服务器进行数据交互来实现,例如服务器控制电子设备实现该定位方法。电子设备和服务器之间可以通过网络连接。
下面以具体的实施例对本申请的技术方案进行详细说明。下面这几个具体的实施例可以相互结合,对于相同或相似的概念或过程可能在某些实施例不再赘述。
图1是本申请提供的伪基站的位置定位方法一实施例的流程示意图。如图1所示,本实施例提供的方法,包括:
步骤101、获取至少一个终端周期性上报的测量报告;测量报告包括:终端的标识和位置;
步骤102、根据各个终端上报的测量报告,生成各个终端的至少两条轨迹记录;
具体的,可以从网络管理平台获取终端周期性上报的测量报告数据,测量报告数据例如包括终端唯一标识、位置,还可以包括数据采集时间,其中位置例如可以通过经纬度表示。
其中,经纬度可以通过A-GPS定位方式得到,或还可以通过其他定位方式得到。
在一实施例中,步骤102可以通过如下方式实现:
针对任一所述终端,将所述终端在连续的至少两个周期内上报的测量报告中的位置,按时间顺序形成所述终端的至少两条轨迹记录;每一条所述轨迹记录中包括至少两个轨迹点,任意相邻的两个所述轨迹点形成第二子轨迹;每个所述轨迹点的位置为对应上报周期的测量报告中的终端的位置。
具体的,将终端唯一标识相同且上报周期连续的位置按时间顺序归纳为一条轨迹记录,轨迹记录可以通过轨迹记录表格的形式表示。相同终端唯一标识下可生成一条或多条轨迹记录。
上报周期连续例如指若有第1周期,第2周期,第3周期的测量报告,则该三个测量报告为上报周期连续的测量报告,若有第1周期,第3周期的测量报告,则该2个测量报告为上报周期不连续的测量报告。
网络管理平台例如包括移动通信***中的无线网络管理平台、核心网管理平台及其他类型的能够记录用户上报的测量报告的综合网络管理平台。终端唯一标识例如包括:国际移动设备识别码(International Mobile Equipment Identity,IMEI)、国际移动用户识别码(International Mobile Subscriber Identification Number,IMSI)、手机号码(Telephone number,Telnum)。在本申请实施例中,以IMEI为例,将上报周期连续的测量报告数据形成一条轨迹记录,记录在表格中。同一终端唯一标识下可以包含多条轨迹记录,示例如下表1,其中,采样点序列表示终端的轨迹点(某一位置)的序号,T为上报周期:
表1
步骤103、针对任一终端,从至少两条轨迹记录中提取满足预设的连续性条件的第一子轨迹,并从任意相邻的两条轨迹记录中的第一子轨迹中提取作为断点的轨迹点。
具体的,根据生成的轨迹记录,将轨迹记录内任意两个上报周期连续的轨迹点之间的连线记为一段第二子轨迹Xn,通过每个轨迹点的位置计算轨迹记录内第二子轨迹的距离,轨迹点的速度、加速度,以及第二子轨迹之间的转角,可以记录在参数表中;
在一实施例中,对于任一所述第二子轨迹Xn,根据所述第二子轨迹Xn中轨迹点n以及轨迹点n+1的位置,确定所述第二子轨迹Xn的距离;n为大于0的整数;所述第二子轨迹Xn为相邻的轨迹点n以及轨迹点n+1形成的;
在一实施例中,根据所述第二子轨迹Xn的距离,确定所述第二子轨迹Xn中轨迹点n+1的速度;
在一实施例中,根据所述第二子轨迹Xn-1的距离,确定所述第二子轨迹Xn-1中第二轨迹点n的速度;所述第二子轨迹Xn-1为相邻的轨迹点n-1以及轨迹点n形成的;
在一实施例中,根据所述轨迹点n+1的速度以及所述轨迹点n的速度,确定所述轨迹点n+1的加速度;
在一实施例中,根据所述第二子轨迹Xn的距离,所述第二子轨迹Xn-1的距离,以及轨迹点n-1和轨迹点n+1的距离,确定所述第二子轨迹Xn和所述第二子轨迹Xn-1之间的转角。
具体的,如图2所示,第二子轨迹Xn表示轨迹点n以及轨迹点n+1之间的子轨迹,距离记为Dn,可以使用球面距离公式计算:
Dn=
R·arc cos[cos(latn)×cos(latn+1)×2×cos(lonn-lonn+1)+sin(latn)sin(latn+1)]
其中,lonn、latn分别表示轨迹点n的经度和纬度,lonn+1、latn+1分别表示轨迹点n+1的经度和纬度。R表示地球的平均半径。
轨迹点n+1速度记为Vn+1,T为周期时间,可以采用如下公式计算:
Vn+1=Dn/T
Vn+1=Dn-1/T
轨迹点n+1的加速度记为αn+1,T为周期时间,计算公式为:
第二子轨迹Xn和第二子轨迹Xn-1之间的转角,记为θn,S为子轨迹向量角,子轨迹转角θn为子轨迹向量角的补角,计算公式为:
其中,d(n-1,n+1)表示轨迹点n-1与轨迹点n+1之间的距离。
例如表2中记录的信息:
表2
在一实施例中,满足预设的连续性条件的第一子轨迹可以通过如下方式形成:
针对任一条所述轨迹记录,将所述轨迹记录中,速度相等的轨迹点形成的任意相邻的第二子轨迹的转角小于第一预设阈值的轨迹点,形成所述第一子轨迹;或,
将所述轨迹记录中,加速度相等的轨迹点形成的任意相邻的第二子轨迹的转角小于第一预设阈值的轨迹点,形成所述第一子轨迹;其中,所述第一子轨迹中任意相邻的轨迹点的上报周期连续。
具体的,对轨迹记录内所有的轨迹点进行进一步整理,按照时间顺序遍历轨迹记录内的所有轨迹点,对于上报周期连续且轨迹点速度相等且第二子轨迹之间的转角小于预设阈值θthresh的采样点,或上报周期连续且轨迹点加速度相等且第二子轨迹之间的转角小于预设阈值θthresh的轨迹点归纳形成一段连续的第一子轨迹,当出现轨迹点的速度以及加速度均不相等或第二子轨迹之间的转角大于预设阈值θthresh的轨迹点,则判定此轨迹点属于另一段连续子轨迹,并以此轨迹点为新的第一子轨迹的第一个轨迹点,继续进行遍历。
如图4所示,对于轨迹记录L1,该轨迹记录中每相邻两个轨迹点,记为一个第二子轨迹,例如X1,X2。对轨迹记录进行第一子轨迹划分,基于速度相等,且第二子轨迹的转角小于第一预设阈值,划分成第一子轨迹P1、P2、P3、P4;基于加速度相等,且第二子轨迹的转角小于第一预设阈值,划分成第一子轨迹S1、S2、S3、S4、S5。
对轨迹记录内的轨迹点进行分段后添加新字段:第一子轨迹序号,代表该轨迹点属哪一段第一子轨迹。如下表3所示:
表3
在提取出满足预设的连续性条件的第一子轨迹之后,可以从相邻的轨迹记录中的第一子轨迹中提取作为断点的轨迹点,例如保留相邻的每条轨迹记录中的第一条第一子轨迹和最后一条第一子轨迹,假设前一条轨迹记录中的最后一条第一子轨迹为a,后一条轨迹记录中的第一条第一子轨迹为b,若第一子轨迹为a和第一子轨迹为b之间满足一定条件,则取第一子轨迹为a的最后一个轨迹点(时间上最后一个)以及第一子轨迹为b的第一个轨迹点(时间上第一个)作为断点。其中,满足一定条件例如是指两个子轨迹之间的时间间隔和/或距离间隔在一定的范围内。
步骤104、对作为断点的轨迹点进行聚类处理,得到至少一个聚类簇;聚类簇包括至少两个轨迹点。
在实施例中,可以通过聚类算法,例如基于密度的聚类算法,如DBSCAN算法,对多个终端的多个作为断点的轨迹点进行聚类处理,得到至少一个聚类簇。一个聚类簇可以表示一个伪基站的范围。
步骤105、将各个聚类簇中的轨迹点的位置进行加权处理,得到各个聚类簇中伪基站的位置。
具体的,可以通过聚类簇内所有作为断点的轨迹点的经纬度进行加权平均,即对于每一个伪基站的范围均能确定唯一一个伪基站Kn位置,其中n为聚类簇的序号,M为聚类簇内轨迹点的数量。例如通过如下公式计算得到伪基站的位置:
其中,lon经表示伪基站的位置的经度,lat伪表示伪基站的位置的纬度。
在本步骤中,将所有的聚类簇进行伪基站的位置确认,确认伪基站位置后可对其展开相关手段处理,提升网络提升和规划建设工作效率。
本实施例中,获取至少一个终端周期性上报的测量报告;所述测量报告包括:所述终端的标识和位置;根据各个所述终端上报的测量报告,生成各个所述终端的至少两条轨迹记录;针对任一所述终端,从至少两条所述轨迹记录中提取满足预设的连续性条件的第一子轨迹,并从任意相邻的两条所述轨迹记录中的所述第一子轨迹中提取作为断点的轨迹点;对所述作为断点的轨迹点进行聚类处理,得到至少一个聚类簇;所述聚类簇包括至少两个轨迹点;将各个所述聚类簇中的轨迹点的位置进行加权处理,得到各个所述聚类簇中伪基站的位置,上述通过终端上报的测量报告生成轨迹记录,进而根据轨迹记录中出现的断点确定伪基站的位置,准确性和效率较高。
在一实施例中,在生成轨迹记录之后,可以对轨迹记录内的异常轨迹点的位置进行修正,若轨迹记录中存在轨迹点的加速度大于预设加速度阈值,则根据轨迹点相邻的轨迹点的位置,对轨迹点的位置进行修正。
具体的,轨迹点的加速度超过预设加速度阈值αthresh的轨迹点可以认为是异常轨迹点,通过该异常轨迹点的前后相邻的轨迹点的经纬度对该异常轨迹点的经纬度进行修正。如图3所示,轨迹点n为异常点。
若轨迹点的加速度αn大于预设加速度阈值αthresh则认为该轨迹点由于定位不准确而导致轨迹偏离度过高,并通过异常轨迹点前后相邻的轨迹点的经纬度进行加权平均对该异常轨迹点的位置进行修正。本步骤中可以通过如下公式修正异常轨迹点的位置,其中n代表异常点序号:
lonn=(lonn-1+lonn+1)/2
latmm=(latn-1+latn+1)/2
在一实施例中,提取作为断点的轨迹点,可以通过如下方式实现:
确定任意相邻的两条所述轨迹记录中最后一条所述第一子轨迹,与第一条所述第一子轨迹之间的时间间隔和距离间隔;
若所述时间间隔处于预设时间范围内,且所述距离间隔处于预设距离范围内,则将最后一条所述第一子轨迹中最后一个轨迹点,以及第一条所述第一子轨迹中第一个轨迹点,确定为所述作为断点的轨迹点。
在一实施例中,最后一条所述第一子轨迹和第一条所述第一子轨迹满足预设的移动性条件。移动性条件在下述实施例中进行详细介绍。满足移动性条件即有移动性,否则无移动性。
具体的,计算相同终端唯一标识下不同轨迹记录的时间间隔与距离间隔。若第一子轨迹的时间间隔与距离间隔满足阈值区间则将该两条第一子轨迹记为一个断点轨迹对,即一个断点轨迹对下包含两段第一子轨迹;
在该实施例中:第一步,基于轨迹记录的首尾第一子轨迹,计算相同终端唯一标识下,不同轨迹记录之间的时间间隔与距离间隔。
例如,计算相同终端唯一标识下不同轨迹记录的时间间隔与距离间隔:计算3515xxx2903号码下第一条轨迹记录Ln与第二条轨迹记录Ln+1的时间间隔与距离间隔,第一条轨迹记录下包含两条第一子轨迹Sn,1,Sn,k;第二条轨迹记录下包含两条第一子轨迹S(n+1),1,S(n+1),k,以第一条轨迹记录中的最后一段第一子轨迹Sn,k的最后一个轨迹点与第二条轨迹记录中的第一段第一子轨迹S(n+1),1的第一个轨迹点根据采集时间进行时间间隔与距离的计算,计算方法参见前述实施例,在此不做赘述。采样点的选择如下表4中第4、第5行所示。
表4
当出现某条轨迹记录Ln的上一条用户轨迹Ln-1或下一条用户轨迹记录Ln+1内的第一子轨迹无移动性时,则可以不进行计算。如轨迹记录Ln内经过筛选与移动性判断后包含两条第一子轨迹Sn,1,Sn,k,用户轨迹记录Ln-1内经过筛选与移动性判断后只包含一条第一子轨迹Sn-1,k,用户轨迹记录Ln+1内经过筛选与移动性判断后只包含一条连续子轨迹Sn+1,w,其中,轨迹记录Ln+1内有w条第一子轨迹,Sn+1,w为轨迹记录Ln+1内的最后一条第一子轨迹。如图5所示,计算相邻的轨迹记录Ln-1和Ln中,第一子轨迹Sn-1,k与Sn,1之间的距离间隔和时间间隔,从而提取作为断点的轨迹点。
完成轨迹记录间的时间间隔与距离间隔计算后,形成IMET号码、轨迹记录序号、第一子轨迹序号、轨迹记录间距离的表格,如下表5所示:
表5
第二步:对轨迹记录之间时间间隔与距离间隔进行判断,若轨迹记录之间的时间间隔与距离间隔分别满足预设的时间范围和距离范围,则将对应的两个第一子轨迹标记为一个断点轨迹对。即
其中,sn表示第一子轨迹Sn,k,s(n+1)表示第一子轨迹Sn+1,1。
通过上述步骤形成IMEI号码、轨迹记录序号、第一子轨迹序号与断点轨迹对序号一一对应的断点轨迹对表格,如下表6所示:
表6
进一步,根据断点轨迹对,筛选出作为断点的轨迹点。将断点轨迹对中前一个第一子轨迹的最后一个轨迹点,以及后一个第一子轨迹的第一个轨迹点,并将保留的轨迹点形成断点轨迹点集合;其中,前一个第一子轨迹和后一个第一子轨迹例如分别为S1,k和S2,1。
可以结合轨迹记录以及断点轨迹对表格两张表格,将断点轨迹对表格内的各断点轨迹对内作为断点的轨迹点提取出来。分别提取第一段第一子轨迹的最后一个轨迹点以及第二段第一子轨迹的第一个轨迹点。如下表7中的1号断点轨迹对,则提取表8中轨迹记录中IMEI号码=3515xxx2903,轨迹记录L1下,第一子轨迹S1,k下的最后一个轨迹点;提取轨迹记录中IMEI号码=3515xxx2903,用户轨迹记录L2下,连续子轨迹S2,1下的第一个轨迹点。
表7
表8
以上,结合轨迹记录以及断点轨迹对表格将所有作为断点的轨迹点提取后,最终形成IMEI号码、轨迹序号、采集时间、经度、纬度以及第一子轨迹序号的断点轨迹点表格,如下表9所示:
表9
在一实施例中,若所述第一子轨迹中各个轨迹点的速度的加权平均值大于第一阈值,所述第一子轨迹包括的上报周期之和大于第二阈值,则确定所述第一子轨迹满足预设的移动性条件。
具体的,在确定断点轨迹对时,可以保留轨迹记录中第一条第一子轨迹和最后一条第一子轨迹,并对第一子轨迹的移动性进行判断,移动性判断可以基于以下两个:1、第一子轨迹内所有轨迹点速度的加权平均值是否超过第一阈值;2、第一子轨迹内所有轨迹点上报周期之和是否超过第二阈值,满足上述两个条件的第一子轨迹可以认为是具有移动性的第一子轨迹,否则舍弃该段第一子轨迹,不进行断点轨迹对的判断。
例如,某条轨迹记录Lm中包含多段第一子轨迹S1,S2,S3,S4...Sk,则保留第一子轨迹S1,Sk。随后对于第一子轨迹移动性判断:该条轨迹记录中的第一条第一子轨迹S1内所有轨迹点的速度加权平均值表示为VS1,所有轨迹点的上报周期之和为TS1,当VS1>第一阈值Vthresh,且TS1>第二阈值Tthresh时,认为该第一子轨迹S1具有移动性。
在本步骤中,对轨迹记录内的第一子轨迹进行筛选与移动性判断后添加新字段:连续子轨迹移动性,如下表10所示:
表10
在一实施例中,聚类处理流程如下:
根据S7所确定的断点采样点表格,对作为断点的轨迹点的集合使用基于DBSCAN的聚类方法,得到核心点、边界点和噪声点,并得到不同聚类簇,一个聚类簇对应一个伪基站的范围;
具体实施流程如下:
1、设定半径E=Ddbscan,MinPts=Ndbscan,初始化聚类簇数K=0;MinPts表示某一样本的距离为E的邻域中样本个数的阈值。
2、给所有作为断点的轨迹点赋一个unvisited的属性,表示这个轨迹点未被访问;
3、在所有属性为unvisited的轨迹点的集合中,随机找一个点P标记为visited,并生成邻域Eps(P);
4、核查点P类型:若在邻域Eps(P)内含有超过MinPts数目的点,则将点P设定为核心点。若P为核心点则执行步骤5;若不是核心点标记为噪声点,返回至步骤2重新选点;
5、建立一个以P为核心点的簇C,将P点邻域Eps(P)内的所有点加入候选集合N;
6、选择候选集合N中的点P’,判断P’是否已被访问,若已被访问且无类别,则将P’加入簇C;若P’未被访问则执行步骤7;
7、计算P’邻域集合Eps(P’),若P’判断为核心点则加入簇C,且将Eps(P’)中所有点的候选集合N’加入候选集合N中;
8、遍历候选集N中所有点直到所有点均被分类则完成一簇的聚类,K=K+1;
9、遍历轨迹点的集合中所有unvisited直到所有轨迹点均被访问完成遍历。
本实施例中,由于终端受到伪基站干扰可能无法上报测量报告从而出现中断,因此可以通过终端上报的测量报告生成轨迹记录,进而根据轨迹记录中出现的断点确定伪基站的位置,准确性和效率较高。
图6为本申请提供的伪基站的位置定位装置一实施例的结构图,如图6所示,本实施例的伪基站的位置定位装置,包括:
获取模块601,用于获取至少一个终端周期性上报的测量报告;所述测量报告包括:所述终端的标识和位置;
预处理模块602,用于根据各个所述终端上报的测量报告,生成各个所述终端的至少两条轨迹记录;
所述预处理模块602,还用于针对任一所述终端,从至少两条所述轨迹记录中提取满足预设的连续性条件的第一子轨迹,并从任意相邻的两条所述轨迹记录中的所述第一子轨迹中提取作为断点的轨迹点;
所述预处理模块602,还用于对所述作为断点的轨迹点进行聚类处理,得到至少一个聚类簇;所述聚类簇包括至少两个轨迹点;
处理模块603,用于将各个所述聚类簇中的轨迹点的位置进行加权处理,得到各个所述聚类簇中伪基站的位置。
在一种可能的实现方式中,所述预处理模块602,具体用于:
针对任一条所述轨迹记录,将所述轨迹记录中,速度相等的轨迹点形成的任意相邻的第二子轨迹的转角小于第一预设阈值的轨迹点,形成所述第一子轨迹;或,
将所述轨迹记录中,加速度相等的轨迹点形成的任意相邻的第二子轨迹的转角小于第一预设阈值的轨迹点,形成所述第一子轨迹;其中,所述第一子轨迹中任意相邻的轨迹点的上报周期连续。
在一种可能的实现方式中,所述预处理模块602,具体用于:
确定任意相邻的两条所述轨迹记录中最后一条所述第一子轨迹,与第一条所述第一子轨迹之间的时间间隔和距离间隔;
若所述时间间隔处于预设时间范围内,且所述距离间隔处于预设距离范围内,则将最后一条所述第一子轨迹中最后一个轨迹点,以及第一条所述第一子轨迹中第一个轨迹点,确定为所述作为断点的轨迹点,其中,最后一条所述第一子轨迹和第一条所述第一子轨迹满足预设的移动性条件。
在一种可能的实现方式中,所述预处理模块602,还用于:
若所述第一子轨迹中各个轨迹点的速度的加权平均值大于第一阈值,所述第一子轨迹包括的上报周期之和大于第二阈值,则确定所述第一子轨迹满足预设的移动性条件。
在一种可能的实现方式中,所述预处理模块602,还用于:
对于任一所述第二子轨迹Xn,根据所述第二子轨迹Xn中轨迹点n以及轨迹点n+1的位置,确定所述第二子轨迹Xn的距离;n为大于0的整数;所述第二子轨迹Xn为相邻的轨迹点n以及轨迹点n+1形成的;
根据所述第二子轨迹Xn的距离,确定所述第二子轨迹Xn中轨迹点n+1的速度;
根据所述第二子轨迹Xn-1的距离,确定所述第二子轨迹Xn-1中第二轨迹点n的速度;所述第二子轨迹Xn-1为相邻的轨迹点n-1以及轨迹点n形成的;
根据所述轨迹点n+1的速度以及所述轨迹点n的速度,确定所述轨迹点n+1的加速度;
根据所述第二子轨迹Xn的距离,所述第二子轨迹Xn-1的距离,以及轨迹点n-1和轨迹点n+1的距离,确定所述第二子轨迹Xn和所述第二子轨迹Xn-1之间的转角。
在一种可能的实现方式中,所述预处理模块602,还用于:
若所述轨迹记录中存在轨迹点的加速度大于预设加速度阈值,则根据所述轨迹点相邻的轨迹点的位置,对所述轨迹点的位置进行修正。
在一种可能的实现方式中,所述预处理模块602,具体用于:
对所述作为断点的轨迹点用基于密度的聚类算法进行聚类处理,得到至少一个聚类簇,所述聚类簇包括多个所述作为断点的轨迹点。
在一种可能的实现方式中,所述处理模块603,具体用于:
针对任一所述聚类簇,将所述聚类簇中作为断点的轨迹点的经度和纬度分别进行加权处理,得到所述聚类簇中伪基站的位置。
本实施例的装置,可以用于执行上述方法实施例的技术方案,其实现原理和技术效果类似,此处不再赘述。
图7为本申请提供的电子设备实施例的结构图,如图7所示,该电子设备包括:
处理器701,以及,用于存储处理器701的可执行指令的存储器702。
可选的,还可以包括:通信接口703,用于实现与其他设备的通信。
上述部件可以通过一条或多条总线进行通信。
其中,处理器701配置为经由执行所述可执行指令来执行前述方法实施例中对应的方法,其具体实施过程可以参见前述方法实施例,此处不再赘述。
本申请实施例中还提供一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现前述方法实施例中对应的方法,其具体实施过程可以参见前述方法实施例,其实现原理和技术效果类似,此处不再赘述。
本领域技术人员在考虑说明书及实践这里公开的发明后,将容易想到本公开的其它实施方案。本申请旨在涵盖本公开的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本公开的一般性原理并包括本公开未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本公开的真正范围和精神由下面的权利要求书指出。
应当理解的是,本公开并不局限于上面已经描述并在附图中示出的精确结构,并且可以在不脱离其范围进行各种修改和改变。本公开的范围仅由所附的权利要求书来限制。
Claims (11)
1.一种伪基站的位置定位方法,其特征在于,包括:
获取至少一个终端周期性上报的测量报告;所述测量报告包括:所述终端的标识和位置;
根据各个所述终端上报的测量报告,生成各个所述终端的至少两条轨迹记录;
针对任一所述终端,从至少两条所述轨迹记录中提取满足预设的连续性条件的第一子轨迹,并从任意相邻的两条所述轨迹记录中的所述第一子轨迹中提取作为断点的轨迹点;
对所述作为断点的轨迹点进行聚类处理,得到至少一个聚类簇;所述聚类簇包括至少两个轨迹点;
将各个所述聚类簇中的轨迹点的位置进行加权处理,得到各个所述聚类簇中伪基站的位置;
其中,
所述从任意相邻的两条所述轨迹记录中的所述第一子轨迹中,提取作为断点的轨迹点,包括:
确定任意相邻的两条所述轨迹记录中最后一条所述第一子轨迹,与第一条所述第一子轨迹之间的时间间隔和距离间隔;
若所述时间间隔处于预设时间范围内,且所述距离间隔处于预设距离范围内,则将最后一条所述第一子轨迹中最后一个轨迹点,以及第一条所述第一子轨迹中第一个轨迹点,确定为所述作为断点的轨迹点,其中,最后一条所述第一子轨迹和第一条所述第一子轨迹满足预设的移动性条件。
2.根据权利要求1所述的方法,其特征在于,所述根据各个所述终端上报的测量报告,生成各个所述终端的至少两条轨迹记录,包括:
针对任一所述终端,将所述终端在连续的至少两个周期内上报的测量报告中的位置信息,按时间顺序形成所述终端的至少两条轨迹记录;每一条所述轨迹记录中包括至少两个轨迹点,任意相邻的两个所述轨迹点形成第二子轨迹;每个所述轨迹点的位置为对应上报周期的测量报告中的终端的位置。
3.根据权利要求2所述的方法,其特征在于,所述从至少两条所述轨迹记录中提取满足预设的连续性条件的第一子轨迹,包括:
针对任一条所述轨迹记录,将所述轨迹记录中,速度相等的轨迹点形成的任意相邻的第二子轨迹的转角小于第一预设阈值的轨迹点,形成所述第一子轨迹;或,
将所述轨迹记录中,加速度相等的轨迹点形成的任意相邻的第二子轨迹的转角小于第一预设阈值的轨迹点,形成所述第一子轨迹;其中,所述第一子轨迹中任意相邻的轨迹点的上报周期连续。
4.根据权利要求1所述的方法,其特征在于,所述将最后一条所述第一子轨迹中最后一个轨迹点,以及第一条所述第一子轨迹中第一个轨迹点,确定为所述作为断点的轨迹点之前,还包括:
若所述第一子轨迹中各个轨迹点的速度的加权平均值大于第一阈值,所述第一子轨迹包括的上报周期之和大于第二阈值,则确定所述第一子轨迹满足预设的移动性条件。
5.根据权利要求3所述的方法,其特征在于,所述从至少两条所述轨迹记录中提取满足预设的连续性条件的第一子轨迹之前,还包括:
对于任一所述第二子轨迹Xn,根据所述第二子轨迹Xn中轨迹点n以及轨迹点n+1的位置,确定所述第二子轨迹Xn的距离;n为大于0的整数;所述第二子轨迹Xn为相邻的轨迹点n以及轨迹点n+1形成的;
根据所述第二子轨迹Xn的距离,确定所述第二子轨迹Xn中轨迹点n+1的速度;
根据所述第二子轨迹Xn-1的距离,确定所述第二子轨迹Xn-1中第二轨迹点n的速度;所述第二子轨迹Xn-1为相邻的轨迹点n-1以及轨迹点n形成的;
根据所述轨迹点n+1的速度以及所述轨迹点n的速度,确定所述轨迹点n+1的加速度;
根据所述第二子轨迹Xn的距离,所述第二子轨迹Xn-1的距离,以及轨迹点n-1和轨迹点n+1的距离,确定所述第二子轨迹Xn和所述第二子轨迹Xn-1之间的转角。
6.根据权利要求1-5任一项所述的方法,其特征在于,所述从至少两条所述轨迹记录中提取满足预设的连续性条件的第一子轨迹之前,还包括:
若所述轨迹记录中存在轨迹点的加速度大于预设加速度阈值,则根据所述轨迹点相邻的轨迹点的位置,对所述轨迹点的位置进行修正。
7.根据权利要求1-5任一项所述的方法,其特征在于,所述对所述作为断点的轨迹点进行聚类处理,得到至少一个聚类簇,包括:
对所述作为断点的轨迹点用基于密度的聚类算法进行聚类处理,得到至少一个聚类簇,所述聚类簇包括多个所述作为断点的轨迹点。
8.根据权利要求1-5任一项所述的方法,其特征在于,所述将各个所述聚类簇中的轨迹点的位置进行加权处理,得到各个所述聚类簇中伪基站的位置,包括:
针对任一所述聚类簇,将所述聚类簇中作为断点的轨迹点的经度和纬度分别进行加权处理,得到所述聚类簇中伪基站的位置。
9.一种伪基站的位置定位装置,其特征在于,包括:
获取模块,用于获取至少一个终端周期性上报的测量报告;所述测量报告包括:所述终端的标识和位置;
预处理模块,用于根据各个所述终端上报的测量报告,生成各个所述终端的至少两条轨迹记录;
所述预处理模块,还用于针对任一所述终端,从至少两条所述轨迹记录中提取满足预设的连续性条件的第一子轨迹,并从任意相邻的两条所述轨迹记录中的所述第一子轨迹中提取作为断点的轨迹点;
所述预处理模块,还用于对所述作为断点的轨迹点进行聚类处理,得到至少一个聚类簇;所述聚类簇包括至少两个轨迹点;
处理模块,用于将各个所述聚类簇中的轨迹点的位置进行加权处理,得到各个所述聚类簇中伪基站的位置;
其中,
所述预处理模块,还用于:
确定任意相邻的两条所述轨迹记录中最后一条所述第一子轨迹,与第一条所述第一子轨迹之间的时间间隔和距离间隔;
若所述时间间隔处于预设时间范围内,且所述距离间隔处于预设距离范围内,则将最后一条所述第一子轨迹中最后一个轨迹点,以及第一条所述第一子轨迹中第一个轨迹点,确定为所述作为断点的轨迹点,其中,最后一条所述第一子轨迹和第一条所述第一子轨迹满足预设的移动性条件。
10.一种电子设备,其特征在于,包括:
处理器;以及
存储器,用于存储所述处理器的可执行指令;
其中,所述处理器配置为经由执行所述可执行指令来执行如权利要求1-8任一项所述的方法。
11.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1-8任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911197361.9A CN111093156B (zh) | 2019-11-29 | 2019-11-29 | 伪基站的位置定位方法、设备和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911197361.9A CN111093156B (zh) | 2019-11-29 | 2019-11-29 | 伪基站的位置定位方法、设备和存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111093156A CN111093156A (zh) | 2020-05-01 |
CN111093156B true CN111093156B (zh) | 2020-12-15 |
Family
ID=70393865
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911197361.9A Active CN111093156B (zh) | 2019-11-29 | 2019-11-29 | 伪基站的位置定位方法、设备和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111093156B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111966729B (zh) * | 2020-08-27 | 2023-05-19 | 平安国际智慧城市科技股份有限公司 | 一种车辆轨迹数据处理方法、装置、设备及存储介质 |
CN112866192B (zh) * | 2020-12-30 | 2022-11-04 | 绿盟科技集团股份有限公司 | 一种识别异常聚集行为的方法及装置 |
CN115373002A (zh) * | 2021-05-19 | 2022-11-22 | 中兴通讯股份有限公司 | 一种道路用户识别方法、装置、存储介质及电子装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102273174A (zh) * | 2009-01-05 | 2011-12-07 | 高通股份有限公司 | 对伪造无线接入点的检测 |
CN105792238A (zh) * | 2014-12-24 | 2016-07-20 | ***通信集团上海有限公司 | 一种伪基站定位方法和装置 |
CN105992210A (zh) * | 2015-02-09 | 2016-10-05 | ***通信集团湖北有限公司 | 移动伪基站定位方法、***、定位设备及定位服务器 |
CN106231597A (zh) * | 2015-06-02 | 2016-12-14 | 中国科学院上海高等研究院 | 一种伪基站的定位方法和*** |
CN107371164A (zh) * | 2017-08-02 | 2017-11-21 | 西北大学 | 一种基于传感器数据与信号差值融合定位伪ap的方法 |
CN109195165A (zh) * | 2018-08-03 | 2019-01-11 | 广州瀚信通信科技股份有限公司 | 一种基于mr数据的移动通信4g伪基站识别方法 |
CN109302706A (zh) * | 2017-07-24 | 2019-02-01 | ***通信有限公司研究院 | 一种伪基站的定位方法、装置及计算机可读存储介质 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6163696A (en) * | 1996-12-31 | 2000-12-19 | Lucent Technologies Inc. | Mobile location estimation in a wireless communication system |
US7286515B2 (en) * | 2003-07-28 | 2007-10-23 | Cisco Technology, Inc. | Method, apparatus, and software product for detecting rogue access points in a wireless network |
DE602006013890D1 (de) * | 2006-06-30 | 2010-06-02 | Ericsson Telefon Ab L M | Verfahren, verwendung des verfahrens und anordnungen in einem elektronischen support-massnahmensystem |
WO2013001884A1 (ja) * | 2011-06-29 | 2013-01-03 | インターナショナル・ビジネス・マシーンズ・コーポレーション | 情報通信端末、エージェント・プログラム、位置推定システムおよび位置推定方法 |
CN102736092B (zh) * | 2012-07-05 | 2013-10-30 | 北京建筑工程学院 | 基于单测距基站与gps的电子标签协同定位方法 |
-
2019
- 2019-11-29 CN CN201911197361.9A patent/CN111093156B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102273174A (zh) * | 2009-01-05 | 2011-12-07 | 高通股份有限公司 | 对伪造无线接入点的检测 |
CN105792238A (zh) * | 2014-12-24 | 2016-07-20 | ***通信集团上海有限公司 | 一种伪基站定位方法和装置 |
CN105992210A (zh) * | 2015-02-09 | 2016-10-05 | ***通信集团湖北有限公司 | 移动伪基站定位方法、***、定位设备及定位服务器 |
CN106231597A (zh) * | 2015-06-02 | 2016-12-14 | 中国科学院上海高等研究院 | 一种伪基站的定位方法和*** |
CN109302706A (zh) * | 2017-07-24 | 2019-02-01 | ***通信有限公司研究院 | 一种伪基站的定位方法、装置及计算机可读存储介质 |
CN107371164A (zh) * | 2017-08-02 | 2017-11-21 | 西北大学 | 一种基于传感器数据与信号差值融合定位伪ap的方法 |
CN109195165A (zh) * | 2018-08-03 | 2019-01-11 | 广州瀚信通信科技股份有限公司 | 一种基于mr数据的移动通信4g伪基站识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111093156A (zh) | 2020-05-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111093156B (zh) | 伪基站的位置定位方法、设备和存储介质 | |
US8618984B2 (en) | Selecting beacons for location inference | |
CN109165527B (zh) | 支持个性化隐私的轨迹敏感数据保护方法 | |
CN108271157B (zh) | 一种伪基站识别方法及装置 | |
CN108243421B (zh) | 伪基站识别方法及*** | |
CN104596507A (zh) | 一种移动终端出行轨迹的确定方法 | |
CN108574934B (zh) | 一种伪基站定位方法和装置 | |
CN111107556B (zh) | 一种移动通信网络的信号覆盖质量评估方法和装置 | |
US20170180940A1 (en) | System and method for geography-based correlation of cellular and wlan identifiers | |
CN110769459B (zh) | 一种移动终端的定位方法、装置、存储介质及服务器 | |
US10425831B2 (en) | Improving coverage in a cellular network using temporary subscriber indentifiers | |
CN108646261B (zh) | 一种车辆gps数据集的数据处理方法和装置 | |
CN111818537A (zh) | 小区覆盖情况确定方法、装置、设备及可存储介质 | |
CN102480784A (zh) | 一种指纹定位误差的评估方法和*** | |
WO2016127880A1 (zh) | 一种确定离线定位数据质量的方法和装置 | |
CN104540138A (zh) | 一种定位伪基站的方法 | |
CN105246157B (zh) | 一种移动终端的定位方法和定位*** | |
CN107155186A (zh) | 一种伪基站定位方法和装置 | |
CN113645625B (zh) | 伪基站定位方法、装置、电子设备和可读介质 | |
CN108289279A (zh) | 定位信息的处理方法、装置和计算机可读存储介质 | |
CN111194001A (zh) | Lte指纹定位校正的方法、装置及*** | |
Chen et al. | A travel mode identification framework based on cellular signaling data | |
CN114697969B (zh) | 一种干扰源定位方法和装置 | |
WO2015084154A1 (en) | A system and method for locating a mobile device | |
CN108260204A (zh) | 一种基于电信定位的方法及装置 |
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 |