How does the function pdegplot(model,'FaceLabels','on') detect the labeled Faces ?

Question

Sleh Eddine Brika 2016-5-17

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/284518-how-does-the-function-pdegplot-model-facelabels-on-detect-the-labeled-faces

回答： Damian Sheehy 2016-5-25

I am using PDE toolbox, when reading an STL file with importGeometry function, and plotted the result with the command pdegplot(model,'FaceLabels','on') I get this kind of result

Do you know how could I get the information of an arbitrary face, let's say I need to manipulate the face F4 how to access to the face F4? I tried to extract the faces, vertices and Edges like this E= [(hc(2).XData)' (hc(2).YData)' (hc(2).ZData)']%Edges V=hc(1).Vertices;%Vertices F=hc(1).Faces;%Faces but the faces that I extracted my self are not the same as the ones plotted. For exemple I have a square, the command pdegplot(model,'FaceLabels','on') will label one face, but the matrix F will have 2. I don't know how it's work ! I think that the fonction find a relation between faces/edges/vertices to detect the labeled faces but I didn't figure out how it's work.

N = 3; model = createpde(N); p=importGeometry(model,'cube.stl');

figure hc=pdegplot(model,'FaceLabels','on'); E= [(hc(2).XData)' (hc(2).YData)' (hc(2).ZData)']%Edge V=hc(1).Vertices; F=hc(1).Faces;

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Damian Sheehy 2016-5-25

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/284518-how-does-the-function-pdegplot-model-facelabels-on-detect-the-labeled-faces#answer_223399

The PDE toolbox does not support the editing of STL geometry at this time. You would need to edit the geometry in the CAD application where you created the file and then import the revision.

For a simple geometry like a cube you could create a triangulation representing the cube and then create geometry from that triangulation. You could then edit the triangulation in MATLAB. For example:

[X, Y, Z] = meshgrid(0:1);

P = [X(:), Y(:), Z(:)];

K = convhull(P);

pdem = createpde;

geometryFromMesh(pdem,P',K')

pdegplot(pdem)

You could then manipulate P.