平滑"空间"数据

Question

我有2个变量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或格子可以完成我的工作,那就没问题了.任何其他更好的可视化的想法/建议也将受到欢迎.

Answer 1

在这里查看@Paul Heimstra的答案.他们建议从字段包中进行薄板样条插值.如果kriging适合您,则包自动化可能很有用.