CN104851286B - 一种利用公交车辆gps数据动态评估道路交通状态的方法 - Google Patents
一种利用公交车辆gps数据动态评估道路交通状态的方法 Download PDFInfo
- Publication number
- CN104851286B CN104851286B CN201410781208.1A CN201410781208A CN104851286B CN 104851286 B CN104851286 B CN 104851286B CN 201410781208 A CN201410781208 A CN 201410781208A CN 104851286 B CN104851286 B CN 104851286B
- Authority
- CN
- China
- Prior art keywords
- road
- gps
- section unit
- mrow
- point
- 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.)
- Expired - Fee Related
Links
Classifications
-
- 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
- G08G1/0133—Traffic data processing for classifying traffic situation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Computer Networks & Wireless Communication (AREA)
- Traffic Control Systems (AREA)
- Navigation (AREA)
Abstract
本发明涉及一种利用公交车辆GPS数据动态评估道路交通状态的方法,属于智能交通领域。配有GPS设备的公交车在道路上行驶时,记录车辆GPS数据,数据含有车辆位置信息和产生的时间信息。计算时,首先将数据中的GPS点转换为空间坐标点,并与行驶道路匹配,获取该条道路在计时周期中,一定范围内的GPS点;然后将这部分GPS点与已知位置信息的公交站点进行匹配,剔除公交站点周边的无效GPS点。之后,利用有效GPS点进行方位角计算,通过与匹配道路方位角的比较,确定公交车辆行车方向。最后,根据有效GPS点位置与时间信息,计算车辆平均行驶速度,利用速度补偿模型修正为道路交通运行速度,与交通状态区间分级进行比较,得到道路交通运行状态。
Description
技术领域
本发明涉及智能交通领域中的交通信息处理技术。可以利用公交车载GPS设备记录的GPS定位数据,分析道路交通运行状态。
背景技术
我国经济和社会的发展,带来机动车保有量高速快速增长,交通基础设施的供给无法满足机动车出行需求,造成严重的交通拥堵问题。提前发现交通拥堵路段,为出行人员提供拥堵路径的信息,可以更好的引导出行人员选择合理出行路径。目前利用公交车载GPS数据进行道路交通状态评估,主要是利用路段边界两侧的GPS数据估计路段边界时刻,从而获取单车路段行程时间,根据路段长度计算得到单车路段行程速度,再根据道路等级、交通流量确定道路速度分级标准,进而到得交通运行状态。但该方法没有考虑公交车辆在站台停靠时对速度的影响,实际中,公交车辆在站台停靠时,乘客上下车花费的时间对平均速度影响很大,进行平均速度估计时应给予考虑。
发明内容
本发明提供的是一种根据公交车辆GPS数据,估计道路交通运行状态的方法。方法如下:
根据本发明,配有车载GPS设备的公交车辆在道路上行驶时,实时记录车辆GPS数据,该数据含有车辆GPS点位置信息和GPS点位产生的时间信息。计算时,首先将信息中的GPS点转换为空间坐标点,并与车辆运行道路进行匹配,获取该条道路在计时周期内的有效范围内的GPS点,计时周期可根据信息发布的要求而定;然后将这部分GPS点与已知位置信息的公交站点进行匹配,剔除公交站点周边的无效GPS点,保留有效GPS点。之后,利用有效GPS点进行方位角计算,通过与匹配道路方位角的比较,确定公交车辆行车方向。最后,根据有效GPS点位信息与时间信息,分方向计算车辆运行速度,利用速度补偿模型修正为道路交通运行速度。
具体方法和步骤如下:
步骤1:划分路段单元。把道路分成上行和下行两个方向。划分路段单元是指为计算道路交通运行速度,把道路的上下行方向划分成多个独立的路段单元,相同的路段单元,认为交通流速度一致,根据交通流的方向计算路段单元的方位角,以正北方向为0度角,路段单元的方位角为路段单元起点和终点的连线与正北方向的顺时针夹角,路段单元i的方位角记为
步骤2:建立路段单元与公交车站对应关系表。根据路段单元和公交车站地理空间位置,建立两者的归属关系,即公交车站属于哪一个路段单元,公交车站j归属于路段单元i,其坐标信息记为
步骤3:过滤路段单元周边的GPS点。过滤路段单元周边dlink范围内的GPS坐标点,dlink取值依据道路宽度,取10-80米,若GPS坐标点n属于路段单元i,则记为包含GPS坐标和时间信息,记为
步骤4:把过滤的路段单元周边GPS数据,按公交车辆编号BIDk、时间t排序,构成路段单元i周边GPS数据集合
步骤5:为消除公交停靠站对速度估计的影响,把公交车站周边dbus范围内的GPS点从集合中剔除,dbus取值依据公交站台的长度,取10-50米,构成集合
步骤6:以BIDk中时间最小的点为起点时间最大的点为终点计算所构成线段的方位角,记为比较与的绝对值大小,若时,认为BIDk的GPS数据是有效数据,否则认为是无效数据,无效数据应剔除。Athreshold的取值为0-45,取值越小,结果精度越高。构成集合
步骤7:在集合中,剔除仅有一条记录的BIDk GPS数据,构成集合
步骤8:计算路段单元i公交车辆平均运行速度集合中,取相同车辆编号最大时间和最小时间的GPS点,计算两点的时间差值和距离差值,进而计算车辆的运行速度。多辆公交车的平均运行速度即是该路段单元公交的平均公交运行速度
步骤9:利用速度补偿模型进行速度补偿,得到道路平均运行速度。道路交通运行速度与公交平均运行速度存在线性关系(见附图2),速度补偿模型为:
步骤10:根据速度的区间分级,把道路平均运行速度落在对应的区间,该区间表示的状态即是道路运行状态。
本发明具有以下有益效果:
该发明利用在公交车上普遍安装的GPS数据采集设备,采集公交车辆的位置点信息,计算公交车的平均运行速度,通过道路交通运行速度大小与公交运行速度存在相关性的特点,估计道路交通运行速度。该方法具有计算效率高,分析结果准确等特点。
附图说明
图1.公交车载GPS数据估计道路交通运行速度流程图
图2.速度补偿模型示意图
图3.京开高速实施过程示意图
具体实施方式
以下结合附图所示,优选实施例作进一步详述,京开高速双向6车道,存在数条公交线路在主线上运营,实施流程如附图3。
步骤1:把京开高速分成上行和下行两个方向,见图3中的AB和BA两个方向,把上下行方向划分成多个独立的路段单元,见图中的link1和link2,计算路段单元的方位角,
步骤2:建立路段单元与公交车站对应关系表。
根据路段单元和公交车站地理空间位置,建立两者的归属关系,即公交车站属于哪一个路段单元,在link1和link2中各有一个公交车站,分别是BID1和BID2,BID1属于路段单元link1,BID2属于路段单元link2,其坐标分别记录为和数据见表1。
表1
步骤3:过滤路段单元周边的GPS点。
选择路段单元周边dlink范围内的GPS坐标点,这里取50米,过滤link1周边的全部GPS点,把过滤的路段单元周边GPS数据,按公交车辆编号BID1、时间t排序,构成路段单元link1周边GPS数据集合见表2所示。
表2
步骤4:消除公交停站对速度估计的影响,把公交车站周边dbus范围内的GPS点从集合中剔除,dbus取值50米,和均被剔除,构成新的集合见表3所示。
表3
步骤5:以BID1中时间最小的点为起点时间最大的点为终点计算所构成线段的方位角,比较与的绝对值大小,取Athreshold=30,在阈值范围内,因此BID1的GPS数据是有效数据,予以保留。构成集合数据与表3相同。
步骤6:在集合中,不存在车辆编号BID1仅有一条记录的GPS数车辆编号BID1的最大时间是12∶15∶50,最小时间是12∶15∶10,对应的GPS点是和计算两点的时间差值和距离差值,进而计算车辆的运行速度,见表4所示,公交车辆平均运行速度(千米/小时)。
表4
步骤8:本实施案例选用的道路属于高速公路,需利用速度补偿模型进行修正,才能得到道路交通运行速度,取速度修正模型参数λ=1/9,得到高速公路交通运行速度
因此:(千米/小时),根据速度区间分级,该速度下的道路处理畅通状态。
最后应说明的是:以上示例仅用以说明本发明而并非限制本发明所描述的技术方案;因此,尽管本说明书参照上述的示例对本发明已进行了详细的说明,但是,本领域的普通技术人员应当理解,仍然可以对本发明进行修改或等同替换;而一切不脱离发明的精神和范围的技术方案及其改进,其均应涵盖在本发明的权利要求范围当中。
Claims (1)
1.一种利用公交车辆GPS数据动态评估道路交通状态的方法,其特征在于包括以下步骤:
步骤1:划分路段单元,把道路分成上行和下行两个方向,划分路段单元是指为计算道路交通运行速度,把道路的上下行方向划分成多个独立的路段单元,相同的路段单元,认为交通流速度一致,根据交通流的方向计算路段单元的方位角,以正北方向为0度角,路段单元的方位角为路段单元起点和终点的连线与正北方向的顺时针夹角,路段单元i的方位角记为
步骤2:建立路段单元与公交车站对应关系表,根据路段单元和公交车站地理空间位置,建立两者的归属关系,即公交车站属于哪一个路段单元,公交车站j归属于路段单元i,其坐标信息记为
步骤3:过滤路段单元周边的GPS点,过滤路段单元周边dlink范围内的GPS坐标点,dlink取值依据道路宽度,取10-80米,若GPS坐标点n属于路段单元i,则记为包含GPS坐标和时间信息,记为
步骤4:把过滤的路段单元周边GPS数据,按公交车辆编号BIDk、时间t排序,构成路段单元i周边GPS数据集合
步骤5:为消除公交停靠站对速度估计的影响,把公交车站周边dbus范围内的GPS点从集合中剔除,dbus取值依据公交站台的长度,取10-50米,构成集合
步骤6:以BIDk中时间最小的点为起点时间最大的点为终点计算所构成线段的方位角,记为比较与的绝对值大小,若时,认为BIDk的GPS数据是有效数据,否则认为是无效数据,无效数据应剔除,Athreshold的取值为0-45,取值越小,结果精度越高,构成集合
步骤7:在集合中,剔除车辆编号BIDk仅有一条记录的GPS数据,构成集合
步骤8:计算路段单元i公交车辆平均运行速度集合中,取相同车辆编号最大时间和最小时间的GPS点,计算两点的时间差值和距离差值,进而计算车辆的运行速度,多辆公交车的平均运行速度即是该路段单元公交的平均公交运行速度
步骤9:利用速度补偿模型进行速度补偿,得到道路平均运行速度;道路交通运行速度与公交平均运行速度存在线性关系,速度补偿模型为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>V</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<msub>
<mi>link</mi>
<mi>i</mi>
</msub>
</mrow>
</msub>
<mo>=</mo>
<mi>&lambda;</mi>
<mo>&CenterDot;</mo>
<msubsup>
<mover>
<mi>V</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<msub>
<mi>link</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mi>b</mi>
<mi>u</mi>
<mi>s</mi>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mover>
<mi>V</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<msub>
<mi>link</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mi>b</mi>
<mi>u</mi>
<mi>s</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo><</mo>
<mi>&lambda;</mi>
<mo><</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
步骤10:根据速度的区间分级,把道路平均运行速度落在对应的区间,该区间表示的状态即是道路运行状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410781208.1A CN104851286B (zh) | 2014-12-18 | 2014-12-18 | 一种利用公交车辆gps数据动态评估道路交通状态的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410781208.1A CN104851286B (zh) | 2014-12-18 | 2014-12-18 | 一种利用公交车辆gps数据动态评估道路交通状态的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104851286A CN104851286A (zh) | 2015-08-19 |
CN104851286B true CN104851286B (zh) | 2017-11-07 |
Family
ID=53850901
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410781208.1A Expired - Fee Related CN104851286B (zh) | 2014-12-18 | 2014-12-18 | 一种利用公交车辆gps数据动态评估道路交通状态的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104851286B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105741545B (zh) * | 2016-03-16 | 2018-07-20 | 山东大学 | 一种基于公交车gnss时空轨迹数据判别交通状态的装置及方法 |
CN107798865B (zh) * | 2016-09-07 | 2020-11-03 | 阿里巴巴(中国)有限公司 | 一种公交线路行驶时间预估方法及装置 |
CN108010357B (zh) * | 2016-11-01 | 2020-11-27 | 武汉四维图新科技有限公司 | 限速信息校验/统计方法、装置及系统 |
CN106530791A (zh) * | 2016-11-03 | 2017-03-22 | 中兴软创科技股份有限公司 | 公交车到站位置匹配方法与系统 |
CN106781506B (zh) * | 2017-02-21 | 2019-02-05 | 济南全通信息科技有限公司 | 基于公交车gps数据的大范围城市公交线网实时运行水平评价方法 |
CN106772355B (zh) * | 2017-02-23 | 2019-10-18 | 四川大学 | 一种计算航迹关联门中心位置的补偿方法 |
CN107590998B (zh) * | 2017-08-16 | 2020-12-29 | 重庆市市政设计研究院 | 一种基于浮动车数据的道路节点运行状态识别系统 |
CN108111173B (zh) * | 2017-12-27 | 2021-03-26 | 东软集团股份有限公司 | 轨迹压缩方法、装置、存储介质及电子设备 |
CN108922178B (zh) * | 2018-07-01 | 2020-05-01 | 北京工业大学 | 基于公交多源数据的公交车辆实时满载率计算方法 |
CN109859515B (zh) * | 2019-01-28 | 2020-05-12 | 武汉元光科技有限公司 | 公交系统中gps定位数据补偿方法及电子设备 |
CN110457413B (zh) * | 2019-07-17 | 2022-05-10 | 中国第一汽车股份有限公司 | 行驶方向的确定方法、装置、设备及存储介质 |
CN110706477B (zh) * | 2019-10-12 | 2021-07-20 | 杭州律橙电子科技有限公司 | 一种公交站点匹配方法 |
CN111308432B (zh) * | 2019-12-03 | 2022-03-22 | 中国人民解放军63921部队 | 一种利用测速数据评估航天器测距数据精度的方法 |
CN111192460A (zh) * | 2019-12-31 | 2020-05-22 | 四川福立盟智慧环保大数据有限公司 | 一种用ai通过gps对渣土车运输行为识别方法 |
CN113611115B (zh) * | 2021-08-06 | 2022-06-21 | 安徽师范大学 | 一种基于路网敏感特征的车辆轨迹聚类方法 |
CN114372663A (zh) * | 2021-11-29 | 2022-04-19 | 北京华录高诚科技有限公司 | 一种公交线网运送效率指数计算方法、装置及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101650876A (zh) * | 2009-08-26 | 2010-02-17 | 重庆大学 | 城市路段交通流平均速度获取方法 |
CN101794507A (zh) * | 2009-07-13 | 2010-08-04 | 北京工业大学 | 基于浮动车数据的宏观路网交通状态评价方法 |
CN103106788A (zh) * | 2013-01-09 | 2013-05-15 | 众联智领科技(北京)有限公司 | 路况采集和服务系统及方法 |
CN103177561A (zh) * | 2011-12-26 | 2013-06-26 | 北京掌城科技有限公司 | 公交实时路况的生成方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105144259B (zh) * | 2013-01-24 | 2018-01-12 | 罗杰·安德烈·埃勒特森 | 交通监察和指导系统 |
-
2014
- 2014-12-18 CN CN201410781208.1A patent/CN104851286B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101794507A (zh) * | 2009-07-13 | 2010-08-04 | 北京工业大学 | 基于浮动车数据的宏观路网交通状态评价方法 |
CN101650876A (zh) * | 2009-08-26 | 2010-02-17 | 重庆大学 | 城市路段交通流平均速度获取方法 |
CN103177561A (zh) * | 2011-12-26 | 2013-06-26 | 北京掌城科技有限公司 | 公交实时路况的生成方法及系统 |
CN103106788A (zh) * | 2013-01-09 | 2013-05-15 | 众联智领科技(北京)有限公司 | 路况采集和服务系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104851286A (zh) | 2015-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104851286B (zh) | 一种利用公交车辆gps数据动态评估道路交通状态的方法 | |
CN106781479B (zh) | 一种基于手机信令数据实时获取高速公路运行状态的方法 | |
CN105865472B (zh) | 一种基于最佳油耗的车载导航方法 | |
CN103895649B (zh) | 一种驾驶员安全驾驶预警方法 | |
CN111862606B (zh) | 一种基于多源数据的非法营运车辆识别方法 | |
CN103794061B (zh) | 基于多种定位数据计算道路融合行程车速的方法 | |
CN106197458A (zh) | 一种基于手机信令数据和导航路线数据的手机用户出行方式识别方法 | |
CN110298500A (zh) | 一种基于出租车数据和城市路网的城市交通轨迹数据集生成方法 | |
CN103632540B (zh) | 基于浮动车数据的城市主干道交通运行信息处理方法 | |
CN105825669A (zh) | 一种识别城市快速路交通瓶颈的系统和方法 | |
CN104916154B (zh) | 一种兼容北斗cors公交精准定位系统及其工作方法 | |
CN105608505A (zh) | 一种基于手机信令数据的居民轨道交通出行方式识别方法 | |
CN104504099A (zh) | 基于位置轨迹的交通出行状态切分方法 | |
CN103177585B (zh) | 基于浮动车数据计算路段各转向平均行程车速方法 | |
CN106846806A (zh) | 基于Isolation Forest的城市道路交通异常检测方法 | |
CN107203945A (zh) | 车险分级评估方法及装置 | |
CN105608895B (zh) | 一种基于局部异常因子的城市交通拥堵路段检测方法 | |
CN105096590B (zh) | 交通信息生成方法和交通信息生成设备 | |
CN103093625B (zh) | 一种基于可信度验证的城市道路交通状态实时估计方法 | |
CN103942952B (zh) | 一种路网功能层次状态等级评估方法 | |
CN106227859A (zh) | 从gps数据中识别交通工具的方法 | |
CN107240264A (zh) | 一种车辆非有效行驶轨迹识别方法和城市道路设施规划方法 | |
WO2023123616A1 (zh) | 一种快速路车辆od位置提取方法及系统 | |
CN105913666A (zh) | 一种快速道路主线可变限速标志优化布设方法 | |
CN105547310A (zh) | 一种基于pm2.5健康出行的路径规划装置及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171107 Termination date: 20181218 |
|
CF01 | Termination of patent right due to non-payment of annual fee |