在ggplot(R)中如何基于距离调整点/图形的填充透明度？

问题描述 #

我正在处理地图相关工作，希望以连续方式呈现某事物的影响力随距中心距离增加而递减的效果（具体递减速率无关紧要），类似图像编辑软件中画笔的「硬度」设置，需要硬度介于0%-100%之间的圆形效果。目前我使用sf包绘制基础地图，代码形式如下：

ggplot() +
  annotation_spatial(mapdata) +
  layer_spatial(mapdata)

想知道能否通过annotate("point",...)或geom_circle()调整填充alpha值来实现？

解决方案 #

直接用annotate("point")或geom_circle()无法实现连续的alpha衰减效果——这两个函数的alpha参数是全局统一值，没法随距离中心的远近动态变化。要实现你要的画笔硬度式渐变效果，推荐以下几种方法：

方法1：嵌套同心圆模拟渐变 #

生成一组半径递增的同心圆，让每个圆的alpha值随半径增大逐步降低，叠加后就能模拟连续衰减的视觉效果。如果需要设置「硬度」，可以让前N个小半径的圆保持alpha=1（对应画笔的硬芯部分），再让外围圆的alpha递减。

示例代码：

library(sf)
library(ggplot2)

# 定义中心坐标（替换成你的实际经纬度/投影坐标）
center <- st_sfc(st_point(c(116.397, 39.908)), crs = 4326)
# 匹配地图数据的坐标系
center <- st_transform(center, st_crs(mapdata))

# 生成同心圆数据：硬度设为50%，即前半段半径保持不透明，后半段衰减
total_radius <- 10  # 最大影响半径，按需调整
hardness <- 0.5  # 0-1之间，对应0%-100%硬度
hard_radii <- seq(0, total_radius * hardness, length.out = 10)
soft_radii <- seq(total_radius * hardness, total_radius, length.out = 10)

# 硬芯部分：alpha=1
hard_circles <- lapply(hard_radii, function(r) {
  st_buffer(center, dist = r) %>% mutate(alpha = 1)
}) %>% do.call(rbind, .)

# 衰减部分：alpha从1线性降到0
soft_circles <- lapply(soft_radii, function(r) {
  st_buffer(center, dist = r) %>% 
    mutate(alpha = 1 - (r - total_radius * hardness)/(total_radius * (1 - hardness)))
}) %>% do.call(rbind, .)

all_circles <- rbind(hard_circles, soft_circles)

# 绘制地图+渐变圆形
ggplot() +
  annotation_spatial(mapdata) +
  layer_spatial(mapdata) +
  geom_sf(data = all_circles, fill = "blue", color = NA, alpha = all_circles$alpha)

方法2：用ggfx包实现径向渐变 #

ggfx包提供了更简洁的径向渐变滤镜，无需手动生成同心圆，直接在指定区域生成从中心到边缘的透明度衰减效果。

示例代码：

library(ggplot2)
library(sf)
library(ggfx)

# 定义中心坐标（需转换为地图的投影坐标）
center_coords <- st_coordinates(st_transform(center, st_crs(mapdata)))
center_x <- center_coords[1]
center_y <- center_coords[2]
radius <- 10

ggplot() +
  annotation_spatial(mapdata) +
  layer_spatial(mapdata) +
  # 绘制径向渐变矩形，裁剪为圆形效果
  annotate(
    "rect",
    xmin = center_x - radius, xmax = center_x + radius,
    ymin = center_y - radius, ymax = center_y + radius,
    fill = radial_gradient(
      colours = c("blue", "transparent"),
      center = c(0.5, 0.5),  # 渐变中心相对于矩形的位置
      radius = c(0.5 * hardness, 0.5)  # 对应硬度：硬芯半径占比，总半径占比
    )
  ) +
  # 用圆形裁剪掉矩形的边角
  geom_circle(
    aes(x0 = center_x, y0 = center_y, r = radius),
    fill = NA, color = NA,
    clip = "on"
  )

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