矢量时空数据变化检测方法及其系统
技术领域
本发明涉及地理信息系统技术领域。具体地说,本发明涉及矢量时空数据变化检测技术。
背景技术
矢量时空数据变化检测技术用于确定实体亲缘关系、实体变化类型,提取时空数据,以实现历史数据和现势数据之间的一体化管理。目前,现有时空数据关联方法的研究还不够深入,还有待于进一步研究,时空变化类型需要进一步进行划分目前用于区分各种时空过程的时空变化类型尚没有形成共识,多种划分方法已被提出,如Christophe Claramunt给出了几种时空过程建模的设计类型,为新变化类型的定义提供可扩展的体系结构。K.Hornsby和M.Egenhofer深入研究了时空变化的语义,提出了9种时空变化类型。但这些变化类型对于时空过程尚未形成一个完整和系统地描述。
在矢量数据变化发现技术方面,目前,国内外相关研究不多,还未形成成熟的矢量数据变化发现技术,现有的矢量数据变化发现技术的不足主要体现在:现有的技术方法主要以矢量要素的对象ID或空间位置为变化发现的匹配基础,不具有通用性,仅适用于某一应用领域;缺乏多时态矢量数据变化的判别方法。
发明内容
因此,根据本发明的优选实施方式,本发明解决了现有技术中在地理信息数据时空变化检测时不能有效地对多时态矢量数据进行变化发现的技术问题。
本发明通过制定时空数据匹配规则的方法,设定好相应的时空数据匹配规则,制定时空数据匹配规则库,以便于用户可选择与应用相适应的时空数据匹配规则,增强了矢量数据变化发现技术的通用性;通过设计时空数据变化规则模版,以便于用户根据应用需求,来确定多时态矢量数据变化判别标准,实现多时态矢量数据变化的判别方法。
根据本发明的一个方面,本发明提供了一种矢量时空数据变化检测方法,所述方法包括以下步骤:
多时态变化检测步骤:针对多时态矢量数据来源不同,可采用多版本数据变化发现方式、增量数据变化发现方式或者人工实时更新变化发现方式来进行所述时空数据的变化发现;
时空数据融合步骤:确定时空数据变化检测结果的正确性,删除没有发生变化的历史数据,将发生了变化的历史数据保留到历史数据中,将时空数据间的变化关系保存在过程数据中,将现势数据保存到现势数据中,对多时态时空数据进行融合,形成具有历史亲缘信息的时空数据。
其中,当采用所述多版本数据变化发现方式或者增量数据变化发现方式时,所述多时态变化检测步骤进一步包括根据用户选择的时空数据匹配规则和时空数据变化规则进行自动变化检测的自动变化检测步骤;
所述自动变化检测步骤进一步包括:
地理要素匹配步骤:基于用户所选的时空数据匹配规则,选取不同时态地理要素作为相关联的粗侯选集,对地理要素关联的粗侯选集中的时态信息进行冲突检查,如果地理要素间的时态信息不发生冲突,将粗侯选集中符合条件的地理要素,进一步选入精细侯选集;
几何变化检测步骤:基于用户选择的时空数据几何变化规则,以及基于在“地理要素匹配步骤”选取的所述精细侯选集中相关联的地理要素,通过自动变化发现算法,检测时空数据变化,对时空数据几何变化进行定性表达,确定时空地理要素的变化关系,及其历史亲缘关系,形成具有历史亲缘信息的时空数据;其中,所述几何变化规则包括点状地理要素变化规则、线状地理要素变化规则、面状地理要素变化规则;
属性变化检测步骤:根据预先设定的时空数据语义变化规则,通过属性值的比较分析,检测时空数据的语义变化,记录属性变化的结果;
其中,所述“几何变化检测步骤”具体包括:
当用户选择点状地理要素变化规则时,执行点状地理要素的变化检测步骤,所述点状地理要素的变化检测步骤包括:
步骤1.1.首先对旧时态点状地理要素和新时态点状地理要素的空间参考和尺度信息进行对比,如果所述空间参考和尺度信息不一致,则退出;如果所述空间参考和尺度信息一致,则进行步骤1.2;
步骤1.2.以旧时态点状地理要素集里面的地理对象为圆心,以阈值d为半径做缓冲区叠加分析;
步骤1.3.判断新时态点状地理要素集里面的点是否落在缓冲区范围内,如果是,则判断为旧时态地理要素和新时态地理要素间没有发生变化,否则,则判断为新时态地理对象的变化类型是新增的,旧时态地理对象的变化类型为消失;
当用户选择线状地理要素变化规则时,执行线状地理要素的变化检测步骤,所述线状地理要素的变化检测步骤包括:
步骤2.1.首先对旧时态线状地理要素和新时态线状地理要素的空间参考和尺度信息进行对比,如果所述空间参考和尺度信息不一致,则退出;如果所述空间参考和尺度信息一致,则进行步骤2.2;
步骤2.2.分别计算旧、新时态线状地理要素的外接矩形;
步骤2.3.判断两个外接矩形的关系,具体包括:
如果旧、新时态线状地理要素的外接矩形的空间关系为相离,则判断为旧线状地理要素的变化类型为消失,新时态线状地理要素的变化类型为新增;
如果旧、新时态线状地理要素的外接矩形的空间关系为相交,则对新、旧线状地理要素进行化简,计算新、旧线状地理要素的特征点,然后对新、旧线状地理要素的特征点之间的距离进行计算;如果新、旧线状地理要素的全部特征点之间的距离在一定的容差范围内,则判断为新、旧线状地理要素没有发生变化;如果新、旧线状地理要素的所有特征点之间的距离超过容差,则判断为旧线状地理要素的变化类型为消失,新时态线状地理要素的变化类型为新增。如果新、旧线状地理要素的部分特征点之间的距离超过容差,则认为旧时态线状地理要素变化为新时态线状地理要素;
当用户选择面状地理要素变化规则时,执行面状地理要素的变化检测步骤,所述面状地理要素的变化检测步骤包括:
步骤3.1.首先对旧时态面状地理要素和新时态面状地理要素的空间参考和尺度信息进行对比,如果所述空间参考和尺度等信息不一致,则退出;如果所述空间参考和尺度信息一致,则进行步骤3.2;
步骤3.2.分别计算旧、新时态面状地理要素的外接矩形;
步骤3.3.判断两个外接矩形之间的关系,具体包括:
如果旧、新时态面状地理要素的外接矩形的空间关系为相离,则判断为旧面状地理要素的变化类型为消失,新时态面状地理要素的变化类型为新增;
如果旧、新时态面状地理要素的外接矩形的空间关系为相交,则计算新、旧面状地理要素的重叠度,如果新、旧面状地理要素的重叠度大于98%,则判断为新、旧面状地理实体没有发生变化;如果新、旧面状地理要素的重叠度小于2%,则判断为旧面状地理要素的变化类型为消失,新时态面状地理要素的变化类型为新增;如果新、旧面状地理要素的重叠度大于2%且小于98%,则判断为旧面状地理要素变为新时态面状实体。
优选的,所述自动变化检测步骤还包括:
人工修正变化步骤:如果地理实体的变化类型不符合实际情况,则根据实际变化情况修改实体的变化类型。
优选的,所述多时态矢量数据变化发现方式为:通过两个不同时态版本数据,进行变化发现;
所述增量数据变化发现方式为:通过导入发生变化的增量数据与旧时态版本数据进行变化发现,确定多时态数据的变化类型;
所述人工实时更新变化发现方式为:通过新时态影像与旧时态矢量时空数据进行叠加,勾绘发生变化的时空地理实体,形成新时态的地理要素,与旧时态版本矢量数据进行变化发现。
本发明提供了矢量时空数据变化检测方法与系统,该发明技术能对多时态矢量数据进行变化发现,并融合形成具有历史信息的时空数据,能更好地实现对时空数据的时空查询、时空分析,能够满足历史数据的保存和恢复、变化的跟踪和预测等要求。可广泛应用于土地利用变化、城市演化等涉及时空变化的管理和应用,直观生动地表达地理现象的变化过程及规律,从而达到对地理现象进行过程推演、过程再现、实时跟踪以及运动模拟,表现地理现象的内在本质和发生规律的目的。
综上所述,本发明的信矢量时空数据变化检测方法与系统能够将多时态数据有效融合起来,形成具有历史信息的时空数据,实现时空数据的一体化管理,能够很好地表达地理现象的变化,生动直观地表现地理现象的变化过程。
附图说明
下面结合附图详细说明本发明的特征和优点,在附图中:
图1示出了根据本发明的优选实施方式的矢量时空数据变化检测流程图。
图2示出了根据本发明的优选实施方式的点的变化检测算法的流程图。
图3示出了根据本发明的优选实施方式的线的变化检测算法的流程图。
图4示出了根据本发明的优选实施方式的面的变化检测算法的流程图。
具体实施方式
下面,参照附图更全面地描述本发明、示出本发明的各方面的优选实施方式。为了更好地突出本发明的实质性特点,在以下的描述中将省略对公知常识的描述。
图1图示了根据本发明的优选实施方式的矢量时空数据变化检测流程图,下面分别按照步骤详细说明:
本发明的矢量时空数据变化检测方法,分为以下几个步骤:
1.制定时空数据匹配规则步骤,本发明的矢量时空数据变化检测方法包含三种匹配规则:a、基于几何的匹配规则,即,通过地理实体的几何位置的相互关系来确定地理实体间的历史亲缘对象;b、基于语义的匹配规则,即,通过地理实体的语义来确定地理实体间的历史亲缘对象;c、几何与属性相结合的匹配规则,即,通过地理实体的几何与语义相结合来确定地理实体间的历史亲缘对象。
2.制定时空数据变化规则步骤,包括:
a、制定时空数据几何变化规则:
根据不同的地理要素类型(点、线、面),设定规则,相同的地理实体在不同时态,由于数字化及变化等原因,其坐标点的个数和状态都可能不一致,因此首先设置判断变化的阈值,针对点状地理实体,设定两个点之间的距离是否超过阈值来判断点发生何种变化;针对线状地理实体,设计基于特征点的方法,采用数据压缩算法提取不同时态中的线状地理实体特征点,设定特征点的阈值,根据不同时态线状地理实体的特征点位置变化情况,来判断线状地理实体发生何种变化;针对面状地理实体,将面进行叠置分析,根据不同时态面状地理实体的面积、周长变化情况,以及不同时态面状地理实体位置等几何位置变化情况,判断发生何种类型的变化。
b、制定时空数据语义变化规则:
制定时空数据语义变化规则,判断两个地理实体间的属性是否发生变化。首先,制定哪些地理实体的属性变化需要监测,设定地理实体的属性字段间的对应关系;再制定属性值间的变化规则,包括数值型属性的变化规则,如数值型属性在一定的误差范围内,则认为没有变化;字符型属性变化规则,如字符型属性,可以采用正则表达式的方式,来判断字型属性是否发生变化。
3.多时态矢量数据变化发现步骤
本发明的多时态矢量数据变化发现方法,针对不同数据来源,可以采用三种方式进行变化发现:
(1)多版本数据变化发现,即,通过两个不同时态版本数据,进行变化发现;
(2)增量数据变化发现,即,通过导入发生变化的增量数据与旧时态版本数据进行变化发现;
(3)人工实时更新变化发现,即,通过新时态影像与旧时态矢量时空数据进行叠加,勾绘发生变化的时空地理实体,形成新时态的地理要素,与旧时态版本矢量数据进行变化发现。
前两种种变化发现方式,可根据用户需要选择采用自动变化检测步骤,或者采用手工变化发现步骤;
1)自动变化检测步骤:
自动变化检测步骤主要包括地理要素匹配步骤、几何变化检测步骤、属性变化检测步骤、人工修正变化关系步骤。自动变化检测,首先进行地理要素匹配步骤、再进行几何变化检测步骤、后进行属性变化检测步骤、最后进行人工修正变化关系步骤。详细的各步骤说明如下:
a.地理要素匹配步骤:
首先用户根据时空数据匹配规则步骤(步骤1)所形成的时空数据匹配规则中选取一种时空数据匹配规则,通过用户所选的时空数据匹配规则,选取不同时态地理要素作为相关联的粗侯选集,对地理要素关联的粗侯选集中的时态信息进行冲突检查,如果地理要素间的时态信息不发生冲突,将粗侯选集中符合条件的地理要素,进一步选入精细侯选集。
b.几何变化检测步骤:
首先用户根据矢量地理要素的类型(点、线、面),选择步骤2制定的与地理要素类型相一致的时空数据几何变化规则,,基于步骤a选取的精细侯选集中相关联的地理要素,通过自动变化发现算法,检测时空数据变化,对时空数据几何变化进行定性表达,确定时空地理要素的变化关系,及其历史亲缘关系,形成具有历史亲缘信息的时空数据。
其中,包括维护时态的一致性步骤:如不同时态的地理要素具有历史亲缘关系,对该时态的地理要素间的时态信息进行自动维护,更改旧时态地理要素的结束时间,将旧时态地理要素的结束时间设置为新时态地理要素的起始时间。
c.属性变化检测步骤:
完成地理要素匹配、几何变化检测步骤后,根据制定时空数据变化规则步骤(步骤2)中所预先设定的时空数据语义变化规则,通过属性值的比较分析,检测时空数据的语义变化,记录属性变化的结果。
d.人工修正变化关系步骤:
基于自动变化检测步骤形成的时空数据,确定实体或实体的变化类型是否符合实际情况,例如,旧时态点状实体、旧时态线状实体、旧时态面状实体、新时态点状实体、新时态线状实体或新时态面状实体的变化类型是否符合实际情况。如果确定某个实体或某些实体的变化类型不符合实际情况,根据实际变化情况修改的实体的变化类型,检查实体变化类型的一致性(实体的变化类型是不能相互矛盾的,用户更改实体变化类型的过程,可能使得实体变化类型发生矛盾,根据关联规则来消除这些矛盾,也称为实体变化类型的一致性检查),修正这些实体的变化类型。
其中几何变化检测方法是按照点状地理要素、线状地理要素、面状地理要素分别进行,下面就点、线、面变化检测分别进行详细说明:
(a)点状地理要素变化检测:
如图2所示,当用户选择点状地理要素变化规则时,执行点状地理要素的变化检测步骤,所述点状地理要素的变化检测步骤包括:
步骤1.1.首先对旧时态点状地理要素和新时态点状地理要素的空间参考和尺度信息进行对比,如果所述空间参考和尺度信息不一致,则退出;如果所述空间参考和尺度信息一致,则进行步骤1.2;
步骤1.2.以旧时态点状地理要素集里面的地理对象为圆心,以阈值d为半径做缓冲区叠加分析;
步骤1.3.判断新时态点状地理要素集里面的点是否落在缓冲区范围内,如果是,则判断为旧时态地理要素和新时态地理要素间没有发生变化,否则,则判断为新时态地理对象是新增的,旧时态地理对象的变化类型为消失。
(b)线状地理要素变化检测:
如图3所示,当用户选择线状地理要素变化规则时,执行线状地理要素的变化检测步骤,所述线状地理要素的变化检测步骤包括:
步骤2.1.首先对旧时态线状地理要素和新时态线状地理要素的空间参考和尺度信息进行对比,如果所述空间参考和尺度信息不一致,则退出;如果所述空间参考和尺度信息一致,则进行步骤2.2;
步骤2.2.分别计算旧、新时态线状地理要素的外接矩形;
步骤2.3.判断两个外接矩形的关系,具体包括:
如果旧、新时态线状地理要素的外接矩形的空间关系为相离,则判断为旧线状地理要素的变化类型为消失,新时态线状地理要素的变化类型为新增;
如果旧、新时态线状地理要素的外接矩形的空间关系为相交,则对新、旧线状地理要素进行化简,计算新、旧线状地理要素的特征点,然后对新、旧线状地理要素的特征点之间的距离进行计算;如果新、旧线状地理要素的全部特征点之间的距离在一定的容差范围内,则判断为新、旧线状地理要素没有发生变化;如果新、旧线状地理要素的所有特征点之间的距离超过容差,则判断为旧线状地理要素的变化类型为消失,新时态线状地理要素的变化类型为新增。如果新、旧线状地理要素的部分特征点之间的距离超过容差,则认为旧时态线状地理要素变化为新时态线状地理要素;
(c)面状地理要素变化检测:
如图4所示当用户选择面状地理要素变化规则时,执行面状地理要素的变化检测步骤,所述面状地理要素的变化检测步骤包括:
步骤3.1.首先对旧时态面状地理要素和新时态面状地理要素的空间参考和尺度信息进行对比,如果所述空间参考和尺度等信息不一致,则退出;如果所述空间参考和尺度信息一致,则进行步骤3.2;
步骤3.2.分别计算旧、新时态面状地理要素的外接矩形;
步骤3.3.判断两个外接矩形之间的关系,具体包括:
如果旧、新时态面状地理要素的外接矩形的空间关系为相离,则判断为旧面状地理要素的变化类型为消失,新时态面状地理要素的变化类型为新增;
如果旧、新时态面状地理要素的外接矩形的空间关系为相交,则计算新、旧面状地理要素的重叠度,如果新、旧面状地理要素的重叠度大于98%,则判断为新、旧面状地理实体没有发生变化;如果新、旧面状地理要素的重叠度小于2%,则判断为旧面状地理要素的变化类型为消失,新时态面状地理要素的变化类型为新增;如果新、旧面状地理要素的重叠度大于2%且小于98%,则判断为旧面状地理要素变为新时态面状实体。
2)手工变化检测步骤:
分别选择不同时态的地理要素,通过软件工具,为选择的不同时态的地理要素,设置不同时态地理要素的变化关系。
3)编辑更新
主要通过对旧时态的矢量数据和新时态影像数据,进行叠加分析,查找旧时态地理要素的变化情况,如果发生变化,则通过编辑工具对新时态地理要素进行勾绘,形成新时态的矢量地理要素,再调用自动变化检测步骤,形成新旧时态地理要素的变化关系。
4)时空数据融合步骤:
确定时空数据变化检测结果的正确性,删除没有发生变化的旧时态数据,将发生了变化的旧时态数据保留到历史数据中,将时空数据间的变化关系保存在过程数据文中,将新时态数据保存到现势数据中,对多时态时空数据融合,形成具有历史亲缘信息的时空数据。
以上内容仅为本发明的较佳实施例,对于本领域的普通技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,本说明书内容不应理解为对本发明的限制。