Auteur ou autrice : Jean-Michel Sarlat.
Mise en ligne le 7 février 2023
Arbres divers décrits par une tortue MetaPost (fichier tortue.mp
). C'est une reprise des L-Systems (https://fr.wikipedia.org/wiki/L-Syst%C3%A8me) codés pour ImageMagick en perl (on les trouve dans la distibution de perlmagick).
L’exemple est extrait de la base MetaPost du site Syracuse : https://syracuse.eu.org/lab/bmp.
Code
%@Auteur : Jean-Michel Sarlat
%@Année : 2005
input tortue.mp;
numeric d,iter;
d = 6 ;
iter = 9;
color marron;
marron = (.5,.1,0);
path gf;
gf = (0,0)--(1,.8)--(0,0)--(-.5,1)--(0,0)--(-.7,-.1)--(0,0)--(0,-1.2)--(0,0)--(.7,-.4)--(0,0);
def pas(expr x) = x * 2.7mm + .5mm enddef;
def epaisseur(expr x) = .15mm * x + .1mm enddef;
def couleur(expr x) = (x/iter)[green,marron] enddef;
vardef trace(expr a,b) =
draw a--b
withpen pencircle scaled e
withcolor c;
enddef;
%% lsys10
%% règle A => S[---LMA][++++B]
%% règle B => S[++LBg][--Cg]
%% règle C => S[-----LB]S[+MC]
%% règle g => ' ' % seulement exécutée au terme
%% règle L => S+S+S+S+S+S--cycle
vardef A (expr n,m) =
save h,e,c;
color c;
h := pas(m);
e := epaisseur(m);
c := couleur(m);
if m>0:
avance(n,h,trace);
blop(n,n+1);
tourne(n+1,-3*d);
L(n+1,m-1);
retourne(n+1);
A(n+1,m-1);
blop(n,n+1);
tourne(n+1,4*d);
B(n+1,m-1);
fi
enddef;
vardef B (expr n,m) =
save h,e,c;
color c;
h := pas(m);
e := epaisseur(m);
c := couleur(m);
if m>0:
avance(n,h,trace);
blop(n,n+1);
tourne(n+1,2*d);
L(n+1,m-1);
B(n+1,m-1);
g(n+1,m-1);
blop(n,n+1);
tourne(n+1,-2*d);
C(n+1,m-1);
g(n+1,m-1);
fi
enddef;
vardef C (expr n,m) =
save h,e,c;
color c;
h := pas(m);
e := epaisseur(m);
c := couleur(m);
if m>0:
avance(n,h,trace);
blop(n,n+1);
tourne(n+1,-5*d);
L(n+1,m-1);
B(n+1,m-1);
avance(n,h,trace);
blop(n,n+1);
tourne(n+1,d);
retourne(n+1);
C(n+1,m-1);
fi
enddef;
vardef L (expr n,m) =
save h,e,c;
color c;
h := pas(m/15);
e := epaisseur(m/10);
c := green;
if m>0:
blop(n,n+1);
for j=1 upto 5:
avance(n+1,h,trace);
tourne(n+1,d);
endfor
avance(n+1,h,trace);
retourne(n+1);
tourne(n+1,180);
avance(n+1,h,trace);
remplis(n+1,.9*green);
fi
enddef;
vardef g (expr n,m) =
if m=0:
draw (gf scaled 3mm) shifted tortue_xy[n]
withpen pencircle scaled 1pt
withcolor .2[red,black];
draw (gf scaled 2.7mm) shifted tortue_xy[n]
withpen pencircle scaled .2pt
withcolor black;
fi
enddef;
beginfig(1);
tortue(1,0,0,90,1);
A(1,iter);
endfig;
end
% package tortue.mp
%% Les tables (chemin,point actuel,orientation,direction)
path tortue_p[];
pair tortue_xy[] ;
numeric tortue_m[];
numeric tortue_a[];
%% L'enregistrement d'une tortue (initialisation)
def tortue(expr n, x , y , a , m ) =
tortue_xy[n] := (x,y) ;
tortue_a[n] := a ;
tortue_m[n] := m ;
tortue_p[n] := (x,y);
enddef;
%% La tortue <n> se dédouble, <m> à cet instant se superpose à <n>
def blop(expr n,m) =
tortue(m,(xpart tortue_xy[n]),(ypart tortue_xy[n]),tortue_a[n],tortue_m[n]);
enddef;
%% La tortue tourne (sur place)
def tourne(expr n, a ) =
tortue_a[n] := tortue_a[n] + a * tortue_m[n];
enddef;
%% La tortue se retourne, sa gauche est notre droite !
def retourne(expr n ) =
tortue_m[n] := -tortue_m[n];
enddef;
%% La tortue avance d'une certaine distance et fait ce qu'on lui demande
vardef avance(expr n,l)(suffix m) =
pair position;
position = tortue_xy[n] + ((l,0) rotated tortue_a[n]);
m(tortue_xy[n],position);
tortue_xy[n] := position;
tortue_p[n] := tortue_p[n]--position;
enddef;
%% Remplissage du polygone circonscrit par le chemin de la tortue
vardef remplis(expr n,c) =
fill tortue_p[n]--cycle withcolor c;
enddef;
%% Une chose que la tortue peut faire en se déplaçant : laisser une trace !
vardef trace(expr a,b) =
draw a--b;
enddef;