用於數據中多個羣集的凸殼/凹殼

Question

我已經在羣集周圍繪製多邊形方面做了大量的閱讀，並且意識到可能是最好的前進方向。基本上我正在尋找一個彈性像多邊形來包裹我的集羣點。

我的數據是由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 

它提供了以下結果：