第4章空间数据的处理.ppt

利用局部范围内的已知采样点的数据内插出未知点的数据。 1、线性内插 将内插点周围的3个数据点的数据值带入多项式，即可解算出系数a0、a1、a2 。 2、双线性多项式内插 将内插点周围的4个数据点的数据值带入多项式，即可解算出系数a0、a1、a2、a3 。 当数据是按正方形格网点布置: 三、局部内插 3、双三次多项式（样条函数）内插 是一种分段函数,每次只用少量的数据点，故内插速度很快；样条函数通过所有的数据点，故可用于精确的内插；可用于平滑处理。 ??? 双三次多项式内插的多项式函数为： ?? 将内插点周围的16个点的数据带入多项式，可计算出所有的系数。 16个点 三、局部内插 4. 移动平均法 在局部范围（或称窗口）内计算个数据点的平均值. 二维平面的移动平均法也可用相同的公式，但位置Xi应被坐标矢量Xi代替。 窗口的大小对内插的结果有决定性的影响。小窗口将增强近距离数据的影响；大窗口将增强远距离数据的影响，减小近距离数据的影响。 当观测点的相互位置越近，其数据的相似性越强；当观测点的相互位置越远，其数据的相似性越低。 加权移动平均法:λi是采样点i对应的权值 加权平均内插的结果随使用的函数及其参数、采样点的分布、窗口的大小等的不同而变化。通常使用的采样点数为6—8点。 对于不规则分布的采样点需要不断地改变窗口的大小、形状和方向，以获取一定数量的采样点。 三、局部内插 4.6 数据的压缩与光滑 一、数据压缩 二、曲线光滑 图形显示输出 数据存储 数据压缩 光滑 矢量数据压缩 通常是对线状实体的压缩。 最常用的是道格拉斯-佩克算法 压缩效果好，但必须在对整条曲线数字化完成后才能进行，且计算量较大 一、数据压缩 栅格数据压缩 可采用游程编码和四叉树编码等方法。 一、数据压缩 是假象曲线为一组离散点，寻找形式较简单、性能良好的曲线解析式。 插值方式：曲线通过给定的离散点。如拉格朗日插值，三次样条曲线 逼近方式：曲线尽量逼近给定离散点。如贝塞尔和B样条曲线。 二、曲线光滑 4.7 空间数据结构的转换 一、矢量向栅格的转换 二、栅格向矢量的转换 点的栅格化 坐标转换 线的栅格化 DDA法(数字微分分析法) Bresenham算法 面的栅格化 线的栅格化+面(多边形)的填充。面的填充法有： 内部点扩散法（种子扩散法） 扫描法 边填充算法 一、矢量向栅格的转换 要求（矢量化过程应保持）： 拓扑转换，即保持实体原有的连通性、邻接性等； 转换实体保持正确的外形。 方法 实际应用中大多数采用人工矢量化法，如扫描矢量化，但工作量大，成为GIS数据输入、更新的瓶颈问题之一。 程序转化转换（全自动或半自动） 1、边界提取 2、二值化 3、二值图像的预处理 4、细化:1）剥皮法 2)骨架法 5、跟踪 6、拓扑化 遥感影象图 栅格分类图 边界提取 二值化 编辑 矢量跟踪 数据压缩 原始线划图 二值化 细化 分类图 扫描 预处理 拓扑化 二、栅格向矢量的转换（矢量化） 第4章 结 束 第4章 空间数据的处理 本章内容 4.1 坐标变换 4.2 图形编辑 4.3 拓扑关系的自动建立 4.4 拓扑关系的编辑 4.5 空间插值 4.6 数据压缩与光滑 4.7 空间数据结构的转换 4.1 坐标变换 一、几何变换（几何纠正） 比例尺变换 乘系数 变形误差改正 通过控制点利用高次变换、二次变换和仿射变换加以改正 坐标旋转和平移 即数字化坐标变换，利用仿射变换改正。 二、投影变换 解析变换 数值变换 解析—数值变换 一、几何纠正  其中A、B代表二次以上高次项之和。上式是高次曲线方程，符合上式的变换称为高次变换。式中有12个未知数，所以在进行高次变换时，需要有6对以上控制点的坐标和理论值，才能求出待定系数。 1、高次变换 2、二次变换 当不考虑高次变换方程中的A和B时，则变成二次曲线方程，称为二次变换。二次变换适用于原图有非线性变形的情况，至少需要5对控制点的坐标及其理论值，才能解算待定系数。 3、仿射变换 实质是两坐标系间的旋转变换。设图纸变形引起x,y两个方向比例尺不同，当x,y比例尺相同时，为相似变换。 特性： 直线变换后仍为直线； 平行线变换后仍为平行线； 不同方向上的长度比发生变化。 求解上式中的6个未知数，需不在一直线上的3对已知控制点，由于误差，需多余观测，所以，用于图幅定向至少需要四对控制点。 一、几何纠正 二、地图投影变换 1、解析变换法 1）反解变换法(又称间接变换法) 假定原图点的坐标为x,y(称为旧坐标)，新图点的坐标为X，Y(称为新坐标)，则由旧坐标变换为新坐标的基本方程式为： 2）正解变换法(又称直接变换法) 2、数值变