How to color triplot faces?

9 次查看(过去 30 天)
Manar Al Asad
Manar Al Asad 2020-8-26
I have a 3D DTM in the form of triangulation (has faces and vertices) where the vertices coordinates are in cartesian coordinates (X,Y,Z). I would like to plot the DTM colored based on the radial position of each point.
In 3D this looks like the following, where r1 is a structure that has 2 fields: facets which a connectivity list, and vertices which is the point coordinates:
llr1(:,1),llr1(:,2),llr1(:,3)] = cart2sph(r1.vertices(:,1),...
r1.vertices(:,2),r1.vertices(:,3)); % convert to spherical coordinates to get radial position
llr1(:,1) = rad2deg(llr1(:,1)); % convert to degrees
llr1(:,2) = rad2deg(llr1(:,2)); % convert to degrees
figure(1)
trisurf(r1.obj.facets,r1.obj.vertices(:,1),r1.obj.vertices(:,2),...
r1.obj.vertices(:,3),llr1(:,3))
Now I would like to do the same thing in 2D simple cylindrical projection
figure(2);
triplot(r1.obj.facets,llr1(:,1),llr1(:,2))
gives me a line map with the right connectivity
but now I would like to color the faces in the 2d plot. Has anyone encountered this problem before and has a simple solution to it?

回答(2 个)

Cris LaPierre
Cris LaPierre 2020-8-27
I believe a triplot is a line plot. You can color the lines and markers (vertices), but there are no faces to color.

Nicolas Hadjittoouli
I managed to go around the problem by saving the points and connectivity list of triplot and I inserted an extra column of ones(m, 1) and use triangulation function to re-create the triangulation. Assuming that m is the length of your points.
For example:
h = triplot(X);
C = h.ConnectivityList;
P = h.Points;
P = [P ones(size(P,1), 1)]
TO = triangulation(C, P);
trisurf(TO);
And after that you can use trisurf as usual to insert face colors.

类别

Help CenterFile Exchange 中查找有关 Delaunay Triangulation 的更多信息

产品


版本

R2019a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by