基于基站定位数据的商圈分析

实验背景随着个人手机和网络的普及，手机已经基本成为所有人必须持有的工具。

根据手机信号再地理空间的覆盖情况结合时间序列的手机定位数据可以完整的还原人群的现实活动轨迹从而得到人口空间分布于活动联系的特征信息。

商圈是现代市场中的重要企业活动空间，商圈划分的目的之一是为了研究潜在的顾客分布，以制定适宜的商业对策。

本次数据，是由通信运营商提供的，特定接口解析得到的用户定位数据。

实验目标对用户的历史定位数据，采用数据挖掘技术，对基站进行分群对不同的商圈分群进行特征分析，比较不同商圈类别的价值，选择合适区域进行针对性的营销活动实验分析方法与过程通过数据是由通信运营商提供，可以推测手机用户在使用短信业务、通话业务和上网业务等信息时产生的定位数据。

由于数据中是由不同基站作为位置的辨别，可以将每个基站当做一个”商圈“，再通过归纳基站范围的人口特征，运用聚类算法，识别不同类别的基站范围，即可等同识别不同类别的商圈。

衡量区域人口特征，可以从人流量和人均停留时间的角度进行分析。

分析流程1）从移动通信运营商提供的特定接口上解析、处理、并滤除用户属性后得到用户定位数据。

2）以单个用户为例，进行数据探索分析，研究在不同基站的停留时间，并进一步地进行预处理，包括数据规约和数据变换。

3）利用步骤2）形成的已完成数据预处理的建模数据，基于基站覆盖范围区域的人流特征进行商圈聚类，对各个商圈分群进行特征分析，选择合适的区域进行运营商的促销活动。

数据抽取分析数据抽取从移动通信运营商提供的特定接口上解析、处理和过滤后，获得位置数据。

时间设定为，2014-1-1为开始时间，2014-6-30为结束时间作为分析的观测窗口，抽取窗口内某市某区域的定位数据形成建模数据。

数据分析为了便于观察数据，先提取用户的ID即EMASI号为”55555“的用户再2014年1月1号的定位数据。

观察数据，可以发现，两条数据可能是同一基站的不同时间。

从表中可以看到，用户在2014年1月1日00:31:48处于36908基站的范围，下一个记录是用户在2014年1月1日00:53:46处于36908基站的范围，这表明了用户从00:31:48到00:53:46都是处于36908基站，共停留了21分58秒，并且在00:53:46进入了36902基站的范围。

判断用户在每个基站的停留时间需要依靠当前记录与下一条记录的对比，发生变化则意味着用户所在基站发生改变，可以记录在上一个基站停留的时间。

数据预处理数据规约原始数据的属性较多，由我们的挖掘目标，网络类型、LOC编号和信令类型这三个属性没有作用，剔除。

衡量用户停留时间，没有必要精确到毫秒，故一同删除。

在计算用户的停留时间，只计算两条记录的时间差，为了减少数据维度，把年月日合并为日期，时分秒合并为时间，得到数据。

数据变换挖掘的目标是寻找高价值的商圈，需要根据用户的定位数据，提取出基站范围内区域的人流特征，如人均停留时间和人流等。

高价值商圈，在人流特征上有人流量大和人均停留时间长的特点。

写字楼的上班族在白天所处基站范围固定，时间也较长，人流量也大。

居住区，也有基站范围固定，时间长，人流量大的特点。

所以，单纯的停留时间无法判断商圈类别。

现代社会工作，以一周为一个工作小周期，分为工作日和周末。

一天中，分为上班时间和下班时间。

综上所述，设计人流特征的四个指标，工作日上班时间人均停留时间、凌晨人均停留时间、周末人均停留时间和日均人流量。

工作日上班时间人均停留时间，意思是所有用户在上班时间9：0018：00处在该基站范围内的平均时间凌晨人均停留时间，意思是所有用户在凌晨时间00：0007：00处在该基站范围的平均时间周末人均停留时间，如上类推。

日均人流量，指的是平均每天曾经在该基站范围内的人数。

这四个指标的计算，直接从原始数据计算比较复杂，需要先处理成中间数据，再从中计算得出四个指标。

对于基站1，有以下公式，再带入所有基站，得出结果。

中间过程数据的计算以单个用户在一天里的定位数据为基础，计算在各个基站范围下的工作日上班时间停留时间、凌晨停留时间、周末停留时间和是否处于基站范围。

假设原始数据所有用户在观测窗口期间（L天）曾经经过的基站有N个，用户有M个，用户i在j天经过的基站有num1、num2和num3，则用户i在j天在num1基站的工作日上班时间停留时间为weekday_num1ij，在num2基站的工作日上班时间停留时间为weekday_num2ij，在num3基站的工作日上班时间停留时间为weekday_num3ij；在num1基站的凌晨停留时间为night_num1ij，在num2基站的凌晨停留时间为night_num2ij，在num3基站的凌晨停留时间为night_num3ij；在num1基站的周末停留时间为weekend_num1ij。

在num1基站是否停留为stay_num1ij。

对于未停留的其他基站，工作日上班时间停留时间、凌晨停留时间、周末停留时间和是否处于基站范围的值均为0。

由于各个属性之间的差异较大。

为了消除数量级数据带来的影响，在聚类之前，需要进行离差标准化处理，离差标准化处理的代码如下，得到建模的样本数据。

Data <- read.csv("./data/business_circle.csv", header = TRUE, encoding = 'utf-8')colnames(Data) <- c("number", "x1", "x2", "x3", "x4")attach(Data)y1 = (x1 - min(x1)) / (max(x1) - min(x1))y2 = (x2 - min(x2)) / (max(x2) - min(x2))y3 = (x3 - min(x3)) / (max(x3) - min(x3))y4 = (x4 - min(x4)) / (max(x4) - min(x4))standardized <- data.frame(Data[, 1], y1, y2, y3, y4)write.csv(standardized, "./tmp/standardizedData.csv", s = TRUE)模型构建——层次聚类算法层次聚类层次聚类方法对给定的数据集进行层次的分解，直到某种条件满足为止。

在已经得到距离值之后，元素间可以被联系起来。

通过分离和融合可以构建一个结构。

传统上，表示的方法是树形数据结构，层次聚类算法，要么是自底向上聚集型的，即从叶子节点开始，最终汇聚到根节点；要么是自顶向下分裂型的，即从根节点开始，递归的向下分裂。

1.凝聚层次聚类：AGNES算法(自底向上)首先将每个对象作为一个簇，然后合并这些原子簇为越来越大的簇，直到某个终结条件被满足2.分裂层次聚类：DIANA算法(自顶向下)首先将所有对象置于一个簇中，然后逐渐细分为越来越小的簇，直到达到了某个终结条件。

数据进行预处理以后，已经形成了建模数据。

这次聚类，采用层次聚类算法，对建模数据进行基于基站数据的商圈聚类，画出谱系聚类图，代码如下。

Data <- read.csv("./data/standardized.csv", header = FALSE)# colnames(Data) <- c("x1", "x2", "x3"," x4")attach(Data)dist <- dist(Data, method = 'euclidean')hc1 <- hclust(dist, "ward.D2")plot(hc1, hang = -1)detach(Data)library(ggplot2)Data <- read.csv("./data/standardized.csv", header = FALSE)attach(Data)dist <- dist(Data, method = 'euclidean')hc1 <- hclust(dist, "ward.D2")plot(hc1)# 分成三类re1 <- rect.hclust(hc1, k = 3)a <- re1[[2]]b <- re1[[3]]c <- re1[[1]]# 商圈类别1matrix <- Data[a,]d <- dim(matrix)y <- as.numeric(t(matrix))row <- factor(rep(1:d[1], each = d[2]))x <- rep(1:d[2], times = d[1])data <- data.frame(y = y, x = x, row = row)ggplot(data = data, aes(x = x, y = y, group = row)) + geom_line() +scale_x_continuous(breaks = c(1, 2, 3, 4),labels = c("工作日人均停留时间", "凌晨人均停留时间","周末人均停留时间", "日均人流量")) +labs(title = "商圈类别1", x = "", y = "")# 商圈类别2matrix <- Data[b, ]d <- dim(matrix)y <- as.numeric(t(matrix))row <- factor(rep(1:d[1], each = d[2]))x <- rep(1:d[2], times = d[1])data <- data.frame(y = y, x = x, row = row)ggplot(data = data, aes(x = x, y = y, group = row)) + geom_line() +scale_x_continuous(breaks = c(1, 2, 3, 4),labels = c("工作日人均停留时间", "凌晨人均停留时间","周末人均停留时间", "日均人流量")) + labs(title="商圈类别2", x="", y="")# 商圈类别3matrix <- Data[c,]d <- dim(matrix)y <- as.numeric(t(matrix))row <- factor(rep(1:d[1], each = d[2]))x <- rep(1:d[2], times = d[1])data <- data.frame(y = y, x = x, row = row)ggplot(data = data, aes(x = x, y = y, group = row)) + geom_line() +scale_x_continuous(breaks = c(1, 2, 3, 4),labels = c("工作日人均停留时间", "凌晨人均停留时间","周末人均停留时间", "日均人流量")) + labs(title = "商圈类别3", x = "", y = "")模型分析针对聚类结果，按不同类别画出了3个特征的折线图：商圈类别1，工作日人均停留时间、凌晨人均停留时间都很低，周末人均停留时间中等，日均人流量极高，这符合商业区的特点。