专利名称：基于网络受限移动对象数据库的交通流数据采集与分析的制作方法
技术领域：
本发明涉及现代交通技术领域中一种网络动态交通流数据的采集以及交通参数的分 析方法。
背景技术：
为了更好地解决交通运输系统所面临的挑战，自20世纪90年代以来，信息技术越来 越多地应用于交通运输部门，其中交通流的分析技术也得到了广泛的研究与应用。为了准 确地获得实时的交通状况信息，人们采取了一系列的技术手段来采集交通数据并获取交通 状况信息，包括使用位置固定的交通流检测设备(如线圈式、微波、光学检测器等)、高 空监视设备，以及浮动车设备等。上述技术手段有着各自的优缺点。其中位置固定的交通流检测设备能够获得交通流 量、平均速度、甚至交通阻塞等信息，但在通常情况下，只能获取设备设置点周围区域的 数据。为了对整个交通网络进行监控，需要设置大量的设备，从而增加了系统的费用。高 空监视设备(如使用直升飞机、卫星等)能够获取交通网络中关于车辆密度、交通拥堵、 平均车流速度等涉及大范围的信息，但同样也存在着费用高昂的缺陷。为了解决上述问题，近年来，越来越多的应用转向了通过浮动车设备来获取交通流数 据，即在整个城市交通网络中， 一些车辆上安装了全球定位系统(GPS)，并通过无线通讯 接口和一定的数据釆集协议，不断地向中心服务器报告自己的位置、速度、行驶方向等原 始信息；服务器对收集到的信息进行分析，得到交通网络的交通状态信息，例如各道路段 的平均通行时间、平均车速、交通阻塞状况等，从而为交通控制与导航服务。这种方法不 受区域和天气条件的限制，并具有经济、灵活的优点，因此得到了人们日益广泛的重视。 然而，现有的基于浮动车法的交通流统计分析技术仍然存在着以下缺陷(1)由于浮动车 系统在进行数据采集时，通常采用的是固定时间间隔或者固定空间距离的方法，当间隔阀 值设得太小时，系统将需要进行频繁的数据通讯，从而带来高昂的通讯代价；反之，当间 隔阀值设得过大时，将造成系统精度的下降。事实上，由于在浮动车系统中，表示和采集 FCD时采用的是直接基于Euclidean空间的方法，没有考虑浮动车与交通路网之间的交互 作用关系，上述矛盾是难以克服的；(2)浮动车系统在进行数据分析时，需要周期性地对 过去某一段时间(如5分钟)之内所采集的浮动车数据进行离线处理，造成了一定程度的 系统延迟，不能实时地得到交通参数；(3)浮动车数据所构成的是Xx YxT空间中大量 离散、无序的点，在与地图进行匹配时，很难精确地还原浮动车的时空轨迹。例如，浮动 车在进行相邻的两次数据采集时，有可能会跨越复杂的道路网络，仅根据两个数据点，系 统无法得到其具体的行驶路线，因此影响了处理的精度。发明内容针对上述现有的三类交通流数据釆集与分析方法所存在的问题与不足，本发明使用网 络受限移动对象数据库技术进行交通流数据的采集与分析，目的是提出一种高效、高精度、 实时的交通流数据的采集与分析方法。为实现上述目的，本发明所采用的技术方案是通过网络受限移动对象数据库技术釆集 和存储移动对象的时空轨迹；然后通过时空统计技术，以原子行驶路段为单位，对移动对 象在运行中所提交的数据进行联机统计，实时地获得交通网络的各项交通参数(如道路网 络中各路段的平均速度、交通流量、交通堵塞状况等)，并保存一段时间周期内交通参数 变化的时态历史数据。为此，需要在每个需要跟踪其位置的移动对象上安装GPS定位设备、实际的交通网络 数据、以及无线通讯设备。其中，上述交通网络数据类似于一张电子地图，可以通过对GIS 地图进行数据转换处理而得到，它包含交通网络中每条道路的信息和每个交叉路口的信息。其中道路的信息包括道路标识、道路的地理几何形状、道路的长度、道路所含原子 行驶路段的集合；交叉路口的信息包括路口的标识、路口的地理位置、路口的区域范围、 路口所连接的道路、路口中交通流的转接关系。道路的标识即道路ID，通常用字母、数字、或者字母与数字的组合来表示，如图l所 示的道路交通网络中共有7条道路，对应的道路ID分别为rl、 r2、 r3、 r4、 r5、 r6、 r7。 每条道路的地理几何形状用一条折线表示，即每条道路都由若干个点来描述，将这些点依 次连接形成的折线即刻画了道路的形状。道路所含的原子行驶路段是道路的某一个行驶方向中(如果道路r是双向道路，则它 有两个行驶方向r+和r-；如果r是单向道路，则它只有一个行驶方向r+或者r-。其中r+ 表示从道路的起点向终点方向行驶，r-表示从道路的终点向起点方向行驶)连接交叉路口 或者道路端点的基本单位，且在其中不再含有移动对象可以驶离该交通流的其它交叉路 口。如图1所示的道路rl是双向道路，因此具有两个交通流^+ (从起点驶向终点)和rl-(从终点驶向起点)。以同一行驶方向所要经过的交叉路口为结点对H+和H-进行分割， 即可得到原子行驶路段。因此rl包含IO个原子路段，且道路双边的原子路段可以不对称。 交叉路口的标识与道路ID的表示方式相同，如jl、j2......，其地理位置用其经纬度坐标表示，其连接的道路用道路ID的集合表示，如在图1中道路rl和r2形成的交叉路口连 接的道路为集合(rl,r2)。另外，用一个数据结构如矩阵(称为交叉路口的连接矩阵)来表 示该交叉路口所连接的各条道路上的交通流之间的转接关系，交通流用道路的ID与正负 号的组合来表示，"+ "表示从道路起点到终点的交通流，而"-"表示从道路终点到起点的 交通流，用不同的数值或符号来表示是否允许移动对象从某一进入交叉路口的交通流(称 为"进入流")切换到某一离开交叉路口的交通流(称为"离开流")。图2给出了一个交 叉路口的连接矩阵的例子。在图2中，道路r8是双向道路，因此它拥有两个交通流，r8+ 和r8-;道路r9是单行道路且只有一个交通流r9-。这三个交通流在连接矩阵中可以有9种 不同的组合，而矩阵元素的值1表示该交叉路口允许移动对象从对应的"进入流"经过该 交叉路口切换到"离开流"，0则表示不允许。数据采集与分析的具体过程包括第一步在交通网络中运行的移动对象通过位置更新操作向中心服务器发送其运行矢 量(包括时间、网络位置、速度和所属交通流方向等信息)。具体的采集方法如下① 移动对象在运行过程中，利用自身的GPS设备不断地测得其最新的运行参数，包括 当前时间t、当前经纬度位置(X， Y)、当前速度v和当前运行方向d等信息。② 将上述GPS测得的位置及方向数据与实际的交通道路网络进行匹配，得到移动对象 在交通网络中的相对数据，即网络位置(rid， pos)和道路交通流方向fd，其中，rid为移 动对象所在的道路的标识；pos是移动对象在道路rid中的相对位置。设每条道路的长度为 1，则道路rid中的任意一个位置可用一个实数pos表示，0￡posSl; fd为移动对象在道 路rid中的运行方向，用"+ "表示由道路rid的起点向终点运行，"-"表示由终点向起点 运行，即fde(+，-)。具体的转换方法如下*计算对应的rid:由于每条道路都是由一条折线来刻画的，通过计算(X， Y)坐 标点位于哪条折线上，即可得出对应的道路ID，即rid值；*计算对应的POS:对于交通网络中每一条道路的两个端点，可以设定X坐标值小的一端为起点，另外一端为终点；当两个端点的X坐标相等时，设定Y坐标值小的一端为起点，另外一端为终点。定义起点的pos值为O，终点的pos值为1， (X，Y)坐标点到起点沿道路的距离与该道路的长度的比值即为该坐标点对应的pos值，显然地，pos的取值范围是0和1之间，即posE[O,l]。 计算移动对象所属的交通流方向fd:比较GPS所测得的移动对象运行方向和道路
rid在pos处的走向，如果运行方向是朝向道路终点的，则fd为"+";若运行方向是朝向道路起点的，则fd为"-"。 经过上述转换后，得到移动对象在该时刻的运行矢量，形如(t,(rid,pos),v， fd)。 ③移动对象将上述转换后的运行矢量与其最后一次位置更新时发送给服务器的运行 矢量相比，考察两个运行矢量是否满足一定的位置更新条件。如果满足以下任何一种条件，就需要向中心服务器发送当前运行矢量(即t，(rid,pOS),V， fd);否则不向中心服务器发送当前运行矢量，并转步骤①进行下一轮处理。移动对象将始终保存最后一次位置更新时的运行矢量数据。具体来说，本发明考察以下三种运行矢量的变化(1) rid发生了变化即该移动对象当前所在的道路不同于上一次位置更新时所在的 道路，即移动对象已经从一条道路运行到另外一条道路；(2) rid没有变化，但是移动对象的运行方向fd发生了变化，或者运行速度v发生了 较大的变化即该移动对象的运行方向由正方向变为负方向或由负方向变为正方向；或者 运行方向没变，但是该移动对象的当前速度与上一次位置更新时的速度相比，其偏差超过 了一个指定的阀值^ (例如，10公里/小时)；(3) rid和v均没有发生上面描述的变化，但是该移动对象的当前位置与预测位置发 生了变化即移动对象的当前位置(记为(rid, posn。w))与根据上次位置更新信息预测的 位置(记为(rid， posevl))相比，距离偏差超过了一个指定的阀值《(例如500米)。其中，posevl的计算方法如下设该移动对象上次位置更新时的运行矢量为(t。， (rid。， pos。),v。,fd。)，道路rid。的长度为len，当前时间为tnow，则pOSevi = pOS。
+fd X -len根据P0S^和p0Sn。w,可以得到预测位置与实际位置之间的距离偏差为|poseV| - pos,l X len如果该值大于《则触发一次位置更新。在判断上述三种位置更新条件时，需要对移动对象的当前位置处于交叉路口区域的特殊情况进行处理。因为在实际的交通道路网络中，交叉路口是一个区域范围(可以在交通网络数据中为每个交叉路口指定一个区域半径"与交叉路口中心的距离Sy的区域都属于该路口的区域范围)，在该区域内移动对象的运行矢量可能发生非常频繁的变化，但实际的位置偏移却并不大，所以在交叉路口区域中没有必要将所有的位置更新信息都发送给服务器。具体措施是a)当移动对象处于交叉路口区域中时，视为其仍然在原来的道路土行
驶，从而保证了移动对象在交叉路口区域中不发生换路位置更新(见上述判断条件(l)); b)当移动对象位于交叉口区域中时，如果其运行矢量发生了以上(2)、 (3)中描述的任意 一种变化，则首先检查该移动对象在当前交叉路口区域范围内是否已经向服务器发送过一 次位置更新信息，如果已经发送过，则不再发送当前的运行矢量，否则将向服务器发送最 新的运行矢量。上述方法保证了移动对象在交叉口区域内最多只发送一次位置更新信息。上述移动对象运行矢量数据采集方法是基于"惯性原理"实现的，即系统假设移动对 象将沿着当前的道路、以大致稳定的速度继续行进一段时间。 一旦该假设条件不成立则向 中心服务器报告一次最新的运行矢量数据。根据上述方法，当中心服务器没有得到一个移动对象的最新位置更新信息时，则可以 认为其仍然在原有道路上近似保持原有的速度运行，且其当前位置可以通过最近一次的位 置更新信息(包括位置、速度、方向和时间)近似地估算出来，而预测偏差在指定的允许 范围内。因此从位置跟踪精度来看，上述方法可以从下面几方面保证所釆集的数据的准确 性a)移动对象在任一时刻所在的道路能够通过所采集的数据唯一确定；b)移动对象在 任意时刻的位置可以通过所釆集的数据进行精确计算，且计算误差必然小于设定的距离阀 值；C)移动对象在任意时刻的速度可以通过该时刻之前采集的运行矢量获得，且误差必然 小于设定的速度阀值。此外，由于本发明提出的位置更新方法考虑了移动对象和道路网络之间的关系，所以 移动对象不需要周期性地报告位置信息，中心服务器也能够通过预测技术较准确地确定出 移动对象在任意时间点的位置，从而在保证数据釆集精度的前提下，减少了数据采集频率， 降低了通信代价。第二步中心服务器接收到移动对象发送来的采集数据之后，将相应的运行矢量(t, (rid, pos), V， fd)附加到对应的移动对象时空轨迹中，从而得到该移动对象的时空轨迹。移动 对象mid的时空轨迹是由该移动对象提交的一组运行矢量按照时间顺序排列所组成的，即((ti，(ridi,posi)，vi， fdi))L (其中，n为正整数，下同)，它可以描述移动对象过去及当前时刻连续的动态位置变化过程。中心服务器在存放移动对象的时空轨迹时是以轨迹片段为单位分片存放的，与各条道 路相对应的轨迹片段存放在相应的道路中，它们将作为统计源数据。轨迹片段是时空轨迹 的一部分，对应于移动对象在某条道路上的行驶过程。移动对象在道路rid上的时空轨迹片段形式为(mid， fd, entry, isactive, (ti, poSi， v，),)，
其中，mid是移动对象的标识，fd是移动对象在道路rid中所属的交通流方向，entry是该 移动对象以同样的方向进入该道路的序号。移动对象可能以相同的方向多次进入同一条道 路，每一次将产生不同的时空轨迹片段，这些片段以entry编号加以区分。isactive标记当
前时刻该时空轨迹片段是否仍处于活动状态。(ti, posi, v,):是移动对象的时空轨迹片段所对 应的时空轨迹曲线，其中ti是一个时间点，posi是ti时刻移动对象在道路rid中的相对位置，
V,是移动对象在tj时刻的速度。(t,, (rid, p0Si)， Vi， fd)就是移动对象在ti时刻的运行矢量，移
动对象所提交的最后一个运行矢量称为活动运行矢量。当一个活动移动对象在路网中行进 时，可能会在不同的道路产生多个时空轨迹片段，这些轨迹片段中只有最后一个处于活动 状态。
从语意上讲，移动对象mid在道路rid上的时空轨迹片段(mid, fd, entry, isactive, (ti, poSi，
Vi);L,)刻画的是该移动对象在道路rid中的连续运行过程，其中任意两次连续的位置更新消
息之间的状态可以通过插值的方法计算出来。此外，如果isactive为真，则表明移动对象 仍然在道路rid中行进，其当前位置可以通过tn、 fd、 posn、 Vn等计算出来。因此我们将时 空轨迹片段上的任意一点(mid， fd, entry, isactive, (t, pos, v))称为时空轨迹点。时空轨迹点反 映移动对象在t时刻瞬间的运行状态。
第三步每次当服务器接收到移动对象发送的位置更新消息时，触发相应道路的交通 参数刷新过程，实时地刷新该道路中各原子行驶路段的交通参数。
在本发明中，每个原子行驶路段与五个基本交通参数相联系，包括移动对象个数 (vmo)、平均点速度(vp)、平均通行时间(Tt)、移动对象流量("mo)、以及交通拥堵 状态(jammed)。
根据统计方式的不同，上述五个基本交通参数可以划分为两大类(1)瞬间交通参数 (包括VIIIO和vp),这类参数反映的是原子行驶路段中的各个移动对象在当前时刻的统计 值，即在进行统计计算时，只有当移动对象目前正处于该原子行驶路段中时，才需要参加
瞬间参数的刷新；(2)时段交通参数(包括Tt、 "mo、 jammed),这类参数值的计算过程 需要对一个时间段(即当前时刻之前的某个时间窗口)内的移动对象时空状态进行统计， 即在统计时，即使该移动对象目前已经驶离了该原子行驶路段，只要其时空轨迹与统计窗 口相交，就需要参加相应的统计计算。
当中心服务器接收到移动对象发送的位置更新消息时，将实时地触发交通参数的刷 新。由于每一条道路中均存放有相关移动对象的轨迹片段，因此参数刷新的过程即是对这些轨迹片段进行统计的过程。
对于瞬间交通参数而言，仅需要对处于活动状态(即isactive为真)的轨迹片段进行 统计，根据其最后一个运行矢量，计算出移动对象的当前位置以及速度，从而刷新相应原 子行驶路段的移动对象个数wmo和平均点速度vp参数值。
时段交通参数的刷新需要统计当前时刻t'，^之前某一段时间〖timeseg之内的轨迹片 段，所有与时间窗口[tn。广5timeseg， t,]相交的时空轨迹片段均需要参加统计计算。由 于轨迹片段包含有相应的移动对象在该道路中行驶时所具体通过的原子行驶路段以及各 路段的通行时间，因此可以方便地计算出对应路段的平均通行时间(Tt)和移动对象流量 ("mo)。另外，轨迹片段中还包含有移动对象的速度信息，因此可以计算出道路各原子 行驶路段的拥堵状态。
在对某个道路r的各原子行驶路段进行交通参数统计的过程中，投影是一个重要的操 作，用n表示。投影操作具体分为以下几种类型(在下面的表示中，TSEGr表示道路r上所 有轨迹片段的集合)
* rUs表示在某个原子行驶路段上的投影，例如，rU(TSEG" ar&)表示计算TSEGr 中所有轨迹片段落在原子行驶路段arsi内的部分，其结果仍然是一个轨迹片段的集合，即 包含所有落在arsi区间内的轨迹片段。在针对某个具体的轨迹片段进行上述投影计算的过 程中，需要对该轨迹片段进行截取操作，即只取其中属于该原子行驶路段的部分；
* EL表示在当前时刻tn。w之前某个时间段上的投影，例如，rU(TSEGr, ltimeseg) 表示计算TSEGr中所有轨迹片段落在时间窗口[t,- ltimeseg， t,]内的部分，其结果仍 然是一个轨迹片段的集合，即包含所有落在该时间区间内的轨迹片段；
* n,表示在某个网络位置(pos)处的投影，例如，n,s(TSEGr， pos)表示计算TSEGr
中所有轨迹片段落在pos处的部分，其结果是一个时空轨迹点的集合，即包含所有轨迹片 段与网络位置pos的交点；
* n'表示在某个时间点上的投影，例如，rL(TSEG"t)表示计算TSEGr中所有轨迹片 段落在时间点t上的值，其结果是一个轨迹点的集合，即包含所有轨迹片段与时间点t的 交点。
另外，求和也是一个重要的操作，用E表示。求和操作具体分为以下几种类型-
* E,表示求集合中元素个数的总和。例如，E (nt(TSEGr, t))表示求集合 rL(TSEG" t)中轨迹片段的个数。由于一个移动对象在任意一个时间点(例如t)上只可 能在道路r上存在最多一个轨迹片段，所以E,(TL(TSEG" t))的结果也就是在时刻t道 路r上的移动对象个数。
* E^表示求集合中点速度的总和。例如，Evp(nt(TSEG" t))表示在时刻t，道路r 上的所有移动对象点速度的总和。由于rit(TSEGr， t)的结果是t时刻各时空轨迹点的集合， 其中的速度即为各移动对象在t时刻的瞬间速度，对它们求和即可得到点速度的总和；
Eu表示求集合中通行时间的总和。例如，UrLrs(TSEG" arsi))表示求集合n ars(TSEGr， arsi)中所有元素所包含的通行时间总和。如前所述，集合rU(TSEGr， arSi)中 的每个元素是一个轨迹片段，该轨迹片段的起始时间和结束时间分别为相应的移动对象进 入和离开原子行驶路段arsi的时间，而两个时间的差即为通行时间，而Eu操作则是求这 些通行时间的总和。
在上述两个操作定义的基础上，下面我们具体给出道路r的各交通参数的统计计算步 骤(设道路r中的所有轨迹片段集合为TSEGr，道路r含有n个原子行驶路段，分别用arSi (1 Si ^n)表示)
(1) 计算在当前时刻t^，原子行驶路段arsi中的移动对象个数(vmo)-
Vino = E,(nt(nars(TSEGr, arsi),t,)) (公式1)
由于rUrU(TSEGr， arSi) ， t,)的结果为在当前时刻t,,原子行驶路段ars!上的所有 轨迹点的集合，根据任意一个移动对象在任意一个时间点只可能在ars,上存在最多一个轨 迹点，可知E,(rL(rU(TSEGr, arsi)， t,,,J)的结果即为在时刻t,原子行驶路段ars;上的 移动对象的个数。
(2) 计算在当前时刻t目，原子行驶路段arsi中的平均点速度(vp):
vp = ￡vp(nt(rL,rs(TSEGr, arsi)， t』/ n (公式2) 其中，n= E um(nt(n,(TSEGr,ar&)，tJ)，即满足计算条件的移动对象的个数。 由于Zw(TL(n,(TSEGr， ars,), t,))表示在当前时刻t,，原子行驶路段ars;上的所有
移动对象的点速度之和，所以它与相应移动对象个数的比值即为在当前时刻原子行驶路段
a,rs,上的平均点速度。
(3) 计算原子行驶路段ars,在给定时间窗口[U,- ^ timeseg， t 。,]内的平均移动对象流 量(como):
"mo= i; u (n P,,s(nts(nars(TSEGr, ars》，Stimeseg), pos)) /《timeseg (公式3) 其中，pos为arsi的终点位置。
nts(n (TSEGr， arSi), I timeseg)表示在时间窗口 [t, - l timeseg， tn。 ]内在ars!行驶的移动对象轨迹片段的集合，n,(rU(rU(TSEGr,arSi),〖timeseg), pos)则表示其 中穿越了arsi终点的部分，结果是轨迹点的集合。对其中的元素计数，即可得到在给定时 间内穿越arsi的移动对象个数，进而除以〖timeseg即可计算出单位时间内通过arsi的移 动对象个数，即移动对象流量。
(4) 计算原子行驶路段arsi在给定时间窗口[t^-〖timeseg， t 。 ]内移动对象的平均通 行时间(Tt):
Tt = E"(nt,(rUs(TSEG"arSi)， ^timeseg)) / n (公式4) 其中，n= H(UTSEG"arSi), "iraeseg))，即满足计算条件的移动对象的 个数。
由于E"(n's(n腦(TSEG" arsi)， S timeseg))表示在时间段[t, - "imeseg， t,]内， 原子行驶路段arsi上的所有移动对象的通行时间之和，所以它与相应移动对象个数的比值 即为在时间段ltimeseg内，原子行驶路段arSi上的平均通行时间。
(5) 计算给定时间窗口 ^ti鹏seg内，原子行驶路段arsi上的交通拥堵状态(jammed):
首先计算arSi上是否存在拥堵区域jaramedarea:
j柳medarea = nsl' (ILs(n,(TSEGr, ars》,《timeseg), vslOT) (公式5) 其中，n&表示求轨迹片段集合中的缓慢轨迹片段(即运行速度低于指定速度阈值
v,^(如10公里/小时)的轨迹片段)所对应的地理区域的交集。
然后，根据ja,edarea是否为空判断ars是否处于拥堵状态，即
jammed = (j柳medarea # NU山 (公式6)
由于交通堵塞是一种集体行为，单个移动对象的缓慢行驶并不能代表道路发生了阻
塞，因此需要进行特殊的处理。方法是通过移动对象所提交的轨迹片段，计算出其行驶缓
慢的路段。如果道路中存在着某个区域，所有的移动对象在过去的一段时间内经过该区域
时均处于缓慢行驶状态，则可以判断该路段为拥堵路段。
第四步将新计算出的交通参数值保存在数据库相应的道路记录中供用户查询。
本发明使用的基于网络受限移动对象数据库的交通流数据采集及分析技术能够高效
准确地获取移动对象完整的时空轨迹，它不仅简化了统计分析的复杂性，提高了处理效率，
而且为数据分析提供了更为丰富的信息，提高了数据处理的精度，有效地加强了交通流分
析的广度和深度。与现有的交通流数据采集及分析方法(如浮动车法)相比，本发明具有 如下优势
(1) 在保证精度的前提下提高了数据采集的效率。由于数据采集是通过移动对象的 位置更新操作完成的，移动对象只有在运行矢量发生变化时才需要触发位置更新过程，因 此当它沿着某条道路平稳行驶时，可以长时间不进行位置更新，从而降低了数据采集的通 讯及计算代价，提高了效率；
(2) 本发明所采集的数据是移动对象的精确时空轨迹数据，其中的运行矢量能够精
确地刻画移动对象的实际运行状态。我们以移动对象提交的最后一个运行矢量(称之为"活
动运行矢量")为例，根据本发明中的位置更新协议能够保证a)移动对象实际所处的道 路必然与活动运行矢量一致；b)移动对象的实际位置与活动运行矢量的计算位置之间的误 差必然在规定的范围之内；C)移动对象的实际速度与活动运行矢量中速度的误差必然在规 定的范围之内；
(3) 本发明中的道路交通流参数统计分析方法是基于时空轨迹进行的，由于时空轨 迹刻画的是移动对象的整个连续的行驶过程，因此克服了其它方法需要根据离散点进行推 断的缺陷(例如浮动车法需要根据两个数据釆样点推断其实际的行车路线，由于两点之间 可以有多条可能的网络连线，从而不可避免地会带来误差)，提高了数据处理的精度；
(4) 本发明的道路交通参数刷新过程是通过位置更新实时触发的，从而消除了统计 参数的信息延迟。
综合以上分析，本发明可以保证交通流统计分析的高效性、高准确性和实时性。该发 明的最终结果可以直接提供给相关的用户使用，例如交通控制中心、特种车辆司机或动态 导航设备用户等，也可用于支持更为复杂的査询，如某个区域或某条道路的平均车速、平 均车流量，以及考虑了当前平均通行时间的动态最短路径等，从而进一步为交通控制及交 通导航服务。
图1是道路、原子行驶路段及交叉路口的示意图。
图2中，(a)是某一交叉路口的交通流转接关系示意图，(b)是对应的连接矩阵。 图3是本发明统计交通流方法的流程图。
图4示意了本发明实施例中移动对象mol某时间段在道路d上的的运行过程及服务器 根据收到的位置更新消息生成的对应的时空轨迹片段。
图5是本发明实施例中移动对象的缓慢路段及拥堵检测的示意图。
具体实施例方式
下面结合附图，通过实例进一步说明本发明，但不以任何方式限制本发明的范围。 根据图2所示的流程进行如下具体操作
1. 利用网络受限移动对象数据库技术采集移动对象的位置更新信息。
例如移动对象mol在某个时间段的运行过程如图4所示，其在时刻t4时发送的位 置更新信息为(t4,(rl，0.12),v4)，接着是(t5,(rl,0.25),v5)， (t6,(rl,0.52),v6)，最新的位置更新信 息为(t7,(rl,0.73)，v7),经过的原子行驶路段依次为arsl、 ars2、 ars3，目前正在ars4上运行。
2. 生成对应的轨迹片段。
如移动对象mol在道路rl上的轨迹片段如图4所示，其中最后一项(t7，0.73,v7)对应于
活动运行矢量。
3. 触发交通参数的刷新。
图4中所示的轨迹片段是道路rl上的轨迹片段之一。交通参数的刷新方法如下
(1) 在计算瞬间变量Vmo和vp时，根据移动对象mol在t7时刻提交的活动运行矢 量，可以计算出mol目前正处于原子行驶路段ars4中的某一位置，因此mol需要参加ars4 的瞬间参数的统计计算，而不需要参加其它原子行驶路段(如arsl、 ars2、 ars3)瞬间参数 的计算。同样地，如果得知移动对象mo2、 mo3目前都处于ars4中，且其活动运行矢量分 别为(ta,(rl,0.85),v8)、 (tb,(rl，0.84),v9)。以原子行驶路段ars4为例，
计算当前移动对象的个数vmo: Vmo-3，即共有mol、 rao2、 mo3三个移动对象。
计算当前的平均点速度vp=(v7+v8+v9)/3，即mol、 mo2、 mo3当前速度的平均 值。
(2) 在计算时段参数"mo和Tt时，设给定的时间窗口[t,- ^timeseg, t哪]为从 时刻t4到tn。w (当前时刻)，则如图4所示的移动对象mol需要参加原子行驶路段arsl、 ars2、.ars3、 ars4的时段参数的计算。以原子行驶路段ars4为例
计算在当前时刻之前的5timeseg (例如5分钟)时间内移动对象的流量"mo:
首先通过公式n"(n ars(TSEGr， ars4), ^ timeseg)计算在时间窗口 [t丽-
^timeseg， t,]内经过了 ars4的所有轨迹片段，然后求这些轨迹片段中经过了
ars4终点(即跨越了ars4)的轨迹片段，对之进行计数，得到满足条件的移动对
象数目，最后除以ltimeseg即可得到所需结果。
计算在当前时刻之前的Stimeseg (例如10分钟)内的平均通行时间计算方法 与"mo相似。在计算出时间窗口 [t 。》 - ^ timeseg， t 。 ]内经过了 ars4的所有轨 迹片段之后，根据轨迹片段的开始和结束时间之间的差值，得到每个移动对象通 过ars4的时间，取它们的平均值，即可得到ars4的平均通行时间Tt. G)计算原子行驶路段arsj的交通拥堵路段状态的方法如图5所示在当前时刻之前 的指定时间范围5 timeseg (例如8分钟)内，有三个移动对象mol、 mo2、 mo3经过了 aw,且它们的运行速度低于指定速度阈值(如10公里/小时)的路段如图所示(见"缓慢 行驶路段")，即:
mol的缓慢行驶路段是从posl到pos2、从pos4到pos6。 mo2的缓慢行驶路段是从pos2到pos5、从pos7到pos8; mo3的缓慢行驶路段是从pos3到pos5、从pos7到pos8;
这些缓慢行驶路段的交集为从pos4到pos5之间的路段，即所有移动对象都在该路段 上都行驶缓慢，所以得知，原子路段ars,的交通拥堵路段(即jammedarea)为从pos4到 pos5的部分，由于jammedarea不为空，因此其拥堵状态jammed值为"真"。
4.将交通参数变化的时态数据保存在动态交通网络相应的道路中，供查询处理使用。
以上通过实施例对本发明进行了详细的描述，本领域的技术人员应当理解，在不超出 本发明的精神和实质的范围内，对本发明做出一定的修改和变动，比如交通参数的刷新过 程可以用定时触发来代替通过移动对象的位置更新实时地进行触发，又如用其他统计方法 代替本发明所用的公式来计算交通参数，仍然可以实现本发明的目的。
权利要求
1.一种基于移动对象数据库的交通流数据采集与分析方法，包括以下步骤1)数据库服务器不断地采集交通网络中各移动对象的运行矢量(t，(rid，pos)，v，fd)，其中t是该运行矢量所对应的时间；(rid，pos)是移动对象在t时刻的位置，其中rid是移动对象所在道路的标识，pos是移动对象在道路rid中的相对位置；v是移动对象在t时刻的运行速度；fd是移动对象所属的交通流方向；2)数据库服务器根据采样数据生成该移动对象在道路rid中的时空轨迹片段(mid，fd，entry，isactive，(ti，posi，vi)i＝1n)，并将其保存在数据库中对应的道路记录中，其中mid是移动对象的标识；fd是该移动对象在道路rid中所属的交通流方向；entry是该移动对象以同样的方向进入该道路的序号；isactive标记当前时刻该时空轨迹片段是否处于活动状态；(ti，posi，vi)i＝1n是移动对象时空轨迹片段所对应的时空轨迹曲线，其中ti是一个时间点，posi是ti时刻移动对象在道路rid中的相对位置，vi是移动对象在ti时刻的速度，n为正整数；3)以原子行驶路段为单位，对所采集的轨迹片段数据进行联机统计分析，实时地获得交通网络的各项交通参数，所述原子行驶路段是道路的某一个行驶方向中连接交叉路口或者道路端点的基本单位，且在其中不再含有移动对象可以驶离该交通流的其它交叉路口；4)将交通参数变化的时态数据保存在数据库相应的道路记录中。
2. 根据权利要求l所述的交通流数据采集与分析方法，其特征在于所述步骤l)采集移动对象运行矢量的方法是，首先在每个移动对象上安装GPS定位设备和无线通讯设备，并配备实际的交通道路网络数据，所述的交通道路网络数据包 含交通网络中每条道路的标识、地理几何形状和长度，以及每个交叉路口的标 识、地理位置、路口的区域范围、路口所连接的道路、路口中交通流的转接关系；移动对象个体在运行过程中，依据下列步骤向数据库服务器报告自身位置 更新信息-1-1.利用GPS设备不断地测得当前运行参数，包括当前时间t、当前经纬度位置(X， Y)、当前速度v及当前方向； 1-2.将GPS测得的位置及方向数据与实际的交通道路网络数据进行匹配，得到 移动对象在交通网络中的相对数据，即网络位置(rid， pos)和道路交通流方向fd,从而得到移动对象当前的运行矢量(t, (rid, pos), v， fd)，其中rid 为移动对象所在道路的标识；pos是移动对象在道路rid中的相对位置，fd 为移动对象在道路rid中所属交通流方向； 1-3.将经过步骤1-2转换的当前运行矢量与其上次位置更新时发送给服务器的运行矢量相比，首先判断是否发生了下列三种情况之一a. rid改变；b. rid没变，但fd改变或者v与上次位置更新时的速度相比，其偏差超过 了一个指定的阀值；c. rid、 fd和v没变，但当前位置(rid， posn。w)与根据上次位置更新时的 运行矢量预测的位置(rid， posevl)相比，距离偏差超过了一个指定的 阀值；如果上述三种情况均未发生，则不向服务器发送当前运行矢量；如果发生 了上述情况a，则直接向服务器发送当前运行矢量；如果发生了上述情况b 或c，则再判断是否当前位置处于交叉路口的区域范围内且上次位置更新也 在该区域内，是则不向服务器发送当前运行矢量，否则发送。
3.根据权利要求2所述的交通流数据采集与分析方法，其特征在于在所述的交通道路网络数据中设定每条道路的两个端点中一个为起点，另一个为终点，定义起点的pos值为0，终点的pos值为1，则在所述步骤1-2中首先通过计算(X， Y)坐标点位于哪条道路上，得到rid值，然后计算(X， Y)坐标点到起点的距 离与该道路的长度的比值得到该坐标点对应的pos值，pose[O,l];移动对象所 属的交通流方向fde(+，-}，其中"+"表示由道路rid的起点向终点运行，"一" 表示由终点向起点运行。
4.根据权利要求2所述的交通流数据采集与分析方法，其特征在于在所述c情 况下，设上次位置更新时的运行矢量为(t。,(rid。，pos。)，v。，fd。)，道路rid。的长度 为len，当前时间为Ww，贝U:<formula>formula see original document page 3</formula>根据p0Sev!和p0Sn。w，可以得到预测位置与实际位置之间的距离偏差为 |poseVi - pos瞎l X len 。
5. 根据权利要求2所述的交通流数据釆集与分析方法，其特征在于服务器一旦接收到移动对象发送来的运行矢量就触发所述步骤3)对交通参数进行实时刷新， 对于瞬间交通参数，根据时空轨迹片段直接查找该移动对象所在的原子行驶路 段，重新计算其相应的参数值；而对于时段交通参数，査找该移动对象与统计窗 口相交的时空轨迹片段，重新计算相应原子行驶路段的参数值。
6. 根据权利要求1 5中任一权利要求所述的交通流数据采集与分析方法，其特征 在于所述步骤3)的交通参数包括在当前时刻t,某原子行驶路段ars,中的移 动对象个数vmo，通过下列计算公式得到vmo = E,(rL(n肌(TSEG" arSi)， t，)) 式中，TSEG,表示原子行驶路段arSi所属道路r上所有轨迹片段的集合，TLb表 示路段投影，rL表示时间点投影，E,表示求集合中元素个数的总和。
7. 根据权利要求1 5中任一权利要求所述的交通流数据采集与分析方法，其特征 在于所述步骤3)的交通参数包括在时刻t,某原子行驶路段arsi中的平均点 速度vp，通过下列计算公式得到<formula>formula see original document page 4</formula> n 式中，n= E,(nLs(TSEGr,ars丄t,)),其中TSEGr表示原子行驶路段arSi所属道路r上所有轨迹片段的集合，n,表示路段投影，rL表示时间点投影，E,表示求集合中元素个数的总和，Evp表示求集合中点速度的总和。
8. 根据权利要求1 5中任一权利要求所述的交通流数据采集与分析方法，其特征 在于所述步骤3)的交通参数包括在给定时间窗口[t,- ^timeseg， t,]内 某原子行驶路段arSi的平均移动对象流量como，通过下列公式计算得到<formula>formula see original document page 4</formula>式中，pos为arsi的终点位置，TSEGr表示原子行驶路段arSi所属道路r上所有轨迹片段的集合，n,表示路段投影，rits表示时间段投影，n,表示位置点投影，E,表示求集合中元素个数的总和。
9. 根据权利要求1 5中任一权利要求所述的交通流数据采集与分析方法，其特征 在于所述步骤3)的交通参数包括在给定时间窗口[t,-〖timeseg， t,]内 某原子行驶路段arSi中移动对象的平均通行时间Tt，通过下列公式计算得到Tt = Hu(nts(n rs(TSEGr, arsj， ltimeseg)) / n 式中，n= H(n,(TSEG"arSi)， Stimeseg))，其中TSEGr表示原子行驶 路段arsi所属道路r上所有轨迹片段的集合，n,表示路段投影，EU表示时间 段投影，E,表示求集合中元素个数的总和，E"表示求集合中通行时间的总和。
10. 根据权利要求1 5中任一权利要求所述的交通流数据采集与分析方法，其特 征在于所述歩骤3)的交通参数包括在给定时间窗口[t。。广UiiTieseg， t,] 内原子行驶路段arsi上的交通拥堵路段j柳medarea，通过下列公式计算得到jaramedarea =门si冊(riu(n,(TSEGr, arsO, l timeseg) , vsk> ) 式中，TSEGr表示原子行驶路段ars,所属道路r上所有轨迹片段的集合，nars表示路段投影，n"表示时间段投影，n^表示求轨迹片段集合中运行速度低于指定速度阈值v.一的地理区域的交集； 根据jammedarea是否为空判断arsi是否处于拥堵状态，艮P: jammed 二 (jammedarea # NULL) 。
全文摘要
本发明提供了一种基于网络受限移动对象数据库的交通流数据采集及分析方法，通过网络受限移动对象数据库技术采集和存储移动对象的时空轨迹，并通过时空统计与分析技术，以原子行驶路段为单位，对所采集的数据进行联机统计，实时地获得交通网络的各项交通参数，同时保存交通参数变化的时态数据。这种方法所采集的数据能够反映移动对象完整的时空轨迹，它不仅简化了统计分析的复杂性，提高了处理效率，而且为数据分析提供了更为丰富的信息，提高了数据处理的精度。另外，交通参数刷新过程可通过移动对象的位置更新实时地进行触发，因此消除了统计参数的信息延迟。与现有的交通流统计分析方法相比，该方法具有高效、高精度、实时的优点。
文档编号H04Q7/38GK101127159SQ20071012193
公开日2008年2月20日 申请日期2007年9月18日 优先权日2007年9月18日
发明者丁治明 申请人:中国科学院软件研究所