如何在ggplot2中利用现成函数轻松绘制椭圆切线？

目前没有专门的geom_tangent类开箱即用的ggplot2几何对象函数，但有两个简便方案，完全不用手动推导切线方程就能实现你的需求：

方案1：使用conicsections包（推荐） #

这个包专门处理圆锥曲线计算，内置了从外部点到椭圆的切线求解功能，只需调用现成方法即可：

示例代码 #

# 首次使用先安装包
# install.packages("conicsections")

library(ggplot2)
library(ggforce)
library(conicsections)

# 定义目标椭圆的参数
target_ellipse <- Ellipse$new(h = 30, k = 6, a = 2, b = 20, angle = 1.5)

# 计算从原点(0,0)出发的切线
origin_tangents <- target_ellipse$get_tangents_from_point(x = 0, y = 0)

# 转换为ggplot可直接使用的线段数据格式
tangent_segments <- do.call(rbind, lapply(origin_tangents, function(t) {
  data.frame(
    x_start = 0, y_start = 0,
    x_end = t$point_on_ellipse$x, y_end = t$point_on_ellipse$y
  )
}))

# 绘制椭圆+切线
ggplot() + 
  geom_ellipse(aes(x0 = 30, y0 = 6, a = 2, b = 20, angle = 1.5)) + 
  geom_segment(data = tangent_segments, aes(x = x_start, y = y_start, xend = x_end, yend = y_end), 
               color = "red", linewidth = 1) +
  scale_x_continuous(limits = c(-20, 55)) +
  scale_y_continuous(limits = c(-30, 30))

如果要指定其他起点，只需修改get_tangents_from_point(x = ?, y = ?)中的坐标，函数会自动计算对应的切线。

方案2：借助ggforce底层工具封装辅助函数 #

ggforce内部有椭圆坐标转换工具，基于这些工具封装函数，同样不用手动推导公式：

示例代码 #

library(ggplot2)
library(ggforce)
library(dplyr)

# 封装函数：输入椭圆参数和起点，返回切线线段数据
get_ellipse_tangents <- function(x0, y0, a, b, angle, start_x, start_y) {
  # 将起点转换到椭圆的局部坐标系
  trans <- ellipse_shape_trans(a = a, b = b, angle = angle)
  local_start <- trans$transform(data.frame(x = start_x - x0, y = start_y - y0))
  
  # 利用椭圆转圆的性质计算切线点（避免手动推导椭圆切线公式）
  r <- max(a, b)
  d <- sqrt(local_start$x^2 + local_start$y^2)
  if (d <= r) stop("起点在椭圆内部，不存在切线")
  theta <- atan2(local_start$y, local_start$x)
  alpha <- asin(r/d)
  tangent_angles <- c(theta - alpha, theta + alpha)
  
  # 将切线点转换回全局坐标系
  local_points <- data.frame(x = r*cos(tangent_angles), y = r*sin(tangent_angles))
  global_points <- trans$inverse_transform(local_points) %>% mutate(x = x + x0, y = y + y0)
  
  # 生成ggplot可用的线段数据
  data.frame(
    x_start = start_x, y_start = start_y,
    x_end = global_points$x, y_end = global_points$y
  )
}

# 生成从原点出发的切线数据
tangent_data <- get_ellipse_tangents(x0 = 30, y0 = 6, a = 2, b = 20, angle = 1.5, start_x = 0, start_y = 0)

# 绘图
ggplot() + 
  geom_ellipse(aes(x0 = 30, y0 = 6, a = 2, b = 20, angle = 1.5)) + 
  geom_segment(data = tangent_data, aes(x = x_start, y = y_start, xend = x_end, yend = y_end), 
               color = "blue", linewidth = 1) +
  scale_x_continuous(limits = c(-20, 55)) +
  scale_y_continuous(limits = c(-30, 30))

这个方案封装后可重复使用，同样支持任意指定起点。

内容的提问来源于stack exchange，提问作者J. Doe