创建许多凸包 MATLAB

Question

我是新来的，这是我的第一篇文章。我想知道是否有一种方法可以从 100 个随机点计算所有可能的凸包。我创建了一个正确的代码，但是它在 convhull 处给我一个错误，并且无法计算该部分下的所有内容..换句话说，我想知道是否有一种方法可以在 for 循环中使用函数 convhull 而不会出错..

这里是代码

hold on
N=100;
Points = zeros(N,2);
for i=1:N
    Points(i,1)= rand(1,1)*100;
    Points(i,2)= rand(1,1)*100;
end

SortedX = sortrows(Points,1);
NewVertices = SortedX;

X = reshape(SortedX(:,1),2,[]);
X = X(:);
Y = reshape(SortedX(:,2),2,[]);
Y = Y(:);
plot(X,Y,'*');

while length(NewVertices)>=3
    NewVertices = NewVertices(NewVertices~=0);
    rowsNV = length(NewVertices(:,1));
    NewVertices = reshape(NewVertices,rowsNV/2,2);
    XNV = reshape(NewVertices(:,1),rowsNV/2,[]);
    XNV = XNV(:);
    YNV = reshape(NewVertices(:,2),rowsNV/2,[]);
    YNV = YNV(:);
    K = convhull(XNV,YNV);
    plot(XNV(K),YNV(K),'r-');
    for i=1:length(K)-1
        NewVertices(K(i),1)=0;
        NewVertices(K(i),2)=0;
    end

end

这是错误

Error using convhull
Error computing the convex hull. Not
enough unique points specified.

Error in test2 (line 28)
    K = convhull(XNV,YNV);

提前致谢

co2ark5

这是在 DAN 的帮助下正确工作的最终代码

hold on
N=100;
Points = rand(N,2);

SortedX = sortrows(Points,1);
NewVertices = SortedX;

plot(SortedX(:,1),SortedX(:,2),'*');

while length(NewVertices)>=3
    X=NewVertices(:,1);
    Y=NewVertices(:,2);
    K = convhull(X,Y);
    plot(X(K),Y(K),'r-');
    for i=1:length(K)-1
        NewVertices(K(i),1)=0;
        NewVertices(K(i),2)=0;
    end
    NewVertices = NewVertices(any(NewVertices,2),:);
end

Answer 1

该错误相当具有表现力，您没有指定足够的点（即至少 3 个非共线点）。请说明您要使用此代码实现的目标。

您何时收到此错误？我假设不在循环的第一次运行？在检查是否有超过 3 个点后，您可能会立即删除更多点，然后调用 convhull 所以很有可能少于 3 个点（这仍然忽略了共线性的机会）。报错后size(NewVertices)是什么？

一些侧面的cmets：请在你绘制之前解释整个重塑业务（也许添加一些cmets）。 Points 也可以更快更简单地初始化：

Points = rand(N, 2);

编辑：

通过阅读下面的 cmets，我很清楚您的循环逻辑是错误的。我无法准确计算出您如何尝试删除点以创建新的子集以在其上找到船体，没有 cmets，我不明白您在使用 reshape 做什么，但我很确定您需要做的就是将循环的第一行移动到末尾，如下所示：

while length(NewVertices)>=3

    rowsNV = length(NewVertices(:,1));
    NewVertices = reshape(NewVertices,rowsNV/2,2);
    XNV = reshape(NewVertices(:,1),rowsNV/2,[]);
    XNV = XNV(:);
    YNV = reshape(NewVertices(:,2),rowsNV/2,[]);
    YNV = YNV(:);
    K = convhull(XNV,YNV);
    plot(XNV(K),YNV(K),'r-');
    for i=1:length(K)-1
        NewVertices(K(i),1)=0;
        NewVertices(K(i),2)=0;
    end

    NewVertices = NewVertices(NewVertices~=0);

end

这样您在检查点数后立即找到您的船体，而不是检查，然后删除点，这就是您的循环在 1 - 2 个点上运行的原因。

Answer 2

非常感谢丹，

首先，我尝试只更改您所说的那一行，但由于我正在使用的重塑功能，它没有再次起作用。然后我决定更改它们，我让它变得更简单并且它确实有效！

这是在 DAN 的帮助下正确工作的最终代码

hold on
N=100;
Points = rand(N,2);

SortedX = sortrows(Points,1);
NewVertices = SortedX;

plot(SortedX(:,1),SortedX(:,2),'*');

while length(NewVertices)>=3
    X=NewVertices(:,1);
    Y=NewVertices(:,2);
    K = convhull(X,Y);
    plot(X(K),Y(K),'r-');
    for i=1:length(K)-1
        NewVertices(K(i),1)=0;
        NewVertices(K(i),2)=0;
    end
    NewVertices = NewVertices(any(NewVertices,2),:);
end