Discussion :

[rhum187]

Bonjour à tous,

j'ai poster un message sur ce forum il y a quelques jours pour un problème de vitesse, et vos réponses m'ont vraiment aidé.

Il se trouve que j'ai une deuxième question :

Comment puis-je faire pour actualiser une figure. Je m'explique :

j'affiche dans une figure un squelette humain (calculs de matrices de translation et rotation) pour une position, ou une autre et ça marche plutôt bien. (je vous joins la figure juste pour voir)

Et j'aimerais pouvoir l'animer, comme si c'était une vidéo, à savoir afficher les frames les une à la suite des autres.

Toutes vos remarques peuvent être intéressantes alors n'hésitez pas (même si vos remarques ne sont pas directement liées à ma question, je pense par exemple à une personne qui me dirait : tu devrais même une autre couleur..., je suis preneur).

Voici mon code même si vous ne pourrez pas l'exécuter car certaines variables dépendent d'autres programmes qui dépendent alors de mon fichier de données.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
 
  data{1,1}
  base=rus1{10,2}
 
for i=1:27
 
    if i==1 % DaveRoot
        q=[0,0,0]';
        s=R(0,0,0,base);
    elseif i==2 % Pelvis
        q=p(:,1);
        s=r(:,:,1);
    elseif i==3 % Spine1
        q=p(:,1);
        s=r(:,:,1); 
    elseif i==4 % Lefthip
        q=p(:,1);
        s=r(:,:,1); 
    elseif i==9 % Righthip
        q=p(:,1);
        s=r(:,:,1);
    elseif i==14 % Spine2
        q=p(:,i-11);
        s=r(:,:,i-11);
   elseif i==17 % Neck
        q=p(:,i-1);
        s=r(:,:,i-1);
    elseif i==20 % L Clavicle
        q=p(:,i-4);
        s=r(:,:,i-4);
    elseif i==24 % R Clavicle
        q=p(:,i-8);
        s=r(:,:,i-8);
 
    else
        q=p(:,i-1);
        s=r(:,:,i-1);
    end
 
    p(:,i)=q+s*[data{1,1}(i,1),data{1,1}(i,2),data{1,1}(i,3)]';
    r(:,:,i)=s*R(data{1,1}(i,4),data{1,1}(i,5),data{1,1}(i,6),base)*R(data{1,200}(i,4),data{1,200}(i,5),data{1,200}(i,6),base);
    v=data{1,1}(i,7)*r(:,2,i);
    line(p(1,i)+[0,v(1)],p(2,i)+[0,v(2)],p(3,i)+[0,v(3)],'color','k')
    pause(0);
 
   % line(p(1,i)+k*[0,r(1,1,i)],p(2,i)+k*[0,r(2,1,i)],p(3,i)+k*[0,r(3,1,i)],'color','r')
   % line(p(1,i)+k*[0,r(1,2,i)],p(2,i)+k*[0,r(2,2,i)],p(3,i)+k*[0,r(3,2,i)],'color','g')
   % line(p(1,i)+k*[0,r(1,3,i)],p(2,i)+k*[0,r(2,3,i)],p(3,i)+k*[0,r(3,3,i)],'color','b')
end
set(findobj('type','line'),'linewidth',2)
set(findobj('type','line','color','k'),'linewidth',2)
view(3)
box on
axis equal
xlabel('x')
ylabel('y')
zlabel('z')

Et quand même la fonction R qui est appelée à de nombreuses reprises (qui est lente en passant)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
 
function R=R(a,b,c,base)
 
 
cosd=inline('cos(v/360*2*pi)');
sind=inline('sin(v/360*2*pi)');
 
Rx=[1, 0, 0; 0, cosd(a), -sind(a); 0, sind(a), cosd(a)];
Ry=[cosd(b), 0, sind(b); 0, 1, 0; -sind(b), 0, cosd(b)];
Rz=[cosd(c), -sind(c), 0; sind(c), cosd(c), 0; 0, 0, 1];
 
 
if (strcmp(base,'XZY'))
R=Rx*Rz*Ry;
end
 
if (strcmp(base,'XYZ'))
R=Rx*Ry*Rz;
end
 
if (strcmp(base,'YXZ'))
R=Ry*Rx*Rz;
end
 
if (strcmp(base,'YZX'))
R=Ry*Rz*Rx;
end
 
if (strcmp(base,'ZYX'))
R=Rz*Ry*Rx;
end
 
if (strcmp(base,'ZXY'))
R=Rz*Rx*Ry;
end

Merci d'avance.

[Invité]

Bonjour,

Une première remarque:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
cosd=inline('cos(v/360*2*pi)');
sind=inline('sin(v/360*2*pi)');

serait plus rapide en utilisant des fonctions anonymes:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
cosd=@(x) cos(x/360*2*pi);
sind=@(x) sin(x/360*2*pi);

Ensuite à la place des innombrables if/elseif, je mettrais plutôt un switch.
Pour actualiser la figure, as-tu essayé la fonction doc drawnow ?

[rhum187]

Je viens d'essayer ce code et c'est pas trop mal mais je visualise toujours sur la figure les anciens squelettes donc j'ai une superposition...

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
 
data{1,1}
base=rus1{10,2}
 
for o=1:20:200
    for i=1:27
 
        if i==1 % DaveRoot
            q=[0,0,0]';
            s=R(0,0,0,base);
        elseif i==2 % Pelvis
            q=p(:,1);
            s=r(:,:,1);
        elseif i==3 % Spine1
            q=p(:,1);
            s=r(:,:,1);
        elseif i==4 % Lefthip
            q=p(:,1);
            s=r(:,:,1);
        elseif i==9 % Righthip
            q=p(:,1);
            s=r(:,:,1);
        elseif i==14 % Spine2
            q=p(:,i-11);
            s=r(:,:,i-11);
        elseif i==17 % Neck
            q=p(:,i-1);
            s=r(:,:,i-1);
        elseif i==20 % L Clavicle
            q=p(:,i-4);
            s=r(:,:,i-4);
        elseif i==24 % R Clavicle
            q=p(:,i-8);
            s=r(:,:,i-8);
 
        else
            q=p(:,i-1);
            s=r(:,:,i-1);
        end
 
        p(:,i)=q+s*[data{1,1}(i,1),data{1,1}(i,2),data{1,1}(i,3)]';
        r(:,:,i)=s*R(data{1,1}(i,4),data{1,1}(i,5),data{1,1}(i,6),base)*R(data{1,o}(i,4),data{1,o}(i,5),data{1,o}(i,6),base);
        v=data{1,1}(i,7)*r(:,2,i);
        line(p(1,i)+[0,v(1)],p(2,i)+[0,v(2)],p(3,i)+[0,v(3)],'color','k')
        pause(0);
 
        % line(p(1,i)+k*[0,r(1,1,i)],p(2,i)+k*[0,r(2,1,i)],p(3,i)+k*[0,r(3,1,i)],'color','r')
        % line(p(1,i)+k*[0,r(1,2,i)],p(2,i)+k*[0,r(2,2,i)],p(3,i)+k*[0,r(3,2,i)],'color','g')
        % line(p(1,i)+k*[0,r(1,3,i)],p(2,i)+k*[0,r(2,3,i)],p(3,i)+k*[0,r(3,3,i)],'color','b')
 
    end
 
    set(findobj('type','line'),'linewidth',2)
    set(findobj('type','line','color','k'),'linewidth',2)
    drawnow expose update
    view(3)
    box on
    axis equal
    xlabel('x')
    ylabel('y')
    zlabel('z')
 
end

[Invité]

Une petite précision, les lignes

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
    view(3)
    box on
    axis equal
    xlabel('x')
    ylabel('y')
    zlabel('z')

peuvent être mises au début en dehors de la boucle.
Ensuite, pour la superposition, tu peux travailler avec le handle de la ligne:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
h = line(x,y,'color','k','linewidth',2);
for o=1:20:200
    for i=1:27
        ...
        set(h,'xdata',p(1,i)+[0,v(1)],'ydata',p(2,i)+[0,v(2)],'zdata',p(3,i)+[0,v(3)])
        ...
    end
end

Ou ajouter un simple doc cla avant la ligne line(...)

[Jerome Briot]

Essaie ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
figure('doublebuffer','on')
 
view(3)
box on
axis equal
xlabel('x')
ylabel('y')
zlabel('z')
 
for i = 1:27   
    h(i) = line(nan,nan,nan);
end
 
set(h,'color','k','linewidth',2);
 
for o = 1:20:200
    for i = 1:27
 
        if i==1 % DaveRoot
            q = [0;0;0];
            s = R(0,0,0,base);
        elseif i==2 || i==3 || i==4 || i==9
            q = p(:,1);
            s = r(:,:,1);
        elseif i==14 % Spine2
            q = p(:,i-11);
            s = r(:,:,i-11);
        elseif i==20 % L Clavicle
            q = p(:,i-4);
            s = r(:,:,i-4);
        elseif i==24 % R Clavicle
            q = p(:,i-8);
            s = r(:,:,i-8);        
        else
            q = p(:,i-1);
            s = r(:,:,i-1);
        end
 
        p(:,i) = q+s*data{1,1}(i,1:3).';
        r(:,:,i) = s*R(data{1,1}(i,4),data{1,1}(i,5),data{1,1}(i,6),base)*R(data{1,o}(i,4),data{1,o}(i,5),data{1,o}(i,6),base);
        v = data{1,1}(i,7)*r(:,2,i);
 
        set(h(i),'xdata',p(1,i)+[0,v(1)],'ydata',p(2,i)+[0,v(2)],'zdata',p(3,i)+[0,v(3)])
 
    end
 
    drawnow
 
end

L'idée est inspirée de la : Comment animer le tracé d'une courbe ?

[rhum187]

Je remercie beaucoup tous les participants de cette conversation. Les réponses que vous m'avez fournis m'ont permis de terminé ce que j'avais entrepris.

Un grand merci !