第3章空间数据处理技巧.ppt

第 3 章 空间数据处理;数据变换指数据从一种数学状态到另一种数学状态的变换，包括几何纠正和地图投影转换等，以实现空间数据的几何配准。 数据重构指数据从一种格式到另一种格式的转换，包括结构转换、类型替换等，以实现空间数据在结构、格式和类型上的统一，多源和异构数据的连接与融合。 数据提取指对数据进行某种条件的取舍，包括类型提取、窗口提取、空间内插等，以适应不同用户对数据的特定要求。;第 3 章 空间数据处理;第1节 空间数据的变换 第2节 空间数据结构的转换 第3节 多元空间数据的融合 第4节 空间数据的压缩与重分类 第5节 空间数据的内插方法 第6节 空间拓扑关系的编辑;第1节 空间数据的变换;一、几何纠正;扭曲 x*=x.Sx y*=y.Sy 其中：Sx,Sy分别为x, y方向 的缩放比例系数，矩阵为： [x*, y*]=[x, y].;比例缩放 x*=x.Sx y*=y.Sy 其中：Sx,Sy分别为x, y方向 的缩放比例系数，矩阵为： [x*, y*]=[x, y].;旋转 将点(x, y)旋转?角 x=A.cos?, y=A.sin? x*=A.cos(?+ ?) =A.(cos?.cos? -sin?.sin?) =x.cos? -y.sin? y*=A.sin(? + ?) =A(sin?.cos? +cos?.sin?) =x.sin? +y.cos? 矩阵为： [x*, y*]=[x, y]. ;平移 x*=x+dx y*=y+dy 其中，dx、dy分别为x, y 方向平移量 相应的向量形式为： [x*, y*]=[x, y]+[dx, dy]; 设x，y为数字化仪坐标，X、Y为理论坐标，m1、m2为地图横向和纵向的长度变化比例，两坐标系夹角为α ，数字化仪原点O’相对于理论坐标系原点平移了a0、b0，则根据图形变换原理，得出坐标变换公式： X = a0 +(m1cosα)x + (m2sinα)y Y = b0–(m1sinα)x + (m2cosα)y;; X = a0 +(m1cosα)x + (m2sinα)y Y = b0–(m1sinα)x + (m2cosα)y 式中，设 a1 = m1cosα ， b1 = -m1sinα a2 = m2sinα ， b2 = m2cosα 则上式可以简化为： X = a0 + a1 x + a2 y Y = b0 – b1x + b2 y 上式中含有6个参数a0、a1、a2、b0、b1、b2，要实现仿射变换，需要知道不在同一直线上的3对控制点的数字化坐标及其理论值，才能求得上述6个待定参数。但在实际应用中，通常利用4个以上的点来进行几何纠正。下面按最小二乘法原理求解待定参数（首先回顾最小二乘法）; 设Qx、Qy表示转换坐标于理论坐标之差，则有： Qx = X – (a0 + a1x + a2y) Qy = Y – (b0 + b1x + b2y) 按照[Qx2] = min和[Qy2] = min的条件，可得到两组法方程： a0n + a1∑x + a2∑y = ∑X a0∑x + a1∑x2 + a2∑xy = ∑xX a0∑y + a1∑xy+ a2∑y2= ∑yX 和 b0n + b1∑x + b2∑y = ∑Y b0∑x + b1∑x2 + b2∑xy = ∑xY b0∑y + b1∑xy+ b2∑y2= ∑yY 式中：n为控制点个数； x、y为控制点的数字化坐标； X、Y为控制点的理论坐标。; 由上述法方程，通过消元法，可求得仿射变换的待定参数a0、a1、a2和b0、b1、b2。 仿射变换是GIS数据处理中使用最多的一种几何纠正方法。它的主要特性为：同时考虑到x和y方向上的变形，因此纠正后的坐标数据在不同方向上的长度比将发生变化。其他方法还有相似变换和二次变换等。 经过仿射变换的空间数据，其精度可用点位中误差表示，即：;二、地图投影及其转换;;第2节 空间数据结构的转换;第2