R语言与GIS：[3]插值与地统计(3)

使用gstat程序包(2)

克里金插值

> lzn.kriged=krige(log(zinc)~1,meuse,meuse.grid,model=lzn.fit)

[using ordinary kriging]

利用krige进行普通克里金插值计算。

> spplot(lzn.kriged["var1.pred"])

利用spplot函数绘制插值的估计结果，如下图所示。

条件模拟

> lzn.condsim=krige(log(zinc)~1,meuse,meuse.grid,model=lzn.fit,nmax=30,nsim=4)

drawing 4 GLS realisations of beta...

[using conditional Gaussian simulation]

利用条件高斯模拟方法完成4个实现。

> spplot(lzn.condsim,main="三个条件模拟")

绘制模拟结果，如下图所示。