根据初始点的颜色值在 MATLAB 中为 voronoi 图着色

Question

我有一个名为包含其 x 和 y 坐标的点矩阵start_coord，以及表示其分类 (1-5) 的列。即第一行看起来像 [75, 100, 4]。

我使用下面的代码计算了该数据的 voronoi 图

[vc_x, vc_y] = voronoi(start_coord(:,1), start_coord(:,2));

我将如何根据每个多边形中包含的点的分类值（即中的第三列）对结果多边形进行着色start_coord？

编辑 要按颜色快速绘制多边形，请参考下面评论中的答案，这有助于告知此编辑。要获取写入可保存为图像的数组的数千个点的 voronoi 多边形，请参阅以下代码：

new_map = zeros(sm_size(1), sm_size(2));
        start_coord = readmatrix(char(join([csv_path, '/', run_types(run), common_name_csv], "")));
        sc_size = size(start_coord);
        
        dt = delaunayTriangulation(start_coord(:,1:2));
        [V,R] = voronoiDiagram(dt);
                
        for i = 1:sc_size(1)
            A=V(R{i},:);
            B=A(any(~isinf(A),2),:); % omit points at infinity
            bw = poly2mask(B(:,1), B(:,2), sm_size(1), sm_size(2));
            new_map(bw == 1) = color_map(start_coord(i,3));
        end

new_map然后可以保存为数组或转换为 RGB 并保存为图像。

Answer 1

用于voronoiDiagram获取多边形。

dt = delaunayTriangulation(start_coord(:,1:2));
[V,R] = voronoiDiagram(dt);

那么 R{i} 将是 start_coord(i,:) 中多边形的顶点，因此将颜色设置为 start_coord(i,3) 的颜色并且：

A=V(R{i},:);
B=A(any(~isinf(A),2),:); % omit points at infinity
plot(polyshape(B));

唯一的问题是无穷远处的顶点被切掉了。但也许这会让你足够接近你想要的。如果您需要填充到边缘，请查看VoronoiLimit函数（我尚未测试过）。

例如：

X = [-1.5 3.2; 1.8 3.3; -3.7 1.5; -1.5 1.3; ...
0.8 1.2; 3.3 1.5; -4.0 -1.0;-2.3 -0.7; ...
0 -0.5; 2.0 -1.5; 3.7 -0.8; -3.5 -2.9; ...
-0.9 -3.9; 2.0 -3.5; 3.5 -2.25];

X(:,3) = [ 1 2 1 3 1 2 2 2 2 3 3 3 3 3 3]';

ccode = ["red","green","blue"];

dt = delaunayTriangulation(X(:,1:2));
[V,R] = voronoiDiagram(dt);
figure
voronoi(X(:,1),X(:,2))
hold on
for i = 1:size(X,1)
    A=V(R{i},:);
    B=A(any(~isinf(A),2),:);
    if(size(B,1)>2)
        plot(polyshape(B),'FaceColor',ccode(X(i,3)));
    end
end

结果：