[英]Increasing 3D surface plot thickness using the Plot3D package in R
我能够通过创建多个表面并使用 plot3D 函数的curtain
参数连接它们来实现这一点。 厚度由每个表面的z
参数控制。 如果有人需要示例代码,这里是;
为简单起见,我使用米来指定表面的厚度/高度
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()
我不知道plot3D
,但rgl
没有特定的功能。 这样做的方法是计算与平板边缘相对应的多边形,并使用polygon3d
绘制这些多边形。 例如:
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.