在R语言中绘制圣诞树农场图形及圣诞树形状图的方法咨询

一、简易绘制圣诞树农场 #

你可以用R的基础绘图系统快速实现多棵圣诞树的农场效果，核心思路是循环生成随机位置，给每个位置绘制由三角形（树冠）和矩形（树干）组成的基础圣诞树：

# 设置画布
plot(0, 0, xlim = c(0, 100), ylim = c(0, 80), type = "n", axes = FALSE, xlab = "", ylab = "")

# 生成20棵随机位置的圣诞树
set.seed(123) # 固定随机种子保证可复现
for (i in 1:20) {
  x <- runif(1, 10, 90)
  y <- runif(1, 5, 60)
  
  # 绘制三层树冠
  polygon(c(x-5, x, x+5), c(y+10, y+20, y+10), col = "darkgreen")
  polygon(c(x-7, x, x+7), c(y+5, y+15, y+5), col = "darkgreen")
  polygon(c(x-9, x, x+9), c(y, y+10, y), col = "darkgreen")
  
  # 绘制树干
  rect(x-2, y-5, x+2, y, col = "saddlebrown")
}

调整循环次数、位置范围和图形尺寸，就能修改农场的规模和树的大小。

二、生成类似示例的圣诞树形状 #

示例中的圣诞树可以用参数方程或分形算法实现，以下是两种实用方案：

1. 参数方程生成平滑圣诞树 #

用三角函数组合模拟树冠的曲线起伏，搭配线性项控制整体锥型：

library(ggplot2)

# 定义参数方程
t <- seq(0, 2*pi, length.out = 1000)
r <- 5*sin(5*t) + 10*(pi - t)/pi 

# 转换为笛卡尔坐标并调整方向
x <- r * cos(t)
y <- -r * sin(t) + max(-r * sin(t))

# 绘制图形
ggplot() +
  geom_path(aes(x, y), color = "darkgreen", linewidth = 2) +
  geom_rect(aes(xmin = -1, xmax = 1, ymin = 0, ymax = 3), fill = "saddlebrown") +
  theme_void()

2. 分形算法生成细节圣诞树 #

基于科赫曲线的递归变种，通过迭代生成带分支的树冠：

# 递归分形圣诞树函数
koch_tree <- function(x1, y1, x2, y2, depth) {
  if (depth == 0) {
    lines(c(x1, x2), c(y1, y2), col = "darkgreen", lwd = 2)
  } else {
    dx <- x2 - x1
    dy <- y2 - y1
    
    # 拆分线段并计算凸起顶点
    x3 <- x1 + dx/3
    y3 <- y1 + dy/3
    x4 <- x1 + 2*dx/3
    y4 <- y1 + 2*dy/3
    angle <- pi/3
    x5 <- x3 + dx/3*cos(angle) - dy/3*sin(angle)
    y5 <- y3 + dx/3*sin(angle) + dy/3*cos(angle)
    
    # 递归绘制各段
    koch_tree(x1, y1, x3, y3, depth-1)
    koch_tree(x3, y3, x5, y5, depth-1)
    koch_tree(x5, y5, x4, y4, depth-1)
    koch_tree(x4, y4, x2, y2, depth-1)
  }
}

# 设置画布并绘制
plot(0, 0, xlim = c(-20, 20), ylim = c(0, 30), type = "n", axes = FALSE, xlab = "", ylab = "")
koch_tree(-15, 0, 15, 0, depth = 3)
rect(-2, 0, 2, -5, col = "saddlebrown")

调整depth参数可控制分形精细度，数值越大细节越丰富。

内容的提问来源于stack exchange，提问作者psysky