发明名称 基于空间聚类分析的城市路网交通小区动态划分方法
摘要 一种城市路网交通小区动态划分方法,从已有的出租车GPS数据出发,对路网的交通状态参数进行估计,考虑城市交通的OD特性,采用路段拥挤度与路段关联度进行空间统计分析,分析路段之间的空间自相关模式,实现对交通小区的自动划分。本发明的有益效果经过交通仿真的实际检验:当出租车辆在路网交通量中的覆盖率到达5%时,路网交通状态的估计正确率在80%以上,在此基础上,诱导车辆避开交通拥堵区域,采取必要的控制管理策略,减少区域拥堵的规模和拥堵时间,加快拥堵的消散。本发明中采用的方法和技术简单易行,运行条件容易满足,易于在全国大中型城市中推广应用。
申请公布号 CN101751777B 申请公布日期 2011.11.16
申请号 CN200810203852.5 申请日期 2008.12.02
申请人 同济大学 发明人 杨晓光;李晓丹;储浩
分类号 G08G1/00(2006.01)I;G08G1/01(2006.01)I 主分类号 G08G1/00(2006.01)I
代理机构 上海智信专利代理有限公司 31002 代理人 吴林松
主权项 1.一种城市路网交通小区动态划分方法,其特征在于:从已有的出租车GPS数据出发,对路网的交通状态参数进行估计,考虑城市交通的OD特性,采用路段拥挤度与路段关联度进行空间统计分析,分析路段之间的空间自相关模式,实现对交通小区的自动划分;动态采集出租车GPS数据,通过数据预处理确定车辆在路网上的逻辑位置;每隔一个估计周期做一次路网交通状态参数估计工作,参数估计结果保存在数据库中,通过空间聚类分析算法实现对交通状态参数和空间关联矩阵的聚类分析而得到对于交通控制和交通诱导不同需求下的若干个交通小区;用实际计算与历史数据结合估计路段平均行程车速,用路段延误时间计算路段拥挤度,结合路段的空间关联矩阵,采用空间聚类分析技术进行交通小区的自动划分;具体包括以下步骤:(1)利用实际GPS数据计算路段行程车速,与历史数据结合,采用加权平均法计算路段当前的平均行程车速;(2)通过计算路段的延误时间,与按照最大限速行驶所用时间的比值,计算路段的当前拥挤度指标;(3)通过分析GIS的空间拓扑关系,计算任意两两路段间的关联关系,量化为空间关联矩阵;(4)采用空间自相关分析Moran’s I指数对路段拥挤度观测值/计算值相似性进行分析,得出路网拥挤度的空间分布模式;(5)采用专家知识库的方法对交通小区进行修正和优化;(6)定时对历史数据进行更新,采用二分平均法减少远期数据对近期数据的影响;对于路段L<sub>j</sub>,其实时平均行程车速<img file="FSB00000563479800011.GIF" wi="49" he="54" />的算法为:<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>V</mi><mi>j</mi><mi>r</mi></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>S</mi><mi>i</mi></msub><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>T</mi><mi>i</mi></msub></mrow></math>]]></maths>其中N为经过路段L<sub>j</sub>的车辆数,S<sub>i</sub>为某车在路段L<sub>j</sub>的行驶距离,T<sub>i</sub>为该车在路段L<sub>j</sub>的行程时间;对路段L<sub>j</sub>,若其平均行程车速为<img file="FSB00000563479800013.GIF" wi="72" he="60" />历史行程车速为<img file="FSB00000563479800014.GIF" wi="73" he="67" />则其行程车速V<sub>j</sub>采用加权调和平均法计算获得:<maths num="0002"><![CDATA[<math><mrow><msub><mi>V</mi><mi>j</mi></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mi>&alpha;</mi><mo>/</mo><msubsup><mi>V</mi><mi>j</mi><mi>r</mi></msubsup><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo></mrow><mo>/</mo><msubsup><mi>V</mi><mi>j</mi><mi>h</mi></msubsup></mrow></mfrac></mrow></math>]]></maths>其中,α∈[0,1],为可调系数,在系统运行的初始阶段历史数据记录少,α取值偏大,取0.8,当系统运行一定时间后,将α取值调低,取0.6,到0.5,长期运行后,根据实测数据将α调到合适值;路段拥挤度λ的计算方法为:<maths num="0003"><![CDATA[<math><mrow><mi>&lambda;</mi><mo>=</mo><mfrac><mrow><msub><mi>T</mi><mi>k</mi></msub><mo>-</mo><msub><mi>T</mi><mrow><mi>k</mi><mn>0</mn></mrow></msub></mrow><msub><mi>T</mi><mrow><mi>k</mi><mn>0</mn></mrow></msub></mfrac></mrow></math>]]></maths>其中:<maths num="0004"><![CDATA[<math><mrow><msub><mi>T</mi><mi>k</mi></msub><mo>=</mo><mfrac><msub><mi>L</mi><mi>k</mi></msub><msub><mover><mi>v</mi><mo>&OverBar;</mo></mover><mi>k</mi></msub></mfrac><mo>=</mo><mfrac><msub><mi>L</mi><mi>k</mi></msub><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><msub><mover><mi>v</mi><mo>&OverBar;</mo></mover><mi>ki</mi></msub><mo>/</mo><mi>n</mi></mrow></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>T</mi><mrow><mi>k</mi><mn>0</mn></mrow></msub><mo>=</mo><mfrac><msub><mi>L</mi><mi>k</mi></msub><msub><mover><mi>v</mi><mo>&OverBar;</mo></mover><mrow><mi>k</mi><mn>0</mn></mrow></msub></mfrac></mrow></math>]]></maths>式中,T<sub>k</sub>为路段k上车辆的实际平均行程时间;T<sub>k0</sub>为路段k上车辆按照最大限速行驶的行程时间;L<sub>k</sub>代表路段k的实际长度;<img file="FSB00000563479800025.GIF" wi="36" he="48" />是路段k上车辆的平均行驶速度;<img file="FSB00000563479800026.GIF" wi="51" he="48" />是路段k上第i辆车的平均行驶速度;<img file="FSB00000563479800027.GIF" wi="54" he="48" />表示路段k上允许或设计的车辆的最大行驶速度;n为路段k上的车辆数;对于任意两路段观测值是否在空间上存在相似性,采用空间自相关统计量Moran’s I按下式计算:<maths num="0006"><![CDATA[<math><mrow><mi>I</mi><mo>=</mo><mfrac><mi>n</mi><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>&CenterDot;</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>W</mi><mi>ij</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>-</mo><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>W</mi><mi>ij</mi></msub></mrow></mfrac></mrow></math>]]></maths>对路段观测值在空间上是否存在空间自相关,根据标准化以后的Moran’s I值或z值:<maths num="0007"><![CDATA[<math><mrow><msub><mi>Z</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><mi>I</mi><mo>-</mo><mi>E</mi></mrow><msqrt><mi>Var</mi><mrow><mo>(</mo><mi>I</mi><mo>)</mo></mrow></msqrt></mfrac></mrow></math>]]></maths>计算,Moran’s I的取值范围为(-1,1),在给定显著性水平下,当Moran’s I显著为正时,表示观测值之间存在显著的正相关,呈现高高积聚或低低积聚分布格局;当Moran’s I显著为负时,表示观测值之间存在显著的负相关,高低相异,呈现空间分散格局;当Moran’s I接近期望值时,表明不存在空间自相关,观测值在空间上随机排列,生成具有不同交通拥挤度的若干个交通小区。
地址 200092 上海市杨浦区四平路1239号