CN110288826B - 基于多源数据融合和milp的交通控制子区聚类划分方法 - Google Patents
基于多源数据融合和milp的交通控制子区聚类划分方法 Download PDFInfo
- Publication number
- CN110288826B CN110288826B CN201910438715.8A CN201910438715A CN110288826B CN 110288826 B CN110288826 B CN 110288826B CN 201910438715 A CN201910438715 A CN 201910438715A CN 110288826 B CN110288826 B CN 110288826B
- Authority
- CN
- China
- Prior art keywords
- road
- milp
- snake
- sequence
- traffic control
- 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 44
- 230000004927 fusion Effects 0.000 title claims abstract description 13
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 14
- 238000010219 correlation analysis Methods 0.000 claims abstract description 6
- 241000270295 Serpentes Species 0.000 claims description 58
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000010606 normalization Methods 0.000 claims description 13
- 238000005457 optimization Methods 0.000 claims description 12
- 238000007499 fusion processing Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000011425 standardization method Methods 0.000 claims description 3
- 230000002159 abnormal effect Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000003064 k means clustering Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000012886 linear function Methods 0.000 description 2
- 238000007500 overflow downdraw method Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 206010039203 Road traffic accident Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/231—Hierarchical techniques, i.e. dividing or merging pattern sets so as to obtain a dendrogram
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
- G06F18/251—Fusion techniques of input or preprocessed data
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0125—Traffic data processing
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/052—Detecting movement of traffic to be counted or controlled with provision for determining speed or overspeed
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/065—Traffic control systems for road vehicles by counting the vehicles in a section of the road or in a parking area, i.e. comparing incoming count with outgoing count
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Traffic Control Systems (AREA)
Abstract
基于多源数据融合和MILP的交通控制子区聚类划分方法,首先,根据实际交通环境,选择两类特征参数数据进行相关性分析,并在数据规范化基础上融合两类数据,形成新的组合特征参数;然后,采用层次聚类算法,为所有路段生成可以表征周围交通流态势的有序序列;接着,基于MILP思想进行初次建模和软件求解,获得初次划分结果;最后,针对需要重新分配的路段进行二次MILP建模和求解,获得最终划分结果。本发明适用于大型城市交通网络,基于多源数据融合的组合特征参数可以获得更符合实际交通状况的交通控制子区划分结果,同时MILP模型具有良好的可解性,可以在有效时间内获得全局最优的交通控制子区划分结果。
Description
技术领域
本发明涉及智能交通信号控制的交通控制子区划分方法,交通控制子区用于城市交通信号的干线协调控制和区域协调控制。
背景技术
在城市交通信号控制中,协调控制可以有效提高整个***的通行效率、减少车均停车延误与行车时间。但是在城市交通路网中,各个交叉口和路段在动态交通流特性上均存在不同程度的差异,如果将其作为同一区域执行统一的控制策略,将达不到较好的控制效果,甚至会加重拥堵或导致交通事故。交通控制子区划分主要是将相邻的交叉口或路段划分成若干个交通控制子区进行协调控制,它是城市交通信号控制***的协调控制功能单元。交通控制子区划分是干线协调控制和区域协调控制的基础,同时也是决定协调控制效果好坏的关键技术。交通控制子区的合理划分能够减少交通***控制与交通管理的复杂性,从而提高交通***的效率和可靠性。
在智能交通领域,主要使用聚类方法来进行交通控制子区的划分。这类方法包括两个步骤,即交通流特征参数的选取和聚类划分方法的设计。在交通流特征参数的选取上,现有的交通控制子区划分方法大多数只选择单一来源的特征参数,如行车速度、车流密度等。由于影响交通流的因素很多,这种来源于单一数据的特征参数不能准确和全面地反映实际的交通流态势,从而导致最终的子区划分结果不符合实际需求。虽然也有少数子区划分方法采用了多个来源的特征参数,但这些方法中涉及的很多特征参数,如旅行时间、车流密度、延误、停车次数和排队长度等,在实际的城市交通环境中,却很难进行采集或采集误差很大。
而在聚类划分方法的设计上,现有方法主要有两类方法。第一类方法是将空间限制条件隐藏在数据中,并由此构建路网图,然后采用基于聚类的图形分割方法进行划分。这类方法主要使用传统的聚类方法,如K均值聚类(即K-means聚类)、谱聚类、非负矩阵分解(No n-negative Matrix Factorization,NMF)等。第二类方法是将空间限制条件显式强加在算法的每一个步骤中,这类方法主要包括一些启发式算法和建模优化类方法。但在实际交通路网中,这两类方法都存在着一定的缺陷和不足。
目前,在实际应用中,现有的交通控制子区聚类划分方法存在以下主要问题:1)大多数方法采用单一来源的特征参数,不能准确和全面地表征交通流态势,从而导致划分结果不符合实际需求;2)少数方法虽然采用了多源的特征参数,但存在数据采集困难或采集误差较大等问题;3)传统聚类划分方法(即第一类方法)存在较多不足。K均值聚类存在初始簇中心选择困难和划分不准确等问题,谱聚类的效果则过度依赖于拉普拉斯矩阵的特征值,非负矩阵分解则要求数据具有较好的线性结构;4)其它的聚类划分方法(即第二类方法)也存在一定的缺陷。启发式算法常常只能获得局部最优的子区划分结果,不符合实际的需求;建模优化类方法可以获得最优结果,但往往求解困难,无法应用于大型城市交通网络。
发明内容
本发明要克服现有技术的上述不足之处,针对大型城市交通网络,提供一种基于多源数据融合和MILP(Mixed Integer Linear Programmi ng:混和整数线性规划)的交通控制子区聚类划分方法。
本发明首先根据实际交通流数据采集的方便性和准确性,选取行车速度和车道流量作为交通流特征参数,并对这两类特征参数进行相关性分析,计算它们的皮尔森相关系数(Pearson correlation coefficie nt)。然后,在数据规范化的基础上,采用多源数据融合的方法对这两类特征参数数据进行线性加权数据融合处理,获得新的组合特征参数。接着,在组合特征参数的基础上,对每一条路段进行聚类分析,获得对应的交通流态势和空间相邻路段的分布情况。最后,在聚类分析的基础上,通过两次MILP(Mixed Integer LinearProgramming:混和整数线性规划)建模和软件求解,获得最终的交通控制子区划分结果。本发明充分利用了多源特征参数之间的相关性,通过皮尔森相关系数和数据规范化进行数据融合,形成新的组合特征参数。新的组合特征参数更加准确和全面地表征了实际的交通流态势。同时,本发明也充分利用了MILP特性和现有求解器软件,可以在有效时间内以离线方式计算获得全局最优的交通控制子区划分结果,适用于大型城市交通网络。
本发明是通过以下技术方案达到上述目的,即基于多源数据融合和MILP的交通控制子区聚类划分方法,具体的实施步骤如下:
(1)通过车载GPS(Global Positioning System:全球定位***)、或手机GPS、或北斗***、或第三方公司(如高德、百度)等方式获得行车速度特征参数数据,同时通过电警或卡口的视频摄像头获得车道流量特征参数数据。行车速度是指单位时间内通过某个路段的车辆行驶速度的平均值,单位为km/h,而车道流量是指单位时间内通过某进口车道停车线的车辆数,单位为pcu/h(pcu,passenger car unit,标准小客车单位,也即标准车当量数)。对这两种特征参数进行相关性分析,计算它们的皮尔森相关系数,具体计算公式如下:
其中,ρ表示行车速度和车道流量的皮尔森相关系数,V和Q分别表示所有路段的行车速度和车道流量,cov(V,Q)表示两者的协方差,σV和σQ分别表示V和Q的方差,vi表示路段i的行车速度(单位:km/h),qi表示路段i的车道流量(单位:pch/h),和则分别表示所有路段行车速度和车道流量的平均值,N为交通路网的路段总数。
(2)对步骤(1)采集的交通数据进行数据规范化处理。数据规范化处理可以消除量纲干扰,同时可以避免异常值和极端值的影响。数据规范化处理采用零均值标准化方法,具体计算公式如下:
(3)基于步骤(1)所得到的皮尔森相关系数,对行车速度和车道流量规范化后的数据进行线性加权数据融合处理,获得新的组合特征参数。其中,根据实际的交通流数据采集状况分析,车道流量数据更为准确,因此主特征定为车道流量,辅特征定为行车速度。数据融合的计算公式如下:
(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”序列生成过程中,加入某相邻路段后,组合特征参数方差的计算公式如下:
(5)基于步骤(4)获得所有路段的“snake”序列,采用MILP思想进行初次建模,并借助现有求解器软件进行优化求解,获得交通控制子区的初步划分结果。MILP建模包括决策变量,目标函数和约束条件。初次MILP决策变量共有四个:xij表示路段j是否属于第i条“snake”,即序列si,x′j表示路段j是否至少属于一个子区,x″j表示路段j是否属于多个子区,ti表示“snake”序列si中选中前Ni条路段的加权方差。初次MILP建模的目标是为了获得具有低重复率、高覆盖率和均匀连通的交通控制子区,所以,初次MILP目标函数是最小化所选“snake”集的组合特征参数的加权方差之和,具体公式如下:
其中,Ns是“snake”的数量,var(si)表示序列si的前Ni条路段的方差,ti是决策变量,即ti=Ni×var(si)。var(si)对于所选路段来说是非线性函数,可以通过Ni作为权重系数,将其改为线性模式,满足MILP优化求解的条件。初次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′表示至少被分配到一个交通控制子区的路段数目占所有路段的比重,该值越大越好,理想值为100%,但会导致最后无解,因此一般设置为90%,而a″表示被分配到多个交通控制子区的路段数目占所有路段的比重,取值越小越好,理想值为0%,一般设置为10%。最后,根据初次MILP的决策变量、目标函数和约束条件,借用现有求解器软件,如CPLEX、SOPLEX等,进行优化求解,获得交通控制子区的初步划分结果。
(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|,表示所有路段的平均相异度,其值当路段l和k属于同一个子区i时,取值为1,目标函数是最小化子区内部的相异度d(l,k),当路段l和k属于不同子区时,目标函数则最小化子区之间的相似度二次MILP约束条件可以确保子区内具有较高的相似度及相邻关系,而子区间具有较低的相似度,具体如公式(10)-(20):
0≤yk≤N,xik∈{0,1} (20)
其中,表示路段k的相邻路段的集合,xik表示路段k是否属于第i条“snake”,即序列si。最后,根据二次MILP的决策变量、目标函数和约束条件,借用现有求解器软件(如CPLEX、SOPLEX等)进行优化求解,获得交通控制子区的最终划分结果。
本发明的有益效果在于:(1)本发明充分考虑了实际交通流数据采集的方便性和准确性,采用多源数据融合方法融合了行车速度和车道流量两个特征参数,融合后的组合特征参数更能准确地表征动态的交通流态势,解决了使用单一交通流特征参数所导致的划分不合理问题;(2)本发明充分利用了MILP的特性,在多源数据融合的基础上,借助现有求解器软件,可以在有效时间内以离线方式计算获得全局最优的交通控制子区划分结果,适用于大型城市交通网络。
附图说明
附图1本发明方法的流程图。
具体实施方式
下面结合附图进一步说明本发明的技术方案。
本发明的基于多源数据融合和MILP的交通控制子区聚类划分方法,具体的实施步骤如下:
(1)通过车载GPS(Global Positioning System:全球定位***)、或手机GPS、或北斗***、或第三方公司(如高德、百度)等方式获得行车速度特征参数数据,同时通过电警或卡口的视频摄像头获得车道流量特征参数数据。行车速度是指单位时间内通过某个路段的车辆行驶速度的平均值,单位为km/h,而车道流量是指单位时间内通过某进口车道停车线的车辆数,单位为pcu/h(pcu,passenger car unit,标准小客车单位,也即标准车当量数)。对这两种特征参数进行相关性分析,计算它们的皮尔森相关系数,具体计算公式如下:
其中,ρ表示行车速度和车道流量的皮尔森相关系数,V和Q分别表示所有路段的行车速度和车道流量,cov(V,Q)表示两者的协方差,σV和σQ分别表示V和Q的方差,vi表示路段i的行车速度(单位:km/h),qi表示路段i的车道流量(单位:pch/h),和则分别表示所有路段行车速度和车道流量的平均值,N为交通路网的路段总数。
(2)对步骤(1)采集的交通数据进行数据规范化处理。数据规范化处理可以消除量纲干扰,同时可以避免异常值和极端值的影响。数据规范化处理采用零均值标准化方法,具体计算公式如下:
(3)基于步骤(1)所得到的皮尔森相关系数,对行车速度和车道流量规范化后的数据进行线性加权数据融合处理,获得新的组合特征参数。其中,根据实际的交通流数据采集状况分析,车道流量数据更为准确,因此主特征定为车道流量,辅特征定为行车速度。数据融合的计算公式如下:
(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”序列生成过程中,加入某相邻路段后,组合特征参数方差的计算公式如下:
(5)基于步骤(4)获得所有路段的“snake”序列,采用MILP思想进行初次建模,并借助现有求解器软件进行优化求解,获得交通控制子区的初步划分结果。MILP建模包括决策变量,目标函数和约束条件。初次MILP决策变量共有四个:xij表示路段j是否属于第i条“snake”,即序列si,x′j表示路段j是否至少属于一个子区,x″j表示路段j是否属于多个子区,ti表示“snake”序列si中选中前Ni条路段的加权方差。初次MILP建模的目标是为了获得具有低重复率、高覆盖率和均匀连通的交通控制子区,所以,初次MILP目标函数是最小化所选“snake”集的组合特征参数的加权方差之和,具体公式如下:
其中,Ns是“snake”的数量,var(si)表示序列si的前Ni条路段的方差,ti是决策变量,即ti=Ni×var(si)。var(si)对于所选路段来说是非线性函数,可以通过Ni作为权重系数,将其改为线性模式,满足MILP优化求解的条件。初次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′表示至少被分配到一个交通控制子区的路段数目占所有路段的比重,该值越大越好,理想值为100%,但会导致最后无解,因此一般设置为90%,而a″表示被分配到多个交通控制子区的路段数目占所有路段的比重,取值越小越好,理想值为0%,一般设置为10%。最后,根据初次MILP的决策变量、目标函数和约束条件,借用现有求解器软件,如CPLEX、SOPLEX等,进行优化求解,获得交通控制子区的初步划分结果。
(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|,表示所有路段的平均相异度,其值当路段l和k属于同一个子区i时,取值为1,目标函数是最小化子区内部的相异度d(l,k),当路段l和k属于不同子区时,目标函数则最小化子区之间的相似度二次MILP约束条件可以确保子区内具有较高的相似度及相邻关系,而子区间具有较低的相似度,具体如公式(10)-(20):
0≤yk≤N,xik∈{0,1} (20)
其中,表示路段k的相邻路段的集合,xik表示路段k是否属于第i条“snake”,即序列si。最后,根据二次MILP的决策变量、目标函数和约束条件,借用现有求解器软件(如CPLEX、SOPLEX等)进行优化求解,获得交通控制子区的最终划分结果。
如附图1,本方法首先通过各种检测设备(GPS设备、视频检测器等)或第三方公司获取到行车速度和车道流量的交通流原始数据;然后,对这两类特征参数数据进行相关性分析,即计算它们的皮尔森相关系数,系数越接近1或-1,表示相关性越强,而越接近0则表示相关性越弱;接着,对行车速度和车道流量的交通原始数据进行数据规范化处理,即零均值标准化,消除量纲和异常值的干扰;随后,根据计算获得的皮尔森相关系数,对规范化处理后的数据进行线性加权数据融合处理,其中,车道流量作为主特征参数,皮尔森相关系数为权重系数,从而获得新的组合特征参数,以此来表征路段的交通流态势;之后,在组合特征参数的基础上,对每条路段执行一次层次聚类算法,生成对应的“snake”有序序列;接着,在“snake”序列的基础上进行MILP初次建模,确定相应的决策变量、目标函数和约束条件,实现对所有路段的初次划分。MILP初次建模的目标函数是最小化所选“snake”集的组合特征参数的加权方差之和,在相应的约束条件下,通过现有求解器软件,对初次MILP模型进行求解,可以获得具有低重复率、高覆盖率和均匀且连通的交通控制子区初步划分结果。但经过初步划分后,仍然有部分路段没有被分配或者分配到多个子区中,因此,需要进行二次MILP建模,处理这些需要重新分配的路段,确定二次MILP的决策变量、目标函数和约束条件。二次MILP建模的目标函数是最小化所有子区间的相似度和子区内的相异度之和,在相应的约束条件下,通过现有求解器软件的优化求解,从而获得交通控制子区的最终划分结果。
实验结果表明,本发明提出的交通控制子区划分方法是有效的,相对于单一特征参数,基于多源数据融合的组合特征参数可以获得更符合实际交通状况的交通控制子区划分结果,同时MILP模型具有良好的可解性,可以借助现有求解器软件,在有效时间内以离线方式计算获得全局最优的交通控制子区划分结果,适用于大型城市交通网络。
本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举,本发明的保护范围不应当被视为仅限于实施例所陈述的具体形式,本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。
Claims (1)
1.基于多源数据融合和MILP的交通控制子区聚类划分方法,包括如下步骤:
(1)通过车载GPS、或手机GPS、或北斗***、或第三方公司获得行车速度特征参数数据,同时通过电警或卡口的视频摄像头获得车道流量特征参数数据;对这两种特征参数进行相关性分析,计算它们的皮尔森相关系数,具体计算公式如下:
其中,ρ表示行车速度和车道流量的皮尔森相关系数,V和Q分别表示所有路段的行车速度和车道流量,cov(V,Q)表示两者的协方差,σV和σQ分别表示V和Q的方差,vi表示路段i的行车速度,qi表示路段i的车道流量,和则分别表示所有路段行车速度和车道流量的平均值,N为交通路网的路段总数;
(2)对步骤(1)采集的交通数据进行数据规范化处理;数据规范化处理采用零均值标准化方法,具体计算公式如下:
(3)基于步骤(1)所得到的皮尔森相关系数,对行车速度和车道流量规范化后的数据进行线性加权数据融合处理,获得新的组合特征参数;其中,根据实际的交通流数据采集状况分析,车道流量数据更为准确,因此主特征定为车道流量,辅特征定为行车速度;数据融合的计算公式如下:
(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”序列生成过程中,加入某相邻路段后,组合特征参数方差的计算公式如下:
(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)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910438715.8A CN110288826B (zh) | 2019-05-24 | 2019-05-24 | 基于多源数据融合和milp的交通控制子区聚类划分方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910438715.8A CN110288826B (zh) | 2019-05-24 | 2019-05-24 | 基于多源数据融合和milp的交通控制子区聚类划分方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110288826A CN110288826A (zh) | 2019-09-27 |
CN110288826B true CN110288826B (zh) | 2020-11-13 |
Family
ID=68002531
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910438715.8A Active CN110288826B (zh) | 2019-05-24 | 2019-05-24 | 基于多源数据融合和milp的交通控制子区聚类划分方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110288826B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112085949B (zh) * | 2020-08-13 | 2023-06-27 | 浙江工业大学 | 一种基于交通运行状况异常的路网脆弱性识别、分析与应对方法 |
CN112652164B (zh) * | 2020-12-02 | 2022-12-30 | 北京北大千方科技有限公司 | 一种交通时段划分方法、装置及设备 |
CN113538908B (zh) * | 2021-07-09 | 2022-05-17 | 大连海事大学 | 一种基于组合特征参数的路况分区*** |
CN113851009B (zh) * | 2021-10-11 | 2023-04-07 | 武汉理工大学 | 一种城市干道时空划分方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012252377A (ja) * | 2011-05-31 | 2012-12-20 | Nec Corp | 交通情報生成システム、交通情報生成方法及びプログラム |
CN103813355A (zh) * | 2014-02-21 | 2014-05-21 | 厦门大学 | 一种分布式网络中协同同步的异常点的识别方法 |
CN103886747A (zh) * | 2014-03-14 | 2014-06-25 | 浙江大学 | 路段交通运行相似度度量方法 |
CN104050809A (zh) * | 2014-06-03 | 2014-09-17 | 北京航天福道高技术股份有限公司 | 交通状况自动分析预测方法 |
CN106781503A (zh) * | 2017-01-22 | 2017-05-31 | 百度在线网络技术(北京)有限公司 | 用于监测驾驶行为的方法和装置 |
CN108762201A (zh) * | 2018-04-18 | 2018-11-06 | 南京工业大学 | 一种基于Pearson相关性的大***图论分解方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103325247B (zh) * | 2012-03-19 | 2015-07-01 | ***通信集团辽宁有限公司 | 一种交通信息处理方法及*** |
CN105118289A (zh) * | 2015-06-17 | 2015-12-02 | 河南理工大学 | 一种交通态势评估方法 |
US10332403B2 (en) * | 2017-01-04 | 2019-06-25 | Honda Motor Co., Ltd. | System and method for vehicle congestion estimation |
-
2019
- 2019-05-24 CN CN201910438715.8A patent/CN110288826B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012252377A (ja) * | 2011-05-31 | 2012-12-20 | Nec Corp | 交通情報生成システム、交通情報生成方法及びプログラム |
CN103813355A (zh) * | 2014-02-21 | 2014-05-21 | 厦门大学 | 一种分布式网络中协同同步的异常点的识别方法 |
CN103886747A (zh) * | 2014-03-14 | 2014-06-25 | 浙江大学 | 路段交通运行相似度度量方法 |
CN104050809A (zh) * | 2014-06-03 | 2014-09-17 | 北京航天福道高技术股份有限公司 | 交通状况自动分析预测方法 |
CN106781503A (zh) * | 2017-01-22 | 2017-05-31 | 百度在线网络技术(北京)有限公司 | 用于监测驾驶行为的方法和装置 |
CN108762201A (zh) * | 2018-04-18 | 2018-11-06 | 南京工业大学 | 一种基于Pearson相关性的大***图论分解方法 |
Non-Patent Citations (2)
Title |
---|
Clustering of heterogeneous networks with directional flows based on "Snake" similarities;Mohammadreza;《Transportation Research PartB:Methodological》;20160930;第254页第1-3段,第255页第1-4段 * |
Dynamic clustering and propagation of congestion in heterogeneously congested urban traffic networks;Mohammadreza Saeedmanesh;《Transportation Research PartB》;20171130;第196页第4-5段,第197页第1-3段,第198页第1-2段,第199页第1-4段,第200页第1段以及第201页第1-3段,第198页表1以及第200页表2 * |
Also Published As
Publication number | Publication date |
---|---|
CN110288826A (zh) | 2019-09-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110288826B (zh) | 基于多源数据融合和milp的交通控制子区聚类划分方法 | |
CN110288825B (zh) | 基于多源数据融合和snmf的交通控制子区聚类划分方法 | |
CN107610464B (zh) | 一种基于高斯混合时间序列模型的轨迹预测方法 | |
CN103295414B (zh) | 一种基于海量历史gps轨迹数据的公交车到站时间预测方法 | |
CN107346460A (zh) | 基于智能网联***下前车运行信息的未来工况预测方法 | |
CN107909247A (zh) | 一种基于空间层次贝叶斯模型的城市宏观道路交通安全影响因素分析方法 | |
CN114463972B (zh) | 一种基于etc门架通信数据的路段区间交通分析预测方法 | |
Wang et al. | Energy consumption characteristics based driving conditions construction and prediction for hybrid electric buses energy management | |
CN109887279B (zh) | 一种交通拥堵预测方法及*** | |
Xie et al. | Impact of connected and autonomous vehicle technology on market penetration and route choices | |
CN110930693B (zh) | 一种用于道路断面的在线短时交通流预测方法 | |
CN108335483A (zh) | 交通拥堵扩散路径的推断方法及其*** | |
Wang et al. | Fusing heterogeneous traffic data by Kalman filters and Gaussian mixture models | |
Shen et al. | An attention-based digraph convolution network enabled framework for congestion recognition in three-dimensional road networks | |
Kodi et al. | Safety evaluation of an adaptive signal control technology using an empirical Bayes approach | |
Schnitzler et al. | Combining a Gauss-Markov model and Gaussian process for traffic prediction in Dublin city center. | |
Zhao et al. | Unifying Uber and taxi data via deep models for taxi passenger demand prediction | |
CN108776650A (zh) | 一种基于混合高斯模型的非机动车道服务水平评价方法 | |
Sobreira et al. | Disaggregated traffic conditions and road crashes in urban signalized intersections | |
Chen et al. | Traffic signal optimization control method based on adaptive weighted averaged double deep Q network | |
Xia et al. | A distributed EMDN-GRU model on Spark for passenger waiting time forecasting | |
CN113553350A (zh) | 一种相似演化模式聚类及动态时区划分的交通流分区模型 | |
DE102018201787A1 (de) | Verfahren und System zum Optimieren und Vorhersagen einer Verkehrssituation | |
Del Ser et al. | Introduction to the special issue on online learning for big-data driven transportation and mobility | |
Cussigh et al. | An all-electric alpine crossing: time-optimal strategy calculation via fleet-based vehicle data |
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 |