CN115206104A - 一种基于手机信令数据的城市居民交通出行方式识别方法 - Google Patents
一种基于手机信令数据的城市居民交通出行方式识别方法 Download PDFInfo
- Publication number
- CN115206104A CN115206104A CN202210878379.0A CN202210878379A CN115206104A CN 115206104 A CN115206104 A CN 115206104A CN 202210878379 A CN202210878379 A CN 202210878379A CN 115206104 A CN115206104 A CN 115206104A
- Authority
- CN
- China
- Prior art keywords
- travel
- uncertain
- moving
- membership
- mobile phone
- 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 47
- 230000011664 signaling Effects 0.000 title claims abstract description 39
- 230000006870 function Effects 0.000 claims abstract description 41
- 238000007781 pre-processing Methods 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 11
- 230000008569 process Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 4
- 238000011835 investigation Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000012850 discrimination method Methods 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 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
-
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/40—Business processes related to the transportation industry
-
- 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/0137—Measuring and analyzing of parameters relative to traffic conditions for specific applications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/20—Services signaling; Auxiliary data signalling, i.e. transmitting data via a non-traffic channel
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Analytical Chemistry (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Chemical & Material Sciences (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Signal Processing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明公开了一种基于手机信令数据的城市居民交通出行方式识别方法,包括:对获取的手机信令数据进行预处理;对手机信令数据进行移动段和停留段的划分,构建出行行程链;构建各种交通出行方式的隶属度函数;计算出每种出行方式的综合隶属度得分,并通过综合隶属度得分阈值,将获得的移动段进一步划分为确定移动段和不确定移动段;分别对确定移动段和不确定移动段的出行方式进行识别。本发明充分考虑了不同交通方式各自的特征和差异构建了适合手机信令数据的交通出行方式识别的隶属度函数,并结合了路网、导航等多源数据,制定了不确定移动段出行方式的识别规则,实现了四种交通方式的识别,有效提高了基于手机信令数据的个体交通方式识别精度。
Description
技术领域
本发明属于交通出行方式识别领域,具体涉及一种基于手机信令数据的城市居民交通出行方式识别方法。
背景技术
通过对居民个体出行方式的调查研究,掌握不同时间、地点的居民出行方式及出行时间等特征信息,一方面可以直接指导城市交通规划,辅助交通决策,有利于制定合理的交通规划方案,优化城市路网结构和交通设施布局。另一方面,有助于交通管理部门合理引导居民选择最佳的出行方式,有效调控城市车流、人流,减少交通拥堵问题的发生。
传统的居民出行方式信息主要通过入户调查(Face to Face)、问卷调查以及计算机辅助电话调查(CATI)等方式获取。但是这些方式受被调查人个人主观认知的影响很大,问卷回收率和收集的数据质量都比较低,同时,需要耗费大量的人力物力和财力。
随着计算机技术、信息与通信技术等现代技术的快速发展,通过手机信令数据来进行个体出行方式的识别成为可能。手机信令数据采集过程简单、便于大规模获取、生成规则统一、实时性好,在人群活动特征分析、居民出行行为分析以及个体出行交通方式识别等方面有及其广泛的应用。
已有的利用手机信令数据来进行城市居民交通出行方式识别的技术,主要可分为基于规则的判别方法、基于机器学习的方法。其中,基于规则的判别方法,例如公开号为CN106197458B、CN105117789A等专利,根据专家的先验知识指定识别规则,识别容易受研究者个人主观因素影响,通常准确率较低。基于机器学习的方法,如公开号为CN112511971A、CN111653096等专利,需要大量的标记样本用于模型的训练,往往存在标记样本获取难度大,模型精度易受样本质量的影响,小数据集模型训练精度不高且易出现过拟合现象等问题。
发明内容
发明目的:针对当前基于手机信令数据的城市居民出现方式识别方法过于简单、没有结合多源数据、识别精度不高等方面的不足,提供一种基于手机信令数据的城市居民交通出行方式识别方法,其充分考虑了不同交通方式各自的特征和差异构建了适合手机信令数据的交通出行方式识别的隶属度函数,并结合了路网、导航等多源数据,制定了不确定移动段出行方式的识别规则,实现了对步行、自行车、公交、小汽车四种交通方式的识别,有效提高了基于手机信令数据的个体交通方式识别精度。
技术方案:为实现上述目的,本发明提供一种基于手机信令数据的城市居民交通出行方式识别方法,包括如下步骤:
S1:对获取的手机信令数据进行预处理,剔除数据中的漂移点和重复记录;
S2:通过SMUoT算法对预处理后的手机信令数据进行移动段和停留段的划分,构建出行行程链;
S3:构建各种交通出行方式的隶属度函数;
S4:根据隶属度函数计算出每种出行方式的综合隶属度得分,并通过综合隶属度得分阈值,将步骤S2获得的移动段进一步划分为确定移动段和不确定移动段;
综合隶属度得分阈值用于将移动段划分为确定移动段和不确定移动段,具体数值由多次试验得出,在0.7时能够获得较好的识别精度;
S5:分别对确定移动段和不确定移动段的出行方式进行识别。
进一步地,所述步骤S1中预处理的过程为:
A1:据获取的手机信令数据的经纬度、时间戳字段,计算相邻轨迹点之间的速度,将速度大于设定值的漂移点剔除;
A2:将手机信令数据中同一时间的重复、错误记录进行剔除。
进一步地,所述步骤S2中出行行程链的构建方法为:
在对居民个体的出行轨迹进行停留段和移动段的识别后,将表示停留状态的连续轨迹点组合成停留段Sk={p1,p2,p3,...,pi},k表示出行链中的第k个停留段,i表示停留段中轨迹点数量,0≤i≤M;
将表示移动状态的连续轨迹点组合成移动段Ml={p1,p2,...,pj},l表示出行链中的第l个移动段,j表示移动段中轨迹点的数量,0≤j≤M;
移动段和停留段交替构成居民个体的出行链T={S1,M1,S2,M2…,Sn},n表示第n个停留段或者移动段。
进一步地,所述步骤S3中交通出行方式包括步行、自行车、公交、小汽车四种交通方式,所述隶属度函数包括出行距离隶属度函数、出行时长隶属度函数、出行速度隶属度函数。
进一步地,所述步骤S4中综合隶属度得分的计算方式为:将计算的每个移动段的移动距离、移动时长、平均移动速度作为自变量分别输入到出行距离隶属度函数、出行时长隶属度函数、出行速度隶属度函数,得到每种出行方式的出行距离隶属度、出行时长隶属度和出行平均速度隶属度,并将三者相乘得到综合隶属度得分。
进一步地,所述步骤S4中确定移动段和不确定移动段的划分方式为:将综合隶属度得分小于综合隶属度得分阈值的出行段划分为不确定移动段,将综合隶属度得分超过综合隶属度得分阈值的出行段划分为确定移动段,并将对应的出行方式作为该确定移动段的出行方式。
进一步地,所述步骤S5中确定移动段的出行方式识别方法为:将综合隶属度得分最高的出行方式作为最终识别的出行方式。
进一步地,所述步骤S5中通过停车频率约束、高速轨迹点占比约束、平均方位角变化约束、时间范围约束、导航时间相似度约束这五个约束条件来实现对不确定移动段出行方式的识别,具体为:
B1:根据停车频率约束识别不确定移动段:计算不确定移动段的停车频率,将停车频率大于规定阈值的识别为公交出行;
B2:根据高速轨迹点占比约束识别不确定移动段:计算不确定移动段中高速轨迹点占比,将占比大的不确定移动段识别为小汽车出行;
B3:根据平均方位角变化约束识别不确定移动段:将平均方位角变化大于阈值的不确定移动段出行方式识别为步行;
B4:根据时间范围约束识别不确定移动段:将晚间及凌晨期间不确定移动段公交出行概率设置为零,即排除居民晚间及凌晨的公交出行;
B5:根据导航时间相似度约束识别不确定移动段:对于使用前四种约束规则都无法识别的不确定移动段,则计算其导航时间相似度以此来判断其出行方式,使用高德地图导航平台API接口,将对应的移动轨迹数据输入可获取平台规划的公交、小汽车、步行、自行车等出行方式的预期花费时间,根据获取的导航预期花费时间与实际花费时间计算导航时间相似度,将导航时间相似度最大的交通方式作为对应不确定移动段的出行方式。
有益效果:本发明与现有技术相比,充分考虑了不同交通方式各自的特征和差异构建了适合手机信令数据的交通出行方式识别的隶属度函数,并结合了路网、导航等多源数据,制定了不确定移动段出行方式的识别规则,实现了对步行、自行车、公交、小汽车四种交通方式的识别,有效提高了基于手机信令数据的个体交通方式识别精度。
附图说明
图1为本发明的基本原理图;
图2为本发明的工作流程图;
图3为本发明中居民个体出行行程链示意图;
图4为本发明中应用隶属度函数进行确定移动段和不确定移动段划分过程示意图;
图5为本发明中的移动段轨迹点方位角示意图。
具体实施方式
下面结合附图和具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
本发明提供一种基于手机信令数据的城市居民交通出行方式识别方法,如图1和图2所示,包括如下步骤:
S1:对获取的手机信令数据进行预处理:
A1:根据获取的某一城市居民个体的出行轨迹,主要包含用户ID、基站经纬度、时间戳等字段,计算相邻轨迹点之间的速度,将速度大于120km/h的漂移点剔除;
A2:将手机信令数据中同一时间的重复、字段不全等错误记录进行剔除。
本发明中居民个体的一次出行轨迹Traj={p1,p2,p3,...,pm},1≤m≤M(M表示轨迹点总数量,Pm表示第m个轨迹点)由一系列手机定位数据记录点p构成,p=<UserID,Time,Longitude,Latitude>为一个四元组,其中UserID表示居民个体编号,Time表示居民个体被定位的时间,Longitude和Latitude表示定位时用户手机连接的基站经纬度。
S2:通过SMUoT算法对预处理后的手机信令数据进行移动段Ml和停留段Sk的划分,构建出行行程链,并且提取出移动段Ml:
居民个体在出行过程中包括两种状态,停留和移动,在对居民个体的出行轨迹进行停留段和移动段的识别后,将表示停留状态的连续轨迹点组合成停留段Sk={p1,p2,p3,...,pi},k表示出行链中的第k个停留段,i表示停留段中轨迹点数量,0≤i≤M;
将表示移动状态的连续轨迹点组合成移动段Ml={p1,p2,...,pj},l表示出行链中的第l个移动段,j表示移动段中轨迹点的数量,0≤j≤M;
移动段和停留段交替构成居民个体的出行行程链T={S1,M1,S2,M2…,Sn},n表示第n个停留段或者移动段,出行行程链T具体如图3所示。
S3:构建各种交通出行方式的隶属度函数:
本发明中交通出行方式包括步行、自行车、公交、小汽车四种交通方式,具体如下表1~表3:
表1:不同出行方式的常见出行距离
表2:不同出行方式的常见出行时长
表3:不同出行方式的常见出行速度
根据上述的关于不同出行方式的出行距离、出行时长、出行速度的先验知识分别构建与步行、自行车、公交、小汽车对应的出行距离隶属度函数、出行时长隶属度函数、出行速度隶属度函数,具体如下表4:
表4:不同出行方式各属性对应的隶属度函数
S4:根据隶属度函数计算出每种出行方式的综合隶属度得分,并通过综合隶属度得分阈值,将步骤S2获得的移动段进一步划分为确定移动段和不确定移动段:
确定移动段和不确定移动段的具体解释如下:
本发明通过手机信令数据识别居民交通出行方式,因此只需要考虑移动段的居民出行方式;通过隶属度函数的计算将最终隶属度得分高于阈值δ的出行段划分为出行方式确定移动段Mcertain,简称确定移动段,确定移动段所对应的交通出行方式为TMcertain;将隶属度得分低于δ的出行段划分为出行方式不确定移动段Muncertain,简称不确定移动段,不确定移动段所对应的交通出行方式为TMunertain。不确定移动段在本发明中属于一种过度性轨迹段,仅借助本发明构建的出行时长隶属度函数、出行距离隶属度函数、出行速度隶属度函数并不能直接识别其交通出行方式。
综合隶属度得分的计算方式为:将计算的每个移动段的移动距离、移动时长、平均移动速度作为自变量分别输入到出行距离隶属度函数、出行时长隶属度函数、出行速度隶属度函数,得到每种出行方式的出行距离隶属度、出行时长隶属度和出行平均速度隶属度,并将三者相乘得到综合隶属度得分。
参照图4,确定移动段和不确定移动段的划分过程为:首先,获取某一移动段,根据表4中的隶属度函数,计算移动段每种出行方式对应的综合隶属度得分Uk,接着根据Uk得到最高隶属度得分Umax,最后将Umax超过阈值δ=0.7的移动段判定为确定移动段Mcertain,反之为不确定移动段Muncertain。
Umax和Uk计算公式如下:
Umax=max(Uk),k∈{walk,bike,bus,car} (1)
Uk=Tk*Dk*Vk,k∈{walk,bike,bus,car} (2)
式中,Tk,Dk,Vk分别表示交通方式为k时所计算的时间隶属度得分,距离隶属度得分,速度隶属度得分,具体公式可参考表4。
S5:对确定移动段Mcertain出行方式识别;令确定移动段的出行方式TMcertain=k,即将确定移动段最高隶属度得分Umax所对应的交通方式k作为确定移动段的交通出行方式。
对不确定移动段Muncertain的出行方式识别:
不确定移动段和确定移动段将构成多种用户换乘方式组合,比如小汽车-不确定移动段-小汽车、公交-不确定移动段、自行车-不确定移动段-自行车等,本发明结合用户出行方式选择规律和出行方式换乘特征提出如下规则来辅助不确定移动段Muncertain出行方式的识别:
规则1:停车频率约束。
公交出行与其他出行一个非常显著的区别就是,公交出行由于存在大量的公交站点,存在大量的乘客上下车行为,那么此时公交出行的停车频率会相交其他出行方式有较大增长,因此本发明将停车频率大于0.8的出行方式不确定移动段的出行方式TMuncertain识别为公交出行,即令TMuncertain=bus;反之,则该不确定移动段待下一步处理。本发明将两个轨迹点之间速度小于等于5km/h且距离小于60m的轨迹段算作一次停留。对于,不确定移动段Muncertain而言,停车频率计算公式如下:
式中:ns表示不确定轨迹段的停车次数,Nb表示不确定轨迹段所经过的公交站总数。
规则2:高速轨迹点占比约束。
城市小汽车行驶的最高速度,明显高于步行、自行车和公交车,因此,本发明将瞬时速度大于15m/s的轨迹点作为高速轨迹点,计算不确定移动段中的高速轨迹点的占比。将高速轨迹点占比超过20%的不确定移动段的出行方式TMuncertain识别为小汽车,即令TMuncertain=car;反之,则该不确定移动段待下一步处理。高速轨迹点占比Rh计算公式如下:
式中,nh表示不确定移动段中高速轨迹点的数量;Nt表示不确定移动段中轨迹点的总数量。
规则3:平均方位角变化约束。
居民步行时可以根据自己的喜好任意的改变步行方向,而自行车、公交和小汽车受到城市路网结构的约束,只能在交叉口等特定区域或有需要的时候才会变换前进方向;所以,步行的轨迹点方位角的变化相对于其他三种出行方式更大。本发明,将平均方向变化Davg大于90°的不确定移动段的出行方式TMuncertain识别为步行,即令TMuncertain=walk;反之,则该不确定移动段待下一步处理。方位角示意图见图5,平均方向变化Davg计算方法如下:
规则4:时间范围约束。
凌晨和晚上公交基本停运,因此,将晚上23点-第二天凌晨6点之间的不确定移动段的公交出行概率设置为零,即排除居民晚间及凌晨的公交出行。
规则5:导航时间相似度约束。
对于使用前四种约束规则都无法识别的不确定移动段,本发明结合导航数据,使用导航时间相似度来识别其最终的出行方式。使用高德地图导航平台API接口,将对应的不确定移动段轨迹数据输入,可获取平台规划的公交、小汽车、步行、自行车等出行方式的预期花费时间,由此计算每种交通方式的导航时间隶属度Rk,navTime,k∈{walk,bike,bus,car},将最大导航时间相似度Rmax=max(Rk,navTime),k∈{walk,bike,bus,car}所对应的出行方式k设置为该不确定移动段的出行方式,即令TMuncertain=k;注意:在计算导航时间相似度过程中,若该不确定移动段符合规则4的时间范围约束,则首先应当令公交对应的导航时间隶属度为零,即令Rbus,navTime=0,然后再计算最大导航时间相似度。时间隶属度计算公式如下:
式中,Treal表示移动段移动的总时间,Tnav表示导航平台返回的不同交通方式的预期花费时间。
S6:至此,已完成某一居民个体出行轨迹中某一个移动段出行方式的识别,重复上述步骤S4到S6完成其他移动段出行方式的识别,最终完成该居民个体出行轨迹中出行方式的识别。
基于上述方案,为了验证本发明方法的有效性,本实施例中将上述方案进行实例应用,具体过程如下:
(1)数据预处理;获取某一居民个体出行轨迹数据Traj,将速度大于120km/h的记录以及同一时间的重复、字段不全等错误记录进行剔除,得到如表5所示的居民个体出行轨迹数据。
表5:预处理后的居民个体出行轨迹数据
(2)行程链构建;使用SMUoT算法将经预处理后的居民个体轨迹数据划分为停留段Si和移动段Mi,并提取其中的移动段Mi,待后续处理,本实施例中该个体轨迹数据共提取6个移动段。如表6所示为提取的移动段。
表6:提取的移动段
(3)综合隶属度计算及确定移动段和不确定移动段划分:
根据移动段中轨迹点的位置和时间戳计算每个移动段的出行时长t(minute)、出行距离d(km)和平均出行速度u(km/h)。将计算出来的出行时长t、出行距离d和出行速度v作为表1.4中隶属度函数的自变量,计算步行、公交、自行车、小汽车四种交通各自的综合隶属度和最大隶属度,将最大隶属度超过阈值δ=0.7的移动段判定为确定移动段Mcertain,反之为不确定移动段Muncertain,划分结果如表7所示。
表7确定移动段和不确定移动段划分
(4)确定移动段Mcertain出行方式识别:
将确定移动段最高隶属度得分Umax所对应的交通方式k作为确定移动段的交通出行方式,如移动段M1中Umax=Ubus,则该移动段的出行方式TMcertain=bus,同理,识别其他确定移动段M2,M3,M6的出行方式依次为:walk、ride、walk。这样依次完成所有确定移动段出行方式的识别。
(5)不确定移动段Muncertain出行方式识别:
不确定移动段的出行方式需依次根据停车频率、高速轨迹点占比、平均方位角变化、时间范围和导航时间相似度这五个约束条件来进行识别。如针对不确定移动段M4首先计算其停车频率不符合停车频率约束;接着计算M4中高速轨迹点占比Rh=0.15,不符合高速轨迹点占比约束;然后计算平均方位角变化Davg=70,不符合平均方位角变化约束;最后将该移动段轨迹输入导航平台API接口获取平台规划的公交、小汽车、步行、自行车等出行方式的预期花费时间,计算不同交通方式的导航时间相似度Rbus,navTime=0.23,Rcar,navTime=0.45,Rwalk,navTime=0.85,Rride,navTime=0.56(如果不确定移动段符合时间范围约束,则应该将Rbus,navTime修正为0,这里M4不符合无需修正),其中Rmax=Rwalk,navTime=0.85,则不确定移动段M4所对应的出行方式为walk。同理,识别不确定移动段M5所对应的出行方式为bus。至此,已完成该居民个体的交通出行方式识别。
为了验证本发明方法的实际效果,本实例将本发明识别方法和现有仅通过距离、时耗、速度三个属性值的映射函数进行识别的方法,如公开号为CN 102136192的专利中提到的方法,进行对比,分别计算两种方法的精确率、召回率。准确率和召回率的计算公式如下:
式中,ACC,REC,Ncorrect,Ndistinguish,Nlabel分别表示:精确率,召回率,算法识别的出行方式与实际的出行方式一致的移动段的数量,算法识别各出行方式的移动段数量,实际的各出行方式的移动段数量。经过计算后可以得到本发明识别方法及公开号为CN102136192的专利中提到的仅凭映射函数的方法对不同交通方式的识别准确率和召回率,计算结果如下:
表8本发明识别方法的识别结果混淆矩阵
表9仅凭映射函数进行识别的结果混淆矩阵
从表8和表9可见,本发明识别方法与现有方法相比,精确率和召回率都有了大幅提升,从而验证了本发明方法的效果。
Claims (8)
1.一种基于手机信令数据的城市居民交通出行方式识别方法,其特征在于,包括如下步骤:
S1:对获取的手机信令数据进行预处理,剔除数据中的漂移点和重复记录;
S2:通过SMUoT算法对预处理后的手机信令数据进行移动段和停留段的划分,构建出行行程链;
S3:构建各种交通出行方式的隶属度函数;
S4:根据隶属度函数计算出每种出行方式的综合隶属度得分,并通过综合隶属度得分阈值,将步骤S2获得的移动段进一步划分为确定移动段和不确定移动段;
S5:分别对确定移动段和不确定移动段的出行方式进行识别。
2.根据权利要求1所述的一种基于手机信令数据的城市居民交通出行方式识别方法,其特征在于,所述步骤S1中预处理的过程为:
A1:据获取的手机信令数据的经纬度、时间戳字段,计算相邻轨迹点之间的速度,将速度大于设定值的漂移点剔除;
A2:将手机信令数据中同一时间的重复、错误记录进行剔除。
3.根据权利要求1所述的一种基于手机信令数据的城市居民交通出行方式识别方法,其特征在于,所述步骤S2中出行行程链的构建方法为:
在对居民个体的出行轨迹进行停留段和移动段的识别后,将表示停留状态的连续轨迹点组合成停留段Sk={p1,p2,p3,…,pi},k表示出行链中的第k个停留段,i表示停留段中轨迹点数量,0≤i≤M;
将表示移动状态的连续轨迹点组合成移动段Ml={p1,p2,…,pj},l表示出行链中的第l个移动段,j表示移动段中轨迹点的数量,0≤j≤M;
移动段和停留段交替构成居民个体的出行链T={S1,M1,S2,M2…,Sn},n表示第n个停留段或者移动段。
4.根据权利要求1所述的一种基于手机信令数据的城市居民交通出行方式识别方法,其特征在于,所述步骤S3中交通出行方式包括步行、自行车、公交、小汽车四种交通方式,所述隶属度函数包括出行距离隶属度函数、出行时长隶属度函数、出行速度隶属度函数。
5.根据权利要求4所述的一种基于手机信令数据的城市居民交通出行方式识别方法,其特征在于,所述步骤S4中综合隶属度得分的计算方式为:将计算的每个移动段的移动距离、移动时长、平均移动速度作为自变量分别输入到出行距离隶属度函数、出行时长隶属度函数、出行速度隶属度函数,得到每种出行方式的出行距离隶属度、出行时长隶属度和出行平均速度隶属度,并将三者相乘得到综合隶属度得分。
6.根据权利要求1所述的一种基于手机信令数据的城市居民交通出行方式识别方法,其特征在于,所述步骤S4中确定移动段和不确定移动段的划分方式为:将综合隶属度得分小于综合隶属度得分阈值的出行段划分为不确定移动段,将综合隶属度得分超过综合隶属度得分阈值的出行段划分为确定移动段,并将对应的出行方式作为该确定移动段的出行方式。
7.根据权利要求1所述的一种基于手机信令数据的城市居民交通出行方式识别方法,其特征在于,所述步骤S5中确定移动段的出行方式识别方法为:将综合隶属度得分最高的出行方式作为最终识别的出行方式。
8.根据权利要求1所述的一种基于手机信令数据的城市居民交通出行方式识别方法,其特征在于,所述步骤S5中通过停车频率约束、高速轨迹点占比约束、平均方位角变化约束、时间范围约束、导航时间相似度约束这五个约束条件来实现对不确定移动段出行方式的识别,具体为:
B1:根据停车频率约束识别不确定移动段:计算不确定移动段的停车频率,将停车频率大于规定阈值的识别为公交出行;
B2:根据高速轨迹点占比约束识别不确定移动段:计算不确定移动段中高速轨迹点占比,将占比大的不确定移动段识别为小汽车出行;
B3:根据平均方位角变化约束识别不确定移动段:将平均方位角变化大于阈值的不确定移动段出行方式识别为步行;
B4:根据时间范围约束识别不确定移动段:将晚间及凌晨期间不确定移动段公交出行概率设置为零,即排除居民晚间及凌晨的公交出行;
B5:根据导航时间相似度约束识别不确定移动段:对于使用前四种约束规则都无法识别的不确定移动段,则计算其导航时间相似度以此来判断其出行方式,使用高德地图导航平台API接口,将对应的移动轨迹数据输入可获取平台规划的出行方式的预期花费时间,根据获取的导航预期花费时间与实际花费时间计算导航时间相似度,将导航时间相似度最大的交通方式作为对应不确定移动段的出行方式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210878379.0A CN115206104B (zh) | 2022-07-25 | 2022-07-25 | 一种基于手机信令数据的城市居民交通出行方式识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210878379.0A CN115206104B (zh) | 2022-07-25 | 2022-07-25 | 一种基于手机信令数据的城市居民交通出行方式识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115206104A true CN115206104A (zh) | 2022-10-18 |
CN115206104B CN115206104B (zh) | 2024-04-23 |
Family
ID=83583076
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210878379.0A Active CN115206104B (zh) | 2022-07-25 | 2022-07-25 | 一种基于手机信令数据的城市居民交通出行方式识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115206104B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116822779A (zh) * | 2023-02-06 | 2023-09-29 | 长安大学 | 一种基于手机信令数据的高速公路机动车碳排放计算方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105117789A (zh) * | 2015-07-29 | 2015-12-02 | 西南交通大学 | 基于手机信令数据的居民出行方式综合判别的方法 |
CN111653094A (zh) * | 2020-05-29 | 2020-09-11 | 南京瑞栖智能交通技术产业研究院有限公司 | 一种基于手机信令数据且含路网修正的城市出行方式综合识别方法 |
CN111653097A (zh) * | 2020-05-29 | 2020-09-11 | 南京瑞栖智能交通技术产业研究院有限公司 | 一种基于手机信令数据且含个人属性修正的城市出行方式综合识别方法 |
CN111653096A (zh) * | 2020-05-29 | 2020-09-11 | 南京瑞栖智能交通技术产业研究院有限公司 | 一种基于手机信令数据的城市出行方式识别方法 |
CN112542045A (zh) * | 2020-12-01 | 2021-03-23 | 江苏欣网视讯软件技术有限公司 | 基于手机信令识别交通出行方式的方法与*** |
CN113573238A (zh) * | 2021-06-11 | 2021-10-29 | 北京交通大学 | 一种基于手机信令的联程旅客出行链识别方法 |
-
2022
- 2022-07-25 CN CN202210878379.0A patent/CN115206104B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105117789A (zh) * | 2015-07-29 | 2015-12-02 | 西南交通大学 | 基于手机信令数据的居民出行方式综合判别的方法 |
CN111653094A (zh) * | 2020-05-29 | 2020-09-11 | 南京瑞栖智能交通技术产业研究院有限公司 | 一种基于手机信令数据且含路网修正的城市出行方式综合识别方法 |
CN111653097A (zh) * | 2020-05-29 | 2020-09-11 | 南京瑞栖智能交通技术产业研究院有限公司 | 一种基于手机信令数据且含个人属性修正的城市出行方式综合识别方法 |
CN111653096A (zh) * | 2020-05-29 | 2020-09-11 | 南京瑞栖智能交通技术产业研究院有限公司 | 一种基于手机信令数据的城市出行方式识别方法 |
WO2021237812A1 (zh) * | 2020-05-29 | 2021-12-02 | 南京瑞栖智能交通技术产业研究院有限公司 | 一种基于手机信令数据且含个人属性修正的城市出行方式综合识别方法 |
CN112542045A (zh) * | 2020-12-01 | 2021-03-23 | 江苏欣网视讯软件技术有限公司 | 基于手机信令识别交通出行方式的方法与*** |
CN113573238A (zh) * | 2021-06-11 | 2021-10-29 | 北京交通大学 | 一种基于手机信令的联程旅客出行链识别方法 |
Non-Patent Citations (1)
Title |
---|
杜亚朋;雒江涛;程克非;唐刚;徐正;罗克韧;余疆;: "基于手机信令和导航数据的出行方式识别方法", 计算机应用研究, no. 08, 21 July 2017 (2017-07-21) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116822779A (zh) * | 2023-02-06 | 2023-09-29 | 长安大学 | 一种基于手机信令数据的高速公路机动车碳排放计算方法 |
CN116822779B (zh) * | 2023-02-06 | 2024-06-04 | 长安大学 | 一种基于手机信令数据的高速公路机动车碳排放计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115206104B (zh) | 2024-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109448370B (zh) | 一种基于车辆轨迹数据的交通控制子区划分方法 | |
CN110298500B (zh) | 一种基于出租车数据和城市路网的城市交通轨迹数据集生成方法 | |
CN110836675B (zh) | 一种基于决策树的自动驾驶搜索决策方法 | |
CN105957342A (zh) | 基于众包时空大数据的车道级道路测图方法及*** | |
CN107330469B (zh) | 一种基于gps轨迹数据的出行方式识别方法 | |
CN112530158B (zh) | 一种基于历史轨迹的路网补充方法 | |
CN106446960A (zh) | 一种密度聚类划分时间段的动态热度路网构建方法 | |
CN109886474B (zh) | 一种面向自动驾驶车辆测试的封闭测试场规划方法 | |
CN113378891A (zh) | 基于轨迹分布表示的城市区域关系可视分析方法 | |
CN106651728A (zh) | 一种综合运输体系客运方式优势运距的确定方法 | |
CN112309118A (zh) | 一种基于时空相似度的车辆轨迹推算方法 | |
CN112884014A (zh) | 一种基于路段拓扑结构分类的交通速度短时预测方法 | |
CN115206104A (zh) | 一种基于手机信令数据的城市居民交通出行方式识别方法 | |
CN109121094A (zh) | 伪码信令数据预处理与出行链识别方法 | |
CN113205213B (zh) | 一种基于互联网地图数据的公交评价方法及*** | |
Huang et al. | Mining crowdsourced trajectory and geo‐tagged data for spatial‐semantic road map construction | |
CN110674226A (zh) | 一种基于轨迹数据的远距离交通节点关联性挖掘方法 | |
CN112801482B (zh) | 一种数据驱动的非机动车骑行适宜性评价*** | |
CN111931968B (zh) | 一种应用手机gps数据来优化公共自行车站点布局的方法 | |
CN113011651A (zh) | 一种以提高骑行质量为目标的城市骑行路径规划方法 | |
CN113723761A (zh) | 一种基于运行数据的多维度城市公交运行服务可靠度评价方法 | |
CN106875688B (zh) | 一种应用gps数据识别公交车和小汽车的方法 | |
Li et al. | A GPS Trajectory Segmentation Method for Transportation Mode Inference | |
CN114973668B (zh) | 一种基于拓扑步数分析的城市道路交通薄弱环节识别方法 | |
CN114896523B (zh) | 一种基于乡村旅游线路的道路规划方法及装置 |
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 |