使用geom_rect在加拿大纽芬兰地图上绘制旋转高亮框标注三一湾

实现斜向标注框的可行方案 #

geom_rect仅支持绘制与坐标轴对齐的矩形，无原生旋转参数，可通过以下两种方法实现需求：

方案1：自定义旋转多边形（推荐，实现成本最低） #

先计算矩形旋转后的四个顶点坐标，再用geom_polygon绘制即可，示例代码如下：

library(ggplot2)
library(maptools) # 配合getData读取GADM数据

# 定义坐标旋转函数：输入原始坐标点、旋转中心、旋转角度（角度制）
rotate_pt <- function(x, y, center_x, center_y, angle = -25) {
  rad <- angle * pi / 180
  x_rel <- x - center_x
  y_rel <- y - center_y
  new_x <- x_rel * cos(rad) - y_rel * sin(rad) + center_x
  new_y <- x_rel * sin(rad) + y_rel * cos(rad) + center_y
  return(data.frame(x = new_x, y = new_y))
}

# 获取底图数据
Provinces <- getData("GADM", country="CAN", level=1)

# 生成原始矩形的四个顶点，按顺序排列
rect_pts <- data.frame(
  x = c(-53.75, -53.75, -53, -53),
  y = c(47.5, 48.5, 48.5, 47.5)
)
# 旋转顶点，中心为原矩形中心，可调整angle参数匹配湾区走向
rotated_pts <- rotate_pt(rect_pts$x, rect_pts$y, 
                         center_x = mean(c(-53.75, -53)), 
                         center_y = mean(c(47.5, 48.5)),
                         angle = -25)

# 绘图
ggplot() +
  geom_polygon(data=Provinces, aes(long, lat, group=group), colour="grey10", 
               fill="cadetblue1", alpha=0.8) +
  coord_map(xlim=c(-60, -52), ylim=c(46, 52)) +
  theme_bw() +
  # 绘制旋转后的矩形框
  geom_polygon(data = rotated_pts, aes(x = x, y = y), 
               alpha = 0, colour = "red", size = 1, linetype = 1)

你可以调整angle参数的数值、原始矩形的长宽范围，直到框体刚好覆盖三一湾区域。

方案2：矢量边界高亮（更精准） #

如果需要完全贴合湾区轮廓标注，可以直接获取三一湾的官方矢量边界数据，导入后用geom_polygon或者geom_sf绘制高亮边框/填充，无需手动调整旋转角度。

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