基于GIS的任意发生元Voronoi图逼近方法

许多地理问题的空间分析中需要采用Voronoi图,但是目前我们尚缺乏一些简单的易于实现的构建任意发生元Voronoi图的方法,也缺乏一个能直接生成任意发生元Voronoi图的软件,为此我们提出了一种基于GIS的构建任意发生元的未加权Voronoi图的逼近方法。首先用有限点来逼近原始发生元,然后构建这些点发生元Voronoi图,最后消除那些属于同一发生元的顶点和边,即得到原始发生元的逼近的Voronoi图。在该算法的具体实现过程中,充分利用了现有GIS软件可以生成点发生元Voronoi图的特性和处理空间数据的能力。试验结果表明,这种方法可以生成未加权的任意形状发生元的逼近Voronoi图,能满足地理问题空间分析的需要。如地理客体可以是点状地理客体(城市、县城、交通枢纽、商业中心和金融中心等)、线状地理客体(交通运输线、经济地带和河系等)、面状地理客体(经济区、公园和绿地等)或者它们的组合,它们的空间影响范围或空间服务范围都可以采用Voronoi图来界定。     

Voronoi tessellation on the ellipsoidal earth for vector data

Voronoi tessellation, and its dual the Delaunay triangulation, provide a cohesive framework for the study and interpretation of phenomena of geographical space in two and three dimensions. The planar and spherical solutions introduce errors in the positional accuracy of both Voronoi vertices and Voronoi edges due to errors in distance computations and the path connecting two locations with planar lines or great circle arcs instead of geodesics. For most geospatial applications the introduction of the above errors is insignificant or tolerable. However, for applications where the accuracy is of utmost importance, the ellipsoidal model of the Earth must be used. Characteristically, the introduction of any positional error in the delimitation of maritime zones and boundaries results in increased maritime space for one state at the expense of another. This is a situation that may, among others, have a serious impact on the financial activities and the relations of the states concerned. In the context of previous work on maritime delimitation we show that the Voronoi diagram constitutes the ideal solution for the development of an automated methodology addressing the problem in its entirety. Due to lack of a vector methodology for the generation of Voronoi diagram on the ellipsoid, the aforementioned solution was constrained by the accuracy of existing approaches. In order to fill this gap, in this paper we deal with the inherent attributes of the ellipsoidal model of the Earth, e.g. the fact that geodesics are open lines, and we elaborate on a methodology for the generation of the Voronoi diagram on the ellipsoid for a set of points in vector format. The resulting Voronoi diagram consists of vertices with positional accuracy that is only bounded by the user needs and edges that are comprised of geodesics densified with vertices equidistant to their generators. Finally, we present the implementation of the proposed algorithm in the Python programming language and the results of two case studies, one on the formation of closest service areas and one on maritime boundaries delimitation, with the positional accuracy set to 1 cm.     

基于网络加权Voronoi图分析的南京市商业中心辐射域研究

道路网络在城市社会和经济活动中发挥重要作用.它是联系商业中心与需求间的桥梁.道路系统的更新与变化必然会对它们的关联产生影响Voronoi图足划分势力空间的有效工具,基于道路网络分析构建的加权Voronoi图可以模拟城市商业中心功能辐射的实际情形.本文给出了基于道路网络分析的Voronoi图和加权Voronoi图构建方法.通过对各种Voronoi图的比较分析,认为基于网络路径时间分析构建的加权Voronoi面域图在商业中心辐射域模拟中最为有效.将空间相互作用和断裂点理论与网络加权Voronoi面域图相结合.模拟不同交通条件下南京市商业中心辐射域空间格局及其需求覆盖变化,并提取了辐射域的覆盖面积、覆盖人口、最远辐射路径距离和平均辐射路径距离等信息.对商业中心的地铁邻近度与其辐射域的扩张面积进行相关分析.采用域邻接矩阵描述空间划分平衡与不平衡之间的相互转变过程,可分析道路系统变更后商业中心之间空间竞争的变化趋势.分析结果表明商业中心的地铁邻近度与辐射域的扩张具有较大的相关性,快速轨道交通在商业中心空间竞争中发挥关键作用.     

基于道路网络的商业网点市场域分析

社会经济活动中人们的空间行为往往是基于道路网络来实现的。但是,目前多数研究都假设地理空间是一个均质空间,采用基于平面欧氏距离的空间分析方法,这是有局限的。本文阐述了基于网络距离的网络Voronoi图基本概念和实现方法,以武汉市商业零售连锁企业为例,分别采用基于欧氏距离的普通Voronoi图方法和基于网络距离的网络Voronoi图方法来确定商业零售网点的市场域,结果表明武汉市主城区商业网点市场域较小、主城区外围市场域较大。市场域大小与城市路网密度呈现一定相关关系,路网密度高、市场域小,路网密度低、市场域大。两种方法的计算结果存在一定差异,差异大小与路网密度有关,路网密度大、差异小,路网密度小、差异大。这表明在路网密度大的情况下,可以采用普通Voronoi图粗略地模拟商业网点的市场域。     

基于局部聚类的网络Voronoi图生成方法研究

提出一种将网络约束下的Voronoi和空间聚类相结合的方法,通过构造局部的聚类分析方法对网络边进行加权,根据实际的点过程性质可以把权重定义为加权或者乘权,进行标准化后与道路段本身长度融合进行计算,依此生成网络Voronoi图,以期理解城市街道的空间特性。以武汉市江汉区为例,对城市网格管理系统产生的城市事件进行算法验证,结果表明,该方法提供了一种灵活的网络约束下的服务区域划分工具,可用于基于网络空间点过程影响下的服务区划分,也可用于系统性地定量刻画城市管理的动态特性。     

黄土丘陵沟壑区农村居民点空间重构——以榆中县为例

农村居民点空间优化研究是进行农村居民点用地空间重构的依据。以兰州市榆中县为例,运用GIS 空间分析能力,计算农村居民点的影响力来划分居民点类型,应用加权Voronoi图的空间分割功能实现宏观上农村居民点的空间重构,确定各搬迁居民点的安置去向,同时计算各个加权Voronoi 图之中的居民点的相对“居住场势”的大小,实现微观上搬迁农户的空间效用增加。结果表明：榆中县共有中心居民点面积6082.60 hm² （1107 个居民点） 占61.25%;Ⅰ类保留型居民点面积2707.55 hm² （4083 个居民点） 占27.26%;Ⅱ类保留型居民点面积235.68 hm² （679 个居民点） 占2.37%;零星居民点面积905.26 hm² （2854 个居民点） 占9.12%。通过加权Voronoi 图的宏观空间分割结合微观的“居住场势”测算,使得居民点空间重构的方案更加合理。     

基于互联网信息的城市商业游憩网点空间竞争域划分研究

提出了一种基于互联网信息的空间竞争域划分方法。在Voronoi图的模型中尝试以网络平台数据作加权因子,使空间竞争域体现互联网中的消费者偏好差异,通过对加权Voronoi图的类型划分和地理学解释,描述竞争环境的空间特征。以北京市健身俱乐部为例,选择大众点评网、爱帮等网络平台信息作为权重因子,构建样本网点的传统Voronoi、互联网信息加权Voronoi计算模型并对比边界线、面积等空间特征分析的结果表明:互联网信息对城市商业游憩网点的空间竞争力产生着不同的影响强度,对Voronoi图空间特征的分类和地理解释可以指示商业网点跨越距离阻碍的消费吸引力及竞争差异,弥补传统加权方法未考虑消费者偏好和购物习惯而对竞争域划分表达不足的缺点。     

Population landscape: a geometric approach to studying spatial patterns of the US urban hierarchy

We present a geometric and graphic approach to studying spatial patterns of urban hierarchy in the US. The multiplicatively weighted Voronoi diagram is found to be effective for visualizing theoretical regions delineated by socio‐economic variables. The population landscape of the continental US demonstrates overall and stepwise patterns reflecting population, neighborhood and distance, with overwhelming influence from huge metropolitan areas. Stepwise exploration and cluster analysis of the spatial pattern reveal an urban hierarchy. Attributes and arrangement are the two important factors of urban hierarchy, with attribute having a stronger local influence and arrangement having a stronger global influence. The study also presents a variation of Zipf's law to visualize the rank‐size distribution from tabular and statistical space to map space.     

Modelling three‐dimensional geoscientific fields with the Voronoi diagram and its dual

Fields as found in the geosciences have properties that are not usually found in other disciplines: the phenomena studied are often three‐dimensional (3D), they tend to change continuously over time, and the collection of samples to study the phenomena is problematic, which often results in highly anisotropic distributions of samples. In the geographical information system (GIS) community, raster structures (voxels or octrees) are the most popular solutions, but, as we show in this paper, they have shortcomings for modelling and analysing 3D geoscientific fields. As an alternative to using rasters, we propose a new spatial model based on the Voronoi diagram (VD) and its dual the Delaunay tetrahedralisation (DT), and argue that they have many advantages over other tessellations. We discuss the main properties of the 3D VD/DT, present some GIS operations that are greatly simplified when the VD/DT is used, and, to analyse two or more fields, we also present a variant of the map algebra framework where all the operations are performed directly on VDs. The usefulness of this Voronoi‐based spatial model is demonstrated with a series of potential applications.     

基于Voronoi图剖分的区域旅游中心城市体系——以苏浙沪地区为例

在区域范围的城市旅游体系中,引入"中心地的中心性"概念,可以在这种旅游性的空间网络中考察城市地位的相对重要性,划分出城市的等级层次,判断出不同级别的中心城市,更可以进一步得到区域空间的城市旅游构成关系.Voronoi图在地理学中是一种利用中心地的中心性广泛用于空间分割、空间邻域查找的空间剖分方法.引入这种方法,以苏浙沪地区为例,在25个地级以上城市中,利用Voronoi图的空间分割原理,找出不仅具有规模性,而且具有空间组织功能的旅游中心城市.     

非均质栅格Voronoi图生成方法研究

提出一个非均质栅格Voronoi图的生成方法.假定在Voronoi领域扩张中,经过每个栅格的距离权值不相等.由于非均质栅格Voronoi图生成复杂,提出了一种基于A*启发式空间搜索确定每个栅格的归属来产生Voronoi区域的方法.讨论了非均质栅格Voronoi图的基本定义、启发式函数确定、具体生成方法等.利用C#语言实现了所提出的方法,并分析了该方法的时间复杂度.但该方法仅适合于简单点发生元情况,在面向复杂发生元的情况下,需进一步研究更有效的方法.     

一种用于界定经济客体空间影响范围的方法——Voronoi图

经济客体的空间影响范围界定十分复杂,但在区域规划和城市规划中有着重要的理论和实际意义。该文提出可采用Voronoi图方法用于经济客体的空间影响范围界定,介绍了Voronoi图的基本原理和其若干扩展,编写了生成Voronoi图的程序,最后以城市为例探讨了Voronoi图在经济客体空间影响范围界定中的应用。     

Development of Voronoi-based cellular automata -an integrated dynamic model for Geographical Information Systems

This paper presents a development of the extended Cellular Automata (CA), a Voronoi-based CA, to model dynamic interactions among spatial objects. Cellular automata are efficient models for representing dynamic spatial interactions. A complex global spatial pattern is generated by a set of simple local transition rules. However, its original definition for a two-dimensional array limits its application to raster spatial data only. This paper presents a newly developed Voronoi-based CA in which the CA is extended by using the Voronoi spatial model as its spatial framework. The Voronoi spatial model offers a ready solution to handling neighbourhood relations among spatial objects dynamically. By implementing this model, we have demonstrated that the Voronoi-based CA can model local interactions among spatial objects to generate complex global patterns. The Voronoi-based CA can further model interactions among point, line and polygon objects with irregular shapes and sizes in a dynamic system. Each of these objects possesses its own set of attributes, transition rules and neighbourhood relationships. The Voronoi-based CA models spatial interactions among real entities, such as shops, residential areas, industries and cities. Compared to the original CA, the Voronoi-based CA is a more natural and efficient representation of human knowledge over space.     

Spatial chromatic model in high-dimensional spaces and the uniqueness of chromatic code: a new perspective of geographic entity–space relationship

The abstraction, representation, and computation of entity–space relationship are keystones of geographic information science (GIS). The newly proposed spatial chromatic tessellation (SCT) provides a novel model to explore this relationship. SCT has demonstrated a variety of potential applications in GIS, such as reasoning spatial topology, point pattern analysis, and Voronoi diagrams. This study aims to theoretically investigate SCT by focusing on two aspects: (1) extending SCT to higher dimensional spaces. Results show that cells missing in lower dimensional spaces are hidden in higher dimensional spaces; (2) exploring the uniqueness of chromatic codes, particularly the chromatic codes of 2-cell and 3-cell clusters: their codes are proved to be unique. In a mathematical perspective, the observed phenomena from the above two aspects bring some new thoughts into the first law of geography and spatial heterogeneity. Based on these new understandings of entity–space relationship, SCT is replaced by spatial chromatic model (SCM) in which spaces are created by entities themselves rather than by partitioning the space preexisted. This makes a change from an absolute geographic space to a relative geographic space.     

Revisiting John Snow's map: network-based spatial demarcation of cholera area

Dr. John Snow's cholera map is known as one of the pioneering examples of an epidemiology map, illustrating the spatial distribution of the victims from the cholera outbreak. This article revisits his map and expands on his attempt at visualizing the distribution of the victims by focusing on spatial demarcation using the sphere of influence along the street network by applying two analytical methods that are designed for analysis of network space. First, the article generates a network-based Voronoi diagram of the water pumps in the map that encompasses Snow's original version of the equidistance line that was drawn around a single pump. The article then presents a new, revised equidistance line that better reflects the circumstances around the time of the outbreak. This is followed by the construction of another set of boundaries derived by the application of a network-based clumping method. Comparing the demarcation lines produced with the two methods shows a sphere of influence which was unclear when using the network Voronoi diagram alone. Results from the analysis using the clumping method also confirm some of Snow's observations on the spatial distribution of the victims.     

新疆城市经济区划分及影响范围

明确且合理的城市经济区和空间组织是区域内各城市实现功能整合的前提,也是实现城乡统筹、有序协调发展的基础。通过因子分析方法来计算城市的中心性指数,在GIS软件的支持下,利用Voronoi图对地理空间进行分割,逐层递归,确定不同等级的中心城市,并构建城市断裂点理论与倍增加权Voronoi图相融合的GIS空间分析方法,划分城市经济影响范围。结果表明:第一,采用简单Voronoi的递归算法,依据城市的中心性指数,选出不同等级的中心城镇,很好地体现了中心职能在地域空间上的分布;第二,基于城市断裂点理论的倍增加权Voronoi图,不仅考虑了城市的实力,而且考虑了城市之间的空间邻近关系,并且很好地反映城市之间的空间竞争关系。第三,这种技术方法结合自然地理条件和行政区划等人文要素的定性分析,具有一定的科学性和现实意义。本文结合新疆独特的地理环境,探讨了新疆各级中心城镇和影响范围,旨在运用新技术和定量化的方法更科学地划分城市经济区,结合现有的城镇现状和发展战略,为未来新疆的城市与区域规划科学决策提供参考。     

Polygon Characterization With the Multiplicatively Weighted Voronoi Diagram*

Many landscape features are represented as polygons in GIS. This paper characterizes polygon shapes with the multiplicatively weighted Voronoi (MW‐Voronoi) diagram and improves its understanding. The MW‐Voronoi diagram's composition is implemented with topological overlay, growth simulation, and vertex calculation methods. The decomposition is done by reversing a polygon to MW‐Voronoi point pairs by segment. It is a new approach to record, characterize, and compare polygons with form and process. The implementation also serves as a geographic education and visualization tool. Applications of the methods are presented with precipitation, fire polygon, and population change data.     

Automated building generalization based on urban morphology and Gestalt theory

Building generalization is a difficult operation due to the complexity of the spatial distribution of buildings and for reasons of spatial recognition. In this study, building generalization is decomposed into two steps, i.e. building grouping and generalization execution. The neighbourhood model in urban morphology provides global constraints for guiding the global partitioning of building sets on the whole map by means of roads and rivers, by which enclaves, blocks, superblocks or neighbourhoods are formed; whereas the local constraints from Gestalt principles provide criteria for the further grouping of enclaves, blocks, superblocks and/or neighbourhoods. In the grouping process, graph theory, Delaunay triangulation and the Voronoi diagram are employed as supporting techniques. After grouping, some useful information, such as the sum of the building's area, the mean separation and the standard deviation of the separation of buildings, is attached to each group. By means of the attached information, an appropriate operation is selected to generalize the corresponding groups. Indeed, the methodology described brings together a number of well-developed theories/techniques, including graph theory, Delaunay triangulation, the Voronoi diagram, urban morphology and Gestalt theory, in such a way that multiscale products can be derived.