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，但我能得到的只是一个不那么花哨的3D散点图 。 我还读到了格子3D线框 ，但我无法理解它。

我应该怎样做才能在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函数工作正常，我能够生成一种颜色的3d图。 如何设置相对于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样式着色的第四个示例。

对于交互式绘图，请考虑在rgl包中使用persp3d ：

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 ，但这smooth=FALSE个人偏好。