一种无线基站与线段型地理地物关联矩阵的建立方法

摘要

本发明属于数据分析技术领域，具体涉及一种无线基站与线段型地理地物关联矩阵的建立方法。包括以下步骤提取训练数据；提取地物信息，建立；对线段型地物进行分段；通过地理关系关联TrainData和PoiData；建立基于训练数据的关系矩阵Rcp_train、Rpc_train；利用传播模型计算每一个无线基站小区到每一个分段点的覆盖强度，推导出覆盖概率；建立所有小区到所有分段点的覆盖概率矩阵Rline_cov；计算训练矩阵的伴随可信矩阵Bline；S9综合Rpc_train、Rline_cov、Bline三个矩阵对线段型地物建立小区-地物的关系矩阵Rline；S10对关系矩阵Rline采用周期性滑动窗口法进行更新。本发明为准确判断手机用户的停留地点之后将为手机用户提供社会化地理信息提供基础数据，为公共安全、防灾减灾、智慧城市、市场营销等基于人员位置的应用提供基础模型。

权利要求

1.一种无线基站与线段型地理地物关联矩阵的建立方法，其特征在于：具体包括以下步骤：S1、从移动运营商大数据中提取训练数据TrainData；

S2、从电子地图中提取地物信息，建立PoiData；

S3、利用线性差分法对线段型地物进行分段；

S4、通过地理关系关联TrainData和PoiData；

S5、建立基于训练数据的关系矩阵Rcp_train、Rpc_train；

S6、利用传播模型计算每一个无线基站小区到每一个分段点的覆盖强度，推导出覆盖概率；

S7、建立所以小区到所以分段点的覆盖概率矩阵Rline_cov；

S8、计算训练矩阵的伴随可信矩阵Bline；

S9、综合Rpc_train、Rline_cov、Bline三个矩阵对线段型地物建立小区-地物的关系矩阵Rline；

S10、对关系矩阵Rline采用周期性滑动窗口法进行更新。

2.根据权利要求1所述的一种无线基站与线段型地理地物关联矩阵的建立方法，其特征在于：所述步骤S2中的地物信息建立PoiData具体为：PoiData＝{(Poiname,Poitype,Centerx,Centery,Points)}其中：Poiname表示地物名称；Poitype表示地物类型；Centerx和Centery分别表示线段中心点的坐标；Points描述多线段的有序点的集合。

3.根据权利要求2所述的一种无线基站与线段型地理地物关联矩阵的建立方法，其特征在于：所述步骤S3具体为：采用线性差分法对于城市环境按照20～50米进行划分，对于郊区环境按照100～5000米进行划分。

4.根据权利要求3所述的一种无线基站与线段型地理地物关联矩阵的建立方法，其特征在于：所述步骤S4中TrainData和PoiData的的关联方法具体为：对于线段型PoiData，对线段沿法线方向扩展为一个矩形，如果TrainData中的经纬度落在该PoiData对应的Points中，则这个TrainData记录可以关联这个线段型PoiData。

5.根据权利要求4所述的一种无线基站与线段型地理地物关联矩阵的建立方法，其特征在于：所述步骤S5中关系矩阵Rcp_train、Rpc_train的具体建立方法为：对所有TrainData进行计算，然后对地理地物进行合并则得到以下的数据其中，Cnt为满足关联条件的训练数据的数量，

把Cnttrain按行进行归一化得到小区-地物的关系矩阵Rcp_train：把Rcp_train进行转置并按行进行归一化得到地物-小区的关系矩阵Rpc_train：

6.根据权利要求5所述的一种无线基站与线段型地理地物关联矩阵的建立方法，其特征在于：所述步骤S6的计算方法具体为：采用以下传播模型计算小区到线段中心点的无线接收功率：Rx＝Tx-(K1+K2log(d)+K3log(hcell)+K4Diff+K5log(d)*log(hcell)+K6(hms)+K7f(clutter))其中，K1为与pi年率相关的因子，取值为23.5；K2为与距离相关的因子，取值为44.9，K3为天线高度相关因子，取值为5.83；K4为与衍射相关的因子，取值为0；K5为与发射天线和距离相关的因子，取值为-6.55；K6为与接收机高度相关因子，取值为0，K7为地物相关因子，取值为0；

Tx为小区的发射功率，从基站的工程参数中获得；

hcell为小区天线高度，从基站的工程参数中获得；

d为小区到线段中心点的距离；

覆盖概率矩阵Rline_cov的建立方法为：

7.根据权利要求6所述的一种无线基站与线段型地理地物关联矩阵的建立方法，其特征在于：所述步骤S8中伴随可信矩阵Bline的计算方法为：假设经过这个地物的用户数usercnt超过N1，并且平均每用户产生的样本数avg_samplecnt超过N2，则认为该样本空间具有统计意义；则训练矩阵的可信度最高，获得最大的加权10x；x取1，N1取100，N2取1000；

当样本空间低于上述空间时，可信度的计算方法如下：

对训练矩阵中的每一个元素进行可信度计算，即得到可信度矩阵Bline：Bline＝[bi,j]n*m。

8.根据权利要求7所述的一种无线基站与线段型地理地物关联矩阵的建立方法，其特征在于：所述步骤S9具体为：根据Rpc_train、Rline_cov、Bline三个矩阵联合计算，最终的合并公式为：R′line＝Rpc_train.*Bline+Rline_cov对R′line按行进行归一化得到最终的Rline

9.根据权利要求8所述的一种无线基站与线段型地理地物关联矩阵的建立方法，其特征在于：所述步骤S10中周期性滑动窗口法具体为：以T为周期，以nT为数据滑动计算窗口的更新方式对关系矩阵Rregion进行更新，其中，T＝1个月，n＝3，表示关系矩阵每个月将重新计算一次，每次计算所采用的数据为过去3个月的数据。

说明书

一种无线基站与线段型地理地物关联矩阵的建立方法

技术领域

[0001] 本发明属于数据分析技术领域，具体涉及一种无线基站与线段型地理地 物关联矩阵的建立方法。

背景技术

[0002] 在应用运营商大数据的过程中，经常需要把基站信息映射到具体的地理 地物，为进一步分析手机用户的社会化位置属性提供数据基础。

[0003] 随着近距离无线通信和移动网络技术的不断发展，基于位置的服务 (location-based service，LBS)越来越受到人们的重视。全球定位系统  (globalpositioning 
system，GPS)在室外已经可以提供位置信息，但是在 室内环境下由于建筑物的遮挡、楼层
密度较大等特点而无法正常工作。在这 种情况下，出现了利用红外线、WiFi、Zigbee等载波
技术的室内定位方案， 但是都没有达到非常好的效果。WiFi技术由于其能耗较大，且容易
受到信号 的干扰，信号覆盖的范围也仅局限于90m以内的空间，所以该技术的优势并 不明
显；Zigbee技术尽管具有低功耗、低成本、通信效率高等优点，但是其 定位结果不稳定，系
统可靠性不强。

[0004] 申请号为201811550838.2的《一种区块特征指纹识别的方法》提出了 一种区块特征指纹识别的方法。包括以下步骤：S1：获取通信运营商提供的 基站工程参数和地图服务
商提供的地理实体实际位置坐标点集合；S2：根据 基站工程参数计算基站的覆盖面；S3：根
据基站的覆盖面和空间区块实际位 置坐标点集合计算地理实体与基站的匹配关系，形成
地理实体特征指纹。本 发明通过基站工程参数结合地理实体实际位置坐标点集合，将每一
个地理区 块实体与基站的位置关系进行匹配，形成地理实体特征指纹。然而该方法存 在
以下缺点：1、无线射频的覆盖是复杂的，有折射、散射和衍射等影响因 素，采用几何的方法
计算覆盖面积重叠，不能准确地计算无线基站和地理地 物实体的关系；2、该方法不能对道
路等线段型地理地物建立无线基站和地 理地物实体的关系。

发明内容

[0005] 本发明要解决的技术问题在于克服现有技术的不足，提供一种无线基站 与线段型地理地物关联矩阵的建立方法，采用的是基于手机用户上报的精准 的GPS和基站编号，
将GPS精确关联到地理实体，建立GPS-小区-地理地物 的关系，通过海量数据可以建立比较
准确的无线基站＝地理地物的关系概率 矩阵。

[0006] 本发明解决上述技术问题的技术方案如下：一种无线基站与线段型地理 地物关联矩阵的建立方法，具体包括以下步骤：

[0007] S1、从移动运营商大数据中提取训练数据TrainData；

[0008] S2、从电子地图中提取地物信息，建立PoiData；

[0009] S3、利用线性差分法对线段型地物进行分段；

[0010] S4、通过地理关系关联TrainData和PoiData；

[0011] S5、建立基于训练数据的关系矩阵Rcp_train、Rpc_train；

[0012] S6、利用传播模型计算每一个无线基站小区到每一个分段点的覆盖强度， 推导出覆盖概率；

[0013] S7、建立所以小区到所以分段点的覆盖概率矩阵Rline_cov；

[0014] S8、计算训练矩阵的伴随可信矩阵Bline；

[0015] S9、综合Rpc_train、Rline_cov、Bline三个矩阵对线段型地物建立小区-地物的 关系矩阵Rline；

[0016] S10、对关系矩阵Rline采用周期性滑动窗口法进行更新。

[0017] 进一步地，所述步骤S2中的地物信息为：高速公路、城市道路、省道、 市道、县道等、出入口、停车场等，并将地物信息建立PoiData具体为：

[0018] PoiData＝{(Poiname,Poitype,Centerx,Centery,Points)}

[0019] 其中：Poiname表示地物名称；Poitype表示地物类型；Centerx和Centery分别 表示线段中心点的坐标；Points描述多线段的有序点的集合。

[0020] 进一步地，所述步骤S3具体为：采用线性差分法对于城市环境按照 20～50米进行划分，对于郊区环境按照100～5000米进行划分。

[0021] 进一步地，所述步骤S4中TrainData和PoiData的的关联方法具体为：对 于线段型PoiData，对线段沿法线方向扩展为一个矩形，如果TrainData中的经 纬度落在该PoiData
对应的Points中，则这个TrainData记录可以关联这个线段 型PoiData。

[0022] 进一步地，所述步骤S5中关系矩阵Rcp_train、Rpc_train的具体建立方法为：

[0023] 对所有TrainData进行计算，然后对地理地物进行合并则得到以下的数据

[0024]

[0025] 其中，Cnt为满足关联条件的训练数据的数量，

[0026] 把Cnttrain按行进行归一化得到小区-地物的关系矩阵Rcp_train：

[0027]

[0028] 把Rcp_train进行转置并按行进行归一化得到地物-小区的关系矩阵Rpc_train：

[0029]

[0030] 进一步地，所述步骤S6的计算方法具体为：采用以下传播模型计算小 区到线段中心点的无线接收功率：

[0031] Rx＝Tx-(K1+K2log(d)+K3log(hcell)+K4Diff+K5log(d)*log(hcell)+K6(hms)

[0032] +K7f(clutter))

[0033] 其中，K1为与pi年率相关的因子，取值为23.5；K2为与距离相关的因 子，取值为44.9，K3为天线高度相关因子，取值为5.83；K4为与衍射相关 的因子，取值为0；K5为与发射
天线和距离相关的因子，取值为-6.55；K6为 与接收机高度相关因子，取值为0，K7为地物相
关因子，取值为0；

[0034] Tx为小区的发射功率，从基站的工程参数中获得；

[0035] hcell为小区天线高度，从基站的工程参数中获得；

[0036] d为小区到线段中心点的距离；

[0037] 覆盖概率矩阵Rline_cov的建立方法为：

[0038]

[0039] 进一步地，所述步骤S8中伴随可信矩阵Bline的计算方法为：假设经过这 个地物的用户数usercnt超过N1，并且平均每用户产生的样本数avg_samplecnt超 过N2，则认为该样本
空间具有统计意义；则训练矩阵的可信度最高，获得最 大的加权10x；x取1，N1取100，N2取
1000；

[0040] 当样本空间低于上述空间时，可信度的计算方法如下：

[0041]

[0042] 对训练矩阵中的每一个元素进行可信度计算，即得到可信度矩阵Bline：

[0043] Bline＝bi,j]n*m。

[0044] 进一步地，所述步骤S9具体为：根据Rcp_train、Rline_cov、Bline三个矩阵联 合计算，最终的合并公式为：

[0045] R'line＝Rpc_train.*Bline+Rline_cov

[0046] 对R'line按行进行归一化得到最终的Rline

[0047]

[0048] 进一步地，所述步骤S10中周期性滑动窗口法具体为：

[0049] 以T为周期，以nT为数据滑动计算窗口的更新方式对关系矩阵Rregion进 行更新，其中，T＝1个月，n＝3，表示关系矩阵每个月将重新计算一次，每次 计算所采用的数据为过去
3个月的数据。

[0050] 本发明采用的是基于手机用户上报的精准的GPS和基站编号，将GPS精 确关联到地理实体，建立GPS-小区-地理地物的关系，通过海量数据可以建 立比较准确的基站＝地
理地物的关系概率矩阵。

[0051] 本发明还将地理地物进行分类，对不同的地物应用不同的计算方法。

[0052] 并发明还提出了一种概率加权合并的方法，把训练得到的关系矩阵与根 据基站和地理实体计算出来的关系矩阵进行合并，从而得到一个更加准确的 关系矩阵。

[0053] 本发明的有益效果是：本发明提出了一种综合利用移动通信网络中的信 令数据、APP数据、基站工程参数以及电子地图中的地物数据联合使用方法。 对一定历史时长的数
据进行提取，采用相对准确的手机上报的GPS数据建立 基站和地理地物的基本关系。在没
有GPS数据或GSP数据不足的地方采用无 线覆盖仿真模型，建立基站和地理地物的关系。并
提出了一个融合方法，把 上述两个关系矩阵融合成统一的关系矩阵。本发明为准确判断手
机用户的停 留地点之后将为手机用户提供社会化地理信息提供基础数据，为公共安全、 
防灾减灾、智慧城市、市场营销等基于人员位置的应用提供基础模型。

附图说明

[0054] 为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面 将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而 易见地，下面描述中
的附图是本发明的一些实施方式，对于本领域普通技术 人员来讲，在不付出创造性劳动的
前提下，还可以根据这些附图获得其他的 附图。

[0055] 图1为本发明具体实施例的步骤S10中采用周期性滑动窗口法对关系矩 阵Rline进行更新的示意图。

具体实施方式

[0056] 下面利用具体实施例对本发明作进一步阐述。在此需要说明的是，对于 这些实施例方式的说明用于帮助理解本发明，但并不构成对本发明的限定。 本文公开的特定结构和
功能细节仅用于描述本发明的示例实施例。然而，可 用很多备选的形式来体现本发明，并
且不应当理解为本发明限制在本文阐述 的实施例中。

[0057] 需要注意的是，除非另有说明，本申请使用的技术术语或者科学术语应 当为本发明所属领域技术人员所理解的通常意义。

[0058] 在本申请的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“长 度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水 平”、“顶”、“底”“内”、“外”、“顺时针”、“逆时针”、“轴向”、“径向”、“周 向”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了 便于描述本发明和简化描述，而不是指示或暗示所指的装置或
元件必须具有 特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。

[0059] 本文使用的术语仅用于描述特定实施例，并且不意在限制本发明的示例 实施例。如本文所使用的，单数形式“一”、“一个”以及“该”意在包括 复数形式，除非上下文明确指
示相反意思。还应当理解术语“包括”、“包 括了”、“包含”、和/或“包含了”当在本文中使用
时，指定所声明的特征、 整数、步骤、操作、单元和/或组件的存在性，并且不排除一个或多
个其他特 征、数量、步骤、操作、单元、组件和/或他们的组合存在性或增加。

[0060] 在下面的描述中提供了特定的细节,以便于对示例实施例的完全理解。 然而,本领域普通技术人员应当理解可以在没有这些特定细节的情况下实现 示例实施例。在其他
实施例中，可以不以非必要的细节来示出众所周知的过 程、结构和技术，以避免使得示例
实施例不清楚。

[0061] 实施例

[0062] 本发明所提供的一种无线基站与线段型地理地物关联矩阵的建立方法， 具体步骤为：

[0063] S1、从移动运营商大数据中提取训练数据TrainData；

[0064] 利用移动运营商管道数据，从移动运营商2G/3G/4G网络的A/Iu/S1-MME 等控制面接口可以采集到信令数据，通过数据转换可以得到位置类数据即每 个用户发生通信事件
所在的基站的位置、时间、时长。

[0065] 从移动运营商2G/3G/4G网络的Gb/IuPS/Gn/S1-u等用户面接口可以采 集到上网数据流，通过数据转换可以得到用户的上网流量数据。随着移动互 联网的发展，大量基于
LBS的APP会上报经纬度，通过解码这些经纬度可以 得到用户的GPS位置，利用各个百度、高
德等地图服务商的坐标转换服务可 以把这些经纬度转换为标准的GPS。

[0066] 通过移动采集数据的控制面和用户面关联方法，比如TMSI，M-TMSI， GTP-ID,APID等方式，经过数据整理，最终可以得到以下记录集。这个记录 集即为关系矩阵的训练数据，
如下所示：

[0067] TrainData＝{(userid,timestart,timeend,cellid,gpsx,gpsy)}

[0068] 其中，userid表示用户ID；timestart表示用户发生通信事件的所在基站的起 始时间；timeend表示用户发生通信事件的所在基站的结束时间；cellid表示手 机ID；gpsx和gpsy
分别表示用户发生通信事件的经纬度。

[0069] S2、从电子地图中提取地物信息，并进行分类，建立PoiData；

[0070] 通过GIS系统获取地区的地理地物，从地物中提取出高速公路、城市道 路、省道、市道、县道等、出入口、停车场等地物信息。把地物信息整理为 以下格式的数据：

[0071] PoiData＝{(Poiname,Poitype,Centerx,Centery,Points)}

[0072] 其中：Poiname表示地物名称，比如：xx大厦；Poitype表示地物类型，比如 楼宇；Centerx和Centery分别表示线段中心点的坐标；Points描述多线段的有 序点的集合。

[0073] S3、利用现行差分法对线段型地物进行分段：

[0074] 高速公路、城市道路往往是线段型地物。而且往往距离比较长。长距离 的地物地理跨度太大，往往与大量的基站小区有关系。因此简单地建立线段 型地物和小区的关系矩
阵，必然不准确，不合理。为了提高准确程度，需要 把这些线段进行划分，按照n米进行划
分，对于城市环境可以按照20～50米 进行划分，对于郊区环境可以按照100～5000米进行
划分。划分方法可以采 用线性差分法。

[0075] S4、通过地理关系关联TrainData和PoiData；

[0076] 对于线段型PoiData，对线段沿法线方向扩展为一个矩形，TrainData和 PoiData的的关联方法为：如果TrainData中的经纬度落在该PoiData对应的 Points中，则这个
TrainData记录可以关联这个区域型PoiData。

[0077] S5、建立基于训练数据的关系矩阵Rcp_train、Rpc_train；

[0078] 对所有TrainData进行计算，然后对地理地物进行合并则得到以下的数据

[0079]

[0080] Cnt为满足关联条件的训练数据的数量。

[0081] 把Cnttrain按行进行归一化得到小区-地物的关系矩阵Rcp_train：

[0082]

[0083] 把Rcp_train进行转置并按行进行归一化得到地物-小区的关系矩阵Rpc_train：

[0084]

[0085] S6、利用传播模型计算每一个无线基站小区到每一个分段点的覆盖强度， 推导出覆盖概率：

[0086] 首先，计算小区到线段中心点的无线路径损耗；

[0087] 采用以下传播模型计算小区到线段中心点的无线接收功率

[0088] Rx＝Tx-(K1+K2log(d)+K3log(hcell)+K4Diff+K5log(d)*log(hcell)+K6(hms)+ K7f(clutter)。

[0089] 其中，K1为与pi年率相关的因子，这里取值为23.5，K2为与距离相关 的因子，这里取值为44.9，K3为天线高度相关因子，这里取值为5.83，K4为与衍射相关的因子，这里可以取
0，K5为与发射天线和距离相关的因子， 这里取值为-6.55，K6接收机高度相关因子，这里可
以取0，K7为地物相关 因子，这里可以取0.

[0090] Tx为小区的发射功率，从基站的工程参数中获得。

[0091] hcell为小区天线高度，从基站的工程参数中获得。

[0092] d为小区到线段中心点的距离。

[0093] 其次，计算小区对线段型地物的覆盖概率：

[0094] 通过上述公式，可以计算每个小区到线段中心点的无线接收功率。

[0095] 为了提高计算速度，可以定义当小区到线段中心点的距离超过小区半径 的3倍时，不参与路径损耗的计算。

[0096] 这样，可以得到各个小区到该线段中心点的接收功率向量(Rx1,Rx2, Rx3…Rxn)，对向量进行归一化，得到归一化的向量，即为小区对该地物的 覆盖概率pline_cov(poii,
cellj)；

[0097] S7、建立基于覆盖的小区-线段型地物的覆盖概率矩阵Rline_cov；

[0098] 遍历每一个线段型地物，即可以得到小区-地物的覆盖概率矩阵：

[0099]

[0100] S8、计算训练矩阵的伴随可信矩阵Bline；

[0101] 由于不是每一个地物都有充足的训练数据，甚至有可能出现训练数据为0的情况，这个时候需要采用覆盖模型对训练模型进行补充和修正。一般来 说基于真实数据的训练
模型的可信度要大于覆盖模型。因此对两个矩阵合并 的时候，要进行权重的设置，由于不
同数据量将产生不同的可信度不一样， 所以针对不同的可信度需要采用不同的权重分配。
由于不同的小区有不同的 训练数据特征，因此需要对每一个小区对应的每一种地物分别
计算融合权重。 这个基于不同小区不同地区的融合权重就是伴随可信矩阵。综合考虑以上
因 素建立针对训练矩阵的伴随可信矩阵B。对于线段型地物，我们把这个矩阵 记为Bline。

[0102] 根据数据分析经验，假设经过这个地物的用户数usercnt超过N1，并且平 均每用户产生的样本数avg_samplecnt超过N2，则认为该样本空间具有统计意 义。则训练矩阵的可信
度最高，获得最大的加权10x。x可以取1，N1可以 取100，N2可以取1000。

[0103] 当样本空间低于上述空间时，可信度的计算方法如下：

[0104]

[0105] 对训练矩阵中的每一个元素进行可信度计算，即得到可信度矩阵Bline：

[0106] Bline＝[bi，j]n*m。

[0107] S9、根据Rcp_train、Rline_cov、Bline三个矩阵联合计算，最终的合并公式为：

[0108] R'line＝Rpc_train.*Bline+Rline_cov

[0109] 对Rl'ine按行进行归一化得到最终的Rline

[0110]

[0111] S10、采用周期性滑动窗口法对关系矩阵Rline进行更新；

[0112] 由于小区的工程参数和无线环境是经常变化的，地理地物也随着城市的 变迁不断地变化，因此需要对Rline矩阵进行不断的更新。本发明综合考虑无 线环境的变化频率和
城市建设的速度，提出以T为周期，以nT为数据滑动 计算窗口的更新方式。如图1所示，比如
T＝1个月，n＝3，则表示关系矩阵 每个月将重新计算一次，每次计算所采用的数据为过去3
个月的数据。

[0113] 最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对 其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通 技术人员应当理解：其
依然可以对前述各实施例所记载的技术方案进行修改， 或者对其中部分或者全部技术特
征进行等同替换；而这些修改或者替换，并 不使相应技术方案的本质脱离本发明各实施例
技术方案的范围。