CN102097002A

CN102097002A - 一种基于ic卡数据获取公交站点od的方法及系统

Info

Publication number: CN102097002A
Application number: CN 201010553014
Authority: CN
Inventors: 陈学武; 陈峥嵘; 杨应翔; 张颂
Original assignee: Southeast University
Current assignee: Southeast University
Priority date: 2010-11-22
Filing date: 2010-11-22
Publication date: 2011-06-15
Anticipated expiration: 2030-11-22
Also published as: CN102097002B

Abstract

本发明公布了一种基于IC卡数据获取公交站点OD的方法及系统。所述的方法为：将采集得到的IC卡刷卡数据进行班次聚类、站点聚类、上车站点匹配、下车站点匹配等步骤，获得公交站点OD矩阵。所述的系统，则包括：数据采集模块、数据存储模块、IC卡刷卡数据处理模块、站点OD信息发布模块。本发明为公交调度和公交建模提供及时的站点OD信息，相比现有站点OD获取技术具有准确、实时的特点。

Description

一种基于IC卡数据获取公交站点OD的方法及系统

技术领域

本发明属于城市公共交通领域，尤其涉及一种基于公交IC卡数据获取公交站点OD的方法及系统。

背景技术

公交站点OD是指公交线网中自上车站点O至下车站点D的乘客量，它在城市公交规划、公交线网优化以及公交运营调度中都有着广泛的应用。传统的获取公交站点OD的方法是通过大规模的抽样调查，如乘客问卷调查、跟车调查等。实施这些调查需要投入大量人力和物力，因此在资金方面往往令人难以接受；由于现场调查工作量非常大，且调查员业务水平参差不齐，获得的公交站点OD矩阵与实际情况有偏差；另外由于从调查实施到公交站点OD矩阵得出需要较长一段时间且不可重复操作，因此得来的调查数据生命周期很短且不能及时反映实际的公交出行情况。

公交数据自动采集设备(如公交IC卡、GPS等)在国内一些大城市的应用，使得公交站点OD的获取方法不再局限于传统的人工调查，通过结合公交车辆GPS轨迹与乘客IC卡刷卡数据便可得到乘客上车站点等信息从而推导出公交站点OD矩阵。然而受诸多因素限制，现阶段大部分城市的大部分公交线路尚未配备GPS设备，只能得到公交IC卡上车刷卡数据。在此前提下，已有的利用IC卡数据获取公交站点OD技术只停留于模糊聚类分析，不能准确判断乘客上车站点、下车站点等信息，由此推导出的公交站点OD矩阵与实际情况存在较大偏差。

发明内容

本发明针对现有公交OD获取技术存在的不足，提出一种基于IC卡数据获取公交站点OD的方法及系统，本发明的技术方案如下：

一种基于IC卡数据获取公交站点OD的方法，其特征在于先用数据采集模块采集一个城市或城市中某个区域内持有公交IC卡乘客的IC卡刷卡数据，所述采集到的刷卡数据先存储在存储单元中；再由IC卡刷卡数据处理模块调取存储单元中的刷卡数据进行处理，得到公交站点OD信息、；最后通过站点OD信息发布模块把公交站点OD信息发布在各终端设备上；

所述IC卡刷卡数据处理模块中，对IC卡刷卡数据进行下述步骤的处理：

1)刷卡数据按班次聚类；

2)上车站点匹配，步骤是，(21)刷卡数据按站点聚类；(22)正交换乘站点刷卡数据匹配；(23)其余站点刷卡数据匹配；

3)下车站点匹配：公交乘客在下车站点没有刷卡记录，用正交换乘和回程进行下车站点匹配；

4)依据匹配出刷卡数据的上车站点、下车站点，生成公交站点OD矩阵，导入数据库后存储。

实现上述方法的系统包括：

(1)数据采集模块：公交车辆上装有IC卡刷卡终端设备；

(2)数据存储模块：数据库服务器；

(3)IC卡刷卡数据处理模块：包括依此连接的班次聚类器、上车站点匹配器、下车站点匹配器以及站点OD生成器；具体实现可以用PC来实现；

(4)站点OD信息发布模块：信息发布服务器；

各个模块采用有线或无线网络介质进行数据传输。

与现有技术相比，本技术方案的效果主要体现是，提高公交OD获取的效率和质量。

附图说明

图1为方法的流程图。

图2为系统的结构图。

图3为关系数据库结构图。

图4为正交换乘示意图。

图5为刷卡时刻分布示意图。

图6为未通过正交换乘判断出的站点匹配过程示意图。

具体实施方式

如图2所示，系统由数据采集模块、数据存储模块、IC卡刷卡数据处理模块和站点OD信息发布模块组成。其中，系统的核心部分----IC卡刷卡数据处理模块包括班次聚类器、上车站点匹配器、下车站点匹配器以及站点OD生成器。

公交车辆上装有IC卡刷卡终端设备，车辆行驶时停靠各个公交站台，持卡乘客上车刷卡，终端设备收集乘客刷卡时间等信息并定时传输到系统数据库；数据库按照设计的格式存储数据，供IC卡刷卡数据处理模块和站点OD信息发布模块调用；IC卡刷卡数据处理模块对存储在数据库里的大量数据按照设计的流程进行处理，并将处理结果上报给数据库服务器端；站点OD信息发布模块根据用户的不同需求采用相应的传输方式将站点OD信息发布在各终端设备上。

采集一个城市或城市中某个区域内持有公交IC卡乘客的刷卡数据，并以此数据建立关系数据库，此数据库包括核心数据信息和参数配置信息，核心数据信息包括IC卡刷卡记录表、公交站点编号表，参数配置信息包括线路参数表、正交换乘参数表、相邻站点行程时间临界值表。

本发明中一种基于IC卡数据获取公交站点OD的方法，主要包括五个步骤：

步骤1：对采集到的数据进行预处理，剔除运行时间范围外的无效数据，并按车辆运行时间先后顺序排列刷卡数据；

步骤2：依据步骤1中得到的数据，识别不同车次的刷卡记录；

步骤3：依据步骤2中得到的数据，判断各条线路各乘客的上车站点，包括：

步骤3-1同车次同站点刷卡数据聚类；

步骤3-2对满足正交换乘条件的站点刷卡数据进行上车站点匹配；

步骤3-3对余下站点刷卡数据进行上车站点匹配；

步骤4：依据步骤1和步骤2得到的数据，判断各条线路各乘客的下车站点，包括：

步骤4-1对满足正交换乘条件的站点刷卡数据进行下车站点匹配；

步骤4-2对满足回程条件的站点刷卡数据进行下车站点匹配；

步骤5：依据上述步骤匹配出的乘客刷卡记录的上车站点、下车站点，生成公交站点OD矩阵，导入数据库。

下面参照附图1，对本发明的具体实施方案作详细描述：

1、数据库设计

本发明采用关系数据库(如Access数据库)管理各类数据，内含两大类数据表：一类是核心数据表，包括IC卡刷卡记录表、公交站点编号表；一类是参数配置表，包括线路参数表、正交换乘参数表。其中，IC卡刷卡记录表可直接从外部数据库导入，其余数据表需进行补充调查获取数据并人工录入。数据库结构如图3所示。

(1)核心数据表

核心数据表的IC卡刷卡记录表包括以下字段：记录编号、IC卡卡号、卡类型、刷卡日期、刷卡时刻、线路编号、车辆编号；公交站点编号表包括以下字段：站点编号、线路编号、行车方向、站点名称。各核心数据表结构如下所示。

表11C卡刷卡记录表结构

字段名称	中文名称	字段类型	字段说明
				ID	记录编号	整型	主键、唯一、非空
Cardno	IC卡卡号	整型
				Cardkind	卡类型	整型
Transdate	刷卡日期	长整型	YYYYMMDD
				Transtime	刷卡时刻	长整型	hhmmss
Lineid	线路编号	整型
				Busid	车辆编号	整型

表2公交站点编号表结构

字段名称	中文名称	字段类型	字段说明
				StationID	站点编号	整型	主键、唯一、非空
Lineid	线路编号	整型

Direction	行车方向	整型	1或2
				Stationname	站点名称	字符串	允许中文，长度不限

(2)参数配置表

参数配置表的线路参数表包括以下字段：记录编号、线路编号、行车方向、车辆末站停留时间、线路行程时间、站点行程时间临界值；正交换乘参数表包括以下字段：换乘站点、正交换乘线路、换乘判断临界时间。其中，相邻站点行程时间临界值是指线路内各相邻站点之间的平均最小行驶时间；正交换乘是指某线路有且仅有一个站点供另一条线路的下车乘客换乘，如图4中线路1与线路2互为对方的正交换乘线路；换乘判断临界时间是指在换乘站点换乘耗时与换乘前所乘线路行程时间之和。

表3线路参数表结构

字段名称	中文名称	字段类型	字段说明
				ID	记录编号	整型	主键、唯一、非空
Lineid	线路编号	整型
				Direction	行车方向	整型	1或2
Depint	车辆末站停留时间	浮点型	补充调查得到
				Tratime	线路行程时间	浮点型	补充调查得到
Stacritime	相邻站点行程时间临界值	浮点型	补充调查得到

表4正交换乘参数表结构

字段名称	中文名称	字段类型	字段说明
				Transtation	换乘站点	整型	外键，来自站点编号表
Transline	正交换乘线路	整型	换乘站点所在线路的正交换乘线路的编号
				Transcritime	换乘判断临界时间	浮点型	换乘耗时与换乘前所乘线路行程时间之和

表5站点标准行程时间表结构

字段名称	中文名称	字段类型	字段说明
				Prestation	前行站点	整型	外键，来自站点编号表

Postation	后行站点	整型	外键，来自站点编号表
				Statratime	行程时间	浮点型	两站点间标准行程时间，补充调查得到

2、数据预处理

(1)抽取某天所有的IC卡刷卡数据；

(2)按照车辆编号、刷卡时刻进行排序，得到新编号刷卡数据；

(3)记新编号的刷卡数据的序号是r；比较该编号记录的刷卡时刻和前一编号记录的刷卡时刻，得到二者的时间差Δt；

(4)最后判断刷卡记录是否处理完。

3、按班次聚类

IC卡刷卡数据的班次聚类由班次聚类器完成，其功能和核心算法描述如下。

公交车辆从始发站出发沿途停靠每个站点，到达终点站后停留一段时间，折返。与之对应的刷卡记录数据有以下特征：抽取某日一条线路的所有刷卡数据，按车辆编号、刷卡时刻顺序排列，以刷卡序列号为X轴、刷卡时刻为Y轴建立二维坐标系，此时刷卡时刻大致呈阶梯状分布，且每隔若干个阶梯出现一次较大时间间隔，记为一个GAP，其时间长度为T_GAP，两相邻GAP间的时间间隔记作t_{GAP(i-l)，GAP(i)}，如图5所示。

按顺序进行以下几步判断，将不同班次的刷卡记录区分开来：

①如果T_GAP(i)≥Depint且t_{GAP(i-l)，GAP(i)}∈[Tratime-Tratime/5，Tratime+Tratime/5]，则GAP(i-l)与GAP(i)之间的刷卡记录作为一个班次的完整刷卡记录；

②如果T_GAP(i)≥Depint/2且t_{GAP(i-l)，GAP(i)}∈[Tratime-Tratime/10，Tratime+Tratime/10]，则GAP(i-l)与GAP(i)之间的刷卡记录作为一个班次的完整刷卡记录；

③如果上述判断均不满足，提取GAP(i-l)后面的所有满足t_{GAP(i-l)，GAP(p)}∈[Tratime-Tratime/10，Tratime+Tratime/10]的刷卡记录，按T_GAP(p)排序，其值最大的GAP(p)与已判断出GAP(i-l)之间的刷卡记录作为一个班次的完整刷卡记录。

4、上车站点匹配

IC卡刷卡数据的上车站点匹配由上车站点匹配器完成，其功能和核心算法描述如下。

(1)刷卡数据按站点聚类

公交车辆从始发站出发沿途停靠每个站点，持有公交IC卡乘客在各站点上车刷卡。与之对应的刷卡记录数据有以下特征：抽取某日一条线路的一个班次的所有刷卡数据，按刷卡时刻顺序排列，以刷卡序列号为X轴、刷卡时刻为Y轴建立二维坐标系，此时刷卡时刻大致呈阶梯状分布，相邻阶梯的时间间隔记为一个STATION，其时间长度为T_STATION，两相邻STATION间的时间间隔记作t_{STATION(j-I)、STATION(j)}，如图5所示。

定义一个站点对应的所有刷卡记录为一个类，由于在某些站点可能没有乘客刷卡，因此通常一个班次包含的类的数量小于等于站点的个数。

按顺序进行以下几步判断，完成刷卡数据的站点聚类：

①如果T_STATION(j)≥60s且t_{STATION(j-l)，STATION(j)}∈[Stacritime-Stacritime/5，Stacritime+Stacritime/5]，则STATION(j-l)与STATION(j)之间的刷卡记录作为一个站点的完整刷卡记录：

②如果T_STATION(j)≥30s且t_{STATION(j-I)，STATION(j)}∈[Stacritime-Stacritime/10，Stacritime+Stacritime/10]，则STATION(j-l)与STATION(j)之间的刷卡记录作为一个站点的完整刷卡记录；

③如果上述判断均不满足，提取STATION(j-1)后面的所有满足t_{STATION(j-I)，STATION(q)}∈[Stacritime-Stacritime/10，Stacritime+Stacritime/10]的刷卡记录，按T_STATION(q)排序，其值最大的STATION(q)的刷卡记录作为一个站点的完整刷卡记录。

(2)正交换乘站点刷卡数据匹配

对于通过步骤(1)处理后的数据，抽取某日某条线路以及可将该线路作为正交换乘线路的其他线路的刷卡数据，提取出在该线路的刷卡时刻晚于其他线路刷卡时刻的对应刷卡记录，判断其刷卡间隔是否小于换乘判断临界时间，如满足则标定该线路的此刷卡记录对应的上车站点为对应换乘站点，完成换乘站点刷卡数据匹配。

(3)其余站点刷卡数据匹配

其余站点刷卡数据的匹配利用最大相似原理。已通过正交换乘完成上车站点匹配的类为基本类，对应的站点作为控制点；其余的类为待定类，对应的站点作为待定站点。比较待定类与基本类间刷卡时刻间隔和待定站点与控制点间标准行程时间的偏差x_i，以完成其余站点刷卡数据匹配。参照图6，具体描述如下：

对于通过步骤(1)聚类后的数据，令t_ij表示类j的第一次刷卡时刻与类i的最后一次刷卡时刻的差值(类i对应站点处于类j对应站点的上游)，T_uv表示实际中站点u至站点v的标准行程时间。已完成匹配的换乘站点的刷卡记录作为基本类，两相邻基本类(或顶点与基本类)之间的部分称为类段，一个类段所包含的类的数量称为该类段的容量。对于容量为M，对应待定站点数量为N的类段，其第m个类的偏差判断值x_i计算如下：

若为第一个类段，则

x_i＝(t_mR-t_iR)²，i＝m，...，N-m

若为中间类段，则

x_i＝(t_Lm-T_Li)²+(t_mR-T_iR)²，i＝m，...，N-m

若为最后一个类段，则

x_i＝(t_Lm-T_Li)²，i＝m，...，N-m

式中，L(R)表示类m的左(右)临基本类；t_Lm(t_mR)表示待定类m与左(右)临基本类的刷卡时刻间隔；T_Li(T_iR)表示待定站点i与左(右)临控制点间的标准行程时间。

求最小的x_i，其对应的类m和站点i满足最大相似匹配条件，即类m的所有刷卡记录对应的上车站点为i。

5、下车站点匹配

公交乘客在下车站点没有刷卡记录，本发明利用正交换乘和回程进行下车站点匹配。IC卡刷卡数据的下车站点匹配由下车站点匹配器完成，其功能和核心算法描述如下。

(1)利用正交换乘匹配

抽取某日某条线路及其正交换乘线路的所有刷卡数据，提取出在该线路的刷卡时刻早于其正交换乘线路刷卡时刻的对应刷卡记录，判断其刷卡间隔是否小于换乘判断临界时间，如满足则标定该线路的此刷卡记录对应的下车站点为对应的换乘站点，完成换乘站点刷卡数据匹配。

(2)利用回程匹配

抽取某条线路一天的所有刷卡数据，按卡号和刷卡时刻排序，提取出同一张卡在同一条线路有偶数条刷卡数据的记录，第奇数次刷卡记录与第偶数次刷卡记录对应的上车站点互为对方的下车站点。

6、公交站点OD矩阵生成

依据上述步骤匹配得到的刷卡记录对应的上车站点、下车站点，由站点OD生成器生成一日一条线路站点OD矩阵P(DAY，LINE)：

其中，P(DAY，LINE)为日期DAY线路LINE的站点OD矩阵，n为线路站点个数，a_od为站点o到站点d的客流量，其值为经过匹配运算得到的“上车站点为o下车站点为d”的记录的条数。

将得到的矩阵导入Access数据库中，汇总各日各线路的站点OD矩阵，得到全部公交站点OD。

按上述步骤处理常州市公交线路3路2008年3月3日全天(6：30～22：30)的IC卡刷卡数据，共计28700条记录，判断出上车站点的记录19084条，判断出下车站点的记录13784条，上车站点、下车站点均判断出的记录8982条。该线路的站点信息和参数配置信息如表6，推导出的公交OD矩阵如表7。

表6常州市公交3路线路参数与站点编号表

下行	线路编号：003	上行
			10min	车辆末站停留时间	10min
60min	线路行程时间	60min
			3.5min	站点行程时间临界值	3.5min
站点编号	站点名称	站点编号
			003101	花园公交中心站	003201
003102	花园新村	003202

003103	丽丰苑	003203
			003104	马公桥	003204
003105	勤业桥自由动力	003205
			003106	第四医院	003206
003107	江南商场	003207
			003108	新世纪商城	003208
003109	小营前路桥市场	003209
			003110	文化宫	003210
003111	天宁寺	003211
			003112	东太平桥	003212
003113	东坡公园(儿童医院)	003213
			003114	牌楼弄	003214
003115	兆丰花园	003215

003116	高家弄	003216
			003117	政成桥	003217

注：如站点编号“003209”中，“003”表示线路为3路，“2”表示上行，“09”表示编号为09的站点。

Claims

1.一种基于IC卡数据获取公交站点OD的方法，其特征在于先用数据采集模块采集一个城市或城市中某个区域内持有公交IC卡乘客的IC卡刷卡数据，所述采集到的刷卡数据先存储在存储单元中；再由IC卡刷卡数据处理模块调取存储单元中的刷卡数据进行处理，得到公交站点OD信息；最后通过站点OD信息发布模块把公交站点OD信息发布在各终端设备上；

1)刷卡数据按班次聚类：

2.根据权利要求1所述的基于IC卡数据获取公交站点OD的方法，其特征在于：所述IC卡刷卡数据包括：记录编号、IC卡卡号、卡类型、刷卡日期、刷卡时刻、线路编号和车辆编号。

3.根据权利要求1所述的基于IC卡数据获取公交站点OD的方法，其特征是把采集得到的IC卡刷卡数据放入相应的关系数据库；

所述IC卡刷卡数据被分为核心数据信息和参数配置信息，它们分别放在所述数据库的核心数据表和参数配置表中；

核心数据信息包括IC卡刷卡数据和公交站点编号，它们分别放在所述数据库的IC卡刷卡记录表和公交站点编号表中；

参数配置信息包括线路参数、正交换乘参数和相邻站点行程时间临界值，它们分别放在所述数据库的线路参数表、正交换乘参数表和相邻站点行程时间临界值表中。

4.根据权利要求1所述的基于IC卡数据获取公交站点OD的方法，其特征是所述步骤1)前，先对IC卡刷卡数据进行预处理，剔除运行时间范围外的无效数据，并按车辆运行时间先后顺序排列刷卡数据，用来识别不同车次的刷卡记录，步骤包括：

a)抽取某天所有的IC卡刷卡数据；

b)按照车辆编号、刷卡时刻进行排序，得到新编号刷卡数据；

c)记新编号的刷卡数据的序号是r；比较该编号记录的刷卡时刻和前一编号记录的刷卡时刻，得到二者的时间差At；

d)最后判断刷卡记录是否处理完。

5.根据权利要求1所述的基于IC卡数据获取公交站点OD的方法，其特征是所述步骤1)中，

公交车辆从始发站出发沿途停靠每个站点、到达终点站后停留一段时间、折返；与之对应的，抽取某日一条线路的所有刷卡数据，按车辆编号、刷卡时刻顺序排列，以刷卡序列号为X轴、刷卡时刻为Y轴建立二维坐标系，刷卡时刻呈类阶梯状分布，且每隔数个阶梯出现一次较大时间间隔，记为一个GAP，其时间长度为T_GAP，两相邻GAP间的时间间隔记作t_{GAP(i-l)，GAP(i)}，进行如下判断，将不同班次的刷卡记录区分开来，完成班次聚类：

①如果T_GAP(i)≥Depint且t_{GAP(i-1)，GAP(i)}∈[Tratime-Tratime/5，Tratime+Tratime/5]，则GAP(i-l)与GAP(i)之间的刷卡记录作为一个班次的完整刷卡记录；

③如果上述判断均不满足，提取GAP(i-l)后面的所有满足t_{GAP(i-l)，GAP(p)}∈[Tratime-Tratime/10，Tratime+Tratime/10]的刷卡记录，按T_GAP(p)排序，其值最大的GAP(p)与已判断出GAP(i-l)之间的刷卡记录作为一个班次的完整刷卡记录；

其中，Depint是车辆末站停留时间、Tratime是线路行程时间。

6.根据权利要求1所述的基于IC卡数据获取公交站点OD的方法，其特征是所述步骤2)中，

(21)刷卡数据按站点聚类：

公交车辆从始发站出发沿途停靠每个站点，持有公交IC卡乘客在各站点上车刷卡，与之对应的，抽取某日一条线路的一个班次的所有刷卡数据，按刷卡时刻顺序排列，以刷卡序列号为X轴、刷卡时刻为Y轴建立二维坐标系，此时刷卡时刻呈类阶梯状分布，相邻阶梯的时间间隔记为一个STATION，其时间长度为T_STATION，两相邻STATION间的时间间隔记作t_{STATION(j-l)，STATION(j)}；定义一个站点对应的所有刷卡记录为一个类，一个班次包含的类的数量小于等于站点的个数；进行如下判断，完成刷卡数据的站点聚类：

①如果T_STATION(j)≥60s且t_{STATION(j-l)，STATION(j)}∈[Stacriti me-Stacritime/5，Stacritime+Stacritime/5]，则STATION(j-l)与STATION(j)之间的刷卡记录作为一个站点的完整刷卡记录；

②如果T_STATION(j)≥30s且t_STATION(j-l)，STATION(j)∈[Stacritime-Stacritime/10，Stacritime+Stacritime/10]，则STATION(j-l)与STATION(j)之间的刷卡记录作为一个站点的完整刷卡记录；

③如果上述判断均不满足，提取STATION(j-l)后面的所有满足t_{STATION(j-l)，STATION(q)}∈[Stacritime-Stacritime/10，Stacritime+Stacritime/10]的刷卡记录，按T_STATION(q)排序，其值最大的STATION(q)的刷卡记录作为一个站点的完整刷卡记录；

其中，Stacritime是相邻站点行程时间I临界值；

(22)正交换乘站点刷卡数据匹配：

抽取某日某条线路以及可将该线路作为正交换乘线路的其他线路的刷卡数据，提取出在该线路的刷卡时刻晚于其他线路刷卡时刻的对应刷卡记录，判断其刷卡间隔是否小于换乘判断临界时间，如满足则标定该线路的此刷卡记录对应的上车站点为对应换乘站点，完成换乘站点刷卡数据匹配；

(23)其余站点刷卡数据匹配：

其余站点刷卡数据的匹配利用最大相似原理，把已通过正交换乘完成上车站点匹配的类为基本类，以及对应的站点作为控制点，其余的类为待定类，以及对应的站点作为待定站点；比较待定类与基本类之间的刷卡时刻间隔，和待定站点与控制点之间标准行程时间的偏差x_i，以完成其余站点刷卡数据匹配。

7.根据权利要求1所述的基于IC卡数据获取公交站点OD的方法，其特征是所述步骤3)中，

(31)利用正交换乘匹配：

抽取某日某条线路及其正交换乘线路的所有刷卡数据，提取出在该线路的刷卡时刻早于其正交换乘线路刷卡时刻的对应刷卡记录，判断其刷卡间隔是否小于换乘判断临界时间，如满足则标定该线路的此刷卡记录对应的下车站点为对应的换乘站点，完成换乘站点刷卡数据匹配；

(32)利用回程匹配：

8.根据权利要求1所述的基于IC卡数据获取公交站点OD的方法，其特征是所述步骤4)中，先依据匹配得到的刷卡数据对应的上车站点、下车站点，创建一日一条线路站点OD矩阵P(DAY，LINE)：

其中，P(DAY，LINE)为日期DAY线路LINE的站点OD矩阵，n为线路站点个数，a_od为站点o到站点d的客流最，其值为经过匹配运算得到的“上车站点为o下车站点为d”的记录的条数；

再将得到的OD矩阵导入数据库中，汇总各日各线路的站点OD矩阵，得到全部公交站点OD矩阵。

9.一种实现权利要求1～8任一所述方法的系统，其特征在于包括：

(1)数据采集模块：公交车辆上装有IC卡刷卡终端设备；

(2)数据存储模块：数据库服务器；

(3)IC卡刷卡数据处理模块：PC机；

(4)站点OD信息发布模块：信息发布服务器；

各个模块采用有线或无线网络介质进行数据传输。