第七讲空间数据处理预览.ppt

二、空间插值方法 滞后距（LAG）的影响 二、空间插值方法 6.4 理论变异函数 块金值（C0； Nuddget): 当h = 0时的非零变异值，由不可解释的原因引起； 基台值（C0+C；sill): 变异曲线变平缓时的变异值，表明在某个距离上样本点不再存在相关性，通常等于数据集的方差； 变程值（range）: 当基台值出现时的h值. 重要的是原点附近变异图的形状，越是最近的点对插值结果的影响越大 二、空间插值方法 kriging插值方法 式中：Z0为待插入点的值； Zx为已知点的值 Wx为每个点的权重值 计算Wi，按采样点数据的变异函数图的统计分析原理来计算 kriging权重计算方法 h11 r(h11) 式中：r(hij)是已知点i和j之间的变异值，r(hi0)是已知点i和待估点之间的变异值， r(hij)和r(hi0)可以通过变异函数求得； Wi为每个已知点的权重，u是拉格朗日系数，是为了将估计误差降低到最小； 根据上式，通过解算联立方程可以计算出Wi 习题 构建泰森多边形 Kriging估值 块金 变程 基台 1.0 110(m) 3.1 泛克里格内插生成预测图 第五节 空间数据的内插方法 主要方法： 叠置法 比重法 区域内插主要解决离散属性数据问题。研究的目标是从已知分区的数据(如社会经济数据)中推出同一地区的另一组分区数据。 如： 已知某地区各县中历年的人口分布数据，但因行政区划分使该地区中某些县的边界线发生了变化，现在需要推算新行政区中历年的人口分布数据，就可用这种插值方法。 二、区域内插 人口 面积 A 35 7 B 30 6 C 10 3 A C B C1 A1 B1 历史区 现实区 第五节 空间数据的内插方法 主要过程： (1)在历史区上叠加满足精度的格网，将历史区内的各格网赋予平均。 (2)按邻域法平滑数据，计算公式可用四邻域或八邻域法。 四邻域法公式：zij=(zi+1, j+ zi-1, j+ zi, j+1 + zi, j -1) ／4 八邻域法公式： zij=(zi+1, j+ zi-1, j+ zi, j+1 + zi, j -1 + zi+1, j+1+ zi+1, j-1+ zi-1, j+1 + zi-1, j -1 ) ／8 区域内插方法-比重法 比重法使用平滑密度函数 第五节 空间数据的内插方法 (3)按四邻域或八邻域法平滑数据，求平滑后区域内各网格值之和。如按四邻域法求U1A，U1B，U1C ，得到数据的变化率，检查是否符合要求。 PA= UA ／ U1A PB= UB ／ U1B PC= UC ／ U1C (4)若数据变化率不符合要求，各格网值乘以变化率，得到调整后的格网值，再进行第二次平滑。如此循环，直到区域数据的变化率满足要求。 第五节 空间数据的内插方法 例证：区域内插方法-比重法 已知某地区有A、B、C三个区域，并已知该三个区域的面积和人口数，人口数分别为UA、UB 、UC ，如下图所示。现将该区域重新划分A1、B1 、C1三个区域，如图，请用区域属性数据的插值法求A1、B1 、C1三个区域的面积和人口数。 人口 面积 A 35 7 B 30 6 C 10 3 A C B C1 A1 B1 历史区 现实区 第五节 空间数据的内插方法 (1)在历史区上叠加满足精度的格网。将历史区内的各格网赋予平均值。假设整个区域分为4x4栅格，求出： A区每个栅格值为 35／7 = 5.0 B区每个橱格值为 30／6 = 5.0 C区每个栅格值为 10／3 = 3.3 如图所示 人口 面积 A 35 7 B 30 6 C 10 3 5.0 5.0 5.0 5.0 5.0 5.0 3.3 3.3 5.0 5.0 5.0 3.3 5.0 5.0 5.0 5.0 图 历史区 A B C 第五节 空间数据的内插方法 (2)对图 历史区数据做第一次平滑后，得图1-1。 5.0 5.0 5.0 5.0 5.0 3.3 3.3 5.0 5.0 5.0 3.3 5.0 5.0 5.0 5.0 5.0 5.0 4.4 4.2 5.0 4.6 4.6 3.9 5.0 5.0 4.2 4.4 5.0 5.0 5.0 4.2 注： 在进行邻域法平滑计算（zij=(zi+1, j+ zi-1, j+ zi, j+1 + zi, j -1) ／4）时，当点在边缘时，应除2或3。 图 历史区数据 图1-1 第一次平滑数据 4.6 5.0 5.0 5.0 3.3 5.0 5.0 5.0 5.0 3.3 4.6 + + + = ( ) / 4 z22 = 第