如何在3D图形中加一个圆或椭圆平面

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

谢谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：3D图 平面 如何

1. # x and y follow independent normal
   
2. sd.x <- 1
   
3. sd.y <- 2
   
4. x <- c(-50:50/10)*sd.x
   
5. y <- c(-50:50/10)*sd.y
   
6. z <- outer(dnorm(x, sd = sd.x), dnorm(y, sd = sd.y))
   
7. # Target z value (half of the maxium)
   
8. z.tgl <- max(z)/2
   
9. # Inverse of dnorm (x > 0 part)
   
10. inv.dnorm <-
    
11.   function(d, mean = 0, sd = 1)
    
12.   {
    
13.     return(mean + sqrt(-2*sd*sd*log(sqrt(2*pi)*sd*d)))
    
14.   }
    
15. # Half range
    
16. hlfrng.x <-
    
17.   seq(0, inv.dnorm(z.tgl/max(dnorm(y, sd = sd.y)), sd = sd.x), length.out = 51)
    
18. hlfrng.y <-
    
19.   inv.dnorm(z.tgl/dnorm(hlfrng.x, sd = sd.x), sd = sd.y)
    
20. # Target x, y value
    
21. x.tgl <-
    
22.   c(hlfrng.x, rev(head(hlfrng.x, -1)), -1*head(hlfrng.x[-1], -1), rev(-1*hlfrng.x[-1]))
    
23. y.tgl <-
    
24.   c(hlfrng.y, rev(-1*hlfrng.y)[-1], -1*head(hlfrng.y[-1], -1), head(rev(hlfrng.y), -1))
    
25. # Assert correct length
    
26. stopifnot(length(x.tgl) == length(y.tgl))
    
27. # Normal surface
    
28. persp(x, y, z, theta = 30, scale = T, box = F, axes = F, zlim = c(0, 0.15), border = "green4") -> res
    
29. # Eclipse
    
30. polygon(trans3d(x.tgl, y.tgl, rep(z.tgl, times = length(x.tgl)), res), border = "Navy", col = "#00007F3F")

复制代码

Result: