如何平滑ggplot2绘制的圆弧？解决路径锯齿问题

问题描述 #

我尝试通过三个点绘制一段圆弧，但得到的圆弧图像存在明显的锯齿状不平滑问题。以下是我的实现代码：

library(ggplot2)
plot_arc <- function(C, r, P1, P2, P3) {
  t <- seq(atan2(P1[2] - C[2], P1[1] - C[1]), atan2(P2[2] - C[2], P2[1] - C[1]), length.out = 100000)
  x <- C[1] + r * cos(t)
  y <- C[2] + r * sin(t)
  data <- data.frame(x = x, y = y)
  
  p <- ggplot() +
    geom_path(data = data, aes(x = x, y = y), color = "blue", size=1.) +
    geom_point(aes(x = P1[1], y = P1[2]), color = "red", size = 3) +
    geom_point(aes(x = P2[1], y = P2[2]), color = "green", size = 3) +
    geom_point(aes(x = P3[1], y = P3[2]), color = "yellow", size = 3) +
    coord_fixed() 
 
  print(p) 
}

plot_arc(C_optimal, r_optimal, P1, P2, P3)

问题分析 #

• 角度序列逻辑缺陷：函数参数包含P3但未实际使用，生成的圆弧仅连接P1和P2，无法保证经过P3；同时若P1与P2的角度跨-π/π边界，seq()会生成绕圆一周的长弧，而非预期的短弧，导致路径异常。
• 过度采样+缺乏抗锯齿：设置100000个采样点属于性能浪费，且未开启渲染层面的抗锯齿优化，加剧了锯齿感。

解决步骤 #

1. 修正圆弧生成逻辑（确保经过三个点） #

重新计算角度范围，生成经过P3的平滑短弧，同时合理控制采样点数量：

library(ggplot2)
plot_arc <- function(C, r, P1, P2, P3) {
  # 计算三个点相对圆心的极角
  theta1 <- atan2(P1[2] - C[2], P1[1] - C[1])
  theta2 <- atan2(P2[2] - C[2], P2[1] - C[1])
  theta3 <- atan2(P3[2] - C[2], P3[1] - C[1])
  
  # 处理角度跨-π/π边界的情况，确定短弧的起止角度
  diff12 <- (theta2 - theta1) %% (2 * pi)
  # 判断是否需要反转角度方向以生成短弧
  if (diff12 > pi) {
    theta_start <- theta2
    theta_end <- theta1
  } else {
    theta_start <- theta1
    theta_end <- theta2
  }
  
  # 生成1000个采样点（足够平滑，避免性能浪费）
  t <- seq(theta_start, theta_end, length.out = 1000)
  # 若角度跨边界，补充完整短弧的采样点
  if ((theta_end - theta_start) %% (2 * pi) > pi) {
    t <- c(t, seq(theta_end + 2 * pi, theta_start, length.out = 1000)[-1])
  }
  
  # 计算圆弧坐标
  x <- C[1] + r * cos(t)
  y <- C[2] + r * sin(t)
  data <- data.frame(x = x, y = y)
  
  p <- ggplot() +
    # 添加抗锯齿参数优化边缘平滑度
    geom_path(data = data, aes(x = x, y = y), color = "blue", size = 1.,
              lineend = "round", linejoin = "round") +
    geom_point(aes(x = P1[1], y = P1[2]), color = "red", size = 3) +
    geom_point(aes(x = P2[1], y = P2[2]), color = "green", size = 3) +
    geom_point(aes(x = P3[1], y = P3[2]), color = "yellow", size = 3) +
    coord_fixed() 
 
  print(p) 
}

# 调用函数（需提前定义C_optimal, r_optimal, P1, P2, P3）
plot_arc(C_optimal, r_optimal, P1, P2, P3)

2. 渲染与输出优化 #

• 启用抗锯齿参数：在geom_path中添加lineend="round"和linejoin="round"，直接优化线条边缘的平滑度。
• 高分辨率保存：使用ggsave保存图像时，指定高DPI并启用Cairo渲染（系统级抗锯齿）：

ggsave("smooth_arc.png", plot = p, dpi = 300, type = "cairo-png")

• RStudio预览优化：打开Tools -> Global Options -> Graphics，将Anti-aliasing设置为2x或4x，提升实时绘图的平滑度。

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