CN111998828B - 一种基于便携式gps的道路坡度估算方法 - Google Patents

一种基于便携式gps的道路坡度估算方法 Download PDF

Info

Publication number
CN111998828B
CN111998828B CN202011006738.0A CN202011006738A CN111998828B CN 111998828 B CN111998828 B CN 111998828B CN 202011006738 A CN202011006738 A CN 202011006738A CN 111998828 B CN111998828 B CN 111998828B
Authority
CN
China
Prior art keywords
road
data
gps
gradient
section
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
Application number
CN202011006738.0A
Other languages
English (en)
Other versions
CN111998828A (zh
Inventor
余飞
郑君瑜
钟庄敏
黄志炯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jinan University
Original Assignee
Jinan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jinan University filed Critical Jinan University
Priority to CN202011006738.0A priority Critical patent/CN111998828B/zh
Publication of CN111998828A publication Critical patent/CN111998828A/zh
Application granted granted Critical
Publication of CN111998828B publication Critical patent/CN111998828B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C9/00Measuring inclination, e.g. by clinometers, by levels

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Navigation (AREA)

Abstract

本发明公开了一种基于便携式GPS的道路坡度估算方法,方法包括以下步骤:将多个GPS接收器放置在车辆上收集数据;对收集得到的数据进行预处理;测试道路划分,考虑数据点的数量,确定路段分割长度并依据确定的长度将测试道路划分为多个分路段;分路段高程数据对齐,对获取的高程数据进行对齐处理;计算路段的道路坡度并定量精度,对每个分路段进行分路段的坡度计算,使用每个分路段对齐后的高程数据进行线性回归来估算道路坡度并根据线性回归的斜率估计的标准误差来推断测量的道路坡度的精度。本发明方法可以方便精确地通过GPS数据解析获取,不受所测试道路类型及车辆行驶状态的影响,通过对累积距离和高程的校正能获得更加精确的坡度结果。

Description

一种基于便携式GPS的道路坡度估算方法
技术领域
本发明属于交通环保技术领域,具体涉及一种基于便携式GPS的道路坡度估算方法。
背景技术
高精度地图提供的信息有利于汽车在行驶中的***控制,有利于减少交通事故的发生,而且自动驾驶技术的成熟发展也离不开高精度地图数据的支持,而道路坡度又是高精度地图的重要组成元素,所以准确高效的测量路面坡度具有重要意义。
目前国内外在路面坡度测量上,如想获得高精度的路面坡度数据需要价格非常昂贵的专业测绘仪器且涉及数据融合和解算等,较为复杂。对于非测绘专业人员在如想获得道路坡度数据则非常困难,仅用简单的GPS接收器或水平仪等设备很难获得高精度的道路坡度数据,且不能满足对路面坡度的动态快速测量需求。
专利CN101838958A提出了一种道路坡度的检测方法,该方法是基于测绘部门的等高线图和道路规划图进行处理计算,主要针对道路设计规划中的坡度获取,需要首先获得从测绘部门的等高线图和道路规划图。该方法所获取的坡度为静态数据,且在道路规划与实际建设后的路面坡度往往存在些偏差,不能准确反映机动车实际行驶中所处道路的坡度信息。而且对于普通使用者也难以从测绘部门获取相关道路规划图再进行专业分析计算。
专利CN102313535A中提出了一种坡度的检测方法,该方法通过加速度传感器和速度信号模块根据经过校正的校正值、加速度信号、机动车加速度值和地球重力加速获得坡度值。该发明依赖于安装在车架上传感器信号的推算,用到不同的车辆时需要分别进行安装调试。该方法更多为了反映车辆的起伏状态变化,无法获取真实的道路坡度值且容易受传感器安装及传感器的精度影响。该方法所获取坡度难以满足机动车尾气排放测量及估算时的高精度和便携测量需求。
发明内容
本发明的主要目的在于克服现有技术的缺点与不足,提出一种基于便携式GPS的道路坡度估算方法,该方法所用设备价格低廉、操作方便、获取的坡度数据精度高、方便普及与应用,在高精度地图坡度数据完善及机动车动态排放估算等领域均有重要的实用价值。
为了达到上述目的,本发明采用以下技术方案:
一种基于便携式GPS的道路坡度估算方法,其特征在于,包括以下步骤:
利用一个或多个便携式GPS接收器收集数据,每个GPS接收器采集的数据包括道路的经纬度和高程数据;
对收集得到的数据进行预处理,根据每个GPS收集到的经纬度数据,计算每条路线轨迹的累积距离并参考实际道路长度进行校正;
对测试道路进行路段划分,考虑路段中数据点的数量,确定路段分割长度并依据确定的长度将测试道路划分为多个分路段,所述数据点为测试中GPS接收器记录的数据;
分路段高程数据对齐,将每个分路段涉及的多组GPS数据合并为数据集,对各组高程数据进行对齐处理,解决平均高程数据误差问题;
计算测试道路的道路坡度并定量精度,对每个分路段进行分路段的坡度计算,使用每个分路段对齐后的高程数据进行线性回归来估算道路坡度并根据线性回归的斜率估计的标准误差来推断测量道路坡度的精度。
进一步的,所述计算路段的累积距离具体为:
假设GPS接收器数量为m,车辆重复n次行驶经过测试道路,则根据GPS接收器的经纬度数据,得到m*n条路线轨迹,基于每条GPS轨迹逐秒速度计算路线轨迹的累积行驶距离,对m*n条路线轨迹累积距离求平均即为路线轨迹累积距离,公式如下:
Figure BDA0002696169180000031
其中,d为路线轨迹累积距离,i表示第i个GPS,1≤i≤n,j表示第j次测试,1≤j≤n,li,j为第i个GPS第j次测试的累积距离。
进一步的,所述对参考实际道路长度进行校正具体为通过对每条路线轨迹累积距离***一致的起点和终点,将其校正为真实距离,每条路线轨迹累积距离的校正因子如下表示:
θi,j=d/li,j
其中,θi,j为第i个GPS第j次测试的累积距离校正因子。
进一步的,所述确定路段分割长度具体为:
数据点数量影响测试道路坡度的精度,所以路段分割长度需保证分割后每个分路段数据点的数量,而数据点的数量取决于数据记录频率、车速以及分路段长度;
分路段的长度越长其数据点数越多获得的坡度精度越高,但是分路段过长会导致路段内的实际变化被平均,从而导致对真实坡度变化的低估;如果分路段长度太短,则分路段内的数据点过少会导致估算坡度不精确。
进一步的,所述数据点数量的计算公式如下:
Figure BDA0002696169180000041
其中,number为分路段中数据点数量,Δli,j为第i个GPS第j次测试的路段长度;Δvi,j为第i个GPS第j次测试的路段平均速度,H为GPS数据记录频率。
进一步的,每个分割后分路段数据点的数量大于20。
进一步的,所述分路段高程数据进行对齐处理具体为:
对每个分路段获取的多组GPS数据进行融合,合并为数据集,对测试道路上所有的高程数据计算平均高程作为参考绝对高程,然后使用每个数据点的高程减去参考绝对高程,每个数据点之间的相对高程差不变。
进一步的,所述计算路段的道路坡度具体为:
对分路段进行路段坡度计算,使用每个分路段对齐后的高程数据点进行一元线性回归拟合来估算道路坡度,拟合直线的斜率乘以100即为测试道路坡度百分比,将全部分路段的坡度进行连接即为测试道路的坡度信息。
进一步的,所述道路坡度的精度具体为拟合直线斜率的标准误差。
本发明与现有技术相比,具有如下优点和有益效果:
1、本发明方法通过对测试道路进行分段并计算坡度,对不同路段的高程数据进行对齐处理,提高了坡度计算准确度;本发明利用低精度的便携式GPS实现了高精度道路坡度的实时快速测量,不受所测试道路类型及车辆行驶状态的影响,解决了普通便携式GPS无法获取高精度坡度数据的问题。
2、本发明方法成本较低、操作方法简单、不依赖于专业测绘工具即可实现高精度道路坡度的测量。
3、本发明方法给出了道路坡度准确度的定量思路,使用者可以根据对坡度精度的需求选择不同的GPS个数进行坡度测量。
附图说明
图1是本发明方法的流程图;
图2a是路段高程数据对齐前数据;
图2b是路段高程数据对齐后数据;
图3a是路段高程数据对齐前拟合坡度结果与精度;
图3b是路段高程数据对齐后拟合坡度结果与精度;
图4是所测试路线的坡度结果。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
实施例
以机动车实际道路排放测试需要知道所测试路线的坡度数据且要求有较高的精度为例,如图1所示,本发明一种基于便携式GPS的道路坡度估算方法,包括以下步骤:
S1、利用多个GPS接收器收集数据,具体为:
将1个或多个便携式GPS接收器放置在车辆上,车辆重复多次正常行驶经过测试道路,采集GPS接收器的所有记录数据,每个GPS接收器采集的数据包括道路的经纬度和高程数据。
在本实施例中,采用4个便携式GPS放置在测试车辆的前端固定,沿机动车排放测试道路正常行驶,重复行驶经过路段4次,共采集了16条路线轨迹数据。
S2、计算路线轨迹的累积距离,对参考实际道路长度进行校正,具体为:
假设GPS接收器数量为m,车辆重复n次行驶经过测试道路,则根据GPS接收器的经纬度数据,得到m*n条路线轨迹,基于每条GPS轨迹逐秒速度计算路线轨迹的累积行驶距离,对m*n条路线轨迹累积距离求平均即为路线轨迹累积距离,路线轨迹累积距离计算公式如下:
Figure BDA0002696169180000061
其中,d为路线轨迹累积距离,i表示第i个GPS,1≤i≤n,j表示第j次测试,1≤j≤n,li,j为第i个GPS第j次测试的累积距离。
参考实际道路长度进行校正具体为:
通过对每条路线轨迹累积距离***一致的起点和终点,将其校正为真实距离,每条路线轨迹累积距离的校正因子如下表示:
θi,j=d/li,j
其中,θi,j为第i个GPS第j次测试的累积距离校正因子。
在本实施例中,对上述16条轨迹进行道路累积距离校正,16条轨迹的累积距离如下表所示,所获得的平均累积距离为35.3km标准偏差为±0.7%。
Figure BDA0002696169180000062
S3、测试道路划分,具体为:
路段内的数据点数量同样决定了所获得道路坡度的精度,为保证每个分路段内所保留的数据点数量,首先要确定划分路段长度。每个分路段内数据点的数量取决于数据记录频率、车速和分别路段长度,在所用GPS为1Hz的数据记录频率下,以速度100km h-1的典型高速公路速度为例,在0.1km的路段内单个GPS仅有3个数据点,但基于m个GPS的n次运行,在100km h-1的速度0.1km的路段长度上则有3*m*n个数据点;如果分路段的长度太长其数据点数越多获得的坡度精度越高,但分路段过长会导致分路段内的实际变化被平均,从而导致对真实坡度变化的低估;如果分路段长度太短,则分路段内的数据点过少可能会导致估算坡度不精确;一般每个分路段的数据点数需要大于20个才能满足坡度精度计算需求。
数据点数量的计算公式如下:
Figure BDA0002696169180000071
其中,number为分路段中数据点数量,Δli,j为第i个GPS第j次测试的路段长度;Δvi,j为第i个GPS第j次测试的路段平均速度,H为GPS数据记录频率。
在本实施例中,划分长度为0.1km,测试道路被划分为353段。因行驶速度越大每个分路段上的数据点越少,以最大行驶速度为120km/h计算则每个分路段上至少有48个数据点。
S4、分路段高程数据对齐,具体为:
首先对每个分路段内的多条数据合并为一个数据集,因不同GPS的绝对高程存在随机误差,需要对所获取的每个分路段的高程数据进行对齐处理。计算沿该路段每次行驶的平均海拔。对每个路段融合之后的数据计算平均高程作为参考绝对高程,然后每个数据点的高程减去绝对高程参考值,其每个数据点之间的相对高程差异不变。
在另一个实施例中,以从GPS B和GPS C获得的高程数据为例,如图2a、图2b所示。
S5、计算测试道路的道路坡度,具体为:
对分路段进行每个分路段坡度的计算,使用每个分路段对齐后的高程数据点进行一元线性回归拟合来估算待测量坡度。拟合直线的斜率乘以100以百分比来表示道路坡度,全部分路段的坡度进行连接即为测试道路的坡度信息;拟合直线斜率的标准偏差为所获取道路坡度的精度。一个分路段中的数据点数量可能会因车速的不同而有所不同,当选择适当的分路段长时该方法可用于计算任何位置的坡度。
在本实施例中,通过对每条分路段上不同轨迹数据点进行校正,然后进行拟合得到该测试道路所有分路段的坡度及其标准偏差,结果如图4所示,获得该测试道路的坡度范围从-8%-6%,精度在±0.5%以内。
在另一个实施例中,如图3a、图3b所示,经过高程数据对齐所获取的坡度精度由±4.6%提高到±0.75%。
另外选择用一台高分辨率厘米级差分GPS的商业坡度仪获取道路坡度数据与本实施施所获得的结果进行对比。对比结果显示通过本发明所提供方法获取的坡度结果与高分辨率厘米级GPS获取坡度结果相近,证实了本发明方法的可靠性。
本发明的坡度测量方法可以方便精确地通过GPS数据解析获取,不受所测试道路类型及车辆行驶状态的影响,通过对累积距离和高程的校正能获得更加精确的坡度结果。
还需要说明的是,在本说明书中,诸如术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其他实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (4)

1.一种基于便携式GPS的道路坡度估算方法,其特征在于,包括以下步骤:
利用一个或多个便携式GPS接收器收集数据,每个GPS接收器采集的数据包括道路的经纬度和高程数据;
对收集得到的数据进行预处理,根据每个GPS收集到的经纬度数据,计算每条路线轨迹的累积距离并参考实际道路长度进行校正;
计算路线轨迹的累积距离具体为:
假设GPS接收器数量为m,车辆重复n次行驶经过测试道路,则根据GPS接收器的经纬度数据,得到m*n条路线轨迹,基于每条GPS轨迹逐秒速度计算路线轨迹的累积行驶距离,对m*n条路线轨迹累积距离求平均即为路线轨迹累积距离,公式如下:
Figure FDA0003533959830000011
其中,d为路线轨迹累积距离,i表示第i个GPS,1≤i≤m,j表示第j次测试,1≤j≤n,li,j为第i个GPS第j次测试的累积距离;
参考实际道路长度进行校正具体为通过对每条路线轨迹累积距离***一致的起点和终点,将其校正为真实距离,每条路线轨迹累积距离的校正因子如下表示:
θi,j=d/li,j
其中,θi,j为第i个GPS第j次测试的累积距离校正因子;
对测试道路进行路段划分,考虑路段中数据点的数量,确定路段分割长度并依据确定的长度将测试道路划分为多个分路段,所述数据点为测试中GPS接收器记录的数据;
确定路段分割长度具体为:
数据点数量影响测试道路坡度的精度,所以路段分割长度需保证分割后每个分路段数据点的数量,而数据点的数量取决于数据记录频率、车速以及分路段长度;
分路段的长度越长其数据点数越多获得的坡度精度越高,但是分路段过长会导致路段内的实际变化被平均,从而导致对真实坡度变化的低估;如果分路段长度太短,则分路段内的数据点过少会导致估算坡度不精确;
分路段高程数据对齐,将每个分路段涉及的多组GPS数据合并为数据集,对各组高程数据进行对齐处理,解决平均高程数据误差问题;具体为:
对每个分路段获取的多组GPS数据进行融合,合并为数据集,对测试道路上所有的高程数据计算平均高程作为参考绝对高程,然后使用每个数据点的高程减去参考绝对高程,每个数据点之间的相对高程差不变;
计算测试道路的道路坡度并定量精度,对每个分路段进行分路段的坡度计算,使用每个分路段对齐后的高程数据进行线性回归来估算道路坡度并根据线性回归的斜率估计的标准误差来推断测量道路坡度的精度;具体为:
对分路段进行路段坡度计算,使用每个分路段对齐后的高程数据点进行一元线性回归拟合来估算道路坡度,拟合直线的斜率乘以100即为测试道路坡度百分比,将全部分路段的坡度进行连接即为测试道路的坡度信息。
2.根据权利要求1所述的一种基于便携式GPS的道路坡度估算方法,其特征在于,所述数据点数量的计算公式如下:
Figure FDA0003533959830000021
其中,number为分路段中数据点数量,Δli,j为第i个GPS第j次测试的路段长度;Δvi,j为第i个GPS第j次测试的路段平均速度,H为GPS数据记录频率。
3.根据权利要求1所述的一种基于便携式GPS的道路坡度估算方法,其特征在于,每个分割后分路段数据点的数量大于20。
4.根据权利要求1所述的一种基于便携式GPS的道路坡度估算方法,其特征在于,所述道路坡度的精度为 拟合直线斜率的标准误差。
CN202011006738.0A 2020-09-23 2020-09-23 一种基于便携式gps的道路坡度估算方法 Active CN111998828B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011006738.0A CN111998828B (zh) 2020-09-23 2020-09-23 一种基于便携式gps的道路坡度估算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011006738.0A CN111998828B (zh) 2020-09-23 2020-09-23 一种基于便携式gps的道路坡度估算方法

Publications (2)

Publication Number Publication Date
CN111998828A CN111998828A (zh) 2020-11-27
CN111998828B true CN111998828B (zh) 2022-05-10

Family

ID=73474732

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011006738.0A Active CN111998828B (zh) 2020-09-23 2020-09-23 一种基于便携式gps的道路坡度估算方法

Country Status (1)

Country Link
CN (1) CN111998828B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112557059B (zh) * 2020-12-15 2023-01-06 潍柴动力股份有限公司 一种坡度采集***、方法和车辆
CN113091708A (zh) * 2021-04-06 2021-07-09 中山大学 基于多样本gps质控数据的道路分段坡度估计方法及***
CN113295826B (zh) * 2021-05-24 2023-06-27 暨南大学 一种基于车联网的机动车尾气排放测试管理***及智能诊断方法
CN113947895A (zh) * 2021-09-14 2022-01-18 中寰卫星导航通信有限公司 一种数据处理方法、装置、设备以及存储介质
CN115331430B (zh) * 2022-07-14 2024-03-08 中通客车股份有限公司 一种智能网联车路谱采集分析方法及***

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006275969A (ja) * 2005-03-30 2006-10-12 Clarion Co Ltd 傾斜角算出装置及びナビゲーション装置
CN110379150B (zh) * 2018-09-13 2021-05-25 北京京东尚科信息技术有限公司 一种评估道路通行性的方法和装置
CN110161513B (zh) * 2018-09-28 2023-11-14 腾讯科技(北京)有限公司 估计道路坡度的方法、装置、存储介质和计算机设备
CN110186428B (zh) * 2019-05-27 2022-03-18 北京泓达九通科技发展有限公司 一种基于大数据曲线拟合的道路坡度确定方法及***

Also Published As

Publication number Publication date
CN111998828A (zh) 2020-11-27

Similar Documents

Publication Publication Date Title
CN111998828B (zh) 一种基于便携式gps的道路坡度估算方法
CN1590965B (zh) 用于在导航***中检测车辆位置的装置和方法
CN107402006B (zh) 基于轨道几何特征信息匹配的列车精密定位方法和***
CN111767354B (zh) 一种高精度地图精度评估方法
CN103162689B (zh) 辅助车载定位***及车辆的辅助定位方法
Boroujeni et al. Road grade measurement using in-vehicle, stand-alone GPS with barometric altimeter
CN100578153C (zh) 车载速度测量仪的校准方法
CN101644756B (zh) 直线型gps动态定位精度检定仪及其检定方法
CN101539589B (zh) 一种机动车gps测速仪检测方法
CN113108808B (zh) 一种车辆里程计在线检定***及检定方法
CN104898139A (zh) 一种车辆定位纠偏方法及装置
CN102175463A (zh) 一种基于改进卡尔曼滤波的汽车路试制动性能检测方法
CN102004246B (zh) 列车车载雷达速度传感器天线夹角偏差故障诊断及速度读数校正方法
CN105882684A (zh) 一种城市轨道交通公里标定标方法
CN109606430B (zh) 轨道车自动定位方法、自动定位***及自动驾驶装置
CN114719884A (zh) 一种惯导***姿态测量精度评估方法及应用
CN107907134A (zh) 一种里程信息辅助地磁匹配的车辆定位***与方法
CN102980592A (zh) 一种通过gps经纬度自动计算车辆脉冲系数的方法及装置
Faghri et al. Application of GPS in traffic management systems
CN109693686A (zh) 一种基于双天线gnss/ins的轨向不平顺测量方法
CN114136275A (zh) 一种轨道线路状态检测装置及路基沉降检测方法
CN110820495A (zh) 非机动车道路面平整度检测方法、***及计算机存储介质
CN108413945B (zh) 轨道坐标点经纬高测量装置及方法
CN110629647B (zh) 一种基于时间触发的任意车速路面高程数据采集的方法
CN113091708A (zh) 基于多样本gps质控数据的道路分段坡度估计方法及***

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