R热图使用插值

Question

这个主题解释了如何使用akima包从x/y/z坐标创建矩阵,但我宁愿不使用新包.毕竟,你可以在gnuplot的一个命令中做同样的事情:gnuplot插值.

使用上述gnuplot线程中的数据点矩阵的图像绘图可以在R中通过调用image()命令来实现.

那么表面图怎么样,我们如何插入x/y/z点来生成2d热图？

Answer 1

好吧,如果你不喜欢akima::interp,也许stats::loess是你的口味？

topo.loess <- loess (z ~ x * y, topo, degree = 2, span = 0.2)
x <- seq (min (topo$x), max (topo$x), .05)
y <- seq (min (topo$y), max (topo$y), .05)
interpolated <- predict (topo.loess, expand.grid (x = x, y = y))
image (x= x, y= y, z = interpolated, asp = 1)
points (topo)

(也沿着MASS的路线)

否则:为什么不使用gnuplot？虽然它可能被认为是一个不同的软件...