在ggplot中平滑XYZ填充等高线：解决锯齿边缘问题

一、ggplot平滑等高线方案 #

原始等高线出现锯齿，核心原因是等高线生成算法基于离散网格点，直接用loess平滑原始数据易导致形状失真。推荐先对数据进行高密度插值，再绘制等高线，既能保留整体形状，又能消除锯齿：

步骤1：安装并加载插值工具包 #

install.packages("akima")
library(akima)
library(tidyverse)

步骤2：对原始数据进行插值处理 #

将原始0.2间隔的网格，插值到更小的间隔（比如0.05），生成平滑的新数据集：

# 从长格式数据中提取x、y、z向量
x_raw <- theta_90_long$x
y_raw <- theta_90_long$y
z_raw <- theta_90_long$z

# 插值生成更密集的网格，间隔可按需调整
interp_data <- interp(x = x_raw, y = y_raw, z = z_raw, 
                      xo = seq(-10, 10, by = 0.05), 
                      yo = seq(0, 10, by = 0.05))  # 统一y轴顺序为0到10，避免矩阵翻转

# 将插值结果转为ggplot可用的长格式
smooth_df <- interp_data %>% 
  as.data.frame() %>% 
  pivot_longer(cols = -y, names_to = "x", values_to = "z") %>% 
  mutate(x = as.numeric(x))

步骤3：绘制平滑等高线 #

替换原始数据为插值后的smooth_df，其余参数保持不变：

k_diff <- c(seq(0, 0.15, 0.01), 0.2, seq(0.3, 1.1, 0.1))

diff_90_smooth <- smooth_df %>% 
  ggplot(aes(x = x, y = y, z = z)) +
  geom_contour_filled(breaks = k_diff) +
  scale_y_continuous(breaks = seq(-10, 10, 2)) +
  scale_x_continuous(breaks = seq(-10, 10, 2)) +
  scale_fill_viridis_d(direction = 1, option = "H") +
  geom_contour(breaks = k_diff, color = "black", linewidth = 0.05) +
  theme(panel.grid.minor = element_blank(),
        legend.text = element_text(size = 9),
        legend.title = element_text(size = 12, face = "bold"),
        legend.key.width = unit(0.3, "cm"),
        legend.key.height = unit(0.3, "cm")) +
  labs(fill = "K") +
  ylab("Y") +
  xlab("X") +
  guides(fill = guide_legend(reverse = TRUE))

diff_90_smooth

二、修复Plotly无输出问题 #

你的Plotly代码无输出，是因为矩阵维度错误：长格式数据中length(theta_90_long$y)是所有行的数量，而非y轴唯一值的数量，导致矩阵形状不匹配。修正如下：

library(plotly)

# 提取x、y的唯一值并排序，确保顺序正确
x_unique <- unique(theta_90_long$x) %>% sort()
y_unique <- unique(theta_90_long$y) %>% sort()

# 将长格式z值转为匹配维度的矩阵，byrow参数需根据原始数据排列方式调整
z_matrix <- matrix(theta_90_long$z, 
                   nrow = length(y_unique), 
                   ncol = length(x_unique),
                   byrow = TRUE)

fig2 <- plot_ly(
  type = "contour", 
  z = z_matrix,
  x = x_unique,  # 指定x轴刻度
  y = y_unique,  # 指定y轴刻度
  autocontour = FALSE,  # 手动指定断点，和ggplot保持一致
  contours = list(
    breaks = k_diff,  # 复用之前的k_diff，保证等高线一致
    line = list(smoothing = 0.85)
  ),
  colorscale = "Viridis"  # 和ggplot的配色风格统一
) 

fig2

关键注意事项 #

• 插值间隔越小，等高线越平滑，但计算量会增加，0.05的间隔在你的数据范围内是均衡选择。
• 若插值后出现NA值，可在interp()函数中添加extrap = TRUE参数进行外推，仅建议在数据边缘少量使用。
• Plotly中关闭自动等高线并复用k_diff，能保证和ggplot的等高线断点完全一致，避免视觉差异。

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