Hello Giometar,
I made the following modifications to your code:
- Removed the call to 'VoronoiLimit' and instead used MATLAB's built-in 'voronoiDiagram' function directly with delaunayTriangulation. This is a direct approach to get the Voronoi vertices V and cells C.
- Simplified the plotting and labeling by using a single loop to plot both the point labels (T1, T2, ..., Tn) and their corresponding area values. This reduces redundancy and makes the code more concise.
- Consolidated the plotting commands with a single hold on and hold off pair, ensuring that all plotting operations happen within the same context and reducing potential errors.
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];
[VX,VY] = voronoi(X(:,1),X(:,2));
h = plot(VX,VY,'-b',X(:,1),X(:,2),'.r');
xlim([-4,4])
ylim([-4,4])
[V,C] = voronoiDiagram(delaunayTriangulation(X));
% Compute area of each Voronoi polygon
A = zeros(length(C),1);
for i = 1:length(C)
v1 = V(C{i},1);
v2 = V(C{i},2);
patch(v1,v2,rand(1,3));
A(i) = polyarea(v1,v2);
end
% Plot points and labels
hold on
plot(X(:,1), X(:,2), '.r');
for i = 1:length(X)
text(X(i,1), X(i,2), sprintf('T%d', i), 'FontWeight', 'bold', 'HorizontalAlignment', 'right');
text(X(i,1), X(i,2), sprintf('%.4f', A(i)), 'FontSize', 8, 'HorizontalAlignment', 'left');
end
hold off
Please use the documentation links below to know more about the functions 'voronoiDiagram' and 'delaunayTriangulation':
- voronoiDiagram: https://www.mathworks.com/help/matlab/ref/delaunaytriangulation.voronoidiagram.html
- delaunayTriangulation:https://www.mathworks.com/help/matlab/ref/delaunaytriangulation.html
