Loxodromie
Loxodromie
Loxodromie
1 Position du problème
Un véhicule se déplace à la surface de la Terre (supposée sphérique de rayon RT = 6370 km) en gardant un
angle α0 constant avec les parallèles (à l’équateur) qu’il rencontre. Cette trajectoire s’appelle une loxodromie.
On repère la position du mobile en coordonnées sphériques (r, θ, ϕ) dans la base (− →r , −
u →θ , −
u u→
ϕ ) où θ est la
colatitude. On définit également la latitude λ telle que θ − λ = π/2.
1. Montrer que l’équation d’une loxodromie partant d’un point de coordonnées (RT , θ0 , ϕ0 ) s’écrit:
tan θ2 = tan θ20 .etan α0 .(ϕ−ϕ0 ) .
2. Tracer la loxodromie entre deux villes: Brest (A: λA = 48, 3◦ N, ϕA = 4, 5◦ W) et Fort de France (B:
λB = 14, 6◦ N, ϕB = 61◦ W).
On peut également repérer la position du mobile en coordonnées cartésiennes dans la base (~ı,~, ~k) où ~k
pointe vers le pôle Nord et ~ı vers l’intersection du méridien de Greenwich avec le plan équatorial.
On appelle orthodromie la trajectoire la plus courte allant du point A au point B. Il s’agit d’un arc de cercle
dans un plan Π où Π est le plan, passant par le centre de la Terre, perpendiculaire à ~uΠ dont les coordonnées
dans la base (~ı,~, ~k) sont données par:
cos λA . cos ϕA cos λB . cos ϕB
~urA cos λA . sin ϕA ∧ ~urB cos λB . sin ϕB
sin λA sin λB
4. Calculer la distance d entre A et B. On montrera que, dans le plan Π, l’orthodromie correspond à un angle
ψ tel que cos ψ = ~urA .~urB . On en déduira que d = RT .ψ.
Solution
1. En coordonnées sphériques:
lorsque θ varie de dθ, on se déplace de RT .dθ vers le Sud;
lorsque ϕ varie de dϕ, on se déplace de RT . sin θ.dϕ vers l’Est.
Si on coupe les parallèles sous un angle α0 , alors: tan α0 = RTR. sin
T .dθ
θ.dϕ .
dθ
En séparant les variables, on a: sin θ = tan α0 .dϕ. On obtient une intégrale classique que l’on résout en
2. cos2 θ2 .du
posant u = tan θ2 ⇒ du = 12 . cos12 θ .dθ. Ainsi, dθ
sin θ = 2. sin
dθ
θ θ = 2. sin θ2 . cos θ2
= du
tan θ2
= du
u s’intègre en
2 . cos 2
ln |u| + cte.
En intégrant entre le point de départ A et un point quelconque de la trajectoire, on a:
θ
ln tan θ2 θ = tan α0 .(ϕ − ϕA ) ⇒ tan θ2 = tan θ20 . exp(tan α0 .(ϕ − ϕ0 ): c’est l’équation de la loxodromie.
A
4. cos ψ = ~urA .~urB = cos λA . cos ϕA . cos λB . cos ϕB + cos λA . sin ϕA . cos λB . sin ϕB + sin λA . sin λB = 0, 54.
On prend la fonction réciproque (arccos) et on obtient: ψ = 57, 07◦ = 0, 996 rad ⇒ d = RT .ψ = 6.345 km.
Equation de la loxodromie
In[1]:= Needs["LinearAlgebra‘CrossProduct‘"]
Needs["Calculus‘VectorAnalysis‘"]
Needs["Graphics‘ParametricPlot3D‘"]
Needs["ProgrammingExamples‘SphericalCurve‘"]
1
ISEN-Brest. Kany. TD: Loxodromie et Orthodromie
In[5]:= alpha0=.;
Simplify[Integrate[1/Sin[theta],{theta,theta0,theta}]
==Integrate[Tan[alpha0],{phi,phi0,phi}]]
Tracé de l’hémisphère
Tracé de la loxodromie
thetaA=Brest[[1]];phiA=Brest[[2]]; thetaB=FortDeFrance[[1]];phiB=FortDeFrance[[2]];
alpha0=N[ArcTan[Log[Abs[Tan[thetaA/2]/Tan[thetaB/2]]]/(phiA-phiB)]];
In[14]:= Trajet={SphericalPoint[1,thetaB,phiB],SphericalPoint[1,thetaA,phiA]}
Orthodromie
In[21]:= thetaC=ArcCos[N[uPi[[3]]]];thetaC*180/N[Pi]
ArcCos[N[uPi[[1]]]/Sin[thetaC]]*180/N[Pi]
ArcSin[N[uPi[[2]]]/Sin[thetaC]]*180/N[Pi]
phiC=-ArcCos[N[uPi[[1]]]/Sin[thetaC]];phiC*180/N[Pi]
Out[21]= 50.2439
Out[22]= 163.514
Out[23]= -16.4855
Out[24]= -163.514
Inverse[M[thetaC,phiC]].{1,0,0}
2
ISEN-Brest. Kany. TD: Loxodromie et Orthodromie
M[thetaC,phiC].uPi
In[31]:= Graph3=SphericalCircle3[thetaC,phiC]
Out[32]= -Graphics3D-
Distance
Out[33]= 6345.66
3
ISEN-Brest. Kany. TD: Loxodromie et Orthodromie
Projection de Mercator
orthodromie
loxodromie
0.6
0.4
0.2
0.0
0.2
0.4
0.6 0.6
0.4
0.2
0.0 0.6
0.2
0.4 0.0 0.2 0.4
0.2
0.6 0.6 0.4