CN113380025A - 一种车辆出行量预测模型构建方法及预测方法和系统 - Google Patents
一种车辆出行量预测模型构建方法及预测方法和系统 Download PDFInfo
- Publication number
- CN113380025A CN113380025A CN202110590338.7A CN202110590338A CN113380025A CN 113380025 A CN113380025 A CN 113380025A CN 202110590338 A CN202110590338 A CN 202110590338A CN 113380025 A CN113380025 A CN 113380025A
- Authority
- CN
- China
- Prior art keywords
- road section
- traffic
- vehicle
- data
- time slot
- 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 32
- 238000010276 construction Methods 0.000 title claims abstract description 15
- 239000011159 matrix material Substances 0.000 claims description 63
- 239000013598 vector Substances 0.000 claims description 48
- 238000012549 training Methods 0.000 claims description 9
- 238000007781 pre-processing Methods 0.000 claims description 8
- AGVAZMGAQJOSFJ-WZHZPDAFSA-M cobalt(2+);[(2r,3s,4r,5s)-5-(5,6-dimethylbenzimidazol-1-yl)-4-hydroxy-2-(hydroxymethyl)oxolan-3-yl] [(2r)-1-[3-[(1r,2r,3r,4z,7s,9z,12s,13s,14z,17s,18s,19r)-2,13,18-tris(2-amino-2-oxoethyl)-7,12,17-tris(3-amino-3-oxopropyl)-3,5,8,8,13,15,18,19-octamethyl-2 Chemical compound [Co+2].N#[C-].[N-]([C@@H]1[C@H](CC(N)=O)[C@@]2(C)CCC(=O)NC[C@@H](C)OP(O)(=O)O[C@H]3[C@H]([C@H](O[C@@H]3CO)N3C4=CC(C)=C(C)C=C4N=C3)O)\C2=C(C)/C([C@H](C\2(C)C)CCC(N)=O)=N/C/2=C\C([C@H]([C@@]/2(CC(N)=O)C)CCC(N)=O)=N\C\2=C(C)/C2=N[C@]1(C)[C@@](C)(CC(N)=O)[C@@H]2CCC(N)=O AGVAZMGAQJOSFJ-WZHZPDAFSA-M 0.000 description 7
- 230000000694 effects Effects 0.000 description 7
- 101001095088 Homo sapiens Melanoma antigen preferentially expressed in tumors Proteins 0.000 description 4
- 102100037020 Melanoma antigen preferentially expressed in tumors Human genes 0.000 description 4
- 230000009471 action Effects 0.000 description 4
- 230000000052 comparative effect Effects 0.000 description 4
- 238000013528 artificial neural network Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000005065 mining Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000011218 segmentation Effects 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003203 everyday effect Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 description 1
- 241000282414 Homo sapiens Species 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000002354 daily effect Effects 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 238000002922 simulated annealing Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
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/0129—Traffic data processing for creating historical data or processing based on historical data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Game Theory and Decision Science (AREA)
- Chemical & Material Sciences (AREA)
- Development Economics (AREA)
- Analytical Chemistry (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明公开了一种车辆出行量预测的模型构建方法及预测方法和系统。所公开方案构建路段出行量OD次数关系图,路段出行量局部关系图、路段出行量全局关系图共同作为空间特征,进行路段车辆出行量预测,并且也考虑外部因素影响,同时采用GCN和TCN混合的深度预测模型进行模型构建与未来出行量的预测。所公开方案提升了车辆出行量预测的准确度,可用于区域内规模车辆出行量的预测。
Description
技术领域
本发明涉及车辆交通流预测技术,具体涉及一种车辆出行量预测模型构建方法及车辆出行量预测方法。
背景技术
实现未来时段内较为准确的车辆(如出租车、网约车等)需求预测可以帮助车辆运营部门进行科学合理的车辆调度与管理,进一步降低城市范围内车辆的空驶率,节约车辆运营成本,同时有效降低乘客打车时由于候车所付出的时间成本,实现司乘双赢。
车辆出行量数据具有复杂的时空依赖关系,现有方法难以有效提取这种高度抽象的非线性关系,并且大多是对经网格划分后的网格内车辆出行量进行预测,这在一定程度上破坏了出行量数据内部潜在的空间关系,在空间粒度上不够精细。在基于图网络的车辆需求预测方法中,构图方式较为单一,不利于充分提取出行量数据的空间关系。
发明内容
针对现有技术的缺陷或不足,本发明的一方面提供了一种车辆出行量预测模型构建方法。
为此,本发明提供的车辆出行量预测模型构建方法包括:
(1)数据预处理,包括:
采集目标区域待预测时隙之前(不包括待预测时隙)多个历史时隙内的车辆轨迹数据,并将车辆轨迹数据与目标区域的路段进行匹配,得到原始数据;
采集每时间段的外部因素获得外部因素数据,所述外部因素包括星期几、空气质量、温度和天气中一种或两种以上因素,所述最后一个历史时隙之前(不包括最后一个历史时隙)的所有历史时隙被划分为多个所述时间段;
(2)利用多个历史时隙内最后一个历史时隙之前的原始数据构建目标区域的路段出行量局部关系图、路段出行量全局关系图和路段出行量OD次数关系图;利用多个历史时隙内最后一个历史时隙之前的外部因素数据构建外部因素特征矩阵;
所述路段出行量局部关系图包括路段出行量矩阵、第一邻接关系矩阵和目标区域的边集;
所述路段出行量矩阵的任一行向量为任一路段i的出行量向量,该路段i的出行量向量由每个历史时隙内该路段i的车辆出行量构成,路段i为目标区域内的任一路段;
所述第一邻接关系矩阵的任一元素aij为路段i与路段j的相邻关系,如路段i与路段j互相交叉或首尾相接,则aij=1,否则aij=0,任一行向量由任一路段i与自身及其他路段的相邻关系构成,路段j为目标区域内的任一路段;
所述路段出行量全局关系图包括所述路段出行量矩阵、第二邻接关系矩阵和所述目标区域的边集;
所述第二邻接关系矩阵的任一元素bij为路段i的出行量向量与路段j的出行量向量相似度,任一行向量由任一路段i的出行量向量与自身的出行量向量及其他路段的出行量向量的相似度构成;
所述路段出行量OD次数关系图包括所述路段出行量矩阵、第三邻接关系矩阵和所述目标区域的边集;
所述第三邻接关系矩阵的任一元素cij为归一化后的路段i到路段j的OD次数,任一行向量由归一化后的路段i到自身的OD次数及归一化后的路段i到其他各路段的OD次数,
xij为目标区域最后一个历史时隙之前路段i与路段j间的OD次数或者xij为所有时间段内的路段i与路段j间OD次数的平均值,R为目标区域的所有路段集合;
所述外部因素特征矩阵由每个时间段的外部因素特征向量构成,任一时间段的外部因素特征向量由该时间段的外部因素构成;
(3)以所述路段出行量局部关系图、路段出行量全局关系图和路段出行量OD次数关系图和外部因素特征矩阵作为输入,以最后一个历史时隙内各路段车辆出行量为标签,对TCN和GCN混合深度网络进行训练,得到预测模型;
所述TCN和GCN混合深度网络由四个通道和卷积层构成,其中三个通道各自均包含依次设置的一层TCN、两层GCN和一层TCN,这三个通道内的一个通道的输入为所述路段出行量全局关系图、一个通道的输入为路段出行量局部关系图、另一个的输入为所述出行量OD次数关系图;第四个通道包括全连接层,该第四个通道的输入为外部因素特征矩阵;将所述四个通道的输出进行融合后,输入所述卷积层。
可选的,所述多个历史时隙的总时长为数周、数个月或数年。
可选的,所述时间段时长为数个小时、半天、一天或多天。
可选的,所述每个历史时隙为数秒或数分钟。
可选的,所述待预测时隙时长为数分钟、数小时、半天或一天。
可选的,对所述各时间段的外部因素中的各因素分别采用独热向量进行编码,之后对各因素的编码进行级联,得到总独热向量,该总独热向量为相应时间段的外部因素特征向量。
本发明还提供了一种车辆出行量预测方法。所提供的预测方法是对上述待预测时隙的车辆出行量进行预测,包括
(一)数据采集与预处理,包括:
采集目标区域待预测时隙之前多个历史时隙内的车辆轨迹数据,并将车辆轨迹数据与目标区域的路段进行匹配,得到原始数据;
采集每时间段的外部因素获得外部因素数据,所述外部因素包括星期几、空气质量、温度和天气中一种或两种以上因素,所述多个历史时隙被划分为多个所述时间段;
(二)采用上述步骤(2)方法利用多个历史时隙内的原始数据构建目标区域的路段出行量局部关系图、路段出行量全局关系图和路段出行量OD次数关系图;利用所述多个历史时隙内的外部因素数据构建外部因素特征矩阵;
(三)将步骤(二)构建的路段出行量局部关系图、路段出行量全局关系图、路段出行量OD次数关系图和外部因素特征矩阵输入上述方法构建的模型中,输出待预测时隙内的车辆出行量预测值。
又一种车辆出行量预测方法,对上述待预测时隙的下一个待预测时隙内的车辆出行量进行预测,包括:
以上述步骤(二)构建的路段出行量局部关系图、路段出行量全局关系图、路段出行量OD次数关系图和外部因素特征矩阵作为输入,以待预测时隙内的车辆出行量为输出对上述方法构建的模型进行更新;
之后采用上步骤(一)、(二)所述方法及更新后的模型对下一个待预测时隙的车辆出行量进行预测。
本发明同时提供了一种车辆出行量预测系统。所提供的预测系统包括数据预处理模块、数据构建模块和预测及更新模块,所述数据预处理模块用于执行上述步骤(一),所述数据构建模块用于执行上所述步骤(二),所述数据预测及更新模块用于对新待预测时隙的预测模型进行更新并利用更新后的模型对新待预测时隙的车辆出行量进行预测。
与现有技术相比,本发明的效果是:
本发明主要是考虑了车辆出行量数据内在的多种空间特征关系,另外考虑了外部因素对车辆需求的影响,同时设计出一种新的车辆出行量预测网络结构,方法整体在充分考虑了车辆出行量数据内在的多种空间特征关系,联合外部因素对深度网络进行训练,具有预测结果更准确的优点。
附图说明
图1为本发明的混合深度网络结构示例图;
图2为实施例6中时隙大小对MFGTN网络的影响图;
图3为实施例7中预测步长对网络的影响(RMSE)图,横坐标1-5分别代表预测步长为5min、10min、15min、20min、25min的对比方案;
图4为实施例7中预测步长对网络的影响(MAE)图,横坐标1-5分别代表预测步长为5min、10min、15min、20min、25min的对比方案;
图5为实施例7中预测步长对网络的影响(MAPE)图,横坐标1-5分别代表预测步长为5min、10min、15min、20min、25min的对比方案。
具体实施方式
除非有特殊说明,本文中的术语或方法根据相关领域普通技术领域人员的常规认识理解或采用已有相应方法实现。
本发明方法中的车辆轨迹数据的采集及数据内容对于本领域技术人员来讲是公知的,所采集的数据内容至少包括车牌号、GPS采样时间、GPS地理位置、载客状态等。本发明所采集的车辆轨迹数据是针对目标区域例如一个城市、一个地区的车辆轨迹数据,并且是其一定历史范围及待预测时隙之前(不包括待预测时隙)多个历史时隙内的数据,该历史时间范围满足预测方法的数据量要求即可,例如一周、多周、半个月、一个月、多个月或一年、多年等。并且,在车辆轨迹数据采集过程中,需对数据进行清洗和校正,数据清洗目的是去除无效数据、格式错误和异常数据、在某个时间段内GPS地理位置没有变化以及不在目标区域内的数据,得到合格数据。校正的目的对数据进行地图匹配,对其中的GPS数据的误差进行校正;数据清洗和校正的方法手段对于本领域技术人员来讲是公知的。
本发明所述历史时隙大小可取数秒、数分钟,也根据待预测时隙确定,例如:一次待预测时隙取5分钟,送入一周的数据进行训练,一次预测可预测未来一个时隙(5min)的出行量信息。时隙取的越小,预测的准确率就越高;时隙增大,预测准确率降低。
本发明的待预测时隙时长即预测步长可根据预测需求来改变,例如:在网络训练好后,可以预测未来一个步长的出行量信息。预测步长取的越小,预测的准确率越高;预测步长取的越大,网络的预测误差均逐渐增大,这是由多次单步预测的误差积累造成的。
本发明的对车辆轨迹数据与目标区域的路段进行匹配是将车辆轨迹数据中的坐标信息与路段的地理坐标信息进行匹配,如可采用地图匹配。进一步,对于一定时间范围内车辆流量小于某一特定次数的路段通常不被视为主要路段,这些路段对于车辆流量预测的影响很小,可将其删除。在具体实施方案中,为方便数据分析,对各路段设置ID,路段ID以取自然数,也可以设置为相应便于数据处理的索引号等。
本发明中采集所述历史时间范围内各时间段的外部因素,所述外部因素主要考虑星期几、天气、空气质量、气温四种因素中一种或几种。其中,星期几指星期一至星期七;天气取常见的天气预报中的天气,如晴、多云、阴、小雨、中雨、大雨;气温根据目标区域常年历史温度确定相应温度范围,为更使预测更精确,根据人类对温度的喜好可将温度范围划分为多个区段,也可等间隔划分为对各区段;空气质量可使用空气质量指数(Air QualityIndex,AQI)来反映,一般常见AQI范围取0~150,根据空气质量相关标准划分为相应区段,也可等间隔划分为多个区段。考虑到外壁因素变动时长的特点,本发明采集一定时间段的外部因素,所述时间段可以取半天、一天、多天,常见如每天的外部因素。在构建外部因素特征矩阵时,可采用独热向量编码方式构建矩阵,例如对于星期几可以使用7位独热向量进行编码,第几位为1就表示当前向量代表星期几;天气使用6位独热向量进行编码,分别表示晴、多云、阴、小雨、中雨、大雨六种天气状况;空气质量AQI范围取0~150,以10为间隔进行分段,则可用15位独热向量对空气质量指数进行编码;气温范围为10℃~28℃,以2℃为间隔进行分段,则可用9位独热向量对气温进行编码。之后将上述四种编码进行级联,得到相应时间段内的37位的独热向量,之后构建外部因素特征矩阵。
本发明所述车辆出行量和OD次数是指,若一辆车的上车点位于路段i上,下车点位于路段j上,则认为发生了一次路段i到路段j的出行,在一定时间范围内路段i与路段j之间的总出行次数即为该时间段两条路段之间的OD次数,某一路段的出行总次数则视为该时间范围内本路段的车辆出行量。具体可通过上客事件和下客事件检测算法,可提取到上车点与下车点,可分时间段计算各路段之间的OD次数。
在对数据相似度度量时,为了对样本间的相似程度进行直观的、数量化的描述,通常通过数据间的距离来反映相似度,数据间的距离大小与样本间的相似度成反比关系。本发明的任意两个路段的出行量向量相似度可采用欧氏距离法、余弦相似度法、DTW算法或FastDTW算法计算。
本发明所述目标区域边集为目标区域内任一路段与其他路段之间的关系构成的矩阵,其中两条路段(如路段i与路段j)之间的关系至少包括连接方式(单行或双行)、距离(欧氏距离或真实路网中的距离)或是路段之间的相似程度。对于任一目标区域的边集本领域技术人员可以通过地图匹配提取或直接采用现有技术中已经提取的结果。
本发明的TCN和GCN混合深度网络由四个通道(channel)和卷积层构成,其中三个通道各自均包含依次设置的一层TCN、两层GCN和一层TCN、构成,这三个通道内的一个通道送入路段出行量全局关系图,一个通道送入路段出行量局部关系图,另一个送入出行量OD次数关系图;第四个通道包含全连接层FC,其中送入的是外部因素特征矩阵;将各个通道的输出矩阵进行融合concat(级联)后,接着送入卷积层conv。
本发明的混合深度网络利用TCN的时间序列处理分析能力,提取路段的出行量的时间序列特征;堆叠两层GCN实现空间特征挖掘,对邻接矩阵中表示的路段之间的空间关系进行特征提取,并且设置所得输出矩阵大小为1×R;最后一层TCN,得到输出矩阵大小仍为1×R,从而完成车辆需求数据时间特征和空间特征的提取;通过使用全连接层进行特征映射,其作用是完成特征空间的变换并且将外部因素的信息提取整合,最后设置输出矩阵大小也为1×R;接着,通过卷积层,将各个通道的输出特征进行提取,综合考虑不同的时空特征的关系,将信息融合。本发明的混合深度网络中三个通道的中各层的操作相同,但作用不同,送入路段出行量全局关系图的通道作用是提取路段出行量序列与路段全局空间特征的关系,送入路段出行量局部关系图的通道作用是提取路段出行量序列与路段局部空间特征的关系,送入出行量OD次数关系图的通道作用是提取路段出行量序列与路段之间的OD次数的关系。
下面结合具体的实施例对本发明做进一步的详细说明,所述是对本发明的解释而不是限定。
实施例1:
本实施例使用的数据集来自2016年10月17日到2016年10月23间采集到的GPS数据共包括7天,其中前6天的数据作为训练集,最后1天的数据作为测试集;该实施例并以5分钟作为一个历史时隙时长,则每天可划分为288个历史时隙,总历史时隙个数为2016;该实施例的待预测时隙时长即预测步长为5min。
本实施例所用车辆轨迹数据集由西安市交通局提供,为西安市车辆真实轨迹数据,数据采集时涉及的车辆超过10000辆,车辆运行区域为西安市三环以内的20000余条路段。每台车辆配备GPS采集设备,每30秒进行一次车辆运行信息采集,每天的运营时间平均超过18个小时。
首先对车辆轨迹数据(包括车辆表示标号、车牌号、GPS采样时间、GPS地理位置、载客状态)进行清洗,删除重复的记录项和无效的GPS记录,并对研究路段的范围进行限定,然后通过地图匹配对采集到的GPS轨迹点进行校正,保证GPS点处于路段上;同时对区域内的路段进行编号,城市路网中连接两个道路交叉口的路为一个路段;该实施例选取日均车辆出行量大于15次的路段进行研究,筛选得到851个路段,进行路段标号;
获取历史时范围段内每天的星期几、天气状况、空气质量、气温四种外部因素;
通过对城市路网内车辆乘客的上车记录点和下车记录点进行挖掘,进而分时隙计算每个时隙内各路段间的OD次数,统计得到各时隙内每个路段的车辆出行量;并且根据原始路网拓扑结构提取每个路段的邻接路段,依据本发明的步骤(3)-(6)构建路段出行量局部关系图、路段出行量全局关系图、路段出行量OD次数关系图和外部因素特征矩阵;
在构建路段出行量全局关系图式,该实施例采用文献Salvador S,Chan P.Towardaccurate dynamic time warping in linear time and space[J].Intelligent DataAnalysis,2007,11(5):561-580.公开的FastDTW算法计算路段i的车辆出行量序列Xi与路段j间的车辆出行量序列Xj相似度ωij计算公式如式:
该实施例的第三邻接关系矩阵的任一元素cij,
其中,xij为目标区域历史时间范围内路段i与路段j间的日均OD次数;
该实施例的目标区域的边集根据文献:殷剑宏,吴开亚.图论及其算法[M].中国科学技术大学出版社,2003.中公开的内容提取;
该实施例所构建的外部因素特征进行提取,主要考虑每天的星期几、天气状况、空气质量和气温四种因素,其中,星期几使用7位独热向量进行编码,第几位为1就表示当前向量代表星期几;天气状况使用6位独热向量进行编码,分别表示晴、多云、阴、小雨、中雨、大雨六种天气状况;空气质量则使用空气质量指数(Air Quality Index,AQI)来反映,表1为研究时段内的空气质量等级示意表,AQI范围取0~150,以10为间隔进行分段,则可用15位独热向量对空气质量指数进行编码;研究时段内气温范围为10℃~28℃,以2℃为间隔进行分段,则可用9位独热向量对气温进行编码;将上述四种编码进行级联,得到37位的独热向量,表示外部因素特征矩阵。
表1
如表1所示,为外部因素特征矩阵中的一个行向量,其表示信息为当前为星期一、晴天、AQI指数为70~80之间,空气质量为70-80(等级为良)、气温为18~20℃。
表2研究时段内的空气质量等级示意表
依照步骤(7)所述方法,进行网络训练,该实施例的TCN和GCN混合深度网络由四个通道和一个卷积层构成,其中第四个通道包含两个全连接层,网络结构如图1所示。
本实施例中,初始学习率设置为0.0001;
本实施例中预测网络基于TensorFlow深度学习框架搭建,采用基于自适应学习率的Adam梯度下降算法进行网络训练批次大小设置为32,迭代次数设置为100,同样采用了提前终止策略避免过拟合。
网络训练的目标是使得各路段车辆出行量的真实值与经过网络处理所得的各路段车辆出行量的预测值间的误差最小化。网络训练过程中所采用的损失函数如式:
模型的评价指标为RMSE、MAE和MAPE。本实施例中,出租出出行量预测对比实验结果如下表所示:
表3对比实验结果
MLP:Guo X,Zhu Q.A Traffic Flow Forecasting Model Based on BP NeuralNetwork[C]//2009 2nd International Conference on Power Electronics andIntelligent Transportation System(PEITS).IEEE,2009:311-314.
SVR:Hong W.Traffic flow forecasting by seasonal SVR with chaoticsimulated annealing algorithm[J].Neurocomputing,2011,74(12–13):2096-2107.
LSTM:Fu R,Zhang Z,Li L.Using LSTM and GRU neural network methods fortraffic flow prediction[C]//2016 31st Youth Academic Annual Conference ofChinese Association of Automation(YAC).IEEE,2016:324-328.
GCN:Wu Z,Pan S,Chen F,et al.A comprehensive survey on graph neuralnetworks[J].arXiv preprint arXiv:1901.00596,2019.
T-GCN:Zhao L,Song Y,Zhang C,et al.T-GCN:A temporal graphconvolutional network for traffic prediction[J].IEEE Transactions onIntelligent Transportation Systems,2019,21(9):3848-3858.
MGLN:Multitask Learning and GCN-Based Taxi Demand Prediction for aTraffic Road Network[J].Sensors,2020,20(13):3776.
实施例2:
该实施例与实施例1不同的是,该实施例所构建的外部因素特征矩阵考虑星期几。
实施例3:
该实施例与实施例1不同的是,该实施例所构建的外部因素特征矩阵考虑天气状况。
实施例4:
该实施例与实施例1不同的是,该实施例所构建的外部因素特征矩阵考虑空气质量。
实施例5:
该实施例与实施例1不同的是,该实施例所构建的外部因素特征矩阵考虑气温。
对比例:
该对比例与实施例1不同的是:不考虑外部因素,且网络训练为七种策略,即路段出行量全局关系图、路段出行量局部关系图、路段出行量OD次数关系图选其一、其二和三者。
上述实施例与对比例预测效果对比如表4所示:
表4
表4中,Glocal、Gglobal和God分别表示路段出行量全局关系图、路段出行量局部关系图、路段出行量OD次数关系图;Xext表示实施例1的外部因素特征矩阵表示为;Xd、Xw、Xa和Xt分别表示星期因素矩阵、天气因素矩阵、空气质量因素矩阵和气温因素矩阵,即对应实施例2-5。
结果显示:当输入包括路段出行量全局关系图、路段出行量局部关系图、路段出行量OD次数关系图三种图时,网络预测的RMSE、MAE和MAPE分别为0.9462、0.6335、7.5553%,可以看出充分考虑车辆出行量数据间的多种空间关系有利于促进网络对空间特征的挖掘,提升网络的预测能力;
此外,可以看出,分别将星期因素矩阵、天气因素矩阵、空气质量因素矩阵和气温因素矩阵逐个加入网络输入也可以降低预测误差,其中加入天气因素矩阵的效果最明显。进一步地,实施例1加入了外部因素特征矩阵作为输入,充分考虑了星期几、天气状况、空气质量、气温等外部因素对车辆出行量的影响,从而有效降低了网络预测误差。
实施例6:
本实施例中,与实施例1的时隙设置大小为5分钟不同的是,分别设置15分钟和30分钟进行时隙划分。图2为时隙大小对本发明网络的影响实验结果图。
实施例7:
本实施例中,使用实施例1网络实现路段车辆出行量预测,以探究预测步长对网络表现的影响,并与使用MGLN网络进行多步预测的结果进行对比,各方案分别设置待预测时隙时长即预测步长为5min、10min、15min、20min、25min进行对比实验。
图3、图4和图5分别为使用MGLN网络和实施例1网络实现多步预测的RMSE、MAE和MAPE结果图,可以看出,由于实施例1网络与MGLN网络相比,在不同的预测步长下,实施例1网络的预测误差均低于MGLN网络。
Claims (9)
1.一种车辆出行量预测模型构建方法,其特征在于,方法包括:
(1)数据预处理,包括:
采集目标区域待预测时隙之前多个历史时隙内的车辆轨迹数据,并将车辆轨迹数据与目标区域的路段进行匹配,得到原始数据;
采集每时间段的外部因素,获得外部因素数据,所述外部因素包括星期几、空气质量、温度和天气中一种或两种以上因素,所述最后一个历史时隙之前的所有历史时隙被划分为多个所述时间段;
(2)利用多个历史时隙内最后一个历史时隙之前的原始数据构建目标区域的路段出行量局部关系图、路段出行量全局关系图和路段出行量OD次数关系图;利用多个历史时隙内最后一个历史时隙之前的外部因素数据构建外部因素特征矩阵;
所述路段出行量局部关系图包括路段出行量矩阵、第一邻接关系矩阵和目标区域的边集;
所述路段出行量矩阵的任一行向量为任一路段i的出行量向量,该路段i的出行量向量由每个历史时隙内该路段i的车辆出行量构成,路段i为目标区域内的任一路段;
所述第一邻接关系矩阵的任一元素aij为路段i与路段j的相邻关系,如路段i与路段j互相交叉或首尾相接,则aij=1,否则aij=0,任一行向量由任一路段i与自身及其他路段的相邻关系构成,路段j为目标区域内的任一路段;
所述路段出行量全局关系图包括所述路段出行量矩阵、第二邻接关系矩阵和所述目标区域的边集;
所述第二邻接关系矩阵的任一元素bij为路段i的出行量向量与路段j的出行量向量相似度,任一行向量由任一路段i的出行量向量与自身的出行量向量及其他路段的出行量向量的相似度构成;
所述路段出行量OD次数关系图包括所述路段出行量矩阵、第三邻接关系矩阵和所述目标区域的边集;
所述第三邻接关系矩阵的任一元素cij为归一化后的路段i到路段j的OD次数,任一行向量由归一化后的路段i到自身的OD次数及归一化后的路段i到其他各路段的OD次数,
xij为目标区域最后一个历史时隙之前路段i与路段j间的OD次数或者xij为所有时间段内的路段i与路段j间OD次数的平均值,R为目标区域的所有路段集合;
所述外部因素特征矩阵由每个时间段的外部因素特征向量构成,任一时间段的外部因素特征向量由该时间段的外部因素构成;
(3)以所述路段出行量局部关系图、路段出行量全局关系图和路段出行量OD次数关系图和外部因素特征矩阵作为输入,以最后一个历史时隙内各路段车辆出行量为标签,对TCN和GCN混合深度网络进行训练,得到预测模型;
所述TCN和GCN混合深度网络由四个通道和卷积层构成,其中三个通道各自均包含依次设置的一层TCN、两层GCN和一层TCN,这三个通道内的一个通道的输入为所述路段出行量全局关系图、一个通道的输入为路段出行量局部关系图、另一个的输入为所述出行量OD次数关系图;第四个通道包括全连接层,该第四个通道的输入为外部因素特征矩阵;将所述四个通道的输出进行融合后,输入所述卷积层。
2.如权利要求1所述的车辆出行量预测模型构建方法,其特征在于,所述多个历史时隙的总时长为数周、数个月或数年。
3.如权利要求1所述的车辆出行量预测模型构建方法,其特征在于,所述时间段时长为数个小时、半天、一天或多天。
4.如权利要求1所述的车辆出行量预测模型构建方法,其特征在于,所述每个历史时隙为数秒或数分钟。
5.如权利要求1所述的车辆出行量预测模型构建方法,其特征在于,所述待预测时隙时长为数分钟、数小时、半天或一天。
6.如权利要求1所述的车辆出行量预测模型构建方法,其特征在于,对所述各时间段的外部因素中的各因素分别采用独热向量进行编码,之后对各因素的编码进行级联,得到总独热向量,该总独热向量为相应时间段的外部因素特征向量。
7.一种车辆出行量预测方法,对权利要求1-6任一权利要求所述待预测时隙的车辆出行量进行预测,其特征在于,包括
(一)数据采集与预处理,包括:
采集目标区域待预测时隙之前多个历史时隙内的车辆轨迹数据,并将车辆轨迹数据与目标区域的路段进行匹配,得到原始数据;
采集每时间段的外部因素获得外部因素数据,所述外部因素包括星期几、空气质量、温度和天气中一种或两种以上因素,所述多个历史时隙被划分为多个所述时间段;
(二)采用权利要求1-6任一权利要求所述步骤(2)方法利用多个历史时隙内的原始数据构建目标区域的路段出行量局部关系图、路段出行量全局关系图和路段出行量OD次数关系图;利用所述多个历史时隙内的外部因素数据构建外部因素特征矩阵;
(三)将步骤(二)构建的路段出行量局部关系图、路段出行量全局关系图、路段出行量OD次数关系图和外部因素特征矩阵输入权利要求1-6任一权利要求所述方法构建的模型中,输出待预测时隙内的车辆出行量预测值。
8.一种车辆出行量预测方法,对权利要求7所述待预测时隙的下一个待预测时隙内的车辆出行量进行预测,其特征在于,包括:
以权利要求7所述步骤(二)构建的路段出行量局部关系图、路段出行量全局关系图、路段出行量OD次数关系图和外部因素特征矩阵作为输入,以待预测时隙内的车辆出行量为输出对权利要求1-6任一权利要求所述方法构建的模型进行更新;
之后采用权利要求7所述步骤(一)、(二)所述方法及更新后的模型对下一个待预测时隙的车辆出行量进行预测。
9.一种车辆出行量预测系统,其特征在于,包括数据预处理模块、数据构建模块和预测及更新模块,所述数据预处理模块用于执行权利要求7所述步骤(一),所述数据构建模块用于执行权利要求7所述步骤(二),所述数据预测及更新模块用于对新待预测时隙的预测模型进行更新并利用更新后的模型对新待预测时隙的车辆出行量进行预测。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110590338.7A CN113380025B (zh) | 2021-05-28 | 2021-05-28 | 一种车辆出行量预测模型构建方法及预测方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110590338.7A CN113380025B (zh) | 2021-05-28 | 2021-05-28 | 一种车辆出行量预测模型构建方法及预测方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113380025A true CN113380025A (zh) | 2021-09-10 |
CN113380025B CN113380025B (zh) | 2022-08-30 |
Family
ID=77574817
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110590338.7A Active CN113380025B (zh) | 2021-05-28 | 2021-05-28 | 一种车辆出行量预测模型构建方法及预测方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113380025B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113514072A (zh) * | 2021-09-14 | 2021-10-19 | 自然资源部第三地理信息制图院 | 一种面向导航数据与大比例尺制图数据的道路匹配方法 |
CN114360254A (zh) * | 2022-03-17 | 2022-04-15 | 湖南工商大学 | 一种基于时空的交通流量预测方法及相关设备 |
CN115376317A (zh) * | 2022-08-22 | 2022-11-22 | 重庆邮电大学 | 一种基于动态图卷积和时序卷积网络的交通流预测方法 |
CN116502539A (zh) * | 2023-05-18 | 2023-07-28 | 淮阴工学院 | 一种VOCs气体浓度预测方法及系统 |
CN116612640A (zh) * | 2023-07-14 | 2023-08-18 | 之江实验室 | 一种车辆调度方法、装置、存储介质及电子设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107230349A (zh) * | 2017-05-23 | 2017-10-03 | 长安大学 | 一种在线式实时短时间交通流预测方法 |
CN110264709A (zh) * | 2019-05-06 | 2019-09-20 | 北京交通大学 | 基于图卷积网络的道路的交通流量的预测方法 |
CN111126680A (zh) * | 2019-12-11 | 2020-05-08 | 浙江大学 | 一种基于时间卷积神经网络的道路断面交通流量预测方法 |
CN111653088A (zh) * | 2020-04-21 | 2020-09-11 | 长安大学 | 一种车辆出行量预测模型构建方法及预测方法和系统 |
CN112419718A (zh) * | 2020-11-17 | 2021-02-26 | 东北大学秦皇岛分校 | 一种基于时空图卷积神经网络的交通拥塞传播预测方法 |
-
2021
- 2021-05-28 CN CN202110590338.7A patent/CN113380025B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107230349A (zh) * | 2017-05-23 | 2017-10-03 | 长安大学 | 一种在线式实时短时间交通流预测方法 |
CN110264709A (zh) * | 2019-05-06 | 2019-09-20 | 北京交通大学 | 基于图卷积网络的道路的交通流量的预测方法 |
CN111126680A (zh) * | 2019-12-11 | 2020-05-08 | 浙江大学 | 一种基于时间卷积神经网络的道路断面交通流量预测方法 |
CN111653088A (zh) * | 2020-04-21 | 2020-09-11 | 长安大学 | 一种车辆出行量预测模型构建方法及预测方法和系统 |
CN112419718A (zh) * | 2020-11-17 | 2021-02-26 | 东北大学秦皇岛分校 | 一种基于时空图卷积神经网络的交通拥塞传播预测方法 |
Non-Patent Citations (2)
Title |
---|
杨云等: "基于卷积网络与支持向量机的云资源预测模型", 《陕西科技大学学报》 * |
樊娜,等: "短时交通流预测模型", 《交通运输工程学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113514072A (zh) * | 2021-09-14 | 2021-10-19 | 自然资源部第三地理信息制图院 | 一种面向导航数据与大比例尺制图数据的道路匹配方法 |
CN114360254A (zh) * | 2022-03-17 | 2022-04-15 | 湖南工商大学 | 一种基于时空的交通流量预测方法及相关设备 |
CN115376317A (zh) * | 2022-08-22 | 2022-11-22 | 重庆邮电大学 | 一种基于动态图卷积和时序卷积网络的交通流预测方法 |
CN115376317B (zh) * | 2022-08-22 | 2023-08-11 | 重庆邮电大学 | 一种基于动态图卷积和时序卷积网络的交通流预测方法 |
CN116502539A (zh) * | 2023-05-18 | 2023-07-28 | 淮阴工学院 | 一种VOCs气体浓度预测方法及系统 |
CN116502539B (zh) * | 2023-05-18 | 2024-02-27 | 淮阴工学院 | 一种VOCs气体浓度预测方法及系统 |
CN116612640A (zh) * | 2023-07-14 | 2023-08-18 | 之江实验室 | 一种车辆调度方法、装置、存储介质及电子设备 |
CN116612640B (zh) * | 2023-07-14 | 2023-12-26 | 之江实验室 | 一种车辆调度方法、装置、存储介质及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN113380025B (zh) | 2022-08-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113380025B (zh) | 一种车辆出行量预测模型构建方法及预测方法和系统 | |
CN111653088B (zh) | 一种车辆出行量预测模型构建方法及预测方法和系统 | |
Diao et al. | A hybrid model for short-term traffic volume prediction in massive transportation systems | |
CN111223301B (zh) | 一种基于图注意力卷积网络的交通流量预测方法 | |
CN110570651A (zh) | 一种基于深度学习的路网交通态势预测方法及系统 | |
CN112561187B (zh) | 一种基于cnn-lstm的网约车目标订单预测方法 | |
CN111242292B (zh) | 基于深度时空网络的od数据预测方法及系统 | |
Han et al. | Short-term traffic prediction based on deepcluster in large-scale road networks | |
CN109376906B (zh) | 基于多维度轨迹的出行时间预测方法、系统及电子设备 | |
CN111242395B (zh) | 用于od数据的预测模型构建方法及装置 | |
CN115204478A (zh) | 一种结合城市兴趣点和时空因果关系的公共交通流量预测方法 | |
CN113780665B (zh) | 一种基于增强循环神经网络的私家车停留位置预测方法和系统 | |
Zhang et al. | Passenger flow forecast of rail station based on multi-source data and long short term memory network | |
CN111598333A (zh) | 客流数据的预测方法及装置 | |
Qiu et al. | Taxi origin-destination demand prediction with contextualized spatial-temporal network | |
CN115169714A (zh) | 城市地铁进出站客流量预测方法 | |
CN110490365B (zh) | 一种基于多源数据融合预测网约车订单量的方法 | |
CN116153125A (zh) | 一种基于ida-rnn模型的空间维度公交到达时间预测方法 | |
CN113192325B (zh) | 一种交通事件相关性分析方法及系统 | |
Xu et al. | A taxi dispatch system based on prediction of demand and destination | |
Wu et al. | A novel dynamically adjusted regressor chain for taxi demand prediction | |
Liu et al. | A MRT daily passenger flow prediction model with different combinations of influential factors | |
Rasaizadi et al. | Short‐Term Prediction of Traffic State for a Rural Road Applying Ensemble Learning Process | |
CN114117903A (zh) | 一种基于bp神经网络的轨道交通短时客流预测方法 | |
Rezazada et al. | What influences passenger’s arrival rate at stops in Melbourne: Wasea-Lstm, a novel deep learning model leveraging multi-source data fusion |
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 |