Arbre (6) dessiné avec une tortue MetaPost (L-Système végétation)

Auteur ou autrice : Jean-Michel Sarlat.

Mise en ligne le 7 février 2023

Image du résultat de l’exemple

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

Mots clés : L-systemsvégétationtortue

Cet exemple fait partie de la collection d’exemples L-Systems (végétation).

Fichiers


lsys10-1.pdf

56.24 KB

Télécharger l’archive complète

Fichiers auxiliaires

tortue.mp


% 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;