发明名称 一种车辆排队长度测算的方法及系统
摘要 一种车辆排队长度测算的方法及系统,通过建立城市道路组合交通流模型,将车流量信息分为稳定流和组合流,可准确模拟现实车流量信息。通过对道路进行数字化图形的采集,可使数字化图形的原始地图中的像素点作为道路的坐标,设定采样时间间隔,每隔一段时间进行采样检测所述道路上车辆的离去率,当检测到所述道路上的车辆在道路通行周期内不再前进时,再次在所述采样位置采集所述道路的数字化图形信息,识别出所述数字图形信息上与所述原始地图上不同像素点,获得拥堵发生的起点位置,根据发生拥堵起点距所述检测点的距离,以及检测到车辆在通行周期内不再前进的时刻通过车辆排队长度计算公式可计算出实时车辆排队长度。
申请公布号 CN103985251B 申请公布日期 2016.03.02
申请号 CN201410161695.1 申请日期 2014.04.21
申请人 东南大学;南京市公安局交通管理局 发明人 顾怀中;陆建;王锦尧;叶坚;苏子毅
分类号 G08G1/01(2006.01)I;G08G1/065(2006.01)I;G06T7/00(2006.01)I;G06K9/46(2006.01)I 主分类号 G08G1/01(2006.01)I
代理机构 北京三聚阳光知识产权代理有限公司 11250 代理人 寇海侠
主权项 一种车辆排队长度测算的方法,其特征在于,包括以下步骤:S1:通过建立城市道路组合交通流模型来描述道路上车辆的到达特征,计算所述道路上车辆的到达率,所述到达率为单位时间内进入所述道路的车辆数;进一步包括以下步骤:S11:将所述组合交通流分为稳定流和离散流,所述稳定流由上游路口直行车道放行的连续驶入车辆构成,离散流由上游路口左转和右转车道驶入的车辆构成,将上游路口直行车道开始放行的时间设为周期开始时间,所述道路上车辆的到达率为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><msub><mi>g</mi><mn>1</mn></msub><msub><mi>c</mi><mn>1</mn></msub></mfrac><mrow><mo>(</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>-</mo><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>&lambda;</mi><mn>2</mn></msub></mrow>]]></math><img file="FDA0000826590380000011.GIF" wi="414" he="127" /></maths>其中,λ<sub>1</sub>为稳定流的到达率,λ<sub>2</sub>为离散流的到达率,g<sub>1</sub>为上游路口直行车道放行的持续时间,c<sub>1</sub>为组合交通流周期,且组合交通流周期等于上游路口信号周期;S2:选取道路中某一位置作为采样位置,采集上下游两个路口间无车辆行驶时道路的数字化图形,将获取到的数字化图形作为所述道路的原始地图,所述数字化图形的起始像素点对应道路下游的起始位置,所述数字化图形结束的像素点对应道路上游的结束位置;S3:在所述道路上选择任意位置作为检测点,每隔采样间隔Δt检测所述道路上车辆的离去率,所述离去率为单位时间内离开所述道路的车辆数;S4:所述检测点在检测到所述道路上的车辆在通行周期内不再前进时,再次在所述采样位置采集所述道路的数字化图形信息;S5:将再次获取到的数字化图形,与所述的原始地图进行比对,比对两幅数字化图形的不同,得到拥堵发生的起点距离检测点的距离l<sub>0</sub>,以及判断拥堵发生的位置位于所述检测点的上游或下游;S6:若发生拥堵的位置位于所述检测点的下游,根据所述的距离l<sub>0</sub>,和发生拥堵的位置,以及所述检测点检测到车辆在通行周期内不再前进的时刻t<sub>s</sub>,通过车辆排队长度测算公式<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>L</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mrow><msub><mi>N</mi><mn>1</mn></msub><mo>+</mo><mrow><mo>&lsqb;</mo><mfrac><mrow><mi>t</mi><mo>-</mo><msub><mi>t</mi><mi>s</mi></msub></mrow><msub><mi>c</mi><mn>1</mn></msub></mfrac><mo>&rsqb;</mo></mrow><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover><msub><mi>c</mi><mn>1</mn></msub><mo>+</mo><msub><mi>N</mi><mn>2</mn></msub><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><msub><mi>q</mi><mi>i</mi></msub><mi>&Delta;</mi><mi>t</mi></mrow></mrow><mo>)</mo></mrow><mo>/</mo><msub><mi>k</mi><mi>j</mi></msub><mo>+</mo><msub><mi>l</mi><mn>0</mn></msub></mrow>]]></math><img file="FDA0000826590380000021.GIF" wi="940" he="158" /></maths>计算实时车辆排队长度;若发生拥堵的位置位于所述检测点的上游,根据所述的距离l<sub>0</sub>,和发生拥堵的位置,以及所述检测点检测到车辆在通行周期内不再前进的时刻t<sub>s</sub>,通过车辆排队长度测算公式<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>L</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mrow><msub><mi>N</mi><mn>1</mn></msub><mo>+</mo><mrow><mo>&lsqb;</mo><mfrac><mrow><mi>t</mi><mo>-</mo><msub><mi>t</mi><mi>s</mi></msub></mrow><msub><mi>c</mi><mn>1</mn></msub></mfrac><mo>&rsqb;</mo></mrow><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover><msub><mi>c</mi><mn>1</mn></msub><mo>+</mo><msub><mi>N</mi><mn>2</mn></msub><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><msub><mi>q</mi><mi>i</mi></msub><mi>&Delta;</mi><mi>t</mi></mrow></mrow><mo>)</mo></mrow><mo>/</mo><msub><mi>k</mi><mi>j</mi></msub><mo>-</mo><msub><mi>l</mi><mn>0</mn></msub></mrow>]]></math><img file="FDA0000826590380000022.GIF" wi="931" he="159" /></maths>计算实时车辆排队长度,其中,t为所述检测点检测到车辆在道路通行周期内不发生移动时或没有检测到车辆通过所述检测点后的任意时刻,N<sub>1</sub>和N<sub>2</sub>分别为车辆到达时组合交通流所在周期剩余时间Δc<sub>1</sub>内到达的车辆数和拥堵扩散结束时组合交通流所在周期剩余时间Δc<sub>2</sub>内到达的车辆数;q<sub>i</sub>为第i个采样时间间隔内检测到的车辆离去率,i为所述检测点检测过的时间间隔的序号,n为到任意时刻t时所述检测点检测过的次数;k<sub>j</sub>为路段阻塞密度,取为100‑150辆车/公里;N<sub>1</sub>:当Δc<sub>1</sub>>c<sub>1</sub>‑g<sub>1</sub>时,N<sub>1</sub>=λ<sub>1</sub>(Δc<sub>1</sub>‑c<sub>1</sub>+g<sub>1</sub>)+λ<sub>2</sub>(c<sub>1</sub>‑g<sub>1</sub>),其中λ<sub>1</sub>(Δc<sub>1</sub>‑c<sub>1</sub>+g<sub>1</sub>)为车辆到达时组合交通流所在周期剩余时间Δc<sub>1</sub>内以稳定流形式到达的车辆数,λ<sub>2</sub>(c<sub>1</sub>‑g<sub>1</sub>)为车辆到达时组合交通流所在周期剩余时间Δc<sub>1</sub>以内离散流形式到达的车辆数;当Δc<sub>1</sub><c<sub>1</sub>‑g<sub>1</sub>时,N<sub>1</sub>=λ<sub>2</sub>c<sub>1</sub>,此时,车辆到达时组合交通流所在周期剩余时间Δc<sub>1</sub>内仅包括离散流,且到达的车辆数为λ<sub>2</sub>Δc<sub>1</sub>;N<sub>2</sub>:当Δc<sub>2</sub><g<sub>1</sub>时,N<sub>2</sub>=λ<sub>1</sub>Δc<sub>2</sub>,此时,拥堵扩散结束时组合交通流所在周期剩余时间Δc<sub>2</sub>内仅包括稳定流,且到达的车辆数为λ<sub>1</sub>Δc<sub>2</sub>;当Δc<sub>2</sub>>g<sub>1</sub>时,N<sub>2</sub>=λ<sub>1</sub>g<sub>1</sub>+λ<sub>2</sub>(Δc<sub>2</sub>‑g<sub>1</sub>),其中,λ<sub>1</sub>g<sub>1</sub>为拥堵扩散结束时组合交通流所在周期剩余时间Δc<sub>2</sub>内以稳定流形式到达的车辆数,λ<sub>2</sub>(Δc<sub>2</sub>‑g<sub>1</sub>)为拥堵扩散结束时组合交通流所在周期剩余时间Δc<sub>2</sub>内以离散流形式到达的车辆数。
地址 210096 江苏省南京市玄武区四牌楼2号