Цвет неограниченных клеток диаграммы ворона в MATLAB

Question

У меня есть два набора точек и plot их в синих звездах и красных точках. Затем I 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 используется для перечисления ячеек диаграммы Вороного.

и это результат:

теперь моя проблема (как вы можете видеть!) О неограниченных клетках. Этот код просто меняет цвет ограниченных ячеек, и я хочу раскрасить неограниченные ячейки с цветом их указанного набора в диапазоне окна оси (т. Е. Поле, которое вы видите на изображении).

Любое предложение?

Answer 1

Ваш пример действительно создает объекты для неограниченных ячеек, но поскольку они содержат вершины с 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);