1.基于多源数据融合和MILP的交通控制子区聚类划分方法，包括如下步骤：

(1)通过车载GPS、或手机GPS、或北斗系统、或第三方公司获得行车速度特征参数数据，同时通过电警或卡口的视频摄像头获得车道流量特征参数数据；对这两种特征参数进行相关性分析，计算它们的皮尔森相关系数，具体计算公式如下：其中，ρ表示行车速度和车道流量的皮尔森相关系数，V和Q分别表示所有路段的行车速度和车道流量，cov(V，Q)表示两者的协方差，σV和σQ分别表示V和Q的方差，vi表示路段i的行车速度，qi表示路段i的车道流量，和 则分别表示所有路段行车速度和车道流量的平均值，N为交通路网的路段总数；

(2)对步骤(1)采集的交通数据进行数据规范化处理；数据规范化处理采用零均值标准化方法，具体计算公式如下：

其中，i＝1，2，...，N，xi表示路段i的某个特征参数，即行车速度或车道流量，μ和σ分别表示这个特征参数的平均值和方差，而 则表示xi进行零均值标准化后的值；

(3)基于步骤(1)所得到的皮尔森相关系数，对行车速度和车道流量规范化后的数据进行线性加权数据融合处理，获得新的组合特征参数；其中，根据实际的交通流数据采集状况分析，车道流量数据更为准确，因此主特征定为车道流量，辅特征定为行车速度；数据融合的计算公式如下：其中，i＝1，2，...，N，ci表示路段i的新的组合特征参数， 表示路段i经过规范化处理后的行车速度， 则表示路段i经过规范化处理后的车道流量；

(4)基于步骤(3)所得到的组合特征参数，对每条路段执行一次层次聚类算法，为每条路段生成一个由路网中所有路段组成的有序序列，即“snake”序列si，其中i＝1，2，...，N；

一个序列即是一条“snake”，这条“snake”代表了起始元素对应路段周围的交通流态势分布情况，以及相邻路段的分布情况；“snake”序列si的生成过程：首先选择一条没有执行过该算法的路段i，获得对应的组合特征值ci，将其作为首个元素加入对应的“snake”序列si；然后，获取序列si中已有路段的所有相邻路段，假设这些相邻路段逐一加入序列si中，并计算加入后的组合特征参数方差，比较所有这些相邻路段逐一加入后的方差，选择方差最小的路段作为新的元素加入序列si中；重复迭代，直到所有路段都加入到序列si中，则路段i的层次聚类算法执行结束，序列si生成完毕；对每一条路段i都重复这个生成过程，直至生成所有的“snake”，即序列s1，s2，...，sN；“snake”序列生成过程中，加入某相邻路段后，组合特征参数方差的计算公式如下：其中， 和 是“snake”序列si前j条路段关于组合特征参数的方差值和平均值，cj是第j条路段的组合特征参数值；

(5)基于步骤(4)获得所有路段的“snake”序列，采用MILP思想进行初次建模，并借助现有求解器软件进行优化求解，获得交通控制子区的初步划分结果；MILP建模包括决策变量，目标函数和约束条件；初次MILP决策变量共有四个：xij表示路段j是否属于第i条“snake”，即序列si，x′j表示路段j是否至少属于一个子区，x″j表示路段j是否属于多个子区，ti表示“snake”序列si中选中前Ni条路段的加权方差；初次MILP目标函数是最小化所选“snake”集的组合特征参数的加权方差之和，具体公式如下：其中，Ns是“snake”的数量，var(si)表示序列si的前Ni条路段的方差，ti是决策变量，即ti＝Ni×var(si)；初次MILP约束条件可以保证所有交通控制子区内部路段具有较高的相似度且具有相邻关系，具体如公式(1)-(9)：xij∈{0，1}，x′j∈{0，1}，x″j∈{0，1}，ti∈R     (9)

其中，表示路网中所有路段的集合，N和Ns分别表示路段的数量和“snake”的数量，Nc和Nmin则分别表示交通控制子区数量以及子区最少包含的路段数量，Ri(j)表示第i条“snake”的第j条路段， 则表示路段j在第i条“snake”中的位置序号，var(S{i，j})表示第i条“snake”的前j条路段的组合特征参数方差值，a′表示至少被分配到一个交通控制子区的路段数目占所有路段的比重，一般取值为90％，而a″表示被分配到多个交通控制子区的路段数目占所有路段的比重，一般取值为10％；最后，根据初次MILP的决策变量、目标函数和约束条件，借用现有求解器软件进行优化求解，获得交通控制子区的初步划分结果；

(6)基于步骤(5)获得的交通控制子区初步划分结果，采用MILP思想进行二次建模和优化求解，把部分分配到多个子区的路段和没有分配到子区的路段进行重新分配，从而获得交通控制子区的最终划分结果；二次MILP的决策变量有四个：yj表示路段j与交通控制子区中心路段的距离，xijk表示路段k是否能通过相邻路段j与第i条“snake”的中心路段相连，表示路段k和l是否属于同一个交通控制子区i，xiok表示路段k是否是交通控制子区i的中心路段；二次MILP目标函数是为了最小化所有子区间的相似度和子区内的相异度之和，具体公式如下：其中，d(l，k)表示路段l和路段k之间的相异度，Dl和Dk分别是路段l和路段k的组合特征参数值，路段之间的相异度d(l，k)＝|Dl-Dk|，表示所有路段的平均相异度，其值二次MILP约束条件可以确保子区内具有较高的相似度及相邻关系，而子区间具有较低的相似度，具体如公式(10)-(20)：

0≤yk≤N，xik∈{0，1}       (20)

其中， 表示路段k的相邻路段的集合，xik表示路段k是否属于第i条“snake”，即序列si；最后，根据二次MILP的决策变量、目标函数和约束条件，借用现有求解器软件进行优化求解，获得交通控制子区的最终划分结果。