在MATLAB中voronoi圖的顏色無限細胞

Question

我有兩組點和plot他們在藍色的星星和紅點。然後我plot沃羅諾伊圖與voronoi(X,Y)功能。我想指定每個單元格的顏色取決於它的網站所屬的集合。我幾乎被使用patch函數來完成這一個是這樣的：

 [v,c]=voronoin(D); 
    for p=1:TheNumberOfSets 
     r=rand()/2+0.5; % random gray color 
     col=[r r r]; 
     for s=1:PointsInSet(p) 
      l=l+1; 
      patch(v(c{l},1),v(c{l},2),col); % color 
      axis([0 10 0 10]); 
     end 
    end 

凡D是集點的座標，TheNumberOfSets表明我們有多少套確實有（在我們剛纔這個特殊的一部分2套），col指定一個隨機的灰色，PointsInSet指定我們在每個集合中有多少個點，並且使用l來枚舉Voronoi圖的單元。

，這是結果：

現在我的問題（你可以看到！）大約是無限的細胞。這段代碼只是改變了有界單元格的顏色，我想用它們在axis box範圍內的指定集合的​​顏色（即可以在圖像中看到的框）對無界單元格着色。

有什麼建議嗎？

Answer 1

您的示例確實爲無界單元創建了patch對象，但由於它們包含的頂點具有表示無界邊的Inf值，因此不會顯示它們。您需要用有限的頂點替換這些頂點來完成修補程序。

由voronoi生成的繪製無界單元邊的頂點對於此目的很有用。繪製Voronoi圖時，您可以獲取這些：

h = voronoi(D); 
v1 = shiftdim(reshape([h(2).XData;h(2).YData],2,3,[]),2); % Arranged one edge per row, one vertex per slice in the third dimension 

的無界邊緣最後繪製，因此您可以通過在c計數無限細胞中分離出這些：

nUnbounded = sum(cellfun(@(ic)ismember(1,ic),c)); 
v1Unbounded = v1(end-(nUnbounded-1):end,:,:); 

這些邊緣的第一家上市的頂點是細胞的有限頂點。這些並不總是完全匹配由於浮點錯誤由voronoin返回的座標，因此確定通過在從pdist2最小成對距離，編號爲頂點，這些對應於：

[~,iBounded] = min(pdist2(v,v1Unbounded(:,:,1))); % Index of the bounded vertex 
vUnbounded = v1Unbounded(:,:,2); % Displayed coordinate of the unbounded end of the cell edge 

要替換這些座標，你可以然後更換patch(v(c{l},1),v(c{l},2),col);下列要求：

cPatch = c{l}; % List of vertex indices 
vPatch = v(cPatch,:); % Vertex coordinates which may contain Inf 
idx = find(cPatch==1); % Check if cell has unbounded edges 
if idx 
    cPatch = circshift(cPatch,-idx); % Move the 1 to the end of the list of vertex indices 
    vPatch = [vPatch(1:idx-1,:) 
       vUnbounded(iBounded == cPatch(end-1),:) 
       vUnbounded(iBounded == cPatch(1),:) 
       vPatch(idx+1:end,:)]; % Replace Inf values at idx with coordinates from the unbounded edges that meet the two adjacent finite vertices 
end 
patch(vPatch(:,1),vPatch(:,2),col);