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

1、2021-11-41 具体转换步骤分如下几步： 1 准备好矢量数据或矢量地图； 2 构架一个与地图等大小的格网，选取好适当的格网密度或分辨率，依据范围则可以算出格网的行列数。 3 根据地图上的点、线、多边形相对这个格网的配置及其属性来确定每一个项元的数值。第1页/共40页2021-11-42 上述三步中，第三步内容较多，需要专门说明：1、点状对象的栅格化：点的变换十分简单，只要这个点落在哪个网格中就是属于那个网格元素，根据该点状对象的特性赋予该像元属性值。点行、列坐标i，j计算公式：第2页/共40页 2、线状对象的栅格化（一）： 假定一线段两端点之间经过若干个网格元素(至少一个)，两端点坐标为

2、(X1，Y1)，(X2，Y2)，则：P82YX（X1，Y1）（X2，Y2）第3页/共40页2021-11-44 判断行差和列差是为了确保单一连接：单一连接定义： 1 、线条不中断，线上每两个像元间必须连接，或邻边连接，或顶角连接。 2 、任何三个像元不得两两间都相邻。第4页/共40页2021-11-453、面的栅格化方法： 一、基于弧段数据的栅格化方法： 按行或列做中心扫描线 求交点栅格坐标的行列值 判断交点左右多边形数值 交点排序并逐段生成栅格数据 特点： 计算量较大，算法比较复杂。第5页/共40页2021-11-46abcdef10110111101110000111101001第6页/共

3、40页2021-11-47000100010000001000000001000011101110001001101100011100000001110110000100111010001100001111111110第7页/共40页2021-11-48二、基于多边形的栅格化方法：、内部点扩散法： 由一个内部的种子点，向其四个方向的邻点扩散，判断新加入的点是否在多边形边界上。如果是，不作为种子点；否则当作新的种子点，直到区域填满。算法特点： 算法设计复杂，而且可能造成阻塞而使扩散不能完成。第8页/共40页2021-11-49、边界代数法：上减下加特点：算法简单，但对于复杂图形，每一像素可能被访

4、问多次，增加了运算量。第9页/共40页2021-11-410 、包含检验法： 点在多边形内判断： 检验夹角之和 射线算法：由待判点向图外某点引射线，判断该射线与某多边形所有边界相交的总次数，如果相交偶数次，则待判点在该多边形外部，如为奇数次，则待判点在该多边形内部如图所示。 n=0内部点外部点n=2n=1n=3n=4n交点个数第10页/共40页2021-11-411三、栅格到矢量： 栅格向矢量转换处理的目的，是为了： 1 将栅格数据分析的结果，通过矢量绘图装置输出； 2 数据压缩的需要，将大量的面状栅格数据转换为由少量数据表示的多边形边界； 3 将自动扫描仪获取的栅格数据加入矢量形式的数据库。

5、第11页/共40页2021-11-412 1、基于图象数据的矢量化方法 图象数据是由不同灰阶的影像或线划，通过自动扫描仪(scanner)，按一定的分辨率进行扫描采样，得到以不同灰度值(0255)表示的数据。目前扫描仪的分辨率可达00125mm，因此对一般粗度(例如0.1mm)的线条，其横断面扫描后平均也有8个像元，而矢量化的要求只能允许横断面保持一个栅格的宽度，因此需要进行从栅格向矢量数据的转换。第12页/共40页2021-11-413 5 9 10141 138 9 5 3 1 0 2 245156 73 144 178 132 23 7 3 212 5 6 8 29 11214 167

6、5 124110 7 6 5 4 7 133 5 192 350 110 135 6 4 7 244 12 2 5 12135 201 166 127155 9 1 1 9 4 8 2112211 43 5 0 （a）扫描前的矢量数据 （b）扫描得到的灰度值第13页/共40页2021-11-414 具体转换的步骤： I 、二值化 线划图形扫描后产生栅格数据，这些数据是按从0255的不同灰度值量度的，设以G(i，j)表示，为了将这种256级不同的灰阶压缩到2个灰阶，即0和1两级，首先要在最大与最小灰阶之间定义一个阈值，设阈值为T，则如果G(i，j)大于等于T，则记此栅格的值为1，如果G(i，j)

7、小于T，则记此栅格的值为0，得到一幅二值图。第14页/共40页2021-11-415根据给定的阈值二值化后得到的栅格数据第15页/共40页2021-11-416 II、细化 细化是消除线划横断面栅格数的差异，使得每一条线只保留代表其轴线或周围轮廓线(对面状符号而言)位置的单个栅格的宽度。 第16页/共40页2021-11-417 对于栅格线划的“细化”方法，常用 “剥皮法” 。剥皮法的实质是从曲线的边缘开始，每次剥掉等于一个栅格宽的一层，直到最后留下彼此连通的由单个栅格点组成的图形。因为一条线在不同位置可能有不同的宽度，故在剥皮过程中必须注意一个条件，即不允许剥去会导致曲线不连通的栅格。 第1

8、7页/共40页2021-11-418 这是这一方法的技术关键所在。其解决办法是，借助一个在计算机中存储着的，由待剥栅格为中心的33栅格组合图来决定。一个33的栅格窗口，其中心栅格有八个邻域，因此组合图有多种不同的排列格式，若将相对位置关系的差异只是转置90、180、270度或互为镜象反射的方法进行归并，则共有51种排列格式。 123456789101112131415161718282930313233343536242526271920212223464748495051424344453738394041第18页/共40页2021-11-41912345678910111213141516

9、1718282930313233343536242526271920212223464748495051424344453738394041 显然，其中只有格式2、3、4、5、10、11、12、16、21、24、28、33、34、35、38、42、43、46和50，可以将中心点剥去。这样，便可确定中间栅格点保留或删除，得到经细化处理后应予保留的栅格系列，并写入数据文件。第19页/共40页2021-11-420III、跟踪 跟踪的目的是将写入数据文件的细化处理后的栅格数据，整理为从结点出发的线段或闭合的线条，并以矢量形式存储于特征栅格点中心的坐标。 第20页/共40页2021-11-421第21

10、页/共40页2021-11-422第三节第三节 空间数据的融合空间数据的融合1、遥感与GISGIS数据的融合1 1）遥感图像与DLGDLG图形的融合 丰富的光谱信息、几何信息、行政界线和属性信息，提高可视化效果。2 2）遥感数据与DEMDEM的融合 有助于遥感影像的几何校正与配准。3 3）遥感与地图扫描图像的融合 发现快速变化的区域，进行自动更新。第22页/共40页2021-11-4232 2、不同格式数据的融合 矢量 交换格式文件 栅格ARC/INFO Coverage E00 GridARCVIEW Shape GridMapinfo Tab Mif AutoCAD DWG DxfMapG

11、IS .wt .wl .wpMicrostation DGN第23页/共40页2021-11-424融合方式： 基于转换器 基于数据标准 基于公共接口 基于直接访问第24页/共40页2021-11-425第四节第四节 压缩与重分类压缩与重分类1 1、数据压缩的意义优化存储空间，减少处理时间2 2、数据压缩：从所取得的数据集合中抽取一个子集，作为一个新的信息源，在规定的精度范围内最好地逼近原集合，而又取得尽可能大的压缩比。 a= a= m/n1 m/n1 1 1）曲线上点的压缩2 2）栅格数据的压缩36578M(1)N(2)曲线上点的筛选第25页/共40页2021-11-426二二 空间数据重分

12、类空间数据重分类第26页/共40页2021-11-427 离散空间：空间具有跳跃特征（土地利用类型），重要变化发生在边界上，边界内的变化则是均匀的，同质的，即在各个方面都是相同的。 连续空间：空间具有渐变特征（地形表面），内插技术必须采用连续的空间渐变模型实现这些连续变化，可用一种平滑的数学表面加以描述。这类技术可分为整体拟合和局部拟合技术两大类。第五节第五节 空间数据内插方法空间数据内插方法 第27页/共40页2021-11-428 整体拟合技术:拟合模型是由研究区域内所有采样点上的全部特征观测值建立的。通常采用的技术是整体趋势面拟合。这种内插技术一般用于模拟大范围内的变化。 局部拟合技术：

13、是仅仅用邻近的数据点来估计未知点的值，而不受局部范围外其它点的影响。这类技术包括双线性多项式内插、样条函数、移动拟合法等等。第28页/共40页2021-11-429 1 分块内插 （1）线性内插 （2）双线性多项式内插 2 逐点内插 （1）移动拟合法 （2）加权平均法 3 数据精度分析第29页/共40页2021-11-430二区域内插1 叠置法2 比重法第30页/共40页2021-11-431123第31页/共40页比重法第32页/共40页第33页/共40页第6节 空间拓扑关系编辑第34页/共40页第35页/共40页作业（2）第36页/共40页一 比重法可以根据平滑密度函数原理将源区人口统计数据从同质性改变为非同质，其步骤如下： 在源区上叠置一张格网 将源区各个分区的平均人口数赋予相应分区的各个格网点。 .、步结果如图，请对 粗线区 进行平滑计算第37页/共40页二 TIN中某三角形如图所示，其坐标（X，Y，Z）如下：A （150,160,20） B（150,210,30） C （180,180,40） P （160,180,?）请用