Matlab中alphaShape的z切片绘图出现异常凸起问题咨询

你遇到的这个z切片异常凸起问题，核心是两个细节没处理到位，我帮你梳理下原因和修复方案：

一、问题根源 #

1. 点集格式错误
   你定义的points是1行24个元素的行向量，但Matlab的alphaShape要求输入必须是N×3的矩阵——每行对应一个三维点的x、y、z坐标。直接传入行向量会让alphaShape错误解析点的结构，生成完全不符合你预期的形状，这必然导致后续切片检测出现奇怪的凸起。

2. alpha值默认设置不匹配
   alphaShape默认会自动计算一个alpha值来拟合点集，但对于你这种切顶金字塔（凸多面体）的场景，默认alpha值可能过小，导致形状生成时出现多余的面或错误连接。而你用convhull能得到预期形状，是因为convhull直接生成凸包，这和把alpha值设为无穷大时的alphaShape效果完全一致。

二、修复后的完整代码 #

我把你的代码调整成正确的格式，同时修复了alpha值的问题：

% 1. 正确构造N×3的三维点集矩阵（每行一个点）
raw_points = [0.3223 0.1761 0.3299;
              0.5752 1.2034 -0.0618;
              -0.0172 1.2034 -0.0618;
              0.2357 0.1761 0.3299;
              0.3223 0.0671 0.2209;
              0.5752 0.4588 -0.8064;
              -0.0172 0.4588 -0.8064;
              0.2357 0.0671 0.2209];
points = 1e3 * raw_points; % 缩放至实际尺寸

% 2. 设置alpha=Inf生成凸包形状，和convhull效果一致
alpha = alphaShape(points, Inf);

% 3. 检测z=198切片上的点是否在形状内
z = 198;
width = 550;
depth = 450;
inout = zeros(width, depth);

for qx = 1:width
    for qy = 1:depth
        if inShape(alpha, qx, qy, z)
            inout(qx, qy) = 1;
        end
    end
end

% 4. 绘图（转置是因为Matlab矩阵索引和绘图轴的对应关系，属于正常操作）
figure;
contourf(inout.');
xlabel('X 坐标');
ylabel('Y 坐标');
title('z=198 处的切片结果');

三、补充说明 #

• 关于contourf的转置：Matlab中矩阵的(行,列)索引对应绘图的(Y,X)轴，而你的inout是按(X方向点数, Y方向点数)创建的，转置后才能让X、Y轴和实际坐标一一对应，所以这个操作是合理的，不用纠结。
• 如果你后续需要处理非凸形状，可以调整alpha值，但对于凸多面体，Inf是最优选择，能保证和convhull完全一致的凸包结果。

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