devoir4/graphique.m at master · racota/devoir4 · GitHub

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
function graphique(poso, xi, yi, zi, face)

% Trouve azimut et elevation pour view()
[~, ~, rb] = getConstantesCylindre();
a = -45;
elevation = acos(norm([(rb(1)-poso(1)) (rb(2)-poso(2))])/(norm(rb-poso)));
e =-rad2deg(elevation);

% Dessin du cylindre
[R, h, rc] = getConstantesCylindre();
[x,y,z] = cylinder;
figure
surf(x*R+rc(1),y*R+rc(2),z*h+(rc(3)-h/2),'FaceAlpha',0,'EdgeAlpha',.3)
axis([0 0.1 0 0.1 0 0.2])

hold on % permet d'ajouter d'autre chose au graphique deja existent (superposition de graphique)

% Ajout de l'observateur
scatter3(poso(1),poso(2),poso(3),'*');

% Ajout de l'image
%Changer les chiffres dans face par les couleurs
color = zeros(length(face),3);
s = zeros(length(face),1);
for i=1:length(face)
    s(i) = 36;
    switch face(i)
        case 1
            color(i,:)=[1 0 0];%'r';
        case 2
            color(i,:)=[0 1 1];%'c';
        case 3
            color(i,:)=[0 1 0];%'g';
        case 4
            color(i,:)=[1 1 0];%'y';
        case 5
            color(i,:)=[0 0 1];%'b';
        case 6
            color(i,:)=[1 0 1];%'m';
        otherwise
            disp('erreur face doivent avoir une valeur de 1 a 6')
    end
end

scatter3(xi,yi,zi,s,color,'.');

view(a,e);

hold off