空间数据结构的转换分析课件

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,计算机图形学演示稿 纪玉波制作(C),*,第2节 空间数据结构的转换,矢量数据到栅格数据的转换，称为矢量栅格化。,许多数据如行政边界、交通干线、土地利用类型、土壤类型等都是用矢量数字化的方法输人计算机或以矢量的方式存在计算机中，表现为点、线、多边形数据。然而，矢量数据直接用于多种数据的复合分析等处理将比较复杂，特别是不同数据要在位置上一一配准。,11/14/2024,1,第2节 空间数据结构的转换 矢量数据到栅格数据的转,具体转换步骤分如下几步,：,1 准备好矢量数据或矢量地图；,2 构架一个与地图等大小的格网，选取好适当的格网密度或分辨率，依据范围则可以算出格网的行列数。,3 根据地图上的点、线、多边形相对这个格网的配置及其属性来确定每一个项元的数值。,11/14/2024,2,具体转换步骤分如下几步：9/21/20232,上述三步中，第三步内容较多，需要专门说明：,1、点状对象的栅格化：,点的变换十分简单，只要这个点落在哪个网格中就是属于那个网格元素，根据该点状对象的特性赋予该像元属性值。,点行、列坐标i，j计算公式：,11/14/2024,3,上述三步中，第三步内容较多，需要专门说明：9/21/,2、线状对象的栅格化（一）：,假定一线段两端点之间经过若干个网格元素(至少一个)，两端点坐标为(X1，Y1)，(X2，Y2)，则：P82,Y,X,（X1，Y1）,（X2，Y2）,11/14/2024,4,2、线状对象的栅格化（一）：YX（X1，Y1）（X2，Y,判断行差和列差是为了确保单一连接：,单一连接定义：,1、线条不中断，线上每两个像元间必须连接，或邻边连接，或顶角连接。,2、任何三个像元不得两两间都相邻。,11/14/2024,5,判断行差和列差是为了确保单一连接：9/21/2023,3、面的栅格化方法：,一、基于弧段数据的栅格化方法：,按行或列做中心扫描线,求交点栅格坐标的行列值,判断交点左右多边形数值,交点排序并逐段生成栅格数据,特点：,计算量较大，算法比较复杂。,11/14/2024,6,3、面的栅格化方法：一、基于弧段数据的栅格化方法：9/21/,1,1,0,1,1,1,1,0,1,1,1,0,0,0,0,1,1,1,1,0,1,0,0,1,11/14/2024,7,1101111011100001111010019/21/2,11/14/2024,8,9/21/20238,二、基于多边形的栅格化方法：,、内部点扩散法：,由一个内部的种子点，向其四个方向的邻点扩散，判断新加入的点是否在多边形边界上。如果是，不作为种子点；否则当作新的种子点，直到区域填满。,算法特点：,算法设计复杂，而且可能造成阻塞而使扩散不能完成。,11/14/2024,9,二、基于多边形的栅格化方法：9/21/20239,、边界代数法：,上减下加,特点：,算法简单，但对于复杂图形，每一像素可能被访问多次，增加了运算量。,11/14/2024,10,、边界代数法：9/21/202310,、包含检验法：,点在多边形内判断：,检验夹角之和,射线算法：,由待判点向图外某点引射线，判断该射线与某多边形所有边界相交的总次数，如果相交偶数次，则待判点在该多边形外部，如为奇数次，则待判点在该多边形内部如图所示。,11/14/2024,11,、包含检验法：9/21/202311,三、栅格到矢量：,栅格向矢量转换处理的目的，是为了：,1 将栅格数据分析的结果，通过矢量绘图装置输出；,2 数据压缩的需要，将大量的面状栅格数据转换为由少量数据表示的多边形边界；,3 将自动扫描仪获取的栅格数据加入矢量形式的数据库。,11/14/2024,12,三、栅格到矢量：9/21/202312,1、基于图象数据的矢量化方法,图象数据是由不同灰阶的影像或线划，通过自动扫描仪(scanner)，按一定的分辨率进行扫描采样，得到以不同灰度值(0255)表示的数据。目前扫描仪的分辨率可达00125mm，因此对一般粗度(例如0.1mm)的线条，其横断面扫描后平均也有8个像元，而矢量化的要求只能允许横断面保持一个栅格的宽度，因此需要进行从栅格向矢量数据的转换。,11/14/2024,13,1、基于图象数据的矢量化方法9/21/202313,（a）扫描前的矢量数据 （b）扫描得到的灰度值,11/14/2024,14,（a）扫描前的矢量数据 （b）,具体转换的步骤：,I、二值化 线划图形扫描后产生栅格数据，这些数据是按从0255的不同灰度值量度的，设以G(i，j)表示，为了将这种256级不同的灰阶压缩到2个灰阶，即0和1两级，首先要在最大与最小灰阶之间定义一个阈值，设阈值为T，则如果G(i，j)大于等于T，则记此栅格的值为1，如果G(i，j)小于T，则记此栅格的值为0，得到一幅二值图。,11/14/2024,15,具体转换的步骤：9/21/202315,根据给定的阈值二值化后得到的栅格数据,11/14/2024,16,根据给定的阈值二值化后得到的栅格数据9/21/202316,II、细化 细化是消除线划横断面栅格数的差异，使得每一条线只保留代表其轴线或周围轮廓线(对面状符号而言)位置的单个栅格的宽度。,11/14/2024,17,II、细化 细化是消除线划横断面栅格数的差异，使得每一条线,对于栅格线划的“细化”方法，常用 “剥皮法”。剥皮法的实质是从曲线的边缘开始，每次剥掉等于一个栅格宽的一层，直到最后留下彼此连通的由单个栅格点组成的图形。因为一条线在不同位置可能有不同的宽度，故在剥皮过程中必须注意一个条件，即不允许剥去会导致曲线不连通的栅格。,11/14/2024,18,对于栅格线划的“细化”方法，常用 “剥皮法”。剥皮法的实,这是这一方法的技术关键所在。其解决办法是，借助一个在计算机中存储着的，由待剥栅格为中心的33栅格组合图来决定。一个33的栅格窗口，其中心栅格有八个邻域，因此组合图有多种不同的排列格式，若将相对位置关系的差异只是转置90、180、270度或互为镜象反射的方法进行归并，则共有51种排列格式。,11/14/2024,19,这是这一方法的技术关键所在。其解决办法是，借助一个在计算机中,显然，其中只有格式2、3、4、5、10、11、12、16、21、24、28、33、34、35、38、42、43、46和50，可以将中心点剥去。这样，便可确定中间栅格点保留或删除，得到经细化处理后应予保留的栅格系列，并写入数据文件。,11/14/2024,20,显然，其中只有格式2、3、4、5、10、11,III、跟踪 跟踪的目的是将写入数据文件的细化处理后的栅格数据，整理为从结点出发的线段或闭合的线条，并以矢量形式存储于特征栅格点中心的坐标。,11/14/2024,21,III、跟踪 跟踪的目的是将写入数据文件的细化处理后的栅格,弧段数据自动生成多边形,11/14/2024,22,弧段数据自动生成多边形9/21/202322,第三节 空间数据的融合,1、,遥感与GIS数据的融合,1）遥感图像与DLG图形的融合,丰富的光谱信息、几何信息、行政界线和属性信息，提高可视化效果。,2）遥感数据与DEM的融合,有助于遥感影像的几何校正与配准。,3）遥感与地图扫描图像的融合,发现快速变化的区域，进行自动更新。,11/14/2024,23,第三节 空间数据的融合1、遥感与GIS数据的融合9/21/,2、不同格式数据的融合,矢量 交换格式文件 栅格,ARC/INFO Coverage E00 Grid,ARCVIEW Shape Grid,Mapinfo Tab Mif,AutoCAD DWG Dxf,MapGIS .wt .wl .wp,Microstation DGN,11/14/2024,24,2、不同格式数据的融合 矢量 交,融合方式：,基于转换器,基于数据标准,基于公共接口,基于直接访问,11/14/2024,25,融合方式：9/21/202325,第四节 压缩与重分类,1、数据压缩的意义,优化存储空间，减少处理时间,2、数据压缩：从所取得的数据集合中抽取一个子集，作为一个新的信息源，在规定的精度范围内最好地逼近原集合，而又取得尽可能大的压缩比。a=m/n1,式中：,m,为曲线的原点数；,n,为曲线经压缩后的点数。,1）曲线上点的压缩,道格拉斯法(Douglas-Peucker)；,2）栅格数据的压缩,游程编码法；四叉树编码压缩法。,3,6,5,7,8,M,(1),N,(2),曲线上点的筛选,11/14/2024,26,第四节 压缩与重分类1、数据压缩的意义36578M(1)N,二 空间数据重分类,11/14/2024,27,二 空间数据重分类9/21/202327,离散空间：,空间具有跳跃特征（土地利用类型），重要变化发生在边界上，边界内的变化则是,均匀的,，,同质的，,即在各个方面都是相同的。,连续空间：,空间具有渐变特征（地形表面），内插技术必须采用连续的空间渐变模型实现这些连续变化，可用一种平滑的数学表面加以描述。这类技术可分为,整体拟合,和,局部拟合技术,两大类。,第五节 空间数据内插方法,11/14/2024,28,离散空间：空间具有跳跃特征（土地利用类型），重要变化发,整体拟合技术:,拟合模型是由研究区域内所有采样点上的全部特征观测值建立的。通常采用的技术是整体趋势面拟合。这种内插技术一般用于模拟大范围内的变化。,局部拟合技术：,是仅仅用邻近的数据点来估计未知点的值，而不受局部范围外其它点的影响。这类技术包括双线性多项式内插、样条函数、移动拟合法等等。,11/14/2024,29,整体拟合技术:拟合模型是由研究区域内所有采样点上的全部,1 分块内插,（1）线性内插,（2）双线性多项式内插,2 逐点内插,（1）移动拟合法,（2）加权平均法,3 数据精度分析,11/14/2024,30,1 分块内插9/21/202330,二区域内插,1 叠置法,2 比重法,11/14/2024,31,二区域内插9/21/202331,1,2,3,11/14/2024,32,1239/21/202332,比重法,11/14/2024,33,计算机图形学演示稿 纪玉波制作(C),比重法9/21/202333计算机图形学演示稿 纪玉波制作,11/14/2024,34,计算机图形学演示稿 纪玉波制作(C),9/21/202334计算机图形学演示稿 纪玉波制作(C),第6节 空间拓扑关系编辑,11/14/2024,35,计算机图形学演示稿 纪玉波制作(C),第6节 空间拓扑关系编辑9/21/202335计算机图形学,11/14/2024,36,计算机图形学演示稿 纪玉波制作(C),9/21/202336计算机图形学演示稿 纪玉波制作(C),作业（2）,11/14/2024,37,计算机图形学演示稿 纪玉波制作(C),作业（2）9/21/202337计算机图形学演示稿 纪玉波,一 比重法可以根据平滑密度函数原理将源区人口统计数据从同质性改变为非同质，其步骤如下：,在源区上叠置一张格网,将源区各个分区的平均人口数赋予相应分区的各个格网点。,.,、步结果如图，请对,粗线区,进行平滑计算,11/14/2024,38,计算机图形学演示稿 纪玉波制作(C),一 比重法可以根据平滑密度函数原理将源区人口统计数据从同质性,二 TIN中某三角形如图所示，其坐标（X，Y，Z）如下：,A（150,160,20）B（150,210,30）,C（180,180,40）P（160,180,?）,请用线性内插法计算P点高程Z。,11/14/2024,39,计算机图形学演示稿 纪玉波制作(C),二 TIN中某三角形如图所示，其坐标（X，Y，Z）如下：9/,三 图1为源