R - 补丁中的3d图

Question

我在数据框中有以下数据：

**x** in (0,1)
**y** in [0,1]
**z** in [0,1]

例如：

X,Y,Z
0.1, 0.2, 0.56
0.1, 0.3, 0.57
...

我想在这种类型的图表上绘制它们：

我试过R，但我能得到的只是not-so-fancy 3d scatterplot。 我也读到了lattice 3d wireframe，但我无法理解它。

我应该怎样做才能在R中获得类似线框的Matlab？ 涉及哪些数据转换？

这是文档中的示例代码：

x <- seq(-pi, pi, len = 20)
y <- seq(-pi, pi, len = 20)
g <- expand.grid(x = x, y = y)
g$z <- sin(sqrt(g$x^2 + g$y^2))
wireframe(z ~ x * y, g, drape = TRUE,
aspect = c(3,1), colorkey = TRUE)

我不觉得特别清楚。

编辑：persp3d功能正常，我可以使用一种颜色生成三维图。如何设置相对于z值的色阶？

感谢任何提示， Mulone

Answer 1

使用outer创建z值，然后使用persp绘制：

z <- outer(x,y, function(x,y) sin(sqrt(x^2+y^2)))
persp(x,y,z)

有着色和设置视角的选项，请参阅?persp。请参阅Matlab样式着色的第四个示例。

对于互动情节，请考虑在persp3d包中使用rgl：

require(rgl)
persp3d(x,y,z,col="blue")

修改

要添加颜色，与persp中的方法略有不同，因为颜色与顶点有关，而不是与构面的中心有关，但它更容易。

jet.colors <- colorRampPalette( c("blue", "green") ) 
pal <- jet.colors(100)
col.ind <- cut(z,100) # colour indices of each point
persp3d(x,y,z,col=pal[col.ind])

帮助文件建议添加参数smooth=FALSE，但这取决于个人偏好。