Delaunay 图中的角度计算

Angle calculation in a Delaunay graph

我在 MATLAB 中绘制了一个像这样的 Delaunay 图:-

我想计算图中所有的角度。我有所有无序形式的点的 x 和 y 值,但我不知道如何对这些点进行排序,因为同一行上的那些点的 x 和 y 值很接近。

一种方法:

x = randn(1,4)*10;
y = randn(1,4)*10;

%calculate the triangulation
tri = delaunay(x,y);

%Plot the graph
triplot(tri,x,y)
hold on

plot(x,y,'ro')
text(x,y,strsplit(num2str(1:length(x))))


% Determine each angle
for i = 1:size(tri,1)
    per = perms(tri(i,:));
    [~, ind] = unique(per(:,2)); %avoid to calculate two time the same angle.
    per = per(ind,:); %the 3 * 3 points that create the angle of each triangle
    for j = 1:3
       P_1 = per(j,1);
       P1 = [x(P_1),y(P_1)];
       P_2 = per(j,2);
       P2 = [x(P_2),y(P_2)];
       P_3 = per(j,3);
       P3 = [x(P_3),y(P_3)];
       ANG = rad2deg(atan2(abs(det([P3-P2;P1-P2])),dot(P3-P2,P1-P2))); %P2 is the point in the middle
       fprintf('Node %d %d %d angle %f\n',P_1, P_2, P_3, ANG)
    end
end