ggplot2中连续图例设置与等高线平滑优化问题

如果用geom_contour_filled默认得到的是离散分箱图例，想要连续渐变的图例，有两种可行方案：

• 方案1：使用栅格插值填充（推荐）
  放弃geom_contour_filled，改用geom_raster实现连续渐变填充，再叠加等高线，代码示例：

  library(ggplot2)
  
  ggplot(data.mi10, aes(x = x, y = y)) +
    # 连续填充，interpolate=TRUE让颜色过渡更平滑
    geom_raster(aes(fill = z), interpolate = TRUE) +
    # 保持15组等高线
    geom_contour(aes(z = z), bins = 15, color = "white", size = 0.5) +
    # 连续型颜色刻度，可替换成你需要的配色（比如scale_fill_gradient）
    scale_fill_viridis_c(option = "plasma")
  

• 方案2：强制geom_contour_filled使用连续图例
  如果你坚持用geom_contour_filled，可以生成足够密集的连续断点，配合连续色条图例：

  z_min <- min(data.mi10$z, na.rm = TRUE)
  z_max <- max(data.mi10$z, na.rm = TRUE)
  # 生成100个连续断点，模拟连续效果
  continuous_breaks <- seq(z_min, z_max, length.out = 100)
  
  ggplot(data.mi10, aes(x = x, y = y, z = z)) +
    geom_contour_filled(breaks = continuous_breaks) +
    geom_contour(breaks = continuous_breaks, color = "white") +
    # 指定使用连续色条图例
    scale_fill_viridis_c(guide = "colorbar")
  

等高线的平滑度核心取决于数据的分辨率，通过插值生成更密集的网格是最有效的办法，步骤如下：

1. 先安装并加载插值工具包interp：

   install.packages("interp")
   library(interp)
   library(tidyr) # 用于数据格式转换
   

2. 对原始数据集进行插值，生成高分辨率网格：

   # 生成500x500的密集网格（可根据需求调整length.out数值）
   interp_result <- interp(
     x = data.mi10$x, y = data.mi10$y, z = data.mi10$z,
     xo = seq(min(data.mi10$x), max(data.mi10$x), length.out = 500),
     yo = seq(min(data.mi10$y), max(data.mi10$y), length.out = 500)
   )
   # 转换为ggplot可用的数据框格式
   interp_df <- as.data.frame(interp_result) %>%
     pivot_longer(cols = -x, names_to = "y", values_to = "z") %>%
     mutate(y = as.numeric(y))
   

3. 用插值后的数据集绘制等高线，保持bins=15：

   ggplot(interp_df, aes(x = x, y = y, z = z)) +
     geom_raster(aes(fill = z), interpolate = TRUE) +
     geom_contour(bins = 15, color = "white", size = 0.5) +
     scale_fill_viridis_c()
   

如果不想额外安装包，也可以尝试给geom_contour添加lineend = "round"参数让线条边缘更圆润，但这只是视觉上的小优化，不如插值的效果明显：

ggplot(data.mi10, aes(x = x, y = y, z = z)) +
  geom_contour_filled(bins = 15) +
  geom_contour(bins = 15, color = "white", size = 0.5, lineend = "round")

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