如何从3D Delaunay三角剖分中获得三角形

我已经有了一个3D三角剖分dt（或者可能是四面体化）的3D点云，我希望找到一个很好的方法从中提取三角形。请注意，我知道如何获得四面体，其简单的dt.ConnectivityList，是否有一种从四面体中获得三角形的有效方法？每个三角形只能在列表中出现一次。如何从3D Delaunay三角剖分中获得三角形

目前我'做以下 - 这不过是痛苦的缓慢：

dt = delaunayTriangulation([X Y Z]); 
tetrahedra = dt.ConnectivityList; 

tris = cell(1, size(tetrahedra, 1)); % contains indices of tris in a tetra 
for tt=1:size(tetrahedra, 1) 
    vertIds = tetrahedra(tt, :); % vertex indices 

    vmask = logical([0 1 1 1]); 
    tris{tt} = [vertIds(circshift(vmask, [0 0 0 0])); 
       vertIds(circshift(vmask, [1 1 1 1])); 
       vertIds(circshift(vmask, [2 2 2 2])); 
       vertIds(circshift(vmask, [3 3 3 3]))]; 
end 

tris = unique(sort(cell2mat(tris'), 2), 'rows');

来源

2014-02-18 brtk

不，我看看你在那里使用的任何函数，但是你是否尝试将'tris'声明为一个数组而不是ac数组？我对大循环中的细胞有过“有趣”的体验。 – McMa

@McMa是的，这是一个可能的优化，但是，我仍然需要对每一行进行排序，然后删除模糊。 – brtk

这里有一个向量化版本：

% take all four subsets of three points from the tets and concatenate 
tris2 = [tetrahedra(:,[1 2 3]); tetrahedra(:,[1 2 4]); tetrahedra(:,[1 3 4]); tetrahedra(:, [2 3 4])]; 
% sort each row 
tris2 = sort(tris2, 2); 
% eliminate duplicates 
tris2 = unique(tris2, 'rows');

拥有10个测试点（随机数）我得到

Cellmethod: time=9.870982, numelements = 1344960 
Vectmethod: time=1.014797, numelements = 1344960

来源

2015-05-15 12:22:51 xenoclast

我知道你的应用程序可能太晚了，但我喜欢这个难题... – xenoclast

如何从3D Delaunay三角剖分中获得三角形

回答

相关问题