CN107545730A - 一种基于公交ic卡数据的站点上下车乘客数估计方法 - Google Patents
一种基于公交ic卡数据的站点上下车乘客数估计方法 Download PDFInfo
- Publication number
- CN107545730A CN107545730A CN201710806799.7A CN201710806799A CN107545730A CN 107545730 A CN107545730 A CN 107545730A CN 201710806799 A CN201710806799 A CN 201710806799A CN 107545730 A CN107545730 A CN 107545730A
- Authority
- CN
- China
- Prior art keywords
- bus
- website
- card
- passenger
- msub
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Traffic Control Systems (AREA)
Abstract
一种基于公交IC卡数据的站点上下车乘客数估计方法,本发明涉及站点上下车乘客数估计方法。本发明的目的是为了解决现有公交出行链法推导站点下车乘客数需要的工作量大,而站点吸引权重法推导站点下车乘客数准确度低的问题。具体过程为:一:根据乘客上车时的刷卡记录时间,求解公交线路的刷卡时间阈值;二:确定公交车未知对应站点时的每个站点上车人数;三:根据公交车未知对应站点时的每个站点上车人数、IC卡信息、公交调度信息中的每一站到站时间、公交线路编号以及公交车辆编号,确定公交车每个对应站点的上车人数;四:根据三确定公交车每个对应站点的下车人数。本发明用于站点上下车乘客数估计领域。
Description
技术领域
本发明涉及站点上下车乘客数估计方法。
背景技术
目前,大部分的公交政策制定的数据来源都是来自人工数据采集的方式,但是该方式 需耗费大量的人力、财力。因此开始有学者将目光聚焦于IC卡数据,试图通过IC卡信息 来获取乘客的出行数据,但是通过该方式获取的下车乘客人数,目前常用的公交出行链以 及站点吸引权重法存在着明显的缺陷,公交出行链法估计站点下车乘客数需要的工作量太 大,而站点吸引权重法估计站点下车乘客数准确度低。
发明内容
本发明的目的是为了解决现有公交出行链法估计站点下车乘客数需要的工作量大,而 站点吸引权重法估计站点下车乘客数准确度低的问题,而提出一种基于公交IC卡数据的 乘客上下车站点估计方法。
一种基于公交IC卡数据的站点上下车乘客数估计方法具体过程为:
步骤一:根据乘客上车时的刷卡记录时间,求解公交线路的刷卡时间阈值tmin;
步骤二:根据公交线路的刷卡时间阈值tmin,确定公交车未知对应站点时的每个站点 上车人数;
步骤三:根据公交车未知对应站点时的每个站点上车人数、IC卡信息、公交调度信息中的每一站到站时间、公交线路编号以及公交车辆编号,确定公交车每个对应站点的上车人数;
步骤四:根据步骤三确定公交车每个对应站点的下车人数。
本发明的有益效果为:
本发明提供了一种基于乘客公交IC卡刷卡时间和公交调度信息,推导乘客上车站点 人数的方法,将公交出行链与站点吸引权重法相结合,通过建立假设条件和相关集合,分 析乘客下一次上车站点满足集合的情况,来推导乘客的下车站点人数,若无法构成闭合的 公交出行链,其高频上车站点,作为本次可参考的下车站点。而若是公交出行链既不闭合, 而又无法获取乘客的高频上车站点,则通过站点的吸引权重确定乘客的下车站点;降低了 推导站点下车乘客数的工作量,提高了准确度,解决了现有公交出行链法推导站点下车乘 客数需要的工作量大,而站点吸引权重法推导站点下车乘客数准确度低的问题。
结合实施例一本发明案例优化10:00~11:00时段37路东山口站发车时刻,因此需分析 10:00~11:00时段从东山口站往赤沙总站的客流需求,从数据库中提取出该时段此运行方向 的公交车辆刷卡记录,其中在该时段共有端号为02560303、09500496、03760206的三趟公 交车辆,发车时刻分别为10:20、10:40、11:00,其中截取端号为02560303的公交车辆部分 刷卡数据,如表3所示,以此为例阐述时间阈值的推导过程,将该部分数据导入上车站点 推导算法的程序中,得到刷卡时间阈值为89s,并且与表1中的公交站点信息相匹配,得到 各站点上车乘客数。
而对于乘客选择何站下车则分下面两种情况进行讨论。
(1)第一种情况
当乘客下一次上车站点位于本线路下游站点或位于下游站点附近站点时,则下一次上 车站点为本次下车站点,而在本案例中02560303端号车辆10:00~11:00的刷卡记录符合该 情况的占42%,以表3中乘客的刷卡记录为例,通过该方法求得乘客下车站点部分刷卡记 录如表4所示。
(2)第二种情况
当不满足第一种情况时,寻找乘客的近日多次上车站点,其作为参考站点的发生时间 为2015.01.02~2015.01.07,其中频次大于等于2次,根据站点近日上车次数求算出乘客在 各站点的下车概率,02560303端号车辆10:00~11:00的刷卡记录符合该情况的占22%,以 表3中乘客的刷卡记录为例,得到部分乘客近日多次上车站点集,见表5,其中/X表示该 站点的上车次数。
而剩余的乘客,其上车站点则根据高站点的客流量作为权重进行分配,高客流量站点 分别为五羊新村、客村立交、珠影、大塘,下车吸引率根据其对应上车人数得到,分别为 0.252、0.260、0.244、0.244。因此,可根据高客流站点的下车吸引率计算得到所对应的下 车人数。
最后得到在10:00~11:00期间,端号为02560303的公交车辆在各公交站点的上下车乘 客数,通过实地调查,各站点刷卡率保持在77%~81%,与城市总体80%的刷卡率接近,忽 略各站点自身刷卡率的浮动情况,将各公交站点按城市总体刷卡率80%的比例换算成站点 实际上下车人数。
附图说明
图1为本发明流程图;
图2为本发明刷卡时间差临界值推导流程图;
图3为本发明确定公交车未知对应站点时的每个站点上车人数流程图;
图4为本发明闭合的公交出行链示意图;
图5为本发明微闭合的公交出行链示意图;
图6为本发明断裂的公交出行链示意图。
具体实施方式
具体实施方式一:结合图1说明本实施方式,本实施方式的一种基于公交IC卡数据的站点上下车乘客数估计方法具体过程为:
步骤一:根据乘客上车时的刷卡记录时间,求解公交线路的刷卡时间阈值tmin;
步骤二:根据公交线路的刷卡时间阈值tmin,确定公交车未知对应站点时的每个站点 上车人数;
步骤三:根据公交车未知对应站点时的每个站点上车人数、IC卡信息、公交调度信息 中的每一站到站时间、公交线路编号、公交车辆编号,确定公交车每个对应站点的上车人 数;
步骤四:根据步骤三确定公交车每个对应站点的下车人数。
具体实施方式二:结合图2说明本实施方式,本实施方式与具体实施方式一不同的是: 所述步骤一中根据乘客上车时的刷卡记录时间,求解公交线路的刷卡时间阈值tmin;具体 过程为:
1):将乘客上车时的刷卡记录按刷卡时间从早到晚排列,计算相邻刷卡记录的刷卡时 间差,并将刷卡时间的差值从小到大排序;
2):选取最小的刷卡时间差值作为刷卡时间阈值,根据公交调度信息确定公交车的停 靠站数N;
3):判断相邻两乘客的刷卡时间差值T是否在刷卡时间阈值范围内,若在刷卡时间阈 值tmin范围内,则判断相邻两乘客为同站上车,并记录该站点为Ps;
若不在刷卡时间阈值范围内,则判断相邻两乘客不为同站上车,并记录站点;
直至记下所有乘客的刷卡上车站点,判断出不同站刷卡上车的站点总数Psn;
4):在确定不同站刷卡上车的站点总数之后,进行不同站刷卡上车的相邻两乘客之间 的投币站点数的计算;
不同站刷卡上车的相邻两乘客之间的投币站点计算公式为:Pt=[T/2tmin],直至记下 所有不同站刷卡上车的相邻两乘客之间的投币站点Ptn;
5):由步骤3)和步骤4)得出上车总站数为p=psn+ptn,而已知判断p与N的关系, 当p<N,刷卡时间阈值即为所选择的刷卡时间差阈值tmin;N为根据公交调度信息确定 公交车的停靠站数,取值为正整数;当p≥N,表示刷卡时间阈值选取的值过小,选取次 小(依次选取第三小、第四小、第…小的刷卡时间差)小的刷卡时间差(刷卡时间差中的 较大值),转至步骤3)继续判断p与N的关系,直至满足p<N。
流程图见图2。
其它步骤及参数与具体实施方式一相同。
具体实施方式三:结合图3说明本实施方式,本实施方式与具体实施方式一或二不同 的是:所述步骤二中根据公交线路的刷卡时间阈值tmin,确定公交车未知对应站点时的每 个站点上车人数;具体过程为:
以tmin为界,判断乘客同站或异站上车,确定了如下各站点上车人数的步骤,用Pathon 语言编写了相应的程序来进一步求解各站点上车人数。
站点上车乘客数计算模型:
乘客上车的刷卡记录按刷卡时间早晚进行排序,记kn为刷卡记录。
(1)设初始刷卡记录次数n=0;初始第m个站点m=0;初始第m个站点的上车人 数Nm=0;
n为刷卡记录次数,m为第m个站点,Nm为第m个站点的上车人数;n、m、Nm取 值为正整数;
(2)令刷卡记录次数n=n+1;
(3)将第n次刷卡记录对应的站点记为m站点;
(4)计算第n次刷卡记录与第n-1次刷卡记录的刷卡时间差值为Tn-1;
(5)当Tn-1≤tmin时,则第n次刷卡记录与第n-1次刷卡记录的相邻两乘客为同站上车,则记第n次刷卡记录对应的站点为m=m+0站点;转至(2);
当Tn-1>tmin时,则第n次刷卡记录与第n-1次刷卡记录的相邻两乘客为异站上车,则 记第n次刷卡记录对应的站点为m=m+1站点;记录下m,Nm;
(6)判断第n次刷卡记录是否为最后一条记录,若不是,转至(2);若是则结束。
求解过程如图3。
其它步骤及参数与具体实施方式一或二相同。
具体实施方式四:本实施方式与具体实施方式一至三之一不同的是:所述步骤三中根 据公交车未知对应站点时的每个站点上车人数、IC卡信息、公交调度信息中的每一站到 站时间、公交线路编号以及公交车辆编号,确定公交车每个对应站点的上车人数;具体过 程为:
将公交车未知对应站点时的每个站点上车人数、IC卡信息、公交调度信息中的到站时 间、公交线路编号以及公交车辆编号进行匹配,若满足公式(1),则站点匹配成功,将步 骤二得到的未知对应站点时的每个站点上车人数按前后顺序依次与匹配成功站点对应,得 到公交车每个对应站点的上车人数;若不满足公式(1),则站点匹配失败;
其中,tmid为每一个站点的第一个人和最后一个人的IC卡刷卡时刻平均值(比如第一 个人IC卡刷卡时刻为11点,最后一个人的IC卡刷卡时刻为11点10分,tmid为11点5 分),Δt为首末站公交IC卡的刷卡时间差与公交调度首末站的发车和到站时间差的差值, N为根据公交调度信息确定公交车的停靠站数N,为每个站点的时间误差,tp为由公 交调度信息得到的公交每一站的停靠站时刻(当第三站tp为9点30分,其中一站 为9点20分,为9点40分,则匹配站点为第三站,当第五 站tp为10点,其中一站为9点50分,为10点10分,则匹 配站点为第五站,当第七站tp为10点40分,其中一站为10点30分, 为10点50分,则匹配站点为第七站,将步骤二得到的未知对应3个站点的 上车人数依次为第三站、第五站、第七站的上车人数)。
其它步骤及参数与具体实施方式一至三之一相同。
具体实施方式五:本实施方式与具体实施方式一至四之一不同的是:所述步骤四中根 据步骤三确定公交车每个对应站点的下车人数;具体过程为:
结合集计层面和非集计层面两方面进行下车站点的推导,联系乘客自身的上下车站点 规律,综合考虑站点本身的吸引率推导乘客的下车站点,并建立假设条件与集合:
假设一:若乘客下一次的上车站点位于本次乘坐公交线路的下游站点,则下一次的上 车站点为本次的下车站点;
假设二:若乘客下一次的上车站点位于本次乘坐公交线路的下游站点100米范围内的 站点,则本次下车站点即为离该上车站点最近的下游站点;
假设三:若乘客下一次的上车站点不是本次乘坐公交线路的下游站点,则乘客会选择 D日内出行的高频上车站点下车;
并且在该高频站点上车的概率与在该高频站点上车的次数成正比。
假设四:若乘客下一次的上车站点不是本次乘坐公交线路的下游站点且无法获取D日 内出行最高频站点时,则乘客会选择下游站点的高客流量站点下车;
并且在该高客流量站点上车的概率与在站点上车的客流量成正比。
集合C1={i+1,i+2…,N},即沿k路公交线路,位于站点i下游方向的公交站点组成的集合;
集合C2={(i+1)',(i+2)',…N'},集合C2为距离集合C1100米范围内的周边公交站点 集合;N'为根据公交调度信息确定公交车的停靠站数;
集合C3:参照乘客D日公交出行链,获得最高频上车站点集;D取值为正整数;
为避免将错误的上车站点归到高频站点集,高频站点集中的站点上车次数应大于一定 的值。
集合C4:参照乘客D日公交出行链,获得下游站点的最高客流量站点;D取值为正整数;若站点即乘客下一次上车站点位于本次乘坐公交线路的下游站点,其公交出行链处于闭合状态,如图4所示,则乘客在站点j的下车概率为:
若站点即认为乘客在距离集合C2100米范围内的下游站点下车,而后换乘至 C2站点中的站点h上车,此时其公交出行链处于微闭合状态,如图5所示,令距离站点j100米范围内的站点为j',则乘客在站点j的下车概率为:
若站点则此时其公交出行链出现断裂,如图6所示,即认为此时需要从集计层面的角度分析在站点j的下车概率,若站点j位于乘客的多日出行高频站点集 中,则其在该站点下车的概率与在该站点上车的次数成正比,若站点j为非高频站点集, 则其在该站点的下车概率与该站点的上车人数呈正比,则乘客在站点j的下车概率为:
其中,Mj为乘客D日出行在站点j上车的次数,Vj为站点j的一天上车人数,d为 集合C1∩C3中的站点个数,a为下游站点个数,l为第l个站点,Ml为乘客D日出行在站 点l上车的次数,Vl为站点l的一天上车人数;D取值为正整数;b为集合C4中的站点个 数;
根据下车概率得出公交车每个对应站点的下车人数。
根据权利要求1所述的刷客时间阈值的求解,其特征在刷卡时间阈值用于判断刷卡时 间对应的上车站点是否为同站。当公交上两乘客之间上车刷卡时间差值大于该阈值时,表 示两者刷卡时间对应的刷卡站点为异站。
而当上车刷卡时间差值小于该阈值时,其停靠站数目偏小;当该阈值选取过小时,使 其停靠站数目偏大。
其它步骤及参数与具体实施方式一至四之一相同。
采用以下实施例验证本发明的有益效果:
实施例一:
本实施例一种基于公交IC卡数据的乘客上下车站点推导方法具体是按照以下步骤制 备的:
广州市居民出行使用公交IC卡已相当普及,据最新的资料显示,目前广州市公交IC 卡的使用率已达到80%,以广州市37路为案例对时刻表的优化进行具体分析,37路的起点为东山口总站,终点为赤沙总站,一共16个站点,在研究时段10:00~11:00内,原发车 时间为10:20、10:40、11:00,全程运行时间为50分钟。以2015.01.01当天为例,37路公 交的IC卡刷卡条数已到达1.5万条,刷卡数量非常大。
本发明案例优化10:00~11:00时段37路东山口站发车时刻,因此需分析10:00~11:00 时段从东山口站往赤沙总站的客流需求,从数据库中提取出该时段此运行方向的公交车辆 刷卡记录,其中在该时段共有端号为02560303、09500496、03760206的三趟公交车辆,发 车时刻分别为10:20、10:40、11:00,其中截取端号为02560303的公交车辆部分刷卡数据, 如表3所示,以此为例阐述时间阈值的推导过程,将该部分数据导入上车站点推导算法的 程序中,得到刷卡时间阈值为89s,并且与表1中的公交站点信息相匹配,得到各站点上车 乘客数。
而对于乘客选择何站下车则分下面两种情况进行讨论。
(1)第一种情况 当乘客下一次上车站点位于本线路下游站点或位于下游站点附近站点时,则下一次上车站点为本次下车站点,而在本案例中02560303端号车辆 10:00~11:00的刷卡记录符合该情况的占42%,以表3中乘客的刷卡记录为例,通过该方法 求得乘客下车站点部分刷卡记录如表4所示。
(2)第二种情况 当不满足第一种情况时,寻找乘客的近日多次上车站点,其作为参考站点的发生时间为2015.01.02~2015.01.07,其中频次大于等于2次,根据站点近日上车次数求算出乘客在各站点的下车概率,02560303端号车辆10:00~11:00的刷卡记录符合该情况的占22%,以表3中乘客的刷卡记录为例,得到部分乘客近日多次上车站点集,见 表5,其中/X表示该站点的上车次数。
而剩余的乘客,其上车站点则根据高站点的客流量作为权重进行分配,高客流量站点 分别为五羊新村、客村立交、珠影、大塘,下车吸引率根据其对应上车人数得到,分别为 0.252、0.260、0.244、0.244。因此,可根据高客流站点的下车吸引率计算得到所对应的下 车人数。
最后得到在10:00~11:00期间,端号为02560303的公交车辆在各公交站点的上下车乘 客数,通过实地调查,各站点刷卡率保持在77%~81%,与城市总体80%的刷卡率接近,忽 略各站点自身刷卡率的浮动情况,将各公交站点按城市总体刷卡率80%的比例换算成站点 实际上下车人数,见表2。
表1 37路公交站点信息
表2各站点上下车乘客数(人)
表2(续表)
表3端号02560206车辆刷卡数据
卡号 | 卡类型 | 终端号 | 交易时间 | 上车站点 | 下车站点 |
’999999****320704 | 普通卡 | ’09****96 | ’2015****112153 | 东山口总站 | 五羊新村 |
’999999****742098 | 老人卡 | ’51****00 | ’2015****122414 | 东山口总站 | 江海大道中 |
’999999****001414 | 普通卡 | ’09****96 | ’2015****105241 | 东山口总站 | 五羊新村 |
’999999****954893 | 普通卡 | ’51****00 | ’2015****121321 | 农林东 | 珠影 |
’999999****102444 | 普通卡 | ’00****04 | ’2015****121211 | 农林东 | 石榴岗 |
’999999****807924 | 老人卡 | ’12****01 | ’2015****113823 | 农林东 | 赤岗 |
’999999****826506 | 普通卡 | ’09****96 | ’2015****110136 | 农林东 | 石羊新村 |
’999999****745213 | 普通卡 | ’05****96 | ’2015****121251 | 梅花村 | 大塘 |
’999999****472926 | 普通卡 | ’60****14 | ’2015****121513 | 梅花村 | 赤岗 |
’999999****306520 | 普通卡 | ’09****96 | ’2015****120251 | 梅花村 | 五羊新村 |
’999999****093870 | 普通卡 | ’60****14 | ’2015****113631 | 杨箕村 | 石榴岗 |
’999999****932813 | 普通卡 | ’00****04 | ’2015****105429 | 杨箕村 | 石榴岗 |
’999999****445990 | 普通卡 | ’05****96 | ’2015****111723 | 南方报社 | 大塘 |
’999999****504634 | 普通卡 | ’02****35 | ’2015****113454 | 南方报社 | 台涌 |
’999999****360316 | 普通卡 | ’04****90 | ’2015****105942 | 五羊新村 | 江海大道中 |
’999999****126203 | 普通卡 | ’00****04 | ’2015****111234 | 五羊新村 | 石榴岗 |
’999999****252950 | 普通卡 | ’06****96 | ’2015****115451 | 五羊新村 | 赤岗 |
’999999****215770 | 普通卡 | ’02****15 | ’2015****114111 | 客村立交 | 珠影 |
’999999****163680 | 普通卡 | ’09****96 | ’2015****10423 | 客村立交 | 五羊新村 |
’999999****620101 | 普通卡 | ’05****96 | ’2015****121836 | 珠影 | 台涌 |
表4 满足情况一的部分乘客站点信息
使用地 | 发卡地 | 卡号 | 卡类型 | 终端号 | 交易时间 |
广州 | 广州 | ’999999****007628 | 普通卡 | ’02****03 | ’2015****102140 |
广州 | 广州 | ’999999****500056 | 普通卡 | ’02****03 | ’2015****102158 |
广州 | 广州 | ’999999****320704 | 普通卡 | ’02****03 | ’2015****102158 |
广州 | 广州 | ’999999****742098 | 老人卡 | ’02****03 | ’2015****102159 |
广州 | 广州 | ’999999****902736 | 普通卡 | ’02****03 | ’2015****102201 |
广州 | 广州 | ’999999****720052 | 普通卡 | ’02****03 | ’2015****102202 |
广州 | 广州 | ’999999****001414 | 普通卡 | ’02****03 | ’2015****102206 |
广州 | 广州 | ’999999****378016 | 普通卡 | ’02****03 | ’2015****102339 |
广州 | 广州 | ’999999****954893 | 普通卡 | ’02****03 | ’2015****102339 |
广州 | 广州 | ’999999****102444 | 普通卡 | ’02****03 | ’2015****102341 |
广州 | 广州 | ’999999****807924 | 老人卡 | ’02****03 | ’2015****102343 |
广州 | 广州 | ’999999****008511 | 普通卡 | ’02****03 | ’2015****102344 |
广州 | 广州 | ’999999****826506 | 普通卡 | ’02****03 | ’2015****102346 |
广州 | 广州 | ’999999****924810 | 老人卡 | ’02****03 | ’2015****102347 |
广州 | 广州 | ’999999****745213 | 普通卡 | ’02****03 | ’2015****102454 |
广州 | 广州 | ’999999****121305 | 老人卡 | ’02****03 | ’2015****102500 |
广州 | 广州 | ’999999****472926 | 普通卡 | ’02****03 | ’2015****102503 |
广州 | 广州 | ’999999****851045 | 老人卡 | ’02****03 | ’2015****102504 |
广州 | 广州 | ’999999****422843 | 普通卡 | ’02****03 | ’2015****102506 |
广州 | 广州 | ’999999****306520 | 普通卡 | ’02****03 | ’2015****102508 |
广州 | 广州 | ’999999****908240 | 老人卡 | ’02****03 | ’2015****102633 |
广州 | 广州 | ’999999****882750 | 老人卡 | ’02****03 | ’2015****102637 |
广州 | 广州 | ’999999****093870 | 普通卡 | ’02****03 | ’2015****102638 |
广州 | 广州 | ’999999****932813 | 普通卡 | ’02****03 | ’2015****102640 |
广州 | 广州 | ’999999****285815 | 老人卡 | ’02****03 | ’2015****102642 |
广州 | 广州 | ’999999****085805 | 老人卡 | ’02****03 | ’2015****102645 |
广州 | 广州 | ’999999****360048 | 普通卡 | ’02****03 | ’2015****102713 |
广州 | 广州 | ’999999****445990 | 普通卡 | ’02****03 | ’2015****102721 |
广州 | 广州 | ’999999****433080 | 普通卡 | ’02****03 | ’2015****102723 |
广州 | 广州 | ’999999****441584 | 普通卡 | ’02****03 | ’2015****102732 |
广州 | 广州 | ’999999****504634 | 普通卡 | ’02****03 | ’2015****102740 |
广州 | 广州 | ’999999****422759 | 普通卡 | ’02****03 | ’2015****102741 |
广州 | 广州 | ’999999****055483 | 普通卡 | ’02****03 | ’2015****102759 |
广州 | 广州 | ’999999****103864 | 普通卡 | ’02****03 | ’2015****102802 |
广州 | 广州 | ’999999****360316 | 普通卡 | ’02****03 | ’2015****102813 |
广州 | 广州 | ’999999****126203 | 普通卡 | ’02****03 | ’2015****102814 |
广州 | 广州 | ’999999****456900 | 普通卡 | ’02****03 | ’2015****102816 |
广州 | 广州 | ’999999****807647 | 老人卡 | ’02****03 | ’2015****102822 |
广州 | 广州 | ’999999****391002 | 学生卡 | ’02****03 | ’2015****102823 |
广州 | 广州 | ’999999****252950 | 普通卡 | ’02****03 | ’2015****102826 |
广州 | 广州 | ’999999****883315 | 学生卡 | ’02****03 | ’2015****102933 |
广州 | 广州 | ’999999****422627 | 普通卡 | ’02****03 | ’2015****102936 |
广州 | 广州 | ’999999****008815 | 老人卡 | ’02****03 | ’2015****102939 |
广州 | 广州 | ’999999****215770 | 普通卡 | ’02****03 | ’2015****102942 |
广州 | 广州 | ’999999****163680 | 普通卡 | ’02****03 | ’2015****102947 |
广州 | 广州 | ’999999****974014 | 普通卡 | ’02****03 | ’2015****102951 |
表5 部分乘客的近日多次上车站点信息
使用地 | 发卡地 | 卡号 | 卡类型 | 多次上车站点1 | 多次站点2 |
广州 | 广州 | ’999999****500056 | 普通卡 | 五羊新村/3 | 赤岗路口/2 |
广州 | 广州 | ’999999****902736 | 普通卡 | 南方报社/4 | |
广州 | 广州 | ’999999****008511 | 普通卡 | 五羊新村/2 | 石榴岗路西/2 |
广州 | 广州 | ’999999****851045 | 普通卡 | 江海大道/3 | |
广州 | 广州 | ’999999****882750 | 老人卡 | 大塘/3 | 杨萁村/2 |
广州 | 广州 | ’999999****055483 | 普通卡 | 石榴岗路西/3 | |
广州 | 广州 | ’999999****103864 | 普通卡 | 石榴岗/2 | |
广州 | 广州 | ’99999****0422627 | 普通卡 | 江海大道中/3 | |
广州 | 广州 | ’999999****529902 | 普通卡 | 珠影/2 | 石榴岗/2 |
本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术 人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发 明所附的权利要求的保护范围。
Claims (5)
1.一种基于公交IC卡数据的站点上下车乘客数估计方法,其特征在于:所述方法具体过程为:
步骤一:根据乘客上车时的刷卡记录时间,求解公交线路的刷卡时间阈值tmin;
步骤二:根据公交线路的刷卡时间阈值tmin,确定公交车未知对应站点时的每个站点上车人数;
步骤三:根据公交车未知对应站点时的每个站点上车人数、IC卡信息、公交调度信息中的每一站到站时间、公交线路编号以及公交车辆编号,确定公交车每个对应站点的上车人数;
步骤四:根据步骤三确定公交车每个对应站点的下车人数。
2.根据权利要求1所述一种基于公交IC卡数据的站点上下车乘客数估计方法,其特征在于:所述步骤一中根据乘客上车时的刷卡记录时间,求解公交线路的刷卡时间阈值tmin;具体过程为:
1):将乘客上车时的刷卡记录按刷卡时间从早到晚排列,计算相邻刷卡记录的刷卡时间差,并将刷卡时间的差值从小到大排序;
2):选取最小的刷卡时间差值作为刷卡时间阈值,根据公交调度信息确定公交车的停靠站数N;
3):判断相邻两乘客的刷卡时间差值T是否在刷卡时间阈值范围内,若在刷卡时间阈值tmin范围内,则判断相邻两乘客为同站上车,并记录该站点为Ps;
若不在刷卡时间阈值范围内,则判断相邻两乘客不为同站上车,并记录站点;
直至记下所有乘客的刷卡上车站点,判断出不同站刷卡上车的站点总数Psn;
4):在确定不同站刷卡上车的站点总数之后,进行不同站刷卡上车的相邻两乘客之间的投币站点数的计算;
不同站刷卡上车的相邻两乘客之间的投币站点计算公式为:Pt=[T/2tmin],直至记下所有不同站刷卡上车的相邻两乘客之间的投币站点Ptn;
5):由步骤3)和步骤4)得出上车总站数为p=psn+ptn,当p<N,刷卡时间阈值即为所选择的刷卡时间差阈值tmin;N为根据公交调度信息确定公交车的停靠站数,取值为正整数;当p≥N,表示刷卡时间阈值选取的值过小,选取次小的刷卡时间差,转至步骤3)继续判断p与N的关系,直至满足p<N。
3.根据权利要求2所述一种基于公交IC卡数据的站点上下车乘客数估计方法,其特征在于:所述步骤二中根据公交线路的刷卡时间阈值tmin,确定公交车未知对应站点时的每个站点上车人数;具体过程为:
(1)设初始刷卡记录次数n=0;初始第m个站点m=0;初始第m个站点的上车人数Nm=0;
n为刷卡记录次数,m为第m个站点,Nm为第m个站点的上车人数;n、m、Nm取值为正整数;
(2)令刷卡记录次数n=n+1;
(3)将第n次刷卡记录对应的站点记为m站点;
(4)计算第n次刷卡记录与第n-1次刷卡记录的刷卡时间差值为Tn-1;
(5)当Tn-1≤tmin时,则第n次刷卡记录与第n-1次刷卡记录的相邻两乘客为同站上车,则记第n次刷卡记录对应的站点为m=m+0站点;转至(2);
当Tn-1>tmin时,则第n次刷卡记录与第n-1次刷卡记录的相邻两乘客为异站上车,则记第n次刷卡记录对应的站点为m=m+1站点;Nm=n-∑Nm-1,记录下m,Nm;
(6)判断第n次刷卡记录是否为最后一条记录,若不是,转至(2);若是则结束。
4.根据权利要求3所述一种基于公交IC卡数据的站点上下车乘客数估计方法,其特征在于:所述步骤三中根据公交车未知对应站点时的每个站点上车人数、IC卡信息、公交调度信息中的每一站到站时间、公交线路编号以及公交车辆编号,确定公交车每个对应站点的上车人数;具体过程为:
将公交车未知对应站点时的每个站点上车人数、IC卡信息、公交调度信息中的每一站到站时间、公交线路编号以及公交车辆编号进行匹配,若满足公式(1),则站点匹配成功,将步骤二得到的未知对应站点时的每个站点上车人数按前后顺序依次与匹配成功站点对应,得到公交车每个对应站点的上车人数;若不满足公式(1),则站点匹配失败;
<mrow>
<msub>
<mi>t</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>-</mo>
<mfrac>
<msub>
<mi>t</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mn>2</mn>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
<mi>N</mi>
</mfrac>
<mo><</mo>
<msub>
<mi>t</mi>
<mi>p</mi>
</msub>
<mo><</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<msub>
<mi>t</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mn>2</mn>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
<mi>N</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,tmid为每一个站点的第一个人和最后一个人的IC卡刷卡时刻平均值,Δt为首末站公交IC卡的刷卡时间差与公交调度首末站的发车和到站时间差的差值,N为根据公交调度信息确定公交车的停靠站数,为每个站点的时间误差,tp为由公交调度信息得到的公交每一站的停靠站时刻。
5.根据权利要求4所述一种基于公交IC卡数据的站点上下车乘客数估计方法,其特征在于:所述步骤四中根据步骤三确定公交车每个对应站点的下车人数;具体过程为:
建立假设条件与集合:
假设一:若乘客下一次的上车站点位于本次乘坐公交线路的下游站点,则下一次的上车站点为本次的下车站点;
假设二:若乘客下一次的上车站点位于本次乘坐公交线路的下游站点100米范围内的站点,则本次下车站点即为离该上车站点最近的下游站点;
假设三:若乘客下一次的上车站点不是本次乘坐公交线路的下游站点,则乘客会选择D日内出行的最高频上车站点下车;
假设四:若乘客下一次的上车站点不是本次乘坐公交线路的下游站点且无法获取D日内出行最高频站点时,则乘客会选择下游站点的最高客流量站点下车;
集合C1={i+1,i+2…,N},即沿k路公交线路,位于站点i下游方向的公交站点组成的集合;
集合C2={(i+1)',(i+2)',…N'},集合C2为距离集合C1100米范围内的公交站点集合;N'为根据公交调度信息确定公交车的停靠站数;
集合C3:参照乘客D日公交出行链,获得最高频上车站点集;D取值为正整数;
集合C4:参照乘客D日公交出行链,获得下游站点的最高客流量站点;D取值为正整数;
若站点即乘客下一次上车站点位于本次乘坐公交线路的下游站点,则乘客在站点j的下车概率为:
<mrow>
<msub>
<mi>P</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
<mo>,</mo>
<mi>h</mi>
<mo>=</mo>
<mi>j</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
<mo>,</mo>
<mi>h</mi>
<mo>&NotEqual;</mo>
<mi>j</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
若站点即认为乘客在距离集合C2100米范围内的下游站点下车,而后换乘至C2站点中的站点h上车,令距离站点j 100米范围内的站点为j',则乘客在站点j的下车概率为:
<mrow>
<msub>
<mi>P</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
<mo>,</mo>
<mi>h</mi>
<mo>=</mo>
<msup>
<mi>j</mi>
<mo>&prime;</mo>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
<mo>,</mo>
<mi>h</mi>
<mo>&NotEqual;</mo>
<msup>
<mi>j</mi>
<mo>&prime;</mo>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
若站点则乘客在站点j的下车概率为:
<mrow>
<msub>
<mi>P</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>M</mi>
<mi>j</mi>
</msub>
<mo>/</mo>
<mstyle>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>d</mi>
</munderover>
<mrow>
<msub>
<mi>M</mi>
<mi>l</mi>
</msub>
<mo>,</mo>
<mi>j</mi>
<mo>&Subset;</mo>
<mo>{</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mo>&cap;</mo>
<msub>
<mi>C</mi>
<mn>3</mn>
</msub>
<mo>}</mo>
</mrow>
</mstyle>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mi>j</mi>
</msub>
<mo>/</mo>
<mstyle>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>b</mi>
</munderover>
<mrow>
<msub>
<mi>V</mi>
<mi>l</mi>
</msub>
<mo>,</mo>
<mi>j</mi>
<mo>&Subset;</mo>
<mo>{</mo>
<msub>
<mi>C</mi>
<mn>4</mn>
</msub>
<mo>}</mo>
</mrow>
</mstyle>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
<mo>,</mo>
<mi>j</mi>
<mo>&NotSubset;</mo>
<mo>{</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mo>&cap;</mo>
<msub>
<mi>C</mi>
<mn>3</mn>
</msub>
<mo>}</mo>
<mo>&cup;</mo>
<mo>{</mo>
<msub>
<mi>C</mi>
<mn>4</mn>
</msub>
<mo>}</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,Mj为乘客D日出行在站点j上车的次数,Vj为站点j的一天上车人数,d为集合C1∩C3中的站点个数,a为下游站点个数,l为第l个站点,Ml为乘客D日出行在站点l上车的次数,Vl为站点l的一天上车人数;D取值为正整数;b为集合C4中的站点个数;
根据下车概率得出公交车每个对应站点的下车人数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710806799.7A CN107545730B (zh) | 2017-09-08 | 2017-09-08 | 一种基于公交ic卡数据的站点上下车乘客数估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710806799.7A CN107545730B (zh) | 2017-09-08 | 2017-09-08 | 一种基于公交ic卡数据的站点上下车乘客数估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107545730A true CN107545730A (zh) | 2018-01-05 |
CN107545730B CN107545730B (zh) | 2020-04-07 |
Family
ID=60958127
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710806799.7A Expired - Fee Related CN107545730B (zh) | 2017-09-08 | 2017-09-08 | 一种基于公交ic卡数据的站点上下车乘客数估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107545730B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109102114A (zh) * | 2018-08-02 | 2018-12-28 | 中南大学 | 一种基于数据融合的公交出行下车站点估计方法 |
CN109409563A (zh) * | 2018-09-07 | 2019-03-01 | 北明软件有限公司 | 一种公交运营车辆实时人数的分析方法、系统和存储介质 |
CN109584555A (zh) * | 2018-12-13 | 2019-04-05 | 昆山市公共交通集团有限公司 | 基于afc数据的公交乘客下车站点推测方法 |
CN110751830A (zh) * | 2019-09-30 | 2020-02-04 | 江苏大学 | 一种集计公交物理站点平均等车时间的方法及系统 |
CN111339159A (zh) * | 2020-02-24 | 2020-06-26 | 交通运输部科学研究院 | 一种一票制公交数据的分析挖掘方法 |
CN111491261A (zh) * | 2020-04-24 | 2020-08-04 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 一种基于智能刷卡数据的个体移动轨迹提取方法 |
CN111723871A (zh) * | 2020-07-09 | 2020-09-29 | 广州市公共交通数据管理中心 | 一种公交车实时车厢满载率的估算方法 |
CN111753217A (zh) * | 2020-06-24 | 2020-10-09 | 武汉元光科技有限公司 | 一种生成公交站台引导词的方法和装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101540098A (zh) * | 2009-04-30 | 2009-09-23 | 北京工业大学 | 基于公交ic卡线路的匹配方法 |
CN105023437A (zh) * | 2015-08-21 | 2015-11-04 | 苏州大学张家港工业技术研究院 | 一种公交od矩阵的构建方法及系统 |
CN105654721A (zh) * | 2016-01-21 | 2016-06-08 | 哈尔滨工业大学 | 一种采用公交ic卡和车载gps数据计算公交od量的方法 |
CN105702035A (zh) * | 2016-03-23 | 2016-06-22 | 武汉大学 | 一种利用历史公交数据评估乘车难易程度的方法 |
JP5986641B2 (ja) * | 2012-10-17 | 2016-09-06 | 株式会社日立製作所 | 交通分析システム |
CN106504525A (zh) * | 2016-10-12 | 2017-03-15 | 电子科技大学 | 基于ic卡数据的od矩阵生成技术及其应用研究 |
-
2017
- 2017-09-08 CN CN201710806799.7A patent/CN107545730B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101540098A (zh) * | 2009-04-30 | 2009-09-23 | 北京工业大学 | 基于公交ic卡线路的匹配方法 |
JP5986641B2 (ja) * | 2012-10-17 | 2016-09-06 | 株式会社日立製作所 | 交通分析システム |
CN105023437A (zh) * | 2015-08-21 | 2015-11-04 | 苏州大学张家港工业技术研究院 | 一种公交od矩阵的构建方法及系统 |
CN105654721A (zh) * | 2016-01-21 | 2016-06-08 | 哈尔滨工业大学 | 一种采用公交ic卡和车载gps数据计算公交od量的方法 |
CN105702035A (zh) * | 2016-03-23 | 2016-06-22 | 武汉大学 | 一种利用历史公交数据评估乘车难易程度的方法 |
CN106504525A (zh) * | 2016-10-12 | 2017-03-15 | 电子科技大学 | 基于ic卡数据的od矩阵生成技术及其应用研究 |
Non-Patent Citations (1)
Title |
---|
闫磊: "基于公交IC卡数据的乘客出行时空推算研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109102114A (zh) * | 2018-08-02 | 2018-12-28 | 中南大学 | 一种基于数据融合的公交出行下车站点估计方法 |
CN109409563B (zh) * | 2018-09-07 | 2021-11-09 | 北明软件有限公司 | 一种公交运营车辆实时人数的分析方法、系统和存储介质 |
CN109409563A (zh) * | 2018-09-07 | 2019-03-01 | 北明软件有限公司 | 一种公交运营车辆实时人数的分析方法、系统和存储介质 |
CN109584555A (zh) * | 2018-12-13 | 2019-04-05 | 昆山市公共交通集团有限公司 | 基于afc数据的公交乘客下车站点推测方法 |
CN110751830A (zh) * | 2019-09-30 | 2020-02-04 | 江苏大学 | 一种集计公交物理站点平均等车时间的方法及系统 |
CN111339159A (zh) * | 2020-02-24 | 2020-06-26 | 交通运输部科学研究院 | 一种一票制公交数据的分析挖掘方法 |
CN111339159B (zh) * | 2020-02-24 | 2023-08-18 | 交通运输部科学研究院 | 一种一票制公交数据的分析挖掘方法 |
CN111491261A (zh) * | 2020-04-24 | 2020-08-04 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 一种基于智能刷卡数据的个体移动轨迹提取方法 |
CN111491261B (zh) * | 2020-04-24 | 2022-03-01 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 一种基于智能刷卡数据的个体移动轨迹提取方法 |
CN111753217A (zh) * | 2020-06-24 | 2020-10-09 | 武汉元光科技有限公司 | 一种生成公交站台引导词的方法和装置 |
CN111753217B (zh) * | 2020-06-24 | 2021-10-08 | 武汉元光科技有限公司 | 一种生成公交站台引导词的方法和装置 |
CN111723871A (zh) * | 2020-07-09 | 2020-09-29 | 广州市公共交通数据管理中心 | 一种公交车实时车厢满载率的估算方法 |
CN111723871B (zh) * | 2020-07-09 | 2022-05-27 | 广州市公共交通数据管理中心有限公司 | 一种公交车实时车厢满载率的估算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107545730B (zh) | 2020-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107545730A (zh) | 一种基于公交ic卡数据的站点上下车乘客数估计方法 | |
CN108269399B (zh) | 一种基于公路网客流od反推技术的高铁客流需求预测方法 | |
CN102592447B (zh) | 一种基于fcm的区域路网的道路交通状态判别方法 | |
CN104574967B (zh) | 一种基于北斗的城市大面积路网交通感知方法 | |
CN101615340A (zh) | 公交动态调度中的实时信息处理方法 | |
CN106384509A (zh) | 考虑出租车运营状态下的城市道路行程时间分布估算方法 | |
CN108151756A (zh) | 一种城市公共交通系统的行程时间计算方法和装置 | |
CN103226801B (zh) | 基于多用户配流模型的空港集疏运交通量确定方法 | |
CN111582750B (zh) | 一种轨道交通反向乘客识别及乘车班次确定方法及系统 | |
CN106898142B (zh) | 一种考虑路段相关性的路径行程时间可靠度计算方法 | |
CN109903553B (zh) | 多源数据挖掘的公交车上下车站点识别和检验方法 | |
CN102324128A (zh) | 基于ic卡记录的公交站间od客流预测方法及装置 | |
CN105206040B (zh) | 一种基于ic卡数据的公交串车预测方法 | |
CN102306243A (zh) | 在轨道交通中的换乘客流清分方法 | |
CN107203945A (zh) | 车险分级评估方法及装置 | |
CN105320997A (zh) | 一种基于大数据的公交客流轨迹实时跟踪方法 | |
CN108665140A (zh) | 一种城际铁路客运枢纽交通衔接系统评价方法 | |
CN106651728A (zh) | 一种综合运输体系客运方式优势运距的确定方法 | |
CN101901293A (zh) | 一种基于多种因素的轨道交通旅行时间计算模型 | |
CN106846214A (zh) | 分析交通枢纽可达性对区域公共运输方式竞争影响的方法 | |
CN112185580B (zh) | 一种基于云端的行程风险信息处理方法 | |
CN107180534A (zh) | 支持向量回归机融合的高速公路路段平均速度估计方法 | |
CN107316096A (zh) | 一种轨道交通一票通乘客进站量预测方法 | |
Huang et al. | Geographic distribution of firm productivity and production: A “, market access” approach | |
CN113408833A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200407 Termination date: 20200908 |
|
CF01 | Termination of patent right due to non-payment of annual fee |