Changeset 37

Timestamp:

03/14/08 22:45:46 (8 months ago)

Author:

christophe

Message:

--

Files:

• trunk/mp-geo/docpdf.tex (modified) (3 diffs)
• trunk/mp-geo/mp-geo.mp (modified) (22 diffs)

trunk/mp-geo/docpdf.tex

@@ -139 +139 @@
-  un problème d'accès direct à des pays comme Cuba ou l'Australie. Un
-  autre souci à régler\ldots}. Enfin les lacs sont regroupés dans un
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
-compilation de la figure \ding{172} a été de 37~s alors que pour la
-figure \ding{173}, la compilation a duré 27~s.
@@ -269 +271 @@
-  \includegraphics[scale=0.7]{polesud}
-  \caption{L'immensité du pôle sud ($\theta=0$; $\phi=90$\degres).}
+\end{figure}
+Ooops ! une petite dernière\ldots\footnote{Ne cherchez pas après les
+  fonctionnalités pour la transparence, elles ne sont pas implémentées
+  dans le package; elles ne sont valables que dans des cas
+  particuliers\ldots}
+\begin{figure}[ht]
+  \centering
+  \includegraphics{globetransparent.png}
+  \caption{Un peu de transparence.}
-\end{figure}
-\section{Projections planes}
@@ -454 +465 @@
-\section{Historique}
-\begin{description}
+\item[14/03/2008] Version \texttt{0.58} : Ajout de fleuves et lacs
+  supplémentaires.
+\item[14/03/2008] Version \texttt{0.56} : Ajout des volcans :).
-\item[18/02/2008] Version \texttt{0.55} : Ajout des projections de Bonne.
-\item[18/02/2008] Version \texttt{0.54} : Modification mineure dans la

trunk/mp-geo/mp-geo.mp

@@ -31 +31 @@
-tan=(0.824,0.705,0.55);
-payscolor=tan;
+
+vardef flamme=reverse(2mm*(0.5,1){dir90}..2mm*(0.75,1.5)..{dir90}2mm*(1,2))--2mm*(0.5,1){dir=-90}..2mm*(1,0.5)..{dir=90}2mm*(1.5,1)--((2mm*(0.5,1){dir90}..2mm*(0.75,1.5)..{dir90}2mm*(1,2)) reflectedabout (2mm*(1,0),2mm*(1,0)+2mm*(0,1)))--cycle
+enddef;
-
-%Les marques
@@ -348 +351 @@
-enddef;
-
-
-vardef Lecturelacs=
-  nomfichiermul:=arborescence&"lacs.dat";
@@ -386 +388 @@
-enddef;
-
+vardef Lecturelacssup=
+  nomfichiermul:=arborescence&"lacssup.dat";
+  nblec:=scantokens readfrom nomfichiermul;
+  for p=1 upto nblec:
+    color Coord[],Pays[];
+    numeric ll;
+    ll:=0;
+    nbpts:=scantokens readfrom nomfichiermul;
+    for k=1 upto nbpts:
+      pair latlon;
+      latlon=scantokens readfrom nomfichiermul;
+      Coord[k]=rayon*(cosd(xpart(latlon))*cosd(ypart(latlon)),cosd(xpart(latlon))*sind(ypart(latlon)),sind(xpart(latlon)));
+      if ((xpart(latlon)>phim) and (xpart(latlon)<phip)):
+        if ProduitScalaire(Coord[k]-pte3,Oeil-pte3)>0:
+          ll:=ll+1;
+          Pays[k]=Coord[k];
+        else:
+          Pays[k]=2*Coord[k];
+        fi;
+      else:
+        Pays[k]=2*Coord[k];
+      fi;
+    endfor;
+    path pays;
+    if ll>0:
+      pays=Projgeo(Pays[1])
+      for l=2 upto nbpts:
+        --Projgeo(Pays[l])
+      endfor;
+      fill pays--cycle withcolor cielfonce;
+      draw pays;
+      clip currentpicture to cercles(pte3,pte1,pte3,pte1,pte4);
+    fi;
+  endfor;
+  closefrom nomfichiermul;
+enddef;
+
+vardef Lecturevolcans=
+  drawoptions(withcolor orange);
+  nomfichiermul:=arborescence&"volcans.dat";
+  nblec:=scantokens readfrom nomfichiermul;
+  show nblec;
+  for p=1 upto nblec:
+    color Coord[],Pays[];
+    pair latlon;
+    latlon=scantokens readfrom nomfichiermul;
+    Coord[p]=rayon*(cosd(xpart(latlon/60))*cosd(ypart(latlon/60)),cosd(xpart(latlon/60))*sind(ypart(latlon/60)),sind(xpart(latlon/60)));
+    if ((xpart(latlon/60)>phim) and (xpart(latlon/60)<phip)):
+      if ProduitScalaire(Coord[p]-pte3,Oeil-pte3)>0:
+        Pays[p]=Coord[p];
+        fill flamme shifted Projgeo(Pays[p]) withcolor 1/2[orange,jaune];
+      %else:
+%       Pays[k]=2*Coord[k];
+      fi;
+%    else:
+%      Pays[k]=2*Coord[k];
+    fi;
+  endfor;
+closefrom nomfichiermul;
+drawoptions();
+enddef;
+
-vardef Lecturerivieres=
-  nomfichiermul:=arborescence&"rivieres.dat";
@@ -399 +463 @@
-      Coord[k]=rayon*(cosd(xpart(latlon/60))*cosd(ypart(latlon/60)),cosd(xpart(latlon/60))*sind(ypart(latlon/60)),sind(xpart(latlon/60)));
-      if ((xpart(latlon/60)>phim) and (xpart(latlon/60)<phip)):
+        if ProduitScalaire(Coord[k]-pte3,Oeil-pte3)>0:
+          ll:=ll+1;
+          Pays[k]=Coord[k];
+        else:
+          Pays[k]=2*Coord[k];
+        fi;
+      else:
+        Pays[k]=2*Coord[k];
+      fi;
+    endfor;
+    path pays;
+    if ll>0:
+      pays=Projgeo(Pays[1])
+      for l=2 upto nbpts:
+        --Projgeo(Pays[l])
+      endfor;
+      draw pays withcolor cielfonce;
+      clip currentpicture to cercles(pte3,pte1,pte3,pte1,pte4);
+    fi;
+  endfor;
+  closefrom nomfichiermul;
+enddef;
+
+vardef Lecturerivieressup=
+  nomfichiermul:=arborescence&"fleuvessup.dat";
+  nblec:=scantokens readfrom nomfichiermul;
+  for p=1 upto nblec:
+    color Coord[],fond,Pays[];
+    numeric ll;
+    ll:=0;
+    nbpts:=scantokens readfrom nomfichiermul;
+    for k=1 upto nbpts:
+      pair latlon;
+      latlon=scantokens readfrom nomfichiermul;
+      Coord[k]=rayon*(cosd(xpart(latlon))*cosd(ypart(latlon)),cosd(xpart(latlon))*sind(ypart(latlon)),sind(xpart(latlon)));
+      if ((xpart(latlon)>phim) and (xpart(latlon)<phip)):
-        if ProduitScalaire(Coord[k]-pte3,Oeil-pte3)>0:
-          ll:=ll+1;
@@ -697 +797 @@
-enddef;
-
+vardef Lecturelacspsup=
+  nomfichiermul:=arborescence&"lacssup.dat";
+  nblec:=scantokens readfrom nomfichiermul;
+  for p=1 upto nblec:
+    pair Coord[];
+    nbpts:=scantokens readfrom nomfichiermul;
+    numeric ll;
+    ll=0;
+    for k=1 upto nbpts:
+      pair latlon;
+      latlon=scantokens readfrom nomfichiermul;
+      if projection="mercator":
+        Coord[k]=mercatorc(xpart(latlon),ypart(latlon));
+      elseif projection="simple":
+        Coord[k]=simplec(xpart(latlon),ypart(latlon));
+      elseif projection="cylindrique":
+        Coord[k]=cylindriquec(xpart(latlon),ypart(latlon));
+      elseif projection="bonne":
+        Coord[k]=bonnec(xpart(latlon),ypart(latlon));
+      elseif projection="coniqueh":
+        if (xpart(latlon)>0) or (xpart(latlon)=0):
+          ll:=ll+1;
+          Coord[ll]=coniquec(xpart(latlon),ypart(latlon),45);
+        fi;
+      elseif projection="coniqueb":
+        if (xpart(latlon)<0):
+          ll:=ll+1;
+          Coord[ll]=coniquec(xpart(latlon),ypart(latlon),-45);
+        fi;
+      fi;
+    endfor;
+    path lac;
+    if (projection="mercator") or (projection="cylindrique") or (projection="simple") or (projection="bonne"):
+      ll:=nbpts
+    fi;
+    if ll>0:
+      lac=Coord[1]
+      for l=2 upto ll:
+        --Coord[l]
+      endfor;
+      fill lac--cycle withcolor cielfonce;
+      draw lac;
+    fi;
+  endfor;
+  closefrom nomfichiermul;
+enddef;
+
-vardef Lecturerivieresp=
-  nomfichiermul:=arborescence&"rivieres.dat";
@@ -725 +872 @@
-          ll:=ll+1;
-          Coord[ll]=coniquec(xpart(latlon/60),ypart(latlon/60),-45);
+        fi;
+      fi;
+    endfor;
+    path riv;
+    if (projection="mercator") or (projection="cylindrique") or (projection="simple") or (projection="bonne"):
+      ll:=nbpts;
+    fi;
+    if ll>0:
+      riv=Coord[1]
+      for l=2 upto ll:
+        --Coord[l]
+      endfor;
+      draw riv withcolor cielfonce;
+    fi;
+  endfor;
+  closefrom nomfichiermul;
+enddef;
+
+vardef Lecturerivierespsup=
+  nomfichiermul:=arborescence&"fleuvessup.dat";
+  nblec:=scantokens readfrom nomfichiermul;
+  for p=1 upto nblec:
+    pair Coord[];
+    nbpts:=scantokens readfrom nomfichiermul;
+    numeric ll;
+    ll:=0;
+    for k=1 upto nbpts:
+      pair latlon;
+      latlon=scantokens readfrom nomfichiermul;
+      if projection="mercator":
+        Coord[k]=mercatorc(xpart(latlon),ypart(latlon));
+      elseif projection="simple":
+        Coord[k]=simplec(xpart(latlon),ypart(latlon));
+      elseif projection="cylindrique":
+        Coord[k]=cylindriquec(xpart(latlon),ypart(latlon));
+      elseif projection="bonne":
+        Coord[k]=bonnec(xpart(latlon),ypart(latlon));
+      elseif projection="coniqueh":
+        if (xpart(latlon)>0) or (xpart(latlon)=0):
+          ll:=ll+1;
+          Coord[ll]=coniquec(xpart(latlon),ypart(latlon),45);
+        fi;
+      elseif projection="coniqueb":
+        if (xpart(latlon)<0):
+          ll:=ll+1;
+          Coord[ll]=coniquec(xpart(latlon),ypart(latlon),-45);
-        fi;
-      fi;
@@ -927 +1120 @@
-  if lacs=true:
-    Lecturelacsp;
+    Lecturelacspsup;
-  fi;
-  Lectureilesp;
@@ -934 +1128 @@
-  if fleuves=true:
-    Lecturerivieresp;
+    Lecturerivierespsup;
-  fi;
-  if maillage=true:
@@ -954 +1149 @@
-  if lacs=true:
-    Lecturelacsp;
+    Lecturelacspsup;
-  fi;
-  Lectureilesp;
@@ -961 +1157 @@
-  if fleuves=true:
-    Lecturerivieresp;
+    Lecturerivierespsup;
-  fi;
-  if maillage=true:
@@ -981 +1178 @@
-  if lacs=true:
-    Lecturelacsp;
+    Lecturelacspsup;
-  fi;
-  Lectureilesp;
@@ -988 +1186 @@
-  if fleuves=true:
-    Lecturerivieresp;
+    Lecturerivierespsup;
-  fi;
-  if maillage=true:
@@ -1009 +1208 @@
-  if lacs=true:
-    Lecturelacsp;
+    Lecturelacspsup;
-  fi;
-  Lectureilesp;
@@ -1016 +1216 @@
-  if fleuves=true:
-    Lecturerivieresp;
+    Lecturerivierespsup;
-  fi;
-  if maillage=true:
@@ -1036 +1237 @@
-  if lacs=true:
-    Lecturelacsp;
+    Lecturelacspsup;
-  fi;
-  Lectureilesp;
@@ -1043 +1245 @@
-  if fleuves=true:
-    Lecturerivieresp;
+    Lecturerivierespsup;
-  fi;
-  drawoptions(withcolor gris);
@@ -1060 +1263 @@
-  if lacs=true:
-    Lecturelacsp;
+    Lecturelacspsup;
-  fi;
-  Lectureilesp;
@@ -1067 +1271 @@
-  if fleuves=true:
-    Lecturerivieresp;
+    Lecturerivierespsup;
-  fi;
-  drawoptions(withcolor gris);
@@ -1077 +1282 @@
-rayon:=2;
-
-
+,volcans
-fleuves=true;
-lacs=true;
@@ -1083 +1288 @@
-noncolore=false;
-maillage=false;
+volcans=false;
-
-vardef Mappemonde(expr longobs,latobs)=
@@ -1116 +1322 @@
-  Lecture("Ceurope.dat");
-  Lecture("Cafrique.dat");
+  if volcans=true:
+    Lecturevolcans;
+  fi;
-  if lacs=true:
-    Lecturelacs;
+    Lecturelacssup;
-  fi;
-  Lectureiles;
@@ -1125 +1335 @@
-  if fleuves=true:
-    Lecturerivieres;
+    Lecturerivieressup;
-  fi;
-  if maillage=true: