Loxodromie

Un article de Wikivisual, l'encyclopédie libre.

Sommaire

1 Présentation générale
2 Navigation loxodromique
3 Démonstration mathématique
4 Voir aussi

[modifier] Présentation générale

Une loxodromie (en grec course oblique) est une courbe qui coupe les méridiens sous un angle constant.

Une route loxodromique est représentée sur une carte marine ou aéronautique en projection de Mercator par une ligne droite mais ne représente pas la distance la plus courte entre deux points. En effet la route la plus courte est appelée route orthodromique ou orthodromie.

La route loxodromique est une route à cap constant.

Image:Loxodromie.svg

Image:Loxo2.svg

[modifier] Navigation loxodromique

Le problème posé est celui de la détermination de la route et de la distance loxodromique entre deux points. Il s'agit donc du problème inverse de la navigation à l'estime.

si les deux points A et B sont peu éloignés, on peut se contenter de formules approchées (latitude moyenne):

<math>M\,</math> étant la distance parcourue à la route <math>R_v\,</math> ; <math>\varphi_A , G_A\,</math> et <math>\varphi_B , G_B\, </math>les coordonnées géographiques (latitude, longitude) des points A et B, et <math>\varphi_m = \frac{\varphi_A + \varphi_B}{2}\,</math> :

<math>\tan R_v = \frac{G_B - G_A}{\varphi_B - \varphi_A} \cos \varphi_m\,</math>

et : <math>M = \frac{\varphi_B - \varphi_A}{\cos R_v}\,</math>

ces formules approchées restent précises à 1 nautique (1' d'arc) près pour <math>M < 375\,</math> nautiques.

formules exactes (latitudes croissantes de la projection de Mercator) :

<math>\tan R_v = - \frac{G_B - G_A}{\lambda_B - \lambda_A}\,</math>

et : <math>M = \frac{\varphi_B - \varphi_A}{\cos R_v}\,</math>

<math>\lambda \,</math>, en ' d'arc, est appelé la latitude croissante (autrefois donné dans les tables de Friocourt) ;

<math>\lambda = 7915,7 . \ln \tan(\frac{\pi}{4} + \frac{\varphi}{2})\, </math>

[modifier] Démonstration mathématique

Sur le globe terrestre, les loxodromies correspondent (lorsqu'elles ne sont pas « dégénérées », c'est-à-dire lorsque l'angle initial donné n'est pas nul) à des spirales s'enroulant autour du pôle (le pôle Nord si l'angle initial et dans <math>]0, \pi[</math> et que le déplacement se fait dans le sens des latitudes croissantes).

Soit à déterminer une équation de la loxodromie et à calculer la longueur parcourue à cap constant <math>\alpha\in\,]0, \pi[</math>.

Considérons les coordonnées sphériques habituelles sur la sphère unité : la longitude <math>\varphi</math> et la colatitude <math>\theta</math>. La loxodromie constitue un arc sur la sphère que l'on suppose de classe <math>C^1</math> : <math>\varphi\mapsto\theta(\varphi)</math> ; soit la fonction <math>f : \varphi\mapsto M(\varphi,\theta(\varphi))</math> qui à la longitude <math>\varphi</math> associe le point courant de la loxodromie de longitude <math>\varphi</math> et de colatitude <math>\theta(\varphi)</math>. Il faut donc bien-sûr se donner au départ une origine des longitudes, puisqu'à un <math>\varphi</math> donné à <math>2\pi</math> près, correspond une infinité de points distincts sur l'arc, de colatitudes différentes. Partons de l'équateur et suivons la loxodromie vers le pôle Nord en nous refusant les classes modulo <math>2\pi</math> pour <math>\varphi</math> : par exemple <math>\theta(\varphi=0)={\pi \over 2}, 0<\theta(\varphi=2\pi)<{\pi \over 2}</math>.

Un vecteur tangent à la loxodromie est ainsi <math>f'(\varphi) = {\partial \vec M \over \partial \varphi}(\varphi,\theta(\varphi)) + \theta'(\varphi)\cdot{\partial \vec M \over \partial \theta}(\varphi,\theta(\varphi))</math>. Ce vecteur, qui dirige la tangente à l'arc, forme donc, par hypothèse, un angle <math>\alpha</math> avec tout vecteur (non nul) dirigeant le parallèle au point considéré. Un vecteur dirigeant le parallèle en <math>M(\varphi,\theta(\varphi))</math> est <math>{\partial \vec M \over \partial \varphi}(\varphi,\theta(\varphi))</math> (tandis qu'un vecteur dirigeant le méridien est bien-sûr <math>{\partial \vec M \over \partial \theta}(\varphi,\theta(\varphi))</math>).

Dans la suite, pour alléger l'écriture, on ne précisera plus le point <math>(\varphi,\theta(\varphi))</math> auquel sont prises les fonctions et leurs dérivées partielles.

En effectuant le produit scalaire d'un vecteur directeur de la tangente à la loxodromie et d'un vecteur directeur du parallèle, on obtient le produit des normes des ces vecteurs par le cosinus de l'angle qu'ils forment :

<math>\left({\partial \vec M \over \partial \varphi}\;|\;{\partial \vec M \over \partial \varphi} + \theta'(\varphi)\cdot{\partial \vec M \over \partial \theta}\right)=\|{\partial \vec M \over \partial \varphi}\|\,\|{\partial \vec M \over \partial \varphi} + \theta'(\varphi)\cdot{\partial \vec M \over \partial \theta}\|\cdot \cos \alpha</math>, en notant <math>(\vec u\;|\;\vec v)</math> le produit scalaire <math>\vec u</math> par <math>\vec v</math>.

En élevant au carré :

<math>\left({\partial \vec M \over \partial \varphi}\;|\;{\partial \vec M \over \partial \varphi} + \theta'(\varphi)\cdot{\partial \vec M \over \partial \theta}\right)^2=\|{\partial \vec M \over \partial \varphi}\|^2\,\|{\partial \vec M \over \partial \varphi} + \theta'(\varphi)\cdot{\partial \vec M \over \partial \theta}\|^2\cdot \cos^2 \alpha</math>.

On a d'autre part clairement : <math>{\partial \vec M \over \partial \varphi}\bot{\partial \vec M \over \partial \theta}</math> (les parallèles et les méridiens sont orthogonaux). Donc, par application du théorème de Pythagore, l'expression se réduit à :

<math>\left({\partial \vec M \over \partial \varphi}\;|\;{\partial \vec M \over \partial \varphi}\right)^2=\|{\partial \vec M \over \partial \varphi}\|^2\,\left(\|{\partial \vec M \over \partial \varphi}\|^2 + \theta'^2(\varphi)\|\cdot{\partial \vec M \over \partial \theta}\|^2\right)\cdot \cos^2 \alpha</math>.

Et en simplifiant :

<math>\|{\partial \vec M \over \partial \varphi}\|^2=\left(\|{\partial \vec M \over \partial \varphi}\|^2 + \theta'(\varphi)^2\cdot\|{\partial \vec M \over \partial \theta}\|^2\right)\cdot \cos^2 \alpha</math>.

D'où, avec « <math>1-\sin^2 =\cos^2</math> »

<math>\sin^2 \alpha \cdot\|{\partial \vec M \over \partial \varphi}\|^2= \theta'(\varphi)^2\cdot\|{\partial \vec M \over \partial \theta}\|^2\cdot \cos^2 \alpha\qquad \mathbf{(1)}</math>.

Calculons les deux normes intervenant dans cette équation :

on sait, d'après le paramétrage sphérique rapporté aux coordonnées cartésiennes dans la base <math>(\vec i, \vec j,\vec k)</math>, que <math>\overrightarrow{OM}(\varphi,\theta)=\cos \theta\;\vec k + \sin \theta\;\vec u_\varphi </math>, où <math>\vec u_\varphi</math> est le vecteur unitaire radial du plan équatorial défini par : <math>\vec u_\varphi = \cos \varphi\; \vec i + \sin\varphi\; \vec j</math>. On définit <math>\vec v_\varphi</math> comme le vecteur dérivé par rapport à <math>\varphi</math> de <math>\vec u_\varphi</math> : <math>\vec v_\varphi = {d\vec u_\varphi\over d\varphi}=-\sin \varphi\; \vec i + \cos\varphi\; \vec j</math>. Alors <math>{\partial \vec M \over \varphi} = \sin \theta\;\vec v_\varphi</math> et <math>{\partial \vec M \over \theta} = -\sin \theta\;\vec k + \cos \theta \vec u_\varphi</math>. Ainsi, <math>\|{\partial \vec M \over \partial \varphi}\|=\sin \theta</math> et <math>\|{\partial \vec M \over \partial \theta}\|=1</math>.

L'équation <math>\mathbf{(1)}</math> se réduit à :

<math>\sin^2 \alpha \sin^2\theta(\varphi)=\theta'^2(\varphi)\cos^2\alpha</math>

et puisque l'on a supposé un trajet vers le pôle Nord, <math>\theta</math> est une fonction décroissante de <math>\varphi</math> et <math>\theta'<0</math>, on suppose là de plus <math>\alpha\in\,]0, \pi/2[</math>(dans les autres cas, on déduit l'arc par une symétrie centrale et/ou une rotation convenable(s), donc on ne perd pas de généralité), par suite :

<math>\sin \alpha \sin\theta(\varphi)=-\theta'(\varphi)\cos\alpha</math>

et <math>\tan \alpha \sin\theta(\varphi)=-{d\theta\over d\varphi}</math>, équation différentielle non linéaire à variables séparables en <math>\theta(\varphi)</math>

En séparant les variables et en intégrant entre 0 et <math>\varphi</math> :

<math>\int_{\pi/2}^{\theta(\varphi)}{d\theta\over \sin\theta}=-\tan\alpha\int_0^\varphi d\varphi</math>

<math>\ln\left(\tan{\theta(\varphi)\over 2}\right)=-\tan \alpha \varphi</math> (cf. Table de primitives)

<math>\theta(\varphi)=2\,\operatorname{Arctan}\left(\mathrm{e}^{-\varphi\,tan\alpha}\right)</math>

La longueur L parcourue vaut alors, par définition :

<math>L=\int_0^{+\infty}\|f'(\varphi)\|d\varphi</math>

où <math>f'(\varphi) = {\partial \vec M \over \partial \varphi}(\varphi,\theta(\varphi)) + \theta'(\varphi)\cdot{\partial \vec M \over \partial \theta}(\varphi,\theta(\varphi))</math> et <math>\|f'(\varphi)\|^2 = \sin^2 \theta + \theta'^2 = \sin^2\theta + \tan^2\alpha\sin^2\theta={\sin^2\theta\over\cos^2\alpha}</math> et pour les mêmes raisons de signe, <math>\|f'(\varphi)\|={\sin\theta\over\cos\alpha}</math>.

<math>L={1\over \cos\alpha}\int_0^{+\infty}\sin\theta(\varphi)d\varphi</math>

En changeant de variable, avec <math>{d\varphi\over d\theta}=-{1\over \tan\alpha\,\sin\theta}</math>, on a <math>L={1\over\cos\alpha\tan\alpha}\int_0^{\pi\over 2}\sin^2\theta\,d\theta</math>