用R语言绘制三维曲面

最新推荐文章于 2024-06-20 11:21:44 发布

zoujiahui_2018

最新推荐文章于 2024-06-20 11:21:44 发布

阅读量4.1k

点赞数 1

分类专栏： # R语言中的绘图文章标签： r语言开发语言

本文链接：https://blog.csdn.net/qq_18055167/article/details/130877198

版权

R语言中的绘图专栏收录该内容

21 篇文章

订阅专栏

该文展示了如何使用R语言的不同图形库，包括graphics,lattice和rgl，来创建3D图形。具体包括使用persp进行3D透视图，wireframe创建3D网格图，rgl库的plot3d和surface3d绘制3D点和表面，以及scatterplot3d绘制3D散点图和拟合平面。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

graphics::persp

x <- seq(-10, 10, length.out = 30)
y <- x
f <- function(x, y) { r <- sqrt(x^2+y^2); 10 * sin(r)/r }
z <- outer(x, y, f)
z[is.na(z)] <- 1
op <- par(bg = "white")
persp(x, y, z, theta = 30, phi = 30, expand = 0.5, col = "lightblue")
persp(x, y, z, theta = 30, phi = 30, expand = 0.5, col = "lightblue",
      ltheta = 120, shade = 0.75, ticktype = "detailed",
      xlab = "X", ylab = "Y", zlab = "Sinc( r )"
) -> res
round(res, 3)

在这里插入图片描述

lattice::wireframe

library(lattice) 
# 构建网格数据
d <- expand.grid(x =10:60, y = 5:40,g=2:3)
# 构建因变量
d$z <-log((d$x^d$g+d$y^2)*d$g)
# 绘制曲面拟合图
wireframe(z~x*y,
          data = d,
          shade=T,  # 曲面颜色渲染
          groups=g, #分组
          drape=T,   # 是否以右侧的色柱区分曲面
          colorkey=T,  # 是否显示色标
          scales=list(arrows=T), #是否显示箭头
          screen=list(z=45,x=-60) #调整图片角度
)

在这里插入图片描述

rgl::plot3d

library(rgl)
d <- expand.grid(x =10:60, y = 5:40,g=2:3)
d$z <-log((d$x^d$g+d$y^2)*d$g)
plot3d(d$x,d$y,d$z,
       type="s",
       size = 0.5,
       col=terrain.colors(length(d$z)))

在这里插入图片描述
plot3d()所绘制的图形是可以用鼠标调整视角的。

plot3d(lm(mpg ~ wt + I(wt^2) + qsec, 
          data = mtcars),
       plane.col =topo.colors(1))

在这里插入图片描述

rgl::surface3d

# 设置x，y，z值
z <- 4 * volcano        
x <- 8 * (1:nrow(z))   
y <- 6 * (1:ncol(z))   
# 设置 颜色
zlim <- range(z)
zlen <- zlim[2] - zlim[1] + 1
colorlut <- topo.colors(zlen) 
col <- colorlut[ z - zlim[1] + 1 ] 
# 绘制曲面图
open3d()
surface3d(x, y, z, color = col, back = "lines")

在这里插入图片描述

scatterplot3d::scatterplot3d

在这里插入图片描述

# 绘制3D 散点图
data(trees)
s3d <- scatterplot3d(trees, 
                     type="p", 
                     highlight.3d=TRUE,
                     angle=60, 
                     #scale.y=0.8, 
                     pch=16, 
                     zlim = c(0,80),
                     col.grid="lightblue",
                     main="scatterplot3d")
                     
# 散点图中添加拟合曲面
tree_lm <- lm(Volume ~ Girth + Height,data = trees)
s3d$plane3d(tree_lm, lty.box = "solid")

在这里插入图片描述