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='三个条件模拟')绘制模拟结果，如下图所示。