“空间”数据的平滑

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

查看@Paul Heimstra 的回答here https://stackoverflow.com/questions/10047870/plotting-interpolated-data-on-map。他们建议从薄板样条插值fields http://cran.r-project.org/web/packages/fields/fields.pdf包裹。如果克里金法适合您，该软件包automap http://cran.r-project.org/web/packages/automap/automap.pdf可能有用。

如果您关心地理值的插值，那么值得一提DG Rossiter 的地统计学和开源统计计算在线/免费课程 http://www.css.cornell.edu/faculty/dgr2/teach/degeostats.html.