How to create a triangulation from a list of edges and list of nodes?

19 views (last 30 days)
I'm beginning with a triangulation so that I know this graph is actually composed of non-intersecting edges that define a set of triangles.
% Start with a random triangulation, get edge list.
dt = delaunayTriangulation(rand(101,2)); % Keep it small for now.
edges = unique(sort(... % Unique mx2 list of
[dt.ConnectivityList(:,1:2);... % all edges in dt.
dt.ConnectivityList(:,[3 1])],...
figure; triplot(dt);
% Convert dt to graph, plot.
g = graph(edges(:,1),edges(:,2));
figure; g.plot;
% We know that g.edges describe dt, how can we reconstruct it?

Answers (2)

Christine Tobler
Christine Tobler on 9 Aug 2018
The problem is that a graph is a more general data structure than a Delaunay triangulation. Many graphs do not represent a triangulation at all, so there is not direct way of doing this.
If you also save the locations of each point in the original triangulation, you could write an algorithm that goes through each point in the graph and finds all the triangles it's a part of.
Could you tell me why you want to retrieve the triangulation from the graph?
  1 Comment
Dominik Mattioli
Dominik Mattioli on 14 Aug 2018
Edited: Dominik Mattioli on 14 Aug 2018
Your suggestion is the most intuitive, however, I imagine that being expensive for very large graphs unless it is done very cleverly. I think that the constraints I've given to the problem help with the problem's complexity.
In a nutshell, I will be using the graph class to selectively removing edges from the triangulation (as a graph) to form quadrilaterals.
Basically: create triangulation -> convert to graph -> remove edges to form some quads -> convert back to tri/quad/mixed 'rangulation'.
P.S. I've coded a triangulation class that accommodates quadrilaterals.

Sign in to comment.

Johannes Korsawe
Johannes Korsawe on 14 Oct 2022
allCyclesWithLengthThree = allcycles(g, 'MaxCycleLength', 3);
connectivityList = cell2mat(allCyclesWithLengthThree);


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!