数据中多个集群的凸包/凹包

Question

我已经阅读了大量关于在集群周围绘制多边形的文章，并意识到 convhull 可能是最好的前进方式。基本上，我正在寻找一个类似弹性的多边形来环绕我的聚类点。

我的数据是由 x（第 1 列）和 y（第 2 列）点组成的矩阵，这些点按簇（第 3 列）分组。我有 700 个这样的集群，因此单独绘制每个集群是不可行的。

有没有办法分别为每个集群执行 convhull，然后将它们中的每一个绘制在一个图表上。

编辑

到目前为止，我编写的代码无法在每个单独的集群上运行凸包...

    [ndata, text, alldata] =  xlsread(fullfile(source_dir));

       [~, y] = sort(ndata(:,end));
       As = ndata(y,:);

       lon = As(:,1); 
       lat = As(:,2);
       cluster = As(:,3);

       %% To find number of points in a cluster (repetitions)
   rep = zeros(size(cluster));

   for j = 1:length(cluster)
       rep(j) = sum(cluster==cluster(j));
   end

   %% Less than 3 points in a cluster are filtered out

   x = lon (rep>3);
   y = lat (rep>3);
   z = cluster (rep>3);

   %% convex hull for each cluster plotted ....hold....then display all.

   figure
   hold on

   clusters = unique(z);

       for i = 1:length(z)

          k=convhull(x(z==clusters(i)), y(z==clusters(i)));

          plot(x, y, 'b.'); %# plot cluster points
          plot(x(k),y(k),'r-'); %# plots only k indices, giving the convex hull


      end

下面是正在显示的图像；

如果这个问题已经被问过，我很抱歉重复，但请引导我找到你认为合适的答案。

请任何人都可以帮助解决这个问题，但我真的很挣扎！

Answer 1

我将遍历所有集群并执行您已经编写的操作，并使用hold on 选项将所有图累积在同一个图中。像这样的：

% Generate three clouds of points in 2D:
c1 = bsxfun(@plus, 0.5 * randn(50,2), [1 3]);
c2 = bsxfun(@plus, 0.6 * randn(20,2), [0 0]);
c3 = bsxfun(@plus, 0.4 * randn(20,2), [1 1]);

data = [c1, ones(50,1); ...
        c2, 2*ones(20,1); ...
        c3, 3*ones(20,1)];

% Plot the data points with different colors

clf
plot(c1(:,1), c1(:,2),'r+', 'LineWidth', 2);
hold on
plot(c2(:,1), c2(:,2),'k+', 'LineWidth', 2);
plot(c3(:,1), c3(:,2),'b+', 'LineWidth', 2);

x = data(:,1);
y = data(:,2);
cluster = data(:,3);

clusters = unique(cluster);

for i = 1:length(clusters)
    px = x(cluster == clusters(i));
    py = y(cluster == clusters(i));
    if length(px) > 2
        k = convhull(px, py);
        plot(px(k), py(k), '-');
    end
end

它给出以下结果：