Orginal Article

Continuous Multi-scale Transformation of Urban Road Networks

  • XIE Kun ,
  • DENG Min , * ,
  • ZHANG Yunfei
Expand
  • School of Geosciences and Info-physics, Central South University, Changsha 410083, China
*Corresponding author: DENG Min, E-mail:

Received date: 2015-12-02

  Request revised date: 2016-03-16

  Online published: 2016-09-27

Copyright

《地球信息科学学报》编辑部 所有

Abstract

Presently, with the emergence of online mapping and Location Based Serveries (LBS), the continuous multi-scale representation of geo-spatial data has become a critical issue in the era of cartography and Geographic Information Science (GIS). This paper proposes a continuous transformation approach for the smooth representation of multi-scale road networks. Firstly, the topological relationships between the multi-scale road networks are reconstructed, and the feature correspondences are identified and classified into three categories (i.e., 1:1 matching features, 1:0 singleton features and M:N complex features) based on a probabilistic relaxation algorithm. For the 1:1 matching features, the corresponding points are identified based on the identical accumulated length ratio calculated from their starting vertex. The middle-scale representation of the 1:1 matching features is interpolated along the straight lines between the corresponding vertexes. For the 1:0 singleton features, considering that some road features are represented as a set of turns having approximately right angles, all the road vertexes containing an approximately right angle are detected and the edges connecting those vertexes are categorized into two characteristic sub-segments regarding as the left and right sides, which are then extended under the continuous direction. Then, the starting and ending nodes of the two characteristic sub-segment sets are selected to split the original road features into a set of divided line segments. Furthermore, the Douglas-Peucker algorithm is introduced to construct a hierarchical tree structure for each of the divided line segment, and all the divided line segments are hierarchically simplified by continuously moving the inner vertexes into the hierarchical tree structure. Moreover, the M:N complex features are visualized in gradually faded colors for a smooth representation. Finally, the continuous transformation results of different matching categories are integrated with respect to a given scale factor t. The experimental results have shown that the proposed approach can achieve a good performance for the continuous representation of the multi-scale urban road networks.

Cite this article

XIE Kun , DENG Min , ZHANG Yunfei . Continuous Multi-scale Transformation of Urban Road Networks[J]. Journal of Geo-information Science, 2016 , 18(9) : 1160 -1166 . DOI: 10.3724/SP.J.1047.2016.01160

1 引言

随着移动互联网技术与地理信息科学技术的交叉发展,导航与移动位置服务已逐渐渗入人们的日常生活中。在信息化大地图时代,用户对电子地图提出了连续多尺度可视化表达的新要求[1]。一方面,用户要求在进行放大、缩小等操作时保持自适应连续LoD(Level of Detail)表达;另一方面,针对不同应用需求,能为用户提供任意比例尺的电子地图。而传统的多版本地图数据库仅在预先定义的固定尺度间进行“跳跃”表达,无法实现不同比例尺数据间的连续过渡。因此,迫切需要研究多比例尺地图数据的平滑连续表达方法[2]
Morphing变换作为一种计算机视觉领域中的连续变形技术,最早用于不同图像间的连续渐变表达,近年来被引入地图制图领域,旨在建立2种比例尺矢量地图间的视觉连续变换模型[3-5]。彭东亮、邓敏等通过构建线状要素的层次弯曲森林,建立弯曲森林结点的对应关系,进行线性Morphing变换[6-7]。李精忠等同样提取线状要素的弯曲特征,利用模拟退火法匹配弯曲特征点并通过线性插值进行Morphing变换,这些方法主要用于河网等自然弯曲要素的连续综合[8]。谢天等建立了多比例尺建筑物多边形转向角函数的Morphing变换方法[9],保证建筑物的直角边界,Guibas等提出一种避免自相交的多边形Morphing变换方法[10]
现有方法主要研究多比例尺河流或建筑物中匹配要素Morphing变换,而对多比例尺道路网的连续变换研究相对较少。相比河网要素,道路要素一般不具有明显的弯曲特征,但存在复杂交叉口、多车道等复杂道路结构[11],由此导致多比例尺道路要素间更加复杂的对应关系。因此,针对多比例尺道路要素间复杂对应关系,本文基于Morphing变换技术提出了一种多比例尺道路网连续变换方法,实现2种不同比例尺道路网的中间尺度连续过渡变换。

2 多比例尺道路网的连续变换方法

基于Morphing的连续变换一般包含2个过程:建立对象间对应关系和确定移位路径以便形状插值。建立对象间对应关系包括建立要素与要素之间的对应关系,以及通过特征匹配建立线状要素中点与点之间的对应关系。本文考虑匹配要素的Morphing变换和未匹配要素的连续表达对多比例尺道路要素进行连续变换,具体流程如图1所示。首先,对2份不同比例尺的道路数据进行拓扑重建和要素匹配,将大小比例尺道路要素匹配关系分为1:1、1:0和M:NM>1或N>1);然后,对给定变换指标t,针对1:1对应线状要素、1:0单独线状要素和M:N复杂线状要素3种情况分别建立线性Morphing变换、连续综合变换和颜色渐变模型进行多尺度连续变换;最后,将3种类型要素的连续变换结果进行合并叠加,得到任一中间比例尺下的道路网可视化 表达。
Fig. 1 The framework of the proposed approach

图1 本文方法流程图

t是一个与比例尺有关的变换参数,不同t值下的变换结果代表了大小比例尺之间某一中间比例尺的地图表达,t愈大,愈接近小比例尺地图表达。本文预先设定统一变换指标t,t∈[0,1],其中t=0代表大比例尺道路网;t=1代表小比例尺道路网。

2.1 要素对应关系的建立

首先,对2份不同比例尺的道路路网数据进行拓扑打断和重新连接,建立结点-路段的图模型。然后,根据文献[12]提出的概率松弛匹配法识别2种比例尺道路要素间的对应关系。该方法首先对2份图结构建立初始匹配矩阵,再综合考虑邻近道路的兼容关系,启发式地更新初始匹配矩阵,使之达到全局最优,最后选取最终的匹配要素。
根据要素对应关系,本文将2份道路数据的线状要素分为3种情况:(1)1:1对应线状要素,即大小比例尺地图中首尾结点相互匹配的1:1匹配要素,主要为综合过程中被简化的线状要素;(2)1:0单独线状要素,即在小比例尺地图中找不到其在大比例尺地图中对应线状要素的独立要素,主要为综合过程中被删除的线状要素;(3)M:N复杂线状要素,即大比例尺地图中多个要素同时与小比例尺地图中一个或多个要素匹配的对应要素集(包括多对一和多对多),主要为综合过程被合并或分解的线状要素。对于上述3种情况,本文分别采用不同的连续变换策略实现不同比例尺道路数据的尺度连续 表达。

2.2 1:1对应要素的Morphing变换

基于复杂弯曲特征的特征匹配方法[6-8]适用于河网、海岸线、行政边界等具有明显弯曲特征的线状要素,而出于安全驾驶考虑,城市道路在设计时一般会避免较大的弯曲构造,因此本文采用简单线性Morphing变换方法[13-14]对1:1对应要素进行连续多尺度变换。该方法分为2步:特征点匹配和确定插值路径。本文根据到对应起点的累积距离与道路长度的比值进行1:1匹配要素的特征匹配。由于Morphing变换幅度较小,内插要素不易产生自相交,本文以对应顶点间的直线作为插值路径。
图2所示,对于大小比例尺地图上的1:1对应线状要素l1={v11,v12,…,v1m}和l2{v21,v22,…,v2n},l1为大比例尺要素,l2为小比例尺要素,v1j,v2j为对应要素的内部顶点。设r1(i)为l1沿起点v11到第i个顶点经过的线路长度占l1长度的比值,r2(j)同理。由于1:1对应要素的首尾顶点已匹配,本文根据到对应起点的累积长度比确定首尾顶点之外的其他匹配特征点。本文定义l1的顶点v1i(i>1)在l2中的对应顶点为l2中到起点v21的累积距离占l2长度的比值等于r1(i)的点,l2的顶点v2j( j>1)在l1上的对应顶点同理。将各要素的原有顶点与新增加的顶点按到起点的距离进行顺序排列,得到一一对应的匹配特征点序列l1={v′11, v′12, …, v′1p}和l2={v′21, v′22, …, v′2p},p=m+n-2。
Fig. 2 The result of linear Morphing method

图2 线性Morphing变换结果

lt={vt,i|i=1,…,p}为给定变换指标t下从要素l1l2的Morphing变换结果,则lt的顶点坐标可通过式(1)的线性插值计算得到。在具体实施过程中,线状要素随着t的增加逐渐从l1变换为l2
v t , i = ( 1 - t ) v 1 i + t v 2 i (1)

2.3 1:0单独线状要素的连续变换

对于1:0的单独线状要素,本文采用保持道路直角特征的Douglas-Peucker算法[15]进行连续变换。道路在进行设计时虽基本遵循连续平滑转向的原则,但在部分路段仍然包含了众多直角转弯(如十字交叉路口、丁字路口和一些延伸的巷子)。在多比例尺连续变换过程中需要尽可能保持这些转角特征。
首先,判断单独要素各顶点的夹角是否接近直角(即夹角是否介于90º-TA与90º+TA之间,TA为预设阈值),若是,则该顶点pi为直角特征顶点。
接着,顺序遍历顶点pi-1, pi-2,…, pj,…, p2,直到遍历到pj所连接的夹角小于180º-TA,则该顶点pj为左直角边的起点,同理可得右直角边终点。若左右直角边的长度都大于一定的阈值 T L ,则该道路转弯被视为直角特征,需在地图连续变换中保持。本文分别在直角顶点、左直角边起点和右直角边终点对该线状要素进行打断。
然后,利用文献[16]提出的弯曲深度二叉树表达方法对打断后道路要素进行层次分割,如图3所示。图3(b)中每个树结点与图3(a)中线状要素的各顶点对应,每个树结点同时存储各顶点到对应基线(图中虚线)的投影距离dpi与移位夹角θi。顶点p2到基线的投影点a需要满足式(2),则点p2的投影距离dp2等于p2与投影点a的距离,点p2的移位夹角θ2为α=∠ p2ap3
p 1 p 2 ¯ p 3 p 2 ¯ = p 1 a ¯ p 3 a ¯ (2)
进行连续简化的过程中,本方法要求各顶点的移位夹角保持不变。即按照层次从高到低的顺序对顶点进行移位。例如,对于第一层顶点p3,根据变换指标t确定该点移位后的投影距离dp′3=(1-tdp3,根据该投影距离在线段p3b之间寻找顶点将p3移动到p′3位置。而对于第二层顶点p2p4的移位,由于顶点p3的移位改变了第二层顶点的基线,因此需要重新计算p2p4投影点a′c′的位置,再执行前面的移位操作。进行下一层顶点的移位时皆需要重新计算该层顶点的基线,以保持该顶点与基线的夹角不变。这样线状要素随着t的增加而逐渐简化。同时,单独道路要素的渲染颜色随着t的增加而逐渐淡化,渲染颜色设定为(255×t,255×t,255×t)。
Fig. 3 The continuous generation of singleton linear features

图3 线状要素连续简化变形

2.4 M:N特殊要素的染色渐变

此外,道路网中还存在其他多对一、多对多的复杂对应关系需要单独处理。如图4(a)、(b)中的道路交叉口是一组相互匹配的线状要素,但无法进一步确定二者的顶点对应关系,从而使用线性Morphing变换。因此,本文综合考虑制图规则和图像Morphing渐变技术,选择染色淡出的方法进行处理。
Fig. 4 Special occasions

图4 特殊情况

具体地,对于变换指标t,大比例尺线状要素采用(255×t, 255×t, 255×t)的染色渐变方法,而小比例尺线状要素采用(255×(1-t), 255×(1-t), 255×(1-t))的染色渐变方式。可以看到,随着t的增加,大比例尺要素逐渐淡化,而小比例尺要素逐渐加深。该方法保留了来自2幅地图的空间信息,同时需要维护要素间的空间关系。

3 实验结果与分析

本文利用实际道路网数据验证方法的可行性和有效性。实验数据采用长沙市不同比例尺的道路网数据,数据来源于国家基础地理信息中心,图5(a)为1:1万比例尺道路数据(共1180条道路),图5(e)为1:5万比例尺道路数据(共696条道路)。通过Microsoft Visual Studio C# 2010与ArcGIS Engine 10.1编程软件实现本文提出的多尺度道路网Morphing变换方法,系统环境为Windows 7,CPU性能是Intel Core 5 Duo T6600 2.20 GHz,内存为4 GB(DDR3 1067 MHz)。
Fig. 5 The experimental results of the practical experimental data

图5 实际算例实验结果

本文将道路网的多尺度尺连续表达分为2部分执行:离线几何计算和在线可视化。离线几何计算包括要素对应关系的建立,1:1对应要素的特征点匹配,1:0单独线状要素的弯曲深度二叉树的建立等,与变换指标t无关的计算过程,可提前进行并建立缓存。在线可视化是指1:1对应要素的线性插值,1:0单独线状要素的连续变换和M:N的颜色渐变显示过程,本文采用了多线程算法对3种情况进行并行处理。据统计,实验数据在线可视化的单次平均运行时间约为120 ms。本文方法对小规模数据具有较高的运算效率,可满足实时可视化需求,但在处理大规模数据时还有一定差距,需要进行分区并行处理。
图5(b)、(c)、(d)分别显示了本文方法对不同变换程度(t=0.25,0.50,0.75)的连续表达结果。实验结果可发现,本文提出的多尺度连续表达方法能够得到良好的视觉连续表达效果,同时有效地保持了路网的整体模式和局部细节特征。
图6所示,本文利用顾及直角特征的Douglas-Peucker改进算法对单独线状要素进行连续变换。本文2.3节的参数TA取值为15°[17],该值能保证大多数情况下对道路的直角转角偏差的容忍,TL根据试验数据实际情况取值80 m。图6 (a)为1:1万比例尺道路数据中某条单独线状要素(点数为21个),图6(b)-(e)分别显示了变换指标t=0.25、0.50、0.75、1.0时的视觉连续变换结果。
Fig. 6 The generalization result of road features considering the right angle feature

图6 顾及直角特征的连续变换方法实验结果

将普通D-P算法与本文提出的顾及直角的连续变换算法进行了比较,结果如图7所示。在t=1时,2个方法保留的点数均为11个,图中红线为变换过程中保留下来的直角特征。分析实验结果可发现,在保持点数相同的情况下,随着变换程度加大,本文方法能够更有效地保持道路要素中的直角特征。
图8所示,对无法使用线性Morphing变换的M:N复杂线状要素,通过自适应改变要素渲染颜色实现了多尺度道路数据的视觉连续渐变效果。图8(a)为1:1万比例尺道路数据的环形交叉口,图8(e)为1:5万比例尺道路数据中十字交叉口,二者相互匹配,图8 (b)-(d)分别为变换程度t=0.25、0.5、0.75的要素渲染结果,视觉连续表达效果显著。
Fig. 7 The comparison between the commonapproach and the proposed approach

图7 普通线状要素连续变换方法与本文方法实验结果对比

Fig. 8 The result of the fade-out method

图8 染色淡出方法结果

4 结论与展望

本文提出了一种城市多比例尺道路网数据的连续表达方法,该方法首先识别不同比例尺道路数据的对应线状要素,然后分别对1:1对应线状要素、1:0单独线状要素和M:N复杂线状要素建立相应的orphing变换,连续综合变换和颜色渐变模型,实现不同比例尺道路网的连续可视化表达。通过实例验证,表明本文提出的道路网连续表达方法能够有效地融合不同比例尺地图的信息,在保持道路网整体模式和局部细节特征基础上,获得了较好的视觉连续表达效果。下一步研究工作将集中在挖掘道路的结构模式特征,进一步改进现有道路视觉连续表达效果。同时,研究大规模数据的分区并行计算方法,以保证运行效率和实时性需求。

The authors have declared that no competing interests exist.

[1]
王家耀. 地图制图学与地理信息工程学科发展趋势[J].测绘学报,2010,39(2):115-119.

[ Wang J Y.Development trends of cartography and geographic information engineering[J]. Acta Geodaetica et Cartographica Sinica, 2010,39(2):115-119. ]

[2]
Li Z, Zhou Q.Integration of linear and areal hierarchies for continuous multi-scale representation of road networks[J]. International Journal of Geographical Information Science, 2012,26(5):855-880.Spatial data can be represented at different scales, and this leads to the issue of multi-scale spatial representation. Multi-scale spatial representation has been widely applied to online mapping products (e.g., Google Maps and Yahoo Maps). However, in most current products, multi-scale representation can only be achieved through a series of maps at fixed scales, resulting in a discontinuity (i.e., with jumps) in the transformation between scales and a mismatch between the available scales and users' desired scales. Therefore, it is very desirable to achieve smoothly continuous multi-scale spatial representations. This article describes an integrated approach to build a hierarchical structure of a road network for continuous multi-scale representation purposes, especially continuous selective omission of roads in a network. In this hierarchical structure, the linear and areal hierarchies are constructed, respectively, using two existing approaches for the linear and areal patterns in a road network. Continuous multi-scale representation of a road network can be achieved by searching in these hierarchies. This approach is validated by applying it to two study areas, and the results are evaluated by both quantitative analysis with two measures (i.e., similarity and average connectivity) and visual inspection. Experimental results show that this integrated approach performs better than existing approaches, especially in terms of preservation of connectivity and patterns of a road network. With this approach, efficient and continuous multi-scale selective omission of road networks becomes feasible.

DOI

[3]
Nöllenburg M, Merrick D, Wolff A, et al.Morphing polylines: a step towards continuous generalization[J]. Computers, Environment and Urban Systems, 2008,32(4):248-260.We study the problem of morphing between two polylines that represent linear geographical features like roads or rivers generalized at two different scales. This problem occurs frequently during continuous zooming in interactive maps. Situations in which generalization operators like typification and simplification replace, for example, a series of consecutive bends by fewer bends are not always handled well by traditional morphing algorithms. We attempt to cope with such cases by modeling the problem as an optimal correspondence problem between characteristic parts of each polyline. A dynamic programming algorithm is presented that solves the matching problem in O(nm) time, where n and m are the respective numbers of characteristic parts of the two polylines. In a case study we demonstrate that the algorithm yields good results when being applied to data from mountain roads, a river and a region boundary at various scales.

DOI

[4]
Efrat A, Har-Peled S, Guibas L J, et al.Morphing between polylines[C]. The 20th Annual ACM-SIAM Symposium on Discrete Algorithms, 2001:680-689.

[5]
Sester M, Brenner C.Continuous generalization for fast and smooth visualization on small displays[C]. The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, 2004,34.

[6]
彭东亮,邓敏,刘慧敏.更充分利用独立弯曲结构的线状要素Morphing变换方法[J].测绘学报,2014,43(6):637-644.<p>本文提出了充分利用独立弯曲结构的线状要素Morphing变换方法。该方法首先对不同比例尺表达的对应线状要素分别构建约束Delaunay三角网并建立弯曲森林,然后进行弯曲匹配以获得对应弯曲。鉴于对应弯曲&ldquo;背面&rdquo;的独立弯曲结构隐藏于更高层次的大弯曲中,对对应弯曲重新构建约束Delaunay三角网进而建立其&ldquo;背面&rdquo;的弯曲森林并进行弯曲匹配得到新的对应弯曲,依此递归充分挖掘对应弯曲结构。在此基础上,将所有对应弯曲的对应始点和对应终点都作为断点切割原线状要素,获得对应线段。最后,采用线性插值算法建立各对应线段之间的对应点关系并以对应点间的直线作为移位路径进行Morphing变换。通过实例分析,验证了本文充分利用独立弯曲结构的方法能够提高对应弯曲特征点的识别能力,从而能够更好地保持弯曲特征点并改善Morphing变换效果。</p>

[ Peng D L, Deng M, Liu H M.Morphing transformationof linear featuresby using independent bend structuresmore sufficiently[J]. Acta Geodaetica et Cartographica Sinica, 2014,43(6):637-644. ]

[7]
彭东亮,邓敏,徐枫.顾及BLG树结构特征的线状要素Morphing变换方法[J].武汉大学学报,2012,37(9):1120-1125.对同一线状要素的不同比例尺表 达,借鉴Douglas-Peucker线状要素简化算法思想分别建立BLG树,通过对两BLG树从根结点到叶子结点进行层次匹配将两线状要素对应分割成 多对线段。在此基础上,借助线性插值算法进行Morphing变换。实验结果证明,此方法有效保持了原线状要素的结构特征,提高了Mor-phing变换 精度,改善了Morphing变换效果。

[ Peng D L, Deng M, Xu F.Morphing linear features considering their BLG-structures[J]. Geomatics and Information Science of Wuhan University, 2012,37(9):1120-1125. ]

[8]
李精忠,吴晨琛,杨泽龙,等.一种利用模拟退火思想的线状要素Morphing方法[J].武汉大学学报·信息科学版,2014,39(12):1446-1451.提出一种基于模拟退火思想的线状要素Morphing方法,针对同名线状要素在大小比例尺下的两种表达,首先利用约束Delaunay三角网提取小比例尺地图上线状要素的弯曲特征点,然后采用模拟退火技术在特征点与大比例尺线状数据顶点之间建立全局最优匹配,匹配结果将两线状要素分割成多对对应线段,最后针对每一对对应线段采用常规线性插值方法进行Morphing插值。模拟算例和实际数据实验证明,该方法较好地顾及了线状要素尺度变换过程中的弯曲化简、删除、夸大、典型化等综合操作,变换结果能有效地保持原线状要素的结构特征,提高了Morphing变换的精度。

DOI

[ Li J, Wu C M, Yang Z.A morphing method for linear features based on simulated annealing[J]. Geomatics and Information Science of Wuhan University, 2014,39(12):1446-1451. ]

[9]
谢天,李精忠.面状居民地Morphing变换的转向角函数法[J].测绘学报,2015,44(7):797-804.<p>提出了一种基于转向角函数的面状居民地morphing方法。针对同名居民地实体在两个不同比例尺下的表达,将传统矢量GIS的坐标串表达形式转为转向角函数表达形式,基于转向角函数分析多边形边特征在不同比例尺下的共性与差异,融合得到了任意中间比例尺下的多边形的转向角函数,最后将中间比例尺下的转向角函数还原为矢量坐标串形式,获得多边形在对应尺度的中间插值形状。试验证明,本文提出的基于转向角函数的面状居民地morphing方法对面状居民地多边形具有较好的适应性,能在保持直角化边界特征的前提下实现连续地图综合与多尺度表达。</p>

DOI

[ Xie T, Li J Z.Steering angle function algorithm of morphing of residential area[J]. Acta Geodaetica et Cartographica Sinica, 2015,44(7):797-804. ]

[10]
Guibas L, Hershberger J, Suri S.Morphing simple polygons[J]. Discrete and Computational Geometry, 2000,21(1):1-34.In this paper we investigate the problem of morphing (i.e., continuously deforming) one simple polygon into another. We assume that our two initial polygons have the same number of sides n , and that corresponding sides are parallel. We show that a morph is always possible through an interpolating simple polygon whose n sides vary but stay parallel to those of the two original ones. If we consider a uniform scaling or translation of part of the polygon as an atomic morphing step, then we show that O(n log n) such steps are sufficient for the morph. Furthermore, the sequence of steps can be computed in O(n log n) time.

DOI

[11]
William A M, Gordon A M.Automating the detection and simplification of junctions in road networks[J]. GeoInformatica, 1999,3(2):185-200.In this paper, the authors present a method for detecting and simplifying road junctions. The junctions are identified using combined spatial clustering and graph theory. The result is a generalization algorithm that can be used to derive generalized products from a single detailed database.

DOI

[12]
张云菲,杨必胜,栾学晨.利用概率松弛法的城市路网自动匹配[J].测绘学报,2012,41(6):933-939.lt;p>多源空间数据匹配是空间数据集成与互操作,变化检测与数据更新的重要前提。路网数据匹配在导航、智能交通和基于位置服务等领域具有重要的研究意义和实用价值。本文提出一种基于概率松弛方法的城市路网自动匹配方法,该方法首先通过路段间几何差异性估算候选路段的初始概率,然后根据邻接候选匹配路段的兼容性不断更新原概率矩阵直到收敛于某一极小值。最后基于收敛的概率矩阵计算各候选路段的结构相似性,并通过设定相应的规则选取和提炼1: 1, 1: M和M: N匹配对。实验选取中国武汉,瑞士苏黎世地区的OpenStreetMap数据与导航数据进行匹配算法的验证。结果表明:本文算法对非刚性偏差较大的路网数据能达到较高精度,不存在匹配方向性问题,且能够识别1: 0, 1: M和M: N匹配。</p>

DOI

[ Zhang Y F, Yang B S, Luan X C.Automated matching urban road networks using probabilistic relaxation[J]. Acta Geodaetica et Cartographica Sinica, 2012,41(6):933-939. ]

[13]
Haunert J H.Link based conflation of geographic datasets[C]. The 8th ICA Workshop on Generalization and Multiple Representation, 2005.

[14]
Oosterom P V, Bos J V D. An object-oriented approach to the design of geographic information systems[C]. The 1st Symposium on Design and Implementation of Large Spatial Databases, 1990:255-269.

[15]
Douglas D H, Peucker T K.Algorithms for the reduction of the number of points required to represent a digitized line or its caricature[J]. The Canadian Cartographer, 1973,10(2):112-122.All digitizing methods, as a general rule, record lines with far more data than is necessary for accurate graphic reproduction or for computer analysis. Two algorithms to reduce the number of points required to represent the line and, if desired, produce caricatures, are presented and compared with the most promising methods so far suggested. Line reduction will form a major part of automated generalization. Règle générale, les méthodes numériques enregistrent des lignes avec beaucoup plus de données qu'il n'est nécessaire à la reproduction graphique précise ou à la recherche par ordinateur. L'auteur présente deux algorithmes pour réduire le nombre de points nécessaires pour représenter la ligne et produire des caricatures si désiré, et les compare aux méthodes les plus prometteuses suggérées jusqu'ici. La réduction de la ligne constituera une partie importante de la généralisation automatique.

DOI

[16]
艾廷华,郭仁忠,刘耀林.曲线弯曲深度层次结构的二叉树表达[J].测绘学报,2001,30(4):343-348.地图综合要顾及目标的几何特征、语义特征和拓扑特征,其中地理意义是控制综合算子选择、参量调整的决定性因素.就线状要素而言,单从角度、距离、矢高等几何特征出发设计的曲线化简算法只能算作对曲线坐标串的几何压缩,不是真正意义上的地图综合.由于曲线的弯曲特征在表达线状地物地理特征上具有重要意义,对弯曲特征的识别、结构描述及操作分析成为目前线要素制图综合的研究热点.本文基于约束Delaunay三角网模型提出一种方法描述曲线弯曲特征在深度上的层次结构,对曲线上的矢量点构建三角网,在三角网覆盖区域里,由外向内进行三角形的“剥皮”操作,根据“剥皮”行进过程中遇到的特征三角形构建二叉树,实现大弯曲套小弯曲层次结构的表达.该方法基于Gestalt对称性、连续性原则,对二叉树结点进行考察,可提取认知意义上的真正弯曲.本文同时给出了弯曲特征二叉树在多边形(闭合曲线)综合化简中的算法设计及实验结果.

DOI

[ Ai T H, Guo R Z, Liu Y L.A binary tree representation of curve hierarchical structure in depth[J]. Acta Geodaetica et Cartographica Sinica, 2001,30(4):343-348. ]

[17]
Zhang M, Shi W, Meng L.A generic matching algorithm for line networks of different resolutions[C]. Proceedings of the 8th ICA Workshop on Generalization and Multiple Representation, 2005:7-8.

Outlines

/