第三章 间数据处理.ppt

分块内插法（双线性多项式内插法） 含义：是在规则分布的已知数据点时，使用最靠近内插点的四个已知数据点组成一个四边形，确定一个双线性多项式来内插其中点的高程。 设所求的函数形式为： Z=a0+a1x+a2y+a3xy 由内插点周围的四个已知数据点的坐标代入计算求得待定系数a0、a1、a2、a3。 分块内插法（二元样条函数内插法—双三次多项式） 含义：在分块插值区用双三次多项式即样条函数模拟地表面。 特点：是一种分段函数，每次只用少量的数据点，故内插速度很快；样条函数通过所有的数据点，故可用于精确的内插；可用于平滑处理。 设所求的函数形式为： Z=a1x3y3+a2x2y3+a3xy3+a4y3+a5x3y2+a6x2y2+a7xy2+a8y2+a9x3y+a10x2y+a11xy+a12x3+a13x3+a14x2+a15x3+a16 将内插点周围的16个点的数据带入多项式，求出待定系数a1-a16。 16个点 以待定点为中心进行插值。 被插值点 逐点内插法 移动拟合法是典型的逐点内插法。对每一个待定点取用一个多项式曲面拟合该点附近的地表面，此时取待定点作平面坐标的原点，并用待定点为圆心，以R为半径的圆内诸数据点来定义函数的待定系数。 反比距离加权法（IDW法） Inverse Distance Weighted 通用，简单 以距离的倒数为权重，距离越小、影响越大。 常用与地形生成。 Kriging法 Kriging法：点与点之间在距离、方向上存在相互影响关系，可以预先知道 计算量大，取样地点有限制， 用于土壤成份、地质构造。 3、数据精度分析 在研究区内随即选取n个数据点，设这些点的图上高程值为zi，其对应的插值数据值为z`i，它们之差用△zi＝ zi – z`i表示。则： 算术平均值为 标准差为 △zi＜原始数据比例尺等高距的1/3 二、区域的内插 区域内插是根据一组多边形分区的已知数据来推求同一地区另一组多边形分区未知数据的内插方法。 1、叠置法 A C B 1 2 3 2、比重法 比重法是根据平滑密度函数的原理，将源区的统计数据从同性质改变为非同性质，而非同性质代表着一般社会经济现象的普遍特点。 第六节 空间拓扑关系的编辑 弧段号 左多边形 右多边形 起始节点 终止节点 a O A 1 8 b A O 1 2 c O B 3 2 d O B 4 3 e O B 5 4 f O C 6 5 g O C 7 6 h O C 8 7 i C A 8 9 j C B 9 5 k C D 12 10 l C D 11 12 m C D 10 11 n B A 9 2 一、多边形连接编辑 设需要对多边形A进行编辑，其算法过程为： 1、从上表的弧段文件中，找出与当前编辑的多边形A相关的所有记录。 弧段号 起始节点 终止节点 左多边形 右多边形 a i n b 1 8 9 1 8 9 2 2 O C B A A A A O 2、在找出的记录中，检查当前编辑的多边形A所处的位置：如果A处在左多边形位置，将之与处于右多边形位置的多边形号相交换，同时也将该记录的节点号位置作对应的交换；反之，如果当前编辑的多边形A处于右多边形位置，则该记录的所有数据项顺序不作改变。 弧段号 起始节点 终止节点 左多边形 右多边形 a i n b 1 8 9 2 8 9 2 1 O C B O A A A A 3、从经过代码位置转换的记录中，任取一个起始节点作为起点，顺序连接各个节点，必要时可对记录的前后顺序作调整，使得连接的节点能自行封闭，即1→8→9→2→1。 弧段号 起始节点 终止节点 左多边形 右多边形 a i n b 1 8 9 2 8 9 2 1 O C B O A A A A 二、节点连接编辑 设需要对节点8进行编辑，其算法过程为： 1、从表的弧段文件中，找出与当前编辑的节点8相关的所有记录。 弧段号 起始节点 终止节点 左多边形 右多边形 a i h 1 8 8 8 9 7 O C O A A C 2、在找出的记录中，检查当前编辑的节点8所在的位置：如果8处在起始节点位置，将之与处在终止节点位置的节点号交换，同时也将该记录的多边形号位置作相应的交换；反之，如果当前编辑的节点8处于终止节点位置，则该记录的所有数据项顺序不作改变。 弧段号 起始节点 终止节点 左多边形 右多边形 a i h 1 9 7 8 8 8 O A C A C O 3、从经过代码位置转换的记录中，任取一个左多边形作为起点，顺序连接各个多边形，同样，必要时可对记录的前后顺序作调整，使得连接的多边形能首尾呼应，即O →A→C→O。 弧段号 起始节点 终止节点 左多边形 右多边形 a i h 1 9 7 8 8