r - "spatial"数据的平滑处理

Question

我有 2 个变量 x 和 y，它们是 [0,1] 处的笛卡尔坐标，z 是 a (连续)变量在这些坐标处。 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 或lattice 可以完成我的工作那就没问题了。任何其他更好的可视化想法/建议也将受到欢迎。

Answer 1

查看 @Paul Heimstra 的回答 here 。他们建议从 fields 进行薄板样条插值包裹。如果克里金法适合您，软件包 automap可能有用。

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