ENSEEIHT/TP-modelisation-images
1
0
Fork 0
Code Issues Pull requests Packages Projects Releases Wiki Activity

feat: merged TP_maillage_suite.m into TP_maillage.m

Browse source
This commit is contained in:
Laurent Fainsin 2022-04-06 23:24:40 +02:00
parent e7a390af90
commit fcd1df58b1
2 changed files with 20 additions and 21 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

20
TP_maillage.m
View file

@ -337,4 +337,24 @@ view(80, -10);
% Sauvegarde des donnees
save donnees;
%% Affichage du maillage, au lieu des surfaces
% Calcul des faces du maillage à garder
FACES = [sort(triBis(:, [2 3 4]), 2); sort(triBis(:, [1 3 4]), 2); sort(triBis(:, [1 2 4]), 2); sort(triBis(:, [1 2 3]), 2)];
FACES = sortrows(FACES);
rep = sum(FACES(1:end - 1, :) == FACES(2:end, :), 2) == 3;
FACES([0; rep] | [rep; 0], :) = [];
fprintf('Calcul du maillage final termine : %d faces. \n', size(FACES, 1));
% Affichage du maillage final
figure(9);
hold on;
for i = 1:size(FACES, 1)
plot3([X(1, FACES(i, 1)) X(1, FACES(i, 2))], [X(2, FACES(i, 1)) X(2, FACES(i, 2))], [X(3, FACES(i, 1)) X(3, FACES(i, 2))], 'r');
plot3([X(1, FACES(i, 1)) X(1, FACES(i, 3))], [X(2, FACES(i, 1)) X(2, FACES(i, 3))], [X(3, FACES(i, 1)) X(3, FACES(i, 3))], 'r');
plot3([X(1, FACES(i, 3)) X(1, FACES(i, 2))], [X(2, FACES(i, 3)) X(2, FACES(i, 2))], [X(3, FACES(i, 3)) X(3, FACES(i, 2))], 'r');
end;
view(80, -10);

21
TP_maillage_suite.m
View file

@ -1,21 +0,0 @@
load donnees;
% Calcul des faces du maillage à garder
FACES = [sort(triBis(:, [2 3 4]), 2); sort(triBis(:, [1 3 4]), 2); sort(triBis(:, [1 2 4]), 2); sort(triBis(:, [1 2 3]), 2)];
FACES = sortrows(FACES);
rep = sum(FACES(1:end - 1, :) == FACES(2:end, :), 2) == 3;
FACES([0; rep] | [rep; 0], :) = [];
fprintf('Calcul du maillage final termine : %d faces. \n', size(FACES, 1));
% Affichage du maillage final
figure;
hold on
for i = 1:size(FACES, 1)
plot3([X(1, FACES(i, 1)) X(1, FACES(i, 2))], [X(2, FACES(i, 1)) X(2, FACES(i, 2))], [X(3, FACES(i, 1)) X(3, FACES(i, 2))], 'r');
plot3([X(1, FACES(i, 1)) X(1, FACES(i, 3))], [X(2, FACES(i, 1)) X(2, FACES(i, 3))], [X(3, FACES(i, 1)) X(3, FACES(i, 3))], 'r');
plot3([X(1, FACES(i, 3)) X(1, FACES(i, 2))], [X(2, FACES(i, 3)) X(2, FACES(i, 2))], [X(3, FACES(i, 3)) X(3, FACES(i, 2))], 'r');
end;
view(80, -10);