如何将带孔洞的散点数据绘制成平滑等高线图？

针对数据中心存在孔洞、使用tricontourf时孔洞被错误填充的问题，提供两个可行的解决思路：

方法一：约束三角剖分排除孔洞区域 #

利用Matlab的delaunayTriangulation添加约束边定义孔洞，让三角剖分自动避开孔洞区域，确保等高线渲染正确：

% 加载数据
Data = readmatrix("result_9_1.txt");
X = Data(:,1);
Y = Data(:,2);
Z = Data(:,3);

% 1. 提取孔洞边界点（示例为中心圆形孔洞，需根据实际数据调整筛选规则）
hole_idx = find((abs(X) <= 0.005) & (abs(Y) <= 0.005));
hole_X = X(hole_idx);
hole_Y = Y(hole_idx);
% 将边界点按极角排序，形成闭合环
[hole_theta, ~] = cart2pol(hole_X, hole_Y);
[~, sort_idx] = sort(hole_theta);
hole_boundary = [hole_X(sort_idx), hole_Y(sort_idx)];

% 2. 创建带约束的三角剖分
dt = delaunayTriangulation(X, Y);
% 补充闭合点，确保约束边形成闭环
hole_boundary = [hole_boundary; hole_boundary(1,:)];
dt = delaunayTriangulation(dt.Points, dt.ConnectivityList, hole_boundary);

% 3. 绘制带孔洞的填充等高线
tricontourf(dt, Z, 20); % 20为等高线层级数，可按需调整
xlim([-0.025,0.025]);
ylim([-0.025,0.025]);
colorbar;
axis equal;

方法二：后处理用Patch遮挡孔洞 #

如果不想修改三角剖分流程，可以先正常绘制等高线，再用与背景色一致的patch覆盖孔洞区域：

Data = readmatrix("result_9_1.txt");
X = Data(:,1);
Y = Data(:,2);
Z = Data(:,3);

% 绘制基础等高线
tricontourf(X, Y, Z, 20);
xlim([-0.025,0.025]);
ylim([-0.025,0.025]);
colorbar;
hold on;

% 创建孔洞形状的遮挡patch（示例为圆形，需匹配实际孔洞形状）
theta = linspace(0, 2*pi, 100);
hole_x = 0.005 * cos(theta);
hole_y = 0.005 * sin(theta);
patch(hole_x, hole_y, get(gca, 'Color'), 'EdgeColor', 'none');
hold off;

关键注意点 #

• 方法一中的孔洞边界筛选规则需根据你的实际数据调整，确保提取到准确的闭合边界；
• 若孔洞为不规则形状，可通过ginput工具手动选取边界点，保存为坐标数组后使用；
• 等高线层级数（代码中的20）可根据Z值范围调整，平衡平滑度与细节展示。

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