圣诞节到了，分享别人种的圣诞树

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

1. part <- list(x0=0,y0=0,x1=0,y1=1,
   
2.     branch1=NULL,branch2=NULL,extend=NULL,
   
3.     lwd=1,depth=0,col='springgreen')
   
4. 
   
5. par(mfrow=c(1,1),mar=c(5, 4, 4, 2) + 0.1)
   
6. segplot <- function(tree) {
   
7.   if (is.null(tree)) return()
   
8.   segments(tree$x0,tree$y0,tree$x1,tree$y1,
   
9.       col=tree$col,
   
10.       lwd=tree$lwd)
    
11.   segplot(tree$branch1)
    
12.   segplot(tree$branch2)
    
13.   segplot(tree$extend)
    
14. }
    
15. #segplot(part)
    
16. 
    
17. grow <- function(tree) {
    
18.   if (is.null(tree) ) return(NULL)
    
19.   
    
20.   tree$lwd=tree$lwd*1.2
    
21.   
    
22.   if (tree$lwd>2.5) tree$col <- 'brown'
    
23.   if (is.null(tree$extend)) {
    
24.     tree$extend <- list(
    
25.         x0=tree$x1,
    
26.         y0=tree$y1,
    
27.         x1=rnorm(1,1,.03)*(2*tree$x1-tree$x0),
    
28.         y1=(rnorm(1,.98,.02)+.02*(tree$x1==tree$x0))*(2*tree$y1-tree$y0),
    
29.         branch1=NULL,
    
30.         branch2=NULL,
    
31.         extend=NULL,
    
32.         lwd=1,
    
33.         depth=tree$depth,
    
34.         col=tree$col
    
35.     )
    
36.     length=sqrt((tree$x1-tree$x0)^2 + (tree$y1-tree$y0)^2)
    
37.     angle <- asin((tree$x1-tree$x0)/length)
    
38.     branch <- list(
    
39.         x0=(tree$x1+tree$x0)/2,
    
40.         y0=(tree$y1+tree$y0)/2,
    
41.         branch1=NULL,
    
42.         branch2=NULL,
    
43.         extend=NULL,
    
44.         lwd=1,
    
45.         depth=tree$depth,
    
46.         col=tree$col
    
47.     )
    
48.     shift <- rnorm(2,.5,.1)
    
49.     branch$x0 <- shift[1]*tree$x1+(1-shift[1])*tree$x0
    
50.     branch$y0 <- shift[1]*tree$y1+(1-shift[1])*tree$y0
    
51.     length=length*rnorm(1,.5,.05)
    
52.     co <- runif(1,.35,.45)
    
53.     branch$x1 <- branch$x0+sin(angle+co)*length
    
54.     branch$y1 <- branch$y0+cos(angle+co)*length
    
55.     tree$branch1 <- branch
    
56.     branch$x0 <- shift[2]*tree$x1+(1-shift[2])*tree$x0
    
57.     branch$y0 <- shift[2]*tree$y1+(1-shift[2])*tree$y0
    
58.     co <- runif(1,.35,.45)
    
59.     branch$x1 <- branch$x0+sin(angle-co)*length
    
60.     branch$y1 <- branch$y0+cos(angle-co)*length
    
61.     tree$branch2 <- branch   
    
62.   } else {
    
63.     tree$branch1 <- grow(tree$branch1)
    
64.     tree$branch2 <- grow(tree$branch2)
    
65.     tree$extend <- grow(tree$extend)
    
66.   }
    
67.   tree$depth <- tree$depth+1
    
68.   if (tree$depth>2)  tree$col <- 'green'
    
69.   if (tree$depth>4)  tree$col <- 'darkgreen'
    
70.   if (tree$depth>6)  tree$col <- 'brown'
    
71.   
    
72.   tree
    
73. }
    
74. tree <- part
    
75. for (i in 1:9) tree <- grow(tree)
    
76. par(mar=c(0,0,0,0))
    
77. plot(x=c(-3,3),y=c(0,9),type='n',axes=FALSE,xlab='',ylab='')
    
78. segplot(tree)

复制代码

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：圣诞节 圣诞树 PART ART 圣诞树 圣诞节

太犀利了 赞一个 学习了

狂赞！！！！！！！！   

这个圣诞树不怎么性感丰满， 比较中国的秀气。
                     

上面再挂些礼物就好了

收藏了，谢谢分享

nuomin 发表于 2013-12-24 11:38
收藏了，谢谢分享

看见你回复了   以为你会画个更好的瞧瞧的。。。

太棒了

哈哈，感谢分享，牛人