Loxodromic curve

A loxodromic curve (of the Greek grc lox (O) - and - grc dromy race obliques ) is a curve which cut the Méridien S under a constant angle .

A loxodromic road is represented on a Sea chart or aeronautics in Projection of Mercator by a straight line but does not represent the the shortest distance between two points. Indeed the shortest road is called orthodromic road or Orthodromie.

The loxodromic road is a road with constant course.

Loxodromic navigation

The problem arising is that of the determination of the road and the loxodromic distance between two points. It is thus about the opposite problem of the Navigation to the regard.

  • if the two points has and B is not very distant, one can be satisfied with approximate formulas (average latitude):

M \, being the distance covered with the road R_v \, ; \ varphi_A, G_A \, and \ varphi_B, G_B \, les coordinated geographical (latitude, longitude) of the points has and B, and \ varphi_m = \ frac {\ varphi_A + \ varphi_B} {2} \, :

\ tan R_v = \ frac {G_B - G_A} {\ varphi_B - \ varphi_A} \ cos \ varphi_m \,

and: M = \ frac {\ varphi_B - \ varphi_A} {\ cos R_v} \,

these approximate formulas remain precise near with 1 nautical (1 ' of arc) for M < 375 \, nautical.

  • exact formulas ( increasing latitudes of the projection of Mercator):

\ tan R_v = - \ frac {G_B - G_A} {\ lambda_B - \ lambda_A} \,

and: M = \ frac {\ varphi_B - \ varphi_A} {\ cos R_v} \,

\ lambda \, , in 'arc, is called the increasing latitude (formerly given in the tables of Friocourt);

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

Mathematical demonstration

On the terrestrial sphere, the loxodromic curves correspond (when they “are not degenerated”, i.e. when the initial angle given is not null) to spirals being rolled up around the pole (the north pole if the initial angle and in ] 0, \ pi and that displacement is done in the direction of [[latitude] S increasing).

Is to be determined an equation of the loxodromic curve and to calculate the length traversed with constant course \ alpha \ in \,] 0, \ pi

Let us consider the usual Coordonnée S spherical on the sphere unit: the Longitude \ varphi and the colatitude \ theta. The loxodromic curve constitutes a arc on the sphere which one supposes of class C^1: \ varphi \ mapsto \ theta (\ varphi) ; that is to say the function f: \ varphi \ mapsto M (\ varphi, \ theta (\ varphi)) which with longitude \ varphi associates the point running of the loxodromic curve of longitude \ varphi and of colatitude \ theta (\ varphi) . It is thus necessary well-sure to give oneself at the beginning an origin longitudes, since to a \ varphi given to 2 \ pi close, corresponds an infinity of distinct points on the arc, of different colatitudes. Let us leave the equator and follow the loxodromic curve towards the north pole by refusing us the classes modulo 2 \ pi for \ varphi: for example \ theta (\ varphi=0) = {\ pi \ over 2}, 0< \ theta (\ varphi=2 \ pi) < {\ pi \ over 2} .

A tangent vector with the loxodromic curve is thus f' (\ varphi) = {\ partial \ vec M \ over \ partial \ varphi} (\ varphi, \ theta (\ varphi)) + \ theta' (\ varphi) \ cdot {\ partial \ vec M \ over \ partial \ theta} (\ varphi, \ theta (\ varphi)). This vector, which directs the tangent to the arc, thus forms, by assumption, an angle \ alpha with any vector (not no one) directing it parallel with the point considered. Vector directing it a parallel in M (\ varphi, \ theta (\ varphi)) is {\ partial \ vec M \ over \ partial \ varphi} (\ varphi, \ theta (\ varphi)) (while a vector directing the meridian line is well-sure {\ partial \ vec M \ over \ partial \ theta} (\ varphi, \ theta (\ varphi))).

In the continuation, to reduce the writing, one will not clarify any more the point (\ varphi, \ theta (\ varphi)) to which are taken the partial functions and their derivative.

By carrying out the scalar product of a directing vector of the tangent to the loxodromic curve and of a directing vector of the parallel, one obtains the product of the standards of these vectors by the cosine of the angle which they form:

\ 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, by noting (\ vec U \; |\; \ vec v) the scalar product \ vec u by \ vec v.

While raising squared:

\ 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.

One has in addition clearly: {\ partial \ vec M \ over \ partial \ varphi} \ club-footed {\ partial \ vec M \ over \ partial \ theta} (the parallels and the meridian lines are orthogonal). Therefore, by application of the theorem of Pythagore, the expression is reduced to:

\ 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.

And while simplifying:

\|{\ 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.

From where, with “1- \ sin^2 = \ cos^2

\ 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)}.

Let us calculate the two standards intervening in this equation:

it is known, according to the spherical parameter setting brought back to the Cartesian coordinates in the base (\ vec I, \ vec J, \ vec K) , that \ overrightarrow {OM} (\ varphi, \ theta) = \ cos \ theta \; \ vec K + \ sin \ theta \; \ vec u_ \ varphi , where \ vec u_ \ varphi is the radial unit vector of the equatorial plan defined by: \ vec u_ \ varphi = \ cos \ varphi \; \ vec I + \ sin \ varphi \; \ vec j. One defines \ vec v_ \ varphi like the vector derived compared to \ varphi from \ vec u_ \ varphi: \ vec v_ \ varphi = {D \ vec u_ \ varphi \ over D \ varphi} = \ sin \ varphi \; \ vec I + \ cos \ varphi \; \ vec j. Then {\ partial \ vec M \ over \ varphi} = \ sin \ theta \; \ vec v_ \ varphi and {\ partial \ vec M \ over \ theta} = - \ sin \ theta \; \ vec K + \ cos \ theta \ vec u_ \ varphi. Thus, \|{\ partial \ vec M \ over \ partial \ varphi} \|= \ sin \ theta and \|{\ partial \ vec M \ over \ partial \ theta} \|=1.

The equation \ mathbf {(1)} is reduced to:

\ sin^2 \ alpha \ sin^2 \ theta (\ varphi) = \ theta'^2 (\ varphi) \ cos^2 \ alpha
and since one supposed a way towards the north pole, \ theta is a decreasing function of \ varphi and there \ theta'<0, one supposes moreover \ alpha \ in \,] 0, \ pi/2 the other cases, one deduces the arc by a central symmetry and/or a suitable rotation (S), therefore one does not lose a general information), consequently: : \ sin \ alpha \ sin \ theta (\ varphi) = \ theta' (\ varphi) \ cos \ alpha
and \ tan \ alpha \ sin \ theta (\ varphi) = {D \ theta \ over D \ varphi} , nonlinear differential equation with separable variables in \ theta (\ varphi)

By separating the variables and while integrating between 0 and \ varphi:

\ int_ {\ pi/2} ^ {\ theta (\ varphi)}{D \ theta \ over \ sin \ theta} = \ tan \ alpha \ int_0^ \ varphi D \ varphi
\ ln \ left (\ tan {\ theta (\ varphi) \ over 2} \ right) = \ tan \ alpha \ varphi (cf Table of primitives)
\ theta (\ varphi) =2 \, \ operatorname {Arctan} \ left (\ mathrm {E} ^ {- \ varphi \, tan \ alpha} \ right)

The length L traversed is worth then, by definition:

L= \ int_0^ {+ \ infty} \|f' (\ varphi) \|D \ varphi
where f' (\ varphi) = {\ partial \ vec M \ over \ partial \ varphi} (\ varphi, \ theta (\ varphi)) + \ theta' (\ varphi) \ cdot {\ partial \ vec M \ over \ partial \ theta} (\ varphi, \ theta (\ varphi)) and \|f' (\ varphi) \|^2 = \ sin^2 \ theta + \ theta'^2 = \ sin^2 \ theta + \ tan^2 \ alpha \ sin^2 \ theta= {\ sin^2 \ theta \ over \ cos^2 \ alpha} and for the same reasons of sign, \|f' (\ varphi) \|= {\ sin \ theta \ over \ cos \ alpha} .
L= {1 \ over \ cos \ alpha} \ int_0^ {+ \ infty} \ sin \ theta (\ varphi) D \ varphi
While changing variable, with {D \ varphi \ over D \ theta} = {1 \ over \ tan \ alpha \, \ sin \ theta} , one has L= {1 \ over \ cos \ alpha \ tan \ alpha} \ int_0^ {\ pi \ over 2} \ sin^2 \ theta \, D \ theta
L= {\ pi \ over 4 \, \ sin \ alpha}

See too

Random links:List catholic ecclesiastical provinces | Frank Furness | Playback | November 11th in sport | Coslett Herbert Waddell | Sarcome