《地理信息系统概论.PPT》由会员分享，可在线阅读，更多相关《地理信息系统概论.PPT（128页珍藏版）》请在金锄头文库上搜索。
1、2024/8/2311第三章第三章 空间数据处理空间数据处理l空间数据处理是GIS的主要功能之一。l空间数据处理是针对空间数据本身完成的操作，不涉及内容的分析。因此，空间数据处理又称为空间数据形式的操作。2024/8/232 2l空间数据处理包括：数据变换：是指数据从一种数学状态到另一种数学状态的变换。包括几何纠正和地图投影转换等，以实现空间数据的几何配准。数据重构：是指数据从一种格式到另一种格式的转换，包括空间数据在结构、格式和类型上的统一，多源和异构数据的连接与融合。数据提取：是指对数据进行某种条件的取舍。包括类型提取、窗口提取、空间内插等，以适应不同用户对数据的特定要求。2024/8/2
2、33 3第一节第一节 空间数据的变换空间数据的变换l空间数据变换：即空间数据坐标系的变换。其实质是建立两个坐标系坐标点之间的一一对应关系。包括几何纠正和投影转换等。他们是空间数据处理的基本内容之一。l空间数据变换的目的：对数字化原图的数据进行坐标转换和变形误差的纠正；不同来源的地图进行地图投影转换与地图比例尺的统一。2024/8/234 4一、几何纠正一、几何纠正l原因：在图形编辑中，只能消除数字化产生的明显误差，而图纸变形产生的误差难以改正，因此要进行几何纠正。l目的：是为了实现对数字化数据的坐标系转换和图纸变形误差的纠正。l（几何纠正）分类：仿射变换相似变换二次变换等2024/8/235
3、5仿射变换仿射变换l仿射变换：可以对坐标数据在x和y方向进行不同比例的缩放，同时进行扭曲、旋转和平移。l仿射变换是使用最多的一种几何纠正方式，只考虑到x和y方向上的变形l仿射变换的特性：直线变换后仍为直线平行线变换后仍为平行线不同方向上的长度比发生变化2024/8/236 6l仿射变换可以不同程度地缩放、倾斜、旋转和平移数据。下图说明了这四种可能的变化。2024/8/237 7仿射变换仿射变换l设x，y为数字化仪坐标，X,Y为理论坐标， m1和m2为地图横向和纵向的长度变化比例，两坐标系的夹角为，数字化仪原点O相对于理论坐标原点平移了A0、B0，则根据图形变换原理，得出坐标变换公式为：l令l则
4、上式简化为2024/8/238 8l对于仿射变换，只需知道不在同一直线上的三对控制点的坐标及其理论值，就可求得待定系数。但在实际使用时，往往利用4个以上的点进行纠正，利用最小二乘法处理，以提高变换的精度。 2024/8/239 9二、地图投影及其转换二、地图投影及其转换一、地图投影的基本原理l1、为什么需要地图投影：将地球椭球体上的空间信息表现到平面地图上，或用GIS的地图图形显示出来，就必须采用某种数学法则，使空间信息在地球表面上的位置和地图平面位置一一对应起来。l2、建立地图投影的目的：采用某种数学法则，使空间信息在地球表面上的位置和地图平面位置一一对应起来，以满足地图制图的要求。2024
5、/8/231010l理解地图投影如何改变空间属性的一种简便方法： 观察光穿过地球投射到表面（称为投影曲面）上。 想像一下，地球表面是透明的，其上绘有经纬网。用一张纸包裹地球。位于地心处的光会将经纬网投影到一张纸上。现在，可以展开这张纸并将其铺平。纸张上的经纬网形状与地球上的形状不同。 地图投影使经纬网发生了变形。2024/8/231111l地球表面是一个不规则曲面，即使把它当做椭球或球表面，数学上也还是一种不能展开成平面的曲面，即把这样的曲面摊开成平面就会发生破裂或褶皱。l理想的地图应该使制图区域不破裂，并保持图形轮廓的完整一致，其边界邻区也应是延续的、对接无误的。2024/8/231212l
6、地图投影在实践中是采用经纬度表示的地球表面上的点位与平面直接坐标或极坐标表示的平面上的点位，通过一定的数学关系式建立起对应的联系。其中，x,y为地球表面经纬度坐标； X，Y为地图平面坐标； f1、f2为对应的数学法则。2024/8/231313二、地图投影的类型l地图投影的主要任务：研究解决曲面如何转化到平面的问题；分析如何解决地图投影中始终存在的变形问题。l地图投影中的变形主要有三种：角度变形、面积变形、长度变形l按照地图投影变形性质，一般把地图投影分为三类：等角投影、等面积投影和任意投影2024/8/231414l投影变形在质和量上的差别： 假设在地球表面上有一个微小的圆形，称为微分圆微分
7、圆。不同投影由于有不同的投影条件，这个微分圆投影后在地图平面上不一定仍为圆形，也可能是一个椭圆，将其称为变形椭圆变形椭圆。 变形椭圆的形状和大小能确切反映出投影变形在质和量上的差别。2024/8/231515根据变形椭圆的原理分析根据变形椭圆的原理分析l等角投影（或称为正形投影）：特点：投影后形状无变化，大小发生变化。优点：等角投影保留局部形状。缺点：由一些弧线围起来的区域将在此过程中发生巨大变形。地图投影无法保留较大区域的形状。l等距离投影：特点：投影后形状发生变化，距离不变。优点：等距地图保留某些点间的距离。任何投影都无法在整幅地图中正确保持比例不变。 多数情况下，地图上总会存在一条或多条
8、这样的线：比例沿着这些线将正确地保持不变。多数等距投影都具有一条或多条这样的线：在此类线中，地图上线的长度（按地图比例尺计算）与地球上同一条线的长度相同，无论它是大圆还是小圆，是直线还是曲线。此类距离被视为真实距离。例如，在正弦投影中，赤道和所有纬线就是其真实长度。在其他等距投影中，赤道和所有经线具有真实长度。而其他投影（例如，两点等距离）仍会显示地图上一点或两点与相隔点间的真实比例。l请记住，任何投影都不能实现地图上的所有点是等距离的。2024/8/231616l等面积投影：特点：投影后面积无变化。优点：等积投影保留所显示要素的面积。 形状、角和比例等其他属性将发生变形。 在等积投影中，经线
9、和纬线可能不垂直相交。 有些情况下，尤其是较小区域的地图，形状不会明显变形,且很难区分等积投影和等角投影，除非加以说明或进行测量。l任意投影：按其性质既不属于等角又不属于等积。l注意：地球表面上的长度、面积和角度经过投影后一般均会发生变化.2024/8/231717l按投影面的形状分为 圆锥投影：投影面为圆锥面 圆柱投影：投影面为圆柱面 方位投影：投影面为平面2024/8/231818地图投影分类地图投影分类l地图投影按投影面与地球的相对位置关系分为：正轴投影：投影面的旋转轴与地球旋转轴重合；横轴投影：投影面的旋转轴与地球旋转轴垂直；斜轴投影：投影面的旋转轴与地球旋转轴既不垂直又不重合；l了解
10、：正轴切圆锥投影和地球表面相切于一条纬线，称为标准纬线。标准纬线投影后保持长度不变。正轴割圆锥投影和地球表面相割于两条标准纬线。l按投影面和地球的空间逻辑关系可以分为： 相切投影 相割投影2024/8/231919l地图主比例尺（或称为普通比例尺）计算地图投影或制作地图时，必须将地球要一定比例缩小表示到平面上，这个比例称为地图主比例尺，或称为普通比例尺。l局部比例尺：由于投影中必定存在某种变形，地图仅能在某些点或线上保持比例尺，其余位置的比例尺都与主比例尺不相同，即大于或小于主比例尺。这个比例尺被称为局部比例尺。l一般地图上注明的比例尺是主比例尺，而对用于测量长度的地图要采用一定的方式设法表示
11、出该图的局部比例尺。这就是在大区域小比例尺地图（小于1:1 000 000）上常见的图解复式比例尺。2024/8/232020三、地理信息系统常用的地图投影l1、高斯-克吕格投影： 是由德国数学家、物理学家和天文学家高斯于19世纪20年代拟定。德国大地测量学家克吕格于1912年对投影公式加以补充，故名高斯-克吕格投影是等角横切椭圆柱投影。 也称为横轴墨卡托，此投影与墨卡托投影类似，不同之处在于圆柱是沿子午线而非赤道纵向排列。通过这种方法生成的等角投影不会保持真实的方向。此投影最适合于南北分布的地块。2024/8/232121l从几何意义上来看，就是假想用一个椭圆柱横向套在地球椭球外面，并于某一
12、子午线相切（此子午线称为中央子午线或中央经线），椭圆柱的中心轴位于地球椭球的赤道上。再按高斯-克吕格投影所规定的条件，将中央经线两侧一定经差范围内的经纬线投影到椭圆柱面上，并将此椭圆柱面展为平面，即得高斯-克吕格投影l1949年后，高斯-克吕格投影被确定为我国地形图系列中1:50万、1:20万、1:10万、1:5万、1:2.5万、1:1万及更大比例尺的数学基础。2024/8/232222l高斯-克吕格投影所规定的条件：中央经线和赤道投影后为互相垂直的直线，且为投影的对称轴；投影具有等角性质；中央经线投影后保持长度不变。2024/8/232323l高斯-克吕格投影的特点：中央经线上没有任何变形，
13、满足中央经线投影后保持长度不变的条件；除中央经线上的长度比为1外，其他任何点上长度比均大于1；在同一条纬线上，离中央经线越远，变形越大，最大值位于投影带的边缘；在同一条经线上，纬度越低，变形越大，变形最大值位于赤道上。投影属于等角性质，故没有角度变形，面积比为长度比的平方；长度比的等变形线平行于中央子午线。2024/8/232424地图上坐标分带地图上坐标分带l高斯-克吕格投影的最大变形处为各投影带在赤道边缘处，为了控制变形，我国地形图采用分带方法：即将地球按一定间隔的经差（6或3）划分为若干个相互不重叠的投影带，各带分别投影。l6分带1:2.5万至1:50万的地形图均采用6分带方案，即从格林
14、尼治零度经线起算，每6为一个投影带，全球共分为60个投影带。我国领土位于东经72到136之间，共包含11个投影带（13-23带）l3 分带：11万的地形图采用3度分带，从东经1.5度的经线开始，每隔3度为一带，用1，2，3，表示，全球共划分120个投影带，即东经1.5 4.5度为第1带，其中央经线的经度为东经3度，东经4.57.5度为第2带，其中央经线的经度为东经6度。地形图上公里网横坐标前2位就是带号，例如：15万地形图上的坐标为（18576000， 293300），其中18即为带号。2024/8/232525当地中央经线经度的计算当地中央经线经度的计算l六度带中央经线经度的计算：当地中央经
15、线经度6当地带号3，例如：地形图上的横坐标为18576000，其所处的六度带的中央经线经度为：6183105。l三度带中央经线经度的计算：中央经线经度3当地带号。l如何判断地形图采用哪种分带一个好记的方法:在中华人民共和国陆地范围内,坐标(Y坐标,8位数,前两位是带号)带号小于等于23的肯定是6度带，大于等于24的肯定是3度带。 2024/8/232626l只知道经纬度时中央经线的计算 将当地经线的整数部分除以6，再取商的整数部分加上1。再将所得结果乘以6后减去3，就可以得到当地的中央经线值。如1061500，用106/6取整得17，（17+1）*6-3=105，即当地的中央经线值为105。毕
16、节，位于东经1053610643，北纬26212746之间， 计算毕节当地的中央经线以及所处的三度带和六度带分别是多少？ 2024/8/232727l2、墨卡托投影（Mercator）是一种等角正切圆柱投影是荷兰地图学家墨卡托在1569年拟定。最初设计该投影的目的是为了精确显示罗盘方位，为海上航行提供保障。此投影的另一功能是能够精确而清晰地定义所有局部形状。许多 Web 制图站点都使用基于球体的墨卡托投影。球体半径等于 WGS 1984 长半轴的长度，即 6378137.0 米。有两种用于仿真 Web 服务所用墨卡托投影的方法。2024/8/232828l基本思想：假设地球被围在一个空的圆柱里
17、，其标准纬线与圆柱相切，然后假想地球轴心有一盏灯，把球面上的形状投影到圆柱面上，再把圆柱面展开。即得墨卡托投影。l特点：角度无变化，每个点向各方向的长度比相等。保持了方向和相互位置关系的正确。2024/8/232929l局限性局限性在墨卡托投影上无法表示极点。可以对所有经线进行投影，但纬度的上下限约为 80 N 和 80 S。大面积变形使得墨卡托投影不适用于常规地理世界地图。l用途和应用用途和应用标准海上航线图（方向）。其他定向使用：航空旅行、风向、洋流。等角世界地图。此投影的等角属性最适合用于赤道附近地区，例如，印尼和太平洋部分地区。2024/8/2330303、UTM投影l通用横轴墨卡托投
18、影（Universal transverse Mercator projection，UTM投影）等角横轴割圆柱投影与高斯-克吕格投影（等角横轴切圆柱投影）相似。属于等角横轴割圆柱投影，圆柱割地球于两条等高圈上，投影后这两条割线上没有变形。lUTM投影可以改善高斯-克吕格投影。两条割线即等高圈上没有任何变形，离开这两条割线越远则变形越大，在两条割线以内长度变形为负值，在两条割线以外长度变形为正值。2024/8/233131lUTM投影分带方法与高斯-克吕格投影相似：从西经180起，每个经差6，自西向东分带，将地球划分为60个投影带，并且每个投影带当地的中央经度的计算公式：当地中央经线经度=6当
19、地带号 例如：地形图上的横坐标为20 345，其所处的6带的中央经线经度=6适用于1:25 000万和1:50 000万地形图。2024/8/233232l用途和应用用途和应用用于比例尺为 1:100,000 的美国地形地图方格。许多国家/地区使用基于现行官方地理坐标系的地方 UTM 区域。前苏联的大比例尺地形制图。我国卫星影像资料常采用。2024/8/233333l4、兰勃特投影兰勃特（Lambert）等角投影在双标准纬线下是一个等角正轴割圆锥投影。由德国数学家兰勃特在1772年拟定。兰勃特等角投影后纬线为同心圆弧，经线为同心圆半径。此投影是最适用于中纬度的一种投影，其描绘形状比描绘面积更准
20、确。美国国家平面坐标系对所有具有较大东西范围的区域均使用此投影。l我国1:100万地形图采用兰勃特投影，其分幅原则与国际地理学会规定的全球统一使用的国际百万分之一地图投影保持一致。l局限性局限性最适合主要为东西范围并且位于北纬或南纬的中间纬度的区域。总的纬度范围不应超过 35。2024/8/233434l5、阿尔伯斯投影阿尔伯斯投影（Albers）是一种正轴等面积割圆锥投影，与兰勃特投影属于同一投影族。这种圆锥投影使用两条标准纬线，相比使用一条标准纬线的投影可在某种程度上减少变形。这种投影最适合于东西方向分布的大陆板块，而不适合南北方向分布的大陆板块。l我国大部分省区图以及大多数这一比例尺的地
21、图也多采用兰勃特投影和阿尔伯斯投影2024/8/233535四、地图投影转换四、地图投影转换l地图投影转换是地图投影和地图编绘的一个重要组成部分当GIS使用的数据取自不同地图投影的图幅时，需要将一种投影的数字化数据转换为所需要投影的坐标数据。l地图投影转换主要研究：从一种地图投影变为另一种地图投影的理论和方法。其实质是建立两平面之间点的一一对应关系。2024/8/233636l投影转换的方式主要有两种：（1）正解变换：通过建立一种投影变换为另一种投影的严密或近似的解析关系式，直接由一种投影的数字化坐标x，y变换到另一种投影的直角坐标X，Y.即：直接求出两种投影点的直角坐标关系式。 （2）反解变
22、换：由一种投影的坐标反解出地理坐标（ x，y B，L），然后将地理坐标带入另一种投影的坐标公式中（B，LX,Y），从而实现由一种投影的坐标到另一种投影坐标的变换（ x，y X，Y ）。2024/8/233737根据转换的方法不同，投影转换可分为：2024/8/2338l（1）解析变换：在获得原数据与新数据的投影公式的情况下，运用正解法或反解法求出原投影坐标x，y与新投影坐标X，Y之间的精确表达式。即：解析出转换表达式38l（2）数值变换：在原数据投影公式未知时（包括投影常数难判别时），或不易求得原数据和新数据两投影间解析关系式的情况下，可以采用多项式来拟合他们间的关系，即利用两投影间已知的若干
23、离散点（也称为共同点）xi，yi和Xi，Yi的数值，根据数值逼近的理论和方法，主要包括插值法、有限差分法、有限元法、待定系数法等，来建立两投影间的关系式的方法。它是地图投影变换中在理论上合实用上的一种较通用的方法。2024/8/233939l（3）解析-数值变换：解析变换与数值变换的一种综合运用。当新数据投影已知，而原数据投影公式不知道的情况下，可以先由原数据投影反解出某些投影点的地理坐标，再代入已知的新数据投影公式中进行计算，便可实现两投影间的变换。2024/8/234040第二节第二节 空间数据结构的转换空间数据结构的转换l回顾矢量数据和栅格数据的基本特点及其优缺点l矢量数据：是面向实体的
24、结构，即对于每一个具体的实体都直接赋有位置和属性的信息，以及目标之间的拓扑关系说明。l栅格数据：是面向位置的结构，平面空间上的任意一点都直接联系到某一个或某一类实体。但对于某一个具体的实体又没有直接聚集所有信息，只能通过遍历栅格矩阵逐一查找，不能完整地建立实体之间的拓扑关系。2024/8/234141矢量与栅格数据结构的比较矢量与栅格数据结构的比较2024/8/234242应用原则应用原则l矢量与栅格数据的应用原则数据采集采用矢量数据结构，有利于保证空间实体的几何精度和拓扑特性的描述；空间分析则主要采用栅格数据结构，有利于加快系统数据的运行速度和分析应用的进程。l在数据处理阶段，经常要进行两种
25、数据结构的相互转换。2024/8/234343一、由矢量向栅格的转换一、由矢量向栅格的转换l栅格化：将点、线或面的矢量数据，转换成对应的栅格数据。也就是矢量向栅格数据转换处理的根本任务。l栅格化可以分别针对点、线和面进行，其中，点的栅格化是线和面的栅格化的基础。l要将矢量数据的平面直角坐标转换成栅格数据以行列表示的栅格坐标，首先要要建立矢量数据的平面直角坐标系和栅格行列坐标系直角的对于关系。2024/8/234444栅格化技术方法栅格化技术方法栅格化技术点的栅格化线的栅格化面的栅格化扫描线算法基于弧段的栅格化基于多边形的栅格化扫描线算法内点填充法边界代数法包含检验法检验夹角之和法铅垂线法（交点
26、个数法）2024/8/234545（一）点的栅格化（一）点的栅格化l设矢量坐标点（x，y），转换后的栅格单元行列值为（I，J），则有：式中：方括号表示取整数运算。2024/8/234646（二）线的栅格化（二）线的栅格化l线的栅格化可以分解成对组成折线的每一个线段的矢量化。l对于线的栅格化：先使用点栅格化的方法，栅格化线段的两个端点，然后再栅格化线段中间的部分。l对于线段中间的部分的栅格化，需要分两种情况来处理：设线段两端点坐标分别为（x1,y1）和（x2,y2），栅格化后的单元行列值分别为（I1，J1）和（I2，J2），则行数差I =|I2-I1|，列数差： J =|J2-J1|一是列数差大
27、于行数差JI的情况二是行数差大于列数差I J ，运用扫描线算法实现。47列数差大于行数差列数差大于行数差l当列数差大于行数差的时候，平行于y轴做每一列的中心线，称为扫描线，求每一条扫描线与线段的交点，按点的栅格化方法将交点转换为栅格坐标。l设：xm为每列的中心扫描线的横坐标， （x1,y1）和（x2,y2）为线段的两个端点坐标，则交点坐标为：X=xm2024/8/234848行数差大于列数差行数差大于列数差l当行数差大于列数差的时候，平行于x轴做每一行的中心扫描线。再求每一条扫描线与线段的交点，按点栅格化方法将交点转换为栅格坐标。l设：ym为每行的中心扫描线的纵坐标，则交点坐标为：Y=ym20
28、24/8/234949（三）面的栅格化（三）面的栅格化l根据多边形矢量数据的不同结构，栅格化有两种不同的方法，即： 基于弧段数据的栅格化拓扑数据结构的矢量数据基于多边形数据的栅格化非拓扑数据结构的矢量 数据。2024/8/2350501、基于弧段数据的栅格化方法、基于弧段数据的栅格化方法l是针对拓扑结构的矢量数据进行栅格化的方法。其方法与线段栅格化的方法类似。l栅格化步骤：对整个要进行栅格化的范围按行或者按列作中心扫描线，对其中的任一条扫描线，求与所有矢量多边形的边界弧段的交点，记录其坐标，并用点栅格化方法求出交点的栅格坐标行列值，再根据弧段的左右多边形信息判断并记录交点多边形的数值。最后通过
29、对一行所有交点按其坐标x值从小到大进行排序，并参照左右多边形配对情况，逐段生成栅格数据。直到全部扫描线都完成从矢量向栅格的转换为止。2024/8/2351512、基于多边形数据的栅格化方法、基于多边形数据的栅格化方法l是针对实体结构的多边形矢量数据栅格化l原理：以非拓扑的实体多边形作为栅格化的处理单元，将一个多边形的内部栅格对于赋予多边形的属性值。l基于多边形数据的栅格化方法包括：内点填充法边界代数法包含检验法2024/8/235252（1）内点填充法： 首先按线的栅格方法把多边形的边界栅格化，然后在多边形的内部找一点作为内点，从该点出发，向外填充多边形区域，直到边界位置。（2）边界代数法：沿
30、着多边形实体的边界环绕多边形一圈，当向上环绕的时候，把边界左边一行中所有的栅格单元的数值都减去属性值，当向下环绕的时候，把边界左边一行中所有的栅格单元的数值都加上属性值，则多边形外部的栅格正负数值抵消，而多边形内部的栅格被赋予属性值。2024/8/235353（3）包含检验法：对每一个栅格单元，逐个判定其是否包含在某个实体多边形之内，若包含在某个多边形之内，则将多边形的属性值赋给该栅格单元。包含检验法实现的关键在于“点在多边形内”的判定，或称为point-in-polygon分析l“点在多边形内”的判定可以通过检查夹角之和的方法或铅垂线法（检验交点数的方法）来实现。2024/8/235454l
31、1.检验夹角之和：对于任一待判别的栅格点P，到一个多边形实体各个顶点连线，构成的夹角之和为0，则点P在多边形之外；如果夹角之和为2，则P在多边形之内l2.铅垂线法（检验交点数法）由任一待判别的栅格点P向下作与y轴平行的射线（称为铅垂线），计算该射线与多边形边界的交点数。若交点数为0或偶数，则栅格点P在多边形之外；若交点数为奇数，则栅格点P在多边形之内。2024/8/235555二、由栅格向矢量的转换二、由栅格向矢量的转换l矢量化：栅格数据结构向矢量数据结构的转换。l矢量化的目的：一是将扫描仪获取的图像栅格数据存入矢量形式的空间数据库；二是将栅格数据进行数据压缩，将面状数据转换为由矢量数据表示的
32、多边形边界。56l栅格向矢量的转换过程比较复杂，它有两种情况：一种是本身为遥感影像或己栅格化的分类图另一种情况通常是从原来的线划图扫描得到的栅格图2024/8/235757（一）基于图像数据的矢量化方法（一）基于图像数据的矢量化方法l图像数据线划图扫描得到的栅格图是不同灰阶的影像，通过扫描仪按一定的分辨率进行扫描采样，得到不同灰度值（0255）表示的数据。l矢量化步骤：由于扫描仪的分辨率一般可达0.0125mm，因此，对一般粗细（如0.1mm）的线条，其横断面扫描后也会有大约8个像元，而矢量的要求只能允许横断面保持一个栅格的宽度，因此，需要进行二值化、细化和跟踪等矢量化步骤。2024/8/23
33、5858矢量化步骤矢量化步骤式中：f(i,j)为原像元灰度。 二值图像中的图形用1表示，背景用0表示。2024/8/235959l（2）细化：细化是消除线划横断面栅格数的差异，使得每一条线只保留代表其轴线或周围轮廓线（对面状符号而言）位置的单个栅格的宽度。细化也称为栅格数据的轴化，就是将占有多个栅格宽的图形要素缩减为只有单栅格宽的图形要素的过程。 l对于栅格线画的细化方法，可以分为： “剥皮法”和“骨架法”2024/8/236060“剥皮法剥皮法”l剥皮的实质就是从曲线的边缘开始，每次删掉外层等于一个栅格宽的一层，直到最后留下彼此连通的由单个栅格组成的图形。l在剥皮过程中必须注意不允许剥去会导
34、致曲线不连通的栅格。解决办法：用一个3*3的栅格窗口，在栅格图上逐个检查每个栅格单元。被查栅格能否删去，由以该栅格为中心的组合图来决定，其原则是不允许剥去会导致图形不连通的栅格，也不能在图形中形成孔。2024/8/236161骨架法骨架法l这种方法就是确定图形的骨架，而将非骨架上的多余栅格删除。l具体做法是扫描全图，凡是像元值为1的栅格都用V值取代。V值是该栅格与北、东和北东三个相邻栅格像元值之和，即l在V值图上保留最大V值的栅格，删去其他栅格，但必须保证连通。因为最大V值的栅格只能分布在图形的中心线上（骨架上），因此选取最大值栅格的过程就是细化的过程2024/8/2362622024/8/2
35、36363l（3）跟踪：跟踪的目的是将细化处理后的栅格数据转换为从节点出发的线段或闭合的线段，并以矢量形式存储线段的坐标。跟踪时，从起点开始，根据八个领域进行搜索下一个相邻点的位置，记录坐标，直到完成全部栅格数据的矢量化。2024/8/236464（二）栅格数据的矢量化方法（二）栅格数据的矢量化方法l针对遥感影像或己栅格化的分类图l栅格数据的矢量化常常针对栅格数据中的多边形进行，其步骤：首先，在栅格数据中搜索多边形边界弧段相交处的节点位置，这些节点通常是相邻栅格单元不相同的属性值个数大于等于3的栅格处。接着，从搜索出的节点里任选一个作为起始跟踪点，顺着栅格单元属性值不同的两个栅格单元之间进行多
36、边形边界弧段的跟踪，记录每一步跟踪的坐标，直到另一个节点为止，则完成一条边界弧段的跟踪矢量化。重复上述过程，做到所有的边界弧段都被生成。最后，将跟踪得到的弧段数据连接组织成多边形，则完成多边形栅格数据的矢量化。2024/8/236565边界提取边界提取l边界提取是遥感图像处理中的一个专门问题。l用一个2x2栅格的窗口，按顺序沿行列方向对栅格图像进行扫描。如果窗口内的四个网格点值相同，它们就属于一个等值区，而无边界通过，否则就存在多边形的边界或边界的结点。如果窗口内有两种栅格值，这四个栅格则均标识为边界点，同时保留原栅格的值。如果窗口内有三个以上不同的值，则标识为结点。2024/8/236666
37、边界点的六种可能结构边界点的六种可能结构 结点的八种可能结构结点的八种可能结构 2024/8/236767第三节第三节 多元空间数据的融合多元空间数据的融合l多元数据的产生的背景：GIS技术经过40多年的发展和应用，已经积累了大量的数据资源。但是由于地理数据的多语意性、多时空性、多尺度性、获取手段的多样性、存储格式的不同以及数据模型与数据结构的差异等，导致多元数据的产生。l多元数据产生的影响：多元数据的产生给数据的集成和信息共享带来了困难。l多源数据融合的目的：为了实现空间数据的共享，特别是随着因特网的发展、数字地球的兴起和GIS应用的日益深入，多源数据的融合已成为GIS设计者和用户的共同要求
38、。2024/8/236868一、遥感与一、遥感与GIS数据的融合数据的融合l遥感数据与GIS数据借助遥感技术获得的信息具有周期动态性、信息丰富、获取效率高等优势；GIS则具有高效的空间数据管理和灵活的空间数据综合分析能力。l遥感数据与GIS数据融合的优点：有利于增强多重数据的复合能力，改善遥感信息提取的及时性和可靠性，便于利用遥感影像辅助GIS空间数据的获取与更新，有效地提供各类数据的使用率。2024/8/236969遥感与遥感与GIS数据的融合，目前最常用的方法数据的融合，目前最常用的方法l（1）遥感影像与数字线画图（DLG）的融合：成果：经过正射纠正后的遥感影像，与数字线画图信息融合，可产
39、生影像地图。优点：这种影像地图具有一定的数学基础，有丰富的光谱信息与几何信息，又有行政界线和属性信息，直接提高了用户的可视化效果。l（2）遥感影像与数字地形模型（DEM）的融合：数字地形模型与遥感数据的融合，有助于实施遥感影像的几何校正与配准，消除遥感影像中因地形起伏所造成的像元位移，提高遥感影像的定位精度，同时，数字地形可参与遥感影像的分类、改善分类精度；2024/8/237070l（3）遥感影像与数字栅格图（DRG）的融合：将数字栅格地图与遥感图像配准叠合，可以从遥感图像中快速发现已发生变化的区域，进而实现空间数据库的自动/半自动更新。2024/8/237171二、不同格式数据的融合二、不
40、同格式数据的融合l必要性：由于GIS软件的多样性，每种GIS软件都有自己特定的数据模型，造成数据存储格式和结构的不同。l常见GIS软件的空间数据格式：ESRI公司的ArcInfo Coverage、Shapefile、E00格式；Autodesk公司的DXF和DWG格式；MapInfo公司的TAB和MIF格式；Intergraph公司的DGN格式；（其中红色加粗格式是用于数据交换格式）2024/8/237272解决不同格式数据之间的融合方法解决不同格式数据之间的融合方法l（1）基于转换器的数据融合：在这种模式下，数据转换一般通过交换格式进行。如，要转换Mapinfo的TAB文件到ArcInfo
41、的Coverage：首先需要使用Mapinfo软件将TAB文件输出为E00或DXF文件，然后运行ArcInfo将E00或DXF文件转换为Arcinfo Coverage。这是目前GIS系统数据融合的主要方法。存在的主要问题： 数据转换过程复杂，系统内部的数据格式需要公开等。2024/8/237373l（2）基于数据标准的数据融合：采用一种空间数据的转换标准来实现多源GIS数据的融合。如：美国国家空间数据协会（NSDI）制定了统一的空间数据格式规范SDTS（Spatial Data Transformation Standard），包括几何坐标、投影、拓扑关系、属性数据、数据字典等，也包括栅格和
42、矢量等不同空间数据格式的转换标准。根据SDTS，目前很多GIS软件提供了标准的空间数据交换格式，如Arcinfo的SDTSIMPORT和SDTSEXPORT模块，可供其他系统调用。优点：这种转换方法能够处理多个数据集、转换次数少，系统内部的数据格式不需公开，只要公开转换采用的技术即可。2024/8/237474l（3）基于公共接口的数据融合：又称为数据互操作模式。接口相当于一种规程，在接口中不仅要考虑数据格式和数据处理，而且还要提供对数据处理应采用的协议，各个系统通过公共接口相互联系，而且允许鸽子系统内部数据结构和数据处理不相同。如OGC(Open GIS Consortium)为数据互操作制
43、定了统一的规范，从而使一个系统同时支持不同的空间数据格式成为可能。主要特点：独立于具体平台，转换技术高度抽象，数据格式不需公开，代表着数据共享技术的发展方向。2024/8/237575l（4）基于直接访问的数据融合：直接数据访问指的是一个GIS软件中实现对其他软件数据格式的直接访问，用户可以使用单个GIS软件存取多种数据格式。直接数据访问不仅避免了繁琐的数据转换，而且在一个GIS软件中访问其他软件的数据格式，不要求用户拥有该数据格式的宿主软件，更不需要该软件的运行，这为多源数据的共享提供了更为实用便捷的支持。如Intergraph公司推出的GeoMedia系列软件提供了这种支持。2024/8/
44、237676第四节第四节 空间数据的压缩与重分类空间数据的压缩与重分类 一、空间数据的压缩1、数据压缩的意义：l优化存储空间，减少处理时间优化存储空间，减少处理时间2、空间数据压缩：从空间坐标数据集合中抽取一个子集，使这个子集在规定的精度范围内最好地逼近原集合，而又取得尽可能大的压缩比。2024/8/237777压缩比表示信息载体减少的程度。 a= m/n1 a= m/n1式中：式中：式中：式中：m m m m为曲线的原点数；为曲线的原点数；为曲线的原点数；为曲线的原点数；n n n n为曲线经压缩后的点数。为曲线经压缩后的点数。为曲线经压缩后的点数。为曲线经压缩后的点数。2024/8/237
45、878（二）基于矢量的压缩（二）基于矢量的压缩l基于矢量的压缩通常是对线状实体上点的数量的压缩。l最常用的是道格拉斯-佩克算法（也被翻译为道格拉斯-普克算法）该算法是D.H.Douglas和T.K.Peucker在1973年提出的。它是基于线状实体的点压缩算法，是对数据的简化处理。2024/8/237979l算法的原理：先拟定一个阈值，然后生成一条连接折线首尾节点的直线段，并计算原始折线上的点到直线段的距离。假如所有折线上的点到直线段的距离都小于预先设定的阈值，这条直线段就被用来代替原来的那条折线；假如有些点的距离大于阈值，距离最远的那一点保留，并将原折线分成两段。对两段折线重复上述过程，最后
46、保留下来的点就是经过数据压缩的折线。 2024/8/238080（三）基于栅格的压缩（三）基于栅格的压缩l栅格数据的压缩可以采用游程编码和四叉树等编码方法。2024/8/238181二、空间数据的重分类二、空间数据的重分类l存储在空间数据库中的数据，是提供为多种目标服务的。当需要进行特定的数据分析时，常常需要先对从数据库中提取的数据做属性的重新分类和空间图像的化简，以构成数据新的使用形式。BBAACBBACB消除具有相同属性相邻多边形的公共边界2024/8/238282第五节第五节 空间数据的内插方法空间数据的内插方法l空间数据的内插：通过已知点或多边形分区的数据，推求出任意点或多边形分区数据
47、的方法称为空间数据的内插。l在已观测点的区域外估算未观测点的数据的过程称为外推称为外推2024/8/238383l根据已知点和已知多边形分区数据的不同，将空间数据内插分为分类：点的内插；多边形分区的内插。l应用：广泛应用于生成等值线、建立数字高程模型、不同区域范围现象的相关分析和比较研究等。2024/8/238484一、点的内插一、点的内插l点的内插：是用来建立具有连续变化特征现象（如地面高程、地形、气温、气压等）的数值方法。l理论基础：空间相关性l空间相关性：对于地理上连续分布的现象，邻近点之间关联性强，较远的点之间关联性弱或者无关。2024/8/238585内插方法内插方法l建立数字高程模
48、型的点的内插方法数字高程模型的建立一般包括包括数据取样、数据内插和数据精度分析数据取样、数据内插和数据精度分析等步骤。2024/8/238686(一一)数据取样数据取样l建立数字高程模型，就要生成按网格形式排列的地面点高程。可以先从现有地形图等高线上进行数据取样。l使用随机取样方案确定取样点：取样点：可以沿着地性线（山脊线、山谷线、坡度变换线），或沿着等高线，或沿着断面线即：数据点应选择在地性线的坡度改变处，或沿等高线在方向改变的地点。优点：通过这种方式进行取样，数据点在地形特征点上，能很好的控制地表面形状。2024/8/238787l高程点格网密度的选择当按照一定间隔的格网取样时，取一个二次
49、曲线来代表相应部分的曲线，而用间隔中点的线性内插值与二次曲线内插值所得的高程的差异来判断格网密度是否恰当。如图3-35l高程差异的数值要在数字高程模型精度要求的限度之内，否则要缩短格网间距。l对地形变化显著的地方，要辅以地形特征点位数据点，以提高数据点控制地面形态特征的能力。2024/8/238888（二）数据内插（二）数据内插l插值原因：由于取样的数据呈离散点分布形式，或者数据虽然按照格网排列，但格网的密度不能满足使用的要求，这样就需要以数据点为基础进行插值运算。l插值原理：插值运算就是选择一个合理的数学模型，利用已知点的数据求出插值函数的待定系数。l插值方法：由于地面形态千变万化，既无规律
50、又无重复性，通常采用局部分块内插法和逐点内插法。2024/8/2389891、分块内插法、分块内插法l原理：把整个内插空间化成若干分块，并对各分块求出各自的曲面函数来刻画曲面形态。l分块内插的关键：是要解决各相邻分块函数间的连续性问题。l分块内插分为：线性内插法、双线性多项式内插法和二元样条函数内插法。2024/8/239090（1）线性内插法）线性内插法l原理：先将所有的已知数据点连接成三角网的形式，使用靠近内插点的三个已知数据点，来确定三角网中的一个三角形形成的空间平面，继而求出该内插点在平面中的高程值。2024/8/239191线性内插函数线性内插函数l所求的线性内插函数形式为l待定系数
51、将内插点周围的3个数据点的数据代入多项式，即可解算出系数a0、a1、a2 。进而带入内插点平面坐标x，y之后，即可求出内插高程Z。2024/8/239292（2 2）双线性多项式内插法）双线性多项式内插法l原理：在规则分布的已知数据点时，使用最靠近内插点的四个已知数据点组成一个四边形，确定一个双线性多项式来内插其中点的高程2024/8/239393双线性多项式函数双线性多项式函数l设确定的函数形式为：l将内插点周围的4个数据点的数据值代入多项式，即可解算出系数a0、a1、a2 、a3 。2024/8/239494（3 3）二次样条函数内插法（双三次多项式）二次样条函数内插法（双三次多项式）l在
52、分块插值区用双三次多项式即样条函数模拟地表面。l二次样条函数内插法：是一种分段函数,每次只用少量的数据点；样条函数通过所有的数据点，故可用于精确的内插；可用于平滑处理。2024/8/239595双三次多项式内插双三次多项式内插l双三次多项式内插的多项式函数为（待定点高程）：l将内插点周围的16个点的数据带入多项式，可计算出所有的系数。2024/8/2396962、逐点内插法、逐点内插法l分块内插法：分块内插法：分块范围在内插过程中一经确定，其形状、大小和位置都保持不变。凡落在分块上的待插值点都用展铺在该分块上的唯一确定的数学面进行内插。l逐点内插法逐点内插法则是以插值点为中心，定义一个局部函数
53、去拟合周围的数据点，数据点的范围随插值点位置的变化而变化，因此又称移动曲面法。移动曲面法。2024/8/239797逐点内插法分类逐点内插法分类l逐点内插法主要有两种基本的插值方法：移动拟合法加权插值方法l注意：克里金法也是一种加权插值方法，只是在计算权重的方法上与加权平均法不同。2024/8/239898（1）移动拟合法：）移动拟合法：l该方法是指对每一个待插值点P，用一个多项式曲面拟合该点附近的表面，从而计算出该点的高程值。此时，取待插值点P为圆心，R为半径的圆（称为搜索圆）内各数据点来计算多项式的待定系数。l取二次多项式来拟合曲面，则待插值点的高程可以写成：z=Ax2+Bxy+Cy2+D
54、x+Ey+F式中：A、B、C、D、E、F为待定系数。l这些待定系数可由落在搜索圆中的n个数据点用最小二乘法进行求解，2024/8/239999搜索圆的半径搜索圆的半径R的确定方法的确定方法l运用移动拟合法的关键：就在于如何确定待插值点的最小邻域范围（搜索圆的半径R），以保证邻近数据点的数量足够计算6个待定系数。l搜索圆的半径R的确定方法，需要考虑两个方面的问题：一是考虑范围，即采用多大面积范围内的数据点来计算待插值点的数值；二是考虑点数，即选择多少数据点参加计算，由此确定需要多大的范围。2024/8/23100100动态搜圆动态搜圆l采用动态搜圆的方法其思路是从数据点的平均密度出发，确定圆内数
55、据点n（如平均n要有12个），以计算搜索圆的半径R，其公式为：R2=n(A/N)=R=l式中N为数据点总数，A为总面积。l特点：1、这种方法实际上综合考虑了点数和范围两个因素。2、先以此半径作搜索圆进行搜索，如果落在搜索圆内的点数大于n，则符合计算要求，可以进行曲面的拟合与插值计算。否则，动态扩大搜索圆的半径，直到落在圆内的数据点个数符合要求为止。2024/8/23101101（2）加权平均法）加权平均法l移动拟合法移动拟合法，往往需要求解复杂的误差方程组来求取曲面函数的待定系数。在实际应用中，更为常用的是加权平均法，可以看作是移动拟合法的特例。l加权平均法加权平均法在使用搜索圆寻找附近数据点
56、的方法上和移动拟合法相同，但加权平均法在计算待插值点的高程时，使用加权平均值代替误差方程求解出的曲面函数2024/8/23102102使用加权平均值求曲面函数使用加权平均值求曲面函数l式中：n为落在搜索圆中的数据点的个数； zi为落在搜索圆中的第i个数据点的高程值；pi为第i个数据点的权重。l权重的计算由于考虑到不同的数据点相对于待插值点的距离不同，对待插值点的高程影响程度不同，所以，一般采用与距离相关的权函数来计算权重，如下式：2024/8/23103103第第i个数据点的权重个数据点的权重l式中：p是数据点的权重；R是搜索圆的半径；d是待插值点到数据点的距离。l这种加权平均法又称为反比距离
57、加权方法反比距离加权方法（inverse distance weighted，IDW）2024/8/23104104(3)克里金法（克里金法（Kriging）l前面介绍的一些空间插值方法对影响插值相关的一些问题没有得到很好的解决，如反比距离加权法中很多人为的因素影响最终的插值效果，即需要人为的设定搜索圆中的数据点的数目、搜索圆的大小、计算权重的方法等，此外，插值的精度或误差难以确定。l为解决这些问题，法国地理数学家G.Matheron和南非矿山工程师D.G.krige研究了一种优化插值方法。Krige首先将这一技术用于更加精确的推测金属矿物的储量。2024/8/23105105Kriging插
58、值方法插值方法lKriging方法是基于一个假设：即被插值的某要素（如地形要素）可以被当做一个区域化的变量来看待。l所谓区域化的变量就是介于完全随机的变量和完全确定的变量之间的一种变量，它随所在区域位置的改变而连续地变化。因此，彼此离得近的点之间有某种程度上的空间相关性，而相隔比较远的点之间在统计上看是相互独立无关的。l通过几十年的实践运用，克里金方法已经成为地统计学（Geostatistics）的基础工具。2024/8/23106106Kriging插值方法插值方法lKriging方法是建立在一个预先定义的协方差模型的基础上，通过线性回归方法把估计值的方差最小化的一种插值方法。lKrigin
59、g方法具体分成许多种，主要有：普通Kriging、简单Kriging和通用Kriging等。2024/8/23107107普通普通KrigingKrigingl首先是利用那些将要用来插值的离散点集合建立一个变量图，变量图通常包括两部分：变量图通常包括两部分：一个是根据实验获得的变量图，另一个是模型变量图。假设要插值的数值用z表示，则通过计算集合中的每一个点相对于其他点的差异，并且用差异和对应的点之间的距离作图，就可以得到根据实验获得的变量图。l通常用来计算的方法是求z差值平方的一半。这样的变量图又称为半方差图。图3-402024/8/23108108l一旦实验获得的变量图计算完成后，就是定义一
60、个模型变量图。l模型变量图是一个简单的数学函数,用来模拟实验获得的变量图的趋势。l如图3-40所示，变量图的形状表明，在小的间距上z的差异是很小的，也就是说，彼此靠近的点有近似z值。在某个距离的间隔之后，z值的差异某种程度上变得随机而没有规律，因此，半方差图在某一个值上变得平直下去，该值是差异的平均值。2024/8/23109109l一旦模型变量图建立后，就被用来计算Kriging方法中的权重，在普通Kriging方法中运用的基本公式如下：l式中：n是集合中离散点的个数； zi是离散点的数值； wi是赋予每个离散点的权重。Kriging方法的这个公式和反比权重插值的公式基本相同，只是权重不是基
61、于一个任意的距离函数，而是基于模型变量图。例如，利用P点周围的三个点P1、P2、P3在P点插值，必须先找到w1、w2、w3，权重可以通过求解联立方程获得（3-32）2024/8/23110110l用变量图来计算权重，则预期的估计误差被用最小二乘方式最小化了。因此，Kriging方法被认为能产生最优的线性无偏估计。lKriging方法的一个重要特点是变量图可以被用来对每一个插值点计算估计的预期误差，因为估计的误差是到周围离散点距离的一个函数。2024/8/23111111简单简单Kriging方法方法l与普通Kriging方法相似，区别在于没有把方程w1+w2+w3=1加入方程组，并且权重相加也
62、不等于一。l简单Kriging是用整个数据集合来平均，而普通Kriging是采用局部平均（对一个插值点的离散点子集的平均）。l简单Kriging不如普通Kriging精确，但通常会产生一个更加平滑的结果。2024/8/23112112通用通用Krigingl在Kriging方法里有一个假设，被估计的数据是固定不变的。这就是说在离散点集合里从一个区域移动到下一个区域，离散点的平均值是相对恒定的。只要在数据值里面存在一个显著的空间趋势。例如一个倾斜的表面或一个局部的平坦区域，这一假设就不成立了。在这种情况下，可以通过运用一个“漂移”项来临时地使用固定条件。l这个漂移是一个简单的多项式函数用来模拟离
63、散点的平均值。l残差是漂移和离散点实际值之间的差，因为残差应该是固定的，Kriging方法带着残差执行，并且插值的残差被加到漂移上来计算估计值。使用一个偏移的这种方式通常称为“通用Kriging”方法。2024/8/23113113（三）数据精度分析（三）数据精度分析l为了有效地判定空间插值结果的精度是否符合实际的要求，必须对建立的数字高程模型进行精度分析。l精度分析的常用方法是，在研究区内随即选n个数据点，设这些点的图上高程值为z i，其对应的插值数值为zi，他们之差用zi=z i-zi表示，则： 算术平均值为2024/8/23114114l标准差为：l根据高程数据精度的评价标准，其允许误差
64、（z）应小于原始数据比例尺等高举的1/3.若符合该条件，即认为所建立的数字高程模型满足精度要求，可以作为数字地形模型分析的依据。2024/8/23115115二、区域的内插二、区域的内插l区域的内插是研究根据一组多边形分区的已知数据来推求同一地区另一组多边形分区未知数据的内插方法；l如图3-41，设一组已知数据的分区称为源区，需要内插的另一组分区称目标区，则根据源区的数据来推求目标区的数据，可以采用以下两种方法:叠置法叠置法,比重法比重法。 ABC 213(a)源(b)目标2024/8/23116116（一）叠置法（一）叠置法l将目标区叠置在源区上，首先确定两者面积的交集ats然后利用下式计算
65、出目标区各个分区t的内插值vt式中：t为目标区各个分区的序号； s为源区各个分区的序号； U为分区s的已知统计数据； ats为t区与s区相交的面积； s为s区的面积2024/8/23117117（二）比重法（二）比重法l是根据平滑密度函数的原理，将源区的统计数据从同性质改变为非同性质，而非同性质代表着一般社会经济现象的普通特点。l如图3-42说明比重法的内插算法过程。2024/8/231181185.05.05.05.05.05.05.05.05.05.05.05.05.05.05.05.05.05.05.05.05.05.05.05.0 5.05.03.33.33.33.33.33.35.0
66、5.05.05.04.64.64.24.24.44.45.05.05.05.05.05.04.24.25.05.05.05.05.05.0 4.24.24.64.63.93.93.33.35.35.35.35.34.84.84.44.44.64.65.35.35.05.05.35.34.54.55.35.35.35.35.35.3 4.54.53.63.63.03.03.43.45.55.55.15.15.15.14.04.04.64.65.35.35.55.55.25.24.54.55.55.55.55.54.94.9 4.54.53.53.53.73.73.33.32024/8/23119
67、119克里格（克里格（Kriging）内插法举例）内插法举例高程数据高程数据2024/8/231201202024/8/23121121插值结果插值结果2024/8/23122122三维显示三维显示2024/8/23123123第第6节节 空间拓扑关系的编辑空间拓扑关系的编辑l拓扑关系拓扑是点、线和多边形要素共享几何的方式的排列布置。l拓扑的主要用途：限制要素共享几何的方式。例如，相邻多边形（如宗地）具有共享边、街道中心线和人口普查区块共享几何以及相邻的土壤多边形共享边。定义并执行数据完整性规则：多边形之间不应存在任何间距、不应有任何叠置要素等。支持拓扑关系查询和导航，如确定要素邻接性和连通性
68、。支持可强制执行数据模型拓扑约束的复杂编辑工具。根据非结构化的几何构造要素，如根据线创建多边形。2024/8/23124124拓扑编辑拓扑编辑l拓扑数据结构最重要的特征是具有拓扑编辑功能，l拓扑编辑功能包括多边形连接编辑多边形连接编辑:指顺序连接组成封闭多边形一组线段的编辑节点连接编辑节点连接编辑:指顺序连接环绕某个节点所有多边形的编辑。2024/8/23125125一、多边形连接编辑一、多边形连接编辑l设需要对多边形P1进行编辑，其算法过程为：1、从弧段文件中，找出与当前编辑的多边形P1相关的所有记录；2、在找出的记录中，检查当前编辑的多边形P1所处的位置；如果P1处在左多边形位置，将之与处
69、于右多边形位置的多边形号相交换，同时也将该记录的节点号位置作相应的交换；反之，如果当前编辑的多边形P1处于右多边形位置，则该记录的所有数据项顺序不作改变。3、从经过代码位置转换的记录中，任取一个起始节点作为起点，顺序连接各个节点，必要时可对记录的前后顺序作调整，使得连接的节点能自行封闭，即N1-N2-N3-N1。如果依照上述顺序连接的节点不能自行闭合，或者出现记录缺损或记录多余等情况，则表示弧段文件有错，必须改正出错的记录。直到所有多边形都经过编辑和改正，再转入节点连接编辑。2024/8/23126126二、节点连接编辑二、节点连接编辑设需要对节点N2进行编辑，其算法过程为1、从弧段文件中，找
70、出与当前编辑的节点N2相关的所有记录；2、在找出的记录中，检查当前编辑的节点N2所在的位置；如果N2处在起始节点位置，将之处于终止节点位置的节点号相交换，同时也将该记录的多边形号位置作相应的交换；反之，如果当前编辑的节点N2处于终止节点位置，则该记录的所有数据项顺序不作改变。3、从经过代码位置转换的记录中，任取一个左多边形作为起点，顺序连接各个多边形，同样必要时可对记录的前后顺序作调整，使得连接的多边形能首尾呼应。即P2-P1-P4-P2如果依照上述顺序连接的多边形不能首尾呼应，或者出现记录缺损或记录多余等情况，同样也表示弧段文件有错，不想改正出错的记录。直到所有节点都经过编辑和改正，才能将该弧段文件应用于数据库的建立。2024/8/231271272024/8/23128128
