&Quot；空间数据的平滑 [英] Smoothing of "spatial" data

问题描述

我有两个变量x和y，它们是在[0，1]处的笛卡尔坐标，并且z是在这些坐标处的(连续)变量的值。z向量有一些重要的异常值

x<-sample(seq(0,1,0.001), replace=F)
y<-sample(seq(0,1,0.001), replace=F)
z<-runif(1001,min=0,max=1)
z[100]<-8;z[400]<-16;z[800]<-4

在filled.contour

中显示这些数据时，我要强调这些异常值

到目前为止我一直在使用

library(akima)
a<-interp(x,y,z)
filled.contour(a$x,a$y,a$z)

但我对这种线性插值法不满意。例如(异常值不能正确显示)。

我在想，我需要的是z的某种最近邻居"空间"平滑(基于x，y位置)。有没有人能帮我或指出可以帮助我的数据/示例/包/代码？我更喜欢基本的R解决方案，但如果ggplot2或晶格可以做我的工作，那就好了。任何其他更好地可视化的想法/建议也将受到欢迎。

推荐答案

查看@Paul Heimstra的答案here。他们建议使用fields程序包中的薄板样条插值法。如果克里金法适合您，automap包可能会很有用。

如果您关心的是地理值的内插，DG Rossiter's online/free course on Geostatistics and Open-source statistical computing值得一提。

这篇关于&Quot；空间数据的平滑的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！