沈阳农业大学土地与环境学院土地信息系统课件 第五章.ppt

半变率图的组成 块金值（c0； Nuddget): 当h = 0时的非零变率，由不可解释的原因引起 ； 基台值（c0+c；sill): 半变率曲线变平缓时的变率值，表明在某个距离上样本点不再存在相关性，通常等于数据集的方差； 变差值（range）: 当基台值出现时的h值 (sill 95% 时的h值). 重要的是原点附近半变率图的形状，越是最近的点对插值结果的影响越大。 空间自相关部分：C/(c0+c) 半方差模型 球面模型 指数模型 线性模型 高斯模型 3、Kriging插值的方法 式中：z0为待插入点的值； Zx为已知点的值 Wx为每个点的权重值 计算Wi，按采样点数据的半方差图的统计分析原理来计算。 权重的计算（三个已知点情况下） 式中：r(hij)是已知点i和j之间的半方差，r(hi0)是已知点i和未知点之间的半方差， r(hij)和r(hi0)可以通过半方差模型来求得； Wi为每个已知点的权重， 是拉格朗日系数，是为了将估计误差降低到最小； 根据上式，通过解算联立方程可以计算出Wi 待插值的计算 Z0=Z1W1+Z2W2+Z3W3 并且通过方差测度来计算误差： 4、ArcGIS下的Kriging插值 5、IDW vs. Kriging Kriging 产生似乎更自然的结果，避免异常值的产生； 同时能给出标准误差。 其它kriging插值 通用克里金插值：要求数据是二阶平稳的或纯平稳的 泛克里金插值：如果数据在空间上存在明显的趋势，那么，应该使用泛克立金方法进行分析 块克里金插值：对中心在X0的小区或块段进行估值 协克里金插值：协同克立格分析是一种空间数据的解释技术，其基本的思想是利用变量之间的空间相关关系进行估计或预测。按照如下的统计准则，它可能是目前最好的空间数据分析方法 要处理的问题是：有两个空间变量，如果第二个变量的分布广，采样密度更高，而第一个变量难以测定或测定的费用较高，那么，可以利用有限样本的变量之间的空间关系来改进对于第一个变量的估计。 ? ESRI ARCGIS中的应用 1.调入“地统计分析“工具 2.参数设置 分析图层及属性设置 插值方式 3.参数设置 误差统计 4.插值结果图 由于插值图像为矩形输出，需重新设置其填充边界。 未填充的区域 填充过多的区域 离散的采样点数据 连续表面 5.裁减边界设置 Click 勾选”启用“ 6.裁减边界设置 选择裁减边界 7.插值图的裁减 8.初步插值图 9.插值图的调整 （1）图例调整 9.插值图的调整 （1）图例调整 选择分级数目 选择手工输入方式 输入级别间距 最后点OK 9.插值图的调整 （1）图例调整 调整分级后的插值图 由于打印中经常采用黑白打印，故需调整图例颜色，使对比度明显。 单击，在出现的对话框中选择颜色或线条 9.插值图的调整 （1）图例调整 成果图 （2）插入指北针 选择指北针 选择指北针，并点击OK。 指北针插入成功。 （3）插入比例尺 点击进行调整修改，直到满意，点OK。 比例尺插入成功。 距离反比权重系数的确定 ARCGIS——IDW 控制反距离加权的参数—权重 权重过高，较近点的影响较大，拟合表面更细致（不光滑）； 权重过低，较远点的影响增加，拟合表面更光滑。缺省值常为 2 。 控制反距离加权的参数—搜索半径 搜索半径－固定 对固定型半径，搜索距离一定，所有在该半径内的样点参与计算。 可预先设定一个阈值，当给定半径内搜索到的点小于该值时可扩大搜索半径，直到达到该阈值为止。 搜索半径－可变 设定参与计算的样点数是固定的，则搜索的半径是可变的。这样对每个插值点的搜索半径可能都不同，因为要达到规定的点数所需要搜索的区域是不一样的。 控制反距离加权的参数—障碍设置 可利用一线状和面状数据集来限制样点的搜索。 线状数据集可作为平坦地表的悬崖或脊状障碍物： 只有位于同侧的样点才符合要求。 权重系数和搜索半径的影响图示 Power = 2, search = 230 Power = 2, search = 150 Power = 2, search = 600 Power = 4, search = 600 距离反比插值评价 优点——简便易行；可为变量值变化很大的数据集提供一个合理的插值结果；不会出现无意义的插值结果而无法解释。 不足——对权重函数的选择十分敏感；易受数据点集群的影响，结果常出现一种孤立点数据明显高于周围数据点的“鸭蛋”分布模式； 全局最大和最小变量值都散布于数据之中。 距离反比很少有预测的特点，内插得到的插值点数据在样点数据取值范围内。 四、高次曲面插值(Multiquadric) 高次曲面插值由 Hardy 于1971年首先提出，随后应用于不同的学科。每个