CN106095802B - 城市路网上基于r树的全时态移动对象索引方法 - Google Patents

Abstract

本发明公开了城市路网上基于R树的全时态移动对象索引方法，首先根据移动对象位置更新信息m查找基于x‑y的上层索引，找到移动对象所处的封闭路段e，并取出其相关信息；然后形成运行矢量mv＝(mid,t,v,rid,pos)，通过rid访问道路索引表，取得道路对应的R‑tree；根据mid访问轨迹单元哈希表，获得轨迹单元链表，访问链表尾结点lu‑tail；根据lu‑tail是否为空，生成链表轨迹单元LU；将连续轨迹单元和活动轨迹单元相继添加到链表尾部，建立索引，同时维护链表信息。本发明引入了动态链表和移动对象哈希表，可直接得出移动对象的历史轨迹信息，无需搜索整个层移动对象树森林。采用了指数平滑法预测移动对象的速度，提高了速度预测的精度，支持移动对象的全时态位置查询。

Description

城市路网上基于R树的全时态移动对象索引方法

技术领域

本发明属于移动对象数据管理技术领域，尤其涉及城市路网上基于R树的全时态移动对象索引方法。

背景技术

随着无线通信与GPS等空间定位技术的发展，通过对移动终端的位置进行追踪管理以提供相关查询服务变得越来越普及。现实生活中人们真实的运动场景往往被限制在一定固定区域内，比如常见的城市路网，因此研究路网中移动对象索引结构具有现实意义。路网中存在规模巨大的移动对象，利用单层的R-tree或者R*-tree对移动对象进行索引会忽略了静态路网自己的特点，得不到很好的查询效果。针对R-tree和R*-tree不能很好地索引移动对象的情况，Saltenis等人基于R*-tree提出了TPR-tree移动对象索引结构，TPR-tree中的最小边界矩形(minimum bounding rectangle，MBR)含有速度等参数信息，能够随着时间参数变化，但是，由于没有存放历史轨迹，TPR-tree及其变种仅能支持移动对象当前以及短暂将来位置的查询，不支持对过去位置信息的查询。

目前基于路网的移动对象位置索引结构的研究还处于起始阶段，大部分结构只能满足部分需求。Frentzos等人提出双层索引结构FNR-tree，由一个2D R-tree和一个1DR-tree森林构成，由于仅存储进入和离开边的数据，导致窗口范围查询结果的准确性比较低，而且只能索引移动对象的历史轨迹，而无法索引移动对象的实时位置。MON-tree是对FNR-tree进行改进而提出来的索引结构，其支持对移动对象的窗口查询和历史轨迹信息查询，但缺乏对时间信息的管理和网络拓扑空间持续性优化，所以时间片查询效率低。丁治明等人提出的NDTR-tree索引结构不仅可以可以快速地出路对移动对象过去可能位置的查询，而且能够对其当前和未来的可能位置进行高效的查询处理。但其在移动对象轨迹查询方面效率不高。

公开号为CN103235831A、名称为"基于路网的移动对象位置索引结构及索引方法"的发明专利提供了一种基于路网的移动对象位置索引结构和索引方法，其优点是在下层R-tree索引中增加了静态对象R-tree，能够很好地支持最近邻查询，其缺点是采用了R-tree构建静态路网索引，查询效率低、构建索引时节点分裂次数多。

发明内容

本发明的目的在于提供一种城市路网上基于R树的索引技术，需要解决的问题是支持对城市路网上移动对象的全时态索引，其次，实现对移动对象全时态的预测，具有良好的查询精度。

为实现以上目的，本发明提出城市路网上基于R树的全时态移动对象索引方法，包括以下步骤：

步骤1：根据移动对象位置更新信息m查找基于x-y的上层索引，找到移动对象所处的封闭路段e，并取出其相关信息；

步骤2：通过e信息以及x，y坐标信息求得移动对象在所属道路上的pos位置，形成运行矢量mv＝(mid,t,v,rid,pos)，其中，pos∈[0，1]，为对象所处道路的相对位置，t为采集该运行矢量的时间，v表示该移动对象在t时刻的速度，rid为道路标识；

步骤3：通过rid访问道路索引表，取得道路对应的R-tree；

步骤4：根据mid访问轨迹单元哈希表，获得轨迹单元链表，访问链表尾结点lu-tail；

步骤5：若lu-tail为空，表示初次收到该对象位置信息，则插入活动轨迹单元U(mv_n)到R-tree中，并记录它在R-tree中的位置uid，通过(U,uid)生成链表轨迹单元LU，并插入到链表；

步骤6：若lu-tail不为空，根据lu-tail中的uid信息找到并删除R-tree中的活动轨迹单元U(mv_n)，插入连续轨迹单元U(mv_s,mv_e)，插入活动轨迹单元U(mv_n)，最后在链表中删除tail；

步骤7：将连续轨迹单元和活动轨迹单元相继添加到链表尾部，建立索引，同时维护链表信息。

进一步，为实现对移动对象全时态的预测，本发明根据上述方法提出一种进行预测移动对象过去、现在和未来的全时态的方法，包括如下步骤：

一、移动对象过去位置索引

步骤1：确认查询时间t_q，判断t_q所属位置更新时间，如果t_q发生在位置更新的时刻执行步骤2；如果t_q发生在两个更新时刻之间执行步骤3，步骤4；

步骤2：t_q＝t_i，根据轨迹单元哈希表和双向链表确定t_q时刻移动对象的运行矢量mv(mid,t,v,rid,pos)，移动对象的位置则为mo_q＝(mid,rid,pos)；

步骤3：t_i<t_q<t_i+1，当前时刻移动对象未改变道路，移动对象在t_q-t_i间的平均速度为两次更新时刻移动对象速度的平均速度，则此时移动对象的位置mo_q，

步骤4：t_i<t_q<t_i+1，当前时刻移动对象已经改变道路，由此可知在时刻t_i+1的更新是因为道路的变化引起的，S_rid是道路的长度，

二、移动对象现在和未来位置的全时态位置预测

移动对象在一段时间内t_q-t_n的近似速度可以用移动对象过去时间内的平均速度与路段平均速度来综合表达，其中t_n为移动对象最后一次更新提交的时间，

步骤1：移动对象过去时间段内的平均速度

n为移动对象更新的次数，v_-i为移动对象距离当前第i次更新的速度，i∈[0,n]，θ为相关系数，θ∈[0,1]；

步骤2：道路中移动对象的平均速度

先求得路段A前路段的平均速度，在此之前必须确定要取路段前路段的数量w，处于同一条道路的移动对象之间的速度影响比较大，因此w值取道路Road中处于道路段A之后的路段数

w＝Road_all-Road_A 公式4；

其中Road_all为Road中路段总数，Road_A为路段A处于Road中的路段位置，对于其中任何一个道路段内的平均速度，为了降低时间开销提高预测准确度，采用t_n时刻的平均速度那么段内移动对象的平均速度可用如下公式得出：

其中ns为该段中移动对象的个数，v_i为移动对象距离t_n时刻最短的更新速度，那么对于个道路段来说，平均速度可以用公式表达为：

步骤3：移动对象的速度预测

综合公式3、4、5、6可得移动对象在未来一段时间内的预测速度，移动对象的预测速度可以由以下公式得到

移动对象的位置信息根据公式8计算得到

又进一步，在预测移动对象现在和未来位置的全时态位置的步骤1中，如果数据波动较大，则θ值应取大一些，可以增加近期数据对预测结果的影响；如果数据波动平稳，则θ值应取小一些。

进一步，本发明还根据上述的城市路网上基于R树的全时态移动对象索引方法提出一种进行索引一定时空范围内移动对象标识集合的方法，包括如下步骤：

步骤1：确认查询区域(△x,△y,△t)；

步骤2：根据查询窗口[(x1,y1),(x2,y2)]，搜索上层R*-tree，遍历到叶子结点，借助封闭路段信息e，求得查询窗口与封闭路段的交点；

步骤3：计算出交点在街道中的位置pos；

步骤4：査询区域与一个道路可能有多个交点，对应于一个偶对(rid,period)其中period包含一个或者多个pos信息，查询区域与多个道路相交，则对应于一组偶对，在rid对应的下层R-tree中查询与(period,△t)相交的轨迹单元，并输出轨迹单元的移动对象标识。

对比现有技术，本发明具有以下优点：

(1)为了提高索引移动对象轨迹信息的效率，引入了动态链表和移动对象哈希表，通过访问哈希表和动态链表可直接得出移动对象的历史轨迹信息，而无需搜索整个层移动对象树森林。

(2)为了应答在两次更新时间之间及最后一次更新之后的查询请求，采用了指数平滑法预测移动对象的速度，提高了移动对象在速度预测方面的精度。

(3)该索引结构支持移动对象的全时态位置查询，包括历史轨迹查询、实时位置查询、将来时刻位置预测和时空范围查询。

附图说明

图1是本发明的移动对象索引结构图。

图2是本发明的路网示意图。

图3是本发明的路网下移动对象三维时空轨迹及二维投影示意图。

图4是本发明的辅助索引结构。

具体实施方案

以下结合附图对本发明做进一步描述。

本发明提供的是支持移动对象全时态查询的索引结构，具体通过以下方案实现的：

一、体系结构

本索引分为上下双层索引结构，对静态路网信息采用R*-tree进行索引，并设计道路索引表来管理下层R-tree森林；对实时更新的移动对象运动轨迹采用结点更新代价较小的R-tree进行索引，设计轨迹单元索引表和双向链表辅助结构对索引协同管理。

下面我们给出几个部分的具体说明：

上层R*-tree索引：基于x-y平面的R*-tree索引，用于索引路网中的封闭路段，其叶子结点的数据结构为<M_xy,id,e>。M_xy表示包含该封闭路段的MBR，id为该叶子节点在索引中的标识，e为表示封闭路段。中间节点用<M_xy,id’,p>表示，M_xy为包含下层节点的MBR，id’为节点在索引中的标识，p为指向下层结点的指针。

下层R-tree索引：下层结构是一组与每条道路对应得R-tree索引，用于索引移动对象轨迹单元。其叶子结点的数据结构为<M_pt,id,mv_s,mv_e,mid>，M_pt为pos-t(位置和时间)平面中包含该轨迹的MBR，id为叶子节点的标识，mv_s和mv_e分别表示轨迹单元起始和结束运行矢量，如果为活动轨迹单元，则mv_e为空。mid表示移动对象的标识。中间结点数据结构为<M_pt,id’,p>，M_pt是包含下层结点的MBR，id’为中间节点标识，p为指向下层结点的指针。

道路索引表：其主要对下层基于道路的R-tree森林进行管理，其中key的值为道路标识rid，value为道路对应的R-tree索引。

轨迹单元索引表和双向链表：通过轨迹单元哈希表来管理移动对象的轨迹单元，key取值为移动对象标识mid，value为一个双向链表，用linkedlist表示，包含具体移动对象的所有轨迹单元LU，其按时间先后顺序加入到链表中，根据移动对象提供的位置信息不断进行更新。当收到一条位置更新信息时，首先形成轨迹单元U，将其插入到R-tree中，并返回在R-tree中的位置标识uid形成链表轨迹单元LU(U,uid)插入到linkedlist中，在维护索引的同时对linkedlist和轨迹单元哈希表进行维护。

本索引的建立和维护的主要步骤如下：

步骤1：根据移动对象位置更新信息m查找基于x-y的上层索引，找到移动对象所处的封闭路段e，并取出其相关信息。

步骤2：通过e信息以及x，y坐标信息求得移动对象在所属道路上的pos

位置，形成运行矢量mv＝(mid,t,v,rid,pos)。其中，pos∈[0，1]为对象所处道路的相对位置，t为采集该运行矢量的时间，v表示该移动对象在t时刻的速度，rid为道路标识。

步骤3：通过rid访问道路索引表，取得道路对应的R-tree；

步骤4：根据mid访问轨迹单元哈希表，获得轨迹单元链表，访问链表尾结点lu-tail；

步骤5：若lu-tail为空，表示初次收到该对象位置信息，则插入活动轨迹单元U(mv_n)到R-tree中，并记录它在R-tree中的位置uid，通过(U,uid)生成链表轨迹单元LU，并插入到链表；

步骤6：若lu-tail不为空，根据lu-tail中的uid信息找到并删除R-tree中的活动轨迹单元U(mv_n)。插入连续轨迹单元U(mv_s,mv_e)，插入活动轨迹单元U(mv_n)，最后在链表中删除lu-tail；

步骤7：将连续轨迹单元和活动轨迹单元相继添加到链表尾部，进而建立索引的同时并维护链表信息。

本发明基于上述的索引结构，采用能够预测移动对象过去，现在和未来的全时态的方法，具体情况如下：

1、移动对象过去位置索引

步骤1：确认查询时间t_q，判断t_q所属位置更新时间，如果t_q发生在位置更新的时刻执行步骤2；如果t_q发生在两个更新时刻之间执行步骤3，步骤4；

步骤2：t_q＝t_i，根据轨迹单元哈希表和双向链表确定t_q时刻移动对象的运行矢量mv(mid,t,v,rid,pos)，移动对象的位置则为mo_q＝(mid,rid,pos)。

步骤3：t_i<t_q<t_i+1，当前时刻移动对象未改变道路，移动对象在t_q-t_i间的平均速度为两次更新时刻移动对象速度的平均速度，则此时移动对象的位置mo_q，

步骤4：t_i<t_q<t_i+1，当前时刻移动对象已经改变道路，由此可知在时刻t_i+1的更新是因为道路的变化引起的，S_rid是道路的长度。

2、移动对象现在和未来位置的全时态位置预测

移动对象在一段时间内t_q-t_n，(其中t_n为移动对象最后一次更新提交的时间)的近似速度可以用移动对象过去时间内的平均速度与路段平均速度来综合表达。本专利采取指数平滑法，指数平滑法符合了移动对象速度预测的条件，既考虑移动对象过去的速度信息，又考虑处在移动对象前面的移动对象的速度信息，对于新数据赋予更大一些的权重而旧数据赋予较小的权重。

步骤1：移动对象过去时间段内的平均速度

n为移动对象更新的次数，v_-i为移动对象距离当前第i次更新的速度，i∈[0,n]。θ为相关系数，θ∈[0,1]，θ的取值容易受主观影响，如果数据波动较大，值应取大一些，可以增加近期数据对预测结果的影响。如果数据波动平稳，值应取小一些。

步骤2：道路中移动对象的平均速度

先求得路段A前路段的平均速度，但在此之前必须确定要取路段前路段的数量w。处于同一条道路的移动对象之间的速度影响比较大，因此w值取道路Road中处于道路段A之后的路段数

w＝Road_all-Road_A 公式4

其中Road_all为Road中路段总数，Road_A为路段A处于Road中的路段位置。对于其中任何一个道路段内的平均速度，为了降低时间开销提高预测准确度，采用t_n时刻的平均速度那么段内移动对象的平均速度可用如下公式得出

其中ns为该段中移动对象的个数，v_i为移动对象距离t_n时刻最短的更新速度，那么对于个道路段来说，平均速度可以用公式表达为

步骤3：移动对象的速度预测

综合公式3、4、5、6可得移动对象在未来一段时间内的预测速度，移动对象的预测速度可以由以下公式得到

移动对象的位置信息根据公式8计算得到

如图1所示，本索引方案分为上下两层，上层索引是基于路网的R*-tree索引结构，下层是基于pos-t平面的R-tree索引森林。下面分别对这两部分进行详细说明。

上层索引结构主要是对如图2所示的路网建立的R*-tree索引，以封闭路段为基本索引单位。图2为路网在上层索引中的对应关系，道路r1包含2个封闭路段(r1.a，r1.b)，于是形成2条记录存储于包含封闭路段MBR信息的叶子节点中。图2矩形框中示例给出了MBR3的范围。地图中的每个封闭路段被一个最小矩形框包含，所有最基本的最小边界矩形被融合为较大的MBR中，进一步迭代，重复上述过程计算更大的MBR，最终路网被一个最大的MBR包含。

如图1下半部分所示，下层索引结构对应R-tree森林，用于索引移动对象的轨迹单元。移动对象基于pos-t平面的轨迹片段被下层森林管理，如图3移动对象在道路上r运动，那么它的所有轨迹单元被道路r对应的R-tree管理，其中图3中的虚线部分表示活动轨迹单元，如果移动对象从一条道路运动到另一条道路，那么它的轨迹单元将被两个道路R-tree管理。

如图4所示的索引辅助结构提高了查询的效率。移动对象的运动轨迹跨越了多条道路，为了查询器历史运动轨迹，需要自上而下遍历整个R-tree森林，而采用哈希表和双向链表构成的辅助结构后，仅需将其加载到缓存进行维护，查找移动对象轨迹不再需要遍历R-tree森林，极大地提高了查询的效率。

本发明提供的是基于上述索引结构的移动对象位置的索引方法，该方法用于索引时空范围内移动对象标识集合，其包括如下步骤：

步骤1：确认查询区域(△x,△y,△t)；

步骤2：根据查询窗口[(x1,y1),(x2,y2)]，搜索上层R*-tree，遍历到叶子结点，借助封闭路段信息e，求得查询窗口与封闭路段的交点；

步骤3：计算出交点在街道中的位置pos；

步骤4：査询区域与一个道路可能有多个交点，对应于一个偶对(rid,period)其中period包含一个或者多个pos信息。查询区域与多个道路相交，则对应于一组偶对，在rid对应的下层R-tree中查询与(period,△t)相交的轨迹单元，并输出轨迹单元的移动对象标识。

以上所述仅为本发明的优先实施例，在本发明的精神和原则之内，所做的任何修改、等同替换、改进，均应在本发明的保护范围之内。

Claims (3)

1.城市路网上基于R树的全时态移动对象索引方法，其特征在于包括以下步骤：

步骤1：根据移动对象位置更新信息m查找基于x-y的上层索引，找到移动对象所处的封闭路段e，并取出其相关信息；

步骤2：通过e信息以及x，y坐标信息求得移动对象在所属道路上的pos位置，形成运行矢量mv＝(mid,t,v,rid,pos)，其中，pos∈[0，1]，为对象所处道路的相对位置，t为采集该运行矢量的时间，v表示该移动对象在t时刻的速度，rid为道路标识；

步骤3：通过rid访问道路索引表，取得道路对应的R-tree；

步骤4：根据mid访问轨迹单元哈希表，获得轨迹单元链表，访问链表尾结点lu-tail；

步骤5：若lu-tail为空，表示初次收到该对象位置信息，则插入活动轨迹单元U(mv_n)到R-tree中，并记录它在R-tree中的位置uid，通过(U,uid)生成链表轨迹单元LU，并插入到链表；

步骤6：若lu-tail不为空，根据lu-tail中的uid信息找到并删除R-tree中的活动轨迹单元U(mv_n)，插入连续轨迹单元U(mv_s,mv_e)，插入活动轨迹单元U(mv_n)，最后在链表中删除tail；

步骤7：将连续轨迹单元和活动轨迹单元相继添加到链表尾部，建立索引，同时维护链表信息。

2.一种根据权利要求1所述的城市路网上基于R树的全时态移动对象索引方法进行预测移动对象过去、现在和未来的全时态的方法，其特征在于包括如下步骤：

一、移动对象过去位置索引

步骤1：确认查询时间t_q，判断t_q所属位置更新时间，如果t_q发生在位置更新的时刻执行步骤2；如果t_q发生在两个更新时刻之间执行步骤3，步骤4；

步骤2：t_q＝t_i，根据轨迹单元哈希表和双向链表确定t_q时刻移动对象的运行矢量mv(mid,t,v,rid,pos)，移动对象的位置则为mo_q＝(mid,rid,pos)；

步骤3：t_i<t_q<t_i+1，当前时刻移动对象未改变道路，移动对象在t_q-t_i间的平均速度为两次更新时刻移动对象速度的平均速度，则此时移动对象的位置mo_q，

步骤4：t_i<t_q<t_i+1，当前时刻移动对象已经改变道路，由此可知在时刻t_i+1的更新是因为道路的变化引起的，S_rid是道路的长度，

二、移动对象现在和未来位置的全时态位置预测

移动对象在一段时间内t_q-t_n的近似速度用移动对象过去时间内的平均速度与路段平均速度来综合表达，其中t_n为移动对象最后一次更新提交的时间，

步骤1：移动对象过去时间段内的平均速度

n为移动对象更新的次数，v_-i为移动对象距离当前第i次更新的速度，i∈[0,n]，θ为相关系数，θ∈[0,1]；

步骤2：道路中移动对象的平均速度

先求得路段A前路段的平均速度，在此之前必须确定要取路段前路段的数量w，处于同一条道路的移动对象之间的速度影响比较大，因此w值取道路Road中处于道路段A之后的路段数

w＝Road_all-Road_A 公式4；

其中Road_all为Road中路段总数，Road_A为路段A处于Road中的路段位置，对于其中任何一个道路段内的平均速度，为了降低时间开销提高预测准确度，采用t_n时刻的平均速度那么段内移动对象的平均速度用如下公式得出：

其中ns为该路段中移动对象的个数，v_i为移动对象距离t_n时刻最短的更新速度，那么对于个道路段来说，平均速度用公式表达为：

步骤3：移动对象的速度预测

综合公式3、4、5、6得到移动对象在未来一段时间内的预测速度，移动对象的预测速度由以下公式得到

移动对象的位置信息根据公式8计算得到

3.一种根据权利要求1所述的城市路网上基于R树的全时态移动对象索引方法进行索引一定时空范围内移动对象标识集合的方法，其特征在于包括如下步骤：

步骤1：确认查询区域(△x,△y,△t)；

步骤2：根据查询窗口[(x1,y1),(x2,y2)]，搜索上层R*-tree，遍历到叶子结点，借助封闭路段信息e，求得查询窗口与封闭路段的交点；

步骤3：计算出交点在街道中的位置pos；

步骤4：査询区域与一个道路有多个交点，对应于一个偶对(rid,period)其中period包含一个或者多个pos信息，查询区域与多个道路相交，则对应于一组偶对，在rid对应的下层R-tree中查询与(period,△t)相交的轨迹单元，并输出轨迹单元的移动对象标识。