等高线平滑方法：如何消除随机扰动引发的等高线棱角？

问题描述 #

用户通过以下Matlab代码生成含随机噪声的数据集P，经griddata插值后绘制contourf图，结果出现明显棱角。增大插值网格数N至2500后，效果无显著改善：

x=0:1:50;
r=3.*rand(1,51);
P=[cos(x')+r',sin(x')+r',sin(x'+r').*cos(x')+r'];
P = sortrows(P,[1,2]);
N = 500;
xv = linspace(min(P(:,1)), max(P(:,1)), N);
yv = linspace(min(P(:,2)), max(P(:,2)), N);
[X,Y] = ndgrid(xv, yv);
Z = griddata(P(:,1), P(:,2), P(:,3), X, Y);
contourf(X, Y, Z, 35)

解决方法 #

核心是针对噪声数据做平滑处理，而非单纯提升插值分辨率，以下是几种有效方案：

1. 更换griddata的插值方法 #

默认的线性插值（'linear'）对噪声敏感，易保留棱角。改用三次样条插值（'cubic'）或v4径向基插值，能生成更平滑的曲面：

% 三次样条插值，平滑度更高
Z = griddata(P(:,1), P(:,2), P(:,3), X, Y, 'cubic');

% 或使用Matlab传统v4插值方法，适合生成连续平滑曲面
Z = griddata(P(:,1), P(:,2), P(:,3), X, Y, 'v4');

2. 对原始数据先做降噪处理 #

在插值前对数据集P的各列做平滑滤波，消除随机噪声的影响，常用方法有：

• 移动平均滤波：适合简单降噪，可调整窗口大小控制平滑程度

  windowSize = 5; % 窗口大小，噪声越强可设越大
  P_smoothed = movmean(P, windowSize, 1);
  Z = griddata(P_smoothed(:,1), P_smoothed(:,2), P_smoothed(:,3), X, Y);
  

• Savitzky-Golay滤波：适合在保留数据趋势的同时降噪（需Signal Processing Toolbox）

  windowSize = 7; % 需为奇数
  order = 2; % 拟合多项式阶数
  P_smoothed = sgolayfilt(P, order, windowSize);
  Z = griddata(P_smoothed(:,1), P_smoothed(:,2), P_smoothed(:,3), X, Y);
  

3. 插值后再做二维平滑 #

若插值后仍有细微棱角，可对生成的Z矩阵做二维高斯滤波：

Z = griddata(P(:,1), P(:,2), P(:,3), X, Y, 'cubic');
Z_smoothed = imgaussfilt(Z, 2); % sigma值越大，平滑程度越高
contourf(X, Y, Z_smoothed, 35)

4. 优化contourf视觉效果 #

通过调整绘图参数弱化棱角感：

contourf(X, Y, Z_smoothed, 35, 'LineStyle', 'none'); % 关闭轮廓线
colormap('viridis'); % 选用过渡更自然的色图

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