MATLAB에서 보로 노이 다이어그램의 무한한 셀

Question

푸른 점과 빨간 점에 두 점 집합이 있고 plot 점이 있습니다. 그렇다면 voronoi(X,Y) 기능을 가진 두 세트의 보로 노이 다이어 plot 다이어그램. 각 셀의 색상을 지정하려면 사이트가 속한 설정에 따라 다릅니다.

 [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은 우리가 가지고이 특정 부분에 (우리가 얼마나 많은 세트를 수행 쇼 : 나는 거의 patch 함수를 사용하여이 방법을이 일을했습니다 2 세트) col 임의의 회색 색상을 지정하면 PointsInSet은 각 세트에 얼마나 많은 포인트가 있는지 지정하고 l은 보로 노이 다이어그램의 셀을 나열하는 데 사용됩니다. (당신이 볼 수있는!) 지금

내 문제가 무제한의 세포에 관한 것입니다 :

는 이것은 결과입니다. 이 코드는 바운드 된 셀의 색상을 변경하기 만하고 축 상자 범위 (즉, 이미지에서 볼 수있는 상자)에서 지정된 세트의 색상으로 묶여 있지 않은 셀에 색상을 지정하려고합니다.

의견이 있으십니까?

Answer 1

예제에서는 실제로 무제한 셀에 대해 patch 개체를 만듭니다. 그러나 테두리가없는 가장자리를 나타내는 Inf 값을 가진 정점이 포함되어 있기 때문에이 개체는 표시되지 않습니다. 패치를 완료하려면이 정점을 유한 패치로 대체해야합니다.

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);