空间预测.PDF

空间预测 金百锁 2017 年 4 月23 日 1. 空间预测 空间预测公式 空间预测是对未知量Z s0 进行预测，它是根据样本数 据Z s 和关于Z 的趋势和方差 以及空间相关性的假设进行预  测。 线性预测公式： ˆ ˆ 1 ˆ Z s0  s0 β  V Z s Xβ 其 中V为Z s的协方差矩阵，是Z s和Z s0 的协方差 向 ˆ 1 1 1 量。β X V X X V Z s是广义最小二乘估计。 2 ˆ 2 预测误差σ s0 EZ s0 Z s0 1 1 1 rZ s0  V  δX V X δ , 这 1 里δ  s0  V X . 3/24 空间预测 简单 （simple ），普通(ordinary)，泛(universal)克 里金 （kriging）决策树 4/24 空间预测 Meuse河堤数据 library (sp) data (meuse) coordinates(meuse) - c (x, y) data (meuse.grid ) coordinates(meuse.grid ) - c (x, y) meuse.grid - as (meuse.grid , +SpatialPixelsDataFrame) 5/24 空间预测 估计变异函数 library (gstat) f - log (zinc) ~ sqrt (dist) vt - variogram(f, meuse) vt.fit - fit.variogram(vt, + vgm(1, Exp, 300, 1)) v - variogram(log (zinc) ~ 1, meuse) v.fit - fit.variogram(v, + vgm(1, Sph, 800, 1)) 6/24 空间预测