How can I calculate the perimeters of Delaunay triangles?
조회 수: 19(최근 30일)
표시 이전 댓글
ZhG 2014년 4월 14일
댓글: Subrata Bhatacharjee 2019년 3월 10일
Assume I have a set of point coordinates in a 2D plane. x = [ 16 16 16 17 18 19 19 20 21 21 23 23 23 ] y = [ 151 369 397 208 84 93 177 326 112 243 164 213 390 ]
Then, I use these points to obtain the Delaunay triangles. I found in a thread which asked about the calculation of area of each triangle. But if I want to calculate the perimeter of each triangle, how can I do for that? If possible, I also want to know the vertices (coordiantes) of each triangle.
댓글 수: 0
Sven 2014년 4월 14일
Here is some code that shows a few different ways of getting the perimeter of triangles. It first just shows the lengths of each edge in the Delaunay Triangulation. Then it sums up the perimeter of the boundary of the triangulation. Then it shows the perimeter of each face (displayed in red).
x = [ 16 16 16 17 18 19 19 20 21 21 23 23 23 ]
y = [ 151 369 397 208 84 93 177 326 112 243 164 213 390 ]
DT = delaunayTriangulation(x',y')
figure, triplot(DT), hold on
edges = DT.edges;
edgeLens = zeros(size(edges,1),1);
for i = 1:size(edges,1)
thisEdgePts = DT.Points(edges(i,:),:);
edgeCpt = mean(thisEdgePts,1);
edgeLen = sqrt(sum(diff(thisEdgePts,,1).^2));
edgeLens(i) = edgeLen;
outerBoundaryPerim = sum(edgeLens(ismember(DT.edges, DT.freeBoundary,'rows')))
for i = 1:size(DT.ConnectivityList,1)
thisVerts = DT.Points(DT.ConnectivityList(i,:),:);
faceEdgeLens = sqrt(sum(diff(thisVerts([1:end 1],:),,1).^2,2));
facePerim = sum(faceEdgeLens);
faceCpt = mean(thisVerts,1);
Did that help you out and answer your question?
댓글 수: 3
Subrata Bhatacharjee 2019년 3월 10일
yes, I can see the perimeter in red color. But I want to save all triangle's perimeter and calculater sum and average. Can you please solve this problem?
Find more on Delaunay Triangulation in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!