Discussion :

[Reya04]

Bonsoir!

Je suis nouveau ici et j'ai cherché pour un problème similaire au mien, sans succès.
Voici mon problème:

J'ai une plaque avec des mailles (Plan X-Y dans mon cas) et je veux tracer un quart de cercle dessus. En fait, ce ne sera pas un cercle parfait puisque ce sont les mailles de la plaque qui vont dicter le passage du cercle.
Pour mieux comprendre (car je pense que mon introduction reste un peu obscure), voici l'ébauche que j'ai effectué:

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
close all
clear all
clc
 
dx = 1;
dy = 2;
r = 10;
 
for x=0:dx:r
 
     y = sqrt(r.^2-x^2);
     a = round(y/dy);
     b = a*dy;
 
     plot(x,b,'o')
     grid on
     hold on
 
end
 
% cercle réel (uniquement a titre informatif)
 
z = 0:dx:r;
d = sqrt(r.^2-z.^2);
plot(z,d,'r')

Voici donc le programme de base.
Un nœud (ou point) ne pouvant donc que se trouver sur une intersection du maillage (de la grille dans ce cas).
Maintenant, mon problème vient du fait que je souhaite rajouter des points lorsqu'il y a des ''sauts'' (je parle de l'axe des Y).
Dans un premier temps, je souhaite placer des points en dessous afin de ''compléter l'escalier''.
mais voilà! comment faire?

J'espère que j'ai été assez explicite, n'hésitez pas à me poser des questions si ce n'est pas clair.

Pour info, j'avais pensé à un if b ~= b-1 dans ma boucle ''for'' mais ça ne fonctionnait pas vraiment...

Merci beaucoup de votre aide!

Bonne soirée

[Kangourou]

salut,

essaie de regarder du côté de l'algorithme de Brensenham pour les cercles, j'ai l'impression que ca correspond à ce que tu cherches à faire.

A+

[Reya04]

Tout d'abord, merci pour ta réponse

Oui, il semblerait que cela corresponde à ce que j'essaie de faire. J'ai trouvé un programme MATLAB qui fait ma première partie mais sur un cercle complet. sans plus...

J'ai progressé dans mon code. Je ne suis pas sûr que cela soit très propre mais il me complète une partie de "l'escalier". Il me reste à trouve comment implémenter les points (8;5), (9;1),(9;2),(9;3). Une idée?

Le nouveau code est:

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
close all
clear all
clc
 
dx = 1;
dy = 1;
r = 10;
 
x = 0;
y = 1;
b = r;
 
 
 
for     x=0:dx:r
 
 
 
            y=sqrt(r.^2-x.^2);
            a=round(y/dy);
            b=a*dy;
 
            plot(x,b,'o')
            grid on
            hold on
 
            if b ~= b-1 
 
            y = sqrt(r.^2-(x+dx)^2);
            a = round(y/dy);
            b = a*dy;
 
            plot(x,b,'or')
            grid on
            hold on
 
            end
 
 
 
 
 
end
 
h = 0:dx:r;
g = sqrt(r.^2-h.^2);
plot(h,g,'g')

Merci de votre aide!

[pseudocode]

J'ai progressé dans mon code. Je ne suis pas sûr que cela soit très propre mais il me complète une partie de "l'escalier". Il me reste à trouve comment implémenter les points (8;5), (9;1),(9;2),(9;3). Une idée?

Tant que la pente est inférieure à 45°, on se déplace toujours de "1" case en horizontal et de "0 ou 1" case en vertical. Quand la pente est supérieure à 45°, c'est l'inverse (vertical <--> horizontal).

Le problème se réduit à trouver le critère qui décide si on se déplace verticalement de 0 ou de 1.

[Reya04]

J'ai toujours pas réussi mais je vous tiens au courant si je trouve quelque chose

[Jean Dumoncel]

Bonjour,

as-tu des contraintes pour le tracé de ta courbe? Ne pourrais-tu pas tracer une courbe de Bézier? Tes mailles seraient alors les points de contrôle de la courbe.