土地信息系统第章土地信息处理与分析技术(“多边形”相关文档)共80张

第5章土地信息处理与分析技术内容提要：1、土地信息录入后的处理2、土地信息的空间分析3、土地信息的空间查询15.1土地信息录入后的处理空间数据不论采用什么方法输入，都会有一些问题，如输入过程中出现意外的错误，输入数据与使用格式不一致，各种来源数据的比例尺、投影不统一，图幅间不匹配等。因此必须对空间数据进行处理，才能得到纯净、统一的数据文件，使存储的空间数据符合规范、标准的要求，满足使用和分析的需要。

2对空间数据的处理坐标变换图幅拼接生成拓扑关系数据核对3（一）坐标变换在地图数据录入完毕后，经常需要进行坐标变换，得到经纬度参照系下的地图。坐标变换的实质是建立两个平面点之间的一一对应的关系，它是空间数据处理的基本内容之一。坐标系一般有三种：用户坐标系规格化数据库坐标系设备坐标系4各种坐标系的定义用户坐标系：指地图采用的坐标系；规格化数据库坐标系：在数据库系统中定义的坐标系；设备坐标系：每一种图形设备都有独特的坐标系，它们使用的坐标都是设备的相对坐标。

5图5-1设备坐标系

61.设备相对坐标到用户坐标的变换一般在图形数字化前后进行，它是利用点的高斯坐标和图幅角点采集坐标建立采集坐标到高斯坐标的变换关系，并利用这些变换参数将图幅中全部采集坐标变换为高斯坐标。2.用户坐标到数据库坐标的变换这个变换用于图形数据的入库。3.数据库坐标到用户坐标的变换这个变换用于空间数据检索，它是2的逆变换。坐标变换的类型74.用户坐标到设备坐标的变换这个变换用于图形显示或绘图仪绘图。用于图形显示时，只要将高斯坐标原点平移至图幅左上角，将坐标顺时针旋转90°，并考虑两种坐标的变换比例，即可实现由高斯坐标到屏幕坐标的变换。5.数据库坐标到屏幕坐标的变换和屏幕坐标到数据库坐标的变换这两个变换是用来进行人机交互编辑并将编辑好的图形数据送回数据库。8

土地信息系统中，如同计算机辅助设计一样存在着对图形缩放、平移、旋转、投影等一系列图形的几何变换问题。这些变换的实质是对组成图形各点进行坐标变换。

图形的平移变换只改变图形位置，不改变形状和大小。图形的旋转变换用来改变图形的视角，变换后图形各部分间线性关系和角度关系不变，变换后直线长度不变。图形的比例变换用来改变图形大小和形状。投影变换比较复杂,是从一种投影方式转成另一种投影方式，其过程是一系列复杂的数学运算。图形的几何变换9图5-2图形的几何变换10（二）图幅接边

在对底图进行数字化时，由于图幅比较大或者使用小型数字化仪，难以将研究区域的底图以整幅的形式来完成，这时需要将整个图幅划分成几部分分别数字化。数字化完成后，又要把相邻图幅拼起来，这就是图幅的接边问题。

接边时应考虑常见的误差不因接边而明显扩大，接边处理的过程如下图：11(a)拼接前；(b)拼接中的边缘不匹配；(c)调整后的拼接结果

图5-3图幅接边过程121、识别和检索相邻图幅

将待拼接的图幅数据按图幅进行编号，编号通常取２位数，其中个位数表示拼幅时横向顺序，十位数表示拼幅时的纵向顺序。在横向拼幅时将十位数相同图，按个位数顺序拼接；纵向拼幅时将个位数相同图，按十位数顺序拼接。图5-4图幅编号及图幅边缘数据的提取范围图幅接边的步骤：132、逻辑一致性的处理

由于人工操作的失误，两个相邻图幅的空间数据库在接合处可能出现逻辑不一致，如一个多边形在一幅图层中具有属性A，而在另一幅图层中属性为B；同一线状地物(铁路，公路，单线河等)在相邻图幅中各部分的属性不一致。此时，必须使用交互编辑的方法，使两相邻图幅各部分的属性相同，取得逻辑一致性。14图幅接边前图幅接边后图5-5逻辑一致性的处理

或153、相邻图幅边界点坐标数据的匹配

常用方法有：①平均法。它是指将图形两边待接点的坐标均值作为接边后点的坐标。该方法简单易行，适用于接边误差在精度允许范围内的各种直线，多义线类的接边处理；②强制法。它是把一条待接边的待接点强制附和到另一条待接边的待接点上。该方法主要用于用户能明显判断出那一条待接边比另一条待接边更准确、可靠、适合交互式的接边处理；

考虑到接边工作的复杂性，一般主要采用强制法进行接边。16图5-6图幅接边前后对照174、相同属性多边形公共边界的删除图5-7相同属性多边形公共边界的删除与属性合并将相同属性的两个或多个相邻图斑组织成一个图斑，即消除公共边界，并对共同属性进行合并。对于合并后的属性表，除多边形的面积和周长重新计算外，其余属性保留其中之一图斑的属性即可。18图形常见错误建立拓扑关系多边形拓扑关系的建立拓扑生成点线拓扑关系的建立（三）建立拓扑关系19拓扑生成

在进行图形数字化错误修改之前，首先要建立拓扑关系，以正确判别地物之间的空间关系。构造拓扑关系可以生成弧的交点、定义构成多边形的弧，将标示号关联到相应的多边形，生成多边形。这样可以帮助我们确定数字化数据中存在的某些错误。20图形常见错误伪节点使一条完整的线变成两段（如图），造成伪节点的原因常常是没有一次录入完毕一条线。一般来说一条线不应该被分割为两条线，但在现实世界中存在两条不同属性的线相连，因此伪结点并非全是错误的。1.伪节点（PseudoNode）21仅与一条弧相连的结点称为悬挂点。与悬挂点相连的弧称为悬挂弧。通常有两种情况：①未及（undershoot）：一种数字化错误类型，导致弧段之间存在缝隙而未接合。②过伸（overshoot）：一种数字化错误类型，它导致弧段过长。2.悬挂点（DangleNode）22碎屑多边形（如图）一般由于重复录入引起，由于前后两次录入同一条线的位置不可能完全一致，造成了“碎屑”多边形。另外，由于用不同比例尺的地图进行数据更新，也可能产生“碎屑”多边形。3.“碎屑”或“条带”多边形（SliverPolygon）23不正规的多边形是由于输入线时，点的次序倒置或者位置不准确引起的。在进行拓扑生成时，同样会产生“碎屑”多边形。

4.不规则的多边形（WeirdPolygon）245.多边形标识点（label）错误

多边形标示点错误通常有两种情况：一是多边形无标示点；一是多边形有多个标示点，如右图中(两个十字点)。

25图5-8数据化常见错误26

在图形修改完毕之后，就意味着可以建立正确的拓扑关系，拓扑关系可以由计算机自动生成，目前大多数GIS软件也都提供了完善的拓扑功能。建立拓扑关系27点线拓扑关系的建立弧段起结点终结点a1N1N2a2N2N3弧段—结点表结点弧段N1a1N2a1a2N3a2结点—弧段表弧段起结点终结点a1N1N2a2N2N3a3N2N4结点弧段N1a1N2a1a2a3N3a2N4a3弧段起结点终结点a1N1N2a2N2N3a3a4N2N4N4N3结点弧段N1a1N2a1a2a3N3a2a4N4a3a4N1N2N3N1N2N3N1N2N3N4N4a1a2a2a1a3a1a2a3a428多边形拓扑关系的建立（1）链的组织：找出在链的中间相关(图5-10(1))而不是在端点相关(图5-10(2))的情况，自动切成新链；把链按一定顺序存储并编号。图5-9链的断开要求29（2）结点匹配：把一定限差内的链的端点作为一个结点，去掉在数据输入过程中形成的结点过交、悬挂线等现象，其坐标值取多个端点的平均值。图5-10结点匹配示意图30（3）检查多边形是否闭合：检查多边形是否闭合可以通过判断一条链的端点是否有与之匹配的端点来进行。多边形不闭合的原因：1）由于结点匹配限差的问题，造成应匹配的端点未匹配；2）由于数字化误差较大，或数字化错误；3）链本身就是悬挂链，不需参加多边形拓扑。图5-11多边形闭合的端点匹配31顺时针方向构多边形：指多边形是在链的右侧。最靠右边的链：指从链的一个端点出发，在这条链的方向上最右边的一条链。

（4）建立多边形基本概念：32建多边形的算法与实例:（1）从Pl结点开始，起始链定为P1P2；从P2点算起，P1P2最右边的链为P2P5；从P5算起，P2P5最右边的链为P5P1。所以，形成的多边形为P1P2P5P1;（2）从P1结点开始,以P1P5为起始链，形成的多边形为P1P5P4Pl;（3）从Pl开始，以PlP4为起始链形成的多边形为PlP4P3P2P1;（4）这时Pl为结点的所有链均被使用了两次，因而转向下一个结点P2，继续进行多边形追踪，直至所有的结点取完。共可追踪出五个多边形，即Al、A2、A3、A4、A5。

顺序取一个结点作为起始结点、过该结点的任一条链为起始链结点取完？Y取该链另一结点，找结点上靠此链最右边的链，作为下一条链N是否回到起点？形成多边形，记录YN取起始结点的、刚才形成多边形的最后一边作为新的起始链链是否被用过两次，即已经成为两个多边形的边YN结束33（5）岛的判断岛的判断即指找出多边形互相包含的情况。

图5-12岛的判断与数据组织34

如果数据是分层多次输入到数据库中，当不同图层数据叠合在一起时，往往会出现该重合的地方不重合的问题。如将道路跟同样位置的地类边界叠合在一起，往往会出现该重合的地方不重合的问题，也可能出现一些很细小的多边形。（四）数据核对防止这一现象的一种方法是在地图数字化跟踪之前，规定不同要素数字化的基本顺序，并且将数字化顺序靠后的其他图层上表示不同要素，而在空间位置上是和数字化顺序靠前的图层的某些要素是重合的点、线不用重复数字化跟踪。35对属性数据的处理

属性数据较为规范，适应于采用表格形式表示，所以许多LIS都采用关系数据库管理系统管理属性数据。文件管理结构设计数据编辑数据查询统计报表创建文件打开文件添加文件存贮文件打印文件图示分析退出系统设计新结构修改老结构拷贝结构删除结构合并结构添加数据插入数据删除数据修改数据合并数据库查询结构信息行拷贝数据列拷贝数据浏览数据显示一条记录显示几列属性条件查询移动数据统计总数据计算平均值统计最大值报表设计多表连接查询打印报表365.2土地信息的空间分析缓冲区分析叠加分析网络分析37缓冲区分析概念

缓冲区是指对点、线或面实体，按指定的条件，在其周围建立一定宽度范围的空间区域作为分析对象，这个区域称为缓冲区。

(1)点的缓冲区(2)线的缓冲区(3)面的缓冲区

图5-13不同类型实体的缓冲区38缓冲区的建立线的缓冲区面的缓冲区点的缓冲区原理：以点为中心，以一定距离为半径绘圆；原理：分别对每个顶点和每条线生成缓冲区，然后对这些缓冲区多边形进行叠置操作；原理：首先生成多边形周长的缓冲区，然后与原始多边形进行叠置操作，多边形缓冲区有内外侧之分。外缓冲区在面状地物的外围形成缓冲区,内缓冲区则在面状地物的内侧形成缓冲区。39图5-14多边形的内外缓冲区40概念叠加分析，就是将具有相同坐标系统的多个空间要素对象的数据层进行叠加，产生一个新数据层面，该数据层面综合了原来两层或多层要素所具有的属性特征。叠加分析41视觉信息叠加视觉信息叠加是将不同侧面的信息内容叠加显示在结果图件或屏幕上，以便研究者判断其相互空间关系，获得更为丰富的空间信息。土地信息系统中视觉信息叠加包括以下几类：点状图，线状图和面状图之间的叠加显示面状图区域边界之间或一个面状图与其他专题区域边界之间的叠加遥感影象与专题地图的叠加专题地图与数字高程模型（DEM）叠加显示立体专题图说明：视觉信息叠加不产生新的数据层面，只是将多层信息复合显示，便于分析42确定一幅图（或数据层）上的点落在哪一幅图（或数据层）的哪个多边形中，这样就可以给相应的点增加新的属性内容。核心算法为判断点是否在多边形内。点与多边形的叠加

C

4

B

3

D

2

A

1属性3多边形号属性2属性1点号1342ABCD4344把一幅图(或一个数据层)中的多边形特征加到另一幅图(或另一个数据层)的线上。线与多边形叠加的算法就是线的多边形裁剪。D34B23A22A11多边形号原线号线号线与多边形的叠加12314ABCD234145例如：如果线状图层为道路网可以得到每个多边形内的道路网密度多边形内部的交通流量进入、离开各个多边形的交通量相邻多边形之间的相互交通量与多边形叠加多边形将穿过它的所有道路打断成弧段如果线状图层为河流与多边形叠加多边形将穿过它的所有河流打断成弧段可以查询任意多边形内的河流长度计算它的河流密度46是将同一地区，同一比例尺的两组或更多的多边形要素的数据层进行叠置，会出现两种成果形式：一种是根据两组多边形边界的交点来建立具有多重属性的多边形，称为合成叠置；另一种是进行多边形范围的属性特征的统计分析，称为统计叠置。多边形与多边形的叠加47合成叠置：是根据两组多边形边界的交点来建立具有多重属性的多边形。合成叠置得到一张新的叠置图，产生了许多新的多边形，每个多边形内都具有两种以上的属性，通过区域多重属性的模拟，寻找和确定同时具有几种地理属性的分布区域。合成叠置48图5-15合成叠加产生的碎屑多边形

49

根据叠加结果最后欲保留空间特征的不同要求，一般的GIS软件都提供了三种类型的合成叠置操作：多边形的不同叠置方式

50Union(并)Identity(并集被第一个输入层边界裁切后剩余部分——叠和)Intersect（交）

ArcGIS中空间多边形叠置51图5-16多边形合成叠加分析52统计叠置统计叠置：是进行多边形范围的属性特征的统计分析。统计叠置的目的是精确地计算一种要素（例如，土地利用）在另一种要素（例如，行政区域）的某个区域多边形内的分布状况和数量特征（包括拥有的类型数、各类型的面积及其所占总面积的百分比等等），或提取某个区域范围内某种专题内容的数据。统计叠置的结果为统计报表或列表输出。53栅格图层叠加

栅格图层叠加有时也称为栅格数据的信息复合,它是指不同层面的栅格数据逐网格按一定的数学法则或逻辑判断进行运算,从而得到新的栅格数据系统的方法。

栅格图层叠加方法：数学运算法：是指不同层面的栅格数据逐网格按一定的数学法则进行运算，从而得到新的栅格数据系统的方法;54图5-17以函数关系进行叠加55逻辑关系分析法：用逻辑表达式来分析处理重合点的非几何特性之间的逻辑关系，实现对空间数据复合、提取、删除等操作。通常是按照两个逻辑子集给定的条件进行逻辑交运算、逻辑并运算、逻辑非运算等。ACBACBACB56图5-2图形的几何变换资源分配模型可用来计算中心地的等时区、等交通距离区、等费用距离区等。悬挂点（DangleNode）用于图形显示时，只要将高斯坐标原点平移至图幅左上角，将坐标顺时针旋转90°，并考虑两种坐标的变换比例，即可实现由高斯坐标到屏幕坐标的变换。用户坐标到设备坐标的变换这个变换用于图形显示或绘图仪绘图。栅格图层叠加有时也称为栅格数据的信息复合,它是指不同层面的栅格数据逐网格按一定的数学法则或逻辑判断进行运算,从而得到新的栅格数据系统的方法。一是多边形有多个标示点，如右图中(两个十字点)。属性数据较为规范，适应于采用表格形式表示，所以许多LIS都采用关系数据库管理系统管理属性数据。线线查询如与某条河流相连的支流有哪停靠点（Stops）：网络路线中资源装、卸的结点点位合成叠置：是根据两组多边形边界的交点来建立具有多重属性的多边形。通常是按照两个逻辑子集给定的条件进行逻辑交运算、逻辑并运算、逻辑非运算等。说明：视觉信息叠加不产生新的数据层面，只是将多层信息复合显示，便于分析ArcView软件中属性查图形的表达方式图形的旋转变换用来改变图形的视角，变换后图形各部分间线性关系和角度关系不变，变换后直线长度不变。该方法简单易行，适用于接边误差在精度允许范围内的各种直线，多义线类的接边处理；例：有土壤厚度（大于和小于50厘米）和土壤类型（红壤和其他类型）两个二值化图层，不同的逻辑运算结果如下：逻辑交运算：结果是将土层厚度大于50厘米，且土壤为红壤的土壤单元显示出来；逻辑并运算：结果将土层厚度大于50厘米，或者土壤为红壤的土壤单元显示出来；逻辑非运算：结果是将土层厚度大于50厘米，但土壤不是红壤的土壤单元显示出来。57网络分析网络：是一个由点和线的二元关系构成的系统，通常用来描述某种资源或物质在空间上的运动。网络分析：是指依据网络拓扑关系（结点与弧段拓扑、弧段的连通性），通过考察网络元素的空间及属性数据，以数学理论模型为基础，对网络的性能特征进行多方面研究的一种分析计算。网络分析的根本目的是研究、筹划一项网络工程如何安排，并使其运行效果最好，如一定资源的最佳分配，从一地到另一地的运输费用最低等。58结点（Nodes）：网络中任意两条线段的交点连通路线或链（Links）：连结两个结点的弧段要素，是网络中资源运移的通道,如街道、河流、水管等转弯（Turns）：资源运移方向从一个链上经结点转向另一个链停靠点（Stops）：网络路线中资源装、卸的结点点位

中心（Centers）：网络中具有接收或发放资源能力，且位于节点处的设施,如水库、商业中心、电站等障碍（Barriers）：禁止网络中链上流动的节点，即资源不能通过的节点网络的基本要素59图5-18网络的基本要素

60主要网络分析功能（1）最短路径分析

最短路径是在网络中从起点经一系列特定的结点至终点的资源运移的最佳路线，即阻力最小的路径。在最短路径选择中，两点之间的距离可以定义为实际的距离，也可定义为两点间的时间、运费、流量等，还可定义为使用这条边所需付出的代价。61（2）地址匹配

地址匹配实质是对地理位置的查询，它涉及到地址的编码。地址匹配与其它网络分析功能结合起来，可以满足实际工作中非常复杂的分析要求。如物流配送，需要将货物送到多家单位。此时需将路径分析和地理编码结合起来使用，可先通过地理编码进行地址查询，获得各单位的地理位置，再利用最短路径方法确定最短送货线路。62（3）资源分配

资源分配模型中的网络主要由中心点（分配中心）组成，有两种分配方式：一是由分配中心向四周输出资源；另一种是从四周向分配中心集中资源。这种分配功能可以解决资源的有效流动和合理分配。资源分配模型可用来计算中心地的等时区、等交通距离区、等费用距离区等。可用来进行城镇中心、商业中心或港口等地的吸引范围分析，以用来寻找区域中最近的商业中心，进行各种区划和港口腹地的模拟等。

63图形查询属性属性查询图形

空间查询图文互访空间关系查询

地址匹配查询

是按一定的要求对土地信息系统所描述的空间实体及其空间信息进行访问，从众多的空间实体中挑选出满足用户要求的空间实体及其相应的属性。

空间信息查询5.3土地信息的空间查询64MapInfo软件中点目标的几何参数查询（一）图形查询属性65MapInfo软件中线目标的几何参数查询66Mapinfo软件中面状目标的几何参数查询67ArcView软件中图形查属性的表达方式68ArcView软件中属性查图形的表达方式（二）属性查询图形69（三）空间关系查询

面面查询

如与某个多边形相邻的多边形有哪些面线查询

如某个多边形的边界有哪些线面点查询

如某个多边形内有哪些点状地物线面查询

如某条线经过（穿过）的多边形有哪些，某条链的左、右多边形是哪些70线线查询

如与某条河流相连的支流有哪些，某条道路跨过哪些河流。线点查询

如某条道路上有哪些桥梁，某条输电线上有哪些变电站。点面查询

如某个点落在哪个多边形内。点线查询

如某个结点由哪些线相交而成。

71（四）地址匹配查询

根据街道地址来查询事物的空间位置和属性信息是地理信息系统特有的一种查询功能，这种查询利用地理编码，输入街道门牌号码，就可知道大致的位置和所在的街区。它对空间分布的社会、经济调查和统计很有帮助，只要在调查表中添了地址，地理信息系统可以自动地从空间位置的角度来统计分析各种经济社会调查资料。

72(五)缓冲区查询

缓冲区查询与前面介绍的缓冲区分析有一点差别，缓冲区查询不对原有图形进行切割，只是根据用户需要给定一个点缓冲、线缓冲或面缓冲的距离，从而形成一个缓冲区的多边形，再根据多边形检索的原理，检索出该缓冲区多边形内的空间地物。

73距黄河150公里范围内的主要城市74

Select需显示的属性项

From属性表

Where条件

or