ggplot2 scale_fill_stepsn未按预期识别自定义配色方案问题

基于研究主题的出版物数量绘制全球国家热力图，地图主体已完成，但色条效果不符合预期：希望0值区间显示灰色，其余数值区间呈现绿色渐变。该配色逻辑在其他ggplot2场景中可正常生效，但在scale_fill_stepsn及相关函数中无法实现「灰色区间+绿色渐变区间」的目标色条。

所用基础数据代码：

# Load the necessary libraries
library(ggplot2)
library(sf)
library(dplyr)
library(rnaturalearth)
library(rnaturalearthdata)

# Load world map data (countries and continents)
world <- ne_countries(scale = "medium", returnclass = "sf")

注：需额外为world数据框添加Publications列，存储各国出版物数量。

原实现代码：

# Define the color palette for the bins (light grey to dark green)
# Create color palette
green_palette <- colorRampPalette(c("#ADBC9F", "#12372A"))(6)
custom_palette <- c("#B4B4B8", green_palette)
first_color <- "#B4B4B8"  # Gray color for 0 values

# Create a world map with countries color-coded by number of publications using binned colors
map_countries <- ggplot(world_data) +
  geom_sf(aes(fill = Publications), color = "white") +  # Color by publications
  scale_fill_stepsn(
    breaks = c(0, 1, 5, 10, 25, 50, 75, Inf),  # Define bin breaks
    colors = c(first_color, green_palette),  # First color is gray, rest are green gradient
    na.value = first_color  # Handle NA values with the same gray color
  ) +
  
  # Adjust the legend appearance
  guides(
    fill = guide_coloursteps(
      even.steps = TRUE,
      show.limits = TRUE,
      barwidth = unit(15, 'cm'),  # Adjust the width of the colorbar
      title = "Publications"  # Optional: Add title to colorbar
    )
  ) +
  
  theme_minimal() +
  
  theme(
    legend.position = "bottom",
    legend.key.height = unit(1, "cm"),
    legend.title = element_text(face = "bold"),
    legend.spacing = unit(0.25, "cm"),
    legend.justification = "center",
    plot.title = element_text(hjust = 0.5)
  ) +
  
  labs(
    title = "World Map: Publications by Country (Binned Colorbar)",
    subtitle = "Total number of publications per country"
  ) +
  
  # Change projection
  coord_sf(crs = "+proj=robin")  # Use the Robinson projection

# Print country map
print(map_countries)

问题根源 #

1. guide_coloursteps(even.steps = TRUE)会强制所有色阶宽度一致，忽略你设置的不等距断点，导致灰色区间（0-1）和绿色渐变区间的显示比例完全不符合预期。
2. 虽颜色数量与区间数匹配，但未通过配置强化颜色与区间的绑定逻辑，易出现映射偏差。

修改后的代码 #

# 加载依赖库
library(ggplot2)
library(sf)
library(dplyr)
library(rnaturalearth)
library(rnaturalearthdata)

# 加载地图数据并添加出版物数据（替换为你的真实数据）
world <- ne_countries(scale = "medium", returnclass = "sf")
# 示例数据：随机生成出版物数量，实际使用时替换为真实数据
set.seed(123)
world$Publications <- sample(c(0, 1:100), nrow(world), replace = TRUE)
world_data <- world

# 定义配色：灰色（0区间）+6阶绿色渐变
gray_color <- "#B4B4B8"
green_palette <- colorRampPalette(c("#ADBC9F", "#12372A"))(6)
custom_colors <- c(gray_color, green_palette)

# 绘制热力图
map_countries <- ggplot(world_data) +
  geom_sf(aes(fill = Publications), color = "white") +
  scale_fill_stepsn(
    breaks = c(0, 1, 5, 10, 25, 50, 75, Inf),  # 8个断点对应7个区间
    colors = custom_colors,
    limits = c(0, Inf),  # 明确数值范围
    na.value = gray_color,
    # 自定义区间标签，提升可读性
    labels = c("0", "1-4", "5-9", "10-24", "25-49", "50-74", "75+")
  ) +
  guides(
    fill = guide_coloursteps(
      even.steps = FALSE,  # 关闭等宽色阶，按实际区间宽度显示
      show.limits = TRUE,
      barwidth = unit(15, 'cm'),
      title = "出版物数量",
      title.position = "top"
    )
  ) +
  theme_minimal() +
  theme(
    legend.position = "bottom",
    legend.key.height = unit(1, "cm"),
    legend.title = element_text(face = "bold"),
    legend.spacing = unit(0.25, "cm"),
    legend.justification = "center",
    plot.title = element_text(hjust = 0.5)
  ) +
  labs(
    title = "全球各国研究主题出版物数量热力图",
    subtitle = "按国家统计的出版物总数"
  ) +
  coord_sf(crs = "+proj=robin")

print(map_countries)

关键修改说明 #

• 关闭even.steps = TRUE，改为FALSE：让色条按断点定义的实际区间宽度显示，灰色区间（0-1）会呈现为窄条，与绿色渐变区间的比例匹配你的需求。
• 新增labels参数：自定义区间标签，让色条的含义更直观。
• 明确limits = c(0, Inf)：确保所有数值被正确纳入配色范围，避免异常值干扰。

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