插值和空间分析(R语言)
> library(lattice)
> library(sp)
> data(meuse)
> coordinates(meuse) <- c("x","y")
> spplot(meuse, "zinc", do.log=T)
> bubble(meuse, "zinc", do.log=T, key.space="bottom")
> xyplot(log(meuse$zinc)~sqrt(meuse$dist), meuse, main="", xlab="dist", ylab="Zn")
> meuse$fitted.s <- predict(zn.lm, meuse) - mean(predict(zn.lm,meuse)) > meuse$residuals <- residuals(zn.lm) > spplot(meuse, c("fitted.s", "residuals")) > spplot(meuse, c("fitted.s", "residuals"))
> library(gstat) > idw.out <- idw(zinc~1, meuse, meuse.grid, idp=1) [inverse distance weighted interpolation] > spplot(idw.out)
> spplot(idw.out, c("var1.pred"))
郑重声明:本站内容如果来自互联网及其他传播媒体,其版权均属原媒体及文章作者所有。转载目的在于传递更多信息及用于网络分享,并不代表本站赞同其观点和对其真实性负责,也不构成任何其他建议。