CN104900061A - 路段行程时间监测方法及装置 - Google Patents
路段行程时间监测方法及装置 Download PDFInfo
- Publication number
- CN104900061A CN104900061A CN201510286096.7A CN201510286096A CN104900061A CN 104900061 A CN104900061 A CN 104900061A CN 201510286096 A CN201510286096 A CN 201510286096A CN 104900061 A CN104900061 A CN 104900061A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- road section
- prediction model
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000012544 monitoring process Methods 0.000 title claims abstract description 30
- 238000005070 sampling Methods 0.000 claims description 22
- 238000012545 processing Methods 0.000 claims description 7
- 238000012806 monitoring device Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 2
- 238000012935 Averaging Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- BVPWJMCABCPUQY-UHFFFAOYSA-N 4-amino-5-chloro-2-methoxy-N-[1-(phenylmethyl)-4-piperidinyl]benzamide Chemical compound COC1=CC(N)=C(Cl)C=C1C(=O)NC1CCN(CC=2C=CC=CC=2)CC1 BVPWJMCABCPUQY-UHFFFAOYSA-N 0.000 description 1
- 206010039203 Road traffic accident Diseases 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
Landscapes
- Traffic Control Systems (AREA)
Abstract
本发明公开了一种路段行程时间监测方法及装置。该方法包括:在一时段内,对经地图匹配过程匹配到监测路段上的浮动车的位置数据进行采样;根据采样的位置数据,确定监测路段在所述时段内的实际行程时间;将所述时段内的实际行程时间分别输入至第一预测模型和第二预测模型,得出监测路段在下一时段的第一预测行程时间和第二预测行程时间;基于第一预测模型的历史预测误差和第二预测模型的历史预测误差确定第一预测模型和第二预测模型的权重系数;基于第一预测模型和第二预测模型的权重系数,对第一预测行程时间和第二预测行程时间进行加权平均,得出监测路段在所述下一时段内的最终预测行程时间。由此,可提高路段行程时间预测的准确性和可靠性。
Description
技术领域
本发明涉及城市交通领域,具体地,涉及一种路段行程时间监测方法及装置。
背景技术
随着城市人口的增长和现代化的不断推进,交通问题已愈加严重。日常的交通拥挤和堵塞,极其消极地影响着人们的正常生活与工作,而交通事故的发生甚至严重危及到人们的生命。实时有效的交通信息是城市交通管控部门了解交通状况,为出行者提供有效信息、合理进行交通诱导以缓解交通压力,解决交通问题的关键。路段行程时间就是反映交通情况的重要参数之一,也是交通诱导系统中需预测的一项重要指标。
路段行程时间预测是基于当前时段的路段行程时间来预测下一时段的路段行程时间,以为出行者提供交通诱导。通常采用单一预测模型来预测路段行程时间。由于仅依赖于该单一预测模型,因此,预测结果的准确性和可靠性往往不高,无法为出行者提供更为准确、可靠的交通诱导信息。
发明内容
本发明的目的是提供一种路段行程时间监测方法及装置,以提高路段行程时间预测的准确性和可靠性。
为了实现上述目的,本发明提供一种路段行程时间监测方法,该方法包括:在一时段内,对经地图匹配过程匹配到监测路段上的浮动车的位置数据进行采样;根据所采样的位置数据,确定所述监测路段在所述时段内的实际行程时间;将所述时段内的实际行程时间分别输入至第一预测模型和第二预测模型,得出所述监测路段在下一时段内的第一预测行程时间和第二预测行程时间;基于所述第一预测模型的历史预测误差和所述第二预测模型的历史预测误差来确定所述第一预测模型和所述第二预测模型的权重系数;以及基于所述第一预测模型和所述第二预测模型的权重系数,对所述第一预测行程时间和所述第二预测行程时间进行加权平均,得出所述监测路段在所述下一时段内的最终预测行程时间。
本发明还提供一种路段行程时间监测装置,该装置包括:数据采样单元,用于在一时段内,对经地图匹配过程匹配到监测路段上的浮动车的位置数据进行采样;第一处理单元,用于根据所采样的位置数据,确定所述监测路段在所述时段内的实际行程时间;第一预测单元,用于将所述时段内的实际行程时间输入至第一预测模型,得出所述监测路段在下一时段内的第一预测行程时间;第二预测单元,用于将所述时段内的实际行程时间输入至第二预测模型,得出所述监测路段在下一时段内的第二预测行程时间;权重确定单元,用于基于所述第一预测模型的历史预测误差和所述第二预测模型的历史预测误差来确定所述第一预测模型和所述第二预测模型的权重系数;以及第二处理单元,用于基于所述第一预测模型和所述第二预测模型的权重系数,对所述第一预测行程时间和所述第二预测行程时间进行加权平均,得出所述监测路段在所述下一时段内的最终预测行程时间。
在上述技术方案中,通过采用两种预测模型来分别预测监测路段在下一时段的路段行程时间,并对两种预测模型得到的预测结果进行融合,从而得出最终的预测结果。由此,可以提高预测结果的准确性和可靠性。此外,可以根据第一预测模型和第二预测模型的历史预测误差动态调整第一预测模型和第二预测模型的权重系数,以保证精度相对高的预测模型得出的结果能够占有更大的权重,从而可以进一步提高预测结果的准确性(例如,相较于单一预测模型,准确性可以大约提高(预测误差可减少15%~20%,平稳度可提高15%~20%)。通过本发明提供的路段行程时间监测方法及装置,能够为出行者及交通监管部门提供更为准确的路段行程时间的预测结果,从而便于出行者选择最佳行驶路段、以及便于交通监管部门及时采取相应的交通诱导措施。
本发明的其他特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本发明的进一步理解,并且构成说明书的一部分,与下面的具体实施方式一起用于解释本发明,但并不构成对本发明的限制。在附图中:
图1是根据本发明的实施方式的路段行程时间监测方法的流程图;
图2是在一种示例情况下,在一时段内,监测路段的示意图;
图3是根据本发明的实施方式的用于对浮动车进行地图匹配的过程的流程图;
图4是在执行本发明提供的浮动车地图匹配过程时,确定的椭圆区域及网格匹配的示意图;
图5是根据本发明的实施方式的路段行程时间监测装置的框图;以及
图6是基于本发明提供的路段行程时间监测方法及装置得到的监测路段的预测行程时间与实际行程时间、和采用单一预测模型得到的预测行程时间之间的结果对比图。
具体实施方式
以下结合附图对本发明的具体实施方式进行详细说明。应当理解的是,此处所描述的具体实施方式仅用于说明和解释本发明,并不用于限制本发明。
图1示出了根据本发明的实施方式的路段行程时间监测方法的流程图。如图1所示,该方法可以包括:步骤S1,在一时段内,对经地图匹配过程匹配到监测路段上的浮动车的位置数据进行采样;步骤S2,根据所采样的位置数据,确定所述监测路段在所述时段内的实际行程时间;步骤S3,将所述时段内的实际行程时间分别输入至第一预测模型和第二预测模型,得出所述监测路段在下一时段内的第一预测行程时间和第二预测行程时间;步骤S4,基于所述第一预测模型的历史预测误差和所述第二预测模型的历史预测误差来确定所述第一预测模型和所述第二预测模型的权重系数;以及步骤S5,基于所述第一预测模型和所述第二预测模型的权重系数,对所述第一预测行程时间和所述第二预测行程时间进行加权平均,得出所述监测路段在所述下一时段内的最终预测行程时间。
具体地,首先,在步骤S1中,在一时段内,对经地图匹配过程匹配到监测路段上的浮动车的位置数据进行采样。假设该时段(亦可称为当前时段)被标记为tc,则上一时段可以被标记为tc-1,下一时段可以被标记为tc+1,以此类推。在本发明中,每个时段的时间间隔可以是相同的。对监测路段的行程时间进行预测就是根据时段tc内的路段行程时间来预测下一时段tc+1的路段行程时间。
图2示出了在一种示例情况下,在时段tc内,监测路段的示意图。应当理解的是,该图2所示的监测路段示意图仅仅用于说明本发明和便于理解本发明,而不用于限制本发明。
如图2所示,在时段tc内,经地图匹配过程匹配到该监测路段N上的浮动车有两个,分别为浮动车A和浮动车B。在时段tc内,可以对浮动车A和浮动车B的位置数据进行采样。例如,可以以固定的采样间隔进行采样,从而获得浮动车A和浮动车B在该监测路段N上的位置移动信息。例如,如图2所示,在时段tc内,第一次采样得到浮动车A的位置数据L1,1和浮动车B的位置数据L1,2,第二次采样得到浮动车A的位置数据L2,1和浮动车B的位置数据L2,2,第三次采样得到浮动车B的位置数据L3,2。由于在第三次采样时浮动车A已驶出监测路段N,因此,浮动车A在该监测路段N上的位置数据记录为两个,分别为L1,1和L2,1,而浮动车B在该监测路段N上的位置数据记录为三个,分别为L1,2、L2,2和L3,2。
在采样得出浮动车的位置数据之后,进行步骤S2,可以根据所采样的位置数据,确定所述监测路段在所述时段tc内的实际行程时间。例如,可以通过以下等式(1)和等式(2)来确定所述监测路段在所述时段tc内的实际行程时间:
其中,Li,j可以表示在所述时段tc内采样的第j辆浮动车在该监测路段上的第i个位置数据;Li+1,j可以表示在所述时段tc内采样的第j辆浮动车在该监测路段上的第i+1个位置数据;Ti,j可以表示第j辆浮动车到达所述第i个位置数据所表示的位置时的时间;Ti+1,j可以表示第j辆浮动车到达所述第i+1个位置数据所表示的位置时的时间;nj可以表示在所述时段tc内采样的第j辆浮动车在该监测路段上的位置数据总个数,并且nj≥2;vj可以表示在所述时段tc内,第j辆浮动车在该监测路段上的平均车速;T可以表示所述监测路段在所述时段tc内的实际行程时间;L可以表示所述监测路段的总距,例如如图2所示;以及m可以表示在所述时段tc内,所述监测路段上的浮动车总数。
下面以图2为例进行说明。如图2所示,经采样得到浮动车A的位置数据为L1,1和L2,1,经采样得到的浮动车B的位置数据为L1,2、L2,2和L3,2。以浮动车A作为该监测路段N上的第一辆浮动车,以浮动车B作为该监测路段N上的第二辆浮动车,那么存在m=2,以及对于浮动车A而言,在所述时段tc内采样的位置数据总个数n=2,对于浮动车B而言,在所述时段tc内采样的位置数据总个数n=3。之后,可以按照等式(1)分别确定出浮动车A和浮动车B在该监测路段N上的平均车速v1和v2。例如:
其中,T2,1可以表示浮动车A到达位置数据L2,1所表示的位置时的时间,T1,1可以表示浮动车A到达位置数据L1,1所表示的位置时的时间,T3,2可以表示浮动车B到达位置数据L3,2所表示的位置时的时间,T2,2可以表示浮动车B到达位置数据L2,2所表示的位置时的时间,以及T1,2可以表示浮动车B到达位置数据L1,2所表示的位置时的时间。这些时间参数可以根据采样时间来确定。例如,T1,2=T1,1=第一次采样时间,T2,2=T2,1=第一次采样时间,T3,2=第三次采样时间。
在得出浮动车A和浮动车B的平均车速v1、v2之后,就可以根据这两个浮动车的平均车速v1、v2和监测路段N的总距L、使用上述等式(2)确定出在时段tc内,监测路段N的实际行程时间T,例如,该实际行程时间为:
在确定出监测路段在所述时段tc内的实际行程时间T之后,进行步骤S3,将所述时段tc内的实际行程时间T分别输入至第一预测模型和第二预测模型,得出所述监测路段在下一时段tc+1内的第一预测行程时间Tp1和第二预测行程时间Tp2。
在本发明中,第一预测模型与第二预测模型是不同的预测模型,并且两个预测模型可以分别选自现有的用于路段行程时间预测的预测模型。例如,所述第一预测模型可以例如为卡尔曼(Kalman)滤波模型,以及所述第二预测模型可以例如为时间序列(ARIMA)模型。应当理解的是,如何使用卡尔曼(Kalman)滤波模型和时间序列(ARIMA)模型来单独进行行程时间预测的原理与实现方法是本领域的技术人员公知的,因此,本发明在此不进行赘述。
在通过两种预测模型得到两个预测行程时间之后,可以对这两个预测行程时间进行加权平均,得出所述监测路段在所述下一时段tc+1内的最终预测行程时间Tp。例如,可以通过以下等式(3)来确定所述最终预测行程时间Tp:
Tp=ω1Tp1+ω2Tp2 (3)
其中,Tp可以表示所述最终预测行程时间;Tp1可以表示所述第一预测行程时间;Tp2可以表示所述第二预测行程时间;ω1可以表示所述第一预测模型的权重系数;以及ω2可以表示所述第二预测模型的权重系数。
为了保证两个预测模型中精度更高的预测模型所得的结果能够占有更大的权重,以使预测结果更为准确、可靠,在本发明中,可以基于所述第一预测模型的历史预测误差和所述第二预测模型的历史预测误差来确定所述第一预测模型和所述第二预测模型的权重系数。
具体地,可以例如通过以下等式(4)和等式(5)来确定所述第一预测模型和所述第二预测模型的权重系数:
其中,e1(k)可以表示所述第一预测模型的第k个历史预测误差;e2(k)可以表示所述第二预测模型的第k个历史预测误差;K1可以表示所述第一预测模型的历史预测误差的总个数;以及K2可以表示所述第二预测模型的历史预测误差的总个数。
在本发明中,预测误差是指对于某个监测路段,针对某一时段,检测出的该时段的实际行程时间与预测出的该时段的预测行程时间之间的相对误差的绝对值。而历史预测误差是指在预测监测路段在下一时段tc+1的路段行程时间之前,获得的针对该下一时段tc+1之前的每一时段(例如,当前时段tc、上一时段tc-1、再上一时段tc-2等等)的预测误差。并且,随着监测过程的进行,每种预测模型的历史预测误差不断递增。
例如,以时段tc为例,根据本发明提供的路段行程时间监测方法,在上一时段tc-1内利用第一预测模型可以得出针对时段tc的第一预测行程时间,以及利用第二预测模型可以得出针对时段tc的第二预测行程时间。而在时段tc内,可以根据在该时段tc内采样的位置数据来确定出该时段tc内的实际行程时间。此时,针对时段tc,就可以根据确定出的实际行程时间和第一预测行程时间以及第二预测行程时间,分别得出第一预测模型的预测误差和第二预测模型的预测误差,即,第一预测模型的预测误差为实际行程时间与第一预测行程时间之间的相对误差的绝对值,第二预测模型的预测误差为实际行程时间与第二预测行程时间之间的相对误差的绝对值。而在预测下一时段tc+1的路段行程时间时,所得到的针对时段tc的第一预测模型的预测误差可以作为该第一预测模型的一个新的历史预测误差,所得到的针对时段tc的第二预测模型的预测误差可以作为该第二预测模型的一个新的历史预测误差,以用于利用等式(4)和等式(5)动态更新第一预测模型和第二预测模型的权重系数,从而用于预测监测路段在下一时段tc+1的路段行程时间。
当在下一时段tc+1内、要预测再下一时段tc+2的路段行程时间时,按照上述同样的方法,可以首先确定出监测路段在时段tc+1内的实际路段行程时间。之后,根据该实际路段行程时间和在上一轮预测过程中由第一预测模型预测出的针对时段tc+1的第一预测路段行程时间,可以确定出第一预测模型针对时段tc+1的预测误差,以及根据该实际路段行程时间和在上一轮预测过程中由第二预测模型预测出的针对时段tc+1的第二预测路段行程时间,可以确定出第二预测模型针对时段tc+1的预测误差。而在预测时段tc+2的路段行程时间时,所得到的针对时段tc+1的第一预测模型的预测误差可以作为该第一预测模型的一个新的历史预测误差,所得到的针对时段tc+1的第二预测模型的预测误差可以作为该第二预测模型的一个新的历史预测误差,以用于利用等式(4)和等式(5)动态更新第一预测模型和第二预测模型的权重系数,从而用于预测监测路段在时段tc+2的路段行程时间。
针对后续时段的路段行程时间预测,都按照上述过程来动态调整模型的权重系数。这种权重系数的动态调整,可以保证在最终得出的预测结果中,精度高的预测模型得到的预测结果能够占有较大的权重,从而提高预测结果的准确性。
此外,如前所述,在步骤S1中,采样的是经地图匹配过程匹配到监测路段上的浮动车的位置数据。也就是说,在执行本发明提供的路段行程时间监测方法之前,需要对浮动车进行地图匹配,以将浮动车匹配到相应的路段上。可以采用多种地图匹配技术来进行浮动车的地图匹配。在本发明中,如图3所示,可以采用如下地图匹配过程:
首先,步骤S31,采集浮动车的GPS定位结果和电子地图数据。浮动车上均搭载GPS定位装置,因此,可以从浮动车的GPS定位装置采集该浮动车的GPS定位结果。此外,可以例如从GIS系统获取电子地图数据,该电子地图数据可以包括各路段的信息,例如,路段标识、路段名称、路段分布等等。
接下来,步骤S32,可以以所述浮动车的GPS定位结果为中心,在电子地图上确定一椭圆区域,该椭圆区域内的全部路段组成初始候选路段集。例如,如图4所示,假设浮动车的GPS定位结果表示该浮动车当前处于P点,那么,可以以该P点为中心,圈定一椭圆区域。例如,可以通过以下等式(6)~等式(8)来确定椭圆区域的大小:
其中,δx可以表示浮动车的GPS定位装置向东方向测量误差的标准差;δy可以表示浮动车的GPS定位装置向北方向测量误差的标准差;可以表示浮动车的GPS定位装置向东方向测量误差的方差;可以表示浮动车的GPS定位装置向北方向测量误差的方差;δxy可以表示协方差。上述参数δx、δy、和δxy可以从GPS定位装置的输出电文中获得。此外,a可以表示椭圆区域的长半轴;b可以表示椭圆区域的短半轴;可以表示椭圆区域的长半轴与正北方向夹角;可以为预设的扩展因子,可以根据需要的置信度来设定该扩展因子的具体值。
在确定出椭圆区域后,电子地图上、该椭圆区域内的全部路段组成初始候选路段集。通过设定椭圆区域来选择初始候选路段,可以有效减小初始候选路段的数量,降低匹配计算量,提高匹配效率。
接下来,步骤S33,可以根据浮动车的GPS定位结果确定浮动车与所述初始候选路段集中的每个初始候选路段的最短距离和方向角度差。即,确定P点距离每个初始候选路段的最短距离和方向角度差。如何确定所述最短距离和所述方向角度差的方法是本领域的技术人员公知的,对此,本发明不再进行赘述。
之后,步骤S34,可以根据浮动车与初始候选路段集中的每个初始候选路段的最短距离和方向角度差、以及浮动车的车速(例如,该车速可以包括在GPS定位结果中),从所述初始候选路段集中确定优选候选路段集。具体地,在该步骤中,可以首先根据浮动车与初始候选路段集中的每个初始候选路段的最短距离和方向角度差、以及浮动车的车速,确定出浮动车与初始候选路段集中的每个初始候选路段的匹配度。之后,可以按照匹配度由高到低的顺序对每个初始候选路段进行排序。最后,可以筛选出排名较高的若干初始候选路段,其余初始候选路段则被剔除。所筛选出的路段组成所述优选候选路段集。
在本发明中,可以例如采用模糊推理的方式确定浮动车与初始候选路段集中的每个初始候选路段的匹配度。其中,浮动车与初始候选路段的最短距离、方向角度差以及该浮动车的车速可以作为模糊推理模型的输入,而模糊推理模型的输出为该浮动车与初始候选路段的匹配度。可以根据大量历史样本数据来构建模糊推理模型,该构建技术是本领域的技术人员公知的,对此,本发明不再做详细阐述。
在确定出优选候选路段集后,步骤S35,可以利用网格匹配方法从所述优选候选路段集中确定与所述浮动车最匹配的路段,并将所述浮动车匹配到所述最匹配的路段上。在该步骤中,利用网格匹配方法来找出与所述浮动车最匹配的路段,并可以将浮动车映射到该最匹配的路段上,以得到浮动车在该最匹配的路段上的具体位置,由此完成浮动车的地图匹配。
例如,如图4所示,假设优选候选路段集中包括两个优选候选路段,分别表示为N1和N2。可以将电子地图栅格化,并确定浮动车的GPS定位结果(即,P点)所处的网格。有时会存在与P点最匹配的路段没有位于该P点所处网格中、而是位于该P点所处网格周围的网格中的情况,因此,为了防止遗漏并缩小搜索量,可以将用于搜索的网格区域大小从该P点所处的网格扩大到包括围绕该网格的周边网格,即,该网格区域共包括9个网格,如图4所示,P点所处网格和围绕该网格的8个周边网格。之后,搜索此9个网格中的优选候选路段,并计算P点到各网格中各优选候选路段的距离,确定此距离值最小所在网格并在此进行浮动车具体位置匹配。若搜索失败(即一个网格中存在多于两条的优选候选路段无法进行识别)则应细化网格,重新进行搜索匹配。按照这种方式,以图4为例,经网格匹配后,确定出优选候选路段N1为与浮动车的GPS定位结果(即,P点)最匹配的路段,并且该浮动车在该路段N1上的具体位置为P点垂直映射到路段N1的Q点。
在本发明中,先采用模糊推理方式进行候选路段筛选,之后再采用网格匹配方法进行匹配,可以简化路网的复杂度,减少搜索匹配时间,提高匹配速率与精度。
图5是根据本发明的实施方式的路段行程时间监测装置的框图。如图5所示,该装置可以包括:数据采样单元10,可以用于在一时段内,对经地图匹配过程匹配到监测路段上的浮动车的位置数据进行采样;第一处理单元20,可以用于根据所采样的位置数据,确定所述监测路段在所述时段内的实际行程时间;第一预测单元30,可以用于将所述时段内的实际行程时间输入至第一预测模型,得出所述监测路段在下一时段内的第一预测行程时间;第二预测单元40,可以用于将所述时段内的实际行程时间输入至第二预测模型,得出所述监测路段在下一时段内的第二预测行程时间;权重确定单元50,可以用于基于所述第一预测模型的历史预测误差和所述第二预测模型的历史预测误差来确定所述第一预测模型和所述第二预测模型的权重系数;以及第二处理单元60,可以用于基于所述第一预测模型和所述第二预测模型的权重系数,对所述第一预测行程时间和所述第二预测行程时间进行加权平均,得出所述监测路段在所述下一时段内的最终预测行程时间。
此外,虽然在附图中未示出,但是本发明提供的路段行程时间监测装置还可以包括地图匹配单元,可以用于对浮动车进行地图匹配。其中,所述地图匹配单元可以例如包括:数据采集模块,可以用于采集浮动车的GPS定位结果和电子地图数据;初始候选路段确定模块,可以用于以所述浮动车的GPS定位结果为中心,在电子地图上确定一椭圆区域,该椭圆区域内的全部路段组成初始候选路段集;计算模块,可以用于根据所述浮动车的GPS定位结果确定所述浮动车与所述初始候选路段集中的每个初始候选路段的最短距离和方向角度差;筛选模块,可以用于根据所述浮动车与所述初始候选路段集中的每个初始候选路段的最短距离和方向角度差、以及所述浮动车的车速,从所述初始候选路段集中确定优选候选路段集;匹配模块,可以用于利用网格匹配方法从所述优选候选路段集中确定与所述浮动车最匹配的路段,并将所述浮动车匹配到所述最匹配的路段上。
本发明的路段行程时间监测装置对应于路段行程时间监测方法,因此相同的内容不再赘述。
综上所述,通过采用两种预测模型来分别预测监测路段在下一时段的路段行程时间,并对两种预测模型得到的预测结果进行融合,从而得出最终的预测结果。由此,可以提高预测结果的准确性和可靠性。此外,可以根据第一预测模型和第二预测模型的历史预测误差动态调整第一预测模型和第二预测模型的权重系数,以保证精度相对高的预测模型得出的结果能够占有更大的权重,从而可以进一步提高预测结果的准确性。
例如,图6示出了基于本发明提供的路段行程时间监测方法及装置得到的监测路段的预测行程时间与实际行程时间、和采用单一预测模型得到的预测行程时间之间的结果对比图,其中,纵坐标表示行程时间,单位为秒;横坐标表示时段,单位为秒。如图6所示,相比于采用单一预测模型(例如,Kalman滤波模型或ARIMA模型)进行预测,基于本发明提供的路段行程时间监测方法及装置得到的监测路段的预测行程时间更接近于实际行程时间,并且,相较于单一预测模型,准确性可以大约提高(预测误差可减少15%~20%,平稳度可提高15%~20%)。
因此,通过本发明提供的路段行程时间监测方法及装置,能够为出行者及交通监管部门提供更为准确的路段行程时间的预测结果,从而便于出行者选择最佳行驶路段、以及便于交通监管部门及时采取相应的交通诱导措施。
以上结合附图详细描述了本发明的优选实施方式,但是,本发明并不限于上述实施方式中的具体细节,在本发明的技术构思范围内,可以对本发明的技术方案进行多种简单变型,这些简单变型均属于本发明的保护范围。
此外,尽管在附图中以特定顺序描述了本发明方法的操作,但是,这并非要求或者暗示必须按照该特定顺序来执行这些操作,或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地,可以省略某些步骤,将多个步骤合并为一个步骤执行,和/或将一个步骤分解为多个步骤执行。
本领域技术人员知道,本发明的实施方式可以实现为一种系统、装置、设备、方法或计算机程序产品。因此,本发明的实施方式可以具体实现为以下形式,即:完全的硬件、完全的软件(包括固件、驻留软件、微代码等),或者硬件和软件结合的形式。
另外需要说明的是,在上述具体实施方式中所描述的各个具体技术特征,在不矛盾的情况下,可以通过任何合适的方式进行组合。为了避免不必要的重复,本发明对各种可能的组合方式不再另行说明。
此外,本发明的各种不同的实施方式之间也可以进行任意组合,只要其不违背本发明的思想,其同样应当视为本发明所公开的内容。
Claims (10)
1.一种路段行程时间监测方法,其特征在于,该方法包括:
在一时段内,对经地图匹配过程匹配到监测路段上的浮动车的位置数据进行采样;
根据所采样的位置数据,确定所述监测路段在所述时段内的实际行程时间;
将所述时段内的实际行程时间分别输入至第一预测模型和第二预测模型,得出所述监测路段在下一时段内的第一预测行程时间和第二预测行程时间;
基于所述第一预测模型的历史预测误差和所述第二预测模型的历史预测误差来确定所述第一预测模型和所述第二预测模型的权重系数;以及
基于所述第一预测模型和所述第二预测模型的权重系数,对所述第一预测行程时间和所述第二预测行程时间进行加权平均,得出所述监测路段在所述下一时段内的最终预测行程时间。
2.根据权利要求1所述的方法,其特征在于,通过以下方式来确定所述监测路段在所述时段内的实际行程时间:
其中,Li,j表示在所述时段内采样的第j辆浮动车在所述监测路段上的第i个位置数据;
Li+1,j表示在所述时段内采样的第j辆浮动车在所述监测路段上的第i+1个位置数据;
Ti,j表示第j辆浮动车到达所述第i个位置数据所表示的位置时的时间;
Ti+1,j表示第j辆浮动车到达所述第i+1个位置数据所表示的位置时的时间;
nj表示在所述时段内采样的第j辆浮动车在所述监测路段上的位置数据总个数,并且nj≥2;以及
vj表示在所述时段内,第j辆浮动车在该监测路段上的平均车速;
T表示所述监测路段在所述时段内的实际行程时间;
L表示所述监测路段的总距;以及
m表示在所述时段内,所述监测路段上的浮动车总数。
3.根据权利要求1所述的方法,其特征在于,所述第一预测模型为卡尔曼(Kalman)滤波模型,以及所述第二预测模型为时间序列(ARIMA)模型。
4.根据权利要求1所述的方法,其特征在于,通过以下方式确定所述第一预测模型和所述第二预测模型的权重系数:
其中,ω1表示所述第一预测模型的权重系数;
ω2表示所述第二预测模型的权重系数;
e1(k)表示所述第一预测模型的第k个历史预测误差;
e2(k)表示所述第二预测模型的第k个历史预测误差;
K1表示所述第一预测模型的历史预测误差的总个数;以及
K2表示所述第二预测模型的历史预测误差的总个数。
5.根据权利要求1-4中任一权利要求所述的方法,其特征在于,所述地图匹配过程包括:
采集浮动车的GPS定位结果和电子地图数据;
以所述浮动车的GPS定位结果为中心,在电子地图上确定一椭圆区域,该椭圆区域内的全部路段组成初始候选路段集;
根据所述浮动车的GPS定位结果确定所述浮动车与所述初始候选路段集中的每个初始候选路段的最短距离和方向角度差;
根据所述浮动车与所述初始候选路段集中的每个初始候选路段的最短距离和方向角度差、以及所述浮动车的车速,从所述初始候选路段集中确定优选候选路段集;以及
利用网格匹配方法从所述优选候选路段集中确定与所述浮动车最匹配的路段,并将所述浮动车匹配到所述最匹配的路段上。
6.一种路段行程时间监测装置,其特征在于,该装置包括:
数据采样单元,用于在一时段内,对经地图匹配过程匹配到监测路段上的浮动车的位置数据进行采样;
第一处理单元,用于根据所采样的位置数据,确定所述监测路段在所述时段内的实际行程时间;
第一预测单元,用于将所述时段内的实际行程时间输入至第一预测模型,得出所述监测路段在下一时段内的第一预测行程时间;
第二预测单元,用于将所述时段内的实际行程时间输入至第二预测模型,得出所述监测路段在下一时段内的第二预测行程时间;
权重确定单元,用于基于所述第一预测模型的历史预测误差和所述第二预测模型的历史预测误差来确定所述第一预测模型和所述第二预测模型的权重系数;以及
第二处理单元,用于基于所述第一预测模型和所述第二预测模型的权重系数,对所述第一预测行程时间和所述第二预测行程时间进行加权平均,得出所述监测路段在所述下一时段内的最终预测行程时间。
7.根据权利要求6所述的装置,其特征在于,所述第一处理单元通过以下方式来确定所述监测路段在所述时段内的实际行程时间:
其中,Li,j表示在所述时段内采样的第j辆浮动车在所述监测路段上的第i个位置数据;
Li+1,j表示在所述时段内采样的第j辆浮动车在所述监测路段上的第i+1个位置数据;
Ti,j表示第j辆浮动车到达所述第i个位置数据所表示的位置时的时间;
Ti+1,j表示第j辆浮动车到达所述第i+1个位置数据所表示的位置时的时间;
nj表示在所述时段内采样的第j辆浮动车在所述监测路段上的位置数据总个数,并且nj≥2;以及
vj表示在所述时段内,第j辆浮动车在该监测路段上的平均车速;
T表示所述监测路段在所述时段内的实际行程时间;
L表示所述监测路段的总距;以及
m表示在所述时段内,所述监测路段上的浮动车总数。
8.根据权利要求6所述的装置,其特征在于,所述第一预测模型为卡尔曼(Kalman)滤波模型,以及所述第二预测模型为时间序列(ARIMA)模型。
9.根据权利要求6所述的装置,其特征在于,所述权重确定单元通过以下方式确定所述第一预测模型和所述第二预测模型的权重系数:
其中,ω1表示所述第一预测模型的权重系数;
ω2表示所述第二预测模型的权重系数;
e1(k)表示所述第一预测模型的第k个历史预测误差;
e2(k)表示所述第二预测模型的第k个历史预测误差;
K1表示所述第一预测模型的历史预测误差的总个数;以及
K2表示所述第二预测模型的历史预测误差的总个数。
10.根据权利要求6-9中任一权利要求所述的装置,其特征在于,该装置还包括地图匹配单元,用于对浮动车进行地图匹配,其中,所述地图匹配单元包括:
数据采集模块,用于采集浮动车的GPS定位结果和电子地图数据;
初始候选路段确定模块,用于以所述浮动车的GPS定位结果为中心,在电子地图上确定一椭圆区域,该椭圆区域内的全部路段组成初始候选路段集;
计算模块,用于根据所述浮动车的GPS定位结果确定所述浮动车与所述初始候选路段集中的每个初始候选路段的最短距离和方向角度差;
筛选模块,用于根据所述浮动车与所述初始候选路段集中的每个初始候选路段的最短距离和方向角度差、以及所述浮动车的车速,从所述初始候选路段集中确定优选候选路段集;以及
匹配模块,用于利用网格匹配方法从所述优选候选路段集中确定与所述浮动车最匹配的路段,并将所述浮动车匹配到所述最匹配的路段上。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510286096.7A CN104900061B (zh) | 2015-05-29 | 2015-05-29 | 路段行程时间监测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510286096.7A CN104900061B (zh) | 2015-05-29 | 2015-05-29 | 路段行程时间监测方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104900061A true CN104900061A (zh) | 2015-09-09 |
CN104900061B CN104900061B (zh) | 2017-08-08 |
Family
ID=54032701
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510286096.7A Expired - Fee Related CN104900061B (zh) | 2015-05-29 | 2015-05-29 | 路段行程时间监测方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104900061B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105489010A (zh) * | 2015-12-29 | 2016-04-13 | 中国城市规划设计研究院 | 一种快速道路行程时间可靠度监测分析系统及方法 |
CN105702036A (zh) * | 2016-03-23 | 2016-06-22 | 深圳市金溢科技股份有限公司 | 一种车辆行驶时间计算方法、装置及系统 |
CN106096767A (zh) * | 2016-06-07 | 2016-11-09 | 中国科学院自动化研究所 | 一种基于lstm的路段行程时间预测方法 |
CN106682754A (zh) * | 2015-11-05 | 2017-05-17 | 阿里巴巴集团控股有限公司 | 事件发生概率预测方法及装置 |
CN109215346A (zh) * | 2018-10-11 | 2019-01-15 | 平安科技(深圳)有限公司 | 一种交通通行时间的预测方法、存储介质和服务器 |
CN109974735A (zh) * | 2019-04-08 | 2019-07-05 | 腾讯科技(深圳)有限公司 | 到达时间的预估方法、装置及计算机设备 |
CN111667689A (zh) * | 2020-05-06 | 2020-09-15 | 浙江师范大学 | 车辆旅行时间预测的方法、设备和计算机设备 |
CN112415892A (zh) * | 2020-11-09 | 2021-02-26 | 东风汽车集团有限公司 | 一种汽油机起动标定控制参数优化方法 |
CN112489418A (zh) * | 2020-10-22 | 2021-03-12 | 浙江交通职业技术学院 | 一种基于道路路段行程时间预测模型的路段行程时间动态误差修正方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101388143A (zh) * | 2007-09-14 | 2009-03-18 | 同济大学 | 基于浮动车数据的公交到站时间预测方法及系统 |
CN102184453A (zh) * | 2011-05-16 | 2011-09-14 | 上海电气集团股份有限公司 | 基于模糊神经网络和支持向量机的风电功率组合预测方法 |
CN103280109A (zh) * | 2013-06-08 | 2013-09-04 | 北京云星宇交通工程有限公司 | 旅行时间的获取方法及装置、预测系统 |
-
2015
- 2015-05-29 CN CN201510286096.7A patent/CN104900061B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101388143A (zh) * | 2007-09-14 | 2009-03-18 | 同济大学 | 基于浮动车数据的公交到站时间预测方法及系统 |
CN102184453A (zh) * | 2011-05-16 | 2011-09-14 | 上海电气集团股份有限公司 | 基于模糊神经网络和支持向量机的风电功率组合预测方法 |
CN103280109A (zh) * | 2013-06-08 | 2013-09-04 | 北京云星宇交通工程有限公司 | 旅行时间的获取方法及装置、预测系统 |
Non-Patent Citations (4)
Title |
---|
姚恩建,左婷: "基于低频浮动车数据的实时地图匹配算法", 《北京工业大学学报》 * |
宋曰聪,胡伟,张涛: "基于遗传算法的交通流量组合预测研究", 《微计算机信息》 * |
张硕,孙剑,李克平: "路径行程时间的组合预测方法研究", 《交通信息与安全》 * |
邹珍: "基于GPS的浮动车数据与实地图匹配的算法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106682754A (zh) * | 2015-11-05 | 2017-05-17 | 阿里巴巴集团控股有限公司 | 事件发生概率预测方法及装置 |
CN105489010A (zh) * | 2015-12-29 | 2016-04-13 | 中国城市规划设计研究院 | 一种快速道路行程时间可靠度监测分析系统及方法 |
CN105489010B (zh) * | 2015-12-29 | 2019-01-04 | 中国城市规划设计研究院 | 一种快速道路行程时间可靠度监测分析系统及方法 |
CN105702036A (zh) * | 2016-03-23 | 2016-06-22 | 深圳市金溢科技股份有限公司 | 一种车辆行驶时间计算方法、装置及系统 |
CN106096767A (zh) * | 2016-06-07 | 2016-11-09 | 中国科学院自动化研究所 | 一种基于lstm的路段行程时间预测方法 |
CN109215346A (zh) * | 2018-10-11 | 2019-01-15 | 平安科技(深圳)有限公司 | 一种交通通行时间的预测方法、存储介质和服务器 |
CN109974735A (zh) * | 2019-04-08 | 2019-07-05 | 腾讯科技(深圳)有限公司 | 到达时间的预估方法、装置及计算机设备 |
CN111667689A (zh) * | 2020-05-06 | 2020-09-15 | 浙江师范大学 | 车辆旅行时间预测的方法、设备和计算机设备 |
CN112489418A (zh) * | 2020-10-22 | 2021-03-12 | 浙江交通职业技术学院 | 一种基于道路路段行程时间预测模型的路段行程时间动态误差修正方法 |
CN112415892A (zh) * | 2020-11-09 | 2021-02-26 | 东风汽车集团有限公司 | 一种汽油机起动标定控制参数优化方法 |
CN112415892B (zh) * | 2020-11-09 | 2022-05-03 | 东风汽车集团有限公司 | 一种汽油机起动标定控制参数优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104900061B (zh) | 2017-08-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104900061B (zh) | 路段行程时间监测方法及装置 | |
EP3848730A1 (en) | Positioning method, apparatus and device, and computer-readable storage medium | |
CN104574967B (zh) | 一种基于北斗的城市大面积路网交通感知方法 | |
CN105321345B (zh) | 一种基于ARIMA模型和kalman滤波的道路交通流预测方法 | |
CN112927513B (zh) | 一种实时在线交通仿真方法及系统 | |
CN112669594A (zh) | 交通路况预测的方法、装置、设备及存储介质 | |
Stipancic et al. | Impact of congestion and traffic flow on crash frequency and severity: application of smartphone-collected GPS travel data | |
CN104851293A (zh) | 一种基于地点车速的路段交通拥堵指数评价方法 | |
CN105225541A (zh) | 基于空管历史数据挖掘的短时航迹预测方法 | |
CN113284369B (zh) | 一种基于ads-b实测航路数据的预测方法 | |
US11908317B2 (en) | Real-time traffic safety management system | |
CN111667689B (zh) | 车辆旅行时间预测的方法、设备和计算机设备 | |
CN116957331A (zh) | 风险旅客流量范围预测方法和装置 | |
CN111582919A (zh) | 一种航班客座率预测方法和系统 | |
CN113053106A (zh) | 基于多传感器信息融合的交通状态估计方法和装置 | |
CN111915880A (zh) | 一种高速公路异常事件事发点有效通行能力估计方法 | |
Dong et al. | An identification model of urban critical links with macroscopic fundamental diagram theory | |
CN108665705B (zh) | 交通诱导统计模型校对方法及装置 | |
CN112116167A (zh) | 数据处理方法、数据处理装置、存储介质和电子设备 | |
Berzina et al. | Development of a post-processing automation procedure for the GPS-based travel time data collection technique | |
Khattak et al. | TabNet-SHAP: A Framework to Estimate Wind Shear-Induced Aviation Turbulence in the Airport Runway Vicinity | |
WO2023032251A1 (ja) | 情報処理装置、情報処理方法、プログラムおよび移動体 | |
CN118378901B (zh) | 一种高速公路交通事故影响程度动态评价方法及设备 | |
CN111524354B (zh) | 基于语言模型的城市交通网络路径选择行为的预测方法、系统、介质及设备 | |
CN109284869B (zh) | 一种基于浮动车数据城市交叉口流量估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170808 Termination date: 20190529 |