[英]Increasing 3D surface plot thickness using the Plot3D package in R
I was able to get this to work by creating multiple surfaces and connecting them using the curtain
argument of the plot3D function.我能够通过创建多个表面并使用 plot3D 函数的curtain
参数连接它们来实现这一点。 The thickness was controlled by the z
argument of each surface.厚度由每个表面的z
参数控制。 In case anyone needs the sample code, here it is;如果有人需要示例代码,这里是;
For simplicity, I used meters to specify thicknesses/height of the surfaces为简单起见,我使用米来指定表面的厚度/高度
library(plot3D)
library(plot3Drgl)
# volcano data range
clim <- range(volcano)
# Surface 1
persp3D(z = volcano, zlim = c(0, 600), clim = clim,
box = FALSE, plot = FALSE, curtain = TRUE, border = "black")
# Surface 2 - 20m below surface 1
persp3D(z = volcano - 20, clim = clim, colvar = volcano,
add = TRUE, colkey = FALSE, plot = FALSE, curtain = TRUE, border = "black")
# Surface 3 - 40m below surface 1
persp3D(z = volcano - 40, clim = clim, colvar = volcano,
add = TRUE, colkey = FALSE, plot = FALSE, curtain = TRUE, border = "black")
# Surface 4 - 60m below surface 1
persp3D(z = volcano - 60, clim = clim, colvar = volcano,
add = TRUE, colkey = FALSE, plot = TRUE, curtain = TRUE, border = "black")
# Surface 5 - 80m below surface 1
persp3D(z = volcano - 80, clim = clim, colvar = volcano,
add = TRUE, colkey = FALSE, plot = TRUE, curtain = FALSE, border = "black")
# Open rgl for interactive viewing
plotrgl()
I don't know about plot3D
, but rgl
doesn't have a specific function for that.我不知道plot3D
,但rgl
没有特定的功能。 The way to do it is to calculate the polygons corresponding to the edges of the slab, and use polygon3d
to plot those.这样做的方法是计算与平板边缘相对应的多边形,并使用polygon3d
绘制这些多边形。 For example:例如:
data(volcano)
persp3d(volcano, col = "green")
# The bottom of the slab
persp3d(volcano-10, col = "blue", add = TRUE)
minx <- 0
maxx <- 1
miny <- 0
maxy <- 1
m <- nrow(volcano)
n <- ncol(volcano)
# The front edge
edgex <- c(seq(minx, maxx, length.out = m),
seq(maxx, minx, length.out = m))
edgey <- miny
edgez <- c(volcano[,1],rev(volcano[,1] - 10))
polygon3d(cbind(edgex, edgey, edgez), coords = c(1,3),
col = "yellow")
# The back edge
edgey <- maxy
edgez <- c(volcano[,n],rev(volcano[,n] - 10))
polygon3d(cbind(edgex, edgey, edgez), coords = c(1,3), col = "white")
edgex <- minx
edgey <- c(seq(miny, maxy, length.out = n),
seq(maxy, miny, length.out = n))
edgez <- c(volcano[1,],rev(volcano[1,] - 10))
polygon3d(cbind(edgex, edgey, edgez), coords = c(2,3), col = "black")
edgex <- maxx
edgez <- c(volcano[m,],rev(volcano[m,] - 10))
polygon3d(cbind(edgex, edgey, edgez), coords = c(2,3), col = "green")
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.