克里金实验报告

1、实验报告姓名：李丹学号：201521170071 专业:摄影测量与遥感一、作业要求对观测数据进行变差分析，并采用课堂所讲的插值方法对其进行插值运 算，随后进行验证。二、实验过程1 .数据介绍本次试验所采用的是1986年5月8日于瑞士 467个降水观测站的降水数 据，为数据集中的“ SIC97-full.xls”。原始数据在空间中的分布如下图所示。 其中颜色越深表示降水量越多。图1原始数据在空间中的分布2 .数据分析在ArcMap中打开上述数据集，进行数据分析频数直方图是指用长方形的高度表示某一段数据出现的频数，通过频数直 方图可以清楚显示各组频数分布情况，并且能够显示各组之间频数的差别。作 出

2、原始数据的频数直方图如下所示频数，10二二0.B9.0. G7值值值差 数小大均准 总最最平标偏度：0.65勺3曜度:27964U分位薮：99-25 市梅蠢：15i 4分位数:262.754422Z Z. 5E L. 17 L. 7b 2. 342. 933.4. 094. S3 .1 35图2频数直方图从上图中可以发现，原始数据的分布基本符合正态分布，可以进行克里金 插值运算。（2）正态QQ图正态Q-Q图就是由标准正态分布的分位数为纵坐标，样本值为横坐标的散 点图。通过该图可以鉴别样本数据是否近似于正态分布，根据样本数据作图如 下。图3正态QQ图从上图中可以看书，给出的点近似分布在一条直线附

3、近，因此可以判定原 始图像为正态分布。（3）趋势分析“趋势分析”工具提供数据的三维透视图。采样点的位置绘制在x,y平 面上。在每个采样点的上方，值由z维中的杆的高度给定。“趋势分析”工具 的唯一功能是值将会作为散点图投影到x,z平面和y,z平面上，可以将其视 为通过三维数据形成的横向视图，之后会根据投影平面上的散点图拟合多项式， 如图中的绿色和蓝色线条所示。（a） XZ平面投影（b） YZ平面投影从上图中可以看出，无论是在XZ平面投影还是在YZ平面投影，两者都具 有确切上升的模式，因此说明数据中存在某种趋势。（4）半方差和协方差分析半方差函数和协方差函数将邻近事物比远处事物更相似这一假设加以量

4、化。 半方差函数和协方差都将统计相关性的强度作为距离函数来测量。作出对应的 函数图像如下图5半方差函数图:0.43： S3 1.3515.032. 55. 912. 22.隹一霆j10图6协方差函数图从上图中可以看出，在短距离内基本符合距离越近越相似的原则。3.克里金插值克里金插值首先考虑的是空间属性在空间位置上的变异分布，确定对一个 待插点值有影响的距离范围，然后用此范围内的采样点来估计待插点的属性值。 该方法在数学上可对所研究的对象提供一种最佳线性无偏估计（某点处的确定值） 的方法。它是考虑了信息样品的形状、大小及与待估计块段相互间的空间位置等几何特征以及品位的空间结构之后，为达到线性、无

5、偏和最小估计方差的估 计，而对每一个样品赋予一定的系数，最后进行加权平均来估计块段品位的方 法。但它仍是一种光滑的内插方法在数据点多时，其内插的结果可信度较高。2.371. 7771. 185 . 592由于平均值未知，因此采用普通克里金插值方法。对于半方差的拟合采用 球形拟合。如下图所示0 m Model Eiimed Averaged图7半方差拟合方法得到的预测误差如下表所示最终得到的插值结果如下图所示。题I点0-66-66-10410+- 125133图8插值结果表1预测误差表项目取值平均值-0.05均方根48.21标准平均值0.00标准均方根1.12平均标准误差42.214.交叉验证交

6、叉验证使用所有数据对趋势和自相关模型进行估计。它会每次移除一个 数据位置，然后预测关联的数据值，然后对第二个点重复此过程，以此类推。 交叉验证会对所有点的测量值和预测值进行比较，完成交叉验证后，一些数据 位置可能被作为异常搁置，这时需要重新拟合趋势模型和自相关模型。已预测E.S55. 1153.644 2.9082. 173 1.437 0.702-0. 034 0.621.2741.9272.5813.2353.8894.5425. 1965.85加璋-10-2图9预测结果上图给出了预测值与真实值的散点图，可以看出它们分布在1:1线（图中灰色线）附近，但是他的坡度（蓝色线）小于1，这也是克里

7、金插值的一个 特点，即大值预测过低，小值预测过高。回归方程为y = 0.89 x x + 17.782.5210.B51.31. 952.53,药94, 5S S3 5.35已删重-W21.2B1-1.2611.891-1.S91-2.521误差图与预测图相同，只是在预测值基础上减去了测量值，从上图中可以看出，误差在零值附近波动。标准化误差4. ITT2.7S5-1.392-2.7S5-4.ITT-5.5700. 651. 31.952. 63.253. 94.555. 25.S5已贝悝-io-2图11标准化误差图预测图、误差图、标准化误差图都显示了克里金法预测的准确程度。标准化误差图要先用预测值减去测量值，然后再除以估计的克里金标准误差。如 果所有数据都是独立的，那么蓝线应该应该是水平的，从上图中可以看出，蓝 线基本接近水平，因此说明使用的克里金模型较为理想。标准化误差5.574 2042. 83B1. 4720. 106-1.26-2. 626-5. 358-3.07 -2.383 -1.705 -1. 023 -0. 341 0. 3411. 0231. 7052.3033. 07标准值图12正态QQ图该图