使用R绘制具有等高线图叠加的3D表面图

时间:2009-12-13 12:53:44

标签: r 3d statistics plot

我有一个3元组的数据集(X,Y,Z点)我想用R绘图。

我想从数据中创建一个曲面图,并在曲面图上叠加一个等高线图,以便创建轮廓图的印象,即从表面图中“阴影”或投影。等高线图将显示在曲面图下方。

我的数据集看起来有点像这样:

Axis  |  Data Type
-------------------
X     |  Date value
Y     |  Float value
Z     |  Float value

我怎样才能做到这一点?

1 个答案:

答案 0 :(得分:85)

编辑:

我刚看到你指出你的一个尺寸是日期。在这种情况下,have a look at Jeff Ryan's chartSeries3d用于绘制三维时间序列图。他在这里显示了随时间变化的收益率曲线:

chartSeries example

原始答案:

据我了解,您希望将计数地图作为3D曲面图下方平面上的投影。我不相信除了创建两个图然后组合它们之外,还有一种简单的方法可以做到这一点。你可以find the spatial view helpful for this

3D绘图有两个主要的R套餐: rgl (或者您可以使用相关的misc3d套餐)和 scatterplot3d

<强> RGL

rgl包使用OpenGL创建交互式3D图(read more on the rgl website)。以下是使用surface3d函数的示例:

library(rgl)
data(volcano)
z <- 2 * volcano # Exaggerate the relief
x <- 10 * (1:nrow(z)) # 10 meter spacing (S to N)
y <- 10 * (1:ncol(z)) # 10 meter spacing (E to W)
zlim <- range(z)
zlen <- zlim[2] - zlim[1] + 1
colorlut <- terrain.colors(zlen,alpha=0) # height color lookup table
col <- colorlut[ z-zlim[1]+1 ] # assign colors to heights for each point
open3d()
rgl.surface(x, y, z, color=col, alpha=0.75, back="lines")

alpha参数使此表面部分透明。现在,您有一个表面的交互式3D绘图,并且您想在下面创建一个计数地图。 rgl允许您向现有图像添加更多绘图:

colorlut <- heat.colors(zlen,alpha=1) # use different colors for the contour map
col <- colorlut[ z-zlim[1]+1 ] 
rgl.surface(x, y, matrix(1, nrow(z), ncol(z)),color=col, back="fill")

在这个表面中,我设置高度= 1,以便我们在另一个表面下面有一个平面。最终看起来像这样,可以用鼠标旋转:

3D surface plot

<强> scatterplot3d

scatterplot3d与R(read the vignette)中的其他绘图函数更像。这是一个简单的例子:

temp <- seq(-pi, 0, length = 50)
x <- c(rep(1, 50) %*% t(cos(temp)))
y <- c(cos(temp) %*% t(sin(temp)))
z <- c(sin(temp) %*% t(sin(temp)))
scatterplot3d(x, y, z, highlight.3d=TRUE,
 col.axis="blue", col.grid="lightblue",
 main="scatterplot3d - 2", pch=20)

在这种情况下,您需要覆盖图像。 R-Wiki has a nice post on creating a tanslucent background image