| 384 | | \subsection{Projection conique} |
|---|
| 385 | | Une projection conique est disponible : il s'agit de la {\em |
|---|
| 386 | | projection conique équidistante tangente}\footnote{D'après |
|---|
| 387 | | \url{http://thierry.hatt.gps.free.fr/projections/images-proj/coniq-equid-lambert.htm}. Pour un point de la sphère terrestre de coordonnées $(\phi;\lambda)$, les formules de projection sont |
|---|
| 388 | | \[\left\{\begin{array}{l} |
|---|
| 389 | | x=\dfrac{\cos\phi\sin(\lambda\sin(45))}{\sin{45}\cos(\phi-45)}\\ |
|---|
| 390 | | \\ |
|---|
| 391 | | y=\dfrac{\cos\phi\cos(\lambda\sin(45))}{\sin{45}\cos(\phi-45)}\\ |
|---|
| 392 | | \end{array} |
|---|
| 393 | | \right. |
|---|
| 394 | | \] |
|---|
| 395 | | }; le parallèle de contact étant soit le 45\degres\ Nord soit le |
|---|
| 396 | | 45\degres\ Sud. On choisit l'hémisphère à représenter par le code |
|---|
| 397 | | source |
|---|
| 398 | | \begin{verbatim} |
|---|
| 399 | | input mp-geo; |
|---|
| 400 | | |
|---|
| 401 | | figure(-20u,-20u,20u,20u); |
|---|
| 402 | | ConiqueH; |
|---|
| 403 | | fin; |
|---|
| 404 | | |
|---|
| 405 | | end |
|---|
| 406 | | \end{verbatim} |
|---|
| 407 | | pour obtenir |
|---|
| 408 | | \begin{figure}[ht] |
|---|
| 409 | | \centering |
|---|
| 410 | | \includegraphics[scale=0.4]{coniqueHN-1.png} |
|---|
| 411 | | \caption{Projection conique - Hémisphère nord.} |
|---|
| 412 | | \end{figure} |
|---|
| 413 | | ou alors |
|---|
| 414 | | \begin{verbatim} |
|---|
| 415 | | input mp-geo; |
|---|
| 416 | | |
|---|
| 417 | | figure(-20u,-20u,20u,20u); |
|---|
| 418 | | ConiqueB; |
|---|
| 419 | | fin; |
|---|
| 420 | | |
|---|
| 421 | | end |
|---|
| 422 | | \end{verbatim} |
|---|
| 423 | | pour voir l'hémisphère sud. |
|---|
| 424 | | \begin{figure}[ht] |
|---|
| 425 | | \centering |
|---|
| 426 | | \includegraphics[scale=0.4]{coniqueBS-1.png} |
|---|
| 427 | | \caption{Projection conique - Hémisphère sud.} |
|---|
| 428 | | \end{figure} |
|---|
| 429 | | \par\`A noter que cette fois-ci, le maillage n'est pas |
|---|
| 430 | | facultatif\footnote{On pourrait peut-être ajouter des graduations}. |
|---|
| | 416 | \subsection{Nouvelle méthode} |
|---|
| | 417 | Avec le développement de ces projections, on s'aperçoit qu'elles ont |
|---|
| | 418 | toutes le même fonctionnement en commun. Aussi, depuis la version |
|---|
| | 419 | \texttt{0.60}, une autre méthode\footnote{Ne vous inquiétez pas : |
|---|
| | 420 | pour des soucis de compatibilité, les autres moyens de projection |
|---|
| | 421 | ont été conservés.} a été implantée. |
|---|
| | 422 | \par En effet, il me semble plus simple de procéder de la manière |
|---|
| | 423 | suivante : |
|---|
| | 424 | \begin{verbatim} |
|---|
| | 425 | input mp-geo; |
|---|
| | 426 | |
|---|
| | 427 | figure(-10u,-10u,10u,10.5u); |
|---|
| | 428 | projection:="bonne"; |
|---|
| | 429 | Projection(0,41,50); |
|---|
| | 430 | Lecturevillesp("France"); |
|---|
| | 431 | fin; |
|---|
| | 432 | |
|---|
| | 433 | end |
|---|
| | 434 | \end{verbatim} |
|---|
| | 435 | Ceci afin d'alléger un peu les commandes à connaître. Un autre |
|---|
| | 436 | avantage est de choisir simplement le mode de projection choisi. |
|---|
| | 437 | \par Pour la création des cartes, il fallait, jusqu'à cette version |
|---|
| | 438 | \texttt{0.60}, parcourir tous les points de la planète pour ne garder |
|---|
| | 439 | que ceux intéressants.\\Dorénavant, les paramètres \verb!Amnord!, |
|---|
| | 440 | \verb!Amsud!, \verb!Amcentrale!, \verb!Caraibes!, \verb!Asie!, |
|---|
| | 441 | \verb!Europe!, \verb!Afrique! et \verb!All! font leurs |
|---|
| | 442 | apparitions. Par défaut, ils sont tous positionnés à \verb!false! sauf |
|---|
| | 443 | \verb!All!. Autrement, le comportement par défaut est celui qui prend |
|---|
| | 444 | le plus de temps ! |
|---|
| | 445 | Aussi, on prendra soin de changer certains de ces positionnements pour |
|---|
| | 446 | obtenir des temps de compilation moins grands. |
|---|
| | 447 | \subsection{Cartes et échelles} |
|---|
| | 448 | Pour pouvoir lire correctement une carte obtenue par |
|---|
| | 449 | projection\footnote{Uniquement les projections déjà vues, c'est-à-dire |
|---|
| | 450 | \verb!simple!, \verb!cylindrique!, \verb!mercator! et \verb!bonne!}, il |
|---|
| | 451 | m'a semblé nécessaire d'y adjoindre une échelle. On ne peut avoir |
|---|
| | 452 | qu'une échelle moyenne; sinon il faudrait déterminer l'échelle de |
|---|
| | 453 | chaque parallèle ou méridien. |
|---|
| | 454 | \par Aussi, la méthode choisie est la suivante : |
|---|
| | 455 | \begin{itemize} |
|---|
| | 456 | \item Supposons que la projection soit centrée sur le point |
|---|
| | 457 | $(\lambda_0;\theta_0)$; on détermine alors la longueur {\em réelle} |
|---|
| | 458 | de l'arc de parallèle situé entre les points $(\lambda_0;\theta_0)$ |
|---|
| | 459 | et $(\lambda_0;\theta_0+5\mbox{\degres})$; |
|---|
| | 460 | \item on détermine ensuite cette même longueur mais {\em sur le dessin}; |
|---|
| | 461 | \item on termine par une simple relation de proportionnalité. |
|---|
| | 462 | \end{itemize} |
|---|
| | 463 | On utilisera donc la syntaxe \verb!echelle(2,46.5,50)! pour un calcul |
|---|
| | 464 | d'échelle de 1~cm pour 50~km pour une projection centrée sur le point $(46,5;2)$. |
|---|
| | 465 | \begin{verbatim} |
|---|
| | 466 | input mp-geo; |
|---|
| | 467 | |
|---|
| | 468 | All:=false; |
|---|
| | 469 | Europe:=true; |
|---|
| | 470 | |
|---|
| | 471 | figure(-10u,-10u,10u,10.5u); |
|---|
| | 472 | projection:="bonne"; |
|---|
| | 473 | Projection(2,46.5,echelle(2,46.5,50)); |
|---|
| | 474 | fin; |
|---|
| | 475 | |
|---|
| | 476 | end |
|---|
| | 477 | \end{verbatim} |
|---|
| | 478 | pour obtenir la figure |
|---|
| | 479 | \begin{figure}[ht] |
|---|
| | 480 | \centering |
|---|
| | 481 | \includegraphics{projectionbonne} |
|---|
| | 482 | \caption{Projection de Bonne -- Centrée sur Paris -- \'Echelle $1\,:\,500\,000$.} |
|---|
| | 483 | \end{figure} |
|---|
| | 484 | \subsection{Projection conique} |
|---|
| | 485 | Cette projection est un peu spéciale, je n'ai pas réussi à l'inclure |
|---|
| | 486 | avec les autres. La projection conique disponible est la {\em |
|---|
| | 487 | projection conique équidistante tangente}\footnote{D'après |
|---|
| | 488 | \url{http://thierry.hatt.gps.free.fr/projections/images-proj/coniq-equid-lambert.htm}. Pour un point de la sphère terrestre de coordonnées $(\phi;\lambda)$, les formules de projection sont |
|---|
| | 489 | \[\left\{\begin{array}{l} |
|---|
| | 490 | x=\dfrac{\cos\phi\sin(\lambda\sin(45))}{\sin{45}\cos(\phi-45)}\\ |
|---|
| | 491 | \\ |
|---|
| | 492 | y=\dfrac{\cos\phi\cos(\lambda\sin(45))}{\sin{45}\cos(\phi-45)}\\ |
|---|
| | 493 | \end{array} |
|---|
| | 494 | \right. |
|---|
| | 495 | \] |
|---|
| | 496 | }; le parallèle de contact étant soit le 45\degres\ Nord soit le |
|---|
| | 497 | 45\degres\ Sud. On choisit l'hémisphère à représenter par le code |
|---|
| | 498 | source |
|---|
| | 499 | \begin{verbatim} |
|---|
| | 500 | input mp-geo; |
|---|
| | 501 | |
|---|
| | 502 | figure(-20u,-20u,20u,20u); |
|---|
| | 503 | ConiqueH; |
|---|
| | 504 | fin; |
|---|
| | 505 | |
|---|
| | 506 | end |
|---|
| | 507 | \end{verbatim} |
|---|
| | 508 | pour obtenir |
|---|
| | 509 | \begin{figure}[ht] |
|---|
| | 510 | \centering |
|---|
| | 511 | \includegraphics{coniqueHN.png} |
|---|
| | 512 | \caption{Projection conique - Hémisphère nord.} |
|---|
| | 513 | \end{figure} |
|---|
| | 514 | ou alors |
|---|
| | 515 | \begin{verbatim} |
|---|
| | 516 | input mp-geo; |
|---|
| | 517 | |
|---|
| | 518 | figure(-20u,-20u,20u,20u); |
|---|
| | 519 | ConiqueB; |
|---|
| | 520 | fin; |
|---|
| | 521 | |
|---|
| | 522 | end |
|---|
| | 523 | \end{verbatim} |
|---|
| | 524 | pour voir l'hémisphère sud. |
|---|
| | 525 | \begin{figure}[ht] |
|---|
| | 526 | \centering |
|---|
| | 527 | \includegraphics{coniqueBS.png} |
|---|
| | 528 | \caption{Projection conique - Hémisphère sud.} |
|---|
| | 529 | \end{figure} |
|---|
| | 530 | \par\`A noter que cette fois-ci, le maillage n'est pas |
|---|
| | 531 | facultatif\footnote{On pourrait peut-être ajouter des graduations}. |
|---|